Diffusion of isolated DNA molecules: dependence on length and topology

DNA confinement in nanochannels: physics and biological applications

DNA is the central storage molecule of genetic information in the cell, and reading that information is a central problem in biology. While sequencing technology has made enormous advances over the past decade, there is growing interest in platforms that can readout genetic information directly from long single DNA molecules, with the ultimate goal of single-cell, single-genome analysis. Such a capability would obviate the need for ensemble averaging over heterogeneous cellular populations and eliminate uncertainties introduced by cloning and molecular amplification steps (thus enabling direct assessment of the genome in its native state). In this review, we will discuss how the information contained in genomic-length single DNA molecules can be accessed via physical confinement in nanochannels. Due to self-avoidance interactions, DNA molecules will stretch out when confined in nanochannels, creating a linear unscrolling of the genome along the channel for analysis. We will first review the fundamental physics of DNA nanochannel confinement--including the effect of varying ionic strength--and then discuss recent applications of these systems to genomic mapping. Apart from the intense biological interest in extracting linear sequence information from elongated DNA molecules, from a physics view these systems are fascinating as they enable probing of single-molecule conformation in environments with dimensions that intersect key physical length-scales in the 1 nm to 100 µm range.

Biologistics--diffusion coefficients for complete proteome of Escherichia coli

Motivation: Biologistics provides data for quantitative analysis of transport (diffusion) processes and their spatio-temporal correlations in cells. Mobility of proteins is one of the few parameters necessary to describe reaction rates for gene regulation. Although understanding of diffusion-limited biochemical reactions in vivo requires mobility data for the largest possible number of proteins in their native forms, currently, there is no database that would contain the complete information about the diffusion coefficients (DCs) of proteins in a given cell type.

Results: We demonstrate a method for the determination of in vivo DCs for any molecule—regardless of its molecular weight, size and structure—in any type of cell. We exemplify the method with the database of in vivo DC for all proteins (4302 records) from the proteome of K12 strain of Escherichia coli, together with examples of DC of amino acids, sugars, RNA and DNA. The database follows from the scale-dependent viscosity reference curve (sdVRC). Construction of sdVRC for prokaryotic or eukaryotic cell requires ~20 in vivo measurements using techniques such as fluorescence correlation spectroscopy (FCS), fluorescence recovery after photobleaching (FRAP), nuclear magnetic resonance (NMR) or particle tracking. The shape of the sdVRC would be different for each organism, but the mathematical form of the curve remains the same. The presented method has a high predictive power, as the measurements of DCs of several inert, properly chosen probes in a single cell type allows to determine the DCs of thousands of proteins. Additionally, obtained mobility data allow quantitative study of biochemical interactions in vivo.

Contact:rholyst@ichf.edu.pl

Supplementary information:Supplementary data are available at Bioinformatics Online.

Lamination and mixing in three fundamental flow sequences driven by electromagnetic body forces

This article pursues the idea that the degree of striations, called lamination, could be engineered to complement stretching and to design new sequential mixers. It explores lamination and mixing in three new mixing sequences experimentally driven by electromagnetic body forces. To generate these three mixing sequences, Lorentz body forces are dynamically controlled to vary the flow geometry produced by a pair of local jets. The first two sequences are inspired from the "tendril and whorl" and "blinking vortex" flows. The third novel sequence is called the "cat's eyes flip." These three mixing sequences exponentially stretch and laminate material lines representing the interface between two domains to be mixed. Moreover, the mixing coefficient (defined as 1-σ(2)/σ(0)(2) where σ(2)/σ(0)(2) is the rescaled variance) and its rate grow exponentially before saturation. This saturation of the mixing process is related to the departure of the mixing rate from an exponential growth when the striations' thicknesses reach the diffusive length scale of the measurements or species and dyes. Incidentally, in our experiments, for the same energy or forcing input, the cat's eyes flip sequence has higher lamination, stretching, and mixing rates than the tendril and whorl and the blinking vortex sequences. These features show that bakerlike in situ mixers can be conceived by dynamically controlling a pair of local jets and by integrating lamination during stirring stages with persistent geometries. Combined with novel insights provided by the quantification of the lamination, this paper should offer perspectives for the development of new sequential mixers, possibly on all scales.

Mobility of a semiflexible chain confined in a nanochannel

The classic results of de Gennes and Odijk describe the mobility of a semiflexible chain confined in a nanochannel only in the limits of very weak and very strong confinement, respectively. Using Monte Carlo sampling of the Kirkwood diffusivity with full hydrodynamic interactions, we show that the mobility of a semiflexible chain exhibits a broad plateau as a function of extension before transitioning to an Odijk regime, and that the width of the plateau depends on the anisotropy of the monomers. For the particular case of DNA in a high ionic strength buffer, which has highly anisotropic monomers, we predict that this Rouse-like behavior will be observed over most of the measurable chain extensions seen in experiments.

Single-molecule studies of viral DNA packaging

Many double-stranded DNA bacteriophages and viruses use specialized ATP-driven molecular machines to package their genomes into tightly confined procapsid shells. Over the last decade, single-molecule approaches - and in particular, optical tweezers - have made key contributions to our understanding of this remarkable process. In this chapter, we review these advances and the insights they have provided on the packaging mechanisms of three bacteriophages: φ 29, λ, and T4.

Looped star polymers show conformational transition from spherical to flat toroidal shapes

Inspired by the topological organization of the circular Escherichia coli chromosome, which is compacted by separate domains, we study a polymer architecture consisting of a central ring to which either looped or linear side chains are grafted. A shape change from a spherical to a toroidal organization takes place as soon as the inner ring becomes large enough for the attached arms to fit within its circumference. Building up a torus, the system flattens, depending on the effective bending rigidity of the chain induced by entropic repulsion of the attached loops and, to a lesser extent, linear arms. Our results suggest that the natural formation of a toroidal structure with a decreased amount of writhe induced by a specific underlying topology could be one driving force, among others, that nature exploits to ensure proper packaging of the genetic material within a rod-shaped, bacterial envelope.

A microfluidics assisted porous silicon array for optical label-free biochemical sensing

A porous silicon (PSi) based microarray has been integrated with a microfluidic system, as a proof of concept device for the optical monitoring of selective label-free DNA-DNA interaction. A 4 × 4 square matrix of PSi one dimensional photonic crystals, each one of 200 μm diameter and spaced by 600 μm, has been sealed by a polydimethylsiloxane (PDMS) channels circuit. The PSi optical microarray elements have been functionalized by DNA single strands after sealing: the microfluidic circuit allows to reduce significantly the biologicals and chemicals consumption, and also the incubation time with respect to a not integrated device. Theoretical calculations, based on finite element method, taking into account molecular interactions, are in good agreement with the experimental results, and the developed numerical model can be used for device optimization. The functionalization process and the interaction between DNA probe and target has been monitored by spectroscopic reflectometry for each PSi element in the microchannels.

Partial hydrodynamic screening of confined linear and circular double-stranded DNA dynamics

We performed experiments and simulations to investigate the influence of hydrodynamic interaction on the diffusion dynamics of circular and linear λ-DNA confined in nanoslits. Contrary to the common assumption that intrachain hydrodynamic interaction (HI) is completely screened when polymers are confined in channels with height h smaller than the radius of gyration R(g), it is found that the HI is partially screened and approaches complete screening only for R(g)≪h. For λ-DNA, the HI becomes nearly completely screened only when the channel height is smaller than the Kuhn length. In addition, the dynamics of linear and circular λ-DNA in very strong confinement is shown to be independent of the chain topology.

Simulation of DNA Extension in Nanochannels

We have used a realistic model for double stranded DNA and Monte Carlo simulations to compute the extension (mean span) of a DNA molecule confined in a nanochannel over the full range of confinement in a high ionic strength buffer. The simulation data for square nanochannels resolve the apparent contradiction between prior simulation studies and the predictions from Flory theory, demonstrating the existence of two transition regimes between weak confinement (the de Gennes regime) and strong confinement (the Odijk regime). The simulation data for rectangular nanochannels support the use of the geometric mean for mapping data obtained in rectangular channels onto models developed for cylinders. The comparison of our results with experimental data illuminates the challenges in applying models for confined, neutral polymers to polyelectrolytes. Using a Flory-type approach, we also provide an improved scaling result for the relaxation time in the transition regime close to that found in experiments.

Conformation of ring polymers in 2D constrained environments

The combination of ring closure and spatial constraints has a fundamental effect on the statistics of semiflexible polymers such as DNA. However, studies of the interplay between circularity and constraints are scarce and single-molecule experimental data concerning polymer conformations are missing. By means of atomic force microscopy we probe the conformation of circular DNA molecules in two dimensions and in the concentrated regime (above the overlap concentration c*). Molecules in this regime experience a collapse, and their statistical properties agree very well with those of simulated vesicles under pressure. Some circular molecules also create confining regions in which other molecules are trapped. Thus we show further that spatially confined molecules fold into specific conformations close to those found for linear chains, and strongly dependent on the size of the confining box.

RNA dimerization monitored by fluorescence correlation spectroscopy

Fluorescence correlation spectroscopy (FCS) provides a versatile tool to investigate molecular interaction under native conditions, approximating infinite dilution. One precondition for its application is a sufficient difference between the molecular weights of the fluorescence-labelled unbound and bound ligand. In previous studies, an 8-fold difference in molecular weights or correspondingly a 1.6-fold difference in diffusion coefficients was required to accurately distinguish between two diffusion species by FCS. In the presented work, the hybridization of two complementary equally sized RNA single strands was investigated at an excellent signal-to-noise ratio enabled by the highly photostable fluorophore Atto647N. The fractions of ssRNA and dsRNA were quantified by applying multicomponent model analysis of single autocorrelation functions and globally fitting several autocorrelation functions. By introducing a priori knowledge into the fitting procedure, 1.3- to 1.4-fold differences in diffusion coefficients of single- and double-stranded RNA of 26, 41, and 54 nucleotides could be accurately resolved. Global fits of autocorrelation functions of all titration steps enabled a highly accurate quantification of diffusion species fractions and mobilities. At a high signal-to-noise ratio, the median of individually fitted autocorrelation functions allowed a robust representation of heterogeneous data. These findings point out the possibility of studying molecular interaction of equally sized molecules based on their diffusional behavior, which significantly broadens the application spectrum of FCS.

Comparative analysis of viscosity of complex liquids and cytoplasm of mammalian cells at the nanoscale

We present a scaling formula for size-dependent viscosity coefficients for proteins, polymers, and fluorescent dyes diffusing in complex liquids. The formula was used to analyze the mobilities of probes of different sizes in HeLa and Swiss 3T3 mammalian cells. This analysis unveils in the cytoplasm two length scales: (i) the correlation length ξ (approximately 5 nm in HeLa and 7 nm in Swiss 3T3 cells) and (ii) the limiting length scale that marks the crossover between nano- and macroscale viscosity (approximately 86 nm in HeLa and 30 nm in Swiss 3T3 cells). During motion, probes smaller than ξ experienced matrix viscosity: η(matrix) ≈ 2.0 mPa·s for HeLa and 0.88 mPa·s for Swiss 3T3 cells. Probes much larger than the limiting length scale experienced macroscopic viscosity, η(macro) ≈ 4.4 × 10(-2) and 2.4 × 10(-2) Pa·s for HeLa and Swiss 3T3 cells, respectively. Our results are persistent for the lengths scales from 0.14 nm to a few hundred nanometers.

Microfluidic devices for bioapplications

Harnessing the ability to precisely and reproducibly actuate fluids and manipulate bioparticles such as DNA, cells, and molecules at the microscale, microfluidics is a powerful tool that is currently revolutionizing chemical and biological analysis by replicating laboratory bench-top technology on a miniature chip-scale device, thus allowing assays to be carried out at a fraction of the time and cost while affording portability and field-use capability. Emerging from a decade of research and development in microfluidic technology are a wide range of promising laboratory and consumer biotechnological applications from microscale genetic and proteomic analysis kits, cell culture and manipulation platforms, biosensors, and pathogen detection systems to point-of-care diagnostic devices, high-throughput combinatorial drug screening platforms, schemes for targeted drug delivery and advanced therapeutics, and novel biomaterials synthesis for tissue engineering. The developments associated with these technological advances along with their respective applications to date are reviewed from a broad perspective and possible future directions that could arise from the current state of the art are discussed.

The distribution of DNA translocation times in solid-state nanopores

This paper systematically investigates the effects of solution viscosity, applied voltage and DNA chain length on the distribution of DNA translocation times through 8 ± 2 nm diameter silicon nitride nanopores. Linear dsDNA translocation events were selected based on the magnitude of current blockage and accumulated into scatter plots of current blockage and event duration (translocation time). The translocation time distribution was fitted to the solution of a Smoluchowski-type equation for 1D biased diffusion to a sink. The DNA drifting speed under bias and diffusion constant were extracted from the fits as functions of solution viscosity, applied voltage and DNA chain length. Combined with the Einstein-Smoluchowski relation, this model allowed evaluation of the viscous drag force on DNA molecules. This model also allowed estimation of the uncertainty in determining the DNA chain length due to the influence of friction on the spread of translocation times in a nanopore measurement. The data analysis suggests that the simple 1D biased diffusion model fits the experimental data well for a wide range of conditions. Some deviations from predicted behavior were observed and show where additional phenomena are likely to contribute to the distribution of DNA translocation times.

Predicting translational diffusion of evolutionary conserved RNA structures by the nucleotide number

Ribonucleic acids are highly conserved essential parts of cellular life. RNA function is determined to a large extent by its hydrodynamic behaviour. The presented study proposes a strategy to predict the hydrodynamic behaviour of RNA single strands on the basis of the polymer size. By atom-level shell-modelling of high-resolution structures, hydrodynamic radius and diffusion coefficient of evolutionary conserved RNA single strands (ssRNA) were calculated. The diffusion coefficients D of 17–174 nucleotides (nt) containing ssRNA depended on the number of nucleotides N with D = 4.56 × 10⁻¹⁰ N^−0.39 m² s⁻¹. The hydrodynamic radius R_H depended on N with R_H = 5.00 × 10⁻¹⁰N^0.38 m. An average ratio of the radius of gyration and the hydrodynamic radius of 0.98 ± 0.08 was calculated in solution. The empirical law was tested by in solution measured hydrodynamic radii and radii of gyration and was found to be highly consistent with experimental data of evolutionary conserved ssRNA. Furthermore, the hydrodynamic behaviour of several evolutionary unevolved ribonucleic acids could be predicted. Based on atom-level shell-modelling of high-resolution structures and experimental hydrodynamic data, empirical models are proposed, which enable to predict the translational diffusion coefficient and molecular size of short RNA single strands solely on the basis of the polymer size.

Biophysics of knotting

Knots appear in a wide variety of biophysical systems, ranging from biopolymers, such as DNA and proteins, to macroscopic objects, such as umbilical cords and catheters. Although significant advancements have been made in the mathematical theory of knots and some progress has been made in the statistical mechanics of knots in idealized chains, the mechanisms and dynamics of knotting in biophysical systems remain far from fully understood. We report on recent progress in the biophysics of knotting-the formation, characterization, and dynamics of knots in various biophysical contexts.

Precise characterization of the conformation fluctuations of freely diffusing DNA: beyond Rouse and Zimm

We studied the dynamics of single freely diffusing fluorescence-labeled double-stranded lambda-phage DNA molecules using dual-color 3-dimensional feedback tracking microscopy and intramolecular fluorescence correlation spectroscopy. Our technique is independently sensitive to the molecule's diffusion coefficient D and radius of gyration R(g) and is concentration insensitive, providing greater precision for characterizing the molecule's intramolecular motion than other methods. We measured D = 0.80 +/- 0.05 microm(2)/s and R(g) approximately 420 nm, consistent with the Kirkwood-Riseman prediction for a flexible polymer with strong hydrodynamic interactions (HI), but we find the statistics of intramolecular motion inconsistent with the Zimm model for such a polymer. We address a dispute in the experimental literature, finding that previous measurements on double-stranded DNA likely lacked the sensitivity to distinguish between the Zimm model and the HI-free Rouse model. Finally, we observe fluorescence fluctuations with a correlation time of over 2 s that cannot be explained by either model and propose that they may be signatures of excluded volume interactions.

Internal dynamics of supercoiled DNA molecules

The intramolecular diffusive motion within supercoiled DNA molecules is of central importance for a wide array of gene regulation processes. It has recently been shown, using fluorescence correlation spectroscopy, that plasmid DNA exhibits unexpected acceleration of its internal diffusive motion upon supercoiling to intermediate density. Here, we present an independent study that shows a similar acceleration for fully supercoiled plasmid DNA. We have developed a method that allows fluorescent labeling of a 200-bp region, as well as efficient supercoiling by Escherichia coli gyrase. Compared to plain circular or linear DNA, the submicrosecond motion within the supercoiled molecules appears faster by up to an order of magnitude. The mean-square displacement as a function of time reveals an additional intermediate regime with a lowered scaling exponent compared to that of circular DNA. Although this unexpected behavior is not fully understood, it could be explained by conformational constraints of the DNA strand within the supercoiled topology in combination with an increased apparent persistence length.

Steric effects and mass-transfer limitations surrounding amplification reactions on immobilized long and clinically relevant DNA templates

DNA and RNA are commonly captured on solid substrates during purification and isolation, where they can be transferred to downstream amplification and transcription reactions. When compared to the solution phase, however, immobilized DNA- and RNA-directed reactions are less efficient because of a variety of complex factors. Steric inhibition because of the bead surface and neighboring biological polymers, a change in solution chemistry because of the high local concentration of template molecules, and mass transfer to the bead surface could all affect the overall reaction kinetics. Furthermore, these effects may be particularly evident when working with long clinically relevant molecules, such as mRNA, viral RNA, and cDNA. In this paper, we focus on the in vitro transcription reaction (IVT) of both a long and short strand of H5 influenza A RNA (1777 and 465 nt) on both free and immobilized DNA templates to study these phenomena. We found that transcription was less efficient on immobilized beads than in solution, but that it can be dramatically increased with optimal solution chemistry. Using high ribonucleotide concentrations (>6 mM total rNTP), the RNA yield from long immobilized cDNA templates was boosted to 60% of solution control. Surprisingly, we found that steric effects because of surrounding immobilized molecules were only significant when the DNA molecules were short enough to achieve a high density (9x10(-4) microm2/molecule) on the silica substrate, such that the gap between molecules is on the order of the polymerase diameter. Eventually, these findings can be exploited in an automated microreactor, where isolation, purification, amplification, and detection of nucleic acids can be unified into one portable device.

Universality in the diffusion of knots

We have evaluated a universal ratio between diffusion constants of the ring polymer with a given knot K and a linear polymer with the same molecular weight in solution through the Brownian dynamics under hydrodynamic interaction. The ratio is found to be constant with respect to the number of monomers, N , and hence the estimate at some N should be valid practically over a wide range of N for various polymer models. Interestingly, the ratio is determined by the average crossing number (N{AC}) of an ideal conformation of knotted curve K , i.e., that of the ideal knot. The N{AC} of ideal knots should therefore be fundamental in the dynamics of knots.

Atomic hydrodynamics of DNA: coil-uncoil-coil transitions in a wall-bounded shear flow

Extensive experimental work on the response of DNA molecules to externally applied forces and on the dynamics of DNA molecules flowing in microchannels and nanochannels has been carried out over the past two decades, however, there has not been available, until now, any atomic-scale means of analyzing nonequilibrium DNA response dynamics. There has not therefore been any way to investigate how the backbone and side-chain atoms along the length of a DNA molecule interact with the molecules and ions of the flowing solvent and with the atoms of passing boundary surfaces. We report here on the application of the nonequilibrium biomolecular dynamics simulation method that we developed [G. M. Wang and W. C. Sandberg, Nanotechnology 18, 4819 (2007)] to analyze, at the atomic interaction force level, the conformational dynamics of short-chain single-stranded DNA molecules in a shear flow near a surface. This is a direct atomic computational analysis of the hydrodynamic interaction between a biomolecule and a flowing solvent. The DNA molecules are observed to exhibit conformational behaviors including coils, hairpin loops, and figure-eight shapes that have neither been previously measured experimentally nor observed computationally, as far as we know. We relate the conformational dynamics to the atomic interaction forces experienced throughout the length of a molecule as it moves in the flowing solvent past the surface boundary. We show that the DNA conformational dynamics is related to the asymmetry in the molecular environment induced by the motion of the surrounding molecules and the atoms of the passing surface. We also show that while the asymmetry in the environment is necessary, it is not sufficient to produce the observed conformational dynamics. A time variation in the asymmetry, due in our case to a shear flow, must also exist. In order to contrast these results with the usual experimental situation of purely diffusive motion in thermal equilibrium we have also carried out computations with a zero shear rate. We show that in thermal equilibrium there is asymmetry and an atomic hydrodynamic coupling between DNA molecules and the solvent molecules but there is no coil-uncoil transition.

Conformation of circular DNA in two dimensions

The conformation of circular DNA molecules of various lengths adsorbed in a 2D conformation on a mica surface is studied. The results confirm the conjecture that the critical exponent nu is topologically invariant and equal to the self-avoiding walk value (in the present case nu=3/4), and that the topology and dimensionality of the system strongly influence the crossover between the rigid regime and the self-avoiding regime at a scale L approximately 7l{p}. Additionally, the bond correlation function scales with the molecular length L as predicted. For molecular lengths L<or=5l{p}, circular DNA behaves like a stiff molecule with an approximately elliptic shape.

Effects of size and topology of DNA molecules on intracellular delivery with non-viral gene carriers

Background

Efforts to improve the efficiency of non-viral gene delivery require a better understanding of delivery kinetics of DNA molecules into clinically relevant cells. Towards this goal, three DNA molecules were employed to investigate the effects of DNA properties on cellular delivery: a circular plasmid DNA (c-DNA), a linearized plasmid DNA (l-DNA) formulated by single-site digestion of c-DNA, and smaller linear gene cassette generated by PCR (pcr-DNA). Four non-viral gene carriers were investigated for DNA delivery: polyethyleneimine (PEI), poly-L-Lysine (PLL), palmitic acid-grafted PLL (PLL-PA), and Lipofectamine-2000™. Particle formation, binding and dissociation characteristics, and DNA uptake by rat bone marrow stromal cells were investigated.

Results

For individual carriers, there was no discernible difference in the morphology of particles formed as a result of carrier complexation with different DNA molecules. With PEI and PLL carriers, no difference was observed in the binding interaction, dissociation characteristics, and DNA uptake among the three DNA molecules. The presence of serum in cell culture media did not significantly affect the DNA delivery by the polymeric carriers, unlike other lipophilic carriers. Using PEI as the carrier, c-DNA was more effective for transgene expression as compared to its linear equivalent (l-DNA) by using the reporter gene for Enhanced Green Fluorescent Protein. pcr-DNA was the least effective despite being delivered into the cells to the same extent.

Conclusion

We conclude that the nature of gene carriers was the primary determinant of cellular delivery of DNA molecules, and circular form of the DNA was more effectively processed for transgene expression.

Prediction of diffusion coefficients of plasmids

The ability to predict diffusion coefficients is important in the design, analysis, and operation of plasmid downstream processing operations such as membrane and fixed-bed chromatography. A correlation is proposed to predict the diffusion coefficient, D, of supercoiled plasmid DNA molecules in dilute solutions on the basis of the molecular weight, M, or size. Experimental data (18 points) collected from the literature confirmed the proposed variation of D with plasmid molecular weight as D proportional, variant M(-2/3), for molecules within the 1,800-287,100 base-pair range. The correlation was able to estimate the available experimental results with an average error of 6.3%.

Direct measurement of the intermolecular forces confining a single molecule in an entangled polymer solution

We use optical tweezers to directly measure the intermolecular forces acting on a single polymer imposed by surrounding entangled polymers (115 kbp DNA, 1 mg/ml). A tubelike confining field was measured in accord with the key assumption of reptation models. A time-dependent harmonic potential opposed transverse displacement, in accord with recent simulation findings. A tube radius of 0.8 microm was determined, close to the predicted value (0.5 microm). Three relaxation modes (approximately 0.4, 5, and 34 s) were measured following transverse displacement, consistent with predicted relaxation mechanisms.

Static conformation and dynamics of single DNA molecules confined in nanoslits

Nearly thirty years ago, Daoud and de Gennes derived the scaling predictions for the linear polymer chains trapped in a slit with dimension close to the Kuhn length; however, these predictions have yet to be compared with experiments. We have fabricated nanoslits with vertical dimension similar to the Kuhn length of ds-DNA (110nm) using standard photolithography techniques. Fluorescently labeled single DNA molecules with contour lengths L ranging from 4 to 75 microm were successfully injected into the slits and the chain molecules undergoing Brownian motions were imaged by fluorescence microscopy. The distributions of the chain radius of gyration and the two-dimensional asphericity were measured. It is found that the DNA molecules exhibit highly anisotropic shape and the mean asphericity is chain length independence. The shape anisotropy of DNA in our measurements is between two and three dimensions (2D and 3D). The static scaling law of the chain extension and the radius of gyration <R parallel>, <Rg> approximately L(nu) were observed with nuR(parallel)=0.65+/-0.02 and nu(Rg)=0.68+/-0.05. These results are close to the average value between two (nuR parallel,Rg=0.75) and three (nuR parallel,Rg=0.6) -dimensional theoretical value. The scaling of the extensional and rotational relaxation time are between Rouse model in nanoslits and Zimm model in the bulk solution, respectively. We show that the conformation and chain relaxation of DNA confined in a slit close to Kuhn length exhibit the quasi-two-dimensional behavior.

Hybrid method coupling fluctuating hydrodynamics and molecular dynamics for the simulation of macromolecules

We present a hybrid computational method for simulating the dynamics of macromolecules in solution which couples a mesoscale solver for the fluctuating hydrodynamics (FH) equations with molecular dynamics to describe the macromolecule. The two models interact through a dissipative Stokesian term first introduced by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We show that our method correctly captures the static and dynamical properties of polymer chains as predicted by the Zimm model. In particular, we show that the static conformations are best described when the ratio sigma/b=0.6, where sigma is the Lennard-Jones length parameter and b is the monomer bond length. We also find that the decay of the Rouse modes' autocorrelation function is better described with an analytical correction suggested by Ahlrichs and Dunweg. Our FH solver permits us to treat the fluid equation of state and transport parameters as direct simulation parameters. The expected independence of the chain dynamics on various choices of fluid equation of state and bulk viscosity is recovered, while excellent agreement is found for the temperature and shear viscosity dependence of center of mass diffusion between simulation results and predictions of the Zimm model. We find that Zimm model approximations start to fail when the Schmidt number Sc < or approximately 30. Finally, we investigate the importance of fluid fluctuations and show that using the preaveraged approximation for the hydrodynamic tensor leads to around 3% error in the diffusion coefficient for a polymer chain when the fluid discretization size is greater than 50 A.

Plasmid purification using non-porous anion-exchange silica fibres

A new type of fibre-based anion-exchange material for plasmid purification was developed. The basic material consisted of non-porous silica fibres with a mean diameter of 1.5 microm and a surface area of 2.4m(2)g(-1). The fibre surface was provided with several types of ligands, either by adsorption of polymers (chitosan or poly(ethyleneimine)) or by polymerization of amine-containing acrylic monomers onto a propyl methacrylate-silanized surface. The resulting polymer layers contained primary, tertiary or quaternary amines as ion-exchange groups. The packing density could be varied considerably, 9-34% (v/v). The loose packing structure provided excellent flow properties suitable for high-speed operations. The best overall performance was shown by silica fibres provided with tertiary amine polymers, having a plasmid-binding capacity of 0.9 mg ml(-1) (pre-purified plasmid) and a plasmid recovery of 62% (performance data remained stable though several adsorption cycles). The high flow rates possible with the fibre material made it especially useful when large volumes of cleared lysate were processed. The columns could be operated with retention of their adsorption properties at speeds of up to 1800 cm h(-1), equivalent to 0.5 column volumes per minute. The binding capacity was found to be lower than anticipated from the design of the fibres. Fluorescence imaging showing individual plasmid molecules indicated the fibre population to be heterogeneous with respect to plasmid adsorption, some fibres displaying poor binding properties. Possible reasons for this heterogeneity are discussed.

Internal mechanical response of a polymer in solution

We observed single molecules of fluorescently labeled double-stranded (ds) lambda DNA held in an anti-Brownian electrokinetic trap. From the measured density fluctuations we extract the density-density response function of the molecule over times >4.5 ms and distances >250 nm, i.e., how a perturbation in density in one part of the molecule propagates through the rest of the molecule. We find a nonmonotonic radial dependence of the relaxation time. In contrast with earlier measurements on freely diffusing dsDNA, we observe clear signs of internal hydrodynamic interactions.

Strong effects of molecular topology on diffusion of entangled DNA molecules

When long polymers such as DNA are in a highly concentrated state they may become entangled, leading to restricted self-diffusion. Here, we investigate the effect of molecular topology on diffusion in concentrated DNA solutions and find surprisingly large effects, even with molecules of modest length and concentration. We measured the diffusion coefficients of linear and relaxed circular molecules by tracking the Brownian motion of single molecules with fluorescence microscopy. Four possible cases were compared: linear molecules surrounded by linear molecules, circular molecules surrounded by linear molecules, linear molecules surrounded by circles, and circles surrounded by circles. In measurements with 45-kbp DNA at 1 mg/ml, we found that circles diffused approximately 100 times slower when surrounded by linear molecules than when surrounded by circles. In contrast, linear and circular molecules diffused at nearly the same rate when surrounded by circles, and circles diffused approximately 10 times slower than linears when surrounded by linears. Thus, diffusion in entangled DNA solutions strongly depends on topology of both the diffusing molecule and the surrounding molecules. This effect also strongly depends on DNA concentration and length. The differences largely disappeared when the concentration was lowered to 0.1 mg/ml or when the DNA length was lowered to 6 kb. Present theories cannot fully explain these effects.

Detection and identification of IHN and ISA viruses by isothermal DNA amplification in microcapillary tubes

Unique base sequences derived from RNA of both infectious hematopoietic necrosis virus (IHNV) and infectious salmon anemia virus (ISAV) were detected and identified using a combination of surface-associated molecular padlock DNA probes (MPPs) and rolling circle amplification (RCA) in microcapillary tubes. DNA oligonucleotides with base sequences identical to RNA obtained from IHNV or ISAV were recognized by MPPs. Circularized MPPs were then captured on the inner surfaces of glass microcapillary tubes by immobilized DNA oligonucleotide primers. Extension of the immobilized primers by isothermal RCA produced DNA concatamers, which were labeled with fluorescent SYBR Green II nucleic acid stain, and measured by microfluorimetry. Molecular padlock probes, combined with this method of surface-associated isothermal RCA, exhibited high selectivity without the need for thermal cycling. This method is applicable to the design of low-power field sensors capable of multiplex detection of viral, bacterial, and protozoan pathogens within localized regions of microcapillary tubes.