Quantitative characterization of homo- and heteroassociations of muscle phosphofructokinase with aldolase

Dissociation of purified phosphofructokinase accompanied with inactivation was analyzed in the absence and presence of aldolase and the data were compared with those obtained with muscle extract. The kinetics of the decrease in enzymatic activity was highly dependent on the dilution factor in both cases, but the inactivation appeared to be biphasic only with extract. The inactivation of the phosphofructokinase was impeded by addition of excess of aldolase. Time courses of kinase inactivation were fitted by alternative kinetic models to characterize the multiple equilibria of several homo- and hetero-oligomers of phosphofructokinase. The combination of modeling data obtained with purified and extract systems suggests that aldolase binds to an intermediate dimer of phosphofructokinase and within this heterocomplex the kinase is completely active. The intermediate dimer is stabilized by association with microtubules and the kinase activity decreased due to dilution can be recovered by addition of excess aldolase. In extract, the phosphofructokinase is of sigmoidal character (Hill coefficient of 2.3); the addition of excess exogenous aldolase to phosphofructokinase resulted in heterocomplex formation displaying Michaelian kinetics. The possible physiological relevance of heterocomplex formation of phosphofructokinase in muscle extract is discussed.

Advantages and disadvantages of aggregating fluxes into synthetic and degradative fluxes when modelling metabolic pathways

It is now widely accepted that mathematical models are needed to predict the behaviour of complex metabolic networks in the cell, in order to have a rational basis for planning metabolic engineering with biotechnological or therapeutical purposes. The great complexity of metabolic networks makes it crucial to simplify them for analysis, but without violating key principles of stoichiometry or thermodynamics. We show here, however, that models for branched complex systems are sometimes obtained that violate the stoichiometry of fluxes at branch points and as a result give unrealistic metabolite concentrations at the steady state. This problem is especially important when models are constructed with the S-system form of biochemical systems theory. However, the same violation of stoichiometry can occur in metabolic control analysis if control coefficients are assumed to be constant when trying to predict the effects of large changes. We derive the appropriate matrix equations to analyse this type of problem systematically and to assess its extent in any given model.

Oxythiamine and dehydroepiandrosterone induce a G1 phase cycle arrest in Ehrlich's tumor cells through inhibition of the pentose cycle

Transketolase (TK) reactions play a crucial role in tumor cell nucleic acid ribose synthesis utilizing glucose carbons, yet, current cancer treatments do not target this central pathway. Experimentally, a dramatic decrease in tumor cell proliferation after the administration of the TK inhibitor oxythiamine (OT) was observed in several in vitro and in vivo tumor models. Here, we demonstrate that pentose cycle (PC) inhibitors, OT and dehydroepiandrosterone (DHEA), efficiently regulate the cell cycle and tumor proliferation processes. Increasing doses of OT or DHEA were administered by daily intraperitoneal injections to Ehrlich's ascites tumor hosting mice for 4 days. The tumor cell number and their cycle phase distribution profile were determined by DNA flow histograms. Tumors showed a dose dependent increase in their G0-G1 cell populations after both OT and DHEA treatment and a simultaneous decrease in cells advancing to the S and G2-M cell cycle phases. This effect of PC inhibitors was significant, OT was more effective than DHEA, both drugs acted synergistically in combination and no signs of direct cell or host toxicity were observed. Direct inhibition of PC reactions causes a G1 cell cycle arrest similar to that of 2-deoxyglucose treatment. However, no interference with cell energy production and cell toxicity is observed. PC inhibitors, specifically ones targeting TK, introduce a new target site for the development of future cancer therapies to inhibit glucose utilizing pathways selectively for nucleic acid production.

New semisynthetic vinca alkaloids: chemical, biochemical and cellular studies

A new semisynthetic anti-tumour bis-indol compound, KAR-2 [3'-(beta-chloroethyl)-2',4'-dioxo-3,5'-spiro-oxazolidino-4-dea cetoxy-vinblastine] with lower toxicity than vinca alkaloids used in chemotherapy binds to calmodulin but, in contrast to vinblastine, does not exhibit anti-calmodulin activity. To investigate whether the modest chemical modification of bis-indol structure is responsible for the lack of anti-calmodulin potency and for the different pharmacological effects, new derivatives have been synthesized for comparative studies. The synthesis of the KAR derivatives are presented. The comparative studies showed that the spiro-oxazolidino ring and the substitution of a formyl group to a methyl one were responsible for the lack of anti-calmodulin activities. The new derivatives, similar to the mother compounds, inhibited the tubulin assembly in polymerization tests in vitro, however their inhibitory effect was highly dependent on the organization state of microtubules; bundled microtubules appeared to be resistant against the drugs. The maximal cytotoxic activities of KAR derivatives in in vivo mice hosting leukaemia P388 or Ehrlich ascites tumour cells appeared similar to that of vinblastine or vincristine, however significant prolongation of life span could be reached with KAR derivatives only after the administration of a single dose. These studies plus data obtained using a cultured human neuroblastoma cell line showed that KAR compounds displayed their cytotoxic activities at significantly higher concentrations than the mother compounds, although their antimicrotubular activities were similar in vitro. These data suggest that vinblastine/vincristine damage additional crucial cell functions, one of which could be related to calmodulin-mediated processes.

Mass isotopomer study of the nonoxidative pathways of the pentose cycle with [1,2-13C2]glucose

We present a single-tracer method for the study of the pentose phosphate pathway (PPP) using [1,2-13C2]glucose and mass isotopomer analysis. The metabolism of [1,2-13C2]glucose by the glucose-6-phosphate dehydrogenase, transketolase (TK), and transaldolase (TA) reactions results in unique pentose and lactate isotopomers with either one or two 13C substitutions. The distribution of these isotopomers was used to estimate parameters of the PPP using the model of Katz and Rognstad (J. Katz and R. Rognstad. Biochemistry 6: 2227-2247, 1967). Mass and position isotopomers of ribose, and lactate and palmitate (products from triose phosphate) from human hepatoma cells (Hep G2) incubated with 30% enriched [1,2-13C2]glucose were determined using gas chromatography-mass spectrometry. After 24-72 h incubation, 1.9% of lactate molecules in the medium contained one 13C substitution (m1) and 10% contained two 13C substitutions (m2). A similar m1-to-m2 ratio was found in palmitate as expected. Pentose cycle (PC) activity determined from incubation with [1,2-13C2]glucose was 5.73 +/- 0.52% of the glucose flux, which was identical to the value of PC (5.55 +/- 0.73%) determined by separate incubations with [1-13C] and [6-13C]glucose, 13C was found to be distributed in four ribose isotopomers ([1-13C]-, [5-13C]-, [1,2-13C2]-, and [4,5-13C2]ribose). The observed ribose isotopomer distribution was best matched with that provided from simulation by substituting 0.032 for TK and 0.85 for TA activity relative to glucose uptake into the model of Katz and Rognstad. The use of [1,2-13C2]glucose not only permits the determination of PC but also allows estimation of relative rates through the TK and TA reactions.

Thiamine supplementation to cancer patients: a double edged sword

The objectives of this review are to (a) explain the mechanism by which thiamine (vitamin B1) promotes nucleic acid ribose synthesis and tumor cell proliferation via the nonoxidative transketolase (TK) pathway; (b) estimate the thiamine intake of cancer patients and (c) provide background information and to develop guidelines for alternative treatments with antithiamine transketolase inhibitors in the clinical setting. Clinical and experimental data demonstrate increased thiamine utilization of human tumors and its interference with experimental chemotherapy. Analysis of RNA ribose indicates that glucose carbons contribute to over 90% of ribose synthesis in cultured cervix und pancreatic carcinoma cells and that ribose is synthesized primarily through the thiamine dependent TK pathway (> 70%). Antithiamine compounds significantly inhibit nucleic acid synthesis and tumor cell proliferation in vitro and in vivo in several tumor models. The medical literature reveals little information regarding the role of the thiamine dependent TK reaction in tumor cell ribose production which is a central process in de novo nucleic acid synthesis and the salvage pathways for purines. Consequently, current thiamine administration protocols oversupply thiamine by 200% to 20,000% of the recommended dietary allowance, because it is considered harmless and needed by cancer patients. The thiamine dependent TK pathway is the central avenue which supplies ribose phosphate for nucleic acids in tumors and excessive thiamine supplementation maybe responsible for failed therapeutic attempts to terminate cancer cell proliferation. Limited administration of thiamine and concomitant treatment with transketolase inhibitors is a more rational approach to treat cancer.

Comparison of control analysis data using different approaches: modelling and experiments with muscle extract

Experimental and model studies have been performed to characterize the control properties of hexokinase and phosphofructokinase in muscle glycolysis and to examine the nature of error associated with experimental flux control coefficient determinations. Different approaches of metabolic control analysis, classical titration, co-response analysis and kinetic modelling indicated that flux control coefficients could be reliably estimated experimentally for the upper part of glycolysis. The kinetic parameters applied to construct the mathematical model were determined in muscle extract under similar conditions used for flux studies. If the kinetic parameters of commercial enzymes are introduced into the model the control analysis data cannot be trusted. Co-response analysis can also be successfully applied to determination of the flux control coefficients of the system. However, the involvement of a rapid-equilibrium enzyme, such as glucose 6-phosphate isomerase, could result in estimation errors for the relevant co-response coefficients that are propagated into the elasticity matrix. If the co-response coefficients related to isomerase activity are replaced by the values obtained by kinetic modelling, the values of elasticities are correct. Our data also suggest that in the upper part of glycolysis hexokinase mainly controls the pathway flux whereas phosphofructokinase exerts dominant control on the turnover of internal metabolite stocks inside the system.

Oxythiamine and dehydroepiandrosterone inhibit the nonoxidative synthesis of ribose and tumor cell proliferation

This study investigates the significance of the glucose-6-phosphate dehydrogenase (G6PD) catalyzed oxidative and the transketolase (TK) catalyzed nonoxidative pentose cycle (PC) reactions in the tumor proliferation process by characterizing tumor growth patterns and synthesis of the RNA ribose moiety in the presence of respective inhibitors of G6PD and TK. Mass spectra analysis of 13C-labeled carbons revealed that these PC reactions contribute to over 85% of de novo ribose synthesis in RNA from [1,2-(13)C]glucose in cultured Mia pancreatic adenocarcinoma cells, with the fraction synthesized through the TK pathway predominating (85%). Five days of treatment with the TK inhibitor oxythiamine (OT) and the G6PD inhibitor dehydroepiandrosterone-sulfate (0.5 microM each) exerted a 39 and a 23% maximum inhibitory effect on cell proliferation in culture, which was increased to 60% when the two drugs were administered in combination. In vivo testing of 400 mg/kg OT or dehydroepiandrosterone-sulfate in C57BL/6 mice hosting Ehrlich's ascitic tumor cells revealed a 90.4 and a 46% decrease in the final tumor mass after 3 days of treatment. RNA ribose fractional synthesis through the TK reaction using metabolites directly from glycolysis declined by 9.1 and 23.9% after OT or the combined treatment, respectively. Nonoxidative PC reactions play a central regulating role in the carbon-recruiting process toward de novo nucleic acid ribose synthesis and cell proliferation in vitro and in vivo. Therefore, enzymes or substrates regulating the nonoxidative synthesis of ribose could also be the sites to preferentially target tumor cell proliferation by new anticancer drugs.

The metabolic productivity of the cell factory

It is widely accepted that some performance function has been optimized during the evolution of metabolic pathways. One can study the nature of such a function by analogy with the industrial manufacturing world, in which there have been efforts over recent decades to optimize production chains, and in which it is now accepted that fluxes are not the only important system variables that determine process efficiency, because inventory turnover must also be considered. Inspired by the parallels between living cells and manufacturing factories, we propose that fluxes and transit time may have simultaneously been major targets of natural selection in the optimization of the design, structure and kinetic parameters of metabolic pathways. Accordingly we define the ratio of flux to transit time as a performance index of productivity in metabolic systems: it measures the efficiency with which stocks are administered, and facilitates comparison of a pathway in different steady states or in different tissues or organisms. For a linear chain of two enzymes, at a fixed total equilibrium constant, we have analysed the variation of flux, transit time and productivity index as functions of the equilibrium constants of the two steps. The results show that only the productivity index has a maximum, which represents a good compromise in optimizing flux and transit time. We have extended control analysis to the productivity index and derived the summation theorem that applies to it. For linear chains of different length with maximum productivity index values, the distribution of control coefficients with regard to the three parameters has a characteristic profile independent of the length of the chain. Finally, this control profile changes when other variables are optimized, and we compare the theoretical results with the control profile of the first steps of glycolysis in rat liver.

Strong control on the transit time in metabolic channelling

A suite of different characteristic times is used to describe the temporal behavior of a metabolic pathway. Here we focus on the 'transit' time, that is the average time it takes for a molecule, entering the steady-state pathway as a substrate, to exit the pathway as a product. We show that metabolic channelling results in dramatic changes in control exerted by pathway enzymes on the transit time. In an 'ideal' pathway a doubling of the enzyme concentrations halves the transit time. In a dynamic channel such an increase can reduce the transit time by a factor of four or more.

Control in channelled pathways. A matrix method calculating the enzyme control coefficients

The usual equations expressing the enzyme control coefficients (quantitative indicators of 'global' control properties of a pathway) via the elasticity coefficients (reflecting local kinetic properties of an enzyme reaction), cannot be applied to a variety of 'non-ideal' pathways, in particular to pathways with metabolic channelling. Here we show that the relationship between the control and elasticity coefficients can be obtained by considering such a metabolic pathway as a network of elemental chemical conversions (steps). To calculate the control coefficients of enzymes one should first determine the elasticity coefficients of such elemental steps and then take their appropriate combinations. Although the method is illustrated for a channelled pathway it can be used for any non-ideal pathway including those with high enzyme concentrations where the sequestration of metabolites by enzymes cannot be neglected.