Introduction to modeling viral infections and immunity

Evaluation of drug resistance for EGFR-TKIs in lung cancer via multicellular lung-on-a-chip

Drug resistance to irreversible epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) is a primary factor affecting their therapeutic efficacy in human non-small cell lung cancer (NSCLC). NSCLC cells can undergo epithelial-mesenchymal transition (EMT) induced by many factors in the tumour microenvironment (TME), which plays a crucial role in tumour drug resistance. In this study, a multicellular lung-on-a-chip that can realise the cell co-culture of the human non-small cell lung cancer cell line HCC827, human foetal lung fibroblasts (HFL-1), and human umbilical vein endothelial cells (HUVECs) is prepared. The TME was simulated on the chip combined with perfusion and other factors, and the drug evaluation of osimertinib was performed to explore the drug resistance mechanism of EGFR-TKIs. In the early stages, a two-dimensional static cell co-culture was achieved by microchip, and the results showed that HFL-1 cells could be transformed into cancer-associated fibroblasts (CAFs), and HCC827 cells could undergo EMT, both of which were mediated by Interleukin-6 (IL-6). Vimentin (VIM) and Alpha Skeletal Muscle Actin (a-SMA) expression of HFL-1 was upregulated, whereas E-cadherin (E-cad) expression of HCC827 was down-regulated. Further, N-cadherin (N-cad) expression of HCC827 was upregulated. In both the static cell co-culture and multicellular lung-on-a-chip, HCC827 cells with CAFs co-culture or IL-6 treatment developed resistance to osimertinib. Further use of the IL-6 antibody inhibitor tocilizumab could reverse EGFR-TKI resistance to a certain extent. Combination therapy with tocilizumab and EGFR-TKIs may provide a novel therapeutic strategy for overcoming EGFR-TKI resistance caused by EMT in NSCLC. Furthermore, the lung-on-a-chip can simulate complex TME and can be used for evaluating tumour resistance and exploring mechanisms, with the potential to become an important tool for personalised diagnosis, treatment, and biomedical research.

Single-cell RNA sequencing unveils the hidden powers of zebrafish kidney for generating both hematopoiesis and adaptive antiviral immunity

The vertebrate kidneys play two evolutionary conserved roles in waste excretion and osmoregulation. Besides, the kidney of fish is considered as a functional ortholog of mammalian bone marrow that serves as a hematopoietic hub for generating blood cell lineages and immunological responses. However, knowledge about the properties of kidney hematopoietic cells, and the functionality of the kidney in fish immune systems remains to be elucidated. To this end, our present study generated a comprehensive atlas with 59 hematopoietic stem/progenitor cell (HSPC) and immune-cells types from zebrafish kidneys via single-cell transcriptome profiling analysis. These populations included almost all known cells associated with innate and adaptive immunity, and displayed differential responses to viral infection, indicating their diverse functional roles in antiviral immunity. Remarkably, HSPCs were found to have extensive reactivities to viral infection, and the trained immunity can be effectively induced in certain HSPCs. In addition, the antigen-stimulated adaptive immunity can be fully generated in the kidney, suggesting the kidney acts as a secondary lymphoid organ. These results indicated that the fish kidney is a dual-functional entity with functionalities of both primary and secondary lymphoid organs. Our findings illustrated the unique features of fish immune systems, and highlighted the multifaced biology of kidneys in ancient vertebrates.

Advanced lung organoids and lung-on-a-chip for cancer research and drug evaluation: a review

Lung cancer has become the primary cause of cancer-related deaths because of its high recurrence rate, ability to metastasise easily, and propensity to develop drug resistance. The wide-ranging heterogeneity of lung cancer subtypes increases the complexity of developing effective therapeutic interventions. Therefore, personalised diagnostic and treatment strategies are required to guide clinical practice. The advent of innovative three-dimensional (3D) culture systems such as organoid and organ-on-a-chip models provides opportunities to address these challenges and revolutionise lung cancer research and drug evaluation. In this review, we introduce the advancements in lung-related 3D culture systems, with a particular focus on lung organoids and lung-on-a-chip, and their latest contributions to lung cancer research and drug evaluation. These developments include various aspects, from authentic simulations and mechanistic enquiries into lung cancer to assessing chemotherapeutic agents and targeted therapeutic interventions. The new 3D culture system can mimic the pathological and physiological microenvironment of the lung, enabling it to supplement or replace existing two-dimensional culture models and animal experimental models and realize the potential for personalised lung cancer treatment.

Monkeypox: a review of epidemiological modelling studies and how modelling has led to mechanistic insight

Human monkeypox (mpox) virus is a viral zoonosis that belongs to the Orthopoxvirus genus of the Poxviridae family, which presents with similar symptoms as those seen in human smallpox patients. Mpox is an increasing concern globally, with over 80,000 cases in non-endemic countries as of December 2022. In this review, we provide a brief history and ecology of mpox, its basic virology, and the key differences in mpox viral fitness traits before and after 2022. We summarize and critique current knowledge from epidemiological mathematical models, within-host models, and between-host transmission models using the One Health approach, where we distinguish between models that focus on immunity from vaccination, geography, climatic variables, as well as animal models. We report various epidemiological parameters, such as the reproduction number, R0, in a condensed format to facilitate comparison between studies. We focus on how mathematical modelling studies have led to novel mechanistic insight into mpox transmission and pathogenesis. As mpox is predicted to lead to further infection peaks in many historically non-endemic countries, mathematical modelling studies of mpox can provide rapid actionable insights into viral dynamics to guide public health measures and mitigation strategies.

Towards systems immunology of critical illness at scale: from single cell 'omics to digital twins

Single-cell 'omics methodology has yielded unprecedented insights based largely on data-centric informatics for reducing, and thus interpreting, massive datasets. In parallel, parsimonious mathematical modeling based on abstractions of pathobiology has also yielded major insights into inflammation and immunity, with these models being extended to describe multi-organ disease pathophysiology as the basis of 'digital twins' and in silico clinical trials. The integration of these distinct methods at scale can drive both basic and translational advances, especially in the context of critical illness, including diseases such as COVID-19. Here, I explore achievements and argue the challenges that are inherent to the integration of data-driven and mechanistic modeling approaches, highlighting the potential of modeling-based strategies for rational immune system reprogramming.

Engineering organ-on-a-chip systems to model viral infections

Infectious diseases remain a public healthcare concern worldwide. Amidst the pandemic of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 infection, increasing resources have been diverted to investigate therapeutics targeting the COVID-19 spike glycoprotein and to develop various classes of vaccines. Most of the current investigations employ two-dimensional (2D) cell culture and animal models. However, 2D culture negates the multicellular interactions and three-dimensional (3D) microenvironment, and animal models cannot mimic human physiology because of interspecies differences. On the other hand, organ-on-a-chip (OoC) devices introduce a game-changer to model viral infections in human tissues, facilitating high-throughput screening of antiviral therapeutics. In this context, this review provides an overview of thein vitroOoC-based modeling of viral infection, highlighting the strengths and challenges for the future.

Biomimetic lung-on-a-chip to model virus infection and drug evaluation

Viral infectious diseases remain a global public health problem. The rapid and widespread spread of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV‑2) has had a severe impact on the global economy and human activities, highlighting the vulnerability of humans to viral infectious diseases and the urgent need to develop new technologies and effective treatments. Organ-on-a-chip is an emerging technology for constructing the physiological and pathological microenvironment of human organs in vitro and has the advantages of portability, high throughput, low cost, and accurate simulation of the in vivo microenvironment. Indeed, organ-on-a-chip provides a low-cost alternative for investigating human organ physiology, organ diseases, toxicology, and drug efficacy. The lung is a main target organ of viral infection, and lung pathophysiology must be assessed after viral infection and treatment with antiviral drugs. This review introduces the construction of lung-on-a-chip and its related pathophysiological models, focusing on the in vitro simulation of viral infection and evaluation of antiviral drugs, providing a developmental direction for research and treatment of viral diseases.

Estimation of timing of infection from longitudinal SARS-CoV-2 viral load data: mathematical modelling study

BACKGROUND

Multiple waves of the COVID-19 epidemic have hit most countries by the end of 2021. Most of those waves are caused by emergence and importation of new variants. To prevent importation of new variants, combination of border control and contact tracing is essential. However, the timing of infection inferred by interview is influenced by recall bias and hinders the contact tracing process.

METHODS

We propose a novel approach to infer the timing of infection, by employing a within-host model to capture viral load dynamics after the onset of symptoms. We applied this approach to ascertain secondary transmission which can trigger outbreaks. As a demonstration, the 12 initial reported cases in Singapore, which were considered as imported because of their recent travel history to Wuhan, were analyzed to assess whether they are truly imported.

RESULTS

Our approach suggested that 6 cases were infected prior to the arrival in Singapore, whereas other 6 cases might have been secondary local infection. Three among the 6 potential secondary transmission cases revealed that they had contact history to previously confirmed cases.

CONCLUSIONS

Contact trace combined with our approach using viral load data could be the key to mitigate the risk of importation of new variants by identifying cases as early as possible and inferring the timing of infection with high accuracy.

Parameter Estimation in the Age of Degeneracy and Unidentifiability

Parameter estimation from observable or experimental data is a crucial stage in any modeling study. Identifiability refers to one’s ability to uniquely estimate the model parameters from the available data. Structural unidentifiability in dynamic models, the opposite of identifiability, is associated with the notion of degeneracy where multiple parameter sets produce the same pattern. Therefore, the inverse function of determining the model parameters from the data is not well defined. Degeneracy is not only a mathematical property of models, but it has also been reported in biological experiments. Classical studies on structural unidentifiability focused on the notion that one can at most identify combinations of unidentifiable model parameters. We have identified a different type of structural degeneracy/unidentifiability present in a family of models, which we refer to as the Lambda-Omega (Λ-Ω) models. These are an extension of the classical lambda-omega (λ-ω) models that have been used to model biological systems, and display a richer dynamic behavior and waveforms that range from sinusoidal to square wave to spike like. We show that the Λ-Ω models feature infinitely many parameter sets that produce identical stable oscillations, except possible for a phase shift (reflecting the initial phase). These degenerate parameters are not identifiable combinations of unidentifiable parameters as is the case in structural degeneracy. In fact, reducing the number of model parameters in the Λ-Ω models is minimal in the sense that each one controls a different aspect of the model dynamics and the dynamic complexity of the system would be reduced by reducing the number of parameters. We argue that the family of Λ-Ω models serves as a framework for the systematic investigation of degeneracy and identifiability in dynamic models and for the investigation of the interplay between structural and other forms of unidentifiability resulting on the lack of information from the experimental/observational data.

Machine learning for mathematical models of HCV kinetics during antiviral therapy

Mathematical models for hepatitis C virus (HCV) dynamics have provided a means for evaluating the antiviral effectiveness of therapy and estimating treatment outcomes such as the time to cure. Recently, a mathematical modeling approach was used in the first proof-of-concept clinical trial assessing in real-time the utility of response-guided therapy with direct-acting antivirals (DAAs) in chronic HCV-infected patients. Several retrospective studies have shown that mathematical modeling of viral kinetics predicts time to cure of less than 12 weeks in the majority of individuals treated with sofosbuvir-based as well as other DAA regimens. A database of these studies was built, and machine learning methods were evaluated for their ability to estimate the time to cure for each patient to facilitate real-time modeling studies. Data from these studies exploring mathematical modeling of HCV kinetics under DAAs in 266 chronic HCV-infected patients were gathered. Different learning methods were applied and trained on part of the dataset ('train' set), to predict time to cure on the untrained part ('test' set). Our results show that this machine learning approach provides a means for establishing an accurate time to cure prediction that will support the implementation of individualized treatment.

Between-group comparison of area under the curve in clinical trials with censored follow-up: Application to HIV therapeutic vaccines

In clinical trials, longitudinal data are commonly analyzed and compared between groups using a single summary statistic such as area under the outcome versus time curve (AUC). However, incomplete data, arising from censoring due to a limit of detection or missing data, can bias these analyses. In this article, we present a statistical test based on splines-based mixed-model accounting for both the censoring and missingness mechanisms in the AUC estimation. Inferential properties of the proposed method were evaluated and compared to ad hoc approaches and to a non-parametric method through a simulation study based on two-armed trial where trajectories and the proportion of missing data were varied. Simulation results highlight that our approach has significant advantages over the other methods. A real working example from two HIV therapeutic vaccine trials is presented to illustrate the applicability of our approach.

Leveraging Computational Modeling to Understand Infectious Diseases

Purpose of Review

Computational and mathematical modeling have become a critical part of understanding in-host infectious disease dynamics and predicting effective treatments. In this review, we discuss recent findings pertaining to the biological mechanisms underlying infectious diseases, including etiology, pathogenesis, and the cellular interactions with infectious agents. We present advances in modeling techniques that have led to fundamental disease discoveries and impacted clinical translation.

Recent Findings

Combining mechanistic models and machine learning algorithms has led to improvements in the treatment of Shigella and tuberculosis through the development of novel compounds. Modeling of the epidemic dynamics of malaria at the within-host and between-host level has afforded the development of more effective vaccination and antimalarial therapies. Similarly, in-host and host-host models have supported the development of new HIV treatment modalities and an improved understanding of the immune involvement in influenza. In addition, large-scale transmission models of SARS-CoV-2 have furthered the understanding of coronavirus disease and allowed for rapid policy implementations on travel restrictions and contract tracing apps.

Summary

Computational modeling is now more than ever at the forefront of infectious disease research due to the COVID-19 pandemic. This review highlights how infectious diseases can be better understood by connecting scientists from medicine and molecular biology with those in computer science and applied mathematics.

Stochastic dynamics of Francisella tularensis infection and replication

We study the pathogenesis of Francisella tularensis infection with an experimental mouse model, agent-based computation and mathematical analysis. Following inhalational exposure to Francisella tularensis SCHU S4, a small initial number of bacteria enter lung host cells and proliferate inside them, eventually destroying the host cell and releasing numerous copies that infect other cells. Our analysis of disease progression is based on a stochastic model of a population of infectious agents inside one host cell, extending the birth-and-death process by the occurrence of catastrophes: cell rupture events that affect all bacteria in a cell simultaneously. Closed expressions are obtained for the survival function of an infected cell, the number of bacteria released as a function of time after infection, and the total bacterial load. We compare our mathematical analysis with the results of agent-based computation and, making use of approximate Bayesian statistical inference, with experimental measurements carried out after murine aerosol infection with the virulent SCHU S4 strain of the bacterium Francisella tularensis, that infects alveolar macrophages. The posterior distribution of the rate of replication of intracellular bacteria is consistent with the estimate that the time between rounds of bacterial division is less than 6 hours in vivo.

Infecting a host cell is required for the replication of many types of bacteria and viruses. In some cases, infected cells release new infectious agents continuously over their lifetime. In others, such as the Francisella tularensis bacterium studied here, they are released in a single burst that coincides with the cell’s death. We show how a stochastic model, the birth-and-death process with catastrophe, can be used to characterise infection in a single cell, thereby allowing us to account for burst events and quantify the kinetics of pathogenesis in the lung, the initial site of infection, as well as in other organs that the infection spreads to. We learn about the parameters of the mathematical model of Francisella tularensis infection making use of the experimental measurements of bacterial loads, together with approximate Bayesian statistical inference methods. The most important parameter describing the pathogenesis is the rate of replication of intracellular bacteria.

Model Averaging in Viral Dynamic Models

The paucity of experimental data makes both inference and prediction particularly challenging in viral dynamic models. In the presence of several candidate models, a common strategy is model selection (MS), in which models are fitted to the data but only results obtained with the "best model" are presented. However, this approach ignores model uncertainty, which may lead to inaccurate predictions. When several models provide a good fit to the data, another approach is model averaging (MA) that weights the predictions of each model according to its consistency to the data. Here, we evaluated by simulations in a nonlinear mixed-effect model framework the performances of MS and MA in two realistic cases of acute viral infection, i.e., (1) inference in the presence of poorly identifiable parameters, namely, initial viral inoculum and eclipse phase duration, (2) uncertainty on the mechanisms of action of the immune response. MS was associated in some scenarios with a large rate of false selection. This led to a coverage rate lower than the nominal coverage rate of 0.95 in the majority of cases and below 0.50 in some scenarios. In contrast, MA provided better estimation of parameter uncertainty, with coverage rates ranging from 0.72 to 0.98 and mostly comprised within the nominal coverage rate. Finally, MA provided similar predictions than those obtained with MS. In conclusion, parameter estimates obtained with MS should be taken with caution, especially when several models well describe the data. In this situation, MA has better performances and could be performed to account for model uncertainty.

A software package for immunologists to learn simulation modeling

Background

As immunology continues to become more quantitative, increasingly sophisticated computational tools are commonly used. One useful toolset are simulation models. Becoming familiar with such models and their uses generally requires writing computer code early in the learning process. This poses a barrier for individuals who do not have prior coding experience.

Results

To help reduce this barrier, I wrote software that teaches the use of mechanistic simulation models to study infection and immune response dynamics, without the need to read or write computer code. The software, called Dynamical Systems Approach to Immune Response Modeling (DSAIRM), is implemented as a freely available package for the R programming language. The target audience are immunologists and other scientists with no or little coding experience. DSAIRM provides a hands-on introduction to simulation models, teaches the basics of those models and what they can be used for. Here, I describe the DSAIRM R package, explain the different ways the package can be used, and provide a few introductory examples.

Conclusions

Working through DSAIRM will equip individuals with the knowledge needed to critically assess studies using simulation models in the published literature and will help them understand when such a modeling approach might be suitable for their own research. DSAIRM also provides users a potential starting point towards development and use of simulation models in their own research.

Simulation modelling for immunologists

The immune system is inordinately complex with many interacting components determining overall outcomes. Mathematical and computational modelling provides a useful way in which the various contributions of different immunological components can be probed in an integrated manner. Here, we provide an introductory overview and review of mechanistic simulation models. We start out by briefly defining these types of models and contrasting them to other model types that are relevant to the field of immunology. We follow with a few specific examples and then review the different ways one can use such models to answer immunological questions. While our examples focus on immune responses to infection, the overall ideas and descriptions of model uses can be applied to any area of immunology.