1
|
Zhang C, Bartels L, Clansey A, Kloiber J, Bondi D, van Donkelaar P, Wu L, Rauscher A, Ji S. A computational pipeline towards large-scale and multiscale modeling of traumatic axonal injury. Comput Biol Med 2024; 171:108109. [PMID: 38364663 DOI: 10.1016/j.compbiomed.2024.108109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 01/26/2024] [Accepted: 02/04/2024] [Indexed: 02/18/2024]
Abstract
Contemporary biomechanical modeling of traumatic brain injury (TBI) focuses on either the global brain as an organ or a representative tiny section of a single axon. In addition, while it is common for a global brain model to employ real-world impacts as input, axonal injury models have largely been limited to inputs of either tension or compression with assumed peak strain and strain rate. These major gaps between global and microscale modeling preclude a systematic and mechanistic investigation of how tissue strain from impact leads to downstream axonal damage throughout the white matter. In this study, a unique subject-specific multimodality dataset from a male ice-hockey player sustaining a diagnosed concussion is used to establish an efficient and scalable computational pipeline. It is then employed to derive voxelized brain deformation, maximum principal strains and white matter fiber strains, and finally, to produce diverse fiber strain profiles of various shapes in temporal history necessary for the development and application of a deep learning axonal injury model in the future. The pipeline employs a structured, voxelized representation of brain deformation with adjustable spatial resolution independent of model mesh resolution. The method can be easily extended to other head impacts or individuals. The framework established in this work is critical for enabling large-scale (i.e., across the entire white matter region, head impacts, and individuals) and multiscale (i.e., from organ to cell length scales) modeling for the investigation of traumatic axonal injury (TAI) triggering mechanisms. Ultimately, these efforts could enhance the assessment of concussion risks and design of protective headgear. Therefore, this work contributes to improved strategies for concussion detection, mitigation, and prevention.
Collapse
Affiliation(s)
- Chaokai Zhang
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, MA, USA
| | - Lara Bartels
- Department of Pediatrics, University of British Columbia, Vancouver, BC, Canada
| | - Adam Clansey
- Department of Mechanical Engineering, University of British Columbia, Vancouver, BC, Canada
| | - Julian Kloiber
- Department of Pediatrics, University of British Columbia, Vancouver, BC, Canada
| | - Daniel Bondi
- Department of Mechanical Engineering, University of British Columbia, Vancouver, BC, Canada
| | - Paul van Donkelaar
- School of Health and Exercise Sciences, University of British Columbia, Kelowna, BC, Canada
| | - Lyndia Wu
- Department of Mechanical Engineering, University of British Columbia, Vancouver, BC, Canada
| | - Alexander Rauscher
- Department of Pediatrics, University of British Columbia, Vancouver, BC, Canada
| | - Songbai Ji
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, MA, USA; Department of Mechanical Engineering, Worcester Polytechnic Institute, Worcester, MA, USA.
| |
Collapse
|
2
|
Abstract
The brain injury modeling community has recommended improving model subject specificity and simulation efficiency. Here, we extend an instantaneous (< 1 sec) convolutional neural network (CNN) brain model based on the anisotropic Worcester Head Injury Model (WHIM) V1.0 to account for strain differences due to individual morphological variations. Linear scaling factors relative to the generic WHIM along the three anatomical axes are used as additional CNN inputs. To generate training samples, the WHIM is randomly scaled to pair with augmented head impacts randomly generated from real-world data for simulation. An estimation of voxelized peak maximum principal strain of the whole-brain is said to be successful when the linear regression slope and Pearson's correlation coefficient relative to directly simulated do not deviate from 1.0 (when identical) by more than 0.1. Despite a modest training dataset (N = 1363 vs. ∼5.7 k previously), the individualized CNN achieves a success rate of 86.2% in cross-validation for scaled model responses, and 92.1% for independent generic model testing for impacts considered as complete capture of kinematic events. Using 11 scaled subject-specific models (with scaling factors determined from pre-established regression models based on head dimensions and sex and age information, and notably, without neuroimages), the morphologically individualized CNN remains accurate for impacts that also yield successful estimations for the generic WHIM. The individualized CNN instantly estimates subject-specific and spatially detailed peak strains of the entire brain and thus, supersedes others that report a scalar peak strain value incapable of informing the location of occurrence. This tool could be especially useful for youths and females due to their anticipated greater morphological differences relative to the generic model, even without the need for individual neuroimages. It has potential for a wide range of applications for injury mitigation purposes and the design of head protective gears. The voxelized strains also allow for convenient data sharing and promote collaboration among research groups.
Collapse
Affiliation(s)
- Nan Lin
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, Massachusetts, USA
| | - Shaoju Wu
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, Massachusetts, USA
| | - Songbai Ji
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, Massachusetts, USA
- Department of Mechanical Engineering, Worcester Polytechnic Institute, Worcester, Massachusetts, USA
| |
Collapse
|
3
|
Lin N, Wu S, Wu Z, Ji S. Efficient Generation of Pretraining Samples for Developing a Deep Learning Brain Injury Model via Transfer Learning. Ann Biomed Eng 2023:10.1007/s10439-023-03354-3. [PMID: 37642795 DOI: 10.1007/s10439-023-03354-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Accepted: 08/18/2023] [Indexed: 08/31/2023]
Abstract
The large amount of training samples required to develop a deep learning brain injury model demands enormous computational resources. Here, we study how a transformer neural network (TNN) of high accuracy can be used to efficiently generate pretraining samples for a convolutional neural network (CNN) brain injury model to reduce computational cost. The samples use synthetic impacts emulating real-world events or augmented impacts generated from limited measured impacts. First, we verify that the TNN remains highly accurate for the two impact types (N = 100 each; [Formula: see text] of 0.948-0.967 with root mean squared error, RMSE, ~ 0.01, for voxelized peak strains). The TNN-estimated samples (1000-5000 for each data type) are then used to pretrain a CNN, which is further finetuned using directly simulated training samples (250-5000). An independent measured impact dataset considered of complete capture of impact event is used to assess estimation accuracy (N = 191). We find that pretraining can significantly improve CNN accuracy via transfer learning compared to a baseline CNN without pretraining. It is most effective when the finetuning dataset is relatively small (e.g., 2000-4000 pretraining synthetic or augmented samples improves success rate from 0.72 to 0.81 with 500 finetuning samples). When finetuning samples reach 3000 or more, no obvious improvement occurs from pretraining. These results support using the TNN to rapidly generate pretraining samples to facilitate a more efficient training strategy for future deep learning brain models, by limiting the number of costly direct simulations from an alternative baseline model. This study could contribute to a wider adoption of deep learning brain injury models for large-scale predictive modeling and ultimately, enhancing safety protocols and protective equipment.
Collapse
Affiliation(s)
- Nan Lin
- Department of Biomedical Engineering, Worcester Polytechnic Institute, 60 Prescott Street, Worcester, MA, 01605, USA
| | - Shaoju Wu
- Department of Biomedical Engineering, Worcester Polytechnic Institute, 60 Prescott Street, Worcester, MA, 01605, USA
| | - Zheyang Wu
- Mathematical Sciences, Worcester Polytechnic Institute, Worcester, MA, 01609, USA
| | - Songbai Ji
- Department of Biomedical Engineering, Worcester Polytechnic Institute, 60 Prescott Street, Worcester, MA, 01605, USA.
- Department of Mechanical Engineering, Worcester Polytechnic Institute, 100 Institute Road, Worcester, MA, 01609, USA.
| |
Collapse
|
4
|
Ji S, Ghajari M, Mao H, Kraft RH, Hajiaghamemar M, Panzer MB, Willinger R, Gilchrist MD, Kleiven S, Stitzel JD. Use of Brain Biomechanical Models for Monitoring Impact Exposure in Contact Sports. Ann Biomed Eng 2022; 50:1389-1408. [PMID: 35867314 PMCID: PMC9652195 DOI: 10.1007/s10439-022-02999-w] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Accepted: 06/22/2022] [Indexed: 02/03/2023]
Abstract
Head acceleration measurement sensors are now widely deployed in the field to monitor head kinematic exposure in contact sports. The wealth of impact kinematics data provides valuable, yet challenging, opportunities to study the biomechanical basis of mild traumatic brain injury (mTBI) and subconcussive kinematic exposure. Head impact kinematics are translated into brain mechanical responses through physics-based computational simulations using validated brain models to study the mechanisms of injury. First, this article reviews representative legacy and contemporary brain biomechanical models primarily used for blunt impact simulation. Then, it summarizes perspectives regarding the development and validation of these models, and discusses how simulation results can be interpreted to facilitate injury risk assessment and head acceleration exposure monitoring in the context of contact sports. Recommendations and consensus statements are presented on the use of validated brain models in conjunction with kinematic sensor data to understand the biomechanics of mTBI and subconcussion. Mainly, there is general consensus that validated brain models have strong potential to improve injury prediction and interpretation of subconcussive kinematic exposure over global head kinematics alone. Nevertheless, a major roadblock to this capability is the lack of sufficient data encompassing different sports, sex, age and other factors. The authors recommend further integration of sensor data and simulations with modern data science techniques to generate large datasets of exposures and predicted brain responses along with associated clinical findings. These efforts are anticipated to help better understand the biomechanical basis of mTBI and improve the effectiveness in monitoring kinematic exposure in contact sports for risk and injury mitigation purposes.
Collapse
Affiliation(s)
- Songbai Ji
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, MA, USA.
| | - Mazdak Ghajari
- Dyson School of Design Engineering, Imperial College London, London, UK
| | - Haojie Mao
- Department of Mechanical and Materials Engineering, Faculty of Engineering, Western University, London, ON, N6A 5B9, Canada
| | - Reuben H Kraft
- Department of Mechanical and Nuclear Engineering, Department of Biomedical Engineering, The Pennsylvania State University, University Park, PA, USA
| | - Marzieh Hajiaghamemar
- Department of Biomedical Engineering, The University of Texas at San Antonio, San Antonio, TX, USA
| | - Matthew B Panzer
- Department of Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, VA, USA
| | - Remy Willinger
- University of Strasbourg, IMFS-CNRS, 2 rue Boussingault, 67000, Strasbourg, France
| | - Michael D Gilchrist
- School of Mechanical & Materials Engineering, University College Dublin, Belfield, Dublin 4, Ireland
| | - Svein Kleiven
- Division of Neuronic Engineering, KTH Royal Institute of Technology, Hälsovägen 11C, 141 57, Huddinge, Sweden
| | - Joel D Stitzel
- Department of Biomedical Engineering, Wake Forest School of Medicine, Winston-Salem, NC, USA.
| |
Collapse
|
5
|
Fiber orientation downsampling compromises the computation of white matter tract-related deformation. J Mech Behav Biomed Mater 2022; 132:105294. [DOI: 10.1016/j.jmbbm.2022.105294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2021] [Revised: 04/13/2022] [Accepted: 05/21/2022] [Indexed: 11/18/2022]
|
6
|
Wu S, Zhao W, Ji S. Real-time dynamic simulation for highly accurate spatiotemporal brain deformation from impact. COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 2022; 394:114913. [PMID: 35572209 PMCID: PMC9097909 DOI: 10.1016/j.cma.2022.114913] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
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 2 >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 2 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.
Collapse
Affiliation(s)
- Shaoju Wu
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, MA, United States of America
| | - Wei Zhao
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, MA, United States of America
| | - Songbai Ji
- Department of Biomedical Engineering, Worcester Polytechnic Institute, Worcester, MA, United States of America
- Department of Mechanical Engineering, Worcester Polytechnic Institute, Worcester, MA, United States of America
- Correspondence to: Department of Biomedical Engineering, Worcester Polytechnic Institute, 60 Prescott Street, Worcester, MA 01506, USA., (S. Ji)
| |
Collapse
|
7
|
|