HiTRACE: high-throughput robust analysis for capillary electrophoresis

MOTIVATION: Capillary electrophoresis (CE) of nucleic acids is a workhorse technology underlying high-throughput genome analysis and large-scale chemical mapping for nucleic acid structural inference. Despite the wide availability of CE-based instruments, there remain challenges in leveraging their full power for quantitative analysis of RNA and DNA structure, thermodynamics and kinetics. In particular, the slow rate and poor automation of available analysis tools have bottlenecked a new generation of studies involving hundreds of CE profiles per experiment. RESULTS: We propose a computational method called high-throughput robust analysis for capillary electrophoresis (HiTRACE) to automate the key tasks in large-scale nucleic acid CE analysis, including the profile alignment that has heretofore been a rate-limiting step in the highest throughput experiments. We illustrate the application of HiTRACE on 13 datasets representing 4 different RNAs, 3 chemical modification strategies and up to 480 single mutant variants; the largest datasets each include 87 360 bands. By applying a series of robust dynamic programming algorithms, HiTRACE outperforms prior tools in terms of alignment and fitting quality, as assessed by measures including the correlation between quantified band intensities between replicate datasets. Furthermore, while the smallest of these datasets required 7-10 h of manual intervention using prior approaches, HiTRACE quantitation of even the largest datasets herein was achieved in 3-12 min. The HiTRACE method, therefore, resolves a critical barrier to the efficient and accurate analysis of nucleic acid structure in experiments involving tens of thousands of electrophoretic bands.     

ShapeFinder: a software system for high-throughput quantitative analysis of nucleic acid reactivity information resolved by capillary electrophoresis

Analysis of the long-range architecture of RNA is a challenging experimental and computational problem. Local nucleotide flexibility, which directly reports underlying base pairing and tertiary interactions in an RNA, can be comprehensively assessed at single nucleotide resolution using high-throughput selective 2'-hydroxyl acylation analyzed by primer extension (hSHAPE). hSHAPE resolves structure-sensitive chemical modification information by high-resolution capillary electrophoresis and typically yields quantitative nucleotide flexibility information for 300-650 nucleotides (nt) per experiment. The electropherograms generated in hSHAPE experiments provide a wealth of structural information; however, significant algorithmic analysis steps are required to generate quantitative and interpretable data. We have developed a set of software tools called ShapeFinder to make possible rapid analysis of raw sequencer data from hSHAPE, and most other classes of nucleic acid reactivity experiments. The algorithms in ShapeFinder (1) convert measured fluorescence intensity to quantitative cDNA fragment amounts, (2) correct for signal decay over read lengths extending to 600 nts or more, (3) align reactivity data to the known RNA sequence, and (4) quantify per nucleotide reactivities using whole-channel Gaussian integration. The algorithms and user interface tools implemented in ShapeFinder create new opportunities for tackling ambitious problems involving high-throughput analysis of structure-function relationships in large RNAs.     

Comparison of rRNA cleavage by complementary 1,10-phenanthroline-Cu(II)- and EDTA-Fe(II)-derivatized oligonucleotides

The chemical nucleases 1,10-phenanthroline-Cu(II) and EDTA-Fe(II), have proven to be valuable tools for structural analysis of nucleic acids. Both have found applications in footprinting and directed proximity studies of DNA and RNA. Derivatives of each that provide for tethering to nucleic acid or protein are commercially available, allowing their widespread use for structural analysis of macromolecules. Although their applications are somewhat overlapping, differences in their cleavage mechanisms and chemical properties allow them to provide distinct and complementary structural information. The purpose of this study is to compare directly the cleavage patterns of tethered 1,10-phenanthroline-Cu(II) and EDTA-Fe(II) complexes within a similar experimental system. Here, the region surrounding nucleotide 1400 of 16S rRNA from Escherichia coli serves as a substrate for chemical cleavage directed by a derivatized complementary oligonucleotide. This region of rRNA is known to be involved in the decoding of mRNA during translation. The results of this study provide evidence in support of the mechanistic differences previously established for EDTA-Fe(II) and 1,10-phenathroline-Cu(II). The delocalized cleavage envelope produced by EDTA-Fe(II) cleavage suggests the involvement of a diffusible reactive species. On the other hand, rRNA cleavage induced by the tethered 1,10-phenanthroline-Cu(II) complex appears localized to the proximity of the chemical nuclease under normal conditions, although the production of an unknown diffusible species appears to occur during long reaction times.     

Monitoring structural changes in nucleic acids with single residue spatial and millisecond time resolution by quantitative hydroxyl radical footprinting

Hydroxyl radical (.OH) footprinting provides comprehensive site-specific quantitative information about the structural changes associated with macromolecular folding, interactions and ligand binding. 'Fast Fenton' footprinting is a laboratory-based method for time-resolved .OH footprinting capable of millisecond time resolution readily applicable to DNA and RNA. This protocol utilizes inexpensive chemical reagents (H2O2, Fe(NH4)2(SO4)2, EDTA, thiourea or ethanol) and widely available quench-flow mixers to reveal transient, often short-lived, intermediate states of complex biochemical processes. We describe a protocol developed to study RNA folding that can be readily tailored to particular applications. Once familiar with quench-flow mixer operation and its calibration, nucleic acid labeling and the conduct of a dose-response experiment, a single kinetic experiment of 30 time points takes about 1 h to perform. Sample processing and separation of the .OH reaction products takes several hours. Data analysis can take 45 min to several weeks depending on the depth of analysis conducted.     

The proto-Nucleic Acid Builder: a software tool for constructing nucleic acid analogs

The helical structures of DNA and RNA were originally revealed by experimental data. Likewise, the development of programs for modeling these natural polymers was guided by known structures. These nucleic acid polymers represent only two members of a potentially vast class of polymers with similar structural features, but that differ from DNA and RNA in the backbone or nucleobases. Xeno nucleic acids (XNAs) incorporate alternative backbones that affect the conformational, chemical, and thermodynamic properties of XNAs. Given the vast chemical space of possible XNAs, computational modeling of alternative nucleic acids can accelerate the search for plausible nucleic acid analogs and guide their rational design. Additionally, a tool for the modeling of nucleic acids could help reveal what nucleic acid polymers may have existed before RNA in the early evolution of life. To aid the development of novel XNA polymers and the search for possible pre-RNA candidates, this article presents the proto-Nucleic Acid Builder (https://github.com/GT-NucleicAcids/pnab), an open-source program for modeling nucleic acid analogs with alternative backbones and nucleobases. The torsion-driven conformation search procedure implemented here predicts structures with good accuracy compared to experimental structures, and correctly demonstrates the correlation between the helical structure and the backbone conformation in DNA and RNA.     

Fluorescently labeled oligonucleotide extension: a rapid and quantitative protocol for primer extension

Identification of nucleotides used for RNA chain initiation or for contacting DNA binding proteins is basic to our understanding of gene regulation. Normally, a radioactive primer is used to copy RNA or DNA. The polymerase extension stops at free ends of mRNA (as in promoter mapping) or at the position of template cleavage or modification (as in footprinting). The locations of these positions are then analyzed by polyacrylamide gel electrophoresis. These analyses have been improved using fluorescently labeled primers and commonly available DNA sequencing machines. The protocol, which we call fluorescently labeled oligonucleotide extension (FLOE), eliminates the need for handling radioactivity and polyacrylamide. The DNA sequencer delivers data as a "trace" that is ready for quantification, which eliminates the need to trace gels separately. The data analysis is further improved by new software, Scanalyze, which we present here. We demonstrate that by using promoter mapping and footprinting, FLOE shortens experimental time, extends the stretch of analyzable sequence, and simplifies quantification compared to radioactive methods and is as sensitive in terms of detecting templates.     

Radical-generating coordination complexes as tools for rapid and effective fragmentation and fluorescent labeling of nucleic acids for microchip hybridization

DNA microchip technology is a rapid, high-throughput method for nucleic acid hybridization reactions. This technology requires random fragmentation and fluorescent labeling of target nucleic acids prior to hybridization. Radical-generating coordination complexes, such as 1,10-phenanthroline-Cu(II) (OP-Cu) and Fe(II)-EDTA (Fe-EDTA), have been commonly used as sequence nonspecific "chemical nucleases" to introduce single-strand breaks in nucleic acids. Here we describe a new method based on these radical-generating complexes for random fragmentation and labeling of both single- and double-stranded forms of RNA and DNA. Nucleic acids labeled with the OP-Cu and the Fe-EDTA protocols revealed high hybridization specificity in hybridization with DNA microchips containing oligonucleotide probes selected for identification of 16S rRNA sequences of the Bacillus group microorganisms.We also demonstrated cDNA- and cRNA-labeling and fragmentation with this method. Both the OP-Cu and Fe-EDTA fragmentation and labeling procedures are quick and inexpensive compared to other commonly used methods. A column-based version of the described method does not require centrifugation and therefore is promising for the automation of sample preparations in DNA microchip technology as well as in other nucleic acid hybridization studies.     

Polymerizing immobilization of acrylamide-modified nucleic acids and its application

The immobilization of nucleic acids on solid supports has been widely used in the detection of DNA and other biomolecules in sensor technology. Because three dimensional (3-D) hydrogel matrixes offer significant advantages for capturing probes over more conventional two dimensional (2-D) rigid substrates and the ability to provide a solution-mimicking environment, they are becoming increasingly attractive as desired supports for bio-analysis. Acrylamide-modified nucleic acids and acrylamide monomers being polymerized directly to immobilize nucleic acids is only one-step chemical process which is not interfered by exterior surroundings, and the 3-D polyacrylamide gel fabricated by this method is not required to be activated by some labile chemical treatments. Moreover, the attachment is extremely stable to withstand the cycling process involved in the polymerase chain reaction (PCR). In this paper, the development of polymerizing immobilization of acrylamide-modified nucleic acids is reviewed, and its applications in DNA sequence high-throughput analysis including mutation analysis and the whole genome sequencing are summarized.     

Quantitative analysis of nucleic acids--the last few years of progress

DNA and RNA quantifications are widely used in biological and biomedical research. In the last ten years, many technologies have been developed to enable automated and high-throughput analyses. In this review, we first give a brief overview of how DNA and RNA quantifications are carried out. Then, five technologies (microarrays, SAGE, differential display, real time PCR and real competitive PCR) are introduced, with an emphasis on how these technologies can be applied and what their limitations are. The technologies are also evaluated in terms of a few key aspects of nucleic acids quantification such as accuracy, sensitivity, specificity, cost and throughput.     

Sharing and archiving nucleic acid structure mapping data

Nucleic acids are particularly amenable to structural characterization using chemical and enzymatic probes. Each individual structure mapping experiment reveals specific information about the structure and/or dynamics of the nucleic acid. Currently, there is no simple approach for making these data publically available in a standardized format. We therefore developed a standard for reporting the results of single nucleotide resolution nucleic acid structure mapping experiments, or SNRNASMs. We propose a schema for sharing nucleic acid chemical probing data that uses generic public servers for storing, retrieving, and searching the data. We have also developed a consistent nomenclature (ontology) within the Ontology of Biomedical Investigations (OBI), which provides unique identifiers (termed persistent URLs, or PURLs) for classifying the data. Links to standardized data sets shared using our proposed format along with a tutorial and links to templates can be found at http://snrnasm.bio.unc.edu.     

Fast Fenton footprinting: a laboratory-based method for the time-resolved analysis of DNA, RNA and proteins

‘Footprinting’ describes assays in which ligand binding or structure formation protects polymers such as nucleic acids and proteins from either cleavage or modification; footprinting allows the accessibility of individual residues to be mapped in solution. Equilibrium and time-dependent footprinting links site-specific structural information with thermodynamic and kinetic transitions. The hydroxyl radical (·OH) is a particularly valuable footprinting probe by virtue of it being among the most reactive of chemical oxidants; it reports the solvent accessibility of reactive sites on macromolecules with as fine as a single residue resolution. A novel method of millisecond time-resolved ·OH footprinting has been developed based on the Fenton reaction, Fe(II) + H₂O₂ → Fe(III) + ·OH + OH⁻. This method can be implemented in laboratories using widely available three-syringe quench flow mixers and inexpensive reagents to study local changes in the solvent accessibility of DNA, RNA and proteins associated with their biological function.     

RNA footprinting analysis using ion pair reverse phase liquid chromatography

Hydroxyl radical footprinting is a powerful technique often employed in characterization of the tertiary interactions between proteins and nucleic acids. Following the generation of a nucleic acid "ladder" either by chemical or enzymatic reactions, the radiolabeled products are traditionally separated by denaturing gel electrophoresis and further quantified by phosphorimaging techniques. Here we report the use of ion pair reverse phase liquid chromatography to analyze the products of an RNA footprinting reaction using fluorescently labeled RNA molecules. This technique offers several advantages over existing procedures, including rapid analysis, automation, and direct quantification of the cleavage products without the need to employ radiolabeling. To illustrate the resolving power of this technique, we have analyzed the products of base hydrolysis, generated from a fluorescently labeled RNA molecule and have subsequently used this method to define the solvent accessibility of the substrate strand as it docks with the hairpin ribozyme.     

New photoactivatable structural and affinity probes of RNAs: specific features and applications for mapping of spermine binding sites in yeast tRNA(Asp) and interaction of this tRNA with yeast aspartyl-tRNA synthetase.

Aryldiazonium salts are shown to be useful as phototriggered structural probes for RNA mapping as well as for footprinting of RNA/protein interaction. In particular the yeast tRNA(Asp)/aspartyl-tRNA synthetase complex is shown to involve the variable loop face and the concave side of the L-shaped nucleic acid bound to a lipophilic area of the enzyme. When chemically linked to spermine, the photoactive group cleaves RNA at polyamine binding sites; 3-4 spermines have been located in the tRNA(Asp), stabilizing the central part of the molecule in regions where two ribose-phosphate strands are close to each other.     

Mutations in 5S DNA and 5S RNA have different effects on the binding of Xenopus transcription factor IIIA

The effects on TFIIIA binding affinity of a series of substitution mutations in the Xenopus laevis oocyte 5S RNA gene were quantified. These data indicate that TFIIIA binds specifically to 5S DNA by forming sequence-specific contacts with three discrete sites located within the classical A and C boxes and the intermediate element of the internal control region. Substitution of the nucleotide sequence at any of the three sites significantly reduces TFIIIA binding affinity, with a 100-fold reduction observed for substitutions in the box C subregion. These results are consistent with a direct interaction of TFIIIA with specific base pairs within the major groove of the DNA. A comparison of the TFIIIA binding data for the same mutations expressed in 5S RNA indicates that the protein does not make any strong sequence-specific contacts with the RNA. Although the protein footprinting sites on the 5S DNA and 5S RNA are coincident, nucleotide substitutions in 5S RNA which moderately reduce TFIIIA binding affinity do not correspond at all to the three specific TFIIIA interaction sites within the gene. The implications of these results for models which attempt to reconcile the DNA and RNA binding activities of TFIIIA by proposing a common structural motif for the two nucleic acids are discussed.     

Thermophoretic melting curves quantify the conformation and stability of RNA and DNA

Measuring parameters such as stability and conformation of biomolecules, especially of nucleic acids, is important in the field of biology, medical diagnostics and biotechnology. We present a thermophoretic method to analyse the conformation and thermal stability of nucleic acids. It relies on the directed movement of molecules in a temperature gradient that depends on surface characteristics of the molecule, such as size, charge and hydrophobicity. By measuring thermophoresis of nucleic acids over temperature, we find clear melting transitions and resolve intermediate conformational states. These intermediate states are indicated by an additional peak in the thermophoretic signal preceding most melting transitions. We analysed single nucleotide polymorphisms, DNA modifications, conformational states of DNA hairpins and microRNA duplexes. The method is validated successfully against calculated melting temperatures and UV absorbance measurements. Interestingly, the methylation of DNA is detected by the thermophoretic amplitude even if it does not affect the melting temperature. In the described setup, thermophoresis is measured all-optical in a simple setup using a reproducible capillary format with only 250 nl probe consumption. The thermophoretic analysis of nucleic acids shows the technique's versatility for the investigation of nucleic acids relevant in cellular processes like RNA interference or gene silencing.     

Semi-automated, single-band peak-fitting analysis of hydroxyl radical nucleic acid footprint autoradiograms for the quantitative analysis of transitions

Hydroxyl radical footprinting can probe the solvent accessibility of the ribose moiety of the individual nucleotides of DNA and RNA. Semi-automated analytical tools are presented for the quantitative analyses of nucleic acid footprint transitions in which processes such as folding or ligand binding are followed as a function of time or ligand concentration. Efficient quantitation of the intensities of the electrophoretic bands comprising the footprinting reaction products is achieved by fitting a series of Lorentzian curves to line profiles obtained from gels utilizing sequentially relaxed constraints consistent with electrophoretic mobility. An automated process of data ‘standardization’ has been developed that corrects for differences in the loading amounts in the electrophoresis. This process enhances the accuracy of the derived transitions and makes generating them easier. Together with visualization of the processed footprinting in false-color two-dimensional maps, DNA and RNA footprinting data can be accurately, precisely and efficiently processed allowing transitions to be objectively and comprehensively analyzed. The utility of this new analysis approach is illustrated by its application to the ion-meditated folding of a large RNA molecule.