Mathematical modeling of the intracellular protein dynamics: The importance of active transport along microtubules

Modelling the motion of organelles in an elongated cell via the coordination of heterogeneous drift-diffusion and long-range transport

 Cellular distribution of organelles in living cells is achieved via a variety of transport mechanisms, including directed motion, mediated by molecular motors along microtubules (MTs), and diffusion which is predominantly heterogeneous in space. In this paper, we introduce a model for particle transport in elongated cells that couples poleward drift, long-range bidirectional transport and diffusion with spatial heterogeneity in a three-dimensional space. Using stochastic simulations and analysis of a related population model, we find parameter regions where the three-dimensional model can be reduced to a coupled one-dimensional model or even a one-dimensional scalar model. We explore the efficiency with which individual model components can overcome drift towards one of the cell poles to reach an approximately even distribution. In particular, we find that if lateral movement is well mixed, then increasing the binding ability of particles to MTs is an efficient way to overcome a poleward drift, whereas if lateral motion is not well mixed, then increasing the axial diffusivity away from MTs becomes an efficient way to overcome the poleward drift. Our three-dimensional model provides a new tool that will help to understand the mechanisms by which eukaryotic cells organize their organelles in an elongated cell, and in particular when the one-dimensional models are applicable.

Bridging intracellular scales by mechanistic computational models

The impact of intracellular spatial organization beyond classical compartments on processes such as cell signaling is increasingly recognized. A quantitative, mechanistic understanding of cellular systems therefore needs to account for different scales in at least three coordinates: time, molecular abundances, and space. Mechanistic mathematical models may span all these scales, but corresponding multi-scale models need to resolve mechanistic details on small scales while maintaining computational tractability for larger ones. This typically results in models that combine different levels of description: from a microscopic representation of chemical reactions up to continuum dynamics in space and time. We highlight recent progress in bridging these model classes and outline current challenges in multi-scale models such as active transport and dynamic geometries.

Computational Modelling of Cancer Development and Growth: Modelling at Multiple Scales and Multiscale Modelling

In this paper, we present two mathematical models related to different aspects and scales of cancer growth. The first model is a stochastic spatiotemporal model of both a synthetic gene regulatory network (the example of a three-gene repressilator is given) and an actual gene regulatory network, the NF-[Formula: see text]B pathway. The second model is a force-based individual-based model of the development of a solid avascular tumour with specific application to tumour cords, i.e. a mass of cancer cells growing around a central blood vessel. In each case, we compare our computational simulation results with experimental data. In the final discussion section, we outline how to take the work forward through the development of a multiscale model focussed at the cell level. This would incorporate key intracellular signalling pathways associated with cancer within each cell (e.g. p53-Mdm2, NF-[Formula: see text]B) and through the use of high-performance computing be capable of simulating up to [Formula: see text] cells, i.e. the tissue scale. In this way, mathematical models at multiple scales would be combined to formulate a multiscale computational model.

Spatio-temporal models of synthetic genetic oscillators

Signal transduction pathways play a major role in many important aspects of cellular function e.g. cell division, apoptosis. One important class of signal transduction pathways is gene regulatory networks (GRNs). In many GRNs, proteins bind to gene sites in the nucleus thereby altering the transcription rate. Such proteins are known as transcription factors. If the binding reduces the transcription rate there is a negative feedback leading to oscillatory behaviour in mRNA and protein levels, both spatially (e.g. by observing fluorescently labelled molecules in single cells) and temporally (e.g. by observing protein/mRNA levels over time). Recent computational modelling has demonstrated that spatial movement of the molecules is a vital component of GRNs and may cause the oscillations. These numerical findings have subsequently been proved rigorously i.e. the diffusion coefficient of the protein/mRNA acts as a bifurcation parameter and gives rise to a Hopf bifurcation. In this paper we first present a model of the canonical GRN (the Hes1 protein) and show the effect of varying the spatial location of gene and protein production sites on the oscillations. We then extend the approach to examine spatio-temporal models of synthetic gene regulatory networks e.g. n-gene repressilators and activator-repressor systems.

On the interplay between mathematics and biology: hallmarks toward a new systems biology

This paper proposes a critical analysis of the existing literature on mathematical tools developed toward systems biology approaches and, out of this overview, develops a new approach whose main features can be briefly summarized as follows: derivation of mathematical structures suitable to capture the complexity of biological, hence living, systems, modeling, by appropriate mathematical tools, Darwinian type dynamics, namely mutations followed by selection and evolution. Moreover, multiscale methods to move from genes to cells, and from cells to tissue are analyzed in view of a new systems biology approach.