On the significance of microtubule flexural behavior in cytoskeletal mechanics

Analysis of the Adherent Cell Response to the Substrate Stiffness Using Tensegrity

Cell's shape is dependent on the cytoskeleton mechanical properties. Hybrid models were developed that combine the discrete structure for the cytoskeleton and continuum parts for other cell organelles. Tensegrity-based structures that consist of tensile and compression elements are useful models to understand the cytoskeleton mechanical behavior. In this study, we are looking to examine the reaction of the cell to a variety of substrate stiffnesses and explain the relationship between cell behavior and substrate mechanical properties. However, which tensegrity structure is appropriate for modeling a living cell? Is the structure's complexity play a major role? We used two spherical tensegrities with different complexities to assess the impact of the structure on the cell's mechanical response versus substrate's stiffness. Six- and twelve-strut tensegrities together with membrane, cytoplasm, nucleoskeleton, and nucleus envelope were assembled in Abaqus package to create a hybrid cell model. A compressive load was applied to the cell model and the reaction forces versus deflection curves were analyzed for number of substrate stiffness values. By analyzing the difference due to two different tensegrities it became clear that the lower density structure is a better choice for modeling stiffer cells. It was also found that the six-strut tensegrity is sensitive to higher range of substrate stiffness.

Tubulin Post-translational Modifications: Potential Therapeutic Approaches to Heart Failure

In recent decades, advancing insights into the mechanisms of cardiac dysfunction have focused on the involvement of microtubule network. A variety of tubulin post-translational modifications have been discovered to fine-tune the microtubules’ properties and functions. Given the limits of therapies based on conserved structures of the skeleton, targeting tubulin modifications appears to be a potentially promising therapeutic strategy. Here we review the current understanding of tubulin post-translational modifications in regulating microtubule functions in the cardiac system. We also discussed how altered modifications may lead to a range of cardiac dysfunctions, many of which are linked to heart failure.

Biomechanics of cells and subcellular components: A comprehensive review of computational models and applications

Cells are a fundamental structural, functional and biological unit for all living organisms. Up till now, considerable efforts have been made to study the responses of single cells and subcellular components to an external load, and understand the biophysics underlying cell rheology, mechanotransduction and cell functions using experimental and in silico approaches. In the last decade, computational simulation has become increasingly attractive due to its critical role in interpreting experimental data, analysing complex cellular/subcellular structures, facilitating diagnostic designs and therapeutic techniques, and developing biomimetic materials. Despite the significant progress, developing comprehensive and accurate models of living cells remains a grand challenge in the 21st century. To understand current state of the art, this review summarises and classifies the vast array of computational biomechanical models for cells. The article covers the cellular components at multi-spatial levels, that is, protein polymers, subcellular components, whole cells and the systems with scale beyond a cell. In addition to the comprehensive review of the topic, this article also provides new insights into the future prospects of developing integrated, active and high-fidelity cell models that are multiscale, multi-physics and multi-disciplinary in nature. This review will be beneficial for the researchers in modelling the biomechanics of subcellular components, cells and multiple cell systems and understanding the cell functions and biological processes from the perspective of cell mechanics.

Geometrical nonlinear elasticity of axon under tension: A coarse-grained computational study

Axon bundles cross-linked by microtubule (MT) associate proteins and bounded by a shell skeleton are critical for normal function of neurons. Understanding effects of the complexly geometrical parameters on their mechanical properties can help gain a biomechanical perspective on the neurological functions of axons and thus brain disorders caused by the structural failure of axons. Here, the tensile mechanical properties of MT bundles cross-linked by tau proteins are investigated by systematically tuning MT length, axonal cross-section radius, and tau protein spacing in a bead-spring coarse-grained model. Our results indicate that the stress-strain curves of axons can be divided into two regimes, a nonlinear elastic regime dominated by rigid-body like inter-MT sliding, and a linear elastic regime dominated by affine deformation of both tau proteins and MTs. From the energetic analyses, first, the tau proteins dominate the mechanical performance of axons under tension. In the nonlinear regime, tau proteins undergo a rigid-body like rotating motion rather than elongating, whereas in the nonlinear elastic regime, tau proteins undergo a flexible elongating deformation along the MT axis. Second, as the average spacing between adjacent tau proteins along the MT axial direction increases from 25 to 125 nm, the Young's modulus of axon experiences a linear decrease whereas with the average space varying from 125 to 175 nm, and later reaches a plateau value with a stable fluctuation. Third, the increment of the cross-section radius of the MT bundle leads to a decrease in Young's modulus of axon, which is possibly attributed to the decrease in MT numbers per cross section. Overall, our research findings offer a new perspective into understanding the effects of geometrical parameters on the mechanics of MT bundles as well as serving as a theoretical basis for the development of artificial MT complexes potentially toward medical applications.

Finite Element Simulations of Mechanical Behaviour of Endothelial Cells

Biomechanical models based on the finite element method have already shown their potential in the simulation of the mechanical behaviour of cells. For instance, development of atherosclerosis is accelerated by damage of the endothelium, a monolayer of endothelial cells on the inner surface of arteries. Finite element models enable us to investigate mechanical factors not only at the level of the arterial wall but also at the level of individual cells. To achieve this, several finite element models of endothelial cells with different shapes are presented in this paper. Implementing the recently proposed bendotensegrity concept, these models consider the flexural behaviour of microtubules and incorporate also waviness of intermediate filaments. The suspended and adherent cell models are validated by comparison of their simulated force-deformation curves with experiments from the literature. The flat and dome cell models, mimicking natural cell shapes inside the endothelial layer, are then used to simulate their response in compression and shear which represent typical loads in a vascular wall. The models enable us to analyse the role of individual cytoskeletal components in the mechanical responses, as well as to quantify the nucleus deformation which is hypothesized to be the quantity decisive for mechanotransduction.

Migration of the 3T3 Cell with a Lamellipodium on Various Stiffness Substrates—Tensegrity Model

Changes in mechanical stimuli and the physiological environment are sensed by the cell. Thesechanges influence the cell’s motility patterns. The cell’s directional migration is dependent on the substrate stiffness. To describe such behavior of a cell, a tensegrity model was used. Cells with an extended lamellipodium were modeled. The internal elastic strain energy of a cell attached to the substrates with different stiffnesses was evaluated. The obtained results show that on the stiffer substrate, the elastic strain energy of the cell adherent to this substrate decreases. Therefore, the substrate stiffness is one of the parameters that govern the cell’s directional movement.

A Finite Element Bendo-Tensegrity Model of Eukaryotic Cell

Mechanical interaction of cell with extracellular environment affects its function. The mechanisms by which mechanical stimuli are sensed and transduced into biochemical responses are still not well understood. Considering this, two finite element (FE) bendo-tensegrity models of a cell in different states are proposed with the aim to characterize cell deformation under different mechanical loading conditions: a suspended cell model elucidating the global response of cell in tensile test simulation and an adherent cell model explicating its local response in atomic force microscopy (AFM) indentation simulation. The force-elongation curve obtained from tensile test simulation lies within the range of experimentally obtained characteristics of smooth muscle cells (SMCs) and illustrates a nonlinear increase in reaction force with cell stretching. The force-indentation curves obtained from indentation simulations lie within the range of experimentally obtained curves of embryonic stem cells (ESCs) and exhibit the influence of indentation site on the overall reaction force of cell. Simulation results have demonstrated that actin filaments (AFs) and microtubules (MTs) play a crucial role in the cell stiffness during stretching, whereas actin cortex (AC) along with actin bundles (ABs) and MTs are essential for the cell rigidity during indentation. The proposed models quantify the mechanical contribution of individual cytoskeletal components to cell mechanics and the deformation of nucleus under different mechanical loading conditions. These results can aid in better understanding of structure-function relationships in living cells.

Capturing 3D large-strain Euler-bending filament dynamics in fibrous media simulations; sample case of compression collapse in dendritic actin network

Cytoskeletal networks to transmission towers are comprised of slender elements. Slender filaments bend and buckle more easily than stretch. Therefore a deforming network is expected to exhaust all possible bending-based modes before engaging filament stretch. While the large-strain bending critically determines fibrous-media response, simulations use small-strain and jointed approximations. At low resolution, these approximations inflate bending resistance and delay buckling onset. The proposed string-of-continuous-beams (SOCB) approach captures 3D nonlinear Euler bending of filaments with high fidelity at low cost. Bending geometry (i.e. angles and its differentials) is solved as primary variables, to fit a 5th order polynomial of the contour angle. Displacement, solved simultaneously as length conservation, is predicted with C3 and C6 smoothness between and within segments, using only 2 nodes. In the chosen analysis frame, in-plane and out-plane moments can be decoupled for arbitrarily-curved segments. Complex crosslink force-transfers can be specified. Simulations show that when a daughter branch is appended, the buckling resistance of a filament changes from linear to nonlinear before reversible collapse. An actin outcrop with 8 generations of mother-daughter branching produced the linear, nonlinear, and collapse regimes observed in compression experiments. 'Collapse' was a redistribution of outcrop forces following the buckling of few strands.

Effects of the cross-linkers on the buckling of microtubules in cells

In cells, the protein cross-linkers lead to a distinct buckling behavior of microtubules (MTs) different from the buckling of individual MTs. This paper thus aims to examine this issue via the molecular structural mechanics (MSM) simulations. The transition of buckling responses was captured as the two-dimensional-linkers were replaced by the three-dimensional (3D) ones. Then, the effects of the radial orientation and the axial density of the 3D-linkers were examined, showing that more uniform distribution of the radial orientation leads to the higher critical load with 3D buckling modes, while the inhomogeneity of the axial density results in the localized buckling patterns. The results demonstrated the important role of the cross-linker in regulating MT stiffness, revealed the physics of the experimentally observed localized buckling and these results will pave the way to a new multi-component mechanics model for whole cells.

Looking "Under the Hood" of Cellular Mechanotransduction with Computational Tools: A Systems Biomechanics Approach across Multiple Scales

Signal modulation has been developed in living cells throughout evolution to promote utilizing the same machinery for multiple cellular functions. Chemical and mechanical modules of signal transmission and transduction are interconnected and necessary for organ development and growth. However, due to the high complexity of the intercommunication of physical intracellular connections with biochemical pathways, there are many missing details in our overall understanding of mechanotransduction processes, i.e., the process by which mechanical signals are converted to biochemical cascades. Cell-matrix adhesions are mechanically coupled to the nucleus through the cytoskeleton. This modulated and tightly integrated network mediates the transmission of mechanochemical signals from the extracellular matrix to the nucleus. Various experimental and computational techniques have been utilized to understand the basic mechanisms of mechanotransduction, yet many aspects have remained elusive. Recently, in silico experiments have made important contributions to the field of mechanobiology. Herein, computational modeling efforts devoted to understanding integrin-mediated mechanotransduction pathways are reviewed, and an outlook is presented for future directions toward using suitable computational approaches and developing novel techniques for addressing important questions in the field of mechanotransduction.

Detyrosinated microtubules buckle and bear load in contracting cardiomyocytes

The microtubule (MT) cytoskeleton can transmit mechanical signals and resist compression in contracting cardiomyocytes. How MTs perform these roles remains unclear because of difficulties in observing MTs during the rapid contractile cycle. Here, we used high spatial and temporal resolution imaging to characterize MT behavior in beating mouse myocytes. MTs deformed under contractile load into sinusoidal buckles, a behavior dependent on posttranslational "detyrosination" of α-tubulin. Detyrosinated MTs associated with desmin at force-generating sarcomeres. When detyrosination was reduced, MTs uncoupled from sarcomeres and buckled less during contraction, which allowed sarcomeres to shorten and stretch with less resistance. Conversely, increased detyrosination promoted MT buckling, stiffened the myocyte, and correlated with impaired function in cardiomyopathy. Thus, detyrosinated MTs represent tunable, compression-resistant elements that may impair cardiac function in disease.

Buckling behavior of individual and bundled microtubules

As the major structural constituent of the cytoskeleton, microtubules (MTs) serve a variety of biological functions that range from facilitating organelle transport to maintaining the mechanical integrity of the cell. Neuronal MTs exhibit a distinct configuration, hexagonally packed bundles of MT filaments, interconnected by MT-associated protein (MAP) tau. Building on our previous work on mechanical response of axonal MT bundles under uniaxial tension, this study is focused on exploring the compression scenarios. Intracellular MTs carry a large fraction of the compressive loads sensed by the cell and therefore, like any other column-like structure, are prone to substantial bending and buckling. Various biological activities, e.g., actomyosin contractility and many pathological conditions are driven or followed by bending, looping, and buckling of MT filaments. The coarse-grained model previously developed in our lab has been used to study the mechanical behavior of individual and bundled in vivo MT filaments under uniaxial compression. Both configurations show tip-localized, decaying, and short-wavelength buckling. This behavior highlights the role of the surrounding cytoplasm and MAP tau on MT buckling behavior, which allows MT filaments to bear much larger compressive forces. It is observed that MAP tau interconnections improve this effect by a factor of two. The enhanced ability of MT bundles to damp buckling waves relative to individual MT filaments, may be interpreted as a self-defense mechanism because it helps axonal MTs to endure harsher environments while maintaining their function. The results indicate that MT filaments in a bundle do not buckle simultaneously implying that the applied stress is not equally shared among the MT filaments, that is a consequence of the nonuniform distribution of MAP tau proteins along the bundle length. Furthermore, from a pathological perspective, it is observed that axonal MT bundles are more vulnerable to failure in compression than tension.

Molecular control of stress transmission in the microtubule cytoskeleton

In this article, we will summarize recent progress in understanding the mechanical origins of rigidity, strength, resiliency and stress transmission in the MT cytoskeleton using reconstituted networks formed from purified components. We focus on the role of network architecture, crosslinker compliance and dynamics, and molecular determinants of single filament elasticity, while highlighting open questions and future directions for this work.

The talin dimer structure orientation is mechanically regulated

Formation of a stable cell-substrate contact can be regulated by mechanical force, especially at the focal adhesion. Individual proteins that make up the focal adhesions, such as talin, can exhibit mechanosensing. We previously described one mode of talin mechanosensing in which the vinculin-binding site of talin is exposed after force-induced stretch of a single talin rod domain. Here, we describe a second mode of talin mechanosensing in which the talin dimer itself can adopt different orientations in response to mechanical stimulation. Using molecular dynamics models, we demonstrate that the C-terminus region of the talin dimer is flexible mainly at the linker between the dimerization helices and the nearby actin-binding helical bundle. Our molecular dynamics simulations reveal two possible orientations of the talin dimer at its C-terminus. The extracellular matrix (ECM)-bound integrins cross-linked by talin can be forced apart leading to an elongated orientation of the talin dimer, and the ECM-bound integrins can be forced together by the ECM producing a collapsed orientation of the talin dimer. Formation of the elongated orientation is shown to be more favorable. Switching between the two talin dimer orientations constitutes a mode of mechanosensing.

(+)-Naloxone inhibits morphine-induced chemotaxis via prevention of heat shock protein 90 cleavage in microglia

BACKGROUND/PURPOSE

Microglia have a crucial role in maintaining neuronal homeostasis in the central nervous system. Immune factors released from microglia have important roles in nociceptive signal transduction. Activation of microglia seems to be a shared mechanism in pathological pain and morphine tolerance because pharmacological attenuation of microglia activation provides satisfactory management in both situations.

METHODS

In the present study, we investigated the effect of 1nM (+)-naloxone, which is not an opioid receptor antagonist, on morphine-induced activation of microglia EOC13.31 cells.

RESULTS

Our results showed that 1μM morphine enhanced microglia activation and migration, decreased α-tubulin acetylation, and induced heat shock protein 90 (HSP90) fragmentation and histone deacetylase 6 (HDAC6) expression. Morphine-induced α-tubulin deacetylation and HSP90 fragmentation were HDAC6-dependent. Pretreatment with (+)-naloxone (1nM) inhibited morphine-evoked microglia activation and chemotaxis and prevented α-tubulin deacetylation and HSP90 fragmentation by inhibiting HDAC6 expression.

CONCLUSION

Based on the findings of the present study, we suggest that (+)-naloxone inhibits morphine-induced microglia activation by regulating HDAC6-dependent α-tubulin deacetylation and HSP90 fragmentation.

Localized vibration of a microtubule surrounded by randomly distributed cross linkers

Based on finite element simulation, the present work studies free vibration of a microtubule surrounded by 3D randomly distributed cross linkers in living cells. A basic result of the present work is that transverse vibration modes associated with the lowest frequencies are highly localized, in sharp contrast to the through-length modes predicted by the commonly used classic elastic foundation model. Our simulations show that the deflected length of localized modes increases with increasing frequency and approaches the entire length of microtubule when frequency approaches the minimum classic frequency given by the elastic foundation model. In particular, unlike the length-sensitive classic frequencies predicted by the elastic foundation model, the lowest frequencies of localized modes predicted by the present model are insensitive to the length of microtubules and are at least 50% lower than the minimum classic frequency for infinitely long microtubules and could be one order of magnitude lower than the minimum classic frequency for shorter microtubules (only a few microns in length). These results suggest that the existing elastic foundation model may have overestimated the lowest frequencies of microtubules in vivo. Finally, based on our simulation results, some empirical relations are proposed for the critical (lowest) frequency of localized modes and the associated wave length. Compared to the classic elastic foundation model, the localized vibration modes and the associated wave lengths predicted by the present model are in better agreement with some known experimental observations.

Effects of an antimitotic drug on mechanical behaviours of the cytoskeleton in distinct grades of colon cancer cells

Biomechanical behaviours of cells change during cancer progression due to alterations in the main cytoskeletal proteins. Microtubules play a vital role in mitosis and in supporting the integrity of the cell due to their ability to withstand high compressive loads. Accordingly, microtubule-targeting agents (MTAs) have become one of the most promising classes of drugs in cancer therapy. This study evaluated changes in visco-elastic parameters induced by an appropriate concentration of an antimitotic drug in two different grades of colon cancer cells. Actin microfilaments and microtubules contents in the cells were evaluated by Western blot analysis and fluorescence intensity calculation. Micropipette aspiration experiments showed that the MTA had distinct mechanical effects on different cell lines. The more aggressive the cells, the greater the reduction in elasticity and viscosity. Invasive cells had a higher initial instantaneous Young's modulus than primary cells, but this reduced to approximately one half of the values for primary cells after 48 h of drug treatment. A considerable association was seen between the changes in mechanical properties and the microtubule to F-actin microfilament content ratio, which decreased with MTA treatment.

Mechanics of biological networks: from the cell cytoskeleton to connective tissue

From the cell cytoskeleton to connective tissues, fibrous networks are ubiquitous in metazoan life as the key promoters of mechanical strength, support and integrity. In recent decades, the application of physics to biological systems has made substantial strides in elucidating the striking mechanical phenomena observed in such networks, explaining strain stiffening, power law rheology and cytoskeletal fluidisation - all key to the biological function of individual cells and tissues. In this review we focus on the current progress in the field, with a primer into the basic physics of individual filaments and the networks they form. This is followed by a discussion of biological networks in the context of a broad spread of recent in vitro and in vivo experiments.

Localized buckling of a microtubule surrounded by randomly distributed cross linkers

Microtubules supported by surrounding cross linkers in eukaryotic cells can bear a much higher compressive force than free-standing microtubules. Different from some previous studies, which treated the surroundings as a continuum elastic foundation or elastic medium, the present paper develops a micromechanics numerical model to examine the role of randomly distributed discrete cross linkers in the buckling of compressed microtubules. First, the proposed numerical approach is validated by reproducing the uniform multiwave buckling mode predicted by the existing elastic-foundation model. For more realistic buckling of microtubules surrounded by randomly distributed cross linkers, the present numerical model predicts that the buckling mode is localized at one end in agreement with some known experimental observations. In particular, the critical force for localized buckling, predicted by the present model, is insensitive to microtubule length and can be about 1 order of magnitude lower than those given by the elastic-foundation model, which suggests that the elastic-foundation model may have overestimated the critical force for buckling of microtubules in vivo. In addition, unlike the elastic-foundation model, the present model can capture the effect of end conditions on the critical force and wavelength of localized buckling. Based on the known data of spacing and elastic constants of cross linkers available in literature, the critical force and wavelength of the localized buckling mode, predicted by the present model, are compared to some experimental data with reasonable agreement. Finally, two empirical formulas are proposed for the critical force and wavelength of the localized buckling of microtubules surrounded by cross linkers.

Measuring microtubule persistence length using a microtubule gliding assay

The mechanical properties of microtubules have been an area of active research for the past two decades, in part because understanding the mechanics of individual microtubules contributes to modeling whole-cell rigidity and structure and hence to better understanding the processes underlying motility and transport. Moreover, the role of microtubule structure and microtubule-associated proteins (MAPs) in microtubule stiffness remains unclear. In this chapter, we present a kinesin-driven microtubule gliding assay analysis of persistence length that is amenable to simultaneous variation of microtubule parameters such as length, structure, or MAP coverage and determination of persistence length. By combining sparse fluorescent labeling of individual microtubules with single particle tracking of individual fluorophores, microtubule gliding trajectories are tracked with nanometer-level precision. The fluctuations in these trajectories, due to thermal fluctuations in the microtubules themselves, are analyzed to extract the microtubule persistence length. In the following, we describe this gliding assay and analysis and discuss two example microtubule variables, length and diameter, in anticipation that the method may be of wide use for in vitro study of microtubule mechanical properties.

Flexural rigidity measurements of biopolymers using gliding assays

Microtubules are cytoskeletal polymers which play a role in cell division, cell mechanics, and intracellular transport. Each of these functions requires microtubules that are stiff and straight enough to span a significant fraction of the cell diameter. As a result, the microtubule persistence length, a measure of stiffness, has been actively studied for the past two decades(1). Nonetheless, open questions remain: short microtubules are 10-50 times less stiff than long microtubules(2-4), and even long microtubules have measured persistence lengths which vary by an order of magnitude(5-9). Here, we present a method to measure microtubule persistence length. The method is based on a kinesin-driven microtubule gliding assay(10). By combining sparse fluorescent labeling of individual microtubules with single particle tracking of individual fluorophores attached to the microtubule, the gliding trajectories of single microtubules are tracked with nanometer-level precision. The persistence length of the trajectories is the same as the persistence length of the microtubule under the conditions used(11). An automated tracking routine is used to create microtubule trajectories from fluorophores attached to individual microtubules, and the persistence length of this trajectory is calculated using routines written in IDL. This technique is rapidly implementable, and capable of measuring the persistence length of 100 microtubules in one day of experimentation. The method can be extended to measure persistence length under a variety of conditions, including persistence length as a function of length along microtubules. Moreover, the analysis routines used can be extended to myosin-based acting gliding assays, to measure the persistence length of actin filaments as well.

Computational model for the cell-mechanical response of the osteocyte cytoskeleton based on self-stabilizing tensegrity structures

The mechanism by which mechanical stimulation on osteocytes results in biochemical signals that initiate the remodeling process inside living bone tissue is largely unknown. Even the type of stimulation acting on these cells is not yet clearly identified. However, the cytoskeleton of osteocytes is suggested to play a major role in the mechanosensory process due to the direct connection to the nucleus. In this paper, a computational approach to model and simulate the cell structure of osteocytes based on self-stabilizing tensegrity structures is suggested. The computational model of the cell consists of the major components with respect to mechanical aspects: the integrins that connect the cell with the extracellular bone matrix, and different types of protein fibers (microtubules and intermediate filaments) that form the cytoskeleton, the membrane-cytoskeleton (microfilaments), the nucleus and the centrosome. The proposed geometrical cell models represent the cell in its physiological environment which is necessary in order to give a statement on the cell behavior in vivo. Studies on the mechanical response of osteocytes after physiological loading and in particular the mechanical response of the nucleus show that the load acting on the nucleus is rising with increasing deformation applied to the integrins.

Computational modeling of axonal microtubule bundles under tension

Microtubule bundles cross-linked by tau protein serve a variety of neurological functions including maintaining mechanical integrity of the axon, promoting axonal growth, and facilitating cargo transport. It has been observed that axonal damage in traumatic brain injury leads to bundle disorientation, loss of axonal viability, and cognitive impairment. This study investigates the initial mechanical response of axonal microtubule bundles under uniaxial tension using a discrete bead-spring representation. Mechanisms of failure due to traumatic stretch loading and their impact on the mechanical response and stability are also characterized. This study indicates that cross-linked axonal microtubule bundles in tension display stiffening behavior similar to a power-law relationship from nonaffine network deformations. Stretching of cross-links and microtubule bending were the primary deformation modes at low stresses. Microtubule stretch was negligible up to tensile stresses of ∼1 MPa. Bundle failure occurred by failure of cross-links leading to pull-out of microtubules and loss of bundle integrity. This may explain the elongation, undulation, and delayed elasticity of axons following traumatic stretch loading. More extensively cross-linked bundles withstood higher tensile stresses before failing. The bundle mechanical behavior uncovered by these computational techniques should guide future experiments on stretch-injured axons.