- Open Access
Envelope 2 protein phosphorylation sites S75 & 277 of hepatitis C virus genotype 1a and interferon resistance: A sequence alignment approach
Virology Journalvolume 8, Article number: 71 (2011)
Hepatitis C is a major health problem affecting more than 200 million individuals in world including Pakistan. Current treatment regimen consisting of interferon alpha and ribavirin does not always succeed to eliminate virus completely from the patient's body.
Interferon induced antiviral protein kinase R (PKR) has a role in the hepatitis C virus (HCV) treatment as dsRNA activated PKR has the capacity to phosphorylate the serine and threonine of E2 protein and dimerization viral RNA. E2 gene of hepatitis C virus (HCV) genotype 1 has an active role in IFN resistance. E2 protein inhibits and terminates the kinase activity of PKR by blocking it in protein synthesis and cell growth. This brings forward a possible relation of E2 and PKR through a mechanism via which HCV evades the antiviral effect of IFN.
A hybrid in-silico and wet laboratory approach of motif prediction, evolutionary and structural anlysis has pointed out serine 75 and 277 of the HCV E2 gene as a promising candidate for the serine phosphorylation. It is proposed that serine phosphorylation of HCV E2 gene has a significant role in interferon resistance.
Hepatitis C virus is a key universal health issue  affecting approximately 200 million individuals worldwide and over 4 million in the United States alone, where it is the most common blood-borne infection . It is the main reason of persistent liver infection and the most familiar sign for liver transplantation . In 60-85% cases HCV develops to cirrhosis and hepatocellular carcinoma . Presently, in Pakistani population 17 million people are infected with HCV and 8-10% individuals are HCV carriers . HCV is a member of genus, hepacivirus, family Flaviviridae and is a positive sense single stranded RNA virus . HCV genome is about 9.6 kb in length . The large open reading frame of Viral RNA having 5` and 3` untranslated regions which is translated into a single polypeptide of 3010 to 3033 amino acids. Which is processed by host as well as viral proteases to yield 10 mature individual proteins out of which 3 are structural and 7 are nonstructural [8, 9]. Detailed structure of HCV virus is still unclear. However, the infectious viral particles are composed of lipid envelope glycoproteins E1 and E2 .
In spite of much recent advancements, still no vaccine is available against HCV infection. The current therapy for HCV infection is pegylated interferon alpha separately or in combination with ribavirin  but it eradicates the virus in only 50-80% of cases and has serious side effects . IFN system is the first line of protection against viral infection in mammals . Many viruses have developed mechanisms to dodge the IFN-dependent cellular response . Amongst these HCV is a significant example in which (70-80%), cases runs away the host defenses and develops a chronic infection. Pathological outcomes of HCV infection changes from individual to individual, unstable from asymptomatic state to liver fibrosis, steatosis, finally to hepatocellular carcinoma [15, 16]. The factors upon which the success or failure of the antiviral therapy depends are unspoken yet, and their recognition characterizes a main confront in HCV virology . IFN-interacts with cells and modifies the expression of a number of genes [17, 18]. E2 glycoprotein is the first viral factor that gets in touch with the host cell surface receptors thus it has an important role in vaccine designing and drug objective [19, 20]. Primary objective of HCV vaccine is to initiate potent humoral responses against E2 protein .
After the death of infected cells, virus particles are discharged and they can infect the cells present in close proximity. Interferon is released from the infected cells to inform the neighboring cells about the presence of virus. In its rapid response, PKR is produced by these adjoining cells. HCV resistance to IFN-treatment is partially related to inhibition of interferon induced anti-viral protein PKR. Interferon induce many protective mechanisms in cells and amongst these the major role in cell protection from many viruses is illustrated by double-stranded RNA (dsRNA)-activated protein kinase PKR [22, 23].
Two HCV proteins (NS5A and E2) are involved in IFN resistance through inhibition of the IFN-α induced double stranded-RNA (dsRNA)-activated protein kinase (PKR) [24, 25]. PKR is a kinase enzyme that reveals varied activities. PKR exhibits autophosphorylation of many serine and threonine positions and dimerization of dsRNA. It also phosphorylates the translation initiation factor eIF-2 (α subunit) that directs towards blockage of protein synthesis . Because of these properties, PKR is taken as an arbitrator of antiviral and anti-inflammatory role of IFN-α . E2 protein blocks PKR activation to bypass its function [24–28]. Taylor et al., (1999)  reported that HCV E2 protein encloses a 12 amino acid sequence domain that is highly homologous to autophosphorylation positions of PKR and initiation factor eIF2. PKR and eIF2 form a phosphorylation homology domain (PePHD). E2 protein inhibits kinase activity of PKR and terminates its blocking role in protein synthesis and cell growth. This suggests that the relationship of E2 and PKR can be considered as a major mechanism by which HCV evades the antiviral effect of IFN.
HCV RNA Detection and identification of E2 gene
Serum samples were collected at the Division of Molecular Virology, National Centre of Excellence in Molecular Biology, University of the Punjab, Lahore. HCV 1a genotype serum samples were used in the current study. 100 μl serum was used to extract viral RNA by using RNA Isolation Kit (Sacace, Italy). Primer 3 software was used to designed primer for the amplification of E2 gene. At 5' end a start codon was introduced artificially in the primers for expression of gene in mammalian cell culture. The required amplified product was of approximately 1089 bps which was confirmed on ethidium bromide stained 1.2% agarose gel evaluated under UV transilluminator and photographed. The necessary band was cut out from the gel and the DNA was extracted with DNA isolation kit (Fermentas Inc. Germany). Isolated DNA was suspended in depc treated water and used for further studies.
Construction of expression vector
Mammalian expression vector pcDNA 3.1 (Invitrogen Tech USA) was used to clone the amplified product encoding E2 gene between Hindi III and EcoR I sites. The plasmid was transformed in chemically-competent cells. The entire ligation reaction was added to a 100 μl aliquot of TOP10F cells. Then 500 μl of SOC medium was added to the cells and incubated at 37°C for 1 hour. Competent cells having plasmid were then selected by spreading the culture onto a Luria-Bertani (LB) agar plate containing 100 μg/ml ampicillin and 12.4 ug/ml tetracyclin. Colonies selection was made by incubating the plate overnight at 37°C. To identify bacteria harboring cloned E2 gene, individual colonies were used to directly inoculate PCR reactions. The PCR reactions were prepared with 5 pmol each of vector-specific primers i.e. T7 (TAATACGACTCACTATAGGG) and BGH (TAGAAGGCACAGTCGAGG). Each reaction was prepared to a final volume of 20 μl. Following amplification, the amplified product was checked on a 1.2%, ethidium bromide stained agarose gel, a successful cloning reaction being visualized as a product at approximately 1.3 kb. Colonies identified as possessing a desired clone were then used to inoculate 3 ml LB culture containing 100 μg/ml ampicillin, shaking at 225 rpm overnight at 37°C. Plasmid was purified using a Plasmid Miniprep Kit, (Fermentas Life Sciences technologies, USA) according to manufacturer's protocol. Quantification of the plasmid prep. was performed by using spectrophotometer. Successful cloning was confirmed through PCR, restriction digestion of plasmid (pcDNA 3.1/myc vector) with EcoR1 & Hindi III and sequencing reaction. At least three clones were sequenced in both directions according to the sequencing protocol.
Applied biosystems prism dye termination method was used to sequence the DNA (PCR products and plasmids) with specific sense and anti sense primers. The sequences reported in this paper have been deposited in genbank database (Accession no. GU736411). In a 10 μl, of total reaction volume, about 300 ng of plasmid, 10 pmol of primer, 1 μl of big dye and 1 μl of dilution buffer were used for each sequencing reaction. Sequencing was achieved in a thermal cycler with parameters; 94°C for 20 s, 58°C for 20 s, 60°C for 4 min, repeated 25 times. The pellet was air dried and DNA was analyzed using an ABI prism sequencer.
Identification of Phosphorylation sites
To study the possible outcome of predication with help of NetPhos 2.0 software as described by Blom et al. (1999) . This software checks the possible post translation modification in the local 1a local strain primary sequence in the amino acid. As a result of analysis the target motif for different kinases is visualized [29, 30]. Scansite is kept at low stringency in order to determine the maximum number of sites that may participate in phosphorylation, upon which further predication is done.
Protein Structure Analysis
An ab-initio model was designed by using I-TASSER as no template model was available in protein Data Bank . Data was uploaded and models were generated. A 3D structure of predicted phosphorylation sites was generated by using two servers, Chimera  and SWISS PDB viewer  were used to obtain 3D structure. It was checked that potential phosphorylation sites are surface exposed or hidden inside. Five sites were observed at the exposed surface of model (aa: S75, S95, S118, S277, Y211).
Protein sequences of envelope gene for genotype 1a from different countries of the world (Japan, France, USA, and UK) were obtained from NCBI. All sequences (GQ898898, EU482831, EU234064, EU362889, DQ061315, AY958052, AY958057, AY956468, AB520610, and AF529293) were then aligned with local envelope gene by using CLUSTALW . A neighbor joining tree was generated using PHYLIP (Figure 1) .
To evaluate its function in disease succession, amplified E2 gene was cloned in mammalian expression vector pc DNA 3.1. A CMV promoter is present in pc DNA 3.1 vector that can transduce eukaryotic cells very efficiently for transient and stable expression studies. Clone was sequenced by the applied biosystems prism dye termination method using vector specific and gene specific primers (Accession no. GU736411). Documented sequences for 1a genotype from diversified areas of the world were compared with local envelope gene sequence to get the percentage nucleotide identity (PNI) (Table 1) and CLUSTAL W. Sequence dissimilarities in the envelope genes explained the different levels of disease outcomes in patients infected with same genotype. These nucleotide differences can help out to plan the genotype specific therapy to avoid and minimize HCV infections.
Possible phosphorylation sites were calculated in the cytoplasmic domain of the E2 protein that can be implicated in the interferon resistance. Twelve putative sites (Figure 2) were chosen as potential phosphorylation sites (S-6, Thr-4, Tyr-2). Stringency of Phosphorylation site predictors was increased as they have tendency to over-predict. Only those motifs were selected that showed a NetPhos score of 0.8 or greater. After it they were also analyzed by Scansite. Finally six sites were predicted by both servers (Table 2).
An online server NetSurfP was used to find the surface accessibility of local envelope sequence, as the phosphorylation sites should be exposed on the surface of proteins (Table 2). Discovery Studio and SWISS PDB Viewer were applied to visualize 3D protein structure of the putative sites. PDB structure was built up by using the server I-TASEER (Ab-initio protein structure predictor) . After this analysis two phosphorylation sites (S75 and S277) were found to be most reliable sites (Figure 3). Scansite was used for finding the phosphorylation interaction motifs and found that S75 and S277 interact with the CLK2 Kinase and YWHAZ which were then investigated in GeneCards and was confirmed from UniGene and UniProt (Table 3).
IFN usually causes the transcription of some antiviral genes. However, detailed principal mechanisms exploring high IFN-α resistance in HCV genotype 1a infected patients are still ambiguous. A double stranded RNA activated PKR phosphorylates the translation initiation factor (eLF2) and thus inhibits its protein synthesis ability. A straight relationship was found in a study of E2 protein of HCV genotype 1a and IFN inhibition pathway. It was shown that it inhibits PKR through a 12 amino acid sequence E2-PePHD domain, and for PKR this domain can perform the function of pseudo-substrate. The functions of PKR i.e. inhibition of viral protein synthesis and its kinase activity, are blocked by this domain. PKR autophosphorylation sites are highly identical to the PePHD domain in HCV genotype 1a than other genotypes .
It has been demonstrated that E2 protein of HCV genotype 1 shares homology with PKR and eIF2 phosphorylation sites and therefore inhibits PKR by binding to it in vitro in mammalian and yeast cells. HCV therefore has evolved a mechanism in the form of E2-PKR interaction to block IFN activity. The possible outcome of PKR inhibition is not only IFN resistance and infection persistence but also cell growth promotion which ultimately leads to hepatocellular carcinoma (HCC).
Since E2 PePHD domain of HCV genotype 3a share less homology to PKR and eIF2 phosphorylation sites than do genotypes 1a and 1b. Therefore, it can be the most satisfactory reason that patients with 3a genotype react to interferon treatment more proficiently as compared to genotype 1a that shows maximum resistance to IFN therapy [37, 38]. The relationship of E2-PePHD domain and its behavior during interferon therapy is quite indistinct. Some studies also explained the modifications happening in the PePHD domain during treatment. This work reveals that in non-responders, slight variations present in pretreated patients quickly evades during treatment . HCV has developed ways to undo the antiviral effects of PKR [40, 41] and it is assumed that IFN resistance may be due to phosphorylation of envelope proteins particularly of E2 protein. Our study focused on detecting such sites in our local isolate of HCV genotype 1a and confirmed two phosphorylation sites (S75 and S277) at the surface of E2 protein that interacted with CLK2 Kinase and YWHAZ. It is further proposed that this interaction might be responsible for HCV resistance to antiviral effects of IFN which could be confirmed by dephosphorylating these sites and analyzing its effects on PKR binding and inhibition.
The role of envelope protein 2 (E2) of hepatitis C virus (HCV) is found to interact with double stranded RNA-dependent protein kinase (PKR) in previous studies. In the present study, the role of phosphorylation of E2 protein towards interferon sensitivity & PKR response was determined. HCV E2 protein was found to contain a sequence identical to phosphorylation sites of the interferon-inducible protein kinase (PKR). As a result E2 inhibits the kinase activity of PKR. This work will help to identify factors that can favor a successful innate immune response to HCV infection.
Khan S, Rai MA, Khan A, Farooqui A, Kazmi SU, Ali SH: Prevalence of HCV and HIV infections in 2005-Earthquake-affected areas of Pakistan. BMC Infect Dis 2008, 8: 147. 10.1186/1471-2334-8-147
Alter MJ: Epidemiology of hepatitis C virus infection. World J Gastroenterol 2007, 13: 2436-2441.
Brown RS: Hepatitis C and liver transplantation. Nature 2005, 436: 973-978. 10.1038/nature04083
Shepard CW, Finelli L, Alter MJ: Global epidemiology of hepatitis C virus infection. Lancet Infect Dis 2005, 5: 558-567. 10.1016/S1473-3099(05)70216-4
Idrees M, Lal A, Naseem M, Khalid M: High prevalence of hepatitis C virus infection in the largest province of Pakistan. J Dig Dis 2008, 9: 95-103. 10.1111/j.1751-2980.2008.00329.x
Bartenschlager R, Lohmann V: Replication of hepatitis C virus. J Gen Virol 2000, 81: 1631-1648.
Jhaveri R, Kundu P, Shapiro AM, Venkatesan A, Dasgupta A: Effect of Heptitis C Virus Core Protein on Cellular Gene Expression: Specific Inhibition of Cyclooxygenase 2. The Journal of Infectious Diseases 2005, 191: 1498-506. 10.1086/429301
Brass V, Moradpour D, Blum HE: Molecular virology of hepatitis C virus (HCV): 2006 update. Int J Med Sci 2006, 3: 29-34.
Penin F, Dubuisson J, Rey FA, Moradpour D, Pawlotsky JM: Structural biology of hepatitis C virus. Hepatology 2004, 39: 5-19. 10.1002/hep.20032
Budkowska A: Mechanism of cell infection with hepatitis C virus (HCV)--a new paradigm in virus-cell interaction. Pol J Microbiol 2009, 58: 93-98.
Di Bisceglie AM, Hoofnagle JH: Optimal therapy of hepatitis C. Hepatology 2002, 36: S121-127. 10.1002/hep.1840360716
Keam SJ, Cvetkovic RS: Peginterferon-alpha-2a (40 kD) plus ribavirin: a review of its use in the management of chronic hepatitis C mono-infection. Drugs 2008, 68: 1273-1317. 10.2165/00003495-200868090-00006
Grandvaux N, tenOever BR, Servant MJ, Hiscott J: The interferon antiviral response: from viral invasion to evasion. Curr Opin Infect Dis 2002, 15: 259-267.
Katze MG, He Y, Gale M: Viruses and interferon: a fight for supremacy. Nat Rev Immunol 2002, 2: 675-687. 10.1038/nri888
Frese M, Pietschmann T, Moradpour D, Haller O, Bartenschlager R: Interferon-alpha inhibits hepatitis C virus subgenomic RNA replication by an MxA-independent pathway. J Gen Virol 2001, 82: 723-733.
Lauer GM, Walker BD: Hepatitis C virus infection. N Engl J Med 2001, 345: 41-52. 10.1056/NEJM200107053450107
Marcello T, Grakoui A, Barba-Spaeth G, Machlin ES, Kotenko SV, MacDonald MR, Rice CM: Interferons alpha and lambda inhibit hepatitis C virus replication with distinct signal transduction and gene regulation kinetics. Gastroenterology 2006, 131: 1887-1898. 10.1053/j.gastro.2006.09.052
Hayashi J, Stoyanova R, Seeger C: The transcriptome of HCV replicon expressing cell lines in the presence of alpha interferon. Virology 2005, 335: 264-275. 10.1016/j.virol.2005.02.018
Cocquerel L, Kuo CC, Dubuisson J, Levy S: CD81-dependent binding of hepatitis C virus E1E2 heterodimers. J Virol 2003, 77: 10677-10683. 10.1128/JVI.77.19.10677-10683.2003
Cooper S, Erickson AL, Adams EJ, Kansopon J, Weiner AJ, Chien DY, Houghton M, Parham P, Walker CM: Analysis of a successful immune response against hepatitis C virus. Immunity 1999, 10: 439-449. 10.1016/S1074-7613(00)80044-8
Steinmann E, Brohm C, Kallis S, Bartenschlager R, Pietschmann T: Efficient trans-encapsidation of hepatitis C virus RNAs into infectious virus-like particles. J Virol 2008, 82: 7034-7046. 10.1128/JVI.00118-08
Kaufman RJ: The Double-stranded RNA-activated Protein Kinase PKR. In Translational Control of Gene Expression. Volume 39. 2nd edition. Cold Spring Harbor Laboratory Press; 2000:503-527.
Pe'ery T, Mathews MB: Viral translational strategies and host defense mechanisms. In Translational Control of Gene Expression. Volume 39. 2nd edition. Cold Spring Harbor Laboratory Press; 2000:371-424.
Gale M Jr, Korth MJ, Tang NM, Tan SL, Hopkins DA, Dever TE, Polyak SJ, Gretch DR, Katze MG: Evidence that hepatitis C virus resistance to interferon is mediated through repression of the PKR protein kinase by the nonstructural 5A protein. Virology 1997, 230: 217-227. 10.1006/viro.1997.8493
Taylor DR, Shi ST, Romano PR, Barber GN, Lai MM: Inhibition of the interferon-inducible protein kinase PKR by HCV E2 protein. Science 1999, 285: 107-110. 10.1126/science.285.5424.107
MacQuillan GC, Caterina P, de Boer B, Allan JE, Platten MA, Reed WD, Jeffrey GP: Ultra-structural localisation of hepatocellular PKR protein using immuno-gold labelling in chronic hepatitis C virus disease. J Mol Histol 2009, 40: 171-176. 10.1007/s10735-009-9227-0
Stark GR, Kerr IM, Williams BR, Silverman RH, Schreiber RD: How cells respond to interferons. Annu Rev Biochem 1998, 67: 227-264. 10.1146/annurev.biochem.67.1.227
Gale M Jr, Blakely CM, Kwieciszewski B, Tan SL, Dossett M, Tang NM, Korth MJ, Polyak SJ, Gretch DR, Katze MG: Control of PKR protein kinase by hepatitis C virus nonstructural 5A protein: molecular mechanisms of kinase regulation. Mol Cell Biol 1998, 18: 5208-5218.
Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol 1999, 294: 1351-1362. 10.1006/jmbi.1999.3310
Obenauer JC, Cantley LC, Yaffe MB: Scansite 2.0: Proteome-wide prediction of cell signaling interactions using short sequence motifs. Nucleic Acids Res 2003, 31: 3635-3641. 10.1093/nar/gkg584
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera--a visualization system for exploratory research and analysis. J Comput Chem 2004, 25: 1605-1612. 10.1002/jcc.20084
Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis 1997, 18: 2714-2723. 10.1002/elps.1150181505
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22: 4673-4680. 10.1093/nar/22.22.4673
Felsenstein J: Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol 1981, 17: 368-376. 10.1007/BF01734359
Zhang Y: I-TASSER server for protein 3D structure prediction. BMC Bioinformatics 2008, 9: 40. 10.1186/1471-2105-9-40
Gaudy C, Lambele M, Moreau A, Veillon P, Lunel F, Goudeau A: Mutations within the hepatitis C virus genotype 1b E2-PePHD domain do not correlate with treatment outcome. J Clin Microbiol 2005, 43: 750-754. 10.1128/JCM.43.2.750-754.2005
Trepo C LK, Niederau C, et al.: Pegylated interferon alfa-2b (PEG-INTRON) monotherapy is superior to interferon alfa-2b (INTRON A) for the treatment of chronic hepatitis C. J Hepatol 2000, 32: 29. 10.1016/S0168-8278(00)80446-8
Zeuzem S, Feinman SV, Rasenack J, Heathcote EJ, Lai MY, Gane E, O'Grady J, Reichen J, Diago M, Lin A, et al.: Peginterferon alfa-2a in patients with chronic hepatitis C. N Engl J Med 2000, 343: 1666-1672. 10.1056/NEJM200012073432301
Chayama K, Suzuki F, Tsubota A, Kobayashi M, Arase Y, Saitoh S, Suzuki Y, Murashima N, Ikeda K, Takahashi N, et al.: Association of amino acid sequence in the PKR-eIF2 phosphorylation homology domain and response to interferon therapy. Hepatology 2000, 32: 1138-1144. 10.1053/jhep.2000.19364
Guo JT, Bichko VV, Seeger C: Effect of alpha interferon on the hepatitis C virus replicon. J Virol 2001, 75: 8516-8523. 10.1128/JVI.75.18.8516-8523.2001
Mathews MB, Hershey JWB, Sonenberg N: Translational Control. NY.: Cold Spring Harbor Laboratory Press, Cold Spring Harbor; 1996.
The authors declare that they have no competing interests.
MI conceived the study and critically reviewed the manuscript. SA performed, nucleotide sequencing and analyzed the results. SA, MA, AH and MI drafted the manuscript. MA, MA, MI, SB, SS, IR, LA and MS participated in data analysis. All the authors studied and approved the final manuscript.