Damage and Failure of Axonal Microtubule under Extreme High Strain Rate: An In-Silico Molecular Dynamics Study

Can repetitive mechanical motion cause structural damage to axons?

Biological structures have evolved to very efficiently generate, transmit, and withstand mechanical forces. These biological examples have inspired mechanical engineers for centuries and led to the development of critical insights and concepts. However, progress in mechanical engineering also raises new questions about biological structures. The past decades have seen the increasing study of failure of engineered structures due to repetitive loading, and its origin in processes such as materials fatigue. Repetitive loading is also experienced by some neurons, for example in the peripheral nervous system. This perspective, after briefly introducing the engineering concept of mechanical fatigue, aims to discuss the potential effects based on our knowledge of cellular responses to mechanical stresses. A particular focus of our discussion are the effects of mechanical stress on axons and their cytoskeletal structures. Furthermore, we highlight the difficulty of imaging these structures and the promise of new microscopy techniques. The identification of repair mechanisms and paradigms underlying long-term stability is an exciting and emerging topic in biology as well as a potential source of inspiration for engineers.

Beyond symptomatic diagnosis of mild traumatic brain injury

It is commonly assumed that there is no brain injury if there are no noticeable symptoms following a head impact. There is growing evidence that traumatic brain injuries can occur with no outward symptoms and that the damage from these injuries can accumulate over time resulting in disease and impairment later in life. It is time to rethink the role that symptoms play in traumatic brain injury and adopt a quantitative understanding of brain health at the cellular level to improve the way we diagnose, prevent, and ultimately heal brain injury.

Viscoelastic damage evaluation of the axon

In this manuscript, we have studied the microstructure of the axonal cytoskeleton and adopted a bottom-up approach to evaluate the mechanical responses of axons. The cytoskeleton of the axon includes the microtubules (MT), Tau proteins (Tau), neurofilaments (NF), and microfilaments (MF). Although most of the rigidity of the axons is due to the MT, the viscoelastic response of axons comes from the Tau. Early studies have shown that NF and MF do not provide significant elasticity to the overall response of axons. Therefore, the most critical aspect of the mechanical response of axons is the microstructural topology of how MT and Tau are connected and construct the cross-linked network. Using a scanning electron microscope (SEM), the cross-sectional view of the axons revealed that the MTs are organized in a hexagonal array and cross-linked by Tau. Therefore, we have developed a hexagonal Representative Volume Element (RVE) of the axonal microstructure with MT and Tau as fibers. The matrix of the RVE is modeled by considering a combined effect of NF and MF. A parametric study is done by varying fiber geometric and mechanical properties. The Young’s modulus and spacing of MT are varied between 1.5 and 1.9 GPa and 20–38 nm, respectively. Tau is modeled as a 3-parameter General Maxwell viscoelastic material. The failure strains for MT and Tau are taken to be 50 and 40%, respectively. A total of 4 RVEs are prepared for finite element analysis, and six loading cases are inspected to quantify the three-dimensional (3D) viscoelastic relaxation response. The volume-averaged stress and strain are then used to fit the relaxation Prony series. Next, we imposed varying strain rates (between 10/sec to 50/sec) on the RVE and analyzed the axonal failure process. We have observed that the 40% failure strain of Tau is achieved in all strain rates before the MT reaches its failure strain of 50%. The corresponding axonal failure strain and stress vary between 6 and 11% and 5–19.8 MPa, respectively. This study can be used to model macroscale axonal aggregate typical of the white matter region of the brain tissue.

Shock-Induced Damage Mechanism of Perineuronal Nets

The perineuronal net (PNN) region of the brain’s extracellular matrix (ECM) surrounds the neural networks within the brain tissue. The PNN is a protective net-like structure regulating neuronal activity such as neurotransmission, charge balance, and action potential generation. Shock-induced damage of this essential component may lead to neuronal cell death and neurodegenerations. The shock generated during a vehicle accident, fall, or improvised device explosion may produce sufficient energy to damage the structure of the PNN. The goal is to investigate the mechanics of the PNN in reaction to shock loading and to understand the mechanical properties of different PNN components such as glycan, GAG, and protein. In this study, we evaluated the mechanical strength of PNN molecules and the interfacial strength between the PNN components. Afterward, we assessed the PNN molecules’ damage efficiency under various conditions such as shock speed, preexisting bubble, and boundary conditions. The secondary structure altercation of the protein molecules of the PNN was analyzed to evaluate damage intensity under varying shock speeds. At a higher shock speed, damage intensity is more elevated, and hyaluronan (glycan molecule) is most likely to break at the rigid junction. The primary structure of the protein molecules is least likely to fail. Instead, the molecules’ secondary bonds will be altered. Our study suggests that the number of hydrogen bonds during the shock wave propagation is reduced, which leads to the change in protein conformations and damage within the PNN structure. As such, we found a direct connection between shock wave intensity and PNN damage.

Effects of Focal Axonal Swelling Level on the Action Potential Signal Transmission

Focal axon swelling refers to localized swelling in axons that may occur because of trauma (e.g., traumatic brain injury) or neurodegenerative diseases (e.g., Alzheimer's disease). Since the swelling region can be many times larger than its original axon size, many researchers hypothesize that the swelling can alter the action potential (AP) signal. This article discusses the results of a series of newly developed computational studies to elucidate the possible intervention or blockage of AP signals due to swelling in the brain. We argue that the spherical geometry of the swelling site with its enlarged conducting interior causes the entering electric currents to spread evenly over the entire swelled membrane. As such, when the swelled surface becomes larger than the threshold size, the electric current will spread too thin to trigger the AP to spike. In this study, we have used a hybrid membrane model to simulate AP propagation across axons of different radii and swelling radii. We used an integrated model where a cylindrical symmetric 2D model is used to examine the electric current inside a spherical swelling site. In addition, two 1D models are used to capture the current flows along the upstream and downstream stretch before and after the swelling site. The parameters for this model are obtained from literature dedicated to modeling the experimental outcomes of mammal neurons. We observed two factors, which simultaneously affect AP transmission across a swelled axon: a) the axon radius and b) the ratio of the swelled and unswelled axon radii. In general, a thicker axon needs a smaller swelling size and axon ratio to block AP transmission. On the other hand, a thinner axon will reach the threshold at a larger swelling size and axon ratio. When only swelling size is considered, then thinner axons will block AP transmission at a smaller swelling radius. The AP transmission delay inside the swelled region determines whether the AP transmits forward or not. Notably, the blockage is worse if the AP fires at a high frequency. An increase in the charging and reset time due to swelling appears to be the main reason for the variation in axonal response.

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.

Viscoelastic Response of Neurofilaments: An Atomistic Simulation Approach

Existent literature has limitations regarding the mechanical behavior of axonal cytoskeletal components in a high strain rate scenario, which is mainly due to limitations regarding the structure of some components such as tau protein and neurofilaments (NF). This study performs molecular dynamics (MD) simulations on NFs to extract their strain rate-dependent behavior. It is found that they are highly stretchable and show multiple stages of unfolding. Furthermore, NFs show high tensile stiffness. Also, viscoelastic modeling shows that they correspond to simplified viscoelastic models. This study effectively enhances the existent axonal models focusing on axonal injury.

Mechanical behavior of actin and spectrin subjected to high strain rate: A molecular dynamics simulation study

Recent nanoscopy and super-resolution microscopy studies have substantiated the structural contribution of periodic actin-spectrin lattice to the axonal cytoskeleton of neuron. However, sufficient mechanical insight is not present for spectrin and actin-spectrin network, especially in high strain rate scenario. To quantify the mechanical behavior of actin-spectrin cytoskeleton in such conditions, this study determines individual stretching characteristics of actin and spectrin at high strain rate by molecular dynamics (MD) simulation. The actin-spectrin separation criteria are also determined. It is found that both actin and spectrin have high stiffness when susceptible to high strain rate and show strong dependence on applied strain rate. The stretching stiffness of actin and forced unfolding mechanism of spectrin are in harmony with the current literature. Actin-spectrin model provides novel insight into their interaction and separation stretch. It is shown that the region vulnerable to failure is the actin-spectrin interface at lower strain rate, while it is the inter-repeat region of spectrin at higher strain rate.

Coarse-Grained Simulation of Mechanical Properties of Single Microtubules With Micrometer Length

Microtubules are one of the most important components in the cytoskeleton and play a vital role in maintaining the shape and function of cells. Because single microtubules are some micrometers long, it is difficult to simulate such a large system using an all-atom model. In this work, we use the newly developed convolutional and K-means coarse-graining (CK-CG) method to establish an ultra-coarse-grained (UCG) model of a single microtubule, on the basis of the low electron microscopy density data of microtubules. We discuss the rationale of the micro-coarse-grained microtubule models of different resolutions and explore microtubule models up to 12-micron length. We use the devised microtubule model to quantify mechanical properties of microtubules of different lengths. Our model allows mesoscopic simulations of micrometer-level biomaterials and can be further used to study important biological processes related to microtubule function.

Domain focused and residue focused phosphorylation effect on tau protein: A molecular dynamics simulation study

Phosphorylation has been hypothesized to alter the ability of tau protein to bind with microtubules (MT), and pathological level of phosphorylation can incorporate formation of Paired Helical Filaments (PHF) in affected tau. Study of the effect of phosphorylation on different domains of tau (projection domain, microtubule binding sites and N-terminus tail) is important to obtain insight about tau neuropathology. In an earlier study, we have already obtained the mechanical properties and behavior of single tau and dimerized tau and observed tau-MT interaction for normal level of phosphorylation. This study attempts to obtain insights on the effect of phosphorylation on different domains of tau, using molecular dynamics (MD) simulation with the aid of CHARMM force field under high strain rate. It also determines the effect of residue focused phosphorylation on tau-MT interaction and tau accumulation tendency. The results show that for single tau protein, unfolding stiffness does not differ significantly due to phosphorylation, but stretching stiffness can be much higher than the normally phosphorylated protein. For dimerized tau protein, the stretching required to separate the protein forming the dimer is similar for phosphorylation in individual domains but is significantly less in case of phosphorylation in all domains. For tau-MT interaction simulations, it is found that for normal phosphorylation, the tau separation from MT occurs at higher strain for phosphorylation in projection domain and N-terminus tail, and earlier for phosphorylation in all domains altogether than the normal phosphorylation state. The residue focused phosphorylation study also shows that tau separates earlier from MT and shows stronger accumulation tendency at the phosphorylated state, while preserving the highly stretchable and flexible characteristic of tau. This study provides important insight on mechanochemical phenomena relevant to traumatic brain injury (TBI) scenario, where the result of mechanical loading and posttranslational modification as well as conformation decides the mechanical behavior.

On the Molecular Level Cavitation in Soft Gelatin Hydrogel

We have studied the molecular level cavitation mechanisms and bubble growth kinetics in soft gelatin hydrogel and water. The apparent difference in cavitation threshold pressure between that generates in pure water and that in gelatin hydrogel is considered. Gelatin, which is derived from collagen, is frequently used as a brain simulant material. In liquid, cavitation bubble is created when surrounding pressure drops below the saturation vapor pressure. In principle, a cavitation bubble should continue to grow as long as tensile pressure continues to increase in the system. In our study, using molecular dynamics simulation, we have investigated the pressure requirement for a nanoscale cavitation to grow in water and gel. First, we have modeled a gel like structure with a preexisting bubble of 5 nm radius. A control model containing a 5 nm bubble in pure water is also created. Then, we have applied hydrostatic tensile pressure at two different expansion rates in the gel and water models. The results show that a gel-like structure requires higher pressure for the cavitation to grow, and both gel and water models exhibit strain rate effect on the cavitation threshold pressure. We have also found that the cavitation collapse time is dominated by the viscosity of the medium.