Numerical Simulation of Stent Angioplasty with Predilation: An Investigation into Lesion Constitutive Representation and Calcification Influence

Drug coated balloons in percutaneous coronary intervention: how can computational modelling help inform evolving clinical practice?

Drug-coated balloons (DCB) represent an emerging therapeutic alternative to drug-eluting stents (DES) for the treatment of coronary artery disease (CAD). Among the key advantages of DCB over DES are the absence of a permanent structure in the vessel and the potential for fast and homogeneous drug delivery. While DCB were first introduced for treatment of in-stent restenosis (ISR), their potential wider use in percutaneous coronary intervention (PCI) has recently been explored in several randomized clinical trials, including for treatment of de novo lesions. Moreover, new hybrid techniques that combine DES and DCB are being investigated to more effectively tackle complex cases. Despite the growing interest in DCB within the clinical community, the mechanisms of drug exchange and the interactions between the balloon, the polymeric coating and the vessel wall are yet to be fully understood. It is, therefore, perhaps surprising that the number of computational (in silico) models developed to study interventions involving these devices is small, especially given the mechanistic understanding that has been gained from computational studies of DES procedures over the last two decades. In this paper, we discuss the current and emerging clinical approaches for DCB use in PCI and review the computational models that have been developed thus far, underlining the potential challenges and opportunities in integrating in silico models of DCB into clinical practice.

Finite element analysis predicts a major mechanical role of epicardial adipose tissue in atherosclerotic coronary disease and angioplasty

BACKGROUND

Understanding how atherosclerosis and angioplasty biomechanically affect the coronary artery wall is crucial for comprehending the pathophysiology of this disease and advancing potential treatments. However, acquiring this information experimentally or in vivo presents challenges. To overcome this, different computational methods have been employed. This research assessed the impact of atherosclerosis and angioplasty on the strains of each coronary artery tunic using the finite element method.

METHODS

Anatomical data were used to create two three-dimensional models of the left anterior descending coronary artery: one representing a normal artery and the other with concentric atherosclerosis, which included the surrounding epicardial fat tissue (EFT) and the three arterial tunics (e.g., intima, media, and adventitia). Blood pressure was applied to both models, and angioplasty was performed in the atherosclerotic model. The mean maximum principal and minimum principal strains were obtained for each layer in each case, and the impact of EFT was analyzed by comparing the results of including and omitting it. Furthermore, a sensitivity analysis was conducted for EFT stiffness, EFT volume, and blood pressure.

RESULTS

Noteworthy biomechanical alterations were observed in the atherosclerotic model before and after angioplasty, compared to the healthy state. After angioplasty, strains in the media and adventitia layers increased on average by up to fivefold, whereas the intima layer experienced a comparatively lower impact. Similarly, excluding EFT resulted in an average fourfold increase in strains in the tunics of both the healthy and atherosclerotic models. In addition, in both healthy and atherosclerotic models, a rise in blood pressure caused the most significant increase in arterial tunic strains, followed by reduced EFT stiffness and increased EFT volume, in order of impact.

CONCLUSION

Coronary artery wall strains are significantly altered by atherosclerosis and angioplasty, leading to cellular growth in the media and adventitia layers and subsequent reobstruction of the lumen after the procedure. EFT strongly influences coronary wall biomechanics, with low EFT stiffness and high volume predicted as risk factors for the development and severity of atherosclerosis. However, all the above may be modulated through interventions targeting epicardial adipose tissue.

Load-sharing characteristics of stenting and post-dilation in heavily calcified coronary artery

In this work, stenting in non-calcified and heavily calcified coronary arteries was quantified in terms of diameter-pressure relationships and load transfer from the balloon to the artery. The efficacy of post-dilation in non-calcified and heavily calcified coronary arteries was also characterized in terms of load sharing and the changes in tissue mechanics. Our results have shown that stent expansion exhibits a cylindrical shape in non-calcified lesions, while it exhibits a dog bone shape in heavily calcified lesions. Load-sharing analysis has shown that only a small portion of the pressure load (1.4 N, 0.8% of total pressure load) was transferred to the non-calcified lesion, while a large amount of the pressure load (19 N, 12%) was transferred to the heavily calcified lesion. In addition, the increasing inflation pressure (from 10 to 20 atm) can effectively increase the minimal lumen diameter (from 1.48 to 2.82 mm) of the heavily calcified lesion, the stress (from 1.5 to 8.4 MPa) and the strain energy in the calcification (1.77 mJ to 26.5 mJ), which are associated with the potential of calcification fracture. Results indicated that increasing inflation pressure can be an effective way to improve the stent expansion if a dog bone shape of the stenting profile is observed. Considering the risk of a balloon burst, our results support the design and application of the high-pressure balloon for post-dilation. This work also sheds some light on the stent design and choice of stent materials for improving the stent expansion at the dog bone region and mitigating stresses on arterial tissues.

Load-sharing characteristics of stenting and post-dilation in heavily calcified coronary artery

In this work, stenting in non-calcified and heavily calcified coronary arteries was quantified in terms of diameter-pressure relationships and load transfer from the balloon to the artery. The efficacy of post-dilation in non-calcified and heavily calcified coronary arteries was also characterized in terms of load sharing and the changes in tissue mechanics. Our results have shown that stent expansion exhibits a cylindrical shape in non-calcified lesions, while it exhibits a dog bone shape in heavily calcified lesions. Load-sharing analysis has shown that only a small portion of the pressure load (1.4 N, 0.8% of total pressure load) was transferred to the non-calcified lesion, while a large amount of the pressure load (19 N, 12%) was transferred to the heavily calcified lesion. In addition, the increasing inflation pressure (from 10 to 20 atm) can effectively increase the minimal lumen diameter (from 1.48 mm to 2.82 mm) of the heavily calcified lesion, the stress (from 1.5 MPa to 8.4 MPa) the strain energy in the calcification (1.77 mJ to 26.5 mJ), which associated with the potential of calcification fracture. Results indicated that increasing inflation pressure can be an effective way to improve the stent expansion if a dog bone shape of the stenting profile is observed. Considering the risk of a balloon burst, our results support the design and application of the high-pressure balloon for post-dilation.

A parametric study assessing Implicit Solver limits for a generic FEM Simulation of PTA without stent deployment

This study focuses on the robustness of a generic Finite Element Model (FEM) of Percutaneous Transluminal Angioplasty (PTA) procedure with permanent set. The influence of three different parameters on simulation robustness were investigated: the stenosis percent, the stenosis offset and the arterial caliber. Five arterial calibers are modeled by adapting the ratio between the inner diameter and the wall thickness. Overall, forty configurations were tested with the same simulation settings and boundary conditions. Results shows convergence issues caused by excessive deformations of elements for stenosis above 65% blockage. Moreover, an increasing stenosis offset tends to decrease convergence. Simulation of PTA on small calibers and large calibers are less robust than intermediate e.g., iliac calibers.Clinical Relevance- PTA can benefit from numerical tools to improve the procedure outcomes. A FEM simulation of PTA without stent deployment can predict the permanent strain induced by this surgery for various configurations. However, robustness of the simulation is required to consider its transfer to clinics. This work aims to determine the robustness boundaries of an implicit solver for PTA simulation. It shows that an implicit solver is robust for all artery calibers with a stenosis below 50% blockage. Moreover medium-caliber arteries exhibit better robustness with converging solutions for stenosis reaching 60% blockage.

Effect of balloon pre-dilation on performance of self-expandable nitinol stent in femoropopliteal artery

Balloon pre-dilation is usually performed before implantation of a nitinol stent in a femoropopliteal artery in a case of severe blockage or calcified plaque. However, its effect on performance of the nitinol stent in a diseased femoropopliteal artery has not been studied yet. This study compares the outcomes of stenting with pre-dilation and without it by modelling the entire processes of stent deployment. Fatigue deformation of the implanted stent is also modelled under diastolic-systolic blood pressure, repetitive bending, torsion, axial compression and their combination. Reduced level of stress in the stent occurs after stenting with pre-dilation, but causing the increased damage in the media layer, i.e. the middle layer of the arterial wall. Generally, pre-dilation increases the risk of nitinol stent's fatigue failure. Additionally, the development of in-stent restenosis is predicted based on the stenting-induced tissue damage in the media layer, and no severe mechanical irritation is induced to the media layer by pre-dilation, stent deployment or fatigue loading.

Towards a Digital Twin of Coronary Stenting: A Suitable and Validated Image-Based Approach for Mimicking Patient-Specific Coronary Arteries

Considering the field of application involving stent deployment simulations, the exploitation of a digital twin of coronary stenting that can reliably mimic the patient-specific clinical reality could lead to improvements in individual treatments. A starting step to pursue this goal is the development of simple, but at the same time, robust and effective computational methods to obtain a good compromise between the accuracy of the description of physical phenomena and computational costs. Specifically, this work proposes an approach for the development of a patient-specific artery model to be used in stenting simulations. The finite element model was generated through a 3D reconstruction based on the clinical imaging (coronary Optical Coherence Tomography (OCT) and angiography) acquired on the pre-treatment patient. From a mechanical point of view, the coronary wall was described with a suitable phenomenological model, which is consistent with more complex constitutive approaches and accounts for the in vivo pressurization and axial pre-stretch. The effectiveness of this artery modeling method was tested by reproducing in silico the stenting procedures of two clinical cases and comparing the computational results with the in vivo lumen area of the stented vessel.

Investigation of mechanical behaviors and improved design of V-shaped braid stents

V-shaped braid stents (VBSs), as highly retrievable and flexible nitinol stents, are extensively applied in endovascular diseases. They also cause less damage to vessel wall compared to tube-cutting stents. However, poor performance of VBS or suboptimal operation can give rise to unwanted clinical situations such as thrombosis and intimal hyperplasia. Therefore, research on designing factors affecting the performance of these devices is of great significance. Furthermore, simulation of stenting process can help designers understand the interactions of stents and vessel wall to reduce time to market. Thus, finite element analysis (FEA) and bench test are performed taking into account both designing factors and stenting process of VBS, including development of parametric modeling tool, research on the relationships among structural parameters and radial force, exploration of the interactions of VBS and vessel wall and pulsating load effect. This research was performed using a commercial solver Abaqus/standard with a user material subroutine (UMAT/nitinol). Structural parameters of VBS, unit-cell height and wire diameter have significant impacts on radial force, unit-cell number has slight influence on radial force, and arc diameter has almost negligible impact on radial force. Without pulsatile load, maximum stress and strain always occur in arc position; however, in pulsatile load, maximum stress and strain are gradually transformed to strut position. The stress created near vessel wall and VBS interface is higher than interaction stress due to pulsating load. The obtained result provided valuable information on the structural design of stents as well as the effects of stent on vessel wall and that vessel wall on stent deformation.Graphical abstract[Formula: see text].

Acute Stent-Induced Endothelial Denudation: Biomechanical Predictors of Vascular Injury

Recent concern for local drug delivery and withdrawal of the first Food and Drug Administration-approved bioresorbable scaffold emphasizes the need to optimize the relationships between stent design and drug release with imposed arterial injury and observed pharmacodynamics. In this study, we examine the hypothesis that vascular injury is predictable from stent design and that the expanding force of stent deployment results in increased circumferential stress in the arterial tissue, which may explain acute injury poststent deployment. Using both numerical simulations and ex vivo experiments on three different stent designs (slotted tube, corrugated ring, and delta wing), arterial injury due to device deployment was examined. Furthermore, using numerical simulations, the consequence of changing stent strut radial thickness on arterial wall shear stress and arterial circumferential stress distributions was examined. Regions with predicted arterial circumferential stress exceeding a threshold of 49.5 kPa compared favorably with observed ex vivo endothelial denudation for the three considered stent designs. In addition, increasing strut thickness was predicted to result in more areas of denudation and larger areas exposed to low wall shear stress. We conclude that the acute arterial injury, observed immediately following stent expansion, is caused by high circumferential hoop stresses in the interstrut region, and denuded area profiles are dependent on unit cell geometric features. Such findings when coupled with where drugs move might explain the drug–device interactions.

Influence of balloon design, plaque material composition, and balloon sizing on acute post angioplasty outcomes: An implicit finite element analysis

In this work we propose a generic modeling approach for simulating percutaneous transluminal angioplasty (PTA) endovascular treatment, and evaluating the influence of balloon design, plaque composition, and balloon sizing on acute post-procedural outcomes right after PTA, without stent implantation. Clinically-used PTA balloons were classified into two categories according to their compliance characteristics, and were modeled correspondingly. Self-defined elastoplastic constitutive laws were implemented within the plaque and artery models, after calibration based on experimental and clinical data. Finite element method (FEM) implicit solver was used to simulate balloon inflation and deflation. Besides balloon profile at max inflation, results are mainly assessed in terms of the elastic recoil ratio (ERR) and lumen gain ratio (LGR) obtained immediately after PTA. No variations in ERR nor LGR values were detected when the balloon design changed, despite the differences observed in their profile at max inflation. Moreover, LGR and ERR inversely varied with the augmentation of calcification level within the plaque (-11% vs. +4% respectively, from fully lipidic to fully calcified plaque). Furthermore, results showed a direct correlation between balloon sizing and LGR and ERR, with noticeably higher rates of change for LGR (+18% and +2% for LGR and ERR respectively for a calcified plaque and a balloon pressure increasing from 10 to 14 atm). However a larger LGR comes with a higher risk of arterial rupture. This proposed methodology opens the way for evaluation of angioplasty balloon selections towards clinical procedure optimization.

Validation of the computational model of a coronary stent: a fundamental step towards in silico trials

The proof of the reliability of a numerical model is becoming of paramount importance in the era of in silico clinical trials. When dealing with a coronary stenting procedure, the virtual scenario should be able to replicate the real device, passing through the different stages of the procedure, which has to maintain the atherosclerotic vessel opened. Nevertheless, most of the published studies adopted commercially resembling geometries and generic material parameters, without a specific validation of the employed numerical models. In this work, a workflow for the generation and validation of the computational model of a coronary stent was proposed. Possible sources of variability in the results, such as the inter-batches variability in the material properties and the choice of proper simulation strategies, were accounted for and discussed. Then, a group of in vitro tests, representative of the device intended use was used as a comparator to validate the model. The free expansion simulation, which is the most used simulation in the literature, was shown to be only partially useful for stent model validation purposes. On the other hand, the choice of proper additional experiments, as the suggested uniaxial tensile tests on the stent and deployment tests into a deformable tube, could provide further suitable information to prove the efficacy of the numerical approach.

Mechanical performances of balloon post-dilation for improving stent expansion in calcified coronary artery: Computational and experimental investigations

Stent deployment in a calcified coronary artery is often associated with suboptimal outcomes such as stent underexpansion and malapposition. Post-dilation after stent deployment is commonly used for optimal stent implantation. There is no guideline for choosing the post-dilation balloon diameter and inflation pressure. In this work, ex-vivo/in-silico experiments were performed to investigate the efficacy of post-dilation balloon diameter and inflation pressure in improving the stent expansion in a calcified lesion. Post-dilations with three balloon diameters (3 mm, 3.5 mm, and 4 mm) were performed. For each balloon diameter, three inflation pressures (10 atm, 20 atm, and 30 atm) were sequentially applied. In ex-vivo experiments, optical coherence tomography images were acquired during the stenting procedure, i.e., pre- and post-deployment of 3 mm diameter stent, as well as after each post-dilation. The results from in-silico experiments were compared with ex-vivo experiments in terms of lumen area. In addition, stretch ratio analysis was developed to predict the stent-induced lumen area, along with the strain analysis and the in-silico experiments. Results have shown that target lumen area could be achieved with an oversized nominal balloon diameter of +0.5 mm (i.e., 0.5 mm greater than reference lumen diameter) at an inflation pressure of 20 atm. After each post-dilation, fibrotic tissue demonstrated a larger strain, contributing to improved lumen gain. However, minimal changes were observed in calcification. Moreover, a strong correlation (R² = 0.95) between the stretch ratio of fibrotic tissue and lumen area after each post-dilation was observed. This indicated that the morphology of the fibrotic tissue could be a potential marker to predict the lumen gain. The detailed mechanistic quantifications of a single lesion cannot be generalized to all clinical cases. However, this work could be used to provide a fundamental understanding of the post-dilations, to develop experimental protocols for producing generalized guidelines, and to exploit their potential for optimal pre- and post-stent strategies.

The influence of angioplasty balloon sizing on acute post-procedural outcomes: a Finite Element Analysis

Atherosclerosis is one of the most common vascular pathologies in the world. Among the most commonly performed endovascular treatments, percutaneous transluminal angioplasty (PTA) has been showing significantly positive clinical outcomes. Due to the complex geometries, material properties and interactions that characterize PTA procedures, finite element analyses of acute angioplasty balloon deployment are limited. In this work, finite element method (FEM) was used to simulate the inflation and deflation of a semi-compliant balloon within the 3D model of a stenosed artery with two different plaque types (lipid and calcified). Self-defined constitutive models for the balloon and the plaque were developed based on experimental and literature data respectively. Balloon deployment was simulated at three different inflation pressures (10, 12 and 14 atm) within the two plaque types. Balloon sizing influence on the arterial elastic recoil obtained immediately after PTA was then investigated. The simulated results show that calcified plaques may lead to higher elastic recoil ratios compared to lipid stenosis, when the same balloon inflation pressures are applied. Also, elastic recoil increases for higher balloon inflation pressure independent of the plaque type. These findings open the way for a data-driven assessment of angioplasty balloon sizing selection and clinical procedures optimization.Clinical Relevance- The FE model developed in this work aims at providing quantitative evaluation of recoil after balloon angioplasty. It may be useful for both manufacturers and clinicians to improve efficiency of angioplasty balloon device design and sizing selection with respect to plaque geometry and constitution, consequently enhancing clinical outcomes.

In vivo and in vitro evaluation of a biodegradable magnesium vascular stent designed by shape optimization strategy

The performance of biodegradable magnesium alloy stents (BMgS) requires special attention to non-uniform residual stress distribution and stress concentration, which can accelerate localized degradation after implantation. We now report on a novel concept in stent shape optimization using a finite element method (FEM) toolkit. A Mg-Nd-Zn-Zr alloy with uniform degradation behavior served as the basis of our BMgS. Comprehensive in vitro evaluations drove stent optimization, based on observed crimping and balloon inflation performance, measurement of radial strength, and stress condition validation via microarea-XRD. Moreover, a Rapamycin-eluting polymer coating was sprayed on the prototypical BMgS to improve the corrosion resistance and release anti-hyperplasia drugs. In vivo evaluation of the optimized coated BMgS was conducted in the iliac artery of New Zealand white rabbit with quantitative coronary angiography (QCA), optical coherence tomography (OCT) and micro-CT observation at 1, 3, 5-month follow-ups. Neither thrombus or early restenosis was observed, and the coated BMgS supported the vessel effectively prior to degradation and allowed for arterial healing thereafter. The proposed shape optimization framework based on FEM provides an novel concept in stent design and in-depth understanding of how deformation history affects the biomechanical performance of BMgS. Computational analysis tools can indeed promote the development of biodegradable magnesium stents.

Computational modelling of magnesium stent mechanical performance in a remodelling artery: Effects of multiple remodelling stimuli

Significant research has been conducted in the area of coronary stents/scaffolds made from resorbable metallic and polymeric biomaterials. These next-generation bioabsorbable stents have the potential to completely revolutionise the treatment of coronary artery disease. The primary advantage of resorbable devices over permanent stents is their temporary presence which, from a theoretical point of view, means only a healed coronary artery will be left behind following degradation of the stent potentially eliminating long-term clinical problems associated with permanent stents. The healing of the artery following coronary stent/scaffold implantation is crucial for the long-term safety of these devices. Computational modelling can be used to evaluate the performance of complex stent devices in silico and assist in the design and development and understanding of the next-generation resorbable stents. What is lacking in computational modelling literature is the representation of the active response of the arterial tissue in the weeks and months following stent implantation, ie, neointimal remodelling, in particular for the case of biodegradable stents. In this paper, a computational modelling framework is developed, which accounts for two major physiological stimuli responsible for neointimal remodelling and combined with a magnesium corrosion model that is capable of simulating localised pitting (realistic) stent corrosion. The framework is used to simulate different neointimal growth patterns and to explore the effects the neointimal remodelling has on the mechanical performance (scaffolding support) of the bioabsorbable magnesium stent.

Finite element evaluation of artery damage in deployment of polymeric stent with pre- and post-dilation

Using finite element method, this paper evaluates damage in an arterial wall and plaque caused by percutaneous coronary intervention. Hyperelastic damage models, calibrated with experimental results, are used to describe stress–stretch responses of arterial layers and plaque; these models are capable to simulate softening behaviour of the tissue due to damage. Abaqus CAE is employed to create the finite element models for the artery wall (with media and adventitia layers), a symmetric uniform plaque, a bioresorbable polymeric stent and a tri-folded expansion balloon. The effect of percutaneous coronary intervention on vessel damage is investigated by simulating the processes of vessel pre-dilation, stent deployment and post-stenting dilation. Energy dissipation density is used to assess the extent of damage in the tissue. Softening of the plaque and the artery, due to the pre-dilation-induced damage, can facilitate the subsequent stent deployment process. The plaque and the artery experienced heterogeneous damage behaviour after the stent deployment, caused by non-uniform deformation. The post-stenting dilation was effective to achieve a full expansion of the stent, but caused additional damage to the artery. The continuous and discontinuous damage models yielded similar results in the percutaneous coronary intervention simulations, while the incorporation of plaque rupture affected the simulated outcomes of stent deployment. The computational evaluation of the artery damage can be potentially used to assess the risk of in-stent restenosis after percutaneous coronary intervention.

Intravascular optical coherence tomography to validate finite-element simulation of angioplasty balloon inflation

Concrete methods are lacking to examine angioplasty simulation results. For the first time, we explored the application of intravascular optical coherence tomography (IVOCT) to experimentally validate results obtained from finite-element simulation of angioplasty balloon deployment. In order to simulate each experimental scenario, IVOCT images were used to create initial geometrical models for the balloon and the phantoms. The study comprised three scenarios. The first scenario involved experimentally monitoring as well as simulating free expansion of the balloon. The second scenario involved experimentally monitoring as well as simulating balloon inflation inside three artery phantoms with different mechanical properties. The third scenario involved experimentally monitoring as well as simulating balloon unfolding and inflation inside a multilayer phantom. Using the first scenario, a constitutive model was developed for the balloon and was tuned to fit the IVOCT balloon inflation monitoring results. This model was used to simulate the balloon's response in simulations involving phantoms corresponding to the second and third scenarios. Diameter values were calculated both from images and simulation results. These values were then compared for each scenario. The obtained results highlight the potentials of IVOCT monitoring to validate simulation procedures.

The influence of the nylon balloon stiffness on the efficiency of the intra-aortic balloon occlusion

In interventional procedures, the balloon inflation is used to occlude the artery and thus reduce bleeding. There is no practically accepted measure of the procedure efficiency. A finite element method model with state-of-the-art modelling techniques was built in order to predict the occlusion levels under the influence of different balloon inflation and its material stiffness. The geometries of a healthy human thoracic aorta and an occlusion balloon were idealized. The non-linear constitutive material of Gasser-Ogden-Holzapfel model was employed for the thoracic aorta; the balloon was model as the hyperelastic model. The realistic physiological blood pressure and the balloon inflation pressures were applied to simulate the different occlusion levels. The final outcome shows an important influence of the material stiffness on the balloon deformation and thus the occlusion efficiency.