Variability in the Geometry of RNA Double Helices Generated from Dinucleoside Building-Blocks

Abstract

A graphical method is presented for the generation of helical parameters from single-crystal structures of RNA nucleic acid fragments that are minimally dinucleosides. The method is compared with other published procedures, for a number of text examples. The RNA double helices generated from three different salts of the dinucleoside monophosphate GpC are examined in relation to the variations in helix morphology that are produced. It is shown that small differences between these GpC salts can be amplified to very distinct helix characteristics.     

Conformational stability in dinucleoside phosphate crystals. Semiempirical potential energy calculations for uridylyl-3'-5'-adenosine monophosphate (UpA) and guanylyl-3',5'-cytidine monophosphate (GpC)

Classical potential energy calculations were performed for the dinucleoside phosphates UpA and GpC. Two widely accessible low-energy regions of conformation space were found for the ω′, ω pair. That of lowest energy contains conformations similar to helical RNA, with ω′ and ω in the vicinity of 300° and 280°, respectively. All five experimental observations of crystalline GpC, two of ApU, and the helical fragment of ApApA fall in this range. The second lowest region has ω′ and ω at about 20° and 80°, respectively, which is in the general region of one experimentally observed crystalline conformer of UpA, and the nonhelical region of ApApA. It is concluded that GpC and ApU, which were crystallized as either sodium or calcium salts, are shielded from each other in the crystal by the water of hydration and are therefore free to adopt their predicted in vacuo minimum energy helical conformations. By contrast, crystalline UpA had only 1/2 water per molecule, and was forced into higher energy conformations in order to maximize intermolecular hydrogen bonding.     

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.     

RNA double helices generated from crystal structures of double helical dinucleoside phosphates.

The dinucleoside phosphates ApU and GpC form right-handed anti-parallel double helical fragments within their crystal lattices. Using a least squares procedure, we have generated the extended double helices which these fragments represent. ApU corresponds to a double helix with 11.9 residues per turn and a pitch of 28. 1Å. The GpC double helix has 10.4 residues per turn and a pitch of 26. 9Å.     

Subsite interactions of ribonuclease T1: viscosity effects indicate that the rate-limiting step of GpN transesterification depends on the nature of N

We report on the effect of the viscogenic agents glycerol and ficoll on the RNase T1 catalyzed turnover of GpA, GpC, GpU, and Torula yeast RNA. For wild-type enzyme, we find that the kcat/Km values for the transesterification of GpC and GpA as well as for the cleavage of RNA are inversely proportional to the relative viscosity of glycerol-containing buffers; no such effect is observed for the conversion of GpU to cGMP and U. The second-order rate constants for His40Ala and Glu46Ala RNase T1, two mutants with a drastically reduced kcat/km ratio, are independent of the microviscosity, indicating that glycerol does not affect the intrinsic kinetic parameters. Consistent with the notion that molecular diffusion rates are unaffected by polymeric viscogens, addition of ficoll has no effect on the kcat/Km for GpC transesterification by wild-type enzyme. The data indicate that the second-order rate constants for GpC, GpA, and Torula yeast RNA are at least partly limited by the diffusion-controlled association rate of substrate and active site; RNase T1 obeys Briggs-Haldane kinetics for these substrates (Km greater than Ks). Calculations suggest that the equilibrium dissociation constants (Ks) for the various GpN-wild-type enzyme complexes are virtually independent of N whereas the measured kcat values follow the order GpC greater than GpA greater than GpU. This is also revealed by the steady-state kinetic parameters of Tyr38Phe and His40Ala RNase T1, two mutants that follow simple Michaelis-Menten kinetics because of a dramatically reduced kcat value (i.e., Km = Ks).(ABSTRACT TRUNCATED AT 250 WORDS)     

Hydrogen-bonded complexes of the ribodinucleoside monophosphates in aqueous solution. Proton magnetic resonance studies.

A proton magnetic resonance study of the chemical shifts of a series of ribodinucleoside monophosphates in neutral H2O solution has been recorded in the 1-100 mM concentration range. The self-complementary dinucleoside monophosphates CpG and GpC and the complementary mixture GpU + ApC form intermolecular hydrogen-bonded complexes at low temperatures. The amino proton chemical shifts in the CpG and GpC spectra are consistent with the formation of a miniature double helical dimer in neutral aqueous solution at low temperatures (approximately 2 degrees C). The complementary mixture of dinucleosides GpU + ApC formed much less stable complexes than either GpC or CpG, while UpA did not show any indication of the formation of intermolecular hydrogen-bonded complexes. This result is consistent with the well-known observation that the stability of a double helix is proportional to the percent of G-C base pairs present.     

Influence of ribose 2′-O-methylation on GpC conformation by classical potential energy calculations

Potential energy calculations were employed to examine the effect of ribose 2′-O-methylation on the conformation of GpC. Minimum energy conformations and allowed conformational regions were calculated for 2′MeGpC and Gp2′MeC. The two lowest energy conformations of 2′MeGpC and Gp2′MeC are similar to those of GpC itself. The helical RNA conformation (sugar pucker-C(3′)-endo, ω′ and ω,g^?g^?, bases-anti) is the global minimum, and a helix-reversing conformation with ω′, ω in the vicinity of 20°, 80° is next in energy. However, subtle differences between the three molecules are noted. When the substitution is on the 5′ ribose (Gp2′MeC), the energy of the helical conformation is less than that of GpC, due to favorable interactions of the added methyl group. When the substitution is at the 3′ ribose (2′MeGpC) these stabilizing interactions are outweighed by steric restrictions, and the helical conformation is of higher energy than for GpC. Furthermore, the statistical weight of the 2′MeGpC g^? g^? helical region is substantially less than the corresponding weight for Gp2′MeC. In addition, 2′MeGpC′s methoxy group is conformationally restricted to a narrow range centered at 76°. This group has a broadly allowed region between 50 and 175° in Gp2′MeC. These differences occur because the appended methyl group in 2′MeGpC is located in the interior of the helix cylinder, as it would be in polynucleotide, while it hangs unimpeded in Gp2′MeC. These findings suggest that 2′-O-methylation has both stabilizing and destabilizing influences on the helical conformation of RNA. For 2′MeGpC the destabilizing steric hindrance imposed by the nature of the guanine base dominates.     

The A and B conformations of DNA and RNA subunits. Potential energy calculations for dGpdC

In order to obtain a molecular picture of the A and B forms of a DNA subunit, potential energy calculations have been made for dGpdC with C(3′)-endo and C(2′)-endo [or C(3′)-exo] sugar puckerings. These are compared with results for GpC. The global minima for dGpdC and GpC are almost identical. They are like A-form duplex DNA and RNA, respectively, with bases anti, the ω′, ω angle pair near 300°, 280°, and sugar pucker C(3′)-endo. For dGpdC, a B-form helical conformer, with sugar pucker C(2′)-endo and ω′ = 257°, ω = 298°, is found only 0.4 kcal/mol above the global minimum. A second low-energy conformation (2.3 kcal/mol) has ω′ = 263°, ω = 158° and ψ near 180°. This has dihedral angles like the original Watson–Crick model of the double helix. In contrast, for GpC, the C(2′)-endo B form is 6.9 kcal/mol above the global minimum. These theoretical results are consistent with experimental studies on DNA and RNA fibers. DNA fibers exist in both A and B forms, while RNA fibers generally assume only the A form. A low-energy conformation unlike the A or B forms was found for both dGpdC and GpC when the sugars were C(3′)-endo. This conformation—ω′,ω near 20°,80°—was not observed for C(2′)-endo dGpdC. Energy surface maps in the ω′,ω plane showed that C(2′)-endo dGpdC has one low-energy valley. It is in the B-form helical region (ω′ ～ 260°, ω ～ 300). When the sugar pucker is C(3′)-endo, dGpdC has two low-energy regions: the A-form helical region and the region with the minimum at ω′ = 16°, ω = 85°.     

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

Effects of neutral salts on the conformational state of glycogen, interacting with iodine

The effects of salts on helix formation of linear (outer) chains of branched alpha-glucane, muscle glycogen, was judged upon by the state of polyiodine chromophore of the iodine reaction surrounded by a polysaccharide asymmetric helix and recorded by the circular dichroism method. It was demonstrated that apart from the known changes in the absorption of the iodine-glycogen complex the salts induce changes in the ellipticity by affecting the helix formation of the linear chains of the polysaccharide. The nature of these effects depends on the type and concentration of the salt used. Monovalent metal salts produce a more favourable effect on helix formation than the divalent metal salts. Among divalent metal salts transient metal salts produce a weaker effect. The action of neutral salts on helix formation is of dual nature: helix formation and helix destruction. The helix formation is largely impaired by high concentrations of LiCl (greater than 6 M) or CaCl2 (greater than 2M). It was shown that a certain role in the mechanism of neutral salt action on helix formation belongs both to electrostatic and lyotropic effects. The latter determine the effects of salts on helix formation at concentrations above the monomolar one, when the specific effect of a salt is especially well-pronounced. The salts which enhance the orderliness of H2O structure produce a favourable effect on helix formation of the linear polysaccharide chains and, consequently, on the iodine reaction. In terms of molar efficiency of their action on the iodine reaction the salts correspond to the Hofmeister lyotropic sequence.     

Thermodynamic and kinetic properties of short RNA helices: the oligomer sequence AnGCUn

We studied the thermodynamic, kinetic and optical properties of the double helices formed by the series of self-complementary oligonucleotides, (AP)_nGpC(pU)_n, 2 ≤ n ≤ 4, and found that the shortest helix, containing just 6 base pairs, is less stable than would be predicted from the properties of the larger molecules. It also shows a markedly smaller hyperchromism on melting than expected. These anomalous properties of a short helix indicate that one cannot always assume that base pair free energies and extinction coefficient changes are independent of helix size.     

Flexibility and rigidity of left-handed Z-DNA

The local variation of torsional angles and helical parameters in Z-DNA was analyzed. The sugar phosphate backbone is fairly rigid but the angles at GpC linkage are more changeable than those at CpG linkage in order to form a variety of structures. The water channel at minor groove is important to stabilize and retain the novel Z-DNA helix.     

Two different molecular conformations found in chitosan type II salts

The type II structure of chitosan acidic salts prepared from crab tendon in solid state was studied using an X-ray fiber diffraction technique together with the linked-atom least-squares (LALS) technique. The cylindrical Patterson method was applied to confirm the molecular conformation of the chitosan. It was shown that there are two different helical conformations for type II salts. One is the relaxed twofold helix having a tetrasaccharide as an asymmetric unit as found in chitosan.HCl salt, which was previously reported as a conformation of chitosan.HCOOH salt. The other is the fourfold helix having a disaccharide as an asymmetric unit newly found in chitosan.HI salt.     

High salt solution structure of a left-handed RNA double helix

Right-handed RNA duplexes of (CG)_n sequence undergo salt-induced helicity reversal, forming left-handed RNA double helices (Z-RNA). In contrast to the thoroughly studied Z-DNA, no Z-RNA structure of natural origin is known. Here we report the NMR structure of a half-turn, left-handed RNA helix (CGCGCG)₂ determined in 6 M NaClO₄. This is the first nucleic acid motif determined at such high salt. Sequential assignments of non-exchangeable proton resonances of the Z-form were based on the hitherto unreported NOE connectivity path [H6_(n)-H5′/H5″_(n)-H8_(n+1)-H1′_(n+1)-H6_(n+2)] found for left-handed helices. Z-RNA structure shows several conformational features significantly different from Z-DNA. Intra-strand but no inter-strand base stacking was observed for both CpG and GpC steps. Helical twist angles for CpG steps have small positive values (4–7°), whereas GpC steps have large negative values (−61°). In the full-turn model of Z-RNA (12.4 bp per turn), base pairs are much closer to the helix axis than in Z-DNA, thus both the very deep, narrow minor groove with buried cytidine 2′-OH groups, and the major groove are well defined. The 2′-OH group of cytidines plays a crucial role in the Z-RNA structure and its formation; 2′-O-methylation of cytidine, but not of guanosine residues prohibits A to Z helicity reversal.     

Crystal structure of the self-complementary 5''-purine start decamer d(GCGCGCGCGC) in the Z-DNA conformation. I.

Alternating self-complementary oligonucleotides starting with a 5'-pyrimidine usually form left-handed Z-DNA; however, with a 5'-purine start sequence they form the right-handed A-DNA. Here we report the crystal structure of the decamer d(GCGCGCGCGC) with a 5'-purine start in the Z-DNA form. The decamer crystallizes in the hexagonal space group P6(5)22, unit cell dimensions a = b = 18.08 and c = 43.10 A, with one of the following four dinucleotide diphosphates in the asymmetric unit: d(pGpC)/d(GpCp)/d(pCpG)/d(CpGp). The molecular replacement method, starting with d(pGpC) of the isomorphous Z-DNA hexamer d(araC-dG)3 without the 2'-OH group of arabinose, was used in the structure analysis. The method gave the solution only after the sugar-phosphate conformation of the GpC step was manipulated. The refinement converged to a final R value of 18.6% for 340 unique reflections in the resolution range 8.0-1.9 A. A result of the sequence alternation is the alternation in the nucleotide conformation; guanosine is C3'-endo, syn, and cytidine is C2'-endo, anti. The CpG step phosphodiester conformation is the same as ZI or ZII, whereas that of the GpC step phosphodiester is "intermediate" in the sense that zeta (O3'-P bond) is the same as ZII but alpha (P-O5' bond) is the same as ZI. The duplexes generated from the dinucleotide asymmetric unit are stacked one on top of the other in the crystal to form an infinite pseudocontinuous helix. This renders it a quasi-polymerlike structure that has assumed the Z-DNA conformation further strengthened by the long inner Z-forming stretch d(CG)4. An interesting feature of the structure is the presence of water strings in both the major and the minor grooves. In the minor groove the cytosine carbonyl oxygen atoms of the GpC and CpG steps are cross-bridged by water molecules that are not themselves hydrogen bonded but are enclosed by the water rings in the mouth of the minor groove. In the major groove three independent water molecules form a zigzagging continuous water string that runs throughout the duplex.     

Fourth 3D structure of the chitosan molecule: conformation of chitosan in its salts with medical organic acids having a phenyl group

Chitosan salts with two medical organic acids having phenyl groups (salicylic and gentisic acids) exhibited fiber diffraction patterns of a new type of crystal which does not compare with known types I and II. The crystals, called type III salts, showed a fiber repeat of 2.550 nm and a meridional reflection at the 5th layer line. These results coupled with a conformational analysis indicate the chain conformation of chitosan with the salts to be a 5/3 helix, this helix differing from those of type I (an extended two-fold helix) and type II (a relaxed two-fold helix or a 4/1 helix). The fiber patterns of all the type III salts were similar. This observation has also been found with type II salts and is an indication that the acid ions are not arranged in regular positions in the crystals. A comparison of solid-state 13C-NMR spectra of the gentisic acid salt and the aspirin salt, which could not be crystallized, suggests that, in the latter salt, the chitosan molecules also formed a 5/3 helix.     

Solid-phase synthesis of the self-complementary hexamer d(c7GpCpc7GpCpc7GpC) via the O-3''-phosphoramidite of 7-deaza-2''-deoxyguanosine.

The synthesis of the O-3'-phosphoramidite of a suitably protected 7-deaza-2'-deoxyguanosine (c7G) which is an isostere of 2'-deoxyguanosine is described. The phosphoramidite of the modified nucleoside was used in the synthesis of the self-complementary hexamer d(c7GpCpc7GpCpc7GpC) on functionalized silica gel in a mini-reactor. As expected from the parent hexamer d(GpCpGpCpGpC) the isosteric d(c7GpCpc7GpCpc7GpC) exhibits a rigid secondary structure (22% hypochromicity at 280 nm) and forms a duplex in 1 M aqueous sodium chloride solution. Due to the altered pi-electron system of the pyrrolo[2,3-d]pyrimidine nucleobase, which affects base stacking and hydrogen bonding, the Tm of the modified duplex is decreased by 10 degrees C compared to that of the parent purine hexamer. Moreover, it is expected that the incorporation of c7G influences the pitch of the helix.     

Probing functional perfection in substructures of ribonuclease T1: double combinatorial random mutagenesis involving Asn43, Asn44, and Glu46 in the guanine binding loop

Combinatorial random mutageneses involving either Asn43 with Asn44 (set 1) or Glu46 with an adjacent insertion (set 2) were undertaken to explore the functional perfection of the guanine recognition loop of ribonuclease T(1) (RNase T(1)). Four hundred unique recombinants were screened in each set for their ability to enhance enzyme catalysis of RNA cleavage. After a thorough selection procedure, only six variants were found that were either as active or more active than wild type which included substitutions of Asn43 by Gly, His, Leu, or Thr, an unplanned Tyr45Ser substitution and Glu46Pro with an adjacent Glu47 insertion. Asn43His-RNase T(1) has the same loop sequence as that for RNases Pb(1) and Fl(2). None of the most active mutants were single substitutions at Asn44 or double substitutions at Asn43 and Asn44. A total of 13 variants were purified, and these were subjected to kinetic analysis using RNA, GpC, and ApC as substrates. Modestly enhanced activities with GpC and RNA involved both k(cat) and K(M) effects. Mutants having low activity with GpC had proportionately even lower relative activity with RNA. Asn43Gly-RNase T(1) and all five of the purified mutants in set 2 exhibited similar values of k(cat)/K(M) for ApC which were the highest observed and about 10-fold that for wild type. The specificity ratio [(k(cat)/K(M))(GpC)/(k(cat)/K(M))(ApC)] varied over 30 000-fold including a 10-fold increase [Asn43His variant; mainly due to a low (k(cat)/K(M))(ApC)] and a 3000-fold decrease (Glu46Ser/(insert)Gly47 variant; mainly due to a low (k(cat)/K(M))(GpC)) as compared with wild type. It is interesting that k(cat) (GpC) for the Tyr45Ser variant was almost 4-fold greater than for wild type and that Pro46/(insert)Glu47 RNase T(1) is 70-fold more active than the permuted variant (insert)Pro47-RNase T(1) which has a conserved Glu46. In any event, the observation that only 6 out of 800 variants surveyed had wild-type activity supports the view that functional perfection of the guanine recognition loop of RNase T(1) has been achieved.