- Open Access
Modelling the lactate response to short-term all out exercise
Dynamic Medicinevolume 6, Article number: 10 (2007)
The maximum post exercise blood lactate concentration (BLCmax) has been positively correlated with maximal short-term exercise (MSE) performance. However, the moment when BLCmax occurs (TBLCmax) is rather unpredictable and interpretation of BLC response to MSE is therefore difficult.
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 (k1) and disappearance (k2). The 4-parameter model includes BLC at MSE termination and amplitudes and rate constants of increase (A1, y1) and decrease (A2, y2) of post MSE-BLC.
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 A1 (7.1 ± 1.6 vs. 10.9 ± 2.0 mmol l-1) were lower (p < 0.05) in MSE10 than in MSE30. k1 (0.610 ± 0.119 vs. 0.505 ± 0.107 min-1), k2 (4.21 10-2 ± 1.06 10-2 vs. 2.45 10-2 ± 1.04 10-2 min-1), and A2 (-563.8 ± 370.8 vs. -1412.6 ± 868.8 mmol l-1), and y1 (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 y2 (0.41 10-2 ± 0.82 10-2 vs. 0.15 10-2 ± 0.42 10-2 min-1) was found.
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.
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 (BLCmax) 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 BLCmax occurs (TBLCmax) 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 . 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 (BLCend_Ex), the amplitude (A1) and rate constant (y1) of the increase and the amplitude (A2) and the rate constant (y2) of the decrease of the BLC (Eq.1).
Recently a 3-parameter bi-exponential model was put forward as a means of calculating the BLC response to MSE . The 3-parameter model was initially developed to model radioactive transformations . Later it was used as pharmacokinetic model . 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 (BLC0). 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 (k1) and the disappearance (k2) of lactate into and out of the blood compartment (Eq.2).
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.
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 , 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 10th minute, and every second minute up to the 20th 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, k1 and k2, and also the bi-exponential 4-parameter model with the constants A1, A2, y1 and y2. The goodness of the fits of the 3- and the 4-parameter model was compared using the F-test . The interrelationships between selected variables were examined using simple regression analysis. Statistical significance was indicated using an alpha level of 0.05.
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).
Goodness of fits, parameter estimations, TBLCmax and BLCmax
Both models adequately described the BLC response to MSE10 and MSE30, respectively (all R2 ≥ 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).
In both tests, there was no difference between k1 of the 3- and y1 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 k2 was higher (p < 0.05) than y2 (Tab. 3 and 4).
TBLCmax in MSE30 and BLCmax in MSE10 were higher (p < 0.05) using the 3- than the 4-parameter model (Tab. 3 and 4). BLCmax 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, BLCmax and TBLCmax in MSE10 than in MSE30 (Tab. 3). k1 and k2 were higher (p < 0.05) in MSE10 than in MSE30 (Tab. 3). A negative correlation (p < 0.05) between k2 and A (r = 0.51; p < 0.005; y = -0.006 x + 0.056) and k2 and BLCmax was found (Fig. 5).
Use of the 4-parameter model showed lower (p < 0.05) values of A1, and higher (p < 0.05) levels of A2 and y1 in MSE10 than in MSE30 (Tab. 4). No significant difference in y2 was found between MSE10 and MSE30 (Tab. 4).
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 BLCmax 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 BLCmax coincided with an increase of TBLCmax from approx. 5 min to 7 min. Differences between calculated and directly measured values of BLCmax of approximately 1 % reflect time intervals of less than ± 1 min of TBLCmax, which appears practically neglectable under the consideration of usual blood sampling intervals between 1 and 3 min. Both models show that the longer TBLCmax 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 . 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 BLCmax increases [23–25]. The k2-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 y2-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 A2 (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 A2 values not lower than -22.5 mmol l-1 and halftimes of the BLC decrease more or less equivalent with k2-values estimated with the 3-parameter model in the present study and others  (Tab. 1, 2, 3).
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 BLCmax and delays TBLCmax 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.
Boycott AE, Haldane JS: The effects of low atmospheric pressures on respiration. J Physiol (Lond). 1908, 37: 355-377.
Embden G, Lehnartz E, Hentschel H: Der zeitliche Verlauf der Milchsäurebildung bei Muskelkontraktion. Z Physiol Chem. 1927, 165: 255-278.
Hill AV, Long CNH, Lupton H: Muscular exercise, lactic acid, and the supply and utilisation of oxygen: parts IV-VI. Proc Roy Soc B. 1924, 97: 84-138.
Bang O: The lactate content of the blood during and after muscular exercise in man. Contributions to the physiology of severe muscular work. Edited by: Bang O, Boje O, Nielsen M. 1936, Skandinavisches Archiv für Physiologie (Acta Societatis Physiologicae Scandinavicae), 51-82.
Margaria R, Edwards HT: The removal of lactic acid from the body during recovery from muscular exercise. Am J Physiol. 1934, 107: 681-686.
Beneke R, Pollmann Ch, Bleif I, Leithäuser RM, Hütler M: How anaerobic is the Wingate Anaerobic Test for humans?. Eur J Appl Physiol. 2002, 87: 388-392. 10.1007/s00421-002-0622-4.
Beneke R, Hütler M, Jung M, Leithäuser RM: Lactate response to maximal short-term high intensity exercise in children, adolescents and adults. J Appl Physiol. 2005, 99: 499-504. 10.1152/japplphysiol.00062.2005.
di Prampero PE: Energetics of muscular exercise. Rev Physiol Biochem Pharmacol. 1981, 89: 144-222.
Freund H, Zouloumian P: Lactate after exercise in man. I. Evolution kinetics in arterial blood. Eur J Appl Physiol. 1981, 46: 121-133. 10.1007/BF00428865.
Freund H, Zouloumian P: Lactate after exercise in man. IV. Physiological observations and model predictions. Eur J Appl Physiol. 1981, 46: 161-176. 10.1007/BF00428868.
Zouloumian P, Freund H: Lactate after exercise in man. II Mathematical model. Eur J Appl Physiol. 1981, 46: 135-147. 10.1007/BF00428866.
Zouloumian P, Freund H: Lactate after exercise in man. III Properties of the Compartment Model. Eur J Appl Physiol. 1981, 46: 149-160. 10.1007/BF00428867.
Bateman H: The solution of a system of differential equations occurring in the theory of radioactive transformations. Proc Cambridge Phil Soc. 1910, 15: 423-427.
Dost FH: Grundlagen der Pharmakokinetik. 1968, Georg Thime Verlag, Stuttgart, 36-47.
Inbar O, Bar-Or O, Skinner JS: The Wingate Anaerobic Test. 1996, Human Kinetics, Champaign, USA, 8-94.
Motulsky HJ, Ransnas LA: Fitting curves to data using nonlinear regression: a practical and nonmathematical review. FASEB J. 1987, 1: 365-374.
Jacobs I, Tesch PA, Bar-Or O, Karlsson J, Dotan R: Lactate in human skeletal muscle after 10 and 30 s of supramaximal exercise. J Appl Physiol. 1983, 55: 365-367.
Baker J, Brown E, Hill G, Pillips G, Williams R, Davies B: Handgrip contribution to lactate production and leg power during high-intensity exercise. Med Sci Sports Exerc. 2002, 34: 1037-1040. 10.1097/00005768-200205001-00413.
Beneke R: Maximal lactate steady state concentration (MLSS): experimental and modelling approaches. Eur J Appl Physiol. 2003, 88: 361-369.
di Prampero PE, Peeters L, Margaria R: Alactic O2 debt and lactic acid production after exhausting exercise in man. J Appl Physiol. 1973, 34 (5): 628-632.
Jones NL, McCartney N, Graham T, Spriet LL, Kowalchuk JM, Heigenhauser GJF, Sutton JR: Muscle performance and metabolism in maximal isokinetic cycling at slow and fast speeds. J Appl Physiol. 1985, 59: 132-136.
Spriet LL, Soderlund K, Bergstrom M, Hultman E: Anaerobic energy release in skeletal muscle during electrical stimulation in men. J Appl Physiol. 1987, 62: 611-615. 10.1063/1.339788.
Freund H, Oyono-Enguelle S, Heitz A, Marbach J, Ott C, Zouloumian P, Lampert E: Work rate-dependent lactate kinetics after exercise in humans. J Appl Physiol. 1986, 61: 932-939.
Heck H: Laktat in der Leistungsdiagnostik. Wissenschaftliche Schriftenreihe des deutschen Sportbundes. 1990, Verlag Karl Hofmann, Schorndorf
McLellan TM, Skinner JS: Blood lactate removal during recovery related to the aerobic threshold. Int J Sports Med. 1982, 3: 224-229.
Bassett DR, Merrill PW, Nagle FR, Agre JC, Sampedro R: Rate of decline in blood lactate after cycling exercise in endurance-trained and -untrained subjects. J Appl Physiol. 1991, 70: 1816-1820.
Freund H, Oyono-Enguelle S, Heitz A, Marbach J, Ott C, Gartner M: Effect of exercise duration on lactate kinetics after short muscular exercise. Eur J Appl Physiol. 1989, 58: 534-542. 10.1007/BF02330709.
Oyono-Enguelle S, Gartner M, Marbach J, Heitz A, Ott C, Freund H: Comparison of arterial and venous blood lactate kinetics after short exercise. Int J Sports Med. 1989, 10: 16-24.
Messonnier L, Freund H, Denis C, Feasson L, Lacour JR: Effect of training on lactate kinetics parameters and their influence on short high-intensity exercise performance. Int J Sport Med. 2006, 27: 60-66. 10.1055/s-2005-837507.
McLoughlin P, McCaffrey N, Moynihan JB: Gentle exercise with a previously inactive muscle group hastens the decline of blood lactate concentration after strenuous exercise. Eur J Appl Physiol. 1991, 62: 274-278. 10.1007/BF00571552.
The authors would like to thank Oliver Grzegorski for his support for the acquisition of the experimental data.
The author(s) declare that they have no competing interests.
RB designed and coordinated the study, helped to analyze the data and drafted the manuscript. MDJ conducted the experiments, analyzed the data and helped to draft the manuscript. RML helped to develop the study design and to draft the manuscript. All authors read and approved the final manuscript.
All experiments were conducted at the Institute of Sports Medicine, Free University Berlin.
Masen D Jumah and Renate M Leithäuser contributed equally to this work.