Stochastic model of integrin-mediated signaling and adhesion dynamics at the leading edges of migrating cells

Force-dependent focal adhesion assembly and disassembly: A computational study

Cells interact with the extracellular matrix (ECM) via cell-ECM adhesions. These physical interactions are transduced into biochemical signals inside the cell which influence cell behaviour. Although cell-ECM interactions have been studied extensively, it is not completely understood how immature (nascent) adhesions develop into mature (focal) adhesions and how mechanical forces influence this process. Given the small size, dynamic nature and short lifetimes of nascent adhesions, studying them using conventional microscopic and experimental techniques is challenging. Computational modelling provides a valuable resource for simulating and exploring various "what if?" scenarios in silico and identifying key molecular components and mechanisms for further investigation. Here, we present a simplified mechano-chemical model based on ordinary differential equations with three major proteins involved in adhesions: integrins, talin and vinculin. Additionally, we incorporate a hypothetical signal molecule that influences adhesion (dis)assembly rates. We find that assembly and disassembly rates need to vary dynamically to limit maturation of nascent adhesions. The model predicts biphasic variation of actin retrograde velocity and maturation fraction with substrate stiffness, with maturation fractions between 18-35%, optimal stiffness of ∼1 pN/nm, and a mechanosensitive range of 1-100 pN/nm, all corresponding to key experimental findings. Sensitivity analyses show robustness of outcomes to small changes in parameter values, allowing model tuning to reflect specific cell types and signaling cascades. The model proposes that signal-dependent disassembly rate variations play an underappreciated role in maturation fraction regulation, which should be investigated further. We also provide predictions on the changes in traction force generation under increased/decreased vinculin concentrations, complementing previous vinculin overexpression/knockout experiments in different cell types. In summary, this work proposes a model framework to robustly simulate the mechanochemical processes underlying adhesion maturation and maintenance, thereby enhancing our fundamental knowledge of cell-ECM interactions.

On multistability and constitutive relations of cell motion on fibronectin lanes

Cell motility on flat substrates exhibits coexisting steady and oscillatory morphodynamics, the biphasic adhesion-velocity relation, and the universal correlation between speed and persistence (UCSP) as simultaneous observations common to many cell types. Their universality and concurrency suggest a unifying mechanism causing all three of them. Stick-slip models for cells on one-dimensional lanes suggest multistability to arise from the nonlinear friction of retrograde flow. This study suggests a mechanical mechanism controlled by integrin signaling on the basis of a biophysical model and analysis of trajectories of MDA-MB-231 cells on fibronectin lanes, which additionally explains the constitutive relations. The experiments exhibit cells with steady or oscillatory morphodynamics and either spread or moving with spontaneous transitions between the dynamic regimes, spread and moving, and spontaneous direction reversals. Our biophysical model is based on the force balance at the protrusion edge, the noisy clutch of retrograde flow, and a response function of friction and membrane drag to integrin signaling. The theory reproduces the experimentally observed cell states, characteristics of oscillations, and state probabilities. Analysis of experiments with the biophysical model establishes a stick-slip oscillation mechanism, and explains multistability of cell states and the statistics of state transitions. It suggests protrusion competition to cause direction reversal events, the statistics of which explain the UCSP. The effect of integrin signaling on drag and friction explains the adhesion-velocity relation and cell behavior at fibronectin density steps. The dynamics of our mechanism are nonlinear flow mechanics driven by F-actin polymerization and shaped by the noisy clutch of retrograde flow friction, protrusion competition via membrane tension, and drag forces. Integrin signaling controls the parameters of the mechanical system.

Modeling cell protrusion predicts how myosin II and actin turnover affect adhesion-based signaling

Orchestration of cell migration is essential for development, tissue regeneration, and the immune response. This dynamic process integrates adhesion, signaling, and cytoskeletal subprocesses across spatial and temporal scales. In mesenchymal cells, adhesion complexes bound to extracellular matrix mediate both biochemical signal transduction and physical interaction with the F-actin cytoskeleton. Here, we present a mathematical model that offers insight into both aspects, considering spatiotemporal dynamics of nascent adhesions, active signaling molecules, mechanical clutching, actin treadmilling, and nonmuscle myosin II contractility. At the core of the model is a positive feedback loop, whereby adhesion-based signaling promotes generation of barbed ends at, and protrusion of, the cell's leading edge, which in turn promotes formation and stabilization of nascent adhesions. The model predicts a switch-like transition and optimality of membrane protrusion, determined by the balance of actin polymerization and retrograde flow, with respect to extracellular matrix density. The model, together with new experimental measurements, explains how protrusion can be modulated by mechanical effects (nonmuscle myosin II contractility and adhesive bond stiffness) and F-actin turnover.

Exploring the influence of cytosolic and membrane FAK activation on YAP/TAZ nuclear translocation

Membrane binding and unbinding dynamics play a crucial role in the biological activity of several nonintegral membrane proteins, which have to be recruited to the membrane to perform their functions. By localizing to the membrane, these proteins are able to induce downstream signal amplification in their respective signaling pathways. Here, we present a 3D computational approach using reaction-diffusion equations to investigate the relation between membrane localization of focal adhesion kinase (FAK), Ras homolog family member A (RhoA), and signal amplification of the YAP/TAZ signaling pathway. Our results show that the theoretical scenarios in which FAK is membrane bound yield robust and amplified YAP/TAZ nuclear translocation signals. Moreover, we predict that the amount of YAP/TAZ nuclear translocation increases with cell spreading, confirming the experimental findings in the literature. In summary, our in silico predictions show that when the cell membrane interaction area with the underlying substrate increases, for example, through cell spreading, this leads to more encounters between membrane-bound signaling partners and downstream signal amplification. Because membrane activation is a motif common to many signaling pathways, this study has important implications for understanding the design principles of signaling networks.

Discrete mechanical model of lamellipodial actin network implements molecular clutch mechanism and generates arcs and microspikes

Mechanical forces, actin filament turnover, and adhesion to the extracellular environment regulate lamellipodial protrusions. Computational and mathematical models at the continuum level have been used to investigate the molecular clutch mechanism, calculating the stress profile through the lamellipodium and around focal adhesions. However, the forces and deformations of individual actin filaments have not been considered while interactions between actin networks and actin bundles is not easily accounted with such methods. We develop a filament-level model of a lamellipodial actin network undergoing retrograde flow using 3D Brownian dynamics. Retrograde flow is promoted in simulations by pushing forces from the leading edge (due to actin polymerization), pulling forces (due to molecular motors), and opposed by viscous drag in cytoplasm and focal adhesions. Simulated networks have densities similar to measurements in prior electron micrographs. Connectivity between individual actin segments is maintained by permanent and dynamic crosslinkers. Remodeling of the network occurs via the addition of single actin filaments near the leading edge and via filament bond severing. We investigated how several parameters affect the stress distribution, network deformation and retrograde flow speed. The model captures the decrease in retrograde flow upon increase of focal adhesion strength. The stress profile changes from compression to extension across the leading edge, with regions of filament bending around focal adhesions. The model reproduces the observed reduction in retrograde flow speed upon exposure to cytochalasin D, which halts actin polymerization. Changes in crosslinker concentration and dynamics, as well as in the orientation pattern of newly added filaments demonstrate the model's ability to generate bundles of filaments perpendicular (actin arcs) or parallel (microspikes) to the protruding direction.

Towards understanding the messengers of extracellular space: Computational models of outside-in integrin reaction networks

The interactions between cells and their extracellular matrix (ECM) are critically important for homeostatic control of cell growth, proliferation, differentiation and apoptosis. Transmembrane integrin molecules facilitate the communication between ECM and the cell. Since the characterization of integrins in the late 1980s, there has been great advancement in understanding the function of integrins at different subcellular levels. However, the versatility in molecular pathways integrins are involved in, the high diversity in their interaction partners both outside and inside the cell as well as on the cell membrane and the short lifetime of events happening at the cell-ECM interface make it difficult to elucidate all the details regarding integrin function experimentally. To overcome the experimental challenges and advance the understanding of integrin biology, computational modeling tools have been used extensively. In this review, we summarize the computational models of integrin signaling while we explain the function of integrins at three main subcellular levels (outside the cell, cell membrane, cytosol). We also discuss how these computational modeling efforts can be helpful in other disciplines such as biomaterial design. As such, this review is a didactic modeling summary for biomaterial researchers interested in complementing their experimental work with computational tools or for seasoned computational scientists that would like to advance current in silico integrin models.

Cell-substrate mechanics guide collective cell migration through intercellular adhesion: a dynamic finite element cellular model

During the process of tissue formation and regeneration, cells migrate collectively while remaining connected through intercellular adhesions. However, the roles of cell-substrate and cell-cell mechanical interactions in regulating collective cell migration are still unclear. In this study, we employ a newly developed finite element cellular model to study collective cell migration by exploring the effects of mechanical feedback between cell and substrate and mechanical signal transmission between adjacent cells. Our viscoelastic model of cells consists many triangular elements and is of high resolution. Cadherin adhesion between cells is modeled explicitly as linear springs at subcellular level. In addition, we incorporate a mechano-chemical feedback loop between cell-substrate mechanics and Rac-mediated cell protrusion. Our model can reproduce a number of experimentally observed patterns of collective cell migration during wound healing, including cell migration persistence, separation distance between cell pairs and migration direction. Moreover, we demonstrate that cell protrusion determined by the cell-substrate mechanics plays an important role in guiding persistent and oriented collective cell migration. Furthermore, this guidance cue can be maintained and transmitted to submarginal cells of long distance through intercellular adhesions. Our study illustrates that our finite element cellular model can be employed to study broad problems of complex tissue in dynamic changes at subcellular level.

A computational study of amoeboid motility in 3D: the role of extracellular matrix geometry, cell deformability, and cell-matrix adhesion

Amoeboid cells often migrate using pseudopods, which are membrane protrusions that grow, bifurcate, and retract dynamically, resulting in a net cell displacement. Many cells within the human body, such as immune cells, epithelial cells, and even metastatic cancer cells, can migrate using the amoeboid phenotype. Amoeboid motility is a complex and multiscale process, where cell deformation, biochemistry, and cytosolic and extracellular fluid motions are coupled. Furthermore, the extracellular matrix (ECM) provides a confined, complex, and heterogeneous environment for the cells to navigate through. Amoeboid cells can migrate without significantly remodeling the ECM using weak or no adhesion, instead utilizing their deformability and the microstructure of the ECM to gain enough traction. While a large volume of work exists on cell motility on 2D substrates, amoeboid motility is 3D in nature. Despite recent progress in modeling cellular motility in 3D, there is a lack of systematic evaluations of the role of ECM microstructure, cell deformability, and adhesion on 3D motility. To fill this knowledge gap, here we present a multiscale, multiphysics modeling study of amoeboid motility through 3D-idealized ECM. The model is a coupled fluid‒structure and coarse-grain biochemistry interaction model that accounts for large deformation of cells, pseudopod dynamics, cytoplasmic and extracellular fluid motion, stochastic dynamics of cell-ECM adhesion, and microstructural (pore-scale) geometric details of the ECM. The key finding of the study is that cell deformation and matrix porosity strongly influence amoeboid motility, while weak adhesion and microscale structural details of the ECM have secondary but subtle effects.

Rac activation is key to cell motility and directionality: An experimental and modelling investigation

Cell migration is a tightly-regulated process that involves protein gradients formed by the Rho family of GTPases, including Rho and Rac. The front (rear) of cells is generally characterized by higher active Rac (Rho) and lower active Rho (Rac) concentrations. Protein clusters, called adhesions, that anchor cells to their external environment have been shown to be dynamic and small (stable and large) at the cell front (rear), forming the force-transmission points necessary for persistent movement. Differences in adhesion sizes and dynamics have been linked to gradients in Rac and Rho activity. Here, we study the effects of Rac activation and gradients in Rac and Rho concentrations and activities on cellular polarity and adhesion size using mathematical and experimental approaches. The former is accomplished by expanding an existing reaction-diffusion model to a 2D domain utilizing stochastic dynamics. The model revealed that a hysteresis between the induced/uninduced states (corresponding to higher/lower Rac concentrations, respectively) along with Rac and Rho activation gradients, generated by chemical cues, were vital for forming polarity. Experimentally, the induced state was generated by increasing the cellular βPIX (a Rac-GEF) level and/or decreasing ROCK (a Rac-GAP effector protein) activity with Y-27632 (a ROCK-inhibitor). In agreement with the simulations, our results showed that cells with elevated RacGTP migrated faster, indicating more robust cellular polarization. However, the directionality of cells was not changed significantly, suggesting that external and/or internal physical or chemical cues were needed. Complementing the faster migration observed, adhesions were smaller, generating the phenotype expected with the induced state.

Dynamics of Mechanosensitive Nascent Adhesion Formation

Cellular migration is a tightly regulated process that involves actin cytoskeleton, adaptor proteins, and integrin receptors. Forces are transmitted extracellularly through protein complexes of these molecules, called adhesions. Adhesions anchor the cell to its substrate, allowing it to migrate. In Chinese hamster ovary cells, three classes of adhesion can be identified: nascent adhesions (NAs), focal complexes, and focal adhesions, ranked here ascendingly based on size and stability. To understand the dynamics and mechanosensitive properties of NAs, a biophysical model of these NAs as colocalized clusters of integrins and adaptor proteins is developed. The model is then analyzed to characterize the dependence of NA area on biophysical parameters that regulate the number of integrins and adaptor proteins within NAs through a mechanosensitive coaggregation mechanism. Our results reveal that NA formation is triggered beyond a threshold of adaptor protein, integrin, or extracellular ligand densities, with these three factors listed in descending order of their relative influence on NA area. Further analysis of the model also reveals that an increase in coaggregation or reductions in integrin mobility inside the adhesion potentiate NA formation. By extending the model to consider the mechanosensitivity of the integrin bond, we identify mechanical stress, rather than mechanical load, as a permissive mechanical parameter that allows for noise-dependent and independent NA assembly, despite both parameters producing a bistable switch possessing a hysteresis. Stochastic simulations of the model confirm these results computationally. This study thus provides insight into the mechanical conditions defining NA dynamics.

Multivalent Binding of a Ligand-Coated Particle: Role of Shape, Size, and Ligand Heterogeneity

We utilize a multiscale modeling framework to study the effect of shape, size, and ligand composition on the efficacy of binding of a ligand-coated particle to a substrate functionalized with the target receptors. First, we show how molecular dynamics along with steered molecular dynamics calculations can be used to accurately parameterize the molecular-binding free energy and the effective spring constant for a receptor-ligand pair. We demonstrate this for two ligands that bind to the α₅β₁-domain of integrin. Next, we show how these effective potentials can be used to build computational models at the meso- and continuum-scales. These models incorporate the molecular nature of the receptor-ligand interactions and yet provide an inexpensive route to study the multivalent interaction of receptors and ligands through the construction of Bell potentials customized to the molecular identities. We quantify the binding efficacy of the ligand-coated-particle in terms of its multivalency, binding free-energy landscape, and the losses in the configurational entropies. We show that 1) the binding avidity for particle sizes less than 350 nm is set by the competition between the enthalpic and entropic contributions, whereas that for sizes above 350 nm is dominated by the enthalpy of binding; 2) anisotropic particles display higher levels of multivalent binding compared to those of spherical particles; and 3) variations in ligand composition can alter binding avidity without altering the average multivalency. The methods and results presented here have wide applications in the rational design of functionalized carriers and also in understanding cell adhesion.

Axon growth regulation by a bistable molecular switch

For the brain to function properly, its neurons must make the right connections during neural development. A key aspect of this process is the tight regulation of axon growth as axons navigate towards their targets. Neuronal growth cones at the tips of developing axons switch between growth and paused states during axonal pathfinding, and this switching behaviour determines the heterogeneous axon growth rates observed during brain development. The mechanisms controlling this switching behaviour, however, remain largely unknown. Here, using mathematical modelling, we predict that the molecular interaction network involved in axon growth can exhibit bistability, with one state representing a fast-growing growth cone state and the other a paused growth cone state. Owing to stochastic effects, even in an unchanging environment, model growth cones reversibly switch between growth and paused states. Our model further predicts that environmental signals could regulate axon growth rate by controlling the rates of switching between the two states. Our study presents a new conceptual understanding of growth cone switching behaviour, and suggests that axon guidance may be controlled by both cell-extrinsic factors and cell-intrinsic growth regulatory mechanisms.

Computational modeling of synergistic interaction between αVβ3 integrin and VEGFR2 in endothelial cells: Implications for the mechanism of action of angiogenesis-modulating integrin-binding peptides

Cooperation between VEGFR2 and integrin αVβ3 is critical for neovascularization in wound healing, cardiovascular ischemic diseases, ocular diseases, and tumor angiogenesis. In the present study, we developed a rule-based computational model to investigate the potential mechanism by which the Src-induced integrin association with VEGFR2 enhances VEGFR2 activation. Simulations demonstrated that the main function of integrin is to reduce the degradation of VEGFR2 and hence stabilize the activation signal. In addition, receptor synthesis rate and recruitment from internal compartment were found to be sensitive determinants of the activation state of VEGFR2. The model was then applied to simulate the effect of integrin-binding peptides such as tumstatin and cilengitide on VEGFR2 signaling. Further, computational modeling proposed potential molecular mechanisms for the angiogenesis-modulating activity of other integrin-binding peptides. The model highlights the complexity of the crosstalk between αVβ3 integrin and VEGFR2 and the necessity of utilizing models to elucidate potential mechanisms in angiogenesis-modulating peptide therapy.

Paxillin phosphorylation at serine 273 and its effects on Rac, Rho and adhesion dynamics

Focal adhesions are protein complexes that anchor cells to the extracellular matrix. During migration, the growth and disassembly of these structures are spatiotemporally regulated, with new adhesions forming at the leading edge of the cell and mature adhesions disassembling at the rear. Signalling proteins and structural cytoskeletal components tightly regulate adhesion dynamics. Paxillin, an adaptor protein within adhesions, is one of these proteins. Its phosphorylation at serine 273 (S273) is crucial for maintaining fast adhesion assembly and disassembly. Paxillin is known to bind to a GIT1-βPIX-PAK1 complex, which increases the local activation of the small GTPase Rac. To understand quantitatively the behaviour of this system and how it relates to adhesion assembly/disassembly, we developed a mathematical model describing the dynamics of the small GTPases Rac and Rho as determined by paxillin S273 phosphorylation. Our model revealed that the system possesses bistability, where switching between uninduced (active Rho) and induced (active Rac) states can occur through a change in rate of paxillin phosphorylation or PAK1 activation. The bistable switch is characterized by the presence of memory, minimal change in the levels of active Rac and Rho within the induced and uninduced states, respectively, and the limited regime of monostability associated with the uninduced state. These results were validated experimentally by showing the presence of bimodality in adhesion assembly and disassembly rates, and demonstrating that Rac activity increases after treating Chinese Hamster Ovary cells with okadaic acid (a paxillin phosphatase inhibitor), followed by a modest recovery after 20 min washout. Spatial gradients of phosphorylated paxillin in a reaction-diffusion model gave rise to distinct regions of Rac and Rho activities, resembling polarization of a cell into front and rear. Perturbing several parameters of the model also revealed important insights into how signalling components upstream and downstream of paxillin phosphorylation affect dynamics.

Cellular migration is crucial in both physiological and pathological functions. Maintenance of proper migration and development of aberrant migration are effectuated by cellular machinery involving protein complexes, called adhesions, that anchor the cell to its environment. Over time, these adhesions assemble at the leading edge, as the cell extends forward, anchoring the front of the cells to its substrate, while those at the cell rear disassemble, allowing detachment and forward movement. Their dynamics are controlled by a number of regulatory factors, occurring on both cell-wide and adhesion-level scales. The coordination of these regulatory factors is complex, but insights about their dynamics can be gained from the use of mathematical modeling techniques which integrate many of these components together. Here, we developed several molecularly explicit models to explore how local regulation of paxillin, an adhesion protein, interacts with the activities of Rac and Rho to produce cell-wide polarization associated with motility and directionality. By altering paxillin phosphorylation/dephosphorylation within such models, we have advanced our understanding of how a shift from a non-motile state to a highly motile state occurs. Deciphering these key processes quantitatively thus helped us gain insight into the subcellular factors underlying polarity and movement.

Cell protrusion and retraction driven by fluctuations in actin polymerization: A two-dimensional model

Animal cells that spread onto a surface often rely on actin-rich lamellipodial extensions to execute protrusion. Many cell types recently adhered on a two-dimensional substrate exhibit protrusion and retraction of their lamellipodia, even though the cell is not translating. Travelling waves of protrusion have also been observed, similar to those observed in crawling cells. These regular patterns of protrusion and retraction allow quantitative analysis for comparison to mathematical models. The periodic fluctuations in leading edge position of XTC cells have been linked to excitable actin dynamics using a one-dimensional model of actin dynamics, as a function of arc-length along the cell. In this work we extend this earlier model of actin dynamics into two dimensions (along the arc-length and radial directions of the cell) and include a model membrane that protrudes and retracts in response to the changing number of free barbed ends of actin filaments near the membrane. We show that if the polymerization rate at the barbed ends changes in response to changes in their local concentration at the leading edge and/or the opposing force from the cell membrane, the model can reproduce the patterns of membrane protrusion and retraction seen in experiment. We investigate both Brownian ratchet and switch-like force-velocity relationships between the membrane load forces and actin polymerization rate. The switch-like polymerization dynamics recover the observed patterns of protrusion and retraction as well as the fluctuations in F-actin concentration profiles. The model generates predictions for the behavior of cells after local membrane tension perturbations.

Nanoscale mechanobiology of cell adhesions

Proper physiological functions of cells and tissues depend upon their abilities to sense, transduce, integrate, and generate mechanical and biochemical signals. Although such mechanobiological phenomena are widely observed, the molecular mechanisms driving these outcomes are still not fully understood. Cell adhesions formed by integrins and cadherins receptors are key structures that process diverse sources of signals to elicit complex mechanobiological responses. Since the nanoscale is the length scale at which molecules interact to relay force and information, the understanding of cell adhesions at the nanoscale level is important for grasping the inner logics of cellular decision making. Until recently, the study of the biological nanoscale has been restricted by available molecular and imaging tools. Fortunately, rapid technological advances have increasingly opened up the nanoscale realm to systematic investigations. In this review, we discuss current insights and key open questions regarding the nanoscale structure and function relationship of cell adhesions, focusing on recent progresses in characterizing their composition, spatial organization, and cytomechanical operation.

A Computational Model of YAP/TAZ Mechanosensing

In cell proliferation, stem cell differentiation, chemoresistance, and tissue organization, the ubiquitous role of YAP/TAZ continues to impact our fundamental understanding in numerous physiological and disease systems. YAP/TAZ is an important signaling nexus integrating diverse mechanical and biochemical signals, such as ECM stiffness, adhesion ligand density, or cell-cell contacts, and thus strongly influences cell fate. Recent studies show that YAP/TAZ mechanical sensing is dependent on RhoA-regulated stress fibers. However, current understanding of YAP/TAZ remains limited due to the unknown interaction between the canonical Hippo pathway and cell tension. Furthermore, the multiscale relationship connecting adhesion signaling to YAP/TAZ activity through cytoskeleton dynamics remains poorly understood. To identify the roles of key signaling molecules in mechanical signal sensing and transduction, we present a, to our knowledge, novel computational model of the YAP/TAZ signaling pathway. This model converts extracellular-matrix mechanical properties to biochemical signals via adhesion, and integrates intracellular signaling cascades associated with cytoskeleton dynamics. We perform perturbations of molecular levels and sensitivity analyses to predict how various signaling molecules affect YAP/TAZ activity. Adhesion molecules, such as FAK, are predicted to rescue YAP/TAZ activity in soft environments via the RhoA pathway. We also found that changes of molecule concentrations result in different patterns of YAP/TAZ stiffness response. We also investigate the sensitivity of YAP/TAZ activity to ECM stiffness, and compare with that of SRF/MAL, which is another important regulator of differentiation. In addition, the model shows that the unresolved synergistic effect of YAP/TAZ activity between the mechanosensing and the Hippo pathways can be explained by the interaction of LIM-kinase and LATS. Overall, our model provides a, to our knowledge, novel platform for studying YAP/TAZ activity in the context of integrating different signaling pathways. This platform can be used to gain, to our knowledge, new fundamental insights into roles of key molecular and mechanical regulators on development, tissue engineering, or tumor progression.

Adhesion-Dependent Wave Generation in Crawling Cells

Dynamic actin networks are excitable. In migrating cells, feedback loops can amplify stochastic fluctuations in actin dynamics, often resulting in traveling waves of protrusion. The precise contributions of various molecular and mechanical interactions to wave generation have been difficult to disentangle, in part due to complex cellular morphodynamics. Here we used a relatively simple cell type-the fish epithelial keratocyte-to define a set of mechanochemical feedback loops underlying actin network excitability and wave generation. Although keratocytes are normally characterized by the persistent protrusion of a broad leading edge, increasing cell-substrate adhesion strength results in waving protrusion of a short leading edge. We show that protrusion waves are due to fluctuations in actin polymerization rates and that overexpression of VASP, an actin anti-capping protein that promotes actin polymerization, switches highly adherent keratocytes from waving to persistent protrusion. Moreover, VASP localizes both to adhesion complexes and to the leading edge. Based on these results, we developed a mathematical model for protrusion waves in which local depletion of VASP from the leading edge by adhesions-along with lateral propagation of protrusion due to the branched architecture of the actin network and negative mechanical feedback from the cell membrane-results in regular protrusion waves. Consistent with our model simulations, we show that VASP localization at the leading edge oscillates, with VASP leading-edge enrichment greatest just prior to protrusion initiation. We propose that the mechanochemical feedbacks underlying wave generation in keratocytes may constitute a general module for establishing excitable actin dynamics in other cellular contexts.

Spatial and temporal coordination of traction forces in one-dimensional cell migration

Migration of a fibroblast along a collagen fiber can be regarded as cell locomotion in one-dimension (1D). In this process, a cell protrudes forward, forms a new adhesion, produces traction forces, and releases its rear adhesion in order to advance itself along a path. However, how a cell coordinates its adhesion formation, traction forces, and rear release in 1D migration is unclear. Here, we studied fibroblasts migrating along a line of microposts. We found that when the front of a cell protruded onto a new micropost, the traction force produced at its front increased steadily, but did so without a temporal correlation in the force at its rear. Instead, the force at the front coordinated with a decrease in force at the micropost behind the front. A similar correlation in traction forces also occurred at the rear of a cell, where a decrease in force due to adhesion detachment corresponded to an increase in force at the micropost ahead of the rear. Analysis with a bio-chemo-mechanical model for traction forces and adhesion dynamics indicated that the observed relationship between traction forces at the front and back of a cell is possible only when cellular elasticity is lower than the elasticity of the cellular environment.

Linking morphodynamics and directional persistence of T lymphocyte migration

T cells play a central role in the adaptive immune response, and their directed migration is essential for homing to sites of antigen presentation. Like neutrophils, T lymphocytes are rapidly moving cells that exhibit amoeboid movement, characterized by a definitive polarity with F-actin concentrated at the front and myosin II elsewhere. In this study, we used total internal reflection fluorescence (TIRF) microscopy to monitor the cells' areas of contact with a surface presenting adhesive ICAM-1 and the chemokine, CXCL12/SDF-1. Our analysis reveals that T-cell migration and reorientation are achieved by bifurcation and lateral separation of protrusions along the leading membrane edge, followed by cessation of one of the protrusions, which acts as a pivot for cell turning. We show that the distribution of bifurcation frequencies exhibits characteristics of a random, spontaneous process; yet, the waiting time between bifurcation events depends on whether or not the pivot point remains on the same side of the migration axis. Our analysis further suggests that switching of the dominant protrusion between the two sides of the migration axis is associated with persistent migration, whereas the opposite is true of cell turning. To help explain the bifurcation phenomenon and how distinct migration behaviours might arise, a spatio-temporal, stochastic model describing F-actin dynamics is offered.

Methods of myofibrillogenesis modeling

Organization in the heart is important on multiple length scales. Myofibrillogenesis processes control the assembly of this multi-scale architecture. Understanding myofibrillogenesis might allow us to better control self-assembly of cardiac tissues. One approach consists of creating phenomenological models and comparing these models to in vitro data from primary myocytes. In this chapter, we present a method for building these models to recapitulate different aspects of myofibrillogenesis. We present a specific example for a cardiomyocyte model, but the same procedure can be used to model fibrillogenesis with other mechanisms such as motility. In sum, the models allow for a better understanding of mechanisms behind self-assembly.

Modeling the formation of cell-matrix adhesions on a single 3D matrix fiber

Cell-matrix adhesions are crucial in different biological processes like tissue morphogenesis, cell motility, and extracellular matrix remodeling. These interactions that link cell cytoskeleton and matrix fibers are built through protein clutches, generally known as adhesion complexes. The adhesion formation process has been deeply studied in two-dimensional (2D) cases; however, the knowledge is limited for three-dimensional (3D) cases. In this work, we simulate different local extracellular matrix properties in order to unravel the fundamental mechanisms that regulate the formation of cell-matrix adhesions in 3D. We aim to study the mechanical interaction of these biological structures through a three dimensional discrete approach, reproducing the transmission pattern force between the cytoskeleton and a single extracellular matrix fiber. This numerical model provides a discrete analysis of the proteins involved including spatial distribution, interaction between them, and study of the different phenomena, such as protein clutches unbinding or protein unfolding.

Model for adhesion clutch explains biphasic relationship between actin flow and traction at the cell leading edge

Cell motility relies on the continuous reorganization of a dynamic actin-myosin-adhesion network at the leading edge of the cell, in order to generate protrusion at the leading edge and traction between the cell and its external environment. We analyze experimentally measured spatial distributions of actin flow, traction force, myosin density, and adhesion density in control and pharmacologically perturbed epithelial cells in order to develop a mechanical model of the actin-adhesion-myosin self-organization at the leading edge. A model in which the F-actin network is treated as a viscous gel, and adhesion clutch engagement is strengthened by myosin but weakened by actin flow, can explain the measured molecular distributions and correctly predict the spatial distributions of the actin flow and traction stress. We test the model by comparing its predictions with measurements of the actin flow and traction stress in cells with fast and slow actin polymerization rates. The model predicts how the location of the lamellipodium-lamellum boundary depends on the actin viscosity and adhesion strength. The model further predicts that the location of the lamellipodium-lamellum boundary is not very sensitive to the level of myosin contraction.

The spatial-temporal characteristics of type I collagen-based extracellular matrix

Type I collagen abounds in mammalian extracellular matrix (ECM) and is crucial to many biophysical processes. While previous studies have mostly focused on bulk averaged properties, here we provide a comprehensive and quantitative spatial-temporal characterization of the microstructure of type I collagen-based ECM as the gelation temperature varies. The structural characteristics including the density and nematic correlation functions are obtained by analyzing confocal images of collagen gels prepared at a wide range of gelation temperatures (from 16 °C to 36 °C). As temperature increases, the gel microstructure varies from a "bundled" network with strong orientational correlation between the fibers to an isotropic homogeneous network with no significant orientational correlation, as manifested by the decaying of length scales in the correlation functions. We develop a kinetic Monte-Carlo collagen growth model to better understand how ECM microstructure depends on various environmental or kinetic factors. We show that the nucleation rate, growth rate, and an effective hydrodynamic alignment of collagen fibers fully determines the spatiotemporal fluctuations of the density and orientational order of collagen gel microstructure. Also the temperature dependence of the growth rate and nucleation rate follow the prediction of classical nucleation theory.

Formation of transient lamellipodia

Cell motility driven by actin polymerization is pivotal to the development and survival of organisms and individual cells. Motile cells plated on flat substrates form membrane protrusions called lamellipodia. The protrusions repeatedly appear and retract in all directions. If a lamellipodium is stabilized and lasts for some time, it can take over the lead and determine the direction of cell motion. Protrusions traveling along the cell perimeter have also been observed. Their initiation is in some situations the effect of the dynamics of the pathway linking plasma membrane receptors to actin filament nucleation, e.g. in chemotaxis. However, lamellipodia are also formed in many cells incessantly during motion with a constant state of the signaling pathways upstream from nucleation promoting factors (NPFs), or spontaneously in resting cells. These observations strongly suggest protrusion formation can also be a consequence of the dynamics downstream from NPFs, with signaling setting the dynamic regime but not initiating the formation of individual protrusions. A quantitative mechanism for this kind of lamellipodium dynamics has not been suggested yet. Here, we present a model exhibiting excitable actin network dynamics. Individual lamellipodia form due to random supercritical filament nucleation events amplified by autocatalytic branching. They last for about 30 seconds to many minutes and are terminated by filament bundling, severing and capping. We show the relevance of the model mechanism for experimentally observed protrusion dynamics by reproducing in very good approximation the repetitive protrusion formation measured by Burnette et al. with respect to the velocities of leading edge protrusion and retrograde flow, oscillation amplitudes, periods and shape, as well as the phase relation between protrusion and retrograde flow. Our modeling results agree with the mechanism of actin bundle formation during lamellipodium retraction suggested by Burnette et al. and Koestler et al.

Bidirectional coupling between integrin-mediated signaling and actomyosin mechanics explains matrix-dependent intermittency of leading-edge motility

A physicochemical model is used to describe the coupling of adhesion, cytoskeletal, and signaling dynamics during cell migration. Analysis of stochastic simulations predicts relationships between measurable quantities that reflect partitioning of stress between F-actin–bound adhesions, which act as a molecular clutch, and retrograde F-actin flow.

Animal cell migration is a complex process characterized by the coupling of adhesion, cytoskeletal, and signaling dynamics. Here we model local protrusion of the cell edge as a function of the load-bearing properties of integrin-based adhesions, actin polymerization fostered by adhesion-mediated signaling, and mechanosensitive activation of RhoA that promotes myosin II–generated stress on the lamellipodial F-actin network. Analysis of stochastic model simulations illustrates how these pleiotropic functions of nascent adhesions may be integrated to govern temporal persistence and frequency of protrusions. The simulations give mechanistic insight into the documented effects of extracellular matrix density and myosin abundance, and they show characteristic, nonnormal distributions of protrusion duration times that are similar to those extracted from live-cell imaging experiments. Analysis of the model further predicts relationships between measurable quantities that reflect the partitioning of stress between tension on F-actin–bound adhesions, which act as a molecular clutch, and dissipation by retrograde F-actin flow.

Mathematical modeling of eukaryotic cell migration: insights beyond experiments

A migrating cell is a molecular machine made of tens of thousands of short-lived and interacting parts. Understanding migration means understanding the self-organization of these parts into a system of functional units. This task is one of tackling complexity: First, the system integrates numerous chemical and mechanical component processes. Second, these processes are connected in feedback interactions and over a large range of spatial and temporal scales. Third, many processes are stochastic, which leads to heterogeneous migration behaviors. Early on in the research of cell migration it became evident that this complexity exceeds human intuition. Thus, the cell migration community has led the charge to build mathematical models that could integrate the diverse experimental observations and measurements in consistent frameworks, first in conceptual and more recently in molecularly explicit models. The main goal of this review is to sift through a series of important conceptual and explicit mathematical models of cell migration and to evaluate their contribution to the field in their ability to integrate critical experimental data.

Dynamic modeling of cell migration and spreading behaviors on fibronectin coated planar substrates and micropatterned geometries

An integrative cell migration model incorporating focal adhesion (FA) dynamics, cytoskeleton and nucleus remodeling, actin motor activity, and lamellipodia protrusion is developed for predicting cell spreading and migration behaviors. This work is motivated by two experimental works: (1) cell migration on 2-D substrates under various fibronectin concentrations and (2) cell spreading on 2-D micropatterned geometries. These works suggest (1) cell migration speed takes a maximum at a particular ligand density (∼1140 molecules/µm(2)) and (2) that strong traction forces at the corners of the patterns may exist due to combined effects exerted by actin stress fibers (SFs). The integrative model of this paper successfully reproduced these experimental results and indicates the mechanism of cell migration and spreading. In this paper, the mechanical structure of the cell is modeled as having two elastic membranes: an outer cell membrane and an inner nuclear membrane. The two elastic membranes are connected by SFs, which are extended from focal adhesions on the cortical surface to the nuclear membrane. In addition, the model also includes ventral SFs bridging two focal adhesions on the cell surface. The cell deforms and gains traction as transmembrane integrins distributed over the outer cell membrane bond to ligands on the ECM surface, activate SFs, and form focal adhesions. The relationship between the cell migration speed and fibronectin concentration agrees with existing experimental data for Chinese hamster ovary (CHO) cell migrations on fibronectin coated surfaces. In addition, the integrated model is validated by showing persistent high stress concentrations at sharp geometrically patterned edges. This model will be used as a predictive model to assist in design and data processing of upcoming microfluidic cell migration assays.

Traveling waves in actin dynamics and cell motility

Much of current understanding of cell motility arose from studying steady treadmilling of actin arrays. Recently, there have been a growing number of observations of a more complex, non-steady, actin behavior, including self-organized waves. It is becoming clear that these waves result from activation and inhibition feedbacks in actin dynamics acting on different scales, but the exact molecular nature of these feedbacks and the respective roles of biomechanics and biochemistry are still unclear. Here, we review recent advances achieved in experimental and theoretical studies of actin waves and discuss mechanisms and physiological significance of wavy protrusions.

Migrating fibroblasts reorient directionality by a metastable, PI3K-dependent mechanism

Mesenchymal cell migration as exhibited by fibroblasts is distinct from amoeboid cell migration and is characterized by dynamic competition among multiple protrusions, which determines directional persistence and responses to spatial cues. Localization of phosphoinositide 3-kinase (PI3K) signaling is thought to play a broadly important role in cell motility, yet the context-dependent functions of this pathway have not been adequately elucidated. By mapping the spatiotemporal dynamics of cell protrusion/retraction and PI3K signaling monitored by total internal reflection fluorescence microscopy, we show that randomly migrating fibroblasts reorient polarity through PI3K-dependent branching and pivoting of protrusions. PI3K inhibition did not affect the initiation of newly branched protrusions, nor did it prevent protrusion induced by photoactivation of Rac. Rather, PI3K signaling increased after, not before, the onset of local protrusion and was required for the lateral spreading and stabilization of nascent branches. During chemotaxis, the branch experiencing the higher chemoattractant concentration was favored, and, thus, the cell reoriented so as to align with the external gradient.

Differences in adhesion and protrusion properties correlate with differences in migration speed under EGF stimulation

Background

Cell migration plays an essential role in many biological processes, such as cancer metastasis, wound healing and immune response. Cell migration is mediated through protrusion and focal adhesion (FA) assembly, maturation and disassembly. Epidermal growth factor (EGF) is known to enhance migration rate in many cell types; however it is not known how FA maturation, FA dynamics and protrusion dynamics are regulated during EGF-induced migration. Here we use total internal reflection fluorescence (TIRF) microscopy and image analysis to quantify FA properties and protrusion dynamics under different doses of EGF stimulation.

Results

EGF was found to broaden the distribution of cell migration rates, generating more fast and slow cells. Furthermore, groups based on EGF stimulation condition or cell migration speed were marked by characteristic signatures. When data was binned based on EGF stimulation conditions, FA intensity and FA number per cell showed the largest difference among stimulation groups. FA intensity decreased with increasing EGF concentration and FA number per cell was highest under intermediate stimulation conditions. No difference in protrusion behavior was observed. However, when data was binned based on cell migration speed, FA intensity and not FA number per cell showed the largest difference among groups. FA intensity was lower for fast migrating cells. Additionally, waves of protrusion tended to correlate with fast migrating cells.

Conclusions

Only a portion of the FA properties and protrusion dynamics that correlate with migration speed, correlate with EGF stimulation condition. Those that do not correlate with EGF stimulation condition constitute the most sensitive output for identifying why cells respond differently to EGF. The idea that EGF can both increase and decrease the migration speed of individual cells in a population has particular relevance to cancer metastasis where the microenvironment can select subpopulations based on some adhesion and protrusion characteristics, leading to a more invasive phenotype as would be seen if all cells responded like an “average” cell.

Coupling actin flow, adhesion, and morphology in a computational cell motility model

Cell migration is a pervasive process in many biology systems and involves protrusive forces generated by actin polymerization, myosin dependent contractile forces, and force transmission between the cell and the substrate through adhesion sites. Here we develop a computational model for cell motion that uses the phase-field method to solve for the moving boundary with physical membrane properties. It includes a reaction-diffusion model for the actin-myosin machinery and discrete adhesion sites which can be in a "gripping" or "slipping" mode and integrates the adhesion dynamics with the dynamics of the actin filaments, modeled as a viscous network. To test this model, we apply it to fish keratocytes, fast moving cells that maintain their morphology, and show that we are able to reproduce recent experimental results on actin flow and stress patterns. Furthermore, we explore the phase diagram of cell motility by varying myosin II activity and adhesion strength. Our model suggests that the pattern of the actin flow inside the cell, the cell velocity, and the cell morphology are determined by the integration of actin polymerization, myosin contraction, adhesion forces, and membrane forces.

Excitable actin dynamics in lamellipodial protrusion and retraction

Many animal cells initiate crawling by protruding lamellipodia, consisting of a dense network of actin filaments, at their leading edge. We imaged XTC cells that exhibit flat lamellipodia on poly-L-lysine-coated coverslips. Using active contours, we tracked the leading edge and measured the total amount of F-actin by summing the pixel intensities within a 5-μm band. We observed protrusion and retraction with period 130-200 s and local wavelike features. Positive (negative) velocities correlated with minimum (maximum) integrated actin concentration. Approximately constant retrograde flow indicated that protrusions and retractions were driven by fluctuations of the actin polymerization rate. We present a model of these actin dynamics as an excitable system in which a diffusive, autocatalytic activator causes actin polymerization; F-actin accumulation in turn inhibits further activator accumulation. Simulations of the model reproduced the pattern of actin polymerization seen in experiments. To explore the model's assumption of an autocatalytic activation mechanism, we imaged cells expressing markers for both F-actin and the p21 subunit of the Arp2/3 complex. We found that integrated Arp2/3-complex concentrations spike several seconds before spikes of F-actin concentration. This suggests that the Arp2/3 complex participates in an activation mechanism that includes additional diffuse components. Response of cells to stimulation by fetal calf serum could be reproduced by the model, further supporting the proposed dynamical picture.

Numerically bridging lamellipodial and filopodial activity during cell spreading reveals a potentially novel trigger of focal adhesion maturation

We present a novel approach to modeling cell spreading, and use it to reveal a potentially central mechanism regulating focal adhesion maturation in various cell phenotypes. Actin bundles that span neighboring focal complexes at the lamellipodium-lamellum interface were assumed to be loaded by intracellular forces in proportion to bundle length. We hypothesized that the length of an actin bundle (with the corresponding accumulated force at its adhesions) may thus regulate adhesion maturation to ensure cell mechanical stability and morphological integrity. We developed a model to test this hypothesis, implementing a "top-down" approach to simplify certain cellular processes while explicitly incorporating complexity of other key subcellular mechanisms. Filopodial and lamellipodial activities were treated as modular processes with functional spatiotemporal interactions coordinated by rules regarding focal adhesion turnover and actin bundle dynamics. This theoretical framework was able to robustly predict temporal evolution of cell area and cytoskeletal organization as reported from a wide range of cell spreading experiments using micropatterned substrates. We conclude that a geometric/temporal modeling framework can capture the key functional aspects of the rapid spreading phase and resultant cytoskeletal complexity. Hence the model is used to reveal mechanistic insight into basic cell behavior essential for spreading. It demonstrates that actin bundles spanning nascent focal adhesions such that they are aligned to the leading edge may accumulate centripetal endogenous forces along their length, and could thus trigger focal adhesion maturation in a force-length dependent fashion. We suggest that this mechanism could be a central "integrating" factor that effectively coordinates force-mediated adhesion maturation at the lamellipodium-lamellum interface.

A review of models of fluctuating protrusion and retraction patterns at the leading edge of motile cells

A characteristic feature of motile cells as they undergo a change in motile behavior is the development of fluctuating exploratory motions of the leading edge, driven by actin polymerization. We review quantitative models of these protrusion and retraction phenomena. Theoretical studies have been motivated by advances in experimental and computational methods that allow controlled perturbations, single molecule imaging, and analysis of spatiotemporal correlations in microscopic images. To explain oscillations and waves of the leading edge, most theoretical models propose nonlinear interactions and feedback mechanisms among different components of the actin cytoskeleton system. These mechanisms include curvature-sensing membrane proteins, myosin contraction, and autocatalytic biochemical reaction kinetics. We discuss how the combination of experimental studies with modeling promises to quantify the relative importance of these biochemical and biophysical processes at the leading edge and to evaluate their generality across cell types and extracellular environments.

Stochastic models of cell protrusion arising from spatiotemporal signaling and adhesion dynamics

During cell migration, local protrusion events are regulated by biochemical and physical processes that are in turn coordinated with the dynamic properties of cell-substratum adhesion structures. In this chapter, we present a modeling approach for integrating the apparent stochasticity and spatial dependence of signal transduction pathways that promote protrusion in tandem with adhesion dynamics. We describe our modeling framework, as well as its abstraction, parameterization, and validation against experimental data. Analytical techniques for identifying and evaluating the effects of model bistability on simulation simulation results are shown, and implications of this analysis for understanding cell protrusion behavior are offered.

Redundant mechanisms for stable cell locomotion revealed by minimal models

Crawling of eukaryotic cells on flat surfaces is underlain by the protrusion of the actin network, the contractile activity of myosin II motors, and graded adhesion to the substrate regulated by complex biochemical networks. Some crawling cells, such as fish keratocytes, maintain a roughly constant shape and velocity. Here we use moving-boundary simulations to explore four different minimal mechanisms for cell locomotion: 1), a biophysical model for myosin contraction-driven motility; 2), a G-actin transport-limited motility model; 3), a simple model for Rac/Rho-regulated motility; and 4), a model that assumes that microtubule-based transport of vesicles to the leading edge limits the rate of protrusion. We show that all of these models, alone or in combination, are sufficient to produce half-moon steady shapes and movements that are characteristic of keratocytes, suggesting that these mechanisms may serve redundant and complementary roles in driving cell motility. Moving-boundary simulations demonstrate local and global stability of the motile cell shapes and make testable predictions regarding the dependence of shape and speed on mechanical and biochemical parameters. The models shed light on the roles of membrane-mediated area conservation and the coupling of mechanical and biochemical mechanisms in stabilizing motile cells.

Pak1 regulates focal adhesion strength, myosin IIA distribution, and actin dynamics to optimize cell migration

p21-activated kinases are essential for spatial and temporal coordination of cytoskeletal dynamics with cellular adhesion during cell migration.

Cell motility requires the spatial and temporal coordination of forces in the actomyosin cytoskeleton with extracellular adhesion. The biochemical mechanism that coordinates filamentous actin (F-actin) assembly, myosin contractility, adhesion dynamics, and motility to maintain the balance between adhesion and contraction remains unknown. In this paper, we show that p21-activated kinases (Paks), downstream effectors of the small guanosine triphosphatases Rac and Cdc42, biochemically couple leading-edge actin dynamics to focal adhesion (FA) dynamics. Quantitative live cell microscopy assays revealed that the inhibition of Paks abolished F-actin flow in the lamella, displaced myosin IIA from the cell edge, and decreased FA turnover. We show that, by controlling the dynamics of these three systems, Paks regulate the protrusive activity and migration of epithelial cells. Furthermore, we found that expressing Pak1 was sufficient to overcome the inhibitory effects of excess adhesion strength on cell motility. These findings establish Paks as critical molecules coordinating cytoskeletal systems for efficient cell migration.

Signaling pathways that control cell migration: models and analysis

Dissecting the intracellular signaling mechanisms that govern the movement of eukaryotic cells presents a major challenge, not only because of the large number of molecular players involved, but even more so because of the dynamic nature of their regulation by both biochemical and mechanical interactions. Computational modeling and analysis have emerged as useful tools for understanding how the physical properties of cells and their microenvironment are coupled with certain biochemical pathways to actuate and control cell motility. In this focused review, we highlight some of the more recent applications of quantitative modeling and analysis in the field of cell migration. Both in modeling and experiment, it has been prudent to follow a reductionist approach in order to characterize what are arguably the principal modules: spatial polarization of signaling pathways, regulation of the actin cytoskeleton, and dynamics of focal adhesions. While it is important that we 'cut our teeth' on these subsystems, focusing on the details of certain aspects while ignoring or coarse-graining others, it is clear that the challenge ahead will be to characterize the couplings between them in an integrated framework.

An adhesion-dependent switch between mechanisms that determine motile cell shape

Keratocytes are fast-moving cells in which adhesion dynamics are tightly coupled to the actin polymerization motor that drives migration, resulting in highly coordinated cell movement. We have found that modifying the adhesive properties of the underlying substrate has a dramatic effect on keratocyte morphology. Cells crawling at intermediate adhesion strengths resembled stereotypical keratocytes, characterized by a broad, fan-shaped lamellipodium, clearly defined leading and trailing edges, and persistent rates of protrusion and retraction. Cells at low adhesion strength were small and round with highly variable protrusion and retraction rates, and cells at high adhesion strength were large and asymmetrical and, strikingly, exhibited traveling waves of protrusion. To elucidate the mechanisms by which adhesion strength determines cell behavior, we examined the organization of adhesions, myosin II, and the actin network in keratocytes migrating on substrates with different adhesion strengths. On the whole, our results are consistent with a quantitative physical model in which keratocyte shape and migratory behavior emerge from the self-organization of actin, adhesions, and myosin, and quantitative changes in either adhesion strength or myosin contraction can switch keratocytes among qualitatively distinct migration regimes.

Cell migration is important for many biological processes: white blood cells chase down and kill bacteria to guard against infection, epithelial cells crawl across open wounds to promote healing, and embryonic cells move collectively to form organs and tissues during embryogenesis. In all of these cases, migration depends on the spatial and temporal organization of multiple forces, including actin-driven protrusion of the cell membrane, membrane tension, cell-substrate adhesion, and myosin-mediated contraction of the actin network. In this work, we have used a simple cell type, the fish epithelial keratocyte, as a model system to investigate the manner in which these forces are integrated to give rise to large-scale emergent properties such as cell shape and movement. Keratocytes are normally fan-shaped and fast-moving, but we have found that keratocytes migrate more slowly and assume round or asymmetric shapes when cell-substrate adhesion strength is too high or too low. By correlating measurements of adhesion-dependent changes in cell shape and speed with measurements of adhesion and myosin localization patterns and actin network organization, we have developed a mechanical model in which keratocyte shape and movement emerge from adhesion and myosin-dependent regulation of the dynamic actin cytoskeleton.

Biomedical engineering strategies in system design space

Modern systems biology and synthetic bioengineering face two major challenges in relating properties of the genetic components of a natural or engineered system to its integrated behavior. The first is the fundamental unsolved problem of relating the digital representation of the genotype to the analog representation of the parameters for the molecular components. For example, knowing the DNA sequence does not allow one to determine the kinetic parameters of an enzyme. The second is the fundamental unsolved problem of relating the parameters of the components and the environment to the phenotype of the global system. For example, knowing the parameters does not tell one how many qualitatively distinct phenotypes are in the organism's repertoire or the relative fitness of the phenotypes in different environments. These also are challenges for biomedical engineers as they attempt to develop therapeutic strategies to treat pathology or to redirect normal cellular functions for biotechnological purposes. In this article, the second of these fundamental challenges will be addressed, and the notion of a "system design space" for relating the parameter space of components to the phenotype space of bioengineering systems will be focused upon. First, the concept of a system design space will be motivated by introducing one of its key components from an intuitive perspective. Second, a simple linear example will be used to illustrate a generic method for constructing the design space in which qualitatively distinct phenotypes can be identified and counted, their fitness analyzed and compared, and their tolerance to change measured. Third, two examples of nonlinear systems from different areas of biomedical engineering will be presented. Finally, after giving reference to a few other applications that have made use of the system design space approach to reveal important design principles, some concluding remarks concerning challenges and opportunities for further development will be made.

Computational investigation of intrinsic and extrinsic mechanisms underlying the formation of carcinoma

The disruption of tissue epithelial architecture is thought to be involved in the initiation of cancer; however, little is known about the cell biology of early neoplastic lesions. Computational models can facilitate in the study of early cancer stagesby simulating different scenarios of tumour inception and by testing in silico the influence of various intrinsic and extrinsic factors on the emerging multicellular morphologies. Here, we have used a computational model to create 3D morphocharts, collections of multicellular morphologies arising from systematic modification of three model parameters defining cell sensitivity to external cues for cell growth, death and adhesion to the extracellular matrix. This systematic search revealed the ranges of parameter values for which robust epithelial structures formed and those that led to abnormal geometrical forms resembling tumour-like morphologies. We also showed how our model can be used to map morphologies of experimental multicellular systems onto their in silico counterpart via the cell sensitivity parameter space that may in turn allow us to identify genetic or epigenetic changes responsible for these morphological distortions.

Stochastic Dynamics of Membrane Protrusion Mediated by the DOCK180/Rac Pathway in Migrating Cells

Cell migration is regulated by processes that control adhesion to extracellular matrix (ECM) and force generation. While our fundamental understanding of how these control mechanisms are actuated at the molecular level (signal transduction) has been refined over many years, appreciation of their dynamics has grown more recently. Here, we formulate and analyze by stochastic simulation a quantitative model of signaling mediated by the integrin family of adhesion receptors. Nascent adhesions foster the activation of the small GTPase Rac by at least two distinct signaling pathways, one of which involves tyrosine phosphorylation of the scaffold protein paxillin and formation of multiprotein complexes containing the guanine nucleotide exchange factor DOCK180. Active Rac promotes protrusion of the cell's leading edge, which in turn enhances the rate of nascent adhesion nucleation; we call this feedback mechanism the core protrusion cycle. Protrusion is antagonized by stable adhesions, which form by myosin-dependent maturation of nascent adhesions, and we propose here a feedforward mechanism mediated by the tyrosine kinase c-Src by which this antagonism is regulated so as to allow transient protrusion at higher densities of ECM. We show that this "buffering of inhibition" mechanism, coupled with the core protrusion cycle, is capable of tuning the frequencies of protrusion and adhesion maturation events.