Segregation of co-cultured multicellular systems: review and modeling consideration.

Cell segregation caused by collective cell migration (CCM) is crucial for morphogenesis, functional development of tissue parts, and is an important aspect in other diseases such as cancer and its metastasis process. Efficiency of the cell segregation depends on the interplay between: (1) biochemical processes such as cell signaling and gene expression and (2) physical interactions between cells. Despite extensive research devoted to study the segregation of various co-cultured systems, we still do not understand the role of physical interactions in cell segregation. Cumulative effects of these physical interactions appear in the form of physical parameters such as: (1) tissue surface tension, (2) viscoelasticity caused by CCM, and (3) solid stress accumulated in multicellular systems. These parameters primarily depend on the interplay between the state of cell-cell adhesion contacts and cell contractility. The role of these physical parameters on the segregation efficiency is discussed on model systems such as co-cultured breast cell spheroids consisting of two subpopulations that are in contact. This review study aims to: (1) summarize biological aspects related to cell segregation, mechanical properties of cell collectives, effects along the biointerface between cell subpopulations and (2) describe from a biophysical/mathematical perspective the same biological aspects summarized before. So that overall it can illustrate the complexity of the biological systems that translate into very complex biophysical/mathematical equations. Moreover, by presenting in parallel these two seemingly different parts (biology vs. equations), this review aims to emphasize the need for experiments to estimate the variety of parameters entering the resulting complex biophysical/mathematical models.

The dynamics along the biointerface between the epithelial and cancer mesenchymal cells: Modeling consideration.

Epithelial cancer is the one of most lethal cancer type worldwide. Targeting the early stage of disease would allow dramatic improvements in the survival of cancer patients. The early stage of the disease is related to cancer cell spreading across surrounding healthy epithelium. Consequently, deeper insight into cell dynamics along the biointerface between epithelial and cancer (mesenchymal) cells is necessary in order to control the disease as soon as possible. Cell dynamics along this epithelial-cancer biointerface is the result of the interplay between various biological and physical mechanisms. Despite extensive research devoted to study cancer cell spreading across the epithelium, we still do not understand the physical mechanisms which influences the dynamics along the biointerface. These physical mechanisms are related to the interplay between physical parameters such as: (1) interfacial tension between cancer and epithelial subpopulations, (2) established interfacial tension gradients, (3) the bending rigidity of the biointerface and its impact on the interfacial tension, (4) surface tension of the subpopulations, (5) viscoelasticity caused by collective cell migration, and (6) cell residual stress accumulation. The main goal of this study is to review some of these physical parameters in the context of the epithelial/cancer biointerface elaborated on the model system such as the biointerface between breast epithelial MCF-10A cells and cancer MDA-MB-231 cells and then to incorporate these parameters into a new biophysical model that could describe the dynamics of the biointerface. We conclude by discussing three biophysical scenarios for cell dynamics along the biointerface, which can occur depending on the magnitude of the generated shear stress: a smooth biointerface, a slightly-perturbed biointerface and an intensively-perturbed biointerface in the context of the Kelvin-Helmholtz instability. These scenarios are related to the probability of cancer invasion.

The rearrangement of co-cultured cellular model systems via collective cell migration.

Cancer invasion through the surrounding epithelium and extracellular matrix (ECM) is the one of the main characteristics of cancer progression. While significant effort has been made to predict cancer cells response under various drug therapies, much less attention has been paid to understand the physical interactions between cancer cells and their microenvironment, which are essential for cancer invasion. Considering these physical interactions on various co-cultured in vitro model systems by emphasizing the role of viscoelasticity, the tissue surface tension, solid stress, and their inter-relations is a prerequisite for establishing the main factors that influence cancer cell spread and develop an efficient strategy to suppress it. This review focuses on the role of viscoelasticity caused by collective cell migration (CCM) in the context of mono-cultured and co-cultured cancer systems, and on the modeling approaches aimed at reproducing and understanding these biological systems. In this context, we do not only review previously-published biophysics models for collective cell migration, but also propose new extensions of those models to include solid stress accumulated within the spheroid core region and cell residual stress accumulation caused by CCM.

Mathematical investigation into the role of macrophage heterogeneity on the temporal and spatio-temporal dynamics of non-small cell lung cancers.

Non Small Cell Lung Cancer (NSCLC) is the most common type of lung cancer, and represents the leading cause of cancer-related deaths worldwide. Experimental studies have shown that these solid cancers are heavily infiltrated with macrophages: anti-tumour M1 macrophages, pro-tumour M2 macrophages, and macrophage subtypes sharing both M1 and M2 properties. In this study we aim to investigate qualitatively the role of macrophages with different functional phenotypes (especially those with mixed phenotypes) on cancer dynamics and the success of different immunotherapies for cancer. To this end, we start with two time-evolving mathematical models for cancer-immune interactions that consider: (i) the effect of the two extreme phenotypes, M1 and M2 cells; (ii) the effect of M1 and M2 cells, as well as a macrophage sub-population with a mixed phenotype (throughout this theoretical study we call these cells "M12 cells"). We compare the dynamics of the two models using computational approaches, paying particular attention to the effect of different anti-cancer immunotherapies that focus on macrophages. Since data available for NSCLC and macrophage interactions are incomplete, we perform a global sensitivity analysis to see the influence of input parameters on model outcomes. Finally, we consider extensions of the previous two models to include also the spatial movement of cells, and investigate the role of macrophages with extreme phenotypes and with mixed phenotypes, on the invasion of cancer cells into the surrounding extracellular matrix (ECM). We use numerical simulations to investigate the macrophages phenotypes at the tumour center versus the invasive margin. Again, we examine the impact of immunotherapies for cancer on the spatial dynamics of cancers and immune cells, and observe a shift in the phenotype of macrophages distributed at the tumour center and invasive margin.

A multi-strain epidemic model for COVID-19 with infected and asymptomatic cases: Application to French data.

Many SARS-CoV-2 variants have appeared over the last months, and many more will continue to appear. Understanding the competition between these different variants could help make future predictions on the evolution of epidemics. In this study we use a mathematical model to investigate the impact of three different SARS-CoV-2 variants on the spread of COVID-19 across France, between January-May 2021 (before vaccination was extended to the full population). To this end, we use the data from Geodes (produced by Public Health France) and a particle swarm optimisation algorithm, to estimate the model parameters and further calculate a value for the basic reproduction number R0. Sensitivity and uncertainty analysis is then used to better understand the impact of estimated parameter values on the number of infections leading to both symptomatic and asymptomatic individuals. The results confirmed that, as expected, the alpha, beta and gamma variants are more transmissible than the original viral strain. In addition, the sensitivity results showed that the beta/gamma variants could have lead to a larger number of infections in France (of both symptomatic and asymptomatic people).

Mathematical investigation of innate immune responses to lung cancer: The role of macrophages with mixed phenotypes.

Macrophages' role in the evolution of solid tumours is a well accepted fact, with the M1-like macrophages having an anti-tumour role and the M2-like macrophages having a pro-tumour role. Despite the fact that some clinical studies on lung tumours have emphasised also the presence of macrophages with mixed M1 and M2 phenotypes in addition to macrophages with distinct phenotypes, the majority of studies still use the distinct M1-M2 classification to predict the evolution of tumours and patient survival. In this theoretical study we use a mathematical modelling and computational approach to investigate the role of macrophages with mixed phenotype on growth/control/elimination of lung tumours. We show that tumour control in the presence of M2→M1 re-polarising treatments is mainly the result of macrophages with mixed phenotypes (due to the assumption of short half-life of M1-like macrophages). We also show that the half-life of various macrophage phenotypes (distinct M1 or mixed M1/M2 phenotypes) impacts the outcome of various therapeutic strategies targeting tumour-associated macrophages. All these results suggest the need for a better experimental understanding of the kinetics of macrophages inside solid tumours.

Leadership Through Influence: What Mechanisms Allow Leaders to Steer a Swarm?

Collective migration of cells and animals often relies on a specialised set of "leaders", whose role is to steer a population of naive followers towards some target. We formulate a continuous model to understand the dynamics and structure of such groups, splitting a population into separate follower and leader types with distinct orientation responses. We incorporate leader influence via three principal mechanisms: a bias in the orientation of leaders towards the destination (orientation-bias), a faster movement of leaders when moving towards the target (speed-bias), and leaders making themselves more clear to followers when moving towards the target (conspicuousness-bias). Analysis and numerical computation are used to assess the extent to which the swarm is successfully shepherded towards the target. We find that successful leadership can occur for each of these three mechanisms across a broad region of parameter space, with conspicuousness-bias emerging as the most robust. However, outside this parameter space we also find various forms of unsuccessful leadership. Forms of excessive influence can result in either swarm-splitting, where the leaders break free and followers are left rudderless, or a loss of swarm cohesion that leads to its eventual dispersal. Forms of low influence, on the other hand, can even generate swarms that move away from the target direction. Leadership must therefore be carefully managed to steer the swarm correctly.

Modelling the effects of environmental heterogeneity within the lung on the tuberculosis life-cycle.

Progress in shortening the duration of tuberculosis (TB) treatment is hampered by the lack of a predictive model that accurately reflects the diverse environment within the lung. This is important as TB has been shown to produce distinct localisations to different areas of the lung during different disease stages, with the environmental heterogeneity within the lung of factors such as air ventilation, blood perfusion and oxygen tension believed to contribute to the apical localisation witnessed during the post-primary form of the disease. Building upon our previous model of environmental lung heterogeneity, we present a networked metapopulation model that simulates TB across the whole lung, incorporating these notions of environmental heterogeneity across the whole TB life-cycle to show how different stages of the disease are influenced by different environmental and immunological factors. The alveolar tissue in the lung is divided into distinct patches, with each patch representing a portion of the total tissue and containing environmental attributes that reflect the internal conditions at that location. We include populations of bacteria and immune cells in various states, and events are included which determine how the members of the model interact with each other and the environment. By allowing some of these events to be dependent on environmental attributes, we create a set of heterogeneous dynamics, whereby the location of the tissue within the lung determines the disease pathological events that occur there. Our results show that the environmental heterogeneity within the lung is a plausible driving force behind the apical localisation during post-primary disease. After initial infection, bacterial levels will grow in the initial infection location at the base of the lung until an adaptive immune response is initiated. During this period, bacteria are able to disseminate and create new lesions throughout the lung. During the latent stage, the lesions that are situated towards the apex are the largest in size, and once a post-primary immune-suppressing event occurs, it is the uppermost lesions that reach the highest levels of bacterial proliferation. Our sensitivity analysis also shows that it is the differential in blood perfusion, causing reduced immune activity towards the apex, which has the biggest influence of disease outputs.

Directionality of Macrophages Movement in Tumour Invasion: A Multiscale Moving-Boundary Approach.

Invasion of the surrounding tissue is one of the recognised hallmarks of cancer (Hanahan and Weinberg in Cell 100: 57-70, 2000. https://doi.org/10.1016/S0092-8674(00)81683-9 ), which is accomplished through a complex heterotypic multiscale dynamics involving tissue-scale random and directed movement of the population of both cancer cells and other accompanying cells (including here, the family of tumour-associated macrophages) as well as the emerging cell-scale activity of both the matrix-degrading enzymes and the rearrangement of the cell-scale constituents of the extracellular matrix (ECM) fibres. The involved processes include not only the presence of cell proliferation and cell adhesion (to other cells and to the extracellular matrix), but also the secretion of matrix-degrading enzymes. This is as a result of cancer cells as well as macrophages, which are one of the most abundant types of immune cells in the tumour micro-environment. In large tumours, these tumour-associated macrophages (TAMs) have a tumour-promoting phenotype, contributing to tumour proliferation and spread. In this paper, we extend a previous multiscale moving-boundary mathematical model for cancer invasion, by considering also the multiscale effects of TAMs, with special focus on the influence that their directional movement exerts on the overall tumour progression. Numerical investigation of this new model shows the importance of the interactions between pro-tumour TAMs and the fibrous ECM, highlighting the impact of the fibres on the spatial structure of solid tumour.

Mathematical models of transmission dynamics and vaccine strategies in Hong Kong during the 2017-2018 winter influenza season.

Two mathematical models described by simple ordinary differential equations are developed to investigate the Hong Kong influenza epidemic during 2017-2018 winter, based on overall epidemic dynamics and different influenza subtypes. The first model, describing the overall epidemic dynamics, provides the starting data for the second model which different influenza subtypes, and whose dynamics is further investigated. Weekly data from December 2017 to May 2018 are obtained from the data base of the Centre of Health Protection in Hong Kong, and used to parametrise the models. With the help of these models, we investigate the impact of different vaccination strategies and determine the corresponding critical vaccination coverage for different vaccine efficacies. The results suggest that at least 72% of Hong Kong population should have been vaccinated during 2017-2018 winter to prevent the seasonal epidemic by herd immunity (while data showed that only a maximum of 11.6% of the population were vaccinated). Our results also show that the critical vaccination coverage decreases with increasing vaccine efficacy, and the increase in one influenza subtype vaccine efficacy may lead to an increase in infections caused by a different subtype.

Improving cancer detection through combinations of cancer and immune biomarkers: a modelling approach.

BACKGROUND: Early cancer diagnosis is one of the most important challenges of cancer research, since in many cancers it can lead to cure for patients with early stage diseases. For epithelial ovarian cancer (which is the leading cause of death among gynaecologic malignancies) the classical detection approach is based on measurements of CA-125 biomarker. However, the poor sensitivity and specificity of this biomarker impacts the detection of early-stage cancers. METHODS: Here we use a computational approach to investigate the effect of combining multiple biomarkers for ovarian cancer (e.g., CA-125 and IL-7), to improve early cancer detection. RESULTS: We show that this combined biomarkers approach could lead indeed to earlier cancer detection. However, the immune response (which influences the level of secreted IL-7 biomarker) plays an important role in improving and/or delaying cancer detection. Moreover, the detection level of IL-7 immune biomarker could be in a range that would not allow to distinguish between a healthy state and a cancerous state. In this case, the construction of solution diagrams in the space generated by the IL-7 and CA-125 biomarkers could allow us predict the long-term evolution of cancer biomarkers, thus allowing us to make predictions on cancer detection times. CONCLUSIONS: Combining cancer and immune biomarkers could improve cancer detection times, and any predictions that could be made (at least through the use of CA-125/IL-7 biomarkers) are patient specific.

Non-local Parabolic and Hyperbolic Models for Cell Polarisation in Heterogeneous Cancer Cell Populations.

Tumours consist of heterogeneous populations of cells. The sub-populations can have different features, including cell motility, proliferation and metastatic potential. The interactions between clonal sub-populations are complex, from stable coexistence to dominant behaviours. The cell-cell interactions, i.e. attraction, repulsion and alignment, processes critical in cancer invasion and metastasis, can be influenced by the mutation of cancer cells. In this study, we develop a mathematical model describing cancer cell invasion and movement for two polarised cancer cell populations with different levels of mutation. We consider a system of non-local hyperbolic equations that incorporate cell-cell interactions in the speed and the turning behaviour of cancer cells, and take a formal parabolic limit to transform this model into a non-local parabolic model. We then investigate the possibility of aggregations to form, and perform numerical simulations for both hyperbolic and parabolic models, comparing the patterns obtained for these models.

Modelling and investigation of the CD4+ T cells - Macrophages paradox in melanoma immunotherapies.

It is generally accepted that tumour cells can be eliminated by M1 anti-tumour macrophages and CD8+ T cells. However, experimental results over the past 10-15 years have shown that B16 mouse melanoma cells can be eliminated by the CD4+ T cells alone (either Th1 or Th2 sub-types), in the absence of CD8+ T cells. In some studies, elimination of B16 melanoma was associated with a Th1 immune response (i.e., elimination occurred in the presence of cytokines produced by Th1 cells), while in other studies melanoma elimination was associated with a Th2 immune response (i.e., elimination occurred in the presence of cytokines produced by Th2 cells). Moreover, macrophages have been shown to be present inside the tumours, during both Th1 and Th2 immune responses. To investigate the possible biological mechanisms behind these apparently contradictory results, we develop a class of mathematical models for the dynamics of Th1 and Th2 cells, and M1 and M2 macrophages in the presence/absence of tumour cells. Using this mathematical model, we show that depending on the re-polarisation rates between M1 and M2 macrophages, we obtain tumour elimination in the presence of a type-I immune response (i.e., more Th1 and M1 cells, compared to the Th2 and M2 cells), or in the presence of a type-II immune response (i.e., more Th2 and M2 cells). Moreover, tumour elimination is also possible in the presence of a mixed type-I/type-II immune response. Tumour growth always occurs in the presence of a type-II immune response, as observed experimentally. Finally, tumour dormancy is the result of a delicate balance between the pro-tumour effects of M2 cells and the anti-tumour effects of M1 and Th1 cells.

The re-polarisation of M2 and M1 macrophages and its role on cancer outcomes.

The anti-tumour and pro-tumour roles of Th1/Th2 immune cells and M1/M2 macrophages have been documented by numerous experimental studies. However, it is still unknown how these immune cells interact with each other to control tumour dynamics. Here, we use a mathematical model for the interactions between mouse melanoma cells, Th2/Th1 cells and M2/M1 macrophages, to investigate the unknown role of the re-polarisation between M1 and M2 macrophages on tumour growth. The results show that tumour growth is associated with a type-II immune response described by large numbers of Th2 and M2 cells. Moreover, we show that (i) the ratio k of the transition rates k12 (for the re-polarisation M1→M2) and k21 (for the re-polarisation M2→M1) is important in reducing tumour population, and (ii) the particular values of these transition rates control the delay in tumour growth and the final tumour size. We also perform a sensitivity analysis to investigate the effect of various model parameters on changes in the tumour cell population, and confirm that the ratio k alone and the ratio of M2 and M1 macrophage populations at earlier times (e.g., day 7) cannot always predict the final tumour size.

Mathematical Models for Immunology: Current State of the Art and Future Research Directions.

The advances in genetics and biochemistry that have taken place over the last 10 years led to significant advances in experimental and clinical immunology. In turn, this has led to the development of new mathematical models to investigate qualitatively and quantitatively various open questions in immunology. In this study we present a review of some research areas in mathematical immunology that evolved over the last 10 years. To this end, we take a step-by-step approach in discussing a range of models derived to study the dynamics of both the innate and immune responses at the molecular, cellular and tissue scales. To emphasise the use of mathematics in modelling in this area, we also review some of the mathematical tools used to investigate these models. Finally, we discuss some future trends in both experimental immunology and mathematical immunology for the upcoming years.

Memory versus effector immune responses in oncolytic virotherapies.

The main priority when designing cancer immuno-therapies has been to seek viable biological mechanisms that lead to permanent cancer eradication or cancer control. Understanding the delicate balance between the role of effector and memory cells on eliminating cancer cells remains an elusive problem in immunology. Here we make an initial investigation into this problem with the help of a mathematical model for oncolytic virotherapy; although the model can in fact be made general enough to be applied also to other immunological problems. According to this model, we find that long-term cancer control is associated with a large number of persistent effector cells (irrespective of the initial peak in effector cell numbers). However, this large number of persistent effector cells is sustained by a relatively large number of memory cells. Moreover, the results of the mathematical model suggest that cancer control from a dormant state cannot be predicted by the size of the memory population.

Symmetries and pattern formation in hyperbolic versus parabolic models of self-organised aggregation.

The study of self-organised collective animal behaviour, such as swarms of insects or schools of fish, has become over the last decade a very active research area in mathematical biology. Parabolic and hyperbolic models have been used intensively to describe the formation and movement of various aggregative behaviours. While both types of models can exhibit aggregation-type patterns, studies on hyperbolic models suggest that these models can display a larger variety of spatial and spatio-temporal patterns compared to their parabolic counterparts. Here we use stability, symmetry and bifurcation theory to investigate this observation more rigorously, an approach not attempted before to compare and contrast aggregation patterns in models for collective animal behaviors. To this end, we consider a class of nonlocal hyperbolic models for self-organised aggregations that incorporate various inter-individual communication mechanisms, and take the formal parabolic limit to transform them into nonlocal parabolic models. We then discuss the symmetry of these nonlocal hyperbolic and parabolic models, and the types of bifurcations present or lost when taking the parabolic limit. We show that the parabolic limit leads to a homogenisation of the inter-individual communication, and to a loss of bifurcation dynamics (in particular loss of Hopf bifurcations). This explains the less rich patterns exhibited by the nonlocal parabolic models. However, for multiple interacting populations, by breaking the population interchange symmetry of the model, one can preserve the Hopf bifurcations that lead to the formation of complex spatio-temporal patterns that describe moving aggregations.

Modeling and investigating malaria P. Falciparum and P. Vivax infections: Application to Djibouti data.

Malaria is an infectious and communicable disease, caused by one or more species of Plasmodium parasites. There are five species of parasites responsible for malaria in humans, of which two, Plasmodium Falciparum and Plasmodium Vivax, are the most dangerous. In Djibouti, the two species of Plasmodium are present in different proportions in the infected population: 77% of P. Falciparum and 33% of P. Vivax. In this study we present a new mathematical model describing the temporal dynamics of Plasmodium Falciparum and Plasmodium Vivax co-infection. We focus briefly on the well posedness of this model and on the calculation of the basic reproductive numbers for the infections with each Plasmodium species that help us understand the long-term dynamics of this model (i.e., existence and stability of various eqiuilibria). Then we use computational approaches to: (a) identify model parameters using real data on malaria infections in Djibouti; (b) illustrate the influence of different estimated parameters on the basic reproduction numbers; (c) perform global sensitivity and uncertainty analysis for the impact of various model parameters on the transient dynamics of infectious mosquitoes and infected humans, for infections with each of the Plasmodium species. The originality of this research stems from employing the FAST method and the LHS method to identify the key factors influencing the progression of the disease within the population of Djibouti. In addition, sensitivity analysis identified the most influential parameter for Falciparium and Vivax reproduction rates. Finally, the uncertainty analysis enabled us to understand the variability of certain parameters on the infected compartments.

Simulating photodynamic therapy for the treatment of glioblastoma using Monte Carlo radiative transport.

Significance: Glioblastoma (GBM) is a rare but deadly form of brain tumor with a low median survival rate of 14.6 months, due to its resistance to treatment. An independent simulation of the INtraoperative photoDYnamic therapy for GliOblastoma (INDYGO) trial, a clinical trial aiming to treat the GBM resection cavity with photodynamic therapy (PDT) via a laser coupled balloon device, is demonstrated. Aim: To develop a framework providing increased understanding for the PDT treatment, its parameters, and their impact on the clinical outcome. Approach: We use Monte Carlo radiative transport techniques within a computational brain model containing a GBM to simulate light path and PDT effects. Treatment parameters (laser power, photosensitizer concentration, and irradiation time) are considered, as well as PDT's impact on brain tissue temperature. Results: The simulation suggests that 39% of post-resection GBM cells are killed at the end of treatment when using the standard INDYGO trial protocol (light fluence = 200 J/cm2 at balloon wall) and assuming an initial photosensitizer concentration of 5 µM. Increases in treatment time and light power (light fluence = 400 J/cm2 at balloon wall) result in further cell kill but increase brain cell temperature, which potentially affects treatment safety. Increasing the p hotosensitizer concentration produces the most significant increase in cell kill, with 61% of GBM cells killed when doubling concentration to 10 µM and keeping the treatment time and power the same. According to these simulations, the standard trial protocol is reasonably well optimized with improvements in cell kill difficult to achieve without potentially dangerous increases in temperature. To improve treatment outcome, focus should be placed on improving the photosensitizer. Conclusions: With further development and optimization, the simulation could have potential clinical benefit and be used to help plan and optimize intraoperative PDT treatment for GBM.

Patterns of spread of influenza A in Canada.

Understanding spatial patterns of influenza transmission is important for designing control measures. We investigate spatial patterns of laboratory-confirmed influenza A across Canada from October 1999 to August 2012. A statistical analysis (generalized linear model) of the seasonal epidemics in this time period establishes a clear spatio-temporal pattern, with influenza emerging earlier in western provinces. Early emergence is also correlated with low temperature and low absolute humidity in the autumn. For the richer data from the 2009 pandemic, a mechanistic mathematical analysis, based on a transmission model, shows that both school terms and weather had important effects on pandemic influenza transmission.