Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs

Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non‐canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte‐Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER‐parm99bsc0 and CHARMM‐36 force‐fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte‐Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson–Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non‐canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse‐grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson–Crick base pairs, as found often in the non‐coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212–226, 2016.     

Stacking interactions in RNA and DNA: Roll‐slide energy hyperspace for ten unique dinucleotide steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT‐D by ωB97X‐D/6‐31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B‐DNA, whereas, smaller twist values correspond to higher stability to RNA and A‐DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B‐DNA or A‐DNA conformations. Conformational preference of BII sub‐state in B‐DNA is preferentially displayed mainly by pyrimidine–purine steps and partly by purine–purine steps. The purine–pyrimidine steps show largest effect of 5‐methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 134–147, 2015.     

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.     

A comparison of four different conformations adopted by human telomeric G‐quadruplex using computer simulations

The telomeric G‐quadruplexes for their unique structural features are considered as potential anticancer drug targets. These, however, exhibit structural polymorphism as different topology types for the intra‐molecular G‐quadruplexes from human telomeric G‐rich sequences have been reported based on NMR spectroscopy and X‐ray crystallography. These techniques provide detailed atomic‐level information about the molecule but relative conformational stability of the different topologies remains unsolved. Therefore, to understand the conformational preference, we have carried out quantum chemical calculations on G‐quartets; used all‐atom molecular dynamics (MD) simulations and steered molecular dynamics (SMD) simulations to characterize the four human telomeric G‐quadruplex topologies based on its G‐tetrad core‐types, viz., parallel, anti‐parallel, mixed‐(3 + 1)‐form1 and mixed‐(3 + 1)‐form2. We have also studied a non‐telomeric sequence along with these telomeric forms giving a comparison between the two G‐rich forms. The structural properties such as base pairing, stacking geometry and backbone conformations have been analyzed. The quantum calculations indicate that presence of a sodium ion inside the G‐tetrad plane or two potassium ions on both sides of the plane give it an overall planarity which is much needed for good stacking to form a helix. MD simulations indicate that capping of the G‐tetrad core by the TTA loops keep the terminal guanine bases away from water. The SMD simulations along with equilibrium MD studies indicate that the parallel and non‐telomeric forms are comparatively less stable. We could come to the conclusion that the anti‐parallel form and also the mixed‐(3 + 1)‐form1 topology are most likely to represent the major conformation., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 83–99, 2016     

Diverse mutants of HIV RRE IIB recognize wild‐type Rev ARM or Rev ARM R35G‐N40V

The binding of human immunodeficiency virus Rev protein via its arginine‐rich motif (ARM) to an internal loop in the Rev‐response element region IIB (RRE IIB) is necessary for viral replication. Many variant RNAs and ARMs that bind Rev and RRE IIB have been found. Despite the essential role of Rev asparagine 40 in recognition, the Rev ARM double‐mutant R35G‐N40V functions well in a Rev–RRE IIB reporter assay, indicating R35G‐N40V uses a distinct recognition strategy. To examine how RRE IIB may evolve specificity to wild‐type Rev ARM and R35G‐N40V, 10 RRE IIB libraries, each completely randomized in overlapping regions, were screened with wild‐type Rev ARM and R35G‐N40V using a reporter system based on bacteriophage λ N antitermination. Consistent with previous studies, a core element of RRE IIB did not vary, and substitutions occurred at conserved residues only in the presence of other substitutions. Notably, the groove‐widening, non‐canonical base‐pair G48:G71 was mutable to U48:G71 without strong loss of binding to wild‐type Rev ARM, suggesting U48:G71 performs the same role by adopting the nearly isosteric, reverse wobble base pair. Originating from RRE IIB, as few as one or two substitutions are sufficient to confer specificity to wild‐type Rev or Rev R35G‐N40. The diversity of RRE IIB mutants that maintain binding to wild‐type Rev ARM and R35G‐N40V supports neutral theories of evolution and illustrates paths by which viral RNA–protein interactions can evolve new specificities. Rev–RRE offers an excellent model with which to study the fine structure of how specificity evolves. Copyright © 2015 John Wiley & Sons, Ltd.     

Probing the role of sequence in the assembly of three‐dimensional DNA crystals

DNA is a widely used biopolymer for the construction of nanometer‐scale objects due to its programmability and structural predictability. One long‐standing goal of the DNA nanotechnology field has been the construction of three‐dimensional DNA crystals. We previously determined the X‐ray crystal structure of a DNA 13‐mer that forms a continuously hydrogen bonded three‐dimensional lattice through Watson‐Crick and non‐canonical base pairs. Our current study sets out to understand how the sequence of the Watson‐Crick duplex region influences crystallization of this 13‐mer. We screened all possible self‐complementary sequences in the hexameric duplex region and found 21 oligonucleotides that crystallized. Sequence analysis showed that one specific Watson‐Crick pair influenced the crystallization propensity and the speed of crystal self‐assembly. We determined X‐ray crystal structures for 13 of these oligonucleotides and found sequence‐specific structural changes that suggests that this base pair may serve as a structural anchor during crystal assembly. Finally, we explored the crystal self‐assembly and nucleation process. Solution studies indicated that these oligonucleotides do not form base pairs in the absence of cations, but that the addition of divalent cations leads to rapid self‐assembly to higher molecular weight complexes. We further demonstrate that crystals grown from mixtures of two different oligonucleotide sequences contain both oligonucleotides. These results suggest that crystal self‐assembly is nucleated by the formation of the Watson‐Crick duplexes initiated by a simple chemical trigger. This study provides new insight into the role of sequence for the assembly of periodic DNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 618–626, 2015.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

DNA bending: the prevalence of kinkiness and the virtues of normality.

DNA bending in 86 complexes with sequence-specific proteins has been examined using normal vector plots, matrices of normal vector angles between all base pairs in the helix, and one-digit roll/slide/twist tables. FREEHELIX, a new program especially designed to analyze severely bent and kinked duplexes, generates the foregoing quantities plus local roll, tilt, twist, slide, shift and rise parameters that are completely free of any assumptions about an overall helix axis. In nearly every case, bending results from positive roll at pyrimidine-purine base pair steps: C-A (= T-G), T-A, or less frequently C-G, in a direction that compresses the major groove. Normal vector plots reveal three well-defined types of bending among the 86 examples: (i) localized kinks produced by positive roll at one or two discrete base pairs steps, (ii) three-dimensional writhe resulting from positive roll at a series of adjacent base pairs steps, or (iii) continuous curvature produced by alternations of positive and negative roll every 5 bp, with side-to-side zig-zag roll at intermediate position. In no case is tilt a significant component of the bending process. In sequences with two localized kinks, such as CAP and IHF, the dihedral angle formed by the three helix segments is a linear function of the number of base pair steps between kinks: dihedral angle = 36 degrees x kink separation. Twenty-eight of the 86 examples can be described as major bends, and significant elements in the recognition of a given base sequence by protein. But even the minor bends play a role in fine-tuning protein/DNA interactions. Sequence-dependent helix deformability is an important component of protein/DNA recognition, alongside the more generally recognized patterns of hydrogen bonding. The combination of FREEHELIX, normal vector plots, full vector angle matrices, and one-digit roll/slide/twist tables affords a rapid and convenient method for assessing bending in DNA.     

Control of box C/D snoRNP assembly by N6‐methylation of adenine

N⁶‐methyladenine is the most widespread mRNA modification. A subset of human box C/D snoRNA species have target GAC sequences that lead to formation of N⁶‐methyladenine at a key trans Hoogsteen‐sugar A·G base pair, of which half are methylated in vivo. The GAC target is conserved only in those that are methylated. Methylation prevents binding of the 15.5‐kDa protein and the induced folding of the RNA. Thus, the assembly of the box C/D snoRNP could in principle be regulated by RNA methylation at its critical first stage. Crystallography reveals that N⁶‐methylation of adenine prevents the formation of trans Hoogsteen‐sugar A·G base pairs, explaining why the box C/D RNA cannot adopt its kinked conformation. More generally, our data indicate that sheared A·G base pairs (but not Watson–Crick base pairs) are more susceptible to disruption by N⁶mA methylation and are therefore possible regulatory sites. The human signal recognition particle RNA and many related Alu retrotransposon RNA species are also methylated at N6 of an adenine that forms a sheared base pair with guanine and mediates a key tertiary interaction.     

In vitro selection of DNAs with an increased propensity to form small circles

A protocol was devised to select for DNA molecules that efficiently form circles from a library of 126 base pair DNAs containing 90 randomized base pairs. After six rounds of selection, individual molecules from the library showed 20‐ to 100‐fold greater j‐factors compared with the starting library, validating the selection protocol. High‐throughput sequencing revealed a sinusoidal pattern of enrichment and de‐enrichment of A/T dinucleotides in the random region with a 10.4 base pair period associated with the helicity of DNA. A similar, but more moderate pattern of C/G dinucleotides was offset by precisely half a helical turn. While C/G dinucleotide enrichments were evenly distributed, A/T dinucleotide enrichments displayed a preference to cluster in individual DNA molecules. The most highly enriched 10 base pair sequences in the random region contained adjacent blocks of A/T and C/G trinucleotides present in some, but not all, rapidly cyclizing molecules. The phased dinucleotide enrichments closely match those present in accurately mapped yeast nucleosomes, confirming the importance of DNA bending in nucleosome formation. However, at certain sites the nucleosomal DNAs show dinucleotide enrichments that differ substantially from the cyclization data. These discrepancies can often be correlated with sequence specific contacts that form between histones and DNA. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 303–320, 2015.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction

BackgroundMolecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs.MethodsWe employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair.ResultsWe noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear.ConclusionsThe quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs “tug of war” situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity.General significanceWe have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.     

Database proton NMR chemical shifts for RNA signal assignment and validation

The Biological Magnetic Resonance Data Bank contains NMR chemical shift depositions for 132 RNAs and RNA-containing complexes. We have analyzed the ¹H NMR chemical shifts reported for non-exchangeable protons of residues that reside within A-form helical regions of these RNAs. The analysis focused on the central base pair within a stretch of three adjacent base pairs (BP triplets), and included both Watson–Crick (WC; G:C, A:U) and G:U wobble pairs. Chemical shift values were included for all 4³ possible WC-BP triplets, as well as 137 additional triplets that contain one or more G:U wobbles. Sequence-dependent chemical shift correlations were identified, including correlations involving terminating base pairs within the triplets and canonical and non-canonical structures adjacent to the BP triplets (i.e. bulges, loops, WC and non-WC BPs), despite the fact that the NMR data were obtained under different conditions of pH, buffer, ionic strength, and temperature. A computer program (RNAShifts) was developed that enables convenient comparison of RNA ¹H NMR assignments with database predictions, which should facilitate future signal assignment/validation efforts and enable rapid identification of non-canonical RNA structures and RNA-ligand/protein interaction sites.     

Structural basis for a novel mechanism of DNA bridging and alignment in eukaryotic DSB DNA repair

Eukaryotic DNA polymerase mu of the PolX family can promote the association of the two 3′‐protruding ends of a DNA double‐strand break (DSB) being repaired (DNA synapsis) even in the absence of the core non‐homologous end‐joining (NHEJ) machinery. Here, we show that terminal deoxynucleotidyltransferase (TdT), a closely related PolX involved in V(D)J recombination, has the same property. We solved its crystal structure with an annealed DNA synapsis containing one micro‐homology (MH) base pair and one nascent base pair. This structure reveals how the N‐terminal domain and Loop 1 of Tdt cooperate for bridging the two DNA ends, providing a templating base in trans and limiting the MH search region to only two base pairs. A network of ordered water molecules is proposed to assist the incorporation of any nucleotide independently of the in trans templating base. These data are consistent with a recent model that explains the statistics of sequences synthesized in vivo by Tdt based solely on this dinucleotide step. Site‐directed mutagenesis and functional tests suggest that this structural model is also valid for Pol mu during NHEJ.     

An innate twist between Crick’s wobble and Watson-Crick base pairs

Non-Watson-Crick pairs like the G·U wobble are frequent in RNA duplexes. Their geometric dissimilarity (nonisostericity) with the Watson-Crick base pairs and among themselves imparts structural variations decisive for biological functions. Through a novel circular representation of base pairs, a simple and general metric scheme for quantification of base-pair nonisostericity, in terms of residual twist and radial difference that can also envisage its mechanistic effect, is proposed. The scheme is exemplified by G·U and U·G wobble pairs, and their predicable local effects on helical twist angle are validated by MD simulations. New insights into a possible rationale for contextual occurrence of G·U and other non-WC pairs, as well as the influence of a G·U pair on other non-Watson-Crick pair neighborhood and RNA-protein interactions are obtained from analysis of crystal structure data. A few instances of RNA-protein interactions along the major groove are documented in addition to the well-recognized interaction of the G·U pair along the minor groove. The nonisostericity-mediated influence of wobble pairs for facilitating helical packing through long-range interactions in ribosomal RNAs is also reviewed.     

Base interaction and conformation manifestations of repetitive nucleotide sequences

To understand why different nucleotide sequences prefer different double helical conformations and to predict conformational behaviour of definite sequences the base-base interaction energy in regular helices consisting of A:U, A:T, G:C and I:C (hypoxanthine-cytosine) base pairs was calculated. Interaction energy was assumed to be a function of eight conformational parameters: H, the distance between adjacent pairs along helix axes; tau, turn angle of one pair relative to the neighbouring one; angles between base planes in a pair (TW, propeller twist and BL, buckle) and position of pairs with respect to helix axes (D and SL, displacements in the plane normal to helix axes, and TL and RL, inclinations to this plane, tilt and roll, respectively). For H and tau characteristic of A- and B-families of nucleic acid conformations (2.5 A less than H less than or equal to 3.5 A, 30 degrees less than or equal to tau less than or equal to 45 degrees) the ranges of conformational parameters corresponding to energy values close to minimal ones (valleys) and correlations between conformational parameters were revealed. Valleys for different sequences largely coincide but have distinctive characteristics for each sequence. Reasons for base pair planarity distortion in double-stranded helices were considered. The calculations permit to account for A-phility of G:C sequences and B-phility for A:T sequences. The valley for I:C sequence branches. This corresponds to A:T-like behaviour in some cases and G:C-like in the others.     

RNA approaches the B‐form in stacked single strand dinucleotide contexts

Duplex RNA adopts an A‐form structure, while duplex DNA interconverts between the A‐ and B‐forms depending on the environment. The C2′‐endo sugar pucker seen in B‐form DNA can occur infrequently in ribose sugars as well, but RNA is not understood to assume B‐form conformations. Through analysis of over 45,000 stacked single strand dinucleotide (SSD) crystal structure conformations, this study demonstrates that RNA is capable of adopting a wide conformational range between the canonical A‐ and B‐forms at the localized SSD level, including many B‐form‐like conformations. It does so through C2′‐endo ribose conformations in one or both nucleotides, and B‐form‐like neighboring base stacking patterns. As chemical reactions on nucleic acids involve localized changes in chemical bonds, the understanding of how enzymes distinguish between DNA and RNA nucleotides is altered by the energetic accessibility of these rare B‐form‐like RNA SSD conformations. The existence of these conformations also has direct implications in parametrization of molecular mechanics energy functions used extensively to model nucleic acid behavior., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 65–82, 2016     

Opening dynamics of 8‐oxoguanine in DNA

8‐oxoguanine is a major lesion of genomic DNA that results from oxidation of guanine by reactive oxygen species. The repair of this lesion is initiated by 8‐oxoguanine glycosylases, which excise the damaged base by “flipping” it outside the DNA double helix. The molecular mechanisms involved in the specific recognition of the damaged base by the enzyme are not yet fully understood. Several models have proposed that, in DNA, the base pair between 8‐oxoguanine and cytosine may possess altered dynamic properties that could help the enzyme locate the lesion and could favor the selective extra‐helical flipping of the damaged base. To test this proposal, we have characterized the spontaneous opening of the base pair between 8‐oxoguanine and cytosine in a DNA double helix using NMR spectroscopy and proton exchange. The results show that the rate of spontaneous opening of 8‐oxoguanine and the lifetime of the base in the extra‐helical state are the same as those of a canonical guanine‐cytosine base pair, in the same base sequence context. This finding suggests that the opening dynamics of 8‐oxoguanine, when paired with cytosine in DNA, does not play a significant role in the recognition of the lesion by glycosylases. Copyright © 2013 John Wiley & Sons, Ltd.     

Is there a special function for U.G basepairs in ribosomal RNA?

U.G basepairs are well-established elements of RNA structure. The geometry of this pair is different, however, from classical Watson-Crick basepairs. This leads to an unusual stacking of the basepair: overlap with the basepair at the 5' side of the U (and the 3' side of the G) is strong (stacked) while it is weak with the basepair on the other side (destacked). The closure of an RNA helix by a U.G pair will be energetically unfavourable when the U residue occupies the 5' end. In transfer RNA there is a strong selection against a 'destacked' U.G pair at helix ends. In the 16S rRNA model of Escherichia coli there are 72 U.G pairs of which 36 or 22 occupy a helix end, depending on how such an end is defined. There is a slight preference for 'stacked' U.G's in these positions. It is remarkable, however, that of 13 very conserved U.G pairs in the 16S (-like) rRNA, 7 occur at helix ends and that 5 of these have the 'destacked' configuration. It is suggested that these pairs, if they exist at all in a hydrogen-bounded form, are stabilized by co-axial stacking with other helices or by interaction with a protein.     

Sequence effects on the propeller twist of base pairs in DNA helices.

Molecular mechanics calculations are performed on all the ten base pair steps (duplex dimers) and also a number of trimer and tetramer duplexes comprising them, in an attempt to systematically examine the possible base sequence effects on the magnitudes of propeller twists of base pairs at a given step. The analysis reveals that though propeller twist is a base pair property, it behaves very much like other base step parameters such as slide, roll, twist etc., Hence, it may be necessary to monitor the nature and variation of magnitudes of pt at a step. Calculations performed on 45 out of the 136 unique tetramer combinations involving all the ten unique base steps show that the difference in magnitudes of propeller twists of the base pairs of a given step has been found to be either steep or moderate depending on base pairs that flank the base step. These observations compare very well with the available experimental data. Tetramer sequences, wherein a base pair of a base step repeats in the same direction, exhibit a relatively steep difference in propeller twist at the step. Tetramers other than these exhibit moderate difference in propeller twist. Such sequences are broadly classified as type-I and type-II respectively. Practically all the tetrads considered in the study, excepting those with GT step and a few involving CG and GC steps, conform to the above classification.