On the Role of Mechanics in Chronic Lung Disease.

Progressive airflow obstruction is a classical hallmark of chronic lung disease, affecting more than one fourth of the adult population. As the disease progresses, the inner layer of the airway wall grows, folds inwards, and narrows the lumen. The critical failure conditions for airway folding have been studied intensely for idealized circular cross-sections. However, the role of airway branching during this process is unknown. Here, we show that the geometry of the bronchial tree plays a crucial role in chronic airway obstruction and that critical failure conditions vary significantly along a branching airway segment. We perform systematic parametric studies for varying airway cross-sections using a computational model for mucosal thickening based on the theory of finite growth. Our simulations indicate that smaller airways are at a higher risk of narrowing than larger airways and that regions away from a branch narrow more drastically than regions close to a branch. These results agree with clinical observations and could help explain the underlying mechanisms of progressive airway obstruction. Understanding growth-induced instabilities in constrained geometries has immediate biomedical applications beyond asthma and chronic bronchitis in the diagnostics and treatment of chronic gastritis, obstructive sleep apnea and breast cancer

Mathematical modeling of collagen turnover in biological tissue

The final publication is available at Springer via http://dx.doi.org/10.1007/s00285-012-0613-yWe present a theoretical and computational model for collagen turnover in soft biological tissues. Driven by alterations in the mechanical environment, collagen fiber bundles may undergo important chronic changes, characterized primarily by alterations in collagen synthesis and degradation rates. In particular, hypertension triggers an increase in tropocollagen synthesis and a decrease in collagen degradation, which lead to the well-documented overall increase in collagen content. These changes are the result of a cascade of events, initiated mainly by the endothelial and smooth muscle cells. Here, we represent these events collectively in terms of two internal variables, the concentration of growth factor TGF-$\beta$ and tissue inhibitors of metalloproteinases TIMP. The upregulation of TGF-$\beta$ increases the collagen density. The upregulation of TIMP also increases the collagen density through decreasing matrix metalloproteinase MMP. We establish a mathematical theory for mechanically-induced collagen turnover and introduce a computational algorithm for its robust and efficient solution. We demonstrate that our model can accurately predict the experimentally observed collagen increase in response to hypertension reported in literature. Ultimately, the model can serve as a valuable tool to predict the chronic adaptation of collagen content to restore the homeostatic equilibrium state in vessels with arbitrary micro-structure and geometry.Peer ReviewedPostprint (author's final draft

Computational modeling of hypertensive growth in the human carotid artery

The final publication is available at Springer via http://dx.doi.org/10.1007/s00466-013-0959-zArterial hypertension is a chronic medical condition associated with an elevated blood pressure. Chronic arterial hypertension initiates a series of events, which are known to collectively initiate arterial wall thickening. However, the correlation between macrostructural mechanical loading, microstructural cellular changes, and macrostructural adaptation remains unclear. Here, we present a microstructurally motivated computational model for chronic arterial hypertension through smooth muscle cell growth. To model growth, we adopt a classical concept based on the multiplicative decomposition of the deformation gradient into an elastic part and a growth part. Motivated by clinical observations, we assume that the driving force for growth is the stretch sensed by the smooth muscle cells. We embed our model into a finite element framework, where growth is stored locally as an internal variable. First, to demonstrate the features of our model, we investigate the effects of hypertensive growth in a real human carotid artery. Our results agree nicely with experimental data reported in the literature both qualitatively and quantitatively.Peer ReviewedPostprint (author's final draft

Physical biology of human brain development

Neurodevelopment is a complex, dynamic process that involves a precisely orchestrated sequence of genetic, environmental, biochemical, and physical events. Developmental biology and genetics have shaped our understanding of the molecular and cellular mechanisms during neurodevelopment. Recent studies suggest that physical forces play a central role in translating these cellular mechanisms into the complex surface morphology of the human brain. However, the precise impact of neuronal differentiation, migration, and connection on the physical forces during cortical folding remains unknown. Here we review the cellular mechanisms of neurodevelopment with a view toward surface morphogenesis, pattern selection, and evolution of shape. We revisit cortical folding as the instability problem of constrained differential growth in a multi-layered system. To identify the contributing factors of differential growth, we map out the timeline of neurodevelopment in humans and highlight the cellular events associated with extreme radial and tangential expansion. We demonstrate how computational modeling of differential growth can bridge the scales-from phenomena on the cellular level toward form and function on the organ level-to make quantitative, personalized predictions. Physics-based models can quantify cortical stresses, identify critical folding conditions, rationalize pattern selection, and predict gyral wavelengths and gyrification indices. We illustrate that physical forces can explain cortical malformations as emergent properties of developmental disorders. Combining biology and physics holds promise to advance our understanding of human brain development and enable early diagnostics of cortical malformations with the ultimate goal to improve treatment of neurodevelopmental disorders including epilepsy, autism spectrum disorders, and schizophrenia

The multiphysics of prion-like diseases: progression and atrophy

Many neurodegenerative diseases are related to the propagation and accumulation of toxic proteins throughout the brain. The lesions created by aggregates of these toxic proteins further lead to cell death and accelerated tissue atrophy. A striking feature of some of these diseases is their characteristic pattern and evolution, leading to well-codified disease stages visible to neuropathology and associated with various cognitive deficits and pathologies. Here, we simulate the anisotropic propagation and accumulation of toxic proteins in full brain geometry. We show that the same model with different initial seeding zones reproduces the characteristic evolution of different prion-like diseases. We also recover the expected evolution of the total toxic protein load. Finally, we couple our transport model to a mechanical atrophy model to obtain the typical degeneration patterns found in neurodegenerative diseases.Comment: 5 pages/5 figure

The role of mechanics during brain development

Convolutions are a classical hallmark of most mammalian brains. Brain surface morphology is often associated with intelligence and closely correlated with neurological dysfunction. Yet, we know surprisingly little about the underlying mechanisms of cortical folding. Here we identify the role of the key anatomic players during the folding process: cortical thickness, stiffness, and growth. To establish estimates for the critical time, pressure, and the wavelength at the onset of folding, we derive an analytical model using the Föppl–von Kármán theory. Analytical modeling provides a quick first insight into the critical conditions at the onset of folding, yet it fails to predict the evolution of complex instability patterns in the post-critical regime. To predict realistic surface morphologies, we establish a computational model using the continuum theory of finite growth. Computational modeling not only confirms our analytical estimates, but is also capable of predicting the formation of complex surface morphologies with asymmetric patterns and secondary folds. Taken together, our analytical and computational models explain why larger mammalian brains tend to be more convoluted than smaller brains. Both models provide mechanistic interpretations of the classical malformations of lissencephaly and polymicrogyria. Understanding the process of cortical folding in the mammalian brain has direct implications on the diagnostics of neurological disorders including severe retardation, epilepsy, schizophrenia, and autism

Theory and implementation of inelastic Constitutive Artificial Neural Networks

Nature has always been our inspiration in the research, design and development of materials and has driven us to gain a deep understanding of the mechanisms that characterize anisotropy and inelastic behavior. All this knowledge has been accumulated in the principles of thermodynamics. Deduced from these principles, the multiplicative decomposition combined with pseudo potentials are powerful and universal concepts. Simultaneously, the tremendous increase in computational performance enabled us to investigate and rethink our history-dependent material models to make the most of our predictions. Today, we have reached a point where materials and their models are becoming increasingly sophisticated. This raises the question: How do we find the best model that includes all inelastic effects to explain our complex data? Constitutive Artificial Neural Networks (CANN) may answer this question. Here, we extend the CANNs to inelastic materials (iCANN). Rigorous considerations of objectivity, rigid motion of the reference configuration, multiplicative decomposition and its inherent non-uniqueness, restrictions of energy and pseudo potential, and consistent evolution guide us towards the architecture of the iCANN satisfying thermodynamics per design. We combine feed-forward networks of the free energy and pseudo potential with a recurrent neural network approach to take time dependencies into account. We demonstrate that the iCANN is capable of autonomously discovering models for artificially generated data, the response of polymers for cyclic loading and the relaxation behavior of muscle data. As the design of the network is not limited to visco-elasticity, our vision is that the iCANN will reveal to us new ways to find the various inelastic phenomena hidden in the data and to understand their interaction. Our source code, data, and examples are available at doi.org/10.5281/zenodo.10066805Comment: 54 pages, 14 figures, 14 table