Assaying phenothiazine derivatives as trypanothione reductase and glutathione reductase inhibitors by theoretical docking and Molecular Dynamics studies

Analysis of the interaction of antimalarial agents with Plasmodium falciparum glutathione reductase through molecular mechanical calculations

CONTEXT

Malaria remains a significant global health challenge with emerging resistance to current treatments. Plasmodium falciparum glutathione reductase (PfGR) plays a critical role in the defense mechanisms of malaria parasites against oxidative stress. In this study, we investigate the potential of targeting PfGR with conventional antimalarials and dual drugs combining aminoquinoline derivatives with GR inhibitors, which reveal promising interactions between PfGR and studied drugs. The naphthoquinone Atovaquone demonstrated particularly high affinity and potential dual-mode binding with the enzyme active site and cavity. Furthermore, dual drugs exhibit enhanced binding affinity, suggesting their efficacy in inhibiting PfGR, where the aliphatic ester bond (linker) is essential for effective binding with the enzyme's active site. Overall, this research provides important insights into the interactions between antimalarial agents and PfGR and encourages further exploration of its role in the mechanisms of action of antimalarials, including dual drugs, to enhance antiparasitic efficacy.

METHODS

The drugs were tested as PfGR potential inhibitors via molecular docking on AutoDock 4, which was performed based on the preoptimized structures in HF/3-21G-PCM level of theory on ORCA 5. Drug-receptor systems with the most promising binding affinities were then studied with a molecular dynamic's simulation on AMBER 16. The molecular dynamics simulations were performed with a 100 ns NPT ensemble employing GAFF2 forcefield in the temperature of 310 K, integration time step of 2 fs, and non-bond cutoff distance of 6.0 Å.

Navigating drug repurposing for Chagas disease: advances, challenges, and opportunities

Chagas disease is a vector-borne illness caused by the protozoan parasite Trypanosoma cruzi (T. cruzi). It poses a significant public health burden, particularly in the poorest regions of Latin America. Currently, there is no available vaccine, and chemotherapy has been the traditional treatment for Chagas disease. However, the treatment options are limited to just two outdated medicines, nifurtimox and benznidazole, which have serious side effects and low efficacy, especially during the chronic phase of the disease. Collectively, this has led the World Health Organization to classify it as a neglected disease. To address this problem, new drug regimens are urgently needed. Drug repurposing, which involves the use of existing drugs already approved for the treatment of other diseases, represents an increasingly important option. This approach offers potential cost reduction in new drug discovery processes and can address pharmaceutical bottlenecks in the development of drugs for Chagas disease. In this review, we discuss the state-of-the-art of drug repurposing approaches, including combination therapy with existing drugs, to overcome the formidable challenges associated with treating Chagas disease. Organized by original therapeutic area, we describe significant recent advances, as well as the challenges in this field. In particular, we identify candidates that exhibit potential for heightened efficacy and reduced toxicity profiles with the ultimate objective of accelerating the development of new, safe, and effective treatments for Chagas disease.

Identification of novel thiadiazin derivatives as potentially selective inhibitors towards trypanothione reductase from Trypanosoma cruzi by molecular docking using the numerical index poses ratio Pr and the binding mode analysis

Abstract

Chagas disease is a serious health problem in Central and South America for which effective treatment is not currently available. This illness is caused by the protozoa Trypanosoma cruzi, a species that relies on a thiol-based metabolism to regulate oxidative stress. Trypanothione reductase enzyme plays a central role in the metabolic pathway of the parasite. In this work, a virtual screening of a library of novel thiadiazine derivatives against trypanothione reductase using molecular docking was performed. Four different series of hybrid ligands having in the structure one or two peptoid moieties (series I and II) or the tetrazole ring (series III and IV) were considered. An ad hoc numerical index called poses ratio was introduced to interpret the results of the docking analysis and to establish relevant structure-interaction relationships. In addition, six binding modes were found for the ligands with the highest populated conformational clusters after applying contact-based analysis. The most regular and relevant were binding modes I and II, found mainly for ligands from series I. A subsequent molecular docking on human glutathione reductase enzyme allowed to assess the possible cytotoxicity of the ligands towards human cells. A selective binding profile was found for ligands with interactions in the Hydrophobic cleft, the spermidine and the Z subsites inside the active site of trypanothione reductase. At the end of the study, new thiadiazine-based compounds were identified as plausible candidates to selectively inhibit the parasitic enzyme.

Graphic abstract

New insight into the role of glutathione reductase in glutathione peroxidase-like activity determination by coupled reductase assay: Molecular Docking Study

Previously we have shown that among 15 substituted salicyloyl (2-hydroxybenzoyl) 5-seleninic acids (SSAs) 4 compounds with longer side chains or a cyclohexyl group exhibit no glutathione peroxidase (GPx)-like activity in the coupled reductase assay. Experimental inhibition of glutathione reductase (GR) by the selenenylsulfide (a main intermediate in the catalytic cycle for GPx-like activity determination) of one of the inactive compounds led us to assess the interactions between 15 selenenylsulfide compounds and the active site of GR by molecular docking. Docking results showed that S and Se atoms in selenenylsulfides of the compounds with no GPx-like activity were beyond 5 Å from S atom of Cys-58 or N atom of imidazole ring of His-467 (Root Mean Square Distances for general assessment of 3 major distances were over 4.8 Å) in the active site, so that they could not be catalyzed to be reduced by GR. Furthermore, their docking scores over 89 Kcal/mol meant that the selenenylsulfides were bound too strongly to the active site to leave it, leading eventually to inhibition of GR. We also applied the molecular docking to other GPx mimics such as ebselen, cyclic seleninate esters and di(propylaminomethylphenyl) diselenides to explain the differences in their GPx-like activity depending to the assays used. Our results suggest that the reduction of a selenenylsulfide by GR plays a positive role in GPx-like activity of GPx mimics in the coupled assay and recommended the prediction of possibility and strength of GPx-like activity by molecular docking before entering experimental research.

Relevance of Trypanothione Reductase Inhibitors on Trypanosoma cruzi Infection: A Systematic Review, Meta-Analysis, and In Silico Integrated Approach

Due to the rudimentary antioxidant defenses in Trypanosoma cruzi, disruptors of redox balance are promising candidates for new antitrypanosomal drugs. We developed an integrated model based on systematic review, meta-analyses, and molecular modeling to evaluate the effect of trypanothione reductase (TR) inhibitors in T. cruzi infections. Our findings indicated that the TR inhibitors analyzed were effective in reducing parasitemia and mortality due to Trypanosoma cruzi infection in animal models. The most investigated drugs (clomipramine and thioridazine) showed no beneficial effects on the occurrence of infection-related electrocardiographic abnormalities or the affinity and density of cardiac β-adrenergic receptors. The affinity between the tested ligands and the active site of TR was confirmed by molecular docking. However, the molecular affinity score was unable to explain TR inhibition and T. cruzi death in vitro or the antiparasitic potential of these drugs when tested in preclinical models of T. cruzi infection. The divergence of in silico, in vitro, and in vivo findings indicated that the anti-T. cruzi effects of the analyzed drugs were not restricted to TR inhibition. As in vivo studies on TR inhibitors are still scarce and exhibit methodological limitations, mechanistic and highly controlled studies are required to improve the quality of evidence.

Structural analysis and molecular docking of trypanocidal aryloxy-quinones in trypanothione and glutathione reductases: a comparison with biochemical data

A set of aryloxy-quinones, previously synthesized and evaluated against Trypanosoma cruzi epimastigotes cultures, were found more potent and selective than nifurtimox. One of the possible mechanisms of the trypanocidal activity of these quinones could be inhibition of trypanothione reductase (TR). Considering that glutathione reductase (GR) is the equivalent of TR in humans, biochemical, kinetic, and molecular docking studies in TR and GR were envisaged and compared with the trypanocidal and cytotoxic data of a set of aryloxy-quinones. Biochemical assays indicated that three naphthoquinones (Nq-h, Nq-g, and Nq-d) selectively inhibit TR and the TR kinetic analyses indicated that Nq-h inhibit TR in a noncompetitive mechanism. Molecular dockings were performed in TR and GR in the following three putative binding sites: the catalytic site, the dimer interface, and the nicotinamide adenine dinucleotide phosphate-binding site. In TR and GR, the aryloxy-quinones were found to exhibit high affinity for a site near it cognate-binding site in a place in which the noncompetitive kinetics could be justified. Taking as examples the three compounds with TR specificity (TRS) (Nq-h, Nq-g, and Nq-d), the presence of a network of contacts with the quinonic ring sustained by the triad of Lys62, Met400', Ser464' residues, seems to contribute hardly to the TRS. Compound Nq-b, a naphthoquinone with nitrophenoxy substituent, proved to be the best scaffold for the design of trypanocidal compounds with low toxicity. However, the compound displayed only a poor and non-selective effect toward TR indicating that TR inhibition is not the main reason for the antiparasitic activity of the aryloxy-quinones.

Molecular Docking and Binding Mode Analysis of Plant Alkaloids as in Vitro and in silico Inhibitors of Trypanothione Reductase from Trypanosoma cruzi

Trypanothione reductase (TryR) is a key enzyme in the metabolism of Trypanosoma cruzi, the parasite responsible for Chagas disease. The available repertoire of TryR inhibitors relies heavily on synthetic substrates of limited structural diversity, and less on plant-derived natural products. In this study, a molecular docking procedure using a Lamarckian Genetic Algorithm was implemented to examine the protein-ligand binding interactions of strong in vitro inhibitors for which no X-ray data is available. In addition, a small, skeletally diverse, set of natural alkaloids was assessed computationally against T. cruzi TryR in search of new scaffolds for lead development. The preferential binding mode (low number of clusters, high cluster population), together with the deduced binding interactions were used to discriminate among the virtual inhibitors. This study confirms the prior in vitro data and proposes quebrachamine, cephalotaxine, cryptolepine, (22 S,25 S)-tomatidine, (22 R,25 S)-solanidine, and (22 R,25 R)-solasodine as new alkaloid scaffold leads in the search for more potent and selective TryR inhibitors.

Antimony(III) complexes with 2-amino-4,6-dimethoxypyrimidines: Synthesis, characterization and biological evaluation

Novel pyrimidine compound bearing disulfide bridge, 5,5'-disulfanediylbis(2-amino-4,6-dimetoxypyrimidine) (3) was synthesized by reduction of 2-amino-4,6-dimethoxy-5-thiocyanatopyrimidine for the first time, and its structure was confirmed by X-ray crystallographic analysis. Novel binuclear antimony(III) compound of (3), {Sb[5,5'-disulfanediylbis(2-amino-4,6-dimetoxypyrimidine)]Cl3}2 (4) and mononuclear antimony(III) compounds, SbL2Cl3, [L: 2-amino-5-thiol-4,6-dimethoxy pyrimidine (2) and 2-amino-5-(1H-tetrazol-5-ylthio)-4,6-dimethoxypyrimidine (6)] were synthesized and characterized with the help of elemental analysis, molecular conductivity, FT-IR, (1)H-NMR and LC-MS techniques. The geometrical structures optimized by a DFT/B3LYP/LANL2DZ method of the compounds, indicated that monomeric compounds have square pyramidal shape. Both antileishmanial activity against Leishmania tropica promastigote and glutathione reductase inhibitory activity were determined in vitro. The results showed that (3) has the best biological activity.

Indazoles: a new top seed structure in the search of efficient drugs against Trypanosoma cruzi

For years, Chagas disease treatment has been limited to only two drugs of highly questionable and controversial use (Nifurtimox® and Benznidazole®). In the search of effective drugs, many efforts have been made, but only a few structures have emerged as actual candidates. Heading into this, the multitarget-directed approach appears as the best choice. In this framework, indazoles were shown to be potent Trypanosoma cruzi growth inhibitors, being able to lead both the formation of reactive oxygen species and the inhibition of trypanothione reductase. Herein, we discuss the main structural factors that rule the anti-T. cruzi properties of indazoles, and how they would be involved in the biological properties as well as in the action mechanisms, attempting to make parallels between the old paradigms and current evidences in order to outline what could be the next steps to follow in regard to the future drug design for Chagas disease treatment.

The receptor-dependent LQTA-QSAR: application to a set of trypanothione reductase inhibitors

A new Receptor-Dependent LQTA-QSAR approach, RD-LQTA-QSAR, is proposed as a new 4D-QSAR method. It is an evolution of receptor independent LQTA-QSAR. This approach uses the free GROMACS package to carry out molecular dynamics simulations and generates a conformational ensemble profile for each compound. Such an ensemble is used to build molecular interaction field-based QSAR models, as in CoMFA. To show the potential of this methodology, a set of 38 phenothiazine derivatives that are specific competitive T. cruzi trypanothione reductase inhibitors, was chosen. Using a combination of molecular docking and molecular dynamics simulations, the binding mode of the phenotiazine derivatives was evaluated in a simulated induced fit approach. The ligands alignments were performed using both ligand and binding site atoms, enabling unbiased alignment. The models obtained were extensively validated by leave-N-out cross-validation and y-randomization techniques to test for their robustness and absence of chance correlation. The final model presented Q(2) LOO of 0.87 and R² of 0.92 and a suitable external prediction of [Formula: see text]= 0.78. The adapted binding site obtained is useful to perform virtual screening and ligand structure-based design and the descriptors in the final model can aid in the design new inhibitors.