Our overall goal is to develop the computational infrastructure to enable pH to be treated the same level as temperature and pressure in molecular dynamics simulation studies of biological molecules. These tools and methods need to be robust and straightforward to employ. Our work to date in this area has had a significant impact on the field; others have implemented similar methods, and key applications demonstrate the utility of direct incorporation of pH effects into molecular simulations.
Our specific aims are directed toward the establishment of a robust methodology for explicit solvent constant pH simulations and are driven by applications to nucleic acids and to a specific protein system where we will explore pH sensing elements in the chaperon activity of a conditionally disordered protein.
Aim 1 is to develop novel and robust explicit solvent CpHMD methods for nucleic acids. Nucleic acids present special challenges to implicit solvent models because of strong electrostatic fields from the nucleic acid polyanion and the significant influences of their ion atmosphere. Explicit solvent, with explicit ions, provides an alternative to the implicit solvet models. New developments in CpHMD will extend the approach to use an explicit solvent representation and will improve the precision, accuracy and robustness of the methodology. Our developments will be benchmarked against RNA and DNA measured pKa values for key bases, predominately the bases A and C, and through explicit collaborations with experimental colleagues.
Aim 2 will apply CpHMD to explore the role of pH in mediating DNA/RNA conformation and base-base interactions. We will examine a number of specific cases, pH-mediated conformational switching in the SARS-CoV coronavirus, in the U6 intra-molecular stem loop of the splicosome, and the role of pH in the catalytic mechanism of the hairpin ribozyme. These applications address fundamental questions regarding the role of pH-driven protonation equilibrium in structure, dynamics and function in nucleic acids. Collaborations with the Al-Hashimi and Walter groups at Michigan and with Victoria D'Souza from Harvard provide both experimental benchmarks for our calculations and enable additional experiment and analysis to be informed from the calculations.
Aim 3 aims to apply CpHMD to delineate and redesign the pH sensing mechanism in the Bacterial pH stress response protein chaperon HdeA. With the Bardwell group, we have redesigned the pH response elements in HdeA to produce a constitutively active and intrinsically disordered protein at pH 7. This provides a basis for exploration of the interactions within the largely disordered, but chaperon active, protein and the interactions it makes with its substrates. Collaboration with the Bardwell and Al-Hashimi groups will combine, coarse-grained simulations with CpHMD and NMR to explore the nature of the active disordered ensemble and how its members interact with substrate molecules.

Public Health Relevance

The role that pH plays in mediating biological function, including protein folding and misfolding, viral infectivity, biochemical catalysis by proteins and nucleic acids, is ubiquitious and thereby has important consequences in understanding biological processes associated with human disease. This proposed work aims to establish a robust theoretical and computational framework to enable the computational investigation of pH mediated processes of biomedical interest, from the pH sensing and chaperone mechanism of enteric Bacterial proteins to the role of pH mediated RNA dynamics in the SAR-CoV coronavirus frameshifting.

National Institute of Health (NIH)
National Institute of General Medical Sciences (NIGMS)
Research Project (R01)
Project #
Application #
Study Section
Special Emphasis Panel (ZRG1)
Program Officer
Preusch, Peter
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of Michigan Ann Arbor
Schools of Arts and Sciences
Ann Arbor
United States
Zip Code
Vilseck, Jonah Z; Armacost, Kira A; Hayes, Ryan L et al. (2018) Predicting Binding Free Energies in a Large Combinatorial Chemical Space Using Multisite ? Dynamics. J Phys Chem Lett 9:3328-3332
Ding, Xinqiang; Hayes, Ryan L; Vilseck, Jonah Z et al. (2018) CDOCKER and ?-dynamics for prospective prediction in D?R Grand Challenge 2. J Comput Aided Mol Des 32:89-102
Ding, Xinqiang; Vilseck, Jonah Z; Hayes, Ryan L et al. (2017) Gibbs Sampler-Based ?-Dynamics and Rao-Blackwell Estimator for Alchemical Free Energy Calculation. J Chem Theory Comput 13:2501-2510
Arthur, Evan J; Brooks 3rd, Charles L (2016) Efficient implementation of constant pH molecular dynamics on modern graphics processors. J Comput Chem 37:2171-80
Horowitz, Scott; Salmon, Loïc; Koldewey, Philipp et al. (2016) Visualizing chaperone-assisted protein folding. Nat Struct Mol Biol 23:691-7
Dickson, Alex; Ahlstrom, Logan S; Brooks 3rd, Charles L (2016) Coupled folding and binding with 2D Window-Exchange Umbrella Sampling. J Comput Chem 37:587-94
Laricheva, Elena N; Goh, Garrett B; Dickson, Alex et al. (2015) pH-dependent transient conformational states control optical properties in cyan fluorescent protein. J Am Chem Soc 137:2892-900
Frank, Aaron T; Law, Sean M; Ahlstrom, Logan S et al. (2015) Predicting protein backbone chemical shifts from C? coordinates: extracting high resolution experimental observables from low resolution models. J Chem Theory Comput 11:325-31
Panahi, Afra; Brooks 3rd, Charles L (2015) Membrane environment modulates the pKa values of transmembrane helices. J Phys Chem B 119:4601-7
Zeng, Xiancheng; Mukhopadhyay, Suchetana; Brooks 3rd, Charles L (2015) Residue-level resolution of alphavirus envelope protein interactions in pH-dependent fusion. Proc Natl Acad Sci U S A 112:2034-9

Showing the most recent 10 out of 16 publications