Sample MDP file for C-Alpha Gromacs simulations
;Run control: Use Langevin dynamics.
integrator = sd
;time step in reduced units. Since we use mass of a carbon=1, the time step is smaller that many other MD simulations.
dt = 0.0005
;number of integration steps
nsteps = 100000
;frequency to write coordinates to output trajectory .trr file
nstxout = 100000
;frequency to write velocities to output trajectory .trr file
nstvout = 100000
;frequency to write energies to log file
nstlog = 1000
;frequency to write energies to energy file
nstenergy = 1000
;frequency to write coordinates to .xtc trajectory
nstxtcout = 1000
;group(s) to write to .xtc trajectory.
;NOTE: group name "system" is defined if you are NOT using an index file. If you use an index file, all group references need to be consistent with the names given in the index file. Our webpage, by default, names the chains "1", "2", ... You can change these names to whatever you like, just make sure the mdp file uses the same naming scheme. For more information on group usage, please consult the Gromacs manual.
xtc_grps = system
;group(s) to write to energy file
energygrps = system
;Frequency to update the neighbor list (and the long-range forces,
;when using twin-range cut-off's).
nstlist = 20
;Make a grid in the box and only check atoms in neighboring grid cells
;when constructing a new neighbor list every nstlist steps.
ns_type = grid
;cut-off distance for the short-range neighbor list (NOTE: cut-offs should be longer if running a CA model than an AA model). 3.0 nm is about 2.5 times the typical longest interaction in a CA model. If you have longer contacts, then increase these numbers.
rlist = 3.0
;treatment of electrostatic interactions
rcoulomb = 3.0
;treatment of van der waals interactions
rvdw = 3.0
;this tells gromacs to use the tables that we give it, instead of the standard interaction types
coulombtype = User
vdwtype = User
; Periodic boudary conditions in all the directions
pbc = xyz
table-extension = 20 ; (nm) !IMPORTANT! Set this so that rvdw + table-extension equals half of the box's longest diagonal.
;Temperature coupling
tc-grps = system
tau_t = 1.0 ; Temperature coupling time constant. Smaller values = stronger coupling.
ref_t = 50.0 ; In reduced units (see Gromacs Manual for explanation. We find many proteins fold around T=100-120)
;Pressure coupling
Pcoupl = no
;Velocity generation
gen_vel = yes
gen_temp = 50.0
gen_seed = 30