[ cavity | charge | clipping plane | coil | column | Collection | current map | current object | Current | command | comp_matrix | conf | cpk ]
icmCavityFinder macro which identifies and characterizes cavities.
To generate a cavity an
skin is generated and converted into grob .
We then split this grob into individual connected entities.
set charge formal
commands.
Formal charges in functional groups can be changed in chemical arrays with the modify command.
The following variables and commands are related to charges in 3D objects or chemical arrays:
set charge formalset chargeset chargemmffCharge( .. [formal] )GRAPHICS.formalChargeDisplay- modify : find and replace charge and uncharged chemical functional groups
"FOG" mode the image gradually disappears between the front and the back planes.
ribbon
type of representation is used. To assign a residue fragment to coil,
do something like this:
assign sstructure
a_/13:40 "_".
#> a b c
1. 15 "alpha"
2. 18 "beta"
5. 10 "gamma"
It is also useful to show arrays of a
table
(see also
database).
Collection is composed of a array of unique keys and a array of values, where each key is associated with one value. Each value can either be of basic type or other collection or array of collections.
A collection can be created in different ways:
- using
Collectionfunction - from JSON format. (`read-json
Collection) - -- read
xmlfrom XML source (read xmlfile.xml orCollection)
Examples:
c = Collection( "firstName", "John", "lastName", "Snith", "age" 25 ) # from key/value pairs
c = Collection( "{ 'firstName': 'John', 'lastName': 'Smith', 'age': 25 }" ) # from JSON string
c = Collection( "<firstName>John</firstName><lastName>Smith</lastName><age>25</age>" ) # from XML string
c["age"]
c["firstName"]
write xml c "c.xml"
Elements can be added, modified and removed dynamically
c = Collection() # empty collection
c["count"] += 1 # accumulates count (when used the first time, creates a new key
c["name"] = "John" # creates a new element
c["name"] = "Jack" # overwrites an existing element
c["addressBook"] = Array()
c["addressBook",1] = Collection( "name" "Jack" "tel" "858-3452238" "cell" "858-8872676" )
c["addressBook",2] = Collection( "name" "John" "tel" "858-7878795" "cell" "858-6373913" )
# index access
c["addressBook",1] # first entry of the addressbook
c["addressBook",1:2] # first two entries of the addressbook
c["addressBook",1:2,"tel"] # first two 'tel' fields
# deletion
delete c["addressBook",2] # deletes second element of the adderssbook
delete c["addressBook"]
See:
an instruction one can execute in the ICM-shell interactively or from
an ICM script file. Typically a command consists of a verb (like
The word order in the argument list is not important,
if arguments have different types. For example:
If two or more argument of the same type are present the order becomes important, e.g.
Redefining global preferences and parameters temporarily for one command.
Sometimes you want to run a command with different values of global
See also:
Individual properties, such as energies, number of visits and the total
numbers of visits can be re-set with the set stack command.
The stack can also be stored in the object it is associated with. It enables multiple object sessions, each object with its own stack, and ability to save the object with multiple conformations to a file. To store the stack in the object use the store conf a_ command.
If cpk and
The spheres sliced by a clipping plane can be displayed in several ways, see
Collection read xml , read json , write-xml, String( xml drugbank example
make map family of commands.
The current map may be changed by the set map command.
To find out which one is the current, use list maps.
build command.
The current object may be changed by the set object [ os_object ] command.
To find out which is the current one, use list object or show object [ os_object].
You may refer to the current object as a_
list tables.
read or delete ) and a bunch of arguments.
display a_//ca,c,n yellow # is as good as ...
display yellow a_//ca,c,n # ... inverse
write "This line" "filename"
or
superimpose as_select1 as_select2 ,
grid["x"] in the plot command.
In the latter case if the "x" string is not preceded by grid word, it will not be considered.
preference and global shell variables of the following type: integer, real, string and logical .
Their values will be automatically restored after the completion
of the command.
E.g.
show surface area surfaceAccuracy=10 vwExpand=1.4
show energy electroMethod = "Coulomb"
file and can be reloaded later
(if you need another matrix) using the read comp_matrix command.
The current matrix can be shown.
All elements are properly reordered and multiplied by a factor
such that the sum of occurrence-weighted diagonal elements (identities)
of the matrix is 1.0 and stored in this form. You can increase all the
values in the matrix (a useful operation for low sequence identities)
with the set comp_matrix command which has similar effect
to reduction of the gapOpen parameter.
The set comp_matrix command also allows one to modify weight
for specific residue pairs (e.g. set comp_matrix 4. "CC").
Two matrices are recommended: the blosum50 with penalties 2.0/0.15 and the
Gonnet matrix (the default) with penalties 2.4/0.15.
A number of matrices are provided. You can use your own data file too.
Normalized matrices with diagonal elements equal to 1. Most of the matrices efficient for sequence alignment have different Cii
weights. For example, two aligned tryptophanes will get better score than
two aligned alanines. The Cij values can also be normalized as follows:
This normalization is used in the Cij_norm = Cij/Sqrt(Cii*Cjj)Score () and Rarray () functions
calculating conservation values for each position in an alignment.
stack .
You can add a conformation to the conformational stack or
load a conformation from the stack.
Example:
% icm
build string "ala his trp"
montecarlo
show stack
iconf> 1 2 3 4 5 6 7
ener> -15.1 -14.6 -14.6 -14.2 -13.9 -11.4 -1.7
rmsd> 0.3 39.2 48.0 44.1 27.4 56.6 39.3
naft> 1 0 0 1 1 1 0
nvis> 4 1 1 4 4 4 1
print Nof(conf)
7
print Iarray(stack) # returns {4 1 1 4 4 4 1}
t = Table(stack) # more conf parameters
for i=1,Nof(conf)
load conf i
center
pause # hit return
endfor
atoms.
In ICM graphics it means solid representation of van der Waals spheres.
Spheres are drawn at van der Waals radii. See also:
display
xstick are displayed simultaneously the cpk sphere becomes transparent. Undisplay xstick to
make the spk solid again.
GRAPHICS.cpkClipCaps .