![]() | This article may be too technical for most readers to understand.(September 2014) |
Transient kinetic isotope effects (or fractionation) occur when the reaction leading to isotope fractionation does not follow pure first-order kinetics (FOK) and therefore isotopic effects cannot be described with the classical equilibrium fractionation equations or with steady-state kinetic fractionation equations (also known as the Rayleigh equation). [1] In these instances, the general equations for biochemical isotope kinetics (GEBIK) and the general equations for biochemical isotope fractionation (GEBIF) can be used.
The GEBIK and GEBIF equations are the most generalized approach to describe isotopic effects in any chemical, catalytic reaction and biochemical reactions because they can describe isotopic effects in equilibrium reactions, kinetic chemical reactions and kinetic biochemical reactions. [2] In the latter two cases, they can describe both stationary and non-stationary fractionation (i.e., variable and inverse fractionation). In general, isotopic effects depend on the number of reactants and on the number of combinations resulting from the number of substitutions in all reactants and products. Describing with accuracy isotopic effects, however, depends also on the specific rate law used to describe the chemical or biochemical reaction that produces isotopic effects. Normally, regardless of whether a reaction is purely chemical or whether it involves some enzyme of biological nature, the equations used to describe isotopic effects base on FOK. This approach systematically leads to isotopic effects that can be described by means of the Rayleigh equation. In this case, isotopic effects will always be expressed as a constant, hence will not be able to describe isotopic effects in reactions where fractionation and enrichment are variable or inverse during the course of a reaction. Most chemical reactions do not follow FOK; neither biochemical reactions can normally be described with FOK. To properly describe isotopic effects in chemical or biochemical reactions, different approaches must be employed such as the use of Michaelis–Menten reaction order (for chemical reactions) or coupled Michaelis–Menten and Monod reaction orders (for biochemical reactions). However, conversely to Michaelis–Menten kinetics, GEBIK and GEBIF equations are solved under the hypothesis of non-steady state. This characteristic allows GEBIK and GEBIF to capture transient isotopic effects.
The GEBIK and GEBIF equations are introduced here below.
The GEBIK and GEBIF equations describe the dynamics of the following state variables
Both S and P contain at least one isotopic expression of a tracer atom. For instance, if the carbon element is used as a tracer, both S and P contain at least one C atom, which may appear as and . The isotopic expression within a molecule is
where is the number of tracer atoms within S, while is the number of isotopic substitutions in the same molecule. The condition must be satisfied. For example, the product in which 1 isotopic substitution occurs (e.g., ) will be described by .
Substrates and products appear in a chemical reaction with specific stoichiometric coefficients. When chemical reactions comprise combinations of reactants and products with various isotopic expressions, the stoichiometric coefficients are functions of the isotope substitution number. If and are the stoichiometric coefficient for substrate and product, a reaction takes the form
For example, in the reaction , the notation is with for both isotopologue reactants of the same substrate with substitution number and , and with for and because the reaction does not comprise production of and .
For isotopomers, the substitution location is taken into account as and , where and indicate a different expressions of the same isotopologue . Isotopomers only exist when and . The substitution location has to be specifically defined depending on the number of tracer atoms a, number of substitutions b, and molecule structure. For multiatomic molecules that are symmetric with respect to tracer position, there is no need to specify the substitution position when . For example, one substitution of deuterium [lower-alpha 1] in the symmetric methane molecule does not require the use of the right superscript. In the case that , the substitution location has to be specified, while for and it is not required. For example, two 2H substitutions in can occur in adjacent or non-adjacent locations. Using this notation, the reaction can be written as
where in defines only one of the two methane forms (either with adjacent or non-adjacent D atoms). The location of D in the two isotopologue water molecules produced on the right-hand side of the reaction has not been indicated because D is present in only one water molecule at saturation, and because the water molecule is symmetric. For asymmetric and multiatomic molecules with and , definition of the substitution location is always required. For instance, the isotopomers of the (asymmetric) nitrous oxide molecule are and .
Reactions of asymmetric isotopomers can be written using the partitioning coefficient u as
where . For example, using N isotope tracers, the isotopomer reactions
can be written as one reaction in which each isotopomer product is multiplied by its partition coefficient as
with . More generally, the tracer element does not necessarily occur in only one substrate and one product. If substrates react releasing products, each having an isotopic expression of the tracer element, then the generalized reaction notation is
(1) |
For instance, consider the and tracers in the reaction
In this case the reaction can be written as
with two substrates and two products without indication of the substitution location because all molecules are symmetric.
Biochemical kinetic reactions of type ( 1 ) are often catalytic reactions in which one or more substrates, , bind to an enzyme, E, to form a reversible activated complex, C, which releases one or more products, , and free, unchanged enzyme. These reactions belong to the type of reactions that can be described by Michaelis–Menten kinetics. Using this approach for substrate and product isotopologue and isotopomer expressions, and under the prescribed stoichiometric relationships among them, leads to the general reactions of the Michaelis–Menten type
(2) |
with the index , where m depends on the number of possible atomic combinations among all isotopologues and isotopomers. Here, , , and are the rate constants indexed for each of the m reactions.
The reactions
can be written as
The following isotope mass balances must hold
To solve for the concentration of all components appearing in any general biochemical reaction as in ( 2 ), the Michaelis–Menten kinetics for an enzymatic reaction are coupled with the Monod kinetics for biomass dynamics. The most general case is to assume that the enzyme concentration is proportional to the biomass concentration and that the reaction is not in quasi-steady state. These hypotheses lead to the following system of equations
(3a) |
(3b) |
(3c) |
(3d) |
(3e) |
with , and where is the concentration of the most limiting substrate in each reaction i, z is the enzyme yield coefficient, Y is the yield coefficient expressing the biomass gain per unit of released product and is the biomass mortality rate. [4]
The isotopic composition of the components in a biochemical system can be defined in different ways depending on the definition of isotopic ratio. Three definitions are described here:
Isotopic ratio relative to each component in the system, each with its isotopic expression, with respect to the concentration of its most abundant isotopologue
Isotopic ratio relative to the mass of the tracer element in each component;
where, and are the molecular weight of each isotopic expression of the substrate and product.
Isotopic ratio relative to the mass of the tracer element in the accumulated substrates and products
Regardless of the definition of the isotopic ratio, the isotopic composition of substrate and product are expressed as
, | (4a) |
. | (4a) |
where is a standard isotopic ration. Here, definition 3 of isotopic ratio has been used, however, any of the three definitions of isotopic ratio can equally be used.
The isotopic ratio of the product can be used to define the instantaneous isotopic ratio
(5) |
and the time-dependent fractionation factor
(6) |
The time-dependent isotopic enrichment is simply defined as
(7) |
Under specific assumptions, the GEBIK and GEBIF equations become equivalent to the equation for steady-state kinetic isotope fractionation in both chemical and biochemical reactions. Here two mathematical treatments are proposed: (i) under biomass-free and enzyme-invariant (BFEI) hypothesis and (ii) under quasi-steady-state (QSS) hypothesis.
In instances where the biomass and enzyme concentrations are not appreciably changing in time, we can assume that biomass dynamics is negligible and that the total enzyme concentration is constant, and the GEBIK equations become
(8a) |
(8b) |
(8c) |
(8d) |
Eqs. ( 4 ) for isotopic compositions, Eq. ( 6 ) for the fractionation factor and Eq. ( 7 ) for the enrichment factor equally applies to the GEBIK equations under the BFEI hypothesis.
If the quasi-steady-state hypothesis is assumed in addition to BFEI hypothesis, then the complex concentration can be assumed to be in a stationary (steady) state according to the Briggs–Haldane hypothesis, and the GEBIK equations become
(9a) |
(9a) |
which are written in a form similar to the classical Micaelis-Menten equations for any substrate and product. Here, the equations also show that the various isotopologue and isotopomer substrates appear as competing species. Eqs. ( 4 ) for isotopic compositions, Eq. ( 6 ) for the fractionation factor and Eq. ( 7 ) for the enrichment factor equally applies to the GEBIK equations under the BFEI and QSS hypothesis.
An example is shown where GEBIK and GEBIF equations are used to describe the isotopic reactions of consumption into according to the simultaneous set of reactions
These can be rewritten using the notation introduced before as.
The substrate has not been included due to its scarcity. In addition, we have not specified the isotopic substitution in the product of the second and third reactions because is symmetric. Assuming that the second and third reactions have identical reaction rates , , and , the full GEBIK and GEBIF equations are
The same reaction can be described with the GEBIK and GEBIF equations under the BFEI and QSS approximations as
where has been substituted with because the rate constants in the third reaction have been assumed to equal those of the second reaction.
In a chemical reaction, chemical equilibrium is the state in which both the reactants and products are present in concentrations which have no further tendency to change with time, so that there is no observable change in the properties of the system. This state results when the forward reaction proceeds at the same rate as the reverse reaction. The reaction rates of the forward and backward reactions are generally not zero, but they are equal. Thus, there are no net changes in the concentrations of the reactants and products. Such a state is known as dynamic equilibrium.
In probability theory and statistics, the gamma distribution is a versatile two-parameter family of continuous probability distributions. The exponential distribution, Erlang distribution, and chi-squared distribution are special cases of the gamma distribution. There are two equivalent parameterizations in common use:
In thermodynamics, the Onsager reciprocal relations express the equality of certain ratios between flows and forces in thermodynamic systems out of equilibrium, but where a notion of local equilibrium exists.
In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.
Compartmental models are a very general modelling technique. They are often applied to the mathematical modelling of infectious diseases. The population is assigned to compartments with labels – for example, S, I, or R,. People may progress between compartments. The order of the labels usually shows the flow patterns between the compartments; for example SEIS means susceptible, exposed, infectious, then susceptible again.
The equilibrium constant of a chemical reaction is the value of its reaction quotient at chemical equilibrium, a state approached by a dynamic chemical system after sufficient time has elapsed at which its composition has no measurable tendency towards further change. For a given set of reaction conditions, the equilibrium constant is independent of the initial analytical concentrations of the reactant and product species in the mixture. Thus, given the initial composition of a system, known equilibrium constant values can be used to determine the composition of the system at equilibrium. However, reaction parameters like temperature, solvent, and ionic strength may all influence the value of the equilibrium constant.
In mathematics, the Gaussian or ordinary hypergeometric function2F1(a,b;c;z) is a special function represented by the hypergeometric series, that includes many other special functions as specific or limiting cases. It is a solution of a second-order linear ordinary differential equation (ODE). Every second-order linear ODE with three regular singular points can be transformed into this equation.
The principle of detailed balance can be used in kinetic systems which are decomposed into elementary processes. It states that at equilibrium, each elementary process is in equilibrium with its reverse process.
In natural language processing, latent Dirichlet allocation (LDA) is a Bayesian network for modeling automatically extracted topics in textual corpora. The LDA is an example of a Bayesian topic model. In this, observations are collected into documents, and each word's presence is attributable to one of the document's topics. Each document will contain a small number of topics.
In the following we solve the second-order differential equation called the hypergeometric differential equation using Frobenius method, named after Ferdinand Georg Frobenius. This is a method that uses the series solution for a differential equation, where we assume the solution takes the form of a series. This is usually the method we use for complicated ordinary differential equations.
In mathematics, the structure constants or structure coefficients of an algebra over a field are the coefficients of the basis expansion of the products of basis vectors. Because the product operation in the algebra is bilinear, by linearity knowing the product of basis vectors allows to compute the product of any elements . Therefore, the structure constants can be used to specify the product operation of the algebra. Given the structure constants, the resulting product is obtained by bilinearity and can be uniquely extended to all vectors in the vector space, thus uniquely determining the product for the algebra.
The Luttinger–Kohn model is a flavor of the k·p perturbation theory used for calculating the structure of multiple, degenerate electronic bands in bulk and quantum well semiconductors. The method is a generalization of the single band k·p theory.
In continuum mechanics, a compatible deformation tensor field in a body is that unique tensor field that is obtained when the body is subjected to a continuous, single-valued, displacement field. Compatibility is the study of the conditions under which such a displacement field can be guaranteed. Compatibility conditions are particular cases of integrability conditions and were first derived for linear elasticity by Barré de Saint-Venant in 1864 and proved rigorously by Beltrami in 1886.
Equilibrium chemistry is concerned with systems in chemical equilibrium. The unifying principle is that the free energy of a system at equilibrium is the minimum possible, so that the slope of the free energy with respect to the reaction coordinate is zero. This principle, applied to mixtures at equilibrium provides a definition of an equilibrium constant. Applications include acid–base, host–guest, metal–complex, solubility, partition, chromatography and redox equilibria.
In mathematics, Jacobi polynomials are a class of classical orthogonal polynomials. They are orthogonal with respect to the weight on the interval . The Gegenbauer polynomials, and thus also the Legendre, Zernike and Chebyshev polynomials, are special cases of the Jacobi polynomials.
In mathematics, the Kodaira–Spencer map, introduced by Kunihiko Kodaira and Donald C. Spencer, is a map associated to a deformation of a scheme or complex manifold X, taking a tangent space of a point of the deformation space to the first cohomology group of the sheaf of vector fields on X.
In physics, relativistic angular momentum refers to the mathematical formalisms and physical concepts that define angular momentum in special relativity (SR) and general relativity (GR). The relativistic quantity is subtly different from the three-dimensional quantity in classical mechanics.
In theoretical physics, the dual graviton is a hypothetical elementary particle that is a dual of the graviton under electric-magnetic duality, as an S-duality, predicted by some formulations of eleven-dimensional supergravity.
The hyperbolastic functions, also known as hyperbolastic growth models, are mathematical functions that are used in medical statistical modeling. These models were originally developed to capture the growth dynamics of multicellular tumor spheres, and were introduced in 2005 by Mohammad Tabatabai, David Williams, and Zoran Bursac. The precision of hyperbolastic functions in modeling real world problems is somewhat due to their flexibility in their point of inflection. These functions can be used in a wide variety of modeling problems such as tumor growth, stem cell proliferation, pharma kinetics, cancer growth, sigmoid activation function in neural networks, and epidemiological disease progression or regression.