Studying the Effect of Site-Specific Hydrophobicity and Polarization on Hydrogen Bond Energy of Protein Using a Polarizable Method.

Quantification of backbone hydrogen bond energies in protein folding has remained elusive despite extensive theoretical and experimental investigations over the past 70 years. This is due to difficulties in experimental mutagenesis study as well as the lack of quantitatively reliable methods in theoretical calculation. Recent advance in experiment has enabled accurate measurement of site-specific backbone hydrogen bond energy in protein. In the present work, we developed an accurate and practical polarizable method to study site-specific hydrogen bond energies in the PIN WW domain. Excellent quantitative agreement between our calculated hydrogen bonding energy and recent experimental measurement is obtained. The direct comparison between theory and experiment helps uncover the microscopic mechanism of experimentally observed context dependent hydrogen bond contribution to protein stability in beta-sheet. In particular, our study reveals two effects that act in a cooperative manner to impact the strength of a hydrogen bond. One is the dynamic stability of the hydrogen bond determined by nearby solvent molecules, and the other is the polarization state of the hydrogen bond influenced by local electrostatic environment. The polar character of the hydrogen bond results in strong coupling between hydrophobic and polarization interactions in a cooperative manner. This nonadditive character in hydrogen bonding should help us better understand the microscopic mechanism in protein folding. Our study also investigated the possible structural effect of backbone amide to ester mutation which should be helpful to experimentalists using this technique in mutagenesis study.