Computational Model of Drug-Coated Balloon Delivery in a Patient-Specific Arterial Vessel with Heterogeneous Tissue Composition

Balloon angioplasty followed by local delivery of antiproliferative drugs to target tissue is increasingly being considered for the treatment of obstructive arterial disease, and yet there is much to appreciate regarding pharmacokinetics in arteries of non-uniform disease. We developed a computational model capable of simulating drug-coated balloon delivery to arteries of heterogeneous tissue composition comprising healthy tissue, as well as regions of fibrous, fibro-fatty, calcified and necrotic core lesions. Image processing using an unsupervised clustering technique was used to reconstruct an arterial geometry from a single, patient-specific color image obtained from intravascular ultrasound-derived virtual histology. Transport of free drug was modeled using a time-dependent reaction-diffusion model and the bound, immobilized drug using the time-dependent reaction equation. The governing equations representing the transport of free as well as bound drug along with a set of initial settings and boundary conditions were solved numerically using an explicit finite difference scheme that satisfied the Courant-Friedrichs-Lewy stability criterion. Our results support previous findings related to the transport and binding of drug in arteries where tissue retention is strongly dependent on local pharmacologic properties. Additionally, modeling results indicate that non-uniform disease composition leads to heterogeneous arterial drug distribution patterns, although further validation using animal studies is required to fully appreciate pharmacokinetics in disease-laden arteries.