Enhancing Cardiovascular Risk Stratification: Radiomics of Coronary Plaque and Perivascular Adipose Tissue - Current Insights and Future Perspectives.

Radiomics, the quantitative extraction and mining of features from radiological images, has recently emerged as a promising source of non-invasive image-based cardiovascular biomarkers, potentially revolutionizing diagnostics and risk assessment. This review explores its application within coronary plaques and pericoronary adipose tissue, particularly focusing on plaque characterization and cardiac events prediction. By shedding light on the current state-of-the-art, achievements, and prospective avenues, this review contributes to a deeper understanding of the evolving landscape of radiomics in the context of coronary arteries. Finally, open challenges and existing gaps are emphasized to underscore the need for future efforts aimed at ensuring the robustness and reliability of radiomics studies, facilitating their clinical translation.

Effect of variability of mechanical properties on the predictive capabilities of vulnerable coronary plaques.

BACKGROUND AND OBJECTIVE: Coronary plaque rupture is a precipitating event responsible for two thirds of myocardial infarctions. Currently, the risk of plaque rupture is computed based on demographic, clinical, and image-based adverse features. However, using these features the absolute event rate per single higher-risk lesion remains low. This work studies the power of a novel framework based on biomechanical markers accounting for material uncertainty to stratify vulnerable and non-vulnerable coronary plaques. METHODS: Virtual histology intravascular ultrasounds from 55 patients, 29 affected by acute coronary syndrome and 26 affected by stable angina pectoris, were included in this study. Two-dimensional vessel cross-sections for finite element modeling (10 sections per plaque) incorporating plaque structure (medial tissue, loose matrix, lipid core and calcification) were reconstructed. A Montecarlo finite element analysis was performed on each section to account for material variability on three biomechanical markers: peak plaque structural stress at diastolic and systolic pressure, and peak plaque stress difference between systolic and diastolic pressures, together with the luminal pressure. Machine learning decision tree classifiers were trained on 75% of the dataset and tested on the remaining 25% with a combination of feature selection techniques. Performance against classification trees based on geometric markers (i.e., luminal, external elastic membrane and plaque areas) was also performed. RESULTS: Our results indicate that the plaque structural stress outperforms the classification capacity of the combined geometric markers only (0.82 vs 0.51 area under curve) when accounting for uncertainty in material parameters. Furthermore, the results suggest that the combination of the peak plaque structural stress at diastolic and systolic pressures with the maximum plaque structural stress difference between systolic and diastolic pressures together with the systolic pressure and the diastolic to systolic pressure gradient is a robust classifier for coronary plaques when the intrinsic variability in material parameters is considered (area under curve equal to [0.91-0.93]). CONCLUSION: In summary, our results emphasize that peak plaque structural stress in combination with the patient's luminal pressure is a potential classifier of plaque vulnerability as it independently considers stress in all directions and incorporates total geometric and compositional features of atherosclerotic plaques.

Impact of Tissue Damage and Hemodynamics on Restenosis Following Percutaneous Transluminal Angioplasty: A Patient-Specific Multiscale Model.

Multiscale agent-based modeling frameworks have recently emerged as promising mechanobiological models to capture the interplay between biomechanical forces, cellular behavior, and molecular pathways underlying restenosis following percutaneous transluminal angioplasty (PTA). However, their applications are mainly limited to idealized scenarios. Herein, a multiscale agent-based modeling framework for investigating restenosis following PTA in a patient-specific superficial femoral artery (SFA) is proposed. The framework replicates the 2-month arterial wall remodeling in response to the PTA-induced injury and altered hemodynamics, by combining three modules: (i) the PTA module, consisting in a finite element structural mechanics simulation of PTA, featuring anisotropic hyperelastic material models coupled with a damage formulation for fibrous soft tissue and the element deletion strategy, providing the arterial wall damage and post-intervention configuration, (ii) the hemodynamics module, quantifying the post-intervention hemodynamics through computational fluid dynamics simulations, and (iii) the tissue remodeling module, based on an agent-based model of cellular dynamics. Two scenarios were explored, considering balloon expansion diameters of 5.2 and 6.2 mm. The framework captured PTA-induced arterial tissue lacerations and the post-PTA arterial wall remodeling. This remodeling process involved rapid cellular migration to the PTA-damaged regions, exacerbated cell proliferation and extracellular matrix production, resulting in lumen area reduction up to 1-month follow-up. After this initial reduction, the growth stabilized, due to the resolution of the inflammatory state and changes in hemodynamics. The similarity of the obtained results to clinical observations in treated SFAs suggests the potential of the framework for capturing patient-specific mechanobiological events occurring after PTA intervention.

MRI radiomics in head and neck cancer from reproducibility to combined approaches.

The clinical applicability of radiomics in oncology depends on its transferability to real-world settings. However, the absence of standardized radiomics pipelines combined with methodological variability and insufficient reporting may hamper the reproducibility of radiomic analyses, impeding its translation to clinics. This study aimed to identify and replicate published, reproducible radiomic signatures based on magnetic resonance imaging (MRI), for prognosis of overall survival in head and neck squamous cell carcinoma (HNSCC) patients. Seven signatures were identified and reproduced on 58 HNSCC patients from the DB2Decide Project. The analysis focused on: assessing the signatures' reproducibility and replicating them by addressing the insufficient reporting; evaluating their relationship and performances; and proposing a cluster-based approach to combine radiomic signatures, enhancing the prognostic performance. The analysis revealed key insights: (1) despite the signatures were based on different features, high correlations among signatures and features suggested consistency in the description of lesion properties; (2) although the uncertainties in reproducing the signatures, they exhibited a moderate prognostic capability on an external dataset; (3) clustering approaches improved prognostic performance compared to individual signatures. Thus, transparent methodology not only facilitates replication on external datasets but also advances the field, refining prognostic models for potential personalized medicine applications.

Developing a robust two-step machine learning multiclassification pipeline to predict primary site in head and neck carcinoma from lymph nodes.

This study aimed to develop a robust multiclassification pipeline to determine the primary tumor location in patients with head and neck carcinoma of unknown primary using radiomics and machine learning techniques. The dataset included 400 head and neck cancer patients with primary tumor in oropharynx, OPC (n = 162), nasopharynx, NPC (n = 137), oral cavity, OC (n = 63), larynx and hypopharynx, HL (n = 38). Two radiomic-based multiclassification pipelines (P1 and P2) were developed. P1 consisted in a direct identification of the primary sites, whereas P2 was based on a two-step approach: in the first step, the number of classes was reduced by merging the two minority classes which were reclassified in the second step. Diverse correlation thresholds (0.75, 0.80, 0.85), feature selection methods (sequential forwards/backwards selection, sequential floating forward selection, neighborhood component analysis and minimum redundancy maximum relevance), and classification models (neural network, decision tree, naïve Bayes, bagged trees and support vector machine) were assessed. P2 outperformed P1, with the best results obtained with the support vector machine classifier including radiomic and clinical features (accuracies of 75.3 % (HL), 75.4 % (OC), 71.3 % (OPC), 92.9 % (NPC)). These results indicate that the two-step multiclassification pipeline integrating radiomics and clinical information is a promising approach to predict the tumor site of unknown primary.

An agent-based model of cardiac allograft vasculopathy: toward a better understanding of chronic rejection dynamics.

Cardiac allograft vasculopathy (CAV) is a coronary artery disease affecting 50% of heart transplant (HTx) recipients, and it is the major cause of graft loss. CAV is driven by the interplay of immunological and non-immunological factors, setting off a cascade of events promoting endothelial damage and vascular dysfunction. The etiology and evolution of tissue pathology are largely unknown, making disease management challenging. So far, in vivo models, mostly mouse-based, have been widely used to study CAV, but they are resource-consuming, pose many ethical issues, and allow limited investigation of time points and important biomechanical measurements. Recently, agent-based models (ABMs) proved to be valid computational tools for deciphering mechanobiological mechanisms driving vascular adaptation processes at the cell/tissue level, augmenting cost-effective in vivo lab-based experiments, at the same time guaranteeing richness in observation time points and low consumption of resources. We hypothesize that integrating ABMs with lab-based experiments can aid in vivo research by overcoming those limitations. Accordingly, this work proposes a bidimensional ABM of CAV in a mouse coronary artery cross-section, simulating the arterial wall response to two distinct stimuli: inflammation and hemodynamic disturbances, the latter considered in terms of low wall shear stress (WSS). These stimuli trigger i) inflammatory cell activation and ii) exacerbated vascular cell activities. Moreover, an extensive analysis was performed to investigate the ABM sensitivity to the driving parameters and inputs and gain insights into the ABM working mechanisms. The ABM was able to effectively replicate a 4-week CAV initiation and progression, characterized by lumen area decrease due to progressive intimal thickening in regions exposed to high inflammation and low WSS. Moreover, the parameter and input sensitivity analysis highlighted that the inflammatory-related events rather than the WSS predominantly drive CAV, corroborating the inflammatory nature of the vasculopathy. The proof-of-concept model proposed herein demonstrated its potential in deepening the pathology knowledge and supporting the in vivo analysis of CAV.

Investigating the effect of drug release on in-stent restenosis: A hybrid continuum - agent-based modelling approach.

BACKGROUND AND OBJECTIVE: In-stent restenosis (ISR) following percutaneous coronary intervention with drug-eluting stent (DES) implantation remains an unresolved issue, with ISR rates up to 10%. The use of antiproliferative drugs on DESs has significantly reduced ISR. However, a complete knowledge of the mechanobiological processes underlying ISR is still lacking. Multiscale agent-based modelling frameworks, integrating continuum- and agent-based approaches, have recently emerged as promising tools to decipher the mechanobiological events driving ISR at different spatiotemporal scales. However, the integration of sophisticated drug models with an agent-based model (ABM) of ISR has been under-investigated. The aim of the present study was to develop a novel multiscale agent-based modelling framework of ISR following DES implantation. METHODS: The framework consisted of two bi-directionally coupled modules, namely (i) a drug transport module, simulating drug transport through a continuum-based approach, and (ii) a tissue remodelling module, simulating cellular dynamics through an ABM. Receptor saturation (RS), defined as the fraction of target receptors saturated with drug, is used to mediate cellular activities in the ABM, since RS is widely regarded as a measure of drug efficacy. Three studies were performed to investigate different scenarios in terms of drug mass (DM), drug release profiles (RP), coupling schemes and idealized vs. patient-specific artery geometries. RESULTS: The studies demonstrated the versatility of the framework and enabled exploration of the sensitivity to different settings, coupling modalities and geometries. As expected, changes in the DM, RP and coupling schemes illustrated a variation in RS over time, in turn affecting the ABM response. For example, combined small DM - fast RP led to similar ISR degrees as high DM - moderate RP (lumen area reduction of ∼13/17% vs. ∼30% without drug). The use of a patient-specific geometry with non-equally distributed struts resulted in a heterogeneous RS map, but did not remarkably impact the ABM response. CONCLUSION: The application to a patient-specific geometry highlights the potential of the framework to address complex realistic scenarios and lays the foundations for future research, including calibration and validation on patient datasets and the investigation of the effects of different plaque composition on the arterial response to DES.

MRI-based radiomic prognostic signature for locally advanced oral cavity squamous cell carcinoma: development, testing and comparison with genomic prognostic signatures.

BACKGROUND: . At present, the prognostic prediction in advanced oral cavity squamous cell carcinoma (OCSCC) is based on the tumor-node-metastasis (TNM) staging system, and the most used imaging modality in these patients is magnetic resonance image (MRI). With the aim to improve the prediction, we developed an MRI-based radiomic signature as a prognostic marker for overall survival (OS) in OCSCC patients and compared it with published gene expression signatures for prognosis of OS in head and neck cancer patients, replicated herein on our OCSCC dataset. METHODS: For each patient, 1072 radiomic features were extracted from T1 and T2-weighted MRI (T1w and T2w). Features selection was performed, and an optimal set of five of them was used to fit a Cox proportional hazard regression model for OS. The radiomic signature was developed on a multi-centric locally advanced OCSCC retrospective dataset (n = 123) and validated on a prospective cohort (n = 108). RESULTS: The performance of the signature was evaluated in terms of C-index (0.68 (IQR 0.66-0.70)), hazard ratio (HR 2.64 (95% CI 1.62-4.31)), and high/low risk group stratification (log-rank p < 0.001, Kaplan-Meier curves). When tested on a multi-centric prospective cohort (n = 108), the signature had a C-index of 0.62 (IQR 0.58-0.64) and outperformed the clinical and pathologic TNM stage and six out of seven gene expression prognostic signatures. In addition, the significant difference of the radiomic signature between stages III and IVa/b in patients receiving surgery suggests a potential association of MRI features with the pathologic stage. CONCLUSIONS: Overall, the present study suggests that MRI signatures, containing non-invasive and cost-effective remarkable information, could be exploited as prognostic tools.

Combining deep learning and machine learning for the automatic identification of hip prosthesis failure: Development, validation and explainability analysis.

AIM: Revision hip arthroplasty has a less favorable outcome than primary total hip arthroplasty and an understanding of the timing of total hip arthroplasty failure may be helpful. The aim of this study is to develop a combined deep learning (DL) and machine learning (ML) approach to automatically detect hip prosthetic failure from conventional plain radiographs. METHODS: Two cohorts of patients (of 280 and 352 patients) were included in the study, for model development and validation, respectively. The analysis was based on one antero-posterior and one lateral radiographic view obtained from each patient during routine post-surgery follow-up. After pre-processing, three images were obtained: the original image, the acetabulum image and the stem image. These images were analyzed through convolutional neural networks aiming to predict prosthesis failure. Deep features of the three images were extracted for each model and two feature-based pipelines were developed: one utilizing only the features of the original image (original image pipeline) and the other concatenating the features of the three images (3-image pipeline). The obtained features were either used directly or reduced through principal component analysis. Both support vector machine (SVM) and random forest (RF) classifiers were considered for each pipeline. RESULTS: The SVM applied to the 3-image pipeline provided the best performance, with an accuracy of 0.958 ± 0.006 in the internal validation and an F1-score of 0.874 in the external validation set. The explainability analysis, besides identifying the features of the complete original images as the major contributor, highlighted the role of the acetabulum and stem images on the prediction. CONCLUSIONS: This study demonstrated the potentialities of the developed DL-ML procedure based on plain radiographs in the detection of the failure of the hip prosthesis.

Predicting 1-year in-stent restenosis in superficial femoral arteries through multiscale computational modelling.

In-stent restenosis in superficial femoral arteries (SFAs) is a complex, multi-factorial and multiscale vascular adaptation process whose thorough understanding is still lacking. Multiscale computational agent-based modelling has recently emerged as a promising approach to decipher mechanobiological mechanisms driving the arterial response to the endovascular intervention. However, the long-term arterial response has never been investigated with this approach, although being of fundamental relevance. In this context, this study investigates the 1-year post-operative arterial wall remodelling in three patient-specific stented SFA lesions through a fully coupled multiscale agent-based modelling framework. The framework integrates the effects of local haemodynamics and monocyte gene expression data on cellular dynamics through a bi-directional coupling of computational fluid dynamics simulations with an agent-based model of cellular activities. The framework was calibrated on the follow-up data at 1 month and 6 months of one stented SFA lesion and then applied to the other two lesions. The calibrated framework successfully captured (i) the high lumen area reduction occurring within the first post-operative month and (ii) the stabilization of the median lumen area from 1-month to 1-year follow-ups in all the stented lesions, demonstrating the potentialities of the proposed approach for investigating patient-specific short- and long-term responses to endovascular interventions.

Impact of nickel-titanium super-elastic material properties on the mechanical performance of self-expandable transcatheter aortic valves.

Self-expandable transcatheter aortic valves (TAVs) elastically resume their initial shape when implanted without the need for balloon inflation by virtue of the nickel-titanium (NiTi) frame super-elastic properties. Experimental findings suggest that NiTi mechanical properties can vary markedly because of a strong dependence on the chemical composition and processing operations. In this context, this study presents a computational framework to investigate the impact of the NiTi super-elastic material properties on the TAV mechanical performance. Finite element (FE) analyses of TAV implantation were performed considering two different TAV frames and three idealized aortic root anatomies, evaluating the device mechanical response in terms of pullout force magnitude exerted by the TAV frame and peak maximum principal stress within the aortic root. The widely adopted NiTi constitute model by Auricchio and Taylor (1997) was used. A multi-parametric sensitivity analysis and a multi-objective optimization of the TAV mechanical performance were conducted in relation to the parameters of the NiTi constitutive model. The results highlighted that: five NiTi material model parameters (EA, σtLS, σtUS, σtUE and σcLS) are significantly correlated with the FE outputs; the TAV frame geometry and aortic root anatomy have a marginal effect on the level of influence of each NiTi material parameter; NiTi alloy candidates with pareto-optimal characteristics in terms of TAV mechanical performance can be successfully identified. In conclusion, the proposed computational framework supports the TAV design phase, providing information on the relationship between the super-elastic behavior of the supplied NiTi alloys and the device mechanical response.

Multiscale agent-based modeling of restenosis after percutaneous transluminal angioplasty: Effects of tissue damage and hemodynamics on cellular activity.

BACKGROUND: Restenosis following percutaneous transluminal angioplasty (PTA) in femoral arteries is a major cause of failure of the revascularization procedure. The arterial wall response to PTA is driven by multifactorial, multiscale processes, whose complete understanding is lacking. Multiscale agent-based modeling frameworks, simulating the network of mechanobiological events at cell-tissue scale, can contribute to decipher the pathological pathways of restenosis. In this context, the present study proposes a fully-automated multiscale agent-based modeling framework simulating the arterial wall remodeling due to the wall damage provoked by PTA and to the altered hemodynamics in the post-operative months. METHODS: The framework, applied to an idealized femoral artery model, integrated: (i) a PTA module (i.e., structural mechanics simulation), computing the post-PTA arterial morphology and the PTA-induced damage (ii) a hemodynamics module (i.e., computational fluid dynamics simulations), quantifying the near-wall hemodynamics, and (iii) a tissue remodeling module simulating cellular activities through an agent-based model. RESULTS: The framework was able to capture relevant features of the 3-month arterial wall response to PTA, namely (i) the impact of the PTA-induced damage and altered hemodynamics on arterial wall remodeling, including the local intimal growth at the most susceptible regions (i.e., elevated damage levels and low wall shear stress), (ii) the lumen area temporal trend resulting from the interaction of the two inputs, and (iii) a 3-month lumen area restenosis of ∼25%, in accordance with clinical evidence. CONCLUSIONS: The overall results demonstrated the framework potentiality in capturing mechanobiological processes underlying restenosis, thus fostering future application to patient-specific scenarios.

A predictive multiscale model of in-stent restenosis in femoral arteries: linking haemodynamics and gene expression with an agent-based model of cellular dynamics.

In-stent restenosis (ISR) is a maladaptive inflammatory-driven response of femoral arteries to percutaneous transluminal angioplasty and stent deployment, leading to lumen re-narrowing as consequence of excessive cellular proliferative and synthetic activities. A thorough understanding of the underlying mechanobiological factors contributing to ISR is still lacking. Computational multiscale models integrating both continuous- and agent-based approaches have been identified as promising tools to capture key aspects of the complex network of events encompassing molecular, cellular and tissue response to the intervention. In this regard, this work presents a multiscale framework integrating the effects of local haemodynamics and monocyte gene expression data on cellular dynamics to simulate ISR mechanobiological processes in a patient-specific model of stented superficial femoral artery. The framework is based on the coupling of computational fluid dynamics simulations (haemodynamics module) with an agent-based model (ABM) of cellular activities (tissue remodelling module). Sensitivity analysis and surrogate modelling combined with genetic algorithm optimization were adopted to explore the model behaviour and calibrate the ABM parameters. The proposed framework successfully described the patient lumen area reduction from baseline to one-month follow-up, demonstrating the potential capabilities of this approach in predicting the short-term arterial response to the endovascular procedure.

Superficial femoral artery stenting: Impact of stent design and overlapping on the local hemodynamics.

BACKGROUND: Superficial femoral arteries (SFAs) treated with self-expanding stents are widely affected by in-stent restenosis (ISR), especially in case of long lesions and multiple overlapping devices. The altered hemodynamics provoked by the stent is considered as a promoting factor of ISR. In this context, this work aims to analyze the impact of stent design and stent overlapping on patient-specific SFA hemodynamics. METHODS: Through a morphing technique, single or multiple stents were virtually implanted within two patient-specific, post-operative SFA models reconstructed from computed tomography. The stented domains were used to perform computational fluid dynamics simulations, quantifying wall shear stress (WSS) based descriptors including time-averaged WSS (TAWSS), oscillatory shear index (OSI), transverse WSS (transWSS), and WSS ratio (WSSRATIO). Four stent designs (three laser-cut - EverFlex, Zilver and S.M.A.R.T. - and one prototype braided stent), and three typical clinical scenarios accounting for different order of stent implantation and overlapping length were compared. RESULTS: The main hemodynamic differences were found between the two types of stent designs (i.e. laser-cut vs. braided stents). The braided stent presented lower median transWSS and higher median WSSRATIO than the laser-cut stents (p < 0.0001). The laser-cut stents presented comparable WSS-based descriptor values, except for the Zilver, exhibiting a median TAWSS ∼30% higher than the other stents. Stent overlapping provoked an abrupt alteration of the WSS-based descriptors. The overlapping length, rather than the order of stent implantation, highly and negatively impacted the hemodynamics. CONCLUSION: The proposed computational workflow compared different SFA stent designs and stent overlapping configurations, highlighting those providing the most favorable hemodynamic conditions.

Coronary Artery Stenting Affects Wall Shear Stress Topological Skeleton.

Despite the important advancements in the stent technology for the treatment of diseased coronary arteries, major complications still affect the postoperative long-term outcome. The stent-induced flow disturbances, and especially the altered wall shear stress (WSS) profile at the strut level, play an important role in the pathophysiological mechanisms leading to stent thrombosis (ST) and in-stent restenosis (ISR). In this context, the analysis of the WSS topological skeleton is gaining more and more interest by extending the current understanding of the association between local hemodynamics and vascular diseases. This study aims to analyze the impact that a deployed coronary stent has on the WSS topological skeleton. Computational fluid dynamics (CFD) simulations were performed in three stented human coronary artery geometries reconstructed from clinical images. The selected cases presented stents with different designs (i.e., two contemporary drug-eluting stents and one bioresorbable scaffold) and included regions with stent malapposition or overlapping. A recently proposed Eulerian-based approach was applied to analyze the WSS topological skeleton features. The results highlighted that the presence of single or multiple stents within a coronary artery markedly impacts the WSS topological skeleton. In particular, repetitive patterns of WSS divergence were observed at the luminal surface, highlighting a WSS contraction action exerted proximal to the stent struts and a WSS expansion action distal to the stent struts. This WSS action pattern was independent from the stent design. In conclusion, these findings could contribute to a deeper understanding of the hemodynamics-driven processes underlying ST and ISR.

Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models.

The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.

3D modelling of drug-coated balloons for the treatment of calcified superficial femoral arteries.

BACKGROUND/OBJECTIVES: Drug-coated balloon therapy for diseased superficial femoral arteries remains controversial. Despite its clinical relevance, only a few computational studies based on simplistic two-dimensional models have been proposed to investigate this endovascular therapy to date. This work addresses the aforementioned limitation by analyzing the drug transport and kinetics occurring during drug-coated balloon deployment in a three-dimensional geometry. METHODS: An idealized three-dimensional model of a superficial femoral artery presenting with a calcific plaque and treated with a drug-coated balloon was created to perform transient mass transport simulations. To account for the transport of drug (i.e. paclitaxel) released by the device, a diffusion-reaction equation was implemented by describing the drug bound to specific intracellular receptors through a non-linear, reversible reaction. The following features concerning procedural aspects, pathologies and modelling assumptions were investigated: (i) balloon application time (60-180 seconds); (ii) vessel wall composition (healthy vs. calcified wall); (iii) sequential balloon application; and (iv) drug wash-out by the blood stream vs. coating retention, modeled as exponential decay. RESULTS: The balloon inflation time impacted both the free and specifically-bound drug concentrations in the vessel wall. The vessel wall composition highly affected the drug concentrations. In particular, the specifically-bound drug concentration was four orders of magnitude lower in the calcific compared with healthy vessel wall portions, primarily as a result of reduced drug diffusion. The sequential application of two drug-coated balloons led to modest differences (~15%) in drug concentration immediately after inflation, which became negligible within 10 minutes. The retention of the balloon coating increased the drug concentration in the vessel wall fourfold. CONCLUSIONS: The overall findings suggest that paclitaxel kinetics may be affected not only by the geometrical and compositional features of the vessel treated with the drug-coated balloon, but also by balloon design characteristics and procedural aspects that should be carefully considered.

In-Stent Restenosis Progression in Human Superficial Femoral Arteries: Dynamics of Lumen Remodeling and Impact of Local Hemodynamics.

In-stent restenosis (ISR) represents a major drawback of stented superficial femoral arteries (SFAs). Motivated by the high incidence and limited knowledge of ISR onset and development in human SFAs, this study aims to (i) analyze the lumen remodeling trajectory over 1-year follow-up period in human stented SFAs and (ii) investigate the impact of altered hemodynamics on ISR initiation and progression. Ten SFA lesions were reconstructed at four follow-ups from computed tomography to quantify the lumen area change occurring within 1-year post-intervention. Patient-specific computational fluid dynamics simulations were performed at each follow-up to relate wall shear stress (WSS) based descriptors with lumen remodeling. The largest lumen remodeling was found in the first post-operative month, with slight regional-specific differences (larger inward remodeling in the fringe segments, p < 0.05). Focal re-narrowing frequently occurred after 6 months. Slight differences in the lumen area change emerged between long and short stents, and between segments upstream and downstream from stent overlapping portions, at specific time intervals. Abnormal patterns of multidirectional WSS were associated with lumen remodeling within 1-year post-intervention. This longitudinal study gave important insights into the dynamics of ISR and the impact of hemodynamics on ISR progression in human SFAs.

Baseline local hemodynamics as predictor of lumen remodeling at 1-year follow-up in stented superficial femoral arteries.

In-stent restenosis (ISR) is the major drawback of superficial femoral artery (SFA) stenting. Abnormal hemodynamics after stent implantation seems to promote the development of ISR. Accordingly, this study aims to investigate the impact of local hemodynamics on lumen remodeling in human stented SFA lesions. Ten SFA models were reconstructed at 1-week and 1-year follow-up from computed tomography images. Patient-specific computational fluid dynamics simulations were performed to relate the local hemodynamics at 1-week, expressed in terms of time-averaged wall shear stress (TAWSS), oscillatory shear index and relative residence time, with the lumen remodeling at 1-year, quantified as the change of lumen area between 1-week and 1-year. The TAWSS was negatively associated with the lumen area change (ρ = - 0.75, p = 0.013). The surface area exposed to low TAWSS was positively correlated with the lumen area change (ρ = 0.69, p = 0.026). No significant correlations were present between the other hemodynamic descriptors and lumen area change. The low TAWSS was the best predictive marker of lumen remodeling (positive predictive value of 44.8%). Moreover, stent length and overlapping were predictor of ISR at follow-up. Despite the limited number of analyzed lesions, the overall findings suggest an association between abnormal patterns of WSS after stenting and lumen remodeling.

A fully coupled computational fluid dynamics - agent-based model of atherosclerotic plaque development: Multiscale modeling framework and parameter sensitivity analysis.

BACKGROUND: Peripheral Artery Disease (PAD) is an atherosclerotic disorder that leads to impaired lumen patency through intimal hyperplasia and the build-up of plaques, mainly localized in areas of disturbed flow. Computational models can provide valuable insights in the pathogenesis of atherosclerosis and act as a predictive tool to optimize current interventional techniques. Our hypothesis is that a reliable predictive model must include the atherosclerosis development history. Accordingly, we developed a multiscale modeling framework of atherosclerosis that replicates the hemodynamic-driven arterial wall remodeling and plaque formation. METHODS: The framework was based on the coupling of Computational Fluid Dynamics (CFD) simulations with an Agent-Based Model (ABM). The CFD simulation computed the hemodynamics in a 3D artery model, while 2D ABMs simulated cell, Extracellular Matrix (ECM) and lipid dynamics in multiple vessel cross-sections. A sensitivity analysis was also performed to evaluate the oscillation of the ABM output to variations in the inputs and to identify the most influencing ABM parameters. RESULTS: Our multiscale model qualitatively replicated both the physiologic and pathologic arterial configuration, capturing histological-like features. The ABM outputs were mostly driven by cell and ECM dynamics, largely affecting the lumen area. A subset of parameters was found to affect the final lipid core size, without influencing cell/ECM or lumen area trends. CONCLUSION: The fully coupled CFD-ABM framework described atherosclerotic morphological and compositional changes triggered by a disturbed hemodynamics.