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

Fei Xia,Fangchen Hu,Yuwei Zhang,Pengfei Li,Ruibo Wu
DOI: https://doi.org/10.26434/chemrxiv-2024-wfqmg
2024-06-12
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 dihedral coupling effect, accurate van der Waals (vdW) interactions, and the influence of strongly polarized charges on conformational balance. The new parameters 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 protein kinase.
Chemistry
What problem does this paper attempt to address?
The main focus of this paper is the development of force field parameters for complexes formed by important molecules ATP (adenosine triphosphate) and GTP (guanosine triphosphate) with magnesium ions (Mg2+) in biochemistry. The existing AMBER force field parameters are based on the methyl triphosphate model, which cannot accurately describe the conformational distribution of ATP∙Mg2+ or GTP∙Mg2+ complexes in solution with three-coordination or two-coordination. The researchers propose a new method based on the ribose triphosphate model, considering the coupling effect of conformations, accurate van der Waals interactions, and the influence of strongly polarized charges on conformational equilibrium. The new parameters can more accurately describe the C3 and C2 conformations of these complexes in solution and can be used to simulate biological molecules involving ATP and GTP, such as protein kinases. The paper also provides a detailed description of the computational methods, including quantum mechanics calculations, force field parameter fitting, correction of dihedral coupling effects, and improvement of van der Waals parameters from density functional theory calculations to enhance the accuracy of simulations.