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.

Recent Development and Application of Constant Ph Molecular Dynamics.

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

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

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

Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

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

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

Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

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

Best Practices in Constant pH MD Simulations: Accuracy and Sampling

Constant pH Coarse-Grained Molecular Dynamics with Stochastic Charge Neutralization

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

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

Conformational changes of enzymes and DNA in molecular dynamics: influenced by pH, temperature, and ligand.

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

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

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

Force Field Limitations of All-Atom Continuous Constant pH Molecular Dynamics

Constant pH Replica Exchange Molecular Dynamics in Biomolecules Using a Discrete Protonation Model

Biomolecular simulations at constant pH