Use and abuse of the quasi-steady-state approximation

The transient kinetic behaviour of an open single enzyme, single substrate reaction is examined. The reaction follows the Van Slyke–Cullen mechanism, a spacial case of the Michaelis–Menten reaction. The analysis is performed both with and without applying the quasi-steady-state approximation. The analysis of the full system shows conditions for biochemical pathway coupling, which yield sustained oscillatory behaviour in the enzyme reaction. The reduced model does not demonstrate this behaviour. The results have important implications in the analysis of open biochemical reactions and the modelling of metabolic systems

MHC-Linked Syngeneic Developmental Preference in Thymic Lobes Colonized with Bone Marrow Cells: A Mathematical model

Reconstitution of the T-cell compartment after bone marrow transplantation depends on successful colonization of the thymus by bone-marrow-derived progenitor cells. Recent studies compared the development of syngeneic and allogeneic bone-marrow-derived cells in cocultures with lymphoid-depleted fetal thymus explants, leading to the discovery of MHC-linked syngeneic developmental preference (SDP) in the thymus. To determine the nature of cell interactions among the bone marrow and thymic elements that might underlie SDP, we analyzed this phenomenon by mathematical modeling. The results indicate that syngeneic mature T cells, responsible for inducing this preference, probably interfere both with the seeding of allogeneic bone-marrow-derived thymocyte progenitors in the thymic stroma and with their subsequent proliferation. In addition, the possibility of augmented death among the developing allogeneic thymocytes cannot be ruled out

Scale-Free topologies and Activatory-Inhibitory interactions

A simple model of activatory-inhibitory interactions controlling the activity of agents (substrates) through a "saturated response" dynamical rule in a scale-free network is thoroughly studied. After discussing the most remarkable dynamical features of the model, namely fragmentation and multistability, we present a characterization of the temporal (periodic and chaotic) fluctuations of the quasi-stasis asymptotic states of network activity. The double (both structural and dynamical) source of entangled complexity of the system temporal fluctuations, as an important partial aspect of the Correlation Structure-Function problem, is further discussed to the light of the numerical results, with a view on potential applications of these general results.Comment: Revtex style, 12 pages and 12 figures. Enlarged manuscript with major revision and new results incorporated. To appear in Chaos (2006

Selection of the scaling solution in a cluster coalescence model

The scaling properties of the cluster size distribution of a system of diffusing clusters is studied in terms of a simple kinetic mean field model. It is shown that a one parameter family of mathematically valid scaling solutions exists. Despite this, the kinetics reaches a unique scaling solution independent of initial conditions. This selected scaling solution is marginally physical; i.e., it is the borderline solution between the unphysical and physical branches of the family of solutions.Comment: 4 pages, 5 figure

Mode transitions in a model reaction-diffusion system driven by domain growth and noise

Pattern formation in many biological systems takes place during growth of the underlying domain. We study a specific example of a reaction–diffusion (Turing) model in which peak splitting, driven by domain growth, generates a sequence of patterns. We have previously shown that the pattern sequences which are presented when the domain growth rate is sufficiently rapid exhibit a mode-doubling phenomenon. Such pattern sequences afford reliable selection of certain final patterns, thus addressing the robustness problem inherent of the Turing mechanism. At slower domain growth rates this regular mode doubling breaks down in the presence of small perturbations to the dynamics. In this paper we examine the breaking down of the mode doubling sequence and consider the implications of this behaviour in increasing the range of reliably selectable final patterns

Passive Scalar: Scaling Exponents and Realizability

An isotropic passive scalar field $T$ advected by a rapidly-varying velocity field is studied. The tail of the probability distribution $P(\theta,r)$ for the difference $\theta$ in $T$ across an inertial-range distance $r$ is found to be Gaussian. Scaling exponents of moments of $\theta$ increase as $\sqrt{n}$ or faster at large order $n$, if a mean dissipation conditioned on $\theta$ is a nondecreasing function of $|\theta|$. The $P(\theta,r)$ computed numerically under the so-called linear ansatz is found to be realizable. Some classes of gentle modifications of the linear ansatz are not realizable.Comment: Substantially revised to conform with published version. Revtex (4 pages) with 2 postscript figures. Send email to [email protected]

Square Patterns and Quasi-patterns in Weakly Damped Faraday Waves

Pattern formation in parametric surface waves is studied in the limit of weak viscous dissipation. A set of quasi-potential equations (QPEs) is introduced that admits a closed representation in terms of surface variables alone. A multiscale expansion of the QPEs reveals the importance of triad resonant interactions, and the saturating effect of the driving force leading to a gradient amplitude equation. Minimization of the associated Lyapunov function yields standing wave patterns of square symmetry for capillary waves, and hexagonal patterns and a sequence of quasi-patterns for mixed capillary-gravity waves. Numerical integration of the QPEs reveals a quasi-pattern of eight-fold symmetry in the range of parameters predicted by the multiscale expansion.Comment: RevTeX, 11 pages, 8 figure

Overview of mathematical approaches used to model bacterial chemotaxis I: the single cell

Mathematical modeling of bacterial chemotaxis systems has been influential and insightful in helping to understand experimental observations. We provide here a comprehensive overview of the range of mathematical approaches used for modeling, within a single bacterium, chemotactic processes caused by changes to external gradients in its environment. Specific areas of the bacterial system which have been studied and modeled are discussed in detail, including the modeling of adaptation in response to attractant gradients, the intracellular phosphorylation cascade, membrane receptor clustering, and spatial modeling of intracellular protein signal transduction. The importance of producing robust models that address adaptation, gain, and sensitivity are also discussed. This review highlights that while mathematical modeling has aided in understanding bacterial chemotaxis on the individual cell scale and guiding experimental design, no single model succeeds in robustly describing all of the basic elements of the cell. We conclude by discussing the importance of this and the future of modeling in this area