DMD Options in dmdinput.json

dmdinput.json Options

Standard DMD options

The following are dmd options at are similar to the configuration file when running DMD by hand. For a more detailed explanation of these options, see the DMD manual.

• Thermostat (string): Specifies the thermostat to use. Only known thermostat available is the ANDERSON thermostat.
• Initial Temperature (double): initial temperature for the DMD simulation to use.
• Final Temperature (double): final temperature for the DMD simulation to achieve. If this is different than Initial Temperature, then the system temperature will slowly be changed to this final value.
• HEAT_X_C (double): specifies the rate at which energy is absorbed into the system from the bath. Essentially heat exchange rate.
• Echo File (string): specifies the name of the echo file.
• Movie File (string): specifies the name of the movie file.
• Restart File (string): specifies the name of the restart file.
• dt (positive int): Number of DMD time steps before saving in the echo, movie, and restart file.
• Time (double): Number of total DMD time steps to perform.

Titratable DMD options

The titr (dictionary) section specifies parameters for titratable-DMD. Note that this section is experimental as the method is not currently published. The following are the available options:

• titr on (bool): specifies if titrations are allowed to occur.
• dt (positive int): specifies the number of DMD time steps between evaluations of the protonation states.
• Buried Cutoff (double between 0 and 1/string): cutoff for the propka buried % above which a residue is considered solvent inaccessible. This only affects protonation state assessment and not the forcefield.
• Partner Distance (double): distance between which 2 residues can freely exhcange protons in an equilibrated manner in the timeframe of DMD steps. Typically 3.5 is used, consistend with the maximum DMD hydrogen bonding distance potential.
• Fixed States (list of strings): residues for which the protonation states should not change at all. Each element of the list should identify the residue using the standard notation in phd3 ("chain:resNum"). That is, if residue 43, and 54 of chain a should not have their protonation states changed, then the following is valid "Fixed States": ["A:43", ""A:54].

Protein Constraints

• Freeze Non-Residues (bool): specifies if any non-canonical amino acid residues (ie substrate and metals) should be frozen during the course of the simulation.
• Restrict Metal Ligands (bool): specifies if any atom bonded to a metal (heteroatom within 3.5 Ang) should be frozen. Additionally, dampenings between these atoms and the atoms it are bound to are placed to soften the impact of not allowing DMD to move these atoms.
• Frozen atoms (dictionary): specifies chains, residues, and atoms, for which they should be frozen during the course of the simulation. The following are valid options:
  • Chains (list of strings): specifies chains to freeze. All residues within the chain will be frozen. Chains are specified as "chainLetter".
  • Residues (list of strings): specifies residues to freeze. All atoms within these residues will be frozen. Residues are specified as "chainLetter:residueNumber".
  • Atoms (list of strings): specifies atoms to freeze. Atoms are specified as "chainLetter:residueNumber:atomID".
• Restrict Displacement (list of (string, string, double)): specifies constraints between two atoms in the protein such that DMD does not change the distance between them by the value specified per DMD time step. For example, for 2 cysteine residues, A:45 and B:65, which have a disulfide bridge, it is necessary to specify "Restrict Displacement": [["A:45:SG", "B:65:SG", 0.02]] as DMD does not understand disulfide bridges. This will prevent the atoms from moving too far apart during the course of the simulation.