Free energy of imperfect nucleic acid helices. II. Small hairpin loops

Physical studies of enzymically synthesized oligonucleotides of defined sequence are used to evaluate quantitatively the stability of small RNA hairpin loops and helices. The series (Ap)₄G(pC) _N(pU)₄, N = 4, 5 or 6, exists as monomolecular hairpin helices when N ≥ 5, and as imperfect dimer helices when N ≤ 4. In this size range, hairpin loops become more favorable (less destabilizing thermodynamically) as they increase in size from 3 to 4 to 5 unbonded nucleotides. Very small hairpin loops are particularly destabilizing; molecules whose base sequence would imply a hairpin loop of three nucleotides will generally exist with a loop of five, including a broken terminal base pair.Thermodynamic parameters for base pair and loop formation are calculated by a method which makes unnecessary the use of measured enthalpies of polynucleotide melting. Literature data on oligonucleotide double helices yield estimates of the free energy contribution from each of the six types of stacking interactions between three possible neighboring base pairs. The advantage of this approach is that the properties of oligonucleotides are used in predicting the stability of small RNA helices, avoiding the long extrapolation from the properties of high polymers.We provide Tables of temperature-dependent free energies that allow one to predict the stability and thermal transition temperature of many simple RNA secondary structures (applicable to ~1 m-Na⁺ concentration). As an example, we apply the rules to an isolated fragment of tRNA^Ser (yeast) (Coutts, 1971), whose properties were not used in calculating the free-energy parameters. The experimental melting temperature of 88 °C is predicted with an error margin of 5 deg. C.     

Enzyme action at 3' termini of ionizing radiation-induced DNA strand breaks

gamma-Irradiation of DNA in vitro produces two types of single strand breaks. Both types of strand breaks contain 5'-phosphate DNA termini. Some strand breaks contain 3'-phosphate termini, some contain 3'-phosphoglycolate termini (Henner, W.D., Rodriguez, L.O., Hecht, S. M., and Haseltine, W. A. (1983) J. Biol. Chem. 258, 711-713). We have studied the ability of prokaryotic enzymes of DNA metabolism to act at each of these types of gamma-ray-induced 3' termini in DNA. Neither strand breaks that terminate with 3'-phosphate nor 3'-phosphoglycolate are substrates for direct ligation by T4 DNA ligase. Neither type of gamma-ray-induced 3' terminus can be used as a primer for DNA synthesis by either Escherichia coli DNA polymerase or T4 DNA polymerase. The 3'-phosphatase activity of T4 polynucleotide kinase can convert gamma-ray-induced 3'-phosphate but not 3'-phosphoglycolate termini to 3'-hydroxyl termini that can then serve as primers for DNA polymerase. E. coli alkaline phosphatase is also unable to hydrolyze 3'-phosphoglycolate groups. The 3'-5' exonuclease actions of E. coli DNA polymerase I and T4 DNA polymerase do not degrade DNA strands that have either type of gamma-ray-induced 3' terminus. E. coli exonuclease III can hydrolyze DNA with gamma-ray-induced 3'-phosphate or 3'-phosphoglycolate termini or with DNase I-induced 3'-hydroxyl termini. The initial action of exonuclease III at 3' termini of ionizing radiation-induced DNA fragments is to remove the 3' terminal phosphate or phosphoglycolate to yield a fragment of the same nucleotide length that has a 3'-hydroxyl terminus. These results suggest that repair of ionizing radiation-induced strand breaks may proceed via the sequential action of exonuclease, DNA polymerase, and DNA ligase. The possible role of exonuclease III in repair of gamma-radiation-induced strand breaks is discussed.     

DNA sequence of the 16S rRNA/23S rRNA intercistronic spacer of two rDNA operons of the archaebacterium Methanococcus vannielii

The DNA sequence of the spacer (plus flanking) regions separating the 16S rRNA and 23S rRNA genes of two presumptive rDNA operons of the archaebacterium Methanococcus vannielii was determined. The spacers are 156 and 242 base pairs in size and they share a sequence homology of 49 base pairs following the 3' terminus of the 16S rRNA gene and of about 60 base pairs preceding the 5' end of the 23S rRNA gene. The 242 base pair spacer, in addition contains a sequence which can be transcribed into tRNAAla, whereas no tRNA-like secondary structure can be delineated from the 156 base pair spacer region. Almost complete sequence homology was detected between the end of the 16S rRNA gene and the 3' termini of either Escherichia coli or Halobacterium halobium 16S rRNA, whereas the putative 5' terminal 23S rRNA sequence shared partial homology with E. coli 23S rRNA and eukaryotic 5.8S rRNA.     

Free energy of imperfect nucleic acid helices. 3. Small internal loops resulting from mismatches

Physical studies of enzymioally synthesized oligoribonucleotides of defined sequence are used to evaluate quantitatively the destabilizing influence of mismatched bases in a double helix. The series (A-)₄G(-C)_n(-U)₄, N = 1 to 6, exist as imperfect dimer helices when N is equal to or less than 4, and as monomolecular hairpin helices when N is 5 and 6. Internal loops become progressively more destabilizing as their size increases from 2 to 4 to 6 nucleotides resulting from 1, 2 and 3 consecutive mismatched base pairs. However, the stability of a helix will generally be greater if a given number of mismatched pairs occur consecutively rather than in isolation from one another.These data may be used for improved calculations of stability of RNA secondary structure, to estimate the frequency of structural fluctuations in a double helix and to assess the stability of modified polynucleotide helices. An unmodified double helix of one million randomly arranged base pairs should contain on the time average approximately 10 G.C and 500 A.U pairs in non-hydrogen bonded, unstacked conformations at 25 °C. Our estimate of the effect of mismatching on T_m values of high polymers is less precise because of the long temperature extrapolation required. However, we estimate that DNA or RNA treated with mutagens which interrupt up to 20% of the nucleotide pairs will show a drop of about 1.2 deg. C in melting temperature with each unit per cent of modification.     

Structural heterogeneity in intramolecular DNA triple helices

Oligodeoxynucleotides designed to form intramolecular triple helices are widely used as model systems in thermodynamic and structural studies. We now report results from UV, Raman and NMR experiments demonstrating that the strand polarity, which also determines the orientation of the connecting loops, has a considerable impact on the formation and stability of pyr x pur x pyr triple helices. There are two types of monomolecular triplexes that can be defined by the location of their purine tract at either the 5'- or 3'-end of the sequence. We have examined four pairs of oligonucleotides with the same base composition but with reversed polarity that can fold into intramolecular triple helices with seven base triplets and two T4 loops under appropriate conditions. UV spectroscopic monitoring of thermal denaturation indicates a consistently higher thermal stability for the 5'-sequences at pH 5.0 in the absence of Mg2+ ions. Raman spectra provide evidence for the formation of triple helices at pH 5 for oligomers with purine tracts located at either the 5'- or 3'-end of the sequence. However, NMR measurements reveal considerable differences in the secondary structures formed by the two types of oligonucleotides. Thus, at acidic pH significant structural heterogeneity is observed for the 3'-sequences. Employing selectively 15N-labeled oligomers, NMR experiments indicate a folding pattern for the competing structures that at least partially changes both Hoogsteen and Watson-Crick base-base interactions.     

Escherichia coli 5S RNA A and B conformers. Characterisation by enzymatic and chemical methods

The structures of the two stable conformers of Escherichia coli 5 S RNA, the and B form, were compared. Information about the structures were obtained using the methods of limited enzymatic hydrolysis and chemical modification of accessible nucleotides. Base-specific modifications were performed for adenosines and cytidines using diethylpyrocarbonate and dimethylsulfate in combination with a strand-scission reaction at the modified site. Base-specific (RNase T1) as well as conformation-specific (nuclease S1, cobra venom nuclease) enzymes were employed for the limited enzymatic hydrolysis. Clear differences in the accessibility of the two 5 S RNA conformers to the enzymes and the chemical reagents were established and the regions with altered reactivities were localized in the 5 S RNA structure. The results are consistent with the disruption of the secondary structural interactions in helix II and partly in helices III and IV during the transition from the A to the B form. (The numbering of the helices is according to the generally accepted Fox and Woese model.) In addition some regions presumably involved in the tertiary structure are distorted. There is evidence, however, for the new formation of structural regions between two distant sites in the 5 S RNA B form. The results enable us to refine the existing 5 S RNA A-form model and provide insight into the structural dynamics that lead to the formation of the 5 S RNA B form.     

A small angle x-ray scattering study of a fragment derived from E. coli 5S RNA.

The structure of a 62 base nuclease resistant fragment of E. coli 5S RNA (bases 1-11, 69-87, 89-120) has been examined by small angle x-ray scattering. The results obtained are indistinguishable from those expected if this oligonucleotide complex were a perfect RNA double helix of about 30 base pairs. These results indicate that this portion of 5S RNA is in a configuration which is approximately double helical, even though proper base pairing is possible over only half its length.     

Properties of cloned and expressed human RNase H1.

We have characterized cloned His-tag human RNase H1. The activity of the enzyme exhibited a bell-shaped response to divalent cations and pH. The optimum conditions for catalysis consisted of 1 mM Mg(2+) and pH 7-8. In the presence of Mg(2+), Mn(2+) was inhibitory. Human RNase H1 shares many enzymatic properties with Escherichia coli RNase H1. The human enzyme cleaves RNA in a DNA-RNA duplex resulting in products with 5'-phosphate and 3'-hydroxy termini, can cleave overhanging single strand RNA adjacent to a DNA-RNA duplex, and is unable to cleave substrates in which either the RNA or DNA strand has 2' modifications at the cleavage site. Human RNase H1 binds selectively to "A-form"-type duplexes with approximately 10-20-fold greater affinity than that observed for E. coli RNase H1. The human enzyme displays a greater initial rate of cleavage of a heteroduplex-containing RNA-phosphorothioate DNA than an RNA-DNA duplex. Unlike the E. coli enzyme, human RNase H1 displays a strong positional preference for cleavage, i.e. it cleaves between 8 and 12 nucleotides from the 5'-RNA-3'-DNA terminus of the duplex. Within the preferred cleavage site, the enzyme displays modest sequence preference with GU being a preferred dinucleotide. The enzyme is inhibited by single-strand phosphorothioate oligonucleotides and displays no evidence of processivity. The minimum RNA-DNA duplex length that supports cleavage is 6 base pairs, and the minimum RNA-DNA "gap size" that supports cleavage is 5 base pairs.     

Identification and assignment of base pairs in three helical stems of Torulopsis utilis ribosomal 5S RNA and its RNase T1 and RNase T2 cleaved fragments via 500-MHz proton homonuclear overhauser enhancements

Imino proton resonances in the downfield region (10-14 ppm) of the 500-MHz 1H NMR spectrum of Torulopsis utilis 5S RNA are identified (A X U, G X C, or G X U) and assigned to base pairs in helices I, IV, and V via analysis of homonuclear Overhauser enhancements (NOE) from intact T. utilis 5S RNA, its RNase T1 and RNase T2 digested fragments, and a second yeast (Saccharomyces cerevisiae) 5S RNA whose nucleotide sequence differs at only six residues from that of T. utilis 5S RNA. The near-identical chemical shifts and NOE behavior of most of the common peaks from these four RNAs strongly suggest that helices I, IV, and V retain the same conformation after RNase digestion and that both T. utilis and S. cerevisiae 5S RNAs share a common secondary and tertiary structure. Of the four G X U base pairs identified in the intact 5S RNA, two are assigned to the terminal stem (helix I) and the other two to helices IV and V. Seven of the nine base pairs of the terminal stem have been assigned. Our experimental demonstration of a G X U base pair in helix V supports the 5S RNA secondary structural model of Luehrsen and Fox [Luehrsen, K. R., & Fox, G.E. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 2150-2154]. Finally, the base-pair proton peak assigned to the terminal G X U in helix V of the RNase T2 cleaved fragment is shifted downfield from that in the intact 5S RNA, suggesting that helices I and V may be coaxial in intact T. utilis 5S RNA.     

Structural analysis of three prokaryotic 5S rRNA species and selected 5S rRNA--ribosomal-protein complexes by means of Pb(II)-induced hydrolysis.

Lead ions have been applied to the structural analysis of 5S rRNA from Thermus thermophilus, Bacillus stearothermophilus and Escherichia coli. Based on the distribution of Pb(II)-induced cleavages, some minor modifications of the consensus secondary structure model of 5S rRNA are proposed. They include the possible base pairing between nucleotides at positions 11 and 109, as well as changes in secondary interactions within the helix B region. The 'prokaryotic arm' region is completely resistant to hydrolysis in the three RNA species, suggesting that it is a relatively stable, highly ordered structure. Hydrolysis of E. coli 5S rRNA complexed with ribosomal protein L18 shows, besides the shielding effect of the bound protein, a highly enhanced cleavage between A108 and A109. It supports the concept that the major L18-induced conformational change involves the junction of helices A, B and D.     

Studies On The Cross Strand Hydrogen Bonds In DNA Double Helices

Abstract

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross- strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may infact be considered as yet another general structural feature of DNA helices.     

Studies on the cross stand hydrogen bonds in DNA double helices.

A systematic analysis has been carried out to examine all the stereochemically possible bifurcated hydrogen bonds including those of cross strand type between propeller twisted base pairs in DNA double helices by stereochemical considerations involving base pairs alone and by molecular mechanics studies on dimer and trimer duplexes. The results show that there are limited number of combinations of adjacent base pairs that would facilitate bifurcated cross-strand hydrogen bond (CSH). B-type helices concomitant with negative propeller twist seem to be more favored for the occurrence of CSH than canonical A-type helices because of slide in the latter. The results also demonstrate that helices with appropriate sequences may possess continuous run of these propeller twist driven cross strand hydrogen bonds indicating that they may in fact be considered as yet another general structural feature of DNA helices.     

Characterization of the gene encoding glutamine synthetase I (glnA) from Bradyrhizobium japonicum.

We have isolated the Bradyrhizobium japonicum gene encoding glutamine synthetase I (glnA) from a phage lambda library by using a fragment of the Escherichia coli glnA gene as a hybridization probe. The rhizobial glnA gene has homology to the E. coli glnA gene throughout the entire length of the gene and can complement an E. coli glnA mutant when borne on an expression plasmid in the proper orientation to be transcribed from the E. coli lac promoter. High levels of glutamine synthetase activity can be detected in cell-free extracts of the complemented E. coli. The enzyme encoded by the rhizobial gene was identified as glutamine synthetase I on the basis of its sedimentation properties and resistance to heat inactivation. DNA sequence analysis predicts a high level of amino acid sequence homology among the amino termini of B. japonicum, E. coli, and Anabaena sp. strain 7120 glutamine synthetases. S1 nuclease protection mapping indicates that the rhizobial gene is transcribed from a single promoter 131 +/- 2 base pairs upstream from the initiation codon. This glnA promoter is active when B. japonicum is grown both symbiotically and in culture with a variety of nitrogen and carbon sources. There is no detectable sequence homology between the constitutively expressed glnA promoter and the differentially regulated nif promoters of the same B. japonicum strain.     

An estimate of the nearest neighbor base-pair content of 5S RNA using CD and absorption spectroscopy.

We analyzed the CD and uv absorption spectra of 5S RNA from Escherichia coli using the method developed in the preceding paper. The analysis of spectra of 5S RNA at 20 degrees C in 0.1M NaClO4, 2.5 mM Na+ (phosphate), pH 7.0, and 0.5 mM MgSO4 gave 7 +/- 3.6 A.U base pairs, 25 +/- 3.6 G.C base pairs, and 7.5 +/- 3.6 G.U base pairs. Estimates of nearest neighbor base pairs were more consistent with the Pieler-Erdmann and the Gewirth-Moore secondary structure models than with the Fox-Woese or the Burns-Luoma-Marshall models. We also examined the structure of 5S RNA as a function of temperature. The melting profile exhibited two transitions--one at about 35 degrees C and one above 50 degrees C. Our spectral data showed that helices I and II were stable during the first transition, and agreed with other data that helix III was the most likely helix to have melted. The results from this in-depth study of 5S RNA indicate that our method of analysis should be useful for studying the secondary structures of other small, unmodified RNAs.     

Structural clues in the sequences of the aquaporins

The large number of sequences available for the aquaporin family represents a valuable source of information to incorporate into three-dimensional structure determination. Phylogenetic analysis was used to define type sequences to avoid extreme over-representation of some subfamilies, and as a measure of the quality of multiple sequence alignment. Inspection of the sequence alignment suggested eight conserved segments that define the core architecture of six transmembrane helices and two functional loops, B and E, projecting into the plane of the membrane. The sum of the core segments and the minimum lengths of the interlinking loops constitute the 208 residues necessary to satisfy the aquaporin architecture. Analysis of hydrophobic and conservation periodicity and of correlated mutations across the alignment indicated the likely assignment and orientation of the helices in the bilayer. This assignment is examined with respect to the structure of the erythrocyte aquaporin 1 determined by electron crystallography. The aquaporin 1 tetramer is described as three rings of helices, each ring with a different exposure to the lipid environment. The sequence analysis clearly suggests that two helices are exposed along their whole lengths, two helices are exposed only at their N termini, and two helices are not exposed to lipid. It is further proposed that, besides loops B and E, the highly conserved motifs on helices 1 and 4, ExxxTxxF/L, could line the water channel.