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.

Phase Vector Incompressible Registration Algorithm for Motion Estimation From Tagged Magnetic Resonance Images

Tagged magnetic resonance imaging has been used for decades to observe and quantify motion and strain of deforming tissue. It is challenging to obtain 3-D motion estimates due to a tradeoff between image slice density and acquisition time. Typically, interpolation methods are used either to combine 2-D motion extracted from sparse slice acquisitions into 3-D motion or to construct a dense volume from sparse acquisitions before image registration methods are applied. This paper proposes a new phase-based 3-D motion estimation technique that first computes harmonic phase volumes from interpolated tagged slices and then matches them using an image registration framework. The approach uses several concepts from diffeomorphic image registration with a key novelty that defines a symmetric similarity metric on harmonic phase volumes from multiple orientations. The material property of harmonic phase solves the aperture problem of optical flow and intensity-based methods and is robust to tag fading. A harmonic magnitude volume is used in enforcing incompressibility in the tissue regions. The estimated motion fields are dense, incompressible, diffeomorphic, and inverse-consistent at a 3-D voxel level. The method was evaluated using simulated phantoms, human brain data in mild head accelerations, human tongue data during speech, and an open cardiac data set. The method shows comparable accuracy to three existing methods while demonstrating low computation time and robustness to tag fading and noise.

Flexural Rigidity and Shear Stiffness of Flagella Estimated from Induced Bends and Counterbends

Motile cilia and flagella are whiplike cellular organelles that bend actively to propel cells or move fluid in passages such as airways, brain ventricles, and the oviduct. Efficient motile function of cilia and flagella depends on coordinated interactions between active forces from an array of motor proteins and passive mechanical resistance from the complex cytoskeletal structure (the axoneme). However, details of this coordination, including axonemal mechanics, remain unclear. We investigated two major mechanical parameters, flexural rigidity and interdoublet shear stiffness, of the flagellar axoneme in the unicellular alga Chlamydomonas reinhardtii. Combining experiment, theory, and finite element models, we demonstrate that the apparent flexural rigidity of the axoneme depends on both the intrinsic flexural rigidity (EI) and the elastic resistance to interdoublet sliding (shear stiffness, ks). We estimated the average intrinsic flexural rigidity and interdoublet shear stiffness of wild-type Chlamydomonas flagella in vivo, rendered immotile by vanadate, to be EI = 840 ± 280 pN⋅μm(2) and ks = 79.6 ± 10.5 pN/rad, respectively. The corresponding values for the pf3; cnk11-6 double mutant, which lacks the nexin-dynein regulatory complex (N-DRC), were EI = 1011 ± 183 pN·μm(2) and ks = 39.3 ± 6.0 pN/rad under the same conditions. Finally, in the pf13A mutant, which lacks outer dynein arms and inner dynein arm c, the estimates were EI = 777 ± 184 pN·μm(2) and ks = 43.3 ± 7.7 pN/rad. In the two mutant strains, the flexural rigidity is not significantly different from wild-type (p > 0.05), but the lack of N-DRC (in pf3; cnk11-6) or dynein arms (in pf13A) significantly reduces interdoublet shear stiffness. These differences may represent the contributions of the N-DRCs (∼40 pN/rad) and residual dynein interactions (∼35 pN/rad) to interdoublet sliding resistance in these immobilized Chlamydomonas flagella.

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.

A Three-Dimensional Computational Human Head Model That Captures Live Human Brain Dynamics

Diffuse axonal injury (DAI) is a debilitating consequence of traumatic brain injury (TBI) attributed to abnormal stretching of axons caused by blunt head trauma or acceleration of the head. We developed an anatomically accurate, subject-specific, three-dimensional (3D) computational model of the human brain, and used it to study the dynamic deformations in the substructures of the brain when the head is subjected to rotational accelerations. The computational head models use anatomy and morphology of the white matter fibers obtained using MRI. Subject-specific full-field shearing motions in live human brains obtained through a recently developed tagged MRI imaging technique are then used to validate the models by comparing the measured and predicted heterogeneous dynamic mechanical response of the brain. These results are used to elucidate the dynamics of local shearing deformations in the brain substructures caused by rotational acceleration of the head. Our work demonstrates that the rotational dynamics of the brain has a timescale of ∼100 ms as determined by the shearing wave speeds, and thus the injuries associated with rotational accelerations likely occur over these time scales. After subject-specific validation using the live human subject data, a representative subject-specific head model is used to simulate a real life scenario that resulted in a concussive injury. Results suggest that regions of the brain, in the form of a toroid, encompassing the white matter, the cortical gray matter, and outer parts of the limbic system have a higher susceptibility to injury under axial rotations of the head.

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 μ₂ , 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.

Contraction and stress-dependent growth shape the forebrain of the early chicken embryo

During early vertebrate development, local constrictions, or sulci, form to divide the forebrain into the diencephalon, telencephalon, and optic vesicles. These partitions are maintained and exaggerated as the brain tube inflates, grows, and bends. Combining quantitative experiments on chick embryos with computational modeling, we investigated the biophysical mechanisms that drive these changes in brain shape. Chemical perturbations of contractility indicated that actomyosin contraction plays a major role in the creation of initial constrictions (Hamburger-Hamilton stages HH11-12), and fluorescent staining revealed that F-actin is circumferentially aligned at all constrictions. A finite element model based on these findings shows that the observed shape changes are consistent with circumferential contraction in these regions. To explain why sulci continue to deepen as the forebrain expands (HH12-20), we speculate that growth depends on wall stress. This idea was examined by including stress-dependent growth in a model with cerebrospinal fluid pressure and bending (cephalic flexure). The results given by the model agree with observed morphological changes that occur in the brain tube under normal and reduced eCSF pressure, quantitative measurements of relative sulcal depth versus time, and previously published patterns of cell proliferation. Taken together, our results support a biphasic mechanism for forebrain morphogenesis consisting of differential contractility (early) and stress-dependent growth (late).

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.

On the accuracy and fitting of transversely isotropic material models

Fiber reinforced structures are central to the form and function of biological tissues. Hyperelastic, transversely isotropic material models are used widely in the modeling and simulation of such tissues. Many of the most widely used models involve strain energy functions that include one or both pseudo-invariants (I4 or I5) to incorporate energy stored in the fibers. In a previous study we showed that both of these invariants must be included in the strain energy function if the material model is to reduce correctly to the well-known framework of transversely isotropic linear elasticity in the limit of small deformations. Even with such a model, fitting of parameters is a challenge. Here, by evaluating the relative roles of I4 and I5 in the responses to simple loadings, we identify loading scenarios in which previous models accounting for only one of these invariants can be expected to provide accurate estimation of material response, and identify mechanical tests that have special utility for fitting of transversely isotropic constitutive models. Results provide guidance for fitting of transversely isotropic constitutive models and for interpretation of the predictions of these models.

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.

Estimation of material parameters from slow and fast shear waves in an incompressible, transversely isotropic material

This paper describes a method to estimate mechanical properties of soft, anisotropic materials from measurements of shear waves with specific polarization and propagation directions. This method is applicable to data from magnetic resonance elastography (MRE), which is a method for measuring shear waves in live subjects or in vitro samples. Here, we simulate MRE data using finite element analysis. A nearly incompressible, transversely isotropic (ITI) material model with three parameters (shear modulus, shear anisotropy, and tensile anisotropy) is used, which is appropriate for many fibrous, biological tissues. Both slow and fast shear waves travel concurrently through such a material with speeds that depend on the propagation direction relative to fiber orientation. A three-parameter estimation approach based on directional filtering and isolation of slow and fast shear wave components (directional filter inversion, or DFI) is introduced. Wave speeds of each isolated shear wave component are estimated using local frequency estimation (LFE), and material properties are calculated using weighted least squares. Data from multiple finite element simulations are used to assess the accuracy and reliability of DFI for estimation of anisotropic material parameters.

Dynein-deficient flagella respond to increased viscosity with contrasting changes in power and recovery strokes

Changes in the flagellar waveform in response to increased viscosity were investigated in uniflagellate mutants of Chlamydomonas reinhardtii. We hypothesized that the waveforms of mutants lacking different dynein arms would change in different ways as viscosity was increased, and that these variations would illuminate the feedback pathways from force to dynein activity. Previous studies have investigated the effects of viscosity on cell body motion, propulsive force, and power in different mutants, but the effect on waveform has not yet been fully characterized. Beat frequency decreases with viscosity in wild-type uniflagellate (uni1) cells, and outer dynein arm deficient (oda2) mutants. In contrast, the inner dynein arm mutant ida1 (lacking I1/f) maintains beat frequency at high viscosity but alters its flagellar waveform more than either wild-type or oda2. The ida1 waveform is narrower than wild-type, primarily due to an abbreviated recovery stroke; this difference is amplified at high viscosity. The oda2 mutant in contrast, maintains a consistent waveform at high and low viscosity with a slightly longer power stroke than wild-type. Analysis of the delays and shear displacements between bends suggest that direct force feedback in the outer dynein arm system may initiate switching of dynein activity. In contrast, I1/f dynein appears to delay switching, most markedly at the initiation of the power stroke, possibly by controlling inter-doublet separation.

Equations of interdoublet separation during flagella motion reveal mechanisms of wave propagation and instability

The motion of flagella and cilia arises from the coordinated activity of dynein motor protein molecules arrayed along microtubule doublets that span the length of axoneme (the flagellar cytoskeleton). Dynein activity causes relative sliding between the doublets, which generates propulsive bending of the flagellum. The mechanism of dynein coordination remains incompletely understood, although it has been the focus of many studies, both theoretical and experimental. In one leading hypothesis, known as the geometric clutch (GC) model, local dynein activity is thought to be controlled by interdoublet separation. The GC model has been implemented as a numerical simulation in which the behavior of a discrete set of rigid links in viscous fluid, driven by active elements, was approximated using a simplified time-marching scheme. A continuum mechanical model and associated partial differential equations of the GC model have remained lacking. Such equations would provide insight into the underlying biophysics, enable mathematical analysis of the behavior, and facilitate rigorous comparison to other models. In this article, the equations of motion for the flagellum and its doublets are derived from mechanical equilibrium principles and simple constitutive models. These equations are analyzed to reveal mechanisms of wave propagation and instability in the GC model. With parameter values in the range expected for Chlamydomonas flagella, solutions to the fully nonlinear equations closely resemble observed waveforms. These results support the ability of the GC hypothesis to explain dynein coordination in flagella and provide a mathematical foundation for comparison to other leading models.

Force generation by endocytic actin patches in budding yeast

Membrane deformation during endocytosis in yeast is driven by local, templated assembly of a sequence of proteins including polymerized actin and curvature-generating coat proteins such as clathrin. Actin polymerization is required for successful endocytosis, but it is not known by what mechanisms actin polymerization generates the required pulling forces. To address this issue, we develop a simulation method in which the actin network at the protein patch is modeled as an active gel. The deformation of the gel is treated using a finite-element approach. We explore the effects and interplay of three different types of force driving invagination: 1), forces perpendicular to the membrane, generated by differences between actin polymerization rates at the edge of the patch and those at the center; 2), the inherent curvature of the coat-protein layer; and 3), forces parallel to the membrane that buckle the coat protein layer, generated by an actomyosin contractile ring. We find that with optimistic estimates for the stall stress of actin gel growth and the shear modulus of the actin gel, actin polymerization can generate almost enough force to overcome the turgor pressure. In combination with the other mechanisms, actin polymerization can the force over the critical value.

Improved measurement of brain deformation during mild head acceleration using a novel tagged MRI sequence

In vivo measurements of human brain deformation during mild acceleration are needed to help validate computational models of traumatic brain injury and to understand the factors that govern the mechanical response of the brain. Tagged magnetic resonance imaging is a powerful, noninvasive technique to track tissue motion in vivo which has been used to quantify brain deformation in live human subjects. However, these prior studies required from 72 to 144 head rotations to generate deformation data for a single image slice, precluding its use to investigate the entire brain in a single subject. Here, a novel method is introduced that significantly reduces temporal variability in the acquisition and improves the accuracy of displacement estimates. Optimization of the acquisition parameters in a gelatin phantom and three human subjects leads to a reduction in the number of rotations from 72 to 144 to as few as 8 for a single image slice. The ability to estimate accurate, well-resolved, fields of displacement and strain in far fewer repetitions will enable comprehensive studies of acceleration-induced deformation throughout the human brain in vivo.

Strain Localization in an Oscillating Maxwell Viscoelastic Cylinder

The transient rotation responses of simple, axisymmetric, viscoelastic structures are of interest for interpretation of experiments designed to characterize materials and closed structures such as the brain using magnetic resonance techniques. Here, we studied the response of a Maxwell viscoelastic cylinder to small, sinusoidal displacement of its outer boundary. The transient strain field can be calculated in closed form using any of several conventional approaches. The solution is surprising: the strain field develops a singularity that appears when the wavefront leaves the center of the cylinder, and persists as the wavefront reflects to the outer boundary and back to the center of the cylinder. The singularity is alternately annihilated and reinitiated upon subsequent departures of the wavefront from the center of the cylinder until it disappears in the limit of steady state oscillations. We present the solution for this strain field, characterize the nature of this singularity, and discuss its potential role in the mechanical response and evolved morphology of the brain.

Mechanical properties of viscoelastic media by local frequency estimation of divergence-free wave fields

Magnetic resonance elastography (MRE) is an imaging modality with which mechanical properties can be noninvasively measured in living tissue. Magnetic resonance elastography relies on the fact that the elastic shear modulus determines the phase velocity and, hence the wavelength, of shear waves which are visualized by motion-sensitive MR imaging. Local frequency estimation (LFE) has been used to extract the local wavenumber from displacement wave fields recorded by MRE. LFE -based inversion is attractive because it allows material parameters to be estimated without explicitly invoking the equations governing wave propagation, thus obviating the need to numerically compute the Laplacian. Nevertheless, studies using LFE have not explicitly addressed three important issues: (1) tissue viscoelasticity; (2) the effects of longitudinal waves and rigid body motion on estimates of shear modulus; and (3) mechanical anisotropy. In the current study we extend the LFE technique to (1) estimate the (complex) viscoelastic shear modulus in lossy media; (2) eliminate the effects of longitudinal waves and rigid body motion; and (3) determine two distinct shear moduli in anisotropic media. The extended LFE approach is demonstrated by analyzing experimental data from a previously-characterized, isotropic, viscoelastic, gelatin phantom and simulated data from a computer model of anisotropic (transversely isotropic) soft material.

CCDC65 mutation causes primary ciliary dyskinesia with normal ultrastructure and hyperkinetic cilia

BACKGROUND

Primary ciliary dyskinesia (PCD) is a genetic disorder characterized by impaired ciliary function, leading to chronic sinopulmonary disease. The genetic causes of PCD are still evolving, while the diagnosis is often dependent on finding a ciliary ultrastructural abnormality and immotile cilia. Here we report a novel gene associated with PCD but without ciliary ultrastructural abnormalities evident by transmission electron microscopy, but with dyskinetic cilia beating.

METHODS

Genetic linkage analysis was performed in a family with a PCD subject. Gene expression was studied in Chlamydomonas reinhardtii and human airway epithelial cells, using RNA assays and immunostaining. The phenotypic effects of candidate gene mutations were determined in primary culture human tracheobronchial epithelial cells transduced with gene targeted shRNA sequences. Video-microscopy was used to evaluate cilia motion.

RESULTS

A single novel mutation in CCDC65, which created a termination codon at position 293, was identified in a subject with typical clinical features of PCD. CCDC65, an orthologue of the Chlamydomonas nexin-dynein regulatory complex protein DRC2, was localized to the cilia of normal nasal epithelial cells but was absent in those from the proband. CCDC65 expression was up-regulated during ciliogenesis in cultured airway epithelial cells, as was DRC2 in C. reinhardtii following deflagellation. Nasal epithelial cells from the affected individual and CCDC65-specific shRNA transduced normal airway epithelial cells had stiff and dyskinetic cilia beating patterns compared to control cells. Moreover, Gas8, a nexin-dynein regulatory complex component previously identified to associate with CCDC65, was absent in airway cells from the PCD subject and CCDC65-silenced cells.

CONCLUSION

Mutation in CCDC65, a nexin-dynein regulatory complex member, resulted in a frameshift mutation and PCD. The affected individual had altered cilia beating patterns, and no detectable ultrastructural defects of the ciliary axoneme, emphasizing the role of the nexin-dynein regulatory complex and the limitations of certain methods for PCD diagnosis.

Measurements of mechanical anisotropy in brain tissue and implications for transversely isotropic material models of white matter

White matter in the brain is structurally anisotropic, consisting largely of bundles of aligned, myelin-sheathed axonal fibers. White matter is believed to be mechanically anisotropic as well. Specifically, transverse isotropy is expected locally, with the plane of isotropy normal to the local mean fiber direction. Suitable material models involve strain energy density functions that depend on the I4 and I5 pseudo-invariants of the Cauchy-Green strain tensor to account for the effects of relatively stiff fibers. The pseudo-invariant I4 is the square of the stretch ratio in the fiber direction; I5 contains contributions of shear strain in planes parallel to the fiber axis. Most, if not all, published models of white matter depend on I4 but not on I5. Here, we explore the small strain limits of these models in the context of experimental measurements that probe these dependencies. Models in which strain energy depends on I4 but not I5 can capture differences in Young's (tensile) moduli, but will not exhibit differences in shear moduli for loading parallel and normal to the mean direction of axons. We show experimentally, using a combination of shear and asymmetric indentation tests, that white matter does exhibit such differences in both tensile and shear moduli. Indentation tests were interpreted through inverse fitting of finite element models in the limit of small strains. Results highlight that: (1) hyperelastic models of transversely isotropic tissues such as white matter should include contributions of both the I4 and I5 strain pseudo-invariants; and (2) behavior in the small strain regime can usefully guide the choice and initial parameterization of more general material models of white matter.

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.