Calculation of protein ionization equilibria with conformational sampling: pK<SUB><FONT SIZE='-1'>a</FONT></SUB> of a model leucine zipper, GCN4 and barnase

@article{citeulike:2983440, abstract = {The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pKa calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pKa estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pKa calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pKa predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pKa predictions for the model leucine zipper was reduced to as low as approximately 0.25 pKa units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pKa calculations, are discussed. Proteins 2002;46:41-60. {\copyright} 2001 Wiley-Liss, Inc.}, address = {Department of Biochemistry, University of Zurich, Zurich, Switzerland}, author = {Gorfe, Alemayehu A. and Ferrara, Philippe and Caflisch, Amedeo and Marti, Daniel N. and Bosshard, Hans R. and Jelesarov, Ilian }, citeulike-article-id = {2983440}, doi = {10.1002/prot.10027}, journal = {Proteins: Structure, Function, and Genetics}, keywords = {pka}, number = {1}, pages = {41--60}, posted-at = {2008-07-10 02:55:12}, priority = {2}, title = {Calculation of protein ionization equilibria with conformational sampling: pK<SUB><FONT SIZE='-1'>a</FONT></SUB> of a model leucine zipper, GCN4 and barnase}, url = {http://dx.doi.org/10.1002/prot.10027}, volume = {46}, year = {2002} }

TY - JOUR ID - citeulike:2983440 L3 - citeulike-article-id:2983440 TI - Calculation of protein ionization equilibria with conformational sampling: pK<SUB><FONT SIZE='-1'>a</FONT></SUB> of a model leucine zipper, GCN4 and barnase JF - Proteins: Structure, Function, and Genetics VL - 46 IS - 1 SP - 41 EP - 60 AD - Department of Biochemistry, University of Zurich, Zurich, Switzerland N2 - The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pKa calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pKa estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pKa calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pKa predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pKa predictions for the model leucine zipper was reduced to as low as approximately 0.25 pKa units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pKa calculations, are discussed. Proteins 2002;46:41-60. © 2001 Wiley-Liss, Inc. KW - pka AU - Gorfe, Alemayehu AU - Ferrara, Philippe AU - Caflisch, Amedeo AU - Marti, Daniel AU - Bosshard, Hans AU - Jelesarov, Ilian PY - 2002/// UR - http://dx.doi.org/10.1002/prot.10027 ER -