SHAMS: combining chemical modification of RNA with mass spectrometry to examine polypurine tract-containing RNA/DNA hybrids

RNA Chemical Labeling with Site-Specific, Relative Quantification by Mass Spectrometry for the Structural Study of a Neomycin-Sensing Riboswitch Aptamer Domain

High-resolution mass spectrometry was used for the label-free, direct localization and relative quantification of CMC+ -modifications of a neomycin-sensing riboswitch aptamer domain in the absence and presence of the aminoglycoside ligands neomycin B, ribostamycin, and paromomycin. The chemical probing and MS data for the free riboswitch show high exposure to solvent of the uridine nucleobases U7, U8, U13, U14, U18 as part of the proposed internal and apical loops, but those of U10 and U21 as part of the proposed internal loop were found to be far less exposed than expected. Thus, our data are in better agreement with the proposed secondary structure of the riboswitch in complexes with aminoglycosides than with that of free RNA. For the riboswitch in complexes with neomycin B, ribostamycin, and paromomycin, we found highly similar CMC+ -modification patterns and excellent agreement with previous NMR studies. Differences between the chemical probing and MS data in the absence and presence of the aminoglycoside ligands were quantitative rather than qualitative (i. e., the same nucleobases were labeled, but to different extents) and can be rationalized by stabilization of both the proposed bulge and the apical loop by aminoglycoside binding. Our study shows that chemical probing and mass spectrometry can provide important structural information and complement other techniques such as NMR spectroscopy.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Computational modeling of RNA 3D structure based on experimental data

RNA molecules are master regulators of cells. They are involved in a variety of molecular processes: they transmit genetic information, sense cellular signals and communicate responses, and even catalyze chemical reactions. As in the case of proteins, RNA function is dictated by its structure and by its ability to adopt different conformations, which in turn is encoded in the sequence. Experimental determination of high-resolution RNA structures is both laborious and difficult, and therefore the majority of known RNAs remain structurally uncharacterized. To address this problem, predictive computational methods were developed based on the accumulated knowledge of RNA structures determined so far, the physical basis of the RNA folding, and taking into account evolutionary considerations, such as conservation of functionally important motifs. However, all theoretical methods suffer from various limitations, and they are generally unable to accurately predict structures for RNA sequences longer than 100-nt residues unless aided by additional experimental data. In this article, we review experimental methods that can generate data usable by computational methods, as well as computational approaches for RNA structure prediction that can utilize data from experimental analyses. We outline methods and data types that can be potentially useful for RNA 3D structure modeling but are not commonly used by the existing software, suggesting directions for future development.

New windows into retroviral RNA structures

Background

The multiple roles of both viral and cellular RNAs have become increasingly apparent in recent years, and techniques to model them have become significantly more powerful, enabling faster and more accurate visualization of RNA structures.

Main body

Techniques such as SHAPE (selective 2’OH acylation analysed by primer extension) have revolutionized the field, and have been used to examine RNAs belonging to many and diverse retroviruses. Secondary structure probing reagents such as these have been aided by the development of faster methods of analysis either via capillary or next-generation sequencing, allowing the analysis of entire genomes, and of retroviral RNA structures within virions. Techniques to model the three-dimensional structures of these large RNAs have also recently developed.

Conclusions

The flexibility of retroviral RNAs, both structural and functional, is clear from the results of these new experimental techniques. Retroviral RNA structures and structural changes control many stages of the lifecycle, and both the RNA structures themselves and their interactions with ligands are potential new drug targets. In addition, our growing understanding of retroviral RNA structures is aiding our knowledge of cellular RNA form and function.

Structure-specific ribonucleases for MS-based elucidation of higher-order RNA structure

Supported by high-throughput sequencing technologies, structure-specific nucleases are experiencing a renaissance as biochemical probes for genome-wide mapping of nucleic acid structure. This report explores the benefits and pitfalls of the application of Mung bean (Mb) and V1 nuclease, which attack specifically single- and double-stranded regions of nucleic acids, as possible structural probes to be employed in combination with MS detection. Both enzymes were found capable of operating in ammonium-based solutions that are preferred for high-resolution analysis by direct infusion electrospray ionization (ESI). Sequence analysis by tandem mass spectrometry (MS/MS) was performed to confirm mapping assignments and to resolve possible ambiguities arising from the concomitant formation of isobaric products with identical base composition and different sequences. The observed products grouped together into ladder-type series that facilitated their assignment to unique regions of the substrate, but revealed also a certain level of uncertainty in identifying the boundaries between paired and unpaired regions. Various experimental factors that are known to stabilize nucleic acid structure, such as higher ionic strength, presence of Mg(II), etc., increased the accuracy of cleavage information, but did not completely eliminate deviations from expected results. These observations suggest extreme caution in interpreting the results afforded by these types of reagents. Regardless of the analytical platform of choice, the results highlighted the need to repeat probing experiments under the most diverse possible conditions to recognize potential artifacts and to increase the level of confidence in the observed structural information.

Ty3 reverse transcriptase complexed with an RNA-DNA hybrid shows structural and functional asymmetry

Retrotransposons are a class of mobile genetic elements that replicate by converting their single-stranded RNA intermediate to double-stranded DNA through the combined DNA polymerase and ribonuclease H (RNase H) activities of the element-encoded reverse transcriptase (RT). Although a wealth of structural information is available for lentiviral and gammaretroviral RTs, equivalent studies on counterpart enzymes of long terminal repeat (LTR)-containing retrotransposons, from which they are evolutionarily derived, is lacking. In this study, we report the first crystal structure of a complex of RT from the Saccharomyces cerevisiae LTR retrotransposon Ty3 in the presence of its polypurine tract-containing RNA-DNA hybrid. In contrast to its retroviral counterparts, Ty3 RT adopts an asymmetric homodimeric architecture whose assembly is substrate dependent. Moreover, our structure and biochemical data suggest that the RNase H and DNA polymerase activities are contributed by individual subunits of the homodimer.

A polymorphism at position 400 in the connection subdomain of HIV-1 reverse transcriptase affects sensitivity to NNRTIs and RNaseH activity

Reverse transcriptase (RT) plays an essential role in HIV-1 replication, and inhibition of this enzyme is a key component of HIV-treatment. However, the use of RT inhibitors can lead to the emergence of drug-resistant variants. Until recently, most clinically relevant resistance mutations were found in the polymerase domain of RT. Lately, an increasing number of resistance mutations has been identified in the connection and RNaseH domain. To further explore the role of these domains we analyzed the complete RT sequence of HIV-1 subtype B patients failing therapy. Position A/T400 in the connection subdomain is polymorphic, but the proportion of T400 increases from 41% in naïve patients to 72% in patients failing therapy. Previous studies suggested a role for threonine in conferring resistance to nucleoside RT inhibitors. Here we report that T400 also mediates resistance to non-nucleoside RT inhibitors. The susceptibility to NVP and EFV was reduced 5-fold and 2-fold, respectively, in the wild-type subtype B NL4.3 background. We show that substitution A400T reduces the RNaseH activity. The changes in enzyme activity are remarkable given the distance to both the polymerase and RNaseH active sites. Molecular dynamics simulations were performed, which provide a novel atomistic mechanism for the reduction in RNaseH activity induced by T400. Substitution A400T was found to change the conformation of the RNaseH primer grip region. Formation of an additional hydrogen bond between residue T400 and E396 may play a role in this structural change. The slower degradation of the viral RNA genome may provide more time for dissociation of the bound NNRTI from the stalled RT-template/primer complex, after which reverse transcription can resume.

Viral reverse transcriptases show selective high affinity binding to DNA-DNA primer-templates that resemble the polypurine tract

Previous results using a SELEX (Systematic Evolution of Ligands by Exponential Enrichment)-based approach that selected DNA primer-template duplexes binding with high affinity to HIV reverse transcriptase (RT) showed that primers mimicking the 3' end, and in particular the six nt terminal G tract, of the RNA polypurine tract (PPT; HIV PPT: 5'-AAAAGAAAAGGGGGG-3') were preferentially selected. In this report, two viral (Moloney murine leukemia virus (MuLV) and avian myeloblastosis virus (AMV)) and one retrotransposon (Ty3) RTs were used for selection. Like HIV RT, both viral RTs selected duplexes with primer strands mimicking the G tract at the PPT 3' end (AMV PPT: 5'-AGGGAGGGGGA-3'; MuLV PPT: 5'-AGAAAAAGGGGGG-3'). In contrast, Ty3, whose PPT lacks a G tract (5'-GAGAGAGAGGAA-3') showed no selective binding to any duplex sequences. Experiments were also conducted with DNA duplexes (termed DNA PPTs) mimicking the RNA PPT-DNA duplex of each virus and a control duplex with a random DNA sequence. Retroviral RTs bound with high affinity to all viral DNA PPT constructs, with HIV and MuLV RTs showing comparable binding to the counterpart DNA PPT duplexes and reduced affinity to the AMV DNA PPT. AMV RT showed similar behavior with a modest preference for its own DNA PPT. Ty3 RT showed no preferential binding for its own or any other DNA PPT and viral RTs bound the Ty3 DNA PPT with relatively low affinity. In contrast, binding affinity of HIV RT to duplexes containing the HIV RNA PPT was less dependent on the G tract, which is known to be pivotal for efficient extension. We hypothesize that the G tract on the RNA PPT helps shift the binding orientation of RT to the 3' end of the PPT where extension can occur.

Probing Retroviral and Retrotransposon Genome Structures: The "SHAPE" of Things to Come

Understanding the nuances of RNA structure as they pertain to biological function remains a formidable challenge for retrovirus research and development of RNA-based therapeutics, an area of particular importance with respect to combating HIV infection. Although a variety of chemical and enzymatic RNA probing techniques have been successfully employed for more than 30 years, they primarily interrogate small (100–500 nt) RNAs that have been removed from their biological context, potentially eliminating long-range tertiary interactions (such as kissing loops and pseudoknots) that may play a critical regulatory role. Selective 2′ hydroxyl acylation analyzed by primer extension (SHAPE), pioneered recently by Merino and colleagues, represents a facile, user-friendly technology capable of interrogating RNA structure with a single reagent and, combined with automated capillary electrophoresis, can analyze an entire 10,000-nucleotide RNA genome in a matter of weeks. Despite these obvious advantages, SHAPE essentially provides a nucleotide “connectivity map,” conversion of which into a 3-D structure requires a variety of complementary approaches. This paper summarizes contributions from SHAPE towards our understanding of the structure of retroviral genomes, modifications to which technology that have been developed to address some of its limitations, and future challenges.

The mechanisms of RNA SHAPE chemistry

The biological functions of RNA are ultimately governed by the local environment at each nucleotide. Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry is a powerful approach for measuring nucleotide structure and dynamics in diverse biological environments. SHAPE reagents acylate the 2'-hydroxyl group at flexible nucleotides because unconstrained nucleotides preferentially sample rare conformations that enhance the nucleophilicity of the 2'-hydroxyl. The critical corollary is that some constrained nucleotides must be poised for efficient reaction at the 2'-hydroxyl group. To identify such nucleotides, we performed SHAPE on intact crystals of the Escherichia coli ribosome, monitored the reactivity of 1490 nucleotides in 16S rRNA, and examined those nucleotides that were hyper-reactive toward SHAPE and had well-defined crystallographic conformations. Analysis of these conformations revealed that 2'-hydroxyl reactivity is broadly facilitated by general base catalysis involving multiple RNA functional groups and by two specific orientations of the bridging 3'-phosphate group. Nucleotide analog studies confirmed the contributions of these mechanisms to SHAPE reactivity. These results provide a strong mechanistic explanation for the relationship between SHAPE reactivity and local RNA dynamics and will facilitate interpretation of SHAPE information in the many technologies that make use of this chemistry.

Correlating SHAPE signatures with three-dimensional RNA structures

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a facile technique for quantitative analysis of RNA secondary structure. In general, low SHAPE signal values indicate Watson-Crick base-pairing, and high values indicate positions that are single-stranded within the RNA structure. However, the relationship of SHAPE signals to structural properties such as non-Watson-Crick base-pairing or stacking has thus far not been thoroughly investigated. Here, we present results of SHAPE experiments performed on several RNAs with published three-dimensional structures. This strategy allows us to analyze the results in terms of correlations between chemical reactivities and structural properties of the respective nucleotide, such as different types of base-pairing, stacking, and phosphate-backbone interactions. We find that the RNA SHAPE signal is strongly correlated with cis-Watson-Crick/Watson-Crick base-pairing and is to a remarkable degree not dependent on other structural properties with the exception of stacking. We subsequently generated probabilistic models that estimate the likelihood that a residue with a given SHAPE score participates in base-pairing. We show that several models that take SHAPE scores of adjacent residues into account perform better in predicting base-pairing compared with individual SHAPE scores. This underscores the context sensitivity of SHAPE and provides a framework for an improved interpretation of the response of RNA to chemical modification.

Investigation of the Reactivity of Oligodeoxynucleotides with Glyoxal and KMnO(4) Chemical Probes by Electrospray Ionization Mass Spectrometry

The reactions of two well-known chemical probes, glyoxal and potassium permanganate (KMnO(4)), with oligodeoxynucleotides were monitored by electrospray ionization (ESI) mass spectrometry to evaluate the influence of the sequence of DNA, its secondary structure, and interactions with associated ligands on the reactivity of the two probes. Glyoxal, a guanine-reactive probe, incorporated a mass shift of 58 Da, and potassium permanganate (KMnO(4)) is a thymine-reactive probe that resulted in a mass shift of 34 Da. The reactions depended on the accessibility of the nucleobases, and the peak abundances of the adducts in the ESI-mass spectra were used to quantify the extent of the chemical probe reactions. In this study, both mixed-base sequences were studied as well as control sequences in which one reactive site was located at the terminus or center of the oligodeoxynucleotide while the surrounding bases were a second, different nucleobase. In addition, the reactions of the chemical probes with non-covalent complexes formed between DNA and either actinomycin D or ethidium bromide, both known to interact with single strand DNA, were evaluated.

MS analysis of nucleic acids in the post-genomic era

Alternative approaches complementing the existing technologies for analysis of nucleic acids and their assemblies are necessary to take on the new challenges posed by the postgenomic era. The versatility of MS in biopolymer analysis and its ability to reach beyond sequence information are the basis of ever expanding applications aimed at the elucidation of nucleic acid structure-function relationships. This Feature summarizes the current state of MS-based approaches devised to overcome the limitations of traditional techniques and to advance different facets of nucleic acids research.

Biophysical analysis of influenza A virus RNA promoter at physiological temperatures

Each segment of the influenza A virus (IAV) genome contains conserved sequences at the 5′- and 3′-terminal ends, which form the promoter region necessary for polymerase binding and initiation of RNA synthesis. Although several models of interaction have been proposed it remains unclear if these two short, partially complementary, and highly conserved sequences can form a stable RNA duplex at physiological temperatures. First, our time-resolved FRET analysis revealed that a 14-mer 3′-RNA and a 15-mer 5′-RNA associate in solution, even at 42 °C. We also found that a nonfunctional RNA promoter containing the 3′-G3U mutation, as well as a promoter containing the compensatory 3′-G3U/C8A mutations, was able to form a duplex as efficiently as wild type. Second, UV melting analysis demonstrated that the wild-type and mutant RNA duplexes have similar stabilities in solution. We also observed an increase in thermostability for a looped promoter structure. The absence of differences in the stability and binding kinetics between wild type and a nonfunctional sequence suggests that the IAV promoter can be functionally inactivated without losing the capability to form a stable RNA duplex. Finally, using uridine specific chemical probing combined with mass spectrometry, we confirmed that the 5′ and 3′ sequences form a duplex which protects both RNAs from chemical modification, consistent with the previously published panhandle structure. These data support that these short, conserved promoter sequences form a stable complex at physiological temperatures, and this complex likely is important for polymerase recognition and viral replication.

Characterization of oligodeoxynucleotides and modifications by 193 nm photodissociation and electron photodetachment dissociation

Ultraviolet photodissociation (UVPD) at 193 nm is compared to collision induced dissociation (CID) for sequencing and determination of modifications of multideprotonated 6-20-mer oligodeoxynucleotides. UVPD at 193 nm causes efficient charge reduction of the deprotonated oligodeoxynucleotides via electron detachment, in addition to extensive backbone cleavages to yield sequence ions of relatively low abundance, including w, x, y, z, a, a-B, b, c, and d ions. Although internal ions populate UVPD spectra, base loss ions from the precursor are absent. Subsequent CID of the charge-reduced oligodeoxynucleotides formed upon electron detachment, in a net process called electron photodetachment dissociation (EPD), results in abundant sequence ions in terms of w, z, a, a-B, and d products, with a marked decrease in the abundance of precursor base loss ions and internal fragments. Complete sequencing was possible for virtually all oligodeoxynucleotides studied. EPD of three modified oligodeoxynucleotides, a methylated oligodeoxynucleotide, a phosphorothioate-modified oligodeoxynucleotide, and an ethylated-oligodeoxynucleotide, resulted in specific and extensive backbone cleavages, specifically, w, z, a, a-B, and d products, which allowed the modification site(s) to be pinpointed to a more specific location than by conventional CID.

Infrared multiphoton dissociation of small-interfering RNA anions and cations

Infrared multiphoton dissociation (IRMPD) on a linear ion trap mass spectrometer is applied for the sequencing of small interfering RNA (siRNA). Both single-strand siRNAs and duplex siRNA were characterized by IRMPD, and the results were compared with that obtained by traditional ion trap-based collision induced dissociation (CID). The single-strand siRNA anions were observed to dissociate via cleavage of the 5' P-O bonds yielding c- and y-type product ions as well as undergo neutral base loss. Full sequence coverage of the siRNA anions was obtained by both IRMPD and CID. While the CID mass spectra were dominated by base loss ions, accounting for approximately 25% to 40% of the product ion current, these ions were eliminated through secondary dissociation by increasing the irradiation time in the IRMPD mass spectra to produce higher abundances of informative sequence ions. With longer irradiation times, however, internal ions corresponding to cleavage of two 5' P-O bonds began to populate the product ion mass spectra as well as higher abundances of [a - Base] and w-type ions. IRMPD of siRNA cations predominantly produced c- and y-type ions with minimal contributions of [a - Base] and w-type ions to the product ion current; the presence of only two complementary series of product ions in the IRMPD mass spectra simplified spectral interpretation. In addition, IRMPD produced high abundances of protonated nucleobases, [G + H](+), [A + H](+), and [C + H](+), which were not detected in the CID mass spectra due to the low-mass cut-off associated with conventional CID in ion traps. CID and IRMPD using short irradiation times of duplex siRNA resulted in strand separation, similar to the dissociation trends observed for duplex DNA. With longer irradiation times, however, the individual single-strands underwent secondary dissociation to yield informative sequence ions not obtained by CID.

Revisiting plus-strand DNA synthesis in retroviruses and long terminal repeat retrotransposons: dynamics of enzyme: substrate interactions

Although polypurine tract (PPT)-primed initiation of plus-strand DNA synthesis in retroviruses and LTR-containing retrotransposons can be accurately duplicated, the molecular details underlying this concerted series of events remain largely unknown. Importantly, the PPT 3' terminus must be accommodated by ribonuclease H (RNase H) and DNA polymerase catalytic centers situated at either terminus of the cognate reverse transcriptase (RT), and in the case of the HIV-1 enzyme, ∼70Å apart. Communication between RT and the RNA/DNA hybrid therefore appears necessary to promote these events. The crystal structure of the HIV-1 RT/PPT complex, while informative, positions the RNase H active site several bases pairs from the PPT/U3 junction, and thus provides limited information on cleavage specificity. To fill the gap between biochemical and crystallographic approaches, we review a multidisciplinary approach combining chemical probing, mass spectrometry, NMR spectroscopy and single molecule spectroscopy. Our studies also indicate that nonnucleoside RT inhibitors affect enzyme orientation, suggesting initiation of plus-strand DNA synthesis as a potential therapeutic target.