• <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)

    国产精品福利在线免费观看| 99热全是精品| 免费少妇av软件| 欧美激情极品国产一区二区三区 | 18禁在线无遮挡免费观看视频| 少妇人妻精品综合一区二区| 国产大屁股一区二区在线视频| 日韩成人伦理影院| 欧美精品一区二区大全| xxx大片免费视频| 天天躁日日操中文字幕| 亚洲国产成人一精品久久久| 黄色一级大片看看| 精华霜和精华液先用哪个| 亚洲内射少妇av| 91久久精品国产一区二区成人| 日韩 亚洲 欧美在线| 国产深夜福利视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 男人舔奶头视频| 各种免费的搞黄视频| 亚洲欧洲日产国产| 国产老妇伦熟女老妇高清| 国产在线一区二区三区精| 欧美日本视频| 自拍偷自拍亚洲精品老妇| 国产精品免费大片| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 日本黄色日本黄色录像| 日本色播在线视频| 丰满迷人的少妇在线观看| 亚洲国产色片| 成年女人在线观看亚洲视频| 国产av码专区亚洲av| 观看av在线不卡| 国产视频首页在线观看| 亚洲欧美一区二区三区国产| 亚洲人成网站在线观看播放| 国产爱豆传媒在线观看| 欧美另类一区| 建设人人有责人人尽责人人享有的 | 赤兔流量卡办理| 永久免费av网站大全| 色吧在线观看| 女性被躁到高潮视频| 最近最新中文字幕大全电影3| 黄色日韩在线| 亚洲av成人精品一二三区| 久久青草综合色| 啦啦啦啦在线视频资源| 久久鲁丝午夜福利片| 欧美日韩在线观看h| 我要看日韩黄色一级片| 国产精品国产三级国产专区5o| 日韩一区二区视频免费看| 人体艺术视频欧美日本| 在线观看三级黄色| 国产精品熟女久久久久浪| 丝袜喷水一区| av网站免费在线观看视频| 一级片'在线观看视频| 国内精品宾馆在线| 欧美97在线视频| 免费观看在线日韩| 国国产精品蜜臀av免费| 少妇人妻精品综合一区二区| 国产一区二区三区综合在线观看 | 少妇精品久久久久久久| 亚洲欧美一区二区三区黑人 | 亚洲精品日韩av片在线观看| 搡老乐熟女国产| 免费黄色在线免费观看| 国产一区二区三区综合在线观看 | 夜夜爽夜夜爽视频| freevideosex欧美| 舔av片在线| 国产精品国产三级国产专区5o| 国产伦理片在线播放av一区| 最新中文字幕久久久久| 亚洲欧美清纯卡通| 国产av一区二区精品久久 | 丰满少妇做爰视频| 日日啪夜夜爽| 亚洲国产高清在线一区二区三| 色吧在线观看| 在线免费十八禁| 成人午夜精彩视频在线观看| 久久久亚洲精品成人影院| 黑人猛操日本美女一级片| 亚洲人与动物交配视频| 99国产精品免费福利视频| 一区在线观看完整版| 少妇裸体淫交视频免费看高清| 免费看av在线观看网站| 成人亚洲精品一区在线观看 | 国产精品伦人一区二区| 久久久久久久久久人人人人人人| 18禁裸乳无遮挡动漫免费视频| 亚洲四区av| 国产淫语在线视频| 日本午夜av视频| 卡戴珊不雅视频在线播放| 国产深夜福利视频在线观看| av网站免费在线观看视频| av一本久久久久| 青春草国产在线视频| 伦精品一区二区三区| 成人黄色视频免费在线看| 一区二区三区四区激情视频| 久久久a久久爽久久v久久| 大话2 男鬼变身卡| 在线观看美女被高潮喷水网站| 久热久热在线精品观看| 黑丝袜美女国产一区| 国产精品国产三级国产av玫瑰| 97热精品久久久久久| 国产成人a区在线观看| 久久久久久九九精品二区国产| 精品一区在线观看国产| 99热这里只有是精品在线观看| 欧美+日韩+精品| 久久精品国产自在天天线| 嫩草影院入口| 美女中出高潮动态图| 亚洲av国产av综合av卡| 91精品伊人久久大香线蕉| 老师上课跳d突然被开到最大视频| 免费观看av网站的网址| 一级二级三级毛片免费看| 极品教师在线视频| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 大片电影免费在线观看免费| 久久久久性生活片| 色网站视频免费| 又粗又硬又长又爽又黄的视频| 免费黄色在线免费观看| 日本黄色日本黄色录像| 美女国产视频在线观看| 国产精品嫩草影院av在线观看| 80岁老熟妇乱子伦牲交| 日韩不卡一区二区三区视频在线| 国产亚洲精品久久久com| 午夜免费观看性视频| 亚洲精品,欧美精品| 一区二区三区四区激情视频| 精品一品国产午夜福利视频| 美女主播在线视频| 丝袜脚勾引网站| 中国美白少妇内射xxxbb| 3wmmmm亚洲av在线观看| 一本久久精品| 国产色婷婷99| 国产 一区精品| 久久韩国三级中文字幕| 在线观看一区二区三区激情| 亚洲美女黄色视频免费看| 王馨瑶露胸无遮挡在线观看| 欧美一级a爱片免费观看看| 亚洲精品久久久久久婷婷小说| 日本爱情动作片www.在线观看| 亚洲精品乱码久久久v下载方式| 天堂中文最新版在线下载| 欧美国产精品一级二级三级 | 日韩,欧美,国产一区二区三区| 男女下面进入的视频免费午夜| 亚洲精品国产av成人精品| 秋霞伦理黄片| 亚洲美女视频黄频| 男人爽女人下面视频在线观看| 国产精品不卡视频一区二区| 亚洲一级一片aⅴ在线观看| 午夜激情福利司机影院| 黄色一级大片看看| 亚洲精品国产成人久久av| 在线亚洲精品国产二区图片欧美 | 久久国产精品大桥未久av | 中文乱码字字幕精品一区二区三区| 熟女av电影| 成人国产麻豆网| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| 成人毛片a级毛片在线播放| 激情 狠狠 欧美| 美女视频免费永久观看网站| 亚洲成色77777| 99久久中文字幕三级久久日本| 女性被躁到高潮视频| 精品亚洲成a人片在线观看 | 一区二区三区乱码不卡18| 一个人免费看片子| 久久久午夜欧美精品| 日韩av免费高清视频| 久久 成人 亚洲| 久久久久久久久久久丰满| 日韩电影二区| 亚洲精品乱码久久久久久按摩| 欧美一区二区亚洲| 一个人免费看片子| 亚洲精品久久久久久婷婷小说| 亚洲国产精品一区三区| 老师上课跳d突然被开到最大视频| 成人毛片60女人毛片免费| 亚洲av电影在线观看一区二区三区| 亚洲成人一二三区av| 97精品久久久久久久久久精品| 中国美白少妇内射xxxbb| 男女下面进入的视频免费午夜| 日本欧美国产在线视频| 亚洲欧美日韩东京热| 成人国产麻豆网| 一本一本综合久久| 亚洲欧美日韩东京热| 菩萨蛮人人尽说江南好唐韦庄| 一本一本综合久久| 久久鲁丝午夜福利片| 十分钟在线观看高清视频www | 亚洲欧美日韩卡通动漫| 美女中出高潮动态图| 中文欧美无线码| 性色avwww在线观看| 亚洲丝袜综合中文字幕| 日韩中字成人| videossex国产| 成年女人在线观看亚洲视频| 亚洲成人中文字幕在线播放| 丝袜脚勾引网站| 伊人久久精品亚洲午夜| 麻豆国产97在线/欧美| 亚洲欧美中文字幕日韩二区| 十八禁网站网址无遮挡 | 97在线视频观看| 男人狂女人下面高潮的视频| 在线看a的网站| 青春草亚洲视频在线观看| 免费黄网站久久成人精品| 成年美女黄网站色视频大全免费 | 永久免费av网站大全| 欧美日韩视频高清一区二区三区二| 亚洲人成网站在线播| 波野结衣二区三区在线| 2018国产大陆天天弄谢| 国产v大片淫在线免费观看| 肉色欧美久久久久久久蜜桃| 99久久综合免费| av视频免费观看在线观看| 亚洲国产精品专区欧美| 欧美成人精品欧美一级黄| 两个人的视频大全免费| 亚洲综合精品二区| 亚洲av.av天堂| 一区二区三区免费毛片| 尤物成人国产欧美一区二区三区| 午夜福利网站1000一区二区三区| 观看美女的网站| 晚上一个人看的免费电影| 亚洲av二区三区四区| 亚洲伊人久久精品综合| 99久久综合免费| 又黄又爽又刺激的免费视频.| 免费观看av网站的网址| av免费在线看不卡| av国产久精品久网站免费入址| 直男gayav资源| 人人妻人人澡人人爽人人夜夜| av免费观看日本| 最黄视频免费看| 国产av国产精品国产| 久久久久久久国产电影| 国产黄片美女视频| 亚洲美女视频黄频| 日韩欧美精品免费久久| 精品一区二区三卡| 午夜福利在线在线| 成人18禁高潮啪啪吃奶动态图 | 又爽又黄a免费视频| 国产在线视频一区二区| 欧美日韩精品成人综合77777| 啦啦啦中文免费视频观看日本| 美女视频免费永久观看网站| 国产欧美亚洲国产| 亚洲欧美一区二区三区黑人 | 日本黄大片高清| 日韩三级伦理在线观看| 简卡轻食公司| 内地一区二区视频在线| 久久久久久久久久人人人人人人| 久久久久视频综合| 日本一二三区视频观看| 国产在线男女| 免费观看a级毛片全部| 国产免费视频播放在线视频| av黄色大香蕉| 国产有黄有色有爽视频| 成人国产av品久久久| 99久久中文字幕三级久久日本| 国产一级毛片在线| 亚洲真实伦在线观看| 在线免费十八禁| av天堂中文字幕网| 91精品国产国语对白视频| 亚洲国产精品999| 又粗又硬又长又爽又黄的视频| 伊人久久国产一区二区| 99久久精品热视频| 欧美成人午夜免费资源| 亚洲av成人精品一二三区| 亚洲精品中文字幕在线视频 | 美女主播在线视频| 国产v大片淫在线免费观看| 午夜福利在线在线| 麻豆乱淫一区二区| 男女边摸边吃奶| 毛片女人毛片| 亚洲第一区二区三区不卡| 美女福利国产在线 | 精品一区在线观看国产| 国产av一区二区精品久久 | 久久精品国产鲁丝片午夜精品| 精华霜和精华液先用哪个| 建设人人有责人人尽责人人享有的 | 亚洲精品色激情综合| 内地一区二区视频在线| 草草在线视频免费看| 美女视频免费永久观看网站| av免费观看日本| 色视频www国产| 午夜日本视频在线| 久久久久性生活片| 国产精品偷伦视频观看了| 97精品久久久久久久久久精品| 午夜福利网站1000一区二区三区| 亚洲国产色片| 天天躁夜夜躁狠狠久久av| 国产片特级美女逼逼视频| 免费观看av网站的网址| 久久久色成人| 亚洲av日韩在线播放| 一本久久精品| 国模一区二区三区四区视频| 欧美老熟妇乱子伦牲交| 国产精品女同一区二区软件| 日日摸夜夜添夜夜添av毛片| 成人漫画全彩无遮挡| 小蜜桃在线观看免费完整版高清| 伦理电影免费视频| 日日摸夜夜添夜夜添av毛片| 欧美一区二区亚洲| 久久精品国产亚洲av涩爱| 国产精品精品国产色婷婷| 日韩电影二区| 激情五月婷婷亚洲| 免费看不卡的av| 国产亚洲午夜精品一区二区久久| 亚洲人成网站在线观看播放| 我要看日韩黄色一级片| 欧美高清成人免费视频www| 免费久久久久久久精品成人欧美视频 | 色网站视频免费| 精品久久久久久久末码| 国产av码专区亚洲av| 丰满少妇做爰视频| 亚洲第一区二区三区不卡| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 好男人视频免费观看在线| 岛国毛片在线播放| 国产成人a∨麻豆精品| 成人免费观看视频高清| 六月丁香七月| 亚洲成人av在线免费| 久久青草综合色| 日本黄色日本黄色录像| 国产成人a∨麻豆精品| 三级经典国产精品| 欧美xxxx黑人xx丫x性爽| 国产高清不卡午夜福利| 久久国内精品自在自线图片| 人妻夜夜爽99麻豆av| 99久国产av精品国产电影| 日本-黄色视频高清免费观看| 亚洲综合精品二区| 大话2 男鬼变身卡| 最近2019中文字幕mv第一页| 中文字幕av成人在线电影| 欧美xxxx性猛交bbbb| 日韩人妻高清精品专区| 青春草国产在线视频| 国产精品福利在线免费观看| 欧美日韩亚洲高清精品| 一区二区av电影网| tube8黄色片| 国产男人的电影天堂91| 久久久精品免费免费高清| 夫妻午夜视频| 欧美老熟妇乱子伦牲交| 一区二区三区四区激情视频| 夫妻性生交免费视频一级片| 婷婷色综合www| 寂寞人妻少妇视频99o| 男人爽女人下面视频在线观看| 男的添女的下面高潮视频| 亚洲人与动物交配视频| 亚洲美女黄色视频免费看| 91精品国产九色| 婷婷色综合大香蕉| 国产在线一区二区三区精| 日日撸夜夜添| 久久精品熟女亚洲av麻豆精品| 国国产精品蜜臀av免费| 人人妻人人爽人人添夜夜欢视频 | 中文字幕精品免费在线观看视频 | 午夜激情福利司机影院| 亚洲欧美日韩无卡精品| 亚洲精品乱码久久久v下载方式| www.色视频.com| 国产大屁股一区二区在线视频| 久久久久视频综合| 国产在线男女| 精品少妇久久久久久888优播| 国产成人精品福利久久| 午夜福利影视在线免费观看| 国产免费一级a男人的天堂| a 毛片基地| 日韩av不卡免费在线播放| 亚洲欧美日韩另类电影网站 | 人妻少妇偷人精品九色| 欧美亚洲 丝袜 人妻 在线| 成人二区视频| 精品一区二区免费观看| 黄色欧美视频在线观看| 激情 狠狠 欧美| 不卡视频在线观看欧美| 成人无遮挡网站| 天堂中文最新版在线下载| 简卡轻食公司| 日本av免费视频播放| 99热这里只有精品一区| 亚洲精品色激情综合| 久久青草综合色| 黄色配什么色好看| 最后的刺客免费高清国语| 国产成人精品久久久久久| 欧美精品亚洲一区二区| 国产 一区精品| 国产深夜福利视频在线观看| 亚洲自偷自拍三级| 午夜免费男女啪啪视频观看| 校园人妻丝袜中文字幕| 人体艺术视频欧美日本| 精品人妻熟女av久视频| 国产午夜精品久久久久久一区二区三区| 中文字幕亚洲精品专区| 各种免费的搞黄视频| 免费人成在线观看视频色| 激情五月婷婷亚洲| 色吧在线观看| 亚洲av在线观看美女高潮| 男女下面进入的视频免费午夜| 成年av动漫网址| 亚洲美女视频黄频| 尤物成人国产欧美一区二区三区| 国产成人免费观看mmmm| 高清在线视频一区二区三区| 最近最新中文字幕免费大全7| 国产高清三级在线| 在线观看三级黄色| 一级毛片aaaaaa免费看小| 黑人猛操日本美女一级片| 国产一区有黄有色的免费视频| 亚洲美女视频黄频| 色视频www国产| 久久久午夜欧美精品| 国产乱来视频区| 亚洲精品乱码久久久久久按摩| 狂野欧美激情性bbbbbb| 人妻夜夜爽99麻豆av| 18禁裸乳无遮挡免费网站照片| 成人免费观看视频高清| 在线观看美女被高潮喷水网站| 免费人成在线观看视频色| 日韩中字成人| 精品一区二区三区视频在线| 中文字幕av成人在线电影| 99热6这里只有精品| 男女啪啪激烈高潮av片| 久久精品久久久久久噜噜老黄| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 亚洲国产精品成人久久小说| 天堂中文最新版在线下载| 插阴视频在线观看视频| 三级经典国产精品| 91在线精品国自产拍蜜月| 色5月婷婷丁香| 激情 狠狠 欧美| 97超碰精品成人国产| 高清日韩中文字幕在线| 日本欧美视频一区| 欧美国产精品一级二级三级 | av黄色大香蕉| 国产一区二区在线观看日韩| a级毛色黄片| 日本色播在线视频| 亚洲美女搞黄在线观看| av国产免费在线观看| 国产一级毛片在线| 99热国产这里只有精品6| 女的被弄到高潮叫床怎么办| 欧美亚洲 丝袜 人妻 在线| 美女视频免费永久观看网站| 草草在线视频免费看| 特大巨黑吊av在线直播| 少妇被粗大猛烈的视频| 91精品国产九色| 中文欧美无线码| 久久久久久久国产电影| 高清日韩中文字幕在线| av黄色大香蕉| 久久人人爽人人爽人人片va| 国产精品人妻久久久久久| 国产精品一区二区在线不卡| 色5月婷婷丁香| 国产精品精品国产色婷婷| 久久久精品免费免费高清| 综合色丁香网| 亚洲经典国产精华液单| 国产精品久久久久久av不卡| 久久女婷五月综合色啪小说| 欧美成人a在线观看| 黄色日韩在线| 欧美成人一区二区免费高清观看| 亚洲成人av在线免费| 国产淫片久久久久久久久| 黑丝袜美女国产一区| 亚洲欧美清纯卡通| 免费观看av网站的网址| 能在线免费看毛片的网站| 永久网站在线| 简卡轻食公司| 三级国产精品欧美在线观看| 水蜜桃什么品种好| 日日啪夜夜撸| 纯流量卡能插随身wifi吗| 亚洲伊人久久精品综合| 国产精品麻豆人妻色哟哟久久| 亚州av有码| 联通29元200g的流量卡| 精品久久久久久久久av| 亚洲丝袜综合中文字幕| 视频中文字幕在线观看| 中文资源天堂在线| 天美传媒精品一区二区| 大香蕉久久网| 成人无遮挡网站| 成年美女黄网站色视频大全免费 | 黄片无遮挡物在线观看| 成年av动漫网址| 国产成人91sexporn| 老女人水多毛片| 有码 亚洲区| 免费看不卡的av| 中国国产av一级| 国产精品一区二区三区四区免费观看| 亚洲va在线va天堂va国产| 1000部很黄的大片| 国产精品99久久99久久久不卡 | 香蕉精品网在线| 亚洲一级一片aⅴ在线观看| 国产淫语在线视频| 国产精品女同一区二区软件| 久久女婷五月综合色啪小说| 中文乱码字字幕精品一区二区三区| 久久久午夜欧美精品| 久热久热在线精品观看| 亚洲国产av新网站| 国产淫语在线视频| 成人漫画全彩无遮挡| 国产一区二区三区综合在线观看 | 国产大屁股一区二区在线视频| 国产精品三级大全| 久久毛片免费看一区二区三区| 97在线视频观看| 国产在线免费精品| 我要看黄色一级片免费的| av在线蜜桃| 国产精品久久久久久精品古装| www.色视频.com| 久久国产精品男人的天堂亚洲 | 国产成人aa在线观看| 国产精品福利在线免费观看| 极品教师在线视频| 久久久成人免费电影| 国产精品麻豆人妻色哟哟久久| 国产亚洲最大av| 国产伦在线观看视频一区| 好男人视频免费观看在线| 国产乱人视频| 欧美高清性xxxxhd video| 日本猛色少妇xxxxx猛交久久| 欧美激情国产日韩精品一区| 亚洲精品国产av成人精品| 国产精品国产三级专区第一集| 国产乱来视频区| 欧美日韩精品成人综合77777| 尤物成人国产欧美一区二区三区| 国产男女内射视频| 伊人久久精品亚洲午夜| 卡戴珊不雅视频在线播放| .国产精品久久| av女优亚洲男人天堂| 免费高清在线观看视频在线观看|