ODE, RDE and SDE models of cell cycle dynamics and clustering in yeast

A Nonlinear Delay Model for Metabolic Oscillations in Yeast Cells

We introduce two time-delay models of metabolic oscillations in yeast cells. Our model tests a hypothesis that the oscillations occur as multiple pathways share a limited resource which we equate to the number of available ribosomes. We initially explore a single-protein model with a constraint equation governing the total resource available to the cell. The model is then extended to include three proteins that share a resource pool. Three approaches are considered at constant delay to numerically detect oscillations. First, we use a spectral element method to approximate the system as a discrete map and evaluate the stability of the linearized system about its equilibria by examining its eigenvalues. For the second method, we plot amplitudes of the simulation trajectories in 2D projections of the parameter space. We use a history function that is consistent with published experimental results to obtain metabolic oscillations. Finally, the spectral element method is used to convert the system to a boundary value problem whose solutions correspond to approximate periodic solutions of the system. Our results show that certain combinations of total resource available and the time delay, lead to oscillations. We observe that an oscillation region in the parameter space is between regions admitting steady states that correspond to zero and constant production. Similar behavior is found with the three-protein model where all proteins require the same production time. However, a shift in the protein production rates peaks occurs for low available resource suggesting that our model captures the shared resource pool dynamics.

Universality of stable multi-cluster periodic solutions in a population model of the cell cycle with negative feedback

We study a population model where cells in one part of the cell cycle may affect the progress of cells in another part. If the influence, or feedback, from one part to another is negative, simulations of the model almost always result in multiple temporal clusters formed by groups of cells. We study regions in parameter space where periodic 'k-cyclic' solutions are stable. The regions of stability coincide with sub-triangles on which certain events occur in a fixed order. For boundary sub-triangles with order 'rs1', we prove that the k-cyclic periodic solution is asymptotically stable if the index of the sub-triangle is relatively prime with respect to the number of clusters k and neutrally stable otherwise. For negative linear feedback, we prove that the interior of the parameter set is covered by stable sub-triangles, i.e. a stable k-cyclic solution always exists for some k. We observe numerically that the result also holds for many forms of nonlinear feedback, but may break down in extreme cases.

Instability of the steady state solution in cell cycle population structure models with feedback

We show that when cell-cell feedback is added to a model of the cell cycle for a large population of cells, then instability of the steady state solution occurs in many cases. We show this in the context of a generic agent-based ODE model. If the feedback is positive, then instability of the steady state solution is proved for all parameter values except for a small set on the boundary of parameter space. For negative feedback we prove instability for half the parameter space. We also show by example that instability in the other half may be proved on a case by case basis.

Evidence for internuclear signaling in drosophila embryogenesis

BACKGROUND

Syncytial nuclei in Drosophila embryos undergo their first 13 divisions nearly synchronously. In the last several cell cycles, these division events travel across the anterior-posterior axis of the syncytial blastoderm in a wave. The phenomenon is well documented but the underlying mechanisms are not yet understood.

RESULTS

We study timing and positional data obtained from in vivo imaging of Drosophila embryos. We determine the statistical properties of the distribution of division times within and across generations with the null hypothesis that timing of division events is an independent random variable for each nucleus. We also compare timing data with a model of Drosophila cell cycle regulation that does not include internuclear signaling, and to a universal model of phase-dependent signaling to determine the probable form of internuclear signaling in the syncytial embryo.

CONCLUSIONS

The statistical variance of division times is lower than one would expect from uncoordinated activity. In fact, the variance decreases between the 10th and 11th divisions, which demonstrates a contribution of internuclear signaling to the observed synchrony and division waves. Our comparison with a coupled oscillator model leads us to conclude that internuclear signaling must be of Response/Signaling type with a positive impulse.

Noise-induced dispersion and breakup of clusters in cell cycle dynamics

We study the effects of random perturbations on collective dynamics of a large ensemble of interacting cells in a model of the cell division cycle. We consider a parameter region for which the unperturbed model possesses asymptotically stable two-cluster periodic solutions. Two biologically motivated forms of random perturbations are considered: bounded variations in growth rate and asymmetric division. We compare the effects of these two dispersive mechanisms with additive Gaussian white noise perturbations. We observe three distinct phases of the response to noise in the model. First, for weak noise there is a linear relationship between the applied noise strength and the dispersion of the clusters. Second, for moderate noise strengths the clusters begin to mix, i.e. individual cells move between clusters, yet the population distribution clearly continues to maintain a two-cluster structure. Third, for strong noise the clusters are destroyed and the population is characterized by a uniform distribution. The second and third phases are separated by an order-disorder phase transition that has the characteristics of a Hopf bifurcation. Furthermore, we show that for the cell cycle model studied, the effects of bounded random perturbations are virtually indistinguishable from those induced by additive Gaussian noise, after appropriate scaling of the variance of noise strength. We then use the model to predict the strength of coupling among the cells from experimental data. In particular, we show that coupling must be rather strong to account for the observed clustering of cells given experimentally estimated noise variance.

Cell cycle dynamics: clustering is universal in negative feedback systems

We study a model of cell cycle ensemble dynamics with cell-cell feedback in which cells in one fixed phase of the cycle S (Signaling) produce chemical agents that affect the growth and development rate of cells that are in another phase R (Responsive). For this type of system there are special periodic solutions that we call k-cyclic or clustered. Biologically, a k-cyclic solution represents k cohorts of synchronized cells spaced nearly evenly around the cell cycle. We show, under very general nonlinear feedback, that for a fixed k the stability of the k-cyclic solutions can be characterized completely in parameter space, a 2 dimensional triangle T. We show that T is naturally partitioned into k(2) sub-triangles on each of which the k-cyclic solutions all have the same stability type. For negative feedback we observe that while the synchronous solution (k = 1) is unstable, regions of stability of k ≥ 2 clustered solutions seem to occupy all of T. We also observe bi-stability or multi-stability for many parameter values in negative feedback systems. Thus in systems with negative feedback we should expect to observe cyclic solutions for some k. This is in contrast to the case of positive feedback, where we observe that the only asymptotically stable periodic orbit is the synchronous solution.

Cell cycle dynamics in a response/signalling feedback system with a gap

We consider a dynamical model of cell cycles of n cells in a culture in which cells in one specific phase (S for signalling) of the cell cycle produce chemical agents that influence the growth/cell cycle progression of cells in another phase (R for responsive). In the case that the feedback is negative, it is known that subpopulations of cells tend to become clustered in the cell cycle; while for a positive feedback, all the cells tend to become synchronized. In this paper, we suppose that there is a gap between the two phases. The gap can be thought of as modelling the physical reality of a time delay in the production and action of the signalling agents. We completely analyse the dynamics of this system when the cells are arranged into two cell cycle clusters. We also consider the stability of certain important periodic solutions in which clusters of cells have a cyclic arrangement and there are just enough clusters to allow interactions between them. We find that the inclusion of a small gap does not greatly alter the global dynamics of the system; there are still large open sets of parameters for which clustered solutions are stable. Thus, we add to the evidence that clustering can be a robust phenomenon in biological systems. However, the gap does effect the system by enhancing the stability of the stable clustered solutions. We explain this phenomenon in terms of contraction rates (Floquet exponents) in various invariant subspaces of the system. We conclude that in systems for which these models are reasonable, a delay in signalling is advantageous to the emergence of clustering.

Clustering in cell cycle dynamics with general response/signaling feedback

Motivated by experimental and theoretical work on autonomous oscillations in yeast, we analyze ordinary differential equations models of large populations of cells with cell-cycle dependent feedback. We assume a particular type of feedback that we call responsive/signaling (RS), but do not specify a functional form of the feedback. We study the dynamics and emergent behavior of solutions, particularly temporal clustering and stability of clustered solutions. We establish the existence of certain periodic clustered solutions as well as "uniform" solutions and add to the evidence that cell-cycle dependent feedback robustly leads to cell-cycle clustering. We highlight the fundamental differences in dynamics between systems with negative and positive feedback. For positive feedback systems the most important mechanism seems to be the stability of individual isolated clusters. On the other hand we find that in negative feedback systems, clusters must interact with each other to reinforce coherence. We conclude from various details of the mathematical analysis that negative feedback is most consistent with observations in yeast experiments.

Measurement of the volume growth rate of single budding yeast with the MOSFET-based microfluidic Coulter counter

We report on measurements of the volume growth rate of ten individual budding yeast cells using a recently developed MOSFET-based microfluidic Coulter counter. The MOSFET-based microfluidic Coulter counter is very sensitive, provides signals that are immune from the baseline drift, and can work with cell culture media of complex composition. These desirable features allow us to directly measure the volume growth rate of single cells of Saccharomyces cerevisiae LYH3865 strain budding yeast in YNB culture media over a whole cell cycle. Results indicate that all budding yeast follow a sigmoid volume growth profile with reduced growth rates at the initial stage before the bud emerges and the final stage after the daughter gets mature. Analysis of the data indicates that even though all piecewise linear, Gomperitz, and Hill's function models can fit the global growth profile equally well, the data strongly support local exponential growth phenomenon. Accurate volume growth measurements are important for applications in systems biology where quantitative parameters are required for modeling and simulation.