Search

Magnesium, the most abundant divalent cation in cells, catalyzes RNA cleavage but also promotes RNA folding. Because folding can protect RNA from cleavage, we predicted a 'Goldilocks landscape', with local maximum in RNA lifetime at Mg2+ concentrations required for folding. Here, we use simulation and experiment to discover an innate and sophisticated mechanism of control of RNA lifetime. By simulation we characterized RNA Goldilocks landscapes and their dependence on cleavage and folding parameters. Experiments with yeast tRNAPhe and the Tetrahymena ribozyme P4-P6 domain show that structured RNAs can inhabit Goldilocks peaks. The Goldilocks peaks are tunable by differences in folded and unfolded cleavage rate constants, Mg2+ binding cooperativity, and Mg2+ affinity. Different folding and cleavage parameters produce Goldilocks landscapes with a variety of features. Goldilocks behavior allows ultrafine control of RNA chemical lifetime, whereas non-folding RNAs do not display Goldilocks peaks of protection. In sum, the effects of Mg2+ on RNA persistence are expected to be pleomorphic, both protecting and degrading RNA. In evolutionary context, Goldilocks behavior may have been a selectable trait of RNA in an early Earth environment containing Mg2+ and other metals., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

rRNA is the single most abundant polymer in most cells. Mammalian rRNAs are nearly twice as large as those of prokaryotes. Differences in rRNA size are due to expansion segments, which contain extended tentacles in metazoans. Here we show that the terminus of an rRNA tentacle of Homo sapiens contains 10 tandem G-tracts that form highly stable G-quadruplexes in vitro. We characterized rRNA of the H. sapiens large ribosomal subunit by computation, circular dichroism, UV melting, fluorescent probes, nuclease accessibility, electrophoretic mobility shifts, and blotting. We investigated Expansion Segment 7 (ES7), oligomers derived from ES7, intact 28S rRNA, 80S ribosomes, and polysomes. We used mass spectrometry to identify proteins that bind to rRNA G-quadruplexes in cell lysates. These proteins include helicases (DDX3, CNBP, DDX21, DDX17) and heterogeneous nuclear ribonucleoproteins. Finally, by multiple sequence alignments, we observe that G-quadruplex-forming sequences are a general feature of LSU rRNA of Chordata but not, as far as we can tell, of other species. Chordata ribosomes present polymorphic tentacles with the potential to switch between inter- and intramolecular G-quadruplexes. To our knowledge, G-quadruplexes have not been reported previously in ribosomes., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Cisplatin is an important anticancer drug with a complex mode of action, a variety of possible targets, and numerous resistance mechanisms. While genomic DNA has traditionally been considered to be its most critical anticancer target, several lines of evidence suggest that various RNAs and other biomolecules may play a role in its anticancer mode of action. In this report we demonstrate that cisplatin modifies pre-miR-200b, impairs its processing to mature miRNA, and decreases miR-200b expression in ovarian cancer cells. Considering the role of miR-200b in epithelial-to-mesenchymal transition and cancer chemosensitivity, cisplatin-induced modification of pre-miR-200b and subsequent deregulation of mature miR-200b may, depending on cell context, limit anticancer activity of this important anticancer drug. More gener- ally, precursor miRNAs may be important targets of cisplatin and play a role in this drug's anticancer activity or modulate cell responses to this drug.

Many experiments demonstrate that regions with higher GC-content in natural DNAs unwind at higher temperatures adsorbing more heat than equivalently sized regions with lower GC-content. This simple observation implies that normalized calorimetric melting profiles (calorimetric cDMCs) will not be equivalent differential melting curves (DMCs). We propose simple expressions for long natural and random DNA sequences to reciprocally convert DMCs and corresponding calorimetric cDMCs. The expressions are confirmed by the Poland-Fixman-Freire method and an approach based upon mixtures of homopolymeric duplexes. Using these expressions and experimental calorimetric data, we demonstrate that the average relative deviation between DMC and cDMC is proportional to the temperature melting range of the helix-coil transition ΔT. Corresponding difference between melting temperatures is proportional to ΔT 2 . In general, sequence and ionic conditions influence the deviation through their effect on ΔT. On the basis of the developed approach, we propose a method to determine the thermodynamic melting temperature (ratio of calorimetric enthalpy and entropy of the helix-coil transition) for natural DNAs from optical DMCs without calorimetric experiments., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Divergence between prokaryotic and eukaryotic ribosomal RNA (rRNA) and among eukaryotic ribosomal RNAs is focused in expansion segments (ESs). Eukaryotic ribosomes are significantly larger than prokaryotic ribosomes partly because of their ESs. We hypothesize that larger rRNAs of complex organisms could confer increased functionality to the ribosome. Here, we characterize the binding partners of Saccharomyces cerevisiae expansion segment 7 (ES7), which is the largest and most variable ES of the eukaryotic large ribosomal subunit and is located at the surface of the ribosome. In vitro RNA-protein pull-down experiments using ES7 as a bait indicate that ES7 is a binding hub for a variety of non-ribosomal proteins essential to ribosomal function in eukaryotes. ES7-associated proteins observed here cluster into four groups based on biological process, (i) response to abiotic stimulus (e.g., response to external changes in temperature, pH, oxygen level, etc.), (ii) ribosomal large subunit biogenesis, (iii) protein transport and localization, and (iv) transcription elongation. Seven synthetases, Ala-, Arg-, Asp-, Asn-, Leu-, Lys- and TyrRS, appear to associate with ES7. Affinities of AspRS, TyrRS and LysRS for ES7 were confirmed by in vitro binding assays. The results suggest that ES7 in S. cerevisiae could play a role analogous to the multi-synthetase complex present in higher order organisms and could be important for the appropriate function of the ribosome. Thermal denaturation studies and footprinting experiments confirm that isolated ES7 is stable and maintains a near-native secondary and tertiary structure., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The interstitial liquid phase within frozen aqueous solutions is an environment that minimizes RNA degradation and facilitates reactions that may have relevance to the RNA World hypothesis. Previous work has shown that frozen solutions support condensation of activated nucleotides into RNA oligomers, RNA ligation by the hairpin ribozyme, and RNA synthesis by a RNA polymerase ribozyme. In the current study, we examined the activity of a hammerhead ribozyme (HHR) in frozen solution. The Schistosoma mansoni hammerhead ribozyme, which predominantly cleaves RNA, can ligate its cleaved products (P1 and P2) with yields up to ~23 % in single turnover experiments at 25 °C in the presence of Mg(2+). Our studies show that this HHR ligates RNA oligomers in frozen solution in the absence of divalent cations. Citrate and other anions that exhibit strong ion-water affinity enhanced ligation. Yields up to 43 % were observed in one freeze-thaw cycle and a maximum of 60 % was obtained after several freeze-thaw cycles using wild-type P1 and P2. Truncated and mutated P1 substrates were ligated to P2 with yields of 14-24 % in one freeze-thaw cycle. A pool of P2 substrates with mixtures of all four bases at five positions were ligated with P1 in frozen solution. High-throughput sequencing indicated that 70 of the 1024 possible P2 sequences were represented in ligated products at 1000 or more read counts per million reads. The results indicate that the HHR can ligate a range of short RNA oligomers into an ensemble of diverse sequences in ice.

We present a molecular-level model for the origin and evolution of the translation system, using a 3D comparative method. In this model, the ribosome evolved by accretion, recursively adding expansion segments, iteratively growing, subsuming, and freezing the rRNA. Functions of expansion segments in the ancestral ribosome are assigned by correspondence with their functions in the extant ribosome. The model explains the evolution of the large ribosomal subunit, the small ribosomal subunit, tRNA, and mRNA. Prokaryotic ribosomes evolved in six phases, sequentially acquiring capabilities for RNA folding, catalysis, subunit association, correlated evolution, decoding, energy-driven translocation, and surface proteinization. Two additional phases exclusive to eukaryotes led to tentacle-like rRNA expansions. In this model, ribosomal proteinization was a driving force for the broad adoption of proteins in other biological processes. The exit tunnel was clearly a central theme of all phases of ribosomal evolution and was continuously extended and rigidified. In the primitive noncoding ribosome, proto-mRNA and the small ribosomal subunit acted as cofactors, positioning the activated ends of tRNAs within the peptidyl transferase center. This association linked the evolution of the large and small ribosomal subunits, proto-mRNA, and tRNA.

The origins and evolution of the ribosome, 3-4 billion years ago, remain imprinted in the biochemistry of extant life and in the structure of the ribosome. Processes of ribosomal RNA (rRNA) expansion can be "observed" by comparing 3D rRNA structures of bacteria (small), yeast (medium), and metazoans (large). rRNA size correlates well with species complexity. Differences in ribosomes across species reveal that rRNA expansion segments have been added to rRNAs without perturbing the preexisting core. Here we show that rRNA growth occurs by a limited number of processes that include inserting a branch helix onto a preexisting trunk helix and elongation of a helix. rRNA expansions can leave distinctive atomic resolution fingerprints, which we call "insertion fingerprints." Observation of insertion fingerprints in the ribosomal common core allows identification of probable ancestral expansion segments. Conceptually reversing these expansions allows extrapolation backward in time to generate models of primordial ribosomes. The approach presented here provides insight to the structure of pre-last universal common ancestor rRNAs and the subsequent expansions that shaped the peptidyl transferase center and the conserved core. We infer distinct phases of ribosomal evolution through which ribosomal particles evolve, acquiring coding and translocation, and extending and elaborating the exit tunnel.

Accurate secondary structures are important for understanding ribosomes, which are extremely large and highly complex. Using 3D structures of ribosomes as input, we have revised and corrected traditional secondary (2°) structures of rRNAs. We identify helices by specific geometric and molecular interaction criteria, not by co-variation. The structural approach allows us to incorporate non-canonical base pairs on parity with Watson-Crick base pairs. The resulting rRNA 2° structures are up-to-date and consistent with three-dimensional structures, and are information-rich. These 2° structures are relatively simple to understand and are amenable to reproduction and modification by end-users. The 2° structures made available here broadly sample the phylogenetic tree and are mapped with a variety of data related to molecular interactions and geometry, phylogeny and evolution. We have generated 2° structures for both large subunit (LSU) 23S/28S and small subunit (SSU) 16S/18S rRNAs of Escherichia coli, Thermus thermophilus, Haloarcula marismortui (LSU rRNA only), Saccharomyces cerevisiae, Drosophila melanogaster, and Homo sapiens. We provide high-resolution editable versions of the 2° structures in several file formats. For the SSU rRNA, the 2° structures use an intuitive representation of the central pseudoknot where base triples are presented as pairs of base pairs. Both LSU and SSU secondary maps are available (http://apollo.chemistry.gatech.edu/RibosomeGallery). Mapping of data onto 2° structures was performed on the RiboVision server (http://apollo.chemistry.gatech.edu/RiboVision).

We present a de novo re-determination of the secondary (2°) structure and domain architecture of the 23S and 5S rRNAs, using 3D structures, determined by X-ray diffraction, as input. In the traditional 2° structure, the center of the 23S rRNA is an extended single strand, which in 3D is seen to be compact and double helical. Accurately assigning nucleotides to helices compels a revision of the 23S rRNA 2° structure. Unlike the traditional 2° structure, the revised 2° structure of the 23S rRNA shows architectural similarity with the 16S rRNA. The revised 2° structure also reveals a clear relationship with the 3D structure and is generalizable to rRNAs of other species from all three domains of life. The 2° structure revision required us to reconsider the domain architecture. We partitioned the 23S rRNA into domains through analysis of molecular interactions, calculations of 2D folding propensities and compactness. The best domain model for the 23S rRNA contains seven domains, not six as previously ascribed. Domain 0 forms the core of the 23S rRNA, to which the other six domains are rooted. Editable 2° structures mapped with various data are provided (http://apollo.chemistry.gatech.edu/RibosomeGallery).

Mg(2+) is essential for RNA folding and catalysis. However, for the first 1.5 billion years of life on Earth RNA inhabited an anoxic Earth with abundant and benign Fe(2+). We hypothesize that Fe(2+) was an RNA cofactor when iron was abundant, and was substantially replaced by Mg(2+) during a period known as the 'great oxidation', brought on by photosynthesis. Here, we demonstrate that reversing this putative metal substitution in an anoxic environment, by removing Mg(2+) and replacing it with Fe(2+), expands the catalytic repertoire of RNA. Fe(2+) can confer on some RNAs a previously uncharacterized ability to catalyse single-electron transfer. We propose that RNA function, in analogy with protein function, can be understood fully only in the context of association with a range of possible metals. The catalysis of electron transfer, requisite for metabolic activity, may have been attenuated in RNA by photosynthesis and the rise of O2.

Ancient components of the ribosome, inferred from a consensus of previous work, were constructed in silico, in vitro and in vivo. The resulting model of the ancestral ribosome presented here incorporates ∼20% of the extant 23S rRNA and fragments of five ribosomal proteins. We test hypotheses that ancestral rRNA can: (i) assume canonical 23S rRNA-like secondary structure, (ii) assume canonical tertiary structure and (iii) form native complexes with ribosomal protein fragments. Footprinting experiments support formation of predicted secondary and tertiary structure. Gel shift, spectroscopic and yeast three-hybrid assays show specific interactions between ancestral rRNA and ribosomal protein fragments, independent of other, more recent, components of the ribosome. This robustness suggests that the catalytic core of the ribosome is an ancient construct that has survived billions of years of evolution without major changes in structure. Collectively, the data here support a model in which ancestors of the large and small subunits originated and evolved independently of each other, with autonomous functionalities.

The three-dimensional structure of the ribosomal large subunit (LSU) reveals a single morphological element, although the 23S rRNA is contained in six secondary structure domains. Based upon maps of inter- and intra-domain interactions and proposed evolutionary pathways of development, we hypothesize that Domain III is a truly independent structural domain of the LSU. Domain III is primarily stabilized by intra-domain interactions, negligibly perturbed by inter-domain interactions, and is not penetrated by ribosomal proteins or other rRNA. We have probed the structure of Domain III rRNA alone and when contained within the intact 23S rRNA using SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension), in the absence and presence of magnesium. The combined results support the hypothesis that Domain III alone folds to a near-native state with secondary structure, intra-domain tertiary interactions, and inter-domain interactions that are independent of whether or not it is embedded in the intact 23S rRNA or within the LSU. The data presented support previous suggestions that Domain III was added relatively late in ribosomal evolution.

Mg²⁺ shares a distinctive relationship with RNA, playing important and specific roles in the folding and function of essentially all large RNAs. Here we use theory and experiment to evaluate Fe²⁺ in the absence of free oxygen as a replacement for Mg²⁺ in RNA folding and catalysis. We describe both quantum mechanical calculations and experiments that suggest that the roles of Mg²⁺ in RNA folding and function can indeed be served by Fe²⁺. The results of quantum mechanical calculations show that the geometry of coordination of Fe²⁺ by RNA phosphates is similar to that of Mg²⁺. Chemical footprinting experiments suggest that the conformation of the Tetrahymena thermophila Group I intron P4-P6 domain RNA is conserved between complexes with Fe²⁺ or Mg²⁺. The catalytic activities of both the L1 ribozyme ligase, obtained previously by in vitro selection in the presence of Mg²⁺, and the hammerhead ribozyme are enhanced in the presence of Fe²⁺ compared to Mg²⁺. All chemical footprinting and ribozyme assays in the presence of Fe²⁺ were performed under anaerobic conditions. The primary motivation of this work is to understand RNA in plausible early earth conditions. Life originated during the early Archean Eon, characterized by a non-oxidative atmosphere and abundant soluble Fe²⁺. The combined biochemical and paleogeological data are consistent with a role for Fe²⁺ in an RNA World. RNA and Fe²⁺ could, in principle, support an array of RNA structures and catalytic functions more diverse than RNA with Mg²⁺ alone.

The Escherichia coli RNA binding protein Hfq plays an important role in regulating mRNA translation through its interactions with small non-coding RNAs (sRNAs) and specific mRNAs sites. The rpoS mRNA, which codes for a transcription factor, is regulated by several sRNAs. DsrA and RprA enhance translation by pairing to a site on this mRNA, while OxyS represses rpoS mRNA translation. To better understand how Hfq interacts with these sRNAs and rpoS mRNA, the binding of wt Hfq and eleven mutant Hfqs to DsrA, RprA, OxyS and rpoS mRNA was examined. Nine of the mutant Hfq had single-residue mutations located on the proximal, distal, and outer-edge surfaces of the Hfq hexamer, while two Hfq had truncated C-terminal ends. Hfq with outer-edge mutations and truncated C-terminal ends behaved similar to wt Hfq with regard to binding the sRNAs, rpoS mRNA segments, and stimulating DsrA•rpoS mRNA formation. Proximal surface mutations decreased Hfq binding to the three sRNAs and the rpoS mRNA segment containing the translation initiation region. Distal surface mutations lowered Hfq's affinity for the rpoS mRNA segment containing the (ARN)(4) sequence. Strong Hfq binding to both rpoS mRNA segments appears to be needed for maximum enhancement of DsrA•rpoS mRNA annealing. OxyS bound tightly to Hfq but exhibited weak affinity for rpoS mRNA containing the leader region and 75 nt of coding sequence in the absence or presence of Hfq. This together with other results suggest OxyS represses rpoS mRNA translation by sequestering Hfq rather than binding to rpoS mRNA., (2011 Elsevier B.V. All rights reserved.)

The Escherichia coli RNA binding protein Hfq is involved in many aspects of post-transcriptional gene expression. Tight binding of Hfq to polyadenylate sequences at the 3' end of mRNAs influences exonucleolytic degradation, while Hfq binding to small noncoding RNAs (sRNA) and their targeted mRNAs facilitate their hybridization which in turn effects translation. Hfq binding to an A-rich tract in the 5' leader region of the rpoS mRNA and to the sRNA DsrA have been shown to be important for DsrA enhanced translation initiation of this mRNA. The complexes of Hfq-A(18) and Hfq-DsrA provide models for understanding how Hfq interacts with these two RNA sequence/structure motifs. Different methods have reported different values for the stoichiometry of Hfq-A(18) and Hfq-DsrA. In this work, mass spectrometry and analytical ultracentrifugation provide direct evidence that the strong binding mode of the Hfq hexamer (Hfq(6)) for A(18) and domain II of DsrA (DsrA(DII)) involve 1:1 complexes. This stoichiometry was also supported by fluorescence anisotropy and a competition gel mobility shift experiment using wild-type and truncated Hfq. More limited studies of Hfq binding to DsrA as well as to the sRNAs RprA, OxyS, and an 18-nt segment of OxyS were also consistent with 1:1 stoichiometry. Mass spectrometry of cross-linked samples of Hfq(6), A(18), and DsrA(DII) exhibit intensity corresponding to a ternary 1:1:1 complex; however, the small intensity of this peak and fluorescence anisotropy experiments did not provide evidence that this ternary complex is stable in solution.

The RNA-binding protein Hfq has been studied extensively for its function as a modulator of gene expression at the post-transcriptional level. While most Hfq studies have focused on the protein's interaction with sRNAs and mRNAs, Hfq binding to DNA has been observed but is less explored. During the isolation of Hfq from Escherichiacoli, we found genomic DNA fragments associated with the protein after multiple steps of purification. Sequences of 41 amplified segments from the DNA fragments associated with Hfq were determined. A large fraction of the DNA segments were predicted to have significant helical axis curvature and were from genes associated with membrane proteins, characteristics unexpected for non-specific binding. Analysis by analytical ultracentrifugation indicated that rA(18) binding to Hfq disrupts Hfq-DNA interactions. The latter observation suggests Hfq binding to DNA involves its distal surface. This was supported by a gel mobility shift assay that showed single amino acid mutations on the distal surface of Hfq inhibited Hfq binding to duplex DNA, while six of seven mutations on the proximal surface and outer circumference of the hexamer did not prevent Hfq binding. Two mutated Hfq which have portions of their C-terminal domain removed also failed to bind to DNA. The apparent K(d) for binding wild type Hfq to several duplex DNA was estimated from a gel mobility shift assay to be ~400nM., (Copyright © 2010 Elsevier B.V. All rights reserved.)

The conversion of a nucleic acid from single strands to double strands is thought to involve slow nucleation followed by fast double-strand propagation. Here, for RNA double-strand propagation, we propose an atomic resolution reaction mechanism. This mechanism, called the stack-ratchet, is based on data-mining of three-dimensional structures and on available thermodynamic information. The stack-ratchet mechanism extends and adds detail to the classic zipper model proposed by Porschke (Porschke, D. Biophysical Chemistry 1974, 2, pp. 97-101). Porschke's zipper model describes the addition of a base pair to a nucleated helix in terms of a single type of elementary reaction; a concerted process in which the two bases, one from each strand, participate in the transition state. In the stack-ratchet mechanism proposed here a net base-pairing step consists of two elementary reactions. Motions of only one strand are required to achieve a given transition state. One elementary reaction preorganizes and stacks the 3' single-strand, driven by base--base stacking interactions. A second elementary reaction stacks the 5' strand and pairs it with the preorganized 3' strand. In the stack-ratchet mechanism, a variable length 3' stack leads the single-strand/double-strand junction. The stack-ratchet mechanism is not a two-state process. A base can be (i) unstacked and unpaired, (ii) stacked and paired, or (ii) stacked and unpaired (only on the 3' strand). The data suggests that helices of DNA and of RNA do not propagate by similar mechanisms.

It is generally accepted that DNA conserves its secondary structure after a freeze-thaw cycle. A negligible amount of degradation occurs after this procedure. Degradation becomes appreciable only after multiple cycles of freezing and thawing. In this study, we have found that a single freeze-thaw cycle in alkaline medium (pH>or=10.8) gives rise to denaturation of calf thymus DNA, although the melting temperature of intact DNA in the solution used for the freeze-thaw experiments is higher than 60 degrees C. The degree of denaturation is almost independent of the regime of freezing. The melting curve obtained after DNA is frozen at -2 degrees C and then thawed is almost the same as after a freezing carried out in liquid nitrogen (-196 degrees C). However, incubation in the same solution at 0 degrees C for 24 hours without freezing does not give rise to any denaturation. The degree of denaturation caused by freezing increases with pH (if pH>or=10.8) and decreases with Na2CO3 concentration at fixed pH and [Na+], although Na2CO3 decreases the melting temperature of intact DNA. A preliminary treatment of DNA with cisplatin or transplatin (0.01 Pt atoms per nucleotide) gives rise to a full recovery of the DNA secondary structure after freezing and thawing similar to what occurs after heating DNA to 100 degrees C and quick cooling. Possible mechanisms that may cause DNA denaturation during a freeze-thaw cycle in alkaline medium are discussed.

The rpoS mRNA encodes a stress response transcription factor in Escherichia coli. It is one of a growing number of mRNAs found to be regulated by small RNAs (sRNA). Translation initiation of rpoS mRNA is enhanced by two sRNAs, DsrA and RprA, that pair to the same site near the rpoS start codon in the presence of the Hfq protein. In this work, we examine the interaction of E. coli Hfq with RprA and two portions of the rpoS mRNA leader region. One rpoS RNA, rpoS-L, contained the entire 565-nucleotide leader region, while the other, rpoS-S, contained the 199-nucleotide sequence surrounding the start codon. An RNase H assay indicated both rpoS RNAs have similar secondary structures in the translation initiation region. Hfq formed two complexes with RprA in a gel mobility assay with binding parameters similar to values previously determined for DsrA. Unlike DsrA, Hfq binding to RprA was inhibited by poly(A) and influenced by Hfq mutations on both the distal and proximal surfaces. Hfq increased the level of RprA binding to both rpoS RNAs but showed a much larger enhancement when rpoS-L, the entire leader region, was examined. The lower affinity of RprA for rpoS-L versus rpoS-S in the absence of Hfq suggests that Hfq overcomes an inhibitory structure within rpoS-L in stimulating RprA binding. Similar results were obtained with DsrA. The results indicate that the full upstream leader sequence of rpoS mRNA influences Hfq-facilitated annealing of RprA and DsrA and is likely to be involved in its regulation.

MicroRNAs (miRNAs) are known to regulate plant development, but have not been studied in gymnosperm seed tissues. The presence and characteristics of several miRNAs were examined in zygotic embryos (ZEs) and female gametophytes (FGs) of Pinus taeda (loblolly pine). Evidence for miRNAs was obtained using northern analyses and quantitative reverse transcription polymerase chain reaction (qRT-PCR) mediated with poly(A) polymerase. Partial sequences of two miRNAs were verified. Three regions of putative mRNA targets were analyzed by qRT-PCR to monitor the occurrence of stage-dependent miRNA-mediated cleavage. Five miRNAs were identified in ZEs and FGs along with partial sequences of Pta-miR166 and Pta-miR167. Both miRNAs showed differing degrees of tissue-specific and stage-specific modulation. Analysis of HB15L mRNA (a potential Pta-miR166 target) suggested miRNA-guided cleavage in ZEs and FGs. Analysis of ARF8L mRNA (a potential Pta-miR167 target) implied cleavage in ZEs but not in FGs. Argonaute9-like mRNA (ptAGO9L) showed stage-specific modulation of expression in ZEs that appeared to be inverted in the corresponding FGs. MicroRNAs and argonaute genes varied spatiotemporally during seed development. The peak levels of Pta-miR166 in FGs and ptAGO9L in embryos occurred at stage 9.1, a critical transition point during embryo development and a point where somatic embryo maturation often stops. MicroRNAs identified in FG tissue may play a role in embryogenesis.

Previous studies employing a 79-nucleotide (nt) RNA indicated that this RNA could form two bands in a native polyacrylamide gel while one band was observed in a denaturing gel. This report describes an investigation on the nature of the two corresponding structures and the segment responsible for forming the slower mobility band. Sedimentation equilibrium of the 79-nt RNA was consistent with the two gel bands corresponding to monomer and dimer forms. The portion of the RNA required for dimer formation was explored using a secondary structure prediction algorithm of two 79-nt RNAs linked in a head-to-tail fashion. The predicted structure suggested that the first 21-nt at the 5' end of each RNA formed a self-complementary duplex. A ribonuclease H assay carried out with RNA prepared as monomer (M), or a mixture of monomer and dimer (M/D), gave results consistent with the predicted M and D structures. Gel mobility experiments on 5' and 3' segments of the 79-nt RNA also indicated that dimer formation was due to the 21-nt 5' end. Studies on the 21-nt RNA molecule and sequence variants showed that this sequence can form a hairpin and a dimer complex. Unexpectedly, the hairpin to dimer conversion was shown to occur at high efficiency in frozen solution, although little or no conversion was observed above 0 degrees C. The results indicate that a freezing environment can promote formation of intermolecular RNA complexes from stable RNA hairpins, supporting the notion that this environment could have played a role in the evolution of RNA complexity.

Hfq is a RNA-binding protein in Escherichia coli that plays an essential role in post-transcriptional regulation of mRNAs by facilitating pairing of noncoding RNAs (ncRNAs) to mRNA target sites. Recent work has provided evidence that E. coli Hfq has two distinct RNA-binding surfaces. In this study, a comparative sequence-structure analysis of hfq genes in bacterial genomes was employed to identify conserved residues that may be involved in binding RNA. A covariance of residue properties at neighboring positions 12 and 39 and conserved surface residues with high propensities at binding sites of RNA-binding proteins suggested several sites for Hfq-RNA interactions. On the basis of these predictions, eight mutant Hfq proteins were produced and their interactions were examined with the 38 nucleotide (nt) domain II of DsrA ncRNA (DsrA(DII)) and A(18) by a gel-mobility shift assay, fluorescence anisotropy, and fluorescence quenching. Mutations on the proximal surface of Hfq had a small affect on Hfq binding to A(18) (

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.

The thermodynamics of nine hairpin DNAs were evaluated using UV-monitored melting curves and differential scanning calorimetry (DSC). Each DNA has the same five-base loop and a stem with 8-10 base pairs. Five of the DNAs have a tandem mismatch in the stem, while four have all base pairs. The tandem mismatches examined (ga/ga, aa/gc, ca/gc, ta/ac, and tc/tc) spanned the range of stability observed for this motif in a previous study of 28 tandem mismatches. UV-monitored melting curves were obtained in 1.0 M Na(+), 0.1 M Na(+), and 0.1 M Na(+) with 5 mM Mg(2+). DSC studies were conducted in 0.1 M Na(+). Transition T(m) values were unchanged over a 50-fold range of strand concentration. Model-independent enthalpy changes (DeltaH degrees ) evaluated by DSC were in good agreement (+/-8%) with enthalpy values determined by van't Hoff analyses of the melting curves in 0.1 M Na(+). The average heat capacity change (DeltaC(p)) associated with the hairpin to single strands transitions was estimated from plots of DeltaH degrees and DeltaS degrees with T(m) and ln T(m), respectively, and from profiles of DSC curves. The average DeltaC(p) values (113 +/- 9 and 42 +/- 27 cal x K(-1) x mol(-1) of bp), were in the range of values reported in previous studies. Consideration of DeltaC(p) produced large changes in DeltaH degrees and DeltaS degrees extrapolated from the transition region to 37 degrees C and smaller but significant changes to free energies. The loop free energy of the five tandem mismatches at 37 degrees C varied over a range of approximately 4 kcal x mol(-1) for each solvent.

RNase H degradation of two 15 nt RNA target sites was examined in the presence of hairpin DNAs with a 5 nt loop and a 10 bp stem or single-stranded 15 nt DNAs. One target site was a segment of a 79 nt RNA, and the other was part of a 53 nt RNA. Secondary structure predictions indicate that the 53 nt RNA target site is entirely single stranded, while a portion of the 79 nt RNA target site forms an intramolecular duplex. Less RNase H and DNA were needed to cleave the 53 nt RNA target site than the less accessible 79 nt RNA site. The hairpin DNAs had their 5 nt loop and 3' side of the stem fully complementary to the target sites or had sequence changes that produced one to nine mismatched pairs. T(m) values ranged from 57 to 80 degrees C. The stability of the hairpin DNAs relative to the stability of their corresponding RNA-DNA hybrids influenced the extent of RNase H degradation at 37 degrees C. Under the assay conditions employed, the amount of degradation directed by the hairpin DNAs was correlated with their predicted DeltaG(o) (37) of binding to the RNA targets. A DNA hairpin with one mismatch to the target site of the 79 nt RNA did not induce degradation under conditions where fully complementary DNA hairpins produced 50-80% degradation. The in vitro results indicate that DNA hairpins can enhance the stringency of RNase H targeted degradation of the RNA sites.

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this "cross-linking effect" (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the "cross-linking effect" (18 - 20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20- 40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.

Hfq, a bacterial RNA-binding protein, was recently shown to contain the Sm1 motif, a characteristic of Sm and LSm proteins that function in RNA processing events in archaea and eukaryotes. In this report, comparative structural modeling was used to predict a three-dimensional structure of the Hfq core sequence. The predicted structure aligns with most major features of the Methanobacterium thermoautotrophicum LSm protein structure. Conserved residues in Hfq are positioned at the same structural locations responsible for subunit assembly and RNA interaction in Sm proteins. A highly conserved portion of Hfq assumes a structural fold similar to the Sm2 motif of Sm proteins. The evolution of the Hfq protein was explored by conducting a BLAST search of microbial genomes followed by phylogenetic analysis. Approximately half of the 140 complete or nearly complete genomes examined contain at least one gene coding for Hfq. The presence or absence of Hfq closely followed major bacterial clades. It is absent from high-level clades and present in the ancient Thermotogales-Aquificales clade and all proteobacteria except for those that have undergone major reduction in genome size. Residues at three positions in Hfq form signatures for the beta/gamma proteobacteria, alpha proteobacteria and low GC Gram-positive bacteria groups.

The Escherichia coli strain FC40 has frequently been employed to investigate the mechanism of adaptive mutations. The strain cannot utilize lactose due to a +1 frameshift mutation that reduces beta-galactosidase to about 1% of normal levels. Cells undergo a high rate of mutation from Lac- to Lac+ when cells are grown with lactose as the sole energy source. Almost all Lac+ colonies arising 3-6 days after plating result from a base pair deletion in runs of iterated base pairs within a 130-bp target region. In this study we characterized Lac+ colonies arising 3-10 days after plating. Temperature gradient gel electrophoresis (TGGE) was used to detect mutations in the target region as a function of the day a colony appears. TGGE results confirmed the occurrence of mutations within the target region in 36 of 37 FC40 Lac+ colonies arising on days 3-7. However, mutations in this region were not detected in 23 of 37 Lac+ colonies arising from days 8-10. Sequencing data verified the TGGE results. Half of the Lac+ mutants arising on days 8-10 with no base pair change in the target region were unstable and exhibited a Lac- phenotype after successive growth cycles in rich medium. The results suggest that amplification of the lac operon region is a common factor in late arising colonies, and that different characteristics distinguish early and late arising Lac+ colonies.

Forty-eight RNA duplexes were constructed that contained all common single base bulges at six different locations. The stabilities of the RNAs were determined by temperature gradient gel electrophoresis (TGGE). The relative stability of a single base bulge was dependent on both base identity and the nearest neighbor context. The single base bulges were placed into two categories. A bulged base with no identical neighboring base was defined as a Group I base bulge. Group II-bulged bases had at least one neighboring base identical to it. Group II bulges were generally more stable than Group I bulges in the same nearest neighbor environments. This indicates that position degeneracy of an unpaired base enhances stability. Differences in the mobility transition temperatures between the RNA fragments with bulges and the completely base-paired reference RNAs were related to free energy differences. Simple models for estimating the free energy contribution of single base bulges were evaluated from the free energy difference data. The contribution of a Group I bulge 5'-(XNZ)-3'.5'-(Z'-X')-3' where N is the unpaired base and X.X' and Z.Z' the neighboring base pairs, could be well-represented (+/-0.34 kcal/mol) by the equation, DeltaG((X)(N)()(Z))(.)((Z)(')(-)(X)(')()) = 3.11 + 0. 40DeltaG(s)()((XZ))(.)((Z)(')(X)(')()). DeltaG(s)()((XZ))(. )((Z)(')(X)(')()) is the stacking energy of the closing base pair doublet. By adding a constant term, delta = -0.3 kcal/mol, to the right side of the above equation, free energies of Group II bulges could also be predicted with the same accuracy. The term delta represents the stabilizing effect due to position degeneracy. A similar equation/model was applied to previous data from 32 DNA fragments with single base bulges. It predicted the free energy differences with a similar standard deviation.

The ability of single-stranded DNA oligomers to form adjacent triplex and duplex domains with two DNA structural motifs was examined. Helix-coil transition curves and a gel mobility shift assay were used to characterize the interaction of single-stranded oligomers 12-20 nt in length with a DNA hairpin and with a DNA duplex that has a dangling end. The 12 nt on the 5'-ends of the oligomers could form a triplex structure with the 12 bp stem of the hairpin or the duplex portion of the DNA with a dangling end. The 3'-ends of the 17-20 nt strands could form Watson-Crick pairs to the five base loop of the hairpin or the dangling end of the duplex. Complexes of the hairpin DNA with the single-stranded oligomers showed two step transitions consistent with unwinding of the triplex strand followed by hairpin denaturation. Melting curve and gel competition results indicated that the complex of the hairpin and the 12 nt oligomer was more stable than the complexes involving the extended single strands. In contrast, results indicated that the extended single-stranded oligomers formed Watson-Crick base pairs with the dangling end of the duplex DNA and enhanced the stability of the adjacent triplex region.

Base pair stacking free energy parameters in a low ionic strength solvent were determined from an analysis of DNA fragments using temperature gradient gel electrophoresis (TGGE). Transition midpoint temperatures (Tu) were determined for the first melting domain (52 +/- 4 bp) of 16, 339 bp DNAs that differed from each other by single base pair substitutions. The data were combined with previously obtained Tu data from 17 similar DNAs that had single base pair changes at different sites [Ke, S. H., and Wartell, R. M. (1995) Biochemistry 34, 4593-4599]. The Tu values were used to evaluate free energy differences (deltaDeltaG) between 31 pairs of DNAs. Linear equations relating the deltaDeltaG values to changes in base pair stacking were analyzed by singular value decomposition (SVD) to determine the 10 nearest neighbor free energy parameters. The order of stability of the parameters, TA < AT < AA < AG < GT approximately TC approximately TG < CC < GC approximately CG, was essentially the same as the hierarchy determined in 1 M Na+ [Allawi, H. T., and SantaLucia, J., Jr. (1997) Biochemistry 36, 10581-10594]. The experimental free energy differences were in good agreement with predictions made using nearest-neighbor parameters determined from several previous studies conducted in medium or high salt concentrations. Conversely the parameters determined in the current study produced good predictions of free energy differences previously determined from 59 DNA oligomers in 1 M Na+. The results indicate that differences between base pair stacking energies are conserved across a wide range of ionic conditions, and in both oligomer and polymer DNA contexts.

Temperature gradient gel electrophoresis (TGGE) and related methods can separate DNA fragments that differ by a single base pair or defect. This article describes the basic features of TGGE, and reviews the theoretical model of DNA unwinding and its ability to predict DNA mobility in a temperature gradient gel. Recent applications of TGGE and related methods that were directed at detecting point mutations, and evaluating the effects of single site defects are also reported.

Cleavage of a RNA target site by RNase H1 from Escherichia coli was examined in the presence of complementary DNA sequences in the form of single-stranded, duplex, and hairpin structures. The target site was a 15 nt sequence in the middle of a 79 nt RNA transcript. DNA molecules employed included seven single-stranded oligodeoxynucleotides 10 or 15 nt long, and five hairpin DNAs each with a 10 bp stem and 5 nt loop. The loop and 3' side of the stem of two of the hairpin DNAs were fully complementary to the target site, while the other hairpin DNAs had sequence changes. A 10 bp duplex DNA with one strand complementary to the target site was also employed. A gel electrophoresis mobility shift assay examined hybrid formation between the RNA and the single-stranded 15 nt DNA and two hairpin DNAs that contained 15 complementary bases. RNA titration of the 32P-labeled single-stranded DNA produced a shifted band indicative of RNA/DNA complex formation. No RNA/DNA complex was detected when the more stable (Tm = 71 degrees C) hairpin DNA was combined with excess RNA. The less stable hairpin DNA (Tm = 62 degrees C) showed a small amount ( approximately 8%) of hybrid formation. Thermodynamic analysis of RNA binding to the DNAs was in qualitative agreement with the results. Although no RNA/DNA hybrid was expected from thermodynamic calculations, a RNase H assay at 25 degrees C showed that hairpin or duplex DNAs with a 10 nt complementary sequence catalyzed RNA degradation. A complementary loop sequence in the hairpin DNA was not required. Cleavage of the RNA did not occur with hairpin DNAs containing three or four noncomplementary bases in the stem. The results show that RNase H can promote the formation and cleavage of a RNA/DNA hybrid between an RNA site and a base paired strand of a stable hairpin or duplex DNA at temperatures below their Tm.

The thermal stability of RNA duplexes differing by a single base pair (bp) substitution or mismatch were investigated by temperature gradient gel electrophoresis (TGGE). All base pair substitutions and mismatches were examined at six sites, and limited changes were investigated at three other sites. DNA templates for in vitro transcription were generated by the polymerase chain reaction (PCR). Transcribed forward and reverse single stranded RNAs were annealed to form 345 bp dupex RNA. Solution melting curves of selected RNAs were in good agreement with the predicted three step transitions. Parallel TGGE was used to determine the relative stabilities of the RNAs, and perpendicular TGGE was employed to obtain mobility transitions and midpoint transition temperatures (Tmu) of the RNAs' first melting domain. The gel solvent included formamide and urea. The Tmu values of the first melting domain were influenced by the identity of the base pair substitution or mismatch as well as by the site's neighboring base pairs. The difference in the transition temperatures (deltaTmu) between pairs of RNA ranged from 0 to 5 degrees C. deltaTmu values were used to determine free energy differences (deltaDeltaG). For RNA pairs distinguished by a base pair substitution, the deltaDeltaG values were closely correlated with free energy differences calculated from stacking free energies determined from melting studies in 1 M Na+ [Serra, M. J., and Turner, D. H. (1995) Methods Enzymol. 259, 242-261.] An algorithm was developed using the free energies of terminal mismatches [Serra, M. J., and Turner, D. H. (1995) Methods Enzymol. 259, 242-261] that provided very good agreement with experimental free energies for the single internal mismatches.

Temperature-Gradient Gel Electrophoresis (TGGE) was employed to determine the thermal stabilities of 28 DNA fragments, 373 bp long, with two adjacent mismatched base pairs, and eight DNAs with Watson-Crick base pairs at the same positions. Heteroduplex DNAs containing two adjacent mismatches were formed by melting and reannealing pairs of homologous 373 bp DNA fragments differing by two adjacent base pairs. Product DNAs were separated based on their thermal stability by parallel and perpendicular TGGE. The polyacrylamide gel contained 3.36 M urea and 19.2 % formamide to lower the DNA melting temperatures. The order of stability was determined in the sequence context d(CXYG).d(CY'X'G) where X.X' and Y.Y" represent the mismatched or Watson-Crick base pairs. The identity of the mismatched bases and their stacking interactions influence DNA stability. Mobility transition melting temperatures (T u) of the DNAs with adjacent mismatches were 1.0-3.6 degrees C (+/-0.2 degree C) lower than the homoduplex DNA with the d(CCAG).d(CTGG) sequence. Two adjacent G.A pairs, d(CGAG).d(CGAG), created a more stable DNA than DNAs with Watson-Crick A.T pairs at the same sites. The d(GA).d(GA) sequence is estimated to be 0.4 (+/-30%) kcal/mol more stable in free energy than d(AA).d(TT) base pairs. This result confirms the unusual stability of the d(GA).d(GA) sequence previously observed in DNA oligomers. All other DNAs with adjacent mismatched base pairs were less stable than Watson-Crick homoduplex DNAs. Their relative stabilities followed an order expected from previous results on single mismatches. Two homoduplex DNAs with identical nearest neighbor sequences but different next-nearest neighbor sequences had a small but reproducible difference in T u value. This result indicates that sequence dependent next neighbor stacking interactions influence DNA stability.

Temperature-gradient gel electrophoresis (TGGE) was used to determine the relative thermal stabilities of 32 DNA fragments that differ by a single unpaired base (base bulge) and 17 DNAs differing by a base pair. Homologus 373 and 372 bp DNA fragments differing by a single base pair substitution or deletion were employed. Heteroduplexes containing a single base bulge were formed by melting and reannealing pairs of 372 and 373 bp DNAs. Product DNAs were separated on the basis of their thermal stability by parallel and perpendicular TGGE. The order of stability was determined for all single unpaired bases in four different nearest neighbor environments: (GXT).(AYC), (GXG).(CYC), (CXA).(TYG), and (TXT).(AYA) with X = A, T, G, or C, and Y = no base, or visa versa. DNA fragments containing a base bulge were destabilized by 2-3.6 degrees C (+/- 0.2 degrees C) with respect to homologous DNAs with complete Watson-Crick base pairing. Both the identity of the unpaired base and the sequence of the flanking base pairs influenced the degree of destabilization. The range of temperature shift correspond to estimated unfavorable free energies from 2.5 to 4.6 kcal/mol. Purine base bulges were generally not as destabilizing as pyrimidine base bulges. An unpaired base which was identical to one of its adjacent bases generally caused less destabilization than an unpaired base with an identity differing from its nearest neighbors. This implies that positional degeneracy of an unpaired base within a run of two or more identical bases is an important factor effecting stability.(ABSTRACT TRUNCATED AT 250 WORDS)

A computer program was developed to automate the selection of DNA fragments for detecting mutations within a long DNA sequence by denaturing gel electrophoresis methods. The program, MELTSCAN, scans through a user specified DNA sequence calculating the melting behavior of overlapping DNA fragments covering the sequence. Melting characteristics of the fragments are analyzed to determine the best fragment for detecting mutations at each base pair position in the sequence. The calculation also determines the optimal fragment for detecting mutations within a user specified mutational hot spot region. The program is built around the statistical mechanical model of the DNA melting transition. The optimal fragment for a given position is selected using the criteria that its melting curve has at least two steps, the base pair position is in the fragment's lowest melting domain, and the melting domain has the smallest number of base pairs among fragments that meet the first two criteria. The program predicted fragments for detecting mutations in the cDNA and genomic DNA of the human p53 gene.

DNA sequences upstream of the rrnB P1 core promoter (-10, -35 region) increase transcription more than 300-fold in vivo and in vitro. This stimulation results from a cis-acting DNA sequence, the UP element, which interacts directly with the alpha subunit of RNA polymerase, increasing transcription about 30-fold, and from a positively acting transcription factor, FIS, which increases expression another 10-fold. A DNA region exhibiting a high degree of intrinsic curvature has been observed upstream of the rrnB P1 core promoter and has thus been often cited as an example of the effect of bending on transcription. However, the precise position of the curvature has not been determined. We address here whether this bend is in fact related to activation of rRNA transcription. Electrophoretic analyses were used to localize the major bend in the rrnB P1 upstream region to position approximately -100 with respect to the transcription initiation site. Since most of the effect of upstream sequences on transcription results from DNA between the -35 hexamer and position -88, i.e. downstream of the bend center, these studies indicate that the curvature leading to the unusual electrophoretic behavior of the upstream region does not play a major role in activation of rRNA transcription. Minor deviations from normal electrophoretic behavior were associated with the region just upstream of the -35 hexamer and could conceivably influence interactions between the UP element and the alpha subunit of RNA polymerase.

Temperature-gradient gel electrophoresis (TGGE) was employed to determine the thermal stabilities of 48 DNA fragments that differ by single base pair mismatches. The approach provides a rapid way for studying how specific base mismatches effect the stability of a long DNA fragment. Homologous 373 bp DNA fragments differing by single base pair substitutions in their first melting domain were employed. Heteroduplexes were formed by melting and reannealing pairs of DNAs, one of which was 32P-labeled on its 5'-end. Product DNAs were separated based on their thermal stability by parallel and perpendicular temperature-gradient gel electrophoresis. The order of stability was determined for all common base pairs and mismatched bases in four different nearest neighbor environments; d(GXT).d(AYC), d(GXG).d(CYC), d(CXA).d(TYG), and d(TXT).d(AYA) with X,Y = A, T, C, or G. DNA fragments containing a single mismatch were destabilized by 1 to 5 degrees C with respect to homologous DNAs with complete Watson-Crick base pairing. Both the bases at the mismatch site and neighboring stacking interactions influence the destabilization caused by a mismatch. G.T, G.G and G.A mismatches were always among the most stable mismatches for all nearest neighbor environments examined. Purine.purine mismatches were generally more stable than pyrimidine.pyrimidine mispairs. Our results are in very good agreement with data where available from solution studies of short DNA oligomers.

Temperature gradient gel electrophoresis (TGGE) and related methods are widely employed to detect mutations in DNA fragments. DNA melting map calculations and GC clamps have been used to enhance the detection of mutations. While generally successful, these methods have not always revealed base changes within a DNA fragment. Previous work suggested that mutations are detected if they are in a DNA's first melting domain, and the melting domain is well separated from final strand dissociation. Two criteria from the DNA melting theory were established to determine when both of these conditions are met. The criteria involve calculating the derivative melting curve as well as the melting map of a DNA sequence. The approach was applied to the cDNA sequence of the human p53 gene. Mutations in the p53 gene are common in human cancers and are generally located in four 'hot spot' regions. Calculations indicated that three DNA fragments are needed to detect base substitutions in the four hot-spot regions. Predicted melting behavior was experimentally tested with eight single base substitutions distributed among the four hot-spot regions. All mutations tested behaved as predicted and were detected by vertical TGGE. Heteroduplex DNAs formed by melting and reannealing various ratios of wild type and mutant DNA fragments were also examined. Results indicated that point mutations can be detected by ethidium bromide staining from samples containing 10% mutant and 90% wild-type sequences.

The secondary structures of the cAMP receptor protein (CRP), a complex of CRP and cAMP, and a cAMP-independent receptor protein mutant (CRP*141 gln) were examined by using Raman spectroscopy. Spectra were obtained from CRP and CRP*141 gln dissolved in 0.3 M NaCl and 30 mM sodium phosphate at protein concentrations of 30-40 mg/mL. CRP and CRP.cAMP1 were compared at lower protein concentrations (10-12 mg/mL) in a solvent of 0.35 M NaCl and 20 mM sodium phosphate. Raman analysis indicates that CRP structural changes induced by one bound cAMP or by the Gly to Gln mutation at residue 141 are small. Spectra of the three CRP samples are essentially identical from 400 to 1900 cm-1. This result differs from the Raman spectroscopy study of CRP and CRP.cAMP2 cocrystals [DeGrazia et al. (1990) Biochemistry 29, 3557]. The latter work showed spectral differences between CRP and CRP.cAMP2 consistent with alterations in the protein conformation. These studies indicate that CRP and CRP.cAMP1 in solution are similar in structure and differ from CRP.cAMP2 cocrystals. Protease digestion and a DNA binding assay were also employed to characterize the wild-type and mutant proteins. CRP*141 gln exhibited the same conformational characteristics of previously reported cAMP-independent mutant proteins. It was sensitive to proteolytic attack in the absence of cAMP, or upon addition of cGMP. In the absence of cAMP, both wild-type and mutant CRPs bound noncooperatively to a 62 bp lac promoter DNA. The equilibrium constants were approximately 10(6) M-1 in 0.1 M Na+. CRP*141 gln had a 2-4-fold higher affinity for the 62 bp DNA than CRP.(ABSTRACT TRUNCATED AT 250 WORDS)

A vertical gel electrophoresis apparatus is described which can distinguish DNA fragments differing by single base pair substitutions. The system employs a homogenous polyacrylamide gel containing urea-formamide and a temperature gradient which runs either perpendicular or parallel to the direction of electrophoresis. The temperature-gradient system simplifies several features of the denaturant-gradient system (1) and is relatively inexpensive to construct. Eight homologous 373 bp DNAs differing by one, two, or nine base pair substitutions were examined. DNA electrophoretic mobility changed abruptly with the temperature induced unwinding of DNA domains. GC to AT substitutions at different locations within the first melting domain, as well as an AT to TA transversion were separated with temperature gradients parallel to the electrophoretic direction. The relative stabilities of the DNAs observed in the gels were compared to predictions of DNA melting theory. General agreement was observed however complete correspondence was not obtained.

Tm curves, CD spectra, and kinetics results of the self-complementary DNA dodecamers d(A6T6), d(A3T3A3T3), d(A2T2A2T2A2T2), d(ATATATATATAT), and d(T6A6) demonstrate that the thermal transitions of these oligomers at low salt concentration involve a hairpin intermediate. At high salt concentrations (greater than 0.1 M Na+) only a duplex to denatured-strand transition appears to occur. The temperature and salt-concentration regions of the transitions are very sequence dependent. Alternating-type AT sequences have a lower duplex stability and a greater tendency to form hairpins than sequences containing more nonalternating AT base pairs. Of the two nonalternating sequences, d(T6A6) is significantly less stable than d(A6T6). Both oligomers have CD curves that are very similar to the unusual CD spectrum of poly(dA).poly(dT). The Raman spectra of these two oligomers are also quite similar, but at low temperature, small intensity differences in two backbone modes and three nucleoside vibrations are obtained. The hairpin to duplex transition for the AT dodecamers was examined by salt-jump kinetics measurements. The transition is faster than transitions for palindromic-sequence oligomers containing terminal GC base pairs. Stopped-flow kinetics studies indicate that the transition is second order and has a relatively low activation energy. The reaction rate increases with increasing ionic strength. These results are consistent with a three-step mechanism for the hairpin to duplex reaction: (i) fraying of the hairpin oligomers' terminal base pairs, (ii) a rate-determining bimolecular step involving formation of a cruciform-type intermediate from two hairpin oligomers with open terminal base pairs, and (iii) base-pair migration and formation in the intermediate to give the duplex.

Raman spectroscopy was employed to examine the secondary structure of the cAMP receptor protein (CRP). Spectra were obtained over the range 400-1900 cm-1 from solutions of CRP and from CRP-cAMP cocrystals. The spectra of CRP dissolved in 30 mM sodium phosphate and 0.15 M NaCl buffered at either pH 6 or pH 8 or dissolved in 0.15-0.2 M NaCl at protein concentrations of 5, 15, and 30 mg/mL were examined. Estimates of the secondary structure distribution were made by analyzing the amide I region of the spectra (1630-1700 cm-1). CRP secondary structure distributions were essentially the same in either pH and at all protein concentrations examined. The amide I analyses indicated a structural distribution of 44% alpha-helix, 28% beta-strand, 18% turn, and 10% undefined for CRP in solution. Raman spectra of CRP-cAMP cocrystals differed from the spectra of CRP in solution. Some differences were assigned to interfering background bands, whereas other spectral differences were attributed to changes in CRP structure. Differences in the amide III region and in the intensity at 935 cm-1 were consistent with alterations in secondary structure. Analysis of the amide I region of the CRP-cAMP cocrystal spectrum indicated a secondary structure distribution of 37% alpha-helix, 33% beta-strand, 17% turn, and 12% undefined. This result is in agreement with a published secondary structure distribution derived from X-ray analysis of CRP-cAMP cocrystals (37% alpha-helix and 36% beta-strand).(ABSTRACT TRUNCATED AT 250 WORDS)

Gal repressor dimer binds to two gal operator sites, OE and OI, which are 16 bp long similar sequences with hyphenated dyad symmetries (11,12). Repressor occupation hinders the reactivity of the N7 atoms in the major groups of guanines, located at positions 1, 3 and 8, and the rotational 1', 3' and 8' of the symmetries. We have shown that Gal repressor binding to OE or OI DNA fragments increases the circular dichroism (CD) spectral peak in the 270 to 300 nm range. The CD change is similar to that observed for Lac repressor binding to its operator site (14). It is consistent with a DNA conformational change during complex formation between Gal repressor and OE and OI DNA. The CD spectral change was not observed when the central 8,8' G-C base pairs in the DNA-protein complex were replaced by A-T base pairs, whereas substitution of the 1,1' G-C base pairs do show the accompanying increase in the spectra during repressor binding. The absence of CD change of the Gal repressor complex with DNA mutated at the 8,8' base pairs suggest that the central G-C base pairs are required for the repressor induced conformational change.

Theoretical calculations were made to determine the influence of side specific 'melting' and 'stabilizing' proteins on the thermal stability of nearby base pairs (bp). A DNA sequence 999bp. long containing the 123 bp. lactose operon control region in the center was examined. Melting curves of base pairs near the binding sites of the catabolite activator protein, CAP, the lactose repressor, and RNA polymerase were calculated in the absence and presence of each protein. The empirical loop entropy model of the helix-coil transition of DNA was employed. Calculations show that melting and stabilizing proteins alter the tm of base pairs 20 to 100 bp-away. The magnitude and range of the effect is strongly influenced by the base pair composition and sequence of the protein site and the immediately adjacent DNA regions.