The design of Artificial Metalloenzymes (ArMs), which result from the incorporation of organometallic cofactors into biological structures, has grown steadily in the last two decades and important new-to-Nature reactions have been reached. These type of exercises could greatly benefit from an understanding of the structural impact that the inclusion of organometallic moieties may have on the biological host. To date though, our understanding of this phenomenon is highly partial. This lack of knowledge is one of the elements that condition that first-generation ArMs generally display relatively poor catalytic profiles. In this work, we approach this matter by assessing the dynamics and stability of a series of ArMs resulting from the inclusion, via different anchoring strategies, of a variety of organometallic cofactors into the Lactococcal multidrug resistance regulator (LmrR) protein. To this aim, we coupled standard force field-based techniques such as Protein-Ligand Docking and Molecular Dynamics simulations with a variety of trajectory convergence analyses, capable of assessing both the stability and flexibility of the different systems under study upon the binding of cofactors. Together with the experimental evidence obtained in other studies, we provide an overview on how these changes can affect the catalytic outcomes obtained from the different ArMs. Fundamentally, our results show that the convergence analysis used in this work can assess how the inclusion of synthetic metallic cofactors in proteins can condition different structural modulations of their host. Those conformational modifications are key to the success of the desired catalytic activity and their proper identification can be wisely used to improve the quality and the rate of success of the ArMs.
Keywords: artificial metalloenzymes; cofactor binding; interactive analysis; molecular dynamics; molecular modeling; molecular plasticity.