Characterization of the three-dimensional kinematic behavior of axons in central nervous system white matter

An Ogden hyperelastic 3D micromechanical model to depict Poynting effect in brain white matter

Shear and torsional load on soft solids such as brain white matter purportedly exhibits the Poynting Effect. It is a typical nonlinear phenomenon associated with soft materials whereby they tend to elongate (positive Poynting effect) or contract (negative Poynting effect) in a direction perpendicular to the shearing or twisting plane. In this research, a novel 3D micromechanical Finite Element Model (FEM) has been formulated to describe the Poynting effect in bi-phasic modeled brain white matter (BWM) representative volume element (RVE) with axons tracts embedded in surrounding extracellular matrix (ECM) for simulating brain matter's response to pure and simple shear. In the presented BWM 3D FEM, nonlinear Ogden hyper-elastic material model is deployed to interpret axons and ECM material phases. The modeled bi-phasic RVEs have axons tied to the surrounding ECM. In this proof-of-concept (POC) FEM, three simple shear loading configurations and a pure shear case were analyzed. Root mean square deviation (RMSD) was calculated for stress and deformation response plots to understand the effect of axon-ECM orientations and loading conditions on the degree of Poynting behavior. Variations in normal stresses (S11, S22, or S33) perpendicular to the shear plane underscored the significance of axonal fiber-matrix interactions. From the simulated ensemble of cases, a transitional dominance trend was noticed, as simple sheared axons showed pronounced Poynting behavior, but shear deformation build-up in the purely sheared brain model exhibited the highest Poynting behavior at higher strain % limits. At lower strain limits, simple shear imparted across and perpendicular to axonal tract directions emerged as the dominant Poynting effect configurations. At high strains, the stress-strain% plots manifested mild strain stiffening effects and bending stresses in purely sheared axons, substantiated the strong non-linearity in brain tissues' response.

Mathematical Optimisation of Magnetic Nanoparticle Diffusion in the Brain White Matter

Magnetic nanoparticles (MNPs) are a promising drug delivery system to treat brain diseases, as the particle transport trajectory can be manipulated by an external magnetic field. However, due to the complex microstructure of brain tissues, particularly the arrangement of nerve fibres in the white matter (WM), how to achieve desired drug distribution patterns, e.g., uniform distribution, is largely unknown. In this study, by adopting a mathematical model capable of capturing the diffusion trajectories of MNPs, we conducted a pilot study to investigate the effects of key parameters in the MNP delivery on the particle diffusion behaviours in the brain WM microstructures. The results show that (i) a uniform distribution of MNPs can be achieved in anisotropic tissues by adjusting the particle size and magnetic field; (ii) particle size plays a key role in determining MNPs' diffusion behaviours. The magnitude of MNP equivalent diffusivity is reversely correlated to the particle size. The MNPs with a dimension greater than 90 nm cannot reach a uniform distribution in the brain WM even in an external magnitude field; (iii) axon tortuosity may lead to transversely anisotropic MNP transport in the brain WM; however, this effect can be mitigated by applying an external magnetic field perpendicular to the local axon track. This study not only advances understanding to answer the question of how to optimise MNP delivery, but also demonstrates the potential of mathematical modelling to help achieve desired drug distributions in biological tissues with a complex microstructure.

Harmonic viscoelastic response of 3D histology-informed white matter model

White matter (WM) consists of bundles of long axons embedded in a glial matrix, which lead to anisotropic mechanical properties of brain tissue, and this complicates direct numerical simulations of WM viscoelastic response. The detailed axonal geometry contains scales that range from axonal diameter (microscale) to many diameters (mesoscale) imposing an additional challenge to numerical simulations. Here we describe the development of a 3D homogenization model for the central nervous system (CNS) that accounts for the anisotropy introduced by the axon/neuroglia composite, the axonal trace curvature, and the tissue dynamic response in the frequency domain. Homogenized models that allow the incorporation of all the above factors are important for accurately simulating the tissue's mechanical behavior, and this in turn is essential in interpreting non-invasive elastography measurements. Geometric and material parameters affect the material properties and thus the response of the brain tissue. More complex, orthotropic, or anisotropic material properties are to be considered as necessitated by the 3D tissue structure. An assembly of micro-scale 3D representative elemental volumes (REVs) is constructed, leading to an integrated mesoscale WM finite element model. Assemblies of microscopic REVs, with orientations based on geometrical reconstructions driven by confocal microscopy data are employed to form the elements of the WM model. Each REV carries local material properties based on a finite element model of biphasic (axon-glial matrix) unidirectional composite. The viscoelastic response of the microscopic REVs is extracted based on geometric information and fiber volume fractions calculated from the relative distance between the local elements and global axonal trace. The response of the WM tissue model is homogenized by averaging the shear moduli over the total volume (thus deriving effective properties) under realistic external loading conditions. Under harmonic shear loading, it is proven that that the effective transverse shear moduli are higher than the axial moduli when the axon moduli are higher than the glial. Methodologically, the process of using micro-scale 3D REVs to describe more complex axon geometries avoids the partition process in traditional composite finite element methods (based on partition of finite element grids) and constitutes a robust algorithm to automatically build a WM model based on available axonal trace information. Analytically, the model provides unmatched simulation flexibility and computational power as the position, orientation, and the magnitude of each tissue building block is calculated using real tissue data, as are the training and testing processes at each level of the multiscale WM tissue.

Morphological changes in glial cells arrangement under mechanical loading: A quantitative study

The mechanical properties and microstructure of brain tissue, as its two main physical parameters, could be affected by mechanical stimuli. In previous studies, microstructural alterations due to mechanical loading have received less attention than the mechanical properties of the tissue. Therefore, the current study aimed to investigate the effect of ex-vivo mechanical forces on the micro-architecture of brain tissue including axons and glial cells. A three-step loading protocol (i.e., loading-recovery-loading) including eight strain levels from 5% to 40% was applied to bovine brain samples with axons aligned in one preferred direction (each sample experienced only one level of strain). After either the first or secondary loading step, the samples were fixed, cut in planes parallel and perpendicular to the loading direction, and stained for histology. The histological images were analyzed to measure the end-to-end length of axons and glial cell-cell distances. The results showed that after both loading steps, as the strain increased, the changes in the cell nuclei arrangement in the direction parallel to axons were more significant compared to the other two perpendicular directions. Based on this evidence, we hypothesized that the spatial pattern of glial cells is highly affected by the orientation of axonal fibers. Moreover, the results revealed that in both loading steps, the maximum cell-cell distance occurred at 15% strain, and this distance decreased for higher strains. Since 15% strain is close to the previously reported brain injury threshold, this evidence could suggest that at higher strains, the axons start to rupture, causing a reduction in the displacement of glial cells. Accordingly, it was concluded that more attention to glial cells' architecture during mechanical loading may lead to introduce a new biomarker for brain injury.

Elucidating Axonal Injuries Through Molecular Modelling of Myelin Sheaths and Nodes of Ranvier

Around half of the traumatic brain injuries are thought to be axonal damage. Disruption of the cellular membranes, or alternatively cytoskeletal damage has been suggested as possible injury trigger. Here, we have used molecular models to have a better insight on the structural and mechanical properties of axon sub-cellular components. We modelled myelin sheath and node of Ranvier as lipid bilayers at a coarse grained level. We built ex-novo a model for the myelin. Lipid composition and lipid saturation were based on the available experimental data. The model contains 17 different types of lipids, distributed asymmetrically between two leaflets. Molecular dynamics simulations were performed to characterize the myelin and node-of-Ranvier bilayers at equilibrium and under deformation and compared to previous axolemma simulations. We found that the myelin bilayer has a slightly higher area compressibility modulus and higher rupture strain than node of Ranvier. Compared to the axolemma in unmyelinated axon, mechanoporation occurs at 50% higher strain in the myelin and at 23% lower strain in the node of Ranvier in myelinated axon. Combining the results with finite element simulations of the axon, we hypothesizes that myelin does not rupture at the thresholds proposed in the literature for axonal injury while rupture may occur at the node of Ranvier. The findings contribute to increases our knowledge of axonal sub-cellular components and help to understand better the mechanism behind axonal brain injury.

Regional Neurodegeneration in vitro: The Protective Role of Neural Activity

Traumatic brain injury is a devastating public health problem, the eighth leading cause of death across the world. To improve our understanding of how injury at the cellular scale affects neural circuit function, we developed a protocol to precisely injure individual neurons within an in vitro neural network. We used high speed calcium imaging to estimate alterations in neural activity and connectivity that occur followed targeted microtrauma. Our studies show that mechanically injured neurons inactivate following microtrauma and eventually re-integrate into the network. Single neuron re-integration is dependent on its activity prior to injury and initial connections in the network: more active and integrated neurons are more resistant to microtrauma and more likely to re-integrate into the network. Micromechanical injury leads to neuronal death 6 h post-injury in a subset of both injured and uninjured neurons. Interestingly, neural activity and network participation after injury were associated with survival in linear discriminate analysis (77.3% correct prediction, Wilks' Lambda = 0.838). Based on this observation, we modulated neuronal activity to rescue neurons after microtrauma. Inhibition of neuronal activity provided much greater survivability than did activation of neurons (ANOVA, p < 0.01 with post-hoc Tukey HSD, p < 0.01). Rescue of neurons by blocking activity in the post-acute period is partially mediated by mitochondrial energetics, as we observed silencing neurons after micromechanical injury led to a significant reduction in mitochondrial calcium accumulation. Overall, the present study provides deeper insight into the propagation of injury within networks, demonstrating that together the initial activity, network structure, and post-injury activity levels contribute to the progressive changes in a neural circuit after mechanical trauma.

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.