(Received for publication, January 27, 1997, and in revised form, March 19, 1997)
From the Laboratoire d'Enzymologie et de Génie Génétique, Université de Nancy I, URA CNRS 457, B.P. 239, 54506 Vandoeuvre-les-Nancy Cédex, France
GapB-encoded protein of Escherichia
coli and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) share
more than 40% amino acid identity. Most of the amino acids involved in
the binding of cofactor and substrates to GAPDH are conserved in
GapB-encoded protein. This enzyme shows an efficient
non-phosphorylating erythrose-4-phosphate dehydrogenase activity (Zhao,
G., Pease, A. J., Bharani, N., and Winkler, M. E. (1995) J. Bacteriol. 177, 2804-2812) but a low phosphorylating
glyceraldehyde-3-phosphate dehydrogenase activity, whereas GAPDH shows
a high efficient phosphorylating glyceraldehyde-3-phosphate dehydrogenase activity and a low phosphorylating erythrose-4-phosphate dehydrogenase activity. To identify the structural factors responsible for these differences, comparative kinetic and binding studies have
been carried out on both GapB-encoded protein of Escherichia coli and GAPDH of Bacillus stearothermophilus. The
KD constant of GapB-encoded protein for NAD is
800-fold higher than that of GAPDH. The chemical mechanism of erythrose
4-phosphate oxidation by GapB-encoded protein is shown to proceed
through a two-step mechanism involving covalent intermediates with
Cys-149, with rates associated to the acylation and deacylation
processes of 280 s1 and 20 s
1,
respectively. No isotopic solvent effect is observed suggesting that
the rate-limiting step is not hydrolysis. The rate of oxidation of
glyceraldehyde 3-phosphate is 0.12 s
1 and is hydride
transfer limiting, at least 2000-fold less efficient compared with that
of erythrose 4-phosphate. Thus, it can be concluded that it is only the
structure of the substrates that prevails in forming a ternary complex
enzyme-NAD-thiohemiacetal productive (or not) for hydride transfer in
the acylation step. This conclusion is reinforced by the fact that the
rate of oxidation for erythrose 4-phosphate by GAPDH is 0.1 s
1 and is limited by the acylation step, whereas
glyceraldehyde 3-phosphate acylation is efficient and is not
rate-determining (
800 s
1). Substituting Asn for His-176
on GapB-encoded protein, a residue postulated to facilitate hydride
transfer as a base catalyst, decreases 40-fold the
kcat of glyceraldehyde 3-phosphate oxidation. This suggests that the non-efficient positioning of the C-1 atom of
glyceraldehyde 3-phosphate relative to the pyridinium of the cofactor
within the ternary complex is responsible for the low catalytic
efficiency. No phosphorylating activity on erythrose 4-phosphate with
GapB-encoded protein is observed although the Pi site is operative as
proven by the oxidative phosphorylation of glyceraldehyde 3-phosphate.
Thus the binding of inorganic phosphate to the Pi site likely is not
productive for attacking efficiently the thioacyl intermediate formed
with erythrose 4-phosphate, whereas a water molecule is an efficient
nucleophile for the hydrolysis of the thioacyl intermediate. Compared
with glyceraldehyde-3-phosphate dehydrogenase activity, this
corresponds to an activation of the deacylation step by
4.5
kcal·mol
1. Altogether these results suggest subtle
structural differences between the active sites of GAPDH and
GapB-encoded protein that could be revealed and/or modulated by the
structure of the substrate bound. This also indicates that a protein
engineering approach could be used to convert a phosphorylating
aldehyde dehydrogenase into an efficient non-phosphorylating one and
vice versa.
The glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH)1 is a tetrameric enzyme that catalyzes reversibly the oxidative phosphorylation of D-glyceraldehyde 3-phosphate (G3P) to form 1,3-diphosphoglycerate (1,3-dPG) in the presence of NAD and inorganic phosphate (1). The refined structures of several phosphorylating GAPDHs have already been reported (2-7). The currently accepted forward reaction pathway involves two steps, first the formation of a covalent ternary complex GAPDH-NAD-G3P preceding an oxidoreduction step that leads to a thioacyl enzyme intermediate and NADH, and second the phosphorylation that produces 1,3-dPG. The chemical mechanism of catalysis was extensively studied and is now well understood (8, 9). The structural determinants of the nicotinamide and adenosine subsites involved in coenzyme specificity that control the catalytic efficiency have been recently analyzed (10-13). The individual contribution of the amino acids implicated in the two phosphate binding sites named Ps and Pi sites has also been studied at the kinetic, structural, and energetic levels (14, 15).
GAPDH is a key enzyme of the glycolysis and gluconeogenesis pathways.
In the gram Escherichia coli, three distinct
gap genes have been isolated so far. The gapA gene was shown
to encode an efficient active GAPDH (16, 17). A second gap
gene, called gapB, was initially characterized by Alefounder
and Perham (18). This gene belongs to a cluster of genes coding for
several glycolytic enzymes. This gene organization is often found in
eubacteria and archaebacteria studied so far (Ref. 17 and references
cited therein) and is usually responsible for the GAPDH activity.
However, a peculiar situation is observed in E. coli. The
gapB gene, which is located within the glycolytic gene
cluster, does not seem to encode GAPDH activity (18). Indeed, Hillman
and Fraenkel (19) showed that a nonsense mutation in the
gapA gene abolishes GAPDH activity. Thus, the GAPDH activity
is only encoded by the gapA gene that is 22.6 min far from
the phosphoglycerate kinase gene on the E. coli chromosome.
This has led to the proposal of either an eukaryotic origin for the
gapA gene (18, 20, 21) or the presence of both
gapA and gapB genes in ancestor bacteria (22). A
third gene named gapC was recently
characterized.2 The presence of several
stop codons in the putative coding regions, clearly indicates that
gapC gene can not encode a functional GAPDH.
Amino acid sequence comparison of GAPDHs from E. coli and Bacillus stearothermophilus and GapB3-encoded protein shows 41.6 and 43.6% identity, respectively. Nearly all the amino acids essential for the chemical mechanism and for the binding of cofactor and substrates are conserved. In the present paper, the enzymatic properties of the protein encoded by the gapB gene are described and compared with those of the GAPDH from B. stearothermophilus. The results are discussed in relation to the three-dimensional structures of the E. coli (7) and B. stearothermophilus (3) GAPDHs and to the recent data of Zhao et al. (23), who showed that GapB-encoded protein displays an efficient non-phosphorylating erythrose-4-phosphate dehydrogenase activity.
The E. coli strain used
for wild-type and mutant GapB-encoded protein production was HB101
(F, supE44, hsd320 Crs
ms
, recA13, ara-14, proA2
lacY1, galK2, rpsL20, xyl-5, mti-1, straA) transformed with a plasmidic
construction containing the gapB gene under either its own
promoter or the gapA promoter.4
Site-directed mutagenesis was performed using the method of Kunkel et al. (24).
For purification, cells were harvested by centrifugation, resuspended
in buffer A (50 mM Tris-HCl, 2 mM EDTA, pH 8)
containing 3 mM -mercaptoethanol and sonicated. The
GapB-encoded protein was then isolated by a 45%
(NH4)2SO4 precipitation. The
contaminating proteins were removed by applying the enzymatic solution
onto gel filtration, ACA 34 resin at pH 8 (buffer A). The enzyme was then purified on a Q-Sepharose column equilibrated with buffer A,
followed by a linear gradient of KCl (0-0.2 M) using a
fast protein liquid chromatography system (Pharmacia, Uppsala). The GapB-encoded protein was eluted at 140 mM KCl. Purified
fractions were pooled, diluted in buffer A, and
(NH4)2SO4 was added to a final
concentration of 1 M. Final purification was achieved by applying the Q-Sepharose pool onto a phenyl-Sepharose column
equilibrated in buffer A containing
(NH4)2SO4 1 M, followed
by a linear gradient of (NH4)2SO4
(100-0%). The GapB-encoded protein was eluted at 200 mM
of (NH4)2SO4.
Phenyl-Sepharose-purified fractions were pooled, and
-mercaptoethanol was added to 1 mM final. Aliquots were
stored at
20 °C without extensive inactivation for several weeks.
Protein concentration during the purification procedure was estimated
by absorbance measurements at 280 nm using an extinction coefficient of
1.0 mg1·ml·cm
1. Purity of the
GapB-encoded protein was checked by electrophoresis on a 10%
SDS-polyacrylamide gel (25) followed by Coomassie Brilliant Blue R-250
staining. Its molecular weight was determined by mass spectrometry.
Enzyme concentration was determined spectrophotometrically, using an
extinction coefficient at 280 nm of 1.50 × 105
M
1·cm
1 deduced from the
method of Scopes (26).
Production and purification of wild-type GAPDH of B. stearothermophilus were performed as described previously by Corbier et al. (14).
Enzyme KineticsInitial rate measurements were carried out
at 25 °C on a Cary 2200 or a Kontron Uvikon spectrophotometer by
following the absorbance of NADH at 340 nm. For the forward reaction,
the experimental conditions were 40 mM triethanolamine, 2 mM EDTA, 50 mM K2HPO4, pH 8.9 (27), or 40 mM triethanolamine, 2 mM
EDTA, pH 8.9, in the absence of phosphate. For the reverse reaction
with 1,3-dPG, assays were performed in 25 mM imidazole,
EDTA 2 mM buffer, pH 7. The turnover number
(kcat) was calculated using a molar extinction coefficient at 280 nm of 1.50 × 105
M1·cm
1 for the GapB-encoded
protein and of 1.31 × 105
M
1·cm
1 for the GAPDH of
B. stearothermophilus. kcat is
expressed per site.
The initial rate data were fitted to the Michaelis-Menten relationship using least squares regression analysis to determine Vm and KM. All KM values were determined at saturating concentrations of the other substrates.
Purification of E4PCommercially available erythrose 4-phosphate (E4P) (Sigma) contained low levels of G3P, glucose 6-phosphate, and other contaminants (23). The true concentration of E4P was deduced from the amount of NADH formed under experimental conditions in which the E4P dehydrogenase reaction proceeded to completion. The average purity was estimated to be 60% as described by the manufacturer.
To isolate pure E4P, contaminant G3P was oxidized by GAPDH from B. stearothermophilus in the presence of phosphate and NAD, and a coupled system using lactate dehydrogenase was used for NADH recycling to displace the equilibrium of the reaction in favor of the formation of 1,3-dPG. The time of the reaction was chosen to oxidize totally G3P without significant oxidation of E4P. Enzymes were then removed by centrifugation after denaturation by addition of H2SO4. The solution was then diluted in imidazole buffer, pH 7, and applied onto a DEAE column (LKB5 PW) previously equilibrated in 25 mM imidazole buffer, pH 7. Elution was performed with a linear gradient of KCl (0-0.5 M) that allowed easy separation of E4P from 1,3-dPG. The content in E4P of the fractions was measured enzymatically using GapB-encoded protein and NAD in 40 mM triethanolamine buffer, 2 mM EDTA, pH 8.9.
Titration of the Cysteine Residues with 5,5Cysteine content was determined
using DTNB under nondenaturing and denaturing conditions. Progress
curves of thionitrobenzoate production for wild-type and mutant enzymes
were recorded at 412 nm in 50 mM Tris-HCl buffer, 2 mM EDTA, pH 8.2, for nondenaturing conditions. For
denaturing conditions, SDS was added (10% final) and enzyme was heated
10 mn at 70 °C. Enzyme concentration was 14.7 µN (in
sites) and DTNB concentration was 0.3 mM final. The amount
of thionitrobenzoate formed was calculated using an extinction coefficient at 412 nm of 13,600 M1·cm
1.
Iodoacetamide was used as a second order labeling probe of
cysteine 149 by measuring the protection against inactivation afforded by the coenzyme. Inactivation of GapB-encoded protein (5 µN) and GAPDH of B. stearothermophilus (0.2 µN) was studied at 37 °C in 0.1 M TES
buffer, EDTA 0.2 mM, pH 7.3. The concentration of
iodoacetamide was 100 µM. The concentration of NAD was
varied over the range 0-7.5 mM for GapB-encoded protein
and 0-5 µM for GAPDH. The concentration of NADH was
varied over the range 0-750 µM for GapB-encoded protein and 0-5 µM for GAPDH. Aliquots were withdrawn from the
incubation mixture at fixed intervals, and the residual activity was
then determined by dilution into the assay mixture. Assumption was made
that the chemical reactivity of cysteine 149 is similar in the apo- and
holoforms. The pseudo-first order constant kobs
was determined for each NAD or NADH concentration from plots of
log(A0/At) versus
time (A0 and At correspond to
the initial enzymatic activity and the activity at time t,
respectively). Dissociation constant (KD) of NAD or
NADH was determined from plots of (1 kA/k0)/[NAD] versus
kA/k0 (k0
and kA correspond to the pseudo-first order constant
kobs in the absence and presence of NAD or NADH,
respectively).
Fast kinetic measurements
were carried out on a Biologic Instruments (SFM3) stopped-flow
apparatus. The dead time of the apparatus under the conditions of flow
rate and measurement cell employed was determined using reduction of
2,6-dichloroindophenol by ascorbic acid (28). It was estimated at 1.4 ms while the apparent first order rate constant of the reaction was
proportional to the ascorbic acid concentration until about 900 s1, thus setting the upper limit of reliable
measurements. Data collected from absorbance measurements at 25 °C
were analyzed with the Biokine program using non-linear regression
analysis. An average of at least six runs was performed to determine
each constant.
Analysis of the steps leading to thioacyl enzyme formation was performed as described previously (14) except that G3P was replaced by E4P at final concentration of 2 mM.
Analysis of the alkylation of Cys-149 by
5,5-dithiobis-(2-nitrobenzoate) was made under nondenaturing
conditions. Progress curves of thionitrobenzoate production for
wild-type and mutant enzymes were recorded at 412 nm in 50 mM Tris-HCl buffer containing 2 mM EDTA, pH
8.2. One syringe was filled with enzyme (14.7 µN), and
the other one contained DTNB (600 µM).
Oxidation of G3P by GapB-encoded protein was done in the absence and presence of phosphate (100 mM). Oxidized products were isolated and enzymatically characterized as already described (9). The chemical nature of the reaction products was also confirmed by 31P NMR. Purified products were lyophilized and dissolved in D2O, and the pH was adjusted to 7. NMR spectral data were obtained on a Bruker AC250 spectrometer. Chemical shifts were referenced to external phosphoric acid at 85%. The signals of the C-3 phosphate were 4.47, 3.37, and 3.94 ppm for G3P, 3-phosphoglycerate, and 1,3-dPG respectively, and that of the C-1 phosphate of 1,3-dPG was 2.21 ppm. The chemical shift of free inorganic phosphate was 2.60 ppm.
Oxidation of E4P by GapB-encoded protein and GAPDH was also done in the absence and presence of phosphate (100 mM). The oxidized products were isolated using a similar procedure as for 3-phosphoglycerate and 1,3-dPG. 1,4-Diphosphoerythronate was detected using the reverse reaction catalyzed by GAPDH. 4-Phosphoerythronate was characterized by 31P NMR spectroscopy (chemical shift of 4-phosphoerythronate, 4.13 ppm).
Deuterium Isotope EffectsD-G3P and D-[1-2H]G3P were prepared as described previously (11). Their concentrations were enzymatically determined. Isotopic effects were measured in 40 mM triethanolamine buffer, 2 mM EDTA, pH 8.9, by direct comparison of initial velocities observed for oxidative phosphorylation of D-G3P and D-[1-2H]G3P (1 mM) in the presence of 10 mM NAD and 50 mM phosphate.
Justification of the Mutations
Site-directed mutagenesis at positions 32, 149, 176, 179, 206-209, and 311 were done for the following reasons (see also Fig. 1).
Positions 149 and 176
These positions are always occupied by a Cys and a His residue, respectively, in all GAPDHs described so far (see Fig. 1). Cys-149 forms a thioacyl intermediate during the oxidoreduction process, and His-176 is postulated to decrease the pKa of Cys-149 to favor the hydride transfer and to stabilize the different intermediates formed along the catalytic process (Ref. 8 and references cited therein). These positions were mutated to check if they could play a similar role during the E4P oxidation process.
Position 311This position is always occupied by a Tyr residue in all GAPDHs described so far except for the GapB-encoded protein where a Cys residue is present. Tyr-311 is at the hinge between catalytic and cofactor domains and is located near the catalytic amino acids Cys-149 and His-176 and also near Cys-153 (see Fig. 7 in Ref. 7). Modeling cannot exclude the possibility for Cys-311 to form an alternative thioacyl enzyme intermediate with E4P that has an additional CHOH group compared with G3P.
Position 179This position belongs to the Ps site and is always occupied by a Thr in all active GAPDHs described so far (15) except for the GapB-encoded protein where a Met residue is present (18).
Position 32This position is always occupied by an Asp
residue, which forms a hydrogen bond with both 2- and 3
- hydroxyl
groups of the ribose of the adenosine moiety, except for the
GapB-encoded protein where a Glu residue is present. Since the
substitution Asp-32
Glu in B. stearothermophilus GAPDH
increased 9-fold the KM of NAD (12), the presence of
a Glu-32 could be one of the factors responsible for the low affinity
of NAD(H) of the GapB-encoded protein (see "Results").
This sequence contained positions always invariant in GAPDH that are involved in the Pi site, i.e. Thr-208 and Gly-209. This sequence is largely changed in the GapB-encoded protein, where in particular Gly-209 is substituted by a Lys residue (see Fig. 1).
Biochemical Properties of Wild-type and Mutated GapB-encoded Protein
GapB-encoded protein was overexpressed in E. coli strain using a plasmidic construction containing the gapB gene under either its own promoter or the gapA promoter.4 Over 15% of the soluble proteins in the supernatant were GapB-encoded proteins. The protocol used to purify GapB-encoded protein to homogeneity took advantage of the higher hydrophobic character of GapB-encoded protein compared with the GAPDH from E. coli. This allowed easy separation of GapB-encoded protein from GAPDH on a hydrophobic column. This protocol is different from that used by Zhao et al. (23) (see Table I). The molecular weight of 37,170 ± 2 determined by mass spectrometry was in agreement with the 37,169 mass predicted from the gapB DNA sequence but is different to that calculated by Zhao et al. (23). The fact that GapB-encoded protein precipitated at low ammonium sulfate concentration (45% instead of 66% for all GAPDHs) and migrated faster in SDS-polyacrylamide gel than expected from its molecular weight compared with the GAPDH from E. coli and B. stearothermophilus (molecular weight of 35,401 and 35,944, respectively, gel not shown) is certainly related to its hydrophobic character. The protein was isolated as apo form, as judged by the ratio A280/A260 of 2. The presence of a reducing agent along the protocol of purification and during the storage of the purified enzyme had no effect on the activity in contrast to what Zhao et al. (23) have observed.
|
Thiol titrations by DTNB showed four titratable cysteines per monomer under denaturing conditions. This result is expected from the gapB DNA sequence that indicates cysteine residues at positions 95, 149, 153, and 311. Two cysteine residues showing different reactivity were titratable per monomer under native conditions, with one having a high reactivity compared with that of the second one. This suggested that, in addition to the Cys-149 that was shown to be accessible and highly reactive in GAPDH (8), a second cysteine residue was also titratable. The fact that two cysteine residues were titratable for C95S and C153S mutants and that only one cysteine was titratable for C311A, S and Y mutants indicates that Cys-311 is the second titratable cysteine.
By combining stopped-flow and classical spectrophotometry experiments,
kobs values of 136 s1 and
8.10
3 s
1 for wild type and of 116 s
1 for C311A mutant were determined in the absence of
NAD. This proved that Cys-149 is the most reactive titratable cysteine
within the active site. A kobs value of 242 s
1 was found for the titratable Cys-149 of the GAPDH of
B. stearothermophilus. Thus, the Cys-149 of GapB-encoded
protein has a chemical reactivity in the range of that observed for
GAPDH. It is noteworthy that in C149A, G, V mutants, the Cys-311 was
not titratable.
Kinetics
With G3P as a SubstrateThe kinetic parameters
(kcat and KM) of GapB-encoded
protein and its various mutants are summarized in Table
II. A 600-fold decrease in the
kcat value in the forward reaction was observed
for the GapB-encoded protein compared with the GAPDH from B. stearothermophilus. KM value for NAD increased 10-fold. To determine the nature of the limiting step and to define whether the catalytic mechanism proceeded via two distinct chemical steps as shown for GAPDH, transient kinetics in the absence or presence
of inorganic phosphate (50 mM) were carried out at
saturating NAD concentrations. Under these experimental conditions no
burst of NADH production was observed. The rate constant of 0.12 s1 is similar to the kcat value
obtained under steady state conditions regardless of the presence or
absence of inorganic phosphate (curves not shown). This demonstrated
that the limiting step is associated with the formation of the acyl
enzyme intermediate and not with steps occurring after the acylation,
as shown for the GAPDHs from E. coli and B. stearothermophilus (11, 14). The foregoing results indicated an
acylation step at least 7000-fold less efficient compared with the
GAPDH from B. stearothermophilus (15). The presence of an
isotopic effect of 5 with D-[1-2H]G3P as a
substrate (data not shown) demonstrated that the rate-limiting step is
associated with the hydride transfer.
|
The fact that acylation was limiting did not exclude an efficient
phosphorylating step. To investigate this possibility, the chemical
nature of the product formed in the forward direction in the absence
and presence of 0.1 M phosphate was determined enzymatically and characterized chemically by 31P NMR. It
was shown to be 3-phosphoglycerate and 1,3-dPG, respectively (see
"Materials and Methods"). Thus, the GapB-encoded protein exhibits a
phosphorylating activity more efficient than the non-phosphorylating one. Clearly, the Pi site is operational. In the reverse direction, 1,3-dPG was shown to be a substrate with a KM value
increased 16-fold, a KM value of NADH increased
5-fold, and a kcat decreased 880-fold compared
with GAPDH (KM (1,3-dPG) 0.08 and 0.005 mM, KM (NADH) 0.05 and 0.011 mM and kcat 0.074 s1
and 65 s
1 for GapB-encoded protein and GAPDH of B. stearothermophilus, respectively).
As shown in Table II, replacing Met-179 by Thr and Cys-311 by Tyr or Ala did not drastically change the GAPDH activity of GapB-encoded protein, whereas changing His-176 into Asn decreased kcat by a factor of 40. No significant activity was observed for C149G, A or V mutants under the experimental conditions used (data not shown, see "Discussion").
With E4P as a SubstrateThe data confirmed those already
published by Zhao et al. (23) showing an E4P dehydrogenase
activity. A kcat of 20 s1 and a
KM of 800 µM for NAD were found under
optimal conditions. The kcat value is 2.5-fold
smaller than that described previously (23). This discrepancy remains
to be explained. An E4P dehydrogenase activity for C149G, A and V
mutants, purified from a E. coli strain, was observed at
least 2000-fold less compared with wild type. This low activity is in
fact due to the GapB-encoded protein expressed from the chromosomic
gapB gene of E. coli strain used to overexpress
the mutants. Indeed the mutants purified from a
gapB
deleted strain showed no significant activity, at least 106-fold less compared with wild type (data not shown). The
fact that CD spectra of the mutants and wild type perfectly superposed (curves not shown) suggests that the absence of activity of the mutants
is not a consequence of change in their structure. Replacing amino
acids located in the catalytic subsite, i.e. Met-179 by Thr,
Cys-311 by Ala or Tyr, His-176 by Asn reduced
kcat by a factor of 5, 6, 33, and 95, respectively.
Pre-steady state stopped-flow experiments showed a burst of 1 mol of
NADH production per monomer with a kobs of 230 and 280 s1 at pH 8.2 and 8.9, respectively (curves not
shown). This supported a two-step mechanism with formation of a
thioacyl intermediate involving Cys-149 and a rate-limiting step
occurring after NADH formation. This limiting step could either be the
deacylation step or the release of 4-phosphoerythronate or NADH that
could include an isomerization step. The fact that no D2O
isotopic effect (data not shown) was observed rather excluded the first
hypothesis. The kinetics are similar in the presence of up to 100 mM phosphate. The fact that the deacylation step is not
rate-limiting does not exclude that phosphorylation can occur with an
efficiency similar or even higher than deacylation. Two kinds of
results argue, however, against this possibility. First, erythrose
4-phosphate oxidation proceeded to completion in the presence or
absence of 50 mM phosphate, as measured from NADH
production (data not shown), whereas, based on the equilibrium constant
of the oxidative phosphorylation of G3P (Keq = 2.10
8 (29)), the oxidative phosphorylation of E4P was
expected to yield only 56% 1,4-diphosphoerythronate under the
experimental conditions used (0.1 mM E4P, 0.1 mM NAD, 40 mM triethanolamine buffer, 2 mM EDTA, 50 mM phosphate, pH 8.9). Second, only
4-phosphoerythronate is formed as characterized by 31P NMR
(see "Materials and Methods").
The catalytic efficiency of the GapB-encoded protein tested with
erythrose is highly decreased (kcat 0.36 s1, KM (erythrose) 90 mM).
This points out the role of the phosphate at the C-4 position for
revealing the erythrose dehydrogenase activity.
E4P is also a substrate for GAPDH from B. stearothermophilus. kcat is reduced 760- and 200-fold when compared with that of G3P activity of GAPDH and E4P activity of GapB-encoded protein, respectively (see Table II). This activity is of phosphorylating type as judged by the yield of NADH production. 1,4-Diphosphoerythronate was isolated using a method similar to that used for 1,3-dPG (9) and characterized enzymatically using the reverse reaction of GAPDH. Transient kinetics in the absence or presence of inorganic phosphate (50 mM) showed no burst of NADH production, with a kobs similar to the kcat value (curves not shown). This showed that the limiting step is associated to the acyl enzyme formation.
NAD Binding to GapB-encoded Protein
To determine whether the higher KM value for NAD is indicative of a decrease of its affinity, KD value was determined using iodoacetamide as a second order labeling probe of cysteine by measuring the protection against inactivation afforded by the coenzyme binding. The deduced KD value provides a measure of macroscopic binding of NAD. Surprisingly, the NAD and NADH dissociation constants KD increased 800- and 300-fold, respectively, compared with GAPDH (see Table III). This can also be related to the observation that GapB-encoded protein is always isolated as an apo form. The low NAD affinity prevented the use of the Racker band (30) as a probe of the positioning of the nicotinamidium ring within the active site (10, 12).
|
Substituting Asp for Glu-32 did not change drastically kcat and KM values, but allosteric property with a Hill number of 2.3 was revealed with both G3P and E4P substrates (curves not shown) confirming that position 32 is important in revealing cofactor cooperativity in GAPDH (12).
E. coli strains deleted for gapA gene have been shown
to be unable to grow on glucose.5 This
indicated that the high GAPDH activity characterized in E. coli is only encoded by gapA gene and that GapB-encoded
protein, if it is translated from its gene, should provide a very low
GAPDH activity. The fact that the expressed gapB
gene-encoded protein could not suppress the glucose
phenotype of a
gapA strain supported this
conclusion.5 The kcat value of 0.12 s
1 with G3P as a substrate also agreed with this
assumption. Indeed, among all the GAPDH mutants from B. stearothermophilus studied so far, only those having a
kcat higher than 1 s
1 were shown
to complement a
gapA E. coli strain (data not shown). The
activity of GapB-encoded protein toward G3P is of phosphorylating type
and proceeds through a two-step mechanism. First, in the presence of
phosphate, formation of 1,3-dPG was observed. Second, 1,3-dPG is a
substrate for the reverse reaction. No significant activity is observed
when Cys-149 is replaced by a neutral amino acid, whereas replacement
of Cys-311 does not abolish the catalytic efficiency of GapB-encoded
protein (see "Results"). This supports the role of Cys-149 in the
formation of the thioacyl intermediate, as already shown for GAPDH (9).
Analysis of the rate constant in the forward direction indicated that
the limiting step is associated with the formation of the thioacyl
intermediate. These steps include G3P binding to the binary complex
enzyme-NAD, nucleophilic attack of the Cys-149 on the aldehydic
function, hydride transfer, and any potential isomerization step of the
ternary complex preceding the hydride transfer. The fact (a)
the chemical reactivity of Cys-149 was similar to that of GAPDH and
(b) an isotopic effect of 5 was observed when
D-[1-2H] G3P was used indicated a limiting
step associated with the hydride transfer. Thus, the efficiency of
hydride transfer is, at least, a 7000-fold decrease when compared with
B. stearothermophilus GAPDH (14).
Although NAD showed affinity decreased by a factor of 800 when compared with B. stearothermophilus GAPDH, this does not prevent an efficient acylation step from occurring. Indeed the enzyme showed a low acylation efficiency with G3P, whereas an efficient acylation step with E4P was observed, at least 2000-fold faster. Therefore, the structure of the substrate and its ability to form a productive ternary complex enzyme-NAD-substrate rather prevail, thereby allowing an efficient hydride transfer. The fact that E4P is not an efficient substrate for GAPDH from B. stearothermophilus (Table II) supports this hypothesis. Different factors can affect the hydride transfer efficiency, in particular the positioning of the hydrogen atom at C-1 relative to the C-4 atom of the nicotinamidium ring and/or the positioning of His-176 relative to the OH group at C-1 (the imidazole ring is postulated to play a role as a base catalyst facilitating hydride transfer). Changing His-176 into Asn decreased 40-fold the kcat for G3P. Hence, it is tempting to speculate that a non-efficient positioning of the C-1 relative to the pyridinium ring is responsible for the kcat decrease. Comparison of the amino acid sequence alignment of the cofactor subsite, i.e. 1-148 and 311-333, indicated that nearly all the residues shown to be involved in the binding of the adenosine (12) and nicotinamide part (11) are conserved, except Asp-32 (see Fig. 1). However, replacing Glu-32 by Asp did not improve the affinity of NAD. Thus other structural factors are involved to explain the drastic increase of KD. Determination of the three-dimensional structure of GapB-encoded protein in the absence and presence of NAD will help characterize these factors.
Two phosphate anion binding sites, Ps and Pi, were identified in the
three-dimensional structure of several GAPDHs (2, 3, 5, 7). The Ps and
Pi sites were historically described as those binding the C-3 phosphate
of G3P and the inorganic phosphate, respectively (2). The functional
role of the two anion binding sites along the catalytic mechanism has
not been totally elucidated so far. Modeling studies on glycosomal
GAPDH from Leishmania mexicana support G3P binding in the Ps
site and an inorganic phosphate binding in the Pi site (5). Although
there is an agreement for positioning the inorganic phosphate in the Pi
site to attack efficiently the thioester, an alternative mechanism has
been proposed for the steps preceding the thioester formation. The C-3
phosphate of G3P would first interact with the Pi site and then flip to the Ps site after the oxidoreduction step had occurred. The structure of the ternary complex GAPDH-NAD-glycidol-3-phosphate from B. stearothermophilus (3) and recent kinetic data obtained from mutants of the Pi and Ps sites favor this mechanism (14, 15). The Pi
site is composed of the side chain of Ser-148, Thr-150, Thr-208, and of
the peptide NH group of Gly-209, and Thr-179 and Arg-231 contribute to
the formation of the Ps site. These amino acids, which are invariant in
all the GAPDHs described so far, are also present in GapB-encoded
protein except, Thr-179 of the Ps site, which is replaced by a Met
residue. This substitution could therefore be responsible for the
inefficiency of the G3P acylation process by GapB-encoded protein.
However, this is unlikely. Indeed, although the T179M mutant of GAPDH
from B. stearothermophilus exhibited a decrease of its
kcat by a factor of 400, only the efficiency of
the phosphorylating step (and not of the acylation step) was affected
(15). Reintroducing a Thr at position 179 did not improve the catalytic
efficiency of GapB-encoded protein, in agreement with the prediction.
The sequence around the invariant residue Thr-208 of the Pi site (which
is located in E. coli GAPDH (7) on a loop between a
-strand and an
-helix) is also largely changed i.e.
S206V, S207D, G209K, A210L. In particular, the presence of the side
chain of Lys-209 could modify the conformation of the polypeptide main
chain around position 208-209 and thus prevent the hydrogen bonding
interaction of the main chain of Lys-209 and of the side chain of
Thr-208 with inorganic phosphate. However, as for position 179, no
improvement of the kcat was observed with G3P as
a substrate when the sequence of GAPDH of E. coli was
reintroduced in GapB-encoded protein. In fact,
kcat is largely decreased by a factor of 100 (kcat 1.4.10
3 s
1 and
1.1.10
3 s
1 for positions 209-210 and
positions 206-207-209-210, respectively), acylation remained
rate-limiting. Clearly, although phosphorylation occurs (as proved by
the formation of 1,3-dPG) and is more efficient than deacylation, we
cannot conclude whether or not inorganic phosphate is bound less
efficiency to the Pi site of GapB-encoded protein than to the GAPDH
site.
How do the above data, which show that GapB-encoded protein is not an efficient phosphorylating GAPDH, fit in with the recent data describing an efficient E4P non-phosphorylating dehydrogenase activity of GapB-encoded protein (23)? First, this new activity raised the question of the nature of the chemical mechanism involved in the oxidation of E4P. Kinetics clearly showed a two-step mechanism with a limiting step occurring after the acylating step. Molecular modeling of the catalytic domain taking the structure of GAPDH from E. coli as a model (see Fig. 7 in Ref. 7) shows an important modification near the catalytic residues Cys-149 and His-176, i.e. Tyr-311, which is conserved in all efficient GAPDHs described so far, is replaced by a Cys residue. Although the presence of the Cys-311 favored the E4P dehydrogenase activity, it is not essential for revealing the activity. Only C149V, A and G mutants lost activity. Clearly, these results support a common chemical mechanism for the oxidation of G3P and E4P catalyzed by GapB-encoded protein with formation of a thiohemiacetal intermediate via Cys-149.
The only structural difference between G3P and E4P is the presence of a
supplementary CHOH group in E4P. How is it possible to relate this
structural difference to the efficiency of the acylation step that is
at least 2000-fold higher for E4P than for G3P? As already mentioned,
the binding of E4P to the enzyme-NAD complex is postulated to favor the
formation of a ternary complex efficient for the hydride transfer. In
the absence of a known three-dimensional structure, it is difficult to
define on a structural point of view why E4P and not G3P is capable to
favor the formation of an efficient thiohemiacetal ternary complex. The
absence of a phosphorylating activity with E4P is also intriguing,
while a phosphorylating activity is observed with GAPDH. The binding of
inorganic phosphate to the Pi site of the GapB-encoded protein, which
should occur as proved by the oxidative phosphorylation of G3P, is not
productive for attacking efficiently the thioacyl intermediate formed
with E4P. Moreover, kcat is not changed in the
presence of a high concentration of phosphate. Although deacylation is
not rate-limiting, this supports a binding of inorganic phosphate to
the Pi site that does not perturb drastically the positioning and the
nucleophilicity of the water molecule attacking the thioacyl intermediate. The rate of deacylation with the E4P substrate for the
GapB-encoded protein is higher than 20 s1, whereas the
rate of deacylation with the G3P substrate for GAPDH is about
10
2 s
1 (15). From an energetic point of
view, this corresponds to an activation of at least 4.5 kcal·mol
1 of the deacylation step. Several kinds of
explanations (mutually not exclusive) can be proposed. This energy of
activation could correspond to the energy required for the water
molecule to attack the thioacyl intermediate. This implies the presence
in the active site of the GapB-encoded protein of a base catalyst
within the catalytic domain to favor the nucleophilicity of the water
molecule. In that context, it would be informative to characterize this amino acid. Another possibility is a more favorable positioning of the
thioacyl intermediate toward the attacking water molecule and/or a
destabilization of the thioacyl bond within the active site as a
consequence of different modes of binding of both substrates. Whatever
this issue, it is remarkable to note that GAPDH and GapB-encoded protein, which shares more than 40% of amino acid identity, have evolved from a common ancestor to be efficient as either a
phosphorylating or a non-phosphorylating aldehyde dehydrogenase.
Moreover, these activities are modulated by the structure of the
substrate. This indicates that it could be possible, by a protein
engineering approach, to convert a phosphorylating aldehyde
dehydrogenase into a non-phosphorylating one and vice versa.
Determination of the three-dimensional structure of GapB-encoded
protein and of inactive ternary complexes with G3P and E4P will be very
instructive in that regard.
Thanks are due to the Service Commun de Biophysicochimie of the University Henri Poincaré Nancy I for giving us the possibility to realize molecular modeling and to the Service Commun de Résonnance Magnétique Nucléaire of the University Henri Poincaré Nancy I for the NMR analyses. We are very grateful to Drs. N. Potier and A. Van Dorsselaer for determining the molecular weights of wild-type and mutants of GapB-encoded protein and to Dr. M. Hebrant for his technical assistance with stopped-flow experiments. We also thank E. Habermacher and J. P. Decle for their efficient technical help. We thank Drs. B. Charpentier and C. Branlant for providing gapB plasmid constructions.