Heterogeneous behavior of metalloproteins toward metal ion binding and selectivity: insights from molecular dynamics studies

About one-third of the existing proteins require metal ions as cofactors for their catalytic activities and structural complexities. While many of them bind only to a specific metal, others bind to multiple (different) metal ions. However, the exact mechanism of their metal preference has not been deduced to clarity. In this study, we used molecular dynamics (MD) simulations to investigate whether a cognate metal (bound to the structure) can be replaced with other similar metal ions. We have chosen seven different proteins (phospholipase A₂, sucrose phosphatase, pyrazinamidase, cysteine dioxygenase (CDO), plastocyanin, monoclonal anti-CD4 antibody Q425, and synaptotagmin 1 C₂B domain) bound to seven different divalent metal ions (Ca²⁺, Mg²⁺, Zn²⁺, Fe²⁺, Cu²⁺, Ba²⁺, and Sr²⁺, respectively). In total, 49 MD simulations each of 50 ns were performed and each trajectory was analyzed independently. Results demonstrate that in some cases, cognate metal ions can be exchanged with similar metal ions. On the contrary, some proteins show binding affinity specifically to their cognate metal ions. Surprisingly, two proteins CDO and plastocyanin which are known to bind Fe²⁺ and Cu²⁺, respectively, do not exhibit binding affinity to any metal ion. Furthermore, the study reveals that in some cases, the active site topology remains rigid even without cognate metals, whereas, some require them for their active site stability. Thus, it will be interesting to experimentally verify the accuracy of these observations obtained computationally. Moreover, the study can help in designing novel active sites for proteins to sequester metal ions particularly of toxic nature.