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.

CHARMM-GUI PDB Reader and Manipulator: Covalent Ligand Modeling and Simulation

CHARMM-GUI PDB manipulator for advanced modeling and simulations of proteins containing nonstandard residues.

CHARMM-GUI ligand reader and modeler for CHARMM force field generation of small molecules

Computer Modeling for Interaction of Ligands and Receptors and Their Applications in Drug Design

GPCR-ModSim: A comprehensive web based solution for modeling G-protein coupled receptors

CovalentDock: Automated Covalent Docking with Parameterized Covalent Linkage Energy Estimation and Molecular Geometry Constraints.

PDB-CAT: A User-Friendly Tool to Classify and Analyze PDB Protein-Ligand Complexes

CHARMM-GUI PACE CG Builder for Solution, Micelle, and Bilayer Coarse-Grained Simulations

CHARMM-GUI QM/MM Interfacer for a Quantum Mechanical and Molecular Mechanical (QM/MM) Simulation Setup: 1. Semiempirical Methods

Expanding Automated Multiconformer Ligand Modeling to Macrocycles and Fragments

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

MolHyb: A Web Server for Structure-Based Drug Design by Molecular Hybridization

CHARMM: The biomolecular simulation program

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

A Practical Guide to Molecular Docking and Homology Modelling for Medicinal Chemists

MDFit: Automated molecular simulations workflow enables high throughput assessment of ligands-protein dynamics

Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

@TOME 3.0: Interfacing Protein Structure Modeling and Ligand Docking

Molcontroller: a VMD Graphical User Interface for Manipulating Molecules

COACH-D: improved protein-ligand binding sites prediction with refined ligand-binding poses through molecular docking