Efficient quantum mechanical calculation of solvation free energies based on density functional theory, numerical atomic orbitals and Poisson–Boltzmann equation

Abstract We have successfully coupled the Kohn–Sham with Poisson–Boltzmann equations to predict the solvation free energy, where the Kohn–Sham equations were solved by implementing the flexible pseudo atomic orbitals as in S iesta package. It was found that the calculated solvation free energy is in good agreement with experimental results for small neutral molecules, and its standard error is 1.33 kcal/mol, the correlation coefficient is 0.97. Due to its high efficiency and accuracy, the proposed model can be a promising tool for computing solvation free energies in computer aided drug design in future.

[1]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[2]  Darrin M. York,et al.  Parameterization and efficient implementation of a solvent model for linear-scaling semiempirical quantum mechanical calculations of biological macromolecules , 1996 .

[3]  C. Cramer,et al.  General parameterized SCF model for free energies of solvation in aqueous solution , 1991 .

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

[5]  D. Sánchez-Portal,et al.  The SIESTA method for ab initio order-N materials simulation , 2001, cond-mat/0111138.

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

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

[8]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[9]  Richard A. Friesner,et al.  Quantum mechanical geometry optimization in solution using a finite element continuum electrostatics method , 1996 .

[10]  A. Bondi.  van der Waals Volumes and Radii , 1964 .

[11]  Kenneth M. Merz,et al.  Fully Quantum Mechanical Description of Proteins in Solution. Combining Linear Scaling Quantum Mechanical Methodologies with the Poisson−Boltzmann Equation , 1999 .

[12]  D. Case,et al.  Generalized born models of macromolecular solvation effects. , 2000, Annual review of physical chemistry.

[13]  Chung F Wong,et al.  Calculation of solvation free energy from quantum mechanical charge density and continuum dielectric theory. , 2006, The journal of physical chemistry. A.

[14]  A. Klamt,et al.  COSMO : a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient , 1993 .

[15]  Haiyan Liu,et al.  The treatment of solvation by a generalized Born model and a self‐consistent charge‐density functional theory‐based tight‐binding method , 2002, J. Comput. Chem..

[16]  Incorporating Protein Environments in Density Functional Theory: A Self-Consistent Reaction Field Calculation of Redox Potentials of [2Fe2S] Clusters in Ferredoxin and Phthalate Dioxygenase Reductase , 1998 .

[17]  F. J. Luque,et al.  Theoretical Methods for the Description of the Solvent Effect in Biomolecular Systems. , 2000, Chemical reviews.

[18]  J. A. McCammon,et al.  Solving the finite difference linearized Poisson‐Boltzmann equation: A comparison of relaxation and conjugate gradient methods , 1989 .

[19]  Nathan A. Baker,et al.  Electrostatics of nanosystems: Application to microtubules and the ribosome , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[21]  C. Cramer,et al.  Implicit Solvation Models: Equilibria, Structure, Spectra, and Dynamics. , 1999, Chemical reviews.