Dynamics and stability of individual base pairs in two homologous RNA-DNA hybrids

Coarse-grained modeling of DNA-RNA hybrids

We introduce oxNA, a new model for the simulation of DNA-RNA hybrids that is based on two previously developed coarse-grained models-oxDNA and oxRNA. The model naturally reproduces the physical properties of hybrid duplexes, including their structure, persistence length, and force-extension characteristics. By parameterizing the DNA-RNA hydrogen bonding interaction, we fit the model's thermodynamic properties to experimental data using both average-sequence and sequence-dependent parameters. To demonstrate the model's applicability, we provide three examples of its use-calculating the free energy profiles of hybrid strand displacement reactions, studying the resolution of a short R-loop, and simulating RNA-scaffolded wireframe origami.

Shedding light on the base-pair opening dynamics of nucleic acids in living human cells

Base-pair opening is a fundamental property of nucleic acids that plays important roles in biological functions. However, studying the base-pair opening dynamics inside living cells has remained challenging. Here, to determine the base-pair opening kinetics inside living human cells, the exchange rate constant ([Formula: see text]) of the imino proton with the proton of solvent water involved in hairpin and G-quadruplex (GQ) structures is determined by the in-cell NMR technique. It is deduced on determination of [Formula: see text] values that at least some G-C base pairs of the hairpin structure and all G-G base-pairs of the GQ structure open more frequently in living human cells than in vitro. It is suggested that interactions with endogenous proteins could be responsible for the increase in frequency of base-pair opening. Our studies demonstrate a difference in dynamics of nucleic acids between in-cell and in vitro conditions.

Highly sensitive MALDI-MS measurement of active ricin: insight from more potential deoxynucleobase-hybrid oligonucleotide substrates

Herein, we report an improved MALDI-MS method for active ricin to contribute toward countermeasures against its real threat to the public. Compared with commonly used DNA or RNA substrates, the deoxynucleobase-hybrid oligonucleotide (RNA_dA, Rd) substrate containing functional Gd[combining low line]A[combining low line]GA loop was revealed as a substrate with more potential and used for the first time in ricin measurement via MALDI-MS. The Rd sequence greatly prompted ricin to exhibit its catalytic activity as rRNA N-glycosylase in ex vitro condition, which was supported by molecular docking simulation and enzymatic parameters depicted in MALDI-MS. Furthermore, we discovered that a highly pure matrix was the most crucial parameter for enhancing the sensitivity, which addressed the major obstacle encountered in the oligo(deoxy)nucleotide measurement, i.e., the interfering alkali metal ion-adducted signals in MALDI-MS. After the optimization of pH and enzymatic reaction buffer composition in this ex vitro condition, this method can provide a wide linearity of up to three orders of magnitude, i.e., 1-5000 ng mL-1, and a high sensitivity of 1 ng mL-1 without any enrichment. Denatured and active ricin could be distinctly differentiated, and the application to practical samples from one international exercise and a soft drink proved the feasibility of this new method. We believe this MALDI-MS method can contribute to the first response to ricin occurrence events in public safety and security, as well as pave a new way for a deep understanding of ricin and other type II ribosome inactivating proteins involved toxicology.

Region 4 of the RNA polymerase σ subunit counteracts pausing during initial transcription

All cellular genetic information is transcribed into RNA by multisubunit RNA polymerases (RNAPs). The basal transcription initiation factors of cellular RNAPs stimulate the initial RNA synthesis via poorly understood mechanisms. Here, we explored the mechanism employed by the bacterial factor σ in promoter-independent initial transcription. We found that the RNAP holoenzyme lacking the promoter-binding domain σ4 is ineffective in de novo transcription initiation and displays high propensity to pausing upon extension of RNAs 3 to 7 nucleotides in length. The nucleotide at the RNA 3' end determines the pause lifetime. The σ4 domain stabilizes short RNA:DNA hybrids and suppresses pausing by stimulating RNAP active-center translocation. The antipausing activity of σ4 is modulated by its interaction with the β subunit flap domain and by the σ remodeling factors AsiA and RbpA. Our results suggest that the presence of σ4 within the RNA exit channel compensates for the intrinsic instability of short RNA:DNA hybrids by increasing RNAP processivity, thus favoring productive transcription initiation. This "RNAP boosting" activity of the initiation factor is shaped by the thermodynamics of RNA:DNA interactions and thus, should be relevant for any factor-dependent RNAP.

GC-Content Dependence of Elastic and Overstretching Properties of DNA:RNA Hybrid Duplexes

DNA:RNA hybrid duplex plays important roles in various biological processes. Both its structural stability and interactions with proteins are highly sequence dependent. In this study, we utilize homebuilt optical tweezers to investigate how GC contents in the sequence influence the structural and mechanical properties of DNA:RNA hybrid by measuring its contour length, elasticities, and overstretching dynamics. Our results support that the DNA:RNA hybrid adopts a conformation between the A- and B-form helix, and the GC content does not affect its structural and elastic parameters obviously when varying from 40 to 60% before the overstretching transition of DNA:RNA hybrid occurs. In the overstretching transition, however, our study unravels significant heterogeneity and strong sequence dependence, suggesting the presence of a highly dynamic competition between the two processes, namely the S-form duplex formation (nonhysteretic) and the unpeeling (hysteretic). Analyzing the components left in DNA:RNA hybrid after the overstretching transition suggests that the RNA strand is more easily unpeeled than the DNA strand, whereas an increase in the GC content from 40 to 60% can significantly reduce unpeeling. Large hysteresis is observed between the stretching and relaxation processes, which is also quantitatively correlated with the percentage of unpeeling in the DNA:RNA duplex. Increasing in both the salt concentration and GC content can effectively reduce the hysteresis with the latter being more significant. Together, our study reveals that the mechanical properties of DNA:RNA hybrid duplexes are significantly different from double-stranded DNA and double-stranded RNA, and its overstretching behavior is highly sequence dependent. These results should be taken into account in the future studies on DNA:RNA-hybrid-related functional structures and motor proteins.

Structural Flexibility of DNA-RNA Hybrid Duplex: Stretching and Twist-Stretch Coupling

DNA-RNA hybrid (DRH) duplexes play essential roles during the replication of DNA and the reverse transcription of RNA viruses, and their flexibility is important for their biological functions. Recent experiments indicated that A-form RNA and B-form DNA have a strikingly different flexibility in stretching and twist-stretch coupling, and the structural flexibility of DRH duplex is of great interest, especially in stretching and twist-stretch coupling. In this work, we performed microsecond all-atom molecular dynamics simulations with new AMBER force fields to characterize the structural flexibility of DRH duplex in stretching and twist-stretch coupling. We have calculated all the helical parameters, stretch modulus, and twist-stretch coupling parameters for the DRH duplex. First, our analyses on structure suggest that the DRH duplex exhibits an intermediate conformation between A- and B-forms and closer to A-form, which can be attributed to the stronger rigidity of the RNA strand than the DNA strand. Second, our calculations show that the DRH duplex has the stretch modulus of 834 ± 34 pN and a very weak twist-stretch coupling. Our quantitative analyses indicate that, compared with DNA and RNA duplexes, the different flexibility of the DRH duplex in stretching and twist-stretch coupling is mainly attributed to its apparently different basepair inclination in the helical structure.

Non-programmed transcriptional frameshifting is common and highly RNA polymerase type-dependent

Background

The viral or host systems for a gene expression assume repeatability of the process and high quality of the protein product. Since level and fidelity of transcription primarily determines the overall efficiency, all factors contributing to their decrease should be identified and optimized. Among many observed processes, non-programmed insertion/deletion (indel) of nucleotide during transcription (slippage) occurring at homopolymeric A/T sequences within a gene can considerably impact its expression. To date, no comparative study of the most utilized Escherichia coli and T7 bacteriophage RNA polymerases (RNAP) propensity for this type of erroneous mRNA synthesis has been reported. To address this issue we evaluated the influence of shift-prone A/T sequences by assessing indel-dependent phenotypic changes. RNAP-specific expression profile was examined using two of the most potent promoters, P_araBAD of E. coli and φ10 of phage T7.

Results

Here we report on the first systematic study on requirements for efficient transcriptional slippage by T7 phage and cellular RNAPs considering three parameters: homopolymer length, template type, and frameshift directionality preferences. Using a series of out-of-frame gfp reporter genes fused to a variety of A/T homopolymeric sequences we show that T7 RNAP has an exceptional potential for generating frameshifts and is capable of slipping on as few as three adenine or four thymidine residues in a row, in a flanking sequence-dependent manner. In contrast, bacterial RNAP exhibits a relatively low ability to baypass indel mutations and requires a run of at least 7 tymidine and even more adenine residues. This difference comes from involvement of various intrinsic proofreading properties. Our studies demonstrate distinct preference towards a specific homopolymer in slippage induction. Whereas insertion slippage performed by T7 RNAP (but not deletion) occurs tendentiously on poly(A) rather than on poly(T) runs, strong bias towards poly(T) for the host RNAP is observed.

Conclusions

Intrinsic RNAP slippage properties involve trade-offs between accuracy, speed and processivity of transcription. Viral T7 RNAP manifests far greater inclinations to the transcriptional slippage than E. coli RNAP. This possibly plays an important role in driving bacteriophage adaptation and therefore could be considered as beneficial. However, from biotechnological and experimental viewpoint, this might create some problems, and strongly argues for employing bacterial expression systems, stocked with proofreading mechanisms.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1034-4) contains supplementary material, which is available to authorized users.

Role of nucleotide identity in effective CRISPR target escape mutations

Prokaryotes use primed CRISPR adaptation to update their memory bank of spacers against invading genetic elements that have escaped CRISPR interference through mutations in their protospacer target site. We previously observed a trend that nucleotide-dependent mismatches between crRNA and the protospacer strongly influence the efficiency of primed CRISPR adaptation. Here we show that guanine-substitutions in the target strand of the protospacer are highly detrimental to CRISPR interference and interference-dependent priming, while cytosine-substitutions are more readily tolerated. Furthermore, we show that this effect is based on strongly decreased binding affinity of the effector complex Cascade for guanine-mismatched targets, while cytosine-mismatched targets only minimally affect target DNA binding. Structural modeling of Cascade-bound targets with mismatches shows that steric clashes of mismatched guanines lead to unfavorable conformations of the RNA-DNA duplex. This effect has strong implications for the natural selection of target site mutations that lead to effective escape from type I CRISPR-Cas systems.

Pseudouridines have context-dependent mutation and stop rates in high-throughput sequencing

The abundant RNA modification pseudouridine (Ψ) has been mapped transcriptome-wide by chemically modifying pseudouridines with carbodiimide and detecting the resulting reverse transcription stops in high-throughput sequencing. However, these methods have limited sensitivity and specificity, in part due to the use of reverse transcription stops. We sought to use mutations rather than just stops in sequencing data to identify pseudouridine sites. Here, we identify reverse transcription conditions that allow read-through of carbodiimide-modified pseudouridine (CMC-Ψ), and we show that pseudouridines in carbodiimide-treated human ribosomal RNA have context-dependent mutation and stop rates in high-throughput sequencing libraries prepared under these conditions. Furthermore, accounting for the context-dependence of mutation and stop rates can enhance the detection of pseudouridine sites. Similar approaches could contribute to the sequencing-based detection of many RNA modifications.

Dynamic and Energetic Signatures of Adenine Tracts in a rA-dT RNA-DNA Hybrid and in Homologous RNA-DNA, RNA-RNA, and DNA-DNA Double Helices

Nuclear magnetic resonance spectroscopy and proton exchange are being used to characterize the opening reactions of individual base pairs in the RNA-DNA hybrid 5'-rGCGAUAAAAAGGCC-3'/5'-dGGCCTTTTTATCGC-3'. The hybrid contains a central tract of five rA-dT base pairs. The rates and the equilibrium constant of the opening reaction for each base pair are determined from the dependence of the exchange rates of imino protons on ammonia concentration, at 10 °C. The results are compared to those previously obtained by our laboratory for three homologous duplexes of the same base sequence (except for the appropriate T/U substitution), containing tracts of dA-rU, rA-rU, or dA-dT base pairs. The rA-dT tract is distinguished by an enhanced propensity of the base pairs to exist in the extrahelical state. The opening rates of rA-dT base pairs also exhibit a strong dependence on the location of the base pair in the structure; namely, as one advances into the tract, the opening rates of rA-dT base pairs gradually decrease. The local stability of each rA-dT base pair within the tract is the same as that of the corresponding rA-rU base pair in the homologous RNA-only duplex but differs from the stabilities of dA-dT and dA-rU base pairs in the other two duplexes (namely, dA-dT > rA-dT > dA-rU). These results demonstrate that, in nucleic acid double helices with the same base sequence, the opening dynamics and the energetics of individual base pairs are strongly influenced by the nature of the strand and by the structural context of the base pair.

A Biophysical Model of CRISPR/Cas9 Activity for Rational Design of Genome Editing and Gene Regulation

The ability to precisely modify genomes and regulate specific genes will greatly accelerate several medical and engineering applications. The CRISPR/Cas9 (Type II) system binds and cuts DNA using guide RNAs, though the variables that control its on-target and off-target activity remain poorly characterized. Here, we develop and parameterize a system-wide biophysical model of Cas9-based genome editing and gene regulation to predict how changing guide RNA sequences, DNA superhelical densities, Cas9 and crRNA expression levels, organisms and growth conditions, and experimental conditions collectively control the dynamics of dCas9-based binding and Cas9-based cleavage at all DNA sites with both canonical and non-canonical PAMs. We combine statistical thermodynamics and kinetics to model Cas9:crRNA complex formation, diffusion, site selection, reversible R-loop formation, and cleavage, using large amounts of structural, biochemical, expression, and next-generation sequencing data to determine kinetic parameters and develop free energy models. Our results identify DNA supercoiling as a novel mechanism controlling Cas9 binding. Using the model, we predict Cas9 off-target binding frequencies across the lambdaphage and human genomes, and explain why Cas9's off-target activity can be so high. With this improved understanding, we propose several rules for designing experiments for minimizing off-target activity. We also discuss the implications for engineering dCas9-based genetic circuits.

Differential stabilities and sequence-dependent base pair opening dynamics of Watson-Crick base pairs with 5-hydroxymethylcytosine, 5-formylcytosine, or 5-carboxylcytosine

5-Hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), and 5-carboxylcytosine (5caC) form during active demethylation of 5-methylcytosine (5mC) and are implicated in epigenetic regulation of the genome. They are differentially processed by thymine DNA glycosylase (TDG), an enzyme involved in active demethylation of 5mC. Three modified Dickerson-Drew dodecamer (DDD) sequences, amenable to crystallographic and spectroscopic analyses and containing the 5'-CG-3' sequence associated with genomic cytosine methylation, containing 5hmC, 5fC, or 5caC placed site-specifically into the 5'-T(8)X(9)G(10)-3' sequence of the DDD, were compared. The presence of 5caC at the X(9) base increased the stability of the DDD, whereas 5hmC or 5fC did not. Both 5hmC and 5fC increased imino proton exchange rates and calculated rate constants for base pair opening at the neighboring base pair A(5):T(8), whereas 5caC did not. At the oxidized base pair G(4):X(9), 5fC exhibited an increase in the imino proton exchange rate and the calculated kop. In all cases, minimal effects to imino proton exchange rates occurred at the neighboring base pair C(3):G(10). No evidence was observed for imino tautomerization, accompanied by wobble base pairing, for 5hmC, 5fC, or 5caC when positioned at base pair G(4):X(9); each favored Watson-Crick base pairing. However, both 5fC and 5caC exhibited intranucleobase hydrogen bonding between their formyl or carboxyl oxygens, respectively, and the adjacent cytosine N(4) exocyclic amines. The lesion-specific differences observed in the DDD may be implicated in recognition of 5hmC, 5fC, or 5caC in DNA by TDG. However, they do not correlate with differential excision of 5hmC, 5fC, or 5caC by TDG, which may be mediated by differences in transition states of the enzyme-bound complexes.

The presence of an RNA:DNA hybrid that is prone to slippage promotes termination by T7 RNA polymerase

Intrinsic termination signals for multisubunit bacterial RNA polymerases (RNAPs) encode a GC-rich stem-loop structure followed by a polyuridine [poly(U)] tract, and it has been proposed that steric clash of the stem-loop with the exit pore of the RNAP imposes a shearing force on the RNA in the downstream RNA:DNA hybrid, resulting in misalignment of the active site. The structurally unrelated T7 RNAP terminates at a similar type of signal (TΦ), suggesting a common mechanism for termination. In the absence of a hairpin (passive conditions), T7 RNAP slips efficiently in both homopolymeric A and U tracts, and we have found that replacement of the U tract in TΦ with a slippage-prone A tract still allows efficient termination. Under passive conditions, incorporation of a single G residue following a poly(U) tract (which is the situation during termination at TΦ) results in a "locked" complex that is unable to extend the transcript. Our results support a model in which transmission of the shearing force generated by steric clash of the hairpin with the exit pore is promoted by the presence of a slippery tracts downstream, resulting in alterations in the active site and the formation of a locked complex that represents an early step in the termination pathway.

Bacteriophage λ N protein inhibits transcription slippage by Escherichia coli RNA polymerase

Transcriptional slippage is a class of error in which ribonucleic acid (RNA) polymerase incorporates nucleotides out of register, with respect to the deoxyribonucleic acid (DNA) template. This phenomenon is involved in gene regulation mechanisms and in the development of diverse diseases. The bacteriophage λ N protein reduces transcriptional slippage within actively growing cells and in vitro. N appears to stabilize the RNA/DNA hybrid, particularly at the 5′ end, preventing loss of register between transcript and template. This report provides the first evidence of a protein that directly influences transcriptional slippage, and provides a clue about the molecular mechanism of transcription termination and N-mediated antitermination.

DNA–RNA hybrid duplexes with decreasing pyrimidine content in the DNA strand provide structural snapshots for the A- to B-form conformational transition of nucleic acids

A gradual increase in the deoxypyrimidine content in DNA–RNA hybrids leads to B- to A-form nucleic acid transition. Possible factors that govern nuclease activity on hybrid duplexes are presented.

Site-specific stabilization of DNA by a tethered major groove amine, 7-aminomethyl-7-deaza-2'-deoxyguanosine

A cationic 7-aminomethyl-7-deaza-2′-deoxyguanosine (7amG) was incorporated site-specifically into the self-complementary duplex d(G¹A²G³A⁴X⁵C⁶G⁷C⁸T⁹C¹⁰T¹¹C¹²)₂ (X = 7amG). This construct placed two positively charged amines adjacent to the major groove edges of two symmetry-related guanines, providing a model for probing how cation binding in the major groove modulates the structure and stability of DNA. Molecular dynamics calculations restrained by nuclear magnetic resonance (NMR) data revealed that the tethered cationic amines were in plane with the modified base pairs. The tethered amines did not form salt bridges to the phosphodiester backbone. There was also no indication of the amines being capable of hydrogen bonding to flanking DNA bases. NMR spectroscopy as a function of temperature revealed that the X⁵ imino resonance remained sharp at 55 °C. Additionally, two 5′-neighboring base pairs, A⁴:T⁹ and G³:C¹⁰, were stabilized with respect to the exchange of their imino protons with solvent. The equilibrium constant for base pair opening at the A⁴:T⁹ base pair determined by magnetization transfer from water in the absence and presence of added ammonia base catalyst decreased for the modified duplex compared to that of the A⁴:T⁹ base pair in the unmodified duplex, which confirmed that the overall fraction of the A⁴:T⁹ base pair in the open state of the modified duplex decreased. This was also observed for the G³:C¹⁰ base pair, where αK_op for the G³:C¹⁰ base pair in the modified duplex was 3.0 × 10⁶ versus 4.1 × 10⁶ for the same base pair in the unmodified duplex. In contrast, equilibrium constants for base pair opening at the X⁵:C⁸ and C⁶:G⁷ base pairs did not change at 15 °C. These results argue against the notion that electrostatic interactions with DNA are entirely entropic and suggest that major groove cations can stabilize DNA via enthalpic contributions to the free energy of duplex formation.

Conformational flexibility influences degree of hydration of nucleic acid hybrids

Four nucleic acid duplexes-DNA/RNA hybrid, RNA/DNA hybrid, RNA duplex, and DNA duplex-were studied under molecular crowding conditions of osmolytes. Destabilization of duplexes (ΔΔG°(25)) indicated that the ΔΔG°(25) values of hybrids were intermediate between those of DNA and RNA duplexes. In the presence of polyethylene glycol 200, the ΔΔG°(25) values were estimated to be +3.0, +3.5, +3.5, and +4.1 kcal mol(-1) for the DNA duplex, DNA/RNA hybrid, RNA/DNA hybrid, and RNA duplex, respectively. Differences in the number of water molecules taken up (-Δn(w)) upon duplex formations between 0 and 37 °C (Δ(-Δn(w))) were estimated to be 44.8 and 59.7 per duplex structure for the DNA/RNA and RNA/DNA hybrids, respectively. While the Δ(-Δn(w)) value for the DNA/RNA hybrid was intermediate between those of the DNA (26.1) and RNA (59.2) duplexes, the value for RNA/DNA hybrid was close to that of RNA duplex. These differences in the thermodynamic parameters and hydration are probably a consequence of the enhanced global flexibility of the RNA/DNA hybrid structure relative to the DNA/RNA hybrid structure observed in molecular dynamics simulations. This molecular crowding study provides information not only on hydration but also on the flexibility of the conformation of nucleic acid duplexes.

Enhanced base-pair opening in the adenine tract of a RNA double helix

Proton exchange and nuclear magnetic resonance spectroscopy are being used to characterize the kinetics and energetics of base-pair opening in two nucleic acid double helices. One is the RNA duplex 5'-r(GCGAUAAAAAGGCC)-3'/5'-r(GGCCUUUUUAUCGC)-3', which contains a central tract of five AU base pairs. The other is the homologous DNA duplex with a central tract of five AT base pairs. The rates and the equilibrium constants of the opening reaction of each base pair are measured from the dependence of the exchange rates of imino protons on ammonia concentration, at 10 °C. The results reveal that the tract of AU base pairs in the RNA duplex differs from the homologous tract of AT base pairs in DNA in several ways. The rates of opening of AU base pairs in RNA are high and increase progressively along the tract, reaching their largest values at the 3'-end of the tract. In contrast, the opening rates of AT base pairs in DNA are much lower than those of AU base pairs. Within the tract, the largest opening rate is observed for the AT base pair at the 5'-end of the tract. These differences in opening kinetics are paralleled by differences in the stabilities of individual base pairs. All AU base pairs in the RNA are less stable than the AT base pairs in the DNA. The presence of the tract enhances these differences by increasing the stability of AT base pairs in DNA while decreasing the stability of AU base pairs in RNA. Due to these divergent trends, along the tracts, the AU base pairs become progressively less stable than AT base pairs. These findings demonstrate that tracts of AU base pairs in RNA have specific dynamic and energetic signatures that distinguish them from similar tracts of AT base pairs in DNA.

Direct observation of the temperature-induced melting process of the Salmonella fourU RNA thermometer at base-pair resolution

In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5'-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine-Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14-C25 base pair. The mismatch base pair in the wt RNA thermometer (A8-G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.

Structural energetics of the adenine tract from an intrinsic transcription terminator

Intrinsic transcription termination sites generally contain a tract of adenines in the DNA template that yields a tract of uracils at the 3' end of the nascent RNA. To understand how this base sequence contributes to termination of transcription, we have investigated two nucleic acid structures. The first is the RNA-DNA hybrid that contains the uracil tract 5'-rUUUUUAU-3' from the tR2 intrinsic terminator of bacteriophage lambda. The second is the homologous DNA-DNA duplex that contains the adenine tract 5'-dATAAAAA-3'. This duplex is present at the tR2 site when the DNA is not transcribed. The opening and the stability of each rU-dA/dT-dA base pair in the two structures are characterized by imino proton exchange and nuclear magnetic resonance spectroscopy. The results reveal concerted opening of the central rU-dA base pairs in the RNA-DNA hybrid. Furthermore, the stability profile of the adenine tract in the RNA-DNA hybrid is very different from that of the tract in the template DNA-DNA duplex. In the RNA-DNA hybrid, the stabilities of rU-dA base pairs range from 4.3 to 6.5 kcal/mol (at 10 degrees C). The sites of lowest stability are identified at the central positions of the tract. In the template DNA-DNA duplex, the dT-dA base pairs are more stable than the corresponding rU-dA base pairs in the hybrid by 0.9 to 4.6 kcal/mol and, in contrast to the RNA-DNA hybrid, the central base pairs have the highest stability. These results suggest that the central rU-dA/dT-dA base pairs in the adenine tract make the largest energetic contributions to transcription termination by promoting both the dissociation of the RNA transcript and the closing of the transcription bubble. The results also suggest that the high stability of dT-dA base pairs in the DNA provides a signal for the pausing of RNA polymerase at the termination site.

R loops stimulate genetic instability of CTG.CAG repeats

Transcription stimulates the genetic instability of trinucleotide repeat sequences. However, the mechanisms leading to transcription-dependent repeat length variation are unclear. We demonstrate, using biochemical and genetic approaches, that the formation of stable RNA.DNA hybrids enhances the instability of CTG.CAG repeat tracts. In vitro transcribed CG-rich repeating sequences, unlike AT-rich repeats and nonrepeating sequences, form stable, ribonuclease A-resistant structures. These RNA.DNA hybrids are eliminated by ribonuclease H treatment. Mutation in the rnhA1 gene that decreases the activity of ribonuclease HI stimulates the instability of CTG.CAG repeats in E. coli. Importantly, the effect of ribonuclease HI depletion on repeat instability requires active transcription. We also showed that transcription-dependent CTG.CAG repeat instability in human cells is stimulated by siRNA knockdown of RNase H1 and H2. In addition, we used bisulfite modification, which detects single-stranded DNA, to demonstrate that the nontemplate DNA strand at transcribed CTG.CAG repeats remains partially single-stranded in human genomic DNA, thus indicating that it is displaced by an RNA.DNA hybrid. These studies demonstrate that persistent hybrids between the nascent RNA transcript and the template DNA strand at CTG.CAG tracts promote instability of DNA trinucleotide repeats.