Overcoming potential energy distortions in constrained internal coordinate molecular dynamics simulations

Journal Article · Thu Jan 28 04:00:00 EST 2016 · Journal of Chemical Physics

DOI:https://doi.org/10.1063/1.4939532· OSTI ID:22493685

Kandel, Saugat; Salomon-Ferrer, Romelia; Larsen, Adrien B.; Vaidehi, Nagarajan ^[1]; Jain, Abhinandan ^[2]

---

1. Division of Immunology, Beckman Research Institute of the City of Hope, Duarte, California 91010 (United States)
2. Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, California 91109 (United States)

The Internal Coordinate Molecular Dynamics (ICMD) method is an attractive molecular dynamics (MD) method for studying the dynamics of bonded systems such as proteins and polymers. It offers a simple venue for coarsening the dynamics model of a system at multiple hierarchical levels. For example, large scale protein dynamics can be studied using torsional dynamics, where large domains or helical structures can be treated as rigid bodies and the loops connecting them as flexible torsions. ICMD with such a dynamic model of the protein, combined with enhanced conformational sampling method such as temperature replica exchange, allows the sampling of large scale domain motion involving high energy barrier transitions. Once these large scale conformational transitions are sampled, all-torsion, or even all-atom, MD simulations can be carried out for the low energy conformations sampled via coarse grained ICMD to calculate the energetics of distinct conformations. Such hierarchical MD simulations can be carried out with standard all-atom forcefields without the need for compromising on the accuracy of the forces. Using constraints to treat bond lengths and bond angles as rigid can, however, distort the potential energy landscape of the system and reduce the number of dihedral transitions as well as conformational sampling. We present here a two-part solution to overcome such distortions of the potential energy landscape with ICMD models. To alleviate the intrinsic distortion that stems from the reduced phase space in torsional MD, we use the Fixman compensating potential. To additionally alleviate the extrinsic distortion that arises from the coupling between the dihedral angles and bond angles within a force field, we propose a hybrid ICMD method that allows the selective relaxing of bond angles. This hybrid ICMD method bridges the gap between all-atom MD and torsional MD. We demonstrate with examples that these methods together offer a solution to eliminate the potential energy distortions encountered in constrained ICMD simulations of peptide molecules.

• OSTI ID: 22493685
• Journal Information: Journal of Chemical Physics, Journal Name: Journal of Chemical Physics Journal Issue: 4 Vol. 144; ISSN JCPSA6; ISSN 0021-9606
• Country of Publication: United States
• Language: English

Related Subjects

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
ATOMS
BOND ANGLE
BOND LENGTHS
CONFORMATIONAL CHANGES
COUPLING
DIFFUSION BARRIERS
LIMITING VALUES
MOLECULAR DYNAMICS METHOD
MOLECULES
PEPTIDES
PHASE SPACE
POLYMERS
POTENTIAL ENERGY
POTENTIALS
TORSION