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‐difluoreothane (II), fluoroacetic acid (III), 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–III), 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 shown for VI in water. The MST‐PC model which mimics the solute–solvent electrostatic interaction only fails to reproduce conformationally dependent solvation energies in acetonitrile, in which the electrostatic contribution is relatively small compared to van der Waals interactions. Exploratory MD simulations within the continuum model in water are reported on urea and 18‐crown‐6. © 1995 by John Wiley & Sons, Inc.

[1]  C. Ghio,et al.  Towards a unified view of the description of internal and external fields acting on chemical functional groups , 1988 .

[2]  L. Troxler,et al.  Conformation and Dynamics of 18-Crown-6, Cryptand 222, and Their Cation Complexes in Acetonitrile Studied by Molecular Dynamics Simulations , 1994 .

[3]  J. Tomasi.  Description and interpretation of molecular phenomena in solution, using effective Hamiltonian operators related to continuous solvent distributions , 1991 .

[4]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[5]  Thomas Simonson,et al.  Solvation Free Energies Estimated from Macroscopic Continuum Theory: An Accuracy Assessment , 1994 .

[6]  Frank A. Momany,et al.  Determination of partial atomic charges from ab initio molecular electrostatic potentials. Application to formamide, methanol, and formic acid , 1978 .

[7]  F. Javier Luque,et al.  Comparison of 6‐31G*‐based MST/SCRF and FEP evaluations of the free energies of hydration for small neutral molecules , 1993, J. Comput. Chem..

[8]  R. Cimiraglia,et al.  On the ab initio evaluation of the solvent shift of electronic absorption spectra , 1981 .

[9]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[10]  R. Zauhar,et al.  The rigorous computation of the molecular electric potential , 1988 .

[11]  P. Claverie,et al.  Theoretical evaluations of the intermolecular interaction energy of a crystal: application to the analysis of crystal geometry , 1975 .

[12]  Jacopo Tomasi,et al.  Approximate evaluations of the electrostatic free energy and internal energy changes in solution processes , 1982 .

[13]  J. Tomasi,et al.  Additivity and non-additivity effects in the conformational changes of 1,2-ethanediol diformatein vacuo and in solution.Ab-initio and molecular mechanics descriptions , 1992 .

[14]  C. Ghio,et al.  Theoretical calculations on 1,2-ethanediol. 2. Equilibrium of the gauche conformers with and without an intramolecular hydrogen bond in aqueous solution , 1992 .

[15]  E. Westhof,et al.  Molecular dynamics simulations of poly(dA)·poly(dT): Comparisons between implicit and explicit solvent representations , 1993, Biopolymers.

[16]  N. L. Allinger,et al.  Conformational analysis—132: The effects of electrostatic interactions and solvation energies on conformational equilibria in dihalides and haloketonesTetrahedron 34, 627, 1978. , 1978 .

[17]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[18]  G. Wipff,et al.  Dramatic solvent effect on the ligand wrapping around a complexed cation: a MD study of p-tert-butylcalix[4]arene tetramide and its complexes with alkali cations and europium(3+) in acetonitrile , 1993 .

[19]  R. Bonaccorsi,et al.  Conformational energy of glycine in aqueous solutions and relative stability of the zwitterionic and neutral forms. An ab initio study , 1984 .

[20]  Alexander A. Rashin,et al.  A simple method for the calculation of hydration enthalpies of polar molecules with arbitrary shapes , 1987 .

[21]  W. C. Still,et al.  A rapid approximation to the solvent accessible surface areas of atoms , 1988 .

[22]  Norman L. Allinger,et al.  Conformational analysis. 130. MM2. A hydrocarbon force field utilizing V1 and V2 torsional terms , 1977 .

[23]  Irwin D. Kuntz,et al.  A new technique to calculate low-energy conformations of cyclic molecules utilizing the ellipsoid algorithm and molecular dynamics: application to 18-crown-6 , 1988 .

[24]  D. Pérahia,et al.  Hydration scheme of the purine and pyrimidine bases and base-pairs of the nucleic acids , 1979 .

[25]  J. Lehn,et al.  Monte Carlo study of the conformation-dependent hydration of the 18-crown-6 macrocycle , 1985 .

[26]  Barry Robson,et al.  Intermolecular Interactions: from diatomics to biopolymers , 1978 .

[27]  F. J. Luque,et al.  An AM1-SCRF approach to the study of changes in molecular properties induced by solvent , 1993 .

[28]  P. Claverie,et al.  Improvements of the continuum model. 1. Application to the calculation of the vaporization thermodynamic quantities of nonassociated liquids , 1988 .

[29]  P. Claverie,et al.  Studies of solvent effects. 1. Discrete, continuum, and discrete-continuum models and their comparison for some simple cases: ammonium(1+) ion, methanol, and substituted ammonium(1+) ion , 1978 .

[30]  Norman L. Allinger,et al.  Treatment of electrostatic effects within the molecular-mechanics method. 1 , 1983 .

[31]  G. Chang,et al.  Macromodel—an integrated software system for modeling organic and bioorganic molecules using molecular mechanics , 1990 .

[32]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[33]  R. Zauhar,et al.  A new method for computing the macromolecular electric potential. , 1985, Journal of molecular biology.

[34]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[35]  Harold A. Scheraga,et al.  A combined iterative and boundary-element approach for solution of the nonlinear Poisson-Boltzmann equation , 1992 .

[36]  G. W. Schnuelle,et al.  A statistical thermodynamic supermolecule-continuum study of ion hydration: Cell and shell methods , 1978 .

[37]  George P. Ford,et al.  Incorporation of hydration effects within the semiempirical molecular orbital framework. AM1 and MNDO results for neutral molecules, cations, anions, and reacting systems , 1992 .

[38]  M. Huron,et al.  Calculation of the interaction energy of one molecule with its whole surrounding. II. Method of calculating electrostatic energy , 1974 .

[39]  J. A. McCammon,et al.  Calculating electrostatic forces from grid‐calculated potentials , 1990 .

[40]  A. Imberty,et al.  Solvent effect on the stability of isomaltose conformers , 1990 .

[41]  S. Harvey.  Treatment of electrostatic effects in macromolecular modeling , 1989, Proteins.

[42]  W F van Gunsteren,et al.  On the dependence of molecular conformation on the type of solvent environment: A molecular dynamics study of cyclosporin A , 1990, Biopolymers.

[43]  Randy J. Zauhar,et al.  Computing the electric potential of biomolecules: Application of a new method of molecular surface triangulation , 1990 .

[44]  L. Troxler,et al.  MD Simulations on Synthetic Ionophores and Their Cation Complexes: Comparison of Aqueous/Non-Aqueous Solvents , 1994 .

[45]  T. P. Straatsma,et al.  Treatment of rotational isomers in free energy calculations. II. Molecular dynamics simulation study of 18‐crown‐6 in aqueous solution as an example of systems with large numbers of rotational isomeric states , 1989 .

[46]  D. Pérahia,et al.  Hydration scheme of uracil and cytosine , 1978 .

[47]  J. Tomasi,et al.  Dispersion and repulsion contributions to the solvation energy: Refinements to a simple computational model in the continuum approximation , 1991 .

[48]  Donald G. Truhlar,et al.  Correlation and solvation effects on heterocyclic equilibria in aqueous solution , 1993 .

[49]  R. Lavery,et al.  Theoretical Biochemistry & Molecular Biophysics. , 1990, Journal of Biomolecular Structure and Dynamics.

[50]  A. Gamba,et al.  Molecular modelling of organophosphorus podands and their complexes with alkali metal cations , 1992 .

[51]  D. Levitt,et al.  Electrostatic calculations for an ion channel. I. Energy and potential profiles and interactions between ions. , 1978, Biophysical journal.

[52]  Alexander A. Rashin,et al.  Hydration phenomena, classical electrostatics, and the boundary element method , 1990 .

[53]  J. Tomasi,et al.  Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects , 1981 .

[54]  P. Claverie,et al.  Calculation of the interaction energy of one molecule with its whole surrounding. I. Method and application to pure nonpolar compounds , 1972 .

[55]  D. Feil,et al.  Molecular dynamics of 18‐crown‐6 complexes with alkali–metal cations and urea: Prediction of their conformations and comparison with data from the cambridge structural database , 1993, J. Comput. Chem..

[56]  B Pullman,et al.  New Paths in the Molecular Orbital Approach to Solvation of Biological Molecules , 1974, Quarterly Reviews of Biophysics.

[57]  A. Varnek,et al.  Complexation of lithium and sodium cations with β-phosphorylate ethers, modelling terminal groups of organophosphorus podands. An experimental and theoretical study , 1992 .

[58]  R. Bonaccorsi,et al.  On the ab initio geometry optimization of molecular solutes , 1991 .

[59]  William J. Dunn,et al.  Theoretical calculations on 1,2-ethanediol. Gauche-trans equilibrium in gas-phase and aqueous solution , 1991 .

[60]  G. W. Schnuelle,et al.  Free energy of a charge distribution in concentric dielectric continua , 1975 .

[61]  E. Westhof,et al.  Three-center hydrogen bonds in DNA: molecular dynamics of poly(dA).cntdot.poly(dT) , 1991 .

[62]  J. Andrew McCammon,et al.  Dielectric boundary smoothing in finite difference solutions of the poisson equation: An approach to improve accuracy and convergence , 1991 .

[63]  M. L. Connolly.  Analytical molecular surface calculation , 1983 .

[64]  Kim A. Sharp,et al.  Incorporating solvent and ion screening into molecular dynamics using the finite‐difference Poisson–Boltzmann method , 1991 .

[65]  A. Rashin.  Continuum electrostatics and hydration phenomena , 1988 .

[66]  Giuliano Alagona,et al.  The role of electrostatics in solute-solvent interactions with the continuum , 1992 .

[67]  A. Varnek,et al.  Molecular dynamics study of p-tert-butylcalix[4]arenetetraamide and its complexes with neutral and cationic guests. Influence of solvation on structures and stabilities , 1993 .

[68]  R. Lavery,et al.  The flexibility of the nucleic acids: (II). The calculation of internal energy and applications to mononucleotide repeat DNA. , 1986, Journal of biomolecular structure & dynamics.

[69]  Malcolm E. Davis,et al.  Electrostatics in biomolecular structure and dynamics , 1990 .

[70]  J. Andrew McCammon,et al.  Computation of electrostatic forces on solvated molecules using the Poisson-Boltzmann equation , 1993 .

[71]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[72]  H. Berendsen,et al.  A consistent empirical potential for water–protein interactions , 1984 .

[73]  J. Dunitz,et al.  Potassium thiocyanate complex of 1,4,7,10,13,16-hexaoxacyclooctadecane , 1974 .

[74]  A. Pullman.  The Supermolecule Approach to the Solvation Problem , 1980 .

[75]  R. Fèvre,et al.  Molecular Refractivity and Polarizability , 1965 .

[76]  P. Auffinger,et al.  Hydration of the 222 cryptand and 222 cryptates studied by molecular dynamics simulations , 1991 .

[77]  J. C. Contador,et al.  Applications of the basic polarizable continuum model. A study of the vibrational properties of diatomic solutes , 1990 .

[78]  Jacopo Tomasi,et al.  Electrostatic interaction of a solute with a continuum. Improved description of the cavity and of the surface cavity bound charge distribution. , 1987 .

[79]  J. Dunitz,et al.  Hydrated sodium thiocyanate complex of 1,4,7,10,13,16‐hexaoxacyclooctadecane , 1974 .

[80]  Arieh Warshel,et al.  Microscopic and semimicroscopic calculations of electrostatic energies in proteins by the POLARIS and ENZYMIX programs , 1993, J. Comput. Chem..

[81]  R. Zauhar,et al.  The incorporation of hydration forces determined by continuum electrostatics into molecular mechanics simulations , 1991 .

[82]  M. Huron,et al.  Calculation of the interaction energy of one molecule with its whole surrounding. III. Application to pure polar compounds , 1974 .

[83]  George P. Ford,et al.  The optimized ellisoidal cavity and its application to the self‐consistent reaction field calculation of hydration energies of cations and neutral molecules , 1992 .

[84]  J. Dunitz,et al.  1,4,7,10,13,16‐Hexaoxacyclooctadecane: crystal structure at 100 K , 1980 .

[85]  J. Tomasi,et al.  Decomposition of the interaction energy with counterpoise corrections to the basis set superposition error for dimers in solution. Method and application to the hydrogen fluoride dimer , 1988 .