Calculation of protein ionization equilibria with conformational sampling: pK<sub>a</sub> of a model leucine zipper, GCN4 and barnase

TitleCalculation of protein ionization equilibria with conformational sampling: pKa of a model leucine zipper, GCN4 and barnase
Publication TypeJournal Article
Year of Publication2002
AuthorsGorfe A.A, Ferrara P., Caflisch A., Marti D.N, Bosshard H.R, Jelesarov I.
JournalProteins: Structure, Function, and Bioinformatics
Date Published2002 Jan 1
Type of ArticleResearch Article
KeywordsComputer Simulation, Crystallography, X-Ray, DNA-Binding Proteins, Electrochemistry, Fungal Proteins, Leucine Zippers, Magnetic Resonance Imaging, Pliability, Protein Conformation, Protein Kinases, Ribonucleases, Saccharomyces cerevisiae Proteins

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.



Alternate JournalProteins
PubMed ID11746702
Full Text PDF: