Growth and remodelling of living tissues: perspectives, challenges and opportunities

Myocardial Bridge and Atherosclerosis, an Intimal Relationship

PURPOSE OF REVIEW

This review investigates the relationship between myocardial bridges (MBs), intimal thickening in coronary arteries, and Atherosclerotic cardiovascular disease. It focuses on the role of mechanical forces, such as circumferential strain, in arterial wall remodeling and aims to clarify how MBs affect coronary artery pathology.

REVIEW FINDINGS

MBs have been identified as influential in modulating coronary artery intimal thickness, demonstrating a protective effect against thickening within the MB segment and an increase in thickness proximal to the MB. This is attributed to changes in mechanical stress and hemodynamics. Research involving arterial hypertension models and vein graft disease has underscored the importance of circumferential strain in vascular remodeling and intimal hyperplasia. Understanding the complex dynamics between MBs, mechanical strain, and vascular remodeling is crucial for advancing our knowledge of coronary artery disease mechanisms. This could lead to improved management strategies for cardiovascular diseases, highlighting the need for further research into MB-related vascular changes.

Review of cardiac-coronary interaction and insights from mathematical modeling

Cardiac-coronary interaction is fundamental to the function of the heart. As one of the highest metabolic organs in the body, the cardiac oxygen demand is met by blood perfusion through the coronary vasculature. The coronary vasculature is largely embedded within the myocardial tissue which is continually contracting and hence squeezing the blood vessels. The myocardium-coronary vessel interaction is two-ways and complex. Here, we review the different types of cardiac-coronary interactions with a focus on insights gained from mathematical models. Specifically, we will consider the following: (1) myocardial-vessel mechanical interaction; (2) metabolic-flow interaction and regulation; (3) perfusion-contraction matching, and (4) chronic interactions between the myocardium and coronary vasculature. We also provide a discussion of the relevant experimental and clinical studies of different types of cardiac-coronary interactions. Finally, we highlight knowledge gaps, key challenges, and limitations of existing mathematical models along with future research directions to understand the unique myocardium-coronary coupling in the heart. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

Obstructed swelling and fracture of hydrogels

Obstructions influence the growth and expansion of bodies in a wide range of settings-but isolating and understanding their impact can be difficult in complex environments. Here, we study obstructed growth/expansion in a model system accessible to experiments, simulations, and theory: hydrogels swelling around fixed cylindrical obstacles with varying geometries. When the obstacles are large and widely-spaced, hydrogels swell around them and remain intact. In contrast, our experiments reveal that when the obstacles are narrow and closely-spaced, hydrogels fracture as they swell. We use finite element simulations to map the magnitude and spatial distribution of stresses that build up during swelling at equilibrium in a 2D model, providing a route toward predicting when this phenomenon of self-fracturing is likely to arise. Applying lessons from indentation theory, poroelasticity, and nonlinear continuum mechanics, we also develop a theoretical framework for understanding how the maximum principal tensile and compressive stresses that develop during swelling are controlled by obstacle geometry and material parameters. These results thus help to shed light on the mechanical principles underlying growth/expansion in environments with obstructions.

Tumor proliferation and invasion are intrinsically coupled and unraveled through tunable spheroid and physics-based models

Proliferation and invasion are two key drivers of tumor growth that are traditionally considered independent multicellular processes. However, these processes are intrinsically coupled through a maximum carrying capacity, i.e., the maximum spatial cell concentration supported by the tumor volume, total cell count, nutrient access, and mechanical properties of the tissue stroma. We explored this coupling of proliferation and invasion through in vitro and in silico methods where we modulated the mechanical properties of the tumor and the surrounding extracellular matrix. E-cadherin expression and stromal collagen concentration were manipulated in a tunable breast cancer spheroid to determine the overall impacts of these tumor variables on net tumor proliferation and continuum invasion. We integrated these results into a mixed-constitutive formulation to computationally delineate the influences of cellular and extracellular adhesion, stiffness, and mechanical properties of the extracellular matrix on net proliferation and continuum invasion. This framework integrates biological in vitro data into concise computational models of invasion and proliferation to provide more detailed physical insights into the coupling of these key tumor processes and tumor growth. STATEMENT OF SIGNIFICANCE: Tumor growth involves expansion into the collagen-rich stroma through intrinsic coupling of proliferation and invasion within the tumor continuum. These processes are regulated by a maximum carrying capacity that is determined by the total cell count, tumor volume, nutrient access, and mechanical properties of the surrounding stroma. The influences of biomechanical parameters (i.e., stiffness, cell elongation, net proliferation rate and cell-ECM friction) on tumor proliferation or invasion cannot be unraveled using experimental methods alone. By pairing a tunable spheroid system with computational modeling, we delineated the interdependencies of each system parameter on tumor proliferation and continuum invasion, and established a concise computational framework for studying tumor mechanobiology.

Forceful patterning: theoretical principles of mechanochemical pattern formation

Biological pattern formation is essential for generating and maintaining spatial structures from the scale of a single cell to tissues and even collections of organisms. Besides biochemical interactions, there is an important role for mechanical and geometrical features in the generation of patterns. We review the theoretical principles underlying different types of mechanochemical pattern formation across spatial scales and levels of biological organization.

Mechanical modeling of the maturation process for tissue-engineered implants: Application to biohybrid heart valves

The development of tissue-engineered cardiovascular implants can improve the lives of large segments of our society who suffer from cardiovascular diseases. Regenerative tissues are fabricated using a process called tissue maturation. Furthermore, it is highly challenging to produce cardiovascular regenerative implants with sufficient mechanical strength to withstand the loading conditions within the human body. Therefore, biohybrid implants for which the regenerative tissue is reinforced by standard reinforcement material (e.g. textile or 3d printed scaffold) can be an interesting solution. In silico models can significantly contribute to characterizing, designing, and optimizing biohybrid implants. The first step towards this goal is to develop a computational model for the maturation process of tissue-engineered implants. This paper focuses on the mechanical modeling of textile-reinforced tissue-engineered cardiovascular implants. First, an energy-based approach is proposed to compute the collagen evolution during the maturation process. Then, the concept of structural tensors is applied to model the anisotropic behavior of the extracellular matrix and the textile scaffold. Next, the newly developed material model is embedded into a special solid-shell finite element formulation with reduced integration. Finally, our framework is used to compute two structural problems: a pressurized shell construct and a tubular-shaped heart valve. The results show the ability of the model to predict collagen growth in response to the boundary conditions applied during the maturation process. Consequently, the model can predict the implant's mechanical response, such as the deformation and stresses of the implant.

Fetal head-down posture may explain the rapid brain evolution in humans and other primates: An interpretative review

Evolutionary cerebrovascular consequences of upside-down postural verticality of the anthropoid fetus have been largely overlooked in the literature. This working hypothesis-based report provides a literature interpretation from an aspect that the rapid evolution of the human brain has been promoted by fetal head-down position due to maternal upright and semi-upright posture. Habitual vertical torso posture is a feature not only of humans, but also of monkeys and non-human apes that spend considerable time in a sitting position. Consequently, the head-down position of the fetus may have caused physiological craniovascular hypertension that stimulated expansion of the intracranial vessels and acted as an epigenetic physiological stress, which enhanced neurogenesis and eventually, along with other selective pressures, led to the progressive growth of the anthropoid brain and its organization. This article collaterally opens a new insight into the conundrum of high cephalopelvic proportions (i.e., the tight fit between the pelvic birth canal and fetal head) in phylogenetically distant lineages of monkeys, lesser apes, and humans. Low cephalopelvic proportions in non-human great apes could be accounted for by their energetically efficient horizontal nest-sleeping and consequently by their larger body mass compared to monkeys and lesser apes that sleep upright. One can further hypothesize that brain size varies in anthropoids according to the degree of exposure of the fetus to postural verticality. The supporting evidence for this postulation includes a finding that in fossil hominins cerebral blood flow rate increased faster than brain volume. This testable hypothesis opens a perspective for research on fetal postural cerebral hemodynamics.

Langevin picture of subdiffusion in nonuniformly expanding medium

Anomalous diffusion phenomena have been observed in many complex physical and biological systems. One significant advance recently is the physical extension of particle's motion in a static medium to a uniformly and even nonuniformly expanding medium. The dynamic mechanism of the anomalous diffusion in the nonuniformly expanding medium has only been investigated by the approach of continuous-time random walk. To study more physical observables and to supplement the physical models of the anomalous diffusion in the expanding mediums, we characterize the nonuniformly expanding medium with a spatiotemporal dependent scale factor a(x,t) and build the Langevin picture describing the particle's motion in the nonuniformly expanding medium. Besides the existing comoving and physical coordinates, by introducing a new coordinate and assuming that a(x,t) is separable at a long-time limit, we build the relation between the nonuniformly expanding medium and the uniformly expanding one and further obtain the moments of the comoving and physical coordinates. Different forms of the scale factor a(x,t) are considered to uncover the combined effects of the particle's intrinsic diffusion and the nonuniform expansion of medium. The theoretical analyses and simulations provide the foundation for studying more anomalous diffusion phenomena in the expanding mediums.

A multiscale computational model of arterial growth and remodeling including Notch signaling

Blood vessels grow and remodel in response to mechanical stimuli. Many computational models capture this process phenomenologically, by assuming stress homeostasis, but this approach cannot unravel the underlying cellular mechanisms. Mechano-sensitive Notch signaling is well-known to be key in vascular development and homeostasis. Here, we present a multiscale framework coupling a constrained mixture model, capturing the mechanics and turnover of arterial constituents, to a cell-cell signaling model, describing Notch signaling dynamics among vascular smooth muscle cells (SMCs) as influenced by mechanical stimuli. Tissue turnover was regulated by both Notch activity, informed by in vitro data, and a phenomenological contribution, accounting for mechanisms other than Notch. This novel framework predicted changes in wall thickness and arterial composition in response to hypertension similar to previous in vivo data. The simulations suggested that Notch contributes to arterial growth in hypertension mainly by promoting SMC proliferation, while other mechanisms are needed to fully capture remodeling. The results also indicated that interventions to Notch, such as external Jagged ligands, can alter both the geometry and composition of hypertensive vessels, especially in the short term. Overall, our model enables a deeper analysis of the role of Notch and Notch interventions in arterial growth and remodeling and could be adopted to investigate therapeutic strategies and optimize vascular regeneration protocols.

Understanding the form and function in Chinese bound foot from last-generation cases

Purpose: Foot adaptation in the typically developed foot is well explored. In this study, we aimed to explore the form and function of an atypical foot, the Chinese bound foot, which had a history of over a thousand years but is not practised anymore. Methods: We evaluated the foot shape and posture via a statistical shape modelling analysis, gait plantar loading distribution via gait analysis, and bone density adaptation via implementing finite element simulation and bone remodelling prediction. Results: The atypical foot with binding practice led to increased foot arch and vertically oriented calcaneus with larger size at the articulation, apart from smaller metatarsals compared with a typically developed foot. This shape change causes the tibia, which typically acts as a load transfer beam and shock absorber, to extend its function all the way through the talus to the calcaneus. This is evident in the bound foot by i) the reduced center of pressure trajectory in the medial-lateral direction, suggesting a reduced supination-pronation; ii) the increased density and stress in the talus-calcaneus articulation; and iii) the increased bone growth in the bound foot at articulation joints in the tibia, talus, and calcaneus. Conclusion: Knowledge from the last-generation bound foot cases may provide insights into the understanding of bone resorption and adaptation in response to different loading profiles.

Sonomechanobiology: Vibrational stimulation of cells and its therapeutic implications

All cells possess an innate ability to respond to a range of mechanical stimuli through their complex internal machinery. This comprises various mechanosensory elements that detect these mechanical cues and diverse cytoskeletal structures that transmit the force to different parts of the cell, where they are transcribed into complex transcriptomic and signaling events that determine their response and fate. In contrast to static (or steady) mechanostimuli primarily involving constant-force loading such as compression, tension, and shear (or forces applied at very low oscillatory frequencies (≤ 1 Hz) that essentially render their effects quasi-static), dynamic mechanostimuli comprising more complex vibrational forms (e.g., time-dependent, i.e., periodic, forcing) at higher frequencies are less well understood in comparison. We review the mechanotransductive processes associated with such acoustic forcing, typically at ultrasonic frequencies (> 20 kHz), and discuss the various applications that arise from the cellular responses that are generated, particularly for regenerative therapeutics, such as exosome biogenesis, stem cell differentiation, and endothelial barrier modulation. Finally, we offer perspectives on the possible existence of a universal mechanism that is common across all forms of acoustically driven mechanostimuli that underscores the central role of the cell membrane as the key effector, and calcium as the dominant second messenger, in the mechanotransduction process.

Proliferating active matter

The fascinating patterns of collective motion created by autonomously driven particles have fuelled active-matter research for over two decades. So far, theoretical active-matter research has often focused on systems with a fixed number of particles. This constraint imposes strict limitations on what behaviours can and cannot emerge. However, a hallmark of life is the breaking of local cell number conservation by replication and death. Birth and death processes must be taken into account, for example, to predict the growth and evolution of a microbial biofilm, the expansion of a tumour, or the development from a fertilized egg into an embryo and beyond. In this Perspective, we argue that unique features emerge in these systems because proliferation represents a distinct form of activity: not only do the proliferating entities consume and dissipate energy, they also inject biomass and degrees of freedom capable of further self-proliferation, leading to myriad dynamic scenarios. Despite this complexity, a growing number of studies document common collective phenomena in various proliferating soft-matter systems. This generality leads us to propose proliferation as another direction of active-matter physics, worthy of a dedicated search for new dynamical universality classes. Conceptual challenges abound, from identifying control parameters and understanding large fluctuations and nonlinear feedback mechanisms to exploring the dynamics and limits of information flow in self-replicating systems. We believe that, by extending the rich conceptual framework developed for conventional active matter to proliferating active matter, researchers can have a profound impact on quantitative biology and reveal fascinating emergent physics along the way.

Recent Methods for Modifying Mechanical Properties of Tissue-Engineered Scaffolds for Clinical Applications

The limited regenerative capacity of the human body, in conjunction with a shortage of healthy autologous tissue, has created an urgent need for alternative grafting materials. A potential solution is a tissue-engineered graft, a construct which supports and integrates with host tissue. One of the key challenges in fabricating a tissue-engineered graft is achieving mechanical compatibility with the graft site; a disparity in these properties can shape the behaviour of the surrounding native tissue, contributing to the likelihood of graft failure. The purpose of this review is to examine the means by which researchers have altered the mechanical properties of tissue-engineered constructs via hybrid material usage, multi-layer scaffold designs, and surface modifications. A subset of these studies which has investigated the function of their constructs in vivo is also presented, followed by an examination of various tissue-engineered designs which have been clinically translated.

A chemo-mechano-biological modeling framework for cartilage evolving in health, disease, injury, and treatment

BACKGROUND AND OBJECTIVE

Osteoarthritis (OA) is a pervasive and debilitating disease, wherein degeneration of cartilage features prominently. Despite extensive research, we do not yet understand the cause or progression of OA. Studies show biochemical, mechanical, and biological factors affect cartilage health. Mechanical loads influence synthesis of biochemical constituents which build and/or break down cartilage, and which in turn affect mechanical loads. OA-associated biochemical profiles activate cellular activity that disrupts homeostasis. To understand the complex interplay among mechanical stimuli, biochemical signaling, and cartilage function requires integrating vast research on experimental mechanics and mechanobiology-a task approachable only with computational models. At present, mechanical models of cartilage generally lack chemo-biological effects, and biochemical models lack coupled mechanics, let alone interactions over time.

METHODS

We establish a first-of-its kind virtual cartilage: a modeling framework that considers time-dependent, chemo-mechano-biologically induced turnover of key constituents resulting from biochemical, mechanical, and/or biological activity. We include the "minimally essential" yet complex chemical and mechanobiological mechanisms. Our 3-D framework integrates a constitutive model for the mechanics of cartilage with a novel model of homeostatic adaptation by chondrocytes to pathological mechanical stimuli, and a new application of anisotropic growth (loss) to simulate degradation clinically observed as cartilage thinning.

RESULTS

Using a single set of representative parameters, our simulations of immobilizing and overloading successfully captured loss of cartilage quantified experimentally. Simulations of immobilizing, overloading, and injuring cartilage predicted dose-dependent recovery of cartilage when treated with suramin, a proposed therapeutic for OA. The modeling framework prompted us to add growth factors to the suramin treatment, which predicted even better recovery.

CONCLUSIONS

Our flexible framework is a first step toward computational investigations of how cartilage and chondrocytes mechanically and biochemically evolve in degeneration of OA and respond to pharmacological therapies. Our framework will enable future studies to link physical activity and resulting mechanical stimuli to progression of OA and loss of cartilage function, facilitating new fundamental understanding of the complex progression of OA and elucidating new perspectives on causes, treatments, and possible preventions.

Langevin picture of anomalous diffusion processes in expanding medium

The expanding medium is very common in many different fields, such as biology and cosmology. It brings a nonnegligible influence on particle's diffusion, which is quite different from the effect of an external force field. The dynamic mechanism of a particle's motion in an expanding medium has only been investigated in the framework of a continuous-time random walk. To focus on more diffusion processes and physical observables, we build the Langevin picture of anomalous diffusion in an expanding medium, and conduct detailed analyses in the framework of the Langevin equation. With the help of a subordinator, both subdiffusion process and superdiffusion process in the expanding medium are discussed. We find that the expanding medium with different changing rate (exponential form and power-law form) leads to quite different diffusion phenomena. The particle's intrinsic diffusion behavior also plays an important role. Our detailed theoretical analyses and simulations present a panoramic view of investigating anomalous diffusion in an expanding medium under the framework of the Langevin equation.

Mechanochemistry of collagen

The collagen molecular family is the result of nearly one billion years of evolution. It is a unique family of proteins, the majority of which provide general mechanical support to biological tissues. Fibril forming collagens are the most abundant collagens in vertebrate animals and are generally found in positions that resist tensile loading. In animals, cells produce fibril-forming collagen molecules that self-assemble into larger structures known as collagen fibrils. Collagen fibrils are the fundamental, continuous, load-bearing elements in connective tissues, but are often further aggregated into larger load-bearing structures, fascicles in tendon, lamellae in cornea and in intervertebral disk. We know that failure to form fibrillar collagen is embryonic lethal, and excessive collagen formation/growth (fibrosis) or uncontrolled enzymatic remodeling (type II collagen: osteoarthritis) is pathological. Collagen is thus critical to vertebrate viability and instrumental in maintaining efficient mechanical structures. However, despite decades of research, our understanding of collagen matrix formation is not complete, and we know still less about the detailed mechanisms that drive collagen remodeling, growth, and pathology. In this perspective, we examine the known role of mechanical force on the formation and development of collagenous structure. We then discuss a mechanochemical mechanism that has the potential to unify our understanding of collagenous tissue assembly dynamics, which preferentially deposits and grows collagen fibrils directly in the path of mechanical force, where the energetics should be dissuasive and where collagen fibrils are most required. We term this mechanism: Mechanochemical force-structure causality. STATEMENT OF SIGNIFICANCE: Our mechanochemical-force structure causality postulate suggests that collagen molecules are components of mechanochemically-sensitive and dynamically-responsive fibrils. Collagen molecules assemble preferentially in the path of applied strain, can be grown in place by mechanical extension, and are retained in the path of force through strain-stabilization. The mechanisms that drive this behavior operate at the level of the molecules themselves and are encoded into the structure of the biomaterial. The concept might change our understanding of structure formation, enhance our ability to treat injuries, and accelerate the development of therapeutics to prevent pathologies such as fibrosis. We suggest that collagen is a mechanochemically responsive dynamic element designed to provide a substantial "material assist" in the construction of adaptive carriers of mechanical signals.

Computer Model-Driven Design in Cardiovascular Regenerative Medicine

Continuing advances in genomics, molecular and cellular mechanobiology and immunobiology, including transcriptomics and proteomics, and biomechanics increasingly reveal the complexity underlying native tissue and organ structure and function. Identifying methods to repair, regenerate, or replace vital tissues and organs remains one of the greatest challenges of modern biomedical engineering, one that deserves our very best effort. Notwithstanding the continuing need for improving standard methods of investigation, including cell, organoid, and tissue culture, biomaterials development and fabrication, animal models, and clinical research, it is increasingly evident that modern computational methods should play increasingly greater roles in advancing the basic science, bioengineering, and clinical application of regenerative medicine. This brief review focuses on the development and application of computational models of tissue and organ mechanobiology and mechanics for purposes of designing tissue engineered constructs and understanding their development in vitro and in situ. Although the basic approaches are general, for illustrative purposes we describe two recent examples from cardiovascular medicine-tissue engineered heart valves (TEHVs) and tissue engineered vascular grafts (TEVGs)-to highlight current methods of approach as well as continuing needs.

ACL microtrauma: healing through nutrition, modified sports training, and increased recovery time

PURPOSE

Sports injuries among youth and adolescent athletes are a growing concern, particularly at the knee. Based on our current understanding of microtrauma and anterior cruciate ligament (ACL) healing characteristics, this clinical commentary describes a comprehensive plan to better manage ACL microtrauma and mitigate the likelihood of progression to a non-contact macrotraumatic ACL rupture.

METHODS

Medical literature related to non-contact ACL injuries among youth and adolescent athletes, collagen and ACL extracellular matrix metabolism, ACL microtrauma and sudden failure, and concerns related to current sports training were reviewed and synthesized into a comprehensive intervention plan.

RESULTS

With consideration for biopsychosocial model health factors, proper nutrition and modified sports training with increased recovery time, a comprehensive primary ACL injury prevention plan is described for the purpose of better managing ACL microtrauma, thereby reducing the incidence of non-contact macrotraumatic ACL rupture among youth and adolescent athletes.

CONCLUSION

Preventing non-contact ACL injuries may require greater consideration for reducing accumulated ACL microtrauma. Proper nutrition including glycine-rich collagen peptides, or gelatin-vitamin C supplementation in combination with healthy sleep, and adjusted sports training periodization with increased recovery time may improve ACL extracellular matrix collagen deposition homeostasis, decreasing sudden non-contact ACL rupture incidence likelihood in youth and adolescent athletes. Successful implementation will require compliance from athletes, parents, coaches, the sports medicine healthcare team, and event organizers. Studies are needed to confirm the efficacy of these concepts.

LEVEL OF EVIDENCE

V.

Three-dimensional chiral morphodynamics of chemomechanical active shells

Morphogenesis of active shells such as cells is a fundamental chemomechanical process that often exhibits three-dimensional (3D) large deformations and chemical pattern dynamics simultaneously. Here, we establish a chemomechanical active shell theory accounting for mechanical feedback and biochemical regulation to investigate the symmetry-breaking and 3D chiral morphodynamics emerging in the cell cortex. The active bending and stretching of the elastic shells are regulated by biochemical signals like actomyosin and RhoA, which, in turn, exert mechanical feedback on the biochemical events via deformation-dependent diffusion and inhibition. We show that active deformations can trigger chemomechanical bifurcations, yielding pulse spiral waves and global oscillations, which, with increasing mechanical feedback, give way to traveling or standing waves subsequently. Mechanical feedback is also found to contribute to stabilizing the polarity of emerging patterns, thus ensuring robust morphogenesis. Our results reproduce and unravel the experimentally observed solitary and multiple spiral patterns, which initiate asymmetric cleavage in Xenopus and starfish embryogenesis. This study underscores the crucial roles of mechanical feedback in cell development and also suggests a chemomechanical framework allowing for 3D large deformation and chemical signaling to explore complex morphogenesis in living shell-like structures.

Minimizing the Elastic Energy of Growing Leaves by Conformal Mapping

During morphogenesis, the shape of living species results from growth, stress relaxation, and remodeling. When the growth does not generate any stress, the body shape only reflects the growth density. In two dimensions, we show that stress free configurations are simply determined by the time evolution of a conformal mapping which concerns not only the boundary but also the displacement field during an arbitrary period of time inside the sample. Fresh planar leaves are good examples for our study: they have no elastic stress, almost no weight, and their shape can be easily represented by holomorphic functions. The growth factor, isotropic or anisotropic, is related to the metrics between the initial and current conformal maps. By adjusting the mathematical shape function, main characteristics such as tips (convex or concave or sharp-pointed), undulating borders, and veins can be mathematically recovered, which are in good agreement with observations. It is worth mentioning that this flexible method allows us to study complex morphologies of growing leaves such as the fenestration process in Monstera deliciosa, and can also shed light on many other 2D biological patterns.

Microvascular remodeling following skin injury

OBJECTIVE

The aim of this study was to describe possible remodeling (i.e., dilatation and elongation) of papillary capillaries induced by increased oxygen demand for the repair process following a skin wound.

METHODS

Computer-assisted video microscopy was used to examine 10 healthy volunteers before (baseline) and after (≈1 h and ≈24 h) an incision (5 mm long and 1 mm deep) on the forearm, 0-1 mm and 30 mm (control site) from the incision. We defined categories from 0 (low) to 3 (high) to grade dilatation and elongation of the nutritive papillary capillaries, as well as the visibility of the superficial vascular plexus. Approximately 10 000 capillaries from 200 films were scored.

RESULTS

The nutritive papillary capillaries were dilated and elongated (p < 0.01) after ≈24 h; that is, elongation (score 1.9 ± 0.9) vs baseline (score 0.9 ± 0.6), p < 0.01 and dilatation (score 2.2 ± 0.7) vs baseline (score 0.3 ± 0.3), p < 0.01. Superficial plexus visibility increased (p < 0.01) after ≈1 h (score 2.0 ± 0.7) and ≈24 h (score 2.7 ± 0.3) vs baseline (score 0.8 ± 0.4).

CONCLUSION

The superficial vascular skin plexus showed enhanced visibility already ≈1 h after the skin trauma. Morphological remodeling in the nutritive papillary capillaries-dilatation and elongation after ≈24 h-facilitate increased O₂  supply.

Computational analysis of the role of mechanosensitive Notch signaling in arterial adaptation to hypertension

Arteries grow and remodel in response to mechanical stimuli. Hypertension, for example, results in arterial wall thickening. Cell-cell Notch signaling between vascular smooth muscle cells (VSMCs) is known to be involved in this process, but the underlying mechanisms are still unclear. Here, we investigated whether Notch mechanosensitivity to strain may regulate arterial thickening in hypertension. We developed a multiscale computational framework by coupling a finite element model of arterial mechanics, including residual stress, to an agent-based model of mechanosensitive Notch signaling, to predict VSMC phenotypes as an indicator of growth and remodeling. Our simulations revealed that the sensitivity of Notch to strain at mean blood pressure may be a key mediator of arterial thickening in hypertensive arteries. Further simulations showed that loss of residual stress can have synergistic effects with hypertension, and that changes in the expression of Notch receptors, but not Jagged ligands, may be used to control arterial growth and remodeling and to intensify or counteract hypertensive thickening. Overall, we identify Notch mechanosensitivity as a potential mediator of vascular adaptation, and we present a computational framework that can facilitate the testing of new therapeutic and regenerative strategies.

Notch signaling regulates strain-mediated phenotypic switching of vascular smooth muscle cells

Mechanical stimuli experienced by vascular smooth muscle cells (VSMCs) and mechanosensitive Notch signaling are important regulators of vascular growth and remodeling. However, the interplay between mechanical cues and Notch signaling, and its contribution to regulate the VSMC phenotype are still unclear. Here, we investigated the role of Notch signaling in regulating strain-mediated changes in VSMC phenotype. Synthetic and contractile VSMCs were cyclically stretched for 48 h to determine the temporal changes in phenotypic features. Different magnitudes of strain were applied to investigate its effect on Notch mechanosensitivity and the phenotypic regulation of VSMCs. In addition, Notch signaling was inhibited via DAPT treatment and activated with immobilized Jagged1 ligands to understand the role of Notch on strain-mediated phenotypic changes of VSMCs. Our data demonstrate that cyclic strain induces a decrease in Notch signaling along with a loss of VSMC contractile features. Accordingly, the activation of Notch signaling during cyclic stretching partially rescued the contractile features of VSMCs. These findings demonstrate that Notch signaling has an important role in regulating strain-mediated phenotypic switching of VSMCs.

Instabilities and Geometry of Growing Tissues

We present a covariant continuum formulation of a generalized two-dimensional vertexlike model of epithelial tissues which describes tissues with different underlying geometries, and allows for an analytical macroscopic description. Using a geometrical approach and out-of-equilibrium statistical mechanics, we calculate both mechanical and dynamical instabilities of a tissue, and their dependences on various variables, including activity, and cell-shape heterogeneity (disorder). We show how both plastic cellular rearrangements and the tissue elastic response depend on the existence of mechanical residual stresses at the cellular level. Even freely growing tissues may exhibit a growth instability depending on the intrinsic proliferation rate. Our main result is an explicit calculation of the cell pressure in a homeostatic state of a confined growing tissue. We show that the homeostatic pressure can be negative and depends on the existence of mechanical residual stresses. This geometric model allows us to sort out elastic and plastic effects in a growing, flowing, tissue.

Multiscale simulations suggest a protective role of neo-adventitia in abdominal aortic aneurysms

Abdominal aortic aneurysms (AAAs) are a dangerous cardiovascular disease, the pathogenesis of which is not yet fully understood. In the present work a recent mechanopathological theory, which correlates AAA progression with microstructural and mechanical alterations in the tissue, is investigated using multiscale models. The goal is to combine these changes, within the framework of mechanobiology, with possible mechanical cues that are sensed by vascular cells along the AAA pathogenesis. Particular attention is paid to the formation of a 'neo-adventitia' on the abluminal side of the aortic wall, which is characterized by a highly random (isotropic) distribution of collagen fibers. Macro- and micro-scale results suggest that the formation of an AAA, as expected, perturbs the micromechanical state of the aortic tissue and triggers a growth and remodeling (G&R) reaction by mechanosensing cells such as fibroblasts. This G&R then leads to the formation of a thick neo-adventitia that appears to bring the micromechanical state of the tissue closer to the original homeostatic level. In this context, this new layer could act like a protective sheath, similar to the tunica adventitia in healthy aortas. This potential 'attempt at healing' by vascular cells would have important implications on the stability of the AAA wall and thus on the risk of rupture. STATEMENT OF SIGNIFICANCE: Current clinical criteria for risk assessment in AAAs are still empirical, as the causes and mechanisms of the disease are not yet fully understood. The strength of the arterial tissue is closely related to its microstructure, which in turn is remodeled by mechanosensing cells in the course of the disease. In this study, multiscale simulations show a possible connection between mechanical cues at the microscopic level and collagen G&R in AAA tissue. It should be emphasized that these micromechanical cues cannot be visualized in vivo. Therefore, the results presented here will help to advance our current understanding of the disease and motivate future experimental studies, with important implications for AAA risk assessment.

Stress-free morphing by means of compatible distortions

We study the morphing of three-dimensional objects within the framework of nonlinear elasticity with large distortions. A distortion field induces a target metric, and the configuration which is effectively realized by a material body is the one that minimizes the distance, measured through the elastic energy, between the target metric and the actual one. Morphing through distortions might have a paramount feature: the resulting configurations might be stress-free; if this is the case, the distortions field is called compatible. We maintain that the morphing through compatible distortions is a key strategy exploited by many soft biological materials, which can exhibit very large shape-change in response to distortions controlled by stimuli such as chemicals or temperature changes, while keeping their stress state almost null. Thus, the study of compatible distortions, and of the related shape-changes, is quite important. Here, we show a blueprint for stress-free morphing based on the notions of metric tensor and of Riemann curvature which can be used to design large morphing of three-dimensional objects.

In vivo development of tissue engineered vascular grafts: a fluid-solid-growth model

Methods of tissue engineering continue to advance, and multiple clinical trials are underway evaluating tissue engineered vascular grafts (TEVGs). Whereas initial concerns focused on suture retention and burst pressure, there is now a pressing need to design grafts to have optimal performance, including an ability to grow and remodel in response to changing hemodynamic loads. Toward this end, there is similarly a need for computational methods that can describe and predict the evolution of TEVG geometry, composition, and material properties while accounting for changes in hemodynamics. Although the ultimate goal is a fluid-solid-growth (FSG) model incorporating fully 3D growth and remodeling and 3D hemodynamics, lower fidelity models having high computational efficiency promise to play important roles, especially in the design of candidate grafts. We introduce here an efficient FSG model of in vivo development of a TEVG based on two simplifying concepts: mechanobiologically equilibrated growth and remodeling of the graft and an embedded control volume analysis of the hemodynamics. Illustrative simulations for a model Fontan conduit reveal the utility of this approach, which promises to be particularly useful in initial design considerations involving formal methods of optimization which otherwise add considerably to the computational expense.

Symmetry breaking and effects of nutrient walkway in time-dependent bone remodeling incorporating poroelasticity

Bone is an extraordinary biological material that continuously adapts its hierarchical microstructure to respond to static and dynamic loads for offering optimal mechanical features, in terms of stiffness and toughness, across different scales, from the sub-microscopic constituents within osteons—where the cyclic activity of osteoblasts, osteoclasts, and osteocytes redesigns shape and percentage of mineral crystals and collagen fibers—up to the macroscopic level, with growth and remodeling processes that modify the architecture of both compact and porous bone districts. Despite the intrinsic complexity of the bone mechanobiology, involving coupling phenomena of micro-damage, nutrients supply driven by fluid flowing throughout hierarchical networks, and cells turnover, successful models and numerical algorithms have been presented in the literature to predict, at the macroscale, how bone remodels under mechanical stimuli, a fundamental issue in many medical applications such as optimization of femur prostheses and diagnosis of the risk fracture. Within this framework, one of the most classical strategies employed in the studies is the so-called Stanford’s law, which allows uploading the effect of the time-dependent load-induced stress stimulus into a biomechanical model to guess the bone structure evolution. In the present work, we generalize this approach by introducing the bone poroelasticity, thus incorporating in the model the role of the fluid content that, by driving nutrients and contributing to the removal of wastes of bone tissue cells, synergistically interacts with the classical stress fields to change homeostasis states, local saturation conditions, and reorients the bone density rate, in this way affecting growth and remodeling. Through two paradigmatic example applications, i.e. a cylindrical slice with internal prescribed displacements idealizing a tract of femoral diaphysis pushed out by the pressure exerted by a femur prosthesis and a bone element in a form of a bent beam, it is highlighted that the present model is capable to catch more realistically both the transition between spongy and cortical regions and the expected non-symmetrical evolution of bone tissue density in the medium–long term, unpredictable with the standard approach. A real study case of a femur is also considered at the end in order to show the effectiveness of the proposed remodeling algorithm.

A Continuum Mathematical Model of Substrate-Mediated Tissue Growth

We consider a continuum mathematical model of biological tissue formation inspired by recent experiments describing thin tissue growth in 3D-printed bioscaffolds. The continuum model, which we call the substrate model, involves a partial differential equation describing the density of tissue, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{u}}(\hat{{\mathbf {x}}},{\hat{t}})$$\end{document}u^(x^,t^) that is coupled to the concentration of an immobile extracellular substrate, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{s}}(\hat{{\mathbf {x}}},{\hat{t}})$$\end{document}s^(x^,t^). Cell migration is modelled with a nonlinear diffusion term, where the diffusive flux is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{s}}$$\end{document}s^, while a logistic growth term models cell proliferation. The extracellular substrate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{s}}$$\end{document}s^ is produced by cells and undergoes linear decay. Preliminary numerical simulations show that this mathematical model is able to recapitulate key features of recent tissue growth experiments, including the formation of sharp fronts. To provide a deeper understanding of the model we analyse travelling wave solutions of the substrate model, showing that the model supports both sharp-fronted travelling wave solutions that move with a minimum wave speed, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c = c_{\mathrm{min}}$$\end{document}c=cmin, as well as smooth-fronted travelling wave solutions that move with a faster travelling wave speed, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$c > c_{\mathrm{min}}$$\end{document}c>cmin. We provide a geometric interpretation that explains the difference between smooth and sharp-fronted travelling wave solutions that is based on a slow manifold reduction of the desingularised three-dimensional phase space. In addition, we also develop and test a series of useful approximations that describe the shape of the travelling wave solutions in various limits. These approximations apply to both the sharp-fronted and smooth-fronted travelling wave solutions. Software to implement all calculations is available at GitHub.

A macroscopic approach for stress-driven anisotropic growth in bioengineered soft tissues

The simulation of growth processes within soft biological tissues is of utmost importance for many applications in the medical sector. Within this contribution, we propose a new macroscopic approach for modelling stress-driven volumetric growth occurring in soft tissues. Instead of using the standard approach of a-priori defining the structure of the growth tensor, we postulate the existence of a general growth potential. Such a potential describes all eligible homeostatic stress states that can ultimately be reached as a result of the growth process. Making use of well-established methods from visco-plasticity, the evolution of the growth-related right Cauchy–Green tensor is subsequently defined as a time-dependent associative evolution law with respect to the introduced potential. This approach naturally leads to a formulation that is able to cover both, isotropic and anisotropic growth-related changes in geometry. It furthermore allows the model to flexibly adapt to changing boundary and loading conditions. Besides the theoretical development, we also describe the algorithmic implementation and furthermore compare the newly derived model with a standard formulation of isotropic growth.

Mathematical models of neuronal growth

The establishment of a functioning neuronal network is a crucial step in neural development. During this process, neurons extend neurites—axons and dendrites—to meet other neurons and interconnect. Therefore, these neurites need to migrate, grow, branch and find the correct path to their target by processing sensory cues from their environment. These processes rely on many coupled biophysical effects including elasticity, viscosity, growth, active forces, chemical signaling, adhesion and cellular transport. Mathematical models offer a direct way to test hypotheses and understand the underlying mechanisms responsible for neuron development. Here, we critically review the main models of neurite growth and morphogenesis from a mathematical viewpoint. We present different models for growth, guidance and morphogenesis, with a particular emphasis on mechanics and mechanisms, and on simple mathematical models that can be partially treated analytically.

Growth and remodeling in the pulmonary autograft: Computational evaluation using kinematic growth models and constrained mixture theory

Computational investigations of how soft tissues grow and remodel are gaining more and more interest and several growth and remodeling theories have been developed. Roughly, two main groups of theories for soft tissues can be distinguished: kinematic-based growth theory and theories based on constrained mixture theory. Our goal was to apply these two theories on the same experimental data. Within the experiment, a pulmonary artery was exposed to systemic conditions. The change in diameter was followed-up over time. A mechanical and microstructural analysis of native pulmonary artery and pulmonary autograft was conducted. Whereas the kinematic-based growth theory is able to accurately capture the growth of the tissue, it does not account for the mechanobiological processes causing this growth. The constrained mixture theory takes into account the mechanobiological processes including removal, deposition and adaptation of all structural constituents, allowing us to simulate a changing microstructure and mechanical behavior.

Bayesian calibration of a computational model of tissue expansion based on a porcine animal model

Tissue expansion is a technique used clinically to grow skin in situ to correct large defects. Despite its enormous potential, lack of fundamental knowledge of skin adaptation to mechanical cues, and lack of predictive computational models limit the broader adoption and efficacy of tissue expansion. In our previous work, we introduced a finite element model of tissue expansion that predicted key patterns of strain and growth which were then confirmed by our porcine animal model. Here we use the data from a new set of experiments to calibrate the computational model within a Bayesian framework. Four 10×10cm2 patches were tattooed in the dorsal skin of four 12 weeks-old minipigs and a total of six patches underwent successful tissue expander placement and inflation to 60cc for expansion times ranging from 1 h to 7 days. Six patches that did not have expanders implanted served as controls for the analysis. We find that growth can be explained based on the elastic deformation. The predicted area growth rate is k∈[0.02,0.08] [h-1]. Growth is anisotropic and reflects the anisotropic mechanical behavior of porcine dorsal skin. The rostral-caudal axis shows greater deformation than the transverse axis, and the time scale of growth in the rostral-caudal direction is given by rate parameters k1∈[0.04,0.1] [h-1] compared to k2∈[0.01,0.05] [h-1] in the transverse direction. Moreover, the calibration results underscore the high variability in biological systems, and the need to create probabilistic computational models to predict tissue adaptation in realistic settings. STATEMENT OF SIGNIFICANCE: Tissue expansion is a widely used technique in reconstructive surgery because it triggers growth of skin for the correction of large skin lesions and for breast reconstruction after mastectomy. Despite of its potential, complications and undesired outcomes persist due to our incomplete understanding of skin mechanobiology. Here we quantify the deformation and growth fields induced by an expander over 7 days in a porcine animal model and use these data to calibrate a computational model of skin growth using finite element simulations and a Bayesian framework. The calibrated model is a leap forward in our understanding skin growth, we now have quantitative understanding of this process: area growth is anisotropic and it is proportional to stretch with a characteristic rate constant of k∈[0.02,0.08] [h-1].

A Direct Comparison of Node and Element-Based Finite Element Modeling Approaches to Study Tissue Growth

Finite element analysis is a useful tool to model growth of biological tissues and predict how growth can be impacted by stimuli. Previous work has simulated growth using node-based or element-based approaches, and this implementation choice may influence predicted growth, irrespective of the applied growth model. This study directly compared node-based and element-based approaches to understand the isolated impact of implementation method on growth predictions by simulating growth of a bone rudiment geometry, and determined what conditions produce similar results between the approaches. We used a previously reported node-based approach implemented via thermal expansion and an element-based approach implemented via osmotic swelling, and we derived a mathematical relationship to relate the growth resulting from these approaches. We found that material properties (modulus) affected growth in the element-based approach, with growth completely restricted for high modulus values relative to the growth stimulus, and no restriction for low modulus values. The node-based approach was unaffected by modulus. Node- and element-based approaches matched marginally better when the conversion coefficient to relate the approaches was optimized based on the results of initial simulations, rather than using the theoretically predicted conversion coefficient (median difference in node position 0.042 cm versus 0.052 cm, respectively). In summary, we illustrate here the importance of the choice of implementation approach for modeling growth, provide a framework for converting models between implementation approaches, and highlight important considerations for comparing results in prior work and developing new models of tissue growth.

The physical basis of mollusk shell chiral coiling

Snails are model organisms for studying the genetic, molecular, and developmental bases of left-right asymmetry in Bilateria. However, the development of their typical helicospiral shell, present for the last 540 million years in environments as different as the abyss or our gardens, remains poorly understood. Conversely, ammonites typically have a bilaterally symmetric, planispiraly coiled shell, with only 1% of 3,000 genera displaying either a helicospiral or a meandering asymmetric shell. A comparative analysis suggests that the development of chiral shells in these mollusks is different and that, unlike snails, ammonites with asymmetric shells probably had a bilaterally symmetric body diagnostic of cephalopods. We propose a mathematical model for the growth of shells, taking into account the physical interaction during development between the soft mollusk body and its hard shell. Our model shows that a growth mismatch between the secreted shell tube and a bilaterally symmetric body in ammonites can generate mechanical forces that are balanced by a twist of the body, breaking shell symmetry. In gastropods, where a twist is intrinsic to the body, the same model predicts that helicospiral shells are the most likely shell forms. Our model explains a large diversity of forms and shows that, although molluscan shells are incrementally secreted at their opening, the path followed by the shell edge and the resulting form are partly governed by the mechanics of the body inside the shell, a perspective that explains many aspects of their development and evolution.

Local and global growth and remodeling in calcific aortic valve disease and aging

Aging and calcific aortic valve disease (CAVD) are the main factors leading to aortic stenosis. Both processes are accompanied by growth and remodeling pathways that play a crucial role in aortic valve pathophysiology. Herein, a computational growth and remodeling (G&R) framework was developed to investigate the effects of aging and calcification on aortic valve dynamics. Particularly, an algorithm was developed to couple the global growth and stiffening of the aortic valve due to aging and the local growth and stiffening due to calcification with the aortic valve transient dynamics. The aortic valve dynamics during baseline were validated with available data in the literature. Subsequently, the changes in aortic valve dynamic patterns during aging and CAVD progression were studied. The results revealed the patterns in geometric orifice area reduction and an increase in the valve stress during local and global growth and remodeling of the aortic valve. The proposed algorithm provides a framework to couple mechanobiology models of disease growth with tissue-scale transient structural mechanics models to study the biomechanical changes during cardiovascular disease growth and aging.

Physically motivated models of polymer networks with dynamic cross-links: comparative study and future outlook

Polymer networks consisting of a mixture of chemical and physical cross-links are known to exhibit complex time-dependent behaviour due to the kinetics of bond association and dissociation. In this article, we highlight and compare two recent physically based constitutive models that describe the nonlinear viscoelastic behaviour of such transient networks. These two models are developed independently by two groups of researchers using different mathematical formulations. Here, we show that this difference can be attributed to different viewpoints: Lagrangian versus Eulerian. We establish the equivalence of the two models under the special situation where chains obey Gaussian statistics and steady-state bond dynamics. We provide experimental data demonstrating that both models can accurately predict the time-dependent uniaxial behaviour of a poly(vinylalcohol) dual cross-link hydrogel. We review the advantages and disadvantages of both approaches in applications and close by discussing a list of open challenges and questions regarding the mathematical modelling of soft, viscoelastic networks.

Influence of the growth gradient on surface wrinkling and pattern transition in growing tubular tissues

Growth-induced pattern formations in curved film-substrate structures have attracted extensive attention recently. In most existing literature, the growth tensor is assumed to be homogeneous or piecewise homogeneous. In this paper, we aim at clarifying the influence of a growth gradient on pattern formation and pattern evolution in bilayered tubular tissues under plane-strain deformation. In the framework of finite elasticity, a bifurcation condition is derived for a general material model and a generic growth function. Then we suppose that both layers are composed of neo-Hookean materials. In particular, the growth function is assumed to decay linearly either from the inner surface or from the outer surface. It is found that a gradient in the growth has a weak effect on the critical state, compared with the homogeneous growth type where both layers share the same growth factor. Furthermore, a finite-element model is built to validate the theoretical model and to investigate the post-buckling behaviours. It is found that the associated pattern transition is not controlled by the growth gradient but by the ratio of the shear modulus between two layers. Different morphologies can occur when the modulus ratio is varied. The current analysis could provide useful insight into the influence of a growth gradient on surface instabilities and suggests that a homogeneous growth field may provide a good approximation on interpreting complicated morphological formations in multiple systems.

Condensation tendency and planar isotropic actin gradient induce radial alignment in confined monolayers

A monolayer of highly motile cells can establish long-range orientational order, which can be explained by hydrodynamic theory of active gels and fluids. However, it is less clear how cell shape changes and rearrangement are governed when the monolayer is in mechanical equilibrium states when cell motility diminishes. In this work, we report that rat embryonic fibroblasts (REF), when confined in circular mesoscale patterns on rigid substrates, can transition from the spindle shapes to more compact morphologies. Cells align radially only at the pattern boundary when they are in the mechanical equilibrium. This radial alignment disappears when cell contractility or cell-cell adhesion is reduced. Unlike monolayers of spindle-like cells such as NIH-3T3 fibroblasts with minimal intercellular interactions or epithelial cells like Madin-Darby canine kidney (MDCK) with strong cortical actin network, confined REF monolayers present an actin gradient with isotropic meshwork, suggesting the existence of a stiffness gradient. In addition, the REF cells tend to condense on soft substrates, a collective cell behavior we refer to as the 'condensation tendency'. This condensation tendency, together with geometrical confinement, induces tensile prestretch (i.e. an isotropic stretch that causes tissue to contract when released) to the confined monolayer. By developing a Voronoi-cell model, we demonstrate that the combined global tissue prestretch and cell stiffness differential between the inner and boundary cells can sufficiently define the cell radial alignment at the pattern boundary.

Fractal Pennes and Cattaneo-Vernotte bioheat equations from product-like fractal geometry and their implications on cells in the presence of tumour growth

In this study, the Pennes and Cattaneo-Vernotte bioheat transfer equations in the presence of fractal spatial dimensions are derived based on the product-like fractal geometry. This approach was introduced recently, by Li and Ostoja-Starzewski, in order to explore dynamical properties of anisotropic media. The theory is characterized by a modified gradient operator which depends on two parameters: R which represents the radius of the tumour and R0 which represents the radius of the spherical living tissue. Both the steady and unsteady states for each fractal bioheat equation were obtained and their implications on living cells in the presence of growth of a large tumour were analysed. Assuming a specific heating/cooling by a constant heat flux equivalent to the metabolic heat generation in the tissue, it was observed that the solutions of the fractal bioheat equations are robustly affected by fractal dimensions, the radius of the tumour growth and the dimensions of the living cell tissue. The ranges of both the fractal dimensions and temperature were obtained, analysed and compared with recent studies. This study confirms the importance of fractals in medicine.

Mechano-chemo-biological Computational Models for Arteries in Health, Disease and Healing: From Tissue Remodelling to Drug-eluting Devices

This review aims to highlight urgent priorities for the computational biomechanics community in the framework of mechano-chemo-biological models. Recent approaches, promising directions and open challenges on the computational modelling of arterial tissues in health and disease are introduced and investigated, together with in silico approaches for the analysis of drug-eluting stents that promote pharmacological-induced healing. The paper addresses a number of chemo-biological phenomena that are generally neglected in biomechanical engineering models but are most likely instrumental for the onset and the progression of arterial diseases. An interdisciplinary effort is thus encouraged for providing the tools for an effective in silico insight into medical problems. An integrated mechano-chemo-biological perspective is believed to be a fundamental missing piece for crossing the bridge between computational engineering and life sciences, and for bringing computational biomechanics into medical research and clinical practice.

Model-based data analysis of tissue growth in thin 3D printed scaffolds

Tissue growth in three-dimensional (3D) printed scaffolds enables exploration and control of cell behaviour in more biologically realistic geometries than that allowed by traditional 2D cell culture. Cell proliferation and migration in these experiments have yet to be explicitly characterised, limiting the ability of experimentalists to determine the effects of various experimental conditions, such as scaffold geometry, on cell behaviour. We consider tissue growth by osteoblastic cells in melt electro-written scaffolds that comprise thin square pores with sizes that were deliberately increased between experiments. We collect highly detailed temporal measurements of the average cell density, tissue coverage, and tissue geometry. To quantify tissue growth in terms of the underlying cell proliferation and migration processes, we introduce and calibrate a mechanistic mathematical model based on the Porous-Fisher reaction-diffusion equation. Parameter estimates and uncertainty quantification through profile likelihood analysis reveal consistency in the rate of cell proliferation and steady-state cell density between pore sizes. This analysis also serves as an important model verification tool: while the use of reaction-diffusion models in biology is widespread, the appropriateness of these models to describe tissue growth in 3D scaffolds has yet to be explored. We find that the Porous-Fisher model is able to capture features relating to the cell density and tissue coverage, but is not able to capture geometric features relating to the circularity of the tissue interface. Our analysis identifies two distinct stages of tissue growth, suggests several areas for model refinement, and provides guidance for future experimental work that explores tissue growth in 3D printed scaffolds.

Constrained Mixture Models of Soft Tissue Growth and Remodeling - Twenty Years After

Soft biological tissues compromise diverse cell types and extracellular matrix constituents, each of which can possess individual natural configurations, material properties, and rates of turnover. For this reason, mixture-based models of growth (changes in mass) and remodeling (change in microstructure) are well-suited for studying tissue adaptations, disease progression, and responses to injury or clinical intervention. Such approaches also can be used to design improved tissue engineered constructs to repair, replace, or regenerate tissues. Focusing on blood vessels as archetypes of soft tissues, this paper reviews a constrained mixture theory introduced twenty years ago and explores its usage since by contrasting simulations of diverse vascular conditions. The discussion is framed within the concept of mechanical homeostasis, with consideration of solid-fluid interactions, inflammation, and cell signaling highlighting both past accomplishments and future opportunities as we seek to understand better the evolving composition, geometry, and material behaviors of soft tissues under complex conditions.

Vascular Mechanobiology: Homeostasis, Adaptation, and Disease

Cells of the vascular wall are exquisitely sensitive to changes in their mechanical environment. In healthy vessels, mechanical forces regulate signaling and gene expression to direct the remodeling needed for the vessel wall to maintain optimal function. Major diseases of arteries involve maladaptive remodeling with compromised or lost homeostatic mechanisms. Whereas homeostasis invokes negative feedback loops at multiple scales to mediate mechanobiological stability, disease progression often occurs via positive feedback that generates mechanobiological instabilities. In this review, we focus on the cell biology, wall mechanics, and regulatory pathways associated with arterial health and how changes in these processes lead to disease. We discuss how positive feedback loops arise via biomechanical and biochemical means. We conclude that inflammation plays a central role in overriding homeostatic pathways and suggest future directions for addressing therapeutic needs.

Nutations in growing plant shoots as a morphoelastic flutter instability

Growing plant shoots exhibit spontaneous oscillations that Darwin observed, and termed 'circumnutations'. Recently, they have received renewed attention for the design and optimal actuation of bioinspired robotic devices. We discuss a possible interpretation of these spontaneous oscillations as a Hopf-type bifurcation in a growing morphoelastic rod. Using a three-dimensional model and numerical simulations, we analyse the salient features of this flutter-like phenomenon (e.g. the characteristic period of the oscillations) and their dependence on the model details (in particular, the impact of choosing different growth models) finding that, overall, these features are robust with respect to changes in the details of the growth model adopted. This article is part of the theme issue 'Topics in mathematical design of complex materials'.

Patient-Specific Computational Fluid Dynamics Reveal Localized Flow Patterns Predictive of Post-Left Ventricular Assist Device Aortic Incompetence

BACKGROUND

Progressive aortic valve disease has remained a persistent cause of concern in patients with left ventricular assist devices. Aortic incompetence (AI) is a known predictor of both mortality and readmissions in this patient population and remains a challenging clinical problem.

METHODS

Ten left ventricular assist device patients with de novo aortic regurgitation and 19 control left ventricular assist device patients were identified. Three-dimensional models of patients' aortas were created from their computed tomography scans, following which large-scale patient-specific computational fluid dynamics simulations were performed with physiologically accurate boundary conditions using the SimVascular flow solver.

RESULTS

The spatial distributions of time-averaged wall shear stress and oscillatory shear index show no significant differences in the aortic root in patients with and without AI (mean difference, 0.67 dyne/cm2 [95% CI, -0.51 to 1.85]; P=0.23). Oscillatory shear index was also not significantly different between both groups of patients (mean difference, 0.03 [95% CI, -0.07 to 0.019]; P=0.22). The localized wall shear stress on the leaflet tips was significantly higher in the AI group than the non-AI group (1.62 versus 1.35 dyne/cm2; mean difference [95% CI, 0.15-0.39]; P<0.001), whereas oscillatory shear index was not significantly different between both groups (95% CI, -0.009 to 0.001; P=0.17).

CONCLUSIONS

Computational fluid dynamics serves a unique role in studying the hemodynamic features in left ventricular assist device patients where 4-dimensional magnetic resonance imaging remains unfeasible. Contrary to the widely accepted notions of highly disturbed flow, in this study, we demonstrate that the aortic root is a region of relatively stagnant flow. We further identified localized hemodynamic features in the aortic root that challenge our understanding of how AI develops in this patient population.

Towards an understanding of the mechanoreciprocity process in adipocytes and its perturbation with aging

Adipose tissue (AT) is a complex organ, with multiple functions that are essential for maintaining metabolic health. A feature of AT is its capability to expand in response to physiological challenges, such as pregnancy and aging, and during chronic states of positive energy balance occurring throughout life. AT grows through adipogenesis and/or an increase in the size of existing adipocytes. One process that is required for healthy AT growth is the remodeling of the extracellular matrix (ECM), which is a necessary step to restore mechanical homeostasis and maintain tissue integrity and functionality. While the relationship between mechanobiology and adipogenesis is now well recognized, less is known about the role of adipocyte mechanosignaling pathways in AT growth. In this review article, we first summarize evidence linking ECM remodelling to AT expansion and how its perturbation is associated to a metabolically unhealthy phenotype. Subsequently, we highlight findings suggesting that molecules involved in the dynamic, bidirectional process (mechanoreciprocity) enabling adipocytes to sense changes in the mechanical properties of the ECM are interconnected to pathways regulating lipid metabolism. Finally, we discuss processes through which aging may influence the ability of adipocytes to appropriately respond to alterations in ECM composition.

Multiscale modeling meets machine learning: What can we learn?

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learning to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

Machine learning approaches to surrogate multifidelity Growth and Remodeling models for efficient abdominal aortic aneurysmal applications

Computational Growth and Remodeling (G&R) models have been widely used to capture the pathological development of arterial diseases and have shown promise for aiding clinical diagnosis, prognosis prediction, and staging classification. However, due to the high complexity of the arterial adaptation mechanism, high-fidelity arterial G&R simulation usually takes hours or even days, which hinders its application in clinical practice. To remedy this problem, we develop a computationally efficient arterial G&R simulation framework that comprehensively combines the physics-based G&R simulations and data-driven machine learning approaches. The proposed framework greatly enhances the computational efficiency of arterial G&R simulations, thereby enabling more time-consuming arterial applications, including personalized parameter estimation and arterial disease progression prediction. In particular, we achieve significant computational cost reduction mainly through two methods: (1) constructing a Multifidelity Surrogate (MFS) to approximate multifidelity G&R simulations by using a cokriging approach and (2) developing a novel iterative optimization algorithm for personalized parameter estimation. The proposed framework is demonstrated by estimating G&R model parameters and predicting individual aneurysm growth using follow-up CT images of Abdominal Aortic Aneurysms (AAAs) from 21 patients. Results show that the personalized parameters are satisfactorily estimated and the growth of AAAs is predicted within the clinically relevant time frame, i.e., less than 2 h, without a loss of accuracy.

Morphoelasticity of large bending deformations of cell sheets during development

Deformations of cell sheets during morphogenesis are driven by developmental processes such as cell division and cell shape changes. In morphoelastic shell theories of development, these processes appear as variations of the intrinsic geometry of a thin elastic shell. However, morphogenesis often involves large bending deformations that are outside the formal range of validity of these shell theories. Here, by asymptotic expansion of three-dimensional incompressible morphoelasticity in the limit of a thin shell, we derive a shell theory for large intrinsic bending deformations and emphasize the resulting geometric material anisotropy and the elastic role of cell constriction. Taking the invagination of the green alga Volvox as a model developmental event, we show how results for this theory differ from those for a classical shell theory that is not formally valid for these large bending deformations and reveal how these geometric effects stabilize invagination.