Metabolic reconstruction and modeling of nitrogen fixation in Rhizobium etli

Transcriptomic Response of Rhizobium leguminosarum to Acidic Stress and Nutrient Limitation Is Versatile and Substantially Influenced by Extrachromosomal Gene Pool

Multipartite genomes are thought to confer evolutionary advantages to bacteria by providing greater metabolic flexibility in fluctuating environments and enabling rapid adaptation to new ecological niches and stress conditions. This genome architecture is commonly found in plant symbionts, including nitrogen-fixing rhizobia, such as Rhizobium leguminosarum bv. trifolii TA1 (RtTA1), whose genome comprises a chromosome and four extrachromosomal replicons (ECRs). In this study, the transcriptomic responses of RtTA1 to partial nutrient limitation and low acidic pH were analyzed using high-throughput RNA sequencing. RtTA1 growth under these conditions resulted in the differential expression of 1035 to 1700 genes (DEGs), which were assigned to functional categories primarily related to amino acid and carbohydrate metabolism, ribosome and cell envelope biogenesis, signal transduction, and transcription. These results highlight the complexity of the bacterial response to stress. Notably, the distribution of DEGs among the replicons indicated that ECRs played a significant role in the stress response. The transcriptomic data align with the Rhizobium pangenome analysis, which revealed an over-representation of functional categories related to transport, metabolism, and regulatory functions on ECRs. These findings confirm that ECRs contribute substantially to the ability of rhizobia to adapt to challenging environmental conditions.

A genome-scale metabolic reconstruction of soybean and Bradyrhizobium diazoefficiens reveals the cost-benefit of nitrogen fixation

Nitrogen-fixing symbioses allow legumes to thrive in nitrogen-poor soils at the cost of diverting some photoassimilate to their microsymbionts. Effort is being made to bioengineer nitrogen fixation into nonleguminous crops. This requires a quantitative understanding of its energetic costs and the links between metabolic variations and symbiotic efficiency. A whole-plant metabolic model for soybean (Glycine max) with its associated microsymbiont Bradyrhizobium diazoefficiens was developed and applied to predict the cost-benefit of nitrogen fixation with varying soil nitrogen availability. The model predicted a nitrogen-fixation cost of c. 4.13 g C g-1 N, which when implemented into a crop scale model, translated to a grain yield reduction of 27% compared with a non-nodulating plant receiving its nitrogen from the soil. Considering the lower nitrogen content of cereals, the yield cost to a hypothetical N-fixing cereal is predicted to be less than half that of soybean. Soybean growth was predicted to be c. 5% greater when the nodule nitrogen export products were amides versus ureides. This is the first metabolic reconstruction in a tropical crop species that simulates the entire plant and nodule metabolism. Going forward, this model will serve as a tool to investigate carbon use efficiency and key mechanisms within N-fixing symbiosis in a tropical species forming determinate nodules.

Genome-Scale Metabolic Modelling of Lifestyle Changes in Rhizobium leguminosarum

Biological nitrogen fixation in rhizobium-legume symbioses is of major importance for sustainable agricultural practices. To establish a mutualistic relationship with their plant host, rhizobia transition from free-living bacteria in soil to growth down infection threads inside plant roots and finally differentiate into nitrogen-fixing bacteroids. We reconstructed a genome-scale metabolic model for Rhizobium leguminosarum and integrated the model with transcriptome, proteome, metabolome, and gene essentiality data to investigate nutrient uptake and metabolic fluxes characteristic of these different lifestyles. Synthesis of leucine, polyphosphate, and AICAR is predicted to be important in the rhizosphere, while myo-inositol catabolism is active in undifferentiated nodule bacteria in agreement with experimental evidence. The model indicates that bacteroids utilize xylose and glycolate in addition to dicarboxylates, which could explain previously described gene expression patterns. Histidine is predicted to be actively synthesized in bacteroids, consistent with transcriptome and proteome data for several rhizobial species. These results provide the basis for targeted experimental investigation of metabolic processes specific to the different stages of the rhizobium-legume symbioses. IMPORTANCE Rhizobia are soil bacteria that induce nodule formation on plant roots and differentiate into nitrogen-fixing bacteroids. A detailed understanding of this complex symbiosis is essential for advancing ongoing efforts to engineer novel symbioses with cereal crops for sustainable agriculture. Here, we reconstruct and validate a genome-scale metabolic model for Rhizobium leguminosarum bv. viciae 3841. By integrating the model with various experimental data sets specific to different stages of symbiosis formation, we elucidate the metabolic characteristics of rhizosphere bacteria, undifferentiated bacteria inside root nodules, and nitrogen-fixing bacteroids. Our model predicts metabolic flux patterns for these three distinct lifestyles, thus providing a framework for the interpretation of genome-scale experimental data sets and identifying targets for future experimental studies.

Crystal structures of the elusive Rhizobium etli L-asparaginase reveal a peculiar active site

Rhizobium etli, a nitrogen-fixing bacterial symbiont of legume plants, encodes an essential L-asparaginase (ReAV) with no sequence homology to known enzymes with this activity. High-resolution crystal structures of ReAV show indeed a structurally distinct, dimeric enzyme, with some resemblance to glutaminases and β-lactamases. However, ReAV has no glutaminase or lactamase activity, and at pH 9 its allosteric asparaginase activity is relatively high, with Km for L-Asn at 4.2 mM and kcat of 438 s-1. The active site of ReAV, deduced from structural comparisons and confirmed by mutagenesis experiments, contains a highly specific Zn2+ binding site without a catalytic role. The extensive active site includes residues with unusual chemical properties. There are two Ser-Lys tandems, all connected through a network of H-bonds to the Zn center, and three tightly bound water molecules near Ser48, which clearly indicate the catalytic nucleophile.

Reconstruction and analysis of a genome-scale metabolic model for Agrobacterium tumefaciens

The plant pathogen Agrobacterium tumefaciens causes crown gall disease and is a widely used tool for generating transgenic plants owing to its virulence. The pathogenic process involves a shift from an independent to a living form within a host plant. However, comprehensive analyses of metabolites, genes, and reactions contributing to this complex process are lacking. To gain new insights about the pathogenicity from the viewpoints of physiology and cellular metabolism, a genome-scale metabolic model (GSMM) was reconstructed for A. tumefaciens. The model, referred to as iNX1344, contained 1,344 genes, 1,441 reactions, and 1,106 metabolites. It was validated by analyses of in silico cell growth on 39 unique carbon or nitrogen sources and the flux distribution of carbon metabolism. A. tumefaciens metabolic characteristics under three ecological niches were modelled. A high capacity to access and metabolize nutrients is more important for rhizosphere colonization than in the soil, and substantial metabolic changes were detected during the shift from the rhizosphere to tumour environments. Furthermore, by integrating transcriptome data for tumour conditions, significant alterations in central metabolic pathways and secondary metabolite metabolism were identified. Overall, the GSMM and constraint-based analysis could decode the physiological and metabolic features of A. tumefaciens as well as interspecific interactions with hosts, thereby improving our understanding of host adaptation and infection mechanisms.

How rhizobia adapt to the nodule environment

Rhizobia are a phylogenetically diverse group of soil bacteria that engage in mutualistic interactions with legume plants. Although specifics of the symbioses differ between strains and plants, all symbioses ultimately result in the formation of specialized root nodule organs which host the nitrogen-fixing microsymbionts called bacteroids. Inside nodules, bacteroids encounter unique conditions that necessitate global reprogramming of physiological processes and rerouting of their metabolism. Decades of research have addressed these questions using genetics, omics approaches, and more recently computational modelling. Here we discuss the common adaptations of rhizobia to the nodule environment that define the core principles of bacteroid functioning. All bacteroids are growth-arrested and perform energy-intensive nitrogen fixation fueled by plant-provided C₄-dicarboxylates at nanomolar oxygen levels. At the same time, bacteroids are subject to host control and sanctioning that ultimately determine their fitness and have fundamental importance for the evolution of a stable mutualistic relationship.

Understanding the host-microbe interactions using metabolic modeling

The human gut harbors an enormous number of symbiotic microbes, which is vital for human health. However, interactions within the complex microbiota community and between the microbiota and its host are challenging to elucidate, limiting development in the treatment for a variety of diseases associated with microbiota dysbiosis. Using in silico simulation methods based on flux balance analysis, those interactions can be better investigated. Flux balance analysis uses an annotated genome-scale reconstruction of a metabolic network to determine the distribution of metabolic fluxes that represent the complete metabolism of a bacterium in a certain metabolic environment such as the gut. Simulation of a set of bacterial species in a shared metabolic environment can enable the study of the effect of numerous perturbations, such as dietary changes or addition of a probiotic species in a personalized manner. This review aims to introduce to experimental biologists the possible applications of flux balance analysis in the host-microbiota interaction field and discusses its potential use to improve human health. Video abstract.

Quantitative models of nitrogen-fixing organisms

Nitrogen-fixing organisms are of importance to the environment, providing bioavailable nitrogen to the biosphere. Quantitative models have been used to complement the laboratory experiments and in situ measurements, where such evaluations are difficult or costly. Here, we review the current state of the quantitative modeling of nitrogen-fixing organisms and ways to enhance the bridge between theoretical and empirical studies.

The Impacts of Domestication and Agricultural Practices on Legume Nutrient Acquisition Through Symbiosis With Rhizobia and Arbuscular Mycorrhizal Fungi

Legumes are unique among plants as they can obtain nitrogen through symbiosis with nitrogen-fixing rhizobia that form root nodules in the host plants. Therefore they are valuable crops for sustainable agriculture. Increasing nitrogen fixation efficiency is not only important for achieving better plant growth and yield, but it is also crucial for reducing the use of nitrogen fertilizer. Arbuscular mycorrhizal fungi (AMF) are another group of important beneficial microorganisms that form symbiotic relationships with legumes. AMF can promote host plant growth by providing mineral nutrients and improving the soil ecosystem. The trilateral legume-rhizobia-AMF symbiotic relationships also enhance plant development and tolerance against biotic and abiotic stresses. It is known that domestication and agricultural activities have led to the reduced genetic diversity of cultivated germplasms and higher sensitivity to nutrient deficiencies in crop plants, but how domestication has impacted the capability of legumes to establish beneficial associations with rhizospheric microbes (including rhizobia and fungi) is not well-studied. In this review, we will discuss the impacts of domestication and agricultural practices on the interactions between legumes and soil microbes, focusing on the effects on AMF and rhizobial symbioses and hence nutrient acquisition by host legumes. In addition, we will summarize the genes involved in legume-microbe interactions and studies that have contributed to a better understanding of legume symbiotic associations using metabolic modeling.

Computing irreversible minimal cut sets in genome-scale metabolic networks via flux cone projection

MOTIVATION

Minimal cut sets (MCSs) for metabolic networks are sets of reactions which, if they are removed from the network, prevent a target reaction from carrying flux. To compute MCSs different methods exist, which may fail to find sufficiently many MCSs for larger genome-scale networks.

RESULTS

Here we introduce irreversible minimal cut sets (iMCSs). These are MCSs that consist of irreversible reactions only. The advantage of iMCSs is that they can be computed by projecting the flux cone of the metabolic network on the set of irreversible reactions, which usually leads to a smaller cone. Using oriented matroid theory, we show how the projected cone can be computed efficiently and how this can be applied to find iMCSs even in large genome-scale networks.

AVAILABILITY AND IMPLEMENTATION

Software is freely available at https://sourceforge.net/projects/irreversibleminimalcutsets/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-scale metabolic reconstruction of the symbiosis between a leguminous plant and a nitrogen-fixing bacterium

The mutualistic association between leguminous plants and endosymbiotic rhizobial bacteria is a paradigmatic example of a symbiosis driven by metabolic exchanges. Here, we report the reconstruction and modelling of a genome-scale metabolic network of Medicago truncatula (plant) nodulated by Sinorhizobium meliloti (bacterium). The reconstructed nodule tissue contains five spatially distinct developmental zones and encompasses the metabolism of both the plant and the bacterium. Flux balance analysis (FBA) suggests that the metabolic costs associated with symbiotic nitrogen fixation are primarily related to supporting nitrogenase activity, and increasing N2-fixation efficiency is associated with diminishing returns in terms of plant growth. Our analyses support that differentiating bacteroids have access to sugars as major carbon sources, ammonium is the main nitrogen export product of N2-fixing bacteria, and N2 fixation depends on proton transfer from the plant cytoplasm to the bacteria through acidification of the peribacteroid space. We expect that our model, called 'Virtual Nodule Environment' (ViNE), will contribute to a better understanding of the functioning of legume nodules, and may guide experimental studies and engineering of symbiotic nitrogen fixation.

Single-Cell Metabolic Profiling: Metabolite Formulas from Isotopic Fine Structures in Heterogeneous Plant Cell Populations

Characterization of the metabolic heterogeneity in cell populations requires the analysis of single cells. Most current methods in single-cell analysis rely on cell manipulation, potentially altering the abundance of metabolites in individual cells. A small sample volume and the chemical diversity of metabolites are additional challenges in single-cell metabolomics. Here, we describe the combination of fiber-based laser ablation electrospray ionization (f-LAESI) with 21 T Fourier transform ion cyclotron resonance mass spectrometry (21TFTICR-MS) for in situ single-cell metabolic profiling in plant tissue. Single plant cells infected by bacteria were selected and sampled directly from the tissue without cell manipulation through mid-infrared ablation with a fine optical fiber tip for ionization by f-LAESI. Ultrahigh performance 21T-FTICR-MS enabled the simultaneous capture of isotopic fine structures (IFSs) for 47 known and 11 unknown compounds, thus elucidating their elemental compositions from single cells and providing information on metabolic heterogeneity in the cell population.

Metabolic Analyses of Nitrogen Fixation in the Soybean Microsymbiont Sinorhizobium fredii Using Constraint-Based Modeling

Nitrogen is the most limiting macronutrient for plant growth, and rhizobia are important bacteria for agriculture because they can fix atmospheric nitrogen and make it available to legumes through the establishment of a symbiotic relationship with their host plants. In this work, we studied the nitrogen fixation process in the microsymbiont Sinorhizobium fredii at the genome level. A metabolic model was built using genome annotation and literature to reconstruct the symbiotic form of S. fredii. Genes controlling the nitrogen fixation process were identified by simulating gene knockouts. Additionally, the nitrogen-fixing capacities of S. fredii CCBAU45436 in symbiosis with cultivated and wild soybeans were evaluated. The predictions suggested an outperformance of S. fredii with cultivated soybean, consistent with published experimental evidence. The reconstruction presented here will help to understand and improve nitrogen fixation capabilities of S. fredii and will be beneficial for agriculture by reducing the reliance on fertilizer applications.

Rhizobia are soil bacteria able to establish symbiosis with diverse host plants. Specifically, Sinorhizobium fredii is a soil bacterium that forms nitrogen-fixing root nodules in diverse legumes, including soybean. The strain S. fredii CCBAU45436 is a dominant sublineage of S. fredii that nodulates soybeans in alkaline-saline soils in the Huang-Huai-Hai Plain region of China. Here, we present a manually curated metabolic model of the symbiotic form of Sinorhizobium fredii CCBAU45436. A symbiosis reaction was defined to describe the specific soybean-microsymbiont association. The performance and quality of the reconstruction had a 70% score when assessed using a standardized genome-scale metabolic model test suite. The model was used to evaluate in silico single-gene knockouts to determine the genes controlling the nitrogen fixation process. One hundred forty-one of 541 genes (26%) were found to influence the symbiotic process, wherein 121 genes were predicted as essential and 20 others as having a partial effect. Transcriptomic profiles of CCBAU45436 were used to evaluate the nitrogen fixation capacity in cultivated versus in wild soybean inoculated with the microsymbiont. The model quantified the nitrogen fixation activities of the strain in these two hosts and predicted a higher nitrogen fixation capacity in cultivated soybean. Our results are consistent with published data demonstrating larger amounts of ureides and total nitrogen in cultivated soybean than in wild soybean. This work presents the first metabolic network reconstruction of S. fredii as an example of a useful tool for exploring the potential benefits of microsymbionts to sustainable agriculture and the ecosystem.

IMPORTANCE Nitrogen is the most limiting macronutrient for plant growth, and rhizobia are important bacteria for agriculture because they can fix atmospheric nitrogen and make it available to legumes through the establishment of a symbiotic relationship with their host plants. In this work, we studied the nitrogen fixation process in the microsymbiont Sinorhizobium fredii at the genome level. A metabolic model was built using genome annotation and literature to reconstruct the symbiotic form of S. fredii. Genes controlling the nitrogen fixation process were identified by simulating gene knockouts. Additionally, the nitrogen-fixing capacities of S. fredii CCBAU45436 in symbiosis with cultivated and wild soybeans were evaluated. The predictions suggested an outperformance of S. fredii with cultivated soybean, consistent with published experimental evidence. The reconstruction presented here will help to understand and improve nitrogen fixation capabilities of S. fredii and will be beneficial for agriculture by reducing the reliance on fertilizer applications.

MICOM: Metagenome-Scale Modeling To Infer Metabolic Interactions in the Gut Microbiota

The bacterial communities that live within the human gut have been linked to health and disease. However, we are still just beginning to understand how those bacteria interact and what potential interventions to our gut microbiome can make us healthier. Here, we present a mathematical modeling framework (named MICOM) that can recapitulate the growth rates of diverse bacterial species in the gut and can simulate metabolic interactions within microbial communities. We show that MICOM can unravel the ecological rules that shape the microbial landscape in our gut and that a given dietary or probiotic intervention can have widely different effects in different people.

Compositional changes in the gut microbiota have been associated with a variety of medical conditions such as obesity, Crohn’s disease, and diabetes. However, connecting microbial community composition to ecosystem function remains a challenge. Here, we introduce MICOM, a customizable metabolic model of the human gut microbiome. By using a heuristic optimization approach based on L2 regularization, we were able to obtain a unique set of realistic growth rates that corresponded well with observed replication rates. We integrated adjustable dietary and taxon abundance constraints to generate personalized metabolic models for individual metagenomic samples. We applied MICOM to a balanced cohort of metagenomes from 186 people, including a metabolically healthy population and individuals with type 1 and type 2 diabetes. Model results showed that individual bacterial genera maintained conserved niche structures across humans, while the community-level production of short-chain fatty acids (SCFAs) was heterogeneous and highly individual specific. Model output revealed complex cross-feeding interactions that would be difficult to measure in vivo. Metabolic interaction networks differed somewhat consistently between healthy and diabetic subjects. In particular, MICOM predicted reduced butyrate and propionate production in a diabetic cohort, with restoration of SCFA production profiles found in healthy subjects following metformin treatment. Overall, we found that changes in diet or taxon abundances have highly personalized effects. We believe MICOM can serve as a useful tool for generating mechanistic hypotheses for how diet and microbiome composition influence community function. All methods are implemented in an open-source Python package, which is available at https://github.com/micom-dev/micom.

IMPORTANCE The bacterial communities that live within the human gut have been linked to health and disease. However, we are still just beginning to understand how those bacteria interact and what potential interventions to our gut microbiome can make us healthier. Here, we present a mathematical modeling framework (named MICOM) that can recapitulate the growth rates of diverse bacterial species in the gut and can simulate metabolic interactions within microbial communities. We show that MICOM can unravel the ecological rules that shape the microbial landscape in our gut and that a given dietary or probiotic intervention can have widely different effects in different people.

Finding MEMo: minimum sets of elementary flux modes

Metabolic network reconstructions are widely used in computational systems biology for in silico studies of cellular metabolism. A common approach to analyse these models are elementary flux modes (EFMs), which correspond to minimal functional units in the network. Already for medium-sized networks, it is often impossible to compute the set of all EFMs, due to their huge number. From a practical point of view, this might also not be necessary because a subset of EFMs may already be sufficient to answer relevant biological questions. In this article, we study MEMos or minimum sets of EFMs that can generate all possible steady-state behaviours of a metabolic network. The number of EFMs in a MEMo may be by several orders of magnitude smaller than the total number of EFMs. Using MEMos, we can compute generating sets of EFMs in metabolic networks where the whole set of EFMs is too large to be enumerated.

Multidisciplinary approaches for studying rhizobium–legume symbioses

The rhizobium–legume symbiosis is a major source of fixed nitrogen (ammonia) in the biosphere. The potential for this process to increase agricultural yield while reducing the reliance on nitrogen-based fertilizers has generated interest in understanding and manipulating this process. For decades, rhizobium research has benefited from the use of leading techniques from a very broad set of fields, including population genetics, molecular genetics, genomics, and systems biology. In this review, we summarize many of the research strategies that have been employed in the study of rhizobia and the unique knowledge gained from these diverse tools, with a focus on genome- and systems-level approaches. We then describe ongoing synthetic biology approaches aimed at improving existing symbioses or engineering completely new symbiotic interactions. The review concludes with our perspective of the future directions and challenges of the field, with an emphasis on how the application of a multidisciplinary approach and the development of new methods will be necessary to ensure successful biotechnological manipulation of the symbiosis.

Interaction and Regulation of Carbon, Nitrogen, and Phosphorus Metabolisms in Root Nodules of Legumes

Members of the plant family Leguminosae (Fabaceae) are unique in that they have evolved a symbiotic relationship with rhizobia (a group of soil bacteria that can fix atmospheric nitrogen). Rhizobia infect and form root nodules on their specific host plants before differentiating into bacteroids, the symbiotic form of rhizobia. This complex relationship involves the supply of C₄-dicarboxylate and phosphate by the host plants to the microsymbionts that utilize them in the energy-intensive process of fixing atmospheric nitrogen into ammonium, which is in turn made available to the host plants as a source of nitrogen, a macronutrient for growth. Although nitrogen-fixing bacteroids are no longer growing, they are metabolically active. The symbiotic process is complex and tightly regulated by both the host plants and the bacteroids. The metabolic pathways of carbon, nitrogen, and phosphate are heavily regulated in the host plants, as they need to strike a fine balance between satisfying their own needs as well as those of the microsymbionts. A network of transporters for the various metabolites are responsible for the trafficking of these essential molecules between the two partners through the symbiosome membrane (plant-derived membrane surrounding the bacteroid), and these are in turn regulated by various transcription factors that control their expressions under different environmental conditions. Understanding this complex process of symbiotic nitrogen fixation is vital in promoting sustainable agriculture and enhancing soil fertility.

The Role of Soil Microorganisms in Plant Mineral Nutrition-Current Knowledge and Future Directions

In their natural environment, plants are part of a rich ecosystem including numerous and diverse microorganisms in the soil. It has been long recognized that some of these microbes, such as mycorrhizal fungi or nitrogen fixing symbiotic bacteria, play important roles in plant performance by improving mineral nutrition. However, the full range of microbes associated with plants and their potential to replace synthetic agricultural inputs has only recently started to be uncovered. In the last few years, a great progress has been made in the knowledge on composition of rhizospheric microbiomes and their dynamics. There is clear evidence that plants shape microbiome structures, most probably by root exudates, and also that bacteria have developed various adaptations to thrive in the rhizospheric niche. The mechanisms of these interactions and the processes driving the alterations in microbiomes are, however, largely unknown. In this review, we focus on the interaction of plants and root associated bacteria enhancing plant mineral nutrition, summarizing the current knowledge in several research fields that can converge to improve our understanding of the molecular mechanisms underpinning this phenomenon.

A Comprehensively Curated Genome-Scale Two-Cell Model for the Heterocystous Cyanobacterium Anabaena sp. PCC 7120

Anabaena sp. PCC 7120 is a nitrogen-fixing filamentous cyanobacterium. Under nitrogen-limiting conditions, a fraction of the vegetative cells in each filament terminally differentiate to nongrowing heterocysts. Heterocysts are metabolically and structurally specialized to enable O₂-sensitive nitrogen fixation. The functionality of the filament, as an association of vegetative cells and heterocysts, is postulated to depend on metabolic exchange of electrons, carbon, and fixed nitrogen. In this study, we compile and evaluate a comprehensive curated stoichiometric model of this two-cell system, with the objective function based on the growth of the filament under diazotrophic conditions. The predicted growth rate under nitrogen-replete and -deplete conditions, as well as the effect of external carbon and nitrogen sources, was thereafter verified. Furthermore, the model was utilized to comprehensively evaluate the optimality of putative metabolic exchange reactions between heterocysts and vegetative cells. The model suggested that optimal growth requires at least four exchange metabolites. Several combinations of exchange metabolites resulted in predicted growth rates that are higher than growth rates achieved by only considering exchange of metabolites previously suggested in the literature. The curated model of the metabolic network of Anabaena sp. PCC 7120 enhances our ability to understand the metabolic organization of multicellular cyanobacteria and provides a platform for further study and engineering of their metabolism.

Host-Microbiome Interaction and Cancer: Potential Application in Precision Medicine

It has been experimentally shown that host-microbial interaction plays a major role in shaping the wellness or disease of the human body. Microorganisms coexisting in human tissues provide a variety of benefits that contribute to proper functional activity in the host through the modulation of fundamental processes such as signal transduction, immunity and metabolism. The unbalance of this microbial profile, or dysbiosis, has been correlated with the genesis and evolution of complex diseases such as cancer. Although this latter disease has been thoroughly studied using different high-throughput (HT) technologies, its heterogeneous nature makes its understanding and proper treatment in patients a remaining challenge in clinical settings. Notably, given the outstanding role of host-microbiome interactions, the ecological interactions with microorganisms have become a new significant aspect in the systems that can contribute to the diagnosis and potential treatment of solid cancers. As a part of expanding precision medicine in the area of cancer research, efforts aimed at effective treatments for various kinds of cancer based on the knowledge of genetics, biology of the disease and host-microbiome interactions might improve the prediction of disease risk and implement potential microbiota-directed therapeutics. In this review, we present the state of the art of sequencing and metabolome technologies, computational methods and schemes in systems biology that have addressed recent breakthroughs of uncovering relationships or associations between microorganisms and cancer. Together, microbiome studies extend the horizon of new personalized treatments against cancer from the perspective of precision medicine through a synergistic strategy integrating clinical knowledge, HT data, bioinformatics, and systems biology.

Understanding Plant Nitrogen Metabolism through Metabolomics and Computational Approaches

A comprehensive understanding of plant metabolism could provide a direct mechanism for improving nitrogen use efficiency (NUE) in crops. One of the major barriers to achieving this outcome is our poor understanding of the complex metabolic networks, physiological factors, and signaling mechanisms that affect NUE in agricultural settings. However, an exciting collection of computational and experimental approaches has begun to elucidate whole-plant nitrogen usage and provides an avenue for connecting nitrogen-related phenotypes to genes. Herein, we describe how metabolomics, computational models of metabolism, and flux balance analysis have been harnessed to advance our understanding of plant nitrogen metabolism. We introduce a model describing the complex flow of nitrogen through crops in a real-world agricultural setting and describe how experimental metabolomics data, such as isotope labeling rates and analyses of nutrient uptake, can be used to refine these models. In summary, the metabolomics/computational approach offers an exciting mechanism for understanding NUE that may ultimately lead to more effective crop management and engineered plants with higher yields.

Genome-scale reconstruction of the metabolic network in Pseudomonas stutzeri A1501

Pseudomonas stutzeri A1501 is an endophytic bacterium capable of nitrogen fixation. This strain has been isolated from the rice rhizosphere and provides the plant with fixed nitrogen and phytohormones. These interesting features encouraged us to study the metabolism of this microorganism at the systems-level. In this work, we present the first genome-scale metabolic model (iPB890) for P. stutzeri, involving 890 genes, 1135 reactions, and 813 metabolites. A combination of automatic and manual approaches was used in the reconstruction process. Briefly, using the metabolic networks of Pseudomonas aeruginosa and Pseudomonas putida as templates, a draft metabolic network of P. stutzeri was reconstructed. Then, the draft network was driven through an iterative and curative process of gap filling. In the next step, the model was evaluated using different experimental data such as specific growth rate, Biolog substrate utilization data and other experimental observations. In most of the evaluation cases, the model was successful in correctly predicting the cellular phenotypes. Thus, we posit that the iPB890 model serves as a suitable platform to explore the metabolism of P. stutzeri.

Metabolic modelling reveals the specialization of secondary replicons for niche adaptation in Sinorhizobium meliloti

The genome of about 10% of bacterial species is divided among two or more large chromosome-sized replicons. The contribution of each replicon to the microbial life cycle (for example, environmental adaptations and/or niche switching) remains unclear. Here we report a genome-scale metabolic model of the legume symbiont Sinorhizobium meliloti that is integrated with carbon utilization data for 1,500 genes with 192 carbon substrates. Growth of S. meliloti is modelled in three ecological niches (bulk soil, rhizosphere and nodule) with a focus on the role of each of its three replicons. We observe clear metabolic differences during growth in the tested ecological niches and an overall reprogramming following niche switching. In silico examination of the inferred fitness of gene deletion mutants suggests that secondary replicons evolved to fulfil a specialized function, particularly host-associated niche adaptation. Thus, genes on secondary replicons might potentially be manipulated to promote or suppress host interactions for biotechnological purposes.

The genome of some bacteria consists of two or more chromosomes or replicons. Here, diCenzo et al. integrate genome-scale metabolic modelling and growth data from a collection of mutants of the plant symbiont Sinorhizobium meliloti to estimate the fitness contribution of each replicon in three environments.

Metabolic Reconstruction of Setaria italica: A Systems Biology Approach for Integrating Tissue-Specific Omics and Pathway Analysis of Bioenergy Grasses

The urgent need for major gains in industrial crops productivity and in biofuel production from bioenergy grasses have reinforced attention on understanding C4 photosynthesis. Systems biology studies of C4 model plants may reveal important features of C4 metabolism. Here we chose foxtail millet (Setaria italica), as a C4 model plant and developed protocols to perform systems biology studies. As part of the systems approach, we have developed and used a genome-scale metabolic reconstruction in combination with the use of multi-omics technologies to gain more insights into the metabolism of S. italica. mRNA, protein, and metabolite abundances, were measured in mature and immature stem/leaf phytomers, and the multi-omics data were integrated into the metabolic reconstruction framework to capture key metabolic features in different developmental stages of the plant. RNA-Seq reads were mapped to the S. italica resulting for 83% coverage of the protein coding genes of S. italica. Besides revealing similarities and differences in central metabolism of mature and immature tissues, transcriptome analysis indicates significant gene expression of two malic enzyme isoforms (NADP- ME and NAD-ME). Although much greater expression levels of NADP-ME genes are observed and confirmed by the correspondent protein abundances in the samples, the expression of multiple genes combined to the significant abundance of metabolites that participates in C4 metabolism of NAD-ME and NADP-ME subtypes suggest that S. italica may use mixed decarboxylation modes of C4 photosynthetic pathways under different plant developmental stages. The overall analysis also indicates different levels of regulation in mature and immature tissues in carbon fixation, glycolysis, TCA cycle, amino acids, fatty acids, lignin, and cellulose syntheses. Altogether, the multi-omics analysis reveals different biological entities and their interrelation and regulation over plant development. With this study, we demonstrated that this systems approach is powerful enough to complement the functional metabolic annotation of bioenergy grasses.

Reconstruction and Use of Microbial Metabolic Networks: the Core Escherichia coli Metabolic Model as an Educational Guide

Biochemical network reconstructions have become popular tools in systems biology. Metabolicnetwork reconstructions are biochemically, genetically, and genomically (BiGG) structured databases of biochemical reactions and metabolites. They contain information such as exact reaction stoichiometry, reaction reversibility, and the relationships between genes, proteins, and reactions. Network reconstructions have been used extensively to study the phenotypic behavior of wild-type and mutant stains under a variety of conditions, linking genotypes with phenotypes. Such phenotypic simulations have allowed for the prediction of growth after genetic manipulations, prediction of growth phenotypes after adaptive evolution, and prediction of essential genes. Additionally, because network reconstructions are organism specific, they can be used to understand differences between organisms of species in a functional context.There are different types of reconstructions representing various types of biological networks (metabolic, regulatory, transcription/translation). This chapter serves as an introduction to metabolic and regulatory network reconstructions and models and gives a complete description of the core Escherichia coli metabolic model. This model can be analyzed in any computational format (such as MATLAB or Mathematica) based on the information given in this chapter. The core E. coli model is a small-scale model that can be used for educational purposes. It is meant to be used by senior undergraduate and first-year graduate students learning about constraint-based modeling and systems biology. This model has enough reactions and pathways to enable interesting and insightful calculations, but it is also simple enough that the results of such calculations can be understoodeasily.

Sphagnum physiology in the context of changing climate: emergent influences of genomics, modelling and host-microbiome interactions on understanding ecosystem function

Peatlands harbour more than one-third of terrestrial carbon leading to the argument that the bryophytes, as major components of peatland ecosystems, store more organic carbon in soils than any other collective plant taxa. Plants of the genus Sphagnum are important components of peatland ecosystems and are potentially vulnerable to changing climatic conditions. However, the response of Sphagnum to rising temperatures, elevated CO2 and shifts in local hydrology have yet to be fully characterized. In this review, we examine Sphagnum biology and ecology and explore the role of this group of keystone species and its associated microbiome in carbon and nitrogen cycling using literature review and model simulations. Several issues are highlighted including the consequences of a variable environment on plant-microbiome interactions, uncertainty associated with CO2 diffusion resistances and the relationship between fixed N and that partitioned to the photosynthetic apparatus. We note that the Sphagnum fallax genome is currently being sequenced and outline potential applications of population-level genomics and corresponding plant photosynthesis and microbial metabolic modelling techniques. We highlight Sphagnum as a model organism to explore ecosystem response to a changing climate and to define the role that Sphagnum can play at the intersection of physiology, genetics and functional genomics.

Steady-state and instationary modeling of proteinogenic and free amino acid isotopomers for flux quantification

Metabolic flux analysis (MFA) is a powerful tool for exploring and quantifying carbon traffic in metabolic networks. Accurate flux quantification requires (1) high-quality isotopomer measurements, usually of biomass components including proteinogenic/free amino acids or central carbon metabolites, and (2) a mathematical model that relates the unknown fluxes to the measured isotopomers. Modeling requires a thorough knowledge of the structure of the underlying metabolic network, often available from many databases, as well as the ability to make reasonable assumptions that will enable simplification of the model. Here we describe a general methodology underlying computer-aided mathematical modeling of a flux-isotopomer relationship and some of the accompanying data-processing steps. One of two modeling strategies will need to be employed, depending on the type of isotope labeling experiment performed. These strategies-steady-state modeling and instationary modeling-have different experimental and computational demands. We discuss the concepts underlying these two types of modeling and demonstrate steady-state modeling in a step-by-step manner. Our methodology should be applicable to most isotope-assisted MFA applications and should serve as a general framework applicable to many realistic metabolic networks with little modification.

Metabolic cancer biology: structural-based analysis of cancer as a metabolic disease, new sights and opportunities for disease treatment

The cancer cell metabolism or the Warburg effect discovery goes back to 1924 when, for the first time Otto Warburg observed, in contrast to the normal cells, cancer cells have different metabolism. With the initiation of high throughput technologies and computational systems biology, cancer cell metabolism renaissances and many attempts were performed to revise the Warburg effect. The development of experimental and analytical tools which generate high-throughput biological data including lots of information could lead to application of computational models in biological discovery and clinical medicine especially for cancer. Due to the recent availability of tissue-specific reconstructed models, new opportunities in studying metabolic alteration in various kinds of cancers open up. Structural approaches at genome-scale levels seem to be suitable for developing diagnostic and prognostic molecular signatures, as well as in identifying new drug targets. In this review, we have considered these recent advances in structural-based analysis of cancer as a metabolic disease view. Two different structural approaches have been described here: topological and constraint-based methods. The ultimate goal of this type of systems analysis is not only the discovery of novel drug targets but also the development of new systems-based therapy strategies.

Comprehensive classification and perspective for modelling photorespiratory metabolism

Biological processes involved in photorespiratory and photosynthetic metabolism operate concurrently and affect the interplay between carbon and nitrogen assimilation reflected in plant growth. Experimental evidence has indicated that photorespiratory metabolism has a wide-ranging influence not only on other principal metabolic pathways but also on a multitude of signalling cascades. Therefore, accurate quantitative models of photorespiration can provide a means for predicting and in silico probing of plant behaviour at various levels of the system. We first present a comprehensive classification of current models of photorespiratory metabolism developed within the existing carbon-centric modelling paradigm. We then offer a perspective for modelling photorespiratory metabolism by considering the coupling of carbon and nitrogen metabolism in the context of compartmentalised, genome-scale metabolic models of C3 plants. In addition, we outline the challenges stemming from the need to consider plant metabolic and signalling pathways in assessing the still controversial role of photorespiration and to confront the devised models with the ever-increasing amounts of high-throughput data.

Network motifs provide signatures that characterize metabolism

Motifs are repeating patterns that determine the local properties of networks. In this work, we characterized all 3-node motifs using enzyme commission numbers of the International Union of Biochemistry and Molecular Biology to show that motif abundance is related to biochemical function. Further, we present a comparative analysis of motif distributions in the metabolic networks of 21 species across six kingdoms of life. We found the distribution of motif abundances to be similar between species, but unique across cellular organelles. Finally, we show that motifs are able to capture inter-species differences in metabolic networks and that molecular differences between some biological species are reflected by the distribution of motif abundances in metabolic networks.

Systems biology of cancer: moving toward the integrative study of the metabolic alterations in cancer cells

One of the main objectives in systems biology is to understand the biological mechanisms that give rise to the phenotype of a microorganism by using high-throughput technologies (HTs) and genome-scale mathematical modeling. The computational modeling of genome-scale metabolic reconstructions is one systemic and quantitative strategy for characterizing the metabolic phenotype associated with human diseases and potentially for designing drugs with optimal clinical effects. The purpose of this short review is to describe how computational modeling, including the specific case of constraint-based modeling, can be used to explore, characterize, and predict the metabolic capacities that distinguish the metabolic phenotype of cancer cell lines. As we show herein, this computational framework is far from a pure theoretical description, and to ensure proper biological interpretation, it is necessary to integrate high-throughput data and generate predictions for later experimental assessment. Hence, genome-scale modeling serves as a platform for the following: (1) the integration of data from HTs, (2) the assessment of how metabolic activity is related to phenotype in cancer cell lines, and (3) the design of new experiments to evaluate the outcomes of the in silico analysis. By combining the functions described above, we show that computational modeling is a useful methodology to construct an integrative, systemic, and quantitative scheme for understanding the metabolic profiles of cancer cell lines, a first step to determine the metabolic mechanism by which cancer cells maintain and support their malignant phenotype in human tissues.

Transport and metabolism in legume-rhizobia symbioses

Symbiotic nitrogen fixation by rhizobia in legume root nodules injects approximately 40 million tonnes of nitrogen into agricultural systems each year. In exchange for reduced nitrogen from the bacteria, the plant provides rhizobia with reduced carbon and all the essential nutrients required for bacterial metabolism. Symbiotic nitrogen fixation requires exquisite integration of plant and bacterial metabolism. Central to this integration are transporters of both the plant and the rhizobia, which transfer elements and compounds across various plant membranes and the two bacterial membranes. Here we review current knowledge of legume and rhizobial transport and metabolism as they relate to symbiotic nitrogen fixation. Although all legume-rhizobia symbioses have many metabolic features in common, there are also interesting differences between them, which show that evolution has solved metabolic problems in different ways to achieve effective symbiosis in different systems.

Strategy diversity stabilizes mutualism through investment cycles, phase polymorphism, and spatial bubbles

There is continuing interest in understanding factors that facilitate the evolution and stability of cooperation within and between species. Such interactions will often involve plasticity in investment behavior, in response to the interacting partner's investments. Our aim here is to investigate the evolution and stability of reciprocal investment behavior in interspecific interactions, a key phenomenon strongly supported by experimental observations. In particular, we present a comprehensive analysis of a continuous reciprocal investment game between mutualists, both in well-mixed and spatially structured populations, and we demonstrate a series of novel mechanisms for maintaining interspecific mutualism. We demonstrate that mutualistic partners invariably follow investment cycles, during which mutualism first increases, before both partners eventually reduce their investments to zero, so that these cycles always conclude with full defection. We show that the key mechanism for stabilizing mutualism is phase polymorphism along the investment cycle. Although mutualistic partners perpetually change their strategies, the community-level distribution of investment levels becomes stationary. In spatially structured populations, the maintenance of polymorphism is further facilitated by dynamic mosaic structures, in which mutualistic partners form expanding and collapsing spatial bubbles or clusters. Additionally, we reveal strategy-diversity thresholds, both for well-mixed and spatially structured mutualistic communities, and discuss factors for meeting these thresholds, and thus maintaining mutualism. Our results demonstrate that interspecific mutualism, when considered as plastic investment behavior, can be unstable, and, in agreement with empirical observations, may involve a polymorphism of investment levels, varying both in space and in time. Identifying the mechanisms maintaining such polymorphism, and hence mutualism in natural communities, provides a significant step towards understanding the coevolution and population dynamics of mutualistic interactions.

Mutualistic interactions between species are often best understood as gradually adjustable reciprocal investments made continuously or iteratively between participants. Prime examples are the mycorrhizal and rhizobial mutualisms so strongly affecting the productivity of plants. When such interactions are described by continuous reciprocal investment games, participants adjust their investments plastically in response to their mutualistic partner's most recent investment. Although common sense suggests that such conditional or reactive behavior provides a potent defense against exploitation, our comprehensive model analysis reveals that the coevolution of investment strategies will often instead induce instability and decay of mutualistic interactions. We also identify several factors that can prevent this decay. First, mutualisms can be stably maintained if the investment strategies of participants are sufficiently diverse. Second, if participants are limited in their movements, the formation of dynamic spatial mosaic structures promotes strategy diversity and thereby facilitates the maintenance of mutualism. These ecological and evolutionary dynamics result in communities with a diversity of interaction types, ranging from mutually beneficial to exploitative, and varying both in space and in time.

Functional modules, structural topology, and optimal activity in metabolic networks

Modular organization in biological networks has been suggested as a natural mechanism by which a cell coordinates its metabolic strategies for evolving and responding to environmental perturbations. To understand how this occurs, there is a need for developing computational schemes that contribute to integration of genomic-scale information and assist investigators in formulating biological hypotheses in a quantitative and systematic fashion. In this work, we combined metabolome data and constraint-based modeling to elucidate the relationships among structural modules, functional organization, and the optimal metabolic phenotype of Rhizobium etli, a bacterium that fixes nitrogen in symbiosis with Phaseolus vulgaris. To experimentally characterize the metabolic phenotype of this microorganism, we obtained the metabolic profile of 220 metabolites at two physiological stages: under free-living conditions, and during nitrogen fixation with P. vulgaris. By integrating these data into a constraint-based model, we built a refined computational platform with the capability to survey the metabolic activity underlying nitrogen fixation in R. etli. Topological analysis of the metabolic reconstruction led us to identify modular structures with functional activities. Consistent with modular activity in metabolism, we found that most of the metabolites experimentally detected in each module simultaneously increased their relative abundances during nitrogen fixation. In this work, we explore the relationships among topology, biological function, and optimal activity in the metabolism of R. etli through an integrative analysis based on modeling and metabolome data. Our findings suggest that the metabolic activity during nitrogen fixation is supported by interacting structural modules that correlate with three functional classifications: nucleic acids, peptides, and lipids. More fundamentally, we supply evidence that such modular organization during functional nitrogen fixation is a robust property under different environmental conditions.

Biological networks are an inherent concept in systems biology that is useful in elucidating how biological entities—as metabolites or proteins—work together in supporting specific phenotypes in microorganisms. Notably, topological analyses carried out over these networks have shown that modular organization is a ubiquitous property at different levels of biological organization, in such a way that modular organization may serve as an organizing principle governing the metabolic activity in microorganisms. With the aim of elucidating the relationship among functional modules, network topology, and optimal metabolic activity, here we present an integrative study that combines computational modeling and metabolome data for evaluation of the metabolic activity of the soil bacterium Rhizobium etli during symbiotic nitrogen fixation with Phaseolus vulgaris. As a result, we supply experimental and computational evidence supporting the concept that the optimal metabolic activity during this biological process is guided by modular structures in the metabolic network of R. etli. Even more fundamentally, we suggest that these biochemical modules interact among each other to ensure an optimal phenotype during nitrogen fixation. Finally, through the in silico analysis on the genome scale metabolic reconstruction for R.etli, we give some examples that suggest that these modular structures supporting nitrogen fixation are robust to external physiological conditions.

Role of trehalose in heat and desiccation tolerance in the soil bacterium Rhizobium etli

BACKGROUND

The compatible solute trehalose is involved in the osmostress response of Rhizobium etli, the microsymbiont of Phaseolus vulgaris. In this work, we reconstructed trehalose metabolism in R. etli, and investigated its role in cellular adaptation and survival to heat and desiccation stress under free living conditions.

RESULTS

Besides trehalose as major compatible solute, R. etli CE3 also accumulated glutamate and, if present in the medium, mannitol. Putative genes for trehalose synthesis (otsAB/treS/treZY), uptake (aglEFGK/thuEFGK) and degradation (thuAB/treC) were scattered among the chromosome and plasmids p42a, p42c, p42e, and p42f, and in some instances found redundant. Two copies of the otsA gene, encoding trehalose-6-P-synthase, were located in the chromosome (otsAch) and plasmid p42a (otsAa), and the latter seemed to be acquired by horizontal transfer. High temperature alone did not influence growth of R. etli, but a combination of high temperature and osmotic stress was more deleterious for growth than osmotic stress alone. Although high temperature induced some trehalose synthesis by R. etli, trehalose biosynthesis was mainly triggered by osmotic stress. However, an otsAch mutant, unable to synthesize trehalose in minimal medium, showed impaired growth at high temperature, suggesting that trehalose plays a role in thermoprotection of R. etli. Desiccation tolerance by R. etli wild type cells was dependent of high trehalose production by osmotic pre-conditioned cells. Cells of the mutant strain otsAch showed ca. 3-fold lower survival levels than the wild type strain after drying, and a null viability after 4 days storage.

CONCLUSIONS

Our findings suggest a beneficial effect of osmotic stress in R. etli tolerance to desiccation, and an important role of trehalose on the response of R. etli to high temperature and desiccation stress.

Structural and dynamical analysis of biological networks

Biological networks are currently being studied with approaches derived from the mathematical and physical sciences. Their structural analysis enables to highlight nodes with special properties that have sometimes been correlated with the biological importance of a gene or a protein. However, biological networks are dynamic both on the evolutionary time-scale, and on the much shorter time-scale of physiological processes. There is therefore no unique network for a given cellular process, but potentially many realizations, each with different properties as a consequence of regulatory mechanisms. Such realizations provide snapshots of a same network in different conditions, enabling the study of condition-dependent structural properties. True dynamical analysis can be obtained through detailed mathematical modeling techniques that are not easily scalable to full network models.

In silico insights into the symbiotic nitrogen fixation in Sinorhizobium meliloti via metabolic reconstruction

BACKGROUND

Sinorhizobium meliloti is a soil bacterium, known for its capability to establish symbiotic nitrogen fixation (SNF) with leguminous plants such as alfalfa. S. meliloti 1021 is the most extensively studied strain to understand the mechanism of SNF and further to study the legume-microbe interaction. In order to provide insight into the metabolic characteristics underlying the SNF mechanism of S. meliloti 1021, there is an increasing demand to reconstruct a metabolic network for the stage of SNF in S. meliloti 1021.

RESULTS

Through an iterative reconstruction process, a metabolic network during the stage of SNF in S. meliloti 1021 was presented, named as iHZ565, which accounts for 565 genes, 503 internal reactions, and 522 metabolites. Subjected to a novelly defined objective function, the in silico predicted flux distribution was highly consistent with the in vivo evidences reported previously, which proves the robustness of the model. Based on the model, refinement of genome annotation of S. meliloti 1021 was performed and 15 genes were re-annotated properly. There were 19.8% (112) of the 565 metabolic genes included in iHZ565 predicted to be essential for efficient SNF in bacteroids under the in silico microaerobic and nutrient sharing condition.

CONCLUSIONS

As the first metabolic network during the stage of SNF in S. meliloti 1021, the manually curated model iHZ565 provides an overview of the major metabolic properties of the SNF bioprocess in S. meliloti 1021. The predicted SNF-required essential genes will facilitate understanding of the key functions in SNF and help identify key genes and design experiments for further validation. The model iHZ565 can be used as a knowledge-based framework for better understanding the symbiotic relationship between rhizobia and legumes, ultimately, uncovering the mechanism of nitrogen fixation in bacteroids and providing new strategies to efficiently improve biological nitrogen fixation.

MetRxn: a knowledgebase of metabolites and reactions spanning metabolic models and databases

BACKGROUND

Increasingly, metabolite and reaction information is organized in the form of genome-scale metabolic reconstructions that describe the reaction stoichiometry, directionality, and gene to protein to reaction associations. A key bottleneck in the pace of reconstruction of new, high-quality metabolic models is the inability to directly make use of metabolite/reaction information from biological databases or other models due to incompatibilities in content representation (i.e., metabolites with multiple names across databases and models), stoichiometric errors such as elemental or charge imbalances, and incomplete atomistic detail (e.g., use of generic R-group or non-explicit specification of stereo-specificity).

DESCRIPTION

MetRxn is a knowledgebase that includes standardized metabolite and reaction descriptions by integrating information from BRENDA, KEGG, MetaCyc, Reactome.org and 44 metabolic models into a single unified data set. All metabolite entries have matched synonyms, resolved protonation states, and are linked to unique structures. All reaction entries are elementally and charge balanced. This is accomplished through the use of a workflow of lexicographic, phonetic, and structural comparison algorithms. MetRxn allows for the download of standardized versions of existing genome-scale metabolic models and the use of metabolic information for the rapid reconstruction of new ones.

CONCLUSIONS

The standardization in description allows for the direct comparison of the metabolite and reaction content between metabolic models and databases and the exhaustive prospecting of pathways for biotechnological production. This ever-growing dataset currently consists of over 76,000 metabolites participating in more than 72,000 reactions (including unresolved entries). MetRxn is hosted on a web-based platform that uses relational database models (MySQL).

Systems biology of bacterial nitrogen fixation: high-throughput technology and its integrative description with constraint-based modeling

Background

Bacterial nitrogen fixation is the biological process by which atmospheric nitrogen is uptaken by bacteroids located in plant root nodules and converted into ammonium through the enzymatic activity of nitrogenase. In practice, this biological process serves as a natural form of fertilization and its optimization has significant implications in sustainable agricultural programs. Currently, the advent of high-throughput technology supplies with valuable data that contribute to understanding the metabolic activity during bacterial nitrogen fixation. This undertaking is not trivial, and the development of computational methods useful in accomplishing an integrative, descriptive and predictive framework is a crucial issue to decoding the principles that regulated the metabolic activity of this biological process.

Results

In this work we present a systems biology description of the metabolic activity in bacterial nitrogen fixation. This was accomplished by an integrative analysis involving high-throughput data and constraint-based modeling to characterize the metabolic activity in Rhizobium etli bacteroids located at the root nodules of Phaseolus vulgaris (bean plant). Proteome and transcriptome technologies led us to identify 415 proteins and 689 up-regulated genes that orchestrate this biological process. Taking into account these data, we: 1) extended the metabolic reconstruction reported for R. etli; 2) simulated the metabolic activity during symbiotic nitrogen fixation; and 3) evaluated the in silico results in terms of bacteria phenotype. Notably, constraint-based modeling simulated nitrogen fixation activity in such a way that 76.83% of the enzymes and 69.48% of the genes were experimentally justified. Finally, to further assess the predictive scope of the computational model, gene deletion analysis was carried out on nine metabolic enzymes. Our model concluded that an altered metabolic activity on these enzymes induced different effects in nitrogen fixation, all of these in qualitative agreement with observations made in R. etli and other Rhizobiaceas.

Conclusions

In this work we present a genome scale study of the metabolic activity in bacterial nitrogen fixation. This approach leads us to construct a computational model that serves as a guide for 1) integrating high-throughput data, 2) describing and predicting metabolic activity, and 3) designing experiments to explore the genotype-phenotype relationship in bacterial nitrogen fixation.

HepatoNet1: a comprehensive metabolic reconstruction of the human hepatocyte for the analysis of liver physiology

We present HepatoNet1, the first reconstruction of a comprehensive metabolic network of the human hepatocyte that is shown to accomplish a large canon of known metabolic liver functions. The network comprises 777 metabolites in six intracellular and two extracellular compartments and 2539 reactions, including 1466 transport reactions. It is based on the manual evaluation of >1500 original scientific research publications to warrant a high-quality evidence-based model. The final network is the result of an iterative process of data compilation and rigorous computational testing of network functionality by means of constraint-based modeling techniques. Taking the hepatic detoxification of ammonia as an example, we show how the availability of nutrients and oxygen may modulate the interplay of various metabolic pathways to allow an efficient response of the liver to perturbations of the homeostasis of blood compounds.

Modeling core metabolism in cancer cells: surveying the topology underlying the Warburg effect

BACKGROUND

Alterations on glucose consumption and biosynthetic activity of amino acids, lipids and nucleotides are metabolic changes for sustaining cell proliferation in cancer cells. Irrevocable evidence of this fact is the Warburg effect which establishes that cancer cells prefers glycolysis over oxidative phosphorylation to generate ATP. Regulatory action over metabolic enzymes has opened a new window for designing more effective anti-cancer treatments. This enterprise is not trivial and the development of computational models that contribute to identifying potential enzymes for breaking the robustness of cancer cells is a priority.

METHODOLOGY/PRINCIPAL FINDINGS

This work presents a constraint-base modeling of the most experimentally studied metabolic pathways supporting cancer cells: glycolysis, TCA cycle, pentose phosphate, glutaminolysis and oxidative phosphorylation. To evaluate its predictive capacities, a growth kinetics study for Hela cell lines was accomplished and qualitatively compared with in silico predictions. Furthermore, based on pure computational criteria, we concluded that a set of enzymes (such as lactate dehydrogenase and pyruvate dehydrogenase) perform a pivotal role in cancer cell growth, findings supported by an experimental counterpart.

CONCLUSIONS/SIGNIFICANCE

Alterations on metabolic activity are crucial to initiate and sustain cancer phenotype. In this work, we analyzed the phenotype capacities emerged from a constructed metabolic network conformed by the most experimentally studied pathways sustaining cancer cell growth. Remarkably, in silico model was able to resemble the physiological conditions in cancer cells and successfully identified some enzymes currently studied by its therapeutic effect. Overall, we supplied evidence that constraint-based modeling constitutes a promising computational platform to: 1) integrate high throughput technology and establish a crosstalk between experimental validation and in silico prediction in cancer cell phenotype; 2) explore the fundamental metabolic mechanism that confers robustness in cancer; and 3) suggest new metabolic targets for anticancer treatments. All these issues being central to explore cancer cell metabolism from a systems biology perspective.

Integration of metabolic databases for the reconstruction of genome-scale metabolic networks

BACKGROUND

Genome-scale metabolic reconstructions have been recognised as a valuable tool for a variety of applications ranging from metabolic engineering to evolutionary studies. However, the reconstruction of such networks remains an arduous process requiring a high level of human intervention. This process is further complicated by occurrences of missing or conflicting information and the absence of common annotation standards between different data sources.

RESULTS

In this article, we report a semi-automated methodology aimed at streamlining the process of metabolic network reconstruction by enabling the integration of different genome-wide databases of metabolic reactions. We present results obtained by applying this methodology to the metabolic network of the plant Arabidopsis thaliana. A systematic comparison of compounds and reactions between two genome-wide databases allowed us to obtain a high-quality core consensus reconstruction, which was validated for stoichiometric consistency. A lower level of consensus led to a larger reconstruction, which has a lower quality standard but provides a baseline for further manual curation.

CONCLUSION

This semi-automated methodology may be applied to other organisms and help to streamline the process of genome-scale network reconstruction in order to accelerate the transfer of such models to applications.

Mathematical modeling: bridging the gap between concept and realization in synthetic biology

Mathematical modeling plays an important and often indispensable role in synthetic biology because it serves as a crucial link between the concept and realization of a biological circuit. We review mathematical modeling concepts and methodologies as relevant to synthetic biology, including assumptions that underlie a model, types of modeling frameworks (deterministic and stochastic), and the importance of parameter estimation and optimization in modeling. Additionally we expound mathematical techniques used to analyze a model such as sensitivity analysis and bifurcation analysis, which enable the identification of the conditions that cause a synthetic circuit to behave in a desired manner. We also discuss the role of modeling in phenotype analysis such as metabolic and transcription network analysis and point out some available modeling standards and software. Following this, we present three case studies—a metabolic oscillator, a synthetic counter, and a bottom-up gene regulatory network—which have incorporated mathematical modeling as a central component of synthetic circuit design.