Wound healing angiogenesis the clinical implications of a simple mathematical model

Nonhealing wounds are a major burden for health care systems worldwide. In addition, a patient who suffers from this type of wound usually has a reduced quality of life. While the wound healing process is undoubtedly complex, in this paper we develop a deterministic mathematical model, formulated as a system of partial differential equations, that focusses on an important aspect of successful healing: oxygen supply to the wound bed by a combination of diffusion from the surrounding unwounded tissue and delivery from newly-formed blood vessels. While the model equations can be solved numerically, the emphasis here is on the use of asymptotic methods to establish conditions under which new blood vessel growth can be initiated and wound-bed angiogenesis can progress. These conditions are given in terms of key model parameters including the rate of oxygen supply and its rate of consumption in the wound. We use our model to discuss the clinical use of treatments such as hyperbaric oxygen therapy, wound bed debridement, and revascularisation therapy that have the potential to initiate healing in chronic, stalled wounds

The Two Regime method for optimizing stochastic reaction-diffusion simulations

The computer simulation of stochastic reaction-diffusion processes in biology is often done using either compartment-based (spatially discretized) simulations or molecular-based (Brownian dynamics) approaches. Compartment-based approaches can yield quick and accurate mesoscopic results but lack the level of detail that is characteristic of the more computationally intensive molecular-based models. Often microscopic detail is only required in a small region but currently the best way to achieve this detail is to use a resource intensive model over the whole domain. We introduce the Two Regime Method (TRM) in which a molecular-based algorithm is used in part of the computational domain and a compartment-based approach is used elsewhere in the computational domain. We apply the TRM to two test problems including a model from developmental biology. We thereby show that the TRM is accurate and subsequently may be used to inspect both mesoscopic and microscopic detail of reaction diffusion simulations according to the demands of the modeller

Multiscale stochastic reaction-diffusion modelling: application to actin dynamics in filopodia

Two multiscale (hybrid) stochastic reaction-diffusion models of actin dynamics in a filopodium are investigated. Both hybrid algorithms combine compartment-based and molecular-based stochastic reaction-diffusion models. The first hybrid model is based on the models previously\ud developed in the literature. The second hybrid model is based on the application of recently developed two-regime method (TRM) to a fully molecular-based model which is also developed in this paper. The results of hybrid models are compared with the results of the molecular-based model. It is shown that both approaches give comparable results, although the TRM model better agrees quantitatively with the molecular-based model

Diffusive spatio-temporal noise in a first-passage time model for intracellular calcium release

The intracellular release of calcium from the endoplasmic reticulum is controlled by ion channels. The resulting calcium signals exhibit a rich spatio-temporal signature, which originates at least partly from microscopic fluctuations. While stochasticity in the gating transition of ion channels has been incorporated into many models, the distribution of calcium is usually described by deterministic reaction-diffusion equations. Here we test the validity of the latter modeling approach by using two different models to calculate the frequency of localized calcium signals (calcium puffs) from clustered IP3 receptor channels. The complexity of the full calcium system is here limited to the basic opening mechanism of the ion channels and, in the mathematical reduction simplifies to the calculation of a first passage time. Two models are then studied: (i) a hybrid model, where channel gating is treated stochastically, while calcium concentration is deterministic and (ii) a fully stochastic model with noisy channel gating and Brownian calcium ion motion. The second model utilises the recently developed two-regime method [M. B. Flegg, S. J. Chapman, and R. Erban, “The two-regime method for optimizing stochastic reaction-diffusion simulations,” J. R. Soc., Interface9, 859–868 (Year: 2012)]10.1098/rsif.2011.0574 in order to simulate a large domain with precision required only near the Ca2+ absorbing channels. The expected time for a first channel opening that results in a calcium puff event is calculated. It is found that for a large diffusion constant, predictions of the interpuff time are significantly overestimated using the model (i) with a deterministic non-spatial calcium variable. It is thus demonstrated that the presence of diffusive noise in local concentrations of intracellular Ca2+ ions can substantially influence the occurrence of calcium signals. The presented approach and results may also be relevant for other cell-physiological first-passage time problems with small ligand concentration and high cooperativity

Multiscale reaction-diffusiion algorithms: PDE-assisted Brownian dynamics

Two algorithms that combine Brownian dynamics (BD) simulations with mean-field partial differential equations (PDEs) are presented. This PDE-assisted Brownian dynamics (PBD) methodology provides exact particle tracking data in parts of the domain, whilst making use of a mean-field reaction-diffusion PDE description elsewhere. The first PBD algorithm couples BD simulations with PDEs by randomly creating new particles close to the interface, which partitions the domain, and by reincorporating particles into the continuum PDE-description when they cross the interface. The second PBD algorithm introduces an overlap region, where both descriptions exist in parallel. It is shown that the overlap region is required to accurately compute variances using PBD simulations. Advantages of both PBD approaches are discussed and illustrative numerical examples are presented

Surgical outcomes of borderline breast lesions detected by needle biopsy in a breast screening program

<p>Abstract</p> <p>Background</p> <p>The Australian Capital Territory and South East New South Wales branch of BreastScreen Australia (BreastScreen ACT&SENSW) performs over 20,000 screening mammograms annually. This study describes the outcome of surgical biopsies of the breast performed as a result of a borderline lesion being identified after screening mammography and subsequent workup.</p> <p>A secondary aim was to identify any parameters, such as a family history of breast cancer, or radiological findings that may indicate which borderline lesions are likely to be upgraded to malignancy after surgery.</p> <p>Methods</p> <p>From a period of just over eight years, all patients of BreastScreen ACT&SENSW who were diagnosed with a borderline breast lesion were identified. These women had undergone needle biopsy in Breastscreen ACT&SENSW and either atypical ductal hyperplasia (ADH), flat epithelial atypia (FEA), atypical lobular hyperplasia (ALH), radial scar/complex sclerosing lesion, papillary lesion, mucocoele-like lesion (MLL) or lobular carcinoma in situ (LCIS) was found. Final outcomes for each type of borderline lesion after referral for surgical biopsy were recorded and analysed. Results of the surgical biopsy were compared to the type of needle biopsy and its result, radiological findings and family history status.</p> <p>Results</p> <p>Of the 94 surgical biopsies performed due to the presence of a borderline breast lesion, 20% showed benign pathology, 55% remained as borderline lesions, 17% showed non-invasive malignancy and 7% showed invasive malignancy. VALCS biopsy was the most common needle biopsy method used to identify the lesions in this study (76%). Malignant outcomes resulted from 24% of the surgical biopsies, with the most common malignant lesion being non-comedo ductal carcinoma <it>in situ </it>(DCIS). The most common borderline lesion for which women underwent surgical biopsy was ADH (38%). Of these women, 22% were confirmed as ADH on surgical biopsy and 47% with a malignancy.</p> <p>Conclusions</p> <p>Further research is required to determine whether characteristics of the mammographic lesion (particularly calcification patterns), the area targeted for biopsy and number of core samples retrieved, can indicate a closer correlation with eventual pathology. This study identified no findings in the diagnostic assessment that could exclude women with borderline lesions from surgical biopsy.</p

Optimal Interruption of P. vivax Malaria Transmission Using Mass Drug Administration

Plasmodium vivax is the most geographically widespread malaria-causing parasite resulting in significant associated global morbidity and mortality. One of the factors driving this widespread phenomenon is the ability of the parasites to remain dormant in the liver. Known as ‘hypnozoites’, they reside in the liver following an initial exposure, before activating later to cause further infections, referred to as ‘relapses’. As around 79–96% of infections are attributed to relapses from activating hypnozoites, we expect it will be highly impactful to apply treatment to target the hypnozoite reservoir (i.e. the collection of dormant parasites) to eliminate P. vivax. Treatment with radical cure, for example tafenoquine or primaquine, to target the hypnozoite reservoir is a potential tool to control and/or eliminate P. vivax. We have developed a deterministic multiscale mathematical model as a system of integro-differential equations that captures the complex dynamics of P. vivax hypnozoites and the effect of hypnozoite relapse on disease transmission. Here, we use our multiscale model to study the anticipated effect of radical cure treatment administered via a mass drug administration (MDA) program. We implement multiple rounds of MDA with a fixed interval between rounds, starting from different steady-state disease prevalences. We then construct an optimisation model with three different objective functions motivated on a public health basis to obtain the optimal MDA interval. We also incorporate mosquito seasonality in our model to study its effect on the optimal treatment regime. We find that the effect of MDA interventions is temporary and depends on the pre-intervention disease prevalence (and choice of model parameters) as well as the number of MDA rounds under consideration. The optimal interval between MDA rounds also depends on the objective (combinations of expected intervention outcomes). We find radical cure alone may not be enough to lead to P. vivax elimination under our mathematical model (and choice of model parameters) since the prevalence of infection eventually returns to pre-MDA levels

Modeling the growth of multicellular cancer spheroids in a\ud bioengineered 3D microenvironment and their treatment with an\ud anti-cancer drug

A critical step in the dissemination of ovarian cancer cells is the formation of multicellular spheroids from cells shed from the primary tumor. The objectives of this study were to establish and validate bioengineered three-dimensional (3D) microenvironments for culturing ovarian cancer cells in vitro and simultaneously to develop computational models describing the growth of multicellular spheroids in these bioengineered matrices. Cancer cells derived from human epithelial ovarian carcinoma were embedded within biomimetic hydrogels of varying stiffness and cultured for up to 4 weeks. Immunohistochemistry was used to quantify the dependence of cell proliferation and apoptosis on matrix stiffness, long-term culture and treatment with the anti-cancer drug paclitaxel.\ud \ud Two computational models were developed. In the first model, each spheroid was treated as an incompressible porous medium, whereas in the second model the concept of morphoelasticity was used to incorporate details about internal stresses and strains. Each model was formulated as a free boundary problem. Functional forms for cell proliferation and apoptosis motivated by the experimental work were applied and the predictions of both models compared with the output from the experiments. Both models simulated how the growth of cancer spheroids was influenced by mechanical and biochemical stimuli including matrix stiffness, culture time and treatment with paclitaxel. Our mathematical models provide new perspectives on previous experimental results and have informed the design of new 3D studies of multicellular cancer spheroids

Growth of confined cancer spheroids: a combined experimental and mathematical modelling approach

We have integrated a bioengineered three-dimensional platform by generating multicellular cancer spheroids in a controlled microenvironment with a mathematical model to investigate\ud confined tumour growth and to model its impact on cellular processes