Slackness between vessel and myocardium is necessary for coronary flow reserve

A multi-scale computational model for the passive mechanical behavior of right ventricular myocardium

We have previously demonstrated the importance of myofiber-collagen mechanical interactions in modeling the passive mechanical behavior of right ventricle free wall (RVFW) myocardium. To gain deeper insights into these coupling mechanisms, we developed a high-fidelity, micro-anatomically realistic 3D finite element model of right ventricle free wall (RVFW) myocardium by combining high-resolution imaging and supercomputer-based simulations. We first developed a representative tissue element (RTE) model at the sub-tissue scale by specializing the hyperelastic anisotropic structurally-based constitutive relations for myofibers and ECM collagen, and equi-biaxial and non-equibiaxial loading conditions were simulated using the open-source software FEniCS to compute the effective stress-strain response of the RTE. To estimate the model parameters of the RTE model, we first fitted a 'top-down' biaxial stress-strain behavior with our previous structurally based (tissue-scale) model, informed by the measured myofiber and collagen fiber composition and orientation distributions. Next, we employed a multi-scale approach to determine the tissue-level (5 x 5 x 0.7 mm specimen size) RVFW biaxial behavior via 'bottom-up' homogenization of the fitted RTE model, recapitulating the histologically measured myofiber and collagen orientation to the biaxial mechanical data. Our homogenization approach successfully reproduced the tissue-level mechanical behavior of our previous studies in all biaxial deformation modes, suggesting that the 3D micro-anatomical arrangement of myofibers and ECM collagen is indeed a primary mechanism driving myofiber-collagen interactions.

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.

Overview of mathematical modeling of myocardial blood flow regulation

The oxygen consumption by the heart and its extraction from the coronary arterial blood are the highest among all organs. Any increase in oxygen demand due to a change in heart metabolic activity requires an increase in coronary blood flow. This functional requirement of adjustment of coronary blood flow is mediated by coronary flow regulation to meet the oxygen demand without any discomfort, even under strenuous exercise conditions. The goal of this article is to provide an overview of the theoretical and computational models of coronary flow regulation and to reveal insights into the functioning of a complex physiological system that affects the perfusion requirements of the myocardium. Models for three major control mechanisms of myogenic, flow, and metabolic control are presented. These explain how the flow regulation mechanisms operating over multiple spatial scales from the precapillaries to the large coronary arteries yield the myocardial perfusion characteristics of flow reserve, autoregulation, flow dispersion, and self-similarity. The review not only introduces concepts of coronary blood flow regulation but also presents state-of-the-art advances and their potential to impact the assessment of coronary microvascular dysfunction (CMD), cardiac-coronary coupling in metabolic diseases, and therapies for angina and heart failure. Experimentalists and modelers not trained in these models will have exposure through this review such that the nonintuitive and highly nonlinear behavior of coronary physiology can be understood from a different perspective. This survey highlights knowledge gaps, key challenges, future research directions, and novel paradigms in the modeling of coronary flow regulation.

Effects of myocardial function and systemic circulation on regional coronary perfusion

Cardiac-coronary interaction and the effects of its pathophysiological variations on spatial heterogeneity of coronary perfusion and myocardial work are still poorly understood. This hypothesis-generating study predicts spatial heterogeneities in both regional cardiac work and perfusion that offer a new paradigm on the vulnerability of the subendocardium to ischemia, particularly at the apex. We propose a mathematical and computational modeling framework to simulate the interaction of left ventricular mechanics, systemic circulation, and coronary microcirculation. The computational simulations revealed that the relaxation rate of the myocardium has a significant effect whereas the contractility has a marginal effect on both the magnitude and transmural distribution of coronary perfusion. The ratio of subendocardial to subepicardial perfusion density (Q_endo/Q_epi) changed by -12 to +6% from a baseline value of 1.16 when myocardial contractility was varied by +25 and -10%, respectively; Q_endo/Q_epi changed by 37% when sarcomere relaxation rate, b, was faster and increased by 10% from the baseline value. The model predicts axial differences in regional myocardial work and perfusion density across the wall thickness. Regional myofiber work done at the apex is 30-50% lower than at the center region, whereas perfusion density in the apex is lower by only 18% compared with the center. There are large axial differences in coronary flow and myocardial work at the subendocardial locations, with the highest differences located at the apex region. A mismatch exists between perfusion density and regional work done at the subendocardium. This mismatch is speculated to be compensated by coronary autoregulation.NEW & NOTEWORTHY We present a model of left ventricle perfusion based on an anatomically realistic coronary tree structure that includes its interaction with the systemic circulation. Left ventricular relaxation rate has a significant effect on the regional distribution of coronary flow and myocardial work.

Morphometric Reconstruction of Coronary Vasculature Incorporating Uniformity of Flow Dispersion

Experimental limitations in measurements of coronary flow in the beating heart have led to the development of in silico models of reconstructed coronary trees. Previous coronary reconstructions relied primarily on anatomical data, including statistical morphometry (e.g., diameters, length, connectivity, longitudinal position). Such reconstructions are non-unique, however, often leading to unrealistic predicted flow features. Thus, it is necessary to impose physiological flow constraints to ensure realistic tree reconstruction. Since a vessel flow depends on its diameter to fourth power, diameters are the logical candidates to guide vascular reconstructions to achieve realistic flows. Here, a diameter assignment method was developed where each vessel diameter was determined depending on its downstream tree size, aimed to reduce flow dispersion to within measured range. Since the coronary micro-vessels are responsible for a major portion of the flow resistance, the auto regulated coronary flow was analyzed in a morphometry-based reconstructed 400 vessel arterial microvascular sub-tree spanning vessel orders 1–6. Diameters in this subtree were re-assigned based on the flow criteria. The results revealed that diameter re-assignment, while adhering to measured morphometry, significantly reduced the flow dispersion to realistic levels while adhering to measured morphometry. The resulting network flow has longitudinal pressure distribution, flow fractal nature, and near-neighboring flow autocorrelation, which agree with measured coronary flow characteristics. Collectively, these results suggest that a realistic coronary tree reconstruction should impose not only morphometric data but also flow considerations. The work is of broad significance in providing a novel computational framework in the field of coronary microcirculation. It is essential for the study of coronary circulation by model simulation, based on a realistic network structure.

Role of Coronary Myogenic Response in Pressure-Flow Autoregulation in Swine: A Meta-Analysis With Coronary Flow Modeling

Myogenic responses (pressure-dependent contractions) of coronary arterioles play a role in autoregulation (relatively constant flow vs. pressure). Publications on myogenic reactivity in swine coronaries vary in caliber, analysis, and degree of responsiveness. Further, data on myogenic responses and autoregulation in swine have not been completely compiled, compared, and modeled. Thus, it has been difficult to understand these physiological phenomena. Our purpose was to: (a) analyze myogenic data with standard criteria; (b) assign results to diameter categories defined by morphometry; and (c) use our novel multiscale flow model to determine the extent to which ex vivo myogenic reactivity can explain autoregulation in vivo. When myogenic responses from the literature are an input for our model, the predicted coronary autoregulation approaches in vivo observations. More complete and appropriate data are now available to investigate the regulation of coronary blood flow in swine, a highly relevant model for human physiology and disease.

Integrative model of coronary flow in anatomically based vasculature under myogenic, shear, and metabolic regulation

Coronary blood flow is regulated to match the oxygen demand of myocytes in the heart wall. Flow regulation is essential to meet the wide range of cardiac workload. The blood flows through a complex coronary vasculature of elastic vessels having nonlinear wall properties, under transmural heterogeneous myocardial extravascular loading. To date, there is no fully integrative flow analysis that incorporates global and local passive and flow control determinants. Here, we provide an integrative model of coronary flow regulation that considers the realistic asymmetric morphology of the coronary network, the dynamic myocardial loading on the vessels embedded in it, and the combined effects of local myogenic effect, local shear regulation, and conducted metabolic control driven by venous O₂ saturation level. The model predicts autoregulation (approximately constant flow over a wide range of coronary perfusion pressures), reduced heterogeneity of regulated flow, and presence of flow reserve, in agreement with experimental observations. Furthermore, the model shows that the metabolic and myogenic regulations play a primary role, whereas shear has a secondary one. Regulation was found to have a significant effect on the flow except under extreme (high and low) inlet pressures and metabolic demand. Novel outcomes of the model are that cyclic myocardial loading on coronary vessels enhances the coronary flow reserve except under low inlet perfusion pressure, increases the pressure range of effective autoregulation, and reduces the network flow in the absence of metabolic regulation. Collectively, these findings demonstrate the utility of the present biophysical model, which can be used to unravel the underlying mechanisms of coronary physiopathology.

Cap buckling as a potential mechanism of atherosclerotic plaque vulnerability

Plaque rupture in atherosclerosis is the primary cause of potentially deadly coronary events, yet about 40% of ruptures occur away from the plaque cap shoulders and cannot be fully explained with the current biomechanical theories. Here, cap buckling is considered as a potential destabilizing factor which increases the propensity of the atherosclerotic plaque to rupture and which may also explain plaque failure away from the cap shoulders. To investigate this phenomenon, quasistatic 2D finite element simulations are performed, considering the salient geometrical and nonlinear material properties of diverse atherosclerotic plaques over the range of physiological loads. The numerical results indicate that buckling may displace the location of the peak von Mises stresses in the deflected caps. Plaque buckling, together with its deleterious effects is further observed experimentally in plaque caps using a physical model of deformable mock coronary arteries with fibroatheroma. Moreover, an analytical approach combining quasistatic equilibrium equations with the Navier-Bresse formulas is used to demonstrate the buckling potential of a simplified arched slender cap under intraluminal pressure and supported by foundations. This analysis shows that plaque caps - calcified, fibrotic or cellular - may buckle in specific undulated shapes once submitted to critical loads. Finally, a preliminary analysis of intravascular ultrasonography recordings of patients with atherosclerotic coronary arteries corroborates the numerical, experimental and theoretical findings and shows that various plaque caps buckle in vivo. By displacing the sites of high stresses in the plaque cap, buckling may explain the atherosclerotic plaque cap rupture at various locations, including cap shoulders.