Articles | Volume 17, issue 19
Ideas and perspectives
02 Oct 2020
Ideas and perspectives |  | 02 Oct 2020

Ideas and perspectives: The same carbon behaves like different elements – an insight into position-specific isotope distributions

Yuyang He, Xiaobin Cao, and Huiming Bao

It is expected that information on the source, reaction pathway, and reaction kinetics of an organic compound can be obtained from its position-specific isotope compositions or intramolecular isotope distribution (Intra-ID). To retrieve the information, we could use its predicted equilibrium Intra-ID as a reference for understanding the observed Intra-IDs. Historically, observed, apparently close-to-equilibrium carbon Intra-ID has prompted an open debate on the nature of biosystems and specifically the pervasiveness of reversible biochemical reactions. Much of the debate remains unresolved, and the discussion has not clearly distinguished between two states of equilibrium: (1) the equilibrium among the corresponding bond-breaking and bond-forming positions in reactant and product and (2) the equilibrium among all carbon positions within a compound. For an organic molecule with multiple carbon positions, equilibrium carbon Intra-ID can be attained only when a specific reaction is in equilibrium and the sources of each position are also in equilibrium with each other. An observed Intra-ID provides limited information on if the sources and pathways are both unconstrained. Here, we elaborate on this insight using examples of the observed Intra-IDs of hydroxyl-bearing minerals, N2O, and acetic acid. Research effort aiming to calibrate position-specific equilibrium and kinetic isotope fractionation factors for defined processes will help to interpret observed Intra-IDs of a compound accurately and fully.

1 Introduction

Biosystems are dominated by a series of nonequilibrium kinetic processes. The understanding of biosystems is rooted in the study of the biochemical reaction mechanism. However, a majority of the biochemical reaction mechanisms remain elusive since they are difficult to isolate and control in laboratory experiments. Stable isotope effects can be used to examine the transition state structure and reversibility of an elementary reaction. Therefore, they can provide information on reaction mechanisms (Bigeleisen, 1949; Galimov, 2006; Bennet, 2012). However, a big organic molecule produced by an organism is the result of complex biochemical reactions that involve multiple kinetic isotope effects (KIEs) and equilibrium isotope effects (EIEs). KIE and EIE refer to the two intrinsic parameters for interpreting the observed isotope fractionations (Bao et al., 2015). According to the transition state theory (Eyring, 1935a, b), the KIE of an elementary step can be defined as the equilibrium fractionation factor between the transition state and reactant (Jones and Urbauer, 1991; Bao et al., 2015):

(1) KIE = β TS / β R ,

where the β factor denotes the reduced partition function ratio of the transition state (TS) or reactant (R). A β factor is the equilibrium isotope fractionation factor between an atom in a specific bond environment and its atomic form that can be predicted theoretically (Urey, 1947; Bigeleisen and Goeppert-Mayer, 1947). For a unidirectional reaction, the KIE of a reaction can also be defined as

(2) KIE = h k / l k ,

where k denotes the reaction rate constant of heavy (h) or light (l) isotopes. To adapt to the convention of geochemists, we define KIE this way so that the normal KIE is less than 1.000, which is the opposite of what Bigeleisen (1949) initially defined.

EIE is the isotope fractionation among reactant and product, which is determined by the bonding environment of the target position or compound:

(3) EIE = β P / β R ,

where P denotes the product of a target reaction. It can also be defined as

(4) EIE = h K / l K ,

where K denotes the equilibrium constant of a target reaction. At equilibrium, the EIE of a reaction equals the ratio of forward reaction KIEf and backward reaction KIEb (EIE=KIEf/KIEb; Bao et al., 2016).

An organic compound usually contains multiple positions of the same element, such as carbon, hydrogen, oxygen, or nitrogen. Compound-specific isotope composition refers to the bulk isotope composition of an element in an individual compound. Position-specific isotope composition refers to the isotope composition of an element at structurally distinct atomic positions within an individual compound. Information on sources, reaction pathways, and reaction kinetics of an organic molecule is pertinent to each position. The compound-specific isotope composition averages isotope compositions of all different positions of the same element in a compound, where information contained in position-specific isotope compositions could be lost (Elsner, 2010; Piasecki et al., 2018).

We named position-specific isotope compositions in a compound intramolecular isotope distribution or Intra-ID (He et al., 2018, 2020). Carbon Intra-ID in organic compounds has invoked a long-standing debate about its fundamental controls. When faced with the observed diverse Intra-IDs, earlier researchers inferred that the patterns “must be the expression of some logical order” (Schmidt, 2003), which is controlled by the EIE and KIE of biochemical reactions (e.g., Hayes, 2001, 2004; Galimov, 2009; Schmidt et al., 2015; Eiler et al., 2018; Gilbert et al., 2019). The Intra-ID was described as being in a “thermodynamic order” or “statistical isotope pattern” when each position in a molecule reaches equilibrium with each other (Galimov, 1985; Schmidt et al., 2015). Here, we name it equilibrium Intra-ID. The nonequilibrium state is expected to be a norm for a biochemical system since life is a dissipative system. At equilibrium, the difference in isotope composition between two positions depends on temperature only, and therefore the deviation of an observed Intra-ID from its predicted equilibrium state has been considered as an ideal reference for interpreting position-specific isotope compositions (Galimov, 1985; Hayes, 2001, 2004; He et al., 2018, 2020; Rustad, 2009; Piasecki et al., 2016).

It has been reported that different carbon fragments of chlorophyll, different carbon positions in acetoin, malonic acid, citric acid, and purine alkaloid have 13βδ13C correlation with regression coefficients in the range of 0.33–0.51 (Galimov, 1985, 2003, 2004, 2006, and references therein). Such a 13βδ13C correlation is written as δ13C-δ13Cave=χβ-βave×103, where χ is the regression coefficient. Galimov interpreted such observed intramolecular 13βδ13C correlations as “equilibrium-like” Intra-IDs produced from sets of reversible biochemical reactions at steady states which are not far from equilibrium. The 13βδ13C correlations were used as supporting evidence that the theorem of minimum entropy production can be applied in biochemical systems. However, other groups interpreted the fair-to-good correlation as fortuitous regardless of the presence or absence of complete reversibility of enzymatic reactions (Buchachenko, 2003, 2007; Schmidt, 2003; Schmidt et al., 2015). In contrast to these reported observed equilibrium-like Intra-IDs, measured position-specific δ13C values are poorly correlated with their predicted 13β values in organic molecules like glucose, nicotine, and tropine (Rossmann et al., 1991; Gleixner and Schmidt, 1997; Robins et al., 2016; Romek et al., 2016). Such an observed nonequilibrium Intra-ID has been termed a “non-statistical isotope pattern” (Rossmann et al., 1991; Gleixner and Schmidt, 1997; Schmidt, 2003; Robins et al., 2016; Romek et al., 2016). Buchachenko (2003, 2007) and Schmidt et al. (2004, 2015) argued that the observed 13βδ13C correlations are a random Intra-ID that only “simulates” the thermodynamic state, which cannot be used as evidence for biochemical reactions favoring an equilibrium state. The 13βδ13C correlation used an unweighted arithmetic mean isotope composition of all the components as the reference of a system. Strictly, only the mass-weighted isotope composition of all components should represent that of a system (Hayes, 2001). In addition, arbitrarily assigning a reference is not mathematically rigorous (He et al., 2018). Therefore, a 13βδ13C correlation cannot be used as supporting evidence for Galimov's hypothesis that the theorem of minimum entropy production applies in biochemical systems. Nevertheless, the invalidity of 13βδ13C correlations cannot fully quell the controversy on the nature of biosystems.

It should be noted that the debate on isotope equilibrium in biosystems among Galimov, Buchachenko, and Schmidt (and their colleagues) did not clearly distinguish between two states of equilibrium: (1) intermolecular isotope equilibrium among the corresponding bond-breaking and bond-forming positions in reactant and product in a defined process and (2) intramolecular isotope equilibrium among all carbon positions in a defined molecule. Such a difference might also be overlooked when discussing the Intra-ID or the site preference (SP) value, i.e., the isotope composition difference among two positions. A fully reversible reaction is necessary for isotope equilibrium between corresponding active positions or functional groups. Similarly, a fully reversible intramolecular exchange mechanism must exist if different positions within a compound are to attain equilibrium. However, an overwhelming majority of biochemical reactions, especially in cases involving large organic molecules, have very few intramolecular exchange pathways. Here, we propose that the utility of parameters like the SP value in organic molecules could be limited before we obtain sufficient details on the source and pathway, as well as on the KIE and EIE of biochemical reactions. To elaborate this point, we present simple cases starting from hydroxyl-bearing minerals, in which oxygen occupies more than one position, and going on to the case of N2O in which unidirectional and fully reversible reactions can produce similar nitrogen Intra-IDs if there exists a symmetric precursor. After presenting the two inorganic cases, we move to examine measured carbon Intra-IDs from the literature of a simple organic molecule, acetic acid, in which Intra-IDs are pathway-dependent.

2 Intramolecular isotope distribution

2.1 Intracrystalline oxygen isotope difference – a rarely effective single mineral geothermometer

The same element, e.g., carbon, occupying different positions in a compound is not a unique feature of organic compounds. Some oxygen-bearing minerals have two or more position-specific oxygens, where the oxygen atoms occupy different positions in a mineral structure and have different chemical properties. Their isotope composition difference has been proposed as a potential single-mineral geothermometer. For example, it has been proposed that water temperature could be reconstructed from intracrystalline oxygen isotope difference in single-mineral copper sulfate pentahydrate (CuSO45H2O; Götz et al., 1975), kaolinite (Al2Si2O5(OH)4), illite (K0.65Al2.0(Al0.65Si3.35O10)(OH)2; Bechtel and Hoernes, 1990), or alunite (KAl3(SO4)2(OH)6; Arehart et al., 1992). By analyzing the isotope composition difference in different oxygens in the same minerals, the early researchers attempted to reconstruct the precipitation temperatures.

To be a single-mineral geothermometer, different oxygen sites must have attained equilibrium within the single mineral, which can be achieved when different positions in a compound have the same source or initially different sources are in equilibrium with each other. Take alunite precipitation from a solution as an example. Alunite has sulfate and hydroxyl oxygen positions in its structure that precipitate from sulfate and hydroxyl ions in the solution (Fig. 1). Alunite with an equilibrium Intra-ID can be produced from an equilibrium precipitation process only if both the oxygen isotope compositions of sulfate and hydroxyl ions in the solution equilibrated with the same ambient water oxygen at the same temperature. However, the oxygen of sulfate does not readily exchange with that of water; the isotope equilibration time for SO42- and ambient water at Earth surface conditions is greater than 106 to 107 years (Lloyd, 1968; Turchyn and Schrag, 2004; Turchyn et al., 2010), while the oxygen in OH can equilibrate with ambient water instantly and can readily exchange during alunite's later burial and diagenetic processes. Thus, the two oxygen positions in alunite can come from different sources at different temperatures, rendering alunite a flawed single-mineral geothermometer. The same is true for gypsum (CaSO42H2O) in which sulfate oxygen is not in equilibrium with formation water and the crystallization water (2H2O) oxygen may be in equilibrium with a different type of water.

Figure 1Sketch of alunite precipitation from water. The alunite could be a single-mineral geothermometer if three conditions are all fully satisfied: (1) H2OA=H2OB, (2) T1=T2=T3, and (3) all the four reactions are fully reversible and attain equilibrium. White, red, yellow, pink, and purple spheres represent hydrogen, oxygen, sulfur, aluminum, and potassium atoms, respectively.


The use of a single-mineral geothermometer requires that oxygens at two different sites have attained equilibrium by exchanging with each other or with the same source oxygen (e.g., water). Unfortunately, this requirement is difficult to meet for most minerals. It is, therefore, not surprising that few successful single-mineral geothermometers exist if any at all.

2.2 Equilibrium-like Intra-ID produced by a kinetic process

For a compound with two different positions of the same element, a simple way to describe its observed Intra-ID is to report the difference between the two isotope compositions, i.e., the site preference (SP) value. The concept of SP originated from the study of nitrous oxide (βNαNO), which is defined as the nitrogen isotope composition difference between the central nitrogen (δ15Nα) and the terminal nitrogen (δ15Nβ) (Yoshida and Toyoda, 2000). The predicted equilibrium SP value at room temperature in N2O is 45 ‰ (Yung and Miller, 1997; Wang et al., 2004; Cao and Liu, 2012). Although most observations fit the equilibrium prediction that 15N preferentially enriches in the αN position by 30 ‰–40 ‰ (Yoshida and Toyoda, 2000; Toyoda et al., 2002; Sutka et al., 2006), negative SP values have been observed nevertheless (Yamulki et al., 2001; Sutka et al., 2003).

In previous literature, the difference in SP values was explained by the difference in synthetic pathways associated with symmetrical or asymmetrical precursors (Schmidt et al., 2004; Toyoda et al., 2005; Sutka et al., 2006). If the precursor of N2O is symmetrical (e.g., –ONNO–, Fig. 2 left), the two nitrogens in the precursor are positionally equivalent; any prior isotope composition and fractionation difference would be erased by the symmetrical structures of the precursor. When producing N2O from a symmetrical precursor, the βN undergoes N–O bond cleavage and therefore has a primary isotope effect which is large, whereas the αN has only a secondary isotope effect which is negligible (close to 1.000; Bigeleisen and Wolfsberg, 1958). Therefore, 15N depletion is expected only on the βN of N2O produced from a symmetrical precursor and is expected to have a positive SP value. If the precursor is asymmetrical (e.g., –NH(OH)NO, Fig. 2 right), the two nitrogens in the precursor are not positionally equivalent. It is assumed that the two nitrogens in the precursor were produced from different EIEs or KIEs because they went through different reaction pathways and may even have different nitrogen sources. Therefore, during the formation of N2O from an asymmetrical precursor, the difference in the position-specific δ15N values of the precursors and the difference in isotope fractionation during the formation processes will be recorded in the SP value of N2O. Such N2O can have either SP>0 or SP<0.

Figure 2Proposed mechanisms for N2O formation from symmetrical and asymmetrical precursors (modified from Schmidt et al., 2004). Light gray, red, blue, and purple spheres represent hydrogen, oxygen, nitrogen, and iron atoms, respectively.


Nevertheless, the two previously proposed mechanisms cannot distinguish N2O with SP>0 produced from the two mechanisms. In addition, for N2O produced from a symmetric precursor, the SP value cannot provide information on the reaction kinetics, since both fully reversible and unidirectional reactions can produce similar observed SP values. When we state that a compound displays an equilibrium Intra-ID, the underlying assumption is that there exists a mechanism for different positions to exchange isotopes intramolecularly. However, not all observed apparent equilibrium or equilibrium-like Intra-IDs are produced by an intramolecular equilibrium process. For reactions like –ONNO– N2O, two types of processes could produce SP>0. First, the N2O formation reaction is fully reversible and attains an equilibrium. When fully reversible, the two nitrogens in N2O are scrambled when it forms the symmetrical precursor through the reverse reaction. At equilibrium, the terminal nitrogen in a weaker-bond environment is expected to be depleted in heavier isotopes than the central nitrogen by 45 ‰ at surface temperature. Second, the N2O formation reaction is unidirectional. When unidirectional, only the N–O bond-breaking position (βN) undergoes a KIE. Thus, the SP value is approximately equal to the KIE value. In this scenario, if there is a normal KIE, the terminal nitrogen is expected to be depleted in heavier isotopes than the central nitrogen by the extent of the KIE value. Such an observed Intra-ID would be similar to the predicted equilibrium Intra-ID, but it is produced by isotope depletion on the unidirectional bond-breaking process. No intramolecular exchange is involved. Therefore, even if the N2O produced by the unidirectional process has SP≈45 ‰, it is not due to a close-to-equilibrium intramolecular isotope exchange. Therefore, it is necessary to distinguish between the mechanisms and reaction kinetics that can produce an observed Intra-ID.

Here we see that both fully reversible and unidirectional processes can result in a similar SP value, but the underlying mechanisms are entirely different. Furthermore, a positive SP value can also be achieved through a combination of nitrogen sources and isotope fractionations from an asymmetrical precursor. Thus, without knowing the underlying process, we cannot interpret an observed Intra-ID or SP value uniquely.

2.3 Position-specific isotope fractionations between reactant and product

As illustrated above, the observed Intra-ID of a compound can be used to gauge the degree of internal thermodynamic equilibrium only if we can determine the mechanisms involved in isotope fractionation. It does not mean, however, that position-specific isotope composition is useless. Based on the predicted equilibrium Intra-ID, a predicted isotope fractionation factor (α) of corresponding positions between the reactant and product in a process can help to evaluate the thermodynamic state of a system and to decipher reaction pathways. In this section, we use a simple organic molecule, acetic acid (CH3COOH), and its measured Intra-IDs from previous studies as examples to illustrate how position-specific isotope fractionation occurs between reactant and product.

The relative isotope enrichment between the carboxyl and methyl carbon in acetic acid is defined as ln13αcarb-met=ln(13Rcarb/13Rmet). 13R(=13C/12C) denotes the carbon isotope molar abundance ratio in a position. Our calculated equilibrium Intra-ID of acetic acid has the carboxyl carbon being 47.3 ‰ heavier than the methyl carbon at 25 C (ln13αcarb-met(eq)=47.3; He et al., 2020). The measured δ13Cmet values from literature can be lower, higher, or approximately equal to the δ13Ccarb values for acetic acids from biological, artificial, or hydrous pyrolysis samples (Table 1). The position-specific δ13C values of biological, artificial, or hydrous-pyrolysis-produced acetic acid are largely overlapping in δ13Cmetδ13Ccarb space. For the majority of biological acetic acids, the δ13Ccarb values are several per mille higher than the δ13Cmet values (Fig. 3 top; ln13αcarb-met=5.1±4.8, n=29), with two cases of ∼18 ‰ higher and one case of −2.2 ‰ lower in δ13Ccarb values. It is expected that the metabolic and catabolic pathways and carbon sources are limited for most natural acetic acid. Therefore, the ln 13αcarb-met value of 5.1±4.8 ‰ could be characteristic but not necessarily exclusive for biologically produced acetic acid. Artificial acetic acids have a very large range of ln 13αcarb-met values from −30.2 ‰ to 24.2 ‰ (Fig. 3 middle; 7.3±14.3 ‰, n=24). Biological and hydrous-pyrolysis-produced acetic acids do not have such negative ln 13αcarb-met values.

Table 1Measured Intra-ID of acetic acid from the literature.

Note that  ln13αcarb-met=ln(13Rcarb/13Rmet)×1000 is defined as the relative isotope enrichments between carboxyl and methyl carbon in acetic acid. 13R=13A12A, and A denotes carbon isotope abundance in a position. a Meinschein et al. (1974). b Rinaldi et al. (1974). c Gelwicks et al. (1989). d Hattori et al. (2011). e Nimmanwudipong et al. (2015). f Yamada et al. (2002). g Yamada et al. (2014). h Thomas et al. (2009). i Dias et al. (2002a). j Dias et al. (2002b).

Download XLSX

Except for the abovementioned features, the production of artificial and biological acetic acid has too many unconstrained parameters. Thus, our discussion will focus on the acetic acid derived from hydrous pyrolysis of oil-prone source rocks. The acetic acids produced from the hydrous pyrolysis of oil-prone source rocks have a ln 13αcarb-met value of 18.3±7.7 ‰ (n=22; Fig. 3 bottom). At 310–350 C, ln 13αcarb-met values of ∼30 ‰ from Mahogany shale or Black shale with a proposed mechanism of unidirectional pyrolysis of precursor acid forms (R-CH2COOHR+CH3COOH; Fig. 4; Dias et al., 2002b). If we consider only the primary KIE between the methylene carbon in R–*CH2COOH and the methyl carbon in acetic acid (*CH3COOH), it is expected that a unidirectional process would lead to a 13C depletion only at the methyl carbon position in acetic acid. The observed Intra-ID of the produced acetic acid should equal the δ13C value difference between the precursors minus the primary KIE value. The primary KIE value is expected to be more negative than the predicted equilibrium isotope fractionation factor, which is −14 ‰ (He et al., 2020). Thus, as long as the position-specific δ13C value difference between the methylene and carboxyl carbon in R−CH2COOH is greater than −14 ‰, the acetic acid produced from unidirectional pyrolysis of such precursor acid should have a carboxyl carbon with a higher δ13C value than that of the methyl carbon. If the carboxyl carbon in the precursor acid has a higher δ13C value than that of the methylene carbon, the pyrolysis process can easily produce acetic acid with an apparent ln 13αcarb-met value close to the predicted equilibrium Intra-ID. Such apparently “equilibrium-like” Intra-ID does not involve intramolecular exchange, but it is the product of unidirectional precursor acid pyrolysis.

3 Implications

Life sustains itself by feeding on negative entropy. Persistent efforts are devoted to describing living systems by rigorous mathematics. Boltzmann first considered living organisms from a thermodynamic perspective, and Schrödinger later applied equilibrium thermodynamics to living systems (Popovic, 2018). Those attempts were not pursued further, since, as we all know today, a living system is an open system that is not in thermodynamic equilibrium, i.e., a dissipative system. The establishment of nonequilibrium thermodynamics by Prigogine and his coworkers has guided researchers to the theorem of minimum entropy production in biological systems (Prigogine and Wiame, 1946). The theorem of minimum entropy production in biological systems states that although a biosystem has increasing entropy, it is usually stable in a steady state. If the system is displaced from the steady state, it tends to return to its original state since the entropy is minimal at this state. Since then, efforts in applying nonequilibrium thermodynamics to living systems have been continued with mixed success (Stoward, 1962; Schneider and Kay, 1994; Hayflick, 2007; Demirel, 2010; Barbacci et al., 2015; Gerber et al., 2016). The theorem of minimum entropy production applies only to linear thermodynamic systems. Therefore, it is necessary to demonstrate that the magnitude of the reaction rate on the scale of interest in a living system is linearly dependent on the driving force responsible for the reaction system. It is reasonable to assume that a complex interacting and constantly evolving nonlinear system is constructed by a series of synergistic reactions and there should exist local linearity, a local steady state, and even local equilibrium (Galimov, 2006). Local nonequilibrium of biochemical systems is potentially significant for the increasing complexity and ordering in the structure of life (Prigogine and Wiame, 1946; Galimov, 2006). Such a system should consist of a set of reversible reactions but not necessarily equilibrium reactions conjugated with energy supplies that are maintained in a steady state not far from equilibrium. Galimov (1985, 2004, 2006) argued that such a close-to-equilibrium steady state should be expressed as a tendency toward equilibrium inter- and intramolecular stable isotope distributions, i.e., a linear inter- and intramolecular 13βδ13C correlation with a regression coefficient smaller than but close to 1. In addition, such an equilibrium-like Intra-ID in organic molecules was proposed as a “special feature of biological systems”, which could be used as a criterion to identify biologically produced extraterrestrial organic molecules (Galimov, 2003). As we have illustrated above, the observed Intra-ID in organic molecules is the product of a set of equilibrium or nonequilibrium processes as well as their source isotope compositions. An observed Intra-ID itself cannot be used as conclusive evidence for the thermodynamic state of a system. Therefore, even if a compound does have a linear intramolecular 13βδ13C correlation with a slope of 1, it does not constitute supporting evidence for the existence of an equilibrium state among biochemical reactions in organisms. To apply nonequilibrium thermodynamics to living systems, further solid evidence is needed.

Figure 3ln 13αcarb-met values of biological, artificial, and hydrous-pyrolysis-produced acetic acid.


Figure 4Acetic acid produced from pyrolysis of precursor acid forms has an Intra-ID that is depleted in 13C in the methyl position. Dark gray, light gray, and red spheres represent carbon, hydrogen, and oxygen atoms, respectively.


A simple comparison of position-specific isotope compositions in one sample, e.g., ln 13αcarb-met values of one acetic acid sample, offers little information on the reaction mechanisms and reaction kinetics of the reaction it involves. Although the position-specific atoms are the same element, without an exchange mechanism, they behave independently as different elements. It would be helpful if we could consider the position-specific atoms independently. The isotope fractionation relationship of different elements in the same compound, i.e.,  (αA-1)/(αB-1), ln αA∕ln αB, or ΔδA∕ΔδB, (named the bonded isotope effect; He and Bao, 2019), is often used to characterize a reaction pathway, for instance, δD and δ18O in H2O (Dansgaard, 1964; Craig, 1961), δ15N and δ18O in NO3- (Casciotti and McIlvin, 2007; Wankel et al., 2009), δ34S and δ18O in SO42- (Antler et al., 2013), or δ13C and δD in organic compounds (Elsner, 2010; Palau et al., 2017). The isotope composition difference in different elements is only useful if the isotope fractionation relationships are considered and their isotope compositions are normalized; e.g., δ(15,18)=(δ15N-δ15Nm)-(15α-1/18α-1)×(δ18O-δ18Om), where δ15Nm and δ18Om are the average isotope composition in a given ocean water column (Sigman et al., 2005). The normalization procedure is necessary because the source isotope compositions can affect the values of the product. Similarly, if the same element at different positions have different sources, their source isotope composition difference must also be considered. Position-specific isotope research can build upon our understanding of the bonded isotope effect.

4 Conclusions

An organic compound usually has an element, e.g., carbon, at different positions and therefore has an Intra-ID. The deviation of an observed Intra-ID from its equilibrium state has been used to evaluate the thermodynamic state of a system. Our analysis of oxygen-bearing minerals, N2O, and acetic acids shows that both isotope sources and all reaction processes need to be in equilibrium to reach an intramolecular equilibrium state. However, such a condition is rarely satisfied. When different positions of the same element cannot exchange with each other, these different positions behave independently like different elements. Observed Intra-ID that is apparently similar to the equilibrium one can also be produced from a combination of different sources and unidirectional processes. Thus, an observed Intra-ID itself is not conclusive without adequate information on sources and reaction kinetics. Compared to position-specific isotope compositions, position-specific isotope fractionation of a defined process is more informative for identifying bond-breaking and bond-forming positions of a large molecule, for predicting its transition state structure, for evaluating the reversibility of a biochemical process, and for determining and qualifying a process in a complex system. All in all, an understanding of a reaction process at the molecular level will always be the first step required for the later sound and wide application of stable isotope composition.

Code and data availability

All data, models, and code generated or used during the study appear in the text.

Author contributions

All authors conceived of the presented idea. YH collected data from the literature, drafted the manuscript, and designed the figures. XC and HB provided critical feedback and helped shape the manuscript. All authors discussed the results and contributed to the final manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Financial support

This research has been supported by the China Postdoctoral Science Foundation (grant no. 2019M660811), the National Natural Science Foundation of China (grant no. 41490635), and the Chinese Academy of Sciences (grant no. XDB18010104).

Review statement

This paper was edited by Yakov Kuzyakov and reviewed by three anonymous referees.


Antler, G., Turchyn, A. V., Rennie, V., Herut, B., and Sivan, O.: Coupled sulfur and oxygen isotope insight into bacterial sulfate reduction in the natural environment, Geochim. Cosmochim. Ac., 118, 98–117,, 2013. 

Arehart, G. B., Kesler, S. E., O'Neil, J. R., and Foland, K.: Evidence for the supergene origin of alunite in sediment-hosted micron gold deposits, Nevada, Econ. Geol., 87, 263–270,, 1992. 

Bao, H., Cao, X., and Hayles, J. A.: The confines of triple oxygen isotope exponents in elemental and complex mass-dependent processes, Geochim. Cosmochim. Ac., 170, 39–50,, 2015. 

Bao, H., Cao, X., and Hayles, J. A.: Triple Oxygen Isotopes: Fundamental Relationships and Applications, Annu. Rev. Earth Pl. Sc., 44, 463–492,, 2016. 

Barbacci, A., Magnenet, V., and Lahaye, M.: Thermodynamical journey in plant biology, Front. Plant Sci., 6, 481,, 2015. 

Bechtel, A. and Hoernes, S.: Oxygen isotope fractionation between oxygen of different sites in illite minerals: a potential single-mineral thermometer, Contrib. Mineral. Petr., 104, 463–470,, 1990. 

Bennet, A. J.: Kinetic isotope effects for studying post-translational modifying enzymes, Curr. Opin. Chem. Biol., 16, 472–478,, 2012. 

Bigeleisen, J.: The Relative Reaction Velocities of Isotopic Molecules, J. Chem. Phys., 17, 675–678,, 1949. 

Bigeleisen, J. and Goeppert-Mayer, M.: Calculation of equilibrium constants for isotopic exchange reactions, J. Chem. Phys., 15, 261–267,, 1947. 

Bigeleisen, J. and Wolfsberg, M.: Theoretical and experimental aspects of isotope effects in chemical kinetics, Adv. Chem. Phys., 15–76,, 1958. 

Buchachenko, A. L.: Are chemical transformations in an enzyme-substrate complex reversible? The experience of the fractionation of isotopes in enzymatic reactions, Russ. J. Phys. Ch., 77, 1298–1302, 2003. 

Buchachenko, A. L.: On isotope fractionation in enzymatic reaction, Russ. J. Phys. Chem. A, 81, 836–836,, 2007. 

Cao, X. and Liu, Y.: Theoretical estimation of the equilibrium distribution of clumped isotopes in nature, Geochim. Cosmochim. Ac., 77, 292–303,, 2012. 

Casciotti, K. and McIlvin, M.: Isotopic analyses of nitrate and nitrite from reference mixtures and application to Eastern Tropical North Pacific waters, Mar. Chem., 107, 184–201,, 2007. 

Craig, H.: Isotopic variations in meteoric waters, Science, 133, 1702–1703,, 1961. 

Dansgaard, W.: Stable isotopes in precipitation, Tellus, 16, 436–468,, 1964. 

Demirel, Y.: Nonequilibrium thermodynamics modeling of coupled biochemical cycles in living cells, J. Non-Newton. Fluid, 165, 953–972,, 2010. 

Dias, R. F., Freeman, K. H., and Franks, S. G.: Gas chromatography–pyrolysis–isotope ratio mass spectrometry: a new method for investigating intramolecular isotopic variation in low molecular weight organic acids, Org. Geochem., 33, 161–168,, 2002a. 

Dias, R. F., Freeman, K. H., Lewan, M. D., and Franks, S. G.: δ13C of low-molecular-weight organic acids generated by the hydrous pyrolysis of oil-prone source rocks, Geochim. Cosmochim. Ac., 66, 2755–2769,, 2002b. 

Eiler, J. M., Clog, M., Lawson, M., Lloyd, M., Piasecki, A., Ponton, C., Xie, H., Lawson, M., Formolo, M. J., and Eiler, J. M.: The isotopic structures of geological organic compounds, in: From Source to Seep: Geochemical Applications in Hydrocarbon Systems, Geological Society of London, ISBN 9781786203663, 2018. 

Elsner, M.: Stable isotope fractionation to investigate natural transformation mechanisms of organic contaminants: principles, prospects and limitations, J. Environ. Monitor., 12, 2005–2031,, 2010. 

Eyring, H.: The activated complex and the absolute rate of chemical reactions, Chem. Rev., 17, 65–77,, 1935a. 

Eyring, H.: The activated complex in chemical reactions, J. Chem. Phys., 3, 107–115,, 1935b. 

Galimov, E. M.: The Biological Fractionation of Isotopes, Academic Press, Orlando,, 1985. 

Galimov, E. M.: Isotopic Criteria for Identification of Organic Carbon on Earth and Meteorites, Space Sci. Rev., 106, 249–262,, 2003. 

Galimov, E. M.: Phenomenon of life: between equilibrium and non-linearity, Origins Life Evol. B., 34, 599–613,, 2004. 

Galimov, E. M.: Phenomenon of life: between equilibrium and non-linearity, Origin and principles of evolution, Geochem. Int., 44, S1–S95,, 2006. 

Galimov, E. M.: Advances in geochemistry during the last four decades: A personal perspective, Appl. Geochem., 24, 1048–1051,, 2009. 

Gelwicks, J. T., Risatti, J. B., and Hayes, J.: Carbon isotope effects associated with autotrophic acetogenesis, Org. Geochem., 14, 441–446,, 1989. 

Gerber, S., Fröhlich, M., Lichtenberg-Fraté, H., Shabala, S., Shabala, L., and Klipp, E.: A thermodynamic model of monovalent cation homeostasis in the yeast Saccharomyces cerevisiae, PLoS Comput. Biol., 12, e1004703,, 2016. 

Gilbert, A., Sherwood Lollar, B., Musat, F., Giunta, T., Chen, S., Kajimoto, Y., Yamada, K., Boreham, C. J., Yoshida, N., and Ueno, Y.: Intramolecular isotopic evidence for bacterial oxidation of propane in subsurface natural gas reservoirs, P. Natl. A. Sci., 116, 6653–6658,, 2019. 

Gleixner, G. and Schmidt, H.-L.: Carbon isotope effects on the fructose-1, 6-bisphosphate aldolase reaction, origin for non-statistical 13C distributions in carbohydrates, J. Biol. Chem., 272, 5382–5387,, 1997. 

Götz, D., Heinzinger, K., and Klemm, A.: The Dependence of the Oxygen Isotope Fractionation in the Hydration Water of CuSO4 5 H2O on the Crystallization Temperature, Z. Naturforsch. A, 30, 1667–1674,, 1975. 

Hattori, R., Yamada, K., Kikuchi, M., Hirano, S., and Yoshida, N.: Intramolecular carbon isotope distribution of acetic acid in vinegar, J. Agr. Food Chem., 59, 9049–9053,, 2011. 

Hayes, J. M.: Fractionation of carbon and hydrogen isotopes in biosynthetic processes, Rev. Mineral Geochem., 43, 225–277,, 2001. 

Hayes, J. M.: Isotopic order, biogeochemical processes, and earth history – Goldschmidt Lecture, Davos, Switzerland, August 2002, Geochim. Cosmochim. Ac., 68, 1691–1700,, 2004. 

Hayflick, L.: Entropy explains aging, genetic determinism explains longevity, and undefined terminology explains misunderstanding both, PLoS Genetics, 3, e220,, 2007. 

He, Y. and Bao, H.: Predicting High-Dimensional Isotope Relationships from Diagnostic Fractionation Factors in Systems with Diffusional Mass Transfer, ACS Earth and Space Chemistry, 3, 120–128,, 2019. 

He, Y., Cao, X., Wang, J., and Bao, H.: Identifying apparent local stable isotope equilibrium in a complex non-equilibrium system, Rapid Commun. Mass Sp., 32, 306–310,, 2018. 

He, Y., Bao, H., and Liu, Y.: Predicting equilibrium intramolecular isotope distribution within a large organic molecule by the cutoff calculation, Geochim. Cosmochim. Ac., 269, 292–302,, 2020. 

Jones, J. P. and Urbauer, J. L.: Theoretical kinetic isotope effects for the hydride transfer from formate to carbon dioxide: a comparison of theory with experiment, J. Comput. Chem., 12, 1134–1141,, 1991. 

Lloyd, R. M.: Oxygen isotope behavior in the Sulfate-Water System, J. Geophys. Res., 73, 6099–6110,, 1968. 

Meinschein, W., Rinaldi, G., Hayes, J., and Schoeller, D.: Intramolecular isotopic order in biologically produced acetic acid, Biol. Mass. Spectrom., 1, 172–174,, 1974. 

Nimmanwudipong, T., Gilbert, A., Yamada, K., and Yoshida, N.: Analytical method for simultaneous determination of bulk and intramolecular 13C-isotope compositions of acetic acid, Rapid Commun. Mass Sp., 29, 2337–2340,, 2015. 

Palau, J., Shouakar-Stash, O., Hatijah Mortan, S., Yu, R., Rosell, M., Marco-Urrea, E., Freedman, D. L., Aravena, R., Soler, A., and Hunkeler, D.: Hydrogen Isotope Fractionation during the Biodegradation of 1, 2-Dichloroethane: Potential for Pathway Identification Using a Multi-element (C, Cl, and H) Isotope Approach, Environ. Sci. Technol., 51, 10526–10535,, 2017. 

Piasecki, A., Sessions, A., Peterson, B., and Eiler, J.: Prediction of equilibrium distributions of isotopologues for methane, ethane and propane using density functional theory, Geochim. Cosmochim. Ac., 190, 1–12, 2016. 

Piasecki, A., Sessions, A., Lawson, M., Ferreira, A., Neto, E. S., Ellis, G. S., Lewan, M. D., and Eiler, J. M.: Position-specific 13C distributions within propane from experiments and natural gas samples, Geochim. Cosmochim. Ac., 220, 110–124,, 2018. 

Popovic, M.: Living organisms from Prigogine's perspective: an opportunity to introduce students to biological entropy balance, J. Biol. Educ., 52, 294–300,, 2018. 

Prigogine, I. and Wiame, J. M.: Biologie et thermodynamique des phénomènes irréversibles, Experientia, 2, 451–453,, 1946. 

Rinaldi, G., Meinschein, W., and Hayes, J.: Intramolecular carbon isotopic distribution in biologically produced acetoin, Biol. Mass. Spectrom., 1, 415–417,, 1974. 

Robins, R. J., Romek, K. M., Remaud, G. S., and Paneth, P.: Non-statistical isotope fractionation as a novel “retro-biosynthetic” approach to understanding alkaloid metabolic pathways, Phytochem. Lett., 20, 499–506,, 2016. 

Romek, K. M., Remaud, G. S., Silvestre, V., Paneth, P., and Robins, R. J.: Non-statistical 13C fractionation distinguishes co-incident and divergent steps in the biosynthesis of the alkaloids nicotine and tropine, J. Biol. Chem., 291, 16620–16629,, 2016. 

Rossmann, A., Butzenlechner, M., and Schmidt, H. L.: Evidence for a nonstatistical carbon isotope distribution in natural glucose, Plant Physiol., 96, 609–614,, 1991. 

Rustad, J. R.: Ab initio calculation of the carbon isotope signatures of amino acids, Org. Geochem., 40, 720–723,, 2009. 

Schmidt, H. L.: Fundamentals and systematics of the non-statistical distributions of isotopes in natural compounds, Naturwissenschaften, 90, 537–552,, 2003. 

Schmidt, H. L., Werner, R. A., Yoshida, N., and Well, R.: Is the isotopic composition of nitrous oxide an indicator for its origin from nitrification or denitrification? A theoretical approach from referred data and microbiological and enzyme kinetic aspects, Rapid Commun. Mass Sp., 18, 2036–2040,, 2004. 

Schmidt, H. L., Robins, R. J., and Werner, R. A.: Multi-factorial in vivo stable isotope fractionation: causes, correlations, consequences and applications, Isot. Environ. Healt. S., 51, 155–199,, 2015. 

Schneider, E. D. and Kay, J. J.: Life as a manifestation of the second law of thermodynamics, Math. Comput. Model., 19, 25–48,, 1994. 

Sigman, D. M., Granger, J., DiFiore, P. J., Lehmann, M. M., Ho, R., Cane, G., and van Geen, A.: Coupled nitrogen and oxygen isotope measurements of nitrate along the eastern North Pacific margin, Global Biogeochem. Cy., 19, GB4022,, 2005. 

Stoward, P. J.: Thermodynamics of Biological Growth, Nature, 194, 977–978,, 1962. 

Sutka, R., Ostrom, N., Ostrom, P., Gandhi, H., and Breznak, J.: Nitrogen isotopomer site preference of N2O produced by Nitrosomonas europaea and Methylococcus capsulatus Bath, Rapid Commun. Mass Sp., 17, 738–745,, 2003. 

Sutka, R. L., Ostrom, N., Ostrom, P., Breznak, J., Gandhi, H., Pitt, A., and Li, F.: Distinguishing nitrous oxide production from nitrification and denitrification on the basis of isotopomer abundances, Appl. Environ. Microb., 72, 638–644,, 2006. 

Thomas, B., Freeman, K. H., and Arthur, M. A.: Intramolecular carbon isotopic analysis of acetic acid by direct injection of aqueous solution, Org. Geochem., 40, 195–200,, 2009. 

Toyoda, S., Yoshida, N., Miwa, T., Matsui, Y., Yamagishi, H., Tsunogai, U., Nojiri, Y., and Tsurushima, N.: Production mechanism and global budget of N2O inferred from its isotopomers in the western North Pacific, Geophys. Res. Lett., 29,, 2002. 

Toyoda, S., Mutobe, H., Yamagishi, H., Yoshida, N., and Tanji, Y.: Fractionation of N2O isotopomers during production by denitrifier, Soil Biol. Biochem., 37, 1535–1545,, 2005. 

Turchyn, A. V. and Schrag, D. P.: Oxygen isotope constraints on the sulfur cycle over the past 10 million years, Science, 303, 2004–2007,, 2004. 

Turchyn, A. V., Brüchert, V., Lyons, T. W., Engel, G. S., Balci, N., Schrag, D. P., and Brunner, B.: Kinetic oxygen isotope effects during dissimilatory sulfate reduction: a combined theoretical and experimental approach, Geochim. Cosmochim. Ac., 74, 2011–2024,, 2010. 

Urey, H. C.: The thermodynamic properties of isotopic substances, J. Chem. Soc. (Resumed), 562–581,, 1947. 

Wang, Z., Schauble, E. A., and Eiler, J. M.: Equilibrium thermodynamics of multiply substituted isotopologues of molecular gases, Geochim. Cosmochim. Ac., 68, 4779–4797,, 2004.  

Wankel, S. D., Kendall, C., and Paytan, A.: Using nitrate dual isotopic composition (δ15N and δ18O) as a tool for exploring sources and cycling of nitrate in an estuarine system: Elkhorn Slough, California, J. Geophys. Res.-Biogeo., 114, GB01011,, 2009. 

Yamada, K., Tanaka, M., Nakagawa, F., and Yoshida, N.: On-line measurement of intramolecular carbon isotope distribution of acetic acid by continuous-flow isotope ratio mass spectrometry, Rapid Commun. Mass Sp., 16, 1059–1064,, 2002. 

Yamada, K., Kikuchi, M., Gilbert, A., Yoshida, N., Wasano, N., Hattori, R., and Hirano, S.: Evaluation of commercially available reagents as a reference material for intramolecular carbon isotopic measurements of acetic acid, Rapid Commun. Mass Sp., 28, 1821–1828,, 2014. 

Yamulki, S., Toyoda, S., Yoshida, N., Veldkamp, E., Grant, B., and Bol, R.: Diurnal fluxes and the isotopomer ratios of N2O in a temperate grassland following urine amendment, Rapid Commun. Mass Sp., 15, 1263–1269,, 2001. 

Yoshida, N. and Toyoda, S.: Constraining the atmospheric N2O budget from intramolecular site preference in N2O isotopomers, Nature, 405, 330–334,, 2000. 

Yung, Y. L. and Miller, C. E.: Isotopic fractionation of stratospheric nitrous oxide, Science, 278, 1778–1780,, 1997. 

Short summary
Different carbon sites in a large organic molecule have different isotope compositions. Different carbon sites may not have the chance to exchange isotopes at all. The lack of appreciation of this notion might be blamed for an unsettled debate on the thermodynamic state of an organism. Here we demonstrate using minerals, N2O, and acetic acid that the dearth of exchange among different carbon sites renders them as independent as if they were different elements in organic molecules.
Final-revised paper