Is the prediction of pKa values by constant-pH molecular dynamics being hindered by inherited problems?
✍ Scribed by Miguel Machuqueiro; António M. Baptista
- Book ID
- 105358294
- Publisher
- John Wiley and Sons
- Year
- 2011
- Tongue
- English
- Weight
- 749 KB
- Volume
- 79
- Category
- Article
- ISSN
- 0887-3585
No coin nor oath required. For personal study only.
✦ Synopsis
Abstract
In this study, we investigate two factors that can hinder the performance of constant‐pH molecular dynamics methods in predicting protein p__K__~a~ values, using hen egg white lysozyme as a test system. The first factor is related to the molecular definition and p__K__~a~ value of model compounds in the Poisson‐Boltzmann framework. We address this by defining the model compound as a molecular fragment with an associated p__K__~a~ value that is calibrated against experimental data, which results in a decrease of 0.12 units in p__K__~a~ errors. The second addressed factor is the possibility that detrimental structural distortions are being introduced in the simulations by the underlying molecular mechanics force field. This issue is investigated by analyzing how the gradual structural rearrangements affect the predicted p__K__~a~ values. The two GROMOS force fields studied here (43A1 and 53A6) yield good p__K__~a~ predictions, although a time‐dependent performance is observed: 43A1 performs better after a few nanoseconds of structural reorganization (p__K__~a~ errors of ∼0.45), while 53A6 gives the best prediction right at the first nanosecond (p__K__~a~ errors of 0.42). These results suggest that the good performance of constant‐pH molecular dynamics methods could be further improved if these force field limitations were overcome. Proteins 2011;. © 2011 Wiley‐Liss, Inc.