Global structure of a DNA three-way junction by solution NMR: towards prediction of 3H fold

Solution and Gas-Phase Stability of DNA Junctions from Temperature-Controlled Electrospray Ionization and Surface-Induced Dissociation

DNA three-way junction (TWJ) structures transiently form during key cellular processes such as transcription, replication, and DNA repair. Despite their significance, the thermodynamics of TWJs, including the influence of strand length, base pair composition, and ligand binding on TWJ stability and dissociation mechanisms, are poorly understood. To address these questions, we interfaced temperature-controlled nanoelectrospray ionization mass spectrometry (TC-nESI-MS) with a cyclic ion mobility spectrometry (cIMS) instrument that was also equipped with a surface-induced dissociation (SID) stage. This novel combination allowed us to investigate the structural intermediates of three TWJ complexes and examine the effects of GC base pairs on their dissociation pathways. We found that two TWJ-specific ligands, 2,7-tris-naphthalene (2,7-TrisNP) and tris-phenoxybenzene (TrisPOB), lead to TWJ stabilization, revealed by an increase in the melting temperature (Tm) by 13 or 26 °C, respectively. To gain insights into conformational changes in the gas phase, we employed cIMS and SID to analyze TWJs and their complexes with ligands. Analysis of IM arrival distributions suggested a single-step dissociation of TWJs and their intermediates for the three studied TWJ complexes. Upon ligand binding, a higher SID energy by 3 V (2,7-TrisNP) and 5 V (TrisPOB) was required to induce 50% dissociation of TWJ, compared to 38 V in the absence of ligands. Our results demonstrate the power of utilizing TC-nESI-MS in combination with cIMS and SID for thermodynamic characterization of TWJ complexes and investigation of ligand binding. These techniques are essential for the TWJ design and development as drug targets, aptamers, and structural units for functional biomaterials.

Solution Nuclear Magnetic Resonance Structures of ATTTT and ATTTC Pentanucleotide Repeats Associated with SCA37 and FAMEs

Expansions of ATTTT and ATTTC pentanucleotide repeats in the human genome are recently found to be associated with at least seven neurodegenerative diseases, including spinocerebellar ataxia type 37 (SCA37) and familial adult myoclonic epilepsy (FAME) types 1, 2, 3, 4, 6, and 7. The formation of non-B DNA structures during some biological processes is thought as a causative factor for repeat expansions. Yet, the structural basis for these pyrimidine-rich ATTTT and ATTTC repeat expansions remains elusive. In this study, we investigated the solution structures of ATTTT and ATTTC repeats using nuclear magnetic resonance spectroscopy. Here, we reveal that ATTTT and ATTTC repeats can form a highly compact minidumbbell structure at the 5'-end using their first two repeats. The high-resolution structure of two ATTTT repeats was determined, showing a regular TTTTA pentaloop and a quasi TTTT/A pentaloop. Furthermore, the minidumbbell structure could escape from proofreading by the Klenow fragment of DNA polymerase I when it was located at five or more base pairs away from the priming site, leading to a small-scale repeat expansion. Results of this work improve our understanding of ATTTT and ATTTC repeat expansions in SCA37 and FAMEs, and provide high-resolution structural information for rational drug design.

Structures and conformational dynamics of DNA minidumbbells in pyrimidine-rich repeats associated with neurodegenerative diseases

Expansions of short tandem repeats (STRs) are associated with approximately 50 human neurodegenerative diseases. These pathogenic STRs are prone to form non-B DNA structure, which has been considered as one of the causative factors for repeat expansions. Minidumbbell (MDB) is a relatively new type of non-B DNA structure formed by pyrimidine-rich STRs. An MDB is composed of two tetraloops or pentaloops, exhibiting a highly compact conformation with extensive loop-loop interactions. The MDB structures have been found to form in CCTG tetranucleotide repeats associated with myotonic dystrophy type 2, ATTCT pentanucleotide repeats associated with spinocerebellar ataxia type 10, and the recently discovered ATTTT/ATTTC repeats associated with spinocerebellar ataxia type 37 and familial adult myoclonic epilepsy. In this review, we first introduce the structures and conformational dynamics of MDBs with a focus on the high-resolution structural information determined by nuclear magnetic resonance spectroscopy. Then we discuss the effects of sequence context, chemical environment, and nucleobase modification on the structure and thermostability of MDBs. Finally, we provide perspectives on further explorations of sequence criteria and biological functions of MDBs.

Binding-driven reactivity attenuation enables NMR identification of selective drug candidates for nucleic acid targets

NMR methods, and in particular ligand-based approaches, are among the most robust and reliable alternatives for binding detection and consequently, they have become highly popular in the context of hit identification and drug discovery. However, when dealing with DNA/RNA targets, these techniques face limitations that have precluded widespread application in medicinal chemistry. In order to expand the arsenal of spectroscopic tools for binding detection and to overcome the existing difficulties, herein we explore the scope and limitations of a strategy that makes use of a binding indicator previously unexploited by NMR: the perturbation of the ligand reactivity caused by complex formation. The obtained results indicate that ligand reactivity can be utilised to reveal association processes and identify the best binders within mixtures of significant complexity, providing a conceptually different reactivity-based alternative within NMR screening methods.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.

Single-tube isothermal label-free fluorescent sensor for pathogen detection based on genetic signatures

We report on a single-tube biosensor for real-time detection of bacterial pathogens with multiplex capabilities. The biosensor consists of two DNA probes, which bind to the complementary fragment of a bacterial RNA to form a three-way junction (3WJ) nucleic acid structure. One of the probes encodes a fluorescent light-up RNA aptamer under T7 promoter. It allows for generation of multiple aptamer copies due to elongation and transcription of the 3WJ structure in the presence of the complementary target. The aptamer coordinates and thereby enhances fluorescence of a cognate fluorogenic dye, allowing for fluorescent detection of the RNA target. Multiple aptamer copies can be produced from a single target-dependent 3WJ structure allowing for amplification and visual observation of the signal. The limit of detection depended on the assay time and was found to be 1.7 nM or 0.6 nM for 30-min or 60-min assay, respectively, when N-methylmesoporphyrin IX (NMM) was used as a fluorescent indicator. The sensor is excellent in analyzing folded RNA targets and differentiating between closely related sequences due to the multicomponent character of the target-interrogating probe. Response to unamplified samples of total bacterial RNA from Mycobacterium tuberculosis complex or Escherichia coli was observed with excellent selectivity within 30 min under isothermal conditions at 50°C in a one-tube one-step assay. Several bacterial species can be detected in multiplex by utilizing biosensors with the template strands encoding different light-up aptamers. The isothermal one-tube-one-step format of the assay and the possibility to monitor the signal visually makes it amenable to use in a point-of-care scenario.

Targeting Pathogenic DNA and RNA Repeats: A Conceptual Therapeutic Way for Repeat Expansion Diseases

Expansions of short tandem repeats (STRs) in the human genome cause nearly 50 neurodegenerative diseases, which are mostly inheritable, nonpreventable and incurable, posing as a huge threat to human health. Non-B DNAs formed by STRs are thought to be structural intermediates that can cause repeat expansions. The subsequent transcripts harboring expanded RNA repeats can further induce cellular toxicity through forming specific structures. Direct targeting of these pathogenic DNA and RNA repeats has emerged as a new potential therapeutic strategy to cure repeat expansion diseases. In this conceptual review, we first introduce the roles of DNA and RNA structures in the genetic instabilities and pathomechanisms of repeat expansion diseases, then describe structural features of DNA and RNA repeats with a focus on the tertiary structures determined by X-ray crystallography and solution nuclear magnetic resonance spectroscopy, and finally discuss recent progress and perspectives of developing chemical tools that target pathogenic DNA and RNA repeats for curing repeat expansion diseases.

Minidumbbell structures formed by ATTCT pentanucleotide repeats in spinocerebellar ataxia type 10

Spinocerebellar ataxia type 10 (SCA10) is a progressive genetic disorder caused by ATTCT pentanucleotide repeat expansions in intron 9 of the ATXN10 gene. ATTCT repeats have been reported to form unwound secondary structures which are likely linked to large-scale repeat expansions. In this study, we performed high-resolution nuclear magnetic resonance spectroscopic investigations on DNA sequences containing two to five ATTCT repeats. Strikingly, we found the first two repeats of all these sequences well folded into highly compact minidumbbell (MDB) structures. The 3D solution structure of the sequence containing two ATTCT repeats was successfully determined, revealing the MDB comprises a regular TTCTA and a quasi TTCT/A pentaloops with extensive stabilizing loop-loop interactions. We further carried out in vitro primer extension assays to examine if the MDB formed in the primer could escape from the proofreading function of DNA polymerase. Results showed that when the MDB was formed at 5-bp or farther away from the priming site, it was able to escape from the proofreading by Klenow fragment of DNA polymerase I and thus retained in the primer. The intriguing structural findings bring about new insights into the origin of genetic instability in SCA10.

Detection of Mg2+-dependent, coaxial stacking rearrangements in a bulged three-way DNA junction by single-molecule FRET

Three-way helical junctions (3WJs) arise in genetic processing, and they have architectural and functional roles in structured nucleic acids. An internal bulge at the junction core allows the helical domains to become oriented into two possible, coaxially stacked conformers. Here, the helical stacking arrangements for a series of bulged, DNA 3WJs were examined using ensemble fluorescence resonance energy transfer (FRET) and single-molecule FRET (smFRET) approaches. The 3WJs varied according to the GC content and sequence of the junction core as well as the pyrimidine content of the internal bulge. Mg²⁺ titration experiments by ensemble FRET show that both stacking conformations have similar Mg²⁺ requirements for folding. Strikingly, smFRET experiments reveal that a specific junction sequence can populate both conformers and that this junction undergoes continual interconversion between the two stacked conformers. These findings will support the development of folding principles for the rational design of functional DNA nanostructures.

Establishment of a universal and rational gene detection strategy through three-way junction-based remote transduction

The polymerase chain reaction and many isothermal amplifications are able to achieve super gene amplification. Unfortunately, most commonly-used transduction methods, such as dye staining and Taqman-like probing, still suffer from shortcomings including false signals or difficult probe design, or are incompatible with multi-analysis. Here a universal and rational gene detection strategy has been established by translating isothermal amplicons to enzyme-free strand displacement circuits via three-way junction-based remote transduction. An assistant transduction probe was imported to form a partial hybrid with the target single-stranded nucleic acid. After systematic optimization the hybrid could serve as an associative trigger to activate a downstream circuit detector via a strand displacement reaction across the three-way junction. By doing so, the detection selectivity can be double-guaranteed through both amplicon-transducer recognition and the amplicon-circuit reaction. A well-optimized circuit can be immediately applied to a new target detection through simply displacing only 10-12 nt on only one component, according to the target. More importantly, this property for the first time enables multi-analysis and logic-analysis in a single reaction, sharing a single fluorescence reporter. In an applicable model, trace amounts of Cronobacter and Enterobacteria genes have been clearly distinguished from samples with no bacteria or one bacterium, with ultra-high sensitivity and selectivity.

Sequence Effect on the Formation of DNA Minidumbbells

The DNA minidumbbell (MDB) is a recently identified non-B structure. The reported MDBs contain two TTTA, CCTG, or CTTG type II loops. At present, the knowledge and understanding of the sequence criteria for MDB formation are still limited. In this study, we performed a systematic high-resolution nuclear magnetic resonance (NMR) and native gel study to investigate the effect of sequence variations in tandem repeats on the formation of MDBs. Our NMR results reveal the importance of hydrogen bonds, base-base stacking, and hydrophobic interactions from each of the participating residues. We conclude that in the MDBs formed by tandem repeats, C-G loop-closing base pairs are more stabilizing than T-A loop-closing base pairs, and thymine residues in both the second and third loop positions are more stabilizing than cytosine residues. The results from this study enrich our knowledge on the sequence criteria for the formation of MDBs, paving a path for better exploring their potential roles in biological systems and DNA nanotechnology.

Conformational Heterogeneity in a Fully Complementary DNA Three-Way Junction with a GC-Rich Branchpoint

DNA three-way junctions (3WJs) are branched structures that serve as important biological intermediates and as components in DNA nanostructures. We recently derived the global structure of a fully complementary 3WJ and found that it contained unpaired bases at the branchpoint, which is consistent with previous observations of branch flexibility and branchpoint reactivity. By combining high-resolution single-molecule Förster resonance energy transfer, molecular modeling, time-resolved ensemble fluorescence spectroscopy, and the first ¹⁹F nuclear magnetic resonance observations of fully complementary 3WJs, we now show that the 3WJ structure can adopt multiple distinct conformations depending upon the sequence at the branchpoint. A 3WJ with a GC-rich branchpoint adopts an open conformation with unpaired bases at the branch and at least one additional conformation with an increased number of base interactions at the branchpoint. This structural diversity has implications for branch interactions and processing in vivo and for technological applications.

Insights into Watson-Crick/Hoogsteen breathing dynamics and damage repair from the solution structure and dynamic ensemble of DNA duplexes containing m1A

In the canonical DNA double helix, Watson-Crick (WC) base pairs (bps) exist in dynamic equilibrium with sparsely populated (∼0.02-0.4%) and short-lived (lifetimes ∼0.2-2.5 ms) Hoogsteen (HG) bps. To gain insights into transient HG bps, we used solution-state nuclear magnetic resonance spectroscopy, including measurements of residual dipolar couplings and molecular dynamics simulations, to examine how a single HG bp trapped using the N1-methylated adenine (m1A) lesion affects the structural and dynamic properties of two duplexes. The solution structure and dynamic ensembles of the duplexes reveals that in both cases, m1A forms a m1A•T HG bp, which is accompanied by local and global structural and dynamic perturbations in the double helix. These include a bias toward the BI backbone conformation; sugar repuckering, major-groove directed kinking (∼9°); and local melting of neighboring WC bps. These results provide atomic insights into WC/HG breathing dynamics in unmodified DNA duplexes as well as identify structural and dynamic signatures that could play roles in m1A recognition and repair.

Flexibility and conformation of the cocaine aptamer studied by PELDOR

The cocaine aptamer is a DNA three-way junction that binds cocaine at its helical junction. We studied the global conformation and overall flexibility of the aptamer in the absence and presence of cocaine by pulsed electron-electron double resonance (PELDOR) spectroscopy, also called double electron-electron resonance (DEER). The rigid nitroxide spin label Ç was incorporated pairwise into two helices of the aptamer. Multi-frequency 2D PELDOR experiments allow the determination of the mutual orientation and the distances between two Çs. Since Ç is rigidly attached to double-stranded DNA, it directly reports on the aptamer dynamics. The cocaine-bound and the non-bound states could be differentiated by their conformational flexibility, which decreases upon binding to cocaine. We observed a small change in the width and mean value of the distance distribution between the two spin labels upon cocaine binding. Further structural insights were obtained by investigating the relative orientation between the two spin-labeled stems of the aptamer. We determined the bend angle between this two stems. By combining the orientation information with a priori knowledge about the secondary structure of the aptamer, we obtained a molecular model describing the global folding and flexibility of the cocaine aptamer.

Conformationally Gated Charge Transfer in DNA Three-Way Junctions

Molecular structures that direct charge transport in two or three dimensions possess some of the essential functionality of electrical switches and gates. We use theory, modeling, and simulation to explore the conformational dynamics of DNA three-way junctions (TWJs) that may control the flow of charge through these structures. Molecular dynamics simulations and quantum calculations indicate that DNA TWJs undergo dynamic interconversion among "well stacked" conformations on the time scale of nanoseconds, a feature that makes the junctions very different from linear DNA duplexes. The studies further indicate that this conformational gating would control charge flow through these TWJs, distinguishing them from conventional (larger size scale) gated devices. Simulations also find that structures with polyethylene glycol linking groups ("extenders") lock conformations that favor CT for 25 ns or more. The simulations explain the kinetics observed experimentally in TWJs and rationalize their transport properties compared with double-stranded DNA.

Interconverting conformations of slipped-DNA junctions formed by trinucleotide repeats affect repair outcome

Expansions of (CTG)·(CAG) repeated DNAs are the mutagenic cause of 14 neurological diseases, likely arising through the formation and processing of slipped-strand DNAs. These transient intermediates of repeat length mutations are formed by out-of-register mispairing of repeat units on complementary strands. The three-way slipped-DNA junction, at which the excess repeats slip out from the duplex, is a poorly understood feature common to these mutagenic intermediates. Here, we reveal that slipped junctions can assume a surprising number of interconverting conformations where the strand opposite the slip-out either is fully base paired or has one or two unpaired nucleotides. These unpaired nucleotides can also arise opposite either of the nonslipped junction arms. Junction conformation can affect binding by various structure-specific DNA repair proteins and can also alter correct nick-directed repair levels. Junctions that have the potential to contain unpaired nucleotides are repaired with a significantly higher efficiency than constrained fully paired junctions. Surprisingly, certain junction conformations are aberrantly repaired to expansion mutations: misdirection of repair to the non-nicked strand opposite the slip-out leads to integration of the excess slipped-out repeats rather than their excision. Thus, slipped-junction structure can determine whether repair attempts lead to correction or expansion mutations.

Expanding the rule set of DNA circuitry with associative toehold activation

Toehold-mediated strand displacement has proven extremely powerful in programming enzyme-free DNA circuits and DNA nanomachines. To achieve multistep, autonomous, and complex behaviors, toeholds must be initially inactivated by hybridizing to inhibitor strands or domains and then relieved from inactivation in a programmed, timed manner. Although powerful and reasonably robust, this strategy has several drawbacks that limit the architecture of DNA circuits. For example, the combination between toeholds and branch migration (BM) domains is 'hard wired' during DNA synthesis thus cannot be created or changed during the execution of DNA circuits. To solve this problem, I propose a strategy called 'associative toehold activation', where the toeholds and BM domains are connected via hybridization of auxiliary domains during the execution of DNA circuits. Bulged thymidines that stabilize DNA three-way junctions substantially accelerate strand displacement reactions in this scheme, allowing fast strand displacement initiated by reversible toehold binding. To demonstrate the versatility of the scheme, I show (1) run-time combination of toeholds and BM domains, (2) run-time recombination of toeholds and BM domains, which results in a novel operation 'toehold switching', and (3) design of a simple conformational self-replicator.

Small-molecule-triggered manipulation of DNA three-way junctions

DNA three-way junctions are frequently used in nanoarchitectures. Ligand-dependent designs that provide well-characterized building blocks for structure-switching DNA nanodevices are presented.

Precise mapping of RNA tertiary structure via nanometer distance measurements with double electron-electron resonance spectroscopy

Divalent metal (Mg(2+)) ion-dependent folding of the hammerhead ribozyme from Schistosoma mansoni was monitored with double electron-electron resonance (DEER) pulsed electron paramagnetic resonance spectroscopy by measuring nanometer-scale distances between paramagnetic spin-labels attached to the RNA. DEER measurements detect global folding of the ribozyme with excellent agreement between predictions from experimental, modeled, and crystallographic measurements. These measurements demonstrate the use of DEER spectroscopy as a tool for structural analysis of complex RNAs.

The unstable part of the apical stem of duck hepatitis B virus epsilon shows enhanced base pair opening but not pico- to nanosecond dynamics and is essential for reverse transcriptase binding

Hepatitis B virus (HBV) replication starts with binding of reverse transcriptase (RT) to the apical stem-loop region of epsilon, a conserved element of the RNA pregenome. For duck HBV, an in vitro replication system has provided molecular details of this interaction. Further insights can be obtained from the structure and dynamics of the duck and human apical stem-loops. Previously, we reported these for the human apical stem-loop. Here, we present the same for the duck counterpart. Unlike its human counterpart, the duck apical stem is unstable in its middle/upper part and contains noncanonical base pairs. This dynamics study is the first of an unstable RNA-DNA stem. Similar to the human stem, the duck apical stem comprises two helical segments with a bend angle of ca. 10 degrees , separated by a nonpaired mobile U residue. It is capped by a well-structured conserved UGUU loop with two residues mobile on the pico- to nanosecond time scale, one of which is involved in RT binding. Remarkably, the unstable middle/upper part of the stem does not show enhanced pico- to nanosecond time scale dynamics. Instead, adenine dispersion relaxation studies indicate enhanced millisecond time scale dynamics involving base pair opening. It can then be concluded that base pair opening is essential for epsilon-RT binding, because stabilization of the stem abolishes binding. We hypothesize that binding occurs by conformational capture of bases in the base pair open state. The unstable secondary structure of the apical stem-loop makes duck epsilon-RT binding unusual in light of recent classifications of RNA target interactions that assume stable secondary structures.

Stabilization of three-way junctions of DNA under molecular crowding conditions

We examined the effects of molecular crowding conditions on the structures and thermodynamics of three-way junctions (TWJs) of DNA. Structural analysis utilizing gel electrophoresis and circular dichroism spectroscopy showed that the designed DNAs folded into TWJ structures in the presence of Na(+) and Mg(2+) under both dilute and molecular crowding conditions with polyethylene glycol 200 (PEG 200). From the thermodynamic parameters evaluated by UV melting techniques in the absence and presence of 5 mM Mg(2+) under dilute and molecular crowding conditions, it was clear that Mg(2+) stabilized all TWJs under the dilute condition, although the extent of stabilization depended on the stacking partners of TWJs. For example, thermodynamic stability (-DeltaG(o) (37)) of A/B-stacked TWJs (A, B, and C are the three helices of TWJ, and among these helices, A and B are stacked together) increased from 3.7 to 5.6 kcal/mol by the addition of 5 mM Mg(2+), while that of A/C-stacked TWJs (A and C are stacked together) increased only from 3.0 to 3.7 kcal/mol. Molecular crowding with PEG 200 destabilized the whole TWJ consisting of a junction point and three helical duplex arms. Crowding agents such as PEG 200 can affect the stability of DNA by modulating its hydration. To explore the crowding effects on the junction point, we evaluated the number of water molecules associated with the whole TWJ as well as the individual arms, and we found that the number of water molecules taken up by the whole TWJ was significantly smaller than the sum of the individual arms. These results show the dehydration from the junction point of the TWJ structure. Therefore, molecular crowding should be favorable for the junction point of TWJ structure and unfavorable for the duplex structure. To prove this concept, we designed truncated TWJ structures that folded into a bimolecular duplex under the dilute condition. With increasing concentrations of PEG 200 from 0 to 30 wt %, the fraction of truncated TWJ structures gradually increased, and that of the bimolecular duplex structure decreased, even in the absence of Mg(2+). We concluded that a cell-mimicking condition, in which the activity of water decreases and hydration becomes less favorable, might facilitate the formation of junction structures in comparison with duplexes.

Preparation of selective and segmentally labeled single-stranded DNA for NMR by self-primed PCR and asymmetrical endonuclease double digestion

We demonstrate a new, efficient and easy-to-use method for enzymatic synthesis of (stereo-)specific and segmental (13)C/(15)N/(2)H isotope-labeled single-stranded DNA in amounts sufficient for NMR, based on the highly efficient self-primed PCR. To achieve this, new approaches are introduced and combined. (i) Asymmetric endonuclease double digestion of tandem-repeated PCR product. (ii) T4 DNA ligase mediated ligation of two ssDNA segments. (iii) In vitro dNTP synthesis, consisting of in vitro rNTP synthesis followed by enzymatic stereo-selective reduction of the C2' of the rNTP, and a one-pot add-up synthesis of dTTP from dUTP. The method is demonstrated on two ssDNAs: (i) a 36-nt three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC labeled and (ii) a 39-nt triple-repeat three-way junction, selectively (13)C(9)/(15)N(3)/(2)H((1',2'',3',4',5',5''))-dC and (13)C(9)/(15)N(2)/(2)H((1',2'',3',4',5',5''))-dT labeled in segment C20-C39. Their NMR spectra show the spectral simplification, while the stereo-selective (2)H-labeling in the deoxyribose of the dC-residues, straightforwardly provided assignment of their C1'-H2' and C2'-H2' resonances. The labeling protocols can be extended to larger ssDNA molecules and to more than two segments.

Folding of the cocaine aptamer studied by EPR and fluorescence spectroscopies using the bifunctional spectroscopic probe Ç

The cocaine aptamer is a DNA molecule that binds cocaine at the junction of three helices. The bifunctional spectroscopic probe Ç was incorporated independently into three different positions of the aptamer and changes in structure and dynamics upon addition of the cocaine ligand were studied. Nucleoside Ç contains a rigid nitroxide spin label and can be studied directly by electron paramagnetic resonance (EPR) spectroscopy and fluorescence spectroscopy after reduction of the nitroxide to yield the fluoroside Ç(f). Both the EPR and the fluorescence data for aptamer 2 indicate that helix III is formed before cocaine binding. Upon addition of cocaine, increased fluorescence of a fully base-paired Ç(f), placed at the three-way junction in helix III, was observed and is consistent with a helical tilt from a coaxial stack of helices II and III. EPR and fluorescence data clearly show that helix I is formed upon addition of cocaine, concomitant with the formation of the Y-shaped three-way helical junction. The EPR data indicate that nucleotides in helix I are more mobile than nucleotides in regular duplex regions and may reflect increased dynamics due to the short length of helix I.

Recognition of DNA three-way junctions by metallosupramolecular cylinders: gel electrophoresis studies

The interaction of metallosupramolecular cylinders with DNA three-way junctions has been studied by gel electrophoresis. A recent X-ray crystal structure of a palindromic oligonucleotide forming part of a complex with such a cylinder revealed binding at the heart of a three-way junction structure. The studies reported herein confirm that this is not solely an artefact of crystallisation and reveal that this is a potentially very powerful new mode of DNA recognition with wide scope. The cylinders are much more effective at stabilizing three-way junctions than simple magnesium di-cations or organic or metallo-organic tetra-cations, with the M cylinder enantiomer being more effective than P. The recognition is not restricted to three-way junctions formed from palindromic DNA with a central AT step at the junction; non-palindromic three-way junctions and those with GC steps are also stabilised. The cylinder is also revealed to stabilise other Y-shaped junctions, such as that formed at a fraying point in duplex DNA (for example, a replication fork), and other DNA three-way junction structures, such as those containing unpaired nucleotides, perhaps by opening up this structure to access the central cavity.

Thermodynamics and NMR studies on Duck, Heron and Human HBV encapsidation signals

Hepatitis B virus (HBV) replication is initiated by binding of its reverse transcriptase (P) to the apical stem-loop (AL) and primer loop (PL) of epsilon, a highly conserved RNA element at the 5'-end of the RNA pregenome. Mutation studies on duck/heron and human in vitro systems have shown similarities but also differences between their P-epsilon interaction. Here, NMR and UV thermodynamic data on AL (and PL) from these three species are presented. The stabilities of the duck and heron ALs were found to be similar, and much lower than that of human. NMR data show that this low stability stems from an 11-nt internal bulge destabilizing the stem of heron AL. In duck, although structured at low temperature, this region also forms a weak point as its imino resonances broaden to disappearance between 30 and 35 degrees C well below the overall AL melting temperature. Surprisingly, the duck- and heron ALs were both found to be capped by a stable well-structured UGUU tetraloop. All avian ALs are expected to adhere to this because of their conserved sequence. Duck PL is stable and structured and, in view of sequence similarities, the same is expected for heron - and human PL.

Prediction of molecular alignment of nucleic acids in aligned media

We demonstrate--using the data base of all deposited DNA and RNA structures aligned in Pf1-medium and RDC refined--that for nucleic acids in a Pf1-medium the electrostatic alignment tensor can be predicted reliably and accurately via a simple and fast calculation based on the gyration tensor spanned out by the phosphodiester atoms. The rhombicity is well predicted over its full range from 0 to 0.66, while the alignment tensor orientation is predicted correctly for rhombicities up to ca. 0.4, for larger rhombicities it appears to deviate somewhat more than expected based on structural noise and measurement error. This simple analytical approach is based on the Debye-Huckel approximation for the electrostatic interaction potential, valid at distances sufficiently far away from a poly-ionic charged surface, a condition naturally enforced when the charge of alignment medium and solute are of equal sign, as for nucleic acids in a Pf1-phage medium. For the usual salt strengths and nucleic acid sizes, the Debye-Huckel screening length is smaller than the nucleic acid size, but large enough for the collective of Debye-Huckel spheres to encompass the whole molecule. The molecular alignment is then purely electrostatic, but it's functional form is under these conditions similar to that for steric alignment. The proposed analytical expression allows for very fast calculation of the alignment tensor and hence RDCs from the conformation of the nucleic acid molecule. This information provides opportunities for improved structure determination of nucleic acids, including better assessment of dynamics in (multi-domain) nucleic acids and the possibility to incorporate alignment tensor prediction from shape directly into the structure calculation process. The procedures are incorporated into MATLAB scripts, which are available on request.

Role of LrpC from Bacillus subtilis in DNA transactions during DNA repair and recombination

Bacillus subtilis LrpC is a sequence-independent DNA-binding and DNA-bending protein, which binds both single-stranded (ss) and double-stranded (ds) DNA and facilitates the formation of higher order protein–DNA complexes in vitro. LrpC binds at different sites within the same DNA molecule promoting intramolecular ligation. When bound to separate molecules, it promotes intermolecular ligation, and joint molecule formation between a circular ssDNA and a homologous ssDNA-tailed linear dsDNA. LrpC binding showed a higher affinity for 4-way (Holliday) junctions in their open conformation, when compared with curved dsDNA. Consistent with these biochemical activities, an lrpC null mutant strain rendered cells sensitive to DNA damaging agents such as methyl methanesulfonate and 4-nitroquinoline-1-oxide, and showed a segregation defect. These findings collectively suggest that LrpC may be involved in DNA transactions during DNA repair and recombination.

Compensating bends in a 16-base-pair DNA oligomer containing a T(3)A(3) segment: A NMR study of global DNA curvature

In-phase ligated DNA containing T(n)A(n) segments fail to exhibit the retarded polyacrylamide gel electrophoresis (PAGE) migration observed for in-phase ligated A(n)T(n) segments, a behavior thought to be correlated with macroscopic DNA curvature. The lack of macroscopic curvature in ligated T(n)A(n) segments is thought to be due to cancellation of bending in regions flanking the TpA steps. To address this issue, solution-state NMR, including residual dipolar coupling (RDC) restraints, was used to determine a high-resolution structure of [d(CGAGGTTTAAACCTCG)2], a DNA oligomer containing a T3A3 tract. The overall magnitude and direction of bending, including the regions flanking the central TpA step, was measured using a radius of curvature, Rc, analysis. The Rc for the overall molecule indicated a small magnitude of global bending (Rc = 138 +/- 23 nm) towards the major groove, whereas the Rc for the two halves (72 +/- 33 nm and 69 +/- 14 nm) indicated greater localized bending into the minor groove. The direction of bending in the regions flanking the TpA step is in partial opposition (109 degrees), contributing to cancellation of bending. The cancellation of bending did not correlate with a pattern of roll values at the TpA step, or at the 5' and 3' junctions, of the T3A3 segment, suggesting a simple junction/roll model is insufficient to predict cancellation of DNA bending in all T(n)A(n) junction sequence contexts. Importantly, Rc analysis of structures refined without RDC restraints lacked the precision and accuracy needed to reliably measure bending.