The genome-scale metabolic extreme pathway structure in Haemophilus influenzae shows significant network redundancy

Development and applications of genome-scale metabolic network models

The genome-scale metabolic network model is an effective tool for characterizing the gene-protein-response relationship in the entire metabolic pathway of an organism. By combining various algorithms, the genome-scale metabolic network model can effectively simulate the influence of a specific environment on the physiological state of cells, optimize the culture conditions of strains, and predict the targets of genetic modification to achieve targeted modification of strains. In this review, we summarize the whole process of model building, sort out the various tools that may be involved in the model building process, and explain the role of various algorithms in model analysis. In addition, we also summarized the application of GSMM in network characteristics, cell phenotypes, metabolic engineering, etc. Finally, we discuss the current challenges facing GSMM.

Nutritional Interactions between Bacterial Species Colonising the Human Nasal Cavity: Current Knowledge and Future Prospects

The human nasal microbiome can be a reservoir for several pathogens, including Staphylococcus aureus. However, certain harmless nasal commensals can interfere with pathogen colonisation, an ability that could be exploited to prevent infection. Although attractive as a prophylactic strategy, manipulation of nasal microbiomes to prevent pathogen colonisation requires a better understanding of the molecular mechanisms of interaction that occur between nasal commensals as well as between commensals and pathogens. Our knowledge concerning the mechanisms of pathogen exclusion and how stable community structures are established is patchy and incomplete. Nutrients are scarce in nasal cavities, which makes competitive or mutualistic traits in nutrient acquisition very likely. In this review, we focus on nutritional interactions that have been shown to or might occur between nasal microbiome members. We summarise concepts of nutrient release from complex host molecules and host cells as well as of intracommunity exchange of energy-rich fermentation products and siderophores. Finally, we discuss the potential of genome-based metabolic models to predict complex nutritional interactions between members of the nasal microbiome.

Access to highly specialized growth substrates and production of epithelial immunomodulatory metabolites determine survival of Haemophilus influenzae in human airway epithelial cells

Haemophilus influenzae (Hi) infections are associated with recurring acute exacerbations of chronic respiratory diseases in children and adults including otitis media, pneumonia, chronic obstructive pulmonary disease and asthma. Here, we show that persistence and recurrence of Hi infections are closely linked to Hi metabolic properties, where preferred growth substrates are aligned to the metabolome of human airway epithelial surfaces and include lactate, pentoses, and nucleosides, but not glucose that is typically used for studies of Hi growth in vitro. Enzymatic and physiological investigations revealed that utilization of lactate, the preferred Hi carbon source, required the LldD L-lactate dehydrogenase (conservation: 98.8% of strains), but not the two redox-balancing D-lactate dehydrogenases Dld and LdhA. Utilization of preferred substrates was directly linked to Hi infection and persistence. When unable to utilize L-lactate or forced to rely on salvaged guanine, Hi showed reduced extra- and intra-cellular persistence in a murine model of lung infection and in primary normal human nasal epithelia, with up to 3000-fold attenuation observed in competitive infections. In contrast, D-lactate dehydrogenase mutants only showed a very slight reduction compared to the wild-type strain. Interestingly, acetate, the major Hi metabolic end-product, had anti-inflammatory effects on cultured human tissue cells in the presence of live but not heat-killed Hi, suggesting that metabolic endproducts also influence HI-host interactions. Our work provides significant new insights into the critical role of metabolism for Hi persistence in contact with host cells and reveals for the first time the immunomodulatory potential of Hi metabolites.

Accrual of functional redundancy along the lifespan and its effects on cognition

Despite the necessity to understand how the brain endures the initial stages of age-associated cognitive decline, no brain mechanism has been quantitatively specified to date. The brain may withstand the effects of cognitive aging through redundancy, a design feature in engineered and biological systems, which entails the presence of substitute elements to protect it against failure. Here, we investigated the relationship between functional network redundancy and age over the human lifespan and their interaction with cognition, analyzing resting-state functional MRI images and cognitive measures from 579 subjects. Network-wide redundancy was significantly associated with age, showing a stronger link with age than other major topological measures, presenting a pattern of accumulation followed by old-age decline. Critically, redundancy significantly mediated the association between age and executive function, with lower anti-correlation between age and cognition in subjects with high redundancy. The results suggest that functional redundancy accrues throughout the lifespan, mitigating the effects of age on cognition.

Haemophilus influenzae Glucose Catabolism Leading to Production of the Immunometabolite Acetate Has a Key Contribution to the Host Airway-Pathogen Interplay

Chronic obstructive pulmonary disease (COPD) is characterized by abnormal inflammatory responses and impaired airway immunity, which provides an opportunistic platform for nontypeable Haemophilus influenzae (NTHi) infection. Clinical evidence supports that the COPD airways present increased concentrations of glucose, which may facilitate proliferation of pathogenic bacteria able to use glucose as a carbon source. NTHi metabolizes glucose through respiration-assisted fermentation, leading to the excretion of acetate, formate, and succinate. We hypothesized that such specialized glucose catabolism may be a pathoadaptive trait playing a pivotal role in the NTHi airway infection. To find out whether this is true, we engineered and characterized bacterial mutant strains impaired to produce acetate, formate, or succinate by inactivating the ackA, pflA, and frdA genes, respectively. While the inactivation of the pflA and frdA genes only had minimal physiological effects, the inactivation of the ackA gene affected acetate production and led to reduced bacterial growth, production of lactate under low oxygen tension, and bacterial attenuation in vivo. Moreover, bacterially produced acetate was able to stimulate the expression of inflammatory genes by cultured airway epithelial cells. These results back the notion that the COPD lung supports NTHi growth on glucose, enabling production of fermentative end products acting as immunometabolites at the site of infection. Thus, glucose catabolism may contribute not only to NTHi growth but also to bacterially driven airway inflammation. This information has important implications for developing nonantibiotic antimicrobials, given that airway glucose homeostasis modifying drugs could help prevent microbial infections associated with chronic lung disease.

Towards scaling elementary flux mode computation

While elementary flux mode (EFM) analysis is now recognized as a cornerstone computational technique for cellular pathway analysis and engineering, EFM application to genome-scale models remains computationally prohibitive. This article provides a review of aspects of EFM computation that elucidates bottlenecks in scaling EFM computation. First, algorithms for computing EFMs are reviewed. Next, the impact of redundant constraints, sensitivity to constraint ordering and network compression are evaluated. Then, the advantages and limitations of recent parallelization and GPU-based efforts are highlighted. The article then reviews alternative pathway analysis approaches that aim to reduce the EFM solution space. Despite advances in EFM computation, our review concludes that continued scaling of EFM computation is necessary to apply EFM to genome-scale models. Further, our review concludes that pathway analysis methods that target specific pathway properties can provide powerful alternatives to EFM analysis.

In silico predicted transcriptional regulatory control of steroidogenesis in spawning female fathead minnows (Pimephales promelas)

Oocyte development and maturation (or oogenesis) in spawning female fish is mediated by interrelated transcriptional regulatory and steroidogenesis networks. This study integrates a transcriptional regulatory network (TRN) model of steroidogenic enzyme gene expressions with a flux balance analysis (FBA) model of steroidogenesis. The two models were functionally related. Output from the TRN model (as magnitude gene expression simulated using extreme pathway (ExPa) analysis) was used to re-constrain linear inequality bounds for reactions in the FBA model. This allowed TRN model predictions to impact the steroidogenesis FBA model. These two interrelated models were tested as follows: First, in silico targeted steroidogenic enzyme gene activations in the TRN model showed high co-regulation (67-83%) for genes involved with oocyte growth and development (cyp11a1, cyp17-17,20-lyase, 3β-HSD and cyp19a1a). Whereas, no or low co-regulation corresponded with genes concertedly involved with oocyte final maturation prior to spawning (cyp17-17α-hydroxylase (0%) and 20β-HSD (33%)). Analysis (using FBA) of accompanying steroidogenesis fluxes showed high overlap for enzymes involved with oocyte growth and development versus those involved with final maturation and spawning. Second, the TRN model was parameterized with in vivo changes in the presence/absence of transcription factors (TFs) during oogenesis in female fathead minnows (Pimephales promelas). Oogenesis stages studied included: PreVitellogenic-Vitellogenic, Vitellogenic-Mature, Mature-Ovulated and Ovulated-Atretic stages. Predictions of TRN genes active during oogenesis showed overall elevated expressions for most genes during early oocyte development (PreVitellogenic-Vitellogenic, Vitellogenic-Mature) and post-ovulation (Ovulated-Atretic). Whereas ovulation (Mature-Ovulated) showed highest expression for cyp17-17α-hydroxylase only. FBA showed steroid hormone productions to also follow trends concomitant with steroidogenic enzyme gene expressions. General trends predicted by in silico modeling were similar to those observed in vivo. The integrated computational framework presented was capable of mechanistically representing aspects of reproductive function in fish. This approach can be extended to study reproductive effects under exposure to adverse environmental or anthropogenic stressors.

Inhibition of C298S mutant of human aldose reductase for antidiabetic applications: Evidence from in silico elementary mode analysis of biological network model

Human aldose reductase (hAR) is the key enzyme in sorbitol pathway of glucose utilization and is implicated in the etiology of secondary complications of diabetes, such as, cardiovascular complications, neuropathy, nephropathy, retinopathy, and cataract genesis. It reduces glucose to sorbitol in the presence of NADPH and the major cause of diabetes complications could be the change in the osmotic pressure due to the accumulation of sorbitol. An activated form of hAR (activated hAR or ahAR) poses a potential obstacle in the development of diabetes drugs as hAR-inhibitors are ineffective against ahAR. The therapeutic efficacy of such drugs is compromised when a large fraction of the enzyme (hAR) undergoes conversion to the activated ahAR form as has been observed in the diabetic tissues. In the present study, attempts have been made to employ systems biology strategies to identify the elementary nodes of human polyol metabolic pathway, responsible for normal metabolic states, followed by the identification of natural potent inhibitors of the activated form of hAR represented by the mutant C298S for possible antidiabetic applications. Quantum Mechanical Molecular Mechanical docking strategy was used to determine the probable inhibitors of ahAR. Rosmarinic acid was found as the most potent natural ahAR inhibitor and warrants for experimental validation in the near future.

Identification of metabolic pathways using pathfinding approaches: a systematic review

Metabolic pathways have become increasingly available for various microorganisms. Such pathways have spurred the development of a wide array of computational tools, in particular, mathematical pathfinding approaches. This article can facilitate the understanding of computational analysis of metabolic pathways in genomics. Moreover, stoichiometric and pathfinding approaches in metabolic pathway analysis are discussed. Three major types of studies are elaborated: stoichiometric identification models, pathway-based graph analysis and pathfinding approaches in cellular metabolism. Furthermore, evaluation of the outcomes of the pathways with mathematical benchmarking metrics is provided. This review would lead to better comprehension of metabolism behaviors in living cells, in terms of computed pathfinding approaches.

Network-Guided Key Gene Discovery for a Given Cellular Process

Identification of key genes for a given physiological or pathological process is an essential but still very challenging task for the entire biomedical research community. Statistics-based approaches, such as genome-wide association study (GWAS)- or quantitative trait locus (QTL)-related analysis have already made enormous contributions to identifying key genes associated with a given disease or phenotype, the success of which is however very much dependent on a huge number of samples. Recent advances in network biology, especially network inference directly from genome-scale data and the following-up network analysis, opens up new avenues to predict key genes driving a given biological process or cellular function. Here we review and compare the current approaches in predicting key genes, which have no chances to stand out by classic differential expression analysis, from gene-regulatory, protein-protein interaction, or gene expression correlation networks. We elaborate these network-based approaches mainly in the context of immunology and infection, and urge more usage of correlation network-based predictions. Such network-based key gene discovery approaches driven by information-enriched 'omics' data should be very useful for systematic key gene discoveries for any given biochemical process or cellular function, and also valuable for novel drug target discovery and novel diagnostic, prognostic and therapeutic-efficiency marker prediction for a specific disease or disorder.

Identifying all moiety conservation laws in genome-scale metabolic networks

The stoichiometry of a metabolic network gives rise to a set of conservation laws for the aggregate level of specific pools of metabolites, which, on one hand, pose dynamical constraints that cross-link the variations of metabolite concentrations and, on the other, provide key insight into a cell's metabolic production capabilities. When the conserved quantity identifies with a chemical moiety, extracting all such conservation laws from the stoichiometry amounts to finding all non-negative integer solutions of a linear system, a programming problem known to be NP-hard. We present an efficient strategy to compute the complete set of integer conservation laws of a genome-scale stoichiometric matrix, also providing a certificate for correctness and maximality of the solution. Our method is deployed for the analysis of moiety conservation relationships in two large-scale reconstructions of the metabolism of the bacterium E. coli, in six tissue-specific human metabolic networks, and, finally, in the human reactome as a whole, revealing that bacterial metabolism could be evolutionarily designed to cover broader production spectra than human metabolism. Convergence to the full set of moiety conservation laws in each case is achieved in extremely reduced computing times. In addition, we uncover a scaling relation that links the size of the independent pool basis to the number of metabolites, for which we present an analytical explanation.

Finding elementary flux modes in metabolic networks based on flux balance analysis and flux coupling analysis: application to the analysis of Escherichia coli metabolism

Elementary modes (EMs) are steady-state metabolic flux vectors with minimal set of active reactions. Each EM corresponds to a metabolic pathway. Therefore, studying EMs is helpful for analyzing the production of biotechnologically important metabolites. However, memory requirements for computing EMs may hamper their applicability as, in most genome-scale metabolic models, no EM can be computed due to running out of memory. In this study, we present a method for computing randomly sampled EMs. In this approach, a network reduction algorithm is used for EM computation, which is based on flux balance-based methods. We show that this approach can be used to recover the EMs in the medium- and genome-scale metabolic network models, while the EMs are sampled in an unbiased way. The applicability of such results is shown by computing “estimated” control-effective flux values in Escherichia coli metabolic network.

In silico predicted structural and functional robustness of piscine steroidogenesis

Assessments of metabolic robustness or susceptibility are inherently dependent on quantitative descriptions of network structure and associated function. In this paper a stoichiometric model of piscine steroidogenesis was constructed and constrained with productions of selected steroid hormones. Structural and flux metrics of this in silico model were quantified by calculating extreme pathways and optimal flux distributions (using linear programming). Extreme pathway analysis showed progestin and corticosteroid synthesis reactions to be highly participant in extreme pathways. Furthermore, reaction participation in extreme pathways also fitted a power law distribution (degree exponent γ=2.3), which suggested that progestin and corticosteroid reactions act as 'hubs' capable of generating other functionally relevant pathways required to maintain steady-state functionality of the network. Analysis of cofactor usage (O2 and NADPH) showed progestin synthesis reactions to exhibit high robustness, whereas estrogen productions showed highest energetic demands with low associated robustness to maintain such demands. Linear programming calculated optimal flux distributions showed high heterogeneity of flux values with a near-random power law distribution (degree exponent γ≥2.7). Subsequently, network robustness was tested by assessing maintenance of metabolite flux-sum subject to targeted deletions of rank-ordered (low to high metric) extreme pathway participant and optimal flux reactions. Network robustness was susceptible to deletions of extreme pathway participant reactions, whereas minimal impact of high flux reaction deletion was observed. This analysis shows that the steroid network is susceptible to perturbation of structurally relevant (extreme pathway) reactions rather than those carrying high flux.

Metabolic versatility in Haemophilus influenzae: a metabolomic and genomic analysis

Haemophilus influenzae is a host adapted human pathogen known to contribute to a variety of acute and chronic diseases of the upper and lower respiratory tract as well as the middle ear. At the sites of infection as well as during growth as a commensal the environmental conditions encountered by H. influenzae will vary significantly, especially in terms of oxygen availability, however, the mechanisms by which the bacteria can adapt their metabolism to cope with such changes have not been studied in detail. Using targeted metabolomics the spectrum of metabolites produced during growth of H. influenzae on glucose in RPMI-based medium was found to change from acetate as the main product during aerobic growth to formate as the major product during anaerobic growth. This change in end-product is likely caused by a switch in the major route of pyruvate degradation. Neither lactate nor succinate or fumarate were major products of H. influenzae growth under any condition studied. Gene expression studies and enzyme activity data revealed that despite an identical genetic makeup and very similar metabolite production profiles, H. influenzae strain Rd appeared to favor glucose degradation via the pentose phosphate pathway, while strain 2019, a clinical isolate, showed higher expression of enzymes involved in glycolysis. Components of the respiratory chain were most highly expressed during microaerophilic and anaerobic growth in both strains, but again clear differences existed in the expression of genes associated e.g., with NADH oxidation, nitrate and nitrite reduction in the two strains studied. Together our results indicate that H. influenzae uses a specialized type of metabolism that could be termed “respiration assisted fermentation” where the respiratory chain likely serves to alleviate redox imbalances caused by incomplete glucose oxidation, and at the same time provides a means of converting a variety of compounds including nitrite and nitrate that arise as part of the host defence mechanisms.

Comparison on extreme pathways reveals nature of different biological processes

Background

Constraint-based reconstruction and analysis (COBRA) is used for modeling genome-scale metabolic networks (MNs). In a COBRA model, extreme pathways (ExPas) are the edges of its conical solution space, which is formed by all viable steady-state flux distributions. ExPa analysis has been successfully applied to MNs to reveal their phenotypic capabilities and properties. Recently, the COBRA framework has been extended to transcriptional regulatory networks (TRNs) and transcriptional and translational networks (TTNs), so efforts are needed to determine whether ExPa analysis is also effective on these two types of networks.

Results

In this paper, the ExPas resulting from the COBRA models of E.coli's MN, TRN and TTN were horizontally compared from 5 aspects: (1) Total number and the ratio of their amount to reaction amount; (2) Length distribution; (3) Reaction participation; (4) Correlated reaction sets (CoSets); (5) interconnectivity degree. Significant discrepancies in above properties were observed during the comparison, which reveals the biological natures of different biological processes. Besides, by demonstrating the application of ExPa analysis on E.coli, we provide a practical guidance of an improved approach to compute ExPas on COBRA models of TRNs.

Conclusions

ExPas of E.coli's MN, TRN and TTN have different properties, which are strongly connected with various biological natures of biochemical networks, such as topological structure, specificity and redundancy. Our study shows that ExPas are biologically meaningful on the newborn models and suggests the effectiveness of ExPa analysis on them.

Functional annotation of conserved hypothetical proteins from Haemophilus influenzae Rd KW20

Haemophilus influenzae is a Gram negative bacterium that belongs to the family Pasteurellaceae, causes bacteremia, pneumonia and acute bacterial meningitis in infants. The emergence of multi-drug resistance H. influenzae strain in clinical isolates demands the development of better/new drugs against this pathogen. Our study combines a number of bioinformatics tools for function predictions of previously not assigned proteins in the genome of H. influenzae. This genome was extensively analyzed and found 1,657 functional proteins in which function of 429 proteins are unknown, termed as hypothetical proteins (HPs). Amino acid sequences of all 429 HPs were extensively annotated and we successfully assigned the function to 296 HPs with high confidence. We also characterized the function of 124 HPs precisely, but with less confidence. We believed that sequence of a protein can be used as a framework to explain known functional properties. Here we have combined the latest versions of protein family databases, protein motifs, intrinsic features from the amino acid sequence, pathway and genome context methods to assign a precise function to hypothetical proteins for which no experimental information is available. We found these HPs belong to various classes of proteins such as enzymes, transporters, carriers, receptors, signal transducers, binding proteins, virulence and other proteins. The outcome of this work will be helpful for a better understanding of the mechanism of pathogenesis and in finding novel therapeutic targets for H. influenzae.

Decomposition of flux distributions into metabolic pathways

Genome-scale reconstructions are often used for studying relationships between fundamental components of a metabolic system. In this study, we develop a novel computational method for analyzing predicted flux distributions for metabolic reconstructions. Because chemical reactions may have multiple reactants and products, a directed hypergraph where hyperarcs may have multiple tail vertices and head vertices is a more appropriate representation of the metabolic network than a conventional network. We use this view to represent predicted flux distributions by maximum generalized flows on hypergraphs. We then demonstrate that the generalized hyperflow problem may be transformed to an equivalent network flow problem with side constraints. This transformation allows a flux to be decomposed into chains of reactions. Subsequent analysis of these chains helps to characterize active pathways in a flux distribution. Such characterizations facilitate comparisons of flux distributions for different environmental conditions. The proposed method is applied to compare predicted flux distributions for Salmonella typhimurium to study changes in metabolism that cause enhanced virulence during a space flight. The differences between flux distributions corresponding to normal and enhanced virulence states confirm previous observations concerning infection mechanisms and suggest new pathways for exploration.

Comparing methods for metabolic network analysis and an application to metabolic engineering

Bioinformatics tools have facilitated the reconstruction and analysis of cellular metabolism of various organisms based on information encoded in their genomes. Characterization of cellular metabolism is useful to understand the phenotypic capabilities of these organisms. It has been done quantitatively through the analysis of pathway operations. There are several in silico approaches for analyzing metabolic networks, including structural and stoichiometric analysis, metabolic flux analysis, metabolic control analysis, and several kinetic modeling based analyses. They can serve as a virtual laboratory to give insights into basic principles of cellular functions. This article summarizes the progress and advances in software and algorithm development for metabolic network analysis, along with their applications relevant to cellular physiology, and metabolic engineering with an emphasis on microbial strain optimization. Moreover, it provides a detailed comparative analysis of existing approaches under different categories.

An algorithm for rapid computational construction of metabolic networks: a cholesterol biosynthesis example

Alternative pathways of metabolic networks represent the escape routes that can reduce drug efficacy and can cause severe adverse effects. In this paper we introduce a mathematical algorithm and a coding system for rapid computational construction of metabolic networks. The initial data for the algorithm are the source substrate code and the enzyme/metabolite interaction tables. The major strength of the algorithm is the adaptive coding system of the enzyme-substrate interactions. A reverse application of the algorithm is also possible, when optimisation algorithm is used to compute the enzyme/metabolite rules from the reference network structure. The coding system is user-defined and must be adapted to the studied problem. The algorithm is most effective for computation of networks that consist of metabolites with similar molecular structures. The computation of the cholesterol biosynthesis metabolic network suggests that 89 intermediates can theoretically be formed between lanosterol and cholesterol, only 20 are presently considered as cholesterol intermediates. Alternative metabolites may represent links with other metabolic networks both as precursors and metabolites of cholesterol. A possible cholesterol-by-pass pathway to bile acids metabolism through cholestanol is suggested.

Finding minimal generating set for metabolic network with reversible pathways

Elementary flux modes give a mathematical representation of metabolic pathways in metabolic networks satisfying the constraint of non-decomposability. The large cost of their computation shifts attention to computing a minimal generating set which is a conically independent subset of elementary flux modes. When a metabolic network has reversible reactions and also admits a reversible pathway, the minimal generating set is not unique. A theoretical development and computational framework is provided which outline how to compute the minimal generating set in this case. The method is based on combining existing software to compute the minimal generating set for a "pointed cone" together with standard software to compute the Reduced Row Echelon Form.

Relationship between extreme pathways and structurally minimal pathways

The determination of reaction pathways is one of the most important functions that should be performed in exploring the kinetics of catalyzed chemical reactions or biochemical reactions, the latter being generally catalyzed by enzymes. It is proven that the terms, "type-I extreme pathway" and "structurally minimal pathway", both introduced to characterize the kinetics of a catalyzed reaction are equivalent. These two terms are based on two distinct methodologies, one mainly rooted in convex analysis and the other in graph theory. The equivalence promises further even more effective methods for reaction-pathway identification by synergistic integration of existing ones.

Elementary mode analysis: a useful metabolic pathway analysis tool for reprograming microbial metabolic pathways

Elementary mode analysis is a useful metabolic pathway analysis tool to characterize cellular metabolism. It can identify all feasible metabolic pathways known as elementary modes that are inherent to a metabolic network. Each elementary mode contains a minimal and unique set of enzymatic reactions that can support cellular functions at steady state. Knowledge of all these pathway options enables systematic characterization of cellular phenotypes, analysis of metabolic network properties (e.g. structure, regulation, robustness, and fragility), phenotypic behavior discovery, and rational strain design for metabolic engineering application. This chapter focuses on the application of elementary mode analysis to reprogram microbial metabolic pathways for rational strain design and the metabolic pathway evolution of designed strains.

Metabolic networks evolve towards states of maximum entropy production

A metabolic network can be described by a set of elementary modes or pathways representing discrete metabolic states that support cell function. We have recently shown that in the most likely metabolic state the usage probability of individual elementary modes is distributed according to the Boltzmann distribution law while complying with the principle of maximum entropy production. To demonstrate that a metabolic network evolves towards such state we have carried out adaptive evolution experiments with Thermoanaerobacterium saccharolyticum operating with a reduced metabolic functionality based on a reduced set of elementary modes. In such reduced metabolic network metabolic fluxes can be conveniently computed from the measured metabolite secretion pattern. Over a time span of 300 generations the specific growth rate of the strain continuously increased together with a continuous increase in the rate of entropy production. We show that the rate of entropy production asymptotically approaches the maximum entropy production rate predicted from the state when the usage probability of individual elementary modes is distributed according to the Boltzmann distribution. Therefore, the outcome of evolution of a complex biological system can be predicted in highly quantitative terms using basic statistical mechanical principles.

Metabolic network reconstruction and genome-scale model of butanol-producing strain Clostridium beijerinckii NCIMB 8052

BACKGROUND

Solventogenic clostridia offer a sustainable alternative to petroleum-based production of butanol--an important chemical feedstock and potential fuel additive or replacement. C. beijerinckii is an attractive microorganism for strain design to improve butanol production because it (i) naturally produces the highest recorded butanol concentrations as a byproduct of fermentation; and (ii) can co-ferment pentose and hexose sugars (the primary products from lignocellulosic hydrolysis). Interrogating C. beijerinckii metabolism from a systems viewpoint using constraint-based modeling allows for simulation of the global effect of genetic modifications.

RESULTS

We present the first genome-scale metabolic model (iCM925) for C. beijerinckii, containing 925 genes, 938 reactions, and 881 metabolites. To build the model we employed a semi-automated procedure that integrated genome annotation information from KEGG, BioCyc, and The SEED, and utilized computational algorithms with manual curation to improve model completeness. Interestingly, we found only a 34% overlap in reactions collected from the three databases--highlighting the importance of evaluating the predictive accuracy of the resulting genome-scale model. To validate iCM925, we conducted fermentation experiments using the NCIMB 8052 strain, and evaluated the ability of the model to simulate measured substrate uptake and product production rates. Experimentally observed fermentation profiles were found to lie within the solution space of the model; however, under an optimal growth objective, additional constraints were needed to reproduce the observed profiles--suggesting the existence of selective pressures other than optimal growth. Notably, a significantly enriched fraction of actively utilized reactions in simulations--constrained to reflect experimental rates--originated from the set of reactions that overlapped between all three databases (P = 3.52 × 10-9, Fisher's exact test). Inhibition of the hydrogenase reaction was found to have a strong effect on butanol formation--as experimentally observed.

CONCLUSIONS

Microbial production of butanol by C. beijerinckii offers a promising, sustainable, method for generation of this important chemical and potential biofuel. iCM925 is a predictive model that can accurately reproduce physiological behavior and provide insight into the underlying mechanisms of microbial butanol production. As such, the model will be instrumental in efforts to better understand, and metabolically engineer, this microorganism for improved butanol production.

Genome-scale metabolic models of Saccharomyces cerevisiae

Systematic analysis of Saccharomyces cerevisiae metabolic functions and pathways has been the subject of extensive studies and established in many aspects. With the reconstruction of the yeast genome-scale metabolic (GSM) network and in silico simulation of the GSM model, the nature of the underlying cellular processes can be tested and validated with the increasing metabolic knowledge. GSM models are also being exploited in fundamental research studies and industrial applications. In this chapter, the principle concepts for construction, simulation and validation of GSM models, progressive applications of the yeast GSM models, and future perspectives are described. This will support and encourage researchers who are interested in systemic analysis of yeast metabolism and systems biology.

C4GEM, a genome-scale metabolic model to study C4 plant metabolism

Leaves of C(4) grasses (such as maize [Zea mays], sugarcane [Saccharum officinarum], and sorghum [Sorghum bicolor]) form a classical Kranz leaf anatomy. Unlike C(3) plants, where photosynthetic CO(2) fixation proceeds in the mesophyll (M), the fixation process in C(4) plants is distributed between two cell types, the M cell and the bundle sheath (BS) cell. Here, we develop a C(4) genome-scale model (C4GEM) for the investigation of flux distribution in M and BS cells during C(4) photosynthesis. C4GEM, to our knowledge, is the first large-scale metabolic model that encapsulates metabolic interactions between two different cell types. C4GEM is based on the Arabidopsis (Arabidopsis thaliana) model (AraGEM) but has been extended by adding reactions and transporters responsible to represent three different C(4) subtypes (NADP-ME [for malic enzyme], NAD-ME, and phosphoenolpyruvate carboxykinase). C4GEM has been validated for its ability to synthesize 47 biomass components and consists of 1,588 unique reactions, 1,755 metabolites, 83 interorganelle transporters, and 29 external transporters (including transport through plasmodesmata). Reactions in the common C(4) model have been associated with well-annotated C(4) species (NADP-ME subtypes): 3,557 genes in sorghum, 11,623 genes in maize, and 3,881 genes in sugarcane. The number of essential reactions not assigned to genes is 131, 135, and 156 in sorghum, maize, and sugarcane, respectively. Flux balance analysis was used to assess the metabolic activity in M and BS cells during C(4) photosynthesis. Our simulations were consistent with chloroplast proteomic studies, and C4GEM predicted the classical C(4) photosynthesis pathway and its major effect in organelle function in M and BS. The model also highlights differences in metabolic activities around photosystem I and photosystem II for three different C(4) subtypes. Effects of CO(2) leakage were also explored. C4GEM is a viable framework for in silico analysis of cell cooperation between M and BS cells during photosynthesis and can be used to explore C(4) plant metabolism.

Systems biochemistry in practice: experimenting with modelling and understanding, with regulation and control

Biology and medicine have become 'big science', even though we may not always like this: genomics and the subsequent analysis of what the genomes encode has shown that interesting living organisms require many more than 300 gene products to interact. We once thought that somewhere in this jungle of interacting macromolecules was hidden the molecule that constitutes the secret of Life, and therewith of health and disease. Now we know that, somehow, the secret of Life is the jungle of interactions. Consequently, we need to find the Rosetta Stones, i.e. interpretations of this jungle of systems biology. We need to find, perhaps convoluted, paths of understanding and intervention. Systems biochemistry is a good place to start, as it has the foothold that what goes in must come out. In the present paper, we review two strategies, which look at control and regulation. We discuss the difference between control and regulation and prove a relationship between them.

Building and analysing genome-scale metabolic models

Reconstructing a model of the metabolic network of an organism from its annotated genome sequence would seem, at first sight, to be one of the most straightforward tasks in functional genomics, even if the various data sources required were never designed with this application in mind. The number of genome-scale metabolic models is, however, lagging far behind the number of sequenced genomes and is likely to continue to do so unless the model-building process can be accelerated. Two aspects that could usefully be improved are the ability to find the sources of error in a nascent model rapidly, and the generation of tenable hypotheses concerning solutions that would improve a model. We will illustrate these issues with approaches we have developed in the course of building metabolic models of Streptococcus agalactiae and Arabidopsis thaliana.

Applications of genome-scale metabolic reconstructions

The availability and utility of genome-scale metabolic reconstructions have exploded since the first genome-scale reconstruction was published a decade ago. Reconstructions have now been built for a wide variety of organisms, and have been used toward five major ends: (1) contextualization of high-throughput data, (2) guidance of metabolic engineering, (3) directing hypothesis-driven discovery, (4) interrogation of multi-species relationships, and (5) network property discovery. In this review, we examine the many uses and future directions of genome-scale metabolic reconstructions, and we highlight trends and opportunities in the field that will make the greatest impact on many fields of biology.

Systems-level engineering of nonfermentative metabolism in yeast

We designed and experimentally validated an in silico gene deletion strategy for engineering endogenous one-carbon (C1) metabolism in yeast. We used constraint-based metabolic modeling and computer-aided gene knockout simulations to identify five genes (ALT2, FDH1, FDH2, FUM1, and ZWF1), which, when deleted in combination, predicted formic acid secretion in Saccharomyces cerevisiae under aerobic growth conditions. Once constructed, the quintuple mutant strain showed the predicted increase in formic acid secretion relative to a formate dehydrogenase mutant (fdh1 fdh2), while formic acid secretion in wild-type yeast was undetectable. Gene expression and physiological data generated post hoc identified a retrograde response to mitochondrial deficiency, which was confirmed by showing Rtg1-dependent NADH accumulation in the engineered yeast strain. Formal pathway analysis combined with gene expression data suggested specific modes of regulation that govern C1 metabolic flux in yeast. Specifically, we identified coordinated transcriptional regulation of C1 pathway enzymes and a positive flux control coefficient for the branch point enzyme 3-phosphoglycerate dehydrogenase (PGDH). Together, these results demonstrate that constraint-based models can identify seemingly unrelated mutations, which interact at a systems level across subcellular compartments to modulate flux through nonfermentative metabolic pathways.

Abundant indispensable redundancies in cellular metabolic networks

Cellular life is a highly redundant complex system; yet, the evolutionary maintenance of the redundancy remains unexplained. Using a systems biology approach, we infer that 37-47% of metabolic reactions in Escherichia coli and yeast can be individually removed without blocking the production of any biomass component under any nutritional condition. However, the majority of these redundant reactions are preserved because they have differential maximal efficiencies at different conditions or their loss causes an immediate fitness reduction that can only be regained via mutation, drift, and selection in evolution. The remaining redundancies are attributable to pleiotropic effects or recent horizontal gene transfers. We find that E. coli and yeast exhibit opposite relationships between the functional importance and redundancy level of a reaction, which is inconsistent with the conjecture that redundancies are preserved as an adaptation to back up important parts in the system. Interestingly, the opposite relationships can both be recapitulated by a simple model in which the natural environments of the organisms change frequently. Thus, adaptive backup is neither necessary nor sufficient to explain the high redundancy of cellular metabolic networks. Taken together, our results strongly suggest that redundant reactions are not kept as backups and that the genetic robustness of metabolic networks is an evolutionary by-product.

Elementary mode analysis: a useful metabolic pathway analysis tool for characterizing cellular metabolism

Elementary mode analysis is a useful metabolic pathway analysis tool to identify the structure of a metabolic network that links the cellular phenotype to the corresponding genotype. The analysis can decompose the intricate metabolic network comprised of highly interconnected reactions into uniquely organized pathways. These pathways consisting of a minimal set of enzymes that can support steady state operation of cellular metabolism represent independent cellular physiological states. Such pathway definition provides a rigorous basis to systematically characterize cellular phenotypes, metabolic network regulation, robustness, and fragility that facilitate understanding of cell physiology and implementation of metabolic engineering strategies. This mini-review aims to overview the development and application of elementary mode analysis as a metabolic pathway analysis tool in studying cell physiology and as a basis of metabolic engineering.

Genome-scale models of bacterial metabolism: reconstruction and applications

Genome-scale metabolic models bridge the gap between genome-derived biochemical information and metabolic phenotypes in a principled manner, providing a solid interpretative framework for experimental data related to metabolic states, and enabling simple in silico experiments with whole-cell metabolism. Models have been reconstructed for almost 20 bacterial species, so far mainly through expert curation efforts integrating information from the literature with genome annotation. A wide variety of computational methods exploiting metabolic models have been developed and applied to bacteria, yielding valuable insights into bacterial metabolism and evolution, and providing a sound basis for computer-assisted design in metabolic engineering. Recent advances in computational systems biology and high-throughput experimental technologies pave the way for the systematic reconstruction of metabolic models from genomes of new species, and a corresponding expansion of the scope of their applications. In this review, we provide an introduction to the key ideas of metabolic modeling, survey the methods, and resources that enable model reconstruction and refinement, and chart applications to the investigation of global properties of metabolic systems, the interpretation of experimental results, and the re-engineering of their biochemical capabilities.

Modelling and analysis of gene regulatory networks

Gene regulatory networks have an important role in every process of life, including cell differentiation, metabolism, the cell cycle and signal transduction. By understanding the dynamics of these networks we can shed light on the mechanisms of diseases that occur when these cellular processes are dysregulated. Accurate prediction of the behaviour of regulatory networks will also speed up biotechnological projects, as such predictions are quicker and cheaper than lab experiments. Computational methods, both for supporting the development of network models and for the analysis of their functionality, have already proved to be a valuable research tool.

Estimation of the number of extreme pathways for metabolic networks

ABSTRACT: BACKGROUND: The set of extreme pathways (ExPa), {pi}, defines the convex basis vectors used for the mathematical characterization of the null space of the stoichiometric matrix for biochemical reaction networks. ExPa analysis has been used for a number of studies to determine properties of metabolic networks as well as to obtain insight into their physiological and functional states in silico. However, the number of ExPas, p = |{pi}|, grows with the size and complexity of the network being studied, and this poses a computational challenge. For this study, we investigated the relationship between the number of extreme pathways and simple network properties. RESULTS: We established an estimating function for the number of ExPas using these easily obtainable network measurements. In particular, it was found that log [p] had an exponential relationship with log[ summation operatori=1Rd-id+ici] MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacyGGSbaBcqGGVbWBcqGGNbWzdaWadaqaamaaqadabaGaemizaq2aaSbaaSqaaiabgkHiTmaaBaaameaacqWGPbqAaeqaaaWcbeaakiabdsgaKnaaBaaaleaacqGHRaWkdaWgaaadbaGaemyAaKgabeaaaSqabaGccqWGJbWydaWgaaWcbaGaemyAaKgabeaaaeaacqWGPbqAcqGH9aqpcqaIXaqmaeaacqWGsbGua0GaeyyeIuoaaOGaay5waiaaw2faaaaa@4414@, where R = |Reff| is the number of active reactions in a network, d-i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaWgaaWcbaGaeyOeI0YaaSbaaWqaaiabdMgaPbqabaaaleqaaaaa@30A9@ and d+i MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazdaWgaaWcbaGaey4kaSYaaSbaaWqaaiabdMgaPbqabaaaleqaaaaa@309E@ the incoming and outgoing degrees of the reactions ri in Reff, and ci the clustering coefficient for each active reaction. CONCLUSION: This relationship typically gave an estimate of the number of extreme pathways to within a factor of 10 of the true number. Such a function providing an estimate for the total number of ExPas for a given system will enable researchers to decide whether ExPas analysis is an appropriate investigative tool.

Bayesian flux balance analysis applied to a skeletal muscle metabolic model

In this article, the steady state condition for the multi-compartment models for cellular metabolism is considered. The problem is to estimate the reaction and transport fluxes, as well as the concentrations in venous blood when the stoichiometry and bound constraints for the fluxes and the concentrations are given. The problem has been addressed previously by a number of authors, and optimization-based approaches as well as extreme pathway analysis have been proposed. These approaches are briefly discussed here. The main emphasis of this work is a Bayesian statistical approach to the flux balance analysis (FBA). We show how the bound constraints and optimality conditions such as maximizing the oxidative phosphorylation flux can be incorporated into the model in the Bayesian framework by proper construction of the prior densities. We propose an effective Markov chain Monte Carlo (MCMC) scheme to explore the posterior densities, and compare the results with those obtained via the previously studied linear programming (LP) approach. The proposed methodology, which is applied here to a two-compartment model for skeletal muscle metabolism, can be extended to more complex models.

Biochemical and statistical network models for systems biology

The normal and abnormal behavior of a living cell is governed by complex networks of interacting biomolecules. Models of these networks allow us to make predictions about cellular behavior under a variety of environmental cues. In this review, we focus on two broad classes of such models: biochemical network models and statistical inference models. In particular, we discuss a number of modeling approaches in the context of the assumptions that they entail, the types of data required for their inference, and the range of their applicability.

Systematic evaluation of objective functions for predicting intracellular fluxes in Escherichia coli

The in vivo distribution of metabolic fluxes in Escherichia coli can be predicted from optimality principles

At least two different sets of optimality principles govern the operation of the metabolic network under different environmental conditions

Metabolism during unlimited growth on glucose in batch culture is best described by the nonlinear maximization of ATP yield per unit of flux

Based on a long history of biochemical and lately genomic research, metabolic networks, in particular microbial ones, are among the best characterized cellular networks. Most components (genes, proteins and metabolites) and their interactions are known. This topological knowledge of the reaction stoichiometry allows to construct metabolic models up to the level of genome scale (Price et al, 2004). Experimentally, sophisticated ¹³C-tracer-based methodologies were developed that enable tracking of the intracellular flux traffic through the reaction network (Sauer, 2006). With the accumulation of such experimental flux data, the question arises why a particular distribution of flux within the network is realized and not one of many alternatives?

Here, we address the question whether the intracellular flux state can be predicted from optimality principles, with the underlying rational that evolution might have optimized metabolic operation toward particular objectives or combinations of multiple objectives. For this purpose, we performed a systematic and rigorous comparison between computational flux predictions and available experimental flux data (Emmerling et al, 2002; Perrenoud and Sauer, 2005; Nanchen et al, 2006) under six different environmental conditions for the model bacterium E. coli. For computational flux predictions, we used a constraint-based modeling approach that requires a stoichiometric model of metabolism (Stelling, 2004). More specifically, we employed flux balance analysis (FBA) where objective functions are defined that represent optimality principles of network operation (Price et al, 2004). This approach has been applied successfully to predict gene deletion lethality (Edwards and Palsson, 2000a, bEdwards and Palsson, 2000a, b; Forster et al, 2003; Kuepfer et al, 2005), network capacities and feasible network states (Edwards 2001, Ibarra 2002), but in only few cases to predict the intracellular flux state (Beard et al, 2002; Holzhütter, 2004).

While different objective functions were proposed for different biological systems (Holzhütter, 2004; Price et al, 2004; Knorr et al, 2006), by far the most common assumption is that microbial cells maximize their growth. To address this issue more generally, we evaluated the accuracy of FBA-based flux predictions for 11 linear and nonlinear objective functions that were combined with eight adjustable constraints. For this purpose, we constructed a highly interconnected stoichiometric network model with 98 reactions and 60 metabolites of E. coli central carbon metabolism. Based on mathematical analyses, the overall model could be reduced to a set of 10 reactions that summarize the actual systemic degree of freedom.

As a quantitative measure of how accurate the experimental data are predicted, we defined predictive fidelity as a single value to quantify the overall deviation between in silico and in vivo fluxes. By comparing all in silico predictions to ¹³C-based in vivo fluxes, we show that prediction of intracellular steady-state fluxes from network stoichiometry alone is, within limits, possible. An unexpected key result is that no further assumptions on network operation in the form of additional and potentially artificial constraints are necessary, provided the appropriate objective function is chosen for a given condition.

While no single objective was able to describe the flux states under all six conditions, we identified two sets of objectives for biologically meaningful predictions without the need for further constraints. For unlimited growth on glucose in aerobic or nitrate-respiring batch cultures, we find that the most accurate and robust results are obtained with the nonlinear maximization of ATP yield per flux unit (Figure 1). Under nutrient scarcity in glucose- or ammonium-limited continuous cultures, in contrast, linear maximization of the overall ATP or biomass yields achieved the highest predictive accuracy.

Since these identified optimality principles describe the system behavior without preconditioning of the network through further constraints, they reflect, to some extent, the evolutionary selection of metabolic network regulation that realizes the various flux states. For conditions of nutrient scarcity, the maximization of energy or biomass yield objective is consistent with the generally observed physiology (Russell and Cook, 1995). The meaning of the maximization of ATP yield per flux unit objective for unlimited growth, however, is less obvious. Generally, it selects for small networks with yet high, albeit suboptimal ATP formation, which has three biological consequences. Firstly, resources are economically allocated since expenditures for enzyme synthesis are, on average, greater for longer pathways. Secondly, suboptimal ATP yields dissipate more energy and thus enable higher catabolic rates. Thirdly, at a constant catabolic rate, a small network results in shorter residence times of substrate molecules until they generate ATP. The relative contribution of these consequences to the evolution of network regulation is unclear, but simultaneous optimization for ATP yield and catabolic rate under this optimality principle identifies a trade-off between the contradicting objectives of maximum overall ATP yield and maximum rate of ATP formation (Pfeiffer et al, 2001).

To which extent can optimality principles describe the operation of metabolic networks? By explicitly considering experimental errors and in silico alternate optima in flux balance analysis, we systematically evaluate the capacity of 11 objective functions combined with eight adjustable constraints to predict ¹³C-determined in vivo fluxes in Escherichia coli under six environmental conditions. While no single objective describes the flux states under all conditions, we identified two sets of objectives for biologically meaningful predictions without the need for further, potentially artificial constraints. Unlimited growth on glucose in oxygen or nitrate respiring batch cultures is best described by nonlinear maximization of the ATP yield per flux unit. Under nutrient scarcity in continuous cultures, in contrast, linear maximization of the overall ATP or biomass yields achieved the highest predictive accuracy. Since these particular objectives predict the system behavior without preconditioning of the network structure, the identified optimality principles reflect, to some extent, the evolutionary selection of metabolic network regulation that realizes the various flux states.

Integrated energy and flux balance based multiobjective framework for large-scale metabolic networks

Flux balance analysis (FBA) provides a framework for the estimation of intracellular fluxes and energy balance analysis (EBA) ensures the thermodynamic feasibility of the computed optimal fluxes. Previously, these techniques have been used to obtain optimal fluxes that maximize a single objective. Because mammalian systems perform various functions, a multi-objective approach is needed when seeking optimal flux distributions in such systems. For example, hepatocytes perform several metabolic functions at various levels depending on environmental conditions; furthermore, there is a potential benefit to enhance some of these functions for applications such as bioartificial liver (BAL) support devices. Herein we developed a multi-objective optimization approach that couples the normalized Normal Constraint (NC) with both FBA and EBA to obtain multi-objective Pareto-optimal solutions. We investigated the Pareto frontiers in gluconeogenic and glycolytic hepatocytes for various combinations of liver-specific objectives (albumin synthesis, glutathione synthesis, NADPH synthesis, ATP generation, and urea secretion). Next, we evaluated the impact of experimental flux measurements on the Pareto frontiers. We found that measurements induce dramatic changes in Pareto frontiers and further constrain the network fluxes. This multi-objective optimality analysis may help explain certain features of the metabolic control of hepatocytes, which is relevant to the response to hepatocytes and liver to various physiological stimuli and disease states.

Understanding the roadmap of metabolism by pathway analysis

The theoretical investigation of the structure of metabolic systems has recently attracted increasing interest. In this chapter, the basic concepts of metabolic pathway analysis are described and various applications are outlined. In particular, the concepts of nullspace and elementary flux modes are explained. The presentation is illustrated by a simple example from tyrosine metabolism and a system describing lysine production in Corynebacterium glutamicum. The latter system gives rise to 37 elementary modes, 36 of which produce lysine with different molar yields. The examples illustrate that metabolic pathway analysis is a useful tool for better understanding the complex architecture of intracellular metabolism, for determining the pathways on which the molar conversion yield of a substrate-product pair under study is maximal, and for assigning functions to orphan genes (functional genomics). Moreover, problems emerging in the modeling of large networks are discussed. An outlook on current trends in the field concludes the chapter.

Uncovering a hidden distributed architecture behind scale-free transcriptional regulatory networks

Numerous studies in both prokaryotes and eukaryotes have shown that, under standard growth conditions, less than 20% of the protein-coding genes are essential for survival. This suggests that biological systems have evolved to have a high degree of robustness to mutational disruptions that can affect the majority of their genes. This mutational robustness could arise either due to redundancy, i.e. direct backup, or due to distributed architecture, i.e. indirect backup where multiple genes contribute to the functioning of a process in the system. Despite clear evidence for direct backup, the prevalence of indirect backup is poorly understood. In this study, we reveal the existence of a hidden distributed architecture behind the scale-free transcriptional regulatory network of yeast by applying a unique network transformation procedure and show that the network is tolerant even to mutations that disrupt regulatory hubs. Contrary to what is generally accepted, our observation that hubs can be lost or replaced in evolution suggests that this hidden distributed architecture behind scale-free networks protects the overall transcriptional program of the organism from mutations affecting major regulatory hubs. We show that the distributed architecture has been provided by an unexpectedly large number of coordinating partners for any regulatory protein. On the basis of these findings, we propose that the existence of such architecture can allow organisms to explore the adaptive landscape in changing environments by providing the plasticity required to reprogram levels of expression of specific genes that may enhance survival. Thus, an "over-engineered" backup system in the form of distributed architecture is likely to be a major determinant of the "evolvability" of the gene expression in organisms faced with environmental diversity.

Integrative analysis of metabolic networks: from peaks to flux models?

Recent developments in genomic and post-genomic technologies have led to the amassment of data describing genome sequences, transcript, protein and metabolite abundances, protein modifications, and protein-protein and protein-DNA interactions. Such technologies have vastly expanded the inventory of detectable molecular species and can be used to describe their interdependence, but they have yet to fulfill their promise in enhancing our knowledge of how flux through metabolic pathways is regulated. A convergence of traditional reductionistic and novel holistic experimental approaches could aid in elucidating flux control.

Systems analyses characterize integrated functions of biochemical networks

Metabolic, regulatory and signaling pathways have been characterized in detail over the past century. As the amount of genomic, proteomic and metabolic data has increased, and the mathematical and analytical capabilities of interrogating these data have advanced, the overlapping roles of pathway constituents have been described. These developments reflect the truly integrated nature of subcellular biochemical networks. Systems analyses, including the reconstruction of stoichiometric networks, provide a key set of tools for quantifying overlap among the metabolic, regulatory and signaling functions of network components. Accounting for this integration is crucial for accurately describing the function of biochemical networks.