PMC:7321036 JSONTXT

The Global Phosphorylation Landscape of SARS-CoV-2 Infection Abstract The causative agent of the coronavirus disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has infected millions and killed hundreds of thousands of people worldwide, highlighting an urgent need to develop antiviral therapies. Here we present a quantitative mass spectrometry-based phosphoproteomics survey of SARS-CoV-2 infection in Vero E6 cells, revealing dramatic rewiring of phosphorylation on host and viral proteins. SARS-CoV-2 infection promoted casein kinase II (CK2) and p38 MAPK activation, production of diverse cytokines, and shutdown of mitotic kinases, resulting in cell cycle arrest. Infection also stimulated a marked induction of CK2-containing filopodial protrusions possessing budding viral particles. Eighty-seven drugs and compounds were identified by mapping global phosphorylation profiles to dysregulated kinases and pathways. We found pharmacologic inhibition of the p38, CK2, CDK, AXL, and PIKFYVE kinases to possess antiviral efficacy, representing potential COVID-19 therapies. Graphical Abstract Highlights • Phosphoproteomics analysis of SARS-CoV-2-infected cells uncovers signaling rewiring • Infection promotes host p38 MAPK cascade activity and shutdown of mitotic kinases • Infection stimulates CK2-containing filopodial protrusions with budding virus • Kinase activity analysis identifies potent antiviral drugs and compounds Phosphoproteomics analysis of SARS-CoV-2-infected Vero E6 cells reveals host cellular pathways hijacked by viral infection, leading to the identification of small molecules that target dysregulated pathways and elicit potent antiviral efficacy. Introduction Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an enveloped positive-sense RNA virus that belongs to the lineage B Betacoronavirus family. It is closely related to SARS-CoV, the causative agent of SARS, which emerged in the human population in 2002 (79% genetic similarity), and several SARS-related coronaviruses that circulate in bats (up to 98% genetic similarity) (Lai et al., 2020; Zhou et al., 2020). The pathophysiology of severe coronavirus disease 2019 (COVID-19) is similar to that of severe disease caused by SARS-CoV and is characterized by acute respiratory distress and excessive inflammation capable of inducing respiratory failure, multi-organ failure, and death (Wong et al., 2004; Zhang et al., 2020). To enter host cells, the SARS-CoV-2 spike (S) protein binds to an ACE2 receptor on the target cell and is subsequently primed by a serine protease, TMPRSS2, that cleaves the S protein and allows fusion of viral and lysosomal membranes (Hoffmann et al., 2020). Following entry, viral genomic RNA is translated to produce the polyproteins ORF1a and ORF1ab, which are subsequently cleaved by viral proteases into non-structural proteins that form the viral replication/transcription complex (RTC). Extensive remodeling of the host endoplasmic reticulum leads to formation of double-membrane vesicles, within which viral RNA synthesis occurs. The viral RNA genome is replicated by transcription of the negative-strand genomic RNA template, whereas subgenomic mRNAs are transcribed and translated to produce structural and accessory proteins. Structural proteins and viral genomes assemble at the endoplasmic reticulum-Golgi intermediate compartment (ERGIC) (de Haan and Rottier, 2005), followed by transport to the cell surface for release by exocytosis (Fehr and Perlman 2015). Although much about SARS-CoV-2 biology can be inferred based on similarity to SARS-CoV, SARS-CoV-2 is a novel coronavirus with unique properties that contribute to its pandemic-scale spread. Unlike SARS-CoV, SARS-CoV-2 infection is commonly asymptomatic, particularly in the younger population (Guan et al., 2020), and contagious prior to symptom onset (Rothe et al., 2020; Peiris et al., 2003; Bai et al., 2020). These characteristics contribute to the difficulty of containing SARS-CoV-2 spread through public health strategies and amplify the need to develop vaccines and therapies to protect against and treat COVID-19. Clinical management of COVID-19 is largely limited to infection prevention and supportive care. So far, remdesivir, a broad-spectrum antiviral agent, is the only medication approved for emergency use to treat COVID-19 by the US Food and Drug Administration (FDA) (Grein et al., 2020). Although the evidence supporting remdesivir use in patients with advanced COVID-19 is promising, there remains an urgent need for potent SARS-CoV-2 therapeutic agents, especially those that could be given in an outpatient setting, to effectively combat the COVID-19 pandemic. Proteomics approaches that globally quantify changes in protein abundance and phosphorylation are powerful tools to elucidate mechanisms of viral pathogenesis by providing a snapshot of how cellular pathways and processes are rewired upon infection (Johnson et al., 2020). Importantly, the functional outcomes of many phosphorylation events are well annotated, especially for kinases where phosphorylation directly regulates their activity. State-of-the-art bioinformatics approaches can then be employed to readily identify regulated kinases from phosphorylation profiles, many of which are likely to be established drug targets with therapeutic potential (Ochoa et al., 2016, 2020). Here we present a quantitative survey of the global phosphorylation and protein abundance landscape of SARS-CoV-2 infection, map phosphorylation changes to disrupted kinases and pathways, and use these profiles to rapidly prioritize drugs and compounds with the potential to treat SARS-CoV-2 infection. Results Phosphorylation Signaling Represents a Primary Host Response to SARS-CoV-2 Infection To determine how SARS-CoV-2 hijacks host-protein signaling, a global phosphoproteomics experiment was performed in Vero E6 cells, a cell line originating from the kidney of a female African green monkey (Chlorocebus sabaeus) (Osada et al., 2014). This cell line was selected because of its high susceptibility to SARS-CoV-2 infection (Harcourt et al., 2020). Cells were harvested in biological triplicate at 6 time points after SARS-CoV-2 infection (0, 2, 4, 8, 12, or 24 h) or after mock infection at 0 or 24 h (Figure 1 A). Using a data-independent acquisition (DIA) proteomics approach, each sample was then partitioned and analyzed for changes in global protein abundance or phosphorylation (data available in Table S1). Chlorocebus sabaeus and human protein sequences were aligned, and phosphorylation sites and protein identifiers were mapped to their respective human protein orthologs. Phosphorylation fold changes calculated using the 0- or 24-h mock control were highly comparable (correlation coefficient r = 0.77); therefore, the 0-h mock control was used for all subsequent comparisons. Figure 1 Global Proteomics of Phosphorylation and Abundance Changes upon SARS-CoV-2 Infection (A) Vero E6 cells were infected with SARS-CoV-2 (MOI 1.0). After 1 h of viral uptake, cells were harvested (0 h) or, subsequently, after 2, 4, 8, 12, or 24 h. As a control, Vero E6 cells were also mock infected for 1 h and harvested immediately thereafter (0 h) or after 24 h of mock infection. All conditions were performed in biological triplicate. Following cell harvest, cells were lysed, and proteins were digested into peptides. Aliquots of all samples were analyzed by mass spectrometry (MS) to measure changes in protein abundance upon infection, whereas the remaining sample was enriched for phosphorylated peptides and subsequently analyzed to measure changes in phosphorylation signaling. A DIA approach was used for all MS acquisitions. Last, all phosphorylation sites and protein identifiers were mapped to their respective human protein orthologs. (B) Principal-component analysis (PCA) of phosphorylation replicates after removing outliers. See also Figure S1. (C) Correlation of protein abundance (AB) and phosphorylation sites (PHs) between replicates within a biological condition (red) and across biological conditions (black). Boxplots depict median (horizonal lines), interquartile range (boxes), maximum and minimum values (vertical lines), and outliers (solid circles). (D) Median AB of individual SARS-CoV-2 proteins in the protein AB analysis. (E) The number of significantly regulated PH groups across the infection time course. (F) Volcano plot of PH group quantification 24 h after infection. (G) The number of significantly regulated proteins across the infection time course. (H) Volcano plot of protein AB quantification 24 h after infection. (I) Gene Ontology enrichment analysis of all significantly changing proteins in terms of AB divided into two sets: downregulated (blue) and upregulated (red). (J) Proportion of significantly regulated PH groups with a correlated (i.e., same direction, AB match) or anticorrelated (i.e., opposite direction, AB mismatch) significant or insignificant (gray) change in protein AB. In (E)–(H), all infection time points are compared with the mock infection at 0 h, and significantly regulated proteins are defined as (absolute value of log2(infection/mock) > 1 and adjusted p < 0.05 or when only detected in infected or mock based on replicate and MS feature counts; STAR Methods). See also Figure S1. Quality control filtering of the data was performed, and two samples from each of the phosphorylation and protein abundance datasets were removed because of poor correlation with their respective replicates (Figures S1 A and S1B). Principal component analysis (PCA) of the remaining samples revealed good separation of mock and infected samples as well as high quantitative reproducibility between biological replicates (Figures 1B, 1C, and S1C). In total, high-quality quantification of 4,624 human-orthologous phosphorylation sites and 3,036 human-orthologous proteins was obtained (Figure S1D). Successful infection was confirmed by the observation of a dramatic increase in viral protein abundance over the course of a 24-h infection period (Figures 1D and S1E). Figure S1 Proteomics: Quality Control (QC), Orthology, Enrichments, and Viral Proteins, Related to Figure 1 (A) Principal component analysis computed on intensities summarized by MSstats at the level of phosphorylation site groups within (from left to right) all runs, with one outlier run removed, and with two outlier runs removed. Outlier runs are labeled 00Hr.2 and 02.Hr.2. (B) Principal components analysis computed on protein intensities as summarized by MSstats (from left to right) within all runs, with one outlier run removed, and with two outlier runs removed. Outlier runs are both labeled 00Hr.2; one is mock and the other is infected. (C) Coefficient of variance boxplot for each condition. Black lines depict the median and their values are indicated above each boxplot. (D) Mapping detected and quantifiable proteins and phosphorylation sites from the green monkey (Chlorocebus sabaeus) protein sequences to human genes. Proteins and sites were considered quantifiable if MSstats computed a non-infinite fold change for any time point or if an infinite log2 fold change passes criteria for inclusion in any time point. (E) Intensities of viral proteins as summarized over all peptide ion fragments by MSstats, averaged across replicates. The MSstats summarization is based on the median intensity of all fragments after data pre-processing (STAR Methods). (F) Gene Ontology enrichment analysis for proteins significantly regulated in terms of abundance upon infection, separated by time point and direction of phosphorylation regulation. All terms with significant over-representation (adjusted p value < 0.01) in the regulated gene set are kept, and redundant terms are removed (see STAR Methods). Numbers in cells indicate the number of genes that match the term for a given time point and direction. (G) Gene Ontology enrichment analysis for significantly phosphorylated proteins upon infection, separated by time point and direction of protein regulation. Details same as for (F). As expected, an increase was observed in the number of significantly regulated phosphorylation sites and proteins over the infection time course, with the majority of regulation occurring at the level of phosphorylation (Figures 1E and 1F) as opposed to protein abundance (Figures 1G and 1H). Of the few proteins that significantly increased in abundance upon infection, the vast majority were SARS-CoV-2 viral proteins (Figure 1H). In contrast, the majority of host proteins decreased in abundance. This finding is consistent with mechanisms of host mRNA nuclear export and/or host mRNA translation inhibition, which are common in viral infections (Kuss et al., 2013; Walsh and Mohr 2011). Gene Ontology enrichment analysis of significantly downregulated proteins revealed several terms related to platelet regulation (Figures 1I and S1F). Several downregulated host proteins are known to be involved in platelet regulation, thrombosis, and prevention of blood coagulation, including APOH, CD9, TSPAN14, AHSG, SERPINA1, and A2M (Mather et al., 2016; Mangin et al., 2009; Taggart et al., 2000). The downregulation of these proteins suggests that they may mechanistically contribute to symptoms of blood coagulation and stroke in COVID-19 patients (Han et al., 2020). Lastly, the contribution of protein abundance to phosphorylation level changes was evaluated. For nearly all cases of a significantly changed phosphorylation site, no corresponding significant change in protein abundance was observed (Figure 1J), further suggesting that phosphorylation signaling represents a primary host response over this time course of infection as opposed to transcriptional regulation, which would influence protein abundance. Phosphorylation of SARS-CoV-2 Viral Proteins by the Host Proteome Viral protein phosphorylation within the host cell may play a role in sensing and responding to cell state. We detected 25 phosphorylation sites in SARS-CoV-2 viral proteins that we combined with another proteomics dataset (Davidson et al., 2020) to amass a total of 49 sites detected across seven viral proteins (Table S2). Of note, this analysis does not distinguish cleaved from uncleaved viral proteins in the assignment of viral phosphorylation sites. The degree of conservation, indicative of functional constraint, was estimated for each residue position (Figure 2 A; Ng and Henikoff 2003), and the sites were mapped to positions within structured regions for five proteins, with the majority observed in accessible positions (i.e., loops) (Figure 2B). The top kinase families predicted by sequence to regulate these sites included casein kinase II (CK2), cyclin-dependent kinase (CDK), and protein kinase C (PKC), among others (Figure 2C), suggesting that these kinases may contribute to regulation of viral replication. Figure 2 Overview of SARS-CoV-2 Viral Protein PHs in the Host Cell (A) Localization of PHs across viral protein sequences from this study and a previous study (Davidson et al., 2020). Stem height indicates predicted deleteriousness of alanine substitutions. Dot color indicates whether the residue is (true) or is not (false) predicted to form part of an interaction interface based on SPPIDER analysis. Positions with no structural coverage are excluded from interface prediction. (B) Distribution of secondary structure elements in which viral PHs were found, as classified by the define secondary structure of proteins (DSSP) tool. (C) Distribution of top matching host kinases to viral PHs according to NetPhorest tool (Horn et al., 2014). (D) Phosphorylation cluster in the C-terminal tail of the M protein (red residues) structure (Heo and Feig 2020) and associated sequence motif. Asterisks indicate PHs. (E) Alignment of M protein phosphorylation clusters across different coronaviruses. Asterisks indicate PHs. (F) Surface electrostatic potential of non-phosphorylated (left) and phosphorylated (right) RNA-binding domains of the N protein (PDB: 6M3M). Positions of PHs are indicated by arrows. Blue denotes a positive charge potential, and red indicates a negative charge potential. Electrostatic potential was computed with the Advanced Poisson-Boltzmann Solver (APBS) tool after preparation with the PDB2PQR tool. Although it is unlikely that all phosphorylation sites on viral proteins play important functional roles, several sites in membrane (M) protein, Nsp9, and nucleocapsid (N) protein (Figures 2D–2F) suggest potential functionality. Five phosphorylation sites were detected in the M protein cluster within a short C-terminal region of the protein (207–215; Figure 2D). Although these acceptor residues are not predicted to be conserved, several are negatively charged residues in M proteins of other related viruses (Figure 2E). This evolutionary pattern suggests that a negative charge in this region may play a functional role, reminiscent of other multi-site phosphorylation events (Serber and Ferrell 2007). To identify phosphorylation sites that may regulate protein-protein interactions, all sites were mapped to 3D structures, and solvent accessibility based protein-protein interface identification and recognition (SPPIDER) was used to assess whether sites resided within interface regions (Porollo and Meller 2007; Figure 2A; Table S2). The single phosphorylation site in Nsp9 was predicted to be at an interface region (“True”), which was supported by inspection of the homodimer structure (PDB: 6W4B). Additional phosphorylation sites were predicted to be at interface residues within the S protein (Figure 2A). However, inspection of S in complex with the ACE2 receptor (Shang et al., 2020; Lan et al., 2020) reveals some of these phosphorylation sites to be near but not at the interface region. Finally, phosphorylation sites in N protein, a structural protein that binds to and assists with packaging viral RNA, were investigated. Most sites occurred within the N-terminal portion of the protein, at or near the RNA binding region, but avoided the C-terminal dimerization domain. The cluster of phosphorylation sites within an arginine/serine (RS)-dipeptide rich region, C-terminal to the RNA binding region (Figure 2A), is conserved in other coronavirus N proteins. This region is phosphorylated in SARS-CoV by serine-arginine (SR) protein kinases, modulating the role of SARS-CoV N protein in host translation inhibition (Peng et al., 2008). It is likely that phosphorylation of this same region in SARS-CoV-2 plays a similar role. Interestingly, in vitro inhibition of SARS-CoV N protein phosphorylation at the RS-rich region results in reduced viral load and cytopathic effects (Wu et al., 2009), highlighting its importance for viral fitness. In addition, sites spanning the sequence of the RNA binding domain, which forms a claw-like structure, have been observed (Kang et al., 2020). Several phosphorylation sites cluster in the structural model, predicted to affect the surface charge of the so-called acidic wrist region (Figure 2F) but not the positive surface charge of the RNA binding pocket. We hypothesize that this surface charge difference may modulate N protein function, potentially via allosteric regulation of RNA binding capacity. Phosphorylation of SARS-CoV-2 Host-Interacting Proteins during Infection The recently published SARS-CoV-2 virus-human protein-protein interaction map identified 332 human proteins interacting with 27 (26 wild-type and 1 mutant) viral proteins (Gordon et al., 2020). Here we found some of these host proteins (40 of 332) to be significantly differentially phosphorylated upon infection (Figure 3 ). Virus-host protein-protein interactions could drive changes in phosphorylation by affecting host protein subcellular localization or by sterically blocking kinase access. Furthermore, phosphorylation of these proteins upon infection may signify an additional mode of functional control over these potential dependency and restriction factors. Figure 3 Phosphorylation on SARS-CoV-2 Virus-Human Interacting Proteins The SARS-CoV-2 virus-host protein-protein interaction map (Gordon et al., 2020) found 332 human proteins interacting with 27 (26 wild-type and 1 mutant) viral proteins. Here we found 40 of 332 proteins significantly differentially phosphorylated across at least two time points (adjusted p < 0.05 and absolute value of log2 fold change [abs(log2FC)] > 1). Viral proteins are shown as red diamonds. Interacting human proteins are shown as gray circles. PHs emanate from human proteins, colored by their log2 fold change compared with uninfected control samples (red, increase; blue, decrease) at each time point (0, 2, 4, 8, 12, and 24 h after infection) in a clockwise fashion. An interactive version of phosphorylation data can be found at https://kroganlab.ucsf.edu/network-maps. The SARS-CoV-2 N protein is known to interact with several RNA-processing proteins that are differentially phosphorylated during infection, including LARP1 and RRP9. Here LARP1 phosphorylation decreases on several sites, which is known to consequently increase LARP1 affinity for 3′ untranslated regions (UTRs) of mRNAs encoding ribosomal proteins, driving inhibition of protein synthesis (Hong et al., 2017). This mechanism may be utilized by SARS-CoV-2 to prioritize synthesis of viral proteins over host proteins. In addition, ORF6 interacts with the NUP98/RAE complex, and NUP98 phosphorylation was observed to increase at S888, a site within its peptidase domain. NUP98 autocatalytic cleavage is required for localization to the nuclear pore; thus, it is possible that NUP98 interaction with ORF6 and/or its virus-induced phosphorylation prevents host mRNA export through the nuclear pore (Krull et al., 2010; Hodel et al., 2002). A similar mechanism is employed by vesicular stomatitis virus (VSV) matrix protein to block host mRNA export by targeting the NUP98/RAE complex, leading to exclusive translation of cytoplasmic VSV mRNAs (Quan et al., 2014). For Nsp12, the majority of its protein interactors displayed decreased phosphorylation during infection. Because Nsp12 is known to encode the RNA-dependent RNA polymerase, responsible for replicating the viral genome, and several of these interacting proteins are related to RNA processing (LARP4B and CRTC3), their regulation may possess functional implications for Nsp12 in viral RNA replication. In addition, Nsp8 interacts with several proteins whose phosphorylation increases (LARP7 and MPHOSPH10) and decreases (CCDC86) on several sites. Notably, LARP7 and MEPCE are important regulators of RNA polymerase II-mediated transcription elongation as part of the 7SK small nuclear ribonucleoprotein particle (snRNP) complex. Regulation of these phosphorylation sites may contribute to the regulation of positive transcription elongation factor b (P-TEFb [CDK9]) and transcriptional regulation of the virus, similar to how these proteins are regulated during HIV infection (Mbonye et al., 2015). SARS-CoV-2 Infection Regulates Host Kinase Signaling To study global changes in kinase signaling and their effect on host protein phosphorylation, regulated phosphorylation sites were grouped in five clusters based on their dynamics using a data-driven clustering approach (Figure 4 A; STAR Methods). For each of the groups, an enrichment analysis was performed for functions and pathways (Figure 4A; Table S3). The dynamics of these changes can be linked to the viral life cycle: entry (0–2 h), replication (4–12 h), and egress (24 h). Clusters 1 and 2 include phosphorylation sites that are, on average, upregulated during infection. Cluster 1 sites tended to be upregulated within 2 h (i.e., linked to viral entry) and were enriched in mRNA processing, cell cycle, apoptosis, and proteins involved in HIV infection. Cluster 2 included apoptosis proteins with later onset of phosphorylation, associated with replication and/or egress. Phosphorylation sites in clusters 3 and 4 were downregulated and enriched in RNA-processing functions. Sites within cluster 5 possessed a dynamic response to infection, with immediate downregulation followed by a rise during the middle and renewed downregulation at late time points. This cluster was enriched for DNA replication and the cell cycle, among others. These observations are corroborated by standard Gene Ontology (GO) enrichment analyses of biological processes regulated by phosphorylation (Figure S1G; Table S3; STAR Methods). Figure 4 Signaling Changes in Host Cells in Response to SARS-CoV-2 Infection (A) Clusters of significantly changing PHs (abs(log2FC) > 1 and adjusted p < 0.05) across the time course of infection with non-redundant enriched Reactome pathway terms (adjusted p value [q] < 0.01) for each cluster. Horizontal red lines below each pathway term correspond to phosphorylated proteins belonging to the pathway, and a black-bordered rectangle is indicative of a significantly enriched term. (B) Kinases depicting a strong change in activity upon infection (abs(log10(p)) > 2.5) in at least one time point, with predicted activity in at least 5 of 6 time points. (C) Schematic representation of interaction between host kinases and SARS-CoV-2 viral proteins from Gordon et al. (2020). Substrate PHs for each kinase are color-coded as blue (down) and red (up) based on the direction of change during infection. Only PHs corresponding to the kinase activity direction are shown. (D) Correlation of kinase activity profiles of each time point with other biological conditions with at least one significantly changing kinase (abs(log10(p)) > 2.5) and having significant correlation with at least one time-point (false discovery rate [FDR] < 5%). (E) Overall phosphorylation change (−log10(p)) of a protein complex, estimated as the change in phosphorylation on member proteins. Only non-redundant protein complexes with a significant change in phosphorylation (abs(log10(p)) > 2.5) in at least one time point are shown. See also Figure S2. We estimated activity regulation for 97 kinases based on the regulation of their known substrates (Ochoa et al., 2016; Hernandez-Armenta et al., 2017; Table S4), with the strongest regulation linked to viral entry (0–2 h) and late replication/egress (24 h). The kinases predicted to be most strongly activated (Figures 4B and S2A) include several members of the p38 pathway, including p38ɣ (MAPK12), CK2 (CSNK2A1/2), Ca2+/calmodulin-dependent protein kinase (CAMK2G), and the guanosine monophosphate (GMP)-dependent protein kinases PRKG1/2, which can inhibit Rho signaling. Kinases predicted to be downregulated include several cell cycle kinases (CDK1/2/5 and AURKA), cell growth-related signaling pathway kinases (PRKACA, AKT1/2, MAPK1/3, and PIM1), and the cytoskeleton regulators (PAK1), among others. Kinase activity estimates based on the 24-h mock control gave highly correlated results (r = 0.81), identifying the same set of highly regulated kinases (Figure S2D). Some of the changes in kinase activity can be directly linked to host-viral protein interactions (Figure 4C). Among the 10 interacting kinases detected in a virus-host protein-protein interaction map (Gordon et al., 2020), an increase in activity for CK2 and a decrease for MARK2 and PRKACA were observed (Figure 4C). Of note, although we predict decreased activity for PRKACA based on phosphorylation of its substrates, we simultaneously detected a significant increase in T198 phosphorylation (8, 12, and 24 h post infection) within its activation loop, suggesting an increase in PRKACA activity. It is possible that Nsp13 is sequestering active PRKACA away from its typical substrates. Figure S2 Full Kinase Activities, Correlated Conditions, and Regulated Complexes, Related to Figure 4 (A) Changes in predicted kinase activities across different time points post-infection. (B) Correlation of kinase activity profiles of each time point with other biological conditions. Kinase activities were estimated for a wide-range of biological conditions obtained from previously published phosphoproteomics datasets (Ochoa et al., 2016). (C) Changes in phosphorylation in protein complexes. Overall phosphorylation change (-log10 (p value) of a protein complex was derived from change in phosphorylation of sites in member proteins. (D) Kinase activity estimates when using either the 0- or 24-h mock controls for those top regulated kinase activities from the 0-h control comparison. To better understand the signaling states of cells over the course of infection, we compared our data with a compilation of public phosphoproteomics datasets of other conditions (Ochoa et al., 2016; Figures 4D and S2B). The first and last time point of infection resembled a kinase activation state induced by inhibition of mTOR, ERK, AKT, and EGFR, consistent with the estimated kinase activities of these growth-related pathways. Between 2 and 12 h after infection, kinase activity states resembling inhibition of phosphatidylinositol 3-kinase (PI3K), p70S6K, and Rho-associated protein kinases (ROCKs) were observed. Finally, several of the time points resembled S/G2 cell cycle state, suggestive of a cell cycle block. Conversely, some conditions were anticorrelated with kinase activity profiles (Figures 4D and S2B). In line with an S/G2 cell-cycle block, infection signaling appeared opposite to that of a mitotic cell. In addition, inhibitors of histone deacetylases (HDACs) (scriptaid and trichostatin A), the proteasome (bortezomib), Hsp90 (geldanamycin), and voltage-gated sodium channels (valproic acid) were also anticorrelated. These drugs, or drugs targeting these protein activities, could induce a signaling state that inhibits viral replication. To further link kinase activities to downstream protein complexes, enrichment of up- or downregulated phosphorylation sites was determined within a curated set of human protein complexes defined by CORUM (Giurgiu et al., 2019; Figures 4E and S2C). This analysis revealed significant changes in phosphorylation of splicing related complexes (spliceosome), the proteasome (PA700-20S-PA28), and chromatin remodeling complexes (HuCHRAC and MLL2). In addition, a subset of regulated phosphorylation sites were detected that have known regulatory functions or high predicted functional scores (Ochoa et al., 2020) that are linked to regulation of protein activities (Table S5). Consistent with the observed signaling changes described above, these regulatory phosphorylation sites are involved in activation of chaperones (including HSP90), proteasome activity, inhibition of the anaphase-promoting complex (APC), and regulation of HDACs and cytoskeleton proteins, among others. CK2 and N Co-localize at Virus-Induced Filopodial Protrusions The phosphoproteomics data indicated regulation of several kinases and effector proteins related to cytoskeleton organization upon SARS-CoV-2 infection. Kinases downstream of the Rho/Rac/Cdc42 GTPases (PAK1/2 and ROCK1/2) and several well-characterized phosphorylation site targets of PAK1/2 kinase in vimentin (VIM S39 and S56) and stathmin (STMN1 S16 and S25) were found to be downregulated during infection (Figures 5A and 5B). The interaction of Nsp7 with RHOA (Gordon et al., 2020) may contribute to this downregulation. In contrast, signaling via CK2 is strongly upregulated, as determined by the increase in phosphorylation of well-characterized target sites (Figures 5A and 5B). Among the many roles of this kinase, we noted increased phosphorylation of cytoskeleton protein targets such as ɑ-Catenin (CTNNA1 S641) and the heavy chain of the motor protein Myosin IIa (MYH9 S1943). In addition to these kinase-mediated effects, the Nsp2 protein of SARS-CoV-2 also interacts directly with strumpellin (WASHC5), a subunit of the actin assembly-inducing WASH complex (Gordon et al., 2020), further implicating cytoskeleton regulation during infection. To study the relevance of these observations in a human infection model, high-resolution immunofluorescence imaging of fixed Caco-2 human colon epithelial cells was performed 24 h after infection (STAR Methods). Figure 5 Colocalization of CK2 and Viral Proteins at Actin Protrusions (A) Pathway of regulated PHs and SARS-CoV-2 interaction partners involved in cytoskeletal reorganization. Dashed lines indicate downregulation of activity, while solid lines indicate upregulation of activity. (B) Regulation of individual kinase activity or PHs depicted in (A). (C) Caco-2 cells infected with SARS-CoV-2 at an MOI of 0.1 for 24 h prior to immunostaining for F-actin and M protein, as indicated. Shown is a confocal section revealing M protein localization along and to the tip of filopodia (left) and magnification of the dashed box (right). (D) Dot plot quantification of the number and length of filopodia in untreated (mock) or infected Caco-2 cells for 24 h with SARS-CoV-2. Filopodium length was measured from the cortical actin to the tip of the filopodium. Error bars represent SD. Statistical testing by Mann-Whitney test. (E) Caco-2 cells infected with SARS-CoV-2 at an MOI of 0.01 for 24 h prior to immunostaining for F-actin, N protein, and casein kinase II (CK2) as indicated (left). Shown is magnification of the dashed box as single channels (right). (F) Magnification of the dashed box from (E) with quantification of colocalization between CK2 and N protein throughout infected Caco-2 cells. Displayed is the proportion of N protein-positive particles colocalizing with CK2. Error bars represent SD. (G and H) Scanning electron microscopy (G) and transmission electron microscopy (H) images of SARS-CoV-2 budding from Vero E6 cell filopodia (black arrow in H). See also Figure S3. SARS-CoV-2 infected Caco-2 cells were imaged for filamentous actin and the SARS-CoV-2 M protein, revealing prominent M protein clusters, possibly marking assembled SARS-CoV-2 viral particles, localized along the shafts and at the tips of actin-rich filopodia (Figures 5B and S3 B). SARS-CoV-2 infection induced a dramatic increase in filopodial protrusions, which were significantly longer and more branched than in uninfected cells (Figure 5D). Uninfected cells also exhibited filopodial protrusions, but their frequency and shape were dramatically different (Figure S3A). Reorganization of the actin cytoskeleton is a common feature of many viral infections and is associated with different stages of the viral life cycle (Taylor et al., 2011). Figure S3 Microscopy Images Showing Response to SARS-CoV-2 Infection, Related to Figure 5 (A) Non-infected Caco2 cells co-stained for F-actin, CK2 and nuclei (DAPI). Magnification of the indicated area is displayed as a single channel and merged images on the right panels. (B) Caco2 cells infected with SARS-CoV-2 at an MOI of 0.1 for 24 h prior to immunostaining for F-actin and M-protein, as indicated. See lower (1) and right (2) panel for magnification of regions indicated by dashed boxes. (C) Scanning electron microscopy and (D) transmission electron microscopy image of SARS-CoV-2 budding from Vero E6 cell filopodia. (E) N protein was found to physically interact with casein kinase II subunits (cartoon, left), CSNK2B and CSNK2A2 (Gordon et al., 2020). To test whether N protein could directly control CK2 activity, N protein was transduced via lentivirus in Vero E6 cells and stably induced via doxycycline for 48 hours followed by phosphoproteomics analysis. Kinase activities were calculated as before (STAR Methods) and top up- (> 1.5, red) and downregulated (< 1.5, blue) kinases are shown. See Table S1 for full phosphoproteomics data and Table S4 for full list of predicted kinase activities. We hypothesize that induction of virus-containing filopodia could be important for SARS-CoV-2 egress and/or cell-to-cell spread within epithelial monolayers. Given that Rho/PAK/ROCK signaling is downregulated, we next asked whether CK2 could play a role in this process. At 24 h, infected cells showed CK2 expression along the thin filopodial protrusions (Figure 5E), partially co-localized with SARS-CoV-2 N protein (Figure 5F). Scanning and transmission electron microscopy were used (Figures 5G, 5H, S3C, and S3D) to image the cellular protrusions at higher resolution. Assembled viral particles are clearly visible along these filopodia (Figure 5G), with instances where the viral particles appear to be budding from the protrusions (Figure 5H). Finally, we performed a global phosphoproteomics analysis of Vero E6 cells overexpressing N protein and observed CK2 activity to be significantly upregulated (Figure S3E; Tables S1 and S4). Because CK2 activity can promote actin polymerization (D’Amore et al., 2019), we hypothesize that N protein may allosterically control CK2 activity and regulate cytoskeleton organization. SARS-CoV-2 Infection Promotes p38/MAPK Signaling Activity and Cell Cycle Arrest Kinase activity analysis of SARS-CoV-2 phosphorylation profiles predicted upregulation of several components of the p38/mitogen-activated protein kinase (MAPK) signaling pathway, including MAP2K3, MAP2K6, MAPK12, MAPKAPK2 (MK2), and MAPKAPK3 (Figures 6A and 6B). Immunoblotting for activated phospho-p38 (T180/Y182), phospho-MK2 (T334), and phospho-cAMP response element-binding protein (CREB) and phospho-ATF-1 at their respective MAPKAPK2 sites (S133 in both) confirmed activation of the p38/MAPK pathway during SARS-CoV-2 infection in ACE2-expressing A549 human lung carcinoma cells (ACE2-A549) (Figure 6C). Furthermore, phosphoproteomics data depict increased phosphorylation of p38 pathway substrates such as negative elongation factor E (NELFE), heat shock protein beta-1 (HSPB1), and signal transducer and activator of transcription 1-alpha/beta (STAT1), among others (Figure 6D). Regulation of these sites occurs late in the time course (24 h after infection), likely reflecting a more advanced stage of viral infection, replication, and egress. Figure 6 SARS-CoV-2 Activates the p38/MAPK Signaling Pathway and Causes Cell Cycle Arrest (A) Diagram of the p38/MAPK signaling pathway. (B) Kinase activity analysis for kinases in the p38/MAPK pathway. (C) Western blot analysis of phosphorylated p38/MAPK signaling components in mock- and SARS-CoV-2-infected ACE2-A549 cells 24 h after infection. (D) Log2 fold change profiles of indicated p38/MAPK substrates during SARS-CoV-2 infection in Vero E6 cells. (E) Transcription factor activity analysis of SARS-CoV-2-infected A549, Calu-3, and NHBE cells, comparing p38/MAPK transcription factors with transcription factors not associated with the p38/MAPK pathway. Statistical test: Mann-Whitney test. (F) qRT-PCR analysis of the indicated mRNA from ACE2-A549 cells pre-treated with the p38 inhibitor SB203580 at the indicated concentrations for 1 h prior to infection with SARS-CoV-2 for 24 h. Statistical test: Student’s t test. See also Figures S4 and S5. (G) Heatmap of Pearson’s correlation coefficients comparing SARS-CoV-2-infected Vero E6 phosphorylation profiles with profiles of cells with induced DNA damage and cells arrested at the indicated cell cycle stages. (H) Log2 fold change profiles of the indicated cell cycle and DNA damage substrates during SARS-CoV-2 infection in Vero E6 cells. (I) DNA content analysis of cells infected with SARS-CoV-2 for 24 h compared with mock-infected cells. The p38/MAPK pathway mediates the cellular response to environmental stress, pathogenic infection, and pro-inflammatory cytokine stimulation, whereas downstream effectors of the pathway include transcription factors and RNA binding proteins that promote inflammatory cytokine production (Cuadrado and Nebreda 2010; Wen et al., 2010). Analysis of estimated transcription factor activity from gene expression data (STAR Methods; Table S6) derived from the infection of a human lung carcinoma cell line (A549), a human epithelial lung cancer cell line (Calu3), and primary human bronchial epithelial (NHBE) cells demonstrated that transcription factors regulated by the p38/MAPK pathway were among the most highly activated upon infection (Figure 6E; Blanco-Melo et al., 2020). To investigate the contribution of the p38/MAPK pathway to cytokine production, SARS-CoV-2-infected ACE2-A549 cells were treated with the p38 inhibitor SB203580. The mRNA of the inflammatory cytokines interleukin-6 (IL-6), tumor necrosis factor alpha (TNF-α), and others increased during infection and were inhibited by p38 inhibition in a dose-dependent manner (Figures 6F, right, and S4A). Interestingly, p38 inhibition also reduced SARS-CoV-2 subgenomic mRNA (Figure 6F, left) in the absence of major cellular toxicity (Figure S5 ), indicative of reduced viral replication. The SB203580-induced decrease in virus production was further confirmed using an anti-SARS-CoV-2 N protein (anti-NP) antibody-based assay (Figure S5, New York Vero E6). Multiplexed ELISA analysis of supernatants of cells from the same experiment demonstrated strong upregulation of inflammatory cytokines at the protein level, including IL-6, CXCL8, CCL20, and CCL2, which were decreased upon p38 inhibition (Figure S4 B; Table S7). However, because SARS-CoV-2 replication is also inhibited by SB203580, we cannot deconvolve the contributions of p38/MAPK pathway activity and SARS-CoV-2 virus presence on cytokine production. Figure S5 Pharmacological Profiling for Viral Titers and Cell Viability, Related to Figure 7 Dose response of phosphoproteomics-informed drugs and compounds. Assays performed in New York (N; red, anti-NP; blue TCID50) and Paris (P; red, RT-qPCR; purple, plaque assays) across two cell lines (A549-ACE2 and Vero E6). Cell viability shown in black. Mean of three biological replicates is shown. Error bars are SEM. Figure S4 Cytokine Profiling upon Infection, p38 Inhibition, and Cell Cycle Analysis, Related to Figure 6 (A) RT-qPCR analysis of indicated mRNA from A549-ACE2 cells pre-treated with p38 inhibitor SB203580 at indicated concentrations for one hour prior to infection with SARS-CoV-2 for 24 hours. Statistical test is Student’s t test. Error bars are SD. (B) Same as in (A) but a Luminex-based quantification of indicated cytokines. Error bars are SD. (C) Cell cycle analysis of Vero E6 cells (same as in Figure 6I) upon SARS-CoV-2 infection at an MOI of 1. Cell stained with DAPI DNA stain prior to flow cytometry analysis. Statistical test is Mann-Whitney test. Error bars are SD. Comparing phosphoproteomics profiles of SARS-CoV-2-infected cells with a database of phosphorylation profiles collected at specific cell cycle stages, viral infection was most highly correlated with cells arrested at the S/G2 transition and was negatively correlated with profiles of cells in mitosis (Figure 6G). We also observed SARS-CoV-2-dependent regulation of CDK2 T14/Y15 phosphorylation, initially increased in response to SARS-CoV-2 infection at 2 h, followed by a decrease over the remainder of the infection time course (Figure 6H, left). CDK2 activity promotes transition from the G2 phase of the cell cycle into mitosis and is inhibited by phosphorylation at positions T14 and Y15 by kinases WEE1 and MYT1, preventing premature entry into mitosis (Parker and Piwnica-Worms 1992; Mueller et al., 1995). CDK2 can also become phosphorylated when the cell cycle is arrested because of checkpoint failure or DNA damage. In addition, H2AX S140 phosphorylation (i.e., γ-H2AX), a hallmark of the DNA damage response, exhibited a profile similar to CDK2, suggesting that the DNA damage response may become activated early during infection (Rogakou et al., 1998; Figure 6H, right). To more directly test whether SARS-CoV-2 infection affects cell cycle progression, cells were infected with SARS-CoV-2 for 24 h, and their DNA content was measured using DAPI DNA staining and flow cytometry. A significant increase in the fraction of cells in S phase and at the G2/M transition and a decrease in the fraction of cells in G0/G1 phase were observed (Figures 6I and S4C). This observation is consistent with arrest between S and G2 phases of the cell cycle. A relationship between p38 activity and cell cycle arrest has been described previously, and the two could be linked mechanistically during SARS-CoV-2 infection (Lee et al., 2002; Yee et al., 2004). Mapping Kinase Activities to Pharmacological Modulators Identifies SARS-CoV-2 Therapies To identify effective therapies for SARS-CoV-2 infection, kinase inhibitors were mapped to the most differentially regulated kinase activities (Figure 7 A) and to specific phosphorylation sites (Table S8; STAR Methods). This resulted in a list of 87 drugs and compounds: 10 FDA-approved, 53 undergoing clinical testing, and 24 pre-clinical. Many of the drugs and compounds identified were reported to target several host kinases in cell-free assays at a minimum, but many have been observed to hit targets in cellular assays as well (Figure 7A). We reasoned that testing molecules with both overlapping and unique targets would help specify the molecular targets of greatest importance for SARS-CoV-2. Here, 68 total drugs and compounds were tested for antiviral efficacy (via qRT-PCR, anti-NP antibody, plaque assay, and/or TCID50) and cellular toxicity at two different institutions (in New York [Mount Sinai, 25 drugs/compounds] and Paris [Institut Pasteur, 62]) and in two cell lines (Vero E6 [68] and A549-ACE2 [61]). All pharmacological profiling results can be found in Figure S5 and Table S8. Figure 7 Mapping Regulated Kinases to Kinase Inhibitors Identifies SARS-CoV-2 Therapies (A) Kinase inhibitors (left) mapped to kinases (right) whose activity was regulated by SARS-CoV-2 infection. Lines connecting them indicate known kinase targets for each drug/compound. (B) Vero E6 cells pre-treated with remdesivir at the indicated doses, followed by SARS-CoV-2 infection for 48 h. Percent viral titer compared with mock drug treatment (anti-NP antibody; red line, dots, and text) and cell viability (black) is depicted. Error bars represent SD. (C) As in (B). Vero E6 cells were treated with the CK2 inhibitor silmitasertib. Physical interactions between N protein and the CSNK2A2 and CSNK2B CK2 subunits were observed in a prior study (Gordon et al., 2020). (D) Predicted increased kinase activity for the p38 signaling pathway and drugs/compounds targeting pathway members (ralimetinib, MAPK13-IN-1, and ARRY-797) and upstream drivers (gilteritinib). (E–G) As in (B). Vero E6 cells treated with the AXL inhibitor gilteritinib (E), the MAPK11/14 inhibitor ralimetinib (F), or the MAPK13 inhibitor MAPK13-IN-1 (G) prior to SARS-CoV-2 infection. (H) A549-ACE2 lung epithelial cells were treated with the MAPK14 inhibitor ARRY-797 prior to SARS-CoV-2 infection. (I) Small interfering RNA (siRNA) knockdown of p38 pathway genes in A549-ACE2 leads to a significant decrease in SARS-CoV-2 viral replication (red), as assessed by qRT-PCR in the absence of effects on cell viability (black). ACE2 and non-targeting siRNAs are included as positive and negative controls, respectively. (J and K) Vero E6 or A549-ACE2 cells were treated with PIKFYVE inhibitor apilimod (J) or the CDK inhibitor dinaciclib (K) prior to SARS-CoV-2 infection. See also Figure S6. We found pharmacological inhibitors of CK2, p38 MAPK signaling, PIKFYVE, and CDKs to possess strong antiviral efficacy. Cells were pre-treated with inhibitor molecules, followed by SARS-CoV-2 infection (STAR Methods), and virus quantity (anti-NP antibody against SARS-CoV-2) and cell viability were quantified 48 h after infection. As a positive control and for comparison, remdesivir was tested, and the expected favorable antiviral activity was observed (half maximal inhibitory concentration [IC50 ] = 1.28 μM; Figure 7B). Silmitasertib, an inhibitor of CSNK2A1 and CSNK2A2, was found to possess antiviral activity (IC50 = 2.34 μM; Figures 7C and S5). In conjunction with data supporting physical interaction (Gordon et al., 2020) and co-localization with N protein (Figure 5F), as well as a potential role in remodeling extracellular matrix upon infection (Figures 5 and S3), CK2 signaling appears to be an important pathway hijacked by SARS-CoV-2. Furthermore, silmitasertib is currently being considered for human testing as a potential treatment for COVID-19. To probe SARS-CoV-2 dependence on MAPK signaling, SARS-CoV-2 replication was measured in response to pharmacological and genetic perturbation of MAPK components that were upregulated during infection (Figure 7D). Potent antiviral activity was observed for gilteritinib (Figure 7E; IC50 = 0.807 μM), an inhibitor of AXL kinase, upstream of p38; ralimetinib (Figure 7F; IC50 = 0.873 μM), an inhibitor of MAPK11 (p38ɑ) and MAPK14 (p38β); MAPK13-IN-1 (Figure 7G; IC50 = 4.63 μM), an inhibitor of MAPK13 (p38-δ); and ARRY-797 (Figure 7H; IC50 = 0.913 μM) in A549-ACE2 cells, a MAPK14 inhibitor. To further probe the dependence of SARS-CoV-2 on p38 pathway members, small interfering RNA (siRNA)-mediated knockdown of MAP2K3, p38-δ (MAPK13), and p38-ɣ (MAPK12) was performed in A549-ACE2 cells, and a significant decrease in SARS-CoV-2 replication was observed for all three, with little to no effect on cell viability (Figure 7I). In addition, we noted marked regulation of phosphatidylinositol enzyme activities for PIK3CA, PLCB3, and PIKFYVE, suggesting a potential role for the appropriate balance of phosphatidylinositol species. To target this process, apilimod, a small-molecule inhibitor of PIKFYVE, was tested and found to possess strong antiviral activity in two cell lines (Vero E6, IC50 < 0.08 μM; A549-ACE2, IC50 = 0.007 μM), corroborated by a recent study (Ou et al., 2020; Figure 7J). Lastly, we noted pronounced regulation of CDK signaling pathways (Figure 4B) and cell cycle stage (Figure 6I) during viral infection, suggesting that the virus may regulate the cell cycle to enhance viral replication. Accordingly, strong antiviral activity for the CDK inhibitor dinaciclib was observed across two cell lines (Vero E6, IC50 = 0.127 μM; A549-ACE2, IC50 = 0.032 μM) (Figure 7K). Discussion We used a mass spectrometry-based approach to study perturbations in protein abundance and phosphorylation during SARS-CoV-2 infection. Viral proteins increased, starting 8 h after infection, indicative of viral replication, whereas only small changes in host protein abundance were observed within 24 h. In contrast, large changes were observed in protein phosphorylation, highlighting the degree by which the virus makes use of the host post-translational regulatory systems to promote rapid changes in cellular signaling. Changes in phosphorylation reflect altered activities of kinases that are hijacked during the infection. Based on changes in phosphorylation of their annotated substrates, we estimated changes in activity of 97 of the 518 human kinases. The changes in kinase activity offer insights into the biology of viral infection. Kinases represent ideal drug targets; here, we identified kinases and pathways altered by SARS-CoV-2 infection that can be targeted by 87 FDA-approved drugs and compounds in clinical trials or in preclinical development (Table S8). The most strongly regulated kinases fall into a set of signaling pathways that include p38/MAPK signaling, AKT and ERK signaling, Rho GTPase and CK2 cytoskeleton signaling, and cell cycle regulation. The downregulation of ROCK and PAK kinase activity and upregulation of CK2 cytoskeleton-related targets suggest virus-induced changes in cytoskeleton organization. Imaging of infected cells revealed formation of actin-rich filopodia containing viral proteins. Higher-resolution electron microscopy data confirm the presence of assembled viral particles in these structures. Many viruses, including vaccinia, Ebola, and Marburg, hijack the host cell cytoskeleton to promote egress and rapid cell-to-cell spread across epithelial monolayers. Vaccinia promotes Arp2/3-dependent actin assembly, producing a filopodial protrusion with a virus at the tip (Leite and Way 2015). In contrast, Marburg virus hijacks the unconventional motor protein Myosin X, which promotes filopodium formation and trafficks the virus along the filopodium shaft. The SARS-CoV-2 protein clusters peppered throughout the length of filopodial protrusions more closely resemble Marburg than vaccinia, but additional work is required to understand whether SARS-CoV-2 makes use of either Myosin X motor activity or actin filament assembly to move along filopodia. CK2 is known to phosphorylate myosin proteins at endocytic sites to drive actin polymerization (Fernández-Golbano et al., 2014). Furthermore, CK2 has been found to regulate actin tail formation during vaccinia virus infection, enabling efficient cell-to-cell spread of the virus (Alvarez and Agaisse 2012; Smith and Law 2004). Here the CK2 inhibitor silmitasertib displayed robust antiviral activity, suggesting a role of this kinase in regulating the SARS-CoV-2 life cycle. In addition, kinase activity profiling analysis shows that CDK1/2 activities are significantly reduced by SARS-CoV-2 infection, leading to a S/G2 phase arrest that is similar to infectious bronchitis virus (IBV), a prototypical coronavirus (Dove et al., 2006; Li et al., 2007a), and other RNA viruses (Lilley et al., 2007; Ariumi et al., 2008). Arresting cells in S/G2 phase may provide benefits for viral replication and progeny production by ensuring an abundant supply of nucleotides and other essential host DNA repair/replication proteins (Chaurushiya and Weitzman 2009). The predicted increase in p38/MAPK activity led us to investigate the effects of p38/MAPK inhibition on pro-inflammatory cytokine production and viral replication in SARS-CoV-2-infected cells. Recent immunological studies have indicated that increased IL-6, IL-10, and TNF-α and lymphopenia are associated with severe COVID-19 cases (Pedersen and Ho 2020). The p38/MAPK pathway responds to and controls production of potentially harmful pro-inflammatory cytokines. Several pathogenic viral infections induce a p38/MAPK signaling state that exhibits uncontrolled positive feedback regulation, leading to excessive inflammation associated with severe disease. Inhibition of p38/MAPK signaling suppressed the overproduction of inflammatory cytokines induced by several viral infections, including SARS-CoV, Dengue virus, and influenza A virus, improving survival in mice (Fu et al., 2014; Growcott et al., 2018; Jimenez-Guardeño et al., 2014). However, p38/MAPK inhibition did not directly impair the virus in these cases but, instead, the host’s immune response to the infection. In contrast, during SARS-CoV-2 infection, p38/MAPK inhibition suppressed cytokine production and impaired viral replication by a still unknown mechanism, suggesting that p38/MAPK inhibition may target multiple mechanisms related to COVID-19 pathogenesis. We tested 68 drugs and compounds and found antiviral activity for several that are FDA approved, in clinical testing, or under preclinical development for various diseases, including silmitasertib (CK2, phase 2), gilteritinib (AXL, FDA approved), ARRY-797 (p38, phase 2/3), MAPK13-IN-1 (p38, preclinical), SB203580 (p38, preclinical), ralimetinib (p38, phase 2), apilimod (PIKFYVE, phase 1), and dinaciclib (CDK, phase 3), among others (Figure S5; Table S8). Silmitasertib, a small molecule undergoing clinical trials for various cancers, is now being considered for testing in humans to combat COVID-19. Although the effectiveness of CK2 inhibition may be attributed to its regulation of stress granules (Gordon et al., 2020), viral egress and dissemination could be facilitated by CK2-mediated remodeling of the extracellular matrix (Figure 5). Ralimetinib is currently in phase 2 clinical trials for treatment of ovarian cancer (Patnaik et al., 2016), and ARRY-797 is in phase 3 clinical trials for treatment of cardiomyopathy. The antiviral activity observed for gilteritinib, an FDA-approved drug for treatment of acute myeloid leukemia, is supported by involvement of another AXL inhibitor, bencentinib, in the RECOVERY COVID-19 clinical trial in the United Kingdom. AXL is known to regulate various intracellular signaling pathways (Allen et al., 2002; Hafizi and Dahlbäck 2006), including Ras/ERK, PI3K, and p38 (Allen et al., 2002); AXL inhibition here may contribute to the downregulation of p38 signaling. Apilimod, a PIKFYVE inhibitor, has been described in a recent study to have antiviral capacity (Ou et al., 2020). Here we expand this into a mechanism of regulation by phosphorylation of PIKFYVE upon SARS-CoV-2 infection. Similar to successful antiretroviral therapy for HIV, a combinatorial drug cocktail may be a viable treatment option for SARS-CoV-2 infection. Specifically, combining remdesivir with the kinase inhibitors identified in this study as well as with translation inhibitors and/or modulators of sigma-1 receptor (Gordon et al., 2020) warrants further testing. Furthermore, pairing genetic and pharmacological perturbations in a systematic fashion could identify new combination therapy approaches and illuminate disease mechanisms. The unbiased, global phosphoproteomics approaches used here highlight cellular processes hijacked during SARS-CoV-2 infection. To address the need for improved therapeutic strategies to fight COVID-19, we employed a data-driven approach by mapping phosphorylation profiles of dysregulated signaling pathways to drugs and compounds targeting those signaling pathways. We hope this paradigm can be employed in the future to find additional therapies for COVID-19 and other infectious diseases. Limitations of Study A limitation of the current study is the use of a non-human cell line for proteomics analysis upon SARS-CoV-2 infection; here, an African green monkey cell line (Vero E6) was used because it has been shown previously to be highly permissible to SARS-CoV-2 infection (Harcourt et al., 2020). However, pharmacological inhibition of SARS-CoV-2 was assessed in human lung A549-ACE2 cells in addition to Vero E6 cells. The majority of drug effects were found to be replicated between cell lines. STAR★Methods Key Resources Table REAGENT or RESOURCE SOURCE IDENTIFIER Antibodies Phospho-p38 MAPK (Thr180/Tyr182) Cell Signaling Cat#9211S ɑ-Rabbit IgG HRP Conjugate Biorad Cat#1706515 ɑ-Mouse IgG HRP Conjugate Biorad Cat#1706516 SARS-CoV-2 N (recombinant mouse monoclonal) Donated from the laboratory of Thomas Moran (Mount Sinai) N/A GAPDH Cell Signaling Cat#2118S Phospho-MAPKAPK2 (T334) Cell Signaling Cat#3007T Phospho-CREB/ATF-1 (S133) Cell Signaling Cat#9198S CK2α Abcam Cat#ab70774 SARS-CoV membrane (M) Rockland Cat#100-401-A55 SARS-CoV nucleocapsid (N) Rockland Cat#200-401-A50 AF568-labeled goat-anti-rabbit InvitrogenCat# Cat#A11011 AF647-labeled goat-anti-mouse Invitrogen Cat#A21235 AF488-labeled Phalloidin Hypermol Cat#8813-01 viral NP protein anti-sera produced in the Garcia-Sastre lab N/A Bacterial and Virus Strains SARS-CoV-2, isolate BetaCoV/France/IDF0372/2020 European Virus Archive goes Global (EVAg) IDF0372/2020 SARS-CoV-2, isolate USA-WA1/2020 BEI Resources Cat#NR-52281 SARS-CoV-2, Isolate Muc-IMB-1/2020 Bundeswehr Institute of Microbiology, Munich, Germany N/A Chemicals, Peptides, and Recombinant Proteins cOmplete protease inhibitor cocktail tablets mini, EDTA-free Roche Cat#11846170001 PhosStop phosphatase inhibitor cocktail tablets Roche Cat#4906837001 Sep-Pak C18 cartridge Waters Cat#WAT054955 Sequencing-grade modified trypsin Promega Cat#V5111 Fetal bovine serum (FBS) Thermo Fisher Scientific Cat#10082 DMEM Thermo Fisher Scientific Cat#MT10013CV Water, HPLC grade Sigma-Aldrich Cat#270733-4 L; CAS#7732-18-5 Ni-NTA Superflow beads QIAGEN Cat#30210 Igepal Sigma-Aldrich Cat#I3021; CAS#9002-93-1 Minimal Essential Media (MEM) Corning Cat#10-009-CV Opti-MEM Thermo Fisher Scientific Cat#31985062 Paraformaldehyde, 4% solution in PBS (PFA) Thermo Scientific Cat#MFCD00133991 Formalin Fisher Scientific Cat#SF100-20; CAS#50-00-0, 67-56-1, 7558-79-4, 10049-21-5, 7732-18-5 4′,6-diamidino-2-phenylindole (DAPI) Thermo Scientific Cat#62247 ReprosilPur 1.9 μm particles ESI SOURCE SOLUTIONS Cat#R119 Crystal Violet Solution Sigma Aldrich Cat#HT90132-1L Thermanox coverslips Ted Pella Cat#26028 Silicon Chips Ted Pella Cat#16007 Aluminum specimen mounts Ted Pella Cat#16111 Double-sided carbon tape Ted Pella Cat#16084-1 Karnovsky’s EM fixative Electron Microscopy Sciences Cat#15720 Sodium Cacodylate Sigma Cat#C4945-10G; CAS#6131-99-3 Osmium Tetroxide Electron Microscopy Sciences Cat#19190; CAS#20816-12-0 Potassium Ferrocyanide Sigma Cat#P-3289; CAS#14459-95-1 Uranyl Acetate Ted Pella Cat#19481; CAS#6159-44-0 Spurr’s resin Ted Pella Cat#18300-4221 Iridium target Electron Microscopy Sciences Cat#3431 PolyJet SignaGen Cat#SL100688 SuperScript IV Reverse Transcriptase Invitrogen Cat#18090010 ARRY-371797 (ARRY-797) MedKoo Biosciences Cat#555466; CAS#1036404-17-7 Dilmapimod (SB-681323) MedChemExpress Cat#HY-10404; CAS#444606-18-2 Doramapimod (BIRB-796) SelleckChem Cat#S1574; CAS#285983-48-4 Losmapimod SelleckChem Cat#S7215; CAS#585543-15-3 MAPK13-IN-1 MedChemExpress Cat#HY-18850; CAS#229002-10-2 Neflamapimod (VX-745) SelleckChem Cat#S1458; CAS#209410-46-8 P38 MAPK Inhibitor (4) Cayman Chemical Cat#22219; CAS#1638-41-1 PF-3644022 Tocris Bioscience Cat#4483; CAS#142557-61-7 PH-797804 SelleckChem Cat#S2726; CAS#586379-66-0 Ralimetinib (LY2228820) SelleckChem Cat#S1494; CAS#862507-23-1 SB203580 SelleckChem Cat#S1076; CAS#152121-47-6 TAK-715 SelleckChem Cat#S2928; CAS#303162-79-0 Talmapimod (SCIO-469) MedChemExpress Cat#HY-10406; CAS#309913-83-5 VX-702 SelleckChem Cat#S6005; CAS#745833-23-2 4EGI-1 SelleckChem Cat#S7369; CAS#315706-13-9 A-484954 Tocris Bioscience Cat#4279; CAS#1276121-88-0 Abemaciclib mesylate SelleckChem Cat#S7158; CAS#1231930-82-7 Alisertib SelleckChem Cat#S1133; CAS#1028486-01-2 Apilimod SelleckChem Cat#S6414; CAS#541550-19-0 AT-13148 SelleckChem Cat#S7563; CAS#1056901-62-2 AZD-1208 SelleckChem Cat#S7104; CAS#1204144-28-4 AZD-5363 (Capivasertib) SelleckChem Cat#S8019; CAS#1143532-39-1 AZD-7648 SelleckChem Cat#S8843; CAS#2230820-11-6 Bafetinib SelleckChem Cat#S1369; CAS#859212-16-1 BMS-582949 SelleckChem Cat#S8124; CAS#623152-17-0 Bortezomib (PS-341) SelleckChem Cat#S1013; CAS#179324-69-7 CHIR-98014 SelleckChem Cat#S2745; CAS#252935-94-7 Dabrafenib (GSK2118436) SelleckChem Cat#S2807; CAS#1195765-45-7 Dasatinib SelleckChem Cat#S1021; CAS#302962-49-8 Dinaciclib SelleckChem Cat#S2768; CAS#779353-01-4 Dorsomorphin 2HCl SelleckChem Cat#S7306; CAS#1219168-18-9 EHT 1864 2HCl SelleckChem Cat#S7482; CAS#754240-09-0 Enzastaurin SelleckChem Cat#S1055; CAS#170364-57-5 Fludarabine SelleckChem Cat#S1491; CAS#S1491 Fosbretabulin (Combretastatin A4 Phosphate (CA4P)) Disodium SelleckChem Cat#S7204; CAS#168555-66-6 Ganetespib SelleckChem Cat#S1159; CAS#888216-25-9 Gilteritinib SelleckChem Cat#S7754; CAS#1254053-43-4 Idarubicin HCl SelleckChem Cat#S1228; CAS#57852-57-0 Ipatasertib SelleckChem Cat#S2808; CAS#1001264-89-6 Ispinesib (SB-715992) SelleckChem Cat#S1452; CAS#336113-53-2 Ixazomib (MLN2238) SelleckChem Cat#S2180; CAS#1072833-77-2 K-252a Cayman Chemical Cat#11338; CAS#99533-80-9 KenPaullone SelleckChem Cat#S7917; CAS#142273-20-9 KW-2449 SelleckChem Cat#S2158; CAS#1000669-72-6 LJI308 SelleckChem Cat#S7871; CAS#1627709-94-7 LMK-235 SelleckChem Cat#S7569; CAS#1418033-25-6 LY2584702 Tosylate SelleckChem Cat#S7704; CAS#1082949-68-5 MBQ-167 SelleckChem Cat#S8749; CAS#2097938-73-1 Mdivi-1 SelleckChem Cat#S7162; CAS#338967-87-6 Midostaurin SelleckChem Cat#S8064; CAS#120685-11-2 MK-2206 2HCl SelleckChem Cat#S1078; CAS#1032350-13-2 ML141 SelleckChem Cat#S7686; CAS#71203-35-5 NSC319276 SelleckChem Cat#S7149; CAS#71555-25-4 NVP-BEP800 SelleckChem Cat#S1498; CAS#847559-80-2 Palbociclib (PD-0332991) HCl SelleckChem Cat#S1116; CAS#827022-32-2 PF-03814735 SelleckChem Cat#S2725; CAS#942487-16-3 PF-3758309 SelleckChem Cat#S7094; CAS#898044-15-0 PHA-793887 SelleckChem Cat#S1487; CAS#718630-59-2 Pictilisib SelleckChem Cat#S1065; CAS#957054-30-7 PKR Inhibitor Cayman Chemical Cat#15323; CAS#608512-97-6 QNZ (EVP4593) SelleckChem Cat#S4902; CAS#545380-34-5 R-406 SelleckChem Cat#S2194; CAS#841290-81-1 R-547 SelleckChem Cat#S2688; CAS#741713-40-6 Rapamycin SelleckChem Cat#S1039; CAS#53123-88-9 Ravoxertinib SelleckChem Cat#S7554; CAS#1453848-26-4 Remdesivir BioVision Cat#B2997; CAS#1809249-37-3 Ripasudil (K-115) hydrochloride dihydrate SelleckChem Cat#S7995; CAS#887375-67-9 Romidepsin (FK228, Depsipeptide) SelleckChem Cat#S3020; CAS#128517-07-7 Ruboxistaurin HCl(LY333531) SelleckChem Cat#S7663; CAS#169939-93-9 Sapanisertib SelleckChem Cat#S2811; CAS#1224844-38-5 SAR-407899 Cayman Chemical Cat#21717; CAS#923359-38-0 SB743921 HCl SelleckChem Cat#S2182; CAS#940929-33-9 SGI-1776 free base SelleckChem Cat#S2198; CAS#1025065-69-3 Silmitasertib SelleckChem Cat#S2248; CAS#1009820-21-6 Sorafenib Tosylate SelleckChem Cat#S1040; CAS#475207-59-1 Sotrastaurin SelleckChem Cat#S2791; CAS#425637-18-9 Stattic SelleckChem Cat#S7024; CAS#19983-44-9 Staurosporine SelleckChem Cat#S1421; CAS#62996-74-1 Tanzisertib SelleckChem Cat#S8490; CAS#899805-25-5(freebase) Tasquinimod SelleckChem Cat#S7617; CAS#254964-60-8 THZ1 2HCl SelleckChem Cat#S7549; CAS#1604810-83-4(freebase) Tpl2 Kinase Inhibitor 1 Cayman Chemical Cat#19710; CAS#871307-18-5 U73122 SelleckChem Cat#S8011; CAS#112648-68-7 Volasertib SelleckChem Cat#S2235; CAS#755038-65-4 XL413 (BMS-863233) SelleckChem Cat#S7547; CAS#1169562-71-3 Y-39983 HCl SelleckChem Cat#S7935; CAS#173897-44-4 CCL2/MCP-1, Capture Antibody R&D Systems Cat#MAB679 CCL5/RANTES, Capture Antibody R&D Systems Cat#MAB678 CCL11/Eotaxin, Capture Antibody R&D Systems Cat#MAB320 GM-CSF, Capture Antibody R&D Systems Cat#MAB615 IL-1RA, Capture Antibody Peprotech Cat#500-P209 CCL13/MCP-4, Capture Antibody R&D Systems Cat#MAB327 IL-4, Capture Antibody R&D Systems Cat#MAB604 IL-10, Capture Antibody Peprotech Cat#500-P20 IL-15, Capture Antibody R&D Systems Cat#MAB647 IL-1beta, Capture Antibody R&D Systems Cat#MAB601 CCL8/MCP-2, Capture Antibody Peprotech Cat#500-P35 CXCL9/MIG, Capture Antibody Peprotech Cat#500-P50 CCL3/MIP-1alpha, Capture Antibody R&D Systems Cat#MAB670 CXCL2/GRO-Beta, Capture Antibody Peprotech Cat#500-P104 TGF-Beta1, Capture Antibody R&D Systems Cat#MAB240 CXCL16, Capture Antibody Peprotech Cat#500-P200 CXCL8/IL-8, Capture Antibody R&D Systems Cat#MAB208 MMP-9, Capture Antibody R&D Systems Cat#MAB936 CXCL1/GRO-alpha, Capture Antibody Peprotech Cat#500-P92 IL-7, Capture Antibody R&D Systems Cat#MAB207 CCL4/MIP-1beta, Capture Antibody R&D Systems Cat#MAB271 IL-1alpha, Capture Antibody R&D Systems Cat#MAB200 IL-12 (p70), Capture Antibody Peprotech Cat#500-P154G IL-16, Capture Antibody R&D Systems Cat#MAB316 CCL20/MIP-3a, Capture Antibody R&D Systems Cat#MAB360 TNF-alpha, Capture Antibody R&D Systems Cat#MAB610 Trappin-2/Elafin, Capture Antibody R&D Systems Cat#AF1747 IL-17A, Capture Antibody R&D Systems Cat#MAB317 CXCL10/IP-10, Capture Antibody R&D Systems Cat#MAB266 S100A8, Capture Antibody R&D Systems Cat#DS8900 IL-22, Capture Antibody Peprotech Cat#500-P211 CXCL5/ENA-78, Capture Antibody R&D Systems Cat#MAB654 IL-6, Capture Antibody R&D Systems Cat#MAB206 CCL7/MCP-3, Capture Antibody R&D Systems Cat#MAB282 CCL2/MCP-1, Protein std R&D Systems Cat#279-MC CCL5/RANTES, Protein std R&D Systems Cat#278-RN CCL11/Eotaxin, Protein std Peprotech Cat#300-21 GM-CSF, Protein std R&D Systems Cat#215-GM IL-1RA, Protein std Peprotech Cat#200-01RA CCL13/MCP-4, Protein std R&D Systems Cat#327-P4 IL-4, Protein std R&D Systems Cat#204-IL IL-10, Protein std Peprotech Cat#200-10 IL-15, Protein std R&D Systems Cat#247-ILB IL-1beta, Protein std R&D Systems Cat#201-LB CCL8/MCP-2, Protein std Peprotech Cat#300-15 CXCL9/MIG, Protein std R&D Systems Cat#392-MG CCL3/MIP-1alpha, Protein std R&D Systems Cat#270-LD CXCL2/GRO-Beta, Protein std Peprotech Cat#300-39 TGF-Beta1, Protein std Peprotech Cat#100-21 CXCL16, Protein std Peprotech Cat#300-55 CXCL8/IL-8, Protein std R&D Systems Cat#208-IL MMP-9, Protein std R&D Systems Cat#911-MP CXCL1/GRO-alpha, Protein std Peprotech Cat#300-11 IL-7, Protein std R&D Systems Cat#207-IL CCL4/MIP-1beta, Protein std R&D Systems Cat#271-BME IL-1alpha, Protein std Peprotech Cat#200-01A IL-12(p70), Protein std R&D Systems Cat#219-IL IL-16, Protein std Peprotech Cat#200-16 CCL20/MIP-3a, Protein std Peprotech Cat#300-29A TNF-alpha, Protein std R&D Systems Cat#210-TA Trappin-2/Elafin, Protein std R&D Systems Cat#1747-PI IL-17A, Protein std Peprotech Cat#200-17 CXCL10/IP-10, Protein std R&D Systems Cat#266-IP S100A8, Protein std R&D Systems Cat#DS8900 IL-22, Protein std Peprotech Cat#200-22 CXCL5/ENA-78, Protein std R&D Systems Cat#254-XB IL-6, Protein std R&D Systems Cat#206-IL CCL7/MCP-3, Protein std R&D Systems Cat#282-P3 CCL2/MCP-1, Detection Antibody R&D Systems Cat#BAF279 CCL5/RANTES, Detection Antibody R&D Systems Cat#BAF278 CCL11/Eotaxin, Detection Antibody R&D Systems Cat#BAF320 GM-CSF, Detection Antibody R&D Systems Cat#BAM 215 IL-1RA, Detection Antibody Peprotech Cat#500-P209BT CCL13/MCP-4, Detection Antibody R&D Systems Cat#BAF327 IL-4, Detection Antibody R&D Systems Cat#BAF204 IL-10, Detection Antibody Peprotech Cat#500-P20BT IL-15, Detection Antibody R&D Systems Cat#BAM 247 IL-1beta, Detection Antibody R&D Systems Cat#BAF201 CCL8/MCP-2, Detection Antibody Peprotech Cat#500-P35BT CXCL9/MIG, Detection Antibody Peprotech Cat#500-P50BT CCL3/MIP-1alpha, Detection Antibody R&D Systems Cat#BAF270 CXCL2/GRO-Beta, Detection Antibody Peprotech Cat#500P104BT TGF-Beta1, Detection Antibody R&D Systems Cat#BAM 2462 CXCL16, Detection Antibody Peprotech Cat#500-P200BT CXCL8/IL-8, Detection Antibody R&D Systems Cat#BAF208 MMP-9, Detection Antibody R&D Systems Cat#BAF911 CXCL1/GRO-alpha, Detection Antibody Peprotech Cat#500P92BT IL-7, Detection Antibody R&D Systems Cat#BAF207 CCL4/MIP-1beta, Detection Antibody R&D Systems Cat#BAF271 IL-1alpha, Detection Antibody R&D Systems Cat#BAF200 IL-12(p70), Detection Antibody Peprotech Cat#500-P154GBT IL-16, Detection Antibody R&D Systems Cat#BAF316 CCL20/MIP-3a, Detection Antibody R&D Systems Cat#BAF360 TNF-alpha, Detection Antibody R&D Systems Cat#BAF210 Trappin-2/Elafin, Detection Antibody R&D Systems Cat#BAF1747 IL-17A, Detection Antibody R&D Systems Cat#BAF317 CXCL10/IP-10, Detection Antibody Peprotech Cat#500-P93BT S100A8, Detection Antibody R&D Systems Cat#DS8900 IL-22, Detection Antibody Peprotech Cat#500-P211BT CXCL5/ENA-78, Detection Antibody R&D Systems Cat#BAF254 IL-6, Detection Antibody R&D Systems Cat#BAF206 CCL7/MCP-3, Detection Antibody R&D Systems Cat#BAF282 Critical Commercial Assays MTT assay Roche Cat#11465007001 Luna Universal One-Step RT-qPCR Kit New England Biolabs Cat#E3005L CellTiter-Glo Luminescent Cell Viability Assay Promega Cat#G7570 Direct-zol RNA Miniprep Plus kit Zymo Research Cat#R2072 KAPA SYBR FAST qPCR Master Mix Roche Cat#07959389001 Bio-Rad II Protein Assay Kit Bio-Rad Cat#5000002 Deposited Data Phosphoproteomics data of SARS-CoV-2 infected Vero E6 cells This paper PRIDE Project ID: PXD019113 Supplemental Tables in Mendeley Data This paper https://dx.doi.org/10.17632/dpkbh2g9hy.1 CORUM v3.0 Giurgiu et al., 2019 http://mips.helmholtz-muenchen.de/corum/ RNA-seq dataset for transcription factor activities Blanco-Melo et al., 2020 GSE147507 Phosphoproteomics datasets of biological conditions Ochoa et al., 2016 https://www.embopress.org/doi/full/10.15252/msb.20167295 Functional scores for phosphorylation sites Ochoa et al., 2020 https://idp.nature.com/authorize?response_type=cookie&client_id=grover&redirect_uri=https%3A%2F%2Fwww.nature.com%2Farticles%2Fs41587-019-0344-3 Gene identifier mapping database from BioMart Ensembl http://uswest.ensembl.org/biomart/martview/5ccd994c794265409f21c6194730ce5b Molecular Signatures Database for gene set enrichment analysis Subramanian et al., 2005 https://www.gsea-msigdb.org/gsea/msigdb/index.jsp Kinase-substrate relationships Bachman et al., 2019 https://www.biorxiv.org/content/10.1101/822668v3.full Chlorocebus proteome sequences for peptide searches Uniprot https://www.uniprot.org/proteomes/UP000029965 Sars-CoV-2 genomic sequences (hCoV-19/France/IDF0372/2020) GISAID Accession number EPI_ISL_406596 Experimental Models: Cell Lines ACE2-A549 human lung carcinoma Donated from the laboratory of Brad Rosenberg N/A Caco-2 ATCC Cat#HTB-37 VERO C1008 [Vero 76, clone E6, Vero E6] (Vero E6) ATCC Cat#CRL-1586 HEK293T/17 cells ATCC Cat#CRL-11268 Oligonucleotides ON-TARGETplus Human MAPK12 (6300) siRNA - SMARTpool, 5 nmol Horizon Discovery/Dharmacon Cat#L-003590-00-0005 ON-TARGETplus Human MAPK13 (5603) siRNA - SMARTpool, 5 nmol Horizon Discovery/Dharmacon Cat#L-003591-00-0005 ON-TARGETplus Human MAP2K3 (5606) siRNA - SMARTpool, 5 nmol Horizon Discovery/Dharmacon Cat#L-003509-00-0005 ON-TARGETplus Human ACE2 (59272) siRNA - SMARTpool Horizon Discovery/Dharmacon Cat#L-005755-00-0005 N gene region (reverse) 5′-CGAAGGTGTGACTTCCATG-3′ Eurofins N/A N gene region (forward)5′-TAATCAGACAAGGAACTGATTA-3′ Eurofins N/A Recombinant DNA pLVX-TetOne-Puro empty vector Takara Cat#631847 pLVX-TetOne-Puro-SARS-CoV-2-N-2xStrep This paper N/A pMD2.G, Addgene Gift from Didier Trono Addgene Plasmid #12259 Gag-Pol-Tat-Rev packaging construct Gift from Judd Hultquist pJH045 Software and Algorithms artMS Bioconductor https://www.bioconductor.org/packages/release/bioc/html/artMS.html MSstats Bioconductor https://www.bioconductor.org/packages/release/bioc/html/artMS.html The R Project for Statistical Computing R Core Team, 2019 http://www.r-project.org/index.html Spectronaut Biognosys https://biognosys.com/shop/spectronaut ZEN (Blue edition) Carl Zeiss Microscopy https://www.zeiss.com/zen SPPIDER v2 Porollo and Meller, 2007 http://sppider.cchmc.org/ DoRothEA Garcia-Alonso et al., 2019 https://saezlab.github.io/dorothea/ VIPER Alvarez et al., 2016 https://www.bioconductor.org/packages/release/bioc/html/viper.html Biostrings Bioconductor https://bioconductor.org/packages/release/bioc/html/Biostrings.html clusterProfiler Bioconductor https://bioconductor.org/packages/release/bioc/html/clusterProfiler.html DESeq2 package Love et al., 2014 http://bioconductor.org/packages/release/bioc/html/DESeq2.html Imaris 64x 9.5.1 Oxford Instruments https://imaris.oxinst.com/ MetaMorph (Version 7.8) Molecular Devices https://www.moleculardevices.com/ Other Orbitrap Exploris 480 MS with internal calibration option Thermo Fisher Scientific Cat#BRE725533 ZEISS LSM 800 Confocal Laser Scanning Microscope Carl Zeiss Microscopy Model#ZEISS LSM 800 QuantStudio 6 Real-Time PCR System Applied Biosystems Cat#4485697 FACSymphony BD Life Sciences N/A Infinity 2000 Plate Reader Tecan N/A Luminex MAGPIX platform Luminex Corp N/A LightCycler 480 Instrument II Roche N/A Bal-Tec Drier Balzers, Liechtenstein Cat#CPD 030 Quorum sputter coater Electron Microscopy Sciences, Hatfield, PA Cat#EMS300T D Hitachi field emission scanning electron microscope Hitachi, Tokyo, Japan Model#SU-8000 Leica UC7 ultramicrotome Leica Microsystems N/A FEI BT Tecnai transmission electron microscope Thermofisher/FEI N/A Gatan Rio camera Gatan N/A Celigo imaging cytometer Nexcelcom N/A Resource Availability Lead Contact Further information and requests for resources and reagents should be directed to and will be fulfilled by Nevan J. Krogan (Nevan.Krogan@ucsf.edu). Materials Availability Plasmid pLVX-TetOne-Puro-SARS-CoV-2-N-2xStrep is available upon request from the Lead Contact. Data and Code Availability The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD019113 (Perez-Riverol et al., 2019). An interactive version of phosphorylation data can be found at https://kroganlab.ucsf.edu/network-maps. Supplementary tables have been deposited to Mendeley Data: https://dx.doi.org/10.17632/dpkbh2g9hy.1. Experimental Model and Subject Details Cells African green monkey kidney epithelial Vero E6 (Vero 76, clone E6, Vero E6, ATCC® CRL-1586) authenticated by ATCC and tested negative for mycoplasma contamination prior to commencement were maintained in a humidified atmosphere at 37°C with 5% CO2, in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% (v/v) fetal bovine serum (FBS, Invitrogen). ACE2-expressing A549 cells, a human lung epithelial cell line, were a kind gift from Brad Rosenberg. A549-ACE2 cells were cultured in DMEM supplemented with 10% FBS and maintained at 37°C with 5% CO2. Caco-2 (ATCC® HTB-37) human colon epithelial cells were maintained in a humidified atmosphere at 37°C with 5% CO2 in DMEM-containing 20% (v/v) FBS. Viruses (Institut Pasteur, Paris) The SARS-CoV-2 isolate BetaCoV/France/IDF0372/2020 was supplied by the National Reference Centre for Respiratory Viruses hosted by Institut Pasteur (Paris, France) and headed by Pr. Sylvie van der Werf. The isolate originated from a human sample provided by Dr. X. Lescure and Pr. Y. Yazdanpanah from the Bichat Hospital, Paris, France. The isolate was supplied through the European Virus Archive goes Global (EVAg) platform. Viral stocks were prepared by propagation in Vero E6 cells in DMEM supplemented with 2% FBS. Viral titers were determined by plaque assay in Minimum Essential Media (MEM) supplemented with 2% (v/v) FBS (Invitrogen) and 0.05% agarose. All experiments involving live SARS-CoV-2 were performed in compliance with Institut Pasteur Paris’s guidelines for Biosafety Level 3 (BSL-3) containment procedures in approved laboratories. All experiments were performed in at least three biologically independent samples. Viruses (Mount Sinai, New York) For infection experiments in ACE2-A549 cells, SARS-CoV-2, isolate USA-WA1/2020 (NR-52281), which shares 99.983% sequence identity with the BetaCoV/France/IDF0372/2020 isolate, was deposited by the Center for Disease Control and Prevention and obtained through BEI Resources, NIAID, NIH. SARS-CoV-2 was propagated in Vero E6 cells in DMEM supplemented with 2% FBS, 4.5 g/L D-glucose, 4 mM L-glutamine, 10 mM Non-Essential Amino Acids, 1 mM Sodium Pyruvate and 10 mM HEPES as described previously (Blanco-Melo et al., 2020). All work involving live SARS-CoV-2 was performed in the CDC/USDA-approved BSL-3 facility of the Global Health and Emerging Pathogens Institute at the Icahn School of Medicine at Mount Sinai in accordance with institutional biosafety requirements. Viruses (University of Freiburg, Germany) For infection experiments in Caco-2 cells, SARS-CoV-2 isolate Muc-IMB-1/2020, kindly provided by the Bundeswehr Institute of Microbiology, Munich, Germany, was used. SARS-CoV-2 was propagated in Vero E6 cells in DMEM supplemented with 2% FBS. All work involving live SARS-CoV-2 was performed in the BSL-3 facility of the Institute of Virology, University Hospital Freiburg, and was approved according to the German Act of Genetic Engineering by the local authority (Regierungspraesidium Tuebingen, permit UNI.FRK.05.16/05). Methods Details Vero E6 cell infection for proteomic analysis Vero E6 cells were seeded using 2x106 cells in T25 flasks. The following day cells were either mock infected or infected with SARS-CoV-2 at a MOI of 1 in serum-free DMEM at 37°C for 1 hour. After absorption the 0 hour samples were lysed immediately, while the media for other samples was replaced with 2% FBS / DMEM (Invitrogen) and incubated at 37°C for times indicated before lysis. Cell lysis and digestion Cells were lysed using 1% IGEPAL (Sigma) in PBS (Invitrogen) for 20 minutes at room temperature (RT) to inactivate the virus. These specific lysis conditions were used as this was the approved virus inactivation protocol. Proteins contained in the cell lysate were then immediately precipitated using 90% methanol (v/v) (Sigma) by centrifugation at 20,000x g for 10 minutes. The protein pellets were frozen at −80°C. Precipitated proteins were resuspended in lysis buffer (8 M urea, 100 mM ammonium bicarbonate (ABC), 150 mM NaCl, protease inhibitor (mini-cOmplete, Roche) and phosphatase inhibitors (phosSTOP, Roche). Tris-(2-carboxyethyl)phosphine (TCEP) was added to a final concentration of 4 mM. DNA was sheared via probe sonication, on ice, at 20% amplitude for 20 s., followed by 10 s of rest. This process was performed a total of three times. Following sonication, protein concentration was determined using Bradford assay. Iodoacetamide (IAA) was added to each sample to a final concentration of 10 mM, and samples were incubated in the dark at room temperature (RT) for 30 minutes. Excess IAA was quenched by the addition of dithiothreitol (DTT) to 10 mM, followed by incubation in the dark at RT for 30 minutes. Samples were then diluted with 0.1 M ABC (pH = 8.0) to a final urea concentration of 2 M. Trypsin (Promega) was added at a 1:100 (enzyme:protein w:w) ratio and digested overnight at 37°C with rotation. Following digestion, 10% trifluoroacetic acid (TFA) was added to each sample to a final pH ∼2. Samples were desalted under vacuum using Sep Pak tC18 cartridges (Waters). Each cartridge was activated with 1 mL 80% acetonitrile (ACN)/0.1% TFA, then equilibrated with 3 × 1 mL of 0.1% TFA. Following sample loading, cartridges were washed with 4 × 1 mL of 0.1% TFA, and samples were eluted with 4 × 0.5 mL 50% ACN/0.25% formic acid (FA). 20 μg of each sample was kept for protein abundance measurements, and the remainder was used for phosphopeptide enrichment. Samples were dried by vacuum centrifugation. Phosphopeptide enrichment For each sample batch, 400 μL (30 μL per sample) of 50% Superflow bead slurry (QIAGEN) was added to a 2 mL bio-spin column. Beads were incubated with 4 × 500 μL of 100 mM EDTA for 30 s, washed with 2 × 500 μL H2O, incubated 4 × 500 μL with 15 mM FeCl3 for 1 minute, washed 3 × 500 μL H2O, and washed once with 500 μL of 0.5% FA to remove residual Fe. Beads were resuspended in 600 μL of H2O and 60 μL were aliquoted into a C18 NEST column that had been equilibrated with 150 μL of 80% ACN, 0.1% TFA. 1 mg of digested peptides were resuspended in 75% ACN/0.15% TFA and incubated with beads for 2 minutes, mixed by pipetting and incubated again for 2 minutes. Beads were washed 4 × 200 μL with 80% ACN, 0.1% TFA, then washed 3 × 200 uL with 0.5% FA, incubated 2 × 200 μL with 500 mM potassium phosphate buffer pH 7 and incubated 2 × 200 μL with 0.5% FA for 15 seconds. Phosphopeptides were eluted by centrifugation at 3000 RPM for 30 seconds with 2 × 75 uL of 50% ACN, 0.1% FA. Mass spectrometry data acquisition Digested samples were analyzed on an Orbitrap Exploris 480 mass spectrometry system (Thermo Fisher Scientific) equipped with an Easy nLC 1200 ultra-high pressure liquid chromatography system (Thermo Fisher Scientific) interfaced via a Nanospray Flex nanoelectrospray source. For all analyses, samples were injected on a C18 reverse phase column (25 cm x 75 μm packed with ReprosilPur 1.9 μm particles). Mobile phase A consisted of 0.1% FA, and mobile phase B consisted of 0.1% FA/80% ACN. Peptides were separated by an organic gradient from 5% to 30% mobile phase B over 112 minutes followed by an increase to 58% B over 12 minutes, then held at 90% B for 16 minutes at a flow rate of 350 nL/minute. Analytical columns were equilibrated with 6 μL of mobile phase A. To build a spectral library, one sample from each set of biological replicates was acquired in a data dependent manner. Data dependent analysis (DDA) was performed by acquiring a full scan over a m/z range of 400-1000 in the Orbitrap at 60,000 resolving power (@200 m/z) with a normalized AGC target of 300%, an RF lens setting of 40%, and a maximum ion injection time of 60 ms. Dynamic exclusion was set to 60 seconds, with a 10 ppm exclusion width setting. Peptides with charge states 2-6 were selected for MS/MS interrogation using higher energy collisional dissociation (HCD), with 20 MS/MS scans per cycle. For phosphopeptide enriched samples, MS/MS scans were analyzed in the Orbitrap using isolation width of 1.3 m/z, normalized HCD collision energy of 30%, normalized AGC of 200% at a resolving power of 30,000 with a 54 ms maximum ion injection time. Similar settings were used for data dependent analysis of samples used to determine protein abundance, with an MS/MS resolving power of 15,000 and a 22 ms maximum ion injection time. Data-independent analysis (DIA) was performed on all samples. An MS scan at 60,000 resolving power over a scan range of 390-1010 m/z, a normalized AGC target of 300%, an RF lens setting of 40%, and a maximum injection time of 60 ms was acquired, followed by DIA scans using 8 m/z isolation windows over 400-1000 m/z at a normalized HCD collision energy of 27%. Loop control was set to All. For phosphopeptide enriched samples, data were collected using a resolving power of 30,000 and a maximum ion injection time of 54 ms. Protein abundance samples were collected using a resolving power of 15,000 and a maximum ion injection time of 22 ms. Spectral library generation and raw data processing Raw mass spectrometry data from each DDA dataset were used to build separate libraries for DIA searches using the Pulsar search engine integrated into Spectronaut version 13.12.200217.43655 (Bruderer et al., 2015) by searching against a database of Uniprot Chlorocebus sequences (19,136 proteins, downloaded April 3, 2020) and 29 SARS-CoV-2 protein sequences translated from genomic sequence downloaded from GISAID (accession EPI_ISL_406596, downloaded April 7, 2020) with two mutations (G22661T Spike V367F and G26144T ORF3a G251V) detected by RNASeq analysis of virus stocks. For protein abundance samples, data were searched using the default BGS settings, variable modification of methionine oxidation, static modification of carbamidomethyl cysteine, and filtering to a final 1% false discovery rate (FDR) at the peptide, peptide spectrum match (PSM), and protein level (Elias and Gygi 2007). For phosphopeptide enriched samples, BGS settings were modified to include phosphorylation of S, T, and Y as a variable modification. The generated search libraries were used to search the DIA data. For protein abundance samples, default BGS settings were used, with no data normalization performed. For phosphopeptide enriched samples, the Significant PTM default settings were used, with no data normalization performed, and the DIA-specific PTM site localization score in Spectronaut was applied. Immunofluorescence microscopy Caco-2 cells seeded on glass coverslips were infected with SARS-CoV-2 Isolate Muc-IMB-1/2020, second passage on Vero E6 cells (2x106 PFU/mL) at an MOI of 0.1 or 0.01. At 24 hours post-infection cells were washed with PBS and fixed in 4% paraformaldehyde in PBS for 20 minutes at RT, followed by permeabilization with 0.3% Triton X-100 in PBS for 10 minutes at RT and blocking in 5% fetal calf serum in PBS for 1 hour at RT. Incubation with primary antibodies against CK2α (Abcam, ab70774, 1:500), SARS-CoV membrane (M) protein (Rockland, #100-401-A55, 1:500) and SARS-CoV nucleocapsid (N) protein (Rockland, #200-401-A50, 1:1000) was performed for 1 hour at RT. After washing with PBS, cells were incubated with AF568-labeled goat-anti-rabbit (Invitrogen, #A11011) and AF647-labeled goat-anti-mouse (Invitrogen, #A21235) secondary antibodies (1:200) as well as AF488-labeled Phalloidin (Hypermol, #8813-01, 1:250) for 2 hours at RT. Fluorescence images were generated using a LSM800 confocal laser-scanning microscope (Zeiss) equipped with a 63X, 1.4 NA oil objective and Airyscan detector and the Zen blue software (Zeiss) and processed with Zen blue software and ImageJ/Fiji. For 3D-reconstruction, cells were fixed and stained as indicated and imaged as z stack with 0.15 μm sections. Z stack was processed using Imaris software 64x 9.5.1 and displayed as MIP. To quantify colocalization between casein kinase II and N protein colocalization events in filopodia of SARS-CoV-2 infected Caco-2 cells were counted. 42 ± 19 % of the N protein particles detected in filopodia colocalize with CK2. Length of filopodia was measured using Metamorph (Version 7.8). Distance was measured starting at cortical actin to the tip of Filopodia. Electron microscopy Vero E6 cells were seeded overnight and then infected for 24 hours with SARS-CoV-2 isolate nCoV-WA1-2020 on silicon chips for scanning electron microscopy or Thermanox for transmission electron microscopy. For scanning electron microscopy (SEM), cells were fixed with Karnovsky’s formulation of 2% paraformaldehyde/2.5% glutaraldehyde in 0.1 M Sorenson’s phosphate buffer, and then post-fixed with 1.0% osmium tetroxide/0.8% potassium ferrocyanide in 0.1 M sodium cacodylate buffer washed with 0.1 M sodium cacodylate buffer then stained with 1% tannic acid in dH2O. After additional buffer washes, the samples were further osmicated with 2% osmium tetroxide in 0.1M sodium cacodylate, then washed with dH2O. Specimens were dehydrated with a graded ethanol series from 50%, 75%, 100% x 3 for 5 minutes each, critical point dried under CO2 in a Bal-Tec model CPD 030 Drier (Balzers, Liechtenstein), mounted with double sided carbon tape on aluminum specimen mounts (Ted Pella), and sputter coated with 35 Å of iridium in a Quorum EMS300T D sputter coater (Electron Microscopy Sciences, Hatfield, PA) prior to viewing at 5 kV in a Hitachi SU-8000 field emission scanning electron microscope (Hitachi, Tokyo, Japan). For transmission electron microscopy (TEM), specimens were fixed as described above for scanning electron microscopy and additionally stained overnight with 1% uranyl acetate at 4°C after the second osmium staining and then dehydrated with the same graded ethanol series, and embedded in Spurr’s resin. Thin sections were cut with a Leica UC7 ultramicrotome (Buffalo Grove, IL) prior to viewing at 120 kV on a FEI BT Tecnai transmission electron microscope (Thermofisher/FEI, Hillsboro, OR). Digital images were acquired with a Gatan Rio camera (Gatan, Pleasanton, CA). N overexpression in Vero E6 cells The N-protein cassette was subcloned from pLVX-EF1a-SARS-CoV-2-N-2xStrep-IRES-Puro (Gordon et al., 2020) into pLVX-TetOne-Puro (Takara) using the restriction enzymes EcoRI and BamHI to create pLVX-TetOne-Puro-SARS-CoV-2-N-2xStrep. Sequence integrity was confirmed by Sanger sequencing (GeneWiz). To produce lentiviruses, HEK293T cells (50% confluent in a T175 flask) were transfected with 5 μg of each of the pLVX-TetOne-Puro lentiviral plasmids, 3.33 μg Gag-Pol packaging construct, and 1.66 μg VSV-G envelope (pMD2.G, Addgene) using PolyJet (SignaGen). Culture supernatant was precipitated with a final concentration of 8.5% PEG-6000 and 0.3M sodium chloride (NaCl), incubated at 4°C for 2-4 hours. Virions were pelleted at 3500 rpm for 20 minutes in a spin bucket rotor, suspended in 0.5 mL 1xPBS, and aliquoted for storage at 80°C. Vero E6 cells were seeded in T75 flasks at 50% confluence and transduced with 200 μL precipitated lentivirus derived from pLVX-TetOne-Puro-SARS-CoV-2-N-2xStrep or pLVX-TetOne-Puro empty vector. 48 hours post transduction, 10 μg/mL Puromycin was added to cultures to select transduced cells. Polyclonal stable cell lines were seeded into 15 cm dishes in triplicate and treated with 1 μg/mL doxycycline for 48 hours. Cells were washed in ice cold PBS and harvested by scraping in cold PBS. Cell pellets were lysed directly (8 M urea, 100 mM ammonium bicarbonate (ABC), 150 mM NaCl, protease inhibitor (mini-cOmplete, Roche) and phosphatase inhibitors (phosSTOP, Roche), and protein digestion and phosphopeptide enrichment was performed as described above for SARS-CoV-2 infected Vero E6 cells. During the analysis, one of three replicates of the N-overexpressed samples was found to be an outlier based on principal component analysis (as in Figures S1A and S1B) and was removed. N-overexpression (n = 2) was compared to empty vector (EV) controls (n = 3) during analysis (full data available in Table S1). Kinase activity predictions were also performed the same as for viral infection phosphoproteomics (full data available in Table S4). Cell cycle analysis 1x105 Vero E6 cells were seeded per sample. The following day cells were either mock infected or infected with SARS-CoV-2 (isolate BetaCoV/France/IDF0372/2020) at an MOI of 1. The samples were incubated for 24 hours before detaching cells using Trypsin 0.05% (Thermo). Cells were pelleted at 500x g for 2 minutes, followed by fixation using 4% Formalin (Thermo) for 30 minutes. Cells were washed 2x with PBS (Thermo) before permeabilizing using 0.1% Triton X-100 (Thermo) for 20 minutes. Finally, cells were pelleted again and incubated with 1 mg/mL 4′,6-diamidino-2-phenylindole (DAPI; Thermo) in PBS for 30 minutes before analyzing cell cycle using fluorescence by flow cytometry using the BD FACSymphony with a violet laser (405 nm) and 450/50 nm filter. Gating was performed by first selecting singlets by gating highly correlated cells in an SSC-A versus SSC-H plot. Next, typical cellular morphology was gated using a FSC-A versus SSC-A plot. These singles were then used to draw histograms which were gated for G0/G1, S, and G2/M phases. SARS-CoV-2 infections in ACE2-A549 cells and lysis for cytokine analysis Approximately 5 × 105 A549-ACE2 cells were pre-treated with either DMSO or SB203580 (0.01, 0.1, 1, 10 μM final concentration). After 1 hour of pre-treatment, cells were infected with SARS-CoV-2 (isolate USA-WA1/2020) at an MOI of 0.1 while maintaining inhibitor concentrations in the media for 30 hours post infection. The supernatants from infected cells were analyzed by multiplexed ELISA for the presence of secreted cytokines and chemokines. For RNA analysis, cells were lysed in TRIzol and total RNA was extracted and DNase treated using the Direct-zol RNA Miniprep Plus kit (Zymo Research) according to the manufacturer’s instructions. For protein analysis, whole cell lysates were obtained by lysis in RIPA buffer and analyzed by SDS-PAGE and western blot. Cytokine RT-qPCR analysis Gene expression of selected cytokines as well as SARS-CoV-2 replication was quantified by RT-qPCR. Reverse transcription of extracted total RNA samples was performed using oligo (dT) primers and SuperScript IV Reverse Transcriptase according to the manufacturer’s instructions. Quantitative real-time PCR analysis of cDNA samples was performed using KAPA SYBR FAST qPCR Master Mix (2X) Universal on a LightCycler 480 Instrument II (Roche). For viral quantification, primers specifically targeting the subgenomic viral N RNA were used. ΔCT values were determined relative to the ACTB and ΔΔCT values were normalized to the average of mock infected samples (for host genes) or to infected DMSO treated samples (for viral replication). Error bars indicate the standard deviation of the mean from three independent biological replicates. Multiplexed cytokine ELISA Supernatants from infected cells were evaluated for 34 cytokines/chemokines by multiplex ELISA for the following analytes: CCL2/MCP-1, CCL3/MIP-1α, CCL4/MIP-1β, CCL5/RANTES, CCL7/MCP-3, CCL8/MCP-2, CCL11/Eotaxin, CCL13/MCP-4, CCL20/MIP-3α, CXCL1/GROα, CXCL2/GROβ, CXCL5/ENA-78, CXCL8/IL-8, CXCL9/MIG, CXCL10/IP-10, CXCL16, IL-1α, IL-1β, IL-1RA, IL-4, IL-6, IL-7, IL-10, IL-12p70, IL-15, IL-16, IL-17A, IL-22, GM-CSF, MMP-9, S100A8, TNFα, TGFβ, and Trappin-2/Elafin. All antibodies and cytokine standards were purchased as antibody pairs from R&D Systems (Minneapolis, MN) or Peprotech (Rocky Hill, NJ). Individual magnetic Luminex bead sets (Luminex Corp, CA) were coupled to cytokine-specific capture antibodies according to the manufacturer’s recommendations. Samples were analyzed on a Luminex MAGPIX platform and quantified using a standard curve. For each bead region, > 50 beads were collected per analyte. The median fluorescence intensity of these beads was recorded for each bead and was used for analysis using a custom R script and a 5P regression algorithm. SARS-CoV-2 viral quantification via antibody staining in presence of inhibitors (Mount Sinai, New York). Two thousand (2,000) Vero E6 cells were seeded into 96-well plates and incubated for 24 hours. Two hours before infection, the medium was replaced with 100 μL of DMEM (2% FBS) containing the compound of interest at concentrations 50% greater than those indicated, including a DMSO control. Plates were then transferred into the BSL-3 facility and 100 PFU of SARS-CoV-2 (MOI 0.025) was added in 50 μL of DMEM (2% FBS), bringing the final compound concentration to those indicated. Plates were then incubated for 48 hours at 37°C. After infection, supernatants were removed and cells were fixed with 4% formaldehyde for 24 hours prior to being removed from the BSL-3 facility. The cells were then immunostained for the viral NP protein (anti-sera produced in the Garcia-Sastre lab; 1:10,000) with a DAPI counterstain. Infected cells (488 nM) and total cells (DAPI) were quantified using the Celigo (Nexcelcom) imaging cytometer. Infectivity was measured by the accumulation of viral NP protein in the nucleus of the Vero E6 cells (fluorescence accumulation). Percent infection was quantified as ((Infected cells/Total cells) - Background) ∗100 and the DMSO control was then set to 100% infection for analysis. The IC50 for each experiment was determined using the Prism software (GraphPad). For select inhibitors, infected supernatants were assayed for infectious viral titer using the Median Tissue Culture Infectious Dose TCID50 method. For this, infectious supernatants were collected at 48 hours post infection and frozen at −80°C until later use. Infectious titers were quantified by limiting dilution titration on Vero E6 cells. Briefly, Vero E6 cells were seeded in 96-well plates at 20,000 cells/well. The next day, SARS-CoV-2-containing supernatant was applied at serial 10-fold dilutions ranging from 10−1 to 10−6 and, after 5 days, viral CPE was detected by staining cell monolayers with crystal violet. TCID50/mL was calculated using the method of Reed and Muench. Cytotoxicity cell viability assays (Mount Sinai, New York). Cytotoxicity was also performed using the MTT assay (Roche), according to the manufacturer’s instructions. Cytotoxicity was performed in uninfected Vero E6 cells with same compound dilutions and concurrent with viral replication assay. All assays were performed in biologically independent triplicates. siRNA-based knockdown of host kinases (Institut Pasteur, Paris). Host kinases were knocked-down in A549 cells stably expressing ACE2 cells using OnTargetPlus siRNA SMARTpools (Horizon Discovery). Briefly, A549-ACE2 cells seeded at 2x104 cells/well in 96-well plates 18 hours prior to the experiment. The next day, each well was transfected with 6 pmoles of siRNA SMARTpool, using Lipofectamine RNAiMAX (Thermo Fisher Scientific) according to the manufacturer’s instructions. Twenty-four (24) hours post transfection, the cell culture supernatant was removed and replaced with virus inoculum (MOI of 0.1 PFU/cell). Following a 1 hour adsorption at 37°C, the virus inoculum was removed and replaced with fresh Opti-MEM media (Thermo Fisher Scientific) containing 2% FBS (v/v). Seventy-two (72) hours post-infection (p.i.), the cell culture supernatant was heat-inactivated at 95°C for 5 minutes, and detection of viral genomes performed by RT-qPCR. SARS-CoV-2 viral quantification via RT-qPCR in presence of inhibitors (Institut Pasteur, Paris). Detection of viral genomes was performed by RT-qPCR, directly from the inactivated supernatant (Lista et al., 2020). SARS-CoV-2 specific primers targeting the N gene region: 5′-TAATCAGACAAGGAACTGATTA-3′ (Forward) and 5′-CGAAGGTGTGACTTCCATG-3′ (Reverse) (Chu et al., 2020) were used with the Luna Universal One-Step RT-qPCR Kit (NEB) in an Applied Biosystems QuantStudio 7 thermocycler, with the following cycling conditions: 55°C for 10 minutes, 95°C for 1 minute, and 40 cycles of 95°C for 10 s, followed by 60°C for 1 minute. The number of viral genomes is expressed as PFU equivalents/mL, and was calculated by performing a standard curve with RNA derived from a viral stock with a known viral titer. SARS-CoV-2 titration by plaque assay (Institut Pasteur, Paris). Vero E6 cells were seeded in 24-well plates at a concentration of 7.5x104 cells/well. The following day, serial dilutions were performed in serum-free MEM media. After 1 hour absorption at 37°C, 2x overlay media was added to the inoculum to give a final concentration of 2% (v/v) FBS / MEM media and 0.05% (w/v) Agarose (all Thermo Fisher Scientific) to achieve a semi-solid overlay. Plaque assays were incubated at 37°C for 3 days. Samples were fixed using 4% Formalin (Sigma Aldrich) and plaques were visualized using crystal Violet solution (Sigma Aldrich). siRNA Cell viability assays (Institut Pasteur, Paris). Cell viability in siRNA knocked-down cells was measured using the CellTiter Glo luminescent cell viability assay (Promega) following the manufacturer’s instructions, and luminescence measured in a Tecan Infinity 2000 plate reader. Percent viability was calculated relative to untreated cells (100% viability) and cells lysed with 20% ethanol (0% viability), included in each experiment. Quantification and Statistical Analysis Mass spectrometry data pre-processing Quantitative analysis was performed in the R statistical programming language (version 3.6.1, 2019-07-05). Initial quality control analyses, including inter-run clusterings, correlations, principal components analysis, peptide and protein counts and intensities were completed with the R package artMS (version 1.5.3). Based on obvious outliers in intensities, correlations, and clusterings, 2 runs were discarded from the protein abundance data and 2 runs were discarded from the phosphopeptide data (Figure S1). Additionally, the phosphopeptide data were filtered based on feature (i.e., peptide ion) intensity, removing any single feature with intensity less than 214—this decision was made based on apparent lack of correlation between runs for feature intensities below this intensity. Thus, for both phosphopeptides and protein abundance, we had 2 control time points and 6 infected time points, each with 3 biologically distinct replicates, except for infected at 0 and 2 hours in the phosphopeptide data and control at 0 hours and infected at 0 hours in the protein abundance data which only had 2 replicates each. Statistical analysis of protein abundance changes between control and infected runs were computed using peptide ion fragment data output from Spectronaut and processed using a pipeline of three functions from the R package MSstats (version 3.19.5) (Choi et al., 2014): function MSstats::SpectronauttoMSstatsFormat with default settings other than setting “removeProtein_with1Feature = TRUE”; function MSstats::dataProcess with default settings other than setting “censoredInt = 0,” “featureSubset = highQuality,” “remove_uninformative_feature_outlier = TRUE,” “clusters=7”; and function MSstats::groupComparison with all default settings. Phosphopeptide intensity data were summarized at the peptide ion level along with confident localization of phosphorylation as described in the previous section. Quantification of phosphorylation based on peptide ions were processed using artMS as a wrapper around MSstats, via functions artMS::doSiteConversion and artMS::artmsQuantification with default settings. All peptides containing the same set of phosphorylated sites were grouped and quantified together into phosphorylation site groups. For both phosphopeptide and protein abundance MSstats pipelines, MSstats performs normalization by median equalization, imputation of missing values and median smoothing to combine intensities for multiple peptide ions or fragments into a single intensity for their protein or phosphorylation site group, and statistical tests of differences in intensity between infected and control time points. Identifying significantly regulated proteins When not explicitly indicated, we used defaults for MSstats for adjusted p values, even in cases of N = 2. By default, MSstats uses Student’s t test for p value calculation and Benjamini-Hochberg method of FDR estimation to adjust p values. Phosphorylation fold change data were filtered for quality based on consistency of observations between treatment and controls by requiring the MSstats reported value missingPercentage < 0.60. This value is the proportion of features (i.e., peptide ions) that are missing across the treatment and control replicates for a specific fold change computation. MSstats phosphorylation results had to be further simplified to effects at single sites. The results of artMS/MSstats are fold changes of specific phosphorylation site groups detected within peptides, so one phosphorylation site can have multiple measurements if it occurs in different phosphorylation site groups (see Table S1 for examples). This complex dataset was reduced to a single fold change per site and time point by choosing (per time point) the fold change with the lowest p value, favoring those detected in both treatment and control, i.e., non-infinite log2 fold change values. This single-site dataset, further reduced to those with human orthology (Table S1), was used as the input for kinase activity analysis. Significance filters of adjusted p value < 0.05 and absolute log2 fold change (infected/mock) > 1 were applied to both phosphorylation and abundance datasets. Additionally, for protein abundance data, if a protein was detected in only infected but not mock, or vice versa, we required detection in 3 replicates or by 3 peptide ion fragments in at least 2 replicates. These significance-filtered data were used for enrichment analyses and counts of significant effects (Figure 1). Mapping C. sabaeus and H. sapiens proteins Orthologous pairs of gene identifiers between C.sabaeus and H.sapiens were downloaded from Ensembl using the BioMart web interface (http://uswest.ensembl.org/biomart/martview/6c48d59058381e6b2b198a1f91ba5e50) on April 6, 2020. Ensembl gene identifiers were mapped to UniProt identifiers and protein sequences using a table of ID mappings and reference proteome sequences for H.sapiens downloaded from Uniprot on April 6, 2020. Orthologous pairs of sequences were aligned using the Needleman-Wunsch global alignment algorithm implemented in the R package Biostrings (v. 2.52) function pairwiseAlignment with default parameters. The resulting alignments were used to convert the sequence positions of detected phosphorylations in C.sabaeus to positions in H.sapiens protein sequences. GO enrichment analysis Sets of genes with significant up and down effects were tested for enrichment of Gene Ontology (GO Biological Process, Molecular Function and Cellular Component) terms. Sets of genes were either combined across time points (Figure 1I) or collected separately per time point (Figures S1F and S1G). The over-representation analysis (ORA) was performed using the enricher function of clusterProfiler package (version 3.12.0) in R with default parameters. The gene ontology terms were obtained from the c5 category of Molecular Signature Database (MSigDBv7.1) (Subramanian et al., 2005). Significant GO terms (adjusted p value < 0.01) were identified and further refined to select non-redundant terms. In order to select non-redundant gene sets, we first constructed a GO term tree based on distances (1-Jaccard Similarity Coefficients of shared genes in MSigDB) between the significant terms. The GO term tree was cut at a specific level (h = 0.99) to identify clusters of non-redundant gene sets. For results with multiple significant terms belonging to the same cluster, we selected the most significant term (i.e., lowest adjusted p value). Clustering of phosphorylation changes Phosphorylation sites with one-to-one mapping to human S, T or Y amino acids and showing significant change in phosphorylation (abs(log2FC) > 1 and adjusted p value < 0.05) in one or more conditions were selected and clustered using hierarchical clustering (complete-linkage) with pearson correlation (1-r) as the distance measure. The cluster tree was cut into 5 clusters using a dynamic tree cutting method, cutreeHybrid function in dynamicTreeCut (Langfelder et al., 2008) package in R, with a minimum cluster size of 130 sites. Phosphorylated proteins in each cluster were tested for enrichment of Reactome pathways (Jassal et al., 2020). The over-representation analysis (ORA) was based on the hypergeometric distribution and performed using the enricher function of clusterProfiler package in R (Yu et al., 2012). The pathway terms were obtained from the c2 category (Reactome) of Molecular Signature Database (MSigDBv6.1). Estimation of kinase activities in the time-course experiment Kinase activities were estimated using known kinase-substrate relationships in literature (Bachman et al., 2019). The resource comprises of a comprehensive collection of phosphosite annotations of direct substrates of kinases obtained from six databases, PhosphoSitePlus, SIGNOR, HPRD, NCI-PID, Reactome, and the BEL Large Corpus, and using three text-mining tools, REACH, Sparser, and RLIMS-P. Kinase activities were inferred as a -log10(p value) of Z-test from the comparison of fold changes in phosphosite measurements of the known substrates against the overall distribution of fold changes across the sample. This statistical approach has been previously shown to perform well at estimating kinase activities (Hernandez-Armenta et al., 2017; Casado et al., 2013). We collected substrate annotations for 400 kinases with available data. Kinase activities for kinases with 3 or more measured substrates were considered leaving us with 97 kinases with activity estimates in at least one or more infection time points. Comparison of phosphorylation profiles with other conditions In addition to the time-course experiment, we also estimated kinase activities of cell line perturbation in a range of biological conditions using the K-test based approach, from a large resource of previously published phosphoproteomics datasets comprising of 435 biological conditions which include several drug treatments (chlorzoxazone, paclitaxel), inhibitors (EGFRi, AKTi) and cell cycle states (Ochoa et al., 2016). After selecting for kinase activities derived from 3 or more substrates, a total of 214 kinases had activity estimates for at least one biological condition including all 97 kinases which were estimated in the time-course experiment. Of the 435 conditions, 309 conditions with at least one kinase with significantly changing activity (-log10(p value) from Z-test > 2.5) were selected. Individual time points of the infection experiment were correlated with each selected biological condition based on Pearson correlation of their kinase activity profiles. Conditions with significant correlation (5% FDR) with one or more time points were identified. Several conditions corresponded to the same biological treatment (i.e., drug perturbation) in different concentrations. In such cases, the conditions are represented together with strongest correlation. Phosphorylation of protein complexes Annotations for protein complexes were obtained from CORUM v3.0 database (Giurgiu et al., 2019). Overall significance of changes of phosphosites of a protein complex was inferred from the log2 fold change of individual phosphosites of protein members of the complex using the same Z-test based approach described for estimation of kinase activities. Among the 3512 complexes, 1070 complexes with phosphorylation changes obtained from 4 or more phosphosites were considered. Of these, 67 complexes had significant phosphorylation change (-log10(p value) from Z-test > 2.5) in at least one time point. In order to select non-redundant complexes, we first constructed a complex term tree based on distances (1-Jaccard Similarity Coefficients of shared genes) between the significant complexes. The term tree was cut at a specific level (h = 0.99) to identify clusters of non-redundant complexes. The largest complex was selected from each cluster. Transcription factor activity after SARS-CoV-2 infection To evaluate the effect of SARS-CoV-2 infection at the Transcription Factor (TF) level, we applied DoRothEA (Garcia-Alonso et al., 2019) to RNA-seq datasets of different human lung cell lines from a recent study (GSE147507) (Blanco-Melo et al., 2020). DoRothEA is a comprehensive resource containing a curated collection of TF-target interactions. Each TF-target interaction is associated with a confidence level based on the number of supporting evidence. Here we selected the most reliable interactions (A, B, and C levels) and computed TF activities based on the normalized expression of their targets using the VIPER algorithm (Alvarez et al., 2016). For the TF activity enrichment analysis, VIPER was executed with the Wald statistic resulting from the differential expression analysis at the gene level between controls and SARS-CoV-2 infected cells using the DESeq2 package (Love et al. 2014). In VIPER, we set the eset.filter parameter to FALSE and consider five as the minimum number of targets allowed per regulon. Pharmacological profiling dose response analysis A hill function was fit to each dose response curve using the lsqcurvefit function in MATLAB (R2018a). IC50 (virus) and CC50 (cell viability) values were defined as the concentration at which the percent measure (virus or cell viability quantification) crossed the 50% mark. If the fit curve did not begin above 50% and cross to below 50% throughout the dose response, an IC50 or CC50 value was marked as greater than the maximum tested concentration. For TCID50 and plaque assay results, which are not depicted in a percent scale, the IC50 values were extracted directly from the parameterized hill equation (see Figure S5 and Table S8). Supporting Citations The following references appear in the Supplemental Information: Abraham et al. (2000); Alfonso et al. (2011); Ali et al. (2000); Amin et al. (2014); Aoh et al. (2009); Araujo et al. (1993); Archambault and Carmena (2012); Aubrey et al. (2016); Auld et al. (2005); Bai et al. (1996); Bakircioglu et al. (2011); Bandau et al. (2012); Barkley et al. (2012); Bartrons et al. (2018); Bassermann et al. (2005); Bayliss et al. (2012); Beeson (2013); Beli et al. (2012); Ben Djoudi Ouadda et al. (2018); Beurel et al. (2015); Biaoxue et al. (2016); Bieling et al. (2008); Bischoff et al. (1990); Bitko et al. (2008); Blethrow et al. (2008); Böck et al. (2010); Boeing et al. (2016); Bradshaw et al. (2013); Brattsand et al. (1994); Bromberg and Darnell (2000); Bruinsma et al. (2014); Byon et al. (1998); Cahill et al. (2016a); Cahill et al. (2016b); Cai et al. (2006); Cai et al. (2007); Casado et al. (2014); Chang et al. (2004); Chen et al. (1999); Chen et al. (2002); Chen et al. (2004); Chen et al. (2009); Chen et al. (2013); Chen et al. (2017); Cheng et al. (1998); Cheng et al. (2000); Chiu et al. (2012); Cho et al. (2019); Chou et al. (1990); Chu et al. (2012); Corkery et al. (2015); Costa-Mattioli et al. (2004); Courcelles et al. (2013); Cox et al. (2003); Das et al. (2018); Das et al. (2019); Day et al. (2010); de Haan and Rottier (2005); Decker and Kovarik (2000); DeDiego et al. (2011); DeDiego et al. (2014); Degryse et al. (2018); Deora et al. (2004); Dephoure et al. (2008); Díaz-Moreno et al. (2009); Dou et al. (2015); Eguchi et al. (2009); Elia et al. (2015); Ellis et al. (2019); Eriksson et al. (2004); Escobar et al. (2015); Eto et al. (2012); Eulalio et al. (2007); Fallahi et al. (2016); Fan et al. (1997); Fan et al. (2018); Fang et al. (2000); Farquar et al. (2008); Feldman et al. (1997); Fleming et al. (2000); Folker et al. (2005); Foltz et al. (2009); Formstecher et al. (2001); Franchin et al. (2015); Frémin et al. (2016); Fujimoto et al. (2017); Gagnon and Delpire (2010); Gao and Luo (2006); Garcia-Moreno et al. (2018); Geuens et al. (2017); Gherzi et al. (2004); Goel et al. (2015); Goldstein et al. (2000); Grzenda et al. (2009); Guo et al. (2017); Hajjar (1991); Hajjar et al. (1994); Halder and Johnson (2011); Handschick et al. (2014); Harder et al. (2014); He et al. (2010); He et al. (2018); Heald and McKeon (1990); Heaton et al. (2016); Hedglin and Benkovic (2015); Hein and Nilsson (2016); Hensley and Kursula (2016); Holz et al. (2005); Hou et al. (2011); Hu et al. (2010); Hu et al. (2013); Hu et al. (2015); Huang and Lai (2001); Huang et al. (2006); Huang et al. (2014); Hwang et al. (2014); Ichimura et al. (2005); Inoki et al. (2002); Inoue et al. (2011); Intine et al. (2003); Ishii et al. (2013); Ishikawa et al. (2015); Izawa and Inagaki (2006); Jacque and Stevenson (2006); Jang et al. (2015); Jiang et al. (2014); Johnson et al. (2011); Kane et al. (2002); Kanzawa et al. (2006); Karin (1995); Kastenhuber and Lowe (2017); Katayama et al. (2004); Kennedy et al. (2017); Khurana and Bhattacharyya (2015); Killoran et al. (2015); Kim et al. (2006); Kindrachuk et al. (2015); Kishimoto et al. (2003); Kjolby et al. (2010); Koff et al. (1992); Kohl et al. (2019); Koseki et al. (1999); Kostenko et al. (2009); Kotula et al. (2013); Kovacina et al. (2003); Kraft et al. (2003); Kramer et al. (1998); Kremer et al. (2007); Kuang et al. (2016); Kumar et al. (2014); Kumar et al. (2018); Kunapuli et al. (2006); Labit et al. (2012); Lamarche-Vane and Hall (1998); Larance et al. (2010); Lau et al. (2006); Law et al. (2007); Le Sage et al. (2016); Lee and Myung (2008); Lee et al. (2011); Lee et al. (2012); Lee et al. (2013); Leite and Way (2015); Lekmine et al. (2003); Lekmine et al. (2004); Leukel and Jost (1995); Levy and Lee (2002); Li and Martinez (2011); Li and Nagy (2011); Li et al. (2007a); Li et al. (2007b); Li et al. (2008); Li et al. (2016); Li et al. (2017); Liang et al. (2016); Liao et al. (2005); Lieberman et al. (2017); Lim et al. (2013); Lin et al. (2009); Lin et al. (2013); Lingel et al. (2017); Lisnock et al. (2000); Liu and Schneider (2013); Liu et al. (2011); Liu et al. (2013); Liu et al. (2015); Lobry et al. (2007); Lokireddy et al. (2015); Long et al. (2014); Lu et al. (2014); Lu et al. (2017); Lubas et al. (2011); Luciani et al. (2000); Lukas et al. (2011); Luo et al. (2005); Ly et al. (2014); Ly et al. (2017); Machado-Neto et al. (2014); Machowska et al. (2015); Madureira et al. (2012); Magnuson et al. (2012); Magron et al. (2015); Mahony et al. (2017); Mancini et al. (2017); Mandl et al. (2005); Manning and Toker (2017); Marcos-Villar et al. (2016); Marklund et al. (1993); Marklund et al. (1996); Marquette et al. (2011); Matsuhashi et al. (2019); Matsumoto et al. (2015); Matsuoka et al. (2003); McCamphill et al. (2015); McClain et al. (2002); Mehellou et al. (2018); Meier and Blobel (1990); Meier and Blobel (1992); Meier and Blobel 1992); Mendoza et al. (2011); Menon et al. (2014); Mizutani et al. (2004); Mo et al. (2012); Moldave (1985); Morris et al. (2017); Muda et al. (1996); Mueller et al. (2007); Müller et al. (2014); Na et al. (2015); Nairn et al. (2001); Nandi et al. (2007); Needham et al. (2019); Nielsen et al. (2001); Ogiso et al. (2004); Ohtsubo et al. (1995); Old et al. (2009); Olsen et al. (2006); Olsen et al. (2010); Oshiro et al. (2007); Palamarchuk et al. (2005); Pan (2010); Paquette et al. (2018); Park et al. (2000); Parplys et al. (2015); Patel et al. (2015); Penzo et al. (2019); Pfeiffer et al. (2011); Pires et al. (2017); Ponnam et al. (2019); Popov (2012); Poppe et al. (2017); Potter et al. (2002); Potthoff and Olson (2007); Preisinger et al. (2004); Proteau et al. (2005); Qiao et al. (2016); Qiu et al. (2016); Radic (2016); Radtke et al. (2006); Radulescu and Cleveland (2010); Rahim et al. (2018); Rai et al. (2017); Ramm et al. (2006); Raught et al. (2004); Ravichandran et al. (2001); Rebeaud et al. (2008); Renganathan et al. (2005); Revenko et al. (2010); Riad et al. (2018); Robinson et al. (2018); Robinson-White and Stratakis (2002); Robitaille et al. (2013); Rogalla et al. (1999); Rogers et al. (2015); Roy et al. (2014); Roy et al. (2017); Ryan et al. (2016); Salvi et al. (2009); Sanborn et al. (2011); Sánchez-Martín and Komatsu (2018); Santamaria et al. (2003); Saredi et al. (1997); Schmucker and Sumara (2014); Schweppe et al. (2013); Sengupta et al. (2013); Shah et al. (2018); Shahbazian et al. (2006); Shandilya et al. (2014); Sharma et al. (2011); Sharma et al. (2014); Shi and Sun (2015); Shi et al. (2008); Shih et al. (2012); Sinclair et al. (2006); Smith and Enquist (2002); Snyder et al. (2004); Stoilov et al. (2004); Stokoe et al. (1992); Su et al. (2001); Tacke et al. (1998); Tafforeau et al. (2011); Tan et al. (2004); Tan et al. (2012); Tanenbaum et al. (2006); Tantos et al. (2013); Tarakhovsky and Prinjha (2018); Tavares et al. (2015); Tee et al. (2002); Thaker et al. (2019); Thejer et al. (2020); Tiedje et al. (2015); Tokuyama et al. (2001); Toughiri et al. (2013); Trencia et al. (2003); Trigg and Ferguson (2015); Tsai and Seto (2002); Tsai et al. (1993); Tsoi et al. (2014); Tsvetkov and Stern (2005); Ui et al. (2015); Ujike and Taguchi (2015); Valvezan and Klein (2012); Vind et al. (2011); Vishwanatha and Kumble (1993); Wada et al. (2001); Walter et al. (2002); Walter et al. (2016); Wang and Prives (1995); Wang and Zhang (1999); Wang et al. (2000); Wang et al. (2007); Wang et al. (2008); Wang et al. (2013); Wang et al. (2020); Ward et al. (2004); Wen et al. (1995); Wertz et al. (2015); Willems et al. (1996); Williams et al. (1997); Wohlbold et al. (2012); Wolin and Cedervall (2002); Wortzel and Seger (2011); Wu et al. (2014); Wu et al. (2018); Wurm et al. (2001); Xu et al. (2011); Yang et al. (2012a); Yang et al. (2012b); Yang et al (2014); Yao et al. (2017); Ye et al. (2008); Yim et al. (2013); Yoo et al. (2019); Yu et al. (2009); Yu et al. (2011); Yu et al. (2012); Yui et al. (2001); Zhang et al. (1998); Zhang et al. (2009); Zhang et al. (2015); Zhang et al. (2019a); Zhang et al. (2019b); Zhao et al. (2015); Zhou et al. (2016); Zhu et al. (2001); Zhu et al. (2017); Zou et al. (2007); Zuo and Manley (1993). Supplemental Information Table S1 Phosphorylation and Abundance Data, Related to Figure 1 Contains proteomic data of Vero E6 cells upon SARS-CoV-2 infection. Find the full list of unfiltered phosphorylation sites occurring upon SARS-CoV-2 infection (PhosphoDataFull tab), full list of protein abundance measurements (AbundanceDataFull tab), and filtered list of all detected phosphorylation sites collapsed into single-site measurements (PhosphoDataFiltered tab). Lastly, find the full list of unfiltered phosphorylation sites upon N protein overexpression in Vero E6 cells (PhosphoNOverexpressionFull tab). Column descriptions are indicated in the final tab. Table S2 Viral Phosphorylation Sites, Related to Figure 2 List of phosphosites found on viral proteins from this study as well as from Davidson et al. (2020). Table S3 Enrichments, Related to Figures 1, 4, and S1 Gene ontology enrichments for significantly changed phosphorylation sites (Enrichment.Phosphorylation tab), significantly changed protein abundance (Enrichment.Abundance tab), and phosphorylation dynamics clusters (from Figure 4A) upon SARS-CoV-2 infection (Enrichment.Ph_Clusters tab). Column descriptions are indicated in the final tab. Table S4 Predicted Kinase Activities, Related to Figure 4 Full results of predicted kinase activities for each time point post infection with SARS-CoV-2 (Kinase Act. Viral Infection tab) and N protein overexpression (Kinase Act. N Overexpression tab). Kinase activities are inferred as a -log10(p value) of Z-test from the comparison of fold changes in phosphosite measurements of the known substrates against the overall distribution of fold changes across the sample. Kinase activities possessing any absolute value change greater than 1.5 are indicated. Column descriptions are indicated in the final tab. Table S5 Prioritized Phosphorylation Site Review, Related to Figure 7 and Table S8 Significantly regulated phosphorylation sites upon SARS-CoV-2 infection prioritized as either annotated within the PhosphoSitePlus database or possessing a high functional score (> = 0.75) (Ochoa et al. 2020). Includes literature context for prioritized phosphorylation sites, including their known functions and proposed relevance to SARS-CoV-2 pathogenesis. Phosphorylation sites are partitioned into eight biological contexts. Table S6 Predicted Transcription Factor Activities, Related to Figure 6 Full results of computed transcription factor activities from RNA-seq analysis of SARS-CoV-2 infected human lung cell lines (GSE147507) (Blanco-Melo et al. 2020) using DoRothEA (Garcia-Alonso et al. 2019). Gene symbols, NES scores, and cell lines are depicted. Table S7 Cytokine Profiling Data, Related to Figure 6 Results from Luminex profiling of SARS-CoV-2 infected ACE2-A549 cells pre-treated with p38 inhibitor SB203580. Supernatants from infected cells were evaluated for 34 cytokines/chemokines. Units are pg/mL. Table S8 Drugs and Compounds, Related to Figures 7 and S5 and Tables S4 and S5 Drugs and compounds mapped to top kinase activities (Table S4) and prioritized phosphorylation sites (Table S5). DrugInfo tab depicts known protein targets, approval status, SMILES, supplier, catalog numbers, chembl IDs, annotation of test site and cell line in which tests were performed, IC50 (viral inhibition) and CC50 (cell viability) values for pharmacological profiling. FullDrugResponseData tab depicts mean and standard deviation for drug screening experiments depicted in Figure S5. Acknowledgments This research was funded by grants from the 10.13039/100000002National Institutes of Health (P50AI150476, U19AI135990, U19AI135972, R01AI143292, R01AI120694, P01A1063302, and R01AI122747 to N.J.K.; 1R01CA221969 and 1R01CA244550 to K.M.S.; R01GM133981 to D.L.S.; 1F32CA236347-01 to J.E.M.; U19AI118610 to J.R.J.; and F32CA239333 to M.B.), Defense Advance Research Projects Agency HR0011-19-2-0020 (to N.J.K., A.G.S., and K.M.S.); by the Laboratory for Genomics Research (LGR) Excellence in Research Award (ERA) from the 10.13039/100014220Innovative Genomics Institute at UC Berkeley (grant number 133122P); by CRIP (Center for Research for Influenza Pathogenesis), a NIAID-supported Center of Excellence for Influenza Research and Surveillance (CEIRS; contract HHSN272201400008C) (to A.G.S.); by supplements to 10.13039/100000060NIAID grant U19AI135972 and DoD grant W81XWH-19-PRMRP-FPA (to A.G.S.); and by the generous support of the 10.13039/100007457JPB Foundation, the Open Philanthropy Project (research grant 2020-215611 [5384]), and other philanthropic donations (to A.G.S.); by the Laboratoire d’Excellence “Integrative Biology of Emerging Infectious Diseases” grant ANR-10-LABX-62-IBEID (to M.V.); by the DFG under Germany's Excellence Strategy (EXC-2189, project ID 390939984 to R.G.); by a Starting Grant Award from the 10.13039/501100000781European Research Council (ERC-2014-STG 638884 PhosFunc to P.B.); by the 10.13039/501100002347Federal Ministry of Education and Research (BMBF, Computational Life Sciences grant 031L0181B to J.S.R.); by the Intramural Research Program of the 10.13039/100000002NIH, 10.13039/100000060National Institute of Allergy and Infectious Diseases (to E.R.F. and E.D.W.); and by funding from 10.13039/100007013F. Hoffmann-La Roche and Vir Biotechnology and gifts from The Ron Conway Family. K.M.S. is an investigator of the Howard Hughes Medical Institute. The views, opinions, and findings contained in this study are those of the authors and do not represent the official views, policies, or endorsement of the Department of Defense or the U.S. Government. We would like to acknowledge Desiree Ho, Innovative Genomics Institute, for SARS-CoV-2 virus illustrations. We thank Randy Albert for support with the BSL3 facility and procedures at the Icahn School of Medicine at Mount Sinai, New York. Author Contributions Conceptualization, M.B., D.L.S., and N.J.K.; Infection Experiments, B.M., V.V.R., B.E.N., L.M., C.K., Q.D.T., A.H., T.V., K.M.W., and E.M.; Proteomics Sample Coordination and Preparation, E.S., M.S., J.M.F., J.Z.G., and J.X.; Proteomics Data Acquisition, A.L.R. and D.L.S.; Data Analysis, M.C.M., B.J.P., D.M., C.H.-A., A. Dunham, M. Modak, D.Q., Y.Z., J.R.J., D.L.S., P.B., J.K.L., M.G., M.B., and C.J.P.M.; Figure Generation, M.C.M., B.J.P., D.M., C.H.-A., A. Dunham, M. Modak, D.Q., J.R.J., D.L.S., P.B., and M.B.; Cell Cycle Experiments, R.R. and B.M.; Infection Imaging, S.W., J.K., S.U., G.K., and R.G.; Manuscript Preparation, M.B., D.E.G., K.O., J.R.J., D.L.S., P.B., N.J.K., and R.D.M.; Literature Review, by R.M.K., M. Modak, J.B., A.R., T.P., Q.L., R.H., M.C., M. Muralidharan, M.K., G.J., B.T., J.H., D.L.S., and M.B.; RNA-seq Analysis, A. Dugourd, A.V., and J.S.-R.; Drug Curation, Y.S., J.E.M., K.M.S., A.R.L., E.J.M., E.F., and S.B.; Interactive Map, T.M., M.C.O., Y.C., J.C.J.C., D.J.B, S.K., M.B., and R.M.K.; Electron Microscopy, E.R.F. and E.D.w.; Work Supervision, T.K., J.K.L., A.G., B.S., M.O., J.S.-R., G.K., R.G., B.R.t., K.M.S., J.R.J., D.L.S., A.G.-S., M.V., P.B., and N.J.K. Declaration of Interests The Krogan laboratory has received research support from Vir Biotechnology and F. Hoffmann-La Roche. K.M.S. has consulting agreements for the following companies involving cash and/or stock compensation: Black Diamond Therapeutics, BridGene Biosciences, Denali Therapeutics, Dice Molecules, eFFECTOR Therapeutics (zotatifin and tomivosertib), Erasca, Genentech/Roche, Janssen Pharmaceuticals, Kumquat Biosciences, Kura Oncology, Merck, Mitokinin, Petra Pharma, Qulab Inc. Revolution Medicines (WDB002), Type6 Therapeutics, Venthera, and Wellspring Biosciences (Araxes Pharma). Supplemental Information can be found online at https://doi.org/10.1016/j.cell.2020.06.034.