Dr Carina Dunlop
Senior Lecturer
Academic and research departments
Department of Mathematics, Centre for Mathematical and Computational Biology.Biography
My research focuses on the mathematical modelling of biological systems. I joined the University of Surrey as a Lecturer in Mathematics in January 2012. Prior to this I spent four years at Heidelberg University as a postdoctoral fellow, and three and a half years as a postdoctoral researcher at Oxford University. I completed my doctoral research in fluid dynamics in 2004 at the University of Oxford, where I read mathematics as an undergraduate.
Areas of specialism
Biophysics;
Mathematical Biology;
Cancer modelling
University roles and responsibilities
- Admissions Tutor
Affiliations and memberships
Society of Mathematical Biology
Member (Developmental Biology subgroup)
European Society for Mathematical and Theoretical Biology
Member
Research
Research interests
I am interested in using mathematical and computational approaches to solve a broad range of problems in developmental biology, tissue morphogenesis and cancer modelling. A particular research focus is on the cell as a physical object, incorporating an understanding of the role of mechanical forces into models of biological processes. I draw on a diverse range of concepts in pursuing this research ranging from population modelling to the theories of fluid dynamics and elasticity theory. Currently ongoing projects include work on tissue self-organization, mechanical regulation of growth and cellular contractility.
Current students:
PhD - Kieran Boniface, Oct. 2019 - present, Mathematical models for tissue growth and development with applications to organoids and tissue eningeering
PhD - Josephine Solowiej-Wedderburn, Oct. 2017 - present, Mechanical models of cell-substrate interactions
PhD- Adam Nasim, Jan. 2018 - present, Mathematical modelling of pharmacological approaches to cancer treatment (collaboration with Dr James Yates, AstraZeneca)
MMath - Robert Dymott, Importance of interstitial fluid flow in tumour growth
Former students
Euan Littlejohns, PhD, University of Surrey (2014 -18), Continuum Elasticity Models for Tissue Growth and Mechanotransduction
Philip Murray, DPhil, University of Oxford (2004-8): From discrete to continuum models of tumour growth, co-supervisors: Prof. P. K. Maini and Dr M.J. Tindall. now Lecturer, Division of Mathematics, University of Dundee
Matthew Johnston, DPhil,University of Oxford (2004-8): Mathematical modelling of cell population dynamics in the colonic crypt, co-supervisors: Prof. P.K. Maini and Prof. S.J. Chapman
Research projects
Open positions
PhD positions available within the group - please do get in touch.
Mechanobiology
The importance of the cell as a physical object and of mechanical interactions in directing cellular behaviours is now clear. We use mathematical modelling to investigate how physical forces internally generated within cells can direct cellular behaviour through mechanotransduction - and in particular on the role of cell-substrate (in vitro) and cell-ECM interactions in this process.
PKPD modelling for cancer treatments
Quantitative systems pharmacology is being increasingly recognised as being of high value in the pharma sector. We here focus on PKPD models of drug effects in tissue within the context of mathematical oncology. (A collaboration with AstraZeneca)
Mathematical models of morphogenesis for tissue engineering
In many context, healthy tissue function is dependent on structure. The question of how tissues develop their shape and size underpins the field of tissue morphogenesis. By considering how growth and active cellular behaviours can generate force and tissue adaption we investigate this process in a range of tissues looking cellular positioning, branching structures, and size control. This has implications for diseases associated with disregulation of morphogenesis and organoid (synthetic tissue) development.
Modelling for biophysical experiment
The range of techniques available for interrogating the biophysics of cells and tissues has exploded over recent years. These new experiments have illuminated the physical nature of cells. We use mathematical modelling of cellular mechanics to inform the interpretation of what these experimental observations could imply for cell cytoskeletal dynamics and force generation. This is particularly important, where cells are seen to interact with each other in multicellular systems so that isolating individual cell behaviours experimentally is difficult.
My teaching
I am lecturing this year on
MATM040 - Mathematical Biology and Physiology
MAT2050 - Inviscid Fluid Dynamics
Supervision
Postgraduate research supervision
Kieran Boniface, PhD, Oct. 2019 - present, Mathematical models for tissue growth and development with applications to organoids and tissue eningeering
Josephine Solowiej-Wedderburn, PhD, Oct. 2017 - present, Mechanical models of cell-substrate interactions
Adam Nasim, PhD, Jan. 2018 - present, Mathematical modelling of pharmacological approaches to cancer treatment (collaboration with Dr James Yates, AstraZeneca)
My publications
Publications
Murray PJ, Walter A, Fletcher AG, Edwards CM, Tindall MJ, Maini PK (2011) Comparing a discrete and continuum model of the intestinal crypt, PHYSICAL BIOLOGY 8 (2) ARTN 026011
van Leeuwen IMM, Edwards CM, Ilyas M, Byrne HM (2007) Towards a multiscale model of colorectal cancer, World Journal of Gastroenterology 13 (9) pp. 1399-1407
Colorectal cancer (CRC) is one of the best characterised cancers, with extensive data documenting the sequential gene mutations that underlie its development. Complementary datasets are also being generated describing changes in protein and RNA expression, tumour biology and clinical outcome. Both the quantity and the variety of information are inexorably increasing and there is now an accompanying need to integrate these highly disparate datasets. In this article we aim to explain why we believe that mathematical modelling represents a natural tool or language with which to integrate these data and, in so doing, to provide insight into CRC. © 2007 The WJG Press. All rights reserved.
Breward C, Dellar PJ, Edwards CM, Kaouri K, Richardson G, Wilson S (2002) Mathematical modelling of pipe-flow and extrusion of composite materials,
Edwards CM, Howison SD, Ockendon H, Ockendon JR (2008) Non-classical shallow water flows, IMA JOURNAL OF APPLIED MATHEMATICS 73 (1) pp. 137-157 Oxford University Press
This paper deals with violent discontinuities in shallow water flows with large Froude number F. On a horizontal base, the paradigm problem is that of the impact of two fluid layers in situations where the flow can be modelled as two smooth regions joined by a singularity in the flow field. Within the framework of shallow water theory, we show that, over a certain time-scale, this discontinuity may be described by a delta shock, which is a weak solution of the underlying conservation laws in which the depth and mass and momentum fluxes have both delta function and step function components. We also make some conjectures about how this model evolves from the traditional model for jet impacts in which a spout is emitted. For flows on a sloping base, we show that for flow with an aspect ratio of O(F?2) on a base with an O(1) or larger slope, the governing equations admit a new type of discontinuous solution that is also modelled as a delta shock. The physical manifestation of this discontinuity is a small ?tube? of fluid bounding the flow. The delta-shock conditions for this flow are derived and solved for a point source on an inclined plane. This latter delta-shock framework also sheds light on the evolution of the layer impact on a horizontal base
Murray PJ, Edwards CM, Tindall MJ, Maini PK (2009) From a discrete to a continuum model of cell dynamics in one dimension, PHYSICAL REVIEW E 80 (3) ARTN 031912
Kanesaki T, Edwards CM, Schwarz US, Grosshans J (2011) Dynamic ordering of nuclei in syncytial embryos: a quantitative analysis of the role of cytoskeletal networks, INTEGRATIVE BIOLOGY 3 pp. 1112-1119 Royal Society of Chemistry
In syncytial embryos nuclei undergo cycles of division and rearrangement within a common cytoplasm. It is presently unclear to what degree and how the nuclear array maintains positional order in the face of rapid cell divisions. Here we establish a quantitative assay, based on image processing, for analysing the dynamics of the nuclear array. By tracking nuclear trajectories in Drosophila melanogaster embryos, we are able to define and evaluate local and time-dependent measures for the level of geometrical order in the array. We find that after division, order is re-established in a biphasic manner, indicating the competition of different ordering processes. Using mutants and drug injections, we show that the order of the nuclear array depends on cytoskeletal networks organised by centrosomes. While both f-actin and microtubules are required for re-establishing order after mitosis, only f-actin is required to maintain the stability of this arrangement. Furthermore, f-actin function relies on myosin-independent non-contractile filaments that suppress individual nuclear mobility, whereas microtubules promote mobility and attract adjacent nuclei. Actin caps are shown to act to prevent nuclear incorporation into adjacent microtubule baskets. Our data demonstrate that two principal ordering mechanisms thus simultaneously contribute: (1) a passive crowding mechanism in which nuclei and actin caps act as spacers and (2) an active self-organisation mechanism based on a microtubule network.
Koke C, Kanesaki T, Grosshans J, Schwarz US, Dunlop CM (2014) A computational model of nuclear self-organisation in syncytial embryos, Journal of Theoretical Biology 359 pp. 92-100
Syncytial embryos develop through cycles of nuclear division and rearrangement within a common cytoplasm. A paradigm example is Drosophila melanogaster in which nuclei form an ordered array in the embryo surface over cell cycles 10-13. This ordering process is assumed to be essential for subsequent cellularisation. Using quantitative tissue analysis, it has previously been shown that the regrowth of actin and microtubule networks after nuclear division generates reordering forces that counteract its disordering effect (Kanesaki et al., 2011). We present here an individual-based computer simulation modelling the nuclear dynamics. In contrast to similar modelling approaches e.g. epithelial monolayers or tumour spheroids, we focus not on the spatial dependence, but rather on the time-dependence of the interaction laws. We show that appropriate phenomenological inter-nuclear force laws reproduce the experimentally observed dynamics provided that the cytoskeletal network regrows sufficiently quickly after mitosis. Then repulsive forces provided by the actin system are necessary and sufficient to regain the observed level of order in the system, after the strong disruption resulting from cytoskeletal network disassembly and spindle formation. We also observe little mixing of nuclei through cell cycles. Our study highlights the importance of the dynamics of cytoskeletal forces during this critical phase of syncytial development and emphasises the need for real-time experimental data at high temporal resolution. © 2014 Elsevier Ltd.
Edwards CM, Chapman SJ (2007) Biomechanical modelling of colorectal crypt budding and fission, BULLETIN OF MATHEMATICAL BIOLOGY 69 (6) pp. 1927-1942 Springer
This paper presents a biomechanical model for the small pits, called crypts, that line the colon. A continuum approach is adopted, with the crypt epithelium modelled as a growing beam attached to the underlying lamina by cell bonds, which generate tension within the layer. These cell attachments are assumed to be viscoelastic thus allowing for cell progression along the crypt. It is shown that any combination of: an increase in net proliferation (i.e. cell production minus apoptosis), an enlargement of the proliferative compartment, an increase in the strength of the cellular attachment to the underlying lamina, or a change in the rate of cell growth or cell bonding may generate buckling of the tissue. These changes can all be generated by an activating mutation of the Wnt cascade, which is generally accepted to be the first genetic change in colorectal cancer, with subsequent deformation, budding, and crypt fission an observed feature of the adenomatous crypt.
Breward C, Dyson R, Edwards CM, Metcalfe P, Please CP, Zyskin M (2004) Modelling of melt on spinning wheels,
Murray PJ, Edwards CM, Tindall MJ, Maini PK (2012) Classifying general nonlinear force laws in cell-based models via the continuum limit, PHYSICAL REVIEW E 85 (2) ARTN 021921
Johnston MD, Edwards CM, Bodmer WF, Maini PK, Chapman SJ (2007) Mathematical modeling of cell population dynamics in the colonic crypt and in colorectal cancer, Proceedings of the National Academy of Sciences of the United States of America 104 (10) pp. 4008-4013
Johnston MD, Maini PK, Chapman SJ, Edwards CM, Bodmer WF (2010) On the proportion of cancer stem cells in a tumour, JOURNAL OF THEORETICAL BIOLOGY 266 (4) pp. 708-711
Johnston MD, Edwards CM, Bodmer WF, Maini PK, Chapman SJ (2007) Examples of mathematical modeling: Tales from the crypt, Cell Cycle 6 (17) pp. 2106-2112
Mathematical modeling is being increasingly recognized within the biomedical sciences as an important tool that can aid the understanding of biological systems. The heavily regulated cell renewal cycle in the colonic crypt provides a good example of how modeling can be used to find out key features of the system kinetics, and help to explain both the breakdown of homeostasis and the initiation of tumorigenesis. We use the cell population model by Johnston et al.5 to illustrate the power of mathematical modeling by considering two key questions about the cell population dynamics in the colonic crypt. We ask: how can a model describe both homeostasis and unregulated growth in tumorigenesis; and to which parameters in the system is the model most sensitive? In order to address these questions, we discuss what type of modeling approach is most appropriate in the crypt. We use the model to argue why tumorigenesis is observed to occur in stages with long lag phases between periods of rapid growth, and we identify the key parameters. ©2007 Landes Bioscience.
Van Leeuwen IMM, Byrne HM, Johnston MD, Edwards CM, Chapman SJ, Bodmer WF, Maini PK (2007) Modelling multiscale aspects of colorectal cancer, AIP Conference Proceedings 971 pp. 3-7
Colorectal cancer (CRC) is responsible for nearly half a million deaths annually world-wide [11]. We present a series of mathematical models describing the dynamics of the intestinal epithelium and the kinetics of the molecular pathway most commonly mutated in CRC, the Wnt signalling network. We also discuss how we are coupling such models to build a multiscale model of normal and aberrant guts. This will enable us to combine disparate experimental and clinical data, to investigate interactions between phenomena taking place at different levels of organisation and, eventually, to test the efficacy of new drugs on the system as a whole. © 2008 American Institute of Physics.
Edwards CM, Schwarz US (2011) Force localization in contracting cell layers, Physical Review Letters 107 (12) American Physical Society
Epithelial cell layers on soft elastic substrates or pillar arrays are commonly used as model systems for investigating the role of force in tissue growth, maintenance and repair. Here we show analytically that the experimentally observed localization of traction forces to the periphery of the cell layers
does not necessarily imply increased local cell activity, but follows naturally from the elastic problem of a finite-sized contractile layer coupled to an elastic foundation. For homogeneous contractility, the force localization is determined by one dimensionless parameter interpolating between linear and
exponential force profiles for the extreme cases of very soft and very stiff substrates, respectively. If contractility is sufficiently increased at the periphery, outward directed displacements can occur at intermediate positions. We also show that anisotropic extracellular stiffness can lead to force
localization in the stiffer direction, as observed experimentally.
does not necessarily imply increased local cell activity, but follows naturally from the elastic problem of a finite-sized contractile layer coupled to an elastic foundation. For homogeneous contractility, the force localization is determined by one dimensionless parameter interpolating between linear and
exponential force profiles for the extreme cases of very soft and very stiff substrates, respectively. If contractility is sufficiently increased at the periphery, outward directed displacements can occur at intermediate positions. We also show that anisotropic extracellular stiffness can lead to force
localization in the stiffer direction, as observed experimentally.
Schwarz US, Dunlop CM (2012) Developmental Biology: A Growing Role for Computer Simulations, CURRENT BIOLOGY 22 (11) pp. R441-R443 CELL PRESS
Keeping cells separated in well-defined domains is essential for development. A new computational-experimental study elucidates the physical mechanisms that establish and maintain the dorsal-ventral compartment boundary in the Drosophila wing disc and demonstrates the increasing value of computer simulations in developmental biology. © 2012 Elsevier Ltd.
Edwards CM, Ovendon N, Rottschäfer V (2003) Some cracking ideas on egg incubation,
Littlejohns Euan, Dunlop Carina (2018) Mechanotransduction mechanisms in growing spherically structured tissues, New Journal of Physics 20 043041 pp. 043041-1-043041-13 IOP Publishing
There is increasing experimental interest in mechanotransduction in
multi-cellular tissues as opposed to single cells. This is driven by a growing awareness of
the importance of physiologically relevant three-dimensional culture and of cell-cell and
cell-gel interactions in directing growth and development. The paradigm biophysical
technique for investigating tissue level mechanobiology in this context is to grow
model tissues in artificial gels with well-defined mechanical properties. These studies
often indicate that the stiµness of the encapsulating gel can significantly alter cellular
behaviours. We demonstrate here potential mechanisms linking tissue growth with
stiµness-mediated mechanotransduction. We show how tissue growth in gel systems
generates points at which there is a significant qualitative change in the cellular stress
and strain experienced. We show analytically how these potential switching points
depend on the mechanical properties of the constraining gel and predict when they
will occur. Significantly, we identify distinct mechanisms that act separately in each of
the stress and strain fields at diµerent times. These observations suggest growth as a
potential physical mechanism coupling gel stiµness with cellular mechanotransduction
in three-dimensional tissues. We additionally show that non-proliferating areas, in
the case that the constraining gel is soft compared with the tissue, will expand and
contract passively as a result of growth. Central compartment size is thus seen to not
be a reliable indicator on its own for growth initiation or active behaviour.
multi-cellular tissues as opposed to single cells. This is driven by a growing awareness of
the importance of physiologically relevant three-dimensional culture and of cell-cell and
cell-gel interactions in directing growth and development. The paradigm biophysical
technique for investigating tissue level mechanobiology in this context is to grow
model tissues in artificial gels with well-defined mechanical properties. These studies
often indicate that the stiµness of the encapsulating gel can significantly alter cellular
behaviours. We demonstrate here potential mechanisms linking tissue growth with
stiµness-mediated mechanotransduction. We show how tissue growth in gel systems
generates points at which there is a significant qualitative change in the cellular stress
and strain experienced. We show analytically how these potential switching points
depend on the mechanical properties of the constraining gel and predict when they
will occur. Significantly, we identify distinct mechanisms that act separately in each of
the stress and strain fields at diµerent times. These observations suggest growth as a
potential physical mechanism coupling gel stiµness with cellular mechanotransduction
in three-dimensional tissues. We additionally show that non-proliferating areas, in
the case that the constraining gel is soft compared with the tissue, will expand and
contract passively as a result of growth. Central compartment size is thus seen to not
be a reliable indicator on its own for growth initiation or active behaviour.
Hardy Kate, Mora Jocelyn M, Dunlop Carina, Carzaniga Raffaella, Franks Stephen, Fenwick Mark A (2018) Nuclear exclusion of SMAD2/3 in granulosa cells is associated with cell proliferation and follicle activation in the mouse ovary., Journal of Cell Science 131 jcs218123 Company of Biologists
Maintenance and activation of the limited supply of primordial follicles in the ovary are
important determinants of reproductive lifespan. Currently, the molecular programme that
maintains the primordial phenotype and the early events associated with follicle activation
are not well defined. Here we have systematically analysed these events using microscopy
and detailed image analysis. Using the immature mouse ovary as a model, we demonstrate
that the onset of granulosa cell (GC) proliferation results in increased packing density on the
oocyte surface and consequent GC cuboidalisation. These events precede oocyte growth
and nuclear translocation of FOXO3a, a transcription factor important in follicle activation.
Immunolabelling of the TGFb signalling mediators and transcription factors, SMAD2/3,
revealed a striking expression pattern specific to GCs of small follicles. SMAD2/3 was
expressed in the nuclei of primordial GCs but was mostly excluded in early growing follicles.
In activated follicles, GC nuclei lacking SMAD2/3 generally expressed Ki67. These findings
suggest that the first phenotypic changes during follicle activation are observed in GCs, and
that TGFb signalling is fundamental for regulating GC arrest and the onset of proliferation.
important determinants of reproductive lifespan. Currently, the molecular programme that
maintains the primordial phenotype and the early events associated with follicle activation
are not well defined. Here we have systematically analysed these events using microscopy
and detailed image analysis. Using the immature mouse ovary as a model, we demonstrate
that the onset of granulosa cell (GC) proliferation results in increased packing density on the
oocyte surface and consequent GC cuboidalisation. These events precede oocyte growth
and nuclear translocation of FOXO3a, a transcription factor important in follicle activation.
Immunolabelling of the TGFb signalling mediators and transcription factors, SMAD2/3,
revealed a striking expression pattern specific to GCs of small follicles. SMAD2/3 was
expressed in the nuclei of primordial GCs but was mostly excluded in early growing follicles.
In activated follicles, GC nuclei lacking SMAD2/3 generally expressed Ki67. These findings
suggest that the first phenotypic changes during follicle activation are observed in GCs, and
that TGFb signalling is fundamental for regulating GC arrest and the onset of proliferation.
Dunlop Carina (2019) Differential cellular contractility as a mechanism for stiffness sensing, New Journal of Physics 21 063005 pp. 1-9 IOP Publishing
The ability of cells to sense and respond to the mechanical properties of
their environments is fundamental to a range of cellular behaviours, with substrate
stiffness increasingly being found to be a key signalling factor. Although active
contractility of the cytoskeleton is clearly necessary for stiffness sensing in cells,
the physical mechanisms connecting contractility with mechanosensing and molecular
conformational change are not well understood. Here we present a contractility-driven
mechanism for linking changes in substrate stiffness with internal conformational
changes. Cellular contractility is often assumed to imply an associated compressive
strain. We show, however, that where the contractility is non-uniform, localized areas
of internal stretch can be generated as stiffer substrates are encountered. This suggests
a physical mechanism for the stretch-activation of mechanotransductive molecules on
stiffer substrates. Importantly, the areas of internal stretch occur deep within the
cell and not near the cellular perimeter, which region is more traditionally associated
with stiffness sensing through e.g. focal adhesions. This supports recent experimental
results on whole-cell mechanically-driven mechanotransduction. Considering cellular
shape we show that aspect ratio acts as an additional control parameter, so that the
onset of positive strain moves to higher stiffness values in elliptical cells.
their environments is fundamental to a range of cellular behaviours, with substrate
stiffness increasingly being found to be a key signalling factor. Although active
contractility of the cytoskeleton is clearly necessary for stiffness sensing in cells,
the physical mechanisms connecting contractility with mechanosensing and molecular
conformational change are not well understood. Here we present a contractility-driven
mechanism for linking changes in substrate stiffness with internal conformational
changes. Cellular contractility is often assumed to imply an associated compressive
strain. We show, however, that where the contractility is non-uniform, localized areas
of internal stretch can be generated as stiffer substrates are encountered. This suggests
a physical mechanism for the stretch-activation of mechanotransductive molecules on
stiffer substrates. Importantly, the areas of internal stretch occur deep within the
cell and not near the cellular perimeter, which region is more traditionally associated
with stiffness sensing through e.g. focal adhesions. This supports recent experimental
results on whole-cell mechanically-driven mechanotransduction. Considering cellular
shape we show that aspect ratio acts as an additional control parameter, so that the
onset of positive strain moves to higher stiffness values in elliptical cells.
We here consider modelling tissue growth and mechanotransduction, utilising a continuum approach based on the theory of elasticity. Our models are valid for both unhealthy and healthy tissues and assume a spherical symmetry is present. In principal, there are a number of tissue types which can be modelled using our framework, however, we choose to focus on three main tissue paradigms; epithelial cysts, ovarian follicles and avascular tumour spheroids, which all share common geometric features of a central core, surrounded by a proliferating rim of cells. In all cases, the tissues are embedded in a constraining outer gel.
We show that growth within the tissue leads to the build-up of internal stress, and we find potential mechanotransductive mechanisms occurring as a result of this growth. These mechanisms are deemed switching points and we show that they can occur both as a result of the internal stress and of the associated strains. We also find that, for systems with a deformable central core, this core can expand and shrink passively as a result of growth within the surrounding material in parameter regimes that we identify.
We also relax the assumption that the inner core is non-growing and consider that both phases of the tissue undergo growth. This leads to a competition of growth and we show that switching still occurs, but is now dependent upon the growth in both regions. This is considered specifically in the case of ovarian follicles, where we further observe that the cuboidalisation of cells can be produced as a consequence of the mechanics within the system. We next consider the addition of the effect of contractility upon the linear model using two techniques of implementation.
The necessity of the use of nonlinear elasticity is then tested, from which we show that for most parameters within the realms of soft biological tissues, the linear approximation to the nonlinear model is of a sufficient likeness to warrant the use of the linear scheme. We find that the ratio of the Young's modulus of the tissue and the surrounding medium is key in determining the effectiveness of the linear model.
These studies consistently highlight the importance of the mechanical properties of the tissue and surrounding extracellular matrix, specifically stressing the significance of the Young's modulus upon tissue growth dynamics.