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

Plaque heterogeneity influences in-stent restenosis following drug-eluting stent implantation: Insights from patient-specific multiscale modelling

In-stent restenosis represents a major cause of failure of percutaneous coronary intervention with drug-eluting stent implantation. Computational multiscale models have recently emerged as powerful tools for investigating the mechanobiological mechanisms underlying vascular adaptation processes during in-stent restenosis. However, to date, the interplay between intervention-induced inflammation, drug delivery and drug retention has been under-investigated. Here, an original patient-specific multiscale agent-based modelling framework was developed to investigate the interplay between drug release, plaque composition and intervention-induced inflammation on in-stent restenosis following drug-eluting stent implantation. The framework integrated a finite element simulation of stent expansion, with a drug transport simulation and an agent-based model of cellular dynamics. A patient-specific coronary cross-section with heterogeneous diseased tissue was considered and rigorously analyzed through a variety of scenarios, including different plaque compositions and different inflammatory responses. The analysis revealed three significant findings: (i) calcifications substantially impeded drug transport, resulting in drug-depleted regions and reduced stent efficacy; (ii) by impacting drug transport, variations in plaque composition influenced arterial wall response, with the fully-calcific scenario showing the greatest lumen area reduction; (iii) the impact of different drug receptor saturation conditions (obtained with different plaque compositions) was particularly evident under conditions of persistent inflammatory state. This study represents a significant advancement in multiscale modelling of in-stent restenosis following drug-eluting stent implantation. The results obtained provided deeper insights into the complex interactions among patient-specific plaque composition, inflammation and drug retention, suggesting a patient-specific management of the intervention, particularly in cases of complex disease.

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.

Haemodynamics of stent-mounted neural interfaces in tapered and deformed blood vessels

The endovascular neural interface provides an appealing minimally invasive alternative to invasive brain electrodes for recording and stimulation. However, stents placed in blood vessels have long been known to affect blood flow (haemodynamics) and lead to neointimal growth within the blood vessel. Both the stent elements (struts and electrodes) and blood vessel wall geometries can affect the mechanical environment on the blood vessel wall, which could lead to unfavourable vascular remodelling after stent placement. With increasing applications of stents and stent-like neural interfaces in venous blood vessels in the brain, it is necessary to understand how stents affect blood flow and tissue growth in veins. We explored the haemodynamics of a stent-mounted neural interface in a blood vessel model. Results indicated that blood vessel deformation and tapering caused a substantial change to the lumen geometry and the haemodynamics. The neointimal proliferation was evaluated in sheep implanted with an endovascular neural interface. Analysis showed a negative correlation with the mean Wall Shear Stress pattern. The results presented here indicate that the optimal stent oversizing ratio must be considered to minimise the haemodynamic impact of stenting.

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.

Effects of stent shape on focal hemodynamics in intracranial atherosclerotic stenosis: A simulation study with computational fluid dynamics modeling

Background and aims

The shape of a stent could influence focal hemodynamics and subsequently plaque growth or in-stent restenosis in intracranial atherosclerotic stenosis (ICAS). In this preliminary study, we aim to investigate the associations between stent shapes and focal hemodynamics in ICAS, using computational fluid dynamics (CFD) simulations with manually manipulated stents of different shapes.

Methods

We built an idealized artery model, and reconstructed four patient-specific models of ICAS. In each model, three variations of stent geometry (i.e., enlarged, inner-narrowed, and outer-narrowed) were developed. We performed static CFD simulation on the idealized model and three patient-specific models, and transient CFD simulation of three cardiac cycles on one patient-specific model. Pressure, wall shear stress (WSS), and low-density lipoprotein (LDL) filtration rate were quantified in the CFD models, and compared between models with an inner- or outer-narrowed stent vs. an enlarged stent. The absolute difference in each hemodynamic parameter was obtained by subtracting values from two models; a normalized difference (ND) was calculated as the ratio of the absolute difference and the value in the enlarged stent model, both area-averaged throughout the arterial wall.

Results

The differences in focal pressure in models with different stent geometry were negligible (ND<1% for all cases). However, there were significant differences in the WSS and LDL filtration rate with different stent geometry, with ND >20% in a static model. Observable differences in WSS and LDL filtration rate mainly appeared in area adjacent to and immediately distal to the stent. In the transient simulation, the LDL filtration rate had milder temporal fluctuations than WSS.

Conclusions

The stent geometry might influence the focal WSS and LDL filtration rate in ICAS, with negligible effect on pressure. Future studies are warranted to verify the relevance of the changes in these hemodynamic parameters in governing plaque growth and possibly in-stent restenosis in ICAS.

The endothelium permeability after bioresorbable scaffolds implantation caused by the heterogeneous expression of tight junction proteins

As one of the main functions of vascular endothelial cells, Vascular permeability is determined by four tight junction proteins (TJPs): Zonula Occludens-1 (ZO-1), Claudin-5, Occludin and Tricellulin. The barrier function of blood vessels will be reconstructed after they are damaged by endothelial mechanical injuries caused by vascular interventions. In this study, the effects of balloon expansion (transient mechanical injury) on four TJPs and vascular permeability were compared with those of poly-l-lactic acid bioresorbable scaffolds (BRSs) implantation (continuous mechanical stimulation). We found that BRSs do not affect vascular permeability, while the recovery of vascular barrier function was found to be only related to the mechanical injuries and repair of endothelium. Mechanical stimulation affects and accelerates the recovery process of vascular permeability with the heterogeneous expression levels of TJPs induced after BRSs implantation. Different TJPs have different sensitivity to different loyal mechanical stimuli. ZO-1 is more sensitive to shear stress and tension than to static pressure. Occludin is sensitive to static pressure and shear stress. Tricellulin is more sensitive to tension stretching. Compared with the other three TJPs, Claudin-5 can respond to mechanical stimulation, with relatively low sensitivity, though. This difference in sensitivity determines the heterogeneous expression of TJPs. Mechanical stimulation of different kinds and strengths can also cause different cell morphological changes and inflammatory reactions. As an important element affecting endothelial function, the mechanical factors emerging after BRSs implantation are worthy of more attention.

•

The repair of vascular permeability is directly related to the type of vascular injuries, while BRSs implantation has little effect on vascular permeability.

•

Transient and persistent mechanical stimulation is the main reason to influence the expression of TJPs.

•

Heterogeneous expression of TJPs caused by their different sensitivity to the form of mechanical stimuli.

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 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.

Assessment of single and double coronary bifurcation stenting techniques using multimodal imaging and 3D modeling in reanimated swine hearts using Visible Heart® methodologies

Stent implantation in bifurcated coronary lesions is technically challenging so that procedural refinements are continuously investigated. Novel procedure modeling and intracoronary imaging techniques may offer critical insights on stent deformations and stent-wall interactions during bifurcation stenting procedures. Thus, we assessed coronary bifurcation stenting techniques using multimodal imaging and 3D modeling in reanimated swine hearts. Harvested swine hearts were reanimated using Visible Heart® methodologies and (under standard fluoroscopic guidance) used to test 1-stent (provisional and inverted provisional) and 2-stent (culotte, TAP and DK-crush) techniques on bifurcations within various coronary vessels using commercially available devices. Intracoronary angioscopy and frequency-domain optical-coherence-tomography (OCT) were obtained during the procedures. 3D OCT reconstruction and micro-computed tomography 3D modeling (post heart fixations) were used to assess stent deformations and stent-wall interactions. We conducted multiple stenting procedures and collected unique endoscopic and OCT images (and subsequent computational models from micro-CT) to assess stent deformations and device/wall interactions during different steps of bifurcation stenting procedures. Endoscopy, micro-CT and virtual reality processing documented that different 1- and 2-stent techniques, practiced according to experts’ recommended steps, achieve optimal post-intervention stent conformation. As compared with intra-procedural endoscopy, software-generated 3D OCT images accurately depicted stent deformations during 1-stent techniques. On the opposite, during more complex 2-stent techniques, some defects were appreciated at 3D OCT reconstruction despite optimal 2D OCT images. This study provided unique insights regarding both stent deformations occurring in the course of bifurcation stenting and the efficacy of OCT to visualize them.

In-Stent Restenosis Progression in Human Superficial Femoral Arteries: Dynamics of Lumen Remodeling and Impact of Local Hemodynamics

In-stent restenosis (ISR) represents a major drawback of stented superficial femoral arteries (SFAs). Motivated by the high incidence and limited knowledge of ISR onset and development in human SFAs, this study aims to (i) analyze the lumen remodeling trajectory over 1-year follow-up period in human stented SFAs and (ii) investigate the impact of altered hemodynamics on ISR initiation and progression. Ten SFA lesions were reconstructed at four follow-ups from computed tomography to quantify the lumen area change occurring within 1-year post-intervention. Patient-specific computational fluid dynamics simulations were performed at each follow-up to relate wall shear stress (WSS) based descriptors with lumen remodeling. The largest lumen remodeling was found in the first post-operative month, with slight regional-specific differences (larger inward remodeling in the fringe segments, p < 0.05). Focal re-narrowing frequently occurred after 6 months. Slight differences in the lumen area change emerged between long and short stents, and between segments upstream and downstream from stent overlapping portions, at specific time intervals. Abnormal patterns of multidirectional WSS were associated with lumen remodeling within 1-year post-intervention. This longitudinal study gave important insights into the dynamics of ISR and the impact of hemodynamics on ISR progression in human SFAs.

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.

Effects of different positions of intravascular stent implantation in stenosed vessels on in-stent restenosis: An experimental and numerical simulation study

Percutaneous coronary intervention (PCI) has been widely used in the treatment of atherosclerosis, while in-stent restenosis (ISR) has not been completely resolved. Studies have shown that changes in intravascular mechanical environment are related to ISR. Hence, an in-depth understanding of the effects of stent intervention on vascular mechanics is important for clinically optimizing stent implantation and relieving ISR. Nine rabbits with stenotic carotid artery were collected by balloon injury. Intravascular stents were implanted into different longitudinal positions (proximal, middle and distal relative to the stenotic area) of the stenotic vessels for numerical simulations. Optical coherence tomography (OCT) scanning was performed to reconstruct the three-dimensional configuration of the stented carotid artery and blood flow velocity waveforms were collected by Doppler ultrasound. The numerical simulations were performed through direct solution of Naiver-Stokes equation in ANSYS. Results showed that the distributions of time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and relative residual time (RRT) in near-end segment were distinctively different from other regions of the stent which considered to promote restenosis for all three models. Spearman rank-correlation analysis showed a significant correlation between hemodynamic descriptors and the stent longitudinal positions (r_TAWSS = -0.718, r_OSI = 0.898, r_RRT = 0.818, p < 0.01). Histology results of the near-end segment showed neointima thickening deepened with the longitudinal positions of stent which was consistent with the numerical simulations. The results suggest that stent implantation can promote restenosis at the near-end segment. As the stenting position moves to distal end, the impact on ISR is more significant.

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.