Open Access

The kinetics of lactate production and removal during whole-body exercise

Theoretical Biology and Medical Modelling20129:7

DOI: 10.1186/1742-4682-9-7

Received: 12 December 2011

Accepted: 13 March 2012

Published: 13 March 2012

Abstract

Background

Based on a literature review, the current study aimed to construct mathematical models of lactate production and removal in both muscles and blood during steady state and at varying intensities during whole-body exercise. In order to experimentally test the models in dynamic situations, a cross-country skier performed laboratory tests while treadmill roller skiing, from where work rate, aerobic power and blood lactate concentration were measured. A two-compartment simulation model for blood lactate production and removal was constructed.

Results

The simulated and experimental data differed less than 0.5 mmol/L both during steady state and varying sub-maximal intensities. However, the simulation model for lactate removal after high exercise intensities seems to require further examination.

Conclusions

Overall, the simulation models of lactate production and removal provide useful insight into the parameters that affect blood lactate response, and specifically how blood lactate concentration during practical training and testing in dynamical situations should be interpreted.

Keywords

Aerobic power Anaerobic power Blood lactate Cross-country skiing Muscle lactate pH

Background

The metabolic power in humans is based on the production and consumption of adenosine triphosphate (ATP). Despite the approximately 100-fold increase in ATP utilization from rest to maximal-intensity exercise, the energetic demands of the muscles are usually satisfied without depleting the intracellular ATP e.g., [13]. In this connection, three sources for ATP synthesis are available; First, ATP can be produced aerobically in the mitochondria by oxidative phosphorylation. Second, ATP can be produced by anaerobic synthesis due to glycolysis or glycogenolysis. Finally, ATP can be produced by phosphocreatine (PCr) break down to Creatine (Cr) (i.e., ADP + PCr gives ATP + Cr in the Creatine Kinease (CK) reaction) e.g., [13].

The rate of oxygen (O2) consumption can be set to the sum of 1) a constant rate (resting O2 consumption), 2) a rate due to unloaded body movements and 3) a rate proportional to the aerobic energy used to perform work. For moderate constant work rates, the aerobic power increases towards a steady state condition. The concept of maximal lactate steady state (MLSS), that is the highest intensity where a steady state lactate can be obtained, has been regarded as important for endurance performance e.g., [46]. For exercise intensities above MLSS, associated with sustained acidosis, a slow component delays the attainment of a steady state value and causes O2 uptake to increase to values greater than those predicted from aerobic steady state demands. For exercise intensities exceeding the maximal oxygen uptake, the steady state corresponds to the level that would be attained if it was possible to carry out the exercise under pure aerobic conditions [7]. Obviously, this virtual steady state is never reached as the increase in oxygen uptake ends when the maximal oxygen uptake is achieved.

When the rate of ATP production by oxidative sources becomes insufficient, high rates of glycolytic or glycogenolytic ATP production are required. The endpoint of glycolysis is pyruvate, which represents a metabolite that can be reduced to lactate or oxidized to CO2 or H2O. Thus, by increasing exercise intensities, the working muscles and various tissues produce more lactate and release it into the plasma. At the same time, the skeletal muscles, the heart, the liver and the kidney cortex remove lactate from the circulation, and lactate is suggested to act as an intermediate for the shuttling of carbohydrate from cells and tissue with relatively low oxidative capacity to cells and tissues with high oxidative capacity [811]. Thus, it is well established that the blood lactate concentration is the result of the production and the removal of lactate in the blood. During steady state sub-maximal exercise, when lactate production (influx) equals lactate removal (outflux), the lactate concentration in the lactate pool stays constant and the rate of oxygen consumption is the measure of the whole body energy expenditure regardless of the magnitude of lactate production and removal or the absolute blood lactate concentration. At exercise intensities above steady state, a rise in the concentration could be attributed to an increase in the rate of lactate production or result from a decrease in the rate of lactate removal.

Lactate itself does not lead to muscle fatigue at high exercise intensities, and fatigue is most likely a result of decreased muscle pH and the associated reactions e.g., [12, 13]. For the skeletal muscles, the pH is normally around 7.1 but can fall to 6.4 during heavy exercise [14], and numerous experiments suggest a negative relationship between decreased pH and the muscle contractile function [1520]. In the capillary blood pH is normally (in steady state conditions) around 7.45, but can fall to around 7.05 during heavy exercise [21]. Previous researchers have found a close relationship with a time delay between lactate concentrations in muscle and blood [22].

Because muscle groups often work unequally, heterogeneity with regard to lactate concentration is likely in muscles during dynamic exercise. A rigorous estimation would require application of several compartments with production, removal and exchange between compartments. As an approximation, a two-compartment model can be applied, where the muscles and other organs that can remove lactate (such as the heart, the liver and the kidney cortex) are regarded as one compartment and the blood space and other tissues are grouped into a second compartment [2325]. These studies indicate that the flux of lactate into a compartment depends upon the lactate gradient and its permeability. Furthermore, blood lactate recovery curves from muscular exercise can be described by a bi-exponential time function and a two-compartment model consisting of the previously worked muscle and the remainder lactate space. The time constants of the bi-exponential time function fitted to the arterial blood lactate recovery curves reflect the abilities to exchange lactate between the two-compartments and the ability to remove lactate from the total lactate space including the working muscle department [26]. The two time constants are found to decrease with work rate and duration of the preceding exercise [27, 28].

The maximum anaerobic energy that can be utilized is proportional to the sum of Cr and lactate that can be accumulated in the body. PCr is an energy buffer that supports the transient failure of other metabolic pathways to support ATP. The equilibrium constant of the CK reaction is around 20 and the slightest drop in ATP allows the reaction to proceed to ATP. Thus, the ATP concentration stays nearly constant until almost all the PCr is utilized. The PCr levels follow an exponential time course after changes in work rate before approaching a steady state condition at moderate exercise intensities. In such cases, a strong similarity has been reported for the time constants of the O2 kinetics and the PCr consumption [29]. However, for exercise intensities above lactate threshold, the anaerobic glycolytic energy supply becomes significant and the association between PCr and O2 rate has not yet been systematically reported. During recovery, the level of PCr must be recovered, the pH must be re-established and ADP removed. While the PCr recovery is mainly due to oxidative ATP synthesis, the PCr stores may be rebuilt by anaerobic glycolysis [3032].

The current study aimed to derive mathematical models for production and removal of lactate in the blood and muscles during dynamic whole-body exercise. As cross-country skiing is a whole body exercise where athletes both train and compete on varying terrain and at constantly varying speeds and work rates [3335] this locomotion was used for that purpose. Thus, the mathematical models were compared with experimental results from the laboratory during treadmill roller skiing. We hypothesized a mathematical two-compartment model of lactate production and removal to accurately predict blood lactate concentration during steady state and at varying exercise intensities.

Methods

Overall design

Initially, the current study derived mathematical simulation models of lactate production and removal in the blood and muscles by utilizing Mathematica 8 (Wolfram Research Inc., Champaign, IL, USA). Thereafter, the simulations were compared with experimental data from an elite skier performing laboratory tests while treadmill roller ski skating both during a steady state and at varying exercise intensities (see details below).

Steady state aerobic power

Work rate W on the treadmill is = m v v ˙ . 0 + μ m g C o s ( α ) v + m g S i n ( α ) v https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq1_HTML.gif, where v is the treadmill velocity, μ the coefficient of friction, m the mass of the skier and αSin(α) is the treadmill incline in radians. The power due to the change of kinetic energy m v v ˙ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq2_HTML.gif is zero on the treadmill since the velocity is constant. μmgCos(α)v is the power of roller friction, and mgSin(α)v is the power of gravity due to the inclination of the treadmill. We define Q max to be the maximal aerobic power and Q ̄ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq3_HTML.gif to be the steady state aerobic power that has been found while treadmill roller skiing. Q a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq4_HTML.gif is regarded as the virtual steady state aerobic power, with Q ̄ a = m o d M i n Q ̃ ̄ a , Q m a x https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq5_HTML.gif, where Min is the minimum function. The Min function is employed to ensure that the aerobic power does not exceed the maximum aerobic power. Below Q max the virtual steady state power Q ̃ ̄ a = Q ̄ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq6_HTML.gif is found to be linear with work rate for a given cycle rate and incline, and as a hypothesis we apply this linearity also for metabolic powers above Q max and find
Q ̃ ̄ a = m o d Q b . r e s t + Q u l ( f ) . u n l o a d e d + c 2 × Ā ( α ) . l o a d e d × W , Q ̄ a = m o d M i n Q ̃ ̄ a , Q m a x Q b = 80 J / s , Q u l ( f ) 111 J / s , c 2 = 5 . 8 , Ā ( α ) = d e f 0 . 92 ( 1 + 1 . 19 E x p ( - 71 . 8 α ) ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ1_HTML.gif
(1)

Q b is the metabolic power at rest (set to 80 J/s), Q un the metabolic power of unloaded movements (zero work rate) which is dependent of the cycle rate. We define the cycle rate to be constant in this article. Thus, Q un = 111 J/s here.

Thereafter, we let Q a be the aerobic power, Q ̃ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq7_HTML.gif the virtual aerobic power and Q a = m o d M i n ( Q m a x , Q ̃ a ( t ) ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq8_HTML.gif. A first order differential equation of aerobic power with the virtual aerobic power as input was used to account mathematically for the delay in aerobic power with a time lag during steady state work rate. Thus,
Q ̃ ˙ a ( t ) = m o d Q ̃ ̄ a - Q ̃ a ( t ) τ a , ( a ) , Q a = m o d M i n ( Q m a x , Q ̃ a ( t ) ) , ( b ) , τ a = 30 s , Q a ( t 0 ) = Q b https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ2_HTML.gif
(2)

The "dot" means time derivative, τ is a time parameter quantifying the time before the aerobic power reaches steady state during sub-maximal work rates. We uses τ = 30 s according to di Prampero [36].

Figure 1 shows the steady state aerobic power as function of the work rate for the G3 skiing technique.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig1_HTML.jpg
Figure 1

The aerobic power Q(W, α) as a function of work rate (W) for inclines of α = 0.05 (upper line) and α = 0.12 (lower line) during treadmill roller skiing using the skating technique. The curve fittings are based on least square fit to the data. Maximal metabolic rate (Q max = 1886 J/s) is represented by a straight horizontal line. The lactate threshold (QLT = 1755 J/s) is represented by the straight dotted horizontal line. ■: Experimental values for α = 0.05, ▲: Experimental values for α = 0.12.

Anaerobic power and lactate concentration in blood and muscle

Lactate concentration in the lactate pool (C(t)), i.e. the mass of lactate per unit volume of this pool, increases only if the rate of lactate appearance (influx) in the lactate pool is larger than the rate of lactate disappearance (outflux). The current study uses a modified version of Brooks [8] and determines the levels of lactate pool concentration by modeling both the influx and outflux streams of lactate and simulate the blood lactate concentration according to Moxnes and Hausken [37]. Here two pyruvate molecules are produced for each glucose or glycogen molecule during glycolysis/glycogenolysis. One molecule of pyruvate gives one molecule of lactate. The increase in glycogenolysis/glycolysis, due to increased exercise intensity, ends when a maximal rate of glycogenolysis/glycolysis is achieved [37]. Therefore, the rate of pyruvate appearance (P) has a least upper bound, which is denoted as P max . By neglecting rates due to changing pyruvate concentration in the plasma, the rate of lactate appearance in mmol/L is R a = P. However, as pyruvate can be oxidized in the mitochondria we set the influx of pyruvate into the mitochondria to be α0Tanh(β0, P)/(β0, α0 < 1, where α0 and β0 are to parameters fitted to the data and Tanh() is a function that accounts for saturation at high lactate concentrations according to Moxnes and Hausken [37]. Thus, the rate of lactate appearance is R a = P - α0Tanh(β0P)/β0 where we expect that α0 is around 1.

During severe exercise, glycogen re-synthesis by the liver is severely depressed. As a hypothesis, we forecast that the rate of lactate disappearance due to both glycogen re-synthesis and lactate oxidation is R d = d0 × (Tanh(χ C )/χ) × D(Q a ) × (Q max - Q a ), where d 0 and χ are two parameters that are fitted to the data. D(Q a (t)) is an unknown function that is monotonically increasing with the aerobic power. When aerobic power equals the maximum aerobic power (i.e. Q a = Q max ) no lactate disappearance takes place. Altogether, for a one compartment model of lactate we have
Ċ ( t ) = R a - R d = m o d P - α 0 T a n h ( β 0 P / ) . a p p e r a n c e / β 0 - d 0 × ( T a n h ( χ C ) / χ ) × D ( Q a ) × ( Q m a x - Q a ) . d i s a p p e r a n c e https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ3_HTML.gif
(3)
In equation (3) we need a model for the rate of pyruvate appearance (P). To our knowledge, no such model exists in the literature, and we hypothesize that this function is linear with work rate up to P max which is the least upper bound. This simple assumption might constitute a potential weakness in our argument. However, if we have experimental support this assumption also seems like the most reasonable one. Thus, we define a virtual steady state ( ( P ̃ ̄ ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq9_HTML.gif) analogous to the virtual aerobic steady state ( Q ̃ ̄ a ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq10_HTML.gif and hypothesize that
P ̃ ˙ ( t ) = m o d P ̃ ̄ - P ̃ ( t ) τ a n , τ a n = 10 s , P = m o d M i n ( P m a x , P ̃ ( t ) ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ4_HTML.gif
(4)

where τ an is the time constant for full activation of glycolysis/glycogenolysis during muscle contractions, set to 10 s [38]. P ̄ ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq11_HTML.gif is the steady state rate of pyruvate appearance. Note that the model in equation (3) applies only for the chosen type of exercise and a fixed concentration of glycogen in the body.

For aerobic power we assume that Q a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq4_HTML.gif is a linear function of work rate and forecast that the steady state rate of pyruvate appearance is approximately proportional to the virtual steady state power leading to
P ̃ ̄ ( t ) = m o d p 0 Q ̃ ̄ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ5_HTML.gif
(5)
Due to the rather fast response time in equation (4), an approximation is that ( t ) P ̃ ̄ = p 0 Q ̃ ̄ a ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq12_HTML.gif. For equation (3) this gives that
Ċ ( t ) p 0 Q ̃ ̄ a ( t ) - α 0 β 0 p 0 T a n h β 0 p 0 Q ̃ ̄ a ( t ) - d 0 × ( T a n h ( χ C ( t ) ) / χ ) × D ( Q a ( t ) ) × ( Q m a x - Q a ( t ) ) , C ( t 0 ) = 0 . 045 k g / m 3 = 0 . 5 m m o l / L https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ6_HTML.gif
(6)
Finally, a model for D(Q a ) is needed. As a hypothesis we set that
D ( Q a ) = m o d Q a - α 0 β p 0 T a n h ( β 0 p 0 Q a ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ7_HTML.gif
(7)

Here the constant of proportionality can be scaled by d0.

A general solution for constant aerobic power is feasible in (5)-(7) when we let T a n h ( χ C ( t ) ) / χ C ( t ) ( Q a = Q ̄ a = constant ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq13_HTML.gif [37].

Altogether, the steady state lactate concentration for constant aerobic power at intensities below the MLSS from equation (5)-(7) is
C ¯ = A r c T a n h ( p 0 χ d 0 ( Q m a x Q ¯ a ) ) / χ = 1 2 L n ( 1 + χ p 0 d 0 Q m a x ( 1 Q ¯ a / Q m a x ) 1 χ p 0 d 0 Q m a x ( 1 Q ¯ a / Q m a x ) ) / χ L i m C ¯ χ 0 = p 0 d 0 ( Q m a x Q ¯ a ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ8_HTML.gif
(8)
The maximum aerobic power that can be used for steady state concentration is given by
1 - p 0 d 0 ( Q m a x - Q ̄ a ) = 0 Q ̄ a = Q ̄ L T = Q m a x - p 0 χ d 0 Q ̄ L T Q m a x = 1 - p 0 χ d 0 Q m a x https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ9_HTML.gif
(9)

where Q ̄ L T = 0 . 93 Q m a x = MLSS =  , p 0 = 1 0 - 5 k g / ( m 3 s ) / ( J / s ) , d 0 = 7 . 241 0 - 8 / ( J / s ) 2 / s , χ = 0 . 95 m 3 / k g https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq14_HTML.gif and Q max = 1886 J/s.

According to equation (8) all steady state lactate concentrations can be achieved for steady state aerobic powers below Q ̄ L T https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq15_HTML.gif (MLSS). Furthermore, the steady state lactate concentration ( C ̄ ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq16_HTML.gif approaches infinity when the aerobic power approaches Q ̄ L T https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq15_HTML.gif. However, it has been discovered that not all steady state lactate levels are tolerated over time. This means that blood lactate levels above a certain level of exercise can be terminated before steady state is reached. The anaerobic power due to anaerobic glycogenolysis/glycolysis ( Q a n G ( t ) ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq17_HTML.gif can be calculated from the increase lactate concentration when using the relation from di Prampero and Ferretti [7] as
Q a n G ( t ) = m × λ × Ċ ( t ) , λ = 3 × 20 : J / ( k g m m o l / L ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ10_HTML.gif
(10)
Blood lactate concentration usually continues to rise a short period of time after exercise. The model in equation (6) does not capture this phenomenon since we only considered one compartment; the lactate pool. In the next step we consider the different compartments involved: working muscles, blood and other tissues such as the liver, kidney and heart. For exercise powers significantly above resting values, we regard a two-compartment model as sufficient: a) the blood compartment and b) muscles and other tissues as one compartment (denoted muscles in the rest of current manuscript). C b (t) is the lactate concentration in the blood and C m (t) the lactate concentration in the muscles. V b (t) and V m (t) are the volumes of muscles and blood, respectively. The total lactate pools volume is V = V b (t) + V m (t), set to 0.18 L per kg body mass as an approximation for the skier modelled here [39]. The muscle mass is set to be 10 kg, based on an iDexa scan of the skier, and muscle volume as 10 L. We assume that lactate moves between the muscle and blood compartment with some time dynamics. Thus, we propose the following model for the muscle's (Cm) and blood's (Cb) concentrations of lactate
C ˙ m ( t ) = ( V V m ) [ p 0 D ( Q ˜ ¯ a ( t ) d 0 ( T a n h ( χ C m ( t ) ) / χ ) D ( Q a ( t ) ) ( Q m a x Q a ) ] k 1 a C m ( t ) C b ( t ) V m C ˙ b ( t ) = k 1 a C m ( t ) C b ( t ) V b https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ11_HTML.gif
(11)
K1ais a parametric function that scales the movement of lactate into and out of the blood. The rate of lactate appearance in the blood is set proportional to the difference in lactate concentration between blood and muscles. To account for different rates of transport into or out of the blood we let K1abe dependent on C m (t) - C b (t). Estimates of these parameters were found by visual curve fitting. Visual curve fitting gives plausible values for the parameters based on plotting of experimental data that is compared with simulations. The parameters were sought to have biologically trustworthy numerical values and a least square fit of the data was performed to produce best fit estimates. The parameters employed here are
p 0 = 1 0 - 5 k g / ( m 3 s ) / ( J / s ) , d 0 = 7 . 241 0 - 8 / ( J / s ) 2 / s , χ = 0 . 95 m 3 / k g , α 0 = 0 . 9 1 / ( β 0 p 0 ) = 0 . 6 Q m a x , α 0 = 0 . 9 , k 1 a = 0 . 05 / s × L , V b = 4 L , V m = 10 L https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Equ12_HTML.gif
(12)

Experimental tests

The derived mathematical simulations are compared with experimental data from an elite crosscountry skier while roller skiing on a treadmill employing the skating G3 technique. The mass of the skier was m = 77.5 kg and the friction coefficient on the treadmill was μ = 0.024 in all tests. Equipment and procedures were similar to the studies by Sandbakk, Holmberg, Leirdal and Ettema [40, 41]. All treadmill tests were performed on a 6 × 3 m motor-driven treadmill (Bonte Technology, Zwolle, The Netherlands). Inclination and speed were calibrated using the Qualisys Pro Reflex system and the Qualisys Track Manager software (Qualisys AB, Gothenburg, Sweden). The treadmill belt consisted of a non-slip rubber surface that allowed the skier to use his own poles (pole length: 90% of body height) with special carbide tips. The skier used a pair of Swenor skating roller skis with standard wheels (Swenor Roller skis, Troesken, Norway) and the Rottefella binding system (Rottefella AS, Klokkartstua, Norway), and the roller skis were pre-warmed before each test through 20-min of roller skiing on the treadmill. The roller skis were tested for rolling friction force (Ff) before the test period, and the friction coefficient (μ) was determined by dividing Ff by the normal force (N) (μ = Ff · N-1). This was performed in a towing test on three subjects (70, 80 and 90 kg) rolling passively at 3.9, 4.4 and 5.0 m/s for 5 min on a flat treadmill (0%) whilst connected to a strain gauge force transducer (S-type 9363, Revere Transducers Europe, Breda, The Netherlands). The measured μ was independent of speed and body mass, and the mean μ-value (0.0237) was incorporated in the work rate calculations. Gas exchange values were measured by open-circuit indirect calorimetry using an Oxycon Pro apparatus (Jaeger GmbH, Hoechberg, Germany). Before each measurement, the VO2 and VCO2 gas analyzers were calibrated using high-precision gases (16.00 ± 0.04% O2 and 5.00 ± 0.1% CO2, Riessner-Gase GmbH & co, Lichtenfels, Germany), the inspiratory flow meter was calibrated with a 3 L volume syringe (Hans Rudolph Inc., Kansas City, MO). Heart rate (HR) was measured with a heart rate monitor (Polar S610, Polar Electro OY, Kempele, Finland), using a 5-s interval for data storage. Blood lactate concentration (BLa) was measured on 5 μL samples taken from the fingertip by a Lactate Pro LT-1710 t (ArkRay Inc, Kyoto, Japan).

In a first test, the skier performed 5-min constant work rates at 0.05 and 0.12 inclines in radians when treadmill roller skiing in the skating G3 technique. Gas exchange values were determined by the average of the last minute during each stage. The lactate threshold was defined at the metabolic power when blood lactate began to accumulate (OBLA) (defined as a concentration of 4 mmol/L, as calculated by a linearly interpolated point out of the three measurement points of blood lactate concentration at the incline of 0.05. Maximal metabolic power was tested at an incline of 0.05 in the G3 technique with a starting speed of 4.4 m/s. The speed was increased by 0.3 m/s every minute until exhaustion. VO2 was measured continuously, and the average of the three highest 10-s consecutive measurements determined VO2max and used to calculate the maximal metabolic power. The test was considered to be a maximal effort if the following three criteria were met: 1) a plateau in VO2 was obtained with increasing exercise intensity, 2) respiratory exchange ratio above 1.10, and 3) blood lactate concentration exceeding 8 mmol/L.

Thereafter, the velocity and the angle of incline on the treadmill were varied according to Figure 2 with gas exchange measured continuously. Immediately after finishing the protocol in Figure 2, the skier had a recovery period while skiing at 0.05 incline at 2.2 m/s, inducing a work rate of 125 J/s and an aerobic power of approximately 900 J/s.
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig2_HTML.jpg
Figure 2

Velocity and incline as a function of time while treadmill roller skiing using the skating G3 technique.___: The skier's velocity (v) in 102 m/s as a function of time (t) in seconds....:The treadmill's incline (α) in radians as a function of time (t) in seconds.

Results

The simulated data for steady state concentration of blood lactate during 5-min constant work rates show good agreement with the experimental results, with less than 0.5 mmol/L disagreement (Figures 3 and 4). For the protocol showed in Figure 2, the simulated and experimental data for work rate and metabolic power as functions of time show general good agreement (Figure 5). However, we see that Q a (t) is delayed compared to the quasi steady state value Q ̄ a ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq18_HTML.gif. Additionally, Q a (t) is delayed compared to the experimental data because the measuring apparatus of VO2 has an inherent time lag (around 15 s) that is not modeled. Based on the same test (i.e., Figure 2), anaerobic power is simulated in Figure 6. Thereafter simulated and experimental lactate levels as functions of time were compared in Figure 7. There was a good agreement between the simulated and experimental results at these varying sub-maximal intensities, showing less than 0.5 mmol/L disagreements. During the recovery period after performing the protocol in Figure 2, the initial simulation model of the concentrations of blood and muscle lactate did not fit well with experimental results. The experimental data showed a much slower reduction in blood lactate concentration compared to the simulated model. Thus, in a second trial we assumed a maximum rate of flux of lactate from blood to muscles. We achieved this mathematically when C m (t) - C b (t) → Max[-0.1, C m (t) - C b (t)] which accounts for a restriction on the lactate flux from blood to muscles. This means that the flux of lactate from the blood to the muscles has a least upper bound. The value of -0.1 was the value that fitted the experimental data best by visual inspection. Using this assumption, the simulated and experimental data were in much better agreement (see Figure 7).
https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig3_HTML.jpg
Figure 3

Steady state concentrations of blood lactate (C) in mmol/L as a function of the fraction of maximum aerobic power while treadmill roller skiing using the skating G3 technique at a 0.05 incline. ■: Experimental data ------: C = C ̄ = 1 2 Ln 1 + χ p 0 d 0 Q max ( 1 - Q a / Q max ) 1 - χ p 0 d 0 Q max ( 1 - Q ̄ a / Q max ) / χ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq19_HTML.gif.

https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig4_HTML.jpg
Figure 4

The concentration of blood and muscle lactate concentration (C) in mmol/L for a skier while treadmill roller skiing in the skating G3 technique at a 0.05 incline at velocities of 3.89 m/s, 5 m/s and 7.44 m/s, respectively. * = experimental data, ___: Muscle compartment simulation, --------: Blood compartment simulation.

https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig5_HTML.jpg
Figure 5

Calculations of work rate (W) and metabolic powers (Q) as functions of time (t) in J/s for a skier while treadmill roller skiing in the skating G3 technique at the velocities and inclines shown in Figure 2. _____: Upper Q ̄ a ( t ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq20_HTML.gif,____: Lower: Work rate (W), - - - - -: Experimental data, .........:Qa(t).

https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig6_HTML.jpg
Figure 6

Simulated anaerobic power (Q) in J/s as a function of time (t) in seconds for a skier while treadmill roller skiing in the skating G3 technique at the velocities and inclines shown in Figure 2 and a subsequent recovery period at 0.05 incline at 2.2 m/s.

https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_Fig7_HTML.jpg
Figure 7

Simulated blood and muscle lactate concentrations (C) in mmol/L as a function of time (t) in seconds for a skier while treadmill roller skiing in the skating G3 technique at the velocities and inclines shown in Figure 2 and a subsequent recovery period at 0.05 incline at 2.2 m/s. *: Experimental data, ___: Muscle compartment simulation, -------: Blood compartment simulation.

Conclusions

The current study constructed mathematical models of lactate production and removal and compared these with experimental results from treadmill roller skiing. The main findings were as follows: 1) a mathematical two-compartment model of lactate production and removal could accurately predict blood lactate concentration during steady state and at varying exercise intensities and; 2) the understanding of lactate removal after high-intensity exercise during whole-body exercise requires further examination.

In the current study, the simulated lactate production and removal fitted the experimental data well during steady state and at varying sub-maximal exercise intensities, and indicate an uncertainty of less than ± 0.5 mmol/L using these mathematical models. To the best of our knowledge, the current study is the first to compare a simulation model of blood lactate concentration during dynamic whole-body exercise against experimental data. Overall, we propose that the current simulation models provide useful insight into how blood lactate concentration during practical training and testing in dynamical situations should be interpreted.

The initial simulation model of lactate removal after high exercise intensities deviated from the experimental findings. This can be related to a longer time delay between lactate concentrations in muscle and blood than simulated by the model [22] or a slower removal of blood lactate than the bi-exponential time function employed here [2628]. In any case, when assuming a maximum rate of flux of lactate from blood to muscles did the simulation model fit experimental data satisfactorily. A rationale for this asymmetry of the time scale for influx and outflux of lactate from the blood may be related to the transport of lactate across membranes accomplished by monocarboxylate transport proteins (MTC). MCT are used for active transport of lactate and is driven by concentration gradients [9, 13]. However, the physiological role of the two important isoforms MCT1 and MCT4 differ. MCT1 is mainly located in slow twitch muscle fibers and is regarded responsible for the influx of lactate to the high oxidative muscles, whereas MCT4 is mainly located in fast muscle fibers and is responsible for outflux of lactate from cells. The current data indicate that the amount of MTC4 is low in blood cells, which leads to a smaller outflux of lactate from the blood during recovery than the influx of lactate to the blood during high intensity exercise. Nevertheless, the understanding of lactate removal after high-intensity exercise during whole-body exercise requires further examination.

Modeling in human biology is always a challenge since one is confronted with conceiving a simple, yet realistic representation of complex phenomena occurring at different levels (cells, organs, tissue). The parameters used in the current study will, in principle, be dependent on the exercise mode employed; fitness level of the individual tested and is only applicable for fixed concentrations of glycogen e.g., [4244]. Thus, to construct valid simulation models of lactate concentration during dynamic exercise these need to be developed for each individual.

Funding

No sources of funding.

Abbreviations

W: 

Work rate in J/s

v: 

Velocity of the treadmill in m/s

α: 

Angle of inclination in radians

m = 77.5 kg: 

Mass of the skier

μ = 0.024: 

Coefficient of roller friction on the treadmill

g = 9.82 m/s2: 

Acceleration of gravity

Q max = 1886: 

Maximal aerobic power in J/s

Q ̄ L T = 0 . 93 Qmax https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq21_HTML.gif

, Aerobic power at the lactate threshold relative to the maximal aerobic power

Qa: 

Aerobic power in J/s

Q ̄ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq3_HTML.gif, Steady state aerobic power in J/s: 

.

Q ̃ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq7_HTML.gif, Virtual aerobic power in J/s: 

.

Q a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq4_HTML.gif, Virtual steady state aerobic power in J/s: 

.

Qb = 80 J/s: 

Resting aerobic power

Q a (t0) = Q b

Initial aerobic power

Qul = 111 J/s: 

Aerobic power of unloaded movement of arms and legs

f

Cycle frequency of the skier in the G3 technique in 1/s

Ā ( α ) = 0 . 92 ( 1 + 1 . 19 E x p ( - 71 . 8 α ) ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq22_HTML.gif, Function that describes the influence of incline on Q ̃ ̄ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq23_HTML.gif

.

c2 = 5.8 J/s, Parameter that describes the influence of work rate on Q ̃ ̄ a https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq23_HTML.gif

.

τ= 30 s: 

Time parameter in seconds quantifying the time before the aerobic power reaches steady state during sub-maximal work rates

C

Lactate concentration in the lactate pool in kg/ml or in mmol/L

( C ̄ ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq16_HTML.gif, Steady state lactate concentration in the lactate pool in kg/ml or in mmol/L

.

C(t0) = 0.045kg/m3 = 0.5mmol/L

Initial concentration of lactate

R a

Rate of lactate appearance in kg/(sm3) or in mmol/(sL)

R d

Rate of lactate disappearance in kg/(sm3) or in mmol/(sL)

P max

Maximum rate of pyruvate appearance in kg/(sm3) or in mmol/(sL)

P

Rate of pyruvate appearance in kg/(sm3) or in mMol/(sL)

P ̄ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq24_HTML.gif, Steady state rate of pyruvate appearance in kg/(sm3) or in mmol/(sL): 

.

P ̃ https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq25_HTML.gif, Virtual rate of pyruvate appearance in kg/(sm3) or in mmol/(sL): 

.

( P ̃ ̄ ) https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq9_HTML.gif, Virtual steady state pyruvate appearance in kg/(sm3) or in mmol/(sL): 

.

τ an = 10 s: 

Time parameter quantifying the time before the pyruvate appearance reaches steady state during sub-maximal work rates

α0 = 0.9: 

Parameter that describe the rate of pyruvate disappearance due to oxidation in the mitochondria

β0 = 1/(0.6 ×p0 ×Q max ): 

Parameter that describes the rate of pyruvate disappearance due to oxidation in the mitochondria

p0 = 10-5kg/(m3s)/(J/s): 

Parameter that describes the rate of lactate appearance as a function of aerobic power

D(Q a (t)): 

Function that describes the rate of lactate disappearance as a function of aerobic power

d0 = 7.2410-8/(J/s)2/s

Parameter that describes the rate of lactate disappearance as a function of aerobic power

χ= 0.95m3/kg

Parameter that describes the saturation of disappearance of lactate

Q a n G https://static-content.springer.com/image/art%3A10.1186%2F1742-4682-9-7/MediaObjects/12976_2011_Article_328_IEq26_HTML.gif, Power caused by the use of ATP produced anaerobically from glycolysis or glycogenolysis: 

.

λ= 3 × 20J

/(kgmmol/L), Parameter that scales the rate of change in lactate concentration and anaerobic power

V m = 10L

Lactate volume of the muscles

V b = 4L

Lactate volume of the blood

V = V b + V m

Total lactate volume

K1a= 0.05/s ×L

Parameter that scales the movement of lactate into and out of the blood

Declarations

Authors’ Affiliations

(1)
Department for Protection Norwegian Defence Research Establishment,
(2)
Department of Human Movement Science,, Norwegian University of science and Technology,

References

  1. Ganong WF: Review of Medical Physiology. 2001, New York: McGraw-Hill Companies, 20
  2. Schmidt-Nielsen K: Animal Physiology: adaptation and environment. 1997, Cambridge: Cambridge University Press, 5
  3. Johnson AT: Biomechanical and Exercise Physiology. 1991, New York: Wiley- Interscience Publication
  4. Margaria R, Cerretelli P, Mangili F: Balance and kinetics of anaerobic energy release during strenuous exercise in man. J Appl Physiol. 1963, 19: 623-628.
  5. Heck H, Mader A, Hess G, Mucke S, Muller R, Hollmann W: Justification of the 4 - mmol/l lactate threshold. International Journal of Sports Medicine. 1985, 6: 117-130. 10.1055/s-2008-1025824.View ArticlePubMed
  6. Billat VL, Sirvent P, Py G, Koralsztein JP, Mercier J: The concept of maximal lactate steady state. Sports Medicine. 2003, 33: 407-426. 10.2165/00007256-200333060-00003.View ArticlePubMed
  7. di Prampero PE, Ferretti G: The energetics of anaerobic muscle metabolism: a reappraisal of older and recent concepts. Respir Physiol. 1999, 118: 103-115. 10.1016/S0034-5687(99)00083-3.View ArticlePubMed
  8. Brooks GA: Lactate production under fully aerobic conditions: the lactate shuttle during rest and exercise. Fed Proc. 1985, 45: 2924-2929.
  9. Brooks GA: Current concepts in lactate exchange. Medicine and Science in Sports and Exercise. 1991, 23: 895-906.View ArticlePubMed
  10. Brooks GA: Lactate shuttles in nature. Biochem Soc Trans. 2002, 30: 258-264.View ArticlePubMed
  11. Brooks GA: Link between glycolytic and oxidative metabolism. Sports Medicine. 2007, 37: 341-343. 10.2165/00007256-200737040-00017.View ArticlePubMed
  12. Myers J: Dangerous curves: a perspective on exercise, lactate, and the anaerobic threshold. Chest. 1997, 111: 787-795. 10.1378/chest.111.3.787.View ArticlePubMed
  13. Juel C: Lactate-proton cotransport in skeletal muscle. Physiol Rev. 1997, 77: 321-358.PubMed
  14. Hermansen L, Osnes JB: Blood and muscle pH after maximal exercise in man. J Appl Physiol. 1972, 32: 304-308.PubMed
  15. Hogan MC, Welch DC: Effect of varied lactate levels on bicycle performance. J Appl Physiol. 1984, 57: 507-513.PubMed
  16. Westerblad H, Lee JA, Allen DG: Cellular mechanisms of fatigue in skeletal muscle. Am J Physiol. 1991, 261: 195-209.
  17. Sahlin K: Metabolic factors in fatigue. Sports Medicine. 1992, 13: 99-107. 10.2165/00007256-199213020-00005.View ArticlePubMed
  18. Hogan MC, Gladden LB, Kurdak S, Polle DC: Increases (lactate) in working dog muscle reduces tension development independent of pH. Medicine and Science in Sports and Exercise. 1995, 27: 371-377.View ArticlePubMed
  19. Andrews MAW, Godt RE, Nosek TM: Influence of physiological L + -lactate concentrations on contractibility of skinned striated muscle fibers of rabbit. J Appl Physiol. 1996, 80: 2060-2065. 10.1063/1.363099.View ArticlePubMed
  20. Favero TG, Zable AC, Colter D, Abramson JJA: Lactate inhibits Ca2 + -activated Ca2 + - channel activity from skeletal muscle sarcoplasmic reticulum. J Appl Physiol. 1997, 82: 447-452.PubMed
  21. Hermansen L, Saltin B: Oxygen uptake during maximal treadmill and bicycle exercise. J Appl Physiol. 1969, 26: 31-37.PubMed
  22. Robergs RA, Ghiasvand F, Parker D: Biochemistry of exercise-induced metabolic acidosis. Am J Physiol Regul Integr Comp Physiol. 2004, 287: 502-516. 10.1152/ajpregu.00114.2004.View Article
  23. Freund H, Zouloumian P: Lactate after exercise in man I. Evolution kinetics in arterial blood. European Journal of Applied Physiology. 1981, 46: 121-133. 10.1007/BF00428865.View Article
  24. Zouloumian P, Freund H: Lactate after exercise in man II. Mathematical model. European Journal of Applied Physiology. 1981, 46: 135-147. 10.1007/BF00428866.View Article
  25. Zouloumian P, Freund H: Lactate after exercise in man III. Properties of the compartment model. European Journal of Applied Physiology. 1981, 46: 149-160. 10.1007/BF00428867.View Article
  26. Oyono-Enguelle S, Freund H, Lampert E, Lonsdorfer A, Lonsdorfer J: Modeling lactate kinetics during recovery from muscular exercise in human. Influence of some physiological factors. Science & Sports. 1993, 8 (31): 181-187.View Article
  27. 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.PubMed
  28. 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.View Article
  29. Rossiter HB, Ward SA, Doyle VL, Howe FA, Griffiths JR, Whipp BJ: Interference from pulmonary O2 uptake with respect to intramuscular (phosphcreatine) kinetcs during moderate exercise in humans. J Physiol. 1999, 518: 921-932. 10.1111/j.1469-7793.1999.0921p.x.PubMed CentralView ArticlePubMed
  30. Crowther GJ, Kemper WF, Carey MF, Conley KE: Control of glycolysis in contracting skeletal muscle. II. Turning it off. Am J Physiol Endocrinol Metab. 2002, 282: 74-79.
  31. Jubrias SA, Esselman PC, Price LB, Cress ME, Conley KE: Large energetic adaptations of elderly muscle to resistance and endurance training. J Appl Physiol. 2001, 90: 1663-1670. 10.1063/1.1374452.View ArticlePubMed
  32. Lanza IR, Wigmore DM, Befroy DE, Kent-Braun JA: In vivo ATP production during free-flow and ischaemic muscle contractions in humans. J Physiol. 2006, 577: 353-367. 10.1113/jphysiol.2006.114249.PubMed CentralView ArticlePubMed
  33. Norman RW, Komi PV: Mechanical energetics of world-class cross-country skiing. International Journal of Sport Biomechanics. 1987, 3: 353-369.
  34. Norman RW, Ounpuu S, Fraser M, Mitchell R: Mechanical power output and estimated metabolic rates of Nordic skiers during Olympic competition. International Journal of Sport Biomechanics. 1989, 5: 169-184.
  35. Sandbakk Ø, Holmberg HC, Leirdal S, Jakobsen V, Ettema G: Analysis of a sprint ski race and associated laboratory determinants of world-class performance. Eur J Appl Physiol. 2011, 111: 947-957. 10.1007/s00421-010-1719-9.PubMed CentralView ArticlePubMed
  36. di Prampero PE: Factors limiting maximal performance in humans. European Journal of Physiology. 2003, 90: 420-429. 10.1007/s00421-003-0926-z.View Article
  37. Moxnes JF, Hausken K: A Dynamic Model for Training Impulse and Lactate Influx and Outflux During Exercise. International Journal of Modern Physics C. 2009, 20: 147-177. 10.1142/S0129183109013522.View Article
  38. Meyer RA: A linear model of muscle respiration explains monoexponential phosphocreatine changes. J Appl Physiol. 1988, 75: 648-656.
  39. Bergman BC, Wolfel EE, Butterfield GE, Lopaschuk GD, Casazza GA, Horning MA, Brooks GA: Active muscle and whole body lactate kinetics after endurance training in men. J Appl Physiol. 1999, 87: 1684-1696.PubMed
  40. Sandbakk Ø, Holmberg HC, Leirdal S, Ettema G: The physiology of world-class sprint skiers. Scandinavian Journal of Medicine and Science in Sports. 2010
  41. Sandbakk Ø, Holmberg HC, Leirdal S, Ettema G: Metabolic rate and gross efficiency at high work rates in elite and national level sprint skiers. Eur J Appl Physiol. 2010, 109: 473-481. 10.1007/s00421-010-1372-3.View ArticlePubMed
  42. Green HJ, Hughson RL, Orr GW, Ranney DA: Anaerobic threshold, blood lactate, and muscle metabolites in progressive exercise. J Appl Physiol. 1983, 54: 1032-1038.PubMed
  43. Gollnick PD, Armstrong RB, Saubert IV, Sembrowich WL, Shepherd RE, Saltin B: Glycogen depletion patterns in human skeletal muscle fibers during prolonged work. Pflugers Archive - European Journal of Physiology. 1973, 344: 1-12. 10.1007/BF00587437.View Article
  44. Galbo H, Holst JJ, Christensen NJ: The effects of different diets and of insulin on the hormonal response to prolonged exercise. Acta Physiologica Scandandinavica. 1979, 107: 19-32. 10.1111/j.1748-1716.1979.tb06438.x.View Article

Copyright

© Moxnes and Sandbakk; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Advertisement