Instytut Podstawowych Problemów Techniki

Streszczenie:

We present a numerical study of the dynamics of an elastic fibre in a shear flow at low Reynolds number, and seek to understand several aspects of the fibre's motion using the equations for slender-body theory coupled to the elastica. The numerical simulations are performed in the bead-spring framework including hydrodynamic interactions in two theoretical schemes: the generalized Rotne-Prager-Yamakawa model and a multipole expansion corrected for lubrication forces. In general, the two schemes yield similar results, including for the dominant scaling features of the shape that we identify. In particular, we focus on the evolution of an initially straight fibre oriented in the flow direction and show that the time scales of fibre bending, curling and rotation, which depend on its length and stiffness, determine the overall motion and evolution of the shapes. We document several characteristic time scales and curvatures representative of the shape that vary as power laws of the bending stiffness and fibre length. The numerical results are further supported by an interpretation using an elastica model.

Słowa kluczowe:

Stokesian dynamics, particle/fluid flow, slender-body theory

Afiliacje autorów:

Streszczenie:

Inclusion of hydrodynamic interactions is essential for a quantitatively accurate Brownian dynamics simulation of colloidal suspensions or polymer solutions. We use the generalized Rotne–Prager–Yamakawa (GRPY) approximation, which takes into account all long-ranged terms in the hydrodynamic interactions, to derive the complete set of hydrodynamic matrices in different geometries: unbounded space, periodic boundary conditions of Lees–Edwards type, and vicinity of a free surface. The construction is carried out both for non-overlapping as well as for overlapping particles. We include the dipolar degrees of freedom, which allows one to use this formalism to simulate the dynamics of suspensions in a shear flow and to study the evolution of their rheological properties. Finally, we provide an open-source numerical package, which implements the GRPY algorithm in Lees–Edwards periodic boundary conditions.

Afiliacje autorów:

Streszczenie:

Two oppositely charged surfaces separated by a dielectric medium attract each other. In contrast we observe a strong repulsion between two plates of a capacitor that is filled with an aqueous electrolyte upon application of an alternating potential difference between the plates. This long-range force increases with the ratio of diffusion coefficients of the ions in the medium and reaches a steady state after a few minutes, which is much larger than the millisecond timescale of diffusion across the narrow gap. The repulsive force, an order of magnitude stronger than the electrostatic attraction observed in the same setup in air, results from the increase in osmotic pressure as a consequence of the field-induced excess of cations and anions due to lateral transport from adjacent reservoirs.

Afiliacje autorów:

Streszczenie:

The charging of electrical double layers inside a cylindrical pore has applications to supercapacitors, batteries, desalination and biosensors. The charging dynamics in the limit of thin double layers, i.e., when the double layer thickness is much smaller than the pore radius, is commonly described using an effective RC transmission line circuit. Here, we perform direct numerical simulations (DNS) of the Poisson-Nernst-Planck equations to study the double layer charging for the scenario of overlapping double layers, i.e., when the double layer thickness is comparable to the pore radius. We develop an analytical model that accurately predicts the results of DNS. Also, we construct a modified effective circuit for the overlapping double layer limit, and find that the modified circuit is identical to the RC transmission line but with different values and physical interpretation of the capacitive and resistive elements. In particular, the effective surface potential is reduced, the capacitor represents a volumetric current source, and the charging timescale is weakly dependent on the ratio of the pore radius and the double layer thickness.

Afiliacje autorów:

Streszczenie:

Injection of cell‐laden scaffolds in the form of mesoscopic particles directly to the site of treatment is one of the most promising approaches to tissue regeneration. Here, a novel and highly efficient method is presented for preparation of porous microbeads of tailorable dimensions (in the range ≈300–1500 mm) and with a uniform and fully interconnected internal porous texture. The method starts with generation of a monodisperse oil‐in‐water emulsion inside a flow‐focusing microfluidic device. This emulsion is later broken‐up, with the use of electric field, into mesoscopic double droplets, that in turn serve as a template for the porous microbeads. By tuning the amplitude and frequency of the electric pulses, the template droplets and the resulting porous bead scaffolds are precisely produced. Furthermore, a model of pulsed electrodripping is proposed that predicts the size of the template droplets as a function of the applied voltage. To prove the potential of the porous microbeads as cell carries, they are tested with human mesenchymal stem cells and hepatic cells, with their viability and degree of microbead colonization being monitored. Finally, the presented porous microbeads are benchmarked against conventional microparticles with nonhomogenous internal texture, revealing their superior performance.

Afiliacje autorów:

Streszczenie:

Two main problems that arise in the context of hydrodynamic bead modeling are an inaccurate treatment of bead overlaps and the necessity of using volume corrections when calculating intrinsic viscosity. We present a formalism based on the generalized Rotne-Prager-Yamakawa approximation that successfully addresses both of these issues. The generalized Rotne-Prager-Yamakawa method is shown to be highly effective for the calculation of transport properties of rigid biomolecules represented as assemblies of spherical beads of different sizes, both overlapping and nonoverlapping. We test the method on simple molecular shapes as well as real protein structures and compare its performance with other computational approaches.

Afiliacje autorów:

Streszczenie:

In bistable reaction–diffusion systems, transitions between stable states typically occur on timescales orders of magnitude longer than the chemical equilibration time. Estimation of transition rates within explicit Brownian dynamics simulations is computationally prohibitively costly. We present a method that exploits a single trajectory, generated by a prior simulation of diffusive motions of molecules, to sample chemical kinetic processes on timescales several orders of magnitude longer than the duration of the diffusive trajectory. In this approach, we „loop” the diffusive trajectory by transferring chemical states of the molecules from the last to the first time step of the trajectory. Trajectory looping can be applied to enhance sampling of rare events in biochemical systems in which the number of reacting molecules is constant, as in cellular signal transduction pathways. As an example, we consider a bistable system of autophosphorylating kinases, for which we calculate state-to-state transition rates and traveling wave velocities. We provide an open-source implementation of the method.

Afiliacje autorów:

Streszczenie:

We develop a generalized Rotne–Prager–Yamakawa approximation for the dipolar components of the inverse friction matrix and use it for calculating the intrinsic viscosity of rigidly connected bead conglomerates. Such bead models are commonly used in the calculation of hydrodynamic properties of macromolecules. We consider both the case of non-overlapping constituent beads as well as overlapping beads of different sizes. We demonstrate the accuracy of the approximation in two test cases and show that it performs well even if the distances between the beads are small or if the beads overlap. Robust performance of this approximation in the case of overlapping beads stems from its correct limiting behaviour at a complete overlap, with one sphere fully immersed in the other. The generalized Rotne–Prager–Yamakawa approximation is thus well suited for evaluation of intrinsic viscosity, which is a key quantity in characterizing molecular conformations of biological macromolecules.

Słowa kluczowe:

complex fluids, low-Reynolds-number flows, mathematical foundations

Afiliacje autorów:

Streszczenie:

The Rotne–Prager–Yamakawa (RPY) approximation is a commonly used approach to model the hydrodynamic interactions between small spherical particles suspended in a viscous fluid at a low Reynolds number. However, when the particles overlap, the RPY tensors lose their positive definiteness, which leads to numerical problems in the Brownian dynamics simulations as well as errors in calculations of the hydrodynamic properties of rigid macromolecules using bead modelling. These problems can be avoided by using regularizing corrections to the RPY tensors; so far, however, these corrections have only been derived for equal-sized particles. Here we show how to generalize the RPY approach to the case of overlapping spherical particles of different radii and present the complete set of mobility matrices for such a system. In contrast to previous ad hoc approaches, our method relies on the direct integration of force densities over the sphere surfaces and thus automatically provides the correct limiting behaviour of the mobilities for the touching spheres and for a complete overlap, with one sphere immersed in the other one. This approach can then be used to calculate hydrodynamic properties of complex macromolecules using bead models with overlapping, different-sized beads, which we illustrate with an example.

Słowa kluczowe:

complex fluids, computational methods, low-Reynolds-number flows, mathematical foundations, suspensions

Afiliacje autorów:

Streszczenie:

Rotne-Prager-Yamakawa approximation is a commonly used approach to model hydrodynamic interactions between particles suspended in fluid. It takes into account all the long-range contributions to the hydrodynamic tensors, with the corrections decaying at least as fast as the inverse fourth power of the interparticle distances, and results in a positive definite mobility matrix, which is fundamental in Brownian dynamics simulations. In this communication, we show how to construct the Rotne-Prager-Yamakawa approximation for the bulk system under shear flow, which is modeled using the Lees-Edwards boundary conditions.

Słowa kluczowe:

Tensor methods, Hydrodynamics, Shear flows, Brownian dynamics, Boundary value problems

Afiliacje autorów:

Streszczenie:

NF-κB is a key transcription factor that regulates innate immune response. Its activity is tightly controlled by numerous feedback loops, including two negative loops mediated by NF-κB inducible inhibitors, IκBα and A20, which assure oscillatory responses, and by positive feedback loops arising due to the paracrine and autocrine regulation via TNFα, IL-1 and other cytokines. We study the NF-κB system of interlinked negative and positive feedback loops, combining bifurcation analysis of the deterministic approximation with stochastic numerical modeling. Positive feedback assures the existence of limit cycle oscillations in unstimulated wild-type cells and introduces bistability in A20-deficient cells. We demonstrated that cells of significant autocrine potential, i.e., cells characterized by high secretion of TNFα and its receptor TNFR1, may exhibit sustained cytoplasmic–nuclear NF-κB oscillations which start spontaneously due to stochastic fluctuations. In A20-deficient cells even a small TNFα expression rate qualitatively influences system kinetics, leading to long-lasting NF-κB activation in response to a short-pulsed TNFα stimulation. As a consequence, cells with impaired A20 expression or increased TNFα secretion rate are expected to have elevated NF-κB activity even in the absence of stimulation. This may lead to chronic inflammation and promote cancer due to the persistent activation of antiapoptotic genes induced by NF-κB. There is growing evidence that A20 mutations correlate with several types of lymphomas and elevated TNFα secretion is characteristic of many cancers. Interestingly, A20 loss or dysfunction also leaves the organism vulnerable to septic shock and massive apoptosis triggered by the uncontrolled TNFα secretion, which at high levels overcomes the antiapoptotic action of NF-κB. It is thus tempting to speculate that some cancers of deregulated NF-κB signaling may be prone to the pathogen-induced apoptosis.

Afiliacje autorów:

Streszczenie:

The Rotne–Prager–Yamakawa approximation is one of the most commonly used methods of including hydrodynamic interactions in modelling of colloidal suspensions and polymer solutions. The two main merits of this approximation are that it includes all long-range terms (i.e. decaying as R−3 or slower in interparticle distances) and that the diffusion matrix is positive definite, which is essential for Brownian dynamics modelling. Here, we extend the Rotne–Prager–Yamakawa approach to include both translational and rotational degrees of freedom, and derive the regularizing corrections to account for overlapping particles. Additionally, we show how the Rotne–Prager–Yamakawa approximation can be generalized for other geometries and boundary conditions.

Słowa kluczowe:

computational methods, low-Reynolds-number flows, suspensions

Afiliacje autorów:

Streszczenie:

The aim of this study is to demonstrate that in molecular dynamical systems with the underlying bi- or multistability, the type of noise determines the most strongly attracting steady state or stochastic attractor. As an example we consider a simple stochastic model of autoregulatory gene with a nonlinear positive feedback, which in the deterministic approximation has two stable steady state solutions. Three types of noise are considered: transcriptional and translational – due to the small number of gene product molecules and the gene switching noise – due to gene activation and inactivation transitions. We demonstrate that the type of noise in addition to the noise magnitude dictates the allocation of probability mass between the two stable steady states. In particular, we found that when the gene switching noise dominates over the transcriptional and translational noise (which is characteristic of eukaryotes), the gene preferentially activates, while in the opposite case, when the transcriptional noise dominates (which is characteristic of prokaryotes) the gene preferentially remains inactive. Moreover, even in the zero-noise limit, when the probability mass generically concentrates in the vicinity of one of two steady states, the choice of the most strongly attracting steady state is noise type-dependent. Although the epigenetic attractors are defined with the aid of the deterministic approximation of the stochastic regulatory process, their relative attractivity is controlled by the type of noise, in addition to noise magnitude. Since noise characteristics vary during the cell cycle and development, such mode of regulation can be potentially employed by cells to switch between alternative epigenetic attractors.

Słowa kluczowe:

Gene expression, Bistability, Stochastic processes, Epigenetic attractors

Afiliacje autorów:

Streszczenie:

Living cells may be considered as biochemical reactors of multiple steady states. Transitions between these states are enabled by noise, or, in spatially extended systems, may occur due to the traveling wave propagation. We analyze a one-dimensional bistable stochastic birth–death process by means of potential and temperature fields. The potential is defined by the deterministic limit of the process, while the temperature field is governed by noise. The stable steady state in which the potential has its global minimum defines the global deterministic attractor. For the stochastic system, in the low noise limit, the stationary probability distribution becomes unimodal, concentrated in one of two stable steady states, defined in this study as the global stochastic attractor. Interestingly, these two attractors may be located in different steady states. This observation suggests that the asymptotic behavior of spatially extended stochastic systems depends on the substrate diffusivity and size of the reactor. We confirmed this hypothesis within kinetic Monte Carlo simulations of a bistable reaction–diffusion model on the hexagonal lattice. In particular, we found that although the kinase–phosphatase system remains inactive in a small domain, the activatory traveling wave may propagate when a larger domain is considered.

Słowa kluczowe:

multi-stability, Markov process, spatially extended system, kinetic Monte Carlo on the lattice, cell signalling

Afiliacje autorów:

Streszczenie:

The aim of this study is to demonstrate that in dynamical systems with underlying bistability the type of noise qualitatively influences the stationary probability distribution (SPD). Specifically, we consider a simplified model of gene expression with the nonlinear positive feedback, which in the deterministic approximation has two stable steady state solutions. Two types of noise are considered; transcriptional - due to the limited number of protein molecules, and gene switching noise - due to gene activation and inactivation. In the limit of zero noise, the SPD generically concentrates in the decreasing vicinity of one of the two stable steady states. We demonstrated that for a range of parameters the SPD corresponding to the system with transcriptional noise only concentrates around a different steady state than SPD corresponding to the system with gene switching noise only.

Słowa kluczowe:

Gene expression, Bistability, Stochastic processes, Epigenetic attractors

Afiliacje autorów:

Streszczenie:

We consider a kinase auto-activation model in which the number of activated kinases follows the timecontinuous Markov process. In the deterministic approximation the process is described by the single nonlinear ordinary differential equation, which may have two stable steady states. We found that for sufficiently large number of kinases, the stationary probability distribution given by the Markov process concentrates in the vicinity of the two stable steady states of the deterministic approximation. However, if the number of kinases diverges to the infinity (zero noise limit), the stationary probability distribution concentrates (generically) in only one of the two steady states.

Afiliacje autorów: