Structural Basis for Selective Inhibition of Trypanosomatid Glyceraldehyde-3-Phosphate Dehydrogenase: Molecular Docking and 3D QSAR Studies

Metallodrugs for the Treatment of Trypanosomatid Diseases: Recent Advances and New Insights

Trypanosomatid parasites are responsible for many Neglected Tropical Diseases (NTDs). NTDs are a group of illnesses that prevail in low-income populations, such as in tropical and subtropical areas of Africa, Asia, and the Americas. The three major human diseases caused by trypanosomatids are African trypanosomiasis, Chagas disease and leishmaniasis. There are known drugs for the treatment of these diseases that are used extensively and are affordable; however, the use of these medicines is limited by several drawbacks such as the development of chemo-resistance, side effects such as cardiotoxicity, low selectivity, and others. Therefore, there is a need to develop new chemotherapeutic against these tropical parasitic diseases. Metal-based drugs against NTDs have been discussed over the years as alternative ways to overcome the difficulties presented by approved antiparasitic agents. The study of late transition metal-based drugs as chemotherapeutics is an exciting research field in chemistry, biology, and medicine due to the ability to develop multitarget antiparasitic agents. The evaluation of the late transition metal complexes for the treatment of trypanosomatid diseases is provided here, as well as some insights about their mechanism of action.

Design of Inhibitors for Glyceraldehyde-3-phosphate Dehydrogenase (GAPDH) Enzyme of <i>Leishmania mexicana</i>

BACKGROUND

Leishmaniosis is a neglected tropical disease and glyceraldehyde 3- phosphate dehydrogenase (GAPDH) is a key enzyme in the design of new drugs to fight this disease.

OBJECTIVE

The present study aimed to evaluate potential inhibitors of GAPDH enzyme found in Leishmania mexicana (L. mexicana).

METHODS

A search for novel antileishmanial molecules was carried out based on similarities from the pharmacophoric point of view related to the binding site of the crystallographic enzyme using the ZINCPharmer server. The molecules selected in this screening were subjected to molecular docking and molecular dynamics simulations.

RESULTS

Consensual analysis of the docking energy values was performed, resulting in the selection of ten compounds. These ligand-receptor complexes were visually inspected in order to analyze the main interactions and subjected to toxicophoric evaluation, culminating in the selection of three compounds, which were subsequently submitted to molecular dynamics simulations. The docking results showed that the selected compounds interacted with GAPDH from L. mexicana, especially by hydrogen bonds with Cys166, Arg249, His194, Thr167, and Thr226. From the results obtained from molecular dynamics, it was observed that one of the loop regions, corresponding to the residues 195-222, can be related to the fitting of the substrate at the binding site, assisting in the positioning and the molecular recognition via residues responsible for the catalytic activity.

CONCLUSION

The use of molecular modeling techniques enabled the identification of promising compounds as inhibitors of the GAPDH enzyme from L. mexicana, and the results obtained here can serve as a starting point to design new and more effective compounds than those currently available.

2D and 3D QSAR Studies on a Series of Antichagasic Fenarimol Derivatives

Chagas disease is one of the most important neglected tropical diseases. Endemic in Latin America, the disease is a global public health problem, affecting several countries in North America, Europe, Asia and Oceania. The disease affects around 8-10 million people worldwide and the limited treatments available present low efficacy and severe side effects, highlighting the urgent need for new therapeutic options. In this work, the authors developed QSAR models for a series of fenarimol derivatives exhibiting anti-T. cruzi activity. The models were constructed using the Hologram QSAR (HQSAR), Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) methods. The QSAR models presented substantial predictive ability for a series of test set compounds (HQSAR, r2pred = 0.66; CoMFA, r2pred = 0.82; and CoMSIA, r2pred = 0.76), and were valuable to identify key structural features related to the observed trypanocidal activity. The results reported herein are useful for the design of novel derivatives having improved antichagasic properties.

A pharmacophore-based virtual screening approach for the discovery of Trypanosoma cruzi GAPDH inhibitors

BACKGROUND

Chagas disease is a major cause of morbidity and death for millions of people in Latin America. The drugs currently available exhibit poor efficacy and severe side effects. Therefore, there is an urgent need for new, safe and effective drugs against Chagas disease. The vital dependence on glycolysis as energy source makes the glycolytic enzymes of Trypanosoma cruzi, the causative agent of Chagas disease, attractive targets for drug design. In this work, glyceraldehyde-3-phosphate dehydrogenase from T. cruzi (TcGAPDH) was employed as molecular target for the discovery of new inhibitors as hits.

RESULTS

Integrated protein-based pharmacophore and structure-based virtual screening approaches resulted in the identification of three hits from three chemical classes with moderate inhibitory activity against TcGAPDH. The inhibitors showed IC50 values in the high micromolar range.

CONCLUSION

The new chemotypes are attractive molecules for future medicinal chemistry efforts aimed at developing new lead compounds for Chagas disease.

Pattern recognition techniques applied to the study of leishmanial glyceraldehyde-3-phosphate dehydrogenase inhibition

Chemometric pattern recognition techniques were employed in order to obtain Structure-Activity Relationship (SAR) models relating the structures of a series of adenosine compounds to the affinity for glyceraldehyde 3-phosphate dehydrogenase of Leishmania mexicana (LmGAPDH). A training set of 49 compounds was used to build the models and the best ones were obtained with one geometrical and four electronic descriptors. Classification models were externally validated by predictions for a test set of 14 compounds not used in the model building process. Results of good quality were obtained, as verified by the correct classifications achieved. Moreover, the results are in good agreement with previous SAR studies on these molecules, to such an extent that we can suggest that these findings may help in further investigations on ligands of LmGAPDH capable of improving treatment of leishmaniasis.

Interactions of antiparasitic alkaloids with Leishmania protein targets: a molecular docking analysis

Background: Leishmaniasis is a collection of chronic diseases caused by protozoa of the genus Leishmania. Current antileishmanial chemotherapeutics have demonstrated adverse side effects and therefore R&D into new safer alternative treatments are needed. Methods: A molecular docking analysis has been carried out to assess possible Leishmania biochemical targets of antiparasitic alkaloids. A total of 209 antiparasitic alkaloids were docked with 24 Leishmania protein targets. Results: The strongest docking alkaloid ligands were flinderoles A and B and juliflorine with Leishmania major methionyl-tRNA synthetase; juliflorine, juliprosine, prosopilosidine and prosopilosine with Leishmania mexicana glycerol-3-phosphate dehydrogenase; and ancistrogriffithine A with L. major N-myristoyl transferase. Conclusion: This molecular docking study has provided evidence for what classes and structural types of alkaloids may be targeting specific Leishmania protein targets.

In-silico Leishmania target selectivity of antiparasitic terpenoids

Neglected Tropical Diseases (NTDs), like leishmaniasis, are major causes of mortality in resource-limited countries. The mortality associated with these diseases is largely due to fragile healthcare systems, lack of access to medicines, and resistance by the parasites to the few available drugs. Many antiparasitic plant-derived isoprenoids have been reported, and many of them have good in vitro activity against various forms of Leishmania spp. In this work, potential Leishmania biochemical targets of antiparasitic isoprenoids were studied in silico. Antiparasitic monoterpenoids selectively docked to L. infantum nicotinamidase, L. major uridine diphosphate-glucose pyrophosphorylase and methionyl t-RNA synthetase. The two protein targets selectively targeted by germacranolide sesquiterpenoids were L. major methionyl t-RNA synthetase and dihydroorotate dehydrogenase. Diterpenoids generally favored docking to L. mexicana glycerol-3-phosphate dehydrogenase. Limonoids also showed some selectivity for L. mexicana glycerol-3-phosphate dehydrogenase and L. major dihydroorotate dehydrogenase while withanolides docked more selectively with L. major uridine diphosphate-glucose pyrophosphorylase. The selectivity of the different classes of antiparasitic compounds for the protein targets considered in this work can be explored in fragment- and/or structure-based drug design towards the development of leads for new antileishmanial drugs.

Identification of electronic and structural descriptors of adenosine analogues related to inhibition of leishmanial glyceraldehyde-3-phosphate dehydrogenase

Quantitative structure-activity relationship (QSAR) studies were performed in order to identify molecular features responsible for the antileishmanial activity of 61 adenosine analogues acting as inhibitors of the enzyme glyceraldehyde 3-phosphate dehydrogenase of Leishmania mexicana (LmGAPDH). Density functional theory (DFT) was employed to calculate quantum-chemical descriptors, while several structural descriptors were generated with Dragon 5.4. Variable selection was undertaken with the ordered predictor selection (OPS) algorithm, which provided a set with the most relevant descriptors to perform PLS, PCR and MLR regressions. Reliable and predictive models were obtained, as attested by their high correlation coefficients, as well as the agreement between predicted and experimental values for an external test set. Additional validation procedures were carried out, demonstrating that robust models were developed, providing helpful tools for the optimization of the antileishmanial activity of adenosine compounds.

Structure- and ligand-based drug design approaches for neglected tropical diseases

Drug discovery has moved toward more rational strategies based on our increasing understanding of the fundamental principles of protein–ligand interactions. Structure- (SBDD) and ligand-based drug design (LBDD) approaches bring together the most powerful concepts in modern chemistry and biology, linking medicinal chemistry with structural biology. The definition and assessment of both chemical and biological space have revitalized the importance of exploring the intrinsic complementary nature of experimental and computational methods in drug design. Major challenges in this field include the identification of promising hits and the development of high-quality leads for further development into clinical candidates. It becomes particularly important in the case of neglected tropical diseases (NTDs) that affect disproportionately poor people living in rural and remote regions worldwide, and for which there is an insufficient number of new chemical entities being evaluated owing to the lack of innovation and R&D investment by the pharmaceutical industry. This perspective paper outlines the utility and applications of SBDD and LBDD approaches for the identification and design of new small-molecule agents for NTDs.

Drug discovery and the use of computational approaches for infectious diseases

For centuries infectious diseases were the scourge of humanity, overcome only by the discovery of vaccination and penicillin. With an armamentarium of effective antibiotics, vaccines and drugs at hand, infectious diseases for many years were considered to be negligible. With the onset of the AIDS pandemic, the return of tuberculosis and influenza (e.g., swine influenza) this notion has changed in recent years. Drug discovery for infectious diseases, therefore, is again gaining increasing interest. This article discusses the drug-discovery process in this area and introduces major computational approaches used to identify suitable drug targets and to discover and optimize chemical lead compounds towards drug candidates using examples from antiparasitic drug discovery.

Discovery of new inhibitors of Schistosoma mansoni PNP by pharmacophore-based virtual screening

Schistosomiasis is considered the second most important tropical parasitic disease, with severe socioeconomic consequences for millions of people worldwide. Schistosoma mansoni , one of the causative agents of human schistosomiasis, is unable to synthesize purine nucleotides de novo, which makes the enzymes of the purine salvage pathway important targets for antischistosomal drug development. In the present work, we describe the development of a pharmacophore model for ligands of S. mansoni purine nucleoside phosphorylase (SmPNP) as well as a pharmacophore-based virtual screening approach, which resulted in the identification of three thioxothiazolidinones (1-3) with substantial in vitro inhibitory activity against SmPNP. Synthesis, biochemical evaluation, and structure-activity relationship investigations led to the successful development of a small set of thioxothiazolidinone derivatives harboring a novel chemical scaffold as new competitive inhibitors of SmPNP at the low-micromolar range. Seven compounds were identified with IC(50) values below 100 μM. The most potent inhibitors 7, 10, and 17 with IC(50) of 2, 18, and 38 μM, respectively, could represent new potential lead compounds for further development of the therapy of schistosomiasis.

New drugs for neglected infectious diseases: Chagas' disease

Chagas' disease (CD) is caused by the protozoan Trypanosoma cruzi (Tc) and remains an important cause of morbidity and mortality. Most researchers in the field now agree that chronic low grade parasite persistence in tissue drives tissue damage and the autoimmune component of CD. Current therapy relies on two compounds: benznidazole and nifurtimox. Despite their long history in the treatment of CD, both compounds induce significant side-effects. In the current issue of the BJP, two contributions demonstrate that NO-donors are active, especially in combination with benznidazole, against Tc in vitro and in experimental models in vivo. The basic concept used by the authors to develop novel anti-Tc compounds relied on the demonstrated ability of nitric oxide to kill the parasite. There are several issues still to be resolved but the reported studies are a clear advance to the field and should be considered for further pre-clinical development.

Novel ruthenium complexes as potential drugs for Chagas's disease: enzyme inhibition and in vitro/in vivo trypanocidal activity

BACKGROUND AND PURPOSE

The discovery of the pharmacological functions of nitric oxide has led to the development of NO donor compounds as therapeutic agents. A new generation of ruthenium NO donors, cis-[Ru(NO)(bpy)(2)L]X(n), has been developed, and our aim was to show that these complexes are able to lyse Trypanosoma cruzi in vitro and in vivo.

EXPERIMENTAL APPROACH

NO donors were incubated with T. cruzi and their anti-T. cruzi activities evaluated as the percentage of lysed parasites compared to the negative control. In vivo, trypanocidal activity was evaluated by observing the levels of parasitaemia, survival rate and elimination of amastigotes in mouse myocardial tissue. The inhibition of GAPDH was monitored by the biochemical reduction of NAD(+) to NADH.

KEY RESULTS

The NO donors cis-[Ru(NO)(bpy)(2)L]X(n) presented inhibitory effects on T. cruzi GAPDH (IC(50) ranging from 89 to 153 microM). The crystal structure of the enzyme shows that the inhibitory mechanism is compatible with S-nitrosylation of the active cysteine (cys166) site. Compounds cis-[Ru(NO)(bpy)(2)imN](PF(6))(3) and cis-[Ru(NO)(bpy)(2)SO(3)]PF(6), at a dose of 385 nmol.kg(-1), yielded survival rates of 80 and 60%, respectively, in infected mice, and eradicated any amastigotes from their myocardial tissue.

CONCLUSIONS AND IMPLICATIONS

The ruthenium compounds exhibited potent in vitro and in vivo trypanocidal activities at doses up to 1000-fold lower than the clinical dose for benznidazole. Furthermore, one mechanism of action of these compounds is via the S-nitrosylation of Cys166 of T. cruzi GAPDH. Thus, these compounds show huge potential as candidates for the development of new drugs for the treatment of Chagas's disease.

Pharmacophore-based 3D QSAR studies on a series of high affinity 5-HT1A receptor ligands

5-HT(1A) receptor antagonists have been employed to treat depression, but the lack of structural information on this receptor hampers the design of specific and selective ligands. In this study, we have performed CoMFA studies on a training set of arylpiperazines (high affinity 5-HT(1A) receptor ligands) and to produce an effective alignment of the data set, a pharmacophore model was produced using Galahad. A statistically significant model was obtained, indicating a good internal consistency and predictive ability for untested compounds. The information gathered from our receptor-independent pharmacophore hypothesis is in good agreement with results from independent studies using different approaches. Therefore, this work provides important insights on the chemical and structural basis involved in the molecular recognition of these compounds.

Genetic validation of aldolase and glyceraldehyde-3-phosphate dehydrogenase as drug targets in Trypanosoma brucei

Aldolase (ALD) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of Trypanosoma brucei are considered to be promising targets for chemotherapeutic treatment of African sleeping sickness, because glycolysis is the single source of ATP for the parasite when living in the human bloodstream. Moreover, these enzymes appeared to possess distinct kinetic and structural properties that have already been exploited for the discovery of effective and selective inhibitors with trypanocidal activity. Here we present an experimental, quantitative assessment of the importance of these enzymes for the glycolytic pathway. This was achieved by decreasing the concentrations of ALD and GAPDH by RNA interference. The effects of these knockdowns on parasite growth, levels of various enzymes and transcripts, enzyme activities and glucose consumption were studied. A partial depletion of ALD and GAPDH was already sufficient to rapidly kill the trypanosomes. An effect was also observed on the activity of some other glycolytic enzymes.

Structure-based approach for the study of estrogen receptor binding affinity and subtype selectivity

Estrogens exert important physiological effects through the modulation of two human estrogen receptor (hER) subtypes, alpha (hERalpha) and beta (hERbeta). Because the levels and relative proportion of hERalpha and hERbeta differ significantly in different target cells, selective hER ligands could target specific tissues or pathways regulated by one receptor subtype without affecting the other. To understand the structural and chemical basis by which small molecule modulators are able to discriminate between the two subtypes, we have applied three-dimensional target-based approaches employing a series of potent hER-ligands. Comparative molecular field analysis (CoMFA) studies were applied to a data set of 81 hER modulators, for which binding affinity values were collected for both hERalpha and hERbeta. Significant statistical coefficients were obtained (hERalpha, q(2) = 0.76; hERbeta, q(2) = 0.70), indicating the internal consistency of the models. The generated models were validated using external test sets, and the predicted values were in good agreement with the experimental results. Five hER crystal structures were used in GRID/PCA investigations to generate molecular interaction fields (MIF) maps. hERalpha and hERbeta were separated using one factor. The resulting 3D information was integrated with the aim of revealing the most relevant structural features involved in hER subtype selectivity. The final QSAR and GRID/PCA models and the information gathered from 3D contour maps should be useful for the design of novel hER modulators with improved selectivity.

Quantitative structure-activity relationships for a series of inhibitors of cruzain from Trypanosoma cruzi: molecular modeling, CoMFA and CoMSIA studies

Human parasitic diseases are the foremost threat to human health and welfare around the world. Trypanosomiasis is a very serious infectious disease against which the currently available drugs are limited and not effective. Therefore, there is an urgent need for new chemotherapeutic agents. One attractive drug target is the major cysteine protease from Trypanosoma cruzi, cruzain. In the present work, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) studies were conducted on a series of thiosemicarbazone and semicarbazone derivatives as inhibitors of cruzain. Molecular modeling studies were performed in order to identify the preferred binding mode of the inhibitors into the enzyme active site, and to generate structural alignments for the three-dimensional quantitative structure-activity relationship (3D QSAR) investigations. Statistically significant models were obtained (CoMFA, r2=0.96 and q2=0.78; CoMSIA, r2=0.91 and q2=0.73), indicating their predictive ability for untested compounds. The models were externally validated employing a test set, and the predicted values were in good agreement with the experimental results. The final QSAR models and the information gathered from the 3D CoMFA and CoMSIA contour maps provided important insights into the chemical and structural basis involved in the molecular recognition process of this family of cruzain inhibitors, and should be useful for the design of new structurally related analogs with improved potency.

Fragment-based QSAR: perspectives in drug design

Drug design is a process driven by innovation and technological breakthroughs involving a combination of advanced experimental and computational methods. A broad variety of medicinal chemistry approaches can be used for the identification of hits, generation of leads, as well as to accelerate the optimization of leads into drug candidates. Quantitative structure-activity relationship (QSAR) methods are among the most important strategies that can be applied for the successful design of small molecule modulators having clinical utility. Hologram QSAR (HQSAR) is a modern 2D fragment-based QSAR method that employs specialized molecular fingerprints. HQSAR can be applied to large data sets of compounds, as well as traditional-size sets, being a versatile tool in drug design. The HQSAR approach has evolved from a classical use in the generation of standard QSAR models for data correlation and prediction into advanced drug design tools for virtual screening and pharmacokinetic property prediction. This paper provides a brief perspective on the evolution and current status of HQSAR, highlighting present challenges and new opportunities in drug design.