Parameter optimization and sensitivity analysis for large kinetic models using a real-coded genetic algorithm

Dynamic modeling is a powerful tool for predicting changes in metabolic regulation. However, a large number of input parameters, including kinetic constants and initial metabolite concentrations, are required to construct a kinetic model. Therefore, it is important not only to optimize the kinetic parameters, but also to investigate the effects of their perturbations on the overall system. We investigated the efficiency of the use of a real-coded genetic algorithm (RCGA) for parameter optimization and sensitivity analysis in the case of a large kinetic model involving glycolysis and the pentose phosphate pathway in Escherichia coli K-12. Sensitivity analysis of the kinetic model using an RCGA demonstrated that the input parameter values had different effects on model outputs. The results showed highly influential parameters in the model and their allowable ranges for maintaining metabolite-level stability. Furthermore, it was revealed that changes in these influential parameters may complement one another. This study presents an efficient approach based on the use of an RCGA for optimizing and analyzing parameters in large kinetic models.     

Interdependence between cooperativity and control coefficients

Influence of the concentration of internal metabolites on the control coefficient (defined as fractional change in flux per fractional change in enzyme activity) and regulatory properties of a given enzyme have been studied theoretically using a cyclic model of three enzymes. This model is useful to investigate the properties of the flux control coefficient for an enzyme following different rate equations. Enzymes can have high or low values of control coefficient irrespective of the type of kinetic equation, but the results obtained show that the sensitivity of these values to substrate variations is strongly dependent on its rate equation. These results help identify which kinetic equation allows the best control of a given metabolic pathway. These results have been applied to the purine nucleotide cycle. It is demonstrated that the best control of the cycle is reached when the irreversible reaction catalyzed by AMP deaminase follows a rate law that corresponds to a rational function of 2:2 degree with respect to AMP concentration.     

Improved hidden Markov models for molecular motors, part 1: basic theory

Hidden Markov models (HMMs) provide an excellent analysis of recordings with very poor signal/noise ratio made from systems such as ion channels which switch among a few states. This method has also recently been used for modeling the kinetic rate constants of molecular motors, where the observable variable—the position—steadily accumulates as a result of the motor's reaction cycle. We present a new HMM implementation for obtaining the chemical-kinetic model of a molecular motor's reaction cycle called the variable-stepsize HMM in which the quantized position variable is represented by a large number of states of the Markov model. Unlike previous methods, the model allows for arbitrary distributions of step sizes, and allows these distributions to be estimated. The result is a robust algorithm that requires little or no user input for characterizing the stepping kinetics of molecular motors as recorded by optical techniques.     

Involvement of a Che-like signal transduction cascade in regulating cyst cell development in Rhodospirillum centenum

Homologues of the E. coli chemotaxis (Che) signal transduction pathway are present in nearly all motile bacteria. Although E. coli contains only one Che cascade, many other bacteria are known to possess multiple sets of che genes. The role of multiple che-like gene clusters could potentially code for parallel Che-like signal transduction pathways that have distinctly different input and output functions. In this study, we describe a che-like gene cluster in Rhodospirillum centenum that controls a developmental cycle. In-frame deletion mutants of homologues of CheW (DeltacheW(3a)and DeltacheW(3b)), CheR (DeltacheR(3)), CheA (DeltacheA(3)) and a methyl-accepting chemotaxis protein (Deltamcp(3)) are defective in starvation-induced formation of heat and desiccation resistant cyst cells. In contrast, mutants of homologues of CheY (DeltacheY(3)), CheB (DeltacheB(3)), and a second input kinase designated as CheS (DeltacheS(3)) result in cells that are derepressed in the formation of cysts. A model of signal transduction is presented in which there are three distinct Che-like signal transduction cascades; one that is involved in chemotaxis, one that is involved in flagella biosynthesis and the third that is involved in cyst development.     

On the stability of metabolic cycles

We investigate the stability properties of two different classes of metabolic cycles using a combination of analytical and computational methods. Using principles from structural kinetic modeling (SKM), we show that the stability of metabolic networks with certain structural regularities can be studied using a combination of analytical and computational techniques. We then apply these techniques to a class of single input, single output metabolic cycles, and find that the cycles are stable under all conditions tested. Next, we extend our analysis to a small autocatalytic cycle, and determine parameter regimes within which the cycle is very likely to be stable. We demonstrate that analytical methods can be used to understand the relationship between kinetic parameters and stability, and that results from these analytical methods can be confirmed with computational experiments. In addition, our results suggest that elevated metabolite concentrations and certain crucial saturation parameters can strongly affect the stability of the entire metabolic cycle. We discuss our results in light of the possibility that evolutionary forces may select for metabolic network topologies with a high intrinsic probability of being stable. Furthermore, our conclusions support the hypothesis that certain types of metabolic cycles may have played a role in the development of primitive metabolism despite the absence of regulatory mechanisms.     

Modelling and analysis of central metabolism operating regulatory interactions in salt stress conditions in a L-carnitine overproducing E. coli strain

Based on experimental data from E. coli cultures, we have devised a mathematical model in the GMA-power law formalism that describes the central and L-carnitine metabolism in and between two steady states, non-osmotic and hyperosmotic (0.3 M NaCl). A key feature of this model is the introduction of type of kinetic order, the osmotic stress kinetic orders (g(OSn)), derived from the power law general formalism, which represent the effect of osmotic stress in each metabolic process of the model.By considering the values of the g(OSn) linked to each metabolic process we found that osmotic stress has a positive and determinant influence on the increase in flux in energetic metabolism (glycolysis); L-carnitine biosynthesis production; the transformation/excretion of Acetyl-CoA into acetate and ethanol; the input flux of peptone into the cell; the anabolic use of pyruvate and biomass decomposition. In contrast, we found that although the osmotic stress has an inhibitory effect on the transformation flux from the glycolytic end products (pyruvate) to Acetyl-CoA, this inhibition is counteracted by other effects (the increase in pyruvate concentration) to the extent that the whole flux increases. In the same vein, the down regulation exerted by osmotic stress on fumarate uptake and its oxidation and the production and export of lactate and pyruvate are reversed by other factors up to the point that the first increased and the second remained unchanged.The model analysis shows that in osmotic conditions the energy and fermentation pathways undergo substantial rearrangement. This is illustrated by the observation that the increase in the fermentation fluxes is not connected with fluxes towards the tricaboxylic acid intermediates and the synthesis of biomass. The osmotic stress associated with these fluxes reflects these changes. All these observations support that the responses to salt stress observed in E. coli might be conserved in halophiles.Flux evolution during osmotic adaptations showed a hyperbolic (increasing or decreasing) pattern except in the case of peptone and fumarate uptake by the cell, which initially decreased. Finally, the model also throws light on the role of L-carnitine as osmoprotectant.     

Optimal parametric sensitivity control for the estimation of kinetic parameters in bioreactors

In this paper the well-known problem of optimal input design is considered. In particular, the focus is on input design for the estimation of kinetic parameters in bioreactors. The problem is formulated as follows: given the model structure (f,g), which is assumed to be affine in the input, and the specific parameter of interest theta;(k) find a feedback law that maximizes the sensitivity of the model output to the parameter under different flow conditions in the bioreactor and, possibly, minimize the input or state costs. Analytical solutions to these problems are presented. As an example a bioreactor with a biomass that grows according to the well-known Monod kinetics is considered.     

Continuous modeling of metabolic networks with gene regulation in yeast and in vivo determination of rate parameters

A continuous model of a metabolic network including gene regulation to simulate metabolic fluxes during batch cultivation of yeast Saccharomyces cerevisiae was developed. The metabolic network includes reactions of glycolysis, gluconeogenesis, glycerol and ethanol synthesis and consumption, the tricarboxylic acid cycle, and protein synthesis. Carbon sources considered were glucose and then ethanol synthesized during growth on glucose. The metabolic network has 39 fluxes, which represent the action of 50 enzymes and 64 genes and it is coupled with a gene regulation network which defines enzyme synthesis (activities) and incorporates regulation by glucose (enzyme induction and repression), modeled using ordinary differential equations. The model includes enzyme kinetics, equations that follow both mass-action law and transport as well as inducible, repressible, and constitutive enzymes of metabolism. The model was able to simulate a fermentation of S. cerevisiae during the exponential growth phase on glucose and the exponential growth phase on ethanol using only one set of kinetic parameters. All fluxes in the continuous model followed the behavior shown by the metabolic flux analysis (MFA) obtained from experimental results. The differences obtained between the fluxes given by the model and the fluxes determined by the MFA do not exceed 25% in 75% of the cases during exponential growth on glucose, and 20% in 90% of the cases during exponential growth on ethanol. Furthermore, the adjustment of the fermentation profiles of biomass, glucose, and ethanol were 95%, 95%, and 79%, respectively. With these results the simulation was considered successful. A comparison between the simulation of the continuous model and the experimental data of the diauxic yeast fermentation for glucose, biomass, and ethanol, shows an extremely good match using the parameters found. The small discrepancies between the fluxes obtained through MFA and those predicted by the differential equations, as well as the good match between the profiles of glucose, biomass, and ethanol, and our simulation, show that this simple model, that does not rely on complex kinetic expressions, is able to capture the global behavior of the experimental data. Also, the determination of parameters using a straightforward minimization technique using data at only two points in time was sufficient to produce a relatively accurate model. Thus, even with a small amount of experimental data (rates and not concentrations) it was possible to estimate the parameters minimizing a simple objective function. The method proposed allows the obtention of reasonable parameters and concentrations in a system with a much larger number of unknowns than equations. Hence a contribution of this study is to present a convenient way to find in vivo rate parameters to model metabolic and genetic networks under different conditions.     

Development of a kinetic metabolic model: application to Catharanthus roseus hairy root

A kinetic metabolic model describing Catharanthus roseus hairy root growth and nutrition was developed. The metabolic network includes glycolysis, pentose-phosphate pathway, TCA cycle and the catabolic reactions leading to cell building blocks such as amino acids, organic acids, organic phosphates, lipids and structural hexoses. The central primary metabolic network was taken at pseudo-steady state and metabolic flux analysis technique allowed reducing from 31 metabolic fluxes to 20 independent pathways. Hairy root specific growth rate was described as a function of intracellular concentration in cell building blocks. Intracellular transport and accumulation kinetics for major nutrients were included. The model uses intracellular nutrients as well as energy shuttles to describe metabolic regulation. Model calibration was performed using experimental data obtained from batch and medium exchange liquid cultures of C. roseus hairy root using a minimal medium in Petri dish. The model is efficient in estimating the growth rate.     

Model reduction of genome-scale metabolic models as a basis for targeted kinetic models

Constraint-based, genome-scale metabolic models are an essential tool to guide metabolic engineering. However, they lack the detail and time dimension that kinetic models with enzyme dynamics offer. Model reduction can be used to bridge the gap between the two methods and allow for the integration of kinetic models into the Design-Built-Test-Learn cycle. Here we show that these reduced size models can be representative of the dynamics of the original model and demonstrate the automated generation and parameterisation of such models. Using these minimal models of metabolism could allow for further exploration of dynamic responses in metabolic networks.     

Analysis of Chinese hamster ovary cell metabolism through a combined computational and experimental approach

Optimization of cell culture processes can benefit from the systematic analysis of experimental data and their organization in mathematical models, which can be used to decipher the effect of individual process variables on multiple outputs of interest. Towards this goal, a kinetic model of cytosolic glucose metabolism coupled with a population-level model of Chinese hamster ovary cells was used to analyse metabolic behavior under batch and fed-batch cell culture conditions. The model was parameterized using experimental data for cell growth dynamics, extracellular and intracellular metabolite profiles. The results highlight significant differences between the two culture conditions in terms of metabolic efficiency and motivate the exploration of lactate as a secondary carbon source. Finally, the application of global sensitivity analysis to the model parameters highlights the need for additional experimental information on cell cycle distribution to complement metabolomic analyses with a view to parameterize kinetic models.     

Sensitivity and Robustness in Covalent Modification Cycles with a Bifunctional Converter Enzyme

Regulation by covalent modification is a common mechanism to transmit signals in biological systems. The modifying reactions are catalyzed either by two distinct converter enzymes or by a single bifunctional enzyme (which may employ either one or two catalytic sites for its opposing activities). The reason for this diversification is unclear, but contemporary theoretical models predict that systems with distinct converter enzymes can exhibit enhanced sensitivity to input signals whereas bifunctional enzymes with two catalytic sites are believed to generate robustness against variations in system’s components. However, experiments indicate that bifunctional enzymes can also exhibit enhanced sensitivity due to the zero-order effect, raising the question whether both phenomena could be understood within a common mechanistic model. Here, I argue that this is, indeed, the case. Specifically, I show that bifunctional enzymes with two catalytic sites can exhibit both ultrasensitivity and concentration robustness, depending on the kinetic operating regime of the enzyme’s opposing activities. The model predictions are discussed in the context of experimental observations of ultrasensitivity and concentration robustness in the uridylylation cycle of the PII protein, and in the phosphorylation cycle of the isocitrate dehydrogenase, respectively.     

Mathematical model of trans-sarcomere exchange of calcium ions and Ca2+-dependent control of smooth muscle contractile activity

The mathematical model of smooth muscles contractile activity Ca(2+)-dependent control has been proposed on the base of Ca ions trans-sarcomal exchange biochemical mechanisms interpretation in myocytes. While analysing the model the conclusion should be made that kinetic parameters changes (in relation to Ca ions) Mg2+, ATP-dependent calcium pump of plasma membrane--Michaelis constant Km and transport process maximal velocity Vmax-render the effect on the character of the intracellular calcium transients and profile of full mechanokinetic curve. As well one more conclusion has been made that plasma membrane Mg2+, ATP-dependent calcium pump, which kinetic parameters under the physiologic conditions are subjected to modulation as the result of metabolic, pharmacologic and physico-chemical factors fulfills the essential role in supplying Ca(2+)-dependent control of the smooth muscles contractile response full cycle.     

Sensitivity and Robustness in Covalent Modification Cycles with a Bifunctional Converter Enzyme

Regulation by covalent modification is a common mechanism to transmit signals in biological systems. The modifying reactions are catalyzed either by two distinct converter enzymes or by a single bifunctional enzyme (which may employ either one or two catalytic sites for its opposing activities). The reason for this diversification is unclear, but contemporary theoretical models predict that systems with distinct converter enzymes can exhibit enhanced sensitivity to input signals whereas bifunctional enzymes with two catalytic sites are believed to generate robustness against variations in system’s components. However, experiments indicate that bifunctional enzymes can also exhibit enhanced sensitivity due to the zero-order effect, raising the question whether both phenomena could be understood within a common mechanistic model. Here, I argue that this is, indeed, the case. Specifically, I show that bifunctional enzymes with two catalytic sites can exhibit both ultrasensitivity and concentration robustness, depending on the kinetic operating regime of the enzyme’s opposing activities. The model predictions are discussed in the context of experimental observations of ultrasensitivity and concentration robustness in the uridylylation cycle of the PII protein, and in the phosphorylation cycle of the isocitrate dehydrogenase, respectively.     

Multi-random binding molecular interactions: a general kinetic model

A kinetic model is suggested to account for the interactions of several ligands with a target whose molecule possesses several independent equivalent receptor sites for each ligand (multiligand multisite model). To analyse the problem, we shall derive solutions for three elementary situations: (a) interactions of a ligand with a mono-receptor site target molecule (monosite model); (b) interactions of several ligands with a target whose molecule possesses one receptor site for each ligand involved (multiligand model); (c) interactions of a ligand with a target whose molecule possesses several receptor sites of the same kind for this ligand (multisite model). Throughout this study, every ligand molecule is assumed to offer one binding site to the target. The main implications of the corresponding analytical solutions are discussed from a molecular point of view. The results cover a great many well-known aspects of the molecular interactions in various fields such as enzymology, endocrinology, radio-immunology and saturation analysis. As suggested by the inhibition patterns obtained, this model may therefore provide a new point of view to interpret the relevant phenomena. Furthermore, a kinetic approach to the generalized mass action law can be deduced from this model, and experimental conditions in which the isotopic dilution law applies are examined.     

Transfer Functions for Protein Signal Transduction: Application to a Model of Striatal Neural Plasticity

We present a novel formulation for biochemical reaction networks in the context of protein signal transduction. The model consists of input-output transfer functions, which are derived from differential equations, using stable equilibria. We select a set of “source” species, which are interpreted as input signals. Signals are transmitted to all other species in the system (the “target” species) with a specific delay and with a specific transmission strength. The delay is computed as the maximal reaction time until a stable equilibrium for the target species is reached, in the context of all other reactions in the system. The transmission strength is the concentration change of the target species. The computed input-output transfer functions can be stored in a matrix, fitted with parameters, and even recalled to build dynamical models on the basis of state changes. By separating the temporal and the magnitudinal domain we can greatly simplify the computational model, circumventing typical problems of complex dynamical systems. The transfer function transformation of biochemical reaction systems can be applied to mass-action kinetic models of signal transduction. The paper shows that this approach yields significant novel insights while remaining a fully testable and executable dynamical model for signal transduction. In particular we can deconstruct the complex system into local transfer functions between individual species. As an example, we examine modularity and signal integration using a published model of striatal neural plasticity. The modularizations that emerge correspond to a known biological distinction between calcium-dependent and cAMP-dependent pathways. Remarkably, we found that overall interconnectedness depends on the magnitude of inputs, with higher connectivity at low input concentrations and significant modularization at moderate to high input concentrations. This general result, which directly follows from the properties of individual transfer functions, contradicts notions of ubiquitous complexity by showing input-dependent signal transmission inactivation.     

Towards Kinetic Modeling of Global Metabolic Networks: Methylobacterium extorquens AM1 Growth as Validation

Here we report a systematic method for constructing a large scale kinetic metabolic model and its initial application to the modeling of central metabolism of Methylobacterium extorquens AM1, a methylotrophic and environmental important bacterium. Its central metabolic network includes formaldehyde metabolism, serine cycle, citric acid cycle, pentose phosphate pathway, gluconeogensis, PHB synthesis and acetyl-CoA conversion pathway, respiration and energy metabolism. Through a systematic and consistent procedure of finding a set of parameters in the physiological range we overcome an outstanding difficulty in large scale kinetic modeling: the requirement for a massive number of enzymatic reaction parameters. We are able to construct the kinetic model based on general biological considerations and incomplete experimental kinetic parameters. Our method consists of the following major steps: 1) using a generic enzymatic rate equation to reduce the number of enzymatic parameters to a minimum set while still preserving their characteristics; 2) using a set of steady state fluxes and metabolite concentrations in the physiological range as the expected output steady state fluxes and metabolite concentrations for the kinetic model to restrict the parametric space of enzymatic reactions; 3) choosing enzyme constants K’s and K’eqs optimized for reactions under physiological concentrations, if their experimental values are unknown; 4) for models which do not cover the entire metabolic network of the organisms, designing a dynamical exchange for the coupling between the metabolism represented in the model and the rest not included.     

Approximative kinetic formats used in metabolic network modeling

An overview is presented of the different approximative enzyme kinetic formats that have been proposed for use in metabolic modeling studies. It is considered that the following four general properties are relevant for approximative kinetics: the rate must be proportional to enzyme level; at high metabolite concentrations, there is downward concave behavior of rate versus concentration; the number of kinetic parameters should be as small as possible; analytical solutions of steady-state network balances are desirable. Six different approximative kinetic formats are evaluated (linear, logarithmic-linear, power law GMA, power law S-systems, thermokinetic, linear-logarithmic) and it is concluded that only the recently proposed linear-logarithmic approach combines all desired properties and therefore seems a most appropriate approximate kinetic format.     

Metabolic control analysis as a mechanism that conserves genetic variance during advanced cycle breeding

The recycling of elite inbreds (i.e., advanced cycle breeding) has led to significant genetic gains but also to a narrow gene pool in plant breeding programs. Sustained yield improvements in many crops have suggested that genetic variance is not depleted at a rate predicted by an additive genetic model. Unlike the additive model in classical quantitative genetic theory, metabolic control analysis relates the variation in a biochemical process with the genetic variation in a quantitative trait. Our objective was to determine whether metabolic control analysis is a mechanism that slows the decrease in genetic variance during advanced cycle breeding. Three cycles of advanced cycle breeding were simulated with 10, 50, or 100 quantitative trait loci (QTL) controlling a trait. In metabolic control analysis, these QTL coded for enzymes involved in a linear metabolic pathway that converted a substrate into a product. In the absence of selection, both the additive model and the metabolic control analysis model led to about a 50% reduction in genetic variance from cycle to cycle. With selection, the additive model led to a 50–58% reduction in genetic variance, but the metabolic control analysis model generally led to only a 12–54% reduction. We suggest selection in a metabolic control analysis model as a mechanism that slows the decrease in genetic variance during advanced cycle breeding. This conservation of genetic variance would allow breeders to achieve genetic gains for a longer period than expected under the additive model.Communicated by H.C. Becker     

Origin and diversification of a metabolic cycle in oligomer world

Based on the oligomer-world hypothesis we propose an abstract model where the molecular recognition among oligomers is described in the shape space. The origin of life in the oligomer world is regarded as the establishment of a metabolic cycle in a primitive cell. The cycle is sustained by the molecular recognition. If an original cell acquires the ability of the replication of oligomers, the relationship among oligomers changes due to the poor fidelity of the replication. This change leads to the diversification of metabolic cycles. The selection among diverse cycles is the basis of the evolution. The evolvability is one of the essential characters of life. We demonstrate the origin and diversification of the metabolic cycle by the computer simulation of our model. Such a simulation is expected to be the simplified demonstration of what actually occurred in the primordial soup. Our model describes an analog era preceding the digital era based on the genetic code.