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.     

Quantitative analysis of DNA hybridization in a flowthrough microarray for molecular testing

Quantitative information about the nucleic acids hybridization reaction on microarrays is fundamental to designing optimized assays for molecular diagnostics. This study presents the kinetic, equilibrium, and thermodynamic analyses of DNA hybridization in a microarray system designed for fast molecular testing of pathogenic bacteria. Our microarray setup uses a porous, nylon membrane for probe immobilization and flowthrough incubation. The Langmuir model was used to determine the reaction rate constants of hybridization with antisense targets specific to Staphylococcus epidermidis and Staphylococcus aureus strains. The kinetic analysis revealed a sequence-dependent reaction rate, with association rate constants on the order of 10⁵ M⁻¹ s⁻¹ and dissociation rate constants of 10⁻⁴ s⁻¹. We found that by increasing the probe surface density from 10¹¹ to 10¹² molecules/cm², the hybridization rate and efficiency are suppressed while the melting temperature of the DNA duplex increases. The maximum fraction of hybridized capture probes at equilibrium did not exceed 50% for hybridization with antisense sequences and was below 6% for hybridization with long targets obtained from PCR. The van’t Hoff analysis of the temperature denaturation data showed that the DNA hybridization in our porous, flowthrough microarray is thermodynamically less favorable than the hybridization of the same sequences in solution.     

The biophysics of DNA hybridization with immobilized oligonucleotide probes.

A mathematical model based on receptor-ligand interactions at a cell surface has been modified and further developed to represent heterogeneous DNA-DNA hybridization on a solid surface. The immobilized DNA molecules with known sequences are called probes, and the DNA molecules in solution with unknown sequences are called targets in this model. Capture of the perfectly complementary target is modeled as a combined reaction-diffusion limited irreversible reaction. In the model, there are two different mechanisms by which targets can hybridize with the complementary probes: direct hybridization from the solution and hybridization by molecules that adsorb nonspecifically and then surface diffuse to the probe. The results indicate that nonspecific adsorption of single-stranded DNA on the surface and subsequent two-dimensional diffusion can significantly enhance the overall reaction rate. Heterogeneous hybridization depends strongly on the rate constants for DNA adsorption/desorption in the non-probe-covered regions of the surface, the two-dimensional (2D) diffusion coefficient, and the size of probes and targets. The model shows that the overall kinetics of DNA hybridization to DNA on a solid support may be an extremely efficient process for physically realistic 2D diffusion coefficients, target concentrations, and surface probe densities. The implication for design and operation of a DNA hybridization surface is that there is an optimal surface probe density when 2D diffusion occurs; values above that optimum do not increase the capture rate. Our model predicts capture rates in agreement with those from recent experimental literature. The results of our analysis predict that several things can be done to improve heterogeneous hybridization: 1) the solution phase target molecules should be about 100 bases or less in size to speed solution-phase and surface diffusion; 2) conditions should be created such that reversible adsorption and two-dimensional diffusion occur in the surface regions between DNA probe molecules; 3) provided that 2) is satisfied, one can achieve results with a sparse probe coverage that are equal to or better than those obtained with a surface totally covered with DNA probes.     

A high throughput method to investigate oligodeoxyribonucleotide hybridization kinetics and thermodynamics.

We describe a high throughput microtiter-based assay to measure binding of oligodeoxyribonucleotides to nucleic acid targets. The assay utilizes oligodeoxyribonucleotide probes labeled with a highly chemiluminescent acridinium ester (AE). Reaction of AE with sodium sulfite renders it non-chemiluminescent. When an AE-labeled probe hybridizes to a target nucleic acid AE is protected from reaction with sodium sulfite and thus remains chemiluminescent. In contrast, unhybridized probe readily reacts with sodium sulfite and is rendered non-chemiluminescent. Hybridization of an AE-labeled probe to a target nucleic acid can therefore be detected without physical separation of unhybridized probe by treatment of the hybridization reaction with sodium sulfite and measurement of the remaining chemiluminescence. Using this method we measured hybridization rate constants and thermodynamic affinities of oligodeoxyribonucleotide probes binding to simple synthetic targets as well as large complex biological targets. The kinetic and thermodynamic parameters were measured with a high degree of accuracy and were in excellent agreement with values measured by other established techniques.     

Compressive Sensing DNA Microarrays

Compressive sensing microarrays (CSMs) are DNA-based sensors that operate using group testing and compressive sensing (CS) principles. In contrast to conventional DNA microarrays, in which each genetic sensor is designed to respond to a single target, in a CSM, each sensor responds to a set of targets. We study the problem of designing CSMs that simultaneously account for both the constraints from CS theory and the biochemistry of probe-target DNA hybridization. An appropriate cross-hybridization model is proposed for CSMs, and several methods are developed for probe design and CS signal recovery based on the new model. Lab experiments suggest that in order to achieve accurate hybridization profiling, consensus probe sequences are required to have sequence homology of at least 80% with all targets to be detected. Furthermore, out-of-equilibrium datasets are usually as accurate as those obtained from equilibrium conditions. Consequently, one can use CSMs in applications in which only short hybridization times are allowed.     

Hybridization kinetics and thermodynamics of molecular beacons

Molecular beacons are increasingly being used in many applications involving nucleic acid detection and quantification. The stem–loop structure of molecular beacons provides a competing reaction for probe–target hybridization that serves to increase probe specificity, which is particularly useful when single-base discrimination is desired. To fully realize the potential of molecular beacons, it is necessary to optimize their structure. Here we report a systematic study of the thermodynamic and kinetic parameters that describe the molecular beacon structure–function relationship. Both probe and stem lengths are shown to have a significant impact on the binding specificity and hybridization kinetic rates of molecular beacons. Specifically, molecular beacons with longer stem lengths have an improved ability to discriminate between targets over a broader range of temperatures. However, this is accompanied by a decrease in the rate of molecular beacon–target hybridization. Molecular beacons with longer probe lengths tend to have lower dissociation constants, increased kinetic rate constants, and decreased specificity. Molecular beacons with very short stems have a lower signal-to-background ratio than molecular beacons with longer stems. These features have significant implications for the design of molecular beacons for various applications.     

Computer simulation study of molecular recognition in model DNA microarrays

DNA microarrays have been widely adopted by the scientific community for a variety of applications. To improve the performance of microarrays there is a need for a fundamental understanding of the interplay between the various factors that affect microarray sensitivity and specificity. We use lattice Monte Carlo simulations to study the thermodynamics and kinetics of hybridization of single-stranded target genes in solution with complementary probe DNA molecules immobilized on a microarray surface. The target molecules in our system contain 48 segments and the probes tethered on a hard surface contain 8-24 segments. The segments on the probe and target are distinct and each segment represents a sequence of nucleotides ( approximately 11 nucleotides). Each probe segment interacts exclusively with its unique complementary target segment with a single hybridization energy; all other interactions are zero. We examine how the probe length, temperature, or hybridization energy, and the stretch along the target that the probe segments complement, affect the extent of hybridization. For systems containing single probe and single target molecules, we observe that as the probe length increases, the probability of binding all probe segments to the target decreases, implying that the specificity decreases. We observe that probes 12-16 segments ( approximately 132-176 nucleotides) long gave the highest specificity and sensitivity. This agrees with the experimental results obtained by another research group, who found an optimal probe length of 150 nucleotides. As the hybridization energy increases, the longer probes are able to bind all their segments to the target, thus improving their specificity. The hybridization kinetics reveals that the segments at the ends of the probe are most likely to start the hybridization. The segments toward the center of the probe remain bound to the target for a longer time than the segments at the ends of the probe.     

Real-time multiplex PCR assays

The ability to multiplex PCR by probe color and melting temperature (T(m)) greatly expands the power of real-time analysis. Simple hybridization probes with only a single fluorescent dye can be used for quantification and allele typing. Different probes are labeled with dyes that have unique emission spectra. Spectral data are collected with discrete optics or dispersed onto an array for detection. Spectral overlap between dyes is corrected by using pure dye spectra to deconvolute the experimental data by matrix algebra. Since fluorescence is temperature dependent and depends on the dye, spectral overlap and color compensation constants are also temperature dependent. Single-labeled probes are easier to synthesize and purify than more complex probes with two or more dyes. In addition, the fluorescence of single-labeled probes is reversible and depends only on hybridization of the probe to the target, allowing study of the melting characteristics of the probe. Although melting curves can be obtained during PCR, data are usually acquired at near-equilibrium rates of 0.05-0.2 degrees C/s after PCR is complete. Using rapid-cycle PCR, amplification requires about 20 min followed by a 10-min melting curve, greatly reducing result turnaround time. In addition to dye color, melting temperature can be used for a second dimension of multiplexing. Multiplexing by color and T(m) creates a "virtual" two-dimensional multiplexing array without the need for an immobilized matrix of probes. Instead of physical separation along the X and Y axes, amplification products are identified by different fluorescence spectra and melting characteristics.     

Thermodynamically modulated partially double-stranded linear DNA probe design for homogeneous real-time PCR

Real-time PCR assays have recently been developed for diagnostic and research purposes. Signal generation in real-time PCR is achieved with probe designs that usually depend on exonuclease activity of DNA polymerase (e.g. TaqMan probe) or oligonucleotide hybridization (e.g. molecular beacon). Probe design often needs to be specifically tailored either to tolerate or to differentiate between sequence variations. The conventional probe technologies offer limited flexibility to meet these diverse requirements. Here, we introduce a novel partially double-stranded linear DNA probe design. It consists of a hybridization probe 5'-labeled with a fluorophore and a shorter quencher oligo of complementary sequence 3'-labeled with a quencher. Fluorescent signal is generated when the hybridization probe preferentially binds to amplified targets during PCR. This novel class of probe can be thermodynamically modulated by adjusting (i) the length of hybridization probe, (ii) the length of quencher oligo, (iii) the molar ratio between the two strands and (iv) signal detection temperature. As a result, pre-amplification signal, signal gain and the extent of mismatch discrimination can be reliably controlled and optimized. The applicability of this design strategy was demonstrated in the Abbott RealTime HIV-1 assay.     

Mutation detection by stacking hybridization on genosensor arrays

A new strategy for analysis of point mutations using oligonucleotide array (genosensor) hybridization was investigated. In the new approach, a single-stranded target strand is preannealed with a labeled "stacking oligonucleotide," and then the partially duplex labeled target molecule is hybridized to an array of glass-tethered oligonucleotide probes, targeted to the region on the target immediately adjacent to the stacking oligomer. In this configuration, the base-stacking interactions between the "capture probe" and the contiguously stacking oligomer stabilize the binding of the target molecule to its complementary probe on the genosensor array. The temperature of hybridization can be adjusted so that the target molecule will bind to the glass-tethered probe only in the presence of the stacking oligomer, and a single mismatch at or near the terminal position ol the capture probe disrupts the stacking interactions and thereby eliminates or greatly reduces the hybridization. This stacking hybridization approach was investigated using a collection of synthetic targets, probes, and stacking oligonucleotides, which permitted identification of conditions for optimal base mismatch discrimination. The oligonucleotide probes were tethered to the glass using a simple, improved attachment chemistry in which a 3'-aminopropanol function introduced into the probe during chemical synthesis binds covalently to silanol groups on clean, underivatized glass. "Operating parameters" examined in the stacking hybridization system included length of capture probe, position, type and number of mismatches between the probe and the target, temperature of hybridization and length of washing, and the presence of terminal phosphate group in the probe, at its junction with the stacking oligomer. The results suggest that in the stacking hybridization configuration: 1. Optimal mismatch discrimination with 9-mer probes occurs at 45 degrees C, after which little or no improvement in mispair rejection occurred on lengthy continued washing at 45 degrees C. 2. At 25 degrees C optimal mismatch discrimination occurred with 7- or 8-mer probes, or with 9-mer probes containing an additional internal mismatch. 3. The presence of a phosphate group on the 5'-end of the glass-tethered probe had no general effect on mismatch discrimination, but influenced the relative stability of different mismatches in the sequence context studied. These results provide a motivation for continued development of the stacking hybridization technique for nucleic acid sequence analysis. This approach offers several advantages over the traditional allele-specific oligonucleotide hybridization technique, and is distinct from the contiguous stacking hybridization sitrategy that the Mirzabekov laboratory has introduced (Yershov et al. (1996) Proc. Natl. Acad. Sci. USA 93, 4913-4918; Parinov et al. (1996) Nucleic Acids Res. 24, 2998-3004).     

Rapid F508del and F508C assay using fluorescent hybridization probes

Amplification and fluorescent genotyping of the cystic fibrosis F508del locus was achieved from human genomic DNA in less than 30 min. The hybridization of adjacent fluorescent probes at the mutation site was monitored by resonance energy transfer between fluorescein and Cy5 during heating or cooling. Characteristic curves were obtained for each genotype; the first derivative of these fluorescent curves has a maximum at an apparent hybridization temperature (Tm) that is specific for each probe/allele duplex. The direction and rate of temperature change determines the difference between the apparent Tm and the true equilibrium Tm. One hundred and five sample were genotyped for the F508del cystic fibrosis mutation by heating and cooling curve profiles. These genotypes were validated by allele-specific amplification. Two fluorescein hybridization probes were designed to match the wild-type sequence perfectly from either codons 502 to 513 or from 504 to 511 on the cystic fibrosis transconductance regulator gene of chromosome 7. While genotyping for the F508del, an allele with the F508C base change was detected. For both F508del and F508C variants, the Tm shift from wild type was greater with a 24-mer probe than with a 35-mer probe. Fluorescent monitoring of hybridization probes is a versatile technique that can detect unexpected sequence alterations.     

Improved confocal characterization of targets in nuclei of cytogenetic preparations

OBJECTIVE: To show that cellular preparations requiring depth analysis of different domains stained by molecular cytogenetic methods (fluorescence in situ hybridization and primed in situ) can be improved by regularized factor analysis of medical image sequences (FAMIS) to isolate fluorescent probes by means of intensity depth profiles of fluorochromes, to track relevant DNA sequences (cosmids and centromeres) in cell nuclei during interphase and to improve the use of cytogenetic techniques resulting in flat preparations of whole cells that are assumed to preserve probe access to their targets. STUDY DESIGN: 3D sequences of images obtained by depth displacement in a confocal microscope were first analyzed by the FAMIS algorithm, which provides factor curves. Factor images then resulted from regularization methods that improve signal/noise ratio while preserving target contours. RESULTS: Factor curves and regularized factor images helped analyze targets inside nuclei. CONCLUSION: It is possible to process preparations containing numerous spots (even when they are on different planes) to differentiate stained targets, to investigate depth differences and to improve visualization and detection.     

Target concentration dependence of DNA melting temperature on oligonucleotide microarrays

The design of microarrays is currently based on studies focusing on DNA hybridization reaction in bulk solution. However, the presence of a surface to which the probe strand is attached can make the solution‐based approximations invalid, resulting in sub‐optimum hybridization conditions. To determine the effect of surfaces on DNA duplex formation, the authors studied the dependence of DNA melting temperature (T_m) on target concentration. An automated system was developed to capture the melting profiles of a 25‐mer perfect‐match probe–target pair initially hybridized at 23°C. Target concentrations ranged from 0.0165 to 15 nM with different probe amounts (0.03–0.82 pmol on a surface area of 10¹⁸ Ų), a constant probe density (5 × 10¹² molecules/cm²) and spacer length (15 dT). The authors found that T_m for duplexes anchored to a surface is lower than in‐solution, and this difference increases with increasing target concentration. In a representative set, a target concentration increase from 0.5 to 15 nM with 0.82 pmol of probe on the surface resulted in a T_m decrease of 6°C when compared with a 4°C increase in solution. At very low target concentrations, a multi‐melting process was observed in low temperature domains of the curves. This was attributed to the presence of truncated or mismatch probes. © 2012 American Institute of Chemical Engineers Biotechnol. Prog., 2012     

Salt Concentration Effects on Equilibrium Melting Curves from DNA Microarrays

DNA microarrays find applications in an increasing number of domains where more quantitative results are required. DNA being a charged polymer, the repulsive interactions between the surface of the microarray and the targets in solution are increasing upon hybridization. Such electrostatic penalty is generally reduced by increasing the salt concentration. In this article, we present equilibrium-melting curves obtained from dedicated physicochemical experiments on DNA microarrays in order to get a better understanding of the electrostatic penalty incurred during the hybridization reaction at the surface. Various salt concentrations have been considered and deviations from the commonly used Langmuir adsorption model are experimentally quantified for the first time in agreement with theoretical predictions.     

Determination of Herpes simplex virus DNA with the use of solid-phase methods for the detection of PCR products

To determine DNA of herpes simplex virus (HSV), types 1 and 2, the polymerase chain reaction (PCR) method was developed with the subsequent detection of amplification products by means of electrophoresis or the molecular hybridization of nucleic acids (MHNA). Two variants of MHNA have been compared: hybridization in the solution of a biotinylated probe with digoxigenin-labeled PCR with the subsequent sorption of hybridization complexes onto streptavidin-covered plates and solid-phase hybridization of digoxigenin-labeled PCR with a biotinylated probe. Effective hybridization was observed after the denaturation of targets at 95 degrees C in the solution of 50 mM NaOH, but not in neutral solutions. To increase the level of sensitivity of hybridization in solution, the exact selection of the amount of the probe was shown to be necessary, for both its excess and deficiency essentially decreased the method sensitivity. A decrease in the ionic power of hybridization solutions from 6 h SSC to 1 h SSC led to greater specificity of hybridization without a decrease in the method sensitivity. A rise in the temperature of hybridization and subsequent washing to 45 degrees C decreased the sensitivity of the method. The limit of the sensitivity PCR with electrophoretic detection was 30 HSV genome equivalents, and 10 genome equivalents in the presence MHNA in the solution and on the solid phase.     

Beyond Affymetrix arrays: expanding the set of known hybridization isotherms and observing pre-wash signal intensities

Microarray hybridization studies have attributed the nonlinearity of hybridization isotherms to probe saturation and post-hybridization washing. Both processes are thought to distort ‘true’ target abundance because immobilized probes are saturated with excess target and stringent washing removes loosely bound targets. Yet the paucity of studies aimed at understanding hybridization and dissociation makes it difficult to align physicochemical theory to microarray results. To fill the void, we first examined hybridization isotherms generated on different microarray platforms using a ribosomal RNA target and then investigated hybridization signals at equilibrium and after stringent wash. Hybridization signal at equilibrium was achieved by treating the microarray with isopropanol, which prevents nucleic acids from dissolving into solution. Our results suggest that (i) the shape of hybridization isotherms varied by microarray platform with some being hyperbolic or linear, and others following a power-law; (ii) at equilibrium, fluorescent signal of different probes hybridized to the same target were not similar even with excess of target and (iii) the amount of target removed by stringent washing depended upon the hybridization time, the probe sequence and the presence/absence of nonspecific targets. Possible physicochemical interpretations of the results and future studies are discussed.     

Matching base-pair number dependence of the kinetics of DNA-DNA hybridization studied by surface plasmon fluorescence spectroscopy

Two single-stranded DNA oligonucleotides consisting of complementary base-pairs can form double strands. This phenomenon is well studied in solutions, however, in order to clarify the physical mechanism of the hybridization occurring at a solid/solution interface, we studied the kinetics by surface plasmon fluorescence spectroscopy (SPFS): one single-stranded oligo-DNA (probe-DNA) was immobilized on the substrate, the other one (target-DNA) labelled with a fluorescent probe was added to the flow cell. After hybridization, the chromophores could be excited by the surface plasmon mode and their fluorescence detected with high sensitivity. The dependence of the k(on) and k(off) rate constants on the length of the hybridizing oligonucleotides was investigated by using a MM0 series (no mismatch) and the kinetics was found to be well described by a Langmuir adsorption model. From these measurements we found that also in the case of surface hybridization the affinity of the duplexes decreases as the number of matching base-pairs decreases from 15 to 10. In order to show that SPFS is the powerful technique with high sensitivity, the hybridization process for mixed target-oligos was measured by SPFS and analyzed by an expanded Langmuir model in which two components of target-oligo can bind to probe-DNA at the sensor surface competitively. Two sets of the k(on) and k(off) obtained from the experiment are successfully consistent with the k(on) and k(off) obtained from experiments for single (pure) target-DNA.