BACKGROUND
obesity is increasing to epidemic proportions  <cit> , <cit> , starting from young ages  <cit> , <cit> , and is associated with an increase in the prevalence of type- <dig> diabetes mellitus  throughout the world  <cit> , <cit> .
in recent years, active participation of immune cells in obesity and t2dm has become evident  <cit> , <cit> .
this chronic activation of the innate immune system can lead to insulin resistance , impaired tolerance to glucose and, eventually, t2dm  <cit> , <cit> .
adipose tissue produces a vast array of adipocyte-derived factors  that regulate the metabolism, inflammation and body mass  <cit> , <cit> .
one of these adipokines, resistin, belongs to a family of cysteine-rich proteins shown to be involved in inflammation and altered insulin sensitivity in rodents  <cit> , <cit> .
resistin is produced from adipose tissue in rodents whereas, in humans, it is secreted by the mononuclear lymphocytes and stromal cells within adipose tissue  <cit> , <cit> .
studies in humanised resistin mice suggest that, even though the site of resistin production differs between species, human resistin exacerbates inflammation in white adipose tissue and contributes to ir, thereby impairing its normal effects  <cit> , <cit> .
ir is defined as the genetic or learned inability of target tissues to respond normally to the action of circulating hormones.
it has been described in skeletal muscle but also in liver and adipose tissues  <cit> .
some studies have shown that resistin affects glucose transport and insulin-stimulated oxidation of glucose in l <dig> skeletal muscle cells  <cit> - <cit> .
it also decreases the uptake and oxidation of long-chain fatty acids  <cit>  and glycogen synthase kinase-3-β, as well as insulin-stimulated insulin receptor substrate- <dig>  tyrosine phosphorylation  <cit> , <cit>  in the same cell line.
furthermore, high levels of resistin in rats leads to ir involving impaired insulin signalling in skeletal muscle, liver and adipose tissues, resulting in glucose intolerance, hyperinsulinemia and hypertriglyceridemia  <cit> .
these observations link resistin in muscle-cell metabolism to ir.
changes in the metabolome  have been described as the “ultimate” response of an organism to various events, such as genetic alterations as well as disease-based or environmental influences  <cit> .
however, the metabolome is dynamic and metabolites are transformed continuously in the cell.
comprehensive characterisation of the metabolic networks and their functional states requires quantitative knowledge of intracellular metabolic fluxes.
these intracellular fluxes can be quantified by analysing incorporation of labelled substrates in metabolic products using appropriate bioinformatic tools  <cit> .
skeletal muscle is considered to be the main tissue involved in the maintenance of glucose homeostasis because its contribution to glucose uptake is ≈ 75% of the total contribution of peripheral tissues  <cit> , and it is the main tissue responsible for insulin-dependent glucose use.
also, pyruvate homeostasis-related fluxes in skeletal muscle play a critical part in glucose homeostasis  <cit> - <cit> .
a better understanding of how these fluxes work together in muscle cells is needed to understand this process.
the aim of the present study was to examine the effect of resistin on the central carbon metabolic network of rat skeletal muscle cells, and to ascertain how resistin alters the response of muscle cells to insulin at the metabolic level.
we used -glucose as a tracer, analysed isotopomer distributions by gas cromatography  coupled with mass spectrometry , and quantified metabolic fluxes using isodyn software  <cit> , <cit> - <cit> .
in the model of the glucose metabolic network in skeletal muscle we included the glycolytic pathway and pentose phosphate pathway .
pyruvate dehydrogenase complex -, phosphoenolpyruvate carboxykinase - and pyruvate carboxylase -catalysed reactions were also included to account for pyruvate homeostasis-related fluxes.
resistin-induced alterations in the profile of metabolic fluxes of rat skeletal muscle cells were identified, and were complemented by analyses of the expression of certain genes.
the study presented here offers a wider vision than previously published of the metabolic reorganisation that resistin induces in rat skeletal muscle cells and in the metabolic response of these cells to insulin.
our observations provide new data for a more precise understanding of the effect of resistin on metabolic-network adaptations of skeletal muscle cells, and help to explain the role of resistin in the ir described in myocytes.
our work reveals the importance of tackling the study of complex biological systems from a systems-biology point of view to fully understand multifactorial diseases and identify new therapeutic targets.
<dig>
RESULTS
for analyses of the effect of resistin on glucose metabolism and on the metabolic response of muscle cells to insulin, l6e <dig> myotubes were pretreated for 8 h with or without resistin .
after preincubation, myotubes were incubated for a further 6 h with glucose  that was 50% enriched in -glucose in the absence or presence of  <dig> nm insulin and in the continued presence or absence of resistin.
after 6 h of incubation, glucose consumption, lactate production, intracellular levels of glycogen and glucose-6-phosphate  were determined , as were the mass-isotopomer distributions in glucose and lactate from the incubation medium, as well as glycogen glucose and ribose isolated from rna in cell pellets .
analyses of the obtained data by isodyn determined the distribution of the metabolic flux profile of l6e <dig> myotubes under various incubation conditions .
figure  <dig> summarises the qualitative changes in metabolites and fluxes.
l6e <dig> myotubes were treated  or not treated  with  <dig> nm resistin for 8 h and then incubated for 6 h in the absence  or presence  of  <dig> nm insulin.
glucose consumption, lactate production, and intracellular contents of glycogen and glucose-6-phosphate were determined as described in the methods section.
results are the mean ± standard deviation .
units of glucose consumption and lactate production refer to the volume of the incubation medium.
indicates p <  <dig>  and  p <  <dig>  compared with the condition of neither resistin nor insulin.
indicates p <  <dig>  compared with the l6e <dig> myotubes treated only with resistin.
lactate
glycogen
ribose
l6e <dig> myotubes were treated  or not treated  with  <dig> nm resistin for 8 h and then incubated for 6 h with 10 mm glucose, 50%-enriched in -glucose in the absence  or presence  of  <dig> nm insulin.
mass isotopomer distributions at the end of incubations were determined, as described in the methods section.
results are the mean ± standard deviation .
indicates p <  <dig>  compared with the condition of neither resistin nor insulin.
glycolytic and ppp fluxes
glycogen-related fluxes
pyruvate homeostasis-related fluxes
tricarboxylic acid cycle fluxes
principal fluxes of central carbon metabolism calculated in l6e <dig> myotubes treated  or not treated  with  <dig> nm resistin for 8 h and then incubated for 6 h with 10 mm glucose, 50%-enriched in -glucose in the absence  or presence  of  <dig> nm insulin.
glucose phosphorylation flux was fixed according to the glucose consumption measured .
remaining fluxes are expressed by median values as well as minimum and maximum values from the  <dig> best flux sets.
fluxes are expressed in nmol · ml–1 · min– <dig>  where volume is the volume of the incubation medium.
χ <dig> values for the averages were: ins– res–:  <dig> ; ins– res+:  <dig> ; ins + res–:  <dig> ; ins + res+:  <dig> .
indicates that the flux is different to the control condition with a 99% confidence interval ,  indicates that the flux is different to the control condition with a 95% ci, and  indicates that the flux differs from the condition in which cells were treated only with resistin with a 99% ci.
ppp, pentose phosphate pathway; oaa, oxaloacetate.
<dig>  basal metabolic flux profile of l6e9
isodyn quantified the metabolic-flux distribution in l6e <dig> myotubes incubated in the presence of glucose  and in the absence of insulin or resistin .
under this incubation condition, cells had a glucose phosphorylation of  <dig>  nmol · ml–1 · min– <dig>  this flux was maintained throughout the subsequent steps of “upper glycolysis”  until the split of fructose- <dig> -bisphosphate into two triose-phosphates by aldolase.
from this reaction onwards, the glycolytic flux  increased twofold .
additionally, the flux across the oxidative branch of the ppp was around eight-times lower than the upper glycolytic flux.
with regard to glycogen metabolism, glycogen phosphorylase   and glycogen synthase   were active, and glycogen recycling was estimated to be ≈ 2% of the total glucose phosphorylation flux.
however, the most important contribution of isodyn for interpretation of experimental data was the possibility of flux quantification around the final steps of lower glycolysis and the tricarboxylic acid  cycle.
in that regard, isodyn revealed that ≈ 49% of phosphoenolpyruvate  and pyruvate entered the tca cycle through the reactions catalysed by pdc  or pc , with a relative pdc:pc ratio of 8: <dig>  around 44% of the pep–pyruvate flux was diverted to the synthesis and output of lactate towards the cell culture medium , and pyruvate cycling flux .
this finding accounted for the fact that feedback of molecules from the tca cycle back towards the pep–pyruvate pool was around ten-times higher than the flux of the tca cycle towards citrate synthesis .
<dig>  effect of resistin on the metabolic flux profile of l6e9
l6e <dig> myotubes treated with  <dig> nm resistin showed 13% higher flux of glucose phosphorylation .
fitting of the experimental data from this incubation condition using isodyn  suggested that resistin increased the glycolytic flux, with a significant increase in the glucose phosphorylation flux, and a proportional increase of glycolytic fluxes.
resistin induced changes in the non-oxidative branch of the ppp , but the relative flux through the oxidative branch  was unchanged.
glycogen recycling-related fluxes  remained unchanged, as indicated by the fact that the net 13c incorporation in glycogen was very similar to the control condition.
moreover, resistin increased lactate production , which translated to an increase in the lactate dehydrogenase flux , though it was not significant.
the most significant metabolic change after isodyn fitting was observed in the pdc-related flux , which was increased significantly .
if we consider the relative ratios of the two main pyruvate incorporation fluxes , there was an increase of the pdc prevalence to 20: <dig>  this finding was also supported by the level of pyruvate dehydrogenase kinase isozyme  <dig>  mrna , and was observed as a trend in the level of pyruvate dehydrogenase pdha <dig> mrna .
despite the large change in this ratio, the effect on fluxes in the tca cycle was almost negligible even though there was an increase in the acetyl-coa output flux .
the pyruvate cycling flux was unchanged by resistin treatment alone, which was supported by the level of pepck mrna .
<dig>  effect of resistin on the metabolic response of insulin in l6e <dig> myotubes
having investigated the alterations caused by resistin on the metabolic profile of l6e <dig> myotubes, we examined its effect on the metabolic response of myotubes to insulin.
first we ascertained if resistin affected stimulation of insulin-induced glucose transport.
thus, induction of 2-deoxy-d--glucose  transport was assayed 30 min after incubation with  <dig> nm insulin in l6e <dig> myotubes treated or not treated with  <dig> nm resistin.
insulin induced an increase in 2-dg transport in l6e <dig> myotubes of  <dig>  ±  <dig> %.
however, it induced an increase in glucose transport in cells that had been pre-treated with  <dig> nm resistin for 8 h of  <dig>  ±  <dig> %.
subsequently, we analysed changes in the metabolic flux profile induced by insulin in l6e <dig> myotubes treated or not treated with  <dig> nm resistin.
in the absence of resistin, insulin stimulated an increase in glucose consumption and lactate production in l6e <dig> myotubes.
there were no changes in glycogen content, but a 150% increase in intracellular g6p content .
fitting of experimental data using isodyn  revealed that the flux of lactate production  was 48% higher than that in cells not treated with insulin.
insulin also increased the flux of glucose phosphorylation  and glycolytic fluxes  by ≈ 59%, and induced a 25% increase in the oxidative ppp flux .
moreover, fluxes related to glycogen metabolism  did not change significantly after insulin treatment .
insulin increased the flux through pdc  by around  <dig> -fold, but did not affect the remainder of tca cycle-related fluxes significantly, as confirmed by the reduction in levels of pdk <dig> mrna and the increasing tendency of levels of pdha <dig> mrna .
the relative ratio of pdc:pc fluxes increased from 8: <dig> to 13: <dig> after insulin treatment.
insulin also led to a 59% increase in the acetyl-coa output fluxes.
in the presence of resistin, insulin increased the flux of glucose phosphorylation by 16% more than when cells were treated with resistin alone.
glycolysis fluxes did not increase significantly compared with the resistin condition, but the lactate dehydrogenase flux increased by 45% .
the increase induced by resistin treatment in the oxidative ppp pathway was not significantly altered if cells were also treated with insulin.
similarly, the gs and gp fluxes did not change significantly after resistin treatment.
interestingly, there was no change in the pdc flux  or the acetyl-coa production flux , and no change was observed in the remainder of tca cycle-related fluxes.
if the ratio between pdc and pc fluxes is taken into account, it shifted from 20: <dig> towards 70: <dig> after resistin was combined with insulin.
the final metabolic state of insulin-stimulated l6e <dig> did not match the one of cells treated with resistin.
some fluxes were equal in resistin-treated and -untreated cells: oxidative ppp flux, gs flux, gp flux, and tca cycle-related fluxes.
most of the quantified fluxes were lower in resistin, but significant changes were found only in the lactate dehydrogenase flux, and the pdc:pc ratio.
<dig>  effect of resistin on gene expression of the enzymes involved in the pdc, pc and pyruvate cycling fluxes in l6e9
to investigate further if the observed effects in metabolic-flux distribution were prompted by a change in the gene expression of any of the enzymes involved in their activities, we analysed the gene expression of several enzymes involved in pyruvate homeostasis at the end of the four incubation conditions tested.
neither insulin, nor resistin, nor combination of the two hormones had a significant effect on levels of pyruvate carboxylase  mrna or pyruvate kinase muscle isozyme  mrna .
however, insulin  increased levels of pyruvate dehydrogenase e <dig>  and decreased pyruvate dehydrogenase kinase isozyme  <dig>  mrna .
presence of resistin did not alter the decrease in levels of pdk <dig> mrna caused by insulin, but prevented the increase in levels of pdha <dig> mrna.
mrna expression of pyruvate dehydrogenase kinase isozyme  <dig>  decreased after treatment with both hormones , with an additive effect if both were combined.
these results suggest that resistin prevents the insulin-induced increase in pdha <dig> expression, and does not prevent the increase in pdc activity caused by a reduction of expression of pdk <dig> and pdk <dig> genes.
moreover, mrna levels of cytosolic nadp-isocitrate dehydrogenase  and nadp-dependent malic enzyme , both causes of pyruvate cycling increasing nadph production, were not affected by any of the hormones tested under these conditions .
however, insulin significantly decreased  levels of pck <dig> mrna in the absence of resistin, but not in the presence of resistin .
these results suggest that resistin prevents inhibition of expression of the pck <dig> gene.
hence, cytosolic flux in pepck may be active in the presence of resistin independent of the presence of insulin.
<dig> discussion
after the discovery of several adipocyte-secreted molecules , adipocytes were recognized as part of endocrine tissue  <cit> , <cit> , <cit> .
these molecules enable adipocytes to communicate with other tissues and organs, and to regulate: the metabolism of lipids and glucose; energy balance; insulin action; cell proliferation.
resistin was identified as an adipokine, and it was suggested that it might link obesity and ir  <cit> , <cit> .
however, the role of resistin in the pathophysiology of ir in humans and animals, and how it acts in muscle, liver and fat, is controversial  <cit> , <cit> .
here, we investigated the effect of resistin on basal and insulin-stimulated glucose metabolism in l6e <dig> rat myotubes.
we used tracer-based metabolomics and our in-house isodyn software to analyse quantitatively metabolic flux distribution in this cell line under different incubation conditions.
the response to insulin and other hormone stimuli in l6e <dig> cells is similar to that observed in skeletal muscle in vivo <cit> .
hence, l6e <dig> cells are considered to be suitable models for analyses of the effect of resistin on glucose metabolism and on the metabolic response of muscle cells to insulin.
analyses with isodyn data provided an overview of the central carbon metabolism of l6e <dig> skeletal muscle cells.
analyses were undertaken assuming a metabolic steady state , and rapid mixing of isotopic isomers of the same metabolite throughout different intracellular compartments.
incubation with a large excess of glucose ensured the required glucose uptake throughout the incubation.
the other conditions of incubation did not change, so the assumption of a metabolic steady state seemed reasonable.
rapid mixing of isotopic isomers of the same species throughout the cell is accepted for analyses of data of stable isotopes  <cit> .
fitting of experimental data in the framework of these assumptions suggested a set of metabolic fluxes underlying the measured distributions of isotopic isomers.
<dig>  the effect of insulin on l6e <dig> myotubes increases glycolysis, ppp, and incorporation of pyruvate to the tca cycle through pdc and fatty-acid synthesis
our results suggested that l6e <dig> myotubes were highly glycolytic, with active synthesis and degradation of glycogen, and showed tca-cycle activity and acetyl-coa output.
oxidative ppp was weakly active in this cell type.
furthermore, use of the same tools revealed that insulin activated the uptake and phosphorylation of glucose, glycolysis and lactate production, the oxidative branch of the ppp, pyruvate incorporation towards the tca cycle and increased acetyl-coa output .
we observed that insulin altered the ratio between pdc and pc fluxes by increasing the pdc flux through an increase of levels of pdha <dig> and decrease of pdk <dig> and pdk <dig> mrna.
the effects of insulin on the metabolic flux distribution of l6e <dig> myotubes reported in the present study correlate with the known activation of pdc by profound suppression of pdk <dig> expression in skeletal muscle  <cit> , inhibition of expression of the pck <dig> gene, and the increase in glucose uptake and lactate production observed using other approaches  <cit> , <cit> .
<dig>  the effect of resistin on l6e <dig> myotubes increases glycolysis and fatty-acid synthesis and alters the ratio between pdc and pc, tilting it towards direct incorporation of pyruvate through pdc
our results revealed that, in the absence of insulin, resistin activated glucose phosphorylation and the glycolytic pathway.
also, resistin did not affect the glycogen content in skeletal muscle.
even though enrichment of the glycogen label  was not significantly different to the control condition, the tendency was for glycogen-related fluxes to be lower, a result that was in accordance with other works  <cit> .
furthermore, our approach showed that, even though many of the tca-cycle fluxes remained unchanged, the ratio between the pdc and pc fluxes was significantly altered, effectively increasing the incorporation of pyruvate through pdc into mitochondria .
pc activity is finely regulated by mechanisms such as allosteric activation by acetyl carnitine, which was not tested in this study.
the acetyl-coa output for fatty-acid synthesis was also increased after resistin treatment, which would explain the higher flux towards the tca cycle.
<dig>  the effect of resistin impairs normal insulin metabolic actions
we showed that resistin affects the normal metabolic response of l6e <dig> myotubes to insulin.
thus, the increase in glucose uptake and glycolytic flux induced by insulin was lower than that observed in resistin-treated myotubes.
these results are in accordance with those published by other authors for l <dig> muscle cells  <cit> , <cit> , <cit> .
resistin treatment, however, did not affect the insulin-stimulated content and recycling of glycogen.
moreover, in resistin-treated cells, insulin did not increase the pdc flux.
after resistin treatment, levels of pdk <dig> mrna and pdk <dig> mrna were lower, accompanied with a reduction in levels of pdha <dig> mrna.
the reduction in expression of the main enzyme and key regulators of its activity resulted in a null net effect on the final flux.
the ratio between the pdc and pc fluxes was also altered, showing reduction in the pc anaplerotic entrance into the cycle.
fluxes through the tca cycle and acetyl-coa output fluxes were not increased, but instead maintained or reduced, by insulin in resistin-treated cells.
this lack of metabolic activation in response to insulin can be explained by the slight decrease observed in glucose uptake and the increase of lactate dehydrogenase flux directed to lactate production, which impaired the increase in the levels of glycolytic intermediates that could be metabolized through other metabolic pathways.
our results provide mechanistic explanations for the observations of palanivel and co-workers  <cit> , who reported that the insulin-stimulated oxidation of glucose via the tca cycle was reduced by resistin, and that this reduction would compromise oxidative respiration in the cell.
<dig>
CONCLUSIONS
several studies have identified positive correlations between resistin levels and ir in vivo <cit>  and in vitro <cit> - <cit> .
overall, our results showed that resistin significantly alters glucose metabolism in l6e <dig> myotubes and their response to insulin.
of the whole set of metabolic changes induced by resistin in l6e <dig> myotubes, the imbalance in pdc and pc fluxes as well as the higher acetyl-coa output flux for fatty-acid synthesis and other metabolic processes should be emphasized.
resistin impairs the normal metabolism of insulin, and the metabolism is directed to try to maintain the imbalance between the pdc and pc fluxes.
in this way, l6e <dig> reduces tca-cycle fluxes and acetyl-coa output fluxes, which leads to an important metabolic imbalance in the central carbon metabolic pathways of l6e <dig> myotubes.
it has been reported that ir in skeletal muscle, liver and adipose tissue is usually accompanied by glucose intolerance and hypertriglyceridemia  <cit> .
in this way, our results showed the importance of the balance of pdc and pc in glucose metabolism.
our results support the idea that metabolic disorders such as t2dm are caused by complex multi-molecular interactions that cannot be explained readily by an alteration in expression of a single gene or gene product, or even alteration of a single enzymatic cascade.
the robustness of results obtained on the effects of insulin on l6e <dig> metabolism found by combining tracer-based metabolomic data and flux analyses using isodyn showed that this approach is a suitable tool to study the effects of hormones on the central carbon metabolism network in myocytes.
<dig> methods
<dig>  cells and reagents
myoblastic cell line l6e <dig> was kindly provided by dr. nadal-ginard .
rat recombinant resistin was produced in a hek <dig> cell line.
the purity of resistin was controlled by sodium dodecyl sulfate–polyacrylamide gel electrophoresis  and estimated to be >99%.
biological activity was validated as described  <cit> .
resistin was tested for endotoxin using a toxinsensor chromogenic lal endotoxin assay kit .
this assay indicated low levels of endotoxins  at the final concentration used, which is below current usa food and drug administion  limits that require eluates from medical devices to be < <dig>  eu/ml.
-d-glucose and -d-glucose were purchased from isotec  and 2-deoxy-d--glucose  from ge healthcare .
<dig>  conditions for cell culture and incubation
myoblastic cell line l6e <dig> was grown as described  <cit> .
for biochemical and mass isotopomer assays, l6e <dig> myoblasts were seeded at  <dig> cells/cm <dig> in 60-mm diameter dishes with 4 ml dulbecco’s modified eagle medium  without glucose or glutamine  supplemented with 10% heat-inactivated foetal bovine serum  , 10 mm glucose, 4 mm asparagine and 5% penicillin/streptomycin .
cell differentiation was induced by lowering fbs to a final concentration of 2% when myoblasts were pre-confluent  <cit> .
cells at day- <dig> of differentiation were washed four times in phosphate-buffered saline  and pre-treated or not with  <dig> nm resistin dissolved in dmem without fbs and supplemented with 10 mm glucose, 4 mm asparagine, 2% hepes  and 1% penicillin/streptomycin.
after 8 h, the medium was replaced by another with the same composition, but now with 10 mm glucose 50%-enriched in -glucose in the presence or absence of  <dig> nm insulin .
at the end of 6-h incubation, media were collected and dishes washed twice with pbs before freezing in liquid nitrogen.
incubation media and dishes were stored at –80°c until biochemical and mass isotopomer analyses.
<dig>  glucose transport assay
to carry out the glucose transport assay,  <dig> l6e <dig> cells/dish were cultured on six-well plates.
at day- <dig> of differentiation, cells were pre-treated or not with  <dig> nm resistin for 8 h, as described above.
next, uptake of 100 μci 2-dg by myotubes incubated with or without  <dig> nm insulin for 30 min was measured as described  <cit> .
<dig>  biochemical determinations of glucose, lactate and g6p
concentrations of glucose and lactate in incubation media at the beginning and end of 6-h incubations were measured as described previously by spectrophotometric methods coupled to a cobas mira plus chemical analyzer .
intracellular g6p content was determined by enzymatic means as described  <cit> .
<dig>  mass isotopomer analyses of lactate, medium and glycogen glucose, and rna ribose, and determination of glycogen content by gc/ms
mass spectral data were obtained on a qp <dig> shimadzu mass selective detector connected to a gc- <dig> gas chromatograph  using helium as the carrier gas and isobutane  <dig>  pa as reagent gas in chemical ionization analyses.
the settings were: gc inlet, 250°c for glucose, ribose and deoxyribose, and 200°c for lactate; transfer line, 250°c; ms source, 200°c.
a varian vf- <dig> capillary column  was used to analyze all compounds.
glucose from the frozen media was purified using a tandem set of dowex-1x8/dowex-50wx <dig>  ion-exchange columns.
the glucose from media was converted to its glucose aldonitrile pentaacetate derivative as described  <cit> .
ion clusters around m/z  <dig>  were monitored.
lactate from the media was extracted with ethyl acetate after acidification with hydrochloric acid  and derivatized to its propylamide-heptafluorobutyric form.
clusters around m/z  <dig>  were monitored  <cit> .
rna ribose was isolated by acid hydrolysis of cellular rna after trizol purification of cell extracts.
ribose isolated from rna was derivatized to its aldonitrile-acetate form using hydroxylamine in pyridine and acetic anhydride.
we monitored ion clusters around m/z  <dig>  to find the molar enrichment and positional distribution of 13c labels in ribose as described  <cit> .
glycogen in frozen cell monolayers from different analytical conditions was extracted and quantified as described  <cit> .
extraction was done by direct digestion of sonicated extracts with amyloglucosidase.
mass isotopomer analyses of glycogen glucose was undertaken as described above for medium glucose.
measurement of glycogen content was carried out using the isotopomer -glucose as recovery standard and internal standard quantification procedures.
the ion cluster for the -glucose of the glucose aldonitrile-pentaacetate derivative was monitored from m/z  <dig> to m/z  <dig>
spectral data were corrected using regression analysis to extract natural 13c enrichment from results  <cit> .
measurement of 13c label distribution determined the different relative distribution percentages of the mass isotopomers, m <dig> , m <dig> , m <dig> , etc., which were reported as molar fractions.
<dig>  analyses and modelling of data using isodyn
transfer of 13c from -glucose medium into intracellular metabolites was simulated by isodyn, a software program written in c++ and designed for analyses of data from stable isotopic tracers  <cit> , <cit> , <cit> .
latest modifications of this software have been described in detail  <cit> , <cit> .
this software automatically constructs and solves the large system of ordinary differential equations describing the evolution of isotopomer concentrations of metabolites produced in glycolysis, tca cycle, and ppp .
all the reactions of the non-oxidative branch catalyzed by transketolase and transaldolase are accounted for, assuming pentose-5-phosphates to be a single pool.
the reactions of the oxidative branch are passed to one reaction leading from g6p to pentose-phosphates.
at this initial moment, all metabolites  are assumed to be non-labeled, and the latter is set in accordance with the known isotopomer composition of the used labeled substrate.
initial total concentrations of intracellular metabolites are calculated as functions of model parameters assuming a steady state at the initial moment.
functions designed especially for each type of reaction  simulate transformation of the carbon skeleton , and the production and consumption rate for each isotopomer of the considered system.
these transformations redistribute 13c isotopes in all metabolites.
individual rates are computed for each isotopomer, and these rates determine the values of the derivatives for the isotopomers.
this system is solved using a method of numerical integration that can be chosen arbitrarily .
starting from the initial values corresponding to the experimental conditions of incubation, isodyn simulates a real-time course of label propagation.
experimental and computed data for corresponding time points are compared.
in this way, reaching an isotopic steady state is not required.
<dig>  rna extraction and real-time rt-qpcr gene-expression analyses
total rna was isolated from l6e <dig> cells using an rneasy mini kit and a qiacube automat .
rna integrity was checked by capillary electrophoresis .
concentrations were checked with ribogreen  and a victor-x <dig> multilabel reader .
labelled cdna was synthesized from 600 ng of total rna using superscript-ii, reverse transcriptase and random hexamers .
oligonucleotide sequences for the reverse transcriptase reaction, synthesized by eurogentec , were: pyruvate carboxylase , pc, rpc-f, ttccgtgtccgaggtgtaaag, and rpc-r, cgctaggaactgctggttgtt; pyruvate kinase muscle isozyme, pkm <dig>  rpkm-784f, gagcaggacgtggacatggt, and rpkm-864r, ctctcccaggaccttcctaacc; pyruvate dehydrogenase e <dig> alpha  <dig>  pdha <dig>  rnm_ <dig> f, tttgtcttctgtgctgggagactg and rnm_ <dig> r, gtagatgggtggcttcaagtttgc; pyruvate dehydrogenase kinase isozyme  <dig>  pdk <dig>  rnm_ <dig> f, catggctaagctcctgtgtgac, and rnm_ <dig> r, ggacgtagaccatgtgaatggg; pyruvate dehydrogenase kinase isozyme  <dig>  pdk <dig>  rnm_ <dig> f, tgctcatgaaccagcacatcctc and rnm_ <dig> r, tcccaatgtggcttgggtttcc; cytosolic nadp-isocitrate dehydrogenase, idh <dig>  rnm_ <dig> f, agacgtccaccaatcccattgc, and rnm_ <dig> r, tcaagctttgctctgtgggctaac; nadp-dependent malic enzyme, me <dig>  rnm_ <dig> f, tgcaagaccatggttcccagac, and rnm_ <dig> r, aagctaagcccagggacattaggg; cytosolic phosphoenolpyruvate carboxykinase   <dig>  pck <dig>  rpepck1-2228f, tggttccacttcgaggtcact, and rpepck1-2310r, caagtatgttttctgtgcactttagct.
real-time pcr  was conducted with a sybrgreen master mix reagent using sureprime core kit  in an icycler  apparatus.
the standard curve method was used for relative quantification of pcr products, and gene expressions normalized to 18s rna quantification, which has been found to be a reliable internal control gene in our hands and by other scholars  <cit> - <cit> .
<dig>  statistical analyses
experiments relating to glucose transport were carried out using three cultures each time for each treatment regimen and then repeated thrice.
metabolic characterization assays were carried out in triplicate and repeated twice.
rnas for rt-qpcr assays were isolated from three independent extracts.
statistical analyses of glucose transport as well as biochemical and metabolite mass isotopomer distributions were calculated by two-way anova between groups.
data for real-time pcr were analyzed by the one-way anova between groups.
confidence intervals at a confidence level of 95% or 99%  were taken as indicating significant differences in each parameter analyzed.
differences from the mean were considered significant for p <  <dig> .
in the case of the metabolic fluxes estimated using isodyn, the difference between calculated and experimental data  was minimized using a simulated annealing algorithm.
the goodness of best fit  was checked based on the value of χ <dig> and number of degrees of freedom  <cit> .
in the case that, according to this checking, the fit was deemed acceptable, the metabolic fluxes corresponding to the best fit were accepted as consistent with the measured isotopomer distribution.
the confidence intervals for the fluxes were determined as the upper and lower limits of fluxes that produced χ <dig> that was less than the fixed threshold  <cit> .
abbreviations
2-dg: 2-deoxy-d--glucose
g6p: glucose-6-phosphate
gp: glycogen phosphorylase
gs: glycogen synthase
ir: insulin resistance
pc: pyruvate carboxylase
pdc: pyruvate dehydrogenase complex
pep: phosphoenolpyruvate
pepck: phosphoenolpyruvate carboxykinase
ppp: pentose phosphate pathway
t2dm: type- <dig> diabetes mellitus
tca: tricarboxylic acid
competing interests
the authors declare that they have no competing interests.
authors’ contributions
sg, sm, fs, az, pr and mc conceived and designed the research.
sg, sm, am, rh, fs, at and yd undertook the experiments.
sg, sm, vs, am, rh, fs, at, yd and mc analysed data.
sg, sm, am, vs, jc, az, pr, and mc interpreted the results of the experiments.
sg, sm, jc, am, rh, fs, at, yd and pr prepared the figures.
sg, sm, am, vs, pr and mc drafted the manuscript.
sg, sm, am, vs, jc, az, pr and mc edited and revised the manuscript.
all authors read and approved the final manuscript.

 obesity is increasing to epidemic proportions,, starting from young ages,, and is associated with an increase in the prevalence of type- diabetes mellitus throughout the world,. in recent years, active participation of immune cells in obesity and t2dm has become evident,. this chronic activation of the innate immune system can lead to insulin resistance, impaired tolerance to glucose and, eventually, t2dm,. adipose tissue produces a vast array of adipocyte-derived factors that regulate the metabolism, inflammation and body mass,. one of these adipokines, resistin, belongs to a family of cysteine-rich proteins shown to be involved in inflammation and altered insulin sensitivity in rodents,. resistin is produced from adipose tissue in rodents whereas, in humans, it is secreted by the mononuclear lymphocytes and stromal cells within adipose tissue,. studies in humanised resistin mice suggest that, even though the site of resistin production differs between species, human resistin exacerbates inflammation in white adipose tissue and contributes to ir, thereby impairing its normal effects,. ir is defined as the genetic or learned inability of target tissues to respond normally to the action of circulating hormones. it has been described in skeletal muscle but also in liver and adipose tissues. some studies have shown that resistin affects glucose transport and insulin-stimulated oxidation of glucose in l skeletal muscle cells -. it also decreases the uptake and oxidation of long-chain fatty acids and glycogen synthase kinase-3-β, as well as insulin-stimulated insulin receptor substrate- tyrosine phosphorylation, in the same cell line. furthermore, high levels of resistin in rats leads to ir involving impaired insulin signalling in skeletal muscle, liver and adipose tissues, resulting in glucose intolerance, hyperinsulinemia and hypertriglyceridemia. these observations link resistin in muscle-cell metabolism to ir. changes in the metabolome have been described as the “ultimate” response of an organism to various events, such as genetic alterations as well as disease-based or environmental influences. however, the metabolome is dynamic and metabolites are transformed continuously in the cell. comprehensive characterisation of the metabolic networks and their functional states requires quantitative knowledge of intracellular metabolic fluxes. these intracellular fluxes can be quantified by analysing incorporation of labelled substrates in metabolic products using appropriate bioinformatic tools. skeletal muscle is considered to be the main tissue involved in the maintenance of glucose homeostasis because its contribution to glucose uptake is ≈ 75% of the total contribution of peripheral tissues, and it is the main tissue responsible for insulin-dependent glucose use. also, pyruvate homeostasis-related fluxes in skeletal muscle play a critical part in glucose homeostasis -. a better understanding of how these fluxes work together in muscle cells is needed to understand this process. the aim of the present study was to examine the effect of resistin on the central carbon metabolic network of rat skeletal muscle cells, and to ascertain how resistin alters the response of muscle cells to insulin at the metabolic level. we used -glucose as a tracer, analysed isotopomer distributions by gas cromatography coupled with mass spectrometry, and quantified metabolic fluxes using isodyn software, -. in the model of the glucose metabolic network in skeletal muscle we included the glycolytic pathway and pentose phosphate pathway. pyruvate dehydrogenase complex -, phosphoenolpyruvate carboxykinase - and pyruvate carboxylase -catalysed reactions were also included to account for pyruvate homeostasis-related fluxes. resistin-induced alterations in the profile of metabolic fluxes of rat skeletal muscle cells were identified, and were complemented by analyses of the expression of certain genes. the study presented here offers a wider vision than previously published of the metabolic reorganisation that resistin induces in rat skeletal muscle cells and in the metabolic response of these cells to insulin. our observations provide new data for a more precise understanding of the effect of resistin on metabolic-network adaptations of skeletal muscle cells, and help to explain the role of resistin in the ir described in myocytes. our work reveals the importance of tackling the study of complex biological systems from a systems-biology point of view to fully understand multifactorial diseases and identify new therapeutic targets.  for analyses of the effect of resistin on glucose metabolism and on the metabolic response of muscle cells to insulin, l6e myotubes were pretreated for 8 h with or without resistin. after preincubation, myotubes were incubated for a further 6 h with glucose that was 50% enriched in -glucose in the absence or presence of nm insulin and in the continued presence or absence of resistin. after 6 h of incubation, glucose consumption, lactate production, intracellular levels of glycogen and glucose-6-phosphate were determined, as were the mass-isotopomer distributions in glucose and lactate from the incubation medium, as well as glycogen glucose and ribose isolated from rna in cell pellets. analyses of the obtained data by isodyn determined the distribution of the metabolic flux profile of l6e myotubes under various incubation conditions. figure  summarises the qualitative changes in metabolites and fluxes. l6e myotubes were treated or not treated with nm resistin for 8 h and then incubated for 6 h in the absence or presence of nm insulin. glucose consumption, lactate production, and intracellular contents of glycogen and glucose-6-phosphate were determined as described in the methods section. results are the mean ± standard deviation. units of glucose consumption and lactate production refer to the volume of the incubation medium. indicates p <  and p <  compared with the condition of neither resistin nor insulin. indicates p <  compared with the l6e myotubes treated only with resistin. lactate         glycogen           ribose           l6e myotubes were treated or not treated with nm resistin for 8 h and then incubated for 6 h with 10 mm glucose, 50%-enriched in -glucose in the absence or presence of nm insulin. mass isotopomer distributions at the end of incubations were determined, as described in the methods section. results are the mean ± standard deviation. indicates p <  compared with the condition of neither resistin nor insulin. glycolytic and ppp fluxes glycogen-related fluxes pyruvate homeostasis-related fluxes tricarboxylic acid cycle fluxes principal fluxes of central carbon metabolism calculated in l6e myotubes treated or not treated with nm resistin for 8 h and then incubated for 6 h with 10 mm glucose, 50%-enriched in -glucose in the absence or presence of nm insulin. glucose phosphorylation flux was fixed according to the glucose consumption measured. remaining fluxes are expressed by median values as well as minimum and maximum values from the best flux sets. fluxes are expressed in nmol · ml–1 · min– where volume is the volume of the incubation medium. χ values for the averages were: ins– res–:; ins– res+:; ins + res–:; ins + res+:. indicates that the flux is different to the control condition with a 99% confidence interval, indicates that the flux is different to the control condition with a 95% ci, and indicates that the flux differs from the condition in which cells were treated only with resistin with a 99% ci. ppp, pentose phosphate pathway; oaa, oxaloacetate. basal metabolic flux profile of l6e9 isodyn quantified the metabolic-flux distribution in l6e myotubes incubated in the presence of glucose and in the absence of insulin or resistin. under this incubation condition, cells had a glucose phosphorylation of  nmol · ml–1 · min– this flux was maintained throughout the subsequent steps of “upper glycolysis” until the split of fructose- -bisphosphate into two triose-phosphates by aldolase. from this reaction onwards, the glycolytic flux increased twofold. additionally, the flux across the oxidative branch of the ppp was around eight-times lower than the upper glycolytic flux. with regard to glycogen metabolism, glycogen phosphorylase and glycogen synthase were active, and glycogen recycling was estimated to be ≈ 2% of the total glucose phosphorylation flux. however, the most important contribution of isodyn for interpretation of experimental data was the possibility of flux quantification around the final steps of lower glycolysis and the tricarboxylic acid cycle. in that regard, isodyn revealed that ≈ 49% of phosphoenolpyruvate and pyruvate entered the tca cycle through the reactions catalysed by pdc or pc, with a relative pdc:pc ratio of 8: around 44% of the pep–pyruvate flux was diverted to the synthesis and output of lactate towards the cell culture medium, and pyruvate cycling flux. this finding accounted for the fact that feedback of molecules from the tca cycle back towards the pep–pyruvate pool was around ten-times higher than the flux of the tca cycle towards citrate synthesis. effect of resistin on the metabolic flux profile of l6e9 l6e myotubes treated with nm resistin showed 13% higher flux of glucose phosphorylation. fitting of the experimental data from this incubation condition using isodyn suggested that resistin increased the glycolytic flux, with a significant increase in the glucose phosphorylation flux, and a proportional increase of glycolytic fluxes. resistin induced changes in the non-oxidative branch of the ppp, but the relative flux through the oxidative branch was unchanged. glycogen recycling-related fluxes remained unchanged, as indicated by the fact that the net 13c incorporation in glycogen was very similar to the control condition. moreover, resistin increased lactate production, which translated to an increase in the lactate dehydrogenase flux, though it was not significant. the most significant metabolic change after isodyn fitting was observed in the pdc-related flux, which was increased significantly. if we consider the relative ratios of the two main pyruvate incorporation fluxes, there was an increase of the pdc prevalence to 20: this finding was also supported by the level of pyruvate dehydrogenase kinase isozyme mrna, and was observed as a trend in the level of pyruvate dehydrogenase pdha mrna. despite the large change in this ratio, the effect on fluxes in the tca cycle was almost negligible even though there was an increase in the acetyl-coa output flux. the pyruvate cycling flux was unchanged by resistin treatment alone, which was supported by the level of pepck mrna. effect of resistin on the metabolic response of insulin in l6e myotubes having investigated the alterations caused by resistin on the metabolic profile of l6e myotubes, we examined its effect on the metabolic response of myotubes to insulin. first we ascertained if resistin affected stimulation of insulin-induced glucose transport. thus, induction of 2-deoxy-d--glucose transport was assayed 30 min after incubation with nm insulin in l6e myotubes treated or not treated with nm resistin. insulin induced an increase in 2-dg transport in l6e myotubes of  ±  %. however, it induced an increase in glucose transport in cells that had been pre-treated with nm resistin for 8 h of  ±  %. subsequently, we analysed changes in the metabolic flux profile induced by insulin in l6e myotubes treated or not treated with nm resistin. in the absence of resistin, insulin stimulated an increase in glucose consumption and lactate production in l6e myotubes. there were no changes in glycogen content, but a 150% increase in intracellular g6p content. fitting of experimental data using isodyn revealed that the flux of lactate production was 48% higher than that in cells not treated with insulin. insulin also increased the flux of glucose phosphorylation and glycolytic fluxes by ≈ 59%, and induced a 25% increase in the oxidative ppp flux. moreover, fluxes related to glycogen metabolism did not change significantly after insulin treatment. insulin increased the flux through pdc by around -fold, but did not affect the remainder of tca cycle-related fluxes significantly, as confirmed by the reduction in levels of pdk mrna and the increasing tendency of levels of pdha mrna. the relative ratio of pdc:pc fluxes increased from 8: to 13: after insulin treatment. insulin also led to a 59% increase in the acetyl-coa output fluxes. in the presence of resistin, insulin increased the flux of glucose phosphorylation by 16% more than when cells were treated with resistin alone. glycolysis fluxes did not increase significantly compared with the resistin condition, but the lactate dehydrogenase flux increased by 45%. the increase induced by resistin treatment in the oxidative ppp pathway was not significantly altered if cells were also treated with insulin. similarly, the gs and gp fluxes did not change significantly after resistin treatment. interestingly, there was no change in the pdc flux or the acetyl-coa production flux, and no change was observed in the remainder of tca cycle-related fluxes. if the ratio between pdc and pc fluxes is taken into account, it shifted from 20: towards 70: after resistin was combined with insulin. the final metabolic state of insulin-stimulated l6e did not match the one of cells treated with resistin. some fluxes were equal in resistin-treated and -untreated cells: oxidative ppp flux, gs flux, gp flux, and tca cycle-related fluxes. most of the quantified fluxes were lower in resistin, but significant changes were found only in the lactate dehydrogenase flux, and the pdc:pc ratio. effect of resistin on gene expression of the enzymes involved in the pdc, pc and pyruvate cycling fluxes in l6e9 to investigate further if the observed effects in metabolic-flux distribution were prompted by a change in the gene expression of any of the enzymes involved in their activities, we analysed the gene expression of several enzymes involved in pyruvate homeostasis at the end of the four incubation conditions tested. neither insulin, nor resistin, nor combination of the two hormones had a significant effect on levels of pyruvate carboxylase mrna or pyruvate kinase muscle isozyme mrna. however, insulin increased levels of pyruvate dehydrogenase e and decreased pyruvate dehydrogenase kinase isozyme mrna. presence of resistin did not alter the decrease in levels of pdk mrna caused by insulin, but prevented the increase in levels of pdha mrna. mrna expression of pyruvate dehydrogenase kinase isozyme decreased after treatment with both hormones, with an additive effect if both were combined. these results suggest that resistin prevents the insulin-induced increase in pdha expression, and does not prevent the increase in pdc activity caused by a reduction of expression of pdk and pdk genes. moreover, mrna levels of cytosolic nadp-isocitrate dehydrogenase and nadp-dependent malic enzyme, both causes of pyruvate cycling increasing nadph production, were not affected by any of the hormones tested under these conditions. however, insulin significantly decreased levels of pck mrna in the absence of resistin, but not in the presence of resistin. these results suggest that resistin prevents inhibition of expression of the pck gene. hence, cytosolic flux in pepck may be active in the presence of resistin independent of the presence of insulin. discussion after the discovery of several adipocyte-secreted molecules, adipocytes were recognized as part of endocrine tissue,,. these molecules enable adipocytes to communicate with other tissues and organs, and to regulate: the metabolism of lipids and glucose; energy balance; insulin action; cell proliferation. resistin was identified as an adipokine, and it was suggested that it might link obesity and ir,. however, the role of resistin in the pathophysiology of ir in humans and animals, and how it acts in muscle, liver and fat, is controversial,. here, we investigated the effect of resistin on basal and insulin-stimulated glucose metabolism in l6e rat myotubes. we used tracer-based metabolomics and our in-house isodyn software to analyse quantitatively metabolic flux distribution in this cell line under different incubation conditions. the response to insulin and other hormone stimuli in l6e cells is similar to that observed in skeletal muscle in vivo. hence, l6e cells are considered to be suitable models for analyses of the effect of resistin on glucose metabolism and on the metabolic response of muscle cells to insulin. analyses with isodyn data provided an overview of the central carbon metabolism of l6e skeletal muscle cells. analyses were undertaken assuming a metabolic steady state, and rapid mixing of isotopic isomers of the same metabolite throughout different intracellular compartments. incubation with a large excess of glucose ensured the required glucose uptake throughout the incubation. the other conditions of incubation did not change, so the assumption of a metabolic steady state seemed reasonable. rapid mixing of isotopic isomers of the same species throughout the cell is accepted for analyses of data of stable isotopes. fitting of experimental data in the framework of these assumptions suggested a set of metabolic fluxes underlying the measured distributions of isotopic isomers. the effect of insulin on l6e myotubes increases glycolysis, ppp, and incorporation of pyruvate to the tca cycle through pdc and fatty-acid synthesis our results suggested that l6e myotubes were highly glycolytic, with active synthesis and degradation of glycogen, and showed tca-cycle activity and acetyl-coa output. oxidative ppp was weakly active in this cell type. furthermore, use of the same tools revealed that insulin activated the uptake and phosphorylation of glucose, glycolysis and lactate production, the oxidative branch of the ppp, pyruvate incorporation towards the tca cycle and increased acetyl-coa output. we observed that insulin altered the ratio between pdc and pc fluxes by increasing the pdc flux through an increase of levels of pdha and decrease of pdk and pdk mrna. the effects of insulin on the metabolic flux distribution of l6e myotubes reported in the present study correlate with the known activation of pdc by profound suppression of pdk expression in skeletal muscle, inhibition of expression of the pck gene, and the increase in glucose uptake and lactate production observed using other approaches,. the effect of resistin on l6e myotubes increases glycolysis and fatty-acid synthesis and alters the ratio between pdc and pc, tilting it towards direct incorporation of pyruvate through pdc our results revealed that, in the absence of insulin, resistin activated glucose phosphorylation and the glycolytic pathway. also, resistin did not affect the glycogen content in skeletal muscle. even though enrichment of the glycogen label was not significantly different to the control condition, the tendency was for glycogen-related fluxes to be lower, a result that was in accordance with other works. furthermore, our approach showed that, even though many of the tca-cycle fluxes remained unchanged, the ratio between the pdc and pc fluxes was significantly altered, effectively increasing the incorporation of pyruvate through pdc into mitochondria. pc activity is finely regulated by mechanisms such as allosteric activation by acetyl carnitine, which was not tested in this study. the acetyl-coa output for fatty-acid synthesis was also increased after resistin treatment, which would explain the higher flux towards the tca cycle. the effect of resistin impairs normal insulin metabolic actions we showed that resistin affects the normal metabolic response of l6e myotubes to insulin. thus, the increase in glucose uptake and glycolytic flux induced by insulin was lower than that observed in resistin-treated myotubes. these results are in accordance with those published by other authors for l muscle cells,,. resistin treatment, however, did not affect the insulin-stimulated content and recycling of glycogen. moreover, in resistin-treated cells, insulin did not increase the pdc flux. after resistin treatment, levels of pdk mrna and pdk mrna were lower, accompanied with a reduction in levels of pdha mrna. the reduction in expression of the main enzyme and key regulators of its activity resulted in a null net effect on the final flux. the ratio between the pdc and pc fluxes was also altered, showing reduction in the pc anaplerotic entrance into the cycle. fluxes through the tca cycle and acetyl-coa output fluxes were not increased, but instead maintained or reduced, by insulin in resistin-treated cells. this lack of metabolic activation in response to insulin can be explained by the slight decrease observed in glucose uptake and the increase of lactate dehydrogenase flux directed to lactate production, which impaired the increase in the levels of glycolytic intermediates that could be metabolized through other metabolic pathways. our results provide mechanistic explanations for the observations of palanivel and co-workers, who reported that the insulin-stimulated oxidation of glucose via the tca cycle was reduced by resistin, and that this reduction would compromise oxidative respiration in the cell.  several studies have identified positive correlations between resistin levels and ir in vivo and in vitro -. overall, our results showed that resistin significantly alters glucose metabolism in l6e myotubes and their response to insulin. of the whole set of metabolic changes induced by resistin in l6e myotubes, the imbalance in pdc and pc fluxes as well as the higher acetyl-coa output flux for fatty-acid synthesis and other metabolic processes should be emphasized. resistin impairs the normal metabolism of insulin, and the metabolism is directed to try to maintain the imbalance between the pdc and pc fluxes. in this way, l6e reduces tca-cycle fluxes and acetyl-coa output fluxes, which leads to an important metabolic imbalance in the central carbon metabolic pathways of l6e myotubes. it has been reported that ir in skeletal muscle, liver and adipose tissue is usually accompanied by glucose intolerance and hypertriglyceridemia. in this way, our results showed the importance of the balance of pdc and pc in glucose metabolism. our results support the idea that metabolic disorders such as t2dm are caused by complex multi-molecular interactions that cannot be explained readily by an alteration in expression of a single gene or gene product, or even alteration of a single enzymatic cascade. the robustness of results obtained on the effects of insulin on l6e metabolism found by combining tracer-based metabolomic data and flux analyses using isodyn showed that this approach is a suitable tool to study the effects of hormones on the central carbon metabolism network in myocytes. methods cells and reagents myoblastic cell line l6e was kindly provided by dr. nadal-ginard. rat recombinant resistin was produced in a hek cell line. the purity of resistin was controlled by sodium dodecyl sulfate–polyacrylamide gel electrophoresis and estimated to be >99%. biological activity was validated as described. resistin was tested for endotoxin using a toxinsensor chromogenic lal endotoxin assay kit. this assay indicated low levels of endotoxins at the final concentration used, which is below current usa food and drug administion limits that require eluates from medical devices to be < eu/ml. -d-glucose and -d-glucose were purchased from isotec and 2-deoxy-d--glucose from ge healthcare. conditions for cell culture and incubation myoblastic cell line l6e was grown as described. for biochemical and mass isotopomer assays, l6e myoblasts were seeded at cells/cm in 60-mm diameter dishes with 4 ml dulbecco’s modified eagle medium without glucose or glutamine supplemented with 10% heat-inactivated foetal bovine serum, 10 mm glucose, 4 mm asparagine and 5% penicillin/streptomycin. cell differentiation was induced by lowering fbs to a final concentration of 2% when myoblasts were pre-confluent. cells at day- of differentiation were washed four times in phosphate-buffered saline and pre-treated or not with nm resistin dissolved in dmem without fbs and supplemented with 10 mm glucose, 4 mm asparagine, 2% hepes and 1% penicillin/streptomycin. after 8 h, the medium was replaced by another with the same composition, but now with 10 mm glucose 50%-enriched in -glucose in the presence or absence of nm insulin. at the end of 6-h incubation, media were collected and dishes washed twice with pbs before freezing in liquid nitrogen. incubation media and dishes were stored at –80°c until biochemical and mass isotopomer analyses. glucose transport assay to carry out the glucose transport assay, l6e cells/dish were cultured on six-well plates. at day- of differentiation, cells were pre-treated or not with nm resistin for 8 h, as described above. next, uptake of 100 μci 2-dg by myotubes incubated with or without nm insulin for 30 min was measured as described. biochemical determinations of glucose, lactate and g6p concentrations of glucose and lactate in incubation media at the beginning and end of 6-h incubations were measured as described previously by spectrophotometric methods coupled to a cobas mira plus chemical analyzer. intracellular g6p content was determined by enzymatic means as described. mass isotopomer analyses of lactate, medium and glycogen glucose, and rna ribose, and determination of glycogen content by gc/ms mass spectral data were obtained on a qp shimadzu mass selective detector connected to a gc- gas chromatograph using helium as the carrier gas and isobutane  pa as reagent gas in chemical ionization analyses. the settings were: gc inlet, 250°c for glucose, ribose and deoxyribose, and 200°c for lactate; transfer line, 250°c; ms source, 200°c. a varian vf- capillary column was used to analyze all compounds. glucose from the frozen media was purified using a tandem set of dowex-1x8/dowex-50wx ion-exchange columns. the glucose from media was converted to its glucose aldonitrile pentaacetate derivative as described. ion clusters around m/z were monitored. lactate from the media was extracted with ethyl acetate after acidification with hydrochloric acid and derivatized to its propylamide-heptafluorobutyric form. clusters around m/z were monitored. rna ribose was isolated by acid hydrolysis of cellular rna after trizol purification of cell extracts. ribose isolated from rna was derivatized to its aldonitrile-acetate form using hydroxylamine in pyridine and acetic anhydride. we monitored ion clusters around m/z to find the molar enrichment and positional distribution of 13c labels in ribose as described. glycogen in frozen cell monolayers from different analytical conditions was extracted and quantified as described. extraction was done by direct digestion of sonicated extracts with amyloglucosidase. mass isotopomer analyses of glycogen glucose was undertaken as described above for medium glucose. measurement of glycogen content was carried out using the isotopomer -glucose as recovery standard and internal standard quantification procedures. the ion cluster for the -glucose of the glucose aldonitrile-pentaacetate derivative was monitored from m/z to m/z spectral data were corrected using regression analysis to extract natural 13c enrichment from results. measurement of 13c label distribution determined the different relative distribution percentages of the mass isotopomers, m, m, m, etc., which were reported as molar fractions. analyses and modelling of data using isodyn transfer of 13c from -glucose medium into intracellular metabolites was simulated by isodyn, a software program written in c++ and designed for analyses of data from stable isotopic tracers,,. latest modifications of this software have been described in detail,. this software automatically constructs and solves the large system of ordinary differential equations describing the evolution of isotopomer concentrations of metabolites produced in glycolysis, tca cycle, and ppp. all the reactions of the non-oxidative branch catalyzed by transketolase and transaldolase are accounted for, assuming pentose-5-phosphates to be a single pool. the reactions of the oxidative branch are passed to one reaction leading from g6p to pentose-phosphates. at this initial moment, all metabolites are assumed to be non-labeled, and the latter is set in accordance with the known isotopomer composition of the used labeled substrate. initial total concentrations of intracellular metabolites are calculated as functions of model parameters assuming a steady state at the initial moment. functions designed especially for each type of reaction simulate transformation of the carbon skeleton, and the production and consumption rate for each isotopomer of the considered system. these transformations redistribute 13c isotopes in all metabolites. individual rates are computed for each isotopomer, and these rates determine the values of the derivatives for the isotopomers. this system is solved using a method of numerical integration that can be chosen arbitrarily. starting from the initial values corresponding to the experimental conditions of incubation, isodyn simulates a real-time course of label propagation. experimental and computed data for corresponding time points are compared. in this way, reaching an isotopic steady state is not required. rna extraction and real-time rt-qpcr gene-expression analyses total rna was isolated from l6e cells using an rneasy mini kit and a qiacube automat. rna integrity was checked by capillary electrophoresis. concentrations were checked with ribogreen and a victor-x multilabel reader. labelled cdna was synthesized from 600 ng of total rna using superscript-ii, reverse transcriptase and random hexamers. oligonucleotide sequences for the reverse transcriptase reaction, synthesized by eurogentec, were: pyruvate carboxylase, pc, rpc-f, ttccgtgtccgaggtgtaaag, and rpc-r, cgctaggaactgctggttgtt; pyruvate kinase muscle isozyme, pkm rpkm-784f, gagcaggacgtggacatggt, and rpkm-864r, ctctcccaggaccttcctaacc; pyruvate dehydrogenase e alpha pdha rnm_ f, tttgtcttctgtgctgggagactg and rnm_ r, gtagatgggtggcttcaagtttgc; pyruvate dehydrogenase kinase isozyme pdk rnm_ f, catggctaagctcctgtgtgac, and rnm_ r, ggacgtagaccatgtgaatggg; pyruvate dehydrogenase kinase isozyme pdk rnm_ f, tgctcatgaaccagcacatcctc and rnm_ r, tcccaatgtggcttgggtttcc; cytosolic nadp-isocitrate dehydrogenase, idh rnm_ f, agacgtccaccaatcccattgc, and rnm_ r, tcaagctttgctctgtgggctaac; nadp-dependent malic enzyme, me rnm_ f, tgcaagaccatggttcccagac, and rnm_ r, aagctaagcccagggacattaggg; cytosolic phosphoenolpyruvate carboxykinase pck rpepck1-2228f, tggttccacttcgaggtcact, and rpepck1-2310r, caagtatgttttctgtgcactttagct. real-time pcr was conducted with a sybrgreen master mix reagent using sureprime core kit in an icycler apparatus. the standard curve method was used for relative quantification of pcr products, and gene expressions normalized to 18s rna quantification, which has been found to be a reliable internal control gene in our hands and by other scholars -. statistical analyses experiments relating to glucose transport were carried out using three cultures each time for each treatment regimen and then repeated thrice. metabolic characterization assays were carried out in triplicate and repeated twice. rnas for rt-qpcr assays were isolated from three independent extracts. statistical analyses of glucose transport as well as biochemical and metabolite mass isotopomer distributions were calculated by two-way anova between groups. data for real-time pcr were analyzed by the one-way anova between groups. confidence intervals at a confidence level of 95% or 99% were taken as indicating significant differences in each parameter analyzed. differences from the mean were considered significant for p < . in the case of the metabolic fluxes estimated using isodyn, the difference between calculated and experimental data was minimized using a simulated annealing algorithm. the goodness of best fit was checked based on the value of χ and number of degrees of freedom. in the case that, according to this checking, the fit was deemed acceptable, the metabolic fluxes corresponding to the best fit were accepted as consistent with the measured isotopomer distribution. the confidence intervals for the fluxes were determined as the upper and lower limits of fluxes that produced χ that was less than the fixed threshold. abbreviations 2-dg: 2-deoxy-d--glucose g6p: glucose-6-phosphate gp: glycogen phosphorylase gs: glycogen synthase ir: insulin resistance pc: pyruvate carboxylase pdc: pyruvate dehydrogenase complex pep: phosphoenolpyruvate pepck: phosphoenolpyruvate carboxykinase ppp: pentose phosphate pathway t2dm: type- diabetes mellitus tca: tricarboxylic acid competing interests the authors declare that they have no competing interests. authors’ contributions sg, sm, fs, az, pr and mc conceived and designed the research. sg, sm, am, rh, fs, at and yd undertook the experiments. sg, sm, vs, am, rh, fs, at, yd and mc analysed data. sg, sm, am, vs, jc, az, pr, and mc interpreted the results of the experiments. sg, sm, jc, am, rh, fs, at, yd and pr prepared the figures. sg, sm, am, vs, pr and mc drafted the manuscript. sg, sm, am, vs, jc, az, pr and mc edited and revised the manuscript. all authors read and approved the final manuscript. 
these data suggest that resistin impairs the metabolic activation of insulin. resistin significantly increased glucose uptake and glycolysis, altering pyruvate utilisation by the cell. in the presence of resistin, insulin only slightly increased glucose uptake and glycolysis, and did not alter the flux profile around pyruvate induced by resistin. it has been suggested that the adipokine resistin links obesity and insulin resistance, although how resistin acts on muscle metabolism is controversial. we aimed to quantitatively analyse the effects of resistin on the glucose metabolic flux profile and on insulin response in l6e myotubes at the metabolic level using a tracer-based metabolomic approach and our in-house developed software, isodyn. resistin prevented the increase in gene expression in pyruvate dehydrogenase-e and the sharp decrease in gene expression in cytosolic phosphoenolpyruvate carboxykinase- induced by insulin.
 it has been suggested that the adipokine resistin links obesity and insulin resistance, although how resistin acts on muscle metabolism is controversial. we aimed to quantitatively analyse the effects of resistin on the glucose metabolic flux profile and on insulin response in l6e myotubes at the metabolic level using a tracer-based metabolomic approach and our in-house developed software, isodyn.  resistin significantly increased glucose uptake and glycolysis, altering pyruvate utilisation by the cell. in the presence of resistin, insulin only slightly increased glucose uptake and glycolysis, and did not alter the flux profile around pyruvate induced by resistin. resistin prevented the increase in gene expression in pyruvate dehydrogenase-e and the sharp decrease in gene expression in cytosolic phosphoenolpyruvate carboxykinase- induced by insulin.  these data suggest that resistin impairs the metabolic activation of insulin. this impairment cannot be explained by the activity of a single enzyme, but instead due to reorganisation of the whole metabolic flux distribution. resistininsulin resistanceglucose metabolismfluxomicsbioinformaticstracer-based metabolomics 