Characterizing Chemotherapy-Induced Neutropenia and Monocytopenia Through Mathematical Modelling

Dynamically adjusted cell fate decisions and resilience to mutant invasion during steady-state hematopoiesis revealed by an experimentally parameterized mathematical model

A major next step in hematopoietic stem cell (HSC) biology is to enhance our quantitative understanding of cellular and evolutionary dynamics involved in undisturbed hematopoiesis. Mathematical models have been and continue to be key in this respect, and are most powerful when parameterized experimentally and containing sufficient biological complexity. In this paper, we use data from label propagation experiments in mice to parameterize a mathematical model of hematopoiesis that includes homeostatic control mechanisms as well as clonal evolution. We find that nonlinear feedback control can drastically change the interpretation of kinetic estimates at homeostasis. This suggests that short-term HSC and multipotent progenitors can dynamically adjust to sustain themselves temporarily in the absence of long-term HSCs, even if they differentiate more often than they self-renew in undisturbed homeostasis. Additionally, the presence of feedback control in the model renders the system resilient against mutant invasion. Invasion barriers, however, can be overcome by a combination of age-related changes in stem cell differentiation and evolutionary niche construction dynamics based on a mutant-associated inflammatory environment. This helps us understand the evolution of e.g., TET2 or DNMT3A mutants, and how to potentially reduce mutant burden.

A multiscale model of the action of a capsid assembly modulator for the treatment of chronic hepatitis B

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound occurs rapidly following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs). Here, we develop a multiscale mathematical model of CAM treatment in chronic HBV infection. By fitting the model to participant data from a phase I trial of the first-generation CAM vebicorvir, we estimate the drug's dose-dependent effectiveness and identify the physiological mechanisms that drive the observed biphasic decline in HBV DNA and RNA, and mechanistic differences between HBeAg-positive and negative infection. Finally, we demonstrate analytically and numerically that HBV RNA is more sensitive than HBV DNA to increases in CAM effectiveness.

A Continuation Technique for Maximum Likelihood Estimators in Biological Models

Estimating model parameters is a crucial step in mathematical modelling and typically involves minimizing the disagreement between model predictions and experimental data. This calibration data can change throughout a study, particularly if modelling is performed simultaneously with the calibration experiments, or during an on-going public health crisis as in the case of the COVID-19 pandemic. Consequently, the optimal parameter set, or maximal likelihood estimator (MLE), is a function of the experimental data set. Here, we develop a numerical technique to predict the evolution of the MLE as a function of the experimental data. We show that, when considering perturbations from an initial data set, our approach is significantly more computationally efficient that re-fitting model parameters while producing acceptable model fits to the updated data. We use the continuation technique to develop an explicit functional relationship between fit model parameters and experimental data that can be used to measure the sensitivity of the MLE to experimental data. We then leverage this technique to select between model fits with similar information criteria, a priori determine the experimental measurements to which the MLE is most sensitive, and suggest additional experiment measurements that can resolve parameter uncertainty.

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

COVID-19 virtual patient cohort suggests immune mechanisms driving disease outcomes

To understand the diversity of immune responses to SARS-CoV-2 and distinguish features that predispose individuals to severe COVID-19, we developed a mechanistic, within-host mathematical model and virtual patient cohort. Our results suggest that virtual patients with low production rates of infected cell derived IFN subsequently experienced highly inflammatory disease phenotypes, compared to those with early and robust IFN responses. In these in silico patients, the maximum concentration of IL-6 was also a major predictor of CD8+ T cell depletion. Our analyses predicted that individuals with severe COVID-19 also have accelerated monocyte-to-macrophage differentiation mediated by increased IL-6 and reduced type I IFN signalling. Together, these findings suggest biomarkers driving the development of severe COVID-19 and support early interventions aimed at reducing inflammation.

The Intestinal Microbiota: Impacts of Antibiotics Therapy, Colonization Resistance, and Diseases

Trillions of microbes exist in the human body, particularly the gastrointestinal tract, coevolved with the host in a mutually beneficial relationship. The main role of the intestinal microbiome is the fermentation of non-digestible substrates and increased growth of beneficial microbes that produce key antimicrobial metabolites such as short-chain fatty acids, etc., to inhibit the growth of pathogenic microbes besides other functions. Intestinal microbiota can prevent pathogen colonization through the mechanism of colonization resistance. A wide range of resistomes are present in both beneficial and pathogenic microbes. Giving antibiotic exposure to the intestinal microbiome (both beneficial and hostile) can trigger a resistome response, affecting colonization resistance. The following review provides a mechanistic overview of the intestinal microbiome and the impacts of antibiotic therapy on pathogen colonization and diseases. Further, we also discuss the epidemiology of immunocompromised patients who are at high risk for nosocomial infections, colonization and decolonization of multi-drug resistant organisms in the intestine, and the direct and indirect mechanisms that govern colonization resistance to the pathogens.

COVID-19 virtual patient cohort reveals immune mechanisms driving disease outcomes

To understand the diversity of immune responses to SARS-CoV-2 and distinguish features that predispose individuals to severe COVID-19, we developed a mechanistic, within-host mathematical model and virtual patient cohort. Our results indicate that virtual patients with low production rates of infected cell derived IFN subsequently experienced highly inflammatory disease phenotypes, compared to those with early and robust IFN responses. In these in silico patients, the maximum concentration of IL-6 was also a major predictor of CD8 + T cell depletion. Our analyses predicted that individuals with severe COVID-19 also have accelerated monocyte-to-macrophage differentiation that was mediated by increased IL-6 and reduced type I IFN signalling. Together, these findings identify biomarkers driving the development of severe COVID-19 and support early interventions aimed at reducing inflammation.

AUTHOR SUMMARY

Understanding of how the immune system responds to SARS-CoV-2 infections is critical for improving diagnostic and treatment approaches. Identifying which immune mechanisms lead to divergent outcomes can be clinically difficult, and experimental models and longitudinal data are only beginning to emerge. In response, we developed a mechanistic, mathematical and computational model of the immunopathology of COVID-19 calibrated to and validated against a broad set of experimental and clinical immunological data. To study the drivers of severe COVID-19, we used our model to expand a cohort of virtual patients, each with realistic disease dynamics. Our results identify key processes that regulate the immune response to SARS-CoV-2 infection in virtual patients and suggest viable therapeutic targets, underlining the importance of a rational approach to studying novel pathogens using intra-host models.