A hybrid discrete-continuum mathematical model of pattern prediction in the developing retinal vasculature.

Pathological angiogenesis has been extensively explored by the mathematical modelling community over the past few decades, specifically in the contexts of tumour-induced vascularisation and wound healing. However, there have been relatively few attempts to model angiogenesis associated with normal development, despite the availability of animal models with experimentally accessible and highly ordered vascular topologies: for example, growth and development of the vascular plexus layers in the murine retina. The current study aims to address this issue through the development of a hybrid discrete-continuum mathematical model of the developing retinal vasculature in neonatal mice that is closely coupled with an ongoing experimental programme. The model of the functional vasculature is informed by a range of morphological and molecular data obtained over a period of several days, from 6 days prior to birth to approximately 8 days after birth. The spatio-temporal formation of the superficial retinal vascular plexus (RVP) in wild-type mice occurs in a well-defined sequence. Prior to birth, astrocytes migrate from the optic nerve over the surface of the inner retina in response to a chemotactic gradient of PDGF-A, formed at an earlier stage by migrating retinal ganglion cells (RGCs). Astrocytes express a variety of chemotactic and haptotactic proteins, including VEGF and fibronectin (respectively), which subsequently induce endothelial cell sprouting and modulate growth of the RVP. The developing RVP is not an inert structure; however, the vascular bed adapts and remodels in response to a wide variety of metabolic and biomolecular stimuli. The main focus of this investigation is to understand how these interacting cellular, molecular, and metabolic cues regulate RVP growth and formation. In an earlier one-dimensional continuum model of astrocyte and endothelial migration, we showed that the measured frontal velocities of the two cell types could be accurately reproduced by means of a system of five coupled partial differential equations (Aubert et al. in Bull. Math. Biol. 73:2430-2451, 2011). However, this approach was unable to generate spatial information and structural detail for the entire retinal surface. Building upon this earlier work, a more realistic two-dimensional hybrid PDE-discrete model is derived here that tracks the migration of individual astrocytes and endothelial tip cells towards the outer retinal boundary. Blood perfusion is included throughout plexus development and the emergent retinal architectures adapt and remodel in response to various biological factors. The resulting in silico RVP structures are compared with whole-mounted retinal vasculatures at various stages of development, and the agreement is found to be excellent. Having successfully benchmarked the model against wild-type data, the effect of transgenic over-expression of various genes is predicted, based on the ocular-specific expression of VEGF-A during murine development. These results can be used to help inform future experimental investigations of signalling pathways in ocular conditions characterised by aberrant angiogenesis.

Dynamics of angiogenesis during murine retinal development: a coupled in vivo and in silico study.

The manner in which the superficial retinal vascular plexus (RVP) develops in neonatal wild-type mice is relatively well documented and poses an interesting challenge to the mathematical modelling community. Prior to birth, astrocyte sprouting and proliferation begin around the edge of the optic nerve head, and subsequent astrocyte migration in response to a chemotactic gradient of platelet-derived growth factor (PDGF)-A results in the formation of a dense scaffold on the surface of the inner retina. Astrocytes express a variety of chemotactic and haptotactic proteins that subsequently induce endothelial cell sprouting and modulate growth of the RVP. An experimentally informed, two-dimensional hybrid partial differential equation-discrete model is derived to track the outward migration of individual astrocyte and endothelial tip cells in response to the appropriate biochemical cues. Blood perfusion is included throughout the development of the plexus, and the evolving retinal trees are allowed to adapt and remodel by means of several biological stimuli. The resulting wild-type in silico RVP structures are compared with corresponding experimental whole mounts taken at various stages of development, and agreement between the respective vascular morphologies is found to be excellent. Subsequent numerical predictions help elucidate some of the key biological processes underlying retinal development and demonstrate the potential of the virtual retina for the investigation of various vascular-related diseases of the eye.

A continuum mathematical model of the developing murine retinal vasculature.

Angiogenesis, the process of new vessel growth from pre-existing vasculature, is crucial in many biological situations such as wound healing and embryogenesis. Angiogenesis is also a key regulator of pathogenesis in many clinically important disease processes, for instance, solid tumour progression and ocular diseases. Over the past 10-20 years, tumour-induced angiogenesis has received a lot of attention in the mathematical modelling community and there have also been some attempts to model angiogenesis during wound healing. However, there has been little modelling work of vascular growth during normal development. In this paper, we describe an in silico representation of the developing retinal vasculature in the mouse, using continuum mathematical models consisting of systems of partial differential equations. The equations describe the migratory response of cells to growth factor gradients, the evolution of the capillary blood vessel density, and of the growth factor concentration. Our approach is closely coupled to an associated experimental programme to parameterise our model effectively and the simulations provide an excellent correlation with in vivo experimental data. Future work and development of this model will enable us to elucidate the impact of molecular cues upon vasculature development and the implications for eye diseases such as diabetic retinopathy and neonatal retinopathy of prematurity.

Mathematical modeling of tumor-induced angiogenesis.

Angiogenesis, the growth of a network of blood vessels, is a crucial component of solid tumor growth, linking the relatively harmless avascular and the potentially fatal vascular growth phases of the tumor. As a process, angiogenesis is a well-orchestrated sequence of events involving endothelial cell migration and proliferation; degradation of tissue; new capillary vessel formation; loop formation (anastomosis) and, crucially, blood flow through the network. Once there is flow associated with the nascent network, subsequent growth evolves both temporally and spatially in response to the combined effects of angiogenic factors, migratory cues via the extracellular matrix, and perfusion-related hemodynamic forces in a manner that may be described as both adaptive and dynamic. In this article, we first present a review of previous theoretical and computational models of angiogenesis and then indicate how recent developments in flow models are providing insight into antiangiogenic and chemotherapeutic drug treatment of solid tumors.

Mathematical modelling of flow through vascular networks: implications for tumour-induced angiogenesis and chemotherapy strategies.

Angiogenesis, the formation of blood vessels from a pre-existing vasculature, is a process whereby capillary sprouts are formed in response to externally supplied chemical stimuli. The sprouts then grow and develop, driven initially by endothelial cell migration, and organize themselves into a branched, connected network structure. Subsequent cell proliferation near the sprout-tip permits further extension of the capillary and ultimately completes the process. Angiogenesis occurs during embryogenesis, wound healing, arthritis and during the growth of solid tumours. In this paper we initially generate theoretical capillary networks (which are morphologically similar to those networks observed in vivo) using the discrete mathematical model of Anderson and Chaplain. This discrete model describes the formation of a capillary sprout network via endothelial cell migratory and proliferative responses to external chemical stimuli (tumour angiogenic factors, TAF) supplied by a nearby solid tumour, and also the endothelial cell interactions with the extracellular matrix. The main aim of this paper is to extend this work to examine fluid flow through these theoretical network structures. In order to achieve this we make use of flow modelling tools and techniques (specifically, flow through interconnected networks) from the field of petroleum engineering. Having modelled the flow of a basic fluid through our network, we then examine the effects of fluid viscosity, blood vessel size (i.e., diameter of the capillaries), and network structure/geometry, upon: (i) the rate of flow through the network; (ii) the amount of fluid present in the complete network at any one time; and (iii) the amount of fluid reaching the tumour. The incorporation of fluid flow through the generated vascular networks has highlighted issues that may have major implications for the study of nutrient supply to the tumour (blood/oxygen supply) and, more importantly, for the delivery of chemotherapeutic drugs to the tumour. Indeed, there are also implications for the delivery of anti-angiogenesis drugs to the network itself. Results clearly highlight the important roles played by the structure and morphology of the network, which is, in turn, linked to the size and geometry of the nearby tumour. The connectedness of the network, as measured by the number of loops formed in the network (the anastomosis density), is also found to be of primary significance. Moreover, under certain conditions, the results of our flow simulations show that an injected chemotherapy drug may bypass the tumour altogether.