Simulation of conformational preconditioning strategies for electrophoretic stretching of DNA in a microcontraction

We have used Brownian dynamics-finite element method to examine two conformational preconditioning approaches for improving DNA stretching in a microcontraction for the purpose of direct gene analysis. The newly proposed "pre-stretching" strategy is found to significantly improve the degree of DNA extension at the exit of the contraction. On the other hand, applying an oscillating extensional field to DNA yields no preconditioning effect. Detailed analysis of the evolution of DNA extension and conformation reveals that the success of our "pre-stretching" strategy relies on the "non-local" effect that cannot be predicted using simple kinematics analysis. In other words, accurate prediction can only be obtained using detailed simulations. Comparing to the existing preconditioning strategies, our "pre-stretching" method is easy to implement while still providing a very good performance. We hope that the insight gained from this study can be useful for future design of biomicrofluidic devices for DNA manipulation.     

Simulation guided design of a microfluidic device for electrophoretic stretching of DNA

We have used Brownian dynamics-finite element method (BD-FEM) to guide the optimization of a microfluidic device designed to stretch DNA for gene mapping. The original design was proposed in our previous study [C. C. Hsieh and T. H. Lin, Biomicrofluidics 5(4), 044106 (2011)] for demonstrating a new pre-conditioning strategy to facilitate DNA stretching through a microcontraction using electrophoresis. In this study, we examine the efficiency of the original device for stretching DNA with different sizes ranging from 48.5 kbp (λ-DNA) to 166 kbp (T4-DNA). The efficiency of the device is found to deteriorate with increasing DNA molecular weight. The cause of the efficiency loss is determined by BD-FEM, and a modified design is proposed by drawing an analogy between an electric field and a potential flow. The modified device does not only regain the efficiency for stretching large DNA but also outperforms the original device for stretching small DNA.     

Stretching DNA by electric field and flow field in microfluidic devices: An experimental validation to the devices designed with computer simulations

We examined the performance of three microfluidic devices for stretching DNA. The first device is a microchannel with a contraction, and the remaining two are the modifications to the first. The modified designs were made with the help of computer simulations [C. C. Hsieh and T. H. Lin, Biomicrofluidics 5(4), 044106 (2011) and C. C. Hsieh, T. H. Lin, and C. D. Huang, Biomicrofluidics 6, 044105 (2012)] and they were optimized for operating with electric field. In our experiments, we first used DC electric field to stretch DNA. However, the experimental results were not even in qualitative agreement with our simulations. More detailed investigation revealed that DNA molecules adopt a globular conformation in high DC field and therefore become more difficult to stretch. Owing to the similarity between flow field and electric field, we turned to use flow field to stretch DNA with the same devices. The evolution patterns of DNA conformation in flow field were found qualitatively the same as our prediction based on electric field. We analyzed the maximum values, the evolution and the distributions of DNA extension at different Deborah number in each device. We found that the shear and the hydrodynamic interaction have significant influence on the performance of the devices.     

Tilted post arrays for separating long DNA

Recent simulations by Chen and Dorfman [Electrophoresis 35, 405–411 (2014)] suggested that “tilting” the electric field with respect to the lattice vectors of a hexagonal post array would lead to a substantial improvement in electrophoretic DNA separations therein. We constructed such an array where the electric field is applied at an angle equidistant between the two lattice vectors. This tilted array leads to (i) baseline resolution of 20 kbp DNA and λ DNA (48.5 kbp) in a 4 mm channel and (ii) measurable separation resolutions for electric fields up to 50 V/cm, both of which are improvements over untilted post arrays of the same post density. The predicted time required to reach a resolution of unity is approximately 5 min, independent of electric field. The separations are more reproducible at higher fields.     

Onset of channeling during DNA electrophoresis in a sparse ordered post array

The “channeling hypothesis” of DNA electrophoresis in sparse, ordered arrays of posts predicts that the DNA will move through the array relatively unhindered if (i) the spacing between the posts is larger than the DNA coil and (ii) the electric field lines are straight. We tested this hypothesis by studying the electrophoretic separation of a small plasmid DNA (pUC19, 2686 base pairs) and a large, linear DNA (λ-DNA, 48 500 base pairs) in a hexagonal array of 1 μm diameter posts with a pitch of 7 μm. At low electric field strengths, these DNAs are separated due to the long-lived, rope-over-pulley collisions of λ-DNA with the posts. The resolution is lost as the electric field increases due to the onset of channeling by the λ-DNA. Using a diffusive model, we show that channeling arises at low electric fields due to the finite size of the array. This channeling is not intrinsic to the system and is attenuated by increasing the size of the array. Higher electric fields lead to intrinsic channeling, which is attributed to the disparate time scales for a rope-over-pulley collision and transverse diffusion between collisions. The onset of channeling is a gradual process, in agreement with extant Brownian dynamics simulation data. Even at weak electric fields, the electrophoretic mobility of λ-DNA in the array is considerably higher than would be expected if the DNA frequently collided with the posts.     

Fluctuations of DNA mobility in nanofluidic entropic traps

We studied the mobility of DNA molecules driven by an electric field through a nanofluidic device containing a periodic array of deep and shallow regions termed entropic traps. The mobility of a group of DNA molecules was measured by fluorescent video microscopy. Since the depth of a shallow region is smaller than the DNA equilibrium size, DNA molecules are trapped for a characteristic time and must compress themselves to traverse the boundary between deep and shallow regions. Consistent with previous experimental results, we observed a nonlinear relationship between the mobility and electric field strength, and that longer DNA molecules have larger mobility. In repeated measurements under seemingly identical conditions, we measured fluctuations in the mobility significantly larger than expected from statistical variation. The variation was more pronounced for lower electric field strengths where the trapping time is considerable relative to the drift time. To determine the origin of these fluctuations, we investigated the dependence of the mobility on several variables: DNA concentration, ionic strength of the solvent, fluorescent dye staining ratio, electroosmotic flow, and electric field strength. The mobility fluctuations were moderately enhanced in conditions of reduced ionic strength and electroosmotic flow.     

Simulation of single DNA molecule stretching and immobilization in a de-wetting two-phase flow over micropillar-patterned surface

We investigate single DNA stretching dynamics in a de-wetting flow over micropillars using Brownian dynamics simulation. The Brownian dynamics simulation is coupled with transient flow field computation through a numerical particle tracking algorithm. The droplet formation on the top of the micropillar during the de-wetting process creates a flow pattern that allows DNA to stretch across the micropillars. It is found that DNA nanowire forms if DNA molecules could extend across the stagnation point inside the connecting water filament before its breakup. It also shows that DNA locates closer to the top wall of the micropillar has higher chance to enter the flow pattern of droplet formation and thus has higher chance to be stretched across the micropillars. Our simulation tool has the potential to become a design tool for DNA manipulation in complex biomicrofluidic devices.     

Inertia and scaling in deterministic lateral displacement

The ability to separate and analyze chemical species with high resolution, sensitivity, and throughput is central to the development of microfluidics systems. Deterministic lateral displacement (DLD) is a continuous separation method based on the transport of species through an array of obstacles. In the case of force-driven DLD (f-DLD), size-based separation can be modelled effectively using a simple particle-obstacle collision model. We use a macroscopic model to study f-DLD and demonstrate, via a simple scaling, that the method is indeed predominantly a size-based phenomenon at low Reynolds numbers. More importantly, we demonstrate that inertia effects provide the additional capability to separate same size particles but of different densities and could enhance separation at high throughput conditions. We also show that a direct conversion of macroscopic results to microfluidic settings is possible with a simple scaling based on the size of the obstacles that results in a universal curve.     

Magnetophoretic-based microfluidic device for DNA isolation

This paper presents a continuous flow microfluidic device for the separation of DNA from blood using magnetophoresis for biological applications and analysis. This microfluidic bio-separation device has several benefits, including decreased sample handling, smaller sample and reagent volumes, faster isolation time, and decreased cost to perform DNA isolation. One of the key features of this device is the use of short-range magnetic field gradients, generated by a micro-patterned nickel array on the bottom surface of the separation channel. In addition, the device utilizes an array of oppositely oriented, external permanent magnets to produce strong long-range field gradients at the interfaces between magnets, further increasing the effectiveness of the device. A comprehensive simulation is performed using COMSOL Multiphysics to study the effect of various parameters on the magnetic flux within the separation channel. Additionally, a microfluidic device is designed, fabricated, and tested to isolate DNA from blood. The results show that the device has the capability of separating DNA from a blood sample with a purity of 1.8 or higher, a yield of up to 33 μg of polymerase chain reaction ready DNA per milliliter of blood, and a volumetric throughput of up to 50 ml/h.     

全球价值链视角下福建省产业集群的战略升级策略

经过十几年的发展,福建省产业集群初具规模并成为福建省区域经济发展的重要支柱,但从全球价值链视角看也面临着升级阻碍。以全球价值链视角分析了福建省产业集群战略升级中面临的问题,并基于链网互动探讨了福建省产业集群的战略升级策略。     

DNA translocation through short nanofluidic channels under asymmetric pulsed electric field

Investigation of single molecule DNA dynamics in confined environments has led to important applications in DNA analysis, separation, and sequencing. Here, we studied the electrophoretic transport of DNA molecules through nanochannels shorter than the DNA contour length and calculated the associated translocation time curves. We found that the longer T4 DNA molecules required a longer time to traverse a fixed length nanochannel than shorter λ DNA molecules and that the translocation time decreased with increasing electric field which agreed with theoretical predictions. We applied this knowledge to design an asymmetric electric pulse and demonstrate the different responses of λ and T4 DNA to the pulses. We used Brownian dynamics simulations to corroborate our experimental results on DNA translocation behaviour. This work contributes to the fundamental understanding of polymer transport through nanochannels and may help in designing better separation techniques in the future.     

Stretching of DNA confined in nanochannels with charged walls

There is currently a growing interest in control of stretching of DNA inside nanoconfined regions due to the possibility to analyze and manipulate single biomolecules for applications such as DNA mapping and barcoding, which are based on stretching the DNA in a linear fashion. In the present work, we couple Finite Element Methods and Monte Carlo simulations in order to study the conformation of DNA molecules confined in nanofluidic channels with neutral and charged walls. We find that the electrostatic forces become more and more important when lowering the ionic strength of the solution. The influence of the nanochannel cross section geometry is also studied by evaluating the DNA elongation in square, rectangular, and triangular channels. We demonstrate that coupling electrostatically interacting walls with a triangular geometry is an efficient way to stretch DNA molecules at the scale of hundreds of nanometers. The paper reports experimental observations of λ-DNA molecules in poly(dimethylsiloxane) nanochannels filled with solutions of different ionic strength. The results are in good agreement with the theoretical predictions, confirming the crucial role of the electrostatic repulsion of the constraining walls on the molecule stretching.     

企业成长与创新视角下的产业链升级研究

产业链升级是中国企业摆脱跨国公司低端锁定的重要途径。以企业成长与企业创新理论为基础,讨论了企业在不同成长阶段入嵌产业链的程度。研究发现青春期企业创新投入越强入嵌程度越低,吸收能力越强入嵌程度越高。新创期入嵌产业链程度较低,而青春期较高。创新投入与入嵌程度之间关系受到新创期负向调节。同时,吸收能力与入嵌程度之间关系受到新创期的负向调节、青春期的正向调节。研究结论对中国企业产业链升级具有理论贡献与实践启示。     

Evaluation of the Kirkwood approximation for the diffusivity of channel-confined DNA chains in the de Gennes regime

We use Brownian dynamics with hydrodynamic interactions to calculate both the Kirkwood (short-time) diffusivity and the long-time diffusivity of DNA chains from free solution down to channel confinement in the de Gennes regime. The Kirkwood diffusivity in confinement is always higher than the diffusivity obtained from the mean-squared displacement of the center-of-mass, as is the case in free solution. Moreover, the divergence of the local diffusion tensor, which is non-zero in confinement, makes a negligible contribution to the latter diffusivity in confinement. The maximum error in the Kirkwood approximation in our simulations is about 2% for experimentally relevant simulation times. The error decreases with increasing confinement, consistent with arguments from blob theory and the molecular-weight dependence of the error in free solution. In light of the typical experimental errors in measuring the properties of channel-confined DNA, our results suggest that the Kirkwood approximation is sufficiently accurate to model experimental data.     

DNA combing on low-pressure oxygen plasma modified polysilsesquioxane substrates for single-molecule studies

Molecular combing and flow-induced stretching are the most commonly used methods to immobilize and stretch DNA molecules. While both approaches require functionalization steps for the substrate surface and the molecules, conventionally the former does not take advantage of, as the latter, the versatility of microfluidics regarding robustness, buffer exchange capability, and molecule manipulation using external forces for single molecule studies. Here, we demonstrate a simple one-step combing process involving only low-pressure oxygen (O₂) plasma modified polysilsesquioxane (PSQ) polymer layer to facilitate both room temperature microfluidic device bonding and immobilization of stretched single DNA molecules without molecular functionalization step. Atomic force microscopy and Kelvin probe force microscopy experiments revealed a significant increase in surface roughness and surface potential on low-pressure O₂ plasma treated PSQ, in contrast to that with high-pressure O₂ plasma treatment, which are proposed to be responsible for enabling effective DNA immobilization. We further demonstrate the use of our platform to observe DNA-RNA polymerase complexes and cancer drug cisplatin induced DNA condensation using wide-field fluorescence imaging.     

Structural and dynamic properties of linker histone H1 binding to DNA

Found in all eukaryotic cells, linker histones H1 are known to bind to and rearrange nucleosomal linker DNA. In vitro, the fundamental nature of H1∕DNA interactions has attracted wide interest among research communities-from biologists to physicists. Hence, H1∕DNA binding processes and structural and dynamical information about these self-assemblies are of broad importance. Targeting a quantitative understanding of H1 induced DNA compaction mechanisms, our strategy is based on using small-angle x-ray microdiffraction in combination with microfluidics. The usage of microfluidic hydrodynamic focusing devices facilitates a microscale control of these self-assembly processes, which cannot be achieved using conventional bulk setups. In addition, the method enables time-resolved access to structure formation in situ, in particular, to transient intermediate states. The observed time dependent structure evolution shows that the H1∕DNA interaction can be described as a two-step process: an initial unspecific binding of H1 to DNA is followed by a rearrangement of molecules within the formed assemblies. The second step is most likely induced by interactions between the DNA and the H1's charged side chains. This leads to an increase in lattice spacing within the DNA∕protein assembly and induces a decrease in the correlation length of the mesophases, probably due to a local bending of the DNA.     

Finite element simulations of hydrodynamic trapping in microfluidic particle-trap array systems

Computational fluid dynamic (CFD) simulation is a powerful tool in the design and implementation of microfluidic systems, especially for systems that involve hydrodynamic behavior of objects such as functionalized microspheres, biological cells, or biopolymers in complex structures. In this work, we investigate hydrodynamic trapping of microspheres in a novel microfluidic particle-trap array device by finite element simulations. The accuracy of the time-dependent simulation of a microsphere''s motion towards the traps is validated by our experimental results. Based on the simulation, we study the fluid velocity field, pressure field, and force and stress on the microsphere in the device. We further explore the trap array''s geometric parameters and critical fluid velocity, which affect the microsphere''s hydrodynamic trapping. The information is valuable for designing microfluidic devices and guiding experimental operation. Besides, we provide guidelines on the simulation set-up and release an openly available implementation of our simulation in one of the popular FEM softwares, COMSOL Multiphysics. Researchers may tailor the model to simulate similar microfluidic systems that may accommodate a variety of structured particles. Therefore, the simulation will be of particular interest to biomedical research involving cell or bead transport and migration, blood flow within microvessels, and drug delivery.     

An automated microfluidic system for single-stranded DNA preparation and magnetic bead-based microarray analysis

We present an integrated microfluidic device capable of performing single-stranded DNA (ssDNA) preparation and magnetic bead-based microarray analysis with a white-light detection for detecting mutations that account for hereditary hearing loss. The entire operation process, which includes loading of streptavidin-coated magnetic beads (MBs) and biotin-labeled polymerase chain reaction products, active dispersion of the MBs with DNA for binding, alkaline denaturation of DNA, dynamic hybridization of the bead-labeled ssDNA to a tag array, and white-light detection, can all be automatically accomplished in a single chamber of the microchip, which was operated on a self-contained instrument with all the necessary components for thermal control, fluidic control, and detection. Two novel mixing valves with embedded polydimethylsiloxane membranes, which can alternately generate a 3-μl pulse flow at a peak rate of around 160 mm/s, were integrated into the chip for thoroughly dispersing magnetic beads in 2 min. The binding efficiency of biotinylated oligonucleotides to beads was measured to be 80.6% of that obtained in a tube with the conventional method. To critically test the performance of this automated microsystem, we employed a commercial microarray-based detection kit for detecting nine mutation loci that account for hereditary hearing loss. The limit of detection of the microsystem was determined as 2.5 ng of input K562 standard genomic DNA using this kit. In addition, four blood samples obtained from persons with mutations were all correctly typed by our system in less than 45 min per run. The fully automated, “amplicon-in-answer-out” operation, together with the white-light detection, makes our system an excellent platform for low-cost, rapid genotyping in clinical diagnosis.     

Size-dependent trajectories of DNA macromolecules due to insulative dielectrophoresis in submicrometer-deep fluidic channels

In this paper, we demonstrate for the first time that insulative dielectrophoresis can induce size-dependent trajectories of DNA macromolecules. We experimentally use λ (48.5 kbp) and T4GT7 (165.6 kbp) DNA molecules flowing continuously around a sharp corner inside fluidic channels with a depth of 0.4 μm. Numerical simulation of the electrokinetic force distribution inside the channels is in qualitative agreement with our experimentally observed trajectories. We discuss a possible physical mechanism for the DNA polarization and dielectrophoresis inside confining channels, based on the observed dielectrophoresis responses due to different DNA sizes and various electric fields applied between the inlet and the outlet. The proposed physical mechanism indicates that further extensive investigations, both theoretically and experimentally, would be very useful to better elucidate the forces involved at DNA dielectrophoresis. When applied for size-based sorting of DNA molecules, our sorting method offers two major advantages compared to earlier attempts with insulative dielectrophoresis: Its continuous operation allows for high-throughput analysis, and it only requires electric field strengths as low as ∼10 V∕cm.     

Hydrodynamic particle focusing design using fluid-particle interaction

For passive sheathless particles focusing in microfluidics, the equilibrium positions of particles are typically controlled by micro channels with a V-shaped obstacle array (VOA). The design of the obstacles is mainly based on the distribution of flow streamlines without considering the existence of particles. We report an experimentally verified particle trajectory simulation using the arbitrary Lagrangian-Eulerian (ALE) fluid-particle interaction method. The particle trajectory which is strongly influenced by the interaction between the particle and channel wall is systematically analyzed. The numerical experiments show that the streamline is a good approximation of particle trajectory only when the particle locates on the center of the channel in depth. As the advantage of fluid-particle interaction method is achieved at a high computational cost and the streamline analysis is complex, a heuristic dimensionless design objective based on the Faxen''s law is proposed to optimize the VOA devices. The optimized performance of particle focusing is verified via the experiments and ALE method.