MCCE2: improving protein pKa calculations with extensive side chain rotamer sampling - PubMed (original) (raw)
MCCE2: improving protein pKa calculations with extensive side chain rotamer sampling
Yifan Song et al. J Comput Chem. 2009.
Abstract
Multiconformation continuum electrostatics (MCCE) explores different conformational degrees of freedom in Monte Carlo calculations of protein residue and ligand pK(a)s. 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 pK(a)s 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 pK(a)s 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 <1.5 pH units. FULL calculations adding heavy atom rotamers and side chain optimization yield an RMSD of 0.90 with 90% of the errors <1.5 pH unit. Good results are also found for pK(a)s in the membrane protein bacteriorhodopsin. The inclusion of extra side chain positions distorts the dielectric boundary and also biases the calculated pK(a)s 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 pK(a)s. 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.
Figures
Figure 1
Fragment of lysozyme structure 4LZT. (a) Single conformation dielectric boundary used to calculate the reaction field energy ΔΔ_G_rxn; and the reference pairwise interactions ΔGAi,BjES (bold lines) between the only conformer with partial charges (a conformer of Arg128 here) and the native conformer of all other residues; (b) Multiconformer dielectric boundary has more low dielectric boundary material so all pairwise interactions, ΔGAi,BjM have larger absolute values than ΔGAi,BjES. The raw pairwise interaction to each non-native conformer is corrected with eq. (6) to give ΔGAi:BjC.
Figure 2
Comparison of pairwise interactions of 1200 conformers in Barnase (1A2P chain A) at ε_p of 4. Δ_G_ES is calculated with only the interacting conformers present (Fig. 1a) whereas Δ_GM uses the standard multiconformation boundary conditions for calculating pairwise interactions (Fig. 1b, Table 1). ΔGAi,BjM versus ΔGAi,BjES for (a) 1613 charge–charge and (c) 9679 charge–dipole interactions and (e) 16641 dipole–dipole interactions. Lines show slope 1 and best-fit lines through the points. ΔGAi:BjC eq. (6) versus ΔGAi,BjES for (c) charge–charge and (c) charge–dipole interactions. Dipole–dipole interactions are generally small and no corrections are used. Line of slope 1 is shown.
Figure 3
Comparison of calculated p_K_a values using FULL MCCE conformer flexibility with experimentally measured values. The error bars represent the standard deviation of the values for different structures. The thick central line is the ideal where p_K_a (calc) = p_K_a (expt); the solid line bracket errors <1 pH unit and the dashed lines errors <2 pH units. Circled points highlight residues buried in the protein with desolvation energies >2.04 kcal/mol (1.5 pH units) or with p_K_as perturbed by >1.5 pH units from the solution value. (A) 305 averaged p_K_as obtained starting with 86 structures obtained by X-ray crystallography of 33 proteins; (B) 265 p_K_as obtained starting with 696 structures obtained by NMR methods of 24 proteins. The calculated and experimental p_K_as are provided in supporting information Table S2.
Figure 4
Shifts in calculated experimental p_K_as versus those calculated with a protein dielectric constant of 4 (●) or 8 (Δ). The dashed and dotted lines show errors of ±1 and ±2 pH units.
Similar articles
- Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins.
Georgescu RE, Alexov EG, Gunner MR. Georgescu RE, et al. Biophys J. 2002 Oct;83(4):1731-48. doi: 10.1016/S0006-3495(02)73940-4. Biophys J. 2002. PMID: 12324397 Free PMC article. - MCCE analysis of the pKas of introduced buried acids and bases in staphylococcal nuclease.
Gunner MR, Zhu X, Klein MC. Gunner MR, et al. Proteins. 2011 Dec;79(12):3306-19. doi: 10.1002/prot.23124. Epub 2011 Sep 9. Proteins. 2011. PMID: 21910138 - Exploring conformational changes coupled to ionization states using a hybrid Rosetta-MCCE protocol.
Song Y. Song Y. Proteins. 2011 Dec;79(12):3356-63. doi: 10.1002/prot.23146. Epub 2011 Aug 30. Proteins. 2011. PMID: 22072519 - Continuum Electrostatics Approaches to Calculating pKas and Ems in Proteins.
Gunner MR, Baker NA. Gunner MR, et al. Methods Enzymol. 2016;578:1-20. doi: 10.1016/bs.mie.2016.05.052. Epub 2016 Jun 20. Methods Enzymol. 2016. PMID: 27497160 Free PMC article. Review. - The pKa Cooperative: a collaborative effort to advance structure-based calculations of pKa values and electrostatic effects in proteins.
Nielsen JE, Gunner MR, García-Moreno BE. Nielsen JE, et al. Proteins. 2011 Dec;79(12):3249-59. doi: 10.1002/prot.23194. Epub 2011 Oct 15. Proteins. 2011. PMID: 22002877 Free PMC article. Review.
Cited by
- MFSD1 with its accessory subunit GLMP functions as a general dipeptide uniporter in lysosomes.
Jungnickel KEJ, Guelle O, Iguchi M, Dong W, Kotov V, Gabriel F, Debacker C, Dairou J, McCort-Tranchepain I, Laqtom NN, Chan SH, Ejima A, Sato K, Massa López D, Saftig P, Mehdipour AR, Abu-Remaileh M, Gasnier B, Löw C, Damme M. Jungnickel KEJ, et al. Nat Cell Biol. 2024 Jul;26(7):1047-1061. doi: 10.1038/s41556-024-01436-5. Epub 2024 Jun 5. Nat Cell Biol. 2024. PMID: 38839979 Free PMC article. - Titratable residues that drive RND efflux: Insights from molecular simulations.
Clark R, Newman KE, Khalid S. Clark R, et al. QRB Discov. 2024 Apr 1;5:e5. doi: 10.1017/qrd.2024.6. eCollection 2024. QRB Discov. 2024. PMID: 38689873 Free PMC article. Review. - Accurately Predicting Protein p_K_a Values Using Nonequilibrium Alchemy.
Wilson CJ, Karttunen M, de Groot BL, Gapsys V. Wilson CJ, et al. J Chem Theory Comput. 2023 Nov 14;19(21):7833-7845. doi: 10.1021/acs.jctc.3c00721. Epub 2023 Oct 11. J Chem Theory Comput. 2023. PMID: 37820376 Free PMC article. - Comparative Performance of High-Throughput Methods for Protein p_K_a Predictions.
Wei W, Hogues H, Sulea T. Wei W, et al. J Chem Inf Model. 2023 Aug 28;63(16):5169-5181. doi: 10.1021/acs.jcim.3c00165. Epub 2023 Aug 7. J Chem Inf Model. 2023. PMID: 37549424 Free PMC article. - PKAD-2: New entries and expansion of functionalities of the database of experimentally measured pKa's of proteins.
Ancona N, Bastola A, Alexov E. Ancona N, et al. J Comput Biophys Chem. 2023 Aug;22(5):515-524. doi: 10.1142/s2737416523500230. Epub 2023 Apr 25. J Comput Biophys Chem. 2023. PMID: 37520074 Free PMC article.
References
- Warshel A, Russell ST. Q Rev Biophys. 1984;17:283. - PubMed
- Honig B, Nicholls A. Science. 1995;268:1144. - PubMed
- Decoursey TE. Physiol Rev. 2003;83:475. - PubMed
- Garcia-Moreno EB, Fitch CA. Methods Enzymol. 2004;380:20. - PubMed
- Bashford D. Front Biosci. 2004;9:1082. - PubMed
Publication types
MeSH terms
Substances
LinkOut - more resources
Full Text Sources