Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Experimental & Molecular Medicine (2022 )Cite this article
Molecules in living organisms are in a constant state of turnover at varying rates, i.e., synthesis, breakdown, oxidation, and/or conversion to different compounds. Despite the dynamic nature of biomolecules, metabolic research has focused heavily on static, snapshot information such as the abundances of mRNA, protein, and metabolites and/or (in)activation of molecular signaling, often leading to erroneous conclusions regarding metabolic status. Over the past century, stable, non-radioactive isotope tracers have been widely used to provide critical information on the dynamics of specific biomolecules (metabolites and polymers including lipids, proteins, and DNA), in studies in vitro in cells as well as in vivo in both animals and humans. In this review, we discuss (1) the historical background of the use of stable isotope tracer methodology in metabolic research; (2) the importance of obtaining kinetic information for a better understanding of metabolism; and (3) the basic principles and model structures of stable isotope tracer methodology using 13C-, 15N-, or 2H-labeled tracers.
Biological compounds in living organisms were once thought to be static until the pioneering work by Rudolf Schoenheimer and his associates in the 1930s, who first developed and used a stable isotope tracer technique using 15N-labeled amino acids to investigate the in vivo kinetics of protein turnover1. They demonstrated that the body protein pool is in a dynamic state of constant turnover (i.e., protein synthesis and protein breakdown)1. In his posthumously compiled book in 1946 entitled “The Dynamic State of Body Constituents”, Schoenheimer eloquently described the dynamic state of living systems as follows: “all constituents of living matter are in a steady state of rapid flux”2. Indeed, virtually everything in the body turns over at varying rates to achieve overall “dynamic” homeostasis. A change in the pool size of any compound can occur when there is an imbalance between its rates of appearance and disappearance. For example, hyperglycemia (i.e., elevation of plasma glucose concentration above normal) can take place if the rate of appearance (Ra) of glucose into circulation, typically through elevation of hepatic glucose output (sum of glycogenolysis and gluconeogenesis), exceeds the rate of disposal or disappearance (Rd) of glucose from the circulation into various tissues for purposes such as glycogen synthesis, oxidation, and de novo synthesis of other compounds. Hyperglycemia can occur even when Rd glucose is elevated compared to the normal condition, provided that Ra glucose is accelerated to a greater extent. The relationships between the blood glucose concentration and the Ra and Rd of glucose underscore the importance of quantitative assessment of substrate kinetics in addition to measurement of the static concentration. Similarly, elucidation of the role of the activation state of molecular signaling molecules requires simultaneous quantification of the corresponding substrate kinetics. Snapshot measurements of molecular signaling pathways may not be consistent with the predicted substrate kinetics3. We will first discuss the importance of obtaining substrate kinetic data for a better understanding of the dynamic nature of substrate metabolism. We will then review the basic principles and models of stable isotope tracer methodology that can be used to assess in vivo kinetics at the whole-body and tissue levels.
The dynamic nature of the biological compounds in the body can be envisioned by considering the analogy of a water tank with flows of water in and out (i.e., water turnover). As shown in Fig. 1, it is intuitively clear that the volume of water in the tank is determined by the balance between the rate of appearance (Ra) of “clean” water into the tank and the rate of disposal or disappearance (Rd) of water from the tank. A steady state in the pool size achieved provided that the two rates are equal, regardless of the absolute rates (Fig. 1). However, the “quality” of the water in the tank may vary depending on the rate of water turnover despite the volume of water and other conditions being identical. To increase the volume of water in preparation for “drought”, Ra must increase and/or Rd must decrease. Indeed, there are an unlimited number of combinations of Ra and Rd that lead to the same volume of water (i.e., pool size) while achieving varying degrees of water “quality.” The amount of water in the tank (i.e., the pool size) is analogous to static concentrations and molecular activation states in biological systems and fails to provide critical kinetic information.
While the water pool size does not change if the two rates are identical, regardless of their absolute rates, differences in the water turnover rate may affect the quality of the water. The figure was created with BioRender.com.
The balance between Ra (e.g., protein synthesis) and Rd (e.g., protein breakdown) determines the pool size of a specific individual protein (e.g., creatine kinase)4,5,6 in a subcellular organelle (such as mitochondrion)7,8, tissue/organ (e.g., skin, muscle, and liver)9,10, or even the whole body11,12,13,14. For example, in healthy adults, muscle mass is relatively constant because the loss of muscle proteins as a result of protein breakdown is replaced by continuous muscle protein synthesis. In growing children or resistance-trained adults who experience muscle hypertrophy, muscle mass increases when the rate of protein synthesis exceeds the rate of protein breakdown. On the other hand, in catabolic states such as sarcopenia or cancer cachexia, the opposite (i.e., breakdown > synthesis) occurs over time. Not only is protein pool size determined by the balance between the rates of synthesis and breakdown, identical pool sizes do not necessarily mean identical rates of protein turnover. Different rates of protein turnover can affect protein quality even when the protein pool size remains constant. For example, the quality of skeletal muscle, typically defined as muscle strength normalized to regional muscle mass15, is positively related to the rate of protein turnover16. Thus, it is important to maintain protein turnover (i.e., synthesis and breakdown) at or above a certain rate to maintain the optimal functional quality (the analogy of “clean water” in the water tank) of muscle by ensuring a constant supply of newly synthesized functional proteins (i.e., protein synthesis) to replace old, nonfunctional or damaged proteins that must be removed (i.e., protein breakdown). Dynamic homeostasis involves the continuous turnover of all body constituents at varying rates. Dysregulation of dynamic homeostasis will lead to physiological and pathological consequences.
Despite the dynamic nature of metabolism in vivo, metabolic research has historically relied on the analysis of static, snapshot information such as the abundance of metabolites, mRNA, and protein or the (in)activity of molecular signaling cascades17. Sole reliance on “static” information, collectively termed “statomics”, regarding the metabolic status of a living organism can result in misleading conclusions. A significant amount of evidence has documented mismatches between “statomics”, or the measurement of concentrations or molecular activation states, and metabolic dynamics in humans3,18 and animals17,19. For example, it was demonstrated that 48 h of fasting in rats resulted in significant elevation in the key enzyme for the first committed step in hepatic gluconeogenesis (GNG), i.e., phosphoenolpyruvate carboxykinase (PEPCK)19. Based solely on the expression of PEPCK, the conclusion would be that the 48-h fast increases GNG, whereas the corresponding in vivo flux rate (i.e., GNG) was actually reduced compared to control and 8-h fasting conditions19. This example of a mismatch between statomics and dynamic flux measurements is not limited to glucose metabolism3,20. Discrepancies between static information (e.g., enzyme abundance/molecular activation state) and actual metabolic flux rates for a given metabolic pathway occur because the actual fluxes of substrate metabolism are accomplished by complex interactions among various factors, including substrate availability, enzyme activity, and implicated signaling cascades21. In the following section, we will cover the fundamental principles of stable isotope tracer methodology and its use to assess in vivo metabolic flux rates, with a focus on the basic model structures of tracer methodology that apply to a wide range of applications.
Information on in vivo substrate dynamics can be obtained with various tracer methods. Here, we will provide basic principles of tracer methodology with the emphasis on stable isotope tracers due not only to health issues related to the use of radioactive isotope tracers but also to the versatility of stable isotope tracers in assessing various aspects of metabolic status. Stable isotopically labeled tracers are any molecules with one or more heavier stable isotopes (e.g., 13C, 2H, or 15N isotopes) incorporated somewhere in the molecule22,23,24. Stable isotope tracers may be administered in the chemical form of the tracer (e.g., 13C glucose) or as heavy water (deuterium oxide, 2H2O) that will produce the desired metabolic tracer in vivo.
The calculation of substrate kinetics is predicated on two basic tracer models: (1) tracer dilution and (2) tracer incorporation22,23. These two basic models can be further divided into four tracer models depending upon the number of pools (single pool vs. multiple pools) and the number of precursors (single precursor vs. multiple precursors) in steady or non-steady states. These models may be used either independently or in combination to assess in vivo kinetics. We will discuss model structures in general terms because tracer methodology can be applied to a variety of different topics, and we will illustrate concepts with common examples. For simplicity, we will discuss these models in the postabsorptive state and physiological steady-state conditions. More comprehensive treatment of the topic can be found in the other references22,23,25.
Substrates turn over in living organisms in a single pool (e.g., fatty acids and water)23,26,27,28 or, more commonly, in multiple pools (e.g., glucose and amino acids)23,29,30,31. The tracer dilution model is based on the dilution of tracer administered into the system by the appearance in the same pool of unlabeled tracees. When an isotopic steady state is achieved, meaning that the rates of tracer and tracee appearance are constant over time and that there is a steady-state enrichment of the tracer in the body pools of the tracee, substrate kinetics can be calculated by the same method, regardless of the number of metabolic pools.
We will discuss the example of palmitate, with the understanding that the principles of calculation apply to determining the kinetics of any other FFA and indeed any substrate distributed structurally or functionally in a single pool, such as water or plasma. In the case of palmitate, although there is exchange between the plasma and interstitial fluid, palmitate is considered to reside essentially in the plasma, as the exchange process between the two compartments (i.e., plasma and interstitial fluid) is too slow to be reflected in the kinetics of plasma palmitate32.
The basic model structure for single-pool kinetics is depicted in Fig. 2. In a physiological steady state, the tracee pool size (i.e., concentration × volume of distribution) is constant over time, meaning that Ra tracee is equal to Rd tracee. Ra tracee and thus Rd tracee can be determined with either a bolus or continuous infusion of tracer. In the case of tracer infusion, the rate of infusion into plasma (F) is constant. In the initial time after the start of tracer infusion, tracer enrichment (i.e., tracer to tracee ratio, TTR) is relatively small but continuously increases over time in a single-exponential manner until reaching an isotopic equilibrium or plateau (Ep) where both the F and the “Rd” of the tracer are constant and equal. For a given F, TTR is inversely related to the rate at which the naturally occurring unlabeled tracee appears in the plasma and dilutes the infused tracer.
a For a given tracer infusion rate of 2 µmol/min (F), the Ra of the tracee, a naturally occurring compound of interest to be traced, is equal to the Rd of the tracee in steady-state conditions where the pool size is constant and is determined as F divided by isotopic enrichment at plateau. b With time, tracer enrichment gradually increases and reaches isotopic equilibrium or plateau enrichment (Ep), where the ratio of Ra tracee to F is equal to the ratio of Rd tracee to the rate at which tracer leaves the pool. After rearranging the relations, Ra tracee can be calculated as F divided by Ep during a continuous tracer infusion study. There is a linear correlation between the magnitude of tracer dilution (reduction in TTR) and Ra tracee for a given F.
In the isotopic steady state, when a plateau in isotopic enrichment has been reached in the sampled pool, the ratio of Ra tracee to F will be directly reflected by the TTR. In Fig. 2, the ratio of the number of tracers (red circle, n = 2) and the number of tracees (black circle, n = 8), i.e., TTR, is 0.25 (2/8), which is dictated by the ratio of Ra tracer (i.e., F, 2 µmol/min) to Ra tracee (8 µmol/min). In the experimental setting, Ra tracee is unknown of interest, F is known, and TTR is measured. Because the ratio of F to Ra tracee dictates the ratio of tracer to tracee (i.e., TTR) at isotopic equilibrium, the following relation will be established: F/Ra tracee = Tracer/Tracee. By rearranging the relation, Ra tracee can be calculated as Ra tracee = F/TTR. Ra tracee in a single pool can also be determined with a bolus tracer injection. Isotopic enrichment in a single pool follows a single-exponential decay after a bolus injection of tracer, which enables the calculation of enrichment at time t over time using a single exponential equation (E = E0•e−kt); thus, Ra will be calculated as the amount of the bolus of tracer divided by the area under the curve of the decay curve of the enrichment.
The single-pool model of tracer dilution is also applicable to other assessments, such as the direct determination of functional muscle mass, called the D3-creatine dilution method (Fig. 3). Briefly, the dilution of a small dose of oral D3-creatine is a direct function of the size of the endogenous creatine pool, which is primarily distributed in muscle. The ingested D3-creatine is distributed throughout the skeletal muscle pool of creatine, and subsequently both labeled and unlabeled creatine are irreversibly converted to creatinine in proportion to the relative abundances of labeled and unlabeled creatine. Measuring the creatinine pool enrichment as a reflection of the muscle creatine enrichment is advantageous because creatinine is excreted into the urine, meaning that the method is minimally invasive33,34. Therefore, the total creatine pool size can be calculated by dividing the dose of D3-creatine by the urinary enrichment of D3-creatinine at equilibrium35. Since the creatine pool size in skeletal muscle is relatively constant36, muscle mass can be calculated by dividing the creatine pool size by the muscle creatine concentration.
The method is predicated on the principle of tracer dilution: the magnitude of dilution of a given small dose of oral D3-creatine intake relative to the amount of preexisting creatine located predominantly in skeletal muscle, determined from the ratio of labeled creatinine to unlabeled creatinine in a urine sample, reflects total muscle mass. Labeled and unlabeled creatine are irreversibly converted to labeled and unlabeled creatinine in proportion to their relative concentrations and then excreted in urine. For example, the urine creatinine enrichment (i.e., ratio of D3-creatinine to creatinine) resulting from 20 kg of muscle will be half that resulting from 10 kg of muscle.
Multiple-pool models apply to many substrates, including glucose and amino acids. We will discuss the multiple-pool model of a single substrate using glucose as the example.
Glucose resides in multiple pools (i.e., blood and interstitial fluid) and is one of the most commonly traced substrates. To simplify the conceptualization of the theory of the model structure, we will discuss the situation of a physiological steady state, meaning that Ra and Rd are equal and the blood level of glucose is constant over time. The basic principles of glucose kinetics are depicted in Fig. 4. Glucose enters the blood (pool 1) directly from the liver but disappears irreversibly into many tissues from the interstitial fluid (pool 2). There is constant exchange of both tracer and tracee glucose between pool 1 and pool 2. Although labeled and unlabeled glucose leave pool 2 to enter tissue, it is typically possible to sample only pool 1 (blood). Regardless, it is possible to use the same steady-state equations (i.e., Ra tracee = F/Ep) in this model as for single pool kinetics, as isotopic exchange between pool 1 and pool 2 will eventually result in equal enrichment in the two pools (Fig. 4). In the case of glucose, isotopic equilibrium would be achieved much more quickly if it was distributed only in a single pool, but mixing in interstitial fluid is much slower than mixing in blood. Consequently, the observed increase in enrichment in blood (pool 1) occurs at a rate that reflects the time course of enrichment of both pool 1 and pool 2. The important point of Fig. 4 with regard to the calculation of steady-state kinetics is that since all the glucose that appears in pool 2 is derived from pool 1, the glucose isotopic enrichment in pool 2 will eventually reach the same equilibrium value as that in pool 1 (Fig. 4). Therefore, Ra tracee is calculated as F divided by Ep, where Ep is the plateau enrichment measured in pool 1.
Substrates such as glucose reside in multiple pools. a For example, glucose dwells in plasma (pool 1) and interstitial fluid (pool 2), between which gradual glucose exchange occurs. In the postabsorptive state, glucose appears in pool 1 via endogenous glucose production and diffuses into pool 2, from which it is taken up by tissues. To determine Ra glucose (and thus Rd glucose) in steady state conditions, a tracer glucose is (typically) continuously infused into pool 1 at a specific rate (F). b The observed enrichment at the point of sampling (pool 1) reflects a mixture of the two pools. Regardless of the differences in the speed of mixing in the two compartments, an isotopic equilibrium will be reached where the isotopic enrichment levels of the two pools are identical. Under this condition, Ra glucose can be calculated as Ra tracee = F/Ep. Ra, rate of appearance; Rd, rate of disappearance, Ep, plateau enrichment.
In addition to the Ra of a single substrate distributed in multiple pools being determined with a constant tracer infusion, a bolus of tracer can be used. In the case of multiple-pool models, the measured decay curve in blood following tracer injection will be described by a multiexponential function. The number of exponential curves that best fit the decay data in a multiple model is related to, but not necessarily equal to, the number of physical pools. For example, the decay in glucose enrichment after a bolus tracer injection may be described by either a two- or three-exponential curve, even though the system behaves physiologically as a two-pool system37. Despite a multiexponential decay curve following a bolus tracer injection in a multiple pool model of a single substrate, Ra can still be accurately calculated from the area under the decay curve in a manner analogous to calculating Ra from the area under the decay curve in a single-pool system. Furthermore, by using the compartmental modeling approach22,25, additional kinetic parameters, including individual pool sizes, decay constant rates (k), half-life, and mean residence time, can be calculated.
The calculation of Ra of a single substrate distributed in multiple pools becomes much more complicated in a physiologically non-steady state. The calculation of Ra glucose in the non-steady state using a constant tracer infusion has been the most widely examined approach. The fundamental problem arises from different rates of mixing of the tracer and tracee in the different pools. Thus, with a rapid change in Ra in conditions such as the start of exercise38,39,40,41,42 or food intake containing carbohydrates43,44, the changing Ra glucose mixes rapidly with the tracer in pool 1 (i.e., blood), from which blood samples are taken and enrichment is measured for calculations of glucose kinetics. As a result, the blood enrichment in the early response of non-steady state conditions largely reflects pool 1 rather than the entire distribution volume of glucose. The mixing between pool 1 and pool 2 occurs more slowly than mixing throughout pool 1. The more extensive the mixing between pools 1 and 2 is, the more closely the sampled compartment will reflect the entire glucose pool. The most commonly used solution to the problem of the sampled compartment reflecting pool 1 and an unknown portion of pool 2 goes back to the 1950s, when Steele proposed calculating Ra glucose in a non-steady state by assuming that the observed values of glucose concentration and TTR reflected pool 1 plus a fixed fractional contribution from pool 245. The limitation of this approach is that the extent of mixing between pools 1 and 2 is time-dependent and not well represented by a single fractional value, particularly when the change in the magnitude and direction of change in Ra are complex. Regardless of the many published criticisms of the Steele equation for calculating non-steady-state kinetics, it remains the most popular approach in the absence of a better method. Although quantification of changes in Ra in the non-steady state can be problematic, it should be recognized that a basic characteristic of the model dictates that a decrease in enrichment during a constant tracer infusion can occur only because of an increase in Ra, and vice versa. The interpretation of changes in isotopic enrichment in a non-steady state is thus always qualitatively valid.
The traditional A-V difference method involves measurement of the difference between the arterial inflow of a substrate (i.e., arterial concentration × blood flow) to an individual organ or tissue (such as skeletal muscle) and venous outflow (i.e., concentration of venous blood draining the tissue of interest × blood flow). This method has been pivotal in determining the interorgan net flux of compounds46 but provides no information about the underlying mechanisms responsible for various responses of net substrate uptake or release. The inclusion of isotopically labeled tracers can greatly expand the metabolic insight provided by the A-V model. We will use the example of the measurement of muscle tissue synthesis and breakdown using a phenylalanine (Phe) tracer. We will describe a simple one-pool model, with the understanding that the addition of more complex model designs is possible47.
The A-V tracer model relies on the same principles of tracer dilution models described above. The significant difference is the necessity of accounting for the blood flow into and out of the tissue. A schematic of the model is shown in Fig. 5. Since Phe is an essential amino acid that cannot be produced in the body, the appearance of Phe in the intracellular pool of free Phe (pool 1 in Fig. 5) can come only from the blood (Ra1) and from the breakdown of protein (pool 1 in Fig. 5), reflected as Ra2. Ra2 is the difference between the total Ra (RaT) and Ra1. RaT is determined by classic tracer dilution, whereby Ra tracer into pool 1 (blood flow × arterial Phe concentration × arterial Phe enrichment) is divided by the measured tracer/tracee ratio in pool 1. It is optimal to measure the intracellular enrichment directly, but if it is not possible to obtain a muscle sample, it can be assumed that the tracer/tracee ratio in the venous outflow reflects the enrichment of free Phe in the muscle pool. Since Ra1 is measured (blood flow × arterial concentration of Phe), Ra2 can be calculated. Ra2 is a direct reflection of protein breakdown since Phe cannot be produced by any metabolic pathway. Once Ra2 is known, the rate of protein synthesis, denoted as Rd2, can be calculated as the sum of Ra1 and the difference between Ra2 and Rd1. In the case of a net synthesis of muscle protein, Ra1 will be greater than Rd1, and the reverse will be true when Rd1 is greater than Ra1. Ra1 minus Rd2 is the A-V difference of Phe.
The A-V model utilizes tracer dilution principles adapted for the measurement of tracer flux across individual tissues and organs. The example shown is for the measurement of muscle protein synthesis and breakdown using Phe as a tracer. For the purpose of explanation, the intracellular free Phe is represented as pool 1 and the protein as pool 2. More complex models, such as distinguishing the interstitial fluid pool from the intracellular free pool and distinguishing leg skin metabolism and leg muscle metabolism, are based on the same basic tracer dilution principles. Since Phe is an essential amino acid that cannot be produced in the body, the appearance of Phe in the intracellular pool of free Phe can come only from the blood (Ra1) and from protein breakdown (Ra2). Ra2 is the difference between the total Ra (RaT) and Ra1. RaT is determined by the classic tracer dilution method, whereby the tracer Ra in pool 1 (blood flow × arterial Phe concentration × arterial Phe enrichment) is divided by the measured tracer/tracee ratio of Phe in pool 1. Since Ra1 is measured (blood flow × arterial concentration of Phe), Ra2 can be calculated. Once Ra2 is known, the rate of protein synthesis (reflected as Rd2) can be calculated as the sum of Ra1 and the difference between Ra2 and Rd1. Ra1 minus Rd2 is the A-V difference of Phe.
The tracer incorporation model calculates the rate of synthesis of a product from the rate of tracer incorporation. There is a distinction between the calculation of the rate of production when the product has only one precursor and when there are multiple potential precursors. Determining the fractional synthetic rate (FSR) of protein from the rate of incorporation of an amino acid tracer precursor is a commonly used example of a single-precursor model. FSR is considered to be a single-precursor model because there are no other precursors that might dilute the increase in enrichment over time resulting from the incorporation of the tracer (amino acid) precursor. In contrast, calculating the rate of glucose production from a single tracer precursor, such as lactate, is an example of a multiple-precursor model. In a multiple-precursor model, the rate of increase in tracer enrichment in the product is a function not only of the rate of incorporation of the tracee into the product but also of the rate of incorporation of other (unlabeled) precursors. In the example of glucose, the rate of increase in the enrichment of glucose from lactate will be diluted to a variable extent by the simultaneous production of glucose from other precursors that are not labeled, such as glycerol.
The calculation of muscle protein FSR will be used as an example. The FSR represents the fraction or percentage of the total protein pool size that has been newly synthesized during a given period of time (%/time). The FSR of any form of polymer with repeating (mono or multiple) precursors, such as DNA48,49, RNA50, glucose (i.e., gluconeogenesis)51, fatty acids (i.e., de novo lipogenesis, DNL)52, and triglycerides52,53, can be determined. The calculation of FSR is predicated upon the precursor-product relation54. In muscle protein synthesis, amino acids are precursors, and muscle protein is the product. An individual amino acid tracer is used to calculate the FSR of muscle protein derived from that precursor. The total muscle protein FSR is then calculated by taking into account the percentage of the protein comprising the amino acid being traced. As FSR is a “fractional” term, calculation of the total synthetic rate requires the multiplication of FSR by the total pool size of the product. FSRs between groups or conditions can be directly compared if the protein pool size remains essentially constant throughout the experimental time frame or if there are no significant differences in protein pool sizes between groups. For example, FSR is a direct reflection of the response of muscle protein synthesis when a crossover design is used in which the responses to two different nutritional supplements are compared in the same subjects. In contrast, muscle protein FSR must be multiplied by muscle mass to make a reasonable comparison between a group of elderly women with depleted muscle mass and a group of young healthy males with significantly greater muscle protein mass.
The basic principle of determining muscle protein FSR is illustrated schematically in Fig. 6. A stable isotope amino acid tracer is typically administered as a constant infusion or as a bolus injection (F in Fig. 6). Alternatively, labeled amino acid precursors can be produced in vivo from orally administered deuterium oxide, 2H2O. Tracer enrichment is diluted by the endogenous rate of appearance of the unlabeled counterpart of the tracer (i.e., tracee). If an essential amino acid is used as the tracer, the appearance of tracee in the fasted state is due entirely to release of that amino acid by protein breakdown (Ra2 in Fig. 6), since essential amino acids cannot be synthesized in vivo. A nonessential amino acid can also be used as a tracer, but the appearance of tracees from protein breakdown and de novo synthesis must also be considered (Ra1 in Fig. 6). The tracee and tracer will leave pool 1 (blood) for protein synthesis (Rd2) or intracellular pathways other than protein synthesis, including oxidation (Rd1). Tracer and tracee amino acids will leave pool 1 in proportion to their relative abundances in pool 1. There is an underlying assumption with the measurement of protein FSR that the turnover rate of the protein relative to its pool size is so slow relative to the experimental period that there is essentially no tracer recycling (i.e., appearance of tracer coming back from protein breakdown into the precursor pool)47,55. This assumption is valid in slow-turnover proteins such as muscle proteins27 but not in fast-turnover proteins such as apolipoprotein B (Apo B)56. In the case of Apo B, both the infused and recycled labeled amino acid precursor incorporate very rapidly into Apo B protein, and the product enrichment will rise quickly to a plateau, which represents the true precursor enrichment.
a Both the labeled tracer (F) and the unlabeled amino acid tracee (Ra1) enter pool 1 (intracellular precursor pool) in proportion to their relative rates of appearance, from which both tracer and tracee are either oxidized (Rd1) or incorporated into protein (Rd2). When using phenylalanine (Phe) as a tracer, Rd2 Phe directly reflects the rate of protein synthesis, as irreversible loss of Phe by oxidation (more correctly, hydroxylation) does not occur in skeletal muscle. Rd2 is the rate of muscle protein synthesis (MPS), calculated as the product of FSR multiplied by pool size. b FSR is defined as the % change in the product pool size resulting from new synthesis during a specified time interval. For example, for a continuous tracer infusion of essential amino acid, the calculation of FSR is based on the change in product enrichment between two time points (thus, two muscle samples obtained at t1 and t2), typically 2–8 h apart, divided by the precursor enrichment. The dotted line indicates that the rate of increase in product muscle protein enrichment is slow, and thus, during the experimental period, the product enrichment never approaches the precursor enrichment. Typically, the enrichment of both precursor and product are measured in the early phase of product enrichment but after precursor enrichment is plateaued. Ra, rate of appearance; Rd, rate of disappearance; F, infusion rate.
An amino acid tracer (e.g., L-[ring-2H5]phenylalanine) is typically infused intravenously, and muscle protein is sampled (usually by needle biopsy) at the beginning and the end of tracer infusion. Protein FSR is calculated from the change in enrichment in the protein-bound tracer over a measured interval of time, divided by the enrichment of the precursor amino acid enrichment in the intracellular compartment: FSR (%/time) = [(EBP2 – EBP1)/(EPC × t)] × 60 × 100, where EBP1 and EBP2 are the enrichments of bound tracer amino acid (e.g., L-[ring-2H5]phenylalanine) in the first and second biopsies, respectively, and EPC is the precursor tracer enrichment, which is typically measured in the intracellular fluid of the muscle tissue sample. The time elapsed between the first and second muscle biopsies is represented by t. The factor 60 converts the units from minutes to hours, and the factor 100 expresses the result as a percentage rather than a fraction.
The correct, or “true”, precursor enrichment is key to the accurate calculation of FSR. The importance of precursor enrichment can be appreciated by considering two identical physiological circumstances, but in one case, the tracer infusion rate is twice that in the other. Tracer incorporation would be twice as fast with the higher rate of tracer infusion, even though the actual FSRs are identical. In the case of muscle FSR, it is possible to determine a reasonably accurate value for precursor enrichment by measuring TTR in the intracellular fluid of the muscle sample57. The intracellular free amino acid pool represents tracer precursor enrichment for cultured fibroblasts and myocytes57. In other circumstances, it may be more difficult to identify a reasonable surrogate for the true precursor enrichment. For example, intrahepatic acetate is the precursor for the de novo FSR of fatty acids, and that compartment is not available for sampling. Two different potential solutions to this issue were developed: mass isotopomer distribution analysis (MIDA)58 and isotopic spectral analysis (ISA)59. The basic principle is that true precursor enrichment can be back-calculated from the labeling patterns of the product58,59. A detailed description of these approaches can be found in the references26,60,61.
In addition to administered preformed tracer precursors, such as labeled amino acids, enriched precursors can be produced in vivo from “heavy water” (2H2O, deuterium oxide). The heavy water labeling of precursors has several potential advantages62,63 over the substrate-labeled tracer infusion method11,12,13,14. Most notably, heavy water labels a number of different precursors, which enables the simultaneous determination of FSRs of a variety of polymers with repeating (mono or multiple) precursor molecules, such as DNA64, fatty acids (de novo lipogenesis)65, and glucose (i.e., gluconeogenesis)66 in free-living conditions. Water moves rapidly and freely throughout the body, including across cell membranes, and equilibrates with all of the body fluid within ~30 min after administration67,68. However, there are several limitations of the heavy water method24. The time required for the measurement of FSR may be much longer than when a prelabeled precursor tracer is infused intravenously. For example, muscle protein FSR can be measured before and after a single meal when an amino acid tracer is infused intravenously at a constant rate, whereas the time frame of labeling an amino acid precursor from heavy water is too long for that experimental design to be practical11,12,13,14,24,67,69. In addition, there is no method to assess the fractional protein breakdown rate (FBR) using heavy water, although it is possible to infer changes in FBR over time from the measured FSR and changes in the protein pool size70. Finally, the amount of 2H2O that can be administered to a human without an adverse response may produce insufficient precursor labeling to easily detect its rate of incorporation into a slow-turnover compound.
The basic principle of heavy water labeling to measure muscle protein FSR is depicted in Fig. 7a, b, d. Heavy water labeling is an approach to precursor labeling, and the determination of FSR relies on the same principles discussed above. To determine muscle protein FSR, 2H2O is administered (after bolus priming) orally. Heavy water labeling has been used in mice71,72, rats69, and humans24. Consumed 2H2O (dideuterated water, M + 2) rapidly equilibrates with body fluid, which in turn produces monodeuterated water (2H1HO) via exchanges of hydrogen (either 1H or 2H) between 2H2O (M + 2) and 2H1HO (M + 1). Due to the relatively slow turnover rate of body water (half-life: ~10 days)24, the water enrichment established following a dose is maintained for a relatively long time period, during which the 2H of water is incorporated into amino acids via transamination reactions, and the labeled amino acids are precursors for protein synthesis24. There is a “fixed” relation between water enrichment (M + 1) and amino acid enrichment (M + 1) because the number of 2H labeling sites for a specific amino acid is constant26. For example, there are theoretically 4 labeling sites of 2H from water in alanine68, but that value has been empirically found to be ~3.7 for mice, rats, and humans. Enrichment of the product (e.g., protein-bound alanine) is determined from muscle tissue sampled at the end of the experimental period divided by the precursor enrichment, as described above. Precursor amino acid enrichment may be determined directly. Alternatively, amino acid precursor enrichment may be determined indirectly from body water enrichment (reflected by any fluid in equilibrium with body water, such as saliva, plasma, urine, etc.) because of “fixed” relations between water enrichment and alanine enrichment.
a A small dose of heavy water (deuterium oxide, 2H2O) is administered to achieve target enrichment (~1% for humans and ~5% for animals) with which to assess the fractional synthesis rate (FSR) of polymers such as muscle protein. 2H2O administered into the system will equilibrate with body water within a relatively short time (<30 min), forming monodeuterated water (i.e., 2H1HO). The 2H will exchange with H attached to various precursors, including amino acids. The number of exchanges and sites of exchange are known for a variety of compounds, including amino acids. For example, there are 4 sites in the case of alanine. The precursors labeled by exchange reactions with heavy water are incorporated into the products of interest, such as muscle protein. b Muscle protein FSR is typically determined from a muscle sample obtained by needle biopsy. d Specific compartments of a muscle sample can be labeled, such as contractile proteins, mitochondrial proteins, sarcoplasmic proteins, and pooled proteins, as well as the free intracellular precursor (e.g., alanine). c Muscle protein FSR can also be assessed without obtaining muscle tissues using the “virtual biopsy” technique. The virtual biopsy technique involves the measurement of tracer incorporation into specific proteins in blood that primarily reflect muscle protein metabolism. e In addition, FSR can be assessed at the proteome level, reflecting the FSR of several hundred or even thousands of individual proteins. The single-precursor incorporation model applies to the calculation of FSR in all the applications described in Fig. 7.
FSR is calculated using the tracer incorporation model discussed above, wherein the rate of incorporation of labeled precursor over time is divided by the enrichment of the precursor. Because water enrichment dictates the amino acid enrichment (EAA), the precursor enrichment can be expressed as the product of the sampled water enrichment (EH2O) and the number of deuterium labeling sites of the tracee amino acid: FSR(%/t) = [EAA/(EH2O × # labeling site × t)] × 10024. It is assumed that FSR is constant (i.e., the function is linear) over the experimental time period, which may be problematic if sampling occurs over a prolonged interval. In the case of protocols spanning several days, it must be assumed that the various changes that occur in FSR within that time frame are averaged out by sampling at the beginning and end of the experimental period.
The incorporation model is well suited to measuring muscle protein FSR but limited by the invasiveness of muscle biopsies. The “virtual biopsy” method has been developed to circumvent the requirement for biopsies to measure muscle protein FSR (Fig. 7c)5. Instead of measuring product enrichment directly from muscle biopsies, the method measures the enrichment of a circulating protein that is mostly released from muscle, such as creatine kinase muscle isoform (CKM, >90% from skeletal muscle) or carbonic anhydrase 3 isoform (CA-3, virtually specific to skeletal muscle)6. Using this approach, it was demonstrated in humans that there are strong correlations between muscle protein FSR and “virtual” FSR (r > 0.89 for CKM; r > 0.87 for CA-3)5 and between “virtual” FSR and myofibrillar proteins, including actin, myosin-2, and tropomyosin alpha-1 (for all, r > 0.66)5. Therefore, the “virtual biopsy” method may provide practical advantages in assessing muscle dynamics when muscle biopsies are not feasible despite its inability to differentiate among muscle tissues.
Proteomics refers to the large-scale analysis of proteins. The abundance of hundreds or even thousands of individual proteins is determined in the same sample. Proteomics is conceptually linked to the description of the entire set of proteins produced in a biological system, which is called the proteome. Proteomics has become a widely used tool to examine the effects of nutritional or pharmacological treatments73. However, as explained in the example provided above regarding the limited kinetic information provided by measuring only the amount of water in the tank, proteomics is limited by the “static” nature of the data. Exploring dynamic changes in the proteome may represent an important advance in proteomic analysis. Measuring dynamic changes in the proteome for days to months with heavy water labeling in conjunction with proteomics can potentially solve the intrinsic limitation of proteomics5,6. With this kinetic approach, the dynamic changes of hundreds to thousands of proteins can be measured simultaneously. The proteome kinetics technique is categorized into four steps: (1) orally ingested 2H2O labels amino acids with deuterium, which are in turn incorporated into proteins; (2) labeled proteins are digested with trypsin to generate labeled peptides for analysis with high-resolution mass spectrometry, and individual proteins are identified from known specific sequences of peptides through database search; (3) turnover rates of individual proteins are determined based on the abundance and patterns of labeling of the peptides by the tracer incorporation technique; and (4) network analysis of FSR of individual proteins in the proteome can enable the identification of biologically meaningful groups of proteins (Fig. 7e). Note that this process relies on the single-precursor incorporation tracer model described above as applied to a variety of individual products.
When a substrate such as glucose is produced from more than one precursor, the appropriate tracer model accounts for multiple-pool kinetics with multiple precursors for the production of a single product. For example, glucose is produced from multiple precursors, such as lactate, glycerol, pyruvate, and various amino acids42,74,75,76,77,78. Determination of the rate of oxidation from a specific substrate is another example of a multiple-pool, multiple-precursor model. In the latter example, CO2 is produced by oxidative phosphorylation in mitochondria from various substrates, including glucose, FFAs, amino acids, and multiple other compounds79,80,81. The multiple-pool, multiple-precursor model is distinct from the single-precursor model because the kinetics of substrates other than the tracee can influence the calculated rate of synthesis of the product. For example, consider the determination of the rate of glucose production from lactate using labeled lactate as the tracer in two different circumstances in which glucose production from lactate is equal. Calculation of the rate of glucose production from lactate by dividing the tracer enrichment in the product (glucose) by the precursor (lactate enrichment), as described above for a single precursor model, will provide only the fraction of glucose production from lactate in this case. The total rate of glucose production from lactate is not determined with only one tracer. The fraction of glucose production from lactate must be multiplied by the total rate of glucose production to calculate the rate of glucose production from lactate. The rate of total glucose production must be determined independently by tracer dilution of a glucose tracer, such as 6,6 2H2-glucose. The reason that a multiple-precursor model is required to calculate the rate of glucose production from lactate is that variable rates of glucose production from other precursors, such as glycerol, amino acids, and glycogen, will affect the total rate of production of glucose and thus the fraction of glucose production from lactate. In this example, the measured fraction of glucose derived from lactate (i.e., product enrichment divided by precursor enrichment) will be reduced by an elevated total rate of glucose production, even if the absolute rate of glucose production from lactate remains constant. Multiplication of the fraction of glucose derived from lactate by the total rate of glucose production is necessary to obtain meaningful data. Thus, in contrast to the single-precursor model, the labeling of the product in a multiple-precursor model is not necessarily a reflection of the rate of production from the labeled precursor. Rather, it is necessary to consider the total rate of production of the product, which must be determined by the tracer dilution method independently of the incorporation of the tracer of the precursor of interest. In this section, the impact of possible isotopic exchange in the TCA cycle will not be considered, as we have focused on the structure of the tracer model. The principle of the multiple-pool, multiple-precursor model is shown schematically in Fig. 8.
a The 13C-lactate tracer is infused into plasma (pool 1). Lactate tracee enters pool 1 from various tissues, especially muscle (Ra1). Lactate can leave pool 1 by various pathways, including glucose production (Rd2). To determine the rate of glucose production from lactate (i.e., Rd2), the fractional contribution of lactate to total glucose production must be calculated as the ratio of the plasma enrichment of 13C-glucose to that of 13C-lactate, multiplied by the total glucose production rate or Ra glucose, which is determined by a separate glucose tracer (e.g., 2H2-glucose). b 13C-glucose enrichment derived from 13C-lactate will never reach the plateau enrichment of plasma lactate because of the dilution of label incorporation by the incorporation of other unlabeled precursors (e.g., glycerol, amino acids, and glycogen). The difference in enrichment between 13C-glucose and 13C-lactate represents the contribution of other precursors, such as glycogen, various amino acids, and glycerol.
The quantification of intracellular metabolic flux (e.g., glycolysis, the pentose phosphate pathway, the Krebs cycle, etc.) has long been of interest because of the important roles these pathways play in maintaining metabolic homeostasis. However, it is generally not possible to quantify intracellular metabolic flux rates with the traditional application of tracer methodology. The problem is that most components of metabolic pathways can be produced from multiple precursors, in which case the multiple-pool, multiple-precursor model applies. This model requires the measurement of both the rate of incorporation of labeled precursor into the product and the total rate of intracellular production of the product. The latter value is generally not accessible in an in vivo experiment designed to assess multiple intermediates in a metabolic pathway. Metabolic flux analysis (MFA) with the use of stable isotope tracers is an approach to approximate the flux rates of components of metabolic pathways without having all the data necessary to directly measure the rates82,83,84. MFA estimates flux rates through several steps, as represented schematically in Fig. 9. The first step is the selection of the optimal tracer and construction of the stoichiometric matrix model. The selection of optimal tracers is determined via in silico simulations (i.e., computer-based simulation)85, and databases (e.g., the Kyoto Encyclopedia of Genes and Genomes) can be used to construct the stoichiometric matrix model. The next step is to perform tracer labeling experiments and measure isotopic labeling patterns and external rates (i.e., substrate uptake and release rates). Since isotopic labeling pattern data provide only relative information, external rates are required to quantify intracellular metabolic flux rates. The third step is flux estimation, which can be accomplished by several software tools, including INCA and Metran, which use the elementary metabolite unit (EMU)86 framework for the simulation of isotopic labeling and flux. The estimated flux results are modified to minimize the residual sum of squares (RSS) between simulated and observed results by adjusting the metabolic network model87. A lower RSS indicates better agreement between the simulated and observed results. Since the minimized RSS is a random variable with a χ2 distribution, the confidence intervals of the RSS can be calculated to verify the flux result. If the flux result is within the statistically acceptable range (e.g., 0.05 for the 95% confidence intervals), it is considered to represent the absolute rates of intracellular metabolic flux. While MFA can provide quantitative flux rates, it should be recognized that the method has been developed to overcome the fact that the requirements of the multiple-precursor, multiple-product model cannot be satisfied. While the precision of results from MFA can be assessed statistically, empirical validation of its accuracy is usually not possible.
MFA is performed as follows: a Selection of isotopic labeling experiments; b metabolic network model construction; c stable isotope tracer labeling experiments and measurement of isotopic labeling patterns of metabolites derived from the tracer substrate; d estimation of intracellular metabolic fluxes using statistical analysis; and assessment of flux values by using the goodness of fit and confidence intervals. e Results of metabolic flux rates of the pathway. RSS, residual sum of squares.
Living organisms are in a constant state of metabolic turnover. Traditional “omics” and molecular tools provide static, snapshot information, collectively termed “statomics”. Stable isotope tracer methodology, on the other hand, provides dynamic data that can quantify metabolic flux rates. Tracer methods are all predicated on two basic model structures, i.e., tracer dilution and tracer incorporation. The basic model structures can be further divided according to the number of pools (single pool or multiple pool) and/or the number of precursors (i.e., single or multiple precursor). Use of tracer dilution alone or together with tracer incorporation can enable the assessment of in vivo and in vitro dynamics of various biological compounds. The tracer dilution model is also applicable to other usages, including direct assessment of functional muscle mass based on the dilution of consumed D3-creatine by the existing unlabeled creatine pool. A key element of tracer methodology is that the tracer labeling of a product, regardless of the model structure, provides fractional information but does not provide actual rates without knowledge of pool size in the case of a single precursor or total production rate in the case of the multiple-precursor model. The use of heavy water labeling has increased over the past decades due to the minimally invasive nature of the method and the diverse array of potential metabolic precursors and products that can become labeled in vivo. However, when using heavy water labeling, the basic principles of the model structure still apply. On the other hand, traditional tracer methodology generally will not apply to the measurement of multiple flux rates in intracellular metabolic pathways. Metabolic flux analysis (MFA) enables the approximation of flux rates that cannot be directly quantified by traditional tracer methodology.
Schoenheimer, R., Ratker, S. & Rittekberg, D. Studies in protein metabolism VII. The metaboism of tyrosine. J. Biol. Chem. 127, 333–344 (1938).
Schoenheimer, R. The Dynamic State of Body Constituents (Harvard University Express, 1946).
Bukhari, S. S. I. et al. Intake of low-dose leucine-rich essential amino acids stimulates muscle anabolism equivalently to bolus whey protein in older women at rest and after exercise. Am. J. Physiol. Metab. 308, E1056–E1065 (2015).
Hines, K. M. et al. Application of high-resolution mass spectrometry to measure low abundance isotope enrichment in individual muscle proteins. Anal. Bioanal. Chem. 407, 4045–4052 (2015).
CAS PubMed PubMed Central Article Google Scholar
Shankaran, M. et al. Circulating protein synthesis rates reveal skeletal muscle proteome dynamics. J. Clin. Invest. 126, 288–302 (2016).
Murphy, C. H. et al. Effect of resistance training and protein intake pattern on myofibrillar protein synthesis and proteome kinetics in older men in energy restriction. J. Physiol. 596, 2091–2120 (2018).
CAS PubMed PubMed Central Article Google Scholar
Burd, N. A. et al. Low-load high volume resistance exercise stimulates muscle protein synthesis more than high-load low volume resistance exercise in young men. PLoS ONE 5, e12033 (2010).
PubMed PubMed Central Article CAS Google Scholar
Wilkinson, S. B. et al. Differential effects of resistance and endurance exercise in the fed state on signalling molecule phosphorylation and protein synthesis in human muscle. J. Physiol. 586, 3701–3717 (2008).
CAS PubMed PubMed Central Article Google Scholar
Zhang, X. J., Chinkes, D. L., Doyle, D. & Wolfe, R. R. Metabolism of skin and muscle protein is regulated differently in response to nutrition. Am. J. Physiol. 274, E484–E492 (1998).
Sheffield-Moore, M. et al. Mixed muscle and hepatic derived plasma protein metabolism is differentially regulated in older and younger men following resistance exercise. Am. J. Physiol. Metab. 288, E922–E929 (2005).
Kim, I.-Y. et al. Quantity of dietary protein intake, but not pattern of intake, affects net protein balance primarily through differences in protein synthesis in older adults. Am. J. Physiol. Endocrinol. Metab. 308, E21–E28 (2015).
CAS PubMed Article Google Scholar
Kim, I.-Y. et al. The anabolic response to a meal containing different amounts of protein is not limited by the maximal stimulation of protein synthesis in healthy young adults. Am. J. Physiol. Metab. 310, E73–E80 (2016).
Kim, I.-Y. et al. Quality of meal protein determines anabolic response in older adults. Clin. Nutr. 37, 2076–2083 (2018).
CAS PubMed Article Google Scholar
Kim, I.-Y. et al. Protein intake distribution pattern does not affect anabolic response, lean body mass, muscle strength or function over 8 weeks in older adults: A randomized-controlled trial. Clin. Nutr. 37, 488–493 (2018).
CAS PubMed Article Google Scholar
Park, S. W. et al. Decreased muscle strength and quality in older adults with type 2 diabetes: the health, aging, and body composition study. Diabetes 55, 1813–1818 (2006).
CAS PubMed Article Google Scholar
Fitts, R. H. et al. The deleterious effects of bed rest on human skeletal muscle fibers are exacerbated by hypercortisolemia and ameliorated by dietary supplementation. Am. J. Physiol. Physiol. 293, C313–C320 (2007).
Neinast, M. D. et al. Quantitative analysis of the whole-body metabolic fate of branched-chain amino acids. Cell Metab. 29, 417–429.e4 (2019).
CAS PubMed Article Google Scholar
Robinson, M. M. et al. Enhanced protein translation underlies improved metabolic and physical adaptations to different exercise training modes in young and old humans. Cell Metab. 25, 581–592 (2017).
CAS PubMed PubMed Central Article Google Scholar
Perry, R. J. et al. Leptin mediates a glucose-fatty acid cycle to maintain glucose homeostasis in starvation. Cell 172, 234–248.e17 (2018).
CAS PubMed PubMed Central Article Google Scholar
Greenhaff, P. L. et al. Disassociation between the effects of amino acids and insulin on signaling, ubiquitin ligases, and protein turnover in human muscle. Am. J. Physiol. Metab. 295, E595–E604 (2008).
Chubukov, V. et al. Transcriptional regulation is insufficient to explain substrate‐induced flux changes in Bacillus subtilis. Mol. Syst. Biol. 9, 709 (2013).
CAS PubMed PubMed Central Article Google Scholar
Wolfe, R. R. & Chinkes, D. L. Isotope Tracers in Metabolic Research: Principles and Practice of Kinetic Analysis (Wiley-Liss, 2005).
Kim, I.-Y., Suh, S.-H., Lee, I.-K. & Wolfe, R. R. Applications of stable, nonradioactive isotope tracers in in vivo human metabolic research. Exp. Mol. Med. 48, e203 (2016).
CAS PubMed PubMed Central Article Google Scholar
Gasier, H. G., Fluckey, J. D. & Previs, S. F. The application of 2H2O to measure skeletal muscle protein synthesis. Nutr. Metab. 7, 31 (2010).
Cobelli, C., Foster, D. & Toffolo, G. Tracer Kinetics in Biomedical Research (Springer US, 2002).
Holmes, W. E., Angel, T. E., Li, K. W. & Hellerstein, M. K. Dynamic proteomics: in vivo proteome-wide measurement of protein kinetics using metabolic labeling. Methods Enzymol. 561, 219–276 (2015).
CAS PubMed Article Google Scholar
Wilkinson, D. J., Brook, M. S., Smith, K. & Atherton, P. J. Stable isotope tracers and exercise physiology: past, present and future. J. Physiol. 595, 2873–2882 (2017).
CAS PubMed Article Google Scholar
Jacobs, K. A., Casazza, G. A., Suh, S. H., Horning, M. A. & Brooks, G. A. Fatty acid reesterification but not oxidation is increased by oral contraceptive use in women. J. Appl. Physiol. 98, 1720–1731 (2005).
CAS PubMed Article Google Scholar
Volpi, E., Ferrando, A. A., Yeckel, C. W., Tipton, K. D. & Wolfe, R. R. Exogenous amino acids stimulate net muscle protein synthesis in the elderly. J. Clin. Invest. 101, 2000–2007 (1998).
CAS PubMed PubMed Central Article Google Scholar
Mason, A., Engelen, M. P. K. J., Ivanov, I., Toffolo, G. M. & Deutz, N. E. P. A four-compartment compartmental model to assess net whole body protein breakdown using a pulse of phenylalanine and tyrosine stable isotopes in humans. Am. J. Physiol. Metab. 313, E63–E74 (2017).
Song, B. et al. Mitoribosomal defects aggravate liver cancer via aberrant glycolytic flux and T cell exhaustion. J. Immunother. Cancer 10, e004337 (2022).
PubMed PubMed Central Article Google Scholar
Baker, N. & Schotz, M. C. Quantitative aspects of free fatty acid metabolism in the fasted rat. J. Lipid Res. 8, 646–660 (1967).
CAS PubMed Article Google Scholar
Wyss, M. & Kaddurah-Daouk, R. Creatine and creatinine metabolism. Physiol. Rev. 80, 1107–1213 (2000).
CAS PubMed Article Google Scholar
Kim, I.-Y., Park, S., Jang, J. & Wolfe, R. R. Understanding muscle protein dynamics: technical considerations for advancing sarcopenia research. Ann. Geriatr. Med. Res. 24, 157–165 (2020).
PubMed PubMed Central Article Google Scholar
Clark, R. V. et al. Total body skeletal muscle mass: estimation by creatine (methyl-d3) dilution in humans. J. Appl. Physiol. 116, 1605–1613 (2014).
CAS PubMed PubMed Central Article Google Scholar
Balsom, P. D., Söderlund, K. & Ekblom, B. Creatine in humans with special reference to creatine supplementation. Sports Med. 18, 268–280 (1994).
CAS PubMed Article Google Scholar
Gastaldelli, A. et al. Glucose kinetics in interstitial fluid can be predicted by compartmental modeling. Am. J. Physiol. 272, E494–E505 (1997).
Romijn, J. A. et al. Regulation of endogenous fat and carbohydrate metabolism in relation to exercise intensity and duration. Am. J. Physiol. 265, E380–E391 (1993).
Romijn, J. A., Coyle, E. F., Sidossis, L. S., Rosenblatt, J. & Wolfe, R. R. Substrate metabolism during different exercise intensities in endurance-trained women. J. Appl. Physiol. 88, 1707–1714 (2000).
CAS PubMed Article Google Scholar
Suh, S. H., Casazza, G. A., Horning, M. A., Miller, B. F. & Brooks, G. A. Effects of oral contraceptives on glucose flux and substrate oxidation rates during rest and exercise. J. Appl. Physiol. 94, 285–294 (2003).
CAS PubMed Article Google Scholar
Brooks, G. A. The precious few grams of glucose during exercise. Int. J. Mol. Sci. 21, 1–19 (2020).
Mazzeo, R. S., Brooks, G. A., Schoeller, D. A. & Budinger, T. F. Disposal of blood [1-13C]lactate in humans during rest and exercise. J. Appl. Physiol. 60, 232–241 (1986).
CAS PubMed Article Google Scholar
Kim, I. Y. et al. Short term elevation in dietary protein intake does not worsen insulin resistance or lipids in older adults with metabolic syndrome: a randomized-controlled trial. BMC Nutr. 3, 1–9 (2017).
Kim, I.-Y. et al. Acute lysine supplementation does not improve hepatic or peripheral insulin sensitivity in older, overweight individuals. Nutr. Metab. 11, 49 (2014).
Steele, R. Influences of glucose loading and of injected insulin on hepatic glucose output. Ann. N. Y. Acad. Sci. 82, 420–430 (1959).
CAS PubMed Article Google Scholar
Cahill, G. F. Starvation in man. N. Engl. J. Med. 282, 668–675 (2010).
Biolo, G., Gastaldelli, A., Zhang, X. J. & Wolfe, R. R. Protein synthesis and breakdown in skin and muscle: a leg model of amino acid kinetics. Am. J. Physiol. 267, E467–E474 (1994).
Zhang, X., Chinkes, D. L., Wu, Z., Martini, W. Z. & Wolfe, R. R. Fractional synthesis rates of DNA and protein in rabbit skin are not correlated. J. Nutr. 134, 2401–2406 (2004).
CAS PubMed Article Google Scholar
Martini, W. Z., Chinkes, D. L. & Wolfe, R. R. Quantification of DNA synthesis from different pathways in cultured human fibroblasts and myocytes. Metabolism 53, 128–133 (2004).
CAS PubMed Article Google Scholar
Hellerstein, M. K. New stable isotope–mass spectrometric techniques for measuring fluxes through intact metabolic pathways in mammalian systems: introduction of moving pictures into functional genomics and biochemical phenotyping. Metab. Eng. 6, 85–100 (2004).
CAS PubMed Article Google Scholar
Siler, S. Q., Neese, R. A., Christiansen, M. P. & Hellerstein, M. K. The inhibition of gluconeogenesis following alcohol in humans. Am. J. Physiol. Metab. 275, E897–E907 (1998).
Donnelly, K. L. et al. Sources of fatty acids stored in liver and secreted via lipoproteins in patients with nonalcoholic fatty liver disease. J. Clin. Invest. 115, 1343–1351 (2005).
CAS PubMed PubMed Central Article Google Scholar
Barrows, B. R., Timlin, M. T. & Parks, E. J. Spillover of dietary fatty acids and use of serum nonesterified fatty acids for the synthesis of VLDL-triacylglycerol under two different feeding regimens. Diabetes 54, 2668–2673 (2005).
CAS PubMed Article Google Scholar
Baumann, P. Q., Stirewalt, W. S., O’Rourke, B. D., Howard, D. & Nair, K. S. Precursor pools of protein synthesis: a stable isotope study in a swine model. Am. J. Physiol. 267, E203–E209 (1994).
Wolfe, R., Kim, I., Park, S., Moughan, P. & Ferrando, A. Advances in stable isotope tracer methodology part 2. New thoughts about and ‘old’ method: measurement of whole-body protein synthesis and breakdown in the fed state. J. Investig. Med. 68, 11–15 (2019).
Reeds, P. J., Hachey, D. L., Patterson, B. W., Motil, K. J. & Klein, P. D. VLDL apolipoprotein B-100, a potential indicator of the isotopic labeling of the hepatic protein synthetic precursor pool in humans: studies with multiple stable isotopically labeled amino acids. J. Nutr. 122, 457–466 (1992).
CAS PubMed Article Google Scholar
Martini, W. Z., Chinkes, D. L. & Wolfe, R. R. The intracellular free amino acid pool represents tracer precursor enrichment for calculation of protein synthesis in cultured fibroblasts and myocytes. J. Nutr. 134, 1546–1550 (2004).
CAS PubMed Article Google Scholar
Hellerstein, M. K. & Neese, R. A. Mass isotopomer distribution analysis: a technique for measuring biosynthesis and turnover of polymers. Am. J. Physiol. 263, E988–E1001 (1992).
Kelleher, J. K. & Masterson, T. M. Model equations for condensation biosynthesis using stable isotopes and radioisotopes. Am. J. Physiol. 262, E118–E125 (1992).
Hellerstein, M. K. & Neese, R. A. Mass isotopomer distribution analysis at eight years: theoretical, analytic, and experimental considerations. Am. J. Physiol. Metab. 276, E1146–E1170 (1999).
Chinkes, D. L., Aarsland, A., Rosenblatt, J. & Wolfe, R. R. Comparison of mass isotopomer dilution methods used to compute VLDL production in vivo. Am. J. Physiol. 271, E373–E383 (1996).
Wilkinson, D. J. et al. A validation of the application of D2O stable isotope tracer techniques for monitoring day-to-day changes in muscle protein subfraction synthesis in humans. Am. J. Physiol. Metab. 306, E571–E579 (2014).
Brook, M. S. et al. Skeletal muscle hypertrophy adaptations predominate in the early stages of resistance exercise training, matching deuterium oxide-derived measures of muscle protein synthesis and mechanistic target of rapamycin complex 1 signaling. FASEB J. 29, 4485–4496 (2015).
CAS PubMed Article Google Scholar
Voogt, J. N. et al. Measurement of very low rates of cell proliferation by heavy water labeling of DNA and gas chromatography/pyrolysis/isotope ratio–mass spectrometric analysis. Nat. Protoc. 2, 3058–3062 (2007).
CAS PubMed Article Google Scholar
Bederman, I. R., Foy, S., Chandramouli, V., Alexander, J. C. & Previs, S. F. Triglyceride synthesis in epididymal adipose tissue: contribution of glucose and non-glucose carbon sources. J. Biol. Chem. 284, 6101–6108 (2009).
CAS PubMed PubMed Central Article Google Scholar
Katanik, J. et al. Measuring gluconeogenesis using a low dose of 2H2O: advantage of isotope fractionation during gas chromatography. Am. J. Physiol. Metab. 284, E1043–E1048 (2003).
Ussing, H. H. The rate of protein renewal in mice and rats studied by means of heavy hydrogen. Acta Physiol. Scand. 2, 209–221 (1941).
Busch, R. et al. Measurement of protein turnover rates by heavy water labeling of nonessential amino acids. Biochim. Biophys. Acta - Gen. Subj. 1760, 730–744 (2006).
Bederman, I. R., Dufner, D. A., Alexander, J. C. & Previs, S. F. Novel application of the “doubly labeled” water method: measuring CO2 production and the tissue-specific dynamics of lipid and protein in vivo. Am. J. Physiol. Metab. 290, E1048–E1056 (2006).
Camera, D. M., Burniston, J. G., Pogson, M. A., Smiles, W. J. & Hawley, J. A. Dynamic proteome profiling of individual proteins in human skeletal muscle after a high-fat diet and resistance exercise. FASEB J. 31, 5478–5494 (2017).
CAS PubMed Article Google Scholar
Jang, J. et al. Myostatin inhibition-induced increase in muscle mass and strength was amplified by resistance exercise training, and dietary essential amino acids improved muscle quality in mice. Nutrients 13, 1508 (2021).
CAS PubMed PubMed Central Article Google Scholar
Kim, Y. et al. Essential amino acid-enriched diet alleviates dexamethasone-induced loss of muscle mass and function through stimulation of myofibrillar protein synthesis and improves glucose metabolism in mice. Metabolites 12, 84 (2022).
CAS PubMed PubMed Central Article Google Scholar
Shankaran, M. et al. Proteome-wide muscle protein fractional synthesis rates predict muscle mass gain in response to a selective androgen receptor modulator in rats. Am. J. Physiol. Endocrinol. Metab. 310, E405–E417 (2016).
Miller, B. F. et al. Lactate and glucose interactions during rest and exercise in men: effect of exogenous lactate infusion. J. Physiol. 544, 963–975 (2002).
CAS PubMed PubMed Central Article Google Scholar
Wang, Y., Kwon, H., Su, X. & Wondisford, F. E. Glycerol not lactate is the major net carbon source for gluconeogenesis in mice during both short and prolonged fasting. Mol. Metab. 31, 36–44 (2020).
CAS PubMed Article Google Scholar
Liu, S., Dai, Z., Cooper, D. E., Kirsch, D. G. & Locasale, J. W. Quantitative analysis of the physiological contributions of glucose to the TCA cycle. Cell Metab. 32, 619–628.e21 (2020).
CAS PubMed PubMed Central Article Google Scholar
Felig, P. The glucose-alanine cycle. Metabolism 22, 179–207 (1973).
CAS PubMed Article Google Scholar
Brooks, G. A. et al. Muscle accounts for glucose disposal but not blood lactate appearance during exercise after acclimatization to 4,300 m. J. Appl. Physiol. 72, 2435–2445 (1992).
CAS PubMed Article Google Scholar
Courtney, K. D. et al. Isotope tracing of human clear cell renal cell carcinomas demonstrates suppressed glucose oxidation in vivo. Cell Metab. 28, 793–800.e2 (2018).
CAS PubMed PubMed Central Article Google Scholar
Borghouts, L. B., Wagenmakers, A. J. M., Goyens, P. L. L. & Keizer, H. A. Substrate utilization in non-obese Type II diabetic patients at rest and during exercise. Clin. Sci. 103, 559–566 (2002).
Metges, C. C. et al. Kinetics of L-[1-13C]leucine when ingested with free amino acids, unlabeled or intrinsically labeled casein. Am. J. Physiol. Metab. 278, E1000–E1009 (2000).
Heux, S., Bergès, C., Millard, P., Portais, J. C. & Létisse, F. Recent advances in high-throughput 13C-fluxomics. Curr. Opin. Biotechnol. 43, 104–109 (2017).
CAS PubMed Article Google Scholar
Antoniewicz, M. R. A guide to 13C metabolic flux analysis for the cancer biologist. Exp. Mol. Med. 50, 1–13 (2018).
CAS PubMed Article Google Scholar
Xu, J. et al. Metabolic flux analysis and fluxomics-driven determination of reaction free energy using multiple isotopes. Curr. Opin. Biotechnol. 64, 151–160 (2020).
CAS PubMed Article Google Scholar
Antoniewicz, M. R. 13C metabolic flux analysis: optimal design of isotopic labeling experiments. Curr. Opin. Biotechnol. 24, 1116–1121 (2013).
CAS PubMed Article Google Scholar
Antoniewicz, M. R., Kelleher, J. K. & Stephanopoulos, G. Elementary Metabolite Units (EMU): a novel framework for modeling isotopic distributions. Metab. Eng. 9, 68–86 (2007).
CAS PubMed Article Google Scholar
Antoniewicz, M. R., Kelleher, J. K. & Stephanopoulos, G. Determination of confidence intervals of metabolic fluxes estimated from stable isotope measurements. Metab. Eng. 8, 324–337 (2006).
CAS PubMed Article Google Scholar
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MSIT) (2021R1A2C3005801), the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1A6A3A13071529), the Brain Pool program funded by the Ministry of Science and ICT through the National Research Foundation of Korea (2019H1D3A1A01071043), and the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1I1A1A01074380).
Department of Molecular Medicine, Lee Gil Ya Cancer and Diabetes Institute Gachon University School of Medicine, Incheon, 21999, Republic of Korea
Il-Young Kim & Sanghee Park
Korea Mouse Metabolic Phenotyping Center, Lee Gil Ya Cancer and Diabetes Institute, Gachon University, Incheon, 21999, Republic of Korea
Department of Health Sciences and Technology, GAIHST, Gachon University, Incheon, 21999, Republic of Korea
Yeongmin Kim & Hee-Joo Kim
Department of Geriatrics, the Center for Translational Research in Aging & Longevity, Donald W. Reynolds Institute on Aging, University of Arkansas for Medical Sciences, Little Rock, AR, USA
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
I.K. and R.R.W. conceptualized the manuscript; I.K. drafted the manuscript; H.K. prepared the figures; S.P., H.K., Y.K., and R.R.W. helped write the manuscript and I.K., S.P., H.K., Y.K., and R.R.W. discussed, edited, and approved the final manuscript.
The authors declare no competing interests.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
Kim, IY., Park, S., Kim, Y. et al. Tracing metabolic flux in vivo: basic model structures of tracer methodology. Exp Mol Med (2022). https://doi.org/10.1038/s12276-022-00814-z
DOI: https://doi.org/10.1038/s12276-022-00814-z
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
Tagging biomolecules with stable isotopes of specific atoms can reveal details of the molecular inter-conversions of metabolism. The masses of the tracer isotopes used are greater than those of the more common atomic forms. This allows their movement through different metabolic pathways to be detected using mass spectrometry and modeling. Il-Young Kim at Gachon University School of Medicine in South Korea and colleagues focus their review on the use of stable, non-radioactive isotope tracers, especially, of carbon, nitrogen, and hydrogen, to study metabolism in live humans and other animals. They cover the basic model structures of tracer methodology that serve as the fundamental basis for various tracer methods available and the most recent applications. Their procedure is especially useful for monitoring the rates of metabolic inter-conversions, which can reveal aspects of health and disease.
Experimental & Molecular Medicine (Exp Mol Med) ISSN 2092-6413 (online) ISSN 1226-3613 (print)