Simulating Metal-Imidazole Complexes

Zhen Li,Subhamoy Bhowmik,Luca Sagresti,Guiseppe Brancato,Madelyn Smith,David Benson,Pengfei Li,Kenneth Merz

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

2024-05-06

Abstract:One commonly observed binding motif in metalloproteins involves the interaction between a metal ion and histidine's imidazole sidechains. Although previous imidazole-M(II) parameters established the flexibility and reliability of the 12-6-4 Lennard-Jones (LJ)-type nonbonded model by simply tuning the ligating atom’s polarizability, they have not been applied to multiple-imidazole complexes. To fill this gap, we systematically simulate multiple-imidazole complexes (ranging from one to six) for five metal ions (Co(II), Cu(II), Mn(II), Ni(II), and Zn(II)) which commonly appear in metalloproteins. Using extensive (40ns per PMF window) sampling to assemble free energy association profiles (using OPC water and standard HID imidazole charge models from AMBER) and comparing them to DFT calculations, a new set of parameters was developed to focus on energetic and geometric features of multiple-imidazole complexes. The obtained free energy profiles agree with the experimental binding free energy and DFT calculated distances. Further investigation of the simulations revealed the existence of ligand-ligand hydrogen bonds and partial π-stacking while one ligand dissociates from the complex. To validate our model, we show that we can close the thermodynamic cycle for metal-imidazole complexes with up to six imidazole molecules in the first solvation shell. Given the success in closing the thermodynamic cycles, we then used the same extended sampling method for six other metal ions (Ag(I), Ca(II), Cd(II), Cu(I), Fe(II), and Mg(II)) to obtain new parameters. Since these new parameters can reproduce the one-imidazole geometry and energy accurately, we hypothesize that they will reasonably predict the binding free energy of higher-level coordination numbers. Hence, we did not extend the analysis of these ions up to six imidazole complexes. Overall, the results shed light on metal-protein interactions by emphasizing the importance of ligand-ligand interaction and metal-π-stacking within metalloproteins.

Chemistry

What problem does this paper attempt to address?

This paper aims to solve the problem of accurately simulating the interaction between metal ions and imidazole rings, especially the simulation of multi - imidazole complexes commonly found in metalloproteins. Although previous parameter settings have demonstrated the flexibility and reliability of the 12 - 6 - 4 Lennard - Jones (LJ) type non - bonding model in regulating the polarizability of coordinated atoms, these parameters have not yet been applied to multi - imidazole complexes. To fill this gap, the authors systematically simulated the complexes formed by five common metal ions (Co(II), Cu(II), Mn(II), Ni(II), Zn(II)) in metalloproteins and different numbers of imidazole molecules (from one to six).

By assembling free - energy correlation curves through extended sampling (40 nanoseconds per PMF window) and comparing them with the results of density functional theory (DFT) calculations, a new set of parameters was developed, focusing on the energy and geometric characteristics of multi - imidazole complexes. The resulting free - energy curves are consistent with the experimental binding free - energy and the distance calculated by DFT. Further research shows that during the dissociation of one ligand from the complex, there are ligand - ligand hydrogen bonds and partial π - stacking phenomena. To verify the validity of the model, the authors showed that a thermodynamic cycle that can close the first solvent shell containing up to six imidazole molecules can be achieved.

Given the results of successfully closing the thermodynamic cycle, the authors used the same extended sampling method to simulate six other metal ions (Ag(I), Ca(II), Cd(II), Cu(I), Fe(II), Mg(II)) to obtain new parameters. Since these new parameters can accurately reproduce the geometry and energy of single - imidazole, the researchers assume that they can reasonably predict the binding free - energy of higher coordination numbers, so they did not extend the analysis of these ions to six - imidazole complexes. Overall, these results emphasize the importance of ligand - ligand interactions and metal - π stacking in metal - protein interactions.