Modelling of Blood Lactate Time-Courses During Exercise and/or the Subsequent Recovery: Limitations and Few Perspectives

Because lactate is an important metabolic intermediate and a signalling molecule between/within cells/organs, it appears essential to be able to describe the kinetics of this central molecule, during and/or after physical exercise. The present study aimed to confront three models and their approaches [Freund and co-workers (F&co), Beneke and co-workers (B&co), and Quittmann and co-workers (Q&co)] to investigate the lactate exchange (γ₁) and removal (γ₂) abilities (min⁻¹) during and/or after exercise. Nine healthy male subjects performed 3- and 6-min easy, moderate, and heavy exercise. Blood lactate concentration (BLC) was measured every 5 s over the entire period of exercise and recovery. Approaches differ depending on the domain in which the model is applied: considering exercise and part of the recovery (B&co and Q&co) or the entire period of recovery (F&co). The different approaches result in differing γ₁ and γ₂ values. Model fitting is closer to the experimental values following the method (model and approach) of F&co. Complementary analyses show that consideration of (i) exercise drastically impairs the quality of model fitting and therefore the γ₁ and γ₂ values and (ii) the entire period of recovery considerably improves the quality of fits and therefore of the γ₁ and γ₂ values. We conclude that (i) it is neither realistic nor reliable to take into account exercise and recovery in the same model and (ii) the longer the period of recovery studied, the better the quality of the γ₁ and γ₂ values.

Introduction

Muscle contraction requires energy that comes from adenosine triphosphate (ATP) hydrolysis. Because of extremely low levels of ATP stores, several metabolic pathways are activated to resynthesize it. Glycogenolysis and glycolysis are major components of these metabolic pathways. Their activation induces lactate production, which increases with exercise intensity (Stanley et al., 1985; MacRae et al., 1992). When this production cannot be balanced anymore by removal processes, lactate accumulates in muscle (Chwalbinska-Moneta et al., 1989; Juel et al., 1990) and blood (Stanley et al., 1985; MacRae et al., 1992; Bergman et al., 1999; Messonnier et al., 2013). Blood lactate profiles obtained during incremental (Davis et al., 1983; Heck et al., 1985) or constant load exercises (Heck et al., 1985; Beneke, 1995) have been well-described. Different indexes drawn from the blood lactate vs. work rate curve have shown to be closely related to performance in different sporting activities (Faude et al., 2009) or markers of exercise intensity during endurance training in athletes (Yu et al., 2012; Tran et al., 2014) and in patients (Messonnier et al., 2021).

Lactate is well-known as a metabolic intermediate between and within cells and organs (lactate shuttle concept) for cell fuelling (Brooks, 1986, 2000, 2009) inducing and contributing to metabolic flexibility of various cells, such as cancer cells (Brooks, 2018; Brooks et al., 2021). It can also serve as a cell-signalling molecule and has the potential for the regulation of gene expression and epigenetic modifications (Hashimoto et al., 2007; Brooks, 2018; Brooks et al., 2021). From that point of view, it appears essential to be able to describe the kinetics of this central molecule, during and/or after physical exercise.

Tracer studies constitute the reference method to determine lactate kinetics parameters, such as the lactate rates of appearance (R_a) and disappearance (R_d) (Brooks et al., 1991; Miller et al., 2002; Messonnier et al., 2013). However, to be fully applicable and reliable, tracer techniques require almost steady states in blood lactate and tracer concentrations and a near equilibrium in concentrations between the compartments of the body (e.g., between active muscles and blood). These requirements are not fulfilled during high-intensity exercise and its subsequent recovery. Therefore, alternative methods should be found and applied. Freund and co-workers (F&Co) described the blood lactate profiles during recovery following short high-intensity exercises using a biexponential time function referring to a two-compartment lactate distribution space (Freund and Gendry, 1978; Freund and Zouloumian, 1981a,b; Zouloumian and Freund, 1981a,b). The interest of this approach lies in the fact that it allows for the determination of two major components of lactate kinetics, namely, the lactate exchange (γ₁) and removal (γ₂) abilities during recovery (Freund et al., 1986; Chatel et al., 2016). Of note, γ₁ and γ₂ should not be mixed up with R_a and R_d, which refer to different concepts. Using this approach, it has been shown that lactate exchange and removal abilities decrease with exercise duration and intensity (Freund et al., 1986, 1989; Chatel et al., 2016). Therefore, these abilities refer to the physiological state of the subject at the end of the exercise. Furthermore, it has been shown that these abilities were (i) improved by endurance training (Messonnier et al., 2001), (ii) different according to physical ability profiles of athletes (Bret et al., 2003), and (iii) related to performance (Messonnier et al., 1997, 2002). Applications of this model have been proposed, allowing for the estimation of net lactate release rate, the net amount of lactate released, lactate disappearance rate, and lactate metabolic clearance rate during recovery (Bret et al., 2003; Messonnier et al., 2006; Chatel et al., 2016).

The model proposed by Freund and Zouloumian and used later by their successors involved a passive and almost complete recovery, meaning that measurements of blood lactate concentrations (BLCs) during recovery may last 60–90 min if the exercise was heavy. To bypass this inconvenience, Beneke and co-workers (B&co) (Beneke et al., 2005) and more recently Quittmann and co-workers (Q&co) (Quittmann et al., 2018) proposed adapted models and approaches allowing reductions of lactate collection time during recovery to 20 and 10 min, respectively. To obtain realistic predictions and especially a return to lifelike resting blood lactate values, their models start at the pre-exercise BLC and use only one amplitude. If this approach forces the models to return to realistic resting BLC (i.e., the pre-exercise BLC), this implies that their models apply to both exercise and recovery. While the approach proposed by F&co provides information concerning lactate exchange and removal abilities during recovery, the approaches proposed by Beneke et al. (2005) and Quittmann et al. (2018) suggest that the same lactate exchange and removal abilities prevail during both exercise and recovery. This contradicts the fact that these abilities are altered by exercise itself (Freund et al., 1986, 1989). Does it constitute a pitfall?

The aim of the present study was to assess the three approaches using the same set of experimental data and compare the results. The hypothesis is that the adapted models are less precise than the models and approaches used by F&co and less universally applicable regardless of the performed exercise.

Materials and Methods

The data set used in the present study had been obtained previously (Freund et al., 1989). The methods are repeated here for the convenience of the reader. The study was conducted according to the laws and standards at the time of the experiment.

Subjects

Nine healthy male subjects volunteered to participate in this study. Their age, weight, height, maximal oxygen uptake (VO_2max), and maximal aerobic power (power associated with VO_2max) were 21.0 ± 2.4 years old, 67.7 ± 6.4 kg, 175 ± 8 cm, 3.76 ± 0.51 L min⁻¹, and 284 ± 33 W, respectively. Prior to giving their written consent, all subjects were informed of the aim and potential risks or discomforts associated with the experiments.

Experimental Design

All exercise tests were performed on the same ergo cycle (Fleisch ergometer) at a constant pedalling frequency of 60 rpm. Prior to the experiments, subjects underwent a physical examination and an incremental exercise up to exhaustion to determine the maximal oxygen uptake (VO_2max) and the corresponding maximal aerobic power. On two occasions 2 weeks apart, subjects were tested in the morning 2 h after having had a light standard breakfast. Under local anaesthesia, an indwelling catheter was placed in the brachial artery. The tests were performed at the normal room temperature (21–23°C). Each experiment was performed in the following order: a rest of 30 min, an easy exercise intensity [(mean ± SEM) 1.81 ± 0.03 W·kg⁻¹; 39–50% VO_2max] followed by a recovery of at least 30 min, a moderate exercise intensity (2.53 ± 0.08 W·kg⁻¹; 52–67% VO_2max) with subsequent recovery of at least 60 min, and finally, a heavy exercise intensity (3.52 ± 0.17 W·kg⁻¹; 76–82% VO_2max) followed by a recovery period of at least 90 min. On the first occasion, five subjects were randomly assigned to 3-min exercises while the four remaining subjects performed 6-min exercises. On the second occasion (2 weeks later), subjects reiterated the protocol on a crossover basis, switching from 3- to 6-min or to 6- to 3-min exercises. Each exercise and its subsequent recovery were considered separately. Seven curves were missing (three after 3-min exercises and four after 6-min exercise) due to exercises not being performed.

Arterial Blood Sampling and Analysis

To avoid coagulation, the subjects were heparinized (100 IU·kg⁻¹ body mass). Throughout exercise and recovery, arterial blood was sampled via a catheter at a rate of 0.32 mL·min⁻¹ by means of a well-calibrated peristaltic pump. Arterial blood was analysed automatically for lactate concentration (Freund, 1967). Lactate concentration was determined on a continuous flow custom-made analyser using an enzymatic method where the oxidation of lactate to pyruvate is coupled in the presence of lactate dehydrogenase and NAD⁺ to NADH + H⁺ formation. The operations required for the biochemical reactions (such as mixing, dialyses, heating and additions of buffers, lactate dehydrogenase, and NAD⁺) were carried out automatically by the analyser. The changes of absorption were measured at 340 or 365 nm (5–10 nm light path). The electrical signals supplied by the colorimeter were converted in lactate concentrations by means of the two standard curves determined before and after each experiment. It is worth mentioning that the time necessary for the blood to flow through the catheter was precisely measured before each experiment. The response time of the analyser was also rigorously measured during its calibration. Both these times were accounted for the synchronisation of data with the experimental protocol. For further details, we refer the reader to previous publications (Freund, 1967; Freund et al., 1989). Total blood loss during a test amounted to 0.100–0.150 L. Signals from the biochemical analyser allowed recording of BLCs every 5 s. Experimental values of maximal blood lactate concentrations (BLC_max, mmol L⁻¹) and time to reach BLC_max (TBLC_max, min) were recorded.

Mathematical Analysis

Arterial blood lactate (exercise and/or recovery) curves were fitted by Equation 1 (Freund and co-workers), Equation 2 (Beneke and co-workers), and Equation 3 (Quittmann and co-workers) using an iterative non-linear regression technique (Kaleidagraph 3.6, Synergy Software, PA, USA).

For the convenience of the readers, the original equations have been reworded to contain consistent terminology over the three equations.

$[La](t)={[La]}_{(\text{completion})}+A_1(1-e^{-{\gamma }_{\text{1}}.\text{t}})+A_2(1-e^{-{\gamma }_2.\text{t}})(1)$

$[La](t)=[(A^{.}{\gamma }_{\text{1}})/({\gamma }_2-{\gamma }_1)].(e^{-{\gamma }_{\text{1}}.\text{t}}-e^{-{\gamma }_{\text{2}}.\text{t}})+{[La]}_{(0)}(2)$

$\begin{array}{ll}\left[La\right](t)={\left[La\right]}_{(0)}+A(1-e^{-{\gamma }_1.\text{t}})-A(1-e^{-{\gamma }_2.\text{t}})& (3)\end{array}$

where [La] is arterial lactate concentration (mmol·L⁻¹), La_(completion) is the BLC at exercise completion (onset of recovery, mmol·L⁻¹), [La]₍₀₎ is the BLC at the beginning of exercise (resting value, mmol·L⁻¹), t is time (min), A₁ and A₂ are the amplitudes of the experimental terms (mmol·L⁻¹) describing lactate appearance and disappearance, respectively, A is this amplitude of appearance and disappearance (difference from maximal and resting BLC, mmol·L⁻¹), and γ₁ and γ₂ are the velocity constants (min⁻¹) describing lactate appearance and disappearance, respectively. These two latter terms, common to the three approaches, are important because they represent dynamic dimensions describing the blood lactate kinetics. BLC_max and TBLC_max values have been calculated from the parameters of the fits.

During the first series of experiments (Series 1), the three equations were applied over the domain proposed by the respective authors: only during the recovery for F&co, the starting value of physical exercise and the recovery values over 20 or 10 min for B&co or Q&co, respectively. In other words, during this first series, the three equations were applied as performed by Freund et al. (1989), Beneke et al. (2005), and Quittmann et al. (2018).

During the second series of experiments (Series 2), the three equations were fitted to the experimental BLC curves over the entire period of recovery, e.g., the starting value of physical exercise for Equations 2, 3.

During the third series of experiments (Series 3), the equations of F&co were applied to three periods of recovery: to the first 10 min, to the first 20 min, and to the entire period of the recovery.

Statistical Analysis

Values are presented as means and SDs. Statistical analysis was processed using Jasp^® software (JASP Team, Version 0.12.2, Amsterdam, the Netherlands). After verification of variable homogeneity and normality by the Shapiro-Wilk test, a sphericity test was performed to compare the means using repeated-measures ANOVA. When the data respected sphericity, Holm's post-hoc test was applied. Otherwise, a Greenhouse-Keizer correction factor was applied. If the data did not follow normality, a row ANOVA was performed, with the application of a post-hoc Conover test. Correlations between sets of data have been performed. Statistical significance was set at p < 0.05.

Results

Series 1

Figure 1 reports three experimental curves of BLCs obtained in the same subject during 3-min exercises and the following recoveries at three different exercise intensities (easy, moderate, and heavy, respectively), and the model fits from Equations 1–3 following the initial recommendations of F&co, B&co, and Q&co (see methods). Graphically, it can be observed that applying a model during the recovery only (F&co method) describes more closely the recovery experimental values rather than taking into account the pre-exercise value and part of the recovery (B&co or Q&co methods).

FIGURE 1

Figure 1. Three experimental blood lactate concentration (La) curves obtained in the same subject during 3-min exercises and the following recoveries at three different exercise intensities [easy (A), moderate (B), and heavy (C), respectively] and the fits from Equations 1–3 applied as Freund and co-workers (F&co), Beneke and co-workers (B&co), and Quittmann and co-workers (Q&co) did (see methods).

Table 1 reports the mean values of the velocity constants (γ₁ and γ₂) along with the R² values of the fits on the experimental data. The repeated measures of ANOVAs show significant differences between the three approaches (p < 0.001) for all the considered parameters. The post hoc tests highlight significant differences between the data obtained by F&co and B&co (p < 0.001), and F&co and Q&co (p < 0.001). The mean values of γ₁ and γ₂ are not different between the methods applied by B&co and Q&co (Table 1). The coefficients of determination R² show significant differences (p < 0.001) between the three approaches, the R² obtained by F&co being the highest. Parameters obtained with the different approaches are on the other hand significantly correlated (Table 1). By focusing only on the 3- and 6-min heavy intensity exercises, all parameters are different (Table 1), and more importantly, correlations between parameters are not obvious. In that sense, it seems that the longer the exercise, the weaker the correlations (Table 1).

TABLE 1

Table 1. Mean values of the velocity constants γ₁ and γ₂ along with the R² values of the fits on the experimental data of Freund et al., 1989.

Table 2 reports the mean values recorded on the experimental data for BLC_max and TBLC_max and those obtained from the fits. As it can be seen in Table 2, BLC_max and TBLC_max are underestimated by the approaches of B&co and Q&co.

TABLE 2

Table 2. Mean values of maximal blood lactate concentration (BLC_max) and time to reach BLC_max (TBLC_max) obtained from the experimental data (ExpData) or predicted by the fits, some of the parameters (γ₁ and γ₂) being reported in Table 1.

Series 2

Figure 2 reports three curves obtained in the same subject during 3-min exercises and the following recoveries at three different exercise intensities (easy, moderate, and heavy, respectively) applying the three approaches over the entire period of recovery. Graphically, it can be observed that applying a model during the recovery only (F&co) describes more closely the experimental values than taking into account the pre-exercise value and part of the recovery (B&co or Q&co).

FIGURE 2

Figure 2. Three experimental blood lactate concentration (La) curves obtained in the same subject during 3-min exercises and the following recoveries at three different exercise intensities [easy (A), moderate (B), and heavy (C), respectively] and the fits from Equations 1–3 applied over the entire period of recovery (for F&co) and taking also into account the initial blood lactate concentration at the beginning of exercise (for B&co and Q&co).

Table 3 reports the mean values of the velocity constants, γ₁ and γ₂, along with the R² values of the fits on the experimental data. The repeated measures of ANOVA show significant differences between the three approaches (p < 0.001) for all the considered parameters. The post-hoc tests highlight significant differences between the data obtained by F&co and B&co (p < 0.001), and F&co and Q&co (p < 0.001). The mean values of γ₁ and γ₂ are not different between the methods applied by B&co and Q&co (Table 3). The coefficients of determination R² show significant differences (p < 0.001) between the three approaches, the R² obtained by the method of F&co being the highest. Parameters obtained with the different approaches are on the other hand significantly correlated (Table 3).

TABLE 3

Table 3. Mean values of the velocity constants γ₁ and γ₂ along with the R² values of the fits on the experimental data of Freund et al. (1989).

Table 4 reports the mean values recorded on the experimental data for BLC_max and TBLC_max and those obtained from the fits. As it can be seen in Table 4, BLC_max and TBLC_max are underestimated and overestimated, respectively, by the approaches of B&co and Q&co.

TABLE 4

Table 4. Mean values of maximal blood lactate concentration (BLC_max) and time to reach BLC_max (TBLC_max) obtained from the experimental data (ExpData) or predicted by the fits, some of the parameters (γ₁ and γ₂) being reported in Table 2.

Series 3

Table 5 reports the mean γ₁ and γ₂ values using the model of Freund et al. applied on the entire period of recovery, or only during the first 20 or 10 min of recovery. γ₁ decreases significantly with the shortening of the recovery duration taken into account. On the contrary, γ₂ increases significantly with the shortening of the recovery duration taken into account. Taking into account only the heavy exercises (of 3 or 6 min), the γ₁ and γ₂ values are not different but they are not correlated to each other, meaning that hierarchisation of results/subjects on the basis of their γ₁ and γ₂ values is drastically obscured.

TABLE 5

Table 5. Mean values of the velocity constants γ₁ and γ₂ along with the R² values of the fits on the experimental data of Freund et al. (1989).

Discussion

The present study aimed to test different models and approaches in depicting and interpreting blood lactate curves during unsteady-state exercise and/or the subsequent recovery. As a whole, it appears that applying models and approaches (i) during exercise and recovery or (ii) during recovery but only over a restricted part of it present pitfalls. In these conditions, the models do not fit adequately to the experimental curves, challenging the interpretation of the velocity constants describing lactate appearance in and disappearance from the blood. Applying models and approaches only during the recovery and over the entire period of recovery (until returning to resting blood lactate values) still appears as the most accurate way to determine the lactate exchange and removal abilities.

Different Approaches and Different Results

Considering the mean values of parameters obtained by Beneke et al. (2005) after application of their model over the exercise and 20 min of recovery in response to a 0.5-min Wingate test, a return to the near resting BLCs after 200 min of recovery can be predicted. This time-to-return for resting blood lactate values is unrealistic. For instance, this duration is much longer (by ~1 h) than the one (~140 min) predicted from data obtained by Maciejewski et al. (2013) using the F&co model applied on the entire period of recovery after 3-min all-out exercises, while the exercise was totally exhaustive and led to higher BLCs than in the former case (Beneke et al., 2005). This first comparison already casts doubt on the fact that the different models and approaches provide congruent results. To answer the question more clearly, our first analysis compared the different parameters obtained after fitting the models on the same experimental set of data (Freund et al., 1989).

In the conditions of applying the equation proposed by each group (i.e., over recovery for F&co, over exercise, and 20 or 10 min of recovery for B&co or Q&co, respectively), the values of the velocity constants are drastically different (Table 1). Figure 1 shows clearly that the equation proposed by F&co fits almost perfectly the recovery curve. On the other hand, the fits from B&co and Q&co equations are more often aside from the experimental data. The deviations are particularly important in the initial phase of the recovery, rendering the determination of the peak BLC during the recovery difficult with the two latter approaches. To illustrate this latter point, BLC_max and TBLC_max are different from the experimental values using the approaches of B&co and Q&co, but are not different using the approach of F&co (Table 2).

One may argue that even if the results are different, correlations still exist between parameters obtained from the different approaches (Table 1). However, it is important to keep in mind that blood lactate curves have been obtained in very different conditions (3- or 6-min, easy-, moderate-, or heavy-intensity exercises). Therefore, the question arises as to whether the correlations are still observable for a given exercise, rendering hierarchisation between subjects possible and accurate. Table 3 reports stratified data according to the performed exercises. As it can be seen, only very few correlations between parameters obtained using the different approaches are obtained.

From that point of view, the lactate exchange and removal abilities (namely, γ₁ and γ₂) obtained from the three approaches are not identical. Two main differences exist between the approaches of F&co, and B&co and Q&co. The two differences lie (i) on the consideration of exercise and recovery or just recovery and (ii) on the considered period of recovery (partial or complete recovery). Therefore, are the differences observed between the approaches related to the first, the second, or a combination of the two differences? The following analyses will try to answer the question.

Effect of Exercise on Blood Lactate Kinetics Parameters (Series 2)

Comparison of data obtained using the three equations on the entire period of recovery (and taking into account the starting exercise value for B&co and Q&co) allows to apprehend the effect of taking into account the exercise on the parameter values. Table 3 reports the data obtained using this paradigm and underlines very different γ₁ and γ₂ values between the approaches used by F&co and both those used by B&co and Q&co. The fits are also closer to the experimental data using the F&co approach as indicated by the R² values and the correspondence with BLC_max and TBLC_max (Table 4; Figure 2). Together, these results suggest that considering exercise is detrimental for good fitting and determination of the parameters of the blood lactate recovery kinetics. Actually, this is not so surprising, previous studies have shown that γ₁ and γ₂ values decrease with exercise duration and intensity (Freund et al., 1986, 1989). Consistent with the decrease of γ₁ as the exercise progresses, Juel et al. (1990) observed that during ~3-min maximal exercises, the net lactate release rate from the active muscle to the blood increased and then tended to plateau, suggesting that lactate exchanges between active muscles and blood are disturbed when maximal exercise progresses. Therefore, if γ₁ and γ₂ are apparently constant during recovery, it appears unrealistic to apply the same constant γ₁ and γ₂ values for exercise and recovery.

Another way to demonstrate the limits of taking into account exercise and recovery is that approaches, such as exercise and recovery, cannot be applied during and after long-lasting (i.e., 60 min) exercises, while the approach proposed by Freund and co-workers is still valuable (Freund et al., 1990).

Effect of Length of Considered Period During Recovery

We also suspect that the length of the considered period of recovery is critical for an accurate determination of γ₁ and γ₂. A first clue tending to confirm this inference was the lower R² values obtained using the approach of Q&co (who considered only the initial 10 min of recovery) compared to the one used by B&co (who considered the first 20 min of recovery). To further investigate the importance of the length of the considered recovery period, we performed two additional analyses. First, a comparison was performed between the data obtained in Series 1 (considering only 20 or 10 min of recovery) and Series 2 (considering the entire period of recovery) for 3- and 6-min heavy-intensity exercises using the approaches of B&co and Q&co (Tables 1, 3). According to the period of recovery taken into account, values of parameters changed. Second, a comparison was performed between the data obtained after exercise using the F&co approach applied (i) on the full length, (ii) over 20 min, or (iii) over 10 min of recovery (Series 3; Table 5). Clearly, γ₁ and γ₂ values are different according to the length of recovery taken into account (Table 5). The shorter the recovery, the lower the γ₁ and the higher the γ₂. These complementary analyses indicate that it is better to take into account the entire recovery period. Actually, as it can be seen in Figures for moderate or heavy exercise, there are important inflexion points on the curves beyond the 20 min of recovery. These points need to be taken into account for precise γ₂ determinations.

Limitations

In the present report, we criticised the fact that taking into account exercise and recovery are not very precise in so far as the dynamical parameters that drive the blood lactate curve (according to the models depicted here), namely, the lactate exchange and removal abilities, change during exercise. In that regard, applying constant values during exercise and recovery makes little sense. Actually, the same criticism could be made to the F&co model and approach. Indeed, the probability that γ₁ and γ₂ remain constant over the entire period of recovery is very low. Indeed, if it is obvious that at the onset of the recovery, γ₁ and γ₂ reflect the physiological state of the subject at that time, it is also obvious that at the end of the recovery, γ₁ and γ₂ are almost returned to resting conditions. Nevertheless, the fits of the model are very close to the experimental data (R² > 0.999) so that everything appears as if the lactate exchange and removal abilities remain constant all over the recovery.

Perspectives

Although theoretically not perfect, the model and approach proposed by F&co seem for the moment the most accurate to describe the blood lactate curve during recovery and determine accurately two important aspects of its kinetics: the lactate exchange and removal abilities (namely, γ₁ and γ₂, respectively). Therefore, we recommend using this model and its domain of application for the entire period of recovery.

In the present study, the models and approaches proposed by B&co and Q&co have been applied on a set of data obtained during and after 3- or 6-min exercises. In such circumstances, strong limitations of these models and approaches to describe the blood lactate curve and thus to accurately determine γ₁ and γ₂ appeared. However, we acknowledge that the impact of exercise on the fits and on the determination of γ₁ and γ₂ would have been much lower if very short exercises (typically 15 or 30 s) had been used. Therefore, we strongly encourage the future scholars who would like to use the B&co or Q&co models and approaches to strictly apply these models (i) only during and after very short exercises and importantly (ii) on the entire period of recovery. Indeed, whatever the model and approach are, a precise determination of the lactate removal ability (γ₂) requires a long period of recovery.

Conclusion

Taken together, caution should be taken concerning the interpretation of the velocity constants obtained from the models and approaches proposed by Beneke et al. (2005) and Quittmann et al. (2018) if exercise is longer than 30 s and if the recovery is incomplete. At present, it seems that the model and approach proposed by F&co, although tedious due to the long period of recovery studied, remains the most accurate way to describe the blood lactate recovery curve and determine the lactate exchange and removal abilities during recovery.

Data Availability Statement

The data analyzed in this study is subject to the following licences/restrictions: Data have been previously published. Original data have been obtained from the authors. Requests to access these datasets should be directed to Hubert Freund, hubert.freund@wanadoo.fr.

Ethics Statement

We reused data already published. At the time of the experiments, the study complied with the law of the country for studies involving human participants and with the declaration of Helsinki. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

LM and HF designed and conducted this study and wrote the first draft of the article. RD, MG, MC, and DB analyzed data. All authors critically reviewed the draft and approved the final version for publication.

Conflict of Interest

LM receives consulting fees from Abbott Diabetes Care. All are outside the submitted work.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors would like to thank Julia Greenfield for her helpful English reviewing expertise.

References

Beneke, R. (1995). Anaerobic threshold, individual anaerobic threshold, and maximal lactate steady state in rowing. Med. Sci Sports Exerc. 27, 863–867. doi: 10.1249/00005768-199506000-00010

PubMed Abstract | CrossRef Full Text | Google Scholar

Beneke, R., Hutler, M., Jung, M., and Leithauser, R. M. (2005). Modeling the blood lactate kinetics at maximal short-term exercise conditions in children, adolescents, and adults. J. Appl. Physiol. 99, 499–504. doi: 10.1152/japplphysiol.00062.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergman, B. C., Wolfel, E. E., Butterfield, G. E., Lopaschuk, G. D., Casazza, G. A., Horning, M. A., et al. (1999). Active muscle and whole body lactate kinetics after endurance training in men. J. Appl. Physiol. 87, 1684–1696. doi: 10.1152/jappl.1999.87.5.1684

PubMed Abstract | CrossRef Full Text | Google Scholar

Bret, C., Messonnier, L., Nouck Nouck, J. M., Freund, H., Dufour, A. B., and Lacour, J. R. (2003). Differences in lactate exchange and removal abilities in athletes specialised in different track running events (100 to 1500 m). Int. J. Sports Med. 24, 108–113. doi: 10.1055/s-2003-38201

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. A. (1986). The lactate shuttle during exercise and recovery. Med. Sci. Sports Exerc. 18, 360–368. doi: 10.1249/00005768-198606000-00019

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. A. (2000). Intra- and extra-cellular lactate shuttles. Med. Sci. Sports Exerc. 32, 790–799. doi: 10.1097/00005768-200004000-00011

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. A. (2009). Cell-cell and intracellular lactate shuttles. J. Physiol. 587, 5591–5600. doi: 10.1113/jphysiol.2009.178350

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. A. (2018). The science and translation of lactate shuttle theory. Cell Metab. 27, 757–785. doi: 10.1016/j.cmet.2018.03.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. A., Arevalo, J. A., Osmond, A. D., Leija, R. G., Curl, C. C., and Tovar, A. P. (2021). Lactate in contemporary biology: a phoenix risen. J. Physiol. doi: 10.1113/JP280955. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

Brooks, G. A., Butterfield, G. E., Wolfe, R. R., Groves, B. M., Mazzeo, R. S., Sutton, J. R., et al. (1991). Decreased reliance on lactate during exercise after acclimatization to 4,300 m. J. Appl. Physiol. (1985) 71, 333–341. doi: 10.1152/jappl.1991.71.1.333

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatel, B., Bret, C., Edouard, P., Oullion, R., Freund, H., and Messonnier, L. A. (2016). Lactate recovery kinetics in response to high-intensity exercises. Eur. J. Appl. Physiol. 116, 1455–1465. doi: 10.1007/s00421-016-3420-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chwalbinska-Moneta, J., Robergs, R. A., Costill, D. L., and Fink, W. J. (1989). Threshold for muscle lactate accumulation during progressive exercise. J. Appl. Physiol. 66, 2710–2716. doi: 10.1152/jappl.1989.66.6.2710

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, H. A., Bassett, J., Hughes, P., and Gass, G. C. (1983). Anaerobic threshold and lactate turnpoint. Eur. J. Appl. Physiol. Occup. Physiol. 50, 383–392. doi: 10.1007/BF00423244

PubMed Abstract | CrossRef Full Text | Google Scholar

Faude, O., Kindermann, W., and Meyer, T. (2009). Lactate threshold concepts: how valid are they? Sports Med. 39, 469–490. doi: 10.2165/00007256-200939060-00003

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, H. (1967). Automatic continuous determination of pyruvic acid: description of the method and its application to simultaneous determination of blood pyruvic acid, blood lactic aicd and blood sugar. Ann. Biol. Clin. 25, 421–437.

PubMed Abstract | Google Scholar

Freund, H., and Gendry, P. (1978). Lactate kinetics after short strenuous exercise in man. Eur. J Appl. Physiol. Occup. Physiol. 39, 123–135. doi: 10.1007/BF00421717

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, H., Oyono-Enguelle, S., Heitz, A., Marbach, J., Ott, C., and Gartner, M. (1989). Effect of exercise duration on lactate kinetics after short muscular exercise. Eur. J. Appl. Physiol. Occup. Physiol. 58, 534–542. doi: 10.1007/BF02330709

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, H., Oyono-Enguelle, S., Heitz, A., Marbach, J., Ott, C., Zouloumian, P., et al. (1986). Work rate-dependent lactate kinetics after exercise in humans. J. Appl. Physiol. 61, 932–939. doi: 10.1152/jappl.1986.61.3.932

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, H., Oyono-Enguelle, S., Heitz, A., Ott, C., Marbach, J., Gartner, M., et al. (1990). Comparative lactate kinetics after short and prolonged submaximal exercise. Int. J. Sports Med. 11, 284–288. doi: 10.1055/s-2007-1024808

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, H., and Zouloumian, P. (1981a). Lactate after exercise in man: I. Evolution kinetics in arterial blood. Eur. J. Appl. Physiol. Occup. Physiol. 46, 121–133. doi: 10.1007/BF00428865

PubMed Abstract | CrossRef Full Text | Google Scholar

Freund, H., and Zouloumian, P. (1981b). Lactate after exercise in man: IV. Physiological observations and model predictions. Eur. J. Appl. Physiol. Occup. Physiol. 46, 161–176. doi: 10.1007/BF00428868

PubMed Abstract | CrossRef Full Text | Google Scholar

Hashimoto, T., Hussien, R., Oommen, S., Gohil, K., and Brooks, G. A. (2007). Lactate sensitive transcription factor network in L6 cells: activation of MCT1 and mitochondrial biogenesis. FASEB J. 21, 2602–2612. doi: 10.1096/fj.07-8174com

PubMed Abstract | CrossRef Full Text | Google Scholar

Heck, H., Mader, A., Hess, G., Mucke, S., Muller, R., and Hollmann, W. (1985). Justification of the 4-mmol/l lactate threshold. Int. J. Sport Med. 6, 117–130. doi: 10.1055/s-2008-1025824

PubMed Abstract | CrossRef Full Text | Google Scholar

Juel, C., Bangsbo, J., Graham, T., and Saltin, B. (1990). Lactate and potassium fluxes from human skeletal muscle during and after intense, dynamic, knee extensor exercise. Acta Physiol. Scand. 140, 147–159. doi: 10.1111/j.1748-1716.1990.tb08986.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Maciejewski, H., Bourdin, M., Lacour, J. R., Denis, C., Moyen, B., and Messonnier, L. (2013). Lactate accumulation in response to supramaximal exercise in rowers. Scand. J. Med. Sci. Sports. 23, 585–592. doi: 10.1111/j.1600-0838.2011.01423.x

PubMed Abstract | CrossRef Full Text | Google Scholar

MacRae, H. S., Dennis, S. C., Bosch, A. N., and Noakes, T. D. (1992). Effects of training on lactate production and removal during progressive exercise in humans. J. Appl. Physiol. 72, 1649–1656. doi: 10.1152/jappl.1992.72.5.1649

PubMed Abstract | CrossRef Full Text | Google Scholar

Messonnier, L., Freund, H., Bourdin, M., Belli, A., and Lacour, J. R. (1997). Lactate exchange and removal abilities in rowing performance. Med. Sci. Sports Exerc. 29, 396–401. doi: 10.1097/00005768-199703000-00016

PubMed Abstract | CrossRef Full Text | Google Scholar

Messonnier, L., Freund, H., Denis, C., Dormois, D., Dufour, A. B., and Lacour, J. R. (2002). Time to exhaustion at VO(2)max is related to the lactate exchange and removal abilities. Int. J. Sports Med. 23, 433–438. doi: 10.1055/s-2002-33740

PubMed Abstract | CrossRef Full Text | Google Scholar

Messonnier, L., Freund, H., Denis, C., Feasson, L., and Lacour, J. R. (2006). Effects of training on lactate kinetics parameters and their influence on short high-intensity exercise performance. Int. J. Sports Med. 27, 60–66. doi: 10.1055/s-2005-837507

PubMed Abstract | CrossRef Full Text | Google Scholar

Messonnier, L., Freund, H., Feasson, L., Prieur, F., Castells, J., Denis, C., et al. (2001). Blood lactate exchange and removal abilities after relative high-intensity exercise: effects of training in normoxia and hypoxia. Eur. J. Appl. Physiol. 84, 403–412. doi: 10.1007/s004210000378

PubMed Abstract | CrossRef Full Text | Google Scholar

Messonnier, L. A., Emhoff, C. A., Fattor, J. A., Horning, M. A., Carlson, T. J., and Brooks, G. A. (2013). Lactate kinetics at the lactate threshold in trained and untrained men. J. Appl. Physiol. 114, 1593–1602. doi: 10.1152/japplphysiol.00043.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Messonnier, L. A., Riccetti, M., Chatel, B., Galacteros, F., Gellen, B., Rupp, T., et al. (2021). How to implement endurance exercise training in sickle cell disease. Haematologica 106, 1476–1479. doi: 10.3324/haematol.2020.267047

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, B. F., Fattor, J. A., Jacobs, K. A., Horning, M. A., Navazio, F., Lindinger, M., et al. (2002). Lactate and glucose interactions during rest and exercise in men: effect of exogenous lactate infusion. J. Physiol. 544, 963–975. doi: 10.1113/jphysiol.2002.027128

PubMed Abstract | CrossRef Full Text | Google Scholar

Quittmann, O. J., Abel, T., Zeller, S., Foitschik, T., and Struder, H. K. (2018). Lactate kinetics in handcycling under various exercise modalities and their relationship to performance measures in able-bodied participants. Eur. J. Appl. Physiol. 118, 1493–1505. doi: 10.1007/s00421-018-3879-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Stanley, W. C., Gertz, E. W., Wisneski, J. A., Morris, D. L., Neese, R. A., and Brooks, G. A. (1985). Systemic lactate kinetics during graded exercise in man. Am J Physiol. 249, E595–602. doi: 10.1152/ajpendo.1985.249.6.E595

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, J., Rice, A. J., Main, L. C., and Gastin, P. B. (2014). Development and implementation of a novel measure for quantifying training loads in rowing: the T2minute method. J Strength Cond. Res. 28, 1172–1180. doi: 10.1519/JSC.0000000000000248

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Chen, X., Zhu, W., and Cao, C. (2012). A quasi-experimental study of Chinese top-level speed skaters' training load: threshold versus polarized model. Int. J. Sports Physiol. Perform. 7, 103–112. doi: 10.1123/ijspp.7.2.103

PubMed Abstract | CrossRef Full Text | Google Scholar

Zouloumian, P., and Freund, H. (1981a). Lactate after exercise in Man: II. Mathematical model. Eur. J. Appl. Physiol. Occup. Physiol. 46, 135–147. doi: 10.1007/BF00428866

PubMed Abstract | CrossRef Full Text | Google Scholar

Zouloumian, P., and Freund, H. (1981b). Lactate after exercise in man: III. Properties of the compartment model. Eur. J. Appl. Physiol. Occup. Physiol. 46, 149–160. doi: 10.1007/BF00428867

PubMed Abstract | CrossRef Full Text | Google Scholar