Skip to main content

Novel coumarins active against Trypanosoma cruzi and toxicity assessment using the animal model Caenorhabditis elegans

Abstract

Background

Chagas disease (CD) is a tropical parasitic disease. Although the number of people infected is very high, the only drugs available to treat CD, nifurtimox (Nfx) and benznidazole, are highly toxic, particularly in the chronic stage of the disease. Coumarins are a large class of compounds that display a wide range of interesting biological properties, such as antiparasitic. Hence, the aim of this work is to find a good antitrypanosomal drug with less toxicity. The use of simple organism models has become increasingly attractive for planning and simplifying efficient drug discovery. Within these models, Caenorhabditis elegans has emerged as a convenient and versatile tool with significant advantages for the toxicological potential identification for new compounds.

Methods

Trypanocidal activity:

Forty-two 4-methylamino-coumarins were assayed against the epimastigote form of Trypanosoma cruzi (Tulahuen 2 strain) by inhibitory concentration 50% (IC50).

Toxicity assays:

Lethal dose 50% (LD50) and Body Area were determined by Caenorhabditis elegans N2 strain (wild type) after acute exposure.

Structure-activity relationship:

A classificatory model was built using 3D descriptors.

Results

Two of these coumarins demonstrated near equipotency to Nifurtimox (IC50 = 5.0 ± 1 μM), with values of: 11 h (LaSOM 266), (IC50 = 6.4 ± 1 μM) and 11 g (LaSOM 231), (IC50 = 8.2 ± 2.3 μM). In C. elegans it was possible to observe that Nfx showed greater toxicity in both the LD50 assay and the evaluation of the development of worms. It is possible to observe that the efficacy between Nfx and the synthesized compounds (11 h and 11 g) are similar. On the other hand, the toxicity of Nfx is approximately three times higher than that of the compounds. Results from the QSAR-3D study indicate that the volume and hydrophobicity of the substituents have a significant impact on the trypanocidal activities for derivatives that cause more than 50% of inhibition. These results show that the C. elegans model is efficient for screening potentially toxic compounds.

Conclusion

Two coumarins (11 h and 11 g) showed activity against T. cruzi epimastigote similar to Nifurtimox, however with lower toxicity in both LD50 and development of C. elegans assays. These two compounds may be a feasible starting point for the development of new trypanocidal drugs.

Background

Chagas disease (CD), or American trypanosomiasis, is a tropical parasitic disease that affects approximately 6-7 million people worldwide, predominantly in America [1]. CD is caused by the flagellated protozoan Trypanosoma cruzi, often transmitted to humans and to other mammals by blood-sucking insects, triatomines [2]. Although the number of people infected is very high and most cases occur in rural areas in Latin American countries [3], the only drugs available to treat CD are Nifurtimox (Nfx) and Benznidazole, which are highly toxic and often ineffective, particularly in the chronic stage of the disease [4]. At the same time, no drug has emerged as an effective candidate for clinical trials in the last 30 years [5]. Therefore, the identification of lead compounds that can be easily taken forward into lead optimization to yield drug candidates to be tested in human clinical trials is an essential move for the development of safer and more effective drugs for the treatment of CD.

Coumarins are a large class of compounds that display a wide range of interesting biological properties such as anticoagulant [6], antimicrobial activities [7], antioxidant [8], and anticancer [9]. They are considered good examples of “privileged structures”, usually rigid, polycyclic heteroatomic systems that are capable of binding to multiple pharmacological targets, thus providing a viable starting point in the search for new drugs [10, 11]. Considering the use of “privileged structures” as a feasible strategy to design new antitrypanosomal drug candidates with favourable pharmacokinetic/toxicity profiles [12], the 2H-chromen-2-one nucleus of coumarins was elected as the main structural feature for the development of new lead compounds with potential trypanocidal activity. This planar ring system is composed of one aromatic ring, capable of establishing hydrophobic, π-π, CH-π and cation-π interactions, and one lactone ring, which contains two oxygen atoms that may interact via hydrogen bonding with a series of amino acid residues, such as serine, threonine, cysteine, asparagine, glutamine and tyrosine. The C-C double bond in the 2-pyrone ring is also essential for conferring planarity to the 2H-chromen-2-one core, allowing for charge delocalization between the carbonyl group and the aromatic ring [13].

Although several studies report the therapeutic potential of coumarins in parasitic diseases such as leishmaniasis, malaria and amebiasis [14], there are few papers describing bioactive coumarins in Trypanosoma cruzi [15,16,17,18,19]. Chalepin, a coumarin isolated from Ruta Angustifolia L. Pers, was shown to be active against T. cruzi with an IC50 = 64 μM. The authors propose as an action mechanism, the glycolytic enzyme glyceraldehyde 3-phosphate dehydrogenase inhibition [15]. The synthetic coumarin described by Oliaro-Bosso et al., proposes the oxidosqualene cyclases as the anti-trypanosome target, with an IC50 = 0.36 μM [19]. Another natural coumarin, Soulamarin, showed an IC50 = 210.1 μM. The authors proposed that the mitochondrial dysfunction and the modification of the plasma membrane permeability as the mechanisms of action [18].

This article describes a complete workflow with new compound synthesis using an in vitro essay against T. cruzi and an in vivo toxicological test using the C. elegans model (Fig. 1). Highlight is made of the synthesis of a library of forty-two 4-methylamino-coumarins containing aromatic and aliphatic rings attached to the amino group, which had their trypanocidal activity assayed in vitro on T. cruzi epimastigotes. The in vivo toxicity profile was assessed for some of the most promising coumarins using the nematode C. elegans as a model for acute toxicity. Finally, a three-dimensional structure-activity relationship (SAR) study was carried out using Molecular Interaction Fields (MIF) by the GRid Independent Descriptors (GRIND) approach. GRIND program was used to determine groups that increase trypanocidal activity.

Fig. 1
figure1

Workflow of new compounds development with anti-trypanosome activity

Methods

Chemistry

All chemicals were purchased as reagent grade and used without further purification. Solvents were distilled and/or dried in accordance with standard methods [20]. Column chromatography was performed on silica gel 60 (0.040-0.063 mm) using a hexane/ethyl acetate gradient. Melting points were determined on a Fisatom 431 apparatus and are uncorrected. Infrared spectra were recorded on Perkin Elmer model FT-Spectrum BRXII/Miracle ATR spectrometer. MS spectra were recorded on Q-TOF micro Waters high-resolution mass spectrometer, operating in electrospray ionization mode. Nuclear magnetic resonance spectra were recorded on Varian Inova (300 MHz for 1H nuclei), on Bruker Ascend (400 MHz for 1H nuclei) and on Anasazi (60 MHz for 1H nuclei) spectrometers. Chemical shifts (δ) are given in parts per million downfield from tetramethylsilane.

6-substituted-4-chloromethylcoumarins (5, 6 and 7)

To a mixture of ethyl 4-chloroacetoacetate (1 mmol) and the corresponding phenol (1 mmol), H2SO4 (1 mL, 98%) was added dropwise. The reaction mixture was stirred for 24 hours at room temperature, and poured onto ice. The resulting precipitate was filtered off and recrystallized from ethanol-dioxane (9:1 v/v) [21].

General procedure for the preparation of compounds 9a-n, 10a-n and 11a-n

In a round bottom flask, a solution of 1 mmol of the 6-substituted-4-chloro methylcoumarin (5, 6 or 7) in THF (2 mL) was stirred at room temperature for 5 minutes for complete solubilization of the solute. Following this, 2 mmol of the corresponding amine (8a-n), 1.5 mmol of potassium iodide and a few droplets of water were added (the latter to achieve the complete solubilization of KI). The reaction mixture was heated to reflux for 1.5 hours and partitioned between a 10% NaOH solution and ethyl acetate (1:1 v/v). The organic phase was dried over MgSO4, filtered and concentrated, and the residue purified by column chromatography with a hexane/ethyl acetate gradient [22, 23].

Anti-Trypanosoma cruzi activity

Trypanosoma cruzi epimastigotes (Tulahuen 2 strain) were grown at 28 °C in an axenic milieu (BHI-Tryptose) supplemented with 5% fetal bovine serum (FBS), as previously described [6, 24]. Cells from a 10-day-old culture (stationary phase) were inoculated into 50 mL of fresh culture milieu to give an initial concentration of 1 × 106 cells per mL. Cell growth was followed by measuring the absorbance of the culture at 600 nm every day. Before inoculation, the milieu was supplemented with a given amount of the coumarin from a freshly prepared stock solution in DMSO (25 mM). The final concentration of DMSO in the culture medium never exceeded 0.4%, and the control was run in the presence of 0.4% DMSO and in the absence of drugs. No effect on epimastigote growth was observed due to the presence of up to 1% DMSO in the culture milieu. Nifurtimox (Nfx) was used as the reference trypanocidal drug. The percentage of growth inhibition (PGI) was calculated as follows: PGI (%) = {1 − [(Ap − A0p) / (Ac − A0c)]} × 100, where Ap = A600nm of the culture containing the studied compound at day 5; A0p = A600nm of the culture containing the studied compound just after addition of the inocula (day 0); Ac = A600nm of the culture in the absence of the studied compound (control) at day 5; A0c = A600nm in the absence of the studied compound at day 0. To determine the IC50 values (50% inhibitory concentrations), parasite growth was followed in the absence (control) and presence of increasing concentrations of the corresponding compound. At day 5, the absorbance of the culture was measured and related to the control. The IC50 value was taken as the concentration of the compound under study necessary to reduce the absorbance ratio to 50%.

Toxicity assays on Caenorhabditis elegans

The nematode strain used was N2 (wild type), originally obtained from the Caenorhabditis Genetics Center (University of Minnesota, Twin Cities, MN, USA), which was maintained on nematode growth medium (NGM) plates seeded with Escherichia coli OP50 at 20°C. Synchronization of C. elegans cultures at the first larval stage (L1) was achieved by washing off the gravid nematodes from the plates into the centrifuge tubes, which were lysed with a bleaching mixture (1% NaOCl; 0.25 M NaOH), followed by flotation on a sucrose solution 30% (m/v) to separate the eggs from the dissolved worms and bacterial debris. The eggs were washed with M9 buffer (0.02 M KH2PO4, 0.04 M Na2HPO4, 0.08 M NaCl, and 0.001 M MgSO4) and allowed to hatch overnight (16 hours) in NGM agar plates without bacteria. After 16 hours it was obtained worms at L1 stage. Then, 1.500 L1 worms were exposed for 30 minutes at 20°C, by constant agitation in a rotator, at crescent doses ranging from 25 to 100 mM of coumarins in a 0.5% NaCl liquid media. Stock solutions of coumarins were made in DMSO, therefore worms treated with 0.5% NaCl and 0.5% DMSO were used as controls. After exposure, worms were washed 3 times with 0.5% NaCl to remove the treatments. Finally, worms were transferred to NGM recovery plates inoculated with Escherichia coli (OP50) for posterior assays.

Lethal dose 50% (LD50) evaluation

The LD50 of coumarins was determined in C. elegans after exposure. The worms were washed three times with NaCl buffer and placed on OP50-seeded NGM plates. The number of surviving worms on each plate was verified 24 hours after exposure. The lethality data was normalized with the controls and presented as a percentage. All the coumarins were tested in three independent experiments with each concentration tested in two replicates within each experiment.

Body area

For the evaluation of body area, 48 hours after exposure the adult worms (20 nematodes per treatment) were photographed and their body contour was measured. For this, the NGM plates were washed with distilled water and the resulting solution was transferred to a centrifuge tube, where the worm was allowed to settle, separating it from the bacteria in suspension. The process was repeated until the solution became clear. After this procedure, 15 μL of the solution with the worms was deposited on a blade covered by agarose and 15 μL of 2.25% levamisole was added. The pictures were acquired from 20 worms per treatment, which were manually contoured and measured with AxioVision software LE version 4.8.2.0 for Windows.

Structure-activity relationship

Pentacle (Molecular Discovery Ltd) use the GRIND approach, which calculates 3D descriptors in an alignment-independent way [25]. Considering the results obtained in the evaluation of biological activity, it was decided to carry out a binary QSAR study based on the approach of [26] compounds with PGI less than 50%, which was set as inactive (1) and with PGI greater or equal to 50%, which was set as active (2). For each compound, 710 GRIND descriptors were used, based on combinations of the following molecular probes: DRY (hydrophobic interactions), O (hydrogen bond acceptor groups), N1 (hydrogen bond donor groups), and TIP (shape descriptor).

The initial set of MIFs was reduced to 60, using the variable selection method available in the Pentacle, Fractional Factorial Design (FDD). This process is performed by constructing models using Partial Least Squares (PLS), where a reduction of data occurs dimensionally. Most of the variations in the dataset are retained, and new variables (latent variables, LV), mutually orthogonal, are built [27]. In this step and in the subsequent ones, the descriptors had to be autoscaled, the basic process of data pre-processing used in QSAR studies [27]. In order to obtain models with simplified interpretation, the generated set went through a new stage of reduction of variables, where descriptors with absolute Pearson’s correlation coefficient (|r|) values with the binary activities greater than 0.2 were maintained, aiming to maintain only the MIF that presented the maximum amount of information relevant for the classification of the compounds. In the sequence, the Ordered Predictors Selection (OPS) variable selection method, available in the QSAR Modeling program, was used [27, 28]. After the variable selection, which also uses the PLS regression method, the model was constructed and refined using the PLS-DA [29, 30]. The aim of this last stage was to maximize the capacity of classification of the model. The necessary threshold adopted was 0.5, a value generally employed with PLS-DA, where compounds below this value were classified as active, and compounds above this were inactive. One of the advantages of using PLS-DA compared to other classification methods is that validation tools, such as calibration models, can be used. Thus, the quality of the data adjustment was assessed based on its coefficient of determination (R2) and the root mean square error of calibration (RMSEC). The results obtained were from cross-validation (using the values of Q2LOO and RMSECV), and by the visual inspection of the separation of the compounds in the two analyzed classes [31,32,33].

Results

Chemistry

The original and green synthetic route for the 4-methylamino-coumarins is outlined in Fig. 2. Those of the 6-substituted-4-methylchlorochromen-2-ones (compounds 5, 6, and 7) were synthesized using the Pechmann reaction [21, 34], where the p-substituted phenols 1-3, undergo electrophilic aromatic substitution under acid catalysis with an ethyl 4-chloroacetoacetate 4, to give an intermediate that cyclizes by a transesterification reaction to produce the coumarin intermediates 5, 6 and 7 with good yields (71-89%). The N-alkylation reaction of primary and secondary amines (8a-n) was then achieved employing coumarins 5-7 and stoichiometric amounts of potassium iodide in THF/water [22, 23].

Fig. 2
figure2

Reagents and conditions: (a) H2SO4, rt., overnight; (b) KI, THF-Water, 50 °C, 1.5 h

For the reactions, 8 aromatic amines with different substituents and 6 different cyclic aliphatic amines were used in order to structure-activity relationship studies. Additional file 1: Figures S1, S2 and S3, show the formed compounds and their yields.

Biological studies

Anti-Trypanosoma cruzi activity

Considering the potential existence of an intracellular epimastigote stage of T. cruzi as an intermediate between amastigote and trypomastigote forms [35, 36], as well as the readiness of experimentation, the forty-two 4-methylamino-coumarins were assayed against the epimastigote form of T. cruzi, Tulahuen 2 strain. In accordance with (Figs. 3, 4 and 5), in general, from the forty-two coumarins tested, fifteen showed a PGI higher than 50% against T. cruzi epimastigotes. The IC50 values of these 15 molecules ranged from 6.4 to > 25 μM, whereas Nifurtimox showed a IC50 equals to 5.0 ± 1.0 μM.

Fig. 3
figure3

Effect of 4-methylamino-coumarins (9 h series) on the growth of Tripanosoma cruzi epimastigotes

Fig. 4
figure4

Effect of 4-methylamino-coumarins (10 h series) on the growth of Tripanosoma cruzi epimastigotes

Fig. 5
figure5

Effect of 4-methylamino-coumarins (11 h series) on the growth of Tripanosoma cruzi epimastigotes

Toxicity assays on Caenorhabditis elegans

The acute toxicity test was conducted using six different concentrations of coumarins and Nifurtimox. The results are represented by the concentration-response curves for each molecule, and are described in Fig. 6. These LD50 values show that the 4-methylcoumarins 11 h (LD50 = 73.4 mM), 11 g (LD50 = 61.7 mM) and 10 g (LD50 = 42.3 mM) present less toxicity that the Nifurtimox (LD50 = 19.50 mM).

Fig. 6
figure6

Lethality of 10 g (a), 11 g (b), 11 h (c) and Nifurtimox (d) after acute exposure of C. elegans. Data were obtained from at least 3 independent experiments performed in duplicate

Another toxicity endpoint used to verify potential toxic effects of the new coumarins was the evaluation of the worm development, which was monitored by measuring body area. In this study, the normal development of C. elegans was affected by acute exposure to Nifurtimox, 11 h, 11 g, and 10 g (Fig. 7). Compared with the control, Nifurtimox (Fig. 7d), a significant reduction of body area was observed at 50, 75 and 100 mM concentrations (p < 0.01, ANOVA/Bonferroni). 11 h (Fig. 7c) presented a significant reduction in the body area at 50 mM concentrations (p < 0.05, ANOVA/Bonferroni). The worms treated with 11 g (Fig. 7b) showed a significant (p < 0.05, ANOVA/Bonferroni) reduction in body area in the concentrations of 25, 50, 75 and 100 mM. For 10 g (Fig. 7a), a significant (p < 0.05, ANOVA/Bonferroni) reduction of body area was observed for the concentrations of 50, 75 and 100 mM.

Fig. 7
figure7

Body areas of C. elegans after acute treatment (30 min) with 10 g, (a) 11 g (b), 11 h (c) and Nifurtimox (d). Results are expressed as mean ± SEM (n = 3 independent experiments performed in duplicate). *p < 0.05, **p < 0.01, ***p < 0.001 compared to the control group. Statistical comparisons were made using one way ANOVA/Bonferroni post-hoc test

Structure-activity relationship

Using the FDD, a 3D model with 60 MIF descriptors was obtained. Following this and after the exclusion of the descriptors with |r| < 0.2 with the binary vector y, the OPS variable selection method was used, generating a PLS model formed by only 14 descriptors. In the last step, the model was refined in the Pirouette program from the analysis of the regression vector, correlogram and loadings plot obtained using the PLS-DA method. Finally, a model was formed consisting of only 11 MIF descriptors (equation) and four latent variables that encode 60.522% of information (LV1: 16.548%; LV2: 15.633%; LV3: 8.868%; LV4: 19.473%).

$$ \mathbf{Class}=0.405+1.182\ast \left(8\_\mathrm{DRY}-\mathrm{DRY}\right)+1.691\ast \left(20\_\mathrm{DRY}-\mathrm{DRY}\right)+2.574\ast \left(151\_\mathrm{TIP}-\mathrm{TIP}\right)+0.966\ast \left(163\_\mathrm{TIP}-\mathrm{TIP}\right)-0.638\ast \left(169\_\mathrm{TIP}-\mathrm{TIP}\right)-1.098\ast \left(170\_\mathrm{TIP}-\mathrm{TIP}\right)-0.756\ast \left(193\_\mathrm{DRY}-\mathrm{O}\right)-1.086\ast \left(261\_\mathrm{DRY}-\mathrm{N}1\right)+1.047\ast \left(285\_\mathrm{DRY}-\mathrm{TIP}\right)+0.716\ast \left(455\_\mathrm{N}1-\mathrm{TIP}\right)+0.640\ast \left(456\_\mathrm{N}1-\mathrm{TIP}\right)\ {\mathbf{R}}^{\mathbf{2}}=0.737;\mathbf{RMSEC}=0.242;\mathbf{Q}\mathbf{2}\mathbf{LOO}=0.613;\mathbf{RMSEC}\mathbf{V}=0.293. $$

Table 1 shows the degrees of importance of each descriptor. It is possible to note that four descriptors are related to probe TIP-TIP, two to N1-TIP, one to DRY-TIP, two to DRY-DRY, one to DRY-O, and one to DRY-N1.

Table 1 Importance of each GRIND descriptor of the Model

Discussion

The experimental design was performed in order to demonstrate the in vivo alternative model C. elegans applied to a drug discovery context. In this paper, it was not only possible to conceive the C. elegans model, but also it was possible to highlight new compounds with anti-Trypanosoma cruzi activity.

Our workflow started with the synthesis of a focused library from the reaction of 3 coumarins × 14 amines resulting in 42 compounds (9a-9n; 10a-10n; 11a-11n), with good molecular diversity. The choices were: 1. To link non-aromatic (primary and secondary amines) and aromatic (several patterns of substitution) amines to 4-methyl position; 2. To check the importance or not of a tertiary amine in this position; for this purpose, we have compounds with an additional methyl group attached to nitrogen; and 3. Substitution at the C-6 position of coumarin nucleus with groups with different electronegativity or volume (ethyl, methyl or Cl).

In the next step, we subjected the compounds to in vitro inhibition assays using the epimastigote form. The compounds with the highest inhibition percentage (>82%) have aromatic amines bearing the coumarin core and, among them, three compounds with no substitution at the benzene ring (9a, 9g, and 10g) and one with 4-Cl substitution (11h). On the other hand, four coumarins presented IC50 < 10 μM (9a, 10g, 11g (LaSOM 231) and 11h (LaSOM 266). One of these compounds, coumarin 11h, presented an IC50 value (6.4 ± 1 μM) very close to Nifurtimox. Although only Nifurtimox was used as the positive control in this study, it was previously determined that Nifurtimox and Benznidazole have very similar IC50 values when tested in vitro against the epimastigote form of T. cruzi, Tulahuen 2 strain [37]. Indeed, it can be observed that the three most active 11h (IC50 = 6.4 ± 1 μM), 11g (IC50 = 8.2 ± 2.3 μM) and 10g (IC50 = 8.9 μM) have the methyl group attached to nitrogen. With regard to the modifications at position 6 of the coumarin core, compounds with the ethyl group were among the most active, while Cl at this position furnished the less active compounds.

A recent study evaluated 19 tricyclic coumarins through a modeling study and in vitro study similar to our work. The study demonstrated for most coumarins evaluated promising activity against the intracellular forms of T. cruzi, with ten compounds showing IC50 values lower than Benznidazole [38].

Following the current trends in medicinal chemistry for using the alternative C. elegans model in the early stages of the discovery of new drugs [39, 40], our research group recently published the use of this method to evaluate the safety of new potentially antitumor compounds [41].

In this study, the evaluation was made of the toxicity of new compounds with potential activity against Chagas disease, using this alternative method. In order to then check the in vivo toxicity of the three hits 11h, 11g and 10g, selected above, the LD50 investigation of these compounds was performed as studies have shown that C. elegans have a good correlation with rodent oral LD50 ranking. To the best of our knowledge this is the first report of the LD50 of Nifurtimox (Nfx) performed in C. elegans. With this propose in mind, the three new coumarins 11h, 11g, 10g and the standard drug Nifurtimox, were evaluated in C. elegans. According to the literature, Nifurtimox has presented several toxic effects, such as: neurotoxicity, testicular damage, ovarian toxicity and deleterious effects in adrenal, colon, oesophageal and mammary tissue [42]. It appeared to be important to know if the two more active coumarins had the same toxicity of Nifurtimox; therefore, the coumarins and the standard drug were evaluated in order to compare the toxicities among them. Of likely importance for the future is the use of C. elegans as a model system for pre-screening drug discovery. It is to be hoped that this will provide new anthelmintics that less toxic [43, 44]. C. elegans, which is about 1 mm in length as an adult, can be cultured in high-throughput format for multiple generations, allowing the identification of molecules that perturb the worm at any point during its life cycle [45, 46].

A study published by a group that used C. elegans as a model for the discovery of new anthelmintic revealed that by first screening in C. elegans, they may have lost molecules that would be effective in killing parasites but are ineffective in model C. elegans. However, they concluded that the speed and ease with which molecules can be traced using C. elegans can overcome the disadvantages it carries as a primary screening system and that it can be used to evaluate hundreds of thousands of molecules at multiple concentrations at a fraction of the cost and time reduced [47].

A study using albendazole in C. elegans demonstrated the EC50 of 18.43 μM for the worm. This concentration showed damage to the body wall of adults and larvae. It was observed intense desquamation of the cuticle of the larvae and the surface of the eggs, preventing their hatching and development. This study reaffirmed the use of C. elegans as a screening system for compounds with potential anthelmintic activity and showed the effects of albendazole on the different life stages of these worms [48].Considering that the acute toxicity tests were performed with a concentration about 1000 times above those used in the efficacy trials (IC50), these three new coumarins 11 g, 11 h and 10 g can be considered safe in comparison to Nifurtimox. Furthermore, the evaluation of the worm development, expressed by the body area measurement, is a good parameter to evaluate toxic effects in C. elegans, considering that the growth of C. elegans is determined by a conservative genetic regulatory pathway [49].

In order to complement the workflow, we added a QSAR study that, through a mathematical model, emphasizes the importance of certain chemical groups substitutions. The results obtained for the fit (R2) and internal prediction ability (Q2LOO) indicate that the model explains and predicts information at the levels recommended by the literature [31,32,33].

The descriptor analysis emphasizes that the importance of steric characteristics related to TIP descriptors predominates in the model. Despite this, the two most important descriptors (193_DRY-O e 261_DRY-N1) indicate that the presence of hydrophobic groups and groups capable of forming hydrogen bonds at the indicated distances tend to lead to molecules with a degree of inhibition below 50%.

It can be seen that for both descriptors this may mean that the pharmacophoric point formed by hydrophobic groups (R1 position) and electronegative atoms (R2 position) is detrimental to activity as observed in compounds 9c, 9d, 10c, 10d and 11d. On the other hand, volume at the R2 position leads to the most important steric descriptor of the model (151_TIP-TIP), which has a positive influence on the activity, indicating that the substituents present in this position may be filling a bulky cavity at the binding site as can be noticed in compounds 11 g and 11 h. The descriptors 169_TIP-TIP and 170_TIP-TIP show that the positioning of bulky groups at a distance of 12.4 to 13.2 (distances formed by the distances between groups R1 and R2) are unfavourable to the activity. The ideal distance between these two groups is given by the descriptor 163_TIP-TIP (10.0–10.4 Å). The descriptors 455_N1-TIP and 456_N1-TIP show that a relation between the bulk group used in position R2 and the lactonic oxygen in the distance of 16.4 to 17.2 Å is important for the activity. In addition, it is possible to observe in the selected inactive molecule that this pharmacophoric characteristic is absent as it can be seen in compounds 9e, 9i, 9j, 10i and 10j. The selected DRY-DRY descriptors, although encoding different pharmacophoric points, are both related to the hydrophobicity of the groups used in the R2 position. Finally, similar interpretation can be made for the 285_DRY-TIP descriptor, the less important of the models. This shows that the hydrophobicity and volume used in the R1 and R2 position should also be favoured, and for both they seem to indicate that lower volume hydrophobic substituents can be detrimental to the activity (Fig. 8). Each one of the descriptors of model 1 is presented in Fig. 9. In Fig. 9, the graphs of weights and scores obtained by PLS-DA are presented, showing the classification of the data set obtained through latent variables 1 and 2, and the distribution of the descriptors in relation to each of these variables. This result indicates that the obtained model is capable of discriminating with good reliability those compounds having a degree of inhibition below or above 50%, and thus the results obtained herein may be useful to direct the synthesis of new derivatives with an acceptable degree of inhibition.

Fig. 8
figure8

GRIND selected descriptors of Model 1 associated with active compounds

Fig. 9
figure9

Plot of the loading and score vectors

C. elegans has proved to be an extremely useful model organism for toxicity studies for pre-screening new drugs. We can observe as limitations of the model used in our study: C. elegans do not possess some mammalian organs, such as lungs, heart, kidneys and liver. In addition, changes in temperature and humidity may alter test results. Another limitation of the study was to respect the solubility limit of the molecules we evaluated, since if we increase the concentrations we could have problems of solubility. Therefore, higher concentrations than those described in the study were not evaluated.

Although it was constructed from coumarin-like structures, the QSAR model was based on molecular interaction fields, which have established the best distances and interaction groups related to the activity. These observations may guide the synthesis of new compounds with better profile against T. cruzi. Nevertheless, the constructed model is limited by the PGI values and by the cutoff point settled to obtain the classificatory model, which is helpful to select features correlated to the activity. However, it lacks information regarding quantitative values, thus turning the data extrapolation difficult to be predicted.

Conclusions

The synthesis of a library of forty two 4-methylamino-coumarins with aromatic and aliphatic rings attached to the amino groups was accomplished using green chemistry conditions. The library activity was assayed in vitro against T. cruzi epimastigotes and two of these coumarins demonstrated to be nearly equipotent to Nifurtimox (IC50 = 5.0 ± 1 μM), 11 h (IC50 = 6.4 μM) and 11 g (IC50 = 8.2 μM). Also, toxicity assay performed on C. elegans showed that these two compounds 11 h (DL50 = 70.81 mM) and 11 g (DL50 = 61.50 mM) were clearly less toxic than Nifurtimox (DL50 = 19.50 mM). In addition, the structure-activity relationship study showed that hydrophobic groups in R1 position and electronegative atoms in R2 position are detrimental to activity. In conclusion, 11 h and 11 g may be a feasible starting point for the development of new trypanocidal compounds. Further studies will be made in order to determine experimentally the mechanism of action of the coumarins using labeled strains of C. elegans.

Availability of data and materials

Not applicable.

References

  1. 1.

    Organization WH. Chagas disease (American trypanosomiasis). 2017.

    Google Scholar 

  2. 2.

    Urbina JA, Payares G, Sanoja C, Molina J, Lira R, Brener Z, Romanha AJ. Parasitological cure of acute and chronic experimental Chagas disease using the long-acting experimental triazole TAK-187. Activity against drug-resistant Trypanosoma cruzi strains. Int J Antimicrob Agents. 2003;21(1):39–48.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Arias AR. Situación de la enfermidad de Chagas en el Continente Americano, in Enfermidad de Chagas: Estrategias em la búsqueda de nuevos medicamentos. Mexico,DF: RIDIMEDCHAG-CYTED; 2012.

    Google Scholar 

  4. 4.

    Trouiller P, Olliaro P, Torreele E, Orbinski J, Laing R, Ford N. Drug development for neglected diseases: a deficient market and a public-health policy failure. Lancet. 2002;359(9324):2188–94.

    PubMed  Article  Google Scholar 

  5. 5.

    Graebin C, Uchoa F, Bernardes L, Campo V, Carvalho I, Eifler-Lima V. Antiprotozoal Agents: An Overview. AntiInfect Agents Med Chem. 2009;8(4):345–66.

    CAS  Article  Google Scholar 

  6. 6.

    Gomez-Outes A, Luisa Suarez-Gea M, Calvo-Rojas G, Lecumberri R, Rocha E, Pozo-Hernandez C, Isabel Terleira-Fernandez A, Vargas-Castrillon E. Discovery of Anticoagulant Drugs: A Historical Perspective. Curr Drug Discov Technol. 2012;9(2):22.

    Google Scholar 

  7. 7.

    Joao Matos M, Vazquez-Rodriguez S, Santana L, Uriarte E, Fuentes-Edfuf C, Santos Y, Munoz-Crego A. Looking for new targets: simple Coumarins as antibacterial agents. Med Chem. 2012;8(6):1140–5.

    Google Scholar 

  8. 8.

    Vianna DR, Bubols G, Meirelles G, Silva BV, da Rocha A, Lanznaster M, Monserrat JM, Garcia SC, von Poser G, Eifler-Lima VL. Evaluation of the antioxidant capacity of synthesized coumarins. Int J Mol Sci. 2012;13(6):7260–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Vianna DR, Hamerski L, Figueiró F, Bernardi A, Visentin LC, Pires ENS, Teixeira HF, Salbego CG, Eifler-Lima VL, Battastini AMO, Von Poser GL, Pinto AC. Selective cytotoxicity and apoptosis induction in glioma cell lines by 5-oxygenated-6,7-methylenedioxycoumarins from Pterocaulon species. Eur J Med Chem. 2012;57:268–74.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Evans BE, Rittle KE, Bock MG. Methods for drug discovery: development of potent, selective, orally effective cholecystokinin antagonists. J Med. 1988;31(12):2235–46.

    CAS  Google Scholar 

  11. 11.

    Vianna DR, Ruschel L, Dietrich F, Figueiró F, Morrone FB, Canto RFS, Corvello F, Velho A, Crestani A, Teixeira H, von Poser GL, Battastini AMO, Eifler-Lima VL. 4-Methylcoumarins with cytotoxic activity against T24 and RT4 human bladder cancer cell lines. Med Chem Commun. 2015;6(5):905–11.

    CAS  Article  Google Scholar 

  12. 12.

    Cavalli A, Lizzi F, Bongarzone S, Brun R, Luise Krauth-Siegel R, Bolognesi ML. Privileged structure-guided synthesis of quinazoline derivatives as inhibitors of trypanothione reductase. Bioorg Med Chem Lett. 2009;19(11):3031–5.

    CAS  PubMed  Article  Google Scholar 

  13. 13.

    Torres FC, Brucker N, Andrade SF, Kawano DF, Garcia SC, Poser GL, Eifler-Lima VL. New insights into the chemistry and antioxidant activity of coumarins. Curr Top Med Chem. 2014;14(22):2600–23.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Peng X-M, L.V. Damu G, He Zhou C. Current developments of Coumarin compounds in medicinal chemistry. Curr Pharm Des 2013;19(21):3884–3930.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Vieira PC, Mafezoli J, Pupo MT, Fernandes JB, da Silva MFDGF, de Albuquerque S, Oliva G, Pavão F. Strategies for the isolation and identification of trypanocidal compounds from the Rutales. Pure Appl Chem. 2001;73(3):617–22.

    CAS  Article  Google Scholar 

  16. 16.

    Brak K, Doyle PS, McKerrow JH, Ellman JA. Identification of a new class of nonpeptidic inhibitors of cruzain. J Am Chem Soc. 2008;130(20):6404–10.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Oliaro-Bosso S, Taramino S, Viola F, Tagliapietra S, Ermondi G, Cravotto G, Balliano G. Umbelliferone aminoalkyl derivatives as inhibitors of human oxidosqualene-lanosterol cyclase. J Enzyme Inhib Med Chem. 2009;24(2):589–98.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Rea A, Tempone AG, Pinto EG, Mesquita JT, Rodrigues E, Silva LGM, Sartorelli P, Lago JHG. Soulamarin Isolated from Calophyllum brasiliense (Clusiaceae) Induces Plasma Membrane Permeabilization of Trypanosoma cruzi and Mytochondrial Dysfunction. PLoS Neglected Trop Dis. 2013;7(12):e2556.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  19. 19.

    Oliaro-Bosso S, Viola F, Taramino S, Tagliapietra S, Barge A, Cravotto G, Balliano G. Inhibitory effect of Umbelliferone Aminoalkyl derivatives on Oxidosqualene Cyclases from S. cerevisiae, T. cruzi, P. carinii, H. sapiens, and A. thaliana: a structure–activity study. ChemMedChem. 2007;2(2):226–33.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  20. 20.

    Armarego WLF, Chai CLL. Purification of Laboratory Chemicals: Fifth Edition; 2003. p. 1–609.

    Google Scholar 

  21. 21.

    Frasinyuk MS, Vinogradova VI, Bondarenko SP, Khilya VP. Synthesis of cytisine derivatives of coumarins. Chem Nat Compd. 2007;43(5):590–3.

    CAS  Article  Google Scholar 

  22. 22.

    Li X, Mintz EA, Bu XR, Zehnder O, Bosshard C, Günter P. Phase transfer catalysis for tandem alkylation of azo dyes for the synthesis of novel multifunctional molecules. Tetrahedron. 2000;56(32):5785–91.

    CAS  Article  Google Scholar 

  23. 23.

    Loeser E, Prasad K, Repic O. Selective N-alkylation of primary amines with chloroacetamides under pH-controlled aqueous conditions. Synth Commun. 2002;32(3):403–9.

    CAS  Article  Google Scholar 

  24. 24.

    Varela J, Lavaggi ML, Cabrera M, Rodríguez A, Miño P, Chiriboga X, Cerecetto H, González M. Bioactive-guided identification of labdane diterpenoids from aerial parts of Aristeguietia glutinosa as anti-Trypanosoma cruzi agents. Nat Prod Commun. 2012;7(9):1139–42.

    CAS  PubMed  Google Scholar 

  25. 25.

    Pastor M, Cruciani G, McLay I, Pickett S, Clementi S. GRid-INdependent descriptors (GRIND): a novel class of alignment-independent three-dimensional molecular descriptors. J Med Chem. 2000;43(17):3233–43.

    CAS  PubMed  Article  Google Scholar 

  26. 26.

    Ermondi G, Caron G, Pintos IG, Gerbaldo M, Pérez M, Pérez DI, Gándara Z, Martínez A, Gómez G, Fall Y. An application of two MIFs-based tools (Volsurf+ and pentacle) to binary QSAR: the case of a palinurin-related data set of non-ATP competitive glycogen synthase kinase 3β (GSK-3β) inhibitors. Eur J Med Chem. 2011;46(3):860–9.

    CAS  PubMed  Article  Google Scholar 

  27. 27.

    Teófilo RF, Martins JPA, Ferreira MMC. Sorting variables by using informative vectors as a strategy for feature selection in multivariate regression. J Chemom. 2009;23(1):32–48.

    Article  CAS  Google Scholar 

  28. 28.

    Martins JP, Ferreira MMC. Qsar modeling: a new open source computational package to generate and validate Qsar models. Quim Nova. 2013;36(4):554–U250.

    CAS  Article  Google Scholar 

  29. 29.

    Barker M, Rayens W. Partial least squares for discrimination. J Chemom. 2003;17(3):166–73.

    CAS  Article  Google Scholar 

  30. 30.

    Lang KL, Silva IT, Machado VR, Zimmermann LA, Caro MSB, Simões CMO, Schenkel EP, Durán FJ, Bernardes LSC, de Melo EB. Multivariate SAR and QSAR of cucurbitacin derivatives as cytotoxic compounds in a human lung adenocarcinoma cell line. J Mol Graph Model. 2014;48:70–9.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Roy PP, Roy K. On some aspects of variable selection for partial least squares regression models. QSAR Comb Sci. 2008;27(3):302–13.

    CAS  Article  Google Scholar 

  32. 32.

    Camargo LTFM, Sena MM, Camargo AJ. A quantum chemical and chemometrical study of indolo [2,1-b] quinazoline and their analogues with cytotoxic activity against breast cancer cells. SAR QSAR Environ Res. 2009;20(5–6):537–49.

    CAS  PubMed  Article  Google Scholar 

  33. 33.

    Dejaegher B, Dhooghe L, Goodarzi M, Apers S, Pieters L, Vander HY. Classification models for neocryptolepine derivatives as inhibitors of the β-haematin formation. Anal Chim Acta. 2011;705(1–2):98–110.

    CAS  PubMed  Article  Google Scholar 

  34. 34.

    Russell A, Frye JR. 2,6-Dihydroxyacetophenone. Organic Syntheses. 1941;21:22.

  35. 35.

    Faucher JF, Baltz T, Petry KG. Detection of an "epimastigote-like" intracellular stage of Trypanosoma cruzi. Parasitol Res. 1995;81(5):441–3.

    CAS  PubMed  Article  Google Scholar 

  36. 36.

    Almeida-de-Faria M, Freymüller E, Colli W, Alves MJM, Freymu E. Trypanosoma cruzi: characterization of an intracellular epimastigote-like form. Exp Parasitol. 1999;92:263–74.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Gerpe A, Álvarez G, Benítez D, Boiani L, Quiroga M, Hernández P, Sortino M, Zacchino S, González M, Cerecetto H. 5-Nitrofuranes and 5-nitrothiophenes with anti-Trypanosoma cruzi activity and ability to accumulate squalene. Bioorg Med Chem. 2009;17(21):7500–9.

    CAS  PubMed  Article  Google Scholar 

  38. 38.

    Coelho GS, Andrade JS, Xavier VF, Sales Junior PA, Rodrigues de Araujo BC, Fonseca KS, Caetano MS, Murta SMF, Vieira PM, Carneiro CM, Taylor JG. Design, synthesis, molecular modelling, and in vitro evaluation of tricyclic coumarins against Trypanosoma cruzi. Chem Biol Drug Des. 2019;93(3):337–50.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Salgueiro WG, Xavier MCDF, Duarte LFB, Câmara DF, Fagundez DA, Soares ATG, Perin G, Alves D, Avila DS. Direct synthesis of 4-organylsulfenyl-7-chloro quinolines and their toxicological and pharmacological activities in Caenorhabditis elegans. Eur J Med Chem. 2014;75:448–59.

    CAS  PubMed  Article  Google Scholar 

  40. 40.

    Carretero M, Solis GM, Petrascheck M. C. elegans as model for drug discovery. Curr Top Med Chem. 2017;17(18):2067–76.

    CAS  PubMed  Article  Google Scholar 

  41. 41.

    Gonçalves IL, Rockenbach L, das Neves GM, Göethel G, Nascimento F, Porto Kagami L, Figueiró F, Oliveira de Azambuja G, de Fraga Dias A, Amaro A, de Souza LM, da Rocha Pitta I, Avila DS, Kawano DF, Garcia SC, Battastini AMO, Eifler-Lima VL. Effect of N-1 arylation of monastrol on kinesin Eg5 inhibition in glioma cell lines. MedChemComm. 2018;9(6):995–1010.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Ja C, de Mecca MM, Bartel LC. Toxic side effects of drugs used to treat Chagas' disease (American trypanosomiasis). Hum Exp Toxicol. 2006;25(8):471–9.

    Article  Google Scholar 

  43. 43.

    Behm CA, Bendig MM, McCarter JP, Sluder AE. RNAi-based discovery and validation of new drug targets in filarial nematodes. Trends Parasitol. 2005;21(3):97–100.

    CAS  PubMed  Article  Google Scholar 

  44. 44.

    Jones AK, Buckingham SD, Sattelle DB. Chemistry-to-gene screens in Caenorhabditis elegans. Nat Rev Drug Discov. 2005;4(4):321–30.

    CAS  PubMed  Article  Google Scholar 

  45. 45.

    Lemieux GA, Liu J, Mayer N, Bainton RJ, Ashrafi K, Werb Z. A whole-organism screen identifies new regulators of fat storage. Nat Chem Biol. 2011;7:206.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Petrascheck M, Ye X, Buck LB. An antidepressant that extends lifespan in adult Caenorhabditis elegans. Nature. 2007;450:553.

    CAS  PubMed  Article  Google Scholar 

  47. 47.

    Burns AR, Luciani GM, Musso G, Bagg R, Yeo M, Zhang Y, Rajendran L, Glavin J, Hunter R, Redman E, Stasiuk S, Schertzberg M, Angus McQuibban G, Caffrey CR, Cutler SR, Tyers M, Giaever G, Nislow C, Fraser AG, MacRae CA, Gilleard J, Roy PJ. Caenorhabditis elegans is a useful model for anthelmintic discovery. Nat Commun. 2015;6:7485.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  48. 48.

    Sant’anna V, Vommaro RC, de Souza W. Caenorhabditis elegans as a model for the screening of anthelminthic compounds: Ultrastructural study of the effects of albendazole. Exp Parasitol. 2013;135(1):1–8.

    PubMed  Article  CAS  Google Scholar 

  49. 49.

    Wu Q, Nouara A, Li Y, Zhang M, Wang W, Tang M, Ye B, Ding J, Wang D. Comparison of toxicities from three metal oxide nanoparticles at environmental relevant concentrations in nematode Caenorhabditis elegans. Chemosphere. 2013;90(3):1123–31.

    CAS  PubMed  Article  Google Scholar 

Download references

Acknowledgements

The authors wish to thank the Brazilian funding agencies CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior), CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico), FAPERGS (Fundação de Amparo à Pesquisa do Estado do Rio Grande do Sul), and INCT-IF (Instituto Nacional de Ciência e Tecnologia para Inovação Farmacêutica) for their financial support.

About this supplement

This article has been published as part of BMC Pharmacology and Toxicology Volume 20 Supplement 1, 2019: Proceedings of Toxi-Latin 2018. The full contents of the supplement are available online at https://bmcpharmacoltoxicol.biomedcentral.com/articles/supplements/volume-20-supplement-1.

Funding

Publciation costs were funded by the Brazilian agencies: Conselho Nacional de Desenvolvimento Cientifico e Tecnológico (CNPq), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001 and Instituto Nacional de Ciência e Tecnologia para Inovação Farmacêutica (INCT-IF/CNPq).

Author information

Affiliations

Authors

Contributions

FGNS and GG contributed equally to this article. FGNS designed, performed and interpreted the data of all the synthesis and T. cruzi experiments, also prepared the initial draft of the manuscript. GG designed, performed and interpreted the data of all the C. elegans toxicity assays and carried out the statistical analysis, prepared the subsequent drafts and also critically revised the manuscript. LPK and GMdN revised the structure-activity relationship studies in order to link the synthesis and biological activity findings. ES critically revised the toxicity assays and the statistical analysis. EB and JV performed the T. cruzi experiments. ILG critically revised all the synthesis procedures and spectra, also carried out the statistical analysis. FRP and EBdM designed and performed the structure-activity relationship studies. GVP and VLEL mentored the synthesis section. MG and HC mentored the T. cruzi experiments and checked the experimental findings. DFK mentored the structure-activity relationship studies. SCG mentored the toxicity assays on C. elegans. VLEL and SCG critically revised all the drafts for important intellectual content. All authors read and approved of the final manuscript.

Corresponding author

Correspondence to Vera Lucia Eifler-Lima.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1:

Novel coumarins active against Trypanosoma cruzi and toxicity assessment using the animal model Caenorhabditis elegans.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Soares, F.G.N., Göethel, G., Kagami, L.P. et al. Novel coumarins active against Trypanosoma cruzi and toxicity assessment using the animal model Caenorhabditis elegans. BMC Pharmacol Toxicol 20, 76 (2019). https://doi.org/10.1186/s40360-019-0357-z

Download citation

Keywords

  • Caenorhabditis elegans
  • Trypanosoma cruzi
  • Coumarins
  • Chagas disease
  • Structure-activity relationship