MCCE2 : improving protein pKa calculations with extensive side chain rotamer sampling

2009 Wiley Periodicals, Inc.

Bibliographische Detailangaben
Veröffentlicht in:Journal of computational chemistry. - 1984. - 30(2009), 14 vom: 15. Nov., Seite 2231-47
1. Verfasser: Song, Yifan (VerfasserIn)
Weitere Verfasser: Mao, Junjun, Gunner, M R
Format: Online-Aufsatz
Sprache:English
Veröffentlicht: 2009
Zugriff auf das übergeordnete Werk:Journal of computational chemistry
Schlagworte:Journal Article Research Support, N.I.H., Extramural Research Support, U.S. Gov't, Non-P.H.S. Ligands Proteins
LEADER 01000caa a22002652 4500
001 NLM186975031
003 DE-627
005 20240511232019.0
007 cr uuu---uuuuu
008 231223s2009 xx |||||o 00| ||eng c
024 7 |a 10.1002/jcc.21222  |2 doi 
028 5 2 |a pubmed24n1404.xml 
035 |a (DE-627)NLM186975031 
035 |a (NLM)19274707 
040 |a DE-627  |b ger  |c DE-627  |e rakwb 
041 |a eng 
100 1 |a Song, Yifan  |e verfasserin  |4 aut 
245 1 0 |a MCCE2  |b improving protein pKa calculations with extensive side chain rotamer sampling 
264 1 |c 2009 
336 |a Text  |b txt  |2 rdacontent 
337 |a ƒaComputermedien  |b c  |2 rdamedia 
338 |a ƒa Online-Ressource  |b cr  |2 rdacarrier 
500 |a Date Completed 19.01.2010 
500 |a Date Revised 11.05.2024 
500 |a published: Print 
500 |a Citation Status MEDLINE 
520 |a 2009 Wiley Periodicals, Inc. 
520 |a 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 
650 4 |a Journal Article 
650 4 |a Research Support, N.I.H., Extramural 
650 4 |a Research Support, U.S. Gov't, Non-P.H.S. 
650 7 |a Ligands  |2 NLM 
650 7 |a Proteins  |2 NLM 
700 1 |a Mao, Junjun  |e verfasserin  |4 aut 
700 1 |a Gunner, M R  |e verfasserin  |4 aut 
773 0 8 |i Enthalten in  |t Journal of computational chemistry  |d 1984  |g 30(2009), 14 vom: 15. Nov., Seite 2231-47  |w (DE-627)NLM098138448  |x 1096-987X  |7 nnns 
773 1 8 |g volume:30  |g year:2009  |g number:14  |g day:15  |g month:11  |g pages:2231-47 
856 4 0 |u http://dx.doi.org/10.1002/jcc.21222  |3 Volltext 
912 |a GBV_USEFLAG_A 
912 |a SYSFLAG_A 
912 |a GBV_NLM 
912 |a GBV_ILN_350 
951 |a AR 
952 |d 30  |j 2009  |e 14  |b 15  |c 11  |h 2231-47