Plaque heterogeneity and the spatial distributions of its components dictate drug-coated balloon therapy

Drug-coated balloon (DCB) angioplasty is one of the potential approaches to alleviating in-stent restenosis and treating peripheral artery disease. An in-silico model has been developed for sirolimus drug eluted from an inflated balloon in a patient-specific arterial cross-section consisting of fibrous tissue, fibrofatty tissue, dense calcium, necrotic core, and healthy tissue. The convection-diffusion-reaction equation represents the transport of drug, while drug binding, both specific and non-specific, can be modelled as a reaction process. The Brinkman equations describe the interstitial flow in porous tissue. An image processing technique is leveraged for reconstructing the computational domain. The Marker and Cell, and Immersed Boundary Methods are used to solve the set of governing equations. The no-flux interface condition and convection do amplify the tissue content, and the regions of dense calcium and necrotic core limited to or extremely close to the interface pose a clinical threat to DCB therapy. Simulations predict the effects of the positioning and clustering of plaque components in the domain. This study demands extensive intravascular ultrasound-derived virtual histology (VH-IVUS) imaging to understand the plaque morphology and determine the relative positions of different plaque compositions about the lumen-tissue interface, which have a significant impact on arterial pharmacokinetics.