Mathematical modelling of extracellular matrix dynamics using discrete cells: fiber orientation and tissue regeneration.

Matrix orientation plays a crucial role in determining the severity of scar tissue after dermal wounding. We present a model framework which allows us to examine the interaction of many of the factors involved in orientation and alignment. Within this framework, cells are considered as discrete objects, while the matrix is modelled as a continuum. Using numerical simulations, we investigate the effect on alignment of changing cell properties and of varying cell interactions with collagen and fibrin.     

Enzyme kinetics for a two-step enzymic reaction with comparable initial enzyme-substrate ratios

We extend the validity of the quasi-steady state assumption for a model double intermediate enzyme-substrate reaction to include the case where the ratio of initial enzyme to substrate concentration is not necessarily small. Simple analytical solutions are obtained when the reaction rates and the initial substrate concentration satisfy a certain condition. These analytical solutions compare favourably with numerical solutions of the full system of differential equations describing the reaction. Experimental methods are suggested which might permit the application of the quasi-steady state assumption to reactions where it may not have been obviously applicable before.     

A multiscale model of complex endothelial cell dynamics in early angiogenesis

We introduce a hybrid two-dimensional multiscale model of angiogenesis, the process by which endothelial cells (ECs) migrate from a pre-existing vascular bed in response to local environmental cues and cell-cell interactions, to create a new vascular network. Recent experimental studies have highlighted a central role of cell rearrangements in the formation of angiogenic networks. Our model accounts for this phenomenon via the heterogeneous response of ECs to their microenvironment. These cell rearrangements, in turn, dynamically remodel the local environment. The model reproduces characteristic features of angiogenic sprouting that include branching, chemotactic sensitivity, the brush border effect, and cell mixing. These properties, rather than being hardwired into the model, emerge naturally from the gene expression patterns of individual cells. After calibrating and validating our model against experimental data, we use it to predict how the structure of the vascular network changes as the baseline gene expression levels of the VEGF-Delta-Notch pathway, and the composition of the extracellular environment, vary. In order to investigate the impact of cell rearrangements on the vascular network structure, we introduce the mixing measure, a scalar metric that quantifies cell mixing as the vascular network grows. We calculate the mixing measure for the simulated vascular networks generated by ECs of different lineages (wild type cells and mutant cells with impaired expression of a specific receptor). Our results show that the time evolution of the mixing measure is directly correlated to the generic features of the vascular branching pattern, thus, supporting the hypothesis that cell rearrangements play an essential role in sprouting angiogenesis. Furthermore, we predict that lower cell rearrangement leads to an imbalance between branching and sprout elongation. Since the computation of this statistic requires only individual cell trajectories, it can be computed for networks generated in biological experiments, making it a potential biomarker for pathological angiogenesis.     

Pattern formation in generalized Turing systems

Turing's model of pattern formation has been extensively studied analytically and numerically, and there is recent experimental evidence that it may apply in certain chemical systems. The model is based on the assumption that all reacting species obey the same type of boundary condition pointwise on the boundary. We call these scalar boundary conditions. Here we study mixed or nonscalar boundary conditions, under which different species satisfy different boundary conditions at any point on the boundary, and show that qualitatively new phenomena arise in this case. For example, we show that there may be multiple solutions at arbitrarily small lengths under mixed boundary conditions, whereas the solution is unique under homogeneous scalar boundary conditions. Moreover, even when the same solution exists under scalar and mixed boundary conditions, its stability may be different in the two cases. We also show that mixed boundary conditions can reduce the sensitivity of patterns to domain changes.Supported in part by NIH Grant # GM29123     

Analysis of a risk-based model for the growth of AIDS infection.

Several models for the spread of AIDS within a homosexual community have been proposed that incorporate biased mixing of different risk groups. A simple model is presented that captures many of the features of these more complex models. Analytical expressions are derived for the time to the state of maximum infection (SMI) in a particular risk group, the proportion infected at SMI, and the number of infected individuals as the group approaches SMI. These results agree qualitatively with numerical simulations of the model.