Journal
PROTEINS-STRUCTURE FUNCTION AND BIOINFORMATICS
Volume 72, Issue 1, Pages 289-298Publisher
WILEY-BLACKWELL
DOI: 10.1002/prot.21923
Keywords
lysozyme; pK(a); simulation; Poisson-Boltzmann; Monte Carlo; tautomer; dielectric constant; stability
Categories
Ask authors/readers for more resources
In this work, we present the first application to a protein of the stochastic constant-pH molecular dynamics (MD) method with the inclusion of proton tautomerism. The acidic titration of HEWL was performed under different conditions. Both generalized reaction field (GRF) and particle mesh Ewald (PME) methods were used in the treatment of the long range electrostatics and, even though the PME simulations revealed to be more stable, the better results were obtained using GRF (pK(a) RMSD of 0.82 for GRF and 1.13 for PME). The results using PME at different dielectric constants (2, 4, and 8) also revealed that there was no significant improvement in pK(a)'s prediction upon increasing the dielectric constant. The secondary structure analysis of HEWL revealed a remarkably stable protein in the acidic pH range. The beta-sheet strands (unlike the alpha-helices) seem to be destabilized upon pH decrease, suggesting that the beta-domain is less stable than the alpha-domain. The four principal alpha-helices were also ordered according to their stability in the acidic pH range and the results (4 < 1 < 2 approximate to 3) were consistent with the ones obtained in thermal denaturation studies.
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