Mosaic amino acid conservation in 3D-structures of surface protein and polymerase of hepatitis B virus

Hepatitis B virus envelope proteins can serve as therapeutic targets embedded in the host cell plasma membrane

Hepatitis B virus (HBV) infection is a major health threat causing 880,000 deaths each year. Available therapies control viral replication but do not cure HBV, leaving patients at risk to develop hepatocellular carcinoma. Here, we show that HBV envelope proteins (HBs)-besides their integration into endosomal membranes-become embedded in the plasma membrane where they can be targeted by redirected T-cells. HBs was detected on the surface of HBV-infected cells, in livers of mice replicating HBV and in HBV-induced hepatocellular carcinoma. Staining with HBs-specific recombinant antibody MoMab recognising a conformational epitope indicated that membrane-associated HBs remains correctly folded in HBV-replicating cells in cell culture and in livers of HBV-transgenic mice in vivo. MoMab coated onto superparamagnetic iron oxide nanoparticles allowed to detect membrane-associated HBs after HBV infection by electron microscopy in distinct stretches of the hepatocyte plasma membrane. Last but not least, we demonstrate that HBs located on the cell surface allow therapeutic targeting of HBV-positive cells by T-cells either engrafted with a chimeric antigen receptor or redirected by bispecific, T-cell engager antibodies. TAKE AWAYS: HBs become translocated to the plasma membrane. Novel, recombinant antibody confirmed proper conformation of HBs on the membrane. HBs provide an interesting target by T-cell-based, potentially curative therapies.

A fusion peptide in preS1 and the human protein disulfide isomerase ERp57 are involved in hepatitis B virus membrane fusion process

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.

Key mutations in the C-terminus of the HBV surface glycoprotein correlate with lower HBsAg levels in vivo, hinder HBsAg secretion in vitro and reduce HBsAg structural stability in the setting of HBeAg-negative chronic HBV genotype-D infection

Increasing evidences suggest that HBsAg-production varies across HBV-genotypes. HBsAg C-terminus plays a crucial role for HBsAg-secretion. Here, we evaluate HBsAg-levels in different HBV-genotypes in HBeAg-negative chronic infection, the correlation of specific mutations in HBsAg C-terminus with HBsAg-levels in-vivo, their impact on HBsAg-secretion in-vitro and on structural stability in-silico.

HBsAg-levels were investigated in 323 drug-naïve HBeAg-negative patients chronically infected with HBV genotype-D(N = 228), -A(N = 65) and -E(N = 30). Genotype-D was characterized by HBsAg-levels lower than genotype-A and -E (3.3[2.7–3.8]IU/ml; 3.8[3.5–4.2]IU/ml and 3.9[3.7–4.2]IU/ml, P < 0.001). Results confirmed by multivariable analysis correcting for patients’demographics, HBV-DNA, ALT and infection-status.

In genotype-D, specific C-terminus mutations (V190A-S204N-Y206C-Y206F-S210N) significantly correlate with HBsAg<1000IU/ml(P-value from <0.001 to 0.04). These mutations lie in divergent pathways involving other HBsAg C-terminus mutations: V190A + F220L (Phi = 0.41, P = 0.003), S204N + L205P (Phi = 0.36, P = 0.005), Y206F + S210R (Phi = 0.47, P < 0.001) and S210N + F220L (Phi = 0.40, P = 0.006). Notably, patients with these mutational pairs present HBsAg-levels 1log lower than patients without them(P-value from 0.003 to 0.02). In-vitro, the above-mentioned mutational pairs determined a significant decrease in HBsAg secretion-efficiency compared to wt(P-value from <0.001 to 0.02). Structurally, these mutational pairs reduced HBsAg C-terminus stability and determined a rearrangement of this domain.

In conclusion, HBsAg-levels in genotype-D are significantly lower than in genotype-A and -E in HBeAg-negative patients. In genotype-D, specific mutational clusters in HBsAg C-terminus correlate with lower HBsAg-levels in-vivo, hamper HBsAg-release in-vitro and affect its structural stability, supporting their detrimental role on HBsAg-secretion. In this light, genotypic-testing can be a valuable tool to optimize the clinical interpretation of HBsAg in genotype-D and to provide information on HBV-pathogenicity and disease-progression.

HDV Can Constrain HBV Genetic Evolution in HBsAg: Implications for the Identification of Innovative Pharmacological Targets

Chronic HBV + HDV infection is associated with greater risk of liver fibrosis, earlier hepatic decompensation, and liver cirrhosis hepatocellular carcinoma compared to HBV mono-infection. However, to-date no direct anti-HDV drugs are available in clinical practice. Here, we identified conserved and variable regions in HBsAg and HDAg domains in HBV + HDV infection, a critical finding for the design of innovative therapeutic agents. The extent of amino-acid variability was measured by Shannon-Entropy (Sn) in HBsAg genotype-d sequences from 31 HBV + HDV infected and 62 HBV mono-infected patients (comparable for demographics and virological-parameters), and in 47 HDAg genotype-1 sequences. Positions with Sn = 0 were defined as conserved. The percentage of conserved HBsAg-positions was significantly higher in HBV + HDV infection than HBV mono-infection (p = 0.001). Results were confirmed after stratification for HBeAg-status and patients’ age. A Sn = 0 at specific positions in the C-terminus HBsAg were correlated with higher HDV-RNA, suggesting that conservation of these positions can preserve HDV-fitness. Conversely, HDAg was characterized by a lower percentage of conserved-residues than HBsAg (p < 0.001), indicating higher functional plasticity. Furthermore, specific HDAg-mutations were significantly correlated with higher HDV-RNA, suggesting a role in conferring HDV replicative-advantage. Among HDAg-domains, only the virus-assembly signal exhibited a high genetic conservation (75% of conserved-residues). In conclusion, HDV can constrain HBsAg genetic evolution to preserve its fitness. The identification of conserved regions in HDAg poses the basis for designing innovative targets against HDV-infection.

Bioinformatic Identification of Rare Codon Clusters (RCCs) in HBV Genome and Evaluation of RCCs in Proteins Structure of Hepatitis B Virus

BACKGROUND

Hepatitis B virus (HBV) as an infectious disease that has nine genotypes (A - I) and a 'putative' genotype J.

OBJECTIVES

The aim of this study was to identify the rare codon clusters (RCC) in the HBV genome and to evaluate these RCCs in the HBV proteins structure.

METHODS

For detection of protein family accession numbers (Pfam) in HBV proteins, the UniProt database and Pfam search tool were used. Protein family accession numbers is a comprehensive and accurate collection of protein domains and families. It contains annotation of each family in the form of textual descriptions, links to other resources and literature references. Genome projects have used Pfam extensively for large-scale functional annotation of genomic data; Pfam database is a large collection of protein families, each represented by multiple sequence alignments and hidden Markov models (HMMs). The Pfam search tools are databases that identify Pfam of proteins. These Pfam IDs were analyzed in Sherlocc program and the location of RCCs in HBV genome and proteins were detected and reported as translated EMBL nucleotide sequence data library (TrEMBL) entries. The TrEMBL is a computer-annotated supplement of SWISS-PROT that contains all the translations of European molecular biology laboratory (EMBL) nucleotide sequence entries not yet integrated in SWISS-PROT. Furthermore, the structures of TrEMBL entries proteins were studied in the PDB database and 3D structures of the HBV proteins and locations of RCCs were visualized and studied using Swiss PDB Viewer software®.

RESULTS

The Pfam search tool found nine protein families in three frames. Results of Pfams studies in the Sherlocc program showed that this program has not identified RCCs in the external core antigen (PF08290) and truncated HBeAg gene (PF08290) of HBV. By contrast, the RCCs were identified in gene of hepatitis core antigen (PF00906 and the residues 224 - 234 and 251 - 255), large envelope protein S (PF00695 and the residues 53-56 and 70 - 84), X protein (PF00739 and the residues 10 - 24, 29 - 83, 95 - 99. 122 - 129, 139 - 143), DNA polymerase (viral) N-terminal domain (PF00242 and the residues 59 - 62, 214 - 217, 407 - 413) and protein P (Pf00336 and the residues 225 - 228). In HBV genome, seven RCCs were identified in the gene area of hepatitis core antigen, large envelope protein S and DNA polymerase, while protein structures of TrEMBL entries sequences found in Sherlocc program outputs were not complete.

CONCLUSIONS

Based on the location of detected RCCs in the structure of HBV proteins, it was found that these RCCs may have a critical role in correct folding of HBV proteins and can be considered as drug targets. The results of this study provide new and deep perspectives about structure of HBV proteins for further researches and designing new drugs for treatment of HBV.

BAsE-Seq: a method for obtaining long viral haplotypes from short sequence reads

We present a method for obtaining long haplotypes, of over 3 kb in length, using a short-read sequencer, Barcode-directed Assembly for Extra-long Sequences (BAsE-Seq). BAsE-Seq relies on transposing a template-specific barcode onto random segments of the template molecule and assembling the barcoded short reads into complete haplotypes. We applied BAsE-Seq on mixed clones of hepatitis B virus and accurately identified haplotypes occurring at frequencies greater than or equal to 0.4%, with >99.9% specificity. Applying BAsE-Seq to a clinical sample, we obtained over 9,000 viral haplotypes, which provided an unprecedented view of hepatitis B virus population structure during chronic infection. BAsE-Seq is readily applicable for monitoring quasispecies evolution in viral diseases.

BAsE-Seq: a method for obtaining long viral haplotypes from short sequence reads

We present a method for obtaining long haplotypes, of over 3 kb in length, using a short-read sequencer, Barcode-directed Assembly for Extra-long Sequences (BAsE-Seq). BAsE-Seq relies on transposing a template-specific barcode onto random segments of the template molecule and assembling the barcoded short reads into complete haplotypes. We applied BAsE-Seq on mixed clones of hepatitis B virus and accurately identified haplotypes occurring at frequencies greater than or equal to 0.4%, with >99.9% specificity. Applying BAsE-Seq to a clinical sample, we obtained over 9,000 viral haplotypes, which provided an unprecedented view of hepatitis B virus population structure during chronic infection. BAsE-Seq is readily applicable for monitoring quasispecies evolution in viral diseases.

Specific mutations in the C-terminus domain of HBV surface antigen significantly correlate with low level of serum HBV-DNA in patients with chronic HBV infection

BACKGROUND

To define HBsAg-mutations correlated with different serum HBV-DNA levels in HBV chronically-infected drug-naive patients.

METHODS

This study included 187 patients stratified into the following ranges of serum HBV-DNA:12-2000 IU/ml, 2000-100,000 IU/ml, and >100,000 IU/ml. HBsAg-mutations were associated with HBV-DNA levels by applying a Bayesian-Partitional-Model and Fisher-exact test. Mutant and wild-type HBV genotype-D genomes were expressed in Huh7 cells and HBsAg-production was determined in cell-supernatants at 3 days-post-transfection.

RESULTS

Specific HBsAg-mutations (M197T,-S204N-Y206C/H-F220L) were significantly correlated with serum HBV-DNA <2000 IU/ml (posterior-probability>90%, P < 0.05). The presence of Y206C/H and/or F220L was also associated with lower median (IQR) HBsAg-levels and lower median (IQR) transaminases (for HBsAg:250[115-840] IU/ml for Y206C/H and/or F220L versus 4300[640-11,838] IU/ml for wild-type, P = 0.023; for ALT:28[21-40] IU/ml versus 53[34-90] IU/ml, P < 0.001). These mutations were localized in the HBsAg C-terminus, known to be involved in virion and/or HBsAg secretion. The co-occurrence of Y206C + F220L was found significant by cluster-analysis, (P = 0.02). In addition, in an in-vitro model Y206C + F220L determined a 2.8-3.3 fold-reduction of HBsAg-amount released in supernatants compared to single mutants and wt (Y206C + F220L = 5,679 IU/ml; Y206H = 16,305 IU/ml; F220L = 18,368 IU/ml; Y206C = 18,680 IU/ml; wt = 14,280 IU/ml, P < 0.05).

CONCLUSIONS

Specific HBsAg-mutations (compartmentalized in the HBsAg C-terminus) correlated with low-serum HBV-DNA and HBsAg-levels. These findings can be important to understand mechanisms underlying low HBV replicative potential including the inactive-carrier state.

Solution structure and in silico binding of a cyclic peptide with hepatitis B surface antigen

A specific ligand targeting the immunodominant region of hepatitis B virus is desired in neutralizing the infectivity of the virus. In a previous study, a disulfide constrained cyclic peptide cyclo S(1) ,S(9) Cys-Glu-Thr-Gly-Ala-Lys-Pro-His-Cys (S(1) , S(9) -cyclo-CETGAKPHC) was isolated from a phage displayed cyclic peptide library using an affinity selection method against hepatitis B surface antigen. The cyclic peptide binds tightly to hepatitis B surface antigen with a relative dissociation constant (KD (rel) ) of 2.9 nm. The binding site of the peptide was located at the immunodominant region on hepatitis B surface antigen. Consequently, this study was aimed to elucidate the structure of the cyclic peptide and its interaction with hepatitis B surface antigen in silico. The solution structure of this cyclic peptide was solved using (1) H, (13) C, and (15) N NMR spectroscopy and molecular dynamics simulations with NMR-derived distance and torsion angle restraints. The cyclic peptide adopted two distinct conformations due to the isomerization of the Pro residue with one structured region in the ETGA sequence. Docking studies of the peptide ensemble with a model structure of hepatitis B surface antigen revealed that the cyclic peptide can potentially be developed as a therapeutic drug that inhibits the virus-host interactions.

Quasispecies structure, cornerstone of hepatitis B virus infection: Mass sequencing approach

Hepatitis B virus (HBV) is a DNA virus with complex replication, and high replication and mutation rates, leading to a heterogeneous viral population. The population is comprised of genomes that are closely related, but not identical; hence, HBV is considered a viral quasispecies. Quasispecies variability may be somewhat limited by the high degree of overlapping between the HBV coding regions, which is especially important in the P and S gene overlapping regions, but is less significant in the X and preCore/Core genes. Despite this restriction, several clinically and pathologically relevant variants have been characterized along the viral genome. Next-generation sequencing (NGS) approaches enable high-throughput analysis of thousands of clonally amplified regions and are powerful tools for characterizing genetic diversity in viral strains. In the present review, we update the information regarding HBV variability and present a summary of the various NGS approaches available for research in this virus. In addition, we provide an analysis of the clinical implications of HBV variants and their study by NGS.

The 'Red Queen' dilemma--running to stay in the same place: reflections on the evolutionary vector of HBV in humans

BACKGROUND

Estimates for the evolutionary rate of HBV until now have been interpreted as showing that HBV is a relatively recent acquisition for mankind. The existence of defined HBV genotypes is thought to represent past founder effects. We have explored virus mutation in a group of 48 persistently infected blood donors sampled at two points in time and suggest otherwise.

METHODS

HBV-infected donors were detected by hepatitis B surface antigen (HBsAg) screening and staged by hepatitis B e markers. Serum DNA from those persistently infected with HBV was characterized by consensus sequencing and the amino acid sequences inferred. These were compared against consensus genotype sequences and divergence measured at two points in time.

RESULTS

Rates of viral mutation were higher across both HBsAg and hepatitis B core antigen in the group of donors seropositive for hepatitis B e antibody (1.36×10⁻³ and 1.54×10⁻³ changes per residue per year, respectively) than in those seropositive for hepatitis B e antigen (4.59×10⁻⁴ and 6.62×10⁻⁴ changes per residue per year, respectively). Codon mutations reverting to the genotype consensus were commonly seen. Codon changes were clustered close to the C-terminal region of HBsAg and were accommodated in overlapping polymerase by synonymous substitutions.

CONCLUSIONS

It is suggested that in vivo HBV behaves as a self-normalizing meme and mutational rates, although high, do not lead to significant change over time in a persistent infection. This would be compatible with co-evolution within its human host and introduction within humans being an ancient occurrence.

Ultra-deep pyrosequencing detects conserved genomic sites and quantifies linkage of drug-resistant amino acid changes in the hepatitis B virus genome

BACKGROUND

Selection of amino acid substitutions associated with resistance to nucleos(t)ide-analog (NA) therapy in the hepatitis B virus (HBV) reverse transcriptase (RT) and their combination in a single viral genome complicates treatment of chronic HBV infection and may affect the overlapping surface coding region. In this study, the variability of an overlapping polymerase-surface region, critical for NA resistance, is investigated before treatment and under antiviral therapy, with assessment of NA-resistant amino acid changes simultaneously occurring in the same genome (linkage analysis) and their influence on the surface coding region.

METHODOLOGY/PRINCIPAL FINDINGS

Serum samples obtained from chronic HBV-infected patients at pre-treatment and during sequential NA treatment with lamivudine, adefovir, and entecavir were analyzed by ultra-deep pyrosequencing (UDPS) using the GS-FLX platform (454 Life Sciences-Roche). The pre-treatment HBV quasispecies was not enriched with NA-resistant substitutions. The frequencies of this type of substitutions at pre-treatment did not predict the frequencies observed during lamivudine treatment. On linkage analysis of the RT region studied, NA-resistant HBV variants (except for rtA181T) were present in combinations of amino acid substitutions that increased in complexity after viral breakthrough to entecavir, at which time the combined variant rtL180M-S202G-M204V-V207I predominated. In the overlapping surface region, NA-resistant substitutions caused selection of stop codons in a significant percentage of sequences both at pre-treatment and during sequential treatment; the rtA181T substitution, related to sW172stop, predominated during treatment with lamivudine and adefovir. A highly conserved RT residue (rtL155), even more conserved than the essential residues in the RT catalytic motif YMDD, was identified in all samples.

CONCLUSIONS

UDPS methodology enabled quantification of HBV quasispecies variants, even those harboring complex combinations of amino acid changes. The high percentage of potentially defective genomes, especially in the surface region, suggests effective trans-complementation of these variants.

Protein X of hepatitis B virus: origin and structure similarity with the central domain of DNA glycosylase

Orthohepadnavirus (mammalian hosts) and avihepadnavirus (avian hosts) constitute the family of Hepadnaviridae and differ by their capability and inability for expression of protein X, respectively. Origin and functions of X are unclear. The evolutionary analysis at issue of X indicates that present strains of orthohepadnavirus started to diverge about 25,000 years ago, simultaneously with the onset of avihepadnavirus diversification. These evolutionary events were preceded by a much longer period during which orthohepadnavirus developed a functional protein X while avihepadnavirus evolved without X. An in silico generated 3D-model of orthohepadnaviral X protein displayed considerable similarity to the tertiary structure of DNA glycosylases (key enzymes of base excision DNA repair pathways). Similarity is confined to the central domain of MUG proteins with the typical DNA-binding facilities but without the capability of DNA glycosylase enzymatic activity. The hypothetical translation product of a vestigial X reading frame in the genome of duck hepadnavirus could also been folded into a DNA glycosylase-like 3D-structure. In conclusion, the most recent common ancestor of ortho- and avihepadnavirus carried an X sequence with orthology to the central domain of DNA glycosylase.

The genetic backbone modulates the phenotype of hepatitis B surface antigen mutants

Hepatitis B virus (HBV) vaccine and diagnostic escape mutants are a growing concern. The principle target of detection, hepatitis B surface antigen (HBsAg), encoded by S, is completely overlapped by the reverse transcriptase encoding P. With the increased incidence of nucleos(t)ide analogue resistance altering P, the concurrent impact on S must be assessed. HBV DNA from 59 HBsAg-positive plasma samples was sequenced across the polymerase/surface region and the amino acid sequence of HBsAg was inferred. ELISAs were formatted containing individually bound monoclonal antibodies directed against three discrete epitopes on HBsAg. Similar point mutations occurring in different genotypes were shown to influence epitope conformation differently, indicating that the genetic backbone is a major factor in predicting phenotype. C-terminal changes associated with antiviral resistance were found to modulate epitope profiles of HBsAg. Treatment options which may promote drug resistance should be avoided to both protect antiviral treatment and prevent facilitation of vaccine and diagnostic escape mutants.

A function essential to viral entry underlies the hepatitis B virus "a" determinant

The hepatitis B virus (HBV) particles bear a receptor-binding site located in the pre-S1 domain of the large HBV envelope protein. Using the hepatitis delta virus (HDV) as a surrogate of HBV, a second infectivity determinant was recently identified in the envelope proteins antigenic loop (AGL), and its activity was shown to depend upon cysteine residues that are essential for the structure of the HBV immunodominant "a" determinant. Here, an alanine-scanning mutagenesis approach was used to precisely map the AGL infectivity determinant to a set of conserved residues, which are predicted to cluster together with cysteines in the AGL disulfide bridges network. Several substitutions suppressed both infectivity and the "a" determinant, whereas others were infectivity deficient with only a partial impact on antigenicity. Interestingly, G145R, a substitution often arising under immune pressure selection and detrimental to the "a" determinant, had no effect on infectivity. Altogether, these findings indicate that the AGL infectivity determinant is closely related to, yet separable from, the "a" determinant. Finally, a selection of HDV entry-deficient mutations were introduced at the surface of HBV virions and shown to also abrogate infection in the HBV model. Therefore, a function can at last be assigned to the orphan "a" determinant, the first-discovered marker of HBV infection. The characterization of the AGL functions at viral entry may lead to novel approaches in the development of antivirals against HBV.

Occult hepatitis B infection: an evolutionary scenario

Background

Occult or latent hepatitis B virus (HBV) infection is defined as infection with detectable HBV DNA and undetectable surface antigen (HBsAg) in patients' blood. The cause of an overt HBV infection becoming an occult one is unknown. To gain insight into the mechanism of the development of occult infection, we compared the full-length HBV genome from a blood donor carrying an occult infection (d4) with global genotype D genomes.

Results

The phylogenetic analysis of polymerase, core and X protein sequences did not distinguish d4 from other genotype D strains. Yet, d4 surface protein formed the evolutionary outgroup relative to all other genotype D strains. Its evolutionary branch was the only one where accumulation of substitutions suggests positive selection (dN/dS = 1.3787). Many of these substitutiions accumulated specifically in regions encoding the core/surface protein interface, as revealed in a 3D-modeled protein complex. We identified a novel RNA splicing event (deleting nucleotides 2986-202) that abolishes surface protein gene expression without affecting polymerase, core and X-protein related functions. Genotype D strains differ in their ability to perform this 2986-202 splicing. Strains prone to 2986-202 splicing constitute a separate clade in a phylogenetic tree of genotype D HBVs. A single substitution (G173T) that is associated with clade membership alters the local RNA secondary structure and is proposed to affect splicing efficiency at the 202 acceptor site.

Conclusion

We propose an evolutionary scenario for occult HBV infection, in which 2986-202 splicing generates intracellular virus particles devoid of surface protein, which subsequently accumulates mutations due to relaxation of coding constraints. Such viruses are deficient of autonomous propagation and cannot leave the host cell until it is lysed.