Quantitative Biology

The number of ribosomes in a cell is considered as limiting, and gene expression is thus largely determined by their cellular concentration. In this work we develop a toy model to study the trade-off between the ribosomal supply and the demand of the translation machinery, dictated by the composition of the transcript pool. Our equilibrium framework is useful to highlight qualitative behaviours and new means of gene expression regulation determined by the fine balance of this trade-off. We also speculate on the possible impact of these mechanisms on cellular physiology.

Mathematical modeling is now used commonly in the analysis of signaling networks. With advances in high resolution microscopy, the spatial location of different signaling molecules and the spatio-temporal dynamics of signaling microdomains are now widely acknowledged as key features of biochemical signal transduction. Reaction-diffusion mechanisms are commonly used to model such features, often with a heavy reliance on numerical simulations to obtain results. However, simulations are parameter dependent and may not be able to provide an understanding of the full range of the system responses. Analytical approaches on the other hand provide a framework to study the entire phase space. In this tutorial, we provide a largely analytical method for studying reaction-diffusion models and analyzing their stability properties. Using two representative biological examples, we demonstrate how this approach can guide experimental design.

Binding, lateral diffusion and exchange are fundamental dynamic processes involved in protein association with cellular membranes. In this study, we developed numerical simulations of lateral diffusion and exchange of fluorophores in membranes with arbitrary bleach geometry and exchange of the membrane localized fluorophore with the cytosol during Fluorescence Recovery after Photobleaching (FRAP) experiments. The model simulations were used to design FRAP experiments with varying bleach region sizes on plasma-membrane localized wild type GFP-Ras2 with a dual lipid anchor and mutant GFP-Ras2C318S with a single lipid anchor in live yeast cells to investigate diffusional mobility and the presence of any exchange processes operating in the time scale of our experiments. Model parameters estimated using data from FRAP experiments with a 1 micron x 1 micron bleach region-of-interest (ROI) and a 0.5 micron x 0.5 micron bleach ROI showed that GFP-Ras2, single or dual lipid modified, diffuses as single species with no evidence of exchange with a cytoplasmic pool. This is the first report of Ras2 mobility in yeast plasma membrane. The methods developed in this study are generally applicable for studying diffusion and exchange of membrane associated fluorophores using FRAP on commercial confocal laser scanning microscopes.

We propose and analyze a mathematical model of cargo transport by non-processive molecular motors. In our model, the motors change states by random discrete events (corresponding to stepping and binding/unbinding), while the cargo position follows a stochastic differential equation (SDE) that depends on the discrete states of the motors. The resulting system for the cargo position is consequently an SDE that randomly switches according to a Markov jump process governing motor dynamics. To study this system we (1) cast the cargo position in a renewal theory framework and generalize the renewal reward theorem and (2) decompose the continuous and discrete sources of stochasticity and exploit a resulting pair of disparate timescales. With these mathematical tools, we obtain explicit formulas for experimentally measurable quantities, such as cargo velocity and run length. Analyzing these formulas then yields some predictions regarding so-called non-processive clustering, the phenomenon that a single motor cannot transport cargo, but that two or more motors can. We find that having motor stepping, binding, and unbinding rates depend on the number of bound motors, due to geometric effects, is necessary and sufficient to explain recent experimental data on non-processive motors.

Microtubule (MT) radial arrays or asters establish the internal topology of a cell by interacting with organelles and molecular motors. We proceed to understand the general pattern forming potential of aster-motor systems using a computational model of multiple MT asters interacting with motors in a cellular confinement. In this model dynein motors are attached to the cell cortex and plus-ended motors resembling kinesin-5 diffuse in the cell interior. The introduction of 'noise' in the form of MT length fluctuations spontaneously results in the emergence of coordinated, achiral vortex-like rotation of asters. The coherence and persistence of rotation requires a threshold density of both cortical dyneins and coupling kinesins, while the onset of rotation if diffusion-limited with relation to cortical dynein mobility. The coordinated rotational motion arises due to the resolution of the 'tug-of-war' of the rotational component due to cortical motors by 'noise' in the form of MT dynamic instability, and such transient symmetry breaking is amplified by local coupling by kinesin-5 complexes. The lack of widespread aster rotation across cell types suggests biophysical mechanisms that suppress such intrinsic dynamics may have evolved. This model is analogous to more general models of locally coupled self-propelled particles (SPP) that spontaneously undergo collective transport in presence of 'noise' that have been invoked to explain swarming in birds and fish. However, the aster-motor system is distinct from SPP models with regard to particle density and 'noise' dependence, providing a set of experimentally testable predictions for a novel sub-cellular pattern forming system.

The asymmetrical inheritance of plasmid DNA, as well as other cellular components, has been shown to be involved in replicative aging. In Saccharomyces cerevisiae, there is an ongoing debate regarding the mechanisms underlying this important asymmetry. Currently proposed models suggest it is established via diffusion, but differ on whether a diffusion barrier is necessary or not. However, no study so far incorporated key aspects to segregation, such as dynamic morphology changes throughout anaphase or plasmids size. Here, we determine the distinct effects and contributions of individual cellular variability, plasmid volume and moving boundaries in the asymmetric segregation of plasmids. We do this by measuring cellular nuclear geometries and plasmid diffusion rates with confocal microscopy, subsequently incorporating this data into a growing domain stochastic spatial simulator. Our modelling and simulations confirms that plasmid asymmetrical inheritance does not require an active barrier to diffusion, and provides a full analysis on plasmid size effects.

The first of N identical independently distributed (i.i.d.) Brownian trajectories that arrives to a small target, sets the time scale of activation, which in general is much faster than the arrival to the target of only a single trajectory. Analytical asymptotic expressions for the minimal time is notoriously difficult to compute in general geometries. We derive here asymptotic laws for the probability density function of the first and second arrival times of a large number of i.i.d. Brownian trajectories to a small target in 1,2, and 3 dimensions and study their range of validity by stochastic simulations. The results are applied to activation of biochemical pathways in cellular transduction.

Due to the stochastic nature of biochemical processes, the copy number of any given type of molecule inside a living cell often exhibits large temporal fluctuations. Here, we develop analytic methods to investigate how the noise arising from a bursting input is reshaped by a transport reaction which is either linear or of the Michaelis-Menten type. A slow transport rate smoothes out fluctuations at the output end and minimizes the impact of bursting on the downstream cellular activities. In the context of gene expression in eukaryotic cells, our results indicate that transcriptional bursting can be substantially attenuated by the transport of mRNA from nucleus to cytoplasm. Saturation of the transport mediators or nuclear pores contributes further to the noise reduction. We suggest that the mRNA transport should be taken into account in the interpretation of relevant experimental data on transcriptional bursting.

Symmetry breaking, a central principle of physics, has been hailed as the driver of self-organization in biological systems in general and biogenesis of cellular organelles in particular, but the molecular mechanisms of symmetry breaking only begin to become understood. Centrioles, the structural cores of centrosomes and cilia, must duplicate every cell cycle to ensure their faithful inheritance through cellular divisions. Work in model organisms identified conserved proteins required for centriole duplication and found that altering their abundance affects centriole number. However, the biophysical principles that ensure that, under physiological conditions, only a single procentriole is produced on each mother centriole remain enigmatic. Here we propose a mechanistic biophysical model for the initiation of procentriole formation in mammalian cells. The model faithfully recapitulates the experimentally observed transition from PLK4 uniformly distributed around the mother centriole, the "ring", to a unique PLK4 focus, the "spot", that triggers the assembly of a new procentriole. This symmetry breaking requires a dual positive feedback based on autocatalytic activation of PLK4 and enhanced centriolar anchoring of PLK4-STIL complexes by phosphorylated STIL. We find that, contrary to previous proposals, in situ degradation of active PLK4 is insufficient to break symmetry. Instead, the model predicts that competition between transient PLK4 activity maxima for PLK4-STIL complexes explains both the instability of the PLK4 ring and formation of the unique PLK4 spot. In the model, strong competition at physiologically normal parameters robustly produces a single procentriole, while increasing overexpression of PLK4 and STIL weakens the competition and causes progressive addition of procentrioles in agreement with experimental observations.

Staphylococcus aureus is a coccal bacterium, which divides by binary fission. After division the cells remain attached giving rise to small clusters, with a characteristic 'bunch of grapes' morphology. S. aureus is an important human pathogen and this, combined with the increasing prevalence of antibiotic-resistant strains, such as Methicillin Resistant S. aureus (MRSA), make it an excellent subject for studies of new methods of antimicrobial action. Many antibiotics, such as penicillin, prevent S. aureus cell division and so an understanding of this fundamental process may pave the way to the identification of novel drugs. We present here a novel image analysis framework for automated detection and segmentation of cells in S. aureus clusters, and identification of their cell division planes. We demonstrate the technique on GFP labelled EzrA, a protein that localises to a mid-cell plane during division and is involved in regulation of cell size and division. The algorithms may have wider applicability in detecting morphologically complex structures of fluorescently-labelled proteins within cells in other cell clusters.