A Comparison of Fully-Coupled 3D In-Stent Restenosis Simulations to In-vivo Data

Impact of geometric and hemodynamic changes on a mechanobiological model of atherosclerosis

BACKGROUND AND OBJECTIVE

In this work, the analysis of the importance of hemodynamic updates on a mechanobiological model of atheroma plaque formation is proposed.

METHODS

For that, we use an idealized and axisymmetric model of carotid artery. In addition, the behavior of endothelial cells depending on hemodynamical changes is analyzed too. A total of three computational simulations are carried out and their results are compared: an uncoupled model and two models that consider the opposite behavior of endothelial cells caused by hemodynamic changes. The model considers transient blood flow using the Navier-Stokes equation. Plasma flow across the endothelium is determined with Darcy's law and the Kedem-Katchalsky equations, considering the three-pore model, which is also employed for the flow of substances across the endothelium. The behavior of the considered substances in the arterial wall is modeled with convection-diffusion-reaction equations, and the arterial wall is modeled as a hyperelastic Yeoh's material.

RESULTS

Significant variations are noted in both the morphology and stenosis ratio of the plaques when comparing the uncoupled model to the two models incorporating updates for geometry and hemodynamic stimuli. Besides, the phenomenon of double-stenosis is naturally reproduced in the models that consider both geometric and hemodynamical changes due to plaque growth, whereas it cannot be predicted in the uncoupled model.

CONCLUSIONS

The findings indicate that integrating the plaque growth model with geometric and hemodynamic settings is essential in determining the ultimate shape and dimensions of the carotid plaque.

Impact of Tissue Damage and Hemodynamics on Restenosis Following Percutaneous Transluminal Angioplasty: A Patient-Specific Multiscale Model

Multiscale agent-based modeling frameworks have recently emerged as promising mechanobiological models to capture the interplay between biomechanical forces, cellular behavior, and molecular pathways underlying restenosis following percutaneous transluminal angioplasty (PTA). However, their applications are mainly limited to idealized scenarios. Herein, a multiscale agent-based modeling framework for investigating restenosis following PTA in a patient-specific superficial femoral artery (SFA) is proposed. The framework replicates the 2-month arterial wall remodeling in response to the PTA-induced injury and altered hemodynamics, by combining three modules: (i) the PTA module, consisting in a finite element structural mechanics simulation of PTA, featuring anisotropic hyperelastic material models coupled with a damage formulation for fibrous soft tissue and the element deletion strategy, providing the arterial wall damage and post-intervention configuration, (ii) the hemodynamics module, quantifying the post-intervention hemodynamics through computational fluid dynamics simulations, and (iii) the tissue remodeling module, based on an agent-based model of cellular dynamics. Two scenarios were explored, considering balloon expansion diameters of 5.2 and 6.2 mm. The framework captured PTA-induced arterial tissue lacerations and the post-PTA arterial wall remodeling. This remodeling process involved rapid cellular migration to the PTA-damaged regions, exacerbated cell proliferation and extracellular matrix production, resulting in lumen area reduction up to 1-month follow-up. After this initial reduction, the growth stabilized, due to the resolution of the inflammatory state and changes in hemodynamics. The similarity of the obtained results to clinical observations in treated SFAs suggests the potential of the framework for capturing patient-specific mechanobiological events occurring after PTA intervention.

Computational modeling of in-stent restenosis: Pharmacokinetic and pharmacodynamic evaluation

Persistence of the pathology of in-stent restenosis even with the advent of drug-eluting stents warrants the development of highly resolved in silico models. These computational models assist in gaining insights into the transient biochemical and cellular mechanisms involved and thereby optimize the stent implantation parameters. Within this work, an already established fully-coupled Lagrangian finite element framework for modeling the restenotic growth is enhanced with the incorporation of endothelium-mediated effects and pharmacological influences of rapamycin-based drugs embedded in the polymeric layers of the current generation drug-eluting stents. The continuum mechanical description of growth is further justified in the context of thermodynamic consistency. Qualitative inferences are drawn from the model developed herein regarding the efficacy of the level of drug embedment within the struts as well as the release profiles adopted. The framework is then intended to serve as a tool for clinicians to tune the interventional procedures patient-specifically.

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.

A multiphysics modeling approach for in-stent restenosis: Theoretical aspects and finite element implementation

Development of in silico models that capture progression of diseases in soft biological tissues are intrinsic in the validation of the hypothesized cellular and molecular mechanisms involved in the respective pathologies. In addition, they also aid in patient-specific adaptation of interventional procedures. In this regard, a fully-coupled high-fidelity Lagrangian finite element framework is proposed within this work which replicates the pathology of in-stent restenosis observed post stent implantation in a coronary artery. Advection-reaction-diffusion equations are set up to track the concentrations of the platelet-derived growth factor, the transforming growth factor-β, the extracellular matrix, and the density of the smooth muscle cells. A continuum mechanical description of volumetric growth involved in the restenotic process, coupled to the evolution of the previously defined vessel wall constituents, is presented. Further, the finite element implementation of the model is discussed, and the behavior of the computational model is investigated via suitable numerical examples. Qualitative validation of the computational model is presented by emulating a stented artery. Patient-specific data are intended to be integrated into the model to predict the risk of in-stent restenosis, and thereby assist in the tuning of stent implantation parameters to mitigate the risk.

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.

A predictive multiscale model of in-stent restenosis in femoral arteries: linking haemodynamics and gene expression with an agent-based model of cellular dynamics

In-stent restenosis (ISR) is a maladaptive inflammatory-driven response of femoral arteries to percutaneous transluminal angioplasty and stent deployment, leading to lumen re-narrowing as consequence of excessive cellular proliferative and synthetic activities. A thorough understanding of the underlying mechanobiological factors contributing to ISR is still lacking. Computational multiscale models integrating both continuous- and agent-based approaches have been identified as promising tools to capture key aspects of the complex network of events encompassing molecular, cellular and tissue response to the intervention. In this regard, this work presents a multiscale framework integrating the effects of local haemodynamics and monocyte gene expression data on cellular dynamics to simulate ISR mechanobiological processes in a patient-specific model of stented superficial femoral artery. The framework is based on the coupling of computational fluid dynamics simulations (haemodynamics module) with an agent-based model (ABM) of cellular activities (tissue remodelling module). Sensitivity analysis and surrogate modelling combined with genetic algorithm optimization were adopted to explore the model behaviour and calibrate the ABM parameters. The proposed framework successfully described the patient lumen area reduction from baseline to one-month follow-up, demonstrating the potential capabilities of this approach in predicting the short-term arterial response to the endovascular procedure.

Uncertainty quantification of a three-dimensional in-stent restenosis model with surrogate modelling

In-stent restenosis is a recurrence of coronary artery narrowing due to vascular injury caused by balloon dilation and stent placement. It may lead to the relapse of angina symptoms or to an acute coronary syndrome. An uncertainty quantification of a model for in-stent restenosis with four uncertain parameters (endothelium regeneration time, the threshold strain for smooth muscle cell bond breaking, blood flow velocity and the percentage of fenestration in the internal elastic lamina) is presented. Two quantities of interest were studied, namely the average cross-sectional area and the maximum relative area loss in a vessel. Owing to the high computational cost required for uncertainty quantification, a surrogate model, based on Gaussian process regression with proper orthogonal decomposition, was developed and subsequently used for model response evaluation in the uncertainty quantification. A detailed analysis of the uncertainty propagation is presented. Around 11% and 16% uncertainty is observed on the two quantities of interest, respectively, and the uncertainty estimates show that a higher fenestration mainly determines the uncertainty in the neointimal growth at the initial stage of the process. The uncertainties in blood flow velocity and endothelium regeneration time mainly determine the uncertainty in the quantities of interest at the later, clinically relevant stages of the restenosis process.

Multiscale Modeling of Vascular Remodeling Induced by Wall Shear Stress

Objective

Hemodynamics-induced low wall shear stress (WSS) is one of the critical reasons leading to vascular remodeling. However, the coupling effects of WSS and cellular kinetics have not been clearly modeled. The aim of this study was to establish a multiscale modeling approach to reveal the vascular remodeling behavior under the interaction between the macroscale of WSS loading and the microscale of cell evolution.

Methods

Computational fluid dynamics (CFD) method and agent-based model (ABM), which have significantly different characteristics in temporal and spatial scales, were adopted to establish the multiscale model. The CFD method is for the second/organ scale, and the ABM is for the month/cell scale. The CFD method was used to simulate blood flow in a vessel and obtain the WSS in a vessel cross-section. The simulations of the smooth muscle cell (SMC) proliferation/apoptosis and extracellular matrix (ECM) generation/degradation in a vessel cross-section were performed by using ABM. During the simulation of the vascular remodeling procedure, the damage index of the SMC and ECM was defined as deviation from the obtained WSS. The damage index decreased gradually to mimic the recovery of WSS-induced vessel damage.

Results

(1) The significant wall thickening region was consistent with the low WSS region. (2) There was no evident change of wall thickness in the normal WSS region. (3) When the damage index approached to 0, the amount and distribution of SMCs and ECM achieved a stable state, and the vessel reached vascular homeostasis.

Conclusion

The established multiscale model can be used to simulate the vascular remodeling behavior over time under various WSS conditions.

Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models

The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.

Do we really understand how drug eluted from stents modulates arterial healing?

The advent of drug-eluting stents (DES) has revolutionised the treatment of coronary artery disease. These devices, coated with anti-proliferative drugs, are deployed into stenosed or occluded vessels, compressing the plaque to restore natural blood flow, whilst simultaneously combating the evolution of restenotic tissue. Since the development of the first stent, extensive research has investigated how further advancements in stent technology can improve patient outcome. Mathematical and computational modelling has featured heavily, with models focussing on structural mechanics, computational fluid dynamics, drug elution kinetics and subsequent binding within the arterial wall; often considered separately. Smooth Muscle Cell (SMC) proliferation and neointimal growth are key features of the healing process following stent deployment. However, models which depict the action of drug on these processes are lacking. In this article, we start by reviewing current models of cell growth, which predominantly emanate from cancer research, and available published data on SMC proliferation, before presenting a series of mathematical models of varying complexity to detail the action of drug on SMC growth in vitro. Our results highlight that, at least for Sodium Salicylate and Paclitaxel, the current state-of-the-art nonlinear saturable binding model is incapable of capturing the proliferative response of SMCs across a range of drug doses and exposure times. Our findings potentially have important implications on the interpretation of current computational models and their future use to optimise and control drug release from DES and drug-coated balloons.

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

Healing in soft biological tissues is a chain of events on different time and length scales. This work presents a computational framework to capture and couple important mechanical, chemical and biological aspects of healing. A molecular-level damage in collagen, i.e., the interstrand delamination, is addressed as source of plastic deformation in tissues. This mechanism initiates a biochemical response and starts the chain of healing. In particular, damage is considered to be the stimulus for the production of matrix metalloproteinases and growth factors which in turn, respectively, degrade and produce collagen. Due to collagen turnover, the volume of the tissue changes, which can result either in normal or pathological healing. To capture the mechanisms on continuum scale, the deformation gradient is multiplicatively decomposed in inelastic and elastic deformation gradients. A recently proposed elasto-plastic formulation is, through a biochemical model, coupled with a growth and remodeling description based on homogenized constrained mixtures. After the discussion of the biological species response to the damage stimulus, the framework is implemented in a mixed nonlinear finite element formulation and a biaxial tension and an indentation tests are conducted on a prestretched flat tissue sample. The results illustrate that the model is able to describe the evolutions of growth factors and matrix metalloproteinases following damage and the subsequent growth and remodeling in the respect of equilibrium. The interplay between mechanical and chemo-biological events occurring during healing is captured, proving that the framework is a suitable basis for more detailed simulations of damage-induced tissue response.

Effects of local coronary blood flow dynamics on the predictions of a model of in-stent restenosis

Computational models are increasingly used to study cardiovascular disease. However, models of coronary vessel remodelling usually make some strong assumptions about the effects of a local narrowing on the flow through the narrowed vessel. Here, we test the effects of local flow dynamics on the predictions of an in-stent restenosis (ISR) model. A previously developed 2D model of ISR is coupled to a 1D model of coronary blood flow. Then, two different assumptions are tested. The first assumption is that the vasculature is always able to adapt, and the volumetric flow rate through the narrowed vessel is kept constant. The second, alternative, assumption is that the vasculature does not adapt at all, and the ratio of the pressure drop to the flow rate (hydrodynamic resistance) stays the same throughout the whole process for all vessels unaffected by the stenosis, and aortic or venous blood pressure does not change either. Then, the dynamics are compared for different locations in coronary tree for two different reendothelization scenarios. The assumptions of constant volumetric flow rate (absolute vascular adaptation) versus constant aortic pressure drop and no adaptation do not significantly affect the growth dynamics for most locations in the coronary tree, and the differences can only be observed at the locations where a strong alternative flow pathway is present. On the other hand, the difference between locations is significant, which is consistent with small vessel size being a risk factor for restenosis. These results suggest that the assumption of a constant flow is a good approximation for ISR models dealing with the typical progression of ISR in the most often stented locations such as the proximal parts of left anterior descending (LAD) and left circumflex (LCX) arteries.

[Review of studies on the biomechanical modelling of the coupling effect between stent degradation and blood vessel remodeling]

The dynamic coupling of stent degradation and vessel remodeling can influence not only the structural morphology and material property of stent and vessel, but also the development of in-stent restenosis. The research achievements of biomechanical modelling and analysis of stent degradation and vessel remodeling were reviewed; several noteworthy research perspectives were addressed, a stent-vessel coupling model was developed based on stent damage function and vessel growth function, and then concepts of matching ratio and risk factor were established so as to evaluate the treatment effect of stent intervention, which may lay the scientific foundation for the structure design, mechanical analysis and clinical application of biodegradable stent.

Finite element simulation and testing of cobalt-chromium stent: a parametric study on radial strength, recoil, foreshortening, and dogboning

The effectiveness of cardiovascular stenting procedure depends on the crimping and expansion characteristics of a stent, influenced by its design parameters. In this study, CoCr stents are fabricated, crimped on a tri-folded balloon, and expanded using manual inflation device. Similarly, in the finite element model, a tri-folded balloon is used to expand the stent. The length and diameter are measured to evaluate the radial strength, recoil, foreshortening, and dogboning. The simulation and experimental results match satisfactorily. The validated FE model can be used with confidence to optimize future stent designs, thus reducing the number of testing and product development time.

Computational modeling for cardiovascular tissue engineering: the importance of including cell behavior in growth and remodeling algorithms

Understanding cardiovascular growth and remodeling (G&R) is fundamental for designing robust cardiovascular tissue engineering strategies, which enable synthetic or biological scaffolds to transform into healthy living tissues after implantation. Computational modeling, particularly when integrated with experimental research, is key for advancing our understanding, predicting the in vivo evolution of engineered tissues, and efficiently optimizing scaffold designs. As cells are ultimately the drivers of G&R and known to change their behavior in response to mechanical cues, increasing efforts are currently undertaken to capture (mechano-mediated) cell behavior in computational models. In this selective review, we highlight some recent examples that are relevant in the context of cardiovascular tissue engineering and discuss the current and future biological and computational challenges for modeling cell-mediated G&R.

Mechanistic evaluation of long-term in-stent restenosis based on models of tissue damage and growth

Development and application of advanced mechanical models of soft tissues and their growth represent one of the main directions in modern mechanics of solids. Such models are increasingly used to deal with complex biomedical problems. Prediction of in-stent restenosis for patients treated with coronary stents remains a highly challenging task. Using a finite element method, this paper presents a mechanistic approach to evaluate the development of in-stent restenosis in an artery following stent implantation. Hyperelastic models with damage, verified with experimental results, are used to describe the level of tissue damage in arterial layers and plaque caused by such intervention. A tissue-growth model, associated with vessel damage, is adopted to describe the growth behaviour of a media layer after stent implantation. Narrowing of lumen diameter with time is used to quantify the development of in-stent restenosis in the vessel after stenting. It is demonstrated that stent designs and materials strongly affect the stenting-induced damage in the media layer and the subsequent development of in-stent restenosis. The larger the artery expansion achieved during balloon inflation, the higher the damage introduced to the media layer, leading to an increased level of in-stent restenosis. In addition, the development of in-stent restenosis is directly correlated with the artery expansion during the stent deployment. The correlation is further used to predict the effect of a complex clinical procedure, such as stent overlapping, on the level of in-stent restenosis developed after percutaneous coronary intervention.

A particle-based model for endothelial cell migration under flow conditions

Endothelial cells (ECs) play a major role in the healing process following angioplasty to inhibit excessive neointima. This makes the process of EC healing after injury, in particular EC migration in a stented vessel, important for recovery of normal vessel function. In that context, we present a novel particle-based model of EC migration and validate it against in vitro experimental data. We have developed a particle-based model of EC migration under flow conditions in an in vitro vessel with obstacles. Cell movement in the model is a combination of random walks and directed movement along the local flow velocity vector. For model calibration, a set of experimental data for cell migration in a similarly shaped channel has been used. We have calibrated the model for a baseline case of a channel with no obstacles and then applied it to the case of a channel with ridges on the bottom surface, representative of stent strut geometry. We were able to closely reproduce the cell migration speed and angular distribution of their movement relative to the flow direction reported in vitro. The model also reproduces qualitative aspects of EC migration, such as entrapment of cells downstream from the flow-disturbing ridge. The model has the potential, after more extensive in vitro validation, to study the effect of variation in strut spacing and shape, through modification of the local flow, on EC migration. The results of this study support the hypothesis that EC migration is strongly affected by the direction and magnitude of local wall shear stress.

Detection of arterial wall abnormalities via Bayesian model selection

Patient-specific modelling of haemodynamics in arterial networks has so far relied on parameter estimation for inexpensive or small-scale models. We describe here a Bayesian uncertainty quantification framework which makes two major advances: an efficient parallel implementation, allowing parameter estimation for more complex forward models, and a system for practical model selection, allowing evidence-based comparison between distinct physical models. We demonstrate the proposed methodology by generating simulated noisy flow velocity data from a branching arterial tree model in which a structural defect is introduced at an unknown location; our approach is shown to accurately locate the abnormality and estimate its physical properties even in the presence of significant observational and systemic error. As the method readily admits real data, it shows great potential in patient-specific parameter fitting for haemodynamical flow models.

Location-Specific Comparison Between a 3D In-Stent Restenosis Model and Micro-CT and Histology Data from Porcine In Vivo Experiments

Background

Coronary artery restenosis is an important side effect of percutaneous coronary intervention. Computational models can be used to better understand this process. We report on an approach for validation of an in silico 3D model of in-stent restenosis in porcine coronary arteries and illustrate this approach by comparing the modelling results to in vivo data for 14 and 28 days post-stenting.

Methods

This multiscale model includes single-scale models for stent deployment, blood flow and tissue growth in the stented vessel, including smooth muscle cell (SMC) proliferation and extracellular matrix (ECM) production. The validation procedure uses data from porcine in vivo experiments, by simulating stent deployment using stent geometry obtained from micro computed tomography (micro-CT) of the stented vessel and directly comparing the simulation results of neointimal growth to histological sections taken at the same locations.

Results

Metrics for comparison are per-strut neointimal thickness and per-section neointimal area. The neointimal area predicted by the model demonstrates a good agreement with the detailed experimental data. For 14 days post-stenting the relative neointimal area, averaged over all vessel sections considered, was 20 ± 3% in vivo and 22 ± 4% in silico. For 28 days, the area was 42 ± 3% in vivo and 41 ± 3% in silico.

Conclusions

The approach presented here provides a very detailed, location-specific, validation methodology for in silico restenosis models. The model was able to closely match both histology datasets with a single set of parameters. Good agreement was obtained for both the overall amount of neointima produced and the local distribution. It should be noted that including vessel curvature and ECM production in the model was paramount to obtain a good agreement with the experimental data.

Electronic supplementary material

The online version of this article (10.1007/s13239-019-00431-4) contains supplementary material, which is available to authorized users.

Mathematical modelling of the restenosis process after stent implantation

The stenting procedure has evolved to become a highly successful technique for the clinical treatment of advanced atherosclerotic lesions in arteries. However, the development of in-stent restenosis remains a key problem. In this work, a novel two-dimensional continuum mathematical model is proposed to describe the complex restenosis process following the insertion of a stent into a coronary artery. The biological species considered to play a key role in restenosis development are growth factors, matrix metalloproteinases, extracellular matrix, smooth muscle cells and endothelial cells. Diffusion-reaction equations are used for modelling the mass balance between species in the arterial wall. Experimental data from the literature have been used in order to estimate model parameters. Moreover, a sensitivity analysis has been performed to study the impact of varying the parameters of the model on the evolution of the biological species. The results demonstrate that this computational model qualitatively captures the key characteristics of the lesion growth and the healing process within an artery subjected to non-physiological mechanical forces. Our results suggest that the arterial wall response is driven by the damage area, smooth muscle cell proliferation and the collagen turnover among other factors.

Semi-intrusive multiscale metamodelling uncertainty quantification with application to a model of in-stent restenosis

We explore the efficiency of a semi-intrusive uncertainty quantification (UQ) method for multiscale models as proposed by us in an earlier publication. We applied the multiscale metamodelling UQ method to a two-dimensional multiscale model for the wound healing response in a coronary artery after stenting (in-stent restenosis). The results obtained by the semi-intrusive method show a good match to those obtained by a black-box quasi-Monte Carlo method. Moreover, we significantly reduce the computational cost of the UQ. We conclude that the semi-intrusive metamodelling method is reliable and efficient, and can be applied to such complex models as the in-stent restenosis ISR2D model. This article is part of the theme issue 'Multiscale modelling, simulation and computing: from the desktop to the exascale'.

Uncertainty Quantification of a Multiscale Model for In-Stent Restenosis

Purpose

Coronary artery stenosis, or abnormal narrowing, is a widespread and potentially fatal cardiac disease. After treatment by balloon angioplasty and stenting, restenosis may occur inside the stent due to excessive neointima formation. Simulations of in-stent restenosis can provide new insight into this process. However, uncertainties due to variability in patient-specific parameters must be taken into account.

Methods

We performed an uncertainty quantification (UQ) study on a complex two-dimensional in-stent restenosis model. We used a quasi-Monte Carlo method for UQ of the neointimal area, and the Sobol sensitivity analysis (SA) to estimate the proportions of aleatory and epistemic uncertainties and to determine the most important input parameters.

Results

We observe approximately 30% uncertainty in the mean neointimal area as simulated by the model. Depending on whether a fast initial endothelium recovery occurs, the proportion of the model variance due to natural variability ranges from 15 to 35%. The endothelium regeneration time is identified as the most influential model parameter.

Conclusion

The model output contains a moderate quantity of uncertainty, and the model precision can be increased by obtaining a more certain value on the endothelium regeneration time. We conclude that the quasi-Monte Carlo UQ and the Sobol SA are reliable methods for estimating uncertainties in the response of complicated multiscale cardiovascular models.

A versatile hybrid agent-based, particle and partial differential equations method to analyze vascular adaptation

Peripheral arterial occlusive disease is a chronic pathology affecting at least 8–12 million people in the USA, typically treated with a vein graft bypass or through the deployment of a stent in order to restore the physiological circulation. Failure of peripheral endovascular interventions occurs at the intersection of vascular biology, biomechanics, and clinical decision making. It is our hypothesis that the majority of endovascular treatment approaches share the same driving mechanisms and that a deep understanding of the adaptation process is pivotal in order to improve the current outcome of the procedure. The postsurgical adaptation of vein graft bypasses offers the perfect example of how the balance between intimal hyperplasia and wall remodeling determines the failure or the success of the intervention. Accordingly, this work presents a versatile computational model able to capture the feedback loop that describes the interaction between events at cellular/tissue level and mechano-environmental conditions. The work here presented is a generalization and an improvement of a previous work by our group of investigators, where an agent-based model uses a cellular automata principle on a fixed hexagonal grid to reproduce the leading events of the graft’s restenosis. The new hybrid model here presented allows a more realistic simulation both of the biological laws that drive the cellular behavior and of the active role of the membranes that separate the various layers of the vein. The novel feature is to use an immersed boundary implementation of a highly viscous flow to represent SMC motility and matrix reorganization in response to graft adaptation. Our implementation is modular, and this makes us able to choose the right compromise between closeness to the physiological reality and complexity of the model. The focus of this paper is to offer a new modular implementation that combines the best features of an agent-based model, continuum mechanics, and particle-tracking methods to cope with the multiscale nature of the adaptation phenomena. This hybrid method allows us to quickly test various hypotheses with a particular attention to cellular motility, a process that we demonstrated should be driven by mechanical homeostasis in order to maintain the right balance between cells and extracellular matrix in order to reproduce a distribution similar to histological experimental data from vein grafts.

A cell-based mechanical model of coronary artery tunica media

A three-dimensional cell-based mechanical model of coronary artery tunica media is proposed. The model is composed of spherical cells forming a hexagonal close-packed lattice. Tissue anisotropy is taken into account by varying interaction forces with the direction of intercellular connection. Several cell-centre interaction potentials for repulsion and attraction are considered, including the Hertz contact model and its neo-Hookean extension, the Johnson-Kendall-Roberts model of adhesive contact, and a wormlike chain model. The model is validated against data from in vitro uni-axial tension tests performed on dissected strips of tunica media. The wormlike chain potential in combination with the neo-Hookean Hertz contact model produces stress-stretch curves which represent the experimental data very well.