Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential
Yongguang Zhang,Yang Jiang,Jiarong Peng,Haiyang Zhang
DOI: https://doi.org/10.1021/acs.jcim.1c00580
IF: 6.162
2021-07-27
Journal of Chemical Information and Modeling
Abstract:Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C4 term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M2+) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/Eb, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion–oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M2+ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion–water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jcim.1c00580.HFE, IOD, and CN values for different R–ε combinations for 11 water models, performance validation of the design models, comparison of R parameters with Shannon's revised effective ion radii, metal coordination distances in simulated E. coli glyoxalase I variants, relative errors of the HFE, IOD, and CN for SPC/E, TIP4P-Ew, OPC3, OPC, TIP3P-FB, and TIP4P-FB, IE between the M2+ ion and TIP3P water, metal-binding center of a dimeric E. coli glyoxalase I variant in the crystal structure, and RMSDs of the protein backbone and metal-binding residues (PDF)Optimized LJ parameters of M2+ models for 11 water models: TIP4P-FB, TIP4P-Ew, TIP4P-D, TIP4P/2005, TIP3P-FB, TIP3P, SPC/Eb, SPC/E, OPC3, OPC, and a99SB-disp (ZIP)This article has not yet been cited by other publications.
chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems