Structure-Based Discovery of Thiosemicarbazone ... - ACS Publications

Sep 26, 2017 - Instituto de Genética e Bioquímica, UFU, Uberlândia, MG, Brazil. •S Supporting Information. ABSTRACT: The venoms of snakes are com...
0 downloads 0 Views 4MB Size
Letter Cite This: ACS Med. Chem. Lett. 2017, 8, 1136-1141

pubs.acs.org/acsmedchemlett

Structure-Based Discovery of Thiosemicarbazone Metalloproteinase Inhibitors for Hemorrhage Treatment in Snakebites Francis B. Ferreira,§ Thiago M. Pereira,§ Dayane L. N. Souza,‡ Daiana S. Lopes,‡ Vitor Freitas,‡ Veridiana M. R. Á vila,‡ Arthur E. Kümmerle,§ and Carlos Mauricio R. Sant’Anna*,§ §

Departamento de Química, Instituto de Ciências Exatas, UFRRJ, Seropédica, RJ, Brazil Instituto de Genética e Bioquímica, UFU, Uberlândia, MG, Brazil



S Supporting Information *

ABSTRACT: The venoms of snakes are composed by many toxins, which are responsible for various toxic effects including intense pain, bleeding disorders, and local tissue damage caused by hemorrhage and necrosis. The snake venom metalloproteinases (SVMPs) are proteolytic zinc-dependent enzymes acting in different hemostatic mechanisms. In this work, a structure-based molecular modeling strategy was used for the rational design, by means of a homology 3D model of an SVMP isolated from Bothrops pauloensis venom (BpMP-I), followed by synthesis and in vitro evaluation of new thiosemicarbazones as the first inhibitors of the B. pauloensis SVMP. Besides being effective for the SVMP inhibition, two molecules were shown to be effective also in vivo, inhibiting hemorrhage caused by the B. pauloensis whole venom. Docking studies on metalloproteinases from other snake species suggest that the thiosemicarbazones activity is not confined to BpMP-I, but seems to be a common feature of metzincins. KEYWORDS: Snakebites, metalloproteinase inhibitors, thiosemicarbazones, molecular docking

S

main effect induced by SVMPs is hemorrhage, an effect caused by degradation and disruption of the capillary basement membrane leading to edema, myonecrosis, and reduced ability to regenerate muscle tissue.9 SVMPs have been exhaustively studied, and their structure and mechanism of action have been determined.9 They are zinc-dependent proteins characterized by a large molecular mass spectrum (20 to 110 kDa), which are classified into three groups, PI to PIII, according to the number of structural domains present in their structures.10 The development of effective SVMP inhibitors may be beneficial for the treatment of victims of snakebites, resulting in a significant overall reduction in local tissue damage after the envenoming. Furthermore, the development of inhibitors to specific toxins present in snake venoms could be the first step toward an alternative to principal treatment for ophidian envenoming, the animal-derived antivenom sera. Besides the low effectiveness on the snakebite area, there are many other problems with the antivenom sera, ranging from the inadequate distribution11 to hypersensitivity reactions caused by them.12,13 New strategies are necessary for the development of molecules with inhibitory activity against toxins from snakes

nakebites, a largely ignored public health problem affecting especially poor populations in rural areas of tropical countries, are considered by WHO a neglected disease, with frequently devastating consequences to these populations.1 In Brazil, snakes from the bothropic and bothropoic genus (Viperidae) are the major causative agents of snakebites, accounting for about 80% of them.2 The venoms of snakes are composed by many toxins, which are used for immobilizing, killing, and digesting their prey.3 These toxins are responsible for various toxic effects including intense pain, bleeding disorders, and local tissue damage caused by hemorrhage and necrosis that can ultimately lead to permanent disability and may even result in limb amputation.4 To date, snake antivenom immunoglobulins (antivenoms) are the only specific treatment for envenoming by snakebites. However, although antivenoms have proved highly effective in the neutralization of systemic effects, they are only partially effective in reverting the local pathological alterations induced by snake venoms. Consequently, it has been suggested that specific inhibitors of enzymes present in the venoms may represent an improvement for the treatment of snakebites.5 In fact, the peptidomimetic hydroxamate batimastat, a known metalloproteinase inhibitor, has been proven to be effective in preventing local tissue damage and some systemic effects induced by the Bothrops asper venom.6,7 Snake venom metalloproteinases (SVMPs) (EC 3.4.24) are proteolytic enzymes acting in different hemostatic mechanisms.8 The © 2017 American Chemical Society

Received: May 1, 2017 Accepted: September 26, 2017 Published: September 26, 2017 1136

DOI: 10.1021/acsmedchemlett.7b00186 ACS Med. Chem. Lett. 2017, 8, 1136−1141

ACS Medicinal Chemistry Letters

Letter

BaP-I crystallographic structure complexed with a pipetidomimetic inhibitor (PDB code 2W12)21 was found as an adequate template. The BaP-I resolution is 1.46 Å, and it has 79% identity and 89% similarity with BpMP-I. The Swiss-MODEL server in alignment mode22 was used for generation of the BpMP-I 3D model. The root-mean-square-deviation (rmsd) between both 3D structures after superposition of their alpha carbons was 0.24 Å. Only 2% of the amino acid residues of the BpMP-I model were not in the allowed regions of the Ramachandran plot, and none of them take part in the catalytic site. Ligands structures were constructed and energy minimized with the semiempirical molecular orbital method PM623 available in Spartan ′14 (Wave function, Inc.). For ligand docking into the active site of the BpMP-I model, GOLD 5.2 (CCDC) was used. The fitness score function GoldScore24 was selected after a redocking study of the peptidomimetic inhibitor into the BaP-I crystallographic structure. Table 1 summarizes the docking results and the inhibition data.

with little or no side effects, as well for the elucidation of their mechanism of action.14 Villalta-Romero and colleagues were able to identify new SVMP inhibitors by means of compound screening and rational peptide design.15 Therefore, we explored in the present work the use of molecular modeling strategies for the rational design, followed by synthesis and evaluation of new substances for the inhibition of metalloproteinase toxins present in Bothrops pauloensis venom, a species responsible for snakebites in Brazil’s Southeast and Midwest regions. Giving the importance of the ZnII ion for the catalytic mechanism of metalloproteinases, our initial proposal was the synthesis of a series of low molecular weight molecules containing semicarbazone and thiosemicarbazone groups, which could be useful as Zn-chelating candidates, as previously described by others.16,17 Condensation reactions of thiosemicarbazide or semicarbazide with the corresponding substituted aromatic aldehydes (1), dissolved in ethanol containing a catalytic amount of hydrochloric acid, were used for the synthesis of the initial series, according to Scheme 1. After ca. 10 h of reaction

Table 1. Docking and Activity Data

Scheme 1. Synthesis of Compounds 2a−e and 3a−g

time followed by precipitation in an ice/water mixture, the solid formed was filtered off and dried furnishing the desired products (2a−e and 3a−g). 1H NMR spectroscopy unequivocally confirmed the product structures and purity of compounds. For the enzyme inhibition assay, BpMP-I, a metalloproteinase from B. pauloensis,18 was dissolved in 0.9% NaCl for protein quantification.19 Azocaseinolytic activity was determined according to Gomes et al. with some modifications.20 Briefly, the azocasein (1 mg/mL) was dissolved in Tris−HCl 0.05 M buffer (pH 7.8) containing 0.15 M NaCl. This solution was added to polystyrene 96-well plates (NUNC MaxiSorp) (160 μL/well). The reaction was initiated by the addition of 45 μL of enzyme solution containing 7.5 μg of BpMP-I. After incubation at 37 °C for 30 min, the reaction was interrupted by the addition of 45 μL of 20% (w/v) trichloroacetic acid. The plate was incubated at room temperature for 30 min and then centrifuged at 3000g for 20 min. The absorbance of the supernatant at 405 nm was determined. One unit (U) of azocaseinolytic activity was defined as an increase of 0.01 absorbance units at 405 nm under standard assay conditions. The inhibition assay was conducted by previously incubating the synthetic compounds with the toxin at 37 °C for 30 min. Immediately after, the azocaseinolytic activity determination was carried out. The IC50 was defined as the concentration capable of inhibiting 50% of azocaseinolytic activity promoted by BpMP-I. The BpMP-I sequence was used to search for structures in the Protein Data Bank to be used as templates for the construction of a 3D homology model for the enzyme. The P-I

a

compd

X

Ar

fitness scorea

IC50 (μM)

2b 2c 2d 2a 2e 3a 3b 3c 3d 3e 3f 3g

S S S S S O O O O O O O

4-OMe-Ph 4-Br-Ph 4-Cl-Ph Ph 3-OH,4-OMe-Ph 4-OMe-Ph 4-COOH-Ph 4-Cl-Ph 2,4-OH-Ph 4-OH-Ph 3-OH,4-OMe-Ph 3-OH -Ph

66.0 62.3 60.8 56.3 53.5 52.2 51.8 50.0 49.0 47.8 47.2 47.0

3.01 × 103 3.14 × 103 4.25 × 103 4.22 × 103 3.41 × 103 inactive inactive inactive inactive inactive inactive inactive

GoldScore fitness function.

According to the scores related with the best docking poses (Table 1), the thiosemicarbazones (2a−2e) are predicted as better ligands than the semicarbazones (3a−3g). It is wellknown that the performance of fitness score data as a predictive tool for activities is generally limited. In this case, however, the enzymatic inhibition assay results show that the docking predictions worked quite well, as the semicarbazones presented no activity, and the best-ranked thiosemicarbazones, 2b and 2c, were, in fact, the two most active inhibitors. Nevertheless, the thiosemicarbazones’ activities were quite modest, so we decided to explore their docking geometries as a starting point for the design of structural modifications aiming at the improvement of the BpMP-I inhibitory potency. As the p-methoxy-substituted thiosemicarbazone (2b) was the most active compound for the in vitro BpMP-I inhibition, we selected its molecular docking poses to analyze the predicted interactions with the catalytic site. In the best-scored docking pose, it was observed that the sulfur atom was coordinated to the ZnII ion of the BpMP-I, as initially planned. The terminal amino group of the thiosemicarbazone group was H bonded to the Gly109 peptidic carbonyl group, while the phenyl ring was interacting by π-stacking with the His142 side chain and the methoxy group was inserted into a hydrophobic pocket next to the catalytic site. It was also observed that below the iminic carbon there was a gorge containing some polar groups, including some backbone peptide bonds, including 1137

DOI: 10.1021/acsmedchemlett.7b00186 ACS Med. Chem. Lett. 2017, 8, 1136−1141

ACS Medicinal Chemistry Letters

Letter

those from Gln106, Ile108, Gly109, and Ile169 (Figure 1). We decided to explore these potential interaction points by

Figure 1. Best scored docking pose for compound 2b inside BpMP-I catalytic site. Figure generated with PyMol software.

introducing a carboxylate group connected to the iminic carbon by variable sized methylene spacers. It would be also expected that the presence of a carboxylate group would have as a second beneficial effect the improvement of the compounds’ water solubility. Molecules containing the carboxylate connected to the iminic carbon by two and three methylene groups were constructed with the Spartan ′14 software, and docking experiments were conducted with these molecules into the BpMP-I catalytic site, with the same methodology described before. In order to get a better evaluation of the energetic effect of the proposed structural modifications, the best scored BpMP-I poses were reoptimized with MOPAC 2012 software, employing the semiempirical molecular orbital method PM623 for the calculation of the interaction enthalpy between the new ligands and the enzyme in the “gas phase” (ΔHint), according to eq 1: ΔHint = ΔHEL + ΔHW − (ΔHEW + ΔHL)

Figure 2. Best scored docking poses (GoldScore) for compounds 5a (A) and 5b (B) inside BpMP-I catalytic site. Figure generated with PyMol software.

compounds the carboxylate group is involved in an H bond with the peptidic NH group of a residue located in the gorge, Ile169. Both compounds also make another H bond with a peptidic carbonyl group, 5a with that of Ser167 and 5b with Glu109. Based on these results, we synthesized both compounds as depicted in Scheme 2. In the first step, we obtained the χ- and δ-ketoacids intermediates by means of a Friedel−Crafts acylation of anisole, using succinic and glutaric anhydride, respectively, for 4a and 4b, and aluminum chloride. After about 3−4 h of reaction times, TLC analyses indicated the formation of the keto acid compounds, which were isolated, after precipitation in a cold solution of HCl and washing with ethylic ether, in yields of 55 and 48% for 4a and 4b, respectively. In the second step, the intermediates 4a and 4b were reacted with thiosemicarbazide in ethanol with a catalytic amount of acetic acid. The reactions were refluxed for approximately 6 h, when TLC analyses indicated the total consumption of starting materials. Compound 5a was isolated as a pure compound in 59% yield, just pouring the reaction onto ice, filtrating, and washing the precipitate formed with cold ethanol. However, precipitation of 5b led to some impurities that were eliminated after performing a purification in a chromatographic column (DCM−methanol), furnishing 5b in 74% yield. HPLC, mass

(1)

where ΔHEL, ΔHW, ΔHEW, and ΔHL are the calculated enthalpies of formation of the enzyme−ligand complex, of a water molecule, of the enzyme−water complex, and of the ligand, respectively, assuming that the ZnII ion was initially coordinated to a water molecule. Both molecular docking and semiempirical results predicted that the proposed structural modifications would result in a large improvement of the interaction profile of the compounds with BpMP-I (see Table 2). Figure 2 presents the best poses obtained for the modified compounds. As can be seen, for both Table 2. Docking, Semiempirical Calculations, and BpMP-I Inhibitory Activity Data compd

fitness scorea

ΔHint (kcal/mol)b

IC50 (μM)

2b 5a 5b

66.0 71.5 73.2

0.0 −65.2 −86.4

3.01 × 103 79.12 1.77

GoldScore fitness function. semiempirical method.

a

b

Relative interaction enthalpy; PM6 1138

DOI: 10.1021/acsmedchemlett.7b00186 ACS Med. Chem. Lett. 2017, 8, 1136−1141

ACS Medicinal Chemistry Letters

Letter

with a reduction of 100% occurring with a ratio of 1:10 Bp/ inhibitor (w/w).

Scheme 2. Two-Step Synthesis of 5a and 5b

spectrometry, and 1H and 13C NMR spectroscopy results unequivocally confirmed the structure and purity of both compounds. The final products 5a and 5b were evaluated in the BpMP-I inhibition assay with the same procedure described before. The results show that, as expected from the theoretical results, the designed structural modifications caused a large increase in the inhibitory activity of both compounds: 5a and 5b were 38 and 1700 times more active than the original compound 2b, respectively (Table 2). As a reference, when BpMP-I was incubated with the known zinc chelator and metalloproteinase inhibitor 1,10-phenanthroline, under the same experimental conditions used for the determination of activity of our compounds, an IC50 higher than those of compounds 5a and 5b was obtained (4.7 × 102 μM). The in vitro results prompted us to evaluate the in vivo effect of the compounds. Hemorrhagic activity was assessed according to Nikai and collaborators,25 with modifications. Initially, groups of Swiss male mice (n = 4, 18−22 g) were injected intradermally in the dorsal region with two minimum hemorrhagic doses (MHD) of B. pauloensis crude venom (Bp) (16.3 μg).26 After 150 min, the animals were anesthetized (ketamine 10% (0.05 mL/kg) + xylazine 2% (0.025 mL/kg)) and sacrificed. Their skins were removed and the hemorrhagic halo was measured with a digital caliper DIGMESS 100.174BL. Control mice received only saline or Bp. For hemorrhagic activity, another protocol was developed in order to simulate a treatment of envenomation (treatment condition). For this purpose, animals received solutions of Bp (16.3 μg), then 10 min later the animals were inoculated with 5b compound by the same route at a 1:10 ratio (w/w). Each experiment was expressed as the mean ± SD. The antihemorrhagic potential of the compounds was determined as percent inhibition of hemorrhagic activity. In accordance with the in vitro enzyme inhibition results, compounds 5a and 5b were able to significantly reduce the hemorrhagic halo induced by B. pauloensis venom in 43% and 100%, respectively, with a ratio of 1:20 Bp/inhibitor (w/w). Compound 5b was then evaluated to determine its dosedependent inhibition of the hemorrhagic halo induced by B. pauloensis venom. As can be seen in Figure 3, compound 5b reduced the venom effect in vivo in a dose-dependent manner

Figure 3. Dose-dependent inhibition by 5b of hemorrhagic activity induced by B. pauloensis venom (Bp). Control mice received saline solution (NaCl 0.9%; m/v) or Bp. (A,C) Inhibition percentage of hemorrhagic activity. (B,D) Representative images of hemorrhagic halos. The same sequence was used for A and B. The ratios of 1:1, 1:5, and 1:10 Bp/5b (w/w). Treatment condition for C and D: injection of Bp, and after 10 min, the animals were inoculated by the same route with the inhibitor 5b at ratio 1:10 (Bp/5b w/w). Results are reported as mean ± SD (n = 4). The difference in inhibition was significant (p < 0.05).

The inhibition was also observed when 5b was inoculated 10 min after venom inoculation (treatment condition), when 5b significantly but not completely suppressed the hemorrhagic process. Figure 3C,D showed that 5b was able to reduce the hemorrhage induced by B. pauloensis venom in approximately 31% (Figure 3C,D). The antihemorrhagic effect of 5b (treatment condition) was similar to the action of BaltMPI, an SVMP inhibitor from B. alternatus serum,27 and triacontyl pcoumarate, isolated from root bark of Bombacopsis glabra extract.28 Preincubation of BaltMPI with Batroxase, a SVMP from B. jararaca, improved the antihemorrhagic effect, when compared with the sequential administration without previous incubation. In this study, the authors related this effect probably to the formation of inactive metalloprotease−inhibitor complexes during the incubation period.27 The treatment 1139

DOI: 10.1021/acsmedchemlett.7b00186 ACS Med. Chem. Lett. 2017, 8, 1136−1141

ACS Medicinal Chemistry Letters

Letter

5b in the presence of increasing amounts of Zn(OAc)2 as a source of ZnII. 1H NMR spectra were recorded and showed unequivocally that an interaction occurred between ZnII and the ligand’s thiosemicarbazone group, in accordance with the complexation mode predicted by our theoretical study (1H NMR data presented as Supporting Information).

condition is an interesting method for evaluating the neutralizing ability of an inhibitor during experimental envenomation. Our results showed that the protective effect of compound 5b against the hemorrhagic effect of B. pauloensis venom can open up new perspectives for its use as a complement for serum therapy. Even though the design of compound 5b was based on the catalytic site of BpMP-I, a nonhemorrhagic (since it has only the metal-binding domain) P-I type metalloprotease, it was able to bind to the catalytic site and effectively inhibit hemorrhagic P-II and P-III types of SVMPs present in B. pauloensis venom.10 Inhibition of other SVMPs present in the venom can be explained by the fact that this class of toxins has a highly conserved catalytic site.29 This probably allows 5b to establish interactions with the catalytic site of other SVMPs similar to those occurring in BpMP-I, from which it was designed. The good performance of GoldScore fitness function in the correlation between scores and activity data (Tables 1 and 2) prompted us to investigate by molecular docking, with the same procedure described before, if compounds 5a and 5b were also able to interact favorably with SVMP structures from other snake species. For this study, we selected crystal structures from seven SVMPs deposited in PDB, using as criteria for selection the crystallographic resolution, sequence identity, and degree of relatedness with the B. pauloensis enzyme: 4Q1L (B. leucurus), 2W12 (B. asper), 3GBO (B. moojeni), 1HTD (Crotalus atrox), 3AIG (C. adamanteus), 1BUD (Agkistrodon acutus), and 4J4M (Protobothrops mucrosquamatus). The docking results are indicative that the inhibitors would perform quite well when binding in the catalytic site of the selected SVMPs, presenting fitness scores very similar to those obtained in the BpMP-I catalytic site (Table 3).



S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsmedchemlett.7b00186. General synthesis procedure and NMR spectra, HPLC, and exact mass analysis of compounds 5a and 5b; effect of zinc acetate additions on 1H NMR of inhibitor 5b (PDF)



*E-mail: [email protected]. Phone: 55 21 2682 1872. ORCID

Carlos Mauricio R. Sant’Anna: 0000-0003-1989-5038 Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The authors would like to thank Brazilian agencies CNPq, CAPES, and FAPERJ, and INCT-INOFAR for fellowships and financial support.

a



ABBREVIATIONS BaP-I, Bothrops asper type I metalloproteinase; Bp, B. pauloensis crude venom; BpMP-I, B. pauloensis type I metalloproteinase; MHD, minimum hemorrhagic doses; SVMP, snake venom metalloproteinase; TLC, thin layer chromatography

fitness scorea

code

2a

5a

5b

BpMP-I 2W12 4Q1L 3GBO 1HTD 3AIG 1BUD 4J4M

66.01 61.70 63.94 63.16 61.89 62.00 54.63 58.87

71.50 66.28 75.49 70.57 71.84 72.9 62.40 65.19

73.17 69.68 79.76 73.88 72.28 72.6 64.87 66.83

AUTHOR INFORMATION

Corresponding Author

Table 3. Docking Fitness Scores for Compounds 2a, 5a, and 5b in SVMPs SVMP

ASSOCIATED CONTENT



REFERENCES

(1) Warrell, D. A. Snake bite. Lancet 2010, 375, 77−88. (2) Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância Epidemiológica. Guia de Vigilância Epidemiológica, 7th ed. Brası ́lia, 2009. (3) Sajevic, T.; Leonardi, A.; Krizaj, I. Haemostatically active proteins in snake venoms. Toxicon 2011, 57, 627−645. (4) Gutierrez, J. M.; Williams, D.; Fan, H. W.; Warrell, D. A. Snakebite envenoming from a global perspective: Towards an integrated approach. Toxicon 2010, 56, 1223−1235. (5) Gutiérrez, J. M.; Theakston, R. D. G.; Warrell, D. A. Confronting the Neglected Problem of Snake Bite Envenoming: The Need for a Global Partnership. PLoS Med. 2006, 3, e150. (6) Rucavado, A.; Escalante, T.; Franceschi, A.; Chaves, F.; León, G.; et al. Inhibition of local hemorrhage and dermonecrosis induced by Bothrops asper snake venom: Effectiveness of early in situ administration of the peptidomimetic metalloproteinase inhibitor batimastat and the chelating agent CaNa2EDTA. Am. J. Trop. Med. Hyg. 2000, 63, 313−319. (7) Rucavado, A.; Escalante, T.; Gutiérrez, J. M. Effect of the metalloproteinase inhibitor batimastat in the systemic toxicity induced by Bothrops asper snake venom: Understanding the role of metalloproteinases in envenomation. Toxicon 2004, 43, 417−424.

GoldScore fitness function.

The values of the scores followed the same trend observed for the values obtained for BPMP-I, with modified compounds 5a and 5b obtaining higher scores when compared with the original inhibitor 2a and compound 5b presenting the highest scores, except for C. adamanteus toxin, where the score was higher for compound 5a. These results are suggestive that compounds 5a and 5b are also candidates to be effective inhibitors of the hemorrhagic effects caused by SVMPs present in the venom of other snake species. Finally, as compounds 5a and 5b contain two groups that could potentially coordinate with the ZnII ion, i.e., the thiosemicarbazone and the carboxylate groups, a question remained about which one was in fact responsible for the interaction. We implemented a model study with solutions of 1140

DOI: 10.1021/acsmedchemlett.7b00186 ACS Med. Chem. Lett. 2017, 8, 1136−1141

ACS Medicinal Chemistry Letters

Letter

the venom of Crotalus atrox (western diamondback rattlesnake). Arch. Biochem. Biophys. 1984, 231, 309−319. (26) Rodrigues, V. M.; Soares, A. M.; Andrião-Escarso, S. H.; Franceschi, A. M.; Rucavado, A.; Gutiérrez, J. M.; Giglio, J. R. Pathological alterations induced by neuwiedase, a metalloproteinase isolated from Bothrops neuwiedi snake venom. Biochimie 2001, 83, 471−479. (27) Palacio, T. Z.; Santos-Filho, N. A.; Rosa, J. C.; Ferreira, R. S.; Barraviera, B.; Sampaio, S. V. Isolation and characterization of a novel metalloprotease inhibitor from Bothrops alternatus snake serum. Int. J. Biol. Macromol. 2017, 98, 436−446. (28) Mendes, M. M.; Vieira, S. A. P. B.; Gomes, M. S. R.; Paula, V. F.; Alcântara, T. M.; Homsi-Brandeburgo, M. I.; dos Santos, J. I.; Magro, A. J.; Fontes, M. R. M.; Rodrigues, V. M. Triacontyl p-coumarate: An inhibitor of snake venom metalloproteinases. Phytochemistry 2013, 86, 72−82. (29) Stöcker, W.; Bode, W. Structural features of a superfamily of zinc-endopeptidases: the metzincins. Curr. Opin. Struct. Biol. 1995, 5, 383−390.

(8) Tanjoni, I.; Evangelista, K.; Della-Casa, M. S.; Butera, D.; Magalhães, G. S.; Baldo, C.; Clissa, P. B.; Fernandes, I.; Eble, J.; Moura-da-Silva, A. M. Different regions of the class P-III snake venom metalloproteinase jararhajin are involved in binding to α2β1 integrin and collagen. Toxicon 2010, 55, 1093−1099. (9) Markland, F. S., Jr.; Swenson, S. Snake venom metalloproteinases. Toxicon 2013, 62, 3−18. (10) Fox, J. W.; Serrano, S. M. T. Insights into and speculations about snake venom metalloproteinase (SVMP) synthesis, folding and disulfide bond formation and their contribution to venom complexity. FEBS J. 2008, 275, 3016−3030. (11) Gutiérrez, J. M.; Rucavado, A.; Escalante, T.; Lomonte, B.; Angulo, Y.; Fox, J. W. Tissue pathology induced by snake venoms: how to understand a complex pattern of alterations from a systems biology perspective. Toxicon 2010, 55, 166−170. (12) Bucaretchi, F.; Herrera, S. R. F.; Hyslop, S.; Bacarat, E. C. E.; Vieira, R. J. Snakebites by Crotalus durissus ssp in children en Campinas, São Paulo, Brazil. Rev. Inst. Med. Trop. Sao Paulo 2002, 44, 133−138. (13) Thiansookon, A.; Rojnuckarin, P. Low incidence of early reactions to horse-derived F(ab)2 antivenom for snakebites in Thailand. Acta Trop. 2008, 105, 203−205. (14) Perales, J.; Neves-Ferreira, A. G. C.; Valente, R. H.; Domont, G. B. Natural inhibitors of snake venom hemorrhagic metalloproteinases. Toxicon 2005, 45, 1013−1020. (15) Villalta-Romero, F.; Gortat, A.; Herrera, A. E.; Arguedas, R.; Quesada, J.; de Melo, R. L.; Calvete, J. J.; Montero, M.; Murillo, R.; Rucavado, A.; Gutiérrez, J. M.; Pérez-Payá, E. Identification of New Snake Venom Metalloproteinase Inhibitors Using Compound Screening and Rational Peptide Design. ACS Med. Chem. Lett. 2012, 3, 540− 543. (16) Ali, M. A.; Mirza, A. H.; Chartres, J. D.; Bernhardt, P. V. Synthesis, characterization and X-ray crystal structures of sevencoordinate pentagonal-bipyramidal zinc(II), cadmium(II) and tin(IV) complexes of a pentadentate N3S2 thiosemicarbazone. Polyhedron 2011, 30, 299−306. (17) Gambino, D. The Wide Pharmacological Versatility of Semicarbazones, Thiosemicarbazones and Their Metal Complexes. Mini-Rev. Med. Chem. 2004, 4, 31−39. (18) Naves-de-Souza, D. L.; Gomes, M. S. R.; Ferreira, F. B.; Rodrigues, R. S.; Achê, D. C.; Richardson, M.; Borges, M. H.; Rodrigues, V. M. Biochemical and enzymatic characterization of BpMP-I, a fibrinogenolytic metalloproteinase isolated from Bothropoides pauloensis snake venom. Comp. Biochem. Physiol., Part B: Biochem. Mol. Biol. 2012, 161, 102−109. (19) Bradford, M. M. Rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem. 1976, 72, 248−254. (20) Gomes, M. S. R.; de Queiroz, M. R.; Mamede, C. C. N.; Mendes, M. M.; Hamaguchi, A.; Homsi-Brandeburgo, M. I.; Sousa, M. V.; Aquino, E. N.; Castro, M. S.; de Oliveira, F.; Rodrigues, V. M. Purification and functional characterization of a new metalloproteinase (BleucMP) from Bothrops leucurus snake venom. Comp. Biochem. Physiol., Part C: Toxicol. Pharmacol. 2011, 153, 290−300. (21) Lingott, T.; Schleberger, C.; Gutiérrez, J. M.; Merfort, I. Highresolution crystal structure of the snake venom metalloproteinase BaP1 complexed with a peptidomimetic: insight into inhibitor binding. Biochemistry 2009, 48, 6166−6174. (22) Arnold, K.; Bordoli, L.; Kopp, J.; Schwede, T. The SWISSMODEL Workspace: A web-based environment for protein structure homology modelling. Bioinformatics 2006, 22, 195−201. (23) Stewart, J. J. P. Optimization of parameters for semiempirical methods VI: more modifications to the NDDO approximations and re-optimization of parameters. J. Mol. Model. 2013, 19, 1−32. (24) Liebeschuetz, J. W.; Cole, J. C.; Korb, O. Pose prediction and virtual screening performance of GOLD scoring functions in a standardized test. J. Comput.-Aided Mol. Des. 2012, 26, 737−748. (25) Nikai, T.; Mori, N.; Kishida, M.; Sugihara, H.; Tu, A. T. Isolation and biochemical characterization of hemorrhagic toxin from 1141

DOI: 10.1021/acsmedchemlett.7b00186 ACS Med. Chem. Lett. 2017, 8, 1136−1141