title: Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent.
abstract: Metal ions play significant roles in biological systems. Accurate molecular dynamics (MD) simulations on these systems require a validated set of parameters. Although there are more detailed ways to model metal ions, the nonbonded model, which employs a 12-6 Lennard-Jones (LJ) term plus an electrostatic potential is still widely used in MD simulations today due to its simple form. However, LJ parameters have limited transferability due to different combining rules, various water models and diverse simulation methods. Recently, simulations employing a Particle Mesh Ewald (PME) treatment for long-range electrostatics have become more and more popular owing to their speed and accuracy. In the present work we have systematically designed LJ parameters for 24 +2 metal (M(II)) cations to reproduce different experimental properties appropriate for the Lorentz-Berthelot combining rules and PME simulations. We began by testing the transferability of currently available M(II) ion LJ parameters. The results showed that there are differences between simulations employing Ewald summation with other simulation methods and that it was necessary to design new parameters specific for PME based simulations. Employing the thermodynamic integration (TI) method and performing periodic boundary MD simulations employing PME, allowed for the systematic investigation of the LJ parameter space. Hydration free energies (HFEs), the ion-oxygen distance in the first solvation shell (IOD) and coordination numbers (CNs) were obtained for various combinations of the parameters of the LJ potential for four widely used water models (TIP3P, SPC/E, TIP4P and TIP4PEW). Results showed that the three simulated properties were highly correlated. Meanwhile, M(II) ions with the same parameters in different water models produce remarkably different HFEs but similar structural properties. It is difficult to reproduce various experimental values simultaneously because the nonbonded model underestimates the interaction between the metal ions and water molecules at short range. Moreover, the extent of underestimation increases successively for the TIP3P, SPC/E, TIP4PEWand TIP4P water models. Nonetheless, we fitted a curve to describe the relationship between ε (the well depth) and radius (Rmin/2) from experimental data on noble gases to facilitate the generation of the best possible compromise models. Hence, by targeting different experimental values, we developed three sets of parameters for M(II) cations for three different water models (TIP3P, SPC/E and TIP4PEW). These parameters we feel represent the best possible compromise that can be achieved using the nonbonded model for the ions in combination with simple water models. From a computational uncertainty analysis we estimate that the uncertainty in our computed HFEs is on the order of ±1kcal/mol. Further improvements will require more advanced non-bonded models likely with inclusion of polarization.
authors.affiliation: 2328 New Physics Building, PO Box 118435, University of Florida, Gainesville, Florida 32611-8435, .
authors.forename: Pengfei
authors.initials: P
authors.lastname: Li
authors.forename: Benjamin P
authors.initials: BP
authors.lastname: Roberts
authors.forename: Dhruva K
authors.initials: DK
authors.lastname: Chakravorty
authors.forename: Kenneth M
authors.initials: KM
authors.lastname: Merz
keywords: Ewald summation; Lennard-Jones parameters; M(II) metal ion
references.citation: Curr Opin Chem Biol. 2008 Apr;12(2):250-5
references.pmid: 18261473
references.citation: Nat Rev Microbiol. 2009 Jan;7(1):25-35
references.pmid: 19079350
references.citation: J Phys Chem B. 2008 Jul 31;112(30):9020-41
references.pmid: 18593145
references.citation: J Mol Graph Model. 2010 Sep;29(2):116-25
references.pmid: 20675161
references.citation: J Phys Chem A. 2011 May 12;115(18):4798-803
references.pmid: 21504174
references.citation: J Phys Chem B. 2010 Dec 2;114(47):15590-7
references.pmid: 21053931
references.citation: Am J Clin Nutr. 1996 May;63(5):797S-811S
references.pmid: 8615367
references.citation: J Phys Chem B. 2006 Sep 21;110(37):18553-9
references.pmid: 16970483
references.citation: Chem Rev. 1996 Nov 7;96(7):2239-2314
references.pmid: 11848828
references.citation: J Biomol NMR. 2013 Jun;56(2):125-37
references.pmid: 23609042
references.citation: FEMS Immunol Med Microbiol. 2005 Feb 1;43(2):105-14
references.pmid: 15681139
references.citation: J Am Chem Soc. 2002 Mar 6;124(9):1958-67
references.pmid: 11866609
references.citation: J Comput Chem. 2005 Dec;26(16):1668-88
references.pmid: 16200636
references.citation: Chem Biol. 2002 Jun;9(6):673-7
references.pmid: 12079778
references.citation: Proteins. 1995 Sep;23(1):12-31
references.pmid: 8539245
references.citation: Prog Biophys Mol Biol. 1997;67(2-3):217-52
references.pmid: 9446936
references.citation: J Am Chem Soc. 2003 Dec 17;125(50):15671-82
references.pmid: 14664617
references.citation: J Am Chem Soc. 2001 Jan 24;123(3):431-41
references.pmid: 11456545
references.citation: J Am Chem Soc. 2005 Apr 6;127(13):4921-9
references.pmid: 15796557
references.citation: J Chem Phys. 2012 Mar 21;136(11):114507
references.pmid: 22443777
references.citation: J Phys Chem A. 2006 Jan 19;110(2):691-9
references.pmid: 16405342
references.citation: Proteins. 2001 Nov 15;45(3):183-9
references.pmid: 11599021
references.citation: J Chem Theory Comput. 2010 Sep 14;6(9):2935-2947
references.pmid: 20856692
references.citation: J Chem Theory Comput. 2011 Feb 8;7(2):433-443
references.pmid: 21552372
references.citation: Annu Rev Phys Chem. 2012;63:401-18
references.pmid: 22404593
references.citation: J Chem Phys. 2004 May 22;120(20):9665-78
references.pmid: 15267980
references.citation: J Phys Chem B. 2011 Aug 25;115(33):10079-85
references.pmid: 21761909
references.citation: Nucleic Acids Res. 2012 Jul;40(Web Server issue):W352-7
references.pmid: 22693222
references.citation: J Chem Phys. 2011 Jan 7;134(1):014511
references.pmid: 21219011
references.citation: Biochemistry. 2012 Mar 27;51(12):2619-29
references.pmid: 22394357
references.citation: J Am Chem Soc. 2004 Feb 4;126(4):1285-9
references.pmid: 14746502
references.citation: Nature. 1998 Oct 15;395(6703):645-8
references.pmid: 9790183
references.citation: J Am Chem Soc. 2002 Mar 6;124(9):1968-76
references.pmid: 11866610
references.citation: Curr Opin Chem Biol. 1998 Apr;2(2):155-8
references.pmid: 9667942
references.citation: Inorg Chem. 2011 May 16;50(10):4572-9
references.pmid: 21495628
references.citation: J Comput Chem. 2009 Apr 30;30(6):864-72
references.pmid: 19191337
references.citation: J Am Chem Soc. 2012 Feb 22;134(7):3367-76
references.pmid: 22007899
references.citation: J Chem Theory Comput. 2012 Apr 10;8(4):1493-502
references.pmid: 26596759
references.citation: J Phys Chem B. 2008 Apr 17;112(15):4694-702
references.pmid: 18366205
references.citation: J Biol Inorg Chem. 2008 Nov;13(8):1205-18
references.pmid: 18604568
references.citation: J Chem Phys. 2012 Jul 28;137(4):044511
references.pmid: 22852635
references.citation: J Am Chem Soc. 2013 Jan 9;135(1):30-3
references.pmid: 23214972
references.citation: Chem Rev. 1996 Nov 7;96(7):2375-2434
references.pmid: 11848831
references.citation: Annu Rev Biochem. 1992;61:897-946
references.pmid: 1497326
references.citation: J Am Chem Soc. 2004 Apr 7;126(13):4167-80
references.pmid: 15053606
references.citation: J Chem Phys. 2006 Aug 7;125(5):054511
references.pmid: 16942230
references.citation: J Am Chem Soc. 2011 Dec 7;133(48):19330-3
references.pmid: 22029374
references.citation: J Comput Chem. 2007 Dec;28(16):2618-40
references.pmid: 17894371
references.citation: J Chem Phys. 2007 Dec 7;127(21):214108
references.pmid: 18067350
references.citation: Inorg Chem. 2011 Sep 5;50(17):8509-15
references.pmid: 21800865
references.citation: J Chem Phys. 2012 Mar 28;136(12):124103
references.pmid: 22462831
references.citation: J Chem Theory Comput. 2010;6(3):774-786
references.pmid: 20300554
references.citation: Annu Rev Biophys Biophys Chem. 1990;19:405-21
references.pmid: 2114117