1. Introduction
Human papillomavirus (HPV) is a small, non-enveloped DNA virus that infects skin or mucosal cells. The circular, double-stranded viral genome is approximately 8 kb in length. It is estimated that 20% of human cancers are caused by infections, and that approximately 12% of the total human cancer burden has a viral association [
1]. HPV causes 5% of all human cancers and approximately 100% of all cervical cancers [
2,
3]. HPV is also responsible for a very high percentage of other anogenital cancers, and up to 70% of head-and-neck cancers, which includes those of the oral cavity, oropharynx, sinus, tonsil, and larynx [
4,
5]. Overall, over 600,000 new cancer cases every year are caused by high-risk HPV types, with types 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 being defined by WHO as cancer-causing in the cervix [
6]. In addition, HPV infection is the most common viral sexually transmitted infection worldwide [
7,
8,
9].
HPV virions are composed of two structural proteins L1 and L2, both of which are required for a successful viral entry, L1 is required for cell attachment and L2 is required for successful viral genome trafficking to the cell nucleus [
10,
11,
12]. From previous studies we know that L1 is the major capsid protein [
13] and L2 is the minor capsid proteins [
14,
15], Previous research has shown that HPV can enter cells through clathrin-dependent or non-clathrin-dependent pathways [
16,
17]; it then travels to the early endosome, and from the early endosome to the late endosome [
18]. During the acidification process the viral capsid begins a process of disassembly [
19] and this allows a small portion of L2 to become exposed and span the endosomal membrance into the cytosol [
20]. This exposure of L2 sequences allows the viral protein to recruit various cargo sorting components [
21]. This in turn promotes endosomal tubulation and trafficking of the L2 DNA complex away from the late endosome [
22,
23], whilst the majority of L1 is left behind and is subsequently degraded within the lysozome [
24,
25]. During this process the L2/DNA complex is transported to the trans-Golgi network (TGN) [
26,
27], and stays there until mitosis [
28,
29], when the nuclear envelope breaks down and thereby allows access to nuclear components and subsequent nuclear retention [
25,
30]. Whilst a great deal is already known about the different cellular partners which facilitate this entry process and trafficking steps within the cell, such as SNX17 [
31], CCT-3 [
32], SNX27 and proteins of the retromer machinery [
12,
33], we were nonetheless interested in investigating other possible L2 interacting partners which may play a role in the infectious entry process. To do this we performed a proteomic analysis of L2 interacting proteins. Amongst these, we found that, we have found that COPB1, a component of the COP1 cargo protein complex is an interacting partner of L2 during HPV infectious entry; furthermore HPV requires. We have found that HPV requires COPI cargo proteins for a successful infection.
As from previous studies we know that there are two types of COP protein COPI and COPII [
34,
35,
36], they are composed of five subunits COPI is composed of (α-COP), (β’-COP), (ε-COP (β-COP), (γ-COP), (δ-COP), (ζ-COP) and COPII is composed of Sar1, Sec23, Sec24, Sec13, and Sec31 [
34,
37], each subunit perform a specific function [
38], COPII proteins transport newly synthesized proteins from ER endoplasmic reticulum to the TGN, [
39,
40,
41] While COPI proteins transports recycling proteins from cis Golgi networks to the ER [
42,
43], Recently, it has also been shown that COPI protein is involved in trafficking of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [
44,
45] and Nairoviruses [
46]. Basically both COP proteins are required for transportation of large molecules from one enclosed compartment to another, formation of membrane domains within the cells [
47,
48], being large transport carrier proteins, their connection with the cytoskeleton and interaction with the target organelles are the basic functions of COP proteins [
49,
50,
51].
2. Materials and Methods
2.1. Mass Spectrometry Analysis
FLAG-HA-tagged HPV-16 L2 expression plasmid or, as a control, empty plasmid were transfected into HEK293 cells. After 48 hours, cell extracts were produced using mass spectrometry lysis buffer (50 mM HEPES [pH 7.4], 150 mM NaCl, 50 mM NaF, 1 mM EDTA, 0.5% NP-40), and they were then incubated for one hour at 4°C on a rotating wheel with EZview Red Anti-HA affinity gel beads (Sigma). The beads were washed using phosphate-buffered saline (PBS) after first being washed with wash buffer (25 mM HEPES [pH 7.0], 0.1% NP-40, and 150 mM NaCl), and then proteome analysis was performed as previously described [
52]
In short, 50 ng of sequencing-grade trypsin in 20 mM triethanolamine bicarbonate (pH 8.5) was used for 12 hours at room temperature to elute proteins from the affinity beads. Following digestion, the beads’ supernatant was extracted, and they were then quickly incubated with another aliquot of 20 mM triethanolamine bicarbonate (pH 8.5). This aliquot was then extracted and combined with the first supernatant. Formic acid was added to 0.1% to halt the reactions. The resultant mixture was then desalted using C18 STAGE (Stop and Go Extraction) tips and lyophilized until it was completely dry. Nanobore columns were constructed using in-house constructed columns packed with 15 cm of 3-μm Ascentis RPA particles (Sigma) using a high-pressure column loader. After being desalted, the samples were reconstituted in 0.1% formic acid and put onto the column. A discontinuous gradient of 0% to 80% acetonitrile in 0.1% formic acid was used to create the column 60 minutes and poured straight into the Amazon ETD ion trap mass spectrometer’s opening (Bruker). During the liquid chromatography (LC) separation, eight data-dependent tandem mass spectrometry (MS/MS) scans were conducted in succession, after one complete scan (375 to 1,700 m/z). Using the Global Proteome Machine interfaced with X!Tandem, raw data files from the Amazon ETD instrument were converted to mgf files by a data analysis software package (Bruker). Then, the files were searched against the Ensembl human protein database and the NCBInr Viral database. The matches were filtered at a 1% false discovery rate.
2.2. Cell Culture and Transfections
The human embryonic kidney cells (HEK293 and HEK293TT-21), and human skin keratinocytes (HaCaT) were cultured in Dulbecco’s modified Eagle medium, supplemented with 10% foetal bovine serum, 100 ug/ml penicillin-streptomycin, and 300 mg/ml glutamine.
HEK293 and HEK293TT cells were transfected using the calcium phosphate precipitation technique (22). For siRNA experiments, HaCaT cells were transfected using Lipofectamine RNAiMAX (Invitrogen) with an ON-TARGETplus SMARTpool of siRNAs against COPB1 (Dharmacon): (5′-CCAAGAUUGCAUUGCGCUA), (5′-CAUAUAAGAAUUCGUGCAA), (5′-AUUAUUAAGGAGAGCGACA), and (5′-CCUCAUGACUUCGCAAAUA). In parallel, a scrambled siSTABLE nontargeting siRNA (Dharmacon) was employed as a control.
2.3. Antibodies
The following primary antibodies were employed in this study: Mouse anti-HPV-16 L1 ; Santa Cruz), Cat #sc-47699) 1:1000 dilutions, with rabbit anti-COPB1 (Abcam) Cat # ab205020]1:1500 dilutions, In Western blotting tests, the secondary antibodies horseradish peroxidase-conjugated swine anti-rabbit (Cat #PO399), and rabbit anti-mouse (Dako), Cat #DK25 (Dako) 1:5000 dilutions, as well as the anti-HA mouse (Santa Cruz, Cat # sc-7392), 1:2000 dilution), HA rabbit (Thermofisher, Cat #MA5-27915), 1:2000 dilution, were used.
2.4. Plasmids
The plasmids expressing FLAG-HA-tagged HPV-16 L2 have been previously described [
22], as have the expression vectors for the HPV-16 L1 and HPV-16 L2 full-length GST fusion proteins [
52]. DNA sequencing was used to authenticate the plasmids produced by the GENEART site-directed mutagenesis system (Invitrogen) for the production of shortened versions of GST-tagged HPV-16 L2, which were made using custom-designed oligonucleotides (Eurofins MWG). HPV-16 L1 and L2 were expressed in a codon-optimized bicistronic fashion using plasmid p16shell L2-3×FLAG-HA, in which L2 is tagged with FLAG and hemagglutinin (HA) tags [
53]. The firefly luciferase gene-carrying pGL3 Luci construct was obtained from Promega. Genscript was used to create the murine Papilloma virus L2 protein construct. The plasmids expressing HPV-2, HPV-5, BPV-1, SfPV-1, Mus musculus papillomavirus type 1 (MmuPV-1), and Merkel Cell Polyoma Virus (MCPyV) capsid proteins were the gift of Dr Chris B. Buck; HPV-16 and HPV-31 capsid protein expression plasmids were the gift of Dr Michelle A. Ozbun; HPV-18 capsid protein expression plasmid was the gift of Dr Samuel K. Campos, the GST BPV-1 L2 construct was the gift of Dr Patricio I. Meneses, and the GST HPV-11 L2 construct was the gift of Dr Bob Garcea, all of which were sequenced for verification.
2.5. Fusion Protein Purification and Binding Assays
As previously reported [
54], GST-tagged fusion proteins or GST alone were produced, purified, and immobilized on glutathione agarose (Sigma). After thorough washing, the GST bead-bound proteins were identified by western blotting with selected primary and secondary antibodies and by Ponceau staining. HaCat cell lysates were treated with GST-tagged fusion proteins for one hour at 4°C.
2.6. Infectivity Assays
In order to assess the infectivity in HaCaT cells, they were treated for 48 hours with either a control siRNA (siCTRL) or a siRNA specifically targeting COPB1. HPV-16 PsVs were added to these cells at concentrations of 150 vge/cell. Using a luciferase assay equipment (Promega), firefly luciferase activity was measured 48 hours later as an indicator of infection. To ensure that the luciferase assays used equivalent protein inputs, the total cell protein was measured. A one-tailed t test was used to assess the statistical significance of the decrease in luciferase activity upon si-COPB1-treated cells, relative to the luciferase readings for siCTRL-treated cells in each experiment. Additionally, the cellular lysates were utilised to validate COPB1’s knockdown effectiveness by western blot. 25*10^5 HEK-293 TT cells were seeded in 15 cm dishes (10 dishes) and incubated overnight at 37°C, the cells were overexpressed with (15ug/15cm plate) HPV p16 Shell (L1and L2 capsid proteins of HPV) with (30ug/ 15cm plate) Luci gene (reporter of capsid proteins), for 48 hours, then the cells were harvested and the pellet were transferred in to siliconized tubes, and 9.5M magnesium chloride were added the same as the volume of pellet, 10% Triton in PBS were added 1/10 volume of pellet, exonuclease and endonuclease enzymes were also added, then we vortexed it and incubate it for overnight at 37°C. then 5M NaCl were added, and the tube were freeze-thawed three times in dry ice, and then the supernatant were collected and we washed the pellet with 1ml HBS, two times, and mixed the washed supernatant with the lysates supernatant. Then a gradient of CSCL were prepared in ultracentrifuge tubes, and we added the lysate on the top and spin for overnight (16 hrs) at 20,000-RPM at 4C. then we got a band in the middle of the tube between heavy and light CSCL in the middle of the tube which contained the HPV PsVs. The band were extracted with a syringe and transferred in to filter containing falcon tubes and washed three times with HBS buffer, and then we runned a 2 ul 2ul of the PsVs in a SDS-PAGE and the gel were stained with Coomassie brilliant blue for 15 minutes and distained over-night with de-stained buffer, and L1 and L2 proteins were detected. And we measured the protein concentration with nanodrop.
Pseudovirions Quantification
LUCI DNA was extracted from PsVs by incubating 10 ul of each PsVs preparation with 5 ul of extraction buffer (20 mM Tris (pH 8.0), 20 mM DTT, 20 mM EDTA, 2.0% SDS, and 0.2% Proteinase K) at 50°C for 15 min. Then, 1.5 ul of the extract was run on an agarose gel to confirm extraction of reporter genome (Luci).
The extracted DNA was purified using a PCR purification kit and eluted by adding 30 ul of nuclease free water.
A standard curve was made for qPCR, as follows: standard stock 2*109 plasmid/µl was used. To make this, 1-2 µg of PGL3 plasmid was digested with BglII enzyme and run on an agarose gel. After excision and purification of the linear plasmid, it was quantified by nanodrop and the copy number (plasmid/ng) determined.
The qPCR was set up with different concentrations using a standard and DNA from the PsVs; we used single digested PGL3 plasmid with BglII as a positive control and water as negative control and every sample was analysed in duplicate.
Next, the VGE and MOI were calculated, based on the qPCR results and the cell numbers that were used; the calculations were done as below. 10: represents the amount of PsV preparation [in microliters] used for the genome extraction. The “ul of prep. used” is the amount of PsV per well used in the infection (0.2 in this case) and the “number of cells” refers to the cells present in the well when the infection is performed. Then, we obtained the mean value (second column) and calculated the copy number per microliter of genome purification (third column, in this case 5 in 500 dilution). Next, to obtain the VGE for each experiment we followed the formula: vge = copy number/ul (considering dilution factor)* ul of prep used / (number of cells *10).
2.7. Synthetic HPV Particle Production and Labeling
HEK293TT cells were used to create HPV-2,5,16,18,31, PMU, BPV, and PsVs with a luciferase reporter plasmid (pGL3 Luci), as previously reported [
55]. Using solutions of bovine serum albumin (BSA) at various concentrations as standards, the purity and capsid protein content of the PsV preparations were assessed using SDS-PAGE and Coomassie brilliant blue staining. Using a standard curve of the published plasmid DNA, the encapsidated DNA was subjected to real-time PCR analysis, allowing the copy number of each PsVs preparation to be determined. The HPV-16 L1-only VLPs were produced by employing a PsV-like procedure.
2.8. Cell Viability Assay
A Trypan Blue exclusion assay was used to count the number of live cells. After being seeded in a 6-well plate, the HaCaT cells were treated with either a control siRNA (siCTRL) or siRNA that selectively targeted COPB1 for 48 hours. The cells were trypsinized, collected, and stained with 0.2% trypan blue after 48 hours. Using a typical hemocytometer, cells were counted directly under a microscope to measure the quantity of viable cells. A one-tailed t test was used to assess statistical significance after normalising the number of viable cells in siCOPB1-treated samples to the number of viable cells obtained for siCTRL-treated samples in each experiment. Western blotting of the cellular lysates was used to verify siCOPB1’s effectiveness.
4. Discussion
From previous studies we know that HPV entry is a multistep process, and that HPV uses the cellular machinery for entry and trafficking to the cell nucleus [
56]. Using a proteomic screen we have found several interesting cargo proteins as potential interacting partners of HPV-16 L2. One of these cargo proteins was COPB1, which is a COPI protein subunit that is responsible for cargo binding and has the function of transporting cargo from TGN to ER [
57]. After confirming the binding between COPBI and HPV-11, -16, and -18 L2 proteins, through GST pulldown and co-immunoprecipitation assays, we investigated whether COPI has any role in the infectious entry process of HPV. Upon knockdown of COPB1 we found that there is more than 80% reduction in HPV-16 PsVs infectivity. We also investigated the effect of COPB1 knockdown on the infectivity of both high-risk and low-risk HPV types, and found a similar requirement for COPB1 for infectivity of diverse HPV types. Finally we also analysed the effects on infectivity of other papillomaviruses - Bovine papillomavirus (BPV) and Murine papillomavirus (MmuPV) - and found that COPB1 knockdown reduces BPV and MmuPV infectivity, Interestingly, we also confirmed that MmuPV L2 can bind COPB1 in GST pulldown assays. As we can see in
Figure 5 that some of HPVs are more sensitive upon COPB1 knockdown and some are less sensitive, which could be due to binding strength, as from GST pulldown in
Figure 2, we can see that COPB1 interacts strong with HPV-16 L2 as compared to HPV-18 L2 and the same goes with infectivity (HPV-16 is more sensitive upon COPB1 knockdown as compared to HPV-18). Thus, from this study, we have understood that COPB1 is a common cargo from TGN towards cell Nucleus for many HPV types, as it is localized in specific parts of the cell that are traversed by most viruses
en route to the cell nucleus through this path. This suggests that COPB1 could be a good target for blocking virus entry into the cell nucleus. Recently it has been also found that COPBI interacts with HPV-16 L2 through the R302/5A L2 [
58], It would also be important to check the infectivity of newly synthesized viruses from nucleus upon COPII protein knockdown and to check their interaction with virus cargo proteins.