Mapping nucleic acid structure by hydroxyl radical cleavage

Systematic Fe(II)-EDTA Method of Dose-Dependent Hydroxyl Radical Generation for Protein Oxidative Footprinting

Correlating the structure and dynamics of proteins with biological function is critical to understanding normal and dysfunctional cellular mechanisms. We describe a quantitative method of hydroxyl radical generation via Fe(II)-ethylenediaminetetraacetic acid (EDTA)-catalyzed Fenton chemistry that provides ready access to protein oxidative footprinting using equipment commonly found in research and process control laboratories. Robust and reproducible dose-dependent oxidation of protein samples is observed and quantitated by mass spectrometry with as fine a single residue resolution. An oxidation analysis of lysozyme provides a readily accessible benchmark for our method. The efficacy of our oxidation method is demonstrated by mapping the interface of a RAS-monobody complex, the surface of the NIST mAb, and the interface between PRC2 complex components. These studies are executed using standard laboratory tools and a few pennies of reagents; the mass spectrometry analysis can be streamlined to map the protein structure with single amino acid residue resolution.

Systematic detection of tertiary structural modules in large RNAs and RNP interfaces by Tb-seq

Compact RNA structural motifs control many aspects of gene expression, but we lack methods for finding these structures in the vast expanse of multi-kilobase RNAs. To adopt specific 3-D shapes, many RNA modules must compress their RNA backbones together, bringing negatively charged phosphates into close proximity. This is often accomplished by recruiting multivalent cations (usually Mg2+), which stabilize these sites and neutralize regions of local negative charge. Coordinated lanthanide ions, such as terbium (III) (Tb3+), can also be recruited to these sites, where they induce efficient RNA cleavage, thereby revealing compact RNA 3-D modules. Until now, Tb3+ cleavage sites were monitored via low-throughput biochemical methods only applicable to small RNAs. Here we present Tb-seq, a high-throughput sequencing method for detecting compact tertiary structures in large RNAs. Tb-seq detects sharp backbone turns found in RNA tertiary structures and RNP interfaces, providing a way to scan transcriptomes for stable structural modules and potential riboregulatory motifs.

"Bind, cleave and leave": multiple turnover catalysis of RNA cleavage by bulge-loop inducing supramolecular conjugates

Antisense sequence-specific knockdown of pathogenic RNA offers opportunities to find new solutions for therapeutic treatments. However, to gain a desired therapeutic effect, the multiple turnover catalysis is critical to inactivate many copies of emerging RNA sequences, which is difficult to achieve without sacrificing the sequence-specificity of cleavage. Here, engineering two or three catalytic peptides into the bulge-loop inducing molecular framework of antisense oligonucleotides achieved catalytic turnover of targeted RNA. Different supramolecular configurations revealed that cleavage of the RNA backbone upon sequence-specific hybridization with the catalyst accelerated with increase in the number of catalytic guanidinium groups, with almost complete demolition of target RNA in 24 h. Multiple sequence-specific cuts at different locations within and around the bulge-loop facilitated release of the catalyst for subsequent attacks of at least 10 further RNA substrate copies, such that delivery of only a few catalytic molecules could be sufficient to maintain knockdown of typical RNA copy numbers. We have developed fluorescent assay and kinetic simulation tools to characterise how the limited availability of different targets and catalysts had restrained catalytic reaction progress considerably, and to inform how to accelerate the catalytic destruction of shorter linear and larger RNAs even further.

Hydroxyl Radical Protein Footprinting: A Mass Spectrometry-Based Structural Method for Studying the Higher Order Structure of Proteins

Hydroxyl radical protein footprinting (HRPF) coupled to mass spectrometry has been successfully used to investigate a plethora of protein-related questions. The method, which utilizes hydroxyl radicals to oxidatively modify solvent-accessible amino acids, can inform on protein interaction sites and regions of conformational change. Hydroxyl radical-based footprinting was originally developed to study nucleic acids, but coupling the method with mass spectrometry has enabled the study of proteins. The method has undergone several advancements since its inception that have increased its utility for more varied applications such as protein folding and the study of biotherapeutics. In addition, recent innovations have led to the study of increasingly complex systems including cell lysates and intact cells. Technological advances have also increased throughput and allowed for better control of experimental conditions. In this review, we provide a brief history of the field of HRPF and detail recent innovations and applications in the field.

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.

RNA Ensembles from Solvent Accessibility Data: Application to the SAM-I Riboswitch Aptamer Domain

Riboswitches are regulatory ribonucleic acid (RNA) elements that act as ligand-dependent conformational switches that recognize their cognate ligand via a binding pocket located in their aptamer domain. In the apo form, the aptamer domain is dynamic, requiring an ensemble representation of its structure. Here, as a proof-of-concept, we used solvent accessibility information to construct a pair of dynamical ensembles of the aptamer domain of the well-studied S-adenosylmethionine (SAM) class-I riboswitch in the absence (-SAM) and presence (+SAM) of SAM. To achieve this, we first generated a large conformational library and then reweighted conformers in the library using solvent-accessible surface area (SASA) data derived from recently reported light-activated structural examination of RNA (LASER) reactivities measured in the -SAM and +SAM states of the riboswitch. The differences in the resulting -SAM and +SAM ensembles are consistent with a SAM-dependent reshaping of the free-energy landscape of the aptamer domain. Within our -SAM ensemble, we identified a "transient" state that is missing a critical long-range contact, leading us to speculate that it may be representative of a folding intermediate. Further structural analysis also revealed that the transient state harbors a hidden binding pocket that is distinct from the SAM-binding pocket and is predicted by docking calculations to selectively bind small-molecule ligands. The SASA-based method we applied to the SAM-I riboswitch aptamer domain is general and could be used to construct dynamical ensembles for other riboswitch aptamer domains and, more broadly, other classes of structured RNAs.

Mass Spectrometry-Based Protein Footprinting for Higher-Order Structure Analysis: Fundamentals and Applications

Proteins adopt different higher-order structures (HOS) to enable their unique biological functions. Understanding the complexities of protein higher-order structures and dynamics requires integrated approaches, where mass spectrometry (MS) is now positioned to play a key role. One of those approaches is protein footprinting. Although the initial demonstration of footprinting was for the HOS determination of protein/nucleic acid binding, the concept was later adapted to MS-based protein HOS analysis, through which different covalent labeling approaches "mark" the solvent accessible surface area (SASA) of proteins to reflect protein HOS. Hydrogen-deuterium exchange (HDX), where deuterium in D2O replaces hydrogen of the backbone amides, is the most common example of footprinting. Its advantage is that the footprint reflects SASA and hydrogen bonding, whereas one drawback is the labeling is reversible. Another example of footprinting is slow irreversible labeling of functional groups on amino acid side chains by targeted reagents with high specificity, probing structural changes at selected sites. A third footprinting approach is by reactions with fast, irreversible labeling species that are highly reactive and footprint broadly several amino acid residue side chains on the time scale of submilliseconds. All of these covalent labeling approaches combine to constitute a problem-solving toolbox that enables mass spectrometry as a valuable tool for HOS elucidation. As there has been a growing need for MS-based protein footprinting in both academia and industry owing to its high throughput capability, prompt availability, and high spatial resolution, we present a summary of the history, descriptions, principles, mechanisms, and applications of these covalent labeling approaches. Moreover, their applications are highlighted according to the biological questions they can answer. This review is intended as a tutorial for MS-based protein HOS elucidation and as a reference for investigators seeking a MS-based tool to address structural questions in protein science.

Binding interface and impact on protease cleavage for an RNA aptamer to HIV-1 reverse transcriptase

RNA aptamers that bind HIV-1 reverse transcriptase (RT) inhibit RT in enzymatic and viral replication assays. Some aptamers inhibit RT from only a few viral clades, while others show broad-spectrum inhibition. Biophysical determinants of recognition specificity are poorly understood. We investigated the interface between HIV-1 RT and a broad–spectrum UCAA-family aptamer. SAR and hydroxyl radical probing identified aptamer structural elements critical for inhibition and established the role of signature UCAA bulge motif in RT-aptamer interaction. HDX footprinting on RT ± aptamer shows strong contacts with both subunits, especially near the C-terminus of p51. Alanine scanning revealed decreased inhibition by the aptamer for mutants P420A, L422A and K424A. 2D proton nuclear magnetic resonance and SAXS data provided constraints on the solution structure of the aptamer and enable computational modeling of the docked complex with RT. Surprisingly, the aptamer enhanced proteolytic cleavage of precursor p66/p66 by HIV-1 protease, suggesting that it stabilizes the productive conformation to allow maturation. These results illuminate features at the RT-aptamer interface that govern recognition specificity by a broad-spectrum antiviral aptamer, and they open new possibilities for accelerating RT maturation and interfering with viral replication.

The evolution of RNA structural probing methods: From gels to next-generation sequencing

RNA molecules are important players in all domains of life and the study of the relationship between their multiple flexible states and the associated biological roles has increased in recent years. For several decades, chemical and enzymatic structural probing experiments have been used to determine RNA structure. During this time, there has been a steady improvement in probing reagents and experimental methods, and today the structural biologist community has a large range of tools at its disposal to probe the secondary structure of RNAs in vitro and in cells. Early experiments used radioactive labeling and polyacrylamide gel electrophoresis as read-out methods. This was superseded by capillary electrophoresis, and more recently by next-generation sequencing. Today, powerful structural probing methods can characterize RNA structure on a genome-wide scale. In this review, we will provide an overview of RNA structural probing methodologies from a historical and technical perspective. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Methods > RNA Analyses in vitro and In Silico RNA Methods > RNA Analyses in Cells.

Extracting information from RNA SHAPE data: Kalman filtering approach

RNA SHAPE experiments have become important and successful sources of information for RNA structure prediction. In such experiments, chemical reagents are used to probe RNA backbone flexibility at the nucleotide level, which in turn provides information on base pairing and therefore secondary structure. Little is known, however, about the statistics of such SHAPE data. In this work, we explore different representations of noise in SHAPE data and propose a statistically sound framework for extracting reliable reactivity information from multiple SHAPE replicates. Our analyses of RNA SHAPE experiments underscore that a normal noise model is not adequate to represent their data. We propose instead a log-normal representation of noise and discuss its relevance. Under this assumption, we observe that processing simulated SHAPE data by directly averaging different replicates leads to bias. Such bias can be reduced by analyzing the data following a log transformation, either by log-averaging or Kalman filtering. Application of Kalman filtering has the additional advantage that a prior on the nucleotide reactivities can be introduced. We show that the performance of Kalman filtering is then directly dependent on the quality of that prior. We conclude the paper with guidelines on signal processing of RNA SHAPE data.

Core-Shell Nanostructured Fe3O4-Poly(styrene-co-vinylbenzyl chloride) Grafted PPI Dendrimers Stabilized with AuNPs/PdNPs for Efficient Nuclease Activity

Four different novel magnetic core-shell nanocomposites stabilized with Au/Pd nanoparticles (NPs) were prepared by a simple procedure and demonstrated their catalytic activity for effective cleavage of pBR322 DNA. Initially, the Fe₃O₄-poly(styrene-divinylbenzene-vinylbenzyl chloride) (ST-DVB-VBC) matrix functionalized with 3-aminobenzoic acid was prepared and grafted with PPI-G(2) and PPI-G(3) dendrimers. Each core-shell matrix was immobilized with AuNPs and PdNPs separately. The resulting composites were characterized by FT-IR, UV-vis, SEM, TEM, XRD, VSM, XPS, Raman, and TGA. The magnetic core-shell nanocomposites at concentrations from 30 to 50 μM were employed separately to study DNA cleavage by agarose gel electrophoresis. Among the four magnetic core-shell nanocomposites, Fe₃O₄-poly(ST-DVB-VBC)-PPI-G(3)-AuNPs showed higher activity than others for DNA cleavage, and formed Form-II and -III DNA. When the concentration of Fe₃O₄-poly(ST-DVB-VBC)-PPI-G(3)-AuNPs was increased from 40 to 45 and 45 to 50 μM, Form-III (linear) DNA was observed with 10 and 22%, respectively, in addition to Form-II. This observation suggests formation of linear DNA from the supercoiled DNA via nicked DNA-intermediated consecutive cleaving process. The magnetic core-shell nanocomposites were stable and monodispersed, and exhibited rapid magnetic response. These properties are crucial for their application in biomolecular separations and targeted drug-delivery in the future.

Covalent labeling-mass spectrometry with non-specific reagents for studying protein structure and interactions

Using mass spectrometry (MS) to obtain information about a higher order structure of protein requires that a protein's structural properties are encoded into the mass of that protein. Covalent labeling (CL) with reagents that can irreversibly modify solvent accessible amino acid side chains is an effective way to encode structural information into the mass of a protein, as this information can be read-out in a straightforward manner using standard MS-based proteomics techniques. The differential reactivity of proteins under two or more conditions can be used to distinguish protein topologies, conformations, and/or binding sites. CL-MS methods have been effectively used for the structural analysis of proteins and protein complexes, particularly for systems that are difficult to study by other more traditional biochemical techniques. This review provides an overview of the non-specific CL approaches that have been combined with MS with a particular emphasis on the reagents that are commonly used, including hydroxyl radicals, carbenes, and diethylpyrocarbonate. We describe the reagent and protein factors that affect the reactivity of amino acid side chains. We also include details about experimental design and workflow, data analysis, recent applications, and some future prospects of CL-MS methods.

PATTERNA: transcriptome-wide search for functional RNA elements via structural data signatures

Establishing a link between RNA structure and function remains a great challenge in RNA biology. The emergence of high-throughput structure profiling experiments is revolutionizing our ability to decipher structure, yet principled approaches for extracting information on structural elements directly from these data sets are lacking. We present PATTERNA, an unsupervised pattern recognition algorithm that rapidly mines RNA structure motifs from profiling data. We demonstrate that PATTERNA detects motifs with an accuracy comparable to commonly used thermodynamic models and highlight its utility in automating data-directed structure modeling from large data sets. PATTERNA is versatile and compatible with diverse profiling techniques and experimental conditions.

Light-activated chemical probing of nucleobase solvent accessibility inside cells

The discovery of functional RNAs that are critical for normal and disease physiology continues to expand at a breakneck pace. Many RNA functions are controlled by the formation of specific structures, and an understanding of each structural component is necessary to elucidate its function. Measuring solvent accessibility intracellularly with experimental ease is an unmet need in the field. Here, we present a novel method for probing nucleobase solvent accessibility, Light Activated Structural Examination of RNA (LASER). LASER depends on light activation of a small molecule, nicotinoyl azide (NAz), to measure solvent accessibility of purine nucleobases. In vitro, this technique accurately monitors solvent accessibility and identifies rapid structural changes resulting from ligand binding in a metabolite-responsive RNA. LASER probing can further identify cellular RNA-protein interactions and unique intracellular RNA structures. Our photoactivation technique provides an adaptable framework to structurally characterize solvent accessibility of RNA in many environments.

Analyzing the structure of macromolecules in their native cellular environment using hydroxyl radical footprinting

Hydroxyl radical footprinting (HRF) has been successfully used to study the structure of both nucleic acids and proteins in live cells.

A structural map of oncomiR-1 at single-nucleotide resolution

The miR-17-92a cluster, also known as 'oncomiR-1', is an RNA transcript that plays a pivotal regulatory role in cellular processes, including the cell cycle, proliferation and apoptosis. Its dysregulation underlies the development of several cancers. Oncomir-1 comprises six constituent miRNAs, each processed with different efficiencies as a function of both developmental time and tissue type. The structural mechanisms that regulate such differential processing are unknown, and this has impeded our understanding of the dysregulation of oncomiR-1 in pathophysiology. By probing the sensitivity of each nucleotide in oncomiR-1 to reactive small molecules, we present a secondary structural map of this RNA at single-nucleotide resolution. The secondary structure and solvent accessible regions of oncomiR-1 reveal that most of its primary microRNA domains are suboptimal substrates for Drosha-DGCR8, and therefore resistant to microprocessing. The structure indicates that the binding of trans-acting factors is required to remodel the tertiary organization and unmask cryptic primary microRNA domains to facilitate their processing into pre-microRNAs.

Dawn of the in vivo RNA structurome and interactome

RNA is one of the most fascinating biomolecules in living systems given its structural versatility to fold into elaborate architectures for important biological functions such as gene regulation, catalysis, and information storage. Knowledge of RNA structures and interactions can provide deep insights into their functional roles in vivo For decades, RNA structural studies have been conducted on a transcript-by-transcript basis. The advent of next-generation sequencing (NGS) has enabled the development of transcriptome-wide structural probing methods to profile the global landscape of RNA structures and interactions, also known as the RNA structurome and interactome, which transformed our understanding of the RNA structure-function relationship on a transcriptomic scale. In this review, molecular tools and NGS methods used for RNA structure probing are presented, novel insights uncovered by RNA structurome and interactome studies are highlighted, and perspectives on current challenges and potential future directions are discussed. A more complete understanding of the RNA structures and interactions in vivo will help illuminate the novel roles of RNA in gene regulation, development, and diseases.

Properties and reactivity of nucleic acids relevant to epigenomics, transcriptomics, and therapeutics

Developments in epigenomics, toxicology, and therapeutic nucleic acids all rely on a precise understanding of nucleic acid properties and chemical reactivity. In this review we discuss the properties and chemical reactivity of each nucleobase and attempt to provide some general principles for nucleic acid targeting or engineering. For adenine-thymine and guanine-cytosine base pairs, we review recent quantum chemical estimates of their Watson-Crick interaction energy, π-π stacking energies, as well as the nuclear quantum effects on tautomerism. Reactions that target nucleobases have been crucial in the development of new sequencing technologies and we believe further developments in nucleic acid chemistry will be required to deconstruct the enormously complex transcriptome.

Genetic variation and RNA structure regulate microRNA biogenesis

MiRNA biogenesis is highly regulated at the post-transcriptional level; however, the role of sequence and secondary RNA structure in this process has not been extensively studied. A single G to A substitution present in the terminal loop of pri-mir-30c-1 in breast and gastric cancer patients had been previously described to result in increased levels of mature miRNA. Here, we report that this genetic variant directly affects Drosha-mediated processing of pri-mir-30c-1 in vitro and in cultured cells. Structural analysis of this variant revealed an altered RNA structure that facilitates the interaction with SRSF3, an SR protein family member that promotes pri-miRNA processing. Our results are compatible with a model whereby a genetic variant in pri-mir-30c-1 leads to a secondary RNA structure rearrangement that facilitates binding of SRSF3 resulting in increased levels of miR-30c. These data highlight that primary sequence determinants and RNA structure are key regulators of miRNA biogenesis.

Osmium tetroxide as a probe of RNA structure

Structured RNAs have a central role in cellular function. The capability of structured RNAs to adopt fixed architectural structures or undergo dynamic conformational changes contributes to their diverse role in the regulation of gene expression. Although numerous biophysical and biochemical tools have been developed to study structured RNAs, there is a continuing need for the development of new methods for the investigation of RNA structures, especially methods that allow RNA structure to be studied in solution close to its native cellular conditions. Here we use osmium tetroxide (OsO₄) as a chemical probe of RNA structure. In this method, we have used fluorescence-based sequencing technologies to detect OsO₄ modified RNA. We characterized the requirements for OsO₄ modification of RNA by investigating three known structured RNAs: the M-box, glycine riboswitch RNAs, and tRNA^asp Our results show that OsO₄ predominantly modifies RNA at uracils that are conformationally exposed on the surface of the RNA. We also show that changes in OsO₄ reactivity at flexible positions in the RNA correlate with ligand-driven conformational changes in the RNA structure. Osmium tetroxide modification of RNA will provide insights into the structural features of RNAs that are relevant to their underlying biological functions.

Comparative and integrative analysis of RNA structural profiling data: current practices and emerging questions

BACKGROUND

Structure profiling experiments provide single-nucleotide information on RNA structure. Recent advances in chemistry combined with application of high-throughput sequencing have enabled structure profiling at transcriptome scale and in living cells, creating unprecedented opportunities for RNA biology. Propelled by these experimental advances, massive data with ever-increasing diversity and complexity have been generated, which give rise to new challenges in interpreting and analyzing these data.

RESULTS

We review current practices in analysis of structure profiling data with emphasis on comparative and integrative analysis as well as highlight emerging questions. Comparative analysis has revealed structural patterns across transcriptomes and has become an integral component of recent profiling studies. Additionally, profiling data can be integrated into traditional structure prediction algorithms to improve prediction accuracy.

CONCLUSIONS

To keep pace with experimental developments, methods to facilitate, enhance and refine such analyses are needed. Parallel advances in analysis methodology will complement profiling technologies and help them reach their full potential.

Synthesis and characterization of FeIII(μ-OH)ZnII complexes: effects of a second coordination sphere and increase in the chelate ring size on the hydrolysis of a phosphate diester and DNA

Effects of a second coordination sphere and of the chelate ring size in Fe^III(μ-OH)Zn^II complexes properties and catalysis.

Data-directed RNA secondary structure prediction using probabilistic modeling

Structure dictates the function of many RNAs, but secondary RNA structure analysis is either labor intensive and costly or relies on computational predictions that are often inaccurate. These limitations are alleviated by integration of structure probing data into prediction algorithms. However, existing algorithms are optimized for a specific type of probing data. Recently, new chemistries combined with advances in sequencing have facilitated structure probing at unprecedented scale and sensitivity. These novel technologies and anticipated wealth of data highlight a need for algorithms that readily accommodate more complex and diverse input sources. We implemented and investigated a recently outlined probabilistic framework for RNA secondary structure prediction and extended it to accommodate further refinement of structural information. This framework utilizes direct likelihood-based calculations of pseudo-energy terms per considered structural context and can readily accommodate diverse data types and complex data dependencies. We use real data in conjunction with simulations to evaluate performances of several implementations and to show that proper integration of structural contexts can lead to improvements. Our tests also reveal discrepancies between real data and simulations, which we show can be alleviated by refined modeling. We then propose statistical preprocessing approaches to standardize data interpretation and integration into such a generic framework. We further systematically quantify the information content of data subsets, demonstrating that high reactivities are major drivers of SHAPE-directed predictions and that better understanding of less informative reactivities is key to further improvements. Finally, we provide evidence for the adaptive capability of our framework using mock probe simulations.

Probing the structure of ribosome assembly intermediates in vivo using DMS and hydroxyl radical footprinting

The assembly of the Escherichia coli ribosome has been widely studied and characterized in vitro. Despite this, ribosome biogenesis in living cells is only partly understood because assembly is coupled with transcription, modification and processing of the pre-ribosomal RNA. We present a method for footprinting and isolating pre-rRNA as it is synthesized in E. coli cells. Pre-rRNA synthesis is synchronized by starvation, followed by nutrient upshift. RNA synthesized during outgrowth is metabolically labeled to facilitate isolation of recent transcripts. Combining this technique with two in vivo RNA probing methods, hydroxyl radical and DMS footprinting, allows the structure of nascent RNA to be probed over time. Together, these can be used to determine changes in the structures of ribosome assembly intermediates as they fold in vivo.

Interactions of two large antiviral polyamides with the long control region of HPV16

PA1 and PA25 are large hairpin polyamides that are effective in nearly eliminating HPV16 episomes (DNA) in cell culture, and PA25 has broad spectrum activity against three cancer-causing forms of HPV (Edwards, T. G., Koeller, K. J., Slomczynska, U., Fok, K., Helmus, M., Bashkin, J. K., Fisher, C., Antiviral Res. 91 (2011) 177-186). Described here are the interactions of these PAs with sequences in the long control region (LCR) of HPV16 (7348-122). Using an FeEDTA conjugate of PA1 (designed to recognize 5'-W2GW7-3'; W = A or T), 34 affinity cleavage (AC) patterns were detected for this fragment. These sites can be rationalized with sequences featuring perfect, single, double, triple and quadruple mismatches. Quantitative DNase I footprinting analysis indicates that perfect sites bind PA1 with Kds between 0.7 and 2.2 nM. Kds for single, double, triple and quadruple mismatch sites range from 1-3 nM-20 nM. Using AC and EDTA conjugates, we report that unlike smaller 8-ring hairpin PAs, introduction of a chiral turn in this large polyamide has no effect on binding orientation (forward vs. reverse). Despite its design to recognize 5'-W2GW5GW4-3' via two Im residues, a motif not represented in this HPV sequence, a PA25-EDTA conjugate yielded 31 affinity cleavage sites on the region. Low nM Kds for PA25 without EDTA indicates a high tolerance for triple and quadruple mismatches. While there is extensive coverage of the sequence examined, AC cleavage patterns for the two PAs show discrete binding events and do not overlap significantly. This indicates that within the context of A/T rich sequences, these PAs do not recognize a simple shared sequence-related feature of the DNA. These insights continue to inform the complex nature of large hairpin PA-DNA interactions and antiviral behavior.

Predicting RNA secondary structures from sequence and probing data

RNA secondary structures have proven essential for understanding the regulatory functions performed by RNA such as microRNAs, bacterial small RNAs, or riboswitches. This success is in part due to the availability of efficient computational methods for predicting RNA secondary structures. Recent advances focus on dealing with the inherent uncertainty of prediction by considering the ensemble of possible structures rather than the single most stable one. Moreover, the advent of high-throughput structural probing has spurred the development of computational methods that incorporate such experimental data as auxiliary information.

Structure Prediction: New Insights into Decrypting Long Noncoding RNAs

Long noncoding RNAs (lncRNAs), which form a diverse class of RNAs, remain the least understood type of noncoding RNAs in terms of their nature and identification. Emerging evidence has revealed that a small number of newly discovered lncRNAs perform important and complex biological functions such as dosage compensation, chromatin regulation, genomic imprinting, and nuclear organization. However, understanding the wide range of functions of lncRNAs related to various processes of cellular networks remains a great experimental challenge. Structural versatility is critical for RNAs to perform various functions and provides new insights into probing the functions of lncRNAs. In recent years, the computational method of RNA structure prediction has been developed to analyze the structure of lncRNAs. This novel methodology has provided basic but indispensable information for the rapid, large-scale and in-depth research of lncRNAs. This review focuses on mainstream RNA structure prediction methods at the secondary and tertiary levels to offer an additional approach to investigating the functions of lncRNAs.

Rsite2: an efficient computational method to predict the functional sites of noncoding RNAs

Noncoding RNAs (ncRNAs) represent a big class of important RNA molecules. Given the large number of ncRNAs, identifying their functional sites is becoming one of the most important topics in the post-genomic era, but available computational methods are limited. For the above purpose, we previously presented a tertiary structure based method, Rsite, which first calculates the distance metrics defined in Methods with the tertiary structure of an ncRNA and then identifies the nucleotides located within the extreme points in the distance curve as the functional sites of the given ncRNA. However, the application of Rsite is largely limited because of limited RNA tertiary structures. Here we present a secondary structure based computational method, Rsite2, based on the observation that the secondary structure based nucleotide distance is strongly positively correlated with that derived from tertiary structure. This makes it reasonable to replace tertiary structure with secondary structure, which is much easier to obtain and process. Moreover, we applied Rsite2 to three ncRNAs (tRNA (Lys), Diels-Alder ribozyme, and RNase P) and a list of human mitochondria transcripts. The results show that Rsite2 works well with nearly equivalent accuracy as Rsite but is much more feasible and efficient. Finally, a web-server, the source codes, and the dataset of Rsite2 are available at http://www.cuialb.cn/rsite2.

Reactivity of Nucleic Acid Radicals

Nucleic acid oxidation plays a vital role in the etiology and treatment of diseases, as well as aging. Reagents that oxidize nucleic acids are also useful probes of the biopolymers' structure and folding. Radiation scientists have contributed greatly to our understanding of nucleic acid oxidation using a variety of techniques. During the past two decades organic chemists have applied the tools of synthetic and mechanistic chemistry to independently generate and study the reactive intermediates produced by ionizing radiation and other nucleic acid damaging agents. This approach has facilitated resolving mechanistic controversies and lead to the discovery of new reactive processes.

Progress and Current Challenges in Modeling Large RNAs

Recent breakthroughs in next-generation sequencing technologies have led to the discovery of several classes of non-coding RNAs (ncRNAs). It is now apparent that RNA molecules are not only just carriers of genetic information but also key players in many cellular processes. While there has been a rapid increase in the number of ncRNA sequences deposited in various databases over the past decade, the biological functions of these ncRNAs are largely not well understood. Similar to proteins, RNA molecules carry out a function by forming specific three-dimensional structures. Understanding the function of a particular RNA therefore requires a detailed knowledge of its structure. However, determining experimental structures of RNA is extremely challenging. In fact, RNA-only structures represent just 1% of the total structures deposited in the PDB. Thus, computational methods that predict three-dimensional RNA structures are in high demand. Computational models can provide valuable insights into structure-function relationships in ncRNAs and can aid in the development of functional hypotheses and experimental designs. In recent years, a set of diverse RNA structure prediction tools have become available, which differ in computational time, input data and accuracy. This review discusses the recent progress and challenges in RNA structure prediction methods.

Correlation of Thermal Stability and Structural Distortion of DNA Interstrand Cross-Links Produced from Oxidized Abasic Sites with Their Selective Formation and Repair

C4'-oxidized (C4-AP) and C5'-oxidized abasic sites (DOB) that are produced following abstraction of a hydrogen atom from the DNA backbone reversibly form cross-links selectively with dA opposite a 3'-adjacent nucleotide, despite the comparable proximity of an opposing dA. A previous report on UvrABC incision of DNA substrates containing stabilized analogues of the ICLs derived from C4-AP and DOB also indicated that the latter is repaired more readily by nucleotide excision repair [Ghosh, S., and Greenberg, M. M. (2014) Biochemistry 53, 5958-5965]. The source for selective cross-link formation was probed by comparing the reactivity of ICL analogues of C4-AP and DOB that mimic the preferred and disfavored cross-links with that of reagents that indirectly detect distortion by reacting with the nucleobases. The disfavored C4-AP and DOB analogues were each more reactive than the corresponding preferred cross-link substrates, suggesting that the latter are more stable, which is consistent with selective ICL formation. In addition, the preferred DOB analogue is more reactive than the respective C4-AP ICL, which is consistent with its more efficient incision by UvrABC. The conclusions drawn from the chemical probing experiments are corroborated by UV melting studies. The preferred ICLs exhibit melting temperatures higher than those of the corresponding disfavored isomers. These studies suggest that oxidized abasic sites form reversible interstrand cross-links with dA opposite the 3'-adjacent thymidine because these products are more stable and the thermodynamic preference is reflected in the transition states for their formation.

The ins and outs of lncRNA structure: How, why and what comes next?

The field of structural biology has the unique advantage of being able to provide a comprehensive picture of biological mechanisms at the molecular and atomic level. Long noncoding RNAs (lncRNAs) represent the new frontier in the molecular biology of complex organisms yet remain the least characterised of all the classes of RNA. Thousands of new lncRNAs are being reported each year yet very little structural data exists for this rapidly expanding field. The length of lncRNAs ranges from 200 nt to over 100 kb in length and they generally exhibit low cellular abundance. Therefore, obtaining sufficient quantities of lncRNA to use for structural analysis is challenging. However, as technologies develop structures of lncRNAs are starting to emerge providing important information regarding their mechanism of action. Here we review the current methods used to determine the structure of lncRNA and lncRNA:protein complexes and describe the significant contribution structural biology has and will make to the field of lncRNA research. This article is part of a Special Issue entitled: Clues to long noncoding RNA taxonomy1, edited by Dr. Tetsuro Hirose and Dr. Shinichi Nakagawa.

Recent trends in the development of vitamin B12 derivatives for medicinal applications

This Feature Article highlights recent developments in the field of vitamin B12 derivatives for medicinal applications. The following topics are emphasized: (1) the development of aquacorrinoids for cyanide detection and detoxification, (2) the use of vitamin B12 conjugates and (3) antivitamins B12 for therapy and diagnosis, and (4) the design of corrinoids as activators of soluble guanylyl cyclase (sGC).

Predicting RNA 3D structure using a coarse-grain helix-centered model

A 3D model of RNA structure can provide information about its function and regulation that is not possible with just the sequence or secondary structure. Current models suffer from low accuracy and long running times and either neglect or presume knowledge of the long-range interactions which stabilize the tertiary structure. Our coarse-grained, helix-based, tertiary structure model operates with only a few degrees of freedom compared with all-atom models while preserving the ability to sample tertiary structures given a secondary structure. It strikes a balance between the precision of an all-atom tertiary structure model and the simplicity and effectiveness of a secondary structure representation. It provides a simplified tool for exploring global arrangements of helices and loops within RNA structures. We provide an example of a novel energy function relying only on the positions of stems and loops. We show that coupling our model to this energy function produces predictions as good as or better than the current state of the art tools. We propose that given the wide range of conformational space that needs to be explored, a coarse-grain approach can explore more conformations in less iterations than an all-atom model coupled to a fine-grain energy function. Finally, we emphasize the overarching theme of providing an ensemble of predicted structures, something which our tool excels at, rather than providing a handful of the lowest energy structures.

Reproducible Analysis of Sequencing-Based RNA Structure Probing Data with User-Friendly Tools

RNA structure-probing data can improve the prediction of RNA secondary and tertiary structure and allow structural changes to be identified and investigated. In recent years, massive parallel sequencing has dramatically improved the throughput of RNA structure probing experiments, but at the same time also made analysis of the data challenging for scientists without formal training in computational biology. Here, we discuss different strategies for data analysis of massive parallel sequencing-based structure-probing data. To facilitate reproducible and standardized analysis of this type of data, we have made a collection of tools, which allow raw sequencing reads to be converted to normalized probing values using different published strategies. In addition, we also provide tools for visualization of the probing data in the UCSC Genome Browser and for converting RNA coordinates to genomic coordinates and vice versa. The collection is implemented as functions in the R statistical environment and as tools in the Galaxy platform, making them easily accessible for the scientific community. We demonstrate the usefulness of the collection by applying it to the analysis of sequencing-based hydroxyl radical probing data and comparing different normalization strategies.

Computational analysis of RNA structures with chemical probing data

RNAs play various roles, not only as the genetic codes to synthesize proteins, but also as the direct participants of biological functions determined by their underlying high-order structures. Although many computational methods have been proposed for analyzing RNA structures, their accuracy and efficiency are limited, especially when applied to the large RNAs and the genome-wide data sets. Recently, advances in parallel sequencing and high-throughput chemical probing technologies have prompted the development of numerous new algorithms, which can incorporate the auxiliary structural information obtained from those experiments. Their potential has been revealed by the secondary structure prediction of ribosomal RNAs and the genome-wide ncRNA function annotation. In this review, the existing probing-directed computational methods for RNA secondary and tertiary structure analysis are discussed.

The pattern of DNA cleavage intensity around indels

Indels (insertions and deletions) are the second most common form of genetic variations in the eukaryotic genomes and are responsible for a multitude of genetic diseases. Despite its significance, detailed molecular mechanisms for indel generation are still unclear. Here we examined 2,656,597 small human and mouse germline indels, 16,742 human somatic indels, 10,599 large human insertions, and 5,822 large chimpanzee insertions and systematically analyzed the patterns of DNA cleavage intensities in the 200 base pair regions surrounding these indels. Our results show that DNA cleavage intensities close to the start and end points of indels are significantly lower than other regions, for both small human germline and somatic indels and also for mouse small indels. Compared to small indels, the patterns of DNA cleavage intensity around large indels are more complex, and there are two low intensity regions near each end of the indels that are approximately 13 bp apart from each other. Detailed analyses of a subset of indels show that there is slight difference in cleavage intensity distribution between insertion indels and deletion indels that could be contributed by their respective enrichment of different repetitive elements. These results will provide new insight into indel generation mechanisms.

Antimalarial, antimicrobial, cytotoxic, DNA interaction and SOD like activities of tetrahedral copper(II) complexes

The mononuclear copper(II) complexes with P, O-donor ligand and different fluoroquinolones have been synthesized and characterized by elemental analysis, electronic spectra, TGA, EPR, FT-IR and LC-MS spectroscopy. An antimicrobial efficiency of the complexes has been tested against five different microorganisms in terms of minimum inhibitory concentration (MIC) and displays very good antimicrobial activity. The binding strength and binding mode of the complexes with Herring Sperm DNA (HS DNA) have been investigated by absorption titration and viscosity measurement studies. The studies suggest the classical intercalative mode of DNA binding. Gel electrophoresis assay determines the ability of the complexes to cleave the supercoiled form of pUC19 DNA. Synthesized complexes have been tested for their SOD mimic activity using nonenzymatic NBT/NADH/PMS system and found to have good antioxidant activity. All the complexes show good cytotoxic and in vitro antimalarial activities.

Synthesis, DNA binding profile and DNA cleavage pathway of divalent metal complexes

Divalent metal complexes of dipyridylamine ligand with an anthracene moiety induced considerable oxidative DNA cleavage in the presence hydrogen peroxide and dioxygen.

Progress and outlook in structural biology of large viral RNAs

The field of viral molecular biology has reached a precipice for which pioneering studies on the structure of viral RNAs are beginning to bridge the gap. It has become clear that viral genomic RNAs are not simply carriers of hereditary information, but rather are active players in many critical stages during replication. Indeed, functions such as cap-independent translation initiation mechanisms are, in some cases, primarily driven by RNA structural determinants. Other stages including reverse transcription initiation in retroviruses, nuclear export and viral packaging are specifically dependent on the proper 3-dimensional folding of multiple RNA domains to recruit necessary viral and host factors required for activity. Furthermore, a large-scale conformational change within the 5'-untranslated region of HIV-1 has been proposed to regulate the temporal switch between viral protein synthesis and packaging. These RNA-dependent functions are necessary for replication of many human disease-causing viruses such as severe acute respiratory syndrome (SARS)-associated coronavirus, West Nile virus, and HIV-1. The potential for antiviral development is currently hindered by a poor understanding of RNA-driven molecular mechanisms, resulting from a lack of structural information on large RNAs and ribonucleoprotein complexes. Herein, we describe the recent progress that has been made on characterizing these large RNAs and provide brief descriptions of the techniques that will be at the forefront of future advances. Ongoing and future work will contribute to a more complete understanding of the lifecycles of retroviruses and RNA viruses and potentially lead to novel antiviral strategies.

Chemical probing of RNA with the hydroxyl radical at single-atom resolution

While hydroxyl radical cleavage is widely used to map RNA tertiary structure, lack of mechanistic understanding of strand break formation limits the degree of structural insight that can be obtained from this experiment. Here, we determine how individual ribose hydrogens of sarcin/ricin loop RNA participate in strand cleavage. We find that substituting deuterium for hydrogen at a ribose 5'-carbon produces a kinetic isotope effect on cleavage; the major cleavage product is an RNA strand terminated by a 5'-aldehyde. We conclude that hydroxyl radical abstracts a 5'-hydrogen atom, leading to RNA strand cleavage. We used this approach to obtain structural information for a GUA base triple, a common tertiary structural feature of RNA. Cleavage at U exhibits a large 5' deuterium kinetic isotope effect, a potential signature of a base triple. Others had noted a ribose-phosphate hydrogen bond involving the G 2'-OH and the U phosphate of the GUA triple, and suggested that this hydrogen bond contributes to backbone rigidity. Substituting deoxyguanosine for G, to eliminate this hydrogen bond, results in a substantial decrease in cleavage at G and U of the triple. We conclude that this hydrogen bond is a linchpin of backbone structure around the triple.

Independent generation and reactivity of uridin-2'-yl radical

The uridin-2'-yl radical (1) has been proposed as an intermediate during RNA oxidation. However, its reactivity has not been thoroughly studied due to the complex conditions under which it is typically generated. The uridin-2'-yl radical was independently generated from a benzyl ketone (2a) via Norrish type I photocleavage upon irradiation at λmax = 350 nm. Dioxygen and β-mercaptoethanol are unable to compete with loss of uracil from 1 in phosphate buffer. Thiol trapping competes with uracil fragmentation in less polar solvent conditions. This is ascribed mostly to a reduction in the rate constant for uracil elimination in the less polar solvent. Hydrogen atom transfer to 1 from β-mercaptoethanol occurs exclusively from the α-face to produce arabinouridine. Mass balances range from 72 to 95%. Furthermore, the synthesis of 2a is amenable to formation of the requisite phosphoramidite for solid-phase oligonucleotide synthesis. This and the fidelity with which the urdin-2'-yl radical is generated from 2a suggest that this precursor should be useful for studying the radical's reactivity in synthetic oligonucleotides.

RNA tertiary structure analysis by 2'-hydroxyl molecular interference

We introduce a melded chemical and computational approach for probing and modeling higher-order intramolecular tertiary interactions in RNA. 2'-Hydroxyl molecular interference (HMX) identifies nucleotides in highly packed regions of an RNA by exploiting the ability of bulky adducts at the 2'-hydroxyl position to disrupt overall RNA structure. HMX was found to be exceptionally selective for quantitative detection of higher-order and tertiary interactions. When incorporated as experimental constraints in discrete molecular dynamics simulations, HMX information yielded accurate three-dimensional models, emphasizing the power of molecular interference to guide RNA tertiary structure analysis and fold refinement. In the case of a large, multidomain RNA, the Tetrahymena group I intron, HMX identified multiple distinct sets of tertiary structure interaction groups in a single, concise experiment.

Interstrand cross-link and bioconjugate formation in RNA from a modified nucleotide

RNA oligonucleotides containing a phenyl selenide derivative of 5-methyluridine were chemically synthesized by solid-phase synthesis. The phenyl selenide is rapidly converted to an electrophilic, allylic phenyl seleneate under mild oxidative conditions. The phenyl seleneate yields interstrand cross-links when part of a duplex and is useful for synthesizing oligonucleotide conjugates. Formation of the latter is illustrated by reaction of an oligonucleotide containing the phenyl selenide with amino acids in the presence of mild oxidant. The products formed are analogous to those observed in tRNA that are believed to be formed posttranslationally via a biosynthetic intermediate that is chemically homologous to the phenyl seleneate.

Prediction of nucleosome rotational positioning in yeast and human genomes based on sequence-dependent DNA anisotropy

BACKGROUND

An organism's DNA sequence is one of the key factors guiding the positioning of nucleosomes within a cell's nucleus. Sequence-dependent bending anisotropy dictates how DNA is wrapped around a histone octamer. One of the best established sequence patterns consistent with this anisotropy is the periodic occurrence of AT-containing dinucleotides (WW) and GC-containing dinucleotides (SS) in the nucleosomal locations where DNA is bent in the minor and major grooves, respectively. Although this simple pattern has been observed in nucleosomes across eukaryotic genomes, its use for prediction of nucleosome positioning was not systematically tested.

RESULTS

We present a simple computational model, termed the W/S scheme, implementing this pattern, without using any training data. This model accurately predicts the rotational positioning of nucleosomes both in vitro and in vivo, in yeast and human genomes. About 65 - 75% of the experimentally observed nucleosome positions are predicted with the precision of one to two base pairs. The program is freely available at http://people.rit.edu/fxcsbi/WS_scheme/. We also introduce a simple and efficient way to compare the performance of different models predicting the rotational positioning of nucleosomes.

CONCLUSIONS

This paper presents the W/S scheme to achieve accurate prediction of rotational positioning of nucleosomes, solely based on the sequence-dependent anisotropic bending of nucleosomal DNA. This method successfully captures DNA features critical for the rotational positioning of nucleosomes, and can be further improved by incorporating additional terms related to the translational positioning of nucleosomes in a species-specific manner.

Finding and characterizing the complexes of drug like molecules with quadruplex DNA: combined use of an enhanced hydroxyl radical cleavage protocol and NMR

Structural information on the complexes of drug like molecules with quadruplex DNAs can aid the development of therapeutics and research tools that selectively target specific quadruplex DNAs. Screening can identify candidate molecules that require additional evaluation. An enhanced hydroxyl radical cleavage protocol is demonstrated that can efficiently provide structural information on the complexes of the candidate molecules with quadruplex DNA. NMR methods have been used to offer additional structural information about the complexes as well as validate the results of the hydroxyl radical approach. This multi-step protocol has been demonstrated on complexes of the chair type quadruplex formed by the thrombin binding aptamer, d(GGTTGGTGTGGTTGG). The hydroxyl radical results indicate that NSC 176319, Cain’s quinolinium that was found by screening, exhibits selective binding to the two TT loops. The NMR results are consistent with selective disruption of the hydrogen bonding between T4 and T13 as well as unstacking of these residues from the bottom quartet. Thus, the combination of screening, hydroxyl radical footprinting and NMR can find new molecules that selectively bind to quadruplex DNAs as well as provide structural information about their complexes.

An RNA aptamer possessing a novel monovalent cation-mediated fold inhibits lysozyme catalysis by inhibiting the binding of long natural substrates

RNA aptamers are being developed as inhibitors of macromolecular and cellular function, diagnostic tools, and potential therapeutics. Our understanding of the physical nature of this emerging class of nucleic acid-protein complexes is limited; few atomic resolution structures have been reported for aptamers bound to their protein target. Guided by chemical mapping, we systematically minimized an RNA aptamer (Lys1) selected against hen egg white lysozyme. The resultant 59-nucleotide compact aptamer (Lys1.2minE) retains nanomolar binding affinity and the ability to inhibit lysozyme's catalytic activity. Our 2.0-Å crystal structure of the aptamer-protein complex reveals a helical stem stabilizing two loops to form a protein binding platform that binds lysozyme distal to the catalytic cleft. This structure along with complementary solution analyses illuminate a novel protein-nucleic acid interface; (1) only 410 Å(2) of solvent accessible surface are buried by aptamer binding; (2) an unusually small fraction (∼18%) of the RNA-protein interaction is electrostatic, consistent with the limited protein phosphate backbone contacts observed in the structure; (3) a single Na(+) stabilizes the loops that constitute the protein-binding platform, and consistent with this observation, Lys1.2minE-lysozyme complex formation takes up rather than displaces cations at low ionic strength; (4) Lys1.2minE inhibits catalysis of large cell wall substrates but not catalysis of small model substrates; and (5) the helical stem of Lys1.2minE can be shortened to four base pairs (Lys1.2minF) without compromising binding affinity, yielding a 45-nucleotide aptamer whose structure may be an adaptable protein binding platform.