Parameterization of a Fluctuating Charge Model for Complexes Containing 3d Transition Metals

Pengfei Li,Luke Landry

DOI: https://doi.org/10.26434/chemrxiv-2024-lpsvj

2024-05-16

Abstract:Metalloproteins widely exist in biology, playing pivotal roles in diverse life processes. Meanwhile, molecular dynamics (MD) simulations, which are driven by classical force fields, have emerged as indispensable tools in scientific research. Among the critical parameters within classical force fields are partial charges, which are traditionally derived from quantum mechanical (QM) calculations. However, QM calculations are often time-insensitive and prone to basis set dependence. Alternatively, fluctuating charge (FQ) models offer a promising avenue for partial charge derivation, boasting significant speed advantages conducive to large-scale screening. Building upon our previous work, which tailored an FQ model for zinc-containing complexes, we have extended this model to include additional 3d transition metals integral to the life sciences, namely chromium, manganese, iron, cobalt, and nickel. Employing CM5 charges as reference points for parameterization, our FQ model accurately reproduces partial charges for 3d metal complexes featuring biologically relevant ligands. Furthermore, by using atomic charges derived by our FQ model, MD simulations have been performed, and the results show excellent performance in simulating proteomic metal sites housing multiple metal ions. Specifically, our simulations accurately capture the behavior of a metalloprotein containing an iron-sulfur cluster and another containing a di-manganese metal site. These results showcase comparable performance to those of RESP charges. We anticipate that our study can accelerate the parameterization of atomic charges, thereby facilitating simulations of metalloproteins featuring 3d transition metals.

Chemistry

What problem does this paper attempt to address?

The main goal of this paper is to address the issue of partial charge parameterization in molecular dynamics (MD) simulations related to 3d transition metals (such as chromium, manganese, iron, cobalt, and nickel). Specifically, the authors aim to develop a fluctuating charge (FQ) model to improve the accuracy and efficiency of simulating these metal ion complexes in biological systems.

### Problems the Paper Attempts to Solve:
1. **Limitations of Traditional Quantum Mechanics Methods**: Traditional quantum mechanics (QM) computational methods, while capable of generating partial charge parameters, are time-consuming and susceptible to basis set dependency.
2. **Need for Efficient Computational Methods**: There is a need for a more computationally efficient method to replace traditional quantum mechanics calculations for application in large-scale screening processes.
3. **Extending Previous Research**: Building on previous research, the goal is to extend the FQ model, which was applicable to zinc complexes, to be applicable to a wider variety of 3d transition metal complexes.

### Main Contributions:
- **Extension of the FQ Model**: The FQ model is extended to include 3d transition metal complexes such as chromium, manganese, iron, cobalt, and nickel.
- **Accurate Partial Charge Prediction**: By parameterizing with reference to CM5 charges, the FQ model can accurately predict the partial charges of these metal complexes.
- **Application in MD Simulations**: The atomic charges obtained from the FQ model perform well in MD simulations, especially in simulations of protein metal sites containing multiple metal ions.

Through this work, the study aims to accelerate the parameterization process of 3d transition metal complexes and promote molecular modeling of related biological systems.