Influence of loop residues on the relative stabilities of DNA hairpin structures

Force-induced melting of DNA hairpin: Unfolding pathways and phase diagrams

Using the exact enumeration technique, we have studied the force-induced melting of a DNA hairpin on the face centered cubic lattice for two different sequences which differ in terms of loop closing base pairs. The melting profiles obtained from the exact enumeration technique is consistent with the Gaussian network model and Langevin dynamics simulations. Probability distribution analysis based on the exact density of states revealed the microscopic details of the opening of the hairpin. We showed the existence of intermediate states near the melting temperature. We further showed that different ensembles used to model single-molecule force spectroscopy setups may give different force-temperature diagrams. We delineate the possible reasons for the observed discrepancies.

Effects of Local Sequence, Reaction Conditions, and Various Additives on the Formation and Stability of Interstrand Cross-Links Derived from the Reaction of an Abasic Site with an Adenine Residue in Duplex DNA

The experiments described here examined the effects of reaction conditions, various additives, and local sequence on the formation and stability interstrand cross-links (ICLs) derived from the reaction of an apurinic/apyrimidinic (AP) site with the exocyclic amino group of an adenine residue on the opposing strand in duplex DNA. Cross-link formation was observed in a range of different buffers, with faster formation rates observed at pH 5. Inclusion of the base excision repair enzyme alkyladenine DNA glycosylase (hAAG) which binds tightly to AP-containing duplexes decreased, but did not completely prevent, formation of the dA-AP ICL. Formation of the dA-AP ICL was not altered by the presence of the biological metal ion Mg²⁺ or the biological thiol, glutathione. Several organocatalysts of imine formation did not enhance the rate of dA-AP ICL formation. Duplex length did not have a large effect on dA-AP yield, so long as the melting temperature of the duplex was not significantly below the reaction temperature (the duplex must remain hybridized for efficient ICL formation). Formation of the dA-AP ICL was examined in over 40 different sequences that varied the neighboring and opposing bases at the cross-linking site. The results indicate that ICL formation can occur in a wide variety of sequence contexts under physiological conditions. Formation of the dA-AP ICL was strongly inhibited by the aldehyde-trapping agents methoxyamine and hydralazine, by NaBH₃CN, by the intercalator ethidium bromide, and by the minor groove-binding agent netropsin. ICL formation was inhibited to some extent in bicarbonate and Tris buffers. The dA-AP ICL showed substantial inherent stability under a variety of conditions and was not a substrate for AP-processing enzymes APE1 or Endo IV. Finally, we characterized cross-link formation in a small (11 bp) stem-loop (hairpin) structure and in DNA-RNA hybrid duplexes.

Can one detect intermediate denaturation states of DNA sequences by following the equilibrium open-close dynamic fluctuations of a single base pair?

Melting of DNA sequences may occur through a few major intermediate states, whose influence on the melting curve has been discussed previously, while their effect on the kinetics has not been explored thoroughly. Here, we chose a simple DNA sequence, forming a hairpin in its native (zipped) state, and study it using molecular dynamic (MD) simulations and a model integrating the Gaussian network model with bond-binding energies-the Gaussian binding energy (GBE) model. We find two major partial denaturation states, a bubble state and a partial unzipping state. We demonstrate the influence of these two states on the closing-opening base pair dynamics, as probed by a tagged bond auto-correlation function (ACF). We argue that the latter is measured by fluorescence correlation spectroscopy experiments, in which one base of the pair is linked to a fluorescent dye, while the complementary base is linked to a quencher, similar to the experiment reported by Altan-Bonnet et al. [Phys. Rev. Lett. 90, 138101 (2003)]. We find that tagging certain base pairs at temperatures around the melting temperature results in a multi-step relaxation of the ACF, while tagging other base pairs leads to an effectively single-step relaxation, albeit non-exponential. Only the latter type of relaxation has been observed experimentally, and we suggest which of the other base pairs should be tagged in order to observe multi-step relaxation. We demonstrate that this behavior can be observed with other sequences and argue that the GBE can reliably predict these dynamics for very long sequences, where MD simulations might be limited.

Artificial nondirectional site-specific recombination systems

Site-specific recombination systems (SRSs) are widely used in studies on synthetic biology and related disciplines. Nondirectional SRSs can randomly trigger excision, integration, reversal, and translocation, which are effective tools to achieve large-scale genome recombination. In this study, we designed 6 new nondirectional SRSs named Vika/voxsym1-4 and Dre/roxsym1-2. All 6 artificial nondirectional SRSs were able to generate random deletion and inversion in Saccharomyces cerevisiae. Moreover, all six SRSs were orthogonal to Cre/loxPsym. The pairwise orthogonal nondirected SRSs can simultaneously initiate large-scale and independent gene recombination in two different regions of the genome, which could not be accomplished using previous orthogonal systems. These SRSs were found to be robust while working in the cells at different growth stages, as well as in the different spatial structure of the chromosome. These artificial pairwise orthogonal nondirected SRSs offer newfound potential for site-specific recombination in synthetic biology.

•

Designed six new artificial nondirectional site-specific recombination systems

•

Pairwise orthogonal nondirected recombination systems in yeast

•

The deletion efficiency of systems is far greater than the inversion efficiency

•

These nondirectional recombination systems were found to be robust

The effect of hairpin loop on the structure and gene expression activity of the long-loop G-quadruplex

G-quadruplex (G4), a four-stranded DNA or RNA structure containing stacks of guanine tetrads, plays regulatory roles in many cellular functions. So far, conventional G4s containing loops of 1–7 nucleotides have been widely studied. Increasing experimental evidence suggests that unconventional G4s, such as G4s containing long loops (long-loop G4s), play a regulatory role in the genome by forming a stable structure. Other secondary structures such as hairpins in the loop might thus contribute to the stability of long-loop G4s. Therefore, investigation of the effect of the hairpin-loops on the structure and function of G4s is required. In this study, we performed a systematic biochemical investigation of model G4s containing long loops with various sizes and structures. We found that the long-loop G4s are less stable than conventional G4s, but their stability increased when the loop forms a hairpin (hairpin-G4). We also verified the biological significance of hairpin-G4s by showing that hairpin-G4s present in the genome also form stable G4s and regulate gene expression as confirmed by in cellulo reporter assays. This study contributes to expanding the scope and diversity of G4s, thus facilitating future studies on the role of G4s in the human genome.

Decoding the Structural Dynamics and Conformational Alternations of DNA Secondary Structures by Single-Molecule FRET Microspectroscopy

In addition to the canonical double helix form, DNA is known to be extrapolated into several other secondary structural patterns involving themselves in inter- and intramolecular type hydrogen bonding. The secondary structures of nucleic acids go through several stages of multiple, complex, and interconvertible heterogeneous conformations. The journey of DNA through these conformers has significant importance and has been monitored thoroughly to establish qualitative and quantitative information about the transition between the unfolded, folded, misfolded, and partially folded states. During this structural interconversion, there always exist specific populations of intermediates, which are short-lived or sometimes even do not accumulate within a heterogeneous population and are challenging to characterize using conventional ensemble techniques. The single-molecule FRET(sm-FRET) microspectroscopic method has the advantages to overcome these limitations and monitors biological phenomena transpiring at a measurable high rate and balanced stochastically over time. Thus, tracing the time trajectory of a particular molecule enables direct measurement of the rate constant of each transition step, including the intermediates that are hidden in the ensemble level due to their low concentrations. This review is focused on the advantages of the employment of single-molecule Forster's resonance energy transfer (sm-FRET), which is worthwhile to access the dynamic architecture and structural transition of various secondary structures that DNA adopts, without letting the donor of one molecule to cross-talk with the acceptor of any other. We have emphasized the studies performed to explore the states of folding and unfolding of several nucleic acid secondary structures, for example, the DNA hairpin, Holliday junction, G-quadruplex, and i-motif.

Use of anion-exchange HPLC to study DNA conformational polymorphism

Anion-exchange chromatography carried out under non-denaturing conditions is a versatile tool to differentiate DNA conformations. In this work, the utility of this form of HPLC was demonstrated in four examples. The hairpin and duplex forms of d(CG)₉ were readily resolved, which allowed for the studies of the influence of salt on the equilibrium of these two forms of secondary structures. Similarly, the minimum size of Tn in the loop region required for the sequence 5'-d(CCCAA-(T)n-TTGGG)-3' to form hairpin was established to be two nucleotides using anion-exchange HPLC and fluorescence resonance energy transfer. Furthermore, the efficiency of hybridization of partially self-complementary sequences d[(CG)₆N_x] was readily monitored by non-denaturing anion-exchange HPLC. Finally, different structures adopted by quadruplex-forming sequences were resolved in the same manner.

Interaction of a photosensitizer methylene blue with various structural forms (cruciform, bulge duplex and hairpin) of designed DNA sequences

Functionally important, local structural transitions in DNA generate various alternative conformations. Cruciform is one of such alternative DNA structures, usually targeted in genomes by various proteins. Symmetry elements in sequence as inverted repeats are the key factor for cruciform formation, facilitated by the presence of the AT-rich regions. Here, we used biophysical and biochemical techniques such as Gel electrophoresis, Circular dichroism (CD), and UV-thermal melting analysis to explore the structural status of the designed DNA sequences, which had potential to form cruciform structures under physiological conditions. The gel electrophoresis analysis revealed that the designed 53-mer DNA oligonucleotide sequence CR forms an intermolecular bulge duplex with flanking ends, while another sequence CRC adopts an intramolecular hairpin structure with flanking ends. Their equimolar complex (CRCRC) bestowed much-retarded migration due to the formation of a quite intriguing cruciform structure. CD studies confirmed that all the alternative structures (cruciform, bulge duplex, and hairpin with flanking ends) exhibit characteristics of B-DNA type conformation. A triphasic UV-thermal melting curve displayed by the complex formed by the equimolar ratio (CRCRC) is also suggestive of the formation of the cruciform structure. The interaction studies of CR, CRC, and their equimolar complex (1:1) with a photosensitizer methylene blue (MB) indicated that MB could not stabilize the discrete structures formed by CR and CRC sequences, however, the cruciform structure showed a quite significant increment in the melting temperature. Such studies facilitate our understanding of various secondary structures possibly present inside the cell and their interactions with drug/dye molecules.

Interstrand DNA Cross-Links Derived from Reaction of a 2-Aminopurine Residue with an Abasic Site

Efficient methods for the site-specific installation of structurally defined interstrand cross-links in duplex DNA may be useful in a wide variety of fields. The work described here developed a high-yield synthesis of chemically stable interstrand cross-links resulting from a reductive amination reaction between an abasic site and the noncanonical nucleobase 2-aminopurine in duplex DNA. Results from footprinting, liquid chromatography-mass spectrometry, and stability studies support the formation of an N²-alkylamine attachment between the 2-aminopurine residue and the Ap site. The reaction performs best when the 2-aminopurine residue on the opposing strand is offset 1 nt to the 5'-side of the abasic site. The cross-link confers substantial resistance to thermal denaturation (melting). The cross-linking reaction is fast (complete in 4 h), employs only commercially available reagents, and can be used to generate cross-linked duplexes in sufficient quantities for biophysical, structural, and DNA repair studies.

Investigation of hairpin DNA and chelerythrine interaction by a single bio-nanopore sensing interface

Chelerythrine (CHE) is one of the potential drugs for cancer treatments. The interaction between hairpin DNA and CHE has been investigated by spectral and mass spectrometry methods. In this paper, the stability of hairpin DNA with different loop bases and its interaction with CHE were explored with a single α-hemolysin (α-HL) nanopore sensing interface. The results showed that the characteristic current pulses not only relate to the loop composition changes of the hairpin DNA, but also provide interaction information between CHE and the hairpin DNA molecules. The dwell time of current pulses for hairpin DNA was significantly increased (hundreds of ms) due to the addition of CHE, and two characteristic current distributions were recognized for the hairpin with T₃ and C₃ loops. The two characteristic current groups could be ascribed to the hairpin DNA and the ones with CHE. This study indicates that it is possible to study the interaction between single CHE and single hairpin DNA molecules by the single-nanopore sensing interface as an alternative method to conventional spectrometric methods for therapeutic mechanism and drug screening purposes.

Short inverted repeats contribute to localized mutability in human somatic cells

Selected repetitive sequences termed short inverted repeats (SIRs) have the propensity to form secondary DNA structures called hairpins. SIRs comprise palindromic arm sequences separated by short spacer sequences that form the hairpin stem and loop respectively. Here, we show that SIRs confer an increase in localized mutability in breast cancer, which is domain-dependent with the greatest mutability observed within spacer sequences (∼1.35-fold above background). Mutability is influenced by factors that increase the likelihood of formation of hairpins such as loop lengths (of 4-5 bp) and stem lengths (of 7-15 bp). Increased mutability is an intrinsic property of SIRs as evidenced by how almost all mutational processes demonstrate a higher rate of mutagenesis of spacer sequences. We further identified 88 spacer sequences showing enrichment from 1.8- to 90-fold of local mutability distributed across 283 sites in the genome that intriguingly, can be used to inform the biological status of a tumor.

Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets

We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C₅ loop closing the triplex stem was replaced with ^5'-CT_nC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a ^5'-GCAA or ^5'-AACG tetraloop. For the triplexes with a ^5'-CT_nC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.

The collective behavior of spring-like motifs tethered to a DNA origami nanostructure

Dynamic DNA nanotechnology relies on the integration of small switchable motifs at suitable positions of DNA nanostructures, thus enabling the manipulation of matter with nanometer spatial accuracy in a trigger-dependent fashion. Typical examples of such motifs are hairpins, whose elongation into duplexes can be used to perform long-range, translational movements. In this work, we used temperature-dependent FRET spectroscopy to determine the thermal stabilities of distinct sets of hairpins integrated into the central seam of a DNA origami structure. We then developed a hybrid spring model to describe the energy landscape of the tethered hairpins, combining the thermodynamic nearest-neighbor energy of duplex DNA with the entropic free energy of single-stranded DNA estimated using a worm-like chain approximation. We show that the organized scaffolding of multiple hairpins enhances the thermal stability of the device and that the coordinated action of the tethered motors can be used to mechanically unfold a G-quadruplex motif bound to the inner cavity of the origami structure, thus surpassing the operational capabilities of freely diffusing motors. Finally, we increased the complexity of device functionality through the insertion of two sets of parallel hairpins, resulting in four distinct states and in the reversible localization of desired molecules within the reconfigurable regions of the origami architecture.

Sufficient minimal model for DNA denaturation: Integration of harmonic scalar elasticity and bond energies

We study DNA denaturation by integrating elasticity - as described by the Gaussian network model - with bond binding energies, distinguishing between different base pairs and stacking energies. We use exact calculation, within the model, of the Helmholtz free-energy of any partial denaturation state, which implies that the entropy of all formed "bubbles" ("loops") is accounted for. Considering base pair bond removal single events, the bond designated for opening is chosen by minimizing the free-energy difference for the process, over all remaining base pair bonds. Despite of its great simplicity, for several known DNA sequences our results are in accord with available theoretical and experimental studies. Moreover, we report free-energy profiles along the denaturation pathway, which allow to detect stable or meta-stable partial denaturation states, composed of bubble, as local free-energy minima separated by barriers. Our approach allows to study very long DNA strands with commonly available computational power, as we demonstrate for a few random sequences in the range 200-800 base-pairs. For the latter, we also elucidate the self-averaging property of the system. Implications for the well known breathing dynamics of DNA are elucidated.

High-throughput thermal stability assessment of DNA hairpins based on high resolution melting

On the basis of high-resolution melting, a high-throughput approach to measure melting temperatures (T_ms) of short DNA hairpins was developed. With this method, T_ms of thousands of triloop, tetraloop, and pentaloop hairpins involving various loop sequences and various closing base pairs (cbp) were obtained in hours. The stability of triloop hairpins decreased with the change of cbp (5'-3') in the order of c-g > g-c > t-a ≥ a-t, showing that the cbp of 5'-Pyr-Pur-3' (Pyr = pyrimidine, Pur = purine) contributed more stability than 5'-Pur-Pyr-3'. For tetraloop hairpins, GNNA, GNAB, and CNNG (N = A, G, C, or T; B = G, C, or T) were found to be highly stable irrespective of the cbp type. TNNA was also stable in both g-c and a-t families, while CGNA only in the c-g family. Pentaloop hairpins of cTGNAGg, cGNYNAg (Y = T or C) and cCGNNAg were exceptionally stable motifs. In most cases, pyrimidine-rich loops were more favorable to stabilize the whole structure than purine-rich ones. The present approach showed a good performance in assessing the thermal stability of large amounts of DNA hairpins comprehensively. These data are useful to understand the sequence dependence of the stability of DNA secondary structures and promising to improve the structure simulation by consummating basic databases.

Assembly of supramolecular DNA complexes containing both G-quadruplexes and i-motifs by enhancing the G-repeat-bearing capacity of i-motifs

The single-step assembly of supramolecular complexes containing both i-motifs and G-quadruplexes (G4s) is demonstrated. This can be achieved because the formation of four-stranded i-motifs appears to be little affected by certain terminal residues: a five-cytosine tetrameric i-motif can bear ten-base flanking residues. However, things become complex when different lengths of guanine-repeats are added at the 3′ or 5′ ends of the cytosine-repeats. Here, a series of oligomers d(XG_iXC₅X) and d(XC₅XG_iX) (X = A, T or none; i < 5) are designed to study the impact of G-repeats on the formation of tetrameric i-motifs. Our data demonstrate that tetramolecular i-motif structure can tolerate specific flanking G-repeats. Assemblies of these oligonucleotides are polymorphic, but may be controlled by solution pH and counter ion species. Importantly, we find that the sequences d(TG_iAC₅) can form the tetrameric i-motif in large quantities. This leads to the design of two oligonucleotides d(TG₄AC₇) and d(TG^BrGG^BrGAC₇) that self-assemble to form quadruplex supramolecules under certain conditions. d(TG₄AC₇) forms supramolecules under acidic conditions in the presence of K⁺ that are mainly V-shaped or ring-like containing parallel G4s and antiparallel i-motifs. d(TG^BrGG^BrGAC₇) forms long linear quadruplex wires under acidic conditions in the presence of Na⁺ that consist of both antiparallel G4s and i-motifs.

DNA binding of Pd(TC3), a conformable cationic porphyrin with a long-lived triplet state

The goal of this work has been to synthesize and investigate Pd(TC3), an intercalating porphyrin that has conformable substituents capable of groove binding to B-form DNA. (TC3 denotes the doubly deprotonated form of 5,10,15,20-tetra[3-(3'-methylimidazolium-1'-yl)prop-1-yl]porphyrin.) Palladium(ii) is an apt choice for the central metal ion because it remains strictly four-coordinate and provides for a luminescent triplet excited state with a long lifetime. The DNA hosts are hairpin-forming sequences programmed to differ in base composition. Luminescence, absorbance, and circular dichroism results are consistent with the idea that congruent structural reorganization takes place at the host and ligand during uptake. Photoexcitation of DNA-bound Pd(TC3) generates a comparatively modest steady state concentration of singlet oxygen, due to a relatively slow reaction with molecular oxygen in solution. The sheer size of the substituent groups disfavors quenching, but groove-binding interactions compound the problem by inhibiting mobility. The results show how ligand design affects adduct structure as well as function.

DNA meter: Energy tunable, quantitative hybridization assay

We describe a novel hybridization assay that employs a unique class of energy tunable, bulge loop-containing competitor strands (C*) that hybridize to a probe strand (P). Such initial "pre-binding" of a probe strand modulates its effective "availability" for hybridizing to a target site (T). More generally, the assay described here is based on competitive binding equilibria for a common probe strand (P) between such tunable competitor strands (C*) and a target strand (T). We demonstrate that loop variable, energy tunable families of C*P complexes exhibit enhanced discrimination between targets and mismatched targets, thereby reducing false positives/negatives. We refer to a C*P complex between a C* competitor single strand and the probe strand as a "tuning fork," since the C* strand exhibits branch points (forks) at the duplex-bulge interfaces within the complex. By varying the loop to create families of such "tuning forks," one can construct C*P "energy ladders" capable of resolving small differences within the target that may be of biological/functional consequence. The methodology further allows quantification of target strand concentrations, a determination heretofore not readily available by conventional hybridization assays. The dual ability of this tunable assay to discriminate and quantitate targets provides the basis for developing a technology we refer to as a "DNA Meter." Here we present data that establish proof-of-principle for an in solution version of such a DNA Meter. We envision future applications of this tunable assay that incorporate surface bound/spatially resolved DNA arrays to yield enhanced discrimination and sensitivity.

DNA hairpins destabilize duplexes primarily by promoting melting rather than by inhibiting hybridization

The effect of secondary structure on DNA duplex formation is poorly understood. Using oxDNA, a nucleotide level coarse-grained model of DNA, we study how hairpins influence the rate and reaction pathways of DNA hybridzation. We compare to experimental systems studied by Gao et al. (1) and find that 3-base pair hairpins reduce the hybridization rate by a factor of 2, and 4-base pair hairpins by a factor of 10, compared to DNA with limited secondary structure, which is in good agreement with experiments. By contrast, melting rates are accelerated by factors of ∼100 and ∼2000. This surprisingly large speed-up occurs because hairpins form during the melting process, and significantly lower the free energy barrier for dissociation. These results should assist experimentalists in designing sequences to be used in DNA nanotechnology, by putting limits on the suppression of hybridization reaction rates through the use of hairpins and offering the possibility of deliberately increasing dissociation rates by incorporating hairpins into single strands.

Breathing dynamics based parameter sensitivity analysis of hetero-polymeric DNA

We study the parameter sensitivity of hetero-polymeric DNA within the purview of DNA breathing dynamics. The degree of correlation between the mean bubble size and the model parameters is estimated for this purpose for three different DNA sequences. The analysis leads us to a better understanding of the sequence dependent nature of the breathing dynamics of hetero-polymeric DNA. Out of the 14 model parameters for DNA stability in the statistical Poland-Scheraga approach, the hydrogen bond interaction ε(hb)(AT) for an AT base pair and the ring factor ξ turn out to be the most sensitive parameters. In addition, the stacking interaction ε(st)(TA-TA) for an TA-TA nearest neighbor pair of base-pairs is found to be the most sensitive one among all stacking interactions. Moreover, we also establish that the nature of stacking interaction has a deciding effect on the DNA breathing dynamics, not the number of times a particular stacking interaction appears in a sequence. We show that the sensitivity analysis can be used as an effective measure to guide a stochastic optimization technique to find the kinetic rate constants related to the dynamics as opposed to the case where the rate constants are measured using the conventional unbiased way of optimization.

Impact of bulge loop size on DNA triplet repeat domains: Implications for DNA repair and expansion

Repetitive DNA sequences exhibit complex structural and energy landscapes, populated by metastable, noncanonical states, that favor expansion and deletion events correlated with disease phenotypes. To probe the origins of such genotype-phenotype linkages, we report the impact of sequence and repeat number on properties of (CNG) repeat bulge loops. We find the stability of duplexes with a repeat bulge loop is controlled by two opposing effects; a loop junction-dependent destabilization of the underlying double helix, and a self-structure dependent stabilization of the repeat bulge loop. For small bulge loops, destabilization of the underlying double helix overwhelms any favorable contribution from loop self-structure. As bulge loop size increases, the stabilizing loop structure contribution dominates. The role of sequence on repeat loop stability can be understood in terms of its impact on the opposing influences of junction formation and loop structure. The nature of the bulge loop affects the thermodynamics of these two contributions differently, resulting in unique differences in repeat size-dependent minima in the overall enthalpy, entropy, and free energy changes. Our results define factors that control repeat bulge loop formation; knowledge required to understand how this helix imperfection is linked to DNA expansion, deletion, and disease phenotypes.

Effect of pressure on thermal stability of g-quadruplex DNA and double-stranded DNA structures

Pressure is a thermodynamic parameter that can induce structural changes in biomolecules due to a volumetric decrease. Although most proteins are denatured by pressure over 100 MPa because they have the large cavities inside their structures, the double-stranded structure of DNA is stabilized or destabilized only marginally depending on the sequence and salt conditions. The thermal stability of the G-quadruplex DNA structure, an important non-canonical structure that likely impacts gene expression in cells, remarkably decreases with increasing pressure. Volumetric analysis revealed that human telomeric DNA changed by more than 50 cm3 mol-1 during the transition from a random coil to a quadruplex form. This value is approximately ten times larger than that for duplex DNA under similar conditions. The volumetric analysis also suggested that the formation of G-quadruplex DNA involves significant hydration changes. The presence of a cosolute such as poly(ethylene glycol) largely repressed the pressure effect on the stability of G-quadruplex due to alteration in stabilities of the interactions with hydrating water. This review discusses the importance of local perturbations of pressure on DNA structures involved in regulation of gene expression and highlights the potential for application of high-pressure chemistry in nucleic acid-based nanotechnology.

Design and evaluation of an i-motif-based allosteric control mechanism in DNA-hairpin molecular devices

Molecular devices designed to assess and manipulate biologically relevant conditions with required accuracy and precision play an essential role in life sciences research. Incorporating allosteric regulation mechanism is an attractive strategy toward more efficient artificial sensing and switching systems. Herein, we report on a new principle of regulating switching parameters of a DNA-based molecular device based on allosteric interaction between spatially separated hairpin stem and a tetraplexed fragment (i.e., i-motif). We characterized thermodynamic and kinetic effects arising from interaction between functional domains of the device and demonstrated the potential of applying the allosteric control principle for rational design of sensors and switches with precisely defined operational characteristics.

Hydration changes upon DNA folding studied by osmotic stress experiments

The thermal stability of nucleic acid structures is perturbed under the conditions that mimic the intracellular environment, typically rich in inert components and under osmotic stress. We now describe the thermodynamic stability of DNA oligonucleotide structures in the presence of high background concentrations of neutral cosolutes. Small cosolutes destabilize the basepair structures, and the DNA structures consisting of the same nearest-neighbor composition show similar thermodynamic parameters in the presence of various types of cosolutes. The osmotic stress experiments reveal that water binding to flexible loops, unstable mismatches, and an abasic site upon DNA folding are almost negligible, whereas the binding to stable mismatch pairs is significant. The studies using the basepair-mimic nucleosides and the peptide nucleic acid suggest that the sugar-phosphate backbone and the integrity of the basepair conformation make important contributions to the binding of water molecules to the DNA bases and helical grooves. The study of the DNA hydration provides the basis for understanding and predicting nucleic acid structures in nonaqueous solvent systems.

Determining the DNA stability parameters for the breathing dynamics of heterogeneous DNA by stochastic optimization

We suggest that the thermodynamic stability parameters (nearest neighbor stacking and hydrogen bonding free energies) of double-stranded DNA molecules can be inferred reliably from time series of the size fluctuations (breathing) of local denaturation zones (bubbles). On the basis of the reconstructed bubble size distribution, this is achieved through stochastic optimization of the free energies in terms of simulated annealing. In particular, it is shown that even noisy time series allow the identification of the stability parameters at remarkable accuracy. This method will be useful to obtain the DNA stacking and hydrogen bonding free energies from single bubble breathing assays rather than equilibrium data.

Target site recognition by a diversity-generating retroelement

Diversity-generating retroelements (DGRs) are in vivo sequence diversification machines that are widely distributed in bacterial, phage, and plasmid genomes. They function to introduce vast amounts of targeted diversity into protein-encoding DNA sequences via mutagenic homing. Adenine residues are converted to random nucleotides in a retrotransposition process from a donor template repeat (TR) to a recipient variable repeat (VR). Using the Bordetella bacteriophage BPP-1 element as a prototype, we have characterized requirements for DGR target site function. Although sequences upstream of VR are dispensable, a 24 bp sequence immediately downstream of VR, which contains short inverted repeats, is required for efficient retrohoming. The inverted repeats form a hairpin or cruciform structure and mutational analysis demonstrated that, while the structure of the stem is important, its sequence can vary. In contrast, the loop has a sequence-dependent function. Structure-specific nuclease digestion confirmed the existence of a DNA hairpin/cruciform, and marker coconversion assays demonstrated that it influences the efficiency, but not the site of cDNA integration. Comparisons with other phage DGRs suggested that similar structures are a conserved feature of target sequences. Using a kanamycin resistance determinant as a reporter, we found that transplantation of the IMH and hairpin/cruciform-forming region was sufficient to target the DGR diversification machinery to a heterologous gene. In addition to furthering our understanding of DGR retrohoming, our results suggest that DGRs may provide unique tools for directed protein evolution via in vivo DNA diversification.

Diversity-generating retroelements function through a unique, reverse transcriptase–mediated “copy and replace” mechanism that enables repeated rounds of protein diversification, selection, and optimization. The ability of DGRs to introduce targeted diversity into protein-coding DNA sequences has the potential to dramatically accelerate the evolution of adaptive traits. The utility of these elements in nature is underscored by their widespread distribution throughout the bacterial domain. Here we define DNA sequences and structures that are necessary and sufficient to direct the diversification machinery to specified target sequences. In addition to providing mechanistic insights into conserved features of DGR activity, our results provide a blueprint for the use of DGRs for a broad range of protein engineering applications.

Tandem GA residues on opposite sides of the loop in molecular beacon-like DNA hairpins compact the loop and increase hairpin stability

The free solution electrophoretic mobilities and thermal stabilities of hairpins formed by two complementary 26-nucleotide oligomers have been measured by capillary electrophoresis. The oligomers are predicted to form molecular beacon-like hairpins with 5 bp stems and 16 nucleotides in the loop. One hairpin, called hairpin2 (hp2), migrates with a relatively fast free solution mobility and exhibits melting temperatures that are reasonably well predicted by the popular structure-prediction program Mfold. Its complement, called hairpin1 (hp1), migrates with a slower free solution mobility and forms a stable hairpin only in solutions containing ≥200 mM Na(+). The melting temperatures observed for hp1 are ~18 °C lower than those observed for hp2 and ~20 °C lower than those predicted by Mfold. The greater thermal stability of hp2 is due to the presence of tandem GA residues on opposite sides of the loop. If the corresponding TC residues in the hp1 loop are replaced by tandem GA residues, the melting temperatures of the modified hairpin are close to those observed for hp2. Eliminating the tandem GA residues in the hp2 loop significantly decreases the thermal stability of hp2. If the loops are replaced by a loop of 16 thymine residues, the free solution mobilities and thermal stabilities of the T-loop hairpin are equal to those observed for hp1. Hence, the loop of hp1 appears to be relatively unstructured, with few base-base stacking interactions. Interactions between tandem GA residues on opposite sides of the hp2 loop appear to compact the loop and increase hairpin stability.

Influence of a hairpin loop on the thermodynamic stability of a DNA oligomer

DSC was used to evaluate the mechanism of the thermally induced unfolding of the single-stranded hairpin HP = 5′-CGGAATTCCGTCTCCGGAATTCCG-3′ and its core duplex D (5′-CGGAATTCCG-3′)₂. The DSC melting experiments performed at several salt concentrations were successfully described for HP and D in terms of a three-state transition model HP↔I (intermediate state) ↔ S (unfolded single-stranded state) and two state transition model D↔2S, respectively. Comparison of the model-based thermodynamic parameters obtained for each HP and D transition shows that in unfolding of HP only the HP↔I transition is affected by the TCTC loop. This observation suggests that in the intermediate state its TCTC loop part exhibits significantly more flexible structure than in the folded state while its duplex part remains pretty much unchanged.

Non-B conformations of CAG repeats using 2-aminopurine

Repetition of trinucleotide sequences is the molecular basis of ~30 hereditary neurological and neurodegenerative diseases, and alternate structures adopted by these sequences are implicated in the etiology of such diseases. Elucidating these structures is important for advancing mechanistic understanding and ultimately treatment. Studies of (CAG) repeats are motivated by their involvement in a number of these diseases, and the structures favored by (CAG)₈ are discussed in this contribution. Utilizing the strong effect of base stacking on fluorescence quantum yield, 2-aminopurine is used in place of adenine to determine the secondary structures adopted by such repeated sequences. Alone, (CAG)₈ folds into a hairpin comprised of a duplex stem and a single-stranded loop. Energetic studies indicate that the hairpin is anchored by the interactions in the stem and has a strained loop environment. As a model for intermediates that form during repeat expansion, (CAG)₈ was also incorporated into a duplex to form a three-way junction. In contrast to the isolated (CAG)₈, this integrated repeat adopts an open, unfolded loop. Enthalpy and entropy changes associated with denaturation indicate that the stability of the three-way junction is dominated by interactions in the duplex arms and that the repeated sequence tracks global unfolding. Because 2-aminopurine provides both structural and energetic information via fluorescence and also is an innocuous substitution for adenine, significant progress in elucidating the secondary structures of (CAG) repeats will be achieved.

Structural studies of a trinucleotide repeat sequence using 2-aminopurine

The secondary structure of repeated trinucleotide sequences results in the development of several neurodegenerative diseases, and these studies consider the (CAG)(8) sequence that forms a stem-loop hairpin. The structural and thermodynamic properties of this hairpin are assessed using 2-aminopurine substitutions for adenine at six positions in this repeated sequence. Circular dichroism spectra and thermal denaturation experiments show that the secondary structure is not disturbed by the modifications. The local structure of the hairpin was monitored using the fluorescence intensities of 2-aminopurines, the changes in the intensity relative to the denatured state, and the sensitivity of the fluorescence to quenching by acrylamide. To establish the stem and loop characteristics in (CAG)(8), known reference points for stem, loop, and exposed base motifs were used. In the vicinity of the loop, the bases become more solvent exposed, which suggests that the instability associated with this repeated hairpin influences the global secondary structure. These results provide the basis to interpret the structures adopted by other repeated (CAG) structures.

Thermodynamics of single strand DNA base stacking

The thermodynamics of the stacking to unstacking transitions of 24 single-stranded DNA sequences (ssDNA), 10-12 bases in length, in sodium phosphate buffer were determined from 10 to 95 degrees C, using differential scanning calorimetry (DSC). An additional 22 ssDNA sequences did not exhibit an S<=>U transition in this temperature range. The transition properties of the ssDNA sequences with <or=60% self-complementarily in the reverse direction were independent of concentration with transition temperatures ranging from 15 to 70 degrees C, van't Hoff transition enthalpies from 92 to 201 kJ mol(-1) and transition enthalpies from 5 to 75% of the corresponding van't Hoff transition enthalpies. Since all the 16 doublets and 60 of the 64 triplets are present in both the transition and the non-transition ssDNA sequences, it is unlikely that the nucleation subset initiating stacking of the sequence is a specific doublet or triplet subset. Of the 141 quadruplet subsets of the 46 sequences, each transition ssDNA sequence contained at least one or more quadruplets not found in the non-transition ssDNA sequences. It could be concluded that the thermal stability of the stacked conformation was dependent on the presence of a possible nucleation quadruplet and the length of the ssDNA sequence and not on the G or C content of the ssDNA sequence, nor on the number of purine bases in the sequence.

An interlocked dimer of the protelomerase TelK distorts DNA structure for the formation of hairpin telomeres

The termini of linear chromosomes are protected by specialized DNA structures known as telomeres that also facilitate the complete replication of DNA ends. The simplest type of telomere is a covalently closed DNA hairpin structure found in linear chromosomes of prokaryotes and viruses. Bidirectional replication of a chromosome with hairpin telomeres produces a catenated circular dimer that is subsequently resolved into unit-length chromosomes by a dedicated DNA cleavage-rejoining enzyme known as a hairpin telomere resolvase (protelomerase). Here we report a crystal structure of the protelomerase TelK from Klebsiella oxytoca phage varphiKO2, in complex with the palindromic target DNA. The structure shows the TelK dimer destabilizes base pairing interactions to promote the refolding of cleaved DNA ends into two hairpin ends. We propose that the hairpinning reaction is made effectively irreversible by a unique protein-induced distortion of the DNA substrate that prevents religation of the cleaved DNA substrate.

The energetics of i-DNA tetraplex structures formed intermolecularly by d(TC5) and intramolecularly by d[(C5T3)3C5]

Cytosine-rich DNA at low pH adopts an antiparallel tetraplex structure via the intercalation of two partially protonated, parallel stranded duplexes. This intriguing structural motif has been named i-DNA. We have used a combination of spectroscopic and calorimetric techniques to characterize the properties of an intermolecular i-DNA formed by d(TC(5)) and an intramolecular i-DNA formed by d[(C(5)T(3))(3)C(5)]. Our measurements reveal that both i-DNA complexes are enthalpically stabilized by 6.5-7.0 kcal/mol(base) and entropically destabilized by 20 cal/mol(base)/K. These values are about 50% larger than the corresponding enthalpy and entropy values per base for Watson and Crick duplexes and for Hoogsteen triplexes, while being similar to per base enthalpy and entropy values reported for G-quadruplexes. Our data also reveal a positive heat capacity change between 20 and 30 cal/mol(base)/K, values similar to that reported for polymeric Watson & Crick DNA duplexes. Solution-dependent studies reveal the overall thermal and thermodynamic stability of i-DNA complexes to be dictated by an interplay between pH and ionic strength. Based on the thermodynamic data measured, we discuss the feasibility of i-DNA formation in the context of conventional DNA sequences, while commenting on potential roles for this structural motif in biological regulatory mechanisms.

Capillary electrophoresis is a sensitive monitor of the hairpin-random coil transition in DNA oligomers

Capillary electrophoresis (CE) has been used to characterize the hairpin-random coil transition of four octamers in the GCxxxxGC minihairpin family, where xxxx is GAAA, TTTC, TTTT, or AAAA. The transition can be monitored by CE because differences in the frictional coefficients of the hairpin and coil forms of each octamer lead to a difference of approximately 9% in the free solution mobilities of the two conformations. The GAAA octamer is unusually stable, with a melting temperature of 65 degrees C. The TTTT octamer forms a minihairpin with a melting temperature of 29 degrees C, the TTTC octamer has a melting temperature of 16 degrees C, and the AAAA octamer has a melting temperature below 0 degrees C. The thermal transitions of the TTTT, TTTC, and AAAA octamers are well fitted by a structure prediction algorithm; however, the GAAA minihairpin is considerably more stable than predicted. The melting temperature of the GAAA minihairpin is reduced to 47 degrees C in aqueous buffers containing 7.2M urea and to 33 degrees C in buffers containing 7.2M urea plus 40% (v/v) formamide. The combined results indicate that CE is a sensitive technique for monitoring conformational transitions in small DNA molecules.

Increase in hybridization rates with oligodeoxyribonucleotides containing locked nucleic acids

Locked nucleic acids (LNAs) incorporated into either stable single stranded oligonucleotides containing tetraloops or their complements have been found to increase second order hybridization rate constants by an order of magnitude compared to the all-DNA hybridization rate constants. Model sequences composed of 20 bases in length that can form hairpins due to a stable GAAA tetraloop were used where LNAs were substituted for the nucleotides in the loop, stem, or end regions of the strand and in the complementary strand. Substitution of the LNAs to the loop predictably raised the melting temperatures of the duplex however, the hybridization rates between the tetraloop and the complementary sequence also increased. In contrast, when LNAs were substituted in the stem, the hybridization rate decreased implying the formation of a more stable hairpin. Substitution of LNAs into the end region of the sequence had little effect on the hybridization rate constants although melting temperatures still showed a predictable increase. Rates also increased when LNAs were substituted into complementary strands of DNA tetraloops. The increase in hybridization rate constant is being attributed to changes in the structure of the stable single strands.

Sequence-dependent gating of an ion channel by DNA hairpin molecules

DNA hairpins produce ionic current signatures when captured by the alpha-hemolysin nano-scale pore under conditions of single molecule electrophoresis. Gating patterns produced by individual DNA hairpins when captured can be used to distinguish differences of a single base pair or even a single nucleotide [Vercoutere,W.A. et al. (2003) Nucleic Acids Res., 31, 1311–1318]. Here we investigate the mechanism(s) that may account for the ionic current gating signatures. The ionic current resistance profile of conductance states produced by DNA hairpin molecules with 3–12 bp stems showed a plateau in resistance between 10 and 12 bp, suggesting that hairpins with 10–12 bp stems span the pore vestibule. DNA hairpins with 9–12 bp stems produced gating signatures with the same relative conductance states. Systematic comparison of the conductance state dwell times and apparent activation energies for a series of 9–10 bp DNA hairpins suggest that the 3′ and 5′ ends interact at or near the limiting aperture within the vestibule of the alpha-hemolysin pore. The model presented may be useful in predicting and interpreting DNA detection using nanopore detectors. In addition, this well-defined molecular system may prove useful for investigating models of ligand-gated channels in biological membranes.

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND

Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis.

RESULTS

HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.

A mass spectrometry-based approach for identifying novel DNA polymerase substrates from a pool of dNTP analogues

There has been a long-standing interest in the discovery of unnatural nucleotides that can be incorporated into DNA by polymerases. However, it is difficult to predict which nucleotide analogs will prove to have biological relevance. Therefore, we have developed a new screening method to identify novel substrates for DNA polymerases. This technique uses the polymerase itself to select a dNTP from a pool of potential substrates via incorporation onto a short oligonucleotide. The unnatural nucleotide(s) is then identified by high-resolution mass spectrometry. By using a DNA polymerase as a selection tool, only the biologically relevant members of a small nucleotide library can be quickly determined. We have demonstrated that this method can be used to discover unnatural base pairs in DNA with a detection threshold of ≤10% incorporation.