Generalized image charge solvation model for electrostatic interactions in molecular dynamics simulations of aqueous solutions

This paper extends the image charge solvation model (ICSM) [J. Chem. Phys. 131, 154103 (2009)], a hybrid explicit/implicit method to treat electrostatic interactions in computer simulations of biomolecules formulated for spherical cavities, to prolate spheroidal and triaxial ellipsoidal cavities, designed to better accommodate non-spherical solutes in molecular dynamics (MD) simulations. In addition to the utilization of a general truncated octahedron as the MD simulation box, central to the proposed extension is an image approximation method to compute the reaction field for a point charge placed inside such a non-spherical cavity by using a single image charge located outside the cavity. The resulting generalized image charge solvation model (GICSM) is tested in simulations of liquid water, and the results are analyzed in comparison with those obtained from the ICSM simulations as a reference. We find that, for improved computational efficiency due to smaller simulation cells and consequently a less number of explicit solvent molecules, the generalized model can still faithfully reproduce known static and dynamic properties of liquid water at least for systems considered in the present paper, indicating its great potential to become an accurate but more efficient alternative to the ICSM when bio-macromolecules of irregular shapes are to be simulated.

[1]  C. Brooks,et al.  Recent advances in the development and application of implicit solvent models in biomolecule simulations. , 2004, Current opinion in structural biology.

[2]  George Dassios,et al.  Magnetoencephalography in ellipsoidal geometry , 2003 .

[3]  Shaozhong Deng,et al.  Electrostatic potential of point charges inside dielectric prolate spheroids. , 2008, Journal of electrostatics.

[4]  R. Abagyan,et al.  Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. , 1994, Journal of molecular biology.

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

[6]  Johan C.-E. Sten.  Ellipsoidal harmonics and their application in electrostatics , 2006 .

[7]  Ronald M. Levy,et al.  Dielectric and thermodynamic response of a generalized reaction field model for liquid state simulations , 1993 .

[8]  Wei Cai,et al.  A comparable study of image approximations to the reaction field , 2007, Comput. Phys. Commun..

[9]  D. A. Dunnett.  Classical Electrodynamics , 2020, Nature.

[10]  J. W. Perram,et al.  On the application of ellipsoidal harmonics to potential problems in molecular electrostatics and magnetostatics , 1976, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[11]  Shaozhong Deng.  Electrostatic potential of point charges inside dielectric oblate spheroids. , 2008, Journal of electrostatics.

[12]  P. Harbury,et al.  Tanford-Kirkwood electrostatics for protein modeling. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[13]  A. Wallqvist.  On the Implementation of Friedman Boundary Conditions in Liquid Water Simulations , 1993 .

[14]  D. Redžić.  An electrostatic problem: A point charge outside a prolate dielectric spheroid , 1994 .

[15]  P. Koehl.  Electrostatics calculations: latest methodological advances. , 2006, Current opinion in structural biology.

[16]  Bo Zhang,et al.  FMM-Yukawa: An adaptive fast multipole method for screened Coulomb interactions , 2009, Comput. Phys. Commun..

[17]  Mark A Olson,et al.  Evaluation of Poisson solvation models using a hybrid explicit/implicit solvent method. , 2005, The journal of physical chemistry. B.

[18]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

[19]  T. Darden,et al.  Molecular dynamics simulations of biomolecules: long-range electrostatic effects. , 1999, Annual review of biophysics and biomolecular structure.

[20]  W. R. Smythe.  Static and Dynamic Electricity , 1989 .

[21]  D. Redz̆ić,et al.  Image charge inclusions in the prolate dielectric spheroid , 2005 .

[22]  Gregory A. Voth,et al.  Simple reversible molecular dynamics algorithms for Nosé-Hoover chain dynamics , 1997 .

[23]  G. Romain,et al.  Ellipsoidal Harmonic expansions of the gravitational potential: Theory and application , 2001 .

[24]  Nathan A. Baker,et al.  Improving implicit solvent simulations: a Poisson-centric view. , 2005, Current opinion in structural biology.

[25]  Changfeng Xue,et al.  Three-layer dielectric models for generalized Coulomb potential calculation in ellipsoidal geometry. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[27]  L. Greengard,et al.  A new version of the fast multipole method for screened Coulomb interactions in three dimensions , 2002 .

[28]  J. Kirkwood,et al.  Statistical Mechanics of Liquid Solutions. , 1936 .

[29]  P. Morse,et al.  Methods of theoretical physics , 1955 .

[30]  J. A. Barker,et al.  Monte Carlo studies of the dielectric properties of water-like models , 1973 .

[31]  D. V. Redzć.  COMMENT ON "ELECTROSTATIC IMAGE THEORY FOR THE DIELECTRIC PROLATE SPHEROID" BY I. V. LINDELL AND K. I. NIKOSKINEN , 2003 .

[32]  J. Kirkwood,et al.  Theory of Solutions of Molecules Containing Widely Separated Charges with Special Application to Zwitterions , 1934 .

[33]  Keijo Nikoskinen,et al.  Electrostatic image theory for the conducting prolate spheroid , 2000 .

[34]  W. T. Norris.  Charge images in a dielectric sphere , 1995 .

[35]  The Ultimate Image Singularities for External Ellipsoidal Harmonics , 1974 .

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

[37]  M. Parrinello,et al.  Nonperiodic boundary conditions for solvated systems. , 2005, The Journal of chemical physics.

[38]  L. Greengard.  The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[39]  J. Hermans,et al.  REACTION FIELD MOLECULAR DYNAMICS SIMULATION WITH FRIEDMAN'S IMAGE CHARGE METHOD , 1995 .

[40]  Bo Zhang,et al.  FMM-Yukawa: An adaptive fast multipole method for screened Coulomb interactions , 2009, Comput. Phys. Commun..

[41]  Mark A Olson,et al.  An efficient hybrid explicit/implicit solvent method for biomolecular simulations , 2004, J. Comput. Chem..

[42]  P. Hünenberger,et al.  Alternative schemes for the inclusion of a reaction-field correction into molecular dynamics simulations: Influence on the simulated energetic, structural, and dielectric properties of liquid water , 1998 .

[43]  Keijo Nikoskinen,et al.  Electrostatic Image Theory for the Dielectric Prolate Spheroid , 2001 .

[44]  M. Willatzen,et al.  On triaxial ellipsoidal quantum dots , 2004 .

[45]  Benzhuo Lu,et al.  Order N algorithm for computation of electrostatic interactions in biomolecular systems , 2006, Proceedings of the National Academy of Sciences.

[46]  H.-J. Dobner,et al.  Verified computation of Lamé functions with high accuracy , 2007, Computing.

[47]  Wei Cai,et al.  Extending the fast multipole method to charges inside or outside a dielectric sphere , 2007, J. Comput. Phys..

[48]  J. Rullmann,et al.  Analysis of discrete and continuum dielectric models : Application to the calculation of protonation energies in solution , 1987 .

[49]  E. Hobson.  The Theory of Spherical and Ellipsoidal Harmonics , 1955 .

[50]  Ismo V. Lindell,et al.  Electrostatic image solution for the conducting prolate spheroid , 1993 .

[51]  D. Jacobs,et al.  Ionic solvation studied by image-charge reaction field method. , 2011, The Journal of chemical physics.

[52]  W. Im,et al.  Generalized solvent boundary potential for computer simulations , 2001 .

[53]  Wei Cai,et al.  ICSM: An order N method for calculating electrostatic interactions added to TINKER , 2013, Comput. Phys. Commun..

[54]  Donald J. Jacobs,et al.  An image-based reaction field method for electrostatic interactions in molecular dynamics simulations of aqueous solutions. , 2009, The Journal of chemical physics.

[55]  M. Washizu,et al.  Equivalent image charges of a prolate spheroid under an external electric field , 2011 .

[56]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[57]  Stefan Ritter.  On the Computation of Lamé Functions, of Eigenvalues and Eigenfunctions of Some Potential Operators , 1998 .

[58]  WALTER GAUTSCHI.  Algorithm 726: ORTHPOL–a package of routines for generating orthogonal polynomials and Gauss-type quadrature rules , 1994, TOMS.

[59]  D. Zorin,et al.  A kernel-independent adaptive fast multipole algorithm in two and three dimensions , 2004 .

[60]  Carlos Simmerling,et al.  Chapter 6 Hybrid Explicit/Implicit Solvation Methods , 2006 .

[61]  R. Levy,et al.  Computer simulations with explicit solvent: recent progress in the thermodynamic decomposition of free energies and in modeling electrostatic effects. , 1998, Annual review of physical chemistry.

[62]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[63]  L. Greengard,et al.  Regular Article: A Fast Adaptive Multipole Algorithm in Three Dimensions , 1999 .

[64]  L. Greengard,et al.  A new version of the Fast Multipole Method for the Laplace equation in three dimensions , 1997, Acta Numerica.

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

[66]  Harold L. Friedman,et al.  Image approximation to the reaction field , 1975 .