1. Introduction
Nearly 200 million people are infected with HCV, and many of these people will develop severe liver disease, with or without hepatocellular carcinoma, because of persistent HCV infection. At present, HCV treatments include pegylated interferon α along with ribavirin and possibly a protease inhibitor. However, this treatment has problems with toxicity, cost, delivery, and efficacy. If severe liver disease develops and treatment fails, the only option is a liver transplant. Following liver transplant, an immunosuppressant is generally administered to help prevent graft rejection. In the majority of these cases, patients receive a calcineurin inhibitor, either cyclosporine (CsA) or tacrolimus. CsA and its nonimmunosuppressive analogs, DEBIO-25 (Alisporivir), SCY-635, and NIM811, are inhibitors of a class of cellular prolyl-peptidyl isomerases called cyclophilins (Cyps), and cyclophilin A is critical for HCV replication [
1,
2,
3,
4].
HCV has a large amount of genetic diversity and is classified in six different genotypes. The virus expresses at least 10 different proteins and several different HCV proteins have been shown to interact with cyclophilins (CypA, CypB, and possibly other cyclophilins) and modulate HCV replication [
2,
5,
6,
7,
8,
9]. The nonstructural protein 5A (NS5A) is a multifunctional, genetically diverse protein with three distinct hypothesized domains (DI, II and III) [
10]. Whether naturally occurring variation in NS5A alters cyclophilin susceptibility has been disputed [
11], but this variation could explain why an antiviral effect of cyclosporine in transplant patients may be detectable only in selected cases, if at all.
Earlier, we and others reported mutations within the DII region of NS5A which conferred resistance to CsA treatments [
2,
9,
12,
13,
14]. In this study we observed varying levels of CsA susceptibility for replicons containing the carboxy terminal regions of NS5A variants of similar lengths from different HCV genotypes. The replicon containing the H77 genotype 1a sequence exhibited the least susceptibility, while the replicons containing genotype 4a and 2a chimera, as well as the genotype 1b replicon, were all more susceptible. The CsA susceptibility of these C-tail NS5A variants correlated with a CypA binding assay. Data from replicons with mutations at amino acid 310, as well as CypA binding of a peptide which contains the 310 proline, but not the defined CypA binding target PDYN, still binds CypA. Our findings are consistent with CypA acting on multiple prolines outside of the previously identified CypA binding site and suggests that conformations of NS5A adopted in viral replication can vary in cyclophilin inhibitor susceptibility.
2. Results and Discussion
Several HCV subgenomic replicons have been described. However, replicons vary in the robustness of replication, and obviously vary genetically at several different places in the genome. To avoid the lower replication capacity associated with certain HCV genotypes and isolate the NS5A variation between genotypes, we constructed NS5A chimeric replicons in the backbone of the commonly used genotype 1b Con1 replicon which displays robust replication capacity. The
in vitro transcribed RNAs derived from the Con1bLN-wt (wild type), and chimeric replicons containing other NS5A genotypic sequences from amino acid 312 to the NS5A-NS5B cleavage site (Con1bLN-5A1a, Con1bLN-5A2a and Con1bLN-5A4a), were electroporated into Huh7.5 cells and luciferase activity was monitored over a period of five days in the presence and absence of CsA. As shown in
Figure 1A, all the replicons exhibited similar replication kinetics in the absence of CsA, thus indicating that the replaced polypeptide derived from genotypes 1a, 2a and 4a did not have deleterious effects on viral replication (red, blue, black and green lines). However the same replicons displayed contrasting susceptibility upon CsA treatment. The Con1bLN-5A4a replicon was found to be most susceptible (solid green
vs. dotted green lines, almost 100-fold less replication,
Figure 1B) to CsA treatment among all replicons. Although the Con1bLN-5A1a replicon had slightly lower replication capacity than the Con1bLN-wt replicon, the Con1bLN-5A1a replicon displayed the least susceptibility to CsA treatment (solid black
vs. dotted black lines, only 10-fold less replication compared to no CsA treatment). The Con1LN-wt and Con1LN-5A2a replicons had slightly better replication capacity than the Con1LN-5A1a and Con1LN-5A4a replicons in the absence of CsA, and showed less inhibition to CsA treatment compared to Con1LN-5A1a replicon (red and blue lines).
By making NS5A chimeras, we directly compared the cyclosporine susceptibility of specific NS5A sequences without the confounding effects of other parts of the genome. Due to the diversity of each subtype, our results do not imply that every genotype 1a HCV is less susceptible than every genotype 1b, only that there is a difference between H77 1a and Con1 in this carboxy terminal region of NS5A. Previous studies indicate the NS5A derived from different genotypes and the strong conservation of the PDYN binding site for CypA identified by NMR has been used to argue that cyclophilin inhibitors are “pangenotypic” and the heterogeneity of NS5A does not correlate with cyclophilin inhibition [
11]. Our data suggests most, if not all, HCV are susceptible to cyclophilin inhibitors. To our knowledge, this is the first demonstration that genotype 4 HCV is likely CsA susceptible, but our data suggests that NS5A polymorphisms outside the conserved PDYN sequence previously associated with CypA binding can also influence the degree of CsA susceptibility.
To test the interaction between NS5A genotypic carboxy terminal regions and CypA, we expressed NS5A polypeptides derived from different genotypes in a cell-free translation system in the presence of
35S cysteine/methionine and performed the CypA binding assay as previously described [
12]. We observed the polypeptide derived from genotype 1a bound more efficiently than the corresponding polypeptides of 1b, 2a and 4a genotypes (
Figure 1C). The polypeptides derived from genotype 1b and 2a bound to CypA but with less efficiency compared to genotype 1a (input lanes; 4 and 7, pull-down lanes; 6 and 9). The apparent migration of protein derived from 2a genotype of similar length in SDS-PAGE was slower compared to 1b genotype due to the presence of an additional 23 amino acids. The genotype 4a polypeptide displayed the least CypA binding in similar experimental set up (input lane 10 with pull-down lane 12). In general, these CypA binding patterns correlated well with their respective replicons’ susceptibility towards CsA (
Figure 1A). In all the pull-down assays, an active site mutant protein CypA55/60 was used as negative control. Overall, the CypA binding data indicates that the carboxy terminal regions of NS5A both interact with CypA as expected and correlate to some degree with the ability of the replicon to replicate in the presence of CsA.
Figure 1.
Role of HCV NS5A C-tails in CsA susceptibility and CypA binding. (A) The Con1bLN replicon was digested with XhoI and BstZ17I restriction enzymes (New England Biolabs) and a corresponding fragment from HCV genotype 1a genotype (aa 311–448; ARALPVWARP to TEDVVCC, accession no. AF009606) was cloned into the replicon, termed Con1bLN-5A1a. A similar strategy was used to clone genotype 2a fragment (aa 307–466; FRRPLPAWARP to EEDDTTVCC, accession no. AB047639) and genotype 4a (aa 313–449; RALPIWARPDYN to VSGSEDVVCC, accession no. Y11604.1), and termed Con1bLN-5A2a and Con1bLN-5A4a. The Huh7.5 cells were electroporated with in vitro synthesized RNA derived from Con1bLN-5A1a, Con1bLN-5A2a, Con1bLN-5A4a and Con1bLN-wt replicons. Equal numbers of electroporated cells were plated. The cells were either untreated (solid lines) or treated with CsA (dotted lines) for 120 hrs and luciferase activity was monitored every 24 hrs and presented. (B) The percent inhibition of respective replicons in (A) were calculated and presented. (C) The CypA binding capacity of NS5A regions derived from different genotypes. The 35S labeled proteins were incubated with either GST-CypA or GST-CypA55/60. The pulled down complexes were resolved by SDS-PAGE and signal was detected after autoradiography. The arrows indicate expected size of polyprotein (I = input (1/20th loaded); P = pull-down).
Figure 1.
Role of HCV NS5A C-tails in CsA susceptibility and CypA binding. (A) The Con1bLN replicon was digested with XhoI and BstZ17I restriction enzymes (New England Biolabs) and a corresponding fragment from HCV genotype 1a genotype (aa 311–448; ARALPVWARP to TEDVVCC, accession no. AF009606) was cloned into the replicon, termed Con1bLN-5A1a. A similar strategy was used to clone genotype 2a fragment (aa 307–466; FRRPLPAWARP to EEDDTTVCC, accession no. AB047639) and genotype 4a (aa 313–449; RALPIWARPDYN to VSGSEDVVCC, accession no. Y11604.1), and termed Con1bLN-5A2a and Con1bLN-5A4a. The Huh7.5 cells were electroporated with in vitro synthesized RNA derived from Con1bLN-5A1a, Con1bLN-5A2a, Con1bLN-5A4a and Con1bLN-wt replicons. Equal numbers of electroporated cells were plated. The cells were either untreated (solid lines) or treated with CsA (dotted lines) for 120 hrs and luciferase activity was monitored every 24 hrs and presented. (B) The percent inhibition of respective replicons in (A) were calculated and presented. (C) The CypA binding capacity of NS5A regions derived from different genotypes. The 35S labeled proteins were incubated with either GST-CypA or GST-CypA55/60. The pulled down complexes were resolved by SDS-PAGE and signal was detected after autoradiography. The arrows indicate expected size of polyprotein (I = input (1/20th loaded); P = pull-down).
We noticed a limited number of differences in the H77 1a genotype compared to the Con1b in the region, approximately 50 amino acids N terminal to WARPDYN. We therefore decided to explore the role of this region towards CsA susceptibility. We constructed two additional chimeric replicons encompassing amino acids 267–312, and the 267–448 region derived from genotype 1a in the backbone of Con1LN-wt replicon, and compared their susceptibility to that of Con1bLN-5A1a (here after known as Con1bLN-5A1a312-448). These replicons replicated well in the absence of CsA but had different degrees of susceptibility to CsA (
Figure 2A). The replicon containing the 1a region from 267–312 (purple lines) was the most sensitive to CsA (~3 log, ~99% reduction,
Figure 2B) suggesting that the 267–312 1a stretch conferred additional susceptibility to Con1bLN-wt replicon. While the Con1bLN-wt (dotted red) and the Con1bLN-5A267-448 replicon (dotted green lines) displayed similar CsA-like susceptibility to the Con1bLN-5A1a267-312 at 96 hours in
Figure 2B, notice the clear separation at an earlier time point in
Figure 2A (compare dotted purple to dotted red/green). This suggested to us the relatively increased susceptibility of Con1bLN-5A1a267-448 replicon, as compared to Con1bLN-5A1a312-448, which was due to the loss of the 267–3121b region. Since there were fewer consistent differences between 1b and 1a amino terminal to 312 compared to the carboxy terminal, we attempted to isolate a subtype “1a susceptibility” feature amino terminal to 312, rather than the “1a relative resistance” feature demonstrated in
Figure 1. Therefore, an effort was made to define the particular residue that alters CsA susceptibility in cell culture in the context of the replicon Con1bLN-5A1a312-448 that exhibits resistance. Shown in
Figure 3A is the amino acids sequence homology between genotype 1a and 1b amino acids N-terminal to an extremely high conserved region, WARPDYN, but within the NS5A region that contributes most to CypA binding as observed in our previous findings. We observed two distinct clusters that have noticeable amino acids changes, named Cluster1 (C1) and Cluster2 (C2) (
Figure 3A). These mutations were incorporated into Con1bLN-5A1a312-448 chimeric replicon termed Con1bLN-1aC1 and Con1bLN-1aC2 such that the genotype 1a residues at C1 and C2 replaced the residues present in genotype 1b. The Con1bLN-5A-1aC1 replicon (yellow lines) was competent enough for replication in both the presence and absence of CsA, thereby indicating that mutating amino acids EQ to DV did not affect genome replication following CsA treatment (
Figure 3B) and maintained similar CsA susceptibility compare to Con1bLN-5A1a312-448. On the other hand, when a 1a stretch of C2 was present in the Con1bLN-1aC2 replicon, there was increased susceptibility in the presence of CsA and produced ~3 log reduction, suggesting that a single or a combination of amino acids within C2 (RKPMI) are capable of altering CsA susceptibility. To our knowledge, this is the first observation that variation ~10 amino acids N-terminal to WARPDYN in genotype 1 alters cyclosporine susceptibility and is consistent with multiple prolines in this region being influenced by cyclophilin, as shown before [
13,
15]. Only mutations in the C2 cluster resulted in increased sensitivity to CsA, but not in cluster 1. Further analysis of amino acid residues present in C2 among genotype 1a reveals this stretch (RKSRRF
ARALPV) of amino acids are fairly conserved, except the alanine (
Figure 3A). The H77 1a genotype (AF009606) has alanine (bold underlined) at this particular position, compared to a proline in 1b. Logo analysis [
16] of amino acids in C2 cluster of 1a and 1b sequences demonstrates that while proline is well conserved in the genotype 1b lineage, genotype 1a commonly has either an alanine or a threonine (
Figure 4A).
To examine the role of amino acid 310 in genotype 1b (
Figure 4A) in its native genotype 1b context, we mutated this amino acid to either an alanine or threonine. The resulting replicons Con1bLN-P310A and Con1bLN-P310T were tested for CsA susceptibility as described above. Both the replicons carrying mutations Ala or Thr became more sensitive to CsA treatment than the Con1bLN-wt replicon indicating the role of proline at position 310 in CypA regulation in genotype 1b (
Figure 4B). We next tested if CypA could bind to this stretch of amino acids and, if so, whether or not a proline at 310 and/or 314 altered binding. A 15 amino acid long peptide representing this region was engineered as an N-terminal fusion protein with GFP and GST-CypA binding assay was performed as above. The GFP alone did not bind to either CypA55/60 or GST-CypA in a pull-down assay (
Figure 5A, lane 2 and 3). The peptide tagged GFP carrying 310P/314P amino acids bound well (lane 6), whereas a peptide tagged GFP containing 310A/314A amino acids exhibited little to no binding to CypA (lane 12), thus indicating that one or both of the prolines may contribute to CypA binding. We then expressed GFP tagged with peptides carrying mutations at the 310 and 314 positions to confirm which prolines contribute to CsA susceptibility. The peptide tagged GFP carrying mutation at 310P/314A (lane 9) and 310A/314P (lane 15) bound to CypA equally well as 310P/314P, thereby indicating that both the prolines contribute to CypA binding. Furthermore, mutating proline at position 310 to Ala or Thr did not abolish the CypA binding completely (lanes 15 and 18) partly due to the fact that both peptides 310A/314P and 310T/314P contain proline at 314 positions. However, a single point mutation at position 310 to either Ala or Thr rendered the 1b replicon more sensitive to CsA, indicating the proline at 310 has a role in CypA binding. Although we and others [
13] observed a contribution of 314P to CypA binding, our data indicates that 310 can also participate in CypA binding and contributes significantly directly or indirectly towards CsA susceptibility in tissue culture. The amino acid analysis of the C2 region also indicates that proline (P310) is highly conserved in non‑genotype 1a HCV. Interestingly, this proline residue, along with two other prolines (P310 and P315) in genotype 2a, have been found to be in the direct vicinity of NS5A:CypA interaction region when tested by gel filtration, circular dichroism, and NMR spectroscopy [
15]. Our transient replication data demonstrate mutating residues in the C2 region, and, in particular, genotype 1b P310 to alanine or threonine residue led to increased susceptibility of 1b replicon (pink and green lines), thus indicating the region’s critical involvement in CsA susceptibility. Our data is consistent with the genotype 2a data on P310 (homologous to P314 in genotype 1b), as well as P342 of genotype 2a in CsA regulation [
13], indicating critical roles of amino acids N-terminal as well as C-terminal to the DYN motif. The role of residue P306 in genotype 2a (corresponding to 310 in genotype 1b identified in this study) in the cyclophilin inhibitor Alisporivir susceptibility to our knowledge has not been investigated (
Figure 5B).
Figure 2.
Con1LN-wt replicons containing different lengths of the carboxy terminus of NS5A genotype 1a have different susceptibility to CsA. (
A) The CsA susceptibility replication assay of Con1bLN-5A1a-267-312, Con1bLN-5A1a-267-448 and Con1bLN-5A1a-312-448 (referring to same construct in
Figure 1A; Con1LN-5A1a) replicons was performed as described above. (
B) The percent inhibition of respective replicons in (A) were calculated and presented.
Figure 2.
Con1LN-wt replicons containing different lengths of the carboxy terminus of NS5A genotype 1a have different susceptibility to CsA. (
A) The CsA susceptibility replication assay of Con1bLN-5A1a-267-312, Con1bLN-5A1a-267-448 and Con1bLN-5A1a-312-448 (referring to same construct in
Figure 1A; Con1LN-5A1a) replicons was performed as described above. (
B) The percent inhibition of respective replicons in (A) were calculated and presented.
Figure 3.
Mutational analysis of Con1bLN-5A1a chimericreplicons reveal linear NS5A regions that alter CsA susceptibility. (A) The HCV 1b (358 isolates) and 1a (224 isolates) amino acid sequences were retrieved from the European HCV Database and subjected to amino acid homology analysis using web based program . The amino acid alignment analysis was performed on 53 amino acids N-terminal to the DYN region of HCV NS5A from 1b and 1a genotypes. The highly conserved DYN region is underlined. Site-directed mutagenesis was performed in two cluster regions, C1 and C2, in Con1bLN-wt replicon to make it similar to genotype 1a. The boxed residues were substituted in 1b replicon with genotype 1a amino acids. Amino acids of interest are numbered. The amino acids of interest are numbered. (B) The CsA susceptibility replication assay was performed on replicons Con1bLN-1aC1, Con1bLN-1aC2, Con1bLN-5A1a and Con1bLN-wt replicons. (C) The percent inhibition of respective replicons in (B) were calculated and presented.
Figure 3.
Mutational analysis of Con1bLN-5A1a chimericreplicons reveal linear NS5A regions that alter CsA susceptibility. (A) The HCV 1b (358 isolates) and 1a (224 isolates) amino acid sequences were retrieved from the European HCV Database and subjected to amino acid homology analysis using web based program . The amino acid alignment analysis was performed on 53 amino acids N-terminal to the DYN region of HCV NS5A from 1b and 1a genotypes. The highly conserved DYN region is underlined. Site-directed mutagenesis was performed in two cluster regions, C1 and C2, in Con1bLN-wt replicon to make it similar to genotype 1a. The boxed residues were substituted in 1b replicon with genotype 1a amino acids. Amino acids of interest are numbered. The amino acids of interest are numbered. (B) The CsA susceptibility replication assay was performed on replicons Con1bLN-1aC1, Con1bLN-1aC2, Con1bLN-5A1a and Con1bLN-wt replicons. (C) The percent inhibition of respective replicons in (B) were calculated and presented.
Figure 4.
Mutational analysis of Con1bLN-wt replicon at position 310 and analysis of CsA susceptibility. (A) Logo analysis of the C2 region of 358 isolates from genotype 1b and 224 from genotype 1a of HCV NS5A. The boxed proline residue at position 310 appears to be highly conserved among most HCV 1b viruses. Sequences derived from genotype 1b and 1a are marked. (B) The CsA susceptibility of Con1bLN-P310A, Con1bLN-P310T was compared to Con1bLN-wt replicons. (C) The percent inhibition of respective replicons in (B) were calculated and presented.
Figure 4.
Mutational analysis of Con1bLN-wt replicon at position 310 and analysis of CsA susceptibility. (A) Logo analysis of the C2 region of 358 isolates from genotype 1b and 224 from genotype 1a of HCV NS5A. The boxed proline residue at position 310 appears to be highly conserved among most HCV 1b viruses. Sequences derived from genotype 1b and 1a are marked. (B) The CsA susceptibility of Con1bLN-P310A, Con1bLN-P310T was compared to Con1bLN-wt replicons. (C) The percent inhibition of respective replicons in (B) were calculated and presented.
Figure 5.
CypA binding analysis of a peptide containing Ala, Pro and Thr at position 310. (
A)
35S labeled proteins derived from NS5A peptide fused to the GFP coding sequences were incubated with either GST-CypA55/60 or with GST-CypA. The pulled-down complexes were resolved by SDS-PAGE and signal was detected by autoradiography. The arrow indicates expected size proteins (I = input (1/20th loaded); P = pull-down). (
B) Alignment CypA binding region from genotype 1b and 2a. The highly conserved DYN region is underlined and amino acids of interests are numbered. In addition to P310 described in this study, the boxed residues were demonstrated previously to regulate the Alisporivir susceptibility in genotype 2a [
13].
Figure 5.
CypA binding analysis of a peptide containing Ala, Pro and Thr at position 310. (
A)
35S labeled proteins derived from NS5A peptide fused to the GFP coding sequences were incubated with either GST-CypA55/60 or with GST-CypA. The pulled-down complexes were resolved by SDS-PAGE and signal was detected by autoradiography. The arrow indicates expected size proteins (I = input (1/20th loaded); P = pull-down). (
B) Alignment CypA binding region from genotype 1b and 2a. The highly conserved DYN region is underlined and amino acids of interests are numbered. In addition to P310 described in this study, the boxed residues were demonstrated previously to regulate the Alisporivir susceptibility in genotype 2a [
13].