Integrated modeling and simulation of recruitment of myocardial perfusion and oxygen delivery in exercise.

While exercise-mediated vasoregulation in the myocardium is understood to be governed by autonomic, myogenic, and metabolic-mediated mechanisms, we do not yet understand the spatial heterogeneity of vasodilation or its effects on microvascular flow patterns and oxygen delivery. This study uses a simulation and modeling approach to explore the mechanisms underlying the recruitment of myocardial perfusion and oxygen delivery in exercise. The simulation approach integrates model components representing: whole-body cardiovascular hemodynamics, cardiac mechanics and myocardial work; myocardial perfusion; and myocardial oxygen transport. Integrating these systems together, model simulations reveal: (1.) To match expected flow and transmural flow ratios at increasing levels of exercise, a greater degree of vasodilation must occur in the subendocardium compared to the subepicardium. (2.) Oxygen extraction and venous oxygenation are predicted to substantially decrease with increasing exercise level preferentially in the subendocardium, suggesting that an oxygen-dependent error signal driving metabolic mediated recruitment of flow would be operative only in the subendocardium. (3.) Under baseline physiological conditions approximately 4% of the oxygen delivered to the subendocardium may be supplied via retrograde flow from coronary veins.

Unraveling the Gordian knot of coronary pressure-flow autoregulation.

The coronary circulation has the inherent ability to maintain myocardial perfusion constant over a wide range of perfusion pressures. The phenomenon of pressure-flow autoregulation is crucial in response to flow-limiting atherosclerotic lesions which diminish coronary driving pressure and increase risk of myocardial ischemia and infarction. Despite well over half a century of devoted research, understanding of the mechanisms responsible for autoregulation remains one of the most fundamental and contested questions in the field today. The purpose of this review is to highlight current knowledge regarding the complex interrelationship between the pathways and mechanisms proposed to dictate the degree of coronary pressure-flow autoregulation. Our group recently likened the intertwined nature of the essential determinants of coronary flow control to the symbolically unsolvable "Gordian knot". To further efforts to unravel the autoregulatory "knot", we consider recent challenges to the local metabolic and myogenic hypotheses and the complicated dynamic structural and functional heterogeneity unique to the heart and coronary circulation. Additional consideration is given to interrogation of putative mediators, role of K+ and Ca2+ channels, and recent insights from computational modeling studies. Improved understanding of how specific vasoactive mediators, pathways, and underlying disease states influence coronary pressure-flow relations stands to significantly reduce morbidity and mortality for what remains the leading cause of death worldwide.

Oxygen-sensing pathways below autoregulatory threshold act to sustain myocardial oxygen delivery during reductions in perfusion pressure.

The coronary circulation has an innate ability to maintain constant blood flow over a wide range of perfusion pressures. However, the mechanisms responsible for coronary autoregulation remain a fundamental and highly contested question. This study interrogated the local metabolic hypothesis of autoregulation by testing the hypothesis that hypoxemia-induced exaggeration of the metabolic error signal improves the autoregulatory response. Experiments were performed on open-chest anesthetized swine during stepwise changes in coronary perfusion pressure (CPP) from 140 to 40 mmHg under normoxic (n = 15) and hypoxemic (n = 8) conditions, in the absence and presence of dobutamine-induced increases in myocardial oxygen consumption (MVO2) (n = 5-7). Hypoxemia (PaO2 < 40 mmHg) decreased coronary venous PO2 (CvPO2) ~ 30% (P < 0.001) and increased coronary blood flow ~ 100% (P < 0.001), sufficient to maintain myocardial oxygen delivery (P = 0.14) over a wide range of CPPs. Autoregulatory responsiveness during hypoxemia-induced reductions in CvPO2 were associated with increases of autoregulatory gain (Gc; P = 0.033) but not slope (P = 0.585) over a CPP range of 120 to 60 mmHg. Preservation of autoregulatory Gc (P = 0.069) and slope (P = 0.264) was observed during dobutamine administration ± hypoxemia. Reductions in coronary resistance in response to decreases in CPP predominantly occurred below CvPO2 values of ~ 25 mmHg, irrespective of underlying vasomotor reserve. These findings support the presence of an autoregulatory threshold under which oxygen-sensing pathway(s) act to preserve sufficient myocardial oxygen delivery as CPP is reduced during increases in MVO2 and/or reductions in arterial oxygen content.

Multiscale model of the physiological control of myocardial perfusion to delineate putative metabolic feedback mechanisms.

Coronary blood flow is tightly regulated to ensure that myocardial oxygen delivery meets local metabolic demand via the concurrent action of myogenic, neural and metabolic mechanisms. Although several competing hypotheses exist, the specific nature of the local metabolic mechanism(s) remains poorly defined. To gain insights into the viability of putative metabolic feedback mechanisms and into the co-ordinated action of parallel regulatory mechanisms, we applied a multiscale modelling framework to analyse experimental data on coronary pressure, flow and myocardial oxygen delivery in the porcine heart in vivo. The modelling framework integrates a previously established lumped-parameter model of myocardial perfusion used to account for transmural haemodynamic variations and a simple vessel mechanics model used to simulate the vascular tone in each of three myocardial layers. Vascular tone in the resistance vessel mechanics model is governed by input stimuli from the myogenic, metabolic and autonomic control mechanisms. Seven competing formulations of the metabolic feedback mechanism are implemented in the modelling framework, and associated model simulations are compared with experimental data on coronary pressures and flows under a range of experimental conditions designed to interrogate the governing control mechanisms. Analysis identifies a maximally probable metabolic mechanism among the seven tested models, in which production of a metabolic signalling factor is proportional to myocardial oxygen consumption and delivery is proportional to flow. Finally, the identified model is validated based on comparisons of simulations with data on the myocardial perfusion response to conscious exercise that were not used for model identification. KEY POINTS: Although several competing hypotheses exist, we lack knowledge of specific nature of the metabolic mechanism(s) governing regional myocardial perfusion. Moreover, we lack an understanding of how parallel myogenic, adrenergic/autonomic and metabolic mechanisms work together to regulatory oxygen delivery in the beating heart. We have developed a multiscale modelling framework to test competing hypotheses against experimental data on coronary pressure, flow and myocardial oxygen delivery in the porcine heart in vivo. The analysis identifies a maximally probable metabolic mechanism among seven tested models, in which the production of a metabolic signalling factor is proportional to myocardial oxygen consumption and delivery is proportional to flow.

CRIMSON: An open-source software framework for cardiovascular integrated modelling and simulation.

In this work, we describe the CRIMSON (CardiovasculaR Integrated Modelling and SimulatiON) software environment. CRIMSON provides a powerful, customizable and user-friendly system for performing three-dimensional and reduced-order computational haemodynamics studies via a pipeline which involves: 1) segmenting vascular structures from medical images; 2) constructing analytic arterial and venous geometric models; 3) performing finite element mesh generation; 4) designing, and 5) applying boundary conditions; 6) running incompressible Navier-Stokes simulations of blood flow with fluid-structure interaction capabilities; and 7) post-processing and visualizing the results, including velocity, pressure and wall shear stress fields. A key aim of CRIMSON is to create a software environment that makes powerful computational haemodynamics tools accessible to a wide audience, including clinicians and students, both within our research laboratories and throughout the community. The overall philosophy is to leverage best-in-class open source standards for medical image processing, parallel flow computation, geometric solid modelling, data assimilation, and mesh generation. It is actively used by researchers in Europe, North and South America, Asia, and Australia. It has been applied to numerous clinical problems; we illustrate applications of CRIMSON to real-world problems using examples ranging from pre-operative surgical planning to medical device design optimization.

Noninvasive quantification of cerebrovascular pressure changes using 4D Flow MRI.

PURPOSE: Hemodynamic alterations are indicative of cerebrovascular disease. However, the narrow and tortuous cerebrovasculature complicates image-based assessment, especially when quantifying relative pressure. Here, we present a systematic evaluation of image-based cerebrovascular relative pressure mapping, investigating the accuracy of the routinely used reduced Bernoulli (RB), the extended unsteady Bernoulli (UB), and the full-field virtual work-energy relative pressure ( ν WERP) method. METHODS: Patient-specific in silico models were used to generate synthetic cerebrovascular 4D Flow MRI, with RB, UB, and ν WERP performance quantified as a function of spatiotemporal sampling and image noise. Cerebrovascular relative pressures were also derived in 4D Flow MRI from healthy volunteers ( n=8 ), acquired at two spatial resolutions (dx = 1.1 and 0.8 mm). RESULTS: The in silico analysis indicate that accurate relative pressure estimations are inherently coupled to spatial sampling: at dx = 1.0 mm high errors are reported for all methods; at dx = 0.5 mm ν WERP recovers relative pressures at a mean error of 0.02 ± 0.25 mm Hg, while errors remain higher for RB and UB (mean error of -2.18 ± 1.91 and -2.18 ± 1.87 mm Hg, respectively). The dependence on spatial sampling is also indicated in vivo, albeit with higher correlative dependence between resolutions using ν WERP (k = 0.64, R2 = 0.81 for dx = 1.1 vs. 0.8 mm) than with RB or UB (k = 0.04, R2 = 0.03, and k = 0.07, R2 = 0.07, respectively). CONCLUSION: Image-based full-field methods such as ν WERP enable cerebrovascular relative pressure mapping; however, accuracy is directly dependent on utilized spatial resolution.

Endovascular ascending aortic repair in type A dissection: A systematic review.

PURPOSE: Up to 10% of acute type A aortic dissection (TAAD) patients are deemed unfit for open surgical repair, exposing these patients to high mortality rates. In recent years, thoracic endovascular aortic repair has proven to be a promising alternative treatment modality in specific cases. This study presents a comprehensive overview of the current state of catheter-based interventions in the setting of primary TAAD. METHODS: A literature search was conducted, using MEDLINE and PubMed databases according to PRISMA guidelines, updated until January 2020. Articles were selected if they reported on the endovascular repair of DeBakey Type I and II aortic dissections. The exclusion criteria were retrograde type A dissection, hybrid procedures, and combined outcome reporting of mixed aortic pathologies (e.g., pseudoaneurysm and intramural hematoma). RESULTS: A total of 31 articles, out of which 19 were case reports and 12 case series, describing a total of 92 patients, were included. The median follow-up was 6 months for case reports and the average follow-up was 14 months for case series. Overall technical success was 95.6% and 30-day mortality of 9%. Stroke and early endoleak rates were 6% and 18%, respectively. Reintervention was required in 14 patients (15%). CONCLUSION: This review not only demonstrates that endovascular repair in the setting of isolated TAAD is feasible with acceptable outcomes at short-term follow-up, but also underlines a lack of mid-late outcomes and reporting consistency. Studies with longer follow-up and careful consideration of patient selection are required before endovascular interventions can be widely introduced.

Practical considerations for territorial perfusion mapping in the cerebral circulation using super-selective pseudo-continuous arterial spin labeling.

PURPOSE: This paper discusses several challenges faced by super-selective pseudo-continuous arterial spin labeling, which is used to quantify territorial perfusion in the cerebral circulation. The effects of off-resonance, pulsatility, vessel movement, and label rotation scheme are investigated, and methods to maximize labeling efficiency and overall image quality are evaluated. A strategy to calculate the territorial perfusion fractions of individual vessels is proposed. METHODS: The effects of off-resonance, label rotation scheme, and vessel movement on labeling efficiency were simulated. Two off-resonance compensation strategies (multiphase prescan, field map), cardiac triggering, and vessel movement were studied in vivo in a group of 10 subjects. Subsequently, a territorial perfusion fraction map was acquired in 2 subjects based on the mean vessel labeling efficiency. RESULTS: Multiphase calibration provided the highest labeling efficiency (P = .002) followed by the field map compensation (P = .037) compared with the uncompensated acquisition. Cardiac triggering resulted in a qualitative improvement of the image and an increase in signal contrast between the perfusion territory and the surrounding tissue (P = .010) but failed to show a significant change in temporal and spatial SNR. The constant clockwise label rotation scheme yielded the highest labeling efficiency. Significant vessel movement (>2 mm according to simulations) was observed in 50% of subjects. The measured territorial perfusion fractions showed good agreement with anatomical data. CONCLUSION: Optimized labeling efficiency resulted in increased image quality and accuracy of territorial perfusion fraction maps. Labeling efficiency depends critically on off-resonance calibration, cardiac triggering, optimal label rotation scheme, and vessel location tracking.

Computational analysis of renal artery flow characteristics by modeling aortoplasty and aortic bypass interventions for abdominal aortic coarctation.

OBJECTIVE: Suprarenal abdominal aortic coarctation (SAAC) alters flow and pressure patterns to the kidneys and is often associated with severe angiotensin-mediated hypertension refractory to drug therapy. SAAC is most often treated by a thoracoabdominal bypass (TAB) or patch aortoplasty (PA). It is currently unclear what effect these interventions have on renal flow and pressure waveforms. This study, using retrospective data from a patient with SAAC subjected to a TAB, undertook computational modeling to analyze aortorenal blood flow preoperatively as well as postoperatively after a variety of TAB and PA interventions. METHODS: Patient-specific anatomic models were constructed from preoperative computed tomography angiograms of a 9-year-old child with an isolated SAAC. Fluid-structure interaction (FSI) simulations of hemodynamics were performed to analyze preoperative renal flow and pressure waveforms. A parametric study was then performed to examine the hemodynamic impact of different bypass diameters and patch oversizing. RESULTS: Preoperative FSI results documented diastole-dominated renal perfusion with considerable high-frequency disturbances in blood flow and pressure. The postoperative TAB right and left kidney volumes increased by 58% and 79%, respectively, reflecting the increased renal artery blood flows calculated by the FSI analysis. Postoperative increases in systolic flow accompanied decreases in high-frequency disturbances, aortic pressure, and collateral flow after all surgical interventions. In general, lesser degrees of high-frequency disturbances followed PA interventions. High-frequency disturbances were eliminated with the 0% PA in contrast to the 30% and 50% PA oversizing and TAB interventions, in which these flow disturbances remained. CONCLUSIONS: Both TAB and PA dramatically improved renal artery flow and pressure waveforms, although disturbed renal waveforms remained in many of the surgical scenarios. Importantly, only the 0% PA oversizing scenario eliminated all high-frequency disturbances, resulting in nearly normal aortorenal blood flow. The study also establishes the relevance of patient-specific computational modeling in planning interventions for the midaortic syndrome.

Comparative Analysis of Porcine and Human Thoracic Aortic Stiffness.

OBJECTIVES: To compare porcine and human thoracic aortic stiffness using the available literature. METHODS: The available literature was searched for studies reporting data on porcine or human thoracic aortic mechanical behaviour. A four fibre constitutive model was used to transform the data from included studies. Thus, equi-biaxial stress stretch curves were generated to calculate circumferential and longitudinal aortic stiffness. Analysis was performed separately for the ascending and descending thoracic aorta. Data on human aortic stiffness were divided by age <60 or ≥60 years. Porcine and human aortic stiffness were compared. RESULTS: Eleven studies were included, six reported on young porcine aortas, four on human aortas of various ages, and one reported on both. In the ascending aorta, circumferential and longitudinal stiffness were 0.42±0.08 MPa and 0.37±0.06 MPa for porcine aortas (4-9 months) versus 0.55±0.15 MPa and 0.45±0.08 MPa for humans <60 years, and 1.02±0.59 MPa and 1.03±0.54 MPa for humans ≥60 years. In the descending aorta, circumferential and longitudinal stiffness were 0.46±0.03 MPa and 0.44±0.01 MPa for porcine aortas (4-10 months) versus 1.04±0.70 MPa and 1.24±0.76 MPa for humans <60 years, and 3.15±3.31 MPa and 1.17±0.31 MPa for humans ≥60 years. CONCLUSIONS: The stiffness of young porcine aortic tissue shows good correspondence with human tissue aged <60 years, especially in the ascending aorta. Young porcine aortic tissue is less stiff than human aortic tissue aged ≥60 years.

Functional assessment of thoracic aortic aneurysms - the future of risk prediction?

Introduction: Treatment guidelines for the thoracic aorta concentrate on size, yet acute aortic dissection or rupture can occur when aortic size is below intervention criteria. Functional imaging and computational techniques are a means of assessing haemodynamic parameters involved in aortic pathology. Sources of data: Original articles, reviews, international guidelines. Areas of agreement: Computational fluid dynamics and 4D flow MRI allow non-invasive assessment of blood flow parameters and aortic wall biomechanics. Areas of controversy: Aortic valve morphology (particularly bicuspid aortic valve) is associated with aneurysm of the ascending aorta, although the exact mechanism of aneurysm formation is not yet established. Growing points: Haemodynamic assessment of the thoracic aorta has highlighted parameters which are linked with both clinical outcome and protein changes in the aortic wall. Wall shear stress, flow displacement and helicity are elevated in patients with bicuspid aortic valve, particularly at locations of aneurysm formation. Areas timely for developing research: With further validation, functional assessment of the aorta may help identify patients at risk of aortic complications, and introduce new haemodynamic indices into management guidelines.

Aortic length measurements for pulse wave velocity calculation: manual 2D vs automated 3D centreline extraction.

BACKGROUND: Pulse wave velocity (PWV) is a biomarker for the intrinsic stiffness of the aortic wall, and has been shown to be predictive for cardiovascular events. It can be assessed using cardiovascular magnetic resonance (CMR) from the delay between phase-contrast flow waveforms at two or more locations in the aorta, and the distance on CMR images between those locations. This study aimed to investigate the impact of different distance measurement methods on PWV. We present and evaluate an algorithm for automated centreline tracking in 3D images, and compare PWV calculations using distances derived from 3D images to those obtained from a conventional 2D oblique-sagittal image of the aorta. METHODS: We included 35 patients from a twin cohort, and 20 post-coarctation repair patients. Phase-contrast flow was acquired in the ascending, descending and diaphragmatic aorta. A 3D centreline tracking algorithm is presented and evaluated on a subset of 30 subjects, on three CMR sequences: balanced steady-state free precession (SSFP), black-blood double inversion recovery turbo spin echo, and contrast-enhanced CMR angiography. Aortic lengths are subsequently compared between measurements from a 2D oblique-sagittal plane, and a 3D geometry. RESULTS: The error in length of automated 3D centreline tracking compared with manual annotations ranged from 2.4 [1.8-4.3] mm (mean [IQR], black-blood) to 6.4 [4.7-8.9] mm (SSFP). The impact on PWV was below 0.5m/s (<5%). Differences between 2D and 3D centreline length were significant for the majority of our experiments (p < 0.05). Individual differences in PWV were larger than 0.5m/s in 15% of all cases (thoracic aorta) and 37% when studying the aortic arch only. Finally, the difference between end-diastolic and end-systolic 2D centreline lengths was statistically significant (p < 0.01), but resulted in small differences in PWV (0.08 [0.04 - 0.10]m/s). CONCLUSIONS: Automatic aortic centreline tracking in three commonly used CMR sequences is possible with good accuracy. The 3D length obtained from such sequences can differ considerably from lengths obtained from a 2D oblique-sagittal plane, depending on aortic curvature, adequate planning of the oblique-sagittal plane, and patient motion between acquisitions. For accurate PWV measurements we recommend using 3D centrelines.

A mathematical model of coronary blood flow control: simulation of patient-specific three-dimensional hemodynamics during exercise.

This work presents a mathematical model of the metabolic feedback and adrenergic feedforward control of coronary blood flow that occur during variations in the cardiac workload. It is based on the physiological observations that coronary blood flow closely follows myocardial oxygen demand, that myocardial oxygen debts are repaid, and that control oscillations occur when the system is perturbed and so are phenomenological in nature. Using clinical data, we demonstrate that the model can provide patient-specific estimates of coronary blood flow changes between rest and exercise, requiring only the patient's heart rate and peak aortic pressure as input. The model can be used in zero-dimensional lumped parameter network studies or as a boundary condition for three-dimensional multidomain Navier-Stokes blood flow simulations. For the first time, this model provides feedback control of the coronary vascular resistance, which can be used to enhance the physiological accuracy of any hemodynamic simulation, which includes both a heart model and coronary arteries. This has particular relevance to patient-specific simulation for which heart rate and aortic pressure recordings are available. In addition to providing a simulation tool, under our assumptions, the derivation of our model shows that ß-feedforward control of the coronary microvascular resistance is a mathematical necessity and that the metabolic feedback control must be dependent on two error signals: the historical myocardial oxygen debt, and the instantaneous myocardial oxygen deficit.

Update in the management of type B aortic dissection.

Stanford type B aortic dissection (TBAD) is a life-threatening aortic disease. The initial management goal is to prevent aortic rupture, propagation of the dissection, and symptoms by reducing the heart rate and blood pressure. Uncomplicated TBAD patients require prompt medical management to prevent aortic dilatation or rupture during subsequent follow-up. Complicated TBAD patients require immediate invasive management to prevent death or injury caused by rupture or malperfusion. Recent developments in diagnosis and management have reduced mortality related to TBAD considerably. In particular, the introduction of thoracic stent-grafts has shifted the management from surgical to endovascular repair, contributing to a fourfold increase in early survival in complicated TBAD. Furthermore, endovascular repair is now considered in some uncomplicated TBAD patients in addition to optimal medical therapy. For more challenging aortic dissection patients with involvement of the aortic arch, hybrid approaches, combining open and endovascular repair, have had promising results. Regardless of the chosen management strategy, strict antihypertensive control should be administered to all TBAD patients in addition to close imaging surveillance. Future developments in stent-graft design, medical therapy, surgical and hybrid techniques, imaging, and genetic screening may improve the outcomes of TBAD patients even further. We present a comprehensive review of the recommended management strategy based on current evidence in the literature.

Biomechanical Changes After Thoracic Endovascular Aortic Repair in Type B Dissection: A Systematic Review.

Thoracic endovascular aortic repair (TEVAR) has evolved into an established treatment option for type B aortic dissection (TBAD) since it was first introduced 2 decades ago. Morbidity and mortality have decreased due to the minimally invasive character of TEVAR, with adequate stabilization of the dissection, restoration of true lumen perfusion, and subsequent positive aortic remodeling. However, several studies have reported severe setbacks of this technique. Indeed, little is known about the biomechanical behavior of implanted thoracic stent-grafts and the impact on the vascular system. This study sought to systematically review the performance and behavior of implanted thoracic stent-grafts and related biomechanical aortic changes in TBAD patients in order to update current knowledge and future perspectives.

A comparative study of altered hemodynamics in iliac vein compression syndrome.

Introduction: Iliac vein compression syndrome (IVCS) is present in over 20% of the population and is associated with left leg pain, swelling, and thrombosis. IVCS symptoms are thought to be induced by altered pelvic hemodynamics, however, there currently exists a knowledge gap on the hemodynamic differences between IVCS and healthy patients. To elucidate those differences, we carried out a patient-specific, computational modeling comparative study. Methods: Computed tomography and ultrasound velocity and area data were used to build and validate computational models for a cohort of IVCS (N = 4, Subject group) and control (N = 4, Control group) patients. Flow, cross-sectional area, and shear rate were compared between the right common iliac vein (RCIV) and left common iliac vein (LCIV) for each group and between the Subject and Control groups for the same vessel. Results: For the IVCS patients, LCIV mean shear rate was higher than RCIV mean shear rate (550 ± 103 s-1 vs. 113 ± 48 s-1, p = 0.0009). Furthermore, LCIV mean shear rate was higher in the Subject group than in the Control group (550 ± 103 s-1 vs. 75 ± 37 s-1, p = 0.0001). Lastly, the LCIV/RCIV shear rate ratio was 4.6 times greater in the Subject group than in the Control group (6.56 ± 0.9 vs. 1.43 ± 0.6, p = 0.00008). Discussion: Our analyses revealed that IVCS patients have elevated shear rates which may explain a higher thrombosis risk and suggest that their thrombus initiation process may share aspects of arterial thrombosis. We have identified hemodynamic metrics that revealed profound differences between IVCS patients and Controls, and between RCIV and LCIV in the IVCS patients. Based on these metrics, we propose that non-invasive measurement of shear rate may aid with stratification of patients with moderate compression in which treatment is highly variable. More investigation is needed to assess the prognostic value of shear rate and shear rate ratio as clinical metrics and to understand the mechanisms of thrombus formation in IVCS patients.