Cooperative response and clustering: Consequences of membrane-mediated interactions among mechanosensitive channels.

Mechanosensitive channels are ion channels which act as cells' safety valves, opening when the osmotic pressure becomes too high and making cells avoid damage by releasing ions. They are found on the cellular membrane of a large number of organisms. They interact with each other by means of deformations they induce in the membrane. We show that collective dynamics arising from the interchannel interactions lead to first- and second-order phase transitions in the fraction of open channels in equilibrium relating to the formation of channel clusters. We show that this results in a considerable delay of the response of cells to osmotic shocks, and to an extreme cell-to-cell stochastic variations in their response times, despite the large numbers of channels present in each cell. We discuss how our results are relevant for E. coli.

Oxygen-Driven Tumour Growth Model: A Pathology-Relevant Mathematical Approach.

Xenografts--as simplified animal models of cancer-differ substantially in vasculature and stromal architecture when compared to clinical tumours. This makes mathematical model-based predictions of clinical outcome challenging. Our objective is to further understand differences in tumour progression and physiology between animal models and the clinic. To achieve that, we propose a mathematical model based upon tumour pathophysiology, where oxygen--as a surrogate for endocrine delivery--is our main focus. The Oxygen-Driven Model (ODM), using oxygen diffusion equations, describes tumour growth, hypoxia and necrosis. The ODM describes two key physiological parameters. Apparent oxygen uptake rate (k'R) represents the amount of oxygen cells seem to need to proliferate. The more oxygen they appear to need, the more the oxygen transport. k'R gathers variability from the vasculature, stroma and tumour morphology. Proliferating rate (kp) deals with cell line specific factors to promote growth. The KH,KN describe the switch of hypoxia and necrosis. Retrospectively, using archived data, we looked at longitudinal tumour volume datasets for 38 xenografted cell lines and 5 patient-derived xenograft-like models. Exploration of the parameter space allows us to distinguish 2 groups of parameters. Group 1 of cell lines shows a spread in values of k'R and lower kp, indicating that tumours are poorly perfused and slow growing. Group 2 share the value of the oxygen uptake rate (k'R) and vary greatly in kp, which we interpret as having similar oxygen transport, but more tumour intrinsic variability in growth. However, the ODM has some limitations when tested in explant-like animal models, whose complex tumour-stromal morphology may not be captured in the current version of the model. Incorporation of stroma in the ODM will help explain these discrepancies. We have provided an example. The ODM is a very simple -and versatile- model suitable for the design of preclinical experiments, which can be modified and enhanced whilst maintaining confidence in its predictions.

Integrative Model of Oxidative Stress Adaptation in the Fungal Pathogen Candida albicans.

The major fungal pathogen of humans, Candida albicans, mounts robust responses to oxidative stress that are critical for its virulence. These responses counteract the reactive oxygen species (ROS) that are generated by host immune cells in an attempt to kill the invading fungus. Knowledge of the dynamical processes that instigate C. albicans oxidative stress responses is required for a proper understanding of fungus-host interactions. Therefore, we have adopted an interdisciplinary approach to explore the dynamical responses of C. albicans to hydrogen peroxide (H2O2). Our deterministic mathematical model integrates two major oxidative stress signalling pathways (Cap1 and Hog1 pathways) with the three major antioxidant systems (catalase, glutathione and thioredoxin systems) and the pentose phosphate pathway, which provides reducing equivalents required for oxidative stress adaptation. The model encapsulates existing knowledge of these systems with new genomic, proteomic, transcriptomic, molecular and cellular datasets. Our integrative approach predicts the existence of alternative states for the key regulators Cap1 and Hog1, thereby suggesting novel regulatory behaviours during oxidative stress. The model reproduces both existing and new experimental observations under a variety of scenarios. Time- and dose-dependent predictions of the oxidative stress responses for both wild type and mutant cells have highlighted the different temporal contributions of the various antioxidant systems during oxidative stress adaptation, indicating that catalase plays a critical role immediately following stress imposition. This is the first model to encapsulate the dynamics of the transcriptional response alongside the redox kinetics of the major antioxidant systems during H2O2 stress in C. albicans.

The dynamics of genome replication using deep sequencing.

Eukaryotic genomes are replicated from multiple DNA replication origins. We present complementary deep sequencing approaches to measure origin location and activity in Saccharomyces cerevisiae. Measuring the increase in DNA copy number during a synchronous S-phase allowed the precise determination of genome replication. To map origin locations, replication forks were stalled close to their initiation sites; therefore, copy number enrichment was limited to origins. Replication timing profiles were generated from asynchronous cultures using fluorescence-activated cell sorting. Applying this technique we show that the replication profiles of haploid and diploid cells are indistinguishable, indicating that both cell types use the same cohort of origins with the same activities. Finally, increasing sequencing depth allowed the direct measure of replication dynamics from an exponentially growing culture. This is the first time this approach, called marker frequency analysis, has been successfully applied to a eukaryote. These data provide a high-resolution resource and methodological framework for studying genome biology.

High-resolution replication profiles define the stochastic nature of genome replication initiation and termination.

Eukaryotic genome replication is stochastic, and each cell uses a different cohort of replication origins. We demonstrate that interpreting high-resolution Saccharomyces cerevisiae genome replication data with a mathematical model allows quantification of the stochastic nature of genome replication, including the efficiency of each origin and the distribution of termination events. Single-cell measurements support the inferred values for stochastic origin activation time. A strain, in which three origins were inactivated, confirmed that the distribution of termination events is primarily dictated by the stochastic activation time of origins. Cell-to-cell variability in origin activity ensures that termination events are widely distributed across virtually the whole genome. We propose that the heterogeneity in origin usage contributes to genome stability by limiting potentially deleterious events from accumulating at particular loci.

Stochastic association of neighboring replicons creates replication factories in budding yeast.

Inside the nucleus, DNA replication is organized at discrete sites called replication factories, consisting of DNA polymerases and other replication proteins. Replication factories play important roles in coordinating replication and in responding to replication stress. However, it remains unknown how replicons are organized for processing at each replication factory. Here we address this question using budding yeast. We analyze how individual replicons dynamically organized a replication factory using live-cell imaging and investigate how replication factories were structured using super-resolution microscopy. Surprisingly, we show that the grouping of replicons within factories is highly variable from cell to cell. Once associated, however, replicons stay together relatively stably to maintain replication factories. We derive a coherent genome-wide mathematical model showing how neighboring replicons became associated stochastically to form replication factories, which was validated by independent microscopy-based analyses. This study not only reveals the fundamental principles promoting replication factory organization in budding yeast, but also provides insight into general mechanisms by which chromosomes organize sub-nuclear structures.

Integrated stress response of Escherichia coli to methylglyoxal: transcriptional readthrough from the nemRA operon enhances protection through increased expression of glyoxalase I.

Methylglyoxal (MG) elicits activation of K(+) efflux systems to protect cells against the toxicity of the electrophile. ChIP-chip targeting RNA polymerase, supported by a range of other biochemical measurements and mutant creation, was used to identify genes transcribed in response to MG and which complement this rapid response. The SOS DNA repair regulon is induced at cytotoxic levels of MG, even when exposure to MG is transient. Glyoxalase I alone among the core MG protective systems is induced in response to MG exposure. Increased expression is an indirect consequence of induction of the upstream nemRA operon, encoding an enzyme system that itself does not contribute to MG detoxification. Moreover, this induction, via nemRA only occurs when cells are exposed to growth inhibitory concentrations of MG. We show that the kdpFABCDE genes are induced and that this expression occurs as a result of depletion of cytoplasmic K(+) consequent upon activation of the KefGB K(+) efflux system. Finally, our analysis suggests that the transcriptional changes in response to MG are a culmination of the damage to DNA and proteins, but that some integrate specific functions, such as DNA repair, to augment the allosteric activation of the main protective system, KefGB.

A systems biology analysis of long and short-term memories of osmotic stress adaptation in fungi.

BACKGROUND: Saccharomyces cerevisiae senses hyperosmotic conditions via the HOG signaling network that activates the stress-activated protein kinase, Hog1, and modulates metabolic fluxes and gene expression to generate appropriate adaptive responses. The integral control mechanism by which Hog1 modulates glycerol production remains uncharacterized. An additional Hog1-independent mechanism retains intracellular glycerol for adaptation. Candida albicans also adapts to hyperosmolarity via a HOG signaling network. However, it remains unknown whether Hog1 exerts integral or proportional control over glycerol production in C. albicans. RESULTS: We combined modeling and experimental approaches to study osmotic stress responses in S. cerevisiae and C. albicans. We propose a simple ordinary differential equation (ODE) model that highlights the integral control that Hog1 exerts over glycerol biosynthesis in these species. If integral control arises from a separation of time scales (i.e. rapid HOG activation of glycerol production capacity which decays slowly under hyperosmotic conditions), then the model predicts that glycerol production rates elevate upon adaptation to a first stress and this makes the cell adapts faster to a second hyperosmotic stress. It appears as if the cell is able to remember the stress history that is longer than the timescale of signal transduction. This is termed the long-term stress memory. Our experimental data verify this. Like S. cerevisiae, C. albicans mimimizes glycerol efflux during adaptation to hyperosmolarity. Also, transient activation of intermediate kinases in the HOG pathway results in a short-term memory in the signaling pathway. This determines the amplitude of Hog1 phosphorylation under a periodic sequence of stress and non-stressed intervals. Our model suggests that the long-term memory also affects the way a cell responds to periodic stress conditions. Hence, during osmohomeostasis, short-term memory is dependent upon long-term memory. This is relevant in the context of fungal responses to dynamic and changing environments. CONCLUSIONS: Our experiments and modeling have provided an example of identifying integral control that arises from time-scale separation in different processes, which is an important functional module in various contexts.

Combinatorial stresses kill pathogenic Candida species.

Pathogenic microbes exist in dynamic niches and have evolved robust adaptive responses to promote survival in their hosts. The major fungal pathogens of humans, Candida albicans and Candida glabrata, are exposed to a range of environmental stresses in their hosts including osmotic, oxidative and nitrosative stresses. Significant efforts have been devoted to the characterization of the adaptive responses to each of these stresses. In the wild, cells are frequently exposed simultaneously to combinations of these stresses and yet the effects of such combinatorial stresses have not been explored. We have developed a common experimental platform to facilitate the comparison of combinatorial stress responses in C. glabrata and C. albicans. This platform is based on the growth of cells in buffered rich medium at 30°C, and was used to define relatively low, medium and high doses of osmotic (NaCl), oxidative (H(2)O(2)) and nitrosative stresses (e.g., dipropylenetriamine (DPTA)-NONOate). The effects of combinatorial stresses were compared with the corresponding individual stresses under these growth conditions. We show for the first time that certain combinations of combinatorial stress are especially potent in terms of their ability to kill C. albicans and C. glabrata and/or inhibit their growth. This was the case for combinations of osmotic plus oxidative stress and for oxidative plus nitrosative stress. We predict that combinatorial stresses may be highly significant in host defences against these pathogenic yeasts.

Optimal placement of origins for DNA replication.

DNA replication is an essential process in biology and its timing must be robust so that cells can divide properly. Random fluctuations in the formation of replication starting points, called origins, and the subsequent activation of proteins lead to variations in the replication time. We analyze these stochastic properties of DNA and derive the positions of origins corresponding to the minimum replication time. We show that under some conditions the minimization of replication time leads to the grouping of origins, and relate this to experimental data in a number of species showing origin grouping.

Are the fractal skeletons the explanation for the narrowing of arteries due to cell trapping in a disturbed blood flow?

We show that common circulatory diseases, such as stenoses and aneurysms, generate chaotic advection of blood particles. This phenomenon has major consequences on the way the biochemical particles behave. Chaotic advection leads to a peculiar filamentary particle distribution, which in turn creates a favorable environment for particle reactions. Furthermore, we argue that the enhanced stretching dynamics induced by chaos can lead to the activation of platelets, particles involved in the thrombus formation. In particular, we vary the size of both stenoses and aneurysms, and model them under resting and exercising conditions. We show that the filamentary particle distribution, governed by the fractal skeleton, depends on the size of the vessel wall irregularity, and investigate how it varies under resting or exercising conditions.

A matter of life or death: modeling DNA damage and repair in bacteria.

DNA damage is a hazard all cells must face, and evolution has created a number of mechanisms to repair damaged bases in the chromosome. Paradoxically, many of these repair mechanisms can create double-strand breaks in the DNA molecule which are fatal to the cell. This indicates that the connection between DNA repair and death is far from straightforward, and suggests that the repair mechanisms can be a double-edged sword. In this report, we formulate a mathematical model of the dynamics of DNA damage and repair, and we obtain analytical expressions for the death rate. We predict a counterintuitive relationship between survival and repair. We can discriminate between two phases: below a critical threshold in the number of repair enzymes, the half-life decreases with the number of repair enzymes, but becomes independent of the number of repair enzymes above the threshold. We are able to predict quantitatively the dependence of the death rate on the damage rate and other relevant parameters. We verify our analytical results by simulating the stochastic dynamics of DNA damage and repair. Finally, we also perform an experiment with Escherichia coli cells to test one of the predictions of our model.

Reacting particles in open chaotic flows.

We study the collision probability p of particles advected by open flows with chaotic advection. We show that p scales with the particle size (or, alternatively, reaction distance) δ as a power law whose coefficient is determined by the fractal dimensions of the invariant sets defined by the advection dynamics. We also argue that this same scaling also holds for the reaction rate of active particles in the low-density regime. These analytical results are compared to numerical simulations, and they are found to agree very well.

Fractal structures in stenoses and aneurysms in blood vessels.

Recent advances in the field of chaotic advection provide the impetus to revisit the dynamics of particles transported by blood flow in the presence of vessel wall irregularities. The irregularity, being either a narrowing or expansion of the vessel, mimicking stenoses or aneurysms, generates abnormal flow patterns that lead to a peculiar filamentary distribution of advected particles, which, in the blood, would include platelets. Using a simple model, we show how the filamentary distribution depends on the size of the vessel wall irregularity, and how it varies under resting or exercise conditions. The particles transported by blood flow that spend a long time around a disturbance either stick to the vessel wall or reside on fractal filaments. We show that the faster flow associated with exercise creates widespread filaments where particles can get trapped for a longer time, thus allowing for the possible activation of such particles. We argue, based on previous results in the field of active processes in flows, that the non-trivial long-time distribution of transported particles has the potential to have major effects on biochemical processes occurring in blood flow, including the activation and deposition of platelets. One aspect of the generality of our approach is that it also applies to other relevant biological processes, an example being the coexistence of plankton species investigated previously.

Random fluctuation leads to forbidden escape of particles.

A great number of physical processes are described within the context of Hamiltonian scattering. Previous studies have rather been focused on trajectories starting outside invariant structures, since the ones starting inside are expected to stay trapped there forever. This is true though only for the deterministic case. We show however that, under finitely small random fluctuations of the field, trajectories starting inside Kolmogorov-Arnold-Moser (KAM) islands escape within finite time. The nonhyperbolic dynamics gains then hyperbolic characteristics due to the effect of the random perturbed field. As a consequence, trajectories which are started inside KAM curves escape with hyperboliclike time decay distribution, and the fractal dimension of a set of particles that remain in the scattering region approaches that for hyperbolic systems. We show a universal quadratic power law relating the exponential decay to the amplitude of noise. We present a random walk model to relate this distribution to the amplitude of noise, and investigate these phenomena with a numerical study applying random maps.

Mathematical modelling of whole chromosome replication.

All chromosomes must be completely replicated prior to cell division, a requirement that demands the activation of a sufficient number of appropriately distributed DNA replication origins. Here we investigate how the activity of multiple origins on each chromosome is coordinated to ensure successful replication. We present a stochastic model for whole chromosome replication where the dynamics are based upon the parameters of individual origins. Using this model we demonstrate that mean replication time at any given chromosome position is determined collectively by the parameters of all origins. Combining parameter estimation with extensive simulations we show that there is a range of model parameters consistent with mean replication data, emphasising the need for caution in interpreting such data. In contrast, the replicated-fraction at time points through S phase contains more information than mean replication time data and allowed us to use our model to uniquely estimate many origin parameters. These estimated parameters enable us to make a number of predictions that showed agreement with independent experimental data, confirming that our model has predictive power. In summary, we demonstrate that a stochastic model can recapitulate experimental observations, including those that might be interpreted as deterministic such as ordered origin activation times.

Escape from attracting sets in randomly perturbed systems.

The dynamics of escape from an attractive state due to random perturbations is of central interest to many areas in science. Previous studies of escape in chaotic systems have rather focused on the case of unbounded noise, usually assumed to have Gaussian distribution. In this paper, we address the problem of escape induced by bounded noise. We show that the dynamics of escape from an attractor's basin is equivalent to that of a closed system with an appropriately chosen "hole." Using this equivalence, we show that there is a minimum noise amplitude above which escape takes place, and we derive analytical expressions for the scaling of the escape rate with noise amplitude near the escape transition. We verify our analytical predictions through numerical simulations of two well-known two-dimensional maps with noise.

Emerging attractors and the transition from dissipative to conservative dynamics.

The topological structure of basin boundaries plays a fundamental role in the sensitivity to the final state in chaotic dynamical systems. Herewith we present a study on the dynamics of dissipative systems close to the Hamiltonian limit, emphasizing the increasing number of periodic attractors, and on the structural changes in their basin boundaries as the dissipation approaches zero. We show numerically that a power law with nontrivial exponent describes the growth of the total number of periodic attractors as the damping is decreased. We also establish that for small scales the dynamics is governed by effective dynamical invariants, whose measure depends not only on the region of the phase space but also on the scale under consideration. Therefore, our results show that the concept of effective invariants is also relevant for dissipative systems.

Transition to chaotic scattering: signatures in the differential cross section.

We show that bifurcations in chaotic scattering manifest themselves through the appearance of an infinitely fine-scale structure of singularities in the cross section. These "rainbow singularities" are created in a cascade, which is closely related to the bifurcation cascade undergone by the set of trapped orbits (the chaotic saddle). This cascade provides a signature in the differential cross section of the complex pattern of bifurcations of orbits underlying the transition to chaotic scattering. We show that there is a power law with a universal coefficient governing the sequence of births of rainbow singularities and we verify this prediction by numerical simulations.

Onset of chaotic advection in open flows.

In this paper we investigate the transition to chaos in the motion of particles advected by open flows with obstacles. By means of a topological argument, we show that the separation points on the surface of the obstacle imply the existence of a saddle point downstream from the obstacle, with an associated heteroclinic orbit. We argue that as soon as the flow becomes time periodic, these orbits give rise to heteroclinic tangles, causing passively advected particles to experience transient chaos. The transition to chaos thus coincides with the onset of time dependence in open flows with stagnant points, in contrast with flows with no stagnant points. We also show that the nonhyperbolic nature of the dynamics near the walls causes anomalous scalings in the vicinity of the transition. These results are confirmed by numerical simulations of the two-dimensional flow around a cylinder.