Computational model of damage-induced growth in soft biological tissues considering the mechanobiology of healing

On the Impact of Residual Strains in the Stress Analysis of Patient-Specific Atherosclerotic Carotid Vessels: Predictions Based on the Homogenous Stress Hypothesis

The identification of carotid atherosclerotic lesion at risk for plaque rupture, eventually resulting in cerebral embolism and stroke, is of paramount clinical importance. High stress in the fibrous plaque cap has been proposed as risk factor. However, among others, residual strains influence said stress predictions, but quantitative and qualitative implications of residual strains in this context are not well explored. We therefore propose a multiplicative kinematics-based Growth and Remodeling (G&R) framework to predict residual strains from homogenizing tissue stress and then investigate its implication on plaque stress. Carotid vessel morphology of four patients was reconstructed from clinical Computed Tomography-Angiography (CT-A) images and equipped with heterogeneous tissue constitutive properties assigned through a histology-based artificial intelligence image segmentation tool. As compared to a purely elastic analysis and depending on patient-specific morphology and tissue distributions, the incorporation of residual strains reduced the maximum wall stress by up to 30 % and resulted in a fundamentally different distribution of stress across the atherosclerotic wall. Regardless residual strains homogenized tissue stresses, the fibrous plaque cap may persistently be exposed to spots of high stress. In conclusion, the incorporation of residual strains in biomechanical studies of atherosclerotic carotids may be important for a reliable assessment of fibrous plaque cap stress.

Mechano-chemo-biological model of atherosclerosis formation based on the outside-in theory

Atherosclerosis is a disease in blood vessels that often results in plaque formation and lumen narrowing. It is an inflammatory response of the tissue caused by disruptions in the vessel wall nourishment. Blood vessels are nourished by nutrients originating from the blood of the lumen. In medium-sized and larger vessels, nutrients are additionally provided from outside through a network of capillaries called vasa vasorum. It has recently been hypothesized (Haverich in Circulation 135:205-207, 2017) that the root of atherosclerotic diseases is the malfunction of the vasa vasorum. This, so-called outside-in theory, is supported by a recently developed numerical model (Soleimani et al. in Arch Comput Methods Eng 28:4263-4282, 2021) accounting for the inflammation initiation in the adventitial layer of the blood vessel. Building on the previous findings, this work proposes an extended material model for atherosclerosis formation that is based on the outside-in theory. Beside the description of growth kinematics and nutrient diffusion, the roles of monocytes, macrophages, foam cells, smooth muscle cells and collagen are accounted for in a nonlinear continuum mechanics framework. Cells are activated due to a lack of vessel wall nourishment and proliferate, migrate, differentiate and synthesize collagen, leading to the formation of a plaque. Numerical studies show that the onset of atherosclerosis can qualitatively be reproduced and back the new theory.

Post-angioplasty remodeling of coronary arteries investigated via a chemo-mechano-biological in silico model

This work presents the application of a chemo-mechano-biological constitutive model of soft tissues for describing tissue inflammatory response to damage in collagen constituents. The material model is implemented into a nonlinear finite element formulation to follow up a coronary standard balloon angioplasty for one year. Numerical results, compared with available in vivo clinical data, show that the model reproduces the temporal dynamics of vessel remodeling associated with subintimal damage. Such dynamics are bimodular, being characterized by an early tissue resorption and lumen enlargement, followed by late tissue growth and vessel constriction. Applicability of the modeling framework in retrospective studies is demonstrated, and future extension towards prospective applications is discussed.

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.

A biochemomechanical model of collagen turnover in arterial adaptations to hemodynamic loading

The production, removal, and remodeling of fibrillar collagen is fundamental to mechanical homeostasis in arteries, including dynamic morphological and microstructural changes that occur in response to sustained changes in blood flow and pressure under physiological conditions. These dynamic processes involve complex, coupled biological, chemical, and mechanical mechanisms that are not completely understood. Nevertheless, recent simulations using constrained mixture models with phenomenologically motivated constitutive relations have proven able to predict salient features of the progression of certain vascular adaptations as well as disease processes. Collagen turnover is modeled, in part, via stress-dependent changes in collagen half-life, typically within the range of 10-70 days. By contrast, in this work we introduce a biochemomechanical approach to model the cellular synthesis of procollagen as well as its transition from an intermediate state of assembled microfibrils to mature cross-linked fibers, with mechano-regulated removal. The resulting model can simulate temporal changes in geometry, composition, and stress during early vascular adaptation (weeks to months) for modest changes in blood flow or pressure. It is shown that these simulations capture salient features from data presented in the literature from different animal models.

Arterial tissues and their inflammatory response to collagen damage: A continuum in silico model coupling nonlinear mechanics, molecular pathways, and cell behavior

Damage in soft biological tissues causes an inflammatory reaction that initiates a chain of events to repair the tissue. This work presents a continuum model and its in silico implementation that describe the cascade of mechanisms leading to tissue healing, coupling mechanical as well as chemo-biological processes. The mechanics is described by means of a Lagrangian nonlinear continuum mechanics framework and follows the homogenized constrained mixtures theory. Plastic-like damage, growth and remodeling as well as homeostasis are taken into account. The chemo-biological pathways account for two molecular and four cellular species, and are activated by damage of collagen molecules in fibers. To consider proliferation, differentiation, diffusion and chemotaxis of species, diffusion-advection-reaction equations are employed. To the best of authors' knowledge, the proposed model combines for the first time such high number of chemo-mechano-biological mechanisms in a consistent continuum biomechanical framework. The resulting set of coupled differential equations describe balance of linear momentum, evolution of kinematic variables as well as mass balance equations. They are discretized in time according to a backward Euler finite difference scheme, and in space through a finite element Galerkin discretization. The features of the model are firstly demonstrated presenting the species dynamics and highlighting the influence of damage intensities on the growth outcome. In terms of a biaxial test, the chemo-mechano-biological coupling and the model's applicability to reproduce normal as well as pathological healing are shown. A last numerical example underlines the model's applicability to complex loading scenarios and inhomogeneous damage distributions. Concluding, the present work contributes towards comprehensive in silico models in biomechanics and mechanobiology.

Computer simulation of the wound process (review of literature)

Relevance. Computer simulation is a mathematical modeling process performed on a computer that is designed to predict the behavior or results of a real or physical system. Computer simulation has a number of advantages over classical models of animal experiments: the cheapness of the method (the need to acquire and maintain animals disappears by itself), the speed of obtaining results, the absence of bioethical problems, the ability to change the conditions of the experiment, etc.

he purpose of this study is to review the methods of computer simulation of the wound process, to identify the shortcomings of the models and propose ways to solve them, as well as to select the best existing model for describing wound regeneration.

Material and methods. In the course of this work, an analysis was made of foreign and domestic literature on the problem of computer modeling of the wound process.

Results. After analyzing the relevant literature on this topic, the problem is seen precisely in the insufficiently studied process of wound regeneration, since many different cells, cytokines, growth factors, enzymes, fibrillar proteins, etc. take part in it. The models that currently exist describe wound regeneration only in an extremely generalized way, which does not allow us to apply them in clinical situations. Analyzing literature sources, we came to the conclusion that both numerical approaches, both cellular-biochemical (the first type of models) and phenomenological (the second type) are applicable in the case of wound modeling and can be used very successfully. The problem is that on the basis of one approach it is impossible to display a complete picture of wound healing, in this way it is possible to predict only individual regeneration parameters necessary for certain purposes due to the complexity and versatility of this typical pathophysiological process.

Conclusion. Computer modeling of wounds is still a controversial and complex topic. Existing models are not intended to describe all the processes occurring in a healing wound. It is much more productive to describe the various phenomena during healing separately. This is due to the fact that many elements are involved in the regeneration of the skin, which are almost impossible to take into account in full. The available models are of exclusively scientific value, consisting in attempts to understand all complex processes and interactions. Practical application is difficult, since existing models require specific input data that require highly specialized equipment. If we abstract from all this, then the best existing model of the first type is the model of the authors Yangyang Wang, Christian F. Guerrero-Juarez, Yuchi Qiu and co-authors, in addition to it, any of the described phenomenological models will do. 

An intricate interplay between stent drug dose and release rate dictates arterial restenosis

Since the introduction of percutaneous coronary intervention (PCI) for the treatment of obstructive coronary artery disease (CAD), patient outcomes have progressively improved. Drug eluting stents (DES) that employ anti-proliferative drugs to limit excess tissue growth following stent deployment have proved revolutionary. However, restenosis and a need for repeat revascularisation still occurs after DES use. Over the last few years, computational models have emerged that detail restenosis following the deployment of a bare metal stent (BMS), focusing primarily on contributions from mechanics and fluid dynamics. However, none of the existing models adequately account for spatiotemporal delivery of drug and the influence of this on the cellular processes that drive restenosis. In an attempt to fill this void, a novel continuum restenosis model coupled with spatiotemporal drug delivery is presented. Our results indicate that the severity and time-course of restenosis is critically dependent on the drug delivery strategy. Specifically, we uncover an intricate interplay between initial drug loading, drug release rate and restenosis, indicating that it is not sufficient to simply ramp-up the drug dose or prolong the time course of drug release to improve stent efficacy. Our model also shows that the level of stent over-expansion and stent design features, such as inter-strut spacing and strut thickness, influence restenosis development, in agreement with trends observed in experimental and clinical studies. Moreover, other critical aspects of the model which dictate restenosis, including the drug binding site density are investigated, where comparisons are made between approaches which assume this to be either constant or proportional to the number of smooth muscle cells (SMCs). Taken together, our results highlight the necessity of incorporating these aspects of drug delivery in the pursuit of optimal DES design.