Modelling the lactate response to short-term all out exercise

Abstract

Background

The maximum post exercise blood lactate concentration (BLC_max) has been positively correlated with maximal short-term exercise (MSE) performance. However, the moment when BLC_max occurs (TBLC_max) is rather unpredictable and interpretation of BLC response to MSE is therefore difficult.

Methods

We compared a 3- and a 4-parameter model for the analysis of the dynamics of BLC response to MSEs lasting 10 (MSE10) and 30 s (MSE30) in eleven males (24.6 ± 2.3 yrs; 182.4 ± 6.8 cm; 75.1 ± 9.4 kg). The 3-parameter model uses BLC at MSE-start, extra-vascular increase (A) and rate constants of BLC appearance (k₁) and disappearance (k₂). The 4-parameter model includes BLC at MSE termination and amplitudes and rate constants of increase (A₁, y₁) and decrease (A₂, y₂) of post MSE-BLC.

Results

Both models consistently explained 93.69 % or more of the variance of individual BLC responses. Reduction of the number of parameters decreased (p < 0.05) the goodness of the fit in every MSE10 and in 3 MSE30. A (9.1 ± 2.1 vs. 15.3 ± 2.1 mmol l^-1) and A₁ (7.1 ± 1.6 vs. 10.9 ± 2.0 mmol l^-1) were lower (p < 0.05) in MSE10 than in MSE30. k₁ (0.610 ± 0.119 vs. 0.505 ± 0.107 min^-1), k₂ (4.21 10^-2 ± 1.06 10^-2 vs. 2.45 10^-2 ± 1.04 10^-2 min^-1), and A₂ (-563.8 ± 370.8 vs. -1412.6 ± 868.8 mmol l^-1), and y₁ (0.579 ± 0.137 vs. 0.489 ± 0.076 min^-1) were higher (p < 0.05) in MSE10 than in MSE30. No corresponding difference in y₂ (0.41 10^-2 ± 0.82 10^-2 vs. 0.15 10^-2 ± 0.42 10^-2 min^-1) was found.

Conclusion

The 3-parameter model estimates of lactate appearance and disappearance were sensitive to differences in test duration and support an interrelation between BLC level and halftime of lactate elimination previously found. The 4-parameter model results support the 3-parameter model findings about lactate appearance; however, parameter estimates for lactate disappearance were unrealistic in the 4-parameter model. The 3-parameter model provides useful information about the dynamics of the lactate response to MSE.

Background

In the early part of the last century, lactate was identified as an indicator of glycolytic activity [1–3]. Soon afterwards, it was observed that blood lactate concentration (BLC) continues to increase for a significant period after termination of maximal short-term high-intensity exercise (MSE) before it starts to decrease [4, 5]. The maximum of the post exercise blood lactate concentration (BLC_max) has been positively correlated with the ability to tolerate high levels of BLC, required to perform well in MSE events [6–8]. An increase in lactate enables calculation of the energy derived from anaerobic metabolism via glycolysis [6, 8]. However, the moment when BLC_max occurs (TBLC_max) is rather unpredictable and interpretation of BLC response to MSE is difficult.

In the early 1980s, a bi-exponential 4-parameter model was developed to analyse the dynamics of lactate after exercise [9–12]. It describes the post exercise BLC based on two compartments, the working muscle and the non-muscular lactate space, and assumes that the intra-vascular lactate concentration represents the average lactate concentration in the non-muscular space [11]. The two exponential terms describe the post exercise flux of lactate from the muscle into the non-muscular space and the subsequent disappearance of lactate out of the non-muscular space. It includes the BLC at the end of exercise (BLC_{end_Ex}), the amplitude (A₁) and rate constant (y₁) of the increase and the amplitude (A₂) and the rate constant (y₂) of the decrease of the BLC (Eq.1).

$$\text{BLC(t)}={\text{BLC}}_{\text{end\_Ex}}+{\text{A}}_{\text{1}}\cdot \text{(1}-{\text{e}}^{-{\text{y}}_{\text{1}}\cdot \text{t}}\text{)}+{\text{A}}_{\text{2}}\cdot \text{(1}-{\text{e}}^{-{\text{y}}_{\text{2}}\cdot \text{t}}\text{)}$$

(1)

Recently a 3-parameter bi-exponential model was put forward as a means of calculating the BLC response to MSE [7]. The 3-parameter model was initially developed to model radioactive transformations [13]. Later it was used as pharmacokinetic model [14]. It describes concentrations in the one compartment model as a function of time with first-order invasion and first-order elimination. It requires a BLC value at the start of exercise (BLC₀). It approximates an increase of lactate (A) in the extra-vascular water space generated by exercise metabolism mainly in working muscles during MSE, which is equivalent to the subsequent decrease to the pre-exercise lactate concentration. It furthermore estimates two velocity constants describing the corresponding kinetics of the appearance (k₁) and the disappearance (k₂) of lactate into and out of the blood compartment (Eq.2).

$$BLC(t)=\frac{A\cdot k_1}{k_2-k_1}\cdot (e^{-k_1\cdot t}-e^{-k_2\cdot t})+BL{C}_0$$

(2)

Differentiation of both models allows for the determination of TBLC_max, insertion of TBLC_max in Eq.1 and Eq.2 give the corresponding BLC_max.

The aims of the present study were 1) to test whether the proposed bi-exponential 3-parameter model sufficiently describes the changes in BLC compared with the previously used 4-parameter model, and 2) to analyse how parameters are changed by the duration of the MSE.

Methods

Eleven male subjects (mean ± SD age: 24.6 ± 2.3 yrs; height: 182.4 ± 6.8 cm; body mass: 75.1 ± 9.4 kg; body mass index: 22.6 ± 2.5 kg m^-2; peak oxygen uptake: 4080 ± 228 ml min^-1) participated in the present study. All subjects were healthy non-smokers, physically active but not specifically trained. None of whom were receiving any pharmacological or specific dietetic treatment. Informed consent was obtained from all subjects after explanation of the nature and risks involved in participation in the experiments, which conformed to internationally accepted policy statements regarding the use of human subjects as approved by the local ethics committee.

Each subject performed two MSE-tests lasting 10 s (MSE10) and 30 s (MSE30) on a mechanically braked cycle ergometer (834 E, Monark) in a randomised order. All tests were carried out at similar times in the morning at least two hours after a light meal. There was a recovery period of one week between testing sessions. The subjects were instructed not to engage in strenuous activity during the day before an exercise test.

In accordance with accepted recommendations for anaerobic performance testing [15], subjects performed a standardised five minutes cycling warm-up at 50 W which included two sprints lasting three seconds performed at the end of the third and the fourth minute as coordinative preparation for MSE10 and MSE30 tests. After a further 10 minutes of rest, the subjects were instructed to pedal as fast as possible. A resistance corresponding to 7.5 % of the body weight was applied after an acceleration phase of three seconds. After termination of each test, the subjects were supervised during a 20 min rest period, where they maintained a seated position.

The BLC was determined from capillary blood samples drawn from the hyperaemic ear lobe. Samples were taken immediately before each test, within 15 s after each test, minute by minute up to the 10^th minute, and every second minute up to the 20^th minute post-test. Samples were haemolysed and analysed utilizing the enzymatic amperometric method (Ebio Plus, Eppendorf).

Data are reported as mean and standard deviation (mean ± SD). Differences within subjects were analysed using the paired t-test. A multiple nonlinear regression analysis was used for the approximation of the BLC time courses using the bi-exponential 3-parameter model to determine the constants A, k₁ and k₂, and also the bi-exponential 4-parameter model with the constants A₁, A₂, y₁ and y₂. The goodness of the fits of the 3- and the 4-parameter model was compared using the F-test [16]. The interrelationships between selected variables were examined using simple regression analysis. Statistical significance was indicated using an alpha level of 0.05.

Results

Performance and BLC

Average mechanical power output was higher (10.1 ± 1.0 vs. 8.3 ± 0.7 W kg^-1, p < 0.05), whilst average mechanical work was lower (100.6 ± 9.8 vs. 248.0 ± 21.5 J kg^-1, p < 0.05) in MSE10 than in MSE30. At any post MSE time, the BLC was consistently higher (p < 0.05) under MSE30- than under MSE10-conditions (Fig. 1, 2, 3, 4).

Figure 1

BLC, individual 3-parameter model approximations for MSE10.

Figure 2

BLC, individual 4-parameter model approximations for MSE10.

Figure 3

BLC, individual 3-parameter model approximations for MSE30.

Figure 4

BLC, individual 4-parameter model approximations for MSE30.

Goodness of fits, parameter estimations, TBLC_max and BLC_max

Both models adequately described the BLC response to MSE10 and MSE30, respectively (all R² ≥ 0.9369). The reduction of the number of parameters decreased the goodness of the fit in every MSE10 and in 3 tests under MSE30-conditions (Tab. 1 and 2 and Fig. 1, 2, 3, 4).

Table 1 3-parameter vs. 4-parameter model to describe the BLC-response to MSE10

Table 2 3-parameter vs. 4-parameter model to describe the BLC-response to MSE30

In both tests, there was no difference between k₁ of the 3- and y₁ of the 4-parameter model. Furthermore, all lactate appearance parameters were highly correlated between both models (r > 0.84, p < 0.001). Contrary to the latter in both tests k₂ was higher (p < 0.05) than y₂ (Tab. 3 and 4).

Table 3 Kinetics of the blood lactate response to given short-term maximal cycling tests based on the bi-exponential 3-parameter model

Table 4 Kinetics of the blood lactate response to given short-term maximal cycling tests based on the bi-exponential 4-parameter model

TBLC_max in MSE30 and BLC_max in MSE10 were higher (p < 0.05) using the 3- than the 4-parameter model (Tab. 3 and 4). BLC_max in MSE10 was also higher (p < 0.05) than the corresponding highest directly measured BLC of 8.4 ± 1.8 mmol l^-1. All maxima estimated with different models were highly correlated (all r > 0.91, p < 0.001).

Effect of test duration

The 3-parameter model revealed lower (p < 0.05) values of A, BLC_max and TBLC_max in MSE10 than in MSE30 (Tab. 3). k₁ and k₂ were higher (p < 0.05) in MSE10 than in MSE30 (Tab. 3). A negative correlation (p < 0.05) between k₂ and A (r = 0.51; p < 0.005; y = -0.006 x + 0.056) and k₂ and BLC_max was found (Fig. 5).

Figure 5

Interrelationship between BLC_max and k₂.

Use of the 4-parameter model showed lower (p < 0.05) values of A₁, and higher (p < 0.05) levels of A₂ and y₁ in MSE10 than in MSE30 (Tab. 4). No significant difference in y₂ was found between MSE10 and MSE30 (Tab. 4).

Discussion

Both MSE tests generated normal to high values of power output, work and BLC values for maximal exercise tests of corresponding durations in non-specifically trained healthy male subjects aged 20 to 30 years [6, 7, 15, 17].

Under these conditions, the bi-exponential 3-parameter and the 4-parameter model adequately described the changes in BLC as a cumulative effect of all related factors irrespective of whether any of these factors' specific behaviour may be described adequately using exponential models as well. Reduction of the number of parameters from 4 to 3 did clearly impair the goodness of the fit of the BLC under MSE10- but only to a minor extent under MSE30-conditions.

The observed differences in A and BLC_max between MSE10 and MSE30 clearly reproduced the frequent observation that after high intensity short term exercise the BLC is not only a result of exercise intensity but also of the duration of the event requiring net lactate production in the extra-vascular compartment. Higher values of A and BLC_max coincided with an increase of TBLC_max from approx. 5 min to 7 min. Differences between calculated and directly measured values of BLC_max of approximately 1 % reflect time intervals of less than ± 1 min of TBLC_max, which appears practically neglectable under the consideration of usual blood sampling intervals between 1 and 3 min. Both models show that the longer TBLC_max results from a combined decrease of the rate constants of lactate invasion and lactate evasion.

Neither the 3- nor the 4-parameter model does directly predict the real amount of intra-muscular lactate increase of any specific muscle. MSE10 and MSE30 leg-cycle ergometries are whole-body exercises involving significant contributions from lean tissue masses throughout the entire body [18]. Assuming a muscle mass of 40 % of total body mass and typical dilution factors for total body water and blood water content [7, 8, 19, 20], the A-values estimated for MSE10 and MSE30 do represent a net increase in muscular lactate of 16.9 mmol kg^-1 (MSE10) and 27.6 mmol kg^-1 (MSE30), equivalent to a rate of net increase in muscular lactate of 1.69 mmol kg^-1 s^-1 and 0.92 mmol kg^-1 s^-1. These values are well within the range of the magnitude of directly measured intra-muscular lactates and equivalents of approx. 1.18 kJ kg^-1 and 1.93 kJ kg^-1 of anaerobic lactic energy production found in previous studies that used comparable durations of maximal cycling or electrical stimulation [17, 21, 22].

Both models provided almost identical results about the invasion of lactate into the blood compartment. This was not the case with respect to the lactate elimination. The 3-parameter model supported previous results that the rate constant decreases if the BLC_max increases [23–25]. The k₂-values found in the present experiments using the 3-parameter model represent rate constants of lactate elimination equivalent to BLC half times similar to others determined previously [23–25]. Contrary, the 4-parameter model provided y₂-estimates equivalent to halftimes of the BLC decrease of approx. 1650 min under MSE10- and approx. 3450 min under MSE30-conditions combined with extreme values of A₂ (Tab. 1, 2 and 4). These extreme parameter estimations of the 4-parameter model seem to indicate that the use of a relatively short post exercise period combined with the limited number of data points were insufficient for the most successful application of the 4-parameter model, despite an excellent goodness of the fits.

The latter may be indirectly supported by other previous studies, which used the 4-parameter model with recovery periods of up to 120 minutes [10, 26–29]. Using the 4-parameter model they estimated A₂ values not lower than -22.5 mmol l^-1 and halftimes of the BLC decrease more or less equivalent with k₂-values estimated with the 3-parameter model in the present study and others [30] (Tab. 1, 2, 3).

Conclusion

The 3-parameter model and the 4-parameter model adequately described the changes in BLC under MSE conditions. Reduction of the number of parameters did impair the goodness of the fit of the BLC. However, under the given testing conditions the 3-parameter model seems to provide more realistic parameter estimations with respect to the elimination of lactate from the blood compartment even at relative short periods of blood sampling.

The proposed 3-parameter model seems to provide useful information about the dynamics of lactate in the blood and in the extra-vascular compartment and it is sensitive to changes in test duration.

Increase of the duration of MSE from 10 to 30 s increases BLC_max and delays TBLC_max from approx. 5 to approx. 7 min. The latter delay is a combined effect of decreased rate constants of lactate invasion and particularly lactate evasion from the blood compartment.