Abstract:Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

On the Transferability of CHARMM36m Protein Force Field with LJ-PME: Hydrogen Bonding Dynamics under Elevated Pressures

Validating the CHARMM36m protein force field with LJ-PME reveals altered hydrogen bonding dynamics under elevated pressures

The Melting Curves and Entropy of Iron under High Pressure

Ab Initioand Molecular Dynamics Studies of Solid Β-Hmx: Effects of Hydrostatic Pressure and High Temperature

Structure Phase Transition from Α to Ε in Fe under Non-Hydrostatic Pressure: an Ab Initio Study

Properties and Phase Transitions of the Solid Β-HMX: Different Force Fields

A Transferrable Range-Separated Force Field for Water: Combining the Power of Both Physically-Motivated Models and Machine Learning Techniques

A Tale of Two Desolvation Potentials: An Investigation of Protein Behavior Under High Hydrostatic Pressure

A Force Balanced Fragmentation Method for Ab Initio Molecular Dynamic Simulation of Protein

Pressure Dependence of the Dynamic Crossover Temperatures in Protein and its Hydration Water

A Transferrable Coarse-Grained Force Field for Simulations of Polyethers and Polyether Blends

Achieving Pressure Consistency in Mechanochemical Simulations of Chemical Reactions Under Pressure

Quantum Fragment Based ab Initio Molecular Dynamics for Proteins.

MACE-OFF23: Transferable Machine Learning Force Fields for Organic Molecules

Pressure-Induced Hydrogen Transfer in 2-Butyne via a Double CH···π Aromatic Transition State

Hybrid Particle-Field Molecular Dynamics Under Constant Pressure

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

A theoretical study on behaviours of H2 on Fe(100) surface under environmental pressure and temperature

Likely Overstabilization of Charge-Charge Interactions in CHARMM36m(w): A Case for a99SB-disp Water

A Comparative Study Using State-of-the-art Electronic Structure Theories on Solid Hydrogen Phases under High Pressures

Building Machine Learning Force Fields of Proteins with Fragment-Based Approach and Transfer Learning