Computational study of human head response to primary blast waves of five levels from three directions

Double Blast Wave Primary Effect on Synaptic, Glymphatic, Myelin, Neuronal and Neurovascular Markers

Explosive blasts are associated with neurological consequences as a result of blast waves impact on the brain. Yet, the neuropathologic and molecular consequences due to blast waves vs. blunt-TBI are not fully understood. An explosive-driven blast-generating system was used to reproduce blast wave exposure and examine pathological and molecular changes generated by primary wave effects of blast exposure. We assessed if pre- and post-synaptic (synaptophysin, PSD-95, spinophilin, GAP-43), neuronal (NF-L), glymphatic (LYVE1, podoplanin), myelin (MBP), neurovascular (AQP4, S100β, PDGF) and genomic (DNA polymerase-β, RNA polymerase II) markers could be altered across different brain regions of double blast vs. sham animals. Twelve male rats exposed to two consecutive blasts were compared to 12 control/sham rats. Western blot, ELISA, and immunofluorescence analyses were performed across the frontal cortex, hippocampus, cerebellum, and brainstem. The results showed altered levels of AQP4, S100β, DNA-polymerase-β, PDGF, synaptophysin and PSD-95 in double blast vs. sham animals in most of the examined regions. These data indicate that blast-generated changes are preferentially associated with neurovascular, glymphatic, and DNA repair markers, especially in the brainstem. Moreover, these changes were not accompanied by behavioral changes and corroborate the hypothesis for which an asymptomatic altered status is caused by repeated blast exposures.

Cerebral Vasculature Influences Blast-Induced Biomechanical Responses of Human Brain Tissue

Multiple finite-element (FE) models to predict the biomechanical responses in the human brain resulting from the interaction with blast waves have established the importance of including the brain-surface convolutions, the major cerebral veins, and using non-linear brain-tissue properties to improve model accuracy. We hypothesize that inclusion of a more detailed network of cerebral veins and arteries can further enhance the model-predicted biomechanical responses and help identify correlates of blast-induced brain injury. To more comprehensively capture the biomechanical responses of human brain tissues to blast-wave exposure, we coupled a three-dimensional (3-D) detailed-vasculature human-head FE model, previously validated for blunt impact, with a 3-D shock-tube FE model. Using the coupled model, we computed the biomechanical responses of a human head facing an incoming blast wave for blast overpressures (BOPs) equivalent to 68, 83, and 104 kPa. We validated our FE model, which includes the detailed network of cerebral veins and arteries, the gyri and the sulci, and hyper-viscoelastic brain-tissue properties, by comparing the model-predicted intracranial pressure (ICP) values with previously collected data from shock-tube experiments performed on cadaver heads. In addition, to quantify the influence of including a more comprehensive network of brain vessels, we compared the biomechanical responses of our detailed-vasculature model with those of a reduced-vasculature model and a no-vasculature model for the same blast-loading conditions. For the three BOPs, the predicted ICP values matched well with the experimental results in the frontal lobe, with peak-pressure differences of 4-11% and phase-shift differences of 9-13%. As expected, incorporating the detailed cerebral vasculature did not influence the ICP, however, it redistributed the peak brain-tissue strains by as much as 30% and yielded peak strain differences of up to 7%. When compared to existing reduced-vasculature FE models that only include the major cerebral veins, our high-fidelity model redistributed the brain-tissue strains in most of the brain, highlighting the importance of including a detailed cerebral vessel network in human-head FE models to more comprehensively account for the biomechanical responses induced by blast exposure.

Mild blast wave exposure produces intensity-dependent changes in MMP2 expression patches in rat brains - Findings from different blast severities

Matrix metalloproteinase 2 (MMP2) is a gelatinase with multiple functions at the neurovascular interface, including local modification of the glia limitans to facilitate access of immune cells into the brain and amyloid-beta degradation during responses to injury or disease. This study examines regional changes in immunoreactive MMP2 in the rat brain after a single mild (2.7-7.9 psi peak) or moderate (13-17.5 psi peak) blast overpressure (BOP) exposure. Immunopositive MMP2 expression was examined quantitatively in histological sections of decalcified rat heads as a marker at 2, 24, and 72 h after BOP. The MMP2 immunoreactivity was isolated to patchy deposits in brain parenchyma surrounding blood vessels. Separate analyses were conducted for the cerebellum, brain stem caudal to the thalamo-mesencephalic junction, and the cerebrum (including diencephalon). The deposits varied in number, size, staining homogeneity (standard deviation of immunopositive region), and a cumulative measure, the product of size, average intensity and number, as a function of blast intensity and time. The sequences of changes in MMP2 spots from sham control animals suggested that the mild BOP exposure differences normalized within 72 h. However, the responses to moderate exposure revealed a delayed response at 72 h in the subtentorial brain stem and the cerebrum, but not the cerebellum. Hence, local MMP2 responses may be a contextual biomarker for locally regulated responses to widely distributed brain injury foci.

Correlating Tissue Mechanics and Spinal Cord Injury: Patient-Specific Finite Element Models of Unilateral Cervical Contusion Spinal Cord Injury in Non-Human Primates

Non-human primate (NHP) models are the closest approximation of human spinal cord injury (SCI) available for pre-clinical trials. The NHP models, however, include broader morphological variability that can confound experimental outcomes. We developed subject-specific finite element (FE) models to quantify the relationship between impact mechanics and SCI, including the correlations between FE outcomes and tissue damage. Subject-specific models of cervical unilateral contusion SCI were generated from pre-injury MRIs of six NHPs. Stress and strain outcomes were compared with lesion histology using logit analysis. A parallel generic model was constructed to compare the outcomes of subject-specific and generic models. The FE outcomes were correlated more strongly with gray matter damage (0.29 < R2 < 0.76) than white matter (0.18 < R2 < 0.58). Maximum/minimum principal strain, Von-Mises and Tresca stresses showed the strongest correlations (0.31 < R2 < 0.76) with tissue damage in the gray matter while minimum principal strain, Von-Mises stress, and Tresca stress best predicted white matter damage (0.23 < R2 < 0.58). Tissue damage thresholds varied for each subject. The generic FE model captured the impact biomechanics in two of the four models; however, the correlations between FE outcomes and tissue damage were weaker than the subject-specific models (gray matter [0.25 < R2 < 0.69] and white matter [R2 < 0.06] except for one subject [0.26 < R2 < 0.48]). The FE mechanical outputs correlated with tissue damage in spinal cord white and gray matters, and the subject-specific models accurately mimicked the biomechanics of NHP cervical contusion impacts.

A mesoscale finite element modeling approach for understanding brain morphology and material heterogeneity effects in chronic traumatic encephalopathy

Chronic Traumatic Encephalopathy (CTE) affects a significant portion of athletes in contact sports but is difficult to quantify using clinical examinations and modeling approaches. We use an in silico approach to quantify CTE biomechanics using mesoscale Finite Element (FE) analysis that bridges with macroscale whole head FE analysis. The sulci geometry produces complex stress waves that interact with one another to create increased shear stresses at the sulci depth that are significantly larger than in analyses without sulci (from 0.5 to 18.0 kPa). Sulci peak stress concentration regions coincide with experimentally observed CTE sites documented in the literature. HighlightsSulci introduce stress localizations at their depth in the gray matterSulci stress fields interact to produce stress concentration sites in white matterDifferentiating brain tissue properties did not significantly affect peak stresses.

Simulation of blast lung injury induced by shock waves of five distances based on finite element modeling of a three-dimensional rat

Blast lung injury (BLI) caused by both military and civilian explosions has become the main cause of death for blast injury patients. By building three-dimensional (3D) models of rat explosion regions, we simulated the surface pressure of the skin and lung. The pressure distributions were performed at 5 distances from the detonation center to the center of the rat. When the distances were 40 cm, 50 cm, 60 cm, 70 cm and 80 cm, the maximum pressure of the body surface were 634.77kPa, 362.46kPa, 248.11kPa, 182.13kPa and 109.29kPa and the surfaces lung pressure ranges were 928-2916 Pa, 733-2254 Pa, 488-1236 Pa, 357-1189 Pa and 314-992 Pa. After setting 6 virtual points placed on the surface of each lung lobe model, simulated pressure measurement and corresponding pathological autopsies were then conducted to validate the accuracy of the modeling. For the both sides of the lung, when the distance were 40 cm, 50 cm and 60 cm, the Pearson's values showed strong correlations. When the distances were 70 cm and 80 cm, the Pearson's values showed weak linear correlations. This computational simulation provided dynamic anatomy as well as functional and biomechanical information.

Astrocyte Mechano-Activation by High-Rate Overpressure Involves Alterations in Structural and Junctional Proteins

Primary blast neurotrauma represents a unique injury paradigm characterized by high-rate overpressure effects on brain tissue. One major hallmark of blast neurotrauma is glial reactivity, notably prolonged astrocyte activation. This cellular response has been mainly defined in primary blast neurotrauma by increased intermediate filament expression. Because the intermediate filament networks physically interface with transmembrane proteins for junctional support, it was hypothesized that cell junction regulation is altered in the reactive phenotype as well. This would have implications for downstream transcriptional regulation via signal transduction pathways like nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB). Therefore, a custom high-rate overpressure simulator was built for in vitro testing using mechanical conditions based on intracranial pressure measurements in a rat model of blast neurotrauma. Primary rat astrocytes were exposed to isolated high-rate mechanical stimulation to study cell junction dynamics in relation to their mechano-activation. First, a time course for "classical" features of reactivity was devised by evaluation of glial fibrillary acidic protein (GFAP) and proliferating cell nuclear antigen (PCNA) expression. This was followed by gene and protein expression for both gap junction (connexins) and anchoring junction proteins (integrins and cadherins). Signal transduction analysis was carried out by nuclear localization of two molecules, NF-κB p65 and mitogen-activated protein kinase (MAPK) p38. Results indicated significant increases in connexin-43 expression and PCNA first at 24 h post-overpressure (p < 0.05), followed by structural reactivity (via increased GFAP, p < 0.05) corresponding to increased anchoring junction dynamics at 48 h post-overpressure (p < 0.05). Moreover, increased phosphorylation of focal adhesion kinase (FAK) was observed in addition to increased nuclear localization of both p65 and p38 (p < 0.05) during the period of structural reactivity. To evaluate the transcriptional activity of p65 in the nucleus, electrophoretic mobility shift assay was conducted for a binding site on the promoter region for intracellular adhesion molecule-1 (ICAM-1), an antagonist of tight junctions. A significant increase in the interaction of nuclear proteins with the NF-κB site on the ICAM-1 corresponded to increased gene and protein expression of ICAM-1 (p < 0.05). Altogether, these results indicate multiple targets and corresponding signaling pathways which involve cell junction dynamics in the mechano-activation of astrocytes following high-rate overpressure.

Do blast induced skull flexures result in axonal deformation?

Subject-specific computer models (male and female) of the human head were used to investigate the possible axonal deformation resulting from the primary phase blast-induced skull flexures. The corresponding axonal tractography was explicitly incorporated into these finite element models using a recently developed technique based on the embedded finite element method. These models were subjected to extensive verification against experimental studies which examined their pressure and displacement response under a wide range of loading conditions. Once verified, a parametric study was developed to investigate the axonal deformation for a wide range of loading overpressures and directions as well as varying cerebrospinal fluid (CSF) material models. This study focuses on early times during a blast event, just as the shock transverses the skull (< 5 milliseconds). Corresponding boundary conditions were applied to eliminate the rotation effects and the resulting axonal deformation. A total of 138 simulations were developed– 128 simulations for studying the different loading scenarios and 10 simulations for studying the effects of CSF material model variance–leading to a total of 10,702 simulation core hours. Extreme strains and strain rates along each of the fiber tracts in each of these scenarios were documented and presented here. The results suggest that the blast-induced skull flexures result in strain rates as high as 150–378 s^-1. These high-strain rates of the axonal fiber tracts, caused by flexural displacement of the skull, could lead to a rate dependent micro-structural axonal damage, as pointed by other researchers.

Primary blast waves induced brain dynamics influenced by head orientations

There is controversy regarding the directional dependence of head responses subjected to blast loading. The goal of this work is to characterize the role of head orientation in the mechanics of blast wave-head interactions as well as the load transmitting to the brain. A three-dimensional human head model with anatomical details was reconstructed from computed tomography images. Three different head orientations with respect to the oncoming blast wave, i.e., front-on with head facing blast, back-on with head facing away from blast, and side-on with right side exposed to blast, were considered. The reflected pressure at the blast wave-head interface positively correlated with the skull curvature. It is evidenced by the maximum reflected pressure occurring at the eye socket with the largest curvature on the skull. The reflected pressure pattern along with the local skull areas could further influence the intracranial pressure distributions within the brain. We did find out that the maximum coup pressure of 1.031 MPa in the side-on case as well as the maximum contrecoup pressure of -0.124 MPa in the back-on case. Moreover, the maximum principal strain (MPS) was also monitored due to its indication to diffuse brain injury. It was observed that the peak MPS located in the frontal cortex region regardless of the head orientation. However, the local peak MPS within each individual function region of the brain depended on the head orientation. The detailed interactions between blast wave and head orientations provided insights for evaluating the brain dynamics, as well as biomechanical factors leading to traumatic brain injury.

Blast Scaling Parameters: Transitioning from Lung to Skull Base Metrics

Neurotrauma from blast exposure is one of the single most characteristic injuries of modern warfare. Understanding blast traumatic brain injury is critical for developing new treatment options for warfighters and civilians exposed to improvised explosive devices. Unfortunately, the pre-clinical models that are widely utilized to investigate blast exposure are based on archaic lung based parameters developed in the early 20th century. Improvised explosive devices produce a different type of injury paradigm than the typical mortar explosion. Protective equipment for the chest cavity has also improved over the past 100 years. In order to improve treatments, it is imperative to develop models that are based more on skull-based parameters. In this mini-review, we discuss the important anatomical and biochemical features necessary to develop a skull-based model.

Biomechanical Responses of the Brain in Swine Subject to Free-Field Blasts

Blast-induced traumatic brain injury (bTBI) is a signature wound of modern warfare. The current incomplete understanding of its injury mechanism impedes the development of strategies for effective protection of bTBI. Despite a considerable amount of experimental animal studies focused on the evaluation of brain neurotrauma caused by blast exposure, there is very limited knowledge on the biomechanical responses of the gyrenecephalic brain subjected to primary free-field blast waves imposed in vivo. This study aims to evaluate the external and internal mechanical responses of the brain against different levels of blast loading with Yucatan swine in free field. The incident overpressure (IOP) was generated using 3.6 kg of C4 charge placed at three standoff distances from the swine. Five swine were exposed to a total of 19 blasts. The three average peak IOP pressure levels in this study were 148.8, 278.9, and 409.2 kPa as measured by a pencil probe. The duration of the first positive wave was in the range of 2.1–3 ms. Pressure changes in the brain and head kinematics were recorded with intracranial pressure (ICP) sensors, linear accelerometers, and angular rate sensors. The corresponding average peak ICPs were in the range of 79–143, 210–281, and 311–414 kPa designated as low, medium, and high blast level, respectively. Peak head linear accelerations were in the range of 120–412 g. A positive correlation between IOP and its corresponding biomechanical responses of the brain was also observed. These experimental data can be used to validate computer models of bTBI.

Posttraumatic dizziness and vertigo

Traumatic brain injury is an increasingly common public health issue, with the mild variant most clinically relevant for this chapter. Common causes of mild traumatic brain injury (mTBI) include motor vehicle accidents, athletics, and military training/deployment. Despite a range of clinically available testing platforms, diagnosis of mTBI remains challenging. Symptoms are primarily neurosensory, and include dizziness, hearing problems, headaches, cognitive, and sleep disturbances. Dizziness is nearly universally present in all mTBI patients, and is the easiest symptom to objectify for diagnosis. Aside from a thorough history and physical exam, in the near future specialized vestibular function tests will be key to mTBI diagnosis. A battery of oculomotor (antisaccade, predictive saccade) and vestibular tasks (head impulse test) has been demonstrated to sensitively and specifically identify individuals with acute mTBI. Vestibular therapy and rehabilitation have shown improvements for mTBI patients in cognitive function, ability to return to activities of daily living, and ability to return to work. Dizziness, as a contributor to short- and long-term disability following mTBI, is ultimately crucial not only for diagnosis but also for treatment.

Relevance of Blood Vessel Networks in Blast-Induced Traumatic Brain Injury

Cerebral vasculature is a complex network that circulates blood through the brain. However, the role of this networking effect in brain dynamics has seldom been inspected. This work is to study the effects of blood vessel networks on dynamic responses of the brain under blast loading. Voronoi tessellations were implemented to represent the network of blood vessels in the brain. The brain dynamics in terms of maximum principal strain (MPS), shear strain (SS), and intracranial pressure (ICP) were monitored and compared. Results show that blood vessel networks significantly affected brain responses. The increased MPS and SS were observed within the brain embedded with vessel networks, which did not exist in the case without blood vessel networks. It is interesting to observe that the alternation of the ICP response was minimal. Moreover, the vessel diameter and density also affected brain dynamics in both MPS and SS measures. This work sheds light on the role of cerebral vasculature in blast-induced traumatic brain injury.