Using All-Atom Potentials to Refine RNA Structure Predictions of SARS-CoV-2 Stem Loops.

A considerable amount of rapid-paced research is underway to combat the SARS-CoV-2 pandemic. In this work, we assess the 3D structure of the 5' untranslated region of its RNA, in the hopes that stable secondary structures can be targeted, interrupted, or otherwise measured. To this end, we have combined molecular dynamics simulations with previous Nuclear Magnetic Resonance measurements for stem loop 2 of SARS-CoV-1 to refine 3D structure predictions of that stem loop. We find that relatively short sampling times allow for loop rearrangement from predicted structures determined in absence of water or ions, to structures better aligned with experimental data. We then use molecular dynamics to predict the refined structure of the transcription regulatory leader sequence (TRS-L) region which includes stem loop 3, and show that arrangement of the loop around exchangeable monovalent potassium can interpret the conformational equilibrium determined by in-cell dimethyl sulfate (DMS) data.

Mechanistic analysis of the hepatitis delta virus (HDV) ribozyme: methods for RNA preparation, structure mapping, solvent isotope effects, and co-transcriptional cleavage.

Small ribozymes such as the hairpin, hammerhead, VS, glm S, and hepatitis delta virus (HDV) are self-cleaving RNAs that are typically characterized by kinetics and structural methods. Working with these RNAs requires attention to numerous experimental details. In this chapter we focus on four different experimental aspects of ribozyme studies: preparing the RNA, mapping its structure with reverse transcription and end-labeled techniques, solvent isotope experiments, and co-transcriptional cleavage assays. Although the focus of these methods is the HDV ribozyme, the methods should be applicable to other ribozymes.

Structural characterization of the viral and cRNA panhandle motifs from the infectious salmon anemia virus.

Infectious salmon anemia virus (ISAV) has emerged as a virus of great concern to the aquaculture industry since it can lead to highly contagious and lethal infections in farm-raised salmon populations. While little is known about the transcription/replication cycle of ISAV, initial evidence suggests that it follows molecular mechanisms similar to those found in other orthomyxoviruses, which include the highly pathogenic influenza A (inf A) virus. During the life cycle of orthomyxoviruses, a panhandle structure is formed by the pairing of the conserved 5' and 3' ends of each genomic RNA. This structural motif serves both as a promoter of the viral RNA (vRNA)-dependent RNA polymerase and as a regulatory element in the transcription/replication cycle. As a first step toward characterizing the structure of the ISAV panhandle, here we have determined the secondary structures of the vRNA and the cRNA panhandles on the basis of solution nuclear magnetic resonance (NMR) and thermal melting data. The vRNA panhandle is distinguished by three noncanonical U · G pairs and one U · U pair in two stem helices that are linked by a highly stacked internal loop. For the cRNA panhandle, a contiguous stem helix with a protonated C · A pair near the terminus and tandem downstream U · U pairs was found. The observed noncanonical base pairs and base stacking features of the ISAV RNA panhandle motif provide the first insight into structural features that may govern recognition by the viral RNA polymerase.

Mechanistic characterization of the HDV genomic ribozyme: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid.

The hepatitis delta virus (HDV) ribozyme uses the nucleobase C75 and a hydrated Mg(2+) ion as the general acid-base catalysts in phosphodiester bond cleavage at physiological salt. A mechanistic framework has been advanced that involves one Mg(2+)-independent and two Mg(2+)-dependent channels. The rate-pH profile for wild-type (WT) ribozyme in the Mg(2+)-free channel is inverted relative to the fully Mg(2+)-dependent channel, with each having a near-neutral pKa. Inversion of the rate-pH profile was used as the crux of a mechanistic argument that C75 serves as general acid both in the presence and absence of Mg(2+). However, subsequent studies on a double mutant (DM) ribozyme suggested that the pKa observed for WT in the absence of Mg(2+) arises from ionization of C41, a structural nucleobase. To investigate this further, we acquired rate-pH/pD profiles and proton inventories for WT and DM in the absence of Mg(2+). Corrections were made for effects of ionic strength on hydrogen ion activity and pH meter readings. Results are accommodated by a model wherein the Mg(2+)-free pKa observed for WT arises from ionization of C75, and DM reactivity is compromised by protonation of C41. The Brønsted base appears to be water or hydroxide ion depending on pH. The observed pKa's are related to salt-dependent pH titrations of a model oligonucleotide, as well as electrostatic calculations, which support the local environment for C75 in the absence of Mg(2+) being similar to that in the presence of Mg(2+) and impervious to bulk ions. Accordingly, the catalytic role of C75 as the general acid does not appear to depend on divalent ions or the identity of the Brønsted base.

Mechanistic characterization of the HDV genomic ribozyme: the cleavage site base pair plays a structural role in facilitating catalysis.

The hepatitis delta virus (HDV) ribozyme occurs in the genomic and antigenomic strands of the HDV RNA and within mammalian transcriptomes. Previous kinetic studies suggested that a wobble pair (G*U or A(+)*C) is preferred at the cleavage site; however, the reasons for this are unclear. We conducted sequence comparisons, which indicated that while G*U is the most prevalent combination at the cleavage site, G-C occurs to a significant extent in genomic HDV isolates, and G*U, G-C, and A-U pairs are present in mammalian ribozymes. We analyzed the folding of genomic HDV ribozymes by free energy minimization and found that variants with purine-pyrimidine combinations at the cleavage site are predicted to form native structures while pyrimidine-purine combinations misfold, consistent with earlier kinetic data and sequence comparisons. To test whether the cleavage site base pair contributes to catalysis, we characterized the pH and Mg(2+)-dependence of reaction kinetics of fast-folding genomic HDV ribozymes with cleavage site base pair purine-pyrimidine combinations: G*U, A-U, G-C, and A(+)*C. Rates for these native-folding ribozymes displayed highly similar pH and Mg(2+) concentration dependencies, with the exception of the A(+)*C ribozyme, which deviated at high pH. None of the four ribozymes underwent miscleavage. These observations support the A(+)*C ribozyme as being more active with a wobble pair at the cleavage site than with no base pair at all. Overall, the data support a model in which the cleavage site base pair provides a structural role in catalysis and does not need to be a wobble pair.

Wild-type is the optimal sequence of the HDV ribozyme under cotranscriptional conditions.

RNA viruses are responsible for a variety of human diseases, and the pathogenicity of RNA viruses is often attributed to a high rate of mutation. Self-cleavage activity of the wild-type hepatitis delta virus (HDV) ribozyme as measured in standard divalent ion renaturation assays is biphasic and mostly slow and can be improved by multiple rational changes to ribozyme sequence or by addition of chemical denaturants. This is unusual in the sense that wild type is the most catalytically active sequence for the majority of protein enzymes, and RNA viruses are highly mutable. To see whether the ribozyme takes advantage of fast-reacting sequence changes in vivo, we performed alignment of 76 genomic and 269 antigenomic HDV isolates. Paradoxically, the sequence for the ribozyme was found to be essentially invariant in nature. We therefore tested whether three ribozyme sequence changes that improve self-cleavage under standard divalent ion renaturation assays also improve self-cleavage during transcription. Remarkably, wild type was as fast, or faster, than these mutants under cotranscriptional conditions. Slowing the rate of transcription or adding the hepatitis delta antigen protein only further stimulated cotranscriptional self-cleavage activity. Thus, the relative activity of HDV ribozyme mutants depends critically on whether the reaction is assayed under in vivo-like conditions. A model is presented for how wild-type ribozyme sequence and flanking sequence work in concert to promote efficient self-cleavage during transcription. Wild type being the optimal ribozyme sequence under in vivo-like conditions parallels the behavior of most protein enzymes.

Insight into the functional versatility of RNA through model-making with applications to data fitting.

The RNA World hypothesis posits that life emerged from self-replicating RNA molecules. For any single biopolymer to be the basis for life, it must both store information and perform diverse functions. It is well known that RNA can store information. Advances in recent years have revealed that RNA can exhibit remarkable functional versatility as well. In an effort to judge the functional versatility of RNA and thereby the plausibility that RNA was at one point the basis for life, a statistical chemical approach is adopted. Essential biological functions are reduced to simple molecular models in a minimalist, biopolymer-independent fashion. The models dictate requisite states, populations of states, and physical and chemical changes occurring between the states. Next, equations are derived from the models, which lead to complex phenomenological constants such as observed and functional constants that are defined in terms of familiar elementary chemical descriptors: intrinsic rate constants, microscopic ligand equilibrium constants, secondary structure stability, and ligand concentration. Using these equations, simulations of functional behavior are performed. These functional models provide practical frameworks for fitting and organizing real data on functional RNAs such as ribozymes and riboswitches. At the same time, the models allow the suitability of RNA as a basis for life to be judged. We conclude that RNA, while inferior to extant proteins in most, but not all, functional respects, may be more versatile than proteins, performing a wider range of elementary biological functions at a tolerable level. Inspection of the functional models and various RNA structures uncovers several surprising ways in which the nucleobases can conspire to afford chemical catalysis and evolvability. These models support the plausibility that RNA, or a closely related informational biopolymer, could serve as the basis for a fairly simple form of life.