Brownian dynamics simulation of DNA in complex geometries.

摘要

This dissertation is concerned with the dynamics of a long DNA molecule in complex geometries, driven by either electrostatic field or flow field. This is accomplished primarily through the use of Brownian dynamics simulation, which captures the essential physics at mesoscopic length scale, and allows us to simulate events happening on long time scale, such as DNA pore translocation and cyclic dynamics of a tethered DNA molecule in shear flow. General methods are developed for both electric field-driven and flow-driven DNA dynamics in complex geometries. In the electric field-driven case, we propose a novel class of electric field-actuated soft mechanical control element for microfluidics. Brownian dynamics (BD)/Finite Element Method (FEM) simulation efficiently explore the design space, and results demonstrate that the On/Off switching could be achieved within a proper parameter space. In the flow-driven case, we first examine the cyclic dynamics of a single DNA molecule tethered to a hard wall in shear flow. Extensive simulation results suggest a classical fluctuation-dissipation stochastic process and question the existence of periodicity of the cyclic dynamics, as previously claimed. We support our numerical findings with a simple analytical calculation for a harmonic dimer in shear flow. In the case of flow-driven DNA molecule in complex geometries, one big challenge is that an efficient algorithm is required to calculate fluctuating hydrodynamic interactions (HI) in complex geometries. We have developed an accelerated immersed boundary method that allows fast calculation of Brownian motion of polymer chains and other particles in complex geometries with HI. With this new method, the first detailed analysis of a recent set of interesting nanofluidic experiments involving DNA dynamics in a complex flow geometry is performed. This analysis explains the observed dynamics over a wide range of parameter values (flow rate, molecular wieght) and illustrates the important quantitative effect of the hydrodynamic interactions on the behavior of the system.