Molecular Dynamics/Free Energy Perturbation Studies of the Thermostable V74I Mutant of Ribonuclease HI

Abstract

Recent site-directed mutagenesis and thermodynamic studies have shown that the V74I mutant of Escherichia coli ribonuclease HI (RNase HI) is more stable than the wild type protein [Ishikawa et al., Biochemistry 32, 6171 (1993)]. In order to clarify the stabilization mechanism of this mutant, we calculated the free energy change due to the mutation Val 74→Ile in both the native and denatured states by free energy perturbations based on molecular dynamics (MD) simulations. We carried out inclusive MD simulations for the protein in water; i.e., fully solvated, no artificial constraints applied, and all long-range Coulomb interactions included. We found that the free energy of the mutant increased slightly relative to the wild type, in the native state by 1.60 kcal/mol, and in the denatured state by 2.25 kcal/mol. The unfolding free energy increment of the mutant (0.66 ± 0.19 kcal/mol) was in good agreement with the experimental value (0.6 kcal/mol). The hysteresis error in the free energy calculations, i.e., forward and reverse perturbations, was only ±0.19 kcal/mol. These results show that the V74I mutant is stabilized relative to the wild type by the increased free energy of the denatured state and not by a decrease in the free energy of the native state as had been proposed earlier based on the mutant X-ray structure. It was found that the stabilization was caused by a loss of solvation energy in the mutant denatured state and not by improved packing interactions inside the native protein.

Key Words:

• Cavity-filling mutation
• free energy perturbation method
• thermal stability
• particle-particle and particle-cell (PPPC) method
• acceptance ratio method (ARM)