Evaluating single-base-pair discriminating capability of planar oligonucleotide microchips using a non-equilibrium dissociation approach

The capability of planar rRNA-based oligonucleotide microarrays for single-base-pair discrimination was evaluated using an approach that compares the non-equilibrium dissociation profiles and dissociation temperatures (Tds) of all probe-target duplexes simultaneously. Three sets of 16S rRNA gene specific probes at different levels of specificity were used along with their counter probes for individual sets having either one or two mismatches (MM) to their targets at specific external (next to terminus) and various internal positions. Criteria based on the Td approach and a discrimination index (DI) were proven to be competent in discriminating PM from internal MM duplexes, but not always for external MM duplexes. Maximal DI for separating PM duplexes from ones with two and one internal MM usually occurred at temperatures approximately 5-10 degrees C and 10-15 degrees C, respectively, higher than the Tds of the PM duplexes. Washing buffer type and salt concentration, and MM number and position were shown statistically to affect dissociation profiles, Td, and single-base-pair discriminating capability. The reusability potential of the planar microchip was further demonstrated.     

Concentration dependency of nonequilibrium thermal dissociation curves in complex target samples

The nonequilibrium thermal dissociation (NTD) methodology has been proposed to provide a superior discrimination between specific and nonspecific hybridizations than the commonly used array techniques involving hybridization followed by a single stringent wash. Multiple studies have used this method on gel-pad, planar, and nylon membrane arrays to identify specific microbial targets in complex target mixtures. A recent physicochemical study revealed several problems, particularly when the method was used to examine complex target samples. In the present study, we investigated the effect of target concentration on NTD of complex target samples obtained from an anaerobic bioreactor. Our purpose was to experimentally demonstrate that variation in the concentrations of both specific and nonspecific targets determines the course of dissociation, which was not evaluated in initial microbiological studies. We also present an approach for analyzing the dissociation curves that is less error prone compared to those used in the previous studies. Our results show that: (i) a specific target in a mixture, at a certain concentration, may have a higher dissociation temperature/time than that of the same pure target, and (ii) the concentration dependence of the dissociation precludes usage of reference curves for identifying a target. Contrary to the previous studies, an explicit calibration is required, which makes the NTD approach impractical for high throughput analysis.     

Specific and nonspecific hybridization of oligonucleotide probes on microarrays

Gene expression analysis by means of microarrays is based on the sequence-specific binding of RNA to DNA oligonucleotide probes and its measurement using fluorescent labels. The binding of RNA fragments involving sequences other than the intended target is problematic because it adds a chemical background to the signal, which is not related to the expression degree of the target gene. The article presents a molecular signature of specific and nonspecific hybridization with potential consequences for gene expression analysis. We analyzed the signal intensities of perfect match (PM) and mismatch (MM) probes of GeneChip microarrays to specify the effect of specific and nonspecific hybridization. We found that these events give rise to different relations between the PM and MM intensities as function of the middle base of the PM, namely a triplet-like (C > G approximately T > A > 0) and a duplet-like (C approximately T > 0 > G approximately A) pattern of the PM-MM log-intensity difference upon binding of specific and nonspecific RNA fragments, respectively. The systematic behavior of the intensity difference can be rationalized on the level of basepairings of DNA/RNA oligonucleotide duplexes in the middle of the probe sequence. Nonspecific binding is characterized by the reversal of the central Watson-Crick (WC) pairing for each PM/MM probe pair, whereas specific binding refers to the combination of a WC and a self-complementary (SC) pairing in PM and MM probes, respectively. The Gibbs free energy contribution of WC pairs to duplex stability is asymmetric for purines and pyrimidines of the PM and decreases according to C > G approximately T > A. SC pairings on the average only weakly contribute to duplex stability. The intensity of complementary MM introduces a systematic source of variation which decreases the precision of expression measures based on the MM intensities.     

Optimization of an oligonucleotide microchip for microbial identification studies: a non-equilibrium dissociation approach

The utility of a high-density oligonucleotide microarray (microchip) for identifying strains of five closely related bacilli ( Bacillus anthracis , Bacillus cereus , Bacillus mycoides , Bacillus medusa and Bacillus subtilis ) was demonstrated using an approach that compares the non-equilibrium dissociation rates ('melting curves') of all probe–target duplexes simultaneously. For this study, a hierarchical set of 30 oligonucleotide probes targeting the 16S ribosomal RNA of these bacilli at multiple levels of specificity (approximate taxonomic ranks of domain, kingdom, order, genus and species) was designed and immobilized in a high-density matrix of gel pads on a glass slide. Reproducible melting curves for probes with different levels of specificity were obtained using an optimized salt concentration. Clear discrimination between perfect match (PM) and mismatch (MM) duplexes was achieved. By normalizing the signals to an internal standard (a universal probe), a more than twofold discrimination (> 2.4×) was achieved between PM and 1-MM duplexes at the dissociation temperature at which 50% of the probe–target duplexes remained intact. This provided excellent differentiation among representatives of different Bacillus species, both individually and in mixtures of two or three. The overall pattern of hybridization derived from this hierarchical probe set also provided a clear 'chip fingerprint' for each of these closely related Bacillus species.     

Bleomycin fragmentation of duplex DNA occurs as staggered single-strand scissions.

Electron microscopy of purified full-length linear duplex molecules produced by bleomycin reaction with PM2 DNA revealed low frequencies of closed circular duplex molecules as well as linear duplex molecules with opposed ends (cyclized molecules which have dissociated to yield a gap between the termini). The occurrence of these latter forms indicates that double-strand scissions produced by bleomycin reaction consist of two single-strand scissions which are physically staggered on the complementary strands. Analysis of the temperature dependence for cyclization led to the estimate that an average of 1.7 +/- 0.44 base-pairs (2.6 +/- 0.5 base pairs without base-stacking energies) occur between the staggered breaks. The reassociated termini cannot be ligated with T4 ligase. When PM2 DNA was fragmented at several sites within each molecule, circular duplexes and linear duplexes with opposed ends with a range of sizes from 350 base pairs up to full-length PM2 DNA were observed. Analysis of the frequency distribution of lengths of these fragments indicates that most, if not all, of the specific sites for bleomycin-directed double-strand scissions in PM2 DNA contain representatives of the same two base single-stranded termini.     

On-chip non-equilibrium dissociation curves and dissociation rate constants as methods to assess specificity of oligonucleotide probes

Nucleic acid hybridization serves as backbone for many high-throughput systems for detection, expression analysis, comparative genomics and re-sequencing. Specificity of hybridization between probes and intended targets is always critical. Approaches to ensure and evaluate specificity include use of mismatch probes, obtaining dissociation curves rather than single temperature hybridizations, and comparative hybridizations. In this study, we quantify effects of mismatch type and position on intensity of hybridization signals and provide a new approach based on dissociation rate constants to evaluate specificity of hybridized signals in complex target mixtures. Using an extensive set of 18mer oligonucleotide probes on an in situ synthesized biochip platform, we demonstrate that mismatches in the center of the probe are more discriminating than mismatches toward the extremities of the probe and mismatches toward the attached end are less discriminating than those toward the loose end. The observed destabilizing effect of a mismatch type agreed in general with predictions using the nearest neighbor model. Use of a new parameter, specific dissociation temperature (T_d-w, temperature of maximum specific dissociation rate constant), obtained from probe–target duplex dissociation profiles considerably improved the evaluation of specificity. These results have broad implications for hybridization data obtained from complex mixtures of nucleic acids.     

Simultaneous quantification of multiple nucleic acid targets in complex rRNA mixtures using high density microarrays and nonspecific hybridization as a source of information

To date, it has been problematic to accurately quantify multiple nucleic acid sequences, representing microbial targets, in multi-target mixtures using oligonucleotide microarrays, primarily due to nonspecific target binding (i.e., cross-hybridization). While some studies ignore the effects of nonspecific binding, other studies have developed approaches to minimize nonspecific binding, such as physical modeling to design highly specific probes, subtracting nonspecific signal using mismatch probes, and/or removing nonspecific duplexes by scanning through a range of wash stringencies. We have developed an alternative approach that, in contrast to previous approaches, uses nonspecific target binding as a source of information. Specifically, the new approach uses hybridization patterns (fingerprints) to quantify specific nucleic acid targets in complex target mixtures. We evaluated the approach by mixing together in vitro transcribed 28S rRNA targets at varying concentrations (up to 1.0 nM), and hybridizing the 24 mixtures to microarrays (n=3160 probes, in duplicate). Three independent Latin-square-designed experiments revealed accurate quantification of the targets. The regression between actual concentration of targets and those determined by the approach were highly positively correlated with high R(2) values (e.g., R(2)=0.90, n=6 targets; R(2)=0.84, n=8 targets; R(2)=0.82, n=10 targets).     

HU binding to DNA: evidence for multiple complex formation and DNA bending

HU, a nonspecific histone-like DNA binding protein, participates in a number of genomic events as an accessory protein and forms multiple complexes with DNA. The HU-DNA binding interaction was characterized by fluorescence, generated with the guanosine analogue 3-methyl-8-(2-deoxy-beta-D-ribofuranosyl)isoxanthopterin (3-MI) directly incorporated into DNA duplexes. The stoichiometry and equilibrium binding constants of complexes formed between HU and 13 and 34 bp DNA duplexes were determined using fluorescence anisotropy and analytical ultracentrifugation. These measurements reveal that three HU molecules bind to the 34 bp duplexes, while two HU molecules bind to the 13 bp duplex. The data are well described by an independent binding site model, and the association constants for the first binding event for both duplexes are similar (approximately 1 x 10(6) M(-1)), indicating that HU binding affinity is independent of duplex length. Further analysis of the binding curves in terms of a nonspecific binding model is indicative that HU binding to DNA exhibits little to no cooperativity. The fluorescence intensity also increases upon HU binding, consistent with decreased base stacking and increased solvent exposure of the 3-MI fluorescence probe. These results are suggestive of a local bending or unwinding of the DNA. On the basis of these results we propose a model in which bending of DNA accompanies HU binding. Up to five complex bands are observed in gel mobility shift assays of HU binding to the 34 bp duplexes. We suggest that protein-induced bending of the DNA leads to the observation of complexes in the gel, which have the same molecular weight but different relative mobilities.     

Development of a statistically robust quantification method for microorganisms in mixtures using oligonucleotide microarrays

High-density oligonucleotide arrays can be extremely useful for identifying and quantifying specific targets (i.e., ribosomal RNA of microorganisms) in mixtures. However, current array identification schemes are severely compromised by nonspecific hybridization, resulting in numerous false-positive and false-negative calls, they lack an adequate internal control for assessing the quality of identification, and are dependent on amplification of specific target sequences which introduce biases. We have developed a novel approach for the routine quantification and identification of metabolically active microorganisms in mixed samples. The advantage of our approach over conventional ones is that it avoids designing, optimizing, validating, and selecting oligonucleotide probes for arrays; also, nonspecific hybridization is no longer a problem. The basic principle of the approach is that a fluorescence pattern of a mixed sample is a superposition of the fluorescent patterns for each target. The superposition can be quantitatively deconvoluted in terms of concentrations of each microbe. We demonstrated the utility of our approach by extracting rRNA from three microorganisms, making test mixtures, labeling the rRNA, and hybridizing each test mixture to DNA oligonucleotide (20-mers, n=346,608) arrays. Comparison of known concentrations of individual targets in mixtures to those estimated by the solution revealed highly consistent results. The goodness-of-fit of the solution revealed that about 90% of the variability in the data could be explained. A new analytical approach for microbial identification and quantification has been presented in this report. Our findings demonstrate that including signal intensity values from all duplexes on the array, which are essentially nonspecific to the target organisms, significantly improved predictions of known microbial targets. To our knowledge, this is the first study to report this phenomenon. In addition, we demonstrate that the method is a self-sufficient analytical procedure since it provides statistical confidence of the quantification.     

B-S transition in short oligonucleotides

Stretching experiments with long double-stranded DNA molecules in physiological ambient revealed a force-induced transition at a force of 65 pN. During this transition between B-DNA and highly overstretched S-DNA the DNA lengthens by a factor of 1.7 of its B-form contour length. Here, we report the occurrence of this so-called B-S transition in short duplexes consisting of 30 basepairs. We employed atomic-force-microscope-based single molecule force spectroscopy to explore the unbinding mechanism of two short duplexes containing 30 or 20 basepairs by pulling at the opposite 5' termini. For a 30-basepair-long DNA duplex the B-S transition is expected to cause a length increase of 6.3 nm and should therefore be detectable. Indeed 30% of the measured force-extension curves exhibit a region of constant force (plateau) at 65 pN, which corresponds to the B-S transition. The observed plateaus show a length between 3 and 7 nm. This plateau length distribution indicates that the dissociation of a 30-basepair duplex mainly occurs during the B-S transition. In contrast, the measured force-extension curves for a 20-basepair DNA duplex exhibited rupture forces below 65 pN and did not show any evidence of a B-S transition.     

Effects of site-specific substitution of 5-fluorouridine on the stabilities of duplex DNA and RNA.

The effects of 5-fluorouridine (FUrd) and 5-fluorodeoxyuridine (FdUrd) substitution on the stabilities of duplex RNA and DNA have been studied to determine how FUrd substitution in nucleic acids may alter the efficiency of biochemical processes that require complementary base pairing for molecular recognition. The parent sequence, 5'-GCGAAUUCGC, contains two non-equivalent uridines. Eight oligonucleotides (four RNA and four DNA) were prepared with either zero, one or two Urd substituted by FUrd. The stability of each self-complementary duplex was determined by measuring the absorbance at 260 nm as a function of temperature. Tm values were calculated from the first derivative of the absorbance versus temperature profiles and values for delta H0 and delta S0 were calculated from the concentration dependence of the Tm. Individual absorbance versus temperature curves were also analyzed by a parametric approach to calculate thermodynamic parameters for the duplex to single-stranded transition. Analysis of the thermodynamic parameters for each oligonucleotide revealed that FUrd substitution had sequence-dependent effects in both A-form RNA and B-form DNA duplexes. Conservation of helix geometry in FUrd-substituted duplexes was determined by CD spectroscopy. FUrd substitution at a single site in RNA stabilized the duplex (delta delta G37 = 0.8 kcal/mol), largely due to more favorable stacking interactions. FdUrd substitution at a single site in DNA destabilized the duplex (delta delta G37 = 0.3 kcal/mol) as a consequence of less favorable stacking interactions. All duplexes melt via single cooperative transitions.     

Biaryl modification of the 5'-terminus of one strand of a microRNA duplex induces strand specificity

MicroRNAs (miRNAs) are single-stranded non-coding RNAs composed of 20-23 nucleotides. They are initially transcribed in the nucleus as pri-miRNAs. After processing, one strand from the miRNA duplex (miR-5p/miR-3p duplex) is loaded onto the RNA-induced silencing complex (RISC) to produce a functional, mature miRNA that inhibits the expression of multiple target genes. In the case of some miRNAs, both strands can be equally incorporated into the RISC as single strands, and both strands can function as mature miRNAs. Thus, a technique for selective expression of miR-5p and miR-3p strands is required to identify distinct targets of miRNAs. In this Letter, we report the synthesis and properties of miRNA duplexes carrying biaryl units at the 5'-terminus of one strand. We found that incorporation of biaryl units at the 5'-terminus of one strand of miRNA duplexes induced strand specificity in these duplexes. Further, we succeeded in identifying endogenous mRNA targets for each strand of the duplex by using the biaryl-modified miRNA duplexes.     

Effect of locked nucleic acid modifications on the thermal stability of noncanonical DNA structure

We studied the kinetic and thermodynamic effects of locked nucleic acid (LNA) modifications on parallel and antiparallel DNA duplexes. The LNA modifications were introduced at cytosine bases of the pyrimidine strand. Kinetic parameters evaluated from melting and annealing curves showed that the association and dissociation rate constants for the formation of the LNA-modified parallel duplex at 25.0 °C were 3 orders of magnitude larger and 6 orders of magnitude smaller, respectively, than that of the unmodified parallel duplex. The activation energy evaluated from the temperature-dependent rate constants was largely altered by the LNA modifications, suggesting that the LNA modifications affected a prenucleation event in the folding process. Moreover, thermodynamic parameters showed that the extent of stabilization by the LNA modification for parallel duplexes (3.6 kcal mol(-1) per one modification) was much more significant than that of antiparallel duplexes (1.6 kcal mol(-1)). This large stabilization was due to the decrease in ΔH° that was more favorable than the decrease in TΔS°. These quantitative parameters demonstrated that LNA modification specifically stabilized the noncanonical parallel duplex. On the basis of these observations, we succeeded to stabilize the parallel duplex by LNA modification at the physiological pH. These results can be useful in the rational design of functional molecules such as more effective antisense and antigene strands, more sensitive strands for detection of target DNA and RNA strands, and molecular switches responding to solution pH.     

Analysis of perfect and mismatched DNA duplexes by a generic hexanucleotide microchip

The thermodynamic analysis was done for the duplexes formed by fluorescently labeled oligonucleotide targets on a genetic hexanucleotide microchip. All 4096 different hexanucleotide chains were immobilized as probes in individual gel pads of the microchip. To strengthen the hybridization, each hexamer was extended at both ends by one nucleotide from the equimolar mixture of all four nucleotides to serve as nonselective linkers. It has been shown that the melting curves for oligonucleotide duplexes formed on the microchip and in a solution are quite similar. The influence of ionic surrounding has been studied in terms of the hybridization efficiency and discrimination between the mismatched and perfect duplexes. Different approaches have been tested to compensate the dependence of duplex stability on the GC content. It has been demonstrated that the use of chaotropic agents, addition of nonlabeled GC-rich competitor oligonucleotides, as well as creation of a temperature gradient along the microchip reproducing the distribution of melting temperatures, efficiently level out the AT/GC differences. The use of tetramethylammonium chloride for the same purpose was accompanied by weakening to some extent the discrimination between the mismatched duplexes and the perfect ones.     

Analysis of Perfect and Mismatched DNA Duplexes with a Generic Hexanucleotide Microchip

Thermodynamic analysis was performed for the duplexes formed by fluorescently labeled oligonucleotide targets on a generic hexanucleotide microchip. All 4096 different hexanucleotide chains were immobilized as probes in individual gel pads of the microchip. To strengthen the hybridization, each hexamer was extended at both ends by one nucleotide from the equimolar mixture of all four nucleotides to serve as nonselective linkers. It has been shown that the melting curves for oligonucleotide duplexes formed on the microchip and in a solution are quite similar. The influence of ionic surrounding has been studied in terms of the hybridization efficiency and discrimination between the mismatched and perfect duplexes. Different approaches have been tested to compensate the dependence of duplex stability on the GC content. It has been demonstrated that the use of chaotropic agents, addition of nonlabeled GC-rich competitor oligonucleotides, as well as creation of a temperature gradient along the microchip reproducing the distribution of melting temperatures, efficiently level out the AT/GC differences. The use of tetramethylammonium chloride for the same purpose was accompanied by weakening to some extent the discrimination between the mismatched duplexes and the perfect ones.     

Hydration pattern of A4T4 and T4A4 DNA: a molecular dynamics study

Hydration pattern and energetics of 'A-tract' containing duplexes have been studied using molecular dynamics on 12-mer self-complementary sequences 5'-d(GCA4T4GC)-3' and 5'-d(CGT4A4CG)-3'. The structural features for the simulated duplexes showed correlation with the corresponding experimental structures. Analysis of the hydration pattern confirmed that water network around the simulated duplexes is more conformation specific rather than sequence specific. The calculated heat capacity change upon duplex formation showed that the process is entropically driven for both the sequences. Furthermore, the theoretical free energy estimates calculated using MMPBSA approach showed a higher net electrostatic contribution for A4T4 duplex formation than for T4A4, however, energetically both the duplexes are observed to be equally stable.     

Long-time stretched exponential kinetics in single DNA duplex dissociation

We probe DNA hybridization kinetics by measuring the lifetime distribution of single 16-bp duplexes under thermal dissociation. Our unique approach, based on two DNA-coated microspheres in an extended optical tweezer, allows the study of single duplex DNA molecules under negligible molecular tension. In contrast to earlier experiments, we find a stretched exponential lifetime distribution, which is likely due to dissociation proceeding via a number of competing pathways between highly force-sensitive intermediate states. Similar measurements of microspheres linked by multiple DNA bridges find they have unexpected short bound lifetimes, also consistent with force sensitivity.     

Use of hybridization melting kinetics for detecting Phytophthora species using three-dimensional microarrays: demonstration of a novel concept for the differentiation of detection targets

Microarray-based detection is limited by variable and inconsistent hybridization intensities across the diversity of probes used in each array. In this paper, we introduce a novel concept for the differentiation of detection targets using duplex melting kinetics. A microarray assay was developed on a PamChip microarray enabling the differentiation of target Phytophthora species using the melting kinetics of probe-target duplexes. In the majority of cases the hybridization kinetics of target and non-target duplexes differed significantly. Analysis of the melting kinetics of duplexes formed by probes with target and non-target DNA was found to be an effective method for determining specific hybridization and was independent of fluctuations in hybridization signal intensity. This form of analysis was more robust than the traditional approach based on hybridization intensity, and enabled the detection of individual Phytophthora species and mixtures thereof.     

Considering the oligonucleotides secondary structures at thermodynamic and kinetic analysis of the DNA-duplexes formation

The influence of secondary structures of DNA oligonucleotides on thermodynamics and kinetics at the formation of their bimolecular complexes (duplexes) has been studied. The models considering inherent secondary structures of duplex components and their influence on quantitative thermodynamic and kinetic characteristics of the duplexes have been developed. The values of thermodynamic impacts given by individual structural elements of the double helix have been shown to depend on hairpin structuring of the duplex components. The "concentration" method to consider oligonucleotides intramolecular structure with thermodynamic parameters of bimolecular duplex formation has been proposed. According to stop-flow measurements, the observed values of association and dissociation constants are influenced by the presence of inherent structures in duplex components. The influence observed is increased with the lowering of the sample temperature. The analysis of experimental data involving the developed models provides the possibility to determine "proper" kinetic constants for the helix-to-coil transition. The difference between observed and calculated rate constants can amount up to two or more orders of magnitude.     

N(4)C-alkyl-N(4)C cross-linked DNA: bending deformations in duplexes that contain a -CNG- interstrand cross-link

Short DNA duplexes containing a 1,3-N(4)C-alkyl-N(4)C interstrand cross-link that joins the two C residues of a -CNG- sequence were prepared using either a phosphoramidite or convertible nucleoside approach. The alkyl cross-link consists of 2, 4, or 7 methylene groups. The duplexes, which contain a seven-base-pair core and A(3)/T(3) complementary 3'-overhanging ends, were characterized by enzymatic digestion and MALDI-TOF mass spectrometry. Ultraviolet thermal denaturation studies showed that the duplexes denature in a cooperative manner and that the length of the cross-link affects the thermal stability. Thus, the transition temperature of the ethyl cross-linked duplex, 42 degrees C, is 16 degrees C higher than the melting temperature of the corresponding non-cross-linked control, whereas the transition temperatures of the butyl and heptyl cross-linked duplexes, 73 and 72 degrees C, respectively, are 46-47 degrees C higher. The reduced molecularity of denaturation of the cross-linked duplexes versus melting of the non-cross-linked duplex most likely accounts for these differences. Examination of molecular models suggests that the ethyl cross-link is too short to span the distance between the two C residues at the site of the cross-link in B-form DNA without causing distortion of the helix, whereas less and no distortion would be expected for the butyl and heptyl cross-links, respectively. The circular dichroism spectra, which show greatest deviation in the ethyl cross-linked duplex from B-form DNA, are consistent with this expectation. Anomalous mobilities on native polyacrylamide gels of multimers produced by self-ligation of each of the cross-linked duplexes suggest that the ethyl and butyl cross-linked duplexes undergo bending deformations, whereas multimers derived from the heptyl cross-linked duplex migrated normally. The bending angle was estimated to be 20 degrees, 13 degrees, and 0 degrees for the ethyl, butyl, and heptyl cross-linked duplexes, respectively. Thus, it appears that the degree of bending in these N(4)C-alkyl-N(4)C cross-linked duplexes is controlled by the length of the cross-link.