Hyperelastic material properties of axonal fibers in brain white matter

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.