Predicting neointimal hyperplasia in stented arteries using time-dependant computational fluid dynamics: A review

In-stent Restenosis After Stenting for Superior Mesenteric Artery Dissection Is Associated With Stent Landing Zone: From Clinical Prediction to Hemodynamic Mechanisms

OBJECTIVE

To identify risk factors for in-stent restenosis (ISR) in patients undergoing stent placement for superior mesenteric artery dissection (SMAD) and to determine the hemodynamic mechanism underlying ISR.

METHODS

For this retrospective study, patients with SMAD who had ISR after stent placement were included in the ISR group, and age- and sex-matched patients with SMAD who did not experience ISR after stent placement were included in the control group. Clinical, imaging, and hemodynamic data were assessed. Multivariable regression was used to identify independent ISR risk factors. Structural and fluid dynamics simulations were applied to determine the hemodynamic mechanism underlying the occurrence of ISR.

RESULTS

The study population included 26 patients with ISR and 26 control patients. Multivariate analysis demonstrated that stent-to-vascular (S/V) ratio (odds ratio [OR], 1.14; 95% confidence interval [CI]: 1.00-1.29; p=0.045), stent proximal position >10 mm away from the SMA root (OR, 108.67; 95% CI: 3.09-3816.42; p=0.010), and high oscillatory shear index (OSI) area (OR, 1.25; 95% CI: 1.02-1.52; p=0.029) were predictors of ISR. In structural and fluid dynamics simulations, a stent proximal position near the abdominal aorta (AA) or entering into the AA reduced the contact area between the proximal struts of the stent and the vascular wall, and alleviated the distal lumen overdilation.

CONCLUSION

The S/V ratio, stent proximal position away from the SMA root (>10 mm), and high OSI area are independent risk factors for ISR in patients with SMAD undergoing stent placement. Deploying the proximal end of the stent near the AA or entering into the AA appears to improve the hemodynamic environment in the SMA lumen and ultimately reduce the risk of ISR.

CLINICAL IMPACT

In-stent restenosis is an uncommon but potentially catastrophic complication after stent placement for the management of superior mesenteric artery dissection. This study identified risk factors for in-stent restenosis and demonstrated that, as long as the stent can fully cover the dissection range, deploying the proximal end of the stent near the abdominal aorta or less entering into the abdominal aorta may reduce the risk of in-stent restenosis in this patient population.

Hemodynamic and Geometric Risk Factors for In-Stent Restenosis in Patients with Intracranial Atherosclerotic Stenosis

Methods

Severe ICAS patients managed with percutaneous transluminal angioplasty and stenting (PTAS) were included in the retrospective cohort study and were divided into two groups according to whether ISR occurred at follow-up (ISR group and no-ISR group). Computational fluid dynamics models were built based on digital subtraction angiography before and after PTAS to simulate blood flow and quantify hemodynamic parameters. The associations between vessel geometry, hemodynamics, and ISR in ICAS patients were investigated.

Results

Among 39 patients, ISR occurred in seven patients (17.95%) after a mean follow-up period of 6.69 ± 3.24 months. Stenting decreased vessel angulation (51.11° [40.07°-67.27°] vs. 15.97° [0.00°-36.16°], P = 0.000) and vessel tortuosity (0.09 [0.06-0.13] vs. 0.01 [0.00-0.03], P = 0.000). Meanwhile, the translational pressure ratio (PR) dramatically increased (0.07 [0.00-0.31] vs. 0.62 [0.41-0.82], P = 0.000) with the wall shear stress ratio decreased (13.93 [8.37-40.30] vs. 2.90 [1.69-4.48], P = 0.000). In the multivariate analysis, smaller Δ tortuosity (P = 0.038) was independently associated with the occurrence of ISR, and smaller post-PTAS translesional PR was also a predictive factor of marginal significance (P = 0.059).

Conclusion

PTAS decreased vessel angulation, vessel tortuosity, and translesional wall shear stress ratio while it increased translesional pressure ratio (PR) dramatically in ICAS patients. Smaller Δ tortuosity was found to be a risk factor for ISR, and smaller post-PTAS translesional PR was also a predictive factor of marginal significance, indicating that both geometric and hemodynamic parameters played important roles in the occurrence of ISR after PTAS.

Hemodynamic alternations following stent deployment and post-dilation in a heavily calcified coronary artery: In silico and ex-vivo approaches

In this work, hemodynamic alterations in a patient-specific, heavily calcified coronary artery following stent deployment and post-dilations are quantified using in silico and ex-vivo approaches. Three-dimensional artery models were reconstructed from OCT images. Stent deployment and post-dilation with various inflation pressures were performed through both the finite element method (FEM) and ex vivo experiments. Results from FEM agreed very well with the ex-vivo measurements, interms of lumen areas, stent underexpansion, and strut malapposition. In addition, computational fluid dynamics (CFD) simulations were performed to delineate the hemodynamic alterations after stent deployment and post-dilations. A pressure time history at the inlet and a lumped parameter model (LPM) at the outlet were adopted to mimic the aortic pressure and the distal arterial tree, respectively. The pressure drop across the lesion, pertaining to the clinical measure of instantaneous wave-free flow ratio (iFR), was investigated. Results have shown that post-dilations are necessary for the lumen gain as well as the hemodynamic restoration towards hemostasis. Malapposed struts induced much higher shear rate, flow disturbances and lower time-averaged wall shear stress (TAWSS) around struts. Post-dilations mitigated the strut malapposition, and thus the shear rate. Moreover, stenting induced larger area of low TAWSS (<0.4 Pa) and lager volume of high shear rate (>2000 s^-1), indicating higher risks of in-stent restenosis (ISR) and stent thrombosis (ST), respectively. Oscillatory shear index (OSI) and relative residence time (RRT) indicated the wall regions more prone to ISR are located near the malapposed stent struts.

Impact of Malapposed and Overlapping Stents on Hemodynamics: A 2D Parametric Computational Fluid Dynamics Study

Despite significant progress, malapposed or overlapped stents are a complication that affects daily percutaneous coronary intervention (PCI) procedures. These malapposed stents affect blood flow and create a micro re-circulatory environment. These disturbances are often associated with a change in Wall Shear Stress (WSS), Time-averaged WSS (TAWSS), relative residence time (RRT) and oscillatory character of WSS and disrupt the delicate balance of vascular biology, providing a possible source of thrombosis and restenosis. In this study, 2D axisymmetric parametric computational fluid dynamics (CFD) simulations were performed to systematically analyze the hemodynamic effects of malapposition and stent overlap for two types of stents (drug-eluting stent and a bioresorbable stent). The results of the modeling are mainly analyzed using streamlines, TAWSS, oscillatory shear index (OSI) and RRT. The risks of restenosis and thrombus are evaluated according to commonly accepted thresholds for TAWSS and OSI. The small malapposition distances (MD) cause both low TAWSS and high OSI, which are potential adverse outcomes. The region of low OSI decrease with MD. Overlap configurations produce areas with low WSS and high OSI. The affected lengths are relatively insensitive to the overlap distance. The effects of strut size are even more sensitive and adverse for overlap configurations compared to a well-applied stent.

Influence of Artery Straightening on Local Hemodynamics in Left Anterior Descending (LAD) Artery after Stent Implantation

Objectives

The study investigates local hemodynamic environment changes caused by straightening phenomenon and the relationship between straightening phenomenon and in-stent restenosis.

Background

Intravascular intervention is an effective treatment in restoring the normal flow conditions and vascular lumen. Unfortunately, in-stent restenosis often occurs in a subset of patients after stent implantation and limits the success of stent implantation outcomes. The implanted stent usually causes artery straightening locally, rather than coinciding and adjusting to the physiological curve exactly. Artery straightening would apparently modify the artery geometry and therefore alter the local hemodynamic environment, which may result in intimal hyperplasia and restenosis after stenting implantation.

Methods

In the current investigation, we verify the hypothesis that the artery straightening influences the local hemodynamic state using the different 3D CT models. Flow analysis for blood in the left anterior descending coronary artery and the straightening model is simulated numerically.

Result

The current results reveal that the straightening phenomenon alters the distribution of wall shear stress and flow patterns, decreases the wall shear stress (WSS), and increases the oscillatory shear index (OSI) and the relative residence time (RRT), especially at the proximal and distal areas of stenting.

Conclusions

The local straightened geometry established after stent implantation was likely to generate portions of the stenting area to a high risk of neointimal hyperplasia and subsequent restenosis.

The Hemodynamic Effect of Enhanced External Counterpulsation Treatment on Atherosclerotic Plaque in the Carotid Artery: A Framework of Patient-Specific Computational Fluid Dynamics Analysis

Long-term enhanced external counterpulsation (EECP) therapy has been recommended for antiatherogenesis in recent clinical observations and trials. However, the precise mechanism underlying the benefits has not been fully clarified. To quantify the effect of EECP intervention on arterial hemodynamic environment, a framework of numerical assessment was introduced using a parallel computing algorithm. A 3D endothelial surface of the carotid artery with mild atherosclerotic plaque was constructed from images of magnetic resonance angiography (MRA). Physiologic boundary conditions were derived from images of the ultrasound flow velocity spectrum measured at the common carotid artery and before and during EECP intervention. Hemodynamic factors relating to wall shear stress (WSS) and its spatial and temporal fluctuations were calculated and analyzed, which included AWSS, OSI, and AWSSG. Measuring and computational results showed that diastole blood pressure, perfusion, and WSS level in carotid bifurcation were significantly increased during EECP intervention. Mean AWSS level throughout the model increased by 16.9%, while OSI level did not show a significant change during EECP. We thus suggested that long-term EECP treatment might inhibit the initiation and development of atherosclerotic plaque via improving the hemodynamic environment in the carotid artery. Meanwhile, EECP performance induced a 19.6% increase in AWSSG level, and whether it would influence the endothelial functions may need a further study. Moreover, the numerical method proposed in this study was expected to be useful for the instant assessment of clinical application of EECP .

Computational modelling of magnesium stent mechanical performance in a remodelling artery: Effects of multiple remodelling stimuli

Significant research has been conducted in the area of coronary stents/scaffolds made from resorbable metallic and polymeric biomaterials. These next-generation bioabsorbable stents have the potential to completely revolutionise the treatment of coronary artery disease. The primary advantage of resorbable devices over permanent stents is their temporary presence which, from a theoretical point of view, means only a healed coronary artery will be left behind following degradation of the stent potentially eliminating long-term clinical problems associated with permanent stents. The healing of the artery following coronary stent/scaffold implantation is crucial for the long-term safety of these devices. Computational modelling can be used to evaluate the performance of complex stent devices in silico and assist in the design and development and understanding of the next-generation resorbable stents. What is lacking in computational modelling literature is the representation of the active response of the arterial tissue in the weeks and months following stent implantation, ie, neointimal remodelling, in particular for the case of biodegradable stents. In this paper, a computational modelling framework is developed, which accounts for two major physiological stimuli responsible for neointimal remodelling and combined with a magnesium corrosion model that is capable of simulating localised pitting (realistic) stent corrosion. The framework is used to simulate different neointimal growth patterns and to explore the effects the neointimal remodelling has on the mechanical performance (scaffolding support) of the bioabsorbable magnesium stent.

An in silico study of the influence of vessel wall deformation on neointimal hyperplasia progression in peripheral bypass grafts

Neointimal hyperplasia (NIH) is a major obstacle to graft patency in the peripheral arteries. A complex interaction of biomechanical factors contribute to NIH development and progression, and although haemodynamic markers such as wall shear stress have been linked to the disease, these have so far been insufficient to fully capture its behaviour. Using a computational model linking computational fluid dynamics (CFD) simulations of blood flow with a biochemical model representing NIH growth mechanisms, we analyse the effect of compliance mismatch, due to the presence of surgical stitches and/or to the change in distensibility between artery and vein graft, on the haemodynamics in the lumen and, subsequently, on NIH progression. The model enabled to simulate NIH at proximal and distal anastomoses of three patient-specific end-to-side saphenous vein grafts under two compliance-mismatch configurations, and a rigid wall case for comparison, obtaining values of stenosis similar to those observed in the computed tomography (CT) scans. The maximum difference in time-averaged wall shear stress between the rigid and compliant models was 3.4 Pa, and differences in estimation of NIH progression were only observed in one patient. The impact of compliance on the haemodynamic-driven development of NIH was small in the patient-specific cases considered.

Early strut protrusion and late neointima thickness in the Absorb bioresorbable scaffold: a serial wall shear stress analysis up to five years

AIMS

The aim of the study was to evaluate the effect of strut protrusion (SP) on wall shear stress (WSS) and neointimal growth (NG) one and five years after implantation of an Absorb bioresorbable vascular scaffold.

METHODS AND RESULTS

Eight patients were selected from a first-in-man study. Following three-dimensional (3D) reconstruction of coronaries, WSS was quantified using Newtonian steady-flow simulation in each cross-section at 5° subunits (sectors) of the circumferential luminal surface. At one year, neointimal thickness (NT) was measured by optical coherence tomography (OCT) and correlated to WSS and SP post procedure. Median SP was 112.9 (90.8, 133.1) µm post implantation. Post procedure, a logarithmic inverse relationship between SP and post-implantation WSS (r=-0.425, p<0.001; correlation coefficients in a range from -0.143 to -0.553) was observed, whereas a correlation between baseline logarithm-transformed WSS (log-WSS) and NT (r=-0.451, p<0.001; correlation coefficients ranged from -0.140 to -0.662) was documented at one year. Mixed-effects analysis between baseline log-WSS and NT at follow-up yielded a slope of 30 µm/ln Pascal (Pa) and a y-intercept of 98 µm. As a result of NG, median flow area decreased from 6.91 (6.53, 7.48) mm2 post implantation to 5.65 (5.47, 6.02) mm2 at one-year follow-up (p=0.01) and to 5.75±1.37 mm2 at five-year follow-up (p=0.024). However, the vessel surface exposed to low WSS (<1 Pa) decreased significantly post procedure (42%) to one year (35.9%) and five years (15.2%) (p-overall <0.0001).

CONCLUSIONS

SP disturbs laminar flow, creates regions of low WSS (<1.0 Pa) that are associated with NG and lumen area reduction. Low WSS post implantation reduced significantly at long-term follow-up. Thin struts with effective embedment would substantially reduce NG and accelerate homogenisation of WSS towards physiological values.

Influence of microvascular sutures on shear strain rate in realistic pulsatile flow

Arterial thrombus formation is directly related to the mechanical shear experienced by platelets within flow. High shear strain rates (SSRs) and large shear gradients cause platelet activation, aggregation and production of thrombus. This study, for the first time, investigates the influence of pulsatile flow on local haemodynamics within sutured microarterial anastomoses. We measured physiological arterial waveform velocities experimentally using Doppler ultrasound velocimetry, and a representative example was applied to a realistic sutured microarterial geometry. Computational geometries were created using measurements taken from sutured chicken femoral arteries. Arterial SSRs were predicted using computational fluid dynamics (CFD) software, to indicate the potential for platelet activation, deposition and thrombus formation. Predictions of steady and sinusoidal inputs were compared to analyse whether the addition of physiological pulse characteristics affects local intravascular flow characteristics. Simulations were designed to evaluate flow in pristine and hand-sutured microarterial anastomoses, each with a steady-state and sinusoidal pulse component. The presence of sutures increased SSR_max in the anastomotic region by factors of 2.1 and 2.3 in steady-state and pulsatile flows respectively, when compared to a pristine vessel. SSR values seen in these simulations are analogous to the presence of moderate arterial stenosis. Steady-state simulations, driven by a constant inflow velocity equal to the peak systolic velocity (PSV) of the measured pulsatile flow, underestimated SSRs by ∼ 9% in pristine, and ∼ 19% in sutured vessels compared with a realistic pulse. Sinusoidal flows, with equivalent frequency and amplitude to a measured arterial waveform, represent a slight improvement on steady-state simulations, but still SSRs are underestimated by 1-2%. We recommend using a measured arterial waveform, of the form presented here, for simulating pulsatile flows in vessels of this nature. Under realistic pulsatile flow, shear gradients across microvascular sutures are high, of the order ∼ 7.9 × 10⁶ m^-1 s^-1, which may also be associated with activation of platelets and formation of aggregates.

Determining the cross-talk between smooth muscle cells and macrophages on a cobalt-chromium stent material surface using an in vitro postimplantation coculture model

Smooth muscle cells (SMCs) and macrophages are important cellular components involved in the development of complications following the implantation of cardiovascular devices. This leads to various disorders such as restenosis, chronic inflammation, and may ultimately result in device failure. In this study, we developed a postimplant stent coculture model using different ratios of SMCs and macrophages seeded on to cobalt-chromium alloy. The macrophages had an increased affinity to the coculture surfaces, which resulted in decreased SMC attachment to the alloy surfaces at the initial time point. Once adhered, the macrophages spread freely and displayed advanced stages of inflammation at 48 h when cocultured with SMCs. This resulted in an increased secretion of proinflammatory cytokines (tumor necrosis factor alpha, monocyte chemotactic protein 1, interleukin [IL]-6, and IL-8) by 48 h in the coculture samples with the greatest increase observed with the high number of macrophages. Therefore, the increased levels of proinflammatory cytokines promoted the growth of SMCs in coculture to a greater extent than when the SMCs were culture alone. Thus, this study demonstrated the constant cross-talk between SMCs and macrophages occurring on the postimplant stent surface. Similar coculture models can be used to test the biocompatibility of drugs and biomaterials at possible postimplantation scenarios. © 2017 Wiley Periodicals, Inc. J Biomed Mater Res Part A: 106A: 673-685, 2018.

Patient-Specific, Multi-Scale Modeling of Neointimal Hyperplasia in Vein Grafts

Neointimal hyperplasia is amongst the major causes of failure of bypass grafts. The disease progression varies from patient to patient due to a range of different factors. In this paper, a mathematical model will be used to understand neointimal hyperplasia in individual patients, combining information from biological experiments and patient-specific data to analyze some aspects of the disease, particularly with regard to mechanical stimuli due to shear stresses on the vessel wall. By combining a biochemical model of cell growth and a patient-specific computational fluid dynamics analysis of blood flow in the lumen, remodeling of the blood vessel is studied by means of a novel computational framework. The framework was used to analyze two vein graft bypasses from one patient: a femoro-popliteal and a femoro-distal bypass. The remodeling of the vessel wall and analysis of the flow for each case was then compared to clinical data and discussed as a potential tool for a better understanding of the disease. Simulation results from this first computational approach showed an overall agreement on the locations of hyperplasia in these patients and demonstrated the potential of using new integrative modeling tools to understand disease progression.

Computational fluid dynamics study of common stent models inside idealised curved coronary arteries

The haemodynamic behaviour of blood inside a coronary artery after stenting is greatly affected by individual stent features as well as complex geometrical properties of the artery including tortuosity and curvature. Regions at higher risk of restenosis, as measured by low wall shear stress (WSS < 0.5 Pa), have not yet been studied in detail in curved stented arteries. In this study, three-dimensional computational modelling and computational fluid dynamics methodologies were used to analyse the haemodynamic characteristics in curved stented arteries using several common stent models. Results in this study showed that stent strut thickness was one major factor influencing the distribution of WSS in curved arteries. Regions of low WSS were found behind struts, particularly those oriented at a large angle relative to the streamwise flow direction. These findings were similar to those obtained in studies of straight arteries. An uneven distribution of WSS at the inner and outer bends of curved arteries was observed where the WSS was lower at the inner bend. In this study, it was also shown that stents with a helical configuration generated an extra swirling component of the flow based on the helical direction; however, this extra swirl in the flow field did not cause significant changes on the distribution of WSS under the current setup.

Local blood flow patterns in stented coronary bifurcations: an experimental and numerical study

PURPOSE

Despite the atheroprone environment of blood flow in coronary bifurcations, limited quantitative information is available on the hemodynamics occurring in these geometries, both before and after their treatment with endovascular stents. Previous studies have focused on computational fluid dynamics (CFD) analyses and have bypassed the challenges associated with experimentally representing the flow environment, providing no means for validation. This study analyzed steady flow conditions in 3 bifurcation angles and 4 different single- and double-stenting procedures, which are used clinically in coronary bifurcations.

METHODS

The numerical aspect of this study utilized geometries derived from CAD models (nonstented cases) and finite element simulations (stented cases). Digital particle image velocimetry (DPIV) testing was conducted within compliant bifurcating models for which an uncertainty analysis was performed at each measurement location for CFD validation purposes. Results were analyzed in terms of velocity magnitude contour maps and axial velocity profiles at several locations in the bifurcated vessels.

RESULTS AND CONCLUSIONS

Qualitatively, the 2 approaches showed agreement in the bulk flow patterns. However, the velocity computed with CFD was outside the DPIV uncertainty estimates, which can be attributed to the intrinsic difference and modeling assumptions of the 2 approaches. The findings reveal that wider bifurcation angles and double-stenting procedures are both characterized by increased areas of low flow and recirculation. Additionally, inferior performance in terms of viscous and wall shear stresses was observed in double-stented cases.

Measurement and modeling of coronary blood flow

Ischemic heart disease that comprises both coronary artery disease and microvascular disease is the single greatest cause of death globally. In this context, enhancing our understanding of the interaction of coronary structure and function is not only fundamental for advancing basic physiology but also crucial for identifying new targets for treating these diseases. A central challenge for understanding coronary blood flow is that coronary structure and function exhibit different behaviors across a range of spatial and temporal scales. While experimental studies have sought to understand this feature by isolating specific mechanisms, in tandem, computational modeling is increasingly also providing a unique framework to integrate mechanistic behaviors across different scales. In addition, clinical methods for assessing coronary disease severity are continuously being informed and updated by findings in basic physiology. Coupling these technologies, computational modeling of the coronary circulation is emerging as a bridge between the experimental and clinical domains, providing a framework to integrate imaging and measurements from multiple sources with mathematical descriptions of governing physical laws. State-of-the-art computational modeling is being used to combine mechanistic models with data to provide new insight into coronary physiology, optimization of medical technologies, and new applications to guide clinical practice.

Numerical investigations of the haemodynamic changes associated with stent malapposition in an idealised coronary artery

The deployment of a coronary stent near complex lesions can sometimes lead to incomplete stent apposition (ISA), an undesirable side effect of coronary stent implantation. Three-dimensional computational fluid dynamics (CFD) calculations are performed on simplified stent models (with either square or circular cross-section struts) inside an idealised coronary artery to analyse the effect of different levels of ISA to the change in haemodynamics inside the artery. The clinical significance of ISA is reported using haemodynamic metrics like wall shear stress (WSS) and wall shear stress gradient (WSSG). A coronary stent with square cross-sectional strut shows different levels of reverse flow for malapposition distance (MD) between 0mm and 0.12 mm. Chaotic blood flow is usually observed at late diastole and early systole for MD=0mm and 0.12 mm but are suppressed for MD=0.06 mm. The struts with circular cross section delay the flow chaotic process as compared to square cross-sectional struts at the same MD and also reduce the level of fluctuations found in the flow field. However, further increase in MD can lead to chaotic flow not only at late diastole and early systole, but it also leads to chaotic flow at the end of systole. In all cases, WSS increases above the threshold value (0.5 Pa) as MD increases due to the diminishing reverse flow near the artery wall. Increasing MD also results in an elevated WSSG as flow becomes more chaotic, except for square struts at MD=0.06 mm.

From Histology and Imaging Data to Models for In-Stent Restenosis

The implantation of stents has been used to treat coronary artery stenosis for several decades. Although stenting is successful in restoring the vessel lumen and is a minimally invasive approach, the long-term outcomes are often compromised by in-stent restenosis (ISR). Animal models have provided insights into the pathophysiology of ISR and are widely used to evaluate candidate drug inhibitors of ISR. Such biological models allow the response of the vessel to stent implantation to be studied without the variation of lesion characteristics encountered in patient studies.

This paper describes the development of complementary in silico models employed to improve the understanding of the biological response to stenting using a porcine model of restenosis. This includes experimental quantification using microCT imaging and histology and the use of this data to establish numerical models of restenosis. Comparison of in silico results with histology is used to examine the relationship between spatial localization of fluid and solid mechanics stimuli immediately post-stenting. Multi-scale simulation methods are employed to study the evolution of neointimal growth over time and the variation in the extent of neointimal hyperplasia within the stented region. Interpretation of model results through direct comparison with the biological response contributes to more detailed understanding of the pathophysiology of ISR, and suggests the focus for follow-up studies.

In conclusion we outline the challenges which remain to both complete our understanding of the mechanisms responsible for restenosis and translate these models to applications in stent design and treatment planning at both population-based and patient-specific levels.

Simulation of oxygen transfer in stented arteries and correlation with in-stent restenosis

Computational models are used to study the combined effect of biomechanical and biochemical factors on coronary in-stent restenosis, which is a postoperative remodeling and regrowth pathology of the stented arteries. More precisely, we address numerical simulations, on the basis of Navier-Stokes and mass transport equations, to study the role of perturbed wall shear stresses and reduced oxygen concentration in a geometrical model reconstructed from a real porcine artery treated with a stent. Joining in vivo and in silico tools of investigation has multiple benefits in this case. On one hand, the geometry of the arterial wall and of the stent closely correspond to a real implanted configuration. On the other hand, the inspection of histological tissue samples informs us on the location and intensity of in-stent restenosis. As a result, we are able to correlate geometrical factors, such as the axial variation of the artery diameter and its curvature; the numerical quantification of biochemical stimuli, such as wall shear stresses; and the availability of oxygen to the inner layers of the artery, with the appearance of in-stent restenosis. This study shows that the perturbation of the vessel curvature could induce hemodynamic conditions that stimulate undesired arterial remodeling.

Construction and expression of retroviral vector pLEGFP-N1-TERT in preparation of seed cells for skin tissue engineering

OBJECTIVE

To construct the retroviral vector pLEGFP-N1-telomerase reverse transcriptase (TERT) and to investigate the expression of TERT in neonatal mouse hypodermal cells.

METHODS

The polymerase chain reaction (PCR)-amplified TERT gene was inserted into plasmid pLEGFP-N1. The positive clone was identified by restriction enzyme digestion and sequencing, then was transfected into packaging cells to produce retrovirus particles. Neonatal mouse hypodermal cells were infected with the virus to generate a stable cell line. The TERT mRNA expression level, telomerase activity, and enhanced green fluorescent protein (EGFP) expression level were analyzed.

RESULTS

Retroviral vector pLEGFP-N1-TERT was constructed successfully, and a stable cell line of neonatal mouse hypodermal cells expressing EGFP was established. Western blot and immunohistochemical assay showed that the expression level of TERT was significantly elevated in the neonatal mouse hypodermal cells.

CONCLUSIONS

A high titer of retrovirus pLEGFP-N1-TERT mediates high-level expression of the exogenous TERT gene in the neonatal mouse hypodermal cells. This protocol has potential applications for skin tissue engineering and cell transplantation therapy.

Computational fluid dynamic simulations of image-based stented coronary bifurcation models

One of the relevant phenomenon associated with in-stent restenosis in coronary arteries is an altered haemodynamics in the stented region. Computational fluid dynamics (CFD) offers the possibility to investigate the haemodynamics at a level of detail not always accessible within experimental techniques. CFD can quantify and correlate the local haemodynamics structures which might lead to in-stent restenosis. The aim of this work is to study the fluid dynamics of realistic stented coronary artery models which replicate the complete clinical procedure of stent implantation. Two cases of pathologic left anterior descending coronary arteries with their bifurcations are reconstructed from computed tomography angiography and conventional coronary angiography images. Results of wall shear stress and relative residence time show that the wall regions more prone to the risk of restenosis are located next to stent struts, to the bifurcations and to the stent overlapping zone for both investigated cases. Considering a bulk flow analysis, helical flow structures are generated by the curvature of the zone upstream from the stent and by the bifurcation regions. Helical recirculating microstructures are also visible downstream from the stent struts. This study demonstrates the feasibility to virtually investigate the haemodynamics of patient-specific coronary bifurcation geometries.

Compound ex vivo and in silico method for hemodynamic analysis of stented arteries

Hemodynamic factors such as low wall shear stress have been shown to influence endothelial healing and atherogenesis in stent-free vessels. However, in stented vessels, a reliable quantitative analysis of such relations has not been possible due to the lack of a suitable method for the accurate acquisition of blood flow. The objective of this work was to develop a method for the precise reconstruction of hemodynamics and quantification of wall shear stress in stented vessels. We have developed such a method that can be applied to vessels stented in or ex vivo and processed ex vivo. Here we stented the coronary arteries of ex vivo porcine hearts, performed vascular corrosion casting, acquired the vessel geometry using micro-computed tomography and reconstructed blood flow and shear stress using computational fluid dynamics. The method yields accurate local flow information through anatomic fidelity, capturing in detail the stent geometry, arterial tissue prolapse, radial and axial arterial deformation as well as strut malapposition. This novel compound method may serve as a unique tool for spatially resolved analysis of the relationship between hemodynamic factors and vascular biology. It can further be employed to optimize stent design and stenting strategies.

Patient-specific simulations of stenting procedures in coronary bifurcations: two clinical cases

Computational simulations of stenting procedures in idealized geometries can only provide general guidelines and their use in the patient-specific planning of percutaneous treatments is inadequate. Conversely, image-based patient-specific tools that are able to realistically simulate different interventional options might facilitate clinical decision-making and provide useful insights on the treatment for each individual patient. The aim of this work is the implementation of a patient-specific model that uses image-based reconstructions of coronary bifurcations and is able to replicate real stenting procedures following clinical indications. Two clinical cases are investigated focusing the attention on the open problems of coronary bifurcations and their main treatment, the provisional side branch approach. Image-based reconstructions are created combining the information from conventional coronary angiography and computed tomography angiography while structural finite element models are implemented to replicate the real procedure performed in the patients. First, numerical results show the biomechanical influence of stents deployment in the coronary bifurcations during and after the procedures. In particular, the straightening of the arterial wall and the influence of two overlapping stents on stress fields are investigated here. Results show that a sensible decrease of the vessel tortuosity occurs after stent implantation and that overlapping devices result in an increased stress state of both the artery and the stents. Lastly, the comparison between numerical and image-based post-stenting configurations proved the reliability of such models while replicating stent deployment in coronary arteries.

Modeling stented coronary arteries: where we are, where to go

In the last two decades, numerical models have become well-recognized and widely adopted tools to investigate stenting procedures. Due to limited computational resources and modeling capabilities, early numerical studies only involved simplified cases and idealized stented arteries. Nowadays, increased computational power allows for numerical models to meet clinical needs and include more complex cases such as the implantation of multiple stents in bifurcations or curved vessels. Interesting progresses have been made in the numerical modeling of stenting procedures both from a structural and a fluid dynamics points of view. Moreover, in the drug eluting stents era, new insights on drug elution capabilities are becoming essential in the stent development. Lastly, image-based methods able to reconstruct realistic geometries from medical images have been proposed in the recent literature aiming to better describe the peculiar anatomical features of coronary vessels and increase the accuracy of the numerical models. In this light, this review provides a comprehensive analysis of the current state-of-the-art in this research area, discussing the main methodological advances and remarkable results drawn from a number of significant studies.

Sequential Structural and Fluid Dynamic Numerical Simulations of a Stented Bifurcated Coronary Artery

Despite their success, stenting procedures are still associated to some clinical problems like sub-acute thrombosis and in-stent restenosis. Several clinical studies associate these phenomena to a combination of both structural and hemodynamic alterations caused by stent implantation. Recently, numerical models have been widely used in the literature to investigate stenting procedures but always from either a purely structural or fluid dynamic point of view. The aim of this work is the implementation of sequential structural and fluid dynamic numerical models to provide a better understanding of stenting procedures in coronary bifurcations. In particular, the realistic geometrical configurations obtained with structural simulations were used to create the fluid domains employed within transient fluid dynamic analyses. This sequential approach was applied to investigate the final kissing balloon (FKB) inflation during the provisional side branch technique. Mechanical stresses in the arterial wall and the stent as well as wall shear stresses along the arterial wall were examined before and after the FKB deployment. FKB provoked average mechanical stresses in the arterial wall almost 2.5 times higher with respect to those induced by inflation of the stent in the main branch only. Results also enlightened FKB benefits in terms of improved local blood flow pattern for the side branch access. As a drawback, the FKB generates a larger region of low wall shear stress. In particular, after FKB the percentage of area characterized by wall shear stresses lower than 0.5 Pa was 79.0%, while before the FKB it was 62.3%. For these reasons, a new tapered balloon dedicated to bifurcations was proposed. The inclusion of the modified balloon has reduced the mechanical stresses in the proximal arterial vessel to 40% and the low wall shear stress coverage area to 71.3%. In conclusion, these results show the relevance of the adopted sequential approach to study the wall mechanics and the hemodynamics created by stent deployment.

Computational flow dynamics of the severe m1 stenosis before and after stenting

Purpose

Computational flow dynamic (CFD) study has not been widely applied in intracranial artery stenosis due to requirement of high resolution in identifying the small intracranial artery. We described a process in CFD study applied to symptomatic severe intracranial (M1) stenosis before and after stenting.

Materials and Methods

Reconstructed 3D angiography in STL format was transferred to Magics (Materialise NV, Leuven, Belgium) for smoothing of vessel surface and trimming of branch vessels and to HyperMesh (Altair Engineering Inc., Auckland, New Zealand) for generating tetra volume mesh from triangular surface-meshed 3D angiogram. Computational analysis of blood flow in the blood vessels was performed using the commercial finite element software ADINA Ver 8.5 (ADINA R & D, Inc., Lebanon, MA). The distribution of wall shear stress (WSS), peak velocity and pressure in a patient was analyzed before and after intracranial stenting.

Results

Computer simulation of wall shear stress, flow velocity and wall pressure before and after stenting could be demonstrated three dimensionally by video mode according to flow vs. time dimension. Such flow model was well correlated with angiographic finding related to maximum degree of stenosis. Change of WSS, peak velocity and pressure at the severe stenosis was demonstrated before and after stenting. There was no WSS after stenting in case without residual stenosis.

Conclusion

Our study revealed that CFD analysis before and after intracranial stenting was feasible despite of limited vessel wall dimension and could reveal change of WSS as well as flow velocity and wall pressure.