Home
JournalsCollections
For Authors For Reviewers For Editorial Board Members
Article Processing Charges Open Access
Ethics Advertising Policy
Editorial Policy Resource Center
Company Information Contact Us
OPEN ACCESS

A Computational Approach for Stent Elution Rate Determined Specific Drug Binding and Receptor-mediated Effects in Arterial Tissue

  • Ramprosad Saha*
Journal of Exploratory Research in Pharmacology   2018;3(4):105-118

doi: 10.14218/JERP.2018.00018

Received:

Revised:

Accepted:

Published online:

 Author information

Citation: Saha R. A Computational Approach for Stent Elution Rate Determined Specific Drug Binding and Receptor-mediated Effects in Arterial Tissue. J Explor Res Pharmacol. 2018;3(4):105-118. doi: 10.14218/JERP.2018.00018.

Abstract

Background and objective

In the present article the effects of drug binding (both specific and nonspecific) in the porous arterial wall following stent-based drug delivery from drug-eluting stents (DESs) are investigated. A three-phase (free, extracellular matrix-bound, and specific receptor-bound) second-order nonlinear saturable reversible binding model is considered in order to describe the binding process with the constituents of the porous arterial wall. Although, there are currently some precise forms of a drug binding model in the arterial tissue in the literature, analyzed by various authors. The specific interest in this present context is in assessing to what extent modelling of specific and nonspecific binding within a single-layered homogeneous porous arterial wall is possible. A novel axi-symmetric model of drug delivery from three stent struts has been developed and is presented.

Methods

The governing equations of motion together with the physiologically realistic boundary conditions are tackled numerically by an explicit finite-difference scheme in staggered grids.

Results

Results include the influence of the significant model parameters, such as Peclet numbers (PeT, Pe1 and Pe2), Damköhler numbers (Da1 and Da2) and time-dependent release kinetics as well as constant release kinetics. Consistency of the proposed approach is shown graphically.

Conclusions

As the porosity (εw) increases, the effective as well as the true diffusivity increases, which eventually leads to expedition of the diffusion process. In a porous media, diffusion takes place in confined tortuous pores and its progression is impeded as the tortuosity increases. The present simulation also demonstrates a decrease in the mean concentration of free as well as extracellular matrix-bound and SR-bound drug with increasing tortuosity. The present observation may be justified in the sense that as the tortuosity increases so too does the effective distance over which diffusion has to take place (i.e. the progression of diffusion is impeded, which eventually lowers the mean concentration of all drug forms).

Keywords

Drug-eluting stent, Free drug, ECM-bound drug, SR-bound drug, Convection-diffusion-reaction, Peclet number, Damköhler number

Introduction

Arterial stents are medical devices that have revolutionized the treatment of coronary artery disease. They serve to reopen the occluded vessel that has become narrowed as a result of atherosclerosis. Atherosclerosis is a common degenerative disease that affects coronary, carotid and other peripheral arteries in the body. Now, it is standard to use the types of stents that gradually release anti-proliferative/anti-inflammatory drugs into the arterial wall to inhibit the cell proliferation that causes the development of restenosis (re-narrowing). The bare metal stents (BMSs), while revolutionary at the time, were soon rendered unsatisfactory due to their inability to prevent in-stent restenosis. The next wave of drug-eluting stents (DESs) consists of a supporting metallic wired scaffold coated with a polymer film that encapsulates the therapeutic drug aimed at preventing hyperplasia of the smooth muscle cells (SMCs) that is responsible for restenosis.

In order to control the release rate, the coating may include a rate-limiting barrier. To ensure effective performance of DES, both the stent geometry and coating design need to be optimized so that therapeutic levels of drug are delivered to the arterial wall for the required period of time.1,2 The success of an anti-proliferative drug therapy from DES depends on the amount of drug eluted from the stent, accumulation of the drug, and drug binding to cells in the arterial wall. Even though DES are now the primary choice of percutaneous coronary interventions (PCIs) for millions of patients, many questions still remain unanswered regarding their longevity and safety.

Although it was not the intention of this article to provide a review of the previously published models, it will describe some experimental studies which have been carried out in the recent past in order to quantify the capability of this device to reduce the in-stent restenosis rate after stent implantation.3–5 The behavior of heparin in explanted arteries allows for the presence of binding site changes along the transmural direction, being higher in the endothelium and lower in the adventitia, as studied by Lovich et al.3 In the experimental study by Migliavacca et al.,6 the pattern of DES drug release in the vascular wall was considered using a single species approach along with a partition coefficient approach to relate the free and the bound drug concentrations. Borgi et al.7 focused on the inclusion of reversible binding leading to delayed release and that the erosion of polymer affects the drug release from a single strut. Horner et al.8 appear to be one of the first groups to provide a three-dimensional reaction-diffusion-convection process of a two-species drug delivery model, including reversible binding sites in a realistic geometry; their model predicted that a single species drug delivery model cannot accurately predict the distribution of the bound drug. Ferreira et al.9 considered a series of nonlinear binding models to describe the degradation of a poly-L-lactic acid stent coating into lactic acid and oligomers.

Although a large number of mathematical models are available to describe drug transport and its binding to arterial tissue sites, only a few (v.i.z. Tzafriri et al.10) considered a nonlinear saturable binding model; later, Bozsak et al.11,12 also considered this. The model from the latter included two phases of drug in the tissue: free and bound. However, it is well established that in addition to binding to specific receptors (SR), there is also the occurrence of nonspecific binding caused by association of drug with membrane constituents or by trapping of the drug in the extracellular medium.13 Most recently, Tzafriri et al.14 and McGinty et al.15 included two equations for drug binding in arterial tissue, namely one for specific binding to receptors and another for nonspecific binding to general extracellular matrix (ECM) sites. Thus, it appears that there are three phases in the tissue, comprising two bound (SR and ECM) and one free.

The main aim of this investigation was to extend the aforementioned work14,15 with a two-species model of specific and nonspecific saturable binding in the arterial wall at different phases following drug transport eluted from three struts, where the transport of free drug is governed by a convection-diffusion-reaction process and that of bound drugs (SR and ECM) by a reaction process only. A simple time-dependent release kinetics is implemented on the surface of the struts.16 The transport of drugs within the arterial tissue is controlled by arterial properties like porosity and tortuosity. At the time of implantation of an endovascular DES, its major impact is on the structure of the arterial wall, which eventually influences the overall rates of diffusion through tissues. The effective diffusivity of a porous wall is supposed to depend on two factors, such as porosity and tortuosity—these parameters regulate the free diffusivity of the drug eluted from struts.17 As such, the present study also deals with the effects of porosity and tortuosity on the diffusivity of drug.

Materials and methods

Please refer to Table 1 for nomenclature.

Table 1

Nomenclature

Nomenclature
LDimensionless length of the artery
csInitial drug concentration on the stent
VwallTransmural filtration velocity
rDimensional radial coordinate
rDimensionless radial coordinate
zDimensional axial coordinate
zDimensionless axial coordinate
tDimensional time
tDimensionless time
cfDimensional concentration of free drug
cfDimensionless concentration of free drug
cbECMDimensional concentration of ECM-bound drug
cbECMDimensionless concentration of ECM-bound drug
cbSRDimensional concentration of SR-bound drug
cbSRDimensionless concentration of SR-bound drug
cbECMmaxECM binding site density
cbSRmaxReceptor density
BMTotal tissue binding capacity
KdECMECM binding on-rate
konSRReceptor binding on-rate
KdEquilibrium association constant
KdECMECM dissociation constant
KdSRReceptor dissociation constant
DfreeFree drug diffusivity
DeffEffective diffusivity of free drug
DTTrue drug diffusivity of arterial wall
DTECMTrue drug diffusivity of ECM sites
DTSRTrue drug diffusivity of SR sites
PeTNon-dimensional Peclet number
Pe1Non-dimensional Peclet number
Pe2Non-dimensional Peclet number
Da1Dimensionless Damköhler number
Da2Dimensionless Damköhler number

Geometric model

The computational domain is comprised of a long axial section of length L and the wall thickness is taken to be 10 times the strut height (δ). The axis of symmetry is taken along the centerline of the artery (cf. Fig. 1). The volume-averaged molar concentration of free drug is denoted by cf, the volume-averaged molar concentration of bound drug that is bound to nonspecific general ECM sites in the tissue is referred to as ECM-bound drug and denoted by cbECM, and the volume-averaged molar concentration of bound drug that is bound to specific receptors is referred to as SR-bound drug and denoted by cbSR. The inter-conversion of drug between the unbound plasma phase and the bound phase of tissue binding sites is controlled by a second-order nonlinear reversible saturable chemical reaction. The transport of free drug eluted from struts is governed by unsteady convection-diffusion-reaction process (Eq. 1),18–20 the ECM-bound drug is represented by unsteady reaction process (Eq. 2),21 and that of the SR-bound drug is represented by unsteady reaction process (Eq. 3).22 Symmetry boundary conditions for both the free and both of the bound drugs are applied at the proximal (Γti) and the distal (Γto) walls (Eq. 4).23 Impermeable boundary condition for both of the bound drugs is assumed at the perivascular wall (Γtp), lumen-tissue (Γbt) and strut-tissue (Γst) interfaces (Eq. 5). For the free drug, a perfectly sink condition is imposed at the perivascular end (Eq. 6). Since a proper boundary condition for the free drug at lumen-tissue interface (Γbt) is not readily apparent, two opposing extremes consider either that flowing blood is extremely efficient at washing out mural-adhered drug, modelled as a zero-concentration interface condition, or mural-adhered drug is insensitive to flowing blood, modelled as a zero-flux boundary condition (Eq. 7).24 Instead of modelling a uniform release of drug from struts, a simple time-dependent release kinetics (Eq. 8)25,26 is assumed.

Schematic of the computational model used for the study.
Fig. 1  Schematic of the computational model used for the study.

Governing equations and boundary conditions

Therefore, the governing equations of the drug transport of free, ECM-bound (saturable binding to general ECM sites) and SR-bound (saturable binding to specific receptors) in the arterial wall are respectively represented in a two-dimensional (2D) Cartesian coordinate system in the following manner:

c¯ft¯+γwεw(Vwallc¯f)r¯=DT[2c¯fr¯2+1r¯c¯fr¯+2c¯fz¯2]c¯bECMt¯c¯bSRt¯,
c¯bECMt¯=[konECMc¯f(c¯bECMmaxc¯bECM)konECMKdECMc¯bECM],
c¯bSRt¯=[konSRc¯f(c¯bSRmaxc¯bSR)konSRKdSRc¯bSR],

Where t denotes time since stent implantation, r is the distance from the intima, cf is the molar concentration of free drug per unit tissue volume, cbECM and cbSR are the molar concentrations of ECM- and receptor-bound drug respectively, cbECMmax and cbSRmax denote the local molar concentration of ECM and receptor drug binding sites respectively, konECM and KonSR are the respective binding on-rate constants, KdECM and KdSR are the respective equilibrium dissociation constants, and Vwall is its transmural convective velocity. Here, DT is the transmural true diffusivity of the drug which can be written as 17,27

DT=[1+BMKd]×Deff,
where
Deff=εwτw×Dfree.

Here, εw and τw are the porosity and the tortuosity of the wall material respectively, Dfree and Deff are the coefficients of free and effective diffusivity respectively, BM and Kd are the net tissue binding capacity and equilibrium association constant respectively.

The symmetry boundary conditions (free, ECM-bound and SR-bound drug) are imposed on the proximal (Γti) and the distal (Γto) walls as

c¯fz¯=0=c¯bECMz¯=c¯bSRz¯onΓtiandΓto,

The impermeable boundary condition for both of the bound drugs (ECM-bound and SR-bound) is assumed at the perivascular wall (Γtp), lumen-tissue (Γbt) and strut-tissue (Γst) interfaces as

c¯bECMr¯=0=c¯bSRr¯onΓtl(=ΓbtΓst)andΓtp,

At the perivascular wall (Γtp), a perfectly sink condition4,11 is imposed for the free drug (cf) as

c¯f=0onΓtp,

At the lumen-tissue interface (Γbt) a proper boundary condition for the free drug is not readily apparent. Considering two opposing extremes, either the flowing blood is extremely efficient at washing out mural-adhered drug, modelled as a zero-concentration interface condition, or mural-adhered drug is insensitive to flowing blood, modelled as zero-flux boundary condition as

c¯f=0orc¯fr¯=0onΓbt,

As the dimension of stent struts is thin with respect to the arterial wall thickness, so their actual geometry is neglected and approximate the drug eluting stent using an equivalent phantom surface that elutes a defined drug load to the arterial wall (Fig. 1). Elution from the phantom surface is modelled by a simple time-dependent release kinetics as follows

c¯f=csexp(λ¯t¯);t¯0onΓst,
where cs is the initial drug concentration on the stent, and λ is the stent drug decay rate.

All the variables and parameters are now made dimensionless to obtain well-behaved computations in the following manner:

r=r¯δ,z=z¯δ,t=t¯.Vwallδ,cf=c¯fcs,cbECM=c¯bECMcbECMmax,cbSR=c¯bSRcbSRmax

Under these assumptions, the above equations (18) take their respective non-dimensional forms as follows:

cft+γwεwcfr=1PeT[2cfr2+1rcfr+2cfz2]1αcbECMt1βcbSRt,
cbECMt=αDa1Pe1[cf(1cbECM)α1cbECM],
cbSRt=βDa2Pe2[cf(1cbSR)β1cbSR],
cfz=0=cbECMz=cbSRzonΓtiandΓto,
cbECMr=0=cbSRronΓtl(=ΓbtΓst)andΓtp,
cf=0onΓtp,
cf=0orcfr=0onΓbt,
cf=csexp(λt);t0onΓst,

where the Peclet numbers (PeT, Pe1 and Pe2), the Damköhler numbers (Da1 and Da2), the scaling parameters (α,α1, β and β1) and dimensionless drug release rate (λ) are defined respectively as:

PeT=Vwall.δDT,Pei=Vwall.δDT(i),Dai=kon(i).cb(i)max.δ2DT(i)(i=1,2;(1)ECM and (2)SR)α=cscbECMmax,α1=KdECMcs,β=cscbSRmax,β1=KdSRcs,λ=λ¯Vwall,
where KdECM(krECM/konECM) and KdSR(krSR/konSR) are, respectively, the equilibrium dissociation constants of the ECM-bound drug and the SR-bound drug, krECM and krSR are the dissociation (backward) rate constant of ECM and receptor drug binding sites in the arterial tissue respectively. Here, DTECM and DTSR are, respectively, the true diffusivities of the ECM-bound drug and the SR-bound drug.

Solution procedure

The governing equations (911) representing the transport of free, ECM-bound and SR-bound drugs were solved numerically using a finite-difference scheme by applying the time-dependent release kinetics and boundary conditions (12–16). For this type of grid alignment, the concentrations for free, ECM-bound and SR-bound drugs were calculated at the cell centers (Fig. 2). The discretization of the time derivative term was based on the first order accurate two-level forward time-differencing formula, while the convective term in the equations is accorded with a hybrid formula, consisting of central differencing and second order up-winding. The diffusive terms are, however, discretized by a second order accurate three-point central difference formula. In order to have second order spatial accuracy of the boundary conditions, some fictitious grid points outside the physical domain were considered. To achieve the steady state criterion, one needs to perform at least 8,00,000 iteration steps for the marker and cell (MAC) scheme. Steady state is achieved when the convergence criterion for concentration is 10−6 for each drug form. Interested readers are referred to Saha et al.28 for the detailed numerical procedure.

A typical combined MAC cell for tissue. Abbreviation: MAC, marker and cell.
Fig. 2  A typical combined MAC cell for tissue. Abbreviation: MAC, marker and cell.

Results and discussion

For the purpose of numerical computation of the quantities of physiological significance, the computational domain has been confined to a finite non-dimensional arterial length of 50. For this computational domain, solutions are computed with grid sizes 501 × 101 for δt = 0.0001.

Based on the numerical values in Table 210,14,17,19,22,23,27,29–39 of the model parameters, an extensive quantitative analysis has been performed through graphical representations, the results are presented in Figures 314. These are representative of a first-generation DES (Cypher) which elutes the drug sirolimus. The radial locations-variant normalized concentration profiles for free drug, ECM-bound drug and SR-bound drug concentrations in the tissue for three different times have been shown in Figure 3a–c respectively. It was observed from the figures that with increasing time, the drug masses (free, ECM-bound and SR-bound) decrease (at z = 21.5). The rate of decrease for the free drug is faster than the ECM-bound drug, whereas the SR-bound drug is slower than the ECM-bound drug. As a result, the drug enters the arterial wall (at r = 15) in the free phase and is rapidly bound to both ECM and SR binding sites. The free and ECM-bound drug concentration profiles rise to a peak (shown in Fig. 3d) before decaying with time as the drug traverses through the tissue, becomes bound to SR in the binding phase, and is absorbed at the adventitial boundary (r = 25).

Drug concentration in the vessel wall versus transmural positions for different times at <italic>z</italic> = 21.5.
Fig. 3  Drug concentration in the vessel wall versus transmural positions for different times at z = 21.5.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Non-dimensional drug mass in each phase versus time. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Drug concentration in the vessel wall versus time for different radial positions at <italic>z</italic> = 21.5.
Fig. 4  Drug concentration in the vessel wall versus time for different radial positions at z = 21.5.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different Peclet numbers <italic>Pe<sub>T</sub></italic>.
Fig. 5  Mean drug concentration in the vessel wall versus time for different Peclet numbers PeT.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different values of <italic>Da</italic><sub>1</sub>/<italic>Pe</italic><sub>1</sub>.
Fig. 6  Mean drug concentration in the vessel wall versus time for different values of Da1/Pe1.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different values of <italic>Da</italic><sub>2</sub>/<italic>Pe</italic><sub>2</sub>.
Fig. 7  Mean drug concentration in the vessel wall versus time for different values of Da2/Pe2.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter <italic>α</italic>.
Fig. 8  Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter α.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter <italic>α</italic><sub>1</sub>.
Fig. 9  Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter α1.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter <italic>β</italic>.
Fig. 10  Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter β.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter <italic>β</italic><sub>1</sub>.
Fig. 11  Mean drug concentration in the vessel wall versus time for different values of the non-dimensional parameter β1.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug; (d) Total drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Visual representation of drug concentration in the vessel wall for time-dependent release kinetics with zero-concentration interface condition in <italic>Γ<sub>bt</sub></italic> at <italic>t</italic> = 300.
Fig. 12  Visual representation of drug concentration in the vessel wall for time-dependent release kinetics with zero-concentration interface condition in Γbt at t = 300.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Visual representation of drug concentration in the vessel wall for constant release kinetics with zero-concentration interface condition in <italic>Γ<sub>bt</sub></italic> at <italic>t</italic> = 300.
Fig. 13  Visual representation of drug concentration in the vessel wall for constant release kinetics with zero-concentration interface condition in Γbt at t = 300.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug. Other model parameters are set as in Table 2. Abbreviations: ECM, extracellular matrix; SR, specific receptors.

Visual representation of drug concentration in the vessel wall for time-dependent release kinetics with zero-flux interface condition in <italic>Γ<sub>bt</sub></italic> at <italic>t</italic> = 300.
Fig. 14  Visual representation of drug concentration in the vessel wall for time-dependent release kinetics with zero-flux interface condition in Γbt at t = 300.

(a) Free drug; (b) ECM-bound drug; (c) SR-bound drug. Other model parameters are set as in Table 2. Abbreviatiosn: ECM, extracellular matrix; SR, specific receptors.

Table 2

Plausible values of involved parameters

ParameterValue with unitReference
δ10−4 m19,29,37
CS10−2 mol m−335
Vwall5.8 × 10−8 m s−131,36
cbECMmax3.63 × 10−1 mol m−310
cbSRmax3.3 × 10−3 mol m−314
BM1.3 mol m−327
KonECM2.0 [mol m−3 s]−110
KonSR8.0 × 102 [mol m−3 s]−138
Kd0.136 mol m−327
KdECM2.6 × 10−3 mol m−310
KdSR2.0 × 10−5 mol m−322
krECM5.2 × 10−3 s−114
krSR1.6 × 10−3 s−114
Dfree3.65 × 10−12 m2 s−123,32
λ10−5 s−133,34
εw0.78717
γw130,39
τw1.33317

Although the free and ECM-bound phase profile shapes are similar, drug concentrations within the SR-bound phase are greater than in the ECM-bound phase, which in turn is greater than the free drug concentrations. Within the time t = 10, the SR-bound drug spanning 30% the thickness of the tissue is saturated; these remain saturated for the duration of the time t = 200 studied (Fig. 3c). The remaining SR sites become saturated in the subsequent times and they too remain at saturation levels for the duration of the times t = 300. In Figure 3d, the temporal variation of normalized mean drug concentration profiles reveals that the drug binding in the SR-bound phase is higher than the drug binding in the ECM-bound phase. Therefore, the results of this study demonstrate that drug delivered to the arterial wall from the stent is too low to occupy a large proportion of ECM binding sites, yet is high enough to saturate SR binding sites. This agrees with Tzafriri et al.14

The time-variant concentration profiles for free drug, ECM-bound drug and SR-bound drug concentrations in the tissue for different radial locations are shown in Figure 4a–c respectively. It was observed from these figures that with increasing radial locations (from lumen-tissue interface), the drug masses (free, ECM-bound and SR-bound) decrease (at z = 21.5). The characteristics of the graphs are quite similar, as anticipated, with the binding and unbinding processes taking place simultaneously. Moreover, at the interface, it has been observed that the concentration of drug attains its maximum value for all time, as expected.

Distributions of normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations over the entire period of time for different values of the Peclet number PeT are presented in Figure 5a–d respectively. It is observed that in each case the drug mass is first increasing, up to some upper bound, and then decreasing asymptotically. Evidently, PeT, depending on DT, again DT depends on Deff, increases with a decrease of the porosity (εw) of the arterial wall (as the porosity decreases, the effective as well as true diffusivity does decrease), and also with an increase of the tortuosity (τw) of the arterial wall (as the tortuosity increases, the effective diffusivity as well as true diffusivity does decrease) (keeping Dfree fixed). It is observed from these figures that all the mean drug (free, ECM-bound, SR-bound and total) concentrations decrease with decreasing porosity and increasing tortuosity of the arterial wall (i.e. increase of the Peclet number (PeT)).

The influence of scaling parameter Da1/Pe1 on the normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations in the arterial tissue are displayed in Figure 6a–d respectively over a stipulated period of time. Evidently, Da1/Pe1, depending on konECM and cbECMmax (keeping δ and Vwall fixed), increases/decreases with an increase/decrease in the ECM binding site density (cbECMmax) and also with an increase/decrease of ECM binding on-rate (konECM). It has been observed that with increasing ECM binding site density and ECM binding on-rate (i.e. increase of Da1/Pe1), the mean drug (free, SR-bound and total) concentrations decrease. But, the mean ECM-bound drug concentration increases with increasing ECM binding site density and ECM binding on-rate (i.e. increase of Da1/Pe1).

The influence of scaling parameter Da2/Pe2 on the normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations in the arterial tissue are displayed in Figure 7a–d respectively over a stipulated period of time. Evidently, Da2/Pe2, depending on konSR and cbSRmax (keeping δ and Vwall fixed), increases/decreases with an increase/decrease in the receptor binding site density (cbSRmax) and also with an increase/decrease in receptor binding on-rate (konSR). It has been observed that with increasing receptor binding site density and receptor binding on-rate (i.e. increase of Da2/Pe2), the mean drug (free and ECM-bound) concentrations decrease. But, the mean SR-bound and total drug concentrations increase with increasing receptor binding site density and receptor binding on-rate (i.e. increase of Da2/Pe2). It was also observed in Figure 7a–d that the effects of Da2/Pe2 on SR-bound and total drugs are more sensitive than free and ECM-bound drugs.

The results in Figure 8a–d, respectively, exhibit the influence of scaling parameter α on the normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations in the arterial tissue over the entire period of time. Evidently, α, depending on cbECMmax, increases with a decrease in the ECM binding site density (cbECMmax) (keeping cs fixed). It is an obvious observation that if the ECM binding site density decreases (i.e. the scaling parameter α increases) then all the mean drug (free, ECM-bound, SR-bound and total) concentrations increase.

The results in Figure 9a–d, respectively, exhibit the influence of scaling parameter α1 on the normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations in the arterial tissue over the entire period of time. Evidently, α1, depending on KdECM, increases with an increase of the dissociation (backward) rate constant (KrECM) in the ECM binding site and also with a decrease in ECM binding on-rate (konECM) (keeping cs fixed). It is an obvious observation that if the ECM binding on-rate decreases and dissociation (backward) rate constant in the ECM binding site increases (i.e. the scaling parameter α1 increases) then the mean drug (free, SR-bound and total) concentrations increase. But, the mean ECM-bound drug concentration decreases with decreasing ECM binding on-rate and increasing dissociation (backward) rate constant in the ECM binding site (i.e. increase of scaling parameter α1).

The results, shown in Figure 10a–d, respectively, project the influence of scaling parameter β on the normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations in the arterial tissue over the entire period of time. Evidently, β, depending on cbSRmax, increases with a decrease in the receptor binding site density (cbSRmax) (keeping cs fixed). It is an obvious observation that if the receptor binding site density decreases (i.e. the scaling parameter β increases) then all of the mean drug (free, ECM-bound, SR-bound and total) concentrations increase.

The influence of scaling parameter β1 on the normalized mean free drug, mean ECM-bound drug, mean SR-bound drug and mean total drug concentrations in the arterial tissue are displayed in Figure 11a–d respectively over a stipulated period of time. Evidently, β1, depending on KdSR, increases with an increase in the dissociation (backward) rate constant (KdSR) in the receptor binding site and also with a decrease in receptor binding on-rate (KonSR) (keeping cs fixed). It was observed that with decreasing receptor binding on-rate and increasing dissociation (backward) rate constant in the receptor binding site (i.e. increase of β1), the mean drug (free and ECM-bound) concentrations increase. But, the mean SR-bound and total drug concentrations decrease with decreasing receptor binding on-rate and increasing dissociation (backward) rate constant in the receptor binding site (i.e. increase of scaling parameter β1). It was also observed from Figure 11a–d that the effects of scaling parameter β1 on SR-bound and total drugs are more sensitive than free and ECM-bound drugs.

Finally, in Figure 12a–c, respectively, the displayed spatial distribution of free, ECM-bound and SR-bound drug concentration may again clearly justify the reduction of late lumen loss at the distal part of the arterial tissue, which again validates the findings of Balakrishnan et al.29 The spatial patterns for free, ECM-bound and SR-bound drug concentration in Figure 13a–c, respectively, clearly establish our findings further for the zero-concentration interface condition and constant release kinetics. In the Figure 14a–c, respectively, the spatial distribution of free, ECM-bound and SR-bound drug concentration for the zero-flux interface condition and time-dependent release kinetics is displayed.

Conclusions

In the present work, a novel analytical closed-form solution of a 2D axi-symmetric model of drug transport eluted from a coronary DES is proposed and focused on the reversible and saturable binding processes in the vascular tissue. The model is based on a single-layered homogeneous multiple-phase system where a system of partial differential equations describes both the dissolution and diffusion processes in the polymeric layer as well as diffusion, convection and reaction in the tissue. The closed-form solution has been established by using the MAC method.

The salient observations of the above findings are the following:

Future research directions

Application of this framework to idealized configurations of arteries stented with a DES yielded some general guidelines for future DES design, especially concerning strategies for the effective elution of the anti-proliferative drug from the stent and its efficacy. Despite these very useful results, many improvements can be envisaged in future models. Limitations of the current model are that the model geometry is 2D axi-symmetric, the considered stent is idealized having three stent struts, and the model is limited to a straight vessel segment. Although the target zone is an atherosclerotic plaque, this consideration has been disregarded here. Inclusion of realistic plaque with anisotropic tissue properties, together with the varying diffusivity of eluted drug within the target lesion, may be the scope of future research.

It will be important to eventually perform the simulations on more realistic three-dimensional geometries where the detailed structure of the stent design is to be taken into account. Based upon the present model validation, the model can be evolved further with the incorporation of several factors, depending on the objectives of the drug release phenomena in various situations, for designing future research in this direction by utilizing the present knowledge of the system.

Medical benefits

Mathematical modelling and numerical simulation are indispensable tools when clinical investigation and/or animal studies are expensive, and in some cases, cumbersome as well. A mathematical model can give us an idea of how an underlying mechanism plays a surrogative role on cardiovascular intervention. The present study gives an overview on the release kinetics and also the binding of sirolimus drug eluted from a coronary stent which helps, certainly, the clinician to estimate the effectiveness of delivery and also the efficacy of drug. Although the present study is an idealized one, it will give an idea to the manufacturer in designing next generation of DES. We point out that the results presented here are for the simulated case of sirolimus release and absorption from the Cypher stent.

The significant influence of estimated parameters on the drug masses has been shown graphically, which establishes the strong fact that by altering the parameters, various means of drug release control can be achieved according to the patients’ needs. Various conclusions can be drawn from the dynamical behavior of the present model study. It is worth noting that as the drug needs a much longer time in the tissue to get absorbed completely, its effect will certainly persist for a long time before the same drug is administered subsequently. This approach may also be applied in any form to other parts of the human body, provided the system does not have major clinical complexity.

Abbreviations

2D: 

two-dimensional

BMS: 

bare metal stents

DES: 

drug-eluting stents

ECM: 

extracellular matrix

MAC: 

marker and cell

PCI: 

percutaneous coronary interventions

SDR: 

specific drug binding

SER: 

stent elution rate

SMC: 

smooth muscle cells

SR: 

specific receptors

Declarations

Acknowledgement

The author gratefully acknowledges the careful scrutiny and suggestions of the learned reviewers. The author also expresses his gratitude to Professor Prashanta Kumar Mandal, Department of Mathematics, Visva-Bharati University, Santiniketan-731235, West Bengal, India, for fruitful discussions about this paper and for useful support in some numerical tasks.

Conflict of interest

The author declares that no conflict of interest exists.

Authors’ contributions

RS is the sole author.

References

  1. Aragon J, Kar S, Tio F, Trauthen B, Parisky A, Watanabe C. The effect of variable release kinetics on paclitaxel efficacy from a drug eluting stent in a porcine model. EuroIntervention 2005;1:228-235 View Article PubMed/NCBI
  2. Guagliumi G, Ikejima H, Sirbu V, Bezerra H, Musumeci G, Lortkipanidze N. Impact of drug release kinetics on vascular response to different zotarolimus-eluting stents implanted in patients with long coronary stenoses: the longoct study (optical coherence tomography in long lesions). JACC Cardiovasc Interv 2011;4:778-785 View Article PubMed/NCBI
  3. Lovich MA, Philbrook M, Sawyer S, Weselcouch E, Edelman ER. Arterial heparin deposition: role of diffusion, convection, and extravascular space. Am J Physiol Heart Circ Physiol 1998;275:H2236-H2242 View Article PubMed/NCBI
  4. Moses JW, Stone GW, Nikolsky E, Mintz GS, Dangas G, Grube E. Drug-eluting stents in the treatment of intermediate lesions: pooled analysis from four randomized trials. J Am Coll Cardiol 2006;47:2164-2171 View Article PubMed/NCBI
  5. Serruys PW, Sianos G, Abizaid A, Aoki J, den Heijer P, Bonnier H. The effect of variable dose and release kinetics on neointimal hyperplasia using a novel paclitaxel-eluting stent platform: the paclitaxel in-stent controlled elution study (pisces). J Am Coll Cardiol 2005;46:253-260 View Article PubMed/NCBI
  6. Migliavacca F, Gervaso F, Prosi M, Zunino P, Minisini S, Formaggia L. Expansion and drug elution model of a coronary stent. Comput Methods Biomech Biomed Engin 2007;10:63-73 View Article PubMed/NCBI
  7. Borghi A, Foa E, Balossino R, Migliavacca F, Dubini G. Modelling drug elution from stents: effects of reversible binding in the vascular wall and degradable polymeric matrix. Comput Methods Biomech Biomed Engin 2008;11:367-377 View Article PubMed/NCBI
  8. Horner M, Joshi S, Dhruva V, Sett S, Stewart S. A two-species drug delivery model is required to predict deposition from drug-eluting stents. Cardiovasc Eng Technol 2010;1:225-234 View Article PubMed/NCBI
  9. Ferreira JA, Naghipoor J, Oliveira Pd. A coupled non-fickian model of a cardiovascular drug delivery system. Math Med Biol 2016;33:329-357 View Article PubMed/NCBI
  10. Tzafriri AR, Levin AD, Edelman ER. Diffusion-limited binding explains binary dose response for local arterial and tumour drug delivery. Cell Proliferation 2009;42:348-363 View Article PubMed/NCBI
  11. Bozsak F, Chomaz JM, Barakat AI. Modeling the transport of drugs eluted from stents: physical phenomena driving drug distribution in the arterial wall. Biomech Model Mechanobiol 2014;13:327-347 View Article PubMed/NCBI
  12. Bozsak F, Gonzalez-Rodriguez D, Sternberger Z, Belitz P, Bewley T, Chomaz JM. Optimization of drug delivery by drug-eluting stents. PLoS one 2015;10:e0130,182 View Article PubMed/NCBI
  13. Lauffenburger DA, Linderman J. Receptors: models for binding, trafficking, and signaling. Oxford University Press; 1993 View Article PubMed/NCBI
  14. Tzafriri AR, Groothuis A, Price GS, Edelman ER. Stent elution rate determines drug deposition and receptor-mediated effects. J Control Release 2012;161:918-926 View Article PubMed/NCBI
  15. McGinty S, Pontrelli G. On the role of specific drug binding in modelling arterial eluting stents. J Math Chem 2016;54:967-976 View Article PubMed/NCBI
  16. McGinty S, Pontrelli G. A general model of coupled drug release and tissue absorption for drug delivery devices. J Control Release 2015;217:327-336 View Article PubMed/NCBI
  17. O’Connell BM, Walsh MT. Demonstrating the influence of compression on artery wall mass transport. Ann Biomed Eng 2010;38:1354-1366 View Article PubMed/NCBI
  18. Balakrishnan B, Dooley JF, Kopia G, Edelman ER. Intravascular drug release kinetics dictate arterial drug deposition, retention, and distribution. J Control Release 2007;123:100-108 View Article PubMed/NCBI
  19. Balakrishnan B, Dooley J, Kopia G, Edelman ER. Thrombus causes fluctuations in arterial drug delivery from intravascular stents. J Control Release 2008;131:173-180 View Article PubMed/NCBI
  20. Hwang CW, Wu D, Edelman ER. Physiological transport forces govern drug distribution for stent-based delivery. Circulation 2001;104:600-605 View Article PubMed/NCBI
  21. Levin AD, Jonas M, Hwang CW, Edelman ER. Local and systemic drug competition in drug-eluting stent tissue deposition properties. J Control Release 2005;109:236-243 View Article PubMed/NCBI
  22. Bierer BE, Mattila PS, Standaert RF, Herzenberg LA, Burakoff SJ, Crabtree G. Two distinct signal transmission pathways in t lymphocytes are inhibited by complexes formed between an immunophilin and either fk506 or rapamycin. Proceedings of the National Academy of Sciences 1990;87:9231-9235 View Article PubMed/NCBI
  23. Kolachalama VB, Levine EG, Edelman ER. Luminal flow amplifies stent-based drug deposition in arterial bifurcations. PloS one 2009;4:e8105 View Article PubMed/NCBI
  24. Kolandaivelu K, O’Brien CC, Shazly T, Edelman ER, Kolachalama VB. Enhancing physiologic simulations using supervised learning on coarse mesh solutions. 2015;12:20141,073 View Article PubMed/NCBI
  25. Higuchi T. Mechanism of sustained-action medication. theoretical analysis of rate of release of solid drugs dispersed in solid matrices. J Pharm Sci 1963;52:1145-1149 View Article PubMed/NCBI
  26. O’Brien CC, Kolachalama VB, Barber TJ, Simmons A, Edelman ER. Impact of flow pulsatility on arterial drug distribution in stent-based therapy. J Control Release 2013;168:115-124 View Article PubMed/NCBI
  27. Kolachalama VB, Pacetti SD, Franses JW, Stankus JJ, Zhao HQ, Shazly T. Mechanisms of tissue uptake and retention in zotarolimus-coated balloon therapy. Circulation 2013;127:2047-2055 View Article PubMed/NCBI
  28. Saha R, Mandal PK. Modelling time-dependent release kinetics in stent-based delivery. J Explor Res Pharmacol 2018;3:61-70 View Article PubMed/NCBI
  29. Balakrishnan B, Tzafriri AR, Seifert P, Groothuis A, Rogers C, Edelman ER. Strut position, blood flow, and drug deposition implications for single and overlapping drug-eluting stents. Circulation 2005;111:2958-2965 View Article PubMed/NCBI
  30. Ferdous J, Chong C. Effect of atherosclerotic plaque on drug delivery from drug-eluting stent. , ; View Article PubMed/NCBI
  31. Huang Z, Tarbell J. Numerical simulation of mass transfer in porous media of blood vessel walls. Am J Physiol 1997;273:H464-H477 View Article PubMed/NCBI
  32. Levin AD, Vukmirovic N, Hwang CW, Edelman ER. Specific binding to intracellular proteins determines arterial transport properties for rapamycin and paclitaxel. Proc Natl Acad Sci U S A 2004;101:9463-9467 View Article PubMed/NCBI
  33. McGinty S. View Article PubMed/NCBI
  34. McGinty S, McKee S, Wadsworth RM, McCormick C. Modelling drug-eluting stents. Math Med Biol 2011;28:1-29 View Article PubMed/NCBI
  35. Sakharov DV, Kalachev LV, Rijken DC. Numerical simulation of local pharmacokinetics of a drug after intravascular delivery with an eluting stent. J Drug Target 2002;10:507-513 View Article PubMed/NCBI
  36. Tedgui A, Lever M. Filtration through damaged and undamaged rabbit thoracic aorta. Am J Physiol Heart Circ Physiol 1984;247:H784-H791 View Article PubMed/NCBI
  37. Vairo G, Cioffi M, Cottone R, Dubini G, Migliavacca F. Drug release from coronary eluting stents: a multidomain approach. J Biomech 2010;43:1580-1589 View Article PubMed/NCBI
  38. Wear MA, Walkinshaw MD. Determination of the rate constants for the fk506 binding protein/rapamycin interaction using surface plasmon resonance: An alternative sensor surface for ni 2+–nitrilotriacetic acid immobilization of his-tagged proteins. Anal Biochem 2007;371:250-252 View Article PubMed/NCBI
  39. Zunino P. Multidimensional pharmacokinetic models applied to the design of drug-eluting stents. Cardiovascular Engineering: An International Journal 2004;4:181-191 View Article PubMed/NCBI