Project-based learning course on metabolic network modelling in computational systems biology.

Project-based learning (PBL) is a dynamic student-centred teaching method that encourages students to solve real-life problems while fostering engagement and critical thinking. Here, we report on a PBL course on metabolic network modelling that has been running for several years within the Master in Integrated Systems Biology (MISB) at the University of Luxembourg. This 2-week full-time block course comprises an introduction into the core concepts and methods of constraint-based modelling (CBM), applied to toy models and large-scale networks alongside the preparation of individual student projects in week 1 and, in week 2, the presentation and execution of these projects. We describe in detail the schedule and content of the course, exemplary student projects, and reflect on outcomes and lessons learned. PBL requires the full engagement of students and teachers and gives a rewarding teaching experience. The presented course can serve as a role model and inspiration for other similar courses.

Towards improved genome-scale metabolic network reconstructions: unification, transcript specificity and beyond.

Genome-scale metabolic network reconstructions provide a basis for the investigation of the metabolic properties of an organism. There are reconstructions available for multiple organisms, from prokaryotes to higher organisms and methods for the analysis of a reconstruction. One example is the use of flux balance analysis to improve the yields of a target chemical, which has been applied successfully. However, comparison of results between existing reconstructions and models presents a challenge because of the heterogeneity of the available reconstructions, for example, of standards for presenting gene-protein-reaction associations, nomenclature of metabolites and reactions or selection of protonation states. The lack of comparability for gene identifiers or model-specific reactions without annotated evidence often leads to the creation of a new model from scratch, as data cannot be properly matched otherwise. In this contribution, we propose to improve the predictive power of metabolic models by switching from gene-protein-reaction associations to transcript-isoform-reaction associations, thus taking advantage of the improvement of precision in gene expression measurements. To achieve this precision, we discuss available databases that can be used to retrieve this type of information and point at issues that can arise from their neglect. Further, we stress issues that arise from non-standardized building pipelines, like inconsistencies in protonation states. In addition, problems arising from the use of non-specific cofactors, e.g. artificial futile cycles, are discussed, and finally efforts of the metabolic modelling community to unify model reconstructions are highlighted.

IDARE2-Simultaneous Visualisation of Multiomics Data in Cytoscape.

Visual integration of experimental data in metabolic networks is an important step to understanding their meaning. As genome-scale metabolic networks reach several thousand reactions, the task becomes more difficult and less revealing. While databases like KEGG and BioCyc provide curated pathways that allow a navigation of the metabolic landscape of an organism, it is rather laborious to map data directly onto those pathways. There are programs available using these kind of databases as a source for visualization; however, these programs are then restricted to the pathways available in the database. Here, we present IDARE2 a cytoscape plugin that allows the visualization of multiomics data in cytoscape in a user-friendly way. It further provides tools to disentangle highly connected network structures based on common properties of nodes and retains structural links between the generated subnetworks, offering a straightforward way to traverse the splitted network. The tool is extensible, allowing the implementation of specialised representations and data format parsers. We present the automated reproduction of the original IDARE nodes using our tool and show examples of other data being mapped on a network of E. coli. The extensibility is demonstrated with two plugins that are available on github. IDARE2 provides an intuitive way to visualise data from multiple sources and allows one to disentangle the often complex network structure in large networks using predefined properties of the network nodes.

A dynamic multi-tissue model to study human metabolism.

Metabolic modeling enables the study of human metabolism in healthy and in diseased conditions, e.g., the prediction of new drug targets and biomarkers for metabolic diseases. To accurately describe blood and urine metabolite dynamics, the integration of multiple metabolically active tissues is necessary. We developed a dynamic multi-tissue model, which recapitulates key properties of human metabolism at the molecular and physiological level based on the integration of transcriptomics data. It enables the simulation of the dynamics of intra-cellular and extra-cellular metabolites at the genome scale. The predictive capacity of the model is shown through the accurate simulation of different healthy conditions (i.e., during fasting, while consuming meals or during exercise), and the prediction of biomarkers for a set of Inborn Errors of Metabolism with a precision of 83%. This novel approach is useful to prioritize new biomarkers for many metabolic diseases, as well as for the integration of various types of personal omics data, towards the personalized analysis of blood and urine metabolites.

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.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0.

Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.

The intertwined metabolism during symbiotic nitrogen fixation elucidated by metabolic modelling.

Genome-scale metabolic network models can be used for various analyses including the prediction of metabolic responses to changes in the environment. Legumes are well known for their rhizobial symbiosis that introduces nitrogen into the global nutrient cycle. Here, we describe a fully compartmentalised, mass and charge-balanced, genome-scale model of the clover Medicago truncatula, which has been adopted as a model organism for legumes. We employed flux balance analysis to demonstrate that the network is capable of producing biomass components in experimentally observed proportions, during day and night. By connecting the plant model to a model of its rhizobial symbiont, Sinorhizobium meliloti, we were able to investigate the effects of the symbiosis on metabolic fluxes and plant growth and could demonstrate how oxygen availability influences metabolic exchanges between plant and symbiont, thus elucidating potential benefits of inter organism amino acid cycling. We thus provide a modelling framework, in which the interlinked metabolism of plants and nodules can be studied from a theoretical perspective.

Constraint Based Modeling Going Multicellular.

Constraint based modeling has seen applications in many microorganisms. For example, there are now established methods to determine potential genetic modifications and external interventions to increase the efficiency of microbial strains in chemical production pipelines. In addition, multiple models of multicellular organisms have been created including plants and humans. While initially the focus here was on modeling individual cell types of the multicellular organism, this focus recently started to switch. Models of microbial communities, as well as multi-tissue models of higher organisms have been constructed. These models thereby can include different parts of a plant, like root, stem, or different tissue types in the same organ. Such models can elucidate details of the interplay between symbiotic organisms, as well as the concerted efforts of multiple tissues and can be applied to analyse the effects of drugs or mutations on a more systemic level. In this review we give an overview of the recent development of multi-tissue models using constraint based techniques and the methods employed when investigating these models. We further highlight advances in combining constraint based models with dynamic and regulatory information and give an overview of these types of hybrid or multi-level approaches.

Benchmarking Procedures for High-Throughput Context Specific Reconstruction Algorithms.

Recent progress in high-throughput data acquisition has shifted the focus from data generation to processing and understanding of how to integrate collected information. Context specific reconstruction based on generic genome scale models like ReconX or HMR has the potential to become a diagnostic and treatment tool tailored to the analysis of specific individuals. The respective computational algorithms require a high level of predictive power, robustness and sensitivity. Although multiple context specific reconstruction algorithms were published in the last 10 years, only a fraction of them is suitable for model building based on human high-throughput data. Beside other reasons, this might be due to problems arising from the limitation to only one metabolic target function or arbitrary thresholding. This review describes and analyses common validation methods used for testing model building algorithms. Two major methods can be distinguished: consistency testing and comparison based testing. The first is concerned with robustness against noise, e.g., missing data due to the impossibility to distinguish between the signal and the background of non-specific binding of probes in a microarray experiment, and whether distinct sets of input expressed genes corresponding to i.e., different tissues yield distinct models. The latter covers methods comparing sets of functionalities, comparison with existing networks or additional databases. We test those methods on several available algorithms and deduce properties of these algorithms that can be compared with future developments. The set of tests performed, can therefore serve as a benchmarking procedure for future algorithms.

Insight into metabolic pathways of the potential biofuel producer, Paenibacillus polymyxa ICGEB2008.

BACKGROUND: Paenibacillus polymyxa is a facultative anaerobe known for production of hydrolytic enzymes and various important biofuel molecules. Despite its wide industrial use and the availability of its genome sequence, very little is known about metabolic pathways operative in the Paenibacillus system. Here, we report metabolic insights of an insect gut symbiont, Paenibacillus polymyxa ICGEB2008, and reveal pathways playing an important role in the production of 2,3-butanediol and ethanol. RESULT: We developed a metabolic network model of P. polymyxa ICGEB2008 with 133 metabolites and 158 reactions. Flux balance analysis was employed to investigate the importance of redox balance in ICGEB2008. This led to the detection of the Bifid shunt, a pathway previously not described in Paenibacillus, which can uncouple the production of ATP from the generation of reducing equivalents. Using a combined experimental and modeling approach, we further studied pathways involved in 2,3-butanediol and ethanol production and also demonstrated the production of hydrogen by the organism. We could further show that the nitrogen source is critical for metabolite production by Paenibacillus, and correctly quantify the influence on the by-product metabolite profile of ICGEB2008. Both simulations and experiments showed that metabolic flux is diverted from ethanol to acetate production when an oxidized nitrogen source is utilized. CONCLUSION: We have created a predictive model of the central carbon metabolism of P. polymyxa ICGEB2008 and could show the presence of the Bifid shunt and explain its role in ICGEB2008. An in-depth study has been performed to understand the metabolic pathways involved in ethanol, 2,3-butanediol and hydrogen production, which can be utilized as a basis for further metabolic engineering efforts to improve the efficiency of biofuel production by this P. polymyxa strain.

Systems approaches to modelling pathways and networks.

It has become commonly accepted that systems approaches to biology are of outstanding importance to gain understanding from the vast amount of data which is presently being generated by advancing high-throughput technologies. The diversity of methods to model pathways and networks has significantly expanded over the past two decades. Modern and traditional approaches are equally important and recent activities aim at integrating the advantages of both. While traditional methods, based on differential equations, are useful to study the dynamics of small systems, modern constraint-based models can be applied to genome-scale systems, but are not able to capture dynamic features. Integrating different approaches is important to develop consistent theoretical descriptions encompassing various scales of biological information. The rapid progress of the field of theoretical systems biology, however, demonstrates how our fundamental theoretical understanding of biology is gaining momentum. The scientific community has apparently accepted the challenge to truly understand the principles of life.