The SARS-Coronavirus-Host Interactome: Identification of Cyclophilins as Target for Pan-Coronavirus Inhibitors
Coronaviruses (CoVs) are important human and animal pathogens that induce fatal respiratory, gastrointestinal and neurological disease. The outbreak of the severe acute respiratory syndrome (SARS) in 2002/2003 has demonstrated human vulnerability to (Coronavirus) CoV epidemics. Neither vaccines nor therapeutics are available against human and animal CoVs. Knowledge of host cell proteins that take part in pivotal virus-host interactions could define broad-spectrum antiviral targets. In this study, we used a systems biology approach employing a genome-wide yeast-two hybrid interaction screen to identify immunopilins (PPIA, PPIB, PPIH, PPIG, FKBP1A, FKBP1B) as interaction partners of the CoV non-structural protein 1 (Nsp1). These molecules modulate the Calcineurin/NFAT pathway that plays an important role in immune cell activation. Overexpression of NSP1 and infection with live SARS-CoV strongly increased signalling through the Calcineurin/NFAT pathway and enhanced the induction of interleukin 2, compatible with late-stage immunopathogenicity and long-term cytokine dysregulation as observed in severe SARS cases. Conversely, inhibition of cyclophilins by cyclosporine A (CspA) blocked the replication of CoVs of all genera, including SARS-CoV, human CoV-229E and -NL-63, feline CoV, as well as avian infectious bronchitis virus. Non-immunosuppressive derivatives of CspA might serve as broad-range CoV inhibitors applicable against emerging CoVs as well as ubiquitous pathogens of humans and livestock.
Published in the journal:
. PLoS Pathog 7(10): e32767. doi:10.1371/journal.ppat.1002331
Category:
Research Article
doi:
https://doi.org/10.1371/journal.ppat.1002331
Summary
Coronaviruses (CoVs) are important human and animal pathogens that induce fatal respiratory, gastrointestinal and neurological disease. The outbreak of the severe acute respiratory syndrome (SARS) in 2002/2003 has demonstrated human vulnerability to (Coronavirus) CoV epidemics. Neither vaccines nor therapeutics are available against human and animal CoVs. Knowledge of host cell proteins that take part in pivotal virus-host interactions could define broad-spectrum antiviral targets. In this study, we used a systems biology approach employing a genome-wide yeast-two hybrid interaction screen to identify immunopilins (PPIA, PPIB, PPIH, PPIG, FKBP1A, FKBP1B) as interaction partners of the CoV non-structural protein 1 (Nsp1). These molecules modulate the Calcineurin/NFAT pathway that plays an important role in immune cell activation. Overexpression of NSP1 and infection with live SARS-CoV strongly increased signalling through the Calcineurin/NFAT pathway and enhanced the induction of interleukin 2, compatible with late-stage immunopathogenicity and long-term cytokine dysregulation as observed in severe SARS cases. Conversely, inhibition of cyclophilins by cyclosporine A (CspA) blocked the replication of CoVs of all genera, including SARS-CoV, human CoV-229E and -NL-63, feline CoV, as well as avian infectious bronchitis virus. Non-immunosuppressive derivatives of CspA might serve as broad-range CoV inhibitors applicable against emerging CoVs as well as ubiquitous pathogens of humans and livestock.
Introduction
Five distinct CoVs (SARS-CoV, hCoV-NL63, hCoV-HKU-1, hCoV-OC43, hCoV-229E) cause respiratory tract illness in humans, ranging from mild common cold to deadly virus-associated pneumonia [1]. At least seven different animal CoVs cause economically significant epizootics in livestock, and deadly disease in companion animals [1]. The agent of SARS was a novel CoV introduced into the human population from an animal reservoir, resulting in a highly lethal epidemic in 2002/2003 [1], [2]. A tremendous diversity of CoVs exists in complex mammalian and avian reservoirs [1], [3], [4]. Host switching is a common feature in CoV evolution, and novel epidemic CoV can emerge anytime [1], [3], [5]. Because the large diversity of CoVs complicates the design of vaccines, the identification of broad-range anti-CoV drug targets might indicate alternative approaches against CoV epidemics [1]. Broad range anti-CoV drugs would also be desirable to treat severe infections caused by known human and animal CoVs.
The SARS-CoV genome is predicted to encode 14 functional open reading frames, leading to the expression of up to 29 structural and non-structural protein products [1]. The functions of many of these proteins are poorly understood or unknown. To study the interplay of viral proteins with the host cell and to identify new targets involved in viral replication we have performed a genome-wide analysis of protein - protein interactions between the SARS-CoV and human proteins via a High-Throughput Yeast Two Hybrid Screen (HTY2H) [6], [7]. Within this framework we identified redundant interactions between SARS-CoV non-structural protein Nsp1 and a group of host proteins with peptidyl-prolyl cis-trans-isomerase activity, including the cyclophilins/immunophilins PPIA, PPIG, PPIH and FKBP1A, FKBP1B. These modulate the Calcineurin/NFAT pathway that plays an important role in immune cell activation [8], [9]. The NFAT family of transcription factors encodes four calcium-regulated proteins of which three (NFAT1, -2, -3) are expressed in a variety of cell types including T-cells, B-cells, mast cells, natural killer cells and eosinophils [8], [9]. NFAT activation regulates pivotal immune processes like apoptosis, anergy, and T-cell development. An essential activation step for NFAT is its dephosphorylation by the phospatase calcineurin A (CnA), resulting in the translocation of NFAT into the nucleus. Cyclosporin A (CspA) forms complexes with cyclophilins that bind to CnA, preventing its activity. This effect is used in transplant patients to prevent organ rejection by suppression of the immune system. Here we show that SARS-CoV nonstructural protein Nsp1, as well as full replicating SARS-CoV, enhance the CnA/NFAT pathway and induce NFAT-responsive promoters. Because interactions with upstream elements of the pathway were redundantly identified in a hypothesis-free virus-host interaction screen, the pathway is likely to play a significant role for virus replication. Indeed, an extensive panel of CoVs covering all three relevant virus genera was strongly inhibited by manipulation of cyclophilins using CspA.
Results
Interaction screening of the SARS-CoV ORFeome and host proteins
All SARS-CoV ORFs and a number of subfragments lacking transmembrane regions were cloned into eukaryotic expression vectors. Using HTY2H, these were screened against a cDNA library of very high complexity (1.4×107) derived from human brain, as well as an additional library of individually-cloned full-length ORFs encoding 5000 human proteins. Inserts from positive yeast clones were sequenced and compared against GenBank. BLAST searches on 2287 DNA sequences yielded 942 different human gene hits. These were divided into four confidence categories: category A (highly confident interaction partners found more than once in one or several screens), category B (single hits), category C (sticky preys interacting with several to many bait proteins) and category D (3′-UTR cDNA regions or inserts in reverse orientation coding for unnatural peptides). We found 132, 383, 245, and 282 hits in categories A – D, respectively. For validation, the cDNAs of 86 category A and category B interaction candidates were cloned in-frame with the Renilla reniformis luciferase and overexpressed in HEK 293 cells. SARS-CoV ORFs were cloned in-frame with N-terminal protein A domains and co-expressed in the same cells. Protein A-directed immunoprecipitates retained on IgG-coated magnetic beads were identified by measuring in-vitro Luciferase activity. About 48% of category A candidates and 36% of category B candidates were confirmed positive with a Z-score >1 (Figure 1, see Materials and Methods for definition), corresponding to previous observations [10]. A list of validated category A and B HTY2H interactor candidates is provided in Table S1.
For an overall estimate of plausibility, more than 5,000 Medline abstracts mentioning “SARS” or “Coronavirus” were screened using the text mining program syngrep, scanning for the mentioning of human protein designations and synonyms. Abstracts mentioning YTH or co-immunoprecipitation assays were specifically sought. Twenty-eight CoV-/host protein interactions were identified in the literature, as listed in Table S2. It was then determined how these literature hits overlapped with the lists of candidate interactors as identified by HTY2H screening in different confidence levels. Using a hypergeometric test, the fractions of overlap were compared to the fraction of literature hits in the list of search terms (31,941 human proteins used for text-mining). Abstracts were enriched for proteins identified as SARS-CoV interaction partners both in the high confidence and the complete data sets (Table S3 and Table S4). Figure 2 summarizes highly confident interactions identified in the overall screen and GO [11], [12] analysis. SARS-CoV proteins were found to preferentially target protein complex subunits (Table S5 and Table S6). Of 9 complexes which were targeted through ≥4 subunits, 4 complexes were found to be significantly enriched: The respiratory chain complex I (7 subunits targeted by SARS-CoV, p-value ∼0.036), the cytoplasmic ribosome (10 subunits targeted by SARS-CoV, p-value ∼0.036), in particular the 60S ribosomal subunit (7 subunits targeted by SARS-CoV, p-value ∼0.036) and the LCR-associated remodeling complex which is involved in DNA conformation modification (4 subunits targeted by SARS-CoV, p-value ∼ 0.039). Furthermore, the analysis of the centrality of SARS targets within the human interaction network (Figure S1) indicated that SARS-CoV proteins target both highly interactive proteins (hubs) as well as so-called bottleneck proteins which are central to many of the shortest paths in their networks [13] (Figure S2). Table S4 summarizes GO results for SARS-CoV nonstructural protein Nsp1, a protein yielding particularly interesting candidate interaction networks. Interactions between Nsp1 and several members of the class of immunophilins (PPIA, PPIG, PPIH, FK506-binding proteins FKBP1A and -B) and calcipressins (RCAN1 and -3) were selected for experimental confirmation.
The N-terminal part of SARS-CoV Nsp1 influences NFAT activation by interacting with several immunophilins and a calcipressin
The immunophilin proteins (cyclophilins and FK506-binding proteins) are all known to bind to CnA in combination with inhibitory molecules, and to influence the CnA/NFAT pathway that plays a major role in the establishment of T-cell immune response [14]. For a more detailed mapping of HTY2H hits, PPIA, PPIB, PPIG, PPIH, FKBP1A and RCAN3, and three versions of Nsp1 [Nsp1(aa 1–180), Nsp1(aa 1–93) and Nsp1(aa 119–180)] were cloned into LUMIER assay vectors to yield luciferase and protein A fusion proteins, respectively. Although PPIB was not identified as an interactor by Y2H it was included in the experiment as it is known to bind to the HIV-1 gag protein [15] and to the HCV NSB5 protein [16]. All tested proteins interacted with Nsp1(aa 1–93), suggesting redundant interactions of SARS-CoV with the CnA/NFAT pathway via the N-terminal part of Nsp1 (Figure 3). To examine the functional consequences of Nsp1 expression on NFAT activity, NFAT and CnA cDNAs were overexpressed in HEK 293 cells. Parallel experiments in Jurkat cells were done without overexpression due to their constitutive activity of the CnA/NFAT pathway. The CnA/NFAT pathway was stimulated by addition of PMA (40 ng/ml) and ionomycin (2 µM) to the culture medium. In both cell lines treated this way, expression of Nsp1 did not induce NFAT activity directly, but increased significantly the stimulatory effect of PMA/ionomycin on NFAT activation (Figure 4A). The increase in NFAT activity could be blocked by CspA, an inhibitor of the NFAT pathway (Figure 4A and B). Coexpression of the calcipressin RCAN3 as shown in Figure 4B attenuated the overall stimulating effect on the NFAT activity. In contrast, coexpression of other CoV proteins or coexpression of PPIA, PPIH, FKBP1A did not impact NFAT activity (data not shown). Experiments up to this point employed overexpression of NFAT3. As different NFAT species are expressed depending on cell type [17], NFAT1 and NFAT2 were alternatively expressed and compared in the same assay. For both species essentially the same influence of Nsp1 on PMA/ionomycin-dependent stimulation was seen (Figure 4C and D). Altogether this suggested a broad effect of Nsp1 on NFAT activation that is mediated via the canonical NFAT activation pathway including CnA.
In order to determine the extent of PMA/ionomycin-dependent NFAT activation during virus infection, HEK 293 lp cells (lp = low passage) with a short passage history were infected with SARS-CoV at an MOI = 1. These cells had been previously demonstrated to support SARS-CoV replication, in contrast to common HEK 293 cells [18], [19]. Figure 5 shows that the CnA/NFAT pathway was induced in the context of SARS-CoV infection at considerable extent, and in a PMA/ionomycin-dependent way.
Effects on cytokine induction
Viruses may interfere with cytokine induction, but on the other hand, may also induce cytokine genes directly. To examine whether the Nsp1-mediated, PMA/ionomycin-dependent activation of NFAT may cause specific induction of relevant cytokines, HEK293 cells were co-transfected with the plasmids described above, except that the NFAT reporter plasmid was replaced by luciferase reporter plasmids carrying the IL-2, IL-4 and IL-8 promoters, respectively (Figure 6). Expression of Nsp1 induced the IL-2 promoter significantly by a factor of about 2.5 (Figure 6A). This effect was inhibited by CspA and RCAN-3, suggesting dependence on the CnA/NFAT pathway. The IL-4 promoter activity was not significantly elevated by Nsp1 expression in the presence of PMA/ionomycin (Figure 6B). Its activity was decreased in the presence of CspA but not RCAN-3. The IL-8 promoter was induced by PMA/ionomycin alone, but significantly downregulated by a factor of about 1.8 in additional presence of Nsp1 (Figure 6C). Expression of RCAN-3 reduced IL-8 promoter activity levels to about half, while CspA inhibited the promoter completely. In Jurkat cells, which express endogenous NFAT3 and CnA, the Nsp1 protein did not induce the IL-2 promoter. The slight induction of IL-4 and the downregulation of the IL-8 promoter activities in presence of Nsp1 (about twofold) were similar to effects seen in HEK 293 cells. These results suggested that Nsp1 expression had the strongest influence on the IL-2 promoter.
Next to NFAT, transcription factors NFκB and Activating Protein 1 (AP-1) determine IL-2 regulation [20]. NFAT, AP-1 and NFκB binding sites are juxtaposed in the IL-2 promoter, and it has been shown that NFAT and AP-1 act in a cooperative manner on the promoter while NFkB has enhancing function [17]. Simultaneously, NFκB induces the IFN-beta gene by binding to the PRDII DNA element [21]. The latter is a more sensitive assay of NFκB nuclear translocation upon viral infection. To examine potential direct effects of Nsp1 on NFκB nuclear translocation and AP-1, HEK 293 and Jurkat cells were cotransfected with SARS-CoV Nsp1fl and p55A2luc containing repeated PRDII elements or pAP-1-luc containing the AP-1 binding site of the IL-2 promoter (Figure 7). Overexpression of Nsp1fl as well as treatment with PMA/ionomycin, respectively, caused small but significant luciferase increases in both cell lines. The combined expression of Nsp1fl with PMA/ionomycin treatment led to significant induction of PRDII by a factor of about 6 in both cell lines. The AP-1 promotor was only slightly upregulated in HEK293 and downregulated in Jurkat cells. This indicated a co-involvement of NFκB but not of AP-1 in the induction of IL-2 by Nsp1, suggesting dependence mainly on the NFAT pathway. In summary, SARS-CoV caused relevant and specific induction of IL-2 by activating the NFAT pathway via Nsp1.
CspA inhibits replication of CoVs
CspA is a highly efficient antagonist of NFAT activation, interacting with cyclophilins. Due to the high specificity of Nsp1-dependent activation of NFAT and due to the high redundancy of SARS-CoV interactions with upstream elements of the CnA/NFAT pathway, we suspected an essential function for the virus. It was therefore investigated whether CspA might influence viral replication (Figure 8). Vero cells were inoculated with a low dose of SARS-CoV (MOI = 0.0001) and growth of virus replication was determined by real-time RT-PCR and plaque titration. In parallel cell cultures treated with the same concentrations of CspA, cell viability was measured with a highly sensitive assay based on ATP provision in metabolically active cells. A profound and dose-dependent inhibition of replication of SARS-CoV strain Frankfurt-1 in Vero E6 cells was seen in absence of cytopathic effects conferred by the compound. The 50% effective inhibitory concentration was 3.3 µM.
Because Nsp1 proteins of group I and SARS coronaviruses share structural and functional similarities [22], it was tested whether the inhibitory effect of CspA could be extended to other pathogenic CoVs. These included members of the genera Alphacoronavirus (human CoV-NL63 and -229E, Feline CoV serotypes I and II [strains Black and 791146], porcine transmissible gastroenteritis virus [TGEV]), Betacoronavirus (SARS-CoV isolates Frankfurt and Hongkong) and Gammacoronavirus (avian infectious bronchitis virus [IBV]). All tested CoVs were inhibited by CspA; replication of TGEV and IBV in the tested range (up to 25 µM) was diminished close to background by CspA. HCoV-NL63 and -229E and the two Feline CoV serotypes were completely inhibited, with 50% effective concentrations of 2.3 µM, 2.3 µM and 2.7 µM, respectively (Figure 8). Figure S3 shows reduction of virus replication in a log scale.
CspA inhibits a SARS-CoV replicon
In order to determine the principal stage of the CoV replication cycle inhibited by CspA, a novel SARS-CoV replicon carrying a secreted Metridia luciferase reporter construct instead of the major structural proteins S, E, and M was used (Figure 9A). The replicon RNA together with an mRNA for the nucleocapsid protein was electroporated in BHK cells. Replicon activity in parallel reactions was controlled to be at the same level after 16 h of incubation (data not shown), and increasing amounts of CspA were added to cells after repeated washing. As shown in Figure 9B, accumulated luciferase activities in supernatants were decreased in a CspA dose-dependent manner after 24 h. Two different specific inhibitors of the CoV main protease, Cinanserin [23] and XP17 (R. H., own unpublished observations), inhibited replicon activity to a comparable extent as CspA, at comparable substance concentrations (Figure 9B). To control against any influence of the nucleocapsid protein that is co-electroporated for maximal replicon efficiency [24], [25], [26] and that is also contained in the replicon RNA, this protein was expressed from a eukaryotic expression vector in the same cells and an NFAT induction assay was conducted as described above. No N-dependent effect on the assay was seen (Figure 9C). These results suggest an action of CspA on genome replication and/or transcription, rather than other stages such as virus entry or egress.
Discussion
Various genomic and proteomic methods have been utilized to identify protein-protein interactions in the context of viral replication [6], [27], [28], [29], [30], [31]. HTY2H is among the most direct approaches to identify interactions between members of viral ORFeomes and large host cDNA libraries. Major advantages of the method include its potential for high throughput testing and automation, as well as its high sensitivity. The latter facilitates investigation of proteins expressed at low levels and of those causing weak and transient interactions [32], [33]. Drawbacks include the inability to control and confirm expression of genes of interest, other than by positive selection of yeast expression clones containing nutritional markers. Moreover, some proteins need posttranslational modifications not provided by the yeast cell, in order to interact with binding partners. Also, since interactions of bait and prey proteins take place in the nucleus of the yeast cell, the assay is influenced by hydrophobic and transmembrane regions affecting the nuclear membrane. It is well known that a considerable fraction of interacting proteins in HTY2H represent false positive findings, making it absolutely necessary to validate interactions by independent eukaryotic assays. We have implemented a version of the Lumier assay that is amenable for screening in mammalian cells at a medium scale of parallelity [34]. Our modified version using a protein A tag instead of a Flag tag enables automated capture of precipitates on IgG Fc-coated magnetic beads. Throughput is mainly limited by the requirement to subclone Y2H plasmid inserts, as the assay does not involve any cell-based imaging or other readouts going beyond in-vitro Luciferase assays. In our analysis of 86 Y2H-positive interaction partners we achieved a positive confirmation rate of about 42% in category A and B interactors, which is in good agreement with a recent standardized comparison of five different interaction assays in which the LUMIER pull-down assay showed the highest sensitivity (36%) on a positive reference set of human proteins [10]. It has to be mentioned that all interaction assay systems carry intrinsic limitations. In the case of our modified Lumier method the Renilla and protein A tags are rather long as compared to His or HA tags. Therefore, true interactions might be missed, and it is possible that the interaction of PPIA with full-length nsp1 is sterically prevented by the length of these tags, as compared to the N-terminal fragment of Nsp1. But the chief attraction of our method is its applicability in high thoughput assays.
The range of interactors identified in this study defines an unprecedented resource for future investigations into pathogenetic mechanisms and antiviral applications against CoVs. In order to demonstrate that HTY2H can afford a direct identification of novel antiviral targets, we have chosen one promising group of interactors for further investigation in the present study. The interaction of Nsp1 with the cyclophilins PPIA, PPIB, PPIG, PPIH, the FK506-binding proteins FKBP1A/B, and the CnA (calcipressin) regulators RCAN1 and RCAN3 represented a highly redundant virus-host interaction involving critical elements of the same regulatory network immediately upstream of the CnA/NFAT pathway. Nsp1 is a virulence factor in-vivo whose action has been linked with early stages of the immune response, including antagonism against IFN signaling and inhibition of host protein synthesis [35], [36], [37], [38]. Our findings add an important new dimension to Nsp1′s role in pathogenicity, identifying this protein as a strong and specific activator of NFAT enhancing the induction of the IL-2 promoter. The increase of NFAT activation extended to all three major NFAT species, suggesting a potential for induction of broad and systemic cytokine dysregulation by affecting several types of immune cells. The pattern of cytokine dysregulation in severe SARS cases differed from the cytokine burst seen in other acute viral diseases in its delayed occurrence, manifesting beyond the second week of symptoms. Interestingly, it was noted upon clinical observations that late aggravation was correlated with severe clinical outcome, triggering parallel efforts in several centers to treat patients empirically with steroids [39], [40], [41]. Our results suggest an influence on fundamental triggers of immune cell activation, contributing an explanation for the cytokine dysregulation and immune-dependent pathogenesis observed in severe cases of SARS.
The detected interactions with immunophilins caught our attention because peptide inhibitors such as CspA and Tacrolimus (FK506) are available that bind to immunophilins and cause dramatic effects on the cellular phosphatase CnA, suppressing the CnA/NFAT immune-regulatory pathway [42]. Cyclophilins are essential cofactors for replication of HCV, HIV and of some parasites [43]. Incorporation of PPIA into HIV-1 particles via binding to gag is paralleled by its binding to the N protein of SARS-CoV by an educated guess finding using surface plasmon resonance biosensor technology [44]. Furthermore, in a recent proteomics study in which viral and cellular proteins incorporated into SARS-CoV virions were spectrometrically profiled PPIA was also found in purified virus particles [31]. We have indeed shown that the inhibition of cyclophilins by CspA exhibited strong and specific inhibitory effects on members of all three genera of CoVs. This broad antiviral action is supported by a recent characterization of Nsp1 structural conservation that extends beyond the limits of CoV genera [22]. Inhibition took place in the low micromolar range, indicating prospects for future investigation of similar (non-immunosuppressive) drugs as broad-range antivirals. Such drugs have already been employed against HCV, whose nonstructural NS5A and NS2 proteins interact with cyclophilins and whose replication can be inhibited by CspA and non-immunosuppressive derivatives thereof [45].
Our study is limited in that it does not clarify the biological functions of the interaction between Nsp1 and its cellular partners, nor does it confirm the involvement of Nsp1 in full virus context by knock out experiments in recombinant viruses. Brockway and Denison showed that deletion of residues in the amino-terminal half of nsp1 is not tolerated for a productive infection [46]. Thus it seems not possible to construct viruses with mutations in the N- terminal half of nsp1 which do not bind to cyclophilins anymore. It will thereofore be difficult to exactly delineate the role of nsp1 in virus-host cooperation. It should also be appreciated that CspA-inhibitable NFAT induction by SARS-CoV nsp1 may be independent of the replication inhibition by CspA. CspA specifically binds to cyclophilins and this complex binds to calcineurin phosphatase preventing the dephosphorylation of NFAT. Binding of nsp1 to cyclophilins and the induction of NFAT is obviously inhibited by Csp A. On an independent level, formation of the Cyclophilin-CspA complex might prevent those cyclophilin functions required for virus replication. Complex further studies involving multiple CoV systems will be required to delineate these functions.
Nevertheless, we have been able to shed more light on the principal stage of the virus replication cycle that is subject to CspA-dependent inhibition. Our experiments using a SARS-CoV replicon suggest strongly that the processes afforded by the replicative proteins rather than stages of virus entry and egress, are affected. Further experiments need to be done in the future in order to investigate on a mechanistic level the potential breadth of the identified antiviral effect. In particular, the large diversity of CoVs in animal reservoirs generates interest in studying how conserved this particular virus-host interaction might be between non-human CoVs and human cells, and whether this could be exploited as a truly broad-range antiviral target that covers epidemic and reservoir-borne viruses alike.
Materials and Methods
Cells
HEK293, HEK293lp (low passage), Vero E6, CaCo-2, HRT 18, Huh 7, FCWF- and St-cells were grown in Dulbecco’s modified Eagle medium (DMEM) containing 10% FBS, 1% L-glutamine, 1% penicillin/streptomycin and 1% non essential amino acids. C8166 SEAP cells were cultured in RPMI medium with 10% FBS and 1% penicillin/streptomycin. Jurkat cells were propagated in RPMI-medium containing 10% FBS, 1% L-glutamine and 1% penicillin/streptomycin.
Y2H library screening
Cloning of the SARS ORFeome into destination bait vector pGBKT7-DEST for Y2H screening was described previously [6]. A series of SARS-CoV (Frankfurt isolate) bait clones containing ORF fragments depleted of transmembrane regions were generated in addition (primers and cloning procedures available on request). Automated yeast two-hybrid screens were essentially done as described previously [7], [47], with the following modifications. Human cDNA libraries from human brain and fetal brain (Clontech) as well as a library of individually cloned full–length open reading frames from cDNAs of 5000 different genes were screened to a minimal coverage of 5 million clones per library. To mate yeast strains, exponentially growing cultures at an OD600nm of 1 were combined, pelleted by centrifugation, and resuspended in an equal volume of YPDA (Yeast Extract Peptone Dextrose Adenine) containing 20% PEG 6000. For the generation of a high-confidence dataset, interaction pairs were selected which were isolated at least twice, or where the bait interacted with two highly related preys, and which did not involve promiscuous preys.
Modified LUMIER assays
For LUMIER assays, proteins were transiently expressed in HEK293 cells as N-terminal fusion proteins with the Staphylococcus aureus protein A tag or Renilla reniformis luciferase. 20 ng of each expression construct were transfected into 10,000 HEK293 cells using 0.05 µl of lipofectamine 2,000 (Invitrogen) in 96 well plates. After 40 hours, medium was removed and cells were lysed on ice in 10 µl of ice-cold lysis buffer (20 mM Tris pH 7.5, 250 mM NaCl, 1% TritonX-100, 10 mM EDTA, 10 mM DTT, Protease Inhibitor Cocktail [Roche # 1 836 170], Phosphatase Inhibitor Cocktail [Roche, # 4 906 837], Benzonase [Novagen #70746, 0,0125 units per µl final concentration]) containing sheep-anti-rabbit IgG-coated magnetic beads (Invitrogen, Dynabeads M280, 2 mg/ml final concentration). Lysates were incubated on ice for 15 minutes. 100 µl of washing buffer (PBS, 1 mM DTT) were added per well, and 10% of the diluted lysate was removed to determine the luciferase activity present in each sample before washing. The rest of the sample was washed 6 times in washing buffer in a Tecan Hydroflex plate washer. Luciferase activity was measured in the lysate as well as in washed beads. Negative controls were transfected with the plasmid expressing the luciferase fusion protein and a vector expressing a dimer of protein A.
For each sample, four values were measured: the luciferase present in 10% of the sample before washing, the luciferase activity present on the beads after washing, and the same values for the negative controls. Normalised interaction signals were calculated as follows: Log(bound)/log(input) – log(bound nc)/log(input nc). Z-scores were calculated by subtracting the mean and dividing by the standard deviation. The mean and standard deviation were calculated from large datasets of protein pairs which were not expected to interact, i.e. from negative reference sets. Normalised signal to noise ratios were calculated as follows: (bound/input)/(bound nc/input nc)[10].
Reporter gene assays
For reporter gene expression HEK293 cells were transiently transfected with 2 µg DNA containing 460 ng of the respective expression plasmids encoding NFAT, calcineurin, SARS-CoV Nsp1 and reporter genes as indicated (six well plates). DNA was transfected using FuGENE HD reagent (Roche Applied Science). For reporter gene expression in Jurkat cells (1×106 cells) 1 µg reporter plasmid and 1 µg expression plasmid encoding SARS-CoV Nsp1 were transiently cotransfected using the Amaxa Cell Line Nucleofector Kit (Lonza). For reporter gene assays in the viral context HEK293lp cells (24-well plates) were transiently transfected with 500 ng reporter plasmid using Lipofectamine LTX (Invitrogen) according to the manufacturer’s instructions. 24 h post transfection cells were infected with SARS-CoV HK at an MOI = 1. 19 h after transfection cells were harvested and Promega’s dual luciferase assays were performed according to the manufacturer’s instructions. All results were normalized to a simultaneously transfected Renilla luciferase (pRL-null, Promega). Cells were stimulated by PMA (40 ng/ml) and Ionomycin (2 µM) in culture medium. The NFAT-pathway was inhibited by Cyclosporin A (50 ng/ml). Virus infected cells were lysed 17 h post infection with Promega passive lysis buffer (PLB) and 1% Igepal (Sigma-Aldrich). All experiments were repeated at least three times.
Virus inhibition experiments and plaque titration
SARS - CoV
Vero E6 cells were seeded in 6-well plates. At a confluence of 80%, cells were washed with phosphate buffered saline (PBS) and infected with a MOI = 0.0001. After 1 h adsorption cells were washed twice with PBS and supplemented with 2 ml fresh media and inhibitor in different concentrations. Supernatants were tested after 24 h p.i. and 48 h p.i.. For quantitative real time RT-PCR 140 µl supernatant were taken and analysed as described previously [48]. Plaque titrations on Vero E6 cells were performed using Avicel overlays (RC581, FMC BioPolymer, Belgium) as described [18].
HCoV-NL63
CaCo-2 cells were seeded in 24-well plates. For quantitative real time RT-PCR cells were infected with a MOI = 0.004 at a confluence of 100%. After 1 h adsorption the viral inoculum was removed, cells were washed twice with PBS and 1 ml fresh medium supplemented with different inhibitor concentrations was added to the cells. 70 µl supernatant were tested on day 2 and day 4 p.i. Viral RNA was extracted from cell culture supernatant with the QIAamp Viral RNA mini Kit (QIAGEN, Hilden, Germany). Real time RT-PCR and determination of viral replication by plaque assay was done as described previously [49]. Various inhibitor concentrations were directly added to the overlay medium. Cells were fixed and plaques were stained with 0.2% crystal violet, 11% formaldehyde, 2% ethanol, and 2% paraformaldehyde.
HCoV-229E
Huh-7 cells were seeded in 24-well plates and incubated until the monolayer was 70–80% confluent. Cells were infected with HCoV229E-luc (V. Thiel unpublished data), in which ORF 4 is replaced by Renilla luciferase, at an MOI = 0.1. After 1 h adsorption viral inoculum was removed, cells were washed with PBS and incubated with 1 ml fresh medium containing different inhibitor concentrations. Viral replication was determined 24 h and 48 h p.i. by Renilla Luciferase Assay System (Promega, Madison, USA) according to the manufacturer’s instructions.
IBV Beaudette
Vero cells grown on 24-well plates were infected with IBV Beaudette at an MOI = 0.00025. Methylcellulose was added after 1 h incubation at 37°C. Different concentrations of CspA at a range from 0 to 25 µM were mixed with the methylcellulose. After 24 h, the cells were fixed with 3% paraformaldehyde. Plaques were stained with an anti-IBV polyclonal serum raised in rabbits and a FITC-labelled secondary antibody.
TGEV PUR46
St cells were seeded in 6-well plates. At a confluence of 100%, cells were washed twice with PBS and infected with a MOI = 0.000066. Plaque assays were performed as described previously [50]. Different inhibitor concentrations were directly added to overlay media.
FCoV Black (Serotype I) and 791146 (Serotype II)
One day before the experiment 1.5 * 106 FCWF cells were seeded in 12-well plates. Cells were infected with respective virus dilutions. 1 h after infection cells were washed with 1 ml PBS and 2 ml fresh medium with different inhibitor concentrations were added. 48 h p.i. supernatants were analyzed via plaque assay on freshly seeded FCWF cells with 100% confluence. Again, cells were infected for 1 h, after washing with 1 ml PBS, 1 ml 1% Carboxymethylcellulose medium was added to cells. After 72 h the plaques were analyzed.
EMCV
Vero E6 cells were seeded in 6-well plates and infected with an MOI = 0.01. After 1 h infection medium was removed and replaced by fresh medium supplemented with different inhibitor concentrations. 24 h p.i. supernatants were analyzed by plaque titration on freshly seeded Vero E6 cells. Cells were infected for 1 h, washed with 1 ml PBS and incubated with 0.4% Noble Agar (Difco) overlay medium for 72 h. After fixation with 2 ml 5% tri-chloro acetic acid (TCA) for 5 min, plaques were stained with 1% crystal violet in 3.6% formaldehyde, 1% methanol, and 20% ethanol.
HIV-1
C8166 SEAP cells were seeded at a density of 30,000 cells per 96-well and CspA was added to a final concentration of 0.2, 0.4, 0.8, 2, 4, and 8 µM. Cultures were inoculated with HIV-1 NL4-3 at low MOI that allowed complete inhibition of replication in the presence of the HIV protease inhibitor Nelfinavir (NFV) to ensure multiple rounds of infection. Half of the culture medium was discarded on day 3 and replaced by fresh medium ± inhibitor. HIV-1 replication was monitored microscopically via cytopathic effects and measured on day 5 via the activity of secreted alkaline phosphatase (SEAP) in the cell culture supernatant employing the Phospha-Light SEAP Reporter Gene Assay (Applied Biosystems).
In all experiments the outcome of the mock-treated cells were set as 100%. Data shown are the mean values of at least three experiments. HIV-1 data show one representative experiment out of three, values are averages of triplicates. DMSO was always added as a control corresponding to the highest inhibitor concentration.
Cytotoxicity assay
Cytotoxicity tests of all cell lines were carried out in a 96-well format with the CellTiter-Glo Luminescent Cell Viability Assay (Promega, Madison,USA) according to the manufacturer’s instructions.
SARS replicon and inhibitor assay
SARS-CoV replicon
The SARS replicon was based on our previously pubished infectious SARS-CoV cDNA clone [51]. By PCR mutagenesis and classical cloning techniques ORFs 2–8 were deleted or replaced by marker genes, respectively. In ORF 2 (spike gene) a secreted Metridia Luciferase (Promega) was inserted, using the original transcription-regulating sequence. Furthermore, for potential selection, a GFP-Zeocin fusion protein was placed in open reading frame 3a. Cloning details will be given by authors upon request.
Inhibitor assays
All assays were performed in triplicates. Replicon RNA and N RNA were transcribed and electroporated in BHK cells as described [51]. For inhibitor assays, cells were seeded in 6 well plates after electroporation using 2×104 cells per well and incubated for 16 hours at 37°C. Hereafter, supernatant was removed entirely and cells were washed twice with PBS. Then 2 ml of fresh medium with indicated concentration of inhibitors were added. For each timepoint 50 µl of supernatant were removed for Luciferase assay (Promega). RLUs were measured as recommended by the manufacturer.
Gene Ontology (GO) over-representation analysis
To determine which cellular pathways are targeted by SARS-CoV, functional categories enriched among SARS cellular interaction partners were identified using a Gene Ontology (GO) over-representation analysis. For this purpose, p-values were determined with the hypergeometric test implemented in the Ontologizer software [12]. P-Values were corrected for multiple testing using the FDR-method by Benjamini and Hochberg [52] and significant terms were identified at a threshold of 0.05.
Analysis of protein complexes
We analyzed whether proteins involved in protein complexes were preferentially targeted by SARS-CoV proteins and which complexes were preferentially targeted. For this purpose, protein complexes for humans were extracted from the CORUM database[53]. After removing complexes which were identical to another complex, we obtained a data set of 1184 complexes containing 2079 distinct proteins. P-values for the enrichment of protein complex subunits among SARS targets were calculated with the hypergeometric test assuming a background of ∼ 25,000 proteins (Table S4). We furthermore analyzed 9 complexes which targeted at least four subunits by SARS proteins. Using a hypergeometric test, we determined p-values for the enrichment of SARS target proteins among the subunits of each complex. P-values were corrected for multiple testing using the FDR- method by Benjamini and Hochberg (Table S5).
Centrality of virus targets
Interactions between SARS proteins and human proteins were connected to a network of human protein-protein interactions taken from the Human Protein Reference Database (HPRD, Release 7) [54] and the Biological General Repository for Interaction Datasets (BioGRID) database [55]. We then compared the distribution of degree (number of interactions) and betweenness centrality [13] for the viral targets against all other proteins in the human networks with the Kolmogorov-Smirnov test in R [56]. P-Values were again corrected for multiple testing using the FDR-method by Benjamini and Hochberg. Degree and betweenness centrality are alternative measures of network centrality for individual proteins. High degree characterizes so-called hubs which are highly interactive while high betweenness centrality characterizes so-called bottlenecks which are central to many connections between proteins.
Supporting Information
Zdroje
1. PerlmanSNetlandJ 2009 Coronaviruses post-SARS: update on replication and pathogenesis. Nat Rev Microbiol 7 439 450
2. PeirisJSYuenKYOsterhausADStohrK 2003 The severe acute respiratory syndrome. N Engl J Med 349 2431 2441
3. WangLFShiZZhangSFieldHDaszakP 2006 Review of bats and SARS. Emerg Infect Dis 12 1834 1840
4. DrexlerJFGloza-RauschFGlendeJCormanVMMuthD 2010 Genomic characterization of severe acute respiratory syndrome-related coronavirus in European bats and classification of coronaviruses based on partial RNA-dependent RNA polymerase gene sequences. J Virol 84 11336 11349
5. WooPCLauSKHuangYYuenKY 2009 Coronavirus diversity, phylogeny and interspecies jumping. Exp Biol Med (Maywood) 234 1117 1127
6. von BrunnATeepeCSimpsonJCPepperkokRFriedelCC 2007 Analysis of intraviral protein-protein interactions of the SARS coronavirus ORFeome. PLoS One 2 e459
7. UetzPDongYAZeretzkeCAtzlerCBaikerA 2006 Herpesviral protein networks and their interaction with the human proteome. Science 311 239 242
8. FeskeSGiltnaneJDolmetschRStaudtLMRaoA 2001 Gene regulation mediated by calcium signals in T lymphocytes. Nat Immunol 2 316 324
9. HoganPGChenLNardoneJRaoA 2003 Transcriptional regulation by calcium, calcineurin, and NFAT. Genes Dev 17 2205 2232
10. BraunPTasanMDrezeMBarrios-RodilesMLemmensI 2009 An experimentally derived confidence score for binary protein-protein interactions. Nat Methods 6 91 97
11. AlexaARahnenfuhrerJLengauerT 2006 Improved scoring of functional groups from gene expression data by decorrelating GO graph structure 3. Bioinformatics 22 1600 1607
12. BauerSGrossmannSVingronMRobinsonPN 2008 Ontologizer 2.0—a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics 24 1650 1651
13. YuHKimPMSprecherETrifonovVGersteinM 2007 The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput Biol 3 e59
14. Hermann-KleiterNBaierG 2010 NFAT pulls the strings during CD4+ T helper cell effector functions. Blood 115 2989 2997
15. LubanJBossoltKLFrankeEKKalpanaGVGoffSP 1993 Human immunodeficiency virus type 1 Gag protein binds to cyclophilins A and B. Cell 73 1067 1078
16. WatashiKIshiiNHijikataMInoueDMurataT 2005 Cyclophilin B is a functional regulator of hepatitis C virus RNA polymerase. Mol Cell 19 111 122
17. MacianF 2005 NFAT proteins: key regulators of T-cell development and function 7. Nat Rev Immunol 5 472 484
18. PfefferleSKraehlingVDittVGrywnaKMuehlbergerE 2009 Reverse genetic characterization of the natural genomic deletion in SARS-Coronavirus strain Frankfurt-1 open reading frame 7b reveals an attenuating function of the 7b protein in-vitro and in-vivo. Virol J 6 131
19. SpiegelMWeberF 2006 Inhibition of cytokine gene expression and induction of chemokine genes in non-lymphatic cells infected with SARS coronavirus. Virol J 3 17
20. JainJLohCRaoA 1995 Transcriptional regulation of the IL-2 gene. Curr Opin Immunol 7 333 342
21. YoneyamaMSuharaWFukuharaYFukudaMNishidaE 1998 Direct triggering of the type I interferon system by virus infection: activation of a transcription factor complex containing IRF-3 and CBP/p300. EMBO J 17 1087 1095
22. WangYShiHRigoletPWuNZhuL 2010 Nsp1 proteins of group I and SARS coronaviruses share structural and functional similarities. Infect Genet Evol 10 919 924
23. ChenLGuiCLuoXYangQGuntherS 2005 Cinanserin is an inhibitor of the 3C-like proteinase of severe acute respiratory syndrome coronavirus and strongly reduces virus replication in vitro. J Virol 79 7095 7103
24. AlmazanFDediegoMLGalanCEscorsDAlvarezE 2006 Construction of a severe acute respiratory syndrome coronavirus infectious cDNA clone and a replicon to study coronavirus RNA synthesis. J Virol 80 10900 10906
25. YountBCurtisKMFritzEAHensleyLEJahrlingPB 2003 Reverse genetics with a full-length infectious cDNA of severe acute respiratory syndrome coronavirus. Proc Natl Acad Sci U S A 100 12995 13000
26. ErikssonKKMakiaDThielV 2008 Generation of recombinant coronaviruses using vaccinia virus as the cloning vector and stable cell lines containing coronaviral replicon RNAs. Methods Mol Biol 454 237 254
27. PanJPengXGaoYLiZLuX 2008 Genome-wide analysis of protein-protein interactions and involvement of viral proteins in SARS-CoV replication. PLoS One 3 e3299
28. ImbertISnijderEJDimitrovaMGuillemotJCLecineP 2008 The SARS-Coronavirus PLnc domain of nsp3 as a replication/transcription scaffolding protein. Virus Res 133 136 148
29. NgLFHibberdMLOoiEETangKFNeoSY 2004 A human in vitro model system for investigating genome-wide host responses to SARS coronavirus infection. BMC Infect Dis 4 34
30. PangRTPoonTCChanKCLeeNLChiuRW 2006 Serum proteomic fingerprints of adult patients with severe acute respiratory syndrome. Clin Chem 52 421 429
31. NeumanBWJosephJSSaikatenduKSSerranoPChatterjeeA 2008 Proteomics analysis unravels the functional repertoire of coronavirus nonstructural protein 3 2. J Virol 82 5279 5294
32. WalhoutAJVidalM 2001 High-throughput yeast two-hybrid assays for large-scale protein interaction mapping. Methods 24 297 306
33. WalhoutAJTempleGFBraschMAHartleyJLLorsonMA 2000 GATEWAY recombinational cloning: application to the cloning of large numbers of open reading frames or ORFeomes. Methods Enzymol 328 575 592
34. Barrios-RodilesMBrownKROzdamarBBoseRLiuZ 2005 High-throughput mapping of a dynamic signaling network in mammalian cells. Science 307 1621 1625
35. ZustRCervantes-BarraganLKuriTBlakqoriGWeberF 2007 Coronavirus non-structural protein 1 is a major pathogenicity factor: implications for the rational design of coronavirus vaccines. PLoS Pathog 3 e109
36. KamitaniWHuangCNarayananKLokugamageKGMakinoS 2009 A two-pronged strategy to suppress host protein synthesis by SARS coronavirus Nsp1 protein. Nat Struct Mol Biol 16 1134 1140
37. WatheletMGOrrMFriemanMBBaricRS 2007 Severe acute respiratory syndrome coronavirus evades antiviral signaling: role of nsp1 and rational design of an attenuated strain. J Virol 81 11620 11633
38. NarayananKHuangCLokugamageKKamitaniWIkegamiT 2008 Severe acute respiratory syndrome coronavirus nsp1 suppresses host gene expression, including that of type I interferon, in infected cells. J Virol 82 4471 4479
39. JonesBMMaESPeirisJSWongPCHoJC 2004 Prolonged disturbances of in vitro cytokine production in patients with severe acute respiratory syndrome (SARS) treated with ribavirin and steroids. Clin Exp Immunol 135 467 473
40. LoutfyMRBlattLMSiminovitchKAWardSWolffB 2003 Interferon alfacon-1 plus corticosteroids in severe acute respiratory syndrome: a preliminary study. JAMA 290 3222 3228
41. PeirisJSChuCMChengVCChanKSHungIF 2003 Clinical progression and viral load in a community outbreak of coronavirus-associated SARS pneumonia: a prospective study. Lancet 361 1767 1772
42. LiuJFarmerJDJrLaneWSFriedmanJWeissmanI 1991 Calcineurin is a common target of cyclophilin-cyclosporin A and FKBP-FK506 complexes. Cell 66 807 815
43. GalatABuaJ 2010 Molecular aspects of cyclophilins mediating therapeutic actions of their ligands. Cell Mol Life Sci 67 3467 3488
44. LuoCLuoHZhengSGuiCYueL 2004 Nucleocapsid protein of SARS coronavirus tightly binds to human cyclophilin A. Biochem Biophys Res Commun 321 557 565
45. FischerGGallayPHopkinsS 2010 Cyclophilin inhibitors for the treatment of HCV infection. Curr Opin Investig Drugs 11 911 918
46. BrockwaySMDenisonMR 2005 Mutagenesis of the murine hepatitis virus nsp1-coding region identifies residues important for protein processing, viral RNA synthesis, and viral replication. Virology 340 209 223
47. AlbersMKranzHKoberIKaiserCKlinkM 2005 Automated yeast two-hybrid screening for nuclear receptor-interacting proteins. Mol Cell Proteomics 4 205 213
48. DrostenCPreiserWGuntherSSchmitzHDoerrHW 2003 Severe acute respiratory syndrome: identification of the etiological agent. Trends Mol Med 9 325 327
49. HerzogPDrostenCMullerMA 2008 Plaque assay for human coronavirus NL63 using human colon carcinoma cells. Virol J 5 138
50. KremplCBallesterosMLZimmerGEnjuanesLKlenkHD 2000 Characterization of the sialic acid binding activity of transmissible gastroenteritis coronavirus by analysis of haemagglutination-deficient mutants. J GenVirol 81 489 496
51. PfefferleSKrahlingVDittVGrywnaKMuhlbergerE 2009 Reverse genetic characterization of the natural genomic deletion in SARS-Coronavirus strain Frankfurt-1 open reading frame 7b reveals an attenuating function of the 7b protein in-vitro and in-vivo. Virol J 6 131
52. BenjaminiYHochbergY 1995 Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Stat Methodol 57 289 300
53. RueppABraunerBDunger-KaltenbachIFrishmanGMontroneC 2008 CORUM: the comprehensive resource of mammalian protein complexes 1. Nucleic Acids Res 36 D646 D650
54. PeriSNavarroJDKristiansenTZAmanchyRSurendranathV 2004 Human protein reference database as a discovery resource for proteomics. Nucleic Acids Res 32 D497 D501
55. BreitkreutzBJStarkCRegulyTBoucherLBreitkreutzA 2008 The BioGRID Interaction Database: 2008 update. Nucleic Acids Res 36 D637 D640
56. Team RDC 2007 R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria
57. MarraMAJonesSJAstellCRHoltRABrooks-WilsonA 2003 The Genome sequence of the SARS-associated coronavirus. Science 300 1399 1404
Štítky
Hygiena a epidemiologie Infekční lékařství LaboratořČlánek vyšel v časopise
PLOS Pathogens
2011 Číslo 10
Nejčtenější v tomto čísle
- Severe Acute Respiratory Syndrome Coronavirus Envelope Protein Regulates Cell Stress Response and Apoptosis
- The SARS-Coronavirus-Host Interactome: Identification of Cyclophilins as Target for Pan-Coronavirus Inhibitors
- Biochemical and Structural Insights into the Mechanisms of SARS Coronavirus RNA Ribose 2′-O-Methylation by nsp16/nsp10 Protein Complex
- Evolutionarily Divergent, Unstable Filamentous Actin Is Essential for Gliding Motility in Apicomplexan Parasites