Anion binding at the active site of trypanosomal triosephosphate isomerase. Monohydrogen phosphate does not mimic sulphate

Hydrophobicity-Based Force Field In Enzymes

The biocatalysis process takes place with the participation of enzymes, which, depending on the reaction carried out, require, apart from the appropriate arrangement of catalytic residues, an appropriate external force field. It is generated by the protein body. The relatively small size of the part directly involved in the process itself is supported by the presence of an often complex structure of the protein body, the purpose of which is to provide an appropriate local force field, eliminating the influence of water. Very often, the large size of the enzyme is an expression of the complex form of this field. In this paper, a comparative analysis of arbitrarily selected enzymes, representatives of different enzyme classes, was carried out, focusing on the measurement of the diversity of the force field provided by a given protein. This analysis was based on the fuzzy oil drop model (FOD) and its modified version (FOD-M), which takes into account the participation of nonaqueous external factors in shaping the structure and thus the force field within the protein. The degree and type of ordering of the hydrophobicity distribution in the protein molecule is the result of the influence of the environment but also the supplier of the local environment for a given process, including the catalysis process in particular. Determining the share of a nonaqueous environment is important due to the ubiquity of polar water, whose participation in processes with high specificity requires control. It can be assumed that some enzymes in their composition have a permanently built-in part, the role of which is reduced to that of a permanent chaperone. It provides a specific external force field needed for the process. The proposed model, generalized to other types of proteins, may also provide a form of recording the environment model for the simulation of the in silico protein folding process, taking into account the impact of its differentiation.

Active-Site Glu165 Activation in Triosephosphate Isomerase and Its Deprotonation Kinetics

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate (DHAP) and d-glyceraldehyde 3-phosphate (GAP) via an enediol(ate) intermediate. The active-site residue Glu165 serves as the catalytic base during catalysis. It abstracts a proton from C1 carbon of DHAP to form the reaction intermediate and donates a proton to C2 carbon of the intermediate to form product GAP. Our difference Fourier transform infrared spectroscopy studies on the yeast TIM (YeTIM)/phosphate complex revealed a C═O stretch band at 1706 cm^-1 from the protonated Glu165 carboxyl group at pH 7.5, indicating that the p K_a of the catalytic base is increased by >3.0 pH units upon phosphate binding, and that the Glu165 carboxyl environment in the complex is still hydrophilic in spite of the increased p K_a. Hence, the results show that the binding of the phosphodianion group is part of the activation mechanism which involves the p K_a elevation of the catalytic base Glu165. The deprotonation kinetics of Glu165 in the μs to ms time range were determined via infrared (IR) T-jump studies on the YeTIM/phosphate and ("heavy enzyme") [U-¹³C,-¹⁵N]YeTIM/phosphate complexes. The slower deprotonation kinetics in the ms time scale is due to phosphate dissociation modulated by the loop motion, which slows down by enzyme mass increase to show a normal heavy enzyme kinetic isotope effect (KIE) ∼1.2 (i.e., slower rate in the heavy enzyme). The faster deprotonation kinetics in the tens of μs time scale is assigned to temperature-induced p K_a decrease, while phosphate is still bound, and it shows an inverse heavy enzyme KIE ∼0.89 (faster rate in the heavy enzyme). The IR static and T-jump spectroscopy provides atomic-level resolution of the catalytic mechanism because of its ability to directly observe the bond breaking/forming process.

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.

The ‘CαNN’ motif: an intrinsic lover of sulfate and phosphate ions

The ‘C^αNN’ motif has an intrinsic affinity for the anions and can recognize anion through local interactions along with augmentation of the helical conformation at the motif segment.

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.

Sulfate ion interaction with 'anion recognition' short peptide motif at the N-terminus of an isolated helix: A conformational landscape

Anion-binding motifs in proteins are typically conserved in sequence and conformation. Crystal structural studies have shown that such motifs often occur in loop regions preceding a helix and interaction with the anions can induce their well defined conformational changes. In order to understand the properties of such motifs in isolation, we have synthesized an 18-residue chimeric polypeptide whose C-terminal part is a designed helix and its N-terminal consists of a C(alpha)NN anion binding structural motif containing residues Leu-Gly-Lys-Gln (residues 107-110 of protein DNA-glycosylase). We present evidence for the interaction of a sulfate (SO(4)(2-)) ion with the L-G-K-Q segment using complementary spectroscopic techniques. Moreover, upon interaction with SO(4)(2-) ion the N-terminal L-G-K-Q segment undergoes a non-helical to helical transition similar to what is observed in protein crystal structure. This work clearly demonstrates the "local" nature of anion binding and the accompanying conformational change that helps in understanding the influence of sequence/structural context of anion binding in proteins.

Dimeric quaternary structure of the prototypical dual specificity phosphatase VH1

The Vaccinia virus H1 gene product, VH1, is a dual specificity phosphatase that down-regulates the cellular antiviral response by dephosphorylating STAT1. The crystal structure of VH1, determined at 1.32 A resolution, reveals a novel dimeric quaternary structure, which exposes two active sites spaced approximately 39 A away from each other. VH1 forms a stable dimer via an extensive domain swap of the N-terminal helix (residues 1-20). In vitro, VH1 can dephosphorylate activated STAT1, in a reaction that is competed by the nuclear transport adapter importin alpha5. Interestingly, VH1 is inactive with respect to STAT1 bound to DNA, suggesting that the viral phosphatase acts predominantly on the cytoplasmic pool of activated STAT1. We propose that the dimeric quaternary structure of VH1 is essential for specific recognition of activated STAT1, which prevents its nuclear translocation, thus blocking interferon-gamma signal transduction and antiviral response.

Kinetic and structural properties of triosephosphate isomerase from Helicobacter pylori

Triosephosphate isomerase (TIM) catalyzes the interconversion between dihydroxyacetone phosphate and D-glyceraldehyde-3-phosphate in the glycolysis-gluconeogenesis metabolism pathway. The Helicobacter pylori TIM gene (HpTIM) was cloned, and HpTIM was expressed and purified. The enzymatic activity of HpTIM for the substrate GAP was determined (K(m) = 3.46 +/- 0.23 mM and k(cat) = 8.8 x 10(4) min(-1)). The crystal structure of HpTIM was determined by molecular replacement at 2.3 A resolution. The overall structure of HpTIM was (beta/alpha)beta(beta/alpha)(6), which resembles the common TIM barrel fold, (beta/alpha)(8); however, a helix is missing after the second beta-strand. The conformation of loop 6 and binding of phosphate ion suggest that the determined structure of HpTIM was in the "closed" state. A highly conserved Arg-Asp salt bridge in the "DX(D/N)G" motif of most TIMs is absent in HpTIM because the sequence of this motif is "(211)SVDG(214)." To determine the significance of this salt bridge to HpTIM, four mutants, including K183S, K183A, D213Q, and D213A, were constructed and characterized. The results suggest that this conserved salt bridge is not essential for the enzymatic activity of HpTIM; however, it might contribute to the conformational stability of HpTIM.

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.

Active site loop motion in triosephosphate isomerase: T-jump relaxation spectroscopy of thermal activation

As for many enzymes, the enzymatic pathway of triosephosphate isomerase (TIM) includes the partially rate determining motion of an active site loop (loop 6, residues 166-176), which must remain closed during chemistry but must open just before product release. The motion of this loop was monitored using laser induced temperature-jump relaxation spectroscopy at nanosecond to millisecond time resolution. Trp168 in the hinge of the mobile loop served as a fluorophore reporter in a mutant of the yeast enzyme. The opening rate was studied as a function of the concentration of glycerol 3-phosphate, a substrate surrogate. Monoexponential kinetics were observed; assuming a simple two-step ligand release mechanism involving an encounter complex intermediate, the time scales of loop opening and closing were derived. The opening rate of the loop at 25 degrees C was determined to be 2500 +/- 1000 s(-1), in remarkable agreement with solution and solid state NMR measurements. The closing rate at the same temperature was 46,700 +/- 1800 s(-1). The rates were also studied as a function of the sample temperature following the jump. Enthalpies of activation of the loop motion, DeltaH(close) and DeltaH(open), were estimated to be 13.8 and 14.1 kcal/mol, respectively. The enthalpy of dissociation estimated from the kinetic studies is in reasonable agreement with steady-state values. Moreover, the enthalpy was dissected, for the first time, into components associated with ion binding and with protein conformational change. The enthalpy of the release reaction appeared to have a substantial contribution from the dissociation of the ligand from the encounter complex, found to be endothermic at 6 kcal/mol. In contrast, the population ratio of the open to closed loop conformations is found to favor the closed conformation but to be substantially less temperature dependent than the release step. Preliminary data of other ligands show that G3P behavior resembles that of the substrate but differs from 2-phosphoglycolate, a tight binding inhibitor, and phosphate. This study represents one of the first detailed comparisons between NMR and fluorescence based probes of protein motion and results in good agreement between the methods. The data in aggregate support a model in which the rate of the loop opening for TIM is dependent on the ligand and results in opening rates in the presence of the product that are comparable to enzymatic throughput, kcat.

Prediction of binding constants of protein ligands: a fast method for the prioritization of hits obtained from de novo design or 3D database search programs

A dataset of 82 protein-ligand complexes of known 3D structure and binding constant Ki was analysed to elucidate the important factors that determine the strength of protein-ligand interactions. The following parameters were investigated: the number and geometry of hydrogen bonds and ionic interactions between the protein and the ligand, the size of the lipophilic contact surface, the flexibility of the ligand, the electrostatic potential in the binding site, water molecules in the binding site, cavities along the protein-ligand interface and specific interactions between aromatic rings. Based on these parameters, a new empirical scoring function is presented that estimates the free energy of binding for a protein-ligand complex of known 3D structure. The function distinguishes between buried and solvent accessible hydrogen bonds. It tolerates deviations in the hydrogen bond geometry of up to 0.25 A in the length and up to 30 degrees in the hydrogen bond angle without penalizing the score. The new energy function reproduces the binding constants (ranging from 3.7 x 10(-2) M to 1 x 10(-14) M, corresponding to binding energies between -8 and -80 kJ/mol) of the dataset with a standard deviation of 7.3 kJ/mol corresponding to 1.3 orders of magnitude in binding affinity. The function can be evaluated very fast and is therefore also suitable for the application in a 3D database search or de novo ligand design program such as LUDI. The physical significance of the individual contributions is discussed.

Protein crystallography and infectious diseases

The current rapid growth in the number of known 3-dimensional protein structures is producing a database of structures that is increasingly useful as a starting point for the development of new medically relevant molecules such as drugs, therapeutic proteins, and vaccines. This development is beautifully illustrated in the recent book, Protein structure: New approaches to disease and therapy (Perutz, 1992). There is a great and growing promise for the design of molecules for the treatment or prevention of a wide variety of diseases, an endeavor made possible by the insights derived from the structure and function of crucial proteins from pathogenic organisms and from man. We present here 2 illustrations of structure-based drug design. The first is the prospect of developing antitrypanosomal drugs based on crystallographic, ligand-binding, and molecular modeling studies of glycolytic glycosomal enzymes from Trypanosomatidae. These unicellular organisms are responsible for several tropical diseases, including African and American trypanosomiases, as well as various forms of leishmaniasis. Because the target enzymes are also present in the human host, this project is a pioneering study in selective design. The second illustrative case is the prospect of designing anti-cholera drugs based on detailed analysis of the structure of cholera toxin and the closely related Escherichia coli heat-labile enterotoxin. Such potential drugs can be targeted either at inhibiting the toxin's receptor binding site or at blocking the toxin's intracellular catalytic activity. Study of the Vibrio cholerae and E. coli toxins serves at the same time as an example of a general approach to structure-based vaccine design. These toxins exhibit a remarkable ability to stimulate the mucosal immune system, and early results have suggested that this property can be maintained by engineered fusion proteins based on the native toxin structure. The challenge is thus to incorporate selected epitopes from foreign pathogens into the native framework of the toxin such that crucial features of both the epitope and the toxin are maintained. That is, the modified toxin must continue to evoke a strong mucosal immune response, and this response must be directed against an epitope conformation characteristic of the original pathogen.

Crystal structure of recombinant human triosephosphate isomerase at 2.8 A resolution. Triosephosphate isomerase-related human genetic disorders and comparison with the trypanosomal enzyme

The crystal structure of recombinant human triosephosphate isomerase (hTIM) has been determined complexed with the transition-state analogue 2-phosphoglycolate at a resolution of 2.8 A. After refinement, the R-factor is 16.7% with good geometry. The asymmetric unit contains 1 complete dimer of 53,000 Da, with only 1 of the subunits binding the inhibitor. The so-called flexible loop, comprising residues 168-174, is in its "closed" conformation in the subunit that binds the inhibitor, and in the "open" conformation in the other subunit. The tips of the loop in these 2 conformations differ up to 7 A in position. The RMS difference between hTIM and the enzyme of Trypanosoma brucei, the causative agent of sleeping sickness, is 1.12 A for 487 C alpha positions with 53% sequence identity. Significant sequence differences between the human and parasite enzymes occur at about 13 A from the phosphate binding site. The chicken and human enzymes have an RMS difference of 0.69 A for 484 equivalent residues and about 90% sequence identity. Complementary mutations ensure a great similarity in the packing of side chains in the core of the beta-barrels of these 2 enzymes. Three point mutations in hTIM have been correlated with severe genetic disorders ranging from hemolytic disorder to neuromuscular impairment. Knowledge of the structure of the human enzyme provides insight into the probable effect of 2 of these mutations, Glu 104 to Asp and Phe 240 to Ile, on the enzyme. The third mutation reported to be responsible for a genetic disorder, Gly 122 to Arg, is however difficult to explain. This residue is far away from both catalytic centers in the dimer, as well as from the dimer interface, and seems unlikely to affect stability or activity. Inspection of the 3-dimensional structure of trypanosomal triosephosphate isomerase, which has a methionine at position 122, only increased the mystery of the effects of the Gly to Arg mutation in the human enzyme.

Structures of the "open" and "closed" state of trypanosomal triosephosphate isomerase, as observed in a new crystal form: implications for the reaction mechanism

The structure of trypanosomal triosephosphate isomerase (TIM) has been solved at a resolution of 2.1A in a new crystal form grown at pH 8.8 from PEG6000. In this new crystal form (space group C2, cell dimensions 94.8 A, 48.3 A, 131.0 A, 90.0 degrees, 100.3 degrees, 90.0 degrees), TIM is present in a ligand-free state. The asymmetric unit consists of two TIM subunits. Each of these subunits is part of a dimer which is sitting on a crystallographic twofold axis, such that the crystal packing is formed from two TIM dimers in two distinct environments. The two constituent monomers of a given dimer are, therefore, crystallographically equivalent. In the ligand-free state of TIM in this crystal form, the two types of dimer are very similar in structure, with the flexible loops in the "open" conformation. For one dimer (termed molecule-1), the flexible loop (loop-6) is involved in crystal contacts. Crystals of this type have been used in soaking experiments with 0.4 M ammonium sulphate (studied at 2.4 A resolution), and with 40 microM phosphoglycolohydroxamate (studied at 2.5 A resolution). It is found that transfer to 0.4 M ammonium sulphate (equal to 80 times the Ki of sulphate for TIM), gives rise to significant sulphate binding at the active site of one dimer (termed molecule-2), and less significant binding at the active site of the other. In neither dimer does sulphate induce a "closed" conformation. In a mother liquor containing 40 microM phosphoglycolohydroxamate (equal to 10 times the Ki of phosphoglycolohydroxamate for TIM), an inhibitor molecule binds at the active site of only that dimer of which the flexible loop is free from crystal contacts (molecule-2). In this dimer, it induces a closed conformation. These three structures are compared and discussed with respect to the mode of binding of ligand in the active site as well as with respect to the conformational changes resulting from ligand binding.

Structure of the complex between trypanosomal triosephosphate isomerase and N-hydroxy-4-phosphono-butanamide: binding at the active site despite an "open" flexible loop conformation

The structure of triosephosphate isomerase from Trypanosoma brucei complexed with the competitive inhibitor N-hydroxy-4-phosphono-butanamide was determined by X-ray crystallography to a resolution of 2.84 A. Full occupancy binding of the inhibitor is observed only at one of the active sites of the homodimeric enzyme where the flexible loop is locked in a completely open conformation by crystal contacts. There is evidence that the inhibitor also binds to the second active site of the enzyme, but with low occupancy. The hydroxamyl group of the inhibitor forms hydrogen bonds to the side chains of Asn 11, Lys 13, and His 95, whereas each of its three methylene units is involved in nonpolar interactions with the side chain of the flexible loop residue Ile 172. Interactions between the hydroxamyl and the catalytic base Glu 167 are absent. The binding of this phosphonate inhibitor exhibits three unusual features: (1) the flexible loop is open, in contrast with the binding mode observed in eight other complexes between triosephosphate isomerase and various phosphate and phosphonate compounds; (2) compared with these complexes the present structure reveals a 1.5-A shift of the anion-binding site; (3) this is the first phosphonate inhibitor that is not forced by the enzyme into an eclipsed conformation about the P-CH2 bond. The results are discussed with respect to an ongoing drug design project aimed at the selective inhibition of glycolytic enzymes of T. brucei.

Crystallographic binding studies with triosephosphate isomerases: conformational changes induced by substrate and substrate-analogues

TIM catalyses the interconversion of a triosephosphate aldehyde into a triosephosphate ketone. This is a simple chemical reaction in which only protons are transferred. The crystallographic studies of TIM from chicken, yeast and trypanosome complexed with substrate and substrate analogues are discussed. The substrate binds in a deep pocket. On substrate binding, large conformational changes are induced in three loops. As a result of these conformational changes in the liganded structure, the active site pocket is sealed off from bulk solvent and the sidechain of the catalytic glutamate becomes optimally positioned for catalysis.

In search of new lead compounds for trypanosomiasis drug design: a protein structure-based linked-fragment approach

A modular method for pursuing structure-based inhibitor design in the framework of a design cycle is presented. The approach entails four stages: (1) a design pathway is defined in the three-dimensional structure of a target protein; (2) this pathway is divided into subregions; (3) complementary building blocks, also called fragments, are designed in each subregion; complementarity is defined in terms of shape, hydrophobicity, hydrogen bond properties and electrostatics; and (4) fragments from different subregions are linked into potential lead compounds. Stages (3) and (4) are qualitatively guided by force-field calculations. In addition, the designed fragments serve as entries for retrieving existing compounds from chemical databases. This linked-fragment approach has been applied in the design of potentially selective inhibitors of triosephosphate isomerase from Trypanosoma brucei, the causative agent of sleeping sickness.