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 γ̇.

Replica-averaging: An algorithm to study mechano-reactive processes for polymers under flow conditions.

A new method based on quasi-independent parallel simulation approach, replica-averaging, has been developed to study the influence of flow on mechanical force-mediated polymer processes, such as denaturation and breaking of bonds. This method considerably mitigates the unphysical prediction of force-mediated events inherent in Brownian dynamics (BD) polymer chain simulations that employ instantaneous force profile-based criteria to identify the occurrence of such events. This inaccuracy in predicting force-mediated event kinetics is due to high fluctuations of the instantaneous force profile around the average force. Replica-averaging reduces such high fluctuation effects by computing a force profile that faithfully represents the average force profile of the polymer chain conformation, which is then used to predict reactive events. For transient conformation conditions, the replica-averaged method more accurately predicts the mechano-reactive kinetics than the time-averaged method typically employed to reduce the unphysical prediction of force-mediated events in BD simulations. Furthermore, the influence of the proposed replica-averaging method parameters on the accuracy of predicting the true average force profile along the polymer is discussed.

Asymptotic method for entropic multiple relaxation time model in lattice Boltzmann method.

To improve the numerical stability of the lattice Boltzmann method, Karlin et al. [Phys. Rev. E 90, 031302(R) (2014)10.1103/PhysRevE.90.031302] proposed the entropic multiple relaxation time (EMRT) collision model. The idea behind EMRT is to construct an optimal postcollision state by maximizing its local entropy value. The critical step of the EMRT model is to solve the entropy maximization problem under certain constraints, which is often computationally expensive and even not feasible. In this paper, we propose to employ perturbation theory and obtain an asymptotic solution to the maximum entropy state. With mathematical analysis of particular cases under relaxed constraints, we obtain the unperturbed form of the original problem and derive the asymptotic solution. We show that the asymptotic solution well approximates the optimal states; thus, our approach provides an efficient way to solve the constrained maximum entropy problem in the EMRT model. Also, we use the same idea of the EMRT model for the initial condition of the distribution function and propose to leave the entropy function to determine the missing information at the initial nodes. Finally, we numerically verify that the simulation results of the EMRT model obtained via the perturbation theory agree well with the exact solution to the Taylor-Green vortex problem. Furthermore, we also demonstrate that the EMRT model exhibits excellent stability performance for under-resolved simulations in the doubly periodic shear layer flow problem.

Label-free focusing of viral particles under a temperature gradient coupled with continuous swirling flow.

The advances of biomedicine and biotechnology demand new approaches to enrich biological nanoparticles, such as viruses, viral vectors and nanovesicles, in an easy-to-operate fashion. Conventional methods, such as ultracentrifugation and ultrafiltration, require bulky instruments and extensive manual operation. Inspired by recent research of thermophoresis of biomolecules and bio-nanoparticles in aqueous solutions, we present a microfluidic design that directly focuses nanoparticles in a label-free and flow-through process by coupling an engineered swirling flow and a moderate, one-dimensional temperature gradient. Enrichment of polystyrene particles, HIV and bacteriophage samples was quantitatively determined, indicating the compatibility of the microfluidic approach with synthetic and biological samples. The focusing results are well predicted using a numerical model. As thermophoresis is ubiquitous, the microfluidic approach can be applied broadly to bio-nanoparticle enrichment without the necessity of labeling, buffer exchange, or sheath fluids, permitting continuous retrieval of concentrated species in a simple, controlled flow with little infrastructure needs.

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.

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.

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.

Prediction of Sub-Monomer A2 Domain Dynamics of the von Willebrand Factor by Machine Learning Algorithm and Coarse-Grained Molecular Dynamics Simulation.

We develop a machine learning tool useful for predicting the instantaneous dynamical state of sub-monomer features within long linear polymer chains, as well as extracting the dominant macromolecular motions associated with sub-monomer behaviors of interest. We employ the tool to better understand and predict sub-monomer A2 domain unfolding dynamics occurring amidst the dominant large-scale macromolecular motions of the biopolymer von Willebrand Factor (vWF) immersed in flow. Results of coarse-grained Molecular Dynamics (MD) simulations of non-grafted vWF multimers subject to a shearing flow were used as input variables to a Random Forest Algorithm (RFA). Twenty unique features characterizing macromolecular conformation information of vWF multimers were used for training the RFA. The corresponding responses classify instantaneous A2 domain state as either folded or unfolded, and were directly taken from coarse-grained MD simulations. Three separate RFAs were trained using feature/response data of varying resolution, which provided deep insights into the highly correlated macromolecular dynamics occurring in concert with A2 domain unfolding events. The algorithm is used to analyze results of simulation, but has been developed for use with experimental data as well.

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.

Long-ranged Protein-glycan Interactions Stabilize von Willebrand Factor A2 Domain from Mechanical Unfolding.

von Willebrand Factor (vWF) is a large multimeric protein that binds to platelets and collagen in blood clotting. vWF A2 domain hosts a proteolytic site for ADAMTS13 (A Disintegrin and Metalloprotease with a ThromboSpondin type 1 motif, member 13) to regulate the size of vWF multimers. This regulation process is highly sensitive to force conditions and protein-glycan interactions as the process occurs in flowing blood. There are two sites on A2 domain (N1515 and N1574) bearing various N-linked glycan structures. In this study, we used molecular dynamics (MD) simulation to study the force-induced unfolding of A2 domain with and without a single N-linked glycan type on each site. The sequential pullout of ß-strands was used to represent a characteristic unfolding sequence of A2. This unfolding sequence varied due to protein-glycan interactions. The force-extension and total energy-extension profiles also show differences in magnitude but similar characteristic shapes between the systems with and without glycans. Systems with N-linked glycans encountered higher energy barriers for full unfolding and even for unfolding up to the point of ADAMTS13 cleavage site exposure. Interestingly, there is not much difference observed for A2 domain structure itself with and without glycans from standard MD simulations, suggesting roles of N-glycans in A2 unfolding through long-ranged protein-glycan interactions.

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.

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.

Measuring the Soret coefficient of nanoparticles in a dilute suspension.

Thermophoresis is an efficient process for the manipulation of molecules and nanoparticles due to the strong force it generates on the nanoscale. Thermophoresis is characterized by the Soret coefficient. Conventionally, the Soret coefficient of nanosized species is obtained by fitting the concentration profile under a temperature gradient at the steady state to a continuous phase model. However, when the number density of the target is ultralow and the dispersed species cannot be treated as a continuous phase, the bulk concentration fluctuates spatially, preventing extraction of temperature-gradient induced concentration profile. The present work demonstrates a strategy to tackle this problem by superimposing snapshots of nanoparticle distribution. The resulting image is suitable for the extraction of the Soret coefficient through the conventional data fitting method. The strategy is first tested through a discrete phase model that illustrates the spatial fluctuation of the nanoparticle concentration in a dilute suspension in response to the temperature gradient. By superimposing snapshots of the stochastic distribution, a thermophoretic depletion profile with low standard error is constructed, indicative of the Soret coefficient. Next, confocal analysis of nanoparticle distribution in response to a temperature gradient is performed using polystyrene nanobeads down to 1e-5% (v/v). The experimental results also reveal that superimposing enhances the accuracy of extracted Soret coefficient. The critical particle number density in the superimposed image for predicting the Soret coefficient is hypothesized to depend on the spatial resolution of the image. This study also demonstrates that the discrete phase model is an effective tool to study particle migration under thermophoresis in the liquid phase.

Gravity-induced swirl of nanoparticles in microfluidics.

Parallel flows of two fluids in microfluidic devices are used for miniaturized chemistry, physics, biology and bioengineering studies, and the streams are often considered to remain parallel. However, as the two fluids do not always have the same density, interface reorientation induced by density stratification is unavoidable. In this paper, flow characteristics of an aqueous polystyrene nanofluid and a sucrose-densified aqueous solution flowing parallel in microchannels are examined. Nanoparticles 100 nm in diameter are used in the study. The motion of the nanoparticles is simulated using the Lagrangian description and directly observed by a confocal microscope. Matched results are obtained from computational and empirical analysis. Although solution density homogenizes rapidly resulting from a fast diffusion of sucrose in water, the nanofluid is observed to rotate for an extended period. Angular displacement of the nanofluid depends on the ratio of gravitational force to viscous force, Re/Fr2, where Re is the Reynolds number and Fr is the Froude number. In the developing region at the steady state, the angular displacement is related to y/Dh, the ratio between distance from the inlet and the hydraulic diameter of the microfluidic channel. The development of nanofluid flow feature also depends on h/w, the ratio of microfluidic channel's height to width. The quantitative description of the angular displacement of nanofluid will aid rational designs of microfluidic devices utilizing multistream, multiphase flows.