Metabolic modeling of microbial strains in silico

Genome-scale metabolic modelling reveals interactions and key roles of symbiont clades in a sponge holobiont

Sponges harbour complex microbiomes and as ancient metazoans and important ecosystem players are emerging as powerful models to understand the evolution and ecology of symbiotic interactions. Metagenomic studies have previously described the functional features of sponge symbionts, however, little is known about the metabolic interactions and processes that occur under different environmental conditions. To address this issue, we construct here constraint-based, genome-scale metabolic networks for the microbiome of the sponge Stylissa sp. Our models define the importance of sponge-derived nutrients for microbiome stability and discover how different organic inputs can result in net heterotrophy or autotrophy of the symbiont community. The analysis further reveals the key role that a newly discovered bacterial taxon has in cross-feeding activities and how it dynamically adjusts with nutrient inputs. Our study reveals insights into the functioning of a sponge microbiome and provides a framework to further explore and define metabolic interactions in holobionts.

Refining microbial community metabolic models derived from metagenomics using reference-based taxonomic profiling

Characterization of microbial community metabolic output is crucial to understanding their functions. Construction of genome-scale metabolic models from metagenome-assembled genomes (MAG) has enabled prediction of metabolite production by microbial communities, yet little is known about their accuracy. Here, we examined the performance of two approaches for metabolite prediction from metagenomes, one that is MAG-guided and another that is taxonomic reference-guided. We applied both on shotgun metagenomics data from human and environmental samples, and validated findings in the human samples using untargeted metabolomics. We found that in human samples, where taxonomic profiling is optimized and reference genomes are readily available, when number of input taxa was normalized, the reference-guided approach predicted more metabolites than the MAG-guided approach. The two approaches showed significant overlap but each identified metabolites not predicted in the other. Pathway enrichment analyses identified significant differences in inferences derived from data based on the approach, highlighting the need for caution in interpretation. In environmental samples, when the number of input taxa was normalized, the reference-guided approach predicted more metabolites than the MAG-guided approach for total metabolites in both sample types and non-redundant metabolites in seawater samples. Nonetheless, as was observed for the human samples, the approaches overlapped substantially but also predicted metabolites not observed in the other. Our findings report on utility of a complementary input to genome-scale metabolic model construction that is less computationally intensive forgoing MAG assembly and refinement, and that can be applied on shallow shotgun sequencing where MAGs cannot be generated.IMPORTANCELittle is known about the accuracy of genome-scale metabolic models (GEMs) of microbial communities despite their influence on inferring community metabolic outputs and culture conditions. The performance of GEMs for metabolite prediction from metagenomes was assessed by applying two approaches on shotgun metagenomics data from human and environmental samples, and validating findings in the human samples using untargeted metabolomics. The performance of the approach was found to be dependent on sample type, but collectively, the reference-guided approach predicted more metabolites than the MAG-guided approach. Despite the differences, the predictions from the approaches overlapped substantially but each identified metabolites not predicted in the other. We found significant differences in biological inferences based on the approach, with some examples of uniquely enriched pathways in one group being invalidated when using the alternative approach, highlighting the need for caution in interpretation of GEMs.

Modeling lignin biosynthesis: a pathway to renewable chemicals

Plant biomass contains lignin that can be converted into high-value-added chemicals, fuels, and materials. The precise genetic manipulation of lignin content and composition in plant cells offers substantial environmental and economic benefits. However, the intricate regulatory mechanisms governing lignin formation challenge the development of crops with specific lignin profiles. Mathematical models and computational simulations have recently been employed to gain fundamental understanding of the metabolism of lignin and related phenolic compounds. This review article discusses the strategies used for modeling plant metabolic networks, focusing on the application of mathematical modeling for flux network analysis in monolignol biosynthesis. Furthermore, we highlight how current challenges might be overcome to optimize the use of metabolic modeling approaches for developing lignin-engineered plants.

Characterization of the enzyme kinetics of EMP and HMP pathway in Corynebacterium glutamicum: reference for modeling metabolic networks

The model of intracellular metabolic network based on enzyme kinetics parameters plays an important role in understanding the intracellular metabolic process of Corynebacterium glutamicum, and constructing such a model requires a large number of enzymological parameters. In this work, the genes encoding the relevant enzymes of the EMP and HMP metabolic pathways from Corynebacterium glutamicum ATCC 13032 were cloned, and engineered strains for protein expression with E.coli BL21 and P.pastoris X33 as hosts were constructed. The twelve enzymes (GLK, GPI, TPI, GAPDH, PGK, PMGA, ENO, ZWF, RPI, RPE, TKT, and TAL) were successfully expressed and purified by Ni2+ chelate affinity chromatography in their active forms. In addition, the kinetic parameters (V max, K m, and K cat) of these enzymes were measured and calculated at the same pH and temperature. The kinetic parameters of enzymes associated with EMP and the HMP pathway were determined systematically and completely for the first time in C.glutamicum. These kinetic parameters enable the prediction of key enzymes and rate-limiting steps within the metabolic pathway, and support the construction of a metabolic network model for important metabolic pathways in C.glutamicum. Such analyses and models aid in understanding the metabolic behavior of the organism and can guide the efficient production of high-value chemicals using C.glutamicum as a host.

From genotype to phenotype: computational approaches for inferring microbial traits relevant to the food industry

When selecting microbial strains for the production of fermented foods, various microbial phenotypes need to be taken into account to achieve target product characteristics, such as biosafety, flavor, texture, and health-promoting effects. Through continuous advances in sequencing technologies, microbial whole-genome sequences of increasing quality can now be obtained both cheaper and faster, which increases the relevance of genome-based characterization of microbial phenotypes. Prediction of microbial phenotypes from genome sequences makes it possible to quickly screen large strain collections in silico to identify candidates with desirable traits. Several microbial phenotypes relevant to the production of fermented foods can be predicted using knowledge-based approaches, leveraging our existing understanding of the genetic and molecular mechanisms underlying those phenotypes. In the absence of this knowledge, data-driven approaches can be applied to estimate genotype-phenotype relationships based on large experimental datasets. Here, we review computational methods that implement knowledge- and data-driven approaches for phenotype prediction, as well as methods that combine elements from both approaches. Furthermore, we provide examples of how these methods have been applied in industrial biotechnology, with special focus on the fermented food industry.

"Minimal metabolism": A key concept to investigate the origins and nature of biological systems

The systems view on life and its emergence from complex chemistry has remarkably increased the scientific attention on metabolism in the last two decades. However, during this time there has not been much theoretical discussion on what constitutes a metabolism and what role it actually played in biogenesis. A critical and updated review on the topic is here offered, including some references to classical models from last century, but focusing more on current and future research. Metabolism is considered as intrinsically related to the living but not necessarily equivalent to it. More precisely, the idea of "minimal metabolism", in contrast to previous, top-down conceptions, is formulated as a heuristic construct, halfway between chemistry and biology. Thus, rather than providing a complete or final characterization of metabolism, our aim is to encourage further investigations on it, particularly in the context of life's origin, for which some concrete methodological suggestions are provided. Also see the video abstract here: https://youtu.be/DP7VMKk2qpA.

Towards a Systems Biology Approach to Understanding the Lichen Symbiosis: Opportunities and Challenges of Implementing Network Modelling

Lichen associations, a classic model for successful and sustainable interactions between micro-organisms, have been studied for many years. However, there are significant gaps in our understanding about how the lichen symbiosis operates at the molecular level. This review addresses opportunities for expanding current knowledge on signalling and metabolic interplays in the lichen symbiosis using the tools and approaches of systems biology, particularly network modelling. The largely unexplored nature of symbiont recognition and metabolic interdependency in lichens could benefit from applying a holistic approach to understand underlying molecular mechanisms and processes. Together with ‘omics’ approaches, the application of signalling and metabolic network modelling could provide predictive means to gain insights into lichen signalling and metabolic pathways. First, we review the major signalling and recognition modalities in the lichen symbioses studied to date, and then describe how modelling signalling networks could enhance our understanding of symbiont recognition, particularly leveraging omics techniques. Next, we highlight the current state of knowledge on lichen metabolism. We also discuss metabolic network modelling as a tool to simulate flux distribution in lichen metabolic pathways and to analyse the co-dependence between symbionts. This is especially important given the growing number of lichen genomes now available and improved computational tools for reconstructing such models. We highlight the benefits and possible bottlenecks for implementing different types of network models as applied to the study of lichens.

Strategies for Enhancing in vitro Degradation of Linuron by Variovorax sp. Strain SRS 16 Under the Guidance of Metabolic Modeling

Phenyl urea herbicides are being extensively used for weed control in both agricultural and non-agricultural applications. Linuron is one of the key herbicides in this family and is in wide use. Like other phenyl urea herbicides, it is known to have toxic effects as a result of its persistence in the environment. The natural removal of linuron from the environment is mainly carried through microbial biodegradation. Some microorganisms have been reported to mineralize linuron completely and utilize it as a carbon and nitrogen source. Variovorax sp. strain SRS 16 is one of the known efficient degraders with a recently sequenced genome. The genomic data provide an opportunity to use a genome-scale model for improving biodegradation. The aim of our study is the construction of a genome-scale metabolic model following automatic and manual protocols and its application for improving its metabolic potential through iterative simulations. Applying flux balance analysis (FBA), growth and degradation performances of SRS 16 in different media considering the influence of selected supplements (potential carbon and nitrogen sources) were simulated. Outcomes are predictions for the suitable media modification, allowing faster degradation of linuron by SRS 16. Seven metabolites were selected for in vitro validation of the predictions through laboratory experiments confirming the degradation-promoting effect of specific amino acids (glutamine and asparagine) on linuron degradation and SRS 16 growth. Overall, simulations are shown to be efficient in predicting the degradation potential of SRS 16 in the presence of specific supplements. The generated information contributes to the understanding of the biochemistry of linuron degradation and can be further utilized for the development of new cleanup solutions without any genetic manipulation.

Modelling Cell Metabolism: A Review on Constraint-Based Steady-State and Kinetic Approaches

Studying cell metabolism serves a plethora of objectives such as the enhancement of bioprocess performance, and advancement in the understanding of cell biology, of drug target discovery, and in metabolic therapy. Remarkable successes in these fields emerged from heuristics approaches, for instance, with the introduction of effective strategies for genetic modifications, drug developments and optimization of bioprocess management. However, heuristics approaches have showed significant shortcomings, such as to describe regulation of metabolic pathways and to extrapolate experimental conditions. In the specific case of bioprocess management, such shortcomings limit their capacity to increase product quality, while maintaining desirable productivity and reproducibility levels. For instance, since heuristics approaches are not capable of prediction of the cellular functions under varying experimental conditions, they may lead to sub-optimal processes. Also, such approaches used for bioprocess control often fail in regulating a process under unexpected variations of external conditions. Therefore, methodologies inspired by the systematic mathematical formulation of cell metabolism have been used to address such drawbacks and achieve robust reproducible results. Mathematical modelling approaches are effective for both the characterization of the cell physiology, and the estimation of metabolic pathways utilization, thus allowing to characterize a cell population metabolic behavior. In this article, we present a review on methodology used and promising mathematical modelling approaches, focusing primarily to investigate metabolic events and regulation. Proceeding from a topological representation of the metabolic networks, we first present the metabolic modelling approaches that investigate cell metabolism at steady state, complying to the constraints imposed by mass conservation law and thermodynamics of reactions reversibility. Constraint-based models (CBMs) are reviewed highlighting the set of assumed optimality functions for reaction pathways. We explore models simulating cell growth dynamics, by expanding flux balance models developed at steady state. Then, discussing a change of metabolic modelling paradigm, we describe dynamic kinetic models that are based on the mathematical representation of the mechanistic description of nonlinear enzyme activities. In such approaches metabolic pathway regulations are considered explicitly as a function of the activity of other components of metabolic networks and possibly far from the metabolic steady state. We have also assessed the significance of metabolic model parameterization in kinetic models, summarizing a standard parameter estimation procedure frequently employed in kinetic metabolic modelling literature. Finally, some optimization practices used for the parameter estimation are reviewed.

An insight to flux-balance analysis for biochemical networks

Systems biology is one of the integrated ways to study biological systems and is more favourable than the earlier used approaches. It includes metabolic pathway analysis, modelling, and regulatory as well as signal transduction for getting insights into cellular behaviour. Among the various techniques of modelling, simulation, analysis of networks and pathways, flux-based analysis (FBA) has been recognised because of its extensibility as well as simplicity. It is widely accepted because it is not like a mechanistic simulation which depends on accurate kinetic data. The study of fluxes through the network is informative and can give insights even in the absence of kinetic data. FBA is one of the widely used tools to study biochemical networks and needs information of reaction stoichiometry, growth requirements, specific measurement parameters of the biological system, in particular the reconstruction of the metabolic network for the genome-scale, many of which have already been built previously. It defines the boundaries of flux distributions which are possible and achievable with a defined set of genes. This review article gives an insight into FBA, from the extension of flux balancing to mathematical representation followed by a discussion about the formulation of flux-balance analysis problems, defining constraints for the stoichiometry of the pathways and the tools that can be used in FBA such as FASIMA, COBRA toolbox, and OptFlux. It also includes broader areas in terms of applications which can be covered by FBA as well as the queries which can be addressed through FBA.

Diagnosing and Predicting Mixed-Culture Fermentations with Unicellular and Guild-Based Metabolic Models

Microbiomes are vital to human health, agriculture, and protecting the environment. Predicting behavior of self-assembled or synthetic microbiomes, however, remains a challenge. In this work, we used unicellular and guild-based metabolic models to investigate production of medium-chain fatty acids by a mixed microbial community that is fed multiple organic substrates. Modeling results provided insights into metabolic pathways of three medium-chain fatty acid-producing guilds and identified potential strategies to increase production of medium-chain fatty acids. This work demonstrates the role of metabolic models in augmenting multi-omic studies to gain greater insights into microbiome behavior.

Multispecies microbial communities determine the fate of materials in the environment and can be harnessed to produce beneficial products from renewable resources. In a recent example, fermentations by microbial communities have produced medium-chain fatty acids (MCFAs). Tools to predict, assess, and improve the performance of these communities, however, are limited. To provide such tools, we constructed two metabolic models of MCFA-producing microbial communities based on available genomic, transcriptomic, and metabolomic data. The first model is a unicellular model (iFermCell215), while the second model (iFermGuilds789) separates fermentation activities into functional guilds. Ethanol and lactate are fermentation products known to serve as substrates for MCFA production, while acetate is another common cometabolite during MCFA production. Simulations with iFermCell215 predict that low molar ratios of acetate to ethanol favor MCFA production, whereas the products of lactate and acetate coutilization are less dependent on the acetate-to-lactate ratio. In simulations of an MCFA-producing community fed a complex organic mixture derived from lignocellulose, iFermGuilds789 predicted that lactate was an extracellular cometabolite that served as a substrate for butyrate (C4) production. Extracellular hexanoic (C6) and octanoic (C8) acids were predicted by iFermGuilds789 to be from community members that directly metabolize sugars. Modeling results provide several hypotheses that can improve our understanding of microbial roles in an MCFA-producing microbiome and inform strategies to increase MCFA production. Further, these models represent novel tools for exploring the role of mixed microbial communities in carbon recycling in the environment, as well as in beneficial reuse of organic residues.

IMPORTANCE Microbiomes are vital to human health, agriculture, and protecting the environment. Predicting behavior of self-assembled or synthetic microbiomes, however, remains a challenge. In this work, we used unicellular and guild-based metabolic models to investigate production of medium-chain fatty acids by a mixed microbial community that is fed multiple organic substrates. Modeling results provided insights into metabolic pathways of three medium-chain fatty acid-producing guilds and identified potential strategies to increase production of medium-chain fatty acids. This work demonstrates the role of metabolic models in augmenting multi-omic studies to gain greater insights into microbiome behavior.

Information Theory in Computational Biology: Where We Stand Today

"A Mathematical Theory of Communication" was published in 1948 by Claude Shannon to address the problems in the field of data compression and communication over (noisy) communication channels. Since then, the concepts and ideas developed in Shannon's work have formed the basis of information theory, a cornerstone of statistical learning and inference, and has been playing a key role in disciplines such as physics and thermodynamics, probability and statistics, computational sciences and biological sciences. In this article we review the basic information theory based concepts and describe their key applications in multiple major areas of research in computational biology-gene expression and transcriptomics, alignment-free sequence comparison, sequencing and error correction, genome-wide disease-gene association mapping, metabolic networks and metabolomics, and protein sequence, structure and interaction analysis.

A Systems-Based Approach for Cyanide Overproduction by Bacillus megaterium for Gold Bioleaching Enhancement

With the constant accumulation of electronic waste, extracting precious metals contained therein is becoming a major challenge for sustainable development. Bacillus megaterium is currently one of the microbes used for the production of cyanide, which is the main leaching agent for gold recovery. The present study aimed to propose a strategy for metabolic engineering of B. megaterium to overproduce cyanide, and thus ameliorate the bioleaching process. For this, we employed constraint-based modeling, running in silico simulations on iJA1121, the genome-scale metabolic model of B. megaterium DSM319. Flux balance analysis (FBA) was initially used to identify amino acids to be added to the culture medium. Considering cyanide as the desired product, we used growth-coupled methods, constrained minimal cut sets (cMCSs) and OptKnock to identify gene inactivation targets. To identify gene overexpression targets, flux scanning based on enforced objective flux (FSEOF) was performed. Further analysis was carried out on the identified targets to determine compounds with beneficial regulatory effects. We have proposed a chemical-defined medium for accelerating cyanide production on the basis of microplate assays to evaluate the components with the greatest improving effects. Accordingly, the cultivation of B. megaterium DSM319 in a chemically-defined medium with 5.56 mM glucose as the carbon source, and supplemented with 413 μM cysteine, led to the production of considerably increased amounts of cyanide. Bioleaching experiments were successfully performed in this medium to recover gold and copper from telecommunication printed circuit boards. The results of inductively coupled plasma (ICP) analysis confirmed that gold recovery peaked out at around 55% after 4 days, whereas copper recovery continued to increase for several more days, peaking out at around 85%. To further validate the bioleaching results, FESEM, XRD, FTIR, and EDAX mapping analyses were performed. We concluded that the proposed strategy represents a viable route for improving the performance of the bioleaching processes.

Alternative Strategies for Microbial Remediation of Pollutants via Synthetic Biology

Continuous contamination of the environment with xenobiotics and related recalcitrant compounds has emerged as a serious pollution threat. Bioremediation is the key to eliminating persistent contaminants from the environment. Traditional bioremediation processes show limitations, therefore it is necessary to discover new bioremediation technologies for better results. In this review we provide an outlook of alternative strategies for bioremediation via synthetic biology, including exploring the prerequisites for analysis of research data for developing synthetic biological models of microbial bioremediation. Moreover, cell coordination in synthetic microbial community, cell signaling, and quorum sensing as engineered for enhanced bioremediation strategies are described, along with promising gene editing tools for obtaining the host with target gene sequences responsible for the degradation of recalcitrant compounds. The synthetic genetic circuit and two-component regulatory system (TCRS)-based microbial biosensors for detection and bioremediation are also briefly explained. These developments are expected to increase the efficiency of bioremediation strategies for best results.

Multiscale Multiobjective Systems Analysis (MiMoSA): an advanced metabolic modeling framework for complex systems

In natural environments, cells live in complex communities and experience a high degree of heterogeneity internally and in the environment. Even in 'ideal' laboratory environments, cells can experience a high degree of heterogeneity in their environments. Unfortunately, most of the metabolic modeling approaches that are currently used assume ideal conditions and that each cell is identical, limiting their application to pure cultures in well-mixed vessels. Here we describe our development of Multiscale Multiobjective Systems Analysis (MiMoSA), a metabolic modeling approach that can track individual cells in both space and time, track the diffusion of nutrients and light and the interaction of cells with each other and the environment. As a proof-of concept study, we used MiMoSA to model the growth of Trichodesmium erythraeum, a filamentous diazotrophic cyanobacterium which has cells with two distinct metabolic modes. The use of MiMoSA significantly improves our ability to predictively model metabolic changes and phenotype in more complex cell cultures.

Genome-Scale Metabolic Network Models of Bacillus Species Suggest that Model Improvement is Necessary for Biotechnological Applications

Background

A genome-scale metabolic network model (GEM) is a mathematical representation of an organism’s metabolism. Today, GEMs are popular tools for computationally simulating the biotechnological processes and for predicting biochemical properties of (engineered) strains.

Objectives

In the present study, we have evaluated the predictive power of two GEMs, namely iBsu1103 (for Bacillus subtilis 168) and iMZ1055 (for Bacillus megaterium WSH002).

Materials and Methods

For comparing the predictive power of Bacillus subtilis and Bacillus megaterium GEMs, experimental data were obtained from previous wet-lab studies included in PubMed. By using these data, we set the environmental, stoichiometric and thermodynamic constraints on the models, and FBA is performed to predict the biomass production rate, and the values of other fluxes. For simulating experimental conditions in this study, COBRA toolbox was used.

Results

By using the wealth of data in the literature, we evaluated the accuracy of in silico simulations of these GEMs. Our results suggest that there are some errors in these two models which make them unreliable for predicting the biochemical capabilities of these species. The inconsistencies between experimental and computational data are even greater where B. subtilis and B. megaterium do not have similar phenotypes.

Conclusions

Our analysis suggests that literature-based improvement of genome-scale metabolic network models of the two Bacillus species is essential if these models are to be successfully applied in biotechnology and metabolic engineering.

Boolean network analysis through the joint use of linear algebra and algebraic geometry

Among the various phenomena that can be modeled by Boolean networks, i.e., discrete-time dynamical systems with binary state variables, gene regulatory interactions are especially well known. Therefore, the analysis of Boolean networks is critical, e.g., to identify genetic pathways and to predict the effects of mutations on the cell functionality. Two methodologies (i.e., the semi-tensor product and the Gröbner bases over finite fields) have recently been proposed to tackle the problem of determining cycles and attractors (with the corresponding basin of attraction) for such systems. Here, it is shown that, by suitably coupling methodologies taken from these two fields (i.e., linear algebra and algebraic geometry), it is not only possible to determine cycles and attractors, but also to find closed-form solutions of the Boolean network. Such a goal is pursued by finding an immersion that recasts the Boolean dynamics in a linear form and by computing the closed-form solution of the latter system. The effectiveness of this technique is demonstrated by fully computing the solutions of the Boolean network modeling the differentiation of the Th-lymphocyte, a type of white blood cells involved in the human adaptive immune system.

Potential Natural Food Preservatives and Their Sustainable Production in Yeast: Terpenoids and Polyphenols

Terpenoids and polyphenols are high-valued plant secondary metabolites. Their high antimicrobial activities demonstrate their huge potential as natural preservatives in the food industry. With the rapid development of metabolic engineering, it has become possible to realize large-scale production of non-native terpenoids and polyphenols by using the generally recognized as safe (GRAS) strain, Saccharomyces cerevisiae, as a cell factory. This review will summarize the major terpenoid and polyphenol compounds with high antimicrobial properties, describe their native metabolic pathways as well as antimicrobial mechanisms, and highlight current progress on their heterologous biosynthesis in S. cerevisiae. Current challenges and perspectives for the sustainable production of terpenoid and polyphenol as natural food preservatives via S. cerevisiae will also be discussed.

Bacillus velezensis: phylogeny, useful applications, and avenues for exploitation

Some members of the Bacillus velezensis (Bv) group (e.g., Bv FZB42T and AS3.43) were previously assigned grouping with B. subtilis and B. amyloliquefaciens, based on the fact that they shared a 99% DNA-DNA percentage phylogenetic similarity. However, hinging on current assessments of the pan-genomic reassignments, the differing phylogenomic characteristics of Bv from B. subtilis and B. amyloliquefaciens are now better understood. Within this re-grouping/reassignment, the various strains within the Bv share a close phylogenomic resemblance, and a number of these strains have received a lot of attention in recent years, due to their genomic robustness, and the growing evidence for their possible utilization in the agricultural industry for managing plant diseases. Only a few applications for their use medicinally/pharmaceutically, environmentally, and in the food industry have been reported, and this may be due to the fact that the majority of those strains investigated are those typically occurring in soil. Although the intracellular unique biomolecules of Bv strains have been revealed via in silico genome modeling and investigated using transcriptomics and proteomics, a further inquisition into the Bv metabolome using newer technologies such as metabolomics could elucidate additional applications of this economically relevant Bacillus species, beyond that of primarily the agricultural sector.

Quantitative flux coupling analysis

Flux coupling analysis (FCA) aims to describe the functional dependencies among reactions in a metabolic network. Currently studied coupling relations are qualitative in the sense that they identify pairs of reactions for which the activity of one reaction necessitates the activity of the other one, but without giving any numerical bounds relating the possible activity rates. The potential applications of FCA are heavily investigated, however apart from some trivial cases there is no clue of what bottleneck in the metabolic network causes each dependency. In this article, we introduce a quantitative approach to the same flux coupling problem named quantitative flux coupling analysis (QFCA). It generalizes the current concepts as we show that all the qualitative information provided by FCA is readily available in the quantitative flux coupling equations of QFCA, without the need for any additional analysis. Moreover, we design a simple algorithm to efficiently identify these flux coupling equations which scales up to the genome-scale metabolic networks with thousands of reactions and metabolites in an effective way. Furthermore, this framework enables us to quantify the "strength" of the flux coupling relations. We also provide different biologically meaningful interpretations, including one which gives an intuitive certificate of precisely which metabolites in the network enforce each flux coupling relation. Eventually, we conclude by suggesting the probable application of QFCA to the metabolic gap-filling problem, which we only begin to address here and is left for future research to further investigate.

Metabolic modeling and response surface analysis of an Escherichia coli strain engineered for shikimic acid production

Background

Classic metabolic engineering strategies often induce significant flux imbalances to microbial metabolism, causing undesirable outcomes such as suboptimal conversion of substrates to products. Several mathematical frameworks have been developed to understand the physiological and metabolic state of production strains and to identify genetic modification targets for improved bioproduct formation. In this work, a modeling approach was applied to describe the physiological behavior and the metabolic fluxes of a shikimic acid overproducing Escherichia coli strain lacking the major glucose transport system, grown on complex media.

Results

The obtained flux distributions indicate the presence of high fluxes through the pentose phosphate and Entner-Doudoroff pathways, which could limit the availability of erythrose-4-phosphate for shikimic acid production even with high flux redirection through the pentose phosphate pathway. In addition, highly active glyoxylate shunt fluxes and a pyruvate/acetate cycle are indicators of overflow glycolytic metabolism in the tested conditions. The analysis of the combined physiological and flux response surfaces, enabled zone allocation for different physiological outputs within variant substrate conditions. This information was then used for an improved fed-batch process designed to preserve the metabolic conditions that were found to enhance shikimic acid productivity. This resulted in a 40% increase in the shikimic acid titer (60 g/L) and 70% increase in volumetric productivity (2.45 gSA/L*h), while preserving yields, compared to the batch process.

Conclusions

The combination of dynamic metabolic modeling and experimental parameter response surfaces was a successful approach to understand and predict the behavior of a shikimic acid producing strain under variable substrate concentrations. Response surfaces were useful for allocating different physiological behavior zones with different preferential product outcomes. Both model sets provided information that could be applied to enhance shikimic acid production on an engineered shikimic acid overproducing Escherichia coli strain.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0632-4) contains supplementary material, which is available to authorized users.

A systematic reconstruction and constraint-based analysis of Leishmania donovani metabolic network: identification of potential antileishmanial drug targets

Visceral leishmaniasis, a lethal parasitic disease, is caused by the protozoan parasite Leishmania donovani. The absence of an effective vaccine, drug toxicity and parasite resistance necessitates the identification of novel drug targets. Reconstruction of genome-scale metabolic models and their simulation has been established as an important tool for systems-level understanding of a microorganism's metabolism. In this work, amalgamating the tools and techniques of computational systems biology with rigorous manual curation, a constraint-based metabolic model for Leishmania donovani BPK282A1 has been developed. New functional annotations for 18 formerly hypothetical or erroneously annotated genes (encountered during iterative refinement of the model) have been proposed. Further, to formulate an accurate biomass objective function, experimental determination of previously uncharacterized biomass constituents was performed. The developed model is a highly compartmentalized metabolic model, comprising 1159 reactions, 1135 metabolites and 604 genes. The model exhibited around 76% accuracy for the prediction of experimental phenotypes of gene knockout studies and drug inhibition assays. Employing in silico gene knockout studies, we identified 28 essential genes with negligible sequence identity to the human proteins. Moreover, by dissecting the functional interdependencies of metabolic pathways, 70 synthetic lethal pairs were identified. Finally, in order to delineate stage-specific metabolism, gene-expression data of the amastigote stage residing in human macrophages were integrated into the model. By comparing the flux distribution, we illustrated the stage-specific differences in metabolism and environmental conditions that are in good agreement with the experimental findings. The developed model can serve as a highly enriched knowledgebase of legacy data and an important tool for generating experimentally verifiable hypotheses.

Genome-scale modeling of yeast: chronology, applications and critical perspectives

Over the last 15 years, several genome-scale metabolic models (GSMMs) were developed for different yeast species, aiding both the elucidation of new biological processes and the shift toward a bio-based economy, through the design of in silico inspired cell factories. Here, an historical perspective of the GSMMs built over time for several yeast species is presented and the main inheritance patterns among the metabolic reconstructions are highlighted. We additionally provide a critical perspective on the overall genome-scale modeling procedure, underlining incomplete model validation and evaluation approaches and the quest for the integration of regulatory and kinetic information into yeast GSMMs. A summary of experimentally validated model-based metabolic engineering applications of yeast species is further emphasized, while the main challenges and future perspectives for the field are finally addressed.

Integrative analysis of human omics data using biomolecular networks

High-throughput '-omics' technologies have given rise to an increasing abundance of genome-scale data detailing human biology at the molecular level. Although these datasets have already made substantial contributions to a more comprehensive understanding of human physiology and diseases, their interpretation becomes increasingly cryptic and nontrivial as they continue to expand in size and complexity. Systems biology networks offer a scaffold upon which omics data can be integrated, facilitating the extraction of new and physiologically relevant information from the data. Two of the most prevalent networks that have been used for such integrative analyses of omics data are genome-scale metabolic models (GEMs) and protein-protein interaction (PPI) networks, both of which have demonstrated success among many different omics and sample types. This integrative approach seeks to unite 'top-down' omics data with 'bottom-up' biological networks in a synergistic fashion that draws on the strengths of both strategies. As the volume and resolution of high-throughput omics data continue to grow, integrative network-based analyses are expected to play an increasingly important role in their interpretation.

Genome scale metabolic reconstruction of Chlorella variabilis for exploring its metabolic potential for biofuels

A compartmentalized genome scale metabolic network was reconstructed for Chlorella variabilis to offer insight into various metabolic potentials from this alga. The model, iAJ526, was reconstructed with 1455 reactions, 1236 metabolites and 526 genes. 21% of the reactions were transport reactions and about 81% of the total reactions were associated with enzymes. Along with gap filling reactions, 2 major sub-pathways were added to the model, chitosan synthesis and rhamnose metabolism. The reconstructed model had reaction participation of 4.3 metabolites per reaction and average lethality fraction of 0.21. The model was effective in capturing the growth of C. variabilis under three light conditions (white, red and red+blue light) with fair agreement. This reconstructed metabolic network will serve an important role in systems biology for further exploration of metabolism for specific target metabolites and enable improved characteristics in the strain through metabolic engineering.

Differential RNA-seq, Multi-Network Analysis and Metabolic Regulation Analysis of Kluyveromyces marxianus Reveals a Compartmentalised Response to Xylose

We investigated the transcriptomic response of a new strain of the yeast Kluyveromyces marxianus, in glucose and xylose media using RNA-seq. The data were explored in a number of innovative ways using a variety of networks types, pathway maps, enrichment statistics, reporter metabolites and a flux simulation model, revealing different aspects of the genome-scale response in an integrative systems biology manner. The importance of the subcellular localisation in the transcriptomic response is emphasised here, revealing new insights. As was previously reported by others using a rich medium, we show that peroxisomal fatty acid catabolism was dramatically up-regulated in a defined xylose mineral medium without fatty acids, along with mechanisms to activate fatty acids and transfer products of β-oxidation to the mitochondria. Notably, we observed a strong up-regulation of the 2-methylcitrate pathway, supporting capacity for odd-chain fatty acid catabolism. Next we asked which pathways would respond to the additional requirement for NADPH for xylose utilisation, and rationalised the unexpected results using simulations with Flux Balance Analysis. On a fundamental level, we investigated the contribution of the hierarchical and metabolic regulation levels to the regulation of metabolic fluxes. Metabolic regulation analysis suggested that genetic level regulation plays a major role in regulating metabolic fluxes in adaptation to xylose, even for the high capacity reactions, which is unexpected. In addition, isozyme switching may play an important role in re-routing of metabolic fluxes in subcellular compartments in K. marxianus.

Designing minimal microbial strains of desired functionality using a genetic algorithm

Background

The rational, in silico prediction of gene-knockouts to turn organisms into efficient cell factories is an essential and computationally challenging task in metabolic engineering. Elementary flux mode analysis in combination with constraint minimal cut sets is a particularly powerful method to identify optimal engineering targets, which will force an organism into the desired metabolic state. Given an engineering objective, it is theoretically possible, although computationally impractical, to find the best minimal intervention strategies.

Results

We developed a genetic algorithm (GA-MCS) to quickly find many (near) optimal intervention strategies while overcoming the above mentioned computational burden. We tested our algorithm on Escherichia coli metabolic networks of three different sizes to find intervention strategies satisfying three different engineering objectives.

Conclusions

We show that GA-MCS finds all practically relevant targets for any (non)-linear engineering objective. Our algorithm also found solutions comparable to previously published results. We show that for large networks optimal solutions are found within a fraction of the time used for a complete enumeration.

Shikimic Acid Production in Escherichia coli: From Classical Metabolic Engineering Strategies to Omics Applied to Improve Its Production

Shikimic acid (SA) is an intermediate of the SA pathway that is present in bacteria and plants. SA has gained great interest because it is a precursor in the synthesis of the drug oseltamivir phosphate (OSF), an efficient inhibitor of the neuraminidase enzyme of diverse seasonal influenza viruses, the avian influenza virus H5N1, and the human influenza virus H1N1. For the purposes of OSF production, SA is extracted from the pods of Chinese star anise plants (Illicium spp.), yielding up to 17% of SA (dry basis content). The high demand for OSF necessary to manage a major influenza outbreak is not adequately met by industrial production using SA from plants sources. As the SA pathway is present in the model bacteria Escherichia coli, several "intuitive" metabolically engineered strains have been applied for its successful overproduction by biotechnological processes, resulting in strains producing up to 71 g/L of SA, with high conversion yields of up to 0.42 (mol SA/mol Glc), in both batch and fed-batch cultures using complex fermentation broths, including glucose as a carbon source and yeast extract. Global transcriptomic analyses have been performed in SA-producing strains, resulting in the identification of possible key target genes for the design of a rational strain improvement strategy. Because possible target genes are involved in the transport, catabolism, and interconversion of different carbon sources and metabolic intermediates outside the central carbon metabolism and SA pathways, as genes involved in diverse cellular stress responses, the development of rational cellular strain improvement strategies based on omics data constitutes a challenging task to improve SA production in currently overproducing engineered strains. In this review, we discuss the main metabolic engineering strategies that have been applied for the development of efficient SA-producing strains, as the perspective of omics analysis has focused on further strain improvement for the production of this valuable aromatic intermediate.

White biotechnology: State of the art strategies for the development of biocatalysts for biorefining

White biotechnology is a term that is now often used to describe the implementation of biotechnology in the industrial sphere. Biocatalysts (enzymes and microorganisms) are the key tools of white biotechnology, which is considered to be one of the key technological drivers for the growing bioeconomy. Biocatalysts are already present in sectors such as the chemical and agro-food industries, and are used to manufacture products as diverse as antibiotics, paper pulp, bread or advanced polymers. This review proposes an original and global overview of highly complementary fields of biotechnology at both enzyme and microorganism level. A certain number of state of the art approaches that are now being used to improve the industrial fitness of biocatalysts particularly focused on the biorefinery sector are presented. The first part deals with the technologies that underpin the development of industrial biocatalysts, notably the discovery of new enzymes and enzyme improvement using directed evolution techniques. The second part describes the toolbox available by the cell engineer to shape the metabolism of microorganisms. And finally the last part focuses on the 'omic' technologies that are vital for understanding and guide microbial engineering toward more efficient microbial biocatalysts. Altogether, these techniques and strategies will undoubtedly help to achieve the challenging task of developing consolidated bioprocessing (i.e. CBP) readily available for industrial purpose.

Systems and photosystems: cellular limits of autotrophic productivity in cyanobacteria

Recent advances in the modeling of microbial growth and metabolism have shown that growth rate critically depends upon the optimal allocation of finite proteomic resources among different cellular functions and that modeling growth rates becomes more realistic with the explicit accounting for the costs of macromolecular synthesis, most importantly, protein expression. The "proteomic constraint" is considered together with its application to understanding photosynthetic microbial growth. The central hypothesis is that physical limits of cellular space (and corresponding solvation capacity) in conjunction with cell surface-to-volume ratios represent the underlying constraints on the maximal rate of autotrophic microbial growth. The limitation of cellular space thus constrains the size the total complement of macromolecules, dissolved ions, and metabolites. To a first approximation, the upper limit in the cellular amount of the total proteome is bounded this space limit. This predicts that adaptation to osmotic stress will result in lower maximal growth rates due to decreased cellular concentrations of core metabolic proteins necessary for cell growth owing the accumulation of compatible osmolytes, as surmised previously. The finite capacity of membrane and cytoplasmic space also leads to the hypothesis that the species-specific differences in maximal growth rates likely reflect differences in the allocation of space to niche-specific proteins with the corresponding diminution of space devoted to other functions including proteins of core autotrophic metabolism, which drive cell reproduction. An optimization model for autotrophic microbial growth, the autotrophic replicator model, was developed based upon previous work investigating heterotrophic growth. The present model describes autotrophic growth in terms of the allocation protein resources among core functional groups including the photosynthetic electron transport chain, light-harvesting antennae, and the ribosome groups.

An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models

Constraint-based models are currently the only methodology that allows the study of metabolism at the whole-genome scale. Flux balance analysis is commonly used to analyse constraint-based models. Curiously, the results of this analysis vary with the software being run, a situation that we show can be remedied by using exact rather than floating-point arithmetic. Here we introduce MONGOOSE, a toolbox for analysing the structure of constraint-based metabolic models in exact arithmetic. We apply MONGOOSE to the analysis of 98 existing metabolic network models and find that the biomass reaction is surprisingly blocked (unable to sustain non-zero flux) in nearly half of them. We propose a principled approach for unblocking these reactions and extend it to the problems of identifying essential and synthetic lethal reactions and minimal media. Our structural insights enable a systematic study of constraint-based metabolic models, yielding a deeper understanding of their possibilities and limitations.

Current tools to analyse constraint-based models of metabolic networks have limited accuracy due to their use of floating-point arithmetic. Here the authors present MONGOOSE, a new computational tool that analyses such models in exact arithmetic, providing improved accuracy and reproducibility.

[Mechanisms of electron transfer to insoluble terminal acceptors in chemoorganotrophic bacteria]

The mechanisms of electron transfer of association of chemoorganotrophic bacteria to the anode in microbial fuel cells are summarized in the survey. These mechanisms are not mutually exclusive and are divided into the mechanisms of mediator electron transfer, mechanisms of electron transfer with intermediate products of bacterial metabolism and mechanism of direct transfer of electrons from the cell surface. Thus, electron transfer mediators are artificial or synthesized by bacteria riboflavins and phenazine derivatives, which also determine the ability of bacteria to antagonism. The microorganisms with hydrolytic and exoelectrogenic activity are involved in electron transfer mechanisms that are mediated by intermediate metabolic products, which are low molecular carboxylic acids, alcohols, hydrogen etc. The direct transfer of electrons to insoluble anode is possible due to membrane structures (cytochromes, pili, etc.). Association of microorganisms, and thus the biochemical mechanisms of electron transfer depend on the origin of the inoculum, substrate composition, mass transfer, conditions of aeration, potentials and location of electrodes and others, that are defined by technological and design parameters.

Genome scale modeling in systems biology: algorithms and resources

In recent years, in silico studies and trial simulations have complemented experimental procedures. A model is a description of a system, and a system is any collection of interrelated objects; an object, moreover, is some elemental unit upon which observations can be made but whose internal structure either does not exist or is ignored. Therefore, any network analysis approach is critical for successful quantitative modeling of biological systems. This review highlights some of most popular and important modeling algorithms, tools, and emerging standards for representing, simulating and analyzing cellular networks in five sections. Also, we try to show these concepts by means of simple example and proper images and graphs. Overall, systems biology aims for a holistic description and understanding of biological processes by an integration of analytical experimental approaches along with synthetic computational models. In fact, biological networks have been developed as a platform for integrating information from high to low-throughput experiments for the analysis of biological systems. We provide an overview of all processes used in modeling and simulating biological networks in such a way that they can become easily understandable for researchers with both biological and mathematical backgrounds. Consequently, given the complexity of generated experimental data and cellular networks, it is no surprise that researchers have turned to computer simulation and the development of more theory-based approaches to augment and assist in the development of a fully quantitative understanding of cellular dynamics.

An efficient graph theory based method to identify every minimal reaction set in a metabolic network

Background

Development of cells with minimal metabolic functionality is gaining importance due to their efficiency in producing chemicals and fuels. Existing computational methods to identify minimal reaction sets in metabolic networks are computationally expensive. Further, they identify only one of the several possible minimal reaction sets.

Results

In this paper, we propose an efficient graph theory based recursive optimization approach to identify all minimal reaction sets. Graph theoretical insights offer systematic methods to not only reduce the number of variables in math programming and increase its computational efficiency, but also provide efficient ways to find multiple optimal solutions. The efficacy of the proposed approach is demonstrated using case studies from Escherichia coli and Saccharomyces cerevisiae. In case study 1, the proposed method identified three minimal reaction sets each containing 38 reactions in Escherichia coli central metabolic network with 77 reactions. Analysis of these three minimal reaction sets revealed that one of them is more suitable for developing minimal metabolism cell compared to other two due to practically achievable internal flux distribution. In case study 2, the proposed method identified 256 minimal reaction sets from the Saccharomyces cerevisiae genome scale metabolic network with 620 reactions. The proposed method required only 4.5 hours to identify all the 256 minimal reaction sets and has shown a significant reduction (approximately 80%) in the solution time when compared to the existing methods for finding minimal reaction set.

Conclusions

Identification of all minimal reactions sets in metabolic networks is essential since different minimal reaction sets have different properties that effect the bioprocess development. The proposed method correctly identified all minimal reaction sets in a both the case studies. The proposed method is computationally efficient compared to other methods for finding minimal reaction sets and useful to employ with genome-scale metabolic networks.

Comparison on extreme pathways reveals nature of different biological processes

Background

Constraint-based reconstruction and analysis (COBRA) is used for modeling genome-scale metabolic networks (MNs). In a COBRA model, extreme pathways (ExPas) are the edges of its conical solution space, which is formed by all viable steady-state flux distributions. ExPa analysis has been successfully applied to MNs to reveal their phenotypic capabilities and properties. Recently, the COBRA framework has been extended to transcriptional regulatory networks (TRNs) and transcriptional and translational networks (TTNs), so efforts are needed to determine whether ExPa analysis is also effective on these two types of networks.

Results

In this paper, the ExPas resulting from the COBRA models of E.coli's MN, TRN and TTN were horizontally compared from 5 aspects: (1) Total number and the ratio of their amount to reaction amount; (2) Length distribution; (3) Reaction participation; (4) Correlated reaction sets (CoSets); (5) interconnectivity degree. Significant discrepancies in above properties were observed during the comparison, which reveals the biological natures of different biological processes. Besides, by demonstrating the application of ExPa analysis on E.coli, we provide a practical guidance of an improved approach to compute ExPas on COBRA models of TRNs.

Conclusions

ExPas of E.coli's MN, TRN and TTN have different properties, which are strongly connected with various biological natures of biochemical networks, such as topological structure, specificity and redundancy. Our study shows that ExPas are biologically meaningful on the newborn models and suggests the effectiveness of ExPa analysis on them.

Systems medicine: a new approach to clinical practice

Most respiratory diseases are considered complex diseases as their susceptibility and outcomes are determined by the interaction between host-dependent factors (genetic factors, comorbidities, etc.) and environmental factors (exposure to microorganisms or allergens, treatments received, etc.) The reductionist approach in the study of diseases has been of fundamental importance for the understanding of the different components of a system. Systems biology or systems medicine is a complementary approach aimed at analyzing the interactions between the different components within one organizational level (genome, transcriptome, proteome), and then between the different levels. Systems medicine is currently used for the interpretation and understanding of the pathogenesis and pathophysiology of different diseases, biomarker discovery, design of innovative therapeutic targets, and the drawing up of computational models for different biological processes. In this review we discuss the most relevant concepts of the theory underlying systems medicine, as well as its applications in the various biological processes in humans.

Plant genome-scale modeling and implementation

Considerable progress has been made in plant genome-scale metabolic reconstruction and modeling in recent years. Such reconstructions made it possible to explore metabolic phenotypes through appropriate model formulation and optimization methods. As a result, plant genome-scale modeling has increasingly attracted interest from the plant research community. In this chapter, the first generation of plant genome-scale metabolic reconstructions is presented, along with the important concepts behind model and constraint formulation. A brief protocol describing the use of constraint-based reconstruction and analysis (COBRA) Toolbox in flux simulation and model modification is provided. This is followed by a presentation of metabolic constraints required to generate fluxes in AraGEM using COBRA that describe photosynthesis, photorespiration, and respiration, respectively. Overall, plant genome-scale modeling is a powerful approach that is accessible and readily adopted.

Incorporation of flexible objectives and time-linked simulation with flux balance analysis

We present two modifications of the flux balance analysis (FBA) metabolic modeling framework which relax implicit assumptions of the biomass reaction. Our flexible flux balance analysis (flexFBA) objective removes the fixed proportion between reactants, and can therefore produce a subset of biomass reactants. Our time-linked flux balance analysis (tFBA) simulation removes the fixed proportion between reactants and byproducts, and can therefore describe transitions between metabolic steady states. Used together, flexFBA and tFBA model a time scale shorter than the regulatory and growth steady state encoded by the biomass reaction. This combined short-time FBA method is intended for integrated modeling applications to enable detailed and dynamic depictions of microbial physiology such as whole-cell modeling. For example, when modeling Escherichia coli, it avoids artifacts caused by low-copy-number enzymes in single-cell models with kinetic bounds. Even outside integrated modeling contexts, the detailed predictions of flexFBA and tFBA complement existing FBA techniques. We show detailed metabolite production of in silico knockouts used to identify when correct essentiality predictions are made for the wrong reason.

Semi-automated curation of metabolic models via flux balance analysis: a case study with Mycoplasma gallisepticum

Primarily used for metabolic engineering and synthetic biology, genome-scale metabolic modeling shows tremendous potential as a tool for fundamental research and curation of metabolism. Through a novel integration of flux balance analysis and genetic algorithms, a strategy to curate metabolic networks and facilitate identification of metabolic pathways that may not be directly inferable solely from genome annotation was developed. Specifically, metabolites involved in unknown reactions can be determined, and potentially erroneous pathways can be identified. The procedure developed allows for new fundamental insight into metabolism, as well as acting as a semi-automated curation methodology for genome-scale metabolic modeling. To validate the methodology, a genome-scale metabolic model for the bacterium Mycoplasma gallisepticum was created. Several reactions not predicted by the genome annotation were postulated and validated via the literature. The model predicted an average growth rate of 0.358±0.12, closely matching the experimentally determined growth rate of M. gallisepticum of 0.244±0.03. This work presents a powerful algorithm for facilitating the identification and curation of previously known and new metabolic pathways, as well as presenting the first genome-scale reconstruction of M. gallisepticum.

Flux balance analysis (FBA) is a powerful approach for genome-scale metabolic modeling. It provides metabolic engineers with a tool for manipulating, predicting, and optimizing metabolism for biotechnological and biomedical purposes. However, we posit that it can also be used as tool for fundamental research in understanding and curating metabolic networks. Specifically, by using a genetic algorithm integrated with FBA, we developed a curation approach to identify missing reactions, incomplete reactions, and erroneous reactions. Additionally, it was possible to take advantage of the ensemble information from the genetic algorithm to identify the most critical reactions for curation. We tested our strategy using Mycoplasma gallisepticum as our model organism. Using the genome annotation as the basis, the preliminary genome-scale metabolic model consisted of 446 metabolites involved in 380 reactions. Carrying out our analysis, we found over 80 incorrect reactions and 16 missing reactions. Based upon the guidance of the algorithm, we were able to curate and resolve all discrepancies. The model predicted an average bacterial growth rate of 0.358±0.12 h⁻¹ compared to the experimentally observed 0.244±0.03 h⁻¹. Thus, our approach facilitated the curation of a genome-scale metabolic network and generated a high quality metabolic model.

From top-down to bottom-up: computational modeling approaches for cellular redoxin networks

SIGNIFICANCE

Thioredoxin, glutaredoxin, and peroxiredoxin systems play critical roles in a large number of redox-sensitive cellular processes. These systems are linked to each other by coupled redox cycles and common reaction intermediates into a larger network. Given the scale and connectivity of this network, computational approaches are required to analyze its dynamics and organization.

RECENT ADVANCES

Theoretical advances, as well as new redox proteomic methods, have led to the development of both top-down and bottom-up systems biology approaches to analyze the these systems and the network as a whole. Top-down approaches have been based on modifications to the Nernst equation or on graph theoretical approaches, while bottom-up approaches have been based on kinetic or stoichiometric modeling techniques.

CRITICAL ISSUES

This review will consider the rationale behind these approaches and focus on their advantages and limitations. Further, the review will discuss modeling standards to ensure model accuracy and availability.

FUTURE DIRECTIONS

Top-down and bottom-up approaches have distinct strengths and limitations in describing cellular redoxin networks. The availability of methods to overcome these limitations, together with the adoption of common modeling standards, is expected to increase the pace of model-led discovery within the redox biology field.

Molecular tools for chemical biotechnology

Biotechnological production of high value chemical products increasingly involves engineering in vivo multi-enzyme pathways and host metabolism. Recent approaches to these engineering objectives have made use of molecular tools to advance de novo pathway identification, tunable enzyme expression, and rapid pathway construction. Molecular tools also enable optimization of single enzymes and entire genomes through diversity generation and screening, whole cell analytics, and synthetic metabolic control networks. In this review, we focus on advanced molecular tools and their applications to engineered pathways in host organisms, highlighting the degree to which each tool is generalizable.

Bioprocess systems engineering: transferring traditional process engineering principles to industrial biotechnology

The complexity of the regulatory network and the interactions that occur in the intracellular environment of microorganisms highlight the importance in developing tractable mechanistic models of cellular functions and systematic approaches for modelling biological systems. To this end, the existing process systems engineering approaches can serve as a vehicle for understanding, integrating and designing biological systems and processes. Here, we review the application of a holistic approach for the development of mathematical models of biological systems, from the initial conception of the model to its final application in model-based control and optimisation. We also discuss the use of mechanistic models that account for gene regulation, in an attempt to advance the empirical expressions traditionally used to describe micro-organism growth kinetics, and we highlight current and future challenges in mathematical biology. The modelling research framework discussed herein could prove beneficial for the design of optimal bioprocesses, employing rational and feasible approaches towards the efficient production of chemicals and pharmaceuticals.