The Atheroprotective Nature of Helical Flow in Coronary Arteries

Investigating the relationship between geometry and hemodynamics in an experimentally derived murine coronary computational model

Despite the critical role of coronary morphology and hemodynamics in the development of coronary artery disease (CAD), comprehensive analyses of these factors in murine models are limited. Our study integrates in vivo approaches with computational methods to yield a complete set of precise and reliable morphologic and hemodynamic measurements and to investigate their interrelationship in the left coronary artery of healthy C57BL/6 mice. The work utilizes advanced micro-computed tomography imaging, enhanced with Microfil® coronary perfusion, complemented by morphometric analysis and computational fluid dynamic simulation. Our results in murine coronary arteries show: i) bifurcations are the most geometrically complex regions, susceptible to disturbed hemodynamics and, consequently, endothelial dysfunction; ii) vascular endothelial cells experience wall shear stress (WSS) an order of magnitude greater than in humans, primarily due to their smaller size, although minimal WSS multi-directionality is noted in both species; iii) intravascular flow exhibits reduced helical patterns compared to human coronaries, indicating a need for further investigation into their potential protective role against disease onset; and iv) strong correlations between geometric and hemodynamic indices highlight the need to integrate these factors for a comprehensive understanding of CAD initiation and progression in preclinical models. Thus, to optimize research based on murine models, it is essential not only to move beyond idealized geometries, but also to avoid uncritically relying on hemodynamic measurements from different species. This study grounds future development of mouse-specific predictive models of CAD, a critical step toward advancing translational research to understand and prevent CAD in humans.

The Effect of Arterial Elongation on Isolated Common Iliac Artery Pathologies

PURPOSE

to investigate the effects of vessel geometry on steno-occlusive and dilatative common iliac artery (CIA) pathologies.

METHODS

this single-center, retrospective study included 100 participants, namely 60 participants with a unilateral, isolated CIA pathology who were divided into three pathology-based groups (a stenosis group, n = 20, an occlusion group, n = 20, and an aneurysm group, n = 20) and 40 participants without a CIA pathology (control group). All participants underwent abdominal and pelvic computed tomography angiography. The aortoiliac region of the participants was reconstructed into three-dimensional models. Elongation parameters (tortuosity index (TI) and absolute average curvature (AAC)) and bifurcation parameters (iliac take-off angle, iliac planarity angle, and bifurcation angle) were determined using an in-house-written piece of software. Demographic data, anthropometric data, cardiovascular risk factor data, and medical history data were obtained from participants' electronic health records. The following statistical methods were used: one-way ANOVA, chi-square test, t-tests, Wilcoxon test, Kruskal-Wallis test, and multivariate linear regression.

RESULTS

in the occlusion group, both TI and AAC values were significantly higher on the contralateral side than on the ipsilateral side (both p < 0.001), whereas in the aneurysm group the AAC values were significantly higher on the ipsilateral side than on the contralateral side (p = 0.001). The ipsilateral and contralateral TI and AAC values of the iliac arteries were significantly higher in the aneurysm group than in the other three groups (all p < 0.001). Age significantly affected all of the elongation parameters except for the TI of the infrarenal aorta (all p < 0.010 except the TI of the infrarenal aorta). In addition, the AAC values for the iliac arteries were significantly associated with obesity (ipsilateral iliac artery, p = 0.045; contralateral iliac artery, p = 0.047). Aortic bifurcation parameters did not differ significantly either within each group (ipsilateral versus contralateral side) or between the individual groups.

CONCLUSIONS

occlusions tend to develop in relatively straight iliac arteries, whereas unilateral, isolated CIA aneurysms are more likely to occur in elongated aortoiliac systems.

Numerical investigation on the impact of different coronary aneurysms morphologies on thrombus formation and hemodynamics: a comparative study

Coronary artery aneurysms (CAAs) are morphologically classified as saccular and fusiform. There is still a great deal of clinical controversy as to which types of CAA are more likely to cause thrombosis. Therefore, the main objective of this study was to evaluate the trend of thrombus growth in CAAs with different morphologies and to assess the risk of possible long-term complications based on hemodynamic parameters. Utilizing computed tomography angiography (CTA) data from eight healthy coronary arteries, two distinct morphologies of coronary artery aneurysms (CAAs) were reconstructed. Distribution of four wall shear stress (WSS)-based indicators and three helicity indicators was analyzed in this study. Meanwhile, a thrombus growth model was introduced to analyze the thrombus formation in CAAs with different morphologies. The research results showed the distribution of most WSS indicators between saccular and fusiform CAAs was not statistically significant. However, due to the presence of a more pronounced helical flow pattern, irregular helical flow structure and longer time of flow stagnation in saccular CAAs during the cardiac cycle, the mean and maximum relative residence time (RRT) were significantly higher in saccular CAAs than in fusiform CAAs (P < 0.05). This may increase the risk of saccular coronary arteries leading to aneurysmal dilatation or even rupture. Although the two CAAs had similar rates of thrombosis, fusiform CAAs may more early cause obstruction of the main coronary flow channel where the aneurysm is located due to thrombosis growth. Thus, the risk of thrombosis in fusiform coronary aneurysms may warrant greater clinical concern.

Advances in the Computational Assessment of Disturbed Coronary Flow and Wall Shear Stress: A Contemporary Review

Coronary artery blood flow is influenced by various factors including vessel geometry, hemodynamic conditions, timing in the cardiac cycle, and rheological conditions. Multiple patterns of disturbed coronary flow may occur when blood flow separates from the laminar plane, associated with inefficient blood transit, and pathological processes modulated by the vascular endothelium in response to abnormal wall shear stress. Current simulation techniques, including computational fluid dynamics and fluid-structure interaction, can provide substantial detail on disturbed coronary flow and have advanced the contemporary understanding of the natural history of coronary disease. However, the clinical application of these techniques has been limited to hemodynamic assessment of coronary disease severity, with the potential to refine the assessment and management of coronary disease. Improved computational efficiency and large clinical trials are required to provide an incremental clinical benefit of these techniques beyond existing tools. This contemporary review is a clinically relevant overview of the disturbed coronary flow and its associated pathological consequences. The contemporary methods to assess disturbed flow are reviewed, including clinical applications of these techniques. Current limitations and future opportunities in the field are also discussed.

Numerical simulation and in vitro experimental study of the hemodynamic performance of vena cava filters with helical forms

Inferior vena cava filter (IVCF) implantation is a common method of thrombus capture. By implanting a filter in the inferior vena cava (IVC), microemboli can be effectively blocked from entering the pulmonary circulation, thereby avoiding acute pulmonary embolism (PE). Inspired by the helical flow effect in the human arterial system, we propose a helical retrievable IVCF, which, due to the presence of a helical structure inducing a helical flow pattern of blood in the region near the IVCF, can effectively avoid the deposition of microemboli in the vicinity of the IVCF while promoting the cleavage of the captured thrombus clot. It also reduces the risk of IVCF dislodging and slipping in the vessel because its shape expands in the radial direction, allowing its distal end to fit closely to the IVC wall, and because its contact structure with the inner IVC wall is curved, increasing the contact area and reducing the risk of the vessel wall being punctured by the IVCF support structure. We used ANSYS 2023 software to conduct unidirectional fluid-structure coupling simulation of four different forms of IVCF, combined with microthrombus capture experiments in vitro, to explore the impact of these four forms of IVCF on blood flow patterns and to evaluate the risk of IVCF perforation and IVCF dislocation. It can be seen from the numerical simulation results that the helical structure does have the function of inducing blood flow to undergo helical flow dynamics, and the increase in wall shear stress (WSS) brought about by this function can improve the situation of thrombosis accumulation to a certain extent. Meanwhile, the placement of IVCF will change the flow state of blood flow and lead to the deformation of blood vessels. In in vitro experiments, we found that the density of the helical support rod is a key factor affecting the thrombus trapping efficiency, and in addition, the contact area between the IVCF and the vessel wall has a major influence on the risk of IVCF displacement.

Shape-Setting of Self-Expanding Nickel-Titanium Laser-Cut and Wire-Braided Stents to Introduce a Helical Ridge

PURPOSE

Altered hemodynamics caused by the presence of an endovascular device may undermine the success of peripheral stenting procedures. Flow-enhanced stent designs are under investigation to recover physiological blood flow patterns in the treated artery and reduce long-term complications. However, flow-enhanced designs require the development of customised manufacturing processes that consider the complex behaviour of Nickel-Titanium (Ni-Ti). While the manufacturing routes of traditional self-expanding Ni-Ti stents are well-established, the process to introduce alternative stent designs is rarely reported in the literature, with much of this information (especially related to shape-setting step) being commercially sensitive and not reaching the public domain, as yet.

METHODS

A reliable manufacturing method was developed and improved to induce a helical ridge onto laser-cut and wire-braided Nickel-Titanium self-expanding stents. The process consisted of fastening the stent into a custom-built fixture that provided the helical shape, which was followed by a shape-setting in air furnace and rapid quenching in cold water. The parameters employed for the shape-setting in air furnace were thoroughly explored, and their effects assessed in terms of the mechanical performance of the device, material transformation temperatures and surface finishing.

RESULTS

Both stents were successfully imparted with a helical ridge and the optimal heat treatment parameters combination was found. The settings of 500 °C/30 min provided mechanical properties comparable with the original design, and transformation temperatures suitable for stenting applications (Af = 23.5 °C). Microscopy analysis confirmed that the manufacturing process did not alter the surface finishing. Deliverability testing showed the helical device could be loaded onto a catheter delivery system and deployed with full recovery of the expanded helical configuration.

CONCLUSION

This demonstrates the feasibility of an additional heat treatment regime to allow for helical shape-setting of laser-cut and wire-braided devices that may be applied to further designs.

Aortic geometry and long-term outcome in patients with a repaired coarctation

OBJECTIVE

This study aims to compare aortic morphology between repaired coarctation patients and controls, and to identify aortic morphological risk factors for hypertension and cardiovascular events (CVEs) in coarctation patients.

METHODS

Repaired coarctation patients with computed tomography angiography (CTA) or magnetic resonance angiography (MRA) were included, followed-up and compared with sex-matched and age-matched controls. Three-dimensional aortic shape was reconstructed using patients' CTA or MRA, or four-dimensional flow cardiovascular magnetic resonance in controls, and advanced geometrical characteristics were calculated and visualised using statistical shape modelling. In patients, we examined the association of geometrical characteristics with (1) baseline hypertension, using multivariable logistic regression; and (2) cardiovascular events (CVE, composite of aortic complications, coronary artery disease, ventricular arrhythmias, heart failure hospitalisation, stroke, transient ischaemic attacks and cardiovascular death), using multivariable Cox regression. The least absolute shrinkage and selection operator (LASSO) method selected the most informative multivariable model.

RESULTS

Sixty-five repaired coarctation patients (23 years (IQR 19-38)) were included, of which 44 (68%) patients were hypertensive at baseline. After a median follow-up of 8.7 years (IQR 4.8-15.4), 27 CVEs occurred in 20 patients. Aortic arch dimensions were smaller in patients compared with controls (diameter p<0.001, wall surface area p=0.026, volume p=0.007). Patients had more aortic arch torsion (p<0.001) and a higher curvature (p<0.001). No geometrical characteristics were associated with hypertension. LASSO selected left ventricular mass, male sex, tortuosity and age for the multivariable model. Left ventricular mass (p=0.014) was independently associated with CVE, and aortic tortuosity showed a trend towards significance (p=0.070).

CONCLUSION

Repaired coarctation patients have a smaller aortic arch and a more tortuous course of the aorta compared with controls. Besides left ventricular mass index, geometrical features might be of importance in long-term risk assessment in coarctation patients.

Predicting Lipid-Rich Plaque Progression in Coronary Arteries Using Multimodal Imaging and Wall Shear Stress Signatures

BACKGROUND

Plaque composition and wall shear stress (WSS) magnitude act as well-established players in coronary plaque progression. However, WSS magnitude per se does not completely capture the mechanical stimulus to which the endothelium is subjected, since endothelial cells experience changes in the WSS spatiotemporal configuration on the luminal surface. This study explores WSS profile and lipid content signatures of plaque progression to identify novel biomarkers of coronary atherosclerosis.

METHODS

Thirty-seven patients with acute coronary syndrome underwent coronary computed tomography angiography, near-infrared spectroscopy intravascular ultrasound, and optical coherence tomography of at least 1 nonculprit vessel at baseline and 1-year follow-up. Baseline coronary artery geometries were reconstructed from intravascular ultrasound and coronary computed tomography angiography and combined with flow information to perform computational fluid dynamics simulations to assess the time-averaged WSS magnitude (TAWSS) and the variability in the contraction/expansion action exerted by WSS on the endothelium, quantifiable in terms of topological shear variation index (TSVI). Plaque progression was measured as intravascular ultrasound-derived percentage plaque atheroma volume change at 1-year follow-up. Plaque composition information was extracted from near-infrared spectroscopy and optical coherence tomography.

RESULTS

Exposure to high TSVI and low TAWSS was associated with higher plaque progression (4.00±0.69% and 3.60±0.62%, respectively). Plaque composition acted synergistically with TSVI or TAWSS, resulting in the highest plaque progression (≥5.90%) at locations where lipid-rich plaque is exposed to high TSVI or low TAWSS.

CONCLUSIONS

Luminal exposure to high TSVI, solely or combined with a lipid-rich plaque phenotype, is associated with enhanced plaque progression at 1-year follow-up. Where plaque progression occurred, low TAWSS was also observed. These findings suggest TSVI, in addition to low TAWSS, as a potential biomechanical predictor for plaque progression, showing promise for clinical translation to improve patient prognosis.

Residence time in complex left main bifurcation disease after stenting

BACKGROUND

Data regarding the mean resident time (RT) after left main (LM) bifurcation stenting are scant. In the present study we performed a patient-specific computational fluid dynamic (CFD) analysis to investigate the different post-stenting mean RT values in LM patients treated with single-or double stenting techniques.

METHODS

Patients were identified after reviewing the local Optical Coherence Tomography (OCT) scans database. Overall, 27 patients (mean age 65.5 ± 12.4, 21 males) [10 patients treated with provisional cross-over stenting, 7 with the double kissing crush (DK crush) and 10 with the nano-inverted T (NIT) technique, respectively] with isolated and significant LM bifurcation disease were analyzed.

RESULTS

After LM bifurcation stenting, the NIT showed a higher averages WSS values at all bifurcation sites compared to DK crush and provisional cross-over stenting. Moreover, the mean RT resulted lower after NIT compared to provisional or DK crush. During the diastolic phase, the average RT of the entire LM bifurcation was 0.46 s, 0.38 s and 0.33 s after using the provisional stenting, DK crush and NIT, respectively. Moreover, the average RT in the LM bifurcation decreased by 17.1 % using the DK crush and by 28.2 % using the NIT compared to the Provisional.

CONCLUSION

The present OCT-derived CFD analysis revealed that, in patients with complex bifurcation LM disease, the provisional approach resulted in lower WSS values, while double stenting techniques, especially the NIT technique, resulted in a marked reduction of average RT compared to the provisional approach.

CONDENSED ABSTRACT

In the present study we performed a patient-specific Optical coherence tomography (OCT)-based computational fluid dynamic (CFD) analysis to investigate the different post-stenting mean RT values in 27 patients treated with provisional cross-over stenting, DK crush and Nano-inverted-T (NIT) stenting. The NIT showed a higher averages WSS values at all bifurcation sites compared to DK crush and Provisional. The mean RT resulted lower in NIT compared to Provisional or DK crush. During the entire diastolic phase, the average RT of the entire LM bifurcation was 0.46 s, 0.38 s and 0.33 s after using the provisional stenting, DK crush and NIT, respectively. Moreover, the average RT in the entire LM bifurcation decreased by 17.1 % using the DK crush and by 28.2 % using the NIT compared to the Provisional.

Impact of Cross-Limb Stent-Graft Configuration on Hemodynamics in Abdominal Aortic Aneurysm Interventional Therapy

PURPOSE

The cross-limb (CL) technique is a commonly used endovascular treatment for addressing unfavorable anatomical features in patients with abdominal aortic aneurysm (AAA). The configuration of CL stent-graft plays a critical role in determining the postoperative hemodynamic properties and physiological behaviors, which ultimately impact the efficacy and safety of endovascular AAA treatment. This study aims to investigate the relationship between hemodynamics and CL stent-graft configuration from a hemodynamic perspective.

METHODS

Five distinct geometric models of cross-limb (CL) stent-graft configurations were constructed by optimizing the real clinical computed tomography angiography (CTA) data. These models varied in main body lengths and cross angles and were used to perform numerical simulations to analyze various hemodynamic parameters. Flow pattern, distribution of wall shear stress (WSS)-related parameters, localized normalized helicity (LNH), pressure drop, and the displacement force of all models were examined in this paper.

RESULTS

In patient-specific cases, helical flow and WSS increase with the main body. However, it also generated secondary flow in localized areas, leading to increased oscillation in the WSS direction. Notably, increasing the stent graft's main body length or decreasing the cross angle reduced the displacement force exerted on the stent-graft. Reducing the cross angle did not significantly alter the hemodynamic characteristics.

CONCLUSION

In the clinical practice of CL deployment, it is crucial to carefully consider the stent-graft configuration and the patient specific to achieve optimal postoperative outcomes. This study provides valuable insights for guiding stent selection and treatment planning in patients with abdominal aortic aneurysm undergoing CL techniques, from a hemodynamic perspective.

Analyzing the effects of helical flow in blood vessels using acoustofluidic-based dynamic flow generator

The effects of helical flow in a blood vessel are investigated in a dynamic flow generator using surface acoustic wave (SAW) in the microfluidic device. The SAW, generated by an interdigital transducer (IDT), induces acoustic streaming, resulting in a stable and consistent helical flow pattern in microscale channels. This approach allows rapid development of helical flow within the channel without directly contacting the medium. The precise design of the window enables the creation of distinct unidirectional vortices, which can be controlled by adjusting the amplitude of the SAW. Within this device, optimal operational parameters of the dynamic flow generator to preserve the integrity of endothelial cells are found, and in such settings, the actin filaments within the cells are aligned to the desired state. Our findings reveal that intracellular Ca2+ concentrations vary in response to flow conditions. Specifically, comparable maximum intensity and graphical patterns were observed between low-flow rate helical flow and high-flow rate Hagen-Poiseuille flow. These suggest that the cells respond to the helical flow through mechanosensitive ion channels. Finally, adherence of monocytes is effectively reduced under helical flow conditions in an inflammatory environment, highlighting the atheroprotective role of helical flow. STATEMENT OF SIGNIFICANCE: Helical flow in blood vessels is well known to prevent atherosclerosis. However, despite efforts to replicate helical flow in microscale channels, there is still a lack of in vitro models which can generate helical flow for analyzing its effects on the vascular system. In this study, we developed a method for generating steady and constant helical flow in microfluidic channel using acoustofluidic techniques. By utilizing this dynamic flow generator, we were able to observe the atheroprotective aspects of helical flow in vitro, including the enhancement of calcium ion flux and reduction of monocyte adhesion. This study paves the way for an in vitro model of dynamic cell culture and offers advanced investigation into helical flow in our circulatory system.

Design of experiment approach to identify the dominant geometrical feature of left coronary artery influencing atherosclerosis

Background and Objective. Coronary artery geometry heavily influences local hemodynamics, potentially leading to atherosclerosis. Consequently, the unique geometrical configuration of an individual by birth can be associated with future risk of atherosclerosis. Although current researches focus on exploring the relationship between local hemodynamics and coronary artery geometry, this study aims to identify the order of influence of the geometrical features through systematic experiments, which can reveal the dominant geometrical feature for future risk assessment.Methods. According to Taguchi's method of design of experiment (DoE), the left main stem (LMS) length (lLMS), curvature (kLMS), diameter (dLMS) and the bifurcation angle between left anterior descending (LAD) and left circumflex (LCx) artery (αLAD-LCx) of two reconstructed patient-specific left coronary arteries (LCA) were varied in three levels to create L9 orthogonal array. Computational fluid dynamic (CFD) simulations with physiological boundary conditions were performed on the resulting eighteen LCA models. Average helicity intensity (h2) and relative atheroprone area (RAA) of near-wall hemodynamic descriptors were analyzed.Results. The proximal LAD (LADproximal) was identified to be the most atheroprone region of the left coronary artery due to higherh2,large RAA of time averaged wall shear stress (TAWSS < 0.4 Pa), oscillatory shear index (OSI ∼ 0.5) and relative residence time (RRT > 4.17 Pa-1). In both patient-specific cases, based onh2and TAWSS,dlmsis the dominant geometric parameter while based on OSI and RRT,αLAD-LCxis the dominant one influencing hemodynamic condition in proximal LAD (p< 0.05). Based on RRT, the rank of the geometrical factors is:αLAD-LCx>dLMS>lLMS>kLMS, indicating thatαLAD-LCxis the most dominant geometrical factor affecting hemodynamics at proximal LAD which may influence atherosclerosis.Conclusion. The proposed identification of the rank of geometrical features of LCA and the dominant feature may assist clinicians in predicting the possibility of atherosclerosis, of an individual, long before it will occur. This study can further be translated to be used to rank the influence of several arterial geometrical features at different arterial locations to explore detailed relationships between the arterial geometrical features and local hemodynamics.

Blood Flow Energy Identifies Coronary Lesions Culprit of Future Myocardial Infarction

The present study establishes a link between blood flow energy transformations in coronary atherosclerotic lesions and clinical outcomes. The predictive capacity for future myocardial infarction (MI) was compared with that of established quantitative coronary angiography (QCA)-derived predictors. Angiography-based computational fluid dynamics (CFD) simulations were performed on 80 human coronary lesions culprit of MI within 5 years and 108 non-culprit lesions for future MI. Blood flow energy transformations were assessed in the converging flow segment of the lesion as ratios of kinetic and rotational energy values (KER and RER, respectively) at the QCA-identified minimum lumen area and proximal lesion sections. The anatomical and functional lesion severity were evaluated with QCA to derive percentage area stenosis (%AS), vessel fractional flow reserve (vFFR), and translesional vFFR (ΔvFFR). Wall shear stress profiles were investigated in terms of topological shear variation index (TSVI). KER and RER predicted MI at 5 years (AUC = 0.73, 95% CI 0.65-0.80, and AUC = 0.76, 95% CI 0.70-0.83, respectively; p < 0.0001 for both). The predictive capacity for future MI of KER and RER was significantly stronger than vFFR (p = 0.0391 and p = 0.0045, respectively). RER predictive capacity was significantly stronger than %AS and ΔvFFR (p = 0.0041 and p = 0.0059, respectively). The predictive capacity for future MI of KER and RER did not differ significantly from TSVI. Blood flow kinetic and rotational energy transformations were significant predictors for MI at 5 years (p < 0.0001). The findings of this study support the hypothesis of a biomechanical contribution to the process of plaque destabilization/rupture leading to MI.

Modelling blood flow in coronary arteries: Newtonian or shear-thinning non-Newtonian rheology?

BACKGROUND

The combination of medical imaging and computational hemodynamics is a promising technology to diagnose/prognose coronary artery disease (CAD). However, the clinical translation of in silico hemodynamic models is still hampered by assumptions/idealizations that must be introduced in model-based strategies and that necessarily imply uncertainty. This study aims to provide a definite answer to the open question of how to properly model blood rheological properties in computational fluid dynamics (CFD) simulations of coronary hemodynamics.

METHODS

The geometry of the right coronary artery (RCA) of 144 hemodynamically stable patients with different stenosis degree were reconstructed from angiography. On them, unsteady-state CFD simulations were carried out. On each reconstructed RCA two different simulation strategies were applied to account for blood rheological properties, implementing (i) a Newtonian (N) and (ii) a shear-thinning non-Newtonian (non-N) rheological model. Their impact was evaluated in terms of wall shear stress (WSS magnitude, multidirectionality, topological skeleton) and helical flow (strength, topology) profiles. Additionally, luminal surface areas (SAs) exposed to shear disturbances were identified and the co-localization of paired N and non-N SAs was quantified in terms of similarity index (SI).

RESULTS

The comparison between paired N vs. shear-thinning non-N simulations revealed remarkably similar profiles of WSS-based and helicity-based quantities, independent of the adopted blood rheology model and of the degree of stenosis of the vessel. Statistically, for each paired N and non-N hemodynamic quantity emerged negligible bias from Bland-Altman plots, and strong positive linear correlation (r > 0.94 for almost all the WSS-based quantities, r > 0.99 for helicity-based quantities). Moreover, a remarkable co-localization of N vs. non-N luminal SAs exposed to disturbed shear clearly emerged (SI distribution 0.95 [0.93, 0.97]). Helical flow topology resulted to be unaffected by blood rheological properties.

CONCLUSIONS

This study, performed on 288 angio-based CFD simulations on 144 RCA models presenting with different degrees of stenosis, suggests that the assumptions on blood rheology have negligible impact both on WSS and helical flow profiles associated with CAD, thus definitively answering to the question "is Newtonian assumption for blood rheology adequate in coronary hemodynamics simulations?".

Computational fluid dynamics as supporting technology for coronary artery disease diagnosis and treatment: an international survey

Background

Computational fluid dynamics (CFD) is emerging as an effective technology able to improve procedural outcomes and enhance clinical decision-making in patients with coronary artery disease (CAD). The present study aims to assess the state of knowledge, use and clinical acceptability of CFD in the diagnosis and treatment of CAD.

Methods

We realized a 20-questions international, anonymous, cross-sectional survey to cardiologists to test their knowledge and confidence on CFD as a technology applied to patients suffering from CAD. Responses were recorded between May 18, 2022, and June 12, 2022.

Results

A total of 466 interventional cardiologists (mean age 48.4 ± 8.3 years, males 362), from 42 different countries completed the survey, for a response rate of 45.9%. Of these, 66.6% declared to be familiar with the term CFD, especially for optimization of existing interventional techniques (16.1%) and assessment of hemodynamic quantities related with CAD (13.7%). About 30% of respondents correctly answered to the questions exploring their knowledge on the pathophysiological role of some CFD-derived quantities such as wall shear stress and helical flow in coronary arteries. Among respondents, 85.9% would consider patient-specific CFD-based analysis in daily interventional practice while 94.2% declared to be interested in receiving a brief foundation course on the basic CFD principles. Finally, 87.7% of respondents declared to be interested in a cath-lab software able to conduct affordable CFD-based analyses at the point-of-care.

Conclusions

Interventional cardiologists reported to be profoundly interested in adopting CFD simulations as a technology supporting decision making in the treatment of CAD in daily practice.

In vivo performance of a tri-leaflet mechanical heart valve prosthesis in an ovine model

OBJECTIVES

We present the long-term results of a trileaflet (Triflo) versus bileaflet (On-X) mechanical valve in both aortic and pulmonary positions in a sheep model.

METHODS

The Triflo valve was implanted in 21 female sheep in aortic (n = 8) and pulmonary position (n = 13). The On-X valve was implanted in 7 female sheep in aortic (n = 1) and pulmonary (n = 6) positions. No antithrombotic medication of any kind was given postoperatively. In the aortic group, survival cohorts were 3 and 5 months. In the pulmonary group, survival cohorts were 10 and 20 weeks. Valve performance was assessed using haematology, echocardiography and acoustic measurements combined with post-mortem pathology analysis of the downstream organs.

RESULTS

The mean gradients were lower for the Triflo valve in both pulmonary [4.30 mmHg (3.70-5.73) vs 6.80 mmHg (4.63-7.96), P = 0.012] and aortic [5.1 mmHg (4.2-7.7) vs 10.7 mmHg (8.7-12.9), P = 0.007] positions. Peak gradients were lower for the Triflo valve in both pulmonary [8.05 mmHg (6.75-10.23) vs 13.15 mmHg (9.20-14.76), P = 0.005] and aortic [8.7 mmHg (7.5-12.5) vs 16.5 mmHg (14.2-19.6), P = 0.009] positions. In both positions, leaflets and housing surface were free from any deposits macro- and microscopically and comparable to nonimplanted control valves. Peripheral organs showed no signs of thrombo-embolic damage. Biochemical and haematological were comparable to preoperative. The closing click sound pressure level of the Triflo was significantly lower in both aortic [108.4 sound pressure level (102.0-115.7) vs 111.7 sound pressure level (105.5-117.0), P < 0.001] and pulmonary [103.6 sound pressure level (99.1-108.9) vs 118.5 sound pressure level (116.7-120.2), P < 0.001] position.

CONCLUSIONS

Preliminary in vivo results of the Triflo valve are promising in both aortic and pulmonary positions in an ovine model. Excellent haemodynamics, stable long-term function, low valve noise and no thrombo-embolic events in the absence of antithrombotic medication lay the foundation to a future clinical first-in-man trial.

Risk evaluation of adverse aortic events in patients with non-circular aortic annulus after transcatheter aortic valve implantation: a numerical study

Transcatheter aortic valve implantation (TAVI) is a micro-invasive surgery used to treat patients with aortic stenosis (AS) efficiently. However, the uneven valve expansion can cause a non-circular annulus, which is one of the main factors leading to complications after TAVI. As a preliminary work, the main purpose of this study was to evaluate the risk of adverse aortic events in patients with a non-circular aortic annulus after TAVI. This study numerically investigated the distribution of four wall shear stress (WSS)-based indicators and three helicity-based indicators in eight patient-specific aortas with different annulus including circular, type I elliptical and type II elliptical shapes. Both elliptical annulus features can significantly enhance the intensity of the helicity (h2) in the ascending aorta (p < 0.001). However, for the type I elliptical annulus, the spiral flow structure was changed into low-velocity and disturbed flow pattern close to the inner side of the aortic arch. For the type II elliptical annulus, the spiral flow remained but became skewed in distribution. The elliptical annulus feature could increase the general level WSS-based indicators, especially in the ascending aorta. However, due to the disturbance of spiral flow or second helical flow in ascending aortas, areas with low TAWSS accompanied by high oscillatory shear index (OSI) and cross flow index (CFI) were observed in all the ascending aortas with non-circular annulus. The elliptical annulus feature can change the hemodynamic environment in the aortic arch, especially in the ascending aorta. Although both elliptical annulus features enhanced the strength of helicity, the uniform distribution of the helical flow was disturbed, especially in the ascending aorta, indicating the potential risk of adverse aortic events may increase. Therefore, for the patients without paravalvular leak but elliptical annulus shape after TAVI treatment, surgeons may be needed to consider further dilatation to make the non-circular annulus become circular.

On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis

BACKGROUND AND OBJECTIVE

In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied.

METHODS

We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology.

RESULTS

22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately.

CONCLUSIONS

These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.

Spiral Laminar Flow is Associated with a Reduction in Disturbed Shear in Patient-Specific Models of an Arteriovenous Fistula

PURPOSE

Areas of disturbed shear that arise following arteriovenous fistula (AVF) creation are believed to contribute to the development of intimal hyperplasia (IH). The presence of helical flow can suppress areas of disturbed shear, which may protect the vasculature from IH. Therefore, the aim of this study is to determine if helical flow, specifically spiral laminar flow (SLF), is present in patient-specific AVF models and is associated with a reduction in exposure to disturbed shear.

METHODS

Four AVF were imaged using MRI within the first two weeks following fistula creation. Patient-specific boundary conditions were obtained using phase-contrast MRI and applied at the inlet and outlets of each model. Computational fluid dynamics was used to analyse the hemodynamics in each model and compare the helical content of the flow to the distribution of disturbed shear.

RESULTS

BC-1 and RC-2 are characterised by the presence of SLF, which coincides with the lowest distribution of disturbed shear. Contrastingly, SLF is absent from BC-2 and RC-1 and experience the largest amount of disturbed shear. Interestingly, BC-2 and RC-1 developed an anastomosis stenosis, while BC-1 and RC-2 remained stenosis free.

CONCLUSION

These findings are in agreement with previous clinical studies and further highlight the clinical potential of SLF as a prognostic marker for a healthy AVF, as its presence correlates with an overall reduction in exposure to disturbed shear and a decrease in the incidence of AVF dysfunction, albeit in a small sample size.

Helical flow in tortuous aortas and its relationship to turbulence: A whole-aorta 4D flow MRI study

Background

Increased vascular tortuosity is a hallmark of ageing of the vascular system, including the aorta. However, the impact of tortuosity on aortic blood flow is unknown. We hypothesized that increased tortuosity would be associated with increased blood flow helicity and with decreased degree of blood flow turbulence as measured by the turbulent kinetic energy (TKE).

Methods

4D Flow MR images covering the entire aorta from the aortic valve to the iliac bifurcation were acquired in 23 normal volunteers aged 18-30 years ("Young") and 23 normal volunteers aged 66-76 years ("Old") without aortic disease. The aorta was segmented and divided into four regions: the ascending, descending, suprarenal abdominal and infrarenal abdominal aorta. Tortuosity, helicity, TKE, flow velocity, and Reynolds number were computed for the whole aorta and for each section.

Results

Tortuosity and helicity were higher whereas TKE, velocity, and Reynolds number were lower in Old than in Young, for all aortic regions (p < 0.05) except for helicity in the descending aorta. Tortuosity correlated positively with helicity and negatively with TKE for all aortic regions (Spearman rho=±0.45-±0.72, p < =0.002) except for TKE in the ascending aorta. Further, helicity correlated with TKE in the descending, suprarenal abdominal and infrarenal abdominal aorta (Spearman rho=-0.56--0.77).

Conclusion

Tortuosity increases with age and blood flow in tortuous aortas is more helical. Increasing helicity, in turn, is associated with decreasing TKE.

Wall shear stress-related plaque growth of lipid-rich plaques in human coronary arteries: an near-infrared spectroscopy and optical coherence tomography study

AIMS

Low wall shear stress (WSS) is acknowledged to play a role in plaque development through its influence on local endothelial function. Also, lipid-rich plaques (LRPs) are associated with endothelial dysfunction. However, little is known about the interplay between WSS and the presence of lipids with respect to plaque progression. Therefore, we aimed to study the differences in WSS-related plaque progression between LRPs, non-LRPs, or plaque-free regions in human coronary arteries.

METHODS AND RESULTS

In the present single-centre, prospective study, 40 patients who presented with an acute coronary syndrome successfully underwent near-infrared spectroscopy intravascular ultrasound (NIRS-IVUS) and optical coherence tomography (OCT) of at least one non-culprit vessel at baseline and completed a 1-year follow-up. WSS was computed applying computational fluid dynamics to a three-dimensional reconstruction of the coronary artery based on the fusion of the IVUS-segmented lumen with a CT-derived centreline, using invasive flow measurements as boundary conditions. For data analysis, each artery was divided into 1.5 mm/45° sectors. Plaque growth based on IVUS-derived percentage atheroma volume change was compared between LRPs, non-LRPs, and plaque-free wall segments, as assessed by both OCT and NIRS. Both NIRS- and OCT-detected lipid-rich sectors showed a significantly higher plaque progression than non-LRPs or plaque-free regions. Exposure to low WSS was associated with a higher plaque progression than exposure to mid or high WSS, even in the regions classified as a plaque-free wall. Furthermore, low WSS and the presence of lipids had a synergistic effect on plaque growth, resulting in the highest plaque progression in lipid-rich regions exposed to low shear stress.

CONCLUSION

This study demonstrates that NIRS- and OCT-detected lipid-rich regions exposed to low WSS are subject to enhanced plaque growth over a 1-year follow-up. The presence of lipids and low WSS proves to have a synergistic effect on plaque growth.

Helical Flow in Healthy and Diseased Patient-specific Coronary Bifurcations

Helical flow (HF) exists in healthy and diseased coronary bifurcations and was found to have a protective atherosclerotic vascular effect in other vessels. However, the role of HF in patient-specific human coronary arteries still needs further study, and is therefore the objective of this study in both healthy and diseased bifurcations. Computational studies were conducted on 16 patient-specific coronary bifurcations, including eight healthy and eight identical cases with idealized narrowing to represent disease. In general, higher HF intensity may have a favorable effect as it corelated to the reduction of the percentage vessel area exposed to adverse time averaged wall shear stress (TAWSS%) in both healthy and diseased models. The HF intensity and distribution of each model varies due to the complex shape of patient-specific models. The presence of disease appears to have an important impact on the downstream HF patterns and the TAWSS distributions. Clinical Relevance- By understanding the relationship between HF and hemodynamics, HF may be used as a predictor for the formation and progression of atherosclerotic plaque in coronary arteries instead of near-wall WSS measures, which can be determined with higher accuracy in vivo.

Modelling coronary flows: impact of differently measured inflow boundary conditions on vessel-specific computational hemodynamic profiles

BACKGROUND AND OBJECTIVES

The translation of hemodynamic quantities based on wall shear stress (WSS) or intravascular helical flow into clinical biomarkers of coronary atherosclerotic disease is still hampered by the assumptions/idealizations required by the computational fluid dynamics (CFD) simulations of the coronary hemodynamics. In the resulting budget of uncertainty, inflow boundary conditions (BCs) play a primary role. Accordingly, in this study we investigated the impact of the approach adopted for in vivo coronary artery blood flow rate assessment on personalized CFD simulations where blood flow rate is used as inflow BC.

METHODS

CFD simulations were carried out on coronary angiograms by applying personalized inflow BCs derived from four different techniques assessing in vivo surrogates of flow rate: continuous thermodilution, intravascular Doppler, frame count-based 3D contrast velocity, and diameter-based scaling law. The impact of inflow BCs on coronary hemodynamics was evaluated in terms of WSS- and helicity-based quantities.

RESULTS

As main findings, we report that: (i) coronary flow rate values may differ based on the applied flow derivation technique, as continuous thermodilution provided higher flow rate values than intravascular Doppler and diameter-based scaling law (p = 0.0014 and p = 0.0023, respectively); (ii) such intrasubject differences in flow rate values lead to different surface-averaged values of WSS magnitude and helical blood flow intensity (p<0.0020); (iii) luminal surface areas exposed to low WSS and helical flow topological features showed robustness to the flow rate values.

CONCLUSIONS

Although the absence of a clinically applicable gold standard approach prevents a general recommendation for one coronary blood flow rate derivation technique, our findings indicate that the inflow BC may impact computational hemodynamic results, suggesting that a standardization would be desirable to provide comparable results among personalized CFD simulations of the coronary hemodynamics.

Non-Newtonian Endothelial Shear Stress Simulation: Does It Matter?

Patient-specific coronary endothelial shear stress (ESS) calculations using Newtonian and non-Newtonian rheological models were performed to assess whether the common assumption of Newtonian blood behavior offers similar results to a more realistic but computationally expensive non-Newtonian model. 16 coronary arteries (from 16 patients) were reconstructed from optical coherence tomographic (OCT) imaging. Pulsatile CFD simulations using Newtonian and the Quemada non-Newtonian model were performed. Endothelial shear stress (ESS) and other indices were compared. Exploratory indices including local blood viscosity (LBV) were calculated from non-Newtonian simulation data. Compared to the Newtonian results, the non-Newtonian model estimates significantly higher time-averaged ESS (1.69 (IQR 1.36)Pa versus 1.28 (1.16)Pa, p < 0.001) and ESS gradient (0.90 (1.20)Pa/mm versus 0.74 (1.03)Pa/mm, p < 0.001) throughout the cardiac cycle, under-estimating the low ESS (<1Pa) area (37.20 ± 13.57% versus 50.43 ± 14.16%, 95% CI 11.28-15.18, p < 0.001). Similar results were also found in the idealized artery simulations with non-Newtonian median ESS being higher than the Newtonian median ESS (healthy segments: 0.8238Pa versus 0.6618Pa, p < 0.001 proximal; 0.8179Pa versus 0.6610Pa, p < 0.001 distal; stenotic segments: 0.8196Pa versus 0.6611Pa, p < 0.001 proximal; 0.2546Pa versus 0.2245Pa, p < 0.001 distal) On average, the non-Newtonian model has a LBV of 1.45 times above the Newtonian model with an average peak LBV of 40-fold. Non-Newtonian blood model estimates higher quantitative ESS values than the Newtonian model. Incorporation of non-Newtonian blood behavior may improve the accuracy of ESS measurements. The non-Newtonian model also allows calculation of exploratory viscosity-based hemodynamic indices, such as local blood viscosity, which may offer additional information to detect underlying atherosclerosis.

Carotid Geometry and Wall Shear Stress Independently Predict Increased Wall Thickness-A Longitudinal 3D MRI Study in High-Risk Patients

Introduction: Carotid geometry and wall shear stress (WSS) have been proposed as independent risk factors for the progression of carotid atherosclerosis, but this has not yet been demonstrated in larger longitudinal studies. Therefore, we investigated the impact of these biomarkers on carotid wall thickness in patients with high cardiovascular risk.

Methods: Ninety-seven consecutive patients with hypertension, at least one additional cardiovascular risk factor and internal carotid artery (ICA) plaques (wall thickness ≥ 1.5 mm and degree of stenosis ≤ 50%) were prospectively included. They underwent high-resolution 3D multi-contrast and 4D flow MRI at 3 Tesla both at baseline and follow-up. Geometry (ICA/common carotid artery (CCA)-diameter ratio, bifurcation angle, tortuosity and wall thickness) and hemodynamics [WSS, oscillatory shear index (OSI)] of both carotid bifurcations were measured at baseline. Their predictive value for changes of wall thickness 12 months later was calculated using linear regression analysis for the entire study cohort (group 1, 97 patients) and after excluding patients with ICA stenosis ≥10% to rule out relevant inward remodeling (group 2, 61 patients).

Results: In group 1, only tortuosity at baseline was independently associated with carotid wall thickness at follow-up (regression coefficient = −0.52, p < 0.001). However, after excluding patients with ICA stenosis ≥10% in group 2, both ICA/CCA-ratio (0.49, p < 0.001), bifurcation angle (0.04, p = 0.001), tortuosity (−0.30, p = 0.040), and WSS (−0.03, p = 0.010) at baseline were independently associated with changes of carotid wall thickness at follow-up.

Conclusions: A large ICA bulb and bifurcation angle and low WSS seem to be independent risk factors for the progression of carotid atherosclerosis in the absence of ICA stenosis. By contrast, a high carotid tortuosity seems to be protective both in patients without and with ICA stenosis. These biomarkers may be helpful for the identification of patients who are at particular risk of wall thickness progression and who may benefit from intensified monitoring and treatment.

The definition of low wall shear stress and its effect on plaque progression estimation in human coronary arteries

Wall shear stress (WSS), the frictional force of the blood on the vessel wall, plays a crucial role in atherosclerotic plaque development. Low WSS has been associated with plaque growth, however previous research used different approaches to define low WSS to investigate its effect on plaque progression. In this study, we used four methodologies to allocate low, mid and high WSS in one dataset of human coronary arteries and investigated the predictive power of low WSS for plaque progression. Coronary reconstructions were based on multimodality imaging, using intravascular ultrasound and CT-imaging. Vessel-specific flow was measured using Doppler wire and computational fluid dynamics was performed to calculate WSS. The absolute WSS range varied greatly between the coronary arteries. On the population level, the established pattern of most plaque progression at low WSS was apparent in all methodologies defining the WSS categories. However, for the individual patient, when using measured flow to determine WSS, the absolute WSS values range so widely, that the use of absolute thresholds to determine low WSS was not appropriate to identify regions at high risk for plaque progression.

Numerical modelling of blood rheology and platelet activation through a stenosed left coronary artery bifurcation

Coronary bifurcations are prone to atherosclerotic plaque growth, experiencing regions of reduced wall shear stress (WSS) and increased platelet adhesion. This study compares effects across different rheological approaches on hemodynamics, combined with a shear stress exposure history model of platelets within a stenosed porcine bifurcation. Simulations used both single/multiphase blood models to determine which approach best predicts phenomena associated with atherosclerosis and atherothrombosis. A novel Lagrangian platelet tracking model was used to evaluate residence time and shear history of platelets indicating likely regions of thrombus formation. Results show a decrease in area of regions with pathologically low time-averaged WSS with the use of multiphase models, particularly in a stenotic bifurcation. Significant non-Newtonian effects were observed due to low-shear and varying hematocrit levels found on the outer walls of the bifurcation and distal to the stenosis. Platelet residence time increased 11% in the stenosed artery, with exposure times to low-shear sufficient for red blood cell aggregation (>1.5 s). increasing the risk of thrombosis. This shows stenotic artery hemodynamics are inherently non-Newtonian and multiphase, with variations in hematocrit (0.163-0.617) and elevated vorticity distal to stenosis (+15%) impairing the function of the endothelium via reduced time-averaged WSS regions, rheological properties and platelet activation/adhesion.

The Hydrodynamics of a Swirling Blood Flow in the Left Heart and Aorta

This paper proposes a new approach to the quantitative analysis of the hydrodynamic structure of a blood flow in the flow channel running from the left atrium to the end of the aorta. This approach is based on the concept of the structural organization of tornado-like swirling jets in channels with a given geometric configuration. Considering the large amount of experimental data in our possession, it was shown that along the entire length of the flow channel, conditions exist for the generation and maintenance of a swirling structure of the jet throughout the entire cardiac cycle. This study has given rise to a new direction in research in fundamental physiology and medicine, which is of great practical importance for diagnosing and treating circulatory disorders accompanied by changes in the geometric configuration and biomechanical characteristics of the heart and great vessels.

The effects of plaque morphological characteristics on the post-stenotic flow in left main coronary artery bifurcation

Local post-stenotic hemodynamics has critical influence in the atherosclerotic plaque progression occurring in susceptible arterial sites, in particular the left main coronary artery (LMCA) bifurcation. Understanding the effects of plaque morphological characteristics: stenosis severity (SS), eccentricity index (EI) and lesion length (LL) on the post-stenotic flow behavior can significantly improve treatment planning. In order to investigate these effects, we have employed computational fluid dynamics (CFD) simulations in twenty computer-generated and five patient-specific LMCA models and the hemodynamic parameters: velocity, pressure (P), wall pressure gradient (WPG), wall shear stress (WSS), time averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT) and helicity intensity (h₂) were analyzed. Our results revealed that the effect of stenosis eccentricity varied significantly for different values of stenosis severity and lesion length. Regions with low WSS, low TAWSS and high RRT were more prominent in models having higher stenosis severity. For smaller lesion length, at low and moderate stenosis severity, surface area with low TAWSS and high RRT decreased with increasing eccentricity index, whereas for high stenosis severity models, low TAWSS region and average RRT values increased with eccentricity. However, for models with longer lesion length, regions with high OSI and RRT overall increased gradually with eccentricity. The helicity intensity (h₂) of all models remained very low except at the most eccentric model with longer lesion length. The presence of very high helical flow in this model suggests the possibility of atheroprotective flow. It can be concluded that all plaque morphological characteristics covered under this investigation play an important role in plaque progression.

Comparison of Swine and Human Computational Hemodynamics Models for the Study of Coronary Atherosclerosis

Coronary atherosclerosis is a leading cause of illness and death in Western World and its mechanisms are still non completely understood. Several animal models have been used to 1) study coronary atherosclerosis natural history and 2) propose predictive tools for this disease, that is asymptomatic for a long time, aiming for a direct translation of their findings to human coronary arteries. Among them, swine models are largely used due to the observed anatomical and pathophysiological similarities to humans. However, a direct comparison between swine and human models in terms of coronary hemodynamics, known to influence atherosclerotic onset/development, is still lacking. In this context, we performed a detailed comparative analysis between swine- and human-specific computational hemodynamic models of coronary arteries. The analysis involved several near-wall and intravascular flow descriptors, previously emerged as markers of coronary atherosclerosis initiation/progression, as well as anatomical features. To do that, non-culprit coronary arteries (18 right–RCA, 18 left anterior descending–LAD, 13 left circumflex–LCX coronary artery) from patients presenting with acute coronary syndrome were imaged by intravascular ultrasound and coronary computed tomography angiography. Similarly, the three main coronary arteries of ten adult mini-pigs were also imaged (10 RCA, 10 LAD, 10 LCX). The geometries of the imaged coronary arteries were reconstructed (49 human, 30 swine), and computational fluid dynamic simulations were performed by imposing individualized boundary conditions. Overall, no relevant differences in 1) wall shear stress-based quantities, 2) intravascular hemodynamics (in terms of helical flow features), and 3) anatomical features emerged between human- and swine-specific models. The findings of this study strongly support the use of swine-specific computational models to study and characterize the hemodynamic features linked to coronary atherosclerosis, sustaining the reliability of their translation to human vascular disease.

In Vivo Based Fluid-Structure Interaction Biomechanics of the Left Anterior Descending Coronary Artery

A fluid-structure interaction-based biomechanical model of the entire left anterior descending coronary artery is developed from in vivo imaging via the finite element method in this paper. Included in this investigation is ventricle contraction, three-dimensional motion, all angiographically visible side branches, hyper/viscoelastic artery layers, non-Newtonian and pulsatile blood flow, and the out-of-phase nature of blood velocity and pressure. The fluid-structure interaction model is based on in vivo angiography of an elite athlete's entire left anterior descending coronary artery where the influence of including all alternating side branches and the dynamical contraction of the ventricle is investigated for the first time. Results show the omission of side branches result in a 350% increase in peak wall shear stress and a 54% decrease in von Mises stress. Peak von Mises stress is underestimated by up to 80% when excluding ventricle contraction and further alterations in oscillatory shear indices are seen, which provide an indication of flow reversal and has been linked to atherosclerosis localization. Animations of key results are also provided within a video abstract. We anticipate that this model and results can be used as a basis for our understanding of the interaction between coronary and myocardium biomechanics. It is hoped that further investigations could include the passive and active components of the myocardium to further replicate in vivo mechanics and lead to an understanding of the influence of cardiac abnormalities, such as arrythmia, on coronary biomechanical responses.

Early Atherosclerotic Changes in Coronary Arteries are Associated with Endothelium Shear Stress Contraction/Expansion Variability

Although unphysiological wall shear stress (WSS) has become the consensus hemodynamic mechanism for coronary atherosclerosis, the complex biomechanical stimulus affecting atherosclerosis evolution is still undetermined. This has motivated the interest on the contraction/expansion action exerted by WSS on the endothelium, obtained through the WSS topological skeleton analysis. This study tests the ability of this WSS feature, alone or combined with WSS magnitude, to predict coronary wall thickness (WT) longitudinal changes. Nine coronary arteries of hypercholesterolemic minipigs underwent imaging with local WT measurement at three time points: baseline (T1), after 5.6 ± 0.9 (T2), and 7.6 ± 2.5 (T3) months. Individualized computational hemodynamic simulations were performed at T1 and T2. The variability of the WSS contraction/expansion action along the cardiac cycle was quantified using the WSS topological shear variation index (TSVI). Alone or combined, high TSVI and low WSS significantly co-localized with high WT at the same time points and were significant predictors of thickening at later time points. TSVI and WSS magnitude values in a physiological range appeared to play an atheroprotective role. Both the variability of the WSS contraction/expansion action and WSS magnitude, accounting for different hemodynamic effects on the endothelium, (1) are linked to WT changes and (2) concur to identify WSS features leading to coronary atherosclerosis.

Numerical analysis of the hemodynamics of rat aorta based on magnetic resonance imaging and fluid-structure interaction

Murine models have been widely used to investigate the mechanobiology of aortic atherosclerosis and dissections, which develop preferably at different anatomic locations of aorta. Based MRI and finite element analysis with fluid-structure interaction, we numerically investigated factors that may affect the blood flow and structural mechanics of rat aorta. The results indicated that aortic root motion greatly increases time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), displacement of the aorta, and enhances helical flow pattern but has limited influence on effective stress, which is highly modulated by blood pressure. Moreover, the influence of the motion component on these indicators is different with axial motion more obvious than planar motion. Surrounding fixation of the intercostal arteries and the branch vessels on aortic arch would reduce the influence of aortic root motion. The compliance of the aorta has different influences at different regions, leading to decrease in TAWSS and helical flow, increase in OSI, RRT at the aortic arch, but has reversed effects on the branch vessels. When compared with the steady flow, the pulsatile blood flow would obviously increase the WSS, the displacement, and the effective stress in most regions. In conclusion, to accurately quantify the blood flow and structural mechanics of rat aorta, the motion of the aortic root, the compliance of aortic wall, and the pulsation of blood flow should be considered. However, when only focusing on the effective stress in rat aorta, the motion of the aortic root may be neglected.

Go with the flow: modeling unique biological flows in engineered in vitro platforms

Interest in recapitulating in vivo phenomena in vitro using organ-on-a-chip technology has grown rapidly and with it, attention to the types of fluid flow experienced in the body has followed suit. These platforms offer distinct advantages over in vivo models with regards to human relevance, cost, and control of inputs (e.g., controlled manipulation of biomechanical cues from fluid perfusion). Given the critical role biophysical forces play in several tissues and organs, it is therefore imperative that engineered in vitro platforms capture the complex, unique flow profiles experienced in the body that are intimately tied with organ function. In this review, we outline the complex and unique flow regimes experienced by three different organ systems: blood vasculature, lymphatic vasculature, and the intestinal system. We highlight current state-of-the-art platforms that strive to replicate physiological flows within engineered tissues while introducing potential limitations in current approaches.

Toward an accurate estimation of wall shear stress from 4D flow magnetic resonance downstream of a severe stenosis

PURPOSE

First, to investigate the agreement between velocity, velocity gradient, and Reynolds stress obtained from four-dimensional flow magnetic resonance (4D flow MRI) measurements and direct numerical simulation (DNS). Second, to propose and optimize based on DNS, 2 alternative methods for the accurate estimation of wall shear stress (WSS) when the resolution of the flow measurements is limited. Thirdly, to validate the 2 methods based on 4D flow MRI data.

METHODS

In vitro 4D MRI has been conducted in a realistic rigid stenosed aorta model under a constant flow rate of 12 L/min. A DNS of transitional stenotic flow has been performed using the same geometry and boundary conditions.

RESULTS

Time-averaged velocity and Reynolds stresses are in good agreement between in vitro 4D MRI data and DNS (errors between 2% and 8% of the reference downsampled data). WSS estimation based on the 2 proposed methods applied to MRI data provide good agreement with DNS for slice-averaged values (maximum error is less than 15% of the mean reference WSS for the first method and 25% for the second method). The performance of both models is not strongly sensitive to spatial resolution up to 1.5 mm voxel size. While the performance of model 1 deteriorates appreciably at low signal-to-noise ratios, model 2 remains robust.

CONCLUSIONS

The 2 methods for WSS magnitude give an overall better agreement than the standard approach used in the literature based on direct calculation of the velocity gradient close to the wall (relative error of 84%).

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.

The Story of Wall Shear Stress in Coronary Artery Atherosclerosis: Biochemical Transport and Mechanotransduction

Coronary artery atherosclerosis is a local, multifactorial, complex disease, and the leading cause of death in the US. Complex interactions between biochemical transport and biomechanical forces influence disease growth. Wall shear stress (WSS) affects coronary artery atherosclerosis by inducing endothelial cell mechanotransduction and by controlling the near-wall transport processes involved in atherosclerosis. Each of these processes is controlled by WSS differently and therefore has complicated the interpretation of WSS in atherosclerosis. In this paper, we present a comprehensive theory for WSS in atherosclerosis. First, a short review of shear stress-mediated mechanotransduction in atherosclerosis was presented. Next, subject-specific computational fluid dynamics (CFD) simulations were performed in ten coronary artery models of diseased and healthy subjects. Biochemical-specific mass transport models were developed to study low-density lipoprotein, nitric oxide, adenosine triphosphate, oxygen, monocyte chemoattractant protein-1, and monocyte transport. The transport results were compared with WSS vectors and WSS Lagrangian coherent structures (WSS LCS). High WSS magnitude protected against atherosclerosis by increasing the production or flux of atheroprotective biochemicals and decreasing the near-wall localization of atherogenic biochemicals. Low WSS magnitude promoted atherosclerosis by increasing atherogenic biochemical localization. Finally, the attracting WSS LCS's role was more complex where it promoted or prevented atherosclerosis based on different biochemicals. We present a summary of the different pathways by which WSS influences coronary artery atherosclerosis and compare different mechanotransduction and biotransport mechanisms.

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.

Possible Early Generation of Physiological Helical Flow Could Benefit the Triflo Trileaflet Heart Valve Prosthesis Compared to Bileaflet Valves

Background—Physiological helical flow in the ascending aorta has been well documented in the last two decades, accompanied by discussions on possible physiological benefits of such axial swirl. Recent 4D-MRI studies on healthy volunteers have found indications of early generation of helical flow, early in the systole and close to the valve plane. Objectives—Firstly, the aim of the study is to investigate the hypothesis of premature swirl existence in the ventricular outflow tract leading to helical flow in the valve plane, and second to investigate the possible impact of two different mechanical valve designs on the preservation of this early helical flow and its subsequent hemodynamic consequences. Methods—We use a pulse duplicator with an aortic arch and High-Speed Particle Image Velocimetry to document the flow evolution in the systolic cycle. The pulse-duplicator is modified with a swirl-generating insert to generate early helical flow in the valve plane. Special focus is paid to the interaction of such helical flow with different designs of mechanical prosthetic heart valves, comparing a classical bileaflet mechanical heart valve, the St. Jude Medical Regent valve (SJM Regent BMHV), with the Triflo trileaflet mechanical heart valve T2B version (Triflo TMHV). Results—When the swirl-generator is inserted, a vortex is generated in the core flow, demonstrating early helical flow in the valve plane, similar to the observations reported in the recent 4D-MRI study taken for comparison. For the Triflo trileaflet valve, the early helical flow is not obstructed in the central orifice, similar as in the case of the natural valve. Conservation of angular momentum leads to radial expansion of the core flow and flattening of the axial flow profile downstream in the arch. Furthermore, the early helical flow helps to overcome separation at the outer and inner curvature. In contrast, the two parallel leaflets for the bileaflet valve impose a flow straightener effect, annihilating the angular momentum, which has a negative impact on kinetic energy of the flow. Conclusion—The results imply better hemodynamics for the Triflo trileaflet valve based on hydrodynamic arguments under the discussed hypothesis. In addition, it makes the Triflo valve a better candidate for valve replacements in patients with a pathological generation of nonaxial velocity in the ventricle outflow tract.

Lipid-rich Plaques Detected by Near-infrared Spectroscopy Are More Frequently Exposed to High Shear Stress

High wall shear stress (WSS) and near-infrared spectroscopy (NIRS) detected lipid-rich plaque (LRP) are both known to be associated with plaque destabilization and future adverse cardiovascular events. However, knowledge of spatial co-localization of LRP and high WSS is lacking. This study investigated the co-localization of LRP based on NIRS and high WSS. Fifty-three patients presenting acute coronary syndrome underwent NIRS-intravascular-ultrasound (NIRS-IVUS) imaging of a non-culprit coronary artery. WSS was obtained using WSS profiling in 3D-reconstructions of the coronary arteries based on fusion of IVUS-segmented lumen and CT-derived 3D-centerline. Thirty-eight vessels were available for final analysis and divided into 0.5 mm/45° sectors. LRP sectors, as identified by NIRS, were more often colocalized with high WSS than sectors without LRP. Moreover, there was a dose-dependent relationship between lipid content and high WSS exposure. This study is a first step in understanding the evolution of LRPs to vulnerable plaques.

Does the inflow velocity profile influence physiologically relevant flow patterns in computational hemodynamic models of left anterior descending coronary artery?

Patient-specific computational fluid dynamics is a powerful tool for investigating the hemodynamic risk in coronary arteries. Proper setting of flow boundary conditions in computational hemodynamic models of coronary arteries is one of the sources of uncertainty weakening the findings of in silico experiments, in consequence of the challenging task of obtaining in vivo 3D flow measurements within the clinical framework. Accordingly, in this study we evaluated the influence of assumptions on inflow velocity profile shape on coronary artery hemodynamics. To do that, (1) ten left anterior descending coronary artery (LAD) geometries were reconstructed from clinical angiography, and (2) eleven velocity profiles with realistic 3D features such as eccentricity and differently shaped (single- and double-vortex) secondary flows were generated analytically and imposed as inflow boundary conditions. Wall shear stress and helicity-based descriptors obtained prescribing the commonly used parabolic velocity profile were compared with those obtained with the other velocity profiles. Our findings indicated that the imposition of idealized velocity profiles as inflow boundary condition is acceptable as long the results of the proximal vessel segment are not considered, in LAD coronary arteries. As a pragmatic rule of thumb, a conservative estimation of the length of influence of the shape of the inflow velocity profile on LAD local hemodynamics can be given by the theoretical entrance length for cylindrical conduits in laminar flow conditions.

Exploring wall shear stress spatiotemporal heterogeneity in coronary arteries combining correlation-based analysis and complex networks with computational hemodynamics

Atherosclerosis at the early stage in coronary arteries has been associated with low cycle-average wall shear stress magnitude. However, parallel to the identification of an established active role for low wall shear stress in the onset/progression of the atherosclerotic disease, a weak association between lesions localization and low/oscillatory wall shear stress has been observed. In the attempt to fully identify the wall shear stress phenotype triggering early atherosclerosis in coronary arteries, this exploratory study aims at enriching the characterization of wall shear stress emerging features combining correlation-based analysis and complex networks theory with computational hemodynamics. The final goal is the characterization of the spatiotemporal and topological heterogeneity of wall shear stress waveforms along the cardiac cycle. In detail, here time-histories of wall shear stress magnitude and wall shear stress projection along the main flow direction and orthogonal to it (a measure of wall shear stress multidirectionality) are analyzed in a representative dataset of 10 left anterior descending pig coronary artery computational hemodynamics models. Among the main findings, we report that the proposed analysis quantitatively demonstrates that the model-specific inlet flow-rate shapes wall shear stress time-histories. Moreover, it emerges that a combined effect of low wall shear stress magnitude and of the shape of the wall shear stress–based descriptors time-histories could trigger atherosclerosis at its earliest stage. The findings of this work suggest for new experiments to provide a clearer determination of the wall shear stress phenotype which is at the basis of the so-called arterial hemodynamic risk hypothesis in coronary arteries.

The role of oxygen transport in atherosclerosis and vascular disease

Atherosclerosis and vascular disease of larger arteries are often associated with hypoxia within the layers of the vascular wall. In this review, we begin with a brief overview of the molecular changes in vascular cells associated with hypoxia and then emphasize the transport mechanisms that bring oxygen to cells within the vascular wall. We focus on fluid mechanical factors that control oxygen transport from lumenal blood flow to the intima and inner media layers of the artery, and solid mechanical factors that influence oxygen transport to the adventitia and outer media via the wall's microvascular system-the vasa vasorum (VV). Many cardiovascular risk factors are associated with VV compression that reduces VV perfusion and oxygenation. Dysfunctional VV neovascularization in response to hypoxia contributes to plaque inflammation and growth. Disturbed blood flow in vascular bifurcations and curvatures leads to reduced oxygen transport from blood to the inner layers of the wall and contributes to the development of atherosclerotic plaques in these regions. Recent studies have shown that hypoxia-inducible factor-1α (HIF-1α), a critical transcription factor associated with hypoxia, is also activated in disturbed flow by a mechanism that is independent of hypoxia. A final section of the review emphasizes hypoxia in vascular stenting that is used to enlarge vessels occluded by plaques. Stenting can compress the VV leading to hypoxia and associated intimal hyperplasia. To enhance oxygen transport during stenting, new stent designs with helical centrelines have been developed to increase blood phase oxygen transport rates and reduce intimal hyperplasia. Further study of the mechanisms controlling hypoxia in the artery wall may contribute to the development of therapeutic strategies for vascular diseases.

Histopathologic and physiologic effect of bifurcation stenting: current status and future prospects

Introduction: Coronary bifurcation lesions are involved in up to 20% of all percutaneous coronary interventions (PCI). However, bifurcation lesion intervention is associated with a high complication rate, and optimal treatment of coronary bifurcation is an ongoing debate.Areas covered: Both different stenting techniques and a variety of devices have been suggested for bifurcation treatment, including the use of conventional coronary stents, bioresorbable vascular scaffolds (BVS), drug-eluting balloons (DEB), and stents dedicated to bifurcations. This review will summarize different therapeutic approaches with their advantages and shortcomings, with special emphasis on histopathologic and physiologic effects of each treatment strategy.Expert opinion: Histopathology and clinical data have shown that a more simple treatment strategy is beneficial in bifurcation lesions, achieving superior results. Bifurcation interventions through balloon angioplasty or placement of stents can importantly alter the bifurcation's geometry and accordingly modify local flow conditions. Computational fluid dynamics (CFD) studies have shown that the outcome of bifurcation interventions is governed by local hemodynamic shear conditions. Minimizing detrimental flow conditions as much as possible should be the ultimate strategy to achieve long-term success of bifurcation interventions.