Quantifying the stabilizing energy of the intraprotein hydrogen bond due to local mutation

MD simulation of the WW domain of PIN based on a dynamically adjusted polarized protein-specific force field from quantum fragment calculations is carried out in both wild and VAL22ALA mutant states. The result shows that the geometry of the Arg14-TYR23 hydrogen bond is conserved upon mutation of VAL22 to ALA. However, the electrostatic energy of this hydrogen bond in the mutant is found to be 0.6 kcal/mol weaker than in the wild state, in close agreement with the experimentally measured upper limit of 1.2 kcal/mol. Analysis shows that the weakened energy of this hydrogen bond in the mutant is due to its dynamically changed polarization resulting from an altered local electrostatic environment near the hydrogen bond which becomes more exposed to the solvent than in the wild.