An in silico study on the role of smooth muscle cell migration in neointimal formation after coronary stenting

Exploring the Potential of Saphenous Vein Grafts Ex Vivo: A Model for Intimal Hyperplasia and Re-Endothelialization

Coronary artery bypass grafting (CABG) utilizing saphenous vein grafts (SVGs) stands as a fundamental approach to surgically treating coronary artery disease. However, the long-term success of CABG is often compromised by the development of intimal hyperplasia (IH) and subsequent graft failure. Understanding the mechanisms underlying this pathophysiology is crucial for improving graft patency and patient outcomes. Objectives: This study aims to explore the potential of an ex vivo model utilizing SVG to investigate IH and re-endothelialization. Methods: A thorough histological examination of 15 surplus SVG procured from CABG procedures at Hospital Canselor Tuanku Muhriz, Malaysia, was conducted to establish their baseline characteristics. Results: SVGs exhibited a mean diameter of 2.65 ± 0.93 mm with pre-existing IH averaging 0.42 ± 0.13 mm in thickness, alongside an observable lack of luminal endothelial cell lining. Analysis of extracellular matrix components, including collagen, elastin, and glycosaminoglycans, at baseline and after 7 days of ex vivo culture revealed no significant changes in collagen but demonstrated increased percentages of elastin and glycosaminoglycans. Despite unsuccessful attempts at re-endothelialization with blood outgrowth endothelial cells, the established ex vivo SVG IH model underscores the multifaceted nature of graft functionality and patency, characterized by IH presence, endothelial impairment, and extracellular matrix alterations post-CABG. Conclusions: The optimized ex vivo IH model provides a valuable platform for delving into the underlying mechanisms of IH formation and re-endothelialization of SVG. Further refinements are warranted, yet this model holds promise for future research aimed at enhancing graft durability and outcomes for CAD patients undergoing CABG.

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.

Computational modeling of in-stent restenosis: Pharmacokinetic and pharmacodynamic evaluation

Persistence of the pathology of in-stent restenosis even with the advent of drug-eluting stents warrants the development of highly resolved in silico models. These computational models assist in gaining insights into the transient biochemical and cellular mechanisms involved and thereby optimize the stent implantation parameters. Within this work, an already established fully-coupled Lagrangian finite element framework for modeling the restenotic growth is enhanced with the incorporation of endothelium-mediated effects and pharmacological influences of rapamycin-based drugs embedded in the polymeric layers of the current generation drug-eluting stents. The continuum mechanical description of growth is further justified in the context of thermodynamic consistency. Qualitative inferences are drawn from the model developed herein regarding the efficacy of the level of drug embedment within the struts as well as the release profiles adopted. The framework is then intended to serve as a tool for clinicians to tune the interventional procedures patient-specifically.

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.

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.

A novel co-culture assay to evaluate the effects of sympathetic innervation on vascular smooth muscle differentiation

Dedifferentiation of vascular smooth muscle cells (VSMCs) from a functional phenotype to an inverse synthetic phenotype is a symptom of cardiovascular disorders, such as atherosclerosis and hypertension. The sympathetic nervous system (SNS) is an essential regulator of the differentiation of vascular smooth muscle cells (VSMCs). In addition, numerous studies suggest that SNS also stimulates VSMCs to retain their contractile phenotype. However, the molecular mechanisms for this stimulation have not been thoroughly studied. In this study, we used a novel in vitro co-culture method to evaluate the effective cellular interactions and stimulatory effects of sympathetic neurons on the differentiation of VSMCs. We co-cultured rat neural-like pheochromocytoma cells (PC12) and rat aortic VSMCs with this method. Expression of VSMCs contractile genes, including smooth muscle actin (acta2), myosin heavy chain (myh11), elastin (eln), and smoothelin (smtn), were determined by quantitative real-time-PCR analysis as an indicator of VSMCs differentiation. Fold changes for specific contractile genes in VSMCs grown in vitro for seven days in the presence (innervated) and absence (non-innervated) of sympathetic neurons were 3.5 for acta2, 6.5 for myh11, 4.19 for eln, and 4 for smtn (normalized to Tata Binding Protein (TBP)). As a result, these data suggest that sympathetic innervation promotes VSMCs' contractile gene expression and also maintains VSMCs' functional phenotype.

An intricate interplay between stent drug dose and release rate dictates arterial restenosis

Since the introduction of percutaneous coronary intervention (PCI) for the treatment of obstructive coronary artery disease (CAD), patient outcomes have progressively improved. Drug eluting stents (DES) that employ anti-proliferative drugs to limit excess tissue growth following stent deployment have proved revolutionary. However, restenosis and a need for repeat revascularisation still occurs after DES use. Over the last few years, computational models have emerged that detail restenosis following the deployment of a bare metal stent (BMS), focusing primarily on contributions from mechanics and fluid dynamics. However, none of the existing models adequately account for spatiotemporal delivery of drug and the influence of this on the cellular processes that drive restenosis. In an attempt to fill this void, a novel continuum restenosis model coupled with spatiotemporal drug delivery is presented. Our results indicate that the severity and time-course of restenosis is critically dependent on the drug delivery strategy. Specifically, we uncover an intricate interplay between initial drug loading, drug release rate and restenosis, indicating that it is not sufficient to simply ramp-up the drug dose or prolong the time course of drug release to improve stent efficacy. Our model also shows that the level of stent over-expansion and stent design features, such as inter-strut spacing and strut thickness, influence restenosis development, in agreement with trends observed in experimental and clinical studies. Moreover, other critical aspects of the model which dictate restenosis, including the drug binding site density are investigated, where comparisons are made between approaches which assume this to be either constant or proportional to the number of smooth muscle cells (SMCs). Taken together, our results highlight the necessity of incorporating these aspects of drug delivery in the pursuit of optimal DES design.

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.

Mechano-chemo-biological Computational Models for Arteries in Health, Disease and Healing: From Tissue Remodelling to Drug-eluting Devices

This review aims to highlight urgent priorities for the computational biomechanics community in the framework of mechano-chemo-biological models. Recent approaches, promising directions and open challenges on the computational modelling of arterial tissues in health and disease are introduced and investigated, together with in silico approaches for the analysis of drug-eluting stents that promote pharmacological-induced healing. The paper addresses a number of chemo-biological phenomena that are generally neglected in biomechanical engineering models but are most likely instrumental for the onset and the progression of arterial diseases. An interdisciplinary effort is thus encouraged for providing the tools for an effective in silico insight into medical problems. An integrated mechano-chemo-biological perspective is believed to be a fundamental missing piece for crossing the bridge between computational engineering and life sciences, and for bringing computational biomechanics into medical research and clinical practice.

A homogenized constrained mixture model of restenosis and vascular remodelling after balloon angioplasty

Restenosis is one of the main adverse effects of the treatment of atherosclerosis through balloon angioplasty or stenting. During the intervention, the arterial wall is overstretched, causing a cascade of cellular events and subsequent neointima formation. This mechanical stimulus and its mechanobiological effects can be reproduced in biomechanical simulations. The aim of these models is to predict the long-term outcome of these procedures, to help increase the understanding of restenosis formation and to allow for in silico optimization of the treatment. We propose a predictive finite-element model of restenosis, using the homogenized constrained mixture modelling framework designed to model growth and remodelling in soft tissues. We compare the results with clinical observations in human coronary arteries and experimental findings in non-human primate models. We also explore the model's clinical relevance by testing its response to different balloon loads and to the use of drug-eluting balloons. The comparison of the results with experimental data shows the relevance of the model. We show its ability to predict both inward and outward remodelling as observed in vivo and we show the importance of an improved understanding of restenosis formation from a biomechanical point of view.

Do we really understand how drug eluted from stents modulates arterial healing?

The advent of drug-eluting stents (DES) has revolutionised the treatment of coronary artery disease. These devices, coated with anti-proliferative drugs, are deployed into stenosed or occluded vessels, compressing the plaque to restore natural blood flow, whilst simultaneously combating the evolution of restenotic tissue. Since the development of the first stent, extensive research has investigated how further advancements in stent technology can improve patient outcome. Mathematical and computational modelling has featured heavily, with models focussing on structural mechanics, computational fluid dynamics, drug elution kinetics and subsequent binding within the arterial wall; often considered separately. Smooth Muscle Cell (SMC) proliferation and neointimal growth are key features of the healing process following stent deployment. However, models which depict the action of drug on these processes are lacking. In this article, we start by reviewing current models of cell growth, which predominantly emanate from cancer research, and available published data on SMC proliferation, before presenting a series of mathematical models of varying complexity to detail the action of drug on SMC growth in vitro. Our results highlight that, at least for Sodium Salicylate and Paclitaxel, the current state-of-the-art nonlinear saturable binding model is incapable of capturing the proliferative response of SMCs across a range of drug doses and exposure times. Our findings potentially have important implications on the interpretation of current computational models and their future use to optimise and control drug release from DES and drug-coated balloons.

[Review of studies on the biomechanical modelling of the coupling effect between stent degradation and blood vessel remodeling]

The dynamic coupling of stent degradation and vessel remodeling can influence not only the structural morphology and material property of stent and vessel, but also the development of in-stent restenosis. The research achievements of biomechanical modelling and analysis of stent degradation and vessel remodeling were reviewed; several noteworthy research perspectives were addressed, a stent-vessel coupling model was developed based on stent damage function and vessel growth function, and then concepts of matching ratio and risk factor were established so as to evaluate the treatment effect of stent intervention, which may lay the scientific foundation for the structure design, mechanical analysis and clinical application of biodegradable stent.

Reversal of Endothelial Extracellular Vesicle-Induced Smooth Muscle Phenotype Transition by Hypercholesterolemia Stimulation: Role of NLRP3 Inflammasome Activation

Recent studies reported that vascular endothelial cells (ECs) secrete NLR family pyrin domain-containing 3 (NLRP3) inflammasome products such as interleukin-1β (IL-1β) via extracellular vesicles (EVs) under various pathological conditions. EVs represent one of the critical mechanisms mediating the cell-to-cell communication between ECs and vascular smooth muscle cells (VSMCs). However, whether or not the inflammasome-dependent EVs directly participate in the regulation of VSMC function remains unknown. In the present study, we found that in cultured carotid ECs, atherogenic stimulation by oxysterol 7-ketocholesterol (7-Ket) induced NLRP3 inflammasome formation and activation, reduced lysosome-multivesicular bodies (MVBs) fusion, and increased secretion of EVs that contain inflammasome product IL-1β. These EC-derived IL-1β-containing EVs promoted synthetic phenotype transition of co-cultured VSMCs, whereas EVs from unstimulated ECs have the opposite effects. Moreover, acid ceramidase (AC) deficiency or lysosome inhibition further exaggerated the 7-Ket-induced release of IL-1β-containing EVs in ECs. Using a Western diet (WD)-induced hypercholesterolemia mouse model, we found that endothelial-specific AC gene knockout mice (Asah1^fl/fl/EC^Cre) exhibited augmented WD-induced EV secretion with IL-1β and more significantly decreased the interaction of MVBs with lysosomes in the carotid arterial wall compared to their wild-type littermates (WT/WT). The endothelial AC deficiency in Asah1^fl/fl/EC^Cre mice also resulted in enhanced VSMC phenotype transition and accelerated neointima formation. Together, these results suggest that NLRP3 inflammasome-dependent IL-1β production during hypercholesterolemia promotes VSMC phenotype transition to synthetic status via EV machinery, which is controlled by lysosomal AC activity. Our findings provide novel mechanistic insights into understanding the pathogenic role of endothelial NLRP3 inflammasome in vascular injury through EV-mediated EC-to-VSMC regulation.

Computational modeling for cardiovascular tissue engineering: the importance of including cell behavior in growth and remodeling algorithms

Understanding cardiovascular growth and remodeling (G&R) is fundamental for designing robust cardiovascular tissue engineering strategies, which enable synthetic or biological scaffolds to transform into healthy living tissues after implantation. Computational modeling, particularly when integrated with experimental research, is key for advancing our understanding, predicting the in vivo evolution of engineered tissues, and efficiently optimizing scaffold designs. As cells are ultimately the drivers of G&R and known to change their behavior in response to mechanical cues, increasing efforts are currently undertaken to capture (mechano-mediated) cell behavior in computational models. In this selective review, we highlight some recent examples that are relevant in the context of cardiovascular tissue engineering and discuss the current and future biological and computational challenges for modeling cell-mediated G&R.

Mechanistic evaluation of long-term in-stent restenosis based on models of tissue damage and growth

Development and application of advanced mechanical models of soft tissues and their growth represent one of the main directions in modern mechanics of solids. Such models are increasingly used to deal with complex biomedical problems. Prediction of in-stent restenosis for patients treated with coronary stents remains a highly challenging task. Using a finite element method, this paper presents a mechanistic approach to evaluate the development of in-stent restenosis in an artery following stent implantation. Hyperelastic models with damage, verified with experimental results, are used to describe the level of tissue damage in arterial layers and plaque caused by such intervention. A tissue-growth model, associated with vessel damage, is adopted to describe the growth behaviour of a media layer after stent implantation. Narrowing of lumen diameter with time is used to quantify the development of in-stent restenosis in the vessel after stenting. It is demonstrated that stent designs and materials strongly affect the stenting-induced damage in the media layer and the subsequent development of in-stent restenosis. The larger the artery expansion achieved during balloon inflation, the higher the damage introduced to the media layer, leading to an increased level of in-stent restenosis. In addition, the development of in-stent restenosis is directly correlated with the artery expansion during the stent deployment. The correlation is further used to predict the effect of a complex clinical procedure, such as stent overlapping, on the level of in-stent restenosis developed after percutaneous coronary intervention.

Location-Specific Comparison Between a 3D In-Stent Restenosis Model and Micro-CT and Histology Data from Porcine In Vivo Experiments

BACKGROUND

Coronary 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.

METHODS

This 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.

RESULTS

Metrics 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.

CONCLUSIONS

The 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.

InSilc: 3D Reconstruction and plaque characterization tool

Coronary artery disease (CAD) is the leading cause of mortality in Europe and worldwide. Atherosclerosis is the most common pathologic process that is highly related with CAD, while the implantation of drug-eluting Bioresorbable Vascular Scaffolds (BVS) is the most promising procedure for treating patients with CAD. InSilc is an textbfin silico clinical trial (ISCT) platform for the development and assessment of drugeluting BVS. The InSilc platform provides insight in the performance of drug-eluting BVS in their short term and medium/long term through the Mechanical Modelling Module, the Deployment Module, the Fluid Dynamics Module, the Myocardial Perfusion Module, the Drug-delivery Module and the Degradation Module. In order for the aforementioned modules to be developed, the utilization of the reconstructed patient specific arterial segment and the BVS design are required, which is achieved through the 3D reconstruction and plaque characterization tool.In this study, the overall architecture of the InSilc platform is presented with special emphasis on the 3D reconstruction and plaque characterization tool. The tool will be able to implement different medical image processing workflows. The workflows will require minimum user intervention in order to be used in large scale clinical trials.

Mathematical modelling of the restenosis process after stent implantation

The stenting procedure has evolved to become a highly successful technique for the clinical treatment of advanced atherosclerotic lesions in arteries. However, the development of in-stent restenosis remains a key problem. In this work, a novel two-dimensional continuum mathematical model is proposed to describe the complex restenosis process following the insertion of a stent into a coronary artery. The biological species considered to play a key role in restenosis development are growth factors, matrix metalloproteinases, extracellular matrix, smooth muscle cells and endothelial cells. Diffusion-reaction equations are used for modelling the mass balance between species in the arterial wall. Experimental data from the literature have been used in order to estimate model parameters. Moreover, a sensitivity analysis has been performed to study the impact of varying the parameters of the model on the evolution of the biological species. The results demonstrate that this computational model qualitatively captures the key characteristics of the lesion growth and the healing process within an artery subjected to non-physiological mechanical forces. Our results suggest that the arterial wall response is driven by the damage area, smooth muscle cell proliferation and the collagen turnover among other factors.

Emerging regulators of vascular smooth muscle cell migration

Vascular smooth muscle cells (VSMCs) are the predominant cell type in the blood vessel wall and normally adopt a quiescent, contractile phenotype. VSMC migration is tightly controlled, however, disease associated changes in the soluble and insoluble environment promote VSMC migration. Classically, studies investigating VSMC migration have described the influence of soluble factors. Emerging data has highlighted the importance of insoluble factors, including extracellular matrix stiffness and porosity. In this review, we will recap on the important signalling pathways that regulate VSMC migration and reflect on the potential importance of emerging regulators of VSMC function.

A versatile hybrid agent-based, particle and partial differential equations method to analyze vascular adaptation

Peripheral arterial occlusive disease is a chronic pathology affecting at least 8–12 million people in the USA, typically treated with a vein graft bypass or through the deployment of a stent in order to restore the physiological circulation. Failure of peripheral endovascular interventions occurs at the intersection of vascular biology, biomechanics, and clinical decision making. It is our hypothesis that the majority of endovascular treatment approaches share the same driving mechanisms and that a deep understanding of the adaptation process is pivotal in order to improve the current outcome of the procedure. The postsurgical adaptation of vein graft bypasses offers the perfect example of how the balance between intimal hyperplasia and wall remodeling determines the failure or the success of the intervention. Accordingly, this work presents a versatile computational model able to capture the feedback loop that describes the interaction between events at cellular/tissue level and mechano-environmental conditions. The work here presented is a generalization and an improvement of a previous work by our group of investigators, where an agent-based model uses a cellular automata principle on a fixed hexagonal grid to reproduce the leading events of the graft’s restenosis. The new hybrid model here presented allows a more realistic simulation both of the biological laws that drive the cellular behavior and of the active role of the membranes that separate the various layers of the vein. The novel feature is to use an immersed boundary implementation of a highly viscous flow to represent SMC motility and matrix reorganization in response to graft adaptation. Our implementation is modular, and this makes us able to choose the right compromise between closeness to the physiological reality and complexity of the model. The focus of this paper is to offer a new modular implementation that combines the best features of an agent-based model, continuum mechanics, and particle-tracking methods to cope with the multiscale nature of the adaptation phenomena. This hybrid method allows us to quickly test various hypotheses with a particular attention to cellular motility, a process that we demonstrated should be driven by mechanical homeostasis in order to maintain the right balance between cells and extracellular matrix in order to reproduce a distribution similar to histological experimental data from vein grafts.

A two-phase model of early fibrous cap formation in atherosclerosis

Atherosclerotic plaque growth is characterised by chronic, non-resolving inflammation that promotes the accumulation of cellular debris and extracellular fat in the inner artery wall. This material is highly thrombogenic, and plaque rupture can lead to the formation of blood clots that occlude major arteries and cause myocardial infarction or stroke. In advanced plaques, vascular smooth muscle cells (SMCs) are recruited from deeper in the artery wall to synthesise a cap of fibrous tissue that stabilises the plaque and sequesters the thrombogenic plaque content from the bloodstream. The fibrous cap provides crucial protection against the clinical consequences of atherosclerosis, but the mechanisms of cap formation are poorly understood. In particular, it is unclear why certain plaques become stable and robust while others become fragile and dangerously vulnerable to rupture. We develop a multiphase model with non-standard boundary conditions to investigate early fibrous cap formation in the atherosclerotic plaque. The model is parameterised using data from a range of in vitro and in vivo studies, and includes highly nonlinear mechanisms of SMC proliferation and migration in response to an endothelium-derived chemical signal. We demonstrate that the model SMC population naturally evolves towards a steady-state, and predict a rate of cap formation and a final plaque SMC content consistent with experimental observations in mice. Parameter sensitivity simulations show that SMC proliferation makes a limited contribution to cap formation, and demonstrate that stable cap formation relies primarily on a critical balance between the rates of SMC recruitment to the plaque, chemotactic SMC migration within the plaque and SMC loss by apoptosis or phenotype change. This model represents the first detailed in silico study of fibrous cap formation in atherosclerosis, and establishes a multiphase modelling framework that can be readily extended to investigate many other aspects of plaque development.

A cell-based mechanical model of coronary artery tunica media

A three-dimensional cell-based mechanical model of coronary artery tunica media is proposed. The model is composed of spherical cells forming a hexagonal close-packed lattice. Tissue anisotropy is taken into account by varying interaction forces with the direction of intercellular connection. Several cell-centre interaction potentials for repulsion and attraction are considered, including the Hertz contact model and its neo-Hookean extension, the Johnson-Kendall-Roberts model of adhesive contact, and a wormlike chain model. The model is validated against data from in vitro uni-axial tension tests performed on dissected strips of tunica media. The wormlike chain potential in combination with the neo-Hookean Hertz contact model produces stress-stretch curves which represent the experimental data very well.

A Comparison of Fully-Coupled 3D In-Stent Restenosis Simulations to In-vivo Data

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.

Prelamin A Accumulation Attenuates Rac1 Activity and Increases the Intrinsic Migrational Persistence of Aged Vascular Smooth Muscle Cells

Vascular smooth muscle cell (VSMC) motility is essential during both physiological and pathological vessel remodeling. Although ageing has emerged as a major risk factor in the development of cardiovascular disease, our understanding of the impact of ageing on VSMC motility remains limited. Prelamin A accumulation is known to drive VSMC ageing and we show that presenescent VSMCs, that have accumulated prelamin A, display increased focal adhesion dynamics, augmented migrational velocity/persistence and attenuated Rac1 activity. Importantly, prelamin A accumulation in proliferative VSMCs, induced by depletion of the prelamin A processing enzyme FACE1, recapitulated the focal adhesion, migrational persistence and Rac1 phenotypes observed in presenescent VSMCs. Moreover, lamin A/C-depleted VSMCs also display reduced Rac1 activity, suggesting that prelamin A influences Rac1 activity by interfering with lamin A/C function at the nuclear envelope. Taken together, these data demonstrate that lamin A/C maintains Rac1 activity in VSMCs and prelamin A disrupts lamin A/C function to reduce Rac1 activity and induce migrational persistence during VSMC ageing.

Towards the virtual artery: a multiscale model for vascular physiology at the physics-chemistry-biology interface

This discussion paper introduces the concept of the Virtual Artery as a multiscale model for arterial physiology and pathologies at the physics-chemistry-biology (PCB) interface. The cellular level is identified as the mesoscopic level, and we argue that by coupling cell-based models with other relevant models on the macro- and microscale, a versatile model of arterial health and disease can be composed. We review the necessary ingredients, both models of arteries at many different scales, as well as generic methods to compose multiscale models. Next, we discuss how this can be combined into the virtual artery. Finally, we argue that the concept of models at the PCB interface could or perhaps should become a powerful paradigm, not only as in our case for studying physiology, but also for many other systems that have such PCB interfaces.This article is part of the themed issue 'Multiscale modelling at the physics-chemistry-biology interface'.