1
|
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
|
2
|
Wang X, Clegg SL, Di Tommaso D. Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions. J Chem Phys 2022; 156:024502. [PMID: 35032987 DOI: 10.1063/5.0074970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Chemical thermodynamic models of solvent and solute activities predict the equilibrium behavior of aqueous solutions. However, these models are semi-empirical. They represent micro-scale ion and solvent behaviors controlling the macroscopic properties using small numbers of parameters whose values are obtained by fitting to activities and other partial derivatives of the Gibbs energy measured for the bulk solutions. We have conducted atomistic simulations of aqueous electrolyte solutions (MgCl2 and CaCl2) to determine the parameters of thermodynamic hydration models. We have implemented a cooperative hydration model to categorize the water molecules in electrolyte solutions into different subpopulations. The value of the electrolyte-specific parameter, k, was determined from the ion-affected subpopulation with the lowest absolute value of the free energy of removing the water molecule. The other equilibrium constant parameter, K1, associated with the first degree of hydration, was computed from the free energy of hydration of hydrated clusters. The hydration number, h, was determined from a reorientation dynamic analysis of the water subpopulations compared to bulk-like behavior. The reparameterized models [R. H. Stokes and R. H. Robinson, J. Solution Chem. 2, 173 (1973) and Balomenos et al., Fluid Phase Equilib. 243, 29 (2006)] using the computed values of the parameters lead to the osmotic coefficients of MgCl2 solutions that are consistent with measurements. Such an approach removes the dependence on the availability of experimental data and could lead to aqueous thermodynamic models capable of estimating the values of solute and solvent activities as well as thermal and volumetric properties for a wide range of compositions and concentrations.
Collapse
Affiliation(s)
- Xiangwen Wang
- Department of Chemistry, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| | - Simon L Clegg
- School of Environmental Sciences, University of East Anglia, Norwich NR4 7TJ, United Kingdom
| | - Devis Di Tommaso
- Department of Chemistry, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| |
Collapse
|
3
|
Nde J, Zhang P, Ezerski JC, Lu W, Knapp K, Wolynes PG, Cheung MS. Coarse-Grained Modeling and Molecular Dynamics Simulations of Ca 2+-Calmodulin. Front Mol Biosci 2021; 8:661322. [PMID: 34504868 PMCID: PMC8421859 DOI: 10.3389/fmolb.2021.661322] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Accepted: 07/21/2021] [Indexed: 12/21/2022] Open
Abstract
Calmodulin (CaM) is a calcium-binding protein that transduces signals to downstream proteins through target binding upon calcium binding in a time-dependent manner. Understanding the target binding process that tunes CaM’s affinity for the calcium ions (Ca2+), or vice versa, may provide insight into how Ca2+-CaM selects its target binding proteins. However, modeling of Ca2+-CaM in molecular simulations is challenging because of the gross structural changes in its central linker regions while the two lobes are relatively rigid due to tight binding of the Ca2+ to the calcium-binding loops where the loop forms a pentagonal bipyramidal coordination geometry with Ca2+. This feature that underlies the reciprocal relation between Ca2+ binding and target binding of CaM, however, has yet to be considered in the structural modeling. Here, we presented a coarse-grained model based on the Associative memory, Water mediated, Structure, and Energy Model (AWSEM) protein force field, to investigate the salient features of CaM. Particularly, we optimized the force field of CaM and that of Ca2+ ions by using its coordination chemistry in the calcium-binding loops to match with experimental observations. We presented a “community model” of CaM that is capable of sampling various conformations of CaM, incorporating various calcium-binding states, and carrying the memory of binding with various targets, which sets the foundation of the reciprocal relation of target binding and Ca2+ binding in future studies.
Collapse
Affiliation(s)
- Jules Nde
- Department of Physics, University of Houston, Houston, TX, United States.,Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Pengzhi Zhang
- Department of Physics, University of Houston, Houston, TX, United States
| | - Jacob C Ezerski
- Department of Physics, University of Houston, Houston, TX, United States
| | - Wei Lu
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Kaitlin Knapp
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| | - Margaret S Cheung
- Department of Physics, University of Houston, Houston, TX, United States.,Center for Theoretical Biological Physics, Rice University, Houston, TX, United States
| |
Collapse
|
4
|
Peng J, Zhang Y, Jiang Y, Zhang H. Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References. J Chem Inf Model 2021; 61:2981-2997. [PMID: 34080414 DOI: 10.1021/acs.jcim.1c00281] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg2+ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl2 solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg2+ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl2 solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.
Collapse
Affiliation(s)
- Jiarong Peng
- 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
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
5
|
Klein F, Cáceres D, Carrasco MA, Tapia JC, Caballero J, Alzate-Morales J, Pantano S. Coarse-Grained Parameters for Divalent Cations within the SIRAH Force Field. J Chem Inf Model 2020; 60:3935-3943. [DOI: 10.1021/acs.jcim.0c00160] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Florencia Klein
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo 11400, Uruguay
| | - Daniela Cáceres
- Escuela de Medicina, Universidad de Talca, 1 Poniente 1141, Talca 3460000, Chile
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingenierı́a, Universidad de Talca, Campus Talca, 1 Poniente No. 1141, Talca 3460000, Chile
| | - Mónica A. Carrasco
- Escuela de Medicina, Universidad de Talca, 1 Poniente 1141, Talca 3460000, Chile
| | - Juan Carlos Tapia
- Escuela de Medicina, Universidad de Talca, 1 Poniente 1141, Talca 3460000, Chile
| | - Julio Caballero
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingenierı́a, Universidad de Talca, Campus Talca, 1 Poniente No. 1141, Talca 3460000, Chile
| | - Jans Alzate-Morales
- Centro de Bioinformática, Simulación y Modelado (CBSM), Facultad de Ingenierı́a, Universidad de Talca, Campus Talca, 1 Poniente No. 1141, Talca 3460000, Chile
| | - Sergio Pantano
- Institut Pasteur de Montevideo, Mataojo 2020, Montevideo 11400, Uruguay
- Shanghai Institute for Advanced Immunochemical Studies, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| |
Collapse
|
6
|
Cai HQ, Li HJ, Wang M, Wang CC, Yi HB. Microscopic insight into the ion aggregation characteristics in aqueous MgCl2 and MgCl2–LiCl solutions: Implications for Mg2+/Li+ separation. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2018.10.059] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
7
|
He Z, Feng G, Yang B, Yang L, Liu CW, Xu HG, Xu XL, Zheng WJ, Gao YQ. Molecular dynamics simulation, ab initio calculation, and size-selected anion photoelectron spectroscopy study of initial hydration processes of calcium chloride. J Chem Phys 2018; 148:222839. [DOI: 10.1063/1.5024279] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Zhili He
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Gang Feng
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- School of Chemistry and Chemical Engineering, Chongqing University, Chongqing 401331, China
| | - Bin Yang
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Lijiang Yang
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Cheng-Wen Liu
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Hong-Guang Xu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xi-Ling Xu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wei-Jun Zheng
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory of Molecular Reaction Dynamics, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yi Qin Gao
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|