Biomedical modeling: the role of transport and mechanics

This issue contains a series of papers that were invited following a workshop held in July 2011 at the University of Notre Dame London Center. The goal of the workshop was to present the latest advances in theory, experimentation, and modeling methodologies related to the role of mechanics in biological systems. Growth, morphogenesis, and many diseases are characterized by time dependent changes in the material properties of tissues—affected by resident cells—that, in turn, affect the function of the tissue and contribute to, or mitigate, the disease. Mathematical modeling and simulation are essential for testing and developing scientific hypotheses related to the physical behavior of biological tissues, because of the complex geometries, inhomogeneous properties, rate dependences, and nonlinear feedback interactions that it entails

Interpreting two-photon imaging data of lymphocyte motility

Recently, using two-photon imaging it has been found that the movement of B and T cells in lymph nodes can be described by a random walk with persistence of orientation in the range of 2 minutes. We interpret this new class of lymphocyte motility data within a theoretical model. The model considers cell movement to be composed of the movement of subunits of the cell membrane. In this way movement and deformation of the cell are correlated to each other. We find that, indeed, the lymphocyte movement in lymph nodes can best be described as a random walk with persistence of orientation. The assumption of motility induced cell elongation is consistent with the data. Within the framework of our model the two-photon data suggest that T and B cells are in a single velocity state with large stochastic width. The alternative of three different velocity states with frequent changes of their state and small stochastic width is less likely. Two velocity states can be excluded

Cutting edge: back to "one-way" germinal centers

The present status of germinal center (GC) research is revisited using in silico simulations based on recent lymphocyte motility data in mice. The generally adopted view of several rounds of somatic hypermutations and positive selection is analyzed with special emphasis on the spatial organization of the GC reaction. We claim that the development of dark zones is not necessary for successful GC reactions to develop. We find that a recirculation of positively selected centrocytes to the dark zone is rather unlikely. Instead we propose a scenario that combines a multiple-step mutation and selection concept with a "one-way" GC in the sense of cell migration

Towards whole-organ modelling of tumour growth

Multiscale approaches to modelling biological phenomena are growing rapidly. We present here some recent results on the formulation of a theoretical framework which can be developed into a fully integrative model for cancer growth. The model takes account of vascular adaptation and cell-cycle dynamics. We explore the effects of spatial inhomogeneity induced by the blood flow through the vascular network and of the possible effects of p27 on the cell cycle. We show how the model may be used to investigate the efficiency of drug-delivery protocols

Parameter domains for Turing and stationary flow-distributed waves: I. The influence of nonlinearity

new type of instability in coupled reaction-diffusion-advection systems is analysed in a one-dimensional domain. This instability, arising due to the combined action of flow and diffusion, creates spatially periodic stationary waves termed flow and diffusion-distributed structures (FDS). Here we show, via linear stability analysis, that FDS are predicted in a considerably wider domain and are more robust (in the parameter domain) than the classical Turing instability patterns. FDS also represent a natural extension of the recently discovered flow-distributed oscillations (FDO). Nonlinear bifurcation analysis and numerical simulations in one-dimensional spatial domains show that FDS also have much richer solution behaviour than Turing structures. In the framework presented here Turing structures can be viewed as a particular instance of FDS. We conclude that FDS should be more easily obtainable in chemical systems than Turing (and FDO) structures and that they may play a potentially important role in biological pattern formation

A design principle for vascular beds: the effects of complex blood rheology

We propose a design principle that extends Murray's original optimization principle for vascular architecture to account for complex blood rheology. Minimization of an energy dissipation function enables us to determine how rheology affects the morphology of simple branching networks. The behavior of various physical quantities associated with the networks, such as the wall shear stress and the flow velocity, is also determined. Our results are shown to be qualitatively and quantitatively compatible with independent experimental observations and simulations

Speed of reaction diffusion in embryogenesis

Reaction diffusion systems have been proposed as mechanisms for patterning during many stages of embryonic development. While much attention has been focused on the study of the steady state patterns formed and the robustness of pattern selection, much less is known about the time scales required for pattern formation. Studies of gradient formation by the diffusion of a single morphogen from a localized source have shown that patterning can occur on realistic time scales over distances of a millimeter or less. Reaction diffusion has the potential to give rise to patterns on a faster time scale, since all points in the domain can act as sources of morphogen. However, the speed at which patterning can occur has hitherto not been explored in depth. In this paper, we investigate this issue in specific reaction diffusion models and address the question of whether patterning via reaction diffusion is fast enough to be applicable to morphogenesis

A mathematical model of the effects of hypoxia on the cell-cycle of normal and cancer cells

The evolution of the cell-cycle is known to be influenced by environmental conditions, including lack of extracellular oxygen (hypoxia). Notably, hypoxia appears to have different effects on normal and cancer cells. Whereas both experience hypoxia-induced arrest of the G1 phase of the cell-cycle (i.e. delay in the transition through the restriction point), experimental evidence suggests that only cancer cells undergo hypoxia-induced quiescence (i.e. the transition of the cell to a latent state in which most of the cell functions, including proliferation, are suspended). Here, we extend a model for the cell-cycle due to Tyson and Novak (J. Theor. Biol. 210 (2001) 249) to account for the action of the protein p27. This protein, whose expression is upregulated under hypoxia, inhibits the activation of the cyclin dependent kinases (CDKs), thus preventing DNA synthesis and delaying the normal progression through the cell-cycle. We use a combination of numerical and analytic techniques to study our model. We show that it reproduces many features of the response to hypoxia of normal and cancer cells, as well as generating experimentally testable predictions. For example our model predicts that cancer cells can undergo quiescence by increasing their levels of p27, whereas for normal cells p27 expression decreases when the cellular growth rate increases

Turing instabilities in general systems

We present necessary and sufficient conditions on the stability matrix of a general n(S2)-dimensional reaction-diffusion system which guarantee that its uniform steady state can undergo a Turing bifurcation. The necessary (kinetic) condition, requiring that the system be composed of an unstable (or activator) and a stable (or inhibitor) subsystem, and the sufficient condition of sufficiently rapid inhibitor diffusion relative to the activator subsystem are established in three theorems which form the core of our results. Given the possibility that the unstable (activator) subsystem involves several species (dimensions), we present a classification of the analytically deduced Turing bifurcations into p (1 h p h (n m 1)) different classes. For n = 3 dimensions we illustrate numerically that two types of steady Turing pattern arise in one spatial dimension in a generic reaction-diffusion system. The results confirm the validity of an earlier conjecture [12] and they also characterise the class of so-called strongly stable matrices for which only necessary conditions have been known before [23, 24]. One of the main consequences of the present work is that biological morphogens, which have so far been expected to be single chemical species [1-9], may instead be composed of two or more interacting species forming an unstable subsystem

A cellular automaton model for tumour growth in inhomogeneous environment

Most of the existing mathematical models for tumour growth and tumour-induced angiogenesis neglect blood flow. This is an important factor on which both nutrient and metabolite supply depend. In this paper we aim to address this shortcoming by developing a mathematical model which shows how blood flow and red blood cell heterogeneity influence the growth of systems of normal and cancerous cells. The model is developed in two stages. First we determine the distribution of oxygen in a native vascular network, incorporating into our model features of blood flow and vascular dynamics such as structural adaptation, complex rheology and red blood cell circulation. Once we have calculated the oxygen distribution, we then study the dynamics of a colony of normal and cancerous cells, placed in such a heterogeneous environment. During this second stage, we assume that the vascular network does not evolve and is independent of the dynamics of the surrounding tissue. The cells are considered as elements of a cellular automaton, whose evolution rules are inspired by the different behaviour of normal and cancer cells. Our aim is to show that blood flow and red blood cell heterogeneity play major roles in the development of such colonies, even when the red blood cells are flowing through the vasculature of normal, healthy tissue