Visualizing metabolic network dynamics through time-series metabolomic data
BMC Bioinformatics volume 21, Article number: 130 (2020)
New technologies have given rise to an abundance of -omics data, particularly metabolomic data. The scale of these data introduces new challenges for the interpretation and extraction of knowledge, requiring the development of innovative computational visualization methodologies. Here, we present GEM-Vis, an original method for the visualization of time-course metabolomic data within the context of metabolic network maps. We demonstrate the utility of the GEM-Vis method by examining previously published data for two cellular systems—the human platelet and erythrocyte under cold storage for use in transfusion medicine.
The results comprise two animated videos that allow for new insights into the metabolic state of both cell types. In the case study of the platelet metabolome during storage, the new visualization technique elucidates a nicotinamide accumulation that mirrors that of hypoxanthine and might, therefore, reflect similar pathway usage. This visual analysis provides a possible explanation for why the salvage reactions in purine metabolism exhibit lower activity during the first few days of the storage period. The second case study displays drastic changes in specific erythrocyte metabolite pools at different times during storage at different temperatures.
The new visualization technique GEM-Vis introduced in this article constitutes a well-suitable approach for large-scale network exploration and advances hypothesis generation. This method can be applied to any system with data and a metabolic map to promote visualization and understand physiology at the network level. More broadly, we hope that our approach will provide the blueprints for new visualizations of other longitudinal -omics data types. The supplement includes a comprehensive user’s guide and links to a series of tutorial videos that explain how to prepare model and data files, and how to use the software SBMLsimulator in combination with further tools to create similar animations as highlighted in the case studies.
Over the last few decades, new technological developments have enabled the generation of vast amounts of “-omics” data . These various -omic data types have helped bring new insights to a vast array of biological questions [2–4]. As more and more data are generated, however, researchers are faced with the enormous challenge of integrating, interpreting, and visualizing these data. The community has recognized these needs, focusing efforts on data visualization as a way to maximize the utility of biological data . Data visualization is particularly crucial for a systems-level perspective of metabolic networks and pathways. Several excellent software tools were made available for drawing and exploring biological network graphs [6–10]. These tools provide impressive descriptions of the network and support for diverse analyses, including the mapping of omics data to networks. In this study, we present GEM-Vis as a new approach for the visualization of time-course metabolomic data in the context of large-scale metabolic network maps.
Metabolomic data provide snapshots of cellular biochemistry, presenting essential insights into a cell’s metabolic state [11, 12]. Visualization tools often allow users to overlay pathway maps with static data sets . Recently, time-course metabolomic data sets that detail cellular changes over time are becoming more prevalent [13–16], leading to the need for dynamic visualizations that can capture the aspect of time —an essential aspect of understanding complex processes such as changes in metabolic activity, concentration, or availability. Many visualization tools [18–21], however, do not yet provide support for the representation of dynamic content. Those visualization tools whose features do include time series visualization [5, 17, 20, 22–24] only provide static depictions of the data. Some progress has been made to provide a stepwise temporal representation of metabolomic data , but a robust and smooth dynamic solution for mapping time series data to networks has yet to be presented.
One reason for the current lack of convincing visual analysis methods for dynamically changing data sets is that time-dependent data add additional layers of complexity to the already difficult problem of visual network exploration. First of all, genome-scale metabolic networks (GEMs) can have enormous sizes: Some published metabolic network maps comprise several thousand biochemical reactions [26, 27], of which human beholders can simultaneously only grasp a very small fraction .
With the development of new experimental technologies and the subsequent generation of -omics data sets, life scientists are faced with the challenge of extracting actionable knowledge. New visualization methods are a critical way that the community can make strides toward making the most of complex data. Here, we present a new method for the visualization of longitudinal metabolomic data in the context of the metabolic network. We provide two case studies that examine (1) a baseline characterization of a physiological process and (2) a set of experimental perturbations that allowed for a side-by-side comparison of different experimental conditions. The introduction of this new visualization method has two significant implications.
The method introduced in this article provides a dynamic visualization of cellular processes. Tools such as Cytoscape  provide visual analysis of networks and supports plugins like TiCoNE  and CyDataSeries  for the visualization of time-course data. However, tools such as these or VANTED  only offer static representations of dynamic data. To our knowledge, only KEGGanim  offers a dynamic visualization of time-course data. The method presented here builds on KEGGanim by offering a smooth interpolation between time points and offers the further advantage of customization concerning the display of both data and the network itself. The method presented outlines an original development for visualizing complex biological data based in a way that a cognition study has found to be useful and support perception .
With a steadily increasing number of carefully prepared metabolic network layouts being published, we here assume a map to be available for the system of interest. If this is not yet the case, a map can be easily drawn using software such as Escher . This paper focuses on the problem of displaying dynamically changing quantitative data of network components. The aim is to answer the question: How to create expressive visual displays of dynamic metabolic networks? Needed are strategies to visually present the data in a way that beholders can best perceive and estimate quantities of network individual components and that at the same time enable them to conceptually narrow down parts of interest even within large networks.
In the next sections, we present the new method GEM-Vis (Genome-Scale Metabolic model Visualization) for the visualization and contextualization of longitudinal metabolomic data in metabolic networks. We developed three different graphical representations of metabolic concentration that allow for different interpretations of metabolomic data through a smooth animation. The method is implemented in the freely available software SBMLsimulator. The supplementary material of this article includes links to a series of four short tutorial videos that explain all aspects needed for creating a GEM-Vis: (1) where to obtain SBMLsimulator, how to run it, (2) where to obtain systems biology models and how to load them into the application, and how to create a simulated time-course data set from the model that can be mapped to an automatically generated pathway map, (3) where to obtain a manually drawn pathway map as well as published time-course data of the yeast Saccharomyces cerevisiae and how to prepare the data set for the import to SBMLsimulator and how to embed this layout information  in an SBML file , (4) how to load model and data into SBMLsimulator to create a GEM-Vis animation video including variation of several visual attribues and to save it to a movie file. Finally, we present two case studies using this method that examine two different cellular systems—the human platelet and the human red blood cell (RBC)—to show how visualizing existing data can provide new insights into cellular metabolism. The result are two animated videos that give detailed information about the systems under study and highlight new insights that were not previously apparent. We hope that this method will aid researchers in visualizing, perceiving, and interpreting complex data sets.
The idea of GEM-Vis is that time series can be adequately observed in the form of an animated sequence of a dynamically changing network map when using an appropriate representation of metabolic quantities. To this end, our technique exploits the repeatedly observed ability of human beholders to estimate quantities most precisely when these are mapped to a lengths scale . Since metabolic maps commonly represent nodes with circles [6, 34], we suggest using the fill level of each node as a visual element to represent its amount at each time point. We experimented with visualization of data in several different ways, based on node size, color, a combination of size and color, or fill level (Supplementary Figure S1). Each of these visual representations provides some advantages over the others, but according to  the notion of the fill level of a node can be the most intuitive as it allows for the user to understand and gauge its minimum or maximum value quickly (see Discussion).
Using this technique, we created such an animation for given longitudinal metabolomic data and a metabolic network map that corresponds to the observed cell type (Fig. 1). To provide a smooth animation, additional time points are interpolated in the provided time series. Further details regarding the development and use of the implementation of the method can be found in the Supplementary Information.
To demonstrate the utility of this method, we applied these visualization methods to four different cellular systems—human hepatocytes , platelets  and RBCs , as well as to yeast . For the two human blood cell types and for yeast longitudinal quantitative data sets were available in the literature [13, 14, 39]. Consequently, all four models provide very different use-case scenarios. Since the hepatocyte model  is a fully-specified kinetic model and available in SBML format from BioModels database , it is well suitable to demonstrate how simulated data can be generated and visualized in the context of an algorithmically generated network (see Additional file 11). The genome-scale model of yeast  can be downloaded in SBML and JSON format from BiGG Models Database , where it comes with a manually drawn network of the organism’s central carbon metabolism. It is, therefore, usable to demonstrate mapping a published time-course metabolomic data set  in the context of a hand-made pathway map (see Additional file 13).
After gaining experience in working with the visualization method, the focus will be applying the GEM-Vis method to study human blood cells in more detail. Transfusion medicine plays a vital role in modern healthcare, making the storage of different blood components important physiological processes to understand. In particular, platelets and RBCs represent relatively simple human cell types that can be intensely studied in the well-defined, static environment provided by blood storage (packed in plastic bags and stored at 22∘C and 4∘C for platelets and RBCs, respectively). While the cells are stored in these conditions, biochemical and morphological changes occur (the “storage lesion”) that are well-studied through the use of metabolomic data [12, 42]. Metabolic models were previously available for both the platelet  and RBC , enabling the creation of network maps for both reconstructions. Thus, these data could be visualized in the context of the entire metabolic network.
Case study: human platelets under storage conditions Our first case study examined the storage of platelets. We manually created a metabolic map for the complete metabolic network of the platelet using Escher . We then overlaid metabolomic data which characterized the baseline storage conditions with eight time points over ten days of storage  to produce a network-level visualization of the data (Fig. 2). Using this network-level visualization, we examined the dynamics of the platelet metabolome.
During the first part of storage, stress due to the non-physiological conditions of storage (i.e., packed in a plastic bag at 22∘C) slows metabolic activity through glycolysis, the pentose phosphate pathway, and purine salvage pathways . Several metabolites are secreted by the cells and accumulate in the storage media, such as hypoxanthine. The metabolite 5-Phospho- α-D-ribose 1-diphosphate (PRPP) is produced from the pentose phosphate pathway and is a cofactor in the salvage reactions that break down hypoxanthine. Because flux through the pentose phosphate pathway is lower, the cells have less capacity to recycle hypoxanthine using the salvage pathways.
By viewing all of the data simultaneously at the network level, we were able to discover that the concentration profile of nicotinamide mirrors that of hypoxanthine. This observation suggests a similar rationale for the accumulation of nicotinamide, providing a hypothesis as to why the salvage pathway within purine metabolism has lower activity during the first few days of storage. These findings are demonstrated in the video highlighted in (Table 1), helping show how network-level visualization allows for improved extraction of biological insight from large, complex data sets.
Case study: human red blood cells under storage conditions Our second case study examined the storage of RBCs. A metabolic map was already available for the RBC  and captures the complete metabolic network . Here, we sought to examine a data set that provided the opportunity to visualize different experimental conditions for the same network. Recently, a study was published  that used quantitative longitudinal metabolomic data to examine the state of the RBC metabolome under four different storage temperatures: 4∘C (storage temperature), 13∘C, 22∘C, and 37∘C (body temperature). For this system, we opted to visualize the dynamics of the metabolite concentrations as nodes with variable size where smaller nodes represent smaller pool sizes, and larger nodes represent larger pool sizes (Fig. 3).
To highlight the differences between the experimental conditions, we examined two of the conditions side-by-side (see the video highlighted in Table 2). This visualization helps supplement the type of statistical and modeling analyses performed previously and helps contextualize the effects of the temperature change across different parts of the network. In particular, it is obvious from a network-level view of the system that certain parts of the network are more active at different points in the time-course. A side-by-side comparison helped emphasize that the availability of reduced glutathione is different with increased temperature, an important physiological feature due to the role of glutathione in neutralizing reactive oxygen species  that accumulate during storage and contribute to the storage lesion . Finally, it can be seen that hypoxanthine—a known toxic metabolite whose concentration has been shown to inversely correlate with the post-transfusion recovery rates of transfusion patients —accumulates faster at higher temperatures. Like in the other case study presented above, the new insights into complex processes (which are not yet fully understood) provide evidence that this method can be beneficial for the simplification and understanding of large, complex data analyses.
In this article, we proposed GEM-Vis as a new method for visualizing time-course metabolomic data in the context of large-scale metabolic networks. The approach was evaluated with a range of different use-cases, ranging from the display of simulated data on automatically generated network layouts to experimentally obtained metabolite concentration data on manually drawn network maps. All experiments were described in elaborate tutorial videos (see supplementary material). Subsequently, the method was applied to study two different cases of human blood cells (platelets and erythrocytes) in more detail.
As a result, a network-level representation of large metabolomic data sets presents a more holistic view of the data than does statistical analysis alone. While visual inspection of data is indeed not a replacement for more detailed statistical or modeling analyses, this method provides an important supplement to existing data analysis pipelines. We demonstrate its utility in such an analysis pipeline by highlighting findings from existing data sets [13, 14]. Visualizing the metabolomic data in the context of the full metabolic network allowed for new insights into existing data sets. A potential explanation why the salvage pathway lowers its activity during the first few days of platelet storage could be deduced for the network of the human platelet. In the RBC network, it could easily be seen that concentrations in certain parts of the network (e.g., nucleotide metabolism) accumulated or depleted together. These findings illustrate the promising potential of visualized time-course data and—combined with in-depth computational data analysis—can help perceiving information and elucidate physiological processes.
The simplification of experimental data interpretation became extremely relevant in the age of high-throughput technologies. The visualization concept presented here offers a systems-level interpretation of metabolomic data. Combined with other data analytics, this method helps provide a holistic view of a data set, moving us closer to being able to realize the full potential of a given data set. More broadly, we hope that the method presented here will provide the starting point for further visualization improvements not only for metabolomic data but for the visualization and contextualization of other data types. Future work may include combining a dynamic representation with static concentration graphs that will continue to improve the capabilities of such software to fully meet the needs of life science researchers.
The method described in this paper utilizes existing software libraries to visually represent metabolomic data in the context of a metabolic network map.
In brief, the metabolic map must be embedded as SBML Layout extension (version 1) into an SBML Level 3 Version 1 file that is provided to the software SBMLsimulator . In this study, the design of metabolic network maps was created using the web-based software Escher  within the web browser Safari 11 and stored in JSON format, resulting in Additional file 4 for iAT-PLT-636 and Additional file 8 for iAB-RBC-283. Subsequently, the generated maps have been converted to SBML using the software EscherConverter (available at https://github.com/draeger-lab/EscherConverter/) and embedded into the metabolic model using a custom Java™ program (Additional file 14). The resulting SBML file with embedded layout for iAT-PLT-636 can be found in Additional file 5, and in Additional file 9 for iAB-RBC-283 (both files are compressed using GZIP). The metabolomic time-course data are provided to SBMLsimulator in a *.csv file format with identifiers matching those of the map (Additional file 3 for iAT-PLT-636 and Additional file 7 for iAB-RBC-283, the latter is a compressed ZIP archive). SBMLsimulator reads in the SBML file with embedded layout and the time-course data. Subsequently, SBMLsimulator uses splines to interpolate the data over time with input from the user. Other features are selected, such as the speed of animation and how metabolite concentrations are represented (e.g., fill level). An optional *.csv file can be provided to SBMLsimulator to define a moving camera animation. To this end, this CSV file needs to contain as the first value the zoom level of the camera followed by a tab-separated list of corner points along the way of the moving camera (these points are the top left corners of the camera’s view port). The result is a smooth animation that allows features such as zooming and panning across different areas of the map, which the user can safe to a video file in one of the supported formats, e.g., AVI, MP4, MPG, WMV, FLV, or MOV. The procedure is depicted in Fig. 4 and demonstrated in detail in Additional files 10, 11, 12 and 13.
This tutorial video (36.3 MB) demonstrates how to download, installation, and launch the software SBMLsimulator. The video is available at https://youtu.be/Eu4uSPmNXVI.
This tutorial video (61.8 MB) demonstrates how to load model files in SBML format and how to run a simulation using the software SBMLsimulator. The video is available at https://youtu.be/CVzp_XtIaHU.
This tutorial video (66.7 MB) demonstrates how to embedding a model layout in an SBML file and how to prepare experimental data for loading it into the software SBMLsimulator. The video is available at https://youtu.be/CoeOh2sFFSQ.
This tutorial video (52.9 MB) demonstrates how to visualize manually created layouts using the software SBMLsimulator. The video is available at https://youtu.be/qv3qPyzofhI.
SBMLsimulator is implemented in Java™ SE 8 under macOS High Sierra on a MacBook Pro 15”, 2016. All computation for the animation videos has been performed under macOS High Sierra version 10.13.2. The animation videos for the two case studies were created using Windows 10.
Audio recording was performed using a ZOOM Handy Recorder H4 and a Steinberg UR22 mkII USB Audio Interface 24 bit/192 kHz (Steinberg Media Technologies GmbH, Hamburg, Germany) in combination with a Røde NT1-A (Røde, Silverwater, NSW, Australia), and the recording software Quicktime (Apple Inc., Cuppertino, CA, USA). Sony VEGAS Pro (version 12) was used for video post-processing, resulting in Additional files 2 and 6.
Availability of data and materials
Data and materials required to reproduce the findings in this article are freely available. See the appendix in Additional file 1 for details. For post-processing generated animation files, a commercial video editing software, such as Sony VEGAS Pro (version 12) that was used for creating Additional files 2 and 6, can be obtained for purchase. Computers and audio recording devices need to be purchased.
Audio video interleave
Biochemically, genetically, genomically structured
Genome-scale metabolic model
QuickTime file format
Moving pictures expert group
Multimedia file format
5-phospho- α-d-ribose 1-diphosphate
Red blood cell
Systems biology markup language
Windows media video
Robinson JL, Nielsen J. Integrative analysis of human omics data using biomolecular networks. Mol BioSyst. 2016; 12(10):2953–64. https://doi.org/10.1039/C6MB00476H.
Österlund T, Cvijovic M, Kristiansson E. Integrative Analysis of Omics Data In: Nielsen J, Hohmann S, editors. Systems Biology. Weinheim: Wiley-VCH Verlag GmbH & Co. KGaA: 2017. p. 1–24. Chap. 1. https://doi.org/10.1002/9783527696130.ch1.
Nielsen J. Systems Biology of Metabolism. Ann Rev Biochem. 2017; 86(1):245–75. https://doi.org/10.1146/annurev-biochem-061516-044757.
Yurkovich JT, Palsson BO. Quantitative -omic data empowers bottom-up systems biology. Curr Opin Biotechnol. 2018; 51:130–6. https://doi.org/10.1016/j.copbio.2018.01.009.
Callaway E. The visualizations transforming biology. Nature. 2016; 535:187–8. https://doi.org/10.1038/535187a.
King ZA, Dräger A, Ebrahim A, Sonnenschein N, Lewis NE, Palsson BO. Escher: A web application for building, sharing, and embedding Data-Rich visualizations of biological pathways. PLoS Comput Biol. 2015; 11(8):1004321.
Smoot ME, Ono K, Ruscheinski J, Wang P-L, Ideker T. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011; 27(3):431–2.
Droste P, Nöh K, Wiechert W. Omix – a visualization tool for metabolic networks with highest usability and customizability in focus. Chem Ing Tech. 2013; 85(6):849–62.
Funahashi A, Matsuoka Y, Jouraku A, Morohashi M, Kikuchi N, Kitano H. CellDesigner 3.5: A versatile modeling tool for biochemical networks. Proc IEEE. 2008; 96(8):1254–65.
Kelder T, van Iersel MP, Hanspers K, Kutmon M, Conklin BR, Evelo CT, Pico AR. WikiPathways: building research communities on biological pathways. Nucleic Acids Res. 2012; 40(Database issue):1301–7.
Patti GJ, Yanes O, Siuzdak G. Innovation: Metabolomics: the apogee of the omics trilogy. Nat Rev Mol Cell Biol. 2012; 13:263–9. https://doi.org/10.1038/nrm3314.
Nemkov T, Hansen KC, Dumont LJ, D’Alessandro A. Metabolomics in transfusion medicine. Transfusion. 2016; 56(4):980–93. https://doi.org/10.1111/trf.13442.
Paglia G, Sigurjónsson OE, Rolfsson O, Valgeirsdottir S, Hansen MB, Brynjólfsson S, Gudmundsson S, Palsson BO. Comprehensive metabolomic study of platelets reveals the expression of discrete metabolic phenotypes during storage. Transfusion. 2014; 54:2911–23. https://doi.org/10.1111/trf.12710.
Yurkovich JT, Zielinski DC, Yang L, Paglia G, Rolfsson O, Sigurjónsson OE, Broddrick JT, Bordbar A, Wichuk K, Brynjólfsson S, Palsson S, Gudmundsson S, Palsson BO. Quantitative time-course metabolomics in human red blood cells reveal the temperature dependence of human metabolic networks. J Biol Chem. 2017; 292:117–804914. https://doi.org/10.1074/jbc.M117.804914.
Bordbar A, Yurkovich JT, Paglia G, Rolfsson O, Sigurjónsson OE, Palsson BO. Elucidating dynamic metabolic physiology through network integration of quantitative timecourse metabolomics. Sci Rep. 2017; 7(46249):1–12. https://doi.org/10.1038/srep46249.
Paglia G, Sigurjónsson ÓE, Rolfsson Ó., Hansen MB, Brynjólfsson S, Gudmundsson S, Palsson BO. Metabolomic analysis of platelets during storage: a comparison between apheresis- and buffy coat-derived platelet concentrates. Transfusion. 2015; 55(2):301–13.
Secrier M, Schneider R. Visualizing time-related data in biology, a review. Brief Bioinforma. 2014; 15(5):771–82. https://doi.org/10.1093/bib/bbt021.
Wiwie C, Rauch A, Haakonsson A, Barrio-Hernandez I, Blagoev B, Mandrup S, Röttger R, Baumbach J. Elucidation of time-dependent systems biology cell response patterns with time course network enrichment. 2017. https://arxiv.org/abs/1710.10262.
Černý M. Improve handling of time series and similar in Cytoscape. 2017. https://github.com/nrnb/GoogleSummerOfCode/issues/76. Accessed 8 Apr 2019.
Rohn H, Junker A, Hartmann A, Grafahrend-Belau E, Treutler H, Klapperstück M, Czauderna T, Klukas C, Schreiber F. VANTED v2: a framework for systems biology applications. BMC Syst Biol. 2012; 6(139). https://doi.org/10.1186/1752-0509-6-139.
Ma DKG, Stolte C, Kaur S, Bain M, O’Donoghue SI. Visual analytics of phosphorylation time-series data on insulin response. AIP Conf Proc. 2013; 1559:185–96. https://doi.org/10.1063/1.4825010.
Pavlopoulos GA, Wegener A-L, Schneider R. A survey of visualization tools for biological network analysis. BioData Min. 2008; 1(12). https://doi.org/10.1186/1756-0381-1-12.
König M, Dräger A, Holzhütter H-G. CySBML: a Cytoscape plugin for SBML. Bioinformatics. 2012; 28(18):2402–3. https://doi.org/10.1093/bioinformatics/bts432.
Gehlenborg N, O’Donoghue SI, Baliga NS, Goesmann A, Hibbs H, Kitano MA, Kohlbacher O, Neuweger H, Schneider R, Tenenbaum D, Gavin A-C. Visualization of omics data for systems biology. Nat Methods. 2010; 7(3 Suppl):56–68. https://doi.org/10.1038/nmeth.1436.
Adler P, Reimand J, Jänes J, Kolde R, Peterson H, Vilo J. KEGGanim: pathway animations for high-throughput data. Bioinformatics. 2008; 24(4):588–90. https://doi.org/10.1093/bioinformatics/btm581.
Brunk E, Sahoo S, Zielinski DC, Altunkaya A, Dräger A, Mih N, Gatto F, Nilsson A, Preciat Gonzalez GA, Aurich MK, Prlić A, Sastry A, Danielsdottir AD, Heinken A, Noronha A, Rose PW, Burley SK, Fleming RMT, Nielsen J, Thiele I, Palsson BO. Recon3D enables a three-dimensional view of gene variation in human metabolism. Nat Biotechnol. 2018:1–37. https://doi.org/10.1038/nbt.4072.
Noronha A, Daníelsdóttir AD, Gawron P, Jóhannsson F, Jónsdóttir S, Jarlsson S, Gunnarsson JP, Brynjólfsson S, Schneider R, Thiele I, Fleming RMT. Reconmap: an interactive visualization of human metabolism. Bioinformatics. 2017; 33(4):605–7. https://doi.org/10.1093/bioinformatics/btw667.
Halford GS, Baker R, McCredden JE, Bain JD. How many variables can humans process?,. Psychol Sci. 2005; 16(1):70–76. https://doi.org/10.1111/j.0956-7976.2005.00782.x.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski BB, Ideker T. Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 2003; 13(11):2498–504. https://doi.org/10.1101/gr.1239303.
Černý M. CyDataSeries - Add time series and the like to your networks. 2017. http://apps.cytoscape.org/apps/cydataseries. Accessed 8 Apr 2019.
Cleveland WS, McGrill R. Graphical Perception: Theory, Experimentation and Application to the Development of Graphical Methods. J Am Stat Assoc. 1984; 79(387):531–54.
Gauges R, Rost U, Sahle S, Wengler K, Bergmann FT. The Systems Biology Markup Language (SBML) Level 3 Package: Layout, Version 1 Core. J Integr Bioinforma. 2015. https://doi.org/10.2390/biecoll-jib-2015-267.
Hucka M, Bergmann FT, Dräger A, Hoops S, Keating SM, Le Novère N, Myers CJ, Olivier BG, Sahle S, Schaff JC, Smith LP, Waltemath D, Wilkinson DJ. Systems Biology Markup Language (SBML) Level 3 Version 1 Core. J Integr Bioinforma. 2018; 15(1):1. https://doi.org/10.1515/jib-2017-0080.
Rougny A, Touré V, Moodie S, Balaur I, Czauderna T, Borlinghaus H, Dogrusoz U, Mazein A, Dräger A, Blinov ML, Villéger AC, Haw R, Demir E, Mi H, Sorokin A, Schreiber F, Luna A. Systems Biology Graphical Notation: Process Description language Level 1 Version 2.0. J Integr Bioinforma. 2019; 16(2). https://doi.org/10.1515/jib-2019-0022.
Bucher J, Riedmaier S, Schnabel A, Marcus K, Vacun G, Weiss TS, Thasler WE, Nüssler AK, Zanger UM, Reuss M. A systems biology approach to dynamic modeling and inter-subject variability of statin pharmacokinetics in human hepatocytes. BMC Syst Biol. 2011; 5(1):66. https://doi.org/10.1186/1752-0509-5-66.
Thomas A, Rahmanian S, Bordbar A, Palsson BO, Jamshidi N. Network reconstruction of platelet metabolism identifies metabolic signature for aspirin resistance. Sci Rep. 2014; 4(3925):1–10. https://doi.org/10.1038/srep03925.
Bordbar A, Jamshidi N, Palsson BO. iAB-RBC-283: A proteomically derived knowledge-base of erythrocyte metabolism that can be used to simulate its physiological and patho-physiological states. BMC Syst Biol. 2011; 5(1):110. https://doi.org/10.1186/1752-0509-5-110.
Mo ML, Palsson B, Herrgård MJ. Connecting extracellular metabolomic measurements to intracellular flux states in yeast. BMC Syst Biol. 2009. https://doi.org/10.1186/1752-0509-3-37.
Bergdahl B, Heer D, Sauer U, Hahn-Hägerdal, B author=Van Niel EW. Dynamic metabolomics differentiates between carbon and energy starvation in recombinant saccharomyces cerevisiae fermenting xylose. Biotechnol Biofuels. 2012. https://doi.org/10.1186/1754-6834-5-34.
Chelliah V, Juty N, Ajmera I, Ali R, Dumousseau M, Glont M, Hucka M, Jalowicki G, Keating S, Knight-Schrijver V, Lloret-Villas A, Nath Natarajan K, Pettit J-B, Rodriguez N, Schubert M, Wimalaratne SM, Zhao Y, Hermjakob H, Le Novère N, Laibe C. BioModels: ten-year anniversary. Nucl Acids Res. 2015. https://doi.org/10.1093/nar/gku1181.
King ZA, Lu JS, Dräger A, Miller PC, Federowicz S, Lerman JA, Ebrahim A, Palsson BO, Lewis NE. BiGG Models: A platform for integrating, standardizing and sharing genome-scale models. Nucleic Acids Res. 2016; 44(D1):515–22. https://doi.org/10.1093/nar/gkv1049.
Yurkovich JT, Bordbar A, Sigurjónsson ÓE, Palsson BO. Systems biology as an emerging paradigm in transfusion medicine. BMC Syst Biol. 2018; 12(1). https://doi.org/10.1186/s12918-018-0558-x.
Le Novère N, Hucka M, Mi H, Moodie S, Schreiber F, Sorokin A, Demir E, Wegner K, Aladjem MI, Wimalaratne SM, Bergmann FT, Gauges R, Ghazal P, Kawaji H, Li L, Matsuoka Y, Villéger A, Boyd SE, Calzone L, Courtot M, Dogrusoz U, Freeman TC, Funahashi A, Ghosh S, Jouraku A, Kim S, Kolpakov F, Luna A, Sahle S, Schmidt E, Watterson S, Wu G, Goryanin I, Kell DB, Sander C, Sauro H, Snoep JL, Kohn K, Kitano H. The Systems Biology Graphical Notation. Nat Biotechnol. 2009; 27(8):735–41. https://doi.org/10.1038/nbt.1558.
Yurkovich JT, Yurkovich BJ, Dräger A, Palsson BO, King ZA. A padawan programmer’s guide to developing software libraries. Cell Syst. 2017; 5(5):431–7. https://doi.org/10.1016/j.cels.2017.08.003.
Whillier S, Raftos JE, Sparrow RL, Kuchel PW. The effects of long-term storage of human red blood cells on the glutathione synthesis rate and steady-state concentration. Transfusion. 2011; 51(7):1450–9. https://doi.org/10.1111/j.1537-2995.2010.03026.x..
D’Alessandro A, D’Amici GM, Vaglio S, Zolla L. Time-course investigation of SAGM-stored leukocyte-filtered red bood cell concentrates: from metabolism to proteomics. Haematologica. 2012; 97(1):107–15.
Nemkov T, Sun K, Reisz JA, Song A, Yoshida T, Dunham A, Wither MJ, Francis RO, Roach RC, Dzieciatkowska M, Rogers SC, Doctor A, Kriebardis A, Antonelou M, Papassideri I, Young C, Thomas T, Hansen KC, Spitalnik SL, Xia Y, Zimring JC, Hod EA, D’Alessandro A. Hypoxia modulates the purine salvage pathway and decreases red blood cell and supernatant levels of hypoxanthine during refrigerated storage. Haematologica. 2017. https://doi.org/10.3324/haematol.2017.178608.
Dörr A, Keller R, Zell A, Dräger A. SBMLsimulator: A Java Tool for Model Simulation and Parameter Estimation in Systems Biology. Computation. 2014; 2(4):246–57. https://doi.org/10.3390/computation2040246.
The authors would like to thank Prof. Dr. Robert Feil, Katrin Keppler for her music contribution, Jan D. Rudolph, and Jakob Matthes for implementing a basic Java™ interpreter of SBML layout.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information about details of the method and implementation can be found in file 12859_2020_3415_MOESM1_ESM.pdf.
Additional file 2 The animated movie iAT-PLT-636.mp4 (547.4 MB) about iAT-PLT-636 available at https://youtu.be/GQuT7R-ldS4.
The data set iAT-PLT-636_Data.csv used in the iAT-PLT-636 animation (15 kB).
The pathway map in Escher format iAT-PLT-636_Map.json.gz compressed with gzip (200 kB).
The pathway map in SBML Level 3 Version 1 format with Layout and FBC (flux balance constraints) packages iAT-PLT-636_Map.xml.gz compressed with gzip (457 kB).
Additional file 6 The animated movie iAB-RBC-283.mp4 (118.9 MB) about iAB-RBC-283 available at https://youtu.be/0INItST4FQc.
The data set iAB-RBC-283_Data.zip used in the iAB-RBC-283 animation (13 kB).
The pathway map in Escher format iAB-RBC-283_Map_cell-outline.json.gz (94 kB).
The pathway map iAB-RBC-283_Map_cell-outline.xml.gz in SBML Level 3 Version 1 format with Layout and FBC packages compressed with gzip (153 kB). Please note that the core model in this file has been reduced to the content of the pathway map and does not comprise all reactions and metabolites of the original model by .
A brief Java™ program that embeds a metabolic map in the format of SBML Layout extensions into an SBML Level 3 Version 1 file.
This document describes the features of the application SBMLsimulator and explains how to use them.
About this article
Cite this article
Buchweitz, L.F., Yurkovich, J.T., Blessing, C. et al. Visualizing metabolic network dynamics through time-series metabolomic data. BMC Bioinformatics 21, 130 (2020). https://doi.org/10.1186/s12859-020-3415-z