Predicting stability of DNA duplexes in solutions containing magnesium and monovalent cations

Sequence-dependent base pair stepping dynamics in XPD helicase unwinding

Helicases couple the chemical energy of ATP hydrolysis to directional translocation along nucleic acids and transient duplex separation. Understanding helicase mechanism requires that the basic physicochemical process of base pair separation be understood. This necessitates monitoring helicase activity directly, at high spatio-temporal resolution. Using optical tweezers with single base pair (bp) resolution, we analyzed DNA unwinding by XPD helicase, a Superfamily 2 (SF2) DNA helicase involved in DNA repair and transcription initiation. We show that monomeric XPD unwinds duplex DNA in 1-bp steps, yet exhibits frequent backsteps and undergoes conformational transitions manifested in 5-bp backward and forward steps. Quantifying the sequence dependence of XPD stepping dynamics with near base pair resolution, we provide the strongest and most direct evidence thus far that forward, single-base pair stepping of a helicase utilizes the spontaneous opening of the duplex. The proposed unwinding mechanism may be a universal feature of DNA helicases that move along DNA phosphodiester backbones. DOI:http://dx.doi.org/10.7554/eLife.00334.001.

The illusion of specific capture: surface and solution studies of suboptimal oligonucleotide hybridization

Background

Hybridization based assays and capture systems depend on the specificity of hybridization between a probe and its intended target. A common guideline in the construction of DNA microarrays, for instance, is that avoiding complementary stretches of more than 15 nucleic acids in a 50 or 60-mer probe will eliminate sequence specific cross-hybridization reactions. Here we present a study of the behavior of partially matched oligonucleotide pairs with complementary stretches starting well below this threshold complementarity length – in silico, in solution, and at the microarray surface. The modeled behavior of pairs of oligonucleotide probes and their targets suggests that even a complementary stretch of sequence 12 nt in length would give rise to specific cross-hybridization. We designed a set of binding partners to a 50-mer oligonucleotide containing complementary stretches from 6 nt to 21 nt in length.

Results

Solution melting experiments demonstrate that stable partial duplexes can form when only 12 bp of complementary sequence are present; surface hybridization experiments confirm that a signal close in magnitude to full-strength signal can be obtained from hybridization of a 12 bp duplex within a 50mer oligonucleotide.

Conclusions

Microarray and other molecular capture strategies that rely on a 15 nt lower complementarity bound for eliminating specific cross-hybridization may not be sufficiently conservative.

Magnesium-free self-assembly of multi-layer DNA objects

Molecular self-assembly with DNA offers a route for building user-defined nanoscale objects, but an apparent requirement for magnesium in solution has limited the range of conditions for which practical utility of such objects may be achieved. Here we report conditions for assembling templated multi-layer DNA objects in the presence of monovalent ions, showing that neither divalent cations in general or magnesium in particular are essential ingredients for the successful assembly of such objects. The percentage of DNA strands in an object that do not form thermally stable double-helical DNA domains (T_m>45 °C) with the template molecule correlated with the sodium requirements for obtaining folded objects. Minimizing the fraction of such weakly binding strands by rational design choices enhanced the yield of folding. The results support the view that DNA-based nanodevices may be designed and produced for a variety of target environments.

Self-assembly of DNA can provide access to a range of nanoscale structures, but assembly using magnesium has been considered essential. Martin and Dietz report conditions that allow the assembly of templated, multi-layer DNA structures in the presence of monovalent ions, rather than magnesium.

Screening interaction between ochratoxin A and aptamers by fluorescence anisotropy approach

By taking advantage of the intrinsic fluorescence of ochratoxin A (OTA), we present a fluorescence anisotropy approach for rapid analysis of the interactions between OTA and aptamers. The specific binding of OTA with a 36-mer aptamer can induce increased fluorescence anisotropy (FA) of OTA as the result of the freedom restriction of OTA and the increase of molecular volume, and the maximum FA change is about 0.160. This FA approach enables an easy way to investigate the effects of buffer compositions like metal ions on the affinity binding. FA analysis shows the interaction between OTA and aptamer is greatly enhanced by the simultaneous presence of Ca(2+) and Na(+), while the binding affinity of aptamer decreases more than 18-fold when only Ca(2+) exists, and the binding is completely lost when Ca(2+) is absent. Crucial region of the aptamer for binding can be mapped through FA analysis and aptamer mutation. The demonstrated FA approach maintains the advantages of FA in simplicity, rapidity, and robustness. This investigation will help the development of aptamer-based assays for OTA detection in optimizing the binding conditions, modification of aptamers, and rational design.

Y(III) interactions with guanine oligonucleotides covalently attached to aqueous/solid interfaces

The binding of Y(III) ions to surface-immobilized single-stranded 20-mers of guanine was studied using the Eisenthal χ((3)) technique and AFM. The free energy of binding for Y(III) to the G(20) sequence was found to be -39.5(8) kJ/mol. Furthermore, yttrium binds much more strongly to surface-immobilized oligonucleotides than the divalent metals previously reported. At maximum surface coverage, Y(III) ion densities range between one to three ions bound per strand. Comparatively, Mg(II) binds to the G(20)-functionalized interface in much higher ion densities. This result may be explained, in part, by the larger hydration sphere radius of Y(III) compared to that of Mg(II). The ion loading and binding free energy results, in conjunction with other surface and bulk aqueous phase studies, suggest that a fully hydrated +2 or +3 yttrium ion binds to the oligonucleotides through an outer-sphere mechanism. Tapping mode AFM results indicate that oligonucleotide height does not appreciably decrease following Y(III) binding. These results, together with the low ion densities for Y(III) ions, indicate that Y(III) strand loading may not significantly decrease the intrastrand Coulombic repulsions in order to cause a significant decrease in oligomer height.

Effect of thiazole orange doubly labeled thymidine on DNA duplex formation

Nucleic acid oligonucleotides are widely used in hybridization experiments for specific detection of complementary nucleic acid sequences. For design and application of oligonucleotides, an understanding of their thermodynamic properties is essential. Recently, exciton-controlled hybridization-sensitive fluorescent oligonucleotides (ECHOs) were developed as uniquely labeled DNA oligomers containing commonly one thymidine having two covalently linked thiazole orange dye moieties. The fluorescent signal of an ECHO is strictly hybridization-controlled, where the dye moieties have to intercalate into double-stranded DNA for signal generation. Here we analyzed the hybridization thermodynamics of ECHO/DNA duplexes, and thermodynamic parameters were obtained from melting curves of 64 ECHO/DNA duplexes measured by ultraviolet absorbance and fluorescence. Both methods demonstrated a substantial increase in duplex stability (ΔΔG°(37) ~ -2.6 ± 0.7 kcal mol(-1)) compared to that of DNA/DNA duplexes of the same sequence. With the exception of T·G mismatches, this increased stability was mostly unaffected by other mismatches in the position opposite the labeled nucleotide. A nearest neighbor model was constructed for predicting thermodynamic parameters for duplex stability. Evaluation of the nearest neighbor parameters by cross validation tests showed higher predictive reliability for the fluorescence-based than the absorbance-based parameters. Using our experimental data, a tool for predicting the thermodynamics of formation of ECHO/DNA duplexes was developed that is freely available at http://genome.gsc.riken.jp/echo/thermodynamics/. It provides reliable thermodynamic data for using the unique features of ECHOs in fluorescence-based experiments.

Ambient temperature detection of PCR amplicons with a novel sequence-specific nucleic acid lateral flow biosensor

In the field of diagnostics, molecular amplification targeting unique genetic signature sequences has been widely used for rapid identification of infectious agents, which significantly aids physicians in determining the choice of treatment as well as providing important epidemiological data for surveillance and disease control assessment. We report the development of a rapid nucleic acid lateral flow biosensor (NALFB) in a dry-reagent strip format for the sequence-specific detection of single-stranded polymerase chain reaction (PCR) amplicons at ambient temperature (22-25°C). The NALFB was developed in combination with a linear-after-the-exponential PCR assay and the applicability of this biosensor was demonstrated through detection of the cholera toxin gene from diarrheal-causing toxigenic Vibrio cholerae. Amplification using the advanced asymmetric PCR boosts the production of fluorescein-labeled single-stranded amplicons, allowing capture probes immobilized on the NALFB to hybridize specifically with complementary targets in situ on the strip. Subsequent visual formation of red lines is achieved through the binding of conjugated gold nanoparticles to the fluorescein label of the captured amplicons. The visual detection limit observed with synthetic target DNA was 0.3 ng and 1 pg with pure genomic DNA. Evaluation of the NALFB with 164 strains of V. cholerae and non-V. cholerae bacteria recorded 100% for both sensitivity and specificity. The whole procedure of the low-cost NALFB, which is performed at ambient temperature, eliminates the need for preheated buffers or additional equipment, greatly simplifying the protocol for sequence-specific PCR amplicon analysis.

MELTING, a flexible platform to predict the melting temperatures of nucleic acids

Background

Computing accurate nucleic acid melting temperatures has become a crucial step for the efficiency and the optimisation of numerous molecular biology techniques such as in situ hybridization, PCR, antigene targeting, and microarrays. MELTING is a free open source software which computes the enthalpy, entropy and melting temperature of nucleic acids. MELTING 4.2 was able to handle several types of hybridization such as DNA/DNA, RNA/RNA, DNA/RNA and provided corrections to melting temperatures due to the presence of sodium. The program can use either an approximative approach or a more accurate Nearest-Neighbor approach.

Results

Two new versions of the MELTING software have been released. MELTING 4.3 is a direct update of version 4.2, integrating newly available thermodynamic parameters for inosine, a modified adenine base with an universal base capacity, and incorporates a correction for magnesium. MELTING 5 is a complete reimplementation which allows much greater flexibility and extensibility. It incorporates all the thermodynamic parameters and corrections provided in MELTING 4.x and introduces a large set of thermodynamic formulae and parameters, to facilitate the calculation of melting temperatures for perfectly matching sequences, mismatches, bulge loops, CNG repeats, dangling ends, inosines, locked nucleic acids, 2-hydroxyadenines and azobenzenes. It also includes temperature corrections for monovalent ions (sodium, potassium, Tris), magnesium ions and commonly used denaturing agents such as formamide and DMSO.

Conclusions

MELTING is a useful and very flexible tool for predicting melting temperatures using approximative formulae or Nearest-Neighbor approaches, where one can select different sets of Nearest-Neighbor parameters, corrections and formulae. Both versions are freely available at http://sourceforge.net/projects/melting/and at http://www.ebi.ac.uk/compneur-srv/melting/under the terms of the GPL license.

Non-specific binding of Na+ and Mg2+ to RNA determined by force spectroscopy methods

RNA duplex stability depends strongly on ionic conditions, and inside cells RNAs are exposed to both monovalent and multivalent ions. Despite recent advances, we do not have general methods to quantitatively account for the effects of monovalent and multivalent ions on RNA stability, and the thermodynamic parameters for secondary structure prediction have only been derived at 1M [Na(+)]. Here, by mechanically unfolding and folding a 20 bp RNA hairpin using optical tweezers, we study the RNA thermodynamics and kinetics at different monovalent and mixed monovalent/Mg(2+) salt conditions. We measure the unfolding and folding rupture forces and apply Kramers theory to extract accurate information about the hairpin free energy landscape under tension at a wide range of ionic conditions. We obtain non-specific corrections for the free energy of formation of the RNA hairpin and measure how the distance of the transition state to the folded state changes with force and ionic strength. We experimentally validate the Tightly Bound Ion model and obtain values for the persistence length of ssRNA. Finally, we test the approximate rule by which the non-specific binding affinity of divalent cations at a given concentration is equivalent to that of monovalent cations taken at 100-fold concentration for small molecular constructs.

Metal ion-promoted conformational changes of oligonucleotides

The review will discuss the influence of metal ions on conformational changes of oligonucleotides. First, a short definition of the torsion angles is given, followed by a concise yet critical overview of the commonly applied experimental techniques. Finally, the possible role of metals upon the following conformational changes of oligonucleotides is discussed: (i) the denaturation of double-strands, (ii) the transition from B- to A-DNA, (iii) the transition from right- to left-handed DNA and RNA, (iv) the condensation, (v) and other conformational changes. We conclude with a summary and outlook.

Enantioselective and label-free detection of oligopeptide via fluorescent indicator displacement

In this work, a simple and label-free fluorescent method via fluorescent indicator displacement (FID) was proposed for enantioselectively determining d-enantiomer of arginine vasopressin (DV) using DV-specific DNA aptamer (V-apt) and one guanidiniophthalocyanine dye (Zn-DIGP). Zn-DIGP that preferentially binds to single-stranded DNA with fluorescence enhancement rather than duplexes occupies the long internal loop of V-apt and generates intensive fluorescence. Then DV is introduced into the solution containing Zn-DIGP and V-apt, and displaces the Zn-DIGP from the binding site of internal loop, leading to fluorescence decrease. But l-enantiomer cannot induce any fluorescence change due to the selectivity of V-apt. This established FID technique can detect DV with a detection limit of 100 nM and exhibits a broad linear range, and is able to discriminate enantiomers of arginine vasopressin unambiguously. Moreover chiral separation by chromatography, complicated experimental procedures and covalent modification of tags (such as organic dyes, redox-active metal complexes) are avoided in our strategy. This simple and label-free method is promising for fabricating diverse aptasensors to determine other biomolecules and drugs.

Electrostatics of nucleic acid folding under conformational constraint

RNA folding is enabled by interactions between the nucleic acid and its ion atmosphere, the mobile sheath of aqueous ions that surrounds and stabilizes it. Understanding the ion atmosphere requires the interplay of experiment and theory. However, even an apparently simple experiment to probe the ion atmosphere, measuring the dependence of DNA duplex stability upon ion concentration and identity, suffers from substantial complexity, because the unfolded ensemble contains many conformational states that are difficult to treat accurately with theory. To minimize this limitation, we measured the unfolding equilibrium of a DNA hairpin using a single-molecule optical trapping assay, in which the unfolded state is constrained to a limited set of elongated conformations. The unfolding free energy increased linearly with the logarithm of monovalent cation concentration for several cations, such that smaller cations tended to favor the folded state. Mg(2+) stabilized the hairpin much more effectively at low concentrations than did any of the monovalent cations. Poisson-Boltzmann theory captured trends in hairpin stability measured for the monovalent cation titrations with reasonable accuracy, but failed to do so for the Mg(2+) titrations. This finding is consistent with previous work, suggesting that Poisson-Boltzmann and other mean-field theories fail for higher valency cations where ion-ion correlation effects may become significant. The high-resolution data herein, because of the straightforward nature of both the folded and the unfolded states, should serve as benchmarks for the development of more accurate electrostatic theories that will be needed for a more quantitative and predictive understanding of nucleic acid folding.

A DNA-based nanomechanical device used to characterize the distortion of DNA by Apo-SoxR protein

DNA-based nanomechanical devices can be used to characterize the action of DNA-distorting proteins. Here, we have constructed a device wherein two DNA triple-crossover (TX) molecules are connected by a shaft, similar to a previous device that measured the binding free energy of integration host factor. In our case, the binding site on the shaft contains the sequence recognized by SoxR protein, the apo form of which is a transcriptional activator. Another active form is oxidized [2Fe-2S] SoxR formed during redox sensing, and previous data suggest that activated Fe-SoxR distorts its binding site by localized DNA untwisting by an amount that corresponds to ~2 bp. A pair of dyes report the fluorescence resonance energy transfer (FRET) signal between the two TX domains, reflecting changes in the shape of the device upon binding of the protein. The TX domains are used to amplify the signal expected from a relatively small distortion of the DNA binding site. From FRET analysis of apo-SoxR binding, the effect of apo-SoxR on the original TX device is similar to the effect of shortening the TX device by 2 bp. We estimate that the binding free energy of apo-SoxR on the DNA target site is 3.2-6.1 kcal/mol.

A coarse-grain three-site-per-nucleotide model for DNA with explicit ions

The "three sites per nucleotide" (3SPN) model provides a coarse-grained representation of nucleic acids for simulation of molecular processes. Previously, this model has relied on an implicit representation of the surrounding ionic environment at the level of Debye-Hückel theory. In this work, we eliminate this limitation and present an explicit representation of ions, both monovalent and divalent. The coarse-grain ion-ion and ion-phosphate potential energy functions are inferred from all-atom simulations and parameterized to reproduce key features of the local structure and organization of ions in bulk water and in the presence of DNA. The resulting model, 3SPN.1-I, is capable of reproducing the local structure observed in detailed atomistic simulations, as well as the experimental melting temperature of DNA for a range of DNA oligonucleotide lengths, CG-content, Na(+) concentration, and Mg(2+) concentration.

Mechanism of translocation of uracil-DNA glycosylase from Escherichia coli between distributed lesions

Uracil-DNA glycosylase (Ung) is a DNA repair enzyme that excises uracil bases from DNA, where they appear through deamination of cytosine or incorporation from a cellular dUTP pool. DNA repair enzymes often use one-dimensional diffusion along DNA to accelerate target search; however, this mechanism remains poorly investigated mechanistically. We used oligonucleotide substrates containing two uracil residues in defined positions to characterize one-dimensional search of DNA by Escherichia coli Ung. Mg(2+) ions suppressed the search in double-stranded DNA to a higher extent than K(+) likely due to tight binding of Mg(2+) to DNA phosphates. Ung was able to efficiently overcome short single-stranded gaps within double-stranded DNA. Varying the distance between the lesions and fitting the data to a theoretical model of DNA random walk, we estimated the characteristic one-dimensional search distance of ~100 nucleotides and translocation rate constant of ~2×10(6) s(-1).

Enhanced annealing of mismatched oligonucleotides using a novel melting curve assay allows efficient in vitro discrimination and restriction of a single nucleotide polymorphism

Background

Many SNP discrimination strategies employ natural restriction endonucleases to discriminate between allelic states. However, SNPs are often not associated with a restriction site and therefore, a number of attempts have been made to generate sequence-adaptable restriction endonucleases. In this study, a simple, sequence-adaptable SNP discrimination mechanism between a 'wild-type' and 'mutant' template is demonstrated. This model differs from other artificial restriction endonuclease models as cis- rather than trans-orientated regions of single stranded DNA were generated and cleaved, and therefore, overcomes potential issues of either inefficient or non-specific binding when only a single variant is targeted.

Results

A series of mismatch 'bubbles' that spanned 0-5-bp surrounding a point mutation was generated and analysed for sensitivity to S1 nuclease. In this model, generation of oligonucleotide-mediated ssDNA mismatch 'bubbles' in the presence of S1 nuclease resulted in the selective degradation of the mutant template while maintaining wild-type template integrity. Increasing the size of the mismatch increased the rate of mutant sequence degradation, until a threshold above which discrimination was lost and the wild-type sequence was degraded. This level of fine discrimination was possible due to the development of a novel high-resolution melting curve assay to empirically determine changes in Tm (~5.0°C per base-pair mismatch) and to optimise annealing conditions (~18.38°C below Tm) of the mismatched oligonucleotide sets.

Conclusions

The in vitro 'cleavage bubble' model presented is sequence-adaptable as determined by the binding oligonucleotide, and hence, has the potential to be tailored to discriminate between any two or more SNPs. Furthermore, the demonstrated fluorometric assay has broad application potential, offering a rapid, sensitive and high-throughput means to determine Tm and annealing rates as an alternative to conventional hybridisation detection strategies.

Divalent metal cation speciation and binding to surface-bound oligonucleotide single strands studied by second harmonic generation

The binding of Sr(II), Ca(II), Mg(II), Ba(II), Mn(II), Zn(II), and Cd(II) to silica/water interfaces functionalized with A(15)T(6) oligonucleotides was quantified at pH 7 and 10 mM NaCl using the Eisenthal χ((3)) technique. The binding free energies range from -31.1(6) kJ/mol for Ba(II) to -33.8(4) kJ/mol for Ca(II). The ion densities were found to range from 2(1) ions/strand for Zn(II) to 11(1) ions/strand for Cd(II). Additionally, we quantified Mg(II) binding in the presence of varying background electrolyte concentrations which showed that the binding free energies changed in a linear fashion from -39.3(8) to -27(1) kJ/mol over the electrolyte concentration range of 1-80 mM, respectively. An adsorption free energy versus interfacial potential analysis allowed us to elucidate the speciation of the bound Mg(II) ions and to identify three possible binding pathways. Our findings suggest that Mg(II) binds as a fully hydrated divalent cation, most likely displacing DNA-bound Na ions. These measurements will serve as a benchmark for computer simulations of divalent metal cation/DNA interactions for geochemical and biosensing applications.

1,2-propanediol-trehalose mixture as a potent quantitative real-time PCR enhancer

Background

Quantitative real-time PCR (qPCR) is becoming increasingly important for DNA genotyping and gene expression analysis. For continuous monitoring of the production of PCR amplicons DNA-intercalating dyes are widely used. Recently, we have introduced a new qPCR mix which showed improved amplification of medium-size genomic DNA fragments in the presence of DNA dye SYBR green I (SGI). In this study we tested whether the new PCR mix is also suitable for other DNA dyes used for qPCR and whether it can be applied for amplification of DNA fragments which are difficult to amplify.

Results

We found that several DNA dyes (SGI, SYTO-9, SYTO-13, SYTO-82, EvaGreen, LCGreen or ResoLight) exhibited optimum qPCR performance in buffers of different salt composition. Fidelity assays demonstrated that the observed differences were not caused by changes in Taq DNA polymerase induced mutation frequencies in PCR mixes of different salt composition or containing different DNA dyes. In search for a PCR mix compatible with all the DNA dyes, and suitable for efficient amplification of difficult-to-amplify DNA templates, such as those in whole blood, of medium size and/or GC-rich, we found excellent performance of a PCR mix supplemented with 1 M 1,2-propanediol and 0.2 M trehalose (PT enhancer). These two additives together decreased DNA melting temperature and efficiently neutralized PCR inhibitors present in blood samples. They also made possible more efficient amplification of GC-rich templates than betaine and other previously described additives. Furthermore, amplification in the presence of PT enhancer increased the robustness and performance of routinely used qPCRs with short amplicons.

Conclusions

The combined data indicate that PCR mixes supplemented with PT enhancer are suitable for DNA amplification in the presence of various DNA dyes and for a variety of templates which otherwise can be amplified with difficulty.

Characterization of the adsorption of oligonucleotides on mercaptopropionic acid-coated CdSe/ZnS quantum dots using fluorescence resonance energy transfer

Semiconductor quantum dots (QDs) coated with thioalkyl acid ligands are often used as probes and reporters for nucleic acid sensing, or protein sensing using aptamers, and are also potential vectors for gene delivery. In such applications, the interactions that potentially lead to the adsorption of oligonucleotides onto the surface of colloidal QDs are an important consideration. To explore such interactions, fluorescence resonance energy transfer (FRET) between QDs and oligonucleotides labeled with a fluorescent dye was used to identify and characterize a set of conditions that favor strong adsorption on 3-mercaptopropionic acid (MPA)-coated CdSe/ZnS QDs. Adsorption curves and competitive binding experiments were used to determine that the order of affinity for nucleobase adsorption was dC>dA≥dG≫dT. The length of the oligonucleotide sequence was also important, with an 80-mer sequence adsorbing more strongly than its 20-mer analog. Adsorption decreased with increasing pH and corresponded to the ionization of the carboxylic acid groups of the MPA ligands. Increased ionic strength partially offsets ligand ionization and increased the extent of adsorption. The interaction between QDs and oligonucleotides was labile, with increases in adsorption at lower concentrations of oligonucleotide and with an increasing number of oligonucleotides per QD. The results were consistent with a hydrogen-bonding model for adsorption, where neutral thioalkyl acid ligands interact favorably with nucleobases and ionized ligands resist adsorption.

Controlling forces and pathways in self-assembly using viruses and DNA

The ability of both viruses and DNA to self-assemble in solution has continues to enable numerous applications at the nanoscale. Here we review the relevant interactions dictating the assembly of these structures, as well as discussing how they can be exploited experimentally. Because self-assembly is a process, we discuss various strategies for achieving spatial and temporal control. Finally, we highlight a few examples of recent advances that exploit the features of these nanostructures.

Specific and nonspecific metal ion-nucleotide interactions at aqueous/solid interfaces functionalized with adenine, thymine, guanine, and cytosine oligomers

This article reports nonlinear optical measurements that quantify, for the first time directly and without labels, how many Mg(2+) cations are bound to DNA 21-mers covalently linked to fused silica/water interfaces maintained at pH 7 and 10 mM NaCl, and what the thermodynamics are of these interactions. The overall interaction of Mg(2+) with adenine, thymine, guanine, and cytosine is found to involve -10.0 ± 0.3, -11.2 ± 0.3, -14.0 ± 0.4, and -14.9 ± 0.4 kJ/mol, and nonspecific interactions with the phosphate and sugar backbone are found to contribute -21.0 ± 0.6 kJ/mol for each Mg(2+) ion bound. The specific and nonspecific contributions to the interaction energy of Mg(2+) with oligonucleotide single strands is found to be additive, which suggests that within the uncertainty of these surface-specific experiments, the Mg(2+) ions are evenly distributed over the oligomers and not isolated to the most strongly binding nucleobase. The nucleobases adenine and thymine are found to bind only three Mg(2+) ions per 21-mer oligonucleotide, while the bases cytosine and guanine are found to bind eleven Mg(2+) ions per 21-mer oligonucleotide.

Monovalent and divalent salt correction algorithms for Tm prediction--recommendations for Primer3 usage

Primer3 is a widely used program for selection of oligonucleotide primers for PCR. The websites used for implementation of Primer3 have recently been updated. PCR requires Mg(2+)(,) which has a significant dsDNA stabilizing effect that must be taken into account when designing PCR primers. The data sets and formulas used to correct for salt concentrations have been updated in Primer3 to give better prediction of melting temperature (T(m)). The liberal combination of different formulas for monovalent and divalent salt correction can lead to different results, depending on the formula chosen by the user. Using published T(m) for 475 different oligonucleotides, it is shown that the combination of the implemented conversion of divalent to monovalent cation concentration works well with one salt correction formula but not with an alternative one. Use of a more recently described alternative formula would lead to equally good T(m) predictions if divalent cations are present. The proper selection of compatible primer pairs depends on the choice of a good combination of salt correction formulas. Currently the SantaLucia salt correction formula should be used if Mg(2+) is present. The alternative formula should be updated to its recent form for future releases.

Cooperative hybridization of oligonucleotides

Nucleic acids have been demonstrated to be versatile nanoscale engineering materials with the construction of dynamic DNA structures, motors, and circuits. These constructions generally rely on the clever use and integration of relatively few reaction mechanisms and design primitives. Here, cooperative hybridization is introduced as a mechanism in which two oligonucleotides of independent sequence can stoichiometrically, simultaneously, and cooperatively hybridize to a DNA complex. Cooperative hybridization is rigorously characterized and modeled and is shown to implement digital concentration comparison with amplification, as well as digital Boolean logic. These designs, based on cooperative hybridization, excel in being robust to impurities and not requiring oligonucleotide purification.

Comparing the Predictions of the Nonlinear Poisson-Boltzmann Equation and the Ion Size-Modified Poisson-Boltzmann Equation for a Low-Dielectric Charged Spherical Cavity in an Aqueous Salt Solution

The ion size-modified Poisson Boltzmann equation (SMPBE) is applied to the simple model problem of a low-dielectric spherical cavity containing a central charge, in an aqueous salt solution to investigate the finite ion size effect upon the electrostatic free energy and its sensitivity to changes in salt concentration. The SMPBE is shown to predict a very different electrostatic free energy than the nonlinear Poisson-Boltzmann equation (NLPBE) due to the additional entropic cost of placing ions in solution. Although the energy predictions of the SMPBE can be reproduced by fitting an appropriatelysized Stern layer, or ion-exclusion layer to the NLPBE calculations, the size of the Stern layer is difficult to estimate a priori. The SMPBE also produces a saturation layer when the central charge becomes sufficiently large. Ion-competition effects on various integrated quantities such the total number of ions predicted by the SMPBE are qualitatively similar to those given by the NLPBE and those found in available experimental results.

Direct detection and quantification of methylation in nucleic acid sequences using high-resolution melting analysis

High-resolution melting (HRM) analysis exploits the reduced thermal stability of DNA fragments that contain base mismatches to detect single nucleotide polymorphisms (SNPs). However, the capacity of HRM to reveal other features of DNA chemistry remains unexplored. DNA methylation plays a key role in regulating gene expression and is essential for normal development in many higher organisms. The presence of methylated bases perturbs the double-stranded DNA structure, although its effect on DNA thermal stability is largely unknown. Here, we reveal that methylated DNA has enhanced thermal stability and is sufficiently divergent from nonmethylated DNA to allow detection and quantification by HRM analysis. This approach reliably distinguishes between sequence-identical DNA differing only in the methylation of one base. The method also provides accurate discrimination between mixes of methylated and nonmethylated DNAs, allowing discrimination between DNA that is 1% and 0% methylated and also between 97.5% and 100% methylated. Thus, the method provides a new means of adjusting thermal optima for DNA hybridization and PCR-based techniques and to empirically measure the impact of DNA methylation marks on the thermostability of regulatory regions. In the longer term, it could enable the development of new techniques to quantify methylated DNA.

Human DHX9 helicase unwinds triple-helical DNA structures

Naturally occurring poly(purine.pyrimidine) rich regions in the human genome are prone to adopting non-canonical DNA structures such as intramolecular triplexes (i.e., H-DNA). Such structure-forming sequences are abundant and can regulate the expression of several disease-linked genes. In addition, the use of triplex-forming oligonucleotides (TFOs) to modulate gene structure and function has potential as an approach to targeted gene therapy. Previously, we found that endogenous H-DNA structures can induce DNA double-strand breaks and promote genomic rearrangements. Herein, we find that the DHX9 helicase co-immunoprecipitates with triplex DNA structures in mammalian cells, suggesting a role in the maintenance of genome stability. We tested this postulate by assessing the helicase activity of purified human DHX9 on various duplex and triplex DNA substrates in vitro. DHX9 displaced the third strand from a specific triplex DNA structure and catalyzed the unwinding with a 3' --> 5' polarity with respect to the displaced third strand. Helicase activity required a 3'-single-stranded overhang on the third strand and was dependent on ATP hydrolysis. The reaction kinetics consisted of a pre-steady-state burst phase followed by a linear, steady-state pseudo-zero-order reaction. In contrast, very little if any helicase activity was detected on blunt triplexes, triplexes with 5'-overhangs, blunt duplexes, duplexes with overhangs, or forked duplex substrates. Thus, triplex structures containing a 3'-overhang represent preferred substrates for DHX9, where it removes the strand with Hoogsteen hydrogen-bonded bases. Our results suggest the involvement of DHX9 in maintaining genome integrity by unwinding mutagenic triplex DNA structures.

Chemical primer extension at submillimolar concentration of deoxynucleotides

Template-directed primer extension usually requires a polymerase, nucleoside triphosphates, and magnesium ions as cofactors. Enzyme-free, chemical primer extensions are known for preactivated nucleotides at millimolar concentrations. Based on a screen of carbodiimides, heterocyclic catalysts, and reactions conditions, we now show that near-quantitative primer conversion can be achieved at submillimolar concentration of any of the four deoxynucleotides (dAMP, dCMP, dGMP and dTMP). The new protocol relies on in situ activation with EDC and 1-methylimidazole and a magnesium-free buffer that was tested successfully for different sequence motifs. The method greatly simplifies chemical primer extension assays, further reduces the cost of such assays, and demonstrates the potential of the in situ activation approach.

Shielding effect of monovalent and divalent cations on solid-phase DNA hybridization: surface plasmon resonance biosensor study

Solid-phase hybridization, i.e. the process of recognition between DNA probes immobilized on a solid surface and complementary targets in a solution is a central process in DNA microarray and biosensor technologies. In this work, we investigate the simultaneous effect of monovalent and divalent cations on the hybridization of fully complementary or partly mismatched DNA targets to DNA probes immobilized on the surface of a surface plasmon resonance sensor. Our results demonstrate that the hybridization process is substantially influenced by the cation shielding effect and that this effect differs substantially for solid-phase hybridization, due to the high surface density of negatively charged probes, and hybridization in a solution. In our study divalent magnesium is found to be much more efficient in duplex stabilization than monovalent sodium (15 mM Mg²⁺ in buffer led to significantly higher hybridization than even 1 M Na⁺). This trend is opposite to that established for oligonucleotides in a solution. It is also shown that solid-phase duplex destabilization substantially increases with the length of the involved oligonucleotides. Moreover, it is demonstrated that the use of a buffer with the appropriate cation composition can improve the discrimination of complementary and point mismatched DNA targets.

Rationally designed aptamer-based fluorescence polarization sensor dedicated to the small target analysis

A direct fluorescence polarization (FP) assay strategy, dedicated to the small molecule sensing and based on the unique induced-fit binding mechanism of end-labelled nucleic acid aptamers, has been recently developed by our group. Small target binding has been successfully converted into a significant increase of the fluorescence anisotropy signal presumably produced by the reduction of the local motional freedom of the dye. In order to generalize the approach, a rational FP sensor methodology was established herein, by engineering instability in the secondary structure of an aptameric recognition element. The anti-adenosine DNA aptamer, labelled by a single fluorescein dye at its 3' extremity, was employed as a model functional nucleic acid probe. The terminal stem of the stem-loop structure was shortened to induce a destabilized/denatured conformation which promoted the local segmental mobility of the dye and then a significant depolarization process. Upon target binding, the structural change of the aptamer induced the formation of a stable stem-loop structure, leading to the reduction of the dye mobility and the increase in the fluorescence anisotropy signal. This reasoned approach was applied to the sensing of adenosine and adenosine monophosphate and their chiral analysis.

Noncompetitive fluorescence polarization aptamer-based assay for small molecule detection

In this paper, a new fluorescence polarization (FP) assay strategy is described reporting the first demonstration of a noncompetitive FP technique dedicated to the small molecule sensing. This approach was based on the unique induced-fit binding mechanism of nucleic acid aptamers which was exploited to convert the small target binding event into a detectable fluorescence anisotropy signal. An anti-L-tyrosinamide DNA aptamer, labeled by a single fluorescent dye at its extremity, was employed as a model functional nucleic acid probe. The DNA conformational change generated by the L-tyrosinamide binding was able to induce a significant increase in the fluorescence anisotropy signal. The method allowed enantioselective sensing of tyrosinamide and analysis in practical samples. The methodology was also applied to the L-argininamide detection, suggesting the potential generalizability of the direct FP-based strategy. Such aptamer-based assay appeared to be a sensitive analytical system of remarkable simplicity and ease of use.

Zip Nucleic Acids: new high affinity oligonucleotides as potent primers for PCR and reverse transcription

Most nucleic acid-based technologies rely upon sequence recognition between an oligonucleotide and its nucleic acid target. With the aim of improving hybridization by decreasing electrostatic repulsions between the negatively charged strands, novel modified oligonucleotides named Zip nucleic acids (ZNAs) were recently developed. ZNAs are oligonucleotide–oligocation conjugates whose global charge is modulated by the number of cationic spermine moieties grafted on the oligonucleotide. It was demonstrated that the melting temperature of a hybridized ZNA is easily predictable and increases linearly with the length of the oligocation. Furthermore, ZNAs retain the ability to discriminate between a perfect match and a single base-pair-mismatched complementary sequence. Using quantitative PCR, we show here that ZNAs are specific and efficient primers displaying an outstanding affinity toward their genomic target. ZNAs are particularly efficient at low magnesium concentration, low primer concentrations and high annealing temperatures, allowing to improve the amplification in AT-rich sequences and potentially multiplex PCR applications. In reverse transcription experiments, ZNA gene-specific primers improve the yield of cDNA synthesis, thus increasing the accuracy of detection, especially for genes expressed at low levels. Our data suggest that ZNAs exhibit faster binding kinetics than standard and locked nucleic acid-containing primers, which could explain why their target recognition is better for rare targets.

Probing the limits of liquid droplet laser desorption mass spectrometry in the analysis of oligonucleotides and nucleic acids

In the present work we demonstrate the advantages of LILBID mass spectrometry (laser-induced liquid bead ion desorption) in the analysis of nucleic acids and large oligonucleotides. For established methods like matrix-assisted laser desorption/ionization (MALDI) and electrospray ionization (ESI), the mass analysis of oligonucleotides or of noncovalent oligonucleotide-protein complexes, in particular of very large ones, still represents a considerable challenge either due to the lack of native solutions or nonspecific adduct formation or due to a reduced salt tolerance or a high charge state of the ions. With LILBID, oligonucleotides, solvated in micro-droplets of aqueous buffer at certain pH and ion strength, are brought into the gas phase by laser ablation. We show that our method is able to detect single- and double-stranded oligonucleotides with high softness, demonstrated by the buffer dependence of the melting of a duplex. The absolute sensitivity is in the attomole range concomitant with a total analyte consumption in the femtomole region. The upper mass limit of oligonucleotides still detected with good signal-to-noise ratio with LILBID is the 1.66 MDa plasmid pUC19. With DNA ladders from short duplexes with sticky ends, we show that LILBID correctly reflects the relative thermodynamic stabilities of the ladders. Moreover, as an example for a specific DNA-protein complex we show that a NF-kappaB p50 homodimer binds sequence specifically to its match DNA. In summary we demonstrate that LILBID, although presently performed only with low mass resolution, due to these advantages, is an alternative mass spectrometric method for the analysis of oligonucleotides in general and of specific noncovalent nucleic acid-protein complexes in particular.

TmPrime: fast, flexible oligonucleotide design software for gene synthesis

Herein we present TmPrime, a computer program to design oligonucleotide sets for gene assembly by both ligase chain reaction (LCR) and polymerase chain reaction (PCR). TmPrime offers much flexibility with no constraints on the gene and oligonucleotide lengths. The program divides the long input DNA sequence based on the input desired melting temperature, and dynamically optimizes the length of oligonucleotides to achieve homologous melting temperatures. The output reports the melting temperatures, oligonucleotide sequences and potential formation of secondary structures. Our program also provides functions on sequence pooling to separate long genes into smaller pieces for multi-pool assembly and codon optimization for expression. The software has been successfully used in the design and synthesis of green fluorescent protein fragment (GFPuv) (760 bp), human protein kinase B-2 (PKB2) (1446 bp) and the promoter of human calcium-binding protein A4 (S100A4) (752 bp) using real-time PCR assembly with LCGreen I, which offers a novel approach to compare the efficiency of gene synthesis. The purity of assembled products is successfully estimated with the use of melting curve analysis, which would potentially eliminate the necessity for agarose gel electrophoresis. This program is freely available at http://prime.ibn.a-star.edu.sg.

Molecular simulation studies of monovalent counterion-mediated interactions in a model RNA kissing loop

A kissing loop is a highly stable complex formed by loop-loop base-pairing between two RNA hairpins. This common structural motif is utilized in a wide variety of RNA-mediated processes, including antisense recognition, substrate recognition in riboswitches, and viral replication. Recent work has shown that the Tar-Tar(*) complex, an archetypal kissing loop, can form without Mg(2+), so long as high concentrations of alkali chloride salts are present. Interestingly, the stability of the complex is found to decrease with increasing cation size. In this work, we used molecular simulations to develop a molecular-level understanding of the origins of the observed counterion specificity. The ionic atmosphere of the Tar-Tar(*) complex was examined in the presence of 800 mm (where m denotes molality) NaCl, KCl, or CsCl. We used spatial free-energy density profiles to analyze differences in counterion accumulation at different spatial extents from the RNA molecule. We found that the lowest free-energy levels, which are situated in the vicinity of the loop-loop interface, can accommodate roughly two counterions, irrespective of counterion type. However, as we moved into higher free-energy levels, which are farther from the loop-loop interface, we observed increased differences in the numbers of accumulated counterions, with this number being largest for Na(+) and smallest for Cs(+). We analyzed the source of these differences and were able to attribute these to two distinct features: The extent of partial dehydration varies based on cation type; the smaller the cation, the greater the degree of dehydration. While smaller ions bind their first-hydration-shell water molecules more tightly than larger ions, they are also able to shed these water molecules for stronger electrostatic interactions with the RNA molecule. Secondly, we observed a distinct asymmetry in the numbers of accumulated cations around each hairpin in the Tar-Tar(*) complex. We were able to ascribe this asymmetry to the presence of a guanine tract in the Tar hairpin, which facilitates partial dehydration of the counterions. However, the smaller ions compensate for this asymmetry by forming a belt around the loop-loop interface in intermediate free-energy levels. As a result, the degree of asymmetry in counterion accumulation around individual hairpins shows an inverse correlation with the experimentally observed cation specificity for the stability of Tar-Tar(*) (i.e., the smaller the asymmetry, the greater the experimentally observed stability). This in turn provides a plausible explanation for why the smaller cations help stabilize the Tar-Tar(*) complex better than the larger cations. These findings suggest that the specific sequence and structural features of the Tar-Tar(*) complex may be the source of the experimentally observed cation specificity in Tar-Tar(*) stability. Our results lead to testable predictions for how changes in sequence might alter the observed counterion specificity in kissing loop stability.

IDT SciTools: a suite for analysis and design of nucleic acid oligomers

DNA and RNA oligomers are used in a myriad of diverse biological and biochemical experiments. These oligonucleotides are designed to have unique biophysical, chemical and hybridization properties. We have created an integrated set of bioinformatics tools that predict the properties of native and chemically modified nucleic acids and assist in their design. Researchers can select PCR primers, probes and antisense oligonucleotides, find the most suitable sequences for RNA interference, calculate stable secondary structures, and evaluate the potential for two sequences to interact. The latest, most accurate thermodynamic algorithms and models are implemented. This free software is available at http://www.idtdna.com/SciTools/SciTools.aspx.