Next Article in Journal
Cellulose Modified with Polyethylenimine (PEI) Using Microwave Methodology for Adsorption of Chromium from Aqueous Solutions
Previous Article in Journal
Odor Perception and Descriptions of Rose-Scented Geranium Pelargonium graveolens ‘Dr. Westerlund’–Sensory and Chemical Analyses
Previous Article in Special Issue
Green Synthesis of Silver Nanoparticles from Allium cepa L. Peel Extract, Their Antioxidant, Antipathogenic, and Anticholinesterase Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Piplartine-Inspired 3,4,5-Trimethoxycinnamates: Trypanocidal, Mechanism of Action, and In Silico Evaluation

by
Carlos S. M. B. Filho
1,
Ramon R. P. P. B. de Menezes
2,
Emanuel P. Magalhães
2,
Yunierkis P. Castillo
3,
Alice M. C. Martins
2 and
Damião P. de Sousa
1,*
1
Department of Pharmaceutical Sciences, Federal University of Paraíba, João Pessoa 58051-900, PB, Brazil
2
Department of Clinical and Toxicological Analysis, Federal University of Ceará, Fortaleza 60020-181, CE, Brazil
3
Escuela de Ciencias Físicas y Matemáticas, Universidad de Las Américas, Quito 170125, Ecuador
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(11), 4512; https://doi.org/10.3390/molecules28114512
Submission received: 24 April 2023 / Revised: 25 May 2023 / Accepted: 29 May 2023 / Published: 2 June 2023

Abstract

:
Chagas disease (CD) is one of the main neglected tropical diseases that promote relevant socioeconomic impacts in several countries. The therapeutic options for the treatment of CD are limited, and parasite resistance has been reported. Piplartine is a phenylpropanoid imide that has diverse biological activities, including trypanocidal action. Thus, the objective of the present work was to prepare a collection of thirteen esters analogous to piplartine (113) and evaluate their trypanocidal activity against Trypanosoma cruzi. Of the tested analogues, compound 11 ((E)-furan-2-ylmethyl 3-(3,4,5-trimethoxyphenyl)acrylate) showed good activity with IC50 values = 28.21 ± 5.34 μM and 47.02 ± 8.70 μM, against the epimastigote and trypomastigote forms, respectively. In addition, it showed a high rate of selectivity to the parasite. The trypanocidal mechanism of action occurs through the induction of oxidative stress and mitochondrial damage. In addition, scanning electron microscopy showed the formation of pores and leakage of cytoplasmic content. Molecular docking indicated that 11 probably produces a trypanocidal effect through a multi-target mechanism, including affinity with proteins CRK1, MPK13, GSK3B, AKR, UCE-1, and UCE-2, which are important for the survival of the parasite. Therefore, the results suggest chemical characteristics that can serve for the development of new trypanocidal prototypes for researching drugs against Chagas disease.

1. Introduction

Chagas disease (CD) is an anthropozoonosis caused by Trypanosoma cruzi that has arthropods of the subfamily Triatominae as vectors [1]. Currently, according to the World Health Organization, CD is considered one of the main neglected tropical diseases that promotes relevant socioeconomic impacts in several countries, mainly in Latin America [2]. Data from the Pan American Health Organization (PAHO) estimate that 70 million people in Latin America are at risk of infection with T. cruzi and that 6–7 million individuals are infected, causing 14,000 deaths annually as a result of the complications of the disease [3].
The therapeutic options for the treatment of Chagas disease are limited, and there are no new drugs for the treatment since the treatment consists of the use of drugs developed more than 50 years ago (benznidazole and nifurtimox), which are more effective in the acute phase of the disease [4]. Furthermore, these drugs are highly toxic, and the use of nifurtimox has been suspended in some countries, and parasite resistance to these drugs has been reported [5,6]. Thus, the search for new drugs with low toxicity and preferably with new pharmacological targets becomes relevant [7,8].
Natural bioactive compounds are important sources of drug candidates [9]. Piplartine, also called piperlongumine (Figure 1), is a phenylpropanoid in the form of an alkamide typically found in plants of the genus Piper and that has diverse biological activities, including cytotoxic, antitumor, schistoicidal, leishmanicidal, and trypanocidal actions [10,11,12,13,14,15,16].
Cell viability studies revealed that piplartine showed greater inhibition of the growth of epimastigote forms of T. cruzi (IC50 = 21.0 μM) than benznidazole (42.7 μM) [17]. In addition, piplartine reduces the expression of antioxidant enzymes involved in the defense of the parasite, such as tryparedoxine reductase and methionine sulfoxide reductase, which is a possible mechanism of action [16]. Furthermore, esters derived from cinnamic acid derivatives, the precursor of piplartine, showed trypanocidal activity against the epimastigote and trypomastigote forms of T. cruzi, especially the compound pentyl p-coumarate, which presented IC50 of 5.16 and 61.63 µM against the epimastigote and trypomastigote forms, respectively [18]. Considering the trypanocidal potential of piplartine and the fact that there are few data about the trypanocidal action of its derivatives, the objective of this work was to prepare a collection of compounds analogous to piplartine for the evaluation of trypanocidal activity.

2. Results and Discussion

2.1. Preparation of Esters 113

In this study, thirteen esters derived from 3,4,5-trimethoxycinnamic acid were prepared by maintaining the 3,4,5-trimethoxycinnamate substructure and differing the R moiety with alkyl and aryl substituents. These esters were obtained in a single step and through three methodologies: Fischer esterification [19] for esters: 1 (R = methyl), 2 (R = ethyl), 3 (R = propyl), and 4 (R = isopropyl), which showed yields from 58.8–91.1%; esterification with alkyl and aryl halides [11]: 5 (R = pentyl), 6 (R = decyl), and 7 (R = 4-methoxybenzyl) with a yield of 38.0–47.5%; and Steglich esterification [20]: 8 (R = 3-methoxybenzyl), 9 (R = 4-methylbenzyl), 10 (R = benzyl), 11 (R = furfuryl), 12 (R = (−)-perillyl) and 13 (R = (−)-bornyl) with yields ranging from 26.7–62.8%. Scheme 1 illustrates these reactions.
All compounds were characterized by infrared spectroscopy (IR) and nuclear magnetic resonance of hydrogen (1H NMR) and carbon thirteen (13C NMR). Based on the 1HNMR spectrum, the compounds showed an approximate purity greater than 95%. In addition, compounds 8, 9, and 12 are unpublished in the literature, and analysis was performed using high resolution mass spectrometry (HRMS). The results and references used in the synthesis of the compounds already described in the literature are available in the Supplementary Materials.

2.2. Trypanocidal Evaluation of Esters against T. cruzi

Initially, the prepared esters were submitted to in vitro tests using the microdilution technique in 96-well plates to evaluate the trypanocidal activity against the epimastigote evolutionary form at different concentrations (200–3.12 μg/mL). Then, the compounds with the best activity were tested against the trypomastigote form of T. cruzi. In both assays, IC50 values were calculated; a compound is considered to have trypanocidal activity when IC50 <40 μM [21]. The results are expressed in Table 1.
In addition, in vitro toxicity tests were performed on renal epithelial cells obtained from monkey LLC-MK2 (Rhesus monkey kidney) using different concentrations of the compounds (200–25 μg/mL), and the selectivity index (SI) was calculated using the ratio between the CC50 (concentration capable of promoting cytotoxicity in 50% of the LLC-MK2 host cells) and the IC50. The compounds were considered selective for the parasite when SI ≥ 10 [21]. The results are expressed in Table 2.
Compound 1 presented an IC50 of 124.9 ± 36.6 μM against the epimastigote form of the parasite, and comparing it with the analogue with 2 (ethyl radical), it is noted that there was a decrease in potency (321.83 ± 72.89 μM) in contrast to the results obtained by Lima et al., in which ethyl caffeate had an IC50 of 18.27 μM against the amastigote form of the parasite [22]. The presence of phenolic hydroxyls in this derivative may have contributed to the better antiparasitic profile compared to 1, in addition to possible differences in the methodological approach of the biological test.
In relation to the analogues 3 (propyl group) and 4 (isopropyl group), which presented IC50 of 321.42 ± 105.95 μM and 265.77 ± 59.22 μM, respectively, comparing these compounds with 1, it is noticed that the increase of the carbonic side chain by up to three carbons drastically reduces the trypanocidal activity of 3,4,5-trimethoxycinnamates. In addition, the use of the decyl group (6) suppressed the trypanocidal activity.
We hypothesize that the decyl group makes the compound highly lipophilic. Consequently, it can remain trapped in the biological membrane and insoluble in the aqueous medium. To further test this hypothesis, the logarithm of the partition coefficient (logP) was predicted for the 13 compounds evaluated in our research with RDKit (RDKit: Open-source cheminformatics. Version: 2023_03_1; RDKit, 2023. http://www.rdkit.org (accessed on 23 May 2023)). The results of the logP predictions are provided as Table S1 in Supplementary Materials and show that compound 6 is the most lipophilic one. In fact, compound 6 is the only one exceeding a logP value of 5.0, which is the limit set by Lipinski’s rule of five for drug-likeness. The atomic contributions to logP were also explored for each compound using RDKit. The results of the later analysis are depicted in Figure S43 of the Supplementary Material. These predictions reveal that the contribution of the 3,4,5-trimethoxycinnamate ester scaffold to logP is conserved among compounds. On the other hand, the aggregated contribution of longer alkyl substituents is related to larger lipophilic contributions, disfavoring the drug-likeness of compound 6 and making it less soluble and more susceptible to getting trapped in biological membranes [23].
However, the increase in the carbon chain with the introduction of the pentyl radical (5) promoted an increase in the antiparasitic action (IC50 = 65.50 ± 7.13 μM) in relation to 1. Similar results were obtained for pentyl p-coumarate (IC50 of 5.16 μM) against the epimastigote form of T. cruzi, compared to methyl p-coumarate (IC50 = 601 μM), as found by Lopes et al. (2019) [18]. Thus, the replacement of the trimethoxylated aromatic ring by a hydroxylated aromatic ring potentiates the trypanocidal activity.
Compound 10 showed an IC50 of 121.66 ± 42.73 μM against the epimastigote form, and when comparing it with substituted compounds 7 (4-methoxybenzyl group), 8 (3-methoxybenzyl group), and 9 (4-methylbenzyl group) (IC50 = 157.37 ± 38.50 μM, 112.36 ± 11.72 μM and 180.49 ± 30.73 μM, respectively), it is noted that there were no significant differences in bioactivity against the epimastigote form of the parasite. However, against the trypomastigote form, the absence of substituents in the aromatic ring resulted in good trypanocidal activity since compound 10 exhibited an IC50 of 40.75 ± 12.36 μM and showed selectivity for the parasite (SI = 13.4).
Compound 11, which has a furfuryl radical, showed the best biological activity against the epimastigote form (IC50 = 28.21 ± 5.34 μM) and trypanocidal activity similar to piplartine (21.0 μM). It is worth noting that both substances have an electronegative atom as a heteroatom in a cyclic system, which can contribute to the potentiation of bioactivity. In this study, the reference drug, benznidazole, had an IC50 of 42.7 μM [17]. Furthermore, 11 showed a high selectivity index (SI > 22.2), and it was not possible to accurately determine the CC50 as the compound did not show toxicity at the highest concentration tested. In addition, compared to trypomastigote cells, the IC50 was 47.02 ± 8.70 μM, with no significant difference compared to 10, also presenting a similar SI. Therefore, further tests were performed with compound 11 to investigate the possible mechanism of trypanocidal activity. Additionally, regarding the investigations of the solubility, stability, and pharmacokinetics of compound 11, there are no studies available that have addressed these aspects. However, due to its being in the chemical class of esters, it is expected to undergo rapid hydrolysis into 3,4,5-trimethoxycinnamic acid and furfuryl alcohol by plasma esterases [24,25]. Furthermore, previous studies conducted with 3,4,5-trimethoxycinnamic acid in rat plasma have indicated that it has rapid oral absorption and high bioavailability [26]. In general, the collection of compounds has a low polarity, so they should cross biological membranes easily due to their lipophilicity.

2.3. Analysis of Cytoplasmic Reactive Oxygen Species (ROS)

To elucidate the mechanism of action of 11, the increase in the concentration of ROS inside epimastigote cells of T. cruzi was evaluated through fluorescent techniques using 2′,7′–dichlorofluorescein diacetate (DCFH-DA) [27]. The results are illustrated in Figure 2 and Figure 3.
From these results, it was observed that there was an increase in relative fluorescence in both treated groups when compared to the control (CT). The result indicates that 11 induces oxidative stress in the parasite through increased levels of ROS, similar to piplartine, which has a trypanocidal effect through the induction of oxidative stress [16].

2.4. Evaluation of Mitochondrial Transmembrane Potential

In addition, the mitochondrial damage exerted by 11 on the parasite was evaluated through the mitochondrial transmembrane potential assay using the fluorescent dye rhodamine 123 (Rho 123) [22,23]. The results are shown in Figure 4 and Figure 5.
It is concluded that 11 promotes a reduction in the accumulation of rhodamine in the intermembrane space of the mitochondria, indicating mitochondrial injury [28,29]. These data explain the possible trypanocidal mechanisms promoted by 11. Similar results were obtained through the treatment of the epimastigotes with pentyl p-coumarate, both cinnamic derivatives [18].

2.5. Structural Evaluation of Epimastigote Forms of T. cruzi Treated with 11

In addition, scanning electron microscopy was used to analyze possible morphological alterations of epimastigote forms treated with compound 11. After 24 h of treatment, ultrastructural alterations were observed in the treated groups in relation to the control, such as pore formation and content leakage in the cytoplasm (Figure 6). Similar results were obtained by treating epimastigote forms with pentyl p-coumarate [18].

2.6. Analysis of the Cell Death Profile of T. cruzi

Epimastigote forms were used to analyze the cell death profile through labeling with AxPE and 7-AAD. Initially, the cells were divided into four cell populations: viable cells—with a low level of labeling (7AAD−/AxPE−); necrotic cells—labeled with 7-AAD (7AAD+/AxPE−); apoptotic cells—labeled with annexin V-PE (7AAD−/AxPE+); and cells in late apoptosis—doubly labeled (7AAD+/AxPE+).
The cells were treated with 11 for 24 h at concentrations of 157 and 314 µM. According to the results (Figure 7), there was a reduction in the percentage of viable cells, an increase in the percentage of necrotic cells, an increase in apoptotic cells, and mainly an increase in cells in late apoptosis.

2.7. Molecular Modeling

Modeling studies were performed with the objective of identifying the most probable molecular targets of compound 11 in T. cruzi. In brief, potential targets for the compound were first identified through computational target fishing predictions. Next, the compound was docked to the proteins retrieved in the previous step, and binding hypotheses for all of them were generated. Finally, the free energies of binding for the docking-predicted complexes were estimated from molecular dynamics (MD) simulations. It must be stressed that the ranking of the most probable targets of compound 11 was performed based on the MD-predicted binding energies. In addition to providing more accurate results than molecular docking in terms of binding energies, free energies calculated from MD trajectories avoid the biases associated with comparing docking scores across different targets.
The objective of computational target-fishing methods is to identify potential targets for chemical compounds. They rely on the similarity principle and use databases containing confirmed ligand-receptor interactions as references. Unfortunately, these databases are biased against human targets. For this reason, the homology-based approach described in the methods section was implemented to identify potential T. cruzi targets for 11. The list of potential molecular targets obtained following this procedure is presented in Table 3. Among the proteins listed in Table 3, TUB-A and TUB-B form a heterodimer that contains the inhibitors’ binding site at the interface of both subunits [30]. Thus, the tubulin (TUB) heterodimer was assembled, taking as reference the crystal structure of the tubulin from Rattus norvegicus in complex with an inhibitor (PDB code 6H9B). 11 was then docked into the dimer interface ligand binding pocket of TUB. In addition, it has been described that the ubiquitin-conjugating enzyme E2 can be inhibited by directly blocking the interaction of ubiquitin with the catalytic cysteine residue as well as by binding to an allosteric site located distant from the catalytic site [31]. For this reason, the binding of 11 to both the catalytic and the allosteric sites of the ubiquitin-conjugating enzyme E2 was explored.
Molecular docking of 11 of the 19 targets listed in Table 3 was performed as described in the Methods section, and its results are summarized in Table 4. In 13 out of the 21 studied systems, more than one probable binding mode is predicted for 11, leading to 42 ligand-receptor complexes for further analyses. According to the obtained results, the best docking scores are observed for the AKR, MetAP2, SIR2, and MPK13 proteins, while the worst-scored targets are MPK4, PK-2, and UCE-1. Despite these differences in scores across targets, the visual inspection of the 42 complexes reveals interactions that could stabilize all these complexes and that, in all cases, the ligand occupies the receptors’ binding sites. Overall, this data is not enough to prioritize a subset of the 11 potential targets studied.
We have previously shown that the additional study of ligand-receptor complexes through MD simulations for the estimation of their free energy of binding can aid the prioritization of the potential molecular targets of chemical compounds [15,18]. Taking this into account, MD simulations and MM-PBSA calculations were performed for the 42 predicted ligand-receptor complexes as described in the methods section, and their results are presented in Table 5.
The obtained results show that for 7 targets plus the catalytic domain of UCE-1, the predicted complexes are unstable (ΔG > 0). Interestingly, there are 12 targets yielding negative estimations of the free energy of binding that can be interpreted as feasible 11-receptor complexes. Among them, the PGFS, TOP2, LGL, MetAP2, SIR2, PPIase, and the catalytic site of UCE-2 present ΔG values higher than −3 kcal/mole, indicating low stability of the complexes predicted for them. On the other hand, better free energies of binding are predicted for AKR, CRK1, MPK13, GSK3B, and the allosteric sites of UCE-1 and UCE-2, being the ubiquitin-conjugating enzymes E2 and those with the lowest predicted ΔG. To get more insights into the possible mechanism of action of 11, the predicted complexes with the later listed targets were analyzed more in detail. CRK1, GSK3B, and MPK13 are protein kinases sharing the same overall folding and containing one ATP binding site. The predicted binding poses of 11 to these protein kinase receptors as well as the network of interactions that it forms with the receptors are presented in Figure 8. Figures representing molecular structures were prepared with UCSF Chimera [32], and interaction networks were produced with Cytoscape [33].
All in all, protein kinase receptor 11 occupies the ATP binding sites. In CRK1 and GSK3, the ligand orients its trimethoxyphenyl moiety toward the bottom of the binding pockets, while the furan substituent occupies the entrance of the binding site, orienting toward the solvent. Both complexes are mainly stabilized by hydrophobic and Van der Waals interactions. The predicted complex with GSK3 shows the possibility of hydrogen bonds between the ligand and the side chains of R110 and C170 through the compound’s carbonyl group and one of the methoxy substituents of the phenyl ring, respectively. Similarly, a hydrogen bond is predicted between one of the methoxy groups of the ligand and K34 of CRK1 in 97% of the MD snapshots used for MM-PBSA calculations.
In the case of MPK13, 11 is predicted to bind in an orientation rotated 180° relative to that observed in GSK3 and CRK1. That is, the furan substituent occupies the bottom of the binding site, while the trimethoxyphenyl moiety is located at the mouth of the pocket and exposed to the solvent. As for GSK3 and CRK1, the ligand is mainly stabilized through hydrophobic and Van der Waals interactions. The predicted complex is further stabilized by hydrogen bonds between the carboxyl group of the compound and the side chains of K33 and D146 and between one of the methoxy substituents and the side chain of N88. The slightly lower free energy of binding estimated for MPK13 relative to GSK3 and CRK1 could be explained by the larger frequency of hydrogen bond formation observed in the first.
The predicted binding poses of 11 to the predicted non-protein kinase receptors are presented in Figure 9. The predicted complex of 11 with AKR shows that the ligand does not directly interact with the NADP cofactor. Instead, it binds at the tunnel that provides access to the cofactor, potentially preventing the enzyme’s substrates access to NADP and their reduction. In this complex, the ligand positions itself favorably to form hydrogen bonds with either the side chain of R242 or the backbone of I325 through its carboxyl oxygen. Additional stability is provided to this complex by the stacking of the furan ring and the trimethoxyphenyl substituent in front of Y39 and Y147, respectively.
According to the predicted free energies of binding, the most stable binding of 11 is obtained at the allosteric sites of UCE-1 and UCE-2. Although these two ubiquitin-conjugating enzymes E2 share only 34% identity, the predicted binding poses of the chemical to them are highly similar. In both cases, the furan group positions itself toward the interior of the cavity, while the trimethoxyphenyl moiety is located at the entrance of the binding pocket. The furan ring is located in both receptors in regions populated with aromatic residues, which favors its stacking with them. Furthermore, the presence of the bulkier part of the ligand (the trimethoxyphenyl group) at the mouth of the site gives stability to the distorted conformations of the enzymes required for the inhibitory effect on ubiquitination previously observed in one human homolog protein [31].
There are several residues at these receptors that could hydrogen bond the ligand; however, this type of interaction is only observed between the ligand’s carboxyl oxygen and either the backbone of V48 or the side chain of Q115 of UCE-1. In contrast, no hydrogen bond is predicted between the ligand and UCE-2. The lack of hydrogen bonds in the UCE-2 complex, despite the highly similar predicted binding modes of 11 to both targets, can be explained by the V48P and Q115L mutations present in UCE-2. This result agrees with the only hydrogen bond observed in the X-ray complex of an allosteric inhibitor with the human Cdc34 ubiquitin-conjugating enzyme E2 [31]. As previously reported in the aforementioned study, the complexes herein reported with UCE-1 and UCE-2 are mainly stabilized by hydrophobic and Van der Waals interactions.
Considering all the obtained results, we postulate that 11 might produce its trypanocidal effect through a multi-target mechanism interfering with different cellular processes. This hypothesis is further supported by the fact that one of the more probable targets of the compound predicted in this research has been shown to be effective for trypanosome inhibition. For example, CRK1, MPK13, and GSK3B are protein kinases that share less than 30% identity between them and play different roles in T. cruzi. GSK3B is essential for the life cycle of T. brucei, and its inhibition leads to a trypanocidal effect [34]. Specifically, indirubin derivatives have demonstrated both GSK3B inhibition and antiparasitic activity against T. brucei and T. cruzi [35,36]. The second protein kinase identified as a potential target of 11 is CRK1. This protein is involved in the cell cycle regulation of trypanosomes, specifically in the G1/S transition, and in protein transport [37,38,39]. The last protein kinase identified as a potential target of the studied compound is MPK13, which is a member of the mitogen-activated protein kinase family. Although no specific role has been reported for MPK13 in T. cruzi, several studies point to a central role for this family of proteins in the regulation and development of the parasite [40,41,42,43].
Among the other probable targets of 11, AKR is annotated as a putative aldo-keto reductase protein in the UniProt database. This superfamily of proteins is related to the maintenance of the redox balance in the cell through the reduction of a variety of substrates [44]. Finally, the ubiquitin conjugating enzymes E2 have been demonstrated to be essential for trypanosome survival [45]. The results obtained by Rojas et al. in [45] also show that the ubiquitin conjugating enzyme E2 from T. brucei is inhibited by an allosteric human enzyme inhibitor. This last result supports the use of the analog allosteric binding site in T. cruzi for UCE-1 and UCE-2 inhibition.
Previous studies on the trypanocidal activity of piplartine have shown that it induces the expression of proteins related to the control of oxidative stress in T. cruzi [16]. On the other hand, several studies demonstrate that interference with the ubiquitination pathway in human cells leads to oxidative stress [46,47]. More importantly, it was reported that the antitumoral activity of piplartine is due to its interference with the ubiquitination pathway, which causes oxidative stress in tumoral cells [48,49]. These precedents and our modeling results strongly point to the ubiquitin-conjugating enzymes E2 of T. cruzi here studied as the preferred targets of the assayed piplartine derivatives.
The presented computational models can guide new experiments to fully clarify the trypanocidal mode of action of this series of chemical compounds and optimize their bioactivity. Our results allow for narrowing the space of molecular targets to be experimentally evaluated as well as prioritizing the proteins with the highest probabilities of serving as targets for these compounds. Future experimental efforts will be directed toward the evaluation of the predicted inhibitory activity of compound 11 against its most probable targets. These assays will include enzymatic inhibition experiments. Afterwards, site-directed mutagenesis experiments can be designed to validate the predicted binding mode of the compound to its confirmed targets. The results of such experiments will enable new rounds of structure-based drug design focusing on the improvement of the potency of piplartine derivatives as trypanocidal agents.

3. Materials and Methods

3.1. Chemical Characterization of Compounds

The infrared was performed by FTIR spectrophotometry, Cary 630 FTIR Agilent Technologies, potassium bromide pills, and frequency measurement in cm−1. The 1H and 13C NMR spectra were obtained on a Mercury-Varian machine operating at 400 MHz and 100 MHz and a Varian-NMR-System machine operating at 500 and 125 MHz. CDCl3 was used as the deuterated solvent. High-resolution mass spectra were performed on Ultraflex II TOF/TOF equipment with a high-performance solid-state laser (λ = 355 nm) and a reflector, using the MALDI technique. Samples were loaded on a steel plate (MTP 384 steel base; Bruker Daltonics Gmbh, Bremen, Germany). Finally, for the analysis of the melting point of the compounds, the heating plate of a device from Microquímica Equipamentos LTDA, model MQAPF 302, was used.

3.2. Preparation of 14

In a 50 mL flat-bottomed flask, 3,4,5-trimethoxycinnamic acid (0.1 g; 0.42 mmol) and 20 mL of the respective alcohol (methanol, ethanol, n-propanol, or isopropyl alcohol) were added. Then, 0.2 mL of concentrated H2SO4 was slowly added. Thus, the reaction was carried out under reflux and magnetic stirring for 6 to 24 h. All reaction steps were monitored by thin-layer chromatography (TLC). After completion of the reaction, the alcohol was partially removed under reduced pressure, then the product was transferred to a separatory funnel with 10 mL of distilled water and extracted using 3 × 10 mL of ethyl acetate. The organic phase was treated with 5% sodium bicarbonate (3 × 10 mL), washed with 10 mL of distilled water, and dried with anhydrous sodium sulfate. The solvent was evaporated under reduced pressure. For the purification of 4, it was necessary to perform adsorption chromatography on a silica gel 60 column, using as the mobile phase a mixture of hexane and ethyl acetate with an increasing polarity gradient [19].

3.3. Preparation of 57

In a 50 mL flat-bottomed flask, 0.1 g (0.42 mmol) of 3,4,5-trimethoxycinnamic acid, 0.43 mmol of the corresponding halide (1-bromopentane, 1-bromodecane, or 4-methoxybenzyl chloride), 5 mL of anhydrous acetone, and 1.68 mmol of triethylamine were added. The reaction was carried out under reflux and magnetic stirring for 24 h and monitored using TLC. The acetone was removed under reduced pressure, then the reaction product was transferred to a separatory funnel and extracted with 10 mL of distilled water and dichloromethane (3 × 10 mL). The organic phase was treated with 5% sodium bicarbonate (3 × 10 mL), 10 mL of distilled water, dried with anhydrous sodium sulfate, and filtered. After extraction, the solvent was evaporated under reduced pressure, and the products were purified by column adsorption chromatography on silica gel 60 using a mixture of hexane and ethyl acetate with an increasing polarity gradient as the mobile phase [11,44].

3.4. Preparation of 813

The 3,4,5-trimethoxycinnamic acid (0.1 g; 0.42 mmol), 0.54 mmol of the respective alcohols (3-methoxybenzyl alcohol, 4-methylbenzyl alcohol, benzyl alcohol, furfuryl alcohol, (−)-perillyl alcohol, or (−)-borneol), 0.015 g of 4-dimethylaminopyridine (DMAP), and 6.3 mL of dichloromethane were added to a 50 mL flat bottom flask. After 5 min 0.11 g (0.54 mmol) of N,N′-dicyclohexylcarbodiimide (DCC) was added, and the reaction was carried out under magnetic stirring and at room temperature for 24–48 h. The solvent was partially removed under reduced pressure, and the reaction product was filtered three times to remove the 1,3-dicyclohexyl urea (DCU). The filtrate was transferred to a separatory funnel and extracted with 10 mL of distilled water and dichloromethane (3 × 10 mL). The organic phase was treated with 10 mL of 5% sodium bicarbonate, 10 mL of distilled water, and dried with anhydrous sodium sulfate. The reaction product was purified on a silica gel 60 column using a mixture of hexane and ethyl acetate [45,46,50].

3.5. In Vitro Cytotoxicity Assay in LLC-MK2 Cells

For the cytotoxicity analysis, the (3-[4,5-dimethylthiazol-2-yl]-2,5 diphenyl tetrazolium bromide) (MTT) assay was used, in which the LLC-MK2 cells were cultured in Dulbecco’s Modified Eagle Medium (DMEM), then exposed to the piplartine analogues at different concentrations (25–200 μg/mL) and incubated for 24 h at 37 °C. Then MTT was added, and the cells were incubated for 4 h, after which 10% sodium dodecyl sulfate was added. Cell viability measurements were performed in a microplate reader (Biochrom® Asys Expert Plus, Cambridge, UK) at a wavelength of 570 nm, and the CC50 was defined. To calculate the relative cell viability, cells treated only with sterile phosphate-buffered saline were used as a negative control. The experiments were conducted in triplicates [51].

3.6. Evaluation of In Vitro Trypanocidal Activity in Y Strains of T. cruzis

3.6.1. Test on Epimastigote Forms of T. cruzi

Initially, epimastigotes were cultured and distributed in 96-well plates at a density of 106 cells/mL in liver infusion tryptose medium supplemented with 10% fetal bovine serum (FBS) and 1% antibiotic. Then, the cells were treated with piplartine analogues at different concentrations and incubated at 28 °C for 24 h. To quantify the growth inhibition of the parasites, a Neubauer chamber was used, and the IC50 was defined [52].

3.6.2. Assay on Trypomastigote Forms of T. cruzi

Trypomastigote forms were obtained by infecting LLC-MK2 cells with the parasite; in this way, the cells were incubated for 24 h in T-25/75 cm2 flasks at 37 °C in an atmosphere with 5% CO2 in DMEM medium supplemented with 2% FBS and 1% antibiotic. Subsequently, the trypomastigote forms were distributed in 96-well plates treated with piplartine analogues at different concentrations (200; 100; 50; 25; 12.5; 6.25; 3.12 μg/mL) and incubated at 37 °C for 24 h. The assessment of cell viability was performed by counting, and the IC50 was defined [53].

3.6.3. The Trypanocidal Mechanism of Action Assays

Thirteen compounds were tested against the epimastigotes form of T. cruzi. The compound that presented a lower IC50 of epimastigotes compared to the other compounds and benznidazole, the standard drug, was chosen to investigate the trypanocidal mechanism of action.

Analysis of Reactive Cytoplasmic Oxygen Species

The DCFH-DA assay was used to evaluate the increase in ROS concentration in the cytoplasm of epimastigote cells treated with 11. Initially, 106 epimastigote cells/mL were added to 24-well plates in the presence of 11 (157 μM and 314 µM). After 3 h, 10 µL of a solution of DCFH-DA (20 µM in DMSO) was added. After 24 h, the cells were centrifuged (2800 RPM for 7 min), washed, and resuspended in 500 µL of phosphate buffered saline for flow cytometer analysis. Thus, the level of labeling in each group was analyzed by the geometric mean of the fluorescence intensity [54], and the results were expressed as relative fluorescence intensity according to the following formula (Equation (1)):
Relative   fluorescence = mTest mControl
mTEST = geometric mean of the interest group; mCONTROL = geometric mean of the control group.

3.6.4. Evaluation of Mitochondrial Transmembrane Potential

Initially, 11 (157 and 314 µM) and epimastigote forms (106 cells/mL) were added in 24-well plates. After 24 h, the cells were centrifuged and washed with phosphate-buffered saline. Subsequently, the cells were resuspended in 100 µL of phosphate-buffered saline, and 5 µL of Rho123 solution (10 µg/mL) was added. Thus, the samples were incubated in the dark for 30 min, washed, and resuspended with phosphate-buffered saline (500 µL/tube) for analysis by flow cytometry [55]. Results were evaluated by relative fluorescence intensity, calculated as described above.

3.6.5. 7-AAD Assay

The epimastigote forms of T. cruzi were cultivated at a density of 106 cells/mL in the presence of 11 (157 and 314 µM) for 24 h in a stove. Then, the cells were transferred to cytometry tubes, centrifuged at 4000 RPM for 5 min and washed twice with phosphate-buffered saline and binding buffer (pH 7.4; a compound of HEPES solution 10 mM; NaCl 140 mM; and CaCl2 2.5 mm). After the last centrifugation, the pellet was resuspended in 100 µL binding buffer, 5 µL 7-AAD (0.5 mg/mL), and 5 µL Ax/PE. The procedure was performed using a commercial kit (Annexin V PE Apoptosis Detection Kit I, BD Biosciences) and following the manufacturer’s guidelines. After 15 min of incubation in the dark, 400 µL of binding buffer was added to each tube, and the cells were analyzed on the FACSCalibur (BD Biosciences). For each tube, a minimum of 104 cells were analyzed to quantify the percentage of unlabeled cells and were either singly or double labeled with 7-AAD and Ax/PE. For this, an argon laser (488 nm) was used to excite the fluorochromes. Ax/PE-labeled cells were read by the FL2 detector (563 to 606 nm, yellow fluorescence), while 7-AAD-labeled cells were read at 615 to 645 nm wavelength (red fluorescence, FL3 detector).

3.6.6. Evaluation of Morphological Changes in T. cruzi Induced by 11

To evaluate the changes induced by 11 in the ultrastructure of T. cruzi, the epimastigote forms were incubated with them for 24 h. After this period, the samples were centrifuged (5000 rpm; 10 min), the supernatant discarded, and the pellet resuspended in glutaraldehyde solution (2.5%) and incubated for 2 h. The fixed samples were again centrifuged and exposed to increasing concentrations of ethanol (30, 50, 70, 90, and 100%) for dehydration, followed by centrifugation (5000 rpm; 5 min). Finally, the cells were transferred to the surface of circular coverslips (15 mm) and dried in a CO2 oven. The coverslips were covered with a gold layer (20 nm thick) using the QT150 ES-Quorum Metallizer and analyzed in a Quanta 450 FEG-FEI scanning electron microscope to observe changes in the three-dimensional structure of the epimastigote forms of T. cruzi.

3.7. Statistical Analysis

For statistical analysis, the GraphPad Prism 5 program (GraphPad Software, San Diego, CA, USA) was used. Values were expressed as mean ± mean standard error (S.P.M), and data were calculated using one-way analysis of variance (ANOVA), followed by Dunnett’s post-test. Significance was defined as * p < 0.05.

3.8. Molecular Docking: Targets Selection

Computational target fishing was performed for 11 using the Similarity Ensemble Approach (SEA) [56]. The identified potential targets were used as input to a NCBI Blast (Altschul, 1997) search against the T. cruzi proteome (taxid: 5693) on the Reference proteins (refseq_protein) database. From the Blast results, T. cruzi proteins sharing at least 35% sequence identity with the SEA predicted potential targets and with alignment coverage higher than or equal to 70% were selected as potential targets of 11. To avoid duplications in the T. cruzi list of potential targets for 11, it was manually curated, and duplicate sequences (e.g., the same protein with different UniProt Accession Codes) were removed.
OpenEye’s Omega [57] and Molcharge [58] software were employed to generate the most stable three-dimensional conformation of 11 and to add atomic partial charges to it, respectively. The X-ray structure of PGFS, PDB code 4GIE, was retrieved from the Protein Data Bank [59]. The rest of the potential targets of the compound had no solved structure, and homology models were generated for them using the SWISS-MODEL server [60]. If necessary, catalytically important cofactors, including metal ions, were added to the homology models, taking their positioning in the models’ templates as a reference.
Molecular docking was performed following the same consensus ranking approach as in our previous publications [15,18]. The Gold [61] software was selected for molecular docking. The receptor binding sites were defined as any residue within 6 Å from the ligands co-crystallized with the homology model templates or with T. cruzi PGFS. The CHEMPLP scoring function was selected for primary docking and was set up to explore 30 different possible binding modes of the ligand to each receptor. The search efficiency parameter of Gold was set to 200%, and the 30 possible binding modes obtained for each receptor were rescored with the GoldScore, ChemScore, and ASP scoring functions of Gold.
The selection of the most probable binding modes of 11 to every receptor was made using the same consensus scoring strategy described in our previous publications [15,18]. For this, the scores predicted for binding pose I according to each scoring function j (Si,j) are aggregated according to Equation (2).
Z i = 1 4 j S i , j S j ¯ s t d S j
where S j ¯ and s t d S j represent the average and standard deviation of scoring function j, respectively, across the 30 predicted binding poses of 11. Either any predicted binding pose with an aggregated score of Zi > 1 or the highest-scored conformation of the ligand were considered for further calculations.

3.9. Molecular Dynamics Simulations and MM-PBSA Calculations

MD simulations and MM-PBSA calculations took place with Amber 18 [62]. All systems were subject to the same simulation procedure that included energy minimization, heating, equilibration, MD production runs, and free energy of binding estimation from MD snapshots. The MD simulations and MM-PBSA calculations protocol were the same as those employed in our previous publication [18].
In brief, the ff14SB force field was used to parametrize proteins, while the gaff force field was employed for non-amino acid molecules. The complexes were embedded in truncated octahedron boxes and solvated with TIP3P water molecules. Excess charges were neutralized by the addition of either Na+ or Cl counterions.
The solvated and neutralized systems were minimized in two stages, the first of which consisted of 500 steps of the steepest descent method followed by 500 cycles of conjugate gradient with all atoms except solvent constrained with a force constant of 500 kcal/mol·Å2. The second stage of energy minimization sequentially implemented 500 steps of the steepest descent algorithm and 1000 cycles of conjugate gradient. No constraint was applied during the second minimization stage. During both energy minimizations, long-range electrostatic interactions were treated with the PME method using a cutoff distance of 10 Å.
Next, the systems were heated from 0 K to 300 K at constant volume along 10,000 steps with a time step of 2 fs. The PME cutoff for heating was set to 12 Å and the solute was constrained with a force constant of 10 kcal/mol·Å2. The temperature was controlled using a Langevin thermostat with a collision frequency of 1.0 ps−1. The SHAKE algorithm was employed to constrain the bonds involving hydrogen atoms, and their interactions were omitted during heating and all subsequent MD steps.
Afterward, the solvated complexes were equilibrated for 100 ps using a time step of 2 fs at a constant pressure of 1 bar and 300 K constant temperature. The temperature was controlled during heating and pressure with isotropic position scaling, setting a relaxation time of 2 ps. The PME cutoff for equilibration was set to 12 Å. Then, 20 different production runs were performed for each equilibrated complex. Different random initial atomic velocities were set for each of these production runs. The parameters for the production runs were the same as those used during the equilibration stage.
MM-PBSA calculations were employed to estimate the free energy of binding 11 to its potential receptors. These were performed with the MMPBSA.py script provided by Amber 18 [63]. For each system, MM-PBSA calculations were performed over 200 MD snapshots (one every 200 ps) selected from the 20 production runs with the ionic strength set to 100 mM.

4. Conclusions

In the present study, it was possible to prepare thirteen esters analogous to piplartine (113), which are structurally related. In the evaluation of trypanocidal activity against T. cruzi, compound 11 showed good activity with IC50 values = 28.21 ± 5.34 μM and 47.02 ± 8.70 μM, against the epimastigote and trypomastigote forms, respectively. In addition, it showed a high rate of selectivity for the parasite (SI > 10). The trypanocidal mechanism of action occurs through the induction of oxidative stress and mitochondrial damage in the parasite’s cells. It was also noted that 11 promotes parasite cell death through necrosis processes and, mainly, late apoptosis. These results could be verified through scanning electron microscopy, which showed that there was a formation of pores and leakage of cytoplasmic content from epimastigote cells. Molecular docking indicated that 11 probably produces a trypanocidal effect through a multi-target mechanism interfering with different cellular processes, including affinity with the protein kinases CRK1, MPK13, and GSK3B, essential for the parasite’s life cycle; the AKR protein involved in the control of cellular oxidative stress; and the ubiquitin E2 conjugating enzymes (UCE-1 and UCE-2), important for the survival of the parasite. Therefore, it is concluded that it was possible to establish chemical characteristics that can serve as a reference for the development of new trypanocidal prototypes with a better biological profile, and other studies can be carried out with piplartine analogous compounds to search for new drugs with trypanocidal activity.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28114512/s1. Chemical characterization compounds 113; Atomic contributions to logP for compounds 113; Predicted logarithm of the partition coefficients for compounds 113. References [11,15,50,64,65,66,67] are cited in supplementary materials.

Author Contributions

Conducting the research, analyzing the data, and writing the original draft preparation were contributed by C.S.M.B.F. Trypanocidal tests were contributed by R.R.P.P.B.d.M., E.P.M. and A.M.C.M. Molecular docking was contributed by Y.P.C. Study planning and coordination were contributed by D.P.d.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Public Call n. 04/2021 Produtividade em Pesquisa PROPESQ/PRPG/Federal University of Paraíba-UFPB, proposal code PIG14859-2021, the National Council for Scientific and Technological Development (CNPq)—grant 306729/2019-9, and the Coordination for the Improvement of Higher Education Personnel (CAPES).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article and in the Supplementary Materials.

Acknowledgments

This work was supported by the Brazilian agencies: the National Council for Scientific and Technological Development (CNPq) and the Coordination for the Improvement of Higher Education Personnel (CAPES). The authors would like to thank the Multi-User Characterization and Analysis Laboratory (LMCA-UFPB) for the analysis of the nuclear magnetic resonance spectroscopic samples and the Northeast Strategic Technologies Center (CETENE-UFPE) and the staff, Julia Campos, for the high-resolution mass spectroscopic analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pérez-Molina, J.A.; Molina, I. Chagas Disease. Lancet 2018, 391, 82–94. [Google Scholar] [CrossRef]
  2. World Heath Organization. Chagas Disease (Also Known as American Trypanosomiasis). Available online: https://www.who.int/news-room/fact-sheets/detail/chagas-disease-(american-trypanosomiasis) (accessed on 14 November 2022).
  3. Pan American Heath Organization. Control, Interrupción de la Transmisión y Eliminación de la Enfermedad de Chagas Como Problema de Salud Pública. Guía de Evaluación, Verificación y Validación; Organización Panamericana de la Salud: Washington, DC, USA, 2019. [Google Scholar] [CrossRef]
  4. Pérez-Molina, J.A.; Crespillo-Andújar, C.; Bosch-Nicolau, P.; Molina, I. Trypanocidal Treatment of Chagas Disease. Enferm. Infecc. Microbiol. Clin. 2021, 39, 458–470. [Google Scholar] [CrossRef] [PubMed]
  5. Arrúa, E.C.; Seremeta, K.P.; Bedogni, G.R.; Okulik, N.B.; Salomon, C.J. Nanocarriers for Effective Delivery of Benznidazole and Nifurtimox in the Treatment of Chagas Disease: A Review. Acta Trop. 2019, 198, 105080. [Google Scholar] [CrossRef]
  6. Müller Kratz, J.; Garcia Bournissen, F.; Forsyth, C.J.; Sosa-Estani, S. Clinical and Pharmacological Profile of Benznidazole for Treatment of Chagas Disease. Expert. Rev. Clin. Pharmacol. 2018, 11, 943–957. [Google Scholar] [CrossRef] [Green Version]
  7. García-Huertas, P.; Cardona-Castro, N. Advances in the Treatment of Chagas Disease: Promising New Drugs, Plants and Targets. Biomed. Pharmacother. 2021, 142, 112020. [Google Scholar] [CrossRef] [PubMed]
  8. Pandey, R.P.; Nascimento, M.S.; Moore, C.E.; Raj, V.S.; Kalil, J.; Cunha-Neto, E. New Approaches for the Treatment of Chagas Disease. Curr. Drug. Targets 2021, 22, 835–841. [Google Scholar] [CrossRef] [PubMed]
  9. Newman, D.J.; Cragg, G.M. Natural Products as Sources of New Drugs over the Nearly Four Decades from 01/1981 to 09/2019. J. Nat. Prod. 2020, 83, 770–803. [Google Scholar] [CrossRef] [Green Version]
  10. Cheng, X.; Tian, P.; Zheng, W.; Yan, X. Piplartine Attenuates the Proliferation of Hepatocellular Carcinoma Cells via Regulating Hsa_circ_100338 Expression. Cancer Med. 2020, 9, 4265–4273. [Google Scholar] [CrossRef] [Green Version]
  11. da Nóbrega, F.; Ozdemir, O.; Nascimento Sousa, S.; Barboza, J.; Turkez, H.; de Sousa, D. Piplartine Analogues and Cytotoxic Evaluation against Glioblastoma. Molecules 2018, 23, 1382. [Google Scholar] [CrossRef] [Green Version]
  12. da Silva, M.A.; Fokoue, H.H.; Fialho, S.N.; Santos, A.P.D.A.D.; Rossi, N.R.D.L.P.; Gouveia, A.d.J.; Ferreira, A.S.; Passarini, G.M.; Garay, A.F.G.; Alfonso, J.J.; et al. Antileishmanial Activity Evaluation of a Natural Amide and Its Synthetic Analogs against Leishmania (V.) Braziliensis: An Integrated Approach in Vitro and in Silico. Parasitol. Res. 2021, 120, 2199–2218. [Google Scholar] [CrossRef]
  13. Mengarda, A.C.; Mendonça, P.S.; Morais, C.S.; Cogo, R.M.; Mazloum, S.F.; Salvadori, M.C.; Teixeira, F.S.; Morais, T.R.; Antar, G.M.; Lago, J.H.G.; et al. Antiparasitic Activity of Piplartine (Piperlongumine) in a Mouse Model of Schistosomiasis. Acta Trop. 2020, 205, 105350. [Google Scholar] [CrossRef]
  14. Ticona, J.C.; Bilbao-Ramos, P.; Flores, N.; Dea-Ayuela, M.A.; Bolás-Fernández, F.; Jiménez, I.A.; Bazzocchi, I.L. (E)-Piplartine Isolated from Piper Pseudoarboreum, a Lead Compound against Leishmaniasis. Foods 2020, 9, 1250. [Google Scholar] [CrossRef] [PubMed]
  15. Turkez, H.; da Nóbrega, F.R.; Ozdemir, O.; Filho, C.d.S.M.B.; de Almeida, R.N.; Tejera, E.; Perez-Castillo, Y.; de Sousa, D.P. NFBTA: A Potent Cytotoxic Agent against Glioblastoma. Molecules 2019, 24, 2411. [Google Scholar] [CrossRef] [Green Version]
  16. Vieira, G.A.L.; da Silva, M.T.A.; Regasini, L.O.; Cotinguiba, F.; Laure, H.J.; Rosa, J.C.; Furlan, M.; Cicarelli, R.M.B. Trypanosoma Cruzi: Analysis of Two Different Strains after Piplartine Treatment. Braz. J. Infect. Dis. 2018, 22, 208–218. [Google Scholar] [CrossRef] [PubMed]
  17. Cotinguiba, F.; Regasini, L.O.; da Silva Bolzani, V.; Debonsi, H.M.; Duó Passerini, G.; Cicarelli, R.M.B.; Kato, M.J.; Furlan, M. Piperamides and Their Derivatives as Potential Anti-Trypanosomal Agents. Med. Chem. Res. 2009, 18, 703–711. [Google Scholar] [CrossRef]
  18. Lopes, S.P.; Castillo, Y.P.; Monteiro, M.L.; de Menezes, R.R.; Almeida, R.N.; Martins, A.; Sousa, D.P. Trypanocidal Mechanism of Action and in Silico Studies of P-Coumaric Acid Derivatives. Int. J. Mol. Sci. 2019, 20, 5916. [Google Scholar] [CrossRef] [Green Version]
  19. Silva, R.H.N.; Andrade, A.C.M.; Nóbrega, D.F.; de Castro, R.D.; Pessôa, H.L.F.; Rani, N.; de Sousa, D.P. Antimicrobial Activity of 4-Chlorocinnamic Acid Derivatives. Biomed. Res. Int. 2019, 2019, 3941242. [Google Scholar] [CrossRef]
  20. Neises, B.; Steglich, W. Simple Method for the Esterification of Carboxylic Acids. Angew. Chem. Int. 1978, 17, 522–524. [Google Scholar] [CrossRef]
  21. Jones, A.; Grkovic, T.; Sykes, M.; Avery, V. Trypanocidal Activity of Marine Natural Products. Mar. Drugs 2013, 11, 4058–4082. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lima, T.C.; Souza, R.J.; Santos, A.D.C.; Moraes, M.H.; Biondo, N.E.; Barison, A.; Steindel, M.; Biavatti, M.W. Evaluation of Leishmanicidal and Trypanocidal Activities of Phenolic Compounds from Calea uniflora Less. Nat. Prod. Res. 2016, 30, 551–557. [Google Scholar] [CrossRef]
  23. Liu, X.; Testa, B.; Fahr, A. Lipophilicity and Its Relationship with Passive Drug Permeation. Pharm. Res. 2011, 28, 962–977. [Google Scholar] [CrossRef] [PubMed]
  24. Hamada, Y. Recent Progress in Prodrug Design Strategies Based on Generally Applicable Modifications. Bioorg. Med. Chem. Lett. 2017, 27, 1627–1632. [Google Scholar] [CrossRef] [PubMed]
  25. Shabbir, A.; Haider, K.; Rehman, K.; Akash, M.S.H.; Chen, S. Biochemical Activation and Functions of Drug-Metabolizing Enzymes. In Biochemistry of Drug Metabolizing Enzymes; Elsevier: Amsterdam, The Netherlands, 2022; pp. 1–27. [Google Scholar] [CrossRef]
  26. Wang, Q.; Xiao, B.-X.; Pan, R.-L.; Liu, X.-M.; Liao, Y.-H.; Feng, L.; Cao, F.-R.; Chang, Q. An LC-MS/MS Method for Simultaneous Determination of Three Polygala Saponin Hydrolysates in Rat Plasma and Its Application to a Pharmacokinetic Study. J. Ethnopharmacol. 2015, 169, 401–406. [Google Scholar] [CrossRef]
  27. Chen, X.; Zhong, Z.; Xu, Z.; Chen, L.; Wang, Y. 2′,7′-Dichlorodihydrofluorescein as a Fluorescent Probe for Reactive Oxygen Species Measurement: Forty Years of Application and Controversy. Free Radic. Res. 2010, 44, 587–604. [Google Scholar] [CrossRef]
  28. Esteras, N.; Adjobo-Hermans, M.J.W.; Abramov, A.Y.; Koopman, W.J.H. Visualization of Mitochondrial Membrane Potential in Mammalian Cells. Methods Cell Biol. 2020, 155, 221–245. [Google Scholar] [CrossRef]
  29. Pokorný, J.; Pokorný, J.; Kobilková, J.; Jandová, A.; Vrba, J.; Vrba, J. Targeting Mitochondria for Cancer Treatment—Two Types of Mitochondrial Dysfunction. Prague Med. Rep. 2014, 115, 104–119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Lama, R.; Sandhu, R.; Zhong, B.; Li, B.; Su, B. Identification of Selective Tubulin Inhibitors as Potential Anti-Trypanosomal Agents. Bioorg. Med. Chem. Lett. 2012, 22, 5508–5516. [Google Scholar] [CrossRef] [Green Version]
  31. Ceccarelli, D.F.; Tang, X.; Pelletier, B.; Orlicky, S.; Xie, W.; Plantevin, V.; Neculai, D.; Chou, Y.-C.; Ogunjimi, A.; Al-Hakim, A.; et al. An Allosteric Inhibitor of the Human Cdc34 Ubiquitin-Conjugating Enzyme. Cell 2011, 145, 1075–1087. [Google Scholar] [CrossRef] [Green Version]
  32. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera?A Visualization System for Exploratory Research and Analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [Green Version]
  33. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  34. Urich, R.; Grimaldi, R.; Luksch, T.; Frearson, J.A.; Brenk, R.; Wyatt, P.G. The Design and Synthesis of Potent and Selective Inhibitors of Trypanosoma brucei Glycogen Synthase Kinase 3 for the Treatment of Human African Trypanosomiasis. J. Med. Chem. 2014, 57, 7536–7549. [Google Scholar] [CrossRef] [PubMed]
  35. Efstathiou, A.; Meira, C.S.; Gaboriaud-Kolar, N.; Bastos, T.M.; Rocha, V.P.C.; Vougogiannopoulou, K.; Skaltsounis, A.-L.; Smirlis, D.; Soares, M.B.P. Indirubin Derivatives Are Potent and Selective Anti-Trypanosoma cruzi Agents. Virulence 2018, 9, 1658–1668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Efstathiou, A.; Gaboriaud-Kolar, N.; Myrianthopoulos, V.; Vougogiannopoulou, K.; Subota, I.; Aicher, S.; Mikros, E.; Bastin, P.; Skaltsounis, A.-L.; Soteriadou, K. Indirubin Analogues Inhibit Trypanosoma brucei Glycogen Synthase Kinase 3 Short and T. brucei Growth. Antimicrob. Agents Chemother. 2019, 63, e02065-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Hammarton, T.C. Cell Cycle Regulation in Trypanosoma brucei. Mol. Biochem. Parasitol. 2007, 153, 1–8. [Google Scholar] [CrossRef]
  38. Hu, H.; Gourguechon, S.; Wang, C.C.; Li, Z. The G1 Cyclin-Dependent Kinase CRK1 in Trypanosoma brucei Regulates Anterograde Protein Transport by Phosphorylating the COPII Subunit Sec31. J. Biol. Chem. 2016, 291, 15527–15539. [Google Scholar] [CrossRef] [Green Version]
  39. Santori, M.I.; Laría, S.; Gómez, E.B.; Espinosa, I.; Galanti, N.; Téllez-Iñón, M.T. Evidence for CRK3 Participation in the Cell Division Cycle of Trypanosoma Cruzi. Mol. Biochem. Parasitol. 2002, 121, 225–232. [Google Scholar] [CrossRef]
  40. Bao, Y.; Weiss, L.M.; Ma, Y.F.; Lisanti, M.P.; Tanowitz, H.B.; Das, B.C.; Zheng, R.; Huang, H. Molecular Cloning and Characterization of Mitogen-Activated Protein Kinase 2 in Trypanosoma cruzi. Cell Cycle 2010, 9, 2960–2968. [Google Scholar] [CrossRef] [Green Version]
  41. Domenicali Pfister, D.; Burkard, G.; Morand, S.; Renggli, C.K.; Roditi, I.; Vassella, E. A Mitogen-Activated Protein Kinase Controls Differentiation of Bloodstream Forms of Trypanosoma brucei. Eukaryot. Cell 2006, 5, 1126–1135. [Google Scholar] [CrossRef] [Green Version]
  42. Kugeratski, F.G.; Batista, M.; de Lima, C.V.P.; Neilson, L.J.; da Cunha, E.S.; de Godoy, L.M.; Zanivan, S.; Krieger, M.A.; Marchini, F.K. Mitogen-Activated Protein Kinase Kinase 5 Regulates Proliferation and Biosynthetic Processes in Procyclic Forms of Trypanosoma brucei. J. Proteome Res. 2018, 17, 108–118. [Google Scholar] [CrossRef] [Green Version]
  43. Müller, I.B.; Domenicali-Pfister, D.; Roditi, I.; Vassella, E. Stage-Specific Requirement of a Mitogen-Activated Protein Kinase by Trypanosoma brucei. Mol. Biol. Cell. 2002, 13, 3787–3799. [Google Scholar] [CrossRef] [Green Version]
  44. Roberts, A.J.; Dunne, J.; Scullion, P.; Norval, S.; Fairlamb, A.H. A Role for Trypanosomatid Aldo-Keto Reductases in Methylglyoxal, Prostaglandin and Isoprostane Metabolism. Biochem. J. 2018, 475, 2593–2610. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Rojas, F.; Koszela, J.; Búa, J.; Llorente, B.; Burchmore, R.; Auer, M.; Mottram, J.C.; Téllez-Iñón, M.T. The Ubiquitin-Conjugating Enzyme CDC34 Is Essential for Cytokinesis in Contrast to Putative Subunits of a SCF Complex in Trypanosoma brucei. PLoS Negl. Trop. Dis. 2017, 11, e0005626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hyun, D.-H.; Lee, M.; Hattori, N.; Kubo, S.-I.; Mizuno, Y.; Halliwell, B.; Jenner, P. Effect of Wild-Type or Mutant Parkin on Oxidative Damage, Nitric Oxide, Antioxidant Defenses, and the Proteasome. J. Biol. Chem. 2002, 277, 28572–28577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Hyun, D.-H.; Gray, D.A.; Halliwell, B.; Jenner, P. Interference with Ubiquitination Causes Oxidative Damage and Increased Protein Nitration: Implications for Neurodegenerative Diseases. J. Neurochem. 2004, 90, 422–430. [Google Scholar] [CrossRef] [PubMed]
  48. Jarvius, M.; Fryknäs, M.; D’Arcy, P.; Sun, C.; Rickardson, L.; Gullbo, J.; Haglund, C.; Nygren, P.; Linder, S.; Larsson, R. Piperlongumine Induces Inhibition of the Ubiquitin–Proteasome System in Cancer Cells. Biochem. Biophys. Res. Commun. 2013, 431, 117–123. [Google Scholar] [CrossRef] [Green Version]
  49. Li, W.; Li, N.; Tang, Y.; Li, B.; Liu, L.; Zhang, X.; Fu, H.; Duan, J. Biological Activity Evaluation and Structure–Activity Relationships Analysis of Ferulic Acid and Caffeic Acid Derivatives for Anticancer. Bioorg. Med. Chem. Lett. 2012, 22, 6085–6088. [Google Scholar] [CrossRef]
  50. Nóbrega, F.R.; Silva, L.v.; da Bezerra Filho, C.S.M.; Lima, T.C.; Castillo, Y.P.; Bezerra, D.P.; Lima, T.K.S.; de Sousa, D.P. Design, Antileishmanial Activity, and QSAR Studies of a Series of Piplartine Analogues. J. Chem. 2019, 2019, 4785756. [Google Scholar] [CrossRef]
  51. Mosmann, T. Rapid Colorimetric Assay for Cellular Growth and Survival: Application to Proliferation and Cytotoxicity Assays. J. Immunol. Methods 1983, 65, 55–63. [Google Scholar] [CrossRef]
  52. da Rodrigues, J.H.S.; Ueda-Nakamura, T.; Corrêa, A.G.; Sangi, D.P.; Nakamura, C.V. A Quinoxaline Derivative as a Potent Chemotherapeutic Agent, Alone or in Combination with Benznidazole, against Trypanosoma Cruzi. PLoS ONE 2014, 9, e85706. [Google Scholar] [CrossRef] [Green Version]
  53. Meira, C.S.; Guimarães, E.T.; dos Santos, J.A.F.; Moreira, D.R.M.; Nogueira, R.C.; Tomassini, T.C.B.; Ribeiro, I.M.; de Souza, C.V.C.; Ribeiro dos Santos, R.; Soares, M.B.P. In Vitro and in Vivo Antiparasitic Activity of Physalis Angulata L. Concentrated Ethanolic Extract against Trypanosoma Cruzi. Phytomedicine 2015, 22, 969–974. [Google Scholar] [CrossRef]
  54. Kessler, R.L.; Soares, M.J.; Probst, C.M.; Krieger, M.A. Trypanosoma Cruzi Response to Sterol Biosynthesis Inhibitors: Morphophysiological Alterations Leading to Cell Death. PLoS ONE 2013, 8, e55497. [Google Scholar] [CrossRef] [Green Version]
  55. Sampaio, T.L.; de Menezes, R.R.P.P.B.; da Costa, M.F.B.; Meneses, G.C.; Arrieta, M.C.V.; Chaves Filho, A.J.M.; de Morais, G.B.; Libório, A.B.; Alves, R.S.; Evangelista, J.S.A.M. Nephroprotective Effects of (−)-α-Bisabolol against Ischemic-Reperfusion Acute Kidney Injury. Phytomedicine 2016, 23, 1843–1852. [Google Scholar] [CrossRef] [PubMed]
  56. Keiser, M.J.; Roth, B.L.; Armbruster, B.N.; Ernsberger, P.; Irwin, J.J.; Shoichet, B.K. Relating Protein Pharmacology by Ligand Chemistry. Nat. Biotechnol. 2007, 25, 197–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hawkins, P.C.D.; Skillman, A.G.; Warren, G.L.; Ellingson, B.A.; Stahl, M.T. Conformer Generation with OMEGA: Algorithm and Validation Using High Quality Structures from the Protein Databank and Cambridge Structural Database. J. Chem. Inf. Model. 2010, 50, 572–584. [Google Scholar] [CrossRef]
  58. OpenEye Scientific Software. QUACPAC; OpenEye Scientific Software: Santa Fe, NM, USA, 2020. [Google Scholar]
  59. Berman, H.M. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [Green Version]
  60. Bienert, S.; Waterhouse, A.; de Beer, T.A.P.; Tauriello, G.; Studer, G.; Bordoli, L.; Schwede, T. The SWISS-MODEL Repository—New Features and Functionality. Nucleic Acids Res. 2017, 45, D313–D319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and Validation of a Genetic Algorithm for Flexible Docking 1 1Edited by F. E. Cohen. J. Mol. Biol. 1997, 267, 727–748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Case, D.A.; Ben-Shalom, I.Y.; Brozell, S.R.; Cerutti, D.S.; Cheatham, T.E., III; Cruzeiro, V.W.D.; Darden, T.A.; Duke, R.E.; Ghoreishi, D.; Gilson, M.K.; et al. AMBER; University of California: San Francisco, CA, USA, 2018. [Google Scholar]
  63. Miller, B.R.; McGee, T.D.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA.Py: An Efficient Program for End-State Free Energy Calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef]
  64. Kumar, S.; Arya, P.; Mukherjee, C.; Singh, B.K.; Singh, N.; Parmar, V.S.; Prasad, A.K.; Ghosh, B. Novel Aromatic Ester from Piper longum and Its Analogues Inhibit Expression of Cell Adhesion Molecules on Endothelial Cells. Biochemistry 2005, 44, 15944–15952. [Google Scholar] [CrossRef]
  65. Jung, J.-C.; Moon, S.; Min, D.; Park, W.K.; Jung, M.; Oh, S. Synthesis and Evaluation of a Series of 3,4,5-Trimethoxycinnamic Acid Derivatives as Potential Antinarcotic Agents. Chem. Biol. Drug Des. 2013, 81, 389–398. [Google Scholar] [CrossRef]
  66. Zhao, Z.; Bai, Y.; Xie, J.; Chen, X.; He, X.; Sun, Y.; Bai, Y.; Zhang, Y.; Wu, S.; Zheng, X. Excavating Precursors from the Traditional Chinese Herb Polygala Tenuifolia and Gastrodia Elata: Synthesis, Anticonvulsant Activity Evaluation of 3,4,5-Trimethoxycinnamic Acid (TMCA) Ester Derivatives. Bioorg. Chem. 2019, 88, 102832. [Google Scholar] [CrossRef] [PubMed]
  67. Katritzky, A.R.; Zhang, S.; Soares, A.; Wang, M. A Facile Synthesis of Benzyl-α, β-Unsaturated Carboxylic Esters. ARKIVOC 2001, 10, 54–57. [Google Scholar] [CrossRef]
Figure 1. Chemical structure of piplartine.
Figure 1. Chemical structure of piplartine.
Molecules 28 04512 g001
Scheme 1. Preparation of esters analogous to piplartine (113). (a) ROH, H2SO4, reflux; (b) alkyl or aryl halide, acetone, Et3N, reflux; (c) ROH, DCC, DMAP, CH2Cl2.
Scheme 1. Preparation of esters analogous to piplartine (113). (a) ROH, H2SO4, reflux; (b) alkyl or aryl halide, acetone, Et3N, reflux; (c) ROH, DCC, DMAP, CH2Cl2.
Molecules 28 04512 sch001
Figure 2. Histogram of staining intensity with DCFH-DA. Experimental groups: control—dark gray; 11 157 μM—light gray; 11 314 µM—pale gray.
Figure 2. Histogram of staining intensity with DCFH-DA. Experimental groups: control—dark gray; 11 157 μM—light gray; 11 314 µM—pale gray.
Molecules 28 04512 g002
Figure 3. Effect of 11 on the production of cytoplasmic reactive species. * p < 0.05 in relation to the control group.
Figure 3. Effect of 11 on the production of cytoplasmic reactive species. * p < 0.05 in relation to the control group.
Molecules 28 04512 g003
Figure 4. Histogram of staining intensity with Rho 123. Experimental groups: control—dark green; 11 157 µM—light green; 11 314 µM—pale green.
Figure 4. Histogram of staining intensity with Rho 123. Experimental groups: control—dark green; 11 157 µM—light green; 11 314 µM—pale green.
Molecules 28 04512 g004
Figure 5. Effect of 11 on mitochondrial transmembrane potential in T. cruzi epimastigotes. * p < 0.05 in relation to the control group.
Figure 5. Effect of 11 on mitochondrial transmembrane potential in T. cruzi epimastigotes. * p < 0.05 in relation to the control group.
Molecules 28 04512 g005
Figure 6. Scanning electron microscopy images of T. cruzi epimastigote forms. (A) Untreated epimastigotes (control), (B) epimastigotes treated with of 11 (157 µM).
Figure 6. Scanning electron microscopy images of T. cruzi epimastigote forms. (A) Untreated epimastigotes (control), (B) epimastigotes treated with of 11 (157 µM).
Molecules 28 04512 g006
Figure 7. Evaluation of the cell death profile of epimastigote cells treated with 11. Data were expressed as mean ± SEM and evaluated by two-way ANOVA with Dunnet’s post-test. * p < 0.05 in relation to the control group.
Figure 7. Evaluation of the cell death profile of epimastigote cells treated with 11. Data were expressed as mean ± SEM and evaluated by two-way ANOVA with Dunnet’s post-test. * p < 0.05 in relation to the control group.
Molecules 28 04512 g007
Figure 8. Predicted binding poses of 11 to the protein kinases CRK1 (top), GSK3 (center), and MPK13 (bottom) and networks of interactions between the ligand and the receptors. The interaction networks were constructed from the 200 MD snapshots used for MM-PBSA calculations. Darker lines indicate more frequent interactions. Receptors are represented in gray and the ligand in cyan, while non-carbon atoms use the following coloring scheme: red for oxygen, blue for nitrogen, and yellow for sulfur. Only residues interacting with the receptors in more than 50% of the analyzed MD snapshots are labeled.
Figure 8. Predicted binding poses of 11 to the protein kinases CRK1 (top), GSK3 (center), and MPK13 (bottom) and networks of interactions between the ligand and the receptors. The interaction networks were constructed from the 200 MD snapshots used for MM-PBSA calculations. Darker lines indicate more frequent interactions. Receptors are represented in gray and the ligand in cyan, while non-carbon atoms use the following coloring scheme: red for oxygen, blue for nitrogen, and yellow for sulfur. Only residues interacting with the receptors in more than 50% of the analyzed MD snapshots are labeled.
Molecules 28 04512 g008
Figure 9. Predicted binding poses of 11 to AKR (top) and the allosteric binding sites of UCE-1 (center) and UCE-2 (bottom). The networks of interactions between the ligand and the receptors are also represented on the left. The interaction networks were constructed from the 200 MD snapshots used for MM-PBSA calculations. Darker lines indicate more frequent interactions. Receptors are represented in gray and the ligand in cyan, while non-carbon atoms use the following coloring scheme: red for oxygen, blue for nitrogen, orange for phosphorus, and yellow for sulfur. Only residues interacting with the receptors in more than 50% of the analyzed MD snapshots are labeled.
Figure 9. Predicted binding poses of 11 to AKR (top) and the allosteric binding sites of UCE-1 (center) and UCE-2 (bottom). The networks of interactions between the ligand and the receptors are also represented on the left. The interaction networks were constructed from the 200 MD snapshots used for MM-PBSA calculations. Darker lines indicate more frequent interactions. Receptors are represented in gray and the ligand in cyan, while non-carbon atoms use the following coloring scheme: red for oxygen, blue for nitrogen, orange for phosphorus, and yellow for sulfur. Only residues interacting with the receptors in more than 50% of the analyzed MD snapshots are labeled.
Molecules 28 04512 g009
Table 1. Trypanocidal activity of compounds 113 against the Y strain of T. cruzi.
Table 1. Trypanocidal activity of compounds 113 against the Y strain of T. cruzi.
EpimastigoteTrypomastigote
EsterSubstituentIC50 (μg/mL)IC50 (μM)IC50 (μg/mL)IC50 (μM)
1Methyl31.5 ± 9.23124.9 ± 36.6--
2Ethyl85.68 ± 19.41321.83 ± 72.89--
3Propyl90.1 ± 29.7321.42 ± 105.95--
4Isopropyl74.5 ± 16.6265.77 ±59.22--
5Pentyl20.2 ± 2.265.50 ±7.13--
6Decyln*n*--
74-methoxybenzyl56.4 ± 13.8157.37 ± 38.50--
83-methoxybenzyl40.27 ± 4.2112.36 ± 11.7264.44 ± 19.01179.80 ± 53.04
94-methylbenzyl61.80 ± 10.52180.49± 30.7339.92 ± 7.34123.06 ± 21.44
10Benzyl39.99 ± 14.03121.66 ± 42.7313.38 ± 4.0640.75 ± 12.36
11Furfuryl8.98 ± 1.7028.21 ± 5.3414.97 ± 2.7747.02 ± 8.70
12(−)-Perillyl70.08 ± 25.56188.16 ± 68.6338.72 ± 5.03103.96 ± 13.51
13
Positive control
(−)-Bornyl
Benznidazole
155.1 ± 43.8
23.0 ± 3.3
414.19 ± 116.97
115.1 ± 16.3
-
32.3 ± 6.4
-
161.4 ± 31.8
n* = could not determine IC50. Benznidazole was used as a positive control.
Table 2. Cytotoxicity of compounds 113 against the LLC-MK2 cells.
Table 2. Cytotoxicity of compounds 113 against the LLC-MK2 cells.
LLCK-MK2Selectivity Index (SI)
EsterSubstituentCC50 (μg/mL)CC50 (μM)EpimastigoteTrypomastigote
1Methyl126.4 ± 25.4501.1 ± 100.74.01-
2Ethyl>200.0>751.1>2.3-
3Propyl147.5 ± 15.6526.2 ± 55.61.6-
4Isopropyl>200.0>713.5>2.7-
5Pentyl>200.0>648.6>9.9-
6Decyl>200.0>528.4n*-
74-methoxybenzyl196.2 ± 26.6547.4 ± 74.23.5-
83-methoxybenzyl100.3 ± 21.7279.9 ± 60.52.51.6
94-methylbenzyl182.8 ± 23.7533.9 ± 69.23.04.6
10Benzyl183.0 ± 47.8557.3 ± 125.04.613.7
11Furfuryl>200.0>628.3>22.2>13.4
12(−)-Perillyl105.9 ± 13.8284.3 ± 37.11.52.7
13(−)-Bornyl71.2 ± 12.7190.1 ± 33.90.5-
n* = could not determine CC50.
Table 3. Potential molecular targets of 11 identified through computational target fishing.
Table 3. Potential molecular targets of 11 identified through computational target fishing.
T. cruzi Target (a)DescriptionID
PGFS_TRYCC9,11-endoperoxide prostaglandin H2 reductase (Prostaglandin F2-alpha synthase)PGFS
Q4D831_TRYCCAldo/keto reductaseAKR
Q4CRI4_TRYCCCell division protein kinase 2 (Cell division protein kinase 2)CRK1
Q4CPK2_TRYCCCell division related protein kinase 2CRK3
Q4DF38_TRYCCCystathionine beta-synthaseCBS
Q4CST7_TRYCCCysteine synthaseCS
Q4DE53_TRYCCDNA topoisomerase 2TOP2
Q4D7B4_TRYCCLactoylglutathione lyase-like proteinLGL
Q4DRC6_TRYCCMethionine aminopeptidase 2 (MAP 2, Peptidase M)MetAP2
Q4D3A0_TRYCCMitogen-activated protein kinaseMPK4 (c)
Q4CNV0_TRYCCNAD-dependent protein deacylase (Regulatory protein SIR2 homolog 5)SIR2
Q4DQ47_TRYCCPeptidylprolyl isomerasePPIase
Q4E4I5_TRYCCProtein kinaseMPK13 (c)
Q4D8B2_TRYCCProtein kinasePK-2
Q4DJM4_TRYCCProtein kinasePK-3
Q4E1X0_TRYCCProtein kinaseGSK3B
Q4CLA1_TRYCCTubulin alpha chainTUB-A (b)
Q4DQP2_TRYCCTubulin beta chainTUB-B (b)
Q4DI08_TRYCCUbiquitin-conjugating enzyme E2UCE2-1
Q4DDU2_TRYCCUbiquitin-conjugating enzyme E2UCE2-2
(a) UniProt entry name. (b) The TUB-A and TUB-B heterodimer is identified as TUB from here on. (c) IDs inferred from homology with Leishmania homolog proteins.
Table 4. Results of the docking of 11 to its potential targets.
Table 4. Results of the docking of 11 to its potential targets.
TargetConformerCHEMPLPGoldScoreChemScoreASPConsensus Z-Score
ScoreZ-ScoreScoreZ-ScoreScoreZ-ScoreScoreZ-Score
PGFS155.380.8719.870.1824.931.9539.551.341.08
256.131.1124.850.2024.171.5839.521.331.05
AKR153.632.0836.752.6925.772.2536.001.252.07
251.841.6230.311.2121.830.8237.501.511.29
353.642.0828.110.7021.560.7336.481.331.21
451.271.4829.070.9223.331.3732.020.551.08
CRK1149.441.4624.231.0219.090.4726.281.931.22
250.321.9026.121.5619.370.6423.260.531.16
349.031.2523.590.8419.941.0025.001.341.11
CRK3157.192.9121.340.8020.571.7931.601.651.79
248.991.1024.531.6219.110.9531.771.701.34
355.012.4325.611.8920.191.5831.031.511.85
CBS150.532.9923.971.4917.832.4331.512.082.25
246.891.6424.581.6515.36-0.0631.292.031.32
TOP2147.311.7232.392.8116.300.5426.992.631.93
LGL155.630.9715.120.4127.542.0336.551.971.35
MetAP2160.841.7830.120.6924.990.7541.402.071.32
261.592.0319.17-1.6629.692.5639.581.371.08
359.041.1432.471.2024.440.5439.121.191.02
MPK4150.543.0022.440.5817.340.6421.510.561.20
SIR2160.772.3130.550.9222.180.9939.461.641.46
257.061.3529.650.7423.361.4936.441.081.17
353.410.4234.401.6822.691.2136.831.151.12
PPIase150.320.8530.103.1722.731.2732.760.281.39
253.582.4023.441.0222.611.1830.56-0.391.05
MPK13171.062.9026.690.7027.891.9733.871.701.82
259.721.1526.640.6926.051.4730.020.811.03
360.421.2628.421.0623.690.8230.460.911.01
PK-2146.082.4122.620.2817.351.4024.651.351.36
238.930.3120.270.2719.702.2925.511.561.11
PK-3146.691.2227.521.7619.561.9529.603.302.06
248.902.1325.681.2918.731.3623.661.111.47
GSK3B150.142.0126.511.2118.851.8628.112.401.87
251.232.3423.630.5618.101.4526.161.481.46
TUB (a)153.551.9129.681.5419.221.7226.380.751.48
CS147.311.7232.392.8116.300.5426.992.631.93
UCE-1 (b)143.012.1923.540.5710.700.3119.280.440.88
UCE-2 (b)151.212.2627.261.3615.181.8025.843.032.11
UCE-1 (Allosteric) (b)146.931.7027.960.2217.770.7326.640.910.89
UCE-2 (Allosteric) (b)160.802.3629.532.2724.611.2231.772.162.00
257.561.4823.810.7225.681.6128.801.261.27
356.571.2127.011.5823.440.7927.790.961.14
(a) Tubulin dimer. (b) Ubiquitin-conjugating enzymes E2 were modeled considering both the catalytic and allosteric binding sites.
Table 5. Predicted free energies of binding of 11 to the studied targets and their components. All energy values are presented in kcal/mole.
Table 5. Predicted free energies of binding of 11 to the studied targets and their components. All energy values are presented in kcal/mole.
TargetConformerMM-PBSA ComponentΔG
of Binding
VDWAALSEELEPBENPOLAREDISPERΔG GasΔG Solv
PGFS1−26.81−8.7019.75−20.4734.45−35.5133.73−1.78
2−23.22−6.7818.11−17.9531.13−30.0031.301.29
AKR1−35.42−16.3830.94−26.8442.96−51.8147.07−4.74
2−35.46−6.2821.80−27.1544.13−41.7438.78−2.96
3−35.89−14.9730.91−27.1844.40−50.8648.13−2.73
4−35.49−14.9731.93−26.9943.63−50.4648.57−1.89
CRK11−39.02−12.9831.65−28.9048.23−52.0050.98−1.02
2−41.25−18.4835.57−30.0049.90−59.7355.46−4.26
3−38.32−12.7232.55−28.0547.09−51.0551.590.55
CRK31−34.41−13.3037.46−25.2145.51−47.7057.7510.05
2−35.85−11.8137.11−26.9246.36−47.6556.558.90
3−30.46−6.2724.34−21.9840.79−36.7343.166.43
CBS1−41.85−14.3139.71−29.6252.03−56.1662.125.96
2−43.60−18.9142.39−29.2452.38−62.5165.533.02
TOP21−52.25−20.6647.10−34.0659.24−72.9172.28−0.63
LGL1−34.52−37.7550.92−29.5448.32−72.2769.70−2.57
MetAP2 1−35.64−7.5930.17−25.3445.31−43.2350.136.90
2−41.01−16.0439.30−29.0951.58−57.0561.794.74
3−41.68−25.6743.35−30.0053.26−67.3566.60−0.75
MPK41−36.97−6.8030.09−26.6046.31−43.7749.796.02
SIR21−36.53−6.0721.95−26.8645.76−42.5940.86−1.73
2−46.74−8.7231.61−32.3055.14−55.4654.45−1.01
3−47.25−7.0929.33−32.6155.42−54.3352.14−2.20
PPIase 1−27.05−4.0715.06−21.3435.67−31.1329.40−1.73
2−26.16−4.9516.80−20.3334.49−31.1130.96−0.15
MPK131−46.26−23.0942.13−33.5855.34−69.3563.90−5.45
2−46.10−20.1238.41−33.0055.24−66.2260.66−5.57
3−45.34−19.6938.34−32.8355.02−65.0360.53−4.50
PK-21−32.10−7.0924.91−23.1040.63−39.1942.443.25
2−33.47−9.7430.99−24.0243.00−43.2149.976.76
PK-31−42.66−13.3437.95−29.4450.88−55.9959.393.40
2−38.93−11.2339.12−28.0749.97−50.1761.0210.85
GSK3B1−39.33−16.5933.76−31.2549.74−55.9352.26−3.67
2−36.30−6.2225.25−26.9845.99−42.5244.261.74
TUB-B1−47.14−6.7831.54−32.2655.49−53.9254.770.85
CS1−38.85−21.0342.10−28.7350.21−59.8863.583.70
UCE-11−30.01−5.4921.02−21.1336.14−35.5136.040.53
UCE-21−32.57−9.2124.57−23.0537.62−41.7839.14−2.64
UCE-1 (Allosteric)1−47.88−3.6924.61−32.9053.44−51.5745.16−6.42
UCE-2 (Allosteric)1−41.24−1.9715.27−28.6845.65−43.2132.25−10.97
2−36.78−3.2116.38−27.0142.67−39.9932.03−7.96
3−40.42−3.1716.67−28.6145.09−43.5933.16−10.43
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Filho, C.S.M.B.; de Menezes, R.R.P.P.B.; Magalhães, E.P.; Castillo, Y.P.; Martins, A.M.C.; de Sousa, D.P. Piplartine-Inspired 3,4,5-Trimethoxycinnamates: Trypanocidal, Mechanism of Action, and In Silico Evaluation. Molecules 2023, 28, 4512. https://doi.org/10.3390/molecules28114512

AMA Style

Filho CSMB, de Menezes RRPPB, Magalhães EP, Castillo YP, Martins AMC, de Sousa DP. Piplartine-Inspired 3,4,5-Trimethoxycinnamates: Trypanocidal, Mechanism of Action, and In Silico Evaluation. Molecules. 2023; 28(11):4512. https://doi.org/10.3390/molecules28114512

Chicago/Turabian Style

Filho, Carlos S. M. B., Ramon R. P. P. B. de Menezes, Emanuel P. Magalhães, Yunierkis P. Castillo, Alice M. C. Martins, and Damião P. de Sousa. 2023. "Piplartine-Inspired 3,4,5-Trimethoxycinnamates: Trypanocidal, Mechanism of Action, and In Silico Evaluation" Molecules 28, no. 11: 4512. https://doi.org/10.3390/molecules28114512

APA Style

Filho, C. S. M. B., de Menezes, R. R. P. P. B., Magalhães, E. P., Castillo, Y. P., Martins, A. M. C., & de Sousa, D. P. (2023). Piplartine-Inspired 3,4,5-Trimethoxycinnamates: Trypanocidal, Mechanism of Action, and In Silico Evaluation. Molecules, 28(11), 4512. https://doi.org/10.3390/molecules28114512

Article Metrics

Back to TopTop