Detailed kinetics and regulation of mammalian 2-oxoglutarate dehydrogenase
© Qi et al; licensee BioMed Central Ltd. 2011
Received: 7 June 2011
Accepted: 26 September 2011
Published: 26 September 2011
Mitochondrial 2-oxoglutarate (α-ketoglutarate) dehydrogenase complex (OGDHC), a key regulatory point of tricarboxylic acid (TCA) cycle, plays vital roles in multiple pathways of energy metabolism and biosynthesis. The catalytic mechanism and allosteric regulation of this large enzyme complex are not fully understood. Here computer simulation is used to test possible catalytic mechanisms and mechanisms of allosteric regulation of the enzyme by nucleotides (ATP, ADP), pH, and metal ion cofactors (Ca2+ and Mg2+).
A model was developed based on an ordered ter-ter enzyme kinetic mechanism combined with con-formational changes that involve rotation of one lipoic acid between three catalytic sites inside the enzyme complex. The model was parameterized using a large number of kinetic data sets on the activity of OGDHC, and validated by comparison of model predictions to independent data.
The developed model suggests a hybrid rapid-equilibrium ping-pong random mechanism for the kinetics of OGDHC, consistent with previously reported mechanisms, and accurately describes the experimentally observed regulatory effects of cofactors on the OGDHC activity. This analysis provides a single consistent theoretical explanation for a number of apparently contradictory results on the roles of phosphorylation potential, NAD (H) oxidation-reduction state ratio, as well as the regulatory effects of metal ions on ODGHC function.
OGDHC was first purified from the pig heart mitochondria by Sanadi et al.  and subsequently studied by many researchers to examine its catalytic and regulatory properties within permeabilized, un-coupled, and coupled mitochondria from a variety of mammalian tissues [5–11]. A catalytic mechanism for the overall reaction of the enzyme complex was also first proposed by Sanadi et al.  which suggested that the coenzyme, NAD-, and 2-oxoglutaric acid participate in the reaction with the help of the cofactors thiamine pyrophosphate (TPP), lipoic acid, and FAD2-[12, 13]. Their proposed mechanism is a Hexa-Uni-Ping-Pong mechanism in Cleland's terminology  where it is assumed that the first product (CO2) is released before the second substrate (CoASH4-) binds, and the second product (Succinyl-CoA4-) is released before the third substrate (NAD-) binds to the enzyme. Subsequently, Koike et al.  postulated another mechanism in which, the lipoic acids transfer intermediates by rotating between the three catalytic sites. Furthermore, experimental results of fluorescence resonance energy transfer and dynamic anisotropy showed that the lipoic acids in the E2 component undergo motion where they rotate between different catalytic sites [16–18]. The results of steady-state kinetic studies done by Hamada et al.  and Smith et al.  contradict each other, and not all results are compatible with the Sanadi mechanism . This issue was addressed by Mcminn and Ottaway  with kinetic studies based on the Fromm method . Mcminn and Ottaway  explained the observed nonlinearity in the reciprocal plots of the results and proposed a phenomenological mechanism with semi-random characteristic. A recent study by Aevarsson et al.  on the crystal structure and architecture of 2-oxo acid dehydrogenase multi-enzyme complexes, provides interesting insights into the plausible kinetic mechanism of 2-oxo acid dehydrogenase family which includes OGDHC.
It has been consistently shown that the activity of OGDHC is controlled by various factors, including the variations of the NAD oxidation-reduction state, the state of phosphorylation of the nucleotide systems, and the ratio of succinyl-CoA to CoA-SH. Regulation by reversible phosphorylation has not been demonstrated. Experiments in isolated mitochondria of heart, liver, and kidney have shown that the OGDHC is regulated by Ca2+ ions with a marked decrease in the apparent K m for 2-oxoglutarate in the presence of adenine nucleotides and minimal effect of Ca2+ at saturating concentration of 2-oxoglutarate [9, 23–25]. Moreover, the apparent K m for 2-oxoglutarate is lowered by a decrease in the ATP/ADP ratio, which can significantly increase the sensitivity of the enzyme to Ca2+ions [23, 25]. It has been reported that the maximum activity of OGDHC is unaffected by changes in pH, while the apparent K m of the enzyme for 2-oxoglutarate is greatly altered by changes in pH over the range of 6.5 - 7.5 . A number of studies have also demonstrated the possible role of Mg2+ ions in the regulation of OGDHC either by directly affecting the activity of the enzyme or by modulating the Ca2+ effect on the enzyme. Mg2+ ion has been shown to increase [26–28] or to have no effect on the activity of OGDHC [29, 30]. McCormack and Denton  studied isolated OGDHC from pig heart mitochondria and found that there is no effect of EDTA and 1 mM Mg2+ on the activity of OGDHC when Ca2+ concentration was effectively less that 1 μM. Panov and Scarpa  concluded that the effects of Mg2+ and Ca2+ ions on the OGDHC activity are additive only at relative low concentration of free cations which suggested that at high concentrations, each ion may compete each other for binding sites. It is also evident that, in the presence of low Ca2+ concentration, Mg2+ ion can strongly modify the enzyme's affinities for 2-oxoglutarate and NAD-. However, the kinetic mechanisms by which these divalent metal ions regulate the properties of mitochon-drial OGDHC are not understood.
Although a number of attempts have been made to understand the catalytic mechanisms of OGDHC, both experimentally and theoretically, there is no mechanistic model that consistently explains the available experimental data on the kinetics of this enzyme complex and adequately describes the regulatory roles of nucleotides and other metal ion cofactors (Ca2+, Mg2+, etc.). Therefore a mechanistic model of OGDHC is needed to understand the orchestrated controlling of OGDHC by cofactors inside mitochondria under different physiological conditions. In the present work, a kinetic model of OGDHC is introduced to quantitatively understand the catalytic properties and regulation of OGDHC, based on the observations from a large number of independent experimental studies in mammalian tissues. The model accurately describes the catalytic properties of this enzyme complex observed experimentally, and clarifies many contradictory results reported in earlier studies.
In this section, we first present a general kinetic model for conformational changes in OGDHC, based on a presumed ter-ter enzyme mechanism via substrate channeling. The model is then used to characterize the kinetics of the 2-oxoglutarate dehydrogenase reaction (Equation 1) and further extended to describe the regulatory roles of cofactors, i.e., nucleotides and various metal ions. The kinetic parameters of the model are estimated using a wide variety of experimental data, available in the literature.
Kinetic scheme for conformational changes in a ter-ter enzyme mechanism
The kinetic equation of the proposed model for OGDHC reaction is derived from a ter-ter enzyme me-chanism combined with a model of conformational changes that represent the rotation of the single lipoic acid between different catalytic sites [15, 18]. The derivation is inspired by a previously developed model for trans-carboxylase . The assumption of the model is that the enzyme complex is composed of three sub-enzyme (E1, E2, and E3), each with one binding site: site 1 binds to 2-oxoglutarate (α KG2-) or corresponding product CO2, site 2 binds to CoASH4- or corresponding product Succinyl-CoA4-, site 3 binds toNAD- or corresponding product NADH2- (Figure 1A). Furthermore, the basic mechanism involves conformational changes, where the rotation of one lipoic acid between three catalytic sites leads to transfer of succinyl from E1 to E2 and proton from E2 to E3. In the catalytic cycle, the disulfide at the tip of the lipoyl can be in oxidized, reduced or semi-reduced li-poate forms, the semi-reduced form is bound with succinyl residue transferred from 2-oxoglutarate (Figure 1B).
where denotes the fraction of free enzyme complex that binds to the reactants; K A , K B , K C , K P , K Q and K R are the dissociation constants associated with the binding of reactants (A: α KG2-, B: CoASH4- and C: NAD-) and products (P: CO2, Q: Succinyl-CoA4- and R: NADH2-) to the enzyme complex. Here, we assumed that these constants do not depend on conformational states of the enzyme complex and the bound reactant at one site does not influence the binding reaction at another site: all binding interactions are independent of one other. These assumptions are necessary to make the model tractable, and are validated by comparing the model predictions to the available experimental data.
Here we used the shorthand notation similar to that of Segel  where num = k1k2k3k4k5k 6 K P K Q K R , num = k-1k-2k-3k-4k-5k-6K A K B K C , and Coef AB = K C K P K Q K R (k1k2k3k6 (k4 + k-4)), and so on. V f and V r have units of mass per unit time per mass of protein. Other kinetic parameters associated with the binding of reactants and products have the units of concentration (mass per unit volume).
reducing the total number of independent unknown kinetic parameter to thirteen.
Kinetic model of OGDHC using a ter-ter enzyme mechanism
We apply the above general form of the ter-ter enzyme mechanism for the analysis of available expe-rimental data on the kinetic of OGDHC to estimate the unknown kinetic parameters and to elucidate whether the proposed mechanism is able to explain the available kinetic data.
Here biochemical reactants, such as α KG, correspond to ensemble chemical species, such as α KG2-, H α KG-, etc. The chemical reaction in Equation (1) is unambiguously balanced in terms of mass and charge, whereas this biochemical reaction is not. In this reaction, the reactant CO2, tot represents the sum of aqueous carbon dioxide and bicarbonate species (, and H2CO3).
where is the standard Gibbs free energy of the reference reaction which is computed using the basic thermodynamic data (298.15 K, I = 0.15 M) listed in Li et al. .
The kinetic expressions for the regulatory effects of various cofactors on OGDHC (which depend on 10 additional kinetic parameters) are parameterized in the Results section, and are estimated separately (see below). Because CO2, tot dependent terms are not included in the denominator in Equation (11), CO2, tot dependent product inhibition is not accounted for in Equation (12).
Kinetic parameter values for 2-oxoglutarate dehydrogenase complex
Basic kinetic parameters
Kinetic parameters for regulatory cofactors
Vmax μmol mg-1 min-1a
K aCa (μM)
K mA (mM)
K iATP (mM)
K mB (μM)
K iADP (mM)
K mC (μM)
K aMg (μM)
K ia (mM)
K ir (mM)
K ic (mM)
K iq (μM)
K aH (pH)
where E* is the least square difference between model simulations and experimental data, and x i is the optimized value of the i th parameter.
Since all kinetic parameters in the model are measured relative to species concentration, we performed a composition analysis to estimate the concentration of all ionic species all experiments analyzed here .
Parameterization of basic kinetic mechanism of OGDHC
In this section, we present the detailed parameterization and validation of the proposed kinetic model based on the available experimental data on the kinetics of OGDHC, measured in a wide variety of experimental conditions. To study the catalytic mechanism of OGDHC, McMinn and Ottaway  investigated the kinetic properties of the OGDHC system, which was prepared from fresh pig heart mitochondria. Following the method of Fromm  in which it was concluded from initial velocity studies that the catalytic mechanism of OGDHC is not consistent with the Hexa-Uni-Ping-Pong mechanism. While their observations suggest a random order kinetic mechanism with respect to the binding of NAD- and CoASH4- and release of Succinyl-CoA4-, the binding of 2-oxoglutarate and release of CO2 is described as a Ping-Pong mechanism. Initial velocity kinetics measured by Smith et al.  with purified pig heart mitochondria OGDHC showed that Succinyl-CoA4- and NADH2- were inhibitors, but no inhibitory effects were observed with GTP or ATP. Their results also show that Succinyl-CoA4- inhibition was competitive with CoASH4- and independent of the NAD- oxidation-reduction state. These data are used here to identify the kinetic parameters of our OGDHC model.
Parameterization of the cofactor-dependent regulatory mechanisms
Denton et al.  conducted a number of experiments to study the effects of Ca2+, pH, and adenine nucleotides on the activity of OGDHC from pig heart mitochondria. Their data are used here to identify the kinetic parameters that characterize the activation/inhibition mechanism of Ca2+, pH, and adenine nucleotides (See Figure caption for details). Enzyme activities are expressed as units of enzyme activity per mg of protein. (One unit of activity is the amount of enzyme which transforms 1 μmol of substrate per minute at 30°C).
The modified flux expression for OGDHC is obtained by substituting Equation (14) into Equation (11). Three adjustable parameters (α Ca , β Ca , and K aCa ) are estimated based on the data from Denton et al. ; the model fits are shown in Figure 3A-B. For model simulations, the kinetic constants for substrates: 2-oxoglutarate, NAD-, and CoASH4- are fixed at their previously estimated values, obtained from data of McMinn and Ottaway . Based on the fits to these data we find that β Ca is close to one. For simplicity, we fixed β Ca = 1, meaning that Ca2+ affects only the K m of 2-oxoglurate, not the Vmax. This mechanism is also consistent with the conclusions of Denton et al. . (Statistical analysis of different model formulizations supports the validity of the null hypothesis that adding the extra adjustable parameters, β Ca , does not lead to any significant improvement in fitting results.) The estimated values of the kinetic parameters are summarized in Table 1.
Experimentally it has been shown that the maximal activity of OGDHC is largely unaffected by changes in pH over the range 6.6-7.4, whereas the Km of the enzyme is markedly altered by pH in this range . In our model, the effect of pH on the OGDHC activity was described based on the observations studies of Denton and colleagues . Here, protons are treated as the essential activators of OGDHC which increase the binding affinity of the enzyme to 2-oxoglutarate. Therefore, the K m of 2-oxoglurate is modified by multiplying the term K aH /[H+] such that KmA, 1= (KmA, 1K aH )/[H+] (in Equation (15)). Figure 3C illustrates the model fits to the data obtained from Table 2 in McCormack and Denton  where the activity of ODGHC was studied under varying pH in both presence and absence of Ca2+ in the buffer.
Mg2+ is known to regulate the activity of OGDHC. In a recent study, Rodriguez-Zavala et al.  examined the effects of ligands, such as ATP, ADP, Ca2+, and Mg2+ on the activity of OGDHC in both isolated pig heart enzyme complex and mitochondrial extracts. These data facilitate the characterization of the regulatory effect of Mg2+ on the OGDHC activity and are used here to estimate the Mg2+ associated kinetic parameters. Enzyme activity is measured in nmol NADH2+ formed per minute per mg protein.
The data shown in Figure 5 are used to estimate the adjustable kinetic parameters related to Mg2+ ions in our kinetic model for OGDHC. Fits to the data are plotted in Figure 5 and the parameter values summarized in Table 1. The developed model is able to satisfactorily explain the effect of Mg2+ ions on the enzyme activity. These results, combined with those shown in Figures 2 and 3, imply that the matrix free Ca2+ and Mg2+ ions concentrations exert significant and distinct effects on the OGDHC activity.
Complete flux expression for the 2-oxoglutarate dehydrogenase complex
Estimated values of K aH , K iATP , K aADP , K aMg , α ADP , α ATP , α Ca , α Mg , and β Mg are listed in Table 1.
Independent validation of the developed kinetic model of OGDHC
Finally, the model is independently validated (corroborated) by comparing the model predictions to the initial rate data of Rutter and Denton  on the kinetics of OGDHC obtained from permeabilized mitochondria and mitochondrial extracts (see Figure 6(A-C)). They studied the regulations of NAD-linked isocitrate dehydrogenase and 2-oxoglutarate dehydrogenase by Ca2+, nucleotide and nicotinamide nucleotides in permeabilized rat heart mitochondria and in mitochondria extracts. Data from their study were not used for estimation of model parameters and used here to further validate the proposed mechanisms and regulation of OGDHC. Therefore, the flux expression of Equation (17) was used for simulations with the values of the kinetic parameters the same as estimated before (see Table 1). The model accurately describes the kinetics and regulation of OGDHC, observed experimentally, without having to re-estimate the model kinetic parameters, signifying the accuracy of the model and the associated model parameters.
To determine the degree to which the model simulations are sensitive to the estimated parameter values, the relative sensitivities are computed and listed in Table 1. A high sensitivity value indicates that a small change in a given parameter can lead to significant changes in model outputs, used to identify the parameter values. All of our adjustable parameters of the model have sensitivities over 30%. Two parameter estimates (K aCa and K aMg ) show relatively low sensitivity compared to the others, indicating that predictions of the developed model are less sensitive to these two values. This implies that these two parameters may not be identified accurately by the present analysis, given the sparseness of the data sets analyzed in this work regarding the regulation of OGDHC by Ca2+ and Mg2+. Further experiments are required to adequately establish the appropriate regulatory mechanisms and the robustness of each model parameters.
A number of kinetic models have been previously developed to explain the basic catalytic mechanisms and regulations by cofactors of OGDHC. Sanadi et al.  first proposed a Hexa-Uni-Ping-Pong mechanism for the overall reaction by studying various roles and locations of the cofactors: thiamine pyrophosphate, lipoic acid, and FAD2- within this complex. Hemada et al.  conducted kinetic studies and proposed a similar mechanism to that of Sandi et al. and suggested that NADH2- is a competitive inhibitor of NAD-. Whereas Smith et al.  suggested a noncompetitive inhibition of NADH2- with NAD-, the catalytic mechanism was not consistent with Sandi et al. Later, McMinn and Ottaway  tested a series of possible alternate mechanisms using computer optimization techniques and initial velocity studies and concluded that the binding of NAD- and CoASH4- and the release of Succinyl-CoA4- is a random order, whereas the binding of the substrate 2-oxoglutarate and release of the product CO2 still follows a Ping-Pong mechanism. Besides above experimental studies, a number of integrated models of mitochondrial bioenergetics have been developed which used different type of OGDHC models. Cortassa et al.  describe the activity of OGDHC as a function of Ca2+, Mg2+ and substrate concentrations using phenomenological terms. Wu et al.  used a simple kinetic model of OGDHC from Kohn and Garfinkel  in their integrated model of TCA cycle that does not incorporate the regulatory effect of metal ion cofactors. In a recent integrated study of mitochondrial bioenergetics, Bazil et al.  developed a kinetic model of OGDHC based on a Hexa-Uni-Ping-Pong mechanism with a general description of the cofactor dependency of OGDHC activity. In summary, there have been a wide variety of kinetic models of OGDHC with contrasting kinetic mechanism and cofactor regulations.
In this paper, we developed a unified mechanistic model of OGDHC, in which Ca2+, Mg2+, ADP3-, and pH are treated as activators and ATP4- as inhibitor of the OGDHC activity. The present model offers more realistic and meaningful explanations on the catalytic properties and regulation mechanisms of OGDHC than previous attempts. The analysis also provides a unique set of kinetic parameters that consistently describe a wide variety of experimental data sets on OGDHC function, obtained from diverse sources. Based on the assumed ter-ter mechanism and associated conformational changes, we are able to consistently reproduce the observed kinetics of OGDHC with a minimal number of model parameters. Thus, the proposed mechanism is found to be more appropriate compared to other alternate kinetic models .
Effects of nucleotides on the OGDHC activity
Energy-linked regulators, ADP and ATP, as well as inorganic phosphate, have been investigated for over two decades for their profound effects on kinetic properties of OGDHC. Kinetic studies of mammalian OGDHC, isolated from varied sources, have shown that ADP causes activation of OGDHC [44–47]. This enzyme complex is sensitive to ADP, where ADP significantly decreases the K m for 2-oxogluterate without affecting the maximum rate (Vmax) of the reaction via allosteric interactions. For example, studies on the OGDHC of rat heart mitochondria show a seven fold decrease in the K m value for 2-oxogluterate by ADP and thereby strongly increases the affinity of OGDHC for the substrate . Other studies of OGDHC from human heart make similar conclusions on the activating effect of ADP . It has also been shown that, at subsaturating concentration of 2-oxoglutarate the relationship between initial reaction rates of OGDHC and concentration of ADP is sigmoidal, suggesting a positive cooperativity in binding of ADP to the enzyme complex . In contrast, at a suboptimal 2-oxoglutarate concentration, ATP was shown to inhibit OGDHC activity in pig heart and bovine kidney mitochondria [23, 45]. In addition direct inhibition effect, recent investigations have shown the possible indirect inhibition OGDHC by ATP because of chelation of divalent ions which activate OGDHC, such as Ca2+, Mg2+. The activating action of ADP and inhibiting action of ATP are in competitive opposition (Figure 3E). It is still unclear if these two effectors bind on the same site on the complex or not. Model analysis based on the available kinetic data cannot exclude either possibility. In our mechanistic model, the regulatory effect of ADP and ATP is incorporated by assuming different binding sites for ADP and ATP in the enzyme complex and the model satisfactorily describes the activating effect of ADP and inhibitory effect ATP observed in many experiments .
While Zavala et al. interpret the data of Figure 5A to indicate that MgADP- is the effective activator of OGDHC activity, our model analysis of the available data sets on Mg2+, ADP, and ATP dependent kinetics (Figures 4 and Figure 5B and 5C) reveals that, magnesium and ADP have independent parallel effects on the OGDHC activity, the most parsimonious explanation of the data. However, al-ternative, more complex, models cannot be ruled out.
Like ADP, Pi has also been shown to decrease the K m value for 2-oxogutarate, without affecting the Vmax of OGDHC reaction [44, 50]. In a recent report, the Pi activation showed biphasic behavior, with pH dependence . In the physiological concentrations range Pi exerts monophasic activation of OGDHC , which can be descried by Equation (15) with three extra parameters. Due to the lack of consistent kinetic data for Pi effects, we do not integrate a Pi dependent regulation mechanism in our current model. However, this energy linked effector may be physiologically important. The overall rate of oxidative phosphorylation is largely determined by phosphorylation potential [51, 52]. In cells when ATP utilization increases, the production of ADP and Pi increase. Therefore, activation of OGDHC by ADP and Pi may represent a compensating effect.
Effect of pH on the OGDHC activity
Mitochondrial matrix proton (H+ ion) concentration is known to affect the OGDHC properties. Specifically, studies on pig heart OGDHC showed that the change in pH in the range 6.6 to 7.4 can significantly alter K m of the enzyme for 2-oxogluterate, without affecting its maximal activity. McCormack and Denton illustrated the effect of pH on OGDHC activity both in the absence and presence of Ca2+ in their assay mediums . Other experimental observations have shown that hydrogen ions favor the higher affinity of OGDHC for 2-oxoglutarate [45, 53]. In our model we hypothesize that hydrogen ions are essential activators of OGDHC activity to describe the observed pH dependency of the OGDHC kinetics.
Effects of Ca2+, Mg2+, and EGTA on the OGDHC activity
Studies by McCormack and Denton demonstrate the activating effects of Ca2+ ions on intra-mitochondrial dehydrogenases: pyruvate (PDH), NAD-isocitrate (NAD-ICDH), and 2-oxoglutarate (OGDHC) [24, 54, 55]. Specifically, the rise in cytosolic Ca2+ concentration in response to extrinsic stimuli, such as hormones can enhance mitochondrial oxidative metabolism via direct activation of these three Ca2+ sensitive dehydrogenases. Such mechanisms may serve as a complementary way to stimulate ATP-synthesis to meet the increased energy demand of the cell [24, 54, 55].
Mg2+ ion has also been shown to regulate the OGDHC activity either by directly activating the enzyme or by modulating the Ca2+ effects on the enzyme. In some studies, Mg2+ shows no effects on OGDHC activity . However, in other studies, Mg2+ is shown to increase the maximal activity of the enzyme and the affinity of OGDHC for 2-oxogluterate by enhancing the Ca2+ stimulatory effects on the enzyme complex [9, 25, 37]. These different observations could be accounted for the different levels of endogenous Ca2+ and Mg2+ present in the purified enzyme complex prepared by different methods. Another possible explanation is that the stimulatory effects of Mg2+ is TPP-dependent, which is not explicitly considered in our model. Panov and Scarpa  found that Mg2+ only exerts its stimulatory effects in the presence of TPP, though exclusion of TPP from the reaction medium has no effect on the initial enzyme activity in the absence of Mg2+. Also, it has been clearly shown that Mg2+ may affect the rate of oxidative phosphorylation in isolated mitochondria primarily via modulating the OGDHC activity . The site of action of Mg2+ ion on OGDHC is unknown. In the present model, we hypothesized a general scheme of nonessential activation of Ca2+, by considering two Ca2+ and Mg2+ binding sites on OGDHC. The Mg2+ effect is incorporated in our model by exclusively modifying the enzyme activity and 2-oxoglutarate binding step. So the parameters Vmax and K mA are accordingly expressed as functions of Mg2+ (Equation (16)). Currently, our model assumes that the turn-over rate of E1 is modified to same value for binding either one ion or two ions. And to make it simple, our model does not include possible interaction between Mg2+ and Ca2+ at high concentrations either. The effects are not additive  at high concentration, suggesting that Mg2+ and Ca2+ may compete for the binding site. Additional kinetic data set are necessary to test different mechanisms and refine our model to more accurately describe the nature of cation dependent kinetic of OGDHC.
EGTA, which is used in many studies to control Ca2+ ion concentration in reaction media, has been shown in experiments and theoretical analysis to inhibit the NAD-linked isocitrate dehydrogenases (ICDH) through the binding complex, MgEGTA . To date, it is still not clear if there is similar inhibition effect of EGTA or EDTA on the activity of OGDHC. McCormack and Denton  concluded that the OGDHC sensitivity toEGTA is very similar to that observed with ICDH , because addition of calcium chelators EGTA or EDTA is associated with a marked decrease in the activity of OGDHC at 0.2 mM 2-oxoglutarate. Panov and Scarpa , in ascribing the inhibition effect of EGTA to the complex formation between Ca2+ and chelators, concluded that the effected of Ca2+ and chelators is associated with different endogenous cation levels in different preparations. But this explanation cannot account for McCormack and Denton's observation that EGTA or EDTA causes a 40% decrease of activity of OGDHC after using Chelex remove much of the endogenous Ca2+ in the buffer. It is also noted that the K M for 2-oxoglutarate in the absence of Ca2+ is 4 ± 1.1 mM measured by Panov and Scarpa  for commercially available enzyme (Sigma, St. Louis, lot 44H80801), which is almost 15 times the estimated value used in our model (Table 1) and that found by previous workers [19, 27]. Only by using the reported K M of Panov and Scarpa , can we reproduce their data using the same mechanism (Equation (17)). In the absence of clear experimental evidence and sufficient data set, our model does not explicitly account for an inhibition effect of EGTA or EDTA.
Our mechanistic OGDHC model based on a detailed catalytic mechanism successfully provides a single consistent theoretical explanation for many previously unresolved experimental observations on the kinetics and regulations of OGDHC. In particular, it suggests the most plausible physiologically regulations of OGDHC by NAD(H) oxidation-reduction state, the nucleotide phosphorylation potential, pH and various metal ions (Mg2+ and Ca2+). As a rise in NADH can reduce the OGDHC flux and thereby provides feedback regulation through the electron transport chain, it is important to ask how NAD oxidation-reduction state and oxidative phosphorylation state exert a coherent regulation of OGDHC in physiological context. Furthermore, how does the OGDHC respond to stimuli via the mitochondrial Ca2+ transport system? Such questions may be addressed by applying the present model in an integrated framework  along with other dehydrogenases , the oxidative phosphorylation system , electron transfer system , and cation transport systems [59–61].
The authors thank Doug Disabato for participating in the initial part of the model development and gathering available experimental data for the analysis, Fan Wu, Kalyan Vinnakota, and Jason Bazil for helpful discussion. This work was supported by NIH grants R01-HL072011 and R01-HL095122.
- Tretter L, Adam-Vizi V: Alpha-ketoglutarate dehydrogenase: a target and generator of oxidative stress. Philos Trans R Soc Lond B Biol Sci. 2005, 360 (1464): 2335-2345. 10.1098/rstb.2005.1764.PubMedPubMed CentralGoogle Scholar
- Perham RN: Domains, motifs, and linkers in 2-oxo acid dehydrogenase multienzyme complexes: a paradigm in the design of a multifunctional protein. Biochemistry. 1991, 30 (35): 8501-8512. 10.1021/bi00099a001.PubMedGoogle Scholar
- Strumilo S: Often ignored facts about the control of the 2-oxoglutarate dehydrogenase complex. Biochemistry and Molecular Biology Education. 2005, 33 (4): 284-287. 10.1002/bmb.2005.49403304284.Google Scholar
- Sanadi DR, Littlefield JW, Bock RM: Studies on alpha-ketoglutaric oxidase. II. Purification and properties. J Biol Chem. 1952, 197 (2): 851-862.PubMedGoogle Scholar
- Denton RM, McCormack JG, Edgell NJ: Role of calcium ions in the regulation of intramitochondrial metabolism. Effects of Na+, Mg2+ and ruthenium red on the Ca2+-stimulated oxidation of oxoglutarate and on pyruvate dehydrogenase activity in intact rat heart mitochondria. Biochem J. 1980, 190 (1): 107-117.PubMedPubMed CentralGoogle Scholar
- Hansford RG, Castro F: Role of Ca2+ in pyruvate dehydrogenase interconversion in brain mitochondria and synaptosomes. Biochem J. 1985, 227 (1): 129-136.PubMedPubMed CentralGoogle Scholar
- McCormack JG, Bromidge ES, Dawes NJ: Characterization of the effects of Ca2+ on the intramitochondrial Ca2+-sensitive dehydrogenases within intact rat-kidney mitochondria. Biochim Biophys Acta. 1988, 934 (3): 282-292. 10.1016/0005-2728(88)90088-6.PubMedGoogle Scholar
- McCormack JG, Denton RM: The role of Ca2+ ions in the regulation of intramitochondrial metabolism and energy production in rat heart. Mol Cell Biochem. 1989, 89 (2): 121-125.PubMedGoogle Scholar
- Panov A, Scarpa A: Independent modulation of the activity of alpha-ketoglutarate dehy-drogenase complex by Ca2+ and Mg2+. Biochemistry. 1996, 35 (2): 427-432. 10.1021/bi952101t.PubMedGoogle Scholar
- Smith CM, Bryla J, Williamson JR: Regulation of mitochondrial alpha-ketoglutarate metabolism by product inhibition at alpha-ketoglutarate dehydrogenase. J Biol Chem. 1974, 249 (5): 1497-1505.PubMedGoogle Scholar
- Starkov AA, Fiskum G, Chinopoulos C, Lorenzo BJ, Browne SE, Patel MS, Beal MF: Mito-chondrial alpha-ketoglutarate dehydrogenase complex generates reactive oxygen species. J Neurosci. 2004, 24 (36): 7779-7788. 10.1523/JNEUROSCI.1899-04.2004.PubMedGoogle Scholar
- Sanadi DR: Enzymes. 1963, 307-344. 2, 7Google Scholar
- Massey V: The composition of the ketoglutarate dehydrogenase complex. Biochim Biophys Acta. 1960, 38: 447-460.PubMedGoogle Scholar
- Cleland WW: The kinetics of enzyme-catalyzed reactions with two or more substrates or products. I. Nomenclature and rate equations. Biochim Biophys Acta. 1963, 67: 104-137.PubMedGoogle Scholar
- Koike M, Reed LJ, Carroll WR: alpha-Keto acid dehydrogenation complexes. IV. Resolution and reconstitution of the Escherichia coli pyruvate dehydrogenation complex. J Biol Chem. 1963, 238: 30-39.PubMedGoogle Scholar
- Angelides KJ, Hammes GG: Mechanism of action of the pyruvate dehydrogenase multienzyme complex from Escherichia coli. Proc Natl Acad Sci USA. 1978, 75 (10): 4877-4880. 10.1073/pnas.75.10.4877.PubMedPubMed CentralGoogle Scholar
- Waskiewicz DE, Hammes GG: Fluorescence polarization study of the alpha-ketoglutarate dehydrogenase complex from Escherichia coli. Biochemistry. 1982, 21 (25): 6489-6496. 10.1021/bi00268a026.PubMedGoogle Scholar
- Angelides KJ, Hammes GG: Structural and mechanistic studies of the alpha-ketoglutarate dehydrogenase multienzyme complex from Escherichia coli. Biochemistry. 1979, 18 (25): 5531-5537. 10.1021/bi00592a001.PubMedGoogle Scholar
- Hamada M, Koike K, Nakaula Y, Hiraoka T, Koike M: A kinetic study of the alpha-keto acid dehydrogenase complexes from pig heart mitochondria. J Biochem. 1975, 77 (5): 1047-1056.PubMedGoogle Scholar
- McMinn CL, Ottaway JH: Studies on the mechanism and kinetics of the 2-oxoglutarate dehydrogenase system from pig heart. Biochem J. 1977, 161 (3): 569-581.PubMedPubMed CentralGoogle Scholar
- Fromm HJ: The use of competitive inhibitors in studying the mechanism of action of some enzyme systems utilizing three substrates. Biochim Biophys Acta. 1967, 139 (2): 221-230.PubMedGoogle Scholar
- Aevarsson A, Seger K, Turley S, Sokatch JR, Hol WG: Crystal structure of 2-oxoisovalerate and dehydrogenase and the architecture of 2-oxo acid dehydrogenase multienzyme complexes. Nat Struct Biol. 1999, 6 (8): 785-792. 10.1038/11563.PubMedGoogle Scholar
- McCormack JG, Denton RM: The effects of calcium ions and adenine nucleotides on the activity of pig heart 2-oxoglutarate dehydrogenase complex. Biochem J. 1979, 180 (3): 533-544.PubMedPubMed CentralGoogle Scholar
- McCormack JG, Denton RM: Role of Ca2+ ions in the regulation of intramitochondrial metabolism in rat heart. Evidence from studies with isolated mitochondria that adrenaline activates the pyruvate dehydrogenase and 2-oxoglutarate dehydrogenase complexes by increasing the intramitochondrial concentration of Ca2+. Biochem J. 1984, 218 (1): 235-247.PubMedPubMed CentralGoogle Scholar
- Rodriguez-Zavala JS, Moreno-Sanchez R: Modulation of oxidative phosphorylation by Mg2+ in rat heart mitochondria. J Biol Chem. 1998, 273 (14): 7850-7855. 10.1074/jbc.273.14.7850.PubMedGoogle Scholar
- Hayakawa T, Kanzaki T, Kitamura T, Fukuyoshi Y, Sakurai Y, Koike K, Suematsu T, Koike M: Mammalian alpha-keto acid dehydrogenase complexes. V. Resolution and reconstitution studies of the pig heart pyruvate dehydrogenase complex. J Biol Chem. 1969, 244 (13): 3660-3670.PubMedGoogle Scholar
- Hirashima M, Hayakawa T, Koike M: Mammalian alpha-keto acid dehydrogenase complexes. II. An improved procedure for the preparation of 2-oxoglutarate dehydrogenase complex from pig heart muscle. J Biol Chem. 1967, 242 (5): 902-907.PubMedGoogle Scholar
- Patel MS: Inhibition by the branched-chain 2-oxo acids of the 2-oxoglutarate dehydrogenase complex in developing rat and human brain. Biochem J. 1974, 144 (1): 91-97.PubMedPubMed CentralGoogle Scholar
- Lai JC, Cooper AJ: Brain alpha-ketoglutarate dehydrogenase complex: kinetic properties, regional distribution, and effects of inhibitors. J Neurochem. 1986, 47 (5): 1376-1386. 10.1111/j.1471-4159.1986.tb00768.x.PubMedGoogle Scholar
- Shylaja N, Maehara M, Watanabe K: Measurement of alpha-ketoglutarate dehydrogenase activity in tissue extracts and human platelets using reversed-phase high-performance liquid chromatography. Anal Biochem. 1990, 191 (2): 223-227. 10.1016/0003-2697(90)90211-Q.PubMedGoogle Scholar
- Northrop DB: Transcarboxylase. VI. Kinetic analysis of the reaction mechanism. J Biol Chem. 1969, 244 (21): 5808-5819.PubMedGoogle Scholar
- Qi F, Dash RK, Han Y, Beard DA: Generating rate equations for complex enzyme systems by a computer-assisted systematic method. BMC Bioinformatics. 2009, 10: 238-10.1186/1471-2105-10-238.PubMedPubMed CentralGoogle Scholar
- Segel IH: Enzyme kinetics: behavior and analysis of rapid equilibrium and steady-state enzyme systems. 1975, New York: WileyGoogle Scholar
- Li X, Dash RK, Pradhan RK, Qi F, Thompson M, Vinnakota KC, Wu F, Yang F, Beard DA: A database of thermodynamic quantities for the reactions of glycolysis and the tricarboxylic acid cycle. J Phys Chem B. 2010, 114 (49): 16068-16082. 10.1021/jp911381p.PubMedPubMed CentralGoogle Scholar
- Beard DA, Qian H: Relationship between thermodynamic driving force and one-way fluxes in reversible processes. PLoS ONE. 2007, 2 (1): e144-10.1371/journal.pone.0000144.PubMedPubMed CentralGoogle Scholar
- Qi F, Chen X, Beard DA: Detailed kinetics and regulation of mammalian NAD-linked isocitrate dehydrogenase. Biochim Biophys Acta. 2008, 1784 (11): 1641-1651.PubMedPubMed CentralGoogle Scholar
- Rodriguez-Zavala JS, Pardo JP, Moreno-Sanchez R: Modulation of 2-oxoglutarate dehydrogenase complex by inorganic phosphate, Mg(2+), and other effectors. Arch Biochem Biophys. 2000, 379 (1): 78-84. 10.1006/abbi.2000.1856.PubMedGoogle Scholar
- Rutter GA, Denton RM: Regulation of NAD+-linked isocitrate dehydrogenase and 2-oxoglutarate dehydrogenase by Ca2+ ions within toluene-permeabilized rat heart mitochondria. Interactions with regulation by adenine nucleotides and NADH/NAD+ ratios. Biochem J. 1988, 252 (1): 181-189.PubMedPubMed CentralGoogle Scholar
- Cortassa S, Aon MA, Marban E, Winslow RL, O'Rourke B: An integrated model of cardiac mitochondrial energy metabolism and calcium dynamics. Biophys J. 2003, 84 (4): 2734-2755. 10.1016/S0006-3495(03)75079-6.PubMedPubMed CentralGoogle Scholar
- Wu F, Yang F, Vinnakota KC, Beard DA: Computer modeling of mitochondrial tricarboxylic acid cycle, oxidative phosphorylation, metabolite transport, and electrophysiology. J Biol Chem. 2007, 282 (34): 24525-24537. 10.1074/jbc.M701024200.PubMedGoogle Scholar
- Kohn MC, Garfinkel D: Computer simulation of metabolism in palmitate-perfused rat heart. II. Behavior of complete model. Ann Biomed Eng. 1983, 11 (6): 511-531. 10.1007/BF02364082.PubMedGoogle Scholar
- Bazil JN, Buzzard GT, Rundell AE: Modeling mitochondrial bioenergetics with integrated volume dynamics. PLoS Comput Biol. 2010, 6 (1): e1000632-10.1371/journal.pcbi.1000632.PubMedPubMed CentralGoogle Scholar
- Beard DA: Simulation of cellular biochemical system kinetics. Wiley Interdiscip Rev Syst Biol Med. 2010, 3 (2): 136-146.PubMedPubMed CentralGoogle Scholar
- Lawlis VB, Roche TE: Regulation of bovine kidney alpha-ketoglutarate dehydrogenase complex by calcium ion and adenine nucleotides. Effects on S0.5 for alpha-ketoglutarate. Biochemistry. 1981, 20 (9): 2512-2518. 10.1021/bi00512a023.PubMedGoogle Scholar
- Lawlis VB, Roche TE: Inhibition of bovine kidney alpha-ketoglutarate dehydrogenase complex by reduced nicotinamide adenine dinucleotide in the presence or absence of calcium ion and effect of adenosine 5'-diphosphate on reduced nicotinamide adenine dinucleotide inhibition. Biochemistry. 1981, 20 (9): 2519-2524. 10.1021/bi00512a024.PubMedGoogle Scholar
- McCormack JG, Denton RM: A comparative study of the regulation of Ca2+ of the activities of the 2-oxoglutarate dehydrogenase complex and NAD+-isocitrate dehydrogenase from a variety of sources. Biochem J. 1981, 196 (2): 619-624.PubMedPubMed CentralGoogle Scholar
- Nichols BJ, Rigoulet M, Denton RM: Comparison of the effects of Ca2+, adenine nucleo-tides and pH on the kinetic properties of mitochondrial NAD(+)-isocitrate dehydrogenase and oxoglutarate dehydrogenase from the yeast Saccharomyces cerevisiae and rat heart. Biochem J. 1994, 303 (Pt 2): 461-465.PubMedPubMed CentralGoogle Scholar
- Ostrovtsova SA, Strumilo SA: Participation of adenosine diphosphate in regulation of the 2-oxoglutarate dehydrogenase complex from human heart. Biomed Biochim Acta. 1990, 49 (6): 515-517.PubMedGoogle Scholar
- Strumilo SA, Taranda NI, Vinogradov VV: Peculiarities of the regulation of adrenal oxoglutarate dehydrogenase complex by NADH and adenosine diphosphate. Biokhimiia. 1982, 47 (5): 724-732.PubMedGoogle Scholar
- Strumilo SA, Taranda NI, Vinogradov VV: Role of phosphate and divalent metal ions in regulation of the activity of the alpha-ketoglutarate dehydrogenase complex from adrenal cortex. Biokhimiia. 1981, 46 (1): 156-161.PubMedGoogle Scholar
- Beard DA, Kushmerick MJ: Strong inference for systems biology. PLoS Comput Biol. 2009, 5 (8): e1000459-10.1371/journal.pcbi.1000459.PubMedPubMed CentralGoogle Scholar
- Wu F, Zhang EY, Zhang J, Bache RJ, Beard DA: Phosphate metabolite concentrations and ATP hydrolysis potential in normal and ischaemic hearts. J Physiol. 2008, 586 (17): 4193-4208. 10.1113/jphysiol.2008.154732.PubMedPubMed CentralGoogle Scholar
- Pawelczyk T, Angielski S: Cooperation of Ca2+ and pH in regulation of the activity of the 2-oxoglutarate dehydrogenase complex and its components from bovine kidney cortex. Acta Biochim Pol. 1984, 31 (3): 289-305.PubMedGoogle Scholar
- McCormack JG, Denton RM: Influence of calcium ions on mammalian intramitochondrial dehydrogenases. Methods Enzymol. 1989, 174: 95-118.PubMedGoogle Scholar
- Denton RM, McCormack JG: Ca2+ transport by mammalian mitochondria and its role in hormone action. Am J Physiol. 1985, 249 (6 Pt 1): E543-554.PubMedGoogle Scholar
- Denton RM, Richards DA, Chin JG: Calcium ions and the regulation of NAD+-linked isocitrate dehydrogenase from the mitochondria of rat heart and other tissues. The Biochemical journal. 1978, 176 (3): 899-906.PubMedPubMed CentralGoogle Scholar
- Beard DA: A Biophysical Model of the Mitochondrial Respiratory System and Oxidative Phosphorylation. PLoS Comput Biol. 2005, 1 (4): e36-10.1371/journal.pcbi.0010036.PubMedPubMed CentralGoogle Scholar
- Chen X, Qi F, Dash RK, Beard DA: Kinetics and regulation of mammalian NADH-ubiquinone oxidoreductase (Complex I). Biophys J. 2010, 99 (5): 1426-1436. 10.1016/j.bpj.2010.06.063.PubMedPubMed CentralGoogle Scholar
- Dash RK, Qi F, Beard DA: A biophysically based mathematical model for the kinetics of mitochondrial calcium uniporter. Biophys J. 2009, 96 (4): 1318-1332. 10.1016/j.bpj.2008.11.005.PubMedPubMed CentralGoogle Scholar
- Pradhan RK, Beard DA, Dash RK: A biophysically based mathematical model for the kinetics of mitochondrial Na+-Ca2+ antiporter. Biophys J. 2010, 98 (2): 218-230. 10.1016/j.bpj.2009.10.005.PubMedPubMed CentralGoogle Scholar
- Pradhan RK, Qi F, Beard DA, Dash RK: Characterization of membrane potential dependency of mitochondrial Ca2+ uptake by an improved biophysical model of mitochondrial Ca2+ uniporter. PLoS One. 2010, 5 (10): e13278-10.1371/journal.pone.0013278.PubMedPubMed CentralGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.