Localized thermodynamic coupling between hydrogen bonding and microenvironment polarity substantially stabilizes proteins

The energetic contributions of hydrogen bonding to protein folding are still unclear, despite over 70 years of study. This is due partly to the difficulty of extracting thermodynamic information about specific interactions from protein mutagenesis data, and partly to the context dependence of hydrogen bond strengths. Herein, we test the hypothesis that hydrogen bond strengths depend on the polarity of their microenvironment, with stronger hydrogen bonds forming in non-polar surroundings. Double mutant thermodynamic cycle analysis using a combination of amide-to-ester backbone mutagenesis and traditional side chain mutagenesis revealed that hydrogen bonds can be stronger by up to 1.2 kcal mol−1 when they are sequestered in hydrophobic surroundings than when they are solvent exposed. Such large coupling energies between hydrogen bond strengths and local polarity suggest that the context dependence of hydrogen bond strengths must be accounted for in any comprehensive account of the forces responsible for protein folding