Langevin dynamics of proteins at constant pH
Abstract
An application of the Langevin dynamics algorithm for simulation of protein conformational equilibria at constant pH is presented. The algorithm is used to compute average protonation of titratable groups in ovomucoid third domain, as functions of pH, resulting in data, basically equivalent to the pH dependencies of chemical shifts obtained from multidimensional nuclear magnetic resonance (NMR) spectroscopy, for the protein titratable residues. The pKa values obtained from the simulation are in reasonable agreement with experimental data. Possible improvements of this methodology, using achievements from other fields of mesoscopic biomolecular simulations, are also discussed.
- Publication:
-
Physical Review E
- Pub Date:
- November 2002
- DOI:
- 10.1103/PhysRevE.66.051911
- Bibcode:
- 2002PhRvE..66e1911W
- Keywords:
-
- 87.15.Aa;
- Theory and modeling;
- computer simulation