|
This script is used to generate structure ensemble structures
through minimization of PDB files. Many options are similar
to minCHARMM.pl, the script for single
minimization runs, and enslatsim.pl, the
script for generating ensembles from lattice simulations.
Two tags are required as argument to ensmin.pl:
The first one identifies the ensemble structures that are used for
input at each run and have to be present already, the second
one is used for the output after the minimization run is complete.
The ensemble directory can be given with -dir if it is
not the current directory. The number of desired runs is specified
using -run either followed by a single number
for the total number of runs or by an explicit interval.
A reference PDB file for comparison with the native conformation
may be given through -natpdb but is not required. For loop or fragment
modeling the range of residues that are being modeled
and a template PDB file with the structure outside the
modeled residues are specified using -l. Like the
native PDB structure this data may already be present in the
ensemble configuration file and does not need to be given in
that case.
Although this script will work in a serial fashion, one would
probably want to take advantage of parallel architectures for
a larger number of simulations. The available parallel options are
listed here.
The remaining options are used to control minimization runs
and are kept in a file tag.options for later reference
and subsequent executions of this script on the same ensemble set.
They are very similar to the options in
minCHARMM.pl with the exception of restraints
for loop modeling:
In minCHARMM.pl the option -l switches
on gradual restraints outside a given fragment list based on distance
in sequence space. An alternative is to use -cons with a
restraint list generated by vicinity.pl. Such a
list will restrain residues surrounding the loop region in two shells
with different force constants to allow some readjustment of residues in
direct contact with the loop. In ensmin.pl a fragment list,
given with -l, is not just used for restraints but also to limit
RMSD and other calculations to the loop region. In typical loop modeling
applications this option would always be specified
even if the restraints are given through -cons. A second
option -conslim is therefore necessary to request explicitly
that the fragment list is also used to generate restraints. However,
if a fragment list is given and no other
restraint list is given with -cons this flag is turned on
automatically, unless -noconslim is given.
The options -limsel and -limforce to specify the
atom selection and maximum force constant used for automatically
restraints from the fragment list. If any of these options is given
the conslim is also turned on automatically.
|