Pichia pastoris growth-coupled heme biosynthesis analysis using metabolic modelling.

Soy leghemoglobin is one of the most important and key ingredients in plant-based meat substitutes that can imitate the colour and flavour of the meat. To improve the high-yield production of leghemoglobin protein and its main component-heme in the yeast Pichia pastoris, glycerol and methanol cultivation conditions were studied. Additionally, in-silico metabolic modelling analysis of growth-coupled enzyme quantity, suggests metabolic gene up/down-regulation strategies for heme production. First, cultivations and metabolic modelling analysis of P. pastoris were performed on glycerol and methanol in different growth media. Glycerol cultivation uptake and production rates can be increased by 50% according to metabolic modelling results, but methanol cultivation-is near the theoretical maximum. Growth-coupled metabolic optimisation results revealed the best feasible upregulation (33 reactions) (1.47% of total reactions) and 66 downregulation/deletion (2.98% of total) reaction suggestions. Finally, we describe reaction regulation suggestions with the highest potential to increase heme production yields.

Metabolic Modeling of Hermetia illucens Larvae Resource Allocation for High-Value Fatty Acid Production.

All plant and animal kingdom organisms use highly connected biochemical networks to facilitate sustaining, proliferation, and growth functions. While the biochemical network details are well known, the understanding of the intense regulation principles is still limited. We chose to investigate the Hermetia illucens fly at the larval stage because this stage is a crucial period for the successful accumulation and allocation of resources for the subsequent organism's developmental stages. We combined iterative wet lab experiments and innovative metabolic modeling design approaches to simulate and explain the H. illucens larval stage resource allocation processes and biotechnology potential. We performed time-based growth and high-value chemical compound accumulation wet lab chemical analysis experiments on larvae and the Gainesville diet composition. We built and validated the first H. illucens medium-size, stoichiometric metabolic model to predict the effects of diet-based alterations on fatty acid allocation potential. Using optimization methods such as flux balance and flux variability analysis on the novel insect metabolic model, we predicted that doubled essential amino acid consumption increased the growth rate by 32%, but pure glucose consumption had no positive impact on growth. In the case of doubled pure valine consumption, the model predicted a 2% higher growth rate. In this study, we describe a new framework for researching the impact of dietary alterations on the metabolism of multi-cellular organisms at different developmental stages for improved, sustainable, and directed high-value chemicals.

Integrative Gene Expression and Metabolic Analysis Tool IgemRNA.

Genome-scale metabolic modeling is widely used to study the impact of metabolism on the phenotype of different organisms. While substrate modeling reflects the potential distribution of carbon and other chemical elements within the model, the additional use of omics data, e.g., transcriptome, has implications when researching the genotype-phenotype responses to environmental changes. Several algorithms for transcriptome analysis using genome-scale metabolic modeling have been proposed. Still, they are restricted to specific objectives and conditions and lack flexibility, have software compatibility issues, and require advanced user skills. We classified previously published algorithms, summarized transcriptome pre-processing, integration, and analysis methods, and implemented them in the newly developed transcriptome analysis tool IgemRNA, which (1) has a user-friendly graphical interface, (2) tackles compatibility issues by combining previous data input and pre-processing algorithms in MATLAB, and (3) introduces novel algorithms for the automatic comparison of different transcriptome datasets with or without Cobra Toolbox 3.0 optimization algorithms. We used publicly available transcriptome datasets from Saccharomyces cerevisiae BY4741 and H4-S47D strains for validation. We found that IgemRNA provides a means for transcriptome and environmental data validation on biochemical network topology since the biomass function varies for different phenotypes. Our tool can detect problematic reaction constraints.

Kinetic and Stoichiometric Modeling-Based Analysis of Docosahexaenoic Acid (DHA) Production Potential by Crypthecodinium cohnii from Glycerol, Glucose and Ethanol.

Docosahexaenoic acid (DHA) is one of the most important long-chain polyunsaturated fatty acids (LC-PUFAs), with numerous health benefits. Crypthecodinium cohnii, a marine heterotrophic dinoflagellate, is successfully used for the industrial production of DHA because it can accumulate DHA at high concentrations within the cells. Glycerol is an interesting renewable substrate for DHA production since it is a by-product of biodiesel production and other industries, and is globally generated in large quantities. The DHA production potential from glycerol, ethanol and glucose is compared by combining fermentation experiments with the pathway-scale kinetic modeling and constraint-based stoichiometric modeling of C. cohnii metabolism. Glycerol has the slowest biomass growth rate among the tested substrates. This is partially compensated by the highest PUFAs fraction, where DHA is dominant. Mathematical modeling reveals that glycerol has the best experimentally observed carbon transformation rate into biomass, reaching the closest values to the theoretical upper limit. In addition to our observations, the published experimental evidence indicates that crude glycerol is readily consumed by C. cohnii, making glycerol an attractive substrate for DHA production.

IMFLer: A Web Application for Interactive Metabolic Flux Analysis and Visualization.

Increasing genome-wide data in biological sciences and medicine has contributed to the development of a variety of visualization tools. Several automatic, semiautomatic, and manual visualization tools have already been developed. Some even have integrated flux balance analysis (FBA), but in most cases, it depends on separately installed third party software that is proprietary and does not allow customization of its functionality and has many restrictions for easy data distribution and analysis. In this study, we present an interactive metabolic flux analyzer and visualizer (IMFLer)-a static single-page web application that enables the reading and management of metabolic model layout maps, as well as immediate visualization of results from both FBA and flux variability analysis (FVA). IMFLer uses the Escher Builder tool to load, show, edit, and save metabolic pathway maps. This makes IMFLer an attractive and easily applicable tool with a user-friendly interface. Moreover, it allows to faster interpret results from FBA and FVA and improves data interoperability by using a standardized file format for the genome-scale metabolic model. IMFLer is a fully open-source tool that enables the rapid visualization and interpretation of the results of FBA and FVA with no time setup and no programming skills required, available at https://lv-csbg.github.io/IMFLer/.

Improvement of Acetaldehyde Production in Zymomonas mobilis by Engineering of Its Aerobic Metabolism.

Acetaldehyde is a valuable product of microbial biosynthesis, which can be used by the chemical industry as the entry point for production of various commodity chemicals. In ethanologenic microorganisms, like yeast or the bacterium Zymomonas mobilis, this compound is the immediate metabolic precursor of ethanol. In aerobic cultures of Z. mobilis, it accumulates as a volatile, inhibitory byproduct, due to the withdrawal of reducing equivalents from the alcohol dehydrogenase reaction by respiration. The active respiratory chain of Z. mobilis with its low energy-coupling efficiency is well-suited for regeneration of NAD+ under conditions when acetaldehyde, but not ethanol, is the desired catabolic product. In the present work, we sought to improve the capacity Z. mobilis to synthesize acetaldehyde, based on predictions of a stoichiometric model of its central metabolism developed herein. According to the model analysis, the main objectives in the course of engineering acetaldehyde producer strains were determined to be: (i) reducing ethanol synthesis via reducing the activity of alcohol dehydrogenase (ADH), and (ii) enhancing the respiratory capacity, either by overexpression of the respiratory NADH dehydrogenase (NDH), or by mutation of other components of respiratory metabolism. Several mutants with elevated respiration rate, decreased alcohol dehydrogenase activity, or a combination of both, were obtained. They were extensively characterized by determining their growth rates, product yields, oxygen consumption rates, ADH, and NDH activities, transcription levels of key catabolic genes, as well as concentrations of central metabolites under aerobic culture conditions. Two mutant strains were selected, with acetaldehyde yield close to 70% of the theoretical maximum value, almost twice the previously published yield for Z. mobilis. These strains can serve as a basis for further development of industrial acetaldehyde producers.

Model-based metabolism design: constraints for kinetic and stoichiometric models.

The implementation of model-based designs in metabolic engineering and synthetic biology may fail. One of the reasons for this failure is that only a part of the real-world complexity is included in models. Still, some knowledge can be simplified and taken into account in the form of optimization constraints to improve the feasibility of model-based designs of metabolic pathways in organisms. Some constraints (mass balance, energy balance, and steady-state assumption) serve as a basis for many modelling approaches. There are others (total enzyme activity constraint and homeostatic constraint) proposed decades ago, but which are frequently ignored in design development. Several new approaches of cellular analysis have made possible the application of constraints like cell size, surface, and resource balance. Constraints for kinetic and stoichiometric models are grouped according to their applicability preconditions in (1) general constraints, (2) organism-level constraints, and (3) experiment-level constraints. General constraints are universal and are applicable for any system. Organism-level constraints are applicable for biological systems and usually are organism-specific, but these constraints can be applied without information about experimental conditions. To apply experimental-level constraints, peculiarities of the organism and the experimental set-up have to be taken into account to calculate the values of constraints. The limitations of applicability of particular constraints for kinetic and stoichiometric models are addressed.

Total enzyme activity constraint and homeostatic constraint impact on the optimization potential of a kinetic model.

The application of biologically and biochemically relevant constraints during the optimization of kinetic models reduces the impact of suggested changes in processes not included in the scope of the model. This increases the probability that the design suggested by model optimization can be carried out by an organism after implementation of design in vivo. A case study was carried out to determine the impact of total enzyme activity and homeostatic constraints on the objective function values and the following ranking of adjustable parameter combinations. The application of constraints on the model of sugar cane metabolism revealed that a homeostatic constraint caused heavier limitations of the objective function than a total enzyme activity constraint. Both constraints changed the ranking of adjustable parameter combinations: no "universal" constraint-independent top-ranked combinations were found. Therefore, when searching for the best subset of adjustable parameters, a full scan of their combinations is suggested for a small number of adjustable parameters, and evolutionary search strategies are suggested for a large number. Simultaneous application of both constraints is suggested.

SpaceScanner: COPASI wrapper for automated management of global stochastic optimization experiments.

MOTIVATION: Due to their universal applicability, global stochastic optimization methods are popular for designing improvements of biochemical networks. The drawbacks of global stochastic optimization methods are: (i) no guarantee of finding global optima, (ii) no clear optimization run termination criteria and (iii) no criteria to detect stagnation of an optimization run. The impact of these drawbacks can be partly compensated by manual work that becomes inefficient when the solution space is large due to combinatorial explosion of adjustable parameters or for other reasons. RESULTS: SpaceScanner uses parallel optimization runs for automatic termination of optimization tasks in case of consensus and consecutively applies a pre-defined set of global stochastic optimization methods in case of stagnation in the currently used method. Automatic scan of adjustable parameter combination subsets for best objective function values is possible with a summary file of ranked solutions. AVAILABILITY AND IMPLEMENTATION: https://github.com/atiselsts/spacescanner . CONTACT: egils.stalidzans@lu.lv. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Modeling of Zymomonas mobilis central metabolism for novel metabolic engineering strategies.

Mathematical modeling of metabolism is essential for rational metabolic engineering. The present work focuses on several types of modeling approach to quantitative understanding of central metabolic network and energetics in the bioethanol-producing bacterium Zymomonas mobilis. Combined use of Flux Balance, Elementary Flux Mode, and thermodynamic analysis of its central metabolism, together with dynamic modeling of the core catabolic pathways, can help to design novel substrate and product pathways by systematically analyzing the solution space for metabolic engineering, and yields insights into the function of metabolic network, hardly achievable without applying modeling tools.

Biotechnological potential of respiring Zymomonas mobilis: a stoichiometric analysis of its central metabolism.

The active, yet energetically inefficient electron transport chain of the ethanologenic bacterium Zymomonas mobilis could be used in metabolic engineering for redox-balancing purposes during synthesis of certain products. Although several reconstructions of Z. mobilis metabolism have been published, important aspects of redox balance and aerobic catabolism have not previously been considered. Here, annotated genome sequences and metabolic reconstructions have been combined with existing biochemical evidence to yield a medium-scale model of Z. mobilis central metabolism in the form of COBRA Toolbox model files for flux balance analysis (FBA). The stoichiometric analysis presented here suggests the feasibility of several metabolic engineering strategies for obtaining high-value products, such as glycerate, succinate, and glutamate that would use the electron transport chain to oxidize the excess NAD(P)H, generated during synthesis of these metabolites. Oxidation of the excess NAD(P)H would also be needed for synthesis of ethanol from glycerol. Maximum product yields and the byproduct spectra have been estimated for each product, with glucose, xylose, or glycerol as the carbon substrates. These novel pathways represent targets for future metabolic engineering, as they would exploit both the rapid Entner-Doudoroff glycolysis, and the energetically uncoupled electron transport of Z. mobilis.