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

Spatial and spatiotemporal pattern formation in generalised Turing systems

Reaction-diffusion, or Turing, models have been proposed to account for a number of pattern formation phenomena in early development. However, there are a number of crucial morphogenetic phenomena that contradict the standard Turing model. Here we review three generalisations of the Turing model and show how they can be applied to two such phenomena. We discuss how these generalisations can provide insight to the processes underlying patterning in these cases

Lattice Gas Automata for Reactive Systems

Reactive lattice gas automata provide a microscopic approachto the dynamics of spatially-distributed reacting systems. After introducing the subject within the wider framework of lattice gas automata (LGA) as a microscopic approach to the phenomenology of macroscopic systems, we describe the reactive LGA in terms of a simple physical picture to show how an automaton can be constructed to capture the essentials of a reactive molecular dynamics scheme. The statistical mechanical theory of the automaton is then developed for diffusive transport and for reactive processes, and a general algorithm is presented for reactive LGA. The method is illustrated by considering applications to bistable and excitable media, oscillatory behavior in reactive systems, chemical chaos and pattern formation triggered by Turing bifurcations. The reactive lattice gas scheme is contrasted with related cellular automaton methods and the paper concludes with a discussion of future perspectives.Comment: to appear in PHYSICS REPORTS, 81 revtex pages; uuencoded gziped postscript file; figures available from [email protected] or [email protected]

A multiscale hybrid model for pro-angiogenic calcium signals in a vascular endothelial cell

Cytosolic calcium machinery is one of the principal signaling mechanisms by which endothelial cells (ECs) respond to external stimuli during several biological processes, including vascular progression in both physiological and pathological conditions. Low concentrations of angiogenic factors (such as VEGF) activate in fact complex pathways involving, among others, second messengers arachidonic acid (AA) and nitric oxide (NO), which in turn control the activity of plasma membrane calcium channels. The subsequent increase in the intracellular level of the ion regulates fundamental biophysical properties of ECs (such as elasticity, intrinsic motility, and chemical strength), enhancing their migratory capacity. Previously, a number of continuous models have represented cytosolic calcium dynamics, while EC migration in angiogenesis has been separately approached with discrete, lattice-based techniques. These two components are here integrated and interfaced to provide a multiscale and hybrid Cellular Potts Model (CPM), where the phenomenology of a motile EC is realistically mediated by its calcium-dependent subcellular events. The model, based on a realistic 3-D cell morphology with a nuclear and a cytosolic region, is set with known biochemical and electrophysiological data. In particular, the resulting simulations are able to reproduce and describe the polarization process, typical of stimulated vascular cells, in various experimental conditions.Moreover, by analyzing the mutual interactions between multilevel biochemical and biomechanical aspects, our study investigates ways to inhibit cell migration: such strategies have in fact the potential to result in pharmacological interventions useful to disrupt malignant vascular progressio

Oscillating focus of SopA associated with filamentous structure guides partitioning of F plasmid

The F plasmid is actively partitioned to daughter cells by the sopABC gene. To elucidate the partitioning mechanisms, we simultaneously analysed movements of the plasmid and the SopA ATPase in single living cells. SopA, which is a putative motor protein assembled densely near nucleoid borders and formed a single discrete focus associated with less dense filamentous distribution along the long axis of the cell. The dense SopA focus oscillates between cell poles. The direction of the plasmid motion switches as the SopA focus switches its position. The velocity of the plasmid motion stays constant while it oscillates moving towards the SopA focus. The low density filamentous distribution of SopA persisted throughout the SopA oscillation. The focus associated with filamentous distribution of SopA was also observed in a cell without nucleoid. The SopA filament may guide the movement of the plasmid as a railway track and lead it to cell quarters

Filament Depolymerization Can Explain Chromosome Pulling during Bacterial Mitosis

Chromosome segregation is fundamental to all cells, but the force-generating mechanisms underlying chromosome translocation in bacteria remain mysterious. Caulobacter crescentus utilizes a depolymerization-driven process in which a ParA protein structure elongates from the new cell pole, binds to a ParB-decorated chromosome, and then retracts via disassembly, pulling the chromosome across the cell. This poses the question of how a depolymerizing structure can robustly pull the chromosome that disassembles it. We perform Brownian dynamics simulations with a simple, physically consistent model of the ParABS system. The simulations suggest that the mechanism of translocation is “self-diffusiophoretic”: by disassembling ParA, ParB generates a ParA concentration gradient so that the ParA concentration is higher in front of the chromosome than behind it. Since the chromosome is attracted to ParA via ParB, it moves up the ParA gradient and across the cell. We find that translocation is most robust when ParB binds side-on to ParA filaments. In this case, robust translocation occurs over a wide parameter range and is controlled by a single dimensionless quantity: the product of the rate of ParA disassembly and a characteristic relaxation time of the chromosome. This time scale measures the time it takes for the chromosome to recover its average shape after it is has been pulled. Our results suggest explanations for observed phenomena such as segregation failure, filament-length-dependent translocation velocity, and chromosomal compaction

Diffusive coupling can discriminate between similar reaction mechanisms in an allosteric enzyme system

<p>Abstract</p> <p>Background</p> <p>A central question for the understanding of biological reaction networks is how a particular dynamic behavior, such as bistability or oscillations, is realized at the molecular level. So far this question has been mainly addressed in well-mixed reaction systems which are conveniently described by ordinary differential equations. However, much less is known about how molecular details of a reaction mechanism can affect the dynamics in diffusively coupled systems because the resulting partial differential equations are much more difficult to analyze.</p> <p>Results</p> <p>Motivated by recent experiments we compare two closely related mechanisms for the product activation of allosteric enzymes with respect to their ability to induce different types of reaction-diffusion waves and stationary Turing patterns. The analysis is facilitated by mapping each model to an associated complex Ginzburg-Landau equation. We show that a sequential activation mechanism, as implemented in the model of Monod, Wyman and Changeux (MWC), can generate inward rotating spiral waves which were recently observed as glycolytic activity waves in yeast extracts. In contrast, in the limiting case of a simple Hill activation, the formation of inward propagating waves is suppressed by a Turing instability. The occurrence of this unusual wave dynamics is not related to the magnitude of the enzyme cooperativity (as it is true for the occurrence of oscillations), but to the sensitivity with respect to changes of the activator concentration. Also, the MWC mechanism generates wave patterns that are more stable against long wave length perturbations.</p> <p>Conclusions</p> <p>This analysis demonstrates that amplitude equations, which describe the spatio-temporal dynamics near an instability, represent a valuable tool to investigate the molecular effects of reaction mechanisms on pattern formation in spatially extended systems. Using this approach we have shown that the occurrence of inward rotating spiral waves in glycolysis can be explained in terms of an MWC, but not with a Hill mechanism for the activation of the allosteric enzyme phosphofructokinase. Our results also highlight the importance of enzyme oligomerization for a possible experimental generation of Turing patterns in biological systems.</p