Abstract
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.