Thermodynamics of DNA hairpins: contribution of loop size to hairpin stability and ethidium binding

Amino Acid Stabilization of Nucleic Acid Secondary Structure: Kinetic Insights from Single-Molecule Studies

Amino acid and nucleic acid interactions are central in biology and may have played a role in the evolutionary development of protein-based life from an early "RNA Universe." To explore the possible role of single amino acids in promoting nucleic acid folding, single-molecule Förster resonance energy transfer experiments have been implemented with a DNA hairpin construct (7 nucleotide double strand with a 40A loop) as a simple model for secondary structure formation. Exposure to positively charged amino acids (arginine and lysine) is found to clearly stabilize the secondary structure. Kinetically, each amino acid promotes folding by generating a large increase in the folding rate with little change in the unfolding rate. From analysis as a function of temperature, arginine and lysine are found to significantly increase the overall exothermicity of folding while imposing only a small entropic penalty on the folding process. Detailed investigations into the kinetics and thermodynamics of this amino acid-induced folding stability reveal arginine and lysine to interact with nucleic acids in a manner reminiscent of monovalent cations. Specifically, these observations are interpreted in the context of an ion atmosphere surrounding the nucleic acid, in which amino acid salts stabilize folding qualitatively like small monovalent cations but also exhibit differences because of the composition of their side chains.

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

Rapid detection of exosomal microRNA biomarkers by electrokinetic concentration for liquid biopsy on chip

An ion concentration polarization (ICP)-based electrokinetic concentration device is used for accelerating the surface hybridization reaction between exosomal microRNAs (miRNAs) and morpholinos (MOs) as a synthetic oligo capture probe in the nanomolar concentration range in a microfluidic channel. Compared with standard hybridization at the same concentration, the hybridization time of the miRNA target on MO capture probes could be reduced from ∼24 h to 30 min, with an increase in detection speed by 48 times. This ICP-enhanced hybridization method not only significantly decreases the detection time but also makes workflow simple to use, circumventing use of quantitative reverse transcription polymerase chain reaction or other conventional enzyme-based amplification methods that can cause artifacts.

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.

Investigation of the Melting Behavior of DNA Three-Way Junctions in the Closed and Open States

Intramolecular three-way junctions are commonly found in both DNA and RNA. These structures are functionally relevant in ribozymes, riboswitches, rRNA, and during replication. In this work, we present a thermodynamic description of the unfolding of DNA intramolecular three-way junctions. We used a combination of spectroscopic and calorimetric techniques to investigate the folding/unfolding thermodynamics of two three-way junctions with a closed (Closed-J) or open (Open-J) junction and their appropriate control stem-loop motifs (GAAATT-Hp, CTATC-Hp, and Dumbbell). The overall results show that both junctions are stable over a wide range of salt concentrations. However, Open-J is more stable due to a higher enthalpy contribution from the formation of a higher number of basepair stacks whereas Closed-J has a defined structure and retains the basepair stacking of all three stems. The comparison of the experimental results of Closed-J and Open-J with those of their component stem-loop motifs allowed us to be more specific about their cooperative unfolding. For instance, Closed-J sacrifices thermal stability of the Dumbbell structure to maintain an overall folded state. At higher salt concentration, the simultaneous unfolding of the above domains is lost, resulting in the unfolding of the three separate stems. In contrast, the junction of Open-J in low salt retains the thermal and enthalpic stability of the Dumbbell structure although sacrificing stability of the CTATC stem. The relative stability of Dumbbell is the primary reason for the higher ΔG°₍₅₎, or free energy, value seen for Open-J at low salt. Higher salt not only maintains thermal stability of the Dumbbell structure in Open-J but causes the CTATC stem to fully fold.

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.

Biofunction-assisted DNA detection through RNase H-enhanced 3' processing of a premature tRNA probe in a wheat germ extract

We have developed a novel type of biofunction-assisted, signal-turn-on sensor for simply and homogenously detecting DNA. This sensor system is composed of two types of in vitro-transcribed label-free RNAs (a 3' premature amber suppressor tRNA probe and an amber-mutated mRNA encoding a reporter protein), RNase H, and a wheat germ extract (WGE). A target DNA induces the 3' end maturation of the tRNA probe, which is enhanced by RNase H and leads to the expression of a full-length reporter protein through amber suppression in WGE, while there is almost no expression without the target due to the inactivity of the premature probe. Therefore, the target can be readily detected with the activity of the translated reporter. The catalytic reuse of the target with the help of RNase H in addition to various bioprocesses in WGE enables this sensor system to exhibit relatively high selectivity and sensitivity.

DNA based arithmetic function: a half adder based on DNA strand displacement

Biomolecular programming utilizes the reactions and information stored in biological molecules, such as proteins and nucleic acids, for computational purposes. DNA has proven itself an excellent candidate for building logic operating systems due to its highly predictable molecular behavior. In this work we designed and realized an XOR logic gate and an AND logic gate based on DNA strand displacement reactions. These logic gates utilize ssDNA as input and output signals. The XOR gate and the AND gate were used as building blocks for constructing a half adder logic circuit, which is a primary step in constructing a full adder, a basic arithmetic unit in computing. This work provides the field of DNA molecular programming with a potential universal arithmetic tool.

Effect of interior loop length on the thermal stability and pK(a) of i-motif DNA

The four-stranded i-motif (iM) conformation of cytosine-rich DNA is important in a wide variety of biochemical systems ranging from its use in nanomaterials to a potential role in oncogene regulation. An iM is stabilized by acidic pH that allows hemiprotonated cytidines to form a C·C(+) base pair. Fundamental studies that aim to understand how the lengths of loops connecting the protonated C·C(+) pairs affect intramolecular iM physical properties are described here. We characterized both the thermal stability and the pK(a) of intramolecular iMs with differing loop lengths, in both dilute solutions and solutions containing molecular crowding agents. Our results showed that intramolecular iMs with longer central loops form at pHs and temperatures higher than those of iMs with longer outer loops. Our studies also showed that increases in thermal stability of iMs when molecular crowding agents are present are dependent on the loop that is lengthened. However, the increase in pK(a) for iMs when molecular crowding agents are present is insensitive to loop length. Importantly, we also determined the proton activity of solutions containing high concentrations of molecular crowding agents to ascertain whether the increase in pK(a) of an iM is caused by alteration of this activity in buffered solutions. We determined that crowding agents alone increase the apparent pK(a) of a number of small molecules as well as iMs but that increases to iM pK(a) were greater than that expected from a shift in proton activity.

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.

Identification and characterization of human apurinic/apyrimidinic endonuclease-1 inhibitors

The repair of abasic sites that arise in DNA from hydrolytic depurination/depyrimidination of the nitrogenous bases from the sugar-phosphate backbone and the action of DNA glycosylases on deaminated, oxidized, and alkylated bases are critical to cell survival. Apurinic/apyrimidinic endonuclease-1/redox effector factor-1 (APE-1; aka APE1/ref-1) is responsible for the initial removal of abasic lesions as part of the base excision repair pathway. Deletion of APE-1 activity is embryonic lethal in animals and is lethal in cells. Potential inhibitors of the repair function of APE-1 were identified based upon molecular modeling of the crystal structure of the APE-1 protein. We describe the characterization of several unique nanomolar inhibitors using two complementary biochemical screens. The most active molecules all contain a 2-methyl-4-amino-6,7-dioxolo-quinoline structure that is predicted from the modeling to anchor the compounds in the endonuclease site of the protein. The mechanism of action of the selected compounds was probed by fluorescence and competition studies, which indicate, in a specific case, direct interaction between the inhibitor and the active site of the protein. It is demonstrated that the inhibitors induce time-dependent increases in the accumulation of abasic sites in cells at levels that correlate with their potency to inhibit APE-1 endonuclease excision. The inhibitor molecules also potentiate by 5-fold the toxicity of a DNA methylating agent that creates abasic sites. The molecules represent a new class of APE-1 inhibitors that can be used to probe the biology of this critical enzyme and to sensitize resistant tumor cells to the cytotoxicity of clinically used DNA damaging anticancer drugs.

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.

Hexamer oligonucleotide topology and assembly under solution phase NMR and theoretical modeling scrutiny

The entire family of noncomplementary hexamer oligodeoxyribonucleotides d(GCXYGC) (X and Y = A, G, C, or T) were assessed for topological indicators and equilibrium thermodynamics using a priori molecular modeling and solution phase NMR spectroscopy. Feasible modeled hairpin structures formed a basis from which solution structure and equilibria for each oligonucleotide were considered. ¹H and ³¹P variable temperature-dependent (VT) and concentration-dependent NMR data, NMR signal assignments, and diffusion parameters led to d(GCGAGC) and d(GCGGGC) being understood as exceptions within the family in terms of self-association and topological character. A mean diffusion coefficient D(298 K) = (2.0 ± 0.07) × 10⁻¹⁰ m² s⁻¹ was evaluated across all hexamers except for d(GCGAGC) (D(298 K) = 1.7 × 10⁻¹⁰ m² s⁻¹) and d(GCGGGC) (D(298 K) = 1.2 × 10⁻¹⁰ m² s⁻¹). Melting under VT analysis (T(m) = 323 K) combined with supporting NMR evidence confirmed d(GCGAGC) as the shortest tandem sheared GA mismatched duplex. Diffusion measurements were used to conclude that d(GCGGGC) preferentially exists as the shortest stable quadruplex structure. Thermodynamic analysis of all data led to the assertion that, with the exception of XY = GA and GG, the remaining noncomplementary oligonucleotides adopt equilibria between monomer and duplex, contributed largely by monomer random-coil forms. Contrastingly, d(GCGAGC) showed preference for tandem sheared GA mismatch duplex formation with an association constant K = 3.9 × 10⁵M⁻¹. No direct evidence was acquired for hairpin formation in any instance although its potential existence is considered possible for d(GCGAGC) on the basis of molecular modeling studies.

RNAstructure: software for RNA secondary structure prediction and analysis

Background

To understand an RNA sequence's mechanism of action, the structure must be known. Furthermore, target RNA structure is an important consideration in the design of small interfering RNAs and antisense DNA oligonucleotides. RNA secondary structure prediction, using thermodynamics, can be used to develop hypotheses about the structure of an RNA sequence.

Results

RNAstructure is a software package for RNA secondary structure prediction and analysis. It uses thermodynamics and utilizes the most recent set of nearest neighbor parameters from the Turner group. It includes methods for secondary structure prediction (using several algorithms), prediction of base pair probabilities, bimolecular structure prediction, and prediction of a structure common to two sequences. This contribution describes new extensions to the package, including a library of C++ classes for incorporation into other programs, a user-friendly graphical user interface written in JAVA, and new Unix-style text interfaces. The original graphical user interface for Microsoft Windows is still maintained.

Conclusion

The extensions to RNAstructure serve to make RNA secondary structure prediction user-friendly. The package is available for download from the Mathews lab homepage at http://rna.urmc.rochester.edu/RNAstructure.html.

Binding the mammalian high mobility group protein AT-hook 2 to AT-rich deoxyoligonucleotides: enthalpy-entropy compensation

HMGA2 is a DNA minor-groove binding protein. We previously demonstrated that HMGA2 binds to AT-rich DNA with very high binding affinity where the binding of HMGA2 to poly(dA-dT)(2) is enthalpy-driven and to poly(dA)poly(dT) is entropy-driven. This is a typical example of enthalpy-entropy compensation. To further study enthalpy-entropy compensation of HMGA2, we used isothermal-titration-calorimetry to examine the interactions of HMGA2 with two AT-rich DNA hairpins: 5'-CCAAAAAAAAAAAAAAAGCCCCCGCTTTTTTTTTTTTTTTGG-3' (FL-AT-1) and 5'-CCATATATATATATATAGCCCCCGCTATATATATATATATGG-3' (FL-AT-2). Surprisingly, we observed an atypical isothermal-titration-calorimetry-binding curve at low-salt aqueous solutions whereby the apparent binding-enthalpy decreased dramatically as the titration approached the end. This unusual behavior can be attributed to the DNA-annealing coupled to the ligand DNA-binding and is eliminated by increasing the salt concentration to approximately 200 mM. At this condition, HMGA2 binding to FL-AT-1 is entropy-driven and to FL-AT-2 is enthalpy-driven. Interestingly, the DNA-binding free energies for HMGA2 binding to both hairpins are almost temperature independent; however, the enthalpy-entropy changes are dependent on temperature, which is another aspect of enthalpy-entropy compensation. The heat capacity change for HMGA2 binding to FL-AT-1 and FL-AT-2 are almost identical, indicating that the solvent displacement and charge-charge interaction in the coupled folding/binding processes for both binding reactions are similar.

Monovalent ion dependence of neomycin B binding to an RNA aptamer characterized by spectroscopic methods

Understanding the interactions of small molecules like antibiotics with RNA is a prerequisite for the development of novel drugs. In this study we address structural and thermodynamic features of such interactions by using a simple model system: the binding of the highly charged antibiotic neomycin B to a short hairpin RNA molecule. Nucleotide A16, which acts as a flap over the neomycin B binding pocket, was substituted by the fluorescent adenine analogue 2-aminopurine (2-AP). Steady-state and time-resolved fluorescence measurements were complemented by UV-melting and circular dichroism studies. The binding of neomycin B at three sites was found to have a strong inverse correlation with Na(+) concentration. For the highest-affinity site, both fluorescence and UV absorption experiments were consistent with a model assuming at least three neomycin NH(3) (+) groups participating in addition to hydrogen bonds in electrostatic interactions with the RNA. The variation of fluorescence intensity and lifetime upon neomycin B binding indicated unstacking of 2-AP16 from neighbouring bases as it flipped over the binding pocket. RNA conformational changes upon binding of the antibiotic were confirmed by circular dichroism. The two weaker binding sites were characterized as unspecific binding to the aptamer, while the high-affinity binding event was shown to be highly specific even at high ionic concentration. In addition, 2-AP was confirmed to be a noninvasive fluorescent probe; it serves as a sensitive spectroscopic tool to investigate details of the interactions between small molecules and RNA.

Unfolding of G-quadruplexes: energetic, and ion and water contributions of G-quartet stacking

It has been shown that DNA oligonucleotides composed, in part, of G repeat sequences can adopt G-quadruplex structures in the presence of specific metal ions. In this work, we use a combination of spectroscopic and calorimetric techniques to determine the spectral and thermodynamic characteristics of two DNA aptamers, d(G2T2G2TGTG2T2G2), G2, and d(G3T2G3TGTG3T2G3), G3; a sequence in the promoter region of the c-MYC oncogene, d(TG4AG3TG4AG3TG4A2G2), NHE-III; and the human telomere sequence d(AG3T2AG3T2AG3T2AG3), 22GG. The circular dichroism spectra of these oligonucleotides in the presence of K+ indicate that all form G-quadruplexes with G-quartets in an antiparallel arrangement (G2), in a parallel arrangement (NHE-III and 22GG), or in a mixed parallel and antiparallel G-quartet arrangement (G3). Melting profiles show transition temperatures, TM, above 45 degrees C that are independent of strand concentration, consistent with the formation of very stable intramolecular G-quadruplexes. We used differential scanning calorimetry to obtain complete thermodynamic profiles for the unfolding of each quadruplex. Subtracting the thermodynamic folding profiles of G2 from those of G3 yielded the following thermodynamic profile for the formation of a G-quartet stack: DeltaG degrees 20 = -2.2 kcal/mol, DeltaHcal = -14.6 kcal/mol, TDeltaScal = -12.4 kcal/mol, DeltanK+ = -0.3 mol of K+/mol, and DeltanW = 13 mol of H2O/mol. Furthermore, we used this profile to estimate the thermodynamic contributions of the loops and/or extra base sequences of each oligonucleotide in the G-quadruplex state. The average free energy contributions of the latter indicate that the incorporation of loops and base overhangs stabilizes quadruplex structures. This stabilization is enthalpy-driven and is due to base-stacking contributions.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30δ DNA

Upon increasing the concentration of ethylene glycol (EG) at 37 degrees C, the twist energy parameter, E(T), which governs the supercoiling free energy, was recently found to undergo a decreasing (or reverse) sigmoidal transition with a midpoint near 20 w/v % EG. In this study, the effects of adding 20 w/v % EG on the torsion elastic constant (alpha) of linear p30delta DNA and on the hydrodynamic radius (R(H)) of a synthetic 24 bp duplex DNA were examined at both 40 and 20 degrees C. The time-resolved fluorescence intensity and fluorescence polarization anisotropy (FPA) of intercalated ethidium were measured in order to assess the effects of 20 w/v % EG on: (1) alpha; (2) R(H); (3) the lifetimes of intercalated and non-intercalated dye; (4) the amplitude of dye wobble in its binding site; and (5) the binding constant for intercalation. The effects of 20 w/v % EG on the circular dichroism (CD) spectrum of the DNA and on the emission spectrum of the free dye were also measured. At 40 degrees C, addition of 20 w/v % EG caused a substantial (1.27- to 1.35-fold) increase in alpha, a significant change in the CD spectrum, and a very small, marginally significant increase in R(H), but little or no change in the amplitude of dye wobble in its binding site or the lifetime of intercalated dye. Together with previously reported measurements of E(T), these results imply that the bending elastic constant of DNA is significantly decreased by 20 w/v % EG at 40 degrees C. At 20 degrees C, addition of 20 w/v % EG caused a marginally significant decrease in alpha and very little change in any other measured properties. Also at 20 degrees C, addition of 30 w/v % betaine caused a marginally significant increase in alpha and significant but modest change in the CD spectrum, but very little change in any other properties.

The initial step of DNA hairpin folding: a kinetic analysis using fluorescence correlation spectroscopy

Conformational fluctuations of single-stranded DNA (ssDNA) oligonucleotides were studied in aqueous solution by monitoring contact-induced fluorescence quenching of the oxazine fluorophore MR121 by intrinsic guanosine residues (dG). We applied fluorescence correlation spectroscopy as well as steady-state and time-resolved fluorescence spectroscopy to analyze kinetics of DNA hairpin folding. We first characterized the reporter system by investigating bimolecular quenching interactions between MR121 and guanosine monophosphate in aqueous solution estimating rate constants, efficiency and stability for formation of quenched complexes. We then studied the kinetics of complex formation between MR121 and dG residues site-specifically incorporated in DNA hairpins. To uncover the initial steps of DNA hairpin folding we investigated complex formation in ssDNA carrying one or two complementary base pairs (dC–dG pairs) that could hybridize to form a short stem. Our data show that incorporation of a single dC–dG pair leads to non-exponential decays for opening and closing kinetics and reduces rate constants by one to two orders of magnitude. We found positive activation enthalpies independent of the number of dC–dG pairs. These results imply that the rate limiting step of DNA hairpin folding is not determined by loop dynamics, or by mismatches in the stem, but rather by interactions between stem and loop nucleotides.

Determination of the reopening temperature of a DNA hairpin structure in vitro

A novel method, based upon primer extension, has been developed for measuring the reopening temperature of a single type of DNA hairpin structure. Two DNA oligonucleotides have been utilized and designated as primers 1 and 2. Primer 1, with its 5- and 3'-termini fully complementary to the hairpin flanking sequences, was used to evaluate primer extension conditions, and primer 2, with its 3'-end competing with the DNA hairpin stem, was used to detect the DNA hairpin reopening temperature. A single DNA hairpin structure was formed on the DNA template by thermal denaturation and renaturation, and this hairpin structure was predicted to prevent the annealing of the 3'-end of primer 2 with the template DNA, which leads to no primer extension. By incubating at different temperatures, the DNA hairpin structure can be reopened at a particular temperature where the primer extension can be carried out. This resulted in the appearance of double-stranded DNA that was detected on an agarose gel. This temperature is defined here as the hairpin reopening temperature.

Protonation of deoxycytidine residues in dC4 tetraloops: UV spectrophotometric study of dC10 and d(A14C4T14)

It is shown that component analysis could be applied to study the UV difference spectra of cytidine oligomers and hairpin oligonucleotides with cytidines in the loop region in order to account for the melting and titration results in terms of cytidine stacking and protonation. Upon acid titration, the dC(10) oligomer undergoes cooperative conformational transition at pH 6.3 accompanied by protonation and formation of the i-structure with half of the residues protonated. The stability of the hemiprotonated structure increases with decreasing pH, the i-structure persisting still in the region of pH<pK of cytidine. An UV difference spectrum that reflects the stacking/unstacking of hemiprotonated cytidine residues was acquired from the melting and titration experiments of the dC(10) oligomer and used to describe the behavior of the dC(4) loop of the hairpin oligonucleotide d(A(14)C(4)T(14)). It is shown that upon titration, the 50% level of protonation of the deoxycytidine tetraloop is attained at pH 5.0. Simultaneously, the stacking interactions of cytidine residues reach the maximum at this pH with two residues stacked, and thereafter decline again. Only marginal stabilization of the oligomer hairpin (DeltaT(m)=1.5 degrees C) is found to accompany the formation of this single hemiprotonated dC.dC(+) base pair. We propose that at pH 5 the cytidines of the dC(4) loop form a hemiprotonated dC.dC(+) pair stacked with the last dA.dT base pair of the hairpin stem.

Determination of melting temperature for variant detection using dHPLC: a comparison between an empirical approach and DNA melting prediction software

Detection of DNA sequence variants by the use of denaturing high-performance liquid chromatography (dHPLC) is a relatively new method (Underhill et al., 1997) and has distinct advantages over other methods such as single-strand conformation polymorphism (SSCP), direct sequencing, and DNA chip hybridization. The dHPLC-based single-nucleotide polymorphism (SNP) screening relies on different DNA thermodynamic properties between perfectly matched base pairs in homoduplex molecules and single base-pair mismatches in heteroduplex DNAs. Separation of the two forms of duplex DNAs by dHPLC is based on ionic forces between the negatively charged DNA and the hydrophobic stationary phase, which consists of C(18) chains on PSDVB (polystyrene-divinylbenzene) beads coated with a positively charged ion-pairing agent (TEAA, triethylammonium acetate). Removal of the DNA from the TEAA-coated beads is dependent upon a mobile organic phase, in the form of a linear acetonitrile gradient. The major factor that influences the success of dHPLC to detect sequence variation is the thermal stability of the duplex DNA, which is determined by the melting temperature (TM(50)), where 50% of the DNA strand is single stranded and 50% is double stranded. The TM(50) predicts the best probability of detecting a single base-pair change based on the altered thermodynamics it imparts to the DNA duplex. Generally, there are two ways to determine this melting temperature, either empirically or with the aid of predictive DNA melting analysis software. Such programs include the DNAMelt program located on the Stanford University DNA Sequencing and Technology Center website, MeltCalc (Schutz and von Ahsen 1999), and WAVEMAKER, the proprietary melting analysis software provided with the Transgenomic WAVE dHPLC system. The goal of the current study was to determine whether currently available predictive DNA melting programs could be used to increase efficiency and throughput of SNP detection. A wide range of amplicons, differing in both size and GC composition, were selected for analysis to simulate the broad spectrum of PCR products that may be encountered during a large-scale dHPLC screening project.

Characterization of nucleic acids by nanopore analysis

Single-stranded DNA and RNA molecules in solution can be driven through a nanoscopic pore by an applied electric field. As each molecule occupies the pore, a characteristic blockade of ionic current is produced. Information about length, composition, structure, and dynamic motion of the molecule can be deduced from modulations of the current blockade.

Effect of secondary structure on the thermodynamics and kinetics of PNA hybridization to DNA hairpins

The binding of a series of PNA and DNA probes to a group of unusually stable DNA hairpins of the tetraloop motif has been observed using absorbance hypochromicity (ABS), circular dichroism (CD), and a colorimetric assay for PNA/DNA duplex detection. These results indicate that both stable PNA-DNA and DNA-DNA duplexes can be formed with these target hairpins, even when the melting temperatures for the resulting duplexes are up to 50 degrees C lower than that of the hairpin target. Both hairpin/single-stranded and hairpin/hairpin interactions are considered in the scope of these studies. Secondary structures in both target and probe molecules are shown to depress the melting temperatures and free energies of the probe-target duplexes. Kinetic analysis of hybridization yields reaction rates that are up to 160-fold slower than hybridization between two unstructured strands. The thermodynamic and kinetic obstacles to hybridization imposed by both target and probe secondary structure are significant concerns for the continued development of antisense agents and especially diagnostic probes.

A semiflexible polymer model applied to loop formation in DNA hairpins

A statistical mechanical "zipper" model is applied to describe the equilibrium melting of short DNA hairpins with poly(dT) loops ranging from 4 to 12 bases in the loop. The free energy of loop formation is expressed in terms of the persistence length of the chain. This method provides a new measurement of the persistence length of single-stranded DNA, which is found to be approximately 1.4 nm for poly(dT) strands in 100 mM NaCl. The free energy of the hairpin relative to the random coil state is found to scale with the loop size with an apparent exponent of > or = 7, much larger than the exponent of approximately 1.5-1.8 expected from considerations of loop entropy alone. This result indicates a strong dependence of the excess stability of the hairpins, from stacking interactions of the bases within the loop, on the size of the loop. We interpret this excess stability as arising from favorable hydrophobic interactions among the bases in tight loops and which diminish as the loops get larger. Free energy profiles along a generalized reaction coordinate are calculated from the equilibrium zipper model. The transition state for hairpin formation is identified as an ensemble of looped conformations with one basepair closing the loop, and with a lower enthalpy than the random coil state. The equilibrium model predicts apparent activation energy of approximately -11 kcal/mol for the hairpin closing step, in remarkable agreement with the value obtained from kinetics measurements.

Binding of ethidium to DNA measured using a 2D diffusion-modulated gradient COSY NMR experiment

The binding of ethidium bromide to a DNA hairpin (dU(5)-hairpin) was investigated using a novel 2D diffusion-modulated gradient correlation spectroscopy (DMG-COSY) experiment to evaluate the applicability of this technique for studying the binding of drugs to DNA. The DMG-COSY experiment includes a preparation period during which coherent magnetization is attenuated due to molecular self-diffusion. Magnetization then evolves due to scalar coupling during an evolution delay, and is detected using gradient pulses for coherence selection. The time-domain data are processed in an analogous manner as for gradient-selected COSY experiments. The diffusion coefficient for uridine in DMSO solution was determined from the H5-H6 crosspeak intensities for a series of 2D DMG-COSY experiments that differed in the magnitude of the gradient pulses applied during the preparation period of the DMG-COSY experiment. The diffusion coefficient for uridine calculated from the DMG-COSY experiments was identical (within experimental error) to that determined from 1D diffusion experiments (5.24x10(-6) cm(2)/s at 26 degrees C). The diffusion coefficients for ethidium bromide and for the dU(5)-hairpin were first measured separately using the DMG-COSY experiment, and then measured in the putative complex. The diffusion coefficient for free ethidium bromide (4.15x10(-6) cm(2)/s at 26 degrees C) was considerably larger than for the dU(5)-hairpin (1. 60x10(-6) cm(2)/s at 26 degrees C), as expected for the smaller molecule. The diffusion coefficient for ethidium was markedly decreased upon addition of the dU(5)-hairpin, consistent with complex formation (1.22x10(-6) cm(2)/s at 26 degrees C). Complex formation of 1:1 stoichiometry between ethidium and the stem of the dU(5)-hairpin was verified independently by fluorescence spectroscopy. These results demonstrate the utility of the DMG-COSY experiment for investigating the binding of drugs to DNA in aqueous solution.