Agent-based modeling of host-pathogen systems: The successes and challenges

Towards early diagnosis of Alzheimer's disease: advances in immune-related blood biomarkers and computational approaches

Alzheimer's disease has an increasing prevalence in the population world-wide, yet current diagnostic methods based on recommended biomarkers are only available in specialized clinics. Due to these circumstances, Alzheimer's disease is usually diagnosed late, which contrasts with the currently available treatment options that are only effective for patients at an early stage. Blood-based biomarkers could fill in the gap of easily accessible and low-cost methods for early diagnosis of the disease. In particular, immune-based blood-biomarkers might be a promising option, given the recently discovered cross-talk of immune cells of the central nervous system with those in the peripheral immune system. Here, we give a background on recent advances in research on brain-immune system cross-talk in Alzheimer's disease and review machine learning approaches, which can combine multiple biomarkers with further information (e.g. age, sex, APOE genotype) into predictive models supporting an earlier diagnosis. In addition, mechanistic modeling approaches, such as agent-based modeling open the possibility to model and analyze cell dynamics over time. This review aims to provide an overview of the current state of immune-system related blood-based biomarkers and their potential for the early diagnosis of Alzheimer's disease.

Engineering innovations in medicine and biology: Revolutionizing patient care through mechanical solutions

The overlap between mechanical engineering and medicine is expanding more and more over the years. Engineers are now using their expertise to design and create functional biomaterials and are continually collaborating with physicians to improve patient health. In this review, we explore the state of scientific knowledge in the areas of biomaterials, biomechanics, nanomechanics, and computational fluid dynamics (CFD) in relation to the pharmaceutical and medical industry. Focusing on current research and breakthroughs, we provide an overview of how these fields are being used to create new technologies for medical treatments of human patients. Barriers and constraints in these fields, as well as ways to overcome them, are also described in this review. Finally, the potential for future advances in biomaterials to fundamentally change the current approach to medicine and biology is also discussed.

Emerging advances in biosecurity to underpin human, animal, plant, and ecosystem health

One Biosecurity is an interdisciplinary approach to policy and research that builds on the interconnections between human, animal, plant, and ecosystem health to effectively prevent and mitigate the impacts of invasive alien species. To support this approach requires that key cross-sectoral research innovations be identified and prioritized. Following an interdisciplinary horizon scan for emerging research that underpins One Biosecurity, four major interlinked advances were identified: implementation of new surveillance technologies adopting state-of-the-art sensors connected to the Internet of Things, deployable handheld molecular and genomic tracing tools, the incorporation of wellbeing and diverse human values into biosecurity decision-making, and sophisticated socio-environmental models and data capture. The relevance and applicability of these innovations to address threats from pathogens, pests, and weeds in both terrestrial and aquatic ecosystems emphasize the opportunity to build critical mass around interdisciplinary teams at a global scale that can rapidly advance science solutions targeting biosecurity threats.

Quantitative approaches for decoding the specificity of the human T cell repertoire

T cell receptor (TCR)-peptide-major histocompatibility complex (pMHC) interactions play a vital role in initiating immune responses against pathogens, and the specificity of TCRpMHC interactions is crucial for developing optimized therapeutic strategies. The advent of high-throughput immunological and structural evaluation of TCR and pMHC has provided an abundance of data for computational approaches that aim to predict favorable TCR-pMHC interactions. Current models are constructed using information on protein sequence, structures, or a combination of both, and utilize a variety of statistical learning-based approaches for identifying the rules governing specificity. This review examines the current theoretical, computational, and deep learning approaches for identifying TCR-pMHC recognition pairs, placing emphasis on each method's mathematical approach, predictive performance, and limitations.

Verifiable biology

The formalization of biological systems using computational modelling approaches as an alternative to mathematical-based methods has recently received much interest because computational models provide a deeper mechanistic understanding of biological systems. In particular, formal verification, complementary approach to standard computational techniques such as simulation, is used to validate the system correctness and obtain critical information about system behaviour. In this study, we survey the most frequently used computational modelling approaches and formal verification techniques for computational biology. We compare a number of verification tools and software suites used to analyse biological systems and biochemical networks, and to verify a wide range of biological properties. For users who have no expertise in formal verification, we present a novel methodology that allows them to easily apply formal verification techniques to analyse their biological or biochemical system of interest.

Systems biology approaches to investigate the role of granulomas in TB-HIV coinfection

The risk of active tuberculosis disease is 15-21 times higher in those coinfected with human immunodeficiency virus-1 (HIV) compared to tuberculosis alone, and tuberculosis is the leading cause of death in HIV+ individuals. Mechanisms driving synergy between Mycobacterium tuberculosis (Mtb) and HIV during coinfection include: disruption of cytokine balances, impairment of innate and adaptive immune cell functionality, and Mtb-induced increase in HIV viral loads. Tuberculosis granulomas are the interface of host-pathogen interactions. Thus, granuloma-based research elucidating the role and relative impact of coinfection mechanisms within Mtb granulomas could inform cohesive treatments that target both pathogens simultaneously. We review known interactions between Mtb and HIV, and discuss how the structure, function and development of the granuloma microenvironment create a positive feedback loop favoring pathogen expansion and interaction. We also identify key outstanding questions and highlight how coupling computational modeling with in vitro and in vivo efforts could accelerate Mtb-HIV coinfection discoveries.

Deterministic Response Threshold Models of Reproductive Division of Labor Are More Robust Than Probabilistic Models in Artificial Ants

We implement an agent-based simulation of the response threshold model of reproductive division of labor. Ants in our simulation must perform two tasks in their environment: forage and reproduce. The colony is capable of allocating ant resources to these roles using different division of labor strategies via genetic architectures and plasticity mechanisms. We find that the deterministic allocation strategy of the response threshold model is more robust than the probabilistic allocation strategy. The deterministic allocation strategy is also capable of evolving complex solutions to colony problems like niche construction and recovery from the loss of the breeding caste. In addition, plasticity mechanisms had both positive and negative influence on the emergence of reproductive division of labor. The combination of plasticity mechanisms has an additive and sometimes emergent impact.

A global method for fast simulations of molecular dynamics in multiscale agent-based models of biological tissues

Agent-based models (ABMs) are a natural platform for capturing the multiple time and spatial scales in biological processes. However, these models are computationally expensive, especially when including molecular-level effects. The traditional approach to simulating this type of multiscale ABM is to solve a system of ordinary differential equations for the molecular events per cell. This significantly adds to the computational cost of simulations as the number of agents grows, which contributes to many ABMs being limited to around10 5 cells. We propose an approach that requires the same computational time independent of the number of agents. This speeds up the entire simulation by orders of magnitude, allowing for more thorough explorations of ABMs with even larger numbers of agents. We use two systems to show that the new method strongly agrees with the traditionally used approach. This computational strategy can be applied to a wide range of biological investigations.

An Agent-Based Interpretation of Leukocyte Chemotaxis in Cancer-on-Chip Experiments

The present paper was inspired by recent developments in laboratory experiments within the framework of cancer-on-chip technology, an immune-oncology microfluidic chip aiming at studying the fundamental mechanisms of immunocompetent behavior. We focus on the laboratory setting where cancer is treated with chemotherapy drugs, and in this case, the effects of the treatment administration hypothesized by biologists are: the absence of migration and proliferation of tumor cells, which are dying; the stimulation of the production of chemical substances (annexin); the migration of leukocytes in the direction of higher concentrations of chemicals. Here, following the physiological hypotheses made by biologists on the phenomena occurring in these experiments, we introduce an agent-based model reproducing the dynamics of two cell populations (agents), i.e., tumor cells and leukocytes living in the microfluidic chip environment. Our model aims at proof of concept, demonstrating that the observations of the biological phenomena can be obtained by the model on the basis of the explicit assumptions made. In this framework, close adherence of the computational model to the biological results, as shown in the section devoted to the first calibration of the model with respect to available observations, is successfully accomplished.

Computational models in inflammatory bowel disease

Inflammatory bowel disease (IBD) is a chronic and relapsing disease with multiple underlying influences and notable heterogeneity among its clinical and response-to-treatment phenotypes. There is no cure for IBD, and none of the currently available therapies have demonstrated clinical efficacies beyond 40%-60%. Data collected about its omics, pathogenesis, and treatment strategies have grown exponentially with time making IBD a prime candidate for artificial intelligence (AI) mediated discovery support. AI can be leveraged to further understand or identify IBD features to improve clinical outcomes. Various treatment candidates are currently under evaluation in clinical trials, offering further approaches and opportunities for increasing the efficacies of treatments. However, currently, therapeutic plans are largely determined using clinical features due to the lack of specific biomarkers, and it has become necessary to step into precision medicine to predict therapeutic responses to guarantee optimal treatment efficacy. This is accompanied by the application of AI and the development of multiscale hybrid models combining mechanistic approaches and machine learning. These models ultimately lead to the creation of digital twins of given patients delivering on the promise of precision dosing and tailored treatment. Interleukin-6 (IL-6) is a prominent cytokine in cell-to-cell communication in the inflammatory responses' regulation. Dysregulated IL-6-induced signaling leads to severe immunological or proliferative pathologies, such as IBD and colon cancer. This mini-review explores multiscale models with the aim of predicting the response to therapy in IBD. Modeling IL-6 biology and generating digital twins enhance the credibility of their prediction.

Spatially distributed infection increases viral load in a computational model of SARS-CoV-2 lung infection

A key question in SARS-CoV-2 infection is why viral loads and patient outcomes vary dramatically across individuals. Because spatial-temporal dynamics of viral spread and immune response are challenging to study in vivo, we developed Spatial Immune Model of Coronavirus (SIMCoV), a scalable computational model that simulates hundreds of millions of lung cells, including respiratory epithelial cells and T cells. SIMCoV replicates viral growth dynamics observed in patients and shows how spatially dispersed infections can lead to increased viral loads. The model also shows how the timing and strength of the T cell response can affect viral persistence, oscillations, and control. By incorporating spatial interactions, SIMCoV provides a parsimonious explanation for the dramatically different viral load trajectories among patients by varying only the number of initial sites of infection and the magnitude and timing of the T cell immune response. When the branching airway structure of the lung is explicitly represented, we find that virus spreads faster than in a 2D layer of epithelial cells, but much more slowly than in an undifferentiated 3D grid or in a well-mixed differential equation model. These results illustrate how realistic, spatially explicit computational models can improve understanding of within-host dynamics of SARS-CoV-2 infection.

A runtime alterable epidemic model with genetic drift, waning immunity and vaccinations

We present methods for building a Java Runtime-Alterable-Model Platform (RAMP) of complex dynamical systems. We illustrate our methods by building a multivariant SEIR (epidemic) RAMP. Underlying our RAMP is an individual-based model that includes adaptive contact rates, pathogen genetic drift, waning and cross-immunity. Besides allowing parameter values, process descriptions and scriptable runtime drivers to be easily modified during simulations, our RAMP can used within R-Studio and other computational platforms. Process descriptions that can be runtime altered within our SEIR RAMP include pathogen variant-dependent host shedding, environmental persistence, host transmission and within-host pathogen mutation and replication. They also include adaptive social distancing and adaptive application of vaccination rates and variant-valency of vaccines. We present simulation results using parameter values and process descriptions relevant to the current COVID-19 pandemic. Our results suggest that if waning immunity outpaces vaccination rates, then vaccination rollouts may fail to contain the most transmissible variants, particularly if vaccine valencies are not adapted to deal with escape mutations. Our SEIR RAMP is designed for easy use by others. More generally, our RAMP concept facilitates construction of highly flexible complex systems models of all types, which can then be easily shared as stand-alone application programs.

Iterative community-driven development of a SARS-CoV-2 tissue simulator

The 2019 novel coronavirus, SARS-CoV-2, is a pathogen of critical significance to international public health. Knowledge of the interplay between molecular-scale virus-receptor interactions, single-cell viral replication, intracellular-scale viral transport, and emergent tissue-scale viral propagation is limited. Moreover, little is known about immune system-virus-tissue interactions and how these can result in low-level (asymptomatic) infections in some cases and acute respiratory distress syndrome (ARDS) in others, particularly with respect to presentation in different age groups or pre-existing inflammatory risk factors. Given the nonlinear interactions within and among each of these processes, multiscale simulation models can shed light on the emergent dynamics that lead to divergent outcomes, identify actionable "choke points" for pharmacologic interventions, screen potential therapies, and identify potential biomarkers that differentiate patient outcomes. Given the complexity of the problem and the acute need for an actionable model to guide therapy discovery and optimization, we introduce and iteratively refine a prototype of a multiscale model of SARS-CoV-2 dynamics in lung tissue. The first prototype model was built and shared internationally as open source code and an online interactive model in under 12 hours, and community domain expertise is driving regular refinements. In a sustained community effort, this consortium is integrating data and expertise across virology, immunology, mathematical biology, quantitative systems physiology, cloud and high performance computing, and other domains to accelerate our response to this critical threat to international health. More broadly, this effort is creating a reusable, modular framework for studying viral replication and immune response in tissues, which can also potentially be adapted to related problems in immunology and immunotherapy.

Thymic B Cells as a New Player in the Type 1 Diabetes Response

Type 1 diabetes (T1d) results from a sustained autoreactive T and B cell response towards insulin-producing β cells in the islets of Langerhans. The autoreactive nature of the condition has led to many investigations addressing the genetic or cellular changes in primary lymphoid tissues that impairs central tolerance- a key process in the deletion of autoreactive T and B cells during their development. For T cells, these studies have largely focused on medullary thymic epithelial cells (mTECs) critical for the effective negative selection of autoreactive T cells in the thymus. Recently, a new cellular player that impacts positively or negatively on the deletion of autoreactive T cells during their development has come to light, thymic B cells. Normally a small population within the thymus of mouse and man, thymic B cells expand in T1d as well as other autoimmune conditions, reside in thymic ectopic germinal centres and secrete autoantibodies that bind selective mTECs precipitating mTEC death. In this review we will discuss the ontogeny, characteristics and functionality of thymic B cells in healthy and autoimmune settings. Furthermore, we explore how in silico approaches may help decipher the complex cellular interplay of thymic B cells with other cells within the thymic microenvironment leading to new avenues for therapeutic intervention.

SBMLWebApp: Web-Based Simulation, Steady-State Analysis, and Parameter Estimation of Systems Biology Models

In systems biology, biological phenomena are often modeled by Ordinary Differential Equations (ODEs) and distributed in the de facto standard file format SBML. The primary analyses performed with such models are dynamic simulation, steady-state analysis, and parameter estimation. These methodologies are mathematically formalized, and libraries for such analyses have been published. Several tools exist to create, simulate, or visualize models encoded in SBML. However, setting up and establishing analysis environments is a crucial hurdle for non-modelers. Therefore, easy access to perform fundamental analyses of ODE models is a significant challenge. We developed SBMLWebApp, a web-based service to execute SBML-based simulation, steady-state analysis, and parameter estimation directly in the browser without the need for any setup or prior knowledge to address this issue. SBMLWebApp visualizes the result and numerical table of each analysis and provides a download of the results. SBMLWebApp allows users to select and analyze SBML models directly from the BioModels Database. Taken together, SBMLWebApp provides barrier-free access to an SBML analysis environment for simulation, steady-state analysis, and parameter estimation for SBML models. SBMLWebApp is implemented in Java™ based on an Apache Tomcat® web server using COPASI, the Systems Biology Simulation Core Library (SBSCL), and LibSBMLSim as simulation engines. SBMLWebApp is licensed under MIT with source code freely available. At the end of this article, the Data Availability Statement gives the internet links to the two websites to find the source code and run the program online.

A Review of Using Mathematical Modeling to Improve Our Understanding of Bacteriophage, Bacteria, and Eukaryotic Interactions

Phage therapy, the therapeutic usage of viruses to treat bacterial infections, has many theoretical benefits in the ‘post antibiotic era.’ Nevertheless, there are currently no approved mainstream phage therapies. One reason for this is a lack of understanding of the complex interactions between bacteriophage, bacteria and eukaryotic hosts. These three-component interactions are complex, with non-linear or synergistic relationships, anatomical barriers and genetic or phenotypic heterogeneity all leading to disparity between performance and efficacy in in vivo versus in vitro environments. Realistic computer or mathematical models of these complex environments are a potential route to improve the predictive power of in vitro studies for the in vivo environment, and to streamline lab work. Here, we introduce and review the current status of mathematical modeling and highlight that data on genetic heterogeneity and mutational stochasticity, time delays and population densities could be critical in the development of realistic phage therapy models in the future. With this in mind, we aim to inform and encourage the collaboration and sharing of knowledge and expertise between microbiologists and theoretical modelers, synergising skills and smoothing the road to regulatory approval and widespread use of phage therapy.

Epidemiological models and COVID-19: a comparative view

Epidemiological models have played a central role in the COVID-19 pandemic, particularly when urgent decisions were required and available evidence was sparse. They have been used to predict the evolution of the disease and to inform policy-making. In this paper, we address two kinds of epidemiological models widely used in the pandemic, namely, compartmental models and agent-based models. After describing their essentials-some real examples are invoked-we discuss their main strengths and weaknesses. Then, on the basis of this analysis, we make a comparison between their respective merits concerning three different goals: prediction, explanation, and intervention. We argue that there are general considerations which could favour any of those sorts of models for obtaining the aforementioned goals. We conclude, however, that preference for particular models must be grounded case-by-case since additional contextual factors, as the peculiarities of the target population and the aims and expectations of policy-makers, cannot be overlooked.

Step-by-step comparison of ordinary differential equation and agent-based approaches to pharmacokinetic-pharmacodynamic models

Mathematical models in oncology aid in the design of drugs and understanding of their mechanisms of action by simulation of drug biodistribution, drug effects, and interaction between tumor and healthy cells. The traditional approach in pharmacometrics is to develop and validate ordinary differential equation models to quantify trends at the population level. In this approach, time‐course of biological measurements is modeled continuously, assuming a homogenous population. Another approach, agent‐based models, focuses on the behavior and fate of biological entities at the individual level, which subsequently could be summarized to reflect the population level. Heterogeneous cell populations and discrete events are simulated, and spatial distribution can be incorporated. In this tutorial, an agent‐based model is presented and compared to an ordinary differential equation model for a tumor efficacy model inhibiting the pERK pathway. We highlight strengths, weaknesses, and opportunities of each approach.

Examining the interplay between face mask usage, asymptomatic transmission, and social distancing on the spread of COVID-19

COVID-19's high virus transmission rates have caused a pandemic that is exacerbated by the high rates of asymptomatic and presymptomatic infections. These factors suggest that face masks and social distance could be paramount in containing the pandemic. We examined the efficacy of each measure and the combination of both measures using an agent-based model within a closed space that approximated real-life interactions. By explicitly considering different fractions of asymptomatic individuals, as well as a realistic hypothesis of face masks protection during inhaling and exhaling, our simulations demonstrate that a synergistic use of face masks and social distancing is the most effective intervention to curb the infection spread. To control the pandemic, our models suggest that high adherence to social distance is necessary to curb the spread of the disease, and that wearing face masks provides optimal protection even if only a small portion of the population comply with social distance. Finally, the face mask effectiveness in curbing the viral spread is not reduced if a large fraction of population is asymptomatic. Our findings have important implications for policies that dictate the reopening of social gatherings.

An agent-based model for simulating COVID-19 transmissions on university campus and its implications on mitigation interventions: a case study

Purpose

Universities across the USA are facing challenging decision-making problems amid the COVID-19 pandemic. The purpose of this study is to facilitate universities in planning disease mitigation interventions as they respond to the pandemic.

Design/methodology/approach

An agent-based model is developed to mimic the virus transmission dynamics on campus. Scenario-based experiments are conducted to evaluate the effectiveness of various interventions including course modality shift (from face-to-face to online), social distancing, mask use and vaccination. A case study is performed for a typical US university.

Findings

With 10%, 30%, 50%, 70% and 90% course modality shift, the number of total cases can be reduced to 3.9%, 20.9%, 35.6%, 60.9% and 96.8%, respectively, comparing against the baseline scenario (no interventions). More than 99.9% of the total infections can be prevented when combined social distancing and mask use are implemented even without course modality shift. If vaccination is implemented without other interventions, the reductions are 57.1%, 90.6% and 99.6% with 80%, 85% and 90% vaccine efficacies, respectively. In contrast, more than 99% reductions are found with all three vaccine efficacies if mask use is combined.

Practical implications

This study provides useful implications for supporting universities in mitigating transmissions on campus and planning operations for the upcoming semesters.

Originality/value

An agent-based model is developed to investigate COVID-19 transmissions on campus and evaluate the effectiveness of various mitigation interventions.

A modular framework for multiscale, multicellular, spatiotemporal modeling of acute primary viral infection and immune response in epithelial tissues and its application to drug therapy timing and effectiveness

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

COVID-19 Intervention Scenarios for a Long-term Disease Management

BACKGROUND

The first outbreak of coronavirus disease 2019 (COVID-19) was successfully restrained in many countries around the world by means of a severe lockdown. Now, we are entering the second phase of the pandemics in which the spread of the virus needs to be contained within the limits that national health systems can cope with. This second phase of the epidemics is expected to last until a vaccination is available or herd immunity is reached. Long-term management strategies thus need to be developed.

METHODS

In this paper we present a new agent-based simulation model "COVID-19 ABM" with which we simulate 4 alternative scenarios for the second "new normality" phase that can help decision-makers to take adequate control and intervention measures.

RESULTS

The scenarios resulted in distinctly different outcomes. A continued lockdown could regionally eradicate the virus within a few months, whereas a relaxation back to 80% of former activity-levels was followed by a second outbreak. Contact-tracing as well as adaptive response strategies could keep COVID-19 within limits.

CONCLUSION

The main insights are that low-level voluntary use of tracing apps shows no relevant effects on containing the virus, whereas medium or high-level tracing allows maintaining a considerably higher level of social activity. Adaptive control strategies help in finding the level of least restrictions. A regional approach to adaptive management can further help in fine-tuning the response to regional dynamics and thus minimise negative economic effects.

Modeling the efficiency of filovirus entry into cells in vitro: Effects of SNP mutations in the receptor molecule

Interaction between filovirus glycoprotein (GP) and the Niemann-Pick C1 (NPC1) protein is essential for membrane fusion during virus entry. Some single-nucleotide polymorphism (SNPs) in two surface-exposed loops of NPC1 are known to reduce viral infectivity. However, the dependence of differences in entry efficiency on SNPs remains unclear. Using vesicular stomatitis virus pseudotyped with Ebola and Marburg virus GPs, we investigated the cell-to-cell spread of viruses in cultured cells expressing NPC1 or SNP derivatives. Eclipse and virus-producing phases were assessed by in vitro infection experiments, and we developed a mathematical model describing spatial-temporal virus spread. This mathematical model fit the plaque radius data well from day 2 to day 6. Based on the estimated parameters, we found that SNPs causing the P424A and D508N substitutions in NPC1 most effectively reduced the entry efficiency of Ebola and Marburg viruses, respectively. Our novel approach could be broadly applied to other virus plaque assays.

Ebola (EBOV) and Marburg (MARV) viruses, which are included viruses of the family Filoviridae, cause severe hemorrhagic fever in humans. Filovirus particles is adsorbed to the cell through glycoprotein (GP), which is the only viral surface protein. Interaction between the filovirus sugar protein (GP) and the Niemann-Pick C1 (NPC1) protein plays a key role in membrane fusion during virus entry. Although some single-nucleotide polymorphism (SNPs) in two surface-exposed loops of NPC1 are known to reduce viral infectivity, the dependence of differences in entry efficiency on SNPs has not been studied. We therefore investigated the cell-to-cell spread of viruses in cultured cells expressing NPC1 or SNP derivatives. Using a mathematical model describing spatial-temporal virus spread, we quantitatively analyze viral entry efficiency and how this affected cell-to-cell spread. Our approach may be applied to not only understanding the roles of genetic polymorphisms in human susceptibility to filoviruses, but also other virus plaque assays.

A modular framework for multiscale, multicellular, spatiotemporal modeling of acute primary viral infection and immune response in epithelial tissues and its application to drug therapy timing and effectiveness: A multiscale model of viral infection in epithelial tissues

Simulations of tissue-specific effects of primary acute viral infections like COVID-19 are essential for understanding disease outcomes and optimizing therapies. Such simulations need to support continuous updating in response to rapid advances in understanding of infection mechanisms, and parallel development of components by multiple groups. We present an open-source platform for multiscale spatiotemporal simulation of an epithelial tissue, viral infection, cellular immune response and tissue damage, specifically designed to be modular and extensible to support continuous updating and parallel development. The base simulation of a simplified patch of epithelial tissue and immune response exhibits distinct patterns of infection dynamics from widespread infection, to recurrence, to clearance. Slower viral internalization and faster immune-cell recruitment slow infection and promote containment. Because antiviral drugs can have side effects and show reduced clinical effectiveness when given later during infection, we studied the effects on progression of treatment potency and time-of-first treatment after infection. In simulations, even a low potency therapy with a drug which reduces the replication rate of viral RNA greatly decreases the total tissue damage and virus burden when given near the beginning of infection. Many combinations of dosage and treatment time lead to stochastic outcomes, with some simulation replicas showing clearance or control (treatment success), while others show rapid infection of all epithelial cells (treatment failure). Thus, while a high potency therapy usually is less effective when given later, treatments at late times are occasionally effective. We illustrate how to extend the platform to model specific virus types (e.g., hepatitis C) and add additional cellular mechanisms (tissue recovery and variable cell susceptibility to infection), using our software modules and publicly-available software repository.

Tissue-resident T cell-derived cytokines eliminate herpes simplex virus-2-infected cells

The mechanisms underlying rapid elimination of herpes simplex virus-2 (HSV-2) in the human genital tract despite low CD8+ and CD4+ tissue-resident T cell (Trm cell) density are unknown. We analyzed shedding episodes during chronic HSV-2 infection; viral clearance always predominated within 24 hours of detection even when viral load exceeded 1 × 107 HSV DNA copies, and surges in granzyme B and IFN-γ occurred within the early hours after reactivation and correlated with local viral load. We next developed an agent-based mathematical model of an HSV-2 genital ulcer to integrate mechanistic observations of Trm cells in in situ proliferation, trafficking, cytolytic effects, and cytokine alarm signaling from murine studies with viral kinetics, histopathology, and lesion size data from humans. A sufficiently high density of HSV-2-specific Trm cells predicted rapid elimination of infected cells, but our data suggest that such Trm cell densities are relatively uncommon in infected tissues. At lower, more commonly observed Trm cell densities, Trm cells must initiate a rapidly diffusing, polyfunctional cytokine response with activation of bystander T cells in order to eliminate a majority of infected cells and eradicate briskly spreading HSV-2 infection.

Coupling the Macroscale to the Microscale in a Spatiotemporal Context to Examine Effects of Spatial Diffusion on Disease Transmission

There are many challenges to coupling the macroscale to the microscale in temporal or spatial contexts. In order to examine effects of an individual movement and spatial control measures on a disease outbreak, we developed a multiscale model and extended the semi-stochastic simulation method by linking individual movements to pathogen’s diffusion, linking the slow dynamics for disease transmission at the population level to the fast dynamics for pathogen shedding/excretion at the individual level. Numerical simulations indicate that during a disease outbreak individuals with the same infection status show the property of clustering and, in particular, individuals’ rapid movements lead to an increase in the average reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document}R0, the final size and the peak value of the outbreak. It is interesting that a high level of aggregation the individuals’ movement results in low new infections and a small final size of the infected population. Further, we obtained that either high diffusion rate of the pathogen or frequent environmental clearance lead to a decline in the total number of infected individuals, indicating the need for control measures such as improving air circulation or environmental hygiene. We found that the level of spatial heterogeneity when implementing control greatly affects the control efficacy, and in particular, an uniform isolation strategy leads to low a final size and small peak, compared with local measures, indicating that a large-scale isolation strategy with frequent clearance of the environment is beneficial for disease control.

Environmental Restrictions: A New Concept Governing HIV-1 Spread Emerging from Integrated Experimental-Computational Analysis of Tissue-Like 3D Cultures

HIV-1 can use cell-free and cell-associated transmission modes to infect new target cells, but how the virus spreads in the infected host remains to be determined. We recently established 3D collagen cultures to study HIV-1 spread in tissue-like environments and applied iterative cycles of experimentation and computation to develop a first in silico model to describe the dynamics of HIV-1 spread in complex tissue. These analyses (i) revealed that 3D collagen environments restrict cell-free HIV-1 infection but promote cell-associated virus transmission and (ii) defined that cell densities in tissue dictate the efficacy of these transmission modes for virus spread. In this review, we discuss, in the context of the current literature, the implications of this study for our understanding of HIV-1 spread in vivo, which aspects of in vivo physiology this integrated experimental-computational analysis takes into account, and how it can be further improved experimentally and in silico.

A mathematical model describing the localization and spread of influenza A virus infection within the human respiratory tract

Within the human respiratory tract (HRT), virus diffuses through the periciliary fluid (PCF) bathing the epithelium. But virus also undergoes advection: as the mucus layer sitting atop the PCF is pushed along by the ciliated cell's beating cilia, the PCF and its virus content are also pushed along, upwards towards the nose and mouth. While many mathematical models (MMs) have described the course of influenza A virus (IAV) infections in vivo, none have considered the impact of both diffusion and advection on the kinetics and localization of the infection. The MM herein represents the HRT as a one-dimensional track extending from the nose down towards the lower HRT, wherein stationary cells interact with IAV which moves within (diffusion) and along with (advection) the PCF. Diffusion was found to be negligible in the presence of advection which effectively sweeps away IAV, preventing infection from disseminating below the depth at which virus first deposits. Higher virus production rates (10-fold) are required at higher advection speeds (40 μm/s) to maintain equivalent infection severity and timing. Because virus is entrained upwards, upper parts of the HRT see more virus than lower parts. As such, infection peaks and resolves faster in the upper than in the lower HRT, making it appear as though infection progresses from the upper towards the lower HRT, as reported in mice. When the spatial MM is expanded to include cellular regeneration and an immune response, it reproduces tissue damage levels reported in patients. It also captures the kinetics of seasonal and avian IAV infections, via parameter changes consistent with reported differences between these strains, enabling comparison of their treatment with antivirals. This new MM offers a convenient and unique platform from which to study the localization and spread of respiratory viral infections within the HRT.

A case for the reuse and adaptation of mechanistic computational models to study transplant immunology

Computational mechanistic models constitute powerful tools for summarizing our knowledge in quantitative terms, providing mechanistic understanding, and generating new hypotheses. The present review emphasizes the advantages of reusing publicly available computational models as a way to capitalize on existing knowledge, reduce the number of parameters that need to be adjusted to experimental data, and facilitate hypothesis generation. Finally, it includes a step-by-step example of the reuse and adaptation of an existing model of immune responses to tuberculosis, tumor growth, and blood pathogens, to study donor-specific antibody (DSA) responses. This review aims to illustrate the benefit of leveraging the currently available computational models in immunology to accelerate the study of alloimmune responses, and to encourage modelers to share their models to further advance our understanding of transplant immunology.

Trends in mathematical modeling of host-pathogen interactions

Pathogenic microorganisms entail enormous problems for humans, livestock, and crop plants. A better understanding of the different infection strategies of the pathogens enables us to derive optimal treatments to mitigate infectious diseases or develop vaccinations preventing the occurrence of infections altogether. In this review, we highlight the current trends in mathematical modeling approaches and related methods used for understanding host-pathogen interactions. Since these interactions can be described on vastly different temporal and spatial scales as well as abstraction levels, a variety of computational and mathematical approaches are presented. Particular emphasis is placed on dynamic optimization, game theory, and spatial modeling, as they are attracting more and more interest in systems biology. Furthermore, these approaches are often combined to illuminate the complexities of the interactions between pathogens and their host. We also discuss the phenomena of molecular mimicry and crypsis as well as the interplay between defense and counter defense. As a conclusion, we provide an overview of method characteristics to assist non-experts in their decision for modeling approaches and interdisciplinary understanding.

A personalized computational model of edema formation in myocarditis based on long-axis biventricular MRI images

Background

Myocarditis is defined as the inflammation of the myocardium, i.e. the cardiac muscle. Among the reasons that lead to this disease, we may include infections caused by a virus, bacteria, protozoa, fungus, and others. One of the signs of the inflammation is the formation of edema, which may be a consequence of the interaction between interstitial fluid dynamics and immune response. This complex physiological process was mathematically modeled using a nonlinear system of partial differential equations (PDE) based on porous media approach. By combing a model based on Biot’s poroelasticity theory with a model for the immune response we developed a new hydro-mechanical model for inflammatory edema. To verify this new computational model, T2 parametric mapping obtained by Magnetic Resonance (MR) imaging was used to identify the region of edema in a patient diagnosed with unspecific myocarditis.

Results

A patient-specific geometrical model was created using MRI images from the patient with myocarditis. With this model, edema formation was simulated using the proposed hydro-mechanical mathematical model in a two-dimensional domain. The computer simulations allowed us to correlate spatiotemporal dynamics of representative cells of the immune systems, such as leucocytes and the pathogen, with fluid accumulation and cardiac tissue deformation.

Conclusions

This study demonstrates that the proposed mathematical model is a very promising tool to better understand edema formation in myocarditis. Simulations obtained from a patient-specific model reproduced important aspects related to the formation of cardiac edema, its area, position, and shape, and how these features are related to immune response.

Modelling cell surface dynamics and cell-cell interactions using Cell Studio: a three-dimensional visualization tool based on gaming technology

Predictive modelling of complex biological systems and biophysical interactions requires the inclusion of multiple nano- and micro-scale events. In many scenarios, however, numerical solutions alone do not necessarily enhance the understanding of the system. Instead, this work explores the use of an agent-based model with visualization capabilities to elucidate interactions between single cells. We present a model of juxtacrine signalling, using Cell Studio, an agent-based modelling system, based on gaming and three-dimensional visualization tools. The main advantages of the system are its ability to apply any cell geometry and to dynamically visualize the diffusion and interactions of the molecules within the cells in real time. These provide an excellent tool for obtaining insight about different biological scenarios, as the user may view the dynamics of a system and observe its emergent behaviour as it unfolds. The agent-based model was validated against the results of a mean-field model of Notch receptors and ligands in two neighbouring cells. The conversion to an agent-based model is described in detail. To demonstrate the advantages of the model, we further created a filopodium-mediated signalling model. Our model revealed that diffusion and endocytosis alone are insufficient to produce significant signalling in a filopodia scenario. This is due to the bottleneck at the cell-filopodium contact region and the long distance to the end of the filopodium. However, allowing active transport of ligands into filopodia enhances the signalling significantly compared with a face-to-face scenario. We conclude that the agent-based approach can provide insights into mechanisms underlying cell signalling. The open-source model can be found in the Internet hosting service GitHub.

On the efficiency of chemotactic pursuit - Comparing blind search with temporal and spatial gradient sensing

In chemotaxis, cells are modulating their migration patterns in response to concentration gradients of a guiding substance. Immune cells are believed to use such chemotactic sensing for remotely detecting and homing in on pathogens. Considering that immune cells may encounter a multitude of targets with vastly different migration properties, ranging from immobile to highly mobile, it is not clear which strategies of chemotactic pursuit are simultaneously efficient and versatile. We tackle this problem theoretically and define a tunable response function that maps temporal or spatial concentration gradients to migration behavior. The seven free parameters of this response function are optimized numerically with the objective of maximizing search efficiency against a wide spectrum of target cell properties. Finally, we reverse-engineer the best-performing parameter sets to uncover strategies of chemotactic pursuit that are efficient under different biologically realistic boundary conditions. Although strategies based on the temporal or spatial sensing of chemotactic gradients are significantly more efficient than unguided migration, such ‘blind search’ turns out to work surprisingly well, in particular if the immune cells are fast and directionally persistent. The resulting simulated data can be used for the design of chemotaxis experiments and for the development of algorithms that automatically detect and quantify goal oriented behavior in measured immune cell trajectories.

Data analysis to modeling to building theory in NK cell biology and beyond: How can computational modeling contribute?

The use of mathematical and computational tools in investigating Natural Killer (NK) cell biology and in general the immune system has increased steadily in the last few decades. However, unlike the physical sciences, there is a persistent ambivalence, which however is increasingly diminishing, in the biology community toward appreciating the utility of quantitative tools in addressing questions of biological importance. We survey some of the recent developments in the application of quantitative approaches for investigating different problems in NK cell biology and evaluate opportunities and challenges of using quantitative methods in providing biological insights in NK cell biology.

Multi-agent simulation regulated by microbe-oriented thermodynamics and kinetics equations for exploiting interspecies dynamics and evolution between methanogenesis, sulfidogenesis, hydrogenesis and exoelectrogenesis

Multi-agent simulation (MAS) regulated by microbe-oriented thermodynamics and kinetics equations were performed for exploiting the interspecies dynamics and evolution in anaerobic respiration and bioelectrochemical systems. A newly-defined kinetically thermodynamic parameter is recognized microbes as agents in various conditions, including electron donors and acceptors, temperature, pH, etc. For verification of the MAS, the treatment of synthetic wastewater containing glucose and acetate was evaluated in four 25°C laboratory-scale reactors with different electron acceptors and cathode materials that had potential for methanogenesis, hydrogenesis, sulfidogenesis and exoelectrogenesis. Within 1000 h operation, the reactors performance and microbial structures using 16S rRNA sequencing matched with the MAS, suggesting acetoclastic exoelectrogenesis predominance (Geobacter). After 2400 h, MAS observed the co-existence of acetoclastic methanogenesis and acetoclastic and propionate exoelectrogenesis, as was reported previously. Such microbial evolution from the short-term to long-term operation likely resulted from the glucose-driven propionate. The MAS developed is applicable in a wide range of complex engineering and natural ecosystems.

Review of the systems biology of the immune system using agent-based models

The immune system is an inherent protection system in vertebrate animals including human beings that exhibit properties such as self-organisation, self-adaptation, learning, and recognition. It interacts with the other allied systems such as the gut and lymph nodes. There is a need for immune system modelling to know about its complex internal mechanism, to understand how it maintains the homoeostasis, and how it interacts with the other systems. There are two types of modelling techniques used for the simulation of features of the immune system: equation-based modelling (EBM) and agent-based modelling. Owing to certain shortcomings of the EBM, agent-based modelling techniques are being widely used. This technique provides various predictions for disease causes and treatments; it also helps in hypothesis verification. This study presents a review of agent-based modelling of the immune system and its interactions with the gut and lymph nodes. The authors also review the modelling of immune system interactions during tuberculosis and cancer. In addition, they also outline the future research directions for the immune system simulation through agent-based techniques such as the effects of stress on the immune system, evolution of the immune system, and identification of the parameters for a healthy immune system.

Deep Reinforcement Learning and Simulation as a Path Toward Precision Medicine

Traditionally, precision medicine involves classifying patients to identify subpopulations that respond favorably to specific therapeutics. We pose precision medicine as a dynamic feedback control problem, where treatment administered to a patient is guided by measurements taken during the course of treatment. We consider sepsis, a life-threatening condition in which dysregulation of the immune system causes tissue damage. We leverage an existing simulation of the innate immune response to infection and apply deep reinforcement learning (DRL) to discover an adaptive personalized treatment policy that specifies effective multicytokine therapy to simulated sepsis patients based on systemic measurements. The learned policy achieves a dramatic reduction in mortality rate over a set of 500 simulated patients relative to standalone antibiotic therapy. Advantages of our approach are threefold: (1) the use of simulation allows exploring therapeutic strategies beyond clinical practice and available data, (2) advances in DRL accommodate learning complex therapeutic strategies for complex biological systems, and (3) optimized treatments respond to a patient's individual disease progression over time, therefore, capturing both differences across patients and the inherent randomness of disease progression within a single patient. We hope that this work motivates both considering adaptive personalized multicytokine mediation therapy for sepsis and exploiting simulation with DRL for precision medicine more broadly.

Causes and Consequences of Spatial Within-Host Viral Spread

The spread of viral pathogens both between and within hosts is inherently a spatial process. While the spatial aspects of viral spread at the epidemiological level have been increasingly well characterized, the spatial aspects of viral spread within infected hosts are still understudied. Here, with a focus on influenza A viruses (IAVs), we first review experimental studies that have shed light on the mechanisms and spatial dynamics of viral spread within hosts. These studies provide strong empirical evidence for highly localized IAV spread within hosts. Since mathematical and computational within-host models have been increasingly used to gain a quantitative understanding of observed viral dynamic patterns, we then review the (relatively few) computational modeling studies that have shed light on possible factors that structure the dynamics of spatial within-host IAV spread. These factors include the dispersal distance of virions, the localization of the immune response, and heterogeneity in host cell phenotypes across the respiratory tract. While informative, we find in these studies a striking absence of theoretical expectations of how spatial dynamics may impact the dynamics of viral populations. To mitigate this, we turn to the extensive ecological and evolutionary literature on range expansions to provide informed theoretical expectations. We find that factors such as the type of density dependence, the frequency of long-distance dispersal, specific life history characteristics, and the extent of spatial heterogeneity are critical factors affecting the speed of population spread and the genetic composition of spatially expanding populations. For each factor that we identified in the theoretical literature, we draw parallels to its analog in viral populations. We end by discussing current knowledge gaps related to the spatial component of within-host IAV spread and the potential for within-host spatial considerations to inform the development of disease control strategies.

Physiological factors leading to a successful vaccination: A computational approach

The immune system mounts a response to an infection by activating T cells. T cell activation occurs when dendritic cells, which have already interacted with the pathogen, scan a T cell that is cognate for (responsive to) the pathogen. This often occurs inside lymph nodes. The time it takes for this scanning event to occur, indeed the probability that it will occur at all, depends on many factors, including the rate that T cells and dendritic cells enter and leave the lymph node as well as the geometry of the lymph node and of course other cellular and molecular parameters. In this paper, we develop a hybrid stochastic-deterministic mathematical model at the tissue scale of the lymph node and simulate dendritic cells and cognate T cells to investigate the most important physiological factors leading to a successful and timely immune response after a vaccination. We use an agent-based model to describe the small population of cognate naive T cells and a partial differential equation description for the concentration of mature dendritic cells. We estimate the model parameters based on the known literature and measurements previously taken in our lab. We perform a parameter sensitivity analysis to quantify the sensitivity of the model results to the parameters. The results show that increasing T cell inflow through high endothelial venules, restricting cellular egress via the efferent lymph and increasing the total dendritic cell count by improving vaccinations are the among the most important physiological factors leading to an improved immune response. We also find that increasing the physical size of lymph nodes improves the overall likelihood that an immune response will take place but has a fairly weak effect on the response rate. The nature of dendritic cell trafficking through the LN (either passive or active transport) seems to have little effect on the overall immune response except if a change in overall egress time is observed.

The effect of delay in viral production in within-host models during early infection

Delay in viral production may have a significant impact on the early stages of infection. During the eclipse phase, the time from viral entry until active production of viral particles, no viruses are produced. This delay affects the probability that a viral infection becomes established and timing of the peak viral load. Deterministic and stochastic models are formulated with either multiple latent stages or a fixed delay for the eclipse phase. The deterministic model with multiple latent stages approaches in the limit the model with a fixed delay as the number of stages approaches infinity. The deterministic model framework is used to formulate continuous-time Markov chain and stochastic differential equation models. The probability of a minor infection with rapid viral clearance as opposed to a major full-blown infection with a high viral load is estimated from a branching process approximation of the Markov chain model and the results are confirmed through numerical simulations. In addition, parameter values for influenza A are used to numerically estimate the time to peak viral infection and peak viral load for the deterministic and stochastic models. Although the average length of the eclipse phase is the same in each of the models, as the number of latent stages increases, the numerical results show that the time to viral peak and the peak viral load increase.

Modeling epidemics: A primer and Numerus Model Builder implementation

Epidemiological models are dominated by compartmental models, of which SIR formulations are the most commonly used. These formulations can be continuous or discrete (in either the state-variable values or time), deterministic or stochastic, or spatially homogeneous or heterogeneous, the latter often embracing a network formulation. Here we review the continuous and discrete deterministic and discrete stochastic formulations of the SIR dynamical systems models, and we outline how they can be easily and rapidly constructed using Numerus Model Builder, a graphically-driven coding platform. We also demonstrate how to extend these models to a metapopulation setting using NMB network and mapping tools.

Mathematical Analysis of Viral Replication Dynamics and Antiviral Treatment Strategies: From Basic Models to Age-Based Multi-Scale Modeling

Viral infectious diseases are a global health concern, as is evident by recent outbreaks of the middle east respiratory syndrome, Ebola virus disease, and re-emerging zika, dengue, and chikungunya fevers. Viral epidemics are a socio-economic burden that causes short- and long-term costs for disease diagnosis and treatment as well as a loss in productivity by absenteeism. These outbreaks and their socio-economic costs underline the necessity for a precise analysis of virus-host interactions, which would help to understand disease mechanisms and to develop therapeutic interventions. The combination of quantitative measurements and dynamic mathematical modeling has increased our understanding of the within-host infection dynamics and has led to important insights into viral pathogenesis, transmission, and disease progression. Furthermore, virus-host models helped to identify drug targets, to predict the treatment duration to achieve cure, and to reduce treatment costs. In this article, we review important achievements made by mathematical modeling of viral kinetics on the extracellular, intracellular, and multi-scale level for Human Immunodeficiency Virus, Hepatitis C Virus, Influenza A Virus, Ebola Virus, Dengue Virus, and Zika Virus. Herein, we focus on basic mathematical models on the population scale (so-called target cell-limited models), detailed models regarding the most important steps in the viral life cycle, and the combination of both. For this purpose, we review how mathematical modeling of viral dynamics helped to understand the virus-host interactions and disease progression or clearance. Additionally, we review different types and effects of therapeutic strategies and how mathematical modeling has been used to predict new treatment regimens.

A Novel Stochastic Multi-Scale Model of Francisella tularensis Infection to Predict Risk of Infection in a Laboratory

We present a multi-scale model of the within-phagocyte, within-host and population-level infection dynamics of Francisella tularensis, which extends the mechanistic one proposed by Wood et al. (2014). Our multi-scale model incorporates key aspects of the interaction between host phagocytes and extracellular bacteria, accounts for inter-phagocyte variability in the number of bacteria released upon phagocyte rupture, and allows one to compute the probability of response, and mean time until response, of an infected individual as a function of the initial infection dose. A Bayesian approach is applied to parameterize both the within-phagocyte and within-host models using infection data. Finally, we show how dose response probabilities at the individual level can be used to estimate the airborne propagation of Francisella tularensis in indoor settings (such as a microbiology laboratory) at the population level, by means of a deterministic zonal ventilation model.

Cell studio: A platform for interactive, 3D graphical simulation of immunological processes

The field of computer modeling and simulation of biological systems is rapidly advancing, backed by significant progress in the fields of experimentation techniques, computer hardware, and programming software. The result of a simulation may be delivered in several ways, from numerical results, through graphs of the simulated run, to a visualization of the simulation. The vision of an in-silico experiment mimicking an in-vitro or in-vivo experiment as it is viewed under a microscope is appealing but technically demanding and computationally intensive. Here, we report “Cell Studio,” a generic, hybrid platform to simulate an immune microenvironment with biological and biophysical rules. We use game engines—generic programs for game creation which offer ready-made assets and tools—to create a visualized, interactive 3D simulation. We also utilize a scalable architecture that delegates the computational load to a server. The user may view the simulation, move the “camera” around, stop, fast-forward, and rewind it and inject soluble molecules into the extracellular medium at any point in time. During simulation, graphs are created in real time for a broad view of system-wide processes. The model is parametrized using a user-friendly Graphical User Interface (GUI). We show a simple validation simulation and compare its results with those from a “classical” simulation, validated against a “wet” experiment. We believe that interactive, real-time 3D visualization may aid in generating insights from the model and encourage intuition about the immunological scenario.

Computational modeling of immune system of the fish for a more effective vaccination in aquaculture

Motivation

A computational model equipped with the main immunological features of the sea bass (Dicentrarchus labrax L.) immune system was used to predict more effective vaccination in fish. The performance of the model was evaluated by using the results of two in vivo vaccinations trials against L. anguillarum and P. damselae.

Results

Tests were performed to select the appropriate doses of vaccine and infectious bacteria to set up the model. Simulation outputs were compared with the specific antibody production and the expression of BcR and TcR gene transcripts in spleen. The model has shown a good ability to be used in sea bass and could be implemented for different routes of vaccine administration even with more than two pathogens. The model confirms the suitability of in silico methods to optimize vaccine doses and the immune response to them. This model could be applied to other species to optimize the design of new vaccination treatments of fish in aquaculture.

Availability and implementation

The method is available at http://www.iac.cnr.it/∼filippo/c-immsim/.

Contact

nromano@unitus.it.

Supplementary information

Supplementary data are available at Bioinformatics online.

Kinetic Modeling of Virus Growth in Cells

When a virus infects a host cell, it hijacks the biosynthetic capacity of the cell to produce virus progeny, a process that may take less than an hour or more than a week. The overall time required for a virus to reproduce depends collectively on the rates of multiple steps in the infection process, including initial binding of the virus particle to the surface of the cell, virus internalization and release of the viral genome within the cell, decoding of the genome to make viral proteins, replication of the genome, assembly of progeny virus particles, and release of these particles into the extracellular environment. For a large number of virus types, much has been learned about the molecular mechanisms and rates of the various steps. However, in only relatively few cases during the last 50 years has an attempt been made-using mathematical modeling-to account for how the different steps contribute to the overall timing and productivity of the infection cycle in a cell. Here we review the initial case studies, which include studies of the one-step growth behavior of viruses that infect bacteria (Qβ, T7, and M13), human immunodeficiency virus, influenza A virus, poliovirus, vesicular stomatitis virus, baculovirus, hepatitis B and C viruses, and herpes simplex virus. Further, we consider how such models enable one to explore how cellular resources are utilized and how antiviral strategies might be designed to resist escape. Finally, we highlight challenges and opportunities at the frontiers of cell-level modeling of virus infections.