Soft constraints-based multiobjective framework for flux balance analysis

Models and molecular mechanisms for trade-offs in the context of metabolism

Accumulating evidence for trade-offs involving metabolic traits has demonstrated their importance in the evolution of organisms. Metabolic models with different levels of complexity have already been considered when investigating mechanisms that explain various metabolic trade-offs. Here we provide a systematic review of modelling approaches that have been used to study and explain trade-offs between: (i) the kinetic properties of individual enzymes, (ii) rates of metabolic reactions, (iii) the rate and yield of metabolic pathways and networks, (iv) different metabolic objectives in single organisms and in metabolic communities, and (v) metabolic concentrations. In providing insights into the mechanisms underlying these five types of metabolic trade-offs obtained from constraint-based metabolic modelling, we emphasize the relationship of metabolic trade-offs to the classical black box Y-model that provides a conceptual explanation for resource acquisition-allocation trade-offs. In addition, we identify several pressing concerns and offer a perspective for future research in the identification and manipulation of metabolic trade-offs by relying on the toolbox provided by constraint-based metabolic modelling for single organisms and microbial communities.

A smart RBS library and its prediction model for robust and accurate fine-tuning of gene expression in Bacillus species

Bacillus species, such as Bacillus subtilis and Bacillus licheniformis, are important industrial bacteria. However, there is a lack of standardized and predictable genetic tools for convenient and reproducible assembly of genetic modules in Bacillus species to realize their full potential. In this study, we constructed a Ribosome Binding Site (RBS) library in B. licheniformis, which provides incremental regulation of expression levels over a 104-fold range. Additionally, we developed a model to quantify the resulting translation rates. We successfully demonstrated the robust expression of various target genes using the RBS library and showed that the model accurately predicts the translation rates of arbitrary coding genes. Importantly, we also extended the use of the RBS library and prediction model to B. subtilis, B. thuringiensis, and B. amyloliquefacie. The versatility of the RBS library and its prediction model enables quantification of biological behavior, facilitating reliable forward engineering of gene expression.

Uses of Multi-Objective Flux Analysis for Optimization of Microbial Production of Secondary Metabolites

Secondary metabolites are not essential for the growth of microorganisms, but they play a critical role in how microbes interact with their surroundings. In addition to this important ecological role, secondary metabolites also have a variety of agricultural, medicinal, and industrial uses, and thus the examination of secondary metabolism of plants and microbes is a growing scientific field. While the chemical production of certain secondary metabolites is possible, industrial-scale microbial production is a green and economically attractive alternative. This is even more true, given the advances in bioengineering that allow us to alter the workings of microbes in order to increase their production of compounds of interest. This type of engineering requires detailed knowledge of the "chassis" organism's metabolism. Since the resources and the catalytic capacity of enzymes in microbes is finite, it is important to examine the tradeoffs between various bioprocesses in an engineered system and alter its working in a manner that minimally perturbs the robustness of the system while allowing for the maximum production of a product of interest. The in silico multi-objective analysis of metabolism using genome-scale models is an ideal method for such examinations.

Metabolic collateral lethal target identification reveals MTHFD2 paralogue dependency in ovarian cancer

Recurrent loss-of-function deletions cause frequent inactivation of tumour suppressor genes but often also involve the collateral deletion of essential genes in chromosomal proximity, engendering dependence on paralogues that maintain similar function. Although these paralogues are attractive anticancer targets, no methodology exists to uncover such collateral lethal genes. Here we report a framework for collateral lethal gene identification via metabolic fluxes, CLIM, and use it to reveal MTHFD2 as a collateral lethal gene in UQCR11-deleted ovarian tumours. We show that MTHFD2 has a non-canonical oxidative function to provide mitochondrial NAD+, and demonstrate the regulation of systemic metabolic activity by the paralogue metabolic pathway maintaining metabolic flux compensation. This UQCR11-MTHFD2 collateral lethality is confirmed in vivo, with MTHFD2 inhibition leading to complete remission of UQCR11-deleted ovarian tumours. Using CLIM's machine learning and genome-scale metabolic flux analysis, we elucidate the broad efficacy of targeting MTHFD2 despite distinct cancer genetic profiles co-occurring with UQCR11 deletion and irrespective of stromal compositions of tumours.

Comparison of Optimization-Modelling Methods for Metabolites Production in Escherichia coli

The metabolic network is the reconstruction of the metabolic pathway of an organism that is used to represent the interaction between enzymes and metabolites in genome level. Meanwhile, metabolic engineering is a process that modifies the metabolic network of a cell to increase the production of metabolites. However, the metabolic networks are too complex that cause problem in identifying near-optimal knockout genes/reactions for maximizing the metabolite’s production. Therefore, through constraint-based modelling, various metaheuristic algorithms have been improvised to optimize the desired phenotypes. In this paper, PSOMOMA was compared with CSMOMA and ABCMOMA for maximizing the production of succinic acid in E. coli. Furthermore, the results obtained from PSOMOMA were validated with results from the wet lab experiment.

A non-dominated sorting Differential Search Algorithm Flux Balance Analysis (ndsDSAFBA) for in silico multiobjective optimization in identifying reactions knockout

Metabolic engineering is defined as improving the cellular activities of an organism by manipulating the metabolic, signal or regulatory network. In silico reaction knockout simulation is one of the techniques applied to analyse the effects of genetic perturbations on metabolite production. Many methods consider growth coupling as the objective function, whereby it searches for mutants that maximise the growth and production rate. However, the final goal is to increase the production rate. Furthermore, they produce one single solution, though in reality, cells do not focus on one objective and they need to consider various different competing objectives. In this work, a method, termed ndsDSAFBA (non-dominated sorting Differential Search Algorithm and Flux Balance Analysis), has been developed to find the reaction knockouts involved in maximising the production rate and growth rate of the mutant, by incorporating Pareto dominance concepts. The proposed ndsDSAFBA method was validated using three genome-scale metabolic models. We obtained a set of non-dominated solutions, with each solution representing a different mutant strain. The results obtained were compared with the single objective optimisation (SOO) and multi-objective optimisation (MOO) methods. The results demonstrate that ndsDSAFBA is better than the other methods in terms of production rate and growth rate.

Seeing the wood for the trees: a forest of methods for optimization and omic-network integration in metabolic modelling

Metabolic modelling has entered a mature phase with dozens of methods and software implementations available to the practitioner and the theoretician. It is not easy for a modeller to be able to see the wood (or the forest) for the trees. Driven by this analogy, we here present a 'forest' of principal methods used for constraint-based modelling in systems biology. This provides a tree-based view of methods available to prospective modellers, also available in interactive version at http://modellingmetabolism.net, where it will be kept updated with new methods after the publication of the present manuscript. Our updated classification of existing methods and tools highlights the most promising in the different branches, with the aim to develop a vision of how existing methods could hybridize and become more complex. We then provide the first hands-on tutorial for multi-objective optimization of metabolic models in R. We finally discuss the implementation of multi-view machine learning approaches in poly-omic integration. Throughout this work, we demonstrate the optimization of trade-offs between multiple metabolic objectives, with a focus on omic data integration through machine learning. We anticipate that the combination of a survey, a perspective on multi-view machine learning and a step-by-step R tutorial should be of interest for both the beginner and the advanced user.

Human Systems Biology and Metabolic Modelling: A Review-From Disease Metabolism to Precision Medicine

In cell and molecular biology, metabolism is the only system that can be fully simulated at genome scale. Metabolic systems biology offers powerful abstraction tools to simulate all known metabolic reactions in a cell, therefore providing a snapshot that is close to its observable phenotype. In this review, we cover the 15 years of human metabolic modelling. We show that, although the past five years have not experienced large improvements in the size of the gene and metabolite sets in human metabolic models, their accuracy is rapidly increasing. We also describe how condition-, tissue-, and patient-specific metabolic models shed light on cell-specific changes occurring in the metabolic network, therefore predicting biomarkers of disease metabolism. We finally discuss current challenges and future promising directions for this research field, including machine/deep learning and precision medicine. In the omics era, profiling patients and biological processes from a multiomic point of view is becoming more common and less expensive. Starting from multiomic data collected from patients and N-of-1 trials where individual patients constitute different case studies, methods for model-building and data integration are being used to generate patient-specific models. Coupled with state-of-the-art machine learning methods, this will allow characterizing each patient's disease phenotype and delivering precision medicine solutions, therefore leading to preventative medicine, reduced treatment, and in silico clinical trials.

System-level analysis of metabolic trade-offs during anaerobic photoheterotrophic growth in Rhodopseudomonas palustris

BACKGROUND

Living organisms need to allocate their limited resources in a manner that optimizes their overall fitness by simultaneously achieving several different biological objectives. Examination of these biological trade-offs can provide invaluable information regarding the biophysical and biochemical bases behind observed cellular phenotypes. A quantitative knowledge of a cell system's critical objectives is also needed for engineering of cellular metabolism, where there is interest in mitigating the fitness costs that may result from human manipulation.

RESULTS

To study metabolism in photoheterotrophs, we developed and validated a genome-scale model of metabolism in Rhodopseudomonas palustris, a metabolically versatile gram-negative purple non-sulfur bacterium capable of growing phototrophically on various carbon sources, including inorganic carbon and aromatic compounds. To quantitatively assess trade-offs among a set of important biological objectives during different metabolic growth modes, we used our new model to conduct an 8-dimensional multi-objective flux analysis of metabolism in R. palustris. Our results revealed that phototrophic metabolism in R. palustris is light-limited under anaerobic conditions, regardless of the available carbon source. Under photoheterotrophic conditions, R. palustris prioritizes the optimization of carbon efficiency, followed by ATP production and biomass production rate, in a Pareto-optimal manner. To achieve maximum carbon fixation, cells appear to divert limited energy resources away from growth and toward CO2 fixation, even in the presence of excess reduced carbon. We also found that to achieve the theoretical maximum rate of biomass production, anaerobic metabolism requires import of additional compounds (such as protons) to serve as electron acceptors. Finally, we found that production of hydrogen gas, of potential interest as a candidate biofuel, lowers the cellular growth rates under all circumstances.

CONCLUSIONS

Photoheterotrophic metabolism of R. palustris is primarily regulated by the amount of light it can absorb and not the availability of carbon. However, despite carbon's secondary role as a regulating factor, R. palustris' metabolism strives for maximum carbon efficiency, even when this increased efficiency leads to slightly lower growth rates.

Optimization of Multi-Omic Genome-Scale Models: Methodologies, Hands-on Tutorial, and Perspectives

Genome-scale metabolic models are valuable tools for assessing the metabolic potential of living organisms. Being downstream of gene expression, metabolism is increasingly being used as an indicator of the phenotypic outcome for drugs and therapies. We here present a review of the principal methods used for constraint-based modelling in systems biology, and explore how the integration of multi-omic data can be used to improve phenotypic predictions of genome-scale metabolic models. We believe that the large-scale comparison of the metabolic response of an organism to different environmental conditions will be an important challenge for genome-scale models. Therefore, within the context of multi-omic methods, we describe a tutorial for multi-objective optimization using the metabolic and transcriptomics adaptation estimator (METRADE), implemented in MATLAB. METRADE uses microarray and codon usage data to model bacterial metabolic response to environmental conditions (e.g., antibiotics, temperatures, heat shock). Finally, we discuss key considerations for the integration of multi-omic networks into metabolic models, towards automatically extracting knowledge from such models.

Flux balance analysis predicts Warburg-like effects of mouse hepatocyte deficient in miR-122a

The liver is a vital organ involving in various major metabolic functions in human body. MicroRNA-122 (miR-122) plays an important role in the regulation of liver metabolism, but its intrinsic physiological functions require further clarification. This study integrated the genome-scale metabolic model of hepatocytes and mouse experimental data with germline deletion of Mir122a (Mir122a^–/–) to infer Warburg-like effects. Elevated expression of MiR-122a target genes in Mir122a^–/–mice, especially those encoding for metabolic enzymes, was applied to analyze the flux distributions of the genome-scale metabolic model in normal and deficient states. By definition of the similarity ratio, we compared the flux fold change of the genome-scale metabolic model computational results and metabolomic profiling data measured through a liquid-chromatography with mass spectrometer, respectively, for hepatocytes of 2-month-old mice in normal and deficient states. The Ddc gene demonstrated the highest similarity ratio of 95% to the biological hypothesis of the Warburg effect, and similarity of 75% to the experimental observation. We also used 2, 6, and 11 months of mir-122 knockout mice liver cell to examined the expression pattern of DDC in the knockout mice livers to show upregulated profiles of DDC from the data. Furthermore, through a bioinformatics (LINCS program) prediction, BTK inhibitors and withaferin A could downregulate DDC expression, suggesting that such drugs could potentially alter the early events of metabolomics of liver cancer cells.

For almost a century, researchers have known that cancer cells have an abnormal metabolism and utilize glucose differently than normal cells do. Aerobic glycolysis or the Warburg effect in cancer cells involves elevated glucose uptake with lactic acid production in the presence of oxygen. MicroRNAs have recently been discovered to be key metabolic regulators that mediate the fine tuning of genes that are involved directly or indirectly in cancer metabolism. MicroRNA-122 (miR-122) plays an important role in the regulation of liver metabolism, but its intrinsic physiological functions require further clarification. This study integrated the genome-scale metabolic modeling (GSMM) of hepatocytes and mouse experimental data with germline deletion of Mir122a (Mir122a^–/–) to infer Warburg-like effects. In silico and in vivo observations indicated that DDC overexpression induced Warburg effect in hepatocyte. Furthermore, through a bioinformatics prediction, BTK inhibitors and withaferin A could downregulate DDC expression, suggesting that such drugs could potentially alter the early events of metabolomics of liver cancer cells.

Exo-MFA - A 13C metabolic flux analysis framework to dissect tumor microenvironment-secreted exosome contributions towards cancer cell metabolism

Dissecting the pleiotropic roles of tumor micro-environment (TME) on cancer progression has been brought to the foreground of research on cancer pathology. Extracellular vesicles such as exosomes, transport proteins, lipids, and nucleic acids, to mediate intercellular communication between TME components and have emerged as candidates for anti-cancer therapy. We previously reported that cancer-associated fibroblast (CAF) derived exosomes (CDEs) contain metabolites in their cargo that are utilized by cancer cells for central carbon metabolism and promote cancer growth. However, the metabolic fluxes involved in donor cells towards packaging of metabolites in extracellular vesicles and exosome-mediated metabolite flux upregulation in recipient cells are still not known. Here, we have developed a novel empirical and computational technique, exosome-mediated metabolic flux analysis (Exo-MFA) to quantify flow of cargo from source cells to recipient cells via vesicular transport. Our algorithm, which is based on ¹³C metabolic flux analysis, successfully predicts packaging fluxes to metabolite cargo in CAFs, dynamic changes in rate of exosome internalization by cancer cells, and flux of cargo release over time. We find that cancer cells internalize exosomes rapidly leading to depletion of extracellular exosomes within 24h. However, metabolite cargo significantly alters intracellular metabolism over the course of 24h by regulating glycolysis pathway fluxes via lactate supply. Furthermore, it can supply up to 35% of the TCA cycle fluxes by providing TCA intermediates and glutamine. Our algorithm will help gain insight into (i) metabolic interactions in multicellular systems (ii) biogenesis of extracellular vesicles and their differential packaging of cargo under changing environments, and (iii) regulation of cancer cell metabolism by its microenvironment.

Mapping the landscape of metabolic goals of a cell

Genome-scale flux balance models of metabolism provide testable predictions of all metabolic rates in an organism, by assuming that the cell is optimizing a metabolic goal known as the objective function. We introduce an efficient inverse flux balance analysis (invFBA) approach, based on linear programming duality, to characterize the space of possible objective functions compatible with measured fluxes. After testing our algorithm on simulated E. coli data and time-dependent S. oneidensis fluxes inferred from gene expression, we apply our inverse approach to flux measurements in long-term evolved E. coli strains, revealing objective functions that provide insight into metabolic adaptation trajectories.

MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models

Comparative metabolic modelling is emerging as a novel field, supported by the development of reliable and standardized approaches for constructing genome-scale metabolic models in high throughput. New software solutions are needed to allow efficient comparative analysis of multiple models in the context of multiple cellular objectives. Here, we present the user-friendly software framework Multi-Metabolic Evaluator (MultiMetEval), built upon SurreyFBA, which allows the user to compose collections of metabolic models that together can be subjected to flux balance analysis. Additionally, MultiMetEval implements functionalities for multi-objective analysis by calculating the Pareto front between two cellular objectives. Using a previously generated dataset of 38 actinobacterial genome-scale metabolic models, we show how these approaches can lead to exciting novel insights. Firstly, after incorporating several pathways for the biosynthesis of natural products into each of these models, comparative flux balance analysis predicted that species like Streptomyces that harbour the highest diversity of secondary metabolite biosynthetic gene clusters in their genomes do not necessarily have the metabolic network topology most suitable for compound overproduction. Secondly, multi-objective analysis of biomass production and natural product biosynthesis in these actinobacteria shows that the well-studied occurrence of discrete metabolic switches during the change of cellular objectives is inherent to their metabolic network architecture. Comparative and multi-objective modelling can lead to insights that could not be obtained by normal flux balance analyses. MultiMetEval provides a powerful platform that makes these analyses straightforward for biologists. Sources and binaries of MultiMetEval are freely available from https://github.com/PiotrZakrzewski/MetEval/downloads.

Metabolic flux ratio analysis and multi-objective optimization revealed a globally conserved and coordinated metabolic response of E. coli to paraquat-induced oxidative stress

The ability of a microorganism to adapt to changes in the environment, such as in nutrient or oxygen availability, is essential for its competitive fitness and survival. The cellular objective and the strategy of the metabolic response to an extreme environment are therefore of tremendous interest and, thus, have been increasingly explored. However, the cellular objective of the complex regulatory structure of the metabolic changes has not yet been fully elucidated and more details regarding the quantitative behaviour of the metabolic flux redistribution are required to understand the systems-wide biological significance of this response. In this study, the intracellular metabolic flux ratios involved in the central carbon metabolism were determined by fractional (13)C-labeling and metabolic flux ratio analysis (MetaFoR) of the wild-type E. coli strain JM101 at an oxidative environment in a chemostat. We observed a significant increase in the flux through phosphoenolpyruvate carboxykinase (PEPCK), phosphoenolpyruvate carboxylase (PEPC), malic enzyme (MEZ) and serine hydroxymethyltransferase (SHMT). We applied an ε-constraint based multi-objective optimization to investigate the trade-off relationships between the biomass yield and the generation of reductive power using the in silico iJR904 genome-scale model of E. coli K-12. The theoretical metabolic redistribution supports that the trans-hydrogenase pathway should not play a direct role in the defence mounted by E. coli against oxidative stress. The agreement between the measured ratio and the theoretical redistribution established the significance of NADPH synthesis as the goal of the metabolic reprogramming that occurs in response to oxidative stress. Our work presents a framework that combines metabolic flux ratio analysis and multi-objective optimization to investigate the metabolic trade-offs that occur under varied environmental conditions. Our results led to the proposal that the metabolic response of E. coli to paraquat-induced oxidative stress is globally conserved and coordinated.

Pyruvate uptake is increased in highly invasive ovarian cancer cells under anoikis conditions for anaplerosis, mitochondrial function, and migration

Anoikis resistance, or the ability for cells to live detached from the extracellular matrix, is a property of epithelial cancers. The "Warburg effect," or the preference of cancer cells for glycolysis for their energy production even in the presence of oxygen, has been shown to be evident in various tumors. Since a cancer cell's metastatic ability depends on microenvironmental conditions (nutrients, stromal cells, and vascularization) and is highly variable for different organs, their cellular metabolic fluxes and nutrient demand may show considerable differences. Moreover, a cancer cell's metastatic ability, which is dependent on the stage of cancer, may further create metabolic alterations depending on its microenvironment. Although recent studies have aimed to elucidate cancer cell metabolism under detached conditions, the nutrient demand and metabolic activity of cancer cells under nonadherent conditions remain poorly understood. Additionally, less is known about metabolic alterations in ovarian cancer cells with varying invasive capability under anoikis conditions. We hypothesized that the metabolism of highly invasive ovarian cancer cells in detachment would differ from less invasive ovarian cancer cells and that ovarian cancer cells will have altered metabolism in detached vs. attached conditions. To assess these metabolic differences, we integrated a secretomics-based metabolic footprinting (MFP) approach with mitochondrial bioenergetics. Interestingly, MFP revealed higher pyruvate uptake and oxygen consumption in more invasive ovarian cancer cells than their less invasive counterparts. Furthermore, ATP production was higher in more invasive vs. less invasive ovarian cancer cells in detachment. We found that pyruvate has an effect on highly invasive ovarian cancer cells' migration ability. Our results are the first to demonstrate that higher mitochondrial activity is related to higher ovarian cancer invasiveness under detached conditions. Importantly, our results bring insights regarding the metabolism of cancer cells under nonadherent conditions and could lead to the development of therapies for modulating cancer cell invasiveness.

Flux Analysis of the Trypanosoma brucei Glycolysis Based on a Multiobjective-Criteria Bioinformatic Approach

Trypanosoma brucei is a protozoan parasite of major of interest in discovering new genes for drug targets. This parasite alternates its life cycle between the mammal host(s) (bloodstream form) and the insect vector (procyclic form), with two divergent glucose metabolism amenable to in vitro culture. While the metabolic network of the bloodstream forms has been well characterized, the flux distribution between the different branches of the glucose metabolic network in the procyclic form has not been addressed so far. We present a computational analysis (called Metaboflux) that exploits the metabolic topology of the procyclic form, and allows the incorporation of multipurpose experimental data to increase the biological relevance of the model. The alternatives resulting from the structural complexity of networks are formulated as an optimization problem solved by a metaheuristic where experimental data are modeled in a multiobjective function. Our results show that the current metabolic model is in agreement with experimental data and confirms the observed high metabolic flexibility of glucose metabolism. In addition, Metaboflux offers a rational explanation for the high flexibility in the ratio between final products from glucose metabolism, thsat is, flux redistribution through the malic enzyme steps.

Advanced stoichiometric analysis of metabolic networks of mammalian systems

Metabolic engineering tools have been widely applied to living organisms to gain a comprehensive understanding about cellular networks and to improve cellular properties. Metabolic flux analysis (MFA), flux balance analysis (FBA), and metabolic pathway analysis (MPA) are among the most popular tools in stoichiometric network analysis. Although application of these tools into well-known microbial systems is extensive in the literature, various barriers prevent them from being utilized in mammalian cells. Limited experimental data, complex regulatory mechanisms, and the requirement of more complex nutrient media are some major obstacles in mammalian cell systems. However, mammalian cells have been used to produce therapeutic proteins, to characterize disease states or related abnormal metabolic conditions, and to analyze the toxicological effects of some medicinally important drugs. Therefore, there is a growing need for extending metabolic engineering principles to mammalian cells in order to understand their underlying metabolic functions. In this review article, advanced metabolic engineering tools developed for stoichiometric analysis including MFA, FBA, and MPA are described. Applications of these tools in mammalian cells are discussed in detail, and the challenges and opportunities are highlighted.

Stoichiometry based steady-state hepatic flux analysis: computational and experimental aspects

: The liver has many complex physiological functions, including lipid, protein and carbohydrate metabolism, as well as bile and urea production. It detoxifies toxic substances and medicinal products. It also plays a key role in the onset and maintenance of abnormal metabolic patterns associated with various disease states, such as burns, infections and major traumas. Liver cells have been commonly used in in vitro experiments to elucidate the toxic effects of drugs and metabolic changes caused by aberrant metabolic conditions, and to improve the functions of existing systems, such as bioartificial liver. More recently, isolated liver perfusion systems have been increasingly used to characterize intrinsic metabolic changes in the liver caused by various perturbations, including systemic injury, hepatotoxin exposure and warm ischemia. Metabolic engineering tools have been widely applied to these systems to identify metabolic flux distributions using metabolic flux analysis or flux balance analysis and to characterize the topology of the networks using metabolic pathway analysis. In this context, hepatic metabolic models, together with experimental methodologies where hepatocytes or perfused livers are mainly investigated, are described in detail in this review. The challenges and opportunities are also discussed extensively.

Optimality and thermodynamics determine the evolution of transcriptional regulatory networks

Transcriptional motifs are small regulatory interaction patterns that regulate biological functions in highly-interacting cellular networks. Recently, attempts have been made to explain the significance of transcriptional motifs through dynamic function. However, fundamental questions remain unanswered. Why are certain transcriptional motifs with similar dynamic function abundant while others occur rarely? What are the criteria for topological generalization of these motifs into complex networks? Here, we present a novel paradigm that combines non-equilibrium thermodynamics with multiobjective-optimality for network analysis. We found that energetic cost, defined herein as specific dissipation energy, is minimal at the optimal environmental conditions and it correlates inversely with the abundance of the network motifs obtained experimentally for E. coli and S. cerevisiae. This yields evidence that dissipative energetics is the underlying criteria used during evolution for motif selection and that biological systems during transcription tend towards evolutionary selection of subgraphs which produces minimum specific heat dissipation under optimal conditions, thereby explaining the abundance/rare occurrence of some motifs. We show that although certain motifs had similar dynamical functionality, they had significantly different energetic cost, thus explaining the abundance/rare occurrence of these motifs. The presented insights may establish global thermodynamic analysis as a backbone in designing and understanding complex networks systems, such as metabolic and protein interaction networks.

Metabolic profiling based quantitative evaluation of hepatocellular metabolism in presence of adipocyte derived extracellular matrix

The elucidation of the effect of extracellular matrices on hepatocellular metabolism is critical to understand the mechanism of functional upregulation. We have developed a system using natural extracellular matrices [Adipogel] for enhanced albumin synthesis of rat hepatocyte cultures for a period of 10 days as compared to collagen sandwich cultures. Primary rat hepatocytes isolated from livers of female Lewis rats recover within 4 days of culture from isolation induced injury while function is stabilized at 7 days post-isolation. Thus, the culture period can be classified into three distinct stages viz. recovery stage [day 0–4], pre-stable stage [day 5–7] and the stable stage [day 8–10]. A Metabolic Flux Analysis of primary rat hepatocytes cultured in Adipogel was performed to identify the key metabolic pathways modulated as compared to collagen sandwich cultures. In the recovery stage [day 4], the collagen-soluble Adipogel cultures shows an increase in TriCarboxylic Acid [TCA] cycle fluxes; in the pre-stable stage [day 7], there is an increase in PPP and TCA cycle fluxes while in the stable stage [day 10], there is a significant increase in TCA cycle, urea cycle fluxes and amino acid uptake rates concomitant with increased albumin synthesis rate as compared to collagen sandwich cultures throughout the culture period. Metabolic analysis of the collagen-soluble Adipogel condition reveals significantly higher transamination reaction fluxes, amino acid uptake and albumin synthesis rates for the stable vs. recovery stages of culture. The identification of metabolic pathways modulated for hepatocyte cultures in presence of Adipogel will be a useful step to develop an optimization algorithm to further improve hepatocyte function for Bioartificial Liver Devices. The development of this framework for upregulating hepatocyte function in Bioartificial Liver Devices will facilitate the utilization of an integrated experimental and computational approach for broader applications of Adipogel in tissue e engineering and regenerative medicine.