Stochastic simulations of DNA in flow: Dynamics and the effects of hydrodynamic interactions

Brownian dynamics simulations of bottlebrush polymers in dilute solution under simple shear and uniaxial extensional flows

We study the dynamics of bottlebrush polymer molecules in dilute solutions subjected to shear and uniaxial extensional flows using Brownian dynamics simulations with hydrodynamic interaction (HI). Bottlebrush polymers are modeled using a coarse-grained representation, consisting of a set of beads interacting pairwise via a purely repulsive potential and connected by finitely extensible nonlinear springs. We present the results for molecular stretching, stress, and solution viscosity during the startup of flow as well as under steady state as a function of side chain length while keeping the backbone length fixed. In extensional flow, the backbone fractional extension and the first normal stress difference decrease with an increase in side chain length at a fixed Weissenberg number (Wi). Using simulation results both in the presence of and in the absence of HI, we show that this is primarily a consequence of steric interaction resulting from the dense grafting of side chains. In shear flow, we observe a shear-thinning behavior in all cases, although it becomes less pronounced with increasing side chain length. Furthermore, nonmonotonicity in the backbone fractional extension is observed under shear, particularly at high Wi. We contextualize our simulation results for bottlebrush polymers with respect to existing studies in the literature for linear polymers and show that the unique dynamical features characterizing bottlebrush polymers arise on account of their additional molecular thickness due to the presence of densely grafted side chains.

Constructing custom thermodynamics using deep learning

One of the most exciting applications of artificial intelligence is automated scientific discovery based on previously amassed data, coupled with restrictions provided by known physical principles, including symmetries and conservation laws. Such automated hypothesis creation and verification can assist scientists in studying complex phenomena, where traditional physical intuition may fail. Here we develop a platform based on a generalized Onsager principle to learn macroscopic dynamical descriptions of arbitrary stochastic dissipative systems directly from observations of their microscopic trajectories. Our method simultaneously constructs reduced thermodynamic coordinates and interprets the dynamics on these coordinates. We demonstrate its effectiveness by studying theoretically and validating experimentally the stretching of long polymer chains in an externally applied field. Specifically, we learn three interpretable thermodynamic coordinates and build a dynamical landscape of polymer stretching, including the identification of stable and transition states and the control of the stretching rate. Our general methodology can be used to address a wide range of scientific and technological applications.

The power of weak, transient interactions across biology: A paradigm of emergent behavior

A growing list of diverse biological systems and their equally diverse functionalities provides realizations of a paradigm of emergent behavior. In each of these biological systems, pervasive ensembles of weak, short-lived, spatially local interactions act autonomously to convey functionalities at larger spatial and temporal scales. In this article, a range of diverse systems and functionalities are presented in a cursory manner with literature citations for further details. Then two systems and their properties are discussed in more detail: yeast chromosome biology and human respiratory mucus.

Characteristic features of self-avoiding active Brownian polymers under linear shear flow

We present Brownian dynamics simulation results of a flexible linear polymer with excluded-volume interactions under shear flow in the presence of active noise. The active noise strongly affects the polymer's conformational and dynamical properties, such as the stretching in the flow direction and compression in the gradient direction, shear-induced alignment, and shear viscosity. In the asymptotic limit of large activities and shear rates, the power-law scaling exponents of these quantities differ significantly from those of passive polymers. The chain's shear-induced stretching at a given shear rate is reduced by active noise, and it displays a non-monotonic behavior, where an initial polymer compression is followed by its stretching with increasing active force. The compression of the polymer in the gradient direction follows the relation ∼WiPe-3/4 as a function of the activity-dependent Weissenberg number WiPe, which differs from the scaling observed in passive systems ∼WiPe-1/2. The flow-induced alignment at large Péclet numbers Pe ≫ 1, where Pe is the Péclet number, and large shear rates WiPe ≫ 1 displays the scaling behavior WiPe-1/2, with an exponent differing from the passive value -1/3. Furthermore, the polymer's zero-shear viscosity displays a non-monotonic behavior, decreasing in an intermediate activity regime due to excluded-volume interactions and increasing again for large Pe. Shear thinning appears with increasing Weissenberg number with the power-laws WiPe-1/2 and WiPe-3/4 for passive and active polymers, respectively. In addition, our simulation results are compared with the results of an analytical approach, which predicts quantitatively similar behaviors for the various aforementioned physical quantities.

A computational model for structural dynamics and reconfiguration of DNA assemblies

Recent advances in constructing a structured DNA assembly whose configuration can be dynamically changed in response to external stimuli have demanded the development of an efficient computational modeling approach to expedite its design process. Here, we present a computational framework capable of analyzing both equilibrium and non-equilibrium dynamics of structured DNA assemblies at the molecular level. The framework employs Langevin dynamics with structural and hydrodynamic finite element models that describe mechanical, electrostatic, base stacking, and hydrodynamic interactions. Equilibrium dynamic analysis for various problems confirms the solution accuracy at a near-atomic resolution, comparable to molecular dynamics simulations and experimental measurements. Furthermore, our model successfully simulates a long-time-scale close-to-open-to-close dynamic reconfiguration of the switch structure in response to changes in ion concentration. We expect that the proposed model will offer a versatile way of designing responsive and reconfigurable DNA machines.

A Continuum Model for the Unfolding of von Willebrand Factor

von Willebrand Factor is a mechano-sensitive protein circulating in blood that mediates platelet adhesion to subendothelial collagen and platelet aggregation at high shear rates. Its hemostatic function and thrombogenic effect, as well as susceptibility to enzymatic cleavage, are regulated by a conformational change from a collapsed globular state to a stretched state. Therefore, it is essential to account for the conformation of the vWF multimers when modeling vWF-mediated thrombosis or vWF degradation. We introduce a continuum model of vWF unfolding that is developed within the framework of our multi-constituent model of platelet-mediated thrombosis. The model considers two interconvertible vWF species corresponding to the collapsed and stretched conformational states. vWF unfolding takes place via two regimes: tumbling in simple shear and strong unfolding in flows with dominant extensional component. These two regimes were demonstrated in a Couette flow between parallel plates and an extensional flow in a cross-slot geometry. The vWF unfolding model was then verified in several microfluidic systems designed for inducing high-shear vWF-mediated thrombosis and screening for von Willebrand Disease. The model predicted high concentration of stretched vWF in key regions where occlusive thrombosis was observed experimentally. Strong unfolding caused by the extensional flow was limited to the center axis or middle plane of the channels, whereas vWF unfolding near the channel walls relied upon the shear tumbling mechanism. The continuum model of vWF unfolding presented in this work can be employed in numerical simulations of vWF-mediated thrombosis or vWF degradation in complex geometries. However, extending the model to 3-D arbitrary flows and turbulent flows will pose considerable challenges.

Wall-anchored semiflexible polymer under large amplitude oscillatory shear flow

The properties of semiflexible polymers tethered by one end to an impenetrable wall and exposed to oscillatory shear flow are investigated by mesoscale simulations. A polymer, confined in two dimensions, is described by a linear bead-spring chain, and fluid interactions are incorporated by the Brownian multiparticle collision dynamics approach. At small strain, the polymers follow the applied flow field. However, at high strain, we find a strongly nonlinear response with major conformational changes. Polymers are stretched along the flow direction and exhibit U-shaped conformations while following the flow. As a consequence of confinement in the half-space, frequency doubling in the time-dependent polymer properties appears along the direction normal to the wall.

Polymer perspective of genome mobilization

Chromosome motion is an intrinsic feature of all DNA-based metabolic processes and is a particularly well-documented response to both DNA damage and repair. By using both biological and polymer physics approaches, many of the contributing factors of chromatin motility have been elucidated. These include the intrinsic properties of chromatin, such as stiffness, as well as the loop modulators condensin and cohesin. Various biological factors such as external tethering to nuclear domains, ATP-dependent processes, and nucleofilaments further impact chromatin motion. DNA damaging agents that induce double-stranded breaks also cause increased chromatin motion that is modulated by recruitment of repair and checkpoint proteins. Approaches that integrate biological experimentation in conjunction with models from polymer physics provide mechanistic insights into the role of chromatin dynamics in biological function. In this review we discuss the polymer models and the effects of both DNA damage and repair on chromatin motion as well as mechanisms that may underlie these effects.

A mechano-reactive coarse-grained model of the blood-clotting agent von Willebrand factor

The von Willebrand Factor (vWF) is a large blood glycoprotein that aids in hemostasis. Within each vWF monomer, the A2 domain hosts a cleavage site for enzyme ADAMTS13, which regulates the size of vWF multimers. This cleavage site can only be exposed when an A2 domain unfolds, and the unfolding reaction energy landscape is highly sensitive to the force conditions on the domain. Based on previous optical tweezer experimental results, we advance here a new activated A2 monomer model (AA2MM) for coarse-grained modeling of vWF that accurately represents the force-based probabilistic change between the unfolded/refolded states. A system of springs is employed to mimic the complex mechanical response of vWF monomers subject to pulling forces. AA2MM was validated by comparing monomer scale simulation results to data from prior pulling experiments on vWF monomer fragments. The model was further validated by comparing multimer scale Brownian dynamics simulation results to experiments using microfluidic chamber microscopy to visualize tethered vWF proteins subject to flow. The A2 domain unfolding reaction was studied in bulk flow simulations (pure shear and elongation flow), giving evidence that elongational flow drives the vWF size regulation process in blood. The mechanoreactive, coarse-grained AA2MM accurately describes the complex mechanical coupling between human blood flow conditions and vWF protein reactivity.

Simulation of semidilute polymer solutions in planar extensional flow via conformationally averaged Brownian noise

The dynamics and rheology of semidilute polymer solutions in strong flows are of great practical relevance. Processing applications can in principle be designed utilizing the relationship between nonequilibrium polymer conformations and the material properties of the solution. However, the interplay between concentration, flow, hydrodynamic interactions (HIs), and topological interactions which govern semidilute polymer dynamics is challenging to characterize. Brownian dynamics (BD) simulations are particularly valuable as a way to directly visualize how molecular interactions arise in these systems and are quantitatively comparable to single-molecule experiments. However, such simulations are often computationally intractable and are limited by the need to calculate the correlated Brownian noise via decomposition of the diffusion tensor. Previously, we have introduced an iterative conformational averaging (CA) method for BD simulations which bypasses these limitations by preaveraging the HI and Brownian noise in an iterative procedure. In this work, we generalize the CA method to flowing semidilute solutions by introducing a conformation dependent diffusion tensor and a strain dependent approximation to the conformationally averaged Brownian noise. We find that this approach nearly quantitatively reproduces both transient and steady state polymer dynamics and rheology while achieving an order of magnitude computational acceleration. We then utilize the CA method to investigate the concentration and flow rate dependence of polymer dynamics in planar extensional flows. Our results are consistent with previous experimental and simulation studies and provide a detailed view of broad conformational distributions in the semidilute regime. We observe interconversion between stretched and coiled states at steady state, which we conjecture occur due to the effect of concentration on the conformation dependent polymer drag. Additionally, we observe transient flow-induced intermolecular hooks in the startup of flow which lead to diverse and unique stretching pathways.

Shear-Induced Extensional Response Behaviors of Tethered von Willebrand Factor

We perform single-molecule flow experiments using confocal microscopy and a microfluidic device for shear rates up to 20,000 s^-1 and present results for the shear-induced unraveling and elongation of tethered von Willebrand factor (VWF) multimers. Further, we employ companion Brownian dynamics simulations to help explain details of our experimental observations using a parameterized coarse-grained model of VWF. We show that global conformational changes of tethered VWF can be accurately captured using a relatively simple mechanical model. Good agreement is found between experimental results and computational predictions for the threshold shear rate of extension, existence of nonhomogenous fluorescence distributions along unraveled multimer contours, and large variations in extensional response behaviors. Brownian dynamics simulations reveal the strong influence of varying chain length, tethering point location, and number of tethering locations on the underlying unraveling response. Through a complex molecule like VWF that naturally adopts a wide distribution of molecular size and has multiple binding sites within each molecule, this work demonstrates the power of tandem experiment and simulation for understanding flow-induced changes in biomechanical state and global conformation of macromolecules.

Tethered Semiflexible Polymer under Large Amplitude Oscillatory Shear

The properties of a semiflexible polymer with fixed ends exposed to oscillatory shear flow are investigated by simulations. The two-dimensionally confined polymer is modeled as a linear bead-spring chain, and the interaction with the fluid is described by the Brownian multiparticle collision dynamics approach. For small shear rates, the tethering of the ends leads to a more-or-less linear oscillatory response. However, at high shear rates, we found a strongly nonlinear reaction, with a polymer (partially) wrapped around the fixation points. This leads to an overall shrinkage of the polymer. Dynamically, the location probability of the polymer center-of-mass position is largest on a spatial curve resembling a limaçon, although with an inhomogeneous distribution. We found shear-induced modifications of the normal-mode correlation functions, with a frequency doubling at high shear rates. Interestingly, an even-odd asymmetry for the Cartesian components of the correlation functions appears, with rather similar spectra for odd x- and even y-modes and vice versa. Overall, our simulations yielded an intriguing nonlinear behavior of tethered semiflexible polymers under oscillatory shear flow.

Structure and proton conduction in sulfonated poly(ether ether ketone) semi-permeable membranes: a multi-scale computational approach

The design of polymeric membranes for proton or ionic exchange highly depends on the fundamental understanding of the physical and molecular mechanisms that control the formation of the conduction channels. There is an inherent relation between the dynamical structure of the polymeric membrane and the electrostatic forces that drive membrane segregation and proton transport. Here, we used a multi-scale computational approach to analyze the morphology of sulfonated poly(ether ether ketone) membranes at the mesoscale. A self-consistent description of the electrostatic phenomenon was adopted, where discrete polymer chains and a continuum proton field were embedded in a continuum fluid. Brownian dynamics was used for the evolution of the suspended polymer molecules, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the proton concentration field. We varied the polymer concentration, the degree of sulfonation and the level of confinement to find relationships between membrane structure and proton conduction. Our results indicate that the reduced mobility of polymer chains, at concentrations above overlap, and a moderate degree of sulfonation - i.e., 30% - are essential elements for membrane segregation and proton domain connectivity. These conditions also ensure that the membrane structure is not affected by size or by potential gradients. Importantly, our analysis shows that membrane conductivity and current are linearly dependent on polymer concentration and quadratically dependent on the degree of sulfonation. We found that the optimal polymeric membrane design requires a polymer concentration above overlap and a degree of sulfonation around 50%. These conditions promote a dynamical membrane morphology with a constant density of proton channels. Our results and measurements agree with previous experimental works, thereby validating our model and observations.

Polymer margination in uniform shear flows

We address the issue of polymer margination (migration towards surfaces) in uniform shear flows through extensive LBMD (lattice-Boltzmann molecular dynamics) simulations. In particular we consider the effect of monomer size, a on the chain's overall margination tendency for chains of length N = 16, 32 monomers in flows at multiple shear rates [small gamma, Greek, dot above]. We observed higher margination of chains with larger radii monomers in comparison to smaller radii monomer chains of the same length N. We quantify this effect by considering various measures such as the distribution of the maximum extent of the chain into the channel bulk, zm, distribution of its center of mass in the direction normal to the surface, zc and the distributions of the chain's radius of gyration in directions parallel and perpendicular to the surface i.e. Rx, Ry and Rz respectively.

Modelling DNA extension and fragmentation in contractive microfluidic devices: a Brownian dynamics and computational fluid dynamics approach

Fragmenting DNA into short pieces is an essential manipulation in many biological studies, ranging from genome sequencing to molecular diagnosis. Among various DNA fragmentation methods, microfluidic hydrodynamic DNA fragmentation has huge advantages especially in terms of handling small-volume samples and being integrated into automatic and all-in-one DNA analysis equipment. Despite the fast progress in experimental studies and applications, a systematic understanding of how DNA molecules are distributed, stretched and fragmented in a confined microfluidic field is still lacking. In this work, we investigate the extension and fragmentation of DNA in a typical contractive microfluidic field, which consists of a shear flow-dominated area and an elongational flow-dominated area, using the Brownian dynamics-computational fluid dynamics method. Our results show that the shear flow at the straight part of the microfluidic channel and the elongational flow at the contractive bottleneck together determine the performance of DNA fragmentation. The average fragment size of DNA decreases with the increase of the strain rate of the elongational flow, and the upstream shear flow can significantly precondition the conformation of DNA to produce shorter and more uniform fragments. A systematic study of the dynamics of DNA fragmentation shows that DNA tends to break at the mid-point when the strain rate of elongational flow is small, and the breakage point largely deviates from the midpoint as the strain rate increases. Our simulation of the thorough DNA fragmentation process in a realistic microfluidic field agrees well with experimental results. We expect that our study can shed new light on the development of future microfluidic devices for DNA fragmentation and integrated DNA analysis devices.

Internal Tensile Force and A2 Domain Unfolding of von Willebrand Factor Multimers in Shear Flow

Using Brownian molecular dynamics simulations, we examine the internal dynamics and biomechanical response of von Willebrand factor (vWF) multimers subject to shear flow. The coarse grain multimer description employed here is based on a monomer model in which the A2 domain of vWF is explicitly represented by a nonlinear elastic spring whose mechanical response was fit to experimental force/extension data from vWF monomers. This permits examination of the dynamic behavior of hydrodynamic forces acting on A2 domains as a function of shear rate and multimer length, as well as position of an A2 domain along the multimer contour. Force/position data reveal that collapsed multimers exhibit a force distribution with two peaks, one near each end of the chain; unraveled multimers, however, show a single peak in A2 domain force near the center of multimers. Guided further by experimental data, significant excursions of force acting on a domain are associated with an increasing probability for A2 domain unfolding. Our results suggest that the threshold shear rate required to induce A2 domain unfolding is inversely proportional to multimer length. By examining data for the duration and location of significant force excursions, convincing evidence is advanced that unfolding of A2 domains, and therefore scission of vWF multimers by the size-regulating blood enzyme ADAMTS13, happen preferentially near the center of unraveled multimers.

Optimization of flow-focusing devices for homogeneous extensional flow

We present a methodology for the shape optimization of flow-focusing devices with the purpose of creating a wide region of homogeneous extensional flow, characterized by a uniform strain-rate along the centerline of the devices. The numerical routines employed include an optimizer, a finite-volume solver, and a mesh generator operating on geometries with the walls parameterized by Bézier curves. The optimizations are carried out for devices with different geometric characteristics (channel aspect ratio and length). The performance of the optimized devices is assessed for varying Reynolds numbers, velocity ratio between streams, and fluid rheology. Brownian dynamics simulations are also performed to evaluate the stretching and relaxation of λ-DNA molecules in the devices. Overall, the optimized flow-focusing devices generate a homogeneous extensional flow over a range of conditions typically found in microfluidics. At high Weissenberg numbers, the extension of λ-DNA molecules in the optimized flow-focusing devices is close to that obtained in an ideal planar extensional flow with an equivalent Hencky strain. The devices presented in this study can be useful in microfluidic applications taking advantage of homogeneous extensional flows and easy control of the Hencky strain and strain-rate.

Active Brownian Filamentous Polymers under Shear Flow

The conformational and rheological properties of active filaments/polymers exposed to shear flow are studied analytically. Using the continuous Gaussian semiflexible polymer model extended by the activity, we derive analytical expressions for the dependence of the deformation, orientation, relaxation times, and viscosity on the persistence length, shear rate, and activity. The model yields a Weissenberg-number dependent shear-induced deformation, alignment, and shear thinning behavior, similarly to the passive counterpart. Thereby, the model shows an intimate coupling between activity and shear flow. As a consequence, activity enhances the shear-induced polymer deformation for flexible polymers. For semiflexible polymers/filaments, a nonmonotonic deformation is obtained because of the activity-induced shrinkage at moderate and swelling at large activities. Independent of stiffness, activity-induced swelling facilitates and enhances alignment and shear thinning compared to a passive polymer. In the asymptotic limit of large activities, a polymer length- and stiffness-independent behavior is obtained, with universal shear-rate dependencies for the conformations, dynamics, and rheology.

Multi-resolution polymer Brownian dynamics with hydrodynamic interactions

A polymer model given in terms of beads, interacting through Hookean springs and hydrodynamic forces, is studied. A Brownian dynamics description of this bead-spring polymer model is extended to multiple resolutions. Using this multiscale approach, a modeller can efficiently look at different regions of the polymer in different spatial and temporal resolutions with scalings given for the number of beads, statistical segment length, and bead radius in order to maintain macro-scale properties of the polymer filament. The Boltzmann distribution of a Gaussian chain for differing statistical segment lengths gives a diffusive displacement equation for the multi-resolution model with a mobility tensor for different bead sizes. Using the pre-averaging approximation, the translational diffusion coefficient is obtained as a function of the inverse of a matrix and then in closed form in the long-chain limit. This is then confirmed with numerical experiments.

Application of the Eckart frame to soft matter: rotation of star polymers under shear flow

The Eckart co-rotating frame is used to analyze the dynamics of star polymers under shear flow, either in melt or solution and with different types of bonds. This formalism is compared with the standard approach used in many previous studies on polymer dynamics, where an apparent angular velocity ω is obtained from the relation between the tensor of inertia and angular momentum. A common mistake is to interpret ω as the molecular rotation frequency, which is only valid for rigid-body rotation. The Eckart frame, originally formulated to analyze the infrared spectra of small molecules, dissects different kinds of displacements: vibrations without angular momentum, pure rotation, and vibrational angular momentum (leading to a Coriolis cross-term). The Eckart frame co-rotates with the molecule with an angular frequency Ω obtained from the Eckart condition for minimal coupling between rotation and vibration. The standard and Eckart approaches are compared with a straight description of the star's dynamics taken from the time autocorrelation of the monomer positions moving around the molecule's center of mass. This is an underdamped oscillatory signal, which can be described by a rotation frequency ω_R and a decorrelation rate Γ. We consistently find that Ω coincides with ω_R, which determines the characteristic tank-treading rotation of the star. By contrast, the apparent angular velocity ω < Ω does not discern between pure rotation and molecular vibrations. We believe that the Eckart frame will be useful to unveil the dynamics of semiflexible molecules where rotation and deformations are entangled, including tumbling, tank-treading motions and breathing modes.

Deciphering the dynamics of star molecules in shear flow

This work analyses the rotation of star polymers under shear flow, in melts, and in good solvent dilute solution. The latter is modeled by single molecule Brownian hydrodynamics, while melts are modeled using non-equilibrium molecular dynamics in closed (periodic) boxes and in open boundaries. A Dissipative Particle Dynamics (DPD) thermostat introduces pairwise monomer friction in melts at will, in directions normal and tangent to the monomer-monomer vectors. Although tangential friction is seldom modeled, we show that it is essential to control hydrodynamic effects in melts. We analyze the different sources of molecular angular momentum in solution and melts and distinguish three dynamic regimes as the shear rate [small gamma, Greek, dot above] is increased. These dynamic regimes are related with the disruption of the different relaxation mechanisms of the star in equilibrium. Although strong differences are found between harmonic springs and finitely extensible bonds, above a critical shear rate the star molecule has a "breathing" mode with successive elongations and contractions in the flow direction with frequency Ω. The force balance in the flow direction unveils a relation between Ω and the orientation angle. Using literature results for the tumbling of rings and linear chains, either in melt or in solution, we show that the relation is general. A different "tank-treading" dynamics determines the rotation of monomers around the center of mass of the molecule. We show that the tank-treading frequency does not saturate but keeps increasing with [small gamma, Greek, dot above]. This is at odds with previous studies which erroneously calculated the molecular angular frequency, used as a proxy for tank-treading.

Modeling the relaxation of internal DNA segments during genome mapping in nanochannels

We have developed a multi-scale model describing the dynamics of internal segments of DNA in nanochannels used for genome mapping. In addition to the channel geometry, the model takes as its inputs the DNA properties in free solution (persistence length, effective width, molecular weight, and segmental hydrodynamic radius) and buffer properties (temperature and viscosity). Using pruned-enriched Rosenbluth simulations of a discrete wormlike chain model with circa 10 base pair resolution and a numerical solution for the hydrodynamic interactions in confinement, we convert these experimentally available inputs into the necessary parameters for a one-dimensional, Rouse-like model of the confined chain. The resulting coarse-grained model resolves the DNA at a length scale of approximately 6 kilobase pairs in the absence of any global hairpin folds, and is readily studied using a normal-mode analysis or Brownian dynamics simulations. The Rouse-like model successfully reproduces both the trends and order of magnitude of the relaxation time of the distance between labeled segments of DNA obtained in experiments. The model also provides insights that are not readily accessible from experiments, such as the role of the molecular weight of the DNA and location of the labeled segments that impact the statistical models used to construct genome maps from data acquired in nanochannels. The multi-scale approach used here, while focused towards a technologically relevant scenario, is readily adapted to other channel sizes and polymers.

Entropy gives rise to topologically associating domains

We investigate chromosome organization within the nucleus using polymer models whose formulation is closely guided by experiments in live yeast cells. We employ bead-spring chromosome models together with loop formation within the chains and the presence of nuclear bodies to quantify the extent to which these mechanisms shape the topological landscape in the interphase nucleus. By investigating the genome as a dynamical system, we show that domains of high chromosomal interactions can arise solely from the polymeric nature of the chromosome arms due to entropic interactions and nuclear confinement. In this view, the role of bio-chemical related processes is to modulate and extend the duration of the interacting domains.

The polymer physics of single DNA confined in nanochannels

In recent years, applications and experimental studies of DNA in nanochannels have stimulated the investigation of the polymer physics of DNA in confinement. Recent advances in the physics of confined polymers, using DNA as a model polymer, have moved beyond the classic Odijk theory for the strong confinement, and the classic blob theory for the weak confinement. In this review, we present the current understanding of the behaviors of confined polymers while briefly reviewing classic theories. Three aspects of confined DNA are presented: static, dynamic, and topological properties. The relevant simulation methods are also summarized. In addition, comparisons of confined DNA with DNA under tension and DNA in semidilute solution are made to emphasize universal behaviors. Finally, an outlook of the possible future research for confined DNA is given.

Mathematical model suitable for efficient simulation of thin semi-flexible polymers in complex environments

We present an alternative approach to simulations of semi-flexible polymers. In contrast with the usual bead-rod compromise between bead-spring and rigid rod models, we use deformable cylindrical segments as basic units of the polymer. The length of each segment is not preserved with end points diffusing under constraints keeping the polymer chain nature intact. The model allows the simulation of tension transport and elasticity properties. In particular we describe a new cooperative regime in the relaxation of the polymer from its fully elongated configuration.

Passive trapping of rigid rods due to conformation-dependent electrophoretic mobility

We present computer simulations of a rigid rod in a combination of an extensional fluid flow and extensional electric field. The electrophoretic mobility of the rod is different parallel or perpendicular to the rod. The dependence of the mobility on the conformation (orientation) leads to a new phenomenon where the rods can be passively trapped in all directions at the stagnation point. This contrasts with the behavior in either fluid flow or electric field alone, in which an object can be pushed towards the stagnation point along some directions but is pushed away in others. We have determined the state space where trapping occurs and have developed a model that describes the strength of trapping when it does occur. This new phenomenon could be used in the future to separate objects based on a coupling between their mobility and ability to be oriented.

Dimensional reduction of duplex DNA under confinement to nanofluidic slits

There has been much interest in the dimensional properties of double-stranded DNA (dsDNA) confined to nanoscale environments as a problem of fundamental importance in both biological and technological fields. This has led to a series of measurements by fluorescence microscopy of single dsDNA molecules under confinement to nanofluidic slits. Despite the efforts expended on such experiments and the corresponding theory and simulations of confined polymers, a consistent description of changes of the radius of gyration of dsDNA under strong confinement has not yet emerged. Here, we perform molecular dynamics (MD) simulations to identify relevant factors that might account for this inconsistency. Our simulations indicate a significant amplification of excluded volume interactions under confinement at the nanoscale due to the reduction of the effective dimensionality of the system. Thus, any factor influencing the excluded volume interaction of dsDNA, such as ionic strength, solution chemistry, and even fluorescent labels, can greatly influence the dsDNA size under strong confinement. These factors, which are normally less important in bulk solutions of dsDNA at moderate ionic strengths because of the relative weakness of the excluded volume interaction, must therefore be under tight control to achieve reproducible measurements of dsDNA under conditions of dimensional reduction. By simulating semi-flexible polymers over a range of parameter values relevant to the experimental systems and exploiting past theoretical treatments of the dimensional variation of swelling exponents and prefactors, we have developed a novel predictive relationship for the in-plane radius of gyration of long semi-flexible polymers under slit-like confinement. Importantly, these analytic expressions allow us to estimate the properties of dsDNA for the experimentally and biologically relevant range of contour lengths that is not currently accessible by state-of-the-art MD simulations.

Coarse-grained model of conformation-dependent electrophoretic mobility and its influence on DNA dynamics

The electrophoretic mobility of molecules such as λ-DNA depends on the conformation of the molecule. It has been shown that electrohydrodynamic interactions between parts of the molecule lead to a mobility that depends on conformation and can explain some experimental observations. We have developed a new coarse-grained model that incorporates these changes of mobility into a bead-spring chain model. Brownian dynamics simulations have been performed using this model. The model reproduces the cross-stream migration that occurs in capillary electrophoresis when pressure-driven flow is applied parallel or antiparallel to the electric field. The model also reproduces the change of mobility when the molecule is stretched significantly in an extensional field. We find that the conformation-dependent mobility can lead to a new type of unraveling of the molecule in strong fields. This occurs when different parts of the molecule have different mobilities and the electric field is large.

Hydrodynamic interactions of deformable polymeric nanocarriers and the effect of crosslinking

We report theoretical as well as numerical investigations of deformable nanocarriers (NCs) under physiologically relevant flow conditions. Specifically, to model the deformable lysozyme-core/dextran-shell crosslinked polymer based NC with internal nanostructure and subject it to external hydrodynamic shear, we have introduced a coarse-grained model for the NC and have adopted a Brownian dynamics framework, which incorporates hydrodynamic interactions, in order to describe the static and dynamic properties of the NC. In order to represent the fluidity of the polymer network in the dextran brush-like corona, we coarse-grain the structure of the NC based on the hypothesis that Brownian motion, polymer melt reptations, and crosslinking density dominate their structure and dynamics. In our model, we specify a crosslinking density and employ the simulated annealing protocol to mimic the experimental synthesis steps in order to obtain the appropriate internal structure of the core-shell polymer. We then compute the equilibrium as well as steady shear rheological properties as functions of the Péclet number and the crosslinking density, in the presence of hydrodynamic interactions. We find that with increasing crosslinking, the stiffness of the nanocarrier increases, the radius of gyration decreases, and as a consequence the self-diffusivity increases. The nanocarrier under shear deforms and orients along the direction of the applied shear and we find that the orientation and deformation under shear are dependent on the shear rate and the crosslinking density. We compare various dynamic properties of the NC as a function of the shear force, such as orientation, deformation, intrinsic stresses etc., with previously reported computational and experimental results of other model systems. The computational approach described here serves as a powerful tool for the rational design of NCs by taking both the physiological as well as the hydrodynamic environments into consideration. Development of such models is essential in order to gain useful insights that may be translated into the optimal design of NCs for diagnostic as well as targeted drug delivery applications.

Effects of excluded volume and hydrodynamic interaction on the deformation, orientation and motion of ring polymers in shear flow

A ring polymer is a classical model to explore the behaviors of biomacromolecules. Compared with its linear counterpart in shear flow, the ring polymer should be more sensitive to excluded volume and hydrodynamic interaction attributed to the absence of chain ends. We carried out multiparticle collision dynamics combined with molecular dynamics simulation to study the effects of excluded volume and hydrodynamic interaction on the behaviors of ring polymers in shear flow. The results show that in the absence of the strong excluded volume interaction, the ring polymer prefers a two-strand linear conformation with high deformation and orientation in the flow-gradient plane, and the tank-treading motion is nearly negligible. Ring polymers without excluded volume show no significant difference from linear polymers in the scaling exponents for the deformation, orientation and tumbling motion. We also observed that the hydrodynamic interaction could efficiently slow down the relaxation of ring polymers while the scaling exponents against the Weissenberg number have rarely been affected.

Self-consistent description of electrokinetic phenomena in particle-based simulations

A new computational method is presented for study suspensions of charged particles undergoing fluctuating hydrodynamic and electrostatic interactions. The proposed model is appropriate for polymers, proteins, and porous particles embedded in a continuum electrolyte. A self-consistent Langevin description of the particles is adopted in which hydrodynamic and electrostatic interactions are included through a Green's function formalism. An Ewald-like split is adopted in order to satisfy arbitrary boundary conditions for the Stokeslet and Poisson Green functions, thereby providing a formalism that is applicable to any geometry and that can be extended to deformable objects. The convection-diffusion equation for the continuum ions is solved simultaneously considering Nernst-Planck diffusion. The method can be applied to systems at equilibrium and far from equilibrium. Its applicability is demonstrated in the context of electrokinetic motion, where it is shown that the ionic clouds associated with individual particles can be severely altered by the flow and concentration, leading to intriguing cooperative effects.

Effects of thermal noise on the transitional dynamics of an inextensible elastic filament in stagnation flow

We investigate the dynamics of a single inextensible elastic filament subject to anisotropic friction in a viscous stagnation-point flow, by employing both a continuum model represented by Langevin type stochastic partial differential equations (SPDEs) and a dissipative particle dynamics (DPD) method. Unlike previous works, the filament is free to rotate and the tension along the filament is determined by the local inextensible constraint. The kinematics of the filament is recorded and studied with normal modes analysis. The results show that the filament displays an instability induced by negative tension, which is analogous to Euler buckling of a beam. Symmetry breaking of normal modes dynamics and stretch-coil transitions are observed above the threshold of the buckling instability point. Furthermore, both temporal and spatial noise are amplified resulting from the interaction of thermal fluctuations and nonlinear filament dynamics. Specifically, the spatial noise is amplified with even normal modes being excited due to symmetry breaking, while the temporal noise is amplified with increasing time correlation length and variance.

Efficient Brownian dynamics simulation of DNA molecules with hydrodynamic interactions in linear flows

The coarse-grained molecular dynamics (MD) or Brownian dynamics (BD) simulation is a particle-based approach that has been applied to a wide range of biological problems that involve interactions with surrounding fluid molecules or the so-called hydrodynamic interactions (HIs). In this paper, an efficient algorithm is proposed to simulate the motion of a single DNA molecule in linear flows. The algorithm utilizes the integrating factor to cope with the effect of the linear flow of the surrounding fluid and applies the Metropolis method (MM) by Bou-Rabee, Donev, and Vanden-Eijnden [Multiscale Model. Simul. 12, 781 (2014)] to achieve more efficient BD simulation. Thus our method permits much larger time step size than previous methods while still maintaining the stability of the BD simulation, which is advantageous for long-time BD simulation. Our numerical results on λ-DNA agree very well with both experimental data and previous simulation results. Finally, when combined with fast algorithms such as the fast multipole method which has nearly optimal complexity in the total number of beads, the resulting method is parallelizable, scalable to large systems, and stable for large time step size, thus making the long-time large-scale BD simulation within practical reach. This will be useful for the study of membranes, long-chain molecules, and a large collection of molecules in the fluids.

Effect of hydrodynamic correlations on the dynamics of polymers in dilute solution

We analyze the effect of time-dependent hydrodynamic interactions on the dynamics of flexible polymers in dilute solution. In analytical calculations, the fluctuating hydrodynamics approach is adopted to describe the fluid, and a Gaussian model to represented the polymer. Simulations are performed exploiting the multiparticle collision dynamics approach, a mesoscale hydrodynamic simulation technique, to explicitly describe the fluid. Polymer center-of-mass velocity correlation functions are calculated for various polymer lengths. Similarly, segment mean square displacements are discussed and polymer diffusion coefficients are determined. Particular attention is paid to the influence of sound propagation on the various properties. The simulations reveal a strong effect of hydrodynamic interactions. Specifically, the time dependence of the center-of-mass velocity correlation functions is determined by polymer properties over a length-dependent time window, but are asymptotically solely governed by fluid correlations, with a long-time tail decaying as t(-3/2). The correlation functions are heavily influenced by sound modes for short polymers, an effect which gradually disappears with increasing polymer length. We find excellent agreement between analytical and simulation results. This allows us to provide a theory-based asymptotic value for the polymer diffusion coefficient in the limit of large system sizes, which is based on a single finite-system-size simulation.

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 Duplex Formation with a Coarse-Grained Model

A middle-resolution coarse-grained model of DNA is proposed. The DNA chain is built of spherical and planar rigid bodies connected by elastic virtual bonds. The bonded part of the potential energy function is fit to potentials of mean force of model systems. The rigid bodies are sets of neutral, charged, and dipolar beads. Electrostatic and van der Waals interactions are parametrized by our recently developed procedure [Maciejczyk, M.; Spasic, A.; Liwo, A.; Scheraga, H.A. J. Comp. Chem.2010, 31, 1644]. Interactions with the solvent and an ionic cloud are approximated by a multipole-multipole Debye-Hückel model. A very efficient R-RATTLE algorithm, for integrating the movement of rigid bodies, is implemented. It is the first coarse-grained model, in which both bonded and nonbonded interactions were parametrized ab initio and which folds stable double helices from separated complementary strands, with the final conformation close to the geometry of experimentally determined structures.

Nonequilibrium thermodynamics of dilute polymer solutions in flow

Modern materials processing applications and technologies often occur far from equilibrium. To this end, the processing of complex materials such as polymer melts and nanocomposites generally occurs under strong deformations and flows, conditions under which equilibrium thermodynamics does not apply. As a result, the ability to determine the nonequilibrium thermodynamic properties of polymeric materials from measurable quantities such as heat and work is a major challenge in the field. Here, we use work relations to show that nonequilibrium thermodynamic quantities such as free energy and entropy can be determined for dilute polymer solutions in flow. In this way, we determine the thermodynamic properties of DNA molecules in strong flows using a combination of simulations, kinetic theory, and single molecule experiments. We show that it is possible to calculate polymer relaxation timescales purely from polymer stretching dynamics in flow. We further observe a thermodynamic equivalence between nonequilibrium and equilibrium steady-states for polymeric systems. In this way, our results provide an improved understanding of the energetics of flowing polymer solutions.

An experimentally-informed coarse-grained 3-Site-Per-Nucleotide model of DNA: structure, thermodynamics, and dynamics of hybridization

A new 3-Site-Per-Nucleotide coarse-grained model for DNA is presented. The model includes anisotropic potentials between bases involved in base stacking and base pair interactions that enable the description of relevant structural properties, including the major and minor grooves. In an improvement over available coarse-grained models, the correct persistence length is recovered for both ssDNA and dsDNA, allowing for simulation of non-canonical structures such as hairpins. DNA melting temperatures, measured for duplexes and hairpins by integrating over free energy surfaces generated using metadynamics simulations, are shown to be in quantitative agreement with experiment for a variety of sequences and conditions. Hybridization rate constants, calculated using forward-flux sampling, are also shown to be in good agreement with experiment. The coarse-grained model presented here is suitable for use in biological and engineering applications, including nucleosome positioning and DNA-templated engineering.