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.

On the Validity of Constant pH Simulations

On Validity of Constant pH Simulations

A new method for reactive constant pH simulations

Comment on "Simulations of ionization equilibria in weak polyelectrolyte solutions and gels" by J. Landsgesell, L. Nová, O. Rud, F. Uhlík, D. Sean, P. Hebbeker, C. Holm and P. Ko\v sovan, Soft Matter, 2019,15, 1155-1185

How to Correct Systematic Errors in Constant-pH Ensemble Simulations

How to Correct Erroneous Symmetry-Breaking in Coarse-Grained Constant-pH Simulations

Biomolecular simulations at constant pH

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

Constant pH Coarse-Grained Molecular Dynamics with Stochastic Charge Neutralization

Canonical titration simulations

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

Application of Explicit-Solvent Constant Ph Molecular Dynamics to Proteins

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

Simulation of weak polyelectrolytes: A comparison between the constant pH and the reaction ensemble method

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

Charge Regulation of Colloidal Particles: Theory and Simulations

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

Best Practices in Constant pH MD Simulations: Accuracy and Sampling

Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

Scalable Constant pH Molecular Dynamics in GROMACS