Identification of novel tetrahydroquinoxaline derived phenyl ureas as modulators of the hepatitis B virus nucleocapsid assembly.

A key step of hepatitis B virus (HBV) replication is the selective packaging of pregenomic RNA (pgRNA) by core protein (Cp) dimers, forming a nucleocapsid where the reverse transcriptional viral DNA replication takes place. One approach in the development of new anti-HBV drugs is to disrupt the assembly of HBV nucleocapsids by misdirecting Cp dimers to assemble morphologically normal capsids devoid of pgRNA. In this study, we built upon our previous discovery of benzamide-derived HBV capsid assembly modulators by exploring fused bicyclic scaffolds with an exocyclic amide that is ß, γ to the fused ring, and identified 1,2,3,4-tetrahydroquinoxaline derived phenyl ureas as a novel scaffold. Structure-activity relationship studies showed that a favorable hydrophobic substitution can be tolerated at the 2-position of the 1,2,3,4-tetrahydroquinoxaline core, and the resulting compound 88 demonstrated comparable or improved antiviral potencies in mouse and human hepatocyte-derived HBV-replicating cell lines compared to our previously reported benzamide compound, 38017 (8). In addition, a novel bis-urea series based on 1,2,3,4-tetrahydroquinoxaline was also found to inhibit HBV DNA replication with sub-micromolar EC50 values. The mode of action of these compounds is consistent with specific inhibition of pgRNA encapsidation into nucleocapsids in hepatocytes.

BMaps: A Web Application for Fragment-Based Drug Design and Compound Binding Evaluation.

Fragment-based drug design uses data about where, and how strongly, small chemical fragments bind to proteins, to assemble new drug molecules. Over the past decade, we have been successfully using fragment data, derived from thermodynamically rigorous Monte Carlo fragment-protein binding simulations, in dozens of preclinical drug programs. However, this approach has not been available to the broader research community because of the cost and complexity of doing simulations and using design tools. We have developed a web application, called BMaps, to make fragment-based drug design widely available with greatly simplified user interfaces. BMaps provides access to a large repository (>550) of proteins with 100s of precomputed fragment maps, druggable hot spots, and high-quality water maps. Users can also employ their own structures or those from the Protein Data Bank and AlphaFold DB. Multigigabyte data sets are searched to find fragments in bondable orientations, ranked by a binding-free energy metric. The designers use this to select modifications that improve affinity and other properties. BMaps is unique in combining conventional tools such as docking and energy minimization with fragment-based design, in a very easy to use and automated web application. The service is available at https://www.boltzmannmaps.com.

Structure-based design of selective, orally available salt-inducible kinase inhibitors that stimulate bone formation in mice.

Osteoporosis is a major public health problem. Currently, there are no orally available therapies that increase bone formation. Intermittent parathyroid hormone (PTH) stimulates bone formation through a signal transduction pathway that involves inhibition of salt-inducible kinase isoforms 2 and 3 (SIK2 and SIK3). Here, we further validate SIK2/SIK3 as osteoporosis drug targets by demonstrating that ubiquitous deletion of these genes in adult mice increases bone formation without extraskeletal toxicities. Previous efforts to target these kinases to stimulate bone formation have been limited by lack of pharmacologically acceptable, specific, orally available SIK2/SIK3 inhibitors. Here, we used structure-based drug design followed by iterative medicinal chemistry to identify SK-124 as a lead compound that potently inhibits SIK2 and SIK3. SK-124 inhibits SIK2 and SIK3 with single-digit nanomolar potency in vitro and in cell-based target engagement assays and shows acceptable kinome selectivity and oral bioavailability. SK-124 reduces SIK2/SIK3 substrate phosphorylation levels in human and mouse cultured bone cells and regulates gene expression patterns in a PTH-like manner. Once-daily oral SK-124 treatment for 3 wk in mice led to PTH-like effects on mineral metabolism including increased blood levels of calcium and 1,25-vitamin D and suppressed endogenous PTH levels. Furthermore, SK-124 treatment increased bone formation by osteoblasts and boosted trabecular bone mass without evidence of short-term toxicity. Taken together, these findings demonstrate PTH-like effects in bone and mineral metabolism upon in vivo treatment with orally available SIK2/SIK3 inhibitor SK-124.

4-Oxooctahydroquinoline-1(2H)-carboxamides as hepatitis B virus (HBV) capsid core protein assembly modulators.

Hepatitis B virus (HBV) core protein, the building block of the HBV capsid, plays multiple roles in viral replication, and is an attractive target for development of antiviral agents with a new mechanism of action. In addition to the heteroaryldihydropyrimidines (HAPs), sulfamoylbenzamides (SBAs), dibenzothiazepine derivatives (DBTs), and sulfamoylpyrrolamides (SPAs) that inhibit HBV replication by modulation of viral capsid assembly and are currently under clinical trials for the treatment of chronic hepatitis B (CHB), other chemical structures with activity to modulate HBV capsid assembly have also been explored. Here we describe our continued optimization of a benzamide originating from our high throughput screening. A new bicyclic carboxamide lead featuring an electron deficient non-planar core structure was discovered. Evaluations of its ADMET (absorption, distribution, metabolism, excretion and toxicity) and pharmacokinetic (PK) profiles demonstrate improved metabolic stability and good bioavailability.

Amino acid residues at core protein dimer-dimer interface modulate multiple steps of hepatitis B virus replication and HBeAg biogenesis.

The core protein (Cp) of hepatitis B virus (HBV) assembles pregenomic RNA (pgRNA) and viral DNA polymerase to form nucleocapsids where the reverse transcriptional viral DNA replication takes place. Core protein allosteric modulators (CpAMs) inhibit HBV replication by binding to a hydrophobic "HAP" pocket at Cp dimer-dimer interfaces to misdirect the assembly of Cp dimers into aberrant or morphologically "normal" capsids devoid of pgRNA. We report herein that a panel of CpAM-resistant Cp with single amino acid substitution of residues at the dimer-dimer interface not only disrupted pgRNA packaging, but also compromised nucleocapsid envelopment, virion infectivity and covalently closed circular (ccc) DNA biosynthesis. Interestingly, these mutations also significantly reduced the secretion of HBeAg. Biochemical analysis revealed that the CpAM-resistant mutations in the context of precore protein (p25) did not affect the levels of p22 produced by signal peptidase removal of N-terminal 19 amino acid residues, but significantly reduced p17, which is produced by furin cleavage of C-terminal arginine-rich domain of p22 and secreted as HBeAg. Interestingly, p22 existed as both unphosphorylated and phosphorylated forms. While the unphosphorylated p22 is in the membranous secretary organelles and the precursor of HBeAg, p22 in the cytosol and nuclei is hyperphosphorylated at the C-terminal arginine-rich domain and interacts with Cp to disrupt capsid assembly and viral DNA replication. The results thus indicate that in addition to nucleocapsid assembly, interaction of Cp at dimer-dimer interface also plays important roles in the production and infectivity of progeny virions through modulation of nucleocapsid envelopment and uncoating. Similar interaction at reduced p17 dimer-dimer interface appears to be important for its metabolic stability and sensitivity to CpAM suppression of HBeAg secretion.

Identification of hepatitis B virus core protein residues critical for capsid assembly, pgRNA encapsidation and resistance to capsid assembly modulators.

Assembly of hepatitis B virus (HBV) capsids is driven by the hydrophobic interaction of core protein (Cp) at dimer-dimer interface. Binding of core protein allosteric modulators (CpAMs) to a hydrophobic "HAP" pocket formed between the inter-dimer interface strengths the dimer-dimer interaction and misdirects the assembly of Cp dimers into non-capsid Cp polymers or morphologically normal capsids devoid of viral pregenomic (pg) RNA and DNA polymerase. In this study, we performed a systematic mutagenesis analysis to identify Cp amino acid residues at Cp dimer-dimer interface that are critical for capsid assembly, pgRNA encapsidation and resistance to CpAMs. By analyzing 70 mutant Cp with a single amino acid substitution of 25 amino acid residues around the HAP pocket, our study revealed that residue W102 and Y132 are critical for capsid assembly. However, substitution of many other residues did not significantly alter the amount of capsids, but reduced the amount of encapsidated pgRNA, suggesting their critical roles in pgRNA packaging. Interestingly, several mutant Cp with a single amino acid substitution of residue P25, T33 or I105 supported high levels of DNA replication, but conferred strong resistance to multiple chemotypes of CpAMs. In addition, we also found that WT Cp, but not the assembly incompetent Cp, such as Y132A Cp, interacted with HBV DNA polymerase (Pol). This later finding implies that encapsidation of viral DNA polymerase may depend on the interaction of Pol with a capsid assembly intermediate, but not free Cp dimers. Taking together, our findings reported herein shed new light on the mechanism of HBV nucleocapsid assembly and mode of CpAM action.

A putative amphipathic alpha helix in hepatitis B virus small envelope protein plays a critical role in the morphogenesis of subviral particles.

Hepatitis B virus (HBV) small (S) envelope protein has the intrinsic ability to direct the formation of small spherical subviral particles (SVPs) in eukaryotic cells. However, the molecular mechanism underlying the morphogenesis of SVPs from the monomeric S protein initially synthesized at the endoplasmic reticulum (ER) membrane remains largely elusive. Structure prediction and extensive mutagenesis analysis suggested that the amino acid residues spanning W156 to R169 of S protein form an amphipathic alpha helix and play essential roles in SVP production and S protein metabolic stability. Further biochemical analyses showed that the putative amphipathic alpha helix was not required for the disulfide-linked S protein oligomerization, but was essential for SVP morphogenesis. Pharmacological disruption of vesicle trafficking between the ER and Golgi complex in SVP producing cells supported the hypothesis that S protein-directed SVP morphogenesis takes place at the ER-Golgi intermediate compartment (ERGIC). Moreover, it was demonstrated that S protein is degraded in hepatocytes via a 20S proteasome-dependent, but ubiquitination-independent non-classic ER-associated degradation (ERAD) pathway. Taken together, the results reported herein favor a model in which the amphipathic alpha helix at the antigenic loop of S protein attaches to the lumen leaflet to facilitate SVP budding from the ERGIC compartment, whereas the failure of budding process may result in S protein degradation by 20S proteasome in an ubiquitination-independent manner.Importance Subviral particles are the predominant viral product produced by HBV-infected hepatocytes. Their levels exceed the virion particles by 10,000 to 100,000-fold in the blood of HBV infected individuals. The high levels of SVPs, or HBV surface antigen (HBsAg), in the circulation induces immune tolerance and contributes to the establishment of persistent HBV infection. The loss of HBsAg, often accompanied by appearance of anti-HBs antibodies, is the hallmark of durable immune control of HBV infection. Therapeutic induction of HBsAg loss is, therefore, considered to be essential for the restoration of host antiviral immune response and functional cure of chronic hepatitis B. Our findings on the mechanism of SVP morphogenesis and S protein metabolism will facilitate the rational discovery and development of antiviral drugs to achieve this therapeutic goal.

Targeting the multifunctional HBV core protein as a potential cure for chronic hepatitis B.

The core (capsid) protein of hepatitis B virus (HBV) is the building block of nucleocapsids where viral DNA reverse transcriptional replication takes place and mediates virus-host cell interaction important for the persistence of HBV infection. The pleiotropic role of core protein (Cp) in HBV replication makes it an attractive target for antiviral therapies of chronic hepatitis B, a disease that affects more than 257 million people worldwide without a cure. Recent clinical studies indicate that core protein allosteric modulators (CpAMs) have a great promise as a key component of hepatitis B curative therapies. Particularly, it has been demonstrated that modulation of Cp dimer-dimer interactions by several chemical series of CpAMs not only inhibit nucleocapsid assembly and viral DNA replication, but also induce the disassembly of double-stranded DNA-containing nucleocapsids to prevent the synthesis of cccDNA. Moreover, the different chemotypes of CpAMs modulate Cp assembly by interaction with distinct amino acid residues at the HAP pocket between Cp dimer-dimer interfaces, which results in the assembly of Cp dimers into either non-capsid Cp polymers (type I CpAMs) or empty capsids with distinct physical property (type II CpAMs). The different CpAMs also differentially modulate Cp metabolism and subcellular distribution, which may impact cccDNA metabolism and host antiviral immune responses, the critical factors for the cure of chronic HBV infection. This review article highlights the recent research progress on the structure and function of core protein in HBV replication cycle, the mode of action of CpAMs, as well as the current status and perspectives on the discovery and development of core protein-targeting antivirals. This article forms part of a symposium in Antiviral Research on "Wide-ranging immune and direct-acting antiviral approaches to curing HBV and HDV infections."

Profiling the circulating mRNA transcriptome in human liver disease.

The human circulation contains cell-free DNA and non-coding microRNA (miRNA). Less is known about the presence of messenger RNA (mRNA). This report profiles the human circulating mRNA transcriptome in people with liver cirrhosis (LC) and hepatocellular carcinoma (HCC) to determine whether mRNA analytes can be used as biomarkers of liver disease. Using RNAseq and RT-qPCR, we investigate circulating mRNA in plasma from HCC and LC patients and demonstrate detection of transcripts representing more than 19,000 different protein coding genes. Remarkably, the circulating mRNA expression levels were similar from person to person over the 21 individuals whose samples were analyzed by RNAseq. Liver derived circulating transcripts such as albumin (ALB), apolipoprotein (APO) A1, A2 & H, serpin A1 & E1, ferritin light chain (FTL) and fibrinogen like 1 (FGL1) were significantly upregulated in HCC patient samples. Higher levels of some of these liver-specific transcripts in the plasma of HCC patients were confirmed by RT-qPCR in another cohort of 20 individuals. Several less abundant circulating transcripts associated with cancer were detected in most HCC samples, but not in healthy subjects. Liver specificity of circulating transcripts was confirmed by investigating their expression in HCC tumor and liver cancer cell lines. Liver specific mRNA sequences in the plasma were predominantly present outside circulating extracellular vesicles. Conclusions: The circulating "mRNA" transcriptome is remarkably consistent in diversity and expression from person to person. Detection of transcripts corresponding to disease selective polypeptides suggests the possibility that circulating mRNA can work as a biomarker analyte for cancer detection.

Discovery and Mechanistic Study of a Novel Human-Stimulator-of-Interferon-Genes Agonist.

Stimulator of interferon genes (STING) is an integral ER-membrane protein that can be activated by 2'3'-cGAMP synthesized by cyclic guanosine monophosphate-adenosine monophosphate synthase (cGAS) upon binding of double-stranded DNA. It activates interferon (IFN) and inflammatory cytokine responses to defend against infection by microorganisms. Pharmacologic activation of STING has been demonstrated to induce an antiviral state and boost antitumor immunity. We previously reported a cell-based high-throughput-screening assay that allowed for identification of small-molecule cGAS-STING-pathway agonists. We report herein a compound, 6-bromo-N-(naphthalen-1-yl)benzo[d][1,3]dioxole-5-carboxamide (BNBC), that induces a proinflammatory cytokine response in a human-STING-dependent manner. Specifically, we showed that BNBC induced type I and III IFN dominant cytokine responses in primary human fibroblasts and peripheral-blood mononuclear cells (PBMCs). BNBC also induced cytokine response in PBMC-derived myeloid dendritic cells and promoted their maturation, suggesting that STING-agonist treatment could potentially regulate the activation of CD4+ and CD8+ T lymphocytes. As anticipated, treatment of primary human fibroblast cells with BNBC induced an antiviral state that inhibited the infection of several kinds of flaviviruses. Taken together, our results indicate that BNBC is a human-STING agonist that not only induces innate antiviral immunity against a broad spectrum of viruses but may also stimulate the activation of adaptive immune responses, which is important for the treatment of chronic viral infections and tumors.

CpAMs induce assembly of HBV capsids with altered electrophoresis mobility: Implications for mechanism of inhibiting pgRNA packaging.

Native agarose gel electrophoresis-based particle gel assay has been commonly used for examination of hepatitis B virus (HBV) capsid assembly and pregenomic RNA encapsidation in HBV replicating cells. Interestingly, treatment of cells with several chemotypes of HBV core protein allosteric modulators (CpAMs) induced the assembly of both empty and DNA-containing capsids with faster electrophoresis mobility. In an effort to determine the physical basis of CpAM-induced capsid mobility shift, we found that the surface charge, but not the size, of capsids is the primary determinant of electrophoresis mobility. Specifically, through alanine scanning mutagenesis analysis of twenty-seven charged amino acids in core protein assembly domain and hinge region, we showed that except for K7 and E8, substitution of glutamine acid (E) or aspartic acid (D) on the surface of capsids reduced their mobility, but substitution of lysine (K) or arginine (R) on the surface of capsids increased their mobility in variable degrees. However, alanine substitution of the charged amino acids that are not exposed on the surface of capsid did not apparently alter capsid mobility. Hence, CpAM-induced electrophoresis mobility shift of capsids may reflect the global alteration of capsid structure that changes the exposure and/or ionization of charged amino acid side chains of core protein. Our findings imply that CpAM inhibition of pgRNA encapsidation is possibly due to the assembly of structurally altered nucleocapsids. Practically, capsid electrophoresis mobility shift is a diagnostic marker of compounds that target core protein assembly and predicts sensitivity of HBV strains to specific CpAMs.

Identification of a novel inhibitor of dengue virus protease through use of a virtual screening drug discovery Web portal.

We report the discovery of a novel small-molecule inhibitor of the dengue virus (DENV) protease (NS2B-NS3pro) using a newly constructed Web-based portal (DrugDiscovery@TACC) for structure-based virtual screening. Our drug discovery portal, an extension of virtual screening studies performed using IBM's World Community Grid, facilitated access to supercomputer resources managed by the Texas Advanced Computing Center (TACC) and enabled druglike commercially available small-molecule libraries to be rapidly screened against several high-resolution DENV NS2B-NS3pro crystallographic structures. Detailed analysis of virtual screening docking scores and hydrogen-bonding interactions between each docked ligand and the NS2B-NS3pro Ser135 side chain were used to select molecules for experimental validation. Compounds were ordered from established chemical companies, and compounds with established aqueous solubility were tested for their ability to inhibit DENV NS2B-NS3pro cleavage of a model substrate in kinetic studies. As a proof-of-concept, we validated a small-molecule dihydronaphthalenone hit as a single-digit-micromolar mixed noncompetitive inhibitor of the DENV protease. Since the dihydronaphthalenone was predicted to interact with NS2B-NS3pro residues that are largely conserved between DENV and the related West Nile virus (WNV), we tested this inhibitor against WNV NS2B-NS3pro and observed a similar mixed noncompetitive inhibition mechanism. However, the inhibition constants were ∼10-fold larger against the WNV protease relative to the DENV protease. This novel validated lead had no chemical features or pharmacophores associated with adverse toxicity, carcinogenicity, or mutagenicity risks and thus is attractive for additional characterization and optimization.

Modeling nonaqueous proton wires built from helical peptides: biased proton transfer driven by helical dipoles.

We report gas-phase electronic structure calculations on helical peptides that act as scaffolds for imidazole-based hydrogen-bonding networks (proton wires). We have modeled various 21-residue polyalanine peptides substituted at regular intervals with histidines (imidazole-bearing amino acids), using a hybrid approach with a semiempirical method (AM1) for peptide scaffolds and density functional theory (B3LYP) for proton wires. We have computed energy landscapes including barriers for Grotthuss-shuttling-type proton motions though wires supported on 3(10)-, α- and π-helical structures, showing the 3(10)- and α-helices to be attractive targets in terms of high proton affinities, low Grotthuss shuttling barriers, and high stabilities. Moreover, bias forces provided by the helical dipole moments were found to promote unidirectional proton translocation.

Modeling energy landscapes of proton motion in nonaqueous, tethered proton wires.

We have modeled structures and energetics of anhydrous proton-conducting wires: tethered hydrogen-bonded chains of the form ···HX···HX···HX···, with functional groups HX = imidazole, triazole, and formamidine; formic, sulfonic, and phosphonic acids. We have applied density functional theory (DFT) to model proton wires up to 19 units long, where each proton carrier is linked to an effective backbone to mimic polymer tethering. This approach allows the direct calculation of hydrogen bond strengths. The proton wires were found to be stabilized by strong hydrogen bonds (up to 50 kJ/mol) whose strength correlates with the proton affinity of HX [related to pK(b)(HX)] and not to pK(a)(HX) as is often assumed. Geometry optimizations and ab initio molecular dynamics near 400 K on imidazole-based proton wires both predict that adding a proton to the end of such wires causes the excess charge to embed into the interior segments of these wires. Proton translocation energy landscapes for imidazole-based wires are sensitive to the imidazole attachment point (head or feet) and to wire architecture (linear or interdigitated). Linear imidazole wires with head-attachment exhibit low barriers for intrawire proton motion, rivaling proton diffusion in liquid imidazole. Excess charge relaxation from the edge of wires is found to be dominated by long-range Grotthuss shuttling for distances as long as 42 Å, especially for interdigitated wires. For imidazole, we predict that proton translocation is controlled by the energetics of desorption from the proton wire, even for relatively long wires (600 imidazole units). Proton desorption energies show no correlation with functional group properties, suggesting that proton desorption is a collective process in proton wires.

Importance of dynamic hydrogen bonds and reorientation barriers in proton transport.

The dynamic nature of hydrogen bonds in phenolic polymers, where the hydrogen bond donor/acceptor reorientation can occur in a single site, presents lower barriers for proton transport.