Enhancing cerebral vasculature analysis with pathlength-corrected 2D angiographic parametric imaging: A feasibility study

BACKGROUND

2D angiographic parametric imaging (API) quantitatively extracts imaging biomarkers related to contrast flow and is conventionally applied to 2D digitally subtracted angiograms (DSA's). In the interventional suite, API is typically performed using 1-2 projection views and is limited by vessel overlap, foreshortening, and depth-integration of contrast motion.

PURPOSE

This work explores the use of a pathlength-correction metric to overcome the limitations of 2D-API: the primary objective was to study the effect of converting 3D contrast flow to projected contrast flow using a simulated angiographic framework created with computational fluid dynamics (CFD) simulations, thereby removing acquisition variability.

METHODS

The pathlength-correction framework was applied to in-silico angiograms, generating a reference (i.e., ground-truth) volumetric contrast distribution in four patient-specific intracranial aneurysm geometries. Biplane projections of contrast flow were created from the reference volumetric contrast distributions, assuming a cone-beam geometry. A Parker-weighted reconstruction was performed to obtain a binary representation of the vessel structure in 3D. Standard ray tracing techniques were then used to track the intersection of a ray from the focal spot with each voxel of the reconstructed vessel wall to a pixel in the detector plane. The lengths of each ray through the 3D vessel lumen were then projected along each ray-path to create a pathlength-correction map, where the pixel intensity in the detector plane corresponds to the vessel width along each source-detector ray. By dividing the projection sequences with this correction map, 2D pathlength-corrected in-silico angiograms were obtained. We then performed voxel-wise (3D) API on the ground-truth contrast distribution and compared it to pixel-wise (2D) API, both with and without pathlength correction for each biplane view. The percentage difference (PD) between the resultant API biomarkers in each dataset were calculated within the aneurysm region of interest (ROI).

RESULTS

Intensity-based API parameters, such as the area under the curve (AUC) and peak height (PH), exhibited notable changes in magnitude and spatial distribution following pathlength correction: these now accurately represent conservation of mass of injected contrast media within each arterial geometry and accurately reflect regions of stagnation and recirculation in each aneurysm ROI. Improved agreement was observed between these biomarkers in the pathlength-corrected biplane maps: the maximum PD within the aneurysm ROI is 3.3% with pathlength correction and 47.7% without pathlength correction. As expected, improved agreement with ROI-averaged ground-truth 3D counterparts was observed for all aneurysm geometries, particularly large aneurysms: the maximum PD for both AUC and PH was 5.8%. Temporal parameters (mean transit time, MTT, time-to-peak, TTP, time-to-arrival, TTA) remained unaffected after pathlength correction.

CONCLUSIONS

This study indicates that the values of intensity-based API parameters obtained with conventional 2D-API, without pathlength correction, are highly dependent on the projection orientation, and uncorrected API should be avoided for hemodynamic analysis. The proposed metric can standardize 2D API-derived biomarkers independent of projection orientation, potentially improving the diagnostic value of all acquired 2D-DSA's. Integration of a pathlength correction map into the imaging process can allow for improved interpretation of biomarkers in 2D space, which may lead to improved diagnostic accuracy during procedures involving the cerebral vasculature.

An arteriovenous mock circulatory loop and accompanying bond graph model for in vitro study of peripheral intravascular bioartificial organs

BACKGROUND

Silicon nanopore membrane-based implantable bioartificial organs are dependent on arteriovenous implantation of a mechanically robust and biocompatible hemofilter. The hemofilter acts as a low-resistance, high-flow network, with blood flow physiology similar to arteriovenous shunts commonly created for hemodialysis access. A mock circulatory loop (MCL) that mimics shunt physiology is an essential tool for refinement and durability testing of arteriovenous implantable bioartificial organs and silicon blood-interfacing membranes. We sought to develop a compact and cost-effective MCL to replicate flow conditions through an arteriovenous shunt and used data from the MCL and swine to inform a bond graph mathematical model of the physical setup.

METHODS

Flow physiology through bioartificial organ prototypes was obtained in the MCL and during extracorporeal attachment to swine for biologic comparison. The MCL was tested for stability overtime by measuring pressurewave variability over a 48-h period. Data obtained in vitro and extracorporeally informed creation of a bond graph model of the MCL.

RESULTS

The arteriovenous MCL was a cost-effective, portable system that reproduced flow rates and pressures consistent with a pulsatile arteriovenous shunt as measured in swine. MCL performance was stable over prolonged use, providing a cost-effective simulator for enhanced testing of peripherally implanted bioartificial organ prototypes. The corresponding bond graph model recapitulates MCL and animal physiology, offering a tool for further refinement of the MCL system.

Encouraging Regular Aortic Valve Opening for EVAHEART 2 LVAD Support Using Virtual Patient Hemodynamic Speed Modulation Analysis

This study focuses on investigating the EVAHEART 2 left ventricular assist device (LVAD) toward designing optimal pump speed modulation (PSM) algorithms for encouraging aortic valve (AV) flow. A custom-designed virtual patient hemodynamic model incorporating the EVAHEART 2 pressure-flow curves, cardiac chambers, and the systemic and pulmonary circulations was developed and used in this study. Several PSM waveforms were tested to evaluate their influence on the mean arterial pressure (MAP), cardiac output (CO), and AV flow for representative heart failure patients. Baseline speeds were varied from 1,600 to 2,000 rpm. For each baseline speed, the following parameters were analyzed: 1) PSM ratio (reduced speed/baseline speed), 2) PSM duration (3-7 seconds), 3) native ventricle contractility, and 4) patient MAP of 70 and 80 mm Hg. More than 2,000 rpm virtual patient scenarios were explored. A lower baseline speed (1,600 and 1,700 rpm) produced more opportunities for AV opening and more AV flow. Higher baseline speeds (1,800 and 2,000 rpm) had lower or nonexistent AV flow. When analyzing PSM ratios, a larger reduction in speed (25%) over a longer PSM (5+ seconds) duration produced the most AV flow. Lower patient MAP and increased native ventricle contractility also contributed to improving AV opening frequency and flow. This study of the EVAHEART 2 LVAD is the first to focus on leveraging PSM to enhance pulsatility and encourage AV flow. Increased AV opening frequency can benefit aortic root hemodynamics, thereby improving patient outcomes.

Modeling Gasotransmitter Availability to Brain Capillary Endothelial Cells with Ultrasound-sensitive Microbubbles

BACKGROUND

Vascular cognitive impairment and dementia results from blood components passing through disrupted blood brain barriers (BBBs). Current treatments can reduce further progress of neuronal damage but do not treat the primary cause. Instead, these treatments typically aim to temporarily disrupt the BBB. Alternatively, this study computationally assessed the feasibility of delivering carbon monoxide (CO) from ultrasound-sensitive microbubbles (MBs) as a strategy to promote BBB repair and integrity. CO can interact with heme-containing compounds within cells and promote cell growth. However, careful dose control is critical for safety and efficacy because CO also binds at high affinity to hemoglobin (Hb).

METHODS

Ultrasound activation was simulated at the internal carotid artery, and CO released from the resulting MB rupture was tracked along the shortest path to the BBB for several activation times and doses. The CO dose available to brain capillary endothelial cells (BCECs) was predicted by considering hemodynamics, mass transport, and binding kinetics.

RESULTS

The half-life of CO binding to Hb indicated that CO is available to interact with BCECs for several cardiac cycles. Further, MB and COHb concentrations would not be near toxic levels and free Hb would be available. The axisymmetric model indicated that biologically-relevant CO concentrations will be available to BCECs, and these levels can be sustained with controlled ultrasound activation. A patient-specific geometry shows that while vessel tortuosity provides a heterogeneous response, a relevant CO concentration could still be achieved.

CONCLUSIONS

This computational study demonstrates feasibility of the CO / MB strategy, and that controlled delivery is important for viability of this strategy.

Application of 1,000 fps High-Speed Angiography to In-Vitro Hemodynamic Evaluation of Left Ventricular Assist Device Outflow Graft Configurations

Left ventricular assist device (LVAD)-induced hemodynamics are characterized by fast-moving flow with large variations in velocity, making quantitative assessments difficult with existing imaging methods. This study demonstrates the ability of 1,000 fps high-speed angiography (HSA) to quantify the effect of the surgical implantation angle of a LVAD outflow graft on the hemodynamics within the ascending aorta in vitro . High-speed angiography was performed on patient-derived, three-dimensional-printed optically opaque aortic models using a nonsoluble contrast media, ethiodol, as a flow tracer. Outflow graft configuration angles of 45° and 90° with respect to the central aortic axis were considered. Projected velocity distributions were calculated from the high-speed experimental sequences using two methods: a physics-based optical flow algorithm and tracking of radio-opaque particles. Particle trajectories were also used to evaluate accumulated shear stress. Results were then compared with computational fluid dynamics (CFD) simulations to confirm the results of the high-speed imaging method. Flow patterns derived from HSA coincided with the impingement regions and recirculation zones formed in the aortic root as seen in the CFD for both graft configurations. Compared with the 45° graft, the 90° configuration resulted in 81% higher two-dimensional-projected velocities (over 100 cm/s) along the contralateral wall of the aorta. Both graft configurations suggest elevated accumulated shear stresses along individual trajectories. Compared with CFD simulations, HSA successfully characterized the fast-moving flow and hemodynamics in each LVAD graft configuration in vitro , demonstrating the potential utility of this technology as a quantitative imaging modality.

A Computational Hemodynamics Approach to Left Ventricular Assist Device (LVAD) Optimization Validated in a Large Patient Cohort

With increasing use of left ventricular assist devices (LVAD) it is critical to devise strategies to optimize LVAD speed while controlling mean arterial pressure (MAP) and flow according to patient physiology. The complex interdependency between LVAD speed, MAP, and flow frequently makes optimization difficult under clinical conditions. We propose a method to guide this procedure in silico, narrowing the conditions to test clinically. A computational model of the circulatory network that simulates HF and LVAD support, incorporating LVAD pressure-flow curves was applied retrospectively to anonymized patient hemodynamics data from the University of Washington Medical Center. MAP management on 61 patient-specific computational models with a target of 70 mm Hg, resulting flow for a given LVAD speed was analyzed, and compared to a target output of 5 L/min. Before performing virtual MAP management, 51% had a MAP>70 mm Hg and CO>5 L/min, and 33% had a MAP>70 mm Hg and CO<5 L/min. After changing systemic resistance to meet the MAP target (without adjusting LVAD speed), 84% of cases resulted in CO higher than 5 L/min, with a median CO of 6.79 L/min, using the computational predictive model. Blood pressure management alone is insufficient in meeting both MAP and CO targets, due to the risk of hypervolemia, and requires appropriate LVAD speed optimization to achieve both targets, while preserving right heart health. Such computational tools can narrow down conditions to be tested for each patient, providing significant insight into the pump-patient interplay. LVAD hemodynamic optimization has the potential to reduce complications and improve outcomes.

Lagrangian Trajectory Simulation of Platelets and Synchrotron Microtomography Augment Hemodynamic Analysis of Intracranial Aneurysms Treated With Embolic Coils

As frequency of endovascular treatments for intracranial aneurysms increases, there is a growing need to understand the mechanisms for coil embolization failure. Computational fluid dynamics (CFD) modeling often simplifies modeling the endovascular coils as a homogeneous porous medium (PM), and focuses on the vascular wall endothelium, not considering the biomechanical environment of platelets. These assumptions limit the accuracy of computations for treatment predictions. We present a rigorous analysis using X-ray microtomographic imaging of the coils and a combination of Lagrangian (platelet) and Eulerian (endothelium) metrics. Four patient-specific, anatomically accurate in vitro flow phantoms of aneurysms are treated with the same patient-specific endovascular coils. Synchrotron tomography scans of the coil mass morphology are obtained. Aneurysmal hemodynamics are computationally simulated before and after coiling, using patient-specific velocity/pressure measurements. For each patient, we analyze the trajectories of thousands of platelets during several cardiac cycles, and calculate residence times (RTs) and shear exposure, relevant to thrombus formation. We quantify the inconsistencies of the PM approach, comparing them with coil-resolved (CR) simulations, showing the under- or overestimation of key hemodynamic metrics used to predict treatment outcomes. We fully characterize aneurysmal hemodynamics with converged statistics of platelet RT and shear stress history (SH), to augment the traditional wall shear stress (WSS) on the vascular endothelium. Incorporating microtomographic scans of coil morphology into hemodynamic analysis of coiled intracranial aneurysms, and augmenting traditional analysis with Lagrangian platelet metrics improves CFD predictions, and raises the potential for understanding and clinical translation of computational hemodynamics for intracranial aneurysm treatment outcomes.

Small Left Ventricular Size Is an Independent Risk Factor for Ventricular Assist Device Thrombosis

The prevalence of ventricular assist device (VAD) therapy has continued to increase due to a stagnant donor supply and growing advanced heart failure (HF) population. We hypothesize that left ventricular (LV) size strongly influences biocompatibility and risk of thrombosis. Unsteady computational fluid dynamics (CFD) was used in conjunction with patient-derived computational modeling and virtual surgery with a standard, apically implanted inflow cannula. A dual-focus approach of evaluating thrombogenicity was employed: platelet-based metrics to characterize the platelet environment and flow-based metrics to investigate hemodynamics. Left ventricular end-diastolic dimensions (LVEDds) ranging from 4.5 to 6.5 cm were studied and ranked according to relative thrombogenic potential. Over 150,000 platelets were individually tracked in each LV model over 15 cardiac cycles. As LV size decreased, platelets experienced markedly increased shear stress histories (SHs), whereas platelet residence time (RT) in the LV increased with size. The complex interplay between increased SH and longer RT has profound implications on thrombogenicity, with a significantly higher proportion of platelets in small LVs having long RT times and being subjected to high SH, contributing to thrombus formation. Our data suggest that small LV size, rather than decreased VAD speed, is the primary pathologic mechanism responsible for the increased incidence of thrombosis observed in VAD patients with small LVs.

Platelet Dynamics and Hemodynamics of Cerebral Aneurysms Treated with Flow-Diverting Stents

Flow-diverting stents (FDS) are used to treat cerebral aneurysms. They promote the formation of a stable thrombus within the aneurysmal sac and, if successful, isolate the aneurysmal dome from mechanical stresses to prevent rupture. Platelet activation, a mechanism necessary for thrombus formation, is known to respond to biomechanical stimuli, particularly to the platelets' residence time and shear stress exposure. Currently, there is no reliable method for predicting FDS treatment outcomes, either a priori or after the procedure. Eulerian computational fluid dynamic (CFD) studies of aneurysmal flow have searched for predictors of endovascular treatment outcome; however, the hemodynamics of thrombus formation cannot be fully understood without considering the platelets' trajectories and their mechanics-triggered activation. Lagrangian analysis of the fluid mechanics in the aneurysmal vasculature provides novel metrics by tracking the platelets' residence time (RT) and shear history (SH). Eulerian and Lagrangian parameters are compared for 19 patient-specific cases, both pre- and post-treatment, to assess the degree of change caused by the FDS and subsequent treatment efficacy.

Left Ventricular Assist Device Inflow Cannula Angle and Thrombosis Risk

BACKGROUND

As heart failure prevalence continues to increase in the setting of a static donor supply, left ventricular assist device (LVAD) therapy for end-stage heart failure continues to grow. Anecdotal evidence suggests that malalignment of the LVAD inflow cannula may increase thrombosis risk, but this effect has not been explored mechanistically or quantified statistically. Our objective is to elucidate the impact of surgical angulation of the inflow cannula on thrombogenicity.

METHODS AND RESULTS

Unsteady computational fluid dynamics is used in conjunction with computational modeling and virtual surgery to model flow through the left ventricle for 5 different inflow cannula angulations. We use a holistic approach to evaluate thrombogenicity: platelet-based (Lagrangian) metrics to evaluate the platelet mechanical environment, combined with flow-based (Eulerian) metrics to investigate intraventricular hemodynamics. The thrombogenic potential of each LVAD inflow cannula angulation is quantitatively evaluated based on platelet shear stress history and residence time. Intraventricular hemodynamics are strongly influenced by LVAD inflow cannula angulation. Platelet behavior indicates elevated thrombogenic potential for certain inflow cannula angles, potentially leading to platelet activation. Our analysis demonstrates that the optimal range of inflow angulation is within 0±7° of the left ventricular apical axis.

CONCLUSIONS

Angulation of the inflow cannula >7° from the apical axis (axis connecting mitral valve and ventricular apex) leads to markedly unfavorable hemodynamics as determined by computational fluid dynamics. Computational hemodynamic simulations incorporating Lagrangian and Eulerian metrics are a powerful tool for studying optimization of LVAD implantation strategies, with the long-term potential of improving outcomes.

Accuracy of Doppler blood pressure measurement in continuous-flow left ventricular assist device patients

Aims

Accurate blood pressure (BP) measurement in continuous‐flow ventricular assist device (CF‐VAD) patients is imperative to reduce stroke risk. This study assesses the accuracy of the Doppler opening pressure method compared with the gold standard arterial line method in CF‐VAD patients.

Methods and results

In a longitudinal cohort of HeartMate II and HVAD patients, arterial line BP and simultaneously measured Doppler opening pressure were obtained. Overall correlation, agreement between Doppler opening pressure and arterial line mean vs. systolic pressure, and the effect of arterial pulsatility on the accuracy of Doppler opening pressure were analysed. A total of 1933 pairs of Doppler opening pressure and arterial line pressure readings within 1 min of each other were identified in 154 patients (20% women, mean age 55 ± 15, 50% HeartMate II and 50% HVAD). Doppler opening pressure had good correlation with invasive mean arterial pressure (r = 0.742, P < 0.0001) and more closely approximated mean than systolic BP (mean error 2.4 vs. −8.4 mmHg). Arterial pulsatility did not have a clinically significant effect on the accuracy of the Doppler opening pressure method.

Conclusions

Doppler opening pressure should be the standard non‐invasive method of BP measurement in CF‐VAD patients.

Blood damage in Left Ventricular Assist Devices: Pump thrombosis or system thrombosis?

INTRODUCTION:

Despite significant technical advancements in the design and manufacture of Left Ventricular Assist Devices, post-implant thrombotic and thromboembolic complications continue to affect long-term outcomes. Previous efforts, aimed at optimizing pump design as a means of reducing supraphysiologic shear stresses generated within the pump and associated prothrombotic shear-mediated platelet injury, have only partially altered the device hemocompatibility.

METHODS:

We examined hemodynamic mechanisms that synergize with hypershear within the pump to contribute to the thrombogenic potential of the overall Left Ventricular Assist Device system.

RESULTS:

Numerical simulations of blood flow in differing regions of the Left Ventricular Assist Device system, that is the diseased native left ventricle, the pump inflow cannula, the impeller, the outflow graft and the anastomosed downstream aorta, reveal that prothrombotic hemodynamic conditions might occur at these specific sites. Furthermore, we show that beyond hypershear, additional hemodynamic abnormalities exist within the pump, which may elicit platelet activation, such as recirculation zones and stagnant platelet trajectories. We also provide evidences that particular Left Ventricular Assist Device implantation configurations and specific post-implant patient management strategies, such as those allowing aortic valve opening, are more hemodynamically favorable and reduce the thrombotic risk.

CONCLUSION:

We extend the perspective of pump thrombosis secondary to the supraphysiologic shear stress environment of the pump to one of Left Ventricular Assist Device system thrombosis, raising the importance of comprehensive characterization of the different prothrombotic risk factors of the total system as the target to achieve enhanced hemocompatibility and improved clinical outcomes.

Blood damage in Left Ventricular Assist Devices: Pump thrombosis or system thrombosis?

Introduction:

Despite significant technical advancements in the design and manufacture of Left Ventricular Assist Devices, post-implant thrombotic and thromboembolic complications continue to affect long-term outcomes. Previous efforts, aimed at optimizing pump design as a means of reducing supraphysiologic shear stresses generated within the pump and associated prothrombotic shear-mediated platelet injury, have only partially altered the device hemocompatibility.

Methods:

We examined hemodynamic mechanisms that synergize with hypershear within the pump to contribute to the thrombogenic potential of the overall Left Ventricular Assist Device system.

Results:

Numerical simulations of blood flow in differing regions of the Left Ventricular Assist Device system, that is the diseased native left ventricle, the pump inflow cannula, the impeller, the outflow graft and the anastomosed downstream aorta, reveal that prothrombotic hemodynamic conditions might occur at these specific sites. Furthermore, we show that beyond hypershear, additional hemodynamic abnormalities exist within the pump, which may elicit platelet activation, such as recirculation zones and stagnant platelet trajectories. We also provide evidences that particular Left Ventricular Assist Device implantation configurations and specific post-implant patient management strategies, such as those allowing aortic valve opening, are more hemodynamically favorable and reduce the thrombotic risk.

Conclusion:

We extend the perspective of pump thrombosis secondary to the supraphysiologic shear stress environment of the pump to one of Left Ventricular Assist Device system thrombosis, raising the importance of comprehensive characterization of the different prothrombotic risk factors of the total system as the target to achieve enhanced hemocompatibility and improved clinical outcomes.

Hyperglycemia Alters the Structure and Hemodynamics of the Developing Embryonic Heart

Congenital heart defects (CHDs) represent the most common form of human birth defects; approximately one-third of heart defects involve malformations of the outflow tract (OFT). Maternal diabetes increases the risk of CHD by 3-5 fold. During heart organogenesis, little is known about the effects of hyperglycemia on hemodynamics, which are critical to normal heart development. Heart development prior to septation in the chick embryo was studied under hyperglycemic conditions. Sustained hyperglycemic conditions were induced, raising the average plasma glucose concentration from 70 mg/dL to 180 mg/dL, akin to the fasting plasma glucose of a patient with diabetes. The OFTs were assessed for structural and hemodynamic alterations using optical coherence tomography (OCT), confocal microscopy, and microcomputed tomography. In hyperglycemic embryos, the endocardial cushions of the proximal OFT were asymmetric, and the OFTs curvature and torsion were significantly altered. The blood flow velocity through the OFT of hyperglycemic embryos was significantly decreased, including flow reversal in 30% of the cardiac cycle. Thus, hyperglycemia at the onset of gestation results in asymmetric proximal endocardial cushions, abnormal OFT curvature, and altered hemodynamics in the developing heart. If present in humans, these results may identify early developmental alterations that contribute to the increased risk for cardiac malformations in babies from diabetic mothers.

LVAD Outflow Graft Angle and Thrombosis Risk

This study quantifies thrombogenic potential (TP) of a wide range of left ventricular assist device (LVAD) outflow graft anastomosis angles through state-of-the-art techniques: 3D imaged-based patient-specific models created via virtual surgery and unsteady computational fluid dynamics with Lagrangian particle tracking. This study aims at clarifying the influence of a single parameter (outflow graft angle) on the thrombogenesis associated with flow patterns in the aortic root after LVAD implantation. This is an important and poorly-understood aspect of LVAD therapy, because several studies have shown strong inter and intrapatient thrombogenic variability and current LVAD implantation strategies do not incorporate outflow graft angle optimization. Accurate platelet-level investigation, enabled by statistical treatment of outliers in Lagrangian particle tracking, demonstrates a strong influence of outflow graft anastomoses angle on thrombogenicity (platelet residence times and activation state characterized by shear stress accumulation) with significantly reduced TP for acutely-angled anastomosed outflow grafts. The methodology presented in this study provides a device-neutral platform for conducting comprehensive thrombogenicity evaluation of LVAD surgical configurations, empowering optimal patient-focused surgical strategies for long-term treatment and care for advanced heart failure patients.

Blood flow through the embryonic heart outflow tract during cardiac looping in HH13-HH18 chicken embryos

Blood flow is inherently linked to embryonic cardiac development, as haemodynamic forces exerted by flow stimulate mechanotransduction mechanisms that modulate cardiac growth and remodelling. This study evaluated blood flow in the embryonic heart outflow tract (OFT) during normal development at each stage between HH13 and HH18 in chicken embryos, in order to characterize changes in haemodynamic conditions during critical cardiac looping transformations. Two-dimensional optical coherence tomography was used to simultaneously acquire both structural and Doppler flow images, in order to extract blood flow velocity and structural information and estimate haemodynamic measures. From HH13 to HH18, peak blood flow rate increased by 2.4-fold and stroke volume increased by 2.1-fold. Wall shear rate (WSR) and lumen diameter data suggest that changes in blood flow during HH13-HH18 may induce a shear-mediated vasodilation response in the OFT. Embryo-specific four-dimensional computational fluid dynamics modelling at HH13 and HH18 complemented experimental observations and indicated heterogeneous WSR distributions over the OFT. Characterizing changes in haemodynamics during cardiac looping will help us better understand the way normal blood flow impacts proper cardiac development.

Effect of Outflow Tract Banding on Embryonic Cardiac Hemodynamics

We analyzed heart wall motion and blood flow dynamics in chicken embryos using in vivo optical coherence tomography (OCT) imaging and computational fluid dynamics (CFD) embryo-specific modeling. We focused on the heart outflow tract (OFT) region of day 3 embryos, and compared normal (control) conditions to conditions after performing an OFT banding intervention, which alters hemodynamics in the embryonic heart and vasculature. We found that hemodynamics and cardiac wall motion in the OFT are affected by banding in ways that might not be intuitive a priori. In addition to the expected increase in ventricular blood pressure, and increase blood flow velocity and, thus, wall shear stress (WSS) at the band site, the characteristic peristaltic-like motion of the OFT was altered, further affecting flow and WSS. Myocardial contractility, however, was affected only close to the band site due to the physical restriction on wall motion imposed by the band. WSS were heterogeneously distributed in both normal and banded OFTs. Our results show how banding affects cardiac mechanics and can lead, in the future, to a better understanding of mechanisms by which altered blood flow conditions affect cardiac development leading to congenital heart disease.

4D subject-specific inverse modeling of the chick embryonic heart outflow tract hemodynamics

Blood flow plays a critical role in regulating embryonic cardiac growth and development, with altered flow leading to congenital heart disease. Progress in the field, however, is hindered by a lack of quantification of hemodynamic conditions in the developing heart. In this study, we present a methodology to quantify blood flow dynamics in the embryonic heart using subject-specific computational fluid dynamics (CFD) models. While the methodology is general, we focused on a model of the chick embryonic heart outflow tract (OFT), which distally connects the heart to the arterial system, and is the region of origin of many congenital cardiac defects. Using structural and Doppler velocity data collected from optical coherence tomography, we generated 4D ([Formula: see text]) embryo-specific CFD models of the heart OFT. To replicate the blood flow dynamics over time during the cardiac cycle, we developed an iterative inverse-method optimization algorithm, which determines the CFD model boundary conditions such that differences between computed velocities and measured velocities at one point within the OFT lumen are minimized. Results from our developed CFD model agree with previously measured hemodynamics in the OFT. Further, computed velocities and measured velocities differ by [Formula: see text]15 % at locations that were not used in the optimization, validating the model. The presented methodology can be used in quantifications of embryonic cardiac hemodynamics under normal and altered blood flow conditions, enabling an in-depth quantitative study of how blood flow influences cardiac development.

Alterations in cancer cell mechanical properties after fluid shear stress exposure: a micropipette aspiration study

Over 90% of cancer deaths result not from primary tumor development, but from metastatic tumors that arise after cancer cells circulate to distal sites via the circulatory system. While it is known that metastasis is an inefficient process, the effect of hemodynamic parameters such as fluid shear stress (FSS) on the viability and efficacy of metastasis is not well understood. Recent work has shown that select cancer cells may be able to survive and possibly even adapt to FSS in vitro. The current research seeks to characterize the effect of FSS on the mechanical properties of suspended cancer cells in vitro. Nontransformed prostate epithelial cells (PrEC LH) and transformed prostate cancer cells (PC-3) were used in this study. The Young's modulus was determined using micropipette aspiration. We examined cells in suspension but not exposed to FSS (unsheared) and immediately after exposure to high (6,400 dyn/cm²) and low (510 dyn/cm²) FSS. The PrEC LH cells were ~140% stiffer than the PC-3 cells not exposed to FSS. Post-FSS exposure, there was an increase of ~77% in Young's modulus after exposure to high FSS and a ~47% increase in Young's modulus after exposure to low FSS for the PC-3 cells. There was no significant change in the Young's modulus of PrEC LH cells post-FSS exposure. Our findings indicate that cancer cells adapt to FSS, with an increased Young's modulus being one of the adaptive responses, and that this adaptation is specific only to PC-3 cells and is not seen in PrEC LH cells. Moreover, this adaptation appears to be graded in response to the magnitude of FSS experienced by the cancer cells. This is the first study investigating the effect of FSS on the mechanical properties of cancer cells in suspension, and may provide significant insights into the mechanism by which some select cancer cells may survive in the circulation, ultimately leading to metastasis at distal sites. Our findings suggest that biomechanical analysis of cancer cells could aid in identifying and diagnosing cancer in the future.

Red blood cell flow in the cardiovascular system: a fluid dynamics perspective

The dynamics of red blood cells (RBCs) is one of the major aspects of the cardiovascular system that has been studied intensively in the past few decades. The dynamics of biconcave RBCs are thought to have major influences in cardiovascular diseases, the problems associated with cardiovascular assistive devices, and the determination of blood rheology and properties. This article provides an overview of the works that have been accomplished in the past few decades and aim to study the dynamics of RBCs under different flow conditions. While significant progress has been made in both experimental and numerical studies, a detailed understanding of the behavior of RBCs is still faced with many challenges. Experimentally, the size of RBCs is considered to be a major limitation that allows measurements to be performed under conditions similar to physiological conditions. In numerical computations, researchers still are working to develop a model that can cover the details of the RBC mechanics as it deforms and moves in the bloodstream. Moreover, most of reported computational models have been confined to the behavior of a single RBC in 2-dimensional domains. Advanced models are yet to be developed for accurate description of RBC dynamics under physiological flow conditions in 3-dimensional regimes.