Evolutionary arms race between SARS-CoV-2 and interferon signaling via dynamic interaction with autophagy.

SARS-CoV-2 emerged, and is evolving to efficiently infect humans worldwide. SARS-CoV-2 evades early innate recognition, interferon signaling activated only in bystander cells. This balance of innate activation and viral evasion has important consequences, but the pathways involved are incompletely understood. Here we find that autophagy genes regulate innate immune signaling, impacting the basal set point of interferons, and thus permissivity to infection. Mechanistically, autophagy genes negatively regulate MAVS, and this low basal level of MAVS is efficiently antagonized by SARS-CoV-2 ORF9b, blocking interferon activation in infected cells. However, upon loss of autophagy increased MAVS overcomes ORF9b-mediated antagonism suppressing infection. This has led to the evolution of SARS-CoV-2 variants to express higher levels of ORF9b, allowing SARS-CoV-2 to replicate under conditions of increased MAVS signaling. Altogether, we find a critical role of autophagy in the regulation of innate immunity and uncover an evolutionary trajectory of SARS-CoV-2 ORF9b to overcome host defenses.

Single cell susceptibility to SARS-CoV-2 infection is driven by variable cell states.

The ability of a virus to infect a cell type is at least in part determined by the presence of host factors required for the viral life cycle. However, even within cell types that express known factors needed for infection, not every cell is equally susceptible, suggesting that our knowledge of the full spectrum of factors that promote infection is incomplete. Profiling the most susceptible subsets of cells within a population may reveal additional factors that promote infection. However, because viral infection dramatically alters the state of the cell, new approaches are needed to reveal the state of these cells prior to infection with virus. Here, we used single-cell clone tracing to retrospectively identify and characterize lung epithelial cells that are highly susceptible to infection with SARS-CoV-2. The transcriptional state of these highly susceptible cells includes markers of retinoic acid signaling and epithelial differentiation. Loss of candidate factors identified by our approach revealed that many of these factors play roles in viral entry. Moreover, a subset of these factors exert control over the infectable cell state itself, regulating the expression of key factors associated with viral infection and entry. Analysis of patient samples revealed the heterogeneous expression of these factors across both cells and patients in vivo. Further, the expression of these factors is upregulated in particular inflammatory pathologies. Altogether, our results show that the variable expression of intrinsic cell states is a major determinant of whether a cell can be infected by SARS-CoV-2.

SARS-CoV-2 ORF6 protein does not antagonize interferon signaling in respiratory epithelial Calu-3 cells during infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused millions of deaths, posing a substantial threat to global public health. Viruses evolve different strategies to antagonize or evade host immune responses. While ectopic expression of SARS-CoV-2 accessory protein ORF6 blocks interferon (IFN) production and downstream IFN signaling, the role of ORF6 in IFN signaling during bona fide viral infection of respiratory cells is unclear. By comparing wild-type (WT) and ORF6-deleted (ΔORF6) SARS-CoV-2 infection and IFN signaling in respiratory cells, we found that ΔORF6 SARS-CoV-2 replicates more efficiently than WT virus and, thus, stimulates more robust immune signaling. Loss of ORF6 does not alter innate signaling in infected cells: both WT and ΔORF6 virus induce delayed IFN responses only in bystander cells. Moreover, expression of ORF6 in the context of SARS-CoV-2 infection has no effect on Sendai virus-stimulated IFN induction: robust translocation of IRF3 is observed in both SARS-CoV-2 infected and bystander cells. Furthermore, IFN pretreatment potently blocks WT and ΔORF6 virus replication similarly, and both viruses fail to suppress the induction of interferon-stimulated genes (ISGs) upon IFN-ß treatment. However, upon treatment with IFN-ß, only bystander cells induce STAT1 translocation during infection with WT virus, whereas ΔORF6 virus-infected cells now show translocation. This suggests that under conditions of high IFN activation, ORF6 can attenuate STAT1 activation. These data provide evidence that ORF6 is not sufficient to antagonize IFN production or IFN signaling in SARS-CoV-2-infected respiratory cells but may impact the efficacy of therapeutics that stimulate innate immune pathways. IMPORTANCE Previous studies identified several SARS-CoV-2 proteins, including ORF6, that antagonize host innate immune responses in the context of overexpression of viral proteins in non-respiratory cells. We set out to determine the role of ORF6 in IFN responses during SARS-CoV-2 infection of respiratory cells. Using a deletion strain, we observed no reduction of infection and no difference in evasion of IFN signaling, with responses limited to bystander cells. Moreover, stimulation of Sendai virus-induced IFN production or IFN-ß-stimulated ISG expression was comparable between SARS-CoV-2 virus and SARS-CoV-2 lacking ORF6 virus, suggesting that ORF6 is not sufficient to counteract IFN induction or IFN signaling during viral infection.

Pyrimidine inhibitors synergize with nucleoside analogues to block SARS-CoV-2.

The SARS-CoV-2 virus has infected more than 261 million people and has led to more than 5 million deaths in the past year and a half1 ( https://www.who.org/ ). Individuals with SARS-CoV-2 infection typically develop mild-to-severe flu-like symptoms, whereas infection of a subset of individuals leads to severe-to-fatal clinical outcomes2. Although vaccines have been rapidly developed to combat SARS-CoV-2, there has been a dearth of antiviral therapeutics. There is an urgent need for therapeutics, which has been amplified by the emerging threats of variants that may evade vaccines. Large-scale efforts are underway to identify antiviral drugs. Here we screened approximately 18,000 drugs for antiviral activity using live virus infection in human respiratory cells and validated 122 drugs with antiviral activity and selectivity against SARS-CoV-2. Among these candidates are 16 nucleoside analogues, the largest category of clinically used antivirals. This included the antivirals remdesivir and molnupiravir, which have been approved for use in COVID-19. RNA viruses rely on a high supply of nucleoside triphosphates from the host to efficiently replicate, and we identified a panel of host nucleoside biosynthesis inhibitors as antiviral. Moreover, we found that combining pyrimidine biosynthesis inhibitors with antiviral nucleoside analogues synergistically inhibits SARS-CoV-2 infection in vitro and in vivo against emerging strains of SARS-CoV-2, suggesting a clinical path forward.

Pyrimidine biosynthesis inhibitors synergize with nucleoside analogs to block SARS-CoV-2 infection.

The ongoing COVID-19 pandemic has highlighted the dearth of approved drugs to treat viral infections, with only ∼90 FDA approved drugs against human viral pathogens. To identify drugs that can block SARS-CoV-2 replication, extensive drug screening to repurpose approved drugs is underway. Here, we screened ∼18,000 drugs for antiviral activity using live virus infection in human respiratory cells. Dose-response studies validate 122 drugs with antiviral activity and selectivity against SARS-CoV-2. Amongst these drug candidates are 16 nucleoside analogs, the largest category of clinically used antivirals. This included the antiviral Remdesivir approved for use in COVID-19, and the nucleoside Molnupirivir, which is undergoing clinical trials. RNA viruses rely on a high supply of nucleoside triphosphates from the host to efficiently replicate, and we identified a panel of host nucleoside biosynthesis inhibitors as antiviral, and we found that combining pyrimidine biosynthesis inhibitors with antiviral nucleoside analogs synergistically inhibits SARS-CoV-2 infection in vitro and in vivo suggesting a clinical path forward.

Drug repurposing screens reveal cell-type-specific entry pathways and FDA-approved drugs active against SARS-Cov-2.

There is an urgent need for antivirals to treat the newly emerged severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). To identify new candidates, we screen a repurposing library of ∼3,000 drugs. Screening in Vero cells finds few antivirals, while screening in human Huh7.5 cells validates 23 diverse antiviral drugs. Extending our studies to lung epithelial cells, we find that there are major differences in drug sensitivity and entry pathways used by SARS-CoV-2 in these cells. Entry in lung epithelial Calu-3 cells is pH independent and requires TMPRSS2, while entry in Vero and Huh7.5 cells requires low pH and triggering by acid-dependent endosomal proteases. Moreover, we find nine drugs are antiviral in respiratory cells, seven of which have been used in humans, and three are US Food and Drug Administration (FDA) approved, including cyclosporine. We find that the antiviral activity of cyclosporine is targeting Cyclophilin rather than calcineurin, revealing essential host targets that have the potential for rapid clinical implementation.

Kinetic Characterization of ASXL1/2-Mediated Allosteric Regulation of the BAP1 Deubiquitinase.

BAP1 is an ubiquitin hydrolase whose deubiquitinase activity is mediated by polycomb group-like protein ASXL2. Cancer-related BAP1 mutations/deletions lead to loss-of-function by targeting the catalytic ubiquitin C-terminal hydrolase (UCH) or UCH37-like domain (ULD) domains of BAP1, and the latter disrupts binding to ASXL2, an obligate partner for BAP1 enzymatic activity. However, the biochemical and biophysical properties of domains involved in forming the enzymatically active complex are unknown. Here, we report the molecular dynamics, kinetics, and stoichiometry of these interactions. We demonstrate that interactions between BAP1 and ASXL2 are direct, specific, and stable to biochemical and biophysical manipulations as detected by isothermal titration calorimetry (ITC), GST association, and optical biosensor assays. Association of the ASXL2-AB box greatly stimulates BAP1 activity. A stable ternary complex is formed, comprised of the BAP1-UCH, BAP1-ULD, and ASXL2-AB domains. Stoichiometric analysis revealed that one molecule of the ULD domain directly interacts with one molecule of the AB box. Real-time kinetic analysis of the ULD/AB protein complex to the BAP1-UCH domain, based on surface plasmon resonance, indicated that formation of the ULD/AB complex with the UCH domain is a single-step event with fast association and slow dissociation rates. In vitro experiments validated in cells that the ASXL-AB box directly regulates BAP1 activity. IMPLICATIONS: Collectively, these data elucidate molecular interactions between specific protein domains regulating BAP1 deubiquitinase activity, thus establishing a foundation for small-molecule approaches to reactivate latent wild-type BAP1 catalytic activity in BAP1-mutant cancers.

Cell-cycle-dependent EBNA1-DNA crosslinking promotes replication termination at oriP and viral episome maintenance.

Epstein-Barr virus (EBV) is an oncogenic human herpesvirus that persists as a multicopy episome in proliferating host cells. Episome maintenance is strictly dependent on EBNA1, a sequence-specific DNA-binding protein with no known enzymatic activities. Here, we show that EBNA1 forms a cell cycle-dependent DNA crosslink with the EBV origin of plasmid replication oriP. EBNA1 tyrosine 518 (Y518) is essential for crosslinking to oriP and functionally required for episome maintenance and generation of EBV-transformed lymphoblastoid cell lines (LCLs). Mechanistically, Y518 is required for replication fork termination at oriP in vivo and for formation of SDS-resistant complexes in vitro. EBNA1-DNA crosslinking corresponds to single-strand endonuclease activity specific to DNA structures enriched at replication-termination sites, such as 4-way junctions. These findings reveal that EBNA1 forms tyrosine-dependent DNA-protein crosslinks and single-strand cleavage at oriP required for replication termination and viral episome maintenance.

The LIM protein Ajuba recruits DBC1 and CBP/p300 to acetylate ERα and enhances ERα target gene expression in breast cancer cells.

Estrogen/ERα signaling is critical for breast cancer progression and therapeutic treatments. Thus, identifying new regulators of this pathway will help to develop new therapeutics to overcome chemotherapy resistance of the breast cancer cells. Here, we report Ajuba directly interacts with ERα to potentiate ERα target gene expression, and biologically Ajuba promotes breast cancer cell growth and contributes to tamoxifen resistance of these cells. Ajuba constitutively binds the DBD and AF2 regions of ERα, and these interactions can be markedly enhanced by estrogen treatment. Mechanistically, Ajuba recruits DBC1 and CBP/p300 and forms a ternary complex to co-activate ERα transcriptional activity and concomitantly enhances ERα acetylation. Moreover, components of this complex can be found at endogenous promoters containing functional ERα responsive elements. Taken together, these data demonstrate that Ajuba functions as a novel co-activator of ERα and that Ajuba/DBC1/CBP/p300 ternary complex may be a new target for developing therapeutics to treat breast cancer.

An evolutionarily conserved DNA architecture determines target specificity of the TWIST family bHLH transcription factors.

Basic helix-loop-helix (bHLH) transcription factors recognize the canonical E-box (CANNTG) to regulate gene transcription; however, given the prevalence of E-boxes in a genome, it has been puzzling how individual bHLH proteins selectively recognize E-box sequences on their targets. TWIST is a bHLH transcription factor that promotes epithelial-mesenchymal transition (EMT) during development and tumor metastasis. High-resolution mapping of TWIST occupancy in human and Drosophila genomes reveals that TWIST, but not other bHLH proteins, recognizes a unique double E-box motif with two E-boxes spaced preferentially by 5 nucleotides. Using molecular modeling and binding kinetic analyses, we found that the strict spatial configuration in the double E-box motif aligns two TWIST-E47 dimers on the same face of DNA, thus providing a high-affinity site for a highly stable intramolecular tetramer. Biochemical analyses showed that the WR domain of TWIST dimerizes to mediate tetramer formation, which is functionally required for TWIST-induced EMT. These results uncover a novel mechanism for a bHLH transcription factor to recognize a unique spatial configuration of E-boxes to achieve target specificity. The WR-WR domain interaction uncovered here sets an example of target gene specificity of a bHLH protein being controlled allosterically by a domain outside of the bHLH region.

Pharmacological protection of retinal pigmented epithelial cells by sulindac involves PPAR-α.

The retinal pigmented epithelial (RPE) layer is one of the major ocular tissues affected by oxidative stress and is known to play an important role in the etiology of age-related macular degeneration (AMD), the major cause of blinding in the elderly. In the present study, sulindac, a nonsteroidal antiinflammatory drug (NSAID), was tested for protection against oxidative stress-induced damage in an established RPE cell line (ARPE-19). Besides its established antiinflammatory activity, sulindac has previously been shown to protect cardiac tissue against ischemia/reperfusion damage, although the exact mechanism was not elucidated. As shown here, sulindac can also protect RPE cells from chemical oxidative damage or UV light by initiating a protective mechanism similar to what is observed in ischemic preconditioning (IPC) response. The mechanism of protection appears to be triggered by reactive oxygen species (ROS) and involves known IPC signaling components such as PKG and PKC epsilon in addition to the mitochondrial ATP-sensitive K(+) channel. Sulindac induced iNOS and Hsp70, late-phase IPC markers in the RPE cells. A unique feature of the sulindac protective response is that it involves activation of the peroxisome proliferator-activated receptor alpha (PPAR-α). We have also used low-passage human fetal RPE and polarized primary fetal RPE cells to validate the basic observation that sulindac can protect retinal cells against oxidative stress. These findings indicate a mechanism for preventing oxidative stress in RPE cells and suggest that sulindac could be used therapeutically for slowing the progression of AMD.

Snail/Gfi-1 (SNAG) family zinc finger proteins in transcription regulation, chromatin dynamics, cell signaling, development, and disease.

The Snail/Gfi-1 (SNAG) family of zinc finger proteins is a group of transcriptional repressors that have been intensively studied in mammals. SNAG family members are similarly structured with an N-terminal SNAG repression domain and a C-terminal zinc finger DNA binding domain, however, the spectrum of target genes they regulate and the ranges of biological functions they govern vary widely between them. They play active roles in transcriptional regulation, formation of repressive chromatin structure, cellular signaling and developmental processes. They can also result in disease states due to deregulation. We have performed a thorough investigation of the relevant literature and present a comprehensive mini-review. Based on the available information, we also propose a mechanism by which SNAG family members may function.

Global analysis of target genes of 21 members of the ZAD transcription factor family in Drosophila melanogaster.

The zinc-finger associated domain (ZAD) family is the largest transcription factor family in dipteran insects. Still, their functional significance is barely recognized in the literature due in part to their resistance to mutagenesis screens in genetic studies. Therefore, we employed in vitro techniques to identify the DNA-binding characteristics of several members of the Drosophila melanogaster ZAD family in an effort to study their target genes. In this comprehensive investigation, we constructed a panel of GST-Zinc finger (ZnF) array chimera from 21 selected ZAD proteins and used them to select binding sites from an oligonucleotide library by employing electrophoretic mobility shift assays (EMSA). Samples of the binding population were sequenced and used to derive DNA-binding consensus sequence for each member. These consensus sequences were tested for complex formation with their respective protein chimera and the specificity of binding ascertained by competition EMSA. Bioinformatics tools were used to identify potential genetic targets. The identified consensus sequences were distinct for each member and the putative genomic targets were clustered in the regulatory regions of specific genes. This appears to be consistent with a conservation of function between members and also suggests that the overlapping functions of ZAD proteins are the result of positive selection to maintain redundancy and not simply artifacts of recent expansion. Putative target genes suggest a major role of the ZAD family members in the regulation of several early developmental genes including homeobox transcription factors.

An efficient and cost-effective protocol for selecting transcription factor binding sites that reduces isotope usage.

To function, transcription factors must position themselves by binding to DNA in a sequence-specific manner. Knowing the binding sites of these factors is a necessary step in understanding their activity. The standard protocols used for selecting a consensus-binding sequence for a DNA binding domain often require the use of radioisotopes to attain the necessary level of power in the assay. Alternatives are often less sensitive and may require an expensive apparatus for visualizing. We have created a modified binding site selection (BSS) protocol to improve efficiency and decrease the use of radioisotope. A GST affinity-tagged DNA binding domain construct was immobilized on a GSH affinity column and used to select from a randomized oligonucleotide library identical to those typically used in a radiolabeled BSS protocol. This produced a library specifically pre-enriched for use in a standard sequential EMSA selection. Use of a pre-enriched library reduced the total number of labeled rounds required for selection, decreasing the use of radioisotope while maintaining efficacy. The protocol was used to select for the binding sequence for several Drosophila melanogaster transcription factors. The consensus sequence was then shown by competitive binding experiments to associate with the protein in a sequence-dependent manner.

Combination of sulindac and dichloroacetate kills cancer cells via oxidative damage.

Sulindac is an FDA-approved non-steroidal anti-inflammatory drug with documented anticancer activities. Our recent studies showed that sulindac selectively enhanced the killing of cancer cells exposed to oxidizing agents via production of reactive oxygen species (ROS) resulting in mitochondrial dysfunction. This effect of sulindac and oxidative stress on cancer cells could be related to the defect in respiration in cancer cells, first described by Warburg 50 years ago, known as the Warburg effect. We postulated that sulindac might enhance the selective killing of cancer cells when combined with any compound that alters mitochondrial respiration. To test this hypothesis we have used dichloroacetate (DCA), which is known to shift pyruvate metabolism away from lactic acid formation to respiration. One might expect that DCA, since it stimulates aerobic metabolism, could stress mitochondrial respiration in cancer cells, which would result in enhanced killing in the presence of sulindac. In this study, we have shown that the combination of sulindac and DCA enhances the selective killing of A549 and SCC25 cancer cells under the conditions used. As predicted, the mechanism of killing involves ROS production, mitochondrial dysfunction, JNK signaling and death by apoptosis. Our results suggest that the sulindac-DCA drug combination may provide an effective cancer therapy.

The LIM protein AJUBA recruits protein arginine methyltransferase 5 to mediate SNAIL-dependent transcriptional repression.

The SNAIL transcription factor contains C-terminal tandem zinc finger motifs and an N-terminal SNAG repression domain. The members of the SNAIL family have recently emerged as major contributors to the processes of development and metastasis via the regulation of epithelial-mesenchymal transition events during embryonic development and tumor progression. However, the mechanisms by which SNAIL represses gene expression are largely undefined. Previously we demonstrated that the AJUBA family of LIM proteins function as corepressors for SNAIL and, as such, may serve as a platform for the assembly of chromatin-modifying factors. Here, we describe the identification of the protein arginine methyltransferase 5 (PRMT5) as an effector recruited to SNAIL through an interaction with AJUBA that functions to repress the SNAIL target gene, E-cadherin. PRMT5 binds to the non-LIM region of AJUBA and is translocated into the nucleus in a SNAIL- and AJUBA-dependent manner. The depletion of PRMT5 in p19 cells stimulates E-cadherin expression, and the SNAIL, AJUBA, and PRMT5 ternary complex can be found at the proximal promoter region of the E-cadherin gene, concomitant with increased arginine methylation of histones at the locus. Together, these data suggest that PRMT5 is an effector of SNAIL-dependent gene repression.

The Ajuba LIM domain protein is a corepressor for SNAG domain mediated repression and participates in nucleocytoplasmic Shuttling.

The SNAG repression domain is comprised of a highly conserved 21-amino acid sequence, is named for its presence in the Snail/growth factor independence-1 class of zinc finger transcription factors, and is present in a variety of proto-oncogenic transcription factors and developmental regulators. The prototype SNAG domain containing oncogene, growth factor independence-1, is responsible for the development of T cell thymomas. The SNAIL proteins also encode the SNAG domain and play key roles in epithelial mesenchymal differentiation events during development and metastasis. Significantly, these oncogenic functions require a functional SNAG domain. The molecular mechanisms of SNAG domain-mediated transcriptional repression are largely unknown. Using a yeast two-hybrid strategy, we identified Ajuba, a multiple LIM domain protein that can function as a corepressor for the SNAG domain. Ajuba interacts with the SNAG domain in vitro and in vivo, colocalizes with it, and enhances SNAG-mediated transcriptional repression. Ajuba shuttles between the cytoplasm and the nucleus and may form a novel intracellular signaling system. Using an integrated reporter gene combined with chromatin immunoprecipitation, we observed rapid, SNAG-dependent assembly of a multiprotein complex that included Ajuba, SNAG, and histone modifications consistent with the repressed state. Thus, SNAG domain proteins may bind Ajuba, trapping it in the nucleus where it functions as an adapter or molecular scaffold for the assembly of macromolecular repression complexes at target promoters.

Regulated recruitment of HP1 to a euchromatic gene induces mitotically heritable, epigenetic gene silencing: a mammalian cell culture model of gene variegation.

Heterochromatin protein 1 (HP1) is a key component of constitutive heterochromatin in Drosophila and is required for stable epigenetic gene silencing classically observed as position effect variegation. Less is known of the family of mammalian HP1 proteins, which may be euchromatic, targeted to expressed loci by repressor-corepressor complexes, and retained there by Lys 9-methylated histone H3 (H3-MeK9). To characterize the physical properties of euchromatic loci bound by HP1, we developed a strategy for regulated recruitment of HP1 to an expressed transgene in mammalian cells by using a synthetic, hormone-regulated KRAB repression domain. We show that its obligate corepressor, KAP1, can coordinate all the machinery required for stable gene silencing. In the presence of hormone, the transgene is rapidly silenced, spatially recruited to HP1-rich nuclear regions, assumes a compact chromatin structure, and is physically associated with KAP1, HP1, and the H3 Lys 9-specific methyltransferase, SETDB1, over a highly localized region centered around the promoter. Remarkably, silencing established by a short pulse of hormone is stably maintained for >50 population doublings in the absence of hormone in clonal-cell populations, and the silent transgenes in these clones show promoter hypermethylation. Thus, like variegation in Drosophila, recruitment of mammalian HP1 to a euchromatic promoter can establish a silenced state that is epigenetically heritable.

SETDB1: a novel KAP-1-associated histone H3, lysine 9-specific methyltransferase that contributes to HP1-mediated silencing of euchromatic genes by KRAB zinc-finger proteins.

Posttranslational modification of histones has emerged as a key regulatory signal in eukaryotic gene expression. Recent genetic and biochemical studies link H3-lysine 9 (H3-K9) methylation to HP1-mediated heterochromatin formation and gene silencing. However, the mechanisms that target and coordinate these activities to specific genes is poorly understood. Here we report that the KAP-1 corepressor for the KRAB-ZFP superfamily of transcriptional silencers binds to SETDB1, a novel SET domain protein with histone H3-K9-specific methyltransferase activity. Although acetylation and phosphorylation of the H3 N-terminal tail profoundly affect the efficiency of H3-K9 methylation by SETDB1, we found that methylation of H3-K4 does not affect SETDB1-mediated methylation of H3-K9. In vitro methylation of the N-terminal tail of histone H3 by SETDB1 is sufficient to enhance the binding of HP1 proteins, which requires both an intact chromodomain and chromoshadow domain. Indirect immunofluoresence staining of interphase nuclei localized SETDB1 predominantly in euchromatic regions that overlap with HP1 staining in nonpericentromeric regions of chromatin. Moreover, KAP-1, SETDB1, H3-MeK9, and HP1 are enriched at promoter sequences of a euchromatic gene silenced by the KRAB-KAP-1 repression system. Thus, KAP-1 is a molecular scaffold that is targeted by KRAB-ZFPs to specific loci and coordinates both histone methylation and the deposition of HP1 proteins to silence gene expression.