Forced transport of deformable containers through narrow constrictions

Acoustothermal transfection for cell therapy

Transfected stem cells and T cells are promising in personalized cell therapy and immunotherapy against various diseases. However, existing transfection techniques face a fundamental trade-off between transfection efficiency and cell viability; achieving both simultaneously remains a substantial challenge. This study presents an acoustothermal transfection method that leverages acoustic and thermal effects on cells to enhance the permeability of both the cell membrane and nuclear envelope to achieve safe, efficient, and high-throughput transfection of primary T cells and stem cells. With this method, two types of plasmids were simultaneously delivered into the nuclei of mesenchymal stem cells (MSCs) with efficiencies of 89.6 ± 1.2%. CXCR4-transfected MSCs could efficiently target cerebral ischemia sites in vivo and reduce the infarct volume in mice. Our acoustothermal transfection method addresses a key bottleneck in balancing the transfection efficiency and cell viability, which can become a powerful tool in the future for cellular and gene therapies.

Closed Formula for Transport across Constrictions

In the last decade, the Fick-Jacobs approximation has been exploited to capture transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that, for fore-aft symmetric channels, the flux of noninteracting systems is fully captured by its linear response regime. For this case, we derive a very simple formula that captures the correct trends and can be exploited as a simple tool to design experiments or simulations. Lastly, we show that higher-order corrections in the flux may appear for nonsymmetric channels.

Multi-scale model of clogging in microfluidic devices with grid-like geometries

We propose a coarse-grained theoretical model to capture the ageing of microfluidic devices under different conditions including constant applied flow rate and constant applied pressure gradient. Microfluidic devices that sort cells by their deformability hold significant promise for medical applications. However, clogging in these microfluidic systems causes their properties to change over time and potentially limits their reliability. We compare the results of the coarse-grained model with those of stochastic simulations and with existing theoretical studies. Lastly, we apply the model to experimental data on the clogging of sickle red blood cells and discuss its wider applicability.

Impacts of Frazil Ice on the Effectiveness of Oil Dispersion and Migration of Dispersed Oil

Oil spills in the Arctic have drawn dramatic attention in recent years. Frazil ice, as the essential formation of sea ice, may affect the effectiveness of dispersants during oil spill response and the associated behaviors of dispersed oil. However, these impacts remain poorly understood, limiting the appropriate usage of dispersants in ice-covered regions. Herein this work explored the effects of frazil ice on the dispersion effectiveness of two dispersants (Corexit 9500A and hydrolyzed shrimp waste) and the migration of dispersed oil within frazil ice. We discovered that frazil ice inhibited dispersion effectiveness by attenuating water velocity. Permeable frazil ice encapsulated 11-30% of dispersed oil, implying a lower oil bioavailability. We thus proposed and verified a microscopic mechanism to unravel the migration of dispersed oil toward permeable constrictions in frazil ice. We predicted the concentration of dispersed oil encapsulated in frazil ice using bed filtration theory and verified the prediction through experiments. Furthermore, the presence of frazil ice can lead to the breakup and coalescence of dispersed oil. Overall, our findings would facilitate the appropriate planning and decision-making of dispersant-based oil spill response and a better understanding of the fate of dispersed oil in the frazil ice-infested ocean.

Squeezing multiple soft particles into a constriction: Transition to clogging

We study numerically how multiple deformable capsules squeeze into a constriction. This situation is largely encountered in microfluidic chips designed to manipulate living cells, which are soft entities. We use fully three-dimensional simulations based on the lattice Boltzmann method to compute the flow of the suspending fluid and on the immersed boundary method to achieve the two-way fluid-structure interaction. The mechanics of the capsule membrane elasticity is computed with the finite-element method. We obtain two main states: continuous passage of the particles and their blockage that leads to clogging the constriction. The transition from one state to another is dictated by the ratio between the size of the capsules and the constriction width and by the capsule membrane deformability. The latter is found to enhance particle passage through narrower constrictions, where rigid particles with similar diameter are blocked and lead to clogging.

Nanoscale Sub-Compartmentalization of the Dendritic Spine Compartment

Compartmentalization of the membrane is essential for cells to perform highly specific tasks and spatially constrained biochemical functions in topographically defined areas. These membrane lateral heterogeneities range from nanoscopic dimensions, often involving only a few molecular constituents, to micron-sized mesoscopic domains resulting from the coalescence of nanodomains. Short-lived domains lasting for a few milliseconds coexist with more stable platforms lasting from minutes to days. This panoply of lateral domains subserves the great variety of demands of cell physiology, particularly high for those implicated in signaling. The dendritic spine, a subcellular structure of neurons at the receiving (postsynaptic) end of central nervous system excitatory synapses, exploits this compartmentalization principle. In its most frequent adult morphology, the mushroom-shaped spine harbors neurotransmitter receptors, enzymes, and scaffolding proteins tightly packed in a volume of a few femtoliters. In addition to constituting a mesoscopic lateral heterogeneity of the dendritic arborization, the dendritic spine postsynaptic membrane is further compartmentalized into spatially delimited nanodomains that execute separate functions in the synapse. This review discusses the functional relevance of compartmentalization and nanodomain organization in synaptic transmission and plasticity and exemplifies the importance of this parcelization in various neurotransmitter signaling systems operating at dendritic spines, using two fast ligand-gated ionotropic receptors, the nicotinic acetylcholine receptor and the glutamatergic receptor, and a second-messenger G-protein coupled receptor, the cannabinoid receptor, as paradigmatic examples.

Dynamics of Vesicles Driven Into Closed Constrictions by Molecular Motors

We study the dynamics of a model of membrane vesicle transport into dendritic spines, which are bulbous intracellular compartments in neurons driven by molecular motors. We reduce the lubrication model proposed in Fai et al. (Phys Rev Fluids 2:113601, 2017) to a fast-slow system, yielding an analytically and numerically tractable equation equivalent to the original model in the overdamped limit. The model's key parameters include: (1) the ratio of motors that prefer to push toward the head of the dendritic spine to the motors that prefer to push in the opposite direction, and (2) the viscous drag exerted on the vesicle by the spine constriction. We perform a numerical bifurcation analysis in these parameters and find that steady-state vesicle velocities appear and disappear through several saddle-node bifurcations. This process allows us to identify the region of parameter space in which multiple stable velocities exist. We show by direct calculations that there can only be unidirectional motion for sufficiently close vesicle-to-spine diameter ratios. Our analysis predicts the critical vesicle-to-spine diameter ratio, at which there is a transition from unidirectional to bidirectional motion, consistent with experimental observations of vesicle trajectories in the literature.

Cooperative size sorting of deformable particles in porous media

Diverse applications-ranging from enhanced oil recovery, filtration, and lab on a chip sorting-rely on the flow-induced transport of deformable particles in porous media. However, how fluid flow can force such particles to squeeze through pore constrictions of complex geometries is poorly understood. Here, we study the transport of model deformable particles in millifluidic porous media with constrictions of tunable aspect ratio. We find that multiple particles can unexpectedly squeeze through large-aspect ratio constrictions, even when isolated particles cannot. This phenomenon arises from pairwise flow-mediated interactions between the particles: when one particle is trapped at a constriction, the increased fluid flow around it enables a second to squeeze past due to locally increased hydrodynamic stresses. This cooperative mechanism causes the particles to ultimately sort themselves by size through the pore space. By revealing a new mode of deformable particle transport in porous media, our work helps to inform real-world applications and provides a straightforward way to sort particles based on size.

Transport mechanism of deformable micro-gel particle through micropores with mechanical properties characterized by AFM

Deformable micro-gel particles (DMP) have been used to enhanced oil recovery (EOR) in reservoirs with unfavourable conditions. Direct pore-scale understanding of the DMP transport mechanism is important for further improvements of its EOR performance. To consider the interaction between soft particle and fluid in complex pore-throat geometries, we perform an Immersed Boundary-Lattice Boltzmann (IB-LB) simulation of DMP passing through a throat. A spring-network model is used to capture the deformation of DMP. In order to obtain appropriate simulation parameters that represent the real mechanical properties of DMP, we propose a procedure via fitting the DMP elastic modulus data measured by the nano-indentation experiments using Atomic Force Microscope (AFM). The pore-scale modelling obtains the critical pressure of the DMP for different particle-throat diameter ratios and elastic modulus. It is found that two-clog particle transport mode is observed in a contracted throat, the relationship between the critical pressure and the elastic modulus/particle-throat diameter ratio follows a power law. The particle-throat diameter ratio shows a greater impact on the critical pressure difference than the elastic modulus of particles.

Mechanical instability and percolation of deformable particles through porous networks

The transport of micron-sized particles such as bacteria, cells, or synthetic lipid vesicles through porous spaces is a process relevant to drug delivery, separation systems, or sensors, to cite a few examples. Often, the motion of these particles depends on their ability to squeeze through small constrictions, making their capacity to deform an important factor for their permeation. However, it is still unclear how the mechanical behavior of these particles affects collective transport through porous networks. To address this issue, we present a method to reconcile the pore-scale mechanics of the particles with the Darcy scale to understand the motion of a deformable particle through a porous network. We first show that particle transport is governed by a mechanical instability occurring at the pore scale, which leads to a binary permeation response on each pore. Then, using the principles of directed bond percolation, we are able to link this microscopic behavior to the probability of permeating through a random porous network. We show that this instability, together with network uniformity, are key to understanding the nonlinear permeation of particles at a given pressure gradient. The results are then summarized by a phase diagram that predicts three distinct permeation regimes based on particle properties and the randomness of the pore network.

A Low-Backpressure Single-Cell Point Constriction for Cytosolic Delivery Based on Rapid Membrane Deformations

Mechanically deforming biological cells through microfluidic constrictions is a recently introduced technique for the intracellular delivery of macromolecules possibly through transient membrane pores induced in the process. The technique is attractive for research and clinical applications mainly because it is simple, fast, and effective while being free of adverse effects often associated with well-known techniques that rely on field- or vector-based delivery. In this nascent approach, an utmost and crucial role is played by the constriction, often in rectangular profile, and it squeezes cells only in one dimension. The results achieved suggest that the longer the constriction is the higher the delivery performance. Contrary to this view, we demonstrate here a unique constriction profile that is highly localized (point) and yet returns comparably effective delivery. Point constrictions are of a semiround geometry, forcing cells in both dimensions while introducing very little backpressure to the system, which is a silicon-glass platform wherein constrictions are arranged in series along an array of channels. The influence of the constriction size and count as well as treatment pressure on delivery performance is presented on the basis of the flow-cytometric analyses of HCT116 cells treated using dextran as model molecules. Delivery performance is also presented for common mammalian cell lines including NIH 3T3, HEK293, and MDCK. Moreover, the versatility of the platform is demonstrated in gene knockdown experiments using synthetic siRNA as well as on the delivery of proteins. Target proteins in some cells exhibit nondiffusive distribution profile raising the plausibility of mechanisms other than transient membrane pores.

Off-center motion of a trapped elastic capsule in a microfluidic channel with a narrow constriction

Owing to their significance in capsule-related engineering and biomedical applications, a number of studies have considered the dynamics of elastic capsules flowing in constricted microchannels. However, these studies have focused on capsules moving along the channel centerline. In the present study, we numerically investigate the transient motion of an elastic capsule in a microfluidic channel with a rectangular constriction, which is initially trapped at the constriction inlet while off the channel centerline (i.e., on the channel bottom-wall). Under the push of the surrounding flow, the capsule can squeeze into the constriction, but only if the capsule deformability or the constriction size is sufficiently large. We find that the critical capillary number leading to the penetration of the capsule into the constriction is larger for off-centerline capsules compared to centered capsules. The centered capsule is stationary at the steady state when it remains stuck at the constriction; in contrast, the off-centerline capsule is not stationary but exhibits a tank-treading motion, i.e., its overall shape maintains a nonspherical shape with a protrusion into the constriction while its membrane exhibits a continuous rotation. Further, we examine the dependence of the capsule motion type, capsule deformation degree and membrane tension distribution on the capillary number (measuring the effects of flow strength and membrane mechanics) and constriction geometries (including the constriction height and width). Finally, we discuss the mechanism governing the capsule motion by analyzing the hydrodynamic forces acting on the capsule. The shear force acting on the capsule top owing to the fluid flow in the gap between the capsule top and the channel top-wall is the main source inducing the membrane tank-treading rotation.

Squeezing bio-capsules into a constriction: deformation till break-up

We study experimentally the deformation and break-up of liquid-filled capsules trapped at an axisymmetric step constriction, and subjected to increasing pressure drops. We considered biological (trout fish eggs) and bioartificial (made of ovalbumin and alginate) ones, with the objective to characterize the transition to break-up. We find that both capsule populations behave as a brittle material. They do not exhibit any plastic deformation prior to break-up. Moreover critical pressure drop exhibits a stochastic behavior as known for the fracture of disordered media. The break-up probability follows a three-parameter Weibull distribution, from which one can deduce the capsule rupture characteristics.

Effect of body deformability on microswimming

In this work we consider the following question: given a mechanical microswimming mechanism, does increased deformability of the swimmer body hinder or promote the motility of the swimmer? To answer this we run immersed-boundary-lattice-Boltzmann simulations of a microswimmer composed of deformable beads connected with springs. We find that the same deformations in the beads can result in different effects on the swimming velocity, namely an enhancement or a reduction, depending on the other parameters. To understand this we determine analytically the velocity of the swimmer, starting from the forces driving the motion and assuming that the deformations in the beads are known as functions of time and are much smaller than the beads themselves. We find that to the lowest order, only the driving frequency mode of the surface deformations contributes to the swimming velocity, and comparison to the simulations shows that both the velocity-promoting and velocity-hindering effects of bead deformability are reproduced correctly by the theory in the limit of small bead deformations. For the case of active deformations we show that there are critical values of the spring constant - which for a general swimmer corresponds to its main elastic degree of freedom - which decide whether the body deformability is beneficial for motion or not.

Theory and algorithms to compute Helfrich bending forces: a review

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

Motion of an elastic capsule in a constricted microchannel

We study the motion of an elastic capsule through a microchannel characterized by a localized constriction. We consider a capsule with a stress-free spherical shape and impose its steady-state configuration in an infinitely long straight channel as the initial condition for our calculations. We report how the capsule deformation, velocity, retention time, and maximum stress of the membrane are affected by the capillary number, Ca , and the constriction shape. We estimate the deformation by measuring the variation of the three-dimensional surface area and a series of alternative quantities easier to extract from experiments. These are the Taylor parameter, the perimeter and the area of the capsule in the spanwise plane. We find that the perimeter is the quantity that best reproduces the behavior of the three-dimensional surface area. This is maximum at the centre of the constriction and shows a second peak after it, whose location depends on the Ca number. We observe that, in general, area-deformation-correlated quantities grow linearly with Ca , while velocity-correlated quantities saturate for large Ca but display a steeper increase for small Ca . The velocity of the capsule divided by the velocity of the flow displays, surprisingly, two different qualitative behaviors for small and large capillary numbers. Finally, we report that longer constrictions and spanwise wall bounded (versus spanwise periodic) domains cause larger deformations and velocities. If the deformation and velocity in the spanwise wall bounded domains are rescaled by the initial equilibrium deformation and velocity, their behavior is undistinguishable from that in a periodic domain. In contrast, a remarkably different behavior is reported in sinusoidally shaped and smoothed rectangular constrictions indicating that the capsule dynamics is particularly sensitive to abrupt changes in the cross section. In a smoothed rectangular constriction larger deformations and velocities occur over a larger distance.

Barriers in the brain: resolving dendritic spine morphology and compartmentalization

Dendritic spines are micron-sized protrusions that harbor the majority of excitatory synapses in the central nervous system. The head of the spine is connected to the dendritic shaft by a 50-400 nm thin membrane tube, called the spine neck, which has been hypothesized to confine biochemical and electric signals within the spine compartment. Such compartmentalization could minimize interspinal crosstalk and thereby support spine-specific synapse plasticity. However, to what extent compartmentalization is governed by spine morphology, and in particular the diameter of the spine neck, has remained unresolved. Here, we review recent advances in tool development - both experimental and theoretical - that facilitate studying the role of the spine neck in compartmentalization. Special emphasis is given to recent advances in microscopy methods and quantitative modeling applications as we discuss compartmentalization of biochemical signals, membrane receptors and electrical signals in spines. Multidisciplinary approaches should help to answer how dendritic spine architecture affects the cellular and molecular processes required for synapse maintenance and modulation.