2.1. Isolation and structural elucidation of bufadienolides
The crude extract from oocytes of
R. alata was prepared by maceration with a mixture of methanol and chloroform and then submitted to C-18 SPE purification. The fraction eluted at 80 % CH
3OH showed anti-trypanosomal activity and was further analyzed by semi-preparative RP-HPLC (
Figure S1). Compounds
1-
7 were isolated and their structures, depicted in
Figure 1, were determined based on extensive spectroscopic analysis and comparison of their data with those reported in the literature. While compounds 1, 2 and 4-7 corresponded to known compounds, the bufadienolide
3 is reported as a new natural product.
Compound
3 was isolated as a white solid. The molecular formula was determined as C
27H
34O
8 based on its HR-MALDI-MS data which showed the molecular ion at
m/z 487.2283 [M + H]
+ (calculated for C
27H
35O
8, 487.2326). Eleven degrees of unsaturation were inferred from the molecular formula. The
1H-NMR and
13C-NMR spectra of compound
3 showed the feature signals for a 2H-pyran-2-one moiety [δ
H 6.24 (1H, dd,
J = 0.8, 9.7 Hz), 7.36 (1H, m), 8.02 (1H, m); δ
C 114.2, 118.2, 150.9, 153.6, 164.0], a tertiary methyl [δ
H 0.81 (3H, s); δ
C 17.3], a secondary methyl [δ
H 0.94 (3H, t,
J = 7.8 Hz); δ
C 9.3], three oxygenated methines [(δ
H 3.75, m; δ
C 60.7), (δ
H 4.15, m; δ
C 68.0), (δ
H 5.47, m; δ
C 76.3)], two oxygenated quaternary carbons [δ
C 75.8, 73.1] and two carbonyl carbons characteristic of aldehyde and ester [δ
C 210.1, 174.9]. Comparison of the NMR signals (
Table 1) of compound
3 with those of the reported bufadienolide dyscinobufotalin [
10]. revealed the data of the steroidal nucleus of the two compounds were very similar, indicating that compound
3 had the same bufadienolide core with two carbinols at C-3 and C-5, an epoxy group at C-14 and C-15, and the esterification of the hydroxyl group at C-16 with propionic acid. The main difference consisted in the presence of an aldehyde group at C-19 in compound
3 instead of the angular methyl group at C-19 of dyscinobufotalin [
10]. HMBC correlations between aldehyde proton at δ
H 10.00 and the quaternary carbon at C-10 allowed for the assignment of the aldehyde group at C-19 (
Figure 1). On the basis of spectral data, the compound
3 was identified as 19-formyl-dyscinobufotalin.
Furthermore, six steroids of the bufadienolide family were identified as 16β-hydroxyl-hellebrigenin (
1) [
11], desacetyl-bufotalin (
2) [
5], bufotalin (
4) [
5], cinobufotalin (
5) [
12], dyscinobufotalin (
6) [
10], and bufalin (
7) [
5].
2.3. Molecular docking of bufadienolides 1, 3 and 7
A docking protocol was validated and its reproducibility confirmed by the root mean square deviation (RMSD) values obtained. The RMSD for each cruzipain-co-crystallized inhibitor complex were as follows: 2.3586 for 1EWO-VSC, 1.9060 for 2AIM-ZRA, 1.6817 for 1U9Q-186 and 2.9319 for 4XUI-2VC. In addition, relative docking scores for each complex resulted in -7.0484, -6.9388, -8.2353 and -8.1351, respectively. After validation of the protocol accuracy, interactions of anti-trypanosomal bufadienolides
1,
3 and
7 with the active site of proteins were studied via docking calculations.
Table 3 shows the calculated parameters for the best pose obtained for each bufadienolide tested. The energy variation for the predictions obtained by simulation ranged from -2.2 to +1.3 kcal/mol. The isolated bufadienolides
1,
3 and
7 displayed interactions with cruzipain. Bufadienolide
1 showed H-bonds with ASP-158 and GLN-19 at a distance of 2.10 Å and 2.32 Å, respectively. On the other hand, bufadienolide
3 presents three H-bonds with amino acid residues GLN-19, HIS-159 and TRP-177. Bufadienolide
7 formed one H-bond with ASP-158 at a distance of 1.95 Å (
Figure 2).
Currently, there is no effective chemotherapeutic treatment available for American trypanosomiasis [
4]. Poisonous terrestrial vertebrates are rarely considered for natural product research programs; however, it has been reported that bioactive materials derived from animals, have resulted in effective medicines, highlighting the potential of vertebrate animals’ toxins in pharmacology. For instance, teprotide a nonapeptide isolated from the venom of the snake
Bothrops jararaca, has served as template to produce captopril and other angiotensin-converted enzyme inhibitors that currently are routinely medicated for hypertension patients [
13]. The alkaloid epibatidine obtained from the skin of the dendrobatid frog
Epipedobates tricolor presents analgesic properties 200 times more potent than morphine. Epibatidine has served as template for producing anti-nociceptive derivatives such as tebanicline (ABT-594) that reached clinical phase II [
14,
15].
In this study, through the analysis of the spectroscopic data, we elucidated the chemical structure of the new bufadienolide
3 from the oocytes of
R. alata. Moreover, five bufadienolides (
2,
4-7) were previously reported in amphibians and one in plants (
1). Bufadienolide
2 was isolated from the poison of the toad
Duttaphrynus melanostictus as a minor constituent [
5]. Compounds
4 and
5 have been obtained from skin and parotoid gland secretions of bufonids, but not in toads of the genus
Rhinella [
16]. The bufadienolide
6 was isolated from Chan Su, a traditional chinese remedy prepared with the poison of Asian toads [
10]. The bufadienolide
7 has been identified in bufonid extracts including oocytes [
17].
Due to the marked inhibitory activity against the pump (Na
+-K
+)ATPase, the antiparasitic potential of bufadienolides against human pathogens has currently received some attention [
18]. In this study, bufadienolides
1-
7 were evaluated in vitro for inhibition of the growth of the trypomastigote form of
T. cruzi. According to our data, substitution of methyl at C-19 by formyl group increases the anti-trypanosomal activity as seen in bufadienolide
3 in comparison with bufadienolide
6. Similarly, in a previous study hellebrigenin showed trypanocidal activity at 91.7 µg/mL while telocinobufagin, which contains a C-19 methyl group instead, was inactive [
19]. Likewise, comparing the activities of bufadienolides
2 and
7, it can be inferred that the lack of a hydroxyl group at C-16 is required for the anti-trypanosomal activity, however if this hydroxyl group is esterified, a moderate bioactivity is observed (compounds
4-
6,
Table 2).
Also, previous evaluations on cancer cell lines with bufadienolides revealed a decrease in growth inhibition by introduction of a hydroxyl group at C-16 position [
20]. Bufadienolide
7 presents higher anti-trypanosomal activity than its 19-OH analogue (IC
50 = 19.4 µM) [
21]. On the basis of these results, we can suggest that the presence of a formyl group at C-19 and the esterification of the hydroxyl group at C-16 are essential elements required for the effective biomolecular interactions with components present in the trypomastigote form of
T. cruzi, which led to the anti-protozoal activity observed.
In mammals, bufadienolides block the transport activity of the (Na
+-K
+)ATPase pump by binding to the α-subunit. This inhibition provokes that intracellular sodium concentration increases until cells get depolarized [
22]. In
T. cruzi, there is no report on the identification of this protein; maintenance of the ionic steady-state is carried out by K
+ channels, Na
+ efflux pumps and (H
+)ATPases [
23]. Kyoichi, I.
et al. [
24] cloned and characterized a (Na
+)ATPase from a gene encoding of
T. cruzi. After biochemical evaluations it was evidenced that the enzymatic activity of this pump was not inhibited by the steroid (cardenolide) ouabain. In the search for anti-trypanosomal drugs, the most considered targets for Chagas disease are cruzipain, sterol-14α-demethylase, trans-syalidase, trypanothione reductase and tubulin [
4,
24]. Currently, the mechanism by which bufadienolides cause inhibition of
T. cruzi is unknown, despite of previous findings suggest that bufadienolides inhibit the growth of
T. cruzi parasites by a biochemical pathway other than the blockage of the (Na
+-K
+)ATPase pump. In this study, in silico analysis by molecular docking was carried out, using the
T. cruzi protease cruzipain. Redocking with previously reported co-crystallized inhibitors was employed for protocol validation and active site prediction (
Figure S9). Molecular docking with the anti-trypanosomal bufadienolides
1,
3 and
7 revealed strong interactions with cruzipains. These interactions correspond to H-bonds between amino acid residues GLN-19, ASP-158, HIS-159 and TRP-177 with hydroxyl groups at C-5 and C-14 position; and formyl at C-19 position of bufadienolides. The most potent anti-trypanosomal bufadienolide (
3), presented three interactions with cruzipain while bufadienolides
1 and
7 showed two and one H-bonds, respectively. Bufadienolide
7 interacts with cruzipain by the hydroxyl at C-14, a functional group essential for bioactivity [
25]. No interaction with polar groups at C-16 such as hydroxyl and propionyloxy were observed for bufadienolides
1 and
3. Cruzipain contains, at the active site, a catalytic triad formed by GLY-23, CYS-25 and GLY-65 [
26]. Furthermore, some amino acid residues have been revealed to be involved in enzymatic catalysis of cruzipain such as residues GLN-19, GLY-66, ASP-158, HIS-159 and TRP-177 that keep peptidyl inhibitors anchored to the active site [
27]. A molecular docking analysis with 173 compounds reported to be inhibitors of cruzipain showed that in addition to the catalytic triad, the amino acid residues GLN-19 and ASP-158 are involved in ligand interactions with a frequency of 47 % and 44 %, respectively [
26].
Selectivity is a well-accepted parameter of the safety and pharmacological activity of a compound [
28]. Bufadienolides from oocytes of
R. alata were evaluated against kidney monkey Vero cells. The bufadienolides
1 (EC
50 = 22.80 µM, S.I. > 5) and
7 (EC
50 = 3.93 µM, S.I. = 2.7) presented significant selectivity (
Table 2). Although bufadienolide
7 shows a promising S.I., less cytotoxicity is required. On the contrary, the low cytotoxicity showed by the bufadienolide
1 highlights its potential for further studies. Recently, 19-hydroxy-bufalin was evaluated against Vero cells and its low cytotoxicity, in comparison with bufadienolide
7, suggest that hydroxyl at position 19 decrease cytotoxic effects against mammals [
21]. Despite of the potent anti-trypanosomal effect shown by compounds
3,
5 and
6, these bufadienolides were toxic to Vero cells, resulting in a reduced selectivity for the assayed parasite. In a previous study, hellebrigenin and telocinobufagin, two bufadienolides with antiprotozoal activity were tested for cytotoxicity on mouse macrophages and erythrocytes, showing no toxic effects [
19]. However, it has been found that rodent cells are over 1000-fold more resistant than human cells to cytotoxic effects of cardiotonic sterols [
29]. Differential cytotoxicity and water solubility are great attributes for potential anti-trypanosomals. In this context, microbial transformations carried out with bufadienolides
4,
5 and
7 have produced nontoxic derivatives to human cells [
30,
31,
32]. Nanoparticles containing bufadienolide
7 were injected in mice xenografted with human colon cancer cells, and their therapeutic potential was tested. The particles were able to inhibit the growth of tumor cells with efficient tumor targeting and cellular uptake [
33].
Some bufadienolides have been found to be active against amphibian pathogens. Arenobufagin, γ-bufotalin and telocinobufagin isolated from the skin secretion of the boreal toad
Anaxyrus boreas inhibit the harmful fungus
Batrachochytrium dendrobatidis in vitro [
34]. Amphibians are also affected by parasitic diseases caused by protozoa [
35]. In fact, different families of anurans, including Bufonidae, are hosts of trypanosomes and are frequently parasitized by more than one species [
36]. Blood samples from
R. alata have not been examined for the presence of trypanosomes; however, congeneric species, as well as syntopic species such as the frog
Engystomops pustulosus, are known to be parasitized by trypanosomes [
37,
38]. Trypanosomiasis in amphibians is characterized by anemia, food refusal, listlessness and localized hemorrhages with swollen lymph glands. Species of leeches from the family Hirudinea are considered to be the aquatic vectors of amphibian trypanosomes, while mosquitos, sand flies and midges, their main terrestrial vectors [
37,
38,
39]. The presence of anti-trypanosomal bufadienolides in oocytes of toads is intriguing. Bufonids biosynthesize bufadienolides, as shown by studies using cholesterol marked isotopically [
40]. These compounds have been isolated from parotoid glands and skin of toads, as well as from their plasma, bile, ovaries and oocytes [
17,
41,
42,
43]. Because they inhibit (Na
+-K
+)ATPase pump, bufadienolides are considered as possible regulators of ionic equilibrium in anurans, and the high levels of bufadienolides present in some bufonid species presumably serve as defense against predation [
44]. However, a possible role as endogenous antiprotozoal substances should not be discarded.