We describe our fully-coupled 3D multiscale model of in-stent restenosis, with blood flow simulations coupled to smooth muscle cell proliferation, and report results of numerical simulations performed with this model. This novel model is based on several previously reported 2D models. We study the effects of various parameters on the process of restenosis and compare with in vivo porcine data where we observe good qualitative agreement. We study the effects of stent deployment depth (and related injury score), reendothelization speed, and simulate the effect of stent width. Also we demonstrate that we are now capable to simulate restenosis in real-sized (18 mm long, 2.8 mm wide) vessel geometries.

BackgroundCoronary artery restenosis is an important side effect of percutaneous coronary intervention. Computational models can be used to better understand this process. We report on an approach for validation of an in silico 3D model of in-stent restenosis in porcine coronary arteries and illustrate this approach by comparing the modelling results to in vivo data for 14 and 28 days post-stenting.MethodsThis multiscale model includes single-scale models for stent deployment, blood flow and tissue growth in the stented vessel, including smooth muscle cell (SMC) proliferation and extracellular matrix (ECM) production. The validation procedure uses data from porcine in vivo experiments, by simulating stent deployment using stent geometry obtained from micro computed tomography (micro-CT) of the stented vessel and directly comparing the simulation results of neointimal growth to histological sections taken at the same locations.ResultsMetrics for comparison are per-strut neointimal thickness and per-section neointimal area. The neointimal area predicted by the model demonstrates a good agreement with the detailed experimental data. For 14 days post-stenting the relative neointimal area, averaged over all vessel sections considered, was 20 ± 3% in vivo and 22 ± 4% in silico. For 28 days, the area was 42 ± 3% in vivo and 41 ± 3% in silico.ConclusionsThe approach presented here provides a very detailed, location-specific, validation methodology for in silico restenosis models. The model was able to closely match both histology datasets with a single set of parameters. Good agreement was obtained for both the overall amount of neointima produced and the local distribution. It should be noted that including vessel curvature and ECM production in the model was paramount to obtain a good agreement with the experimental data.Electronic supplementary materialThe online version of this article (10.1007/s13239-019-00431-4) contains supplementary material, which is available to authorized users.