• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    An Integrated Systems Biology Approach Identifies the Proteasome as A Critical Host Machinery for ZIKV and DENV Replication

    2021-12-03 09:02:22GuangSongEmilyLeeJianboPanMiaoXuHeeSoolRhoYichenChengNadiaWhittShuYangJenniferKouznetsovaCarleenKlumppThomasSamuelMichaelCedricMooreKiJunYoonKimberlyChristianAntonSimeonovWenweiHuangMenghangXiaRuiliHuangMadhuLalNag
    Genomics,Proteomics & Bioinformatics 2021年1期

    Guang Song, Emily M. Lee, Jianbo Pan, Miao Xu, Hee-Sool Rho,Yichen Cheng, Nadia Whitt, Shu Yang, Jennifer Kouznetsova,Carleen Klumpp-Thomas, Samuel G. Michael, Cedric Moore, Ki-Jun Yoon,Kimberly M. Christian, Anton Simeonov, Wenwei Huang, Menghang Xia,Ruili Huang, Madhu Lal-Nag,*, Hengli Tang,*,Wei Zheng,*,Jiang Qian*,Hongjun Song,8,9,0,*,Guo-li Ming,8,9,*,Heng Zhu,*

    1 Department of Pharmacology & Molecular Sciences, Johns Hopkins School of Medicine, Baltimore, MD 21205, USA

    2 Department of Biological Science, Florida State University, Tallahassee, FL 32306, USA

    3 Department of Ophthalmology, Johns Hopkins School of Medicine, Baltimore, MD 21205, USA

    4 National Center for Advancing Translational Sciences, National Institutes of Health, Bethesda, MD 20892, USA

    5 Sir Run Run Shaw Hospital, Zhejiang University School of Medicine, Hangzhou 310016, China

    6 Institute for Cell Engineering, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA

    7 Department of Neuroscience and Mahoney Institute for Neurosciences, Perelman School for Medicine, University of Pennsylvania,Philadelphia, PA 19104, USA

    8 Department of Cell and Developmental Biology,Perelman School for Medicine,University of Pennsylvania,Philadelphia,PA 19104,USA

    9 Institute for Regenerative Medicine, Perelman School for Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA

    10 The Epigenetics Institute, Perelman School for Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA

    KEYWORDS

    Abstract The Zika virus(ZIKV) and dengue virus (DENV) flaviviruses exhibit similar replicative processes but have distinct clinical outcomes.A systematic understanding of virus–host protein–protein interaction networks can reveal cellular pathways critical to viral replication and disease pathogenesis. Here we employed three independent systems biology approaches toward this goal. First,protein array analysis of direct interactions between individual ZIKV/DENV viral proteins and 20,240 human proteins revealed multiple conserved cellular pathways and protein complexes,including proteasome complexes. Second, an RNAi screen of 10,415 druggable genes identified the host proteins required for ZIKV infection and uncovered that proteasome proteins were crucial in this process.Third,high-throughput screening of 6016 bioactive compounds for ZIKV inhibition yielded 134 effective compounds,including six proteasome inhibitors that suppress both ZIKV and DENV replication. Integrative analyses of these orthogonal datasets pinpoint proteasomes as critical host machinery for ZIKV/DENV replication. Our study provides multi-omics datasets for further studies of flavivirus–host interactions, disease pathogenesis, and new drug targets.

    Introduction

    The dengue virus(DENV)and Zika virus(ZIKV)are two closely related pathogens of the Flaviviridae family[1]. Although dengue disease has been recognized in the Americas since the 1600’s, DENV was only isolated in 1943 and is still one of the most widespread global mosquito-borne viruses. DENV contributes to symptoms in 96 million people and results in over 20,000 deaths each year [2,3]. ZIKV was first discovered as a mild, obscure human pathogen in 1947, but has emerged as a major public health concern in the past few years. This is largely due to its role as an etiological agent in several neurological pathologies, including congenital microcephaly and Guillain–Barre syndrome [4].

    The genomes of both DENV and ZIKV are composed of a single positive-strand RNA, which is directly translated into a polyprotein and subsequently processed to generate components necessary for viral replication and assembly [1]. Because of the limited number of proteins encoded by viral genomes,these viruses are obligate intracellular pathogens. This means that they are completely dependent on their hosts for survival and reproduction, which is mediated by direct interactions between the virus and host cellular components[5–7].A better understanding of virus–host interactions can reveal critical cellular pathways that are necessary for viral replication and pathogenesis. In turn, this could be used to identify effective treatment regimens targeting host proteins [5,7–9]. Advancements in high-throughput technologies over the last decade have made it possible to systematically analyze the protein–protein interactome between a virus and its host [10–14]. Previous studies have identified several new host pathways that are essential to the life cycles of several pathogens, including Kaposi’s sarcoma-associated herpesvirus (KSHV) [15–18],influenza virus [19], HIV [20], and Epstein–Barr virus [21].

    Most antiviral drugs are classified as direct-acting antivirals(DAAs).DAAs directly target specific viral proteins critical for infection. While there are many successful DAAs currently in use for viral infections (e.g., hepatitis C virus), it is wellknown that many RNA viruses rapidly develop drug resistance. This is due to the selective stress imparted by targeting essential viral proteins and the high mutation rate in their RNA-based genomes[22,23].For this reason,a drug targeting critical host proteins would provide a higher genetic barrier for a virus to develop drug resistance [6].

    Genetic similarities between DENV and ZIKV, together with recent findings about the host cell dependency factors they share,suggest that these two related flaviviruses likely utilize a similar replicative strategy in the host [24,25]. Consequently, characterization of conserved flavivirus–human protein–protein interactions (PPIs) can reveal critical cellular pathways that are essential for flavivirus infection [5,25,26].On the other hand, differences in PPIs between ZIKV and DENV may provide insight into how these two viruses lead to different pathological outcomes,for example,microcephaly induced by ZIKV[27].Here,we comprehensively surveyed the human proteome with individual ZIKV and DENV proteins to identify virus–host PPI networks. Bioinformatic analyses revealed multiple cellular pathways and protein complexes,including the proteasome complex. In parallel, a RNAi screen targeting druggable genes in combination with a highthroughput chemical genetics approach also revealed overlapping cellular pathways and protein complexes. Through integrative analysis of these three omics datasets, we identified several conserved cellular machineries important for ZIKV and DENV infection, including the proteasome complex.Cell-based assays confirmed that proteasome inhibitors effectively suppressed both ZIKV and DENV replication.Together, our study not only provides a valuable multiomics dataset resource for the field, but also suggests new strategies for understanding the molecular mechanisms of virus–host interactions and pathogenesis and for identifying cellular host-based targets to develop antiviral therapeutics.

    Results

    ZIKV and DENV recombinant proteins

    The flavivirus genome encodes a total of ten proteins, including three structural proteins and seven non-structural proteins.The three structural proteins are the capsid protein (C); the pre-membrane protein (prM), which is subsequently cleaved upon viral maturation into a Pr peptide and a mature membrane protein (M); and the envelope protein (E), which mediates fusion between viral and cellular membranes. The seven non-structural proteins are non-structural protein 1 (NS1),which is required for formation of the replication complex and recruitment of other non-structural proteins to the ERderived membrane structures; NS2A, which is involved in virion assembly and antagonizes the host alpha/beta interferon antiviral response; serine protease NS2B; serine protease NS3; NS4A, which regulates the ATPase activity of NS3;NS4B, which induces the formation of ER-derived membrane vesicles; RNA-directed RNA polymerase NS5; and the short peptide 2K [28]. To construct ZIKV– and DENV–host PPI networks, we cloned the genes from the ZIKV MR766 strain(African strain) and DENV serotype 1 (Figure S1A). We confirmed cloning fidelity by Sanger sequencing (Figure S1B–D;Table S1). Using a previously reported protocol [21], the viral proteins were individually purified from yeast as N-terminal tagged GST fusion proteins and fluorescently labeled (Figure S1E). The quality and quantity of these labeled proteins were evaluated using SDS-PAGE (Figure S1E).

    Considering the various post-translational modifications(e.g.,glycosylation)catalyzed by yeast cells and the importance of correct disulfide bond formation on a protein’s function and binding activity, we decided to focus on the six homologous non-structural proteins (i.e., NS2A, NS2B, NS3, NS4A,NS4B, and NS5) and two variants with the signal peptide 2K (i.e., NS4A + 2K and 2K + NS4B) encoded by ZIKV MR766 strain and DENV serotype 1 to construct virus–host PPI networks.

    Construction of ZIKV– and DENV–human PPIs

    The Human Proteome Microarray v3.0 (HuProtTMarray),comprised of 20,240 immobilized human proteins from>15,000 full-length genes,was used to identify virus–human PPI networks [29]. Each viral protein was fluorescently labeled and individually probed on the HuProt array.Fluorescent signals indicating viral protein bound to immobilized human protein were acquired, normalized, and quantified[30].We used a very stringent cut-off(Z-score ≥15)to identify positive hits for each viral protein(Figure 1A).The assays performed in duplicate showed high reproducibility as measured by Pearson correlation coefficients.An example of binding signals obtained with DENV-NS5 is shown in Figure 1B.

    We identified a total of 1708 ZIKV–host PPIs and 1408 DENV–host PPIs, involving 581 human proteins (Table S2).The majority of host proteins were found to interact with specific individual viral proteins.For example,152 human proteins only interacted with a single viral protein, whereas 75 human proteins bound to two viral proteins (Figure S2A).We found 24 human proteins that interacted with all viral proteins tested,which is possibly a consequence of the common Nterminal GST tag. These 24 human proteins were removed from further analysis.

    We recently used the NS2A PPI dataset to investigate how ZIKV-NS2A causes microcephaly-like phenotypes in the embryonic mouse cortex and human forebrain organoid models [27]. Using a co-immunoprecipitation (Co-IP) assay, we confirmed interactions between ZIKV-NS2A and three endogenous PPI targets (ARPC3, SMAD7, and NUMBL) in neural stem cells [27]. We also evaluated the ability of our approach to recover human proteins known to be targeted by viruses.We acquired a total of 754 human targets from VirusMINT [31] and Virhostome [32]. Of the 557 host proteins identified in our PPI analysis, 54 overlapped with the 754 targets from VirusMINT and Virhostome (hypergeometric P value = 1.9E-5; Figure S2B). The result that only around 10%of the identified hits were shared with the other two databases may be due to the use of different technologies,or due to the fact that different varieties of the virus have specific binding proteins to maintain their replication. Furthermore, we noted that a recently published paper identified 701 and 688 human binding proteins by IP-MS and BioID, respectively(Figure S2C), both of which were based on MS [33]. Of these,48 overlapped with our 581 host proteins (hypergeometric P value = 0.004).

    Host cellular machineries involved in PPIs

    To compare PPIs between ZIKV and DENV, we assembled a global PPI network involving 557 human, 8 ZIKV, and 8 DENV protein nodes(Figure S2D;Table S2).From this data,110 and 42 host proteins were exclusively connected to a single ZIKV or DENV protein, respectively, suggesting that these virus-specific PPIs could contribute to ZIKV- or DENVspecific infection outcomes or pathological effects. For example, PLEK was connected only to ZIKV-NS2A; DDX49 and TTR only to ZIKV-NS4B; and 75 proteins only to ZIKVNS4A. Our recent study also confirmed the interactions of ARPC3 and NUMBL to ZIKV-NS2A but not to DENVNS2A, which was identified in our previous study by using a Co-IP method in HEK293 cells [27].

    In the PPI networks, 368 (66.1%) human proteins were connected to both ZIKV and DENV homologous proteins,supporting the notion that these two related viruses exploit similar cellular machineries. Statistical analysis showed a significant overlap between human proteins recognized among each viral homologous protein pair(Figure 1C).For example,ZIKV-NS3 and DENV-NS3 proteins were found to interact with 238 and 240 human proteins, respectively, of which 187 were shared (77.9%-78.6%; hypergeometric P value = 3.7E-324). Conversely, ZIKV-NS3 and ZIKV-NS4A, two unrelated proteins, interacted with 238 and 401 human proteins,respectively, of which only 168 overlapped (41.9%-70.6%).Similarly, only 127 proteins (52.9%-67.6%) overlapped between 240 DENV-NS3 bound and 188 DENV-NS4A bound human proteins.

    Gene Ontology(GO)analysis for human proteins that were targeted by each individual viral protein revealed several interesting GO features (Figure 2A; Table S3). First, host proteins connected to homologous ZIKV and DENV proteins were often enriched for the same GO terms. This is consistent with the result that a large number of shared host proteins interacted with homologous viral proteins. Second, host proteins targeted by different viral proteins were enriched for diverse biological processes and protein complexes. Third, many different viral proteins interacted with different components of the same enriched biological processes and protein complexes.

    Figure 1 Identification of ZIKV– and DENV–host PPIsA. Sample images of HuProt arrays showing human proteins bound by individual viral proteins. Each human protein was printed in duplicate.The orange,blue,and green boxes represent shared,ZIKV-specific,and DENV-specific interactions,respectively.B.Duplicate experiments performed for each viral protein probe showed high reproducibility. Pearson correlation coefficient analysis showed that signals of duplicate experiments based on DENV-NS5 have a high linear relationship(r=0.961).C.Summary of numbers of unique and conserved virus–host interactions between each ZIKV and DENV homologous pair.

    These observations raised the question of whether the conserved and virus-specific PPIs reflected different biological processes. Indeed, GO analysis of ZIKV/DENV conserved PPIs and ZIKV- or DENV-specific PPIs demonstrated distinct enrichments (Figure 2A). For instance, the GO term of cell–cell adhesion was enriched mainly in human proteins specifically targeted by ZIKV proteins.On the other hand,GO terms of proteasome complex and NIK/NF-kappaB signaling were enriched in PPIs shared by ZIKV and DENV, suggesting that these virus-relevant protein complexes and biological processes may be important for flavivirus infections.For example,six of the non-structural ZIKV proteins (NS2B, NS3, NS4A,NS4A+2K,2K+NS4B,and NS5)interact with eight components in the proteasome complex (Figure 2B). Similar phenomena were observed for the spliceosome complex(Figure 2C). Furthermore, Co-IP performed in HEK293FT cells verified the physical binding of proteasome subunits PSMA3/PSMB4 to ZIKV-NS3 and PSMA1/PSMA3/PSMB4 to ZIKV-NS5 (Figure 2D). Consistent with our finding, a recent study reported that DENV-NS5 protein interfered with host mRNA splicing through direct binding to proteins in the spliceosome complex [34].

    RNAi screening identified critical host proteins

    To validate whether host proteins enriched in the aforementioned biological processes and protein complexes were functionally involved in ZIKV infection, we carried out a siRNA-knockdown assay similar to those used for other viruses [12,35–38]. Specifically, 10,415 druggable target genes were individually knocked down and ZIKV NS1 protein level was measured using a high-throughput homogenous timeresolved fluorescence (HTRF) assay as a surrogate for viral load in ZIKV-infected HEK293 cells.Among the 10,415 target genes,knockdown of 120(1.2%)genes resulted in significantly reduced NS1 levels (> 30%; Table S4). GO analysis revealed that proteasome,spliceosome,RNA polymerase,COPI vesicle coat, and Eukaryotic 43S preinitiation complex were significantly enriched among those 120 genes, with proteasome showing the lowest P value (P = 3.8E-25; Figure 3A).

    Of the 10,415 target genes, protein products of 327 genes were found to interact with ZIKV proteins during our PPI analysis, and individual siRNA-knockdown of three of them(i.e., PSMC3, PSMA1, and OVOL2) resulted in > 30%reduction of NS1 levels. Notably, a significant increase in the success rate of the knockdown assays was observed for those genes whose protein products were found in the enriched GO terms identified by the PPI analysis (Figure 3B). For example, individual knockdown of 20 of the 47 members(42.5%)in the proteasome complex showed>30%reduction of NS1 levels (Figure 3B).

    High-throughput drug screening identified small molecule inhibitors

    Figure 2 GO analyses of host proteins in the PPI networksA. Enriched GO terms in the categories of molecular function, biological process, and cellular component are found in both shared and virus-specific PPI networks. The folds of enrichment are color-coded by P value. As examples, interactions of six non-structural ZIKV proteins(NS2B,NS3,NS4A,NS4A+2K,2K+NS4B,and NS5)with proteasome complex(B)and interactions of eight non-structural ZIKV proteins (NS2A, NS2B, NS3, NS4A, NS4A + 2K, 2K + NS4B, SN4B, and NS5) with spliceosome complex (C) were shown respectively. Here, only the subunits capable of binding with ZIKV proteins were included. Circles with bright blue outlines indicate previously reported virus binding proteins. D. Co-IP of overexpressed FLAG-tagged ZIKV proteins (NS3 and NS5) and V5-tagged human proteasome subunits (PSMA1, PSMA3, and PSMB4) in HEK293FT cells. IP assays were performed with anti-FLAG mAb magnetic beads and eluted fractions were analyzed by Western blot using mouse anti-V5 antibodies.Mouse IgG magnetic beads were used as a negative control to evaluate any non-specific binding on the beads. Inputs correspond to 2% of total lysate incubating with anti-FLAG mAb magnetic beads.

    To further substantiate our results, we employed an independent chemical genetics approach to screen for and validate chemical compounds that target host proteins essential for viral replication and that interfere with the ZIKV life cycle.A total of 6016 compounds, including the Library of Pharmacologically Active Compounds (LOPAC, 1280 compounds),the NIH Chemical Genomics Center (NCGC) pharmaceutical collection of 2816 approved drugs, and 1920 bioactive compounds [39], were screened for antiviral activity against ZIKV infection of HEK293 cells. ZIKV infection was quantified by ZIKV-NS1 antibody-based Time-Resolved Fluorescence Resonance Energy Transfer (ZIKV-NS1 TR-FRET) detection[40]. The ZIKV-NS1 TR-FRET assay measures the total amount of intra- and extracellular NS1 protein levels in infected cell culture,which was used as an indicator for ZIKV replication levels in cells(Figure 4A).Of the 6016 compounds,256 were identified as preliminary hits and selected for secondary validation by the ZIKV-NS1 TR-FRET assay and cytotoxicity evaluation in the same cells(Figure 4A,Figure S3).Viral inhibition was confirmed for 217 of the preliminary hits,and 134 compounds exhibited greater than four-fold selectivity of ZIKV-NS1 inhibition over compound cytotoxicity(Table S5), which included the 24 compounds previously reported [41].

    Figure 3 Critical host proteins for ZIKV replicationA. STRING analysis of genes that significantly affected ZIKV replication in RNAi screening. B. Percentage of genes with over 30%reduction of NS1 levels by siRNA-knockdown among all genes in a specific category. The‘‘All” group indicates the collection of 10,415 siRNA-targeted druggable genes;the‘‘ZIKV-host PPI”group indicates the 327 genes whose protein products were found to interact with ZIKV proteins in the PPI dataset. Note the high success rate (20 out of 47 members) in the category of ‘‘Proteasome complex.”

    Based on the reported mechanisms of action(https://tripod.nih.gov/npc/),ZIKV inhibition exhibited by 92 of 134 effective compounds was mainly mediated by 12 biological categories[39]: proteasome inhibitors, antibacterials/antifungals, CDK inhibitors, PI3K inhibitors, Chk1 inhibitors, antiprotozoal agents, DNA topoisomerase-I inhibitors, kinase inhibitors,inosine 5′-monophosphate dehydrogenase inhibitors, mTOR inhibitors, tubulin inhibitors, and anticancer agents (Figure 4B).Then,a ZIKV virus titer assay was performed to further confirm the anti-ZIKV activity of these compounds.Among these compounds, the antiviral activity of emetine was confirmed in the mouse models of ZIKV infection [42],which validated our compound screening approach.

    Integrative analysis of the omics data

    We compared data deposited in Drugbank [43], Therapeutic Target Database(TTD)[44],and STITCH 5.0[45]and identified 1065 human proteins as targets of the 134 effective anti-ZIKV compounds from our screen. Of the 1065 protein targets, 45 were found to interact with ZIKV proteins in our PPI analysis (Figure 5A). STRING analysis revealed that the majority of these proteins (80.0%, 36/45) were highly connected via functional associations, such as physical interactions, co-expression, tissue specificity, and functional similarity [46]. Indeed, 46 connections were found among 45 proteins, compared to only 18 expected connections (PPI enrichment P value =1.3E-8). GO analysis of these proteins revealed significant enrichment in proteasome, vesicle, and in the regulation of cell death (Figure 5A; Table S6).

    Among the 6016 tested compounds, 3671 have known targets. Of these 3671 compounds, 98 (2.67%) showed selective inhibition against ZIKV infection. For the 766 drugs that are known to target proteins in our PPI analysis,29(3.79%)were effective, demonstrating a 1.42-fold enrichment. Individual pathways and complexes also showed enrichment for identified effective drugs, except for peroxisome and oxidation–reduction process (Figure 5B).

    Proteasome inhibitors suppress ZIKV and DENV replication

    The integration of the three orthogonal datasets presented strong evidence that the same conserved cellular machineries play an important role in ZIKV and DENV replication. The proteasome complex stood out for several reasons. First, the PPI network analysis revealed that six ZIKV and six DENV proteins interacted with eight and seven proteasome subunits,respectively, most of which are part of the 20S core particle(Figure 6A and B). Second, individual knockdown of 20 proteasome genes resulted in substantially reduced ZIKV replication in the RNAi screen (Figure 3B). Third, the proteasome complex was the most significantly enriched pathway targeted by the 134 effective compounds identified by the chemical genetics approach to inhibit ZIKV.

    To further validate our results, we selected six proteasome inhibitors (MLN-2238, carfilzomib, bortezomib, delanzomib,oprozomib, and MG-115) for further evaluation of their inhibitory activities on ZIKV and DENV in the human glioblastoma cell line SNB-19. We used a recent clinical isolate of the Puerto Rico PRVABC59 ZIKV strain for this analysis.The cultures were infected with ZIKV or DENV at a multiplicity of infection(MOI)of 1 in the presence of these compounds at a concentration of 1 μM.DMSO served as the negative control. In this assay, all of the proteasome inhibitors tested suppressed both of ZIKV and DENV envelope expression as compared to the DMSO control (Figure 6C, Figure S4A and B).

    Figure 4 Small molecule inhibitors against ZIKV replicationA. Flowchart of compound screening and confirmation with the ZIKV-NS1 TR-FRET assay. Precultured cells in 1536-well plates were treated with 6016 compounds for 1 h, and then infected with virus for 1 day, followed by the ZIKV-NS1 TR-FRET assay. Of the 6016 compounds, 256 were identified as preliminary hits and selected for secondary validation by the ZIKV-NS1 TR-FRET assay and cytotoxicity evaluation with the same cells.217 of the preliminary hits were confirmed and 134 compounds exhibited greater than four-fold selectivity of ZIKV-NS1 inhibition over compound cytotoxicity. B. Summary of behaviors and IC50 values of 12 groups of potent compounds categorized based upon their reported mechanisms of action.Values represent mean±SD(n=3 cultures).Curves represent best fits for calculating IC50.

    Figure 5 Integrative analysis of PPI and chemical genetics screenA.STRING analysis of 45 anti-ZIKV drug target human proteins that were found to interact with ZIKV proteins in our PPI analysis.B.Functional association networks among the proteins that interact with viral proteins and are targeted by effective compounds.

    Finally, we used a colorimetric focus-forming unit (FFU)assay to determine the dose response and IC50of these compounds on ZIKV production.Consistent with the intracellular antigen expression assay,all six proteasome inhibitors reduced infectious ZIKV production, with IC50values for MLN-2238,carfilzomib, bortezomib, delanzomib, and oprozomib in the nanomolar range (Figure 6D and E).

    Discussion

    Figure 6 Experimental validation of the proteasome inhibitorsA. PPI network analysis of virus proteins and human proteasome subunits reveals that most of the interacting proteasome subunits are part of the 20S core particle. B. Percentage of the ZIKV-binding subunits in 26S proteasome and its two sub-complexes, the 20S core particle and the 19S regulatory particle. C. Inhibition of ZIKV expression in human glioblastoma cell line SNB-19 by a panel of proteasome inhibitors. The SNB-19 cells were infected by ZIKV PRVABC59 (MOI = 1) in the presence of 1 μM of each inhibitor and then incubated for 48 h before the cultures were analyzed for ZIKV-E protein expression by immunostaining.Mock indicates cells without ZIKV infection. Scale bar: 100 μm. D. and E. Sample images (D) and quantification (E) of titer assay to assess the potency of the proteasome inhibitors against infectious ZIKV production in SNB-19 cells.All data were normalized to that of 0 μM for each compound.Dose-dependent antiviral activity is presented as fluorescent focus-forming units per ml(FFU/ml)and data are represented as mean±SD(n=6).Curves represent best fits for calculating IC50 values(listed to the right).MOI,multiplicity of infection;ZIKV-E,ZIKV envelope.

    In this study,we employed three high-throughput platforms to investigate host cellular machineries that are critical for ZIKV and DENV replication. First, HuProt arrays were used to screen for direct PPIs between each ZIKV/DENV protein and 20,240 human proteins. Next, an RNAi screen targeting 10,415 druggable genes was adapted to identify the critical human genes required for ZIKV replication.Lastly,a chemical genetics approach was employed to screen 6016 bioactive compounds for their ability to inhibit ZIKV replication. We have confirmed the anti-ZIKV activities of 217 compounds, with 134 of them having a selectivity index greater than 4-fold,which represents a comprehensive list of approved drugs and bioactive compounds with anti-ZIKV activity. Integration of the three independent omics datasets identified several host machineries, including the proteasome complex, the spliceosome complex, and regulation of mRNA stability. The integrated data, including PPIs, RNAi screening, and compound screening focused on ZIKV and DENV, provide useful resources for further studies to understand viral biology and disease pathogenesis and to identify new drug targets. Moreover, the systematic screening illustrated by our approach can be readily implemented to study other virus–host interactions in order to uncover the nuances of disease pathogenesis and discover novel therapeutic strategies.

    Our multi-omics datasets could have many applications.As an example, we recently took advantage of the PPI dataset to understand molecular mechanisms underlying the differential pathogenic impact on host cells induced by ZIKV and DENV[27]. Consistent with the clinical phenotype that ZIKV infection, but not DENV infection, could lead to microcephaly,our functional screen showed that expression of ZIKVNS2A, but not DENV-NS2A, leads to reduced proliferation and accelerated depletion of cortical neural stem cells in both embryonic mouse cortex in vivo and cultured human forebrain organoids.To understand how these two very similar proteins lead to different consequences in the same host cells,we mined the PPI dataset (Table S2) and found differential interactions of ZIKV-NS2A and DENV-NS2A with adherens junction proteins. We further validated this finding in neural stem cells with endogenous proteins[27].This critical information generated the hypothesis that the differential impact of ZIKV-NS2A and DENV-NS2A on adherens junctions may underlie their differential impact on neural stem cell properties. We tested and confirmed this hypothesis in both in vivo embryonic mouse cortex and in vitro human brain organoid models [27]. Other viral proteins have also been implicated in the pathogenesis of virus infection; for example, ZIKV-NS4A and ZIKVNS4B cooperatively suppressed the Akt-mTOR pathway to inhibit neurogenesis and induce autophagy in human fetal neural stem cell [47]. Additionally, we found that targeting multiple components of the same protein complexes/signaling pathways seems to be a reoccurring event in pathogen–host interactions. Using the spliceosome complex as an example(Figure 2C), host proteins API5 and HNRNPDL were found to interact with ZIKV proteins NS4A and NS5, respectively.It is an intriguing finding that the same process/complex can be targeted by a pathogen at different times. It is conceivable that such ‘‘multivalency” interactions could serve as an effective means to ensure the robust hijacking of the host cell machinery by a pathogen.For example,in one of our previous studies using in vitro phosphorylation assays on human protein arrays, we observed that four conserved viral protein kinases encoded by four different herpesviruses could all phosphorylate 14 components of the DNA damage response pathways,such as TIP60, RAD51, RPA1, and RPA2 [17]. In-depth in vivo studies confirmed that these phosphorylation events played an important role in promoting viral DNA replication in all four viruses. In another study, we observed that a secreted protein kinase ROP18, encoded by Toxoplasma gondii, could phosphorylate multiple components in the MAPK pathway [29]. A third example is the observation that the KSHV-encoded LANA protein could bind to all three components of the NER damage recognition/verification complex XPA–RPA (i.e., XPA, RPA1, and RPA2) [48]. Therefore,our virus–host PPI database can be used to explore both conserved and unique pathogenic processes induced by ZIKV and DENV in different cellular contexts in the future.

    In this study, we focused the investigation of our datasets on viral replication to identify critical cellular machineries as candidate drug targets [19]. Using high-throughput drug screening to reveal hijacked host machinery, we identified potential antiviral compounds with a higher genetic barrier for the virus to develop drug-resistance. In addition, we could potentially use these host-targeting drugs as broad-acting antivirals for closely related viruses,such as DENV and ZIKV,because of their substantially overlapping PPI networks with the host. Integrative analysis of independently identified pathways and PPI networks presents a strong case for the proteasome as a conserved,critical machinery for ZIKA and DENV replication.The proteasome complex is a part of the ubiquitin–proteasome pathway and regulates many fundamental cellular processes[49].Emerging evidence implicates the proteasome as a critical player in viral pathogenesis by modulating the function of viral proteins to favor viral propagation and evade the host immune response[50–52].Until now,there have been few FDA approved antiviral drugs targeting intracellular host proteins,due to the potential side effects[9].Notably,Maraviroc,a CCR5 receptor antagonist, has been approved as an antiretroviral drug for the treatment of HIV infection, which could prevent viral entry by blocking binding of viral envelope gp120 to CCR5[53].Several proteasome inhibitor drugs tested in this study,including carfilzomib and bortezomib,have been approved by the FDA for the therapy of various cancers,such as breast cancer,multiple myeloma,and Hodgkin’s lymphoma[54–57]. Consequently, these drugs could potentially be repurposed to further evaluate their efficacy and tolerance in a clinical setting as novel therapies for ZIKV and DENV infection.

    In summary, we discovered a multitude of cellular pathways and protein complexes related to ZIKV and DENV infection by integrating three high-throughput systems biology methods, i.e., ZIKA/DENV–human PPIs, a druggable gene screen, and a high-throughput chemical genetics screen. We identified the human proteasome as a conserved, critical machinery for ZIKV and DENV replication with functional confirmation by pharmacological proteasome inhibitors. We also found a comprehensive list of 134 selective ZIKV inhibitors that span over 12 cellular pathways and mechanisms.Our study provides a rich resource of multi-omics datasets for future investigation of viral pathogenesis and drug development and highlights a systematic biological approach to investigate virus–host interactions.

    Materials and methods

    Viral cDNA preparation

    The African prototype ZIKV strain MR766 and DENV serotype 1 Hawaii strain were used to infect mosquito cells,as previously described [58]. Lysates of virus-infected mosquito cells were prepared,and 1 μg of the total RNA was used to prepare cDNA by Superscript III (Catalog No. 18080044, Thermo-Fisher Scientific, Waltham, MA) for PCR templates.

    Gateway cloning and protein expression

    Gateway cloning and protein expression were performed using the method as in our previous publication[59].In short,primer sets with the attB1 or attB2 sequences at the 5′- and 3′-ends respectively (Table S1) were designed to amplify the fulllength viral genes,which were then cloned into Gateway Entry vector pDONR221 using the Gateway recombination reaction(Catalog No. 11789021, ThermoFisher Scientific). Each initial cloning was examined by BsrGI (Catalog No. R0575S, New England Biolabs,Ipswich,MA)digestion and Sanger sequencing. Then, each insert viral gene was shuttled into the yeast expression vector pEGH-A to carry out the protein expression.

    Protein labeling

    The quality of each ZIKV and DENV protein was determined using SDS-PAGE, followed by Coomassie staining. Proteins that passed this quality control test were then labeled directly with NHS-tethered Cy5 dye (Catalog No. GEPA15101,Sigma-Aldrich, St. Louis, MO) on glutathione beads. After quenching the dye molecules, the labeled protein was eluted and the quantity of these purified proteins was examined on SDS-PAGE gels.

    Identification of virus-binding host proteins on HuProt arrays

    PPI assays on the Huprot array and signal extraction of each spot were performed using the same methods described previously[27].In short,the signal intensity(Rij)of a given protein spot (i,j) was generated as foreground signal (Fij) minus the corresponding background signal (Bij). The averaged Rijfrom duplicate spots was defined as the signal intensity of the protein probe (Rp). For the replicate samples, the signal profiles were quantile normalized to a merged profile. Using a similar method as described in our previous study[60], the Z-score of each binding assay with a virus protein was computed based on the distribution of Rp.SD andrepresent the standard deviation and mean of the noise distribution on the array,respectively.A stringent cutoff(Z ≥15) was used to determine the positive hits in this study.The proteins determined as positives in all assays were removed from further analysis.

    Comparison to other datasets

    The statistical significance of the overlap between our set of identified virus-binding human proteins and those deposited in VirusMINT and Virhostome was calculated using a hypergeometric test implemented in R[31,32].The number of background proteins was defined as the number of unique wellannotated human proteins detected in our HuProt Array(n = 13,816).

    Functional annotation

    Database for Annotation, Visualization, and Integrated Discovery (DAVID) was used to identify the enriched functional terms(molecular function,cellular component,biological process, and KEGG pathway) for virus-binding proteins [61].Some enriched terms (P value < 0.05) were selected and represented in a heat map by the fold change.

    PPI network

    Virus–human PPIs identified in this study were input into Cytoscape to construct flavivirus–host PPI networks [62].Human PPIs were extracted and drawn from STRING 10.0[46].The significance of functional terms and interaction numbers were also calculated and provided by STRING.

    Drug–target interaction

    Drug targets were collected from three resources, Drugbank,TTD, and STITCH 5.0 [43–45]. Drugbank and TTD include known targets of experimental drugs and FDA-approved drugs. STITCH combines chemical–protein interactions from experimental chemical screens, prediction, known databases,and text mining.For chemical–protein interactions in STITCH 5.0, only those with greater than 0.7 of high combined confidence score and with experimental or database scores were chosen for analysis.Those targets not identified as positive hits on HuProt arrays were removed from this study.

    Propagation of ZIKV

    ZIKV stocks were generated in Aedes albopictus clone C6/36 cells as previously described [40]. The ZIKV PRVABC59 strain was purchased from American Type Culture Collection(ATCC; Manassas, VA). The ZIKV MR766 stock was purchased from Zeptomatrix, Buffalo, NY. Briefly, a T-75 flask of C6/36 cells (90%–95% confluency) was inoculated with 1 × 106ZIKV virions in low volume (3 ml) for 1 h, rocking it every 15 min. After 1 h, 17 ml of media was added and C6/36 cells were maintained at 28°C in 5%CO2.At 7–8 days post-viral inoculation, supernatants were harvested, filtered,and stored at -80 °C. ZIKV titer was determined by an FFU assay.

    Viral infection

    For SNB-19 and HEK293 cell infections,cells were seeded into 12-or 96-well plates 1 day prior to viral infection.For SNB-19 cells, compounds were added 1 h before addition of ZIKV at MOI = 1. SNB-19 cells were harvested at 48 h after infection for analysis by immunofluorescence. For viral production assays, infected SNB-19 cell supernatant was harvested 24 h post infection. ZIKV and DENV titers in cell supernatants were measured by an FFU assay (FFU/ml), as previously described [40].

    Co-IP and Western blot

    FLAG-tagged ZIKV proteins (NS3 and NS5) and V5-tagged human proteasome subunits (PSMA1, PSMA3, and PSMB4)were overexpressed in HEK293FT cells. For Co-IP assays,HEK293FT cells were incubated in the lysis buffer (Catalog No. 9803, Cell Signaling Technology, Danvers, MA) for 30 min on ice. After sonication and centrifugation, the supernatants were subjected to IP with anti-FLAG mAb magnetic beads(Catalog No.M8823,Sigma-Aldrich)at 4°C overnight.Then, the beads were washed six times using lysis buffer and used to perform an immunoblot assay with mouse anti-V5 antibodies (Catalog No. R960-25, ThermoFisher Scientific).Mouse IgG magnetic beads (Catalog No. 5873, Cell Signaling Technology) were used as a negative control to evaluate any non-specific binding on the beads. After incubating with Alex647 labeled Rabbit anti-mouse IgG secondary antibody(Catalog No.A-21239,ThermoFisher Scientific,)and washing,the membranes were visualized with Odyssey? CLx Imaging System.

    Immunocytochemistry

    SNB-19 cells were seeded onto coverslips in 12-well plates 1 day prior to infection.At 24-h post infection,cells were fixed with 4% paraformaldehyde for 15 min at roon temperature,followed by three 10-min washes in PBS at room temperature and permeabilization in PBT (PBS with 0.1% Triton X-100)for 10 min at room temperature. Cells were blocked for 1 h at room temperature in PBTG, incubated with anti-flavivirus group antigen 4G2 (1:1000; catalog No. ATCC? VR-1852TM,ATCC) at 4 °C, washed three times with PBS, and incubated with goat anti-mouse-FITC (1:500; catalog No. AP127F,Sigma-Aldrich)for 1 h at room temperature,followed by three 15-min washes with PBS.Coverslips were mounted and nuclei stained using VECTASHIELD (Catalog No. H-1200, Vector Labs, Burlingame, CA).

    Compound screening using the ZIKV-NS1 TR-FRET assay

    The primary compound screen was performed in 1536-well plates with the ZIKV-NS1 TR-FRET assay as described previously [42]. In total, there are 6016 compounds, including the LOPAC (1280 compounds; catalog No. LO1280, Sigma-Aldrich), the NCGC pharmaceutical collection of 2816 approved drugs, and 1920 bioactive compounds [39].

    For compound screening, HEK293 cells were seeded at 1000 cells/well and incubated at 37 °C with 5% CO2for 16 h. Then, the compounds were transferred to cells in assay plates at 23 nl/well using a pintool workstation (Catalog No.NX-TR pintool station, Wako Automation, San Diego, CA)and incubated for 1 h. ZIKV (MOI = 1) was added to the assay plates at 2 μl/well followed by a 24-h incubation. For detection of NS1 protein levels, 2.5 μl/well of TR-FRET NS1 reagent mixture was added and incubated overnight at 4 °C. The plates were measured in the TR-FRET mode in an EnVision plate reader (Catalog No. 2105-0010, PerkinElmer,Waltham, MA). The experiment for hit compound confirmation was carried out in the same assay as the primary screen except that all the compounds were diluted at a 1:3 ratio for 11 concentrations. The primary screening data and the curve fitting were analyzed as in a previous publication [63]. For the the concentration-response curves and IC50 values of compounds, the confirmation data were analyzed using Prism software (https://www.graphpad.com/, GraphPad Software San Diego, CA).

    Compound cytotoxicity assay

    To eliminate the false positive compounds due to compound cytotoxicity, an ATP content assay [40] was used to measure cell viability after cells were treated with compounds in the absence of ZIKV MR766 infection. Briefly, cells were plated in 1536-well white assay plates in the same way as described above. After a 24-h incubation with compounds, 3.5 μl ATP content reagent mixture (Catalog No. 6016941, PerkinElmer)was added to each well in the assay plates and incubated for 30 min at room temperature.Luminescence signals were determined in a ViewLux plate reader (Catalog No. ViewLuxTMultraHTS microplate imager, PerkinElmer). Compounds with cytotoxicity were eliminated from hit compound list as false positive compounds.

    RNAi screening

    RNAi screening was conducted using the Ambion Silencer?Select Human Druggable Genome siRNA Library Version 4 as described previously[64]and the HTRF assay for NS1 antigen was performed as described above. The HTRF signal for each unique non-overlapping siRNA against the target genes was normalized to a negative control targeting siRNA. The value for each siRNA was divided by the median negative control value and multiplied by 100 to generate the negative normalized metric for each well/siRNA. The median value of negative controls in each plate was used for normalization,while the positive control was set to assess the assay performance and transfection efficiency.

    Data availability

    The primary screening data of all 6016 compounds and the concentration-response curves and IC50 values of compounds with confirmation data were deposited into the PubChem database (PubChem: 1347053), and are publicly accessible at https://pubchem.ncbi.nlm.nih.gov/assay/assay.cgi?aid=1347053.

    CRediT author statement

    Guang Song: Methodology, Validation, Investigation,Resources,Writing-original draft,Writing-review&editing,Visualization. Emily M.Lee: Methodology,Validation, Investigation, Resources, Writing - original draft, Visualization.Jianbo Pan:Methodology,Software,Data curation,Investigation,Visualization,Writing-original draft.Miao Xu:Methodology, Validation, Investigation, Resources, Writing - original draft. Hee-Sool Rho: Investigation, Resources. Yichen Cheng:Investigation, Resources. Nadia Whitt: Investigation,Resources.Shu Yang:Investigation,Resources.Jennifer Kouznetsova: Investigation, Resources. Carleen Klumpp-Thomas:Investigation, Resources. Samuel G. Michael: Investigation,Resources.Cedric Moore:Resources.Ki-Jun Yoon:Resources.Kimberly M. Christian: Resources. Anton Simeonov:Resources. Wenwei Huang: Resources. Menghang Xia:Resources. Ruili Huang: Resources. Madhu Lal-Nag: Investigation, Resources, Validation, Writing - original draft. Hengli Tang:Conceptualization,Writing-original draft,Supervision,Project administration,Funding acquisition.Wei Zheng:Conceptualization, Writing - original draft, Supervision, Project administration, Funding acquisition. Jiang Qian: Conceptualization,Software,Writing-original draft,Supervision,Project administration, Funding acquisition. Hongjun Song: Conceptualization, Writing - original draft, Supervision, Project administration, Funding acquisition. Guo-li Ming: Conceptualization,Writing-original draft,Supervision, Project administration, Funding acquisition. Heng Zhu: Conceptualization,Supervision, Project administration, Funding acquisition,Writing-original draft,Writing-review&editing.All authors read and approved the final manuscript.

    Competing interests

    The authors declare no competing interests.

    Acknowledgments

    The authors would like to thank Paul Shinn and the compound management group at National Center for Advancing Translational Sciences (NCATS) for their professional support,and Jordan Schnoll for coordination.This work was supported by the National Institutes of Health (NIH), USA(Grant Nos. U19AI131130, R01GM111514, R21AI131706,R35NS097370, and R37NS047344) and the Intramural Research Program of the NCATS/NIH, USA. Finally, we would like to thank Jessica Dunn for proofreading this manuscript.

    Supplementary material

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.gpb.2020.06.016.

    ORCID

    0000-0002-7630-716X (Guang Song)

    0000-0001-8451-978X (Emily M. Lee)

    0000-0001-6014-8160 (Jianbo Pan)

    0000-0001-9968-733X (Miao Xu)

    0000-0002-1688-3237 (Hee-Sool Rho)

    0000-0003-3127-0102 (Yichen Cheng)

    0000-0002-2886-9109 (Nadia Whitt)

    0000-0001-8909-1962 (Shu Yang)

    0000-0003-0210-0796 (Jennifer Kouznetsova)

    0000-0002-6646-6132 (Carleen Klumpp-Thomas)

    0000-0002-1070-6563 (Samuel G. Michael)

    0000-0002-3459-392X (Cedric Moore)

    0000-0003-2985-2541 (Ki-Jun Yoon)

    0000-0003-4315-6710 (Kimberly M. Christian)

    0000-0002-4523-9977 (Anton Simeonov)

    0000-0002-7727-9287 (Wenwei Huang)

    0000-0001-7285-8469 (Menghang Xia)

    0000-0001-8886-8311 (Ruili Huang)

    0000-0003-4924-6493 (Madhu Lal-Nag)

    0000-0002-9358-410X (Hengli Tang)

    0000-0003-1034-0757 (Wei Zheng)

    0000-0002-0476-3596 (Jiang Qian)

    0000-0002-8720-5310 (Hongjun Song)

    0000-0002-2517-6075 (Guo-li Ming)

    0000-0002-8426-2889 (Heng Zhu)

    久久午夜福利片| 国产成人一区二区在线| 91久久精品电影网| 婷婷六月久久综合丁香| 最近视频中文字幕2019在线8| 日韩亚洲欧美综合| 亚洲欧美日韩高清在线视频| 一区二区三区激情视频| 国产男人的电影天堂91| 久久午夜亚洲精品久久| 99热这里只有精品一区| 成人高潮视频无遮挡免费网站| 国产av在哪里看| 国产成人福利小说| 国产成人aa在线观看| 极品教师在线视频| 女人被狂操c到高潮| 乱系列少妇在线播放| 天天躁日日操中文字幕| aaaaa片日本免费| 美女被艹到高潮喷水动态| 国产精品亚洲一级av第二区| 美女黄网站色视频| 久久九九热精品免费| 中文在线观看免费www的网站| 老熟妇乱子伦视频在线观看| 丰满的人妻完整版| 久久国内精品自在自线图片| 欧美最黄视频在线播放免费| 男人舔奶头视频| 精品久久国产蜜桃| 久久久久性生活片| netflix在线观看网站| 亚洲精品亚洲一区二区| 国产一区二区三区在线臀色熟女| bbb黄色大片| 亚洲三级黄色毛片| 午夜福利高清视频| 极品教师在线免费播放| 99久国产av精品| 国产免费男女视频| 日韩精品有码人妻一区| 亚洲aⅴ乱码一区二区在线播放| 欧美丝袜亚洲另类 | 日本五十路高清| a级毛片免费高清观看在线播放| av女优亚洲男人天堂| av在线老鸭窝| 女同久久另类99精品国产91| 国产高清视频在线观看网站| 无人区码免费观看不卡| 色5月婷婷丁香| 国产三级在线视频| 亚洲av熟女| 极品教师在线视频| 国产白丝娇喘喷水9色精品| 午夜福利欧美成人| 亚洲va在线va天堂va国产| 国产在视频线在精品| 婷婷精品国产亚洲av| 婷婷六月久久综合丁香| 好男人在线观看高清免费视频| 男女那种视频在线观看| 色哟哟哟哟哟哟| 久久这里只有精品中国| 午夜精品在线福利| 久久久久久久久中文| 国产三级中文精品| 91在线精品国自产拍蜜月| 如何舔出高潮| 日本成人三级电影网站| 91麻豆精品激情在线观看国产| 色精品久久人妻99蜜桃| 联通29元200g的流量卡| 2021天堂中文幕一二区在线观| 欧美一区二区亚洲| 别揉我奶头 嗯啊视频| 亚洲精品一卡2卡三卡4卡5卡| 精品国产三级普通话版| 毛片女人毛片| 免费av观看视频| av在线观看视频网站免费| av在线天堂中文字幕| xxxwww97欧美| 看十八女毛片水多多多| 国产亚洲91精品色在线| 成年免费大片在线观看| 欧美最新免费一区二区三区| 长腿黑丝高跟| 精品人妻熟女av久视频| 天天一区二区日本电影三级| 在线免费观看的www视频| 亚洲av免费高清在线观看| 国产色婷婷99| 色在线成人网| 久久久久九九精品影院| 免费人成在线观看视频色| 欧美一区二区精品小视频在线| 男女边吃奶边做爰视频| 如何舔出高潮| 午夜福利在线观看免费完整高清在 | 午夜日韩欧美国产| 久久精品国产自在天天线| 男女视频在线观看网站免费| 精品久久久噜噜| 欧美一区二区精品小视频在线| 欧美性猛交黑人性爽| 桃色一区二区三区在线观看| 国语自产精品视频在线第100页| 99热网站在线观看| 51国产日韩欧美| 中文资源天堂在线| avwww免费| www.www免费av| 精品久久久噜噜| 又粗又爽又猛毛片免费看| 俄罗斯特黄特色一大片| 高清在线国产一区| 日日啪夜夜撸| 亚洲av熟女| 日韩欧美国产在线观看| 波多野结衣高清无吗| 我要搜黄色片| 老女人水多毛片| 亚洲中文字幕一区二区三区有码在线看| 亚洲熟妇熟女久久| 亚洲欧美激情综合另类| 免费av不卡在线播放| 久久久久久伊人网av| 日本欧美国产在线视频| 国产精品亚洲美女久久久| 精品久久久久久久久亚洲 | 老熟妇仑乱视频hdxx| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 成年女人看的毛片在线观看| 成年女人永久免费观看视频| 99热这里只有精品一区| 动漫黄色视频在线观看| 天美传媒精品一区二区| 国内精品久久久久精免费| 成人国产一区最新在线观看| 一区二区三区四区激情视频 | 91麻豆av在线| 人妻夜夜爽99麻豆av| 最近视频中文字幕2019在线8| 老司机福利观看| 欧美日韩综合久久久久久 | 国产av一区在线观看免费| 国产一区二区三区视频了| 淫秽高清视频在线观看| 精品人妻一区二区三区麻豆 | 亚洲精品成人久久久久久| 欧美成人一区二区免费高清观看| 小蜜桃在线观看免费完整版高清| www.色视频.com| 99精品久久久久人妻精品| 中文资源天堂在线| 亚洲国产欧洲综合997久久,| 91麻豆av在线| 中文字幕av成人在线电影| 亚洲七黄色美女视频| 国产aⅴ精品一区二区三区波| 久久久久久九九精品二区国产| 听说在线观看完整版免费高清| 亚洲男人的天堂狠狠| 欧美日韩黄片免| 嫩草影院精品99| 亚洲人成网站在线播| 亚洲国产高清在线一区二区三| 日韩强制内射视频| 一本精品99久久精品77| 色在线成人网| 最近在线观看免费完整版| 久久久久免费精品人妻一区二区| 国产精品,欧美在线| 热99re8久久精品国产| 亚洲欧美日韩东京热| 99久久中文字幕三级久久日本| 亚洲人与动物交配视频| 亚洲av二区三区四区| 99在线人妻在线中文字幕| 又爽又黄无遮挡网站| 欧美区成人在线视频| 男女边吃奶边做爰视频| 精品久久久久久久久久久久久| 两个人的视频大全免费| 欧美在线一区亚洲| 人人妻人人看人人澡| 午夜福利欧美成人| 亚洲人成伊人成综合网2020| 国产老妇女一区| 亚洲av日韩精品久久久久久密| 精品久久久久久成人av| 波多野结衣高清无吗| 中亚洲国语对白在线视频| 99精品在免费线老司机午夜| 国产精品av视频在线免费观看| 亚洲成人久久爱视频| 成人高潮视频无遮挡免费网站| 国产高清视频在线观看网站| 91精品国产九色| 国产老妇女一区| 嫩草影院精品99| 黄色丝袜av网址大全| 91午夜精品亚洲一区二区三区 | 1000部很黄的大片| 久久亚洲真实| 中文资源天堂在线| 床上黄色一级片| 极品教师在线视频| 久久久午夜欧美精品| 色精品久久人妻99蜜桃| 伊人久久精品亚洲午夜| 国产精品爽爽va在线观看网站| 欧美国产日韩亚洲一区| 国产一区二区三区av在线 | 免费人成视频x8x8入口观看| 91在线观看av| 天堂√8在线中文| 美女高潮的动态| 免费av不卡在线播放| 狂野欧美激情性xxxx在线观看| 久久久久精品国产欧美久久久| 97碰自拍视频| 日韩欧美精品免费久久| 欧美精品啪啪一区二区三区| 成人特级av手机在线观看| 婷婷丁香在线五月| 老熟妇仑乱视频hdxx| 国产精品久久久久久av不卡| 欧美一级a爱片免费观看看| 精品久久久久久,| 久久久久国内视频| 色综合色国产| 深爱激情五月婷婷| 国产伦精品一区二区三区四那| 黄片wwwwww| 午夜爱爱视频在线播放| 99视频精品全部免费 在线| 国产精品爽爽va在线观看网站| 久久精品国产鲁丝片午夜精品 | 直男gayav资源| 大型黄色视频在线免费观看| 给我免费播放毛片高清在线观看| 亚洲专区国产一区二区| 一级毛片久久久久久久久女| 久久精品国产鲁丝片午夜精品 | 久久精品国产鲁丝片午夜精品 | 亚洲人成网站在线播| 国产乱人视频| 午夜老司机福利剧场| 尾随美女入室| a级毛片a级免费在线| 国内久久婷婷六月综合欲色啪| 在线观看av片永久免费下载| 嫁个100分男人电影在线观看| 国产免费男女视频| 露出奶头的视频| 亚洲国产精品久久男人天堂| 色在线成人网| 久久精品国产亚洲网站| 999久久久精品免费观看国产| 嫩草影院新地址| 淫妇啪啪啪对白视频| 少妇丰满av| 亚洲最大成人手机在线| 一个人看的www免费观看视频| 变态另类丝袜制服| 国产精品一区二区免费欧美| 国产精品av视频在线免费观看| 51国产日韩欧美| .国产精品久久| 噜噜噜噜噜久久久久久91| 精品久久久久久久久久久久久| 国产精品福利在线免费观看| 99在线视频只有这里精品首页| 日韩国内少妇激情av| 国产精品人妻久久久影院| 校园春色视频在线观看| 黄片wwwwww| 亚洲人成伊人成综合网2020| 1000部很黄的大片| 免费观看的影片在线观看| 日韩国内少妇激情av| 亚洲av成人av| 乱码一卡2卡4卡精品| 午夜福利欧美成人| 一级av片app| 日韩人妻高清精品专区| 麻豆国产av国片精品| av福利片在线观看| 成人午夜高清在线视频| 18+在线观看网站| 久久久久久伊人网av| 亚洲av不卡在线观看| 国产爱豆传媒在线观看| 久久国内精品自在自线图片| 在线观看免费视频日本深夜| 一区二区三区免费毛片| 看免费成人av毛片| 最近视频中文字幕2019在线8| 午夜精品一区二区三区免费看| 日韩欧美一区二区三区在线观看| 如何舔出高潮| 久久人人精品亚洲av| 国产真实乱freesex| 女同久久另类99精品国产91| 久久久久久久精品吃奶| 看片在线看免费视频| 久久中文看片网| 成熟少妇高潮喷水视频| 国产不卡一卡二| 日韩亚洲欧美综合| 国产精品电影一区二区三区| 亚洲三级黄色毛片| 国产免费男女视频| 午夜免费成人在线视频| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 少妇丰满av| 亚洲精品久久国产高清桃花| 国产成人a区在线观看| 久久热精品热| 观看美女的网站| 亚洲不卡免费看| 悠悠久久av| 中文字幕熟女人妻在线| 又爽又黄无遮挡网站| h日本视频在线播放| 免费不卡的大黄色大毛片视频在线观看 | 国产精品福利在线免费观看| 欧美色欧美亚洲另类二区| av.在线天堂| 性插视频无遮挡在线免费观看| 熟妇人妻久久中文字幕3abv| 少妇裸体淫交视频免费看高清| 日日摸夜夜添夜夜添小说| 国产高清激情床上av| 欧美3d第一页| 免费看a级黄色片| 三级毛片av免费| 精品一区二区三区人妻视频| 精品午夜福利在线看| 国产aⅴ精品一区二区三区波| 久久亚洲精品不卡| 亚洲美女搞黄在线观看 | а√天堂www在线а√下载| 91在线观看av| 深爱激情五月婷婷| 伊人久久精品亚洲午夜| 日日干狠狠操夜夜爽| 91狼人影院| 亚洲性夜色夜夜综合| 看黄色毛片网站| 精品久久国产蜜桃| 波多野结衣巨乳人妻| 九色成人免费人妻av| 高清在线国产一区| 久久久久久伊人网av| 精品久久久久久久久久久久久| 国产主播在线观看一区二区| 尤物成人国产欧美一区二区三区| 久久久久久伊人网av| 国产在视频线在精品| 男女做爰动态图高潮gif福利片| 特大巨黑吊av在线直播| 午夜a级毛片| 99久国产av精品| 18禁黄网站禁片免费观看直播| 麻豆av噜噜一区二区三区| 午夜福利在线在线| 日本黄色片子视频| 亚洲欧美日韩无卡精品| 久久草成人影院| 99国产极品粉嫩在线观看| 国产成人av教育| 国产精品久久久久久亚洲av鲁大| 免费在线观看成人毛片| 人人妻人人澡欧美一区二区| 一级av片app| 无遮挡黄片免费观看| 97超级碰碰碰精品色视频在线观看| 久久久午夜欧美精品| 精品无人区乱码1区二区| 22中文网久久字幕| 中文字幕精品亚洲无线码一区| 久久欧美精品欧美久久欧美| 色播亚洲综合网| 91在线观看av| 美女免费视频网站| 嫩草影视91久久| 日日干狠狠操夜夜爽| 看黄色毛片网站| 成人国产一区最新在线观看| 麻豆国产97在线/欧美| 人妻久久中文字幕网| 一本精品99久久精品77| 嫩草影视91久久| 俺也久久电影网| 国产免费男女视频| 99久久中文字幕三级久久日本| 国内精品久久久久久久电影| 一级黄色大片毛片| 免费不卡的大黄色大毛片视频在线观看 | 在线国产一区二区在线| 香蕉av资源在线| 久久久久性生活片| 亚洲最大成人中文| 乱人视频在线观看| 天堂av国产一区二区熟女人妻| 男女下面进入的视频免费午夜| 午夜免费激情av| 亚洲,欧美,日韩| 99在线人妻在线中文字幕| 日韩欧美国产在线观看| 精华霜和精华液先用哪个| 亚洲精品久久国产高清桃花| netflix在线观看网站| 免费高清视频大片| 亚洲图色成人| 成人国产一区最新在线观看| 制服丝袜大香蕉在线| av天堂在线播放| 十八禁网站免费在线| 三级毛片av免费| av福利片在线观看| 极品教师在线免费播放| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线 | 国产精品98久久久久久宅男小说| 国产男人的电影天堂91| 亚洲欧美精品综合久久99| 亚洲va日本ⅴa欧美va伊人久久| 黄片wwwwww| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 一区二区三区四区激情视频 | 最后的刺客免费高清国语| 日韩,欧美,国产一区二区三区 | 禁无遮挡网站| 亚洲最大成人中文| or卡值多少钱| 女的被弄到高潮叫床怎么办 | 丝袜美腿在线中文| 真实男女啪啪啪动态图| 麻豆久久精品国产亚洲av| 美女cb高潮喷水在线观看| av天堂中文字幕网| 日韩亚洲欧美综合| 波多野结衣高清无吗| 极品教师在线视频| 91麻豆精品激情在线观看国产| 免费在线观看日本一区| 嫩草影院精品99| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久视频播放| 国产免费一级a男人的天堂| 久99久视频精品免费| av中文乱码字幕在线| 欧美成人免费av一区二区三区| 黄色女人牲交| 联通29元200g的流量卡| 18禁在线播放成人免费| 国产高清不卡午夜福利| 国产亚洲av嫩草精品影院| 精品人妻偷拍中文字幕| 蜜桃亚洲精品一区二区三区| 最新中文字幕久久久久| 一本精品99久久精品77| 国产精品爽爽va在线观看网站| 一个人看视频在线观看www免费| 亚洲欧美日韩卡通动漫| 国产精品1区2区在线观看.| 亚洲欧美日韩高清专用| 久久热精品热| 国产69精品久久久久777片| 黄色配什么色好看| 亚洲国产精品成人综合色| 九九在线视频观看精品| 国产成人av教育| 亚洲乱码一区二区免费版| 亚洲久久久久久中文字幕| 亚洲国产欧洲综合997久久,| 亚洲国产日韩欧美精品在线观看| 性插视频无遮挡在线免费观看| 亚洲电影在线观看av| 国内精品久久久久久久电影| 一区福利在线观看| 啦啦啦观看免费观看视频高清| 嫩草影院精品99| 国产精品98久久久久久宅男小说| 麻豆国产av国片精品| 日韩欧美精品v在线| 亚洲av不卡在线观看| 一级a爱片免费观看的视频| 久久草成人影院| 久久欧美精品欧美久久欧美| 久久午夜亚洲精品久久| h日本视频在线播放| 最近最新免费中文字幕在线| h日本视频在线播放| 欧美激情在线99| 伦精品一区二区三区| 国产不卡一卡二| 亚洲欧美精品综合久久99| 国产真实伦视频高清在线观看 | 人妻夜夜爽99麻豆av| 一本一本综合久久| 色哟哟哟哟哟哟| 午夜精品在线福利| 久久精品91蜜桃| 午夜精品在线福利| 999久久久精品免费观看国产| 国产精品嫩草影院av在线观看 | АⅤ资源中文在线天堂| 国产aⅴ精品一区二区三区波| 久久久成人免费电影| 国产在视频线在精品| 久久婷婷人人爽人人干人人爱| 校园春色视频在线观看| 草草在线视频免费看| 久久精品国产亚洲av天美| 免费av毛片视频| 又爽又黄a免费视频| 国产精品久久久久久精品电影| 日韩欧美在线二视频| 搡女人真爽免费视频火全软件 | 亚洲成a人片在线一区二区| 国内揄拍国产精品人妻在线| 精品久久久噜噜| 我要搜黄色片| 国产午夜福利久久久久久| 欧美色欧美亚洲另类二区| 国产91精品成人一区二区三区| 色综合婷婷激情| 性插视频无遮挡在线免费观看| 国产精品综合久久久久久久免费| 91麻豆av在线| 国产亚洲91精品色在线| 欧美激情久久久久久爽电影| 国产欧美日韩一区二区精品| 97人妻精品一区二区三区麻豆| 乱码一卡2卡4卡精品| 18禁在线播放成人免费| 熟女电影av网| 一级av片app| 88av欧美| 成人高潮视频无遮挡免费网站| 久久久色成人| 黄片wwwwww| 日日夜夜操网爽| 免费人成视频x8x8入口观看| 国产爱豆传媒在线观看| 别揉我奶头 嗯啊视频| 色综合站精品国产| 亚洲精品一卡2卡三卡4卡5卡| 性色avwww在线观看| 亚洲一区二区三区色噜噜| 久9热在线精品视频| 亚洲一区二区三区色噜噜| 麻豆成人午夜福利视频| 麻豆久久精品国产亚洲av| 色综合站精品国产| 内射极品少妇av片p| 毛片一级片免费看久久久久 | av国产免费在线观看| 又爽又黄a免费视频| 精品午夜福利在线看| 国产av麻豆久久久久久久| 日韩欧美免费精品| 亚洲国产精品合色在线| 夜夜爽天天搞| 亚洲电影在线观看av| 国产麻豆成人av免费视频| 国产精品福利在线免费观看| 精品日产1卡2卡| 日韩欧美一区二区三区在线观看| 老师上课跳d突然被开到最大视频| 美女黄网站色视频| 亚洲av中文字字幕乱码综合| 欧美日韩乱码在线| 婷婷亚洲欧美| 午夜福利视频1000在线观看| 97碰自拍视频| 亚洲久久久久久中文字幕| 搡女人真爽免费视频火全软件 | 亚洲不卡免费看| 亚州av有码| 亚洲成av人片在线播放无| 91在线精品国自产拍蜜月| 亚洲国产高清在线一区二区三| 能在线免费观看的黄片| 色在线成人网| 99热这里只有精品一区| 国产精品一及| 真人做人爱边吃奶动态| 免费大片18禁| 国产毛片a区久久久久| 国内精品宾馆在线| 亚洲av中文av极速乱 | 国产三级在线视频| 久久精品国产鲁丝片午夜精品 | av专区在线播放| 欧美日韩乱码在线| 一级黄色大片毛片| 三级毛片av免费| 3wmmmm亚洲av在线观看| 亚洲国产精品久久男人天堂| 久久精品影院6| 伦精品一区二区三区| 亚洲va在线va天堂va国产| 男人舔奶头视频| 精品久久久久久成人av|