Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models

Arkajyoti Sengupta,Zhen Li,Lin Frank Song,Pengfei Li,Kenneth M. Merz

DOI: https://doi.org/10.1021/acs.jcim.0c01390

IF: 6.162

2021-02-04

Journal of Chemical Information and Modeling

Abstract:Monovalent ions play significant roles in various biological and material systems. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB), with significantly improved descriptions of condensed phase water, have been developed. The pairwise interaction between the metal ion and water necessitates the development of ion parameters specifically for these water models. Herein, we parameterized the 12-6 and the 12-6-4 nonbonded models for 12 monovalent ions with the respective four new water models. These monovalent ions contain eight cations including alkali metal ions (Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺), transition-metal ions (Cu⁺ and Ag⁺), and Tl⁺ from the boron family, along with four halide anions (F^-, Cl^-, Br^-, I^-). Our parameters were designed to reproduce the target hydration free energies (the 12-6 hydration free energy (HFE) set), the ion-oxygen distances (the 12-6 ion-oxygen distance (IOD) set), or both of them (the 12-6-4 set). The 12-6-4 parameter set provides highly accurate structural features overcoming the limitations of the routinely used 12-6 nonbonded model for ions. Specifically, we note that the 12-6-4 parameter set is able to reproduce experimental hydration free energies within 1 kcal/mol and experimental ion-oxygen distances within 0.01 Å simultaneously. We further reproduced the experimentally determined activity derivatives for salt solutions, validating the ion parameters for simulations of ion pairs. The improved performance of the present water models over our previous parameter sets for the TIP3P, TIP4P, and SPC/E water models (Li, P. et al <i>J. Chem. Theory Comput.</i> 2015 11 1645 1657) highlights the importance of the choice of water model in conjunction with the metal ion parameter set.

chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop monovalent ion parameters applicable to four new water models (OPC3, OPC, TIP3P - FB, and TIP4P - FB). Specifically, the authors aim to parameterize 12 monovalent ions through the 12 - 6 and 12 - 6 - 4 non - bonding models in order to accurately reproduce the experimental hydration free energy (HFE) and the ion - oxygen distance (IOD) in the first solvent shell. This work is of great significance for improving the accuracy of simulating the interactions between metal ions and water molecules, especially in various systems involving monovalent ions in biology and materials science.

### Main Objectives:
1. **Parameterize 12 - 6 and 12 - 6 - 4 Non - Bonding Models**: Develop parameters for 12 monovalent ions (including 8 cations and 4 anions) under four new water models.
2. **Reproduce Experimental Data**: Ensure that the parameterized models can accurately reproduce the experimentally measured hydration free energy (HFE) and ion - oxygen distance (IOD).
3. **Improve Model Performance**: By introducing the 12 - 6 - 4 model, overcome the limitations of the traditional 12 - 6 model in simultaneously reproducing HFE and structural features (such as IOD).

### Background:
- **Importance of Monovalent Ions**: Monovalent ions play important roles in biological, geological, and materials sciences, especially in biological ion - exchange processes and RNA folding.
- **Limitations of Existing Models**: The traditional 12 - 6 Lennard - Jones (LJ) non - bonding model cannot simultaneously and accurately reproduce the experimentally measured HFE and IOD in some cases, especially for high - charge systems, because of the lack of ion - induced dipole interactions.
- **Advantages of the New Model**: The 12 - 6 - 4 model can better describe ion - water interactions by introducing an additional attraction term \(C_4\), thereby improving the accuracy and reliability of the model.

### Methods:
- **Potential Function and Combination Rules**: Use the 12 - 6 and 12 - 6 - 4 non - bonding models and combine with the Berthelot - Lorentz (LB) combination rules for parameterization.
- **Thermodynamic Integration**: Calculate the free - energy difference through the thermodynamic integration (TI) method to ensure the accuracy of parameterization.
- **Parameter Scanning**: Optimize the parameters to match the experimental data by scanning the \(R_{\text{min}}/2\) values.

### Results:
- **Performance of the 12 - 6 Model**: For Na+, K+, Rb+, and Cs+ ions, the 12 - 6 model can simultaneously reproduce the experimentally measured HFE and IOD. However, for Li+, Tl+, Cu+, and Ag+ ions, the 12 - 6 model has relatively large errors, especially in terms of hydration free energy.
- **Performance of the 12 - 6 - 4 Model**: The 12 - 6 - 4 model can more accurately reproduce HFE and IOD simultaneously, especially when dealing with high - charge systems.

### Discussion:
- **Importance of Model Selection**: Selecting appropriate water models and ion parameter sets is crucial for simulating the interactions between metal ions and water molecules.
- **Future Directions**: Further optimize the parameters to improve the applicability and transferability of the model in different environments.

Through this research, the authors provide a reliable parameter set for the accurate simulation of monovalent ions under four new water models, which is of great significance for understanding the ion behavior in complex biological and materials systems.