Cells as Active Particles in Asymmetric Potentials: Motility under External Gradients

Cell migration is a crucial event during development and in disease. Mechanical constraints and chemical gradients can contribute to the establishment of cell direction, but their respective roles remain poorly understood. Using a microfabricated topographical ratchet, we show that the nucleus dictates the direction of cell movement through mechanical guidance by its environment. We demonstrate that this direction can be tuned by combining the topographical ratchet with a biochemical gradient of fibronectin adhesion. We report competition and cooperation between the two external cues. We also quantitatively compare the measurements associated with the trajectory of a model that treats cells as fluctuating particles trapped in a periodic asymmetric potential. We show that the cell nucleus contributes to the strength of the trap, whereas cell protrusions guided by the adhesive gradients add a constant tunable bias to the direction of cell motion.     

Analytical aspects of the Brownian motor effect in randomly flashing ratchets

The muscle contraction, operation of ATP synthase, maintaining the shape of a cell are believed to be secured by motor proteins, which can be modelled using the Brownian ratchet mechanism. We consider the randomly flashing ratchet model of a Brownian motor, where the particles can be in two states, only one of which is sensitive the applied spatially periodic potential (the mathematical setting is a pair of weakly coupled reaction-diffusion and Fokker–Planck equations). We prove that this mechanism indeed generates unidirectional transport by showing that the amount of mass in the wells of the potential decreases/increases from left to right. The direction of transport is unambiguously determined by the location of each minimum of the potential with respect to the so-called diffusive mean of its adjacent maxima. The transport can be generated not only by an asymmetric potential, but also by a symmetric potential and asymmetric transition rates, and as a consequence of the general result we derive explicit conditions when the latter happens. When the transitions are localized on narrow active sites in the protein conformation space, we find a more explicit characterization of the bulk transport direction, and infer that some common preconditions of the motor effect are redundant.     

Automated topographical cell proliferation analysis

BACKGROUND: Cell proliferation is often studied using the incorporation of bromodeoxyuridine (BrdU). Immunohistochemical staining is then used to detect BrdU in the nucleus. To circumvent the observer bias and labor-intensive nature of manually counting BrdU-labeled nuclei, an automated topographical cell proliferation analysis method is developed. METHODS: Sections stained with fluorescein-labeled anti-BrdU and counterstained with To-Pro-3 are scanned using confocal laser scanning microscopy (CLSM). For every point in the image, the nucleus density of BrdU-labeled nuclei and the total nucleus density of the neighborhood of that point are calculated from the BrdU and the To-Pro-3 signal, respectively. The ratio of these densities gives an indication of the amount of cell proliferation at that point. The automated measure is validated by comparing it with the ratio of BrdU-stained nuclei to the total number of nuclei obtained from a manual count. RESULTS: A positive correlation is found between the automated measure and the ratios calculated from the manual counting (r = 0.86, P < 0.001). Calculating the topographical cell proliferation using the automated method is faster and does not suffer from interobserver variability. CONCLUSIONS: Automated topographical cell proliferation analysis is a fast method to objectively find differences in cell proliferation within a tissue. This can be visualized by a topographical map that corresponds to the tissue under study.     

Effects of Surface Asymmetry on Neuronal Growth

Detailed knowledge of how the surface physical properties, such as mechanics, topography and texture influence axonal outgrowth and guidance is essential for understanding the processes that control neuron development, the formation of functional neuronal connections and nerve regeneration. Here we synthesize asymmetric surfaces with well-controlled topography and texture and perform a systematic experimental and theoretical investigation of axonal outgrowth on these substrates. We demonstrate unidirectional axonal bias imparted by the surface ratchet-based topography and quantify the topographical guidance cues that control neuronal growth. We describe the growth cone dynamics using a general stochastic model (Fokker-Planck formalism) and use this model to extract two key dynamical parameters: diffusion (cell motility) coefficient and asymmetric drift coefficient. The drift coefficient is identified with the torque caused by the asymmetric ratchet topography. We relate the observed directional bias in axonal outgrowth to cellular contact guidance behavior, which results in an increase in the cell-surface coupling with increased surface anisotropy. We also demonstrate that the disruption of cytoskeletal dynamics through application of Taxol (stabilizer of microtubules) and Blebbistatin (inhibitor of myosin II activity) greatly reduces the directional bias imparted by these asymmetric surfaces. These results provide new insight into the role played by topographical cues in neuronal growth and could lead to new methods for stimulating neuronal regeneration and the engineering of artificial neuronal tissue.     

Getting in shape and swimming: the role of cortical forces and membrane heterogeneity in eukaryotic cells

Recent research has shown that motile cells can adapt their mode of propulsion to the mechanical properties of the environment in which they find themselves—crawling in some environments while swimming in others. The latter can involve movement by blebbing or other cyclic shape changes, and both highly-simplified and more realistic models of these modes have been studied previously. Herein we study swimming that is driven by membrane tension gradients that arise from flows in the actin cortex underlying the membrane, and does not involve imposed cyclic shape changes. Such gradients can lead to a number of different characteristic cell shapes, and our first objective is to understand how different distributions of membrane tension influence the shape of cells in an inviscid quiescent fluid. We then analyze the effects of spatial variation in other membrane properties, and how they interact with tension gradients to determine the shape. We also study the effect of fluid–cell interactions and show how tension leads to cell movement, how the balance between tension gradients and a variable bending modulus determine the shape and direction of movement, and how the efficiency of movement depends on the properties of the fluid and the distribution of tension and bending modulus in the membrane.

     

Protrusion Fluctuations Direct Cell Motion

Many physiological phenomena involve directional cell migration. It is usually attributed to chemical gradients in vivo. Recently, other cues have been shown to guide cells in vitro, including stiffness/adhesion gradients or micropatterned adhesive motifs. However, the cellular mechanism leading to these biased migrations remains unknown, and, often, even the direction of motion is unpredictable. In this study, we show the key role of fluctuating protrusions on ratchet-like structures in driving NIH3T3 cell migration. We identified the concept of efficient protrusion and an associated direction index. Our analysis of the protrusion statistics facilitated the quantitative prediction of cell trajectories in all investigated conditions. We varied the external cues by changing the adhesive patterns. We also modified the internal cues using drug treatments, which modified the protrusion activity. Stochasticity affects the short- and long-term steps. We developed a theoretical model showing that an asymmetry in the protrusion fluctuations is sufficient for predicting all measures associated with the long-term motion, which can be described as a biased persistent random walk.     

Direction reversal of fluctuation-induced biased Brownian motion on distorted ratchets

The biased movement of Brownian particles on a fluctuating two-state periodic potential made of identical distorted ratchets is studied. The purpose is to investigate how the direction of the particle movement is related to the asymmetry of the potential. In general, distorting one of the two linear arms of a regular symmetric ratchet (with equal arm lengths) can create a driving force for the Brownian particle to execute biased movement. The direction of the induced biased movement depends on the type of the distortion. It has been found that if one linear arm is kinked into two linear sub-arms, the direction of the movement can be either positive or negative depending on the frequency of the fluctuation and the location and the degree of the kink. In contrast, if one arm of the symmetric ratchet is replaced by a continuous nonlinear sinusoidal function, the movement is always unidirectional. Thus, for the latter case to generate the direction reversal phenomenon, the ratchets have to have an additional asymmetry. We also have found that two potentials with different distorted ratchets can generate identical fluxes if the distortions are polar symmetric about the mid-point of the arm(s) of the basic linear two-arm ratchet. The results are useful for designing experimental apparatuses for the separation of protein particles based on their sizes and charges and the viscosity of the medium.     

A stochastic model for chemotaxis based on the ordered extension of pseudopods

Many amoeboid cells move by extending pseudopods. Here I present a new stochastic model for chemotaxis that is based on pseudopod extensions by Dictyostelium cells. In the absence of external cues, pseudopod extension is highly ordered with two types of pseudopods: de novo formation of a pseudopod at the cell body in random directions, and alternating right/left splitting of an existing pseudopod that leads to a persistent zig-zag trajectory. We measured the directional probabilities of the extension of splitting and de novo pseudopods in chemoattractant gradients with different steepness. Very shallow cAMP gradients can bias the direction of splitting pseudopods, but the bias is not perfect. Orientation of de novo pseudopods require much steeper cAMP gradients and can be more precise. These measured probabilities of pseudopod directions were used to obtain an analytical model for chemotaxis of cell populations. Measured chemotaxis of wild-type cells and mutants with specific defects in these stochastic pseudopod properties are similar to predictions of the model. These results show that combining splitting and de novo pseudopods is a very effective way for cells to obtain very high sensitivity to stable gradient and still be responsive to changes in the direction of the gradient.     

Protrusion Fluctuations Direct Cell Motion

Many physiological phenomena involve directional cell migration. It is usually attributed to chemical gradients in vivo. Recently, other cues have been shown to guide cells in vitro, including stiffness/adhesion gradients or micropatterned adhesive motifs. However, the cellular mechanism leading to these biased migrations remains unknown, and, often, even the direction of motion is unpredictable. In this study, we show the key role of fluctuating protrusions on ratchet-like structures in driving NIH3T3 cell migration. We identified the concept of efficient protrusion and an associated direction index. Our analysis of the protrusion statistics facilitated the quantitative prediction of cell trajectories in all investigated conditions. We varied the external cues by changing the adhesive patterns. We also modified the internal cues using drug treatments, which modified the protrusion activity. Stochasticity affects the short- and long-term steps. We developed a theoretical model showing that an asymmetry in the protrusion fluctuations is sufficient for predicting all measures associated with the long-term motion, which can be described as a biased persistent random walk.     

Protein-protein ratchets: stochastic simulation and application to processive enzymes

Interaction between a protein and a series of binding sites on a cytoskeletal substrate can create a resistance, or "protein friction," as the protein is moved along the substrate. If attachment and detachment rates are specified asymmetrically, this resistance can depend on the direction of movement, and the binding interaction acts as a ratchet. Stochastic computer simulations have been used to examine this type of protein-protein interaction. The performance of a protein-protein ratchet in the piconewton and nanometer range is significantly limited by thermal fluctuations, which in experimental measurements with single molecules are evident as Brownian motion. Simulations with a two-component model combining a conventional motor enzyme model with a protein-protein ratchet confirm previous suggestions that the processive movement of a single motor enzyme molecule against a load, as seen in experiments with inner arm dynein molecules, might be made possible by an accessory protein interaction that prevents backward slippage. When this accessory protein interaction is defined so that it acts as a ratchet, backward slippage can be prevented with minimal interference with forward progression.     

Fluctuations of Fitness Distributions and the Rate of Muller's Ratchet

The accumulation of deleterious mutations is driven by rare fluctuations that lead to the loss of all mutation free individuals, a process known as Muller's ratchet. Even though Muller's ratchet is a paradigmatic process in population genetics, a quantitative understanding of its rate is still lacking. The difficulty lies in the nontrivial nature of fluctuations in the fitness distribution, which control the rate of extinction of the fittest genotype. We address this problem using the simple but classic model of mutation selection balance with deleterious mutations all having the same effect on fitness. We show analytically how fluctuations among the fittest individuals propagate to individuals of lower fitness and have dramatically amplified effects on the bulk of the population at a later time. If a reduction in the size of the fittest class reduces the mean fitness only after a delay, selection opposing this reduction is also delayed. This delayed restoring force speeds up Muller's ratchet. We show how the delayed response can be accounted for using a path-integral formulation of the stochastic dynamics and provide an expression for the rate of the ratchet that is accurate across a broad range of parameters.     

Temporal variation in selection accelerates mutational decay by Muller's ratchet

Asexual species accumulate deleterious mutations through an irreversible process known as Muller's ratchet. Attempts to quantify the rate of the ratchet have ignored the role of temporal environmental heterogeneity even though it is common in nature and has the potential to affect overall ratchet rate. Here we examine Muller's ratchet in the context of conditional neutrality (i.e., mutations that are deleterious in some environmental conditions but neutral in others) as well as more subtle changes in the strength (but not sign) of selection. We find that temporal variation increases the rate of the ratchet (mutation accumulation) and the rate of fitness decline over that of populations experiencing constant selection of equivalent average strength. Temporal autocorrelation magnifies the effects of temporal heterogeneity and can allow the ratchet to operate at large population sizes in which it would be halted under constant selection. Classic studies of Muller's ratchet show that the rate of fitness decline is maximized when selection is of a low but intermediate strength. This relationship changes quantitatively with all forms of temporal heterogeneity studied and changes qualitatively when there is temporal autocorrelation in selection. In particular, the rate of fitness decline can increase indefinitely with the strength of selection with some forms of temporal heterogeneity. Our finding that temporal autocorrelation in selection dramatically increases ratchet rate and rate of fitness decline may help to explain the paucity of asexual taxa.     

Adaptive evolution of asexual populations under Muller's ratchet

We study the population genetics of adaptation in nonequilibrium haploid asexual populations. We find that the accumulation of deleterious mutations, due to the operation of Muller's ratchet, can considerably reduce the rate of fixation of advantageous alleles. Such reduction can be approximated reasonably well by a reduction in the effective population size. In the absence of Muller's ratchet, a beneficial mutation can only become fixed if it creates the best possible genotype; if Muller's ratchet operates, however, mutations initially arising in a nonoptimal genotype can also become fixed in the population, since the loss of the least-loaded class implies that an initially nonoptimal background can become optimal. We show that, while the rate at which adaptive mutations become fixed is reduced, the rate of fixation of deleterious mutations due to the ratchet is not changed by the presence of beneficial mutations as long as the rate of their occurrence is low and the deleterious effects of mutations (s(d)) are higher than the beneficial effects (s(a)). When s(a) > s(d), the advantage of a beneficial mutation can outweigh the deleterious effects of associated mutations. Under these conditions, a beneficial allele can drag to fixation deleterious mutations initially associated with it at a higher rate than in the absence of advantageous alleles. We propose analytical approximations for the rates of accumulation of deleterious and beneficial mutations. Furthermore, when allowing for the possible occurrence of interference between beneficial alleles, we find that the presence of deleterious mutations of either very weak or very strong effect can marginally increase the rate of accumulation of beneficial mutations over that observed in the absence of such deleterious mutations.     

Local and spatially coordinated movements in Dictyostelium discoideum amoebae during chemotaxis

We have studied chemotaxis by individual Dictyostelium discoideum amoebae using strong, local gradients of the chemoattractant cyclic AMP. Gradients were provided by diffusion of cyclic AMP from a microneedle, which could be positioned at various points around the cell. Responses to changes in the gradient indicate how the cell is structurally organized for chemotactic movement. There is a polarity in the responsiveness of the surface to stimulation by cyclic AMP along the length of the amoeba. Furthermore, two aspects of chemotactic movement can be distinguished. The first response to cyclic AMP is a locally generated extension of a hyaline pseudopod from the region of the surface nearest the stimulus. The second response, the flow of cytoplasm in the direction of the stimulus, is coordinated and separate from the first response. The coordination appears to depend on the nucleus or on the microtubule-organizing center.     

Loss of least-loaded class in asexual populations due to drift and epistasis

We consider the dynamics of a nonrecombining haploid population of finite size that accumulates deleterious mutations irreversibly. This ratchet-like process occurs at a finite speed in the absence of epistasis, but it has been suggested that synergistic epistasis can halt the ratchet. Using a diffusion theory, we find explicit analytical expressions for the typical time between successive clicks of the ratchet for both nonepistatic and epistatic fitness functions. Our calculations show that the interclick time is of a scaling form that in the absence of epistasis gives a speed that is determined by size of the least-loaded class and the selection coefficient. With synergistic interactions, the ratchet speed is found to approach zero rapidly for arbitrary epistasis. Our analytical results are in good agreement with the numerical simulations.     

A multiscale whole-cell theory for mechanosensitive migration on viscoelastic substrates

Increasing experimental evidence validates that both the elastic stiffness and viscosity of the extracellular matrix regulate mesenchymal cell behavior, such as the rational switch between durotaxis (cell migration to stiffer regions), anti-durotaxis (migration to softer regions), and adurotaxis (stiffness-insensitive migration). To reveal the mechanisms underlying the crossover between these motility regimes, we have developed a multiscale chemomechanical whole-cell theory for mesenchymal migration. Our framework couples the subcellular focal adhesion dynamics at the cell-substrate interface with the cellular cytoskeletal mechanics and the chemical signaling pathways involving Rho GTPase proteins. Upon polarization by the Rho GTPase gradients, our simulated cell migrates by concerted peripheral protrusions and contractions, a hallmark of the mesenchymal mode. The resulting cell dynamics quantitatively reproduces the experimental migration speed as a function of the uniform substrate stiffness and explains the influence of viscosity on the migration efficiency. In the presence of stiffness gradients and absence of chemical polarization, our simulated cell can exhibit durotaxis, anti-durotaxis, and adurotaxis respectively with increasing substrate stiffness or viscosity. The cell moves toward an optimally stiff region from softer regions during durotaxis and from stiffer regions during anti-durotaxis. We show that cell polarization through steep Rho GTPase gradients can reverse the migration direction dictated by the mechanical cues. Overall, our theory demonstrates that opposing durotactic behaviors emerge via the interplay between intracellular signaling and cell-medium mechanical interactions in agreement with experiments, thereby elucidating complex mechanosensing at the single-cell level.     

The Chromatin Extrusion Phenomenon in Amoeba proteus Cell Cycle

Amoeba proteus is possibly the best known of all unicellular eukaryotes. At the same time, several quintessential issues of its biology, including some aspects of the cell cycle, remain unsolved. Here, we show that this obligate agamic amoebae and related species have a special type of cyclic polyploidy. Their nucleus has an euploid status only for a small fraction of the cell cycle, during metaphase and telophase. The rest of the time it has an aneuploid status, which is a consequence of polyploidization. Extrusion of “excess” chromatin from the nucleus in late interphase and during prophase results in depolyploidization. Such a strategy of life cycle in unicellular eukaryotes is thought to be the main mechanism of “resetting” the Muller's ratchet and is a satisfactory alternative to the meiotic recombination for agamic protists.     

Kick-starting the ratchet: the fate of mutators in an asexual population

Muller's ratchet operates in asexual populations without intergenomic recombination. In this case, deleterious mutations will accumulate and population fitness will decline over time, possibly endangering the survival of the species. Mutator mutations, i.e., mutations that lead to an increased mutation rate, will play a special role for the behavior of the ratchet. First, they are part of the ratchet and can come to dominance through accumulation in the ratchet. Second, the fitness-loss rate of the ratchet is very sensitive to changes in the mutation rate and even a modest increase can easily set the ratchet in motion. In this article we simulate the interplay between fitness loss from Muller's ratchet and the evolution of the mutation rate from the fixation of mutator mutations. As long as the mutation rate is increased in sufficiently small steps, an accelerating ratchet and eventual extinction are inevitable. If this can be countered by antimutators, i.e., mutations that reduce the mutation rate, an equilibrium can be established for the mutation rate at some level that may allow survival. However, the presence of the ratchet amplifies fluctuations in the mutation rate and, even at equilibrium, these fluctuations can lead to dangerous bursts in the ratchet. We investigate the timescales of these processes and discuss the results with reference to the genome degradation of the aphid endosymbiont Buchnera aphidicola.