Search

In this paper we use a hybrid multiscale mathematical model that incorporates both individual cell behaviour through the cell-cycle and the effects of the changing microenvironment through oxygen dynamics to study the multiple effects of radiation therapy. The oxygenation status of the cells is considered as one of the important prognostic markers for determining radiation therapy, as hypoxic cells are less radiosensitive. Another factor that critically affects radiation sensitivity is cell-cycle regulation. The effects of radiation therapy are included in the model using a modified linear quadratic model for the radiation damage, incorporating the effects of hypoxia and cell-cycle in determining the cell-cycle phase-specific radiosensitivity. Furthermore, after irradiation, an individual cell's cell-cycle dynamics are intrinsically modified through the activation of pathways responsible for repair mechanisms, often resulting in a delay/arrest in the cell-cycle. The model is then used to study various combinations of multiple doses of cell-cycle dependent chemotherapies and radiation therapy, as radiation may work better by the partial synchronisation of cells in the most radiosensitive phase of the cell-cycle. Moreover, using this multi-scale model, we investigate the optimum sequencing and scheduling of these multi-modality treatments, and the impact of internal and external heterogeneity on the spatio-temporal patterning of the distribution of tumour cells and their response to different treatment schedules. [ABSTRACT FROM AUTHOR]

Pneumococcal conjugate vaccination has proved highly effective in eliminating vaccine-type pneumococcal carriage and disease. However, the potential adverse effects of serotype replacement remain a major concern when implementing routine childhood pneumococcal conjugate vaccination programmes. Applying a concise predictive model, we present a ready-to-use quantitative tool to investigate the implications of serotype replacement on the net effectiveness of vaccination against invasive pneumococcal disease (IPD) and to guide in the selection of optimal vaccine serotype compositions. We utilise pre-vaccination data on pneumococcal carriage and IPD and assume partial or complete elimination of vaccine-type carriage, its replacement by non-vaccine-type carriage, and stable case-to-carrier ratios (probability of IPD per carriage episode). The model predicts that the post-vaccination IPD incidences in Finland for currently available vaccine serotype compositions can eventually decrease among the target age group of children <5 years of age by 75%. However, due to replacement through herd effects, the decrease among the older population is predicted to be much less (20–40%). We introduce a sequential algorithm for the search of optimal serotype compositions and assess the robustness of inferences to uncertainties in data and assumptions about carriage and IPD. The optimal serotype composition depends on the age group of interest and some serotypes may be highly beneficial vaccine types in one age category (e.g. 6B in children), while being disadvantageous in another. The net effectiveness will be improved only if the added serotype has a higher case-to-carrier ratio than the average case-to-carrier ratio of the current non-vaccine types and the degree of improvement in effectiveness depends on the carriage incidence of the serotype. The serotype compositions of currently available pneumococcal vaccines are not optimal and the effectiveness of vaccination in the population at large could be improved by including new serotypes in the vaccine (e.g. 22 and 9N). [ABSTRACT FROM AUTHOR]

Since the discovery of tumour initiating cells (TICs) in solid tumours, studies focussing on their role in cancer initiation and progression have abounded. The biological interrogation of these cells continues to yield volumes of information on their pro-tumourigenic behaviour, but actionable generalised conclusions have been scarce. Further, new information suggesting a dependence of tumour composition and growth on the microenvironment has yet to be studied theoretically. To address this point, we created a hybrid, discrete/continuous computational cellular automaton model of a generalised stem-cell driven tissue with a simple microenvironment. Using the model we explored the phenotypic traits inherent to the tumour initiating cells and the effect of the microenvironment on tissue growth. We identify the regions in phenotype parameter space where TICs are able to cause a disruption in homeostasis, leading to tissue overgrowth and tumour maintenance. As our parameters and model are non-specific, they could apply to any tissue TIC and do not assume specific genetic mutations. Targeting these phenotypic traits could represent a generalizable therapeutic strategy across cancer types. Further, we find that the microenvironmental variable does not strongly affect the outcomes, suggesting a need for direct feedback from the microenvironment onto stem-cell behaviour in future modelling endeavours. [ABSTRACT FROM AUTHOR]

Interactions among multiple infectious agents are increasingly recognized as a fundamental issue in the understanding of key questions in public health regarding pathogen emergence, maintenance, and evolution. The full description of host-multipathogen systems is, however, challenged by the multiplicity of factors affecting the interaction dynamics and the resulting competition that may occur at different scales, from the within-host scale to the spatial structure and mobility of the host population. Here we study the dynamics of two competing pathogens in a structured host population and assess the impact of the mobility pattern of hosts on the pathogen competition. We model the spatial structure of the host population in terms of a metapopulation network and focus on two strains imported locally in the system and having the same transmission potential but different infectious periods. We find different scenarios leading to competitive success of either one of the strain or to the codominance of both strains in the system. The dominance of the strain characterized by the shorter or longer infectious period depends exclusively on the structure of the population and on the the mobility of hosts across patches. The proposed modeling framework allows the integration of other relevant epidemiological, environmental and demographic factors, opening the path to further mathematical and computational studies of the dynamics of multipathogen systems. [ABSTRACT FROM AUTHOR]

The field of medical systems biology aims to advance understanding of molecular mechanisms that drive disease progression and to translate this knowledge into therapies to effectively treat diseases. A challenging task is the investigation of long-term effects of a (pharmacological) treatment, to establish its applicability and to identify potential side effects. We present a new modeling approach, called Analysis of Dynamic Adaptations in Parameter Trajectories (ADAPT), to analyze the long-term effects of a pharmacological intervention. A concept of time-dependent evolution of model parameters is introduced to study the dynamics of molecular adaptations. The progression of these adaptations is predicted by identifying necessary dynamic changes in the model parameters to describe the transition between experimental data obtained during different stages of the treatment. The trajectories provide insight in the affected underlying biological systems and identify the molecular events that should be studied in more detail to unravel the mechanistic basis of treatment outcome. Modulating effects caused by interactions with the proteome and transcriptome levels, which are often less well understood, can be captured by the time-dependent descriptions of the parameters. ADAPT was employed to identify metabolic adaptations induced upon pharmacological activation of the liver X receptor (LXR), a potential drug target to treat or prevent atherosclerosis. The trajectories were investigated to study the cascade of adaptations. This provided a counter-intuitive insight concerning the function of scavenger receptor class B1 (SR-B1), a receptor that facilitates the hepatic uptake of cholesterol. Although activation of LXR promotes cholesterol efflux and -excretion, our computational analysis showed that the hepatic capacity to clear cholesterol was reduced upon prolonged treatment. This prediction was confirmed experimentally by immunoblotting measurements of SR-B1 in hepatic membranes. Next to the identification of potential unwanted side effects, we demonstrate how ADAPT can be used to design new target interventions to prevent these. [ABSTRACT FROM AUTHOR]

Malaria and lymphatic filariasis (LF) continue to cause a considerable public health burden globally and are co-endemic in many regions of sub-Saharan Africa. These infections are transmitted by the same mosquito species which raises important questions about optimal vector control strategies in co-endemic regions, as well as the effect of the presence of each infection on endemicity of the other; there is currently little consensus on the latter. The need for comprehensive modelling studies to address such questions is therefore significant, yet very few have been undertaken to date despite the recognised explanatory power of reliable dynamic mathematical models. Here, we develop a malaria-LF co-infection modelling framework that accounts for two key interactions between these infections, namely the increase in vector mortality as LF mosquito prevalence increases and the antagonistic Th1/Th2 immune response that occurs in co-infected hosts. We consider the crucial interplay between these interactions on the resulting endemic prevalence when introducing each infection in regions where the other is already endemic (e.g. due to regional environmental change), and the associated timescale for such changes, as well as effects on the basic reproduction number R0 of each disease. We also highlight potential perverse effects of vector controls on human infection prevalence in co-endemic regions, noting that understanding such effects is critical in designing optimal integrated control programmes. Hence, as well as highlighting where better data are required to more reliably address such questions, we provide an important framework that will form the basis of future scenario analysis tools used to plan and inform policy decisions on intervention measures in different transmission settings. [ABSTRACT FROM AUTHOR]

Epithelial cells in the colon are arranged in cylindrical structures called crypts in which cellular proliferation and migration are tightly regulated. We hypothesized that the proliferation patterns of cells may determine the stability of crypts as well as the rates of somatic evolution towards colorectal tumorigenesis. Here, we propose a linear process model of colonic epithelial cells that explicitly takes into account the proliferation kinetics of cells as a function of cell position within the crypt. Our results indicate that proliferation kinetics has significant influence on the speed of cell movement, kinetics of mutation propagation, and sensitivity of the system to selective effects of mutated cells. We found that, of all proliferation curves tested, those with mitotic activities concentrated near the stem cell, including the actual proliferation kinetics determined in in vivo labeling experiments, have a greater ability of delaying the rate of mutation accumulation in colonic stem cells compared to hypothetical proliferation curves with mitotic activities focused near the top of the crypt column. Our model can be used to investigate the dynamics of proliferation and mutation accumulation in spatially arranged tissues. [ABSTRACT FROM AUTHOR]

Breast cancer is the most common malignancy in women and is responsible for hundreds of thousands of deaths annually. As with most cancers, it is a heterogeneous disease and different breast cancer subtypes are treated differently. Understanding the difference in prognosis for breast cancer based on its molecular and phenotypic features is one avenue for improving treatment by matching the proper treatment with molecular subtypes of the disease. In this work, we employed a competition-based approach to modeling breast cancer prognosis using large datasets containing genomic and clinical information and an online real-time leaderboard program used to speed feedback to the modeling team and to encourage each modeler to work towards achieving a higher ranked submission. We find that machine learning methods combined with molecular features selected based on expert prior knowledge can improve survival predictions compared to current best-in-class methodologies and that ensemble models trained across multiple user submissions systematically outperform individual models within the ensemble. We also find that model scores are highly consistent across multiple independent evaluations. This study serves as the pilot phase of a much larger competition open to the whole research community, with the goal of understanding general strategies for model optimization using clinical and molecular profiling data and providing an objective, transparent system for assessing prognostic models. [ABSTRACT FROM AUTHOR]