Tracer diffusion in a crowded cylindrical channel

Structure and dynamics of an active polymer chain inside a nanochannel grafted with polymers

We use computer simulations to investigate the complex dynamics of a polymer, made of active Brownian particles, inside a channel grafted internally with passive polymer chains. Our simulations reveal that this probe-polymer, if passive, exhibits a compact structure when its interaction is repulsive with the grafted chains as it tends to stay within the hollow space created along the axis of the channel. On increasing the attractive interaction, the passive probe-polymer is pulled towards the grafted polymeric region and adopts an extended structure. By contrast, switching on the activity helps the probe-polymer to escape from the local traps caused by the sticky grafted chains. The interplay between the activity of the probe-polymer and its sticky interaction with the grafted chains results in shrinking, followed by swelling as the activity is increased. To elucidate the dynamics we compute the mean square displacement (MSD) of the center of mass of the probe-polymer, which increases monotonically with activity and displays superdiffusive behavior at an intermediate time and enhanced diffusion at a long time period. In addition, compared with the attractive interaction, the active probe-polymer shows faster dynamics when the interaction is repulsive to the grafted polymers. We believe that our current study will provide insights into the structural changes and dynamics of active polymers in heterogeneous media and will be useful in designing polymer-based drug delivery vehicles.

Specific Ion and Electric Field Controlled Diverse Ion Distribution and Electroosmotic Transport in a Polyelectrolyte Brush Grafted Nanochannel

Controlling ion distribution inside a charged nanochannel is central to using such channels in diverse applications. Here, we show the possibility of using a charged polyelectrolyte (PE) brush-grafted nanochannel for triggering diverse nanoscopic ion distribution and nanofluidic electroosmotic transport by controlling the valence and size of the counterions (that screen the charges of the PE brushes) and the strength of an externally applied axial electric field. We atomistically simulate separate cases of fully charged polyacrylic acid (PAA) brush functionalized nanochannels with Na⁺, Cs⁺, Ca²⁺, Ba²⁺, and Y³⁺ counterions screening the PE charges. Four key findings emerge from our simulations. First, we find that the counterions with a greater valence and a smaller size prefer to remain localized inside the brush layer. Second, for the case where there is an added chloride salt with the same cation (as the screening counterions), there are more coions (Cl^- ions) in the brush-free bulk than counterions (for counterions Na⁺, Ca²⁺, Ba²⁺, Y³⁺): this is a manifestation of the overscreening (OS) of the PE brush layer. Contrastingly, the number of Cs⁺ ions remain higher than the Cl^- ions inside the brush-free bulk, ensuring that there is no OS effect for this case. Third, large applied electric field enables a few Na⁺, Cs⁺, and Ba²⁺ counterions to leave the brush layer and to go to the bulk: this makes the OS of the PE brush layer disappear for the cases of PE brushes being screened by the Na⁺ and Ba²⁺ ions. On the other hand, no such electric-field-mediated disappearance of OS is observed for the cases of Ca²⁺ and Y³⁺ screening counterions; we attribute this to the firm attachment of these counterions to the negatively charged monomers. Free energy associated with a counterion binding to a PE chain corroborates this diversity in the counterion-specific response to the applied electric field. Finally, we demonstrate that such diverse ion distributions, along with specific electric-field-strength-dependent ion properties, lead to (1) electroosmotic (EOS) transport in nanochannels grafted with PAA brushes screened with Cs⁺ ions to be always counterion dominated, (2) EOS transport in nanochannels grafted with PAA brushes screened with Ca²⁺ and Y³⁺ ions to be always coion-dominated, and (3) EOS transport in nanochannels grafted with PAA brushes screened with Na⁺ and Ba²⁺ ions to be coion dominated for smaller electric fields and counterion dominated for larger electric fields.

Transport of a self-propelled tracer through a hairy cylindrical channel: interplay of stickiness and activity

Active transport of biomolecules assisted by motor proteins is imperative for the proper functioning of cellular activities. Inspired by the diffusion of active agents in crowded cellular channels, we computationally investigate the transport of an active tracer through a polymer grafted cylindrical channel by varying the activity of the tracer and stickiness of the tracer to the polymers. Our results reveal that the passive tracer exhibits profound subdiffusion with increasing stickiness by exploring deep into the grafted polymeric zone, while purely repulsive one prefers to diffuse through the pore-like space created along the cylindrical axis of the channel. In contrast, the active tracer shows faster dynamics and intermediate superdiffusion even though the tracer preferentially stays close to the dense polymeric region. This observation is further supported by the sharp peaks in the density profile of the probability of radial displacement of the tracer. We discover that the activity plays an important role in deciding the pathway that the tracer takes through the narrow channel. Interestingly, increasing the activity washes out the effect of stickiness. Adding to this, van-Hove functions manifest that the active tracer dynamics deviates from Gaussianity, and the degree of deviation grows with the activity. Our work has direct implications on how effective transportation and delivery of cargo can be achieved through a confined medium where activity, interactions, and crowding are interplaying. Looking ahead, these factors will be crucial for understanding the mechanism of artificial self-powered machines navigating through the cellular channels and performing in vivo challenging tasks.

Translational and rotational dynamics of a self-propelled Janus probe in crowded environments

We computationally investigate the dynamics of a self-propelled Janus probe in crowded environments. The crowding is caused by the presence of viscoelastic polymers or non-viscoelastic disconnected monomers. Our simulations show that the translational as well as rotational mean square displacements have a distinctive three-step growth for fixed values of self-propulsion force, and steadily increase with self-propulsion, irrespective of the nature of the crowder. On the other hand, in the absence of crowders, the rotational dynamics of the Janus probe is independent of self-propulsion force. On replacing the repulsive polymers with sticky ones, translational and rotational mean square displacements of the Janus probe show a sharp drop. Since different faces of a Janus particle interact differently with the environment, we show that the direction of self-propulsion also affects its dynamics. The ratio of long-time translational and rotational diffusivities of the self-propelled probe with a fixed self-propulsion, when plotted against the area fraction of the crowders, passes through a minimum and at higher area fraction merges to its value in the absence of the crowder. This points towards the decoupling of the translational and rotational dynamics of the self-propelled probe at an intermediate area fraction of the crowders. However, such translational-rotational decoupling is absent for passive probes.

Tracer diffusion in crowded solutions of sticky polymers

Macromolecular diffusion in strongly confined geometries and crowded environments is still to a large extent an open subject in soft matter physics and biology. In this paper, we employ large-scale Langevin dynamics simulations to investigate how the diffusion of a tracer is influenced by the combined action of excluded-volume and weak attractive crowder-tracer interactions. We consider two species of tracers, standard hard-core particles described by the Weeks-Chandler-Andersen (WCA) repulsive potential and core-softened (CS) particles, which model, e.g., globular proteins, charged colloids, and nanoparticles covered by polymeric brushes. These systems are characterized by the presence of two length scales in the interaction and can show waterlike anomalies in their diffusion, stemming from the inherent competition between different length scales. Here we report a comprehensive study of both diffusion and structure of these two tracer species in an environment crowded by quenched configurations of polymers at increasing density. We analyze in detail how the tracer-polymer affinity and the system density affect transport as compared to the emergence of specific static spatial correlations. In particular, we find that, while hardly any differences emerge in the diffusion properties of WCA and CS particles, the propensity to develop structural order for large crowding is strongly frustrated for CS particles. Surprisingly, for large enough affinity for the crowding matrix, the diffusion coefficient of WCA tracers display a nonmonotonic trend as their density is increased when compared to the zero affinity scenario. This waterlike anomaly turns out to be even larger than what observed for CS particle and appears to be rooted in a similar competition between excluded-volume and affinity effects.

Field-driven tracer diffusion through curved bottlenecks: fine structure of first passage events

Using scaling arguments and extensive numerical simulations, we study the dynamics of a tracer particle in a corrugated channel represented by a periodic sequence of broad chambers and narrow funnel-like bottlenecks enclosed by a hard-wall boundary.

Transport of probe particles in a polymer network: effects of probe size, network rigidity and probe-polymer interaction

Fundamental understanding of the effect of microscopic parameters on the dynamics of probe particles in different complex environments has wide implications. Examples include diffusion of proteins in biological hydrogels, porous media, polymer matrix, etc. Here, we use extensive molecular dynamics simulations to investigate the dynamics of the probe particle in a polymer network on a diamond lattice, which provides substantial crowding to mimic the cellular environment. Our simulations show that the dynamics of the probe increasingly becomes restricted, non-Gaussian and subdiffusive on increasing the network rigidity, binding affinity and probe size. In addition, the velocity autocorrelation functions show negative dips owing to the viscoelasticity and caging due to the surrounding network. These observations go with the general experimental findings. Importantly, for a probe particle of size comparable to the mesh size, unrestricted motion engulfing large length scales has been observed. This happens with a more flexible polymer network, which is easily pushed by the bigger probe. On increasing the rigidity of the network, the bigger probe can not efficiently push the network and as a result the long tail disappears. Our study gives a general qualitative picture of the transport of probes in a gel-like medium, as encountered in different contexts.

Polymer translocation through a hairy channel mimicking the inner plug of a nuclear pore complex

A microscopic transport model of a polymer translocating through a nuclear pore complex (NPC) is presented based on self-consistent field theory (SCFT), with the NPC and its nucleoporins mimicked by a hairy channel. Multiple cell environment effects (electrolyte effect, excluded volume effect, NPC drag effect, and hydrophobic effect) are all considered in this hairy channel model. The influence of various parameters (polymer chain length, length of NPC, strength of hydrophobic effect, and excluded volume effect) on translocation time is studied through theoretical analysis and numerical calculation. Numerical simulation results show that an area of low nucleoporin number density exists in the NPC, which facilitates the translocation of the polymer. The results also show that the translocation time curves with increasing NPC length and polymer charge number are concave. In addition, there are critical values for NPC length and polymer charge number for which the translocation time has a minimal value. The translocation time decreases with the increasing strength of the hydrophobic effect and excluded volume effect.

Looking "Under the Hood" of Cellular Mechanotransduction with Computational Tools: A Systems Biomechanics Approach across Multiple Scales

Signal modulation has been developed in living cells throughout evolution to promote utilizing the same machinery for multiple cellular functions. Chemical and mechanical modules of signal transmission and transduction are interconnected and necessary for organ development and growth. However, due to the high complexity of the intercommunication of physical intracellular connections with biochemical pathways, there are many missing details in our overall understanding of mechanotransduction processes, i.e., the process by which mechanical signals are converted to biochemical cascades. Cell-matrix adhesions are mechanically coupled to the nucleus through the cytoskeleton. This modulated and tightly integrated network mediates the transmission of mechanochemical signals from the extracellular matrix to the nucleus. Various experimental and computational techniques have been utilized to understand the basic mechanisms of mechanotransduction, yet many aspects have remained elusive. Recently, in silico experiments have made important contributions to the field of mechanobiology. Herein, computational modeling efforts devoted to understanding integrin-mediated mechanotransduction pathways are reviewed, and an outlook is presented for future directions toward using suitable computational approaches and developing novel techniques for addressing important questions in the field of mechanotransduction.

Diffusing diffusivity: Rotational diffusion in two and three dimensions

We consider the problem of calculating the probability distribution function (pdf) of angular displacement for rotational diffusion in a crowded, rearranging medium. We use the diffusing diffusivity model and following our previous work on translational diffusion [R. Jain and K. L. Sebastian, J. Phys. Chem. B 120, 3988 (2016)], we show that the problem can be reduced to that of calculating the survival probability of a particle undergoing Brownian motion, in the presence of a sink. We use the approach to calculate the pdf for the rotational motion in two and three dimensions. We also propose new dimensionless, time dependent parameters, αrot,2D and αrot,3D, which can be used to analyze the experimental/simulation data to find the extent of deviation from the normal behavior, i.e., constant diffusivity, and obtain explicit analytical expressions for them, within our model.

Tracer diffusion in a sea of polymers with binding zones: mobile vs. frozen traps

We use molecular dynamics simulations to investigate the tracer diffusion in a sea of polymers with specific binding zones for the tracer. These binding zones act as traps. Our simulations show that the tracer can undergo normal yet non-Gaussian diffusion under certain circumstances, e.g., when the polymers with traps are frozen in space and the volume fraction and the binding strength of the traps are moderate. In this case, as the tracer moves, it experiences a heterogeneous environment and exhibits confined continuous time random walk (CTRW) like motion resulting in a non-Gaussian behavior. Also the long time dynamics becomes subdiffusive as the number or the binding strength of the traps increases. However, if the polymers are mobile then the tracer dynamics is Gaussian but could be normal or subdiffusive depending on the number and the binding strength of the traps. In addition, with increasing binding strength and number of polymer traps, the probability of the tracer being trapped increases. On the other hand, removing the binding zones does not result in trapping, even at comparatively high crowding. Our simulations also show that the trapping probability increases with the increasing size of the tracer and for a bigger tracer with the frozen polymer background the dynamics is only weakly non-Gaussian but highly subdiffusive. Our observations are in the same spirit as found in many recent experiments on tracer diffusion in polymeric materials and question the validity of using Gaussian theory to describe diffusion in a crowded environment in general.

Diffusion in a Crowded, Rearranging Environment

It has been found in many experiments that the mean square displacement of a Brownian particle x(T) diffusing in a rearranging environment is strictly Fickian, obeying ⟨(x(T))(2)⟩ ∝ T, but the probability distribution function for the displacement is not Gaussian. An explanation of this is that the diffusivity of the particle itself is changing as a function of time. Models for this diffusing diffusivity have been solved analytically in the limit of small time, but simulations were necessary for intermediate and large times. We show that one of the diffusing diffusivity models is equivalent to Brownian motion in the presence of a sink and introduce a class of models for which it is possible to find analytical solutions. Our solution gives ⟨(x(T))(2)⟩ ∝ T for all times and at short times the probability distribution function of the displacement is exponential which crosses over to a Gaussian in the limit of long times and large displacements.

Non-universal tracer diffusion in crowded media of non-inert obstacles

For tracer motion in an array of attractive obstacles we observe transient, non-ergodic anomalous diffusion depending on the obstacle density.