Structural distinctions between NAD+ riboswitch domains 1 and 2 determine differential folding and ligand binding.

Riboswitches are important gene regulatory elements frequently encountered in bacterial mRNAs. The recently discovered nadA riboswitch contains two similar, tandemly arrayed aptamer domains, with the first domain possessing high affinity for nicotinamide adenine dinucleotide (NAD+). The second domain which comprises the ribosomal binding site in a putative regulatory helix, however, has withdrawn from detection of ligand-induced structural modulation thus far, and therefore, the identity of the cognate ligand and the regulation mechanism have remained unclear. Here, we report crystal structures of both riboswitch domains, each bound to NAD+. Furthermore, we demonstrate that ligand binding to domain 2 requires significantly higher concentrations of NAD+ (or ADP retaining analogs) compared to domain 1. Using a fluorescence spectroscopic approach, we further shed light on the structural features which are responsible for the different ligand affinities, and describe the Mg2+-dependent, distinct folding and pre-organization of their binding pockets. Finally, we speculate about possible scenarios for nadA RNA gene regulation as a putative two-concentration sensor module for a time-controlled signal that is primed and stalled by the gene regulation machinery at low ligand concentrations (domain 1), and finally triggers repression of translation as soon as high ligand concentrations are reached in the cell (domain 2).

Effects of flanking regions on HDV cotranscriptional folding kinetics.

Hepatitis delta virus (HDV) ribozyme performs the self-cleavage activity through folding to a double pseudoknot structure. The folding of functional RNA structures is often coupled with the transcription process. In this work, we developed a new approach for predicting the cotranscriptional folding kinetics of RNA secondary structures with pseudoknots. We theoretically studied the cotranscriptional folding behavior of the 99-nucleotide (nt) HDV sequence, two upstream flanking sequences, and one downstream flanking sequence. During transcription, the 99-nt HDV can effectively avoid the trap intermediates and quickly fold to the cleavage-active state. It is different from its refolding kinetics, which folds into an intermediate trap state. For all the sequences, the ribozyme regions (from 1 to 73) all fold to the same structure during transcription. However, the existence of the 30-nt upstream flanking sequence can inhibit the ribozyme region folding into the active native state through forming an alternative helix Alt1 with the segments 70-90. The longer upstream flanking sequence of 54 nt itself forms a stable hairpin structure, which sequesters the formation of the Alt1 helix and leads to rapid formation of the cleavage-active structure. Although the 55-nt downstream flanking sequence could invade the already folded active structure during transcription by forming a more stable helix with the ribozyme region, the slow transition rate could keep the structure in the cleavage-active structure to perform the activity.

Evidence That Nucleophile Deprotonation Exceeds Bond Formation in the HDV Ribozyme Transition State.

Steric constraints imposed by the active sites of protein and RNA enzymes pose major challenges to the investigation of structure-function relationships within these systems. As a strategy to circumvent such constraints in the HDV ribozyme, we have synthesized phosphoramidites from propanediol derivatives and incorporated them at the 5'-termini of RNA and DNA oligonucleotides to generate a series of novel substrates with nucleophiles perturbed electronically through geminal fluorination. In nonenzymatic, hydroxide-catalyzed intramolecular transphosphorylation of the DNA substrates, pH-rate profiles revealed that fluorine substitution reduces the maximal rate and the kinetic p Ka, consistent with the expected electron-withdrawing effect. In HDV ribozyme reactions, we observed that the RNA substrates undergo transphosphorylation relatively efficiently, suggesting that the conformational constraints imposed by a ribofuranose ring are not strictly required for ribozyme catalysis. In contrast to the nonenzymatic reactions, however, substrate fluorination modestly increases the ribozyme reaction rate, consistent with a mechanism in which (1) the 2'-hydroxyl nucleophile exists predominantly in its neutral, protonated form in the ground state and (2) the 2'-hydroxyl bears some negative charge in the rate-determining step, consistent with a transition state in which the extent of 2'-OH deprotonation exceeds the extent of P-O bond formation.

Complexity in pH-Dependent Ribozyme Kinetics: Dark pKa Shifts and Wavy Rate-pH Profiles.

Charged bases occur in RNA enzymes, or ribozymes, where they play key roles in catalysis. Cationic bases donate protons and perform electrostatic catalysis, while anionic bases accept protons. We previously published simulations of rate-pH profiles for ribozymes in terms of species plots for the general acid and general base that have been useful for understanding how ribozymes respond to pH. In that study, we did not consider interaction between the general acid and general base or interaction with other species on the RNA. Since that report, diverse small ribozyme classes have been discovered, many of which have charged nucleobases or metal ions in the active site that can either directly interact and participate in catalysis or indirectly interact as "influencers". Herein, we simulate experimental rate-pH profiles in terms of species plots in which reverse protonated charged nucleobases interact. These analyses uncover two surprising features of pH-dependent enzyme kinetics. (1) Cooperativity between the general acid and general base enhances population of the functional forms of a ribozyme and manifests itself as hidden or "dark" pKa shifts, real pKa shifts that accelerate the reaction but are not readily observed by standard experimental approaches, and (2) influencers favorably shift the pKas of proton-transferring nucleobases and manifest themselves as "wavy" rate-pH profiles. We identify parallels with the protein enzyme literature, including reverse protonation and wavelike behavior, while pointing out that RNA is more prone to reverse protonation. The complexities uncovered, which arise from simple pairwise interactions, should aid deconvolution of complex rate-pH profiles for RNA and protein enzymes and suggest veiled catalytic devices for promoting catalysis that can be tested by experiment and calculation.

Development and Characterization of Recombinant Virus Generated from a New World Zika Virus Infectious Clone.

Zika virus (ZIKV; family Flaviviridae, genus Flavivirus) is a rapidly expanding global pathogen that has been associated with severe clinical manifestations, including devastating neurological disease in infants. There are currently no molecular clones of a New World ZIKV available that lack significant attenuation, hindering progress toward understanding determinants of transmission and pathogenesis. Here we report the development and characterization of a novel ZIKV reverse genetics system based on a 2015 isolate from Puerto Rico (PRVABC59). We generated a two-plasmid infectious clone system from which infectious virus was rescued that replicates in human and mosquito cells with growth kinetics representative of wild-type ZIKV. Infectious clone-derived virus initiated infection and transmission rates in Aedes aegypti mosquitoes comparable to those of the primary isolate and displayed similar pathogenesis in AG129 mice. This infectious clone system provides a valuable resource to the research community to explore ZIKV molecular biology, vaccine development, antiviral development, diagnostics, vector competence, and disease pathogenesis. IMPORTANCE: ZIKV is a rapidly spreading mosquito-borne pathogen that has been linked to Guillain-Barré syndrome in adults and congenital microcephaly in developing fetuses and infants. ZIKV can also be sexually transmitted. The viral molecular determinants of any of these phenotypes are not well understood. There is no reverse genetics system available for the current epidemic virus that will allow researchers to study ZIKV immunity, develop novel vaccines, or develop antiviral drugs. Here we provide a novel infectious clone system generated from a recent ZIKV isolated from a patient infected in Puerto Rico. This infectious clone produces virus with in vitro and in vivo characteristics similar to those of the primary isolate, providing a critical tool to study ZIKV infection and disease.

The Third Signal Cytokine Interleukin 12 Rather Than Immune Checkpoint Inhibitors Contributes to the Functional Restoration of Hepatitis D Virus-Specific T Cells.

BACKGROUND: Hepatitis D virus (HDV) infection affects 15-20 million individuals worldwide and causes severely progressive hepatitis. It is unknown to what extent cellular immune responses contribute to liver disease and control of viral replication in HDV infection. METHODS: Immune cell frequencies and phenotypes were determined in 49 HDV-infected patients, 25 individuals with hepatitis B virus monoinfection and 18 healthy controls. T-cell proliferative and cytokine-producing capacities were analyzed by stimulation with overlapping peptides spanning the large HDV antigen. To restore T-cell responses, blocking antibodies (anti-cytotoxic T-lymphocyte antigen 4, anti-programmed death ligand 1) or proinflammatory cytokines (interleukin [IL] 12) were used. RESULTS: Immune cell frequencies and phenotypes did not vary between the groups. Exclusively, the senescence marker CD57 was significantly up-regulated in CD8+ T cells from patients with hepatitis delta. HDV-specific T-cell proliferation and cytokine production were weak and could only partly be rescued by blockade of the programmed death 1 pathway. However, a more robust and consistent increase in HDV-specific CD4+ and CD8+ T-cell responses was evident when the third signal cytokine IL-12 was added, which also affected cytomegalovirus- and Epstein-Barr virus-specific T cells. CONCLUSIONS: This investigation of virus-specific T-cell immunity in patients with HDV infection, the largest to date, revealed premature aging of immune cells and impaired T-cell functionality. This could be restored by blocking inhibitory pathways and, in particular, by supplementing with IL-12.

Hepatitis delta virus: insights into a peculiar pathogen and novel treatment options.

Chronic hepatitis D is the most severe form of viral hepatitis, affecting ∼20 million HBV-infected people worldwide. The causative agent, hepatitis delta virus (HDV), is a unique human pathogen: it is the smallest known virus; it depends on HBV to disseminate its viroid-like RNA; it encodes only one protein (HDAg), which has both structural and regulatory functions; and it replicates using predominantly host proteins. The failure of HBV-specific nucleoside analogues to suppress the HBV helper function, and the limitations of experimental systems to study the HDV life cycle, have impeded the development of HDV-specific drugs. Thus, the only clinical regimen for HDV is IFNα, which shows some efficacy but long-term virological responses are rare. Insights into the receptor-mediated entry of HDV, and the observation that HDV assembly requires farnesyltransferase, have enabled novel therapeutic strategies to be developed. Interference with entry, for example through blockade of the HBV-HDV-specific receptor sodium/taurocholate cotransporting polypeptide NTCP by Myrcludex B, and inhibition of assembly by blockade of farnesyltransferase using lonafarnib or nucleic acid polymers such as REP 2139-Ca, have shown promising results in phase II studies. In this Review, we summarize our knowledge of HDV epidemiology, pathogenesis and molecular biology, with a particular emphasis on possible future developments.

Transition State Features in the Hepatitis Delta Virus Ribozyme Reaction Revealed by Atomic Perturbations.

Endonucleolytic ribozymes constitute a class of non-coding RNAs that catalyze single-strand RNA scission. With crystal structures available for all of the known ribozymes, a major challenge involves relating functional data to the physically observed RNA architecture. In the case of the hepatitis delta virus (HDV) ribozyme, there are three high-resolution crystal structures, the product state of the reaction and two precursor variants, with distinct mechanistic implications. Here, we develop new strategies to probe the structure and catalytic mechanism of a ribozyme. First, we use double-mutant cycles to distinguish differences in functional group proximity implicated by the crystal structures. Second, we use a corrected form of the Brønsted equation to assess the functional significance of general acid catalysis in the system. Our results delineate the functional relevance of atomic interactions inferred from structure, and suggest that the HDV ribozyme transition state resembles the cleavage product in the degree of proton transfer to the leaving group.

HDVDB: a data warehouse for hepatitis delta virus.

Hepatitis Delta Virus (HDV) is an RNA virus and causes delta hepatitis in humans. Although a lot of data is available for HDV, but retrieval of information is a complicated task. Current web database 'HDVDB' provides a comprehensive web-resource for HDV. The database is basically concerned with basic information about HDV and disease caused by this virus, genome structure, pathogenesis, epidemiology, symptoms and prevention, etc. Database also supplies sequence data and bibliographic information about HDV. A tool 'siHDV Predict' to design the effective siRNA molecule to control the activity of HDV, is also integrated in database. It is a user friendly information system available at public domain and provides annotated information about HDV for research scholars, scientists, pharma industry people for further study.

The role of an active site Mg(2+) in HDV ribozyme self-cleavage: insights from QM/MM calculations.

The hepatitis delta virus (HDV) ribozyme is a catalytic RNA motif embedded in the human pathogenic HDV RNA. It catalyzes self-cleavage of its sugar-phosphate backbone with direct participation of the active site cytosine C75. Biochemical and structural data support a general acid role of C75. Here, we used hybrid quantum mechanical/molecular mechanical (QM/MM) calculations to probe the reaction mechanism and changes in Gibbs energy along the ribozyme's reaction pathway with an N3-protonated C75H(+) in the active site, which acts as the general acid, and a partially hydrated Mg(2+) ion with one deprotonated, inner-shell coordinated water molecule that acts as the general base. We followed eight reaction paths with a distinct position and coordination of the catalytically important active site Mg(2+) ion. For six of them, we observed feasible activation barriers ranging from 14.2 to 21.9 kcal mol(-1), indicating that the specific position of the Mg(2+) ion in the active site is predicted to strongly affect the kinetics of self-cleavage. The deprotonation of the U-1(2'-OH) nucleophile and the nucleophilic attack of the resulting U-1(2'-O(-)) on the scissile phosphodiester are found to be separate steps, as deprotonation precedes the nucleophilic attack. This sequential mechanism of the HDV ribozyme differs from the concerted nucleophilic activation and attack suggested for the hairpin ribozyme. We estimate the pKa of the U-1(2'-OH) group to range from 8.8 to 11.2, suggesting that it is lowered by several units from that of a free ribose, comparable to and most likely smaller than the pKa of the solvated active site Mg(2+) ion. Our results thus support the notion that the structure of the HDV ribozyme, and particularly the positioning of the active site Mg(2+) ion, facilitate deprotonation and activation of the 2'-OH nucleophile.

Kinetic partitioning mechanism of HDV ribozyme folding.

RNA folding kinetics is directly tied to RNA biological functions. We introduce here a new approach for predicting the folding kinetics of RNA secondary structure with pseudoknots. This approach is based on our previous established helix-based method for predicting the folding kinetics of RNA secondary structure. In this approach, the transition rates for an elementary step: (1) formation, (2) disruption of a helix stem, and (3) helix formation with concomitant partial melting of an incompatible helix, are calculated with the free energy landscape. The folding kinetics of the Hepatitis delta virus (HDV) ribozyme and the mutated sequences are studied with this method. The folding pathways are identified by recursive searching the states with high net flux-in(out) population starting from the native state. The theory results are in good agreement with that of the experiments. The results indicate that the bi-phasic folding kinetics for the wt HDV sequence is ascribed to the kinetic partitioning mechanism: Part of the population will quickly fold to the native state along the fast pathway, while another part of the population will fold along the slow pathway, in which the population is trapped in a non-native state. Single mutation not only changes the folding rate but also the folding pathway.

Sodium taurocholate cotransporting polypeptide is a functional receptor for human hepatitis B and D virus.

Human hepatitis B virus (HBV) infection and HBV-related diseases remain a major public health problem. Individuals coinfected with its satellite hepatitis D virus (HDV) have more severe disease. Cellular entry of both viruses is mediated by HBV envelope proteins. The pre-S1 domain of the large envelope protein is a key determinant for receptor(s) binding. However, the identity of the receptor(s) is unknown. Here, by using near zero distance photo-cross-linking and tandem affinity purification, we revealed that the receptor-binding region of pre-S1 specifically interacts with sodium taurocholate cotransporting polypeptide (NTCP), a multiple transmembrane transporter predominantly expressed in the liver. Silencing NTCP inhibited HBV and HDV infection, while exogenous NTCP expression rendered nonsusceptible hepatocarcinoma cells susceptible to these viral infections. Moreover, replacing amino acids 157-165 of nonfunctional monkey NTCP with the human counterpart conferred its ability in supporting both viral infections. Our results demonstrate that NTCP is a functional receptor for HBV and HDV.DOI:http://dx.doi.org/10.7554/eLife.00049.001.

Programming a highly structured ribozyme into complex allostery using RNA oligonucleotides.

RNA possesses great potential for expanding the toolbox currently available to synthetic biologists. Here, the modulation of the Hepatitis Delta Virus ribozyme's activity with a series of rationally designed aptamers and effector RNA oligonucleotides is described. The ribozyme was initially fused with an 18-nucleotide hairpin structure that abolished its self-cleaving activity. The binding of a 14-mer oligonucleotide to the hairpin rescued the self-cleavage in a concentration-dependent manner. This modified ribozyme was inserted into the 5' UTR of a reporter gene, and the resulting construct was used to demonstrate that it is possible to modulate the ribozyme activity in cellulo with the oligonucleotide. Subsequently, ribozymes possessing specialized aptamers respecting other logic gates were also successfully designed and found to be functional in vitro. To our knowledge, this is the first example of HDV ribozyme regulation by oligonucleotides, as well as the first allosteric regulation of HDV ribozyme in mammalian cells.

Significantly improved rescue of rabies virus from cDNA plasmids.

The rescue of recombinant rabies virus (RV) from cloned cDNA is an inefficient process because it relies on the de novo formation within cells of functional ribonucleoprotein (RNP) complexes from plasmid-expressed viral-like antigenome RNAs and three helper proteins. In the standard RV reverse genetics systems, bacteriophage T7 RNA polymerase drives the transcription of virus antigenome-like RNAs containing three nonviral G residues at the 5'-end and a correct 3'-end generated by the autocatalytic activity of an 85 nucleotides long hepatitis delta virus antigenomic "core" ribozyme (HDVagrz). Here, we show that employing optimized ribozyme sequences significantly improves RV rescue. Substitution of the "core" HDVagrz by a ribozyme with an enhanced cleavage activity resulted in an approximately 10-fold higher number of rescue events and faster initiation of an infectious cycle. The alternative use of a hammerhead ribozyme for the generation of an exact 5'-end similarly enhanced rescue efficiency. Notably, RV cDNA clones containing the combination of optimized 3'- and 5'-ribozymes were rescued at an at least 100-fold increase. In addition to virus rescue, reporter gene expression from transfected minigenome cDNAs was significantly enhanced by the novel ribozymes. The improved RV reverse genetics system greatly facilitates recovery of strongly attenuated viruses and vectors for biomedical applications.

Developing three-dimensional models of putative-folding intermediates of the HDV ribozyme.

Both the role and the interacting partners of an RNA molecule can change depending on its tertiary structure. Consequently, it is important to be able to accurately predict the complete folding pathway of an RNA molecule. The hepatitis delta virus (HDV) ribozyme is a small catalytic RNA with the greatest number of folding intermediates making it the model of choice with which to address this problem. The tertiary structures of the known putative intermediates along the folding pathway of the HDV ribozyme were predicted using the Macromolecular Conformations Symbolic programming (MC-Sym) software. The structures obtained by this method received physical support from Selective 2'-Hydroxyl Acylation analyzed by Primer Extension (SHAPE). The analysis of these structures elucidated several features of the HDV ribozyme. In addition, this report represents an application for MC-Sym that permits progression one step further toward the computer prediction of an RNA molecule-folding pathway.

Nondenaturing purification of co-transcriptionally folded RNA avoids common folding heterogeneity.

Due to the energetic frustration of RNA folding, tertiary structured RNA is typically characterized by a rugged folding free energy landscape where deep kinetic barriers separate numerous misfolded states from one or more native states. While most in vitro studies of RNA rely on (re)folding chemically and/or enzymatically synthesized RNA in its entirety, which frequently leads into kinetic traps, nature reduces the complexity of the RNA folding problem by segmental, co-transcriptional folding starting from the 5' end. We here have developed a simplified, general, nondenaturing purification protocol for RNA to ask whether avoiding denaturation of a co-transcriptionally folded RNA can reduce commonly observed in vitro folding heterogeneity. Our protocol bypasses the need for large-scale auxiliary protein purification and expensive chromatographic equipment and involves rapid affinity capture with magnetic beads and removal of chemical heterogeneity by cleavage of the target RNA from the beads using the ligand-induced glmS ribozyme. For two disparate model systems, the Varkud satellite (VS) and hepatitis delta virus (HDV) ribozymes, we achieve >95% conformational purity within one hour of enzymatic transcription, without the need for any folding chaperones. We further demonstrate that in vitro refolding introduces severe conformational heterogeneity into the natively-purified VS ribozyme but not into the compact, double-nested pseudoknot fold of the HDV ribozyme. We conclude that conformational heterogeneity in complex RNAs can be avoided by co-transcriptional folding followed by nondenaturing purification, providing rapid access to chemically and conformationally pure RNA for biologically relevant biochemical and biophysical studies.

Intrinsic disorder and oligomerization of the hepatitis delta virus antigen.

The 195 amino acid basic protein (δAg) of hepatitis delta virus (HDV) is essential for replication of the HDV RNA genome. Numerous properties have been mapped to full-length δAg and attempts made to link these to secondary, tertiary and quaternary structures. Here, for the full-size δAg, extensive intrinsic disorder was predicted using PONDR-FIT, a meta-predictor of intrinsic disorder, and evidenced by circular dichroism measurements. Most δAg amino acids are in disordered configurations with no more than 30% adopting an α-helical structure. In addition, dynamic light scattering studies indicated that purified δAg assembled into structures of as large as dodecamers. Cross-linking followed by denaturing polyacrylamide gel electrophoresis revealed hexamers to octamers for this purified δAg and at least this size for δAg found in virus-like particles. Oligomers of purified δAg were resistant to elevated NaCl and urea concentrations, and bound without specificity to RNA and single- and double-stranded DNAs.

Development of mathematical models for the analysis of hepatitis delta virus viral dynamics.

BACKGROUND: Mathematical models have shown to be extremely helpful in understanding the dynamics of different virus diseases, including hepatitis B. Hepatitis D virus (HDV) is a satellite virus of the hepatitis B virus (HBV). In the liver, production of new HDV virions depends on the presence of HBV. There are two ways in which HDV can occur in an individual: co-infection and super-infection. Co-infection occurs when an individual is simultaneously infected by HBV and HDV, while super-infection occurs in persons with an existing chronic HBV infection. METHODOLOGY/PRINCIPAL FINDINGS: In this work a mathematical model based on differential equations is proposed for the viral dynamics of the hepatitis D virus (HDV) across different scenarios. This model takes into consideration the knowledge of the biology of the virus and its interaction with the host. In this work we will present the results of a simulation study where two scenarios were considered, co-infection and super-infection, together with different antiviral therapies. Although, in general the predicted course of HDV infection is similar to that observed for HBV, we observe a faster increase in the number of HBV infected cells and viral load. In most tested scenarios, the number of HDV infected cells and viral load values remain below corresponding predicted values for HBV. CONCLUSIONS/SIGNIFICANCE: The simulation study shows that, under the most commonly used and generally accepted therapy approaches for HDV infection, such as lamivudine (LMV) or ribavirine, peggylated alpha-interferon (IFN) or a combination of both, LMV monotherapy and combination therapy of LMV and IFN were predicted to more effectively reduce the HBV and HDV viral loads in the case of super-infection scenarios when compared with the co-infection. In contrast, IFN monotherapy was found to reduce the HDV viral load more efficiently in the case of super-infection while the effect on the HBV viral load was more pronounced during co-infection. The results suggest that there is a need for development of high efficacy therapeutic approaches towards the specific inhibition of HDV replication. These approaches may additionally be directed to the reduction of the half-life of infected cells and life-span of newly produced circulating virions.