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

Modelling Time-dependent Release Kinetics in Stent-based Delivery

  • Ramprosad Saha1 and
  • Prashanta Kumar Mandal2,*
Journal of Exploratory Research in Pharmacology   2018;3(2):61-70

doi: 10.14218/JERP.2018.00001

Received:

Revised:

Accepted:

Published online:

 Author information

Citation: Saha R, Mandal PK. Modelling Time-dependent Release Kinetics in Stent-based Delivery. J Explor Res Pharmacol. 2018;3(2):61-70. doi: 10.14218/JERP.2018.00001.

Abstract

Background and objective

The present study deals with a computational model of the transport and retention of drug within the arterial wall eluted from a drug-eluting stent, to enhance our understanding of the performance of this device. We considered a two-species model (free and bound) incorporating a reversible reaction to describe drug interactions with the constituents of the arterial wall. An axisymmetric model of drug delivery from a pair of stent struts has been developed, where the transport of free drug is modelled as an unsteady reaction-diffusion process, while the bound drug, assuming complete immobilization in the tissue, is modeled as an unsteady reaction process. The model also took into account a second-order binding process that describes a saturating reversible binding and time-dependent release kinetics of the drug-eluting stent. Considering that diffusion takes place over a tortuous path in a porous media, the effects of porosity and tortuosity on diffusion cannot be ruled out from this present investigation.

Methods

An explicit finite-difference scheme is leveraged to tackle numerically the governing equations of motion, together with the physiologically realistic boundary conditions.

Results

The quantitative effects of significant factors, such as Peclet number (PeT), Damköhler number (Da), tortuosity, porosity, interstrut distance and time-dependent release kinetics of drug-eluting stent, on the distribution and retention of drug are determined graphically.

Conclusions

Predicted results are consistent with several existing results in the literature, which certainly validates the applicability of the model considered.

Keywords

Two-species model, Drug-eluting stent, Reaction-diffusion process, Time-dependent release kinetics, Interstrut distance, Porosity, Tortuosity

Introduction

Arterial stenting, acting as a supporting scaffold, has revolutionized the treatment of coronary artery disease, as it reopens the occluded vessel and restores the normal flow of blood. 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 arterial stents coated by drug–the so-called drug-eluting stents (DESs) raised all sorts of questions by releasing antiproliferative agents in a controlled manner into the injured site to reduce restenosis rates.1–9 Drug-eluting stents are now the primary choice of percutaneous coronary interventions (PCIs) in millions of patients, but questions regarding their longevity and safety still arise.

Some experimental studies have been carried out in the recent past, with the aim of quantifying the capability of this device to reduce in-stent restenosis after stent implantation.10–12 Lovich et al.10 studied the behaviour of heparin in explanted arteries and concluded that the presence of binding sites changes along the transmural direction, being higher in the endothelium and lower in the adventitia. Lovich and Edelman studied the effects of specific binding sites inside the arterial wall on the drug uptake,13 where the presence of specific binding site action was modelled using the reversible chemical reaction. Sakharov et al.14 disregarded the convective effects on the transport of free drug. Hwang et al.15 predicted the free as well as bound drug concentrations by solving for distribution of free drug and then using a multiplicative factor (partition approach) to predict the concentration of bound drug. Migliavacca et al.16 studied the drug release pattern in vascular wall from drug-eluting stents using a single species approach along with a partition coefficient approach to relate the free and the bound drug concentrations. Borghi et al.17 opined that the inclusion of reversible binding leads to delayed release and that the erosion of polymer affects the drug release from a single strut. Horner et al.18 considered a two-species drug delivery model including reversible binding sites, and their model predicted that a single species drug delivery model cannot accurately predict the distribution of bound drug. They also concluded that a two-species approach that includes reversible binding is the way forward for future stent-based drug delivery systems.

The main objective of this study was to advance the aforementioned work with a two-species model of drug transport eluted from a pair of struts, where the transport of free drug is governed by reaction-diffusion process and that of bound drug, assuming complete immobilization in the tissue, by reaction process. Following Tzafriri et al.,19 a second-order dynamic model that describes a saturating reversible binding process by treating bound drug as a dynamic variable has been taken into account to explore drug interaction with cells of the arterial wall. In most of the studies cited above, transient drug release has been modelled as a uniform release, which is unrealistic and not representative of actual stent-based delivery. Instead, a simple time-dependent Dirichlet boundary condition is applied on the surface of the struts.20–24 Arterial properties, such as porosity and tortuosity, dictate the transport of drugs within the arterial tissue.

When an endovascular drug-eluting stent is implanted, it has major impact on the structure of the arterial wall, eventually influencing the overall rates of diffusion through tissues.25 For diffusion in a porous material, the effective diffusion coefficient is assumed to depend on two factors viz porosity (a dimensionless parameter, which is the ratio of pore volume to the total material volume) and diffusion path tortuosity (ratio of the actual pore length to the distance between its ends; i.e. arc-chord ratio)26—these parameters change the free diffusivity of the drug eluted from a pair of struts.27 Thus, the effects of porosity and tortuosity on the diffusivity of drug have also been accounted for in this present investigation.

Governing equations and boundary conditions

The freely-transported drug, referred to as ‘free drug’, is denoted by cf and the drug, completely immobilized in the tissue, is referred to as ‘bound drug’, which is denoted by cb. The interconversion of drug between the unbound plasma phase and the bound phase of tissue binding sites is controlled by a second-order chemical reaction. 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 (Fig. 1). The transport of free drug eluted from a pair of struts is governed by unsteady reaction-diffusion process (Eq. 1) and that of the bound drug is represented by unsteady reaction process (Eq. 2).

Schematic diagram of the computational domain.
Fig. 1  Schematic diagram of the computational domain.

Symmetry boundary conditions for both the free and the bound drugs are applied at the proximal (Γti) and the distal (Γto) walls (Eq. 3).28,29 Impermeable boundary condition for bound drug is assumed at the perivascular wall (Γtp), lumen-tissue (Γbt) and strut-tissue (Γst) interfaces (Eq. 4). For the free drug, perfect sink condition is imposed at the perivascular end (Eq. 5). Since a proper boundary condition for the free drug at lumen-tissue interface (Γbt) is not readily apparent, we considered two opposing extremes, either that flowing blood is extremely efficient at washing out mural-adhered drug, modelled as a zero-concentration interface condition, or that mural-adhered drug is insensitive to flowing blood, modelled as zero-flux boundary condition (Eq. 6).30–31 Instead of modelling uniform release of drug from a pair of struts, we assumed simple time-dependent release kinetics with a flux condition (Eq. 7).13,20 Therefore, the governing equations representing the transport of free as well as bound drugs together with their appropriate boundary conditions are as follows:

cft=DT[2cfr2+1rcfr+2cfz2]     [kacf(BMcb)kdcb],
cbt=[kacf(BMcb)kdcb],
cfz=0=cbz on Γti and Γto,
cbr=0 on Γtp and Γtl(=ΓbtΓst),
cf=0 on Γtp,
cf=0 or cfr=0 on Γbt,
Jb(t)=Dcc02πt on Γst,
where z and r coordinates are along the axial and the radial directions respectively. ka and kd are the rates of association and dissociation constants respectively. BM is the net tissue binding capacity. Dc and c0 are the diffusion coefficient and initial concentration in coating respectively.

Here, DT, the true diffusivity of the free drug can be written as:27,32

DT=(1+BMRd)×Deff,
where
Deff=ετ×Dfree.
ε and τ are the porosity and the tortuosity of the wall material respectively. Dfree and Deff are the coefficients of free and effective diffusivity respectively. Rd(= kd/ka) is the equilibrium dissociation constant.

We now introduce nondimensional variables as

x=rδ,z¯=zδ,t¯=tVyδ,c¯f=cfc0,c¯b=cbBM,
where Vy is the transmural filtration velocity.

Under these assumptions, the equations (17) take their respective nondimensional forms (dropping tilde) as

cft=1PeT[2cfx2+1xcfx+2cfz2]DaPeT[cf(1cb)ε1cb],
cbt=ε2DaPeT[cf(1cb)ε1cb],
cfz=0=cbz on Γti and Γto,
cbx=0 on Γtp and Γtl(=ΓbtΓst),
cf=0 on Γtp,
cf=0 or cfx=0 on Γbt,
Jb(t)=1Pecπt on Γst,
where PeT[= (Vyδ)/(DT)] and Da[= (kaBMδ2)/(DT)] are the Peclet number and Damköhler number in the tissue respectively. Here, ε1(= Rd/c0) and ε2(= c0/BM) are two scaling parameters. Pec{= [Vy(h2/δ)]/Dc} is the Peclet number in the coating of struts and h is the thickness of the coating of the strut.

Solution procedure

The governing equations (89) representing the transport of free and bound drugs along with the boundary conditions (10–14) are solved numerically by finite-difference scheme. Forward-time centered-space discretization technique is leveraged in this explicit scheme. Following, we describe our finite-difference scheme in more detail. We denote xj = jδx, zi = iδz, tn = nδt, where n refers to the time level and δt is the time increment. Here, δx and δz stand for the step sizes along the radial and the axial directions respectively. Steady state is achieved when the convergence criterion for concentration was 10−7 for both drug forms.

The finite-difference approximation of (8) for the transport of free drug is given by

cfi,jn+1cfi,jnδt=[diffcf]i,jn[reac]i,jn,
where [diffcf]i,jn=(1/PeT) {[(2cf)/(x2)]i,jn+ 1/xj [(cf)/(x)]i,jn + [(2cf)/(z2)]i,jn} and [reac]i,jn= (Da/PeT) [cfi,j(1cbi,j)ε1cbi,j]

Likewise, the discretized version of Eq. 9 may be written as

cbi,jn+1cbi,jnδt=ε2[reac]i,jn.

Result and discussion

For the purpose of numerical computation of the desired quantities of major physiological significance by using the baseline values included in Table 17,13,23,28,33–41, solutions are computed through the generation of grids, with a size of 301 × 101 for δt = 0.00001. The simulation concerning the grid independence study was performed for the purpose of examining the error associated with the grid sizes used, and is presented in Figure 2a. One may notice from this figure that the transmural variation of normalized free drug concentration concerning the three grid sizes viz 151 × 51, 301 × 101, 601 × 201 overlap one another at t = 100 for PeT = 2 and Da = 40. Thus, the grid independence study in the present context of numerical simulation has its own importance to establish the correctness of the results obtained. Figure 2b displays the results of normalized free drug concentration for different time steps δt viz 0.0001, 0.00001, 0.00005. Here too, the concentration profiles get overlapped for three distinct time steps, which also establish the correctness of the time step used.

Table 1

Plausible values of involved parameters

DescriptionParameterValueReference
Strut dimension, mδ0.0001[13]
Strut coating thickness, mh5.0 × 10−5[33]
Mean wall thickness, mAw (= 10δ)0.001[34]
Interstrut distance, mΔ (= 7δ)0.0007[35]
Transmural filtration velocity, m s−1Vy4 × 10−8[36]
Porosity of the arterial wallε0.787[37]
Tortuosity of the arterial wallτ1.333[37]
Coating drug diffusivity, m2 s−1Dc1 × 10−15[38]
Coefficient of free diffusivity, m2 s−1Dfree3.65 × 10−12[28,39]
Coefficient of effective diffusivity, m2 s−1Deff2.15 × 10−12Our study
True diffusivity of the free drug, m2 s−1DT24 × 10−12Our study
Initial drug concentration in the coating, mol m−3c01[23]
Net tissue binding capacity, mol m−3BM0.01[7]
Association rate constant, m3 mol−1 s−1ka10[40]
Dissociation rate constant, s−1kd0.01[41]
Equilibrium dissociation constant, mol m−3Rd0.001Our study
Dimensionless Peclet number in the coatingPec100Our study
Dimensionless Peclet number in the tissuePeT2Our study
Dimensionless Damköhlar number in the tissueDa40Our study
Dimensionless scaling parameterε10.001Our study
Dimensionless scaling parameterε2100Our study
Transmural variation of normalized free drug concentration for <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40 and <italic>t</italic> = 100 at <italic>z</italic> = 11.
Fig. 2  Transmural variation of normalized free drug concentration for PeT = 2, Da = 40 and t = 100 at z = 11.

a: Different grid sizes; b: different time steps.

Distributions of mean free as well as bound drug concentrations for different values of scaling parameter ε1(= Rd/c0) are presented in Figure 3a and 3b respectively, and the same for different values of ε2(= c0/BM), which are depicted in Figure 4a and 4b respectively. Evidently, ε1, depending on Rd (= kd/ka), decreases with a decrease in the dissociation rate constant kd and also with an increase in the association rate constant ka. However, ε2 increases with decreasing BM (keeping c0 as fixed). Figure 3a shows that the normalized mean free drug concentration decreases with decreasing ε1 for PeT = 2, Da = 40, ε2 = 100, up to a certain time and, thereafter, no significant changes occurred. It may be justified in the sense that, as ε1 decreases, the rate of reversible binding (kd) decreases and/or the rate of forward binding increases, which is lowering the mean concentration of free drug.

Distribution of normalized mean drug concentration for different values of ε<sub>1</sub> at <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40, ε<sub>2</sub> = 100 and <italic>t</italic> = 500.
Fig. 3  Distribution of normalized mean drug concentration for different values of ε1 at PeT = 2, Da = 40, ε2 = 100 and t = 500.

a: Free drug; b: bound drug.

Distribution of normalized mean drug concentration for different values of ε<sub>2</sub> at <italic>Pe<sub>T</sub></italic> = 2, <italic>Da =</italic> 40, ε<sub>1</sub> = 0.001 and <italic>t</italic> = 500.
Fig. 4  Distribution of normalized mean drug concentration for different values of ε2 at PeT = 2, Da = 40, ε1 = 0.001 and t = 500.

a: Free drug; b: bound drug.

An interesting phenomenon may be observed when the transported free drug is insensitive to flowing blood [i.e. modelled as zero-flux lumen-tissue interface condition, the mean concentration of free drug is always higher (zoom inset) as compared to zero-concentration interface condition and approaches a quasi-steady state faster]. Figure 3b shows how the rates of forward as well as reversible binding affect the mean concentration of bound drug within the arterial tissue. It is to be observed that the mean concentration is increased with the decrease of ε1, owing to the increase in the rate of forward binding and/or to the decrease in the rate of reversible binding. Here too, the mean concentration of bound drug appears to be higher in case of zero-flux lumen-tissue interface condition, and saturation of binding sites takes place very rapidly.

Effects of ε2 (i.e. net tissue binding potential on the mean concentrations of free as well as bound drug) are displayed in Figure 4a and 4b respectively. It may be recalled that ε2 increases with decreasing binding potential. The results of these figures indicate that the mean concentration of free drug increases with decreasing binding potential up. to ε2 = 100, but the concentration reaches a quasi-steady state for weaker binding capacity (ε2 = 1,000) than the other cases. As anticipated, the mean concentration of free drug appears to be higher for zero-flux lumen-tissue interface condition (zoom inset) (Fig. 4a). Similar observations can be made for normalized bound drug concentration but with a distinction. In the case of that free drug for ε2 ≤ 100, the quasi-equilibrium is not fully established until t = 500. On the other hand, in the case that bound drug for ε2 ≥ 10, the quasi-equilibrium is attained very rapidly (Fig. 4b). Furthermore, the effects of zero-flux lumen-tissue interface condition on the spatiotemporal distribution of free and bound drug can be visualized in Figure 5, panels a and b.

Visual representation of drug concentration at <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40 and <italic>t =</italic> 500.
Fig. 5  Visual representation of drug concentration at PeT = 2, Da = 40 and t = 500.

a: Free drug (Jb = 0 on Γbt for cf); b: bound drug (Jb = 0 on Γbt for cf).

The time-evolution of the normalized free and bound drug concentrations at positions P1 and P2 (near the left and the perivascular boundaries in Fig. 1) are shown in Figure 6, panels a, b and c, d respectively. The drug distribution profiles indicate that quasi-equilibrium has not yet been fully established until t = 500 for small values of Da and also that the normalized free drug concentration at P1 and P2 decreases with increasing Da (Fig. 6a and 6b). Moreover, the reduced local free drug concentration at P2 with increasing Da further reduces the drug dissipation at the perivascular boundary. Similar results on the distribution of bound drug have been shown in Figure 6c and 6d. All of the above observations are in good agreement with those of Hwang et al.15 who studied the distribution of free drug only. The lack of drug far away from the strut (i.e. at P2) can be a serious factor for in-stent restenosis occurrence, as high drug concentration at the perivascular end is more important in suppressing restenosis.42–44

Temporal variation of normalized drug concentration at different positions of the domain for different <italic>Da</italic> at <italic>Pe<sub>T</sub></italic> = 2 and <italic>t</italic> = 500.
Fig. 6  Temporal variation of normalized drug concentration at different positions of the domain for different Da at PeT = 2 and t = 500.

a: Free drug at P1; b: free drug at P2, c bound drug at P1; d: bound drug at P2.

The influences of diffusivity (Peclet number) and time-dependent release kinetics of a drug on the temporal variation of free as well as bound drug concentrations within the arterial tissue are portrayed in Figure 7, panels a, b and c, d respectively. Figure 7a and 7b indicate that the free drug concentrations at both P1 and P2 do not attain quasi-equilibrium state until t = 500, due to the time-dependent release-kinetics. However, if one assumes constant release of drug from a well apposed strut, a quasi-steady state has been fully established at the positions considered. On the contrary, the bound drug concentration attains a quasi-steady state irrespective of time-dependent and constant release of drug. The effect of release kinetics on the spatial distribution of free as well as bound drug can also be visualized clearly in Figure 8, panels a–d, in which heterogeneous distribution and retention of drug are found to be observed throughout the domain and higher concentration of free drug is observed in case of time-dependent release of drug from struts.

Temporal variation of normalized drug concentration at different positions of the domain for different <italic>Pe<sub>T</sub></italic> at <italic>Da</italic> = 40 and <italic>t</italic> = 500.
Fig. 7  Temporal variation of normalized drug concentration at different positions of the domain for different PeT at Da = 40 and t = 500.

a: Free drug at P1; b: free drug at P2, c bound drug at P1; d: bound drug at P2.

Visual representation of drug concentration at <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40 and <italic>t</italic> = 500.
Fig. 8  Visual representation of drug concentration at PeT = 2, Da = 40 and t = 500.

a: Free (time-dependent release); b: bound (time-dependent release); c: free (constant release); d: bound (constant release).

Figure 9a shows the concentration profile at a height of one strut thickness for free drug concentration depending on interstrut distance for PeT = 2 and Da = 40 at t = 500. A single peak is noted when the struts are placed one-strut width apart. As the interstrut distance increases, the peak concentration falls and distinct peaks over each strut are observed. The distributions of bound drug concentration are analogous to those of the distributions of free drug qualitatively (Fig. 9b).

Distribution of normalized drug concentration for various interstrut distances at one-strut height within tissue for <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40 and <italic>t</italic> = 500.
Fig. 9  Distribution of normalized drug concentration for various interstrut distances at one-strut height within tissue for PeT = 2, Da = 40 and t = 500.

a: Free drug; b: bound drug.

When a stent is implanted within an artery, it has a major impact on the structures of the wall and subsequently alters the transport and retention of the drug within the arterial tissue. Any porous material under compression will demonstrate smaller pore ratio than that of its relaxed state, which in turn influences the tortuosity. Arterial property such as porosity (ε) dictates the transport and retention of drug eluted from a pair of struts, as depicted in Figure 10a and 10b respectively. It is observed from these figures that the mean concentrations of free as well as bound drug increases with increasing ε. However, the enhancement is maximum when ε is allowed to change from 0.2 to 0.78, up to t = 250 for both drug forms. It is interesting to note that the mean concentration of bound drug attains a quasi-steady state for larger ε (i.e. saturation of binding sites takes places very rapidly).

Temporal variation of normalized mean drug concentration for different porosities (ε) at <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40 and <italic>t</italic> = 500.
Fig. 10  Temporal variation of normalized mean drug concentration for different porosities (ε) at PeT = 2, Da = 40 and t = 500.

a: Free drug; b: bound drug.

The above observation may be justified in the sense that as the porosity increases, the effective as well as true diffusivity do increase, 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. Our simulation also demonstrates the fact that a decrease in the mean concentration of free drug with increasing tortuosity (τ) (Fig. 11a) (i.e. an inverse relationship between free drug concentration and tortuosity is revealed). A similar pattern is also observed for bound drug (Fig. 11b). The above 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 eventually lowering the mean concentration of both drug forms is impeded).

Normalized mean drug concentration for varying tortuosity (τ) at <italic>Pe<sub>T</sub></italic> = 2, <italic>Da</italic> = 40 and <italic>t</italic> = 500.
Fig. 11  Normalized mean drug concentration for varying tortuosity (τ) at PeT = 2, Da = 40 and t = 500.

a: Free drug; b: bound drug.

Conclusions and future work

In this numerical study, a two-dimensional axisymmetric model of drug transport eluted from drug-eluting stents has been developed. We consider a two-species model, capable of predicting the delivery of drug and its retention in the arterial tissue where the transport of free drug is governed by reaction-diffusion process and that of bound drug, assuming complete immobilization in the tissue, by reaction process. Following Tzafriri et al.19 a second-order dynamic model that describes a saturating reversible binding process has been taken into account. A time-dependent release kinetic is applied on the surface of the struts.20,23 As the arterial structures such as porosity and tortuosity affect the diffusion process, their influence on the transport and retention of drug eluted from a pair of struts has been investigated successfully.

Though experimental studies can provide information on release kinetics and other histological information, computational studies can provide detailed predictions of the drug distribution and its retention within the arterial tissue. With the rapid ascent of stent-based drug delivery in the treatment of vascular disease, many important issues concerning drug delivery and its retention in the arterial tissue need to be addressed. Indeed, arterial vessels with disease are the target sites for endovascular intervention, and quantifying the pharmacokinetics for this scenario is needed towards full appreciation of drug-eluting stents and like technologies.

Abbreviations

BMS: 

bare metal stent

CAD: 

coronary artery disease

DES: 

drug-eluting stent

ISR: 

in-stent restenosis

PCI: 

percutaneous coronary intervention

Declarations

Acknowledgement

The authors gratefully acknowledge the careful scrutiny and suggestions of the learned reviewers. The authors gratefully acknowledge the partial financial support from Special Assistance Programme (SAP-III) sponsored by University Grants Commission (UGC), New Delhi, India (Grant No. F.510/3/DRS-III/2015(SAP-I)).

Conflict of interest

The authors have no conflict of interests related to this publication.

Authors’ contributions

Formulating the problem (RS, PKM); running the in-house code (RS); writing the manuscript (RS, PKM).

References

  1. 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
  2. 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
  3. Costa MA, Simon DI. Molecular basis of restenosis and drug-eluting stents. Circulation 2005;111:2257-2273 View Article PubMed/NCBI
  4. Daemen J, Serruys PW. Drug-eluting stent update 2007 part I. A survey of current and future generation drug-eluting stents: meaningful advances or more of the same?. Circulation 2007;116:316-328 View Article PubMed/NCBI
  5. Hose D, Narracott A, Griffiths B, Mahmood S, Gunn J, Sweeney D. A thermal analogy for modelling drug elution from cardiovascular stents. Comput Methods Biomech Biomed Engin 2004;7:257-264 View Article PubMed/NCBI
  6. Hwang CW, Levin AD, Jonas M, Li PH, Edelman ER. Thrombosis modulates arterial drug distribution for drug-eluting stents. Circulation 2005;111:1619-1626 View Article PubMed/NCBI
  7. 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
  8. McGinty S, McKee S, Wadsworth RM, McCormick C. Modelling drug-eluting stents. Math Med Biol 2011;28:1-29 View Article PubMed/NCBI
  9. McGinty S, McKee S, Wadsworth RM, McCormick C. Modeling arterial wall drug concentrations following the insertion of a drug-eluting stent. SIAM J Appl Math 2014;73:2004-2028 View Article PubMed/NCBI
  10. Lovich MA, Philbrook M, Sawyer S, Weselcouch E, Edelman ER. Arterial heparin deposition: role of diffusion, convection, and extravascular space. Am J Physiol 1998;275:H2236-H2242 View Article PubMed/NCBI
  11. 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
  12. 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
  13. Lovich MA, Edelman ER. Computational simulations of local vascular heparin deposition and distribution. Am J Physiol 1996;271:H2014-H2024 View Article PubMed/NCBI
  14. 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
  15. Hwang C, Wu D, Edelman E. Physiological transport forces govern drug distribution for stent-based delivery. Circulation 2001;104:600-605 View Article PubMed/NCBI
  16. 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
  17. 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
  18. 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. Cardiovascular Engineering and Technology 2010;1:225-234 View Article PubMed/NCBI
  19. Tzafriri AR, Levin AD, Edelman ER. Diffusion-limited binding explains binary dose response for local arterial and tumour drug delivery. Cell Prolif 2009;42:348-363 View Article PubMed/NCBI
  20. 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
  21. Mandal AP, Mandal PK. An unsteady analysis of arterial drug transport from half-embedded drug-eluting stent. Appl Math Comput 2015;266:968-981 View Article PubMed/NCBI
  22. Mandal AP, Mandal PK. Computational modelling of three-phase stent-based delivery. J Explor Res Pharmacol 2017;2:31-40 View Article PubMed/NCBI
  23. O’Brien C, Kolachalama V, Barber T, Simmons A, Edelman E. Impact of flow pulsatility on arterial drug distribution in stent-based therapy. J Control Release 2013;168:115-124 View Article PubMed/NCBI
  24. Saha R, Sarifuddin, Misra JC, Mandal PK. Impact of luminal flow on mass transport through coronary arteries: a study relevant to drug-eluting stent. Int J Math Comput 2016;27:40-58 View Article PubMed/NCBI
  25. Saltzman WM. Drug delivery: engineering principles for drug therapy. USA: Oxford University Press; 2001 View Article PubMed/NCBI
  26. Satterfield CN. Mass transfer in heterogeneous catalysis. The MIT Press; 1970 View Article PubMed/NCBI
  27. 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
  28. 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
  29. Sarifuddin, Mandal PK. Effect of diffusivity on the transport of drug eluted from drug-eluting stent. Int J Appl Comput Math 2016;2:291-301 View Article PubMed/NCBI
  30. Kolandaivelu K, O’Brien CC, Shazly T, Edelman ER, Kolachalama VB. Enhancing physiologic simulations using supervised learning on coarse mesh solutions. J R Soc Interface 2015;12:20141073 View Article PubMed/NCBI
  31. Pontrelli G, De Monte F. Modeling of mass dynamics in arterial drug-eluting stents. Journal of Porous Media 2009;12:19-28 View Article PubMed/NCBI
  32. 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
  33. Mongrain R, Faik I, Leask RL, Rodes-Cabau J, Larose E, Bertrand OF. Effects of diffusion coefficients and struts apposition using numerical simulations for drug eluting coronary stents. J Biomech Eng 2007;129:733-742 View Article PubMed/NCBI
  34. Moore JE, Ku DN. Pulsatile velocity measurements in a model of the human abdominal aorta under resting conditions. J Biomech Eng 1994;116:337-346 View Article PubMed/NCBI
  35. 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
  36. Khakpour M, Vafai K. Critical assessment of arterial transport models. Int J Heat Mass Transf 2008;51:807-822 View Article PubMed/NCBI
  37. O’Connell BM, Walsh MT. Biomedical Science, Engineering and Technology. London: InTechOpen; 2012, 301-324 View Article PubMed/NCBI
  38. Kamath KR, Barry JJ, Miller KM. The taxus drug-eluting stent: A new paradigm in controlled drug delivery. Adv Drug Deliv Rev 2006;58:412-436 View Article PubMed/NCBI
  39. Creel CJ, Lovich MA, Edelman ER. Arterial paclitaxel distribution and deposition. Circ Res 2000;86:879-884 View Article PubMed/NCBI
  40. Zhang F, Fath M, Marks R, Linhardt RJ. A highly stable covalent conjugated heparin biochip for heparin–protein interaction studies. Anal Biochem 2002;304:271-273 View Article PubMed/NCBI
  41. Zhu X, Pack DW, Braatz RD. Modelling intravascular delivery from drug-eluting stents with biodurable coating: investigation of anisotropic vascular drug diffusivity and arterial drug distribution. Comput Methods Biomech Biomed Engin 2014;17:187-198 View Article PubMed/NCBI
  42. 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
  43. McGinty S, McKee S, McCormick C, Wheel M. Release mechanism and parameter estimation in drug-eluting stent systems: analytical solutions of drug release and tissue transport. Math Med Biol 2015;32:163-186 View Article PubMed/NCBI
  44. Wessely R, Schömig A, Kastrati A. Sirolimus and paclitaxel on polymer-based drug-eluting stents: similar but different. J Am Coll Cardiol 2006;47:708-714 View Article PubMed/NCBI