MMTSB
Tool Set Documentation

Difference between revisions of "hlamc.pl"

From MMTSB
Jump to: navigation, search
Line 1: Line 1:
__NOTOC__
+
dummy
QqMPan  <a href="http://wciwqozvaykk.com/">wciwqozvaykk</a>, [url=http://zorhhsphbbmh.com/]zorhhsphbbmh[/url], [link=http://iqlxelzaugym.com/]iqlxelzaugym[/link], http://vpyqztywhdob.com/
 
 
 
== Description ==
 
 
 
This utility is used for running single hybrid lattice/all-atom simulations. A replica exchange
 
version <docmark>hlamcrex.pl</docmark> is also available.<BR>
 
Two Monte Carlo simulation schemes are available for hybrid simulations that combine lattice and all-atom models.
 
In both modes lattice models are effectively simulated on an all-atom surface combining the advantages of
 
reduced lattice models with a more accurate all-atom potential energy function.<BR>
 
In the first mode, called the "free" mode, a Monte Carlo simulation cycle consists of a lattice simulation and all-atom
 
reconstruction of the final lattice chain followed by all-atom minimization and molecular dynamics.
 
The all-atom energy at the end of each cycle is used for accepting or rejecting conformations according to Metropolis
 
criteria.<BR>
 
In the second mode, called "coupled", lattice and all-atom simulations are more tightly coupled.
 
Here, final conformations from very short lattice simulations are used for generating side chain
 
center of mass restraints in continuous all-atom minimization or molecular dynamics runs.
 
As in the first mode, all-atom energies are used for accepting or rejecting conformations in a Monte
 
Carlo scheme. The advantage of the second method is that the all-atom structures remain mostly
 
relaxed throughout the simulation while following the guiding restraints from the lattice simulations.
 
However, one faces the same barriers due to steric hindrance as in conventional all-atom simulations
 
that are easily overcome with the reduced lattice model in the first mode.<BR>
 
The default mode is "free". It can be selected with the <font color=#508060>mode</font> keyword of the <B>-par</B> option.
 
Other parameters are the Metropolis 'temperature' in units of kcal/mol that is set with <font color=#508060>mctemp</font>,
 
the force constant that is used for C-alpha restraints in the all-atom part during the "free" mode,
 
the force constant used in the "coupled" mode for side chain center restraints,
 
and two additional options <font color=#508060>[no]keepchain</font> and <font color=#508060>[no]eval</font>. The first one sets whether lattice simulations
 
are continued from the final chain of the previous lattice run (<font color=#508060>keepchain</font>, default)
 
or from a chain generated at the end of the all-atom simulation section (<font color=#508060>nokeepchain</font>).
 
The second option controls whether a separate instance of CHARMM is used for evaluating all-atom energies at
 
the end of the simulation cycle(<font color=#508060>eval</font> vs. <font color=#508060>noeval</font>). This is useful when a different,
 
more costly function, e.g. using an implicit solvation term, should be used for evaluating energies while
 
a simple vacuum energy function may be sufficient for restrained all-atom minimizations and short molecular dynamics
 
runs. The number of total Monte Carlo cycles is set with <B>-n</B>.<BR>
 
Lattice simulations require a MONSSTER sequence file as input with <B>-seq</B>. Other lattice simulation
 
parameters such as the number of lattice simulation Monte Carlo cycles, the internal temperature or
 
potential component scaling factors as listed in <docmark>geninput.pl</docmark> are given with <B>-latpar</B>.<BR>
 
All-atom simulation parameters are set with <B>-aapar</B>. The complete list of possible parameters is
 
given <A HREF="http://mmtsb.scripps.edu/cgi-bin/charmmpar" TARGET=cpar>here</A>. Minimization runs are
 
performed if <font color=#508060>sdsteps</font> and/or <font color=#508060>minsteps</font> are set to non-zero values. Molecular dynamics
 
is turned off by default but will be performed if <font color=#508060>dynsteps</font> is set. Parameters for a separate
 
energy evaluation are set with <B>-evalpar</B>.<BR>
 
All output is written to the current directory unless a different directory is specified with <B>-dir</B>.
 
Structures that are accepted at the end of a Monte Carlo cycle are saved into an ensemble under the same
 
directory. The tag that is used can be set with <B>-enstag</B>.<BR>
 
Log files for the all-atom simulation part, all-atom evaluation, and overall hybrid sampling are generated
 
if the options <B>-aalog</B>, <B>-evallog</B>, or <B>-log</B> are used.<BR>
 
 
 
== Options ==
 
 
 
; -help : usage information
 
 
 
 
 
== Examples ==
 
 
 
<mmtsbToolExample cmd="hlamc.pl" set="/apps/mmtsb/bench/hlamc.pl-test/test1"></mmtsbToolExample>
 
 
 
<mmtsbToolExample cmd="hlamc.pl" set="/apps/mmtsb/bench/hlamc.pl-test/test2"></mmtsbToolExample>
 

Revision as of 16:59, 6 December 2008

dummy