Getting Deeper into the Molecular Events of Heme Binding Mechanisms: A Comparative Multi-level Computational Study of HasAsm and HasAyp Hemophores

[Image: see text] Many biological systems obtain their activity by the inclusion of metalloporphyrins into one or several binding pockets. However, decoding the molecular mechanism under which these compounds bind to their receptors is something that has not been widely explored and is a field with open questions. In the present work, we apply computational techniques to unravel and compare the mechanisms of two heme-binding systems, concretely the HasA hemophores from Gram negative bacteria Serratiamarcescens (HasAsm) and Yersinia pestis (HasAyp). Despite the high sequence identity between both systems, the comparison between the X-ray structures of their apo and holo forms suggests different heme-binding mechanisms. HasAyp has extremely similar structures for heme-free and heme-bound forms, while HasAsm presents a very large displacement of a loop that ultimately leads to an additional coordination to the metal with respect to HasAyp. We combined Gaussian accelerated molecular dynamics simulations (GaMDs) in explicit solvent and protein–ligand docking optimized for metalloligands. GaMDs were first carried out on heme-free forms of both hemophores. Then, protein–ligand dockings of the heme were performed on cluster representatives of these simulations and the best poses were then subjected to a new series of GaMDs. A series of analyses reveal the following: (1) HasAyp has a conformational landscape extremely similar between heme-bound and unbound states with no to limited impact on the binding of the cofactor, (2) HasAsm presents as a slightly broader conformational landscape in its apo state but can only visit conformations similar to the X-ray of the holo form when the heme has been bound. Such behavior results from a complex cascade of changes in interactions that spread from the heme-binding pocket to the flexible loop previously mentioned. This study sheds light on the diversity of molecular mechanisms of heme-binding and discusses the weight between the pre-organization of the receptor as well as the induced motions resulting in association.