Statistical analysis of RNA backbone

Local conformation is an important determinant of RNA catalysis and binding. The analysis of RNA conformation is particularly difficult due to the large number of degrees of freedom (torsion angles) per residue. Proteins, by comparison, have many fewer degrees of freedom per residue. In this work, we use and extend classical tools from statistics and signal processing to search for clusters in RNA conformational space. Results are reported both for scalar analysis, where each torsion angle is separately studied, and for vectorial analysis, where several angles are simultaneously clustered. Adapting techniques from vector quantization and clustering to the RNA structure, we find torsion angle clusters and RNA conformational motifs. We validate the technique using well-known conformational motifs, showing that the simultaneous study of the total torsion angle space leads to results consistent with known motifs reported in the literature and also to the finding of new ones.     

Characterization of the solution conformations of unbound and Tat peptide-bound forms of HIV-1 TAR RNA.

Basic peptides from the carboxy terminus of the HIV-1 Tat protein bind to the apical stem-loop region of TAR RNA with high affinity and moderate specificity. The conformations of the unbound and 24 residue Tat peptide (Tfr24)-bound forms of TAR RNA have been characterized by NMR spectroscopy. The unbound form of TAR exists in major and minor forms having different trinucleotide bulge conformations. A specific TAR RNA conformational change is observed upon complex formation with Tfr24, consisting of coaxial stacking of helical stems and base triple formation. A U23-A27-U38 base triple is proposed based on exchangeable proton NMR data, where U23 forms a base pair with A27 in the major groove. No evidence for base triple formation was found for Tat peptides in which lysine residues are extensively substituted for arginine.     

''A'' forms of RNAs in single strands, duplexes and RNA-DNA hybrids.

Helical parameters have been calculated for the 'A' form minimum energy conformations of ApA, CpC, GpG, UpU, GpC and UpA. The helix geometries are base sequence dependent. The single strands are narrower and more tightly wound than that duplex RNA-11 form. 9-12 kcal./mole are needed to convert these single strands to the RNA-11 conformation. However, in some sequences other 'A' type conformers capable of complementary base pairing may be formed at lower energetic cost. There is substantially more base stacking in the calculated single strands than in the RNA-11 conformation. Calculated intrastrand base stacking energies reflect these differences, and also are sequence dependent. The 'A' form RNA subunits differ from the analogous DNAs in possessing a larger rise per residue, needed to accomodate the 2'-OH. RNA-DNA hybrids are consequently more likely to be in the 'A-RNA than in the 'A'-DNA conformation, although the base sequence determines the extent of the preference.     

Low temperature solution structures and base pair stacking of double helical d(CGTACG)(2)

Solution structures and base pair stacking of a self- complementary DNA hexamer d(CGTACG)(2) have been studied at 5, 10 and 15 degrees C, respectively. The stacking interactions among the center base pair steps of the DNA duplex are found to improve when the terminal base pairs became less stable due to end fraying. A new structural quantity, the stacking sum (Sigma(s)), is introduced to indicate small changes in the stacking overlaps between base pairs. The improvements in the stacking overlaps to maintain the double helical conformation are probably the cause for the observed temperature dependent structural changes in double helical DNA molecule. A detailed analysis of the helical parameters, backbone torsion angles, base orientations and sugar conformations of these structures has been performed.     

The dynamic structural basis of differential enhancement of conformational stability by 5'- and 3'-dangling ends in RNA

Unpaired bases at the end of an RNA duplex (dangling ends) can stabilize the core duplex in a sequence-dependent manner and are important determinants of RNA folding, recognition, and functions. Using 2-aminopurine as a dangling end purine base, we have employed femtosecond time-resolved fluorescence spectroscopy, combined with UV optical melting, to quantitatively investigate the physical and structural nature of the stacking interactions between the dangling end bases and the terminal base pairs. A 3'-dangling purine base has a large subpopulation that stacks on the guanine base of the terminal GC or UG pair, either intrastrand or cross-strand depending on the orientation of the pair, thus providing stabilization of different magnitudes. On the contrary, a 5'-dangling purine base only has a marginal subpopulation that stacks on the purine of the same strand (intrastrand) but has little cross-strand stacking. Thus a 5'-dangling purine does not provide significant stabilization. These stacking structures are not static, and a dangling end base samples a range of stacked and unstacked conformations with respect to the terminal base pair. Femtosecond time-resolved anisotropy decay reveals certain hindered base conformational dynamics that occur on the picosecond to nanosecond time scales, which allow the dangling base to sample these substates. When the dangling purine is opposite to a U and is able to form a potential base pair at the end of the duplex, there is an interplay of base stacking and hydrogen-bonding interactions that depends on the orientation of the base pair relative to the adjacent GC pair. By resolving these populations that are dynamically exchanging on fast time scales, we elucidated the correlation between dynamic conformational distributions and thermodynamic stability.     

Raman studies of nucleic acids. XII. Conformations of oligonucleotides and deuterated polynucleotides

Laser Raman spectra of the trinucleoside diphoshate ApApA and dinucleoside phosphates ApU, UpA, GpC, CpG, and GpU are reported and discussed. Assignments of conformationally sensitive frequencies are-facilitated by comparison with spectra reported here of poly(rA), poly(rC), and poly(rU) in deuterium oxide solutions. The significant spectral differences between ApU and UpA, and between GpC and CpG, reveal that the sequence isomers have nonidentical conformations in aqueous solution. In UpA at low temperature the bases are stacked and the backbone conformation is similar to that found in ordered polynucleotide structures and RNA. In ApU no base stacking can be detected and the backbone conformation differs from that found in UpA, both in the orientation of phosphodiester linkages and in the internal conformation of ribose. At the conditions employed neither ApU nor UpA exhibits base pairing in aqueous solutions. In both GpC and CpG the bases are stacked and the phosphodiester conformations are similar to those encountered for UpA and RNA. However, major differences between spectra of GpC and CpG indicate that the geometries of stacking and ribosyl conformations are different. In GpC the Raman data favor the formation of hydrogen bonded dimers containing GC pairs. Protonation of C in GpC is sufficient to eliminate the ordered conformation detected by Raman spectroscopy. Despite the ordered backbone conformation evident in GpU, this dinucleoside apparently contains neither stacked nor hydrogen bonded bases at the conditions employed here. The Raman data also confirm the stacking interactions in ApApA, poly(rA), and poly(rC) but suggest that the backbone conformation in poly(rC) differs qualitatively from that found in most ordered polynucleotide structures and is thermally more stable. The present results demonstrate the sensitivity of the Raman technique to sequence-related structural differences in oligonucleotides and provide additional spectra–structure correlations for future conformational studies of RNA by laser Raman spectroscopy.     

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model.

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A -form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.     

Solution structure of duplex DNA containing the mutagenic lesion 1,N(2)-etheno-2'-deoxyguanine

We have used high-resolution NMR spectroscopy and molecular dynamics simulations to determine the solution structure of DNA containing the genotoxic lesion 1, N (2)-etheno-2'-deoxyguanosine (epsilonG), paired to dC. The NMR data suggest the presence of a major, minimally perturbed structure at neutral pH. NOESY spectra indicate the presence of a right-handed helix with all nucleotides in anti, 2'-deoxyribose conformations within the C2'-endo/C1'-exo range and proper Watson-Crick base pair alignments outside the lesion site. The epsilonG residue remains deeply embedded inside the helix and stacks between the flanking base pairs. The lesion partner dC is extrahelical and is located in the minor groove of the duplex, where it is highly exposed to solvent. Upon acidification of the sample, a second conformation at the lesion site of the duplex emerges, with protonation of the lesion partner dC and possible formation of a Hoogsteen base pair. Restrained molecular dynamics simulations of the neutral-pH structure generated a set of three-dimensional models that show epsilonG inside the helix, where the lesion is stabilized by stacking interactions with flanking bases but without participating in hydrogen bonding. The lesion counterbase dC is displaced in the minor groove of the duplex where it can form a hydrogen bond with the sugar O4' atom of a residue 2 bp away.     

Efficient in silico exploration of RNA interhelical conformations using Euler angles and WExplore

HIV-1 TAR RNA is a two-helix bulge motif that plays a critical role in HIV viral replication and is an important drug target. However, efforts at designing TAR inhibitors have been challenged by its high degree of structural flexibility, which includes slow large-amplitude reorientations of its helices with respect to one another. Here, we use the recently introduced algorithm WExplore in combination with Euler angles to achieve unprecedented sampling of the TAR conformational ensemble. Our ensemble achieves similar agreement with experimental NMR data when compared with previous TAR computational studies, and is generated at a fraction of the computational cost. It clearly emerges from configuration space network analysis that the intermittent formation of the A22-U40 base pair acts as a reversible switch that enables sampling of interhelical conformations that would otherwise be topologically disallowed. We find that most previously determined ligand-bound structures are found in similar location in the network, and we use a sample-and-select approach to guide the construction of a set of novel conformations which can serve as the basis for future drug development efforts. Collectively, our findings demonstrate the utility of WExplore in combination with suitable order parameters as a method for exploring RNA conformational space.     

Nearest neighbor parameters for Watson-Crick complementary heteroduplexes formed between 2'-O-methyl RNA and RNA oligonucleotides

Results from optical melting studies of Watson–Crick complementary heteroduplexes formed between 2′-O-methyl RNA and RNA oligonucleotides are used to determine nearest neighbor thermodynamic parameters for predicting the stabilities of such duplexes. The results are consistent with the physical model assumed by the individual nearest neighbor-hydrogen bonding model, which contains terms for helix initiation, base pair stacking and base pair composition. The sequence dependence is similar to that for Watson–Crick complementary RNA/RNA duplexes, which suggests that the sequence dependence may also be similar to that for other backbones that favor A-form RNA conformations.     

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Abstract

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.     

Database of non-canonical base pairs found in known RNA structures

Atomic resolution RNA structures are being published at an increasing rate. It is common to find a modest number of non-canonical base pairs in these structures in addition to the usual Watson-Crick pairs. This database summarizes the occurrence of these rare base pairs in accordance with standard nomenclature. The database, http://prion.bchs.uh.edu/, contains information such as sequence context, sugar pucker conformation, anti / syn base conformations, chemical shift, p K (a)values, melting temperature and free energy. Of the 29 anticipated pairs with two or more hydrogen bonds, 20 have been encountered to date. In addition, four unexpected pairs with two hydrogen bonds have been reported bringing the total to 24. Single hydrogen bond versions of five of the expected geometries have been encountered among the single hydrogen bond interactions. In addition, 18 different types of base triplets have been encountered, each of which involves three to six hydrogen bonds. The vast majority of the rare base pairs are antiparallel with the bases in the anti configuration relative to the ribose. The most common are the GU wobble, the Sheared GA pair, the Reverse Hoogsteen pair and the GA imino pair.     

Energy hyperspace for stacking interaction in AU/AU dinucleotide step: Dispersion‐corrected density functional theory study

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence‐directed features, such as small roll values for the purine–pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near‐zero slide to be most favorable for the purine–pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar–phosphate backbone with C3′‐endo sugars and this demands C1′–C1′ distance of about 5.4 Å along the chains. Consideration of an energy penalty term for deviation of C1′–C1′ distance from the mean value, to the recent DFT‐D functionals, specifically ωB97X‐D appears to predict reliable energy contour for AU/AU step. Such distance‐based penalty improves energy contours for the other purine–pyrimidine sequences also. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 107–120, 2014.     

Theoretical studies on protein-nucleic acid interactions. III. Stacking of aromatic amino acids with bases and base pairs of nucleic acids

Stacking of aromatic amino acids tryptophan (Trp), tyrosine (Tyr), phenylalanine (Phe), and histidine (His) with bases and base pairs of nucleic acids has been studied. Stacking energies of the amino acid–base (or base pair) complexes have been calculated by second-order perturbation theory. Our results show that, in general, the predominant contribution to the total stacking energy comes from the dispersion terms. In these cases, repulsion energy is greater than the sum of electrostatic and polarization energies. In contrast to this, interaction of histidine with the bases and base pairs is largely Coulombic in nature. The complexes of guanine with aromatic amino acids are more stable than the corresponding complexes of adenine. Among pyrimidines, cytosine forms the most stable complexes with the aromatic amino acids. The G · C base pair has the highest affinity with aromatic amino acids among various sets of base pairs. Optimized geometries of the stacked complexes show that the aromatic moieties overlap only partially. The heteroatom of one residue generally overlaps with the other aromatic moiety. There is a considerable degree of configurational freedom in the stacked geometries. The role of stacking in specific recognition of base sequences by proteins is discussed.     

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure

The crystal structure of d(CATGGGCCCATG)₂ shows unique stacking patterns of a stable B↔A-DNA intermediate. We evaluated intrinsic base stacking energies in this crystal structure using an ab initio quantum mechanical method. We found that all crystal base pair steps have stacking energies close to their values in the standard and crystal B-DNA geometries. Thus, naturally occurring stacking geometries were essentially isoenergetic while individual base pair steps differed substantially in the balance of intra-strand and inter-strand stacking terms. Also, relative dispersion, electrostatic and polarization contributions to the stability of different base pair steps were very sensitive to base composition and sequence context. A large stacking flexibility is most apparent for the CpA step, while the GpG step is characterized by weak intra-strand stacking. Hydration effects were estimated using the Langevin dipoles solvation model. These calculations showed that an aqueous environment efficiently compensates for electrostatic stacking contributions. Finally, we have carried out explicit solvent molecular dynamics simulation of the d(CATGGGCCCATG)₂ duplex in water. Here the DNA conformation did not retain the initial crystal geometry, but moved from the B↔A intermediate towards the B-DNA structure. The base stacking energy improved in the course of this simulation. Our findings indicate that intrinsic base stacking interactions are not sufficient to stabilize the local conformational variations in crystals.     

X-Ray analysis of d(CGCGAATTXGCG)(2) containing a 2' or minute-deoxy-N(4)-methoxycytosine residue at X: a characteristic pattern of sugar puckers in the crystalline state of the Dickerson-Drew type DNA dodecamers

In a series of structural studies on damaged DNA, a modified Dickerson-Drew dodecamer with the sequence d(CGCGAATTmo(4)CGCG), where mo(4)C is 2'-deoxy-N(4)-methoxycytidine, was synthesized and its structure in a new crystal form has been determined by the X-ray diffraction method. The two dodecamers form a B-form duplex, in which the two mo(4)C residues, respectively, form a wobble pair and a Watson-Crick type pair with the guanine residues of the opposite strand. A comparison of the sugar conformations with those of the other related Dickerson-Drew dodecamers indicates a common feature of their puckering patterns. The sugar pucker of the third residue always adopts an intermediate state (C4'-exo-O4'-endo) between the A-form and B-form. This deviation is ascribed to the stacking interaction of the ribose ring at the third residue with the guanine base at the 12th residue, which is brought about by an extra G12:G2 interaction between two duplexes related by a crystallographic 2(1) symmetry.     

Stacking geometry for non‐canonical G:U wobble base pair containing dinucleotide sequences in RNA: dispersion‐corrected DFT‐D study

Emergence of thousands of crystal structures of noncoding RNA molecules indicates its structural and functional diversity. RNA function is based upon a large variety of structural elements which are specifically assembled in the folded molecules. Along with the canonical Watson‐Crick base pairs, different orientations of the bases to form hydrogen‐bonded non‐canonical base pairs have also been observed in the available RNA structures. Frequencies of occurrences of different non‐canonical base pairs in RNA indicate their important role to maintain overall structure and functions of RNA. There are several reports on geometry and energetic stabilities of these non‐canonical base pairs. However, their stacking geometry and stacking stability with the neighboring base pairs are not well studied. Among the different non‐canonical base pairs, the G:U wobble base pair (G:U W:WC) is most frequently observed in the RNA double helices. Using quantum chemical method and available experimental data set we have studied the stacking geometry of G:U W:WC base pair containing dinucleotide sequences in roll‐slide parameters hyperspace for different values of twist. This study indicates that the G:U W:WC base pair can stack well with the canonical base pairs giving rise to large interaction energy. The overall preferred stacking geometry in terms of roll, twist and slide for the eleven possible dinucleotide sequences is seen to be quite dependent on their sequences. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 328–338, 2015.     

RNABC: forward kinematics to reduce all-atom steric clashes in RNA backbone

Although accurate details in RNA structure are of great importance for understanding RNA function, the backbone conformation is difficult to determine, and most existing RNA structures show serious steric clashes (>or= 0.4 A overlap) when hydrogen atoms are taken into account. We have developed a program called RNABC (RNA Backbone Correction) that performs local perturbations to search for alternative conformations that avoid those steric clashes or other local geometry problems. Its input is an all-atom coordinate file for an RNA crystal structure (usually from the MolProbity web service), with problem areas specified. RNABC rebuilds a suite (the unit from sugar to sugar) by anchoring the phosphorus and base positions, which are clearest in crystallographic electron density, and reconstructing the other atoms using forward kinematics. Geometric parameters are constrained within user-specified tolerance of canonical or original values, and torsion angles are constrained to ranges defined through empirical database analyses. Several optimizations reduce the time required to search the many possible conformations. The output results are clustered and presented to the user, who can choose whether to accept one of the alternative conformations. Two test evaluations show the effectiveness of RNABC, first on the S-motifs from 42 RNA structures, and second on the worst problem suites (clusters of bad clashes, or serious sugar pucker outliers) in 25 unrelated RNA structures. Among the 101 S-motifs, 88 had diagnosed problems, and RNABC produced clash-free conformations with acceptable geometry for 71 of those (about 80%). For the 154 worst problem suites, RNABC proposed alternative conformations for 72. All but 8 of those were judged acceptable after examining electron density (where available) and local conformation. Thus, even for these worst cases, nearly half the time RNABC suggested corrections suitable to initiate further crystallographic refinement. The program is available from http://kinemage.biochem.duke.edu .     

Structural Dynamics of Human Telomeric G-Quadruplex Loops Studied by Molecular Dynamics Simulations

Loops which are linkers connecting G-strands and supporting the G-tetrad core in G-quadruplex are important for biological roles of G-quadruplexes. TTA loop is a common sequence which mainly resides in human telomeric DNA (hTel) G-quadruplex. A series of molecular dynamics (MD) simulations were carried out to investigate the structural dynamics of TTA loops. We found that (1) the TA base pair formed in TTA loops are very stable, the occupied of all hydrogen bonds are more than 0.95. (2) The TA base pair makes the adjacent G-quartet more stable than others. (3) For the edgewise loop and the diagonal loop, most loop bases are stacking with others, only few bases have considerable freedom. (4) The stabilities of these stacking structures are distinct. Part of the loops, especially TA base pairs, and bases stacking with the G-quartet, maintain certain stable conformations in the simulation, but other parts, like TT and TA stacking structures, are not stable enough. For the first time, spontaneous conformational switches of TTA edgewise loops were observed in our long time MD simulations. (5) For double chain reversal loop, it is really hard to maintain a stable conformation in the long time simulation under present force fields (parm99 and parmbsc0), as it has multiple conformations with similar free energies.