FluxMap: a VANTED add-on for the visual exploration of flux distributions in biological networks

A next-generation dynamic programming language Julia: Its features and applications in biological science

BACKGROUND

The advent of Julia as a sophisticated and dynamic programming language in 2012 represented a significant milestone in computational programming, mathematical analysis, and statistical modeling. Having reached its stable release in version 1.9.0 on May 7, 2023, Julia has developed into a powerful and versatile instrument. Despite its potential and widespread adoption across various scientific and technical domains, there exists a noticeable knowledge gap in comprehending its utilization within biological sciences.

THE AIM OF REVIEW

This comprehensive review aims to address this particular knowledge gap and offer a thorough examination of Julia's fundamental characteristics and its applications in biology.

KEY SCIENTIFIC CONCEPTS OF THE REVIEW

The review focuses on a research gap in the biological science. The review aims to equip researchers with knowledge and tools to utilize Julia's capabilities in biological science effectively and to demonstrate the gap. It paves the way for innovative solutions and discoveries in this rapidly evolving field. It encompasses an analysis of Julia's characteristics, packages, and performance compared to the other programming languages in this field. The initial part of this review discusses the key features of Julia, such as its dynamic and interactive nature, fast processing speed, ease of expression manipulation, user-friendly syntax, code readability, strong support for multiple dispatch, and advanced type system. It also explores Julia's capabilities in data analysis, visualization, machine learning, and algorithms, making it suitable for scientific applications. The next section emphasizes the importance of using Julia in biological research, highlighting its seamless integration with biological studies for data analysis, and computational biology. It also compares Julia with other programming languages commonly used in biological research through benchmarking and performance analysis. Additionally, it provides insights into future directions and potential challenges in Julia's applications in biology.

Metabolic Rearrangements Causing Elevated Proline and Polyhydroxybutyrate Accumulation During the Osmotic Adaptation Response of Bacillus megaterium

For many years now, Bacillus megaterium serves as a microbial workhorse for the high-level production of recombinant proteins in the g/L-scale. However, efficient and stable production processes require the knowledge of the molecular adaptation strategies of the host organism to establish optimal environmental conditions. Here, we interrogated the osmotic stress response of B. megaterium using transcriptome, proteome, metabolome, and fluxome analyses. An initial transient adaptation consisted of potassium import and glutamate counterion synthesis. The massive synthesis of the compatible solute proline constituted the second longterm adaptation process. Several stress response enzymes involved in iron scavenging and reactive oxygen species (ROS) fighting proteins showed higher levels under prolonged osmotic stress induced by 1.8 M NaCl. At the same time, the downregulation of the expression of genes of the upper part of glycolysis resulted in the activation of the pentose phosphate pathway (PPP), generating an oversupply of NADPH. The increased production of lactate accompanied by the reduction of acetate secretion partially compensate for the unbalanced (NADH/NAD⁺) ratio. Besides, the tricarboxylic acid cycle (TCA) mainly supplies the produced NADH, as indicated by the higher mRNA and protein levels of involved enzymes, and further confirmed by ¹³C flux analyses. As a consequence of the metabolic flux toward acetyl-CoA and the generation of an excess of NADPH, B. megaterium redirected the produced acetyl-CoA toward the polyhydroxybutyrate (PHB) biosynthetic pathway accumulating around 30% of the cell dry weight (CDW) as PHB. This direct relation between osmotic stress and intracellular PHB content has been evidenced for the first time, thus opening new avenues for synthesizing this valuable biopolymer using varying salt concentrations under non-limiting nutrient conditions.

A Genome-Scale Metabolic Model of Soybean (Glycine max) Highlights Metabolic Fluxes in Seedlings

Until they become photoautotrophic juvenile plants, seedlings depend upon the reserves stored in seed tissues. These reserves must be mobilized and metabolized, and their breakdown products must be distributed to the different organs of the growing seedling. Here, we investigated the mobilization of soybean (Glycine max) seed reserves during seedling growth by initially constructing a genome-scale stoichiometric model for this important crop plant and then adapting the model to reflect metabolism in the cotyledons and hypocotyl/root axis (HRA). A detailed analysis of seedling growth and alterations in biomass composition was performed over 4 d of postgerminative growth and used to constrain the stoichiometric model. Flux balance analysis revealed marked differences in metabolism between the two organs, together with shifts in primary metabolism occurring during different periods postgermination. In particular, from 48 h onward, cotyledons were characterized by the oxidation of fatty acids to supply carbon for the tricarboxylic acid cycle as well as production of sucrose and glutamate for export to the HRA, while the HRA was characterized by the use of a range of imported amino acids in protein synthesis and catabolic processes. Overall, the use of flux balance modeling provided new insight into well-characterized metabolic processes in an important crop plant due to their analysis within the context of a metabolic network and reinforces the relevance of the application of this technique to the analysis of complex plant metabolic systems.

Mental-Map Preserving Visualisation of Partitioned Networks in Vanted

Biological networks can be large and complex, often consisting of different sub-networks or parts. Separation of networks into parts, network partitioning and layouts of overview and sub-graphs are of importance for understandable visualisations of those networks. This article presents NetPartVis to visualise non-overlapping clusters or partitions of graphs in the Vanted framework based on a method for laying out overview graph and several sub-graphs (partitions) in a coordinated, mental-map preserving way.

Novel column generation-based optimization approach for poly-pathway kinetic model applied to CHO cell culture

Mathematical modelling can provide precious tools for bioprocess simulation, prediction, control and optimization of mammalian cell-based cultures. In this paper we present a novel method to generate kinetic models of such cultures, rendering complex metabolic networks in a poly-pathway kinetic model. The model is based on subsets of elementary flux modes (EFMs) to generate macro-reactions. Thanks to our column generation-based optimization algorithm, the experimental data are used to identify the EFMs, which are relevant to the data. Here the systematic enumeration of all the EFMs is eliminated and a network including a large number of reactions can be considered. In particular, the poly-pathway model can simulate multiple metabolic behaviors in response to changes in the culture conditions.

We apply the method to a network of 126 metabolic reactions describing cultures of antibody-producing Chinese hamster ovary cells, and generate a poly-pathway model that simulates multiple experimental conditions obtained in response to variations in amino acid availability. A good fit between simulated and experimental data is obtained, rendering the variations in the growth, product, and metabolite uptake/secretion rates. The intracellular reaction fluxes simulated by the model are explored, linking variations in metabolic behavior to adaptations of the intracellular metabolism.

•

Novel method to model multiple states by a poly-pathway kinetic model.

•

EFMs relevant to data identified by column generation (CG)-based optimization.

•

CG optimization enables use of networks much larger than systematic enumeration.

•

A kinetic model simulates changes in metabolic rates linked to available amino acids.

•

The flux distribution of each metabolic state is visualized in the original network.

Recent developments in synthetic biology and metabolic engineering in microalgae towards biofuel production

In the wake of the uprising global energy crisis, microalgae have emerged as an alternate feedstock for biofuel production. In addition, microalgae bear immense potential as bio-cell factories in terms of producing key chemicals, recombinant proteins, enzymes, lipid, hydrogen and alcohol. Abstraction of such high-value products (algal biorefinery approach) facilitates to make microalgae-based renewable energy an economically viable option. Synthetic biology is an emerging field that harmoniously blends science and engineering to help design and construct novel biological systems, with an aim to achieve rationally formulated objectives. However, resources and tools used for such nuclear manipulation, construction of synthetic gene network and genome-scale reconstruction of microalgae are limited. Herein, we present recent developments in the upcoming field of microalgae employed as a model system for synthetic biology applications and highlight the importance of genome-scale reconstruction models and kinetic models, to maximize the metabolic output by understanding the intricacies of algal growth. This review also examines the role played by microalgae as biorefineries, microalgal culture conditions and various operating parameters that need to be optimized to yield biofuel that can be economically competitive with fossil fuels.

Fluxomics links cellular functional analyses to whole-plant phenotyping

Fluxes through metabolic pathways reflect the integration of genetic and metabolic regulations. While it is attractive to measure all the mRNAs (transcriptome), all the proteins (proteome), and a large number of the metabolites (metabolome) in a given cellular system, linking and integrating this information remains difficult. Measurement of metabolome-wide fluxes (termed the fluxome) provides an integrated functional output of the cell machinery and a better tool to link functional analyses to plant phenotyping. This review presents and discusses sets of methodologies that have been developed to measure the fluxome. First, the principles of metabolic flux analysis (MFA), its 'short time interval' version Inst-MFA, and of constraints-based methods, such as flux balance analysis and kinetic analysis, are briefly described. The use of these powerful methods for flux characterization at the cellular scale up to the organ (fruits, seeds) and whole-plant level is illustrated. The added value given by fluxomics methods for unravelling how the abiotic environment affects flux, the process, and key metabolic steps are also described. Challenges associated with the development of fluxomics and its integration with 'omics' for thorough plant and organ functional phenotyping are discussed. Taken together, these will ultimately provide crucial clues for identifying appropriate target plant phenotypes for breeding.

Isotopically Nonstationary Metabolic Flux Analysis (INST-MFA) of Photosynthesis and Photorespiration in Plants

Photorespiration is a central component of photosynthesis; however to better understand its role it should be viewed in the context of an integrated metabolic network rather than a series of individual reactions that operate independently. Isotopically nonstationary ¹³C metabolic flux analysis (INST-MFA), which is based on transient labeling studies at metabolic steady state, offers a comprehensive platform to quantify plant central metabolism. In this chapter, we describe the application of INST-MFA to investigate metabolism in leaves. Leaves are an autotrophic tissue, assimilating CO₂ over a diurnal period implying that the metabolic steady state is limited to less than 12 h and thus requiring an INST-MFA approach. This strategy results in a comprehensive unified description of photorespiration, Calvin cycle, sucrose and starch synthesis, tricarboxylic acid (TCA) cycle, and amino acid biosynthetic fluxes. We present protocols of the experimental aspects for labeling studies: transient ¹³CO₂ labeling of leaf tissue, sample quenching and extraction, mass spectrometry (MS) analysis of isotopic labeling data, measurement of sucrose and amino acids in vascular exudates, and provide details on the computational flux estimation using INST-MFA.

ReactionFlow: an interactive visualization tool for causality analysis in biological pathways

Background

Molecular and systems biologists are tasked with the comprehension and analysis of incredibly complex networks of biochemical interactions, called pathways, that occur within a cell. Through interviews with domain experts, we identified four common tasks that require an understanding of the causality within pathways, that is, the downstream and upstream relationships between proteins and biochemical reactions, including: visualizing downstream consequences of perturbing a protein; finding the shortest path between two proteins; detecting feedback loops within the pathway; and identifying common downstream elements from two or more proteins.

Results

We introduce ReactionFlow, a visual analytics application for pathway analysis that emphasizes the structural and causal relationships amongst proteins, complexes, and biochemical reactions within a given pathway. To support the identified causality analysis tasks, user interactions allow an analyst to filter, cluster, and select pathway components across linked views. Animation is used to highlight the flow of activity through a pathway.

Conclusions

We evaluated ReactionFlow by providing our application to two domain experts who have significant experience with biomolecular pathways, after which we conducted a series of in-depth interviews focused on each of the four causality analysis tasks. Their feedback leads us to believe that our techniques could be useful to researchers who must be able to understand and analyze the complex nature of biological pathways. ReactionFlow is available at https://github.com/CreativeCodingLab/ReactionFlow.

Flux visualization using VANTED/FluxMap

The calculation of metabolic fluxes has been shown to be a valuable asset in systems biology. Several procedures are commonly used to achieve this. Flux balance analyses or metabolic flux analyses usually result in a list of reaction rates (fluxes) provided in a spreadsheet format. This makes it difficult to quickly assess general characteristics of the solution. A fast and easy mapping of these results to a graphical map template facilitates an easy visual data inspection. Here, we describe a protocol that helps in setting up user-specific network templates, mapping flux results to it, and creating multiple exportable flux maps at one time.

CeCaFDB: a curated database for the documentation, visualization and comparative analysis of central carbon metabolic flux distributions explored by 13C-fluxomics

The Central Carbon Metabolic Flux Database (CeCaFDB, available at http://www.cecafdb.org) is a manually curated, multipurpose and open-access database for the documentation, visualization and comparative analysis of the quantitative flux results of central carbon metabolism among microbes and animal cells. It encompasses records for more than 500 flux distributions among 36 organisms and includes information regarding the genotype, culture medium, growth conditions and other specific information gathered from hundreds of journal articles. In addition to its comprehensive literature-derived data, the CeCaFDB supports a common text search function among the data and interactive visualization of the curated flux distributions with compartmentation information based on the Cytoscape Web API, which facilitates data interpretation. The CeCaFDB offers four modules to calculate a similarity score or to perform an alignment between the flux distributions. One of the modules was built using an inter programming algorithm for flux distribution alignment that was specifically designed for this study. Based on these modules, the CeCaFDB also supports an extensive flux distribution comparison function among the curated data. The CeCaFDB is strenuously designed to address the broad demands of biochemists, metabolic engineers, systems biologists and members of the -omics community.

Plant metabolic modeling: achieving new insight into metabolism and metabolic engineering

Models are used to represent aspects of the real world for specific purposes, and mathematical models have opened up new approaches in studying the behavior and complexity of biological systems. However, modeling is often time-consuming and requires significant computational resources for data development, data analysis, and simulation. Computational modeling has been successfully applied as an aid for metabolic engineering in microorganisms. But such model-based approaches have only recently been extended to plant metabolic engineering, mainly due to greater pathway complexity in plants and their highly compartmentalized cellular structure. Recent progress in plant systems biology and bioinformatics has begun to disentangle this complexity and facilitate the creation of efficient plant metabolic models. This review highlights several aspects of plant metabolic modeling in the context of understanding, predicting and modifying complex plant metabolism. We discuss opportunities for engineering photosynthetic carbon metabolism, sucrose synthesis, and the tricarboxylic acid cycle in leaves and oil synthesis in seeds and the application of metabolic modeling to the study of plant acclimation to the environment. The aim of the review is to offer a current perspective for plant biologists without requiring specialized knowledge of bioinformatics or systems biology.

Computational approaches for microalgal biofuel optimization: a review

The increased demand and consumption of fossil fuels have raised interest in finding renewable energy sources throughout the globe. Much focus has been placed on optimizing microorganisms and primarily microalgae, to efficiently produce compounds that can substitute for fossil fuels. However, the path to achieving economic feasibility is likely to require strain optimization through using available tools and technologies in the fields of systems and synthetic biology. Such approaches invoke a deep understanding of the metabolic networks of the organisms and their genomic and proteomic profiles. The advent of next generation sequencing and other high throughput methods has led to a major increase in availability of biological data. Integration of such disparate data can help define the emergent metabolic system properties, which is of crucial importance in addressing biofuel production optimization. Herein, we review major computational tools and approaches developed and used in order to potentially identify target genes, pathways, and reactions of particular interest to biofuel production in algae. As the use of these tools and approaches has not been fully implemented in algal biofuel research, the aim of this review is to highlight the potential utility of these resources toward their future implementation in algal research.

Flux analysis in plant metabolic networks: increasing throughput and coverage

Quantitative information about metabolic networks has been mainly obtained at the level of metabolite contents, transcript abundance, and enzyme activities. However, the active process of metabolism is represented by the flow of matter through the pathways. These metabolic fluxes can be predicted by Flux Balance Analysis or determined experimentally by (13)C-Metabolic Flux Analysis. These relatively complicated and time-consuming methods have recently seen significant improvements at the level of coverage and throughput. Metabolic models have developed from single cell models into whole-organism dynamic models. Advances in lab automation and data handling have significantly increased the throughput of flux measurements. This review summarizes advances to increase coverage and throughput of metabolic flux analysis in plants.

High-throughput data pipelines for metabolic flux analysis in plants

In this chapter we illustrate the methodology for high-throughput metabolic flux analysis. Central to this is developing an end to end data pipeline, crucial for integrating the wet lab experiments and analytics, combining hardware and software automation, and standardizing data representation providing importers and exporters to support third party tools. The use of existing software at the start, data extraction from the chromatogram, and the end, MFA analysis, allows for the most flexibility in this workflow. Developing iMS2Flux provided a standard, extensible, platform independent tool to act as the "glue" between these end points. Most importantly this tool can be easily adapted to support different data formats, data verification and data correction steps allowing it to be central to managing the data necessary for high-throughput MFA. An additional tool was needed to automate the MFA software and in particular to take advantage of the course grained parallel nature of high-throughput analysis and available high performance computing facilities.In combination these methods show the development of high-throughput pipelines that allow metabolic flux analysis to join as a full member of the omics family.

Isotopically nonstationary MFA (INST-MFA) of autotrophic metabolism

Metabolic flux analysis (MFA) is a powerful approach for quantifying plant central carbon metabolism based upon a combination of extracellular flux measurements and intracellular isotope labeling measurements. In this chapter, we present the method of isotopically nonstationary (13)C MFA (INST-MFA), which is applicable to autotrophic systems that are at metabolic steady state but are sampled during the transient period prior to achieving isotopic steady state following the introduction of (13)CO2. We describe protocols for performing the necessary isotope labeling experiments, sample collection and quenching, nonaqueous fractionation and extraction of intracellular metabolites, and mass spectrometry (MS) analysis of metabolite labeling. We also outline the steps required to perform computational flux estimation using INST-MFA. By combining several recently developed experimental and computational techniques, INST-MFA provides an important new platform for mapping carbon fluxes that is especially applicable to autotrophic organisms, which are not amenable to steady-state (13)C MFA experiments.

Integrated network analysis and effective tools in plant systems biology

One of the ultimate goals in plant systems biology is to elucidate the genotype-phenotype relationship in plant cellular systems. Integrated network analysis that combines omics data with mathematical models has received particular attention. Here we focus on the latest cutting-edge computational advances that facilitate their combination. We highlight (1) network visualization tools, (2) pathway analyses, (3) genome-scale metabolic reconstruction, and (4) the integration of high-throughput experimental data and mathematical models. Multi-omics data that contain the genome, transcriptome, proteome, and metabolome and mathematical models are expected to integrate and expand our knowledge of complex plant metabolisms.

3DScapeCS: application of three dimensional, parallel, dynamic network visualization in Cytoscape

Background

The exponential growth of gigantic biological data from various sources, such as protein-protein interaction (PPI), genome sequences scaffolding, Mass spectrometry (MS) molecular networking and metabolic flux, demands an efficient way for better visualization and interpretation beyond the conventional, two-dimensional visualization tools.

Results

We developed a 3D Cytoscape Client/Server (3DScapeCS) plugin, which adopted Cytoscape in interpreting different types of data, and UbiGraph for three-dimensional visualization. The extra dimension is useful in accommodating, visualizing, and distinguishing large-scale networks with multiple crossed connections in five case studies.

Conclusions

Evaluation on several experimental data using 3DScapeCS and its special features, including multilevel graph layout, time-course data animation, and parallel visualization has proven its usefulness in visualizing complex data and help to make insightful conclusions.

Structure and dynamics of molecular networks: a novel paradigm of drug discovery: a comprehensive review

Despite considerable progress in genome- and proteome-based high-throughput screening methods and in rational drug design, the increase in approved drugs in the past decade did not match the increase of drug development costs. Network description and analysis not only give a systems-level understanding of drug action and disease complexity, but can also help to improve the efficiency of drug design. We give a comprehensive assessment of the analytical tools of network topology and dynamics. The state-of-the-art use of chemical similarity, protein structure, protein-protein interaction, signaling, genetic interaction and metabolic networks in the discovery of drug targets is summarized. We propose that network targeting follows two basic strategies. The "central hit strategy" selectively targets central nodes/edges of the flexible networks of infectious agents or cancer cells to kill them. The "network influence strategy" works against other diseases, where an efficient reconfiguration of rigid networks needs to be achieved by targeting the neighbors of central nodes/edges. It is shown how network techniques can help in the identification of single-target, edgetic, multi-target and allo-network drug target candidates. We review the recent boom in network methods helping hit identification, lead selection optimizing drug efficacy, as well as minimizing side-effects and drug toxicity. Successful network-based drug development strategies are shown through the examples of infections, cancer, metabolic diseases, neurodegenerative diseases and aging. Summarizing >1200 references we suggest an optimized protocol of network-aided drug development, and provide a list of systems-level hallmarks of drug quality. Finally, we highlight network-related drug development trends helping to achieve these hallmarks by a cohesive, global approach.

Isotopically nonstationary 13C metabolic flux analysis

(13)C metabolic flux analysis (MFA) is a powerful approach for quantifying cell physiology based upon a combination of extracellular flux measurements and intracellular isotope labeling measurements. In this chapter, we present the method of isotopically nonstationary (13)C MFA (INST-MFA), which is applicable to systems that are at metabolic steady state, but are sampled during the transient period prior to achieving isotopic steady state following the introduction of a (13)C tracer. We describe protocols for performing the necessary isotope labeling experiments, for quenching and extraction of intracellular metabolites, for mass spectrometry (MS) analysis of metabolite labeling, and for computational flux estimation using INST-MFA. By combining several recently developed experimental and computational techniques, INST-MFA provides an important new platform for mapping carbon fluxes that is especially applicable to animal cell cultures, autotrophic organisms, industrial bioprocesses, high-throughput experiments, and other systems that are not amenable to steady-state (13)C MFA experiments.

VANTED v2: a framework for systems biology applications

BACKGROUND

Experimental datasets are becoming larger and increasingly complex, spanning different data domains, thereby expanding the requirements for respective tool support for their analysis. Networks provide a basis for the integration, analysis and visualization of multi-omics experimental datasets.

RESULTS

Here we present VANTED (version 2), a framework for systems biology applications, which comprises a comprehensive set of seven main tasks. These range from network reconstruction, data visualization, integration of various data types, network simulation to data exploration combined with a manifold support of systems biology standards for visualization and data exchange. The offered set of functionalities is instantiated by combining several tasks in order to enable users to view and explore a comprehensive dataset from different perspectives. We describe the system as well as an exemplary workflow.

CONCLUSIONS

VANTED is a stand-alone framework which supports scientists during the data analysis and interpretation phase. It is available as a Java open source tool from http://www.vanted.org.

Software applications toward quantitative metabolic flux analysis and modeling

Metabolites and their pathways are central for adaptation and survival. Metabolic modeling elucidates in silico all the possible flux pathways (flux balance analysis, FBA) and predicts the actual fluxes under a given situation, further refinement of these models is possible by including experimental isotopologue data. In this review, we initially introduce the key theoretical concepts and different analysis steps in the modeling process before comparing flux calculation and metabolite analysis programs such as C13, BioOpt, COBRA toolbox, Metatool, efmtool, FiatFlux, ReMatch, VANTED, iMAT and YANA. Their respective strengths and limitations are discussed and compared to alternative software. While data analysis of metabolites, calculation of metabolic fluxes, pathways and their condition-specific changes are all possible, we highlight the considerations that need to be taken into account before deciding on a specific software. Current challenges in the field include the computation of large-scale networks (in elementary mode analysis), regulatory interactions and detailed kinetics, and these are discussed in the light of powerful new approaches.