Journal
PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS
Volume 79, Issue 12, Pages 3437-3447Publisher
WILEY
DOI: 10.1002/prot.23115
Keywords
pentapeptides; hen egg white lysozyme; simulation; Poisson-Boltzmann; Monte Carlo; GROMOS; 43A1; 53A6
Categories
Funding
- Fundacao para a Ciencia e Tecnologia, Portugal [SFRH/BPD/29358/2006, PTDC/BIA-PRO/104378/2008]
- Fundação para a Ciência e a Tecnologia [SFRH/BPD/29358/2006, PTDC/BIA-PRO/104378/2008] Funding Source: FCT
Ask authors/readers for more resources
In this study, we investigate two factors that can hinder the performance of constant-pH molecular dynamics methods in predicting protein pK(a) values, using hen egg white lysozyme as a test system. The first factor is related to the molecular definition and pK(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 pK(a) value that is calibrated against experimental data, which results in a decrease of 0.12 units in pK(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 pK(a) values. The two GROMOS force fields studied here (43A1 and 53A6) yield good pK(a) predictions, although a time-dependent performance is observed: 43A1 performs better after a few nanoseconds of structural reorganization (pK(a) errors of similar to 0.45), while 53A6 gives the best prediction right at the first nanosecond (pK(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; 79:3437-3447. (C) 2011 Wiley-Liss, Inc.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available