Conformational changes in Leishmania mexicana glyceraldehyde-3-phosphate dehydrogenase induced by designed inhibitors

Stoichiometry of ligand binding and role of C-terminal lysines in Mycobacterium tuberculosis and human GAPDH multifunctionality

Glyceraldehyde-3-phosphate-dehydrogenase (GAPDH; EC1.2.1.12) has several functions in Mycobacterium tuberculosis (Mtb) and the human host. Apart from its role in glycolysis, it serves both as a cell surface and a secreted receptor for plasmin(ogen) (Plg/Plm), transferrin (Tf), and lactoferrin (Lf). Plg sequestration by Mtb GAPDH facilitates bacterial adhesion and tissue invasion, while an equivalent interaction with host GAPDH regulates immune cell migration. In both, host and microbe, internalization of Tf/Lf-GAPDH complexes serves as a route for iron acquisition. To date, the structure of Mtb GAPDH or the residues involved in these moonlighting interactions have not been identified. This study provides the first known X-ray crystal structure of Mtb GAPDH. Through further mutagenesis and functional assays, we found that the C-terminal lysines of Mtb and human GAPDH affect enzyme activity and ligand binding. We also establish the stoichiometry of Plg, Tf and Lf interactions with the GAPDH tetramer. Lastly, molecular simulation studies reveal the interactions of the C-terminal lysine residues.

Experimental structure based drug design (SBDD) applications for anti-leishmanial drugs: A paradigm shift?

Leishmaniasis is a group of neglected tropical diseases caused by at least 20 species of Leishmania protozoa, which are spread by the bite of infected sandflies. There are three main forms of the disease: cutaneous leishmaniasis (CL, the most common), visceral leishmaniasis (VL, also known as kala-azar, the most serious), and mucocutaneous leishmaniasis. One billion people live in areas endemic to leishmaniasis, with an annual estimation of 30,000 new cases of VL and more than 1 million of CL. New treatments for leishmaniasis are an urgent need, as the existing ones are inefficient, toxic, and/or expensive. We have revised the experimental structure-based drug design (SBDD) efforts applied to the discovery of new drugs against leishmaniasis. We have grouped the explored targets according to the metabolic pathways they belong to, and the key achieved advances are highlighted and evaluated. In most cases, SBDD studies follow high-throughput screening campaigns and are secondary to pharmacokinetic optimization, due to the majoritarian belief that there are few validated targets for SBDD in leishmaniasis. However, some SBDD strategies have significantly contributed to new drug candidates against leishmaniasis and a bigger number holds promise for future development.

Structure of glyceraldehyde-3-phosphate dehydrogenase from Paracoccidioides lutzii in complex with an aldonic sugar acid

The pathogen Paracoccidioides lutzii (Pb01) is found in South America countries Colombia, Ecuador, Venezuela and Brazil, especially in the central, west, and north regions of the latter. It belongs to the Ajellomycetaceae family, Onygenales order, and is typically thermodimorphic, presenting yeast cells when it grows in animal tissues, but mycelia when in the environment, where it produces the infectious propagule. This fungus is one of the etiologic agents of Paracoccidioidomycosis (PCM), the most important endemic fungal infection in Latin America. Investigations on its genome have contributed to a better understanding about its metabolism and revealed the complexity of several metabolic glycolytic pathways. Glyceraldehyde-3-Phosphate Dehydrogenase from Paracoccidioides lutzii (PlGAPDH) is considered a moonlighting protein and participates in several biological processes of this pathogen. The enzyme was expressed and purified, as seen in SDS-PAGE gel, crystallized and had its three dimensional structure (3D) determined in complex with NAD+, a sulphate ion and d-galactonic acid, therefore, a type of 'GAA site'. It is the first GAPDH structure to show this chemical type in this site and how this protein can bind an acid derived from oxidation of a linear hexose.

Virtual Screening of Benzimidazole Derivatives as Potential Triose Phosphate Isomerase Inhibitors with Biological Activity against Leishmania mexicana

Leishmania mexicana (L. mexicana) is a causal agent of cutaneous leishmaniasis (CL), a “Neglected disease”, for which the search for new drugs is a priority. Benzimidazole is a scaffold used to develop antiparasitic drugs; therefore, it is interesting molecule against L. mexicana. In this work, a ligand-based virtual screening (LBVS) of the ZINC15 database was performed. Subsequently, molecular docking was used to predict the compounds with potential binding at the dimer interface of triosephosphate isomerase (TIM) of L. mexicana (LmTIM). Compounds were selected on binding patterns, cost, and commercial availability for in vitro assays against L. mexicana blood promastigotes. The compounds were analyzed by molecular dynamics simulation on LmTIM and its homologous human TIM. Finally, the physicochemical and pharmacokinetic properties were determined in silico. A total of 175 molecules with docking scores between −10.8 and −9.0 Kcal/mol were obtained. Compound E2 showed the best leishmanicidal activity (IC50 = 4.04 µM) with a value similar to the reference drug pentamidine (IC50 = 2.23 µM). Molecular dynamics analysis predicted low affinity for human TIM. Furthermore, the pharmacokinetic and toxicological properties of the compounds were suitable for developing new leishmanicidal agents.

3-Bromopyruvate: A new strategy for inhibition of glycolytic enzymes in Leishmania amazonensis

The compound 3-bromopyruvate (3-BrPA) is well-known and studies from several researchers have demonstrated its involvement in tumorigenesis. It is an analogue of pyruvic acid that inhibits ATP synthesis by inhibiting enzymes from the glycolytic pathway and oxidative phosphorylation. In this work, we investigated the effect of 3-BrPA on energy metabolism of L. amazonensis. In order to verify the effect of 3-BrPA on L. amazonensis glycolysis, we measured the activity level of three glycolytic enzymes located at different points of the pathway: (i) glucose kinases, step 1, (ii) glyceraldehyde 3-phosphate dehydrogenase (GAPDH), step 6, and (iii) enolase, step 9. 3-BrPA, in a dose-dependent manner, significantly reduced the activity levels of all the enzymes. In addition, 3-BrPA treatment led to a reduction in the levels of phosphofruto-1-kinase (PFK) protein, suggesting that the mode of action of 3-BrPA involves the downregulation of some glycolytic enzymes. Measurement of ATP levels in promastigotes of L. amazonensis showed a significant reduction in ATP generation. The O₂ consumption was also significantly inhibited in promastigotes, confirming the energy depletion effect of 3-BrPA. When 3-BrPA was added to the cells at the beginning of growth cycle, it significantly inhibited L. amazonensis proliferation in a dose-dependent manner. Furthermore, the ability to infect macrophages was reduced by approximately 50% when promastigotes were treated with 3-BrPA. Taken together, these studies corroborate with previous reports which suggest 3-BrPA as a potential drug against pathogenic microorganisms that are reliant on glucose catabolism for ATP supply.

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.

Structure determination and analyses of the GAPDH from the parasite Schistosoma mansoni, the first one from a platyhelminth

The enzyme Glyceraldehyde-3-Phosphate Dehydrogenase from Schistosoma mansoni (SmGAPDH) is characterized as a therapeutical target for schistosomiasis. In this context, we report here the experimental structure, structural analyses and comparisons of SmGAPDH, the first one from a Platyhelminth. The enzyme was expressed, purified and assayed for crystallization, what allowed the obtainment of crystals of sufficient quality to collect X-ray diffraction data up to 2.51 Å resolution. SmGAPDH is the only GAPDH to present the sequence NNR (its residues 114-116) which leads to (especially R116) a hydrogen bond network that possibly reflects on the flexibility of residues to interact with the adenine part of NAD⁺, speculated to be important for differential drug design.

Unraveling the conformational dynamics of glycerol 3-phosphate dehydrogenase, a nicotinamide adenine dinucleotide-dependent enzyme of Leishmania mexicana

Allosteric changes modulate the enzymatic activity, leading to activation or inhibition of the molecular target. Understanding the induced fit accommodation mechanism of a ligand in its lowest-free energy state and the subsequent conformational changes induced in the protein are important questions for drug design. In the present study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were applied to analyze the glycerol-3-phosphate dehydrogenase of Leishmania mexicana (LmGPDH) conformational changes induced by its cofactor and substrate binding. GPDH is a nicotinamide adenine dinucleotide (NAD)-dependent enzyme, which has been reported as an interesting target for drug discovery and development against leishmaniasis. Despite its relevance for glycolysis and pentose phosphate pathways, the structural flexibility and conformational motions of LmGPDH in complex with NADH and dihydroxyacetone phosphate (DHAP) remain unexplored. Here, we analyzed the conformational dynamics of the enzyme-NADH complex (cofactor), and the enzyme-NADH-DHAP complex (adduct), mapped the hydrogen-bond interactions for the complexes and pointed some structural determinants of the enzyme that emerge from these contacts to NADH and DHAP. Finally, we proposed a consistent mechanism for the conformational changes on the first step of the reversible redox conversion of dihydroxyacetone phosphate to glycerol 3-phosphate, indicating key residues and interactions that could be further explored in drug discovery.

The Antimicrobials Anacardic Acid and Curcumin Are Not-Competitive Inhibitors of Gram-Positive Bacterial Pathogenic Glyceraldehyde-3-Phosphate Dehydrogenase by a Mechanism Unrelated to Human C5a Anaphylatoxin Binding

The ubiquitous and highly abundant glycolytic enzyme D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is pivotal for the energy and carbon metabolism of most organisms, including human pathogenic bacteria. For bacteria that depend mostly on glycolysis for survival, GAPDH is an attractive target for inhibitor discovery. The availability of high-resolution structures of GAPDH from various pathogenic bacteria is central to the discovery of new antibacterial compounds. We have determined the X-ray crystal structures of two new GAPDH enzymes from Gram-positive bacterial pathogens, Streptococcus pyogenes and Clostridium perfringens. These two structures, and the recent structure of Atopobium vaginae GAPDH, reveal details in the active site that can be exploited for the design of novel inhibitors based on naturally occurring molecules. Two such molecules, anacardic acid and curcumin, have been found to counter bacterial infection in clinical settings, although the cellular targets responsible for their antimicrobial properties remain unknown. We show that both anacardic acid and curcumin inhibit GAPDH from two bacterial pathogens through uncompetitive and non-competitive mechanisms, suggesting GAPDH as a relevant pharmaceutical target for antibacterial development. Inhibition of GAPDH by anacardic acid and curcumin seems to be unrelated to the immune evasion function of pathogenic bacterial GAPDH, since neither natural compound interfere with binding to the human C5a anaphylatoxin.

Identification of common immunodominant antigens of Eimeria tenella, Eimeria acervulina and Eimeria maxima by immunoproteomic analysis

Clinical chicken coccidiosis is mostly caused by simultaneous infection of several Eimeria species, and host immunity against Eimeria is species-specific. It is urgent to identify common immunodominant antigen of Eimeria for developing multivalent anticoccidial vaccines. In this study, sporozoite proteins of Eimeria tenella, Eimeria acervulina and Eimeria maxima were analyzed by two-dimensional electrophoresis (2DE). Western bot analysis was performed on the yielded 2DE gel using antisera of E. tenella E. acervulina and E. maxima respectively. Next, the detected immunodominant spots were identified by comparing the data from MALDI-TOF-MS/MS with available databases. Finally, Eimeria common antigens were identified by comparing amino acid sequence between the three Eimeria species. The results showed that analysis by 2DE of sporozoite proteins detected 629, 626 and 632 protein spots from E. tenella, E. acervulina and E. maxima respectively. Western bot analysis revealed 50 (E. tenella), 64 (E. acervulina) and 57 (E. maxima) immunodominant spots from the sporozoite 2DE gels of the three Eimeria species. The immunodominant spots were identified as 33, 27 and 25 immunodominant antigens of E. tenella, E. acervulina and E. maxima respectively. Fifty-four immunodominant proteins were identified as 18 ortholog proteins among the three Eimeria species. Finally, 5 of the 18 ortholog proteins were identified as common immunodominant antigens including elongation factor 2 (EF-2), 14-3-3 protein, ubiquitin-conjugating enzyme domain-containing protein (UCE) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH). In conclusion, our results not only provide Eimeria sporozoite immunodominant antigen map and additional immunodominant antigens, but also common immunodominant antigens for developing multivalent anticoccidial vaccines.

Crystal Structure of Glyceraldehyde-3-Phosphate Dehydrogenase from the Gram-Positive Bacterial Pathogen A. vaginae, an Immunoevasive Factor that Interacts with the Human C5a Anaphylatoxin

The Gram-positive anaerobic human pathogenic bacterium Atopobium vaginae causes most diagnosed cases of bacterial vaginosis as well as opportunistic infections in immunocompromised patients. In addition to its well-established role in carbohydrate metabolism, D-glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from Streptococcus pyogenes and S. pneumoniae have been reported to act as extracellular virulence factors during streptococcal infections. Here, we report the crystal structure of GAPDH from A. vaginae (AvGAPDH) at 2.19 Å resolution. The refined model has a crystallographic R_free of 22.6%. AvGAPDH is a homotetramer wherein each subunit is bound to a nicotinamide adenine dinucleotide (NAD⁺) molecule. The AvGAPDH enzyme fulfills essential glycolytic as well as moonlight (non-glycolytic) functions, both of which might be targets of chemotherapeutic intervention. We report that AvGAPDH interacts in vitro with the human C5a anaphylatoxin and inhibits C5a-specific granulocyte chemotaxis, thereby suggesting the participation of AvGAPDH in complement-targeted immunoevasion in a context of infection. The availability of high-quality structures of AvGAPDH and other homologous virulence factors from Gram-positive pathogens is critical for drug discovery programs. In this study, sequence and structural differences between AvGAPDH and related bacterial and eukaryotic GAPDH enzymes are reported in an effort to understand how to subvert the immunoevasive properties of GAPDH and evaluate the potential of AvGAPDH as a druggable target.

D-Glyceraldehyde-3-Phosphate Dehydrogenase Structure and Function

Aside from its well-established role in glycolysis, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) has been shown to possess many key functions in cells. These functions are regulated by protein oligomerization , biomolecular interactions, post-translational modifications , and variations in subcellular localization . Several GAPDH functions and regulatory mechanisms overlap with one another and converge around its role in intermediary metabolism. Several structural determinants of the protein dictate its function and regulation. GAPDH is ubiquitously expressed and is found in all domains of life. GAPDH has been implicated in many diseases, including those of pathogenic, cardiovascular, degenerative, diabetic, and tumorigenic origins. Understanding the mechanisms by which GAPDH can switch between its functions and how these functions are regulated can provide insights into ways the protein can be modulated for therapeutic outcomes.

Crystal Structures of Group B Streptococcus Glyceraldehyde-3-Phosphate Dehydrogenase: Apo-Form, Binary and Ternary Complexes

Glyceraldehyde 3-phosphate dehydrogenase or GAPDH is an evolutionarily conserved glycolytic enzyme. It catalyzes the two step oxidative phosphorylation of D-glyceraldehyde 3-phosphate into 1,3-bisphosphoglycerate using inorganic phosphate and NAD+ as cofactor. GAPDH of Group B Streptococcus is a major virulence factor and a potential vaccine candidate. Moreover, since GAPDH activity is essential for bacterial growth it may serve as a possible drug target. Crystal structures of Group B Streptococcus GAPDH in the apo-form, two different binary complexes and the ternary complex are described here. The two binary complexes contained NAD+ bound to 2 (mixed-holo) or 4 (holo) subunits of the tetrameric protein. The structure of the mixed-holo complex reveals the effects of NAD+ binding on the conformation of the protein. In the ternary complex, the phosphate group of the substrate was bound to the new Pi site in all four subunits. Comparison with the structure of human GAPDH showed several differences near the adenosyl binding pocket in Group B Streptococcus GAPDH. The structures also reveal at least three surface-exposed areas that differ in amino acid sequence compared to the corresponding areas of human GAPDH.

Structural analyses to identify selective inhibitors of glyceraldehyde 3-phosphate dehydrogenase-S, a sperm-specific glycolytic enzyme

STUDY HYPOTHESIS

Detailed structural comparisons of sperm-specific glyceraldehyde 3-phosphate dehydrogenase, spermatogenic (GAPDHS) and the somatic glyceraldehyde 3-phosphate dehydrogenase (GAPDH) isozyme should facilitate the identification of selective GAPDHS inhibitors for contraceptive development.

STUDY FINDING

This study identified a small-molecule GAPDHS inhibitor with micromolar potency and >10-fold selectivity that exerts the expected inhibitory effects on sperm glycolysis and motility.

WHAT IS KNOWN ALREADY

Glycolytic ATP production is required for sperm motility and male fertility in many mammalian species. Selective inhibition of GAPDHS, one of the glycolytic isozymes with restricted expression during spermatogenesis, is a potential strategy for the development of a non-hormonal contraceptive that directly blocks sperm function.

STUDY DESIGN, SAMPLES/MATERIALS, METHODS

Homology modeling and x-ray crystallography were used to identify structural features that are conserved in GAPDHS orthologs in mouse and human sperm, but distinct from the GAPDH orthologs present in somatic tissues. We identified three binding pockets surrounding the substrate and cofactor in these isozymes and conducted a virtual screen to identify small-molecule compounds predicted to bind more tightly to GAPDHS than to GAPDH. Following the production of recombinant human and mouse GAPDHS, candidate compounds were tested in dose-response enzyme assays to identify inhibitors that blocked the activity of GAPDHS more effectively than GAPDH. The effects of a selective inhibitor on the motility of mouse and human sperm were monitored by computer-assisted sperm analysis, and sperm lactate production was measured to assess inhibition of glycolysis in the target cell.

MAIN RESULTS AND THE ROLE OF CHANCE

Our studies produced the first apoenzyme crystal structures for human and mouse GAPDHS and a 1.73 Å crystal structure for NAD(+)-bound human GAPDHS, facilitating the identification of unique structural features of this sperm isozyme. In dose-response assays T0501_7749 inhibited human GAPDHS with an IC50 of 1.2 μM compared with an IC50 of 38.5 μM for the somatic isozyme. This compound caused significant reductions in mouse sperm lactate production (P= 0.017 for 100 μM T0501_7749 versus control) and in the percentage of motile mouse and human sperm (P values from <0.05 to <0.0001, depending on incubation conditions).

LIMITATIONS, REASONS FOR CAUTION

The chemical properties of T0501_7749, including limited solubility and nonspecific protein binding, are not optimal for drug development.

WIDER IMPLICATIONS OF THE FINDINGS

This study provides proof-of-principle evidence that GAPDHS can be selectively inhibited, causing significant reductions in sperm glycolysis and motility. These results highlight the utility of structure-based drug design and support further exploration of GAPDHS, and perhaps other sperm-specific isozymes in the glycolytic pathway, as contraceptive targets.

LARGE SCALE DATA

None. Coordinates and data files for three GAPDHS crystal structures were deposited in the RCSB Protein Data Bank (http://www.rcsb.org).

STUDY FUNDING AND COMPETING INTERESTS

This work was supported by grants from the National Institutes of Health (NIH), USA, including U01 HD060481 and cooperative agreement U54 HD35041 as part of the Specialized Cooperative Centers Program in Reproduction and Infertility Research from the Eunice Kennedy Shriver National Institute of Child Health and Human Development, and TW/HD00627 from the NIH Fogarty International Center. Additional support was provided by subproject CIG-05-109 from CICCR, a program of CONRAD, Eastern Virginia Medical School, USA. There are no conflicts of interest.

Computer-Aided Drug Discovery Approaches against the Tropical Infectious Diseases Malaria, Tuberculosis, Trypanosomiasis, and Leishmaniasis

Despite the tremendous improvement in overall global health heralded by the adoption of the Millennium Declaration in the year 2000, tropical infections remain a major health problem in the developing world. Recent estimates indicate that the major tropical infectious diseases, namely, malaria, tuberculosis, trypanosomiasis, and leishmaniasis, account for more than 2.2 million deaths and a loss of approximately 85 million disability-adjusted life years annually. The crucial role of chemotherapy in curtailing the deleterious health and economic impacts of these infections has invigorated the search for new drugs against tropical infectious diseases. The research efforts have involved increased application of computational technologies in mainstream drug discovery programs at the hit identification, hit-to-lead, and lead optimization stages. This review highlights various computer-aided drug discovery approaches that have been utilized in efforts to identify novel antimalarial, antitubercular, antitrypanosomal, and antileishmanial agents. The focus is largely on developments over the past 5 years (2010-2014).

Three-dimensional structures in the design of therapeutics targeting parasitic protozoa: reflections on the past, present and future

Parasitic protozoa cause a range of diseases which threaten billions of human beings. They are responsible for tremendous mortality and morbidity in the least-developed areas of the world. Presented here is an overview of the evolution over the last three to four decades of structure-guided design of inhibitors, leads and drug candidates aiming at targets from parasitic protozoa. Target selection is a crucial and multi-faceted aspect of structure-guided drug design. The major impact of advances in molecular biology, genome sequencing and high-throughput screening is touched upon. The most advanced crystallographic techniques, including XFEL, have already been applied to structure determinations of drug targets from parasitic protozoa. Even cryo-electron microscopy is contributing to our understanding of the mode of binding of inhibitors to parasite ribosomes. A number of projects have been selected to illustrate how structural information has assisted in arriving at promising compounds that are currently being evaluated by pharmacological, pharmacodynamic and safety tests to assess their suitability as pharmaceutical agents. Structure-guided approaches are also applied to incorporate properties into compounds such that they are less likely to become the victim of resistance mechanisms. A great increase in the number of novel antiparasitic compounds will be needed in the future. These should then be combined into various multi-compound therapeutics to circumvent the diverse resistance mechanisms that render single-compound, or even multi-compound, drugs ineffective. The future should also see (i) an increase in the number of projects with a tight integration of structural biology, medicinal chemistry, parasitology and pharmaceutical sciences; (ii) the education of more `medicinal structural biologists' who are familiar with the properties that compounds need to have for a high probability of success in the later steps of the drug-development process; and (iii) the expansion of drug-development capabilities in middle- and low-income countries.

Potential therapeutic targets and the role of technology in developing novel antileishmanial drugs

Leishmaniasis is the most prevalent pathogenic disease in many countries around the world, but there are few drugs available to treat it. Most antileishmanial drugs available are highly toxic, have resistance issues or require hospitalization for their use; therefore, they are not suitable for use in most of the affected countries. Over the past decade, the completion of the genomes of many human pathogens, including that of Leishmania spp., has opened new doors for target identification and validation. Here, we focus on the potential drug targets that can be used for the treatment of leishmaniasis and bring to light how recent technological advances, such as structure-based drug design, structural genomics, and molecular dynamics (MD), can be used to our advantage to develop potent and affordable antileishmanial drugs.

Proteome scale comparative modeling for conserved drug and vaccine targets identification in Corynebacterium pseudotuberculosis

Corynebacterium pseudotuberculosis (Cp) is a pathogenic bacterium that causes caseous lymphadenitis (CLA), ulcerative lymphangitis, mastitis, and edematous to a broad spectrum of hosts, including ruminants, thereby threatening economic and dairy industries worldwide. Currently there is no effective drug or vaccine available against Cp. To identify new targets, we adopted a novel integrative strategy, which began with the prediction of the modelome (tridimensional protein structures for the proteome of an organism, generated through comparative modeling) for 15 previously sequenced C. pseudotuberculosis strains. This pan-modelomics approach identified a set of 331 conserved proteins having 95-100% intra-species sequence similarity. Next, we combined subtractive proteomics and modelomics to reveal a set of 10 Cp proteins, which may be essential for the bacteria. Of these, 4 proteins (tcsR, mtrA, nrdI, and ispH) were essential and non-host homologs (considering man, horse, cow and sheep as hosts) and satisfied all criteria of being putative targets. Additionally, we subjected these 4 proteins to virtual screening of a drug-like compound library. In all cases, molecules predicted to form favorable interactions and which showed high complementarity to the target were found among the top ranking compounds. The remaining 6 essential proteins (adk, gapA, glyA, fumC, gnd, and aspA) have homologs in the host proteomes. Their active site cavities were compared to the respective cavities in host proteins. We propose that some of these proteins can be selectively targeted using structure-based drug design approaches (SBDD). Our results facilitate the selection of C. pseudotuberculosis putative proteins for developing broad-spectrum novel drugs and vaccines. A few of the targets identified here have been validated in other microorganisms, suggesting that our modelome strategy is effective and can also be applicable to other pathogens.

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.

Selection of molecular targets for drug development against trypanosomatids

Trypanosomatid parasites are a group of flagellated protozoa that includes the genera Leishmania and Trypanosoma, which are the causative agents of diseases (leishmaniases, sleeping sickness and Chagas disease) that cause considerable morbidity and mortality, affecting more than 27 million people worldwide. Today no effective vaccines for the prevention of these diseases exist, whereas current chemotherapy is ineffective, mainly due to toxic side effects of current drugs and to the emergence of drug resistance and lack of cost effectiveness. For these reasons, rational drug design and the search of good candidate drug targets is of prime importance. The search for drug targets requires a multidisciplinary approach. To this end, the completion of the genome project of many trypanosomatid species gives a vast amount of new information that can be exploited for the identification of good drug candidates with a prediction of "druggability" and divergence from mammalian host proteins. In addition, an important aspect in the search for good drug targets is the "target identification" and evaluation in a biological pathway, as well as the essentiality of the gene in the mammalian stage of the parasite, which is provided by basic research and genetic and proteomic approaches. In this chapter we will discuss how these bioinformatic tools and experimental evaluations can be integrated for the selection of candidate drug targets, and give examples of metabolic and signaling pathways in the parasitic protozoa that can be exploited for rational drug design.

Pharmacophore binding motifs for nicotinamide adenine dinucleotide analogues across multiple protein families: a detailed contact-based analysis of the interaction between proteins and NAD(P) cofactors

We have analyzed the protein-binding pharmacophore of NAD and its close analogues in all protein-ligand structures available in the RCSB database as of February 2012; this analysis has then been used to assess the novelty of structures emerging after that date. We show that proteins have evolved diverse pharmacophore motifs for binding the adenine moiety, fewer, but still diverse, motifs for nicotinamide, and a very limited set of motifs for binding the pyrophosphate linker. Our exhaustive analysis includes a pharmacophore contact analysis for over 1900 protein-ligand structures containing NAD analogues; we have benchmarked this set of contacts against nearly 27 000 protein-ligand structures to demonstrate that the diversity of interactions seen with NAD is very similar to that seen for all other ligands. Hence, variation in binding motifs for NAD is not distinct from that observed for other ligands and they show significant variation across protein families.

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.

Role of cytosolic glyceraldehyde-3-phosphate dehydrogenase in visceral organ infection by Leishmania donovani

The initial 7 steps of the glycolytic pathway from glucose to 3-phosphoglycerate are localized in the glycosomes in Leishmania, including step 6, catalyzed by the enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH). In L. donovani and L. mexicana, there exists a second GAPDH enzyme present in the cytosol that is absent in L. braziliensis and that has become a pseudogene in L. major. To investigate the role of the cytosolic GAPDH (cGAPDH), an L. donovani cGAPDH-null mutant was generated, and conversely, the functional L. donovani cGAPDH was introduced into L. major and the resulting engineered parasites were characterized. The L. donovani cGAPDH-null mutant was able to proliferate at the same rate as the wild-type parasite in glucose-deficient medium. However, in the presence of glucose, the L. donovani cGAPDH-null mutant consumed less glucose and proliferated more slowly than the wild-type parasite and displayed reduced infectivity in visceral organs of experimentally infected mice. This demonstrates that cGAPDH is functional in L. donovani and is required for survival in visceral organs. Restoration of cGAPDH activity in L. major, in contrast, had an adverse effect on L. major proliferation in glucose-containing medium, providing a possible explanation of why it has evolved into a pseudogene in L. major. This study indicates that there is a difference in glucose metabolism between L. donovani and L. major, and this may represent an important factor in the ability of L. donovani to cause visceral disease.

Anticancer agents that counteract tumor glycolysis

Can we consider cancer to be a "metabolic disease"? Tumors are the result of a metabolic selection, forming tissues composed of heterogeneous cells that generally express an overactive metabolism as a common feature. In fact, cancer cells have increased needs for both energy and biosynthetic intermediates to support their growth and invasiveness. However, their high proliferation rate often generates regions that are insufficiently oxygenated. Therefore, their carbohydrate metabolism must rely mostly on a glycolytic process that is uncoupled from oxidative phosphorylation. This metabolic switch, also known as the Warburg effect, constitutes a fundamental adaptation of tumor cells to a relatively hostile environment, and supports the evolution of aggressive and metastatic phenotypes. As a result, tumor glycolysis may constitute an attractive target for cancer therapy. This approach has often raised concerns that antiglycolytic agents may cause serious side effects toward normal cells. The key to selective action against cancer cells can be found in their hyperbolic addiction to glycolysis, which may be exploited to generate new anticancer drugs with minimal toxicity. There is growing evidence to support many glycolytic enzymes and transporters as suitable candidate targets for cancer therapy. Herein we review some of the most relevant antiglycolytic agents that have been investigated thus far for the treatment of cancer.

Structure and kinetic characterization of human sperm-specific glyceraldehyde-3-phosphate dehydrogenase, GAPDS

hGAPDS (human sperm-specific glyceraldehyde-3-phosphate dehydrogenase) is a glycolytic enzyme essential for the survival of spermatozoa, and constitutes a potential target for non-hormonal contraception. However, enzyme characterization of GAPDS has been hampered by the difficulty in producing soluble recombinant protein. In the present study, we have overexpressed in Escherichia coli a highly soluble form of hGAPDS truncated at the N-terminus (hGAPDSΔN), and crystallized the homotetrameric enzyme in two ligand complexes. The hGAPDSΔN–NAD+–phosphate structure maps the two anion-recognition sites within the catalytic pocket that correspond to the conserved Ps site and the newly recognized Pi site identified in other organisms. The hGAPDSΔN–NAD+–glycerol structure shows serendipitous binding of glycerol at the Ps and new Pi sites, demonstrating the propensity of these anion-recognition sites to bind non-physiologically relevant ligands. A comparison of kinetic profiles between hGAPDSΔN and its somatic equivalent reveals a 3-fold increase in catalytic efficiency for hGAPDSΔN. This may be attributable to subtle amino acid substitutions peripheral to the active centre that influence the charge properties and protonation states of catalytic residues. Our data therefore elucidate structural and kinetic features of hGAPDS that might provide insightful information towards inhibitor development.

Structure of insoluble rat sperm glyceraldehyde-3-phosphate dehydrogenase (GAPDH) via heterotetramer formation with Escherichia coli GAPDH reveals target for contraceptive design

Sperm glyceraldehyde-3-phosphate dehydrogenase has been shown to be a successful target for a non-hormonal contraceptive approach, but the agents tested to date have had unacceptable side effects. Obtaining the structure of the sperm-specific isoform to allow rational inhibitor design has therefore been a goal for a number of years but has proved intractable because of the insoluble nature of both native and recombinant protein. We have obtained soluble recombinant sperm glyceraldehyde-3-phosphate dehydrogenase as a heterotetramer with the Escherichia coli glyceraldehyde-3-phosphate dehydrogenase in a ratio of 1:3 and have solved the structure of the heterotetramer which we believe represents a novel strategy for structure determination of an insoluble protein. A structure was also obtained where glyceraldehyde 3-phosphate binds in the P(s) pocket in the active site of the sperm enzyme subunit in the presence of NAD. Modeling and comparison of the structures of human somatic and sperm-specific glyceraldehyde-3-phosphate dehydrogenase revealed few differences at the active site and hence rebut the long presumed structural specificity of 3-chlorolactaldehyde for the sperm isoform. The contraceptive activity of alpha-chlorohydrin and its apparent specificity for the sperm isoform in vivo are likely to be due to differences in metabolism to 3-chlorolactaldehyde in spermatozoa and somatic cells. However, further detailed analysis of the sperm glyceraldehyde-3-phosphate dehydrogenase structure revealed sites in the enzyme that do show significant difference compared with published somatic glyceraldehyde-3-phosphate dehydrogenase structures that could be exploited by structure-based drug design to identify leads for novel male contraceptives.

An unexpected phosphate binding site in glyceraldehyde 3-phosphate dehydrogenase: crystal structures of apo, holo and ternary complex of Cryptosporidium parvum enzyme

Background

The structure, function and reaction mechanism of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) have been extensively studied. Based on these studies, three anion binding sites have been identified, one 'Ps' site (for binding the C-3 phosphate of the substrate) and two sites, 'Pi' and 'new Pi', for inorganic phosphate. According to the original flip-flop model, the substrate phosphate group switches from the 'Pi' to the 'Ps' site during the multistep reaction. In light of the discovery of the 'new Pi' site, a modified flip-flop mechanism, in which the C-3 phosphate of the substrate binds to the 'new Pi' site and flips to the 'Ps' site before the hydride transfer, was proposed. An alternative model based on a number of structures of B. stearothermophilus GAPDH ternary complexes (non-covalent and thioacyl intermediate) proposes that in the ternary Michaelis complex the C-3 phosphate binds to the 'Ps' site and flips from the 'Ps' to the 'new Pi' site during or after the redox step.

Results

We determined the crystal structure of Cryptosporidium parvum GAPDH in the apo and holo (enzyme + NAD) state and the structure of the ternary enzyme-cofactor-substrate complex using an active site mutant enzyme. The C. parvum GAPDH complex was prepared by pre-incubating the enzyme with substrate and cofactor, thereby allowing free movement of the protein structure and substrate molecules during their initial encounter. Sulfate and phosphate ions were excluded from purification and crystallization steps. The quality of the electron density map at 2Å resolution allowed unambiguous positioning of the substrate. In three subunits of the homotetramer the C-3 phosphate group of the non-covalently bound substrate is in the 'new Pi' site. A concomitant movement of the phosphate binding loop is observed in these three subunits. In the fourth subunit the C-3 phosphate occupies an unexpected site not seen before and the phosphate binding loop remains in the substrate-free conformation. Orientation of the substrate with respect to the active site histidine and serine (in the mutant enzyme) also varies in different subunits.

Conclusion

The structures of the C. parvum GAPDH ternary complex and other GAPDH complexes demonstrate the plasticity of the substrate binding site. We propose that the active site of GAPDH can accommodate the substrate in multiple conformations at multiple locations during the initial encounter. However, the C-3 phosphate group clearly prefers the 'new Pi' site for initial binding in the active site.

Glycolytic enzyme inhibitors in cancer treatment

BACKGROUND

The radio- and chemotherapeutics currently used for the treatment of cancer are widely known to be characterized by a low therapeutic index. An interesting approach to overcoming some of the limits of these techniques is the exploitation of the so-called Warburg effect, which typically characterizes neoplastic cells. Interestingly, this feature has already been utilized with good results, but only for diagnostic purposes (PET and SPECT). From a pharmacological point of view, drugs able to perturb cancer cell metabolism, specifically at the level of glycolysis, may display interesting therapeutic activities in cancer.

OBJECTIVE

The pharmacological actions of these glycolytic enzyme inhibitors, based primarily on ATP depletion, could include: i) amelioration of drug selectivity by exploiting the particular glycolysis addiction of cancer cell; ii) inhibition of energetic and anabolic processes; iii) reduction of hypoxia-linked cancer-cell resistance; iv) reduction of ATP-dependent multi-drug resistance; and v) cytotoxic synergism with conventional cancer treatments.

CONCLUSION

Several glycolytic inhibitors are currently in preclinical and clinical development. Their clinical value as anticancer agents, above all in terms of therapeutic index, strictly depends on a careful reevaluation of the pathophyiological role of the unique metabolism of cancer cells in general and of Warburg effect in particular.

Exploiting 3D structural templates for detection of metal-binding sites in protein structures

High throughput structural genomics efforts have been making the structures of proteins available even before their function has been fully characterized. Therefore, methods that exploit the structural knowledge to provide evidence about the functions of proteins would be useful. Such methods would be needed to complement the sequence-based function annotation approaches. The current study describes generation of 3D-structural motifs for metal-binding sites from the known metalloproteins. It then scans all the available protein structures in the PDB database for putative metal-binding sites. Our analysis predicted more than 1000 novel metal-binding sites in proteins using three-residue templates, and more than 150 novel metal-binding sites using four-residue templates. Prediction of metal-binding site in a yeast protein YDR533c led to the hypothesis that it might function as metal-dependent amidopeptidase. The structural motifs identified by our method present novel metal-binding sites that reveal newer mechanisms for a few well-known proteins.

The Crystal and Solution Structures of Glyceraldehyde-3-phosphate Dehydrogenase Reveal Different Quaternary Structures

The presence of an isoform of glyceraldehyde-3-phosphate dehydrogenase (kmGAPDH1p) associated with the cell wall of a flocculent strain of Kluyveromyces marxianus was the first report of a non-cytosolic localization of a glycolytic enzyme, but the mechanism by which the protein is transported to the cell surface is not known. To identify structural features that could account for the multiple localizations of the protein, the three-dimensional structure of kmGAPDH1p was determined by x-ray crystallography and small angle x-ray scattering. The x-ray crystallographic structure of kmGAPDH1p revealed a dimer, although all GAPDH homologs studied thus far have a tetrameric structure with 222 symmetry. Interestingly, the structure of kmGAPDH1p in solution revealed a tetramer with a 70 degrees tilt angle between the dimers. Moreover, the separation between the centers of the dimers composing the kmGAPDH1p tetramer diminished from 34 to 30 A upon NAD(+) binding, this latter value being similar to the observed in the crystallographic models of GAPDH homologs. The less compact structure of apo-kmGAPDH1p could already be the first image of the transition intermediate between the tetramer observed in solution and the dimeric form found in the crystal structure, which we postulate to exist in vivo because of the protein's multiple subcellular localizations in this yeast species.

Crystal structure of glyceraldehyde-3-phosphate dehydrogenase from Plasmodium falciparum at 2.25 A resolution reveals intriguing extra electron density in the active site

The crystal structure of D-glyceraldehyde-3-phosphate dehydrogenase (PfGAPDH) from the major malaria parasite Plasmodium falciparum is solved at 2.25 A resolution. The structure of PfGAPDH is of interest due to the dependence of the malaria parasite in infected human erythrocytes on the glycolytic pathway for its energy generation. Recent evidence suggests that PfGAPDH may also be required for other critical activities such as apical complex formation. The cofactor NAD(+) is bound to all four subunits of the tetrameric enzyme displaying excellent electron densities. In addition, in all four subunits a completely unexpected large island of extra electron density in the active site is observed, approaching closely the nicotinamide ribose of the NAD(+). This density is most likely the protease inhibitor AEBSF, found in maps from two different crystals. This putative AEBSF molecule is positioned in a crucial location and hence our structure, with expected and unexpected ligands bound, can be of assistance in lead development and design of novel antimalarials.

Selective inhibition of Trypanosoma cruzi GAPDH by “bi-substrate” analogues

A new series of "bi-substrate" analogues have been synthesized as potential inhibitors of the glyceraldehyde-3-phosphate dehydrogenase and one lead compound has been identified that inhibits the enzyme from Trypanosoma cruzi with good affinity and very high (50-fold) specificity.

Structure–activity relationships of novel inhibitors of glyceraldehyde-3-phosphate dehydrogenase

3D QSAR studies were performed on a library of 120 GAPDH inhibitors, including a series of coumarins, flavonoids, and nucleosides. The VolSurf method was successfully used to calculate surface descriptors for protein-ligand affinity and binding site of the enzyme. PCA/PLS analyses have permitted the evaluation of the structural features crucial for potency, selectivity, and favorable pharmacokinetic properties, and are important for the design of new ligands.

Implications of protein flexibility for drug discovery

Proteins are in constant motion between different conformational states with similar energies. This has often been ignored in drug design. However, protein flexibility is fundamental to understanding the ways in which drugs exert biological effects, their binding-site location, binding orientation, binding kinetics, metabolism and transport. Protein flexibility allows increased affinity to be achieved between a drug and its target. This is crucial, because the lipophilicity and number of polar interactions allowed for an oral drug is limited by absorption, distribution, metabolism and toxicology considerations.

Leishmania mexicana glycerol-3-phosphate dehydrogenase showed conformational changes upon binding a bi-substrate adduct

Certain pathogenic trypanosomatids are highly dependent on glycolysis for ATP production, and hence their glycolytic enzymes, including glycerol-3-phosphate dehydrogenase (GPDH), are considered attractive drug targets. The ternary complex structure of Leishmania mexicana GPDH (LmGPDH) with dihydroxyacetone phosphate (DHAP) and NAD(+) was determined to 1.9A resolution as a further step towards understanding this enzyme's mode of action. When compared with the apo and binary complex structures, the ternary complex structure shows an 11 degrees hinge-bending motion of the C-terminal domain with respect to the N-terminal domain. In addition, residues in the C-terminal domain involved in catalysis or substrates binding show significant movements and a previously invisible five-residue loop region becomes well ordered and participates in NAD(+) binding. Unexpectedly, DHAP and NAD(+) appear to form a covalent bond, producing an adduct in the active site of LmGPDH. Modeling a ternary complex glycerol 3-phosphate (G3P) and NAD(+) with LmGPDH identified ten active site residues that are highly conserved among all GPDHs. Two lysine residues, Lys125 and Lys210, that are presumed to be critical in catalysis, were mutated resulting in greatly reduced catalytic activity. Comparison with other structurally related enzymes found by the program DALI suggested Lys210 as a key catalytic residue, which is located on a structurally conserved alpha-helix. From the results of site-directed mutagenesis, molecular modeling and comparison with related dehydrogenases, a catalytic mechanism of LmGPDH and a possible evolutionary scenario of this group of dehydrogenases are proposed.

Anomalous differences of light elements in determining precise binding modes of ligands to glycerol-3-phosphate dehydrogenase

Pathogenic protozoa such as Trypanosome and Leishmania species cause tremendous suffering worldwide. Because of their dependence on glycolysis for energy, the glycolytic enzymes of these organisms, including glycerol-3-phosphate dehydrogenase (GPDH), are considered attractive drug targets. Using the adenine part of NAD as a lead compound, several 2,6-disubstituted purines were synthesized as inhibitors of Leishmania mexicana GPDH (LmGPDH). The electron densities for the inhibitor 2-bromo-6-chloro-purine bound to LmGPDH using a "conventional" wavelength around 1 A displayed a quasisymmetric shape. The anomalous signals from data collected at 1.77 A clearly indicated the positions of the halogen atoms and revealed the multiple binding modes of this inhibitor. Intriguing differences in the observed binding modes of the inhibitor between very similarly prepared crystals illustrate the possibility of crystal-to-crystal variations in protein-ligand complex structures.

Structure of Trypanosoma cruzi glycosomal glyceraldehyde-3-phosphate dehydrogenase complexed with chalepin, a natural product inhibitor, at 1.95 A resolution

The structure of the glycosomal glyceraldehyde-3-phosphate dehydrogenase (gGAPDH) from Trypanosoma cruzi complexed with chalepin, a natural product from Pilocarpus spicatus, has been determined by X-ray crystallography to 1.95 A resolution. The structure is in the apo form without cofactors in the subunits of the tetrameric gGAPDH in the asymmetric unit. Unequivocal density corresponding to the inhibitor was clearly identified in one monomer. The final refined model of the complex shows extensive conformational changes when compared with the native structure. The mode of binding of chalepin to gGAPDH and its implications for inhibitor design are discussed.