Nearest-neighbor effects in the structure and function of nucleic acids

In vitro codon-reading specificities of unmodified tRNA molecules with different anticodons on the sequence background of Escherichia coli tRNASer

The codon-reading properties of wobble-position variants of the unmodified form of Escherichia coli tRNASer1 (the UGA anticodon) were measured in a cell-free translation system. Two variants, with the AGA and CGA anticodons, each exclusively read a single codon, UCU and UCG, respectively. The only case of efficient wobbling occurred with the variant with the GGA anticodon, which reads the UCU codon in addition to the UCC codon. Surprisingly, this wobble reading is more efficient than the Watson-Crick reading by the variant with the AGA anticodon. Furthermore, we prepared tRNA variants with AA, UC, and CU, instead of GA, in the second and third positions and measured their relative efficiencies in the reading of codons starting with UU, GA, and AG, respectively. The specificity concerning the wobble position is essentially the same as that in the case of the codons starting with UC.

Inferring the conformation of RNA base pairs and triples from patterns of sequence variation

The success of comparative analysis in resolving RNA secondary structure and numerous tertiary interactions relies on the presence of base covariations. Although the majority of base covariations in aligned sequences is associated to Watson-Crick base pairs, many involve non-canonical or restricted base pair exchanges (e.g. only G:C/A:U), reflecting more specific structural constraints. We have developed a computer program that determines potential base pairing conformations for a given set of paired nucleotides in a sequence alignment. This program (ISOPAIR) assumes that the base pair conformation is maintained through sequence variation without significantly affecting the path of the sugar-phosphate backbone. ISOPAIR identifies such 'isomorphic' structures for any set of input base pair or base triple sequences. The program was applied to base pairs and triples with known structures and sequence exchanges. In several instances, isomorphic structures were correctly identified with ISOPAIR. Thus, ISOPAIR is useful when assessing non-canonical base pair conformations in comparative analysis. ISOPAIR applications are limited to those cases where unusual base pair exchanges indeed reflect a non-canonical conformation.

Effect of the nucleotide-37 on the interaction of tRNA(Phe) with the P site of Escherichia coli ribosomes

The method of anticodon loop replacement has been used to make derivatives of yeast tRNA(Phe) with the substitution at position 37 (tRNA(Phe)GAAA) and at the anticodon(tRNA(Phe)GCAG). A quantitative study of the interaction of various types of deacylated yeast tRNA(Phe) (tRNA(Phe)+Y, tRNA(Phe)GAAA, tRNA(Phe)-Y) with the P site of the [70S ribosome*poly(U)]-complex was carried out at different Mg2+ concentrations and temperatures. The presence and nature of the nucleotide situated at the 3'-end of the anticodon are essential for such interaction in E coli ribosomes. Replacement of the Y base with the unmodified adenosine decreases the interaction enthalpy from 39 kcal/mol to 24 kcal/mol, whereas its removal reduces the interaction enthalpy to 16 kcal/mol. Replacement of the second anticodon nucleotide, adenosine, with cytosine further reduces the enthalpy to 6 kcal/mol, which is typical of tRNA-P site interaction in the absence of poly(U). In the absence of poly(U) the affinity of tRNA(PheY) for the P site of the 70S ribosome is five times lower than the affinity of tRNA(Phe+Y) or tRNA(Phe)GCAG. Thus, in the ribosome the modified nucleotide stabilizes the codon-anticodon interaction through its stacking interaction with the codon-anticodon base stack. In addition, this decreases the free energy of binding as a result of the interaction of the modified nucleotide itself with the hydrophobic center of the P site.

tRNA anticodons with the modified nucleoside 2-methylthio-N6-(4-hydroxyisopentenyl)adenosine distinguish between bases 3' of the codon

The modified nucleoside 2-methylthio-N6-(4-hydroxyisopentenyl)adenosine (ms2io6A) is present immediately to the 3' side of the anticodon (position 37) in tRNAs that read codons starting with uridine and hence include amber (UAG) suppressor tRNAs. We have used strains of Salmonella typhimurium that differ only in their ability to synthesize ms2io6A in order to determine specifically how this modified nucleoside influences the efficiency of amber suppression in two codon contexts differing by only which base is 3' of the codon. The results show that the presence of the modified nucleoside ms2io6A not only improves the efficiency of the suppressor tRNAs but also allows them to distinguish between at least two bases 3' of the codon. Thus, the presence of ms2io6A reduces the intrinsic codon context sensitivity of the tRNA and specifically counteracts an unfavourable nucleotide on the 3' side of the codon. The possible codon-anticodon interactions responsible for this effect are discussed.

Preferential alkylation by 1,3-bis(2-chloroethyl)-1-nitrosourea (BCNU) of guanines with guanines as neighboring bases in DNA

The base sequence of DNA has been shown to influence the kinds and amounts of alkylation of purine bases by N-methyl-N-nitrosourea [W. T. Briscoe and L-E. Cotter, Chem. Biol. Interact. 56, 321 (1985)]. In the present study, the alkylation of DNA polymers of defined sequence by 1,3-bis(2-chloroethyl)-1-nitrosourea (BCNU) has been investigated. The assay involved treating poly (dG).poly(dC), poly(dG-dC).poly(dG-dC), poly(dA-dC).poly(dG-dT), poly(dA-dG).poly(dC-dT), and calf thymus DNA with BCNU, followed by hydrolysis to release the modified purine bases and separation and quantitation of these by HPLC. Analysis of the results revealed that there was a 24-fold increase of 7-(beta-hydroxyethyl)guanine (HOEtG) in poly(dG).poly(dC) relative to poly(dA-dG).poly(dC-dT). There was also a 3-fold increase in HOEtG in poly(dG-dC).poly(dG-dC), poly(dA-dC).poly(dG-dT) and calf thymus DNA relative to poly(dA-dG).poly(dC-dT). A 2- to 4-fold increase of 7(beta-aminoethyl)guanine (AmEtG) was observed in poly(dG-dC).poly(dG-dC) relative to the other polymers tested. This study has determined that guanines in certain base sequences in polydeoxyribonucleotides are more susceptible to BCNU alkylation at the N-7 position than guanines in other sequences.

Molecular mechanics calculations of dA12.dT12 and of the curved molecule d(GCTCGAAAAA)4.d(TTTTTCGAGC)4

Using the AMBER software package (Weiner and Kollman 1981) substantially modified for electrostatic contributions, the structural energies of the double-stranded oligonucleotides dA12.dT12 and d(GCTCGAAAAA)4.d(TTTTTCGAGC)4 were minimized. Using various starting structures for the molecule dA12.dT12, one final structure is obtained which possesses the experimentally determined properties of poly(dA).poly(dT). This structure is an A-form-B-form-hybrid structure similar to that of Arnott et al. (1983). The dA-strand is similar to an A-form while the dT-strand is similar to normal B-form. This structure and separately optimized B-form sequence stretches were used to construct the double-stranded fragment d(GCTCGAAAAA)4 which again was optimized. This sequence, when imbedded in a DNA fragment as contiguous repeats, shows a gel migration anomaly which has been interpreted as stable curvature of the DNA (Diekmann 1986). The calculated structure of this sequence indeed has a curved helix axis and is discussed as a model for curved DNA. A theoretical formalism is presented which allows one to calculate the structural parameters of any nucleic acid double helix in two different geometrical representations. This formalism is used to determine the parameters of the base-pair orientations of the curved structure in terms of wedge as well as cylindrical parameters. In the structural model presented here, the curvature of the helix axis results from an alternation of two different DNA structures in which the base-pairs possess different angles with the helix axis ('cylinder tilt'). Resulting from geometric restraints, a negative cylinder tilt angle correlates strongly with the closing of the minor groove ('wedge roll'). The blocks with different structure are not exactly coincident with the dA5-blocks and the B-DNA stretches. Within the dA5 block, base-pair tilt and wedge roll adopt large values which proceed into the 3' flanking B-DNA sequence by about one base-pair. These properties of the structure calculated here are discussed in terms of different models explaining DNA curvature.

13C-NMR of ribosyl ApApA, ApApG and ApUpG

The chemical shifts as well as the 13C-31P coupling constants of the carbon-13 nuclei in single-stranded ApApA, ApApG, and ApUpG are sensitive to sequence and temperature. ApApA and ApApG have similar properties with large shielding (up to 1.7 ppm) of many of the base carbons upon decreasing the temperature from 70 degrees C to 11 degrees C; the base carbons have smaller shielding changes in ApUpG. Large shielding and deshielding effects are observed for the 1', 3', 4' and 5'-carbons over this temperature range. Analysis of the 13C-31P couplings measured at the 4' ribose carbons show that the population of the anti rotamer about O5'-C5' varies from 98 to 75%, and is higher in ApApA and ApApG than in ApUpG. The CCOP coupling data at 2' and 4' is consistent with a blend of the -antiperiplanar/-synclinal nonclassical rotamers about the C3'-O3' bond, varying from 89/11% in ApApG to 55/45% in ApUpG. The coupling and chemical shift data support the thesis that ApUpG is stacked much less than the other two molecules. The stacked forms of all three trinucleotides is most easily interpreted by a standard A-RNA model. It is not necessary to invoke the "bulged base" hypothesis [Lee, C.-H. and Tinoco, Jr., I. (1981) Biophysical Chemistry 1, 283-294; Lankhorst, P.P., Wille, G., van Boom., J.H., Altona, C., and Haasnoot, C.A.G. (1983) Nucleic Acids Research 11, 2839-2856] to explain the contrast in 13C spectroscopic properties of ApUpG in comparison to ApApG and ApApA.

DNA sequence has an effect on the extent and kinds of alkylation of DNA by a potent carcinogen

A system has been developed to study the effects of base sequence (neighboring bases) upon the alkylation of guanine (G) and adenine (A) bases in DNA. The study was performed on the synthetic polydeoxyribonucleotides, poly(dG).poly(dC), poly(dG-dC).poly(dG-dC), poly(dA).poly(dT), poly(dA-dT).poly(dA-dT), poly(dA-dC).poly(dG-dT), poly(dA-dG).poly(dC-dT), as well as calf thymus DNA. Each polynucleotide was treated with N-[3H]methyl-N-nitrosourea (MNU), depurinated, and the freed alkylpurines separated by HPLC and quantitated by liquid scintillation counting. The amounts of 3-methylguanine (3-MG), 7-MG, and O6-MG relative to guanine, and 3-methyladenine (3-MA) and 1-MA plus 7-MA relative to adenine, and also the O6-MG/7-MG ratios were highly reproducible for a given polynucleotide. Significant differences were found in the amounts of each of the methylpurines formed when compared among the six synthetic polynucleotides and DNA. This evidence is interpreted as an effect upon alkylation which is ultimately dependent upon the base sequence. These findings may have significance in defining the specificity of chemical carcinogens in terms of the susceptability to modification of nucleotide sequences such as those found in certain oncogenes.

The structure of yeast tRNA(Asp). A model for tRNA interacting with messenger RNA

The anticodon of yeast tRNA(Asp), GUC, presents the peculiarity to be self-complementary, with a slight mismatch at the uridine position. In the orthorhombic crystal lattice, tRNA(Asp) molecules are associated by anticodon-anticodon interactions through a two-fold symmetry axis. The anticodon triplets of symmetrically related molecules are base paired and stacked in a normal helical conformation. A stacking interaction between the anticodon loops of two two-fold related tRNA molecules also exists in the orthorhombic form of yeast tRNA(Phe). In that case however the GAA anticodon cannot be base paired. Two characteristic differences can be correlated with the anticodon-anticodon association: the distribution of temperature factors as determined from the X-ray crystallographic refinements and the interaction between T and D loops. In tRNA(Asp) T and D loops present higher temperature factors than the anticodon loop, in marked contrast to the situation in tRNA(Phe). This variation is a consequence of the anticodon-anticodon base pairing which rigidifies the anticodon loop and stem. A transfer of flexibility to the corner of the tRNA molecule disrupts the G19-C56 tertiary interactions. Chemical mapping of the N3 position of cytosine 56 and analysis of self-splitting patterns of tRNA(Asp) substantiate such a correlation.

Temperature jump relaxation studies on the interactions between transfer RNAs with complementary anticodons. The effect of modified bases adjacent to the anticodon triplet

We have used the temperature-jump relaxation technique to determine the kinetic and thermodynamic parameters for the association between the following tRNAs pairs having complementary anticodons: tRNA(Ser) with tRNA(Gly), tRNA(Cys) with tRNA(Ala) and tRNA(Trp) with tRNA(Pro). The anticodon sequence of E. coli tRNA(Ser), GGA, is complementary to the U*CC anticodon of E. coli tRNA(Gly(2] (where U* is a still unknown modified uridine base) and A37 is not modified in none of these two tRNAs. E. coli tRNA(Ala) has a VGC anticodon (V is 5-oxyacetic acid uridine) while tRNA(Cys) has the complementary GCA anticodon with a modified adenine on the 3' side, namely 2-methylthio N6-isopentenyl adenine (mS2i6A37) in E. Coli tRNA(Cys) and N6-isopentenyl adenine (i6A37) in yeast tRNA(Cys). The brewer yeast tRNA(Trp) (anticodon CmCA) differs from the wild type E. coli tRNA(Trp) (anticodon CCA) in several positions of the nucleotide sequence. Nevertheless, in the anticodon loop, only two interesting differences are present: A37 is not modified while C34 at the first anticodon position is modified into a ribose 2'-O methyl derivative (Cm). The corresponding complementary tRNA is E.coli tRNA(Pro) with the VGG anticodon. Our results indicate a dominant effect of the nature and sequence of the anticodon bases and their nearest neighbor in the anticodon loop (particularly at position 37 on the 3' side); no detectable influence of modifications in the other tRNA stems has been detected. We found a strong stabilizing effect of the methylthio group on i6A37 as compared to isopentenyl modification of the same residue. We have not been able so far to assess the effect of isopentenyl modification alone in comparison to unmodified A37. The results obtained with the complex yeast tRNA(Trp)-E.coli tRNA(Pro) also suggest that a modification of C34 to Cm34 does not significantly increase the stability of tRNA(Trp) association with its complementary anticodon in tRNA(Pro). The observations are discussed in the light of inter- and intra-strand stacking interactions among the anticodon triplets and with the purine base adjacent to them, and of possible biological implications.

Coxsackievirus B3: primary structure of the 5' non-coding and capsid protein-coding regions of the genome

The nucleotide sequence from the 5' terminus to nucleotide 3822 has been determined for the genome of the enterovirus, coxsackievirus B3 (CB3). This region encompasses the 5'-terminal 738 nucleotide non-coding sequence and the region which codes for the four viral capsid proteins and for the initial products of the P2 region. Regions in the 5' non-translated RNA sequence which may be involved with a structural and/or regulatory role are discussed.

Quantitation of base substitutions in eukaryotic 5S rRNA: selection for the maintenance of RNA secondary structure

Eukaryotic 5S rRNA sequences from 34 diverse species were compared by the following method: (1) The sequences were aligned; (2) the positions of substitutions were located by comparison of all possible pairs of sequences; (3) the substitution sites were mapped to an assumed general base pairing model; and (4) the R-Y model of base stacking was used to study stacking pattern relationships in the structure. An analysis of the sequence and structure variability in each region of the molecule is presented. It was found that the degree of base substitution varies over a wide range, from absolute conservation to occurrence of over 90% of the possible observable substitutions. The substitutions are located primarily in stem regions of the 5S rRNA secondary structure. More than 88% of the substitutions in helical regions maintain base pairing. The disruptive substitutions are primarily located at the edges of helical regions, resulting in shortening of the helical regions and lengthening of the adjacent nonpaired regions. Base stacking patterns determined by the R-Y model are mapped onto the general secondary structure. Intrastrand and interstrand stacking could stabilize alternative coaxial structures and limit the conformational flexibility of nonpaired regions. Two short contiguous regions are 100% conserved in all species. This may reflect evolutionary constraints imposed at the DNA level by the requirement for binding of a 5S gene transcription initiation factor during gene expression.

Relative stability of guanosine-cytidine diribonucleotide cores: a 1H NMR assessment

Proton NMR was used to study the secondary structure and melting behavior of six self-complementary oligoribonucleotide tetramers, each containing two guanosine and two cytidine residues (GGCC, CCGG, GCCG, CGGC, GCGC, and CGCG). GGCC and CCGG formed perfect duplexes containing four G.C base pairs with Tms of 54 and 47.8 degrees C, respectively; GCCG and CGGC formed staggered duplexes with two G.C base pairs and four 3' double-dangling bases, with Tms of 35.5 and 29.2 degrees C, respectively; GCGC formed a perfect duplex with a Tm of 49.9 degrees C, while CGCG formed a staggered duplex with a Tm of 36.9 degrees C. From these results, an order of stability of the cores containing two G.C base pairs was proposed: GC:GC is more stable than GG:CC which is more stable than CG:CG. The RY model for secondary structure stability prediction was applied to the above tetramers with reasonable success. Suggestions for refinements are discussed.

Conformation and dynamics of short DNA duplexes: (dC-dG)3 and (dC-dG)4

Natural abundance 13C NMR spectra of duplexed (dC-dG)3 and (dC-dG)4 exhibit resolved resonances for most of the carbons at 0.1M NaCl in aqueous solution. Large transitions in chemical shift for many of the hexamer carbons (up to 1.8 ppm) are observed in variable temperature measurements. Determination of spin-lattice relaxation times and nuclear Overhauser enhancements in 0.1M NaCl indicate that the duplexes tumble almost isotropically, with overall correlation times near 5 nsec; the sugar carbons experience more rapid local motions than do the base carbons. The relaxation data are also consistent with the most rapid local motions occurring at the chain-terminal residues, especially in the Cyd(1) sugar. 4M NaCl causes changes in the 13C chemical shifts of most of the guanine base carbons, and rearrangements in the deoxyribose carbon shifts; this is consistent with changes predicted by a salt-induced B to Z transition, viz. conversion of the guanylates from the anti to syn range about the glycosyl bond, and from the S to N pseudorotational state of the deoxyribose ring.