Bioinformatics Modeling of Host Interferon-Stimulated Gene (ISG) Restriction Factors

By Dr. Zubair Khalid, DVM, MS, PhD · 2026-05-10

Key Takeaways

• Bioinformatics approaches, including transcriptomic screening (e.g., scRNA-seq) and protein-protein interaction networks, are crucial for identifying host Interferon-Stimulated Genes (ISGs) that function as antiviral restriction factors. These methods have revealed cell-type-specific ISG expression patterns correlating with viral susceptibility, such as IRF1 restricting Hepatitis Delta Virus (HDV) in hepatocytes.
• Structural modeling techniques, such as homology modeling, molecular docking, and molecular dynamics simulations, are essential for elucidating the molecular interfaces between restriction factors and viral targets. This allows for the mapping of critical contact residues and the understanding of mechanisms like CCND3 blocking SFTSV nucleoprotein multimerization.
• Viruses have evolved sophisticated escape mechanisms, including direct antagonism of restriction factors (e.g., SARS-CoV-2 ORF6), degradation via ubiquitin-proteasome pathways, and mutations within viral proteins that disrupt binding to host restriction factors. Bioinformatics can predict these escape variants by analyzing interface residues.
• Computational prediction models, including machine learning and protein language models, are increasingly used to identify novel ISGs with antiviral activity based on sequence features and structural homology. These models aid in prioritizing candidates for functional validation and understanding broad-spectrum antiviral activity, such as ZAP's RNA binding and exosome recruitment.
• Comparative transcriptomics across species is instrumental in identifying species-specific restriction barriers and candidate ISGs, as demonstrated by the black flying fox RTP4 orthologue inhibiting Hepatitis C Virus (HCV) replication in mouse cells. This approach is vital for understanding host range and developing veterinary antiviral strategies.

Introduction

Interferon-stimulated genes (ISGs) constitute the primary effector arm of the host innate antiviral response. Upon interferon induction, hundreds of ISG products collectively target viral replication at multiple stages, including entry, transcription, translation, assembly, and egress [<a href="#ref-1">1</a>, <a href="#ref-2">2</a>]. The characterization of these antiviral proteins, collectively termed restriction factors, has been advanced substantially by bioinformatics screening, structural modeling, and interactomics [<a href="#ref-3">3</a>, <a href="#ref-4">4</a>]. In veterinary virology, understanding ISG restriction mechanisms is critical for elucidating host range barriers, viral pathogenesis, and for designing host-directed antiviral strategies [<a href="#ref-5">5</a>, <a href="#ref-6">6</a>].

This review provides a systematic examination of bioinformatics methodologies used to model host ISG restriction factors. We focus on the molecular interfaces between restriction factors and viral targets, the structural basis of antiretroviral activity, and the viral escape mutations that circumvent these defenses. The discussion is contextualized within veterinary species, with comparative references to human pathogens where relevant for host-range parallels.

Bioinformatics Approaches for ISG Discovery and Functional Annotation

The identification of ISG restriction factors has been driven by transcriptomic screening, loss-of-function CRISPR libraries, and protein interaction networks [<a href="#ref-1">1</a>, <a href="#ref-3">3</a>]. Single-cell RNA sequencing (scRNA-seq) has revealed cell-type-specific ISG expression patterns that correlate with susceptibility to viral infection [<a href="#ref-3">3</a>]. For example, analysis of hepatocyte-like cells identified interferon regulatory factor 1 (IRF1) as a potent restriction factor for hepatitis delta virus (HDV) through a mechanism blocking early viral replication [<a href="#ref-3">3</a>]. Similarly, transcriptomic profiling of porcine epidemic diarrhea virus (PEDV) variant strains demonstrated that enhanced immune responses upregulate ISGs such as IFI44 and OASL, which restrict viral replication through positive regulation of type I interferon signaling [<a href="#ref-6">6</a>].

Comparative transcriptomics across species has been instrumental in identifying candidate restriction factors. The black flying fox orthologue of receptor transporter protein 4 (RTP4) was shown to inhibit hepatitis C virus (HCV) replication in mouse cells, illustrating how cross-species comparisons can reveal species-specific restriction barriers [<a href="#ref-7">7</a>]. In veterinary contexts, similar approaches applied to porcine, avian, and equine cells have identified ISGs that limit coronaviruses, influenza A viruses, and arteriviruses [<a href="#ref-6">6</a>, <a href="#ref-8">8</a>].

A key computational pipeline for ISG discovery involves integrating transcriptomic data with protein-protein interaction (PPI) networks. The work by Zhao et al. [<a href="#ref-1">1</a>] screened over 200 ISGs for activity against bandaviruses using a library-based overexpression approach in cell culture, followed by co-immunoprecipitation and structural modeling to map the CCND3-nucleoprotein interaction interface. Such integrative frameworks are now standard for prioritizing ISG candidates for mechanistic validation.

flowchart TD
    A["Interferon induction"] --> B["Transcriptomic profiling (RNA-seq, scRNA-seq)"]
    B --> C["ISG library screening / CRISPR knockout"]
    C --> D["Validation of antiviral activity"]
    D --> E["Identification of viral target"]
    E --> F["Co-IP / proximity ligation to map interaction"]
    F --> G["Structural modeling (AlphaFold2, MD simulations)"]
    G --> H["Mutation analysis of interface residues"]
    H --> I["Viral escape variant identification"]
    I --> J["Functional assays in relevant host cells"]
    J --> D

Structural Modeling of ISG-Virus Interfaces

High-resolution structural characterization of restriction factor-viral protein complexes is essential for understanding the molecular basis of antiviral activity. Computational approaches, including homology modeling, molecular docking, and molecular dynamics (MD) simulations, have been deployed to predict interaction interfaces and identify contact residues [<a href="#ref-1">1</a>, <a href="#ref-9">9</a>].

One paradigmatic example is the interaction between cyclin D3 (CCND3) and the nucleoprotein (NP) of severe fever with thrombocytopenia syndrome virus (SFTSV), a high-pathogenic bandavirus [<a href="#ref-1">1</a>]. Using structural modeling, the CN domain of CCND3 was shown to bind the NP “head” region in an RNA-independent manner, blocking NP multimerization, NP-RNA binding, and NP association with the viral polymerase. The interface residues mapped from this analysis are critical for viral ribonucleoprotein (RNP) function [<a href="#ref-1">1</a>]. Such modeling not only reveals the mechanism of restriction but also predicts residues that, when mutated, may allow escape.

For the interferon-induced transmembrane (IFITM) family, lymphocyte antigen 6 complex locus E (LY6E) was shown to restrict coronavirus entry by interfering with spike protein-mediated membrane fusion [<a href="#ref-9">9</a>]. Structural and mutagenesis studies have identified that LY6E acts at the fusion step, and computational docking simulations have suggested that LY6E interacts with the heptad repeat regions of the spike protein, though high-resolution structures remain limited [<a href="#ref-9">9</a>]. Similarly, the guanylate-binding protein (GBP) family, particularly human GBP1, inhibits SARS-CoV-2 replication through a mechanism distinct from spike processing inhibition, as demonstrated by farrukee et al. [<a href="#ref-2">2</a>]. Structural modeling of GBP1 in complex with viral nonstructural proteins is ongoing.

The TRIM family of E3 ubiquitin ligases includes several ISG restriction factors. TRIM22 is constitutively expressed in respiratory epithelium and restricts influenza A virus (IAV) transcription [<a href="#ref-8">8</a>]. The RING domain of TRIM22 mediates ubiquitination of the viral nucleoprotein, and computational modeling of the TRIM22-NP interface has identified key lysine residues required for degradation [<a href="#ref-8">8</a>]. Tetherin (BST2) is another well-characterized ISG that physically traps enveloped virions at the cell surface. In the context of herpes simplex virus type 1 (HSV-1) corneal infection, tetherin upregulation via STING-dependent signaling is critical for containing viral spread [<a href="#ref-10">10</a>]. Structural studies of tetherin have revealed a coiled-coil domain that spans the viral envelope, and the specific residues involved in virion tethering have been mapped by mutagenesis and cryo-electron microscopy.

Viral Escape Mechanisms and Mutation of Restriction Interfaces

Viruses have evolved multiple strategies to evade ISG restriction. These include direct antagonism of the restriction factor, degradation via the ubiquitin-proteasome pathway, and mutational alteration of the interface recognized by the host factor [<a href="#ref-1">1</a>, <a href="#ref-11">11</a>, <a href="#ref-12">12</a>]. For example, the SARS-CoV-2 ORF6 protein is a potent interferon antagonist that inhibits ISG expression. The host CRL4B E3 ligase recruited by PRPF19 targets ORF6 for ubiquitin-dependent degradation, thereby relieving interferon inhibition [<a href="#ref-11">11</a>]. Viral escape can occur through mutations in ORF6 that reduce its recognition by the ligase complex.

The nonstructural protein NSs of bandaviruses partially antagonizes CCND3 by attenuating its induction and promoting autophagic degradation [<a href="#ref-1">1</a>]. This antagonism represents a direct viral countermeasure against an ISG restriction factor. Similarly, the Vif protein of lentiviruses (not discussed in provided papers) degrades APOBEC3 proteins, but in the context of hepatitis B virus (HBV), the X protein (HBx) downregulates tree shrew APOBEC3 (tsAPOBEC3) via exosomes, although tsAPOBEC3 retains partial antiviral activity [<a href="#ref-13">13</a>].

Mutations in the viral interface that disrupt binding to restriction factors are common. The RNA-dependent RNA polymerase 3Dpol of enterovirus A71 is targeted by the ISG SHFL (shiftless antiviral protein) for ubiquitin-proteasome degradation. The zinc-finger domain and residues 164-199 of SHFL are critical for interaction with 3Dpol, and mutations in 3Dpol that alter this interface confer resistance [<a href="#ref-12">12</a>]. Bioinformatics modeling of these mutations can predict escape variants before they emerge in the field.

The mermaid diagram above illustrates the iterative workflow from ISG discovery to escape variant prediction.

Computational Prediction of Restriction Factor Activity

Machine learning and deep learning approaches have been applied to predict which ISGs are likely to restrict specific viruses based on sequence features, expression patterns, and structural homology. For example, protein language models trained on interferon-induced transcriptomes can identify novel ISGs with high precision [<a href="#ref-14">14</a>]. The zinc finger antiviral protein (ZAP) was initially identified as a retrovirus restriction factor and later shown to inhibit LINE-1 retroelements [<a href="#ref-14">14</a>]. ZAP binds viral RNA via its zinc fingers and recruits the exosome for degradation. Computational models of ZAP-RNA interaction have identified the RNA sequence motifs that confer binding specificity, although ZAP shows broad spectrum activity, suggesting recognition of structural features rather than primary sequence [<a href="#ref-14">14</a>].

The spliceosome factor SART1 regulates ISG expression through alternative splicing, as demonstrated by Lin et al. [<a href="#ref-15">15</a>]. In silico prediction of splicing events influenced by SART1 has identified novel antiviral isoforms of EIF4G3 and GORASP2 that restrict HCV. This highlights that restriction factor activity can be regulated at the post-transcriptional level, and computational modeling must account for splicing variants.

An important computational challenge is the prediction of interface contact residues from sequence alone. The CCND3-NP interface [<a href="#ref-1">1</a>] and the SHFL-3Dpol interface [<a href="#ref-12">12</a>] were initially mapped using co-immunoprecipitation and then refined with AlphaFold2-based docking. Molecular dynamics simulations of the complexes can predict binding free energies and the impact of mutations on affinity [<a href="#ref-9">9</a>, <a href="#ref-12">12</a>].

Tables of Key ISG Restriction Factors and Their Mechanisms

Conclusion

Bioinformatics modeling of host ISG restriction factors has become an indispensable tool for veterinary virology. By integrating transcriptomic data, protein interaction networks, and structural modeling, researchers can identify candidate restriction factors, map their interfaces with viral targets, and predict escape mutations. The examples discussed highlight the diversity of restriction mechanisms, from direct binding and degradation (CCND3, SHFL) to membrane fusion inhibition (LY6E) and RNA degradation (ZAP). The computational prediction of these interactions informs the development of host-directed therapeutics and the assessment of cross-species transmission risk. Future work will likely incorporate advanced artificial intelligence methods to predict interface dynamics and co-evolutionary constraints, further refining our understanding of the ongoing arms race between hosts and viruses.

Related Clinical & Scientific Guides

• A Practical Guide to Detecting Antimicrobial Resistance Genes in Shotgun Metagenomic Data
• Computational Immunology: Modeling the Immune System
• How to Set Hard Filters for Germline Variant Calling: A Practical Guide to GATK Best Practices