Abstract:In cells, adenosine triphosphate (ATP) and guanosine triphosphate (GTP) molecules typically form tri-coordinated or bi-coordinated ATP∙Mg2+ or GTP∙Mg2+ complexes with Mg2+ ions and bind to proteins, participating in and regulating many important cellular functions. The accuracy of their force field parameters plays a crucial role in studying the function-related conformations of ATP∙Mg2+ or GTP∙Mg2+ using molecular dynamics (MD) simulations. The parameters developed based on the methyl triphosphate model in existing AMBER force fields cannot accurately describe the conformational distribution of tri-coordinated or bi-coordinated ATP∙Mg2+ or GTP∙Mg2+ complexes in solution. In this study, we develop force field parameters for the triphosphate group based on the new ribosyl triphosphate model, considering the dihedral coupling effect, accurate van der Waals (vdW) interactions, and the influence of strongly polarized charges on conformational balance. The new force fields can accurately describe the conformational balance of tri-coordinated and bi-coordinated ATP∙Mg2+ or GTP∙Mg2+ conformations in solution and can be applied to simulate biological systems containing ATP∙Mg2+ or GTP∙Mg2+ complexes.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is that in molecular dynamics (MD) simulations, the existing force - field parameters cannot accurately describe the conformational distribution of ATP·Mg²⁺ and GTP·Mg²⁺ complexes in solution, especially the equilibrium between the tricoordinate (C3) and dicoordinate (C2) conformations. Specifically: 1. **Limitations of existing models**: The existing AMBER force - field parameters were developed based on the methyl triphosphate (MTP) model and failed to fully consider the influence of the ribose ring on the conformation of the phosphate group, resulting in an inability to accurately describe the conformational distribution of ATP·Mg²⁺ and GTP·Mg²⁺ complexes. 2. **Conformational energy differences**: There are significant differences in the free - energy differences (ΔG_C23) of the ATP·Mg²⁺ complex from the C3 to C2 conformational transition obtained by different research groups using different force - field parameters. For example, the ΔG_C23 predicted by the AMBER force - field is approximately 7.8 kcal/mol, while that predicted by the CHARMM22 force - field is approximately - 5.0 kcal/mol. These differences indicate the inconsistency of the existing force - field parameters. 3. **Accuracy of non - bonded interactions**: The existing non - bonded interaction parameters (such as van der Waals (vdW) interactions and charge distributions) may not be sufficient to accurately describe the interactions between Mg²⁺ ions and phosphate oxygen atoms. To overcome these problems, the authors propose the following improvement measures: 1. **New ribose triphosphate (RTP) model**: A RTP model including the ribose ring was introduced to more accurately describe the conformational distribution of the phosphate group. 2. **Consideration of dihedral - angle coupling effects**: Through the 2D CMAP method, the coupling effects between dihedral angles were considered to improve the accuracy of conformational description. 3. **Accurate vdW parameters**: Based on high - precision density functional theory (DFT) calculations, the vdW interaction parameters between Mg²⁺ ions and phosphate oxygen atoms were re - parameterized. 4. **Optimization of charge distribution**: The RESP charge distribution of ATP molecules was recalculated using the B3LYP - D3BJ/6 - 311+G** method combined with the PCM model. Through these improvements, the authors aim to develop new force - field parameters that can more accurately describe the conformational equilibrium of ATP·Mg²⁺ and GTP·Mg²⁺ complexes in solution, thereby providing a more reliable basis for the simulation of biological systems.

Developing Force Fields for Accurate Mg2+ and Triphosphate Interactions in ATP∙Mg2+ and GTP∙Mg2+ Complexes

Developing Force Fields for Accurate Mg2+ and Triphosphate Interactions in ATP∙Mg2+ Complexes

ATP–Magnesium Coordination: Protein Structure-Based Force Field Evaluation and Corrections

Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties

Extensive Evaluation of Force Fields for G-Quadruplexes

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties in SPC/E, TIP3P-fb, TIP4P/2005, TIP4P-Ew, and TIP4P-D

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations.

Polarizable Multipole-Based Force Field for Dimethyl and Trimethyl Phosphate

ATP and its N6-substituted analogues: parameterization, molecular dynamics simulation and conformational analysis

Phosphorylation Modification Force Field FB18CMAP Improving Conformation Sampling of Phosphoproteins

Force fields for monovalent and divalent metal cations in TIP3P water based on thermodynamic and kinetic properties

Role of Mg2+ Ions in Protein Kinase Phosphorylation: Insights from Molecular Dynamics Simulations of ATP-kinase Complexes

Study of interactions between metal ions and protein model compounds by energy decomposition analyses and the AMOEBA force field.

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA–Protein Complexes

Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Coupled Two-Dimensional Main-Chain Torsional Potential for Protein Dynamics II: Performance and Validation.

A New Quantum Calibrated Force Field for Zinc-Protein Complex

A Multi-State Coarse-Grained Simulation Model Captures Conformational Cycling in P-Type ATPases

A consistent picture of phosphate–divalent cation binding from models with implicit and explicit electronic polarization