The force generated by biological membranes on a polymer rod and its response: Statics and dynamics

Islands of conformational stability for filopodia

Filopodia are long, thin protrusions formed when bundles of fibers grow outwardly from a cell surface while remaining closed in a membrane tube. We study the subtle issue of the mechanical stability of such filopodia and how this depends on the deformation of the membrane that arises when the fiber bundle adopts a helical configuration. We calculate the ground state conformation of such filopodia, taking into account the steric interaction between the membrane and the enclosed semiflexible fiber bundle. For typical filopodia we find that a minimum number of fibers is required for filopodium stability. Our calculation elucidates how experimentally observed filopodia can obviate the classical Euler buckling condition and remain stable up to several tens of . We briefly discuss how experimental observation of the results obtained in this work for the helical-like deformations of enclosing membrane tubes in filopodia could possibly be observed in the acrosomal reactions of the sea cucumber Thyone, and the horseshoe crab Limulus. Any realistic future theories for filopodium stability are likely to rely on an accurate treatment of such steric effects, as analysed in this work.

The physical foundation of vasoocclusion in sickle cell disease

The pathology of sickle cell disease arises from the occlusion of small blood vessels because of polymerization of the sickle hemoglobin within the red cells. We present measurements using a microfluidic method we have developed to determine the pressure required to eject individual red cells from a capillary-sized channel after the cell has sickled. We find that the maximum pressure is only ∼100 Pa, much smaller than typically found in the microcirculation. This explains why experiments using animal models have not observed occlusion beginning in capillaries. The magnitude of the pressure and its dependence on intracellular concentration are both well described as consequences of sickle hemoglobin polymerization acting as a Brownian ratchet. Given the recently determined stiffness of sickle hemoglobin gels, the observed obstruction seen in sickle cell disease as mediated by adherent cells can now be rationalized, and surprisingly suggests a window of maximum vulnerability during circulation of sickle cells.

Predicting the morphology of sickle red blood cells using coarse-grained models of intracellular aligned hemoglobin polymers

Sickle red blood cells (SS-RBCs) exhibit heterogeneous cell morphologies (sickle, holly leaf, granular, etc.) in the deoxygenated state due to the polymerization of the sickle hemoglobin. Experimental evidence points to a close relationship between SS-RBC morphology and intracellular aligned hemoglobin polymers. Here, we develop a coarse-grained (CG) stochastic model to represent the growth of the intracellular aligned hemoglobin polymer domain. The CG model is calibrated based on the mechanical properties (Young's modulus, bending rigidity) of the sickle hemoglobin fibers reported in experiments. The process of the cell membrane transition is simulated for physiologic aligned hemoglobin polymer configurations and mean corpuscular hemoglobin concentration. Typical SS-RBC morphologies observed in experiments can be obtained from the current model as a result of the intracellular aligned hemoglobin polymer development without introducing any further ad hoc assumptions. It is found that the final shape of SS-RBCs is primarily determined by the angular width of the aligned hemoglobin polymer domain, but it also depends, to a lesser degree, on the polymer growth rate and the cell membrane rigidity. Cell morphologies are quantified by structural shape factors, which agree well with experimental results from medical images.

Dynamic structure factor of a stiff polymer in a glassy solution

We provide a comprehensive overview of the current theoretical understanding of the dynamic structure factor of stiff polymers in semidilute solution based on the wormlike chain (WLC) model. We extend previous work by computing exact numerical coefficients and an expression for the dynamic mean square displacement (MSD) of a free polymer and compare various common approximations for the hydrodynamic interactions, which need to be treated accurately if one wants to extract quantitative estimates for model parameters from experimental data. A recent controversy about the initial slope of the dynamic structure factor is thereby resolved. To account for the interactions of the polymer with a surrounding (sticky) polymer solution, we analyze an extension of the WLC model, the glassy wormlike chain (GWLC), which predicts near power law and logarithmic long-time tails in the dynamic structure factor.

Universal metastability of sickle hemoglobin polymerization

Sickle hemoglobin (HbS) polymerization occurs when the concentration of deoxyHbS exceeds a well-defined solubility. In experiments using sickle hemoglobin droplets suspended in oil, it has been shown that when polymerization ceases the monomer concentration is above equilibrium solubility. We find that the final concentration in uniform bulk solutions (i.e., with negligible boundaries) agrees with the droplet measurements, and both exceed the expected solubility. To measure hemoglobin in uniform solutions, we used modulated excitation of trace amounts of CO in gels of HbS. In this method, a small amount of CO is introduced to a spatially uniform deoxyHb sample, so that less than 2% of the sample is liganded. The liganded fraction is photolyzed repeatedly and the rate of recombination allows the concentration of deoxyHbS in the solution phase to be determined, even if polymers have formed. Both uniform and droplet samples exhibit the same quantitative behavior, exceeding solubility by an amount that depends on the initial concentration of the sample, as well as conditions under which the gel was formed. We hypothesize that the early termination of polymerization is due to the obstruction in polymer growth, which is consistent with the observation that pressing on slides lowers the final monomer concentration, making it closer to solubility. The thermodynamic solubility in free solution is thus achieved only in conditions with low polymer density or under external forces (such as found in sedimentation) that disrupt polymers. Since we find that only about 67% of the expected polymer mass forms, this result will impact any analysis predicated on predicting the polymer fraction in a given experiment.

Monomer depletion, pressure difference, and membrane tube radius reduction due to fiber polymerization in microspikes

In many processes vital to life, the growth of biological fibers outwards from a membrane surface naturally produces membrane tube tethers or microspikes in biological cells. Here, we investigate the novel effect of pressure difference (due to monomer depletion) on the polymerization dynamics of biological fibers within long membrane tubes. We crucially find that fiber monomers become depleted close to the growing tip as the fiber polymerizes, thus reducing the local pressure, and hence decreasing the membrane tube radius at the tip. This process is found to slow the growth of the fiber, a process which becomes important when we go on to construct a dynamical theory for biopolymer growth in long, narrow tubes. Our result is interesting in that it emphasizes how "passive" biological transport mechanisms such as via pressure differences may play an important role in cell movements.

Deforming biological membranes: how the cytoskeleton affects a polymerizing fiber

We give a theoretical treatment of the force exerted by a fluctuating membrane on a polymer rod tip, taking into account the effects of an underlying biological cytoskeleton by way of a simple harmonic dependence on displacement. We also consider theoretically and experimentally the dynamics of a growing fiber tip under the influence of such a fluctuation-induced membrane force, including the effects of an underlying cytoskeletal network. We compare our model with new experimental data for the growth of hemoglobin fibers within red blood cells, revealing a good agreement. We are also able to estimate the force and membrane/cytoskeletal displacement required to stall growth of, or buckle, a growing fiber. We discuss the significance of our results in a biological context, including how the properties of the membrane and cytoskeleton relate to the thermodynamics of rod polymerization.

Metastable polymerization of sickle hemoglobin in droplets

Sickle cell disease arises from a genetic mutation of one amino acid in each of the two hemoglobin beta chains, leading to the polymerization of hemoglobin in the red cell upon deoxygenation, and is characterized by vascular crises and tissue damage due to the obstruction of small vessels by sickled cells. It has been an untested assumption that, in red cells that sickle, the growing polymer mass would consume monomers until the thermodynamically well-described monomer solubility was reached. By photolysing droplets of sickle hemoglobin suspended in oil we find that polymerization does not exhaust the available store of monomers, but stops prematurely, leaving the solutions in a supersaturated, metastable state typically 20% above solubility at 37 degrees C, though the particular values depend on the details of the experiment. We propose that polymer growth stops because the growing ends reach the droplet edge, whereas new polymer formation is thwarted by long nucleation times, since the concentration of hemoglobin is lowered by depletion of monomers into the polymers that have formed. This finding suggests a new aspect to the pathophysiology of sickle cell disease; namely, that cells deoxygenated in the microcirculation are not merely undeformable, but will actively wedge themselves tightly against the walls of the microvasculature by a ratchet-like mechanism driven by the supersaturated solution.

The cellulose synthase complex: a polymerization driven supramolecular motor

We present a biophysical model for the propulsion of the cellulose synthase complex, the motile transmembrane protein complex responsible for the biosynthesis of cellulose microfibrils, the dominant architectural component of the cell walls of higher plants. Our model identifies the polymerization and the crystallization of the cellulose chains as the combined driving forces and elucidates the role of polymer flexibility and membrane elasticity as force transducers. The model is elaborated using both stochastic simulations and a simplified analytical treatment. On the basis of the model and approximate values for the relevant physical constants, we estimate the speed of the cellulose synthase complex to be in the range v(p) approximately 10(-9)-10(-8) m/s, consistent with the recently reported experimental value of 5.8 x 10(-9) m/s.

Stall, spiculate, or run away: The fate of fibers growing towards fluctuating membranes

We study the dynamics of a growing semiflexible fiber approaching a membrane at an angle. At late times we find three regimes: fiber stalling, when growth stops, runaway, in which the fiber bends away from the membrane, and another regime in which spicules form. We discuss which regions of the resulting "phase diagram" are explored by (i) single and bundled actin fibers in living cells, (ii) sickle hemoglobin fibers, and (iii) microtubules inside vesicles. We complement our analysis with 3D stochastic simulations.

Spicules and the effect of rigid rods on enclosing membrane tubes

Membrane tubes (spicules) arise in cells, or artificial membranes, in the nonlinear deformation regime due to, e.g., the growth of microtubules, actin filaments, or sickle hemoglobin fibers towards a membrane. We calculate the axial force f exerted by the tube, and its average radius, taking into account steric interactions between the fluctuating membrane and the enclosed rod. We find a smooth crossover of the axial force between f approximately square root of (sigma) and f approximately sigma as the membrane tension sigma increases and the tube radius shrinks. This crossover occurs around the most physiologically relevant membrane tensions. Our work may be important in (i) interpreting experiments in which axial force is related to the tube radius or membrane tension, and (ii) constructing dynamical theories for biopolymer growth in narrow tubes where these fluctuation effects control the tube radius.

Three-dimensional dynamic Monte Carlo simulations of elastic actin-like ratchets

We present three-dimensional dynamic Monte Carlo simulations of the growth of a semiflexible fiber against a fluctuating obstacle. The natural reference for our numerical study are the elastic and Brownian ratchet models previously analyzed semianalytically. We find that the decay of the velocity versus applied load is exponential to a good degree of accuracy, provided we include in the load the drag force felt by the moving obstacle. If the fiber and obstacle only interact via excluded volume, there are small corrections to the Brownian ratchet predictions which suggest that tip fluctuations play a minor role. If on the other hand fiber and obstacle interact via a soft potential, the corrections are much larger when the obstacle diffuses slowly. This means that microscopic assumptions can profoundly affect the dynamics. We also identify and characterize a novel "pushing catastrophe"--which is distinct from the usual fiber buckling--in which the growth of the fiber decouples from the obstacle movement. The time distribution of catastrophes can be explained via an approximate analytical treatment, and our numerics suggest that the time taken to lose propulsive force is largely dependent on the fiber incidence angle. Our results are a first step in realizing numerical polymer models for the motion of sets or networks of semiflexible fibers close to a fluctuating membrane or obstacle.