Free energy calculations on snake venom metalloproteinase BaP1

Interactions between Triterpenes and a P-I Type Snake Venom Metalloproteinase: Molecular Simulations and Experiments

Small molecule inhibitors of snake venom metalloproteinases (SVMPs) could provide a means to rapidly halt the progression of local tissue damage following viperid snake envenomations. In this study, we examine the ability of candidate compounds based on a pentacyclic triterpene skeleton to inhibit SVMPs. We leverage molecular dynamics simulations to estimate the free energies of the candidate compounds for binding to BaP1, a P-I type SVMP, and compare these results with experimental assays of proteolytic activity inhibition in a homologous enzyme (Batx-I). Both simulation and experiment suggest that betulinic acid is the most active candidate, with the simulations predicting a standard binding free energy of ΔG∘=−11.0±1.4 kcal/mol. The simulations also reveal the atomic interactions that underlie binding between the triterpenic acids and BaP1, most notably the electrostatic interaction between carboxylate groups of the compounds and the zinc cofactor of BaP1. Together, our simulations and experiments suggest that occlusion of the S1′ subsite is essential for inhibition of proteolytic activity. While all active compounds make hydrophobic contacts in the S1′ site, β-boswellic acid, with its distinct carboxylate position, does not occlude the S1′ site in simulation and exhibits negligible activity in experiment.

Unraveling the distinctive features of hemorrhagic and non-hemorrhagic snake venom metalloproteinases using molecular simulations

Snake venom metalloproteinases are important toxins that play fundamental roles during envenomation. They share a structurally similar catalytic domain, but with diverse hemorrhagic capabilities. To understand the structural basis for this difference, we build and compare two dynamical models, one for the hemorrhagic atroxlysin-I from Bothrops atrox and the other for the non-hemorraghic leucurolysin-a from Bothrops leucurus. The analysis of the extended molecular dynamics simulations shows some changes in the local structure, flexibility and surface determinants that can contribute to explain the different hemorrhagic activity of the two enzymes. In agreement with previous results, the long Ω-loop (from residue 149 to 177) has a larger mobility in the hemorrhagic protein. In addition, we find some potentially-relevant differences at the base of the S1' pocket, what may be interesting for the structure-based design of new anti-venom agents. However, the sharpest differences in the computational models of atroxlysin-I and leucurolysin-a are observed in the surface electrostatic potential around the active site region, suggesting thus that the hemorrhagic versus non-hemorrhagic activity is probably determined by protein surface determinants.

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM-GBSA (ΔGBind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.