Improved protein free energy calculation by more accurate treatment of nonbonded energy: application to chymotrypsin inhibitor 2, V57A.

We developed a software package for improved free energy calculation, in which spherical solvent boundary potential, cell multipole method, and Nosé-Hoover equation are employed. The performance of the developed software package is demonstrated in the case of valine to alanine mutation of the 57th residue in chymotrypsin inhibitor 2. By using this package, we obtained results quantitatively comparable to experimental results. By the free energy component analysis, it is shown that leucine 51, arginine 65, arginine 67, and phenylalanine 69 residues contribute significantly to the total free energy shift, deltadeltaG. Among them, contribution from the hydrophilic arginine 67 residue, which is in close contact with the mutation site, is the largest. Structure around the mutation site is largely changed by the mutation. The structure change is caused mainly by two effects, hydrophobic interaction and short-range interaction along the sequence. Effects of Nosé-Hoover algorithm and Kirkwood reaction field are also discussed.