An Application of the Miertus-Scrocco-Tomasi Solvation Model in Molecular Mechanics and Dynamics Simulations

The point-chart approximation of the Miertus-Scrocco-Tomasi solvation model (MST-PC) based on a continuum representation of the solvent has been incorporated in force field calculations. Application in molecular mechanics (MM) involves conformational equilibria in solution: rotational isomers of ethylene glycol (I), 1,2-difluoroethane (10, fluoroacetic acid (110, and representative conformers of macrocyclic receptors such as 18-crown-6 (IV), cryptand 2.2.2 (V), and t-butyl-calix[4]arenetetraamide (VI). Assessment of the MST-PC results is based on the comparison with ab initio reactive field calculations (for I-111), with the continuum model of Still (W. C. Still et al., J. Am. Chem. SOC., 1990, 112, 6127) (for I-VI), and with average solute-solvent interaction energies obtained from molecular dynamics (MD) simulations with explicit solvent in water (I-VI) and in acetonitrile (IV-VI). It is demonstrated that the continuum solvent model qualitatively reproduces the trends in solvation energies in water. The few exceptions may be related to particular topological features of the solute. An improved discrete/continuum approach in which some first-shell solvent molecules are considered as a part of the solute embedded in the dielectric continuum provides more realistic results, as is *Author to wh,om all correspondence should be addressed