Mimicking Titration Experiments with MD Simulations: A Protocol for the Investigation of PH-dependent Effects on Proteins

Overview

Journal Sci Rep

Specialty Science

Date 2016 Mar 4

PMID 26936826

Citations 8

Authors

Eileen Socher

Heinrich Sticht

Affiliations

Soon will be listed here.

Abstract

Protein structure and function are highly dependent on the environmental pH. However, the temporal or spatial resolution of experimental approaches hampers direct observation of pH-induced conformational changes at the atomic level. Molecular dynamics (MD) simulation strategies (e.g. constant pH MD) have been developed to bridge this gap. However, one frequent problem is the sampling of unrealistic conformations, which may also lead to poor pKa predictions. To address this problem, we have developed and benchmarked the pH-titration MD (pHtMD) approach, which is inspired by wet-lab titration experiments. We give several examples how the pHtMD protocol can be applied for pKa calculation including peptide systems, Staphylococcus nuclease (SNase), and the chaperone HdeA. For HdeA, pHtMD is also capable of monitoring pH-dependent dimer dissociation in accordance with experiments. We conclude that pHtMD represents a versatile tool for pKa value calculation and simulation of pH-dependent effects in proteins.

Citing Articles

Cysteine Oxidation in Proteins: Structure, Biophysics, and Simulation.

Ruiz D, Sandoval-Perez A, Rangarajan A, Gunderson E, Jacobson M Biochemistry. 2022; 61(20):2165-2176.

PMID: 36161872 PMC: 9583617. DOI: 10.1021/acs.biochem.2c00349.

LAWS: Local alignment for water sites-Tracking ordered water in simulations.

Klyshko E, Kim J, Rauscher S Biophys J. 2022; 122(14):2871-2883.

PMID: 36116009 PMC: 10397812. DOI: 10.1016/j.bpj.2022.09.012.

Large-Scale Ligand Perturbations of the Protein Conformational Landscape Reveal State-Specific Interaction Hotspots.

Stachowski T, Fischer M J Med Chem. 2022; 65(20):13692-13704.

PMID: 35970514 PMC: 9619398. DOI: 10.1021/acs.jmedchem.2c00708.

A study on the secondary structure of the metalloregulatory protein CueR: effect of pH, metal ions and DNA.

Balogh R, Nemeth E, Jones N, Vronning Hoffmann S, Jancso A, Gyurcsik B Eur Biophys J. 2021; 50(3-4):491-500.

PMID: 33907862 DOI: 10.1007/s00249-021-01539-z.

Identification of a unique Ca-binding site in rat acid-sensing ion channel 3.

Zuo Z, Smith R, Chen Z, Agharkar A, Snell H, Huang R Nat Commun. 2018; 9(1):2082.

PMID: 29802295 PMC: 5970173. DOI: 10.1038/s41467-018-04424-0.

References

Fitch C, Karp D, Lee K, Stites W, Lattman E, Garcia-Moreno E B . Experimental pK(a) values of buried residues: analysis with continuum methods and role of water penetration. Biophys J. 2002; 82(6):3289-304. PMC: 1302117. DOI: 10.1016/s0006-3495(02)75670-1. View

Khandogin J, Brooks 3rd C . Constant pH molecular dynamics with proton tautomerism. Biophys J. 2005; 89(1):141-57. PMC: 1366513. DOI: 10.1529/biophysj.105.061341. View

Swails J, York D, Roitberg A . Constant pH Replica Exchange Molecular Dynamics in Explicit Solvent Using Discrete Protonation States: Implementation, Testing, and Validation. J Chem Theory Comput. 2014; 10(3):1341-1352. PMC: 3985686. DOI: 10.1021/ct401042b. View

Berman H, Westbrook J, Feng Z, Gilliland G, Bhat T, Weissig H . The Protein Data Bank. Nucleic Acids Res. 1999; 28(1):235-42. PMC: 102472. DOI: 10.1093/nar/28.1.235. View

Baptista A, MARTEL P, Petersen S . Simulation of protein conformational freedom as a function of pH: constant-pH molecular dynamics using implicit titration. Proteins. 1997; 27(4):523-44. View

Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C . Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins. 2006; 65(3):712-25. PMC: 4805110. DOI: 10.1002/prot.21123. View

Mongan J, Case D . Biomolecular simulations at constant pH. Curr Opin Struct Biol. 2005; 15(2):157-63. DOI: 10.1016/j.sbi.2005.02.002. View

McNutt M, Mullins L, Raushel F, Pace C . Contribution of histidine residues to the conformational stability of ribonuclease T1 and mutant Glu-58----Ala. Biochemistry. 1990; 29(33):7572-6. DOI: 10.1021/bi00485a005. View

Tournaire-Roux C, Sutka M, Javot H, Gout E, Gerbeau P, Luu D . Cytosolic pH regulates root water transport during anoxic stress through gating of aquaporins. Nature. 2003; 425(6956):393-7. DOI: 10.1038/nature01853. View

10.

Yang F, Gustafson K, Boyd M, Wlodawer A . Crystal structure of Escherichia coli HdeA. Nat Struct Biol. 1998; 5(9):763-4. DOI: 10.1038/1796. View

11.

Vila-Vicosa D, Teixeira V, Santos H, Machuqueiro M . Conformational study of GSH and GSSG using constant-pH molecular dynamics simulations. J Phys Chem B. 2013; 117(25):7507-17. DOI: 10.1021/jp401066v. View

12.

Olsson M, Sondergaard C, Rostkowski M, Jensen J . PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. J Chem Theory Comput. 2015; 7(2):525-37. DOI: 10.1021/ct100578z. View

13.

Lee J, Miller B, Damjanovic A, Brooks B . Enhancing constant-pH simulation in explicit solvent with a two-dimensional replica exchange method. J Chem Theory Comput. 2015; 11(6):2560-74. DOI: 10.1021/ct501101f. View

14.

Williams S, Blachly P, McCammon J . Measuring the successes and deficiencies of constant pH molecular dynamics: a blind prediction study. Proteins. 2011; 79(12):3381-8. PMC: 3227005. DOI: 10.1002/prot.23136. View

15.

Lee M, Salsbury Jr F, Brooks 3rd C . Constant-pH molecular dynamics using continuous titration coordinates. Proteins. 2004; 56(4):738-52. DOI: 10.1002/prot.20128. View

16.

Donnini S, Tegeler F, Groenhof G, Grubmuller H . Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics. J Chem Theory Comput. 2011; 7(6):1962-1978. PMC: 3114466. DOI: 10.1021/ct200061r. View

17.

Garrison M, Crowhurst K . NMR-monitored titration of acid-stress bacterial chaperone HdeA reveals that Asp and Glu charge neutralization produces a loosened dimer structure in preparation for protein unfolding and chaperone activation. Protein Sci. 2013; 23(2):167-78. PMC: 3926742. DOI: 10.1002/pro.2402. View

18.

Roche S, Bressanelli S, Rey F, Gaudin Y . Crystal structure of the low-pH form of the vesicular stomatitis virus glycoprotein G. Science. 2006; 313(5784):187-91. DOI: 10.1126/science.1127683. View

19.

Bashford D, Case D, Dalvit C, Tennant L, Wright P . Electrostatic calculations of side-chain pK(a) values in myoglobin and comparison with NMR data for histidines. Biochemistry. 1993; 32(31):8045-56. DOI: 10.1021/bi00082a027. View

20.

Chu W, Zhang J, Zheng Q, Chen L, Wu Y, Xue Q . Constant pH molecular dynamics (CpHMD) and molecular docking studies of CquiOBP1 pH-induced ligand releasing mechanism. J Mol Model. 2012; 19(3):1301-9. DOI: 10.1007/s00894-012-1680-0. View