Microtubule Polymerization and Cross-Link Dynamics Explain Axonal Stiffness and Damage

Sex Differences in Axonal Dynamic Responses Under Realistic Tension Using Finite Element Models

Existing axonal finite element models do not consider sex morphological differences or the fidelity in dynamic input. To facilitate a systematic investigation into the micromechanics of diffuse axonal injury, we develop a parameterized modeling approach for automatic and efficient generation of sex-specific axonal models according to specified geometrical parameters. Baseline female and male axonal models in the corpus callosum with random microtubule (MT) gap configurations are generated for model calibration and evaluation. They are then used to simulate a realistic tensile loading consisting of both a loading and a recovery phase (to return to an initial undeformed state) generated from dynamic corpus callosum fiber strain in a real-world head impact simulation. We find that MT gaps and the dynamic recovery phase are both critical to successfully reproduce MT undulation as observed experimentally, which has not been reported before. This strengthens confidence in model dynamic responses. A statistical approach is further employed to aggregate axonal responses from a large sample of random MT gap configurations for both female and male axonal models (n = 10,000 each). We find that peak strains in MTs and the Ranvier node and associated neurofilament failures in female axons are substantially higher than those in male axons because there are fewer MTs in the former and also because of the random nature of MT gap locations. Despite limitations in various model assumptions as a result of limited experimental data currently available, these findings highlight the need to systematically characterize MT gap configurations and to ensure a realistic model input for axonal dynamic simulations. Finally, this study may offer fresh and improved insight into the biomechanical basis of sex differences in brain injury, and sets the stage for more systematic investigations at the microscale in the future, both numerically and experimentally.

Dynamics of the axon plasma membrane skeleton

It was recently revealed via super-resolution microscopy experiments that the axon plasma membrane skeleton (APMS) comprises a series of periodically arranged azimuthal actin rings connected via longitudinal spectrin filaments forming an orthotropic network. The common perception is that APMS enhances structural stability of the axon but its impact on axon deformation is unknown. To investigate the response of the APMS to extension, we introduce a coarse-grain molecular dynamics model consisting of actin particles forming rings and chains of particles representing spectrin tetramers with repeats than can unfold. We observe that the shape of force-extension curve is initially linear and the force level depends on the extension rate. Even during the initial deformation stage, unfolding of spectrin repeats occurs, but the saw-tooth shape of the corresponding force-extension curve observed in the case of one spectrin tetramer does not appear in the case of the entire APMS. The reason is that spectrin unfolding is not synchronized across filaments during extension. If actin-spectrin associations remain intact, the force-extension response reaches a perfectly plastic region because of increased spectrin unfolding frequency. However, when actin-spectrin links dissociate, which can happen at moderate and high extension rates, APMS softens and the resistance force decreases linearly as the axon elongates until it reaches a point where the APMS is completely severed. Furthermore, when the ring-to-ring distance is maintained fixed under stretch, the resistance force relaxes exponentially as a function of time due to additional unfolding of spectrin tetramers following the Kelvin-Voigt representation of the Zener model.

Unravelling the mechanotransduction pathways in Alzheimer's disease

Alzheimer's disease (AD) represents one of the most common and debilitating neurodegenerative disorders. By the end of 2040, AD patients might reach 11.2 million in the USA, around 70% higher than 2022, with severe consequences on the society. As now, we still need research to find effective methods to treat AD. Most studies focused on the tau and amyloid hypothesis, but many other factors are likely involved in the pathophysiology of AD. In this review, we summarize scientific evidence dealing with the mechanotransduction players in AD to highlight the most relevant mechano-responsive elements that play a role in AD pathophysiology. We focused on the AD-related role of extracellular matrix (ECM), nuclear lamina, nuclear transport and synaptic activity. The literature supports that ECM alteration causes the lamin A increment in the AD patients, leading to the formation of nuclear blebs and invaginations. Nuclear blebs have consequences on the nuclear pore complexes, impairing nucleo-cytoplasmic transport. This may result in tau hyperphosphorylation and its consequent self-aggregation in tangles, which impairs the neurotransmitters transport. It all exacerbates in synaptic transmission impairment, leading to the characteristic AD patient's memory loss. Here we related for the first time all the evidence associating the mechanotransduction pathway with neurons. In addition, we highlighted the entire pathway influencing neurodegenerative diseases, paving the way for new research perspectives in the context of AD and related pathologies.

Image-based axon model highlights heterogeneity in initiation of damage

Head injury simulations predict the occurrence of traumatic brain injury by placing a threshold on the calculated strains for axon tracts within the brain. However, a current roadblock to accurate injury prediction is the selection of an appropriate axon damage threshold. While several computational studies have used models of the axon cytoskeleton to investigate damage initiation, these models all employ an idealized, homogeneous axonal geometry. This homogeneous geometry with regularly spaced microtubules, evenly distributed throughout the model, overestimates axon strength because, in reality, the axon cytoskeleton is heterogeneous. In the heterogeneous cytoskeleton, the weakest cross section determines the initiation of failure, but these weak spots are not present in a homogeneous model. Addressing one source of heterogeneity in the axon cytoskeleton, we present a new semiautomated image analysis pipeline for using serial-section transmission electron micrographs to reconstruct the microtubule geometry of an axon. The image analysis procedure locates microtubules within the images, traces them throughout the image stack, and reconstructs the microtubule structure as a finite element mesh. We demonstrate the image analysis approach using a C. elegans touch receptor neuron due to the availability of high-quality serial-section transmission electron micrograph data sets. The results of the analysis highlight the heterogeneity of the microtubule structure in the spatial variation of both microtubule number and length. Simulations comparing this image-based geometry with homogeneous geometries show that structural heterogeneity in the image-based model creates significant spatial variation in deformation. The homogeneous geometries, on the other hand, deform more uniformly. Since no single homogeneous model can replicate the mechanical behavior of the image-based model, our results argue that heterogeneity in axon microtubule geometry should be considered in determining accurate axon failure thresholds.

Mathematical models of neuronal growth

The establishment of a functioning neuronal network is a crucial step in neural development. During this process, neurons extend neurites—axons and dendrites—to meet other neurons and interconnect. Therefore, these neurites need to migrate, grow, branch and find the correct path to their target by processing sensory cues from their environment. These processes rely on many coupled biophysical effects including elasticity, viscosity, growth, active forces, chemical signaling, adhesion and cellular transport. Mathematical models offer a direct way to test hypotheses and understand the underlying mechanisms responsible for neuron development. Here, we critically review the main models of neurite growth and morphogenesis from a mathematical viewpoint. We present different models for growth, guidance and morphogenesis, with a particular emphasis on mechanics and mechanisms, and on simple mathematical models that can be partially treated analytically.

Toward a Comprehensive Delineation of White Matter Tract-Related Deformation

Finite element (FE) models of the human head are valuable instruments to explore the mechanobiological pathway from external loading, localized brain response, and resultant injury risks. The injury predictability of these models depends on the use of effective criteria as injury predictors. The FE-derived normal deformation along white matter (WM) fiber tracts (i.e., tract-oriented strain) recently has been suggested as an appropriate predictor for axonal injury. However, the tract-oriented strain only represents a partial depiction of the WM fiber tract deformation. A comprehensive delineation of tract-related deformation may improve the injury predictability of the FE head model by delivering new tract-related criteria as injury predictors. Thus, the present study performed a theoretical strain analysis to comprehensively characterize the WM fiber tract deformation by relating the strain tensor of the WM element to its embedded fiber tract. Three new tract-related strains with exact analytical solutions were proposed, measuring the normal deformation perpendicular to the fiber tracts (i.e., tract-perpendicular strain), and shear deformation along and perpendicular to the fiber tracts (i.e., axial-shear strain and lateral-shear strain, respectively). The injury predictability of these three newly proposed strain peaks along with the previously used tract-oriented strain peak and maximum principal strain (MPS) were evaluated by simulating 151 impacts with known outcome (concussion or non-concussion). The results preliminarily showed that four tract-related strain peaks exhibited superior performance than MPS in discriminating concussion and non-concussion cases. This study presents a comprehensive quantification of WM tract-related deformation and advocates the use of orientation-dependent strains as criteria for injury prediction, which may ultimately contribute to an advanced mechanobiological understanding and enhanced computational predictability of brain injury.

Effect of Strain Rate on Single Tau, Dimerized Tau and Tau-Microtubule Interface: A Molecular Dynamics Simulation Study

Microtubule-associated protein (MAP) tau is a cross-linking molecule that provides structural stability to axonal microtubules (MT). It is considered a potential biomarker for Alzheimer's disease (AD), dementia, and other neurological disorders. It is also a signature protein for Traumatic Brain Injury (TBI) assessment. In the case of TBI, extreme dynamic mechanical energies can be felt by the axonal cytoskeletal members. As such, fundamental understandings of the responses of single tau protein, polymerized tau protein, and tau-microtubule interfaces under high-rate mechanical forces are important. This study attempts to determine the high-strain rate mechanical behavior of single tau, dimerized tau, and tau-MT interface using molecular dynamics (MD) simulation. The results show that a single tau protein is a highly stretchable soft polymer. During deformation, first, it significantly unfolds against van der Waals and electrostatic bonds. Then it stretches against strong covalent bonds. We found that tau acts as a viscoelastic material, and its stiffness increases with the strain rate. The unfolding stiffness can be ~50-500 MPa, while pure stretching stiffness can be >2 GPa. The dimerized tau model exhibits similar behavior under similar strain rates, and tau sliding from another tau is not observed until it is stretched to >7 times of original length, depending on the strain rate. The tau-MT interface simulations show that very high strain and strain rates are required to separate tau from MT suggesting Tau-MT bonding is stronger than MT subunit bonding between themselves. The dimerized tau-MT interface simulations suggest that tau-tau bonding is stronger than tau-MT bonding. In summary, this study focuses on the structural response of individual cytoskeletal components, namely microtubule (MT) and tau protein. Furthermore, we consider not only the individual response of a component, but also their interaction with each other (such as tau with tau or tau with MT). This study will eventually pave the way to build a bottom-up multiscale brain model and analyze TBI more comprehensively.

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.

Modeling links softening of myelin and spectrin scaffolds of axons after a concussion to increased vulnerability to repeated injuries

Damage to the microtubule lattice, which serves as a rigid cytoskeletal backbone for the axon, is a hallmark mechanical initiator of pathophysiology after concussion. Understanding the mechanical stress transfer from the brain tissue to the axonal cytoskeleton is essential to determine the microtubule lattice's vulnerability to mechanical injury. Here, we develop an ultrastructural model of the axon's cytoskeletal architecture to identify the components involved in the dynamic load transfer during injury. Corroborative in vivo studies were performed using a gyrencephalic swine model of concussion via single and repetitive head rotational acceleration. Computational analysis of the load transfer mechanism demonstrates that the myelin sheath and the actin/spectrin cortex play a significant role in effectively shielding the microtubules from tissue stress. We derive failure maps in the space spanned by tissue stress and stress rate to identify physiological conditions in which the microtubule lattice can rupture. We establish that a softer axonal cortex leads to a higher susceptibility of the microtubules to failure. Immunohistochemical examination of tissue from the swine model of single and repetitive concussion confirms the presence of postinjury spectrin degradation, with more extensive pathology observed following repetitive injury. Because the degradation of myelin and spectrin occurs over weeks following the first injury, we show that softening of the myelin layer and axonal cortex exposes the microtubules to higher stress during repeated incidences of traumatic brain injuries. Our predictions explain how mechanical injury predisposes axons to exacerbated responses to repeated injuries, as observed in vitro and in vivo.

Biomimicking Fiber Platform with Tunable Stiffness to Study Mechanotransduction Reveals Stiffness Enhances Oligodendrocyte Differentiation but Impedes Myelination through YAP-Dependent Regulation

A key hallmark of many diseases, especially those in the central nervous system (CNS), is the change in tissue stiffness due to inflammation and scarring. However, how such changes in microenvironment affect the regenerative process remains poorly understood. Here, a biomimicking fiber platform that provides independent variation of fiber structural and intrinsic stiffness is reported. To demonstrate the functionality of these constructs as a mechanotransduction study platform, these substrates are utilized as artificial axons and the effects of axon structural versus intrinsic stiffness on CNS myelination are independently analyzed. While studies have shown that substrate stiffness affects oligodendrocyte differentiation, the effects of mechanical stiffness on the final functional state of oligodendrocyte (i.e., myelination) has not been shown prior to this. Here, it is demonstrated that a stiff mechanical microenvironment impedes oligodendrocyte myelination, independently and distinctively from oligodendrocyte differentiation. Yes-associated protein is identified to be involved in influencing oligodendrocyte myelination through mechanotransduction. The opposing effects on oligodendrocyte differentiation and myelination provide important implications for current work screening for promyelinating drugs, since these efforts have focused mainly on promoting oligodendrocyte differentiation. Thus, the platform may have considerable utility as part of a drug discovery program in identifying molecules that promote both differentiation and myelination.

Beading of injured axons driven by tension- and adhesion-regulated membrane shape instability

The formation of multiple beads along an injured axon will lead to blockage of axonal transport and eventually neuron death, and this has been widely recognized as a hallmark of nervous system degeneration. Nevertheless, the underlying mechanisms remain poorly understood. Here, we report a combined experimental and theoretical study to reveal key factors governing axon beading. Specifically, by transecting well-developed axons with a sharp atomic force microscope probe, significant beading of the axons was triggered. We showed that adhesion was not required for beading to occur, although when present strong axon–substrate attachments seemed to set the locations for bead formation. In addition, the beading wavelength, representing the average distance between beads, was found to correlate with the size and cytoskeleton integrity of axon, with a thinner axon or a disrupted actin cytoskeleton both leading to a shorter beading wavelength. A model was also developed to explain these observations which suggest that axon beading originates from the shape instability of the membrane and is driven by the release of work done by axonal tension as well as the reduction of membrane surface energy. The beading wavelength predicted from this theory was in good agreement with our experiments under various conditions. By elucidating the essential physics behind axon beading, the current study could enhance our understanding of how axonal injury and neurodegeneration progress as well as provide insights for the development of possible treatment strategies.

Localized Axolemma Deformations Suggest Mechanoporation as Axonal Injury Trigger

Traumatic brain injuries are a leading cause of morbidity and mortality worldwide. With almost 50% of traumatic brain injuries being related to axonal damage, understanding the nature of cellular level impairment is crucial. Experimental observations have so far led to the formulation of conflicting theories regarding the cellular primary injury mechanism. Disruption of the axolemma, or alternatively cytoskeletal damage has been suggested mainly as injury trigger. However, mechanoporation thresholds of generic membranes seem not to overlap with the axonal injury deformation range and microtubules appear too stiff and too weakly connected to undergo mechanical breaking. Here, we aim to shed a light on the mechanism of primary axonal injury, bridging finite element and molecular dynamics simulations. Despite the necessary level of approximation, our models can accurately describe the mechanical behavior of the unmyelinated axon and its membrane. More importantly, they give access to quantities that would be inaccessible with an experimental approach. We show that in a typical injury scenario, the axonal cortex sustains deformations large enough to entail pore formation in the adjoining lipid bilayer. The observed axonal deformation of 10–12% agree well with the thresholds proposed in the literature for axonal injury and, above all, allow us to provide quantitative evidences that do not exclude pore formation in the membrane as a result of trauma. Our findings bring to an increased knowledge of axonal injury mechanism that will have positive implications for the prevention and treatment of brain injuries.

The model of local axon homeostasis - explaining the role and regulation of microtubule bundles in axon maintenance and pathology

Axons are the slender, cable-like, up to meter-long projections of neurons that electrically wire our brains and bodies. In spite of their challenging morphology, they usually need to be maintained for an organism's lifetime. This makes them key lesion sites in pathological processes of ageing, injury and neurodegeneration. The morphology and physiology of axons crucially depends on the parallel bundles of microtubules (MTs), running all along to serve as their structural backbones and highways for life-sustaining cargo transport and organelle dynamics. Understanding how these bundles are formed and then maintained will provide important explanations for axon biology and pathology. Currently, much is known about MTs and the proteins that bind and regulate them, but very little about how these factors functionally integrate to regulate axon biology. As an attempt to bridge between molecular mechanisms and their cellular relevance, we explain here the model of local axon homeostasis, based on our own experiments in Drosophila and published data primarily from vertebrates/mammals as well as C. elegans. The model proposes that (1) the physical forces imposed by motor protein-driven transport and dynamics in the confined axonal space, are a life-sustaining necessity, but pose a strong bias for MT bundles to become disorganised. (2) To counterbalance this risk, MT-binding and -regulating proteins of different classes work together to maintain and protect MT bundles as necessary transport highways. Loss of balance between these two fundamental processes can explain the development of axonopathies, in particular those linking to MT-regulating proteins, motors and transport defects. With this perspective in mind, we hope that more researchers incorporate MTs into their work, thus enhancing our chances of deciphering the complex regulatory networks that underpin axon biology and pathology.

The Shrinking Brain: Cerebral Atrophy Following Traumatic Brain Injury

Cerebral atrophy in response to traumatic brain injury is a well-documented phenomenon in both primary investigations and review articles. Recent atrophy studies focus on exploring the region-specific patterns of cerebral atrophy; yet, there is no study that analyzes and synthesizes the emerging atrophy patterns in a single comprehensive review. Here we attempt to fill this gap in our current knowledge by integrating the current literature into a cohesive theory of preferential brain tissue loss and by identifying common risk factors for accelerated atrophy progression. Our review reveals that observations for mild traumatic brain injury remain inconclusive, whereas observations for moderate-to-severe traumatic brain injury converge towards robust patterns: brain tissue loss is on the order of 5% per year, and occurs in the form of generalized atrophy, across the entire brain, or focal atrophy, in specific brain regions. The most common regions of focal atrophy are the thalamus, hippocampus, and cerebellum in gray matter and the corpus callosum, corona radiata, and brainstem in white matter. We illustrate the differences of generalized and focal gray and white matter atrophy on emerging deformation and stress profiles across the whole brain using computational simulation. The characteristic features of our atrophy simulations-a widening of the cortical sulci, a gradual enlargement of the ventricles, and a pronounced cortical thinning-agree well with clinical observations. Understanding region-specific atrophy patterns in response to traumatic brain injury has significant implications in modeling, simulating, and predicting injury outcomes. Computational modeling of brain atrophy could open new strategies for physicians to make informed decisions for whom, how, and when to administer pharmaceutical treatment to manage the chronic loss of brain structure and function.

Tension- and Adhesion-Regulated Retraction of Injured Axons

Damage-induced retraction of axons during traumatic brain injury is believed to play a key role in the disintegration of the neural network and to eventually lead to severe symptoms such as permanent memory loss and emotional disturbances. However, fundamental questions such as how axon retraction progresses and what physical factors govern this process still remain unclear. Here, we report a combined experimental and modeling study to address these questions. Specifically, a sharp atomic force microscope probe was used to transect axons and trigger their retraction in a precisely controlled manner. Interestingly, we showed that the retracting motion of a well-developed axon can be arrested by strong cell-substrate attachment. However, axon retraction was found to be retriggered if a second transection was conducted, albeit with a lower shrinking amplitude. Furthermore, disruption of the actin cytoskeleton or cell-substrate adhesion significantly altered the retracting dynamics of injured axons. Finally, a mathematical model was developed to explain the observed injury response of neural cells in which the retracting motion was assumed to be driven by the pre-tension in the axon and progress against neuron-substrate adhesion as well as the viscous resistance of the cell. Using realistic parameters, model predictions were found to be in good agreement with our observations under a variety of experimental conditions. By revealing the essential physics behind traumatic axon retraction, findings here could provide insights on the development of treatment strategies for axonal injury as well as its possible interplay with other neurodegenerative diseases.

An Integrated Cytoskeletal Model of Neurite Outgrowth

Neurite outgrowth underlies the wiring of the nervous system during development and regeneration. Despite a significant body of research, the underlying cytoskeletal mechanics of growth and guidance are not fully understood, and the relative contributions of individual cytoskeletal processes to neurite growth are controversial. Here, we review the structural organization and biophysical properties of neurons to make a semi-quantitative comparison of the relative contributions of different processes to neurite growth. From this, we develop the idea that neurons are active fluids, which generate strong contractile forces in the growth cone and weaker contractile forces along the axon. As a result of subcellular gradients in forces and material properties, actin flows rapidly rearward in the growth cone periphery, and microtubules flow forward in bulk along the axon. With this framework, an integrated model of neurite outgrowth is proposed that hopefully will guide new approaches to stimulate neuronal growth.

Modeling the Axon as an Active Partner with the Growth Cone in Axonal Elongation

Forces generated by the growth cone are vital for the proper development of the axon and thus brain function. Although recent experiments show that forces are generated along the axon, it is unknown whether the axon plays a direct role in controlling growth cone advance. Here, we use analytic and finite element modeling of microtubule dynamics and the activity of the molecular motors myosin and dynein to investigate mechanical force balance along the length of the axon and its effects on axonal outgrowth. Our modeling indicates that the paradoxical effects of stabilizing microtubules and the consequences of microtubule disassembly on axonal outgrowth can be explained by changes in the passive and active mechanical properties of axons. Our findings suggest that a full understanding of growth cone motility requires a consideration of the mechanical contributions of the axon. Our study not only has potential applications during neurodevelopment but might also help identify strategies to manipulate and promote axonal regrowth to treat neurodegeneration.

Utilizing a Structural Mechanics Approach to Assess the Primary Effects of Injury Loads Onto the Axon and Its Components

Diffuse axonal injury (DAI) occurs as a result of the transmission of rapid dynamic loads from the head to the brain and specifically to its neurons. Despite being one of the most common and most deleterious types of traumatic brain injury (TBI), the inherent cell injury mechanism is yet to be understood. Experimental observations have led to the formulation of different hypotheses, such as mechanoporation of the axolemma and microtubules (MTs) breakage. With the goal of singling out the mechanical aspect of DAI and to resolve the ambiguity behind its injury mechanism, a composite finite element (FE) model of a representative volume of an axon was developed. Once calibrated and validated against published experimental data, the axonal model was used to simulate injury scenarios. The resulting strain distributions along its components were then studied in dependence of strain rate and of typical modeling choices such as the applied MT constraints and tau proteins failure. Results show that oversimplifying the MT bundle kinematic entails a systematic attenuation (cf = 2.33) of the computed maximum MT strain. Nevertheless, altogether, results support the hypothesis of axolemma mechanoporation as a cell-injury trigger. Moreover, for the first time the interconnection between the axolemma and the MT bundle is shown to play a role in damage localization. The proposed FE axonal model is a valuable tool to understand the axonal injury mechanism from a mechanical perspective and could be used in turn for the definition of well-informed injury criteria at the tissue and organ level.

Physical Biology of Axonal Damage

Excessive physical impacts to the head have direct implications on the structural integrity at the axonal level. Increasing evidence suggests that tau, an intrinsically disordered protein that stabilizes axonal microtubules, plays a critical role in the physical biology of axonal injury. However, the precise mechanisms of axonal damage remain incompletely understood. Here we propose a biophysical model of the axon to correlate the dynamic behavior of individual tau proteins under external physical forces to the evolution of axonal damage. To propagate damage across the scales, we adopt a consistent three-step strategy: First, we characterize the axonal response to external stretches and stretch rates for varying tau crosslink bond strengths using a discrete axonal damage model. Then, for each combination of stretch rates and bond strengths, we average the axonal force-stretch response of n = 10 discrete simulations, from which we derive and calibrate a homogenized constitutive model. Finally, we embed this homogenized model into a continuum axonal damage model of [1-d]-type in which d is a scalar damage parameter that is driven by the axonal stretch and stretch rate. We demonstrate that axonal damage emerges naturally from the interplay of physical forces and biological crosslinking. Our study reveals an emergent feature of the crosslink dynamics: With increasing loading rate, the axonal failure stretch increases, but axonal damage evolves earlier in time. For a wide range of physical stretch rates, from 0.1 to 10 /s, and biological bond strengths, from 1 to 100 pN, our model predicts a relatively narrow window of critical damage stretch thresholds, from 1.01 to 1.30, which agrees well with experimental observations. Our biophysical damage model can help explain the development and progression of axonal damage across the scales and will provide useful guidelines to identify critical damage level thresholds in response to excessive physical forces.