Equilibration protocol for globular and membrane proteins It includes a flatbottom potential box to retrain a ligand for example within this box.
Number of steps to run the simulations in units of 4fs
Temperature of the thermostat in Kelvin
Force constant of the flatbottom potential in kcal/mol/A^2. E.g. 5
Reference selection to use as dynamic center of the flatbottom box.
Selection of atoms to apply the flatbottom potential
Position of the flatbottom box in term of the reference center given as [xmin, xmax, ymin, ymax, zmin, zmax]
For membrane protein simulations set it to true so that the barostat does not modify the xy aspect ratio.
A dictionary containing as keys the atomselections of the constraints and as values the constraint scaling factor. 0 factor means no constraint, 1 full constraints and in between values are used for scaling. The order with which the constraints are applied is random, so make atomselects mutually exclusive to be sure you get the correct constraints.
>>> from htmd.protocols.equilibration_v1 import Equilibration >>> md = Equilibration() >>> md.numsteps = 10000000 >>> md.temperature = 300 >>> md.useconstantratio = True # only for membrane sims >>> # this is only needed for setting the flatbottom potential, otherwise remove it >>> md.reference = 'protein and resid 293' >>> md.selection = 'segname L and noh' >>> md.box = [-25, 25, -25, 25, 43, 45] >>> md.k = 5 >>> md.write('./build','./equil')
Convenience function for adding a new constraint to existing constraints.
>>> eq.addConstraint('chain X', 0.3)
Write the equilibration protocol
Writes the equilibration protocol and files into a folder for execution using files inside the inputdir directory
>>> md = Equilibration() >>> md.write('./build','./equil')