AJP - Heart Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 278: H1490-H1499, 2000;
0363-6135/00 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (12)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Cornelissen, A. J. M.
Right arrow Articles by Spaan, J. A. E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Cornelissen, A. J. M.
Right arrow Articles by Spaan, J. A. E.
Vol. 278, Issue 5, H1490-H1499, May 2000

Myogenic reactivity and resistance distribution in the coronary arterial tree: a model study

Annemiek J. M. Cornelissen1, Jenny Dankelman1, Ed VanBavel2, Henk G. Stassen1, and Jos A. E. Spaan2

1 Faculty of Design, Engineering and Production, Mechanical Engineering and Marine Technology, Man Machine Systems and Control Group, Delft University of Technology, 2628 CD Delft; and 2 Department of Medical Physics, Cardiovascular Research Institute Amsterdam, Academic Medical Center, University of Amsterdam, 1100 DE Amsterdam, The Netherlands


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The objectives of this study were to evaluate the myogenic behavior of blood vessels and their interaction within the coronary arterial tree and to evaluate the possible role of the myogenic response in autoregulation. The model consists of 10 compartments in series, each representing a class of vessel sizes. Diameter and resistance in each class are determined by their value at full dilation (dp, Rp) and by the myogenic response. Three distributions of Rp and three distributions of myogenic strength, Mi (slope of pressure-diameter curve, range -0.05 to -0.4%/mmHg) were evaluated (9 cases). It was found that larger vessels attenuate the myogenic activity of smaller vessels and that myogenic responsiveness is sufficient to achieve autoregulation. When Mi has a maximum in vessels of 84 µm, the maximum effect of perfusion pressure on active diameter occurs in vessels between 123 and 181 µm, depending on the distribution of Rp. Distribution of resistance and control mechanisms in the coronary arterial tree are important for interpretation of individual vessel responses as observed in vivo.

coronary circulation; myogenic response; autoregulation; mathematical model


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

MYOGENIC CONSTRICTION of the smooth muscle cells in the arterial wall due to a rise in pressure is thought to be one of the mechanisms involved in coronary autoregulation, i.e., a relative constant flow despite changes in arterial pressure. Myogenic responsiveness in isolated resistance arteries is well documented (6, 16-20, 23, 26). However, the position of a resistance vessel in the arterial tree determines how effective the response can be. As Johnson (9, 10) argued in 1980, myogenic constriction of proximal arterioles due to arterial pressure elevation may limit the increase in pressure that is transmitted to the more distal vessels.

Studies in isolated vessels generally demonstrate heterogeneous myogenic responsiveness with the stronger myogenic reactivity in the smaller arterioles (4, 22). For example, in isolated porcine coronary arteries and arterioles with active diameters ranging from 37 to 180 µm, arterioles of ~62 µm were shown to exhibit the strongest myogenic reactivity (19). This agrees well with the observations of only a moderate responsiveness in porcine arteries of 180 µm (26) and an absence of myogenic responsiveness in porcine arteries of 223 µm (23). Therefore, on the basis of isolated vessel studies, one would expect that the smaller arterioles form the primary site for myogenic flow control.

Several studies report on the diameter changes of small arteries with diameters varying between 50 and 400 µm in the epicardial vascular bed in vivo, supporting the concept of diameter-dependent myogenic response (2, 12). However, since all these vessels are part of a network, it is not clear how the response of a single vessel under observation is determined by the other vessels in the tree.

The first aim of this study was to evaluate myogenic behavior of single coronary segments within a network and their interaction. The second aim was to evaluate the possible role of myogenic tone in autoregulation. We used a simple theoretical model with 10 resistance compartments in series. This is a strong simplification, since it has been shown that heterogeneity of branching has to be taken into account to achieve a good prediction of pressure distribution as a function of the vessel diameter (27). However, the complexities of such a heterogeneous model obscure the role of the distribution of resistance with respect to our specific aims. In our compartmental model, we evaluated three different resistance distributions, which covered the range of pressure distributions that can be found in literature, either as measured (3) or as estimated from the coronary branching structure (13, 14, 27).


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The model consists of 10 resistance compartments in series. The first nine compartments are applied to mimic the resistance distribution of the arterial tree. Nine compartments are used because this is approximately the number of Strahler orders of the vessels in the coronary arterial tree covering diameters ranging from 10 to 500 µm (15, 27). The resistance of capillaries and venules are assumed to be constant and are lumped in the 10th compartment. Passive and active pressure-diameter relations of vessels with certain diameter determine the resistance behaviors of each of the first nine compartments. Therefore, we assume that each of these compartments, i, represent Ni identical vessels in parallel. Hence, according to Poiseuille, the resistance can be written as
<FENCE><AR><R><C><IT>R</IT><SUB><IT>i</IT></SUB> = <FR><NU>128&eegr;<SUB><IT>i</IT></SUB><IT>L</IT><SUB><IT>i</IT></SUB></NU><DE><IT>N</IT><SUB><IT>i</IT></SUB>&pgr;(<IT>d</IT><SUB><IT>i</IT></SUB>)<SUP>4</SUP></DE></FR></C></R><R><C><IT>K</IT><SUB><IT>i</IT></SUB> = <FR><NU>&eegr;<SUB><IT>i</IT></SUB><IT>L</IT><SUB><IT>i</IT></SUB></NU><DE><IT>N</IT><SUB><IT>i</IT></SUB></DE></FR></C></R></AR></FENCE> ⇒ <IT>R</IT><SUB><IT>i</IT></SUB> = <FR><NU>128<IT>K</IT><SUB><IT>i</IT></SUB></NU><DE>&pgr;(<IT>d</IT><SUB><IT>i</IT></SUB>)<SUP>4</SUP></DE></FR>  <IT>i</IT> = 1, 2, … , 9 (1)
where eta i is the viscosity constant; Li is the length of the vessel; Ni is the number of parallel vessels in the compartment; and di is the diameter of the vessel. Since Ki includes the length and number of vessels as well as the viscosity constant in each compartment, it represents anatomical and rheological data. We assume Ki in each compartment independent of pressure and independent of diameter variations induced by myogenic response or pressure variations directly. The Fahraeus-Lindqvist effect is taken into account insofar as the viscosity in the reference condition is considered. However, possible viscosity changes due to this effect and induced by diameter changes are ignored. Local pressure Pi is taken as the mean of the inlet and the outlet pressure of the compartment, Pin,i and Pin,i+1, respectively.
P<SUB><IT>i</IT></SUB> = <FR><NU>P<SUB>in,<IT>i</IT></SUB> + P<SUB>in,<IT>i</IT>+1</SUB></NU><DE>2</DE></FR>  for  <IT>i</IT> = 1, 2, … , 9 (2)
The flow through each compartment, Qi, can be calculated as follows

Q<SUB><IT>i</IT></SUB> = <FR><NU>P<SUB>in,<IT>i</IT></SUB> − P<SUB>in,<IT>i</IT>+1</SUB></NU><DE><IT>R</IT><SUB><IT>i</IT></SUB> (P<SUB><IT>i</IT></SUB>)</DE></FR>  for  <IT>i</IT> = 1, 2, … , 9

Q<SUB>10</SUB> = <FR><NU>P<SUB>in,10</SUB></NU><DE><IT>R</IT><SUB>10</SUB></DE></FR> (3)
Since the flows through all compartments (Q1, Q2, ... , Q10) must be equal, the local pressures in the vessels can be calculated providing the parameters Ki and relations di(Pi) are known as well as the inlet and outlet pressures of the 1st and 10th compartment, respectively.

In this study, the outlet pressure of the 10th compartment is kept constant and equal to zero. Furthermore, it is assumed arbitrarily that the flow equaled 1 ml/s when the inlet pressure equals 90 mmHg at full dilation. All normalized passive pressure-diameter relations (see below) are taken to be identical for all compartments. For the resistance of the most distal compartment (R10), it is assumed that the inflow pressure Pin,10 equals 30 mmHg at a perfusion pressure of 90 mmHg at full vasodilation. Since the present model ignores the heterogeneity in topology of the coronary tree, it is not possible to estimate Ki from observed arteriolar length, density, and viscosity distributions. Therefore, a set of Ki values has been determined from a given pressure distribution in the coronary tree at full vasodilation, as explained below.

Three different pressure distributions are considered, given by the dashed lines in Fig. 1. These three curves in Fig. 1 result in three different resistance distributions at full dilation as given in Fig. 2. We will further refer to these distributions as: 1) distal dominant resistance, 2) gradual resistance distribution, and 3) proximal dominant resistance. As is clear from Fig. 1, our distal dominant and proximal dominant resistance distributions form an upper and lower bound for measured pressure distributions (1, 3), while the gradual resistance distribution forms an intermediate case. Without an evaluation of best fit, this gradual resistance distribution is considered the most realistic case of the three.


View larger version (25K):
[in this window]
[in a new window]
 
Fig. 1.   Definition of pressure distribution as function of vascular diameter. Dashed lines provide the pressure profiles for the fully dilated tree at a perfusion pressure of 90 mmHg. Solid, vertical lines provide the passive pressure-diameter relations for the 9 compartments. Local pressure (Pi) and diameter (di) in compartments with distal dominant (), gradual (black-triangle), and proximal dominant () resistance distribution. Data from Chilian et al. (3) measured in coronary arterioles of the beating cat epicardium, dilated with dipyridamole (open circle ) or in control conditions (). Interpolation of data measured by Chilian (1) in coronary subendocardium (down-triangle) and subepicardium (triangle ) arterioles of an isolated porcine heart.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 2.   Distribution of resistance over compartments in series. Inset at top: the 10 compartments model. The first 9 compartments represent the arterial tree with resistances depending on local pressure. The most distal compartment lumps the capillaries and the venules, and its resistance is kept constant. Bottom: the distal dominant (top), gradual (middle), and proximal dominant (bottom) resistance distribution for the vasodilated bed at a perfusion pressure of 90 mmHg. These resistance distributions correspond to the pressure profiles shown in Fig. 1.

Passive Pressure-Diameter Relations

Figure 3A shows data on passive pressure-diameter relations observed by various authors on isolated coronary vessels. The anatomical diameter of these vessels at 100 mmHg varied from 65 to 260 µm. As can be seen, the normalized passive pressure-diameter data from the different studies are rather similar over this range of vessel diameters. The function as provided by Liao and Kuo (19) is used to fit through these data points

<FR><NU><IT>d</IT><SUB><IT>i</IT></SUB></NU><DE><IT>d</IT><SUB>anat,<IT>i</IT></SUB></DE></FR> = <IT>d</IT><SUB>0</SUB> + (<IT>d</IT><SUB>max</SUB> − <IT>d</IT><SUB>0</SUB>) <FR><NU>P<SUB><IT>i</IT></SUB>/P<SUB>0.5</SUB></NU><DE>1 + P<SUB><IT>i</IT></SUB>/P<SUB>0.5</SUB></DE></FR>  <IT>i</IT> = 1, 2, … , 9

with <IT>d</IT><SUB>max</SUB> = 1.07, <IT>d</IT><SUB>0</SUB> = 0.7 and P<SUB>0.5</SUB> = 22.7 mmHg (4)
where danat,i equals the anatomical diameter defined as the passive diameter at local pressure of 100 mmHg. In contrast to Liao and Kuo (19), the passive pressure-normalized diameter relations are taken equal for all vessel orders. This curve is depicted by the solid line in Fig. 3A. To obtain the different pressure-diameter curves in each compartment, the anatomical diameters danat,i need to be chosen.



View larger version (3322K):
[in this window]
[in a new window]
 
Fig. 3.   Passive and active pressure diameter relations. A: passive pressure-diameter relations. Different symbols show diameters of various isolated coronary porcine vessels (6, 16-19, 23, 26) or human vessels (20), normalized to the diameter at 100 mmHg (danat), in some cases after interpolation of the original data. Solid line shows the passive pressure-diameter relation used for this study (Eq. 4). B: active pressure-diameter relations. Solid line corresponds with the passive pressure-diameter relation shown in A. Two dashed lines represent two extreme active pressure-diameter relations used for this study, the strong one with a decrease in diameter of 0.4%/mmHg and the weak one with 0.05%/mmHg. Below 5 mmHg, the vessels are thought to behave passively (Eq. 5). Data are from porcine coronary vessels, with basal tone, having a passive diameter of 91 µm at 100 mmHg (54 µm in vivo) () and from porcine coronary vessels with a passive diameter of 255 µm at 100 mmHg (180 µm in vivo) (black-diamond ). Data were normalized to the passive diameter at 100 mmHg.

Estimation of K and danat in the Compartments

Since it is assumed that the control of flow takes place in vessels with diameters ranging from 20 to 400 µm, we chose danat,1 = 400 µm and danat,9 = 20 µm. According to the application of Strahler ordering to topological findings, log (di) should be linear related to i (27). This determines the values of danat,i for i = 1, 2, ... , 8. Figure 1 provides the passive pressure-diameter relations for the vessels with different values of danat,i over the relevant pressure range, with Pi on the y-axis and di on the x-axis (vertical running solid curves). The values for di at Pi = 100 mmHg are indicated along the top axis of Fig. 1 and in Table 1. The intercepts of these pressure-diameter relations with the pressure distribution curves provide the actual diameter at full dilation of all compartments when the inlet pressure of the first compartment equals 90 mmHg. From these intercepts the values of Ki for a given pressure distribution can be calculated when the flow is defined. For the three resistance distributions, the three sets of Ki values are reported in Table 1.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Assigned diameters of the dilated vessels at 100 mmHg in the different compartments and the Ki values in each compartment

Active Pressure-Diameter Relations

Besides a standardized passive pressure-diameter relationship we also need a standardized active pressure-diameter relationship. This is implemented as shown in Fig. 3B. Up to a pressure of 5 mmHg the vessel is assumed to behave passively. For higher pressures the active pressure-diameter relationship is assumed to be a straight line with negative slope crossing the passive pressure-diameter relation at 5 mmHg. The magnitude of the slope represents the strength of the myogenic response. In Fig. 3B, the two extremes for myogenic strength used in this study are depicted. The active pressure-diameter relationship is calculated as follows

<IT>Eq. 4</IT> for P<SUB><IT>i</IT></SUB> ≤ 5 mmHg

<IT>d</IT><SUB><IT>i</IT></SUB>(P<SUB><IT>i</IT></SUB>) = <IT>d</IT><SUB><IT>i</IT></SUB>(5) − <FR><NU>M<SUB><IT>i</IT></SUB></NU><DE>100</DE></FR> <IT>d</IT><SUB>anat.<IT>i</IT></SUB>P<SUB><IT>i</IT></SUB> for P<SUB><IT>i</IT></SUB> > 5 mmHg 

<IT>i</IT> = 1, 2, … , 9 (5)
where Mi is the myogenic strength in percent per millimeter Hg.

Three different distributions of myogenic strength over the compartments are studied. Two distributions are homogeneous, but with different strengths; the weak and strong homogeneous distributions correspond to the two extremes in Fig. 3B. In the third distribution, myogenic strength is dependent on the anatomical diameter, with a maximum at 84 µm passive diameter. The three different distributions are depicted in Fig. 4.


View larger version (18K):
[in this window]
[in a new window]
 
Fig. 4.   Assumed distributions of myogenic strength, as defined in Fig. 3B, along the arterial tree; , strong response, homogeneously distributed; , weak response, homogeneously distributed; black-triangle, myogenic strength is varying with diameter between the weak and the strong case, having a maximum at an anatomical diameter of 84 µm.

Simulations

The diameters of the vessels, the local pressures in the vessels, and the resistances of the compartments are calculated for perfusion pressures ranging from 5 to 140 mmHg. For the calculations, an iterative procedure is used, as illustrated in Fig. 5, which is stopped when the flows (Qi values) in all compartments are equal with an accuracy of 10-10%.


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 5.   Flow chart illustrating the iterative simulation process. Diameters of vessels were calculated according to the passive, weak, strong, or diameter-dependent pressure-diameter relation. Depending on the resistance distribution, a set of Ki values was chosen (Table 1) to calculate the resistances. In this iterative procedure, Ri values and Qi values are continuously adjusted until flow in all compartments become equal. A maximum number of 51 iteration steps was needed before this stop criterion was fulfilled. Solutions were not dependent on the initially chosen resistances (Rstart).


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Autoregulation Curves

Figure 6 shows the pressure-flow relations for the three resistance distributions at full dilation and 90 mmHg. The passive pressure-flow relations are nearly identical; at equal pressure, flow differs only by 0.2% maximal between the different resistance distributions. The slight curvature of the passive pressure-flow line reflects the distensibility of the vessels in the first nine compartments. When the myogenic response is distributed homogeneously, the resistance distribution is hardly affecting the autoregulation curves. However, for the diameter-dependent myogenic response, the three different resistance distributions result in quite different autoregulation curves. When resistance is proximal dominant and myogenic tone diameter dependent, the autoregulation is weak due to the fact that distal vessels having the high myogenic strength hardly contribute to the total resistance.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 6.   Simulated coronary pressure-flow relations. For the dilated case, there is practically no difference between the three different resistance distributions. The stronger the myogenic strength, the better the autoregulation. The distribution of resistance only has an influence in cases where myogenic strength is diameter dependent.

Individual Diameter Responses in the Networks

Figure 7 shows the normalized diameter distribution as a function of the perfusion pressure when resistance is gradually distributed, and the different distributions in myogenic strength. Although quantitatively different, comparable results are found when resistance is dominant in the proximal compartments or in the distal compartments (results not shown). In the low perfusion pressure range, the diameters of the vessels in the different compartments increase with pressure when local pressure is smaller than 5 mmHg. At higher perfusion pressures the myogenic reduction in diameter emerges. The slope of the perfusion pressure-diameter relation at 90 mmHg is defined as the sensitivity of the myogenic response to changes in perfusion pressure. This sensitivity is provided for the nine compartments in Fig. 8, which relates to the three resistance distributions. In case of homogeneously distributed myogenic response (Fig. 7, left and middle; and Fig. 8), the strongest diameter reduction is always in the most proximal compartment. However, when myogenic reactivity is diameter dependent (Fig. 7, right; and Fig. 8), the strongest diameter reduction is in smaller vessels with the diameter of highest response dependent on the resistance distribution at full vasodilation. Maximal sensitivity of diameter to perfusion pressure is found at 123, 181, or 181 µm for the distal dominant, gradual distributed, or proximal dominant resistance distribution, respectively.


View larger version (26K):
[in this window]
[in a new window]
 
Fig. 7.   Dependence of individual active normalized diameters on perfusion pressure with gradual resistance distribution, reflecting the different distributions in myogenic strength. Different symbols at top reflect the compartments shown in the scheme at bottom. At low perfusion pressure, the diameters of the different compartments increase with pressure, because the local pressure under these conditions is smaller than 5 mmHg. At higher perfusion pressures, a myogenic reduction in diameter emerges. For the cases of homogeneous distribution of myogenic strength (left and middle), the maximum response is always in the proximal vessels (). With the diameter-dependent myogenic strength, the 4th compartment (181 µm, down-triangle) shows maximal myogenic diameter reduction. All curves were calculated for perfusion pressures of 0, 5, ... , 85, 89, 90, 91, 95, 100, ... , 140 mmHg.



View larger version (15K):
[in this window]
[in a new window]
 
Fig. 8.   Sensitivities of diameter to changes in perfusion pressure for the various compartments. These values are determined by calculating the slopes of the curves in Fig. 7 at a perfusion pressure of 90 mmHg; 90 mmHg is chosen because this is approximately the physiological pressure in vessels with a diameter of 400 µm. Homogeneously distributed myogenic reactivities ( and ) show the strongest diameter reduction in the most proximal compartment, irrespective of the resistance distribution. For the diameter-dependent myogenic tone, the strongest diameter reduction is dependent on the resistance distribution and emerges in vessels with diameters of 123, 181, or 181 µm for the distal dominant, gradual distributed, and proximal dominant resistance, respectively, whereas for the isolated vessels the maximal strength was at 84 µm.

Distributed Response of Resistance

Figure 9 provides the dependence of total resistance (top) and resistances of the compartments (bottom) on perfusion pressure for the case of gradual resistance distribution, reflecting the different distributions in myogenic strength. Despite quantitative differences, the two other resistance distributions resulted in comparable changes in resistance in the compartments as a function of perfusion pressure (results not shown). The sensitivity of resistance to perfusion pressure determined at 90 mmHg is demonstrated in Fig. 10 for all three resistance distributions.


View larger version (32K):
[in this window]
[in a new window]
 
Fig. 9.   Total resistance (Rtot) and individual resistances (Ri) of the compartments as a function of perfusion pressure (Pin,1) for the gradual resistance distribution, showing the different distributions in myogenic strength. The different symbols at top reflect the compartments shown in the scheme at bottom. Myogenic activity becomes clear by an increase in resistance with increasing pressure. Gradual resistance distribution is noticeable by a constant resistance for all compartments at low pressures where vessels are dilated. Myogenic response results in a large range of resistances at higher perfusion pressures. With the gradual resistance distribution and diameter-dependent myogenic strength, the resistance in compartments 3 and 4 becomes more dominant with increasing perfusion pressure. All curves were calculated for perfusion pressures of 0, 5, ... , 85, 89, 90, 91, 95, 100, ... , 140 mmHg.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 10.   Sensitivity of resistance to changes in perfusion pressure for the various compartments. These values are determined by calculating the slopes of the curves in Fig. 9 at a perfusion pressure of 90 mmHg. The highest sensitivities were in compartments 5, 4, and 2, representing vessels with diameters of 84, 123, and 268 µm, for the distal dominant, gradual distributed, and proximal dominant resistance, respectively.

The myogenic activity of a compartment should result in an increase of resistance with perfusion pressure. The importance of resistance distribution becomes apparent when considering the homogeneously distributed cases. When resistance is proximal dominant (Fig. 10, right), the resistance increase occurs mainly in the proximal compartments and annihilates the myogenic rise in resistance of the more distal compartments. When resistance is gradually distributed (Fig. 9, left and middle; and Fig. 10, middle), the myogenic resistance increase becomes more distributed as well, but the effect remains dominant in the proximal compartments. When resistance is distal dominant (Fig. 10, left), the strongest myogenic rise in resistance occurs at intermediate diameters notwithstanding the homogeneous distribution of myogenic strength.

For diameter-dependent myogenic strength (peak strength at 84 µm), maximal sensitivity of active resistance to perfusion pressure is found at 84, 123, or 268 µm for the distal dominant, gradual distributed, or proximal dominant resistance distribution, respectively (Fig. 9, right; Fig. 10).


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Our results demonstrate that myogenic activity of the large resistance vessels in the arterial tree attenuates the myogenic stimulus of the small vessels during changes in perfusion pressure. The more proximal the resistance, the stronger this effect. Hence, in flow control, the distribution of resistance is important for the role of myogenic behavior of individual vessels within the coronary network. This is most clearly demonstrated in the model with diameter-dependent myogenic tone. Peak myogenic tone was assumed here to occur in vessel segments with a passive diameter of 84 µm. However, in the network, peak myogenic diameter alterations to a change in perfusion pressure shifted to vessels with passive diameters in the range of 123 to 181 µm.

The model with homogeneously distributed resistance in combination with a diameter-dependent myogenic response seems to be the most realistic of the cases studied. It describes best the experimental data on pressure distribution in the dilated coronary bed (3) and diameter dependence of myogenic response of isolated vessels (19). However, coronary vascular trees show heterogeneous branching, and therefore resistance distributions along various paths may well be quite different, sometimes more closely resembling either the proximal dominant case or the distal dominant case. Hence, the different resistance distributions chosen provide insight into the importance of heterogeneity of the coronary tree in relation to physiological measurements. Obviously, the more local the measurement, e.g., diameter of individual in vivo vessels vs. pressure-flow lines in arteries, the more important heterogeneity will be.

Limitations of the Model

The model consists of nine compartments, which corresponds to the number of Strahler orders established for porcine coronary subtrees starting with a vessel of 400-500 µm in diameter (13, 27). A model consisting of a series of resistances can be rationalized as the equivalent of a symmetrical dichotomous branching tree. In that case, resistance of each compartment corresponds to the equivalent value of all vessels with a certain diameter and length in parallel. However, when applying the Strahler order vs. vessel diameter relation and the vessel diameter vs. length relation as reported by VanBavel and Spaan (27), the pressure distribution follows the top dashed curve of Fig. 1 (resistance dominant in the distal compartments). This curve deviates from experimental data on pressure-diameter relations in vivo, as is demonstrated in Fig. 1 as well. It has been shown that the artificial steep pressure drop at the smallest vessel diameter disappears when heterogeneity in branching characteristics is taken into account (27). To arrive at more realistic pressure distributions for our serial model, we changed the value of Ki for each compartment, which can be interpreted either as modifying the length (Li) or changing the number of parallel vessels (Ni). Also, changing the viscosity (eta i) in the compartment can modify Ki. Note that we did not alter the passive and active local pressure-diameter relations of the vessels representative for the compartments.

Two rheological effects are not accounted for in our model but may have an impact since viscosity may change with diameter change. As has been described by Pries et al. (24), apparent viscosity may decrease with increasing diameter for diameters lower than 40 µm, whereas for larger diameters, viscosity will increase. The latter effect is known as the Fahraeus-Lindqvist effect, whereas the former effect in apparent viscosity may relate to an endothelial surface layer (25, 29). However, for both rheological phenomena, it holds that their influence depends on the relative diameter variation of the vessels and the sensitivity of the apparent viscosity to these diameter changes. We calculated the variations in apparent viscosity based on diameter variation induced by a change in Pin,1 from 140 to 10 mmHg predicted with constant Ki. The maximal effect occurs in the case of the distal dominant resistance distribution and strong homogeneous myogenic strength. In this case, most of the diameter variation is limited to the first compartment, and apparent viscosity does change by 20%. For the gradual resistance distribution and diameter-dependent myogenic response, the maximal viscosity change is 8.5% and occurs in the 4th compartment. These changes are all due to the Fahraeus-Lindqvist effect. The other effect hardly plays a role, because vessels in the relevant compartments do not change much in diameter.

Especially at low pressures, there may be differences in the normalized passive pressure-diameter relationships between vessels of various caliber. However, for the sake of simplicity, these relations were taken independent of the size of the vessel. We also assumed myogenic characteristics of each compartment to be represented by a linear myogenic pressure-diameter relation above 5 mmHg. The current simplifications do not, however, influence the basic conclusions of our study, since myogenic strengths applied to the model cover the experimentally observed pressure-diameter relations found in isolated vessel studies, as is shown in Fig. 3B. In the case of diameter-dependent myogenic tone, the peak strength was taken for vessels with 84 µm in passive diameter at local pressure of 100 mmHg, which is in agreement with Liao and Kuo (19).

Our model does not contain an element describing metabolic vasodilation. Hence, the predictions by the model should be considered to be at a constant metabolic rate. Furthermore, possible other control mechanisms, including flow-dependent dilation, were purposely not taken into account to study the effect of myogenic response per se. Obviously, these modifications should be taken into account in a more comprehensive model. Griffith et al. (8) reported that the inhibition of endothelium-derived relaxing factor (EDRF) changes the resistance distribution in the microvascular network of a rabbit ear. Therefore, in vivo an interaction between flow-dependent dilation and the myogenic mechanism exists. Since the flow-dependent mechanism is generally assumed to be active in the larger vessels, the effect of myogenic constriction of these vessels on attenuation of local pressure changes in the smaller vessels may have been overestimated in the present study.

Previous Model Studies

Liao and Kuo (19) used a model approach similar to ours to investigate the potentiating effect of shear-dependent dilation by adenosine. That study used four compartments in series. Each compartment exhibited myogenic responsiveness as well. Granger (7) used a series-coupled model with three compartments to investigate the interaction of myogenic control, metabolic control, and flow-dependent dilation in the arterial tree. However, neither study gave special attention to the interaction between myogenic strength and resistance distribution over the compartments.

Model Predictions and Experimental Observations

Autoregulation curves. The pressure-flow lines show quite good autoregulation and support the notion that the myogenic response contributes to this phenomenon (5). To compare the slopes of our predicted autoregulation curves with in vivo observations, we may apply the correlation between flow (Q) on the one hand and pressure (Pin,1) and oxygen consumption (MVO2) on the other hand, as provided by Vergroesen et al. (28)
Q = <IT>a</IT>P<SUB>in,1</SUB> + <IT>b</IT>M<A><AC>V</AC><AC>˙</AC></A><SC>o</SC><SUB>2</SUB> + <IT>c</IT> (6)
where a equals the sensitivity of flow to Pin,1, and b indicates the sensitivity of flow to MVO2; c is the rest term. These authors calculated the relative sensitivity, a*, which equals the ratio of the slopes of the autoregulation curve (a) and of the pressure-flow curve at maximal dilation. They reported values for a* between 0.041 and 0.202. Our simulations are at constant MVO2, and therefore a* can be directly calculated from the slopes of the autoregulation curve and the fully dilated pressure-flow curve. These slopes, determined by linear regression between 60 and 120 mmHg, are shown in Table 2. Table 2 indicates that the pressure-flow lines with distal dominant and gradual resistance distribution and with the strong and diameter-dependent myogenic strength are within the range of the experimental found pressure-flow lines. When resistance is proximal and myogenic strength is diameter dependent, predicted autoregulation is too strong, whereas in the other combinations autoregulation is too weak. It should again be stressed that in in vivo studies, other control mechanisms are active as well, and therefore comparing quantitative results should be done with care.

                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Relative sensitivities of flow to perfusion pressure, calculated from the simulated autoregulation curves

Our model demonstrates flow reserve at pressures lower than actually found in experimental studies, where the dilated and the active pressure-flow lines start to diverge beyond 35 mmHg (21). However, coronary isolated vessel studies demonstrate myogenic tone at low pressures. This discrepancy could be due to either metabolic vasodilation at low perfusion pressure in the whole heart, vessel selection for isolated vessel studies, or alterations in vessel properties due to isolation or cannulation.

The current study demonstrates the difficulty of drawing conclusions with respect to myogenic tone and resistance distribution from coronary arterial pressure-flow lines. For example, for both the strong and the weak myogenic response, the pressure-flow line is independent of the resistance distribution (Fig. 6), notwithstanding quite large differences in the distribution of resistance responses over the compartments (Figs. 9 and 10).

When myogenic strength is assumed to be diameter dependent, autoregulation is poorer when resistance is proximal dominant. This, however, is obvious, since in this case most of the resistance is located in the proximal vessels with lower tone. When resistance is distal dominant or when resistance is gradually distributed, the model predictions differ only marginally. This is caused by the coincidence of strong myogenic tone with the location of the highest resistance.

Distribution of resistance and diameter variation induced by maximal vasodilation or changing perfusion pressure. In case myogenic tone is different in the various compartments, the maximal dilatory response to adenosine is expected to be different as well. Chilian and Layne (2) infused adenosine intracoronarily [Kanatsuka et al. (11) infused adenosine intravenously, and their data are not considered here] and demonstrated that smaller arterioles dilated more than proximal ones. The results are replotted in Fig. 11A: vessels ranging in diameter from 40 to 150 µm dilated on average 20%, and vessels between 150 and 300 µm dilated on average 10%, when perfusion pressure was kept at 96 mmHg. The percentage dilation predicted with the three resistance distributions and diameter-dependent myogenic strength, assuming that adenosine caused full dilation, is also shown in Fig. 11A. The dilation measured by Chilian and Layne (2) is smaller than we predict in our model, also for the most realistic case of gradual distributed resistance with diameter-dependent myogenic tone. This difference may be due to too high a level of vasoconstriction assumed in our model. However, this would then imply that the local spontaneous tone of the isolated vessels (6, 16-20, 26) is stronger than the tone in vivo under the working conditions of the experiments of Chilian and Layne (2). This may be possible since heart rate was high in the experiments of Chilian and Layne (2), ~188 beats/min, which should have resulted in lower coronary resistance at this baseline. This is also reflected in their coronary flow reserve at coronary pressure of 90 mmHg, which was a factor of 3 in the experiments and was a factor of 4 for our prediction applying the gradual resistance distribution.



View larger version (1719K):
[in this window]
[in a new window]
 
Fig. 11.   Diameter variation induced by maximal vasodilation and pressure steps. A: diameter variation induced by maximal vasodilation, with data of Chilian and Layne (2) (open circle ) obtained in anesthetized open-chest dogs by intracoronary infusion of adenosine while perfusion pressure was kept constant at ~96 mmHg. Solid symbols show the simulated results for the different combinations in resistance distribution and myogenic strength. Relative diameter change from myogenic constriction to maximal dilation, [di,pas(90) - di,act(90)]/di,act(90), was calculated at a constant perfusion pressure of 90 mmHg; di,pas(90) values are the passive diameters for each compartment according to Eq. 4, with Pin,1 = 90 mmHg; di,act(90) values are the active diameters for each compartment according to Eq. 5, with Pin,1 = 90 mmHg. In B, open symbols are relative diameter changes measured by Chilian and Layne (2) (open circle ) and Kanatsuka et al. (12) (triangle ) in anesthetized open-chest dogs after reductions in perfusion pressure from approximately 90 mmHg to 40 mmHg, [di,act(40) - di,act(90)]/di,act(90). Solid symbols show the results of simulating this pressure step for the different combinations in resistance distribution and myogenic strength.

Chilian and Layne (2) and Kanatsuka et al. (12) have induced pressure steps by an occluder around a major coronary artery. Pressure reductions at steady state varied from ~90 mmHg to pressures in the range of 90 to 30 mmHg. The major outcome was that smaller vessels dilated more than larger ones, except for a pressure step down to 30 mmHg, in which case the small vessels (<150 µm) did not dilate at all. We have compiled relative diameter changes of the different vessels as induced by pressure reduction from 90 to 40 mmHg from the two studies in Fig. 11B (2, 12). In Fig. 11, we also plotted our model predictions for five cases. Only the cases with diameter-dependent myogenic tone predict the dominant dilation in the smaller vessels (down to 50 µm). In contrast, the simulated homogeneously distributed weak and strong myogenic responses predict a continuous decreasing vessel diameter. The simulations with diameter-dependent myogenic tone further predict that the dilation becomes less in vessels with very small diameters, but experimental data in this range are absent. For the large vessels Kanatsuka et al. (12) found on average constriction of about 10%, in contrast to Chilian et al. (2), who found on average no variation in diameter. This constriction can be explained by a weaker myogenic response than assumed here or by involvement of flow-dependent dilation.

In conclusion, this study shows that interaction between different vessels in the coronary arterial tree influences the effectiveness of myogenic flow control. Myogenic activity of larger vessels attenuates the myogenic activity of smaller vessels. The location for myogenic control depends not only on the myogenic strength of the vessels in isolation but also on the resistance distribution. For a given flow demand, the distributed myogenic response is consistent with both experimentally found diameter changes and the course of the pressure-flow line. Hence, the myogenic response is sufficient to explain autoregulation, notwithstanding the fact that other control mechanisms have their own role in flow control. Therefore, the myogenic response obscures the in vivo study of other mechanisms mediating autoregulation by its potential for flow control. This study emphasizes the importance of knowledge of the distribution of resistance and control mechanisms when studying vessels in the coronary tree.


    ACKNOWLEDGEMENTS

This work is supported by Netherlands Heart Foundation Grant 95020 (to J. Dankelman).


    FOOTNOTES

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. §1734 solely to indicate this fact.

Address for reprint requests and other correspondence: J. A. E. Spaan, Dept. of Medical Physics, Cardiovascular Research Institute Amsterdam, Academic Medical Center, Univ. of Amsterdam, PO Box 22700, 1100 DE Amsterdam, The Netherlands (E-mail: j.a.spaan{at}amc.uva.nl).

Received 2 August 1999; accepted in final form 17 November 1999.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Chilian, WM. Microvascular pressures and resistances in the left ventricular subepicardium and subendocardium. Circ Res 69: 561-570, 1991[Abstract/Free Full Text].

2.   Chilian, WM, and Layne SM. Coronary microvascular responses to reductions in perfusion pressure. Evidence for persistent arteriolar vasomotor tone during coronary hypoperfusion. Circ Res 66: 1227-1238, 1990[Abstract/Free Full Text].

3.   Chilian, WM, Layne SM, Klausner EC, Eastham CL, and Marcus ML. Redistribution of coronary microvascular resistance produced by dipyridamole. Am J Physiol Heart Circ Physiol 256: H383-H390, 1989[Abstract/Free Full Text].

4.   Davis, MJ. Myogenic response gradient in an arteriolar network. Am J Physiol Heart Circ Physiol 264: H2168-H2179, 1993[Abstract/Free Full Text].

5.   Dole, WP. Autoregulation of the coronary circulation. Prog Cardiovasc Dis 29: 293-323, 1987[ISI][Medline].

6.   Giezeman, MJMM, VanBavel E, Grimbergen CA, and Spaan JAE Compliance of isolated porcine coronary small arteries and coronary pressure-flow relations. Am J Physiol Heart Circ Physiol 267: H1190-H1198, 1994[Abstract/Free Full Text].

7.   Granger, HJ. Coordinated microvascular reactions mediated by series coupling of metabolic, myogenic and flow-sensitive resistance elements (Abstract). FASEB J 3: A1387, 1989.

8.   Griffith, TM, Edwards DH, Davies RL, Harrison TJ, and Evans KT. EDRF coordinates the behaviour of vascular resistance vessels. Nature 329: 442-445, 1987[Medline].

9.   Johnson, PC. The myogenic response. In: Handbook of Physiology. The Cardiovascular System. Vascular Smooth Muscle. Bethesda, MD: Am. Physiol. Soc, 1980, sect. 2, vol. II, chapt. 15, p. 409-442.

10.   Johnson, PC. The myogenic response in the microcirculation and its interaction with other control systems. J Hypertens Suppl 7: S33-S39, 1989[Medline].

11.   Kanatsuka, H, Lamping KG, Eastham CL, Dellsperger KC, and Marcus ML. Comparison of the effects of increased myocardial oxygen consumption and adenosine on the coronary microvascular resistance. Circ Res 65: 1296-1305, 1989[Abstract/Free Full Text].

12.   Kanatsuka, H, Lamping KG, Eastham CL, and Marcus ML. Heterogeneous changes in epimyocardial microvascular size during graded coronary stenosis. Evidence of the microvascular site for autoregulation. Circ Res 66: 389-396, 1990[Abstract/Free Full Text].

13.   Kassab, GS, Berkley J, and Fung YC. Analysis of pig's coronary arterial blood flow with detailed anatomical data. Ann Biomed Eng 25: 204-217, 1997[ISI][Medline].

14.   Kassab, GS, Pallencaoe E, Schatz A, and Fung YC. Longitudinal position matrix of the pig coronary vasculature and its hemodynamic implications. Am J Physiol Heart Circ Physiol 273: H2832-H2842, 1997[Abstract/Free Full Text].

15.   Kassab, GS, Rider CA, Tang NJ, and Fung YC. Morphometry of pig coronary arterial trees. Am J Physiol Heart Circ Physiol 265: H350-H365, 1993[Abstract/Free Full Text].

16.   Kuo, L, Chilian WM, and Davis MJ. Coronary arteriolar myogenic response is independent of endothelium. Circ Res 66: 860-866, 1990[Abstract/Free Full Text].

17.   Kuo, L, Chilian WM, and Davis MJ. Interaction of pressure- and flow-induced responses in porcine coronary resistance vessels. Am J Physiol Heart Circ Physiol 261: H1706-H1715, 1991[Abstract/Free Full Text].

18.   Kuo, L, Davis MJ, and Chilian WM. Myogenic activity in isolated subepicardial and subendocardial coronary arterioles. Am J Physiol Heart Circ Physiol 255: H1558-H1562, 1988[Abstract/Free Full Text].

19.   Liao, JC, and Kuo L. Interaction between adenosine and flow-induced dilation in coronary microvascular network. Am J Physiol Heart Circ Physiol 272: H1571-H1581, 1997[Abstract/Free Full Text].

20.   Miller, FJJ, Dellsperger KC, and Gutterman DD. Myogenic constriction of human coronary arterioles. Am J Physiol Heart Circ Physiol 273: H257-H264, 1997[Abstract/Free Full Text].

21.   Mosher, P, Ross J, Jr, McFate PA, and Shaw RF. Control of coronary blood flow by an autoregulatory mechanism. Circ Res 14: 250-259, 1964[Abstract/Free Full Text].

22.   Muller, JM, Davis MJ, and Chilian WM. Integrated regulation of pressure and flow in the coronary microcirculation. Cardiovasc Res 32: 668-678, 1996[ISI][Medline].

23.   Nakayama, K, Osol G, and Halpern W. Differential contractile responses of pressurized porcine coronary resistance-sized and conductance coronary arteries to acetylcholine, histamine and prostaglandin F2 alpha. Blood Vessels 26: 235-245, 1989[ISI][Medline].

24.   Pries, AR, Secomb TW, Gessner T, Sperandio MB, Gross JF, and Gaehtgens P. Resistance to blood flow in microvessels in vivo. Circ Res 75: 904-915, 1994[Abstract/Free Full Text].

25.   Pries, AR, Secomb TW, Jacobs H, Sperandio M, Osterloh K, and Gaehtgens P. Microvascular blood flow resistance: role of endothelial surface layer. Am J Physiol Heart Circ Physiol 273: H2272-H2279, 1997.

26.   Rajagopalan, S, Dube S, and Canty JM, Jr. Regulation of coronary diameter by myogenic mechanisms in arterial microvessels greater than 100 µm in diameter. Am J Physiol Heart Circ Physiol 268: H788-H793, 1995[Abstract/Free Full Text].

27.   VanBavel, E, and Spaan JAE Branching patterns in the porcine coronary arterial tree. Estimation of flow heterogeneity. Circ Res 71: 1200-1212, 1992[Abstract/Free Full Text].

28.   Vergroesen, I, Noble MI, Wieringa PA, and Spaan JA. Quantification of O2 consumption and arterial pressure as independent determinants of coronary flow. Am J Physiol Heart Circ Physiol 252: H545-H553, 1987[Abstract/Free Full Text].

29.  Vink H, Constantinescu A, and Spaan JAE. Oxidized lipoproteins degrade the endothelial surface layer: implications for platelet-endothelial cell adhesion. Circulation In press.


Am J Physiol Heart Circ Physiol 278(5):H1490-H1499
0363-6135/00 $5.00 Copyright © 2000 the American Physiological Society



This article has been cited by other articles:


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
B. E. Carlson, J. C. Arciero, and T. W. Secomb
Theoretical model of blood flow autoregulation: roles of myogenic, shear-dependent, and metabolic responses
Am J Physiol Heart Circ Physiol, October 1, 2008; 295(4): H1572 - H1579.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Renal Physiol.Home page
N. Kleinstreuer, T. David, M. J. Plank, and Z. Endre
Dynamic myogenic autoregulation in the rat kidney: a whole-organ model
Am J Physiol Renal Physiol, June 1, 2008; 294(6): F1453 - F1464.
[Abstract] [Full Text] [PDF]


Home page
Physiol. Rev.Home page
N. Westerhof, C. Boer, R. R. Lamberts, and P. Sipkema
Cross-talk between cardiac muscle and coronary vasculature.
Physiol Rev, October 1, 2006; 86(4): 1263 - 1308.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
F. Spohr, A. J. M. Cornelissen, C. Busch, M. M. Gebhard, J. Motsch, E. O. Martin, and J. Weimann
Role of endogenous nitric oxide in endotoxin-induced alteration of hypoxic pulmonary vasoconstriction in mice
Am J Physiol Heart Circ Physiol, August 1, 2005; 289(2): H823 - H831.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
A. J. M. Cornelissen, J. Dankelman, E. VanBavel, and J. A. E. Spaan
Balance between myogenic, flow-dependent, and metabolic flow control in coronary arterial tree: a model study
Am J Physiol Heart Circ Physiol, June 1, 2002; 282(6): H2224 - H2237.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
G. S. Kassab and S. Molloi
Cross-sectional area and volume compliance of porcine left coronary arteries
Am J Physiol Heart Circ Physiol, August 1, 2001; 281(2): H623 - H628.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
D. Merkus, I. Vergroesen, O. Hiramatsu, H. Tachibana, H. Nakamoto, E. Toyota, M. Goto, Y. Ogasawara, J. A. E. Spaan, and F. Kajiya
Stenosis differentially affects subendocardial and subepicardial arterioles in vivo
Am J Physiol Heart Circ Physiol, April 1, 2001; 280(4): H1674 - H1682.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
A. J. M. Cornelissen, J. Dankelman, E. VanBavel, and J. A. E. Spaan
Balance between myogenic, flow-dependent, and metabolic flow control in coronary arterial tree: a model study
Am J Physiol Heart Circ Physiol, June 1, 2002; 282(6): H2224 - H2237.
[Abstract] [Full Text] [PDF]


<