Dynamic assembly of malate dehydrogenase–citrate synthase multienzyme complex in the mitochondria

  1. Joy Omini
  2. Inga Krassovskaya
  3. Taiwo Adeolu Dele-Osibanjo
  4. Connor Pedersen
  5. Toshihiro Obata  Is a corresponding author
  1. Department of Biochemistry and Center for Plant Science Innovation, University of Nebraska-Lincoln, United States

eLife Assessment

This important study provides novel information on multi-enzyme complexes, known as metabolons, that form between sequential enzymes in a metabolic pathway. Using an innovative NanoBiT split-luciferase system, the authors present compelling evidence that malate dehydrogenase (MDH1) and citrate synthase (CIT1) dynamically associate under different metabolic conditions in Saccharomyces cerevisiae. The findings suggest the dynamic MDH1-CIT1 interaction facilitates control of TCA pathway flux rate.

https://doi.org/10.7554/eLife.107953.3.sa0

Abstract

The tricarboxylic acid (TCA) cycle enzymes malate dehydrogenase (MDH1) and citrate synthase (CIT1) form a multienzyme complex, referred to as a metabolon, that channels intermediate oxaloacetate between their reaction centers. Given that the MDH1–CIT1 metabolon enhances pathway reactions in vitro, its dynamic assembly is hypothesized to contribute to TCA cycle regulation in response to cellular metabolic demands. Here, we demonstrated that yeast mitochondrial MDH1 and CIT1 dissociated when aerobic respiration was suppressed by the Crabtree effect and associated when the respiratory activity was enhanced by acetate. Pharmacological TCA cycle inhibition dissociated the complex, whereas electron transport chain inhibition enhanced the interaction. The multienzyme complex assembly was related to the mitochondrial matrix acidification and oxidation, as well as cellular levels of malate, fumarate, and citrate. These factors significantly affected the MDH1–CIT1 complex affinity in vitro. Especially, variations in buffer pH within the physiological pH range between 6.0 and 7.0 in the mitochondrial matrix significantly impacted the MDH1–CIT1 affinity. These results demonstrate the dynamic association and dissociation of the MDH1–CIT1 metabolon and its relationship with respiratory activity, supporting metabolon dynamics as an integral factor in metabolic regulation governed by multiple factors such as mitochondrial pH and metabolite levels.

Introduction

Enzymes catalyzing sequential reactions can interact to form a multienzyme complex, often called a ‘metabolon’, which channels the metabolic intermediate within the complex. Metabolite channeling can mediate pathway reactions by concentrating the reaction intermediates near the enzyme active site and sequestering them from competing reactions (Obata, 2020; Spivey and Ovádi, 1999; Wheeldon et al., 2016). Thus, the dynamic assembly of multienzyme complex is believed to quickly regulate cellular metabolic flux by changing their degree of association and dissociation without involving time-consuming and resource-demanding protein synthesis, degradation, and modification (Halper and Srere, 1977; Morgunov and Srere, 1998; Obata, 2019; Omini et al., 2024; Srere, 1985; Sumegi et al., 1993). However, limited experimental evidence shows metabolon dynamics and the mechanisms regulating its association.

The tricarboxylic acid (TCA) cycle multienzyme complex composed of malate dehydrogenase (MDH) and citrate synthase (CS) is conserved in all organisms analyzed so far, including animals, bacteria, yeast, and plants (Meyer et al., 2011; Morgunov and Srere, 1998; Robinson et al., 1987; Zhang et al., 2017). MDH and CS catalyze the key steps in respiratory metabolism as their reactions comprise the carbon entry steps from the glycolytic processes as malate and acetyl-CoA, respectively (Sweetlove et al., 2010). MDH–CS complex is considered a dynamic protein complex (Bulutoglu et al., 2016; Omini et al., 2021) that protects the channeled intermediate, oxaloacetate, from degradation and competing pathways in the bulk phase of the cell (Bulutoglu et al., 2016). Importantly, MDH–CS complex formation and oxaloacetate channeling are considered essential for the forward TCA cycle flux to occur since MDH forward reaction to synthesize oxaloacetate is thermodynamically unfavorable in physiological conditions (Fukuda et al., 2008). These findings suggest that the dynamic association and dissociation of the MDH–CS multienzyme complex play a role in coordinating the TCA cycle and associated metabolic pathways. To address this hypothesis, this study focuses on demonstrating the in vivo dynamics of the MDH–CS complex in relation to respiratory states, a crucial aspect of achieving metabolic regulation. Another key requirement, the direct impacts of the metabolon on pathway flux, will be addressed in a separate study.

This study utilized the yeast mitochondrial MDH (MDH1)–CS (CIT1) multienzyme complex as the model. Budding yeast, Saccharomyces cerevisiae, dramatically rearranges its central carbon metabolism in response to changes in nutrient availability and stress conditions (Causton et al., 2001; Rep et al., 1999; Rolland et al., 2002). Remarkably, yeast respiratory metabolism, involving aerobic respiration and fermentation, is highly adaptable to substrate availability (Pfeiffer and Morley, 2014). When oxygen and non-fermentable respiratory substrates, such as acetate, are abundant, aerobic respiration is upregulated, with increased carbon flux through the TCA cycle (Gerstmeir et al., 2003; Lee et al., 2011). In the presence of fermentable sugars, such as glucose and fructose, the fermentation pathway is upregulated, and aerobic respiration is repressed, even when oxygen is available (Zhang et al., 2022). Furthermore, aerobic respiration and fermentation cooperate in the presence of a poorly fermentable carbon source like raffinose (Guaragnella et al., 2013). Application of fermentable substrates to respiring yeast cells induces a substantial shift from aerobic respiration to fermentation, known as the ‘Crabtree effect’, involving a massive transcriptional rearrangement of enzyme genes related to aerobic respiration and fermentation (Gancedo, 1998; Ronne, 1995; Thierie and Penninckx, 2010). This inducible metabolic shift in yeast makes it an ideal model for investigating the dynamic relationship between the MDH1–CIT1 multienzyme complex and TCA cycle flux. Understanding the mechanisms of the Crabtree effect is crucial for metabolic engineering applications to enhance the supply of TCA cycle intermediates for desired product synthesis (Yin et al., 2015) and for gaining insights into metabolic regulation in other eukaryotic systems, including cancer cells, which exhibit metabolic shifts similar to the Crabtree effect (Diaz-Ruiz et al., 2011).

Various allosteric regulators and environmental factors, including pH and redox state, affect MDH and CS enzyme activities (Beeckmans, 1984; Williamson and Cooper, 1980). These factors alter protein conformations, influencing the MDH–CS complex affinities. Previous in vitro studies have demonstrated that NAD+, malate, succinate, acetyl-CoA, α-ketoglutarate, and acidic pH enhance the MDH–CS interaction, while NADH, citrate, and basic pH weaken it (Omini et al., 2024; Omini et al., 2021; Tompa et al., 1987; Wu and Minteer, 2015). These allosteric regulators promote the MDH–CS interaction when the TCA cycle substrates are abundant and products are limited, indicating the role of the MDH–CS multienzyme complex in the TCA cycle feedback regulation. The electron transport chain (ETC) activity affects the mitochondrial matrix pH, redox state, and ATP content and is closely related to cellular respiratory flux distributions (Santo-Domingo and Demaurex, 2012). Therefore, the MDH–CS multienzyme complex may associate and dissociate according to the microenvironment in the mitochondrial matrix, such as metabolite concentrations, pH, redox state, and energy levels.

We adopted a NanoBiT protein–protein interaction assay system (Dixon et al., 2016) to monitor real-time multienzyme complex interactions in living yeast cells using continuous microplate readings. The NanoBiT split-luciferase system is based on the 18 kDa NanoLUC luciferase derived from deep-sea shrimp. The small size of this reporter is designed to minimize structural and behavioral interference with the target proteins, while offering a dynamic range and brightness approximately 150-fold greater than those of conventional luciferases, enabling detection of weak, transient interactions. The NanoBiT subunits exhibit low intrinsic affinity (Kd = 190 µM) and rapid association and dissociation kinetics (kon = 500 M–1 s–1, koff = 0.2 s–1), ensuring that the luminescent complex formation is dictated by the interaction of the fused proteins rather than the affinity of the tags themselves (Dixon et al., 2016). While other fluorescence-based methods, such as FLIM-FRET and BRET, are widely used for in vivo protein–protein interaction studies, they are limited in yeast microplate assays due to high cellular autofluorescence. The NanoBiT system overcomes these constraints, providing semi-quantitative monitoring of relative fluctuations in protein interaction levels over time across a diverse range of metabolic conditions.

In this study, the NanoBiT system examined the MDH1–CIT1 multienzyme complex association in the presence of various carbon sources and respiratory inhibitors. The results showed the relationship between respiratory activity and the dynamic MDH1–CIT1 complex assembly. Analysis of the mitochondrial matrix microenvironment using fluorescent biosensors and cellular metabolite profiles revealed crosstalk among cellular respiratory status, the MDH1–CIT1 interaction, and the mitochondrial matrix microenvironment, including pH. These findings indicate the functional dynamics of the MDH–CS multienzyme complex and how they coordinate with the metabolic states of the TCA cycle and adjacent pathways.

Results

MDH1–CIT1 interaction was detected exclusively in respiratory conditions

We employed a NanoBiT split-luciferase system (Dixon et al., 2016) for semi-quantitative relative detection of the MDH1–CIT1 multienzyme complex association in living yeast cells. The codon-optimized sequences of NanoBiT subunits were integrated into the direct downstream of the MDH1 and CIT1 exons in the yeast genome by gene homologous recombination to fuse the NanoBiT subunits to the native enzymes. The NanoBiT subunits reconstitute the NanoLUC holoenzyme to produce luciferase signals when MDH1 and CIT1 are in proximity. The affinity between the NanoBiT subunits is very low and does not stabilize the interactions of the tagged proteins (Dixon et al., 2016). The integration of NanoBiT subunits had no significant effect on the cellular growth rate or cellular MDH and CS enzyme activities (Figure 1—figure supplement 1). Although the MDH1–CIT1 reporter line did not show a detectable luciferase signal when it was grown in a fermentation condition (glucose-containing SD media; SD-Gluc), it showed a substantial luciferase signal in mixed-respiration (raffinose-containing SD media; SD-Raff) and respiration (acetate-containing SD media; SD-Acet) conditions (Figure 1A). These results indicate that MDH1 and CIT1 interact under respiratory conditions. The MDH1–CIT1 interaction can be related to the oxidative respiration and the TCA cycle activity since oxygen consumption was detected only in SD-Raff and SD-Acet media but not in the SD-Gluc medium (Figure 1B), as previously reported (Klein et al., 1998; Polakis and Bartley, 1965; Rep et al., 1999; Rolland et al., 2002; Strittmatter, 1957). However, MDH1 and CIT1 expressions were downregulated in SD-Gluc growing cells near the detection limit (Figure 1C), making it difficult to conclude if the undetectable MDH1–CIT1 interaction was due to complex dissociation or low protein levels.

Figure 1 with 1 supplement see all
MDH1–CIT1 interaction under respiration, fermentation, and mixed-respiration conditions.

Yeast cells were grown in the minimum media containing acetate (SD-Acet), glucose (SD-Gluc), and raffinose (SD-Raff) to the exponential growth phase. (A) Luciferase signal indicating MDH1–CIT1 complex interaction (N = 20). (B) Cellular oxygen consumption rate (N = 3). (C) MDH1 and CIT1 protein levels detected by western blotting. Numbers on the left indicate the position of the molecular weight markers. DLD1 is a loading control of mitochondrial protein. . In A and B, data are presented as mean ± SEM, and the differences between conditions were tested by Student’s t-test. Asterisks indicate significant differences (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant).

Figure 1—source data 1

PDF file containing the original western blotting images showing MDH1 and CIT1 protein abundance in Figure 1C.

https://cdn.elifesciences.org/articles/107953/elife-107953-fig1-data1-v1.pdf
Figure 1—source data 2

Original files for western blot analysis displayed in Figure 1C.

https://cdn.elifesciences.org/articles/107953/elife-107953-fig1-data2-v1.zip

Crabtree induction reduced the MDH1–CIT1 interaction

To test the relationship between respiratory activity and the MDH1–CIT1 complex association when MDH1 and CIT1 enzymes are abundant in the cell, we monitored the time course of the MDH1–CIT1 interaction after a rapid shift from aerobic to anaerobic respiration via the Crabtree effect. The Crabtree effect was induced by applying 2% glucose to the SD-Raff-grown MDH1–CIT1 NanoBiT reporter line (Figure 2, Figure 2—figure supplement 1). The luciferase signals only slightly changed for 30 min upon glucose application, followed by a steep decline. In contrast, the control cells retained the initial signal level for 100 min (Figure 2A, Figure 2—figure supplement 1D). Decrease in the MDH1–CIT1 interaction was validated by a co-immunoprecipitation assay at 1.5 and 2.5 hr after the glucose application (Figure 2—figure supplement 1A). Other Crabtree inducers, such as fructose and sucrose, also reduced the MDH1–CIT1 interaction, while the addition of non-fermentable sugars, including galactose, caused no significant change in MDH1–CIT1 interaction (Figure 2—figure supplement 2A–D, Figure 2—figure supplement 3), likely due to raffinose in the media suppressing its effects. The signal decline following glucose application was partially reversed by co-application with a fermentation inhibitor, 100 mM phosphate (van Urk et al., 1989; Figure 2A). Phosphate co-application slowed the decrease of luciferase signal compared to the glucose-applied cells, and the signal was not statistically significantly lower than the control until 80 min after application. The oxygen consumption rate significantly decreased when glucose was added to SD-Raff-grown cells but was slightly recovered by co-application of glucose with phosphate (Figure 2—figure supplement 2E). These results indicate that the respiratory suppression by the Crabtree effect is related to MDH1–CIT1 complex disruption.

Figure 2 with 4 supplements see all
MDH1–CIT1 complex association, mitochondrial microenvironments, and cellular metabolite levels during Crabtree effect induction.

Cells were cultured in fresh SD-Raff media in the control condition (black). The Crabtree effect was induced by the 2% glucose application to the SD-Raff-grown cells at 0 min (red). (A) NanoBiT signal indicating MDH1–CIT1 interaction. Relative luminescence unit (RLU) was calculated by normalizing the luciferase signals to the average signals during three pre-treatment time points. SD-Raff-grown cells were also co-treated with 2% glucose and a fermentation inhibitor, 100 mM phosphate, at 0 min (gray). (B) Mitochondrial matrix pH. (C) Mitochondrial matrix redox states as GSH/GSSG equivalent (mV). (D) Mitochondrial matrix ATP level indicated by the ratio between 560 and 510 nm emission signals. All data in A–D are presented as mean ± s.d. (N = 4). (E) Cellular metabolite levels at 80 min. The boxes, lines, error bars, and points indicate interquartile range, median, minimum, and maximum values, and outliers, respectively (N = 7). Statistical differences against the control samples were assessed using the Student’s t-test at each time point. Asterisks indicate significant differences (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant).

MDH1 protein level monitored with the yeast strain expressing MDH1 fused with full-length NanoLUC luciferase showed a 20% decline within 100 min upon glucose addition (Figure 2, Figure 2—figure supplement 1B). Western blotting also showed a slight decrease in the MDH1 and CIT1 levels (Figure 2—figure supplement 1C). The time course of the MDH1–CIT1 NanoBiT signal was inconsistent with that of MDH1-NanoLUC, and the 20% decline in protein content cannot fully explain the decrease of MDH1–CIT1 signal by over 50% within 100 min after glucose addition (Figure 2A, Figure 2—figure supplement 1B). The reduction of MDH1–CIT1 interaction was not accompanied by the MDH1 and CIT1 abundance in a repeated experiment (Figure 2—figure supplement 1D, E). The interaction index, calculated by normalizing the NanoBiT signal to MDH1 and CIT1 abundances measured with MDH1-NanoLUC and CIT1-NanoLUC (see Materials and methods), showed reduced interaction after 60 min of glucose application. Therefore, the decline in the MDH1–CIT1 NanoBiT signal reflects the MDH1–CIT1 complex dissociation co-occurring with the Crabtree effect, even though we cannot exclude a partial effect of protein degradation (Litsios et al., 2018).

Oxidative respiration significantly influences the microenvironments in the mitochondrial matrix (Beeckmans, 1984; Williamson and Cooper, 1980). To assess if the MDH1–CIT1 interaction decline is related to the changes in the microenvironments of the mitochondrial matrix by the Crabtree effect, we monitored the mitochondrial matrix pH, ATP concentration, and redox state by expressing fluorescent biomarkers, pHluorin, mito-GoAteam2, and mito-roGFP1, respectively, in the mitochondrial matrix of the MDH1–CIT1 reporter line (Figure 2—figure supplement 4A–D). Mitochondrial matrix pH in the SD-Raff-grown cells was 7.2 and temporally declined to 6.8 in the first 25 min of glucose application. The control cell matrix pH gradually decreased to 6.2 and remained lower than in glucose-applied cells (Figure 2B). The alterations of the MDH1–CIT1 interaction and mitochondrial matrix parameters in the control condition are not due to carbon starvation, since 2% raffinose application had no significant effect (Figure 2—figure supplement 4E). The redox state stayed around –288.0 mV in both the control and glucose cells (Figure 2C). The ATP level slightly increased during the experimental period of 100 min in both the control and glucose cells, with slightly but significantly lower levels in the glucose than in the control cells from 55 to 100 min (Figure 2D). Thus, mitochondrial matrix pH and ATP levels were maintained after glucose was applied to cells, while they decreased and increased in control cells during the experimental period, indicating that mitochondrial pH and ATP levels are associated with MDH–CS complex dissociation during the Crabtree induction.

Our previous in vitro study showed that metabolite abundance influences the MDH–CS complex interaction (Omini et al., 2021). To assess the possible roles of metabolite accumulation in controlling MDH1–CIT1 complex interaction, we determined the cellular levels of 38 metabolites (Source data 1) by gas chromatography-mass spectrometry (GC-MS). We focus here on the intermediates of the TCA cycle and related pathways (Figure 2E). At 80 min after 2% glucose application to SD-Raff-grown cells, malate, citrate, fumarate, and glutamate levels decreased, while α-ketoglutarate content increased significantly compared to the control cells (Figure 2G). Thus, the glucose-induced shift from aerobic respiration to fermentation altered the intracellular metabolite profile in SD-Raff-grown cells, potentially influencing the MDH1–CIT1 affinity.

MDH1–CIT1 complex interacts in relation to the TCA cycle activity

We further assessed the relationship between the MDH1–CIT1 complex association and TCA cycle activity. Acetate supplies acetyl-CoA to the TCA cycle and enhances its activity (Cavero et al., 2003). The application of 1% sodium acetate gradually increased the MDH1–CIT1 interaction while MDH1 and CIT1 protein abundance slightly reduced (Figure 3A, Figure 3—figure supplement 1A–D). The TCA cycle activation by ethanol (Schüller, 2003) also enhanced the complex association (Figure 3—figure supplement 1E). However, it is accompanied by an increase in MDH1 and CIT1 abundance (Figure 3—figure supplement 1F, G), making the effects of ethanol on the complex association unclear (Figure 3—figure supplement 1H). On the other hand, TCA cycle inhibition by 5 mM sodium arsenite, which impedes α-ketoglutarate dehydrogenase (Lee et al., 2011), led to a decline in MDH1–CIT1 interaction within 10 min (Figure 3B, Figure 3—figure supplement 1I–L). Aminooxyacetate (AOA) is an aminotransferase inhibitor and reduces TCA cycle activity by blocking the malate-aspartate NADH shuttle, which intersects with the cycle (Borst, 2020; Eto et al., 1999; Molinié et al., 2022). The MDH1–CIT1 interaction gradually reduced following 0.5 mM AOA application (Figure 3C, Figure 3—figure supplement 1M–P). These results support that the MDH1–CIT1 multienzyme complex association is related to TCA cycle activity.

Figure 3 with 1 supplement see all
MDH1–CIT1 complex association, mitochondrial matrix microenvironments, and cellular metabolite levels following tricarboxylic acid (TCA) cycle activation and inhibition.

Cells were cultured in SD-Raff media in the control condition (black). The TCA cycle activator (acetate, dark red) and inhibitors (arsenite, blue, and aminooxyacetate (AOA), pink) were applied at 0 min. (A–C) NanoBiT signal indicating MDH1–CIT1 interaction. Relative luminescence unit (RLU) was calculated by normalizing the luciferase signals by the average signals during three pre-treatment time points. (D–F) Mitochondrial matrix pH in control cells (black) and cells treated with acetate (dark red), arsenite (blue), and AOA (pink). (G–I) Mitochondrial matrix redox states as GSH/GSSG equivalent (mV). (D, H, L) Mitochondrial matrix ATP level indicated by the ratio between 560 and 510 nm emission signals of mito-GoATeam2 sensor. All data in A-L are presented as mean ± s.d. (N = 4). (M) Cellular metabolite levels after 80 min of treatment. The boxes, lines, error bars, and points indicate interquartile range, median, minimum, and maximum values, and outliers, respectively (N = 7). Statistical differences against the control samples were assessed using the Student’s t-test at each time point. Asterisks indicate significant differences (*p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001; ns, not significant).

Mitochondrial matrix pH declined to 6.1 immediately after acetate application (Figure 3D). Mitochondrial pH went slightly lower than that of the control cells after 40 min of TCA cycle inhibition by arsenite (Figure 3E), while AOA application had no effect on mitochondrial pH (Figure 3F). The mitochondrial matrix redox state was gradually increased in the control cells during the experimental period, while it was kept constant in the cells treated by acetate, arsenite, and AOA (Figure 3G–I), resulting in the significantly lower mitochondrial matrix redox state in the treated cells in the latter half of the experimental period. Mitochondrial ATP level was decreased by acetate treatment after 40 min (Figure 3J), although arsenite and AOA showed no effect on the mitochondrial ATP levels within the experimental period (Figure 3J, K). These changes in mitochondrial microenvironments did not show clear, direct relationships with the MDH1–CIT1 interaction.

We evaluated the cellular metabolite profile after treatment of SD-Raff grown cells for 80 min with 1% acetate, 5 mM arsenite, and 0.5 mM AOA to determine the relationship between change in metabolite levels and the MDH1–CIT1 complex interaction (Figure 3M). Acetate application significantly reduced malate, α-ketoglutarate, fumarate, succinate, and glutamate levels, while aspartate abundance significantly increased. Inhibition of the TCA cycle by arsenite significantly reduced the levels of all TCA cycle metabolites other than α-ketoglutarate. Aspartate and glutamate were also decreased, and glutamate was almost depleted in the cell (Figure 3M). AOA significantly decreased citrate, fumarate, and succinate levels while it increased α-ketoglutarate abundance (Figure 3M). Reduction of TCA cycle organic acids by the TCA cycle inhibition may be related to the reduction in the MDH1–CIT1 interaction.

ETC inhibition enhanced the MDH1–CIT1 interaction, coinciding with matrix acidification

Mitochondrial ETC activity is directly related to the mitochondrial microenvironments (Selivanov et al., 2008). To further assess the effects of mitochondrial matrix microenvironments on MDH1–CIT1 complex interaction, we treated SD-Raff-grown cells with ETC inhibitors. Inhibition of ETC complexes II, III, IV, and V with 20 mM malonate, 10 µM antimycin, 0.5 mM cyanide, and 1 mM oligomycin, respectively, significantly reduced oxygen consumption of SD-Raff-grown cells (Figure 4—figure supplement 1), showing their effects on respiration. MDH1–CIT1 interaction was enhanced by ETC inhibitors with minor alterations of MDH1 and CIT1 protein levels (Figure 4A–C, Figure 4—figure supplement 2), whereas the complex V (ATP synthase) inhibitor exerted no effect on MDH1–CIT1 interaction and mitochondrial matrix microenvironment in this experimental condition (Figure 4—figure supplement 3A–D). Complex II inhibition slightly and temporally enhanced the MDH1–CIT1 interaction, although the increase was not statistically significant compared to the control (Figure 4A, Figure 4—figure supplement 1B–E). The MDH1–CIT1 interaction was enhanced in the first 20 min of complex IV inhibition and reduced to the basal level afterward (Figure 4B, Figure 4—figure supplement 2E–H). Complex III inhibition by antimycin slowly increased MDH1–CIT1 interaction, and the signal increase peaked at 60 min of treatment (Figure 4C; Figure 4—figure supplement 2I–L).

Figure 4 with 4 supplements see all
MDH1–CIT1 complex association, mitochondria microenvironments, and cellular metabolite levels following mitochondrial electron transport chain (ETC) inhibition.

Cells were cultured in SD-Raff media in the control condition (black). The ETC inhibitors for complex II (malonate, purple, A–D), complex IV (cyanide, green, E–H), and complex III (antimycin, orange, I–L) were applied at 0 min. (A, E, I) NanoBiT signal indicating MDH1–CIT1 interaction. Relative luminescence unit (RLU) was calculated by normalizing the luciferase signals by the average signals during three pre-treatment time points. (B, F, J) Mitochondrial matrix pH. (C, G, K) Mitochondrial matrix redox states as GSH/GSSG equivalent (mV). (D, H, L) Mitochondrial matrix ATP level indicated by the ratio between 560 and 510 nm emission signals of mito-GoATeam2 sensor. All data in A-L are presented as mean ± s.d. (M) Cellular metabolite levels after 30 min for malonate and cyanide and after 80 min for antimycin treatment. The boxes, lines, error bars, and points indicate interquartile range, median, minimum, and maximum values, and outliers, respectively. Statistical differences against the control samples were assessed using the Student’s t-test at each time point. Asterisks indicate significant differences with p < 0.05.

Mitochondrial matrix pH was reduced by all ETC inhibitor treatments (Figure 4D–F), although the effect was very small in malonate-treated cells (Figure 4D). The matrix pH was transiently decreased and exhibited an inverse relationship with the MDH1–CIT1 interaction (Figure 4E). The mitochondrial matrix pH was reduced to 5.5 immediately after antimycin application and remained at the same level for 100 min (Figure 4F). These treatments show minor effects on the redox state of the mitochondrial matrix (Figure 4G–I). The matrix was gradually oxidized following antimycin treatment, but the difference relative to control cells was not statistically significant (Figure 4I). ETC inhibitions also showed relatively minor effects on mitochondrial ATP levels (Figure 4J–L), except for the reduction of ATP levels at the first 20 min of complex IV inhibition (Figure 4K), which is accompanied by the increase of MDH1–CIT1 interaction (Figure 4B). These results indicate the relationship between MDH1–CIT1 interaction and mitochondrial matrix acidification.

Since the results indicate a relationship between mitochondrial matrix pH and the MDH1–CIT1 interaction, we have tested the effect of an uncoupler, carbonyl cyanide 3-chlorophenylhydrazone (CCCP), on this interaction (Figure 4—figure supplement 4), which is expected to equilibrate the mitochondrial pH with the medium pH of 5.8 (Rosado et al., 2008). Treatment with 2 µM CCCP increased the MDH1–CIT1 interaction signal, suggesting that acidic pH enhances complex association (Figure 4—figure supplement 4A, D). However, this result should be interpreted with caution, as CCCP treatment compromised the luciferase signal from MDH1 and CIT1 fused to full-length NanoLUC (Figure 4—figure supplement 4B, C).

Cellular metabolite profiles were determined 30 min after malonate and cyanide treatments and 80 min after antimycin and oligomycin treatments. Complex II and V inhibition did not significantly alter cellular metabolite levels other than the malate and aspartate accumulation in malonate-treated cells (Figure 4M, Figure 4—figure supplement 3E). Complex IV inhibition by cyanide significantly reduced citrate and glutamate levels, while malate, α-ketoglutarate, fumarate, and succinate levels significantly increased (Figure 4M). Complex III inhibition significantly decreased citrate and α-ketoglutarate levels and increased the malate level (Figure 4M). These results support the relationship between the levels of TCA cycle intermediates and MDH1–CIT1 interaction.

MDH1–CIT1 interaction is affected by pH, ATP concentration, and metabolite availability in vitro

To evaluate the specific effects of microenvironments and metabolite availability on the MDH1–CIT1 interaction, we investigated the MDH1–CIT1 affinity in vitro using recombinant proteins. Acidic pH significantly enhanced the binding affinity of the MDH1–CIT1 complex (Figure 5A) within the physiological range observed in this study (Figures 2B, 3D). The Kd of the MDH1–CIT1 interaction was 3.48 µM at pH 7.2, while it was one magnitude lower at acidic pH (0.223 µM at pH 6.4 and 0.033 µM at pH 6.0). We also assessed the effects of the metabolites on the MDH1–CIT1 interaction (Figure 5B, Figure 5—figure supplement 1). We tested malate, fumarate, citrate, α-ketoglutarate, succinate, glutamate, and aspartate since their cellular levels changed when the MDH1–CIT1 interaction altered (Figures 2E, 3M). Citrate destroyed the MDH1–CIT1 interaction, while malate and fumarate significantly enhanced complex affinity (Figure 5B). α-ketoglutarate, succinate, and aspartate slightly enhanced MDH1–CIT1 complex affinity, although glutamate did not affect the interaction (Figure 5—figure supplement 1A). ATP also enhanced the MDH1–CIT1 complex affinity (Figure 5—figure supplement 1B) as observed in porcine enzymes (Omini et al., 2021). Thus, malate, fumarate, citrate, ATP levels, and pH potentially influence the MDH1–CIT1 interaction in yeast mitochondria.

Figure 5 with 1 supplement see all
Effects of pH and metabolites on the yeast MDH1–CIT1 multienzyme complex affinity.

The affinity of the MDH–CS multienzyme complex was analyzed by microscale thermophoresis (MST) using fluorescently labeled MDH1 as the target and CIT1 as the ligand. Curves represent the response (fraction bound) against CIT1 concentration. Data is presented as mean ± s.d. (N = 3). (A) Effects of pH. The MDH1–CIT1 interaction was determined in the buffer with pH 7.2 (pink), 6.8 (orange), 6.4 (olive green), 6.0 (green), and 5.8 (blue). (B) Effects of 10 mM malate (red), α-ketoglutarate (green), succinate (brown), citrate (blue), aspartate (purple), glutamate (pink), and fumarate (orange). The Kd values of MDH1–CIT1 interaction were shown next to the legend.

Discussion

Here, we report the relationships between respiratory activity and mitochondrial MDH1–CIT1 multienzyme complex interaction in living yeast cells. MDH1–CIT1 multienzyme complex dynamically associated and dissociated in response to respiratory status. This study is the first in vivo evidence in any organism showing the dynamic association and dissociation of the TCA cycle metabolon in real time. The MDH–CS metabolon, specifically, is considered essential for the forward TCA cycle flux due to its ability to overcome the unfavorable thermodynamics of MDH reaction under physiological oxaloacetate concentrations (Beeckmans and Kanarek, 1981; Huang et al., 2018). The results of this study demonstrate that the MDH1–CIT1 multienzyme complex associates when the TCA cycle is active and dissociates when it is suppressed, supporting the idea that metabolon dynamics are related to pathway regulation. Given that the MDH–CS metabolon is proposed to enhance pathway reactions (Bulutoglu et al., 2016; Fukuda et al., 2008), its dynamic assembly and disassembly would influence pathway activity. Together, these observations align with the proposed theory that metabolon dynamics act as a regulatory mechanism of metabolic flux (Noor et al., 2014; Obata, 2020).

The MDH1–CIT1 complex dissociated under conditions suppressing the TCA cycle. TCA cycle inhibition by arsenite and AOA-induced dissociation of the MDH1–CIT1 complex (Figure 3B, C). In a physiological condition, the MDH1–CIT1 complex dissociated when Crabtree-inducing fermentable sugars were supplied to the media (Figures 1 and 2, Figure 2—figure supplements 13). The Crabtree effect redistributes respiratory flux to stimulate fermentation and repress aerobic respiratory pathways, including the TCA cycle (Gancedo, 1998; Ronne, 1995). The metabolic shift by the Crabtree effect is a concerted action of multiple mechanisms, such as the suppression of enzyme gene expressions and the inhibition of mitochondrial transporters (Ahmadzadeh et al., 1996). Interestingly, glucose-induced rapid repression of mitochondrial respiration occurs at a rate that is too great to be explained only by the inhibition of enzyme synthesis, and additional forms of regulation are anticipated to be involved (Ahmadzadeh et al., 1996). The dissociation of the TCA cycle metabolon coincided with this repression (Figure 2A), highlighting metabolon dynamics as a potential additional mechanism to swiftly downregulate mitochondrial respiration during the Crabtree effect.

On the other hand, the MDH1–CIT1 complex interaction was enhanced under conditions that activate the TCA cycle. An increase in TCA cycle activity by the addition of acetate (Orlandi et al., 2014; Vilela-Moura et al., 2011) enhanced the MDH1–CIT1 complex interaction (Figure 3A), suggesting the role of the MDH–CS metabolon in facilitating TCA cycle activity.

However, ETC inhibitions also increased MDH1–CIT1 interaction (Figure 4A–C), despite their inhibitory effects on oxidative respiration (Figure 4—figure supplement 1) and the TCA cycle. This apparent contradiction can be explained by the transient nature of the response. The rise in MDH1–CIT1 interaction following ETC inhibition was transient, suggesting that it does not represent steady-state changes in respiratory activity. Instead, it is more likely linked to acute changes in the mitochondrial matrix microenvironment. Notably, the timing of the increased interaction coincided with changes in matrix pH, supporting a role for pH-dependent regulation (see below). These findings align with our model, which posits that the MDH1–CIT1 interaction is controlled by microenvironmental changes accompanying shifts in respiration. In the case of ETC inhibition, transient perturbations of the matrix microenvironment appear to be disconnected from steady-state respiratory output, thereby promoting MDH1–CIT1 interaction despite overall inhibition of respiration.

We also investigated the molecular mechanisms regulating MDH1–CIT1 complex dynamics. We previously reported that protein conformation changes induced by the solution pH and allosteric regulators affect porcine MDH–CS complex interaction in vitro (Omini et al., 2021). The recombinant yeast MDH1 and CIT1 showed similar responses to the porcine enzymes; MDH1–CIT1 affinity was higher in low pH and in the presence of upstream TCA cycle intermediates (malate and fumarate) but disrupted in the presence of reaction product (citrate; Figure 5). These factors can affect the TCA cycle metabolon formation in vivo. The relationship between these factors and the MDH1–CIT1 interaction observed in this study is summarized in Figure 6A.

Relationship between the respiratory metabolism and the MDH1–CIT1 metabolon association.

(A) Summary of the effects of metabolic treatments on MDH1–CIT1 complex association, mitochondrial matrix parameters (pH, redox state, and ATP levels), and cellular metabolite levels. Each row corresponds to a specific metabolic perturbation. Blue and red labels indicate treatments that decrease or increase the complex association, respectively. Upward and downward arrows indicate increases and decreases in each parameter, respectively. Bold arrows denote changes that are consistent with the observed in vivo alterations in MDH1–CIT1 interaction. The color of the parameter and metabolite indicates its effect on MDH1–CIT1 interaction in vitro (blue, inhibitory; red, promotive). αKG, α-ketoglutarate; fum, fumarate; succ, succinate. (B) A diagram depicting the proposed regulatory mechanism of the MDH1–CIT1 metabolon association. In conditions with low respiratory flux, the MDH1–CIT1 multienzyme complex dissociates, and the tricarboxylic acid (TCA) cycle flux reduces. Reduced electron transport chain (ETC) flux results in higher mitochondrial matrix pH, which reduces MDH1–CIT1 affinity. When the respiratory flux and the TCA cycle flux are high, MDH1–CIT1 metabolon associates and likely channels the intermediate oxaloacetate (OXA). High ETC flux lowers mitochondrial matrix pH and enhances the MDH1–CIT1 interaction. The TCA cycle intermediates affect MDH1–CIT1 metabolon formation; fumarate and malate enhance (red arrows with dotted lines), and citrate inhibits (blue arrows with dotted lines) the interaction. The arrow thickness represents the metabolic fluxes.

We observed an enhanced MDH1–CIT1 interaction under conditions that lower mitochondrial matrix pH in vivo. The mitochondrial matrix pH decreased after adding acetate, corresponding to the enhanced MDH1–CIT1 interaction (Figure 3A, D). The time courses of the temporal increase in MDH1–CIT1 interaction and pH decrease were inversely related following the complex III and IV inhibitions (Figures 4B, C, E, F, 6A). The decrease in the pH from 7.5 to 6.0 ranges is remarkable since the Kd of MDH1–CIT1 association decreases by an order of magnitude within this pH range (Figure 5A). The mitochondrial matrix acidification under respiratory conditions likely favors the MDH1–CIT1 interaction. On the other hand, the mitochondria matrix pH was maintained around 7.2 in the presence of glucose (Figure 2D), which likely weakens the MDH1–CIT1 complex interaction. These results indicate that mitochondrial matrix pH, which is related to proton transport activity by the ETC, can stabilize or destabilize the MDH1–CIT1 complex, thereby regulating the MDH1–CIT1 interaction in response to respiratory activity (Figure 6B). Enhancing the MDH1–CIT1 interaction with the uncoupler CCCP (Figure 4—figure supplement 4) also supports the prominent effect of pH.

Nevertheless, mitochondrial pH may influence the MDH1–CIT1 association, but it is not always the predominant factor regulating the interaction. The observed changes in pH and interaction do not consistently occur simultaneously following the Crabtree induction and acetate supplementation (Figures 2 and 3). TCA cycle inhibitions reduce the MDH1–CIT1 interaction without significantly affecting mitochondrial matrix pH (Figure 3B, C, E, F). These findings suggest that mitochondrial matrix pH can be a contributing factor, either favoring or disfavoring the interaction, rather than the primary regulator, which appears to be the direct induction and dissociation of the MDH1–CIT1 complex.

Little evidence in this study supports a primary role for mitochondrial matrix redox state or ATP levels in controlling MDH1–CIT1 interaction. Our in vitro analyses demonstrate that NAD(H) and ATP can modulate MDH1–CIT1 interaction (Figure 5—figure supplement 1B; Omini et al., 2021), suggesting that these factors may contribute to complex dynamics. Although the lower matrix ATP level following glucose supplementation (Figure 2D) and matrix reduction upon arsenite and AOA-induced TCA cycle inhibition (Figure 3H, I) associated with decreased interaction (Figures 2A and 3B and C), these relationships were not consistently observed across conditions (Figure 6A), indicating that neither redox state nor ATP levels predict complex association in vivo. However, we cannot rule out their contribution under specific metabolic contexts not captured in the current study.

The TCA cycle intermediates and cofactors can also regulate the MDH1–CIT1 interaction, considering their effects on the interaction in vitro (Omini et al., 2021). The yeast enzymes showed responses to malate, α-ketoglutarate, succinate, and citrate (Figure 5B, Figure 5—figure supplement 1) similar to the enzymes of other organisms (Omini et al., 2021; Tompa et al., 1987; Wu et al., 2015), while fumarate is newly identified as an effector of the MDH–CS interaction. Especially, fumarate, malate, and citrate showed significant influences on yeast enzymes (Figure 5B). Cellular levels of these effector metabolites significantly altered under the conditions tested in this study (Figures 2E, 3M). Increased and decreased levels of malate and fumarate following the TCA cycle and ETC inhibitions (Figures 3M and 4M), respectively, are likely related to the MDH1–CIT1 interaction since malate enhances the interaction (Figure 5D; Omini et al., 2021). These results indicate the involvement of metabolite effectors, such as malate and fumarate, in the regulation of the MDH1–CIT1 interaction (Figure 6). However, their precise effects must be evaluated through site-specific, time-dependent perturbation and quantification of metabolite levels in the mitochondrial matrix, as a whole-cell metabolite profile may not reflect the metabolite concentrations that directly influence the complex association.

Aside from MDH1–CIT1 interaction dynamics, our results highlight the complex regulation of TCA cycle metabolism. The activation and inhibition of respiratory activity resulted in diverse metabolic phenotypes (Figure 6A, Source data 1), where intermediate levels did not simply reflect overall pathway activity. This complexity stems from the distinct mechanisms of each inhibitor, such as arsenite affecting α-ketoglutarate dehydrogenase and AOA disrupting the malate-aspartate shuttle (Cavero et al., 2003; Eto et al., 1999; Lee et al., 2011; Liu and Butow, 2006), off-target effects of the inhibitors, and from the adaptive reorganization of intersecting metabolic networks to bypass local blockades (Herrgård et al., 2008; Lehtinen et al., 2013; Rogers et al., 2021). These diverse metabolic phenotypes allow us to assess the relationships between metabolites and metabolon assembly independently of respiratory activity.

This study demonstrates that the TCA cycle MDH1–CIT1 multienzyme complex dynamically interacts in response to the cellular respiratory status. Cues of cellular respiratory state may be transmitted to the multienzyme complex, at least in part, by mitochondrial matrix microenvironment and metabolite levels. In particular, mitochondrial matrix pH and malate and fumarate levels likely have significant effects on the stability of the MDH1–CIT1 complex, since these factors strongly influence complex affinity, and their in vivo changes coincide with complex dynamics (Figure 6). However, none of these factors consistently correlates with the MDH1–CIT1 interaction (Figure 6A), suggesting that none of them is a predominant regulator, but multiple factors work together in a cooperative manner to regulate the formation of the multienzyme complex. This coordinated regulation probably serves to fine-tune the flux of the TCA cycle, allowing it to adapt efficiently to varying metabolic demands and maintain cellular homeostasis. Although we focused on allosteric regulators in this study, further factors are potentially involved in the MDH1–CIT1 complex regulation. For example, 44 and 33 post-translational modifications have been identified in CIT1 and MDH1, respectively (Bhagwat et al., 2021; Henriksen et al., 2012; Holt et al., 2009; Lanz et al., 2021; Reinders et al., 2007; Swaney et al., 2013; Weinert et al., 2013), some of which likely affect the MDH1–CIT1 complex affinity. Various scaffolding molecules, such as long noncoding RNAs, lipid layers, and scaffolding proteins, have been shown to stabilize the multienzyme complexes in other systems (Zhu et al., 2022). Future studies should investigate the effects of these factors to understand the regulatory mechanisms for MDH1–CIT1 interaction.

Changes in environmental conditions and nutrient availability occur swiftly, and a corresponding change in metabolic flux must occur at a similar rate for the successful adaptation of living cells. Dynamic metabolons can be a system to regulate and fine-tune metabolic network flux quickly, allowing cells to maintain metabolic homeostasis in rapidly fluctuating environments (Obata, 2019). Considering the thermodynamic unfavorability of the forward MDH reaction and that substrate channeling overcomes this thermodynamic barrier in vitro, the MDH1–CIT1 complex formation likely enhances the forward TCA cycle flux (Bulutoglu et al., 2016; Fukuda et al., 2008; Sweetlove et al., 2010). The dynamic assembly of the metabolon probably provides fine-tuning of TCA cycle activity in response to respiratory demand. While the present study delineates the environmental and metabolic factors governing this dynamic assembly, future studies utilizing targeted genetic tools to independently modulate specific variables will be critical for dissecting the contributions of individual regulatory factors in vivo. Further studies assessing the impacts of metabolon formation on metabolic pathway flux in living cells are essential to understand the functions of metabolons in metabolic network regulation. Elucidating the regulatory system of the TCA cycle metabolon can lead to a novel strategy to manipulate the TCA cycle flux in metabolic engineering to achieve efficient industrial production of various molecules requiring TCA cycle intermediates as substrates or to control the Warburg effect to suppress cancer cells.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (S. cerevisiae)CIT1SGDSGD:S000005284YNR001C
Gene (S. cerevisiae)MDH1SGDSGD:S000001568YKL085W
Strain, strain background (S. cerevisiae)BY4741SGDMATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0Background strain
Genetic reagent (S. cerevisiae)MDH1/CIT1-BiTThis paperMDH1-SmBiT-cMyc and CIT1-LgBiT-HA
Genetic reagent (S. cerevisiae)MDH1-nLUCThis paperMDH1 fused C-terminally to full-length NanoLUC-HA
Genetic reagent (S. cerevisiae)CIT1-nLUC strainThis paperCIT1 fused C-terminally to full-length NanoLUC-HA
Genetic reagent (S. cerevisiae)MDH1/CIT1-BiT + mito GO ATeamThis paperMDH1/CIT1 reporter with mitochondrial ATP sensor
Genetic reagent (S. cerevisiae)MDH1/CIT1-BiT + mito roGFPThis paperMDH1/CIT1 reporter with mitochondrial redox sensor
Genetic reagent (S. cerevisiae)MDH1/CIT1-BiT + pHluorinThis paperMDH1/CIT1 reporter with mitochondrial pH sensor
Antibodyanti-cMyc (mouse monoclonal)Thermo ScientificRRID:AB_55847WB (1:10,000)
Antibodyanti-cMyc (mouse monoclonal)Thermo ScientificRRID:AB_2533049WB (1:5000)
Antibodyanti-PGK1 (mouse monoclonal)Thermo ScientificRRID:AB_2532235WB (1:5000)
Antibodyanti-DLD1 (rabbit polyclonal)Baile et al., 2013WB (1:5000)
Recombinant DNA reagentpDGB2Ω2 vectorSarrion-Perdigones et al., 2013Golden Braid vector
Recombinant DNA reagentpGDBα1 vectorSarrion-Perdigones et al., 2013Golden Braid vector
Recombinant DNA reagentpGDB2Ω2-1; pGDB2Ω2-2; pGDB2Ω2-3This paperModified cloning-site vectors
Recombinant DNA reagentpGDB2Ω2-5′LgBiT; pGDB2Ω2-3′LgBiT; pGDB2Ω2-5′SmBiT; pGDB2Ω2-3′SmBiT; pGDB2Ω2-5′nLUC; pGDB2Ω2-3′nLUC; pGDB2Ω2-URA3This paperIntermediate Golden Braid vectors
Recombinant DNA reagentpGDBα1-LgBiT; pGDBα1-SmBiT; pGDBα1-nLUCThis paperIntegration-cassette plasmids
Recombinant DNA reagentnLUC-HA; LgBiT-HA; SmBiT-cMycThis paperGenewizSynthetic DNA fragments
Recombinant DNA reagentpET21bNovagen69741-3 (Millipore Sigma)E. coli expression plasmid
Recombinant DNA reagentpET-cit1; pET-mdh1This paperVectors for recombinant MDH1 and CIT1 expression
Recombinant DNA reagentp415-GPDpro-mito GO ATeam; p416-GPDpro-mito roGFPVevea et al., 2013Mitochondrial ATP and redox biosensor plasmids
Recombinant DNA reagentnLUC-HA; pAG416-COX4-pHluorinAyer et al., 2013Mitochondrial matrix pH biosensor plasmid
Sequence-based reagent40 PCR and sequencing primersMillipore SigmaFull nucleotide sequences are provided in Appendix 1
Commercial assay or kitNano-Glo Luciferase Assay SystemPromegaN1110Furimazine substrate for NanoBiT/NanoLUC luminescence assays
Commercial assay or kitProtein Labeling Kit RED-NHS 2nd GenerationNanoTemperMO-L011Labeling kit for microscale thermophoresis assay
Software, algorithmMassHunter Software suitesAgilentRRID:SCR_016657;
RRID:SCR_015040;
RRID:SCR_019081
GC–MS data acquisition and analysis

Strains, media, and culture conditions

Request a detailed protocol

S. cerevisiae BY4741 (MATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0) was used as the background strain. Cells were grown in synthetic complete medium (SD) containing 0.67% yeast nitrogen base lacking amino acids (Research Products International, Mt. Prospect, IL, USA) with 2% (wt/vol) D-raffinose and 1% amino acid complete drop-out mix. The complete amino acid drop-out mix was replaced with the amino acid drop-out mix lacking leucine for GoAteam expressing cells or the amino acid drop-out mix lacking uracil for RoGFP and pHluorin expressing cells. Cell growth cultures were incubated with shaking at 220 rpm at 28°C in an incubator shaker. Cell cultures were grown to exponential phase with an OD600 = 0.5–1.0 in all analyses.

Oxygen consumption rate measurement

Request a detailed protocol

Oxygen consumption was measured using a Clark-type electrode (Oxygraph, Hansatech Instruments, Norfolk, UK) as described previously (Agrimi et al., 2011). Cells were grown to exponential phase OD600=0.5, harvested, centrifuged at 3000 x g for 5 min at 4°C, and resuspended in growth medium to obtain a density of OD600 = 5.0. A total of 1 ml reaction volume consisting of equal parts of 445 µl cell culture, 445 µl growth media, and 10 µl of each inhibitor was added to the oxygraph chamber. Sodium malonate (20 mM), antimycin A (10 µM), sodium cyanide (0.5 mM), and oligomycin (1 mM) were used as specific inhibitors of ETC complexes II, III, IV, and V, respectively. The change in oxygen concentration was followed subsequently. Oxygen consumption rates were determined before and after the addition of each inhibitor from the slope of a plot of O2 concentration versus time. All measurements were conducted in triplicates.

Generation of the S. cerevisiae reporter strains with tagged MDH1 and CIT1

Request a detailed protocol

The NanoBiT split NanoLUC luciferase complementation system (Dixon et al., 2016) was adopted to monitor MDH1–CIT1 interaction in yeast cells. The yeast lines expressing MDH1 and CIT1 proteins fused with the small (SmBiT) and large (LgBiT) NanoBiT subunits, respectively, were generated by inserting the yeast codon-optimized tag-coding sequences into the BY4741 genome following the scarless C-terminal tagging procedure (Landgraf et al., 2016). The full-length NanoLUC, LgBiT, and SmBiT coding sequences were integrated to direct downstream of the Mdh1 (chrIV:3300230) and Cit1 (chrX:303993) genes on the BY4741 genome for C-terminal fusion. The LgBiT and SmBiT sequences include flexible linkers (ACKIPNDLKQKVMNH; Hu et al., 2002) with HA and cMyc epitope sequences, respectively, for immunological detections. Details of the vector construction procedure, primers, and vectors are described in Appendix 1 (Supplementary method).

Briefly, we generated the plasmids with integration cassettes composed of the linker, 5′ half of tag-coding sequence, URA3 selection marker, and 3′ half of tag-coding sequence in this order using the Golden Braid technology (Sarrion-Perdigones et al., 2013). The halves of the tag-coding sequences have an overlapping sequence. The integration cassettes were amplified by PCR using the primers with around 40 bases 5′ extension with sequences homologous to the franking region of the tag insertion sites. BY4741 cells were transformed with the purified PCR products to insert the cassettes into the target sites directly downstream of the Mdh1 and Cit1 genes by homologous recombination due to the 5′ extension of the primers. The transformants were selected on uracil-deficient plates. To pop-out the URA3 marker and reconstruct NanoLUC and NanoBiT subunits, the transformants were cultured overnight in uracil-containing YPD media. The URA3 gene was excluded from the genome by recombination based on the overlapping sequences in the cassette. The cell suspension was spread on SD plates containing 1 mg ml–1 5-fluoroorotic acid to select the cell lines without the URA3 gene. SmBiT was initially fused with the Mdh1 gene, and the resulting strain was further transformed to fuse LgBiT with Cit1 (MDH1/CIT1-BiT strain). Another strain expressing MDH1 fused with full-length NanoLUC luciferase was also generated following the same procedure to monitor MDH1 protein levels (MDH1-nLUC strain).

In vivo MDH1–CIT1 interaction measurement

Request a detailed protocol

The MDH1/CIT1-BiT strain expressing MDH1 and CIT1 enzymes tagged with split halves of luciferase enzyme was grown to OD600 = 0.35–0.45. Cells were collected and resuspended to obtain a cell density of OD600 = 2.0 in fresh SD-Raff media. Each sample consisted of 80 µl of media, 10 µl of cells, and 10 µl of 50x furimazine luciferase substrate (Promega, Madison, WI, USA). The luminescence signal was measured every five min with a microplate reader (CLARIOSTAR Plus, BMG LABTECH, Ortenberg, Germany) at 28°C with 200 rpm shaking at the beginning of each cycle. Baseline luminescence was measured for 20 min before the treatment was applied. The time-dependent luminescence was measured for an additional 80 min after treatment. Relative luminescence unit (RLU) was calculated by normalizing the luciferase signals by the average signals during three pre-treatment time points. The protein levels of MDH1 and CIT1 were assessed using yeast lines expressing these proteins tagged with full-length NanoLUC luciferase. Using the luminescence data of these relative protein levels, we calculated the normalized interaction index by dividing the NanoBiT interaction signal by the product of the relative abundances of both proteins:

NormalizedInteractionIndex=NanoBiTMDH1CIT1

In this formula, NanoBiT, MDH1, and CIT1 are the relative luminescence levels at each time point.

Western blotting

Request a detailed protocol

Cells were grown to OD600 = 0.5 and cell pellet was collected. Cells were first pretreated with 2 M LiAc and 0.4 M NaOH to permeabilize the cells and then treated with SDS–PAGE sample buffer to extract proteins according to the method described by Zhang et al., 2011. The cell lysate was centrifuged at 27,000 × g for 10 min at 4°C and the proteins were detected by SDS–PAGE and western blotting following the method described previously (Rajasekaran et al., 2022). cMyc Tag Monoclonal Antibody (MA1213161MG, Thermo Scientific) and HA Tag Monoclonal Antibody (LSG26183, Thermo Fisher Scientific) were used as primary antibodies to detect MDH1 and CIT1, respectively. DLD1 antibody generated by S. Claypool laboratory at the Johns Hopkins University (Baile et al., 2013), and PGK1 Monoclonal Antibody (Invitrogen 459250, Thermo Fisher Scientific) was used to detect phosphoglycerate kinase as a housekeeping protein and internal standard.

In vitro MDH1–CIT1 interaction measurement

Request a detailed protocol

The recombinant MDH1 and CIT1 were produced according to the method described in Appendix 1. The interaction between recombinant MDH1 and CIT1 enzymes was analyzed by microscale thermophoresis (MST) according to the method described by Omini et al., 2021 with slight modifications. Base MST buffer contained 50 mM Tris-HCl (pH 8), 150 mM NaCl, 10 mM MgCl2, 5 mM DTT, and 0.05% Tween-20. Recombinant MDH1 (10 µM) was labeled with Protein Labeling Kit RED-NHS 2nd Generation (NanoTemper, München, Germany) and used as the target. CIT1 was used as the ligand. Two times serial dilution of 80 µM CIT1 was conducted for 16 concentrations. A total of 10 µl of CIT1, 10 µl MST buffer, and 10 µl labeled MDH1 were mixed and loaded to Monolith NT.115 Capillaries (NanoTemper). Capillaries were incubated at room temperature for 1 min, and the interaction was analyzed by MST using Monolith NT.115 (NanoTemper). To test the effect of pH on MDH1–CIT1 interaction, 50 mM Tris buffers with pH 6.0, 7.5, and 8.0 were prepared and used as the MST buffer for sample preparation. To test the effect of reducing environment on interaction, various concentrations (5, 2.5, and 1.25 mM) of DTT were prepared with the MST buffer and used for sample preparation. To test the effect of ATP level on interaction, an MST buffer containing different concentrations of ATP (5, 2.5, and 1.25 mM) was made and used for sample preparation. To test the effect of metabolite availability on the interaction, MST buffer containing respective metabolites at 10 mM concentration was made and used for sample preparation and experiments.

Enzyme activity assays

Request a detailed protocol

Cells were grown to an exponential phase with an OD600 of 0.5, and 2 ml of the cells were harvested for the enzyme activity assay. Yeast cell lysates were prepared by disruption with glass beads as described previously (Mukherjee et al., 2020), with the lysis buffer omitted butylhydroxytoluene. The protein concentration of the lysate was determined using the Pierce BCA Protein Assay Kit (Thermo Fisher Scientific). CS activity was determined using a method described previously (Srere, 1969) that measures free thiols by coupling the CS reaction to thiol reaction with 5,5′-dithiobis-(2-nitrobenzoate) (DTNB). The CS enzyme activity assay mixture contained 154 µl of distilled water, 20 µl of 1 mM DTNB, and 6 µl of 10 mM acetyl CoA. The CS reaction was initiated by the addition of 10 µl of 10 mM oxaloacetate. The absorption at 412 nm was followed to measure CS activity. MDH activity assay mixture contained 50 mM TES (pH 7.2), 5 mM MgCl2, 0.2 mM NADH, and 0.05% Triton X-100. To obtain a total reaction volume of 300 µl, 285 µl of the assay mixture, and 5 µl cell lysate were added to the wells, and the reaction was initiated with 10 µl of 30 mM oxaloacetate. Reduction of NADH was followed at 340 nm to determine MDH enzyme activity. Enzyme activity was measured using a microplate reader absorbance function (CLARIOSTAR Plus, BMG LABTECH).

Metabolite profiling

Request a detailed protocol

Yeast cells were grown to the exponential phase (OD600 ~ 0.5), and treatments were applied and incubated for the time described in the figure legends. Cellular metabolites were extracted following the protocol described by Obata et al., 2013. The cells in 1 ml culture were harvested by vacuum filtration using a membrane filter (0.45 µm HV Durapore 25 mm diameter; Millipore Sigma, Burlington, MA, USA). The filter was put into a 2-ml microcentrifuge tube, flash-frozen in liquid nitrogen, and stored at −80°C. The metabolites were extracted from the filtered cells with methanol:water:chloroform, and a 50-µl aliquot was dried down by vacuum centrifugation. Dried metabolites were derivatized with methoxyamine hydrochloride in pyridine and further trimethylsilylated by N-methyl-N-(trimethylsilyl) trifluoroacetamide (Millipore Sigma). Derivatized samples were analyzed by 7200 GC-QTOF system (Agilent, Santa Clara, CA, USA) exactly as described in Wase et al., 2022. Each metabolite’s peak height was normalized by the peak height of the internal standard (ribitol) to represent relative levels of metabolite.

Expression of mitochondrial biosensors

Request a detailed protocol

The pH, redox state, and ATP levels in the mitochondrial matrix were measured using the mito-Go Ateam2, mito-ROGFP1, and pHluorin (pAG416-COX4-pHluorin, URA selection marker) fluorescence biosensors, respectively, specifically localizing in the mitochondrial matrix. The mito-Go Ateam2 (p415-GPDpro-mito GO ATeam) and mito-roGFP1 (p416-GPDpro-mito roGFP) encoding plasmids (Vevea et al., 2013) were generous gifts from Dr. Liza Pon at the Department of Pathology and Cell Biology, Columbia University. The pHluorin encoded plasmid (pAG416-COX4-pHluorin; Ayer et al., 2013) was a generous gift from Dr. Anita Ayer at the University of New South Wales, Sydney, Australia. The plasmids encoding these biosensors were transformed into split-Luc tagged MDH1/CIT1 strain using the lithium acetate method (Chen et al., 1992). The p416-GPDpro-mito roGFP and pAG416-COX4-pHluorin harboring cells were selected on URA- media, and the cells with p415-GPDpro-mito GO ATeam were selected on Leu- media.

The localization of biosensors within the mitochondria was confirmed by confocal microscopy. Cells were grown in SD selection media to OD600 = 0.50. Cells were stained with MitoTracker Deep Red FM (Thermo Fisher Scientific) in a 100-nM dye solution for 30 min with shaking at 28°C to visualize the mitochondria. The stained cells were resuspended in fresh 10 mM HEPES buffer (pH 7.4) with 2% raffinose to obtain a final concentration of OD600 = 10. Confocal imaging was performed on an A1R-Ti2 Confocal Laser Scanning Microscope (Nikon, Tokyo, Japan) with a Plan Apo 60x 1.40 Oil lens 0.17 WD 0.13 (Nikon) and a 2x digital zoom for a total magnification of ×1200. Brightfield and fluorescent images were captured. Imaging was done sequentially with excitation at 405 nm and emission at 425–475 nm for pHluorin and mito-roGFP1, excitation at 488 nm and emission at 500–550 nm for mito-GoAteam2, and excitation at 640 nm and emission at 663–738 nm for MitoTracker. Images were collected with NIS Elements software (Nikon) and processed with ImageJ analysis software (Vevea et al., 2013).

Measurement of redox state, pH, and ATP level in the mitochondrial matrix

Request a detailed protocol

Cells expressing pHluorin, mito-roGFP1, and mito-GoAteam2 were grown to OD600 = 0.50 in their respective selection media. Cells were prepared for fluorescence measurement according to the method described by Morgan et al., 2011. Time-based fluorescence intensity was measured using a microplate reader (CLARIOSTAR Plus, BMG LABTECH) at 28°C with shaking at 200 rpm with 20 min of baseline measurement before treatments. In all experiments, the strains harboring empty vectors were grown simultaneously as a reference for background fluorescence at the different excitation wavelengths. Background fluorescence was subtracted from fluorescence intensity from cells expressing biosensors.

For the pH measurement, the ratio of emission intensity at 510 nm resulting from excitation of pHluorin at 390 and 470 nm was calculated (R390/470) using the pHluorin-expressing strain. The mitochondrial matrix pH was calculated from the R390/470 ratio using a calibration curve generated according to the method described by Orij et al., 2009. Briefly, pHluorin-expressing cells were permeabilized by resuspension in PBS containing 100 µg ml−1 digitonin for 10 min and washed with PBS. Cells were resuspended in citric acid/Na2HPO4 buffer of pH values ranging from 5.0 to 9.0. The ratio of pHluorin emission at 510 nm upon excitation at 390 and 470 nm (R390/470) was plotted against buffer pH to obtain a calibration curve.

The mitochondrial matrix redox state was measured using the mito-roGFP1-expressing strain. The ratio of emission intensity at 510 nm from excitation at 365 and 470 nm (R470/365) was measured. Mito-roGFP1 in situ calibration was performed following the method described by Vevea et al., 2013. Digitonin-treated mito-roGFP1-expressing cells were incubated in 0, 5, and 10 mM H2O2 and DTT for 20 min at 28°C with shaking at 200 rpm. The R470/365 was plotted against the redox potential of the solutions calculated using the formula described previously (Morgan et al., 2011).

Relative ATP level was determined by measuring the ratio of the emission intensity of mito-GoAteam2 at 510 and 560 nm with excitation at 470 nm (R560/510). The ATP levels were analyzed as the relative value without calibration to absolute concentrations.

Statistical analysis

Request a detailed protocol

The differences between the control and test samples were evaluated by a two-tailed unpaired Student’s t-test. p < 0.05 was considered a statistically significant difference. For the time course data, the test was applied at each time point. All data were obtained from triplicated independent experiments.

Appendix 1

Supplementary method

General procedures

PCR reactions for the construction of vectors and integration cassettes were conducted with Phusion High-Fidelity DNA polymerase (Thermo Fisher Scientific, Waltham, MA, USA). DreamTaq Green DNA polymerase (Thermo Fisher Scientific) was used for the amplification of vector cloning sites and yeast genome to verify the integration of desired sequences. The Sanger sequencing of the PCR products was conducted by Eurofins Scientific (Luxembourg City, Luxembourg). Restriction digestion and ligation were conducted by manufacturer’s protocols using the enzymes from New England Biolabs (NEB, Ipswich, MA, USA). Chemically competent E. coli DH5α was used for vector construction. Frozen-EZ Yeast Transformation Kit (Zymo Research Corporation, CA, USA) was used for yeast cell transformation. All oligo DNAs were synthesized by Millipore Sigma (Burlington, MA, USA) and Eurofins Scientific. Sequences of the primers used in this study are in Appendix 1—table 1.

Generation of the nanoBiT and nanoLUC reporter yeast strains

Vector construction

The DNA fragments for genomic DNA recombination were constructed in the pGDBα1 vector using the Golden Braid technology (Sarrion-Perdigones et al., 2013). The procedure is summarized in Figure. First, the cloning site in the pDGB2Ω2 vector (Sarrion-Perdigones et al., 2013) was replaced to generate three vectors containing desired 4-base overhang by BsaI digestion. The pDGB2Ω2 cloning site was amplified using the primers with modified overhang sequences and EcoRV recognition sequences on both ends (pGDB2Ω2-1, pGDB2Ω2-2, and pGDB2Ω2-3 primer sets; Appendix 1—table 1). Purified PCR products were cloned into pGDB2Ω2 at the EcoRV site by T4 ligase, resulting in the pGDB2Ω2-1, pGDB2Ω2-2, and pGDB2Ω2-3 vectors including GGAG/TCCT, TCCT/CGTC, and CGTC/CGCT BsaI restriction sites, respectively (Figure). The inserted sequences were amplified using pDGB1seq primer sets and the sequences were verified.

Second, the vectors harboring 5′ half of tag-coding sequence, URA3 selection marker, and 3′ half of tag-coding sequence were generated. The DNA sequence coding linker-HA-nanoBiT large subunit (LgBiT-HA), linker-cMyc-nanoBiT small subunit (SmBiT-cMyc), and linker-HA-full-length nanoLUC luciferase (nLUC-HA) were synthesized by Genewiz (Azenta Life Sciences, South Plainfield, NJ, USA). The sequences of the synthesized DNA are in Appendix 1—table 2. The 5′ and 3′ halves of tag-coding sequences were amplified using the primers with 5′ extensions containing BsaI recognition sequences to generate GGAG/TCCT and CGTC/CGCT overhangs to be cloned into the BsaI sites of pGDB2Ω2-1 and pGDB2Ω2-3, resulting in pGDB2Ω2-5′ and pGDB2Ω2-3′ vectors, respectively. pDGB2Ω2-LgBiT5′ and pDGB2Ω2-LgBiT3′ primer sets were used to amplify the 5′ and 3′ halves of the large nanoBiT subunits using LgBit-HA as the template to construct pDGB2Ω2-5′LgBiT and pDGB2Ω2-3′LgBiT vectors, respectively. pDGB2Ω2- SmBiT 5′ and pDGB2Ω2-SmBiT 3′ primer sets were used to amplify the 5′ and 3′ halves of the small nanoBiT subunits using SmBit-cMyc as the template to construct pDGB2Ω2-5′SmBiT and pDGB2Ω2-3′SmBiT vectors, respectively. Full-length nanoLUC luciferase sequence was amplified using pDGB2Ω2-nLUC5′ and pDGB2Ω2-nLUC3′ primer sets with nLUC-HA as the template to construct pDGB2Ω2-5′nLUC and pDGB2Ω2-3′nLUC vectors, respectively. URA3 gene was amplified from genomic DNA of Kluyveromyces lactis using the primers with BsaI sites with TCCT/CGTC sequences (pDGB2Ω2-Ura primer set) and cloned into pGDB2Ω2-2 to construct the pGDB2Ω2-URA3 vector. The inserted sequences were amplified using the pDGB1seq primer set and their sequences were verified.

Third, the DNA fragments were assembled into the pGDBα1 vector using the Golden Braid procedure. The full-length nanoLUC, LgBiT, and SmBiT coding sequences were integrated directly downstream of the Mdh1 (chrIV:3300230) and Cit1 (chrX:303993) genes on the BY4741 genome for C-terminal fusion. pGDB2Ω2-5′, pGDB2Ω2-URA3, pGDB2Ω2-3′ (300 ng each), and 200 ng pGDBα1 were incubated with 20U BsaI-HF (NEB) and 400U T4DNA ligase (NEB) at 37°C for 4 min, 16°C for 5 min, and 50°C for 5 min for 28 cycles. The reaction was terminated by heating at 80°C for 5 min after the cycles. The reaction mixtures were used to transform E. coli and the colonies were selected on LB plates containing 50 µg/ml kanamycin, 20 µg/ml Xgal (GoldBio), and 25 µg/ml isopropyl-β-D-thiogalactopyranoside (IPTG; GoldBio). The plasmids were extracted from the white colonies. Their inserted sequences were amplified using the pDGB1seq primer set and their sequences were verified by sequencing. pGDB2Ω2-5′LgBiT, pGDB2Ω2-URA3, and pGDB2Ω2-3′LgBiT were used to construct pGDBα1-LgBiT vector. pGDB2Ω2-5′SmBiT, pGDB2Ω2-URA3, and pGDB2Ω2-3′SmBiT vectors were used to construct pGDBα1-SmBiT vector. pGDB2Ω2-5′SmBiT, pGDB2Ω2-URA3, and pGDB2Ω2-3′SmBiT vectors were used to construct pGDBα1-SmBiT vector. pGDB2Ω2-5′nLUC, pGDB2Ω2-URA3, and pGDB2Ω2-3′nLUC vectors were used to construct pGDBα1-nLUC vector.

Scarless integration of tagging sequences into yeast genome

The integration cassettes were PCR-amplified using the primers annealing directly downstream of the BsaI sites with 37-40 base 5′ extension complementary to the franking sequences of the genomic insertion sites. The fragment to integrate small nanoBiT subunit into the downstream of MDH1 gene was amplified using MDH1-SmBiT primer set with pGDBα1-SmBiT vector as the template. The fragment to integrate large nanoBiT subunit into the downstream of CIT1 gene was amplified using CIT1-LgBiT primer set with pGDBα1-LgBiT vector as the template. The fragment to integrate full-length nanoLUC coding sequence into the downstream of MDH1 gene was amplified using MDH1-nLUC primer set with pGDBα1-nLUC vector as the template. Yeast cells were transformed with the PCR products and the strains with integrated sequences were selected on URA- plates. The genomic DNA was isolated from the positive clones and the tag sequence insertion was assessed by PCR using the primers annealing the adjacent region of the insertion sites (MDH1seq and CIT1seq primer sets) and the URA3 coding sequence (Ura3intFw/MDH1seqRv and Ura3intFw/CIT1seqRv primer sets).

The URA3 marker was ejected by endogenous homologous recombination between the repeats in the inserted sequence. The verified transformants were cultured in liquid YPD media containing uracil for 2 days. The strains with removed URA3 marker were selected on the URA-containing SD plate with 1 mg ml–1 5-fluoroorotic acid. The lack of URA3 marker was assessed by PCR using MDH1seq and CIT1seq primer sets. Correct integration of the inserted fragments was confirmed by sequencing.

Recombinant MDH1 and CIT1 protein production

Vector construction

The Cit1 and Mdh1 coding sequences without the mitochondrial targeting signal were amplified from S. cerevisiae BY4741 cDNA library generated by SuperScript III First-Strand Synthesis System (Invitrogen, Thermo Fisher Scientific) using the cit1-Fw/cit1-Rv and mdh1-Fw/mdh1-Rv primer sets (Appendix 1—table 1). The cDNA fragments were cloned into the NdeI/XhoI restriction site in the expression plasmid pET21b with a C-terminal 6x His tag, using the ClonExpress II One Step Cloning kit (Nanjing Vazyme Biotech Co. Ltd., China) following the manufacturer’s instructions. The resulting plasmids were extracted, and the inserts were amplified by cit1-verFw/ctmd-verRv and mdh1-verFw/ctmd-verRv primer sets (Appendix 1—table 1) to verify the inserted sequences. The plasmids for the expression of CIT1 (pET-cit1) and MDH1 (pET-mdh1) were transformed into chemically competent E. coli BL21 (DE3) cells.

Protein expression and purification

The E. coli BL21 (DE3) strains harboring pET-cit1 and pET-mdh1 were cultivated in terrific broth medium at 37°C with shaking at 220 rpm to OD600 = 0.6–0.8, and 0.8 and 0.2 mM of IPTG was applied for CIT1 and MDH1 expression induction, respectively. The cultures were incubated for an additional 20 hr at 18°C, and the cells were harvested by centrifugation at 6000 × g and 4°C for 20 min. The cell pellets were resuspended in cell lysis buffer (50 mM Tris-HCl buffer, pH 8.0, 500 mM NaCl, 10% glycerol, 20 mM imidazole, 1 tablet of cocktail of EDTA-free protease inhibitors/100 ml, 200 µg/ml lysozyme, 20 µg/ml DNase, 1 mM MgCl2, 0.5% Tween 20, and 10 mM DTT) using fivefold volume of the weight of cell pellets. The cells were then disrupted by sonication using Sonic Dismembrator 550 (Fisher Scientific, Hampton, NH, USA) with the amplitude set at 6, pulser on 5 s, pulser off 15 s and a total process time of 7 min. Following the centrifugation at 16,000 × g at 4°C for 70 min, recombinant proteins in the supernatants were purified using Ni-NTA His SpinTrap columns (GE Healthcare, Chicago, IL, USA), and the proteins were eluted with concentration-gradient imidazole in protein purification buffer (50 mM Tris-HCl buffer, pH 8.0, 500 mM NaCl, 10 mM β-mercaptoethanol, 20 mM imidazole (binding buffer), 40, 60, and 80 mM imidazole (wash buffers), and 500 mM imidazole (elution buffer)). The purified enzymes were desalted using an Amicon Ultra-4 centrifugal concentrator (10 kDa MWCO, Millipore Sigma) at 4°C to exchange to buffer A (50 mM Tris-HCl buffer, pH 8.0). Enzyme purity was estimated using SDS–PAGE with Coomassie blue staining. Glycerol was added to the purified proteins to a final concentration of 10%. The purified proteins were then flash-frozen with liquid nitrogen and stored at −80°C until further use. Protein concentration was determined using the Bio-Rad Protein Assay Kit (Bio-Rad Laboratories, Hercules, CA, USA) following the manufacturer’s instructions. The enzyme activity was verified as described in the Method section.

Appendix 1—figure 1
Schematics of vector construction and genomic integration procedures.

Colored texts indicate the overhang sequences of the BsaI digestion sites. The ‘tag’ indicates the integrated tag-coding sequences, namely large nanoBiT subunit, small nanoBiT subunit, and nanoLUC luciferase. GOI, gene of interest; CDS, coding sequence; gDS, genomic downstream sequence; Genome, genomic sequence.

Appendix 1—table 1
Primers used in this study.

Fw and Rv indicate forward and reverse primers, respectively. Bold letters indicate the 4-base overhang in the BsaI digestion sites. Underlined sequences indicate plasmid-annealing sequences.

Primer namePrimer sequence
pGDB2Ω2-1FwTCATGATATCGGTCTCAGGAGCACAGCTTGTCTGTAAGCG
pGDB2Ω2-1RvCATTGATATCGGTCTCAAGGACAGCTGGCACGACAGGTTTC
pGDB2Ω2-2FwGATATCTCCTAGAGACCCACAGCTTGTCTGTAAGCGG
pGDB2Ω2-2RvCATTGATATCGACGAGAGACCCAGCTGGCACGACAGGTTTC
pGDB2Ω2-3FwCATTGATATCGGTCTCACGTCCACAGCTTGTCTGTAAGCGG
pGDB2Ω2-3RvCATTGATATCGGTCTCAAGCGCAGCTGGCACGACAGGTTTC
pDGB2Ω2-LgBiT5′FwCATTCATTGGTCTCAGGAGCAGCTTGTAAGATCCCAAACG
pDGB2Ω2-LgBiT5′RvGTTGTATTGGGTCTCAAGGAGTAACACCATCAATAACTAAAGTACCA
pDGB2Ω2-UraFwGGTGTTACTCCTTGAGACCCAATACAACAGATCACGTG
pDGB2Ω2-UraRvAGCAGTTTGACGTGAGACCCGTTTTATTTAGGTTCTATCGAGG
pDGB2Ω2-LgBiT3′FwCATTACGGTCTCACGTCCAAACTGCTGCT
pDGB2Ω2-LgBiT3′RvCTAGAGTAGGTCTCTAGCGGCTAAAATACGTTCACATAAACGCCAACC
pDGB2Ω2-SmBiT5′FwATAGGTCTCAGGAGGCATGTAAAATTCCTAATGATTTAA
pDGB2Ω2-SmBiT5′RvGAACTAGGTCTCAAGGATAAAATTTCTTCAAATAAACGATAACCAG
pDGB2Ω2-SmBiT3′FwCATTACGGTCTCACGTCGCATGTAAAATTCC
pDGB2Ω2-SmBiT3′RvCATAAGGGTCTCTAGCGTTATAAAATTTCTTCAAATAAACGATAA
pDGB2Ω2-nLUC5′FwCATTCATTGGTCTCAGGAGCAGCTTGTAAGATCCCAAACG
pDGB2Ω2-nLUC5′RvGTTGTATTGGGTCTCAAGGAGTAACACCATCAATAACTAAAGTACCA
pDGB2Ω2-nLUC3′FwATAAAACGGGTCTCACGTCAAACTGCTGGTTATAATTTAGATCAAG
pDGB2Ω2-nLUC3′RvCTAGAGTAGGTCTCTAGCGGCTAAAATACGTTCACATAAACGCCAACC
MDH1-SmBiT-FwAAGAATATCGAAAAGGGTGTCAACTTTGTTGCTAGTAAAGCATGTAAAATTCCTAATGA
MDH1-SmBiT-RvTTTTTTTCCCTATTTTTCACTCTATTTCTGATCTTATAAAATTTCTTCAAATAAACG
CIT1-LgBiT-FwAATACAAGGAGTTGGTAAAGAAAATCGAAAGTAAGAACGCTTGTAAGATCCCAAACGAC
CIT1-LgBiT-RvTGAAAATACGTGTTTGAATAGTCGCATACCCTGAATCTCGTGTTACTATTAATTCTTA
MDH1-nLUC-FwAAGAATATCGAAAAGGGTGTCAACTTTGTTGCTAGTAAAAGACCAGCTTGTAAGATCCCA
MDH1-nLUC-RvTTCCCTATTTTTCACTCTATTTCTGATCTTGAACAATTTAAGCTAAAATACGTTCACA
MDH1seqFwGCATCTCCGGTCACTTTGGG
MDH1seqRvCTAGTTGATTTTTGGCAGTTTCCTTCCTTTC
CIT1seqFwGCCAGAGCTATTGGTGTGTTACC
CIT1seqRvCGGTAGGCATAGGGGACTCAAAG
Ura3intFwGCGTTACCACCATCCAATGCAGAC
pDGB1seqFwCGCCAGCAACGCGGCCTT
pDGB1seqRvGCAAGCGGTTGCCCACCGTC
cit1-FwCTTTAAGAAGGAGATATACATATGAGTAGCGCCTCCGAACAAACG
cit1-RvAGTGGTGGTGGTGGTGGTGCTCGAGGTTCTTACTTTCGATTTTCTTTACCAACTCCTTG
mdh1-FwCTTTAAGAAGGAGATATACATATGTATAAAGTGACTGTTTTGGGTGCAGGC
mdh1-RvAGTGGTGGTGGTGGTGGTGCTCGAGTTTACTAGCAACAAAGTTGACACCCTTTTC
cit1-verFwGGGCTACGAAAACAAGGATTTTATTGAC
mdh1-verFwCATCAACGCAAGCATCGTTC
ctmd-verRvGTTATGCTAGTTATTGCTCAGCGGTG
Appendix 1—table 2
Sequences of the synthesized DNA.
DNA nameSequences
nLUC-HACCCGGGAGACCAGCTTGTAAGATCCCAAACGACTTGAAGCAAAAGGTTATGAACCACTACCCATACGACGTACCAGATTACGCTATGGTTTTTACTTTAGAAGATTTTGTTGGTGATTGGCGTCAAACTGCTGGTTATAATTTAGATCAAGTTTTAGAACAAGGTGGTGTTTCTTCTTTATTTCAAAATTTAGGTGTTTCTGTTACTCCTATTCAACGTATTGTTTTATCTGGTGAAAATGGTTTAAAAATTGATATTCATGTTATTATTCCTTATGAAGGTTTATCTGGTGATCAAATGGGTCAAATTGAAAAAATTTTTAAAGTTGTTTATCCTGTTGATGATCATCATTTTAAAGTTATTTTACATTATGGTACTTTAGTTATTGATGGTGTTACTCCTAATATGATTGATTATTTTGGTCGTCCTTATGAAGGTATTGCTGTTTTTGATGGTAAAAAAATTACTGTTACTGGTACTTTATGGAATGGTAATAAAATTATTGATGAACGTTTAATTAATCCTGATGGTTCTTTATTATTTCGTGTTACTATTAATGGTGTTACTGGTTGGCGTTTATGTGAACGTATTTTAGCTGTCGAC
LgBiT-HAGGATCCAGACCAGCTTGTAAGATCCCAAACGACTTGAAGCAAAAGGTTATGAACCACTACCCATACGACGTACCAGATTACGCTATGGTTTTTACTTTAGAAGATTTTGTTGGTGATTGGGAACAAACTGCTGCTTATAATTTAGATCAAGTTTTAGAACAAGGTGGTGTTTCTTCTTTATTGCAAAATTTAGCTGTTTCTGTTACTCCTATTCAACGTATTGTTAGATCTGGTGAAAATGCTTTAAAAATTGATATTCATGTTATTATTCCTTATGAAGGTTTATCTGCTGATCAAATGGCTCAAATTGAAGAAGTTTTTAAAGTTGTTTATCCTGTTGATGATCATCATTTTAAAGTTATTTTACCATATGGTACTTTAGTTATTGATGGTGTTACTCCTAATATGTTGAATTATTTTGGTCGTCCTTATGAAGGTATTGCTGTTTTTGATGGTAAAAAAATTACTGTTACTGGTACTTTATGGAATGGTAATAAAATTATTGATGAACGTTTAATTACACCTGATGGTTCTATGTTATTTCGTGTTACTATTAATTCTTAAGATATC
SmBiT-cMycGATATCAGACCAGCATGTAAAATTCCTAATGATTTAAAACAAAAAGTTATGAATCATATGGAACAAAAGTTGATTTCTGAAGAAGATTTGGTTACTGGTTATCGTTTATTTGAAGAAATTTTATAAAAGCTT

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files. The dataset for the metabolite profiling is provided in Source data 1.

References

  1. Book
    1. Omini J
    (2024)
    Doctoral Dissertation: Role of the Malate Dehydrogenase-Citrate Synthase Metabolon in Metabolic Flux Regulation
    University of Nebraska-Lincoln.
    1. Robinson JB
    2. Inman L
    3. Sumegi B
    4. Srere PA
    (1987)
    Further characterization of the Krebs tricarboxylic acid cycle metabolon
    The Journal of Biological Chemistry 262:1786–1790.
    1. Srere PA
    (1985) The metabolon
    Trends in Biochemical Sciences 10:109–110.
    https://doi.org/10.1016/0968-0004(85)90266-X
  2. Book
    1. Thierie J
    2. Penninckx M
    (2010) Crabtree effect
    In: Flickinger MC, editors. Encyclopedia of Industrial Biotechnology. John Wiley & Sons, Ltd. pp. 1–384.
    https://doi.org/10.1002/9780470054581.eib243

Article and author information

Author details

  1. Joy Omini

    Department of Biochemistry and Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, United States
    Contribution
    Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
  2. Inga Krassovskaya

    Department of Biochemistry and Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, United States
    Contribution
    Investigation, Methodology, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
  3. Taiwo Adeolu Dele-Osibanjo

    Department of Biochemistry and Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, United States
    Contribution
    Validation, Investigation, Visualization, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2537-6998
  4. Connor Pedersen

    Department of Biochemistry and Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, United States
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  5. Toshihiro Obata

    Department of Biochemistry and Center for Plant Science Innovation, University of Nebraska-Lincoln, Lincoln, United States
    Contribution
    Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Visualization, Writing – original draft, Project administration, Writing – review and editing
    For correspondence
    tobata2@unl.edu
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8931-7722

Funding

U.S. National Science Foundation (1845451)

  • Toshihiro Obata

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

This study is supported by the NSF CAREER Award to TO under Grant No. 1845451. Dr. Liza Pon at Columbia University kindly gifted us the mito-Go Ateam2 and mito-ROGFP1 biosensors, and Dr. Anita Ayer at the University of New South Wales, Sydney, Australia, gifted us the pHluorin mitochondrial biosensor. Dr. Oleh Khalimonchuk at the University of Nebraska-Lincoln provided us with the anti-DLD1 antibody. The Microscopy work assisted by Terri Fangman was carried out at the Microscopy Research Core Facility of the Center for Biotechnology at the University of Nebraska-Lincoln, which is partially funded by the Nebraska Center for Integrated Biomolecular Communication COBRE grant (P20 GM113126 and NIGMS) and the Nebraska Research Initiative. This paper is based on a dissertation submitted by Joy Omini to fulfill the requirements for the degree of Doctor of Philosophy, University of Nebraska-Lincoln (Omini, 2024).

Version history

  1. Sent for peer review:
  2. Preprint posted:
  3. Reviewed Preprint version 1:
  4. Reviewed Preprint version 2:
  5. Version of Record published:

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.107953. This DOI represents all versions, and will always resolve to the latest one.

Copyright

© 2025, Omini et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 814
    views
  • 29
    downloads
  • 1
    citation

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Citations by DOI

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Joy Omini
  2. Inga Krassovskaya
  3. Taiwo Adeolu Dele-Osibanjo
  4. Connor Pedersen
  5. Toshihiro Obata
(2026)
Dynamic assembly of malate dehydrogenase–citrate synthase multienzyme complex in the mitochondria
eLife 14:RP107953.
https://doi.org/10.7554/eLife.107953.3

Share this article

https://doi.org/10.7554/eLife.107953