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.

A new method for reactive constant pH simulations

On Validity of Constant pH Simulations

Canonical titration simulations

On the Validity of Constant pH Simulations

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

Reactive Monte Carlo simulations for charge regulation of colloidal particles

Theory of Charge Regulation of Colloidal Particles in Electrolyte Solutions

Constant pH Coarse-Grained Molecular Dynamics with Stochastic Charge Neutralization

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

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

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

A Generalized Grand-Reaction Method for Modelling the Exchange of Weak (Polyprotic) Acids between a Solution and a Weak Polyelectrolyte Phase

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

Biomolecular simulations at constant pH

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

Charge regulation of colloidal particles in aqueous solutions

Titration in Canonical and Grand-Canonical Ensembles

How to Correct Systematic Errors in Constant-pH Ensemble Simulations

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

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

Charge Regulation of Colloidal Particles: Theory and Simulations