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

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

Standardizing the measurement of parasite clearance in falciparum malaria: the parasite clearance estimator

<p>Abstract</p> <p>Background</p> <p>A significant reduction in parasite clearance rates following artesunate treatment of falciparum malaria, and increased failure rates following artemisinin combination treatments (ACT), signaled emergent artemisinin resistance in Western Cambodia. Accurate measurement of parasite clearance is therefore essential to assess the spread of artemisinin resistance in <it>Plasmodium falciparum</it>. The slope of the log-parasitaemia <it>versus </it>time relationship is considered to be the most robust measure of anti-malarial effect. However, an initial lag phase of numerical instability often precedes a steady exponential decline in the parasite count after the start of anti-malarial treatment. This lag complicates the clearance estimation, introduces observer subjectivity, and may influence the accuracy and consistency of reported results.</p> <p>Methods</p> <p>To address this problem, a new approach to modelling clearance of malaria parasites from parasitaemia-time profiles has been explored and validated. The methodology detects when a lag phase is present, selects the most appropriate model (linear, quadratic or cubic) to fit log-transformed parasite data, and calculates estimates of parasite clearance adjusted for this lag phase. Departing from previous approaches, parasite counts below the level of detection are accounted for and not excluded from the calculation.</p> <p>Results</p> <p>Data from large clinical studies with frequent parasite counts were examined. The effect of a lag phase on parasite clearance rate estimates is discussed, using individual patient data examples. As part of the World Wide Antimalarial Resistance Network's (WWARN) efforts to make innovative approaches available to the malaria community, an automated informatics tool: the parasite clearance estimator has been developed.</p> <p>Conclusions</p> <p>The parasite clearance estimator provides a consistent, reliable and accurate method to estimate the lag phase and malaria parasite clearance rate. It could be used to detect early signs of emerging resistance to artemisinin derivatives and other compounds which affect ring-stage clearance.</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

On the mathematical modeling of wound healing angiogenesis in skin as a reaction-transport process

Over the last 30 years, numerous research groups have attempted to provide mathematical descriptions of the skin wound healing process. The development of theoretical models of the interlinked processes that underlie the healing mechanism has yielded considerable insight into aspects of this critical phenomenon that remain difficult to investigate empirically. In particular, the mathematical modeling of angiogenesis, i.e., capillary sprout growth, has offered new paradigms for the understanding of this highly complex and crucial step in the healing pathway. With the recent advances in imaging and cell tracking, the time is now ripe for an appraisal of the utility and importance of mathematical modeling in wound healing angiogenesis research. The purpose of this review is to pedagogically elucidate the conceptual principles that have underpinned the development of mathematical descriptions of wound healing angiogenesis, specifically those that have utilized a continuum reaction-transport framework, and highlight the contribution that such models have made toward the advancement of research in this field. We aim to draw attention to the common assumptions made when developing models of this nature, thereby bringing into focus the advantages and limitations of this approach. A deeper integration of mathematical modeling techniques into the practice of wound healing angiogenesis research promises new perspectives for advancing our knowledge in this area. To this end we detail several open problems related to the understanding of wound healing angiogenesis, and outline how these issues could be addressed through closer cross-disciplinary collaboration

Crystalline Order on a Sphere and the Generalized Thomson Problem

We attack generalized Thomson problems with a continuum formalism which exploits a universal long range interaction between defects depending on the Young modulus of the underlying lattice. Our predictions for the ground state energy agree with simulations of long range power law interactions of the form 1/r^{gamma} (0 < gamma < 2) to four significant digits. The regime of grain boundaries is studied in the context of tilted crystalline order and the generality of our approach is illustrated with new results for square tilings on the sphere.Comment: 4 pages, 5 eps figures Fig. 2 revised, improved Fig. 3, reference typo fixe

Breast milk-derived human milk oligosaccharides promote Bifidobacterium interactions within a single ecosystem

Diet-microbe interactions play an important role in modulating the early-life microbiota, with Bifidobacterium strains and species dominating the gut of breast-fed infants. Here, we sought to explore how infant diet drives distinct bifidobacterial community composition and dynamics within individual infant ecosystems. Genomic characterisation of 19 strains isolated from breast-fed infants revealed a diverse genomic architecture enriched in carbohydrate metabolism genes, which was distinct to each strain, but collectively formed a pangenome across infants. Presence of gene clusters implicated in digestion of human milk oligosaccharides (HMOs) varied between species, with growth studies indicating that within single infants there were differences in the ability to utilise 2'FL and LNnT HMOs between strains. Cross-feeding experiments were performed with HMO degraders and non-HMO users (using spent or 'conditioned' media and direct co-culture). Further H-NMR analysis identified fucose, galactose, acetate, and N-acetylglucosamine as key by-products of HMO metabolism; as demonstrated by modest growth of non-HMO users on spend media from HMO metabolism. These experiments indicate how HMO metabolism permits the sharing of resources to maximise nutrient consumption from the diet and highlights the cooperative nature of bifidobacterial strains and their role as 'foundation' species in the infant ecosystem. The intra- and inter-infant bifidobacterial community behaviour may contribute to the diversity and dominance of Bifidobacterium in early life and suggests avenues for future development of new diet and microbiota-based therapies to promote infant health

Spatio-temporal spread of artemisinin resistance in Southeast Asia

Current malaria elimination targets must withstand a colossal challenge–resistance to the current gold standard antimalarial drug, namely artemisinin derivatives. If artemisinin resistance significantly expands to Africa or India, cases and malaria-related deaths are set to increase substantially. Spatial information on the changing levels of artemisinin resistance in Southeast Asia is therefore critical for health organisations to prioritise malaria control measures, but available data on artemisinin resistance are sparse. We use a comprehensive database from the WorldWide Antimalarial Resistance Network on the prevalence of non-synonymous mutations in the Kelch 13 (K13) gene, which are known to be associated with artemisinin resistance, and a Bayesian geostatistical model to produce spatio-temporal predictions of artemisinin resistance. Our maps of estimated prevalence show an expansion of the K13 mutation across the Greater Mekong Subregion from 2000 to 2022. Moreover, the period between 2010 and 2015 demonstrated the most spatial change across the region. Our model and maps provide important insights into the spatial and temporal trends of artemisinin resistance in a way that is not possible using data alone, thereby enabling improved spatial decision support systems on an unprecedented fine-scale spatial resolution. By predicting for the first time spatio-temporal patterns and extents of artemisinin resistance at the subcontinent level, this study provides critical information for supporting malaria elimination goals in Southeast Asia

A model for malaria treatment evaluation in the presence of multiple species

Plasmodium (P.) falciparum and P. vivax are the two most common causes of malaria. While the majority of deaths and severe morbidity are due to P. falciparum, P. vivax poses a greater challenge to eliminating malaria outside of Africa due to its ability to form latent liver stage parasites (hypnozoites), which can cause relapsing episodes within an individual patient. In areas where P. falciparum and P. vivax are co-endemic, individuals can carry parasites of both species simultaneously. These mixed infections complicate dynamics in several ways; treatment of mixed infections will simultaneously affect both species, P. falciparum can mask the detection of P. vivax, and it has been hypothesised that clearing P. falciparum may trigger a relapse of dormant P. vivax. When mixed infections are treated for only blood-stage parasites, patients are at risk of relapse infections due to P. vivax hypnozoites. We present a stochastic mathematical model that captures interactions between P. falciparum and P. vivax, and incorporates both standard schizontocidal treatment (which targets blood-stage parasites) and radical treatment (which additionally targets liver-stage parasites). We apply this model to assess the implications of different treatment coverage of radical cure for mixed and P. vivax infections and a so-called "unified radical cure" treatment strategy for P. falciparum, P. vivax and mixed infections. We find that a unified radical cure strategy, with G6PD screening, leads to a substantially lower incidence of malaria cases and deaths overall. We perform a one-way sensitivity analysis to highlight important model parameters