Computational study for the effects of coil configuration on blood flow characteristics in coil-embolized cerebral aneurysm

A comprehensive numerical approach to coil placement in cerebral aneurysms: mathematical modeling and in silico occlusion classification

Endovascular coil embolization is one of the primary treatment techniques for cerebral aneurysms. Although it is a well-established and minimally invasive method, it bears the risk of suboptimal coil placement which can lead to incomplete occlusion of the aneurysm possibly causing recurrence. One of the key features of coils is that they have an imprinted natural shape supporting the fixation within the aneurysm. For the spatial discretization, our mathematical coil model is based on the discrete elastic rod model which results in a dimension-reduced 1D system of differential equations. We include bending and twisting responses to account for the coils natural curvature and allow for the placement of several coils having different material parameters. Collisions between coil segments and the aneurysm wall are handled by an efficient contact algorithm that relies on an octree based collision detection. In time, we use a standard symplectic semi-implicit Euler time stepping method. Our model can be easily incorporated into blood flow simulations of embolized aneurysms. In order to differentiate optimal from suboptimal placements, we employ a suitable in silico Raymond-Roy-type occlusion classification and measure the local packing density in the aneurysm at its neck, wall region and core. We investigate the impact of uncertainties in the coil parameters and embolization procedure. To this end, we vary the position and the angle of insertion of the micro-catheter, and approximate the local packing density distributions by evaluating sample statistics.

Influence of framing coil orientation and its shape on the hemodynamics of a basilar aneurysm model

Aneurysms are bulges of an artery, which require clinical management solutions. Due to the inherent advantages, endovascular coil filling is emerging as the treatment of choice for intracranial aneurysms (IAs). However, after successful treatment of coil embolization, there is a serious risk of recurrence. It is well known that optimal packing density will enhance treatment outcomes. The main objective of endovascular coil embolization is to achieve flow stasis by enabling significant reduction in intra-aneurysmal flow and facilitate thrombus formation. The present study numerically investigates the effect of framing coil orientation on intra-aneurysmal hemodynamics. For the purpose of analysis, actual shape of the embolic coil is used, instead of simplified ideal coil shape. Typically used details of the framing coil are resolved for the analysis. However, region above the framing coil is assumed to be filled with a porous medium. Present simulations have shown that orientation of the framing coil loop (FCL) greatly influences the intra-aneurysmal hemodynamics. The FCLs which were placed parallel to the outlets of basilar tip aneurysm (Coil A) were found to reduce intra-aneurysmal flow velocity that facilitates thrombus formation. Involving the coil for the region is modeled using a porous medium model with a packing density of 20 % . The simulations indicate that the framing coil loop (FCL) has a significant influence on the overall outcome.

Modeling the Mechanical Microenvironment of Coiled Cerebral Aneurysms

Successful occlusion of cerebral aneurysms using coil embolization is contingent upon stable thrombus formation, and the quality of the thrombus depends upon the biomechanical environment. The goal of this study was to investigate how coil embolization alters the mechanical micro-environment within the aneurysm dome. Inertialess particles were injected in three-dimensional, computational simulations of flow inside patient aneurysms using patient-specific boundary conditions. Coil embolization was simulated as a homogenous porous medium of known permeability and inertial constant. Lagrangian particle tracking was used to calculate the residence time and shear stress history for particles in the flow before and after treatment. The percentage of particles entering the aneurysm dome correlated with the neck surface area before and after treatment (pretreatment: R2 = 0.831, P < 0.001; post-treatment: R2 = 0.638, P < 0.001). There was an inverse relationship between the change in particles entering the dome and coil packing density (R2 = 0.600, P < 0.001). Following treatment, the particles with the longest residence times tended to remain within the dome even longer while accumulating lower shear stress. A significant correlation was observed between the treatment effect on residence time and the ratio of the neck surface area to porosity (R2 = 0.390, P = 0.007). The results of this study suggest that coil embolization triggers clot formation within the aneurysm dome via a low shear stress-mediated pathway. This hypothesis links independently observed findings from several benchtop and clinical studies, furthering our understanding of this treatment strategy.

Tailored Sac Embolization During EVAR for Preventing Persistent Type II Endoleak

BACKGROUND

Persistent type II endoleaks (ELIIp) occur in 8-23% of patients submitted to endovascular aneurysm repair (EVAR) and may lead to aneurysm progression and rupture. Intraoperative embolization of the abdominal aortic aneurysm (AAA) sac is effective to prevent their occurrence, however a method to achieve complete sac thrombosis has not been standardized yet. Aim of our study was to identify factors associated with prevention of ELIIp after intraoperative embolization, in order to optimize technical details.

METHODS

Patients at high risk for ELIIp, who underwent EVAR with AAA - sac coil embolization were prospectively collected into a dedicated database from January 2012 to March 2015. The endoluminal residual sac volume (ERV), not occupied by the endograft [ERV= AAA total volume (TV) - (AAA-thrombus volume (THV) + endograft volume (EgV)] was calculated on preoperative computed tomography and the concentration of coils implanted (CCoil= n coils implanted/ERV) for each patient was evaluated. AAA volumetric evaluation was conducted by dedicated vessels analysis software (3Mensio). ELIIp presence was evaluated by contrast-enhanced ultrasound at 6 and 12-month. Patients with ELIIp at 12 months (Group 1) were clustered and compared to patients without ELIIp (Group 2), in order to evaluate the incidence of ELIIp in patients undergone to preventive AAA-sac embolization, and identify the predictors of ELIIp prevention. Morphological potential risk factors for ELIIp such as TV, THV, VR% and EgV were also considered in all patients. Statistical correlation was assessed by Fisher Exact Test.

RESULTS

Among 326 patients undergone to standard EVAR, 61 (19% - M: 96.7%, median age: 72 [IQR: 8] years, median AAA diameter: 57 [IQR: 7] mm) were considered at high risk for ELIIp and were submitted to coil embolization. The median AAA total volume (TV) and median ERV were 156 (IQR: 59) cc and 46 (IQR: 26) cc, respectively. The median number and concentration of coils (IMWCE-38-16-45 Cook M-Ray) positioned in AAA-sac were 5 (IQR: 1) coils and 0.17 coil/cm³ (range 0.02-1.20). Among this high-risk population, the incidence of ELIIp was 29.5% and 23% at 6 and 12-month, respectively. Fourteen patients (23%) were clustered in Group1 and 47 (77%) in Group 2. Both groups were homogeneous for clinical characteristics and preoperative morphological risk factors. There were no differences in the preoperative median TV, AAA-thrombus volume (THV), %VR, EgV and number of implanted coils between Group1 and Group2. Patients in Group1 had a significantly higher ERV (59 [IQR: 13] cm³ vs. 42 [IQR: 27] cm³, P = 0.002) and lower CCoil (0.09 [IQR: 0.03] vs. 0.18 [IQR: 0.21], P = 0.006) than patients of Group2. ELIIp was significantly related to the presence of ERV > 49 cm³ (86 % vs. 42 %, Group1 and Group2 respectively, P = 0.006) and CCoil < 0.17coil/ cm³ (100% vs. 68%, Group1 e Group2 respectively, P = 0.014).

CONCLUSION

According with our results, Coil concentration and endoluminal residual volume can affect the efficacy of the AAA - sac embolization in the prevention of ELIIp, moreover CCoil ≥0.17coil/ cm³ maight be considered to determine the tailored number of coils.

Cerebral arterial architectonics and CFD simulation in mice with type 1 diabetes mellitus of different duration

Type 1 diabetes is a chronic autoimmune disease that affects tens of millions of people. Diabetes mellitus is one of the strongest factors in the development of cerebrovascular diseases. In this study we used NOD.CB17 Prkdcscid mice and the pharmacological model of type 1 diabetes mellitus of different duration to study changes in the cerebral vasculature. We used two combined approaches using magnetic resonance angiography both steady and transient CFD blood flow modeling. We identified the influence of type 1 diabetes on the architectonics and hemodynamics of the large blood vessels of the brain as the disease progresses. For the first time, we detected a statistically significant change in angioarchitectonics (the angles between the vessels of the circle of Willis, cross-sections areas of vessels) and hemodynamic (maximum blood flow rate, hydraulic resistance) in animals with diabetes duration of 2 months, that is manifested by the development of asymmetry of cerebral blood flow. The result shows the negative effect of diabetes on cerebral circulation as well as the practicability of CFD modeling. This may be of extensive interest, in pharmacological and preclinical studies.

Effect of Local Coil Density on Blood Flow Stagnation in Densely Coiled Cerebral Aneurysms: A Computational Study Using a Cartesian Grid Method

Aneurysm recurrence is the most critical concern following coil embolization of a cerebral aneurysm. Adequate packing density (PD) and coil uniformity are believed necessary to achieve sufficient flow stagnation, which decreases the risk of aneurysm recurrence. The effect of coil distribution on the extent of flow stagnation, however, especially in cases of dense packing (high PD), has received less attention. Thus, the cause of aneurysm recurrence despite dense packing is still an open question. The primary aim of this study is to evaluate the effect of local coil density on the extent of blood flow stagnation in densely coiled aneurysms. For this purpose, we developed a robust computational framework to determine blood flow using a Cartesian grid method, by which the complex fluid pathways in coiled aneurysms could be flexibly treated using an implicit function. This tool allowed us to conduct blood flow analyses in two patient-specific geometries with 50 coil distribution patterns in each aneurysm at clinically adequate PD. The results demonstrated that dense packing in the aneurysm may not necessarily block completely the inflow into the aneurysm and local flow that formed in the neck region, whose strength was inversely related to this local PD. This finding suggests that local coil density in the neck region still plays an important role in disturbing the remaining local flow, which possibly prevents thrombus formation in a whole aneurysm sac, increasing the risk of aneurysm regrowth and subsequent recurrence.

In vitro measurement of platelet adhesion to intact endothelial cells under low shear conditions

BACKGROUND

Prediction of thrombus formation at intact arterial walls under low shear flow conditions is clinically important particularly for better prognoses of embolisation in cerebral aneurysms. Although a new mathematical model for this purpose is necessary, little quantitative information has been known about platelet adhesion to intact endothelial cells.

OBJECTIVE

The objective of this study is to measure the number of platelets adhering to intact endothelial cells with a focus upon the influence of the shear rate.

METHODS

Endothelial cells disseminated in μ-slides were exposed to swine whole blood at different shear rates. Adenosine diphosphate (ADP) was used as an agonist. Adherent platelets were counted by means of scanning electron microscopy.

RESULTS

At an ADP concentration of 1 µM, 20.8 ± 3.1 platelets per 900 µm2 were observed after 30-minute perfusion at a shear rate of 0.8 s-1 whereas only 3.0 ± 1.4 per 900 µm2 at 16.8 s-1.

CONCLUSIONS

The number of adherent platelets is determined by a balance between the shear and the degree of stimulation by the agonist. At an ADP concentration of 1 µM, a limit to the shear rate at which platelets can adhere to intact endothelial cells is considered to be slightly higher than 16.8 s-1.

In vitro measurement of the permeability of endovascular coils deployed in cerebral aneurysms

BACKGROUND AND PURPOSE

Aneurysm recurrence is the primary limitation of endovascular coiling treatment for cerebral aneurysms. Coiling is currently quantified by a volumetric porosity measure called packing density (pd). Blood flow through a coil mass depends on the permeability of the coil mass, and not just its pd. The permeability of coil masses has not yet been quantified. Here we measure coil permeability with a traditional falling-head permeameter modified to incorporate idealized aneurysms.

METHODS

Silicone replicas of idealized aneurysms were manufactured with three different aneurysm diameters (4, 5, and 8 mm). Four different coil types (Codman Trufill Orbit, Covidien Axium, Microvention Microplex 10, and Penumbra 400) were deployed into the aneurysms with a target pd of 35%. Coiled replicas were installed on a falling-head permeameter setup and the time taken for a column of fluid above the aneurysm to drop a certain height was recorded. Permeability of the samples was calculated based on a simple modification of the traditional permeameter equation to incorporate a spherical aneurysm.

RESULTS

The targeted 35% pd was achieved for all samples (35%±1%, P=0.91). Coil permeabilities were significantly different from each other (P<0.001) at constant pd. Microplex 10 coils had the lowest permeability of all coil types. Data suggest a trend of increasing permeability with thicker coil wire diameter (not statistically significant).

CONCLUSIONS

A simple in vitro setup was developed to measure the permeabilities of coil masses based on traditional permeametry. Coil permeability should be considered when evaluating the hemodynamic efficacy of coiling instead of just packing density. Coils made of thicker wires may be more permeable, and thus less effective, than coils made from thinner wires. Whether aneurysm recurrence is affected by coil wire diameter or permeability needs to be confirmed with clinical trials.