Quantum mechanical geometry optimization in solution using a finite element continuum electrostatics method

We present a new algorithm for performing ab initio solution phase geometry optimizations. The procedure is based on the self consistent‐reaction‐field method developed in our laboratory which combines electronic structure calculations with a finite element formulation of the continuum electrostatics problem. A gradient for the total solution phase free energy is obtained by combining different contributions from the gradient of the classical polarization free energy and the derivatives of the quantum mechanical energy. The method used in obtaining the classical gradient is based on exact linear algebra relations and a Green function formalism due to Handy and Schaefer. Both the classical and quantum mechanical gradients are validated by comparison with energy finite differences. The result of applications to a number of small organic compounds are discussed. Comparisons between the predicted location and depth of the various solution phase minima of the Ramachandran map for the alanine dipeptide and thos...

[1]  T. J. Baker,et al.  Developments and trends in three-dimensional mesh generation , 1989 .

[2]  Ian R. Gould,et al.  A quantum Mechanical Investigation of the Conformational Energetics of the Alanine and Glycine Dipeptides in the Gas Phase and in Aqueous Solution , 1994 .

[3]  B. Honig,et al.  Accurate First Principles Calculation of Molecular Charge Distributions and Solvation Energies from Ab Initio Quantum Mechanics and Continuum Dielectric Theory , 1994 .

[4]  Arieh Ben-Naim,et al.  Solvation thermodynamics of nonionic solutes , 1984 .

[5]  N. L. Allinger,et al.  Molecular Mechanics (MM3). Calculations of Furan, Vinyl Ethers, and Related Compounds , 1993 .

[6]  Charles L. Brooks,et al.  Theoretical study of blocked glycine and alanine peptide analogs , 1991 .

[7]  D. Rinaldi,et al.  Fast geometry optimizationin self‐cosistent reaction field computations on solvated molecules , 1992 .

[8]  Barry Honig,et al.  Focusing of electric fields in the active site of Cu‐Zn superoxide dismutase: Effects of ionic strength and amino‐acid modification , 1986, Proteins.

[9]  J Kraut,et al.  Theoretical studies on the dihydrofolate reductase mechanism: electronic polarization of bound substrates. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[10]  H. Schaefer,et al.  Generalization of analytic energy third derivatives for the RHF closed‐shell wave function: Derivative energy and integral formalisms and the prediction of vibration–rotation interaction constants , 1986 .

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

[12]  Daniel Rinaldi,et al.  Lamé's functions and ellipsoidal harmonics for use in chemical physics , 1982, Comput. Chem..

[13]  Jacopo Tomasi,et al.  Analytical derivatives for molecular solutes. II. Hartree–Fock energy first and second derivatives with respect to nuclear coordinates , 1994 .

[14]  Z. Maksić,et al.  Modelling of structure and properties of molecules , 1987 .

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

[16]  H A Scheraga,et al.  Ab initio self‐consistent field and potential‐dependent partial equalization of orbital electronegativity calculations of hydration properties of N‐acetyl‐N′‐methyl‐alanineamide , 1990, Biopolymers.

[17]  B. Honig,et al.  New Model for Calculation of Solvation Free Energies: Correction of Self-Consistent Reaction Field Continuum Dielectric Theory for Short-Range Hydrogen-Bonding Effects , 1996 .

[18]  Jacopo Tomasi,et al.  Analytical derivatives for molecular solutes. I. Hartree–Fock energy first derivatives with respect to external parameters in the polarizable continuum model , 1994 .

[19]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[20]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[21]  R. Friesner,et al.  Pseudospectral Hartree-Fock gradient calculations , 1991 .

[22]  K. Kopple,et al.  Solvent-dependent conformational distributions of some dipeptides , 1980 .

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