An introduction to mathematical models of coagulation-fragmentation processes: a discrete deterministic mean-field approach

We summarise the properties and the fundamental mathematical results associated with basic models which describe coagulation and fragmentation processes in a deterministic manner and in which cluster size is a discrete quantity (an integer multiple of some basic unit size). In particular, we discuss Smoluchowski's equation for aggregation, the Becker-Döring model of simultaneous aggregation and fragmentation, and more general models involving coagulation and fragmentation.

Elucidating the interactions between the adhesive and transcriptioanl functions of beta-catenin in normal and cancerous cells

Wnt signalling is involved in a wide range of physiological and pathological processes. The presence of an extracellular Wnt stimulus induces cytoplasmic stabilisation and nuclear translocation of beta-catenin, a protein that also plays an essential role in cadherin-mediated adhesion. Two main hypotheses have been proposed concerning the balance between beta-catenin's adhesive and transcriptional functions: either beta-catenin's fate is determined by competition between its binding partners, or Wnt induces folding of beta-catenin into a conformation allocated preferentially to transcription. The experimental data supporting each hypotheses remain inconclusive. In this paper we present a new mathematical model of the Wnt pathway that incorporates beta-catenin's dual function. We use this model to carry out a series of in silico experiments and compare the behaviour of systems governed by each hypothesis. Our analytical results and model simulations provide further insight into the current understanding of Wnt signalling and, in particular, reveal differences in the response of the two modes of interaction between adhesion and signalling in certain in silico settings. We also exploit our model to investigate the impact of the mutations most commonly observed in human colorectal cancer. Simulations show that the amount of functional APC required to maintain a normal phenotype increases with increasing strength of the Wnt signal, a result which illustrates that the environment can substantially influence both tumour initiation and phenotype.

A bidomain threshold model of propagating calcium waves

We present a bidomain fire-diffuse-fire model that facilitates mathematical analysis of propagating waves of elevated intracellular calcium (Ca) in living cells. Modelling Ca release as a threshold process allows the explicit construction of travelling wave solutions to probe the dependence of Ca wave speed on physiologically important parameters such as the threshold for Ca release from the endoplasmic reticulum (ER) to the cytosol, the rate of Ca resequestration from the cytosol to the ER, and the total [Ca] (cytosolic plus ER). Interestingly, linear stability analysis of the bidomain fire-diffuse-fire model predicts the onset of dynamic wave instabilities leading to the emergence of Ca waves that propagate in a back-and-forth manner. Numerical simulations are used to confirm the presence of these so-called "tango waves" and the dependence of Ca wave speed on the total [Ca]. The original publication is available at www.springerlink.com (Journal of Mathematical Biology)

Coarse-graining and renormalisation group methods for the elucidation of the kinetics of complex nucleation and growth processes

We review our work on generalisations of the Becker-Doring model of cluster-formation as applied to nucleation theory, polymer growth kinetics, and the formation of upramolecular structures in colloidal chemistry. One valuable tool in analysing mathematical models of these systems has been the coarse-graining approximation which enables macroscopic models for observable quantities to be derived from microscopic ones. This permits assumptions about the detailed molecular mechanisms to be tested, and their influence on the large-scale kinetics of surfactant self-assembly to be elucidated. We also summarise our more recent results on Becker-Doring systems, notably demonstrating that cross-inhibition and autocatalysis can destabilise a uniform solution and lead to a competitive environment in which some species flourish at the expense of others, phenomena relevant in models of the origins of life.

Frequency of elemental events of intracellular Ca²⁺ dynamics

The dynamics of intracellular Ca²⁺ is driven by random events called Ca²⁺ puffs, in which Ca²⁺ is liberated from intracellular stores. We show that the emergence of Ca²⁺ puffs can be mapped to an escape process. The mean first passage times that correspond to the stochastic fraction of puff periods are computed from a novel master equation and two Fokker-Planck equations. Our results demonstrate that the mathematical modeling of Ca²⁺ puffs has to account for the discrete character of the Ca²⁺ release sites and does not permit a continuous description of the number of open channels.