MCCE2: Improving protein pKa calculations with extensive side chain rotamer sampling
- 24 August 2009
- journal article
- research article
- Published by Wiley in Journal of Computational Chemistry
- Vol. 30 (14), 2231-2247
- https://doi.org/10.1002/jcc.21222
Abstract
Multiconformation continuum electrostatics (MCCE) explores different conformational degrees of freedom in Monte Carlo calculations of protein residue and ligand pKas. Explicit changes in side chain conformations throughout a titration create a position dependent, heterogeneous dielectric response giving a more accurate picture of coupled ionization and position changes. The MCCE2 methods for choosing a group of input heavy atom and proton positions are described. The pKas calculated with different isosteric conformers, heavy atom rotamers and proton positions, with different degrees of optimization are tested against a curated group of 305 experimental pKas in 33 proteins. QUICK calculations, with rotation around Asn and Gln termini, sampling His tautomers and torsion minimum hydroxyls yield an RMSD of 1.34 with 84% of the errors being Kas in the membrane protein bacteriorhodopsin. The inclusion of extra side chain positions distorts the dielectric boundary and also biases the calculated pKas by creating more neutral than ionized conformers. Methods for correcting these errors are introduced. Calculations are compared with multiple X-ray and NMR derived structures in 36 soluble proteins. Calculations with X-ray structures give significantly better pKas. Results with the default protein dielectric constant of 4 are as good as those using a value of 8. The MCCE2 program can be downloaded from http://www.sci.ccny.cuny.edu/∼mcce. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009Keywords
This publication has 114 references indexed in Scilit:
- Analysis of the electrochemistry of hemes with Ems spanning 800 mVProteins-Structure Function and Bioinformatics, 2008
- Electrostatic basis for the unidirectionality of the primary proton transfer in cytochrome c oxidaseProceedings of the National Academy of Sciences of the United States of America, 2008
- Factors influencing the energetics of electron and proton transfers in proteins. What can be learned from calculationsBiochimica et Biophysica Acta (BBA) - Bioenergetics, 2006
- Calculated Proton Uptake on Anaerobic Reduction of CytochromecOxidase: Is the Reaction Electroneutral?Biochemistry, 2006
- Electrostatic Environment of Hemes in Proteins: pKas of Hydroxyl LigandsBiochemistry, 2006
- Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsinJournal of Molecular Biology, 1992
- On the calculation of electrostatic interactions in proteinsJournal of Molecular Biology, 1985
- CHARMM: A program for macromolecular energy, minimization, and dynamics calculationsJournal of Computational Chemistry, 1983
- Calculation of the electric potential in the active site cleft due to α-helix dipolesJournal of Molecular Biology, 1982
- Refinement of protein conformations using a macromolecular energy minimization procedureJournal of Molecular Biology, 1969