Ab initio quantum mechanical calculation of protein in solution is carried out to generate polarized protein-specific charge(s) (PPC) for molecular dynamics (MD) stimulation of protein. The quantum calculation of protein is made possible by developing a fragment-based quantum chemistry approach in combination with the implicit continuum solvent model. The computed electron density of protein is utilized to derive PPCs that represent the polarized electrostatic state of protein near the native structure. These PPCs are atom-centered like those in the standard force fields and are thus computationally attractive for molecular dynamics simulation of protein. Extensive MD simulations have been carried out to investigate the effect of electronic polarization on the structure and dynamics of thioredoxin. Our study shows that the dynamics of thioredoxin is stabilized by electronic polarization through explicit comparison between MD results using PPC and AMBER charges. In particular, MD free-energy calculation using PPCs accurately reproduced the experimental value of pK(a) shift for ionizable residue Asp(26) buried inside thioredoxin, whereas previous calculations using standard force fields overestimated pK(a) shift by twice as much. Accurate prediction of pK(a) shifts by rigorous MD free energy simulation for ionizable residues buried inside protein has been a significant challenge in computational biology for decades. This study presented strong evidence that electronic polarization of protein plays an important role in protein dynamics.

A new method that incorporates the conductorlike polarizable continuum model (CPCM) with the recently developed molecular fractionation with conjugate caps (MFCC) approach is developed for ab initio calculation of electrostatic solvation energy of protein. The application of the MFCC method makes it practical to apply CPCM to calculate electrostatic solvation energy of protein or other macromolecules in solution. In this MFCC-CPCM method, calculation of protein solvation is divided into calculations of individual solvation energies of fragments (residues) embedded in a common cavity defined with respect to the entire protein. Besides computational efficiency, the current approach also provides additional information about contribution to protein solvation from specific fragments. Numerical studies are carried out to calculate solvation energies for a variety of peptides including α helices and β sheets. Excellent agreement between the MFCC-CPCM result and those from the standard full system CPCM calculation is obtained. Finally, the MFCC-CPCM calculation is applied to several real proteins and the results are compared to classical molecular mechanics Poisson-Boltzmann (MM/PB) and quantum Divid-and-Conque Poisson-Boltzmann (D&C-PB) calculations. Large wave function distortion energy (solute polarization energy) is obtained from the quantum calculation which is missing in the classical calculation. The present study demonstrates that the MFCC-CPCM method is readily applicable to studying solvation of proteins.