Computational model of stent-based delivery from a half-embedded two-layered coating

Plaque heterogeneity and the spatial distributions of its components dictate drug-coated balloon therapy

Drug-coated balloon (DCB) angioplasty is one of the potential approaches to alleviating in-stent restenosis and treating peripheral artery disease. An in-silico model has been developed for sirolimus drug eluted from an inflated balloon in a patient-specific arterial cross-section consisting of fibrous tissue, fibrofatty tissue, dense calcium, necrotic core, and healthy tissue. The convection-diffusion-reaction equation represents the transport of drug, while drug binding, both specific and non-specific, can be modelled as a reaction process. The Brinkman equations describe the interstitial flow in porous tissue. An image processing technique is leveraged for reconstructing the computational domain. The Marker and Cell, and Immersed Boundary Methods are used to solve the set of governing equations. The no-flux interface condition and convection do amplify the tissue content, and the regions of dense calcium and necrotic core limited to or extremely close to the interface pose a clinical threat to DCB therapy. Simulations predict the effects of the positioning and clustering of plaque components in the domain. This study demands extensive intravascular ultrasound-derived virtual histology (VH-IVUS) imaging to understand the plaque morphology and determine the relative positions of different plaque compositions about the lumen-tissue interface, which have a significant impact on arterial pharmacokinetics.

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.

An unsteady analysis of two-phase binding of drug in an asymmetric stenosed vessel

In this paper, we investigate endovascular delivery to get a step ahead of the pharmacological limitations it has due to the complexity of dealing with a patient-specific vessel through a mathematical model. We divide the domain of computation into four sub-domains: the lumen, the lumen-tissue interface, the upper tissue and the lower tissue which are extracted from an asymmetric atherosclerotic image derived by the intravascular ultrasound (IVUS) technique. The injected drug at the luminal inlet is transported with the streaming blood which is considered Newtonian. An irreversible uptake kinetics of the injected drug at the lumen-tissue interface from the luminal side to the tissue domains is assumed. Subsequently, the drug is dispersed within the tissue followed by its retention in the extracellular matrix (ECM) and by receptor-mediated binding. The Marker and Cell (MAC) method has been leveraged to get a quantitative insight into the model considered. The effect of the wall absorption parameter on the concentration of all drug forms (free as well as two-phase bound) has been thoroughly investigated, and some other important factors, such as the averaged concentration, the tissue content, the fractional effect, the concentration variance and the effectiveness of drug have been graphically analyzed to gain a clear understanding of endovascular delivery. The simulated results predict that with increasing values of the absorption parameter, the averaged concentrations of all drug forms do decrease. An early saturation of binding sites takes place for smaller values of the absorption parameter, and also rapid saturation of ECM binding sites occurs as compared to receptor binding sites. Results also predict the influence of surface roughness as well as asymmetry of the domain about the centerline on the distribution and retention of drug. A thorough sensitivity analysis has been carried out to determine the influence of some parameters involved.

Specific and nonspecific binding of drug eluted from a half-embedded stent in presence of atherosclerotic plaque

This study is dealt with the two-phase binding (specific and nonspecific) of drug eluted from a half- embedded drug-eluting stent in presence of atherosclerotic plaque. The specific binding due to the interaction of drug molecules with specific receptors and nonspecific binding caused by the trapping of drug in the extra-cellular matrix have been paid due attention. An idealised wall consisting of a plaque and a healthy tissue region has been considered. Moreover, a Dirichlet release condition is imposed on the strut surface. In this investigation, a two-dimensional model governing drug transport and its two-phase binding in cylindrical polar coordinate system has been solved numerically by a finite-difference method. Our simulation predicts that plaque behaves like a physical barrier in two types of the binding process and there is an inverse relationship between bound drug concentration and plaque thickness. Simulations show that a single peak profile of drug is noted when the struts are situated one-strut radius apart and as the inter-strut distance increases, the peak concentration falls and distinct peak profiles over each strut are visualised. The model also reveals that in the region downstream of a strut, the concentration of both bound drug forms in the plaque and healthy regions increases, and eventually, the saturation length of binding sites increases. Predicted results show for smaller Damköhler number, the rapid saturation of binding sites takes place and the stent having thinner strut may perform well in terms of effectiveness as well as efficacy in the stent-based delivery.