Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community’s traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators

Kinetic Intersections as a Source of Information on Rate Coefficients

C

Regulation-Structured Dynamic Metabolic Model Provides a Potential Mechanism for Delayed Enzyme Response in Denitrification Process

In a recent study of denitrification dynamics in hyporheic zone sediments, we observed a significant time lag (up to several days) in enzymatic response to the changes in substrate concentration. To explore an underlying mechanism and understand the interactive dynamics between enzymes and nutrients, we developed a trait-based model that associates a community’s traits with functional enzymes, instead of typically used species guilds (or functional guilds). This enzyme-based formulation allows to collectively describe biogeochemical functions of microbial communities without directly parameterizing the dynamics of species guilds, therefore being scalable to complex communities. As a key component of modeling, we accounted for microbial regulation occurring through transcriptional and translational processes, the dynamics of which was parameterized based on the temporal profiles of enzyme concentrations measured using a new signature peptide-based method. The simulation results using the resulting model showed several days of a time lag in enzymatic responses as observed in experiments. Further, the model showed that the delayed enzymatic reactions could be primarily controlled by transcriptional responses and that the dynamics of transcripts and enzymes are closely correlated. The developed model can serve as a useful tool for predicting biogeochemical processes in natural environments, either independently or through integration with hydrologic flow simulators

Tropical geometries and dynamics of biochemical networks. Application to hybrid cell cycle models

We use the Litvinov-Maslov correspondence principle to reduce and hybridize networks of biochemical reactions. We apply this method to a cell cycle oscillator model. The reduced and hybridized model can be used as a hybrid model for the cell cycle. We also propose a practical recipe for detecting quasi-equilibrium QE reactions and quasi-steady state QSS species in biochemical models with rational rate functions and use this recipe for model reduction. Interestingly, the QE/QSS invariant manifold of the smooth model and the reduced dynamics along this manifold can be put into correspondence to the tropical variety of the hybridization and to sliding modes along this variety, respectivelyComment: conference SASB 2011, to be published in Electronic Notes in Theoretical Computer Scienc

Approximations and their consequences for dynamic modelling of signal transduction pathways

Signal transduction is the process by which the cell converts one kind of signal or stimulus into another. This involves a sequence of biochemical reactions, carried out by proteins. The dynamic response of complex cell signalling networks can be modelled and simulated in the framework of chemical kinetics. The mathematical formulation of chemical kinetics results in a system of coupled differential equations. Simplifications can arise through assumptions and approximations. The paper provides a critical discussion of frequently employed approximations in dynamic modelling of signal transduction pathways. We discuss the requirements for conservation laws, steady state approximations, and the neglect of components. We show how these approximations simplify the mathematical treatment of biochemical networks but we also demonstrate differences between the complete system and its approximations with respect to the transient and steady state behavior

On the dynamics of the adenylate energy system: homeorhesis vs homeostasis.

Biochemical energy is the fundamental element that maintains both the adequate turnover of the biomolecular structures and the functional metabolic viability of unicellular organisms. The levels of ATP, ADP and AMP reflect roughly the energetic status of the cell, and a precise ratio relating them was proposed by Atkinson as the adenylate energy charge (AEC). Under growth-phase conditions, cells maintain the AEC within narrow physiological values, despite extremely large fluctuations in the adenine nucleotides concentration. Intensive experimental studies have shown that these AEC values are preserved in a wide variety of organisms, both eukaryotes and prokaryotes. Here, to understand some of the functional elements involved in the cellular energy status, we present a computational model conformed by some key essential parts of the adenylate energy system. Specifically, we have considered (I) the main synthesis process of ATP from ADP, (II) the main catalyzed phosphotransfer reaction for interconversion of ATP, ADP and AMP, (III) the enzymatic hydrolysis of ATP yielding ADP, and (IV) the enzymatic hydrolysis of ATP providing AMP. This leads to a dynamic metabolic model (with the form of a delayed differential system) in which the enzymatic rate equations and all the physiological kinetic parameters have been explicitly considered and experimentally tested in vitro. Our central hypothesis is that cells are characterized by changing energy dynamics (homeorhesis). The results show that the AEC presents stable transitions between steady states and periodic oscillations and, in agreement with experimental data these oscillations range within the narrow AEC window. Furthermore, the model shows sustained oscillations in the Gibbs free energy and in the total nucleotide pool. The present study provides a step forward towards the understanding of the fundamental principles and quantitative laws governing the adenylate energy system, which is a fundamental element for unveiling the dynamics of cellular life

Novel strategies for process control based on hybrid semi-parametric mathematical systems

Tese de doutoramento. Engenharia Química. Universidade do Porto. Faculdade de Engenharia. 201

Optimal Control of a Nonlinear Time-Delay System in Batch Fermentation Process

The main control goal in batch process is to get a high yield of products. In this paper, to maximize the yield of 1,3-propanediol (1,3-PD) in bioconversion of glycerol to 1,3-PD, we consider an optimal control problem involving a nonlinear time-delay system. The control variables in this problem include the initial concentrations of biomass and glycerol and the terminal time of the batch process. By a time-scaling transformation, we transcribe the optimal control problem into a new one with fixed terminal time, which yields a new nonlinear system with variable time-delay. The gradients of the cost and constraint functionals with respect to the control variables are derived using the costate method. Then, a gradient-based optimization method is developed to solve the optimal control problem. Numerical results show that the yield of 1,3-PD at the terminal time is increased considerably compared with the experimental data

Dynamic energy budget approach to evaluate antibiotic effects on biofilms

Quantifying the action of antibiotics on biofilms is essential to devise therapies against chronic infections. Biofilms are bacterial communities attached to moist surfaces, sheltered from external aggressions by a polymeric matrix. Coupling a dynamic energy budget based description of cell metabolism to surrounding concentration fields, we are able to approximate survival curves measured for different antibiotics. We reproduce numerically stratified distributions of cell types within the biofilm and introduce ways to incorporate different resistance mechanisms. Qualitative predictions follow that are in agreement with experimental observations, such as higher survival rates of cells close to the substratum when employing antibiotics targeting active cells or enhanced polymer production when antibiotics are administered. The current computational model enables validation and hypothesis testing when developing therapies.Comment: to appear in Communications in Nonlinear Science and Numerical Simulatio