A fusion peptide in preS1 and the human protein disulfide isomerase ERp57 are involved in hepatitis B virus membrane fusion process
Abstract
Cell entry of enveloped viruses relies on the fusion between the viral and plasma or endosomal membranes, through a mechanism that is triggered by a cellular signal. Here we used a combination of computational and experimental approaches to unravel the main determinants of hepatitis B virus (HBV) membrane fusion process. We discovered that ERp57 is a host factor critically involved in triggering HBV fusion and infection. Then, through modeling approaches, we uncovered a putative allosteric cross-strand disulfide (CSD) bond in the HBV S glycoprotein and we demonstrate that its stabilization could prevent membrane fusion. Finally, we identified and characterized a potential fusion peptide in the preS1 domain of the HBV L glycoprotein. These results underscore a membrane fusion mechanism that could be triggered by ERp57, allowing a thiol/disulfide exchange reaction to occur and regulate isomerization of a critical CSD, which ultimately leads to the exposition of the fusion peptide.
Introduction
Hepatitis B is a major public health problem; it affects over 250 million people worldwide, and 850,000 deaths occur each year as a result of hepatitis B complications. The structure of its etiological agent, the hepatitis B virus (HBV), features a nucleocapsid that is surrounded by a lipid bilayer containing the HBV surface antigen (HBsAg) consiting in three envelope glycoproteins (GPs) designated as small (S), medium (M), and large (L), which are the products of a single open reading frame. They share the C-terminal S domain that contains four putative transmembrane (TM) domains. The L and M proteins have N-terminal extensions (preS1/preS2 and preS2, respectively) that mediate diverse functions in nucleocapsid binding and receptor recognition (Baumert et al., 2014). The first 2- to 75-amino-acid sequence of the preS1 domain of the L protein (Blanchet and Sureau, 2007; Bremer et al., 2011; Le Seyec et al., 1999) and the antigenic loop (AGL) of the S domain (Le Duff et al., 2009; Salisse and Sureau, 2009; Schulze et al., 2007) have been identified as essential determinants for infectivity of HBV and hepatitis delta virus (HDV), a pathogen that depends on HBV GPs for its propagation.
Entry of enveloped viruses into cells can be defined as the sequence of events occurring from the attachment of the virus to the host cell until the release of the genome into the cytoplasm, via fusion between viral and cellular membranes. Like for most enveloped viruses, HBV entry into cells is a finely regulated and complex process consisting of different steps, in which several viral and cellular factors are involved. Its first step involves low-affinity binding to heparan sulfate proteoglycans (HSPGs) residing on the hepatocytes’ surface (Leistner et al., 2008; Schulze et al., 2007). This attachment is mediated by the preS1 region of the L protein and/or the AGL of the S protein (Ni et al., 2014; Schulze et al., 2007). Afterwards, the virus interacts with its high-affinity receptor, the sodium taurocholate-cotransporting polypeptide (NTCP) (Ni et al., 2014; Yan et al., 2012), through the amino-terminal end of the L protein preS1 domain (Glebe et al., 2005; Gripon et al., 2005; Yan et al., 2012). NTCP is an integral membrane protein expressed at the basolateral membrane of hepatocytes, which explains the tropism of HBV for the liver.
The post-binding entry steps of HBV occur through endocytosis; however, the exact mechanism is still unclear and somehow controversial. One early study in HepaRG cells showed that HBV is internalized via caveolin-mediated endocytosis (Macovei et al., 2010). Nevertheless, inhibition of caveolin-mediated endocytosis or silencing of caveolin-1 did not impair HBV infection in tupaia hepatocytes (Bremer et al., 2009) or in HepaG2-NTCP cells (Herrscher et al., 2020). Contrastingly, several other studies presented evidence that HBV endocytosis is clathrin-dependent (Herrscher et al., 2020; Huang et al., 2012; Umetsu et al., 2018). Recent studies have reported that HBV infection of HepaRG cells depends on Rab5 and Rab7 (Macovei et al., 2013), which are GTPases involved in the biogenesis of endosomes, and that the epidermal growth factor receptor (EGFR) is a host-entry cofactor that interacts with NTCP and mediates HBV internalization (Iwamoto et al., 2019). These findings support the hypothesis that HBV is transported from early to mature endosomes. After the early endosome stage, translocation is associated with a gradually decreasing pH, from about 6.2 in early endosomes to close to 5.5 in late endosomes, which allows fusion of many enveloped viruses with the endosomal membrane. However, in the case of HBV, pharmacological agents that raise or neutralize the pH in the endocytic pathway do not affect infection (Macovei et al., 2010; Macovei et al., 2013; Rigg and Schaller, 1992). Furthermore, treatments with protease inhibitors have no effect on infection (Macovei et al., 2013), suggesting that HBV transport into the degradative branch of the endocytic pathway is not required per se to initiate this process.
Virus entry by membrane fusion involves interactions between viral fusion proteins and host receptors, which result in conformational changes of the virus envelope proteins. However, the molecular determinants and mechanism of membrane fusion of HBV remain to be defined. Previous results have indicated the essential role of the cysteine residues of the AGL, as shown by the reduction of virus entry levels by inhibitors of thiol/disulfide exchange reaction (Abou-Jaoudé and Sureau, 2007), hence suggesting a redox state responsible for conformational changes that can have a role during the fusion step.
Here, using a combination of computational and experimental approaches, we sought to better understand how HBV induces the fusion of its lipid membrane with that of the infected cell. Specifically, using a coevolution analysis of HBV GPs and molecular modeling combined with experimental investigations ex vivo in molecular virology and in vivo in liver humanized mice, we provide evidence that the mechanism triggering HBV membrane fusion involves ERp57, a cellular protein disulfide isomerase (PDI). Furthermore, our results highlight the role of specific cysteines in the AGL determinant as well as a sequence (aa 48–66) in the preS1 determinant that could ultimately act as a fusion peptide mediating HBV membrane fusion.
Results
HBV membrane fusion is independent of acidic pH and receptor expression
To investigate the fusion activation mechanism and to identify the fusion determinants of HBV, we designed a cell-cell fusion assay whereby Huh7 ‘donor’ cells, expressing a luciferase reporter gene under control of the HIV-1 promoter, were co-cultured with either Huh7-tat or Huh7-NTCP-tat ‘indicator’ cells, expressing the HIV-1 transactivator of transcription (Tat) protein, which induces luciferase expression only in fused donor and indicator cells (Lavillette et al., 2007). We transfected donor cells with pT7HB2.7 (Sureau et al., 1994), an expression plasmid encoding the wild-type HBV glycoproteins L, M, and S. The transfected donor cells were then co-cultivated with Huh7-tat or Huh7-NTCP-tat indicator cells for 1 day. The medium of the co-cultures was then acidified at pH 4 for 3 min to trigger fusion and the next day, the luciferase activity in the lysates of co-cultured cells was measured as a read-out of membrane fusion (Figure 1A). The GPs of vesicular stomatitis virus (VSV) or of Crimean-Congo hemorrhagic fever virus (CCHFV) were used as controls for viruses that need acidic pH to promote membrane fusion. We found that HBV GPs induced similar levels of fusion in co-cultures that were exposed to either acidic or neutral pH, as well as in co-cultures lacking or expressing the NTCP receptor (Figure 1A; see raw data in Figure 1—figure supplement 1). Since HBV entry requires HSPGs to mediate the capture of its viral particles through HBsAg (Leistner et al., 2008; Schulze et al., 2007), we addressed whether blocking of HBsAg/HSPG interaction could inhibit cell-cell fusion using heparin as the competitor. Yet, while the applied doses of heparin could prevent cell-free entry, as shown previously (Schulze et al., 2007), addition of soluble heparin to the co-cultures did not prevent HBsAg-mediated fusion, whether the indicator cells expressed NTCP or not (Figure 1B). We confirmed these results by using CHO and CHO-pgsB618 (Richard et al., 1995) cells as donor and/or indicator cells. While both cell types do not express NTCP, only the former expresses HSPGs. We found that cell-cell fusion could be detected for either indicator cell type to the same extent as for Huh7 cells (Figure 1C).
Altogether, these results indicated that cell-cell fusion mediated by HBV GPs is independent of acidic pH and requires neither HSPG nor NTCP receptor, which underscores an alternative fusion trigger.
The preS1 domain of HBV L protein harbors a critical determinant of membrane fusion
The L, M, and S GPs of HBV are produced by a single open reading frame and share a common C-terminal S domain. M and L proteins harbor additional N-terminal extensions (preS2 and preS1/preS2, respectively), with preS1 harboring the NTCP-binding determinant (Glebe et al., 2005; Gripon et al., 2005). Noteworthy, the fusion determinants of HBV GPs and, particularly, the fusion peptide that could induce merging of viral and endosomal membranes, have not yet been functionally identified in infection or cell-cell fusion assays.
First, to address which GP is responsible for HBV membrane fusion, we evaluated the role of either protein in cell-cell fusion assays (Figure 1D). Huh7 cells were transfected with plasmids encoding wild-type (wt) HBV GPs, that is, L, M, and S (pT7HB2.7 plasmid) vs only L, M, or S (using pCiL, pCiM, and pCiS plasmids, respectively) (Komla-Soukha and Sureau, 2006). To analyze the expression of either protein at cell surface, transfected cells were labeled with sulfo-NHS-SS-biotin, a chemical compound that is unable to penetrate biological membranes. After lysis and immunoprecipitation of biotinylated proteins, we found that the individually expressed L, M, or S proteins were detected at similar levels as compared to HBV GPs (L, M, and S) expressed simultaneously, as in cells transfected with the wt pT7HB2.7 plasmid (Figure 1E and F). Then, to determine the fusion activity of either protein, we performed cell-cell fusion assays as described above. We found that none of the L, M, or S proteins expressed alone were able to induce membrane fusion (Figure 1D). Furthermore, when we tested the pT7HB2.7Mless plasmid, which induces co-expression of S and L only (‘noM’ in Figure 1D–F), we detected a cell-cell fusion activity at the same level than for wt HBV GPs (Figure 1D). This indicated that M is not necessary for membrane fusion, in agreement with previous results (Ni et al., 2010; Sureau et al., 1994) showing that M is dispensable for infectivity of viral particles (Figure 1—figure supplement 3).
Altogether, these results suggested that the determinants of membrane fusion are harbored within L and S GPs.
Next, aiming to identify a fusion peptide in either protein, we used a computational approach to pinpoint regions of the HBV GPs that may potentially interact with membrane bilayers. Using Membrane Protein Explorer (MPEx), a tool based on the Wimley-White interfacial hydrophobicity scale (Snider et al., 2009), five regions of high interfacial hydrophobicity were identified (Figure 2—figure supplement 1A). Two out of the five hydrophobic regions did not correspond to HBV GP transmembrane regions (TM1, TM2, and TM3/TM4) and, therefore, were considered as candidate fusion peptides (Figure 2A and B). The first predicted segment comprised amino acids 48–66 that overlap with the preS1 domain. The second segment, which includes amino acids 127–145, is included in the preS2 region. Our prediction analyses indicated that the first segment (ΔG = -3.38) was more likely to be a fusion peptide than the second one (ΔG = -0.85) (Figure 2—figure supplement 1B). Considering the Wimley-White scale, a set of mutants was designed to alter the hydrophobicity of the two predicted segments (Figure 2B and Figure 2—figure supplement 1B). In the first segment, three mutants were studied by changing the aromatic residues to an alanine or glutamate: F52A, F56A, W66A, F52A/W66A (FW/AA), and F52E/W66E (FW/EE), or a glycine to an alanine (G53A). In the second segment, four mutants were considered: Y129A, F130A, S136E, and L144A; while the first two mutants targeted aromatic residues, S136 and L144 were also considered important because they are at the center of the predicted region and have a relatively high hydrophobicity.
-
Figure 2—source data 1
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig2-data1-v2.xlsx
-
Figure 2—source data 2
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig2-data2-v2.xlsx
-
Figure 2—source data 3
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig2-data3-v2.xlsx
To evaluate the role of these two sequences in HBV fusion, we introduced these single or double mutations in both regions and inserted them in the pT7HB2.7 HBV GP expression plasmid. Each mutant was compared to wt HBV GPs in both infection assays, using HDV particles (Sureau, 2010; Perez-Vargas et al., 2019), and cell-cell fusion assays, as described above. We found that HDV particles carrying these mutant GPs were produced by Huh7 cells at levels similar to those produced with wt GPs (Figure 2C and D), hence ruling out gross misfolding induced by the mutations that would otherwise prevent HBV GP incorporation into viral particles (Abou-Jaoudé and Sureau, 2007). Interestingly, no infectivity could be detected for most of the mutations introduced in the preS1 peptide (Figure 2E), whereas the HDV particles with mutations in the preS2 peptide showed levels of infectivity that were similar to those obtained with the wt GPs (Figure 2F). Correlating with the results of these infection assays, we found that the mutants in the preS1 peptide that prevented HDV infectivity also abrogated cell-cell fusion activity (Figure 2G) in a manner unrelated to the levels of GP cell-surface expression (Figure 2I and Figure 2—figure supplement 2). In contrast, mutations in the preS2 peptide displayed the same levels of cell-cell fusion activity as compared to wt (Figure 2H and J).
Altogether, these results indicated that the preS1 region harbors a potential fusion peptide.
Stabilizing cross-strand disulfide exchanges in HBV S protein prevents membrane fusion
Next, we sought to investigate the mechanisms that could induce fusion-activating conformational changes in the HBV GPs, leading to exposure of the fusion peptide. As neither the HBV receptor interaction nor the acidic pH could trigger membrane fusion (Figure 1), we thought that conformational rearrangement of HBV GPs might involve reshuffling of their disulfide bonds. Indeed, previous studies have shown that cysteine residues of the HBV S antigenic loop are essential for HDV infectivity and that viral entry is blocked by inhibitors of thiol/disulfide exchange reactions, such as Tris(2-carboxyethyl)phosphine hydrochloride (TCEP), dithiothreitol (DTT), 5,5-dithiobis(2-nitrobenzoic acid) (DTNB), or 4-acetamido-4'-maleimidyl-stilbene-2,2'-disulfonate (AMS) (Abou-Jaoudé and Sureau, 2007). Thus, to extend the notion that thiol/disulfide exchange reactions are implicated during membrane fusion and entry, we performed HBV infection and fusion assays in the presence of DTNB, an alkylating agent. First, using different DTNB concentrations that were added either at the onset of infection or at 16 hr post-infection, we confirmed that DTNB could block HDV infection in a dose-dependent manner, but only when it was added at the onset of infection (Abou-Jaoudé and Sureau, 2007; Figure 3—figure supplement 1). Second, using time-of-addition experiments, we found that DTNB could inhibit infection only if added within the first 2 hr after inoculation with HDV particles (Figure 3A). These results suggested that DTNB blocks a thiol/disulfide exchange reaction that could be necessary at an early step of infection, such as a trigger of the fusion mechanism, though not at a later stage of the entry process. Third, to evaluate the effect of DTNB on membrane fusion, we performed cell-cell fusion assays in presence of DTNB, which was added at the onset of cell co-cultures vs at 16 hr after seeding the cell co-cultures. We showed that DTNB added during the co-culture neither induced cytotoxicity (Figure 1—figure supplement 2) nor affected expression of HBV glycoproteins on the cell surface (Figure 3C and D). Yet, we found a dose-dependent reduction in the level of cell-cell fusion when DTNB was added immediately after cell-cell contact, whereas we detected a much lower effect in fusion activity when DTNB was added at 16 hr after cell contact (Figure 3B).
-
Figure 3—source data 1
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig3-data1-v2.xlsx
-
Figure 3—source data 2
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig3-data2-v2.xlsx
-
Figure 3—source data 3
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig3-data3-v2.xlsx
Altogether, these results suggested a role of the disulfide bond network during HBV membrane fusion steps, perhaps at the level of the fusion trigger.
To address this possibility and to identify potential mechanisms involved in fusion triggering, we focused on the ‘a’ determinant of protein S that exhibits eight conserved Cys, which, for some of them, are in strong proximity in the sequence (Figure 4A). To avoid trivial contact predictions between consecutive Cys, we defined four Cys-containing regions in a way that Cys pairs that are potentially in contact should have a sequence separation of at least four amino acids. The first Cys-containing region includes C270, the second, C284 and C287, the third, C300, C301, and C302, and the last one, C310 and C312 (Figure 4A). We applied secondary and tertiary structure prediction methods together with the contact prediction method RaptorX (Ma et al., 2015; Wang et al., 2017), based on coevolution signals, to predict disulfide connectivity in the ‘a’ determinant, which may identify which Cys forms disulfide bonds. Notably, RaptorX predicted structural contacts between either region (Figure 4—figure supplement 1) and we highlighted pairs of residues in contact in the four Cys regions, with the strongest signal detected between the third and fourth regions (Figure 4B). Next, applying the JPred secondary structure prediction method (Cole et al., 2008), we predicted two β-strands in the Cys-rich regions delimited by the S segments 298–303 and 310–313 (Figure 4A). Then, considering the secondary structure prediction and the contact prediction, we built a three-dimensional model for the region 294–317 (Figure 4C), which indicated that this sequence is compatible with a β-hairpin structural motif containing a cross-strand disulfide (CSD) bond between C301 and C310. Finally, through the analysis of its five χ dihedral angles (Figure 4—figure supplement 2), this disulfide bond was classified in a ‘-H Staple conformation’, which is a particular type of disulfide geometry associated with allosteric functions by triggering a conformational change upon switching between the reduced and oxidized states (Chiu and Hogg, 2019; Hogg, 2003).
We, therefore, hypothesized that the redox state of this disulfide may act as an allosteric switch that could contribute to control conformational rearrangements of the S protein. Thus, we used our structural model of the C301-C310 disulfide bond to design a mutant of S that could disrupt this hypothetical allosteric function, that is, the T303C/G308C double mutant that induces an additional C303-C308 disulfide bond (Figure 4C). Further molecular dynamic (MD) simulations (1000 frames per MD trajectory) carried out to differentiate between allosteric and structurally stabilizing disulfides, where the disulfides can be classified based on their angles (Figure 4—figure supplement 2), showed that the T303C/G308C mutant predominantly forms a structural disulfide bond.
Aiming to validate our prediction that an additional disulfide bond between the two β-strands could, by stabilizing the 298–313 β-hairpin motif, prevent membrane fusion from occurring, we produced HDV particles carrying the individual (T303C or G308C) and double (T303C/G308C) mutations in HBV GPs. By measuring HDV RNAs in cell supernatants, we found that all mutants could produce comparable levels of viral particles relative to wt virus (Figure 5A), suggesting absence of gross alterations of HBV GP conformation that would otherwise preclude virion assembly (Abou-Jaoudé and Sureau, 2007). Importantly, we found that while HDV particles generated with GPs harboring the individual mutations were as infectious as wt, those that were produced with the T303C/G308C double mutation (noted TG/CC in Figure 5) and the putative additional C303-C308 CSD bond were not infectious (Figure 5B). Moreover, we found that HDV particles harboring GPs with this T303C/G308C mutation had similar binding levels on Huh7 cells than those generated with wt GPs (Figure 6A), underscoring a post-binding defect. Then, to address this possibility, we performed cell-cell fusion assays with either HBV GP mutant, which was readily expressed at the cell surface (Figure 5C). We found that whereas the single mutations displayed similar fusion activity as compared to wt HBV GPs, the T303C/G308C double mutation completely prevented HBV GP-induced cell-cell fusion activity (Figure 5D).
-
Figure 5—source data 1
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig5-data1-v2.xlsx
-
Figure 5—source data 2
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig5-data2-v2.xlsx
-
Figure 5—source data 3
- https://cdn.elifesciences.org/articles/64507/elife-64507-fig5-data3-v2.xlsx
Altogether, these results suggested that the putative C303-C308 additional disulfide bond stabilizing the loop containing the C301-C310 CSD bond inhibited HBV entry and fusion, perhaps by preventing conformational rearrangements of HBV GPs that are required for promoting membrane fusion.
ERp57 is a protein disulfide isomerase that promotes HBV entry and infectivity in vivo
We reasoned that isomerization of the C301-C310 CSD (Figure 4) or of another CSD of the AGL determinant with allosteric functions could facilitate some conformational rearrangements required to promote membrane fusion. We, therefore, hypothesized that such an isomerization could be induced by a host factor from the PDI family, which are enzymes that can both reduce and oxidize disulfide bonds.
To address if PDIs are involved in HBV entry, we tested the effect of inhibitors (NTZ, EGCG, rutin, bacitracin, PX-12) that target different PDI species (PDIA1, ERp5, ERp57, TMX1) for their effect in cell entry of viral particles. First, through binding assays of viral particles to Huh7 or Huh7-NTCP cells performed in the presence of either inhibitor, we found that none of these inhibitors affected binding of HDV particles generated with either wt or T303C/G308C mutant GPs (Figure 6A). Second, using infection assays with Huh7-NTCP cells pre-incubated with either inhibitor, we found that HDV particles had strongly reduced infectivity in presence of nitazoxanide (NTZ) and (−)-epigallocatechin 3-gallate (EGCG) inhibitors that both target ERp57 (Figure 6B). Third, we confirmed these results using infection assays with authentic HBV particles (Figure 6C). Finally, to demonstrate that the inhibitors acted at the level of membrane fusion, we performed cell-cell fusion assays, as described above, whereby either inhibitor was added at the onset of co-cultures of HBV GP-expressing Huh7 donor and Huh7-NTCP-tat indicator cells and was kept throughout the assay period. Remarkably, we found a strong reduction in the levels of cell-cell fusion with the same drugs that inhibited HDV infection (Figure 6D). Hence, these results suggested a potential role of ERp57 in HBV membrane fusion.
Next, aiming to confirm and extend these findings, we selected a subset of the above PDIs, that is, ERp46, ERp57, and ERp72, which displayed low but significant expression at the surface of Huh7 cells (Figure 7A), in agreement with a previous report (Turano et al., 2002). We down-regulated either PDI in target cells via transduction of Huh7-NTCP cells with short hairpin RNA (shRNA)-expressing lentiviral vectors (Figure 7—figure supplement 1 and Figure 7—figure supplement 2). We found that down-regulation of ERp57, though not ERp46 or ERp72, strongly reduced the levels of HDV (Figure 7B) and HBV infection (Figure 7C) and of cell-cell fusion (Figure 7D). Finally, through confocal microscopy analysis of Huh7-NTCP cells (Figure 8A), we investigated the colocalization between ERp57 and Rab5 (early endosomes), Rab7 (late endosomes), Rab11 (recycling endosomes), or Lamp1 (lysosomes). The quantifications of these results showed that ERp57 could be detected in late endosomes but poorly in the other above-tested locations (Figure 8B), in line with the notion that fusion of HBV particles occurs in late endosomes (Macovei et al., 2013). Thus, ERp57 can be found at locations compatible for both cell-cell fusion and cell-free entry by internalization.
Altogether, these results indicated that ERp57 is likely a PDI that promotes HBV entry at a membrane fusion step.
Finally, we sought to demonstrate that ERp57 inhibition may prevent HBV propagation in vivo using NTZ, which has a short half-life of about 1.5 hr (Ruiz-Olmedo et al., 2017; Stockis et al., 1996). We generated a cohort of liver humanized mice (HuHep-mice) derived from the NFRG mouse model (Azuma et al., 2007; Figure 9A). We retained the animals that displayed >15 mg/ml of human serum albumin (HSA), which corresponds to 40–70% of human hepatocytes in the liver (Calattini et al., 2015). In agreement with previous reports (Perez-Vargas et al., 2019), these animals supported HBV infection (Group #1) for several weeks/months (Figure 9B; see Figure 9—figure supplement 1 for individual mice). The second group of HuHep mice was treated with NTZ 30 min prior to inoculation with HBV and, then, treated again with NTZ 1 hr later. We found that viremia in this group was delayed by about 4 weeks, as compared to Group #1 for which HBV could disseminate immediately after inoculation. This indicated that the blocking of ERp57 could prevent HBV infection in vivo.
Discussion
The entry process of enveloped viruses into cells is the series of all events that take place from the attachment of the virus to the host cell until the release of its genome into the cytoplasm. It is a finely regulated and complex process with several steps, in which many viral and cellular factors are involved. The first interaction often occurs with HSPGs. It may lack specificity but serves to give a virus an initial catch hold from which it can recruit specific receptors and entry co-factors that drive the reactions leading to entry. Fusion is the last step of enveloped virus entry and allows the release of the viral capsid into the cytoplasm following the merging of the viral membrane with the membrane of the infected cell. The interactions with the target cell that trigger conformational changes of the viral surface glycoproteins, ultimately leading to the insertion of their fusion peptide into the cell membrane, vary widely for enveloped viruses and can be divided into different scenarios. In the first one (e.g., human immunodeficiency virus), fusion is triggered directly by the interaction of the viral glycoprotein with its cellular receptor, through allosteric conformational rearrangements. In some cases, a sequential interaction with additional host factors is required to trigger the conformational changes required for fusion. In the second scenario (e.g., influenza virus), the interactions with the receptor at the cell surface leads to the endocytosis of viral particles, which is followed by GP protonation in the low-pH environment of the intracellular endosomal organelles that triggers the fusogenic conformational change. In the third scenario (e.g., Ebola virus), the initial interactions of the virion with the cell surface trigger its endocytosis followed by a second interaction with an internal receptor, often found in late endosomes, which is preceded by proteolytic cleavage of the fusion protein by an endosomal protease, leading to fusion activation (Harrison, 2015; White and Whittaker, 2016). Finally, for certain viruses (e.g., Sindbis virus), the fusion protein requires an activating redox reaction involving disulfide bonds of their glycoproteins to induce membrane fusion (Key et al., 2015; Rey and Lok, 2018).
Using a cell-cell fusion assay, we found that HBV fusion activity reached similar levels whether indicator cells expressed or not HSPG and/or NTCP but was not increased when the cell co-cultures were exposed at low pH, in contrast to bona fide pH-dependent GPs such as VSV-G or CCHFV Gn/Gc (Figure 1). That both HSPG and NTCP, which are respectively HBV virion membrane capture molecules (Leistner et al., 2008; Schulze et al., 2007) and specific entry factors (Ni et al., 2014; Yan et al., 2015), are not required for cell-cell fusion highlights that this fusion assay reveals late entry events, such as those occurring after virus interaction with either factor. Similarly, for other viruses such as influenza virus or hepatitis C virus (HCV), binding to their cell entry receptor is not a requirement in both cell-cell fusion (Lin and Cannon, 2002) and liposome fusion (Lavillette et al., 2006) assays triggered by a low-pH treatment. Thus, while it is clear that cell-cell fusion does not recapitulate per se all the events required to promote cell entry of viral particles since it bypasses the step of internalization that subsequently allows membrane fusion in endosomes that are required for entry of the above-mentioned viruses and of HBV (Macovei et al., 2013; Iwamoto et al., 2019), it is a suitable experimental tool for investigating some of the structural and functional determinants that promote envelope glycoprotein membrane-fusion activity (Earp et al., 2004). Accordingly, our results indicate that the trigger for the HBV membrane fusion mechanism not only is independent of an allosteric interaction of its GPs with the NTCP receptor but also is independent of GP protonation that is induced by the low-pH environment of endosomes. That low pH does not increase HBV cell-cell fusion agrees with previous results indicating that pharmacological agents that raise or neutralize the pH of the endocytic pathway had no effect on HBV infection (Macovei et al., 2010; Macovei et al., 2013; Rigg and Schaller, 1992).
Previous results from Abou-Jaoudé and Sureau, 2007 showed that cysteine residues of the HBV antigenic loop are essential for HDV infectivity and that viral entry is blocked by inhibitors of thiol/disulfide exchange reaction. Our results extend this notion as they indicate that such reactions seem to be necessary to mediate a critical early post-binding event but not at a later stage of the infection process since no effect in virus infectivity could be detected when DTNB was added at 4 hr post-infection (Figure 3). Since isomerization of disulfide bonds has been shown to be crucial for conformational rearrangements of GPs from other enveloped viruses leading to fusion (Fenouillet et al., 2007), we sought to investigate if and how such reactions could be implicated during the membrane fusion step of HBV entry. Here, using our cell-cell fusion assay, we found that DTNB blocked HBV GP-mediated membrane fusion (Figure 3B). Altogether, these results indicated a role of disulfide bond network of S GP during HBV membrane fusion.
Capitalizing on the above-mentioned experimental information that inhibitors of thiol/disulfide exchange reactions alter virus entry, we sought to examine how disulfide bonds of the HBV GPs, or rather, how a potential reshuffling of its disulfide bond profile, could be important for HBV entry. Indeed, cross-strand disulfides occurring in some viral surface GPs are believed to play a role in virus entry (Barbouche et al., 2003; Jain et al., 2007; Rosenthal et al., 1998; Wouters et al., 2004). Particularly, allosteric disulfide bonds can modulate the activity of the proteins in which they reside by mediating a structural change when they are reduced or oxidized (Hogg, 2003; Schmidt et al., 2006). Allosteric control of protein function is defined as a change in one site (the allosteric site) that influences another site by exploiting the protein’s flexibility; an allosteric disulfide bond represents the ‘allosteric site’ and the conformational change triggered by cleavage of such bonds alters protein function. For the HBV S protein, we used the contact prediction method RaptorX (Ma et al., 2015; Wang et al., 2017) to predict contacts between four Cys-rich regions of the AGL determinant (Figure 4), which highlighted that two of these regions may likely interact, that is, the Cys-rich regions III and IV (Figure 4—figure supplement 1). Using the secondary structure prediction method JPred (Cole et al., 2008), we proposed that these regions organize in two β-strands and we constructed a three-dimensional model of the region 294–317 of the HBV S GP, which indicated that this sequence is compatible with a β-hairpin structural motif containing a CSD bond between C301 and C310 (Figure 4). Interestingly, the analysis of the signs of the five χ dihedral angles defined by the Cys residues allowed to classify this particular disulfide bond in a -H Staple conformation, which is a particular type of disulfide geometry associated with allosteric functions that is known to trigger conformational changes upon switching between the reduced and oxidized states (Chiu and Hogg, 2019; Schmidt et al., 2006). Hence, we hypothesized that the redox state of the C301-C310 disulfide bond may act as an allosteric switch controlling conformational rearrangements of the HBV GP leading, ultimately, to exposure of the fusion peptide. Of note, the β-hairpin region with the predicted CSD lies at the surface of the S protein according to a three-dimensional in silico model (van Hemert et al., 2008), which may allow interactions with other HBsAg subunits and/or cellular factors. Yet, to test our structural and dynamic model involving a C301-C310 CSD bond in S GP (Figure 4), we reasoned that creating an additional, neighboring disulfide bond between positions 303 and 308 may stabilize the β-hairpin motif (Figure 4), which may prevent molecular rearrangements and, thus, allow membrane fusion to occur. The in silico analysis indicates that the T303C/G308C double mutant most probably generates two structural CSD according to our JPred-based (Cole et al., 2008) structural modeling, which affects the structural conformation of the C301-C310 CSD that is no longer classified as an allosteric bond. When we tested the T303C/G308C mutation in functional assays, we found that the mutant HBV GPs induced an almost complete loss of infection and fusion activity (Figure 5), hence suggesting that by stabilizing cross-strand disulfide exchange, the putative additional disulfide bond prevented conformational rearrangements of HBV GPs that are required for promoting membrane fusion. One possibility is that such stabilization could prevent an isomerization of the C301-C310 CSD bond that generates alternative disulfide bond(s), for example, between C284 and C310, which was proposed in a previous study (Mangold et al., 1995). Yet, the antigenic loss of S induced by these mutations did not allow us to design an assay that would detect the additional disulfide bond in the T303C-G308C mutant nor the block of conformational rearrangements that is suggested by its phenotype.
Assuming that the isomerization of the C301-C310 allosteric CSD or of other thiols/disulfides of the AGL determinant could facilitate the conformational rearrangements of HBV GPs required to promote membrane fusion, we hypothesized that such an isomerization could be induced by a host factor from the PDI family, which are enzymes that can both reduce and oxidize disulfide bonds. PDIs consist of a family of 21 structurally related proteins with a thioredoxin-like domain. Most of these isomerases have a CXXC motif that catalyzes formation, reduction, and rearrangement of the disulfide bonds in proteins (Abell and Brown, 1993). These isomerases are primarily involved in the folding of proteins in the endoplasmic reticulum (ER), catalyzing formation of their disulfide bonds, and most of these isomerases have ER retention signals. However, some isomerases from the PDI family have also been shown to be present at the cell surface, both in functional and in biochemical assays (Turano et al., 2002). Accordingly, cell-surface-localized PDIs are involved in processes such as cell adhesion, nitric oxide signaling, and cell entry of different viruses (Diwaker et al., 2013; Fenouillet et al., 2007). In support of the notion that PDIs are involved in HBV entry, we found that inhibitors that target different PDI members could block infection and cell-cell fusion though not the binding of viral particles to the cell surface. Of note, we found that bacitracin, which targets PDIA1, did not inhibit HBV entry and membrane fusion, in line with a previous study showing that it could not inhibit HDV entry (Abou-Jaoudé and Sureau, 2007). While the above ruled out PDIA1 as an entry co-factor of HBV, we found a strong reduction in the levels of HBV and HDV infection as well as of HBV-induced cell-cell fusion when we used the NTZ and EGCG inhibitors (Figure 6), which target ERp57 (Müller et al., 2008; Pacello et al., 2016). Consistently, we detected a low but significant expression of ERp57 as well as of ERp46 and ERp72 at the cell surface (Figure 7), in agreement with a previous study (Turano et al., 2002). Furthermore, we detected ERp57 in late endosomes (Figure 8), which is meaningful since previous reports have shown that HBV infection of HepaRG cells depends on Rab5 and Rab7 (Macovei et al., 2013), which are GTPases involved in the biogenesis of endosomes, and that the epidermal growth factor receptor is a host-entry cofactor that interacts with NTCP and mediates HBV internalization (Iwamoto et al., 2019). Using a gene silencing approach, we confirmed that down-regulation of ERp57 but not of these alternative PDIs could decrease the levels of HDV and HBV infection as well as of cell-cell fusion (Figure 7). Importantly, we showed that a short-time treatment of liver humanized mice with NTZ could delay HBV infection by approximately 2–4 weeks (Figure 9). Since NTZ has a short half-life of about 1.5 hr in vivo (Ruiz-Olmedo et al., 2017; Stockis et al., 1996) and since NTZ was administrated at very short times before and after HBV inoculation, we calculated that less than 10% of the drug was still present in those mice at 7 hr post-infection, which likely precludes an effect on HBV post-entry steps (Korba et al., 2008). Altogether, these results support the role of ERp57 at early steps of HBV infection and validate this PDI as a therapeutic target. Note that our results did not discard the possibility that some other PDIs could also play a role during HBV entry into cells.
The fusion-mediating GPs of enveloped viruses contain a sequence, termed fusion peptide, that interacts with and destabilizes the cellular target membrane. Such an event is finely controlled so as to occur at the appropriate time and location and to prevent fortuitous inactivation of GP fusion activity and virus infectivity. Hence, a conformational change in these GPs is a requirement to induce the accessibility and function of the fusion peptide segment. Candidate fusion peptides are generally identified as hydrophobic sequences, of approximately 16–26 residues in length, that are conserved within a virus family and that may adopt α-helical conformation with strongly hydrophobic faces. They can be internal or located at the amino-terminus of fusion GP subunits (Apellániz et al., 2014; Epand, 2003; Martin et al., 1999). There are a number of criteria that characterize fusion peptide segments and, while none of these criteria taken individually are absolute to define a fusion peptide segment, they are sufficiently restrictive to predict if a given region of a protein presents features of a fusion peptide segment (Delos and White, 2000; Delos et al., 2000), which needs to be further functionally tested.
Previously, a conserved peptide comprising 23 amino acids at the N-terminal end of the HBV S protein and overlapping its TM1 sequence was shown to interact with model membranes, causing liposome destabilization in a pH-dependent manner (Berting et al., 2000; Rodríguez-Crespo et al., 1994; Rodríguez-Crespo et al., 1995; Rodríguez-Crespo et al., 1999). However, it was also demonstrated that hydrophobic residues in TM1 were critical for S protein expression as well as for infectivity (Siegler and Bruss, 2013). An essential role of TM1 in fusion mechanism, albeit in a pH-independent manner, could be shown for the duck hepatitis B virus (DHBV) (Chojnacki et al., 2005; Grgacic and Schaller, 2000), although there is also evidence for the involvement of the preS domain of DHBV at an early step of infection, likely during the fusion process (Delgado et al., 2012).
Here, through a computational hydropathy analysis of the HBV GPs, we identified two potential short sequences within the preS1 and preS2 regions that may potentially interact with membrane bilayers. To validate these predictions, we characterized in both infection and cell-cell fusion assays HBV GP mutants in key positions in either sequence. We found that while none of the mutations in the preS2 segment altered infection or membrane fusion activities, mutations in the preS1 sequence induced an almost complete loss of infectivity and cell-cell fusion (Figure 2). Note that these mutants had similar if not identical levels of cell-surface-expressed L, M, and S proteins and/or capacity to induce the formation of HDV particles. These results suggested that the preS1 region harbors a fusion peptide in addition to the NTCP-binding determinant.
Overall, our study characterizes some crucial determinants of HBV entry and membrane fusion. The mechanism by which fusion proteins are activated and undergo conformational rearrangements or fusion intermediates is a particularly complex process involving several regions of viral surface GPs. Our results suggest that for HBV, this mechanism could be triggered by ERp57, allowing a thiol/disulfide exchange reaction to occur and regulate isomerization of critical CSD(s), which ultimately results in the exposition of the fusion peptide that seems to be located within the preS1 region.
Materials and methods
Plasmids
Plasmid pSVLD3 harboring a trimer of the HDV gt1 genome (accession number M21012.1), pCiL encoding the L protein, pCiS encoding the S protein (Komla-Soukha and Sureau, 2006), and pT7HB2.7 encoding the three HBV envelope proteins were a gift from Sureau, 2010; Sureau et al., 1994. To induce the expression of L and S only, the pT7HB2.7 plasmid was modified at the M start codon and Kozak consensus sequence in order to silence the expression of M protein, resulting in pT7HB2.7Mless construct. The pCiM plasmid encoding the M protein was constructed by deleting the preS1 region from pCiL until the N-terminal methionine of preS2. All mutations in pT7HB2.7 plasmid were introduced by point-directed mutagenesis. The phCMV-VSV-G encoding the G protein from VSV and pCAGGS-GP/wt-M encoding the Gn and Gc glycoproteins from CCHFV were described previously (Freitas et al., 2020). The plasmid encoding the luciferase reporter under control of an HIV-1 long terminal repeat internal promoter (pLTR-luc) was a gift from Olivier Schwartz and was used as described before (Lavillette et al., 2007). Mission shRNA plasmids (Sigma), shRNA sequences, and oligonucleotides used for introducing mutations in HBV GPs are described in Supplementary file 1.
Cells
Huh7 human hepatocarcinoma cells and Huh7-NTCP cells, which were generated by transduction of Huh7 cells with a retroviral vector transducing the NTCP plasmid (pLX304NTCP, DNASU) and selected for blasticidin resistance, were grown in William’s E medium (WME) (Gibco, France) supplemented with non-essential amino acids, 2 mM L-glutamine, 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) buffer, 100 U/ml of penicillin, 100 µg/ml of streptomycin, and 10% fetal bovine serum. 293T human kidney cells (ATCC CRL-1573), CHO-K1 (CHO) Chinese hamster ovary cells (ATCC CCL-61), and CHO-pgsB-618 cells (ATCC CRL-2241), which do not produce glycosaminoglycans, were grown in Dulbecco’s modified minimal essential medium (DMEM, Gibco) supplemented with 100 U/ml of penicillin, 100 µg/ml of streptomycin, and 10% fetal calf serum. Huh7-tat and Huh7-NTCP-tat indicator cells expressing HIV Tat were generated by transduction of Huh7 and Huh7-NTCP cells, respectively, with the LXSN-tat retroviral vector and selected for G418 resistance. HepG2.2.15 human hepatoma cells were used to produce HBV virus and were maintained in WME medium complemented with 10% fetal bovine serum. Authentication of purchased cell lines was performed by ATCC. Authentication of Huh7 cells was based on expression of human transferrin and serum albumin. Authentication of HepG2.2.15 cells was based on titration of released infectious HBV particles. All cell lines were certified mycoplasma-free, as per our monthly contamination testing.
PDI inhibitors
Request a detailed protocolDTNB, NTZ, EGCG, rutin, bacitracin, and PX-12 were purchased from Sigma-Aldrich and dissolved in dimethyl sulfoxide (DMSO), ethanol, or water according to the manufacturer’s instructions.
Antibodies
For western blot analysis, HBs antigen and calnexin were detected with goat anti-HBV polyclonal antibody (Murex, DiaSorin) coupled to horseradish peroxidase (HRP) and rabbit calnexin polyclonal antibody (Enzo), respectively. The rabbit anti-ERp46 (Abcam), mouse anti-ERp57 (Abcam), and mouse anti-ERp72 (Santa Cruz Biotechnology) antibodies were used for detecting PDI proteins by flow cytometry and western blot. NTCP was detected with polyclonal NTCP/SLC10A1 antibody (Bioss Antibodies) coupled to Phycoerythrin (PE) for flow cytometry. Mouse anti-ERp57 (Abcam), Rabbit anti-Rab5, anti-Rab7, anti-Rab11, and anti-Lamp1 (Cell Signaling Technology), and Donkey anti-Rabbit-Alexa-Fluor-488 and Donkey anti-Mouse-Alexa-Fluor-568 (Thermo Fisher) antibodies were used for immunofluorescence (IF) studies.
shRNA-expressing stable cell lines
Request a detailed protocol293 T cells were seeded 24 hr prior to transfection with VSV-G plasmid, psPAX2 packaging plasmid, and pLKO.1 expression vector carrying shRNA against ERp46, ERp57, or ERp72 using calcium phosphate precipitation. Medium was replaced 16 hr post-transfection. Vector supernatants were harvested 24 hr later, filtered through a 0.45 μm filter. Stable knockdown of ERp72, ERp57, or ERp46 in Huh7-NTCP, Huh7-tat, and Huh7-NTCP-tat cells was performed by selection with puromycin after lentiviral transduction. The knockdown was validated by flow cytometry and western blot using antibodies against ERp46, ERp57, or ERp72.
Cell-cell fusion assays
Request a detailed protocolHuh7 ‘donor’ cells (2.5x105 cells/well seeded in six-well tissue culture dishes 24 hr prior to transfection) were co-transfected using FuGENE six transfection reagent (Promega) with 3 μg of pT7HB2.7- wt or mutated glycoproteins and 50 ng of pLTR-luc reporter plasmid. For a positive control, cells were co-transfected with 3 μg of either pCAGGS-GP/wt-M, expressing CCHFV GPs, or 1 μg of phCMV-VSV-G and with 50 ng of the pLTR-luc plasmid. For negative controls, cells were co-transfected with 2 μg of an empty phCMV plasmid and 50 ng of the pLTR-luc plasmid. 12 hr later, transfected cells were detached with Versene (0.53 mM ethylenediaminetetraacetic acid (EDTA); Gibco), counted, and reseeded at the same concentration (105 cells/well) in 12-well plates. Huh7-tat or Huh7-NTCP-tat indicator cells, detached with EDTA and washed, were then added to the transfected cells (3x105 cells per well). After 24 hr of cocultivation, the cells were washed with phosphate-buffered saline (PBS), incubated for 3 min in fusion buffer (130 mM NaCl, 15 mM sodium citrate, 10 mM MES [2-(N-morpholino)ethanesulfonic acid], 5 mM HEPES) at pH 4, pH 5, or pH 7, and then washed three times with normal medium. The luciferase activity was measured 24 hr later using a luciferase assay kit according to the manufacturer’s instructions (Promega).
HDV particle production and infection
Request a detailed protocolHuh7 cells were seeded in 10 cm plates at a density of 106 cells per plate and were transfected with a mixture of 2.5 µg of pSVLD3 plasmid and 10 µg of plasmid allowing the expression of surface envelope glycoproteins of VSV or HBV using FuGENE six transfection reagent (Promega), as described previously (Perez-Vargas et al., 2019). Transfected cells were grown for up to 9 days in primary hepatocyte maintenance medium containing 2% DMSO to slow cell growth.
The supernatants of virus producer cells were filtrated through 0.45-nm-pore filters and were analyzed by quantitative reverse transcription PCR (RTqPCR) for detection of HDV RNA, using the primers described below. These supernatants were also used for infection experiments in Huh7-NTCP cells or PDI-down-regulated Huh7-NTCP cells, which were seeded in 48-well plates at a density of 1.5x104 cells per well. Infected cells were cultured in primary hepatocyte maintenance medium containing 2% DMSO following infection. RTqPCR assays were used to assess infectivity of viral particles at 7 days post-infection.
For inhibition assays, drugs were incubated with cells for 2 hr at 37°C before virus addition or at different times post-infection and the infectivity was assessed 7 days post-infection by RT-qPCR.
Binding assays
Request a detailed protocolHDV wt particles (107genome equivalent [GE]) were added to Huh7-NTCP cells and incubated for 1 hr at 4°C. Unbound virions were removed by three washes with cold PBS, and RTqPCR was used to assess the amount of bound viral particles.
RTqPCR detection of HDV RNAs in virus producer and infected cells
Request a detailed protocolCells were washed with PBS and total RNA was extracted with TRI Reagent according to the manufacturer’s instructions (Molecular Research Center). RNAs were reverse transcribed using random oligonucleotide primers with iScript (Bio-Rad). The following specific primers were used: for HDV RNA quantification, forward primer 5’-GGACCCCTTCAGCGAACA and reverse primer 5’-CCTAGCATCTCCTCCTATCGCTAT. Quantitative PCR (qPCR) was performed using FastStart Universal SYBR Green Master (Roche) on a StepOne Real-Time PCR System (Applied Biosystems). As an internal control of extraction, in vitro-transcribed exogenous RNAs from the linearized Triplescript plasmid pTRI-Xef (Invitrogen) were added to the samples prior to RNA extraction and quantified with specific primers (5′-CGACGTTGTCACCGGGCACG and 5′-ACCAGGCATGGTGGTTACCTTTGC). All values of intracellular HDV RNAs were normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene transcription. For GAPDH mRNA quantification, we used the forward 5’-AGGTGAAGGTCGGAGTCAACG and reverse 5’-TGGAAGATGGTGATGGGATTTC primers.
Western blot analyses
Request a detailed protocolThe proteins from pelleted cell supernatants or extracted from total cell lysates were denatured in Lammeli buffer at 95°C for 5 min and were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS–PAGE) and then transferred to nitrocellulose membranes (GE Healthcare). Membranes were blocked with 5% nonfat dried milk in PBS and incubated at 4°C with a rabbit or mouse antibody diluted in PBS-0.01% milk, followed by incubation with an IRdye secondary antibody (Li-Cor Biosciences). Membrane visualization was performed using an Odyssey infrared imaging system CLx (LI-COR Biosciences).
For cell-surface biotinylation, Huh7 cells were transfected into 10 cm plates with plasmids encoding wt or mutant HBV GPs. After 48 hr, the cell monolayers were rinsed three times with ice-cold PBS and overlaid with 0.5 ml biotin solution (0.5 mg sulpho-N-hydroxysuccinimide–biotin (Pierce) per ml of PBS, pH 7.2). The cells were then labeled for 30 min at 4°C. The biotin solution was removed and the cells were rinsed once with ice-cold 100 mM glycine solution and then incubated for 15 min with 100 mM glycine at 4°C to stop the reaction. The last washing step was performed with ice-cold PBS. Proteins were solubilized by the addition of 1 ml radioimmunoprecipitation assay buffer and equivalent quantities of protein lysates from each sample (Nanodrop quantification; Thermofisher) were immunoprecipitated with biotin-agarose beads. Proteins were electrophoresed under reducing conditions in SDS–PAGE followed by electrophoretic transfer to nitrocellulose. Surface‐biotinylated proteins were detected with anti-HBV antibody (Murex) coupled to HRP and enhanced chemiluminescence (ECL; Roche). The membranes of biotinylated samples were routinely re-probed with anti-calnexin antibody to confirm the absence of the intracellular protein calnexin. In addition, 10% of each lysate was denatured and loaded onto separate gels. Immunoblotting for calnexin on the membranes of lysate was done to confirm uniform protein loading.
Densitometry analysis (Image Lab BioRad software) was used to estimate the relative total amount of L, M, and S mutant proteins, which were expressed relative of the wild-type L, M, and S total proteins.
Flow cytometry
Request a detailed protocolThe surface expression of NTCP, ERp46, ERp57, and ERp72 was quantified by fluorescence activated cell sorting (FACS) analysis from 106 live cells using antibodies added to cells for 1 hr at 4°C. After washing, the binding of antibody to the cell surface was detected using PE (phycoerythrin)-conjugated anti-mouse antibodies.
Immunofluorescence, confocal microscopy imaging, and image analysis
Request a detailed protocolHuh7-NTCP cells were grown on uncoated 14-mm-diameter glass coverslips. 48 hr after seeding, cells were washed with PBS, fixed with 3% paraformaldehyde in PBS for 15 min, quenched with 50 mM NH4Cl, and permeabilized with 0.1% Triton X-100 for 7 min. Fixed cells were then incubated for 1 hr with primary antibodies in 1% bovine serum albumin (BSA)/PBS, and washed and stained for 1 hr with the corresponding fluorescent Alexa Fluor conjugated secondary antibody (Alexa-Fluor-488 and Alexa-Fluor-568; Thermo Fisher) in 1% BSA/PBS. Cells were washed three times with PBS, stained for nuclei with Hoechst 33342 (Thermo Fisher) for 5 min, washed two times with PBS, and mounted on Mowiol 40–88 (Sigma-Aldrich) prior to image acquisition with LSM-710 confocal microscope (Zeiss). Single-section confocal images of 0.6 µm thickness were analyzed with the ImageJ software. The Pearson’s correlation coefficients were calculated by using the JACoP plugin for ImageJ.
Cytotoxicity assay
Request a detailed protocolThe release of lactate dehydrogenase (LDH) from damaged cells was measured with CytoTox-ONE (Promega, MA, USA) homogeneous membrane integrity assay. Cells were grown in a 96-well flat-bottom culture plate at a density of 3x103 cells per well and treated with the different drugs for 2 hr or 24 hr. Maximum LDH release was determined by adding 2 µl of CytoTox-ONE lysis solution to control wells for 10 min. The assay was performed in 96-well plates by adding 100 µl of the sample supernatant and 100 µl of CytoTox-ONE reagent, after which the plate was shaken for 10 s. After 10 min of incubation, 50 µl CytoTox-One stop solution was added and the plate was shaken again for 10 s. The fluorescence signal was measured at λEX = 560 nm, λEM = 590. LDH release was calculated as the percentage of LDH released in the culture media of total LDH (media and lysates).
Fusion peptide prediction
Request a detailed protocolThe HBV surface sequence used was taken from the UniProt database, with accession number P03138. Hydropathy plots were obtained with Membrane Protein eXplorer software (Snider et al., 2009) using as input the reference sequence. Hydropathy plots were also used to evaluate the effect of residue mutations. Sequences with a propensity to partition into the lipid bilayer were identified using interfacial settings and pH = 5.0.
Contact prediction on the Cys-rich region
Request a detailed protocolContact prediction was performed using RaptorX (Wang et al., 2017; Teppa et al., 2020). RaptorX integrates evolutionary coupling and sequence conservation information through an ultra-deep neural network formed by two deep residual neural networks. RaptorX predicts pairs of residues, whose mutations have arisen simultaneously during evolution.
Structural models and molecular dynamic simulation studies
Request a detailed protocolThe HBV surface protein sequence was taken from the UniProt database, with accession number P03138. Secondary structure prediction was performed with Jpred (Cole et al., 2008). The S protein region 294–317 was modeled using MODELLER (Sali and Blundell, 1993). The template crystal structure of the Newcastle disease virus fusion protein (PDB code: 1G5G) was retrieved from the PDB database (Berman et al., 2000). Sequence alignment was generated with Clustal X (Larkin et al., 2007). The model evaluation was conducted using the Ramachandran plot (Ramachandran et al., 1963). The model of the wild-type sequence was further used to create two structural models with mutations using UCSF Chimera package (Pettersen et al., 2004). One model contains the double mutations T303C/G308C, which may create an extra disulfide bond. The overall effect of those mutations would be to ‘shift’ the disulfide bridge of two amino acids toward the turn of the β-hairpin motif. After mutations, the models were energy minimized by applying Molecular Modelling Toolkit (MMTK) with Amber parameters for standard residues and 100 steepest descent minimization steps with a step size of 0.02 Å. To investigate the stability of the disulfide bonds, MD simulations of the three models were carried out by GROMACS version 2020 (Abraham et al., 2015) in conjunction with OPLS-AA/L all-atom force field. The models were immersed in the cubic boxes filled with water molecules with a minimal distance of 1.0 nm between the peptide surface and the box. Each system was equilibrated to the desired temperature through a stepwise heating protocol in NVT ensemble followed by 100.0 ps equilibration in NPT ensemble with position restraints on the protein molecule. The final productive MD was performed for each system for 10 ns under periodic boundary conditions without any restraints on the protein with a time step of 2 fs at constant pressure (1 atm) and temperature (300 K). Coordinates were saved every 10 ps, yielding 1000 frames per MD trajectory. All the frames were further investigated to differentiate between allosteric and structurally stabilizing disulfides. Disulfide bonds were classified based on the five relevant torsion angles (χ1, χ2, χ3, χ2, and χ1') (see Figure 4—figure supplement 2), disulfides being treated as symmetrical. In this system, 20 conformational categories are possible (Marques et al., 2010; Schmidt and Hogg, 2007; Schmidt et al., 2006). The three central angles (χ2, χ3, and χ2) define the basic shapes: Spiral, Hook, and Staple (Eklund et al., 1984). χ3 defines the orientational motif: left-handed (LH) or right-handed (RH) if the sign is negative or positive, respectively (Eklund et al., 1984). χ1 and χ1' determine the signs of the nomenclature (Qi and Grishin, 2005).
In vivo experiments
Request a detailed protocolAll experiments were performed in accordance with the European Union guidelines for approval of the protocols by the local ethics committee (Authorization Agreement C2EA-15, ‘Comité Rhône-Alpes d’Ethique pour l’Expérimentation Animale’, Lyon, France - APAFIS#27316–2020060810332115 v4). Primary human hepatocytes (PHH; Corning, BD Gentest) were intrasplenically injected into NFRG mice (Azuma et al., 2007), a triple mutant mouse knocked out for fumarylacetoacetate hydrolase (fah-/-), recombinase activating gene 2 (rag2-/-), and interleukin two receptor gamma chain (IL2rg-/-), 48 hr after adeno-uPA conditioning (Bissig et al., 2010; Calattini et al., 2015). Mice were subjected to 2-(2-nitro-4-trifluoro-methylbenzoyl)1,3-cyclohexedione (NTBC) cycling during the liver repopulation process, as described previously (Calattini et al., 2015). Mice with HSA levels >15 mg/ml, as determined using a Cobas C501 analyzer (Roche Applied Science), were inoculated with virus preparations by intra-peritoneal injections. Sera were collected at different time points before and after infection. Mice were sacrificed 6 weeks post-infection.
Statistical analysis
Request a detailed protocolStatistical analyses were performed using GraphPad Prism version 5.02 for Windows, GraphPad Software (San Diego, California, USA). The Mann-Whitney or Wilcoxon tests were used for statistical comparisons. A p-value of 0.05 or less was considered as significant. When applicable, data are presented as mean ± standard deviation and results of the statistical analysis are shown as follows: ns, not significant (p>0.05); *p<0.05; **p<0.01; and ***p<0.001.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for figures 1-3, 5-7 and 9.
References
-
The three lives of viral fusion peptidesChemistry and Physics of Lipids 181:40–55.https://doi.org/10.1016/j.chemphyslip.2014.03.003
-
Robust expansion of human hepatocytes in fah−/−/Rag2−/−/Il2rg−/− miceNature Biotechnology 25:903–910.https://doi.org/10.1038/nbt1326
-
Entry of hepatitis B and C viruses - recent progress and future impactCurrent Opinion in Virology 4:58–65.https://doi.org/10.1016/j.coviro.2013.12.002
-
The protein data bank and the challenge of structural genomicsNature Structural Biology 7:957–959.https://doi.org/10.1038/80734
-
Human liver chimeric mice provide a model for hepatitis B and C virus infection and treatmentJournal of Clinical Investigation 120:924–930.https://doi.org/10.1172/JCI40094
-
Hepatitis B virus infection is dependent on cholesterol in the viral envelopeCellular Microbiology 11:249–260.https://doi.org/10.1111/j.1462-5822.2008.01250.x
-
Functional and biochemical characterization of hepatitis C virus (HCV) Particles produced in a humanized liver mouse modelJournal of Biological Chemistry 290:23173–23187.https://doi.org/10.1074/jbc.M115.662999
-
Allosteric disulfides: sophisticated molecular structures enabling flexible protein regulationJournal of Biological Chemistry 294:2949–5908.https://doi.org/10.1074/jbc.REV118.005604
-
The jpred 3 secondary structure prediction serverNucleic Acids Research 36:W197–W201.https://doi.org/10.1093/nar/gkn238
-
Potential roles of protein disulphide isomerase in viral infectionsActa Virologica 57:293–304.
-
BookThe Many Mechanisms of Viral Membrane Fusion ProteinsIn: Marsh M, editors. Membrane Trafficking in Viral Replication. Berlin Heidelberg: Springer. pp. 25–66.
-
Fusion peptides and the mechanism of viral fusionBiochimica Et Biophysica Acta (BBA) - Biomembranes 1614:116–121.https://doi.org/10.1016/S0005-2736(03)00169-X
-
Cell entry by enveloped viruses: redox considerations for HIV and SARS-coronavirusAntioxidants & Redox Signaling 9:1009–1034.https://doi.org/10.1089/ars.2007.1639
-
Hepatitis B virus entry into HepG2-NTCP cells requires clathrin-mediated endocytosisCellular Microbiology 22:13205.https://doi.org/10.1111/cmi.13205
-
Disulfide bonds as switches for protein functionTrends in Biochemical Sciences 28:210–214.https://doi.org/10.1016/S0968-0004(03)00057-4
-
Targeting allosteric disulphide bonds in CancerNature Reviews Cancer 13:425–431.https://doi.org/10.1038/nrc3519
-
Hepatitis B surface antigen assembles in a post-ER, pre-Golgi compartmentJournal of Cell Biology 118:1305–1320.https://doi.org/10.1083/jcb.118.6.1305
-
The p10 FAST protein fusion peptide functions as a cystine noose to induce cholesterol-dependent liposome fusion without liposome tubulationBiochimica Et Biophysica Acta (BBA) - Biomembranes 1848:408–416.https://doi.org/10.1016/j.bbamem.2014.10.020
-
Clustal W and clustal X version 2.0Bioinformatics 23:2947–2948.https://doi.org/10.1093/bioinformatics/btm404
-
Hepatitis C virus glycoproteins mediate low pH-dependent membrane fusion with liposomesJournal of Biological Chemistry 281:3909–3917.https://doi.org/10.1074/jbc.M509747200
-
Role of glycosaminoglycans for binding and infection of hepatitis B virusCellular Microbiology 10:122–133.https://doi.org/10.1111/j.1462-5822.2007.01023.x
-
Regulation of hepatitis B virus infection by Rab5, Rab7, and the endolysosomal compartmentJournal of Virology 87:6415–6427.https://doi.org/10.1128/JVI.00393-13
-
Conformational characterization of disulfide bonds: a tool for protein classificationJournal of Theoretical Biology 267:388–395.https://doi.org/10.1016/j.jtbi.2010.09.012
-
Role of the N-terminal peptides of viral envelope proteins in membrane fusionAdvanced Drug Delivery Reviews 38:233–255.https://doi.org/10.1016/S0169-409X(99)00031-9
-
Enveloped viruses distinct from HBV induce dissemination of hepatitis D virus in vivoNature Communications 10:10117-z.https://doi.org/10.1038/s41467-019-10117-z
-
UCSF chimera--a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Structural classification of thioredoxin-like fold proteinsProteins: Structure, Function, and Bioinformatics 58:376–388.https://doi.org/10.1002/prot.20329
-
Stereochemistry of polypeptide chain configurationsJournal of Molecular Biology 7:95–99.https://doi.org/10.1016/S0022-2836(63)80023-6
-
Fibroblast growth Factor-2 can mediate cell attachment by linking receptors and heparan sulfate proteoglycans on neighboring cellsJournal of Biological Chemistry 270:24188–24196.https://doi.org/10.1074/jbc.270.41.24188
-
Duck hepatitis B virus infection of hepatocytes is not dependent on low pHJournal of Virology 66:2829–2836.https://doi.org/10.1128/jvi.66.5.2829-2836.1992
-
Prediction of a putative fusion peptide in the S protein of hepatitis B virusJournal of General Virology 75:637–639.https://doi.org/10.1099/0022-1317-75-3-637
-
Phospholipid interactions of the putative fusion peptide of hepatitis B virus surface antigen S proteinJournal of General Virology 76:301–308.https://doi.org/10.1099/0022-1317-76-2-301
-
Effect of nitazoxanide on albendazole pharmacokinetics in cerebrospinal fluid and plasma in ratsSaudi Pharmaceutical Journal 25:413–418.https://doi.org/10.1016/j.jsps.2016.09.005
-
Comparative protein modelling by satisfaction of spatial restraintsJournal of Molecular Biology 234:779–815.https://doi.org/10.1006/jmbi.1993.1626
-
A function essential to viral entry underlies the hepatitis B virus "a" determinantJournal of Virology 83:9321–9328.https://doi.org/10.1128/JVI.00678-09
-
Search for allosteric disulfide bonds in NMR structuresBMC Structural Biology 7:49.https://doi.org/10.1186/1472-6807-7-49
-
MPEx: a tool for exploring membrane proteinsProtein Science 18:2624–2628.https://doi.org/10.1002/pro.256
-
Pharmacokinetics of nitazoxanide after single oral dose administration in 6 healthy volunteersInternational Journal of Clinical Pharmacology and Therapeutics 34:349–351.
-
The use of hepatocytes to investigate HDV infection: the HDV/HepaRG modelMethods in Molecular Biology 640:463–473.https://doi.org/10.1007/978-1-60761-688-7_25
-
COVTree: coevolution in OVerlapped sequences by tree analysis serverNucleic Acids Research 48:W558–W565.https://doi.org/10.1093/nar/gkaa330
-
Proteins of the PDI family: unpredicted non-ER locations and functionsJournal of Cellular Physiology 193:154–163.https://doi.org/10.1002/jcp.10172
-
Inhibitory effect of silibinin on hepatitis B virus entryBiochemistry and Biophysics Reports 14:20–25.https://doi.org/10.1016/j.bbrep.2018.03.003
-
Accurate de novo prediction of protein contact map by Ultra-Deep learning modelPLOS Computational Biology 13:e1005324.https://doi.org/10.1371/journal.pcbi.1005324
-
NTCP opens the door for hepatitis B virus infectionAntiviral Research 121:24–30.https://doi.org/10.1016/j.antiviral.2015.06.002
-
Identification of allosteric disulfides from prestress analysisBiophysical Journal 107:672–681.https://doi.org/10.1016/j.bpj.2014.06.025
Article and author information
Author details
Funding
LabEx Ecofect of the Université de Lyon (ANR-11-LABX-0048)
- François-Loïc Cosset
ANR (ANR-11-IDEX-0007)
- François-Loïc Cosset
LabEX CALSIMLAB of Sorbonne Université (ANR-11-LABX-0037–01)
- Alessandra Carbone
ANR (ANR-11-IDEX-0004–02)
- Alessandra Carbone
ANRS|MIE (ECTZ71388)
- François-Loïc Cosset
- Anja Böckmann
ANRS|MIE (ECTZ160643)
- Alessandra Carbone
- François-Loïc Cosset
ANRS|MIE (ECTZ160726)
- Alessandra Carbone
- François-Loïc Cosset
- Natalia Freitas
ANRS|MIE (ECTZ38814)
- Alessandra Carbone
- François-Loïc Cosset
- Elin Teppa
ANRS|MIE (ECTZ41733)
- Alessandra Carbone
- François-Loïc Cosset
- Christophe Combet
ANRS|MIE (ECTZ123278)
- François-Loïc Cosset
ANRS|MIE (ECTZ119828)
- François-Loïc Cosset
Foundation FINOVI
- François-Loïc Cosset
2017 Call for Joriss Projects
- François-Loïc Cosset
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We are grateful to Camille Sureau for the HBV GP expression constructs and for sharing the HDV infection assay. We thank Solène Denolly for helpful discussions.
We thank the ‘Plateforme de Thérapie Génique’ in Nantes (France) for the production of the in vivo-certified lots of adeno-uPA vectors. We thank Jean-François Henry, Nadine Aguilera and Tiphaine Dorel from the animal facility (PBES, Plateau de Biologie Experimental de la Souris, UMS3444/CNRS, US8/Inserm, ENS de Lyon, UCBL), and Véronique Pierre for their technical help in the handling of mice. We acknowledge the contribution of ANIRA-Genetic Analysis and the PLATIM-Microscopy facilities of SFR Biosciences (UMS3444/CNRS, US8/Inserm, ENS de Lyon, UCBL) for image quantifications, technical assistance, and support. We thank Didier Décimo for support with the BSL3 facility. We thank Omran Allatif for guidance with the statistical analysis. This work was supported by the French ‘Agence Nationale de la Recherche sur le SIDA et les Hépatites Virales | Maladie Infectieuses Emergentes’ (ANRS|MIE, grants ECTZ71388, ECTZ160643, ECTZ160726, ECTZ38814, and ECTZ41733), the Foundation FINOVI, the 2017 Call for Joriss Projects, the LabEx Ecofect (ANR-11-LABX-0048) of the ‘Université de Lyon’, within the program ‘Investissements d’Avenir’ (ANR-11-IDEX-0007) operated by the French National Research Agency (ANR), and the LabEX CALSIMLAB (ANR-11-LABX-0037–01 and ANR-11-IDEX-0004–02) of Sorbonne Université.
Ethics
Animal experimentation: All experiments were performed in accordance with the European Union guidelines for approval of the protocols by the local ethics committee (Authorization Agreement C2EA-15, "Comité Rhône-Alpes d'Ethique pour l'Expérimentation Animale", Lyon, France - APAFIS#27316-2020060810332115 v4).
Copyright
© 2021, Pérez-Vargas et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,926
- views
-
- 326
- downloads
-
- 18
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Microbiology and Infectious Disease
Timely and effective use of antimicrobial drugs can improve patient outcomes, as well as help safeguard against resistance development. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) is currently routinely used in clinical diagnostics for rapid species identification. Mining additional data from said spectra in the form of antimicrobial resistance (AMR) profiles is, therefore, highly promising. Such AMR profiles could serve as a drop-in solution for drastically improving treatment efficiency, effectiveness, and costs. This study endeavors to develop the first machine learning models capable of predicting AMR profiles for the whole repertoire of species and drugs encountered in clinical microbiology. The resulting models can be interpreted as drug recommender systems for infectious diseases. We find that our dual-branch method delivers considerably higher performance compared to previous approaches. In addition, experiments show that the models can be efficiently fine-tuned to data from other clinical laboratories. MALDI-TOF-based AMR recommender systems can, hence, greatly extend the value of MALDI-TOF MS for clinical diagnostics. All code supporting this study is distributed on PyPI and is packaged at https://github.com/gdewael/maldi-nn.
-
- Computational and Systems Biology
- Genetics and Genomics
Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.