Towards patient-specific biomechanical human brain models

Abstract

Reliable characterization of spatial variations in brain tissue stiffness is essential for predictive biomechanical modeling, yet most current methods rely on coarse regional parameter assignments based on invasive mechanical testing. In this study, we propose a new approach to obtain subject-specific mechanical properties at voxel resolution from in vivo diffusion tensor magnetic resonance imaging (DTI) based on a linear regression between the fractional anisotropy (FA) from DTI and experimentally measured stiffness values. To assess how such heterogeneity in mechanical properties influences simulated brain deformation, we construct a finite element model based on two material parameterizations of the same human brain: one employing nine anatomically defined regions, each with uniform material parameters, and another in which the shear modulus is assigned voxel-wise on the corresponding FA value. Applying this FA-stiffness mapping yields a smoothly varying mechanical property distribution that better captures local microstructural differences not represented by region-wise parameterizations. Both parameterizations are subjected to an identical atrophy-driven loading scenario. They exhibit comparable overall volume loss, but diverge in regional behavior. The voxel-resolved parameterization predicts more pronounced ventricular expansion and differs in the displacement and stretch distributions, indicating that variability in stiffness can alter local predicted responses even when global outcomes appear similar. This work presents a pipeline for estimating individualized mechanical properties directly from imaging protocols that are routinely performed on patients, with important implications for brain biomechanics. While the approach depends on a simplified linear FA-stiffness relation and assumes isotropic constitutive behavior, it provides a framework for integrating imaging-based microstructure into subject-specific simulations. Future validation against in vivo or experimental deformation data is needed to determine the fidelity and clinical utility of FA-derived stiffness fields.