Pharmacological Modulation of Hemodynamics in Adult Zebrafish In Vivo

INTRODUCTION: Hemodynamic parameters in zebrafish receive increasing attention because of their important role in cardiovascular processes such as atherosclerosis, hematopoiesis, sprouting and intussusceptive angiogenesis. To study underlying mechanisms, the precise modulation of parameters like blood flow velocity or shear stress is centrally important. Questions related to blood flow have been addressed in the past in either embryonic or ex vivo-zebrafish models but little information is available for adult animals. Here we describe a pharmacological approach to modulate cardiac and hemodynamic parameters in adult zebrafish in vivo. MATERIALS AND METHODS: Adult zebrafish were paralyzed and orally perfused with salt water. The drugs isoprenaline and sodium nitroprusside were directly applied with the perfusate, thus closely resembling the preferred method for drug delivery in zebrafish, namely within the water. Drug effects on the heart and on blood flow in the submental vein were studied using electrocardiograms, in vivo-microscopy and mathematical flow simulations. RESULTS: Under control conditions, heart rate, blood flow velocity and shear stress varied less than ± 5%. Maximal chronotropic effects of isoprenaline were achieved at a concentration of 50 μmol/L, where it increased the heart rate by 22.6 ± 1.3% (n = 4; p < 0.0001). Blood flow velocity and shear stress in the submental vein were not significantly increased. Sodium nitroprusside at 1 mmol/L did not alter the heart rate but increased blood flow velocity by 110.46 ± 19.64% (p = 0.01) and shear stress by 117.96 ± 23.65% (n = 9; p = 0.03). DISCUSSION: In this study, we demonstrate that cardiac and hemodynamic parameters in adult zebrafish can be efficiently modulated by isoprenaline and sodium nitroprusside. Together with the suitability of the zebrafish for in vivo-microscopy and genetic modifications, the methodology described permits studying biological processes that are dependent on hemodynamic alterations.

Maximal exercise hemodynamics and risk of mortality in apparently healthy men and women

The use of exercise electrocardiography (ECG) to detect latent coronary heart disease (CHD) is discouraged in apparently healthy populations because of low sensitivity. These recommendations however, are based on the efficacy of evaluation of ischemia (ST segment changes) with little regard for other measures of cardiac function that are available during exertion. The purpose of this investigation was to determine the association of maximal exercise hemodynamic responses with risk of mortality due to all-causes, cardiovascular disease (CVD), and coronary heart disease (CHD) in apparently healthy individuals. Study participants were 20,387 men (mean age = 42.2 years) and 6,234 women (mean age = 41.9 years) patients of a preventive medicine center in Dallas, TX examined between 1971 and 1989. During an average of 8.1 years of follow-up, there were 348 deaths in men and 66 deaths in women. In men, age-adjusted all-cause death rates (per 10,000 person years) across quartiles of maximal systolic blood pressure (SBP) (low to high) were: 18.2, 16.2, 23.8, and 24.6 (p for trend $<$0.001). Corresponding rates for maximal heart rate were: 28.9, 15.9, 18.4, and 15.1 (p trend $<$0.001). After adjustment for confounding variables including age, resting systolic pressure, serum cholesterol and glucose, body mass index, smoking status, physical fitness and family history of CVD, risks (and 95% confidence interval (CI)) of all-cause mortality for quartiles of maximal SBP, relative to the lowest quartile, were: 0.96 (0.70-1.33), 1.36 (1.01-1.85), and 1.37 (0.98-1.92) for quartiles 2-4 respectively. Similar risks for maximal heart rate were: 0.61 (0.44-0.85), 0.69 (0.51-0.93), and 0.60 (0.41-0.87). No associations were noted between maximal exercise rate-pressure product mortality. Similar results were seen for risk of CVD and CHD death. In women, similar trends in age-adjusted all-cause and CVD death rates across maximal SBP and heart rate categories were observed. Sensitivity of the exercise test in predicting mortality was enhanced when ECG results were evaluated together with maximal exercise SBP or heart rate with a concomitant decrease in specificity. Positive predictive values were not improved. The efficacy of the exercise test in predicting mortality in apparently healthy men and women was not enhanced by using maximal exercise hemodynamic responses. These results suggest that an exaggerated systolic blood pressure or an attenuated heart rate response to maximal exercise are risk factors for mortality in apparently healthy individuals. ^

Response of B-type natriuretic pepticle to exercise in hypertensive patients with suspected diastolic heart failure: Correlation with cardiac function, hemodynamics, and workload

Background Diastolic heart failure (DHF) is characterized by dyspnea due to increased left ventricular (LV) filling pressures during stress. We sought the relationship of exercise-induced increases in B-type natriuretic peptide (BNP) to LV filling pressures and parameters of cardiovascular performance in suspected DHF. Methods Twenty-six treated hypertensive patients with suspected DHF (exertional dyspnea, LV ejection fraction >50%, and diastolic dysfunction) underwent maximal exercise echocardiography using the Bruce protocol. BNP, transmitral Doppler, and tissue Doppler for systolic (So) and early (Ea) and late (Aa) diastolic mitral annular velocities were obtained at rest and peak stress. LV filling pressures were estimated with E/Ea ratios. Results Resting BNP correlated with resting pulse pressure (r=0.45, P=0.02). Maximal exercise performance (4.6 +/- 2.5min) was limited by dyspnea. Blood pressure increased with exercise (from 143 +/- 19/88 +/- 8 to 191 +/- 22/90 +/- 10 mm Hg); 13 patients (50%) had a hypertensive response. Peak exercise BNP correlated with peak transmitral E velocity (r = 0.41, P <.05) and peak heart rate (r = -0.40, P <.05). BNP increased with exercise (from 48 57 to 74 97 pg/mL, P =.007), and the increment of BNP with exercise was associated with maximal workload and peak exercise So, Ea, and Aa (P <.01 for all). Filling pressures, approximated by lateral E/Ea ratio, increased with exercise (7.7 +/- 2.0 to 10.0 +/- 4.8, P <.01). BNP was higher in patients with possibly elevated filling pressures at peak exercise (E/Ea >10) compared to those with normal pressures (123 +/- 124 vs 45 +/- 71 pg/mL, P =.027). Conclusions Augmentation of BNP with exercise in hypertensive patients with suspected DHF is associated with better exercise capacity, LV systolic and diastolic function, and left atrial function. Peak exercise BNP levels may identify exercise-induced elevation of filling pressures in DHF.

Influence of Outlet Boundary Conditions on Cerebrovascular Aneurysm Hemodynamics

Computational fluid dynamic (CFD) studies of blood flow in cerebrovascular aneurysms have potential to improve patient treatment planning by enabling clinicians and engineers to model patient-specific geometries and compute predictors and risks prior to neurovascular intervention. However, the use of patient-specific computational models in clinical settings is unfeasible due to their complexity, computationally intensive and time-consuming nature. An important factor contributing to this challenge is the choice of outlet boundary conditions, which often involves a trade-off between physiological accuracy, patient-specificity, simplicity and speed. In this study, we analyze how resistance and impedance outlet boundary conditions affect blood flow velocities, wall shear stresses and pressure distributions in a patient-specific model of a cerebrovascular aneurysm. We also use geometrical manipulation techniques to obtain a model of the patient’s vasculature prior to aneurysm development, and study how forces and stresses may have been involved in the initiation of aneurysm growth. Our CFD results show that the nature of the prescribed outlet boundary conditions is not as important as the relative distributions of blood flow through each outlet branch. As long as the appropriate parameters are chosen to keep these flow distributions consistent with physiology, resistance boundary conditions, which are simpler, easier to use and more practical than their impedance counterparts, are sufficient to study aneurysm pathophysiology, since they predict very similar wall shear stresses, time-averaged wall shear stresses, time-averaged pressures, and blood flow patterns and velocities. The only situations where the use of impedance boundary conditions should be prioritized is if pressure waveforms are being analyzed, or if local pressure distributions are being evaluated at specific time points, especially at peak systole, where the use of resistance boundary conditions leads to unnaturally large pressure pulses. In addition, we show that in this specific patient, the region of the blood vessel where the neck of the aneurysm developed was subject to abnormally high wall shear stresses, and that regions surrounding blebs on the aneurysmal surface were subject to low, oscillatory wall shear stresses. Computational models using resistance outlet boundary conditions may be suitable to study patient-specific aneurysm progression in a clinical setting, although several other challenges must be addressed before these tools can be applied clinically.

Association of myocardial T1-mapping CMR with hemodynamics and RV performance in pulmonary hypertension

Early detection of right ventricular (RV) involvement in chronic pulmonary hypertension (PH) is essential due to prognostic implications. T1 mapping by cardiac magnetic resonance (CMR) has emerged as a noninvasive technique for extracellular volume fraction (ECV) quantification. We assessed the association of myocardial native T1 time and equilibrium contrast ECV (Eq-ECV) at the RV insertion points with pulmonary hemodynamics and RV performance in an experimental model of chronic PH. Right heart catheterization followed by immediate CMR was performed on 38 pigs with chronic PH (generated by surgical pulmonary vein banding) and 6 sham-operated controls. Native T1 and Eq-ECV values at the RV insertion points were both significantly higher in banded animals than in controls and showed significant correlation with pulmonary hemodynamics, RV arterial coupling, and RV performance. Eq-ECV values also increased before overt RV systolic dysfunction, offering potential for the early detection of myocardial involvement in chronic PH.

A compact mock circulation loop for the in vitro testing of cardiovascular devices

In vitro cardiovascular device performance evaluation in a mock circulation loop (MCL) is a necessary step prior to in vivo testing.A MCL that accurately represents the physiology of the cardiovascular system accelerates the assessment of the device’s ability to treat pathological conditions. To serve this purpose, a compact MCL measuring 600 ¥ 600 ¥ 600 mm (L ¥ W¥ H) was constructed in conjunction with a computer mathematical simulation.This approach allowed the effective selection of physical loop characteristics, such as pneumatic drive parameters, to create pressure and flow, and pipe dimensions to replicate the resistance, compliance, and fluid inertia of the native cardiovascular system. The resulting five-element MCL reproduced the physiological hemodynamics of a healthy and failing heart by altering ventricle contractility, vascular resistance/compliance, heart rate, and vascular volume. The effects of interpatient anatomical variability, such as septal defects and valvular disease, were also assessed. Cardiovascular hemodynamic pressures (arterial, venous, atrial, ventricular), flows (systemic, bronchial, pulmonary), and volumes (ventricular, stroke) were analyzed in real time. The objective of this study is to describe the developmental stages of the compact MCL and demonstrate its value as a research tool for the accelerated development of cardiovascular devices.

Reversal of acute coagulopathy during hypotensive resuscitation using small-volume 7.5% NaCl adenocaine and Mg2+ in the rat model of severe hemorrhagic shock

OBJECTIVE: : Acute traumatic coagulopathy occurs early in hemorrhagic trauma and is a major contributor to mortality and morbidity. Our aim was to examine the effect of small-volume 7.5% NaCl adenocaine (adenosine and lidocaine, adenocaine) and Mg on hypotensive resuscitation and coagulopathy in the rat model of severe hemorrhagic shock. DESIGN: : Prospective randomized laboratory investigation. SUBJECTS: : A total of 68 male Sprague Dawley Rats. INTERVENTION: : Post-hemorrhagic shock treatment for acute traumatic coagulopathy. MEASUREMENTS AND METHODS: : Nonheparinized male Sprague-Dawley rats (300-450 g, n = 68) were randomly assigned to either: 1) untreated; 2) 7.5% NaCl; 3) 7.5% NaCl adenocaine; 4) 7.5% NaCl Mg; or 5) 7.5% NaCl adenocaine/Mg. Hemorrhagic shock was induced by phlebotomy to mean arterial pressure of 35-40 mm Hg for 20 mins (~40% blood loss), and animals were left in shock for 60 mins. Bolus (0.3 mL) was injected into the femoral vein and hemodynamics monitored. Blood was collected in Na citrate (3.2%) tubes, centrifuged, and the plasma snap frozen in liquid N2 and stored at -80°C. Coagulation was assessed using activated partial thromboplastin times and prothrombin times. RESULTS: : Small-volume 7.5% NaCl adenocaine and 7.5% NaCl adenocaine/Mg were the only two groups that gradually increased mean arterial pressure 1.6-fold from 38-39 mm Hg to 52 and 64 mm Hg, respectively, at 60 mins (p < .05). Baseline plasma activated partial thromboplastin time was 17 ± 0.5 secs and increased to 63 ± 21 secs after bleeding time, and 217 ± 32 secs after 60-min shock. At 60-min resuscitation, activated partial thromboplastin time values for untreated, 7.5% NaCl, 7.5% NaCl/Mg, and 7.5% NaCl adenocaine rats were 269 ± 31 secs, 262 ± 38 secs, 150 ± 43 secs, and 244 ± 38 secs, respectively. In contrast, activated partial thromboplastin time for 7.5% NaCl adenocaine/Mg was 24 ± 2 secs (p < .05). Baseline prothrombin time was 28 ± 0.8 secs (n = 8) and followed a similar pattern of correction. CONCLUSIONS: : Plasma activated partial thromboplastin time and prothrombin time increased over 10-fold during the bleed and shock periods prior to resuscitation, and a small-volume (~1 mL/kg) IV bolus of 7.5% NaCl AL/Mg was the only treatment group that raised mean arterial pressure into the permissive range and returned activated partial thromboplastin time and prothrombin time clotting times to baseline at 60 mins.

Anatomic fitting of total artificial hearts for in vivo evaluation

Successful anatomic fitting of a total artificial heart (TAH) is vital to achieve optimal pump hemodynamics after device implantation. Although many anatomic fitting studies have been completed in humans prior to clinical trials, few reports exist that detail the experience in animals for in vivo device evaluation. Optimal hemodynamics are crucial throughout the in vivo phase to direct design iterations and ultimately validate device performance prior to pivotal human trials. In vivo evaluation in a sheep model allows a realistically sized representation of a smaller patient, for which smaller third-generation TAHs have the potential to treat. Our study aimed to assess the anatomic fit of a single device rotary TAH in sheep prior to animal trials and to use the data to develop a threedimensional, computer-aided design (CAD)-operated anatomic fitting tool for future TAH development. Following excision of the native ventricles above the atrio-ventricular groove, a prototype TAH was inserted within the chest cavity of six sheep (28–40 kg).Adjustable rods representing inlet and outlet conduits were oriented toward the center of each atrial chamber and the great vessels, with conduit lengths and angles recorded for future analysis. A threedimensional, CAD-operated anatomic fitting tool was then developed, based on the results of this study, and used to determine the inflow and outflow conduit orientation of the TAH. The mean diameters of the sheep left atrium, right atrium, aorta, and pulmonary artery were 39, 33, 12, and 11 mm, respectively. The center-to-center distance and outer-edge-to-outer-edge distance between the atria, found to be 39 ± 9 mm and 72 ± 17 mm in this study, were identified as the most critical geometries for successful TAH connection. This geometric constraint restricts the maximum separation allowable between left and right inlet ports of a TAH to ensure successful alignment within the available atrial circumference.

Computational fluid dynamic analysis of intracranial aneurysmal bleb formation

Background The management of unruptured aneurysms is controversial with the decision to treat influenced by aneurysm characteristics including size and morphology. Aneurysmal bleb formation is thought to be associated with an increased risk of rupture. Objective To correlate computational fluid dynamic (CFD) indices with bleb formation. Methods Anatomical models were constructed from three-dimensional rotational angiogram (3DRA) data in 27 patients with cerebral aneurysms harbouring single blebs. Additional models representing the aneurysm before bleb formation were constructed by digitally removing the bleb. We characterised haemodynamic features of models both with and without the bleb using CFDs. Flow structure, wall shear stress (WSS), pressure and oscillatory shear index (OSI) were analysed. Results There was a statistically significant association between bleb location at or adjacent to the point of maximal WSS (74.1%, p=0.019), irrespective of rupture status. Aneurysmal blebs were related to the inflow or outflow jet in 88.9% of cases (p<0.001) whilst 11.1% were unrelated. Maximal wall pressure and OSI were not significantly related to bleb location. The bleb region attained a lower WSS following its formation in 96.3% of cases (p<0.001) and was also lower than the average aneurysm WSS in 86% of cases (p<0.001). Conclusion Cerebral aneurysm blebs generally form at or adjacent to the point of maximal WSS and are aligned with major flow structures. Wall pressure and OSI do not contribute to determining bleb location. The measurement of WSS using CFD models may potentially predict bleb formation and thus improve the assessment of rupture risk in unruptured aneurysms.

A follow up MRI-based geometry and computational fluid dynamics study of carotid bifurcation

Cardiovascular disease is the leading causes of death in the developed world. Wall shear stress (WSS) is associated with the initiation and progression of atherogenesis. This study combined the recent advances in MR imaging and computational fluid dynamics (CFD) and evaluated the patient-specific carotid bifurcation. The patient was followed up for 3 years. The geometry changes (tortuosity, curvature, ICA/CCA area ratios, central to the cross-sectional curvature, maximum stenosis) and the CFD factors (Velocity distribute, Wall Shear Stress (WSS) and Oscillatory Shear Index (OSI)) were compared at different time points.The carotid stenosis was a slight increase in the central to the cross-sectional curvature, and it was minor and variable curvature changes for carotid centerline. The OSI distribution presents ahigh-values in the same region where carotid stenosis and normal border, indicating complex flow and recirculation.The significant geometric changes observed during the follow-up may also cause significant changes in bifurcation hemodynamics.

Comparison between computational models of abdominal aortic aneurysm

Objective: To compare the differences in the hemodynamic parameters of abdominal aortic aneurysm (AAA) between fluid-structure interaction model (FSIM) and fluid-only model (FM), so as to discuss their application in the research of AAA. Methods: An idealized AAA model was created based on patient-specific AAA data. In FM, the flow, pressure and wall shear stress (WSS) were computed using finite volume method. In FSIM, an Arbitrary Lagrangian-Eulerian algorithm was used to solve the flow in a continuously deforming geometry. The hemodynamic parameters of both models were obtained for discussion. Results: Under the same inlet velocity, there were only two symmetrical vortexes in the AAA dilation area for FSIM. In contrast, four recirculation areas existed in FM; two were main vortexes and the other two were secondary flow, which were located between the main recirculation area and the arterial wall. Six local pressure concentrations occurred in the distal end of AAA and the recirculation area for FM. However, there were only two local pressure concentrations in FSIM. The vortex center of the recirculation area in FSIM was much more close to the distal end of AAA and the area was much larger because of AAA expansion. Four extreme values of WSS existed at the proximal of AAA, the point of boundary layer separation, the point of flow reattachment and the distal end of AAA, respectively, in both FM and FSIM. The maximum wall stress and the largest wall deformation were both located at the proximal and distal end of AAA. Conclusions: The number and center of the recirculation area for both models are different, while the change of vortex is closely associated with the AAA growth. The largest WSS of FSIM is 36% smaller than that of FM. Both the maximum wall stress and largest wall displacement shall increase with the outlet pressure increasing. FSIM needs to be considered for studying the relationship between AAA growth and shear stress.

MRI-based fatigue analysis predicts plaque vulnerability: A study comparing symptomatic and asymptomatic individuals

BACKGROUND: Rupture of atheromatous plaque in the carotid artery often leads to thrombosis and subsequent stroke. The mechanism of plaque rupture is not entirely clear but is thought to be a multi-factorial process involving thinning and weakening of the fibrous cap and biomechanical stress as the trigger leading to plaque rupture. As the cardiovascular system is a classic fatigue environment, the weakening of plaque leading to rupture may be a fatigue process, which is a symptomatically quiescent but potentially progressive failure process. In this study, we used a fatigue analysis based on in vivo magnetic resonance imaging (MRI) to investigate the rupture initiation location, crack propagation path and fatigue life within plaques of asymptomatic and symptomatic individuals. METHODS: Forty non-consecutive subjects (20 symptomatic and 20 asymptomatic) underwent high-resolution multi-sequence in vivo MRI of the carotid bifurcation. Fatigue analysis was performed based on the plaque geometry derived from in vivo MRI of the carotid artery at the point of maximum stenosis. Paris’ Law in fracture mechanics is adopted to determine the fatigue crack growth rate. Incremental crack propagation was dynamically simulated based on stress distributions. Plaque initiation location, crack propagation path and fatigue cycle of symptomatic and asymptomatic individuals were compared. RESULTS: Cracks were often found to begin at the lumen wall at areas of stress concentration. The preferred rupture direction was radial from the lumen center. The crack initially advanced slowly but accelerated as it developed, depending on plaque morphology. The fatigue cycles of symptomatic plaques were significantly less than those in the asymptomatic group (2.3 ± 0.9 vs 3.1 ± 0.7 (x106); p = 0.003). CONCLUSIONS: The number of cycles to rupture in symptomatic patients was higher than those predicted in asymptomatic patients by fatigue analysis, suggesting the possibility that plaques with a less fatigue life may be more prone to be symptomatic and rupture. If further validated by large-scale longitudinal studies, fatigue analysis based on high resolution in vivo MRI could potentially act as a useful tool for risk assessment of carotid atheroma.