Solvent reaction field potential inside an uncharged globular protein: A bridge between implicit and explicit solvent models?

The solvent reaction field potential of an uncharged protein immersed in simple point charge/extended explicit solvent was computed over a series of molecular dynamics trajectories, in total $1560\mathrm{ns}$ of simulation time. A finite, positive potential of 13–24 $k_{b}T{e}_c^{-1}$ (where $T=300\mathrm{K}$ ), dependent on the geometry of the solvent-accessible surface, was observed inside the biomolecule. The primary contribution to this potential arose from a layer of positive charge density $1.0\mathrm{\AA }$ from the solute surface, on average $0.008e_c∕{\mathrm{\AA }}^3$ , which we found to be the product of a highly ordered first solvation shell. Significant second solvation shell effects, including additional layers of charge density and a slight decrease in the short-range solvent-solvent interaction strength, were also observed. The impact of these findings on implicit solvent models was assessed by running similar explicit solvent simulations on the fully charged protein system. When the energy due to the solvent reaction field in the uncharged system is accounted for, correlation between per-atom electrostatic energies for the explicit solvent model and a simple implicit (Poisson) calculation is 0.97, and correlation between per-atom energies for the explicit solvent model and a previously published, optimized Poisson model is 0.99.