Fast high-resolution mapping of long fragments of genomic DNA based on single-molecule detection.

Here we describe bacterial genotyping by direct linear analysis (DLA) single-molecule mapping. DLA involves preparation of restriction digest of genomic DNA labeled with a sequence-specific fluorescent probe and stained nonspecifically with intercalator. These restriction fragments are stretched one by one in a microfluidic device, and the distribution of probes on the fragments is determined by single-molecule measurement of probe fluorescence. Fluorescence of the DNA-bound intercalator provides information on the molecule length. Because the probes recognize short sequences, they encounter multiple cognate sites on 100- to 300-kb-long DNA fragments. The DLA maps are based on underlying DNA sequences of microorganisms; therefore, the maps are unique for each fragment. This allows fragments of similar lengths that cannot be resolved by standard DNA sizing techniques to be readily distinguished. DNA preparation, data collection, and analysis can be carried out in as little as 5h when working with monocultures. We demonstrate the ability to discriminate between two pathogenic Escherichia coli strains, O157:H7 Sakai and uropathogenic 536, and we use DLA mapping to identify microorganisms in mixtures. We also introduce a second color probe to double the information used to distinguish molecules and increase the length range of mapped fragments.

An automated sample preparation system with mini-reactor to isolate and process submegabase fragments of bacterial DNA.

Existing methods for extraction and processing of large fragments of bacterial genomic DNA are manual, time-consuming, and prone to variability in DNA quality and recovery. To solve these problems, we have designed and built an automated fluidic system with a mini-reactor. Balancing flows through and tangential to the ultrafiltration membrane in the reactor, cells and then released DNA can be immobilized and subjected to a series of consecutive processing steps. The steps may include enzymatic reactions, tag hybridization, buffer exchange, and selective removal of cell debris and by-products of the reactions. The system can produce long DNA fragments (up to 0.5 Mb) of bacterial genome restriction digest and perform DNA tagging with fluorescent sequence-specific probes. The DNA obtained is of high purity and floating free in solution, and it can be directly analyzed by pulsed-field gel electrophoresis (PFGE) or used in applications requiring submegabase DNA fragments. PFGE-ready samples of DNA restriction digests can be produced in as little as 2.1 h and require less than 10(8) cells. All fluidic operations are automated except for the injection of the sample and reagents.

Base-stacking and base-pairing contributions into thermal stability of the DNA double helix.

Two factors are mainly responsible for the stability of the DNA double helix: base pairing between complementary strands and stacking between adjacent bases. By studying DNA molecules with solitary nicks and gaps we measure temperature and salt dependence of the stacking free energy of the DNA double helix. For the first time, DNA stacking parameters are obtained directly (without extrapolation) for temperatures from below room temperature to close to melting temperature. We also obtain DNA stacking parameters for different salt concentrations ranging from 15 to 100 mM Na+. From stacking parameters of individual contacts, we calculate base-stacking contribution to the stability of A*T- and G*C-containing DNA polymers. We find that temperature and salt dependences of the stacking term fully determine the temperature and the salt dependence of DNA stability parameters. For all temperatures and salt concentrations employed in present study, base-stacking is the main stabilizing factor in the DNA double helix. A*T pairing is always destabilizing and G*C pairing contributes almost no stabilization. Base-stacking interaction dominates not only in the duplex overall stability but also significantly contributes into the dependence of the duplex stability on its sequence.

Pseudocomplementary PNAs as selective modifiers of protein activity on duplex DNA: the case of type IIs restriction enzymes.

This study evaluates the potential of pseudocomplementary peptide nucleic acids (pcPNAs) for sequence-specific modification of enzyme activity towards double-stranded DNA (dsDNA). To this end, we analyze the ability of pcPNA-dsDNA complexes to site-selectively interfere with the action of four type IIs restriction enzymes. We have found that pcPNA-dsDNA complexes exhibit a different degree of DNA protection against cleaving/nicking activity of various isoschizomeric endonucleases under investigation (PleI, MlyI and N.BstNBI) depending on their type and mutual arrangement of PNA-binding and enzyme recognition/cleavage sites. We have also found that the pcPNA targeting to closely located PleI or BbsI recognition sites on dsDNA generates in some cases the nicking activity of these DNA cutters. At the same time, MlyI endonuclease, a PleI isoschizomer, does not exhibit any DNA nicking/cleavage activity, being completely blocked by the nearby pcPNA binding. Our results have general implications for effective pcPNA interference with the performance of DNA-processing proteins, thus being important for prospective applications of pcPNAs.

Stacked-unstacked equilibrium at the nick site of DNA.

Stability of duplex DNA with respect to separation of complementary strands is crucial for DNA executing its major functions in the cell and it also plays a central role in major biotechnology applications of DNA: DNA sequencing, polymerase chain reaction, and DNA microarrays. Two types of interaction are well known to contribute to DNA stability: stacking between adjacent base-pairs and pairing between complementary bases. However, their contribution into the duplex stability is yet to be determined. Now we fill this fundamental gap in our knowledge of the DNA double helix. We have prepared a series of 32, 300 bp-long DNA fragments with solitary nicks in the same position differing only in base-pairs flanking the nick. Electrophoretic mobility of these fragments in the gel has been studied. Assuming the equilibrium between stacked and unstacked conformations at the nick site, all 32 stacking free energy parameters have been obtained. Only ten of them are essential and they govern the stacking interactions between adjacent base-pairs in intact DNA double helix. A full set of DNA stacking parameters has been determined for the first time. From these data and from a well-known dependence of DNA melting temperature on G.C content, the contribution of base-pairing into duplex stability has been estimated. The obtained energy parameters of the DNA double helix are of paramount importance for understanding sequence-dependent DNA flexibility and for numerous biotechnology applications.

Single molecule DNA intercalation in continuous homogenous elongational flow.

Sequence-nonspecific staining of DNA with intercalating fluorophores is required for fluorescence-based length estimation of elongated DNA in optical mapping techniques. However, the observed length of a DNA molecule is affected by the relative concentrations of DNA and dye. In some applications, predetermination of DNA concentration may not be possible. Here we present a microfluidic approach in which individual DNA molecules are entrained by converging laminar sheath flows containing the intercalating dye PO-PRO-1. This provides uniform staining regardless of DNA concentration, and uniform elastic stretching of DNA in continuous elongational flow. On-chip intercalation provides a unique process for concentration-independent staining of long DNA fragments for the optical mapping method Genome Sequence Scanning (GSS), and normalizes intramolecular elasticity across a broad range of molecule lengths. These advances permit accurate mapping of observed molecules to sequence derived templates, thus improving detection of complex bacterial mixtures using GSS.

High-throughput genome scanning in constant tension fluidic funnels.

Genome Sequence Scanning (GSS) is a bacterial identification technology that detects sparse sequence-specific fluorescent tags on long DNA molecules linearized in a continuous-flow microfunnel. The efficiency and sensitivity of GSS depends on the detection throughput of well-stretched molecules. Previous studies have investigated the fundamental roles of elongational and shear flow on DNA stretching in continuous flow devices. Here, we present a novel compound funnel design that significantly improves effective molecule throughput over previously described devices. First, exploring the relationship between fluid strain rate and molecule tension allows for design of funnel geometries that accommodate high fluid velocities without compromising molecules due to overstretching under high tension. Second, a constant-strain detection channel is utilized to significantly reduce the number of molecules lost to further analysis due to shear-induced molecular tumbling and relaxation. Finally, application of the constant-strain detection channel allows for a priori prediction of spatial resolution bias due to accelerating flow. In all, the refined funnel geometries presented here yield over thirty-fold increase in effective molecule throughput due to increased fluid flow and improved retention of stretched molecules, compared to previously described devices.

Sequence-dependent base pair opening in DNA double helix.

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.

Monitoring of single nicks in duplex DNA by gel electrophoretic mobility-shift assay.

We demonstrate that the gel electrophoretic mobility-shift assay (EMSA) can be used for site-selective and quantitative monitoring of nicks in linear double-stranded DNA (dsDNA) thus allowing to expediently follow the nicking activity of enzymes or other agents targeted to a designated dsDNA site. At elevated temperature and/or in the presence of urea, DNA fragments carrying a single nick produced by the nicking enzyme N.BstNBI exhibit a well-detectable gel retardation effect. On the basis of permutation analysis, the decreased electrophoretic mobility of nicked dsDNA fragments is attributed to a bend (or hinge) in the DNA double helix sequence-specifically generated by a nick. Since nick-induced DNA bending depends on interaction between base pairs adjacent to a nick, the change in mobility is different for nicked DNA sites with different sequences. Therefore, EMSA monitoring of differential mobility change caused by nicks within various DNA sequences could be useful for studying the differential base stacking and nearest-neighbor energetics.

Specific versus nonspecific binding of cationic PNAs to duplex DNA.

Although peptide nucleic acids (PNAs) are neutral by themselves, they are usually appended with positively charged lysine residues to increase their solubility and binding affinity for nucleic acid targets. Thus obtained cationic PNAs very effectively interact with the designated duplex DNA targets in a sequence-specific manner forming strand-invasion complexes. We report on the study of the nonspecific effects in the kinetics of formation of sequence-specific PNA-DNA complexes. We find that in a typical range of salt concentrations used when working with strand-invading PNAs (10-20 mM NaCl) the PNA binding rates essentially do not depend on the presence of nontarget DNA in the reaction mixture. However, at lower salt concentrations (<10 mM NaCl), the rates of PNA binding to DNA targets are significantly slowed down by the excess of unrelated DNA. This effect of nontarget DNA arises from depleting the concentration of free PNA capable of interacting with DNA target due to adhesion of positively charged PNA molecules on the negatively charged DNA duplex. As expected, the nonspecific electrostatic effects are more pronounced for more charged PNAs. We propose a simple model quantitatively describing all major features of the observed phenomenon. This understanding is important for design of and manipulation with the DNA-binding polycationic ligands in general and PNA-based drugs in particular.

Tailoring the activity of restriction endonuclease PleI by PNA-induced DNA looping.

DNA looping is one of the key factors allowing proteins bound to different DNA sites to signal one another via direct contacts. We demonstrate that DNA looping can be generated in an arbitrary chosen site by sequence-directed targeting of double-stranded DNA with pseudocomplementary peptide-nucleic acids (pcPNAs). We designed pcPNAs to mask the DNA from cleavage by type IIs restriction enzyme PleI while not preventing the enzyme from binding to its primary DNA recognition site. Direct interaction between two protein molecules (one bound to the original recognition site and the other to a sequence-degenerated site) results in a totally new activity of PleI: it produces a nick near the degenerate site. The PNA-induced nicking efficiency varies with the distance between the two protein-binding sites in a phase with the DNA helical periodicity. Our findings imply a general approach for the fine-tuning of proteins bound to DNA sites well separated along the DNA chain.

Kinetics and mechanism of the DNA double helix invasion by pseudocomplementary peptide nucleic acids.

If adenines and thymines in two mutually complementary mixed-base peptide nucleic acid (PNA) oligomers are substituted with diaminopurines and thiouracils, respectively, so-called pseudocomplementary PNAs (pcPNAs) are created. Pairs of pcPNAs have recently demonstrated an ability to highly selectively target essentially any designated site on double-stranded DNA (dsDNA) by forming very stable PNA-DNA strand-displacement complexes via double duplex invasion (helix invasion). These properties of pcPNAs make them unique and very promising ligands capable of denying the access of DNA-binding proteins to dsDNA. To elucidate the sequence-unrestricted mechanism of sequence-specific dsDNA recognition by pcPNAs, we have studied the kinetics of formation of corresponding PNA-DNA complexes at various temperatures by the gel-shift assay. In parallel, the conditions for possible self-hybridization of pcPNA oligomers have been assayed by mixing curve (Job plot) and thermal melting experiments. The data indicate that, at physiological temperatures ( approximately 37 degrees C), the equilibrium is shifted toward the pairing of corresponding pcPNAs with each other. This finding explains a linear concentration dependence, within the submicromolar range, of the pcPNA invasion rate into dsDNA at 37 degrees C. At elevated temperatures (>50 degrees C), the rather unstable pcPNA duplexes dissociate, yielding the expected quadratic dependence for the rate of pcPNA invasion on the PNA concentration. The polycationic character of pcPNA pairs, carrying the duplicated number of protonated terminal PNA residues commonly used to increase the PNA solubility and binding affinity, also explains the self-inhibition of pcPNA invasion observed at higher PNA concentrations. Melting of pcPNA duplexes occurs with the integral transition enthalpies ranged from -235 to -280 kJ.mol(-1), contributing to an anomalously high activation energy of approximately 150 kJ.mol(-1) found for the helix invasion of pcPNAs carrying four different nucleobases. A simplified kinetic model for pcPNAs helix invasion is proposed that interprets all unusual features of pcPNAs binding to dsDNA. Our findings have important implications for rational use of pcPNAs.