Related Experiment Video
Updated: Jun 16, 2026

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package
Published on: September 17, 2021
Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics.
Sarah L Williams1, César Augusto F de Oliveira, J Andrew McCammon
1Department of Chemistry & Biochemistry, University of California San Diego, La Jolla, California 92093-0365, Center for Theoretical Biological Physics, University of California San Diego, La Jolla, California 92093, Howard Hughes Medical Institute, University of California San Diego, La Jolla, California 92093-0365, Department of Pharmacology, University of California San Diego, La Jolla, California 92093-0365.
This study enhances the constant pH method by coupling it with accelerated molecular dynamics to improve sampling for molecular simulations. The new approach accelerates convergence and yields more accurate pK(a) values for complex protein systems like hen egg-white lysozyme (HEWL).
Area of Science:
- Computational Chemistry
- Biomolecular Modeling
- Molecular Dynamics
Background:
- Standard constant pH molecular dynamics methods face sampling challenges, leading to slow convergence for certain residues.
- Previous work by Mongan et al. showed promise but noted inconsistencies in predicted pK(a) values for systems like hen egg-white lysozyme (HEWL) due to slow convergence of strongly interacting residues.
Purpose of the Study:
- To propose an extension of the constant pH method by integrating it with accelerated molecular dynamics.
- To overcome sampling limitations inherent in standard constant pH molecular dynamics approaches.
- To improve the accuracy and speed of pK(a) predictions for complex biomolecular systems.
Main Methods:
- Coupling the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics.
- Applying the adapted method to the hen egg-white lysozyme (HEWL) system for pK(a) value prediction.
Main Results:
- The coupled method demonstrates improved sampling efficiency compared to the standard constant pH method.
- Faster convergence was observed for residues that previously converged slowly.
- Predicted pK(a) values for HEWL showed closer agreement with experimental data, particularly for challenging residues.
Conclusions:
- The integration of accelerated molecular dynamics with the constant pH method effectively addresses sampling issues in molecular simulations.
- This enhanced approach provides more reliable and experimentally consistent pK(a) predictions for biomolecules.
- The method offers a valuable advancement for computational studies of protein protonation states.
Related Concept Videos
Spin–Spin Coupling Constant: Overview
Qualitatively, any spin plus-half nucleus polarizes the spins of its electrons to the minus-half state. Consequently, the paired electron in the hydrogen–carbon bond must have a...
Spin–Spin Coupling: Three-Bond Coupling (Vicinal Coupling)
The extent of coupling depends on the C‑C bond length, the two H‑C‑C angles, any electron-withdrawing substituents, and the dihedral angle between the involved orbitals. The...
Spin–Spin Coupling: One-Bond Coupling
¹H NMR: Long-Range Coupling
In alkenes, spin information is communicated via σ–π overlap, as seen in allylic (four-bond) and homoallylic (five-bond) couplings. These coupling interactions are stronger when the σ bond is parallel to the alkene π orbitals.
Double Resonance Techniques: Overview
Spin decoupling is usually achieved by...
Spin–Spin Coupling: Two-Bond Coupling (Geminal Coupling)
The central atom need not be NMR-active because its electrons are affected by the electron polarization of the spin-active atoms. However, spin information is transmitted less effectively than in one-bond coupling, and 2J values are usually weaker than 1J values. The energy of...

