Wall stress of the cervical carotid artery in patients with carotid dissection: a case-control study

Spontaneous internal carotid artery (ICA) dissection (sICAD) results from an intimal tear located around the distal carotid sinus. The mechanisms causing the tear are unknown. This case-control study tested the hypotheses that head movements increase the wall stress in the cervical ICA and that the stress increase is greater in patients with sICAD than in controls. Five patients with unilateral, recanalized, left sICAD and five matched controls were investigated before and after maximal head rotation to the left and neck hyperextension after 45° head rotation to the left. The anatomy of the extracranial carotid arteries was assessed by magnetic resonance imaging and used to create finite element models of the right ICA. Wall stress increased after head movements. Increases above the 80th and 90th percentile were located at the intimal side of the artery wall from 7.4 mm below to 10 mm above the cranial edge of the carotid sinus, i.e., at the same location as histologically confirmed tears in patients with sICAD. Wall stress increase did not differ between patients and controls. The present findings suggest that wall stress increases at the intimal side of the artery wall surrounding the distal edge of the carotid bulb after head movements may be important for the development of carotid dissection. The lack of wall stress difference between the two groups indicates that the carotid arteries of patients with carotid dissection have either distinct functional or anatomical properties or endured unusually heavy wall stresses to initiate dissection.

A mathematical method for the 3D analysis of rotating deformable systems applied on lumen-forming MDCK cell aggregates

Cell motility contributes to the formation of organs and tissues, into which multiple cells self-organize. However such mammalian cellular motilities are not characterized in a quantitative manner and the systemic consequences are thus unknown. A mathematical tool to decipher cell motility, accounting for changes in cell shape, within a three-dimensional (3D) cell system was missing. We report here such a tool, usable on segmented images reporting the outline of clusters (cells) and allowing the time-resolved 3D analysis of circular motility of these as parts of a system (cell aggregate). Our method can analyze circular motility in sub-cellular, cellular, multi-cellular, and also non-cellular systems for which time-resolved segmented cluster outlines are available. To exemplify, we characterized the circular motility of lumen-initiating MDCK cell aggregates, embedded in extracellular matrix. We show that the organization of the major surrounding matrix fibers was not significantly affected during this cohort rotation. Using our developed tool, we discovered two classes of circular motion, rotation and random walk, organized in three behavior patterns during lumen initiation. As rotational movements were more rapid than random walk and as both could continue during lumen initiation, we conclude that neither the class nor the rate of motion regulates lumen initiation. We thus reveal a high degree of plasticity during a developmentally critical cell polarization step, indicating that lumen initiation is a robust process. However, motility rates decreased with increasing cell number, previously shown to correlate with epithelial polarization, suggesting that migratory polarization is converted into epithelial polarization during aggregate development.

Cerebrospinal fluid dynamics in the human cranial subarachnoid space: an overlooked mediator of cerebral disease. II. In vitro arachnoid outflow model

The arachnoid membrane (AM) and granulations (AGs) are important in cerebrospinal fluid (CSF) homeostasis, regulating intracranial pressure in health and disease. We offer a functional perspective of the human AM's transport mechanism to clarify the role of AM in the movement of CSF and metabolites. Using cultures of human AG cells and a specialized perfusion system, we have shown that this in vitro model mimics the in vivo characteristics of unidirectional fluid transport and we present the first report of serum-free permeability values (92.5 microl min(-1) mm Hg(-1) cm(-2)), which in turn are in agreement with the CSF outflow rates derived from a dynamic, in vivo magnetic resonance imaging-based computational model of the subarachnoid cranial space (130.9 microl min(-1) mm Hg(-1) cm(-2)). Lucifer yellow permeability experiments have verified the maintenance of tight junctions by the arachnoidal cells with a peak occurring around 21 days post-seeding, which is when all perfusion experiments were conducted. Addition of ruthenium red to the perfusate, and subsequent analysis of its distribution post-perfusion, has verified the passage of perfusate via both paracellular and transcellular mechanisms with intracellular vacuoles of approximately 1 microm in diameter being the predominant transport mechanism. The comparison of the computational and in vitro models is the first report to measure human CSF dynamics functionally and structurally, enabling the development of innovative approaches to modify CSF outflow and will change concepts and management of neurodegenerative diseases resulting from CSF stagnation.

Cerebrospinal fluid dynamics in the human cranial subarachnoid space: an overlooked mediator of cerebral disease. I. Computational model

Abnormal cerebrospinal fluid (CSF) flow is suspected to be a contributor to the pathogenesis of neurodegenerative diseases such as Alzheimer's through the accumulation of toxic metabolites, and to the malfunction of intracranial pressure regulation, possibly through disruption of neuroendocrine communication. For the understanding of transport processes involved in either, knowledge of in vivo CSF dynamics is important. We present a three-dimensional, transient, subject-specific computational analysis of CSF flow in the human cranial subarachnoid space (SAS) based on in vivo magnetic resonance imaging. We observed large variations in the spatial distribution of flow velocities with a temporal peak of 5 cm s(-1) in the anterior SAS and less than 4 mm s(-1) in the superior part. This could reflect dissimilar flushing requirements of brain areas that may show differences in susceptibility to pathological CSF flow. Our methods can be used to compare the transport of metabolites and neuroendocrine substances in healthy and diseased brains.

Long-term follow-up, computed tomography, and computational fluid dynamics of the Cabrol procedure

OBJECTIVES

The Cabrol procedure is characterized by insertion of an ascending aortic composite graft with reimplantation of the coronary arteries by the interposition of a graft tube. Our purpose is to report the clinical long-term follow-up and computed tomographic findings in patients having undergone the Cabrol procedure and to determine blood flow in the Cabrol graft using computational fluid dynamics.

METHODS

Clinical follow-up (76.6 +/- 16.6 months) and dual-source computed tomographic angiography data of 7 patients (all men, mean age 54.9 +/- 9.6 years) with 12 Cabrol grafts (left main coronary artery, n = 7; right coronary artery, n = 5) were reviewed. In 2 patients, the right coronary artery was directly reattached to the aortic graft. Computational fluid dynamics were calculated using computed tomographic data of a patient with the Cabrol procedure and compared with those in a Valsalva graft and a healthy aortic root.

RESULTS

Computed tomography showed Cabrol graft occlusions to 1 of 7 (14%) left main and of 2 of 5 (40%) right coronary arteries. Six grafts to the left main and 3 to the right coronary artery were fully patent, similar to the 2 directly reattached right coronary arteries to the aortic graft. Computational fluid dynamics results show similar blood flow parameters into the coronaries for the healthy aortic root and Valsalva graft. In the Cabrol graft, a spiraling flow pattern with low flow into the right coronary artery was found (right coronary artery = 1 mL/min at both systole and diastole).

CONCLUSIONS

Our study indicates low flow rates particularly in the right Cabrol graft correlating with a higher incidence of occlusions of the right as compared with the left Cabrol graft at long-term follow-up.

Three-dimensional computational modeling of subject-specific cerebrospinal fluid flow in the subarachnoid space

This study aims at investigating three-dimensional subject-specific cerebrospinal fluid (CSF) dynamics in the inferior cranial space, the superior spinal subarachnoid space (SAS), and the fourth cerebral ventricle using a combination of a finite-volume computational fluid dynamics (CFD) approach and magnetic resonance imaging (MRI) experiments. An anatomically accurate 3D model of the entire SAS of a healthy volunteer was reconstructed from high resolution T2 weighted MRI data. Subject-specific pulsatile velocity boundary conditions were imposed at planes in the pontine cistern, cerebellomedullary cistern, and in the spinal subarachnoid space. Velocimetric MRI was used to measure the velocity field at these boundaries. A constant pressure boundary condition was imposed at the interface between the aqueduct of Sylvius and the fourth ventricle. The morphology of the SAS with its complex trabecula structures was taken into account through a novel porous media model with anisotropic permeability. The governing equations were solved using finite-volume CFD. We observed a total pressure variation from -42 Pa to 40 Pa within one cardiac cycle in the investigated domain. Maximum CSF velocities of about 15 cms occurred in the inferior section of the aqueduct, 14 cms in the left foramen of Luschka, and 9 cms in the foramen of Magendie. Flow velocities in the right foramen of Luschka were found to be significantly lower than in the left, indicating three-dimensional brain asymmetries. The flow in the cerebellomedullary cistern was found to be relatively diffusive with a peak Reynolds number (Re)=72, while the flow in the pontine cistern was primarily convective with a peak Re=386. The net volumetric flow rate in the spinal canal was found to be negligible despite CSF oscillation with substantial amplitude with a maximum volumetric flow rate of 109 mlmin. The observed transient flow patterns indicate a compliant behavior of the cranial subarachnoid space. Still, the estimated deformations were small owing to the large parenchymal surface. We have integrated anatomic and velocimetric MRI data with computational fluid dynamics incorporating the porous SAS morphology for the subject-specific reconstruction of cerebrospinal fluid flow in the subarachnoid space. This model can be used as a basis for the development of computational tools, e.g., for the optimization of intrathecal drug delivery and computer-aided evaluation of cerebral pathologies such as syrinx development in syringomelia.

Patient-specific three-dimensional simulation of LDL accumulation in a human left coronary artery in its healthy and atherosclerotic states

We calculate low-density lipoprotein (LDL) transport from blood into arterial walls in a three-dimensional, patient-specific model of a human left coronary artery. The in vivo anatomy data are obtained from computed tomography images of a patient with coronary artery disease. Models of the artery anatomy in its healthy and diseased states are derived after segmentation of the vessel lumen, with and without the detected plaque, respectively. Spatial shear stress distribution at the endothelium is determined through the reconstruction of the arterial blood flow field using computational fluid dynamics. The arterial endothelium is represented by a shear stress-dependent, three-pore model, taking into account blood plasma and LDL passage through normal junctions, leaky junctions, and the vesicular pathway. Intraluminal pressures of 70 and 120 mmHg are employed as the normal and hypertensive operating pressures, respectively. By applying our model to both the healthy and diseased states, we show that the location of the plaque in the diseased state corresponds to one of the two sites with predicted high-LDL concentration in the healthy state. We further show that, in the diseased state, the site with high-LDL concentration has shifted distal to the plaque, which is in agreement with the clinical observation that plaques generally grow in the downstream direction. We also demonstrate that hypertension leads to increased number of regions with high-LDL concentration, elucidating one of the ways in which hypertension may promote atherosclerosis.

Identification of atherosclerotic lesion-prone sites through patient-specific simulation of low-density lipoprotein accumulation

We present a patient-specific model of low-density lipoprotein (LDL) transport from blood into arterial walls. To this end, the arterial endothelium is represented by a shear-stress dependent three-pore model taking into account blood plasma and LDL passage through the vesicular pathway, normal junctions and leaky junctions. We virtually remove atherosclerotic plaque from an in-vivo left coronary artery computed tomography (CT) dataset to obtain an approximation of the artery anatomy in its healthy state. By applying our model, we show that the location of the plaque in the diseased state corresponds to one of the two sites with predicted high LDL concentration in the healthy state. We further show that in the diseased state, the site with high LDL concentration has shifted distally, which is in agreement with the clinical observation that plaques generally grow in downstream direction.

Mixing and modes of mass transfer in the third cerebral ventricle: a computational analysis

Anatomic, velocimetric, and brain motion MRI scans were combined with a computational fluid dynamics model to investigate cerebrospinal fluid (CSF) mixing in the third cerebral ventricle of a healthy male adult. It was found that advection dominates over diffusion in most of the third ventricle. Three zones where diffusion plays an important role in the mixing process were identified. One of these zones, consisting of recessus infundibulus, recessus opticus and the adjacent regions up to commissura anterior, is likely to exist in the general population. We hypothesize that this zone may act as a buffer to flatten concentration peaks of pituitary gland hormones released into the CSF of the third ventricle. We further hypothesize that this zone may facilitate the communication between hypothalamus and the pituitary gland through the third ventricle cerebrospinal fluid by prolonging residence times of the communicated hormones.

Computational modeling of coupled blood-wall mass transport of LDL: effects of local wall shear stress

The work herein represents a novel approach for the modeling of low-density lipoprotein (LDL) transport from the artery lumen into the arterial wall, taking into account the effects of local wall shear stress (WSS) on the endothelial cell layer and its pathways of volume and solute flux. We have simulated LDL transport in an axisymmetric representation of a stenosed coronary artery, where the endothelium is represented by a three-pore model that takes into account the contributions of the vesicular pathway, normal junctions, and leaky junctions also employing the local WSS to yield the overall volume and solute flux. The fraction of leaky junctions is calculated as a function of the local WSS based on published experimental data and is used in conjunction with the pore theory to determine the transport properties of this pathway. We have found elevated levels of solute flux at low shear stress regions because of the presence of a larger number of leaky junctions compared with high shear stress regions. Accordingly, we were able to observe high LDL concentrations in the arterial wall in these low shear stress regions despite increased filtration velocity, indicating that the increase in filtration velocity is not sufficient for the convective removal of LDL.

Computational investigation of subject-specific cerebrospinal fluid flow in the third ventricle and aqueduct of Sylvius

The cerebrospinal fluid flow in the third ventricle of the brain and the aqueduct of Sylvius was studied using computational fluid dynamics (CFD) based on subject-specific boundary conditions derived from magnetic resonance imaging (MRI) scans. The flow domain geometry was reconstructed from anatomical MRI scans by manual image segmentation. The movement of the domain boundary was derived from MRI brain motion scans. Velocimetric MRI scans were used to reconstruct the velocity field at the inferior end of the aqueduct of Sylvius based on the theory of pulsatile flow in pipes. A constant pressure boundary condition was assigned at the foramina of Monro. Three main flow features were observed: a fluid jet emerging from the aqueduct of Sylvius, a moderately mobile recirculation zone above the jet and a mobile recirculation below the jet. The flow in the entire domain was laminar with a maximum Reynolds number of 340 in the aqueduct. The findings demonstrate that by combining MRI scans and CFD simulations, subject-specific detailed quantitative information of the flow field in the third ventricle and the aqueduct of Sylvius can be obtained.

Computational modeling of the mechanical behavior of the cerebrospinal fluid system

A computational fluid dynamics (CFD) model of the cerebrospinal fluid system was constructed based on a simplified geometry of the brain ventricles and their connecting pathways. The flow is driven by a prescribed sinusoidal motion of the third ventricle lateral walls, with all other boundaries being rigid. The pressure propagation between the third and lateral ventricles was examined and compared to data obtained from a similar geometry with a stenosed aqueduct. It could be shown that the pressure amplitude in the lateral ventricles increases in the presence of aqueduct stenosis. No difference in phase shift between the motion of the third ventricle walls and the pressure in the lateral ventricles because of the aqueduct stenosis could be observed. It is deduced that CFD can be used to analyze the pressure propagation and its phase shift relative to the ventricle wall motion. It is further deduced that only models that take into account the coupling between ventricles, which feature a representation of the original geometry that is as accurate as possible and which represent the ventricle boundary motion realistically, should be used to make quantitative statements on flow and pressure in the ventricular space.

Reconstruction of Cerebrospinal Fluid Flow in the Third Ventricle Based on MRI Data

A finite-volume model of the cerebrospinal fluid (CSF) system encompassing the third ventricle and the aqueduct of Sylvius was used to reconstruct CSF velocity and pressure fields based on MRI data. The flow domain geometry was obtained through segmentation of MRI brain anatomy scans. The movement of the domain walls was interpolated from brain motion MRI scans. A constant pressure boundary condition (BC) was specified at the foramina of Monro. A transient velocity BC reconstructed from velocimetric MRI scans was employed at the inferior end of the aqueduct of Sylvius. It could be shown that a combination of MRI scans and computational fluid dynamics (CFD) simulation can be used to reconstruct the flow field in the third ventricle. Pre-interventional knowledge of patient-specific CSF flow has the potential to improve neurosurgical interventions such as shunt placement in case of hydrocephalus.