Explicit Modeling of White Matter Axonal Fiber Tracts in a Finite Element Brain Model

MRI-based quantification of cardiac-driven brain biomechanics for early detection of neurological disorders

We present a pipeline to quantify biomechanical environment of the brain using solely MRI-derived data in order to elucidate the role of biomechanical factors in neurodegenerative disorders. Neurological disorders, like Alzheimer's and Parkinson's diseases, are associated with physical changes, including the accumulation of amyloid-β and tau proteins, damage to the cerebral vasculature, hypertension, atrophy of the cortical gray matter, and lesions of the periventricular white matter. Alterations in the external mechanical environment of cells can trigger pathological processes, and it is known that AD causes reduced stiffness in the brain tissue during degeneration. However, there appears to be a significant lag time between microscale changes and macroscale obstruction of neurological function in the brain. Here, we present a pipeline to quantify the whole brain biomechanical environment to bridge the gap in understanding how underlying brain changes affect macroscale brain biomechanics. This pipeline enables image-based quantification of subject-specific displacement field of the whole brain to subject-specific strain, strain rate, and stress across 133 labeled functional brain regions. We have focused our development efforts on utilizing solely MRI-derived data to facilitate clinical applicability of our approach and have emphasized automation in all aspects of our methods to reduce operator dependance. Our pipeline has the potential to improve early detection of neurological disorders and facilitate the identification of disease before widespread, irreversible damage has occurred.

Effect of Human Head Shape on the Risk of Traumatic Brain Injury: A Gaussian Process Regression-based Machine Learning Approach

INTRODUCTION

Computational head injury models are promising tools for understanding and predicting traumatic brain injuries. However, most available head injury models are "average" models that employ a single set of head geometry (e.g., 50th-percentile U.S. male) without considering variability in these parameters across the human population. A significant variability of head shapes exists in U.S. Army soldiers, evident from the Anthropometric Survey of U.S. Army Personnel (ANSUR II). The objective of this study is to elucidate the effects of head shape on the predicted risk of traumatic brain injury from computational head injury models.

MATERIALS AND METHODS

Magnetic resonance imaging scans of 25 human subjects are collected. These images are registered to the standard MNI152 brain atlas, and the resulting transformation matrix components (called head shape parameters) are used to quantify head shapes of the subjects. A generative machine learning model is used to generate 25 additional head shape parameter datasets to augment our database. Head injury models are developed for these head shapes, and a rapid injurious head rotation event is simulated to obtain several brain injury predictor variables (BIPVs): Peak cumulative maximum principal strain (CMPS), average CMPS, and the volume fraction of brain exceeding an injurious CMPS threshold. A Gaussian process regression model is trained between head shape parameters and BIPVs, which is then used to study the relative sensitivity of the various BIPVs on individual head shape parameters. We distinguish head shape parameters into 2 types: Scaling components ${T_{xx}}$, ${T_{yy}}$, and ${T_{zz}}$ that capture the breadth, length, and height of the head, respectively, and shearing components (${T_{xy}},{T_{xz}},{T_{yx}},{T_{yz}},{T_{zx}}$, and ${T_{zy}}$) that capture the relative skewness of the head shape.

RESULTS

An overall positive correlation is evident between scaling components and BIPVs. Notably, a very high, positive correlation is seen between the BIPVs and the head volume. As an example, a 57% increase in peak CMPS was noted between the smallest and the largest investigated head volume parameters. The variation in shearing components ${T_{xy}},{T_{xz}},{T_{yx}},{T_{yz}},{T_{zx}}$, and ${T_{zy}}$ on average does not cause notable changes in the BIPVs. From the Gaussian process regression model, all 3 BIPVs showed an increasing trend with each of the 3 scaling components, but the BIPVs are found to be most sensitive to the height dimension of the head. From the Sobol sensitivity analysis, the ${T_{zz}}$ scaling parameter contributes nearly 60% to the total variance in peak and average CMPS; ${T_{yy}}$ contributes approximately 20%, whereas ${T_{xx}}$ contributes less than 5%. The remaining contribution is from the 6 shearing components. Unlike peak and average CMPS, the VF-CMPS BIPV is associated with relatively evenly distributed Sobol indices across the 3 scaling parameters. Furthermore, the contribution of shearing components on the total variance in this case is negligible.

CONCLUSIONS

Head shape has a considerable influence on the injury predictions of computational head injury models. Available "average" head injury models based on a 50th-percentile U.S. male are likely associated with considerable uncertainty. In general, larger head sizes correspond to greater BIPV magnitudes, which point to potentially a greater injury risk under rapid neck rotation for people with larger heads.

Validation of a computational biomechanical mouse brain model for rotational head acceleration

Recent mouse brain injury experiments examine diffuse axonal injury resulting from accelerative head rotations. Evaluating brain deformation during these events would provide valuable information on tissue level thresholds for brain injury, but there are many challenges to imaging the brain's mechanical response during dynamic loading events, such as a blunt head impact. To address this shortcoming, we present an experimentally validated computational biomechanics model of the mouse brain that predicts tissue deformation, given the motion of the mouse head during laboratory experiments. First, we developed a finite element model of the mouse brain that computes tissue strains, given the same head rotations as previously conducted in situ hemicephalic mouse brain experiments. Second, we calibrated the model using a single brain segment, and then validated the model based on the spatial and temporal strain responses of other regions. The result is a computational tool that will provide researchers with the ability to predict brain tissue strains that occur during mouse laboratory experiments, and to link the experiments to the resulting neuropathology, such as diffuse axonal injury.

Dynamic strain fields of the mouse brain during rotation

Mouse models are used to better understand brain injury mechanisms in humans, yet there is a limited understanding of biomechanical relevance, beginning with how the murine brain deforms when the head undergoes rapid rotation from blunt impact. This problem makes it difficult to translate some aspects of diffuse axonal injury from mouse to human. To address this gap, we present the two-dimensional strain field of the mouse brain undergoing dynamic rotation in the sagittal plane. Using a high-speed camera with digital image correlation measurements of the exposed mid-sagittal brain surface, we found that pure rotations (no direct impact to the skull) of 100-200 rad/s are capable of producing complex strain fields that evolve over time with respect to rotational acceleration and deceleration. At the highest rotational velocity tested, the largest tensile strains (≥ 21% elongation) in selected regions of the mouse brain approach strain thresholds previously associated with axonal injury in prior work. These findings provide a benchmark to validate the mechanical response in biomechanical computational models predicting diffuse axonal injury, but much work remains in correlating tissue deformation patterns from computational models with underlying neuropathology.

Multi-scale measurement of stiffness in the developing ferret brain

Cortical folding is an important process during brain development, and aberrant folding is linked to disorders such as autism and schizophrenia. Changes in cell numbers, size, and morphology have been proposed to exert forces that control the folding process, but these changes may also influence the mechanical properties of developing brain tissue. Currently, the changes in tissue stiffness during brain folding are unknown. Here, we report stiffness in the developing ferret brain across multiple length scales, emphasizing changes in folding cortical tissue. Using rheometry to measure the bulk properties of brain tissue, we found that overall brain stiffness increases with age over the period of cortical folding. Using atomic force microscopy to target the cortical plate, we found that the occipital cortex increases in stiffness as well as stiffness heterogeneity over the course of development and folding. These findings can help to elucidate the mechanics of the cortical folding process by clarifying the concurrent evolution of tissue properties.

Inverse identification of region-specific hyperelastic material parameters for human brain tissue

The identification of material parameters accurately describing the region-dependent mechanical behavior of human brain tissue is crucial for computational models used to assist, e.g., the development of safety equipment like helmets or the planning and execution of brain surgery. While the division of the human brain into different anatomical regions is well established, knowledge about regions with distinct mechanical properties remains limited. Here, we establish an inverse parameter identification scheme using a hyperelastic Ogden model and experimental data from multi-modal testing of tissue from 19 anatomical human brain regions to identify mechanically distinct regions and provide the corresponding material parameters. We assign the 19 anatomical regions to nine governing regions based on similar parameters and microstructures. Statistical analyses confirm differences between the regions and indicate that at least the corpus callosum and the corona radiata should be assigned different material parameters in computational models of the human brain. We provide a total of four parameter sets based on the two initial Poisson's ratios of 0.45 and 0.49 as well as the pre- and unconditioned experimental responses, respectively. Our results highlight the close interrelation between the Poisson's ratio and the remaining model parameters. The identified parameters will contribute to more precise computational models enabling spatially resolved predictions of the stress and strain states in human brains under complex mechanical loading conditions.

Rate- and Region-Dependent Mechanical Properties of Göttingen Minipig Brain Tissue in Simple Shear and Unconfined Compression

Traumatic brain injury (TBI), particularly from explosive blasts, is a major cause of casualties in modern military conflicts. Computational models are an important tool in understanding the underlying biomechanics of TBI but are highly dependent on the mechanical properties of soft tissue to produce accurate results. Reported material properties of brain tissue can vary by several orders of magnitude between studies, and no published set of material parameters exists for porcine brain tissue at strain rates relevant to blast. In this work, brain tissue from the brainstem, cerebellum, and cerebrum of freshly euthanized adolescent male Göttingen minipigs was tested in simple shear and unconfined compression at strain rates ranging from quasi-static (QS) to 300 s-1. Brain tissue showed significant strain rate stiffening in both shear and compression. Minimal differences were seen between different regions of the brain. Both hyperelastic and hyper-viscoelastic constitutive models were fit to experimental stress, considering data from either a single loading mode (unidirectional) or two loading modes together (bidirectional). The unidirectional hyper-viscoelastic models with an Ogden hyperelastic representation and a one-term Prony series best captured the response of brain tissue in all regions and rates. The bidirectional models were generally able to capture the response of the tissue in high-rate shear and all compression modes, but not the QS shear. Our constitutive models describe the first set of material parameters for porcine brain tissue relevant to loading modes and rates seen in blast injury.

A sensor-enabled cloud-based computing platform for computational brain biomechanics

BACKGROUND AND OBJECTIVES

Driven by the risk of repetitive head trauma, sensors have been integrated into mouthguards to measure head impacts in contact sports and military activities. These wearable devices, referred to as "instrumented" or "smart" mouthguards are being actively developed by various research groups and organizations. These instrumented mouthguards provide an opportunity to further study and understand the brain biomechanics due to impact. In this study, we present a brain modeling service that can use information from these sensors to predict brain injury metrics in an automated fashion.

METHODS

We have built a brain modeling platform using several of Amazon's Web Services (AWS) to enable cloud computing and scalability. We use a custom-built cloud-based finite element modeling code to compute the physics-based nonlinear response of the intracranial brain tissue and provide a frontend web application and an application programming interface for groups working on head impact sensor technology to include simulated injury predictions into their research pipeline.

RESULTS

The platform results have been validated against experimental data available in literature for brain-skull relative displacements, brain strains and intracranial pressure. The parallel processing capability of the platform has also been tested and verified. We also studied the accuracy of the custom head surfaces generated by Avatar 3D.

CONCLUSION

We present a validated cloud-based computational brain modeling platform that uses sensor data as input for numerical brain models and outputs a quantitative description of brain tissue strains and injury metrics. The platform is expected to generate transparent, reproducible, and traceable brain computing results.

Development, validation and a case study: The female finite element head model (FeFEHM)

BACKGROUND AND OBJECTIVE

Traumatic brain injuries are one of the leading causes of death and disability in the world. To better understand the interactions and forces applied in different constituents of the human head, several finite element head models have been developed throughout the years, for offering a good cost-effective and ethical approach compared to experimental tests. Once validated, the female finite element head model (FeFEHM) will allow a better understanding of injury mechanisms resulting in neuronal damage, which can later evolve into neurodegenerative diseases.

METHODS

This work encompasses the approached methodology starting from medical images and finite element modelling until the validation process using novel experimental data of brain displacements conducted on human cadavers. The material modelling of the brain is performed using an age-specific characterization of the brain using microindentation at dynamic rates and under large deformation, with a similar age to the patient used to model the FeFEHM.

RESULTS

The numerical displacement curves are in good accordance with the experimental data, displaying similar peak times and values, in all three anatomical planes. The case study result shows a similarity between the pressure fields of the FeFEHM compared to another model, highlighting the future potential of the model.

CONCLUSIONS

The initial objective was met, and a new female finite element head model has been developed with biofidelic brain motion. This model will be used for the assessment of repetitive impact scenarios and its repercussions on the female brain.

A Review of Cyclist Head Injury, Impact Characteristics and the Implications for Helmet Assessment Methods

Head injuries are common for cyclists involved in collisions. Such collision scenarios result in a range of injuries, with different head impact speeds, angles, locations, or surfaces. A clear understanding of these collision characteristics is vital to design high fidelity test methods for evaluating the performance of helmets. We review literature detailing real-world cyclist collision scenarios and report on these key characteristics. Our review shows that helmeted cyclists have a considerable reduction in skull fracture and focal brain pathologies compared to non-helmeted cyclists, as well as a reduction in all brain pathologies. The considerable reduction in focal head pathologies is likely to be due to helmet standards mandating thresholds of linear acceleration. The less considerable reduction in diffuse brain injuries is likely to be due to the lack of monitoring head rotation in test methods. We performed a novel meta-analysis of the location of 1809 head impacts from ten studies. Most studies showed that the side and front regions are frequently impacted, with one large, contemporary study highlighting a high proportion of occipital impacts. Helmets frequently had impact locations low down near the rim line. The face is not well protected by most conventional bicycle helmets. Several papers determine head impact speed and angle from in-depth reconstructions and computer simulations. They report head impact speeds from 5 to 16 m/s, with a concentration around 5 to 8 m/s and higher speeds when there was another vehicle involved in the collision. Reported angles range from 10° to 80° to the normal, and are concentrated around 30°-50°. Our review also shows that in nearly 80% of the cases, the head impact is reported to be against a flat surface. This review highlights current gaps in data, and calls for more research and data to better inform improvements in testing methods of standards and rating schemes and raise helmet safety.

Approximating subject-specific brain injury models via scaling based on head-brain morphological relationships

Most human head/brain models represent a generic adult male head/brain. They may suffer in accuracy when investigating traumatic brain injury (TBI) on a subject-specific basis. Subject-specific models can be developed from neuroimages; however, neuroimages are not typically available in practice. In this study, we establish simple and elegant regression models between brain outer surface morphology and head dimensions measured from neuroimages along with age and sex information (N = 191; 141 males and 50 females with age ranging 14-25 years). The regression models are then used to approximate subject-specific brain models by scaling a generic counterpart, without using neuroimages. Model geometrical accuracy is assessed using adjusted [Formula: see text] and absolute percentage error (e.g., 0.720 and 3.09 ± 2.38%, respectively, for brain volume when incorporating tragion-to-top). For a subset of 11 subjects (from smallest to largest in brain volume), impact-induced brain strains are compared with those from "morphed models" derived from neuroimage-based mesh warping. We find that regional peak strains from the scaled subject-specific models are comparable to those of the morphed counterparts but could be considerably different from those of the generic model (e.g., linear regression slope of 1.01-1.03 for gray and white matter regions versus 1.16-1.19, or up to ~ 20% overestimation for the smallest brain studied). These results highlight the importance of incorporating brain morphological variations in impact simulation and demonstrate the feasibility of approximating subject-specific brain models without neuroimages using age, sex, and easily measurable head dimensions. The scaled models may improve subject specificity for future TBI investigations.

Biomechanics of Traumatic Head and Neck Injuries on Women: A State-of-the-Art Review and Future Directions

The biomechanics of traumatic injuries of the human body as a consequence of road crashes, falling, contact sports, and military environments have been studied for decades. In particular, traumatic brain injury (TBI), the so-called "silent epidemic", is the traumatic insult responsible for the greatest percentage of death and disability, justifying the relevance of this research topic. Despite its great importance, only recently have research groups started to seriously consider the sex differences regarding the morphology and physiology of women, which differs from men and may result in a specific outcome for a given traumatic event. This work aims to provide a summary of the contributions given in this field so far, from clinical reports to numerical models, covering not only the direct injuries from inertial loading scenarios but also the role sex plays in the conditions that precede an accident, and post-traumatic events, with an emphasis on neuroendocrine dysfunctions and chronic traumatic encephalopathy. A review on finite element head models and finite element neck models for the study of specific traumatic events is also performed, discussing whether sex was a factor in validating them. Based on the information collected, improvement perspectives and future directions are discussed.

Development of detailed finite element models for in silico analyses of brain impact dynamics

BACKGROUND AND OBJECTIVE

In the last few decades, several studies have been performed to investigate traumatic brain injuries (TBIs) and to understand the biomechanical response of brain tissues, by using experimental and computational approaches. As part of computational approaches, human head finite element (FE) models show to be important tools in the analysis of TBIs, making it possible to estimate local mechanical effects on brain tissue for different accident scenarios. The present study aims to contribute to the computational approach by means of the development of three advanced FE head models for accurately describing the head tissue dynamics, the first step to predict TBIs.

METHODS

We have developed three detailed FE models of human heads from magnetic resonance images of three volunteers: an adult female (32 yrs), an adult male (35 yrs), and a young male (16 yrs). These models have been validated against experimental data of post mortem human subjects (PMHS) tests available in the literature. Brain tissue displacements relative to the skull, hydrostatic intracranial pressure, and head acceleration have been used as the parameters to compare the model response with the experimental response for validation. The software CORAplus (CORrelation and Analysis) has been adopted to evaluate the bio-fidelity level of FE models.

RESULTS

Numerical results from the three models agree with experimental data. FE models presented in this study show a good bio-fidelity for hydrostatic pressure (CORA score of 0.776) and a fair bio-fidelity brain tissue displacements relative to the skull (CORA score of 0.443 and 0.535). The comparison among numerical simulations carried out with the three models shows negligible differences in the mechanical state of brain tissue due to the different morphometry of the heads, when the same acceleration history is considered.

CONCLUSIONS

The three FE models, thanks to their accurate description of anatomical morphology and to their bio-fidelity, can be useful tools to investigate brain mechanics due to different impact scenarios. Therefore, they can be used for different purposes, such as the investigation of the correlation between head acceleration and tissue damage, or the effectiveness of helmet designs. This work does not address the issue to define injury thresholds for the proposed models.

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.

Development and validation of subject-specific 3D human head models based on a nonlinear visco-hyperelastic constitutive framework

Computational head models are promising tools for understanding and predicting traumatic brain injuries. Most available head models are developed using inputs (i.e. head geometry, material properties and boundary conditions) from experiments on cadavers or animals and employ hereditary integral-based constitutive models that assume linear viscoelasticity in part of the rate-sensitive material response. This leads to high uncertainty and poor accuracy in capturing the nonlinear brain tissue response. To resolve these issues, a framework for the development of subject-specific three-dimensional head models is proposed, in which all inputs are derived in vivo from the same living human subject: head geometry via magnetic resonance imaging (MRI), brain tissue properties via magnetic resonance elastography (MRE), and full-field strain-response of the brain under rapid head rotation via tagged MRI. A nonlinear, viscous dissipation-based visco-hyperelastic constitutive model is employed to capture brain tissue response. Head models are validated using quantitative metrics that compare spatial strain distribution, temporal strain evolution, and the magnitude of strain maxima, with the corresponding experimental observations from tagged MRI. Results show that our head models accurately capture the strain-response of the brain. Further, employment of the nonlinear visco-hyperelastic constitutive framework provides improvements in the prediction of peak strains and temporal strain evolution over hereditary integral-based models.

Concussion Prone Scenarios: A Multi-Dimensional Exploration in Impact Directions, Brain Morphology, and Network Architectures Using Computational Models

While individual susceptibility to traumatic brain injury (TBI) has been speculated, past work does not provide an analysis considering how physical features of an individual's brain (e.g., brain size, shape), impact direction, and brain network features can holistically contribute to the risk of suffering a TBI from an impact. This work investigated each of these features simultaneously using computational modeling and analyses of simulated functional connectivity. Unlike the past studies that assess the severity of TBI based on the quantification of brain tissue damage (e.g., principal strain), we approached the brain as a complex network in which neuronal oscillations orchestrate to produce normal brain function (estimated by functional connectivity) and, to this end, both the anatomical damage location and its topological characteristics within the brain network contribute to the severity of brain function disruption and injury. To represent the variations in the population, we analyzed a publicly available database of brain imaging data and selected five distinct network architectures, seven different brain sizes, and three uniaxial head rotational conditions to study the consequences of 74 virtual impact scenarios. Results show impact direction produces the most significant change in connections across brain areas (structural connectome) and the functional coupling of activity across these brain areas (functional connectivity). Axial rotations were more injurious than those with sagittal and coronal rotations when the head kinematics were the same for each condition. When the impact direction was held constant, brain network architecture showed a significantly different vulnerability across axial and sagittal, but not coronal rotations. As expected, brain size significantly affected the expected change in structural and functional connectivity after impact. Together, these results provided groupings of predicted vulnerability to impact-a subgroup of male brain architectures exposed to axial impacts were most vulnerable, while a subgroup of female brain architectures was the most tolerant to the sagittal impacts studied. These findings lay essential groundwork for subject-specific analyses of concussion and provide invaluable guidance for designing personalized protection equipment.

Brain architecture-based vulnerability to traumatic injury

The white matter tracts forming the intricate wiring of the brain are subject-specific; this heterogeneity can complicate studies of brain function and disease. Here we collapse tractography data from the Human Connectome Project (HCP) into structural connectivity (SC) matrices and identify groups of similarly wired brains from both sexes. To characterize the significance of these architectural groupings, we examined how similarly wired brains led to distinct groupings of neural activity dynamics estimated with Kuramoto oscillator models (KMs). We then lesioned our networks to simulate traumatic brain injury (TBI) and finally we tested whether these distinct architecture groups’ dynamics exhibited differing responses to simulated TBI. At each of these levels we found that brain structure, simulated dynamics, and injury susceptibility were all related to brain grouping. We found four primary brain architecture groupings (two male and two female), with similar architectures appearing across both sexes. Among these groupings of brain structure, two architecture types were significantly more vulnerable than the remaining two architecture types to lesions. These groups suggest that mesoscale brain architecture types exist, and these architectural differences may contribute to differential risks to TBI and clinical outcomes across the population.

Assessment of brain injury characterization and influence of modeling approaches

In this study, using computational biomechanics models, we investigated influence of the skull-brain interface modeling approach and the material property of cerebrum on the kinetic, kinematic and injury outputs. Live animal head impact tests of different severities were reconstructed in finite element simulations and DAI and ASDH injury results were compared. We used the head/brain models of Total HUman Model for Safety (THUMS) and Global Human Body Models Consortium (GHBMC), which had been validated under several loading conditions. Four modeling approaches of the skull-brain interface in the head/brain models were evaluated. They were the original models from THUMS and GHBMC, the THUMS model with skull-brain interface changed to sliding contact, and the THUMS model with increased shear modulus of cerebrum, respectively. The results have shown that the definition of skull-brain interface would significantly influence the magnitude and distribution of the load transmitted to the brain. With sliding brain-skull interface, the brain had lower maximum principal stress compared to that with strong connected interface, while the maximum principal strain slightly increased. In addition, greater shear modulus resulted in slightly higher the maximum principal stress and significantly lower the maximum principal strain. This study has revealed that using models with different modeling approaches, the same value of injury metric may correspond to different injury severity.

Cumulative strain-based metrics for predicting subconcussive head impact exposure-related imaging changes in a cohort of American youth football players

OBJECTIVE

Youth football athletes are exposed to repetitive subconcussive head impacts during normal participation in the sport, and there is increasing concern about the long-term effects of these impacts. The objective of the current study was to determine if strain-based cumulative exposure measures are superior to kinematic-based exposure measures for predicting imaging changes in the brain.

METHODS

This prospective, longitudinal cohort study was conducted from 2012 to 2017 and assessed youth, male football athletes. Kinematic data were collected at all practices and games from enrolled athletes participating in local youth football organizations in Winston-Salem, North Carolina, and were used to calculate multiple risk-weighted cumulative exposure (RWE) kinematic metrics and 36 strain-based exposure metrics. Pre- and postseason imaging was performed at Wake Forest School of Medicine, and diffusion tensor imaging (DTI) measures, including fractional anisotropy (FA), and its components (CL, CP, and CS), and mean diffusivity (MD), were investigated. Included participants were youth football players ranging in age from 9 to 13 years. Exclusion criteria included any history of previous neurological illness, psychiatric illness, brain tumor, concussion within the past 6 months, and/or contraindication to MRI.

RESULTS

A total of 95 male athletes (mean age 11.9 years [SD 1.0 years]) participated between 2012 and 2017, with some participating for multiple seasons, resulting in 116 unique athlete-seasons. Regression analysis revealed statistically significant linear relationships between the FA, linear coefficient (CL), and spherical coefficient (CS) and all strain exposure measures, and well as the planar coefficient (CP) and 8 strain measures. For the kinematic exposure measures, there were statistically significant relationships between FA and RWE linear (RWEL) and RWE combined probability (RWECP) as well as CS and RWEL. According to area under the receiver operating characteristic (ROC) curve (AUC) analysis, the best-performing metrics were all strain measures, and included metrics based on tensile, compressive, and shear strain.

CONCLUSIONS

Using ROC curves and AUC analysis, all exposure metrics were ranked in order of performance, and the results demonstrated that all the strain-based metrics performed better than any of the kinematic metrics, indicating that strain-based metrics are better discriminators of imaging changes than kinematic-based measures. Studies relating the biomechanics of head impacts with brain imaging and cognitive function may allow equipment designers, care providers, and organizations to prevent, identify, and treat injuries in order to make football a safer activity.

Physics-Informed Machine Learning Improves Detection of Head Impacts

In this work we present a new physics-informed machine learning model that can be used to analyze kinematic data from an instrumented mouthguard and detect impacts to the head. Monitoring player impacts is vitally important to understanding and protecting from injuries like concussion. Typically, to analyze this data, a combination of video analysis and sensor data is used to ascertain the recorded events are true impacts and not false positives. In fact, due to the nature of using wearable devices in sports, false positives vastly outnumber the true positives. Yet, manual video analysis is time-consuming. This imbalance leads traditional machine learning approaches to exhibit poor performance in both detecting true positives and preventing false negatives. Here, we show that by simulating head impacts numerically using a standard Finite Element head-neck model, a large dataset of synthetic impacts can be created to augment the gathered, verified, impact data from mouthguards. This combined physics-informed machine learning impact detector reported improved performance on test datasets compared to traditional impact detectors with negative predictive value and positive predictive values of 88 and 87% respectively. Consequently, this model reported the best results to date for an impact detection algorithm for American football, achieving an F₁ score of 0.95. In addition, this physics-informed machine learning impact detector was able to accurately detect true and false impacts from a test dataset at a rate of 90% and 100% relative to a purely manual video analysis workflow. Saving over 12 h of manual video analysis for a modest dataset, at an overall accuracy of 92%, these results indicate that this model could be used in place of, or alongside, traditional video analysis to allow for larger scale and more efficient impact detection in sports such as American Football.

Finite element modeling of the complex anisotropic mechanical behavior of the human sclera and pia mater

BACKGROUND AND OBJECTIVE

Accurate finite element (FE) simulation of the optic nerve head (ONH) depends on accurate mechanical properties of the load-bearing tissues. The peripapillary sclera in the ONH exhibits a depth-dependent, anisotropic, heterogeneous collagen fiber distribution. This study proposes a novel cable-in-solid modeling approach that mimics heterogeneous anisotropic collagen fiber distribution, validates the approach against published experimental biaxial tensile tests of scleral patches, and demonstrates its effectiveness in a complex model of the posterior human eye and ONH.

METHODS

A computational pipeline was developed that defines control points in the sclera and pia mater, distributes the depth-dependent circumferential, radial, and isotropic cable elements in the sclera and pia in a pattern that mimics collagen fiber orientation, and couples the cable elements and solid matrix using a mesh-free penalty-based cable-in-solid algorithm. A parameter study was performed on a model of a human scleral patch subjected to biaxial deformation, and computational results were matched to published experimental data. The new approach was incorporated into a previously published eye-specific model to test the method; results were then interpreted in relation to the collagen fibers' (cable elements) role in the resultant ONH deformations, stresses, and strains.

RESULTS

Results show that the cable-in-solid approach can mimic the full range of scleral mechanical behavior measured experimentally. Disregarding the collagen fibers/cable elements in the posterior eye model resulted in ∼20-60% greater tensile and shear stresses and strains, and ∼30% larger posterior deformations in the lamina cribrosa and peripapillary sclera.

CONCLUSIONS

The cable-in-solid approach can easily be implemented into commercial FE packages to simulate the heterogeneous and anisotropic mechanical properties of collagenous biological tissues.

Displacement voxelization to resolve mesh-image mismatch: Application in deriving dense white matter fiber strains

BACKGROUND AND OBJECTIVE

It is common to combine biomechanical modeling and medical images for multimodal analyses. However, mesh-image mismatch may occur that prevents direct information exchange. To eliminate mesh-image mismatch, we develop a simple but elegant displacement voxelization technique based on image voxel corner nodes to achieve voxel-wise strain. We then apply the technique to derive dense white matter fiber strains along whole-brain tractography (∼35 k fiber tracts consisting of ∼3.3 million sampling points) resulting from head impact.

METHODS

Displacements at image voxel corner nodes are first obtained from model simulation via scattered interpolation. Each voxel is then scaled linearly to form a unit hexahedral element. This allows convenient and efficient voxel-wise strain tensor calculation and displacement interpolation at arbitrary fiber sampling points via shape functions. Fiber strains from displacement interpolation are then compared with those from the commonly used strain tensor projection using either voxel- or element-wise strain tensors.

RESULTS

Based on a synthetic displacement field, fiber strains interpolated from voxelized displacement are considerably more accurate than those from strain tensor projection relative to the prescribed ground-truth (determinant of coefficient (R²) of 1.00 and root mean squared error (RMSE) of 0.01 vs. 0.87 and 0.10, respectively). For a set of real-world reconstructed head impacts (N = 53), the strain tensor projection method performs similarly poorly (R² of 0.80-0.90 and RMSE of 0.03-0.07), with overestimation strongly correlated with strain magnitude (Pearson correlation coefficient >0.9). Up to ∼15% of the fiber strains are overestimated by more than the lower bound of a conservative injury threshold of 0.09. The percentage increases to ∼37% when halving the threshold. Voxel interpolation is also significantly more efficient (15 s vs. 40 s for element strain tensor projection, without parallelization).

CONCLUSIONS

Voxelized displacement interpolation is considerably more accurate and efficient in deriving dense white matter fiber strains than strain tensor projection. The latter generally overestimates with overestimation magnitude strongly correlating with fiber strain magnitude. Displacement voxelization is an effective technique to eliminate mesh-image mismatch and generates a convenient image representation of tissue deformation. This technique can be generalized to broadly facilitate a diverse range of image-related biomechanical problems for multimodal analyses. The convenient image format may also promote and facilitate biomechanical data sharing in the future.

Integrating Human and Non-Human Primate Data to Estimate Human Tolerances for Traumatic Brain Injury

Traumatic brain injury (TBI) contributes to a significant portion of the injuries resulting from motor vehicle crashes, falls, and sports collisions. The development of advanced countermeasures to mitigate these injuries requires a complete understanding of the tolerance of the human brain to injury. In this study, we developed a new method to establish human injury tolerance levels using an integrated database of reconstructed football impacts, sub-injurious human volunteer data, and non-human primate data. The human tolerance levels were analyzed using tissue-level metrics determined using harmonized species-specific finite element brain models. Kinematics-based metrics involving complete characterization of angular motion (e.g., DAMAGE) showed better power of predicting tissue-level deformation in a variety of impact conditions and were subsequently used to characterize injury tolerance. The proposed human brain tolerances for mild and severe TBI were estimated and presented in the form of injury risk curves based on selected tissue-level and kinematics-based injury metrics. The application of the estimated injury tolerances was finally demonstrated using real-world automotive crash data.

Relating strain fields with microtubule changes in porcine cortical sulci following drop impact

The biomechanical response of brain tissue to strain and the immediate neural outcomes are of fundamental importance in understanding mild traumatic brain injury (mTBI). The sensitivity of neural tissue to dynamic strain events and the resulting strain-induced changes are considered to be a primary factor in injury. Rodent models have been used extensively to investigate impact-induced injury. However, the lissencephalic structure is inconsistent with the human brain, which is gyrencephalic (convoluted structure), and differs considerably in strain field localization effects. Porcine brains have a similar structure to the human brain, containing a similar ratio of white-grey matter and gyrification in the cortex. In this study, coronal brain slabs were extracted from female pig brains within 2hrs of sacrifice. Slabs were implanted with neutral density radiopaque markers, sealed inside an elastomeric encasement, and dropped from 0.9 m onto a steel anvil. Particle tracking revealed elevated tensile strains in the sulcus. One hour after impact, decreased microtubule associated protein 2 (MAP2) was found exclusively within the sulcus with no increase in cell death. These results suggest that elevated tensile strain in the sulcus may result in compromised cytoskeleton, possibly indicating a vulnerability to pathological outcomes under the right circumstances. The results demonstrated that the observed changes were unrelated to shear strain loading of the tissues but were more sensitive to tensile load.

Instantaneous Brain Strain Estimation for Automotive Head Impacts via Deep Learning

Efficient brain strain estimation is critical for routine application of a head injury model. Lately, a convolutional neural network (CNN) has been successfully developed to estimate spatially detailed brain strains instantly and accurately in contact sports. Here, we extend its application to automotive head impacts, where impact profiles are typically more complex with longer durations. Head impact kinematics (N=458) from two public databases were used to generate augmented impacts (N=2694). They were simulated using the anisotropic Worcester Head Injury Model (WHIM) V1.0, which provided baseline elementwise peak maximum principal strain (MPS). For each augmented impact, rotational velocity (v_rot) and the corresponding rotational acceleration (a_rot) profiles were concatenated as static images to serve as CNN input. Three training strategies were evaluated: 1) "baseline", using random initial weights; 2) "transfer learning", using weight transfer from a previous CNN model trained on head impacts drawn from contact sports; and 3) "combined training", combining previous training data from contact sports (N=5661) for training. The combined training achieved the best performances. For peak MPS, the CNN achieved a coefficient of determination (R²) of 0.932 and root mean squared error (RMSE) of 0.031 for the real-world testing dataset. It also achieved a success rate of 60.5% and 94.8% for elementwise MPS, where the linear regression slope, k, and correlation coefficient, r, between estimated and simulated MPS did not deviate from 1.0 (when identical) by more than 0.1 and 0.2, respectively. Cumulative strain damage measure (CSDM) from the CNN estimation was also highly accurate compared to those from direct simulation across a range of thresholds (R² of 0.899-0.943 with RMSE of 0.054-0.069). Finally, the CNN achieved an average k and r of 0.98±0.12 and 0.90±0.07, respectively, for six reconstructed car crash impacts drawn from two other sources independent of the training dataset. Importantly, the CNN is able to efficiently estimate elementwise MPS with sufficient accuracy while conventional kinematic injury metrics cannot. Therefore, the CNN has the potential to supersede current kinematic injury metrics that can only approximate a global peak MPS or CSDM. The CNN technique developed here may offer enhanced utility in the design and development of head protective countermeasures, including in the automotive industry. This is the first study aimed at instantly estimating spatially detailed brain strains for automotive head impacts, which employs >8.8 thousand impact simulations generated from ~1.5 years of nonstop computations on a high-performance computing platform.

Mechanical characterisation of the human dura mater, falx cerebri and superior sagittal sinus

The cranial meninges have been shown to play a pivotal role in traumatic brain injury mechanopathology. However, while the mechanical response of the brain and its many subregions have been studied extensively, the meninges have conventionally been overlooked. This paper presents the first comparative mechanical analysis of human dura mater, falx cerebri and superior sagittal sinus tissues. Biaxial tensile analysis identified that these tissues are mechanically heterogeneous, in contrast to the assumption that the tissues are mechanically homogeneous which is typically employed in FE model design. A thickness of 0.91 ± 0.05 (standard error) mm for the falx cerebri was also identified. This data can aid in improving the biofidelity of the influential falx structure in FE models. Additionally, the use of a collagen hybridizing peptide on the superior sagittal sinus suggests this structure is particularly susceptible to the effects of circumferential stretch, which may have important implications for clinical treatment of dural venous sinus pathologies. Collectively, this research progresses understanding of meningeal mechanical and structural characteristics and may aid in elucidating the behaviour of these tissues in healthy and diseased conditions. STATEMENT OF SIGNIFICANCE: This study presents the first evaluation of human falx cerebri and superior sagittal sinus mechanical, geometrical and structural properties, along with a comparison to cranial dura mater. To mechanically characterise the tissues, biaxial tensile testing is conducted on the tissues. This analysis identifies, for the first time, mechanical stiffness differences between these tissues. Additionally, geometrical analysis identifies that there are thickness differences between the tissues. The evaluation of human meningeal tissues allows for direct implementation of the novel data to finite element head injury models to enable improved biofidelity of these influential structures in traumatic brain injury simulations. This work also identifies that the superior sagittal sinus may be easily damaged during clinical angioplasty procedures, which may inform the treatment of dural sinus pathologies.

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.

Investigation of the Compressive Viscoelastic Properties of Brain Tissue Under Time and Frequency Dependent Loading Conditions

The mechanical characterization of brain tissue has been generally analyzed in the frequency and time domain. It is crucial to understand the mechanics of the brain under realistic, dynamic conditions and convert it to enable mathematical modelling in a time domain. In this study, the compressive viscoelastic properties of brain tissue were investigated under time and frequency domains with the same physical conditions and the theory of viscoelasticity was applied to estimate the prediction of viscoelastic response in the time domain based on frequency-dependent mechanical moduli through Finite Element models. Storage and loss modulus were obtained from white and grey matter, of bovine brains, using dynamic mechanical analysis and time domain material functions were derived based on a Prony series representation. The material models were evaluated using brain testing data from stress relaxation and hysteresis in the time dependent analysis. The Finite Element models were able to represent the trend of viscoelastic characterization of brain tissue under both testing domains. The outcomes of this study contribute to a better understanding of brain tissue mechanical behaviour and demonstrate the feasibility of deriving time-domain viscoelastic parameters from frequency-dependent compressive data for biological tissue, as validated by comparing experimental tests with computational simulations.

MR Imaging of Human Brain Mechanics In Vivo: New Measurements to Facilitate the Development of Computational Models of Brain Injury

Computational models of the brain and its biomechanical response to skull accelerations are important tools for understanding and predicting traumatic brain injuries (TBIs). However, most models have been developed using experimental data collected on animal models and cadaveric specimens, both of which differ from the living human brain. Here we describe efforts to noninvasively measure the biomechanical response of the human brain with MRI-at non-injurious strain levels-and generate data that can be used to develop, calibrate, and evaluate computational brain biomechanics models. Specifically, this paper reports on a project supported by the National Institute of Neurological Disorders and Stroke to comprehensively image brain anatomy and geometry, mechanical properties, and brain deformations that arise from impulsive and harmonic skull loadings. The outcome of this work will be a publicly available dataset ( http://www.nitrc.org/projects/bbir ) that includes measurements on both males and females across an age range from adolescence to older adulthood. This article describes the rationale and approach for this study, the data available, and how these data may be used to develop new computational models and augment existing approaches; it will serve as a reference to researchers interested in using these data.

Displacement Error Propagation From Embedded Markers to Brain Strain

Head injury model validation has evolved from against pressure to relative brain-skull displacement, and more recently, against marker-based strain. However, there are concerns on strain data quality. In this study, we parametrically investigate how displacement random errors and synchronization errors propagate into strain. Embedded markers from four representative configurations are used to form unique and nonoverlapping tetrahedrons, triangles, and linear elements. Marker displacements are then separately subjected to up to ±10% random displacement errors and up to ±2 ms synchronization errors. Based on 100 random trials in each perturbation test, we find that smaller strain errors relative to the baseline peak strains are significantly associated with larger element sizes (volume, area, or length; p < 0.05). When displacement errors are capped at the two extreme levels, the earlier "column" and "cluster" configurations provide few usable elements with relative strain error under an empirical threshold of 20%, while about 30-80% of elements in recent "repeatable" and "uniform" configurations are considered otherwise usable. Overall, denser markers are desired to provide exhaustive pairwise linear elements with a range of sizes to balance the need for larger elements to minimize strain error but smaller elements to increase the spatial resolution in strain sampling. Their signed strains also provide unique and unambiguous information on tissue tension and compression. This study may provide useful insights into the scrutinization of existing experimental data for head injury model strain validation and to inform how best to design new experiments in the future.

Toward subject-specific evaluation: methods of evaluating finite element brain models using experimental high-rate rotational brain motion

Computational models of the brain have become the gold standard in biomechanics to understand, predict, and mitigate traumatic brain injuries. Many models have been created and evaluated with limited experimental data and without accounting for subject-specific morphometry of the specimens in the dataset. Recent advancements in the measurement of brain motion using sonomicrometry allow for a comprehensive evaluation of brain model biofidelity using a high-rate, rotational brain motion dataset. In this study, four methods were used to determine the best technique to compare nodal displacement to experimental brain motion, including a new morphing method to match subject-specific inner skull geometry. Three finite element brain models were evaluated in this study: the isotropic GHBMC and SIMon models, as well as an anisotropic model with explicitly embedded axons (UVA-EAM). Using a weighted cross-correlation score (between 0 and 1), the anisotropic model yielded the highest average scores across specimens and loading conditions ranging from 0.53 to 0.63, followed by the isotropic GHBMC with average scores ranging from 0.46 to 0.58, and then the SIMon model with average scores ranging from 0.36 to 0.51. The choice of comparison method did not significantly affect the cross-correlation score, and differences of global strain up to 0.1 were found for the morphed geometry relative to baseline models. The morphed or scaled geometry is recommended when evaluating computational brain models to capture the subject-specific skull geometry of the experimental specimens.

Evaluation of Tissue-Level Brain Injury Metrics Using Species-Specific Simulations

Traumatic brain injury (TBI) is a significant public health burden, and the development of advanced countermeasures to mitigate and prevent these injuries during automotive, sports, and military impact events requires an understanding of the intracranial mechanisms related to TBI. In this study, the efficacy of tissue-level injury metrics for predicting TBI was evaluated using finite element reconstructions from a comprehensive, multi-species TBI database. The database consisted of human volunteer tests, laboratory-reconstructed head impacts from sports, in vivo non-human primate (NHP) tests, and in vivo pig tests. Eight tissue-level metrics related to brain tissue strain, axonal strain, and strain-rate were evaluated using survival analysis for predicting mild and severe TBI risk. The correlation between TBI risk and most of the assessed metrics were statistically significant, but when injury data was analyzed by species, the best metric was often inconclusive and limited by the small datasets. When the human and animal datasets were combined, the injury analysis was able to delineate maximum axonal strain as the best predictor of injury for all species and TBI severities, with maximum principal strain as a suitable alternative metric. The current study is the first to provide evidence to support the assumption that brain strain response between human, pig, and NHP result in similar injury outcomes through a multi-species analysis. This assumption is the biomechanical foundation for translating animal brain injury findings to humans. The findings in the study provide fundamental guidelines for developing injury criteria that would contribute towards the innovation of more effective safety countermeasures.

Calibration of a Heterogeneous Brain Model Using a Subject-Specific Inverse Finite Element Approach

Central to the investigation of the biomechanics of traumatic brain injury (TBI) and the assessment of injury risk from head impact are finite element (FE) models of the human brain. However, many existing FE human brain models have been developed with simplified representations of the parenchyma, which may limit their applicability as an injury prediction tool. Recent advances in neuroimaging techniques and brain biomechanics provide new and necessary experimental data that can improve the biofidelity of FE brain models. In this study, the CAB-20MSym template model was developed, calibrated, and extensively verified. To implement material heterogeneity, a magnetic resonance elastography (MRE) template image was leveraged to define the relative stiffness gradient of the brain model. A multi-stage inverse FE (iFE) approach was used to calibrate the material parameters that defined the underlying non-linear deviatoric response by minimizing the error between model-predicted brain displacements and experimental displacement data. This process involved calibrating the infinitesimal shear modulus of the material using low-severity, low-deformation impact cases and the material non-linearity using high-severity, high-deformation cases from a dataset of in situ brain displacements obtained from cadaveric specimens. To minimize the geometric discrepancy between the FE models used in the iFE calibration and the cadaveric specimens from which the experimental data were obtained, subject-specific models of these cadaveric brain specimens were developed and used in the calibration process. Finally, the calibrated material parameters were extensively verified using independent brain displacement data from 33 rotational head impacts, spanning multiple loading directions (sagittal, coronal, axial), magnitudes (20–40 rad/s), durations (30–60 ms), and severity. Overall, the heterogeneous CAB-20MSym template model demonstrated good biofidelity with a mean overall CORA score of 0.63 ± 0.06 when compared to in situ brain displacement data. Strains predicted by the calibrated model under non-injurious rotational impacts in human volunteers (N = 6) also demonstrated similar biofidelity compared to in vivo measurements obtained from tagged magnetic resonance imaging studies. In addition to serving as an anatomically accurate model for further investigations of TBI biomechanics, the MRE-based framework for implementing material heterogeneity could serve as a foundation for incorporating subject-specific material properties in future models.

Instantaneous Whole-Brain Strain Estimation in Dynamic Head Impact

Head injury models are notoriously time consuming and resource demanding in simulations, which prevents routine application. Here, we extend a convolutional neural network (CNN) to instantly estimate element-wise distribution of peak maximum principal strain (MPS) of the entire brain (>36 k speedup accomplished on a low-end computing platform). To achieve this, head impact rotational velocity and acceleration temporal profiles are combined into two-dimensional images to serve as CNN input for training and prediction of MPS. Compared with the directly simulated counterparts, the CNN-estimated responses (magnitude and distribution) are sufficiently accurate for 92.1% of the cases via 10-fold cross-validation using impacts drawn from the real world (n = 5661; range of peak rotational velocity in augmented data extended to 2-40 rad/sec). The success rate further improves to 97.1% for "in-range" impacts (n = 4298). When using the same CNN architecture to train (n = 3064) and test on an independent, reconstructed National Football League (NFL) impact dataset (n = 53; 20 concussions and 33 non-injuries), 51 out of 53, or 96.2% of the cases, are sufficiently accurate. The estimated responses also achieve virtually identical concussion prediction performances relative to the directly simulated counterparts, and they often outperform peak MPS of the whole brain (e.g., accuracy of 0.83 vs. 0.77 via leave-one-out cross-validation). These findings support the use of CNN for accurate and efficient estimation of spatially detailed brain strains across the vast majority of head impacts in contact sports. Our technique may hold the potential to transform traumatic brain injury (TBI) research and the design and testing standards of head protective gears by facilitating the transition from acceleration-based approximation to strain-based design and analysis. This would have broad implications in the TBI biomechanics field to accelerate new scientific discoveries. The pre-trained CNN is freely available online at https://github.com/Jilab-biomechanics/CNN-brain-strains.

An anatomically detailed and personalizable head injury model: Significance of brain and white matter tract morphological variability on strain

Finite element head (FE) models are important numerical tools to study head injuries and develop protection systems. The generation of anatomically accurate and subject-specific head models with conforming hexahedral meshes remains a significant challenge. The focus of this study is to present two developmental works: first, an anatomically detailed FE head model with conforming hexahedral meshes that has smooth interfaces between the brain and the cerebrospinal fluid, embedded with white matter (WM) fiber tracts; second, a morphing approach for subject-specific head model generation via a new hierarchical image registration pipeline integrating Demons and Dramms deformable registration algorithms. The performance of the head model is evaluated by comparing model predictions with experimental data of brain-skull relative motion, brain strain, and intracranial pressure. To demonstrate the applicability of the head model and the pipeline, six subject-specific head models of largely varying intracranial volume and shape are generated, incorporated with subject-specific WM fiber tracts. DICE similarity coefficients for cranial, brain mask, local brain regions, and lateral ventricles are calculated to evaluate personalization accuracy, demonstrating the efficiency of the pipeline in generating detailed subject-specific head models achieving satisfactory element quality without further mesh repairing. The six head models are then subjected to the same concussive loading to study the sensitivity of brain strain to inter-subject variability of the brain and WM fiber morphology. The simulation results show significant differences in maximum principal strain and axonal strain in local brain regions (one-way ANOVA test, p < 0.001), as well as their locations also vary among the subjects, demonstrating the need to further investigate the significance of subject-specific models. The techniques developed in this study may contribute to better evaluation of individual brain injury and the development of individualized head protection systems in the future. This study also contains general aspects the research community may find useful: on the use of experimental brain strain close to or at injury level for head model validation; the hierarchical image registration pipeline can be used to morph other head models, such as smoothed-voxel models.

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.

White Matter Tract-Oriented Deformation Is Dependent on Real-Time Axonal Fiber Orientation

Traumatic axonal injury (TAI) is a critical public health issue with its pathogenesis remaining largely elusive. Finite element (FE) head models are promising tools to bridge the gap between mechanical insult, localized brain response, and resultant injury. In particular, the FE-derived deformation along the direction of white matter (WM) tracts (i.e., tract-oriented strain) has been shown to be an appropriate predictor for TAI. The evolution of fiber orientation in time during the impact and its potential influence on the tract-oriented strain remains unknown, however. To address this question, the present study leveraged an embedded element approach to track real-time fiber orientation during impacts. A new scheme to calculate the tract-oriented strain was proposed by projecting the strain tensors from pre-computed simulations along the temporal fiber direction instead of its static counterpart directly obtained from diffuse tensor imaging. The results revealed that incorporating the real-time fiber orientation not only altered the direction but also amplified the magnitude of the tract-oriented strain, resulting in a generally more extended distribution and a larger volume ratio of WM exposed to high deformation along fiber tracts. These effects were exacerbated with the impact severities characterized by the acceleration magnitudes. Results of this study provide insights into how best to incorporate fiber orientation in head injury models and derive the WM tract-oriented deformation from computational simulations, which is important for furthering our understanding of the underlying mechanisms of TAI.

Integrating material properties from magnetic resonance elastography into subject-specific computational models for the human brain

Advances in brain imaging and computational methods have facilitated the creation of subject-specific computational brain models that aid researchers in investigating brain trauma using simulated impacts. The emergence of magnetic resonance elastography (MRE) as a non-invasive mechanical neuroimaging tool has enabled in vivo estimation of material properties at low-strain, harmonic loading. An open question in the field has been how this data can be integrated into computational models. The goals of this study were to use a novel MRI dataset acquired in human volunteers to generate models with subject-specific anatomy and material properties, and then to compare simulated brain deformations to subject-specific brain deformation data under non-injurious loading. Models of five subjects were simulated with linear viscoelastic (LVE) material properties estimated directly from MRE data. Model predictions were compared to experimental brain deformation acquired in the same subjects using tagged MRI. Outcomes from the models matched the spatial distribution and magnitude of the measured peak strain components as well as the 95^th percentile in-plane peak strains within 0.005 mm/mm and maximum principal strain within 0.012 mm/mm. Sensitivity to material heterogeneity was also investigated. Simulated brain deformations from a model with homogenous brain properties and a model with brain properties discretized with up to ten regions were very similar (a mean absolute difference less than 0.0015 mm/mm in peak strains). Incorporating material properties directly from MRE into a biofidelic subject-specific model is an important step toward future investigations of higher-order model features and simulations under more severe loading conditions.

Variations in User Implementation of the CORA Rating Metric

The CORA rating metric is frequently used in the field of injury biomechanics to compare the similarity of response time histories. However, subjectivity exists within the CORA metric in the form of user-customizable parameters that give the metric the flexibility to be used for a variety of applications. How these parameters are customized is not always reported in the literature, and it is unknown how these customizations affect the CORA scores. Therefore, the purpose of this study was to evaluate how variations in the CORA parameters affect the resulting similarity scores. A literature review was conducted to determine how the CORA parameters are commonly customized within the literature. Then, CORA scores for two datasets were calculated using the most common parameter customizations and the default parameters. Differences between the CORA scores using customized and default parameters were statistically significant for all customizations. Furthermore, most customizations produced score increases relative to the default settings. The use of standard deviation corridors and exclusion of the corridor component were found to produce the largest score differences. The observed differences demonstrated the need for researchers to exercise transparency when using customized parameters in CORA analyses.

A Review of Validation Methods for the Intracranial Response of FEHM to Blunt Impacts

The following is a review of the processes currently employed when validating the intracranial response of Finite Element Head Models (FEHM) against blunt impacts. The authors aim to collate existing validation tools, their applications and findings on their effectiveness to aid researchers in the validation of future FEHM and potential efforts in improving procedures. In this vain, publications providing experimental data on the intracranial pressure, relative brain displacement and brain strain responses to impacts in human subjects are surveyed and key data are summarised. This includes cases that have previously been used in FEHM validation and alternatives with similar potential uses. The processes employed to replicate impact conditions and the resulting head motion are reviewed, as are the analytical techniques used to judge the validity of the models. Finally, publications exploring the validation process and factors affecting it are critically discussed. Reviewing FEHM validation in this way highlights the lack of a single best practice, or an obvious solution to create one using the tools currently available. There is clear scope to improve the validation process of FEHM, and the data available to achieve this. By collecting information from existing publications, it is hoped this review can help guide such developments and provide a point of reference for researchers looking to validate or investigate FEHM in the future, enabling them to make informed choices about the simulation of impacts, how they are generated numerically and the factors considered during output assessment, whilst being aware of potential limitations in the process.

An anatomically detailed and personalizable head injury model: Significance of brain and white matter tract morphological variability on strain

Finite element head (FE) models are important numerical tools to study head injuries and develop protection systems. The generation of anatomically accurate and subject-specific head models with conforming hexahedral meshes remains a significant challenge. The focus of this study is to present two developmental works: first, an anatomically detailed FE head model with conforming hexahedral meshes that has smooth interfaces between the brain and the cerebrospinal fluid, embedded with white matter (WM) fiber tracts; second, a morphing approach for subject-specific head model generation via a new hierarchical image registration pipeline integrating Demons and Dramms deformable registration algorithms. The performance of the head model is evaluated by comparing model predictions with experimental data of brain-skull relative motion, brain strain, and intracranial pressure. To demonstrate the applicability of the head model and the pipeline, six subject-specific head models of largely varying intracranial volume and shape are generated, incorporated with subject-specific WM fiber tracts. DICE similarity coefficients for cranial, brain mask, local brain regions, and lateral ventricles are calculated to evaluate personalization accuracy, demonstrating the efficiency of the pipeline in generating detailed subject-specific head models achieving satisfactory element quality without further mesh repairing. The six head models are then subjected to the same concussive loading to study the sensitivity of brain strain to inter-subject variability of the brain and WM fiber morphology. The simulation results show significant differences in maximum principal strain and axonal strain in local brain regions (one-way ANOVA test, p < 0.001), as well as their locations also vary among the subjects, demonstrating the need to further investigate the significance of subject-specific models. The techniques developed in this study may contribute to better evaluation of individual brain injury and the development of individualized head protection systems in the future. This study also contains general aspects the research community may find useful: on the use of experimental brain strain close to or at injury level for head model validation; the hierarchical image registration pipeline can be used to morph other head models, such as smoothed-voxel models.

Structural Anisotropy vs. Mechanical Anisotropy: The Contribution of Axonal Fibers to the Material Properties of Brain White Matter

Brain's micro-structure plays a critical role in its macro-structure material properties. Since the structural anisotropy in the brain white matter has been introduced due to axonal fibers, considering the direction of axons in the continuum models has been mediated to improve the results of computational simulations. The aim of the current study was to investigate the role of fiber direction in the material properties of brain white matter and compare the mechanical behavior of the anisotropic white matter and the isotropic gray matter. Diffusion tensor imaging (DTI) was employed to detect the direction of axons in white matter samples, and tensile stress-relaxation loads up to 20% strains were applied on bovine gray and white matter samples. In order to calculate the nonlinear and time-dependent properties of white matter and gray matter, a visco-hyperelastic model was used. The results indicated that the mechanical behavior of white matter in two orthogonal directions, parallel and perpendicular to axonal fibers, are significantly different. This difference indicates that brain white matter could be assumed as an anisotropic material and axons have contribution in the mechanical properties. Also, up to 15% strain, white matter samples with axons parallel to the force direction are significantly stiffer than both the gray matter samples and white matter samples with axons perpendicular to the force direction. Moreover, the elastic moduli of white matter samples with axons both parallel and perpendicular to the loading direction and gray matter samples at 15-20% strain are not significantly different. According to these observations, it is suggested that axons have negligible roles in the material properties of white matter when it is loaded in the direction perpendicular to the axon direction. Finally, this observation showed that the anisotropy of brain tissue not only has effects on the elastic behavior, but also has effects on the viscoelastic behavior.