4.7 Article

Systematic Parametrization of Divalent Metal Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Journal

JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 16, Issue 7, Pages 4429-4442

Publisher

AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.0c00194

Keywords

-

Funding

  1. National Institutes of Health [GM056207, GM044974, GM066859]

Ask authors/readers for more resources

Divalent metal ions play important roles in biological and materials systems. Molecular dynamics simulation is an efficient tool to investigate these systems at the microscopic level. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB) have been developed and better represent the physical properties of water than previous models. Metal ion parameters are dependent on the water model employed, making it necessary to develop metal ion parameters for select new water models. In the present work, we performed parameter scanning for the 12-6 Lennard-Jones nonbonded model of divalent metal ions in conjunction with the four new water models as well as four previous water models (TIP3P, SPC/E, TIP4P, and TIP4P-Ew). We found that these new three-point and four-point water models provide comparable or significantly improved performance for the simulation of divalent metal ions when compared to previous water models in the same category. Among all eight water models, the OPC3 water model yields the best performance for the simulation of divalent metal ions in the aqueous phase when using the 12-6 model. On the basis of the scanning results, we independently parametrized the 12-6 model for 24 divalent metal ions with each of the four new water models. As noted previously, the 12-6 model still fails to simultaneously reproduce the experimental hydration free energy (HFE) and ion-oxygen distance (IOD) values even with these new water models. To solve this problem, we parametrized the 12-6-4 model for the 16 divalent metal ions for which we have both experimental HFE and IOD values for each of the four new water models. The final parameters are able to reproduce both the experimental HFE and IOD values accurately. To validate the transferability of our parameters, we carried out benchmark calculations to predict the energies and geometries of ion-water clusters as well as the ion diffusivity coefficient of Mg2+. By comparison to quantum chemical calculations and experimental data, these results show that our parameters are well designed and have excellent transferability. The metal ion parameters for the 12-6 and 12-6-4 models reported herein can be employed in simulations of various biological and materials systems when using the OPC3, OPC, TIP3P-FB, or TIP4P-FB water model.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Chemistry, Physical

Systematic Evaluation of Ion Diffusion and Water Exchange

Zhen Li, Kenneth M. Merz

Summary: Diffusion, as a fundamental property of fluids, plays a crucial role in both scientific research and everyday life. While the diffusion properties of many liquids, such as water, have been extensively studied, there has been a lack of theoretical research on transition metal ions. In this study, a workflow called ISAIAH was developed to accurately simulate the diffusion coefficients of 15 monoatomic ions in water. The results showed good agreement with experimental data, which allowed for the theoretical prediction of the diffusion coefficient of 239Pu4+ ion. The use of an augmented 12-6-4 Lennard-Jones potential parameter set showed better results compared to some other parameter sets, indicating that both hydration free energy and ion-oxygen distance play important roles in diffusion.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2022)

Article Chemistry, Physical

Mechanism of Zinc Transport through the Zinc Transporter YiiP

Gaurav Sharma, Kenneth M. Merz

Summary: Zinc is an essential transition metal ion that plays multiple roles in cellular function. This study used molecular dynamics simulations to explore the structure-function relationship of YiiP protein in a lipid bilayer and proposed a hypothesis on the zinc efflux mechanism. The study revealed the important role of zinc ions in restraining the protein's transmembrane domains and suggested the guiding role of amino acid H153 in transporting zinc ions. Understanding the zinc transport process can provide insights into treating diseases such as diabetes and cancer.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2022)

Article Biochemical Research Methods

In Silico Collision Cross Section Calculations to Aid Metabolite Annotation

Susanta Das, Kiyoto Aramis Tanemura, Laleh Dinpazhoh, Mithony Keng, Christina Schumm, Lydia Leahy, Carter K. Asef, Markace Rainey, Arthur S. Edison, Facundo M. Fernandez, Kenneth M. Merz

Summary: In this study, an efficient CCS computational workflow was developed to accurately predict unknown structures using a machine learning model and standard DFT methods. TWIMS experiments were performed to validate the experimental values and assess uncertainties. The workflow yielded accurate structural predictions and provided unique insights into preferred conformations.

JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY (2022)

Article Biochemistry & Molecular Biology

Computer-aided drug design, quantum-mechanical methods for biological problems

Madushanka Manathunga, Andreas W. Gotz, Kenneth M. Merz

Summary: Quantum chemistry allows for the study of systems with high accuracy, but its computational expense limits its application to a few atoms. To optimize performance, methods such as implementing quantum mechanical methods on modern hardware and using multiscale approaches have been developed. Simplified QM methods, including the use of machine learning, have also been created to enhance the speed and accuracy of calculations.

CURRENT OPINION IN STRUCTURAL BIOLOGY (2022)

Article Chemistry, Medicinal

Quantum Mechanics/Molecular Mechanics Simulations on NVIDIA and AMD Graphics Processing Units

Madushanka Manathunga, Hasan Metin Aktulga, Andreas W. Gotz, Jr Kenneth M. Merz

Summary: The researchers have ported and optimized the GPU-accelerated QUICK and AMBER-based QM/MM implementation on AMD GPUs. This includes the entire Fock matrix build and force calculation in QUICK, along with general performance improvements to the QUICK GPU code. Benchmark tests on NVIDIA V100 and AMD MI100 cards show similar performance for standalone HF/DFT calculations and QM/MM molecular dynamics simulations with QUICK and QUICK/AMBER. In addition, significant speedups are observed for QM/MM molecular dynamics simulations compared to the previous release version.

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2023)

Article Chemistry, Medicinal

Molecular Insights into the Calcium Binding in Troponin C through a Molecular Dynamics Study

Madelyn Smith, Pengfei Li

Summary: Calcium-binding proteins play critical roles in various biological processes. Understanding the ion binding mechanism of EF-hand motif in calcium-binding proteins is important for inhibitor design. This study systematically investigates the ion binding affinity of four calcium-binding peptides containing the EF-hand loop. The results provide molecular insights into the ion binding and the importance of the E-helix motif in this regulation.

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2023)

Editorial Material Chemistry, Medicinal

Editorial: Machine Learning in Bio-cheminformatics

Kenneth M. Merz, Guo-Wei Wei, Feng Zhu

JOURNAL OF CHEMICAL INFORMATION AND MODELING (2023)

Article Chemistry, Physical

AMBER Free Energy Tools: A New Framework for the Design of Optimized Alchemical Transformation Pathways

Hsu-Chun Tsai, Tai-Sung Lee, Abir Ganguly, Timothy J. Giese, Maximilian C. C. J. C. Ebert, Paul Labute, Kenneth M. Merz, Darrin M. York

Summary: We propose a framework for optimized alchemical transformation pathways in free energy simulations using nonlinear mixing and a new functional form for softcore potentials. The framework is implemented and tested in the GPU-accelerated AMBER software suite. The optimized pathways integrate important features such as smoothstep functions, power scaling of interactions, LJ pairwise form, and smoothing of the potential at the nonbonded cutoff boundary. The pathways demonstrate superior numerical stability and minimal variance of free energy estimates compared to traditional methods.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Multidisciplinary

An Inorganic Fluorescent Chemosensor: Rational Design and Selective Mg2+ Detection

Theodore Gerard, Yang Wei, Erwin Weerawardhana, Adriana Lugosan, Matthias Zeller, Diane A. Dickie, Pengfei Li, Wei-Tsung Lee

Summary: A Zn2+ based complex, 3, exhibits enhanced fluorescence emission in the presence of Mg2+. Both fluorescent and computational studies reveal that 3 selectively interacts with Mg2+ due to the optimal cavity size formed by two uncoordinated pyrazole side arms. This research presents a novel approach to the development of fluorescent chemosensors.

ACS OMEGA (2023)

Article Biochemistry & Molecular Biology

13C Electron Nuclear Double Resonance Spectroscopy-Guided Molecular Dynamics Computations Reveal the Structure of the Enzyme-Substrate Complex of an Active, N-Linked Glycosylated Lipoxygenase

Ajay Sharma, Chris Whittington, Mohammed Jabed, S. Gage Hill, Anastasiia Kostenko, Tao Yu, Pengfei Li, Peter E. Doan, Brian M. Hoffman, Adam R. Offenbacher

Summary: We report the active-site structure of the MoLOX substrate complex, which provides insight into reactivity differences across the LOX family and can guide the development of MoLOX inhibitors. The ENDOR-guided MD approach used in this study demonstrates its robustness in describing LOX-substrate structures.

BIOCHEMISTRY (2023)

Article Chemistry, Physical

Consequences of Overfitting the van der Waals Radii of Ions

Madelyn Smith, Zhen Li, Luke Landry, Kenneth M. Merz, Pengfei Li

Summary: Atomic radii are important in scientific research and can be derived from crystal structures. However, determining the van der Waals (VDW) radii of ions is a challenge due to the combination of VDW and electrostatic interactions in crystal structures. In this study, VDW radii were determined based on a wavefunction analysis, showing excellent agreement with experimentally determined VDW radii of noble-gas atoms.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Multidisciplinary

Understanding the Key Roles of pH Buffer in Accelerating Lignin Degradation by Lignin Peroxidase

Wenhan Fang, Shishi Feng, Zhihui Jiang, Wanzhen Liang, Pengfei Li, Binju Wang

Summary: pH buffer plays a critical role in accelerating the degradation of lignin substrate in lignin peroxidase (LiP). The pH buffer of tartaric acid stabilizes the Trp171H·+ cation radical and enhances the oxidizing power of the enzyme, facilitating the thermodynamics of the reaction and increasing the rate of lignin degradation. These findings provide valuable insights into pH-dependent redox reactions and tryptophan-mediated biological electron transfer reactions.

JACS AU (2023)

No Data Available