ABSTRACT
Snakebite envenoming affects millions of people worldwide, being officially considered a neglected tropical disease by the World Health Organization. The antivenom is effective in neutralizing the systemic effects of envenomation, but local effects are poorly neutralized, often leading to permanent disability. The natural resistance of the South American pit viper Bothrops jararaca to its venom is partly attributed to BJ46a, a natural snake venom metalloendopeptidase inhibitor. Upon complex formation, BJ46a binds non-covalently to the metalloendopeptidase, rendering it unable to exert its proteolytic activity. However, the structural features that govern this interaction are largely unknown. In this work, we applied structural mass spectrometry techniques (cross-linking-MS and hydrogen-deuterium exchange MS) and in silico analyses (molecular modeling, docking, and dynamics simulations) to understand the interaction between BJ46a and jararhagin, a metalloendopeptidase from B. jararaca venom. We explored the distance restraints generated from XL-MS experiments to guide the modeling of BJ46a and jararhagin, as well as the protein-protein docking simulations. HDX-MS data pinpointed regions of protection/deprotection at the interface of the BJ46a-jararhagin complex which, in addition to the molecular dynamics simulation data, reinforced our proposed interaction model. Ultimately, the structural understanding of snake venom metalloendopeptidases inhibition by BJ46a could lead to the rational design of drugs to improve anti-snake venom therapeutics, alleviating the high morbidity rates currently observed.