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

A Spatial Interpolation Approach to Assign Magnetic Resonance Imaging-Derived Material Properties for Finite Element Models of Adeno-Associated Virus Infusion Into a Recurrent Brain Tumor

Adeno-associated virus (AAV) is a clinically useful gene delivery vehicle for treating neurological diseases. To deliver AAV to focal targets, direct infusion into brain tissue by convection-enhanced delivery (CED) is often needed due to AAV's limited penetration across the blood-brain-barrier and its low diffusivity in tissue. In this study, computational models that predict the spatial distribution of AAV in brain tissue during CED were developed to guide future placement of infusion catheters in recurrent brain tumors following primary tumor resection. The brain was modeled as a porous medium, and material property fields that account for magnetic resonance imaging (MRI)-derived anatomical regions were interpolated and directly assigned to an unstructured finite element mesh. By eliminating the need to mesh complex surfaces between fluid regions and tissue, mesh preparation was expedited, increasing the model's clinical feasibility. The infusion model predicted preferential fluid diversion into open fluid regions such as the ventricles and subarachnoid space (SAS). Additionally, a sensitivity analysis of AAV delivery demonstrated that improved AAV distribution in the tumor was achieved at higher tumor hydraulic conductivity or lower tumor porosity. Depending on the tumor infusion site, the AAV distribution covered 3.67-70.25% of the tumor volume (using a 10% AAV concentration threshold), demonstrating the model's potential to inform the selection of infusion sites for maximal tumor coverage.

The Impact of Drop Test Conditions on Brain Strain Location and Severity: A Novel Approach Using a Deep Learning Model

In contact sports such as rugby, players are at risk of sustaining traumatic brain injuries (TBI) due to high-intensity head impacts that generate high linear and rotational accelerations of the head. Previous studies have established a clear link between high-intensity head impacts and brain strains that result in concussions. This study presents a novel approach to investigating the effect of a range of laboratory controlled drop test parameters on regional peak and mean maximum principal strain (MPS) predictions within the brain using a trained convolutional neural network (CNN). The CNN is publicly available at https://github.com/Jilab-biomechanics/CNN-brain-strains . The results of this study corroborate previous findings that impacts to the side of the head result in significantly higher regional MPS than forehead impacts. Forehead impacts tend to result in the lowest region-averaged MPS values for impacts where the surface angle was at 0° and 45°, while side impacts tend to result in higher regional peak and mean MPS. The absence of a neck in drop tests resulted in lower regional peak and mean MPS values. The results indicated that the relationship between drop test parameters and resulting regional peak and mean MPS predictions is complex. The study's findings offer valuable insights into how deep learning models can be used to provide more detailed insights into how drop test conditions impact regional MPS. The novel approach used in this paper to predict brain strains can be applied in the development of better methods to reduce the brain strain resulting from head accelerations such as protective sports headgear.

Subject-specific finite element head models for skull fracture evaluation-a new tool in forensic pathology

Post-mortem computed tomography (PMCT) enables the creation of subject-specific 3D head models suitable for quantitative analysis such as finite element analysis (FEA). FEA of proposed traumatic events is an objective and repeatable numerical method for assessing whether an event could cause a skull fracture such as seen at autopsy. FEA of blunt force skull fracture in adults with subject-specific 3D models in forensic pathology remains uninvestigated. This study aimed to assess the feasibility of FEA for skull fracture analysis in routine forensic pathology. Five cases with blunt force skull fracture and sufficient information on the kinematics of the traumatic event to enable numerical reconstruction were chosen. Subject-specific finite element (FE) head models were constructed by mesh morphing based on PMCT 3D models and A Detailed and Personalizable Head Model with Axons for Injury Prediction (ADAPT) FE model. Morphing was successful in maintaining subject-specific 3D geometry and quality of the FE mesh in all cases. In three cases, the simulated fracture patterns were comparable in location and pattern to the fractures seen at autopsy/PMCT. In one case, the simulated fracture was in the parietal bone whereas the fracture seen at autopsy/PMCT was in the occipital bone. In another case, the simulated fracture was a spider-web fracture in the frontal bone, whereas a much smaller fracture was seen at autopsy/PMCT; however, the fracture in the early time steps of the simulation was comparable to autopsy/PMCT. FEA might be feasible in forensic pathology in cases with a single blunt force impact and well-described event circumstances.

Traumatic axonal injury: Clinic, forensic and biomechanics perspectives

Identification of Traumatic axonal injury (TAI) is critical in clinical practice, particularly in terms of long-term prognosis, but also for medico-legal issues, to verify whether the death or the after-effects were attributable to trauma. Multidisciplinary approaches are an undeniable asset when it comes to solving these problems. The aim of this work is therefore to list the different techniques needed to identify axonal lesions and to understand the lesion mechanisms involved in their formation. Imaging can be used to assess the consequences of trauma, to identify indirect signs of TAI, to explain the patient's initial symptoms and even to assess the patient's prognosis. Three-dimensional reconstructions of the skull can highlight fractures suggestive of trauma. Microscopic and immunohistochemical techniques are currently considered as the most reliable tools for the early identification of TAI following trauma. Finite element models use mechanical equations to predict biomechanical parameters, such as tissue stresses and strains in the brain, when subjected to external forces, such as violent impacts to the head. These parameters, which are difficult to measure experimentally, are then used to predict the risk of injury. The integration of imaging data with finite element models allows researchers to create realistic and personalized computational models by incorporating actual geometry and properties obtained from imaging techniques. The personalization of these models makes their forensic approach particularly interesting.

Mechanical and Functional Responses in Astrocytes under Alternating Deformation Modes Using Magneto-Active Substrates

This work introduces NeoMag, a system designed to enhance cell mechanics assays in substrate deformation studies. NeoMag uses multidomain magneto-active materials to mechanically actuate the substrate, transmitting reversible mechanical cues to cells. The system boasts full flexibility in alternating loading substrate deformation modes, seamlessly adapting to both upright and inverted microscopes. The multidomain substrates facilitate mechanobiology assays on 2D and 3D cultures. The integration of the system with nanoindenters allows for precise evaluation of cellular mechanical properties under varying substrate deformation modes. The system is used to study the impact of substrate deformation on astrocytes, simulating mechanical conditions akin to traumatic brain injury and ischemic stroke. The results reveal local heterogeneous changes in astrocyte stiffness, influenced by the orientation of subcellular regions relative to substrate strain. These stiffness variations, exceeding 50% in stiffening and softening, and local deformations significantly alter calcium dynamics. Furthermore, sustained deformations induce actin network reorganization and activate Piezo1 channels, leading to an initial increase followed by a long-term inhibition of calcium events. Conversely, fast and dynamic deformations transiently activate Piezo1 channels and disrupt the actin network, causing long-term cell softening. These findings unveil mechanical and functional alterations in astrocytes during substrate deformation, illustrating the multiple opportunities this technology offers.

Brain Deformation Estimation With Transfer Learning for Head Impact Datasets Across Impact Types

OBJECTIVE

The machine-learning head model (MLHM) to accelerate the calculation of brain strain and strain rate, which are the predictors for traumatic brain injury (TBI), but the model accuracy was found to decrease sharply when the training/test datasets were from different head impacts types (i.e., car crash, college football), which limits the applicability of MLHMs to different types of head impacts and sports. Particularly, small sizes of target dataset for specific impact types with tens of impacts may not be enough to train an accurate impact-type-specific MLHM.

METHODS

To overcome this, we propose data fusion and transfer learning to develop a series of MLHMs to predict the maximum principal strain (MPS) and maximum principal strain rate (MPSR).

RESULTS

The strategies were tested on American football (338), mixed martial arts (457), reconstructed car crash (48) and reconstructed American football (36) and we found that the MLHMs developed with transfer learning are significantly more accurate in estimating MPS and MPSR than other models, with a mean absolute error (MAE) smaller than 0.03 in predicting MPS and smaller than [Formula: see text] in predicting MPSR on all target impact datasets. High performance in concussion detection was observed based on the MPS and MPSR estimated by the transfer-learning-based models.

CONCLUSION

The MLHMs can be applied to various head impact types for rapidly and accurately calculating brain strain and strain rate.

SIGNIFICANCE

This study enables developing MLHMs for the head impact type with limited availability of data, and will accelerate the applications of MLHMs.

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.

Revealing the role of material properties in impact-related injuries: Investigating the influence of brain and skull density variations on head injury severity

Traumatic brain injuries (TBI) resulting from head impacts are a major public health concern, which prompted our research to investigate the complex relationship between the material properties of brain tissue and the severity of TBI. The goal of this research is to investigate how variations in brain and skull density influence the vulnerability of brain tissue to traumatic injury, thereby enhancing our understanding of injury mechanism. To achieve this goal, we employed a well-validated finite element head model (FEHM). The current investigation was divided into two phases: in the first one, three distinct brain viscoelastic materials that had been utilized in prior studies were analyzed. The review of the properties of these three materials has been meticulous, encompassing both the spectrum of mechanical properties and the behaviors that are relevant to the way in which brain tissue reacts to traumatic loading conditions. In the second phase, the material properties of both the brain and skull tissue, alongside the impact conditions, were held constant. After this step, the focus was directed towards the variation of density in the brain and skull, which was consistent with the results obtained from previous experimental investigations, in order to determine the precise impact of these variations in density. This approach allowed a more profound comprehension of the impacts that density had on the simulation results. In the first phase, Material No. 2 exhibited the highest maximum first principal strain value in the frontal region (ε max = 15.41 % ), indicating lower stiffness to instantaneous deformation. This characteristic suggests that Material No. 2 may deform more extensively upon impact, potentially increasing the risk of injury due to its viscoelastic behavior. In contrast, Material No. 1, with a lower maximum first principal strain in the frontal region (ε max = 7.87 % ), displayed greater stiffness to instantaneous deformation, potentially reducing the risk of brain injury upon head impact. The second phase provided quantitative findings revealing a proportional relationship between brain tissue density and the pressures experienced by the brain. A 2 % increase in brain tissue density corresponded to approximately a 1 % increase in pressure on the brain tissue. Similarly, changes in skull density exhibited a similar quantitative relationship, with a 6 % increase in skull density leading to a 2.5 % increase in brain pressure. This preliminary approximate ratio of 2 to 1 between brain and skull density variations provides an initial quantitative framework for assessing the impact of density changes on brain vulnerability. These findings have several implications for the development of protective measures and injury prevention strategies, particularly in contexts where head trauma is a major issue.

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.

Head injury: Importance of the deep brain nuclei in force transmission to the brain

Finite element modeling provides a digital representation of the human body. It is currently the most pertinent method to study the mechanisms of head injury, and is becoming a scientific reference in forensic expert reports. Improved biofidelity is a recurrent aim of research studies in biomechanics in order to improve earlier models whose mechanical properties conformed to simplified elastic behavior and mechanic laws. We aimed to study force transmission to the brain following impacts to the head, using a finite element head model with increased biofidelity. To the model developed by the Laboratory of Applied Biomechanics of Marseille, we added new brain structures (thalamus, central gray nuclei and ventricular systems) as well as three tracts involved in the symptoms of head injury: the corpus callosum, uncinate tracts and corticospinal tracts. Three head impact scenarios were simulated: an uppercut with the prior model and an uppercut with the improved model in order to compare the two models, and a lateral impact with an impact velocity of 6.5 m/s in the improved model. In these conditions, in uppercuts the maximum stress values did not exceed the injury risk threshold. On the other hand, the deep gray matter (thalamus and central gray nuclei) was the region at highest risk of injury during lateral impacts. Even if injury to the deep gray matter is not immediately life-threatening, it could explain the chronic disabling symptoms of even low-intensity head injury.

Prediction of skull fractures in blunt force head traumas using finite element head models

Traumatic head injuries remain a leading cause of death and disability worldwide. Although skull fractures are one of the most common head injuries, the fundamental mechanics of cranial bone and its impact tolerance are still uncertain. In the present study, a strain-rate-dependent material model for cranial bone has been proposed and implemented in subject-specific Finite Element (FE) head models in order to predict skull fractures in five real-world fall accidents. The subject-specific head models were developed following an established image-registration-based personalization pipeline. Head impact boundary conditions were derived from accident reconstructions using personalized human body models. The simulated fracture lines were compared to those visible in post-mortem CT scans of each subject. In result, the FE models did predict the actual occurrence and extent of skull fractures in all cases. In at least four out of five cases, predicted fracture patterns were comparable to ones from CT scans and autopsy reports. The tensile material model, which was tuned to represent rate-dependent tensile data of cortical skull bone from literature, was able to capture observed linear fractures in blunt indentation loading of a skullcap specimen. The FE model showed to be sensitive to modeling parameters, in particular to the constitutive parameters of the cortical tables. Nevertheless, this study provides a currently lacking strain-rate dependent material model of cranial bone that has the capacity to accurately predict linear fracture patterns. For the first time, a procedure to reconstruct occurrences of skull fractures using computational engineering techniques, capturing the all-in-all fracture initiation, propagation and final pattern, is presented.

Evaluating child helmet protection and testing standards: A study using PIPER child head models aged 1.5, 3, 6, and 18 years

The anatomy of children's heads is unique and distinct from adults, with smaller and softer skulls and unfused fontanels and sutures. Despite this, most current helmet testing standards for children use the same peak linear acceleration threshold as for adults. It is unclear whether this is reasonable and otherwise what thresholds should be. To answer these questions, helmet-protected head responses for different ages are needed which is however lacking today. In this study, we apply continuously scalable PIPER child head models of 1.5, 3, and 6 years old (YO), and an upgraded 18YO to study child helmet protection under extensive linear and oblique impacts. The results of this study reveal an age-dependence trend in both global kinematics and tissue response, with younger children experiencing higher levels of acceleration and velocity, as well as increased skull stress and brain strain. These findings indicate the need for better protection for younger children, suggesting that youth helmets should have a lower linear kinematic threshold, with a preliminary value of 150g for 1.5-year-old helmets. However, the results also show a different trend in rotational kinematics, indicating that the threshold of rotational velocity for a 1.5YO is similar to that for adults. The results also support the current use of small-sized adult headforms for testing child helmets before new child headforms are available.

Development of personalizable female and male pedestrian SAFER human body models

OBJECTIVES

Vulnerable road users are globally overrepresented as victims of road traffic injuries. Developing biofidelic male and female pedestrian human body models (HBMs) that represent diverse anthropometries is essential to enhance road safety and propose intervention strategies.

METHODS

In this study, 50th percentile male and female pedestrians of the SAFER HBM were developed via a newly developed image registration-based mesh morphing framework. The performance of the HBMs was evaluated by means of a set of cadaver experiments, involving subjects struck laterally by a generic sedan buck.

RESULTS

In simulated whole-body pedestrian collisions, the personalized HBMs effectively replicate trajectories of the head and lower body regions, as well as head kinematics, in lateral impacts. The results also demonstrate the personalization framework's capacity to generate personalized HBMs with reliable mesh quality, ensuring robust simulations.

CONCLUSIONS

The presented pedestrian HBMs and personalization framework provide robust means to reconstruct and evaluate head impacts in pedestrian-to-vehicle collisions thoroughly and accurately.

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.

An atlas of the heterogeneous viscoelastic brain with local power-law attenuation synthesised using Prony-series

This review addresses the acute need to acknowledge the mechanical heterogeneity of brain matter and to accurately calibrate its local viscoelastic material properties accordingly. Specifically, it is important to compile the existing and disparate literature on attenuation power-laws and dispersion to make progress in wave physics of brain matter, a field of research that has the potential to explain the mechanisms at play in diffuse axonal injury and mild traumatic brain injury in general. Currently, viscous effects in the brain are modelled using Prony-series, i.e., a sum of decaying exponentials at different relaxation times. Here we collect and synthesise the Prony-series coefficients appearing in the literature for twelve regions: brainstem, basal ganglia, cerebellum, corona radiata, corpus callosum, cortex, dentate gyrus, hippocampus, thalamus, grey matter, white matter, homogeneous brain, and for eight different mammals: pig, rat, human, mouse, cow, sheep, monkey and dog. Using this data, we compute the fractional-exponent attenuation power-laws for different tissues of the brain, the corresponding dispersion laws resulting from causality, and the averaged Prony-series coefficients. STATEMENT OF SIGNIFICANCE: Traumatic brain injuries are considered a silent epidemic and finite element methods (FEMs) are used in modelling brain deformation, requiring access to viscoelastic properties of brain. To the best of our knowledge, this work presents 1) the first multi-frequency viscoelastic atlas of the heterogeneous brain, 2) the first review focusing on viscoelastic modelling in both FEMs and experimental works, 3) the first attempt to conglomerate the disparate existing literature on the viscoelastic modelling of the brain and 4) the largest collection of viscoelastic parameters for the brain (212 different Prony-series spanning 12 different tissues and 8 different animal surrogates). Furthermore, this work presents the first brain atlas of attenuation power-laws essential for modelling shear waves in brain.

Machine-learning-based head impact subtyping based on the spectral densities of the measurable head kinematics

BACKGROUND

Traumatic brain injury can be caused by head impacts, but many brain injury risk estimation models are not equally accurate across the variety of impacts that patients may undergo, and the characteristics of different types of impacts are not well studied. We investigated the spectral characteristics of different head impact types with kinematics classification.

METHODS

Data were analyzed from 3262 head impacts from lab reconstruction, American football, mixed martial arts, and publicly available car crash data. A random forest classifier with spectral densities of linear acceleration and angular velocity was built to classify head impact types (e.g., football, car crash, mixed martial arts). To test the classifier robustness, another 271 lab-reconstructed impacts were obtained from 5 other instrumented mouthguards. Finally, with the classifier, type-specific, nearest-neighbor regression models were built for brain strain.

RESULTS

The classifier reached a median accuracy of 96% over 1000 random partitions of training and test sets. The most important features in the classification included both low- and high-frequency features, both linear acceleration features and angular velocity features. Different head impact types had different distributions of spectral densities in low- and high-frequency ranges (e.g., the spectral densities of mixed martial arts impacts were higher in the high-frequency range than in the low-frequency range). The type-specific regression showed a generally higher R2 value than baseline models without classification.

CONCLUSION

The machine-learning-based classifier enables a better understanding of the impact kinematics spectral density in different sports, and it can be applied to evaluate the quality of impact-simulation systems and on-field data augmentation.

Investigating the Effect of Brain Size on Deformation Magnitude Using Subject-Specific Finite Element Models

Abstract In the last decade, computational models of the brain have become the gold standard tool for investigating traumatic brain injury (TBI) mechanisms and developing novel protective equipment and other safety countermeasures. However, most studies utilizing finite element (FE) models of the brain have been conducted using models developed to represent the average neuroanatomy of a target demographic, such as the 50th percentile male. Although this is an efficient strategy, it neglects normal anatomical variations present within the population and their contributions on the brain's deformation response. As a result, the contributions of structural characteristics of the brain, such as brain volume, on brain deformation are not well understood. The objective of this study was to develop a set of statistical regression models relating measures of the size and shape of the brain to the resulting brain deformation. This was performed using a database of 125 subject-specific models, simulated under six independent head kinematic boundary conditions, spanning a range of impact modes (frontal, oblique, side), severity (non-injurious and injurious), and environments (volunteer, automotive, and American football). Two statistical regression techniques were utilized. First, simple linear regression (SLR) models were trained to relate intracranial volume (ICV) and the 95th percentile of maximum principal strain (MPS-95) for each of the impact cases. Second, a partial least squares regression model was constructed to predict MPS-95 based on the affine transformation parameters from each subject, representing the size and shape of their brain, considering the six impact conditions collectively. Both techniques indicated a strong linear relationship between ICV and MPS-95, with MPS-95 varying by approximately 5% between the smallest and largest brains. This difference represented up to 40% of the mean strain across all subjects. This study represents a comprehensive assessment of the relationships between brain anatomy and deformation, which is crucial for the development of personalized protective equipment, identifying individuals at higher risk of injury, and using computational models to aid clinical diagnostics of TBI.

Personalization of human body models and beyond via image registration

Finite element human body models (HBMs) are becoming increasingly important numerical tools for traffic safety. Developing a validated and reliable HBM from the start requires integrated efforts and continues to be a challenging task. Mesh morphing is an efficient technique to generate personalized HBMs accounting for individual anatomy once a baseline model has been developed. This study presents a new image registration-based mesh morphing method to generate personalized HBMs. The method is demonstrated by morphing four baseline HBMs (SAFER, THUMS, and VIVA+ in both seated and standing postures) into ten subjects with varying heights, body mass indices (BMIs), and sex. The resulting personalized HBMs show comparable element quality to the baseline models. This method enables the comparison of HBMs by morphing them into the same subject, eliminating geometric differences. The method also shows superior geometry correction capabilities, which facilitates converting a seated HBM to a standing one, combined with additional positioning tools. Furthermore, this method can be extended to personalize other models, and the feasibility of morphing vehicle models has been illustrated. In conclusion, this new image registration-based mesh morphing method allows rapid and robust personalization of HBMs, facilitating personalized simulations.

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.

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.

Image-based axon model highlights heterogeneity in initiation of damage

Head injury simulations predict the occurrence of traumatic brain injury by placing a threshold on the calculated strains for axon tracts within the brain. However, a current roadblock to accurate injury prediction is the selection of an appropriate axon damage threshold. While several computational studies have used models of the axon cytoskeleton to investigate damage initiation, these models all employ an idealized, homogeneous axonal geometry. This homogeneous geometry with regularly spaced microtubules, evenly distributed throughout the model, overestimates axon strength because, in reality, the axon cytoskeleton is heterogeneous. In the heterogeneous cytoskeleton, the weakest cross section determines the initiation of failure, but these weak spots are not present in a homogeneous model. Addressing one source of heterogeneity in the axon cytoskeleton, we present a new semiautomated image analysis pipeline for using serial-section transmission electron micrographs to reconstruct the microtubule geometry of an axon. The image analysis procedure locates microtubules within the images, traces them throughout the image stack, and reconstructs the microtubule structure as a finite element mesh. We demonstrate the image analysis approach using a C. elegans touch receptor neuron due to the availability of high-quality serial-section transmission electron micrograph data sets. The results of the analysis highlight the heterogeneity of the microtubule structure in the spatial variation of both microtubule number and length. Simulations comparing this image-based geometry with homogeneous geometries show that structural heterogeneity in the image-based model creates significant spatial variation in deformation. The homogeneous geometries, on the other hand, deform more uniformly. Since no single homogeneous model can replicate the mechanical behavior of the image-based model, our results argue that heterogeneity in axon microtubule geometry should be considered in determining accurate axon failure thresholds.

[Mathematical modeling using finite element analysis in forensic medical examination]

The purpose of the work is to develop methods of mathematical modeling using finite element analysis in forensic medical examination. The stages of the methodology for solving problems of deformable body mechanics in forensic medicine are considered, which allows to reliably establish the possibility of formation and morphology of damage under specific conditions and circumstances, to focus the researcher's attention on problem points when creating and evaluating the model. The use of simplified models of the human body makes the expert's conclusion more reasonable, which increases the confidence of law enforcement agencies in the activities of the forensic medical expert service and allows for a new look at solving the problems of forensic medicine and forensic medical examination.

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.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser).

SIGNIFICANCE

There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

Finding the Spatial Co-Variation of Brain Deformation With Principal Component Analysis

OBJECTIVE

Strain and strain rate are effective traumatic brain injury metrics. In finite element (FE) head model, thousands of elements were used to represent the spatial distribution of these metrics. Owing that these metrics are resulted from brain inertia, their spatial distribution can be represented in more concise pattern. Since head kinematic features and brain deformation vary largely across head impact types (Zhan et al., 2021), we applied principal component analysis (PCA) to find the spatial co-variation of injury metrics (maximum principal strain (MPS), MPS rate (MPSR) and MPS × MPSR) in four impact types: simulation, football, mixed martial arts and car crashes, and used the PCA to find patterns in these metrics and improve the machine learning head model (MLHM).

METHODS

We applied PCA to decompose the injury metrics for all impacts in each impact type, and investigate the spatial co-variation using the first principal component (PC1). Furthermore, we developed a MLHM to predict PC1 and then inverse-transform to predict for all brain elements. The accuracy, the model complexity and the size of training dataset of PCA-MLHM are compared with previous MLHM (Zhan et al., 2021).

RESULTS

PC1 explained variance on the datasets. Based on PC1 coefficients, the corpus callosum and midbrain exhibit high variance on all datasets. Finally, the PCA-MLHM reduced model parameters by 74% with a similar MPS estimation accuracy.

CONCLUSION

The brain injury metric in a dataset can be decomposed into mean components and PC1 with high explained variance.

SIGNIFICANCE

The spatial co-variation analysis enables better interpretation of the patterns in brain injury metrics. It also improves the efficiency of MLHM.

Quantifying the Effect of Sex and Neuroanatomical Biomechanical Features on Brain Deformation Response in Finite Element Brain Models

Recent automotive epidemiology studies have concluded that females have significantly higher odds of sustaining a moderate brain injury or concussion than males in a frontal crash after controlling for multiple crash and occupant variables. Differences in neuroanatomical features, such as intracranial volume (ICV), have been shown between male and female subjects, but how these sex-specific neuroanatomical differences affect brain deformation is unknown. This study used subject-specific finite element brain models, generated via registration-based morphing using both male and female magnetic resonance imaging scans, to investigate sex differences of a variety of neuroanatomical features and their effect on brain deformation; additionally, this study aimed to determine the relative importance of these neuroanatomical features and sex on brain deformation metrics for a single automotive loading environment. Based on the Bayesian linear mixed models, sex had a significant effect on ICV, white matter volume and gray matter volume, as well as a section of cortical gray matter regions' thicknesses and volumes; however, after these neuroanatomical features were accounted for in the statistical model, sex was not a significant factor in predicting brain deformation. ICV had the highest relative effect on the brain deformation metrics assessed. Therefore, ICV should be considered when investigating both brain injury biomechanics and injury risk.

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.

American Football Helmet Effectiveness Against a Strain-Based Concussion Mechanism

Brain strain is increasingly being used in helmet design and safety performance evaluation as it is generally considered as the primary mechanism of concussion. In this study, we investigate whether different helmet designs can meaningfully alter brain strains using two commonly used metrics, peak maximum principal strain (MPS) of the whole brain and cumulative strain damage measure (CSDM). A convolutional neural network (CNN) that instantly produces detailed brain strains is first tested for accuracy for helmeted head impacts. Based on N = 144 impacts in 12 impact conditions from three random and representative helmet models, we conclude that the CNN is sufficiently accurate for helmet testing applications, for elementwise MPS (success rate of 98.6%), whole-brain peak MPS and CSDM (coefficient of determination of 0.977 and 0.980, with root mean squared error of 0.015 and 0.029, respectively). We then apply the technique to 23 football helmet models (N = 1104 impacts) to reproduce elementwise MPS. Assuming a concussion would occur when peak MPS or CSDM exceeds a threshold, we sweep their thresholds across the value ranges to evaluate the number of predicted hypothetical concussions that different helmets sustain across the impact conditions. Relative to the 12 impact conditions tested, we find that the "best" and "worst" helmets differ by an average of 22.5% in terms of predicted concussions, ranging from 0 to 42% (the latter achieved at the threshold value of 0.28 for peak MPS and 0.4 for CSDM, respectively). Such a large variation among helmets in strain-based concussion predictions demonstrate that helmet designs can still be optimized in a clinically meaningful way. The robustness and accuracy of the CNN tool also suggest its potential for routine use for helmet design and safety performance evaluation in the future. The CNN is freely available online at https://github.com/Jilab-biomechanics/CNN-brain-strains .

Real-time dynamic simulation for highly accurate spatiotemporal brain deformation from impact

Real-time dynamic simulation remains a significant challenge for spatiotemporal data of high dimension and resolution. In this study, we establish a transformer neural network (TNN) originally developed for natural language processing and a separate convolutional neural network (CNN) to estimate five-dimensional (5D) spatiotemporal brain-skull relative displacement resulting from impact (isotropic spatial resolution of 4 mm with temporal resolution of 1 ms). Sequential training is applied to train (N = 5184 samples) the two neural networks for estimating the complete 5D displacement across a temporal duration of 60 ms. We find that TNN slightly but consistently outperforms CNN in accuracy for both displacement and the resulting voxel-wise four-dimensional (4D) maximum principal strain (e.g., root mean squared error (RMSE) of ~1.0% vs. ~1.6%, with coefficient of determination, R ² >0.99 vs. >0.98, respectively, and normalized RMSE (NRMSE) at peak displacement of 2%-3%, based on an independent testing dataset; N = 314). Their accuracies are similar for a range of real-world impacts drawn from various published sources (dummy, helmet, football, soccer, and car crash; average RMSE/NRMSE of ~0.3 mm/~4%-5% and average R ² of ~0.98 at peak displacement). Sequential training is effective for allowing instantaneous estimation of 5D displacement with high accuracy, although TNN poses a heavier computational burden in training. This work enables efficient characterization of the intrinsically dynamic brain strain in impact critical for downstream multiscale axonal injury model simulation. This is also the first application of TNN in biomechanics, which offers important insight into how real-time dynamic simulations can be achieved across diverse engineering fields.

The Presence of the Temporal Horn Exacerbates the Vulnerability of Hippocampus During Head Impacts

Hippocampal injury is common in traumatic brain injury (TBI) patients, but the underlying pathogenesis remains elusive. In this study, we hypothesize that the presence of the adjacent fluid-containing temporal horn exacerbates the biomechanical vulnerability of the hippocampus. Two finite element models of the human head were used to investigate this hypothesis, one with and one without the temporal horn, and both including a detailed hippocampal subfield delineation. A fluid-structure interaction coupling approach was used to simulate the brain-ventricle interface, in which the intraventricular cerebrospinal fluid was represented by an arbitrary Lagrangian-Eulerian multi-material formation to account for its fluid behavior. By comparing the response of these two models under identical loadings, the model that included the temporal horn predicted increased magnitudes of strain and strain rate in the hippocampus with respect to its counterpart without the temporal horn. This specifically affected cornu ammonis (CA) 1 (CA1), CA2/3, hippocampal tail, subiculum, and the adjacent amygdala and ventral diencephalon. These computational results suggest that the presence of the temporal horn exacerbate the vulnerability of the hippocampus, highlighting the mechanobiological dependency of the hippocampus on the temporal horn.

Modelling glioma progression, mass effect and intracranial pressure in patient anatomy

Increased intracranial pressure is the source of most critical symptoms in patients with glioma, and often the main cause of death. Clinical interventions could benefit from non-invasive estimates of the pressure distribution in the patient's parenchyma provided by computational models. However, existing glioma models do not simulate the pressure distribution and they rely on a large number of model parameters, which complicates their calibration from available patient data. Here we present a novel model for glioma growth, pressure distribution and corresponding brain deformation. The distinct feature of our approach is that the pressure is directly derived from tumour dynamics and patient-specific anatomy, providing non-invasive insights into the patient's state. The model predictions allow estimation of critical conditions such as intracranial hypertension, brain midline shift or neurological and cognitive impairments. A diffuse-domain formalism is employed to allow for efficient numerical implementation of the model in the patient-specific brain anatomy. The model is tested on synthetic and clinical cases. To facilitate clinical deployment, a high-performance computing implementation of the model has been publicly released.

Subject-Specific Head Model Generation by Mesh Morphing: A Personalization Framework and Its Applications

Finite element (FE) head models have become powerful tools in many fields within neuroscience, especially for studying the biomechanics of traumatic brain injury (TBI). Subject-specific head models accounting for geometric variations among subjects are needed for more reliable predictions. However, the generation of such models suitable for studying TBIs remains a significant challenge and has been a bottleneck hindering personalized simulations. This study presents a personalization framework for generating subject-specific models across the lifespan and for pathological brains with significant anatomical changes by morphing a baseline model. The framework consists of hierarchical multiple feature and multimodality imaging registrations, mesh morphing, and mesh grouping, which is shown to be efficient with a heterogeneous dataset including a newborn, 1-year-old (1Y), 2Y, adult, 92Y, and a hydrocephalus brain. The generated models of the six subjects show competitive personalization accuracy, demonstrating the capacity of the framework for generating subject-specific models with significant anatomical differences. The family of the generated head models allows studying age-dependent and groupwise brain injury mechanisms. The framework for efficient generation of subject-specific FE head models helps to facilitate personalized simulations in many fields of neuroscience.

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.

Predictive Factors of Kinematics in Traumatic Brain Injury from Head Impacts Based on Statistical Interpretation

Brain tissue deformation resulting from head impacts is primarily caused by rotation and can lead to traumatic brain injury. To quantify brain injury risk based on measurements of kinematics on the head, finite element (FE) models and various brain injury criteria based on different factors of these kinematics have been developed, but the contribution of different kinematic factors has not been comprehensively analyzed across different types of head impacts in a data-driven manner. To better design brain injury criteria, the predictive power of rotational kinematics factors, which are different in (1) the derivative order (angular velocity, angular acceleration, angular jerk), (2) the direction and (3) the power (e.g., square-rooted, squared, cubic) of the angular velocity, were analyzed based on different datasets including laboratory impacts, American football, mixed martial arts (MMA), NHTSA automobile crashworthiness tests and NASCAR crash events. Ordinary least squares regressions were built from kinematics factors to the 95% maximum principal strain (MPS95), and we compared zero-order correlation coefficients, structure coefficients, commonality analysis, and dominance analysis. The angular acceleration, the magnitude and the first power factors showed the highest predictive power for the majority of impacts including laboratory impacts, American football impacts, with few exceptions (angular velocity for MMA and NASCAR impacts). The predictive power of rotational kinematics about three directions (x: posterior-to-anterior, y: left-to-right, z: superior-to-inferior) of kinematics varied with different sports and types of head impacts.

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.

Time Window of Head Impact Kinematics Measurement for Calculation of Brain Strain and Strain Rate in American Football

Wearable devices have been shown to effectively measure the head's movement during impacts in sports like American football. When a head impact occurs, the device is triggered to collect and save the kinematic measurements during a predefined time window. Then, based on the collected kinematics, finite element (FE) head models can calculate brain strain and strain rate, which are used to evaluate the risk of mild traumatic brain injury. To find a time window that can provide a sufficient duration of kinematics for FE analysis, we investigated 118 on-field video-confirmed football head impacts collected by the Stanford Instrumented Mouthguard. The simulation results based on the kinematics truncated to a shorter time window were compared with the original to determine the minimum time window needed for football. Because the individual differences in brain geometry influence these calculations, we included six representative brain geometries and found that larger brains need a longer time window of kinematics for accurate calculation. Among the different sizes of brains, a pre-trigger time of 40 ms and a post-trigger time of 70 ms were found to yield calculations of brain strain and strain rate that were not significantly different from calculations using the original 200 ms time window recorded by the mouthguard. Therefore, approximately 110 ms is recommended for complete modeling of impacts for football.

Influence of Anisotropic White Matter on Electroosmotic Flow Induced by Direct Current

Treatment of cerebral edema remains a major challenge in clinical practice and new innovative therapies are needed. This study presents a novel approach for mitigating cerebral edema by inducing bulk fluid transport utilizing the brain’s electroosmotic property using an anatomically detailed finite element head model incorporating anisotropy in the white matter (WM). Three representative anisotropic conductivity algorithms are employed for the WM and compared with isotropic WM. The key results are (1) the electroosmotic flow (EOF) is driven from the edema region to the subarachnoid space under an applied electric field with its magnitude linearly correlated to the electric field and direction following current flow pathways; (2) the extent of EOF distribution variation correlates highly with the degree of the anisotropic ratio of the WM regions; (3) the directions of the induced EOF in the anisotropic models deviate from its isotropically defined pathways and tend to move along the principal fiber direction. The results suggest WM anisotropy should be incorporated in head models for more reliable EOF evaluations for cerebral edema mitigation and demonstrate the promise of the electroosmosis based approach to be developed as a new therapy for edema treatment as evaluated with enhanced head models incorporating WM anisotropy.

Effect of Strain Rate on Single Tau, Dimerized Tau and Tau-Microtubule Interface: A Molecular Dynamics Simulation Study

Microtubule-associated protein (MAP) tau is a cross-linking molecule that provides structural stability to axonal microtubules (MT). It is considered a potential biomarker for Alzheimer's disease (AD), dementia, and other neurological disorders. It is also a signature protein for Traumatic Brain Injury (TBI) assessment. In the case of TBI, extreme dynamic mechanical energies can be felt by the axonal cytoskeletal members. As such, fundamental understandings of the responses of single tau protein, polymerized tau protein, and tau-microtubule interfaces under high-rate mechanical forces are important. This study attempts to determine the high-strain rate mechanical behavior of single tau, dimerized tau, and tau-MT interface using molecular dynamics (MD) simulation. The results show that a single tau protein is a highly stretchable soft polymer. During deformation, first, it significantly unfolds against van der Waals and electrostatic bonds. Then it stretches against strong covalent bonds. We found that tau acts as a viscoelastic material, and its stiffness increases with the strain rate. The unfolding stiffness can be ~50-500 MPa, while pure stretching stiffness can be >2 GPa. The dimerized tau model exhibits similar behavior under similar strain rates, and tau sliding from another tau is not observed until it is stretched to >7 times of original length, depending on the strain rate. The tau-MT interface simulations show that very high strain and strain rates are required to separate tau from MT suggesting Tau-MT bonding is stronger than MT subunit bonding between themselves. The dimerized tau-MT interface simulations suggest that tau-tau bonding is stronger than tau-MT bonding. In summary, this study focuses on the structural response of individual cytoskeletal components, namely microtubule (MT) and tau protein. Furthermore, we consider not only the individual response of a component, but also their interaction with each other (such as tau with tau or tau with MT). This study will eventually pave the way to build a bottom-up multiscale brain model and analyze TBI more comprehensively.

Measurement of brain simulant strains in head surrogate under impact loading

Impact-induced traumatic brain injury (TBI) is a major source of disability and mortality. Knowledge of brain strains during impact (accelerative) loading is critical for the overall management of TBI, including the development of injury thresholds, personal protective equipment, and validation of computational models. Despite these needs, the current understanding of brain strains in humans or humanlike surrogates is limited, especially for injury causing loading magnitudes. Toward this end, we measured full-field, in-plane (2D) strains in a brain simulant using the hemispherical head surrogate. The hemispherical head was mounted on the Hybrid-III neck and subjected to impact loading using a linear impactor system. The resulting head kinematics was measured using a triaxial accelerometer and angular rate sensors. Dynamic, 2D strains in a brain simulant were obtained using high-speed imaging and digital image correlation. Concurrent finite element (FE) simulations of the experiment were also performed to gain additional insights. The role of stiff membranes of the head was also studied using experiments. Our results suggest that rotational modes dominate the response of the brain simulant. The wave propagation in the brain simulant as a result of impact has a timescale of ~100 ms. We obtain peak strains of ~20%, ~40%, ~60% for peak rotational accelerations of ~838, ~5170, ~11,860 rad/s², respectively. Further, peak strains in cortical regions are higher than subcortical regions by up to ~70%. The agreement between the experiments and FE simulations is reasonable in terms of spatiotemporal evolution of strain pattern and peak strain magnitudes. Experiments with the addition of falx and tentorium indicate significant strain concentration (up to 115%) in the brain simulant near the interface of falx or tentorium and brain simulant. Overall, this work provides important insights into the biomechanics of strain in the brain simulant during impact loading.

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.

Symmetry of the Human Head—Are Symmetrical Models More Applicable in Numerical Analysis?

The study of symmetrical and non-symmetrical effects in physics, mathematics, mechanics, medicine, and numerical methods is a current topic due to the complexity of the experiments, calculations, and virtual simulations. However, there is a limited number of research publications in computational biomechanics focusing on the symmetry of numerical head models. The majority of the models in the researched literature are symmetrical. Thus, we stated a hypothesis wherever the symmetrical models might be more applicable in numerical analysis. We carried out in-depth studies about head symmetry through clinical data, medical images, materials models, and computer analysis. We concluded that the mapping of the entire geometry of the skull and brain is essential due to the significant differences that affect the results of numerical analyses and the possibility of misinterpretation of the tissue deformation under mechanical load results.

A Methodology to Compare Biomechanical Simulations With Clinical Brain Imaging Analysis Utilizing Two Blunt Impact Cases

According to the US Defense and Veterans Brain Injury Center (DVBIC) and Centers for Disease Control and Prevention (CDC), mild traumatic brain injury (mTBI) is a common form of head injury. Medical imaging data provides clinical insight into tissue damage/injury and injury severity, and helps medical diagnosis. Computational modeling and simulation can predict the biomechanical characteristics of such injury, and are useful for development of protective equipment. Integration of techniques from computational biomechanics with medical data assessment modalities (e.g., magnetic resonance imaging or MRI) has not yet been used to predict injury, support early medical diagnosis, or assess effectiveness of personal protective equipment. This paper presents a methodology to map computational simulations with clinical data for interpreting blunt impact TBI utilizing two clinically different head injury case studies. MRI modalities, such as T1, T2, diffusion-weighted imaging (DWI) and apparent diffusion coefficient (ADC), were used for simulation comparisons. The two clinical cases have been reconstructed using finite element analysis to predict head biomechanics based on medical reports documented by a clinician. The findings are mapped to simulation results using image-based clinical analyses of head impact injuries, and modalities that could capture simulation results have been identified. In case 1, the MRI results showed lesions in the brain with skull indentation, while case 2 had lesions in both coup and contrecoup sides with no skull deformation. Simulation data analyses show that different biomechanical measures and thresholds are needed to explain different blunt impact injury modalities; specifically, strain rate threshold corresponds well with brain injury with skull indentation, while minimum pressure threshold corresponds well with coup–contrecoup injury; and DWI has been found to be the most appropriate modality for MRI data interpretation. As the findings from these two cases are substantiated with additional clinical studies, this methodology can be broadly applied as a tool to support injury assessment in head trauma events and to improve countermeasures (e.g., diagnostics and protective equipment design) to mitigate these injuries.

The relationship between brain injury criteria and brain strain across different types of head impacts can be different

Multiple brain injury criteria (BIC) are developed to quickly quantify brain injury risks after head impacts. These BIC originated from different head impact types (e.g. sports and car crashes) are widely used in risk evaluation. However, the accuracy of using the BIC on brain injury risk estimation across head impact types has not been evaluated. Physiologically, brain strain is often considered the key parameter of brain injury. To evaluate the BIC's risk estimation accuracy across five datasets comprising different head impact types, linear regression was used to model 95% maximum principal strain, 95% maximum principal strain at the corpus callosum and cumulative strain damage (15%) on 18 BIC. The results show significantly different relationships between BIC and brain strain across datasets, indicating the same BIC value may suggest different brain strain across head impact types. The accuracy of brain strain regression is generally decreasing if the BIC regression models are fitted on a dataset with a different type of head impact rather than on the dataset with the same type. Given this finding, this study raises concerns for applying BIC to estimate the brain injury risks for head impacts different from the head impacts on which the BIC was developed.

Rapid Estimation of Entire Brain Strain Using Deep Learning Models

OBJECTIVE

Many recent studies have suggested that brain deformation resulting from a head impact is linked to the corresponding clinical outcome, such as mild traumatic brain injury (mTBI). Even though several finite element (FE) head models have been developed and validated to calculate brain deformation based on impact kinematics, the clinical application of these FE head models is limited due to the time-consuming nature of FE simulations. This work aims to accelerate the process of brain deformation calculation and thus improve the potential for clinical applications.

METHODS

We propose a deep learning head model with a five-layer deep neural network and feature engineering, and trained and tested the model on 2511 total head impacts from a combination of head model simulations and on-field college football and mixed martial arts impacts.

RESULTS

The proposed deep learning head model can calculate the maximum principal strain (Green Lagrange) for every element in the entire brain in less than 0.001s with an average root mean squared error of 0.022, and with a standard deviation of 0.001 over twenty repeats with random data partition and model initialization.

CONCLUSION

Trained and tested using the dataset of 2511 head impacts, this model can be applied to various sports in the calculation of brain strain with accuracy, and its applicability can even further be extended by incorporating data from other types of head impacts.

SIGNIFICANCE

In addition to the potential clinical application in real-time brain deformation monitoring, this model will help researchers estimate the brain strain from a large number of head impacts more efficiently than using FE models.