Search

Within the areas of Computational Organisation Theory and Artificial Intelligence, techniques have been developed to simulate and analyse dynamics within organisations in society. Usually these modelling techniques are applied to factories and to the internal organisation of their process flows, thus obtaining models of complex organisations at various levels of aggregation. The dynamics in living cells are often interpreted in terms of well-organised processes, a bacterium being considered a (micro)factory. This suggests that organisation modelling techniques may also benefit their analysis. Using the example of Escherichia coli it is shown how indeed agent-based organisational modelling techniques can be used to simulate and analyse E.coli's intracellular dynamics. Exploiting the abstraction levels entailed by this perspective, a concise model is obtained that is readily simulated and analysed at the various levels of aggregation, yet shows the cell's essential dynamic patterns. © 2010 Springer-Verlag Berlin Heidelberg.

What is the highest possible replication rate for living organisms? The cellular growth rate is controlled by a variety of processes. Therefore, it is unclear which metabolic process or group of processes should be activated to increase growth rate. An organism that is already growing fast may already have optimized through evolution all processes that could be optimized readily, but may be confronted with a more generic limitation. Here we introduce a method called 'cell selection' to select for highest growth rate, and show how such a cellular site of 'growth control' was identified. By applying pH-auxostat cultivation to the already fast-growing yeast Kluyveromyces marxianus for a sufficiently long time, we selected a strain with a 30% increased growth rate; its cell-cycle time decreased to 52 min, much below that reported to date for any eukaryote. The increase in growth rate was accompanied by a 40% increase in cell surface at a fairly constant cell volume. We show how the increase in growth rate can be explained by a dominant (80%) limitation of growth by the group of membrane processes (a 0.7% increase of specific growth rate to a 1% increase in membrane surface area). Simultaneous activation of membrane processes may be what is required to accelerate growth of the fastest-growing form of eukaryotic life to growth rates that are even faster, and may be of potential interest for single-cell protein production in industrial 'White' biotechnology processes. © 2008 The Authors.

Genomic data allow the large-scale manual or semi-automated assembly of metabolic network reconstructions, which provide highly curated organism-specific knowledge bases. Although several genome-scale network reconstructions describe Saccharomyces cerevisiae metabolism, they differ in scope and content, and use different terminologies to describe the same chemical entities. This makes comparisons between them difficult and underscores the desirability of a consolidated metabolic network that collects and formalizes the 'community knowledge' of yeast metabolism. We describe how we have produced a consensus metabolic network reconstruction for S. cerevisiae. In drafting it, we placed special emphasis on referencing molecules to persistent databases or using database-independent forms, such as SMILES or InChI strings, as this permits their chemical structure to be represented unambiguously and in a manner that permits automated reasoning. The reconstruction is readily available via a publicly accessible database and in the Systems Biology Markup Language (http://www.comp-sys-bio.org/yeastnet). It can be maintained as a resource that serves as a common denominator for studying the systems biology of yeast. Similar strategies should benefit communities studying genome-scale metabolic networks of other organisms. © 2008 Nature Publishing Group.

In both industrial fermenters and in their natural habitats, microorganisms often experience an inhomogeneous and fluctuating environment. In this paper we mimicked one aspect of this nonideal behaviour by imposing a low and oscillating extracellular glucose concentration on nonoscillating suspensions of yeast cells. The extracellular dynamics changed the intracellular dynamics - which was monitored through NADH fluorescence - from steady to equally dynamic; the latter followed the extracellular dynamics at the frequency of glucose pulsing. Interestingly, the amplitude of the oscillation of the NADH fluorescence increased with time. This increase in amplitude was sensitive to inhibition of protein synthesis, and was due to a change in the cells rather than in the medium; the cell population was 'trained' to respond to the extracellular dynamics. To examine the mechanism behind this 'training', we subjected the cells to a low and constant extracellular glucose concentration. Seventy-five minutes of adaptation to a low and constant glucose concentration induced the same increase of the amplitude of the forced NADH oscillations as did the train of glucose pulses. Furthermore, 75 min of adaptation to a low (oscillating or continuous) glucose concentration decreased the K

The effects of pre-incubation with mercury (Hg

The effects of pre-incubation with mercury (Hg2+) and cadmium (Cd2+) on the activities of individual glycolytic enzymes, on the flux and on internal metabolite concentrations of the upper part of glycolysis were investigated in mouse muscle extracts. In the range of metal concentrations analysed we found that only hexokinase and phosphofructokinase, the enzymes that shared the control of the flux, were inhibited by Hg2+ and Cd2+. The concentrations of the internal metabolites glucose-6-phosphate and fructose-6-phosphate did not change significantly when Hg2+ and Cd2+ were added. A mathematical model was constructed to explore the mechanisms of inhibition of Hg2+ and Cd2+ on hexokinase and phosphofructokinase. Equations derived from detailed mechanistic models for each inhibition were fitted to the experimental data. In a concentration-dependent manner these equations describe the observed inhibition of enzyme activity. Under the conditions analysed, the integral model showed that the simultaneous inhibition of hexokinase and phosphofructokinase explains the observation that the concentrations of glucose-6-phosphate and fructose-6-phosphate did not change as the heavy metals decreased the glycolytic flux.

Chronic myeloid leukaemia (CML) is a clonal myeloproliferative disease resulting from an aberrant BCR.ABL gene and protein. To predict BCR.ABL protein abundance and phosphorylation in individual cells in a population of CML cells, we modelled BCR.ABL protein regulation through associated miRNAs using a systems approach. The model rationalizes the level of BCR.ABL protein heterogeneity in CML cells in correlation with the heterogeneous BCR.ABL mRNA levels. We also measured BCR.ABL mRNA and BCR.ABLp phosphorylation in individual cells. The experimental data were consistent with the modelling results, thereby partly validating the model. Provided it is tested further, the model may be used to support effective therapeutic strategies including the combined application of a tyrosine kinase inhibitor and miRNAs targeting BCR.ABL. It appears able to predict different effects of the two types of drug on cells with different expression levels and consequently different effects on the generation of resistance.

Regulation analysis is a methodology that quantifies to what extent a change in the flux through a metabolic pathway is regulated by either gene expression or metabolism. Two extensions to regulation analysis were developed over the past years: (i) the regulation of Vmax can be dissected into the various levels of the gene-expression cascade, such as transcription, translation, protein degradation, etc. and (ii) a time-dependent version allows following flux regulation when cells adapt to changes in their environment. The methodology of the original form of regulation analysis as well as of the two extensions will be described in detail. In addition, we will show what is needed to apply regulation analysis in practice. Studies in which the different versions of regulation analysis were applied revealed that flux regulation was distributed over various processes and depended on time, enzyme, and condition of interest. In the case of the regulation of glycolysis in baker's yeast, it appeared, however, that cells that remain under respirofermentative conditions during a physiological challenge tend to invoke more gene-expression regulation, while a shift between respirofermentative and respiratory conditions invokes an important contribution of metabolic regulation. The complexity of the regulation observed in these studies raises the question what is the advantage of this highly distributed and condition-dependent flux regulation. © 2011 Elsevier Inc. All rights reserved.

Awareness is growing that drug target validation should involve systems analysis of cellular networks. There is less appreciation, though, that the composition of networks may change in response to drugs. If the response is homeostatic (e.g. through upregulation of the target protein), this may neutralize the inhibitory effect. In this scenario the effect on cell growth and survival would be less than anticipated based on affinity of the drug for its target. Glycolysis is the sole free-energy source for the deadly parasite Trypanosoma brucei and is therefore a possible target pathway for anti-trypanosomal drugs. Plasma-membrane glucose transport exerts high control over trypanosome glycolysis and hence the transporter is a promising drug target. Here we show that at high inhibitor concentrations, inhibition of trypanosome glucose transport causes cell death. Most interestingly, sublethal concentrations initiate a domino effect in which network adaptations enhance inhibition. This happens via (i) metabolic control exerted by the target protein, (ii) decreases in mRNAs encoding the target protein and other proteins in the same pathway, and (iii) partial differentiation of the cells leading to (low) expression of immunogenic insect-stage coat proteins. We discuss how these 'anti-homeostatic' responses together may facilitate killing of parasites at an acceptable drug dosage. © 2010 Blackwell Publishing Ltd.

Realistic quantitative models require data from many laboratories. Therefore, standardization of experimental systems and assay conditions is crucial. Moreover, standards should be representative of the in vivo conditions. However, most often, enzyme-kinetic parameters are measured under assay conditions that yield the maximum activity of each enzyme. In practice, this means that the kinetic parameters of different enzymes are measured in different buffers, at different pH values, with different ionic strengths, etc. In a joint effort of the Dutch Vertical Genomics Consortium, the European Yeast Systems Biology Network and the Standards for Reporting Enzymology Data Commission, we have developed a single assay medium for determining enzyme-kinetic parameters in yeast. The medium is as close as possible to the in vivo situation for the yeast Saccharomyces cerevisiae, and at the same time is experimentally feasible. The in vivo conditions were estimated for S. cerevisiae strain CEN. PK113-7D grown in aerobic glucose-limited chemostat cultures at an extracellular pH of 5.0 and a specific growth rate of 0.1 h(-1). The cytosolic pH and concentrations of calcium, sodium, potassium, phosphorus, sulfur and magnesium were determined. On the basis of these data and literature data, we propose a defined in vivo-like medium containing 300 mm potassium, 50 mm phosphate, 245 mm glutamate, 20 mm sodium, 2 mm free magnesium and 0.5 mm calcium, at a pH of 6.8. The V-max values of the glycolytic and fermentative enzymes of S. cerevisiae were measured in the new medium. For some enzymes, the results deviated conspicuously from those of assays done under enzyme-specific, optimal conditions.

Numerous top-down kinetic models have been constructed to describe the cell cycle. These models have typically been constructed, validated and analyzed using model species (molecular intermediates and proteins) and phenotypic observations, and therefore do not focus on the individual model processes (reaction steps). We have developed a method to: (a) quantify the importance of each of the reaction steps in a kinetic model for the positioning of a switch point [i.e. the restriction point (RP)]; (b) relate this control of reaction steps to their effects on molecular species, using sensitivity and co-control analysis; and thereby (c) go beyond a correlation towards a causal relationship between molecular species and effects. The method is generic and can be applied to responses of any type, but is most useful for the analysis of dynamic and emergent responses such as switch points in the cell cycle. The strength of the analysis is illustrated for an existing mammalian cell cycle model focusing on the RP [Novak B, Tyson J (2004) J Theor Biol230, 563-579]. The reactions in the model with the highest RP control were those involved in: (a) the interplay between retinoblastoma protein and E2F transcription factor; (b) those synthesizing the delayed response genes and cyclin D/Cdk4 in response to growth signals; (c) the E2F-dependent cyclin E/Cdk2 synthesis reaction; as well as (d) p27 formation reactions. Nine of the 23 intermediates were shown to have a good correlation between their concentration control and RP control. Sensitivity and co-control analysis indicated that the strongest control of the RP is mediated via the cyclin E/Cdk2:p27 complex concentration. Any perturbation of the RP could be related to a change in the concentration of this complex; apparent effects of other molecular species were indirect and always worked through cyclin E/Cdk2:p27, indicating a causal relationship between this complex and the positioning of the RP. © 2009 FEBS.

Much of enzyme kinetics builds on simplifications enabled by the quasi-steady-state approximation and is highly useful when the concentration of the enzyme is much lower than that of its substrate. However, in vivo, this condition is often violated. In the present study, we show that, under conditions of realistic yet high enzyme concentrations, the quasi-steady-state approximation may readily be off by more than a factor of four when predicting concentrations. We then present a novel extension of the quasi-steady-state approximation based on the zero-derivative principle, which requires considerably less theoretical work than did previous such extensions. We show that the first-order zero-derivative principle, already describes much more accurately the true enzyme dynamics at enzyme concentrations close to the concentration of their substrates. This should be particularly relevant for enzyme kinetics where the substrate is an enzyme, such as in phosphorelay and mitogen-activated protein kinase pathways. We illustrate this for the important example of the phosphotransferase system involved in glucose uptake, metabolism and signaling. We find that this system, with a potential complexity of nine dimensions, can be understood accurately using the first-order zero-derivative principle in terms of the behavior of a single variable with all other concentrations constrained to follow that behavior. © 2009 FEBS.

Time-dependent regulation analysis is a new methodology that allows us to unravel, both quantitatively and dynamically, how and when functional changes in the cell are brought about by the interplay of gene expression and metabolism. In this first experimental implementation, we dissect the initial and late response of baker's yeast upon a switch from glucose-limited growth to nitrogen starvation. During nitrogen starvation, unspecific bulk degradation of cytosolic proteins and small organelles (autophagy) occurs. If this is the primary cause of loss of glycolytic capacity, one would expect the cells to regulate their glycolytic capacity through decreasing simultaneously and proportionally the capacities of the enzymes in the first hour of nitrogen starvation. This should lead to regulation of the flux which is initially dominated by changes in the enzyme capacity. However, metabolic regulation is also known to act fast. To analyse the interplay between autophagy and metabolism, we examined the first 4 h of nitrogen starvation in detail using time-dependent regulation analysis. Some enzymes were initially regulated more by a breakdown of enzyme capacity and only later through metabolic regulation. However, other enzymes were regulated metabolically in the first hours and then shifted towards regulation via enzyme capacity. We conclude that even initial regulation is subtle and governed by different molecular levels. © 2009 FEBS.

In Saccharomyces cerevisiae the HXK2 gene, which encodes the glycolytic enzyme hexokinase II, is involved in the regulatory mechanism known as ‘glucose repression’. Its deletion leads to fully respiratory growth at high glucose concentrations where the wild type ferments profusely. Here we describe that deletion of the HXK2 gene resulted in a 75% reduction in fermentative capacity. Using regulation analysis we found that the fluxes through most glycolytic and fermentative enzymes were regulated cooperatively by changes in their capacities (Vmax) and by changes in the way they interacted with the rest of the metabolism. Glucose transport and phosphofructokinase were regulated purely at the metabolic level. The reduction of fermentative capacity in the mutant was accompanied by a remarkable resilience of the remaining capacity to nutrient starvation. After starvation, the fermentative capacity of the hxk2Δ mutant was similar to that of the wild type. Based on our results and previous reports, we suggest an inverse correlation between glucose repression and the resilience of fermentative capacity towards nutrient starvation. Only a limited number of glycolytic enzyme activities changed upon starvation of the hxk2Δ mutant and we discuss to what extent this could explain the stability of the fermentative capacity.