|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 282, Issue 7, 4859-4867, February 16, 2007
Characterization, Kinetics, and Crystal Structures of Fructose-1,6-bisphosphate Aldolase from the Human Parasite, Giardia lamblia*![]() ![]() ![]() ![]() 1![]() ![]() 2
From the
Received for publication, October 10, 2006 , and in revised form, November 28, 2006.
Class I and class II fructose-1,6-bisphosphate aldolases (FBPA), glycolytic pathway enzymes, exhibit no amino acid sequence homology and utilize two different catalytic mechanisms. The mammalian class I FBPA employs a Schiff base mechanism, whereas the human parasitic protozoan Giardia lamblia class II FBPA is a zinc-dependent enzyme. In this study, we have explored the potential exploitation of the Giardia FBPA as a drug target. First, synthesis of FBPA was demonstrated in Giardia trophozoites by using an antibody-based fluorescence assay. Second, inhibition of FBPA gene transcription in Giardia trophozoites suggested that the enzyme is necessary for the survival of the organism under optimal laboratory growth conditions. Third, two crystal structures of FBPA in complex with the transition state analog phosphoglycolohydroxamate (PGH) show that the enzyme is homodimeric and that its active site contains a zinc ion. In one crystal form, each subunit contains PGH, which is coordinated to the zinc ion through the hydroxamic acid hydroxyl and carbonyl oxygen atoms. The second crystal form contains PGH only in one subunit and the active site of the second subunit is unoccupied. Inspection of the two states of the enzyme revealed that it undergoes a conformational transition upon ligand binding. The enzyme cleaves D-fructose-1,6-bisphosphate but not D-tagatose-1,6-bisphosphate, which is a tight binding competitive inhibitor. The essential role of the active site residue Asp-83 in catalysis was demonstrated by amino acid replacement. Determinants of catalysis and substrate recognition, derived from comparison of the G. lamblia FBPA structure with Escherichia coli FBPA and with a closely related enzyme, E. coli tagatose-1,6-bisphosphate aldolase (TBPA), are described.
Giardia lamblia, a flagellated protozoan, is the most common disease-causing parasite in developed countries. G. lamblia is also responsible for the most frequently diagnosed infective disease in developing countries (for a recent review see Ref. 1). The manifestations of giardiasis disease range from severe diarrhea, weight loss, vomiting, and malnutrition to asymptomatic carriage. G. lamblia is an ancient eukaryote. Its cellular biology and biochemistry are not well defined, and as a result, there is much to learn about its survival and adaptation tactics. Currently, giardiasis is treated by using metronidazole or tinidazole. Although the cure rates using tinidazole and metronidazole drug therapy are presently in excess of 60%, it is clear that alternative drug treatments are needed. Both drugs produce undesirable side effects, and increasing resistance to drug regimes presents a potential problem (2-6). Recurrence is also a problem, with frequency reported to be as high as 90% (4). In pursuit of alternative drug therapies, we have examined the recently reported G. lamblia genome sequence (7) in order to uncover genes that encode optimal targets. This analysis has led to our interest in the G. lamblia fructose-1,6-bisphosphate aldolase (EC 4.1.2.13 [EC] ) (glFBPA),3 an enzyme that catalyzes the reversible cleavage of D-fructose 1,6-bisphosphate (FBP) to dihydroxyacetone phosphate (DHAP) and D-glyceraldehyde 3-phosphate (G3P) (Fig. 1A). This process is a key step in the classical Embden-Meyerhof-Parnas glycolytic pathway. Because G. lamblia lacks mitochondria as well as the components of oxidative phosphorylation, glucose degradation via glycolysis serves as its major source of ATP (8-10). We envisioned that selective inhibition of FBPA might disrupt the functioning of the glycolytic pathway and thereby hinder survival of G. lamblia within its human host. Unlike the mammalian FBPA, which is a class I aldolase, glFBPA belongs to the class II family (11). The former FBPA uses an active site lysine residue to activate the substrate for C(3)-C(4) bond cleavage via Schiff base formation (Fig. 2). In contrast, class II FBPAs use a Zn2+ cofactor to coordinate both the C(2) carbonyl and C(3) hydroxyl oxygens of the substrate, consequently stabilizing the enediolate intermediate formed by C(3)-C(4) bond cleavage (12, 13) (Fig. 2). Importantly, our genomic sequence analysis indicates that FBPA is the only Giardia class II aldolase and that the Giardia genome does not encode a class I aldolase.
The mechanistically dissimilar class I and II aldolases do not share significant amino acid sequence homology. The evolutionary significance of this finding has been discussed previously (12, 14). The divergence in class I and II aldolase structures and catalytic mechanisms suggests a novel strategy for the design of drugs that operate by selectively inhibiting the activity of FBPA in the G. lamblia parasite but not in mammalian FBPA.
To date, the most thoroughly characterized class II FBPA is the enzyme from Escherichia coli (ecFBPA). Crystal structures have been determined for apo-ecFBPA and for the enzyme complexed with a transition state analog, phosphoglycolohydroxamate (PGH) (Fig. 1B) (15, 16). The structures, together with kinetic properties of active site site-directed mutants, have provided significant insight into the class II FBPA catalytic mechanism (17-19). The class II FBPA adopts an ( In this study, we describe our studies of glFBPA, which is the first example of a eukaryotic class II aldolase to be characterized. The x-ray structure of the enzyme in complex with the Zn2+ cofactor and PGH is reported, as are the kinetic constants for glFBPA catalysis of FBP and TBP cleavage. Finally, the essential role of the active site residue, Asp-83, in catalysis was confirmed by amino acid replacement and kinetic analysis of the mutant enzyme.
Protein Cloning, Expression, and PurificationThe trophozoites of the G. lamblia isolate WB, clone 1267, were grown as described previously (21). The genomic DNA preparation was performed with the DNA Stat kit (Stratagene). The PCR product was introduced into the pET100/D-TOPO expression vector (Invitrogen). Recombinant plasmids were isolated from the E. coli TOP10 strain (pET100/FBPAn).
For protein production, the E. coli strain BL21(DE3) Star was transformed with the pET100/FBPAn recombinant plasmid. Cells were grown in LB media at 30 °C to A600 = 0.5, when isopropyl 1-thio- Preparation of D83A glFBPA MutantThe QuickChange mutagenesis kit (Stratagene) was used in combination with the pET100/FBPAn plasmid template and the E. coli strain BL21(DE3) Star cell line. The mutant sequence was confirmed by DNA sequencing. Homogeneous D83A glFBPA was prepared using the same procedure described above for the wild-type enzyme.
Steady-state Kinetic Constant DeterminationInitial velocities were measured at 25 °C using 1-ml reaction solutions containing glFBPA (0.016 µM), 200 µM NADH, 5 units of triose-phosphate isomerase, 2 units of glycerol-3-phosphate dehydrogenase, and varying concentrations of FBP (Km, 0.5-10) in 50 mM K+HEPES (pH 7.5). The absorbance of the reaction solution was monitored at 340 nm (
To determine the kinetic constants, the initial velocity data were fitted to Equation 1 with KinetAsystI (IntelliKinetics, State College, PA),
where [S] is the substrate concentration; V0 is the initial velocity; Vmax is the maximum velocity; and Km is the Michaelis-Menten constant for the substrate. The kcat value was calculated from Vmax and the enzyme concentration [E] (determined employing the Bradford method (22)) using the equation kcat = Vmax/[E].
The competitive inhibition constants (Ki) of PGH (synthesized as the cyclohexylammonium chloride salt using a previously published protocol (23)) and TBP were determined by measuring the initial velocity of the above reactions at 0, 50, 100, and 150 nM PGH ([I]), and 6 and 20 µM TBP, and fitting the data to Equation 2 with KinetAsystI.
pH Rate Profile DeterminationThe initial velocity data were measured as a function of the pH by using the following (50 mM) buffer at the indicated pH values: MES and potassium acetate (4.0-5.0), MES (5.5), BisTris (6.0-6.5), HEPES (7.0-8.0), TAPS (8.5), CHES (9.0-9.5), and CAPS (10.0-10.5). The kcat and kcat/Km values were fitted to Equations 3 and 4, respectively, using the computer program KinetAsyst,
where Y is kcat or kcat/Km; [H] is the hydrogen ion concentration; C is the pH-independent value of kcat or kcat/Km; Ka is the acid dissociation constant; and Kb is the base dissociation constant. Analytical Size Exclusion ChromatographyFor estimation of the oligomeric state of the protein, analytical size exclusion chromatography was performed using a Superdex-200 HR 10/30 column (Amersham Biosciences). Runs were performed at a flow rate of 0.4 ml/min with 50 mM Tris-HCl (pH 7.5) and 0.1 M NaCl serving as eluant. Polyclonal Antibody Production and Immunofluorescence AssayPolyclonal antibodies against purified glFBPA were produced in BALB/c mice (Spring Valley Laboratories, Inc.). Fluorescein isothiocyanate-conjugated, affinity-purified goat anti-mouse immunoglobulin was used as the secondary antibody, and the control consisted of omission of the primary antibody. The specimens were examined with an Axioplan fluorescence microscope (Carl Zeiss, Thornwood, NY). Functional Knock-out Experiments of the FBPA Gene in G. lambliaThe FBPA gene was inserted into the plasmid pdsRNA between two opposite tetracycline-inducible Ras-related nuclear GTP-binding protein (ran) promoters, and knock-out experiments were carried out as described previously (24). Crystallization and Data CollectionCrystals were grown at room temperature in hanging drops using the vapor diffusion method. The protein solution was mixed with an equal volume of mother liquor containing 18-23% polyethylene glycol monomethyl ether 2000, 0.1 M Tris-HCl (pH 8.8), and 0.2 M MgCl2. Hexagonal and orthorhombic crystals were obtained from protein solutions that were incubated on ice for 30 min with 20 mM PGH. The hexagonal crystals diffracted x-rays to a resolution of 2.3 Å. The orthorhombic crystals appeared within 4-8 weeks, and they diffracted x-rays to a resolution of 1.75 Å. For data collection, the crystals were transferred to solutions containing mother liquor and 20% glycerol and flash-cooled with liquid propane cooled in liquid nitrogen. Diffraction data were acquired at the Industrial Macromolecular Crystallography Association-Collaborative Access Team (IMCA-CAT) 17-ID/BM beamlines at the Advanced Photon Source (Argonne National Laboratory, Argonne, IL). For data acquisition, the IMCA-CAT beamline was equipped with a Q210 CCD detector (ADSC, San Diego) or MARCCD 165 detector (Marresearch, Norderstedt, Germany). Data processing was carried out using HKL (25). The asymmetric unit of either crystal form contains two protein molecules. The statistics of data collection are provided in Table 1.
Structure Determination and RefinementThe structure of glFBPA was determined by molecular replacement using the program CNS (26) and the 2.3 Å diffraction data from the hexagonal crystal form. The search model was built based on E. coli tagatose-1,6-bisphosphate aldolase (ecTBPA), which shares 38% amino acid sequence identity with glFBPA (PDB entry 1GVF (20)). Structure refinement was carried out using the CNS program (26). The orthorhombic crystal form structure was obtained by molecular replacement using a partially refined structure of the hexagonal crystal form as the search model. The models were inspected and modified on a graphics work station using the program O (27). Water molecules were added to the model based on the Fo - Fc difference Fourier electron density map (where Fo and Fc are the observed and calculated structure factors, respectively), using peaks with density 3 as the acceptance criteria. Refinement statistics are provided in Table 1. PROCHECK was used for analysis of geometry (28), QUANTA for molecular modeling and structural alignment (Molecular Simulations Inc.), and PyMOL for depiction of the structures (29). Zinc Ion DetectionThe presence of zinc in the crystals was confirmed by x-ray fluorescence at the zinc absorption edge energy measured at the IMCA-CAT beam line 17-ID. Data were measured using a Bicron fluorescent detector while increasing the monochromator energy in small steps from 9600 to 9900 eV. Modeling the Enzyme-Substrate ComplexThe FBP coordinates, taken from the crystal structure of the class I rabbit liver FBPA (PDB entry 1FDJ [PDB] ), were used to place FBP in the active site of glFBPA. The DHAP portion of FBP was placed manually in a position equivalent to that of PGH. The positioning of the G3P portion of FBP was based on information from site-directed mutagenesis experiments, as detailed under "Results and Discussion." The final model was optimized by energy minimization of the docked FBP and residues surrounding the G3P portion of the substrate using the program CNS.
Identification of glFBPA as Potential Drug TargetglFBPA was selected together with a number of other potential Giardia drug target candidates based on genome-scale sequence analysis. The criteria for selection were as follows: 1) the likely absolute requirement for Giardia survival, and 2) either the absence in the human genome or sufficient sequence divergence from the human counterpart (30, 31). The sequencing of the G. lamblia genome is nearly complete (7). At the time of target selection, the list of Saccharomyces cerevisiae essential genes included 878 genes (32), comprising the most comprehensive cell viability data available for a eukaryotic organism. BlastP sequence analysis (33) of this list against the genome sequence of G. lamblia yielded 332 nontransmembrane proteins larger than 100 amino acid residues and that exhibit sequence homology to yeast essential genes with E score <10-3. The gene set was compared with the human genome data base (www.ncbi.nlm.nih.gov/genome/guide/human). One of the promising candidates that arose from this analysis was the gene encoding class II FBPA. In this context, it is noteworthy that Plasmodium falciparum and Trypanosoma brucei aldolases have been proposed as anti-malarial and anti-trypanosomal drug targets even though they belong to the class I enzymes and thus are closely related to the human enzyme (34, 35). Expression and Characterization of glFBPATo confirm that the gene encodes a true FBPA, the recombinant protein was expressed in E. coli. The theoretical mass of glFBPA is 35,376.9 Da, which is consistent with the mobility of the protein on SDS-polyacrylamide gels. Gel filtration analysis indicated that the protein is predominantly a homodimer in solution. Homodimeric association was also observed for ecFBPA (16), in contrast to FBPAs from Mycobacterium tuberculosis (36) and Thermus aquaticus (37), and ecTBPA (20), which associate into homotetramers. The steady-state kinetic constants for the glFBPA-catalyzed cleavage of FBP at pH 7.5 and 25 °C in the absence of added Zn2+ are kcat = 3.55 ± 0.05 s-1, Km = 1.7 ± 0.1 µM, and kcat/Km = 2.1 x 106 s-1 M-1. These values are within the range determined for other FBPAs (38). Addition of Zn2+ to reaction solutions did not increase the catalytic rate. In fact, at a concentration in excess of 30 µM ZnCl2, inhibition of catalysis was observed. This result suggested that the Zn2+ cofactor binds so tightly to glFBPA that the protein isolated from the E. coli cells is completely Zn2+-occupied. Moreover, attempts to remove the Zn2+ cofactor from glFBPA by exhaustive dialysis against EDTA in Zn2+-free buffer and then to restore catalytic activity by addition of Zn2+ were not successful. The dialyzed enzyme retained 50% of its original catalytic activity. The inclusion of ZnCl2 in the activity assay did not restore the lost activity. The pH dependence of glFBPA catalysis was determined by measuring the kcat and kcat/Km values as a function of reaction pH. Both pH rate profiles (Fig. 3) are bell-shaped, indicating that loss of activity occurs at acidic and basic pH. Therefore, all further kinetic determinations were carried out at pH 7.5 where FBPA displays optimal activity. The pH profiles were fitted to define the apparent pKa values reported in Fig. 3 legend. Class II FBPAs are known to be inhibited by PGH (23), an inhibitor that mimics the enediolate intermediate in the catalyzed reaction (Fig. 2). PGH is also known to inhibit triosephosphate isomerase (23), which promotes the ensuing step in the glycolytic pathway. Consequently, PGH has little potential for drug targeting of the glFBPA within the human host. Nevertheless, PGH is an excellent probe for identification of the glFBPA active site residues via x-ray crystal structure determination of the glFBPA-Zn2+-PGH complex. The binding affinity for PGH to glFBPA is defined by Ki = 37 ± 2nM.
An analysis of glFBPA sequence homologs suggests that the enzyme is more closely related to the ecTBPA (TBP is a C(4) epimer of FBP) than it is to the ecFBPA (38 and 23% sequence identity, respectively). This finding prompted the evaluation of the enzymatic activity of glFBPA with TBP. The turnover rate for TBP is below the detection limit of the assay (<5.5 x 10-5 s-1). The binding affinity for TBP to glFBPA is reflected by the measured Ki = 2.4 ± 0.2 µM. In Vivo Enzyme Expression and Essentiality TestPolyclonal antibodies raised against purified glFBPA were used in conjunction with a fluorescein isothiocyanate-conjugated secondary antibody to detect the expression of FBPA within G. lamblia trophozoites (Fig. 4). No fluorescence was observed for the control experiment in which trophozoites were treated with the fluorescein isothiocyanate-conjugated secondary antibody in the absence of primary antibody. These results show that the trophozoites produce FBPA in a form that is recognized by antibodies raised against the native enzyme and that the protein is distributed throughout the organism. The requirement of the enzyme for the Giardia parasite viability was examined by carrying out an in vivo experiment. This involved RNA interference/antisense RNA gene silencing techniques (39), recently adapted to G. lamblia trophozoites (24). Within 10 days after cultivation, trophozoites transfected with the silencing plasmid yielded no viable organisms even prior to induction by tetracycline. In control experiments that employed the nonessential gene encoding adaptor protein 1 involved in protein trafficking (24), viable trophozoites were observed after the transfection. The ran promoter can drive a low level of expression of the reporter gene even in the absence of tetracycline (40). Presumably, the "leakage" of the ran promoter enabled double-stranded RNA synthesis without induction. This triggered transcriptional inhibition and glFBPA gene silencing. The results of the transcriptional inhibition experiments are consistent with the crucial role that the targeted glFBPA plays in the Giardia glycolytic cycle as was suggested by genome analysis. Overall Structure of glFBPAThe structure was determined by molecular replacement using the ecTBPA structure as a search model, which exhibits higher sequence identity to glFBPA than ecFBPA (Table 2). Indeed, the glFBPA structure more closely matches that of ecTBPA than ecFBPA despite the difference in substrate specificity of the former two enzymes (Table 2). This is opposite of what would be predicted based on substrate specificity alone.
Consistent with the results of the analytical gel filtration experiments, glFBPA associates into a tightly packed dimer with approximate noncrystallographic 2-fold symmetry (Fig. 5). Each subunit (denoted A and B) folds into an ( / )8 barrel as seen in structures of other bacterial class II aldolases (15, 16, 20, 37).
In the hexagonal crystal form, the two subunits are similar, as reflected by the low root-mean-square deviation (r.m.s.d.) of their
In the orthorhombic crystal form, the active site of one subunit (molecule A) is occupied by Zn2+ and PGH, whereas the active site of molecule B is not. Superposition of the two subunits yields an r.m.s.d. value between
The relatively large r.m.s.d. value between the ligand-bound and ligand-free molecules is the consequence of the fact that two loops, spanning residues 174-194 and residues 227-237, have undergone conformational transitions. The conformational changes do not appear to be the result of crystal packing. Most notably, the 174-194 loop is largely disordered in the unbound state, whereas, with the exception of residues 186-189, the loop conformation is defined in the bound state. Upon ligand binding, the two loops adopt ordered conformations, with several residues stationed to shape the active site cavity. Key residues on the 174-194 loop are involved in PGH and Zn2+ binding (Lys-182 and His-178) (Fig. 6). The 227-237 loop undergoes substantial rearrangement, with the Pro-231 C
glFBPA PGH-Zn2+-binding SitesZn2+ was not added to any of the solutions during protein purification and crystallization. Yet exposure of a crystal of the PGH-bound enzyme to synchrotron-tunable x-rays revealed a peak at 9674.2 eV corresponding to the absorption edge of zinc. Difference Fourier electron density maps show a bound ligand and Zn2+ in both molecules of the hexagonal crystal form and only in molecule A of the orthorhombic crystal form (Fig. 6A). The structures of the ligand-bound glFBPA are the same in both crystal forms. Because the orthorhombic crystal form was determined at higher resolution (1.75 versus 2.3 Å), the analysis and the following discussion refer to the more accurately refined structure.
The glFBPA active site is a highly polar cavity located at the C-terminal end of the
Zn2+ SiteThe bound Zn2+ has a trigonal bipyramidal coordination geometry similar to that seen in the crystal structures of ecTBPA and ecFBPA. Three Zn2+ ligands are provided by the N The DHAP SiteThe tight binding inhibitor PGH mimics the enediolate intermediate, which, in the ensuing step of catalysis, undergoes protonation to form DHAP (Fig. 2). The inhibitor binds in a deep polar cavity with the phosphate moiety stationed in a site enriched with the positive dipoles of four backbone amide groups (Gly-179, Ser-213, Asp-255, and Ser-256) (Fig. 6B). In addition, the phosphate group interacts with the amino group of Lys-182 on the 174-194 flexible loop and with the hydroxyl groups of Ser-213 and Ser-256. The phosphate-binding site is very similar to that present in ecFBPA and ecTBPA. An interesting difference between the bacterial aldolases and glFBPA is that the monovalent cation that counters the phosphate negative charge in the bacterial aldolases is absent in glFBPA. The cation is replaced by the ammonium substituent of the Lys-182 side chain. In addition to coordination with Zn2+, the hydroxamate moiety of the PGH interacts with the backbone amide group of Gly-211, Asn-253, and Asp-83 (Fig. 6B). Asn-253 and Asp-83 are conserved in ecFBPA (numbered Asn-286 and Asp-109). The G3P SiteThe interaction of groups in glFBPA with the G3P moiety was assessed by using a model of the enzyme-FBP complex (Fig. 7). In constructing this model, the PGH-binding mode was used to define the orientation of the DHAP-forming end of FBP and interactions with the G3P-forming end were optimized. Inspection of the model leads to the prediction that Asp-83 binds the FBP through electrostatic interactions formed between the side chain carboxylate group and the substrate C(3) and C(4) hydroxyl groups. In addition, it suggests that Asp-83 functions to deprotonate the C(4) hydroxyl during catalytic turnover. To test this prediction, the D83A mutant was prepared and characterized. Although this mutant displays native protein behavior during purification, it does not have significant catalytic activity when subjected to forcing conditions (viz. 170 µM D83A FBPA with 2 mM FBP at pH 7.5 and 25 °C for 3 h). This condition defines a maximum turnover rate that is 5 orders of magnitude slower than the kcat = 3.55 ± 0.05 s-1 of the wild-type enzyme. The corresponding ecFBPA mutant D109A is also known to be catalytically impaired (17).
The properties of site-directed mutants of ecFBPA active site residues support the position of the G3P-forming end of the FBP. The studies indicate that Arg-331, Ser-61, Asp-288, and Asn-35 of the E. coli enzyme interact with the C(6) phosphate group and the C(4) and C(5) hydroxyl groups of FBP (17, 18, 42). Inspection of the glFBPA-FBP model leads to the prediction that the FBP C(6) phosphate interacts with the Ser-50 hydroxyl, guanidinium group of Arg-259 and guanidinium group of Arg-280 (located on the opposing subunit). The observation that both subunits contribute to the active site interactions has been noted in earlier studies with ecFBPA, where the counterpart of Arg-280 is Arg-331 (42). In contrast, Arg-259 of glFBPA does not have a counterpart in the ecFBPA. In ecFBPA this position is occupied by a glutamine (Gln-292). Consequently, it appears that the glFBPA active site contains a more extensive charge network. Specifically, Arg-259 is anchored by the two carboxyl groups of Asp-255 and Asp-278 (the latter residing on the partner subunit), and in turn, Asp-278 interacts with Arg-280 (Fig. 7). Interaction of the G3P phosphate group with the hydroxyl group of Ser-50 is conserved in the class II FBP aldolases (Ser-61 in ecFBPA). This is an important interaction as indicated by characterization of the S61A ecFBPA mutant (18). In the glFBPA-FBP model, the amide group of Asn-24 interacts with the G3P phosphate group as well as with the C(5) hydroxyl group of FBP. This is consistent with the observation that replacement of Asn-35, the corresponding residue in ecFBPA, by alanine results in impaired enzyme activity (17). The glFBPA-FBP-Zn2+ active site model described above serves two purposes. First, it provides a starting point for inhibitor design. Second, it provides insight into why TBP is not recognized as a substrate. Aldolase II Substrate SpecificityThe results of this study demonstrate that that TBP, the FBP C(4) hydroxyl epimer, binds tightly at the glFBPA substrate site (Ki = 2.4 µM), but it is not converted to products. As shown by Hunter and co-workers (20) (and demonstrated in Figs. 7 and 8), it is evident that if TBP binds to the glFBPA active site in the same conformation as does FBP, the C(4) hydroxyl group would not be properly aligned for deprotonation by the Asp-83. Thus, from this perspective the basis for the stringent specificity is obvious. Because the TBPA aspartic acid that deprotonates the TBP C(4) hydroxyl is located in exactly the same position as Asp-83 in glFBPA (Fig. 8), rotation about the C(3)-C(4) bond in TBP must occur in order to bring the C(4) hydroxyl close to the aspartic acid. However, the atoms beyond C(4) are forced into different spatial orientations with respect to the surrounding active site residues in this rotamer. This results in bad steric interactions with groups of both TBPA and FBPA (recall that the active site architecture of the two enzymes is very similar). To alleviate these clashes, additional bond rotation must take place, which could lead to nonproductive binding and no turnover.
In E. coli FBP and TBP are converted to G3P and DHAP by two different class II aldolases as follows: ecFBPA and ecTBPA. The ecFBPA substrate discrimination factor, expressed as (kcat/Km(FBP))/(kcat/Km(TBP)), was determined to be 1423 (18). The ecTBPA substrate discrimination factor, expressed as (kcat/Km(TBP))/(kcat/Km(FBP)), is 333. Although the substrate specificity is not as high in the two E. coli aldolases as it is in the glFBPA, it is significant. Hunter and co-workers (20) identified Asp-288 and Gln-59 as residues that play an important role in discriminating between FBP and TBP substrates. In glFBPA, Asp-288 corresponds to Asp-255 and Gln-59 corresponds to Gln-48. In ecTBPA, the positions of these two residues are occupied by alanine (Fig. 8). The fact that glFBPA retains the "ecFBPA Asp-288" as Asp-255 is consistent with the substrate preference for FBP. Hunter and co-workers (20) suggested Gln-59 might influence the sidechain orientation of another conserved residue positioned to contact FBP, Asn-35 (Asn-24 in ecTBPA and Asn-24 in glFBPA), and that the smaller alanine side chain leads to an alternate conformation of Asn-24 in TBPA. This conclusion is inconsistent with the observation that the active site arrangement of glFBPA has Asn-24 oriented in a nearly identical manner as that observed for Asn-24 in ecTBPA. Moreover, site-directed mutagenesis experiments that attempted to "switch" specificity failed to provide a clear-cut answer (18). Clearly, there must be other factors involved in determining substrate specificity. Additional studies (e.g. structure of an enzyme-substrate complex with an ordered active site loop) are required to gain better insight into this issue. ConclusionWe have identified the key glycolytic cycle enzyme FBPA as a potential Giardia drug target. The glFBPA, which belongs to the class II aldolase family, contrasts with the human FBPA that belongs to the evolutionarily unrelated class I aldolase family. We have shown that Giardia trophozoites produce this enzyme in the cytoplasm and that RNA silencing results in death of the parasite. Although both the amino acid sequence and crystal structure of glFBPA are more similar to those of ecTBPA than ecFBPA (38 and 23% sequence identity, respectively), the kinetic characterization shows unequivocally that the true substrate is FBP and that its C(4) epimer, TBP, is a glFBPA inhibitor. This finding underscores the importance of determining function experimentally rather than inferring it solely based on sequence/structure comparison.
The atomic coordinates and structure factors (code 2ISV, 2ISW) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/).
* This work was supported by National Institutes of Health Grant RO1 AI059733. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
1 Present address: Chemistry Dept., Division of Natural Sciences, Mount Marty College, Yankton, SD 57078. 2 To whom correspondence should be addressed: Center for Advanced Research in Biotechnology, 9600 Gudelsky Dr., Rockville MD 20850. Tel.: 240-314-6245; Fax: 240-314-6255; E-mail: osnat{at}carb.nist.gov.
3 The abbreviations used are: glFBPA, Giardia FBPA; FBPA, fructose-1,6-bisphosphate aldolase; ecTBPA, E. coli tagatose-1,6-bisphosphate aldolase; ecFBPA, E. coli FBPA; PGH, phosphoglycolohydroxamate; FBP, D-fructose-1,6-bisphosphate; TBP, D-tagatose-1,6-bisphosphate; DHAP, dihydroxyacetone phosphate; G3P, D-glyceraldehyde 3-phosphate; MES, 2-(N-morpholino)ethanesulfonic acid; BisTris, 2,2-bis(hydroxymethyl)-2,2',2-nitrilotriethanol; TAPS, N-tris(hydroxymethyl)methyl-3-aminopropanesulfonic acid; CHES, N-cyclohexyl-2-aminoethanesulfonic acid; CAPS, N-cyclohexyl-3-aminopropanesulfonic acid; PDB, Protein Data Bank; r.m.s.d., root-mean-square deviation; IMCACAT, Industrial Macromolecular Crystallography Association-Collaborative Access Team.
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||