Regulation of differentiation in a population of cells interacting through a common pool

We consider a model of a suspension of a cell population in a well-mixed medium. There are two chemical substances, say A and H, reacting in each cell of the population and the substance H can only diffuse from the inside of cell to the medium or vice versa across the cell membrane. The medium is well mixed that the concentration of H is kept uniform over the medium. Cells interact indirectly with each other through the medium. The differential equations governing the dynamics of the suspension are analyzed using standard techniques for differential equations. It is shown that the cell population divides into several groups in respect of the chemical concentrations as time elapses. It is also shown how the fraction of the number of cells belonging to each subgroup to the total number of cells is regulated. The results may be used to explain the mechanism for differentiation of multi-cellular organisms.     

Dynamics of chemostat in which one microbial population grows on multiple complementary nutrients

The equations of a chemostat in which one microbial population grows on multiple rate-limiting nutrients are formulated. The dynamics of a chemostat involving growth on complementary nutrients is studied through stability analysis of the system of equations. Some conditions are derived that relate the dynamic behavior of the chemostat to its operating conditions and can be applied to any model for the specific growth rate of the population. It is shown that, if maintenance of the population is neglected, the system exhibits no sustained or damped oscillations. If maintenance of the population is considered, damped oscillations are observed for some operating conditions.     

A 3D Hybrid Model for Tissue Growth: The Interplay between Cell Population and Mass Transport Dynamics

To provide theoretical guidance for the design and in vitro cultivation of bioartificial tissues, we have developed a multiscale computational model that can describe the complex interplay between cell population and mass transport dynamics that governs the growth of tissues in three-dimensional scaffolds. The model has three components: a transient partial differential equation for the simultaneous diffusion and consumption of a limiting nutrient; a cellular automaton describing cell migration, proliferation, and collision; and equations that quantify how the varying nutrient concentration modulates cell division and migration. The hybrid discrete-continuous model was parallelized and solved on a distributed-memory multicomputer to study how transport limitations affect tissue regeneration rates under conditions encountered in typical bioreactors. Simulation results show that the severity of transport limitations can be estimated by the magnitude of two dimensionless groups: the Thiele modulus and the Biot number. Key parameters including the initial seeding mode, cell migration speed, and the hydrodynamic conditions in the bioreactor are shown to affect not only the overall rate, but also the pattern of tissue growth. This study lays the groundwork for more comprehensive models that can handle mixed cell cultures, multiple nutrients and growth factors, and other cellular processes, such as cell death.     

The age structure of populations of Saccharomyces cerevisiae

A discrete deterministic model is described for the growth of an age-structured population of yeast, Saccharomyces cerevisiae, incorporating recent information on the asymmetry of cell division and control of the cell cycle in this species. Solutions are obtained for the age structure of the population at equilibrium, and for the equilibrium distribution of relative frequency of cells through the cell cycle. The model is applied to experimental data on the changing age structure of nonequilibrium populations of yeast. The model predicts well both the transient behavior and the equilibrium structure of such populations. It is shown that the asymmetry of cell division explains (1) the excess of newly formed daughter cells in the population as compared to the frequency of older cells and (2) the damped oscillations in the frequencies of cells of different ages as demographic equilibrium is approached.     

Complementary replication.

Differential equations for the kinetics of complementary replicating macromolecules in a flow reactor are derived. It is shown that such a model has many features in common with the differential equation for direct replication, the replicator equation. Two special cases of replication, and the influence of mutation on them, have been studied in detail. In the case of first-order mass action kinetics--the quasi-species model--complementary replication, like direct replication, exhibits an error threshold for the replication accuracy, below which the genetic information is lost. In turns out that the long-time behavior of many special cases of the second-order kinetics model can be described in terms of second-order replicator equations, although this is not possible in general.     

How Chlamydomonas keeps track of the light once it has reached the right phototactic orientation.

By using a real-time assay that allows measurement of the phototactic orientation of the unicellular alga Chlamydomonas with millisecond time resolution, it can be shown that single photons not only induce transient direction changes but that fluence rates as low as 1 photon cell(-1) s(-1) can already lead to a persistent orientation. Orientation is a binary variable, i.e., in a partially oriented population some organisms are fully oriented while the rest are still at random. Action spectra reveal that the response to a pulsed stimulus follows the Dartnall-nomogram for a rhodopsin while the response to a persistent stimulus falls off more rapidly toward the red end of the spectrum. Thus light of 540 nm, for which chlamy-rhodopsin is equally sensitive as for 440-nm light, induces no measurable persistent orientation while 440-nm light does. A model is presented which explains not only this behavior, but also how Chlamydomonas can track the light direction and switches between a positive and negative phototaxis. According to the model the ability to detect the direction of light, to make the right turn and to stay oriented, is a direct consequence of the helical path of the organism, the orientation of its eyespot relative to the helix-axis, and the special shielding properties of eyespot and cell body. The model places particular emphasis on the fact that prolonged swimming into the correct direction not only requires making a correct turn initially, but also avoiding further turns once the right direction has been reached.     

The role of seasonality in the dynamics of deer tick populations

In this paper, we formulate a nonlinear system of difference equations that models the three-stage life cycle of the deer tick over four seasons. We study the effect of seasonality on the stability and oscillatory behavior of the tick population by comparing analytically the seasonal model with a non-seasonal one. The analysis of the models reveals the existence of two equilibrium points. We discuss the necessary and sufficient conditions for local asymptotic stability of the equilibria and analyze the boundedness and oscillatory behavior of the solutions. A main result of the mathematical analysis is that seasonality in the life cycle of the deer tick can have a positive effect, in the sense that it increases the stability of the system. It is also shown that for some combination of parameters within the stability region, perturbations will result in a return to the equilibrium through transient oscillations. The models are used to explore the biological consequences of parameter variations reflecting expected environmental changes.     

Mechanisms Explaining Transitions between Tonic and Phasic Firing in Neuronal Populations as Predicted by a Low Dimensional Firing Rate Model

Several firing patterns experimentally observed in neural populations have been successfully correlated to animal behavior. Population bursting, hereby regarded as a period of high firing rate followed by a period of quiescence, is typically observed in groups of neurons during behavior. Biophysical membrane-potential models of single cell bursting involve at least three equations. Extending such models to study the collective behavior of neural populations involves thousands of equations and can be very expensive computationally. For this reason, low dimensional population models that capture biophysical aspects of networks are needed. The present paper uses a firing-rate model to study mechanisms that trigger and stop transitions between tonic and phasic population firing. These mechanisms are captured through a two-dimensional system, which can potentially be extended to include interactions between different areas of the nervous system with a small number of equations. The typical behavior of midbrain dopaminergic neurons in the rodent is used as an example to illustrate and interpret our results. The model presented here can be used as a building block to study interactions between networks of neurons. This theoretical approach may help contextualize and understand the factors involved in regulating burst firing in populations and how it may modulate distinct aspects of behavior.     

Pair formation models with maturation period

The standard model for pair formation is generalized to include a maturation period. This model in the form of three coupled delay equations is a special case of the general age-structured model for a two-sex population. The exact conditions for the existence of an exponential (persistent) two-sex solution are derived. It is shown that this solution is unique and locally stable. In order to achieve these results the theory of homogeneous differential equations is extended to a class of homogeneous delay equations.     

Dynamics of yeast metabolism and regulation

The dynamic behavior of yeasts is discussed on the basis of transient experiments, such as pulses and shifts in continuous culture, and of oscillating synchronized cultures. The minimal elements of a structured model are evaluated and extensions by regulatory mechanisms are proposed. Segregation of a population in cell classes reflecting the position of a cell in the cell cycle, its individual age, is shown to be necessary in order to account for the spontaneous synchronization of continuous S. cerevisisae cultures.     

A model of the neuronal network using eye micromovements for distinguishing contrasts on an image

A monochromatic model of a neuron network of foveal retina with the straight way of information stream of the system "cone-diminutive bipolar-diminutive ganglion cell" is proposed. The network can distinguish contrasts on an image. It is shown that cells of the output level have concentric receptive fields formed due to eye micromovements. The paper discusses the functioning of the model in the case when eye movements contain resetting microsaccades as well as drift. A special attention is paid to the compensation of transient processes.     

Outline of a mathematical biosociology of beliefs and prejudices

Behavior based on acceptance on faith is considered as a result of a special conditioning, and is assumed to be mediated through a special center. Behavior based on rational critical analysis is mediated through a different center. The previously developed theory of two mutually exclusive behaviors is applied to the situation. Making simple assumptions about the development of the two centers with age, it is shown how equations can be derived which describe the percentage of individuals exhibiting one type of behavior or another as a function of time and of different biosociological parameters.     

Nonlinear viscoelastic, thermodynamically consistent, models for biological soft tissue

The mechanical behavior of most biological soft tissue is nonlinear viscoelastic rather than elastic. Many of the models previously proposed for soft tissue involve ad hoc systems of springs and dashpots or require measurement of time-dependent constitutive coefficient functions. The model proposed here is a system of evolution differential equations, which are determined by the long-term behavior of the material as represented by an energy function of the type used for elasticity. The necessary empirical data is time independent and therefore easier to obtain. These evolution equations, which represent non-equilibrium, transient responses such as creep, stress relaxation, or variable loading, are derived from a maximum energy dissipation principle, which supplements the second law of thermodynamics. The evolution model can represent both creep and stress relaxation, depending on the choice of control variables, because of the assumption that a unique long-term manifold exists for both processes. It succeeds, with one set of material constants, in reproducing the loading-unloading hysteresis for soft tissue. The models are thermodynamically consistent so that, given data, they may be extended to the temperature-dependent behavior of biological tissue, such as the change in temperature during uniaxial loading. The Holzapfel et al. three-dimensional two-layer elastic model for healthy artery tissue is shown to generate evolution equations by this construction for biaxial loading of a flat specimen. A simplified version of the Shah-Humphrey model for the elastodynamical behavior of a saccular aneurysm is extended to viscoelastic behavior.     

Unsteady-state operation of continuous fermentor for enhancement of cell mass production

The transient behavior of continuous fermentation is studied concentrating on the time scale intrinsic to the system. The time scale is the time required for the fermentorto reach a stable steady state after the disturbance of cell mass is introduced. When the cell concentration is disturbed from the steady-state value, in particular, at the dilution rate near washout, the transient period becomes extended significantly, and the steady state is resumed sluggishly. This sluggish transient behavior could be turned to an advantage for enhancing the cell mass output rate. The proposed transient operation is a continuous fermentation whereby a positive disturbance in the cell mass is introduced, so that the cell concentration is higher than the steady-state value for an extended transient period. It is shown that a significantly higher cell mass production than that from the steady-state continuous fermentation can be achieved. Simple experiments were performed to demonstrate the improvement of cell (Candida utilis) productivity.     

The effects of the functional response on the bifurcation behavior of a mite predator–prey interaction model

 A predator–prey interaction model based on a system of differential equations with temperature-dependent parameters chosen appropriately for a mite interaction on apple trees is analyzed to determine how the type of functional response influences bifurcation and stability behavior. Instances of type I, II, III, and IV functional responses are considered, the last of which incorporates prey interference with predation. It is shown that the model systems with the type I, II, and III functional responses exhibit qualitatively similar bifurcation and stability behavior over the interval of definition of the temperature parameter. Similar behavior is found in the system with the type IV functional response at low levels of prey interference. Higher levels of interference are destabilizing, as illustrated by the prevalence of bistability and by the presence of three attractors for some values of the model parameters. All four systems are capable of modeling population oscillations and outbreaks. Received 12 March 1996; received in revised form 25 October 1996     

Microbial growth dynamics on the basis of individual budgets

The popular theories for microbial dynamics by Monod, Pirt and Droop are shown to be special cases of a model for individual budgets, in which growth and maintenance are on the expense of reserve materials. The dynamics of reserve materials is a first order process with a relaxation time proportional to cell length; maintenance is proportional to cell volume, and uptake, which depends hyperbolically on substrate density, is proportional to cell volume as well. Because of the latter, population dynamics depends on the behaviour of the individuals in a simple way, such that the cell volume distribution has no quantitative effect.When uptake is proportional to the surface area of the cell, which is realistic from a physical point of view, the relation between the individual level and the population one becomes more complicated and the cell size and shape distribution affects population dynamics. It is shown how the changing shape of rods modifies uptake and, consequently, growth.The concept of energy conductance, defined as the ratio, of the maximum surface area specific uptake and the volume specific energy reserve has been introduced in the analysis of microbial dynamics. The first tentative results indicate that the value for E. coli is close to the mean value for a wide variety of animals.Properties of the model for cell suspension at constant substrate densities are analyzed and tested against a variety of experimental data from the literature on both the individual and the population level.     

Mathematical description of a bursting pacemaker neuron by a modification of the Hodgkin-Huxley equations.

Modifications based on experimental results reported in the literature are made to the Hodgkin-Huxley equations to describe the electrophysiological behavior of the Aplysia abdominal ganglion R15 cell. The system is then further modified to describe the effects with the application of the drug tetrodotoxin (TTX) to the cells' bathing medium. Methods of the qualitative theory of differential equations are used to determine the conditions necessary for such a system of equations to have an oscillatory solution. A model satisfying these conditions is shown to preduct many experimental observations of R15 cell behavior. Numerical solutions are obtained for differential equations satisfying the conditions of the model. These solutions are shown to have a form similar to that of the bursting which is characteristic of this cell, and to preduct many results of experiments conducted on this cell. The physiological implications of the model are discussed.     

Stochastic Fluctuations Through Intrinsic Noise in Evolutionary Game Dynamics

A one-step (birth–death) process is used to investigate stochastic noise in an elementary two-phenotype evolutionary game model based on a payoff matrix. In this model, we assume that the population size is finite but not fixed and that all individuals have, in addition to the frequency-dependent fitness given by the evolutionary game, the same background fitness that decreases linearly in the total population size. Although this assumption guarantees population extinction is a globally attracting absorbing barrier of the Markov process, sample trajectories do not illustrate this result even for relatively small carrying capacities. Instead, the observed persistent transient behavior can be analyzed using the steady-state statistics (i.e., mean and variance) of a stochastic model for intrinsic noise that assumes the population does not go extinct. It is shown that there is good agreement between the theory of these statistics and the simulation results. Furthermore, the ESS of the evolutionary game can be used to predict the mean steady state.     

A model of tumor growth based on cell cycle kinetics

This work describes a mathematical model of growth based on the kinetics of the cell cycle. A traditional model of the cell cycle has been used, with the addition of a resting (G₀) state from which cells could reenter the reproductive cycle. The model assumes that a growth regulatory substance regulates the transition of cells to and from the resting state. Other transitions between the phases of the cycle were modeled as a first order process. Cell loss is an important feature of growth kinetics, and has been represented by a general but tractable mathematical form. The resulting model forms a system of ordinary nonlinear differential equations. Analytic methods are employed first in the study of this system. Simplifying assumptions regarding cell loss give rise to special cases for which equilibrium solutions can be found. One special case, which assumes first order loss from all cell cycle phases at equal rates, is presented here. For small time values, approximations corresponding to exponential growth were developed. The equations describing an intrinsic growth rate were derived. Simulation methods were used to further characterize the behavior of this model. Parameter values were chosen based on animal tumor cell cycle kinetic data, resulting in a set of 45 model simulations. Several tumor treatment protocols were simulated which illustrated the importance of the intrinsic growth rate and cell loss concepts. Although the qualitative behavior regarding absolute and relative growth is reasonable, this model awaits data for model fitting, parameter estimation, or revision of the equations.