Unconventional systems analysis problems in molecular biology: a case study in gene regulatory network modeling

A novel kinetic model for a cocoa waste fermentation to ethanol reaction and its experimental validation

A non-segregated kinetic model is proposed to describe a fermentation process of agro-industrial residues derived via cocoa (mucilage juice) by Pichia kudriavzevii. The novel proposed hybrid model is based on a multiple coupling reaction mechanisms (structured) to describe the kinetics of substrate consumption, biomass, carbon dioxide, and ethanol, coupled to an unstructured model for the activity enzyme. The parameters of the kinetic model are estimated by non-linear least-squares curve fitting using the Marquardt-Levenberg algorithm. In addition, numerical simulations were compared with the experimental data via residual graphs. The effectiveness of the model was statistically evaluated using dimensionless efficiency coefficients under different initial conditions. A global sensitivity analysis was applied (Fisher's information matrix). The experimental results of the batch reactor showed a maximum ethanol concentration of 29 g/L, with a yield of 0.48 g-ethanol/g-glucose and a productivity of 0.30 g/L h. The method determined that the cell formation coefficient and the specific substrate consumption rate (θ1 and θ2) directly influence most of the states of our system. The proposed scheme is particularly suitable to assist in the rational design of cell factory properties or fermentation processes because it can represent the complex biochemistry in more detail and under different initial experimental conditions; the above reveals that the generated model is robust and can be considered for control and optimization purposes.

Tryptophan Production Maximization in a Fed-Batch Bioreactor with Modified E. coli Cells, by Optimizing Its Operating Policy Based on an Extended Structured Cell Kinetic Model

Hybrid kinetic models, linking structured cell metabolic processes to the dynamics of macroscopic variables of the bioreactor, are more and more used in engineering evaluations to derive more precise predictions of the process dynamics under variable operating conditions. Depending on the cell model complexity, such a math tool can be used to evaluate the metabolic fluxes in relation to the bioreactor operating conditions, thus suggesting ways to genetically modify the microorganism for certain purposes. Even if development of such an extended dynamic model requires more experimental and computational efforts, its use is advantageous. The approached probative example refers to a model simulating the dynamics of nanoscale variables from several pathways of the central carbon metabolism (CCM) of Escherichia coli cells, linked to the macroscopic state variables of a fed-batch bioreactor (FBR) used for the tryptophan (TRP) production. The used E. coli strain was modified to replace the PTS system for glucose (GLC) uptake with a more efficient one. The study presents multiple elements of novelty: (i) the experimentally validated modular model itself, and (ii) its efficiency in computationally deriving an optimal operation policy of the FBR.

Homeostasis in a feed forward loop gene regulatory motif

The internal state of a cell is affected by inputs from the extra-cellular environment such as external temperature. If some output, such as the concentration of a target protein, remains approximately constant as inputs vary, the system exhibits homeostasis. Special sub-networks called motifs are unusually common in gene regulatory networks (GRNs), suggesting that they may have a significant biological function. Potentially, one such function is homeostasis. In support of this hypothesis, we show that the feed-forward loop GRN produces homeostasis. Here the inputs are subsumed into a single parameter that affects only the first node in the motif, and the output is the concentration of a target protein. The analysis uses the notion of infinitesimal homeostasis, which occurs when the input-output map has a critical point (zero derivative). In model equations such points can be located using implicit differentiation. If the second derivative of the input-output map also vanishes, the critical point is a chair: the output rises roughly linearly, then flattens out (the homeostasis region or plateau), and then starts to rise again. Chair points are a common cause of homeostasis. In more complicated equations or networks, numerical exploration would have to augment analysis. Thus, in terms of finding chairs, this paper presents a proof of concept. We apply this method to a standard family of differential equations modeling the feed-forward loop GRN, and deduce that chair points occur. This function determines the production of a particular mRNA and the resulting chair points are found analytically. The same method can potentially be used to find homeostasis regions in other GRNs. In the discussion and conclusion section, we also discuss why homeostasis in the motif may persist even when the rest of the network is taken into account.

An Evaluation of the Information Technology of Gene Expression Profiles Processing Stability for Different Levels of Noise Components

This paper presents the results of research concerning the evaluation of stability of information technology of gene expression profiles processing with the use of gene expression profiles, which contain different levels of noise components. The information technology is presented as a structural block-chart, which contains all stages of the studied data processing. The hybrid model of objective clustering based on the SOTA algorithm and the technology of gene regulatory networks reconstruction have been investigated to evaluate the stability to the level of the noise components. The results of the simulation have shown that the hybrid model of the objective clustering has high level of stability to noise components and vice versa, the technology of gene regulatory networks reconstruction is rather sensitive to the level of noise component. The obtained results indicate the importance of gene expression profiles preprocessing at the early stage of the gene regulatory network reconstruction in order to remove background noise and non-informative genes in terms of the used criteria.

Quantitative modeling and analytic assessment of the transcription dynamics of the XlnR regulon in Aspergillus niger

Background

Transcription of genes coding for xylanolytic and cellulolytic enzymes in Aspergillus niger is controlled by the transactivator XlnR. In this work we analyse and model the transcription dynamics in the XlnR regulon from time-course data of the messenger RNA levels for some XlnR target genes, obtained by reverse transcription quantitative PCR (RT-qPCR). Induction of transcription was achieved using low (1 mM) and high (50 mM) concentrations of D-xylose (Xyl). We investigated the wild type strain (Wt) and a mutant strain with partial loss-of-function of the carbon catabolite repressor CreA (Mt).

Results

An improved kinetic differential equation model based on two antagonistic Hill functions was proposed, and fitted to the time-course RT-qPCR data from the Wt and the Mt by numerical optimization of the parameters. We show that perturbing the XlnR regulon with Xyl in low and high concentrations results in different expression levels and transcription dynamics of the target genes. At least four distinct transcription profiles were observed, particularly for the usage of 50 mM Xyl. Higher transcript levels were observed for some genes after induction with 1 mM rather than 50 mM Xyl, especially in the Mt. Grouping the expression profiles of the investigated genes has improved our understanding of induction by Xyl and the according regulatory role of CreA.

Conclusions

The model explains for the higher expression levels at 1 mM versus 50 mM in both Wt and Mt. It does not yet fully encapsulate the effect of partial loss-of-function of CreA in the Mt. The model describes the dynamics in most of the data and elucidates the time-dynamics of the two major regulatory mechanisms: i) the activation by XlnR, and ii) the carbon catabolite repression by CreA.

Electronic supplementary material

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