Abstract:The structural dynamics of biological macromolecules, such as proteins, DNA/RNA, or complexes thereof, are strongly influenced by protonation changes of their typically many titratable groups, which explains their sensitivity to pH changes. Conversely, conformational and environmental changes of the biomolecule affect the protonation state of these groups. With few exceptions, conventional force field-based molecular dynamics (MD) simulations do not account for these effects, nor do they allow for coupling to a pH buffer. Here we present a GROMACS implementation of a rigorous Hamiltonian interpolation $\lambda$-dynamics constant pH method, which rests on GPU-accelerated Fast Multipole Method (FMM) electrostatics. Our implementation supports both CHARMM36m and Amber99sb*-ILDN force fields and is largely automated to enable seamless switching from regular MD to constant pH MD, involving minimal changes to the input files. Here, the first of two companion papers describes the underlying constant pH protocol and sample applications to several prototypical benchmark systems such as cardiotoxin V, lysozyme, and staphylococcal nuclease. Enhanced convergence is achieved through a new dynamic barrier height optimization method, and high p$K_a$ accuracy is demonstrated. We use Functional Mode Analysis and Mutual Information to explore the complex intra- and intermolecular couplings between the protonation states of titratable groups as well as those between protonation states and conformational dynamics. We identify striking conformation-dependent p$K_a$ variations and unexpected inter-residue couplings. Conformation-protonation coupling is identified as a primary cause of the slow protonation convergence notorious to constant pH simulations involving multiple titratable groups, suggesting enhanced sampling methods to accelerate convergence.

Biomolecular simulations at constant pH

Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

On the Validity of Constant pH Simulations

Constant pH Simulation with FMM Electrostatics in GROMACS. (A) Design and Applications

Constant-pH Molecular Dynamics Simulations of Closed and Open States of a Proton-gated Ion Channel

Development of constant-pH simulation methods in implicit solvent and applications in biomolecular systems

Introducing titratable water to all-atom molecular dynamics at constant pH.

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

Poisson–Boltzmann continuum-solvation models: applications to pH-dependent properties of biomolecules

Constant pH Coarse-Grained Molecular Dynamics with Stochastic Charge Neutralization

On Validity of Constant pH Simulations

A new method for reactive constant pH simulations

All-Atom Continuous Constant Ph Molecular Dynamics with Particle Mesh Ewald and Titratable Water.

An equilibrium constant pH molecular dynamics method for accurate prediction of pH-dependence in protein systems: Theory and application

Best Practices in Constant pH MD Simulations: Accuracy and Sampling

Scalable Constant pH Molecular Dynamics in GROMACS

Constant-pH Simulations with the Polarizable Atomic Multipole AMOEBA Force Field

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

Mimicking titration experiments with MD simulations: A protocol for the investigation of pH-dependent effects on proteins

Effects Of System Net Charge And Electrostatic Truncation On All-Atom Constant Ph Molecular Dynamics

Biomolecular conformational changes and ligand binding: from kinetics to thermodynamics.