Influence of anatomical dominance and hypertension on coronary conduit arterial and microcirculatory flow patterns: a multiscale modeling study

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.

Review of cardiac-coronary interaction and insights from mathematical modeling

Cardiac-coronary interaction is fundamental to the function of the heart. As one of the highest metabolic organs in the body, the cardiac oxygen demand is met by blood perfusion through the coronary vasculature. The coronary vasculature is largely embedded within the myocardial tissue which is continually contracting and hence squeezing the blood vessels. The myocardium-coronary vessel interaction is two-ways and complex. Here, we review the different types of cardiac-coronary interactions with a focus on insights gained from mathematical models. Specifically, we will consider the following: (1) myocardial-vessel mechanical interaction; (2) metabolic-flow interaction and regulation; (3) perfusion-contraction matching, and (4) chronic interactions between the myocardium and coronary vasculature. We also provide a discussion of the relevant experimental and clinical studies of different types of cardiac-coronary interactions. Finally, we highlight knowledge gaps, key challenges, and limitations of existing mathematical models along with future research directions to understand the unique myocardium-coronary coupling in the heart. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.

Mechanistic insights on age-related changes in heart-aorta-brain hemodynamic coupling using a pulse wave model of the entire circulatory system

Age-related changes in aortic biomechanics can impact the brain by reducing blood flow and increasing pulsatile energy transmission. Clinical studies have shown that impaired cardiac function in patients with heart failure is associated with cognitive impairment. Although previous studies have attempted to elucidate the complex relationship between age-associated aortic stiffening and pulsatility transmission to the cerebral network, they have not adequately addressed the effect of interactions between aortic stiffness and left ventricle (LV) contractility (neither on energy transmission nor on brain perfusion). In this study, we use a well-established and validated one-dimensional blood flow and pulse wave computational model of the circulatory system to address how age-related changes in cardiac function and vasculature affect the underlying mechanisms involved in the LV-aorta-brain hemodynamic coupling. Our results reveal how LV contractility affects pulsatile energy transmission to the brain, even with preserved cardiac output. Our model demonstrates the existence of an optimal heart rate (near the normal human heart rate) that minimizes pulsatile energy transmission to the brain at different contractility levels. Our findings further suggest that the reduction in cerebral blood flow at low levels of LV contractility is more prominent in the setting of age-related aortic stiffening. Maintaining optimal blood flow to the brain requires either an increase in contractility or an increase in heart rate. The former consistently leads to higher pulsatile power transmission, and the latter can either increase or decrease subsequent pulsatile power transmission to the brain.NEW & NOTEWORTHY We investigated the impact of major aging mechanisms of the arterial system and cardiac function on brain hemodynamics. Our findings suggest that aging has a significant impact on heart-aorta-brain coupling through changes in both arterial stiffening and left ventricle (LV) contractility. Understanding the underlying physical mechanisms involved here can potentially be a key step for developing more effective therapeutic strategies that can mitigate the contributions of abnormal LV-arterial coupling toward neurodegenerative diseases and dementia.

Mechanism of exercise intolerance in heart diseases predicted by a computer model of myocardial demand-supply feedback system

BACKGROUND AND OBJECTIVE

The myocardial demand-supply feedback system plays an important role in augmenting blood supply in response to exercise-induced increased myocardial demand. During this feedback process, the myocardium and coronary blood flow interact bidirectionally at many different levels.

METHODS

To investigate these interactions, a novel computational framework that considers the closed myocardial demand-supply feedback system was developed. In the framework coupling the systemic circulation of the left ventricle and coronary perfusion with regulation, myocardial work affects coronary perfusion via flow regulation mechanisms (e.g., metabolic regulation) and myocardial-vessel interactions, whereas coronary perfusion affects myocardial contractility in a closed feedback system. The framework was calibrated based on the measurements from healthy subjects under graded exercise conditions, and then was applied to simulate the effects of graded exercise on myocardial demand-supply under different physiological and pathological conditions.

RESULTS

We found that the framework can recapitulate key features found during exercise in clinical and animal studies. We showed that myocardial blood flow is increased but maximum hyperemia is reduced during exercise, which led to a reduction in coronary flow reserve. For coronary stenosis and myocardial inefficiency, the model predicts that an increase in heart rate is necessary to maintain the baseline cardiac output. Correspondingly, the resting coronary flow reserve is exhausted and the range of heart rate before exhaustion of coronary flow reserve is reduced. In the presence of metabolic regulation dysfunction, the model predicts that the metabolic vasodilator signal is higher at rest, saturates faster during exercise, and as a result, causes quicker exhaustion of coronary flow reserve.

CONCLUSIONS

Model predictions showed that the coronary flow reserve deteriorates faster during graded exercise, which in turn, suggests a decrease in exercise tolerance for patients with stenosis, myocardial inefficiency and metabolic flow regulation dysfunction. The findings in this study may have clinical implications in diagnosing cardiovascular diseases.

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.

A Closed-Loop Modeling Framework for Cardiac-to-Coronary Coupling

The mechanisms by which cardiac mechanics effect coronary perfusion (cardiac-to-coronary coupling) remain incompletely understood. Several coronary models have been proposed to deepen our understanding of coronary hemodynamics, but possibilities for in-depth studies on cardiac-to-coronary coupling are limited as mechanical properties like myocardial stress and strain are most often neglected. To overcome this limitation, a mathematical model of coronary mechanics and hemodynamics was implemented in the previously published multi-scale CircAdapt model of the closed-loop cardiovascular system. The coronary model consisted of a relatively simple one-dimensional network of the major conduit arteries and veins as well as a lumped parameter model with three transmural layers for the microcirculation. Intramyocardial pressure was assumed to arise from transmission of ventricular cavity pressure into the myocardial wall as well as myocardial stiffness, based on global pump mechanics and local myofiber mechanics. Model-predicted waveforms of global epicardial flow velocity, as well as of intramyocardial flow and diameter were qualitatively and quantitatively compared with reported data. Versatility of the model was demonstrated in a case study of aortic valve stenosis. The reference simulation correctly described the phasic pattern of coronary flow velocity, arterial flow impediment, and intramyocardial differences in coronary flow and diameter. Predicted retrograde flow during early systole in aortic valve stenosis was in agreement with measurements obtained in patients. In conclusion, we presented a powerful multi-scale modeling framework that enables realistic simulation of coronary mechanics and hemodynamics. This modeling framework can be used as a research platform for in-depth studies of cardiac-to-coronary coupling, enabling study of the effect of abnormal myocardial tissue properties on coronary hemodynamics.

Transmural Distribution of Coronary Perfusion and Myocardial Work Density Due to Alterations in Ventricular Loading, Geometry and Contractility

Myocardial supply changes to accommodate the variation of myocardial demand across the heart wall to maintain normal cardiac function. A computational framework that couples the systemic circulation of a left ventricular (LV) finite element model and coronary perfusion in a closed loop is developed to investigate the transmural distribution of the myocardial demand (work density) and supply (perfusion) ratio. Calibrated and validated against measurements of LV mechanics and coronary perfusion, the model is applied to investigate changes in the transmural distribution of passive coronary perfusion, myocardial work density, and their ratio in response to changes in LV contractility, preload, afterload, wall thickness, and cavity volume. The model predicts the following: (1) Total passive coronary flow varies from a minimum value at the endocardium to a maximum value at the epicardium transmurally that is consistent with the transmural distribution of IMP; (2) Total passive coronary flow at different transmural locations is increased with an increase in either contractility, afterload, or preload of the LV, whereas is reduced with an increase in wall thickness or cavity volume; (3) Myocardial work density at different transmural locations is increased transmurally with an increase in either contractility, afterload, preload or cavity volume of the LV, but is reduced with an increase in wall thickness; (4) Myocardial work density-perfusion mismatch ratio at different transmural locations is increased with an increase in contractility, preload, wall thickness or cavity volume of the LV, and the ratio is higher at the endocardium than the epicardium. These results suggest that an increase in either contractility, preload, wall thickness, or cavity volume of the LV can increase the vulnerability of the subendocardial region to ischemia.

Diagnostic performance of virtual fractional flow reserve derived from routine coronary angiography using segmentation free reduced order (1-dimensional) flow modelling

Introduction

Fractional flow reserve (FFR) improves assessment of the physiological significance of coronary lesions compared with conventional angiography. However, it is an invasive investigation. We tested the performance of a virtual FFR (1D-vFFR) using routine angiographic images and a rapidly performed reduced order computational model.

Methods

Quantitative coronary angiography (QCA) was performed in 102 with coronary lesions assessed by invasive FFR. A 1D-vFFR for each lesion was created using reduced order (one-dimensional) computational flow modelling derived from conventional angiographic images and patient specific estimates of coronary flow. The diagnostic accuracy of 1D-vFFR and QCA derived stenosis was compared against the gold standard of invasive FFR using area under the receiver operator characteristic curve (AUC).

Results

QCA revealed the mean coronary stenosis diameter was 44% ± 12% and lesion length 13 ± 7 mm. Following angiography calculation of the 1DvFFR took less than one minute. Coronary stenosis (QCA) had a significant but weak correlation with FFR (r = -0.2, p = 0.04) and poor diagnostic performance to identify lesions with FFR <0.80 (AUC 0.39, p = 0.09), (sensitivity - 58% and specificity - 26% at a QCA stenosis of 50%). In contrast, 1D-vFFR had a better correlation with FFR (r = 0.32, p = 0.01) and significantly better diagnostic performance (AUC 0.67, p = 0.007), (sensitivity - 92% and specificity - 29% at a 1D-vFFR of 0.7).

Conclusions

1D-vFFR improves the determination of functionally significant coronary lesions compared with conventional angiography without requiring a pressure-wire or hyperaemia induction. It is fast enough to influence immediate clinical decision-making but requires further clinical evaluation.

Cardiovascular deconditioning during long-term spaceflight through multiscale modeling

Human spaceflight has been fascinating man for centuries, representing the intangible need to explore the unknown, challenge new frontiers, advance technology, and push scientific boundaries further. A key area of importance is cardiovascular deconditioning, that is, the collection of hemodynamic changes-from blood volume shift and reduction to altered cardiac function-induced by sustained presence in microgravity. A thorough grasp of the 0G adjustment point per se is important from a physiological viewpoint and fundamental for astronauts' safety and physical capability on long spaceflights. However, hemodynamic details of cardiovascular deconditioning are incomplete, inconsistent, and poorly measured to date; thus a computational approach can be quite valuable. We present a validated 1D-0D multiscale model to study the cardiovascular response to long-term 0G spaceflight in comparison to the 1G supine reference condition. Cardiac work, oxygen consumption, and contractility indexes, as well as central mean and pulse pressures were reduced, augmenting the cardiac deconditioning scenario. Exercise tolerance of a spaceflight traveler was found to be comparable to an untrained person with a sedentary lifestyle. At the capillary-venous level significant waveform alterations were observed which can modify the regular perfusion and average nutrient supply at the cellular level. The present study suggests special attention should be paid to future long spaceflights which demand prompt physical capacity at the time of restoration of partial gravity (e.g., Moon/Mars landing). Since spaceflight deconditioning has features similar to accelerated aging understanding deconditioning mechanisms in microgravity are also relevant to the understanding of aging physiology on the Earth.

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.

A One-Dimensional Hemodynamic Model of the Coronary Arterial Tree

One-dimensional (1D) hemodynamic models of arteries have increasingly been applied to coronary circulation. In this study, we have adopted flow and pressure profiles in Olufsen's 1D structured tree as coronary boundary conditions, with terminals coupled to the dynamic pressure feedback resulting from the intra-myocardial stress because of ventricular contraction. We model a trifurcation structure of the example coronary tree as two adjacent bifurcations. The estimated results of blood pressure and flow rate from our simulation agree well with the clinical measurements and published data. Furthermore, the 1D model enables us to use wave intensity analysis to simulate blood flow in the developed coronary model. Six characteristic waves are observed in both left and right coronary flows, though the waves' magnitudes differ from each other. We study the effects of arterial wall stiffness on coronary blood flow in the left circumflex artery (LCX). Different diseased cases indicate that distinct pathological reactions of the cardiovascular system can be better distinguished through Wave Intensity analysis, which shows agreement with clinical observations. Finally, the feedback pressure in terminal vessels and measurement deviation are also investigated by changing parameters in the LCX. We find that larger feedback pressure increases the backward wave and decreases the forward one. Although simplified, this 1D model provides new insight into coronary hemodynamics in healthy and diseased conditions. We believe that this approach offers reference resources for studies on coronary circulation disease diagnosis, treatment and simulation.

Association between coronary dominance and acute inferior myocardial infarction: a matched, case-control study

Background

Previous studies have found a connection between left coronary artery dominance and worse prognoses in patient with acute coronary syndrome, which remains a predominant cause of morbidity and mortality globally. The aim of this study was to investigate whether coronary dominance is associated with the incidence of acute inferior myocardial infarction (MI).

Methods

Between January 2011 and November 2014, 265 patients with acute inferior MI and 530 age-matched and sex-matched controls were recruited for a case-control study in the Second Affiliated Hospital of Xi’an Jiaotong University in Xi’an, China. All participants underwent coronary angiography. The exclusion criteria included history of coronary artery bypass graft surgery, chronic or systemic diseases (including hepatic failure, kidney failure, hypothyroidism and Grave’s disease), ventricular fibrillation, and known allergy to iodinated contrast agent. Patients with left- or co-dominant anatomies were placed into the LD group and those with right-dominant anatomy were included in the RD group. The association of acute inferior MI and coronary dominant anatomy were assessed using multivariable conditional logistic regression, and to estimate the odds ratio (OR) and 95% confidence interval (95%CI).

Results

Distributions of right dominance were significantly different between the acute inferior MI group and control group (94.0% vs. 87.9%, P = 0.018). Univariable conditional logistic regression revealed that right dominance may be a risk factor for the incident acute inferior MI (OR: 2.137; 95% CI: 1.210–3.776; P = 0.009). After adjusting for baseline systolic blood pressure, heart rate, smoking status, diabetes mellitus, hypertension, hyperlipidaemia, and family history of coronary artery disease, results of multivariate conditional logistic regression showed that right dominance was associated with the incidence of acute inferior MI (OR: 2.396; 95% CI: 1.328–4.321; P = 0.004).

Conclusions

Right coronary dominance may play a disadvantageous role in the incidence of acute inferior MI. However, further studies are needed to verify our findings, especially with regard to the underlying mechanisms.

A multi-scale model of the coronary circulation applied to investigate transmural myocardial flow

Distribution of blood flow in myocardium is a key determinant of the localization and severity of myocardial ischemia under impaired coronary perfusion conditions. Previous studies have extensively demonstrated the transmural difference of ischemic vulnerability. However, it remains incompletely understood how transmural myocardial flow is regulated under in vivo conditions. In the present study, a computational model of the coronary circulation was developed to quantitatively evaluate the sensitivity of transmural flow distribution to various cardiovascular and hemodynamic factors. The model was further incorporated with the flow autoregulatory mechanism to simulate the regulation of myocardial flow in the presence of coronary artery stenosis. Numerical tests demonstrated that heart rate (HR), intramyocardial tissue pressure (P_im ), and coronary perfusion pressure (P_per ) were the major determinant factors for transmural flow distribution (evaluated by the subendocardial-to-subepicardial (endo/epi) flow ratio) and that the flow autoregulatory mechanism played an important compensatory role in preserving subendocardial perfusion against reduced P_per . Further analysis for HR variation-induced hemodynamic changes revealed that the rise in endo/epi flow ratio accompanying HR decrease was attributable not only to the prolongation of cardiac diastole relative to systole, but more predominantly to the fall in P_im . Moreover, it was found that P_im and P_per interfered with each other with respect to their influence on transmural flow distribution. These results demonstrate the interactive effects of various cardiovascular and hemodynamic factors on transmural myocardial flow, highlighting the importance of taking into account patient-specific conditions in the explanation of clinical observations.

Major influence of a 'smoke and mirrors' effect caused by wave reflection on early diastolic coronary arterial wave intensity

KEY POINTS

Coronary wave intensity analysis (WIA) is an emerging technique for assessing upstream and downstream influences on myocardial perfusion. It is thought that a dominant backward decompression wave (BDW_dia ) is generated by a distal suction effect, while early-diastolic forward decompression (FDW_dia ) and compression (FCW_dia ) waves originate in the aorta. We show that wave reflection also makes a substantial contribution to FDW_dia , FCW_dia and BDW_dia , as quantified by a novel method. In 18 sheep, wave reflection accounted for ∼70% of BDW_dia , whereas distal suction dominated in a computer model representing a hypertensive human. Non-linear addition/subtraction of mechanistically distinct waves (e.g. wave reflection and distal suction) obfuscates the true contribution of upstream and downstream forces on measured waves (the 'smoke and mirrors' effect). The mechanisms underlying coronary WIA are more complex than previously thought and the impact of wave reflection should be considered when interpreting clinical and experimental data.

ABSTRACT

Coronary arterial wave intensity analysis (WIA) is thought to provide clear insight into upstream and downstream forces on coronary flow, with a large early-diastolic surge in coronary flow accompanied by a prominent backward decompression wave (BDW_dia ), as well as a forward decompression wave (FDW_dia ) and forward compression wave (FCW_dia ). The BDW_dia is believed to arise from distal suction due to release of extravascular compression by relaxing myocardium, while FDW_dia and FCW_dia are thought to be transmitted from the aorta into the coronary arteries. Based on an established multi-scale computational model and high-fidelity measurements from the proximal circumflex artery (Cx) of 18 anaesthetized sheep, we present evidence that wave reflection has a major impact on each of these three waves, with a non-linear addition/subtraction of reflected waves obscuring the true influence of upstream and downstream forces through concealment and exaggeration, i.e. a 'smoke and mirrors' effect. We also describe methods, requiring additional measurement of aortic WIA, for unravelling the separate influences of wave reflection versus active upstream/downstream forces on coronary waves. Distal wave reflection accounted for ∼70% of the BDW_dia in sheep, but had a lesser influence (∼25%) in the computer model representing a hypertensive human. Negative reflection of the BDW_dia at the coronary-aortic junction attenuated the Cx FDW_dia (by ∼40% in sheep) and augmented Cx FCW_dia (∼5-fold), relative to the corresponding aortic waves. We conclude that wave reflection has a major influence on early-diastolic WIA, and thus needs to be considered when interpreting coronary WIA profiles.

In silico analysis of the anti-hypertensive drugs impact on myocardial oxygen balance

Hypertension is a very common pathology, and its clinical treatment largely relies on different drugs. Some of these drugs exhibit specific protective functions in addition to those resulting from blood pressure reduction. In this work, we study the impact of commonly used anti-hypertensive drugs (RAAS, [Formula: see text] and calcium channel blockers) on myocardial oxygen supply-consumption balance, which plays a crucial role in type 2 myocardial infarction. To this aim, 42 wash-out hypertensive patients were selected, a number of measured data were used to set a validated multi-scale cardiovascular model to subject-specific conditions, and the administration of different drugs was suitably simulated. Our results ascribe the well-known major cardioprotective efficiency of [Formula: see text] blockers compared to other drugs to a positive change of myocardial oxygen balance due to the concomitant: (1) reduction in aortic systolic, diastolic and pulse pressures, (2) decrease in left ventricular work, diastolic cavity pressure and oxygen consumption, (3) increase in coronary flow and (4) ejection efficiency improvement. RAAS blockers share several positive outcomes with [Formula: see text] blockers, although to a reduced extent. In contrast, calcium channel blockers seem to induce some potentially negative effects on the myocardial oxygen balance.