aarex.pl

Description

This script is used to run replica exchange simulations. In most parallel environments it will start the replica exchange server automatically. The options -n, -par, -temp, -condfile, -ens, -ensdir are available as in <docmark>rexserver.pl</docmark>. Alternatively, it can connect to an external replica exchange server if its address, port, and ID are given with -rserv. In that case -jobs may be used to run only some of the temperature windows instead of all windows if multiple clients are launched on different machines.
For shared memory environments the option -cpus can be used to specify a smaller number of CPUs than temperature windows if necessary due to computational restraints. A host file can be given with -hosts for automatic remote submission in a distributed environment. (see <docmark>calcprop.pl</docmark> for a more detailed explanation).
In order to contact the server from an external program (e.g. for monitoring purposes) the server ID is required. The option -saveid is available to write this server information to a file.
The intial input PDB files for each MD simulation are expected either as the last command line arguments or from an external file that is given with -f.
The remaining parameters are used to control the MD simulations. Parameters that can be given with -mdpar are shown <A HREF="http://mmtsb.scripps.edu/cgi-bin/charmmpar" TARGET=cpar>here</A>.
The option -l is available to provide a list of residues and a template PDB structure for loop modeling. Please note that with this option RMSD values that are automatically calculated if a reference PDB structure is given are also limited to only these residues.
Further options given with -mdopt control whether a trajectory or restart file is written out ([no]trajout and [no]restout), how many MD steps are used (default: 1) for averaging the energy score used in the replica exchange Metropolis criteria (avgener), whether the rest of protein is restrained outside the loop ([no]conslim), and the maximum restraint force (limforce) and type of restrained atoms (limsel) in this case. Other restraints may be specified with -cons.
The simulation protocol can be further customized by providing CHARMM commands through external files. This can be done with -custom which expects a keyword and a file name as arguments. The keyword is used to specify when the custom command sequence should be inserted in the standard protocol. If setup is used the commands will be sent to CHARMM only once during the initial setup phase. If the keywords pre or post the commands will be executed before or after the dynamics command, respectively. These keywords can be further qualified with :init, :equi, and :prod corresponding to the replica exchange cycle modes to allow custom equilibration protocols.
The option -log is available to request a server log file. A CHARMM log file is generated for each client if -charmmlog is given, an energy log file is generated with -elog.

Options

-help 

usage information

Examples

aarex.pl -n 8 -mdpar dynsteps=100,param=22,gb,nocut -par initruns=2,equilruns=2,natpdb=1vii.exp.pdb -temp 4:298:400 1vii.sample.{1,2,3,4}.pdb
runs 8 replica exchange MD simulation cycles with four exponentially spaced temperature windows from 298 to 400K. The first 2 runs are considered initialization runs, the next 2 runs are equilibration runs. The native PDB structure is given as reference for calculation RMSD values. MD parameters are set to run 100 steps for each cycle, use CHARMM22 parameters with GB implicit solvent and no electrostatic cutoffs. Initial conformations are taken from the files 1vii.sample.?.pdb.

aarex.pl -n 4
runs 4 additional cycles continuing a previous replica exchange simulation run

aarex.pl -n 5 -par initruns=2,equilruns=0,nosave -temp 4:298:400 -ensdir data -ens rex -mdpar dynsteps=200,gb,nocut 1vii.sample.{1,2,3,4}.pdb
runs 5 replica exchange simulation runs with 2 initialization and no equilibration runs. The conformation from the lowest temperature at each run is saved under the rex tag in an ensemble in the directory data. No other conformations during the simulation are saved.

aarex.pl -n 6 -par initruns=2,equilruns=2,natpdb=1vii.exp.pdb -temp 4:298:400 -cpus 2 -mdpar dynsteps=200,gb -l 1vii.exp.pdb 10:21 -mdopt conslim,limforce=5.0,limsel=ca -log server.log -charmmlog charmm.log -f init.files
runs a replica exchange simulation runs for loop modeling. The loop residues are located at 10 through 21, the rest of the protein is restrained at C-alpha atoms with a force constant of up to 5 kcal/mol. This run has 4 temperature windows but uses only 2 CPUs at a time. A server log file and for each client a CHARMM log file are written out. The intial input files are taken from init.files.

aarex.pl -n 8 -mdpar dynsteps=100,gb,nocut -par initruns=2,equilruns=2 -temp 4:298:400 -hosts sgi.workstations 1vii.sample.{1,2,3,4}.pdb
runs a replica exchange simulation across distributed workstations

aarex.pl -n 8 -mdpar dynsteps=100,gb,nocut -par initruns=2,equilruns=2 -temp 4:298:400 -mp -hosts sgi.local 1vii.sample.1.pdb
runs a replica exchange simulation across distributed workstations with local directories. All temperature windows are started from the same initial file.

aarex.pl -n 10 -mdpar dynsteps=100,gb,nocut -par initruns=2,equilruns=0,natpdb=1vii.exp.pdb -condfile conditions 1vii.sample.{1,2,3,4}.pdb
replica exchange simulation with all replicas at the same temperature but with different radius of gyration umbrella potentials

The condition file looks the following:

bias rg
298 force=1.25,target=18.0
298 force=0.60,target=9.0
298 force=0.60,target=13.0
298 force=0.75,target=14.0

aarex.pl -n 10 -mdpar dynsteps=100,gb -par initruns=2,equilruns=0,natpdb=1vii.exp.pdb -condfile conditions -log server.log -charmmlog charmm.log 1vii.sample.{1,2,3,4}.pdb
2D replica exchange with two temperatures and two different radius of gyration umbrellas

with the following condition file:

bias rg
298 force=1.00,target=15.0
298 force=0.60,target=12.0
320 force=0.60,target=12.0
320 force=1.00,target=15.0

aarex.pl -charmmlog clog -n 8 -mdpar dynsteps=100,param=19,nogb -custom setup acesetup.inp -custom pre ace.inp -par initruns=2,equilruns=0,natpdb=1vii.exp.pdb -temp 4:298:350 1vii.sample.2.pdb
runs a replica exchange simulation with a customized potential function. Through the files acesetup.inp and ace.inp the ACE solvation model is used instead of the default Generalized Born model.