Next Article in Journal
Revealing the Therapeutic Targets and Mechanism of Ginsenoside Rg1 for Liver Damage Related to Anti-Oxidative Stress Using Proteomic Analysis
Next Article in Special Issue
Molecular Docking and In-Silico Analysis of Natural Biomolecules against Dengue, Ebola, Zika, SARS-CoV-2 Variants of Concern and Monkeypox Virus
Previous Article in Journal
Irrigation with Magnetized Water Alleviates the Harmful Effect of Saline–Alkaline Stress on Rice Seedlings
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 17
10.3390/ijms231710047
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

by
Lenci K. Vázquez-Jiménez
1,
Alfredo Juárez-Saldivar
1,
Rogelio Gómez-Escobedo
2,
Timoteo Delgado-Maldonado
1,
Domingo Méndez-Álvarez
1,
Isidro Palos
3,
Debasish Bandyopadhyay
4,
Carlos Gaona-Lopez
1,
Eyra Ortiz-Pérez
1,
Benjamín Nogueda-Torres
2,
Esther Ramírez-Moreno
5 and
Gildardo Rivera
1,*
1
Laboratorio de Biotecnología Farmacéutica, Centro de Biotecnología Genómica, Instituto Politécnico Nacional, Reynosa 88710, Mexico
2
Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Ciudad de México 11340, Mexico
3
Unidad Académica Multidisciplinaria Reynosa-Rodhe, Universidad Autónoma de Tamaulipas, Reynosa 88779, Mexico
4
Department of Chemistry and SEEMS, University of Texas Rio Grande Valley, Edinburg, TX 78539, USA
5
Escuela Nacional de Medicina y Homeopatía, Instituto Politécnico Nacional, Ciudad de México 07320, Mexico
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(17), 10047; https://doi.org/10.3390/ijms231710047
Submission received: 30 June 2022 / Revised: 25 August 2022 / Accepted: 30 August 2022 / Published: 2 September 2022
(This article belongs to the Special Issue New Avenues in Molecular Docking for Drug Design 2022)
Browse Figures
Versions Notes

Abstract

:
Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC50) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.

1. Introduction

Chagas disease or American trypanosomiasis is a “neglected tropical” disease because it occurs mainly in low-income countries [1]. This disease causes 56,000 new cases annually in Latin America, and it is estimated that it causes approximately 12,000 deaths per year [2]. Nifurtimox (Nfx) and benznidazole (Bzn) are the only drugs available on the market for the treatment of this disease, which are not suitable because they cause severe adverse effects. Additionally, both drugs have poor efficacy in the chronic phase [3].
Chagas disease is caused by the parasite Trypanosoma cruzi (T. cruzi), a protozoan dependent on glycolysis as the only source of supply of adenosine triphosphate (ATP) for energy generation [4]. Therefore, glycolytic enzymes are potential targets in the search for new therapeutic agents [5]. In this context, triosephosphate isomerase (TIM) enzyme has been proposed by different research groups as a target for the design of new drugs against T. cruzi [6,7].
TIM from T. cruzi (TcTIM) catalyzes the reversible interconversion of the products that come from the catalysis of aldolase (glyceraldehyde-3-phosphate and dihydroxyacetone phosphate) through an intermediate enediol (alkenol), in the fifth step of the glycolysis pathway [8]. TcTIM structurally adopts an α/ß barrel folding, where each TIM monomer is made up of eight parallel ß strands that form the inner part of the barrel and are surrounded by eight α helices that are joined by handles [9]. An important site for the binding of various small molecules is the interface [5,10,11,12,13,14]; it is found among the monomers and occupies a significant part of the molecular surface area of each monomer (1490 Å2) in addition to being mainly composed of loops 1, 2, 3, 4 and 8 [15,16,17].
In the search for new molecules with biological activity, benzimidazole derivatives have presented a variety of activities [18,19,20] highlighting the antiparasitic activity against T. cruzi [21,22]. Among these compounds, two 2,3-dihydroimidazo[1,2-a]benzimidazole analogs have been shown to have trypanocidal effects (IC50 = 1.10 and 2.10 μM), with a better activity than the reference drug Bzn (IC50 = 20.7 μM) and without severe cytotoxicity against U2OS cells (CC50 = 36.5 and 18.8 μM) [22]. Other authors reported that benzimidazole derivatives (IC50 = 5 μM) are more active than Bzn (IC50 = 7.5 μM) to inhibit the growth of the parasite in epimastigote form [23]. Furthermore, the benzimidazole scaffold is considered a privileged structure in medicinal chemistry, as it has the potential to interact with different biological targets, including TcTIM [24,25].
In the present work, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole substructure was carried out; the compounds obtained were analyzed by molecular docking, and subsequently, a selection of compounds was made based on the interaction profile and docking score for an in vitro evaluation against trypomastigotes of T. cruzi and a molecular dynamics simulation analysis. Benzimidazole derivatives have been used in drug designing; however, these kinds of compounds have shown some disadvantages such as hepatotoxicity, among others [26,27,28,29,30]. Therefore, in this study, a prediction of the adsorption, distribution, metabolism and excretion (ADME) properties by computational analysis was performed. Finally, to know their potential selectivity, a molecular docking analysis and molecular dynamics simulation of the compounds against human TIM (HsTIM) were performed.

2. Results

2.1. TcTIM Inhibitors Analysis

In a first step, ten compounds with inhibitory activity against TcTIM were analyzed by molecular docking. Table 1 shows the docking score and interactions profile at the TcTIM interface. Six compounds are benzimidazole derivatives (L1–L6), and four presented a different structure. L7 and 3-(2-benzothiazolylthio-1-propanesulfonic acid (BTS) are benzothiazoles; L8, a thiazole; and L9, a cyclohexanone derivative. The control ligands derived from benzimidazole (L1-L6) show inhibition percentages of 48% to 69%, at a concentration of 200 µM (except the ligand L3) [21,25]. Ligands L7-L9 show IC50 values between 0.086 and 8.0 µM (Table 1) [30,31,32,33].
To determine the most favorable binding site of control ligands for molecular docking analysis, the structure of TcTIM protein (PDB 1SUX) was analyzed with the DoGSiteScorer tool. The results showed two possible binding sites (interface and the active site) with a score (drug score) greater than 0.5. This value indicates that both sites could bind with high affinity to a drug. Additionally, blind molecular docking was performed (PyRx software was used), which showed that most of the control ligands bind at the interface, with docking scores ranging from −6.2 to −8.9 Kcal/mol (Table 1). The results of both analyses showed that the best binding site for control ligands is the interface (Figure 1).

2.2. Molecular Docking of Control Ligands at the TcTIM Interface

Once the possible binding site of the control ligands was identified, a molecular docking analysis was carried out at the TcTIM interface (PyRx software was used). Table 1 shows the docking scores of each of the control ligands evaluated, as well as the interactions with the amino acid residues. The compound with the lowest score (−8.9 Kcal/mol) was ligand L1. The compound BTS co-crystallized presented the highest docking score of −6.2 Kcal/mol.

2.3. LBVS from ZINC15 Database and Molecular Docking Analysis

Considering the benzimidazole scaffold, a LBVS of 750 million compounds from ZINC15 database was carried out. A total of 67,141 molecules per substructure was obtained. Subsequently, the Lipinski’s rule of five was applied as an inclusion criterion, obtaining 53,410 molecules. These compounds were analyzed by molecular docking at the TcTIM interface, obtaining 1604 compounds that showed a binding energy value between −10.6 and −8.9 Kcal/mol. Subsequently, these compounds were grouped based on the type of interactions present in each compound (using the scikit-learn library and the DataWarrior program), obtaining ten groups (Table 2) [32,33]. Group three had the largest number of compounds (206), and group six had the fewest number of compounds (103).
Figure 2 shows the interactions of the lead compound of each group. The interactions that predominated in the ten lead compounds were hydrophobic interactions, hydrogen bonds, and π-stacking. The most common interaction is hydrophobic with amino acid residues Tyr102 (A), Ile109 (B), and Ile69 (B) in all ten compounds, followed by hydrogen bonding interactions with amino acid residues Thr70 (B), Tyr103 (B), and Lys113 (B) in six compounds.

2.4. Trypanocidal Activity

Based on availability and accessibility, compounds BP2 and BP5 (Table 2) were purchased to be evaluated against blood trypomastigotes of the NINOA and INC-5 strains of T. cruzi to determine their trypanocidal activity. The results are shown in Table 3. Compound BP2 presented better half-maximal lytic concentration (LC50) against the NINOA strain (LC50 = 155.86 ± 3.4 µM), a similar value to Bzn, however, two times less than Nfx. Compound BP5 showed similar activity against the two strains evaluated (LC50 = 179.55 ± 19.7 and 179.71 ± 19.0 µM, respectively).
Compounds BP2 and BP5 showed trypanocidal activity and were considered as potential TcTIM inhibitors due to their high docking scores (−10.4 and −10.2 Kcal/mol, respectively); therefore, a molecular dynamic simulation at the TcTIM interface was done to predict the stability of ligand–TcTIM complex.

2.5. Molecular Dynamics Simulation

In the molecular dynamics analysis, the TcTIM (free) showed a Root Mean Square Deviation (RMSD) from 0.33 to 2.04 Å, with a difference in oscillation of 1.71 Å and a mean of 1.43 Å (Figure 3). The L1–TcTIM complex showed stability in the first 45 ns, with an RMSD of 2.17 Å and a mean of 6.17 Å. Overall, BP2–TcTIM complex showed a fluctuation from 0.66 to 12.16 Å, with a difference between the oscillations of 11.5 Å, and a mean of 7.55 Å. Meanwhile, the BP5–TcTIM complex showed similar stability as L1 in the first 45 ns. Additionally, in the next 55 ns, it showed good stability without big fluctuations (RMSD from 0.97 Å to 6.07 Å, showing a difference between the oscillations of 5.11 Å and a mean of 3.59 Å). Additionally, the RMSD of the ligand BTS co-crystallized with the protein was analyzed, which presented an RMSD value of 1.02 to 11.36 Å, and the difference in the oscillation is 10.34 Å, with a mean of 8.65 A.
The RMSD obtained could be influenced by changes that occur in other parts of the structure; therefore, a clustering of the RMSD of each ligand-protein versus the simulation time was performed (Figure 4). An RMSD matrix of each complex was made (Appendix A, Figure A1), and the Ward method and the Euclidean clustering distance were applied to obtain more information on the trajectory of molecular dynamics [34,35]. For the RMSD clustering of BP2, three clusters were observed (corresponding to 59.5, 10.4 and 30.1% frames, respectively), whereas for BP5 and L1, two clusters were observed (for BP5 of 51 and 49 % frames and L1 of 45.1 and 54.8%, respectively) [36]. Additionally, the clustering of the RMSF and Rg was formed (Appendix A, Figure A2, and Figure A3, respectively). In addition to the clustering of the RMSD, an interaction frequency analysis by group was performed to observe the rate of interactions present throughout the analyzed time, which is shown in Appendix A (Figure A4).
The Root Mean Square Fluctuation (RMSF) was also analyzed during molecular dynamics simulation (Figure 5). The RMSF showed a very similar fluctuation pattern between the free protein and the protein in complex with the ligands, especially with the BP5 compound. However, the RMSF calculation revealed that the complex with the BP2 ligand showed a high fluctuation in some regions according to the RMSD pattern as well as in the L1 and BTS ligands in variable proportions.
Additionally, the radius of gyration (Rg) of the free protein was analyzed, as well as of the protein in complex with the ligands (L1, BP2, BP5 and BTS). This parameter allows the calculation of the structural variations of the protein during the analysis of the molecular dynamics [37]. The free TcTIM folding maintains an almost constant fluctuation between 24.64 and 25.57 Å, with a difference in the oscillation of 0.93 Å over 100 ns (Figure 6). The L1–TcTIM complex presented fluctuations between 24.60 and 25.44 Å, with a difference of 0.83 Å. The BP2–TcTIM complex showed a fluctuation from 24.51 to 25.36 Å, with a difference of 0.84 Å, whereas the complex with the BP5 compound fluctuated from 24.72 to 25.50 Å, with a difference of 0.77 Å, this being the one with the smallest difference of the three complexes. On the other hand, the BTS–TcTIM complex showed a fluctuation with a minimum of 24.68 Å and a maximum of 25.40 Å, whereas the difference in oscillation was 0.72 Å.
Additionally, the mean run values of the different components of the calculated Molecular Mechanics Poisson–Boltzmann Surface Area (MMGBSA) binding free energies were determined (Table 4).

2.6. Analysis of Molecular Physicochemical Properties

Molecular physicochemical properties were calculated for the compounds BP2 and BP5 (SwissADME website was used). The descriptors obtained (Table 5) were: physicochemical properties (number of hydrogen bond acceptors, number of hydrogen bond donors, number of rotational bonds, polar surface area, molecular weight, and partition coefficient), solubility coefficient, pharmacokinetics (human gastrointestinal absorption, blood-brain permeability, P-glycoprotein substrate, and inhibition of different CYP450 isoforms), and hepatotoxicity (ProTox-II server was used).

2.7. Molecular Docking at the HsTIM Interface

HsTIM and TcTIM interfaces have been described with 52% sequence identity [38]. Therefore, compounds capable of specifically disrupting the TcTIM interface but not the HsTIM interface would be expected to be converted to selective trypanocidal compounds. Thus, the binding affinity of BP2, BP5, and the control ligand (L1) at the HsTIM interface was analyzed by molecular docking (Figure 7). L1 presented a docking score of −7.2 Kcal/mol, BP2 presented a score of −8.0 Kcal/mol, and BP5 presented a docking score of −5.9 Kcal/mol.

2.8. Molecular Dynamics Simulation at the HsTIM Interface

To complement the results obtained in the molecular docking analysis of L1, BP2, and BP5 at the HsTIM interface, molecular dynamics simulations were carried out (Figure 8). The RMSD values of HsTIM (free) remained constant with a minimum fluctuation of 0.3 Å and a maximum of 3.2 Å, and a difference in oscillation of 2.9 Å and a mean of 1.87 Å (Figure 8A). The RMSD value of L1–HsTIM reached its equilibrium at 3 ns with a minimum fluctuation of 0.3 to a maximum fluctuation of 3.2 Å, with a difference in oscillation of 2.9 Å and a mean of 1.87 Å. Subsequently, the RMSD values for the BP2 and BP5 complexes showed a rise in RMSD values, with a mean fluctuation of 24.69 and 27.91 Å, respectively. Figure 8B shows the RMSF plot, which showed high fluctuations in most of the regions according to the RMSD pattern.
In addition, the Rg was determined (Figure 8C) for the protein alone (24.36 to 25.42 Å, difference is 1.06 Å and a mean of 24.89 Å) as well as for the complexes. For L1 and BP2 complex, the Rg was very similar with differences in the oscillation of 1.36 and 1.30 Å, and with means of 24.95 and 24.99 Å, respectively. BP5–HsTIM presented a minimum fluctuation of 24.26 Å and a maximum of 25.29 Å, with a difference of 1.03 Å and a mean of 24.71 Å.
Finally, as for the Ligand–TcTIM complexes, the mean execution values of the different components of the free binding energies of the MMGBSA of the molecular mechanics calculated for the Ligand–HsTIM complexes were also determined, which are shown in Table 4.

3. Discussion

3.1. TcTIM Inhibitors Analysis

Ten compounds reported in the literature with inhibitory activity against TcTIM were selected as control ligands for molecular docking analysis. Compound L1 showed the lowest score (−8.9 Kcal/mol) [22]; this benzimidazole derivative has a 5-nitro-1,3-thiazol-2-yl aminocarbonyl-2-methylsulfanyl fragment at 2-position, and a 1-naphthyloxy group and a chlorine atom at 5- and 6-position, respectively. L4 (−7.3 Kcal/mol), L5 (−7.2 Kcal/mol), and L6 (−7.2 Kcal/mol) have a similar structure as L1. However, the 1-naphthyloxy group was replaced by a methoxycarbonyl group. This replacement could explain the difference in the docking values between these compounds, according to the interactions calculated. L1 has a larger structure that generates a greater number of hydrophobic interactions. Compounds L2, L3, L7 (1,2,4-thiadiazole derivatives), and L8 (benzothiazole derivative) presented higher docking scores (between −6.8 and −6.9 Kcal/mol). On the other hand, the compound BTS co-crystallized with the protein at the interface binding site [39] presented the highest docking score of −6.2 Kcal/mol. These results suggest that the benzimidazole scaffold is key to obtaining a better binding at the TcTIM interface.
All control ligands showed hydrophobic interactions, hydrogen bonds, and π-stacking interactions. Hydrophobic interactions have been reported to have a greater effect on the docking score because this type of ligand-protein interaction is important in the function of a molecule [40]. All control ligands showed interactions with the amino acid residues Tyr102 (A), Tyr103 (A), and Ile69 (B), suggesting a triad catalytic. Various researchers have described interactions with residues important in the binding of compounds at the TcTIM interface. Other interactions with Asn67, Phe75, Thr70, Glu105, and Lys113 have been reported [31,41,42,43].

3.2. LBVS from ZINC15 Database and Molecular Docking Analysis

The LBVS allowed 67,141 molecules with a benzimidazole scaffold to be obtained from the ZINC15 database. After that, 1604 molecules by molecular docking showed better binding energy than ligand L1. The lead compound of the ten groups obtained from molecular docking analysis showed interactions with amino acid residues previously reported as important for the binding on the TcTIM interface, such as Tyr102 (A), Tyr103 (A), and Ile69 (B) [25,26,29,31]. The presence of hydrogen bonds, as well as hydrophobic interactions, also promote the formation of stronger and more robust ligand–protein complexes [32]. There were also some π-stacking interactions that are also related to binding stability [6], these being with the amino acid residues Tyr103 (A), Tyr102 (A), and Tyr103 (B) in all compounds except compound BP2.
On the other hand, only compounds BP4 and BP6 presented a cation-π interaction with Lys113 (B), whereas compound BP1 additionally had interactions with residues Glu115 (B) and Asn67 (B), which contributed to stability through hydrophobic interactions [43]. It is worth mentioning that the control ligands also presented these types of interactions with most of the residues, where Tyr102 (A), Tyr103 (A), and Ile69 (B) stand out in the ten control ligands.

3.3. Trypanocidal Activity

According to low price and accessibility, compounds BP2 and BP5 were purchased for trypanocidal effect analysis. Both compounds showed trypanocidal activity against the trypomastigote form of T. cruzi, the infective form of the parasite in the mammalian host [44].

3.4. Molecular Dynamics Simulation at the TcTIM Interface

The free protein (TcTIM), L1-TcTIM, B2-TcTIM, B5-TcTIM, and BTS-TcTIM were analyzed by molecular dynamics simulations to determine the stability of the complex formed [44,45]. In addition, the complex of the ligand co-crystallized with the protein (BTS-TcTIM) was analyzed to expand the information that has been described in the literature.
L1-TcTIM showed a fluctuation peak in RMSD values, suggesting a change in the initial binding position, outside of the interface [43]. According to the literature, high RMSD values of a complex compared to the initial framework indicate a change in the initial binding position [33,40]. The results obtained suggest that BP5–TcTIM is better than the control ligand (L1–TcTIM) and the co-crystallized ligand (BTS–TcTIM), since it has been described that the complex with the lowest RMSD and with minimal differences between the oscillations is the most stable complex [46,47].
Additionally, RMSD clustering was carried out (L1–TcTIM, B2–TcTIM and B5–TcTIM). A distribution by groups was clearly observed. Notably, two of them could be dynamic events, whereas the third group is best described as a conformational transition from one group to another (Figure 4) [36].
In addition, the frequencies of the interactions by group (see Appendix A) show that L1 is the compound that presented the highest number of different interactions in both groups, highlighting π-stacking interactions with Tyr100 (A) (99.75 and 81.58% presence over time in each group, respectively) and donors of hydrogen bonds with Tyr99 (A), among others such as cation-π and cationic interactions that remain for approximately 100% of the simulation time (Figure A2). In the case of BP5, the interaction frequency also shows hydrogen bond donor with Tyr99 (A) for 95.41% of the time and π-stacking interactions with Tyr100 (A) for 81.75% of the time, in addition to interactions with hydrogen bond acceptors Tyr100 (A) of 98%, which decreases in cluster 2 (to 62.65%) as well as the frequency of the π-stacking interaction with Tyr349 (B) (from 18.24 to 9.03%) to strengthen the interaction with Thr316 (B) (from 1.99 to 37.34%) (Figure A2).
In BP2, only three types of interaction are present in the three groups (hydrophobic, hydrogen bond donors and hydrogen bond acceptors), Tyr99 (A) being the only hydrogen bond donor interaction with an interaction frequency between 96 and 100% over time in the three clusters (Figure A2). It should be noted that in the interaction frequencies of each cluster, the number of interactions and the type of cation-π influence the stability of the union, as is the case with BP2, which does not present it [48].
The RMSF showed a very similar fluctuation pattern between the free protein and the protein in complex with the ligands, especially with the BP5 compound. However, the RMSF calculation revealed that the complex with the BP2 ligand showed a high fluctuation in some regions according to the RMSD pattern. Additionally, L1-TcTIM and BTS-TcTIM complexes showed a fluctuation of variable proportions.
In the calculation of RMSF, we observe that the high fluctuation in some regions ac-cording to the RMSD pattern may be due to mobile residues in both protein monomers that coincide with the regions of loop 3 (between residues 68 to 79), an adjacent moving α-helix (from residues 120 to 140), and loop 6 (residues 165 to 177), also known as the flexible or catalytic loop, which forms a kind of cap that opens and closes on the active site [15,37,49], being a factor that contributed to the elevation of RMSD values.
Additionally, the Rg value for all four complexes was slightly lower than that of free TcTIM. Therefore, it can be suggested that these complexes could be relatively more rigid than the free TcTIM. In addition, the similar Rg values in the four complexes suggest that the interactions of the ligands analyzed do not influence the variation of the protein structure (Figure 6).
MMGBSA analysis calculated the binding free energy of the ligand–TcTIM complex [50]. The van der Waals (ΔEvdw) forces form the greatest contribution to the free energy of binding in the three complexes (L1–TcTIM, BP2–TcTIM and BP5–TcTIM), especially with BP5. These kind of forces are more crucial than the electrostatic to determine binding [43,51]. The BP5–TcTIM complex presented the most favorable binding free energy (ΔGb) (−25.68 ± 0.19 Kcal/mol), followed by L1–TcTIM (−24.37 ± 0.21 Kcal/mol).
On the other hand, a low SASA energy (ΔGSA), in the case of the three complexes (L1–TcTIM, BP2–TcTIM and BP5–TcTIM), also indicates that the nonpolar residues are buried in the solvent, which may favor the stability of the complex through the synergistic effects of hydrogen bonding and hydrophobic interactions [50]. Additionally, the BTS–TcTIM complex was analyzed in the same way, where the ΔGSA was the lowest (−2.93 ± 0.01 Kcal/mol), favoring the stability of the complex and the contributions of ΔEvdw as well as the polar contribution (ΔGpolar), which were the ones that formed the highest contribution to ΔGb (−14.12 ± 0.23 Kcal/mol), which compared to the other compounds was the highest, validating the docking score obtained in molecular docking (−6.2 Kcal/mol).

3.5. Analysis of Molecular Physicochemical Properties

Compounds BP2 and BP5 were analyzed through the SwissADME website and the ProTox-II server. Neither compound violated any of the physicochemical properties, and, as such, they have drug-like characteristics [52]. The in silico analysis predicted high human intestinal absorption and null permeability of the blood–brain barrier for the two compounds, showing the potential for good oral absorption without causing central nervous system damage. In addition, they presented moderate solubility, an important aspect, since if the compound is not soluble in water, it cannot be absorbed [53].
On the other hand, the compounds showed hepatotoxic inactivity and were substrates of P-glycoprotein, an efflux pump related to drug resistance [54]. This glycoprotein takes the substance that is related to it and incorporates it into its structure to expel it— for example, through the intestine. This is a way in which the cell defends itself against toxic substances [55]. Many compounds can be a substrate of this glycoprotein—for example, it has been observed with antibiotics and anticancer agents [56]. In addition, the CYP450 inhibition prediction results showed that BP2 and BP5 compounds are likely to behave as inhibitors of almost all CYP450 isoforms (except isoform 1A2), potentially causing these isoforms to be unable to metabolize other drugs because they will be inhibited, maybe giving rise to drug–drug interactions [40]; therefore, it is necessary to improve these properties in these kinds of compounds.

3.6. Molecular Docking at the HsTIM Interface

To know the selectivity of control ligand L1 and the compounds BP2 and BP5 against TcTIM versus HsTIM, the molecular docking analysis at the HsTIM interface was performed. The results showed that the control ligand L1 exhibited a docking score of −7.2 Kcal/mol at the HsTIM interface. The orientation of the L1 ligand at the HsTIM interface allowed for hydrophobic interactions between the L1 and HsTIM ring systems with amino acid residues Arg17 (B), Asn71 (A), and Leu236 (B). Furthermore, the thiazole of ligand L1 formed a hydrogen bond between the thiazole ring nitrogen and guanidino nitrogen (C(NH2)+3). Compound BP2 showed a higher affinity, with a docking score of −8.0 Kcal/mol, than the control ligand L1 for the HsTIM interface. Due to its structure formed by a ring system, it allowed it to interact with many residues, such as Gly16 (B), Arg17 (B), Gln19 (B), Ser20 (B), Leu21 (B), Leu24 (B), Asn71 (A), Asp85 (A), and Leu236 (B). The putative position of the BP2 compound at the interface suggests that it may alter the stability of HsTIM monomers. Compound BP5 exhibited low affinity with a docking score of −5.9 Kcal/mol towards HsTIM. Likewise, the interactions observed were poor with respect to compound BP2 and the control ligand L1. Only three hydrogen bonding interactions with Gly16 (B) and Asn71 (A), and one hydrophobic contact with Thr70 (A), were observed.
Molecular dynamics simulations for the three complexes (L1–HsTIM, BP2–HsTIM, and BP5–HsTIM) and free–HsTIM were also analyzed (Figure 8). The trajectories were first analyzed by RMSD and RMSF analysis to understand the stability and fluctuations of these structures [50]. Subsequently, the Rg was analyzed to observe the compactness of the system, where BP5 was slightly lower than free-HsTIM, suggesting that this complex may be relatively more rigid. In addition, the free energies of binding were calculated by the MMGBSA method from the molecular dynamics trajectories. The ΔGb values for the three complexes (L1–HsTIM, BP2–HsTIM, and BP5–HsTIM) show lower values than the complexes with TcTIM, supporting the docking scores obtained in molecular docking. In general, these results suggested that BP5 showed lower affinity against HsTIM than LI and BP2.

4. Materials and Methods

4.1. Data Preparation

10 compounds with inhibitory activity against TcTIM were used as control ligands, which were drawn in ChemDraw and saved in SDF format. Afterwards, they were minimized and converted to pdbqt format from OpenBabel.

4.2. Molecular Docking

The crystallographic structure of the TcTIM protein in complex with 3-(2-benzothiazolylthio-1-propanesulfonic acid (BTS) was obtained from the Protein Data Bank (PDB) “http://www.pdb.org (accessed on 12 May 2021)” [57] with the PDB ID 1SUX as well as the HsTIM protein with the accession code 4POC. The protein was prepared for docking using the UCSF Chimera 1.14.1 DockPrep tool [58]. Additionally, the prepare_receptor4.py script from MGLTools 1.5.6 was used to add AutoDock atom types and add Gasteiger charges to the protein structure.
The prediction of potential binding sites was performed first using the DoGSiteScorer tool of the Proteins Plus server “https://proteins.plus/ (accessed on 18 May 2021)”. Subsequently, a blind molecular docking was performed. For blind molecular docking, the receptor was defined as rigid and PyRx software was used, which works with AutoDock vina 1.1.2 (vina) [59]. For docking at the binding site, the conformational search space was determined by establishing the coordinates in the center of residues at the interface (X = 29,077, Y = 101,654, and Z = 62,161) using PyRx software. The binding site on HsTIM was determined by overlap between TcTIM (PDB 1SUX) and the HsTIM apoprotein (PDB 4POC) using UCSF Chimera. Then, docking at HsTIM was centered on interface residues, as previously described [56].

4.3. Virtual Screening

LBVS using benzimidazole scaffold was performed from the ZINC15 database “https://zinc15.docking.org (accessed on 25 May 2021)” [60]. Subsequently, all structures were downloaded in “CSV” format and the Lipinski’s rule of five was applied using the OpenBabel program. Finally, the structures were prepared for molecular docking on the TcTIM interface using the PyRx program. Once molecular docking at the TcTIM interface was concluded, all compounds were selected based on the docking score of −8.9 Kcal/mol obtained by the L1 ligand (see results section). Subsequently, through the PLIP web service, an interaction profile was generated for each of the docking complexes [61]. With the scikit-learn library and the DataWarrior program “https://openmolecules.org/datawarrior (accessed on 8 June 2021)” [33,62], the compounds were grouped according to the interaction profiles. Finally, two compounds were selected based on cost and availability for evaluation against bloodstream trypomastigotes of T. cruzi.

4.4. Trypanocidal Assay

CD1 mice, 6 to 8 weeks old, infected with bloodstream trypomastigotes of INC-5 and NINOA strain, were used for the trypanocidal assay. At the peak of parasitemia (2 to 4 weeks), parasitized blood was collected by cardiac puncture using sodium heparin as an anticoagulant. Blood was adjusted to 1 × 106 trypomastigotes/mL. In a 96-well plate, 90 µL of infected blood and 10 µL of benzimidazole derivatives or dilutions of the reference drugs (Bzn and Nfx) were deposited for a final volume of 100 µL per well. Reference drugs were used as positive lysis control, and wells with untreated blood trypomastigotes were used as negative lysis control; the microplates were incubated at 4 °C for 24 h. Motile trypomastigotes were subsequently quantified using the Pizzi–Brener method [63], for which 5 μL of blood was placed on a slide and covered with a 18 × 18 mm coverslip. Motile protozoa were counted in 15 fields at 40× using an optical microscope. The lysis percentage of each treatment was calculated by comparing the viable trypomastigotes with the negative control [64]. The half-maximum lytic concentration (LC50) was determined for each compound using the Probit statistical tool. Subsequently, the results were converted to micromolar units. Benzimidazole derivatives (BP2 and BP5) were purchased from the commercial house MolPort and were worked without any additional purification. The assay was performed in triplicate.

4.5. Molecular Dynamics Analysis

For analysis of molecular dynamics of the selected compounds, the open-source software suite GROMACS 5.1.2 was used [65]. Protein with access code 1SUX was parameterized in AMBER03 (ff94/ff99 modification from Duan et al., 2003) [66] force field with the GROMACS software suite pdb2gmx. Protein protonation state pH 7 was previously calculated using PROPKA tool implemented in UCSF Chimera. On the other hand, the topology of the compounds was generated with ACPYPE Server “http://webapps.ccpn.ac.uk/acpype (accessed on 14 July 2021)” [67], which is based on the General Amber Force Field (GAFF). The system was a dodecahedron with periodic boundary conditions. In addition to containing the ligand-protein complex, it was filled with TIP3P water molecules and the number of ions (Cl− or Na+) necessary to have a neutral charge in the system. Before performing dynamics, the system was energetically minimized by the steepest descent algorithm. Then, two equilibrium steps were performed with restrictions of 1000 kJ/mol/nm2 on the movement of protein and ligand heavy atoms. The first stage was at constant pressure, implementing the leapfrog method and the v-rescale thermostat to bring the system from 0 to 300 K. The second stage was carried out at constant temperature, again with the frog jump method, but now with the Berendsen barostat to bring the system from 1 to 2 bar. Both stages achieved a duration of 100 ps. Once the system was balanced, molecular dynamics were performed with a 100 ns trajectory, where long-range interactions and forces were calculated with the Particle-Mesh Ewald (PME) method, establishing the Lennard–Jones and Coulomb contributions at 1.2 nm, balancing the system with samples at 100 ps. Finally, the interactions for each complex in the simulation were obtained with the ProLIF tool (Protein–Ligand Interaction Fingerprints) [68]. The determination of the percentage of the frequencies of interactions was carried out according to the grouping obtained by each complex. The stability of the complexes was determined using the GROMACS software tools. First, the RMSD between the α carbons and the ligand was obtained. Then, the pairwise RMSD matrix was calculated during the 100 ns. It was performed with MDAnalysis [69] to make the RMSD clustering with the metric ward and the Euclidean distance from the python3 SciPy library. RMSF of α-carbons along with 2D structure [70] was performed to understand the effect of the compound on the secondary structure of the TIM surface. Finally, the Rg was calculated to corroborate the maintenance of the three-dimensional compactness of the TIM throughout the simulation. The free energy of binding of the complexes was calculated with the MMPBSA using the g_mmpbsa software [71], retrieving 200 stable frames (between 90 and 100 ns).

4.6. Analysis of Molecular Physicochemical Properties

In addition, the compounds BP2 and BP5 were in silico analyzed with the SwissADME website “http://www.swissadme.ch (accessed on 22 August 2021)” [72] to determine their pharmacokinetic properties, such as absorption, bioavailability, permeability, blood–brain barrier penetration, metabolism, and excretion, as well as hepatotoxicity via the ProTox-II server “https://tox-new.charite.de/protox_II (accessed on 23 August 2022)” [73].

5. Conclusions

In this study, a LBVS using a benzimidazole scaffold and molecular docking at the TcTIM interface allowed obtaining 1604 new benzimidazole derivatives as potential TcTIM inhibitors. The biological evaluation in the in vitro model against the trypomastigote form of T. cruzi determined that compounds BP2 and BP5 have trypanocidal activity comparable to reference drugs. Additionally, the molecular dynamics simulation demonstrated that compound BP5 formed the most stable complex with TcTIM during the time analyzed, and van der Waals forces were described as the most important in the binding process. Finally, a low affinity (−5.9 Kcal/mol) of BP5 at the HsTIM interface suggest a selectivity on TcTIM. These results encourage researchers to explore these kinds of compounds to develop new trypanocidal agents.

Author Contributions

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

Funding

This research was funded by Secretaria de Investigacion y Posgrado del Instituto Politecnico Nacional, grant number 20220935.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Ijms 23 10047 g0a1 550
Figure A1. Matrices for RMSD cluster analysis.
Figure A1. Matrices for RMSD cluster analysis.
Ijms 23 10047 g0a1
Ijms 23 10047 g0a2 550
Figure A2. RMSF cluster analysis.
Figure A2. RMSF cluster analysis.
Ijms 23 10047 g0a2
Ijms 23 10047 g0a3 550
Figure A3. Rg cluster analysis.
Figure A3. Rg cluster analysis.
Ijms 23 10047 g0a3
Ijms 23 10047 g0a4 550
Figure A4. Frequency of interactions present in each cluster by RSMD.
Figure A4. Frequency of interactions present in each cluster by RSMD.
Ijms 23 10047 g0a4

References

  1. Zamora-Ledesma, S.; Hernández-Camacho, N.; Sánchez-Moreno, M.; Ruiz-Piña, H.; Villagrán-Herrera, M.E.; Marín-Sánchez, C.; Carrillo-Angeles, I.G.; Jones, R.W.; Camacho-Macías, B. Seropositivity for Trypanosoma cruzi and Leishmania mexicana in dogs from a metropolitan region of Central México. Vet. Parasitol. Reg. Stud. Rep. 2020, 22, 100459. [Google Scholar] [CrossRef] [PubMed]
  2. PAHO Pan-American Health Organization. Available online: http://www.paho.org/hq/index.php?option=com_topics&view=article&id=10&Itemid=40743 (accessed on 15 April 2021).
  3. Nunes, M.C.; Dones, W.; Morillo, C.A.; Encina, J.J.; Ribeiro, A.L. Chagas disease: An overview of clinical and epidemiological aspects. J. Am. Coll. Cardiol. 2013, 62, 767–776. [Google Scholar] [CrossRef] [PubMed]
  4. Olivares-Illana, V.; Riveros-Rosas, H.; Cabrera, N.; de Gómez-Puyou, M.; Pérez-Montfort, R.; Costas, M.; Gómez-Puyou, A. A guide to the effects of a large portion of the residues of triosephosphate isomerase on catalysis, stability, druggability, and human disease. Proteins 2017, 85, 1190–1211. [Google Scholar] [CrossRef] [PubMed]
  5. Vázquez-Raygoza, A.; Cano-González, L.; Velázquez-Martínez, I.; Trejo-Soto, P.J.; Castillo, R.; Hernández-Campos, A.; Hernández-Luis, F.; Oria-Hernández, J.; Castillo-Villanueva, A.; Avitia-Domínguez, C.; et al. Species-Specific Inactivation of Triosephosphate Isomerase from Trypanosoma brucei: Kinetic and Molecular Dynamics Studies. Molecules 2017, 22, 2055. [Google Scholar] [CrossRef] [PubMed]
  6. Olivares-Illana, V.; Pérez-Montfort, R.; López-Calahorra, F.; Costas, M.; Rodríguez-Romero, A.; Tuena de Gómez-Puyou, M.; Gómez Puyou, A. Structural differences in triosephosphate isomerase from different species and discovery of a multitrypanosomatid inhibitor. Biochemistry 2006, 45, 2556–2560. [Google Scholar] [CrossRef] [PubMed]
  7. D’Antonio, E.L.; Deinema, M.S.; Kearns, S.P.; Frey, T.A.; Tanghe, S.; Perry, K.; Roy, T.A.; Gracz, H.S.; Rodriguez, A.; D’Antonio, J. Structure-based Approach to the Identification of a Novel Group of Selective Glucosamine Analogue Inhibitors of Trypanosoma Cruzi Glucokinase. Mol. Biochem. Parasitol. 2015, 204, 64–76. [Google Scholar] [CrossRef]
  8. Maldonado, E.; Soriano-García, M.; Moreno, A.; Cabrera, N.; Garza-Ramos, G.; Tuena de Gómez-Puyou, M.; Gómez-Puyou, A.; Pérez-Montfort, R. Differences in the intersubunit contacts in triosephosphate isomerase from two closely related pathogenic trypanosomes. J. Mol. Biol. 1998, 283, 193–203. [Google Scholar] [CrossRef]
  9. Zomosa-Signoret, V.; Hernández-Alcántara, G.; Reyes-Vivas, H.; Martínez-Martínez, E.; Garza-Ramos, G.; Pérez-Montfort, R.; de Gómez-Puyou, M.T.; Gómez-Puyou, A. Control of the reactivation kinetics of homodimeric triosephosphate isomerase from unfolded monomers. Biochemistry 2003, 42, 3311–3318. [Google Scholar] [CrossRef]
  10. Pérez-Montfort, R.; Garza-Ramos, G.; Alcántara, G.; Reyes-Vivas, H.; Gao, X.G.; Maldonado, E.; de Gómez-Puyou, M.T.; Gómez-Puyou, A. Derivatization of the Interface Cysteine of Triosephosphate Isomerase from Trypanosoma Brucei and Trypanosoma Cruzi as Probe of the Interrelationship Between the Catalytic Sites and the Dimer Interface. Biochemistry 1999, 38, 4114–4120. [Google Scholar] [CrossRef]
  11. Álvarez, G.; Aguirre-López, B.; Varela., J.; Cabrera., M.; Merlino., A.; López., G.V.; Lavaggi., M.L.; Porcal, W.; Di-Maio, R.; González, M.; et al. Massive screening yields novel and selective Trypanosoma cruzi triosephosphate isomerase dimer-interface irreversible inhibitors with anti-trypanosomal activity. Eur. J. Med. Chem. 2010, 45, 5767–5772. [Google Scholar] [CrossRef]
  12. Flores-Sandoval, C.A.; Cuevas-Hernández, R.I.; CorreaBasurto, J.; Beltrán-Conde, H.I.; Padilla-Martínez, I.I.; Farfán-García, J.N.; Nogueda-Torres, B.; Trujillo-Ferrara, J.G. Synthesis and theoretic calculations of benzoxazoles and docking studies of their interactions with triosephosphate isomerase. Med. Chem. Res. 2013, 22, 2768–2777. [Google Scholar] [CrossRef]
  13. Cuevas-Hernández, R.I.; Correa-Basurto, J.; FloresSandoval, C.A.; Padilla-Martínez, I.I.; Nogueda-Torres, B.; Villa-Tanaca, M.L.; Tamay-Cach, F.; Nolasco-Fidencio, J.J.; Trujillo-Ferrara, J.G. Fluorine-containing benzothiazole as a novel trypanocidal agent: Design, in silico study, synthesis and activity evaluation. Med. Chem. Res. 2016, 25, 211–224. [Google Scholar] [CrossRef]
  14. Kurkcuoglu, Z.; Findik, D.; Akten, E.D.; Doruker, P. How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics. Biophys. J. 2015, 109, 1169–1178. [Google Scholar] [CrossRef] [PubMed]
  15. Wierenga, R.K.; Noble, M.E.; Vriend, G.; Nauche, S.; Hol, W.G. Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex. J. Mol. Biol. 1991, 220, 995–1015. [Google Scholar] [CrossRef]
  16. Wierenga, R.K.; Noble, M.E. Comparison of the refined crystal structures of liganded and unliganded chicken, yeast and trypanosomal triosephosphate isomerase. J. Mol. Biol. 1992, 224, 1115–1126. [Google Scholar] [CrossRef]
  17. Álvarez, G.; Martínez, J.; Varela, J.; Birriel, E.; Cruces, E.; Gabay, M.; Leal, S.M.; Escobar, P.; Aguirre-López, B.; Cabrera, N.; et al. Development of bis-thiazoles as inhibitors of triosephosphate isomerase from Trypanosoma cruzi. Identification of new non-mutagenic agents that are active in vivo. Eur. J. Med. Chem. 2015, 100, 246–256. [Google Scholar] [CrossRef]
  18. Shaaban, M.R.; Saleh, T.S.; Mayhoub, A.S.; Mansour, A.; Farag, A.M. Synthesis and analgesic/anti-inflammatory evaluation of fused heterocyclic ring systems incorporating phenylsulfonyl moiety. Bioorg. Med. Chem. 2008, 16, 6344–6352. [Google Scholar] [CrossRef]
  19. El Rashedy, A.A.; Aboul-Enein, H.Y. Benzimidazole derivatives as potential anticancer agents. Mini Rev. Med. Chem. 2013, 13, 399–407. [Google Scholar] [CrossRef]
  20. Song, D.; Shutao, M. Recent Development of Benzimidazole-Containing Antibacterial Agents. Chem. Med. Chem. 2016, 11, 646–659. [Google Scholar] [CrossRef]
  21. Francesconi, V.; Cichero, E.; Schenone, E.; Naesens, L.; Tonelli, M. Synthesis and Biological Evaluation of Novel (thio) semicarbazone-Based Benzimidazoles as Antiviral Agents against Human Respiratory Viruses. Molecules 2020, 25, 1487. [Google Scholar] [CrossRef] [Green Version]
  22. Romo-Mancillas, A.; Téllez-Valencia, A.; Yépez-Mulia, L.; Hernández-Luis, F.; Hernández-Campos, A.; Castillo, R. The design and inhibitory profile of new benzimidazole derivatives against triosephosphate isomerase from Trypanosoma cruzi: A problem of residue motility. J. Mol. Graph. Model. 2011, 30, 90–99. [Google Scholar] [CrossRef]
  23. Oh, S.; Kim, S.; Kong, S.; Yang, G.; Lee, N.; Han, D.; Goo, J.; Siqueira-Neto, J.L.; Freitas-Junior, L.H.; Song, R. Synthesis and biological evaluation of 2,3-dihydroimidazo[1,2-a]benzimidazole derivatives against Leishmania donovani and Trypanosoma cruzi. Eur. J. Med. Chem. 2014, 84, 395–403. [Google Scholar] [CrossRef]
  24. Melchor-Doncel de la Torre, S.; Vázquez, C.; González-Chávez, Z.; Yépez-Mulia, L.; Nieto-Meneses, R.; Jasso-Chávez, R.; Saavedra, M.; Hernández-Luis, F. Synthesis and biological evaluation of 2-methyl-1H-benzimidazole-5-carbohydrazides derivatives as modifiers of redox homeostasis of Trypanosoma cruzi. Bioorg. Med. Chem. Lett. 2017, 27, 3403–3407. [Google Scholar] [CrossRef] [PubMed]
  25. Velázquez-López, J.M.; Hernández-Campos, A.; Yépez-Mulia, L.; Téllez-Valencia, A.; Flores-Carrillo, P.; Nieto-Meneses, R.; Castillo, R. Synthesis and trypanocidal activity of novel benzimidazole derivatives. Eur. J. Med. Chem. Lett. 2016, 26, 4377–4381. [Google Scholar] [CrossRef] [PubMed]
  26. Bezabeh, S.; Mackey, A.C.; Kluetz, P.; Jappar, D.; Korvick, J. Accumulating Evidence for a Drug-Drug Interaction between Methotrexate and Proton Pump Inhibitors. Oncologist 2012, 17, 550–554. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, M.; Han, X.; Zhou, Z. New Substituted Benzimidazole Derivatives: A Patent Review (2013–2014). Expert Opin. Ther. Pat. 2015, 25, 595–612. [Google Scholar] [CrossRef] [PubMed]
  28. Hossain, M.J.; Sultan, M.Z.; Rashid, M.A.; Kuddus, M.R. Does Rabeprazole Sodium Alleviate the Anti-diabetic Activity of Linagliptin? Drug-Drug Interaction Analysis by In Vitro and In Vivo Methods. Drug Res. 2020, 70, 519–527. [Google Scholar] [CrossRef]
  29. Hossain, M.J.; Sultan, M.Z.; Rashid, M.A.; Kuddus, M.R. Interactions of Linagliptin, Rabeprazole Sodium, and Their Formed Complex with Bovine Serum Albumin: Computational Docking and Fluorescence Spectroscopic Methods. Anal. Sci. Adv. 2021, 7, 202000153. [Google Scholar] [CrossRef]
  30. LiverTox: Clinical and Research Information on Drug-Induced Liver Injury [Internet]; Anthelmintic Agents; National Institute of Diabetes and Digestive and Kidney Diseases: Bethesda, MD, USA, 2012.
  31. Kurkcuoglu, Z.; Ural, G.; Demet Akten, E.; Doruker, P. Blind Dockings of Benzothiazoles to Multiple Receptor Conformations of Triosephosphate Isomerase from Trypanosoma cruzi and Human. Mol. Inform. 2011, 30, 986–995. [Google Scholar] [CrossRef]
  32. Sander, T.; Freyss, J.; von Korff, M.; Rufener, C. DataWarrior: An open-source program for chemistry aware data visualization and analysis. J. Chem. Inf. Model. 2015, 55, 460–473. [Google Scholar] [CrossRef]
  33. Abraham, A.; Pedregosa, F.; Eickenberg, M.; Gervais, P.; Mueller, A.; Kossaifi, J.; Gramfort, A.; Thirion, B.; Varoquaux, G. Machine learning for neuroimaging with scikit-learn. Front. Neuroinform. 2014, 8, 14. [Google Scholar] [CrossRef] [PubMed]
  34. De Paris, R.; Quevedo, C.V.; Ruiz, D.D.; Norberto de Souza, O. An effective approach to cluster inhA molecular dynamics pathway using substrate-binding cavity features. PLoS ONE 2015, 10, e0133172. [Google Scholar] [CrossRef] [PubMed]
  35. De Paris, R.; Quevedo, C.V.; Ruiz, D.D.; Norberto de Souza, O.; Barros, R.C. Clustering molecular dynamics trajectories for optimizing docking experiments. Comput. Intell. Neurosci. 2015, 2015, 1–9. [Google Scholar] [CrossRef] [PubMed]
  36. Wolf, A.; Kirschner, K.N. Principal component and clustering analysis on molecular dynamics data of the ribosomal L11·23S subdomain. J. Mol. Model. 2013, 19, 539–549. [Google Scholar] [CrossRef] [PubMed]
  37. Singh, S. Dynamics of heroin molecule inside the lipid membrane: A molecular dynamics study. J. Mol. Model. 2019, 25, 121. [Google Scholar] [CrossRef]
  38. Minini, L.; Álvarez, G.; González, M.; Cerecetto, H.; Merlino, A. Molecular docking and molecular dynamics simulation studies of Trypanosoma cruzi triosephosphate isomerase inhibitors. Insights into the inhibition mechanism and selectivity. J. Mol. Graph. Model. 2015, 58, 40–49. [Google Scholar] [CrossRef]
  39. Téllez-Valencia, A.; Olivares-Illana, V.; Hernandez-Santoyo, A.; Perez-Montfort, R.; Costas, M.; Rodriguez-Romero, A.; López-Calahorra, F.; de Gómez-Puyou, M.T.; Gómez-Puyou, A. Inactivation of triosephosphate isomerase from Trypanosoma cruzi by an agent that perturbs its dimer interface. J. Mol. Biol. 2004, 341, 1355–1365. [Google Scholar] [CrossRef]
  40. Juárez-Saldivar, A.; Barbosa-Cabrera, E.; Lara-Ramírez, E.E.; Paz-González, A.D.; Martínez-Vázquez, A.V.; Bocanegra-García, V.; Palos, I.; Campillo, N.E.; Rivera, G. Virtual Screening of FDA-Approved Drugs against Triose Phosphate Isomerase from Entamoeba histolytica and Giardia lamblia Identifies Inhibitors of Their Trophozoite Growth Phase. Int. J. Mol. Sci. 2021, 22, 5943. [Google Scholar] [CrossRef]
  41. Téllez-Valencia, A.; Avila-Ríos, S.; Pérez-Montfort, R.; Rodríguez-Romero, A.; de Gómez-Puyou, M.T.; López-Calahorra, F.; Gómez-Puyou, A. Highly Specific Inactivation of Triosephosphate Isomerase from Trypanosoma Cruzi. Biochem. Biophys Res. Commun. 2002, 295, 958–963. [Google Scholar] [CrossRef]
  42. Unni, S.; Aouti, S.; Thiyagarajan, S.; Padmanabhan, B. Identification of a repurposed drug as an inhibitor of Spike protein of human coronavirus SARS-CoV-2 by computational methods. J. biosci. 2020, 45, 130. [Google Scholar] [CrossRef]
  43. Espinosa-Bustos, C.; Ortiz-Pérez, M.; González-González, 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]
  44. Perilla, J.R.; Goh, B.C.; Cassidy, C.K.; Liu, B.; Bernardi, R.C.; Rudack, T.; Yu, H.; Wu, Z.; Schulten, K. Molecular dynamics simulations of large macromolecular complexes. Curr. Opin. Struct. Biol. 2015, 31, 64–74. [Google Scholar] [CrossRef] [PubMed]
  45. Méndez-Álvarez, D.; Herrera-Mayorga, D.; Juárez-Saldivar, A.; Paz-González, A.D.; Ortiz-Pérez, E.; Bandyopadhyay, D.; Pérez-Sánchez, H.; Rivera, G. Ligand-based virtual screening, molecular docking, and molecular dynamics of eugenol analogs as potential acetylcholinesterase inhibitors with biological activity against Spodoptera frugiperda. Mol. Divers. 2021, 26, 2025–2037. [Google Scholar] [CrossRef]
  46. De Vivo, M.; Masetti, M.; Bottegoni, G.; Cavalli, A. Role of molecular dynamics and related methods in drug discovery. J. Med. Chem. 2016, 59, 4035–4061. [Google Scholar] [CrossRef] [PubMed]
  47. Zhang, X.; Yan, J.; Wang, H.; Wang, Y.; Wang, J.; Zhao, D. Molecular docking, 3D-QSAR, and molecular dynamics simulations of thieno [3, 2-b] pyrrole derivatives against anticancer targets of KDM1A/LSD1. J. Biomole. Struct. Dyn. 2020, 59, 4035–4061. [Google Scholar] [CrossRef]
  48. Luise, C.; Robaa, D.; Regenass, P.; Maurer, D.; Ostrovskyi, D.; Seifert, L.; Bacher, J.; Burgahn, T.; Wagner, T.; Seitz, J.; et al. Structure-Based Design, Docking and Binding Free Energy Calculations of A366 Derivatives as Spindlin1 Inhibitors. Int. J. Mol. Sci. 2021, 22, 5910. [Google Scholar] [CrossRef]
  49. Williams, J.C.; McDermott, A.E. Dynamics of the flexible loop of triosephosphate isomerase: The loop motion is not ligand gated. Biochemistry 1995, 34, 8309–8319. [Google Scholar] [CrossRef]
  50. Ghosh, R.; Chakraborty, A.; Biswas, A.; Chowdhuri, S. Evaluation of green tea polyphenols as novel corona virus (SARS-CoV-2) main protease (Mpro) inhibitors-an in silico docking and molecular dynamics simulation study. J. Biomol. Struct. Dyn. 2021, 39, 4362–4374. [Google Scholar] [CrossRef]
  51. Niu, Y.; Pan, D.; Shi, D.; Bai, Q.; Liu, H.; Yao, X. Influence of Chirality of Crizotinib on Its MTH1 Protein Inhibitory Activity: Insight from Molecular Dynamics Simulations and Binding Free Energy Calculations. PLoS ONE 2015, 10, e0145219. [Google Scholar] [CrossRef] [PubMed]
  52. Benet, L.Z.; Hosey, C.M.; Ursu, O.; Oprea, T.I. BDDCS, the Rule of 5 and drugability. Adv. Drug Deliv. Rev. 2016, 101, 89–98. [Google Scholar] [CrossRef] [Green Version]
  53. Ibrahim, Z.Y.; Uzairu, A.; Shallangwa, G.A.; Abechi, S.E. Application of QSAR Method in the Design of Enhanced Antimalarial Derivatives of Azetidine-2-carbonitriles, their Molecular Docking, Drug-likeness, and SwissADME Properties. Iran J. Pharm. Res. 2021, 20, 254–270. [Google Scholar] [CrossRef]
  54. Mollazadeh, S.; Sahebkar, A.; Hadizadeh, F.; Behravan, J.; Arabzadeh, S. Structural and functional aspects of P-glycoprotein and its inhibitors. Life Sci. 2018, 214, 118–123. [Google Scholar] [CrossRef] [PubMed]
  55. Glaeser, H. Importance of P-glycoprotein for drug-drug interactions. Handb. Exp. Pharmacol. 2011, 201, 285–297. [Google Scholar] [CrossRef]
  56. Waghray, D.; Zhang, Q. Inhibit or Evade Multidrug Resistance P-Glycoprotein in Cancer Treatment. J. Med. Chem. 2018, 61, 5108–5121. [Google Scholar] [CrossRef]
  57. Berman, H.M.; Kleywegt, G.J.; Nakamura, H.; Markley, J.L. How community has shaped the Protein Data Bank. Structure 2013, 21, 1485–1491. [Google Scholar] [CrossRef]
  58. 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]
  59. 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]
  60. Sterling, T.; Irwin, J.J. ZINC 15--Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [Google Scholar] [CrossRef]
  61. Salentin, S.; Schreiber, S.; Joachim-Haupt, V.; Adasme, M.F.; Schroeder, M. PLIP: Fully automated protein–Ligand interaction profiler. Nucleic Acids Res. 2015, 43, W443–W447. [Google Scholar] [CrossRef]
  62. López-López, E.; Naveja, J.J.; Medina-Franco, J.L. DataWarrior: An evaluation of the open-source drug discovery tool. Expert Opin. Drug Discov. 2019, 14, 335–341. [Google Scholar] [CrossRef]
  63. Chacón-Vargas, K.F.; Nogueda-Torres, B.; Sánchez-Torres, L.E.; Suarez-Contreras, E.; Villalobos-Rocha, J.C.; Torres-Martinez, Y.; Lara-Ramírez, 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]
  64. Becerra, N.A.; Espinosa-Bustos, C.; Vázquez, K.; Rivera, G.; Paulino, M.; Cantero, J.; Nogueda, B.; Chacón-Vargas, F.; Castillo-Velazquez, U.; Elizondo Rodríguez, A.F.; 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]
  65. Van-Der-Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  66. Duan, Y.; Wu, C.; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; Lee, T.; et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comput. Chem. 2003, 24, 1999–2012. [Google Scholar] [CrossRef]
  67. Sousa-da-Silva, A.W.; Wim Vranken, F. ACPYPE AnteChamber PYthon Parser interface. BMC Res. Notes 2012, 5, 367. [Google Scholar] [CrossRef] [PubMed]
  68. Bouysset, C.; Fiorucci, S. ProLIF: A library to encode molecular interactions as fingerprints. J. Cheminform. 2021, 13, 1–9. [Google Scholar] [CrossRef]
  69. Michaud-Agrawal, N.; Denning, E.J.; Woolf, T.B.; Beckstein, O. MDAnalysis: A toolkit for the analysis of molecular dynamics simulations. J. Comput. Chem. 2011, 32, 2319–2327. [Google Scholar] [CrossRef] [PubMed]
  70. Kunzmann, P.; Hamacher, K. Biotite: A unifying open-source computational biology framework in Python. BMC Bioinform. 2018, 19, 1–8. [Google Scholar] [CrossRef]
  71. Kumari, R.; Kumar, R.; Open Source Drug Discovery Consortium; Lynn, A. g_mmpbsa—A GROMACS Tool for High-Throughput MM-PBSA Calculations. J. Chem. Inf. Model. 2014, 54, 1951–1962. [Google Scholar] [CrossRef] [PubMed]
  72. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  73. Banerjee, P.; Eckert, A.O.; Schrey, A.K.; Preissner, R. ProTox-II: A webserver for the prediction of toxicity of chemicals. Nucleic Acids Res. 2018, 46, W257–W263. [Google Scholar] [CrossRef] [Green Version]
Ijms 23 10047 g001 550
Figure 1. Binding site prediction by the DoGSiteScorer web Server and blind molecular docking in TcTIM. The protein is represented in gray, the interface is represented in green, and the ligands in red crosses.
Figure 1. Binding site prediction by the DoGSiteScorer web Server and blind molecular docking in TcTIM. The protein is represented in gray, the interface is represented in green, and the ligands in red crosses.
Ijms 23 10047 g001
Ijms 23 10047 g002 550
Figure 2. Lead compound of each group obtained by interaction profile analysis. Hydrophobic interactions are shown in gray, hydrogen bonds in blue, π-stacking in green, and cation-π interactions in purple.
Figure 2. Lead compound of each group obtained by interaction profile analysis. Hydrophobic interactions are shown in gray, hydrogen bonds in blue, π-stacking in green, and cation-π interactions in purple.
Ijms 23 10047 g002
Ijms 23 10047 g003 550
Figure 3. RMSD of molecular dynamics simulation of free TcTIM, benzimidazole derivatives, and the BTS ligand in complex with TcTIM.
Figure 3. RMSD of molecular dynamics simulation of free TcTIM, benzimidazole derivatives, and the BTS ligand in complex with TcTIM.
Ijms 23 10047 g003
Ijms 23 10047 g004 550
Figure 4. RMSD clustering of molecular dynamics simulation of free TcTIM and benzimidazole derivatives in complex with TcTIM.
Figure 4. RMSD clustering of molecular dynamics simulation of free TcTIM and benzimidazole derivatives in complex with TcTIM.
Ijms 23 10047 g004
Ijms 23 10047 g005 550
Figure 5. RMSF of molecular dynamics simulation of free TcTIM, benzimidazole derivatives, and the BTS ligand in complex with TcTIM.
Figure 5. RMSF of molecular dynamics simulation of free TcTIM, benzimidazole derivatives, and the BTS ligand in complex with TcTIM.
Ijms 23 10047 g005
Ijms 23 10047 g006 550
Figure 6. Rg of molecular dynamics simulation of free TcTIM, benzimidazole derivatives and the BTS ligand in complex with TcTIM.
Figure 6. Rg of molecular dynamics simulation of free TcTIM, benzimidazole derivatives and the BTS ligand in complex with TcTIM.
Ijms 23 10047 g006
Ijms 23 10047 g007 550
Figure 7. Interaction profile of L1, BP2, and BP5 at the HsTIM interface. Hydrogen bonds are shown as blue lines and hydrophobic interaction as gray dashed lines. Yellow dashed lines represent Salt Bridge interaction. Interactions were obtained using the PLIP online server.
Figure 7. Interaction profile of L1, BP2, and BP5 at the HsTIM interface. Hydrogen bonds are shown as blue lines and hydrophobic interaction as gray dashed lines. Yellow dashed lines represent Salt Bridge interaction. Interactions were obtained using the PLIP online server.
Ijms 23 10047 g007
Ijms 23 10047 g008 550
Figure 8. Analysis of molecular dynamics simulations of L1–, BP2– and BP5–HsTIM complex. (A) RMSD diagram of the three complexes. (B) RMSF diagram of the three complexes. (C) Rg diagram of the three complexes.
Figure 8. Analysis of molecular dynamics simulations of L1–, BP2– and BP5–HsTIM complex. (A) RMSD diagram of the three complexes. (B) RMSF diagram of the three complexes. (C) Rg diagram of the three complexes.
Ijms 23 10047 g008
Table
Table 1. Docking scores and molecular interactions of TcTIM inhibitory compounds used as control ligands at the TcTIM interface.
Table 1. Docking scores and molecular interactions of TcTIM inhibitory compounds used as control ligands at the TcTIM interface.
IDCompoundAntiparasitic ActivityEnzymatic
Activity		Docking Score
Kcal/mol		Hydrophobic InteractionsHydrogen Bondsπ-Stacking Interactions
L1Ijms 23 10047 i001-% inhibition =
48 at 200 µM
[20]		−8.9Tyr102 (A)
Ile69 (B)
Arg71 (B)
Phe75 (B)		Glu75 (A)
Tyr102 (A)
Arg71 (B)		Tyr103 (A)
L2Ijms 23 10047 i002-% inhibition =
52 at 200 µM [21]		−6.8Tyr102 (A)
Ile69 (B)		Leu101 (A)
Tyr102 (A)
Tyr103 (A)
Arg71 (B)		Tyr103 (A)
L3Ijms 23 10047 i003Epimastigotes IC50 = 186.23 μM [25]% inhibition =
65 at 100 µM
[21]		−6.9Tyr102 (A)
Ile69 (B)		Tyr103 (A)
Arg71 (B)		Tyr103 (A)
Tyr103 (B)
L4Ijms 23 10047 i004Epimastigotes IC50 = 48.91 μM [25]% inhibition =
59 at 200 µM
[25]		−7.3Tyr102 (A)
Ile69 (B)		Tyr103 (A)
Arg71 (B)		Tyr103 (A)
Tyr103 (B)
L5Ijms 23 10047 i005Epimastigotes IC50 = 42.51 μM [25]% inhibition =
69 at 200 µM [25]		−7.3Tyr102 (A)
Ile69 (B)		Tyr103 (A)
Arg71 (B)		Tyr103 (A)
L6Ijms 23 10047 i006Epimastigotes IC50 = 28.67 μM [25]% de inhibición =
68 at 200 µM [25]		−7.2Tyr102 (A)
Ile69 (B)
Phe75 (B)		Tyr103 (A)
Tyr103 (B)
Gly104 (B)		Tyr103 (A)
L7Ijms 23 10047 i007-IC50 = 8.0 µM;
% inhibition = 95 at 250 µM
[26,27]		−6.9Tyr102 (A)
Ile69 (B)
Glu108 (B)
Ile109 (B)		Tyr103 (A)
Glu112 (B)		Tyr103 (A)
L8Ijms 23 10047 i008Epimastigotes IC50 = 5.8 µM
[28]		IC50 = 3.5 µM
[28]		−6.9Tyr102 (A)
Ile69 (B)
Phe75 (B)		Tyr102 (A)
Tyr103 (B)		Tyr103 (A)
L9Ijms 23 10047 i009Epimastigotes IC50 = 0.6 µM
[29]		IC50 = 0.086 µM
[29]		−7.2Tyr102 (A)
Tyr103 (A)
Ile69 (B)
Arg71 (B)
Phe75 (B)
Tyr103 (B)
Glu105 (B)		Arg71 (B)Tyr102 (A)
BTSIjms 23 10047 i010-semi-maximal inhibition at 35 μM [31]−6.2Tyr102 (A)
Ile69 (B)
Phe75 (B)		Tyr103 (B)Tyr103 (A)
Table
Table 2. Lead compounds from each group obtained by interaction profile analysis of 1604 potential TcTIM inhibitors.
Table 2. Lead compounds from each group obtained by interaction profile analysis of 1604 potential TcTIM inhibitors.
GroupCompounds in the GroupID Best Composite of Each GroupDocking Score (Kcal/mol)Structure
1167BP1 (ZINC000150134991)−10.2Ijms 23 10047 i011
2135BP2 (ZINC000040071949)−10.4Ijms 23 10047 i012
3206BP3 (ZINC000030009841)−10.3Ijms 23 10047 i013
4171BP4 (ZINC000150010278)−10.0Ijms 23 10047 i014
5142BP5 (ZINC000040013445)−10.2Ijms 23 10047 i015
6103BP6 (ZINC000140212311)−10.2Ijms 23 10047 i016
7181BP7 (ZINC000040170214)−10.5Ijms 23 10047 i017
8162BP8 (ZINC000040170215)−10.6Ijms 23 10047 i018
9191BP9 (ZINC000170014382)−10.4Ijms 23 10047 i019
10146BP10 (ZINC000170072839)−10.4Ijms 23 10047 i020
Table
Table 3. Trypanocidal activity of benzimidazole derivatives against NINOA and INC-5 strains of T. cruzi.
Table 3. Trypanocidal activity of benzimidazole derivatives against NINOA and INC-5 strains of T. cruzi.
CompoundLC50 (µM) *
NINOAINC-5
BP2Ijms 23 10047 i021155.86 ± 3.4226.30 ± 15.4
BP5Ijms 23 10047 i022179.55 ± 19.7179.71 ± 19.0
Nfx70.41 ± 8.0139.37 ± 3.0
Bzn130.72 ± 8.8191.28 ± 1.8
* The results are means of three experiments.
Table
Table 4. Free energy of binding (Kcal/mol) of the Ligand-TcTIM and Ligand-HsTIM complexes.
Table 4. Free energy of binding (Kcal/mol) of the Ligand-TcTIM and Ligand-HsTIM complexes.
ComplexesΔEvdw (kcal/mol)ΔEele (kcal/mol)ΔGpolar (kcal/mol)ΔGSA (kcal/mol)ΔGb (kcal/mol)
L1-TcTIM−33.97 ± 0.15−3.85 ± 0.0616.99 ± 0.12−3.89 ± 0.01−24.73 ± 0.15
BP2-TcTIM−29.73 ± 0.24−3.78 ± 0.1117.77 ± 0.20−3.28 ± 0.03−19.03 ± 0.16
BP5-TcTIM−46.58 ± 0.2−5.21 ± 0.130.61 ± 0.13−4.49 ± 0.01−25.68 ± 0.20
BTS-TcTIM−24.53 ± 0.23−14.93 ± 0.3728.28 ± 0.21−2.93 ± 0.01−14.12 ± 0.23
L1-HsTIM−32.19 ± 0.44−11.71 ± 0.3734.74 ± 0.67−3.56 ± 0.03−12.72 ± 0.25
BP2-HsTIM−30.68 ± 0.24−0.09 ± 0.1214.02 ± 0.2−2.98 ± 0.02−19.72 ± 0.21
BP5-HsTIM−23.04 ± 0.18−5.22 ± 0.1815.63 ± 0.83−2.26 ± 0.02−14.89 ± 0.84
ΔEvdw, van der Waals contributions; ΔEele, Electrostatic contributions; ΔGpolar, Polar contributions; ΔGSA, SASA contributions, ΔGb; affinity energy.
Table
Table 5. Molecular physicochemical properties of BP2 and BP5 using SwissADME website and the ProTox-II server.
Table 5. Molecular physicochemical properties of BP2 and BP5 using SwissADME website and the ProTox-II server.
Physicochemical Properties
CompoundMW (g/mol) < 500Nb < 10Nhba < 10Nhbd < 5TPSA (Å2) < 140Log P < 5Log S 1
BP2459.5474298.402.70Moderately soluble
BP5427.5034188.883.20Moderately soluble
Probability of pharmacokinetic properties
CompoundBlood-brain permeabilityHuman gastrointestinal absorption 2P-glycoprotein substrateCYP1A2 inhibitorCYP2C19 inhibitorCYP2C9 inhibitorCYP2D6 inhibitorCYP3A4 inhibitorHepatotoxicity
BP2NoSolubleYesNoYesYesYesYesInactive
BP5NoSolubleYesNoYesYesYesYesInactive
MW: molecular weight, Nb: number of rotatable bonds, Nhba: number of hydrogen bond acceptors, Nhbd: number of hydrogen bond donors, TPSA: polar surface area, Log P: partition coefficient, Log S: solubility coefficient. null.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vázquez-Jiménez, L.K.; Juárez-Saldivar, A.; Gómez-Escobedo, R.; Delgado-Maldonado, T.; Méndez-Álvarez, D.; Palos, I.; Bandyopadhyay, D.; Gaona-Lopez, C.; Ortiz-Pérez, E.; Nogueda-Torres, B.; et al. Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents. Int. J. Mol. Sci. 2022, 23, 10047. https://doi.org/10.3390/ijms231710047

AMA Style

Vázquez-Jiménez LK, Juárez-Saldivar A, Gómez-Escobedo R, Delgado-Maldonado T, Méndez-Álvarez D, Palos I, Bandyopadhyay D, Gaona-Lopez C, Ortiz-Pérez E, Nogueda-Torres B, et al. Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents. International Journal of Molecular Sciences. 2022; 23(17):10047. https://doi.org/10.3390/ijms231710047

Chicago/Turabian Style

Vázquez-Jiménez, Lenci K., Alfredo Juárez-Saldivar, Rogelio Gómez-Escobedo, Timoteo Delgado-Maldonado, Domingo Méndez-Álvarez, Isidro Palos, Debasish Bandyopadhyay, Carlos Gaona-Lopez, Eyra Ortiz-Pérez, Benjamín Nogueda-Torres, and et al. 2022. "Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents" International Journal of Molecular Sciences 23, no. 17: 10047. https://doi.org/10.3390/ijms231710047

APA Style

Vázquez-Jiménez, L. K., Juárez-Saldivar, A., Gómez-Escobedo, R., Delgado-Maldonado, T., Méndez-Álvarez, D., Palos, I., Bandyopadhyay, D., Gaona-Lopez, C., Ortiz-Pérez, E., Nogueda-Torres, B., Ramírez-Moreno, E., & Rivera, G. (2022). Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents. International Journal of Molecular Sciences, 23(17), 10047. https://doi.org/10.3390/ijms231710047

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop