The Chirality of a Four-Way Helical Junction in RNA

CHSalign: A Web Server That Builds upon Junction-Explorer and RNAJAG for Pairwise Alignment of RNA Secondary Structures with Coaxial Helical Stacking

RNA junctions are important structural elements of RNA molecules. They are formed when three or more helices come together in three-dimensional space. Recent studies have focused on the annotation and prediction of coaxial helical stacking (CHS) motifs within junctions. Here we exploit such predictions to develop an efficient alignment tool to handle RNA secondary structures with CHS motifs. Specifically, we build upon our Junction-Explorer software for predicting coaxial stacking and RNAJAG for modelling junction topologies as tree graphs to incorporate constrained tree matching and dynamic programming algorithms into a new method, called CHSalign, for aligning the secondary structures of RNA molecules containing CHS motifs. Thus, CHSalign is intended to be an efficient alignment tool for RNAs containing similar junctions. Experimental results based on thousands of alignments demonstrate that CHSalign can align two RNA secondary structures containing CHS motifs more accurately than other RNA secondary structure alignment tools. CHSalign yields a high score when aligning two RNA secondary structures with similar CHS motifs or helical arrangement patterns, and a low score otherwise. This new method has been implemented in a web server, and the program is also made freely available, at http://bioinformatics.njit.edu/CHSalign/.

Modeling and design by hierarchical natural moves

We develop a unique algorithm implemented in the program MOSAICS (Methodologies for Optimization and Sampling in Computational Studies) that is capable of nanoscale modeling without compromising the resolution of interest. This is achieved by modeling with customizable hierarchical degrees of freedom, thereby circumventing major limitations of conventional molecular modeling. With the emergence of RNA-based nanotechnology, large RNAs in all-atom representation are used here to benchmark our algorithm. Our method locates all favorable structural states of a model RNA of significant complexity while improving sampling accuracy and increasing speed many fold over existing all-atom RNA modeling methods. We also modeled the effects of sequence mutations on the structural building blocks of tRNA-based nanotechnology. With its flexibility in choosing arbitrary degrees of freedom as well as in allowing different all-atom energy functions, MOSAICS is an ideal tool to model and design biomolecules of the nanoscale.

3D maps of RNA interhelical junctions

More than 50% of RNA secondary structure is estimated to be A-form helices, which are linked together by various junctions. Here we describe a protocol for computing three interhelical Euler angles describing the relative orientation of helices across RNA junctions. 5' and 3' helices, H1 and H2, respectively, are assigned based on the junction topology. A reference canonical helix is constructed using an appropriate molecular builder software consisting of two continuous idealized A-form helices (iH1 and iH2) with helix axis oriented along the molecular Z-direction running toward the positive direction from iH1 to iH2. The phosphate groups and the carbon and oxygen atoms of the sugars are used to superimpose helix H1 of a target interhelical junction onto the corresponding iH1 of the reference helix. A copy of iH2 is then superimposed onto the resulting H2 helix to generate iH2'. A rotation matrix R is computed, which rotates iH2' into iH2 and expresses the rotation parameters in terms of three Euler angles α(h), β(h) and γ(h). The angles are processed to resolve a twofold degeneracy and to select an overall rotation around the axis of the reference helix. The three interhelical Euler angles define clockwise rotations around the 5' (-γ(h)) and 3' (α(h)) helices and an interhelical bend angle (β(h)). The angles can be depicted graphically to provide a 'Ramachandran'-type view of RNA global structure that can be used to identify unusual conformations as well as to understand variations due to changes in sequence, junction topology and other parameters.

A quantitative measure of chirality inside nucleic acid databank

We show the capability of a chirality index (Pietropaolo et al., Proteins 2008;70:667-677) to investigate nucleic acid structures because of its high sensitivity to helical conformations. By analyzing selected structures of DNA and RNA, we have found that sequences rich in cytosine and guanine have a tendency to left-handed chirality, in contrast to regions rich in adenine or thymine which show strong negative, right-handed, chirality values. We also analyze RNA structures, where specific loops and hairpin motifs are characterized by a well-defined chirality value. We find that in nucleosome the chirality is exalted, whereas in ribosome it is reduced. Our results illustrate the sensitivity of this descriptor for nucleic acid conformations.

Comparative gel electrophoresis analysis of helical junctions in RNA

Comparative gel electrophoresis provides information on the relative angles subtended between helical arms at a branchpoint in RNA. It is based upon the comparison of electrophoretic mobility in polyacrylamide gels of species containing two long arms, with the remaining one(s) being significantly shorter. Although the method currently lacks a really well-established basis of physical theory, it is very powerful, yet simple to apply. It has had a number of significant successes in RNA, DNA and DNA-protein complexes, and in all cases to date the results have stood the test of time and eventual comparison with crystallographic analysis.

Analysis of branched nucleic acid structure using comparative gel electrophoresis

Electrophoresis in polyacrylamide gels provides a simple yet powerful means of analyzing the relative disposition of helical arms in branched nucleic acids. The electrophoretic mobility of DNA or RNA with a central discontinuity is determined by the angle subtended between the arms radiating from the branchpoint. In a multi-helical branchpoint, comparative gel electrophoresis can provide a relative measure of all the inter-helical angles and thus the shape and symmetry of the molecule. Using the long-short arm approach, the electrophoretic mobility of all the species with two helical arms that are longer than all others is compared. This can be done as a function of conditions, allowing the analysis of ion-dependent folding of branched DNA and RNA species. Notable successes for the technique include the four-way (Holliday) junction in DNA and helical junctions in functionally significant RNA species such as ribozymes. Many of these structures have subsequently been proved correct by crystallography or other methods, up to 10 years later in the case of the Holliday junction. Just as important, the technique has not failed to date. Comparative gel electrophoresis can provide a window on both fast and slow conformational equilibria such as conformer exchange in four-way DNA junctions. But perhaps the biggest test of the approach has been to deduce the structures of complexes of four-way DNA junctions with proteins. Two recent crystallographic structures show that the global structures were correctly deduced by electrophoresis, proving the worth of the method even in these rather complex systems. Comparative gel electrophoresis is a robust method for the analysis of branched nucleic acids and their complexes.