A General Treatment of Solvent Effects Based on Screened Coulomb Potentials

A new implicit solvent model for Monte Carlo simulations of peptides and proteins is presented. The important features of the model are:  electrostatics of polar solvation described by nonlinear, distance-dependent screened Coulomb potentials that are formally connected to classical dielectric theory of polar solvation; the self-energy is derived from the integral form of the Born equation; and the introduction of the Born radius of an atom embedded in a protein modeled as a dielectric continuum. The model was parametrized in the context of the CHARMM PAR22 force field by fitting to experimental solvation free energies of amino acid side chain analogues. Introduction of an intrinsic vacuum screening constant for each amino acid was required to obtain the correct screening behavior in solution. Solvation energies of several peptides in different conformations were calculated and compared with the corresponding quantities calculated with the Poisson−Boltzmann equation. The comparative results were highly co...