Combining mathematical modelling with in vitro experiments to predict in vivo drug-eluting stent performance

3D-printed, citrate-based bioresorbable vascular scaffolds for coronary artery angioplasty

Fully bioresorbable vascular scaffolds (BVSs) aim to overcome the limitations of metallic drug-eluting stents (DESs). However, polymer-based BVSs, such as Abbott's Absorb, the only US FDA-approved BVS, have had limited use due to increased strut thickness (157 μm for Absorb), exacerbated tissue inflammation, and increased risk of major cardiac events leading to inferior clinical performance when compared to metallic DESs. Herein we report the development of a drug-eluting BVS (DE-BVS) through the innovative use of a photopolymerizable, citrate-based biomaterial and a high-precision additive manufacturing process. BVS with a clinically relevant strut thickness of 62 μm can be produced in a high-throughput manner, i.e. one BVS per minute, and controlled release of the anti-restenosis drug everolimus can be achieved by engineering the structure of polymer coatings to fabricate drug-eluting BVS. We achieved the successful deployment of BVSs and DE-BVSs in swine coronary arteries using a custom-built balloon catheter and BVS delivery system and confirmed BVS safety and efficacy regarding maintenance of vessel patency for 28 days, observing an inflammation profile for BVS and DE-BVS that was comparable to the commercial XIENCE™ DES (Abbott Vascular).

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.

Modeling Dual Drug Delivery from Eluting Stents: The Influence of Non-Linear Binding Competition and Non-Uniform Drug Loading

OBJECTIVE

There is increasing interest in simultaneous endovascular delivery of more than one drug from a drug-loaded stent into a diseased artery. There may be an opportunity to obtain a therapeutically desirable uptake profile of the two drugs over time by appropriate design of the initial drug distribution in the stent. Due to the non-linear, coupled nature of diffusion and reversible specific/non-specific binding of both drugs as well as competition between the drugs for a fixed binding site density, a comprehensive numerical investigation of this problem is critically needed.

METHODS

This paper presents numerical computation of dual drug delivery in a stent-artery system, accounting for diffusion as well as specific and non-specific reversible binding. The governing differential equations are discretized in space, followed by integration over time using a stiff numerical solver. Three different cases of initial dual drug distribution are considered.

RESULTS

For the particular case of sirolimus and paclitaxel, results show that competition for a limited non-specific binding site density and the significant difference in the forward/backward reaction coefficients play a key role in determining the nature of drug uptake. The nature of initial distribution of the two drugs in the stent is also found to influence the binding process, which can potentially be used to engineer a desirable dual drug uptake profile.

CONCLUSIONS

These results help improve the fundamental understanding of endovascular dual drug delivery. In addition, the numerical technique and results presented here may be helpful for designing and optimizing other drug delivery problems as well.

How does the Nature of an Excipient and an Atheroma Influence Drug-Coated Balloon Therapy?

The advent of drug-eluting stents and drug-coated balloons have significantly improved the clinical outcome of patients with vascular occlusions. However, ischemic vascular disease remains the most common cause of death worldwide. Improving the current treatment modalities demands a better understanding of the processes which govern drug uptake and retention in blood vessels. In this study, we evaluated the influence of urea and butyryl-trihexyl citrate, as excipients, on the efficacy of drug-coated balloon therapy. An integrated approach, utilizing both in-vitro and in-silico methods, was used to quantify the tracking loss, vessel adhesion, drug release, uptake, and distribution associated with the treatment. Moreover, a parametric study was used to evaluate the potential influence of different types of lesions on drug-coated balloon therapy. Despite the significantly higher tracking loss (urea: 35.5% vs. butyryl-trihexyl citrate: 8.13%) observed in the urea-based balloons, the drug uptake was almost two times greater than with its hydrophobic counterpart. Non-calcified lesions were found to delay the transmural propagation of sirolimus while calcification was shown to limit the retentive potential of lesions. Ultimately this study helps to elucidate how different excipients and types of lesions may influence the efficacy of drug-coated balloon therapy.

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.

Anti-coagulation and anti-hyperplasia coating for retrievable vena cava filters by electrospraying and their performance in vivo

A novel drug eluting retrievable vena cava filter (RVCF) with a heparin-modified poly(ε-caprolactone) (hPCL) coating containing rapamycin was prepared by electrospraying. The in vitro drug release pattern showed that the encapsulated rapamycin in the coating can be sustainably released within one month, whereas activated partial thromboplastin time (APTT) and in vitro cell culture showed that the drug eluting RVCF can effectively extend blood clotting time and inhibit smooth muscle cell (SMC) and endothelial cell (EC) proliferation, respectively. The as-prepared drug eluting RVCF and corresponding commercial RVCF were implanted into the vena cava of sheep. The retrieval operation at a predetermined time point showed that the drug eluting RVCF had a much higher retrieval rate than the commercial RVCF. Comprehensive investigations, including histological, immunohistological and immunofluorescence analyses, on explanted veins were carried out. The results demonstrated that the as-prepared RVCF possessed excellent antihyperplasia properties in vivo, significantly improving the retrieval rate and extending the in vivo dwelling time in sheep. Consequently, the drug eluting RVCF has promising potential for application in the clinic to improve RVCF retrieval rates.

How Does Fluid Flow Influence Drug Release from Drug Filled Implants?

Drug-filled implants (DFIs) have emerged as an innovative approach to control the delivery of drugs. These devices contain the drug within the structure of the implant itself and avoid the need to include additional drug carrier materials such as a polymers, which are often associated with inflammation and delayed healing/tissue regeneration at the implant site. One common feature of in vitro experiments to generate drug release profiles is stirring or agitation of the release medium. However, the influence of the resulting fluid flow on the rate of drug release from DFIs has yet to be quantified. In this paper we consider two DFIs, which although similar in shape and size, employ different strategies to control the release of drug: a porous pin with pores on the order of μm and a pin drilled with orifices of the order of mm. We develop a multiphysics mathematical model of drug release from these DFIs, subject to fluid flow induced through stirring and show that fluid flow greatly influences the drug release profile for the orifice pin, but that the porous pin drug release profile is relatively insensitive to flow. We demonstrate that drug release from the porous pin may adequately be described through a simplified radial 1D dissolution-diffusion model, while a 3D dissolution-advection-diffusion model is required to describe drug release from the orifice pin. A sensitivity analysis reveals that that the balance of reaction-advection-diffusion in terms of key nondimensional numbers governs the overall drug release. Our findings potentially have important implications in terms of devising the most relevant experimental protocol for quantifying drug release from DFIs.

3D modelling of drug-coated balloons for the treatment of calcified superficial femoral arteries

Background/Objectives

Drug-coated balloon therapy for diseased superficial femoral arteries remains controversial. Despite its clinical relevance, only a few computational studies based on simplistic two-dimensional models have been proposed to investigate this endovascular therapy to date. This work addresses the aforementioned limitation by analyzing the drug transport and kinetics occurring during drug-coated balloon deployment in a three-dimensional geometry.

Methods

An idealized three-dimensional model of a superficial femoral artery presenting with a calcific plaque and treated with a drug-coated balloon was created to perform transient mass transport simulations. To account for the transport of drug (i.e. paclitaxel) released by the device, a diffusion-reaction equation was implemented by describing the drug bound to specific intracellular receptors through a non-linear, reversible reaction. The following features concerning procedural aspects, pathologies and modelling assumptions were investigated: (i) balloon application time (60–180 seconds); (ii) vessel wall composition (healthy vs. calcified wall); (iii) sequential balloon application; and (iv) drug wash-out by the blood stream vs. coating retention, modeled as exponential decay.

Results

The balloon inflation time impacted both the free and specifically-bound drug concentrations in the vessel wall. The vessel wall composition highly affected the drug concentrations. In particular, the specifically-bound drug concentration was four orders of magnitude lower in the calcific compared with healthy vessel wall portions, primarily as a result of reduced drug diffusion. The sequential application of two drug-coated balloons led to modest differences (~15%) in drug concentration immediately after inflation, which became negligible within 10 minutes. The retention of the balloon coating increased the drug concentration in the vessel wall fourfold.

Conclusions

The overall findings suggest that paclitaxel kinetics may be affected not only by the geometrical and compositional features of the vessel treated with the drug-coated balloon, but also by balloon design characteristics and procedural aspects that should be carefully considered.

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.

Mechano-chemo-biological Computational Models for Arteries in Health, Disease and Healing: From Tissue Remodelling to Drug-eluting Devices

This review aims to highlight urgent priorities for the computational biomechanics community in the framework of mechano-chemo-biological models. Recent approaches, promising directions and open challenges on the computational modelling of arterial tissues in health and disease are introduced and investigated, together with in silico approaches for the analysis of drug-eluting stents that promote pharmacological-induced healing. The paper addresses a number of chemo-biological phenomena that are generally neglected in biomechanical engineering models but are most likely instrumental for the onset and the progression of arterial diseases. An interdisciplinary effort is thus encouraged for providing the tools for an effective in silico insight into medical problems. An integrated mechano-chemo-biological perspective is believed to be a fundamental missing piece for crossing the bridge between computational engineering and life sciences, and for bringing computational biomechanics into medical research and clinical practice.

Conditioning the microenvironment for soft tissue regeneration in a cell free scaffold

The use of cell-free scaffolds for the regeneration of clinically relevant volumes of soft tissue has been challenged, particularly in the case of synthetic biomaterials, by the difficulty of reconciling the manufacturing and biological performance requirements. Here, we investigated in vivo the importance of biomechanical and biochemical cues for conditioning the 3D regenerative microenvironment towards soft tissue formation. In particular, we evaluated the adipogenesis changes related to 3D mechanical properties by creating a gradient of 3D microenvironments with different stiffnesses using 3D Poly(Urethane-Ester-ether) PUEt scaffolds. Our results showed a significant increase in adipose tissue proportions while decreasing the stiffness of the 3D mechanical microenvironment. This mechanical conditioning effect was also compared with biochemical manipulation by loading extracellular matrices (ECMs) with a PPAR-γ activating molecule. Notably, results showed mechanical and biochemical conditioning equivalency in promoting adipose tissue formation in the conditions tested, suggesting that adequate mechanical signaling could be sufficient to boost adipogenesis by influencing tissue remodeling. Overall, this work could open a new avenue in the design of synthetic 3D scaffolds for microenvironment conditioning towards the regeneration of large volumes of soft and adipose tissue, with practical and direct implications in reconstructive and cosmetic surgery.

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.

Influence of vessel curvature and plaque composition on drug transport in the arterial wall following drug-eluting stent implantation

In the last decade, many computational models have been developed to describe the transport of drug eluted from stents and the subsequent uptake into arterial tissue. Each of these models has its own set of limitations: for example, models typically employ simplified stent and arterial geometries, some models assume a homogeneous arterial wall, and others neglect the influence of blood flow and plasma filtration on the drug transport process. In this study, we focus on two common limitations. Specifically, we provide a comprehensive investigation of the influence of arterial curvature and plaque composition on drug transport in the arterial wall following drug-eluting stent implantation. The arterial wall is considered as a three-layered structure including the subendothelial space, the media and the adventitia, with porous membranes separating them (endothelium, internal and external elastic lamina). Blood flow is modelled by the Navier-Stokes equations, while Darcy's law is used to calculate plasma filtration through the porous layers. Our findings demonstrate that arterial curvature and plaque composition have important influences on the spatiotemporal distribution of drug, with potential implications in terms of effectiveness of the treatment. Since the majority of computational models tend to neglect these features, these models are likely to be under- or over-estimating drug uptake and redistribution in arterial tissue.

Effect of aging on heat transfer, fluid flow and drug transport in anterior human eye: A computational study

There are a lot of geometrical and morphological changes that happen in the human eye with age. Primary open-angle glaucoma, which is caused by the increase in intraocular pressure inside the anterior chamber of the eye is also associated with the physiological aging of the eye. Therefore, it is crucial to understand the effects of aging on drug delivery in the human eye when applied topically. Consequently, a numerical model of topical drug delivery for an aging human eye has been developed using commercial software COMSOL Multiphysics in the current study. Three different age groups (young, middle and old) have been considered and the changes in geometrical and tissue properties of different domains of the eye with age have been included in the numerical model. The effect of aging on heat transfer, aqueous humor flow, intraocular pressure and drug concentration in different domains and orientations of the eye have been investigated. Additionally, an attempt has been made to predict the best class of anti-glaucomatic treatment in silico that should be preferred to treat primary open-angle glaucoma effectively. Results illustrate that there is a decrease in the average corneal temperature and an increase in the temperature deviation across the cornea with age. Further, there is a decrease in the aqueous humor flow magnitude in the anterior chamber of the eye and an increase in intraocular pressure in the anterior chamber of older age groups, which leads to primary open-angle glaucoma. The reduced aqueous humor flow leads to increased drug concentration in the anterior chamber as well as iris and reduced drug concentration in the trabecular mesh of the older age groups, thereby affecting the treatment efficacy. Additionally, our simulated results demonstrate that anti-glaucomatic treatments should be more focused on treating the trabecular mesh rather than the ciliary body of the eye.

Assessing the influence of atherosclerosis on drug coated balloon therapy using computational modelling

Interventional therapies such as drug-eluting stents (DES) and drug-coated balloons (DCB) have significantly improved the clinical outcomes of patients with coronary occlusions in recent years. Despite this marked improvement, ischemic cardiovascular disease remains the most common cause of death worldwide. To address this, research efforts are focused on improving the safety and efficacy of the next generation of these devices. However, current experimental methods are unable to account for the influence of atherosclerotic lesions on drug uptake and retention. Therefore, in this study, we used an integrated approach utilizing both in vitro and in silico methods to assess the performance of DCB therapy. This approach was validated against existing in vivo results before being used to numerically estimate the effect of the atheroma. A bolus release of sirolimus was observed with our coating matrix. This, coupled with the rapid saturation of specific and non-specific binding sites observed in our study, indicated that increasing the therapeutic dose coated onto the balloons might not necessarily result in greater uptake and/or retention. Additionally, our findings alluded to an optimal exposure time, dependent on the coating matrix, for the DCBs to be expanded against the vessel. Moreover, our findings suggest that a biphasic drug release profile might be beneficial for establishing and maintaining the saturation of bindings sites within severely occluded vessels. Ultimately, we have demonstrated that computational methods may be capable of assessing the efficacy of DCB therapy as well as predict the influence of atherosclerotic lesions on said efficacy.

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

An attempt is made in the present investigation to develop a computational model for the purpose of studying the effect of interstitial flow in the porous media on the distribution of drug eluted from a half-embedded drug-eluting stent and its retention in the presence of two-layered coating of the stent. The transport of free drug inside the coatings is considered as an unsteady diffusion process while that in the tissue as an unsteady convection-diffusion-reaction process. The bound drug is governed by an unsteady reaction process only. Immersed boundary method (IBM) in the staggered grid formulation, popularly known as marker and cell (MAC) method, has been leveraged to tackle numerically the governing equations. This model highlights the benefits of consideration of two-layered coating and does predict underlying mechanism for better efficacy by tweaking the kinetics parameters. Comparisons are also made with the results available for stent-based delivery.

Development of a Bioactive Polymeric Drug Eluting Coronary Stent Coating Using Electrospraying

Drug-eluting stents are now routinely used in the treatment of acute coronary syndromes caused by coronary artery disease. Whilst the sustained release of anti-proliferative drugs from these devices has greatly reduced the need for repeat revascularisation procedures, this approach is not suitable for all patients and appears to delay regrowth of the endothelium, necessitating the use of prolonged dual anti-platelet therapy. Although the development of more advanced stent platforms and drug coatings has produced modest improvements in performance, these devices have not fully addressed the limitations experienced with their first-generation counterparts. In the present study, we developed a novel stent coating that provides controlled sirolimus release from a bioactive polymer (accelerate™ AT) that has previously been shown to support endothelial cell growth in vitro. A bespoke electrospray deposition process provided control over the coating thickness, surface roughness, drug load, and release kinetics. The resultant optimised coating combines rapid release of an anti-proliferative agent from a bioactive polymer coating that promotes re-endothelialisation, thereby offering potential protection against in-stent restenosis and thrombosis. This novel, dual-action coating therefore has significant therapeutic potential, with the enhanced control of drug load and release kinetics offered by electrospray deposition also opening up opportunities for more personalised treatment approaches. Further development and evaluation of these technologies in vitro and in vivo is therefore warranted.