1
|
Koca Fındık B, Jafari M, Song LF, Li Z, Aviyente V, Merz KM. Binding of Phosphate Species to Ca 2+ and Mg 2+ in Aqueous Solution. J Chem Theory Comput 2024; 20:4298-4307. [PMID: 38718258 PMCID: PMC11137831 DOI: 10.1021/acs.jctc.4c00218] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 04/27/2024] [Accepted: 05/01/2024] [Indexed: 05/29/2024]
Abstract
Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H2PO4-, HPO42-, PO43- coordinating with Ca2+ and Mg2+. First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca2+ and Mg2+ generally show similar trends at longer metal-phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation-anion interactions.
Collapse
Affiliation(s)
- Basak Koca Fındık
- Department
of Chemistry, Bogazici University, Bebek, Istanbul 34342, Turkey
| | - Majid Jafari
- Department
of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
| | - Lin Frank Song
- Biochemical
and Biophysical Systems Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Zhen Li
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Viktorya Aviyente
- Department
of Chemistry, Bogazici University, Bebek, Istanbul 34342, Turkey
| | - Kenneth M. Merz
- Department
of Biochemistry and Molecular Biology, Michigan
State University, East Lansing, Michigan 48824, United States
- Department
of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
2
|
Goswami A, Peña-Torres A, Jónsson EÖ, Egorov SA, Jónsson H. Evidence of Sharp Transitions between Octahedral and Capped Trigonal Prism States of the Solvation Shell of the Aqueous Fe 3+ Ion. J Phys Chem Lett 2024; 15:4523-4530. [PMID: 38634894 DOI: 10.1021/acs.jpclett.4c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/19/2024]
Abstract
The structure of the solvation shell of the aqueous Fe3+ ion has been a subject of controversy due to discrepancies between experiments and different levels of theory. We address this issue by performing simulations for a wide range of ion concentrations, using various potential energy functions, supplemented by density functional theory calculations of selected configurations. The solvation shell undergoes abrupt transitions between two states: a hexacoordinated octahedral (OH) state and a capped trigonal prism (CTP) state with 7-fold coordination. The lifetime of these states is dependent on concentration. In dilute FeCl3 solutions, the lifetimes of both are similar (≈1 ns). However, the lifetime of the OH state increases with ion concentration, while that of the CTP state decreases slightly. When a uniform negative background charge is used instead of explicit counterions, the lifetime of the OH state is greatly overestimated. These findings underscore the need for further experimental measurements and higher-level simulations.
Collapse
Affiliation(s)
- Amrita Goswami
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Alejandro Peña-Torres
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Elvar Ö Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| | - Sergei A Egorov
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
- Department of Chemistry, University of Virginia, Charlottesville, Virginia 22901, United States
| | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, VR-III, 107 Reykjavík, Iceland
| |
Collapse
|
3
|
Ćeranić K, Milovanović B, Petković M. Density functional theory study of crown ether-magnesium complexes: from a solvated ion to an ion trap. Phys Chem Chem Phys 2023; 25:32656-32665. [PMID: 38010878 DOI: 10.1039/d3cp03991a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Metal ion detection rests on host-guest recognition. We propose a theoretical protocol for designing an optimal trap for a desired metal cation. A host for magnesium ions was sought for among derivatives of crown ethers 12-crown-4, 15-crown-5, and 18-crown-6. Mg-crown complexes and their hydrated counterparts with water molecules bound to the cation were optimized using density functional theory. Based on specific geometric criteria, Interacting quantum atoms analysis and density functional theory-based molecular dynamics of Mg-crown complexes immersed in water, crown ethers for optimal accommodation of Mg2+ in aqueous solution were identified. Selectivity of the chosen crowns towards Na+, K+, and Ca2+ ions is addressed.
Collapse
Affiliation(s)
- Katarina Ćeranić
- Innovative Centre of the Faculty of Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Branislav Milovanović
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| | - Milena Petković
- University of Belgrade - Faculty of Physical Chemistry, Studentski trg 12-16, 11158 Belgrade, Serbia.
| |
Collapse
|
4
|
Zhao X, Meng K, Niu Y, Ming S, Rong J, Yu X, Zhang Y. Surface/interfacial transport through pores control desalination mechanisms in 2D carbon-based membranes. Phys Chem Chem Phys 2023; 25:30296-30307. [PMID: 37930335 DOI: 10.1039/d3cp03133k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2023]
Abstract
The shortage of freshwater is a critical concern for contemporary society, and reverse osmosis desalination technology has gathered considerable attention as a potential solution to this problem. It has been recognized that the desalination process involving water flow through angstrom-sized pores has tremendous potential. However, it is challenging to obtain angstrom-sized pore structures with internal mass transfer and surface/interface properties matching the application conditions. Herein, a two-dimensional (2D) zeolite-like carbon structure (Carzeo-ANG) was constructed with unique angstrom-sized pores in the zeolite structure; then, the surface/interfacial transport behavior and percolation effect of the Carzeo-ANG desalination membrane were evaluated by density functional theory (DFT) calculations and classical molecular dynamics. The first-principles calculations in density functional theory were implemented through the Vienna ab initio simulation package (VASP), which is a commercial package for the simulation of carbon-based materials. The results show that Carzeo-ANG is periodically distributed with angstrom-sized pores (effective diameter = 5.4 Å) of dodecacyclic carbon rings, which ensure structural stability while maintaining sufficient mechanical strength. The remarkable salt-ion adsorption properties and mass transfer activity combined with the reasonable density distribution and free energy barrier for water molecules endow the membrane with superior desalination ability. At the pressure of 80 MPa, the rejection efficiency of Cl- and Na+ were 100% and 96.25%, and the membrane could achieve a water flux of 132.71 L cm-2 day-1 MPa-1. Moreover, the interconnected electronic structure of Carzeo-ANG imparts a self-cleaning effect.
Collapse
Affiliation(s)
- Xiaoyang Zhao
- Faculty of Materials Science and Engineering, Kunming University of Science and Technology, Kunming, 650093, China.
| | - Kun Meng
- Faculty of Materials Science and Engineering, Kunming University of Science and Technology, Kunming, 650093, China.
| | - Yutao Niu
- Faculty of Materials Science and Engineering, Kunming University of Science and Technology, Kunming, 650093, China.
| | - Sen Ming
- Faculty of Materials Science and Engineering, Kunming University of Science and Technology, Kunming, 650093, China.
| | - Ju Rong
- Faculty of Materials Science and Engineering, Kunming University of Science and Technology, Kunming, 650093, China.
| | - Xiaohua Yu
- Faculty of Materials Science and Engineering, Kunming University of Science and Technology, Kunming, 650093, China.
| | - Yannan Zhang
- National and Local Joint Engineering Laboratory for Lithium-ion Batteries and Materials Preparation Technology, Kunming University of Science and Technology, Kunming 650093, China
| |
Collapse
|
5
|
Zheng N, Jiang W, Zhang P, Ma L, Chen J, Zhang H. Repurposing of World-Approved Drugs for Potential Inhibition against Human Carbonic Anhydrase I: A Computational Study. Int J Mol Sci 2023; 24:12619. [PMID: 37628799 PMCID: PMC10454238 DOI: 10.3390/ijms241612619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 08/06/2023] [Accepted: 08/08/2023] [Indexed: 08/27/2023] Open
Abstract
Human carbonic anhydrases (hCAs) have enzymatic activities for reversible hydration of CO2 and are acknowledged as promising targets for the treatment of various diseases. Using molecular docking and molecular dynamics simulation approaches, we hit three compounds of methyl 4-chloranyl-2-(phenylsulfonyl)-5-sulfamoyl-benzoate (84Z for short), cyclothiazide, and 2,3,5,6-tetrafluoro-4-piperidin-1-ylbenzenesulfonamide (3UG for short) from the existing hCA I inhibitors and word-approved drugs. As a Zn2+-dependent metallo-enzyme, the influence of Zn2+ ion models on the stability of metal-binding sites during MD simulations was addressed as well. MM-PBSA analysis predicted a strong binding affinity of -18, -16, and -14 kcal/mol, respectively, for these compounds, and identified key protein residues for binding. The sulfonamide moiety bound to the Zn2+ ion appeared as an essential component of hCA I inhibitors. Vina software predicted a relatively large (unreasonable) Zn2+-sulfonamide distance, although the relative binding strength was reproduced with good accuracy. The selected compounds displayed potent inhibition against other hCA isoforms of II, XIII, and XIV. This work is valuable for molecular modeling of hCAs and further design of potent inhibitors.
Collapse
Affiliation(s)
| | | | | | | | | | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
6
|
Matias PMC, Sousa JFM, Bernardino EF, Vareda JP, Durães L, Abreu PE, Marques JMC, Murtinho D, Valente AJM. Reduced Chitosan as a Strategy for Removing Copper Ions from Water. Molecules 2023; 28:molecules28104110. [PMID: 37241851 DOI: 10.3390/molecules28104110] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 05/12/2023] [Accepted: 05/12/2023] [Indexed: 05/28/2023] Open
Abstract
Toxic heavy metals are priority pollutants in wastewater, commonly present in dangerous concentrations in many places across the globe. Although in trace quantities copper is a heavy metal essential to human life, in excess it causes various diseases, whereby its removal from wastewater is a necessity. Among several reported materials, chitosan is a highly abundant, non-toxic, low-cost, biodegradable polymer, comprising free hydroxyl and amino groups, that has been directly applied as an adsorbent or chemically modified to increase its performance. Taking this into account, reduced chitosan derivatives (RCDs 1-4) were synthesised by chitosan modification with salicylaldehyde, followed by imine reduction, characterised by RMN, FTIR-ATR, TGA and SEM, and used to adsorb Cu(II) from water. A reduced chitosan (RCD3), with a moderate modification percentage (43%) and a high imine reduction percentage (98%), proved to be more efficient than the remainder RCDs and even chitosan, especially at low concentrations under the best adsorption conditions (pH 4, RS/L = 2.5 mg mL-1). RCD3 adsorption data were better described by the Langmuir-Freundlich isotherm and the pseudo-second-order kinetic models. The interaction mechanism was assessed by molecular dynamics simulations, showing that RCDs favour Cu(II) capture from water compared to chitosan, due to a greater Cu(II) interaction with the oxygen of the glucosamine ring and the neighbouring hydroxyl groups.
Collapse
Affiliation(s)
- Pedro M C Matias
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Joana F M Sousa
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Eva F Bernardino
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - João P Vareda
- University of Coimbra, CIEPQPF, Department of Chemical Engineering, 3030-790 Coimbra, Portugal
| | - Luisa Durães
- University of Coimbra, CIEPQPF, Department of Chemical Engineering, 3030-790 Coimbra, Portugal
| | - Paulo E Abreu
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Jorge M C Marques
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Dina Murtinho
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| | - Artur J M Valente
- University of Coimbra, CQC-IMS, Department of Chemistry, 3004-535 Coimbra, Portugal
| |
Collapse
|
7
|
Melse O, Antes I, Kaila VRI, Zacharias M. Benchmarking biomolecular force field-based Zn 2+ for mono- and bimetallic ligand binding sites. J Comput Chem 2023; 44:912-926. [PMID: 36495007 DOI: 10.1002/jcc.27052] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Revised: 11/23/2022] [Accepted: 11/24/2022] [Indexed: 12/14/2022]
Abstract
Zn2+ is one of the most versatile biologically available metal ions, but accurate modeling of Zn2+ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn2+ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn2+ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn2+ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn2+ site geometries that could guide the further development of biomolecular Zn2+ models.
Collapse
Affiliation(s)
- Okke Melse
- Center for Functional Protein Assemblies (CPA), Technical University of Munich, Garching, Germany.,SynBiofoundry@TUM, Technical University of Munich, Straubing, Germany
| | - Iris Antes
- Center for Functional Protein Assemblies (CPA), Technical University of Munich, Garching, Germany.,SynBiofoundry@TUM, Technical University of Munich, Straubing, Germany
| | - Ville R I Kaila
- Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | - Martin Zacharias
- Center for Functional Protein Assemblies (CPA), Technical University of Munich, Garching, Germany
| |
Collapse
|
8
|
Sinha S, Pindi C, Ahsan M, Arantes PR, Palermo G. Machines on Genes through the Computational Microscope. J Chem Theory Comput 2023; 19:1945-1964. [PMID: 36947696 PMCID: PMC10104023 DOI: 10.1021/acs.jctc.2c01313] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/24/2023]
Abstract
Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.
Collapse
|
9
|
Li Z, Song LF, Sharma G, Koca Fındık B, Merz KM. Accurate Metal-Imidazole Interactions. J Chem Theory Comput 2022; 19:619-625. [PMID: 36584400 DOI: 10.1021/acs.jctc.2c01081] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Modeling the interaction between a metal ion and small molecules can provide pivotal information to bridge and close the gap between two types of simulations: metal ions in water and metal ions in metalloproteins. As previously established, the 12-6-4 Lennard-Jones (LJ)-type nonbonded model, because of its ability to account for the induced dipole effect, has been highly successful in simulating metal ion systems. Using the potential of mean force (PMF) method, the polarizability of the metal-chelating nitrogen from two types of imidazole molecules, delta nitrogen protonated (HID) and epsilon nitrogen protonated (HIE), has been parametrized against experiment for 11 metals (Ag(I), Ca(II), Cd(II), Co(II), Cu(I), Cu(II), Fe(II), Mg(II), Mn(II), Ni(II), and Zn(II)) in conjunction with three commonly used water models (TIP3P, SPC/E, and OPC). We show that the standard 12-6 and unmodified 12-6-4 models are not able to accurately model these interactions and, indeed, predict that the complex should be unstable. The resultant parameters further establish the flexibility and the reliability of the 12-6-4 LJ-type nonbonded model, which can correctly describe three-component interactions between a metal, ligand, and solvent by simply tuning the polarizability of the chelating atom. Also, the transferability of this model was tested, showing the capability of describing metal-ligand interactions in various environments.
Collapse
Affiliation(s)
- Zhen Li
- Department of Chemistry, Michigan State University, East Lansing, Michigan48824, United States
| | - Lin Frank Song
- Department of Chemistry, Michigan State University, East Lansing, Michigan48824, United States
| | - Gaurav Sharma
- Department of Chemistry, Michigan State University, East Lansing, Michigan48824, United States
| | - Basak Koca Fındık
- Department of Chemistry, Boğaziçi University, Bebek, 34342Istanbul, Turkey
| | - Kenneth M Merz
- Department of Chemistry, Michigan State University, East Lansing, Michigan48824, United States
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan48824, United States
| |
Collapse
|
10
|
Theoretical investigation on the structure and physicochemical properties of choline chloride-based deep eutectic solvents. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
11
|
Fan K, Zhang Y, Qiu Y, Zhang H. Impacts of targeting different hydration free energy references on the development of ion potentials. Phys Chem Chem Phys 2022; 24:16244-16262. [PMID: 35758314 DOI: 10.1039/d2cp01237e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Hydration free energy (HFE) as the most important solvation parameter is often targeted in ion model development, even though the reported values differ by dozens of kcal mol-1 mainly due to the experimentally undetermined HFE of the proton ΔG°(H+). The choice of ΔG°(H+) obviously affects the hydration of single ions and the relative HFE between the ions with different (magnitude or sign) charges, and the impacts of targeted HFEs on the ion solvation and ion-ion interactions are largely unrevealed. Here we designed point charge models of K+, Mg2+, Al3+, and Cl- ions targeting a variety of HFE references and then investigated the HFE influences on the simulations of dilute and concentrated ion solutions and of the salt ion pairs in gas, liquid, and solid phases. Targeting one more property of ion-water oxygen distances (IOD) leaves the ion-water binding distance invariant, while the binding strength increases with the decreasing (more negative) HFE of ions as a result of a decrease in ΔG°(H+) for the cation and an increase in ΔG°(H+) for the anion. The increase in ΔG°(H+) leads to strengthened cation-anion interactions and thus to close ion-ion contacts, low osmotic pressures, and small activity derivatives in concentrated ion solutions as well as too stable ion pairs of the salts in different phases. The ion diffusivity and water exchange rates around the ions are simply not HFE dependent but rather more complex. Targeting both the aqueous IOD and salt crystal properties of KCl was also attempted and the comparison between different models indicates the complexity and challenge in obtaining a balanced performance between different phases using classical force fields. Our results also support that a real ΔG°(H+) value of -259.8 kcal mol-1 recommended by Hünenberger and Reif guides ion models to reproduce ion-water and ion-ion interactions reasonably at relatively low salt concentrations. Simulations of a metalloprotein show that a relatively more positive ΔG°(H+) for Mg2+ model is better for a reasonable description of the metal binding network.
Collapse
Affiliation(s)
- Kun Fan
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China.
| |
Collapse
|
12
|
Qiu Y, Jiang Y, Zhang Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential. J Phys Chem B 2021; 125:13502-13518. [PMID: 34860517 DOI: 10.1021/acs.jpcb.1c09103] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C4 term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C4 term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C4 term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.
Collapse
Affiliation(s)
- Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
13
|
Zhang Y, Jiang Y, Qiu Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Highly Charged Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4613-4629. [PMID: 34467756 DOI: 10.1021/acs.jcim.1c00723] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Here, we developed nonbonded point charge models using a simple Lennard-Jones (LJ) 12-6 potential for highly charged metal cations (18 trivalent and 6 tetravalent ions) for use with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models simultaneously reproduce the hydration free energy (HFE) and ion-oxygen distance (IOD) in the first hydration shell with an error within 1 kcal/mol and 0.01 Å on average, respectively, and yield reasonable coordination numbers for most cations. Such performance is equivalent to the previously reported point charge models using a more complex 12-6-4 LJ-type potential, while the LJ R parameters of our models are much close to Shannon's revised effective ion radii than that of the 12-6-4 models. Our designed models overestimate the diffusion constants of several trivalent ions by 5-68%. The performance in predicting osmotic coefficients of trivalent chlorides in aqueous solution depends on the salt type. A calibration of cation-anion interacting LJ parameters reproduces the experimental osmotic coefficients of an AlCl3 solution at 0.2-3.0 mol/L. The effectiveness of our new models is further demonstrated by simulating a metalloprotein system with four force field/water combinations. This work facilitates accurate modeling of metal-containing systems by a variety of force fields and water models in aqueous solution.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|