CumoVis is an Omix plug-in for three dimensional visualization of the entire atom transport layer of a metabolic network. The atom transport layer is a graph that explains the fate of all individual atoms of all metabolite molecules involved in a metabolic network. CumoVis uses the third dimension in order to show all atom mappings. Furthermore, CumoVis can visualize the so called isotopomer graph and cumomer graph arising from the atom transport layer. In CumoVis you can perform network decomposition and atom path tracing.
Table of contents
- 1.2. Using CumoVis
- 1.3. Further Reading
The individual molecules consist of a atoms which are transported from the educt molecules to the product molecules of a reaction. The description of this transport is called an atom mapping or atom transition. In an atom mapping, series of atoms can be split, joined as well as recombined. The totality of the atom mappings of all reactions lead to the atomic layer:
Isotope labeling is an inportant tool in metabolic flux analysis. Here, the input substrate given into the biological system is isotopically labeled at specific positions in the molecule. The labeled substrate is metabolized by the system and thus, the isotopic labeling is distributed all over the metabolic network. The isotopic enrichment can be measured and allows inference on metabolic flux rates within the network.
A molecule with n labeling positions (being the atoms of the chemical element used for isotope labeling, e.g. carbon) has 2n different possible labeling states which we call isotopomers. The isotopomer network describes how these individual isotopomers are transported, split, joined and recombined in the biological system:
Shortly spoken, cumomers are cumulative sets of isotopomers. Instead of tracing labeled atoms, the cumomer network describes the fate of labeled molecule fragments. By this, splitting reactions does no longer have to be considered and the entire network can be substructured in a cascade of different levels. Inside of such a level, all cumomers have the same weight, i.e. no splitting or joining steps have to be considered.
A CC can be further decomposed into different stroncly connected components (SCC) by considering the actual edge direction. An SCC is a set of nodes within a CC whereas each nodes can reach all other nodes via a certain path. This, SCCs represent the cyclic subnetworks in a CC.
Atomic Layer Plug-In. Currently, CumoVis can only visualize the mappings of carbon atoms and arising isotopomer/cumomer networks. It is recommended to use atom mappings that base upon chemical structures.
Now, CumoVis opens as an additional window and loads the metabolic network model given in Omix.
The network appears different from the classic design known from Omix. The metabolites are displayed with their respective chemical structure if available. Otherwise, a tube-like item us shown filled with as many spheres as labeling positions the metabolite has. The metabolite's position corresponds to its position in the Omix diagram. A text label on top of the chemical structure item shows the metabolites name.
Biochemical reactions in CumoVis are displayed as red spheres placed at the same position as in the Omix diagram. However, the connections in CumoVis do not use the edge shape from Omix. Instead, a standard curve progression is computed from the positions of the connected nodes which is similar for all types and levels of graphs visualized by CumoVis.
On top, a toolbar provides access to different functionality of the tool. In the center of the window the network is shown in 3D. The left side of CumoVis is the component view for navigating between atomic layer, isotopomer network and cumomer cascade.
- Moving the mouse with pressed left mouse button rotates the camera around the scene center point.
- Moving the mouse with pressed right mouse button moves the camera relative to the scene center (up, down, left, right).
- Using the mouse wheel moves the camera towards or away from the scene center point.
- Arrow keys ↑ and ↓ move the network fore- and backwards.
- Arrow keys ← and → rotate the network scene around its yaw axis.
- Page up and down keys rotate the network scene around its pitch axis.
- + and - keys can be used to zoom in and out.
- Key "u" moves the camera up and "d" moves it down.
- "Cumomer Cascade"
- "Isotopomer Graph"
- "Atom Transition Graph"
- "Metabolic Network"
CumoVis visualizes the atom mappings by showing individual items for the single atoms. Their connections show how the atom is transported in the biochemical reaction:
In order to view the isotopomer network and cumomer cascade click on the corresponding entries in the component view. The internal model of the cumomer and isotopomer network is only generated on demand because depending on the network size it can be very memory and runtime expensive. Thus, CumoVis will ask you to confirm loading the cumomer / isotopomer network.
The nodes in the atomix layer, isotopomer network as well as cumomer levels do not display the metabolite name by default. However, you can cause the metabolite name of individual items to be displayed by moving the mouse curser over the corresponding item in the scene.
The "Cumomer Cascade" item contains one child element for each weight level of the cascade. By clicking on a weight level item, only the selected weight level is displayed in the scene.
For the isotopomer network, no network decomposition is provided.
Here, you have to specify the start metabolite and select the atoms to be traced (start isotopomer). You have to select at lease one atom position up to all available positions. After confirming the dialog, a second dialog appears where you have to configure the path tracing:
Here, you have following options:
- Target - If you speficy a target isotopomer CumoVis searchs for all possible paths between the start and end node through the network.
- Maximal Depth - Since path tracing is a extreme complex computational task highly depending on the network size, the maximal depth restricts the search to the environment of the start isotopomer. Take a small maximal depth and increase it iteratively.
- Tracing - You have the coice betreen path tracing on the cumomer weight level or in the isotopomer network. If you select weight level, the path tracing does not include paths where weight is lost or fragments are recombined.
- "Atom specific" - means that each individual atom of the start node is color coded and the successor nodes display the colors of the source atom.
- "Direction" - you can decide to perform path tracing in forward as well as in backward direction. Furthermore, you can decide if reversible reactions are traced in both directions or if only their nominal direction is taken into account.
The path tracing result usually consists of a certain number of paths through the network. On the left side of the toolbar you find a spinbox where you can select the currently displayed path. By pressing the "Shortest" button you can restrict the available paths to the shortest ones (usually multiple). By pressing the "Longest" button the longest paths are shown. By pressing "Show all" all paths are displayed simultaneously.
Furhermore, the displayed path depth can be specified. This is initially 1 which causes only the start node to be displayed. The slider in the center of the toolbar can be used in order to manually specify the path depth. Furthermore, you can customize the number of simultaneously displayed path subsections by changing the length in the spinbox (center).
The media-player-like controls allow to visualize the path tracing as animation.
The path tracing results can be saved as CSV file in list or matrix form.
By pressing "Continue" you can perform another path tracing step at the end of the paths of the current path tracinf result.
In the subsequently appearing configuration window you can specify following program settings:
- Coordinate scale - specify a scale factor which is applied to the coordinates taken from the Omix diagram layout. This option is useful if the distance of the metabolite items is too high or low.
- Edge shape resolution - the shape of the connecting edges is the most expensive part of the 3D scene. In order to improve memory and runtime performance reduce the resolution. In order to get better visual results increase it.
- Chemical structures - you can force CumoVis not to display chemical structures but to use the tube-like shapes instead.
- Item fonts - you can customize the font of metabolites and reactions.
- Wiechert W, Möllney M, Isermann N, Wurzel M, de Graaf AA (1999) Bidirectional reaction steps in metabolic networks. Part III: Explicit solution and analysis of isotopomer labeling systems. Biotechnol Bioeng 66:69–85
- Wiechert W, Wurzel M (2001) Metabolic isotopomer labeling systems. Part I: Global dynamic behaviour. Math Biosci 169:173–205
- Weitzel, M., Wiechert, W., and Nöh, K. 2007. The topology of metabolic isotope labeling networks. BMC Bioinformatics, 8 (1), 315.
- Nöh, K., Weitzel, M., and Wiechert, W. 2008. From Isotope Labeling Patterns to Metabolic Flux Rates. In U. H. E. Hansmann, J. H. Meinke, S. Mohanty, W. Nadler, and O. Zimmermann (Eds.) From Computational Biophysics to Systems Biology (CBSB08), vol. 40 of NIC Series, (pp. 345–348).
- Droste, P., Weitzel, M., and Wiechert, W. 2008b. Visual Exloration of Isotope Labeling Networks in 3D. Bioprocess and Biosystems Engineering, 31 (3), 227–239.