A potential target enzyme for trypanocidal drugs revealed by the crystal structure of NAD-dependent glycerol-3-phosphate dehydrogenase from Leishmania mexicana

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.

Kinetic Model for the Heterogeneous Biocatalytic Reactions Using Tethered Cofactors

Understanding the mechanism of interfacial enzyme kinetics is critical to the development of synthetic biological systems for the production of value-added chemicals. Here, the interfacial kinetics of the catalysis of β-nicotinamide adenine dinucleotide (NAD+)-dependent enzymes acting on NAD+ tethered to the surface of silica nanoparticles (SiNPs) has been investigated using two complementary and supporting kinetic approaches: enzyme excess and reactant (NAD+) excess. Kinetic models developed for these two approaches characterize several critical reaction steps including reversible enzyme adsorption, complexation, decomplexation, and catalysis of the surface-bound enzyme/NAD+ complex. The analysis reveals a concentrating effect resulting in a very high local concentration of enzyme and cofactor on the particle surface, in which the enzyme is saturated by surface-bound NAD, facilitating a rate enhancement of enzyme/NAD+ complexation and catalysis. This resulted in high enzyme efficiency within the tethered NAD+ system compared to that of the free enzyme/NAD+ system, which increases with decreasing enzyme concentration. The role of enzyme adsorption onto solid substrates with a tethered catalyst (such as NAD+) has potential for creating highly efficient flow biocatalytic systems.

Lysine245 plays a crucial role in stability and function of glycerol 3-phosphate dehydrogenase (Gpd1) in Saccharomyces cerevisiae

Glycerol 3-phosphate dehydrogenase (Gpd1 isoform) catalyzes the rate limiting step of glycerol synthesis and is a critical component of the osmo-responsive machinery in yeast. The three-dimensional structure of the enzyme is similar to the enzyme from many other organisms, including humans. A recent study with the human enzyme has proposed K120 (K152 in yeast) to be in the correct orientation for catalysis; K204 (K245 in yeast) is out of plane and is not a participant in the catalytic cycle. The current work was carried out to establish the role of K245 in the catalytic cycle of yeast Gpd1. K245A mutant was found to possess lower catalytic activity. Osmotically stressed cells expressing Gpd1 (K245A) showed no change in intracellular glycerol as compared with wild-type cells which showed ~60% increase. Fluorescence microscopy, native polyacrylamide gel electrophoresis (PAGE) analysis, fluorescence spectroscopy, and Thioflavin T spectrofluorimetry showed a relatively unstable, aggregation- and degradation-prone conformation for the mutant. In silico studies showed an aggregation "hotspot" around K245. This study establishes the requirement of K245 for conformational stability and functional adaptation of Gpd1 in Saccharomyces cerevisiae.

Cloning and RNAi-mediated three lethal genes that can be potentially used for Chilo suppressalis (Lepidoptera: Crambidae) management

RNA interference (RNAi) has gained attention in recent years as a viable pest control strategy. Here, RNAi assays were performed to screen the potential functionality of genes in Chilo suppressalis, a serious pest of rice, and to determine their potential for developing a highly targeted molecular control approach. Potential homologs of NADH dehydrogenase (ND), glycerol 3-phosphate dehydrogenase (GPDH) and male specific lethal 3 (MSL3) were cloned from C. suppressalis, and their spatiotemporal gene expression evaluated. The expression of all three genes was higher in the pupal and adult stages than the larval stages and largely higher in the larval head compared to other tissues. Newly hatched larvae exhibited high mortalities and suppressed growth when fed bacteria producing double-stranded RNAs (dsRNAs) corresponding to the three target genes. This study provides insights into the function of ND, GPDH and MSL3 during C. suppressalis larval development and suggests that all may be candidate gene targets for C. suppressalis pest management.

Systematic in silico Evaluation of Leishmania spp. Proteomes for Drug Discovery

Leishmaniasis is a group of neglected infectious diseases, with approximately 1. 3 million new cases each year, for which the available therapies have serious limitations. Therefore, it is extremely important to apply efficient and low-cost methods capable of selecting the best therapeutic targets to speed up the development of new therapies against those diseases. Thus, we propose the use of integrated computational methods capable of evaluating the druggability of the predicted proteomes of Leishmania braziliensis and Leishmania infantum, species responsible for the different clinical manifestations of leishmaniasis in Brazil. The protein members of those proteomes were assessed based on their structural, chemical, and functional contexts applying methods that integrate data on molecular function, biological processes, subcellular localization, drug binding sites, druggability, and gene expression. These data were compared to those extracted from already known drug targets (BindingDB targets), which made it possible to evaluate Leishmania proteomes for their biological relevance and treatability. Through this methodology, we identified more than 100 proteins of each Leishmania species with druggability characteristics, and potential interaction with available drugs. Among those, 31 and 37 proteins of L. braziliensis and L. infantum, respectively, have never been tested as drug targets, and they have shown evidence of gene expression in the evolutionary stage of pharmacological interest. Also, some of those Leishmania targets showed an alignment similarity of <50% when compared to the human proteome, making these proteins pharmacologically attractive, as they present a reduced risk of side effects. The methodology used in this study also allowed the evaluation of opportunities for the repurposing of compounds as anti-leishmaniasis drugs, inferring potential interaction between Leishmania proteins and ~1,000 compounds, of which only 15 have already been tested as a treatment for leishmaniasis. Besides, a list of potential Leishmania targets to be tested using drugs described at BindingDB, such as the potential interaction of the DEAD box RNA helicase, TRYR, and PEPCK proteins with the Staurosporine compound, was made available to the public.

Distinct sequence and structural feature of trypanosoma malate dehydrogenase

Glycosomal malate dehydrogenase from Trypanosoma cruzi (tcgMDH) catalyzes the oxidation/reduction of malate/oxaloacetate, a crucial step of the glycolytic process occurring in the glycosome of the human parasite. Inhibition of tcgMDH is considered a druggable trait for the development of trypanocidal drugs. Sequence comparison of MDHs from different organisms revealed a distinct insertion of a prolin rich 9-mer (62-KLPPVPRDP-70) in tcgMDH as compared to other eukaryotic MDHs. Crystal structure of tcgMDH is solved here at 2.6 Å resolution with R_work/R_free values of 0.206/0.216. The tcgMDH forms homo-dimer with the solvation free energy (ΔG^o) gain of -9.77 kcal/mol. The dimeric form is also confirmed in solution by biochemical assays, chemical-crosslinking and dynamic light scattering. The inserted 9-mer adopts a structure of a solvent accessible loop in the vicinity of NAD⁺ binding site. The distinct sequence and structural feature of tcgMDH, revealed in the present report, provides an anchor point for the development of inhibitors specific for tcgMDH, possible trypanocidal agents of the future.

Decoding the Molecular Effects of Atovaquone Linked Resistant Mutations on Plasmodium falciparum Cytb-ISP Complex in the Phospholipid Bilayer Membrane

Atovaquone (ATQ) is a drug used to prevent and treat malaria that functions by targeting the Plasmodium falciparum cytochrome b (PfCytb) protein. PfCytb catalyzes the transmembrane electron transfer (ET) pathway which maintains the mitochondrial membrane potential. The ubiquinol substrate binding site of the protein has heme bL, heme bH and iron-sulphur [2FE-2S] cluster cofactors that act as redox centers to aid in ET. Recent studies investigating ATQ resistance mechanisms have shown that point mutations of PfCytb confer resistance. Thus, understanding the resistance mechanisms at the molecular level via computational approaches incorporating phospholipid bilayer would help in the design of new efficacious drugs that are also capable of bypassing parasite resistance. With this knowledge gap, this article seeks to explore the effect of three drug resistant mutations Y268C, Y268N and Y268S on the PfCytb structure and function in the presence and absence of ATQ. To draw reliable conclusions, 350 ns all-atom membrane (POPC:POPE phospholipid bilayer) molecular dynamics (MD) simulations with derived metal parameters for the holo and ATQ-bound -proteins were performed. Thereafter, simulation outputs were analyzed using dynamic residue network (DRN) analysis. Across the triplicate MD runs, hydrophobic interactions, reported to be crucial in protein function were assessed. In both, the presence and absence of ATQ and a loss of key active site residue interactions were observed as a result of mutations. These active site residues included: Met 133, Trp136, Val140, Thr142, Ile258, Val259, Pro260 and Phe264. These changes to residue interactions are likely to destabilize the overall intra-protein residue communication network where the proteins' function could be implicated. Protein dynamics of the ATQ-bound mutant complexes showed that they assumed a different pose to the wild-type, resulting in diminished residue interactions in the mutant proteins. In summary, this study presents insights on the possible effect of the mutations on ATQ drug activity causing resistance and describes accurate MD simulations in the presence of the lipid bilayer prior to conducting inhibitory drug discovery for the PfCytb-iron sulphur protein (Cytb-ISP) complex.

Discovery of Alternative Chemotherapy Options for Leishmaniasis through Computational Studies of Asteraceae

Leishmaniasis is a complex disease caused by over 20 Leishmania species that primarily affects populations with poor socioeconomic conditions. Currently available drugs for treating leishmaniasis include amphotericin B, paromomycin, and pentavalent antimonials, which have been associated with several limitations, such as low efficacy, the development of drug resistance, and high toxicity. Natural products are an interesting source of new drug candidates. The Asteraceae family includes more than 23 000 species worldwide. Secondary metabolites that can be found in species from this family have been widely explored as potential new treatments for leishmaniasis. Recently, computational tools have become more popular in medicinal chemistry to establish experimental designs, identify new drugs, and compare the molecular structures and activities of novel compounds. Herein, we review various studies that have used computational tools to examine various compounds identified in the Asteraceae family in the search for potential drug candidates against Leishmania.

Deciphering the Binding of Salicylic Acid to Arabidopsis thaliana Chloroplastic GAPDH-A1

Salicylic acid (SA) has an essential role in the responses of plants to pathogens. SA initiates defence signalling via binding to proteins. NPR1 is a transcriptional co-activator and a key target of SA binding. Many other proteins have recently been shown to bind SA. Amongst these proteins are important enzymes of primary metabolism. This fact could stand behind SA’s ability to control energy fluxes in stressed plants. Nevertheless, only sparse information exists on the role and mechanisms of such binding. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was previously demonstrated to bind SA both in human and plants. Here, we detail that the A1 isomer of chloroplastic glyceraldehyde 3-phosphate dehydrogenase (GAPA1) from Arabidopsis thaliana binds SA with a K_D of 16.7 nM, as shown in surface plasmon resonance experiments. Besides, we show that SA inhibits its GAPDH activity in vitro. To gain some insight into the underlying molecular interactions and binding mechanism, we combined in silico molecular docking experiments and molecular dynamics simulations on the free protein and protein–ligand complex. The molecular docking analysis yielded to the identification of two putative binding pockets for SA. A simulation in water of the complex between SA and the protein allowed us to determine that only one pocket—a surface cavity around Asn35—would efficiently bind SA in the presence of solvent. In silico mutagenesis and simulations of the ligand/protein complexes pointed to the importance of Asn35 and Arg81 in the binding of SA to GAPA1. The importance of this is further supported through experimental biochemical assays. Indeed, mutating GAPA1 Asn35 into Gly or Arg81 into Leu strongly diminished the ability of the enzyme to bind SA. The very same cavity is responsible for the NADP⁺ binding to GAPA1. More precisely, modelling suggests that SA binds to the very site where the pyrimidine group of the cofactor fits. NADH inhibited in a dose-response manner the binding of SA to GAPA1, validating our data.

The unusual di-domain structure of Dunaliella salina glycerol-3-phosphate dehydrogenase enables direct conversion of dihydroxyacetone phosphate to glycerol

Dunaliella has been extensively studied due to its intriguing adaptation to high salinity. Its di-domain glycerol-3-phosphate dehydrogenase (GPDH) isoform is likely to underlie the rapid production of the osmoprotectant glycerol. Here, we report the structure of the chimeric Dunaliella salina GPDH (DsGPDH) protein featuring a phosphoserine phosphatase-like domain fused to the canonical glycerol-3-phosphate (G3P) dehydrogenase domain. Biochemical assays confirm that DsGPDH can convert dihydroxyacetone phosphate (DHAP) directly to glycerol, whereas a separate phosphatase protein is required for this conversion process in most organisms. The structure of DsGPDH in complex with its substrate DHAP and co-factor nicotinamide adenine dinucleotide (NAD) allows the identification of the residues that form the active sites. Furthermore, the structure reveals an intriguing homotetramer form that likely contributes to the rapid biosynthesis of glycerol.

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 Potential of Secondary Metabolites from Plants as Drugs or Leads against Protozoan Neglected Diseases-Part III: In-Silico Molecular Docking Investigations

Malaria, leishmaniasis, Chagas disease, and human African trypanosomiasis continue to cause considerable suffering and death in developing countries. Current treatment options for these parasitic protozoal diseases generally have severe side effects, may be ineffective or unavailable, and resistance is emerging. There is a constant need to discover new chemotherapeutic agents for these parasitic infections, and natural products continue to serve as a potential source. This review presents molecular docking studies of potential phytochemicals that target key protein targets in Leishmania spp., Trypanosoma spp., and Plasmodium spp.

Regulatory networks, genes and glycerophospholipid biosynthesis pathway in schistosomiasis: a systems biology view for pharmacological intervention

Understanding network topology through embracing the global dynamical regulation of genes in an active state space rather than traditional one-gene-one trait approach facilitates the rational drug development process. Schistosomiasis, a neglected tropical disease, has glycerophospholipids as abundant molecules present on its surface. Lack of effective clinical solutions to treat pathogens encourages us to carry out systems-level studies that could contribute to the development of an effective therapy. Development of a strategy for identifying drug targets by combined genome-scale metabolic network and essentiality analyses through in silico approaches provides tantalizing opportunity to investigate the role of protein/substrate metabolism. A genome-scale metabolic network model reconstruction represents choline-phosphate cytidyltransferase as the rate limiting enzyme and regulates the rate of phosphatidylcholine (PC) biosynthesis. The uptake of choline was regulated by choline concentration, promoting the regulation of phosphocholine synthesis. In Schistosoma, the change in developmental stage could result from the availability of choline, hampering its developmental cycle. There are no structural reports for this protein. In order to inhibit the activity of choline-phosphate cytidyltransferase (CCT), it was modeled by homology modeling using 1COZ as the template from Bacillus subtilis. The transition-state stabilization and catalytic residues were mapped as 'HXGH' and 'RTEGISTT' motif. CCT catalyzes the formation of CDP-choline from phosphocholine in which nucleotidyltransferase adds CTP to phosphocholine. The presence of phosphocholine permits the parasite to survive in an immunologically hostile environment. This feature endeavors development of an inhibitor specific for cytidyltransferase in Schistosoma. Flavonolignans were used to inhibit this activity in which hydnowightin showed the highest affinity as compared to miltefosine.

Antileishmanial phytochemical phenolics: molecular docking to potential protein targets

A molecular docking analysis has been carried out to examine potential Leishmania protein targets of antiprotozoal plant-derived polyphenolic compounds. A total of 352 phenolic phytochemicals, including 10 aurones, six cannabinoids, 34 chalcones, 20 chromenes, 52 coumarins, 92 flavonoids, 41 isoflavonoids, 52 lignans, 25 quinones, eight stilbenoids, nine xanthones, and three miscellaneous phenolic compounds, were used in the virtual screening study using 24 Leishmania enzymes (52 different protein structures from the Protein Data Bank). Noteworthy protein targets were Leishmania dihydroorotate dehydrogenase, N-myristoyl transferase, phosphodiesterase B1, pteridine reductase, methionyl-tRNA synthetase, tyrosyl-tRNA synthetase, uridine diphosphate-glucose pyrophosphorylase, nicotinamidase, and glycerol-3-phosphate dehydrogenase. Based on in-silico analysis of antiparasitic polyphenolics in this study, two aurones, one chalcone, five coumarins, six flavonoids, one isoflavonoid, three lignans, and one stilbenoid, can be considered to be promising drug leads worthy of further investigation.

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.

Specificity in transition state binding: the Pauling model revisited

Linus Pauling proposed that the large rate accelerations for enzymes are caused by the high specificity of the protein catalyst for binding the reaction transition state. The observation that stable analogues of the transition states for enzymatic reactions often act as tight-binding inhibitors provided early support for this simple and elegant proposal. We review experimental results that support the proposal that Pauling's model provides a satisfactory explanation for the rate accelerations for many heterolytic enzymatic reactions through high-energy reaction intermediates, such as proton transfer and decarboxylation. Specificity in transition state binding is obtained when the total intrinsic binding energy of the substrate is significantly larger than the binding energy observed at the Michaelis complex. The results of recent studies that aimed to characterize the specificity in binding of the enolate oxygen at the transition state for the 1,3-isomerization reaction catalyzed by ketosteroid isomerase are reviewed. Interactions between pig heart succinyl-coenzyme A:3-oxoacid coenzyme A transferase (SCOT) and the nonreacting portions of coenzyme A (CoA) are responsible for a rate increase of 3 × 10(12)-fold, which is close to the estimated total 5 × 10(13)-fold enzymatic rate acceleration. Studies that partition the interactions between SCOT and CoA into their contributing parts are reviewed. Interactions of the protein with the substrate phosphodianion group provide an ~12 kcal/mol stabilization of the transition state for the reactions catalyzed by triosephosphate isomerase, orotidine 5'-monophosphate decarboxylase, and α-glycerol phosphate dehydrogenase. The interactions of these enzymes with the substrate piece phosphite dianion provide a 6-8 kcal/mol stabilization of the transition state for reaction of the appropriate truncated substrate. Enzyme activation by phosphite dianion reflects the higher dianion affinity for binding to the enzyme-transition state complex compared with that of the free enzyme. Evidence is presented that supports a model in which the binding energy of the phosphite dianion piece, or the phosphodianion group of the whole substrate, is utilized to drive an enzyme conformational change from an inactive open form E(O) to an active closed form E(C), by closure of a phosphodianion gripper loop. Members of the enolase and haloalkanoic acid dehalogenase superfamilies use variable capping domains to interact with nonreacting portions of the substrate and sequester the substrate from interaction with bulk solvent. Interactions of this capping domain with the phenyl group of mandelate have been shown to activate mandelate racemase for catalysis of deprotonation of α-carbonyl carbon. We propose that an important function of these capping domains is to utilize the binding interactions with nonreacting portions of the substrate to activate the enzyme for catalysis.

Biochemical and Molecular-Genetic Characterization of SFD1's Involvement in Lipid Metabolism and Defense Signaling

The Arabidopsis thaliana SFD1 (suppressor of fatty acid desaturase deficiency1) gene (also known as GLY1) is required for accumulation of 34:6 (i.e., 18:3-16:3) monogalactosyldiacylglycerol (MGDG) and for the activation of systemic acquired resistance (SAR), an inducible defense mechanism that confers resistance against a broad spectrum of pathogens. SFD1, which has been suggested to be involved in lipid-based signaling in SAR, contains a putative chloroplast transit peptide and has glycerol-3-phosphate synthesizing dihydroxyacetone phosphate (DHAP) reductase (also referred as glycerol-3-phosphate dehydrogenase) activity. The goals of this study were to determine if the DHAP reductase activity and chloroplast localization are required for SFD1's involvement in galactolipid metabolism and SAR signaling. The crystal structure of a Leishmania mexicana glycerol-3-phosphate dehydrogenase was used to model SFD1 structure and identify Lys194, Lys279, and Asp332 as potential catalytic site residues in SFD1. Mutational analysis of SFD1 confirmed that Lys194, Lys279, and Asp332 are critical for SFD1's DHAP reductase activity, and its involvement in SAR. SFD1 proteins with these residues individually substituted by Ala lacked DHAP reductase activity and were unable to complement the SAR defect of the sfd1 mutant. The SFD1-Ala279 protein was also unable to restore 34:6-MGDG content when expressed in the sfd1 mutant. In vivo imaging of a green fluorescent protein-tagged SFD1 protein demonstrated that SFD1 is targeted to the chloroplast. The N-terminal 43 amino acids, which are required for proper targeting of SFD1 to the chloroplast, are also required for SFD1's function in lipid metabolism and SAR. Taken together, these results demonstrate that SFD1's DHAP reductase activity is required in the chloroplast for lipid metabolism and defense signaling.

Molecular cloning, expression, and immunolocalization of the NAD(+)-dependent glycerol 3-phosphate dehydrogenase (GPD) from Clonorchis sinensis

Glycerol 3-phosphate dehydrogenase (GPD) plays an important role in the energy metabolism and nutrition metabolism. In order to know about the biological functions of GPD of Clonorchis sinensis (C. sinensis), we identified a complete gene coding GPD from C. sinensis metacercaria cDNA library. This novel cDNA sequence contains 1,056 bp with a putative open reading frame of 351 amino acids and shares 74% identity with GPD from Schistosoma mansoni. Recombinant CsGPD was expressed and purified from Escherichia coli BL21 (DE3). Western blot analysis displayed that recombinant CsGPD can be recognized by anti-CsGPD serum and C. sinensis-infected serum. RT-PCR and immunolocalization analysis confirmed that GPD expressed both at the stage of adult worm and metacercaria of C. sinensis and immunolocated at the tegument of adult worm, tegument and tegumentary cells of metacercaria. Our current study has paved the way for the further researches about the biological functions involved in the growth of C. sinensis.

Drug targets in Leishmania

Leishmaniasis is a major public health problem and till date there are no effective vaccines available. The control strategy relies solely on chemotherapy of the infected people. However, the present repertoire of drugs is limited and increasing resistance towards them has posed a major concern. The first step in drug discovery is to identify a suitable drug target. The genome sequences of Leishmania major and Leishmania infantum has revealed immense amount of information and has given the opportunity to identify novel drug targets that are unique to these parasites. Utilization of this information in order to come up with a candidate drug molecule requires combining all the technology and using a multi-disciplinary approach, right from characterizing the target protein to high throughput screening of compounds. Leishmania belonging to the order kinetoplastidae emerges from the ancient eukaryotic lineages. They are quite diverse from their mammalian hosts and there are several cellular processes that we are getting to know of, which exist distinctly in these parasites. In this review, we discuss some of the metabolic pathways that are essential and could be used as potential drug targets in Leishmania.

An examination of the relationship between active site loop size and thermodynamic activation parameters for orotidine 5'-monophosphate decarboxylase from mesophilic and thermophilic organisms

Closure of the active site phosphate gripper loop of orotidine 5'-monophosphate decarboxylase from Saccharomyces cerevisiae (ScOMPDC) over the bound substrate orotidine 5'-monophosphate (OMP) activates the bound substrate for decarboxylation by at least 10(4)-fold [Amyes, T. L., Richard, J. P., and Tait, J. J. (2005) J. Am. Chem. Soc. 127, 15708-15709]. The 19-residue phosphate gripper loop of the mesophilic ScOMPDC is much larger than the nine-residue loop at the ortholog from the thermophile Methanothermobacter thermautotrophicus (MtOMPDC). This difference in loop size results in a small decrease in the total intrinsic phosphate binding energy of the phosphodianion group of OMP from 11.9 to 11.6 kcal/mol, along with a modest decrease in the extent of activation by phosphite dianion of decarboxylation of the truncated substrate 1-(beta-D-erythrofuranosyl)orotic acid. The activation parameters DeltaH(double dagger) and DeltaS(double dagger) for k(cat) for decarboxylation of OMP are 3.6 kcal/mol and 10 cal K(-1) mol(-1) more positive, respectively, for MtOMPDC than for ScOMPDC. We suggest that these differences are related to the difference in the size of the active site loops at the mesophilic ScOMPDC and the thermophilic MtOMPDC. The greater enthalpic transition state stabilization available from the more extensive loop-substrate interactions for the ScOMPDC-catalyzed reaction is largely balanced by a larger entropic requirement for immobilization of the larger loop at this enzyme.

A substrate in pieces: allosteric activation of glycerol 3-phosphate dehydrogenase (NAD+) by phosphite dianion

The ratio of the second-order rate constants for reduction of dihydroxyacetone phosphate (DHAP) and of the neutral truncated substrate glycolaldehyde (GLY) by glycerol 3-phosphate dehydrogenase (NAD (+), GPDH) saturated with NADH is (1.0 x 10 (6) M (-1) s (-1))/(8.7 x 10 (-3) M (-1) s (-1)) = 1.1 x 10 (8), which was used to calculate an intrinsic phosphate binding energy of at least 11.0 kcal/mol. Phosphite dianion binds very weakly to GPDH ( K d > 0.1 M), but the bound dianion strongly activates GLY toward enzyme-catalyzed reduction by NADH. Thus, the large intrinsic phosphite binding energy is expressed only at the transition state for the GPDH-catalyzed reaction. The ratio of rate constants for the phosphite-activated and the unactivated GPDH-catalyzed reduction of GLY by NADH is (4300 M (-2) s (-1))/(8.7 x 10 (-3) M (-1) s (-1)) = 5 x 10 (5) M (-1), which was used to calculate an intrinsic phosphite binding energy of -7.7 kcal/mol for the association of phosphite dianion with the transition state complex for the GPDH-catalyzed reduction of GLY. Phosphite dianion has now been shown to activate bound substrates for enzyme-catalyzed proton transfer, decarboxylation, hydride transfer, and phosphoryl transfer reactions. Structural data provide strong evidence that enzymic activation by the binding of phosphite dianion occurs at a modular active site featuring (1) a binding pocket complementary to the reactive substrate fragment which contains all the active site residues needed to catalyze the reaction of the substrate piece or of the whole substrate and (2) a phosphate/phosphite dianion binding pocket that is completed by the movement of flexible protein loop(s) to surround the nonreacting oxydianion. We propose that loop motion and associated protein conformational changes that accompany the binding of phosphite dianion and/or phosphodianion substrates lead to encapsulation of the substrate and/or its pieces in the protein interior, and to placement of the active site residues in positions where they provide optimal stabilization of the transition state for the catalyzed reaction.

New competitive inhibitors of cytosolic (NADH-dependent) rabbit muscle glycerophosphate dehydrogenase

We report the synthesis and biochemical evaluation of new competitive inhibitors of the cytosolic (NADH-dependent) glycerophosphate dehydrogenase. The best tested compound, phosphono-propionohydroxamic acid, with a Ki of 6 microM, might be of interest as an anti-obesity drug.

Evolution of protein fold in the presence of functional constraints

The functional requirement to form and maintain the active site structure probably exerts a strong selective pressure on a protein to adopt just one stable and evolutionarily conserved fold. Nonetheless, new evidence suggests the likelihood of protein fold being neither physically nor biologically invariant. Alternative folds discovered in several proteins are composed of constant and variable parts. The latter display context-dependent conformations and a tendency to form new oligomeric interfaces. In turn, oligomerisation mediates fold evolution without loss of protein function. Gene duplication breaks down homo-oligomeric symmetry and relieves the pressure to maintain the local architecture of redundant active sites; this can lead to further structural changes.

Crystal structures of human glycerol 3-phosphate dehydrogenase 1 (GPD1)

Homo sapiens L-alpha-glycerol-3-phosphate dehydrogenase 1 (GPD1) catalyzes the reversible biological conversion of dihydroxyacetone (DHAP) to glycerol-3-phosphate. The GPD1 protein was expressed in Escherichia coli, and purified as a fusion protein with glutathione S-transferase. Here we report the apoenzyme structure of GPD1 determined by multiwavelength anomalous diffraction phasing, and other complex structures with small molecules (NAD+ and DHAP) by the molecular replacement method. This enzyme structure is organized into two distinct domains, the N-terminal eight-stranded beta-sheet sandwich domain and the C-terminal helical substrate-binding domain. An electrophilic catalytic mechanism by the epsilon-NH3+ group of Lys204 is proposed on the basis of the structural analyses. In addition, the inhibitory effects of zinc and sulfate on GPDHs are assayed and discussed.

Structural and functional analysis of the gpsA gene product of Archaeoglobus fulgidus: a glycerol-3-phosphate dehydrogenase with an unusual NADP+ preference

NAD(+)-dependent glycerol-3-phosphate dehydrogenase (G3PDH) is generally absent in archaea, because archaea, unlike eukaryotes and eubacteria, utilize glycerol-1-phosphate instead of glycerol-3-phosphate for the biosynthesis of membrane lipids. Surprisingly, the genome of the hyperthermophilic archaeon Archaeoglobus fulgidus comprises a G3PDH ortholog, gpsA, most likely due to horizontal gene transfer from a eubacterial organism. Biochemical characterization proved G3PDH-like activity of the recombinant gpsA gene product. However, unlike other G3PDHs, the up to 85 degrees C thermostable A. fulgidus G3PDH exerted a 15-fold preference for NADPH over NADH. The A. fulgidus G3PDH bears the hallmarks of adaptation to halotolerance and thermophilicity, because its 1.7-A crystal structure showed a high surface density for negative charges and 10 additional intramolecular salt bridges compared to a mesophilic G3PDH structure. Whereas all amino acid residues required for dihydroxyacetone phosphate binding and reductive catalysis are highly conserved, the binding site for the adenine moiety of the NAD(P) cosubstrate shows a structural variation that reflects the observed NADPH preference, for example, by a putative salt bridge between R49 and the 2'-phosphate.

4SCOPmap: automated assignment of protein structures to evolutionary superfamilies

Background

Inference of remote homology between proteins is very challenging and remains a prerogative of an expert. Thus a significant drawback to the use of evolutionary-based protein structure classifications is the difficulty in assigning new proteins to unique positions in the classification scheme with automatic methods. To address this issue, we have developed an algorithm to map protein domains to an existing structural classification scheme and have applied it to the SCOP database.

Results

The general strategy employed by this algorithm is to combine the results of several existing sequence and structure comparison tools applied to a query protein of known structure in order to find the homologs already classified in SCOP database and thus determine classification assignments. The algorithm is able to map domains within newly solved structures to the appropriate SCOP superfamily level with ~95% accuracy. Examples of correctly mapped remote homologs are discussed. The algorithm is also capable of identifying potential evolutionary relationships not specified in the SCOP database, thus helping to make it better. The strategy of the mapping algorithm is not limited to SCOP and can be applied to any other evolutionary-based classification scheme as well. SCOPmap is available for download.

Conclusion

The SCOPmap program is useful for assigning domains in newly solved structures to appropriate superfamilies and for identifying evolutionary links between different superfamilies.

Quinoxaline N , N ′-dioxide derivatives and related compounds as growth inhibitors of Trypanosoma cruzi . Structure–activity relationships

Quinoxaline derivatives presented good inhibitor activity of growth of Trypanosoma cruzi in in vitro assays. The 50% inhibitory doses were of the same order of that of Nifurtimox. Derivative 13, a quinoxaline N,N'-dioxide derivative, and the reduced derivatives 19 and 20 were the most cytotoxic compounds against the protozoan. Structural requirements for optimal activity were studied by computational methods. From statistical analysis we could establish a multiple correlation between activity and lipophilic properties and LUMO energy.

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.

Crystal structure of Pseudomonas fluorescens mannitol 2-dehydrogenase: evidence for a very divergent long-chain dehydrogenase family

Mannitol 2-dehydrogenase from Pseudomonas fluorescens (pfMDH) is a secondary alcohol dehydrogenase that catalyzes the reversible NAD(P)-dependent oxidation of D-mannitol to D-fructose, D-arabinitol to D-xylulose, and D-sorbitol to L-sorbose. It is a member of the mostly prokaryotic family of long-chain mannitol dehydrogenases that so far includes 66 members. Unlike other alcohol and polyol dehydrogenases that utilize metal cofactors or a conserved active-site tyrosine for catalysis, an invariant lysine is the general base. The crystal structure of pfMDH in a binary complex with NAD(H) and a ternary complex with NAD(H) and D-mannitol have been determined to 1.7 and 1.8 A resolution respectively. Comparison of secondary structure assignment to sequence alignments suggest the shortest members of this family, mannitol-1-phosphate 5-dehydrogenases, retain core elements but lack secondary structural components found on the surface of pfMDH. The elements predicted to be absent are distributed throughout the primary sequence, implying that a simple truncation or fusion did not occur. The closest structural neighbors are 6-phosphogluconate dehydrogenase, UDP-glucose dehydrogenase, N-(1-D-carboxyethyl)-L-norvaline dehydrogenase, and glycerol-3-phosphate dehydrogenase. Although sequence identity is only a barely recognizable 7-10%, conservation of secondary structural elements as well as homologous residues that are contributed to the active site indicates they may be related by divergent evolution.

Analysis of membrane stereochemistry with homology modeling of sn-glycerol-1-phosphate dehydrogenase

Different enantiomeric isomers, sn-glycerol-1-phosphate and sn-glycerol-3-phosphate, are used as the glycerophosphate backbones of phospholipids in the cellular membranes of Archaea and the remaining two kingdoms, respectively. In Archaea, sn-glycerol-1-phosphate dehydrogenase is involved in the generation of sn-glycerol-1-phosphate, while sn-glycerol-3-phosphate dehydrogenase synthesizes the enantiomer in Eukarya and Bacteria. The coordinates of sn-glycerol-3-phosphate dehydrogenase are available, although neither the tertiary structure nor the reaction mechanism of sn-glycerol-1-phosphate dehydrogenase is known. Database searching revealed that the archaeal enzyme shows sequence similarity to glycerol dehydrogenase, dehydroquinate synthase and alcohol dehydrogenase IV. The glycerol dehydrogenase, with coordinates that are available today, is closely related to the archaeal enzyme. Using the structure of glycerol dehydrogenase as the template, we built a model structure of the Methanothermobacter thermautotrophicus sn-glycerol-1-phosphate dehydrogenase, which could explain the chirality of the product. Based on the model structure, we determined the following: (1) the enzyme requires a Zn(2+) ion for its activity; (2) the enzyme selectively uses the pro-R hydrogen of the NAD(P)H; (3) the putative active site and the reaction mechanism were predicted; and (4) the archaeal enzyme does not share its evolutionary origin with sn-glycerol-3-phosphate dehydrogenase.

Crystal structure of Pseudomonas fluorescens mannitol 2-dehydrogenase binary and ternary complexes. Specificity and catalytic mechanism

Long-chain mannitol dehydrogenases are secondary alcohol dehydrogenases that are of wide interest because of their involvement in metabolism and potential applications in agriculture, medicine, and industry. They differ from other alcohol and polyol dehydrogenases because they do not contain a conserved tyrosine and are not dependent on Zn(2+) or other metal cofactors. The structures of the long-chain mannitol 2-dehydrogenase (54 kDa) from Pseudomonas fluorescens in a binary complex with NAD(+) and ternary complex with NAD(+) and d-mannitol have been determined to resolutions of 1.7 and 1.8 A and R-factors of 0.171 and 0.176, respectively. These results show an N-terminal domain that includes a typical Rossmann fold. The C-terminal domain is primarily alpha-helical and mediates mannitol binding. The electron lone pair of Lys-295 is steered by hydrogen-bonding interactions with the amide oxygen of Asn-300 and the main-chain carbonyl oxygen of Val-229 to act as the general base. Asn-191 and Asn-300 are involved in a web of hydrogen bonding, which precisely orients the mannitol O2 proton for abstraction. These residues also aid in stabilizing a negative charge in the intermediate state and in preventing the formation of nonproductive complexes with the substrate. The catalytic lysine may be returned to its unprotonated state using a rectifying proton tunnel driven by Glu-292 oscillating among different environments. Despite low sequence homology, the closest structural neighbors are glycerol-3-phosphate dehydrogenase, N-(1-d-carboxylethyl)-l-norvaline dehydrogenase, UDP-glucose dehydrogenase, and 6-phosphogluconate dehydrogenase, indicating a possible evolutionary relationship among these enzymes.

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.

Inhibition of plasmodium falciparum triose-phosphate isomerase by chemical modification of an interface cysteine. Electrospray ionization mass spectrometric analysis of differential cysteine reactivities

Plasmodium falciparum triose-phosphate isomerase, a homodimeric enzyme, contains four cysteine residues at positions 13, 126, 196, and 217 per subunit. Among these, Cys-13 is present at the dimer interface and is replaced by methionine in the corresponding human enzyme. We have investigated the effect of sulfhydryl labeling on the parasite enzyme, with a view toward developing selective covalent inhibitors by targeting the interface cysteine residue. Differential labeling of the cysteine residues by iodoacetic acid and iodoacetamide has been followed by electrospray ionization mass spectrometry and positions of the labels determined by analysis of tryptic fragments. The rates of labeling follows the order Cys-196 > Cys-13 Cys-217/Cys-126, which correlates well with surface accessibility calculations based on the enzyme crystal structure. Iodoacetic acid labeling leads to a soluble, largely inactive enzyme, whereas IAM labeling leads to precipitation. Carboxyl methylation of Cys-13 results in formation of monomeric species detectable by gel filtration. Studies with an engineered C13D mutant permitted elucidation of the effects of introducing a negative charge at the interface. The C13D mutant exhibits a reduced stability to denaturants and 7-fold reduction in the enzymatic activity even under the concentrations in which dimeric species are observed.

Homology modeling and molecular dynamics study of NAD-dependent glycerol-3-phosphate dehydrogenase from Trypanosoma brucei rhodesiense, a potential target enzyme for anti-sleeping sickness drug development

Sleeping sickness and Chagas disease are among the most severe diseases in Africa as well as Latin America. These two diseases are caused by Trypanosoma spp. Recently, an enzyme of a glycolytic pathway, NAD-dependent glycerol-3-phosphate dehydrogenase, of Leishmania mexicana was crystallized and its structure determined by x-ray crystallography. This structure has offered an excellent template for modeling of the homologous enzymes from another Trypanosoma species. Here, a homology model of the T. brucei enzyme based on the x-ray structure of LmGPDH has been generated. This model was used as the starting point for molecular dynamics simulation in a water box. The analysis of the molecular dynamics trajectory indicates that the functionally important motifs have both a very stable secondary structure and tertiary arrangement.

Biogenesis and function of peroxisomes and glycosomes

Peroxisomes of higher eukaryotes, glycosomes of kinetoplastids, and glyoxysomes of plants are related microbody organelles that perform differing metabolic functions tailored to their cellular environments. The close evolutionary relationship of these organelles is most clearly evidenced by the conservation of proteins involved in matrix protein import and biogenesis. The glycosome can be viewed as an offshoot of the peroxisomal lineage with additional metabolic functions, specifically glycolysis and purine salvage. Within the parasitic protozoa, only kinetoplastids have been conclusively demonstrated to possess glycosomes or indeed any peroxisome-like organelle. The importance of glycosomal pathways and their compartmentation emphasizes the potential of the glycosome and glycosomal proteins as drug targets.

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

The glycolytic enzymes of trypanosomes are attractive drug targets, since the blood-stream form of Trypanosoma brucei lacks a functional citric acid cycle and is dependent solely on glycolysis for its energy requirements. Glyceraldehyde-3-phosphate dehydrogenases (GAPDH) from the pathogenic trypanosomatids T. brucei, Trypanosoma cruzi and Leishmania mexicana are quite similar to each other, and yet have sufficient structural differences compared to the human enzyme to enable the structure-based design of compounds that selectively inhibit all three trypanosomatid enzymes but not the human homologue. Adenosine analogs with substitutions on N-6 of the adenine ring and on the 2' position of the ribose moiety were designed, synthesized and tested for inhibition. Two crystal structures of L. mexicana glyceraldehyde-3-phosphate dehydrogenase in complex with high-affinity inhibitors that also block parasite growth were solved at a resolution of 2.6 A and 3.0 A. The complexes crystallized in the same crystal form, with one and a half tetramers in the crystallographic asymmetric unit. There is clear electron density for the inhibitor in all six copies of the binding site in each of the two structures. The L. mexicana GAPDH subunit exhibits substantial structural plasticity upon binding the inhibitor. Movements of the protein backbone, in response to inhibitor binding, enlarge a cavity at the binding site to accommodate the inhibitor in a classic example of induced fit. The extensive hydrophobic interactions between the protein and the two substituents on the adenine scaffold of the inhibitor provide a plausible explanation for the high affinity of these inhibitors for trypanosomatid GAPDHs.

Enzymology, structure, and dynamics of acetohydroxy acid isomeroreductase

Acetohydroxy acid isomeroreductase is a key enzyme involved in the biosynthetic pathway of the amino acids isoleucine, valine, and leucine. This enzyme is of great interest in agrochemical research because it is present only in plants and microorganisms, making it a potential target for specific herbicides and fungicides. Moreover, it catalyzes an unusual two-step reaction that is of great fundamental interest. With a view to characterizing both the mechanism of inhibition by potential herbicides and the complex reaction mechanism, various techniques of enzymology, molecular biology, mass spectrometry, X-ray crystallography, and theoretical simulation have been used. The results and conclusions of these studies are described briefly in this paper.

Enzymes of carbohydrate metabolism as potential drug targets

The potential for chemotherapeutic exploitation of carbohydrate metabolism in the Trypanosomatidae is reviewed. This review is based largely on discussions held at a meeting of the COST B9 Action, entitled 'Bioenergetics of Protozoan Parasites'. The major questions posed were: which enzymes are the best to target; what further information is required to allow their use for rational drug development; what compounds would constitute the best inhibitors and which of the enzymes of the pentose-phosphate pathway are present inside the glycosomes, as well? Only partial answers could be obtained in many cases, but the interactive discussion between the multidisciplinary group of participants, comprising chemists, biochemists and molecular biologists, provided thought-provoking ideas and will help direct future research.

Glycolysis as a target for the design of new anti-trypanosome drugs

Glycolysis is perceived as a promising target for new drugs against parasitic trypanosomatid protozoa because this pathway plays an essential role in their ATP supply. Trypanosomatid glycolysis is unique in that it is compartmentalized, and many of its enzymes display unique structural and kinetic features. Structure- and catalytic mechanism-based approaches are applied to design compounds that inhibit the glycolytic enzymes of the parasites without affecting the corresponding proteins of the human host. For some trypanosomatid enzymes, potent and selective inhibitors have already been developed that affect only the growth of cultured trypanosomatids, and not mammalian cells.

Adenosine analogues as inhibitors of Trypanosoma brucei phosphoglycerate kinase: elucidation of a novel binding mode for a 2-amino-N(6)-substituted adenosine

As part of a project aimed at structure-based design of adenosine analogues as drugs against African trypanosomiasis, N(6)-, 2-amino-N(6)-, and N(2)-substituted adenosine analogues were synthesized and tested to establish structure-activity relationships for inhibiting Trypanosoma brucei glycosomal phosphoglycerate kinase (PGK), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), and glycerol-3-phosphate dehydrogenase (GPDH). Evaluation of X-ray structures of parasite PGK, GAPDH, and GPDH complexed with their adenosyl-bearing substrates led us to generate a series of adenosine analogues which would target all three enzymes simultaneously. There was a modest preference by PGK for N(6)-substituted analogues bearing the 2-amino group. The best compound in this series, 2-amino-N(6)- [2''(p-hydroxyphenyl)ethyl]adenosine (46b), displayed a 23-fold improvement over adenosine with an IC(50) of 130 microM. 2-[[2''-(p-Hydroxyphenyl)ethyl]amino]adenosine (46c) was a weak inhibitor of T. brucei PGK with an IC(50) of 500 microM. To explore the potential of an additive effect that having the N(6) and N(2) substitutions in one molecule might provide, the best ligands from the two series were incorporated into N(6),N(2)-disubstituted adenosine analogues to yield N(6)-(2''-phenylethyl)-2-[(2'' -phenylethyl)amino]adenosine (69) as a 30 microM inhibitor of T. brucei PGK which is 100-fold more potent than the adenosine template. In contrast, these series gave no compounds that inhibited parasitic GAPDH or GPDH more than 10-20% when tested at 1.0 mM. A 3.0 A X-ray structure of a T. brucei PGK/46b complex revealed a binding mode in which the nucleoside analogue was flipped and the ribosyl moiety adopted a syn conformation as compared with the previously determined binding mode of ADP. Molecular docking experiments using QXP and SAS program suites reproduced this "flipped and rotated" binding mode.