Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex

Ligand-Based Virtual Screening and Molecular Docking of Benzimidazoles as Potential Inhibitors of Triosephosphate Isomerase Identified New Trypanocidal Agents

Trypanosoma cruzi (T. cruzi) is a parasite that affects humans and other mammals. T. cruzi depends on glycolysis as a source of adenosine triphosphate (ATP) supply, and triosephosphate isomerase (TIM) plays a key role in this metabolic pathway. This enzyme is an attractive target for the design of new trypanocidal drugs. In this study, a ligand-based virtual screening (LBVS) from the ZINC15 database using benzimidazole as a scaffold was accomplished. Later, a molecular docking on the interface of T. cruzi TIM (TcTIM) was performed and the compounds were grouped by interaction profiles. Subsequently, a selection of compounds was made based on cost and availability for in vitro evaluation against blood trypomastigotes. Finally, the compounds were analyzed by molecular dynamics simulation, and physicochemical and pharmacokinetic properties were determined using SwissADME software. A total of 1604 molecules were obtained as potential TcTIM inhibitors. BP2 and BP5 showed trypanocidal activity with half-maximal lytic concentration (LC₅₀) values of 155.86 and 226.30 µM, respectively. Molecular docking and molecular dynamics simulation analyzes showed a favorable docking score of BP5 compound on TcTIM. Additionally, BP5 showed a low docking score (−5.9 Kcal/mol) on human TIM compared to the control ligand (−7.2 Kcal/mol). Both compounds BP2 and BP5 showed good physicochemical and pharmacokinetic properties as new anti-T. cruzi agents.

Highly Similar Sequence and Structure Yet Different Biophysical Behavior: A Computational Study of Two Triosephosphate Isomerases

Homodimeric triosephosphate isomerases (TIMs) from Trypanosoma cruzi (TcTIM) and Trypanosoma brucei (TbTIM) have markedly similar amino-acid sequences and three-dimensional structures. However, several of their biophysical parameters, such as their susceptibility to sulfhydryl agents and their reactivation speed after being denatured, have significant differences. The causes of these differences were explored with microsecond-scale molecular dynamics (MD) simulations of three different TIM proteins: TcTIM, TbTIM, and a chimeric protein, Mut1. We examined their electrostatic interactions and explored the impact of simulation length on them. The same salt bridge between catalytic residues Lys 14 and Glu 98 was observed in all three proteins, but key differences were found in other interactions that the catalytic amino acids form. In particular, a cation-π interaction between catalytic amino acids Lys 14 and His 96 and both a salt bridge and a hydrogen bond between catalytic Glu 168 and residue Arg 100 were only observed in TcTIM. Furthermore, although TcTIM forms less hydrogen bonds than TbTIM and Mut1, its hydrogen bond network spans almost the entire protein, connecting the residues in both monomers. This work provides new insight into the mechanisms that give rise to the different behavior of these proteins. The results also show the importance of long simulations.

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

BACKGROUND

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

OBJECTIVE

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

RESULTS

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

CONCLUSION

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

Structural basis for the modulation of plant cytosolic triosephosphate isomerase activity by mimicry of redox-based modifications

Reactive oxidative species (ROS) and S-glutathionylation modulate the activity of plant cytosolic triosephosphate isomerases (cTPI). Arabidopsis thaliana cTPI (AtcTPI) is subject of redox regulation at two reactive cysteines that function as thiol switches. Here we investigate the role of these residues, AtcTPI-Cys13 and At-Cys218, by substituting them with aspartic acid that mimics the irreversible oxidation of cysteine to sulfinic acid and with amino acids that mimic thiol conjugation. Crystallographic studies show that mimicking AtcTPI-Cys13 oxidation promotes the formation of inactive monomers by reposition residue Phe75 of the neighboring subunit, into a conformation that destabilizes the dimer interface. Mutations in residue AtcTPI-Cys218 to Asp, Lys, or Tyr generate TPI variants with a decreased enzymatic activity by creating structural modifications in two loops (loop 7 and loop 6) whose integrity is necessary to assemble the active site. In contrast with mutations in residue AtcTPI-Cys13, mutations in AtcTPI-Cys218 do not alter the dimeric nature of AtcTPI. Therefore, modifications of residues AtcTPI-Cys13 and AtcTPI-Cys218 modulate AtcTPI activity by inducing the formation of inactive monomers and by altering the active site of the dimeric enzyme, respectively. The identity of residue AtcTPI-Cys218 is conserved in the majority of plant cytosolic TPIs, this conservation and its solvent-exposed localization make it the most probable target for TPI regulation upon oxidative damage by reactive oxygen species. Our data reveal the structural mechanisms by which S-glutathionylation protects AtcTPI from irreversible chemical modifications and re-routes carbon metabolism to the pentose phosphate pathway to decrease oxidative stress.

Enzyme Architecture: Amino Acid Side-Chains That Function To Optimize the Basicity of the Active Site Glutamate of Triosephosphate Isomerase

We report pH rate profiles for k_cat and K_m for the isomerization reaction of glyceraldehyde 3-phosphate catalyzed by wildtype triosephosphate isomerase (TIM) from three organisms and by ten mutants of TIM; and, for K_i for inhibition of this reaction by phosphoglycolate trianion (I^3–). The pH profiles for K_i show that the binding of I^3– to TIM (E) to form EH·I₃^– is accompanied by uptake of a proton by the carboxylate side-chain of E165, whose function is to abstract a proton from substrate. The complexes for several mutants exist mainly as E^–·I₃^– at high pH, in which cases the pH profiles define the pK_a for deprotonation of EH·I₃^–. The linear free energy correlation, with slope of 0.73 (r² = 0.96), between k_cat/K_m for TIM-catalyzed isomerization and the disassociation constant of PGA trianion for TIM shows that EH·I₃^– and the transition state are stabilized by similar interactions with the protein catalyst. Values of pK_a = 10–10.5 were estimated for deprotonation of EH·I₃^– for wildtype TIM. This pK_a decreases to as low as 6.3 for the severely crippled Y208F mutant. There is a correlation between the effect of several mutations on k_cat/K_m and on pK_a for EH·I₃^–. The results support a model where the strong basicity of E165 at the complex to the enediolate reaction intermediate is promoted by side-chains from Y208 and S211, which serve to clamp loop 6 over the substrate; I170, which assists in the creation of a hydrophobic environment for E165; and P166, which functions in driving the carboxylate side-chain of E165 toward enzyme-bound substrate.

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

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

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.

Structural insights from a novel invertebrate triosephosphate isomerase from Litopenaeus vannamei

Triosephosphate isomerase (TIM; EC 5.3.1.1) is a key enzyme involved in glycolysis and gluconeogenesis. Glycolysis is one of the most regulated metabolic pathways, however little is known about the structural mechanisms for its regulation in non-model organisms, like crustaceans. To understand the structure and function of this enzyme in invertebrates, we obtained the crystal structure of triosephosphate isomerase from the marine Pacific whiteleg shrimp (Litopenaeus vannamei, LvTIM) in complex with its inhibitor 2-phosphogyceric acid (2-PG) at 1.7Å resolution. LvTIM assembles as a homodimer with residues 166-176 covering the active site and residue Glu166 interacting with the inhibitor. We found that LvTIM is the least stable TIM characterized to date, with the lowest range of melting temperatures, and with the lowest activation enthalpy associated with the thermal unfolding process reported. In TIMs dimer stabilization is maintained by an interaction of loop 3 by a set of hydrophobic contacts between subunits. Within these contacts, the side chain of a hydrophobic residue of one subunit fits into a cavity created by a set of hydrophobic residues in the neighboring subunit, via a "ball and socket" interaction. LvTIM presents a Cys47 at the "ball" inter-subunit contact indicating that the character of this residue is responsible for the decrease in dimer stability. Mutational studies show that this residue plays a role in dimer stability but is not a solely determinant for dimer formation.

Crystal structures of two monomeric triosephosphate isomerase variants identified via a directed-evolution protocol selecting for L-arabinose isomerase activity

The crystal structures are described of two variants of A-TIM: Ma18 (2.7 Å resolution) and Ma21 (1.55 Å resolution). A-TIM is a monomeric loop-deletion variant of triosephosphate isomerase (TIM) which has lost the TIM catalytic properties. Ma18 and Ma21 were identified after extensive directed-evolution selection experiments using an Escherichia coli L-arabinose isomerase knockout strain expressing a randomly mutated A-TIM gene. These variants facilitate better growth of the Escherichia coli selection strain in medium supplemented with 40 mM L-arabinose. Ma18 and Ma21 differ from A-TIM by four and one point mutations, respectively. Ma18 and Ma21 are more stable proteins than A-TIM, as judged from CD melting experiments. Like A-TIM, both proteins are monomeric in solution. In the Ma18 crystal structure loop 6 is open and in the Ma21 crystal structure loop 6 is closed, being stabilized by a bound glycolate molecule. The crystal structures show only small differences in the active site compared with A-TIM. In the case of Ma21 it is observed that the point mutation (Q65L) contributes to small structural rearrangements near Asn11 of loop 1, which correlate with different ligand-binding properties such as a loss of citrate binding in the active site. The Ma21 structure also shows that its Leu65 side chain is involved in van der Waals interactions with neighbouring hydrophobic side-chain moieties, correlating with its increased stability. The experimental data suggest that the increased stability and solubility properties of Ma21 and Ma18 compared with A-TIM cause better growth of the selection strain when coexpressing Ma21 and Ma18 instead of A-TIM.

The activating oxydianion binding domain for enzyme-catalyzed proton transfer, hydride transfer, and decarboxylation: specificity and enzyme architecture

The kinetic parameters for activation of yeast triosephosphate isomerase (ScTIM), yeast orotidine monophosphate decarboxylase (ScOMPDC), and human liver glycerol 3-phosphate dehydrogenase (hlGPDH) for catalysis of reactions of their respective phosphodianion truncated substrates are reported for the following oxydianions: HPO3(2-), FPO3(2-), S2O3(2-), SO4(2-) and HOPO3(2-). Oxydianions bind weakly to these unliganded enzymes and tightly to the transition state complex (E·S(‡)), with intrinsic oxydianion Gibbs binding free energies that range from -8.4 kcal/mol for activation of hlGPDH-catalyzed reduction of glycolaldehyde by FPO3(2-) to -3.0 kcal/mol for activation of ScOMPDC-catalyzed decarboxylation of 1-β-d-erythrofuranosyl)orotic acid by HOPO3(2-). Small differences in the specificity of the different oxydianion binding domains are observed. We propose that the large -8.4 kcal/mol and small -3.8 kcal/mol intrinsic oxydianion binding energy for activation of hlGPDH by FPO3(2-) and S2O3(2-), respectively, compared with activation of ScTIM and ScOMPDC reflect stabilizing and destabilizing interactions between the oxydianion -F and -S with the cationic side chain of R269 for hlGPDH. These results are consistent with a cryptic function for the similarly structured oxydianion binding domains of ScTIM, ScOMPDC and hlGPDH. Each enzyme utilizes the interactions with tetrahedral inorganic oxydianions to drive a conformational change that locks the substrate in a caged Michaelis complex that provides optimal stabilization of the different enzymatic transition states. The observation of dianion activation by stabilization of active caged Michaelis complexes may be generalized to the many other enzymes that utilize substrate binding energy to drive changes in enzyme conformation, which induce tight substrate fits.

Reflections on the catalytic power of a TIM-barrel

The TIM-barrel fold is described and its propagation throughout the enzyme universe noted. The functions of the individual front loops of the eponymous TIM-barrel of triosephosphate isomerase are presented in a discussion of: (a) electrophilic catalysis, by amino acid side chains from loops 1 and 4, of abstraction of an α-carbonyl hydrogen from substrate dihydroxyacetone phosphate (DHAP) or d-glyceraldehyde 3-phosphate (DGAP). (b) The engineering of loop 3 to give the monomeric variant monoTIM and the structure and catalytic properties of this monomer. (c) The interaction between loops 6, 7 and 8 and the phosphodianion of DHAP or DGAP. (d) The mechanism by which a ligand-gated conformational change, dominated by motion of loops 6 and 7, activates TIM for catalysis of deprotonation of DHAP or DGAP. (e) The conformational plasticity of TIM, and the utilization of substrate binding energy to "mold" the distorted active site loops of TIM mutants into catalytically active enzymes. The features of the TIM-barrel fold that favor effective protein catalysis are discussed.

Enzyme architecture: on the importance of being in a protein cage

Substrate binding occludes water from the active sites of many enzymes. There is a correlation between the burden to enzymatic catalysis of deprotonation of carbon acids and the substrate immobilization at solvent-occluded active sites for ketosteroid isomerase (KSI--small burden, substrate pKa=13), triosephosphate isomerase (TIM, substrate pKa≈18) and diaminopimelate epimerase (DAP epimerase, large burden, substrate pKa≈29) catalyzed reaction. KSI binds substrates at a surface cleft, TIM binds substrate at an exposed 'cage' formed by closure of flexible loops; and, DAP epimerase binds substrate in a tight cage formed by an 'oyster-like' clamping motion of protein domains. Directed evolution of a solvent-occluded active site at a designed protein catalyst of the Kemp elimination reaction is discussed.

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.

High-resolution crystal structure and redox properties of chloroplastic triosephosphate isomerase from Chlamydomonas reinhardtii

Triosephosphate isomerase (TPI) catalyzes the interconversion of glyceraldehyde-3-phosphate to dihydroxyacetone phosphate. Photosynthetic organisms generally contain two isoforms of TPI located in both cytoplasm and chloroplasts. While the cytoplasmic TPI is involved in the glycolysis, the chloroplastic isoform participates in the Calvin-Benson cycle, a key photosynthetic process responsible for carbon fixation. Compared with its cytoplasmic counterpart, the functional features of chloroplastic TPI have been poorly investigated and its three-dimensional structure has not been solved. Recently, several studies proposed TPI as a potential target of different redox modifications including dithiol/disulfide interchanges, glutathionylation, and nitrosylation. However, neither the effects on protein activity nor the molecular mechanisms underlying these redox modifications have been investigated. Here, we have produced recombinantly and purified TPI from the unicellular green alga Chlamydomonas reinhardtii (Cr). The biochemical properties of the enzyme were delineated and its crystallographic structure was determined at a resolution of 1.1 Å. CrTPI is a homodimer with subunits containing the typical (β/α)8-barrel fold. Although no evidence for TRX regulation was obtained, CrTPI was found to undergo glutathionylation by oxidized glutathione and trans-nitrosylation by nitrosoglutathione, confirming its sensitivity to multiple redox modifications.

Magnitude and origin of the enhanced basicity of the catalytic glutamate of triosephosphate isomerase

Glu-167 of triosephosphate isomerase from Trypanosoma brucei brucei (TbbTIM) acts as the base to deprotonate substrate to form an enediolate phosphate trianion intermediate. We report that there is a large ~6 pK unit increase in the basicity of the carboxylate side chain of Glu-167 upon binding of the inhibitor phosphoglycolate trianion (I(3-)), an analog of the enediolate phosphate intermediate, from pKEH ≈ 4 for the protonated free enzyme EH to pK(EHI) ≈ 10 for the protonated enzyme-inhibitor complex EH•I(3-). We propose that there is a similar increase in the basicity of this side chain when the physiological substrates are deprotonated by TbbTIM to form an enediolate phosphate trianion intermediate and that it makes an important contribution to the enzymatic rate acceleration. The affinity of wildtype TbbTIM for I(3-) increases 20,000-fold upon decreasing the pH from 9.3 to 4.9, because TbbTIM exists mainly in the basic form E over this pH range, while the inhibitor binds specifically to the rare protonated enzyme EH. This reflects the large increase in the basicity of the carboxylate side chain of Glu-167 upon binding of I(3-) to EH to give EH•I(3-). The I172A mutation at TbbTIM results in an ~100-fold decrease in the affinity of TbbTIM for I(3-) at pH < 6 and an ~2 pK unit decrease in the basicity of the carboxylate side chain of Glu-167 at the EH•I(3-) complex, to pK(EHI) = 7.7. Therefore, the hydrophobic side chain of Ile-172 plays a critical role in effecting the large increase in the basicity of the catalytic base upon the binding of substrate and/or inhibitors.

Crystal structures of triosephosphate isomerase from methicillin resistant Staphylococcus aureus MRSA252 provide structural insights into novel modes of ligand binding and unique conformations of catalytic loop

Staphylococcus aureus is one of the most dreaded pathogens worldwide and emergence of notorious antibiotic resistant strains have further exacerbated the present scenario. The glycolytic enzyme, triosephosphate isomerase (TIM) is one of the cell envelope proteins of the coccus and is involved in biofilm formation. It also plays an instrumental role in adherence and invasion of the bacteria into the host cell. To structurally characterize this important enzyme and analyze it's interaction with different inhibitors, substrate and transition state analogues, the present article describes several crystal structures of SaTIM alone and in complex with different ligands: glycerol-3-phosphate (G3P), glycerol-2-phosphate (G2P), 3-phosphoglyceric acid (3PG) and 2-phosphoglyceric acid (2PG). Unique conformations of the catalytic loop 6 (L6) has been observed in the different complexes. It is found to be in "almost closed" conformation in both subunits of the structure complexed to G3P. However L6 adopts the open conformation in presence of G2P and 2PG. The preference of the conformation of the catalytic loop can be correlated with the position of the phosphate group in the ligand. Novel modes of binding have been observed for G2P and 3PG for the very first time. The triose moiety is oriented away from the catalytic residues and occupies an entirely different position in some subunits. A completely new binding site for phosphate has also been identified in the complex with 2PG which differs substantially from the conventional phosphate binding site of the ligand in the crystal structures of TIM determined so far.

Mechanism for activation of triosephosphate isomerase by phosphite dianion: the role of a hydrophobic clamp

The role of the hydrophobic side chains of Ile-172 and Leu-232 in catalysis of the reversible isomerization of R-glyceraldehyde 3-phosphate (GAP) to dihydroxyacetone phosphate (DHAP) by triosephosphate isomerase (TIM) from Trypanosoma brucei brucei (Tbb) has been investigated. The I172A and L232A mutations result in 100- and 6-fold decreases in k(cat)/K(m) for the isomerization reaction, respectively. The effect of the mutations on the product distributions for the catalyzed reactions of GAP and of [1-(13)C]-glycolaldehyde ([1-(13)C]-GA) in D(2)O is reported. The 40% yield of DHAP from wild-type Tbb TIM-catalyzed isomerization of GAP with intramolecular transfer of hydrogen is found to decrease to 13% and to 4%, respectively, for the reactions catalyzed by the I172A and L232A mutants. Likewise, the 13% yield of [2-(13)C]-GA from isomerization of [1-(13)C]-GA in D(2)O is found to decrease to 2% and to 1%, respectively, for the reactions catalyzed by the I172A and L232A mutants. The decrease in the yield of the product of intramolecular transfer of hydrogen is consistent with a repositioning of groups at the active site that favors transfer of the substrate-derived hydrogen to the protein or the oxygen anion of the bound intermediate. The I172A and L232A mutations result in (a) a >10-fold decrease (I172A) and a 17-fold increase (L232A) in the second-order rate constant for the TIM-catalyzed reaction of [1-(13)C]-GA in D(2)O, (b) a 170-fold decrease (I172A) and 25-fold increase (L232A) in the third-order rate constant for phosphite dianion (HPO(3)(2-)) activation of the TIM-catalyzed reaction of GA in D(2)O, and (c) a 1.5-fold decrease (I172A) and a larger 16-fold decrease (L232A) in K(d) for activation of TIM by HPO(3)(2-) in D(2)O. The effects of the I172A mutation on the kinetic parameters for the wild-type TIM-catalyzed reactions of the whole substrate and substrate pieces are consistent with a decrease in the basicity of the carboxylate side chain of Glu-167 for the mutant enzyme. The data provide striking evidence that the L232A mutation leads to a ca. 1.7 kcal/mol stabilization of a catalytically active loop-closed form of TIM (E(C)) relative to an inactive open form (E(O)).

A paradigm for enzyme-catalyzed proton transfer at carbon: triosephosphate isomerase

Triosephosphate isomerase (TIM) catalyzes the stereospecific 1,2-proton shift at dihydroxyacetone phosphate (DHAP) to give (R)-glyceraldehyde 3-phosphate through a pair of isomeric enzyme-bound cis-enediolate phosphate intermediates. The chemical transformations that occur at the active site of TIM were well understood by the early 1990s. The mechanism for enzyme-catalyzed isomerization is similar to that for the nonenzymatic reaction in water, but the origin of the catalytic rate acceleration is not understood. We review the results of experimental work that show that a substantial fraction of the large 12 kcal/mol intrinsic binding energy of the nonreacting phosphodianion fragment of TIM is utilized to activate the active site side chains for catalysis of proton transfer. Evidence is presented that this activation is due to a phosphodianion-driven conformational change, the most dramatic feature of which is closure of loop 6 over the dianion. The kinetic data are interpreted within the framework of a model in which activation is due to the stabilization by the phosphodianion of a rare, desolvated, loop-closed form of TIM. The dianion binding energy is proposed to drive the otherwise thermodynamically unfavorable desolvation of the solvent-exposed active site. This reduces the effective local dielectric constant of the active site, to enhance stabilizing electrostatic interactions between polar groups and the anionic transition state, and increases the basicity of the carboxylate side chain of Glu-165 that functions to deprotonate the bound carbon acid substrate. A rebuttal is presented to the recent proposal [Samanta, M., Murthy, M. R. N., Balaram, H., and Balaram, P. (2011) ChemBioChem 12, 1886-1895] that the cationic side chain of K12 functions as an active site electrophile to protonate the carbonyl oxygen of DHAP.

Cloning, overexpression and characterization of Leishmania donovani triosephosphate isomerase

Triosephosphate isomerase (TIM) is a major enzyme in the glycolytic pathway, which catalyzes the interconversion of glyceraldehyde 3-phosphate to dihydroxyacetone phosphate. Here, we report cloning, expression and purification of a catalytically active recombinant TIM of Leishmania donovani (LdTIM). The recombinant LdTIM had a pH optimum in the range of 7.2-9.0, found stable at 25°C for 30 min and K(m) and V(max) for the substrate glyceraldehyde 3-phosphate was 0.328±0.02mM and 10.05mM/min/mg, respectively. The cysteine-reactive agent methylmethane thiosulphonate (MMTS) was used as probe, in order to test its effect on enzyme activity. The MMTS induced 75% enzyme inactivation within 15 min at 250 μM concentration. The biochemical characterization of LdTIM described in this work is the essential step towards deeper understanding of its role in parasite survival. The purification of LdTIM in bioactive form provides important tools for further functional and structural studies.

High resolution crystal structures of triosephosphate isomerase complexed with its suicide inhibitors: the conformational flexibility of the catalytic glutamate in its closed, liganded active site

The key residue of the active site of triosephosphate isomerase (TIM) is the catalytic glutamate, which is proposed to be important (i) as a catalytic base, for initiating the reaction, as well as (ii) for the subsequent proton shuttling steps. The structural properties of this glutamate in the liganded complex have been investigated by studying the high resolution crystal structures of typanosomal TIM, complexed with three suicide inhibitors: (S)-glycidol phosphate ((S)-GOP, at 0.99 Å resolution), (R)-glycidol phosphate, ((R)-GOP, at 1.08 Å resolution), and bromohydroxyacetone phosphate (BHAP, at 1.97 Å resolution). The structures show that in the (S)-GOP active site this catalytic glutamate is in the well characterized, competent conformation. However, an unusual side chain conformation is observed in the (R)-GOP and BHAP complexes. In addition, Glu97, salt bridged to the catalytic lysine in the competent active site, adopts an unusual side chain conformation in these two latter complexes. The higher chemical reactivity of (S)-GOP compared with (R)-GOP, as known from solution studies, can be understood: the structures indicate that in the case of (S)-GOP, Glu167 can attack the terminal carbon of the epoxide in a stereoelectronically favored, nearly linear OCO arrangement, but this is not possible for the (R)-GOP isomer. These structures confirm the previously proposed conformational flexibility of the catalytic glutamate in its closed, liganded state. The importance of this conformational flexibility for the proton shuttling steps in the TIM catalytic cycle, which is apparently achieved by a sliding motion of the side chain carboxylate group above the enediolate plane, is also discussed.

Determining the molecular mechanism of inactivation by chemical modification of triosephosphate isomerase from the human parasite Giardia lamblia: a study for antiparasitic drug design

Giardiasis, the most prevalent intestinal parasitosis in humans, is caused by Giardia lamblia. Current drug therapies have adverse effects on the host, and resistant strains against these drugs have been reported, demonstrating an urgent need to design more specific antigiardiasic drugs. ATP production in G. lamblia depends mainly on glycolysis; therefore, all enzymes of this pathway have been proposed as potential drug targets. We previously demonstrated that the glycolytic enzyme triosephosphate isomerase from G. lamblia (GlTIM), could be completely inactivated by low micromolar concentrations of thiol-reactive compounds, whereas, in the same conditions, the activity of human TIM (HuTIM) was almost unaltered. We found that the chemical modification (derivatization) of at least one Cys, of the five Cys residues per monomer in GlTIM, causes this inactivation. In this study, structural and functional studies were performed to describe the molecular mechanism of GlTIM inactivation by thiol-reactive compounds. We found that the Cys222 derivatization is responsible for GlTIM inactivation; this information is relevant because HuTIM has a Cys residue in an equivalent position (Cys217). GlTIM inactivation is associated with a decrease in ligand affinity, which affects the entropic component of ligand binding. In summary, this work describes a mechanism of inactivation that has not been previously reported for TIMs from other parasites and furthermore, we show that the difference in reactivity between the Cys222 in GlTIM and the Cys217 in HuTIM, indicates that the surrounding environment of each Cys residue has unique structural differences that can be exploited to design specific antigiardiasic drugs.

Wildtype and engineered monomeric triosephosphate isomerase from Trypanosoma brucei: partitioning of reaction intermediates in D2O and activation by phosphite dianion

Product yields for the reactions of (R)-glyceraldehyde 3-phosphate (GAP) in D2O at pD 7.9 catalyzed by wildtype triosephosphate isomerase from Trypanosoma brucei brucei (Tbb TIM) and a monomeric variant (monoTIM) of this wildtype enzyme were determined by (1)H NMR spectroscopy and were compared with the yields determined in earlier work for the reactions catalyzed by TIM from rabbit and chicken muscle [O'Donoghue, A. C., Amyes, T. L., and Richard, J. P. (2005), Biochemistry 44, 2610 - 2621]. Three products were observed from the reactions catalyzed by TIM: dihydroxyacetone phosphate (DHAP) from isomerization with intramolecular transfer of hydrogen, d-DHAP from isomerization with incorporation of deuterium from D2O into C-1 of DHAP, and d-GAP from incorporation of deuterium from D2O into C-2 of GAP. The yield of DHAP formed by intramolecular transfer of hydrogen decreases from 49% for the muscle enzymes to 40% for wildtype Tbb TIM to 34% for monoTIM. There is no significant difference in the ratio of the yields of d-DHAP and d-GAP for wildtype TIM from muscle sources and Trypanosoma brucei brucei, but partitioning of the enediolate intermediate of the monoTIM reaction to form d-DHAP is less favorable ((k(C1))(D)/(k(C2))(D) = 1.1) than for the wildtype enzyme ((k(C1))(D)/(k(C2))(D) = 1.7). Product yields for the wildtype Tbb TIM and monoTIM-catalyzed reactions of glycolaldehyde labeled with carbon-13 at the carbonyl carbon ([1-(13)C]-GA) at pD 7.0 in the presence of phosphite dianion and in its absence were determined by (1)H NMR spectroscopy [Go, M. K., Amyes, T. L., and Richard, J. P. (2009) Biochemistry 48, 5769-5778]. There is no detectable difference in the yields of the products of wildtype muscle and Tbb TIM-catalyzed reactions of [1-(13)C]-GA in D2O. The kinetic parameters for phosphite dianion activation of the reactions of [1-(13)C]-GA catalyzed by wildtype Tbb TIM are similar to those reported for the enzyme from rabbit muscle [Amyes, T. L. and Richard, J. P. (2007) Biochemistry 46, 5841-5854], but there is no detectable dianion activation of the reaction catalyzed by monoTIM. The engineered disruption of subunit contacts at monoTIM causes movement of the essential side chains of Lys-13 and His-95 away from the catalytic active positions. We suggest that this places an increased demand that the intrinsic binding energy of phosphite dianion be utilized to drive the change in the conformation of monoTIM back to the active structure for wildtype TIM.

A role for flexible loops in enzyme catalysis

Triosephosphate isomerase (TIM), glycerol 3-phosphate dehydrogenase, and orotidine 5'-monophosphate decarboxylase each use the binding energy from the interaction of phosphite dianion with a flexible phosphate gripper loop to activate a second, phosphodianion-truncated, substrate towards enzyme-catalyzed proton transfer, hydride transfer, and decarboxylation, respectively. Studies on TIM suggest that the most important general effect of loop closure over the substrate phosphodianion, and the associated conformational changes, is to extrude water from the enzyme active site. This should cause a decrease in the effective active-site dielectric constant, and an increase in transition state stabilization from enhanced electrostatic interactions with polar amino acid side chains. The most important specific effect of these conformational changes is to increase the basicity of the carboxylate side chain of the active site glutamate base by its placement in a 'hydrophobic cage'.

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.

Crystallographic binding studies with an engineered monomeric variant of triosephosphate isomerase

Crystallographic binding studies have been carried out to probe the active-site binding properties of a monomeric variant (A-TIM) of triosephosphate isomerase (TIM). These binding studies are part of a structure-based directed-evolution project aimed towards changing the substrate specificity of monomeric TIM and are therefore aimed at finding binders which are substrate-like molecules. A-TIM has a modified more extended binding pocket between loop-7 and loop-8 compared with wild-type TIM. The A-TIM crystals were grown in the presence of citrate, which is bound in the active site of each of the two molecules in the asymmetric unit. In this complex, the active-site loops loop-6 and loop-7 adopt the closed conformation, similar to that observed in liganded wild-type TIM. Extensive crystal-soaking protocols have been developed to flush the bound citrate out of the active-site pocket of both molecules and the crystal structure shows that the unliganded open conformation of the A-TIM active site is the same as in unliganded wild-type TIM. It is also shown that sulfonate compounds corresponding to the transition-state analogue 2-phosphoglycolate bind in the active site, which has a closed conformation. It is also shown that the new binding pocket of A-TIM can bind 3-phosphoglycerate (3PGA; an analogue of a C4-sugar phosphate) and 4-phospho-D-erythronohydroxamic acid (4PEH; an analogue of a C5-sugar phosphate). Therefore, these studies have provided a rationale for starting directed-evolution experiments aimed at generating the catalytic properties of a C5-sugar phosphate isomerase on the A-TIM framework.

Triosephosphate isomerase: 15N and 13C chemical shift assignments and conformational change upon ligand binding by magic-angle spinning solid-state NMR spectroscopy

Microcrystalline uniformly (13)C,(15)N-enriched yeast triosephosphate isomerase (TIM) is sequentially assigned by high-resolution solid-state NMR (SSNMR). Assignments are based on intraresidue and interresidue correlations, using dipolar polarization transfer methods, and guided by solution NMR assignments of the same protein. We obtained information on most of the active-site residues involved in chemistry, including some that were not reported in a previous solution NMR study, such as the side-chain carbons of His95. Chemical shift differences comparing the microcrystalline environment to the aqueous environment appear to be mainly due to crystal packing interactions. Site-specific perturbations of the enzyme's chemical shifts upon ligand binding are studied by SSNMR for the first time. These changes monitor proteinwide conformational adjustment upon ligand binding, including many of the sites probed by solution NMR and X-ray studies. Changes in Gln119, Ala163, and Gly210 were observed in our SSNMR studies, but were not reported in solution NMR studies (chicken or yeast). These studies identify a number of new sites with particularly clear markers for ligand binding, paving the way for future studies of triosephosphate isomerase dynamics and mechanism.

An internal sequence targets Trypanosoma brucei triosephosphate isomerase to glycosomes

In kinetoplastid protists, glycolysis is compartmentalized in glycosomes, organelles belonging to the peroxisome family. The Trypanosoma brucei glycosomal enzyme triosephosphate isomerase (TPI) does not contain either of the two established peroxisome-targeting signals, but we identified a 22 amino acids long fragment, present at an internal position of the polypeptide, that has the capacity to route a reporter protein to glycosomes in transfected trypanosomes, as demonstrated by cell-fractionation experiments and corroborating immunofluorescence studies. This polypeptide-internal routing information seems to be unique for the sequence of the trypanosome enzyme: a reporter protein fused to a Saccharomyces cerevisiae peptide containing the sequence corresponding to the 22-residue fragment of the T. brucei enzyme, was not targeted to glycosomes. In yeasts, as in most other organisms, TPI is indeed exclusively present in the cytosol. These results suggest that it may be possible to develop new trypanocidal drugs by targeting specifically the glycosome import mechanism of TPI.

Specific sequence combinations at parallel and antiparallel helix-helix interfaces

Orientation of helices at parallel and antiparallel helix-helix interfaces in proteins depends on interacting amino acids from both helices. Particularly important are amino acids at positions analogous to a and d in GCN4 leucine zipper nomenclature, which form hydrophobic core. In this work repeating sequence combinations at a and d positions characteristic for both parallel and antiparallel packing are shown. Layer packing of hydrophobic groups is compared for possible combinations of aliphatic amino acids at all four positions. Correlation between specific position of methyl groups and interhelical angle is found for parallel and antiparallel types of packing.

Perturbation of the dimer interface of triosephosphate isomerase and its effect on Trypanosoma cruzi

BACKGROUND

Chagas disease affects around 18 million people in the American continent. Unfortunately, there is no satisfactory treatment for the disease. The drugs currently used are not specific and exert serious toxic effects. Thus, there is an urgent need for drugs that are effective. Looking for molecules to eliminate the parasite, we have targeted a central enzyme of the glycolytic pathway: triosephosphate isomerase (TIM). The homodimeric enzyme is catalytically active only as a dimer. Because there are significant differences in the interface of the enzymes from the parasite and humans, we searched for small molecules that specifically disrupt contact between the two subunits of the enzyme from Trypanosoma cruzi but not those of TIM from Homo sapiens (HTIM), and tested if they kill the parasite.

METHODOLOGY/PRINCIPAL FINDINGS

Dithiodianiline (DTDA) at nanomolar concentrations completely inactivates recombinant TIM of T. cruzi (TcTIM). It also inactivated HTIM, but at concentrations around 400 times higher. DTDA was also tested on four TcTIM mutants with each of its four cysteines replaced with either valine or alanine. The sensitivity of the mutants to DTDA was markedly similar to that of the wild type. The crystal structure of the TcTIM soaked in DTDA at 2.15 A resolution, and the data on the mutants showed that inactivation resulted from alterations of the dimer interface. DTDA also prevented the growth of Escherichia coli cells transformed with TcTIM, had no effect on normal E. coli, and also killed T. cruzi epimastigotes in culture.

CONCLUSIONS/SIGNIFICANCE

By targeting on the dimer interface of oligomeric enzymes from parasites, it is possible to discover small molecules that selectively thwart the life of the parasite. Also, the conformational changes that DTDA induces in the dimer interface of the trypanosomal enzyme are unique and identify a region of the interface that could be targeted for drug discovery.

Outliers in SAR and QSAR: 2. Is a flexible binding site a possible source of outliers?

Structure-activity relationship (SAR) and/or quantitative structure-activity relationship (QSAR) studies play an important role in a lead optimization of drug discovery research. When there is a lack of ligand-bound protein structural information, one of the assumptions in SAR and QSAR studies is that similar analogs bind to the same binding site in a similar binding mode. In such studies, outliers have often been observed, especially in QSAR. However, most of these studies have focused their attention on the development of QSAR and left outliers unattended. We searched ligand-bound X-ray crystal structures from the protein structure database to find evidences that could indicate a possible source of outliers in SAR or QSAR. Our results showed the possibility of conformational changes in a flexible binding site as one possible source of outliers.

Enzymatic catalysis of proton transfer at carbon: activation of triosephosphate isomerase by phosphite dianion

More than 80% of the rate acceleration for enzymatic catalysis of the aldose-ketose isomerization of (R)-glyceraldehyde 3-phosphate (GAP) by triosephosphate isomerase (TIM) can be attributed to the phosphodianion group of GAP [Amyes, T. L., O'Donoghue, A. C., and Richard, J. P. (2001) J. Am. Chem. Soc. 123, 11325-11326]. We examine here the necessity of the covalent connection between the phosphodianion and triose sugar portions of the substrate by "carving up" GAP into the minimal neutral two-carbon sugar glycolaldehyde and phosphite dianion pieces. This "two-part substrate" preserves both the alpha-hydroxycarbonyl and oxydianion portions of GAP. TIM catalyzes proton transfer from glycolaldehyde in D2O, resulting in deuterium incorporation that can be monitored by 1H NMR spectroscopy, with kcat/Km = 0.26 M-1 s-1. Exogenous phosphite dianion results in a very large increase in the observed second-order rate constant (kcat/Km)obsd for turnover of glycolaldehyde, and the dependence of (kcat/Km)obsd on [HPO32-] exhibits saturation. The data give kcat/Km = 185 M-1 s-1 for turnover of glycolaldehyde by TIM that is saturated with phosphite dianion so that the separate binding of phosphite dianion to TIM results in a 700-fold acceleration of proton transfer from carbon. The binding of phosphite dianion to the free enzyme (Kd = 38 mM) is 700-fold weaker than its binding to the fleeting complex of TIM with the altered substrate in the transition state (Kd = 53 muM); the total intrinsic binding energy of phosphite dianion in the transition state is 5.8 kcal/mol. We propose a physical model for catalysis by TIM in which the intrinsic binding energy of the substrate phosphodianion group is utilized to drive closing of the "mobile loop" and a protein conformational change that leads to formation of an active site environment that is optimally organized for stabilization of the transition state for proton transfer from alpha-carbonyl carbon.

Amino acid composition of parallel helix-helix interfaces

Amino acids at helix-helix parallel interfaces influence arrangement of helices and interhelical angles. Parallel interfaces in 79 proteins were considered. Location of amino acids at the positions analogous to a and d in GCN4 leucine zipper nomenclature shows that certain combinations of amino acids characteristic for parallel packing occur more often than could be expected by chance. Repeating sequence combinations occur at a and d positions of parallel helix-helix interfaces with similar values of interhelical angles not only in homologous proteins but also within the same protein and in nonhomologous proteins. Within each group of observed combinations correlation exists between the size of amino acid and magnitude of the interhelical angle.

Functional role of the conserved active site proline of triosephosphate isomerase

The importance of the fully conserved active site proline, Pro168, for the reaction mechanism of triosephosphate isomerase (TIM) has been investigated by studying the enzymatic and crystallographic properties of the P168A variant of trypanosomal TIM. In TIM, Pro168 follows the key catalytic residue Glu167, situated at the beginning of the flexible active site loop (loop 6). Turnover numbers of the P168A variant for its substrates are reduced approximately 50-fold, whereas the Km values are approximately 2 times lower. The affinity of the P168A variant for the transition state analogue 2-phosphoglycolate (2PG) is reduced 5-fold. The crystal structures of unliganded and liganded (2PG) P168A show that the phosphate moiety of 2PG is bound similarly as in wild-type TIM, whereas the interactions of the carboxylic acid moiety with the side chain of the catalytic Glu167 differ. The unique properties of the proline side chain at position 168 are required to transmit ligand binding to the conformational change of Glu167: the side chain of Glu167 flips from the inactive swung-out to the active swung-in conformation on ligand binding in wild-type TIM, whereas in the mutant this conformational change does not occur. Further structural comparisons show that in the wild-type enzyme the concerted movement of loop 6 and loop 7 from unliganded-open to liganded-closed appears to be facilitated by the interactions of the phosphate moiety with loop 7. Apparently, the rotation of 90 degrees of the Gly211-Gly212 peptide plane of loop 7 plays a key role in this concerted movement.

Disulfide bridges in the mesophilic triosephosphate isomerase from Giardia lamblia are related to oligomerization and activity

Triosephosphate isomerase from the mesophile Giardia lamblia (GlTIM) is the only known TIM with natural disulfide bridges. We previously found that oxidized and reduced thiol states of GlTIM are involved in the interconversion between native dimers and higher oligomeric species, and in the regulation of enzymatic activity. Here, we found that trophozoites and cysts have different oligomeric species of GlTIM and complexes of GlTIM with other proteins. Our data indicate that the internal milieu of G. lamblia is favorable for the formation of disulfide bonds. Enzyme mutants of the three most solvent exposed Cys of GlTIM (C202A, C222A, and C228A) were prepared to ascertain their contribution to oligomerization and activity. The data show that the establishment of a disulfide bridge between two C202 of two dimeric GlTIMs accounts for multimerization. In addition, we found that the establishment of an intramonomeric disulfide bond between C222 and C228 abolishes catalysis. Multimerization and inactivation are both reversed by reducing conditions. The 3D structure of the C202A GlTIM was solved at 2.1 A resolution, showing that the environment of the C202 is prone to hydrophobic interactions. Molecular dynamics of an in silico model of GlTIM when the intramonomeric disulfide bond is formed, showed that S216 is displaced 4.6 A from its original position, causing loss of hydrogen bonds with residues of the active-site loop. This suggests that this change perturb the conformational state that aligns the catalytic center with the substrate, inducing enzyme inactivation.

The planar conformation of a strained proline ring: a QM/MM study

QM and QM/MM energy calculations have been carried out on an atomic resolution structure of liganded triosephosphate isomerase (TIM) that has an active site proline (Pro168) in a planar conformation. The origin of the planarity of this proline has been identified. Steric interactions between the atoms of the proline ring and a tyrosine ring (Tyr166) on one side of the proline prevent the ring from adopting the up pucker (chi1 is approximately -30 degrees), while the side chain of a nearby alanine (Ala171) forbids the down pucker (chi1 is approximately +30 degrees). To obtain a proline conformation that is in agreement with the experimentally observed planar state, a quantum system of sufficient size is required and should at least include the nearby side chains of Tyr166, Ala171, and Glu129 to provide enough stabilization. It is argued that the current force fields for structure optimization do not describe strained protein fragments correctly. The proline is part of a catalytic loop that closes upon ligand binding. Comparison of the proline conformation in different TIM X-ray structures, indicates that in the closed conformation of TIM the proline is planar or nearly planar, while in the open conformation it is down puckered. This suggests that the planarity possibly plays a role in the overall catalytic cycle of TIM, presumable acting as a reservoir of energy that becomes available upon loop opening.

Drosophila model of human inherited triosephosphate isomerase deficiency glycolytic enzymopathy

Heritable mutations, known as inborn errors of metabolism, cause numerous devastating human diseases, typically as a result of a deficiency in essential metabolic products or the accumulation of toxic intermediates. We have isolated a missense mutation in the Drosophila sugarkill (sgk) gene that causes phenotypes analogous to symptoms of triosephosphate isomerase (TPI) deficiency, a human familial disease, characterized by anaerobic metabolic dysfunction resulting from pathological missense mutations affecting the encoded TPI protein. In Drosophila, the sgk gene encodes the glycolytic enzyme TPI. Our analysis of sgk mutants revealed TPI impairment associated with reduced longevity, progressive locomotor deficiency, and neural degeneration. Biochemical studies demonstrate that mutation of this glycolytic enzyme gene does not result in a bioenergetic deficit, suggesting an alternate cause of enzymopathy associated with TPI impairment.

Solution NMR and Computer Simulation Studies of Active Site Loop Motion in Triosephosphate Isomerase†

Solution NMR spin relaxation experiments and classical MD simulations are used to study the dynamics of triosephosphate isomerase (TIM) in complex with glycerol 3-phosphate (G3P). Three regions in TIM exhibit conformational transitions on the micros-ms time scale as detected by chemical exchange broadening effects in NMR spectroscopy: residue Lys 84 on helix C, located at the dimeric interface; active site loop 6; and helix G. The results indicate that the conformational exchange process affecting the residues of loop 6 is the correlated opening and closing of the loop. Distinct processes are responsible for the chemical exchange linebroadening observed in the other regions of TIM. MD simulations confirm that motions of individual residues within the active site loop are correlated and suggest that the chemical exchange processes observed for residues in helix G arise from transitions between 3(10)- and alpha-helical structures. The results of the joint NMR and MD study provide global insight into the role of conformational dynamic processes in the function of TIM.

Structural differences in triosephosphate isomerase from different species and discovery of a multitrypanosomatid inhibitor

We examined the interfaces of homodimeric triosephosphate isomerase (TIM) from eight different species. The crystal structures of the enzymes showed that a portion of the interface is markedly similar in TIMs from Trypanosoma cruzi (TcTIM), Trypanosoma brucei, and Leishmania mexicana and significantly different from that of TIMs from human, yeast, chicken, Plasmodium falciparum, and Entamoeba histolytica. Since this interfacial region is central in the stability of TcTIM, we hypothesized that it would be possible to find agents that selectively affect the stability of TIMs from the three trypanosomatids. We found that 6,6'-bisbenzothiazole-2,2' diamine in the low micromolar range causes a desirable irreversible inactivation of the enzymes from the three trypanosomatids and has no effect on the other five TIMs. Thus, the data indicate that it is possible to find compounds that induce selective inactivation of the enzymes from three different trypanosomatids.

Reversible Equilibrium Unfolding of Triosephosphate Isomerase from Trypanosoma cruzi in Guanidinium Hydrochloride Involves Stable Dimeric and Monomeric Intermediates

The reversible guanidinium hydrochloride-induced unfolding of Trypanosoma cruzi triosephosphate isomerase (TcTIM) was characterized under equilibrium conditions. The catalytic activity was followed as a native homodimeric functional probe. Circular dichroism, intrinsic fluorescence, and size-exclusion chromatography were used as secondary, tertiary, and quaternary structural probes, respectively. The change in ANS fluorescence intensity with increasing denaturant concentrations was also determined. The results show that two stable intermediates exist in the transition from the homodimeric native enzyme to the unfolded monomers: one (N(2*)) is a slightly more expanded, non-native, and active dimer, and the other is a partially expanded monomer (M) that binds ANS. Spectroscopic and activity data were used to reach a thermodynamic characterization. The results indicate that the Gibbs free energies for the partial reactions are 4.5 (N(2) <==> N(2*)), 65.8 (N(2*) <==> 2M), and 17.8 kJ/mol (M <==> U). It appears that TcTIM monomers are more stable than those found for other TIM species (except yeast TIM), where monomer stability is only marginal. These results are compared with those for the guanidinium hydrochloride-induced denaturation of TIM from different species, where despite the functional and three-dimensional similarities, a remarkable heterogeneity exists in the unfolding pathways.

Biogenesis of peroxisomes and glycosomes: trypanosomatid glycosome assembly is a promising new drug target

In trypanosomatids (Trypanosoma and Leishmania), protozoa responsible for serious diseases of mankind in tropical and subtropical countries, core carbohydrate metabolism including glycolysis is compartmentalized in peculiar peroxisomes called glycosomes. Proper biogenesis of these organelles and the correct sequestering of glycolytic enzymes are essential to these parasites. Biogenesis of glycosomes in trypanosomatids and that of peroxisomes in other eukaryotes, including the human host, occur via homologous processes involving proteins called peroxins, which exert their function through multiple, transient interactions with each other. Decreased expression of peroxins leads to death of trypanosomes. Peroxins show only a low level of sequence conservation. Therefore, it seems feasible to design compounds that will prevent interactions of proteins involved in biogenesis of trypanosomatid glycosomes without interfering with peroxisome formation in the human host cells. Such compounds would be suitable as lead drugs against trypanosomatid-borne diseases.

A molecular modeling and QSAR study of suppressors of the growth of Trypanosoma cruzi epimastigotes

In this work, we have used molecular modeling and QSAR tools to study 18 dithiocarbamate suppressors of the growth of Trypanosoma cruzi epimastigotes, which have been reported in the literature as superoxide dismutase (SOD) inhibitors. The principal component analysis (PCA) showed that the descriptors superficial area, heat of formation, logarithm of the partition coefficient, charge of the nitrogen atom from the dithiocarbamate group and Charges of the two carbon atoms adjacent to that nitrogen are responsible for the classification between the higher and lower trypanomicid activity. Using multiple linear regression (MLR) and docking methods it was possible to identify the probable bioactive isomers that suppress of the growth of T. cruzi epimastigotes. Our best partial least square (PLS) model obtained with these six descriptors yields a good correlation between experimental and predicted biological activities and compares two different SODs as possible target for interaction with the dithiocarbamates.