Experimental Physiology
	

Email Content Delivery
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Experimental Physiology 90.5 pp 715-726
DOI: 10.1113/expphysiol.2005.030189
© The Physiological Society 2005
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
90/5/715    most recent
expphysiol.2005.030189v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
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 Google Scholar
Google Scholar
Right arrow Articles by Ferreira, L. F
Right arrow Articles by Barstow, T. J
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ferreira, L. F
Right arrow Articles by Barstow, T. J
Related Collections
Right arrow Human, Environmental & Exercise

Kinetics of estimated human muscle capillary blood flow during recovery from exercise

Leonardo F Ferreira1, Allison J Harper2, Dana K Townsend1, Barbara J Lutjemeier1 and Thomas J Barstow2

Departments of 1 Anatomy and Physiology2 Kinesiology, Kansas State University, Manhattan, KA 66506-0302, USA


    Abstract
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
The kinetic characteristics of muscle capillary blood flow {eph_162_mu1} during recovery from exercise are controversial (e.g. one versus two phases). Furthermore, it is not clear how the overall {eph_162_mu2} kinetics are temporally associated with muscle oxygen uptake {eph_162_mu3} kinetics. To address these issues, we examined the kinetics of {eph_162_mu4} estimated from the rearrangement of the Fick equation {eph_162_mu5} using the kinetics of pulmonary {eph_162_mu6} ({eph_162_mu7}, primary component) and deoxy-haemoglobin concentration ([HHb]) as indices of {eph_162_mu8} and C(a – v)O2 (arterio-venous oxygen difference) kinetics, respectively. {eph_162_mu9} (l min–1) was measured breath by breath and [HHb] (µm) was measured by near infrared spectroscopy during moderate (M; below lactate threshold, LT) and heavy exercise (H, above LT) in nine subjects. The kinetics of {eph_162_mu10} were biphasic, with an initial fast phase ({tau}I; M = 9.3 ± 4.9 s and H = 6.0 ± 3.8 s) followed by a slower phase 2 ({tau}P; M = 29.9 ± 8.6 s and H = 47.7 ± 26.0 s). For moderate exercise, the overall kinetics of {eph_162_mu11} (mean response time [MRT], 36.1 ± 8.6 s) were significantly slower than the kinetics of {eph_162_mu12} ({tau}P; 27.8 ± 5.3 s) and [HHb] (MRT for [HHb]; 16.2 ± 6.3 s). However, for heavy exercise, there was no significant difference between MRT-[HHb] (34.7 ± 10.4 s) and {tau}P for {eph_162_mu13} (32.3 ± 6.7 s), while MRT for {eph_162_mu14} (48.7 ± 21.8 s) was significantly slower than MRT for [HHb] and {tau}P for {eph_162_mu15}. In conclusion, during recovery from exercise the estimated {eph_162_mu16} kinetics were biphasic, showing an early rapid decrease in blood flow. In addition, the overall kinetics of {eph_162_mu17} were slower than the estimated {eph_162_mu18} kinetics.

(Received 31 January 2005; accepted after revision 18 May 2005; first published online 20 May 2005)
Corresponding author T. J. Barstow: Department of Kinesiology, 1A Natatorium, Kansas State University, Manhattan, KA 66506-0302, USA. Email: tbarsto{at}ksu.edu


    Introduction
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Muscle blood flow {eph_162_mu19}, coupled with arterial O2 content, will dictate the dynamics of O2 delivery during the recovery period following cessation of exercise. Hence, the kinetics of {eph_162_mu20} can have profound effects on the recuperation of muscle metabolism to pre-exercise levels depending on the population studied (trained subjects, Haseler et al. 1999; heart failure patients, Hanada et al. 2000). However, the temporal (kinetic) characteristics of {eph_162_mu21} during recovery from exercise remain unclear.

It has been reported that the recovery of {eph_162_mu22} was not exponential in human larger vessels (moderate and heavy exercise, Van Leeuwen et al. 1992) and in rat muscle feed arteries (high-intensity muscle contractions, Lash, 1994). In contrast, other investigators have used exponential equations to describe the {eph_162_mu23} response during recovery from exercise (Shoemaker et al. 1994, 1997; Whipp et al. 1995). However, these studies either did not attempt to distinguish (Shoemaker et al. 1994; Whipp et al. 1995), or did not report (Shoemaker et al. 1997) whether the best description of the data was obtained by mono- or two-exponential functions. On the other hand, Leyk et al. (1999) contended that the recovery of blood flow did not have a fast component such as observed following the onset of exercise. Conversely, Van Leeuwen et al. (1992) reported an initial rapid decrease in {eph_162_mu24}, followed by a slower response after cessation of exercise, suggesting that the {eph_162_mu25} response had at least two distinct phases. Thus, the question of whether {eph_162_mu26} during recovery has one or two phases remains unresolved. This is an important issue because distinct mechanisms regulating muscle blood flow may operate in different phases of the response. For instance, following the onset of exercise {eph_162_mu27} is biphasic and the initial phase (first 15–20 s) appears to be determined by a complex interaction between the muscle pump and vasodilation (Tschakovsky & Sheriff, 2004) while the subsequent increase in {eph_162_mu28} is determined by neural, vascular and metabolic factors (Laughlin & Korzick, 2001). Furthermore, computer simulations of {eph_162_mu29} and muscle oxygen uptake {eph_162_mu30} responses following the onset of exercise suggest that, since {eph_162_mu31} on-kinetics are biphasic, blood flow in the initial 10–30 s of exercise is critical in determining the adequacy of O2 delivery to the contracting muscle (Ferreira et al. 2005a). Since a mechanical factor, such as the muscle pump, is involved in the early phase of {eph_162_mu32} response following the onset of exercise (e.g. Sheriff & Hakeman, 2001), a qualitative symmetry of onset and recovery kinetics (i.e. biphasic recovery of {eph_162_mu33}) might be expected.

The muscle microcirculation is the primary site of O2 exchange. Thus, muscle capillary blood flow {eph_162_mu34} and {eph_162_mu35} kinetics will determine the dynamics of capillary O2 extraction {eph_162_mu36}. In humans, the dynamics of O2 extraction in the microcirculation have been estimated non-invasively by the deoxy-haemoglobin concentration ([HHb]) signal from near infrared spectroscopy (NIRS; DeLorey et al. 2003; Grassi et al. 2003). Although the kinetics of [HHb] following the onset of exercise have been previously described (DeLorey et al. 2003; Grassi et al. 2003; Ferreira et al. 2005b), the time course of [HHb] during recovery from exercise is presently unknown. The study of microvascular O2 pressure (Pmv,O2) recovery kinetics (approximately equivalent to the time course of O2 extraction) in the rat muscle microcirculation has revealed extremely different responses following the onset and cessation of muscle contraction (McDonough et al. 2001), indicating asymmetry of the temporal association between {eph_162_mu37} and {eph_162_mu38} for the onset and recovery phases. In humans, following the transition from unloaded exercise to exercise the kinetics of {eph_162_mu39} were similar to the {eph_162_mu40} kinetics (Grassi et al. 1996; Koga et al. 2005; Ferreira et al. 2005b). However, it is unclear whether blood flow and O2 utilization remain coupled after cessation of exercise. Shoemaker et al. (1994) found similar {eph_162_mu41} (femoral artery) and estimated {eph_162_mu42} kinetics (but see Discussion), while data from other studies indicate that {eph_162_mu43} recovered at a slower rate than {eph_162_mu44} after moderate (Leyk et al. 1999; Van Beekvelt et al. 2001) and heavy exercise (Bangsbo et al. 1994; Radegran & Saltin, 1999), although the kinetic parameters were not determined. Importantly, if the recovery kinetics of {eph_162_mu45} were faster than the {eph_162_mu46} kinetics, such as following rest-to-exercise transitions (MacDonald et al. 1998; Behnke et al. 2002), there would be an O2 delivery limitation to the restoration of muscle metabolism (Barstow et al. 1990). Therefore, it could be expected that the recovery of {eph_162_mu47} would be slower than that of {eph_162_mu48} in normal healthy individuals, suggesting that O2 delivery would not be a limiting factor to the recovery of muscle metabolism, as recently shown in sedentary individuals (Haseler et al. 2004).

Mechanistically, it is important to distinguish whether the kinetics of {eph_162_mu49} are a first-order (i.e. one phase) or higher-order response. If the {eph_162_mu50} kinetics are first order it implies, but does not prove, that a single mechanism is controlling {eph_162_mu51} during recovery from exercise. Likewise, higher-order responses suggest that multiple mechanisms operate to maintain {eph_162_mu52} elevated during recovery, and these vary as a function of time. Insights gained from this analysis would also be useful to understand the dynamics of muscle oxygenation {eph_162_mu53} reported in physiological (Chance et al. 1992; McCully et al. 1994) and pathological conditions (Belardinelli et al. 1997; Hanada et al. 2000).

The purpose of the present study was to examine the estimated kinetics of {eph_162_mu54} to determine whether during recovery from exercise blood flow in the microcirculation has one or two phases. We hypothesized that the putative muscle pump effect observed following the onset of exercise would also be evident after cessation of exercise, leading to a biphasic recovery of {eph_162_mu55} similar to, but a mirror image of, the onset kinetics (Ferreira et al. 2005b). We also tested the hypothesis that the overall recovery kinetics of {eph_162_mu56} would be slower than the {eph_162_mu57} kinetics, which would indicate that different mechanisms are involved in the regulation of blood flow kinetics during recovery from cycling exercise compared to the onset, when the overall kinetics of {eph_162_mu58} were similar to those of {eph_162_mu59} (transition from unloaded exercise to exercise, Grassi et al. 1996; Ferreira et al. 2005b; Koga et al. 2005).


    Methods
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Nine healthy subjects (7 men and 2 women) with mean (± S.D.) age 24.7 ± 6.3 years, body weight 67.9 ± 12.2 kg and height 175.4 ± 13.1 cm participated in this study. All procedures were explained to each subject, who provided signed written informed consent before enrolment in the study. The protocol was approved by the Institutional Review Board for Research Involving Human Subjects at Kansas State University and followed the principles outlined in the Declaration of Helsinki.

The exercise protocol was performed on three separate days. All exercise tests were performed on an electronically braked cycle ergometer (Corival 400, Lode, The Netherlands). On the first test day, the subjects performed an incremental (ramp) exercise test for determination of peak oxygen uptake {eph_162_mu60} and the estimated lactate threshold (LT). From these, the work rates calculated to elicit 90% of the LT {eph_162_mu61} (90% LT) and a {eph_162_mu62} half-way between the LT and {eph_162_mu63} (50% {eph_162_mu64}) were determined. On each of the two consecutive visits to the laboratory, subjects performed four bouts of constant work rate exercise, the first three bouts at 90% LT (6 min each) interspersed by 6 min at 20 W and the last bout at 50% {Delta} (8 min duration). The first bout was preceded by 4 min of baseline pedaling at 20 W. Therefore, a total of six transitions to the 90% LT and two transitions to the 50% {Delta} work rates were performed. Three subjects performed four transitions at the 90% LT work rate due to schedule conflicts. Pulmonary {eph_162_mu65} and muscle oxygenation by NIRS were measured in all subjects and transitions.

Pulmonary gas exchange ({eph_162_mu66} and {eph_162_mu67}) and minute expired ventilation {eph_162_mu68} were measured breath by-breath (CardiO2, MedGraphics, St Paul, MN, USA). Before each exercise test the volume signal was calibrated using a 3 l syringe, while the O2 and CO2 analysers were calibrated with gases of known concentration. Heart rate was recorded from the electrocardiogram using a modified lead I configuration and stored in the breath-by-breath file.

Muscle oxygenation was evaluated by a frequency-domain multidistance (FDMD) near-infrared spectroscopy system (OxiplexTS, ISS, Champaign, IL, USA) operating at two wavelengths (690 and 830 nm). The light-emitting diodes and light detector were connected to a plastic probe by optical fibres. The probe consisted of light source-detector separations of 2.0, 2.5, 3.0 and 3.5 cm for both wavelengths. The frequency modulation of laser intensity was 110 MHz and the heterodyne detection was performed at a 5 kHz cross-correlation frequency. The output frequency for data storage was selected as 31.25 Hz. After the site had been carefully shaved and dried, the probe was positioned longitudinally over the belly of the vastus lateralis muscle ~15 cm above the patella, bound to the skin with a skin cement (Skin-Bond, Smith & Nephew, Largo, FL, USA) and secured with Velcro straps around the thigh. The probe position was marked to check for any sliding and for accurate repositioning on subsequent test days. No movement (sliding) was observed in any exercise test. The near-infrared spectrometer was calibrated on each test day according to the manufacturer's recommendations. To estimate the actual values of [HHb], [HbO2] and total haemoglobin concentration ([THb] = [HHb] + [HbO2]) dynamic measurements of the reduced scattering coefficient were performed continuously throughout the exercise and recovery periods and automatically incorporated into the algorithms utilized by the manufacturer's software.

Kinetics analysis

The breath-by-breath O2 and 31.25 Hz NIRS-oxygenation data were converted to second-by-second values, time aligned to the end of exercise and ensemble averaged for each subject, to generate a single data set for each variable, at each exercise intensity. The kinetics of {eph_162_mu69}, [HHb] and {eph_162_mu70} (see below) were determined by non-linear regression using a least squares technique (Marquadt-Levenberg, SigmaPlot 2001, Systat Software, Point Richmond, CA). The model used for fitting the [HHb] response was:


{eph_162_m1}

(1)
(see definitions below). For responses without ‘undershoots’ or slow decreases in [HHb] after the fast component (e.g. Fig. 1) the parameter AS converged to zero, which is consistent with a monoexponential model. Whenever an undershoot was seen the slow component refers to the increase in [HHb] above the nadir achieved during the recovery period.



View larger version (21K):
[in this window]
[in a new window]
 
Figure 1.  Dynamics of deoxyhaemoglobin concentration ([HHb]) from three representative subjects
Note that the responses are qualitatively similar during the initial 60 s, showing a short time delay and ‘rapid’ decrease (steep slope) thereafter. In general, these responses qualitatively resemble the dynamics of C(a – v)O2 shown by Van Beekvelt et al. (2001). The source of variability between subjects is unclear, but this is also seen for total haemoglobin concentration (Fig. 2) and probably involves differences in adipose tissue thickness, muscle capillarization ({approx} capillary blood volume) and tissue scattering properties. Similarly, changes in [HHb] from moderate to heavy exercise within subjects can be misleading since {Delta}[HHb] from unloaded cycling to peak work rate can range from 2.4 to 19.5 µM (Barstow et al. 2004).

 
The {eph_162_mu71} response during recovery from exercise was fitted with the following model:


{eph_162_m2}

(2)
where, in eqns (1) and (2), EE refers to end of exercise; A, amplitude; TD, time delay; and {tau}, the time constant of each phase of the response (subscripts I, P and S). For recovery from exercise below the LT only the initial and primary component were included in the non-linear regression function, whereas for recovery from exercise above the LT two different models were used. In model I, TDS was substituted for TDP (i.e. primary and slow components of recovery began together at the end of phase 1), whereas for model II, TDS != TDP (i.e. the slow component emerged after an independent time delay). The two models were then compared by a F test to determine whether model II described the data significantly better then model I (Motulsky & Ransnas, 1987). If so, model II was used in the procedures to estimate {eph_162_mu72} kinetics during recovery from heavy exercise, otherwise model I was employed. The time constant of the primary (fast) component, and the physiological interpretation, are different when the data are fitted by models I or II (Barstow & Mole, 1991). However, we utilized this approach for analysing the {eph_162_mu73} instead of defining a priori the model for non-linear regression because we were interested in using the {eph_162_mu74} kinetics to estimate the temporal profile of {eph_162_mu75} during recovery (see below). In this circumstance, the model that better described the {eph_162_mu76} data appeared to be the appropriate choice.

The initial (cardiodynamic) component of {eph_162_mu77} was described up to TDP and the amplitude of the response at TDP (A'I) was calculated as A'I = AI (1 – e–(TDP/{tau}I)). The relevant amplitude of the primary component was calculated as A'P = A'I + AP.

Muscle capillary blood flow

The {eph_162_mu78} response to exercise was derived from the kinetics of {eph_162_mu79} (time constants and amplitude of response) and raw [HHb] data. The methods employed to derive {eph_162_mu80} have been described in detail elsewhere (Ferreira et al. 2005b). Briefly, the kinetics of {eph_162_mu81} (primary component) during recovery from constant work rate exercise closely represents the muscle phosphocreatine and, presumably, {eph_162_mu82} kinetics (Rossiter et al. 2002), whereas [HHb] measured by NIRS is thought to reflect the muscle capillary fractional O2 extraction (DeLorey et al. 2003; Grassi et al. 2003). Therefore, by rearranging the Fick equation, the temporal characteristics of {eph_162_mu83} can be estimated from the ratio of {eph_162_mu84} to [HHb], specifically:


{eph_162_m3}

(3)
Under these circumstances, absolute values of {eph_162_mu85} are uncertain but its temporal (kinetic) characteristics are preserved (Ferreira et al. 2005b). The muscle O2 response was estimated using the kinetic parameters obtained from the curve fitting of {eph_162_mu86}, i.e. by assuming that the recovery of {eph_162_mu87} was exponential, starting at the end of exercise with the time constant determined for {eph_162_mu88} (primary component) and amplitude equal to A'P (Ferreira et al. 2005b). The resulting {eph_162_mu89} was used in the estimation of {eph_162_mu90} kinetics (e.g. Figs 3 and 4). The time course of {eph_162_mu91} was analysed by exponential equations as described above in eqn (2) (model II), where {eph_162_mu92} is substituted for {eph_162_mu93}. This method is consistent with previous studies (Shoemaker et al. 1994, 1997). It is not known, at present, whether the {eph_162_mu94} profile during recovery has an initial fast component as seen following the onset of exercise (Kindig et al. 2002; see Discussion). Therefore, a model not including the initial (fast) phase 1 (monoexponential, moderate exercise and two-exponential, heavy exercise) was also employed for non-linear regression analysis and an F test was used to determine the better description of the data. The mean response time (MRT) for {eph_162_mu95}, which approximates the time to reach 63% of the response, was calculated as:


{eph_162_m4}

(4)
where the parameters used in the equation are defined above.



View larger version (16K):
[in this window]
[in a new window]
 
Figure 3.  Representative data for recovery from moderate exercise (90% LT)
A, muscle {eph_162_mu208} based on kinetic parameters of pulmonary {eph_162_mu209} (primary component); B, deoxyhaemoglobin concentration ([HHb], µM); and C, estimated temporal profile of muscle capillary blood flow ({circ}, {eph_162_mu210}) obtained as {eph_162_mu211} and best regression fit (continuous line) obtained from non-linear regression using eqn (2). Time constant of initial ({tau}I) and primary components ({tau}P) as well as mean response time (MRT) are shown. Dotted line demarcates the Phase 1–2 transition of the estimated {eph_162_mu212} response. Note that the [HHb] response is consistent with an O2 extraction profile resulting from the on-going monoexponential decrease in muscle {eph_162_mu213}, and an initial fast decrease in muscle blood flow that projects to a plateau at an elevated value compared to end of recovery, characteristic of a biphasic response.

 


View larger version (17K):
[in this window]
[in a new window]
 
Figure 4.  Representative data of recovery from heavy exercise (50% {Delta})
See Fig. 3 for abbreviations and further details on panels A, B and C. The value of {tau} shown in A was obtained from fitting the {eph_162_mu214} response corresponding to the primary component of {eph_162_mu215} during recovery. See Methods for further details.

 
Data analysis for comparison of overall kinetics during recovery from heavy exercise

To compare the overall kinetics of the estimated {eph_162_mu96} and {eph_162_mu97}, we reanalysed the pulmonary {eph_162_mu98} response of each subject with a two-exponential model (phases 1–2 of eqn (2)) over the same region of data for the initial and primary components of {eph_162_mu99}, i.e. the {eph_162_mu100} data corresponding in time to the slow component for {eph_162_mu101} (from t {approx} TDS of {eph_162_mu102} to 480 s) were ignored in the non-linear regression. The resulting recovery kinetics of estimated {eph_162_mu103} ({tau}Pb) were then compared to the MRT for both [HHb] and {eph_162_mu104}.

Statistical analysis

To determine significant differences between two means, Student's two-tailed paired t test, or Wilcoxon's test when appropriate, was performed. A repeated-measures analysis of variance was conducted to compare more than two means and the Tukey–Kramer's post hoc test was used for pairwise comparisons. The relationship between two variables was analysed by the Pearson product–moment correlation. Significance was accepted when P ≤ 0.05. All tests were conducted using a commercial statistical software package (NCSS 2000, NCSS Statistical Software, Kaysville, UT). Values were reported as means ± S.D.


    Results
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
{eph_162_mu105} was 48.8 ± 7.0 ml kg–1 min–1 and the estimated LT occurred at 56.3 ± 8.5% {eph_162_mu106}. The work rates for the constant work rate tests were 115 ± 36 W (90% LT) and 206 ± 56 W (50% {Delta}). In five subjects the {eph_162_mu107} response during recovery from heavy exercise was significantly better described by the model with the slow component emerging after an independent time delay (model II). The overall [HHb] profile beyond 60 s of recovery was variable among subjects and between exercise intensities (e.g. Fig. 1). However, in all subjects we observed a time delay before the decrease in [HHb]. The changes in total haemoglobin ([THb]) and oxyhaemoglobin concentration ([HbO2]) are shown in Fig. 2. The [THb] responses suggest that there were changes in blood volume under the NIRS probe during the recovery period. However, these are expected to have minor effects on the [HHb] signal (Ferrari et al. 1997; Grassi et al. 2003).



View larger version (36K):
[in this window]
[in a new window]
 
Figure 2.  Representative responses of total haemoglobin concentration (•, [THb]) and oxyhemoglobin concentration ({circ}, [HbO2])
The data shown for top, middle and bottom panels are, respectively, from subjects depicted in Fig. 1. Although there are substantial changes in [THb], these are expected to have less effect on the [HHb] response, as previously shown (Grassi et al. 2003). Accordingly, in general the temporal profile of [HbO2] resembles that of [THb], further suggesting that changes in blood volume affect primarily the [HbO2] signal.

 
Regarding the temporal profile of {eph_162_mu108}, all recovery responses (for moderate and heavy exercise) were better described by the model which included an initial component, indicating that the primary response of the kinetics of muscle capillary blood flow during recovery was biphasic. The recovery responses of estimated {eph_162_mu109} and [HHb] of a representative subject are shown in Figs 3 (90% LT) and 4 (50% {Delta}). The kinetic parameters of [HHb] and {eph_162_mu110} are shown in Table 1. The TD of [HHb] was shorter and {tau}[HHb] was longer for recovery from moderate compared to heavy exercise (P < 0.05), while the MRT of [HHb] (TD + {tau}) was significantly faster during recovery from moderate exercise (P < 0.01). The results for the overall kinetics of {eph_162_mu111}, [HHb] and {eph_162_mu112} are shown in Fig. 5. For moderate exercise, the overall kinetics of {eph_162_mu113} (MRT) were significantly slower than the kinetics of {eph_162_mu114} (as {tau}P) and [HHb] (as MRT). In addition, [HHb] kinetics were faster than the {eph_162_mu115} kinetics. During recovery from heavy exercise, MRT-{eph_162_mu116} was significantly slower than the ‘lumped’ time constant of {eph_162_mu117} ({tau}Pb), while {tau}Pb of {eph_162_mu118} did not differ significantly from the MRT of [HHb]. There was a significant correlation between MRT-{eph_162_mu119} and {eph_162_mu120} (r = 0.79, P < 0.01; Fig. 6). However, the data points were widely scattered, indicating a dissociation between the overall recovery kinetics of {eph_162_mu121} and muscle {eph_162_mu122}.


View this table:
[in this window]
[in a new window]
 
Table 1. Kinetic parameters of {eph_162_mu232} and [HHb] for recovery from moderate and heavy exercise
 


View larger version (13K):
[in this window]
[in a new window]
 
Figure 5.  Overall kinetics of deoxyhaemoglobin concentration ([HHb]), pulmonary oxygen uptake {eph_162_mu216} and muscle capillary blood flow {eph_162_mu217}
MRT, mean response time. The data for pulmonary {eph_162_mu218} represent the time constant of the primary component ({tau}P) of the response (eqn (2)). For heavy exercise, {tau}P was obtained by fitting the {eph_162_mu219} response over the time window corresponding to the initial and primary components of {eph_162_mu220} during recovery (i.e. from end of exercise to TDS of {eph_162_mu221}). The values shown are the means ± S.D. * Significantly different from MRT-{eph_162_mu222} at the same exercise intensity (P < 0.05); {dagger} significantly different from {eph_162_mu223} at the same exercise intensity (P < 0.05); {ddagger} significantly different from heavy exercise value (P < 0.05); § P = 0.068 compared with heavy exercise value.

 


View larger version (16K):
[in this window]
[in a new window]
 
Figure 6.  Relationship between the mean response time of estimated muscle capillary blood flow {eph_162_mu224} and time constant of pulmonary {eph_162_mu225} ( {eph_162_mu226} ; primary component) for moderate and heavy exercise
Dashed line, regression line; continuous line, line of identity. It is reasonable to interpret this relationship as muscle metabolism {eph_162_mu227} driving {eph_162_mu228}. However, the overall recovery kinetics of {eph_162_mu229} (MRT) appear to become progressively greater (slower) than the kinetics of {eph_162_mu230} for individuals with slower {eph_162_mu231}. See text and Fig. 4 for further details.

 
Comparison with onset of exercise

The results of selected parameters of {eph_162_mu123} kinetics following the onset of exercise obtained from a previous study with the same subjects (Ferreira et al. 2005b) are shown in Table 2. For moderate exercise, the contribution of phase 1 to the overall response (A'I/A'P) was smaller during recovery, but the remaining kinetic parameters were not significantly different from the onset kinetics. The overall effect was a slower MRT-{eph_162_mu124} during recovery. For heavy exercise, the A'I/A'P was similar for the onset and recovery from exercise. However, the {tau} of the initial (phase 1) component ({tau}I) of {eph_162_mu125} was longer during recovery, while both the duration of phase 1 (TDP) and {tau}P were greater for recovery, resulting in a slower MRT-{eph_162_mu126} during recovery compared to the onset of exercise. (N.B. This comparison should be made with caution because the temporal characteristics of {eph_162_mu127} were estimated and not directly measured for both onset and recovery. However, the {eph_162_mu128} kinetics were estimated for the same subjects and exercise bouts for both onset and recovery kinetics, which justifies comparison of the kinetic parameters.)


View this table:
[in this window]
[in a new window]
 
Table 2. Selected kinetic parameters of {eph_162_mu236} following the onset of moderate and heavy exercise
 

    Discussion
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
This investigation is the first, to our knowledge, specifically to address the temporal kinetic characteristics of [HHb] and muscle capillary blood flow (estimated) during recovery from exercise. In addition, we examined the association between the recovery kinetics of estimated {eph_162_mu129} and {eph_162_mu130}. There are two principal new findings in the present study. First, the kinetics of {eph_162_mu131} were biphasic during recovery (similar to the characteristics following exercise onset), with an initial fast phase followed by a second slower response emerging ~15–30 s after the end of exercise. Second, the overall kinetics of {eph_162_mu132} were slower than those of {eph_162_mu133}, indicating a temporal dissociation between the recovery of O2 delivery and utilization. It is important to emphasize that {eph_162_mu134} kinetics were estimated by making two basic assumptions: (a) {eph_162_mu135}; and (b) [HHb](t) {propto} C(a – v)O2(t) (see Study limitations below). Therefore, some error will probably exist in the estimated {eph_162_mu136} kinetics when compared with the ‘true’ {eph_162_mu137} response, as discussed in detail elsewhere (Ferreira et al. 2005b).

Recovery kinetics of muscle blood flow

The temporal profile of {eph_162_mu138} during recovery from exercise has been investigated in several studies (Van Leeuwen et al. 1992; Lash, 1994; Leyk et al. 1999; Van Beekvelt et al. 2001), but the kinetic characteristics of {eph_162_mu139} remain equivocal. Although the recovery of blood flow has been described by exponential equations (similar to eqn (2); Shoemaker et al. 1994, 1997; Whipp et al. 1995), theses studies did not distinguish (or report) whether the response was better described by mono- or two-exponential models. Leyk et al. (1999) observed that leg blood flow (measured using Doppler ultrasound) during recovery from calf exercise did not show an initial fast response (i.e. the decline in {eph_162_mu140} was monoexponential). In contrast, in a similar exercise model the recovery of {eph_162_mu141} demonstrated an early rapid decrease after cessation of exercise (Van Leeuwen et al. 1992). However, the authors (Van Leeuwen et al. 1992) contended that the response was not universally exponential because in 50% of the subjects blood flow increased again for a very short period after the initial rapid decrease. In the present investigation, the estimated {eph_162_mu142} response was better described, in all subjects, by a two- (moderate exercise) or three-exponential model (heavy exercise), indicating that in the microcirculation the recovery of muscle blood flow was biphasic (Figs 3 and 4).

As stated in the Introduction, it is relevant to determine whether the recovery of {eph_162_mu143} is better described by first- or higher-order models because this implies the existence of single versus multiple controlling mechanism(s), respectively, for the variable studied. Therefore, this is an important step to clarify the mechanisms operating to sustain an elevated blood flow (and O2 delivery) during the recovery period. The initial decrease in {eph_162_mu144} (phase 1) observed in our study ({tau}I {approx} 6–9 s) was substantially faster than the estimated {eph_162_mu145} (~30 s) or the subsequent phase 2 of the {eph_162_mu146} response (~29 s for moderate exercise and ~48 s for heavy exercise). It has been suggested that the skeletal muscle vasculature is able to respond (dilate or constrict) very rapidly during the transitional phase (Hamann et al. 2004; Tschakovsky et al. 2004). However, these suggestions are based on investigations of on-responses to a stimulus, while few studies have focused on the time course of response to the removal of a stimulus. In this context, direct recordings of arteriolar diameter after muscle contraction have shown a ~10 s delay to the onset of arteriolar constriction (i.e. arteriolar diameter temporarily remained similar to the exercise level; Gorczynski & Duling, 1978; VanTeeffelen & Segal, 2000). These observations suggest that the initial component of recovery of {eph_162_mu147} observed in the present study was determined primarily by a mechanical factor. Consistent with this, after cessation of calf exercise, the recovery of {eph_162_mu148} (half-time) became progressively faster when contractions were performed in the supine, sitting and standing positions (Van Leeuwen et al. 1992). This speeding of the decline of {eph_162_mu149} under conditions of augmented venous hydrostatic pressure was attributed to removal of the muscle pump (Van Leeuwen et al. 1992). In the present study we examined transitions from loaded to ‘unloaded’ cycling, which could have diminished the effect of the muscle pump when going from exercise to recovery (Sheriff & Hakeman, 2001; Shiotani et al. 2002). However, venous pressure in superficial veins (Shiotani et al. 2002) may not reflect the profile of intramuscular vessels that are subjected to contraction-induced pressure swings that will have substantial effects on the muscle pump-generated blood flow (Laughlin, 1987). Therefore, based on the ~10 s delay to the onset of arteriolar constriction observed in microcirculatory preparations, we speculate that the initial decrease in {eph_162_mu150} during recovery observed herein was largely caused by the removal of the muscle pump effect in the transition to unloaded pedalling. However, this hypothesis must await confirmation with kinetic analysis under experimental conditions that can modify the potential contribution of the muscle pump (e.g. body position).

As mentioned above, the second phase of {eph_162_mu151} was appreciably slower than the initial response. Considering that the putative mechanical factor is mostly evident during phase 1, phase 2 of {eph_162_mu152} would probably be determined by the interaction among vasoactive substances of neural (e.g. noradrenaline), muscular (metabolic byproducts) and vascular (endothelium) origin. For instance, pharmacologic intervention studies have shown that blockade of adenosine (Kille & Klabunde, 1984), prostaglandins (Kilbom & Wennmalm, 1976) and nitric oxide (Radegran & Saltin, 1999) decreased the muscle blood flow during recovery from exercise. Clearly, the mechanisms operating to keep blood flow elevated during recovery remain to be further investigated, and our results emphasize the importance of distinguishing the effects of each putative mechanism on each phase of the {eph_162_mu153} response. The available evidence, however, suggests that the interaction among these factors (neural, muscular and vascular) will probably depend on the type (e.g. isometric or dynamic), intensity (moderate or heavy) and body position (supine versus upright) of exercise (reviewed by Bangsbo & Hellsten, 1998).

Temporal relationship between {eph_162_mu154} and {eph_162_mu155} during recovery

{eph_162_mu156} and {eph_162_mu157} achieve their baseline (pre-exercise) levels within 3–4 min after the end of moderate exercise (Shoemaker et al. 1994; Radegran & Saltin, 1999; van Beekvelt et al. 2001). After heavy exercise {eph_162_mu158} and {eph_162_mu159} remain elevated for a longer period of time (> 10–15 min; Bangsbo et al. 1991; Van Beekvelt et al. 2001). Accumulation of metabolites during the heavy exercise period has been considered as a possible mechanism for the elevated postexercise {eph_162_mu160} (Bangsbo & Hellsten, 1998), whereas the cause of elevated {eph_162_mu161} is more controversial (for reviews see Gaesser & Brooks, 1984; Bangsbo & Hellsten, 1998). Although several studies have investigated the dynamics of both {eph_162_mu162} and {eph_162_mu163} during recovery from exercise (e.g. Bangsbo et al. 1994; Radegran & Saltin, 1999; Van Beekvelt et al. 2001), the kinetic parameters of these variables have been determined in only a few studies (Shoemaker et al. 1994; Whipp et al. 1995). In the study of Shoemaker et al. (1994) the femoral artery blood flow kinetics were similar to the kinetics of {eph_162_mu164} (primary component) after moderate exercise. In the present study, the overall kinetics of {eph_162_mu165} were associated with (Fig. 4), but slower (~8 and 16 s for moderate and heavy exercise, respectively) than the estimated {eph_162_mu166} kinetics. The cause of the discrepancy between our results and those of Shoemaker et al. (1994) is not clear. However, the monoexponential model used by Shoemaker et al. (1994) for analysis of {eph_162_mu167} included a time delay, but only the time constant was reported. Thus, it is possible that the reported {tau} reflected primarily the kinetics of {eph_162_mu168} during phase 2. If true, their results would then be in agreement with the present results, where the kinetics of {eph_162_mu169} (primary component) and {tau}P for {eph_162_mu170} were not significantly different. Therefore, in agreement with results from other studies (van Beekvelt et al. 2001) our data indicate that after moderate exercise the overall kinetics of muscle blood flow were slower than the kinetics of {eph_162_mu171}.

Although not the main focus of our study, the implications of the relationship between {eph_162_mu172} and {eph_162_mu173} kinetics need to be considered. Previous studies have shown that the recovery of muscle metabolism was limited by O2 delivery in exercise-trained subjects (Haseler et al. 1999), possibly by a faster recovery of {eph_162_mu174} compared to the {eph_162_mu175} kinetics (Barstow et al. 1990). In contrast, in sedentary individuals O2 delivery was not the limiting factor to the recovery kinetics of muscle phosphocreatine {eph_162_mu176}, suggesting that the recovery of {eph_162_mu177} was not faster than the kinetics of {eph_162_mu178}. Importantly, the slower recovery dynamics of {eph_162_mu179} compared to {eph_162_mu180} can be considered indirect evidence that, in the population presently studied (healthy, physically active subjects), the recovery kinetics of {eph_162_mu181} were not limited by O2 availability. However, this hypothesis remains to be tested directly with interventions that would increase convective and/or diffusive O2 delivery during recovery from exercise.

It has been shown that after heavy exercise the recovery of {eph_162_mu182} was appreciably slower than that of {eph_162_mu183} (Bangsbo et al. 1994; Van Beekvelt et al. 2001). The difference between the MRT of the primary component of {eph_162_mu184} estimated herein and {eph_162_mu185} for heavy exercise (~16 s) was twice that for moderate exercise (~8 s). This greater temporal dissociation of the restoration of O2 delivery and utilization, compared to moderate exercise, suggests that there is further alteration in the balance between sympathetic vasoconstriction and locally released vasodilators after heavy exercise. This is consistent with the progressively greater ‘functional sympatholysis’ observed at higher exercise intensities (Ruble et al. 2002).

Recovery of muscle oxygenation in the microcirculation

The recovery kinetics of muscle oxygenation have been investigated in studies examining the dynamics of tissue O2 saturation (StO2) in humans by NIRS (Chance et al. 1992; McCully et al. 1994; Belardinelli et al. 1997; Hanada et al. 2000) and microvascular O2 pressure (Pmv,O2) in rat muscle by phosphorescence quenching (McDonough et al. 2001). The overall kinetics of StO2 after moderate cycling exercise ({tau} of 20–25 s; Chance et al. 1992) were, on average, slower than the [HHb] kinetics found in the present study (~17 s). However, changes in vascular volume (i.e. total haemoglobin concentration; [THb] = [HHb] + [HbO2]) under the NIRS probe can distort the kinetics of StO2 (StO2 = [HbO2]/[THb]) compared to the O2 extraction kinetics (MacDonald et al. 1999; for discussion see DeLorey et al. 2003; Grassi et al. 2003; Ferreira et al. 2005b). In contrast, [HHb] is less sensitive to changes in local vascular volume and will probably better reflect the dynamics of O2 extraction (Ferrari et al. 1997; Grassi et al. 2003), although the kinetics of [HHb] and O2 extraction have not been directly compared.

The analysis of individual parameters describing the kinetics of [HHb] showed that the TD to the recovery of [HHb] was ~5 s for moderate exercise and 11 s for heavy exercise (Table 1). After cessation of muscle contractions that did not elicit an increase in blood lactate concentration (relatively light–moderate intensity), a time delay to the onset of recovery of Pmv,O2 (~5 s) has also been observed (McDonough et al. 2001), which is very similar to the TD of [HHb] observed in our study. The overall Pmv,O2 kinetics consisted of two phases, with an initial faster phase lasting ~55 s followed by a sluggish Pmv,O2 response (McDonough et al. 2001). The [HHb] response could also be considered to have at least two phases: an initial steep decrease in [HHb] (after the delay-like phase) followed by minor, slower changes thereafter (Figs 1, 3, 4). This [HHb] profile is probably a consequence of the concurrent monoexponential decrease in {eph_162_mu186} and a biphasic recovery of {eph_162_mu187}, which is qualitatively similar to that observed following the onset of exercise (DeLorey et al. 2003; Grassi et al. 2003; Ferreira et al. 2005b).

One important aspect to consider is the effect of fibre type on the dynamics of skeletal muscle oxygenation. McDonough et al. (2004) demonstrated that Pmv,O2 recovery kinetics of a slow-twitch (ST) muscle (rat soleus muscle) were faster than those of the fast-twitch (FT) counterpart (peroneal muscle). This suggests that the recovery of {eph_162_mu188}, versus that of {eph_162_mu189}, were slower for ST compared to FT muscles, which is consistent with the greater endothelium-dependent vasodilatation of ST muscles (Woodman et al. 2001). This aspect cannot be addressed with the present data because of the mixed composition of human muscles and lack of biopsy information from our subjects. However, we acknowledge that the slight predominance of FT fibres in superficial areas of vastus lateralis muscle (Lexell et al. 1983) sampled by the NIRS will probably lead to a slower mean response time of [HHb] than deeper muscle areas with a greater percentage of ST fibres.

Asymmetry of {eph_162_mu190} kinetics during exercise onset and recovery

Symmetry of on–off responses is a characteristic of linear control systems (Lamarra, 1990). In this context, the adjustment of cardiovascular variables (e.g. cardiac output and {eph_162_mu191}) in the transitional phases of exercise appear to be characteristically non-linear. The recovery kinetics of cardiac output (Yoshida & Whipp, 1994), forearm- (Van Beekvelt et al. 2001) and thigh-muscle blood flow (Whipp et al. 1995) were slower than their respective kinetics following the onset of exercise. Similarly, in the present study the estimated overall kinetics of {eph_162_mu192} were slower during recovery compared to the onset of moderate and heavy exercise.

Assuming that a mechanical factor (muscle pump) is the primary mediator of phase 1 of {eph_162_mu193} following the onset and recovery from exercise, one might expect that phase 1 characteristics would differ between these conditions. Following the onset of exercise the muscle pump is imposed on a ‘non-dilated’ vascular bed, whereas during recovery the removal of the muscle pump takes place under vasodilated conditions. In addition, arterial and venous pressures may be different at the end of exercise compared to unloaded exercise, prior to exercise. Therefore, potential differences in the physical properties of the muscle vascular system could account, at least in part, for the non-linearities of on–off kinetics. Moreover, phase 2 kinetics of {eph_162_mu194} during recovery from heavy exercise were slower than the on-kinetics, suggesting that different mechanisms (possibly metabolites; see above) are involved in the regulation of blood flow following the onset and recovery from exercise. Collectively, it seems reasonable to consider that the haemodynamic adjustments following the onset and recovery from exercise are asymmetrical, from central (cardiac output; Yoshida & Whipp, 1994) to peripheral circulation (larger vessels; Whipp et al. 1995; Van Beekvelt et al. 2001) and the microcirculation (present study).

Study limitations

The limitations of the methods used in our study, including the controversy regarding any potential myoglobin influence on the NIRS signal (Tran et al. 1999), and the assumptions made to estimate {eph_162_mu195} were addressed in detail previously (Ferreira et al. 2005b). Briefly, we assumed that the kinetics of the primary component of pulmonary {eph_162_mu196} reflect the kinetics of muscle {eph_162_mu197}. Based on results from computer models (Barstow et al. 1990) and the close agreement between phosphocreatine {eph_162_mu198} and the primary component of {eph_162_mu199} kinetics during recovery from exercise (Rossiter et al. 2002), this assumption appears to be valid. Moreover, we used the kinetics of [HHb] as a proxy of O2 extraction. The [HHb] kinetics following the onset of exercise in humans (Grassi et al. 2003; Ferreira et al. 2005b) are similar to those of C(a – v)O2 measured in different studies (Grassi et al. 1996, 2002). Although these observations were made for responses following the onset of exercise, there is no evidence to suggest that the relationship between [HHb] and C(a – v)O2 during recovery would be disrupted. However, the kinetics of [HHb] and C(a – v)O2 have not been directly compared. Van Beekvelt et al. (2001) reported that after moderate exercise forearm C(a – v)O2 returned to baseline levels within 90 s, which is similar to the [HHb] response observed herein. In some subjects we observed undershoots in [HHb] during recovery that were qualitatively similar to the C(a – v)O2 measurements of Van Beekvelt et al. (2001) for recovery from moderate and heavy exercise. Therefore, it is reasonable to assume that during recovery from exercise, as for the onset, changes in [HHb] reflect the dynamics of O2 extraction in the microcirculation.

Conclusions

In summary, we have demonstrated that during recovery from exercise the estimated {eph_162_mu200} kinetics were biphasic, with an early rapid decrease followed by a slower phase that emerged after a time delay (~15–30 s). We have suggested that, as for the onset of exercise (Sheriff & Hakeman, 2001), this early rapid phase of {eph_162_mu201} during recovery may be largely determined by a muscle pump effect, but this hypothesis needs to be tested. After cessation of exercise the overall kinetics of {eph_162_mu202} were slower than after the onset of exercise, suggesting an asymmetry of the mechanisms determining the on- versus off-kinetics of blood flow. In this context, following the onset of exercise the estimated {eph_162_mu203} kinetics were similar to the {eph_162_mu204} kinetics (Ferreira et al. 2005b), while the overall recovery kinetics of {eph_162_mu205} were slower than the {eph_162_mu206} recovery kinetics, indicating that the mechanisms controlling {eph_162_mu207} kinetics during recovery from exercise are not temporally associated to O2 uptake.


    References
 Top
 Abstract
 Introduction
 Methods
 Results
 Discussion
 References
 
Bangsbo J, Gollnick PD, Graham TE & Saltin B (1991). Substrates for muscle glycogen synthesis in recovery from intense exercise in man. J Physiol 434, 423–440.[Abstract/Free Full Text]

Bangsbo J, Graham T, Johansen L & Saltin B (1994). Muscle lactate metabolism in recovery from intense exhaustive exercise: impact of light exercise. J Appl Physiol 77, 1890–1895.[Abstract/Free Full Text]

Bangsbo J & Hellsten Y (1998). Muscle blood flow and oxygen uptake in recovery from exercise. Acta Physiol Scand 162, 305–312.[CrossRef][Medline]

Barstow TJ, Ferreira LF, Lutjemeier BJ & Townsend DK (2004). Tissue oxygen saturation by NIRS during 60 and 100 rpm ramp tests. Med Sci Sports Exerc 36, S25.

Barstow TJ, Lamarra N & Whipp BJ (1990). Modulation of muscle and pulmonary O2 uptakes by circulatory dynamics during exercise. J Appl Physiol 68, 979–989.[Abstract/Free Full Text]

Barstow TJ & Mole PA (1991). Linear and nonlinear characteristics of oxygen uptake kinetics during heavy exercise. J Appl Physiol 71, 2099–2106.[Abstract/Free Full Text]

Behnke BJ, Barstow TJ, Kindig CA, McDonough P, Musch TI & Poole DC (2002). Dynamics of oxygen uptake following exercise onset in rat skeletal muscle. Respir Physiol Neurobiol 133, 229–239.[CrossRef][Medline]

Belardinelli R, Barstow TJ, Nguyen P & Wasserman K (1997). Skeletal muscle oxygenation and oxygen uptake kinetics following constant work rate exercise in chronic congestive heart failure. Am J Cardiol 80, 1319–1324.[CrossRef][Medline]

Chance B, Dait MT, Zhang C, Hamaoka T & Hagerman F (1992). Recovery from exercise-induced desaturation in the quadriceps muscles of elite competitive rowers. Am J Physiol 262, C766–C775.