Structural considerations for the rational design of selective anti-trypanosomal agents: the role of the aromatic clusters at the interface of triosephosphate isomerase dimer

Preclinical Studies and Drug Combination of Low-Cost Molecules for Chagas Disease

Chagas disease is caused by the protozoan Trypanosoma cruzi (T. cruzi). It remains the major parasitic disease in Latin America and is spreading worldwide, affecting over 10 million people. Hundreds of new compounds with trypanosomicidal action have been identified from different sources such as synthetic or natural molecules, but they have been deficient in several stages of drug development (toxicology, scaling-up, and pharmacokinetics). Previously, we described a series of compounds with simple structures, low cost, and environmentally friendly production with potent trypanosomicidal activity in vitro and in vivo. These molecules are from three different families: thiazolidenehydrazines, diarylideneketones, and steroids. From this collection, we explored their capacity to inhibit the triosephosphate isomerase and cruzipain of T. cruzi. Then, the mechanism of action was explored using NMR metabolomics and computational molecular dynamics. Moreover, the mechanism of death was studied by flow cytometry. Consequently, five compounds, 314, 793, 1018, 1019, and 1260, were pre-clinically studied and their pharmacologic profiles indicated low unspecific toxicity. Interestingly, synergetic effects of diarylideneketones 793 plus 1018 and 793 plus 1019 were evidenced in vitro and in vivo. In vivo, the combination of compounds 793 plus 1018 induced a reduction of more than 90% of the peak of parasitemia in the acute murine model of Chagas disease.

Recent Advances in the Development of Triose Phosphate Isomerase Inhibitors as Antiprotozoal Agents

BACKGROUND

Parasitic diseases caused by protozoa such as Chagas disease, leishmaniasis, malaria, African trypanosomiasis, amebiasis, trichomoniasis, and giardiasis are considered serious public health problems in developing countries. Drug-resistance among parasites justifies the search for new therapeutic drugs and the identification of new targets becomes a valuable approach. In this scenario, glycolysis pathway which consists of the conversion of glucose into pyruvate plays an important role in the protozoa energy supply and it is therefore considered as a promising target. In this pathway, triose phosphate isomerase (TIM) plays an essential role in efficient energy production. Furthermore, protozoa TIM show structural differences with human enzyme counterparts suggesting the possibility of obtaining selective inhibitors. Therefore, TIM is considered a valid approach to develop new antiprotozoal agents, inhibiting the glycolysis in the parasite.

OBJECTIVE

In this review, we discuss the drug design strategies, structure-activity relationship, and binding modes of outstanding TIM inhibitors against Trypanosoma cruzi, Trypanosoma brucei, Plasmodium falciparum, Giardia lamblia, Leishmania mexicana, Trichomonas vaginalis, and Entamoeba histolytica.

RESULTS

TIM inhibitors showed mainly aromatic systems and symmetrical structure, where the size and type of heteroatom are important for enzyme inhibition. This inhibition is mainly based on the interaction with i) the interfacial region of TIM inducing changes on the quaternary and tertiary structure or ii) with the TIM catalytic region were the main pathways that disabled the catalytic activity of the enzyme.

CONCLUSION

Benzothiazole, benzoxazole, benzimidazole, and sulfhydryl derivatives stand out as TIM inhibitors. In silico and in vitro studies demonstrate that the inhibitors bind mainly at the TIM dimer interface. In this review, the development of new TIM inhibitors as antiprotozoal drugs is demonstrated as an important pharmaceutical strategy that may lead to new therapies for these ancient parasitic diseases.

Phenylpyrrole fungicides act on triosephosphate isomerase to induce methylglyoxal stress and alter hybrid histidine kinase activity

Fludioxonil, a natural product of pyrrolnitrin, is a potent fungicide used on crops worldwide. Drug action requires the presence of a group III hybrid histidine kinase (HHK) and the high osmolarity glycerol (HOG) pathway. We have reported that the drug does not act directly on HHK, but triggers the conversion of the kinase to a phosphatase, which dephosphorylates Ypd1 to constitutively activate HOG signaling. Still, the direct drug target remains unknown and mode of action ill defined. Here, we heterologously expressed a group III HHK, dimorphism-regulating kinase 1 (Drk1) in Saccharomyces cerevisae to delineate fludioxonil’s target and action. We show that the drug interferes with triosephosphate isomerase (TPI) causing release of methylglyoxal (MG). MG activates the group III HHK and thus the HOG pathway. Drug action involved Drk1 cysteine 392, as a C392S substitution increased drug resistance in vivo. Drug sensitivity was reversed by dimedone treatment, indicating Drk1 responds in vivo to an aldehydic stress. Fludioxonil treatment triggered elevated cytosolic methylglyoxal. Likewise, methylglyoxal treatment of Drk1-expressing yeast phenocopied treatment with fludioxonil. Fludioxonil directly inhibited TPI and also caused it to release methylglyoxal in vitro. Thus, TPI is a drug target of the phenylpyrrole class of fungicides, inducing elevated MG which alters HHK activity, likely converting the kinase to a phosphatase that acts on Ypd1 to trigger HOG pathway activation and fungal cell death.

Species-Specific Inactivation of Triosephosphate Isomerase from Trypanosoma brucei: Kinetic and Molecular Dynamics Studies

Human African Trypanosomiasis (HAT), a disease that provokes 2184 new cases a year in Sub-Saharan Africa, is caused by Trypanosoma brucei. Current treatments are limited, highly toxic, and parasite strains resistant to them are emerging. Therefore, there is an urgency to find new drugs against HAT. In this context, T. brucei depends on glycolysis as the unique source for ATP supply; therefore, the enzyme triosephosphate isomerase (TIM) is an attractive target for drug design. In the present work, three new benzimidazole derivatives were found as TbTIM inactivators (compounds 1, 2 and 3) with an I₅₀ value of 84, 82 and 73 µM, respectively. Kinetic analyses indicated that the three molecules were selective when tested against human TIM (HsTIM) activity. Additionally, to study their binding mode in TbTIM, we performed a 100 ns molecular dynamics simulation of TbTIM-inactivator complexes. Simulations showed that the binding of compounds disturbs the structure of the protein, affecting the conformations of important domains such as loop 6 and loop 8. In addition, the physicochemical and drug-like parameters showed by the three compounds suggest a good oral absorption. In conclusion, these molecules will serve as a guide to design more potent inactivators that could be used to obtain new drugs against HAT.

Structure–activity relationship and binding mode studies for a series of diketo-acids as HIV integrase inhibitors by 3D-QSAR, molecular docking and molecular dynamics simulations

We explored the main factors affecting the activity of compounds by different statistical and computational methods.

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)'s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM's catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand's effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop's dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.

Blind Dockings of Benzothiazoles to Multiple Receptor Conformations of Triosephosphate Isomerase from Trypanosoma cruzi and Human

We aim to uncover the binding modes of benzothiazoles, which have been reported as specific inhibitors of triosephosphate isomerase from the parasite Trypanosoma cruzi (TcTIM), by performing blind dockings on both TcTIM and human TIM (hTIM). Detailed analysis of binding sites and specific interactions are carried out based on ensemble dockings to multiple receptor conformers obtained from molecular dynamics simulations. In TcTIM dimer dockings, the inhibitors preferentially bind to the tunnel-shaped cavity formed at the interface of the subunits, whereas non-inhibitors mostly choose other sites. In contrast, TcTIM monomer binding interface and hTIM dimer interface do not present a specific binding site for the inhibitors. These findings point to the importance of the tunnel and of the dimeric form for inhibition of TcTIM. Specific interactions of the inhibitors and their sulfonate-free derivatives with the receptor residues indicate the significance of sulfonate group for binding affinity and positioning on the TcTIM dimer interface. One of the inhibitors also binds to the active site, which may explain its relatively higher inhibition effect on hTIM.

The design and inhibitory profile of new benzimidazole derivatives against triosephosphate isomerase from Trypanosoma cruzi: a problem of residue motility

To develop a new set of compounds with inhibitory activity against the triosephosphate isomerase of Trypanosoma cruzi (TcTIM), a group of benzimidazole derivatives was studied using four different docking procedures. These docking procedures differ in the number and type of mobile residues considered in the analysis. As a result of this methodology, a clustered analysis of plausible candidate structures was produced. A different set of previously synthesized compounds was used to validate this analysis. The validation showed that the best results correspond to the docking procedure in which the residues near the hydrophobic pocket of the protein's interface were considered mobile. A binding site for the best candidates was identified. Residues Tyr103, Glu105 and Lys113, among others, are important for the binding of this kind of compound. Residue Tyr103 is different in the human TIM, thus establishing a key feature for the future design of selective inhibitors.

Massive screening yields novel and selective Trypanosoma cruzi triosephosphate isomerase dimer-interface-irreversible inhibitors with anti-trypanosomal activity

Triosephosphate isomerase from Trypanosoma cruzi (TcTIM), an enzyme in the glycolytic pathway that exhibits high catalytic rates of glyceraldehyde-3-phosphate- and dihydroxyacetone-phosphate-isomerization only in its dimeric form, was screened against an in-house chemical library containing nearly 230 compounds belonging to different chemotypes. After secondary screening, twenty-six compounds from eight different chemotypes were identified as screening positives. Four compounds displayed selectivity for TcTIM over TIM from Homo sapiens and, concomitantly, in vitro activity against T. cruzi.

Toward a rational design of selective multi-trypanosomatid inhibitors: a computational docking study

Compound V7, a benzothiazole which was recently found as selective inhibitor of trypanosomal TIMs, was docked into TIMs from Trypanosoma cruzi, Trypanosoma brucei, Entamoeba histolytica, Plasmodium falciparum, yeast, and human. Structural analyses revealed the importance of the accessibility to the two aromatic clusters located at the dimer's interface for the selective inhibition of trypanosomal TIMs. Thus, it was found that different accessibilities of the protein interface of TIMs plays an important role in the inhibitory activity of benzothiazoles. These findings will contribute to the rational development and improvement of benzothiazoles to be used as multi-trypanosomatid inhibitors.

An intersubunit lock-and-key 'clasp' motif in the dimer interface of Delta class glutathione transferase

Structural investigations of a GST (glutathione transferase), adGSTD4-4, from the malaria vector Anopheles dirus show a novel lock-and-key 'Clasp' motif in the dimer interface of the Delta class enzyme. This motif also appears to be highly conserved across several insect GST classes, but differs from a previously reported mammalian lock-and-key motif. The aromatic 'key' residue not only inserts into a hydrophobic pocket, the 'lock', of the neighbouring subunit, but also acts as part of the 'lock' for the other subunit 'key'. The 'key' residues from both subunits show aromatic ring stacking with each other in a pi-pi interaction, generating a 'Clasp' in the middle of the subunit interface. Enzyme catalytic and structural characterizations revealed that single amino acid replacements in this 'Clasp' motif impacted on catalytic efficiencies, substrate selectivity and stability. Substitutions to the 'key' residue create strong positive co-operativity for glutathione binding, with a Hill coefficient approaching 2. The lock-and-key motif in general and especially the 'Clasp' motif with the pi-pi interaction appear to play a pivotal role in subunit communication between active sites, as well as in stabilizing the quaternary structure. Evidence of allosteric effects suggests an important role for this particular intersubunit architecture in regulating catalytic activity through conformational transitions of subunits. The observation of co-operativity in the mutants also implies that glutathione ligand binding and dimerization are linked. Quaternary structural changes of all mutants suggest that subunit assembly or dimerization basically manipulates subunit communication.