3D RISM theory with fast reciprocal-space electrostatics.

The calculation of electrostatic solute-solvent interactions in 3D RISM ("three-dimensional reference interaction site model") integral equation theory is recast in a form that allows for a computational treatment analogous to the "particle-mesh Ewald" formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.

[1]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[2]  M. Fedorov,et al.  3DRISM Multigrid Algorithm for Fast Solvation Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[3]  Carlos Simmerling,et al.  Three-dimensional molecular theory of solvation coupled with molecular dynamics in Amber. , 2010, Journal of chemical theory and computation.

[4]  Andriy Kovalenko,et al.  Efficient treatment of solvation shells in 3D molecular theory of solvation , 2012, J. Comput. Chem..

[5]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[6]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[7]  Fumio Hirata,et al.  Ab initio study of water. II. Liquid structure, electronic and thermodynamic properties over a wide range of temperature and density , 1999 .

[8]  Yasuomi Kiyota,et al.  The electronic-structure theory of a large-molecular system in solution: Application to the intercalation of proflavine with solvated DNA , 2011 .

[9]  Tohru Morita,et al.  A New Approach to the Theory of Classical Fluids. I , 1960 .

[10]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[11]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[12]  B. Montgomery Pettitt,et al.  A site-site theory for finite concentration saline solutions , 1992 .

[13]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[14]  Andriy Kovalenko,et al.  An MM/3D-RISM approach for ligand binding affinities. , 2010, The journal of physical chemistry. B.

[15]  S. Ten-no,et al.  Ab initio study of water: self-consistent determination of electronic structure and liquid state properties , 1997 .

[16]  Robert C. Harris,et al.  Influence of Grid Spacing in Poisson-Boltzmann Equation Binding Energy Estimation. , 2013, Journal of chemical theory and computation.

[17]  Irwin Oppenheim,et al.  Statistical Mechanical Theory of Transport Processes. VII. The Coefficient of Thermal Conductivity of Monatomic Liquids , 1954 .

[18]  H. Schwalbe,et al.  Model development for the viral Kcv potassium channel. , 2009, Biophysical journal.

[19]  R. Gebauer,et al.  Density functional theory of solvation in a polar solvent: extracting the functional from homogeneous solvent simulations. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  Christian Holm,et al.  How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines , 1998 .

[21]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[22]  H. G. Petersen.  Accuracy and efficiency of the particle mesh Ewald method , 1995 .

[23]  P. Pulay.  Convergence acceleration of iterative sequences. the case of scf iteration , 1980 .

[24]  S. Kast,et al.  Molecular dynamics simulation of the cytosolic mouth in Kcv-type potassium channels. , 2007, Biochemistry.

[25]  M. Levesque,et al.  Molecular Density Functional Theory of Water. , 2013, The journal of physical chemistry letters.

[26]  Stefan M. Kast,et al.  A minimalist model for ion partitioning and competition in a K+ channel selectivity filter , 2011, The Journal of general physiology.

[27]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[28]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[29]  Daniel Potts,et al.  Parallel Three-Dimensional Nonequispaced Fast Fourier Transforms and Their Application to Particle Simulation , 2013, SIAM J. Sci. Comput..

[30]  Communication: An exact bound on the bridge function in integral equation theories. , 2012, The Journal of chemical physics.

[31]  D. Beglov,et al.  Finite representation of an infinite bulk system: Solvent boundary potential for computer simulations , 1994 .

[32]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[33]  Hartmut Löwen,et al.  Highly asymmetric electrolytes in the primitive model: Hypernetted chain solution in arbitrary spatial dimensions , 2013, J. Comput. Chem..

[34]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[35]  Gillian C. Lynch,et al.  Protein solvation from theory and simulation: Exact treatment of Coulomb interactions in three-dimensional theories. , 2010, The Journal of chemical physics.

[36]  Maxim V Fedorov,et al.  Combination of RISM and Cheminformatics for Efficient Predictions of Hydration Free Energy of Polyfragment Molecules: Application to a Set of Organic Pollutants. , 2011, Journal of chemical theory and computation.

[37]  Fumio Hirata,et al.  Solution of three‐dimensional reference interaction site model and hypernetted chain equations for simple point charge water by modified method of direct inversion in iterative subspace , 1999 .

[38]  Stefan M. Kast,et al.  Prediction of tautomer ratios by embedded-cluster integral equation theory , 2010, J. Comput. Aided Mol. Des..

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

[40]  Andriy Kovalenko,et al.  Modeling solvatochromic shifts using the orbital-free embedding potential at statistically mechanically averaged solvent density. , 2010, The journal of physical chemistry. A.

[41]  Stefan M Kast,et al.  Treatment of charged solutes in three-dimensional integral equation theory. , 2008, Journal of Chemical Physics.

[42]  B. Montgomery Pettitt,et al.  A dielectrically consistent interaction site theory for solvent—electrolyte mixtures , 1992 .

[43]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[44]  R. Kjellander,et al.  The chemical potential of simple fluids in a common class of integral equation closures , 1989 .

[45]  J. D. Talman,et al.  Numerical Fourier and Bessel transforms in logarithmic variables , 1978 .

[46]  Fumio Hirata,et al.  Modified Anderson Method for Accelerating 3D-RISM Calculations Using Graphics Processing Unit. , 2012, Journal of chemical theory and computation.

[47]  S. Kast.  Free energies from integral equation theories: enforcing path independence. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[48]  D. Case,et al.  Simple electrolyte solutions: comparison of DRISM and molecular dynamics results for alkali halide solutions. , 2013, The Journal of chemical physics.

[49]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[50]  V. Lebedev,et al.  A QUADRATURE FORMULA FOR THE SPHERE OF THE 131ST ALGEBRAIC ORDER OF ACCURACY , 1999 .

[51]  B. Pettitt,et al.  An Integral Equation Study of the Hydrophobic Interaction between Graphene Plates. , 2008, Journal of chemical theory and computation.

[52]  Jochen Heil,et al.  Quantum chemistry in solution by combining 3D integral equation theory with a cluster embedding approach. , 2008, The journal of physical chemistry. B.

[53]  S. Kast,et al.  Closed-form expressions of the chemical potential for integral equation closures with certain bridge functions. , 2008, The Journal of chemical physics.

[54]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[55]  C Holm,et al.  How to Convert SPME to P3M: Influence Functions and Error Estimates. , 2012, Journal of chemical theory and computation.

[56]  F. Hirata,et al.  Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach , 1998 .

[57]  Fumio Hirata,et al.  Placevent: An algorithm for prediction of explicit solvent atom distribution—Application to HIV‐1 protease and F‐ATP synthase , 2012, J. Comput. Chem..

[58]  David Chandler,et al.  Free energy functions in the extended RISM approximation , 1985 .

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