Biomedical subjects

Marta Cascante

Publications and source records attributed to Marta Cascante.

At least 19 recordsLinked to original sources

PubMed2006

Pentose phosphate cycle oxidative and nonoxidative balance: A new vulnerable target for overcoming drug resistance in cancer.

The metabolic network of cancer cells confers adaptive mechanisms against many chemotherapeutic agents, but also presents critical constraints that make the cells vulnerable to perturbation of the network due to drug therapy. To identify these fragilities, combination therapies based on targeting the nucleic acid synthesis metabolic network at multiple points were tested. Results showed that cancer cells overcome single hit strategies through different metabolic network adaptations, demonstrating the robustness of cancer cell metabolism. Analysis of these adaptations also identified the maintenance of pentose phosphate cycle oxidative and nonoxidative balance to be critical for cancer cell survival and vulnerable to chemotherapeutic intervention. The vulnerability of cancer cells to the imbalance on pentose phosphate cycle was demonstrated by phenotypic phase plane analysis.

Animals↗

PubMed2006

Integration of enzyme kinetic models and isotopomer distribution analysis for studies of in situ cell operation.

A current trend in neuroscience research is the use of stable isotope tracers in order to address metabolic processes in vivo. The tracers produce a huge number of metabolite forms that differ according to the number and position of labeled isotopes in the carbon skeleton (isotopomers) and such a large variety makes the analysis of isotopomer data highly complex. On the other hand, this multiplicity of forms does provide sufficient information to address cell operation in vivo. By the end of last millennium, a number of tools have been developed for estimation of metabolic flux profile from any possible isotopomer distribution data. However, although well elaborated, these tools were limited to steady state analysis, and the obtained set of fluxes remained disconnected from their biochemical context. In this review we focus on a new numerical analytical approach that integrates kinetic and metabolic flux analysis. The related computational algorithm estimates the dynamic flux based on the time-dependent distribution of all possible isotopomers of metabolic pathway intermediates that are generated from a labeled substrate. The new algorithm connects specific tracer data with enzyme kinetic characteristics, thereby extending the amount of data available for analysis: it uses enzyme kinetic data to estimate the flux profile, and vice versa, for the kinetic analysis it uses in vivo tracer data to reveal the biochemical basis of the estimated metabolic fluxes.

Algorithms↗

PubMed2006

(2Alpha,3beta)-2,3-dihydroxyolean-12-en-28-oic acid, a new natural triterpene from Olea europea, induces caspase dependent apoptosis selectively in colon adenocarcinoma cells.

Triterpenoids are known to induce apoptosis and to be anti-tumoural. Maslinic acid, a pentacyclic triterpene, is present in high concentrations in olive pomace. This study examines the response of HT29 and Caco-2 colon-cancer cell lines to maslinic-acid treatment. At concentrations inhibiting cell growth by 50-80% (IC50HT29=61+/-1 microM, IC80HT29=76+/-1 microM and IC50Caco-2=85+/-5 microM, IC80Caco-2=116+/-5 microM), maslinic acid induced strong G0/G1 cell-cycle arrest and DNA fragmentation, and increased caspase-3 activity. However, maslinic acid did not alter the cell cycle or induce apoptosis in the non-tumoural intestine cell lines IEC-6 and IEC-18. Moreover, maslinic acid induced cell differentiation in colon adenocarcinoma cells. These findings support a role for maslinic acid as a tumour suppressant and as a possible new therapeutic tool for aberrant cell proliferation in the colon. In this report, we demonstrate for the first time that, in tumoural cancer cells, maslinic acid exerts a significant anti-proliferation effect by inducing an apoptotic process characterized by caspase-3 activation by a p53-independent mechanism, which occurs via mitochondrial disturbances and cytochrome c release.

Adenocarcinoma↗

PubMed2006

Software for dynamic analysis of tracer-based metabolomic data: estimation of metabolic fluxes and their statistical analysis.

MOTIVATION: Metabolic flux analysis of biochemical reaction networks using isotope tracers requires software tools that can analyze the dynamics of isotopic isomer (isotopomer) accumulation in metabolites and reveal the underlying kinetic mechanisms of metabolism regulation. Since existing tools are restricted by the isotopic steady state and remain disconnected from the underlying kinetic mechanisms, we have recently developed a novel approach for the analysis of tracer-based metabolomic data that meets these requirements. The present contribution describes the last step of this development: implementation of (i) the algorithms for the determination of the kinetic parameters and respective metabolic fluxes consistent with the experimental data and (ii) statistical analysis of both fluxes and parameters, thereby lending it a practical application. RESULTS: The C++ applications package for dynamic isotopomer distribution data analysis was supplemented by (i) five distinct methods for resolving a large system of differential equations; (ii) the 'simulated annealing' algorithm adopted to estimate the set of parameters and metabolic fluxes, which corresponds to the global minimum of the difference between the computed and measured isotopomer distributions; and (iii) the algorithms for statistical analysis of the estimated parameters and fluxes, which use the covariance matrix evaluation, as well as Monte Carlo simulations. An example of using this tool for the analysis of (13)C distribution in the metabolites of glucose degradation pathways has demonstrated the evaluation of optimal set of parameters and fluxes consistent with the experimental pattern, their range and statistical significance, and also the advantages of using dynamic rather than the usual steady-state method of analysis. AVAILABILITY: Software is available free from http://www.bq.ub.es/bioqint/selivanov.htm

Algorithms↗

PubMed2006

Mathematical modeling of polyamine metabolism in mammals.

Polyamines are considered as essential compounds in living cells, since they are involved in cell proliferation, transcription, and translation processes. Furthermore, polyamine homeostasis is necessary to cell survival, and its deregulation is involved in relevant processes, such as cancer and neurodegenerative disorders. Great efforts have been made to elucidate the nature of polyamine homeostasis, giving rise to relevant information concerning the behavior of the different components of polyamine metabolism, and a great amount of information has been generated. However, a complex regulation at transcriptional, translational, and metabolic levels as well as the strong relationship between polyamines and essential cell processes make it difficult to discriminate the role of polyamine regulation itself from the whole cell response when an experimental approach is given in vivo. To overcome this limitation, a bottom-up approach to model mathematically metabolic pathways could allow us to elucidate the systemic behavior from individual kinetic and molecular properties. In this paper, we propose a mathematical model of polyamine metabolism from kinetic constants and both metabolite and enzyme levels extracted from bibliographic sources. This model captures the tendencies observed in transgenic mice for the so-called key enzymes of polyamine metabolism, ornithine decarboxylase, S-adenosylmethionine decarboxylase and spermine spermidine N-acetyl transferase. Furthermore, the model shows a relevant role of S-adenosylmethionine and acetyl-CoA availability in polyamine homeostasis, which are not usually considered in systemic experimental studies.

Acetyl Coenzyme A↗

PubMed2006

K-ras Asp12 mutant neither interacts with Raf, nor signals through Erk and is less tumorigenic than K-ras Val12.

Different mutant amino acids in the Ras proteins lead to distinct transforming capacities and different aggressiveness in human tumors. K-Ras Asp12 (K12D) is more prevalent in benign than in malignant human colorectal tumors, whereas K-Ras Val12 (K12V) associates with more advanced and metastatic carcinomas, higher recurrence and decreased survival. Here, we tested, in a nude mouse xenograft model, whether different human K-Ras oncogenes mutated at codon 12 to Val, Asp or Cys would confer NIH3T3 fibroblasts distinct oncogenic phenotypes. We studied tumor histology and growth, apoptotic and mitotic rates, activation of signal transduction pathways downstream of Ras and regulation of the cell cycle and apoptotic proteins in tumors derived from the implanted transformants. We found that the K12V oncogene induces a more aggressive tumorigenic phenotype than the K12D oncogene, whereas K12C does not induce tumors in this model. Thus, K12V mutant tumors proliferate about seven times faster, and have higher cellularity and mitotic rates than the K12D mutant tumors. A molecular analysis of the induced tumors shows that the K12V mutant protein interacts with Raf-1 and transduces signals mainly through the Erk pathway. Unexpectedly, in tumors induced by the K12D oncogene, the K-Ras mutant protein does not interact with Raf-1 nor activates the Erk canonical pathway. Instead, it transduces signals through the PI3K/Akt, JNK, p38 and FAK pathways. Finally, the higher growth rate of the K12V tumors associates with enhanced Rb phosphorylation, and PCNA and cyclin B upregulation, consistent with faster G1/S and G2/M transitions, without alteration of apoptotic regulation.

Animals↗

PubMed2006

Electron-transfer capacity of catechin derivatives and influence on the cell cycle and apoptosis in HT29 cells.

Galloylated and nongalloylated catechin conjugates with cysteine derivatives have been synthesized and evaluated for their capacity to scavenge free radicals and to influence crucial functions (cell cycle, apoptosis) in HT29 colon carcinoma cells. We show that the nonphenolic part of the molecule modified the capacity of catechins to donate hydrogen atoms and to transfer electrons to free radicals. Nongalloylated derivatives did not significantly influence either the cell cycle or apoptosis. Among the galloylated species, 4beta-[S-(O-ethyl-cysteinyl)]epicatechin 3-O-gallate, which showed a high electron-transfer capacity (5 e- per molecule), arrested the cell cycle and induced apoptosis as expected for galloylated catechins such as tea (-)-epigallocatechin 3-O-gallate. 4beta-[S-(N-Acetyl-O-methyl-cysteinyl)]epicatechin 3-O-gallate, which showed the highest hydrogen-donating capacity (10 H per molecule) while keeping the electron-transfer capacity low (2.9 e- per molecule), did not trigger any significant apoptosis. The gallate moiety did not appear to be sufficient for the pro-apoptotic effect of the catechin derivatives in HT29 cells. Instead, a high electron-transfer capacity is more likely to be behind this effect. The use of stable radicals sensitive exclusively to electron transfer may help to design molecules with either preventive scavenging action (high hydrogen donation, low electron transfer) or therapeutic pro-apoptotic activity (high electron transfer).

Apoptosis↗

PubMed2006

Bistability from double phosphorylation in signal transduction. Kinetic and structural requirements.

Previous studies have suggested that positive feedback loops and ultrasensitivity are prerequisites for bistability in covalent modification cascades. However, it was recently shown that bistability and hysteresis can also arise solely from multisite phosphorylation. Here we analytically demonstrate that double phosphorylation of a protein (or other covalent modification) generates bistability only if: (a) the two phosphorylation (or the two dephosphorylation) reactions are catalyzed by the same enzyme; (b) the kinetics operate at least partly in the zero-order region; and (c) the ratio of the catalytic constants of the phosphorylation and dephosphorylation steps in the first modification cycle is less than this ratio in the second cycle. We also show that multisite phosphorylation enlarges the region of kinetic parameter values in which bistability appears, but does not generate multistability. In addition, we conclude that a cascade of phosphorylation/dephosphorylation cycles generates multiple steady states in the absence of feedback or feedforward loops. Our results show that bistable behavior in covalent modification cascades relies not only on the structure and regulatory pattern of feedback/feedforward loops, but also on the kinetic characteristics of their component proteins.

Animals↗

PubMed2005

Metabolic homeostasis in the human erythrocyte: in silico analysis.

A detailed computer model of human erythrocyte metabolism was shown to predict three steady states, two stable and one unstable. The most extreme steady state is characterized by almost zero concentrations of all the phosphorylated intermediates. The "normal" steady state is remarkably robust in the face of large changes in the activity of most of the enzymes of glycolysis and the pentose phosphate pathway: this steady state can be viewed as an attractor towards which the system returns following a metabolic perturbation. Focus is given to three responses of the system: (1) the 'energy charge' that pertains to the concentration of ATP relative to all purine nucleotides; (2) redox power expressed as the ratio of reduced-to-total glutathione and (3) the concentration of 2,3-bisphosphoglycerate, that directly affects the oxygen affinity of haemoglobin thus affecting the main physiological function of the cell. The collapse of the normal steady state in what can be viewed topologically as a catastrophe is posited as one key element of erythrocyte senescence and it is particularly important for erythrocyte destruction in patients with an inborn enzyme deficiency.

2,3-Diphosphoglycerate↗

PubMed2005

In situ localization of transketolase activity in epithelial cells of different rat tissues and subcellularly in liver parenchymal cells.

Metabolic mapping of enzyme activities (enzyme histochemistry) is an important tool to understand (patho)physiological functions of enzymes. A new enzyme histochemical method has been developed to detect transketolase activity in situ in various rat tissues and its ultrastructural localization in individual cells. In situ detection of transketolase is important because this multifunctional enzyme has been related with diseases such as cancer, diabetes, Alzheimer's disease, and Wernicke-Korsakoff's syndrome. The proposed method is based on the tetrazolium salt method applied to unfixed cryostat sections in the presence of polyvinyl alcohol. The method appeared to be specific for transketolase activity when the proper control reaction is performed and showed a linear increase of the amount of final reaction product with incubation time. Transketolase activity was studied in liver, small intestine, trachea, tongue, kidney, adrenal gland, and eye. Activity was found in liver parenchyma, epithelium of small intestine, trachea, tongue, proximal tubules of kidney and cornea, and ganglion cells in medulla of adrenal gland. To demonstrate transketolase activity ultrastructurally in liver parenchymal cells, the cupper iron method was used. It was shown that transketolase activity was present in peroxisomes and at membranes of granular endoplasmic reticulum. This ultrastructural localization is similar to that of glucose-6-phosphate dehydrogenase activity, suggesting activity of the pentose phosphate pathway at these sites. It is concluded that the method developed for in situ localization of transketolase activity for light and electron microscopy is specific and allows further investigation of the role of transketolase in (proliferation of) cancer cells and other pathophysiological processes.

Animals↗

PubMed2005

Effect of new antioxidant cysteinyl-flavanol conjugates on skin cancer cells.

Novel catechin derivatives obtained from grape procyanidins and l-cysteine scavenge free radicals by hydrogen atom donation, rather than electron transfer, and reduce cell viability in A375 and M21 melanoma cells. In particular, 4beta-(S-cysteinyl)epicatechin 3-O-gallate has a free radical scavenging capacity as strong as that of tea (-)-epigallocatechin gallate and causes a significant S-phase cell-cycle arrest in both cell lines at doses higher than 100 microM. The other cysteinyl compounds do not affect normal cell cycle distribution. The gallate derivative also induces apoptosis in melanoma cells more strongly than the other derivatives and the parent (-)-epicatechin do. The gallate compound seems to trigger nuclear condensation and fragmentation, which is confirmed by DNA laddering. Interestingly, they do not induce apoptosis in keratinocytes (HaCaT).

Antineoplastic Agents, Phytogenic↗

PubMed2005

Improved localization of glucose-6-phosphate dehydrogenase activity in cells with 5-cyano-2,3-ditolyl-tetrazolium chloride as fluorescent redox dye reveals its cell cycle-dependent regulation.

Since the introduction of cyano-ditolyl-tetrazolium chloride (CTC), a tetrazolium salt that gives rise to a fluorescent formazan after reduction, it has been applied to quantify activity of dehydrogenases in individual cells using flow cytometry. Confocal laser scanning microscopy (CLSM) showed that the fluorescent formazan was exclusively localized at the surface of individual cells and not at intracellular sites of enzyme activity. In the present study, the technique has been optimized to localize activity of glucose-6-phosphate dehydrogenase (G6PD) intracellularly in individual cells. Activity was demonstrated in cultured fibrosarcoma cells in different stages of the cell cycle. Cells were incubated for the detection of G6PD activity using a medium containing 6% (w/v) polyvinyl alcohol, 5 mM CTC, magnesium chloride, sodium azide, the electron carrier methoxyphenazine methosulphate, NADP, and glucose-6-phosphate. Before incubation, cells were permeabilized with 0.025% glutaraldehyde. Fluorescent formazan was localized exclusively in the cytoplasm of fibrosarcoma cells. The amount of fluorescent formazan in cells increased linearly with incubation time when measured with flow cytometry and CLSM. When combining the Hoechst staining for DNA with the CTC method for the demonstration of G6PD activity, flow cytometry showed that G6PD activity of cells in S phase and G2/M phase is 27 +/- 4% and 43 +/- 4% higher, respectively, than that of cells in G1 phase. CLSM revealed that cells in all phases of mitosis as well as during apoptosis contained considerably lower G6PD activity than cells in interphase. It is concluded that posttranslational regulation of G6PD is responsible for this cell cycle-dependent activity.

Animals↗

PubMed2005

Rapid simulation and analysis of isotopomer distributions using constraints based on enzyme mechanisms: an example from HT29 cancer cells.

MOTIVATION: Addition of labeled substrates and the measurement of the subsequent distribution of the labels in isotopomers in reaction networks provide a unique method for assessing metabolic fluxes in whole cells. However, owing to insufficiency of information, attempts to quantify the fluxes often yield multiple possible sets of solutions that are consistent with a given experimental pattern of isotopomers. In the study of the pentose phosphate pathways, the need to consider isotope exchange reactions of transketolase (TK) and transaldolase (TA) (which in past analyses have often been ignored) magnifies this problem; but accounting for the interrelation between the fluxes known from biochemical studies and kinetic modeling solves it. The mathematical relationships between kinetic and equilibrium constants restrict the domain of estimated fluxes to the ones compatible not only with a given set of experimental data, but also with other biochemical information. METHOD: We present software that integrates kinetic modeling with isotopomer distribution analysis. It solves the ordinary differential equations for total concentrations (accounting for the kinetic mechanisms) as well as for all isotopomers in glycolysis and the pentose phosphate pathway (PPP). In the PPP the fluxes created in the TK and TA reactions are expressed through unitary rate constants. The algorithms that account for all the kinetic and equilbrium constant constraints are integrated with the previously developed algorithms, which have been further optimized. The most time-consuming calculations were programmed directly in assembly language; this gave an order of magnitude decrease in the computation time, thus allowing analysis of more complex systems. The software was developed as C-code linked to a program written in Mathematica (Wolfram Research, Champaign, IL), and also as a C++ program independent from Mathematica. RESULTS: Implementing constraints imposed by kinetic and equilibrium constants in the isotopomer distribution analysis in the data from the cancer cells eliminated estimates of fluxes that were inconsistent with the kinetic mechanisms of TK and TA. Fluxes measured experimentally in cells can be used to estimate better the kinetics of TK and TA as they operate in situ. Thus, our approach of integrating various methods for in situ flux analysis opens up the possibility of designing new types of experiments to probe metabolic interrelationships, including the incorporation of additional biochemical information. AVAILABILITY: Software is available freely at: http://www.bq.ub.es/bioqint/selivanov.htm CONTACT: martacascante@ub.edu

Algorithms↗

PubMed2005

K-ras codon-specific mutations produce distinctive metabolic phenotypes in NIH3T3 mice [corrected] fibroblasts.

Among K-ras mutations, codon 12 mutations have been identified as those conferring a more aggressive phenotype. This aggressiveness is primarily associated with slow proliferation but greatly increased resistance to apoptosis. Using transfected NIH3T3 fibroblasts with a mutated K-ras minigene either at codon 12 (K12) or at codon 13 (K13), and taking advantage of [1,2-13C2]glucose tracer labeling, we show that codon 12 mutant K-ras (K12)-transformed cells exhibit greatly increased glycolysis with only a slight increase in activity along pathways that produce nucleic acid and lipid synthesis precursors in the oxidative branch of the pentose phosphate pathway and via pyruvate dehydrogenase flux. K13 mutants display a modest increase in anaerobic glycolysis associated with a large increase in oxidative pentose phosphate pathway activity and pyruvate dehydrogenase flux. The distinctive differences in metabolic profiles of K12 and K13 codon mutated cells indicate that a strong correlation exists between the flow of glucose carbons towards either increased anaerobic glycolysis, and resistance to apoptosis (K12), or increased macromolecule synthesis, rapid proliferation, and increased sensitivity to apoptosis.

Anaerobiosis↗

PubMed2005

Procyanidin fractions from pine (Pinus pinaster) bark: radical scavenging power in solution, antioxidant activity in emulsion, and antiproliferative effect in melanoma cells.

Pine (Pinus pinaster) bark is a rich source of procyanidin oligomers. From a total polyphenolic extract, we have generated fractions of different procyanidin composition. The mixtures, devoid of gallate esters, were active as free radical scavengers against ABTS(*+), DPPH, and HNTTM. Pine bark fractions were tested for antioxidant activity in solution (hydrogen donation and electron transfer) and emulsion (inhibition of lipid peroxidation) and compared with their galloylated counterparts from grape origin. While galloylation clearly influenced the free radical scavenging efficiency in solution, it did not seem to play a determinant role in protection against lipid peroxidation in emulsion. The fractions were very mild inhibitors of cell proliferation. Because gallate esters appear to interfere with crucial cell functions, gallate free pine procyanidins may be the innocuous chemopreventative agents of choice for many applications in food and skin protection.

Antioxidants↗

PubMed2005

DNA damage induced by a quinoxaline 1,4-di-N-oxide derivative (hypoxic selective agent) in Caco-2 cells evaluated by the comet assay.

The DNA damage induced by 7-chloro-3-[[(N,N-dimethylamino)propyl]amino]-2-quinoxalinecarbonitrile 1,4-di-N-oxide hydrochloride (Q-85 HCl) in Caco-2 cells under hypoxic and well-oxygenated conditions has been studied by using the comet assay. This compound has shown a good in vitro profile of high selective toxicity in hypoxia, but its mechanism of action is unknown. The DNA damage has been evaluated by performing the comet assay after a 2-h treatment with Q-85 HCl (0.1, 0.2, 0.4 microM in hypoxia; 20, 40 microM in well-oxygenated conditions). The number of cells in apoptosis has also been assessed by flow cytometry analysis of Annexin V-FITC staining. The capability of the cells to repair the DNA damage and the proliferation rate was evaluated at different times after the treatment (24-168 h). Under hypoxic conditions, a clear dose-dependent increase in the number of nuclei with a comet was observed (comet score: 132 +/- 13, 343 +/- 30 and 399 +/- 1; control comet score: 42 +/- 14). Under well-oxygenated conditions, the number of nuclei with comet increased significantly with respect to the control (comet score: 273 +/- 14 and 312 +/- 9; control comet score: 27 +/- 4). Cells in apoptosis were not detected by the comet assay nor by flow cytometry. The recovery from DNA damage was time- and concentration-dependent in hypoxia (cells treated with the highest concentration still showed DNA damage after 72 h) and rather time-dependent in well-oxygenated conditions (DNA was completely repaired after 24 h). In conclusion, Q-85 HCl acts by DNA damage and not only the reduced intermediate is genotoxic but also some other derivatives and Q-85 HCl itself may be acting.

Apoptosis↗

PubMed2005

Selective toxicity of a quinoxaline 1,4-di-N-oxide derivative in human tumour cell lines.

The presence of hypoxic cells in human solid tumours is one of the causes of tumour resistance to conventional therapy, and is also associated with processes that promote the tumour progression. Different chemical agents have been designed in order to take advantage of the particular metabolic characteristics of hypoxic regions. These drugs, called bioreductive agents, are activated inside the hypoxic cells to give active species that, in the presence of oxygen, are oxidised back to the non-toxic parent compound. Several quinoxaline 1,4-di-N-oxides have been described as potential bioreductive agents, and among them, 7-cloro-3-[[(N,N-dimethylamino)propy]amino]-2-quinoxalinecarbonitrile 1,4-di-N-oxide hydrocloride (Q-85 HCl) appeared to be the most promising one. In the present work, the selective cytotoxicity of Q-85 HCl was studied in several human tumour cell lines of different origin (Caco2, MCF-7, HT-29 and Tk-10). Cell viability was calculated after 2 h treatment under hypoxic and well-oxygenated conditions. The potency (the concentration that gives 1% of cell survival) in hypoxia and hypoxia cytotoxicity ratio (HCR = potency in oxygenated conditions/potency in hypoxia) were calculated after a 14-day clonogenic assay. Q-85 HCl was more toxic in hypoxia than in well-oxygenated cells in all the tumour cell lines. The best profile of potency in hypoxia (0.4 micromol/L) and selectivity (HCR=155) was found in CaCo-2 cells. Altogether, these results suggest an in vitro biological profile for Q-85 HCl that makes it an interesting candidate for the development as a bioreductive agent.

Antineoplastic Agents↗

PubMed2004

An optimized algorithm for flux estimation from isotopomer distribution in glucose metabolites.

MOTIVATION: Analysis of the conversion of (13)C glucose within the metabolic network allows the evaluation of the biochemical fluxes in interconnecting metabolic pathways. Such analyses require solving hundreds of equations with respect to individual isotopomer concentrations, and this assumes applying special software even for constructing the equations. The algorithm, proposed by others could be improved. METHOD: A C-code linked to the program written in Mathematica (Wolfram Research Inc.), constructs and solves differential equations for all isotopomer concentrations, using the general enzyme characteristics (K(m), equilibrium constant, etc.). This code uses innovative algorithm of determination for the isotopomers-products, thus essentially decreasing the computation time. Feasible metabolic fluxes are provided by the parameters of enzyme kinetics found from the data fitting. RESULTS: The software effectively evaluates metabolic fluxes based on the measured isotopomer distribution, as was illustrated by the analysis of glycolysis and pentose phosphate cycle. The mechanism of transketolase and transaldolase catalysis was shown to induce a specific kind of isotopomer re-distribution, which, despite the significance of its effect, usually is not taken into account. AVAILABILITY: The software could be freely downloaded from the site: http://bq.ub.es/bioqint/label_distribution/.

Algorithms↗