Next Article in Journal
Selectivity Screening and Structure–Cytotoxic Activity Observations of Selected Oleanolic Acid (OA)-Type Saponins from the Amaranthaceae Family on a Wide Panel of Human Cancer Cell Lines
Next Article in Special Issue
Pilose Antler Protein Relieves UVB-Induced HaCaT Cells and Skin Damage
Previous Article in Journal
Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel
Previous Article in Special Issue
N-(2-Hydroxyphenyl)-2-Propylpentanamide (HO-AAVPA) Induces Apoptosis and Cell Cycle Arrest in Breast Cancer Cells, Decreasing GPER Expression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Docking-Based Virtual Screening of FDA-Approved Drugs Using Trypanothione Reductase Identified New Trypanocidal Agents

by
Rogelio Gómez-Escobedo
1,
Domingo Méndez-Álvarez
2,
Citlali Vázquez
3,
Emma Saavedra
3,
Karina Vázquez
4,
Verónica Alcántara-Farfán
5,
Joaquín Cordero-Martínez
5,
Alonzo Gonzalez-Gonzalez
2,
Gildardo Rivera
2,* and
Benjamín Nogueda-Torres
1,*
1
Departamento de Parasitología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Ciudad de México 07738, Mexico
2
Laboratorio de Biotecnología Farmacéutica, Centro de Biotecnología Genómica, Instituto Politécnico Nacional, Reynosa 88710, Mexico
3
Departamento de Bioquímica, Instituto Nacional de Cardiología Ignacio Chávez, Ciudad de México 14080, Mexico
4
Facultad de Medicina Veterinaria y Zootecnia, Universidad Autónoma de Nuevo León, General Escobedo 66050, Mexico
5
Laboratorio de Bioquímica Farmacológica, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Ciudad de México 07738, Mexico
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(16), 3796; https://doi.org/10.3390/molecules29163796
Submission received: 1 July 2024 / Revised: 5 August 2024 / Accepted: 7 August 2024 / Published: 10 August 2024
(This article belongs to the Special Issue The Design, Synthesis, and Biological Activity of New Drug Candidates)

Abstract

:
American trypanosomiasis or Chagas disease, caused by Trypanosoma cruzi (T. cruzi), affects approximately 6–7 million people worldwide. However, its pharmacological treatment causes several uncomfortable side effects, causing patients’ treatment abandonment. Therefore, there is a need for new and better treatments. In this work, the molecular docking of nine hundred twenty-four FDA-approved drugs on three different sites of trypanothione reductase of T. cruzi (TcTR) was carried out to find potential trypanocidal agents. Finally, biological evaluations in vitro and in vivo were conducted with the selected FDA-approved drugs. Digoxin, alendronate, flucytosine, and dihydroergotamine showed better trypanocidal activity than the reference drugs benznidazole and nifurtimox in the in vitro evaluation against the trypomastigotes form. Further, these FDA-approved drugs were able to reduce 20–50% parasitemia in a short time in an in vivo model, although with less efficiency than benznidazole. Therefore, the results suggest a combined therapy of repurposed and canonical drugs against T. cruzi infection.

1. Introduction

Chagas disease, also known as American trypanosomiasis, is an anthropozoonotic disease caused by the Trypanosoma cruzi (T. cruzi) parasite, which is acquired mainly through the feces of hematophagous Hemiptera vectors of the Triatomine subfamily. Other ways of contracting the parasite are congenital transmission, blood transfusions, organ transplants, oral transmission by food contaminated with infected feces, and laboratory accidents. It is estimated that Chagas Disease affects approximately 6–7 million people, mainly in Latin America and Caribbean endemic countries, although cases in non-endemic countries have been documented due to human immigration [1,2]. Chagas disease presents three phases: the acute phase, with a duration of 2–8 weeks with high parasitemia levels but with nonspecific signs and symptoms; the indetermined phase, with no signs or symptoms, where two-thirds of those infected remain in this phase; and the chronic phase, where the remaining third progresses to and where cardiopathies such as arrhythmias, septal defects, cardiac failure, and death occur [3].
At present, there are only two drugs for its treatment, nifurtimox (Nfx) and benznidazole (Bzn), which have shown a moderately parasitological cure in the acute phase but are ineffective in the chronic phase. In addition to the latter, the long periods of treatment and their severe side effects cause people to abandon the treatment, which means that the parasitological cure is not achieved, with the possibility of the emergence of variants that are resistant to these drugs. Therefore, new drugs are required for the treatment of the disease [4].
In the search for new drugs, targeted drug repositioning is a strategy that has been used in the last decades to obtain new trypanocidal agents; the use of this strategy allows for reducing time and costs in the development and search for treatments against existing diseases [5]. Different approaches have been tested to reposition drugs for the treatment of Chagas disease, the most common being through in silico studies. In this approach, molecular docking on a directed target has been a very useful tool that has allowed for the virtual screening of huge libraries of compounds; the strategy has also permitted the establishment of predictions about the potential ways in which drugs/compounds are able to inhibit enzymes and to understand the molecular interactions that exists between them [6].
Trypanothione reductase (TR) is an enzyme that is central in the trypanothione-dependent redox system of the parasite, which indirectly participates in the detoxification of reactive oxygen species (ROS) that the parasite must deal with once it is inside the mammalian host. This redox system involves the utilization of NADPH reducing equivalents that are fed into a variety of enzymatic detoxification pathways utilizing trypanothione as a vehicle. Trypanothione is oxidized as it yields its electrons into the detoxification paths, later it is then reduced by TR, and it may be now recycled [7]. This enzyme represents an attractive pharmacological target since it is not found in the mammalian host and its identity with human glutathione reductase (which is its analogue) is low, so even inhibiting the parasite enzyme, the host enzyme will not be affected. There are three sites of interest: the Z site, mepacrine, and the catalytic site. The Z site is a hydrophobic barrier located in the proximity of the catalytic site, deep within the interface cavity, in close vicinity to the NADPH binding site; in previous studies, it has been mentioned that by binding at this Z site, it is possible to sequester the catalytic site of TR through non-competitive inhibition [8]. The mepacrine site is a hydrophobic region located at the entrance of the trypanothione disulfide binding site; this site is where most inhibitors bind, acting through competitive inhibition. The catalytic site, the place where the reduction of trypanothione disulfide takes place, was chosen to inhibit through competitive inhibition [9]. The catalytic triad is composed of Cys53, Cys58, and His461.
In this study, a molecular docking-based virtual screening of FDA-approved drugs using the enzyme TR of T. cruzi (TcTR) as a pharmacological target was performed centered at each of these three sites to identify new agents with potential trypanocidal activity. Finally, biological in vitro and in vivo evaluations against trypomastigotes from T. cruzi NINOA and INC-5 strains were carried out with the selected FDA-approved drugs.

2. Results

2.1. Binding Sites

Nine hundred twenty-four FDA-approved drugs were evaluated by molecular docking on three distinct sites of the TcTR enzyme: (a) Z site composed of three residues, Phe396, Pro398, and Leu399, is a hydrophobic region present in the neighborhood of the catalytic site and the NADPH binding site [10]; (b) mepacrine site is a hydrophobic region found at the entrance of the trypanothione disulfide binding cavity, the most prominent residues in this site are Trp21 and Met113, and this site takes relevance since most of the inhibitors that have been tested for this enzyme bind here, thus providing competitive inhibition [9]; and (c) catalytic site, where the reduction of trypanothione disulfide takes place, was selected to search for competitive inhibitors, as the residues involved with interactions at this site are mainly His461, Glu466, and Glu467 [9]. The results of the top 20 FDA-approved drugs according to docking score values for each of the docking sites are shown in Table 1. The most promising interaction profiles were constructed from the docking of several known TR inhibitors at each of the docking sites; the full results may be consulted in the Supplementary Materials.

2.2. Molecular Docking on Z Site

Molecular docking on the Z site of the nine hundred twenty-four FDA-approved drugs allowed us to obtain the binding energy value of each drug. The top 20 are shown in Table 1; these compounds have docking scores (DS) ranging from −8.864 to −10.148 kcal/mol, where most are comparable to control inhibitors tested. The control inhibitors ZINC12151998, 7i, and 7e (previously reported as enzyme inhibitors at the Z site) had binding energy values of −10.3, −9.838, and −9.428 kcal/mol, respectively.
Subsequently, an analysis of the ligand–protein interactions of the top 20 FDA-approved drugs and control compounds was performed. The profile with hydrophobic interactions (HI), hydrogen bonds (HB), ππ–cation interactions (π-c), ππ–π stacking (π-s), halogen interactions (HalB), and salt bridges (SB) is shown in Figure 1. As the reported Z site is a triad of amino acid residues deep within the cavity of the interface, all the calculated docking grid boxes encompass the NADPH binding site as well. The NADPH binding site resides in the vicinity of the catalytic site given its function as the source of reducing equivalents; thus, docking in this deep portion of the cavity resulted in a dual site docking. The interacting residues presented in Figure 1 occur mostly within the NADPH cavity, related to residues that are part of the binding of the natural ligand NADPH: Gly197, Gly198, Phe199, Tyr222, Arg223, Arg229, Asn255, and Ile286.

2.3. Molecular Docking on Mepacrine Site

Of the nine hundred twenty-four FDA drugs evaluated by molecular docking on the mepacrine active site, the top 20 FDA-approved drugs with the best binding energy were obtained, with DS values ranging from −8.717 to −11.136 kcal/mol, where most are better than the control compound JWZ. The control compound JWZ had a binding energy of −8.914 kcal/mol (Table 1). The analysis of the ligand–protein interaction profile is shown in Figure 2.

2.4. Molecular Docking on the Catalytic Site

Like the previous binding sites, the top 20 FDA-approved drugs with the best binding energy on the catalytic site were selected with DS values ranging from −9.536 to −11.725 kcal/mol. The control compound ZINC12151998 presented a binding energy of −10.799 kcal/mol (Table 1). The analysis of the ligand–protein interaction profile is shown in Figure 3.

2.5. Molecular Dynamic Analysis

A molecular dynamic analysis was performed for alendronate (at mepacrine site), digoxin (at Z-site), dihydroergotamine (at catalytic site), and flucytosine (at mepacrine site), the four compounds chosen to be further tested as trypanocidal agents to determine the ligand–TcTR complex stability; the dynamics were analyzed with three measurements: RMSD, RMSF, and the radius of gyration.

2.5.1. Root Mean Square Deviation (RMSD) Analysis

Figure 4 shows the RMDS fluctuations for the alendronate–TcTR complex (orange), digoxin–TcTR (blue), dihydroergotamine–TcTR (green), flucytosine–TcTR (violet), and free TcTR (red). The alendronate–TcTR complex shows a maximum fluctuation of 4.31 Å; still, most of the dynamics maintain a fluctuation in the range of 2–3 Å, with a mean of 2.09 ± 0.33 Å, showing stable behavior throughout most of the 120 ns trajectory. The digoxin–TcTR complex shows a maximum fluctuation of 3.72 Å; still, most of the dynamics maintain a fluctuation in the range of 2–3 Å, with a mean of 2.47 ± 0.23 Å, showing stable behavior throughout most of the 120 ns trajectory. The dihydroergotamine–TcTR complex shows a maximum fluctuation of 9.06 Å, which may be observed at about 35 ns; still, the dynamics remain mostly stable after the first 40 Å, with a fluctuation around 4.5 Å, and the overall fluctuation has a mean of 4.66 ± 0.78 Å, while for the plateau-like portion of the graph, it has a fluctuation with a mean of 4.58 ± 0.32 Å. On the other hand, flucytosine–TcTR is the least stable complex, as it shows a fluctuating behavior throughout except for two plateau-like portions at 72–100 ns and 112–116 ns; the overall average RMSD is 5.61 ± 2.40 Å, while for these two short lapses, it holds an average RMSD of 4.39 ± 0.34 Å and 4.36 ± 0.25 Å, respectively.
Ten evenly spaced frames were extracted from each nanosecond trajectory (every 100 ps) for the second half of the trajectory (60–120 ns), where stability was observed for the RMSD plot, and these were analyzed to determine their interaction profiles, which, altogether, were 600 frames. It was observed that for alendronate, the top five most frequently occurring interactions were Ile339 (HI: 95.67%), Val54 (HI: 78.70%), Pro336 (HI: 77.54%), and Val59 (HB: 63.89%), and it may be added that two less frequent, but important interactions occur with His461 (HI: 59.07%) and Glu466 (HI: 57.90%). For the analysis of the interaction profile, it is observed that digoxin maintains a stable binding mode (over 90% observed frequency) with a relatively high number of residues, that is 7, five of these through a strong type of interaction (HB and SB), His428 (SB: 99.67%), Gly198 (HB: 98.34%), Phe199 (HI: 96.01%), Asp232 (HB: 95.02%), Ala285 (HB: 94.52%), Ile286 (HB: 94.35%), and Phe231 (HI: 92.69%); additionally, two are less frequent but are relevant to NADPH binding Arg229 (SB:88.87%) and Gly197 (HB: 76.58%). As for dihydroergotamine, the analysis of the interaction profile over the 60 ns trajectory shows that five residues are the most common, Glu466 (HB: 99.34%), Ile339 (HI: 94.52%), Val59 (HI: 94.02%), Ile107 (HI: 81.56%), and His461 (HB: 58.64%). For the short stable lapse (72–100 ns), flucytosine maintained its profile that was observed to occur with Asn340 (HB: 90.75%), Tyr455 (HB: 59.43%), Pro336 (HI: 58.72%), and Arg472 (HB: 55.52%).

2.5.2. Root Mean Square Fluctuation (RMSF) Analysis

The RMSF fluctuations for ligand–TcTR complexes were determined; the alendronate–TcTR complex (orange), digoxin–TcTR (blue), dihydroergotamine–TcTR (green), flucytosine–TcTR (violet), and free TcTR (red) are shown in Figure 5. RMSF graphs for TR show minor fluctuations located at loop regions, which are prone to fluctuations; still, proteins remain mostly stable through the dynamics, suggesting that ligand interaction did not considerably affect proteins. Altogether, this RMSF supports the notion that the complexes formed with each ligand are stable at the site each binds to and do not significantly alter the residue fluctuations at any regions.

2.5.3. Radius of Gyration Analysis

Figure 6 shows the radius of gyration for ligand–TcTR complexes: the alendronate–TcTR complex (orange), digoxin–TcTR (blue), dihydroergotamine–TcTR (green), flucytosine–TcTR (violet), and apo-TcTR. Complexes remain stable throughout the 120 ns of molecular dynamics with minimal fluctuation; comparing the receptor to the complexes, there is not a major difference, which suggests that the protein remains compact in its dynamics.

2.6. In Vitro Activity on Blood Trypomastigotes

Four drugs (Figure 7) were selected and evaluated at a single concentration to find those with a lysis percentage equal to or higher than the reference drugs on trypomastigotes present in the blood of infected mice (Table 2). Flucytosine has the lowest percentage of lysis on the NINOA strain and alendronate on the INC-5 strain. However, they present similar trypanocidal values to the reference drugs against the NINOA strain; therefore, the half-maximal lytic concentration (LC50) of all drugs was determined in both strains. Alendronate, digoxin, and dihydroergotamine presented LC50 values similar or better than the reference drugs.

2.7. Short-Term In Vivo Assay in a Murine Model of T. cruzi Infection

The results obtained from the evaluation of alendronate, digoxin, and dihydroergotamine mesylate in a short-term in vivo model of infection with the T. cruzi NINOA strain are shown in Figure 8

3. Discussion

In general, in the molecular docking analysis on three active sites of TcTR (Table 1), some drugs appear as potential ligands in two out of the three docked sites, and four drugs are predicted with a high inhibitory potential on the three active sites: digoxin, flucytosine, regorafenib, and vilazodone. These results could be an advantage since in this way, the enzyme could simultaneously be inhibited at different sites, resulting in a better inhibition of the enzyme, leaving the parasite more vulnerable to ROS. However, only two drugs were selected, digoxin and flucytosine, due to reasons that are mentioned specifically for each site.

3.1. Molecular Docking on Z Site

The interaction profiles observed for the binding site of the control compounds were used for comparison with the FDA-approved drugs to test their potential to binding TcTR. The most frequently occurring amino acids interactions were eight (HI: Val195, Phe199, Tyr222, Arg223, and Ile286. HB: Gly197, Arg223, and Arg229). Ten FDA-approved drugs (digoxin, tannic acid, paclitaxel, dibucaine, nafarelin, conivaptan, flucytosine, edrophonium, tolvaptan, and ganirelix) had the best docking score value with a percentage of similarity (≥40%) of the ligand–protein interaction profile, which is presented at the right edge of Figure 1, where the percentage of similarity between the ligands is presented in comparison with the amino acids mentioned above (those that occur most frequently for control compounds in this site); summary interaction profiles may be consulted in the Supplementary Materials. However, some drugs have a molecular weight (>1000 g/mol) and cause several adverse effects or could be potential pan-assay interference compounds (PAINS), which tend to interact nonspecifically with many biological targets, generating false positives; such criteria are described in Supplementary Table S1. Therefore, only digoxin (Figure 4) was proposed for further in vitro studies. Digoxin is a digitalis drug used to treat arrhythmias; there are no reports of having been tested on T. cruzi.

3.2. Molecular Docking on Mepacrine Site

An analysis of the ligand–protein interactions of the top 20 FDA-approved drugs was carried out, and a comparison was made with the set of ligand–protein interaction profiles of the controls used for this site (Figure 2). Five amino acid residues were the most frequently shared (HI: Leu18, Trp22, Ile339, and Leu399 and HB: Tyr 111), Try111, Ile339, and Leu399 being the most recurrent. The FDA-approved drugs with the best docking score value and PS ≥ 40% were posaconazole, digitoxin, flucytosine, alendronate, digoxin, cetrorelix, vilazodone, metocurine, anidulafungine, and plerixafor. However, from the previous criteria applied (Supplementary Table S2), only flucytosine and alendronate (Figure 4) were selected for further studies. Alendronate is used to treat osteoporosis, and it has been previously tested as an inhibitor of T. cruzi farnesyl diphosphate synthase [11,12], whereas flucytosine, an antifungal drug, has also not been tested against T. cruzi.

3.3. Molecular Docking on the Catalytic Site

Out of the top 20 FDA-approved drugs considering the percentage similarity of the ligand–protein interaction profile with the controls was performed, and the amino acid residues with the highest occurrence are six in number (HI: Leu18, Ile107, Ile339, Leu399, and Pro462 and HB: Glu466). The amino acid residues with the highest frequency were Ile339, Phe396, Pro398, Leu399, and Pro462. The drugs with the best binding score value and PS ≥ 40% were flucytosine, digitoxin, telmisartan, nilotinib, zafirlukast, alendronate, digoxin, and dihydroergotamine (Supplementary Table S3), although only dihydroergotamine (Figure 4) was selected for in vitro studies. Dihydroergotamine is used to treat migraines, and it has not been reported to be tested against T. cruzi.

3.4. Molecular Dynamic Analysis

The behavior for the complexes for alendronate and digoxin with TcTR suggests that these have the highest potential to behave as TcTR inhibitors, as throughout the full 120 ns trajectory, these show only minor fluctuations regarding their RMSD, under 2.5 Å. Considering that the interactions for alendronate, as both Val54 and Val59, which have a frequency over 70%, are neighboring to catalytic Cys53 and Cys58 and, additionally, His461, the remaining catalytic triad residue, participates in nearly 60% of the analyzed frames, it suggests that alendronate, while it maintains a significantly stable pose, could potentially block the access of trypanothione disulfide to the catalytic triad. The most frequent interactions for digoxin occur with residues reported to interact with NADPH (Gly197, Gly198, Phe199, Arg229, Ile286); thus, digoxin may act as a molecule that limits the access for the electron carrier and may inhibit the reductase activity. Second, the behavior observed for dihydroergotamine shows a mild potential to act as a TcTR inhibitor, as though it shows instability initially, after the first third (40 ns) of the dynamic trajectory, it attains a relative stability at about 4.5 Å, which is maintained to the end of the trajectory. Dihydroergotamine maintains the interaction with Val59 over 90% of the frames analyzed; as this is a catalytic neighboring residue, it may be suggested that it could potentially hinder the access of trypanothione disulfide to the catalytic residues. Additionally, dihydroergotamine bears an interaction with His461 (a catalytic residue) for nearly 60% of the analyzed frames, altogether suggesting that the pose it attains after fluctuating during the first 40 ns is an orientation that may result in being detrimental to the reductase activity. Flucytosine shows the least stability for the ligand–TcTR complex, setting it as the least probable candidate to act as a TcTR inhibitor; thus, it is likely to act using a different mechanism. For flucytosine, the interacting residues reside at the cavity entrance, thus potentially blocking the access to the trypanothione disulfide, yet the stability is low; thus, its potential to act as a TcTR inhibitor is low. Both the analysis of RMSF and the radius of gyration show that these molecules do not significantly alter either the residue fluctuations at their respective docking sites nor do they destabilize the protein’s 3D conformation, as its radius of gyration remains compact throughout. Altogether, this analysis supports the initial proposal (proposed using molecular docking protocol) of alendronate, digoxin, and dihydroergotamine as potential TcTR inhibitors.

3.5. In Vitro Activity on Blood Trypomastigotes

The results regarding screening (Table 1) show that all the drugs have a lysis percentage comparable to the reference drugs, indicating that their trypanocidal activity is present; for this reason, a dose-dependent trial was continued using different concentrations to determine the LC50. As shown in Table 2 flucytosine had an LC50 above the reference drugs; for this reason, it was discarded to continue its evaluation in vivo, while the remaining three drugs show results that warrant further studies. These results become relevant since blood trypomastigotes is the morphology present in the mammalian host and thus responsible for the acute form of Chagas disease.

3.6. Short-Term In Vivo Assay in a Murine Model of T. cruzi Infection

The control Bzn was the drug with the lowest blood parasite survival after 8 h drug administration (20.83%), followed by alendronate with 59.83%, dihydroergotamine mesylate with 62.44%, and finally digoxin with 71.60% parasite survival. The Kruskal–Wallis statistical test determined that there is a significant difference between the three drugs tested and the positive control, Bzn (p < 0.01); however, when comparing the FDA-approved drugs, we can observe that there is no significant difference between them (p > 0.05). It is observed that the percentage of parasite clearance by the drugs is lower than the control drug; this could be due to the requirements of the metabolism and absorption of the drugs. Notwithstanding, they did present an effect on parasitemia, suggesting that they can be evaluated in combined therapies with Bzn to determine a possible synergistic effect to reduce the treatment time or the dose of the Bzn to achieve a reduction in the adverse effects.

4. Materials and Methods

4.1. Protein Preparation

The three-dimensional structure of TcTR was obtained from the Protein Data Bank (PDB) database with the accession code 1BZL [13]. The protein was prepared in the open access software UCSF-Chimera 1.16 [14], with which all co-crystallized molecules were removed; subsequently, hydrogens and charges were added using the UCSF-Chimera built-in Dock Prep tool and finally converted to PDBQT format by adding the Gasteiger charges with the MGTools 1.5.7 tool. The force field applied was AMBER [15].

4.2. Ligand Library Preparation

The FDA drugs were obtained from the BindingDB database [16] in a single file in .tsv format, then using the Open Babel 3.1 [17] software, all those repeated were eliminated, and with the same tool, the file was converted to .smi format. Then, each drug was separated into individual files and converted to .sdf format; then, the molecules were minimized and converted to .mol2 format, and the force field was MMFF94 [15].

4.3. Molecular Docking

The software used to perform molecular docking of all ligands was AutoDock Vina 1.1.2 [18]. Molecular docking was performed on three distinct active sites: catalytic site, with coordinates X = 24.663, Z = 9.125, and Y = −4.301, obtained through selection of the center of mass of one of the two trypanothione molecules present in 1BZL with PyMOL 3.0.3 [19]; mepacrine site, with coordinates X = 17.733, Z = 13.499, and Y = −1.197; and the site centered at the Z site, with coordinates X = 41.825, Y = 4.351, and Z = −28.343, reported previously (Espinosa-Bustos et al., 2022) [20]. The box size was calculated for each of the different ligands according to the radius of gyration of each ligand (Feinstein & Brylinski, 2015) [21].

4.4. Molecular Docking Analysis

A series of controls (Biscriptolepine, Komaroviquinone, ZINC12151998, 3-Methoxycarpachromene, JWZ, M9Y, M9J, RD0, RD7, M9S, WPE, JV0, WP7, 2JR, WP6, QUM, WP5, WPF, TS8, RDS, diverse quinone derivatives) (Battista et al., 2020; Espinosa-Bustos et al., 2022; Maamri et al., 2021; Matadamas-Martínez et al., 2019; Saha & Sharma, 2015) [9,20,22,23,24] were gathered from different literature reports of TcTR inhibitors and were docked to TcTR at the three docking sites of interest. The docking score threshold values and the most common interaction profiles were obtained from this initial analysis this interaction profile was used as consensus profile to contrast with the profiles for the FDA-approved ligands; full profiles may be consulted in Supplementary Materials. Further on, molecular docking was analyzed using two criteria: binding energy and protein–ligand interaction profile. The binding energy of the ligands was compared with the controls (Supplementary Materials) that were downloaded from the different databases; from this comparison, the drugs that were either above or very close to the binding energy of the controls were selected, and the top ten drugs were selected from each of the sites mentioned above. For the interaction profile, the comparison was made against the controls (Supplementary Materials), selecting those amino acid residues that appear more frequently among the interactions of the drugs and controls with the proteins presented at the right edge of Figure 1–3 column as percentage of similarity (PS); likewise, those FDA drugs that shared these amino acids in greater proportion were selected.

4.5. Molecular Dynamic Analysis

Once the four drugs were chosen, molecular dynamics was carried out in GROMACS software version 2018 using the AMBER03 force field. The topology of each of the chosen drugs was generated with the ACPYPE 2022.1.3 antechamber module with the General Amber Force Field. Solvation was carried out with water molecules in a dodecahedron with the minimum distance from the wall of 10 Å, and the TIP3P water model was used. Next, chlorine and sodium ions were added to neutralize the system with a 50,000-fold energy minimization. Two equilibrium steps were then carried out at 300 K. In the first, the drug was simulated at NVT conditions (constant number of particles, volume, and temperature). For the second step, the drug was simulated at NTP conditions (constant number of particles, pressure, and temperature). Each step reaches duration of 100 ps. The simulation was performed at a temperature of 300 K for a 120 ns trajectory using a V-rescale thermostat and a Parrinello–Rahman barostat with tau_ and tau_p, respectively. The calculation for long-range electrostatic interactions was performed using the Particle Mesh Ewald method, and the LINCS algorithm was used to know the H-bond length restrictions. To determine the stability of each complex, the root mean square deviation (RMSD), the root mean square fluctuation (RMSF), and the radius of gyration (gyr) were obtained (González-González et al., 2023) [25].

4.6. In Vitro Trypanocidal Assay of Blood Trypomastigotes

Parasite obtention was performed using female mice of the CD1 strain weighing 25–35 g, which were infected with blood trypomastigotes intraperitoneally, one group with the NINOA strain and the other with the INC-5 strain of T. cruzi. Once the maximum peak of parasitemia was reached (approximately 2–4 weeks), infected blood was extracted by cardiac puncture using sodium heparin as anticoagulant, and this blood was adjusted to 1 × 106 trypomastigotes/mL (diluted with saline) (Chacón-Vargas et al., 2017) [26]. First, a screening was performed to select drugs with trypanocidal activity; in a 96-well microplate, 90 µL of the previously adjusted blood and 10 µL of the selected FDA drugs were added at a final concentration of 12.5 µg/mL. In parallel, the reference drugs, Nfx and Bzn, were used as positive controls at the same concentration (12.5 µg/mL); same volume of the vehicle in which the drugs were dissolved (DMSO 1%) was used as negative control. Each drug and reference were mounted in triplicate. The microplate was incubated at 4 °C for 24 h; then, the trypomastigote count was performed using the Pizzi–Brener method, which consists of taking 5 µL of blood from the wells, putting it on a slide, and covering it with an 18 × 18 mm coverslip; 15 random fields were counted under a brightfield optical microscope at a magnification of 40×. The percentage of lysis by each compound was obtained by comparing them with the live trypomastigotes of the negative control (no drug) (Juarez-Saldivar et al., 2024) [27]. Those compounds whose percentage of lysis was close or higher than the reference drugs (Nfx and Bzn) were selected to continue with the tests and calculation of the half-maximal lytic concentration (LC50). For this, the compounds were tested at five concentrations in serial dilutions (100–6.25 µg/mL); the percentage of lysis was determined, and the computation of the LC50 was determined using the Probit tool. Later, the results were converted to micromolar units [28].

4.7. Short-Term In Vivo Assay

For the in vivo assay, mice of the CD1 strain weighing 30–40 g were inoculated intraperitoneally with blood trypomastigotes at a concentration of 1 × 105 trypomastigotes/mL, with modifications to the methodology reported by Romanha et al. [29] and Wong-Baeza et al. [30], one group with the NINOA strain and another group with the INC-5 strain. After the necessary days to reach the necessary parasitemia (12–18 dpi), they were separated into 5 groups, each consisting of three mice; the first group was administered with only 4% gum Arabic (vehicle medium for the reference drugs and medicines), the second group Bzn, the third group dihydroergotamine, the fourth group digoxin, and the fifth group alendronate. The dose administered was 100 mg/kg. From this moment on, counted as zero hour, 5 μL of blood was taken from each mouse from the caudal vein; the trypomastigotes count was performed as described above, and the procedure was repeated at two, four, six, and eight hours. The percentage of lysis of each compound was calculated from the comparison of the number of blood trypomastigotes counted compared to time zero for each of the groups.

5. Conclusions

This study, based on a molecular docking approach to evaluate a library of FDA-approved drugs at three different sites of the pharmacological target TcTR, allowed us to identify alendronate, digoxin, and dihydroergotamine as potential anti-T. cruzi drugs by presenting similar or better LC50 values than the reference drugs in an in vitro model; although it is observed that these drugs have this activity, it is necessary to demonstrate in future studies through enzymatic assays that it is acting on the TcTR enzyme or if these may act through an alternative mechanism. However, the effect was not pronounced in the short-term in vivo assay in an animal model; although they did show a trypanocidal effect, it did not reach that of the canonical antichagasic, Bzn. This opens the possibility of combined therapies to reduce the adverse effects of Bzn.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29163796/s1, Tables S1–S3, Controls used for TR, and Interaction profiles for docked control TcTR inhibitors.

Author Contributions

Conceptualization, G.R., B.N.-T. and E.S; methodology, R.G.-E., D.M.-Á. V.A.-F., J.C.-M., K.V., A.G.-G. and C.V.; formal analysis, D.M.-Á.; investigation, R.G.-E.; data curation, D.M.-Á. and C.V.; writing—original draft preparation, R.G.-E. and G.R.; writing—review and editing, D.M.-Á., C.V., K.V., A.G.-G., V.A.-F., J.C.-M., B.N.-T. and E.S.; project administration and funding acquisition, G.R., B.N.-T. and E.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Secretaria de Investigación y Posgrado del Instituto Politécnico Nacional, grant number 20240460 and Consejo Nacional de Humanidades, Ciencias y Tecnología (CONACHCyT), grant number 282663 and CBF2023-2024-2740 to E.S.

Institutional Review Board Statement

The animal study protocol was approved by the Ethics Committee of Instituto Politécnico Nacional (protocol code ZOO-006-2020 and date of approval 12/20/2020) for studies involving animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Acknowledgments

We want to express our acknowledgement to Instituto Politecnico Nacional and Instituto Nacional de Cardiología Ignacio Chávez for their support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Organización Panamericana de la Salud Enfermedad de Chagas. Available online: https://www.paho.org/es/temas/enfermedad-chagas#:~:text=La%20enfermedad%20de%20Chagas%20es,deficitarios%2C%20consider%C3%A1ndosela%20una%20enfermedad%20desatendida (accessed on 30 April 2024).
  2. World Health Organization La Enfermedad de Chagas (Tripanosomiasis Americana). Available online: https://www.who.int/es/news-room/fact-sheets/detail/chagas-disease-(american-trypanosomiasis) (accessed on 30 April 2024).
  3. Rojo-Medina, J.; Ruiz-Matus, C.; Salazar-Schettino, P.M.; González-Roldán, J.F. Enfermedad de Chagas En México. Gac. Med. Mex. 2018, 154, 605–612. [Google Scholar] [CrossRef] [PubMed]
  4. 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] [PubMed]
  5. Adasme, M.F.; Bolz, S.N.; Adelmann, L.; Salentin, S.; Haupt, V.J.; Moreno-Rodríguez, A.; Nogueda-Torres, B.; Castillo-Campos, V.; Yepez-Mulia, L.; De Fuentes-Vicente, J.A.; et al. Repositioned Drugs for Chagas Disease Unveiled via Structure-Based Drug Repositioning. Int. J. Mol. Sci. 2020, 21, 8809. [Google Scholar] [CrossRef] [PubMed]
  6. Juárez-Saldivar, A.; Schroeder, M.; Salentin, S.; Haupt, V.J.; Saavedra, E.; Vázquez, C.; Reyes-Espinosa, F.; Herrera-Mayorga, V.; Villalobos-Rocha, J.C.; García-Pérez, C.A.; et al. Computational Drug Repositioning for Chagas Disease Using Protein-Ligand Interaction Profiling. Int. J. Mol. Sci. 2020, 21, 4270. [Google Scholar] [CrossRef] [PubMed]
  7. Olin-Sandoval, V.; Moreno-Sanchez, R.; Saavedra, E. Targeting Trypanothione Metabolism in Trypanosomatid Human Parasites. Curr. Drug Targets 2010, 11, 1614–1630. [Google Scholar] [CrossRef] [PubMed]
  8. Khan, M.O.F.; Austin, S.E.; Chan, C.; Yin, H.; Marks, D.; Vaghjiani, S.N.; Kendrick, H.; Yardley, V.; Croft, S.L.; Douglas, K.T. Use of an Additional Hydrophobic Binding Site, the Z Site, in the Rational Drug Design of a New Class of Stronger Trypanothione Reductase Inhibitor, Quaternary Alkylammonium Phenothiazines. J. Med. Chem. 2000, 43, 3148–3156. [Google Scholar] [CrossRef] [PubMed]
  9. Battista, T.; Colotti, G.; Ilari, A.; Fiorillo, A. Targeting Trypanothione Reductase, a Key Enzyme in the Redox Trypanosomatid Metabolism, to Develop New Drugs against Leishmaniasis and Trypanosomiases. Molecules 2020, 25, 1924. [Google Scholar] [CrossRef]
  10. Venkatesan, S.K.; Shukla, A.K.; Dubey, V.K. Molecular Docking Studies of Selected Tricyclic and Quinone Derivatives on Trypanothione Reductase of Leishmania Infantum. J. Comput. Chem. 2010, 31, 2463–2475. [Google Scholar] [CrossRef] [PubMed]
  11. Demoro, B.; Caruso, F.; Rossi, M.; Benítez, D.; González, M.; Cerecetto, H.; Galizzi, M.; Malayil, L.; Docampo, R.; Faccio, R.; et al. Bisphosphonate Metal Complexes as Selective Inhibitors of Trypanosoma Cruzi Farnesyl Diphosphate Synthase. Dalton Trans. 2012, 41, 6468. [Google Scholar] [CrossRef]
  12. Assíria Fontes Martins, T.; de Figueiredo Diniz, L.; Mazzeti, A.L.; da Silva do Nascimento, Á.F.; Caldas, S.; Caldas, I.S.; de Andrade, I.M.; Ribeiro, I.; Bahia, M.T. Benznidazole/Itraconazole Combination Treatment Enhances Anti-Trypanosoma Cruzi Activity in Experimental Chagas Disease. PLoS ONE 2015, 10, e0128707. [Google Scholar] [CrossRef]
  13. Bond, C.S.; Zhang, Y.; Berriman, M.; Cunningham, M.L.; Fairlamb, A.H.; Hunter, W.N. Crystal Structure of Trypanosoma Cruzi Trypanothione Reductase in Complex with Trypanothione, and the Structure-Based Discovery of New Natural Product Inhibitors. Structure 1999, 7, 81–89. [Google Scholar] [CrossRef]
  14. 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] [PubMed]
  15. Méndez-Álvarez, D.; Torres-Rojas, M.F.; Lara-Ramirez, E.E.; Marchat, L.A.; Rivera, G. Ligand-Based Virtual Screening, Molecular Docking, and Molecular Dynamic Simulations of New β-Estrogen Receptor Activators with Potential for Pharmacological Obesity Treatment. Molecules 2023, 28, 4389. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, T.; Lin, Y.; Wen, X.; Jorissen, R.N.; Gilson, M.K. BindingDB: A Web-Accessible Database of Experimentally Determined Protein-Ligand Binding Affinities. Nucleic Acids Res. 2007, 35, D198–D201. [Google Scholar] [CrossRef] [PubMed]
  17. O’boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An Open Chemical Toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef] [PubMed]
  18. Trott, O.; Olson, A.J. AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed]
  19. DeLano, W.L. The PyMOL Molecular Graphics System. 2002. Available online: https://cir.nii.ac.jp/crid/1570572701247244160 (accessed on 20 March 2024).
  20. Espinosa-Bustos, C.; Ortiz Pérez, M.; Gonzalez-Gonzalez, A.; Zarate, A.M.; Rivera, G.; Belmont-Díaz, J.A.; Saavedra, E.; Cuellar, M.A.; Vázquez, K.; Salas, C.O. New Amino Naphthoquinone Derivatives as Anti-Trypanosoma Cruzi Agents Targeting Trypanothione Reductase. Pharmaceutics 2022, 14, 1121. [Google Scholar] [CrossRef] [PubMed]
  21. Feinstein, W.P.; Brylinski, M. Calculating an Optimal Box Size for Ligand Docking and Virtual Screening against Experimental and Predicted Binding Pockets. J. Cheminform. 2015, 7, 18. [Google Scholar] [CrossRef]
  22. Maamri, S.; Benarous, K.; Yousfi, M. Identification of 3-Methoxycarpachromene and Masticadienonic Acid as New Target Inhibitors against Trypanothione Reductase from Leishmania Infantum Using Molecular Docking and ADMET Prediction. Molecules 2021, 26, 3335. [Google Scholar] [CrossRef]
  23. Matadamas-Martínez, F.; Hernández-Campos, A.; Téllez-Valencia, A.; Vázquez-Raygoza, A.; Comparán-Alarcón, S.; Yépez-Mulia, L.; Castillo, R. Leishmania Mexicana Trypanothione Reductase Inhibitors: Computational and Biological Studies. Molecules 2019, 24, 3216. [Google Scholar] [CrossRef]
  24. Saha, D.; Sharma, A. Docking-Based Screening of Natural Product Database in Quest for Dual Site Inhibitors of Trypanosoma Cruzi Trypanothione Reductase (TcTR). Med. Chem. Res. 2015, 24, 316–333. [Google Scholar] [CrossRef]
  25. González-González, A.; Méndez-Álvarez, D.; Vázquez-Jiménez, L.K.; Delgado-Maldonado, T.; Ortiz-Pérez, E.; Paz-González, A.D.; Bandyopadhyay, D.; Rivera, G. Molecular Docking and Dynamic Simulations of Quinoxaline 1,4-Di-N-Oxide as Inhibitors for Targets from Trypanosoma Cruzi, Trichomonas Vaginalis, and Fasciola Hepatica. J. Mol. Model. 2023, 29, 180. [Google Scholar] [CrossRef] [PubMed]
  26. Chacón-Vargas, K.F.; Nogueda-Torres, B.; Sánchez-Torres, L.E.; Suarez-Contreras, E.; Villalobos-Rocha, J.C.; Torres-Martinez, Y.; Lara-Ramirez, E.E.; Fiorani, G.; Krauth-Siegel, R.L.; Bolognesi, M.L.; et al. Trypanocidal Activity of Quinoxaline 1,4 Di-N-Oxide Derivatives as Trypanothione Reductase Inhibitors. Molecules 2017, 22, 220. [Google Scholar] [CrossRef] [PubMed]
  27. Juarez-Saldivar, A.; Gómez-Escobedo, R.; Corral-Ruiz, G.; Chacón-Vargas, K.F.; Horta-Montaño, V.; Sanchez-Torres, L.; Vazquez-Jimenez, L.k.; Nogueda-Torres, B.; Rivera, G. Repositioning FDA-Approved Drug Against Chagas Disease and Cutaneous Leishmaniosis by Structure-Based Virtual Screening. Arch. Med. Res. 2024, 55, 102958. [Google Scholar] [CrossRef]
  28. Becerra, N.A.; Espinosa-Bustos, C.; Vázquez, K.; Rivera, G.; Paulino, M.; Cantero, J.; Nogueda, B.; Chacón-Vargas, F.; Castillo-Velazquez, U.; Rodríguez, A.F.E.; et al. Expanding the Chemical Space of Aryloxy-Naphthoquinones as Potential Anti-Chagasic Agents: Synthesis and Trypanosomicidal Activity. Med. Chem. Res. 2021, 30, 2256–2265. [Google Scholar] [CrossRef]
  29. Romanha, A.J.; Castro, S.L.d.; Soeiro, M. de N.C.; Lannes-Vieira, J.; Ribeiro, I.; Talvani, A.; Bourdin, B.; Blum, B.; Olivieri, B.; Zani, C.; et al. In Vitro and in Vivo Experimental Models for Drug Screening and Development for Chagas Disease. Mem. Inst. Oswaldo Cruz 2010, 105, 233–238. [Google Scholar] [CrossRef]
  30. Wong-Baeza, C.; Nogueda-Torres, B.; Serna, M.; Meza-Toledo, S.; Baeza, I.; Wong, C. Trypanocidal Effect of the Benzyl Ester of N-Propyl Oxamate: A Bi-Potential Prodrug for the Treatment of Experimental Chagas Disease. BMC Pharmacol. Toxicol. 2015, 16, 10. [Google Scholar] [CrossRef]
Figure 1. Interaction profiles for the top 20 FDA-approved drugs docked and centered at the Z site. PS: The percentage of similarity; red boxes highlight the eight amino acid residues that are part of the observed consensus profile in docked inhibitors (full information in Supplementary Materials). Dashed lines separate between interaction types, and the number inside the colored squares indicates the number of interactions for that residue with the ligand.
Figure 1. Interaction profiles for the top 20 FDA-approved drugs docked and centered at the Z site. PS: The percentage of similarity; red boxes highlight the eight amino acid residues that are part of the observed consensus profile in docked inhibitors (full information in Supplementary Materials). Dashed lines separate between interaction types, and the number inside the colored squares indicates the number of interactions for that residue with the ligand.
Molecules 29 03796 g001
Figure 2. Interaction profiles for the top 20 FDA-approved drugs docked and centered at the mepacrine site. PS: The percentage of similarity; red boxes highlight the five amino acid residues that are part of the observed consensus profile in docked inhibitors (full information in Supplementary Materials). Dashed lines separate between interaction types, and the number inside the colored squares indicates the number of interactions for that residue with the ligand.
Figure 2. Interaction profiles for the top 20 FDA-approved drugs docked and centered at the mepacrine site. PS: The percentage of similarity; red boxes highlight the five amino acid residues that are part of the observed consensus profile in docked inhibitors (full information in Supplementary Materials). Dashed lines separate between interaction types, and the number inside the colored squares indicates the number of interactions for that residue with the ligand.
Molecules 29 03796 g002
Figure 3. Interaction profiles for the top 20 FDA-approved drugs docked centered at the catalytic site. PS: Percentage of similarity, red boxes highlight the six amino acid residues that are part of the observed consensus profile in docked inhibitors (full information in Supplementary Materials). Dashed lines separate between interaction types, and the number inside the colored squares indicates the number of interactions for that residue with the ligand.
Figure 3. Interaction profiles for the top 20 FDA-approved drugs docked centered at the catalytic site. PS: Percentage of similarity, red boxes highlight the six amino acid residues that are part of the observed consensus profile in docked inhibitors (full information in Supplementary Materials). Dashed lines separate between interaction types, and the number inside the colored squares indicates the number of interactions for that residue with the ligand.
Molecules 29 03796 g003
Figure 4. RMSD graph for fluctuations over time for TcTR complex fluctuation: alendronate (1.25–4.31 Å), digoxin (1.15–3.72 Å), dihydroergotamine (0.98–9.06 Å), flucytosine (1.08–13.1 Å), and TcTR (0.33–2.22 Å).
Figure 4. RMSD graph for fluctuations over time for TcTR complex fluctuation: alendronate (1.25–4.31 Å), digoxin (1.15–3.72 Å), dihydroergotamine (0.98–9.06 Å), flucytosine (1.08–13.1 Å), and TcTR (0.33–2.22 Å).
Molecules 29 03796 g004
Figure 5. The RMSF graph for fluctuations over time for TcTR complex fluctuation: alendronate, digoxin, dihydroergotamine, flucytosine, and apo-TcTR. Blue spiral (alpha helix), green triangle (beta sheet), in between space (loop).
Figure 5. The RMSF graph for fluctuations over time for TcTR complex fluctuation: alendronate, digoxin, dihydroergotamine, flucytosine, and apo-TcTR. Blue spiral (alpha helix), green triangle (beta sheet), in between space (loop).
Molecules 29 03796 g005
Figure 6. The radius of the gyration graph for molecular dynamics over time for alendronate (orange), digoxin (blue), dihydroergotamine (green), flucytosine (violet), and apo-TcTR (red).
Figure 6. The radius of the gyration graph for molecular dynamics over time for alendronate (orange), digoxin (blue), dihydroergotamine (green), flucytosine (violet), and apo-TcTR (red).
Molecules 29 03796 g006
Figure 7. Four selected FDA-approved drugs for biological evaluation against trypomastigotes of T. cruzi in an in vitro model.
Figure 7. Four selected FDA-approved drugs for biological evaluation against trypomastigotes of T. cruzi in an in vitro model.
Molecules 29 03796 g007
Figure 8. The percentage of parasite survival induced by three FDA drugs and benznidazole treatments in a short-term in vivo model of trypanosomiasis (NINOA strain). MD: Dihydroergotamine, D: digoxin, AS: alendronate. BZN reference drug. The control is infection without treatment.
Figure 8. The percentage of parasite survival induced by three FDA drugs and benznidazole treatments in a short-term in vivo model of trypanosomiasis (NINOA strain). MD: Dihydroergotamine, D: digoxin, AS: alendronate. BZN reference drug. The control is infection without treatment.
Molecules 29 03796 g008
Table 1. The docking score (DS) of the top 20 FDA-approved drugs and reference inhibitors on the three molecular docking sites.
Table 1. The docking score (DS) of the top 20 FDA-approved drugs and reference inhibitors on the three molecular docking sites.
Z SiteMepacrine SiteCatalytic Site
DrugDS (kcal/mol)DrugDS (kcal/mol)DrugDS (kcal/mol)
Digoxin−10.148Dactinomycin−11.136Flucytosine−11.725
Anidulafungine−9.873Irinotecan−10.64Digitoxin−11.012
Tannic acid−9.551Posaconazole−10.381Dactinomycin−10.973
Metocurine−9.238Digitoxin−10.379Irinotecan−10.489
Nilotinib−9.203Tannic acid−10.075Anidulafungin−10.44
Palperidone−9.186Flucytosine−10.034Telmisartan−10.334
Crizotinib−9.169Regorafenib−9.731Nilotinib−10.264
Paclitaxel−9.158Tetracaine−9.58Regorafenib−9.899
Dibucaine−9.085Alendronate−9.437Vilazodone−9.849
Nafarelin−9.051Sorafenib−9.34Dutasteride−9.832
Conivaptan−9.049Telbivudine−9.303Lanreotide−9.802
Ponatinib−9.047Digoxin−9.169Zafirlukast−9.777
Flucytosine−9.031Dutasteride−9.156Lapatinib−9.77
Lomitapide−8.979Temsirolimus−9.082Sorafenib−9.748
Edrophonium−8.969Cetrorelix−9.009Tetracaine−9.718
Itraconazole−8.945Lanreotide−8.936Alendronate−9.714
Vilazodone−8.881Vilazodone−8.817Ponatinib−9.627
Regorafenib−8.878Metocurine−8.801Telbivudine−9.569
Tolvaptan−8.866Anidulafungine−8.758Digoxin−9.55
Ganirelix−8.864Plerixafor−8.717Dihidroergotamine−9.536
ZINC12151998−10.3JWZ−8.914ZINC12151998−10.799
7i−9.838
7e−9.428
Table 2. The percentage of lysis and LC50 values of four FDA drugs on Trypanosoma cruzi strains NINOA and INC-5.
Table 2. The percentage of lysis and LC50 values of four FDA drugs on Trypanosoma cruzi strains NINOA and INC-5.
Drug% Lysis at 12.5 µg/mLLC50 µmol
NINOAINC-5NINOAINC-5
Flucytosine13.4 ± 2.921.9 ± 1.5613 ± 221272 ± 59
Alendronate28.2 ± 617.3 ± 15.3174 ± 11277 ± 13.3
Digoxin36.2 ± 233.6 ± 2.545 ± 2.876 ± 10
Dihydroergotamine29.4 ± 1.525 ± 728.1 ± 3.157 ± 2.6
Nifurtimox28.6 ± 4.525 ± 3.6161 ± 33255 ± 39
Benznidazole32.2 ± 630.3 ± 5.5220 ± 40337 ± 34
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

Gómez-Escobedo, R.; Méndez-Álvarez, D.; Vázquez, C.; Saavedra, E.; Vázquez, K.; Alcántara-Farfán, V.; Cordero-Martínez, J.; Gonzalez-Gonzalez, A.; Rivera, G.; Nogueda-Torres, B. Molecular Docking-Based Virtual Screening of FDA-Approved Drugs Using Trypanothione Reductase Identified New Trypanocidal Agents. Molecules 2024, 29, 3796. https://doi.org/10.3390/molecules29163796

AMA Style

Gómez-Escobedo R, Méndez-Álvarez D, Vázquez C, Saavedra E, Vázquez K, Alcántara-Farfán V, Cordero-Martínez J, Gonzalez-Gonzalez A, Rivera G, Nogueda-Torres B. Molecular Docking-Based Virtual Screening of FDA-Approved Drugs Using Trypanothione Reductase Identified New Trypanocidal Agents. Molecules. 2024; 29(16):3796. https://doi.org/10.3390/molecules29163796

Chicago/Turabian Style

Gómez-Escobedo, Rogelio, Domingo Méndez-Álvarez, Citlali Vázquez, Emma Saavedra, Karina Vázquez, Verónica Alcántara-Farfán, Joaquín Cordero-Martínez, Alonzo Gonzalez-Gonzalez, Gildardo Rivera, and Benjamín Nogueda-Torres. 2024. "Molecular Docking-Based Virtual Screening of FDA-Approved Drugs Using Trypanothione Reductase Identified New Trypanocidal Agents" Molecules 29, no. 16: 3796. https://doi.org/10.3390/molecules29163796

APA Style

Gómez-Escobedo, R., Méndez-Álvarez, D., Vázquez, C., Saavedra, E., Vázquez, K., Alcántara-Farfán, V., Cordero-Martínez, J., Gonzalez-Gonzalez, A., Rivera, G., & Nogueda-Torres, B. (2024). Molecular Docking-Based Virtual Screening of FDA-Approved Drugs Using Trypanothione Reductase Identified New Trypanocidal Agents. Molecules, 29(16), 3796. https://doi.org/10.3390/molecules29163796

Article Metrics

Back to TopTop