Molecular mechanism of capsid disassembly in hepatitis B virus

Mechanisms of Hepatitis B Virus cccDNA and Minichromosome Formation and HBV Gene Transcription

Hepatitis B virus (HBV) is the etiologic agent of chronic hepatitis B, which puts at least 300 million patients at risk of developing fibrosis, cirrhosis, and hepatocellular carcinoma. HBV is a partially double-stranded DNA virus of the Hepadnaviridae family. While HBV was discovered more than 50 years ago, many aspects of its replicative cycle remain incompletely understood. Central to HBV persistence is the formation of covalently closed circular DNA (cccDNA) from the incoming relaxed circular DNA (rcDNA) genome. cccDNA persists as a chromatinized minichromosome and is the major template for HBV gene transcription. Here, we review how cccDNA and the viral minichromosome are formed and how viral gene transcription is regulated and highlight open questions in this area of research.

Impact of DNA on interactions between core proteins of Hepatitis B virus-like particles comprising different C-terminals

Hepatitis B virus (HBV) virus-like particles (VLPs) are promising therapeutic agents derived from HBV core proteins (Cp). This study investigates the assembly dynamics of HBV VLPs, which is crucial for their potential as drug carriers or gene delivery systems. Coarse-grained molecular dynamics simulations explore the impact of C-terminal domain length (in the Cp ranging from Cp149 to wild-type Cp183) on Cp assembly and stability, particularly in the presence of DNA. Our findings reveal that the C-terminal nucleic acid binding region significantly influences Cp assembly and stability of trimers comprising Cp dimers. Shorter C-terminal domains (Cp164, Cp167) enhance stability and protein-protein interactions, while interactions between naturally occurring Cp183 are destabilized in the absence of DNA. Interestingly, DNA addition further stabilizes Cp assemblies, and this effect is influenced by the length of the nucleic acid binding region. Shorter C-terminal domains show less dependency on DNA content. This stabilization is attributed to electrostatic forces between positively charged C-terminal chains and negatively charged nucleic acids. Our study sheds light on the molecular mechanisms governing protein-protein and protein-DNA interactions in HBV VLP assembly, providing insights into Cp processability and informing the development of efficient gene therapy carriers using VLP technology.

Computational approaches to investigate fluoride binding, selectivity and transport across the membrane

The use of molecular dynamics (MD) simulations to study biomolecular systems has proven reliable in elucidating atomic-level details of structure and function. In this chapter, MD simulations were used to uncover new insights into two phylogenetically unrelated bacterial fluoride (F-) exporters: the CLCF F-/H+ antiporter and the Fluc F- channel. The CLCF antiporter, a member of the broader CLC family, has previously revealed unique stoichiometry, anion-coordinating residues, and the absence of an internal glutamate crucial for proton import in the CLCs. Through MD simulations enhanced with umbrella sampling, we provide insights into the energetics and mechanism of the CLCF transport process, including its selectivity for F- over HF. In contrast, the Fluc F- channel presents a novel architecture as a dual topology dimer, featuring two pores for F- export and a central non-transported sodium ion. Using computational electrophysiology, we simulate the electrochemical gradient necessary for F- export in Fluc and reveal details about the coordination and hydration of both F- and the central sodium ion. The procedures described here delineate the specifics of these advanced techniques and can also be adapted to investigate other membrane protein systems.

A computational spatial whole-Cell model for hepatitis B viral infection and drug interactions

Despite a vaccine, hepatitis B virus (HBV) remains a world-wide source of infections and deaths. We develop a whole-cell computational platform combining spatial and kinetic models describing the infection cycle of HBV in a hepatocyte host. We simulate key parts of the infection cycle with this whole-cell platform for 10 min of biological time, to predict infection progression, map out virus-host and virus-drug interactions. We find that starting from an established infection, decreasing the copy number of the viral envelope proteins shifts the dominant infection pathway from capsid secretion to re-importing the capsids into the nucleus, resulting in more nuclear-localized viral covalently closed circular DNA (cccDNA) and boosting transcription. This scenario can mimic the consequence of drugs designed to manipulate viral gene expression. Mutating capsid proteins facilitates capsid destabilization and disassembly at nuclear pore complexes, resulting in an increase in cccDNA copy number. However, excessive destabilization leads to premature cytoplasmic disassembly and does not increase the cccDNA counts. Finally, our simulations can predict the best drug dosage and its administration timing to reduce the cccDNA counts. Our adaptable computational platform can be parameterized to study other viruses and identify the most central viral pathways that can be targeted by drugs.

Study of self-assembling properties of HBc-VLP derivatives aided by molecular dynamic simulations from a thermodynamic perspective

Self-assembling protein nanoparticles showed promise for vaccine design due to efficient antigen presentations and safety. However, the unpredictable formations of epitopes-fused protein assemblies remain challenging in the upstream design. This study suggests employing molecular dynamic (MD) simulations to investigate the assembly properties of Hepatitis B core protein (HBc) from thermodynamic perspectives. Eight HBc derivatives were expressed in E. coli, with their self-assembly properties characterised by high-performance liquid chromatography and transmission electron microscopy. MD simulations on the dimers, based on AlphaFold-predicted 3D structures, analysed the derivative at the atomic level. Results revealed that HBc derivatives can form dissociative polymers or large multi-subunit structures due to assembly failures. The instability of the dimer in aqueous solvents or inappropriate intradimer distances could cause major assembly failures. Polar solvation energies played a vital role too in forming assemble-incompetent dimers. Importantly, our study demonstrated that MD simulations on dimers can provide preliminary predictions on the assembly properties of HBc derivatives, thus aiding vaccine design by lowering the risk of self-assembling failures in engineered proteins.Communicated by Ramaswamy H. Sarma.

Impact of core protein naturally selected mutants associated with HBeAg-negative status in HBV biosynthesis

Hepatitis B e antigen (HBeAg) loss represents a late stage of chronic hepatitis B virus (HBV) infection associated with a drastic decrease in HBV-DNA, a lower risk of disease progression, and the occurrence of several mutations in the preCore/core region. However, the underlying mechanisms supporting the downregulation of viral replication have yet to be elucidated. In the present study, the analysis of the frequency of subgenotype D1 core protein (HBc) mutations associated with HBeAg status revealed a higher mutation rate in HBeAg-negative sequences compared to HBeAg-positive ones. Particularly, 22 amino acids exhibited a higher frequency of mutation in HBeAg-negative sequences, while the remaining residues showed a high degree of conservation. Subsequently, the assessment of HBc mutants derived from HBeAg-negative patients in viral structure and replicative capacity revealed that HBc mutations have the ability to modulate the subcellular localization of the protein (either when the protein was expressed alone or in the context of viral replication), capsid assembly, and, depending on specific mutation patterns, alter covalently closed circular DNA (cccDNA) recycling and up- or downregulate viral replication. In conclusion, HBc mutations associated with HBeAg-negative status impact on various stages of the HBV life cycle modulating viral replication during the HBeAg-negative stage of infection.

Relaxational dynamics of the T-number conversion of virus capsids

We extend a recently proposed kinetic theory of virus capsid assembly based on Model A kinetics and study the dynamics of the interconversion of virus capsids of different sizes triggered by a quench, that is, by sudden changes in the solution conditions. The work is inspired by in vitro experiments on functionalized coat proteins of the plant virus cowpea chlorotic mottle virus, which undergo a reversible transition between two different shell sizes (T = 1 and T = 3) upon changing the acidity and salinity of the solution. We find that the relaxation dynamics are governed by two time scales that, in almost all cases, can be identified as two distinct processes. Initially, the monomers and one of the two types of capsids respond to the quench. Subsequently, the monomer concentration remains essentially constant, and the conversion between the two capsid species completes. In the intermediate stages, a long-lived metastable steady state may present itself, where the thermodynamically less stable species predominate. We conclude that a Model A based relaxational model can reasonably describe the early and intermediate stages of the conversion experiments. However, it fails to provide a good representation of the time evolution of the state of assembly of the coat proteins in the very late stages of equilibration when one of the two species disappears from the solution. It appears that explicitly incorporating the nucleation barriers to assembly and disassembly is crucial for an accurate description of the experimental findings, at least under conditions where these barriers are sufficiently large.

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

Multiscale Modeling of Hepatitis B Virus Capsid Assembly and Its Dimorphism

Hepatitis B virus (HBV) is an endemic, chronic virus that leads to 800000 deaths per year. Central to the HBV lifecycle, the viral core has a protein capsid assembled from many copies of a single protein. The capsid protein adopts different (quasi-equivalent) conformations to form icosahedral capsids containing 180 or 240 proteins: T = 3 or T = 4, respectively, in Caspar-Klug nomenclature. HBV capsid assembly has become an important target for recently developed antivirals; nonetheless, the assembly pathways and mechanisms that control HBV dimorphism remain unclear. We describe computer simulations of the HBV assembly, using a coarse-grained model that has parameters learned from all-atom molecular dynamics simulations of a complete HBV capsid and yet is computationally tractable. Dynamical simulations with the resulting model reproduce experimental observations of HBV assembly pathways and products. By constructing Markov state models and employing transition path theory, we identify pathways leading to T = 3, T = 4, and other experimentally observed capsid morphologies. The analysis shows that capsid polymorphism is promoted by the low HBV capsid bending modulus, where the key factors controlling polymorphism are the conformational energy landscape and protein-protein binding affinities.

Region-Specific Hepatitis B Virus Genome Exposure from Nucleocapsid Modulated by Capsid Linker Sequence and Inhibitor: Implications for Uncoating

Hepatitis B virus (HBV) contains a partially double-stranded, relaxed circular (RC) DNA genome synthesized within a nucleocapsid (NC) in the host cell cytoplasm. The release of RC DNA from the NC, in an ill-defined process called uncoating, to the nucleus is required for its conversion to the covalently closed circular (CCC) DNA, the viral episome serving as the transcriptional template for all viral RNAs necessary for replication and, thus, essential for establishing and sustaining viral infection. In efforts to better understand uncoating, we analyzed HBV core (HBc) mutants that show various levels of nuclear CCC DNA but little to no cytoplasmic RC DNA. We found that RC DNA could be synthesized by these mutants outside the cell, but in contrast to the wild type (wt), the mutant NCs were unable to protect RC DNA from digestion by the endogenous nuclease(s) in cellular lysates or exogenous DNase. Subcellular fractionation suggested that the major RC DNA-degrading activity was membrane associated. Digestion with sequence-specific and nonspecific DNases revealed the exposure of specific regions of RC DNA from the mutant NC. Similarly, treatment of wt NCs with a core inhibitor known to increase CCC DNA by affecting uncoating also led to region-specific exposure of RC DNA. Furthermore, a subpopulation of untreated wild type (wt) mature NCs showed site-specific exposure of RC DNA as well. Competition between RC DNA degradation and its conversion to CCC DNA during NC uncoating thus likely plays an important role in the establishment and persistence of HBV infection and has implications for the development of capsid-targeted antivirals. IMPORTANCE Disassembly of the hepatitis B virus (HBV) nucleocapsid (NC) to release its genomic DNA, in an ill-understood process called uncoating, is required to form the viral nuclear episome in the host cell nucleus, a viral DNA essential for establishing and sustaining HBV infection. The elimination of the HBV nuclear episome remains the holy grail for the development of an HBV cure. We report here that the HBV genomic DNA is exposed in a region-specific manner during uncoating, which is enhanced by mutations of the capsid protein and a capsid-targeted antiviral compound. The exposure of the viral genome can result in its rapid degradation or, alternatively, can enhance the formation of the nuclear episome, thus having a major impact on HBV infection and persistence. These results are thus important for understanding fundamental mechanisms of HBV replication and persistence and for the ongoing pursuit of an HBV cure.

The Mechanism of Action of Hepatitis B Virus Capsid Assembly Modulators Can Be Predicted from Binding to Early Assembly Intermediates

Interfering with the self-assembly of virus nucleocapsids is a promising approach for the development of novel antiviral agents. Applied to hepatitis B virus (HBV), this approach has led to several classes of capsid assembly modulators (CAMs) that target the virus by either accelerating nucleocapsid assembly or misdirecting it into noncapsid-like particles, thereby inhibiting the HBV replication cycle. Here, we have assessed the structures of early nucleocapsid assembly intermediates, bound with and without CAMs, using molecular dynamics simulations. We find that distinct conformations of the intermediates are induced depending on whether the bound CAM accelerates or misdirects assembly. Specifically, the assembly intermediates with bound misdirecting CAMs appear to be flattened relative to those with bound accelerators. Finally, the potency of CAMs within the same class was studied. We find that an increased number of contacts with the capsid protein and favorable binding energies inferred from free energy perturbation calculations are indicative of increased potency.