No Arabic abstract
We have simulated the motion of a bead subjected to a constant force while embedded in a network of semiflexible polymers which can represent actin filaments. We find that the bead displacement obeys the power law x ~ t^alfa. After the initial stage characterized by the exponent alfa=0.75 we find a new regime with alfa=0.5. The response in this regime is linear in force and scales with the polymer concentration as c^(-1.4). We find that the polymers pile up ahead of the moving bead, while behind it the polymer density is reduced. We show that the force resisting the bead motion is due to steric repulsion exerted by the polymers on the front hemisphere of the bead.
Using a recently developed bead-spring model for semiflexible polymers that takes into account their natural extensibility, we report an efficient algorithm to simulate the dynamics for polymers like double-stranded DNA (dsDNA) in the absence of hydrodynamic interactions. The dsDNA is modelled with one bead-spring element per basepair, and the polymer dynamics is described by the Langevin equation. The key to efficiency is that we describe the equations of motion for the polymer in terms of the amplitudes of the polymers fluctuation modes, as opposed to the use of the physical positions of the beads. We show that, within an accuracy tolerance level of $5%$ of several key observables, the model allows for single Langevin time steps of $approx1.6$, 8, 16 and 16 ps for a dsDNA model-chain consisting of 64, 128, 256 and 512 basepairs (i.e., chains of 0.55, 1.11, 2.24 and 4.48 persistence lengths) respectively. Correspondingly, in one hour, a standard desktop computer can simulate 0.23, 0.56, 0.56 and 0.26 ms of these dsDNA chains respectively. We compare our results to those obtained from other methods, in particular, the (inextensible discretised) WLC model. Importantly, we demonstrate that at the same level of discretisation, i.e., when each discretisation element is one basepair long, our algorithm gains about 5-6 orders of magnitude in the size of time steps over the inextensible WLC model. Further, we show that our model can be mapped one-on-one to a discretised version of the extensible WLC model; implying that the speed-up we achieve in our model must hold equally well for the latter. We also demonstrate the use of the method by simulating efficiently the tumbling behaviour of a dsDNA segment in a shear flow.
Mucus is a viscoelastic gel secreted by the pulmonary epithelium in the tracheobronchial region of the lungs. The coordinated beating of cilia moves mucus upwards towards pharynx, removing inhaled pathogens and particles from the airways. The efficacy of this clearance mechanism depends primarily on the rheological properties of mucus. Here we use magnetic wire based microrheology to study the viscoelastic properties of human mucus collected from human bronchus tubes. The response of wires between 5 and 80 microns in length to a rotating magnetic field is monitored by optical time-lapse microscopy and analyzed using constitutive equations of rheology, including those of Maxwell and Kelvin-Voigt. The static shear viscosity and elastic modulus can be inferred from low frequency (from 0.003 to 30 rad s-1) measurements, leading to the evaluation of the mucin network relaxation time. This relaxation time is found to be widely distributed, from one to several hundred seconds. Mucus is identified as a viscoelastic liquid with an elastic modulus of 2.5 +/- 0.5 Pa and a static viscosity of 100 +/- 40 Pa s. Our work shows that beyond the established spatial variations in rheological properties due to microcavities, mucus exhibits secondary inhomogeneities associated with the relaxation time of the mucin network that may be important for its flow properties.
We present a theoretical framework for the linear and nonlinear visco-elastic properties of reversibly crosslinked networks of semiflexible polymers. In contrast to affine models where network strain couples to the polymer end-to-end distance, in our model strain rather serves to locally distort the network structure. This induces bending modes in the polymer filaments, the properties of wich are slaved to the surrounding network structure. Specifically, we investigate the frequency-dependent linear rheology, in particular in combination with crosslink binding/unbinding processes. We also develop schematic extensions to describe the nonlinear response during creep measurements as well as during constant-strainrate ramps.
Motivated by the structure of networks of cross-linked cytoskeletal biopolymers, we study the orientationally ordered phases in two-dimensional networks of randomly cross-linked semiflexible polymers. We consider permanent cross-links which prescribe a finite angle and treat them as quenched disorder in a semi-microscopic replica field theory. Starting from a fluid of un-cross-linked polymers and small polymer clusters (sol) and increasing the cross-link density, a continuous gelation transition occurs. In the resulting gel, the semiflexible chains either display long range orientational order or are frozen in random directions depending on the value of the crossing angle, the crosslink concentration and the stiffness of the polymers. A crossing angle $thetasim 2pi/M$ leads to long range $M$-fold orientational order, e.g., hexatic or tetratic for $theta=60^{circ}$ or $90^{circ}$, respectively. The transition is discontinuous and the critical cross-link density depends on the bending stiffness of the polymers and the cross-link geometry: the higher the stiffness and the lower $M$, the lower the critical number of cross-links. In between the sol and the long range ordered state, we always observe a gel which is a statistically isotropic amorphous solid (SIAS) with random positional and random orientational localization of the participating polymers.
In this work we study the assisted translocation of a polymer across a membrane nanopore, inside which a molecular motor exerts a force fuelled by the hydrolysis of ATP molecules. In our model the motor switches to its active state for a fixed amount of time, while it waits for an ATP molecule binding and triggering the impulse, during an exponentially distributed time lapse. The polymer is modelled as a beads-springs chain with both excluded volume and bending contributions, and moves in a stochastic three dimensional environment modelled with a Langevin dynamics at fixed temperature. The resulting dynamics shows a Michaelis-Menten translocation velocity that depends on the chain flexibility. The scaling behavior of the mean translocation time with the polymer length for different bending values is also investigated.