Long dwell-time passage of DNA through nanometer-scale pores: kinetics and sequence dependence of motion

A detailed understanding of the kinetics of DNA motion though nanometer-scale pores is important for the successful development of many of the proposed next-generation rapid DNA sequencing and analysis methods. Many of these approaches require DNA motion through nanopores to be slowed by several orders of magnitude from its native translocation velocity so that the translocation times for individual nucleotides fall within practical timescales for detection. With the increased dwell time of DNA in the pore, DNA-pore interactions begin to play an increasingly important role in translocation kinetics. In previous work, we and others observed that when the DNA dwell time in the pore is substantial (>1 ms), DNA motion in α-hemolysin (α-HL) pores leads to nonexponential kinetics in the escape of DNA out of the pore. Here we show that a three-state model for DNA escape, involving stochastic binding interactions of DNA with the pore, accurately reproduces the experimental data. In addition, we investigate the sequence dependence of the DNA escape process and show that the interaction strength of adenine with α-HL is substantially lower relative to cytosine. Our results indicate a difference in the process by which DNA moves through an α-HL nanopore when the motion is fast (microsecond timescale) as compared with when it is slow (millisecond timescale) and strongly influenced by DNA-pore interactions of the kind reported here. We also show the ability of wild-type α-HL to detect and distinguish between 5-methylcytosine and cytosine based on differences in the absolute ionic current through the pore in the presence of these two nucleotides. The results we present here regarding sequence-dependent (and dwell-time-dependent) DNA-pore interaction kinetics will have important implications for the design of methods for DNA analysis through reduced-velocity motion in nanopores.     

Competitive displacement of DNA during surface hybridization

Using real-time dual-color fluorescence detection, we have experimentally tracked individual target species during competitive DNA surface hybridization in a two-component sample. Our experimental results demonstrate displacement of the lower affinity species by the higher affinity species and corroborate recent theoretical models describing competitive DNA surface hybridization. Competition at probe sites complementary to one of the two DNA species was monitored in separate experiments for two different target pairs. Each pair differs in sequence by a single nucleotide polymorphism, and one pair includes a folding target. We propose a mechanistic interpretation of the differences between hybridization curves of targets in multi-component and single-component experiments.     

DNA interstrand crosslink formation by mechlorethamine at a cytosine-cytosine mismatch pair: kinetics and sequence dependence

Expansion of the triplet repeat DNA sequence d[CGG]n.d[CCG]n is a characteristic of Fragile X syndrome, a human neurodegenerative disease. Stable intrastrand conformations formed by both d[CGG]n and d[CCG]n, and involving G-G and C-C mismatch pairs, respectively, are believed to be of importance in the development of the disease. We have shown previously that C-C mismatch pairs can be crosslinked covalently by mechlorethamine, a nitrogen mustard alkylating agent, and hence this reaction may be of value as a probe for conformers of d[CCG]n. To characterize the mechlorethamine C-C crosslink reaction further, here we report the kinetics and sequence dependence of formation of the crosslink species, using a series of model duplexes. The rate of reaction depends on the base sequence proximal to the C-C mismatch pair. Hence, in 19mer duplexes containing a central d[M4M3M2M1Cn1n2n3n4].d[N4N3N2N1Cm1m2m3m4] sequence, where M-m and N-n are complementary base pairs, the amount of crosslink increased with increasing G-C content of the eight base pairs neighboring the C-C mismatch and with the proximity of the G-C pairs to the C-C mismatch. Molecular dynamics simulations of the solvated duplexes provided an explanation of these data. Hence, for a C-C pair flanked by G-C base pairs the mismatched cytosine bases remain stacked within the duplex, but for a C-C pair flanked by A-T base pairs, the simulations suggested local opening of the duplex around the C-C pair, making it a less effective target for mechlorethamine.     

DNA reassociation kinetics and hybridization in different angiosperms

Reassociation kinetics ofDaucus carota andPetroselinum crispum (Apiaceae), andDatura innoxia (Solanaceae) are presented. Hybridization of³H-labelled DNA of two carrot cultivars indicate strong qualitative homologies of DNA sequences; nevertheless, certain quantitative differences in some Cotregions seem to exist. However, homologous sequences ofDaucus DNA with DNA ofDatura, and, suprisingly, even with DNA ofPetroselinum are very restricted: between 8% in the repeated regions and ca. 7–9% in the unique regions.     

Modeling of DNA hybridization kinetics for spatially resolved biochips

The marriage of microfluidics with detection technologies that rely on highly selective nucleic acid hybridization will provide improvements in bioanalytical methods for purposes such as detection of pathogens or mutations and drug screening. The capability to deliver samples in a controlled manner across a two-dimensional hybridization detection platform represents a substantial technical challenge in the development of quantitative and reusable biochips. General theoretical and numerical models of heterogeneous hybridization kinetics are required in order to design and optimize such biochips and to develop a quantitative method for online interpretation of experimental results. In this work we propose a general kinetic model of heterogeneous hybridization and develop a technique for estimating the kinetic coefficients for the case of well-spaced, noninteracting surface-bound probes. The experimentally verified model is then incorporated into the BLOCS (biolab-on-a-chip simulation) 3D microfluidics finite element code and used to model the dynamic hybridization on a biochip surface in the presence of a temperature gradient. These simulations demonstrate how such a device can be used to discriminate between fully complementary and single-base-pair mismatched hybridization using fluorescence detection by interpretation of the unique spatially resolved intensity pattern. It is also shown how the dynamic transport of the targets is likely to affect the rate and location of hybridization as well as that, although nonspecific hybridization is present, the change in the concentration of hybridized targets over the sensor platform is sufficiently high to determine if a fully complementary match is present. Practical design information such as the optimum transport speed, target concentration, and channel height is presented. The results presented here will aid in the interpretation of results obtained with such a temperature-gradient biochip.     

Origin of DNA helical structure and its sequence dependence

Conformational analysis of DNA shows that the origin of the B-form double helix can be attributed in large part to the atomic charge pattern in the base pairs. The charge patterns favor specific helical stacking of the base pairs. Base pairs alone--without backbones--have a strong tendency to form helix, indicating that the backbones play a rather passive role in determining the basic helical structure of DNA. It is mainly the electrostatic interactions determined by the charge pattern on base pairs that stabilize a particular helical conformation. The charge pattern in the base pairs appears to be responsible for much of the sequence dependence of DNA conformation, rather than steric clashes.     

Optical study of DNA surface hybridization reveals DNA surface density as a key parameter for microarray hybridization kinetics

We investigate the kinetics of DNA hybridization reactions on glass substrates, where one 22 mer strand (bound-DNA) is immobilized via phenylene-diisothiocyanate linker molecule on the substrate, the dye-labeled (Cy3) complementary strand (free-DNA) is in solution in a reaction chamber. We use total internal reflection fluorescence for surface detection of hybridization. As a new feature we perform a simultaneous real-time measurement of the change of free-DNA concentration in bulk parallel to the total internal reflection fluorescence measurement. We observe that the free-DNA concentration decreases considerably during hybridization. We show how the standard Langmuir kinetics needs to be extended to take into account the change in bulk concentration and explain our experimental results. Connecting both measurements we can estimate the surface density of accessible, immobilized bound-DNA. We discuss the implications with respect to DNA microarray detection.     

Secondary structure effects on DNA hybridization kinetics: a solution versus surface comparison

The hybridization kinetics for a series of designed 25mer probe–target pairs having varying degrees of secondary structure have been measured by UV absorbance and surface plasmon resonance (SPR) spectroscopy in solution and on the surface, respectively. Kinetic rate constants derived from the resultant data decrease with increasing probe and target secondary structure similarly in both solution and surface environments. Specifically, addition of three intramolecular base pairs in the probe and target structure slow hybridization by a factor of two. For individual strands containing four or more intramolecular base pairs, hybridization cannot be described by a traditional two-state model in solution-phase nor on the surface. Surface hybridization rates are also 20- to 40-fold slower than solution-phase rates for identical sequences and conditions. These quantitative findings may have implications for the design of better biosensors, particularly those using probes with deliberate secondary structure.     

Drosophila Rrp1 3'-exonuclease: demonstration of DNA sequence dependence and DNA strand specificity.

Drosophila Rrp1 (recombination repair protein 1) is a DNA repair enzyme whose nuclease activities include AP-endonuclease, 3'-exonuclease, 3'-phosphodiesterase and 3'-phosphatase. This study investigates the sequence specificity of the dsDNA 3'-exonuclease activity of Rrp1. We demonstrate that the activity is more efficient in purine-rich regions of dsDNA than in pyrimidine-rich regions. Rrp1 exonuclease activity is examined at 3'-terminal homopurine or homopyrimidine tracts, at junctions between purine- and pyrimidine-rich sequences and upon encountering repeated dinucleotide runs. The data show that purine-purine and 3'-pyrimidine-5'-purine dinucleotide bonds are cleaved faster than 3'-purine-5'-pyrimidine or pyrimidine-pyrimidine bonds. Thus, the base occupying the penultimate position in the 3'-terminal dinucleotide may be important in determining the relative efficiency of bond cleavage by Rrp1. These findings may reflect upon specific DNA-protein interactions in the enzyme active site.     

On translocation through a membrane channel via an internal binding site: kinetics and voltage dependence

Here we present a model for maltodextrin translocation through maltoporin channels. In a first step, our theoretical analysis does consider the case of a single binding site for a given substrate in a structurally unaffected channel with a possibly different entrance barrier on either side. It is shown how by means of conventional electrical conductance measurements (including current noise analysis) the basic equilibrium and rate constants can be determined as functions of the applied voltage. Then also the net translocation rate of the substrate becomes accessible quantitatively. This most simple model mechanism has been extended to include a voltage-dependent fast conformational change of the channel that prevents the binding process. The so developed approach has been tested with experimental data for a single maltoporin trimer being reconstituted in black lipid membranes when studied in the presence of maltohexaose as the substrate. The experimental results turned out to be clearly incompatible with binding alone. They are, however, very satisfactorily fitted by pertinent theoretical curves if also inhibition of binding by a conformational transition is taken into account. Accordingly, quantitative evaluations of the underlying parameters and eventually of the translocation rate have been carried out successfully. Our analysis reveals a set of parameters necessary for an optimal translocation that nicely corresponds to natural conditions.     

Nuclear ribosomal DNA sequence polymorphism and hybridization in checker mallows (Sidalcea, Malvaceae)

Checker mallows (Sidalcea, Malvaceae) constitute a western North American genus of annuals and perennials that have been regarded as taxonomically difficult because of complex patterns of morphological variation putatively stemming from hybridization and polyploidy. In recent molecular phylogenetic investigations extensive polymorphism was observed in the internal and external transcribed spacers (ITS and ETS) of 18S-26S nuclear ribosomal DNA for some Sidalcea samples. To resolve the evolutionary basis for this polymorphism and to readdress the evolutionary impact of hybridization in Sidalcea we cloned and sequenced the polymorphic DNAs and included the clones in phylogenetic analyses together with direct sequences of non-polymorphic samples. The positions of cloned spacer sequences in the phylogenetic trees suggest that S. reptans and two subspecies of S. malviflora may have been influenced by past hybridization with lineages of the "glaucescens" clade. Polymorphic sequence patterns in other taxa may be a result of extensive interbreeding within young clades, in keeping with the minimal sequence divergence, largely overlapping geographic distributions and morphology, and ploidy variation in these groups. Other possible explanations for polymorphic sequences in members of Sidalcea include slow concerted evolution relative to mutation rates, incomplete lineage sorting, and recent pseudogene formation.     

Influence of secondary structure on kinetics and reaction mechanism of DNA hybridization

Hybridization of nucleic acids with secondary structure is involved in many biological processes and technological applications. To gain more insight into its mechanism, we have investigated the kinetics of DNA hybridization/denaturation via fluorescence resonance energy transfer (FRET) on perfectly matched and single-base-mismatched DNA strands. DNA hybridization shows non-Arrhenius behavior. At high temperature, the apparent activation energies of DNA hybridization are negative and independent of secondary structure. In contrast, when temperature decreases, the apparent activation energies of DNA hybridization change to positive and become structure dependent. The large unfavorable enthalpy of secondary structure melting is compensated for by concomitant duplex formation. Based on our results, we propose a reaction mechanism about how the melting of secondary structure influences the hybridization process. A significant point in the mechanism is that the rate-limiting step switches along with temperature variation in the hybridization process of structured DNA, because the free energy profile of hybridization in structured DNA varies with the variation in temperature.     

Nucleotide-binding kinetics of Na,K-ATPase: cation dependence

Correlation between the Na,K-ATPase affinity to ADP and the cation (its nature and concentration) present in the medium was investigated. In buffer with low ionic strength (I approximately 1 mM) high-affinity ADP binding was not observed, while a stepwise increase in the concentrations of added cation (Na(+), Tris(+), imidazole(+), N-methylglucamine(+), choline(+)) induced an increase in the ADP affinity. The effect was fully saturated at 30-50 mM for all of the cations tested. The maximal affinity for ADP was slightly higher in the presence of Na(+), Tris(+), or imidazole(+) than in the presence of N-methylglucamine(+) or choline(+) (equilibrium dissociation constant K(d) 0.2-0.3 vs 0.7 microM). The ADP dissociation rates from its complex with enzyme in the presence of Na(+) or Tris(+) were similar, implying identity of the nucleotide-binding enzyme conformations, which therefore are assigned to E(1). The ability to compete with K(+) clearly distinguished Na(+) from other cations, which speaks against the sole involvement of the transport sites in the induction of the ADP-binding E(1) conformation. Since the cations are similar in their mode of induction of the high ADP affinity but they demonstrate a pronounced difference in ability to compete with K(+), their effects cannot be combined within any scheme with only one type of cation-binding sites. We suggest that the high affinity toward nucleotide is induced by cation interactions within the protein or lipid and that these nucleotide-domain-related sites coexist with the transport sites, which bind only Na(+) or K(+).     

Steady-state opticohydrodynamic properties of DNA: Molecular weight dependence and the internal viscosity problem

Extinction angles, flow birefringence, and intrinsic viscosities are compared for linear, bihelical DNAs from viral and other sources that span a range in molecular weight from ～10⁵ to 1.3 × 10⁸. This range effectively spans the region over which transition from rigid-rod to expanded-coil hydrodynamic property behavior occurs. All DNAs are in identical, phosphate–EDTA, neutral-pH buffers, 0.1M in NaCl. The extinction angle is a hydrodynamic property only and is thus particularly sensitive to effects of kinetic chain rigidity or internal viscosity. Our extinction angle results cannot be interpreted by any simple, single-function theoretical expression. Rather, they must be divided into distinct high- and low-molecular-weight domains. The low-molecular-weight region is typical of rigid-particle opticohydrodynamic property behaviour characterized primarily by particle orientation. The high-molecular-weight domain shows evidence for a finite internal viscosity effect, however, which can be interpreted as very nearly Kuhnian using Cerf's amplification of the Gaussian subchain model to include internal viscosity. It is found that the high-molecular-weight, monodisperse viral DNAs from T7, T5, and T2 bacteriophage show an internal viscosity contribution to the limiting extinction angle–shear rate ratio of ～3 × 10^?3 s. An effect of this magnitude may be marginally important in interpreting extinction angle and certain other hydrodynamic property data for high-molecular-weight DNA systems. Internal viscosity effects do not appear to be manifest in the ratio of flow birefringence to intrinsic viscosity, however, and the persistence length of the high-molecular-weight DNAs is found to be independent of molecular weight to within estimated experimental uncertainty.     

Efficient identification of DNA hybridization partners in a sequence database

MOTIVATION: The specific hybridization of complementary DNA molecules underlies many widely used molecular biology assays, including the polymerase chain reaction and various types of microarray analysis. In order for such an assay to work well, the primer or probe must bind to its intended target, without also binding to additional sequences in the reaction mixture. For any given probe or primer, potential non-specific binding partners can be identified using state-of-the-art models of DNA binding stability. Unfortunately, these models rely on dynamic programming algorithms that are too slow to apply on a genomic scale. RESULTS: We present an algorithm that efficiently scans a DNA database for short (approximately 20-30 base) sequences that will bind to a query sequence. We use a filtering approach, in which a series of increasingly stringent filters is applied to a set of candidate k-mers. The k-mers that pass all filters are then located in the sequence database using a precomputed index, and an accurate model of DNA binding stability is applied to the sequence surrounding each of the k-mer occurrences. This approach reduces the time to identify all binding partners for a given DNA sequence in human genomic DNA by approximately three orders of magnitude, from two days for the ENCODE regions to less than one minute for typical queries. Our approach is scalable to large DNA sequences. Our method can scan the human genome for medium strength binding sites to a candidate PCR primer in an average of 34.5 minutes. AVAILABILITY: Software implementing the algorithms described here is available at http://noble.gs.washington.edu/proj/dna-binding.     

Radiation-induced DNA breaks in different human satellite DNA sequence areas,analyzed by DNA breakage detection-fluorescence in situ hybridization

Human blood leukocytes were exposed to X rays to analyze the initial level of DNA breakage induced within different satellite DNA sequence areas and telomeres, using the DNA breakage detection-FISH procedure. The satellite DNA families analyzed comprised alphoid sequences, satellite 1, and 5-bp classical satellite DNA sequences from chromosome 1 (D1Z1 locus), from chromosome 9 (D9Z3 locus), and from the Y chromosome (DYZ1 locus). Since the control hybridization signal was quite different in each of the DNA targets, the relative increase in whole fluorescence intensity with respect to unirradiated controls was the parameter used for comparison. Irradiation of nucleoids obtained after protein removal demonstrated that the alkaline unwinding solution generates around half the amount of signal when breaks are present in the 5-bp classical DNA satellites as when the same numbers of breaks are present the genome overall, whereas the signal is slightly stronger when the breaks are within the alphoids or satellite 1 sequences. After correction for differences in sensitivity to the alkaline unwinding-renaturation, DNA housed in chromatin corresponding to 5-bp classical satellites proved to be more sensitive to breakage than the overall genome, whereas DNA in the chromatin corresponding to alphoids or satellite 1 showed a sensitivity similar to that of the whole genome. The minimum detectable dose was 0.1 Gy for the whole genome, 0.2 Gy for alphoids and satellite 1, and 0.4 Gy for the 5-bp classical satellites. Telomeric DNA sequences appeared to be maximally labeled in unirradiated cells. Thus telomeric ends behave like DNA breaks, constituting a source of background in alkaline unwinding assays.     

The effect of sequence heterogeneity on DNA melting kinetics

We consider kinetics of the cooperative melting of DNA sections situated at the edge of the helix. Accurate calculations based on the real sequences of such sections demonstrate that their internal heterogeneity has a drastic effect on the melting kinetics. Allowance for the internal heterogeneity increases the relaxation time by several orders of magnitude as compared with a model based on the assumption of equal base-pair stability within a section. The relaxation times obtained are in good agreement with the experimental data of Suyama and Wada (A. Suyama and A. Wada, Biopolymers, 23, 409 (1984)). An analysis of the melting process revealed some simple sequence characteristics that determine its rate. An examination of the temperature dependence of the relaxation time led to a distinct interpretation of the apparent activation energies of the denaturation and renaturation. The relaxation time proved to reach its maximum near the equilibrium melting point of the section examined.     

Rapid hybridization kinetics of DNA attached to submicron latex particles.

We describe a novel method for attaching any DNA molecule to submicron latex beads and characterize the hybridization kinetic properties of these bead-DNA conjugates. The conjugates hybridize to DNA in solution with rates comparable to homogeneous hybridization reactions, are compatible with common hybridization conditions and are conveniently manipulated. They should thus serve as useful reagents for the fractionation and characterization of DNA and RNA.