skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Mathematical Modeling of RNA Virus Sensing Pathways Reveals Paracrine Signaling as the Primary Factor Regulating Excessive Cytokine Production
Title: Mathematical Modeling of RNA Virus Sensing Pathways Reveals Paracrine Signaling as the Primary Factor Regulating Excessive Cytokine Production
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
Award ID(s):
1943777
NSF-PAR ID:
10215875
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Processes
Volume:
8
Issue:
6
ISSN:
2227-9717
Page Range / eLocation ID:
719
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; ( , PLOS Computational Biology)
    Kosakovsky Pond, Sergei L. (Ed.)
    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
  2. ; ; ; ; ; ; ; ; ; ; et al ( , Journal of Virology)
    Viral Hemorrhagic Septicemia virus (VHSV) is a pathogenic fish rhabdovirus found in discrete locales throughout the northern hemisphere. VHSV infection of fish cells leads to upregulation of the host's virus detection response, but the virus quickly suppresses interferon (IFN) production and antiviral genes expression. By systematically screening each of the six VHSV structural and nonstructural genes, we have identified matrix protein (M) as its most potent anti-host protein. VHSV-IVb M alone suppressed mitochondrial antiviral signaling protein (MAVS) and type I IFN-induced gene expression in a dose-dependent manner. M also suppressed the constitutively active SV40 promoter and globally decreased cellular RNA levels. Chromatin immunoprecipitation (ChIP) studies illustrated that M inhibited RNA polymerase II (RNAP II) recruitment to gene promoters, and decreased RNAP II CTD Ser2 phosphorylation during VHSV infection. However, transcription directed by RNAP I-III was suppressed by M. To identify regions of functional importance, M proteins from a variety of VHSV strains were tested in cell-based transcriptional inhibition assays. M protein of a particular VHSV-Ia strain, F1, was significantly less potent than -IVb M at inhibiting SV40/luc expression, yet differed by just four amino acids. Mutation of D62 to alanine alone, or in combination with an E181 to alanine mutation (D62A/E181A), dramatically reduced the ability of -IVb M to suppress host transcription. Introducing either M D62A or D62A/E181A mutations into VHSV-IVb via reverse genetics resulted in viruses that replicated efficiently but exhibited less cytotoxicity and reduced anti-transcriptional activities, implicating M as a primary regulator of cytopathicity and host transcriptional suppression. Importance: Viruses must suppress host antiviral responses to replicate and spread between hosts. In these studies, we identified the matrix protein of the deadly fish Novirhabdovirus, VHSV, as a critical mediator of host suppression during infection. Our studies indicated that M alone could block cellular gene expression at very low expression levels. We identified several subtle mutations in M that were less potent at suppressing host transcription. When these mutations were engineered back into recombinant viruses, the resulting viruses replicated well but elicited less toxicity in infected cells and activated host innate immune responses more robustly. These data demonstrated that VHSV M plays an important role in mediating both virus-induced cell toxicity and viral replication. Our data suggest that its roles in these two processes can be separated to design effective attenuated viruses for vaccine candidates. 
    more » « less
  3. ; ; ; ( , Seminars in Cell & Developmental Biology)
    Ubiquitin-like proteins (Ubls) share some features with ubiquitin (Ub) such as their globular 3D structure and the ability to attach covalently to other proteins. Interferon Stimulated Gene 15 (ISG15) is an abundant Ubl that similar to Ub, marks many hundreds of cellular proteins, altering their fate. In contrast to Ub, , ISG15 requires interferon (IFN) induction to conjugate efficiently to other proteins. Moreover, despite the multitude of E3 ligases for Ub-modified targets, a single E3 ligase termed HERC5 (in humans) is responsible for the bulk of ISG15 conjugation. Targets include both viral and cellular proteins spanning an array of cellular compartments and metabolic pathways. So far, no common structural or biochemical feature has been attributed to these diverse substrates, raising questions about how and why they are selected. Conjugation of ISG15 mitigates some viral and bacterial infections and is linked to a lower viral load pointing to the role of ISG15 in the cellular immune response. In an apparent attempt to evade the immune response, some viruses try to interfere with the ISG15 pathway. For example, deconjugation of ISG15 appears to be an approach taken by coronaviruses to interfere with ISG15 conjugates. Specifically, coronaviruses such as SARS-CoV, MERS-CoV, and SARS-CoV-2, encode papain-like proteases (PL1pro) that bear striking structural and catalytic similarities to the catalytic core domain of eukaryotic deubiquitinating enzymes of the Ubiquitin-Specific Protease (USP) sub-family. The cleavage specificity of these PLpro enzymes is for flexible polypeptides containing a consensus sequence (R/K)LXGG, enabling them to function on two seemingly unrelated categories of substrates: (i) the viral polyprotein 1 (PP1a, PP1ab) and (ii) Ub- or ISG15-conjugates. As a result, PLpro enzymes process the viral polyprotein 1 into an array of functional proteins for viral replication (termed non-structural proteins; NSPs), and it can remove Ub or ISG15 units from conjugates. However, by de-conjugating ISG15, the virus also creates free ISG15, which in turn may affect the immune response in two opposite pathways: free ISG15 negatively regulates IFN signaling in humans by binding non-catalytically to USP18, yet at the same time free ISG15 can be secreted from the cell and induce the IFN pathway of the neighboring cells. A deeper understanding of this protein-modification pathway and the mechanisms of the enzymes that counteract it will bring about effective clinical strategies related to viral and bacterial infections 
    more » « less
  4. ; ; ( , Viruses)
    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
  5. ; ; ( , mSphere)
    Lee, Benhur (Ed.)
    ABSTRACT Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has infected over 40 million people worldwide, with over 1 million deaths as of October 2020 and with multiple efforts in the development and testing of antiviral drugs and vaccines under way. In order to gain insights into SARS-CoV-2 evolution and drug targets, we investigated how and to what extent the SARS-CoV-2 genome sequence differs from those of other well-characterized human and animal coronavirus genomes, as well as how polymorphic SARS-CoV-2 genomes are generally. We ultimately sought to identify features in the SARS-CoV-2 genome that may contribute to its viral replication, host pathogenicity, and vulnerabilities. Our analyses suggest the presence of unique sequence signatures in the 3′ untranslated region (3′-UTR) of betacoronavirus lineage B, which phylogenetically encompasses SARS-CoV-2 and SARS-CoV as well as multiple groups of bat and animal coronaviruses. In addition, we identified genome-wide patterns of variation across different SARS-CoV-2 strains that likely reflect the effects of selection. Finally, we provide evidence for a possible host-microRNA-mediated interaction between the 3′-UTR and human microRNA hsa-miR-1307-3p based on the results of multiple computational target prediction analyses and an assessment of similar interactions involving the influenza A H1N1 virus. This interaction also suggests a possible survival mechanism, whereby a mutation in the SARS-CoV-2 3′-UTR leads to a weakened host immune response. The potential roles of host microRNAs in SARS-CoV-2 replication and infection and the exploitation of conserved features in the 3′-UTR as therapeutic targets warrant further investigation. IMPORTANCE The coronavirus disease 2019 (COVID-19) outbreak is having a dramatic global effect on public health and the economy. As of October 2020, SARS-CoV-2 has been detected in over 189 countries, has infected over 40 million people, and is responsible for more than 1 million deaths. The genome of SARS-CoV-2 is small but complex, and its functions and interactions with human host factors are being studied extensively. The significance of our study is that, using extensive SARS-CoV-2 genome analysis techniques, we identified potential interacting human host microRNA targets that share similarity with those of influenza A virus H1N1. Our study results will allow the development of virus-host interaction models that will enhance our understanding of SARS-CoV-2 pathogenesis and motivate the exploitation of both the interacting viral and host factors as therapeutic targets. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email