Dynamics of collapsed polymers under the simultaneous influence of elongational and shear flows

Predicting reaction behavior of tethered polymers in shear flow

Kinetics of force-mediated chemical reactions of end-tethered polymers with varying chain length N in varying shear rate flow γ̇ are explored via coarse-grained Brownian dynamics simulations. At fixed γ̇, force F along a polymer increases linearly with N as previously predicted; however, contrary to existing theory, the F(N) slope increases for N above a transition length that exhibits minimal dependence on γ̇. Force profiles are used in a stochastic model of a force-mediated reaction to compute the time for x percent of a polymer population to experience a reaction, tx. Observations are insensitive to the selected value of x in that tx data for varying N and γ̇ can be consistently collapsed onto a single curve via appropriate scaling, with one master curve for systems below the transition N (small N) and another for those above (large N). Different force scaling for small and large N results in orders of magnitude difference in force-mediated reaction kinetics as represented by the population response time. Data presented illustrate the possibility of designing mechano-reactive polymer populations with highly controlled response to flow across a range in γ̇.

Enhanced Pulley Effect for Translocation: The Interplay of Electrostatic and Hydrodynamic Forces

Solid-state nanopore sensors remain a promising solution to the rising global demand for genome sequencing. These single-molecule sensing technologies require single-file translocation for high resolution and accurate detection. In a previous publication, we discovered a hairpin unraveling mechanism, namely, the pulley effect, in a pressure-driven translocation system. In this paper, we further investigate the pulley effect in the presence of pressure-driven fluid flow and an opposing force provided by an electrostatic field as an approach to increase single-file capture probability. A hydrodynamic flow is used to move the polymer forward, and two oppositely charged electrostatic square loops are used to create an opposing force. By optimizing the balance between forces, we show that the single-file capture can be amplified from about 50% to almost 95%. The force location, force strength, and flow rate are used as the optimizing variables.

Statistical Dynamics of Flow-Driven Globular Polymers

The response of collapsed polymers to the effects of linear mixed flow is studied theoretically in this paper using a model of a self-interacting finitely extensible Gaussian chain that evolves stochastically in the presence of random thermal fluctuations and an external fluid velocity gradient. The interactions that produce compact chain configurations are described by a harmonic pair potential of strength κ that acts between nonbonded sites on the chain backbone. Several chain properties are calculated analytically from this model as a function of κ for elongational and shear flows, including the dependence of the chain's steady-state mean-square end-to-end distance on the Weissenberg number of the flow, the time-dependence of the chain's relaxation to equilibrium from a steady-state of given chain extension, and the nature of the force-extension curves that are obtained from the free energy change between unperturbed and flow-stretched states of the chain. For both elongational and shear flows (but to different degrees), it is found that the greater the value of κ (and the more compact the chain), the more difficult it is, in general, for the imposed flow to induce a transition between compact and extended states, in broad agreement with available data from numerical simulations. For the relaxation process, the differences between the two flow types are more marked. The characteristic decay time for relaxation from a state prepared by elongational flow is essentially independent of κ, whereas in the case of a state prepared by shear flow, it is distinctly κ-dependent, the relaxation becoming faster at larger κ.

Uncertainty quantification of a thrombosis model considering the clotting assay PFA-100®

Mathematical models of thrombosis are currently used to study clinical scenarios of pathological thrombus formation. As these models become more complex to predict thrombus formation dynamics high computational cost must be alleviated and inherent uncertainties must be assessed. Evaluating model uncertainties allows to increase the confidence in model predictions and identify avenues of improvement for both thrombosis modeling and anti-platelet therapies. In this work, an uncertainty quantification analysis of a multi-constituent thrombosis model is performed considering a common assay for platelet function (PFA-100®). The analysis is facilitated thanks to time-evolving polynomial chaos expansions used as a parametric surrogate for the full thrombosis model considering two quantities of interest; namely, thrombus volume and occlusion percentage. The surrogate is thoroughly validated and provides a straightforward access to a global sensitivity analysis via computation of Sobol' coefficients. Six out of 15 parameters linked to thrombus consitution, vWF activity, and platelet adhesion dynamics were found to be most influential in the simulation variability considering only individual effects; while parameter interactions are highlighted when considering the total Sobol' indices. The influential parameters are related to thrombus constitution, vWF activity, and platelet to platelet adhesion dynamics. The surrogate model allowed to predict realistic PFA-100® closure times of 300,000 virtual cases that followed the trends observed in clinical data. The current methodology could be used including common anti-platelet therapies to identify scenarios that preserve the hematological balance.

Flow-regulated nucleation protrusion theory for collapsed polymers

The globular-stretch transition of a collapsed polymer in low strain rate elongational flow is studied using polymeric protrusion kinetics scaling laws and numerical simulation. Results demonstrate the influence of fluid flow on the occurrence probability of long-length thermally nucleated polymeric protrusions, which regulate collapsed polymer unfolding in low strain rate flows. Further, we estimate that the globular-stretch transition rate (k_{s}) in low strain rate (∈[over ̇]) elongational flows varies as k_{s}∼e^{-α∈[over ̇]^{-1}}. Results here reveal that the existing approach of neglecting the effects of fluid flow on thermally nucleated protrusions distribution is not valid for analyzing polymer unfolding behavior in low strain rate flows. Neglecting such an effect overestimates the constant α in the scaling law of transition rate (k_{s}∼e^{-α∈[over ̇]^{-1}}) by a factor of 2.

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.

Predicting pathological von Willebrand factor unraveling in elongational flow

The globular-to-unraveled conformation transition of von Willebrand factor (vWF), a large polymeric glycoprotein in human blood plasma, is a crucial step in the process of clotting at sites of vascular injury. However, unraveling of vWF multimers in uninjured vasculature can lead to pathology (i.e., thrombus formation or degradation of vWF proteins by enzyme ADAMTS13, making them nonfunctional). To identify blood flow conditions that might induce pathological unraveling of vWF multimers, here we have computed the globular-to-unraveled transition rate of vWF multimers subjected to varying strain rate elongational flow by employing an enhanced sampling technique, the weighted ensemble method. Weighted ensemble sampling was employed instead of standard brute-force simulations because pathological blood flow conditions can induce undesired vWF unraveling on timescales potentially inaccessible to standard simulation methods. Results here indicate that brief but periodic exposure of vWF to the elongational flow of strain rate greater than or equal to 2500 s^-1 represents a source of possible pathology caused by the undesired unraveling of vWF multimers.

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.

Coarse-Grain Modeling of Shear-Induced Binding between von Willebrand Factor and Collagen

Von Willebrand factor (VWF) is a large multimeric protein that aids in blood clotting. Near injury sites, hydrodynamic force from increased blood flow elongates VWF, exposing binding sites for platelets and collagen. To investigate VWF binding to collagen that is exposed on injured arterial surfaces, Brownian dynamics simulations are performed with a coarse-grain molecular model. Accounting for hydrodynamic interactions in the presence of a stationary surface, shear flow conditions are modeled. Binding between beads in coarse-grain VWF and collagen sites on the surface is described via reversible ligand-receptor-type bond formation, which is governed via Bell model kinetics. For conditions in which binding is energetically favored, the model predicts a high probability for binding at low shear conditions; this is counter to experimental observations but in agreement with what prior modeling studies have revealed. To address this discrepancy, an additional binding criterion that depends on the conformation of a submonomer feature in the model local to a given VWF binding site is implemented. The modified model predicts shear-induced binding, in very good agreement with experimental observations; this is true even for conditions in which binding is significantly favored energetically. Biological implications of the model modification are discussed in terms of mechanisms of VWF activity.

Ring polymer dynamics and tumbling-stretch transitions in planar mixed flows

The properties of dilute polymer solutions are governed by the conformational dynamics of individual polymers which can be perturbed in the presence of an applied flow. Much of our understanding of dilute solutions comes from studying how flows manipulate the molecular features of polymer chains out of equilibrium, primarily focusing on linear polymer chains. Recently there has been an emerging interest in the dynamics of nonlinear architectures, particularly ring polymers, which exhibit surprising out-of-equilibrium dynamics in dilute solutions. In particular, it has been observed that hydrodynamics can couple to topology in planar elongational and shear flows, driving molecular expansion in the nonflow direction that is not observed for linear chains. In this paper, we extend our understanding of dilute ring polymer dynamics to mixed flows, which represent flow profiles intermediate between simple shear or planar elongation. We map the conformational behaviors at a number of flow geometries and strengths, demonstrating transitions between coiled, tumbling, and stretched regimes. Indeed, these observations are consistent with how linear chains respond to mixed flows. For both linear and ring polymers, we observe a marked first-order-like transition between tumbling and stretched polymers that we attribute to a dynamic energy barrier between the two states. This manifests as bimodal extension distributions in a narrow range of flow strengths and geometries, with the primary difference between rings and linear chains being the presence of molecular expansion in the vorticity direction.

Bioinspired Shear-Flow-Driven Layer-by-Layer in Situ Self-Assembly

Layer-by-layer (LbL) assembly is widely applied as a coating technique for the nanoscale control of architecture and related properties. However, its translational applications are limited by the time-consuming and laborious nature of the process. Inspired by the blood-clotting process, herein, we develop a shear-flow-driven LbL (SF-LbL) self-assembly approach that accelerates the adsorption rate of macromolecules by mechanically configuring the polymer chain via a coil-stretch transition, which effectively simplifies and speeds the diffusion-controlled assembly process. The structural characteristics and surface homogeneity of the SF-LbL films are improved, and diverse three-dimensional structures can be achieved. Functional SF-LbL-assembled surfaces for corneal modification are successfully fabricated, and the surface of wounded rat corneas and skin can be directly decorated in situ with SF-LbL nanofilms due to the advantages of this approach. Furthermore, in situ SF-LbL self-assembly has promise as a simple approach for the wound dressing for interventional therapeutics in the clinic, as illustrated by the successful in situ fabrication of drug-free layers consisting of chitosan and heparin on the dorsal skin of diabetic mice to rescue defective wound healing. This bioinspired self-assembly approach is expected to provide a robust and versatile platform with which to explore the surface engineering of nanofilms in science, engineering, and medicine.

Shear-Induced Unfolding and Enzymatic Cleavage of Full-Length VWF Multimers

Proteolysis of the multimeric blood coagulation protein von Willebrand Factor (VWF) by ADAMTS13 is crucial for prevention of microvascular thrombosis. ADAMTS13 cleaves VWF within the mechanosensitive A2 domain, which is believed to open under shear flow. In this study, we combine fluorescence correlation spectroscopy (FCS) and a microfluidic shear cell to monitor real-time kinetics of full-length VWF proteolysis as a function of shear stress. For comparison, we also measure the Michaelis-Menten kinetics of ADAMTS13 cleavage of wild-type VWF in the absence of shear but partially denaturing conditions. Under shear, ADAMTS13 activity on full-length VWF arises without denaturing agent as evidenced by FCS and gel-based multimer analysis. In agreement with Brownian hydrodynamics simulations, we find a sigmoidal increase of the enzymatic rate as a function of shear at a threshold shear rate γ˙1/2 = 5522/s. The same flow-rate dependence of ADAMTS13 activity we also observe in blood plasma, which is relevant to predict hemostatic dysfunction.

Stretching self-entangled DNA molecules in elongational fields

We present experiments of self-entangled DNA molecules stretching under a planar elongational field, and their stretching dynamics are compared to identical molecules without entanglements. Self-entangled molecules stretch in a stage-wise fashion, persisting in an "arrested" state for decades of strain prior to rapidly stretching, slowing down the stretching dynamics by an order of magnitude compared to unentangled molecules. Self-entangled molecules are shown to proceed through a transient state where one or two ends of the molecule are protruding from an entangled, knotted core. This phenomenon sharply contrasts with the wide array of transient configurations shown here and by others for stretching polymers without entanglements. The rate at which self-entangled molecules stretch through this transient state is demonstrably slower than unentangled molecules, providing the first direct experimental evidence of a topological friction. These experimental observations are shown to be qualitatively and semi-quantitatively reproduced by a dumbbell model with two fitting parameters, the values of which are reasonable in light of previous experiments of knotted DNA.

Fluid Mechanics of Blood Clot Formation

Intravascular blood clots form in an environment in which hydrodynamic forces dominate and in which fluid-mediated transport is the primary means of moving material. The clotting system has evolved to exploit fluid dynamic mechanisms and to overcome fluid dynamic challenges to ensure that clots that preserve vascular integrity can form over the wide range of flow conditions found in the circulation. Fluid-mediated interactions between the many large deformable red blood cells and the few small rigid platelets lead to high platelet concentrations near vessel walls where platelets contribute to clotting. Receptor-ligand pairs with diverse kinetic and mechanical characteristics work synergistically to arrest rapidly flowing cells on an injured vessel. Variations in hydrodynamic stresses switch on and off the function of key clotting polymers. Protein transport to, from, and within a developing clot determines whether and how fast it grows. We review ongoing experimental and modeling research to understand these and related phenomena.

Von Willlebrand adhesion to surfaces at high shear rates is controlled by long-lived bonds

Von Willebrand factor (vWF) adsorbs and immobilizes platelets at sites of injury under high-shear-rate conditions. It has been recently demonstrated that single vWF molecules only adsorb significantly to collagen above a threshold shear, and here we explain such counterintuitive behavior using a coarse-grained simulation and a phenomenological theory. We find that shear-induced adsorption only occurs if the vWF-surface bonds are slip-resistant such that force-induced unbinding is suppressed, which occurs in many biological bonds (i.e., catch bonds). Our results quantitatively match experimental observations and may be important to understand the activation and mechanical regulation of vWF activity during blood clotting.

Shear-induced dynamics of polymeric globules at adsorbing homogeneous and inhomogeneous surfaces

The dynamics and adsorption behavior of a single collapsed homopolymer on a surface in shear flow is investigated by means of Brownian hydrodynamics simulations. We study different homogeneous and inhomogeneous surface models and determine dynamic state diagrams as a function of the cohesive strength, the adhesive strength, and the shear rate. We find distinct dynamical adsorbed states that are classified into rolling and slipping states, globular and coil-like states, as well as isotropic and prolate states. We identify two different cyclic processes based on trajectories of the polymer stretching and the polymer separation from the surface. For adsorption on an inhomogeneous surface consisting of discrete binding sites, we observe stick-roll motion for highly corrugated surface potentials. Although the resulting high surface friction leads to low drift velocities and reduced hydrodynamic lift forces on such inhomogeneous surfaces, a shear-induced adsorption is not found in the presence of full hydrodynamic interactions. A hydrodynamically stagnant surface model is introduced for which shear-induced adsorption is observed in the absence of hydrodynamic interactions.