Modelling drug elution from stents: effects of reversible binding in the vascular wall and degradable polymeric matrix

Today the most popular approach for the prevention of the restenosis consists in the use of the drug eluting stents. The stent acts as a source of drug, from a coating or from a reservoir, which is transported into and through the artery wall. In this study, the behaviour of a model of a hydrophilic drug (heparin) released from a coronary stent into the arterial wall is investigated. The presence of the specific binding site action is modelled using a reversible chemical reaction that explains the prolonged presence of drug in the vascular tissue. An axi-symmetric model of a single stent strut is considered. First an advection–diffusion problem is solved using the finite element method. Then a simplified model with diffusion only in the arterial wall is compared with: (i) a model including the presence of reversible binding sites in the vascular wall and (ii) a model featuring a drug reservoir made of a degradable polymeric matrix. The results show that the inclusion of a reversible binding for the drug leads to delayed release curves and that the polymer erosion affects the drug release showing a quicker elution of the drug from the stent.     

Modelling drug elution from stents: effects of reversible binding in the vascular wall and degradable polymeric matrix

Today the most popular approach for the prevention of the restenosis consists in the use of the drug eluting stents. The stent acts as a source of drug, from a coating or from a reservoir, which is transported into and through the artery wall. In this study, the behaviour of a model of a hydrophilic drug (heparin) released from a coronary stent into the arterial wall is investigated. The presence of the specific binding site action is modelled using a reversible chemical reaction that explains the prolonged presence of drug in the vascular tissue. An axi-symmetric model of a single stent strut is considered. First an advection-diffusion problem is solved using the finite element method. Then a simplified model with diffusion only in the arterial wall is compared with: (i) a model including the presence of reversible binding sites in the vascular wall and (ii) a model featuring a drug reservoir made of a degradable polymeric matrix. The results show that the inclusion of a reversible binding for the drug leads to delayed release curves and that the polymer erosion affects the drug release showing a quicker elution of the drug from the stent.     

Geometry parameterization and multidisciplinary constrained optimization of coronary stents

Coronary stents are tubular type scaffolds that are deployed, using an inflatable balloon on a catheter, most commonly to recover the lumen size of narrowed (diseased) arterial segments. A common differentiating factor between the numerous stents used in clinical practice today is their geometric design. An ideal stent should have high radial strength to provide good arterial support post-expansion, have high flexibility for easy manoeuvrability during deployment, cause minimal injury to the artery when being expanded and, for drug eluting stents, should provide adequate drug in the arterial tissue. Often, with any stent design, these objectives are in competition such that improvement in one objective is a result of trade-off in others. This study proposes a technique to parameterize stent geometry, by varying the shape of circumferential rings and the links, and assess performance by modelling the processes of balloon expansion and drug diffusion. Finite element analysis is used to expand each stent (through balloon inflation) into contact with a representative diseased coronary artery model, followed by a drug release simulation. Also, a separate model is constructed to measure stent flexibility. Since the computational simulation time for each design is very high (approximately 24?h), a Gaussian process modelling approach is used to analyse the design space corresponding to the proposed parameterization. Four objectives to assess recoil, stress distribution, drug distribution and flexibility are set up to perform optimization studies. In particular, single objective constrained optimization problems are set up to improve the design relative to the baseline geometry—i.e. to improve one objective without compromising the others. Improvements of 8, 6 and 15% are obtained individually for stress, drug and flexibility metrics, respectively. The relative influence of the design features on each objective is quantified in terms of main effects, thereby suggesting the design features which could be altered to improve stent performance. In particular, it is shown that large values of strut width combined with smaller axial lengths of circumferential rings are optimal in terms of minimizing average stresses and maximizing drug delivery. Furthermore, it is shown that a larger amplitude of the links with minimum curved regions is desirable for improved flexibility, average stresses and drug delivery.     

Hydrodynamic effects of compliance mismatch in stented arteries

Cardiovascular diseases are the number one cause of death in the world, making the understanding of hemodynamics and development of treatment options imperative. The most common modality for treatment of occlusive coronary artery diseases is the use of stents. Stent design profoundly influences the postprocedural hemodynamic and solid mechanical environment of the stented artery. However, despite their wide acceptance, the incidence of stent late restenosis is still high (Zwart et al., 2010, "Coronary Stent Thrombosis in the Current Era: Challenges and Opportunities for Treatment," Current Treatment Options in Cardiovascular Medicine, 12(1), pp. 46-57), and it is most prevailing at the proximal and distal ends of the stent. In this work, we focus our investigation on the localized hemodynamic effects of compliance mismatch due to the presence of a stent in an artery. The compliance mismatch in a stented artery is maximized at the proximal and distal ends of the stent. Hence, it is our objective to understand and reveal the mechanism by which changes in compliance contribute to the generation of nonphysiological wall shear stress (WSS). Such adverse hemodynamic conditions could have an effect on the onset of restenosis. Three-dimensional, spatiotemporally resolved computational fluid dynamics simulations of pulsatile flow with fluid-structure interaction were carried out for a simplified coronary artery with physiologically relevant flow parameters. A model with uniform elastic modulus is used as the baseline control case. In order to study the effect of compliance variation on local hemodynamics, this baseline model is compared with models where the elastic modulus was increased by two-, five-, and tenfold in the middle of the vessel. The simulations provided detailed information regarding the recirculation zone dynamics formed during flow reversals. The results suggest that discontinuities in compliance cause critical changes in local hemodynamics, namely, altering the local pressure and velocity gradients. The change in pressure gradient at the discontinuity was as high as 90%. The corresponding changes in WSS and oscillatory shear index calculated were 9% and 15%, respectively. We demonstrate that these changes are attributed to the physical mechanism associating the pressure gradient discontinuities to the production of vorticity (vorticity flux) due to the presence of the stent. The pressure gradient discontinuities and augmented vorticity flux are affecting the wall shear stresses. As a result, this work reveals how compliance variations act to modify the near wall hemodynamics of stented arteries.     

A Thermal Analogy for Modelling Drug Elution from Cardiovascular Stents

Restriction of blood flow by the narrowing or occlusion of arteries is one of the most common presentations of cardiovascular disease. One treatment involves the introduction of a metal scaffold, or stent, designed to prevent recoil and to provide structural stability to the vessel. On the occasions that this treatment is ineffective, failure is usually associated with re-invasion of tissue. This can be prevented by local delivery of drugs which inhibit tissue growth. The drug might be delivered locally in a polymer coating on the stent. This paper develops and explores the use of a thermal analogue of the drug delivery process and the associated three-dimensional convection–diffusion equation to model the spatial and temporal distribution of drug concentration within the vessel wall. This allows the routine use of commercial finite element analysis software to investigate the dynamics of drug distribution, assist in the understanding of the treatment process and develop improved delivery systems. Two applications illustrate how the model might be used to investigate the effects of controllable or measurable parameters on the progression of the process. It is demonstrated that the geometric characteristics of the stent can have significant impact on the homogeneity of the dosing in the vessel wall.     

Modelling intravascular delivery from drug-eluting stents with biodurable coating: investigation of anisotropic vascular drug diffusivity and arterial drug distribution

In-stent restenosis occurs in coronary arteries after implantation of drug-eluting stents with non-uniform restenosis thickness distribution in the artery cross section. Knowledge of the spatio-temporal drug uptake in the arterial wall is useful for investigating restenosis growth but may often be very expensive/difficult to acquire experimentally. In this study, local delivery of a hydrophobic drug from a drug-eluting stent implanted in a coronary artery is mathematically modelled to investigate the drug release and spatio-temporal drug distribution in the arterial wall. The model integrates drug diffusion in the coating and drug diffusion with reversible binding in the arterial wall. The model is solved by the finite volume method for both high and low drug loadings relative to its solubility in the stent coating with varied isotropic–anisotropic vascular drug diffusivities. Drug release profiles in the coating are observed to depend not only on the coating drug diffusivity but also on the properties of the surrounding arterial wall. Time dependencies of the spatially averaged free- and bound-drug levels in the arterial wall on the coating and vascular drug diffusivities are discussed. Anisotropic vascular drug diffusivities result in slightly different average drug levels in the arterial wall but with very different spatial distributions. Higher circumferential vascular diffusivity results in more uniform drug loading in the upper layers and is potentially beneficial in reducing in-stent restenosis. An analytical expression is derived which can be used to determine regions in the arterial with higher free-drug concentration than bound-drug concentration.     

Analysis of prolapse in cardiovascular stents: a constitutive equation for vascular tissue and finite-element modelling

The effectiveness of a cardiovascular stent depends on many factors, such as its ability to sustain the compression applied by the vessel wall, minimal longitudinal contraction when it is expanded, and its ability to flex when navigating tortuous blood vessels. The long-term reaction of the tissue to the stent is also device dependant; in particular some designs provoke in-stent restenosis (i.e., regrowth of the occlusion around the stent). The mechanism of restenosis is thought to involve injury or damage to the vessel wall due to the high stresses generated around the stent when it expands. Because of this, the deflection of the tissue between the struts of the stent (called prolapse or "draping") has been used as a measure of the potential of a stent to cause restenosis. In this paper, uniaxial and biaxial experiments on human femoral artery and porcine aortic vascular tissue are used to develop a hyperelastic constitive model of vascular tissue suitable for implementation in finite-element analysis. To analyze prolapse, four stent designs (BeStent 2, Medtronic AVE; NIROYAL, Boston Scientific; VELOCITY, Cordis; TETRA, Guidant) were expanded in vitro to determine their repeating-unit dimensions. This geometric data was used to generate a finite element model of the vascular tissue supported within a repeating-unit of the stent. Under a pressure of 450 mm Hg (representing the radial compression of the vessel wall), maximum radial deflection of 0.253 mm, 0.279 mm, 0.348 mm and 0.48 mm were calculated for each of the four stents. Stresses in the vascular wall were highest for the VELOCITY stent. The method is proposed as a way to compare stents relative to their potential for restenosis and as a basis for a biomechanical design of a stent repeating-unit that would minimize restenosis.     

Anomalous hemodynamic effects of a self-expanding intracranial stent: Comparing in-vitro and ex-vivo models using ultra-high resolution MicroCT based CFD

Previous research on the effects of intracranial stents on arterial hemodynamics has involved computational hemodynamics (CHD) simulations applied to artificially generated stent models. In this study, accurate geometric reconstructions of in-vitro (PTFE tube) and ex-vivo (canine artery) deployed stents based on ultra-high resolution MicroCT imaging were used. The primary goal was to compare the hemodynamic effects of deployment in these two different models and to identify flow perturbations due to deployment anomalies such as stent malapposition and strut prolapse, important adverse mechanics occurring in clinical practice, but not considered in studies using idealized stent models.Ultra-high resolution MicroCT data provided detailed visualization of deployment characteristics allowing for accurate in-stent flow simulation. For stent cells that are regularly and symmetrically deployed, the near wall flow velocities and wall shear stresses were similar to previously published results derived from idealized models. In-stent hemodynamics were significantly altered by misaligned or malapposed stent cells, important effects not realistically captured in previous models. This research shows the feasibility and value of an ex-vivo stent model for MicroCT based CHD studies. It validates previous in-vitro studies and further contributes to the understanding of in-stent hemodynamics associated with adverse mechanics of self-expanding intracranial stents.     

Impact of transport and drug properties on the local pharmacology of drug-eluting stents

Drugs released from stents are driven by physiological transport forces, principally solvent-driven flow (convection) and random molecular agitation (diffusion). The relative strength of these two forces determines drug penetration and distribution in the arterial wall. Drug physicochemical factors can induce critical modulations to the primary distribution, both transiently and at steady state. Hydrophobic interactions and nonspecific binding, for example, can both result in tissue drug concentrations severalfold above administered concentration. Drug interaction with native proteins may also interfere with drug transfer at the stent-artery interface. These transport forces and tissue interactions can induce local drug concentrations even at steady state to vary by one or more orders of magnitude over the span of a few cells. To account for significant local variations in drug concentrations following stent-based delivery, rational design of vascular delivery systems requires consideration of drug distribution and tissue interactions on a local, continuum basis. Continuum analysis adapts traditional pharmacokinetics to the local environment by supplementing discrete global parameters of drug content with continuous local values of concentration, transport and binding. The interplay of these parameters with local flux conditions and drug and tissue properties defines the local drug distribution in space and over time. This type of analysis may well become increasingly relevant given the trend toward stent-based drug therapy in cardiovascular care.     

Hemodynamic analysis of a novel stent graft design with slit perforations in thoracic aortic aneurysm

Thoracic endovascular aortic repair (TEVAR) has been introduced as a less invasive approach to the treatment of thoracic aortic aneurysm (TAA). However, the effectiveness of TEVAR in the treatment of TAA is often limited due to the complex anatomy of aortic arch. Flow preservation at the three supra-aortic branches further increases the overall technical difficulty. This study proposes a novel stent graft design with slit perforations that can positively alter the hemodynamics at the aortic arch while maintaining blood flow to supra-aortic branches. We carried out a computational fluid dynamic (CFD) analysis to evaluate flow characteristics near stented aortic arch in simplified TAA models, followed by in-vitro experiments using particle image velocimetry (PIV) in a mock circulatory loop. The hemodynamics result was studied in terms of time-averaged wall shear stress (TAWSS), oscillating shear index (OSI), and endothelial cell action potential (ECAP). The results showed that the stent graft with slit perforations can reduce the disturbed flow region considerably. Furthermore, the effect of the slits on flow preservation to the supra-aortic branches was simulated and compared with experimental results. The effectiveness of the stent graft with slit perforations in preserving flow to the branches was demonstrated by both simulated and experimental results. Low TAWSS and elevated ECAP were observed in the aortic arch aneurysm after the placement of the stent graft with slits, implying the potential of thrombus formation in the aneurysm. On the other hand, the effects of the stent grafts with full-slit design and half-slit design on the shear stress did not differ significantly. The present analysis indicated that not only could the stent graft with slit perforations shield the aneurysm from rupture, but also it resulted in a favorable environment for thrombus that can contribute to the shrinkage of the aneurysm.     

A thermal analogy for modelling drug elution from cardiovascular stents

Restriction of blood flow by the narrowing or occlusion of arteries is one of the most common presentations of cardiovascular disease. One treatment involves the introduction of a metal scaffold, or stent, designed to prevent recoil and to provide structural stability to the vessel. On the occasions that this treatment is ineffective, failure is usually associated with re-invasion of tissue. This can be prevented by local delivery of drugs which inhibit tissue growth. The drug might be delivered locally in a polymer coating on the stent. This paper develops and explores the use of a thermal analogue of the drug delivery process and the associated three-dimensional convection-diffusion equation to model the spatial and temporal distribution of drug concentration within the vessel wall. This allows the routine use of commercial finite element analysis software to investigate the dynamics of drug distribution, assist in the understanding of the treatment process and develop improved delivery systems. Two applications illustrate how the model might be used to investigate the effects of controllable or measurable parameters on the progression of the process. It is demonstrated that the geometric characteristics of the stent can have significant impact on the homogeneity of the dosing in the vessel wall.     

Endovascular stent configuration affects intraluminal flow dynamics and in vitro endothelialization

Neointimal hyperplasia influenced by intravascular hemodynamics is considered partly responsible for restenosis after endovascular stenting. To evaluate the effect of stent configuration on fluid flow behavior, we visualized flow near stents, and measured the proliferation of cultured endothelial cells (ECs). A single-coil stent (coil pitch; CP = 2.5, 5, or 10 mm) was inserted into a glass tube and perfused at 30-90 ml/min, while the flow pattern was determined by particle imaging velocimetry. The reduction of the flow velocity near the wall was correlated with the decrease in the coil interval of the stent. In perfusion cultures with stents, the proliferation of ECs was influenced by the local flow velocity distribution. When a stent with a CP value of 10 mm was used, the doubling time of ECs was 30.7 h, while the doubling time was 38.5 h when the CP was 5 mm. The doubling time of ECs was shorter at sites upstream of the stent wire where the velocity was higher than downstream of the wire. In conclusion, a single-coil stent can be used to modify hemodynamic factors, suggesting that improved stent design may facilitate rapid endothelialization after stent implantation.     

A multiscale theoretical model for diffusive mass transfer in cellular biological media

An integrated methodology is developed for the theoretical analysis of solute transport and reaction in cellular biological media, such as tissues, microbial flocs, and biofilms. First, the method of local spatial averaging with a weight function is used to establish the equation which describes solute conservation at the cellular biological medium scale, starting with a continuum-based formulation of solute transport at finer spatial scales. Second, an effective-medium model is developed for the self-consistent calculation of the local diffusion coefficient in the cellular biological medium, including the effects of the structural heterogeneity of the extra-cellular space and the reversible adsorption to extra-cellular polymers. The final expression for the local effective diffusion coefficient is: D(Abeta)=lambda(beta)D(Aupsilon), where D(Aupsilon) is the diffusion coefficient in water, and lambda(beta) is a function of the composition and fundamental geometric and physicochemical system properties, including the size of solute molecules, the size of extra-cellular polymer fibers, and the mass permeability of the cell membrane. Furthermore, the analysis sheds some light on the function of the extra-cellular hydrogel as a diffusive barrier to solute molecules approaching the cell membrane, and its implications on the transport of chemotherapeutic agents within a cellular biological medium. Finally, the model predicts the qualitative trend as well as the quantitative variability of a large number of published experimental data on the diffusion coefficient of oxygen in cell-entrapping gels, microbial flocs, biofilms, and mammalian tissues.     

Effects of cellular pharmacology on drug distribution in tissues.

The efficacy of targeted therapeutics such as immunotoxins is directly related to both the extent of distribution achievable and the degree of drug internalization by individual cells in the tissue of interest. The factors that influence the tissue distribution of such drugs include drug transport; receptor/drug binding; and cellular pharmacology, the processing and routing of the drug within cells. To examine the importance of cellular pharmacology, previously treated only superficially, we have developed a mathematical model for drug transport in tissues that includes drug and receptor internalization, recycling, and degradation, as well as drug diffusion in the extracellular space and binding to cell surface receptors. We have applied this "cellular pharmacology model" to a model drug/cell system, specifically, transferrin and the well-defined transferrin cycle in CHO cells. We compare simulation results to models with extracellular diffusion only or diffusion with binding to cell surface receptors and present a parameter sensitivity analysis. The comparison of models illustrates that inclusion of intracellular trafficking significantly increases the total transferrin concentration throughout much of the tissue while decreasing the penetration depth. Increasing receptor affinity or tissue receptor density reduces permeation of extracellular drug while increasing the peak value of the intracellular drug concentration, resulting in "internal trapping" of transferrin near the source; this could account for heterogeneity of drug distributions observed in experimental systems. Other results indicate that the degree of drug internalization is not predicted by the total drug profile. Hence, when intracellular drug is required for a therapeutic effect, the optimal treatment may not result from conditions that produce the maximal total drug distribution. Examination of models that include cellular pharmacology may help guide rational drug design and provide useful information for whole body pharmacokinetic studies.     

Optimization of cardiovascular stent design using computational fluid dynamics

Coronary stent design affects the spatial distribution of wall shear stress (WSS), which can influence the progression of endothelialization, neointimal hyperplasia, and restenosis. Previous computational fluid dynamics (CFD) studies have only examined a small number of possible geometries to identify stent designs that reduce alterations in near-wall hemodynamics. Based on a previously described framework for optimizing cardiovascular geometries, we developed a methodology that couples CFD and three-dimensional shape-optimization for use in stent design. The optimization procedure was fully-automated, such that solid model construction, anisotropic mesh generation, CFD simulation, and WSS quantification did not require user intervention. We applied the method to determine the optimal number of circumferentially repeating stent cells (N(C)) for slotted-tube stents with various diameters and intrastrut areas. Optimal stent designs were defined as those minimizing the area of low intrastrut time-averaged WSS. Interestingly, we determined that the optimal value of N(C) was dependent on the intrastrut angle with respect to the primary flow direction. Further investigation indicated that stent designs with an intrastrut angle of approximately 40 deg minimized the area of low time-averaged WSS regardless of vessel size or intrastrut area. Future application of this optimization method to commercially available stent designs may lead to stents with superior hemodynamic performance and the potential for improved clinical outcomes.     

Influence of stent-induced vessel deformation on hemodynamic feature of bloodstream inside ICA aneurysms

One of the effective treatment options for intracranial aneurysms is stent-assisted coiling. Though, previous works have demonstrated that stent usage would result in the deformation of the local vasculature. The effect of simple stent on the blood hemodynamics is still uncertain. In this work, hemodynamic features of the blood stream on four different ICA aneurysm with/without interventional are investigated. To estimate the relative impacts of vessel deformation, four distinctive ICA aneurysm is simulated by the one-way FSI technique. Four hemodynamic factors of aneurysm blood velocity, wall pressure and WSS are compared in the peak systolic stage to disclose the impact of defamation by the stent in two conditions. The stent usage would decrease almost all of the mentioned parameters, except for OSI. Stenting reduces neck inflow rate, while the effect of interventional was not consistent among the aneurysms. The deformation of an aneurysm has a strong influence on the hemodynamics of an aneurysm. This outcome is ignored by most of the preceding investigations, which focused on the pre-interventional state for studying the relationship between hemodynamics and stents. Present results show that the application of stent without coiling would improve most hemodynamic factors, especially when the deformation of the aneurysm is high enough.

     

Fitting density functions and diffusion tensors to three-dimensional drug transport within brain tissue

In the research on drug transport within brain tissue, typical data sets consist of a collection of volume integrals from various locations in the brain, with each datum measuring the total amount of drug in a sample of tissue. These samples contain medically useful information about the extent and orientation of the distribution of the drug, but extracting this information requires some mathematical analysis. The method of analysis presented here performs a nonlinear regression on the data to fit a multivariate density function (for example, a Gaussian density) to model such transport processes as diffusion and dispersion. The principal components of that density then characterize the size and shape of the distribution of drug.     

The release kinetics of drug eluting stents containing sirolimus as coated drug: Role of release media

Prior understanding on the in vitro release profile of the drug from drug eluting devices such as stent (DES) is crucial in designing and optimizing the drug embedded coating or matrices. In fact, assessing in vitro release profile is a mandatory requirement prior to the clinical evaluation of DES. The in vitro release is also employed to estimate parameters such as T1/2. The release profile largely depends on the release medium selected for the studies. Normally PBS with a pH of 7.4 is used for assessing the release kinetics of the drug. Often drug undergoes irreversible changes such as hydrolysis in PBS leading to erroneous assessment of the release profile. This is particularly true in the case of sirolimus, one of the widely used drugs in various applications. We studied the influence of various media on the release profile of sirolimus from DES. The data generated suggested that a release medium consisting of 9:1 (v/v) of normal saline and isopropanol is a most suitable one for assessing in vitro the release kinetics of sirolimus from DES.     

Stent expansion in curved vessel and their interactions: a finite element analysis

Coronary restenosis after angioplasty has been reduced by stenting procedure, but in-stent restenosis (ISR) has not been eliminated yet, especially in tortuous vessels. In this paper, we proposed a finite element method (FEM) to study the expansion of a stent in a curved vessel (the CV model) and their interactions. A model of the same stent in a straight vessel (the SV model) was also studied and mechanical parameters of both models were researched and compared, including final lumen area, tissue prolapse between stent struts and stress distribution. Results show that in the CV model, the vessel was straightened by stenting and a hinge effect can be observed at extremes of the stent. The maximum tissue prolapse of the CV model was more severe (0.079 mm) than the SV model (0.048 mm); and the minimum lumen area of the CV was decreased (6.10 mm(2)), compared to that of the SV model (6.28 mm(2)). Tissue stresses of the highest level were concentrated in the inner curvature of the CV model. The simulations offered some explanations for the clinical results of ISR in curved vessels and gave design suggestions of the stent and balloon for tortuous vessels. This FEM provides a tool to study mechanisms of stents in curved vessels and can improve new stent designs especially for tortuous vessels.