Wall shear stress gradient is independently associated with middle cerebral artery aneurysm development: a case-control CFD patient-specific study based on 77 patients

Hemodynamic factors of spontaneous vertebral artery dissecting aneurysms assessed with numerical and deep learning algorithms: Role of blood pressure and asymmetry

BACKGROUND AND OBJECTIVES

The pathophysiology of spontaneous vertebral artery dissecting aneurysms (SVADA) is poorly understood. Our goal is to investigate the hemodynamic factors contributing to their formation using computational fluid dynamics (CFD) and deep learning algorithms.

METHODS

We have developed software that can use patient imagery as input to recreate the vertebrobasilar arterial system, both with and without SVADA, which we used in a series of three patients. To obtain the kinematic blood flow data before and after the aneurysm forms, we utilized numerical methods to solve the complex Navier-Stokes partial differential equations. This was accomplished through the application of a finite volume solver (OpenFoam/Helyx OS). Additionally, we trained a neural ordinary differential equation (NODE) to learn and replicate the dynamical streamlines obtained from the computational fluid dynamics (CFD) simulations.

RESULTS

In all three cases, we observed that the equilibrium of blood pressure distributions across the VAs, at a specific vertical level, accurately predicted the future SVADA location. In the two cases where there was a dominant VA, the dissection occurred on the dominant artery where blood pressure was lower compared to the contralateral side. The SVADA sac was characterized by reduced wall shear stress (WSS) and decreased velocity magnitude related to increased turbulence. The presence of a high WSS gradient at the boundary of the SVADA may explain its extension. Streamlines generated by CFD were learned with a neural ordinary differential equation (NODE) capable of capturing the system's dynamics to output meaningful predictions of the flow vector field upon aneurysm formation.

CONCLUSION

In our series, asymmetry in the vertebrobasilar blood pressure distributions at and proximal to the site of the future SVADA accurately predicted its location in all patients. Deep learning algorithms can be trained to model blood flow patterns within biological systems, offering an alternative to the computationally intensive CFD. This technology has the potential to find practical applications in clinical settings.

Constrained estimation of intracranial aneurysm surface deformation using 4D-CTA

BACKGROUND AND OBJECTIVE

Intracranial aneurysms are relatively common life-threatening diseases, and assessing aneurysm rupture risk and identifying the associated risk factors is essential. Parameters such as the Oscillatory Shear Index, Pressure Loss Coefficient, and Wall Shear Stress are reliable indicators of intracranial aneurysm development and rupture risk, but aneurysm surface irregular pulsation has also received attention in aneurysm rupture risk assessment.

METHODS

The present paper proposed a new approach to estimate aneurysm surface deformation. This method transforms the estimation of aneurysm surface deformation into a constrained optimization problem, which minimizes the error between the displacement estimated by the model and the sparse data point displacements from the four-dimensional CT angiography (4D-CTA) imaging data.

RESULTS

The effect of the number of sparse data points on the results has been discussed in both simulation and experimental results, and it shows that the proposed method can accurately estimate the surface deformation of intracranial aneurysms when using sufficient sparse data points.

CONCLUSIONS

Due to a potential association between aneurysm rupture and surface irregular pulsation, the estimation of aneurysm surface deformation is needed. This paper proposed a method based on 4D-CTA imaging data, offering a novel solution for the estimation of intracranial aneurysm surface deformation.

Investigation of paraclinoid aneurysm formation by comparing the combined influence of hemodynamic parameters between aneurysmal and non-aneurysmal arteries

Numerous studies have evaluated the effects of hemodynamic parameters on aneurysm formation. However, the reasons why aneurysms do not initiate in intracranial arteries are still unclear. This study aimed to investigate the influence of hemodynamic parameters, wall shear stress (WSS) and strain, on aneurysm formation by comparing between aneurysmal and non-aneurysmal arteries. Fifty-eight patients with paraclinoid aneurysms on one side were enrolled. Based on magnetic resonance angiography, each patient's left and right internal carotid arteries (ICAs) were reconstructed. For a patient having an aneurysm on one side, the ICA with the paraclinoid aneurysm was defined as the aneurysmal artery after eliminating the aneurysm, whereas the opposite ICA without aneurysm was defined as the non-aneurysmal artery. Computational fluid dynamics and fluid-structure interaction analyses were then performed for both aneurysmal and non-aneurysmal arteries. Finally, the relationship between high hemodynamic parameters and aneurysm location was investigated. For aneurysmal arteries, high WSS and strain locations were well-matched with the aneurysm formation site. Also, considerable correlations between high WSS and strain locations were observed. However, there was no significant relationship between high hemodynamic parameters and aneurysm formation for non-aneurysmal arteries. The findings are helpful for understanding aneurysm formation mechanism and encouraging further relevant research.

Association of Arterial Tortuosity with Hemodynamic Parameters-A Computational Fluid Dynamics Study

BACKGROUND

Tortuosity of intracranial arteries has been proven to be associated with the risk of intracranial aneurysm development. We decided to analyze which aspects of tortuosity are correlated with hemodynamics parameters promoting intracranial aneurysm development.

METHODS

We constructed 73 idealized models of tortuous artery (length: 25 mm, diameter: 2.5 mm) with single bifurcation. For each model, on the course of segment before bifurcation, we placed 1-3 angles with measures 15, 30, 45, 60, or 75 degrees and arc lengths 2, 5, 7, 10, or 15 mm. We performed computational fluid dynamics analysis. Blood was modeled as Newtonian fluid. We have set velocity wave of 2 cardiac cycles. After performing simulation we calculated following hemodynamic parameters at the bifurcation: time average wall shear stress (TAWSS), time average wall shear stress gradient (TAWSSG), oscillatory shear index (OSI), and relative residence time (RRT).

RESULTS

We found a significant positive correlation with number of angles and TAWSS (R = 0.329; P < 0.01), TAWSSG (R = 0.317; P < 0.01), and negative with RRT (R = -0.335; P < 0.0.01). Similar results were obtained in terms of arcs lengths. On the other hand, mean angle measure was negatively correlated to TAWSS (R = -0.333; P < 0.01), TAWSSG (R = -0.473 P < 0.01), OSI (R = -0.463; P < 0.01), and positively to RRT (R = 0.332; P < 0.01). On the basis of the obtained results, we developed new tortuosity descriptor, which considered angle measures normalized to its arc length and distance from bifurcation. For such descriptor we found strong negative correlation with TAWSS (R = -0.701; P < 0.01), TAWSSG (R = 0.778; P < 0.01), OSI (R = -0.776; P < 0.01), and positive with RRT (R = 0.747; P < 0.01).

CONCLUSIONS

Hemodynamic parameters promoting aneurysm development are correlated with larger number of smaller angles located on larger arcs.

3D bioprinted multilayered cerebrovascular conduits to study cancer extravasation mechanism related with vascular geometry

Cerebral vessels are composed of highly complex structures that facilitate blood perfusion necessary for meeting the high energy demands of the brain. Their geometrical complexities alter the biophysical behavior of circulating tumor cells in the brain, thereby influencing brain metastasis. However, recapitulation of the native cerebrovascular microenvironment that shows continuities between vascular geometry and metastatic cancer development has not been accomplished. Here, we apply an in-bath 3D triaxial bioprinting technique and a brain-specific hybrid bioink containing an ionically crosslinkable hydrogel to generate a mature three-layered cerebrovascular conduit with varying curvatures to investigate the physical and molecular mechanisms of cancer extravasation in vitro. We show that more tumor cells adhere at larger vascular curvature regions, suggesting that prolongation of tumor residence time under low velocity and wall shear stress accelerates the molecular signatures of metastatic potential, including endothelial barrier disruption, epithelial-mesenchymal transition, inflammatory response, and tumorigenesis. These findings provide insights into the underlying mechanisms driving brain metastases and facilitate future advances in pharmaceutical and medical research.

Characteristics and evaluation of atherosclerotic plaques: an overview of state-of-the-art techniques

Atherosclerosis is an important cause of cerebrovascular and cardiovascular disease (CVD). Lipid infiltration, inflammation, and altered vascular stress are the critical mechanisms that cause atherosclerotic plaque formation. The hallmarks of the progression of atherosclerosis include plaque ulceration, rupture, neovascularization, and intraplaque hemorrhage, all of which are closely associated with the occurrence of CVD. Assessing the severity of atherosclerosis and plaque vulnerability is crucial for the prevention and treatment of CVD. Integrating imaging techniques for evaluating the characteristics of atherosclerotic plaques with computer simulations yields insights into plaque inflammation levels, spatial morphology, and intravascular stress distribution, resulting in a more realistic and accurate estimation of plaque state. Here, we review the characteristics and advancing techniques used to analyze intracranial and extracranial atherosclerotic plaques to provide a comprehensive understanding of atheroma.

Studying the imaging features and infarction mechanism of vertebrobasilar dolichoectasia with high-resolution magnetic resonance imaging

The mechanisms underlying ischemic infarction in patients with vertebrobasilar dolichoectasia (VBD) remain unclear. In this study, we retrospectively analyzed the imaging characteristics of high-resolution magnetic resonance imaging (HR-MRI) in patients with VBD to explore the possible mechanisms of ischemic stroke (IS) in patients with VBD. Patients with VBD were recruited from the HR-MRI database between July 2017 and June 2021. HR-MRI was used to evaluate the diameter, bifurcation height, laterality, arterial dissection, and atherosclerotic plaques of the basilar artery (BA). Transcranial Doppler was used to measure the vertebrobasilar mean velocity (Vm), peak systolic velocity (Vs), end-diastolic velocity (Vd), and pulsatile index. Twenty-six patients with VBD were enrolled, of which 15 had IS and 11 did not. The incidence of classical vascular risk factors, including age, sex, hypertension, diabetes, and hypercholesterolemia, did not differ significantly between the two groups. The BA diameters of the stroke group were significantly higher than that of the nonstroke group (6.57 ± 1.00 mm vs. 5.06 ± 0.50 mm, p = 0.000). The height of the BA bifurcation in the stroke and nonstroke groups was statistically significant (p = 0.002). Compared with the nonstroke group, the Vm, Vs, and Vd of the BA in the stroke group were lower, but the difference was not significant. In the 16 patients with atherosclerotic stenosis, 30 atherosclerotic plaques were found in the BA, 18 (60%) in the greater curvature, and 12 (40%) in the lesser curvature. In addition, one artery dissection (on the lesser curvature) and two dissecting aneurysms (on the greater curvature) were found in the BA of three patients, respectively. The BA diameter and bifurcation height are factors related to IS in patients with VBD. The mechanism of stroke in patients with VBD may involve abnormal hemodynamics, artery dissection, and atherosclerosis. HR-MRI is a useful method for evaluating the risk and underlying mechanism of stroke in patients with VBD.

Advanced cross-sectional imaging of cerebral aneurysms

While the rupture rate of cerebral aneurysms is only 1% per year, ruptured aneurysms are associated with significant morbidity and mortality, while aneurysm treatments have their own associated risk of morbidity and mortality. Conventional markers for aneurysm rupture include patient-specific and aneurysm-specific characteristics, with the development of scoring systems to better assess rupture risk. These scores, however, rely heavily on aneurysm size, and their accuracy in assessing risk in smaller aneurysms is limited. While the individual risk of rupture of small aneurysms is low, due to their sheer number, the largest proportion of ruptured aneurysms are small aneurysms. Conventional imaging techniques are valuable in characterizing aneurysm morphology; however, advanced imaging techniques assessing the presence of inflammatory changes within the aneurysm wall, hemodynamic characteristics of blood flow within aneurysm sacs, and imaging visualization of irregular aneurysm wall motion have been used to further determine aneurysm instability that otherwise cannot be characterized by conventional imaging techniques. The current manuscript reviews conventional imaging techniques and their value and limitations in cerebral aneurysm characterization, and evaluates the applications, value and limitations of advanced aneurysm imaging and post-processing techniques including intracranial vessel wall MRA, 4D-flow, 4D-CTA, and computational fluid dynamic simulations.

A predictive hemodynamic model based on risk factors for ruptured mirror aneurysms

Objectives

To identify hemodynamic risk factors for intracranial aneurysm rupture and establish a predictive model to aid evaluation.

Methods

We analyzed the hemodynamic parameters of 91 pairs of ruptured mirror aneurysms. A conditional univariate analysis was used for the continuous variables. A conditional multivariate logistic regression analysis was performed to identify the independent risk factors. Differences where p < 0.05 were statistically significant. A predictive model was established based on independent risk factors. Odds ratios (ORs) were used to score points. The validation cohort consisted of 189 aneurysms. Receiver operating characteristic curves were generated to determine the cutoff values and area under the curves (AUCs) of the predictive model and independent risk factors.

Results

The conditional multivariate logistic analysis showed that the low shear area (LSA) (OR = 70.322, p = 0.044, CI = 1.112–4,445.256), mean combined hemodynamic parameter (CHP) (>0.087) (OR = 3.171, p = 0.034, CI = 1.089–9.236), and wall shear stress gradient (WSSG) ratio (>893.180) (OR = 5.740, p = 0.003, CI = 1.950–16.898) were independent risk factors. A prediction model was established: 23^*LSA + 1^*CHP mean (>0.087: yes = 1, no = 0) + 2 ^* WSSG ratio (>893.180: yes = 1, no = 0). The AUC values of the predictive model, LSA, mean CHP (>0.087), and WSSG ratio (>893.180) were 0.748, 0.700, 0.654, and 0.703, respectively. The predictive model and LSA cutoff values were 1.283 and 0.016, respectively. In the validation cohort, the predictive model, LSA, CHP (>0.087), and WSSG ratio (>893.180) were 0.736, 0.702, 0.689, and 0.706, respectively.

Conclusions

LSA, CHP (>0.087), and WSSG ratio (>893.180) were independent risk factors for aneurysm rupture. Our predictive model could aid practical evaluation.

Impinging Flow Induces Expression of Monocyte Chemoattractant Protein-1 in Endothelial Cells Through Activation of the c-Jun N-terminal Kinase/c-Jun/p38/c-Fos Pathway

OBJECTIVE

Monocyte chemoattractant protein-1 (MCP-1) is an important regulator of the formation and development of intracranial aneurysms. This study explored the molecular mechanisms underlying the induction of MCP-1 and related inflammatory factors in human umbilical vein endothelial cells (HUVECs) under hemodynamic conditions.

METHODS

A modified T chamber was used to simulate fluid flow at the bifurcation of the artery and wall shear stress on HUVECs in vitro. Changes in HUVECs were analyzed in response to impinging flow. And HUVECs without impinging flow were used as the control group. Protein expression levels of extracellular signal-regulated kinase (ERK), c-Jun N-terminal kinase (JNK), p38, activator protein-1, and MCP-1 were detected by Western blot, and the messenger RNA expression levels of MCP-1, interleukin (IL)-1β, and IL-6 were determined by quantitative reverse transcription polymerase chain reaction.

RESULTS

Under impinging flow, the phosphorylation levels of ERK, JNK, and p38, as well as the protein levels of MCP-1, c-Jun, and c-Fos, increased. The messenger RNA expression of MCP-1, IL-1β, and IL-6 also increased in HUVECs. Pretreatment of the HUVECs with inhibitors of JNK and p38 significantly attenuated the increased expression of MCP-1, IL-1β, and IL-6, while ERK inhibitors had no obvious effect.

CONCLUSIONS

Under impinging flow, MCP-1 and inflammatory factors are regulated through the JNK/c-Jun/p38/c-Fos pathway and participate in EC inflammation.

Medial Gap: A Structural Factor at the Arterial Bifurcation Aggravating Hemodynamic Insult

Previous studies have reported that intracranial aneurysms frequently occur adjacent to the medial gap. However, the role of the medial gap in aneurysm formation is controversial. We designed this study to explore the potential role of the medial gap in aneurysm formation. Widened artery bifurcations with or without medial gaps were microsurgically created and pathologically stained in the carotid arteries of 30 rats. Numerical artery bifurcation models were constructed, and bidirectional fluid-solid interaction analyses were performed. Animal experiments showed that the apexes of widened bifurcations with a medial gap were prone to being insulted by blood flow compared to those without a medial gap. The bidirectional fluid-solid interaction analyses indicated that artery bifurcations with the medial gap exhibited higher wall shear stress (WSS) and von Mises stress (VMS) at the apex of the bifurcation. The disparity of stress between the gap and no-gap model was larger for widened bifurcations, peaking at 180° with a maximum of 1.9 folds. The maximum VMS and relatively high WSS were located at the junction between the medial gap and the adjacent arterial wall. Our results suggest that the medial gap at the widened arterial bifurcation may promote aneurysm formation.

Morphometry of cerebral arterial bifurcations harbouring aneurysms: a case-control study

Background

Conclusions from studies evaluating vessel dimensions and their deviations from values resulting from the principle of minimum work (PMW) on the formation of intracranial aneurysms (IAs) are still inconclusive. Our study aimed to perform a morphometric analysis of cerebral arterial bifurcations harbouring aneurysms.

Methods

The study comprised 147 patients with basilar artery (BA) and middle cerebral artery (MCA) aneurysms and 106 patients constituting the control group. The following morphometric parameters were evaluated: the radii of vessels forming the bifurcation, the junction exponent, the values of the bifurcation angles (Φ₁ and Φ₂ angles between the parent vessel trunk axis and the larger or smaller branches, respectively; α angle, the total bifurcation angle) and the difference between the predicted optimal and observed branch angles.

Results

The analysed parameters for internal carotid artery (ICA) bifurcations were not significantly different among the groups. The MCA and BA bifurcation angles and the radii of the parent MCA and BA vessels with aneurysms were significantly higher than those of the control group. The differences between the predicted optimal and observed branch angles were significantly higher for BA and MCA bifurcations with aneurysms compared to the control group. The mean junction exponent for bifurcations in the circle of Willis (i.e., ICA and BA bifurcations, respectively) and MCA bifurcations with aneurysms was significantly lower than the theoretical optimum and did not significantly differ among the groups. In a multilevel multivariate logistic regression analysis, the branch angles and the radius from the parent vessel were significant independent predictors of the presence of an IA. The ROC analysis indicated that the α angle was the best performer in discriminating between aneurysmal and nonaneurysmal bifurcations.

Conclusions

The dimensions of the arteries forming the circle of Willis do not follow the PMW. Deviation from the energetically optimum geometry for bifurcations beyond the circle of Willis (particularly, a larger radius of the parent artery and a wider total bifurcation angle) may lead to the formation of IAs. Further studies are warranted to investigate the significance of vessel dimensions and the bifurcation angle on the magnitude of shear stress in the walls of arterial bifurcations.