(Received for publication, February 6, 1996, and in revised form, October 22, 1996)
From the Lehrstuhl für Allgemeine Chemie und Biochemie,
TU-München, Vöttingerstrae 40, 85350 Freising, Germany
The kinetic and equilibrium isotope effects on the fructose-1,6-bisphosphate aldolase reaction have been determined using the rabbit muscle enzyme. The natural 13C abundance for both atoms participating in the bond splitting were measured in position C-1 of dihydroxyacetone phosphate and glyceraldehyde 3-P after irreversible conversion to glycerol-3-P and 3-phosphoglycerate, respectively, and chemical degradation. The carbon isotope effects were determined comparing the 13C content of the corresponding positions after partial and complete turnover, and after complete equilibration of the reactants. 13(Vmax/Km) on C-3 was 1.016 ± 0.007 and 0.997 ± 0.009 on position C-4, and the equilibrium isotope effects K12/K13 on these positions were 1.0036 ± 0.0002 and 1.0049 ± 0.0001.
The observed kinetic isotope effect on C-3 is discussed to originate from the formation of the enamine, which comes to equilibrium before the rate determining release of glyceraldehyde 3-P from the ternary complex. The equilibrium isotope effect is seen as the reason for an earlier-found relative 13C enrichment in position C-3 and C-4 of glucose and for varying enrichments in 13C of carbohydrates from different compartments of cells. The kinetic isotope effect is suggested to cause 13C discriminations in the C-3 pool in context with the hexose formation in competition with other dihydroxyacetone phosphate turnover reactions.
The relative enrichment of carbon-13 in the carboxyl group of amino acids, as observed by Abelson and Hoering (1) in 1961, was the first indication for the existence of non-statistical isotope distributions in biological compounds. Later results on acetic acid (2) and on acetoin (3) gave evidence that this observation was just one example of a common phenomenon. In order to find a general explanation for this observation Galimov (4) discussed that even in chemically unequilibrated systems, e.g. biological systems, a microscopic reversibility in enzymatic reactions is the origin for a thermodynamically ordered isotope distribution. The author's calculations could in fact explain some of the isotopic patterns of natural compounds known at that time. However, the presumption of a general thermodynamic equilibrium in biological systems is probably not realistic. In our opinion kinetic isotope effects on enzymatic reactions should be considered as primary causes for isotope discriminations. This has been proven for the primary CO2-fixing reactions (5-8). In secondary metabolism, the isotope effect on the pyruvate dehydrogenase reaction has been made responsible for the general depletion of 13C in metabolites of acetyl-CoA, such as fatty acids or isoprenoids (9-15).
More detailed interpretations were possible when the total isotopic
patterns of primary and secondary metabolites became available. Especially, our corresponding investigations on glucose indicated an
enrichment of 13C in positions C-3 and C-4, and a depletion
of 13C in positions C-1 and C-6 (Fig. 1) of
this important primary metabolite (16). We have also identified the
corresponding pattern in some direct metabolites of glucose (15, 17).
As one of the reasons for this glucose pattern, we have discussed
isotope effects on several reactions of the pentose phosphate cycle,
especially on the fructose-1,6-bisphosphate aldolase reaction. In order
to verify this assumption we have now determined the kinetic and the
equilibrium isotope effects on this reaction.
We used for our experiments the well characterized Class I rabbit
muscle aldolase (18, 19), which belongs to the same class as the
cytosolic plant enzyme (20, 21). Rabbit muscle aldolase is a
homotetramer with a molecular mass of 160 kDa, each subunit occurring
in an /
barrel (22). The enzyme catalyzes the reversible aldol
cleavage of 1-phospho-ketoses (mainly fructose 1,6-bisphosphate
(FBP))1 into dihydroxyacetone phosphate
(DHAP) and aldehydes (mainly glyceraldehyde phosphate (glyceraldehyde
3-P)). In the condensation reaction it is specific for dihydroxyacetone
phosphate with retention of pro-S-configuration of C-3 in
the ketoses. More than 50 different aldehydes are condensed with DHAP
(23, 24).
The currently understood enzyme mechanism consists of several distinct
steps (25-27), as shown in Fig. 2. They are preceded by
the opening of the predominant ring form of the substrate in solution;
step I is building up of the enzyme substrate complex, E·FBP; step II a Schiff base
(E·SBFBP) formation from this complex; in step
III the -C-C-bond cleavage to the "ternary" complex,
E·enamine·glyceraldehyde 3-P, the enamine of DHAP and
glyceraldehyde 3-P, takes place; step IV is the dissociation of
glyceraldehyde 3-P from the ternary complex; step V is the DHAP-Schiff
base formation, E·SBDHAP; step VI is a proton
transfer giving rise to the protonated Schiff base, E·DHAP; step VII is the dissociation of DHAP from its
Michaelis complex, E·DHAP. It is still not clear whether
the velocity of the total reaction rate is limited by the -C-C-bond
cleavage (step III) or by the release of the triose phosphates
(step IV and VII).
The main purpose of this work is to elucidate the role of the enzyme reaction for the non-statistical 13C pattern of glucose, however, in the context of this investigation we also expect to identify the rate-limiting step in the mechanism by the carbon isotope effect on the reaction. Yet, our main interest deals with the influence of the aldolase reaction on the non-statistical isotope distribution in carbohydrate metabolism.
All ordinary chemicals were of analytical grade and purchased from local suppliers. Fructose 1,6-bisphosphate-trisodium salt, C6H11O12P2Na3*8H2O, was from Boehringer Mannheim GmbH, Mannheim, Germany.
EnzymesTriose-phosphate isomerase (EC 5.3.1.1), glyceraldehyde-3-phosphate dehydrogenase (EC 1.2.1.12), and phosphoglycerate kinase (EC 2.7.2.3), used for enzymatic assays of fructose 1,6-bisphosphate, 3-phosphoglycerate, and glycerol 3-phosphate (28), were from local suppliers (Sigma, München, and Boehringer Mannheim GmbH). Fructose-1,6-bisphosphate aldolase (EC 4.2.1.13) from rabbit muscle, specific activity 0.15 µkat/mg, glycerol-3-phosphate dehydrogenase (EC 1.1.1.18), specific activity 2.83 µkat/mg, and glyceraldehyde-3-phosphate dehydrogenase (EC 1.2.1.12), specific activity 1.33 µkat/mg, used for the reaction assay, were purchased from Boehringer Mannheim GmbH. The activity of triose-phosphate isomerase in these enzymes, which is a common impurity of glycolytic enzymes, was below 0.1%.
Determination of the Kinetic Isotope Effect on the Aldolase Reaction and Isolation of the Reaction ProductsThe measurement of the kinetic isotope effect on the aldolase reaction was performed under pH- and temperature-controlled conditions (Dosimat, Methrom, Herisau, CH), with the competitive method after O'Leary (29). 1 mmol of fructose 1,6-bisphosphate, 0.1 mmol of NAD+, 25 µmol of phosphoglycolate, and 2 mmol of NaH2AsO4 were dissolved in 100 ml of water (high performance liquid chromatography grade). The pH was adjusted to 7.6 (25% NaOH). For "partial turnover" (approximately 10% or 1 h), an aliquot of 90 ml of this medium was incubated in the Dosimat vessel at 25 °C with 1.25 µkat of glycerol-3-phosphate dehydrogenase, 834 nkat of glyceraldehyde-3-phosphate dehydrogenase, and 33.3 nkat of fructose-1,6-bisphosphate aldolase. For the "total turnover" assay, to an aliquot of 10 ml of medium, 834 nkat of aldolase were added, while the activities of the other enzymes added were identical as in the partial turnover assay. After turnover, the reaction mixtures were heated for 10 min to 65 °C. The denatured proteins were eliminated by filtration (0.45 µm, Millipore, type HA), and the exact turnover was determined in an aliquot of the filtrate (enzymatic assays for fructose 1,6-bisphosphate, glycerol 3-phosphate, and 3-phosphoglycerate). The bulk of the filtrate was lyophilized overnight, and the residue was dissolved in 5 ml of water. The reaction products 3-phosphoglycerate and glycerol-3-phosphate were separated by anion-exchange chromatography of 1-ml aliquots (Parisil SAX, 10 µm 20 × 250 mm, Grom, Herrenberg, Germany) on a Sykam high performance liquid chromatograph system (Sykam, Gilching, Germany). Before addition of the medium, the column was equilibrated for 30 min with 5 mM KH2PO4 (pH 2.8) (flow rate 5 ml/min.). By a linear gradient between 5 and 150 mM KH2PO4 (pH 2.8) the reaction products were eluted. The fractions containing the reaction products, identified by their UV-absorption (200-350 nm) and enzymatic assays, were pooled, lyophilized, and stored at 4 °C. The column was regenerated with 25 ml of 500 mM KH2PO4 (pH 2.8).
Determination of the Equilibrium Isotope Effect on the Aldolase Reaction and Isolation of the Reaction Products270 mg of fructose 1,6-bisphosphate and 10 mg of 2-phosphoglycolate (for inhibition of the triose-phosphate isomerase) were dissolved in 10 ml of 20 mM sodium arsenate solution (pH 7.6, 25 °C). The solution was divided into 4 aliquots. In each, the reaction was started by the addition of 3.3 nkat of fructose-1,6-bisphosphate aldolase. The reaction was stopped after 8 h (2 aliquots), respectively, 24 h (2 aliquots), by heating (65 °C, 10 min), and the denatured protein was eliminated by filtration (0.45 µm, Millipore type HA).
To the solution, 0.1 µmol of NADH, 0.1 µmol of NAD+, 283 nkat of glycerol-3-phosphate dehydrogenase solution, and 133 nkat of glyceraldehyde-3-phosphate dehydrogenase solution were added, and by 1 h incubation the triose phosphates were converted into glycerol 3-phosphate and 3-phosphoglycerate, respectively. After this conversion the proteins were denatured by heating. The filtrate of the medium was concentrated to 1 ml, and the products were separated as described before, by anion-exchange chromatography. The pooled fractions (about 5 ml) were concentrated under reduced pressure to 500 µl and finally lyophilized directly in 700-µl autosampler vials, in which the chemical degradation was performed.
Degradation of Reaction Products for Positional Isotope AnalysisFor the periodate degradation of glycerol 3-phosphate
into formaldehyde and phosphoglycolaldehyde (30), the compound
(isolation and determination of total amount see above) was dissolved
in 500 mM KH2PO4 buffer (pH
5.8) to a final concentration of at least 100 mM, and a
120% molar excess sodium metaperiodate was added. The reaction was
complete after 5 min at room temperature. In the case of millimole samples, this was controlled by the determination of the formaldehyde formed (31), while for nanomole samples the turnover was determined simultaneously to the isotope ratio measurement in the GC-C-IRMS system. For the determination of the isotope ratio of the
phosphoglycolaldehyde formed, the solution was evaporated to dryness;
after addition of water, the evaporation was repeated twice. The
residue was directly submitted to wet combustion by adding 1 ml of 50%
H2SO4 and a two times molar excess of
KMnO4. The reaction was performed in an evacuated vessel at
room temperature, and was complete after 3 min. The CO2
produced was measured volumetrically and then condensed into a vial for
isotope ratio determination.
3-Phosphoglycerate was oxidized by Ce(IV)(NH4)4(SO4)4 to CO2 and phosphoglycolate (32). In a 10-ml vacutainer 3-phosphoglycerate (10-1000 nmol), isolated as described, was dissolved in 1 ml of 50 mM KH2PO4, and 2-2000 µl of reaction mixture (12.7 g of Ce(IV)(NH4)4(SO4)4 and 15 ml of concentrated H2SO4, diluted in 85 ml of water) were added. The vacutainer was evacuated and the mixture incubated for 1 h at 100 °C. With millimole samples the yield was controlled volumetrically, and the CO2 was transferred to the isotope ratio mass spectrometer for measurement. In the case of nanomole samples, after reaction the vacutainer was filled to a final gas volume of 1 ml with 10% H2SO4, and an aliquot of the gaseous phase was transferred by means of a gas tight syringe to the GC-C-IRMS system, where turnover and isotope ratio were determined.
Determination of Isotope Ratios and Isotope EffectsGas
samples in the millimole range and solid substances were analyzed in a
MM 903 isotope ratio mass spectrometer (VG Isogas, Middelwich,
Cheshire, Great Britan). Prior, solid samples were converted into
CO2 in an elemental analyzer Roboprep CN (Europa Scientific, Crewe, Great Britan). Gaseous samples in the nanomole range
were analyzed in an Isochrome I GC-C-IRMS system (VG-Isogas, Middelwich, Cheshire, Great Britan). The separation of CO2,
formic acid, and formaldehyde in the GC-C-IRMS system was performed on a PoraPLOT U column (0.32 mm × 5 m, Chrompack, Frankfurt,
Germany) using a temperature program (2 min 110 °C, 0.3 °C/min to
140 °C). All isotope ratios were corrected after Craig (33) and were expressed in 13C values2
against PDB. The calculation of isotope effects was done after O'Leary
(29).
The usual enzyme assay (28) for the determination of the enzyme activity of the fructose-1,6-bisphosphate aldolase is based on the isomerization of the trioses formed by triose-phosphate isomerase and the reduction of dihydroxyacetone phosphate using an excess of the subsidiary enzymes.
For the determination of the kinetic isotope effects it was important,
not only to convert the products as fast as possible, but also to do
this by independent reactions, in order to obtain two distinct
compounds as carriers of atoms 3 and 4 of the substrate. Therefore the
triose-phosphate isomerase, if present at all, was blocked by
2-phosphoglycolate (34), and the glyceraldehydephosphate formed was
oxidized by glyceraldehyde phosphate dehydrogenase to
3-phosphoglycerate, while dihydroxyacetone phosphate was independently reduced to glycerol 3-phosphate (Fig. 3). The coupled
reactions had the advantage of recycling NAD+, and thus the
coenzyme concentration could be kept constant. In addition, the
irreversibility of the reaction was attained by incubation in 20 mM arsenate (35).
Control for Isotope Effects on the Developed Method Itself
The evaluation of the separation and degradation procedure
for isotope effects was tested with reference mixtures corresponding to
the reaction medium but without enzymes. These mixtures were separated
as described under "Experimental Procedures," and the products were
degraded as indicated. As an isotope effect could also be implied in
the degradation of the products, 3-phosphoglycerate and glycerol
3-phosphate were degraded in millimole quantities under the conditions
applied for the analysis of the incubation products. The turnover of
the Ce(IV)-fission of 3-phosphoglycerate was controlled gas
volumetrically on the CO2 formed, in the case of the
glycerol 3-phosphate degradation (NaIO4) by colorimetric determination of the formaldehyde formed (32). The 13C
values of CO2, formaldehyde, and phosphoglycolaldehyde were measured as described under "Experimental Procedures." For both reaction products, the
13C values as calculated from
that of the fragments were compared to that of the original compound
(Table I). For position C-1 in both products they did
not differ by more than 0.5
. Hence an isotope effect on the
separation and degradation procedures would be negligible.
|
In addition, an isotope balance of the whole reaction was established
by comparison of the 13C value of fructose
1,6-bisphosphate to that of the sum of all fragments isolated (Table
II). The difference in the
13C value of
1.3
is small enough to exclude an isotope effect in the whole
reaction, even though the errors on each position should not be
neglected. We guess that the reason for these relatively high standard
deviations of the positional
13C value is the
application of an on-line method, permitting the use of the small
amount of only 50-100 nmol of substrate for the determination of the
isotope effects. On the other hand the advantage of this method was the
possibility to supply small incubation volumes for the experiments
(1-5 ml), to use small analytical columns for the separation, and
hence to have the possibility of making more independent experiments
with a limited amount of enzyme.
|
This isotope balance and the value assignment to the fragments also
reveals the specific isotope distribution in the substrate fructose
1,6-bisphosphate used. In contrast to natural glucose (16) the
positions C-3 and C-4 are depleted in 13C. As the origin of
the compound is not known, a discussion of the reasons is not possible.
In the present context it has to be pointed out that the measurement of
the isotope effect is always based on a relative shift in a given
position, and absolute values of the C-atom in question are not
important.
The kinetic isotope effect on the
fructose-1,6-bisphosphate aldolase reaction was determined according to
the competitive method (29). In the present case the 13C
value of the products in positions corresponding to positions C-3 and
C-4 of fructose 1,6-bisphosphate after partial and complete turnover
were determined after,
![]() |
(Eq. 1) |
|
The thermodynamic equilibrium of the fructose-1,6-bisphosphate aldolase reaction should be expressed as,
![]() |
(Eq. 2) |
![]() |
The calculation of the isotope distribution and hence the equilibrium
isotope effect was done by comparing the isotope abundance in the
triose phosphates measured with the isotope content of the
corresponding positions of the substrate fructose 1,6-bisphosphate in
equilibrium (13CSE), derived from its
original pattern (
13CS0) and taking into
account the shift due to the corresponding abundance of the triose
phosphates. Then the isotope content of the positions in question in
the triose phosphates in equilibrium (
13CPE)
was correlated to that of fructose 1,6-bisphosphate in equilibrium (
13CSE) using the equilibrium concentrations
according to the equilibrium constant for [glyceraldehyde 3-P] *
[DHAP]/[FBP] = 8.1 * 10
5 M (36) and the
values as compiled in Table IV.
![]() |
(Eq. 3) |
|
At the beginning, it has to be pointed out that presently in the reaction sequence in question (Fig. 2), the final release of DHAP (step VII) is discussed to be the rate-limiting step of the whole sequence (37, 38), however, on the other side in the performed Vmax/Km,FBP experiments, only the reactions before the first irreversible step, here the release of glyceraldehyde 3-P (step IV), can contribute to the observed isotope effects. Intriguingly we found on C-3 a kinetic isotope effect significantly different from unity, while the one on C-4 was practically unity. This is not compatible with the -C-C- bond fission itself as a reason for the isotope effect, but with the formation of the sp2 configuration between C-2 and C-3 leading to the enamine of DHAP. The measured value 1.016 is reasonable for such a conversion of a secondary alcohol to an sp2 terminal carbon, accompanied by the simultaneous fission of the adjacent C-C bond (39). This process also has to be understood as an equilibration between the secondary alcohol and the enamine structures, in which the equilibration is attained before the release of glyceraldehyde 3-P. The observed effect is hence an equilibrium isotope effect on C-3. On the other hand the formation of the carbonyl sp2 structure in glyceraldehyde 3-P does not affect the isotope abundance on C-4. In agreement, the observation of large secondary tritium isotope effects for the conversion of the substrate with 3H in positions C-3 and C-4 (1.33 and 1.37, respectively, Ref. 40) proves the existence of an equilibrium between intermediate sp2 structures before the release of glyceraldehyde 3-P, which is hence definitely the slowest partial step in Vmax/Km,FDB experiments.
We cannot explain at present the observed 13C equilibrium isotope effects on the aldolase reaction. The predictions of Cleland (39) for a substitution of a C-H bond for a C-C bond at C-3 and the substitution of a C-O bond for the cleaved C-C bond at C-4, as well as the suggestions of Hogg et al. (41) for the hydration of aldehydes predict an inverse isotope effect for C-4, hence a 13C enrichment in the released glyceraldehyde 3-P. On the other hand there are no examples available to which our observation could be compared.
Most important with regard to our main interest, the isotopic pattern
of hexoses and descendants (Fig. 4) is the existence of
the isotope effects in respect to the observed 13C
enrichments in positions C-3 and C-4. The thermodynamic isotope factors
(i values) after Galimov (4) do not generally correspond to the actually found 13C pattern of glucose (16). However,
the 13C enrichments in positions C-3 and C-4 (Fig. 4)
support the assumption that for the pattern in these positions, the
equilibrium isotope effects found in this investigation are mainly
responsible. The positive values,
K12/K13 = 1.004 for C-3
and 1.005 for C-4, respectively, indicate that the dissociation
equilibrium for the 12C isotopomers favors the trioses and
that consequently hexose isotopomers with 13C enrichments
in positions C-3 and C-4 are in equilibrium with these trioses,
explaining different enrichments of 13C in carbohydrates
from various compartments of plants (17). Even more, the found
equilibrium isotope effects also predict a relatively higher
13C enrichment in position C-4, and that is exactly what
has been found in the pattern (Fig. 1).
On the other hand, applying this to a realistic dynamic biological system like a plant we have to take into account that the existing fast transport, metabolism, and turnover of substrates and different conditions of metabolic fluxes prevent all chemical species from being in equilibrium. Triose phosphates are the most important intermediates of carbohydrate metabolism, providing both the carbon and the energy source for the metabolism of all other cell compartments. As we have found that in spite of the fast turnover and scrambling mentioned above, in most metabolites 13C enrichments dominate in the positions originating from C-3 and C-4 of glucose. Maintenance of an isotopic balance demands a corresponding depletion in representatives of minor pathways originating from trioses, and this may be caused by the kinetic isotope effect found.
The existence of a kinetic isotope effect on the aldolase reaction
would, for example, demand a depletion in 13C of position
C-1 of glycerol. As a matter of fact we have found that this compound
from various sources is depleted in 13C relative to glucose
from the same origin3 (Fig. 4). The
depletion is exclusively limited to position C-1, and a
"normalization" of the pattern in positions C-2 and C-3 as
reference shows that it can attain, relative to the expected pattern
from DHAP, 25 (e.g. glycerol in wine, less intense in glycerol from cattle fat). As glycerol is a minor product in the sources investigated, this depletion compensates the enrichments of a
few
in the main products. A quantitative calculation of a metabolic
and isotopic balance, as has been possible in the system
CO2/CH4 for a rumen simulation (42) will, in a
complicated system like a plant, only be possible after the
determination and isotopic analysis of all main compounds. This would
in principle permit the measurement of metabolic fluxes in plants and
cell compartments and will be investigated in the future.
![]() |
We thank Dr. I. A. Rose for the helpful comments on the manuscript.