Roger M. Wartell

Thermal denaturation of DNA molecules: A comparison of theory with experiment

Abstract Experimental and theoretical results on the helix-coil transition of DNA are reviewed. The theoretical model of the transition is described, and the influence of heterogeneous base pair stacking, and strand dissociation on the predicted melting transition is examined. New experimental transition data on seven DNAs, 154–587 base pairs (bp) long, are reported and compared with theoretical calculations. We review and evaluate previous studies on long DNAs (≥1000 bp) as well as previous and recent results on short DNAs. The comparison of theory with equilibrium melting curves of short DNAs indicates that base pair sequence has a relatively small influence on the stacking free energy. Excellent agreement is obtained between theory and equilibrium transitions of 14 out of 15 fragments 80–587 pb. The deviation between theory and experiment for a 516 bp DNA can be attributed to the formation of stem-loop structures. This may provide the explanation for inconsistent results observed with long DNAs. The effect of single base pair changes on DNA transitions is discussed. Current views on fluctuational opening of base pairs at temperatures below the transition are described.

Evolution of the ribosome at atomic resolution

Significance Ribosomes exist in every cell and are responsible for translation from mRNA to protein. The structure of the ribosomal common core is highly conserved in all living species, while the outer regions of the ribosome are variable. Ribosomal RNA of eukaryotes contains expansion segments accreted onto the surface of the core, which is nearly identical in structure to that in prokaryotic ribosomes. Comparing eukaryotic and prokaryotic ribosomes allows us to identify 3D insertion fingerprints of the expansion segments. Similar fingerprints allow us to analyze the common core and detect ancestral expansion segments within it. We construct a molecular model of ribosomal evolution starting from primordial biological systems near the dawn of life, culminating with relatively recent changes specific to metazoans. 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.

Effect of Hfq on RprA-rpoS mRNA pairing: Hfq-RNA binding and the influence of the 5' rpoS mRNA leader region.

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.

Evidence for stage‐specific modulation of specific microRNAs (miRNAs) and miRNA processing components in zygotic embryo and female gametophyte of loblolly pine (Pinus taeda)

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.

Detecting single base substitutions, mismatches and bulges in DNA by temperature gradient gel electrophoresis and related methods

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.

Secondary structure and domain architecture of the 23S and 5S rRNAs

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

RNA Folding and Catalysis Mediated by Iron (II)

Mg2+ 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 Fe2+ in the absence of free oxygen as a replacement for Mg2+ in RNA folding and catalysis. We describe both quantum mechanical calculations and experiments that suggest that the roles of Mg2+ in RNA folding and function can indeed be served by Fe2+. The results of quantum mechanical calculations show that the geometry of coordination of Fe2+ by RNA phosphates is similar to that of Mg2+. Chemical footprinting experiments suggest that the conformation of the Tetrahymena thermophila Group I intron P4–P6 domain RNA is conserved between complexes with Fe2+ or Mg2+. The catalytic activities of both the L1 ribozyme ligase, obtained previously by in vitro selection in the presence of Mg2+, and the hammerhead ribozyme are enhanced in the presence of Fe2+ compared to Mg2+. All chemical footprinting and ribozyme assays in the presence of Fe2+ 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 Fe2+. The combined biochemical and paleogeological data are consistent with a role for Fe2+ in an RNA World. RNA and Fe2+ could, in principle, support an array of RNA structures and catalytic functions more diverse than RNA with Mg2+ alone.

DNA bending induced by the catabolite activator protein allows ring formation of a 144 bp DNA.

The effect of the catabolite activator protein, CAP, on the ligation of a 144 bp DNA was examined. This DNA has EcoRI ends and contains the lac operon CAP site and promoter-operator region. At low DNA concentrations (nM) and 37 degrees C the presence of CAP and cAMP enables T4 ligase to convert the linear duplex to a covalently closed ring. Nuclease digestion and sedimentation equilibrium studies show that the ring is a monomer circle. Ring formation does not occur in the absence of either CAP or cAMP. The kinetics of ring closure, and the bimolecular joining of two fragments were measured. The presence of CAP decreased the rate of bimolecular joining of the EcoRI ends of linear DNAs. Thus the measured rates of ring closure are likely to be a lower limit for this process. Closure reactions carried out with ethidium bromide indicate that CAP induced bending rather than twisting is responsible for ring formation. The all or none nature of the closure reaction suggests that persistence length DNAs may be useful in a simple assay for protein induced DNA bending.

History of the ribosome and the origin of translation

Significance The ribosome, in analogy with a tree, contains a record of its history, spanning 4 billion years of life on earth. The information contained within ribosomes connects us to the prehistory of biology. Details of ribosomal RNA variation, observed by comparing three-dimensional structures of ribosomes across the tree of life, form the basis of our molecular-level model of the origins and evolution of the translational system. We infer many steps in the evolution of translation, mapping out acquisition of structure and function, revealing much about how modern biology originated from ancestral chemical systems. 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.

E. coli DNA associated with isolated Hfq interacts with Hfq's distal surface and C-terminal domain

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