Hyperelastic material properties of axonal fibers in brain white matter

A theoretical framework for predicting the heterogeneous stiffness map of brain white matter tissue

Finding the stiffness map of biological tissues is of great importance in evaluating their healthy or pathological conditions. However, due to the heterogeneity and anisotropy of biological fibrous tissues, this task presents challenges and significant uncertainty when characterized only by single-mode loading experiments. In this study, we propose a new theoretical framework to map the stiffness landscape of fibrous tissues, specifically focusing on brain white matter tissue. Initially, a finite element (FE) model of the fibrous tissue was subjected to six loading cases, and their corresponding stress-strain curves were characterized. By employing multiobjective optimization, the material constants of an equivalent anisotropic material model were inversely extracted to best fit all six loading modes simultaneously. Subsequently, large-scale FE simulations were conducted, incorporating various fiber volume fractions and orientations, to train a convolutional neural network capable of predicting the equivalent anisotropic material properties solely based on the fibrous architecture of any given tissue. The proposed method, leveraging brain fiber tractography, was applied to a localized volume of white matter, demonstrating its effectiveness in precisely mapping the anisotropic behavior of fibrous tissue. In the long-term, the proposed method may find applications in traumatic brain injury, brain folding studies, and neurodegenerative diseases, where accurately capturing the material behavior of the tissue is crucial for simulations and experiments.

Toward understanding the brain tissue behavior due to preconditioning: an experimental study and RVE approach

Brain tissue under preconditioning, as a complex issue, refers to repeated loading-unloading cycles applied in mechanical testing protocols. In previous studies, only the mechanical behavior of the tissue under preconditioning was investigated; However, the link between macrostructural mechanical behavior and microstructural changes in brain tissue remains underexplored. This study aims to bridge this gap by investigating bovine brain tissue responses both before and after preconditioning. We employed a dual approach: experimental mechanical testing and computational modeling. Experimental tests were conducted to observe microstructural changes in mechanical behavior due to preconditioning, with a focus on axonal damage. Concurrently, we developed multiscale models using statistically representative volume elements (RVE) to simulate the tissue's microstructural response. These RVEs, featuring randomly distributed axonal fibers within the extracellular matrix, provide a realistic depiction of the white matter microstructure. Our findings show that preconditioning induces significant changes in the mechanical properties of brain tissue and affects axonal integrity. The RVE models successfully captured localized stresses and facilitated the microscopic analysis of axonal injury mechanisms. These results underscore the importance of considering both macro and micro scales in understanding brain tissue behavior under mechanical loading. This comprehensive approach offers valuable insights into mechanotransduction processes and improves the analysis of microstructural phenomena in brain tissue.

Computational modeling and uncertainty prediction of hyperelastic constitutive responses of damaged brain tissue under different temperature and strain rates

The effect of strain rate and temperature on the hyperelastic material stress-strain characteristics of the damaged porcine brain tissue is evaluated in this present work. The desired constitutive responses are obtained using the commercially available finite element (FE) tool ABAQUS, utilizing 8-noded brick elements. The model's accuracy has been verified by comparing the results from the previously published literature. Further, the stress-strain behavior of the brain tissue is evaluated by varying the damages at various strain rates and temperatures (13, 20, 27, and 37°C) under compression test. Additionally, the sensitivity analysis of the model is computed to check the effect of input parameters, that is, the temperature, strain rate, and damages on the material properties (shear modulus). The modeling and discussion sections enumerate the inclusive features and model capabilities.

Porosity-permeability tensor relationship of closely and randomly packed fibrous biomaterials and biological tissues: Application to the brain white matter

The constitutive model for the porosity-permeability relationship is a powerful tool to estimate and design the transport properties of porous materials, which has attracted significant attention for the advancement of novel materials. However, in comparison with other materials, biomaterials, especially natural and artificial tissues, have more complex microstructures e.g. high anisotropy, high randomness of cell/fibre dimensions/position and very low porosity. Consequently, a reliable microstructure-permeability relationship of fibrous biomaterials has proven elusive. To fill this gap, we start a mathematical derivation from the fundamental brain white matter (WM) formed by nerve fibres. This is augmented by a numerical characterisation and experimental validations to obtain an anisotropic permeability tensor of the brain WM as a function of the tissue porosity. A versatile microstructure generation software (MicroFiM) for fibrous biomaterial with complex microstructure and low porosity was built accordingly and made freely accessible here. Moreover, we propose an anisotropic poro-hyperelastic model enhanced by the newly defined porosity-permeability tensor relationship which precisely captures the tissues macro-scale permeability changes due to the microstructural deformation in an infusion scenario. The constitutive model, theories and protocols established in this study will both provide improved design strategies to tailor the transport properties of fibrous biomaterials and enable the non-invasive characterisation of the transport properties of biological tissues. This will lead to the provision of better patient-specific medical treatments, such as drug delivery. STATEMENT OF SIGNIFICANCE: Due to the microstructural complexity, a reliable microstructure-permeability relationship of fibrous biomaterials has proven elusive, which hinders our way of tuning the fluid transport property of the biomaterials by directly programming their microstructure. The same problem hinders non-invasive characterisations of fluid transport properties in biological tissues, which can significantly improve the efficiency of treatments e.g. drug delivery, directly from the tissues accessible microstructural information, e.g. porosity. Here, we developed a validated mathematical formulation to link the random microstructure to a fibrous material's macroscale permeability tensor. This will advance our capability to design complex biomaterials and make it possible to non-invasively characterise the permeability of living tissues for precise treatment planning. The newly established theory and protocol can be easily adapted to various types of fibrous biomaterials.

Histology-informed multiscale modeling of human brain white matter

In this study, we propose a novel micromechanical model for the brain white matter, which is described as a heterogeneous material with a complex network of axon fibers embedded in a soft ground matrix. We developed this model in the framework of RVE-based multiscale theories in combination with the finite element method and the embedded element technique for embedding the fibers. Microstructural features such as axon diameter, orientation and tortuosity are incorporated into the model through distributions derived from histological data. The constitutive law of both the fibers and the matrix is described by isotropic one-term Ogden functions. The hyperelastic response of the tissue is derived by homogenizing the microscopic stress fields with multiscale boundary conditions to ensure kinematic compatibility. The macroscale homogenized stress is employed in an inverse parameter identification procedure to determine the hyperelastic constants of axons and ground matrix, based on experiments on human corpus callosum. Our results demonstrate the fundamental effect of axon tortuosity on the mechanical behavior of the brain's white matter. By combining histological information with the multiscale theory, the proposed framework can substantially contribute to the understanding of mechanotransduction phenomena, shed light on the biomechanics of a healthy brain, and potentially provide insights into neurodegenerative processes.

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.

Mechanical hierarchy in the formation and modulation of cortical folding patterns

The important mechanical parameters and their hierarchy in the growth and folding of the human brain have not been thoroughly understood. In this study, we developed a multiscale mechanical model to investigate how the interplay between initial geometrical undulations, differential tangential growth in the cortical plate, and axonal connectivity form and regulate the folding patterns of the human brain in a hierarchical order. To do so, different growth scenarios with bilayer spherical models that features initial undulations on the cortex and uniform or heterogeneous distribution of axonal fibers in the white matter were developed, statistically analyzed, and validated by the imaging observations. The results showed that the differential tangential growth is the inducer of cortical folding, and in a hierarchal order, high-amplitude initial undulations on the surface and axonal fibers in the substrate regulate the folding patterns and determine the location of gyri and sulci. The locations with dense axonal fibers after folding settle in gyri rather than sulci. The statistical results also indicated that there is a strong correlation between the location of positive (outward) and negative (inward) initial undulations and the locations of gyri and sulci after folding, respectively. In addition, the locations of 3-hinge gyral folds are strongly correlated with the initial positive undulations and locations of dense axonal fibers. As another finding, it was revealed that there is a correlation between the density of axonal fibers and local gyrification index, which has been observed in imaging studies but not yet fundamentally explained. This study is the first step in understanding the linkage between abnormal gyrification (surface morphology) and disruption in connectivity that has been observed in some brain disorders such as Autism Spectrum Disorder. Moreover, the findings of the study directly contribute to the concept of the regularity and variability of folding patterns in individual human brains.

Origins of brain tissue elasticity under multiple loading modes by analyzing the microstructure-based models

Constitutive behaviors and material properties of brain tissue play an essential role in accurately modeling its mechanical responses. However, the measured mechanical behaviors of brain tissue exhibit a large variability, and the reported elastic modulus can differ by orders of magnitude. Here we develop the micromechanical models based on the actual microstructure of the longitudinally anisotropic plane of brain tissue to investigate the microstructural origins of the large variability. Specifically, axonal fiber bundles with the specified configurations are distributed in an equivalent matrix. All micromechanical models are subjected to multiple loading modes, such as tensile, compressive, and shear loading, under periodic boundary conditions. The predicted results agree well with the experimental results. Furthermore, we investigate how brain tissue elasticity varies with its microstructural features. It is revealed that the large variability in brain tissue elasticity stems from the volume fraction of axonal fiber, the aspect ratio of axonal fiber, and the distribution of axonal fiber orientation. The volume fraction has the greatest impact on the mechanical behaviors of brain tissue, followed by the distribution of axonal fiber orientation, then the aspect ratio. This study provides critical insights for understanding the microstructural origins of the large variability in brain tissue elasticity.

Effect of Particle Size and Surface Charge on Nanoparticles Diffusion in the Brain White Matter

Purpose

Brain disorders have become a serious problem for healthcare worldwide. Nanoparticle-based drugs are one of the emerging therapies and have shown great promise to treat brain diseases. Modifications on particle size and surface charge are two efficient ways to increase the transport efficiency of nanoparticles through brain-blood barrier; however, partly due to the high complexity of brain microstructure and limited visibility of Nanoparticles (NPs), our understanding of how these two modifications can affect the transport of NPs in the brain is insufficient.

Methods

In this study, a framework, which contains a stochastic geometric model of brain white matter (WM) and a mathematical particle tracing model, was developed to investigate the relationship between particle size/surface charge of the NPs and their effective diffusion coefficients (D) in WM.

Results

The predictive capabilities of this method have been validated using published experimental tests. For negatively charged NPs, both particle size and surface charge are positively correlated with D before reaching a size threshold. When Zeta potential (Zp) is less negative than -10 mV, the difference between NPs’ D in WM and pure interstitial fluid (IF) is limited.

Conclusion

A deeper understanding on the relationships between particle size/surface charge of NPs and their D in WM has been obtained. The results from this study and the developed modelling framework provide important tools for the development of nano-drugs and nano-carriers to cure brain diseases.