Natural oscillatory modes of 3D deformation of the human brain in vivo.

Natural modes and frequencies of three-dimensional (3D) deformation of the human brain were identified from in vivo tagged magnetic resonance images (MRI) acquired dynamically during transient mild acceleration of the head. Twenty 3D strain fields, estimated from tagged MRI image volumes in 19 adult subjects, were analyzed using dynamic mode decomposition (DMD). These strain fields represented dynamic, 3D brain deformations during constrained head accelerations, either involving rotation about the vertical axis of the neck or neck extension. DMD results reveal fundamental oscillatory modes of deformation at damped frequencies near 7 Hz (in neck rotation) and 11 Hz (in neck extension). Modes at these frequencies were found consistently among all subjects. These characteristic features of 3D human brain deformation are important for understanding the response of the brain in head impacts and provide valuable quantitative criteria for the evaluation and use of computer models of brain mechanics.

A 3D Computational Head Model Under Dynamic Head Rotation and Head Extension Validated Using Live Human Brain Data, Including the Falx and the Tentorium.

We employ an advanced 3D computational model of the head with high anatomical fidelity, together with measured tissue properties, to assess the consequences of dynamic loading to the head in two distinct modes: head rotation and head extension. We use a subject-specific computational head model, using the material point method, built from T1 magnetic resonance images, and considering the anisotropic properties of the white matter which can predict strains in the brain under large rotational accelerations. The material model now includes the shear anisotropy of the white matter. We validate the model under head rotation and head extension motions using live human data, and advance a prior version of the model to include biofidelic falx and tentorium. We then examine the consequences of incorporating the falx and tentorium in terms of the predictions from the computational head model.

Mechanics of cortical folding: stress, growth and stability.

Cortical folding, or gyrification, coincides with several important developmental processes. The folded shape of the human brain allows the cerebral cortex, the thin outer layer of neurons and their associated projections, to attain a large surface area relative to brain volume. Abnormal cortical folding has been associated with severe neurological, cognitive and behavioural disorders, such as epilepsy, autism and schizophrenia. However, despite decades of study, the mechanical forces that lead to cortical folding remain incompletely understood. Leading hypotheses have focused on the roles of (i) tangential growth of the outer cortex, (ii) spatio-temporal patterns in the birth and migration of neurons, and (iii) internal tension in axons. Recent experimental studies have illuminated not only the fundamental cellular and molecular processes underlying cortical development, but also the stress state, mechanical properties and spatio-temporal patterns of growth in the developing brain. The combination of mathematical modelling and physical measurements has allowed researchers to evaluate hypothesized mechanisms of folding, to determine whether each is consistent with physical laws. This review summarizes what physical scientists have learned from models and recent experimental observations, in the context of recent neurobiological discoveries regarding cortical development. Here, we highlight evidence of a combined mechanism, in which spatio-temporal patterns bias the locations of primary folds (i), but tangential growth of the cortical plate induces mechanical instability (ii) to propagate primary and higher-order folds.This article is part of the Theo Murphy meeting issue 'Mechanics of development'.

modCHIMERA: a novel murine closed-head model of moderate traumatic brain injury.

Traumatic brain injury is a major source of global disability and mortality. Preclinical TBI models are a crucial component of therapeutic investigation. We report a tunable, monitored model of murine non-surgical, diffuse closed-head injury-modCHIMERA-characterized by impact as well as linear and rotational acceleration. modCHIMERA is based on the Closed-Head Impact Model of Engineered Rotational Acceleration (CHIMERA) platform. We tested this model at 2 energy levels: 1.7 and 2.1 Joules-substantially higher than previously reported for this system. Kinematic analysis demonstrated linear acceleration exceeding injury thresholds in humans, although outcome metrics tracked impact energy more closely than kinematic parameters. Acute severity metrics were consistent with a complicated-mild or moderate TBI, a clinical population characterized by high morbidity but potentially reversible pathology. Axonal injury was multifocal and bilateral, neuronal death was detected in the hippocampus, and microglial neuroinflammation was prominent. Acute functional analysis revealed prolonged post-injury unconsciousness, and decreased spontaneous behavior and stimulated neurological scores. Neurobehavioral deficits were demonstrated in spatial learning/memory and socialization at 1-month. The overall injury profile of modCHIMERA corresponds with the range responsible for a substantial portion of TBI-related disability in humans. modCHIMERA should provide a reliable platform for efficient analysis of TBI pathophysiology and testing of treatment modalities.

Frequency-dependent shear properties of annulus fibrosus and nucleus pulposus by magnetic resonance elastography.

Aging and degeneration are associated with changes in mechanical properties in the intervertebral disc, generating interest in the establishment of mechanical properties as early biomarkers for the degenerative cascade. Magnetic resonance elastography (MRE) of the intervertebral disc is usually limited to the nucleus pulposus, as the annulus fibrosus is stiffer and less hydrated. The objective of this work was to adapt high-frequency needle MRE to the characterization of the shear modulus of both the nucleus pulposus and annulus fibrosus. Bovine intervertebral discs were removed from fresh oxtails and characterized by needle MRE. The needle was inserted in the center of the disc and vibrations were generated by an amplified piezoelectric actuator. MRE acquisitions were performed on a 4.7-T small-animal MR scanner using a spin echo sequence with sinusoidal motion encoding gradients. Acquisitions were repeated over a frequency range of 1000-1800 Hz. The local frequency estimation inversion algorithm was used to compute the shear modulus. Stiffness maps allowed the visualization of the soft nucleus pulposus surrounded by the stiffer annulus fibrosus surrounded by the homogeneous gel. A significant difference in shear modulus between the nucleus pulposus and annulus fibrosus, and an increase in the shear modulus with excitation frequency, were observed, in agreement with the literature. This study demonstrates that global characterization of both the nucleus pulposus and annulus fibrosus of the intervertebral disc is possible with needle MRE using a preclinical magnetic resonance imaging (MRI) scanner. MRE can be a powerful method for the mapping of the complex properties of the intervertebral disc. The developed method could be adapted for in situ use by preserving adjacent vertebrae and puncturing the side of the intervertebral disc, thereby allowing an assessment of the contribution of osmotic pressure to the mechanical behavior of the intervertebral disc.

Pre-clinical MR elastography: Principles, techniques, and applications.

Magnetic resonance elastography (MRE) is a method for measuring the mechanical properties of soft tissue in vivo, non-invasively, by imaging propagating shear waves in the tissue. The speed and attenuation of waves depends on the elastic and dissipative properties of the underlying material. Tissue mechanical properties are essential for biomechanical models and simulations, and may serve as markers of disease, injury, development, or recovery. MRE is already established as a clinical technique for detecting and characterizing liver disease. The potential of MRE for diagnosing or characterizing disease in other organs, including brain, breast, and heart is an active research area. Studies involving MRE in the pre-clinical setting, in phantoms and artificial biomaterials, in the mouse, and in other mammals, are critical to the development of MRE as a robust, reliable, and useful modality.

Measurement of anisotropic mechanical properties in porcine brain white matter ex vivo using magnetic resonance elastography.

The mechanical properties of brain tissue, particularly those of white matter (WM), need to be characterized accurately for use in finite element (FE) models of brain biomechanics and traumatic brain injury (TBI). Magnetic resonance elastography (MRE) is a powerful tool for non-invasive estimation of the mechanical properties of soft tissues. While several studies involving direct mechanical tests of brain tissue have shown mechanical anisotropy, most MRE studies of brain tissue assume an isotropic model. In this study, an incompressible transversely isotropic (TI) material model parameterized by minimum shear modulus (µ2), shear anisotropy parameter (ϕ), and tensile anisotropy parameter (ζ) is applied to analyze MRE measurements of ex vivo porcine white matter (WM) brain tissue. To characterize shear anisotropy, "slow" (pure transverse) shear waves were propagated at 100, 200 and 300Hz through sections of ex vivo brain tissue including both WM and gray matter (GM). Shear waves were found to propagate with elliptical fronts, consistent with TI material behavior. Shear wave fields were also analyzed within regions of interest (ROI) to find local shear wavelengths parallel and perpendicular to fiber orientation. FE simulations of a TI material with a range of plausible shear modulus (µ2) and shear anisotropy parameters (ϕ) were run and the results were analyzed in the same fashion as the experimental case. Parameters of the FE simulations which most closely matched each experiment were taken to represent the mechanical properties of that particular sample. Using this approach, WM in the ex vivo porcine brain was found to be mildly anisotropic in shear with estimates of minimum shear modulus (actuation frequencies listed in parenthesis): µ2= 1.04 ± 0.12 kPa (at 100Hz), µ2= 1.94 ± 0.29 kPa (at 200Hz), and µ2= 2.88 ± 0.34 kPa (at 300Hz) and corresponding shear anisotropy factors of ϕ= 0.27 ± 0.09 (at 100Hz), ϕ= 0.29 ± 0.14 (at 200Hz) and ϕ= 0.34 ± 0.13 (at 300Hz). Future MRE studies will focus on tensile anisotropy, which will require both slow and fast shear waves for accurate estimation.

Actin growth profile in clathrin-mediated endocytosis.

Clathrin-mediated endocytosis in yeast is driven by a protein patch containing close to 100 different types of proteins. Among the proteins are 5000-10000 copies of polymerized actin, and successful endocytosis requires growth of the actin network. Since it is not known exactly how actin network growth drives endocytosis, we calculate the spatial distribution of actin growth required to generate the force that drives the process. First, we establish the force distribution that must be supplied by actin growth, by combining membrane-bending profiles obtained via electron microscopy with established theories of membrane mechanics. Next, we determine the profile of actin growth, using a continuum mechanics approach and an iterative procedure starting with an actin growth profile obtained from a linear analysis. The profile has fairly constant growth outside a central hole of radius 45-50 nm, but very little growth in this hole. This growth profile can reproduce the required forces if the actin shear modulus exceeds 80 kPa, and the growing filaments can exert very large polymerization forces. The growth profile prediction could be tested via electron-microscopy or super-resolution experiments in which the turgor pressure is suddenly turned off.

Requirements for accurate estimation of anisotropic material parameters by magnetic resonance elastography: A computational study.

PURPOSE: To establish the essential requirements for characterization of a transversely isotropic material by magnetic resonance elastography (MRE). THEORY AND METHODS: Three methods for characterizing nearly incompressible, transversely isotropic (ITI) materials were used to analyze data from closed-form expressions for traveling waves, finite-element (FE) simulations of waves in homogeneous ITI material, and FE simulations of waves in heterogeneous material. Key properties are the complex shear modulus µ2 , shear anisotropy ϕ=µ1/µ2-1, and tensile anisotropy ζ=E1/E2-1. RESULTS: Each method provided good estimates of ITI parameters when both slow and fast shear waves with multiple propagation directions were present. No method gave accurate estimates when the displacement field contained only slow shear waves, only fast shear waves, or waves with only a single propagation direction. Methods based on directional filtering are robust to noise and include explicit checks of propagation and polarization. Curl-based methods led to more accurate estimates in low noise conditions. Parameter estimation in heterogeneous materials is challenging for all methods. CONCLUSIONS: Multiple shear waves, both slow and fast, with different propagation directions, must be present in the displacement field for accurate parameter estimates in ITI materials. Experimental design and data analysis can ensure that these requirements are met. Magn Reson Med 78:2360-2372, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Steady dynein forces induce flutter instability and propagating waves in mathematical models of flagella.

Cilia and flagella are highly conserved organelles that beat rhythmically with propulsive, oscillatory waveforms. The mechanism that produces these autonomous oscillations remains a mystery. It is widely believed that dynein activity must be dynamically regulated (switched on and off, or modulated) on opposite sides of the axoneme to produce oscillations. A variety of regulation mechanisms have been proposed based on feedback from mechanical deformation to dynein force. In this paper, we show that a much simpler interaction between dynein and the passive components of the axoneme can produce coordinated, propulsive oscillations. Steady, distributed axial forces, acting in opposite directions on coupled beams in viscous fluid, lead to dynamic structural instability and oscillatory, wave-like motion. This 'flutter' instability is a dynamic analogue to the well-known static instability, buckling. Flutter also occurs in slender beams subjected to tangential axial loads, in aircraft wings exposed to steady air flow and in flexible pipes conveying fluid. By analysis of the flagellar equations of motion and simulation of structural models of flagella, we demonstrate that dynein does not need to switch direction or inactivate to produce autonomous, propulsive oscillations, but must simply pull steadily above a critical threshold force.

A longitudinal magnetic resonance elastography study of murine brain tumors following radiation therapy.

An accurate and noninvasive method for assessing treatment response following radiotherapy is needed for both treatment monitoring and planning. Measurement of solid tumor volume alone is not sufficient for reliable early detection of therapeutic response, since changes in physiological and/or biomechanical properties can precede tumor volume change following therapy. In this study, we use magnetic resonance elastography to evaluate the treatment effect after radiotherapy in a murine brain tumor model. Shear modulus was calculated and compared between the delineated tumor region of interest (ROI) and its contralateral, mirrored counterpart. We also compared the shear modulus from both the irradiated and non-irradiated tumor and mirror ROIs longitudinally, sampling four time points spanning 9-19 d post tumor implant. Results showed that the tumor ROI had a lower shear modulus than that of the mirror ROI, independent of radiation. The shear modulus of the tumor ROI decreased over time for both the treated and untreated groups. By contrast, the shear modulus of the mirror ROI appeared to be relatively constant for the treated group, while an increasing trend was observed for the untreated group. The results provide insights into the tumor properties after radiation treatment and demonstrate the potential of using the mechanical properties of the tumor as a biomarker. In future studies, more closely spaced time points will be employed for detailed analysis of the radiation effect.

Magnetic resonance elastography of slow and fast shear waves illuminates differences in shear and tensile moduli in anisotropic tissue.

Mechanical anisotropy is an important property of fibrous tissues; for example, the anisotropic mechanical properties of brain white matter may play a key role in the mechanics of traumatic brain injury (TBI). The simplest anisotropic material model for small deformations of soft tissue is a nearly incompressible, transversely isotropic (ITI) material characterized by three parameters: minimum shear modulus (µ), shear anisotropy (ϕ=µ1µ-1) and tensile anisotropy (ζ=E1E2-1). These parameters can be determined using magnetic resonance elastography (MRE) to visualize shear waves, if the angle between the shear-wave propagation direction and fiber direction is known. Most MRE studies assume isotropic material models with a single shear (µ) or tensile (E) modulus. In this study, two types of shear waves, "fast" and "slow", were analyzed for a given propagation direction to estimate anisotropic parameters µ, ϕ, and ζ in two fibrous soft materials: turkey breast ex vivo and aligned fibrin gels. As expected, the speed of slow shear waves depended on the angle between fiber direction and propagation direction. Fast shear waves were observed when the deformations due to wave motion induced stretch in the fiber direction. Finally, MRE estimates of anisotropic mechanical properties in turkey breast were compared to estimates from direct mechanical tests.

Extension of the Optimized Virtual Fields Method to estimate viscoelastic material parameters from 3D dynamic displacement fields.

In-vivo measurement of the mechanical properties of soft tissues is essential to provide necessary data in biomechanics and medicine (early cancer diagnosis, study of traumatic brain injuries, etc.). Imaging techniques such as Magnetic Resonance Elastography (MRE) can provide 3D displacement maps in the bulk and in vivo, from which, using inverse methods, it is then possible to identify some mechanical parameters of the tissues (stiffness, damping etc.). The main difficulties in these inverse identification procedures consist in dealing with the pressure waves contained in the data and with the experimental noise perturbing the spatial derivatives required during the processing. The Optimized Virtual Fields Method (OVFM) [1], designed to be robust to noise, present natural and rigorous solution to deal with these problems. The OVFM has been adapted to identify material parameter maps from Magnetic Resonance Elastography (MRE) data consisting of 3-dimensional displacement fields in harmonically loaded soft materials. In this work, the method has been developed to identify elastic and viscoelastic models. The OVFM sensitivity to spatial resolution and to noise has been studied by analyzing 3D analytically simulated displacement data. This study evaluates and describes the OVFM identification performances: different biases on the identified parameters are induced by the spatial resolution and experimental noise. The well-known identification problems in the case of quasi-incompressible materials also find a natural solution in the OVFM. Moreover, an a posteriori criterion to estimate the local identification quality is proposed. The identification results obtained on actual experiments are briefly presented.

Analysis of unstable modes distinguishes mathematical models of flagellar motion.

The mechanisms underlying the coordinated beating of cilia and flagella remain incompletely understood despite the fundamental importance of these organelles. The axoneme (the cytoskeletal structure of cilia and flagella) consists of microtubule doublets connected by passive and active elements. The motor protein dynein is known to drive active bending, but dynein activity must be regulated to generate oscillatory, propulsive waveforms. Mathematical models of flagellar motion generate quantitative predictions that can be analysed to test hypotheses concerning dynein regulation. One approach has been to seek periodic solutions to the linearized equations of motion. However, models may simultaneously exhibit both periodic and unstable modes. Here, we investigate the emergence and coexistence of unstable and periodic modes in three mathematical models of flagellar motion, each based on a different dynein regulation hypothesis: (i) sliding control; (ii) curvature control and (iii) control by interdoublet separation (the 'geometric clutch' (GC)). The unstable modes predicted by each model are used to critically evaluate the underlying hypothesis. In particular, models of flagella with 'sliding-controlled' dynein activity admit unstable modes with non-propulsive, retrograde (tip-to-base) propagation, sometimes at the same parameter values that lead to periodic, propulsive modes. In the presence of these retrograde unstable modes, stable or periodic modes have little influence. In contrast, unstable modes of the GC model exhibit switching at the base and propulsive base-to-tip propagation.

Mechanical forces in cerebral cortical folding: a review of measurements and models.

Folding of the cerebral cortical surface is a critical process in human brain development, yet despite decades of indirect study and speculation the mechanics of the process remain incompletely understood. Leading hypotheses have focused on the roles of circumferential expansion of the cortex, radial growth, and internal tension in neuronal fibers (axons). In this article, we review advances in the mathematical modeling of growth and morphogenesis and new experimental data, which together promise to clarify the mechanical basis of cortical folding. Recent experimental studies have illuminated not only the fundamental cellular and molecular processes underlying cortical development, but also the stress state and mechanical behavior of the developing brain. The combination of mathematical modeling and biomechanical data provides a means to evaluate hypothesized mechanisms objectively and quantitatively, and to ensure that they are consistent with physical law, given plausible assumptions and reasonable parameter values.

A cortical folding model incorporating stress-dependent growth explains gyral wavelengths and stress patterns in the developing brain.

In humans and many other mammals, the cortex (the outer layer of the brain) folds during development. The mechanics of folding are not well understood; leading explanations are either incomplete or at odds with physical measurements. We propose a mathematical model in which (i) folding is driven by tangential expansion of the cortex and (ii) deeper layers grow in response to the resulting stress. In this model the wavelength of cortical folds depends predictably on the rate of cortical growth relative to the rate of stress-induced growth. We show analytically and in simulations that faster cortical expansion leads to shorter gyral wavelengths; slower cortical expansion leads to long wavelengths or even smooth (lissencephalic) surfaces. No inner or outer (skull) constraint is needed to produce folding, but initial shape and mechanical heterogeneity influence the final shape. The proposed model predicts patterns of stress in the tissue that are consistent with experimental observations.

Viscoelastic properties of the ferret brain measured in vivo at multiple frequencies by magnetic resonance elastography.

Characterization of the dynamic mechanical behavior of brain tissue is essential for understanding and simulating the mechanisms of traumatic brain injury (TBI). Changes in mechanical properties may also reflect changes in the brain due to aging or disease. In this study, we used magnetic resonance elastography (MRE) to measure the viscoelastic properties of ferret brain tissue in vivo. Three-dimensional (3D) displacement fields were acquired during wave propagation in the brain induced by harmonic excitation of the skull at 400 Hz, 600 Hz and 800 Hz. Shear waves with wavelengths in the order of millimeters were clearly visible in the displacement field, in strain fields, and in the curl of displacement field (which contains no contributions from longitudinal waves). Viscoelastic parameters (storage and loss moduli) governing dynamic shear deformation were estimated in gray and white matter for these excitation frequencies. To characterize the reproducibility of measurements, two ferrets were studied on three different dates each. Estimated viscoelastic properties of white matter in the ferret brain were generally similar to those of gray matter and consistent between animals and scan dates. In both tissue types G' increased from approximately 3 kPa at 400 Hz to 7 kPa at 800 Hz and G″ increased from approximately 1 kPa at 400 Hz to 2 kPa at 800 Hz. These measurements of shear wave propagation in the ferret brain can be used to both parameterize and validate finite element models of brain biomechanics.

Elastic characterization of transversely isotropic soft materials by dynamic shear and asymmetric indentation.

The mechanical characterization of soft anisotropic materials is a fundamental challenge because of difficulties in applying mechanical loads to soft matter and the need to combine information from multiple tests. A method to characterize the linear elastic properties of transversely isotropic soft materials is proposed, based on the combination of dynamic shear testing (DST) and asymmetric indentation. The procedure was demonstrated by characterizing a nearly incompressible transversely isotropic soft material. A soft gel with controlled anisotropy was obtained by polymerizing a mixture of fibrinogen and thrombin solutions in a high field magnet (B = 11.7 T); fibrils in the resulting gel were predominantly aligned parallel to the magnetic field. Aligned fibrin gels were subject to dynamic (20-40 Hz) shear deformation in two orthogonal directions. The shear storage modulus was 1.08 ± 0. 42 kPa (mean ± std. dev.) for shear in a plane parallel to the dominant fiber direction, and 0.58 ± 0.21 kPa for shear in the plane of isotropy. Gels were indented by a rectangular tip of a large aspect ratio, aligned either parallel or perpendicular to the normal to the plane of transverse isotropy. Aligned fibrin gels appeared stiffer when indented with the long axis of a rectangular tip perpendicular to the dominant fiber direction. Three-dimensional numerical simulations of asymmetric indentation were used to determine the relationship between direction-dependent differences in indentation stiffness and material parameters. This approach enables the estimation of a complete set of parameters for an incompressible, transversely isotropic, linear elastic material.