Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

Mitochondrial DNA damage and animal longevity: insights from comparative studies

Chemical reactions in living cells are under strict enzyme control and conform to a tightly regulated metabolic program. However, uncontrolled and potentially deleterious endogenous reactions occur, even under physiological conditions. Aging, in this chemical context, could be viewed as an entropic process, the result of chemical side reactions that chronically and cumulatively degrade the function of biological systems. Mitochondria are a main source of reactive oxygen species (ROS) and chemical sidereactions in healthy aerobic tissues and are the only known extranuclear cellular organelles in animal cells that contain their own DNA (mtDNA). ROS can modify mtDNA directly at the sugar-phosphate backbone or at the bases, producing many different oxidatively modified purines and pyrimidines, as well as single and double strand breaks and DNA mutations. In this scenario, natural selection tends to decrease the mitochondrial ROS generation, the oxidative damage to mtDNA, and the mitochondrial mutation rate in long-lived species, in agreement with the mitochondrial oxidative stress theory of aging.

Identifying the preferred subset of enzymatic profiles in nonlinear kinetic metabolic models via multiobjective global optimization and pareto filters

Optimization models in metabolic engineering and systems biology focus typically on optimizing a unique criterion, usually the synthesis rate of a metabolite of interest or the rate of growth. Connectivity and non-linear regulatory effects, however, make it necessary to consider multiple objectives in order to identify useful strategies that balance out different metabolic issues. This is a fundamental aspect, as optimization of maximum yield in a given condition may involve unrealistic values in other key processes. Due to the difficulties associated with detailed non-linear models, analysis using stoichiometric descriptions and linear optimization methods have become rather popular in systems biology. However, despite being useful, these approaches fail in capturing the intrinsic nonlinear nature of the underlying metabolic systems and the regulatory signals involved. Targeting more complex biological systems requires the application of global optimization methods to non-linear representations. In this work we address the multi-objective global optimization of metabolic networks that are described by a special class of models based on the power-law formalism: the generalized mass action (GMA) representation. Our goal is to develop global optimization methods capable of efficiently dealing with several biological criteria simultaneously. In order to overcome the numerical difficulties of dealing with multiple criteria in the optimization, we propose a heuristic approach based on the epsilon constraint method that reduces the computational burden of generating a set of Pareto optimal alternatives, each achieving a unique combination of objectives values. To facilitate the post-optimal analysis of these solutions and narrow down their number prior to being tested in the laboratory, we explore the use of Pareto filters that identify the preferred subset of enzymatic profiles. We demonstrate the usefulness of our approach by means of a case study that optimizes the ethanol production in the fermentation of Saccharomyces cerevisiae.

Effects of a post-weaning cafeteria diet in young rats: metabolic syndrome, reduced activity and low anxiety-like behaviour.

Among adolescents, overweight, obesity and metabolic syndrome are rapidly increasing in recent years as a consequence of unhealthy palatable diets. Animal models of diet-induced obesity have been developed, but little is known about the behavioural patterns produced by the consumption of such diets. The aim of the present study was to determine the behavioural and biochemical effects of a cafeteria diet fed to juvenile male and female rats, as well as to evaluate the possible recovery from these effects by administering standard feeding during the last week of the study. Two groups of male and female rats were fed with either a standard chow diet (ST) or a cafeteria (CAF) diet from weaning and for 8 weeks. A third group of males (CAF withdrawal) was fed with the CAF diet for 7 weeks and the ST in the 8th week. Both males and females developed metabolic syndrome as a consequence of the CAF feeding, showing overweight, higher adiposity and liver weight, increased plasma levels of glucose, insulin and triglycerides, as well as insulin resistance, in comparison with their respective controls. The CAF diet reduced motor activity in all behavioural tests, enhanced exploration, reduced anxiety-like behaviour and increased social interaction; this last effect was more pronounced in females than in males. When compared to animals only fed with a CAF diet, CAF withdrawal increased anxiety in the open field, slightly decreased body weight, and completely recovered the liver weight, insulin sensitivity and the standard levels of glucose, insulin and triglycerides in plasma. In conclusion, a CAF diet fed to young animals for 8 weeks induced obesity and metabolic syndrome, and produced robust behavioural changes in young adult rats, whereas CAF withdrawal in the last week modestly increased anxiety, reversed the metabolic alterations and partially reduced overweight.

Assessment of metabolic fluxes in the mouse brain in vivo using 1H-[13C] NMR spectroscopy at 14.1 Tesla.

(13)C magnetic resonance spectroscopy (MRS) combined with the administration of (13)C labeled substrates uniquely allows to measure metabolic fluxes in vivo in the brain of humans and rats. The extension to mouse models may provide exclusive prospect for the investigation of models of human diseases. In the present study, the short-echo-time (TE) full-sensitivity (1)H-[(13)C] MRS sequence combined with high magnetic field (14.1 T) and infusion of [U-(13)C6] glucose was used to enhance the experimental sensitivity in vivo in the mouse brain and the (13)C turnover curves of glutamate C4, glutamine C4, glutamate+glutamine C3, aspartate C2, lactate C3, alanine C3, γ-aminobutyric acid C2, C3 and C4 were obtained. A one-compartment model was used to fit (13)C turnover curves and resulted in values of metabolic fluxes including the tricarboxylic acid (TCA) cycle flux VTCA (1.05 ± 0.04 μmol/g per minute), the exchange flux between 2-oxoglutarate and glutamate Vx (0.48 ± 0.02 μmol/g per minute), the glutamate-glutamine exchange rate V(gln) (0.20 ± 0.02 μmol/g per minute), the pyruvate dilution factor K(dil) (0.82 ± 0.01), and the ratio for the lactate conversion rate and the alanine conversion rate V(Lac)/V(Ala) (10 ± 2). This study opens the prospect of studying transgenic mouse models of brain pathologies.

GD3 Synthase Overexpression Sensitizes Hepatocarcinoma Cells to Hypoxia and Reduces Tumor Growth by Suppressing the cSrc/NF-κB Survival Pathway

Abstract Background: Hypoxia-mediated HIF-1a stabilization and NF-kB activation play a key role in carcinogenesis by fostering cancer cell survival, angiogenesis and tumor invasion. Gangliosides are integral components of biological membranes with an increasingly recognized role as signaling intermediates. In particular, ganglioside GD3 has been characterized as a proapoptotic lipid effector by promoting cell death signaling and suppression of survival pathways. Thus, our aim was to analyze the role of GD3 in hypoxia susceptibility of hepatocarcinoma cells and in vivo tumor growth. Methodology/Principal Findings: We generated and characterized a human hepatocarcinoma cell line stably expressing GD3 synthase (Hep3B-GD3), which catalyzes the synthesis of GD3 from GM3. Despite increased GD3 levels (2-3 fold), no significant changes in cell morphology or growth were observed in Hep3B-GD3 cells compared to wild type Hep3B cells under normoxia. However, exposure of Hep3B-GD3 cells to hypoxia (2% O2) enhanced reactive oxygen species (ROS) generation, resulting in decreased cell survival, with similar findings observed in Hep3B cells exposed to increasing doses of exogenous GD3. In addition, hypoxia-induced c-Src phosphorylation at tyrosine residues, NF-kB activation and subsequent expression of Mn-SOD were observed in Hep3B cells but not in Hep3B-GD3 cells. Moreover, MnTBAP, an antioxidant with predominant SOD mimetic activity, reduced ROS generation, protecting Hep3B-GD3 cells from hypoxia-induced death. Finally, lower tumor growth, higher cell death and reduced Mn-SOD expression were observed in Hep3B-GD3 compared to Hep3B tumor xenografts. Conclusion: These findings underscore a role for GD3 in hypoxia susceptibility by disabling the c-Src/NF-kB survival pathway resulting in lower Mn-SOD expression, which may be of relevance in hepatocellular carcinoma therapy.

Metabolic origin of δ15N values in nitrogenous compounds from Brassica napus L. leaves

Nitrogen isotope composition (δ15N) in plant organic matter is currently used as a natural tracer of nitrogen acquisition efficiency. However, the δ15N value of whole leaf material does not properly reflect the way in which N is assimilated because isotope fractionations along metabolic reactions may cause substantial differences among leaf compounds. In other words, any change in metabolic composition or allocation pattern may cause undesirable variability in leaf δ15N. Here, we investigated the δ15N in different leaf fractions and individual metabolites from rapeseed (Brassica napus) leaves. We show that there were substantial differences in δ15N between nitrogenous compounds (up to 30 ) and the content in (15N enriched) nitrate had a clear influence on leaf δ15N. Using a simple steady-state model of day metabolism, we suggest that the δ15N value in major amino acids was mostly explained by isotope fractionation associated with isotope effects on enzyme-catalysed reactions in primary nitrogen metabolism. δ15N values were further influenced by light versus dark conditions and the probable occurrence of alternative biosynthetic pathways. We conclude that both biochemical pathways (that fractionate between isotopes) and nitrogen sources (used for amino acid production) should be considered when interpreting the δ15N value of leaf nitrogenous compounds

Sphingomonas wittichii RW1 gene reporters interrogating the dibenzofuran metabolic network highlight conditions for early successful development in contaminated microcosms.

In order to better understand the fate and activity of bacteria introduced into contaminated material for the purpose of enhancing biodegradation rates, we constructed Sphingomonas wittichii RW1 variants with gene reporters interrogating dibenzofuran metabolic activity. Three potential promoters from the dibenzofuran metabolic network were selected and fused to the gene for enhanced green fluorescent protein (EGFP). The stability of the resulting genetic constructions in RW1 was examined, with plasmids based on the broad-host range vector pME6012 being the most reliable. One of the selected promoters, upstream of the gene Swit_4925 for a putative 2-hydroxy-2,4-pentadienoate hydratase, was inducible by growth on dibenzofuran. Sphingomonas wittichii RW1 equipped with the Swit_4925 promoter egfp fusion grew in a variety of non-sterile sandy microcosms contaminated with dibenzofuran and material from a former gasification site. The strain also grew in microcosms without added dibenzofuran but to a very limited extent, and EGFP expression indicated the formation of consistent small subpopulations of cells with an active inferred dibenzofuran metabolic network. Evidence was obtained for competition for dibenzofuran metabolites scavenged by resident bacteria in the gasification site material, which resulted in a more rapid decline of the RW1 population. Our results show the importance of low inoculation densities in order to observe the population development of the introduced bacteria and further illustrate that the limited availability of unique carbon substrate may be the most important factor impinging growth.

Clinical evolution after a non-reactive hypothermic EEG following cardiac arrest.

BACKGROUND: Lack of electroencephalography (EEG) background reactivity during therapeutic hypothermia (TH) has been associated with poor outcome in post-anoxic comatose patients. However, decision on intensive care withdrawal is based on normothermic (NT) evaluations. This study aims at exploring whether patients showing recovery of EEG reactivity in NT after a non-reactive EEG in TH differ from those remaining non-reactive. METHODS: Patients with non-reactive EEG during TH were identified from our prospective registry of consecutive comatose adults admitted after successful resuscitation from CA between April 2009 and June 2014. Variables including neurological examination, serum neuron-specific enolase (NSE), procalcitonin, and EEG features were compared regarding impact on functional outcome at 3 months. RESULTS: Seventy-two of 197 patients (37 %) had a non-reactive EEG background during TH with thirteen (18 %) evolving towards reactivity in NT. Compared to those remaining non-reactive (n = 59), they showed significantly better recovery of brainstem reflexes (p < 0.001), better motor responses (p < 0.001), transitory consciousness improvement (p = 0.008), and a tendency toward lower NSE (p = 0.067). One patient recovering EEG reactivity survived with good functional outcome at 3 months. CONCLUSIONS: Recovery of EEG reactivity from TH to NT seems to distinguish two patients' subgroups regarding early neurological assessment and transitory consciousness improvement, corroborating the role of EEG in providing information about cerebral functions. Understanding these dynamic changes encourages maintenance of intensive support in selected patients even after a non-reactive EEG background in TH, as a small subgroup may indeed recover with good functional outcome.

Metabolic effects of insulin and IGFs on gilthead sea bream (Sparus aurata) muscle cells.

Primary cultures of gilthead sea bream myocytes were performed in order to examine the relative metabolic function of insulin compared with IGF-I and IGF-II (insulin-like growth factors, IGFs) at different stages in the cell culture. In these cells, the in vitro effects of insulin and IGFs on 2-deoxyglucose (2-DG) and L-alanine uptake were studied in both myocytes (day 4) and small myotubes (day 9). 2-DG uptake in gilthead sea bream muscle cells was increased in the presence of insulin and IGFs in a time dependent manner and along with muscle cell differentiation. On the contrary, L-alanine uptake was also stimulated by insulin and IGFs but showed an inverse pattern, being the uptake higher in small myocytes than in large myotubes. The results of preincubation with inhibitors (PD-98059, wortmannin, and cytochalasin B) on 2-DG uptake indicated that insulin and IGFs stimulate glucose uptake through the same mechanisms, and evidenced that mitogenesis activator protein kinase (MAPK) and PI3K-Akt transduction pathways mediate the metabolic function of these peptides. In the same way, we observed that GLUT4 protein synthesis was stimulated in the presence of insulin and IGFs in gilthead sea bream muscle cells in a different manner at days 4 or 9 of the culture. In summary we describe here, for the first time, the effects of insulin and IGFs on 2-DG and L-alanine uptake in primary culture of gilthead sea bream muscle cells. We show that both MAPK and PI3K-Akt transduction pathways are needed in order to control insulin and IGFs actions in these cells. Moreover, changes in glucose uptake can be explained by the action of the GLUT4 transporter, which is stimulated in the presence of insulin and IGFs throughout the cell culture.