skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Multicellular spatial model of RNA virus replication and interferon responses reveals factors controlling plaque growth dynamics
Respiratory viruses present major public health challenges, as evidenced by the 1918 Spanish Flu, the 1957 H2N2, 1968 H3N2, and 2009 H1N1 influenza pandemics, and the ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Severe RNA virus respiratory infections often correlate with high viral load and excessive inflammation. Understanding the dynamics of the innate immune response and its manifestations at the cell and tissue levels is vital to understanding the mechanisms of immunopathology and to developing strain-independent treatments. Here, we present a novel spatialized multicellular computational model of RNA virus infection and the type-I interferon-mediated antiviral response that it induces within lung epithelial cells. The model is built using the CompuCell3D multicellular simulation environment and is parameterized using data from influenza virus-infected cell cultures. Consistent with experimental observations, it exhibits either linear radial growth of viral plaques or arrested plaque growth depending on the local concentration of type I interferons. The model suggests that modifying the activity of signaling molecules in the JAK/STAT pathway or altering the ratio of the diffusion lengths of interferon and virus in the cell culture could lead to plaque growth arrest. The dependence of plaque growth arrest on diffusion lengths highlights the importance of developing validated spatial models of cytokine signaling and the need for in vitro measurement of these diffusion coefficients. Sensitivity analyses under conditions leading to continuous or arrested plaque growth found that plaque growth is more sensitive to variations of most parameters and more likely to have identifiable model parameters when conditions lead to plaque arrest. This result suggests that cytokine assay measurements may be most informative under conditions leading to arrested plaque growth. The model is easy to extend to include SARS-CoV-2-specific mechanisms or to use as a component in models linking epithelial cell signaling to systemic immune models.  more » « less
Award ID(s):
1943777 1720625
PAR ID:
10316925
Author(s) / Creator(s):
; ; ; ;
Editor(s):
Kosakovsky Pond, Sergei L.
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
17
Issue:
10
ISSN:
1553-7358
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    RNA viruses, such as influenza and Severe Acute Respiratory Syndrome (SARS), invoke excessive immune responses; however, the kinetics that regulate inflammatory responses within infected cells remain unresolved. Here, we develop a mathematical model of the RNA virus sensing pathways, to determine the intracellular events that primarily regulate interferon, an important protein for the activation and management of inflammation. Within the ordinary differential equation (ODE) model, we incorporate viral replication, cell death, interferon stimulated genes’ antagonistic effects on viral replication, and virus sensor protein (TLR and RIG-I) kinetics. The model is parameterized to influenza infection data using Markov chain Monte Carlo and then validated against infection data from an NS1 knockout strain of influenza, demonstrating that RIG-I antagonism significantly alters cytokine signaling trajectory. Global sensitivity analysis suggests that paracrine signaling is responsible for the majority of cytokine production, suggesting that rapid cytokine production may be best managed by influencing extracellular cytokine levels. As most of the model kinetics are host cell specific and not virus specific, the model presented provides an important step to modeling the intracellular immune dynamics of many RNA viruses, including the viruses responsible for SARS, Middle East Respiratory Syndrome (MERS), and Coronavirus Disease (COVID-19). 
    more » « less
  2. The timing and magnitude of the immune response (i.e., the immunodynamics) associated with the early innate immune response to viral infection display distinct trends across influenza A virus subtypes in vivo. Evidence shows that the timing of the type-I interferon response and the overall magnitude of immune cell infiltration are both correlated with more severe outcomes. However, the mechanisms driving the distinct immunodynamics between infections of different virus strains (strain-specific immunodynamics) remain unclear. Here, computational modeling and strain-specific immunologic data are used to identify the immune interactions that differ in mice infected with low-pathogenic H1N1 or high-pathogenic H5N1 influenza viruses. Computational exploration of free parameters between strains suggests that the production rate of interferon is the major driver of strain-specific immune responses observed in vivo, and points towards the relationship between the viral load and lung epithelial interferon production as the main source of variance between infection outcomes. A greater understanding of the contributors to strain-specific immunodynamics can be utilized in future efforts aimed at treatment development to improve clinical outcomes of high-pathogenic viral strains. 
    more » « less
  3. Smith, Amber M (Ed.)
    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. 
    more » « less
  4. At the end of 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a novel human coronavirus, emerged and rapidly caused a global pandemic. SARS-CoV-2 is the causative agent of coronavirus disease 2019 (COVID-19), which affects the respiratory tract and lungs of infected individuals. Due to the increased transmissibility of the SARS-CoV-2 virus compared to its previous versions, determining as fully as possible the various structural aspects of the virus became critical for the development of therapeutics and vaccines to combat this virus. Knowing the structures of viral proteins and their glycosylation is an essential foundation for the understanding of the mechanism of the disease. Glycopeptide analysis has been used to map the glycosylation of viral glycoproteins, including those of influenza and HIV. Thanks to the developments in the field over the last few decades, scientists were able to quickly develop therapeutics against SARS-CoV-2. This chapter discusses the four structural proteins of SARS-CoV-2, their glycosylation and modifications, and the techniques used to map SARS-CoV-2 glycosylation. 
    more » « less
  5. The emergence and mutation of pathogenic viruses have been occurring at an unprecedented rate in recent decades. The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has developed into a global public health crisis due to extensive viral transmission. In situ RNA mapping has revealed angiotensin-converting enzyme 2 (ACE2) expression to be highest in the nose and lower in the lung, pointing to nasal susceptibility as a predominant route for infection and the cause of subsequent pulmonary effects. By blocking viral attachment and entry at the nasal airway using a cyclodextrin-based formulation, a preventative therapy can be developed to reduce viral infection at the site of entry. Here, we assess the safety and antiviral efficacy of cyclodextrin-based formulations. From these studies, hydroxypropyl beta-cyclodextrin (HPBCD) and hydroxypropyl gamma-cyclodextrin (HPGCD) were then further evaluated for antiviral effects using SARS-CoV-2 pseudotypes. Efficacy findings were confirmed with SARS-CoV-2 Delta variant infection of Calu-3 cells and using a K18-hACE2 murine model. Intranasal pre-treatment with HPBCD-based formulations reduced viral load and inflammatory signaling in the lung. In vitro efficacy studies were further conducted using lentiviruses, murine hepatitis virus (MHV), and influenza A virus subtype H1N1. These findings suggest HPBCD may be used as an agnostic barrier against transmissible pathogens, including but not limited to SARS-CoV-2. 
    more » « less