An Integrated Effective Fragment—Polarizable Continuum Approach to Solvation: Theory and Application to Glycine

A new discrete/continuum solvation model has been developed by combining the effective fragment potential (EFP) for the discrete part and the polarizable continuum model (PCM) for the continuum part. The usefulness of this model is demonstrated by applying it to the calculation of the relative energies of the neutral and zwitterionic forms of glycine. These calculations were performed by treating glycine with ab initio wave functions. Water clusters were treated with both ab initio and EFP methods for comparison purposes, and the effect of the continuum was accounted for by the PCM model. The energy barrier connecting the zwitterionic and neutral three-water clusters was also examined. The computationally efficient EFP/PCM model gives results that are in close agreement with the much more expensive full ab initio/PCM calculation. The use of methods that account for electron correlation is necessary to obtain accurate relative energies for the isomers of glycine.

[1]  J. Langlet,et al.  Theoretical study of solvent effect on intramolecular proton transfer of glycine , 2000 .

[2]  Mark S. Gordon,et al.  A combined discrete/continuum solvation model: Application to glycine , 2000 .

[3]  I. Tuñón,et al.  On the tautomerization process of glycine in aqueous solution , 2000 .

[4]  Jacopo Tomasi,et al.  Glycine and alanine: a theoretical study of solvent effects upon energetics and molecular response properties , 2000 .

[5]  Jacopo Tomasi,et al.  Second-Order Møller-Plesset Analytical Derivatives for the Polarizable Continuum Model Using the Relaxed Density Approach , 1999 .

[6]  Seiichiro Ten-no,et al.  Ab initio Monte Carlo simulation using multicanonical algorithm: temperature dependence of the average structure of water dimer , 1999 .

[7]  I. Tuñón,et al.  Aminoacid zwitterions in solution: Geometric, energetic, and vibrational analysis using density functional theory-continuum model calculations , 1998 .

[8]  Benedetta Mennucci,et al.  New applications of integral equations methods for solvation continuum models: ionic solutions and liquid crystals , 1998 .

[9]  Jacopo Tomasi,et al.  Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications , 1997 .

[10]  Jacopo Tomasi,et al.  A new integral equation formalism for the polarizable continuum model: Theoretical background and applications to isotropic and anisotropic dielectrics , 1997 .

[11]  V. Barone,et al.  SOLVENT EFFECTS ON THE CONFORMATIONAL BEHAVIOR OF MODEL PEPTIDES. A COMPARISON BETWEEN DIFFERENT CONTINUUM MODELS , 1996 .

[12]  Iñaki Tuñón,et al.  Why is glycine a zwitterion in aqueous solution? A theoretical study of solvent stabilising factors , 1996 .

[13]  Andreas Klamt,et al.  Incorporation of solvent effects into density functional calculations of molecular energies and geometries , 1995 .

[14]  Jacopo Tomasi,et al.  Remarks on the use of the apparent surface charges (ASC) methods in solvation problems: Iterative versus matrix‐inversion procedures and the renormalization of the apparent charges , 1995, J. Comput. Chem..

[15]  T. Truong,et al.  Analytical first and second energy derivatives of the generalized conductorlike screening model for free energy of solvation , 1995 .

[16]  M. Gordon,et al.  On the Number of Water Molecules Necessary To Stabilize the Glycine Zwitterion , 1995 .

[17]  A. H. de Vries,et al.  Implementation of reaction field methods in quantum chemistry computer codes , 1995, J. Comput. Chem..

[18]  Jacopo Tomasi,et al.  Molecular Interactions in Solution: An Overview of Methods Based on Continuous Distributions of the Solvent , 1994 .

[19]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[20]  R. Levy,et al.  Molecular dynamics simulation of time-resolved fluorescence and nonequilibrium solvation of formaldehyde in water , 1990 .

[21]  W. J. Stevens,et al.  Transferability of molecular distributed polarizabilities from a simple localized orbital based method , 1989 .

[22]  H. Scheraga,et al.  Monte Carlo-minimization approach to the multiple-minima problem in protein folding. , 1987, Proceedings of the National Academy of Sciences of the United States of America.

[23]  A. Warshel,et al.  Calculations of electrostatic energies in proteins. The energetics of ionized groups in bovine pancreatic trypsin inhibitor. , 1985, Journal of molecular biology.

[24]  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 .

[25]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[26]  M. Okina,et al.  On the Ratio of Zwitterion Form to Uncharged Form of Glycine at Equilibrium in Various Aqueous Media , 1982 .

[27]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[28]  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 .

[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]  D. D. Yue,et al.  Theory of Electric Polarization , 1974 .

[31]  M. Plesset,et al.  Note on an Approximation Treatment for Many-Electron Systems , 1934 .

[32]  Mark S. Gordon,et al.  The Effective Fragment Potential Method: A QM-Based MM Approach to Modeling Environmental Effects in Chemistry , 2001 .

[33]  H. Rzepa,et al.  An AM1 and PM3 molecular orbital and self-consistent reaction-field study of the aqueous solvation of glycine, alanine and proline in their neutral and zwitterionic forms , 1991 .

[34]  N. P. Buu‐Hoï,et al.  Carcinogenic nitrogen compounds. Part LXXVII. A novel synthesis of β-carbolines , 1972 .