A coarsegrained NADH redox model enables inference of subcellular metabolic fluxes from fluorescence lifetime imaging
Abstract
Mitochondrial metabolism is of central importance to diverse aspects of cell and developmental biology. Defects in mitochondria are associated with many diseases, including cancer, neuropathology, and infertility. Our understanding of mitochondrial metabolism in situ and dysfunction in diseases are limited by the lack of techniques to measure mitochondrial metabolic fluxes with sufficient spatiotemporal resolution. Herein, we developed a new method to infer mitochondrial metabolic fluxes in living cells with subcellular resolution from fluorescence lifetime imaging of NADH. This result is based on the use of a generic coarsegrained NADH redox model. We tested the model in mouse oocytes and human tissue culture cells subject to a wide variety of perturbations by comparing predicted fluxes through the electron transport chain (ETC) to direct measurements of oxygen consumption rate. Interpreting the fluorescence lifetime imaging microscopy measurements of NADH using this model, we discovered a homeostasis of ETC flux in mouse oocytes: perturbations of nutrient supply and energy demand of the cell do not change ETC flux despite significantly impacting NADH metabolic state. Furthermore, we observed a subcellular spatial gradient of ETC flux in mouse oocytes and found that this gradient is primarily a result of a spatially heterogeneous mitochondrial proton leak. We concluded from these observations that ETC flux in mouse oocytes is not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration.
Editor's evaluation
This paper describes the derivation and validation of a coarsegrained model to measure mitochondrial metabolism at cellular and subcellular resolution by exploiting fluorescence lifetime imaging of NADH. This technique is applied to mouse oocytes subjected to a variety of metabolic stresses and to human tissue culture cells, revealing spatial gradients in mitochondrial NADH oxidation. This method represents an exciting new approach to quantifying mitochondrial electron transport chain rates and provides for the first time a method to study mitochondrial metabolic flux with subcellular resolution.
https://doi.org/10.7554/eLife.73808.sa0Introduction
Cells transduce energy from the environment to power cellular processes. Decades of extensive research have produced a remarkable body of detailed information about the biochemistry of mitochondrial energy metabolism (Salway, 2017). In brief, metabolites, such as pyruvate, are transported into mitochondria, where they are broken down and their products enter the tricarboxylic acid cycle (TCA). The TCA is composed of a number of chemical reactions, which ultimately reduces NAD^{+} to NADH. NADH and oxygen are then utilized by the electron transport chain (ETC) to pump hydrogen ions across the mitochondrial membrane. ATP synthase uses this proton gradient to power the synthesis of ATP from ADP (Mitchell, 1961). The activities of mitochondrial energy metabolism are characterized by the fluxes through these pathways: that is, the number of molecules turned over per unit time (Stephanopoulos, 1999). However, despite the wealth of knowledge concerning mitochondrial biochemistry, the spatiotemporal dynamics of cellular energy usage remains elusive and it is still unclear how cells partition energy across different cellular processes (Dumollard et al., 2007; Van Blerkom, 2011; Yellen, 2018; Yang et al., 2021) and how energy metabolism is misregulated in diseases (Brand and Nicholls, 2011; Lin and Beal, 2006; Wallace, 2012; Bratic and Larsson, 2013; Lowell and Shulman, 2005; Mick et al., 2020). Metabolic heterogeneities, between and within individual cells, are believed to be widespread, but remain poorly characterized (Takhaveev and Heinemann, 2018; Aryaman et al., 2018). Mitochondria have been observed to associate with the cytoskeleton (Lawrence et al., 2016), spindle (Wang et al., 2020), and endoplasmic reticulum (Dumollard et al., 2004) and display subcellular heterogeneities in mtDNA sequence (Morris et al., 2017) and mitochondrial membrane potential (Smiley et al., 1991). These observations suggest the potential existence of subcellular patterning of mitochondrial metabolic fluxes that could be critical in processes such as oocyte maturation (Yu et al., 2010) and embryo development (Sanchez et al., 2019). The limitations of current techniques for measuring mitochondrial metabolic fluxes with sufficient spatiotemporal resolution present a major challenge. In particular, there is a lack of techniques to measure mitochondrial metabolic fluxes with single cell and subcellular resolution.
Bulk biochemical techniques for measuring metabolic fluxes, such as oxygen consumption and nutrient uptake rates (Ferrick et al., 2008; Houghton et al., 1996; Lopes et al., 2005), and isotope tracing by mass spectrometry (Wiechert, 2001), require averaging over large populations of cells. Such techniques cannot resolve cellular, or subcellular, metabolic heterogeneity (Takhaveev and Heinemann, 2018; Aryaman et al., 2018). Biochemical approaches for measuring mitochondrial metabolic fluxes, such as mass spectrometry, are also often destructive (Wiechert, 2001; Saks et al., 1998), and thus cannot be used to observe continual changes in fluxes over time. Fluorescence microscopy provides a powerful means to measure cellular and subcellular metabolic heterogeneity continuously and nondestructively, with high spatiotemporal resolution. However, while fluorescent probes can be used to measure mitochondrial membrane potential (Perry et al., 2011) and the concentration of key metabolites (Imamura et al., 2009; Berg et al., 2009; DíazGarcía et al., 2017; San Martín et al., 2014), it is not clear how to relate those observables to mitochondrial metabolic fluxes.
NADH is an important cofactor that is involved in many metabolic pathways, including the TCA and ETC in mitochondria. NADH binds with enzymes and acts as an electron carrier that facilitates redox reactions. In the ETC, for example, NADH binds to complex I and donates its electron to ubiquinone and ultimately to oxygen, becoming oxidized to NAD^{+}. Endogenous NADH has long been used to noninvasively probe cellular metabolism because NADH is autofluorescent, while NAD^{+} is not (Heikal, 2010). Fluorescence lifetime imaging microscopy (FLIM) of NADH autofluorescence allows quantitative measurements of the concentration of NADH, the fluorescence lifetimes of NADH, and the fraction of NADH molecules bound to enzymes (Becker, 2012; Becker, 2019; Bird et al., 2005; Skala et al., 2007; Heikal, 2010; Sharick et al., 2018; Sanchez et al., 2018; Sanchez et al., 2019; Ma et al., 2019). It has been observed that the fraction of enzymebound NADH and NADH fluorescence lifetimes are correlated with the activity of oxidative phosphorylation, indicating that there is a connection between NADH enzymebinding and mitochondrial metabolic fluxes (Bird et al., 2005; Skala et al., 2007). The mechanistic basis of this empirical correlation has been unclear.
Here, we developed a generic coarsegrained NADH redox model that enables the inference of ETC flux with subcellular resolution from FLIM measurements. We validated this model in mouse oocytes and human tissue culture cells subject to a wide range of perturbations by comparing predicted ETC fluxes from FLIM to direct measurements of oxygen consumption rate (OCR), and by a selfconsistency criterion. Using this method, we discovered that perturbing nutrient supply and energy demand significantly impacts NADH metabolic state but does not change ETC flux. We also discovered a subcellular spatial gradient of ETC flux in mouse oocytes and found that this flux gradient is primarily due to a spatially heterogeneous mitochondrial proton leak. We concluded from these observations that ETC flux in mouse oocytes is not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration. Thus, FLIM of NADH can be used to noninvasively and continuously measure mitochondrial ETC fluxes with subcellular resolution and provides novel insights into the spatiotemporal regulation of metabolic fluxes in cells.
Results
Quantifying response of mitochondrial metabolism to changing oxygen levels and metabolic inhibitors using FLIM of NADH
We used meiosis II arrested mouse oocytes as a model system. MII oocytes are in a metabolic steadystate, which eases interpretations of metabolic perturbations. ATP synthesis in mouse oocytes occurs primarily through oxidative phosphorylation using pyruvate, without an appreciable contribution from glycolysis (Houghton et al., 1996), providing an excellent system to study mitochondrial metabolism. Mouse oocytes can be cultured in vitro using chemically welldefined media (Biggers and Racowsky, 2002). In our work, we used AKSOM as the culturing media (Summers, 2013). The oocytes can directly take up pyruvate supplied to them or derive it from lactate through the activity of lactate dehydrogenase (LDH) (Lane and Gardner, 2000), and they can remain in a steadystate for hours with constant metabolic fluxes. While NADH and NADPH are difficult to distinguish with fluorescence measurements due to their overlapping fluorescence spectrum, the concentration of NADH in mouse oocytes is 40 times greater than the concentration of NADPH for the whole cell (Bustamante et al., 2017) and potentially even greater in mitochondria (Zhao et al., 2011), so the autofluorescence signal from these cells (particularly from mitochondria) can be safely assumed to result from NADH.
To investigate how FLIM measurements vary with mitochondrial activities, we performed quantitative metabolic perturbations. We first continually varied the concentration of oxygen in the media, from 50±2 μM to 0.26±0.04 μM, over a course of 30 min while imaging NADH autofluorescence of oocytes with FLIM (Figure 1a, top, black curve; Video 1). NADH is present in both mitochondria and cytoplasm where it is involved in different metabolic pathways. To specifically study the response of NADH in mitochondria, we used a machine learningbased algorithm to segment mitochondria from the NADH intensity images (Berg et al., 2019; Figure 1b and Figure 1—figure supplement 1). We verified the accuracy of the segmentation with a mitochondrial labeling dye, MitoTracker Red FM, which showed a 80.6±1% (SEM) accuracy of the segmentation (Appendix 1).
Using the segmentation mask, we obtained the intensity of NADH, $I$, in mitochondria by averaging the photon count over all mitochondrial pixels. The intensity increased with decreasing oxygen concentration (Figure 1a, top, red), as is readily seen from the raw images (Figure 1a, middle). Restoring oxygen to its original level caused a recovery of NADH intensity, indicating that the observed changes are reversible (Figure 1a; Video 1). These observations are consistent with the expectation that NADH concentration will increase at low oxygen levels due to oxygen’s role as the electron acceptor in the ETC. In addition to intensity, FLIM can be used to determine the enzyme engagement of NADH by measuring the photon arrival times, from which fluorescence lifetimes can be fitted. Enzymebound NADH has a much longer fluorescence lifetime than free NADH (Sharick et al., 2018), allowing bound and free NADH to be separately resolved, but the precise fluorescence lifetimes of NADH depend on a range of factors, including viscosity, pH, and the identity of the enzyme NADH binds to (Sharick et al., 2018; Ghukasyan and Heikal, 2015). To fit NADH fluorescence lifetimes, we grouped all detected photons from mitochondria to form histograms of photon arrival times from NADH autofluorescence for each time point (Figure 1a, lower). We fitted the histograms using a model in which the NADH fluorescence decay, $F(\tau )$, is described by the sum of two exponentials,
where ${\tau}_{\mathrm{l}}$ and ${\tau}_{\mathrm{s}}$ are long and short fluorescence lifetimes, corresponding to enzymebound NADH and free NADH, respectively, and $f$ is the fraction of enzymebound NADH (Sanchez et al., 2018; Sanchez et al., 2019) (Materials and methods).
We repeated the oxygen drop experiments for a total of 68 oocytes. Since the oxygen drop is much slower than the NADH redox reactions (30 min compared to a timescale of seconds), the oxygen perturbation can be safely assumed to be quasistatic, allowing the FLIM measurements to be determined as a function of oxygen levels. We averaged data from all oocytes to obtain a total of four FLIM parameters: mitochondrial NADH intensity, I, long and short fluorescence lifetimes, ${\tau}_{\mathrm{l}}$ and ${\tau}_{\mathrm{s}}$, and the fraction of enzymebound NADH, $f$. We determined how these parameters varied with oxygen level (Figure 1a and c). All parameters are insensitive to oxygen level until oxygen drops below ~10 μM. This observation is consistent with previous studies that showed mitochondria have a very high apparent affinity for oxygen (Chance and Williams, 1955; Gnaiger et al., 1998).
We next explored the relationship between the measured FLIM parameters and the concentration of NADH. Since bound and free NADH have different fluorescence lifetimes, and hence different molecular brightnesses, the NADH concentration is not generally proportional to NADH intensity. Assuming molecular brightness is proportional to fluorescence lifetime (Lakowicz, 2006), we derived a relation between NADH intensity, fluorescence lifetimes, and concentrations as
where ${c}_{\mathrm{s}}$ is a calibration factor that relates intensities and concentrations (see Appendix 1). We measured the calibration factor by titrating free NADH in vitro and acquiring FLIM data (Figure 1—figure supplement 2, Equation S4). To test the validity of this approach, we used Equations 2a,b to measure concentrations of free and bound NADH in solutions with different concentrations of purified LDH, an enzyme to which NADH can bind. The measured NADH bound concentration increases with LDH concentration while the sum of free and bound NADH concentration remains a constant and equal to the amount of NADH added to the solution (Figure 1—figure supplement 3). This result demonstrates that Equations 2a,b can be used to measure free and bound NADH concentrations from NADH intensity and lifetimes. It is well established that FLIM can be used to distinguish bound and free NADH in vivo based on the large change of fluorescence lifetime when NADH binds to enzymes (Skala et al., 2007; Heikal, 2010). Even though the exact amount that the lifetime changes depend on the specific enzyme NADH binds to (Sharick et al., 2018), enzymebound NADH always has a much longer fluorescence lifetime than free NADH. Therefore, the method to calculate free and bound concentrations of NADH from FLIM measurements is expected to hold in vivo. We next used this method to study NADH in mitochondria in oocytes. We applied Equations 2a,b to our FLIM data from oocytes and determined how the concentrations of free NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, and enzymebound NADH, $\left[\mathrm{N}\mathrm{A}\mathrm{D}{\mathrm{H}}_{\mathrm{b}}\right]$, depended on oxygen level (Figure 1c, lower right). Interestingly, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ increased as oxygen fell below ~10 μM, while $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$ did not vary with oxygen level.
We next explored the impact of metabolic inhibitors on mitochondrial NADH. We first inhibited LDH by adding 9 mM of oxamate to the AKSOM media. This led to a decrease of NADH intensity in the mitochondria (Figure 2a, upper) and significant changes in all FLIM parameters (Figure 2a, lower, p<0.001). We next inhibited complex I of the ETC using 5 μM of rotenone (in the presence of 9 mM of oxamate, to reduce cytoplasmic NADH signal for better mitochondrial segmentation). This resulted in a dramatic increase of NADH intensity in the mitochondria (Figure 2b, upper) and significant changes in NADH bound ratio and long lifetime (Figure 2b, lower, p<0.001). Then we inhibited ATP synthase with 5 μM of oligomycin (in the presence of 9 mM of oxamate), which, similar to rotenone, resulted in an increase of mitochondrial NADH intensity (Figure 2c, upper) and significant changes in all FLIM parameters (Figure 2c, lower, p<0.001). Finally, we subjected the oocytes to 5 μM of FCCP (in the presence of 9 mM of oxamate), which uncouples proton translocation from ATP synthesis, and observed a decrease of mitochondrial NADH intensity (Figure 2d, upper) and significant changes in FLIM parameters (Figure 2d, lower, p<0.001). Interestingly, the direction of change of FLIM parameters under FCCP is opposite to those under rotenone and oligomycin. For each of these conditions, we used Equations 2a,b to calculate the concentrations of free NADH, $\left[\mathrm{N}\mathrm{A}\mathrm{D}{\mathrm{H}}_{\mathrm{f}}\right]$, (Figure 2e) and bound NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$, (Figure 2f) from the measured intensity and FLIM parameters. While rotenone and oligomycin significantly increased $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and decreased $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$, FCCP decreased $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$. It remains unclear how to relate these changes of the free and bound concentrations of NADH to the activities of mitochondrial respiration.
Developing a coarsegrained NADH redox model to relate FLIM measurements of NADH to activities of mitochondrial metabolic pathways
We next developed a mathematical model of NADH redox reactions to relate these quantitative FLIM measurements to activities of mitochondrial metabolic pathways. NADH is a central coenzyme that binds to enzymes and facilitates redox reactions by acting as an electron carrier. There are two categories of enzymes associated with NADH redox reactions, which together form a redox cycle: oxidases that oxidize NADH to NAD^{+} and reductases that reduce NAD^{+} to NADH. The major NADH oxidase in mitochondria is complex I of ETC for mammalian cells. There are many NADH reductases in mitochondria because NADH can be reduced through different pathways depending on the energy substrate. These pathways include the TCA cycle, fatty acid catabolism via beta oxidation, amino acid catabolism such as glutaminolysis and the malateaspartate shuttle (Salway, 2017). A comprehensive NADH redox model will include all the oxidases and reductases. For generality, we consider N oxidases and M reductases.
For convenience, we introduced a reduced notation to describe models of the enzyme kinetics of these oxidases and reductases. We began by illustrating our reduced notation using reversible MichaelisMenten kinetics as an example (Keleti, 1986; Miller and Alberty, 2002; Smith, 1992). The conventional, full notation for these kinetics (Figure 3a, left) explicitly displays all chemical species that are modeled in this reaction scheme — free NADH, ${\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}$, free enzyme, ${\mathrm{O}\mathrm{x}}_{i}$, free NAD^{+}, ${\mathrm{N}\mathrm{A}\mathrm{D}}_{\mathrm{f}}^{+}$, and NADH bound to the enzyme — as well as the forward and reverse reaction rates — ${k}_{1}$, ${k}_{1}$, ${k}_{2}$, and ${k}_{2}$ Our reduced notation for reversible MichaelisMenten kinetics (Figure 3a, right) is an alternative way of representing the same mathematical model. In this reduced notation, only free NADH, free NAD^{+}, and NADH bound to the enzyme are explicitly shown, while the free enzyme concentration is only represented as entering through the effective binding rates ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{b}}$ and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{b}}$. The conventional, full notation, and the reduced notation are alternative ways of representing the same mathematical model, but the reduced notation is convenient to use in the derivation that follows (see Appendix 2).
We next introduced a generalized enzyme kinetics using our reduced notation (Figure 3b), which contains not only free NADH, free NAD^{+}, and NADH bound to the enzyme, but also NAD^{+} bound to the enzyme, and the reaction rates for oxidation and reduction of the bound coenzymes. In this reduced notation, all the binding and unbinding rates, and the reaction rates, can be functions of metabolite concentrations, protein concentrations, and other factors such as pH and membrane potential. As in the reversible MichaelisMenten kinetics example, these rates can depend on the concentration of the free enzyme itself. This dependency on free enzyme concentration can be nonlinear, as could occur if the enzyme oligomerizes. Furthermore, the rates may depend on the concentration of free NADH, free NAD^{+}, and the enzyme complexes. Thus, while the reduced notation for the generalized enzyme might appear to describe a firstorder reaction, it can actually be used to represent reactions of any order, with arbitrary, nonlinear dependencies on the concentration of the enzyme itself, as well as arbitrary, nonlinear dependencies on other factors. In order to model the dynamics of enzymes described by such generalized kinetics, it is necessary to specify the functional form of all the rates, as well as specify mathematical models for all the variables that enter these rates (i.e., free enzyme concentration, membrane potential, pH, etc.) (Appendix 2). However, in what follows, we will derive results that hold true, irrespective of the functional form of the rates or the presence of additional, implicit variables. Thus, remarkably, these quantitative predictions are valid for enzyme kinetics of any order, with arbitrary nonlinearities in the rates.
To begin our derivation, we note that under this generalized enzyme kinetics (Figure 3b), the net flux through the ith oxidase at steadystate is:
where $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet {\mathrm{O}\mathrm{x}}_{i}\right]$, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}}^{+}\bullet {\mathrm{O}\mathrm{x}}_{i}\right]$, and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ are the concentrations of the ith oxidasebound NADH, NAD^{+}, and free NADH, respectively. ${r}_{{\mathrm{o}\mathrm{x}}_{i}}^{+}$ and ${r}_{{\mathrm{o}\mathrm{x}}_{i}}^{}$ are the forward and reverse oxidation rates. ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{b}}$ and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{u}}$ are the binding and unbinding rates. The second equality in Equation 3 results from the steadystate condition, where the net binding and unbinding flux equals the net oxidation flux.
We next considered a redox cycle between NADH and NAD^{+} with multiple oxidases and reductases. To account for all possible NADH redox pathways, we developed a detailed NADH redox model with N oxidases and M reductases described by the generalized enzyme kinetics (Figure 4a and Figure 4—figure supplement 1). In this model, NADH and NAD^{+} can bind and unbind to each oxidase and reductase. Once bound, NADH can be reversibly oxidized to NAD^{+} by the oxidases, and NAD^{+} can be reversibly reduced to NADH by the reductases, forming a redox cycle. The functional dependencies of the binding and unbinding rates, and the reaction rates, can be different for each oxidase and reductase, and each of these rates can be nonlinear functions of free enzyme concentrations, NADH concentration, and other factors such as pH and membrane potential. Modeling the dynamics of this redox cycle requires specifying the precise number of oxidases and reductases, the functional forms of the rates, and mathematical models for all the variables the rates implicitly depend on. However, we will show that quantitative predictions regarding the interpretation of FLIM measurements can be made that generally hold, independent of these modeling choices.
FLIM cannot resolve the association of NADH with individual enzymes in cells, but rather, provides quantitative information on the global states of bound and free NADH. Thus, to facilitate comparison to FLIM experiments, we coarsegrained the detailed redox model by mapping all N oxidases into a single effective oxidase and all M reductases into a single effective reductase (Figure 4b and Appendix 3). This coarsegraining is mathematically exact and involves no approximations or assumptions.
In the coarsegrained redox model, NADH can be bound to the effective oxidase, $\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}$, bound to the effective reductase, $\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{R}\mathrm{e}$, or can be free, ${\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}$. Hence, the concentration of NADH bound to all enzymes is, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]=\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}\right]+[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{R}\mathrm{e}]$, and the total concentration of NADH is, $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\right]=\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]+\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$. The kinetics of the effective oxidase and reductase are represented by the coarsegrained forward, ${r}_{\mathrm{o}\mathrm{x}}^{+}$, and reverse, ${r}_{\mathrm{o}\mathrm{x}}^{}$, oxidation rates, and the forward, ${r}_{\mathrm{r}\mathrm{e}}^{+}$, and reverse, ${r}_{\mathrm{r}\mathrm{e}}^{}$, reduction rates. The global flux through all the oxidases in the detailed redox model equals the global flux through the coarsegrained oxidase, which at steadystate is:
where ${k}_{\mathrm{o}\mathrm{x}}^{\mathrm{b}}$ is the rate that free NADH binds the effective oxidase, ${k}_{\mathrm{o}\mathrm{x}}^{\mathrm{u}}$ is the rate that NADH unbinds the effective oxidase, and the last equality results because the coarsegrained redox loop is a linear pathway so the global oxidative flux must equal the global binding and unbinding flux at steadystate. The conservation of global flux explicitly relates the effective binding and unbinding rates and the reaction rates of the coarsegrained model to those of the detailed model (Appendix 3, Figure 4—figure supplement 1). The binding and unbinding kinetics of NADH and NAD^{+} to the effective oxidase and reductase are described by eight coarsegrained binding and unbinding rates (Figure 4b). The coarsegrained reaction rates and binding and unbinding rates can be arbitrary functions of metabolite concentrations, enzyme concentrations, and other factors (i.e., pH, membrane potential, etc.). These effective rates can even be functions of $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}}_{\mathrm{f}}^{+}\right]$, and the concentration of other variables, and thus can include reactions of arbitrary order. Hence, this coarsegrained model is a generic model of NADH redox reactions. Fully specifying this model would require explicitly choosing the functional form of all the rates and incorporating additional equations to describe the dynamics of all the implicit variables that the rates depend on (Appendix 2). We next demonstrate that quantitative predictions regarding the interpretation of FLIM measurements of NADH can be made that are valid irrespective of the form of the rates or the presence of implicit variables.
Accurately predicting ETC flux from FLIM of NADH using the NADH redox model
At steadystate, the model can be further coarsegrained, without approximation, to consider only free NADH, with a turnover rate of $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, and free NAD^{+}, with a turnover rate of $\stackrel{~}{r}}_{\mathrm{r}\mathrm{e}$ (Figure 4c). Our key prediction is that the steadystate global oxidative flux of NADH is (Appendix 4):
where
and
This prediction results from the steadystate assumption where the net binding and unbinding flux of NADH from the oxidase balances the net oxidative flux through the oxidase (Equation 4 and Appendix 4). The turnover rate of free NADH, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, is proportional to the difference between the NADH bound ratio $\beta $, that is, the ratio between bound and free NADH concentrations, and the equilibrium NADH bound ratio, ${\beta}_{\mathrm{e}\mathrm{q}}$ (i.e., what the bound ratio would be if the global oxidative flux is zero). $\beta}_{\mathrm{e}\mathrm{q}$ and the prefactor $\alpha$ are independent of the reaction rates of the oxidase and reductase and can be explicitly related to the binding and unbinding rates of the coarsegrained model (Appendix 4, Equation S43 and S45).
In mitochondria, the major NADH oxidation pathway is the ETC. Thus, Equations 5ac predict that there is a direct connection between quantities that can be measured by FLIM of NADH in mitochondria (i.e., $\beta $ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$) and the flux through the ETC (i.e., ${J}_{\mathrm{o}\mathrm{x}}$). Equations 5ac suggest a procedure for using FLIM to infer flux through the ETC: if a condition can be found under which there is no net flux through the ETC, then ${\beta}_{\mathrm{e}\mathrm{q}}$ can be measured with FLIM. Once ${\beta}_{\mathrm{e}\mathrm{q}}$ is known, then subsequent FLIM measurements of $\beta $ allows $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, and hence ${J}_{\mathrm{o}\mathrm{x}}$, to be inferred (up to a constant of proportionality $\alpha$) (Appendix 5).
Equations 5ac are valid irrespective of the functional forms of the rate laws, which may have nonlinear dependencies on metabolite concentrations, enzyme concentrations, and other factors. While Equation 5a seems to imply firstorder kinetics in $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, the rates can also be arbitrary functions of $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, so Equations 5a–c hold for kinetics of any order. Equations 5ac are also applicable if the rates depend on additional variables that have their own dynamical equations (as long as the system is at steadystate): as an example, Appendix 7 shows that Equations 5ac result when the N oxidases and M reductases are each described by reversible MichaelisMenten kinetics, a model in which the rates depend on the concentration of free enzymes (which is a dynamical variable). More generally, if detailed biophysical models of the NADH oxidases are available, then parameters of these models can be explicitly mapped to the coarsegrained parameters of the NADH redox model. Appendix 7 and Appendix 7—table 1 contain mappings between the coarsegrained model and a number of previously proposed detailed biophysical models of NADH oxidation in the ETC (Beard, 2005; Korzeniewski and Zoladz, 2001; Hill, 1977; Jin and Bethke, 2002; Chang et al., 2011). However, since Equations 5ac are valid for a broad set of models, they can be used for flux inference without the need to specify the functional form of the rates or the variables they depend on. This is because the rates are coarsegrained into two effective parameters $\alpha $ and ${\beta}_{\mathrm{e}\mathrm{q}}$, which can be experimentally determined. This generality results from the steadystate assumption and the topology of the reactions resulting in the net binding and unbinding flux of NADH from the oxidase balancing the net oxidative flux. Thus, Equations 5a–c provide a general procedure to infer the ETC flux from FLIM measurements of NADH in mitochondria.
We applied this procedure to analyze our oxygen drop experiments (Figure 1) by assuming that there was no net flux through the ETC at the lowest oxygen level achieved for each oocyte (implying that the measured value of $\beta $ at that oxygen concentration corresponds to ${\beta}_{\mathrm{e}\mathrm{q}}$ for that oocyte). We also assumed that $\alpha $ and ${\beta}_{\mathrm{e}\mathrm{q}}$ do not change with oxygen levels, which is reasonable since, as noted above, they are independent of the reaction rates of the oxidase and reductase. The measured value of ${\beta}_{\mathrm{e}\mathrm{q}}$ allowed us to obtain a prediction for ${J}_{\mathrm{o}\mathrm{x}}$ as a function of oxygen concentration for the oocytes (Figure 5a). To test these predictions, we directly determined ${J}_{\mathrm{o}\mathrm{x}}$ as a function of oxygen concentration by measuring the OCR of the oocytes using a nanorespirometer (Lopes et al., 2005) (Materials and methods). The direct measurements of ${J}_{\mathrm{o}\mathrm{x}}$ from OCR quantitatively agree with the predictions of ${J}_{\mathrm{o}\mathrm{x}}$ from FLIM for all oxygen concentrations (Figure 5a), strongly arguing for the validity of the model and the inference procedure. This agreement supports the assumption that $\alpha $ and ${\beta}_{\mathrm{e}\mathrm{q}}$ are independent of oxygen levels.
So far, we have inferred the ETC flux up to a constant of proportionality $\alpha$, allowing the relative changes of ETC flux to be inferred from FLIM of NADH. $\alpha$ cannot be determined by FLIM alone. If an absolute measurement of the ETC flux can be obtained at one condition, then $\alpha$ can be calibrated to predict absolute ETC fluxes for other conditions. OCR measurement provides a means to calibrate $\alpha$ (Appendix 5, Equation S48). We used oocytes cultured in AKSOM media at 50±2 μM oxygen as a reference state, which, from our OCR measurements yielded ${J}_{\mathrm{o}\mathrm{x}}=56.6\pm 2\text{}\mu \mathrm{M}/s$ (SEM) and hence a constant of proportionality of $\alpha =5.4\pm 0.2{\mathrm{}\mathrm{s}}^{1}$. Using this value of $\alpha $, we can predict absolute values of ${J}_{\mathrm{o}\mathrm{x}}$ under various perturbations assuming $\alpha $ remains a constant. We note that ${J}_{\mathrm{o}\mathrm{x}}$ is a flux density with units of concentration per second, an intensive quantity that does not depend on the mitochondrial volume. Multiplying ${J}_{\mathrm{o}\mathrm{x}}$ by the volume of mitochondria in an oocyte gives the total ETC flux, proportional to OCR, in that oocyte. In all subsequent discussions, ETC flux refers to flux density unless otherwise noted.
We next applied the inference procedure and a constant of $\alpha =5.4\pm 0.2{\mathrm{}\mathrm{s}}^{1}$ to analyze the experiments of oxamate, FCCP, rotenone and oligomycin perturbations (Figure 2). We dropped oxygen levels to determine ${\beta}_{\mathrm{e}\mathrm{q}}$ in the presence of oxamate (Figure 5—figure supplement 1h) and applied Equations 5ac to infer the impact of oxamate on ${J}_{\mathrm{o}\mathrm{x}}$ at 50 μM oxygen (i.e., control levels of oxygen). Surprisingly, while the addition of oxamate greatly impacts FLIM parameters, including a 29±2% (SEM) decrease in intensity and a 10±3% increase in bound ratio (Figure 2a), this procedure revealed that the predicted ETC flux with oxamate (${J}_{\mathrm{o}\mathrm{x}}=55.2\pm 3.2\text{}\mu \mathrm{M}/s$) is the same as that without oxamate (${J}_{\mathrm{o}\mathrm{x}}=55.4\pm 1.9\text{}\mu \mathrm{M}/s$) (Figure 5b; p=0.95), which was confirmed by direct measurements of oocytes’ OCR that yielded ${J}_{\mathrm{o}\mathrm{x}}=55.4\pm 1.5\text{}\mu \mathrm{M}/s$ and ${J}_{\mathrm{o}\mathrm{x}}=54.9\pm 0.7\text{}\mu \mathrm{M}/s$ before and after the addition of oxamate, respectively (Figure 5b; p=0.85). We next analyzed the FCCP experiment. We obtained ${\beta}_{\mathrm{e}\mathrm{q}}$ by dropping oxygen in the presence of FCCP (Figure 5—figure supplement 1h) and applied Equations 5ac to infer the impact of FCCP on ${J}_{\mathrm{o}\mathrm{x}}$ at 50 μM oxygen. We predicted that FCCP increased the flux to ${J}_{\mathrm{o}\mathrm{x}}=67.7\pm 1.5\text{}\mu \mathrm{M}/s$, which was confirmed by the directly measured ${J}_{\mathrm{o}\mathrm{x}}=74.0\pm 3.7\text{}\mu \mathrm{M}/s$ from OCR (Figure 5b; p=0.30). Following the same FLIM based inference procedures, we predicted that the addition of rotenone and oligomycin reduced the fluxes to ${J}_{\mathrm{o}\mathrm{x}}=16.7\pm 1.4\text{}\mu \mathrm{M}/s$ and ${J}_{\mathrm{o}\mathrm{x}}=15.7\pm 1.3\text{}\mu \mathrm{M}/s$, respectively, which was again confirmed by corresponding direct measurements of OCR that yielded ${J}_{\mathrm{o}\mathrm{x}}=11.1\pm 0.4\text{}\mu \mathrm{M}/s$ and ${J}_{\mathrm{o}\mathrm{x}}=22.3\pm 0.6\text{}\mu \mathrm{M}/s$ (Figure 5b; p=0.31 and p=0.17). The quantitative agreement between predicted fluxes from FLIM and directly measured fluxes from OCR under a variety of conditions (i.e., varying oxygen tension, sodium oxamate, FCCP, rotenone, and oligomycin) demonstrates that Equations 5ac can be successfully used to infer flux through the ETC in mouse oocytes. This agreement also supports the assumption that $\alpha $ is a constant across these different perturbations.
The work described above used the relation ${\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}}=\alpha \left(\beta {\beta}_{\mathrm{e}\mathrm{q}}\right)$ to predict the flux through the ETC from FLIM measurements. We next show that the model also predicts a relationship between $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ and the fluorescence lifetime of enzymebound NADH, ${\tau}_{\mathrm{l}}$, in mitochondria. This provides a second means to use the model to infer $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, and hence ${J}_{\mathrm{o}\mathrm{x}}$, from FLIM of NADH. Specifically, we assumed that NADH bound to the oxidases have a different average lifetime, ${\tau}_{\mathrm{o}\mathrm{x}}$, than NADH bound to the reductases, ${\tau}_{\mathrm{r}\mathrm{e}}$, which is reasonable because NADH bound to different enzymes do exhibit different fluorescence lifetimes (Sharick et al., 2018). This assumption implies that the experimentally measured long lifetime of NADH in mitochondria, ${\tau}_{\mathrm{l}}$, is a weighted sum of these two lifetimes,
Using the coarsegrained NADH redox model at steadystate, Equation 6 leads to a nontrivial prediction that ${\tau}_{\mathrm{l}}$ is linearly related to $1/\beta $ (Appendix 5):
where the slope $A$ and offset $B$ can be explicitly related to ${\tau}_{\mathrm{o}\mathrm{x}}$, ${\tau}_{\mathrm{r}\mathrm{e}}$, and the coarsegrained binding and unbinding rates. Such a linear relationship is indeed observed in individual oocytes subject to oxygen drops (Figure 6a and Figure 6—figure supplement 1), supporting the assumptions of the model. Combining Equations 7 and 5b leads to a predicted relationship between $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ and NADH long fluorescence lifetime (Appendix 5):
where ${\tau}_{\mathrm{e}\mathrm{q}}$ is the equilibrium NADH long fluorescence lifetime, that is, the value of the long lifetime when the global oxidative flux is zero. This provides a second means to infer $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ from FLIM measurements: dropping oxygen and plotting the relationship between ${\tau}_{\mathrm{l}}$ and $1/\beta $ provides a means to measure A and B from Equation 7, while ${\tau}_{\mathrm{e}\mathrm{q}}$ can be obtained from the NADH long fluorescence lifetime obtained at the lowest oxygen level. Once A, B, and ${\tau}_{\mathrm{e}\mathrm{q}}$ are known, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ can be inferred solely from NADH long fluorescence lifetime ${\tau}_{\mathrm{l}}$, using Equation 8.
We next used the lifetime method (Equation 8) and the bound ratio method (Equation 5b) to separately infer $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ in oocytes subject to a wide variety of conditions (varying oxygen levels, with oxamate, FCCP, rotenone, and oligomycin). We obtained A, B, ${\beta}_{\mathrm{e}\mathrm{q}}$, and ${\tau}_{\mathrm{e}\mathrm{q}}$ for these different conditions (Figure 5—figure supplement 1 and Figure 6—figure supplement 1), and used the two different methods to provide two independent measures of $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ (assuming $\alpha $ is constant across all conditions). The predictions of $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ from these two methods quantitatively agree under all conditions (Figure 6b, p=0.73), which is a strong selfconsistency check that further supports the use of the model to infer ETC flux from FLIM measurements of NADH.
The NADH redox model enables accurate prediction of ETC flux in human tissue culture cells
After thoroughly testing the NADH redox model and the inference procedure in mouse oocytes, we next investigated if it can be used in other cell types. We chose human tissue culture cells for this purpose, since they are widely used as model systems to study metabolic dysfunctions in human diseases including cancer (Vander Heiden et al., 2009) and neuropathology (Lin and Beal, 2006).
While mouse oocytes have a negligible level of NADPH compared to NADH (Bustamante et al., 2017), the concentrations of NADH and NADPH are similar in tissue culture cells (10–100 μM averaged over the whole cell) (Lu et al., 2018; Park et al., 2016; Blacker et al., 2014). Since NADPH and NADH have overlapping fluorescent spectra (Patterson et al., 2000), the presence of NADPH may complicate the interpretation of FLIM experiments. Thus, we investigated the impact of background fluorescence, such as from NADPH, on the flux inference procedure. If the background fluorescence does not change with the perturbations under study, then it can be treated as an additive offset that systematically makes the measured concentrations of free and bound NADH different from their actual values. In this case, a derivation in Appendix 5 demonstrates that the background fluorescence can be incorporated into the equilibrium bound ratio ${\beta}_{\mathrm{e}\mathrm{q}}$ and does not impact the flux inference procedures. In other words, if the modified ${\beta}_{\mathrm{e}\mathrm{q}}$ can be reliably determined, then the measured concentrations of free and bound fluorescent species can be used in place of the true values of NADH in Equations 5ac to infer the ETC flux. An alternative possibility is that the background fluorescence does change with the perturbations under study, but in a manner that is proportional to the change in NADH. In this case, the background fluorescence can be incorporated into the equilibrium bound ratio $\alpha $ and, once more, does not impact the flux inference procedures (Appendix 5). If the background fluorescence changes in some more complicated manner, then the inference procedure may no longer be valid. Thus, depending on the behavior of NADPH, it either might or might not interfere with the inference procedure: no impact if NADPH is either constant or proportional to changes in NADH, a possible impact otherwise. Therefore, the validity of the inference procedure in the presence of significant NADPH fluorescence must be established empirically.
We next tested the inference procedures experimentally in hTERTRPE1 (hTERTimmortalized retinal pigment epithelial cell line) tissue culture cells. We started by exploring the impact of metabolic perturbations on mitochondrial NAD(P)H: the combined signal from NADH and NADPH (which are indistinguishable) from mitochondria. We first cultured the cells in DMEM with 10 mM galactose (Materials and methods). We then inhibited complex I of the ETC by adding 8 μM of rotenone to the media. This resulted in a significant increase of mitochondrial NAD(P)H intensity (Figure 7a, upper). We segmented mitochondria using a machine learningbased algorithm from the intensity images of NAD(P)H, and fitted the fluorescence decay curves of mitochondrial NAD(P)H to obtain changes in FLIM parameters (Materials and methods). All FLIM parameters displayed significant changes (Figure 7a, lower, p<0.001, and Figure 7—figure supplement 1). We next uncoupled proton translocation from ATP synthesis by adding 3.5 μM CCCP to the media. This led to a decrease of NAD(P)H intensity in the mitochondria (Figure 7b, upper) and significant changes in NAD(P)H bound ratio and short lifetime, but in opposite directions as compared to rotenone perturbation (Figure 7b, lower, p<0.01, and Figure 7—figure supplement 1). Finally, we perturbed the nutrient conditions by culturing the cells in DMEM with 10 mM glucose. FLIM imaging revealed an increase of mitochondrial NAD(P)H intensity (Figure 7c, upper) and significant changes in all FLIM parameters as compared to the galactose condition (Figure 7c, lower, p<0.001, and Figure 7—figure supplement 1).
Inference of the ETC flux from FLIM measurements requires a measurement of ${\beta}_{\mathrm{e}\mathrm{q}}$. Since rotenone is known to drastically decrease the OCR of hTERTRPE1 cells to near zero (MacVicar and Lane, 2014), we used the NAD(P)H bound ratio measured in the presence of rotenone as ${\beta}_{\mathrm{e}\mathrm{q}}$. Different values of ${\beta}_{\mathrm{e}\mathrm{q}}$ were obtained for glucose and galactose conditions by adding 8 μM of rotenone to each condition (Figure 7—figure supplement 1). We next calculated the concentrations of free NAD(P)H, $\left[\mathrm{N}\mathrm{A}\mathrm{D}(\mathrm{P}){\mathrm{H}}_{\mathrm{f}}\right]$, from the FLIM parameters using Equation 2a. $\left[\mathrm{N}\mathrm{A}\mathrm{D}(\mathrm{P}){\mathrm{H}}_{\mathrm{f}}\right]$ displayed significant changes for all perturbations (Figure 7d). Using Equation 5b and assuming $\alpha$ is a constant, we calculated the NAD(P)H turnover rate, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, from the FLIM measurements and ${\beta}_{\mathrm{e}\mathrm{q}}$. $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ changed significantly for all perturbations (Figure 7e). Multiplying $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ and $\left[\mathrm{N}\mathrm{A}\mathrm{D}(\mathrm{P}){\mathrm{H}}_{\mathrm{f}}\right]$, we obtained the predicted ETC flux, ${J}_{\mathrm{o}\mathrm{x}}$, which increased under FCCP, decreased under glucose and reduced to zero under rotenone (Figure 7f).
To test the model predictions, we compared the predicted ETC flux with previous direct OCR measurements of the same cell type that we used, under the same conditions (MacVicar and Lane, 2014). Remarkably, the predicted changes in ETC fluxes are in quantitative agreement with the directly measured OCR across all conditions as estimated from Figure 1A of MacVicar and Lane, 2014: CCCP is predicted to increase the ETC flux by 14±3% (SEM), in agreement with the 18±21% increase from OCR measurement (p=0.80); Glucose is predicted to decrease ETC flux by 33±3%, in agreement with the 46±9% decrease from OCR measurements (p=0.30), shifting metabolism from oxidative phosphorylation to aerobic glycolysis. Since we used $\beta $ from rotenone treatment as ${\beta}_{\mathrm{e}\mathrm{q}}$, the predicted decrease in ETC flux after the addition of rotenone is 101±2%, which is in agreement with the 82±2% decrease from OCR measurement (p=0.28). This quantitative agreement between predicted ETC fluxes and measured OCR across all perturbations demonstrated the applicability of the NADH redox model and the flux inference procedures to tissue culture cells, even though they contain substantial levels of NADPH.
Homeostasis of ETC flux in mouse oocytes: perturbations of nutrient supply and energy demand impact NADH metabolic state but do not impact ETC flux
Having established the validity of the NADH redox model and the associated flux inference procedures, we next applied it to study energy metabolism in mouse oocytes. We began by investigating the processes that determine the ETC flux in MII mouse oocytes. Mitochondrialbased energy metabolism can be viewed as primarily consisting of three coupled cycles: the NADH/NAD^{+} redox cycle (which our NADH redox model describes), the proton pumping/dissipation cycle, and the ATP/ADP production/consumption cycle (Figure 8a). At the most upstream portion of this pathway, the reduction of NAD^{+} to NADH is powered by a supply of nutrients, while at the most downstream portion, energydemanding cellular processes hydrolyze ATP to ADP. To test whether nutrient supply and energy demand set ETC flux, we investigated the effect of perturbing these processes. To perturb supply, we first varied the concentration of pyruvate in the media from 181 μM (which is standard for AKSOM) to either 18.1 μM or 1.81 mM, and observed significant changes in NADH intensity and FLIM parameters (Figure 8b, left), demonstrating that the NADH metabolic state is altered. To perturb demand, we began by adding 10 μM nocodazole to the media, which disassembled the meiotic spindle, an energy user, and resulted in significant changes in NADH FLIM parameters (Figure 8b, center). Similarly, the addition of 10 μM latrunculin A disassembled the actin cortex and also produced significant changes in NADH FLIM parameters (Figure 8b, right).
We next performed additional perturbations of nutrient supply, inhibiting the conversion of lactate to pyruvate by LDH (with 9 mM oxamate) and inhibiting the malateaspartate shuttle (with 11 mM AOA). We performed additional perturbations of energy demand by inhibiting protein synthesis (with 1 mM cycloheximide) and ion homeostasis, by varying extracellular potassium concentrations from 0 mM to 15 mM, inhibiting the Na^{+}/K^{+} pump (with 2 mM ouabain), and adding an ionophore (10 μM gramicidin). All perturbations resulted in significant changes in NADH FLIM parameters (Figure 8—figure supplement 1), showing that NADH metabolic state is generally impacted by varying nutrient supply and cellular energy demand. We next used the NADH redox model and the measured FLIM parameters to infer the concentration and effective turnover rate of free NADH for these perturbations. The free NADH concentrations, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, and turnover rates, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, displayed large variations across the perturbations, ranging from 33.5±1 μM (SEM) to 56.0±2.8 μM and from 1.0±0.05 s^{1} to 1.65±0.09 s^{1}, respectively (Figure 8c). Surprisingly, the changes in $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ were highly anticorrelated such that the data points primarily fell within a region where the inferred ETC flux, ${J}_{\mathrm{o}\mathrm{x}}={\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}}\left[\mathrm{N}\mathrm{A}\mathrm{D}{\mathrm{H}}_{\mathrm{f}}\right]$, is a constant 55.5 μM·s^{1} (Figure 8c, solid line, shaded region indicates 5% error). Indeed, ANOVA tests confirmed that perturbing nutrient supplies and cellular energy demand lead to no significant change in either the inferred ETC flux (Figure 8d, p=0.20) or directly measured OCR (Figure 8e, p=0.07). Thus, while nutrient supply and cellular energy demand strongly affect mitochondrial NADH redox metabolism, they do not impact ETC flux. In contrast, ETC flux is impacted by perturbing proton leak and ATP synthesis (Figure 5). Taken together, this suggests that the ETC flux in mouse oocytes is set by the intrinsic properties of their mitochondria, which can adjust their NADH redox metabolism to maintain a constant flux when nutrient supplies and cellular energy demand are varied. The mechanistic basis of this homeostasis of ETC flux is unclear and will be an exciting topic for future research.
Subcellular spatial gradient of ETC flux in mouse oocytes: spatially inhomogeneous mitochondrial proton leak leads to a higher ETC flux in mitochondria closer to cell periphery
Our results presented so far were performed by averaging together FLIM measurements from all mitochondria within an oocyte. However, FLIM data is acquired with optical resolution, enabling detailed subcellular measurements. To see if there are spatial variations in FLIM measurements within individual oocytes, we computed the mean NADH fluorescence decay time for each mitochondrial pixel. The mean NADH fluorescence decay time displays a clear spatial gradient, with higher values closer to the oocyte center (Figure 9a).
To quantify this gradient in more detail, we partitioned mouse oocytes into equally spaced concentric regions (Figure 9b) and fitted the fluorescence decay curves from mitochondrial pixels within each region to obtain FLIM parameters as a function of distance from the oocyte center. NADH intensity, bound ratio, and long lifetime in mitochondria all display significant spatial gradient within oocytes (Figure 9c). Next, using Equations 5ac and ${\beta}_{\mathrm{e}\mathrm{q}}$ obtained at the lowest oxygen level, and confirming that ${\beta}_{\mathrm{e}\mathrm{q}}$ is uniform within the oocyte with complete inhibition of ETC using high concentration of rotenone (Figure 9—figure supplement 1), we predicted the ETC flux, $J}_{\mathrm{o}\mathrm{x}$, as a function of distance from the oocyte’s center. The ETC flux displayed a strong spatial gradient within oocytes, with a higher flux closer to the cell periphery (Figure 9d). Note that, as described above, ${J}_{\mathrm{o}\mathrm{x}}$ is actually a flux density with units of concentration per second. Thus, the measured flux gradient is not merely a reflection of variations in mitochondrial density, but instead indicates the existence of subcellular spatial heterogeneities in mitochondrial activities.
To investigate the origin of this flux gradient, we inhibited ATP synthase using 5 μM of oligomycin and repeated measurements of subcellular spatial variations in inferred fluxes. After inhibition, ${J}_{\mathrm{o}\mathrm{x}}$ decreased at all locations throughout the oocytes and displayed an even more dramatic flux gradient (Figure 9e). If oligomycin completely blocks ATP synthase, then the remaining flux must be the result of proton leak. If it is further assumed that proton leak remains the same with and without oligomycin, then the flux due to ATP synthase in control oocytes can be determined by subtracting the flux after oligomycin inhibition (i.e., the proton leak) from the flux before inhibition. Performing this procedure throughout oocytes indicates that proton leak greatly increases in mitochondria near the periphery of oocytes, where ATP production decreases (Figure 9f). This implies that the subcellular gradient in ETC flux is primarily caused by a gradient in proton leak and that mitochondria near the periphery of oocytes are less active in ATP production than those in the middle of the oocyte.
We hypothesized that a gradient in proton leak would result in a gradient of mitochondrial membrane potential, with lower membrane potential closer to the cell periphery where proton leak is the greatest. To test this, we measured mitochondrial membrane potential using the membrane potentialsensitive dye TMRM, which preferentially accumulates in mitochondria with higher membrane potential (AlZubaidi et al., 2019). We observed a strong spatial gradient of the intensity of TMRM in mitochondria within oocytes, with dimmer mitochondria near the cell periphery (Figure 9g and h), indicating that mitochondria near the periphery of the oocyte have a lower membrane potential. This result is robust to locally normalizing TMRM intensity by mitochondrial mass using a membrane potential insensitive dye (Mitotracker Red FM), or using an alternative membrane potentialsensitive dye, JC1 (Figure 9—figure supplement 2). The predicted flux of proton leak and mitochondrial TMRM intensity shows a strong negative correlation (Figure 9i), confirming our hypothesis.
Taken together, these results show that MII mouse oocytes contain subcellular spatial heterogeneities of mitochondrial metabolic activities. The observation that proton leak is responsible for the gradient of ETC flux suggests that the flux heterogeneity is a result of intrinsic mitochondrial heterogeneity. This is consistent with our conclusion from the homeostasis of ETC flux (Figure 8) that it is the intrinsic rates of mitochondrial respiration, not energy demand or supply, that controls the ETC flux. The causes and consequences of the subcellular spatial variation in mitochondrial activity remain unclear and are an exciting topic for future research.
Discussion
The NADH redox model is a general model to relate FLIM measurements of NADH to ETC fluxes
Despite extensive studies and applications of FLIM in metabolic research (Bird et al., 2005; Skala et al., 2007; Heikal, 2010; Sharick et al., 2018; Sanchez et al., 2018; Liu et al., 2018; Sanchez et al., 2019; Ma et al., 2019), it remains a challenge to relate FLIM measurements to the activities of the underlying metabolic pathways in cells. We overcame this challenge by developing a coarsegrained NADH redox model that leads to quantitative predictions for the relationship between FLIM measurements and the flux through the ETC. The model was constructed by explicitly coarsegraining a detailed NADH redox model with an arbitrary number of oxidases and reductases that represent all the possible enzymes involved in NADH redox reactions. The reactions in the detailed NADH redox model can be of arbitrary order and depend on implicit variables (i.e., free enzyme concentration, membrane potential, pH, etc.), which obey their own dynamical equations. The dynamics of the redox model will, of course, depend on the precise number of oxidases and reductases, the functional forms of the rates, and specific mathematical models for all the variables the rates implicitly depend on. However, the quantitative predictions relating FLIM measurements and ETC flux are independent of these modeling choices. Coarsegraining the detailed NADH redox model reduces all oxidases to an effective oxidase and all reductases to an effective reductase. The kinetic rates of the coarsegrained model can be related to those of the detailed model by keeping the global fluxes through the oxidases and the reductases the same in both models. The coarsegrained model predicts that the flux through the ETC is a product of the turnover rate and the concentration of free NADH (Equation 5a). The turnover rate is proportional to the difference between the nonequilibrium and the equilibrium NADH bound ratio (Equation 5b), which are measurable by FLIM of NADH (Equation 5c). Thus, this model provides a generic framework to relate FLIM measurements of NADH to the flux through the ETC in mitochondria.
The central assumption required for the validity of Equations 5ac is that the redox reactions, and binding and unbinding processes, can be approximated as being at steadystate (i.e., undergoing only quasistatic changes over perturbations or development). At steadystate, the net binding and unbinding flux balances the oxidative flux of NADH. Therefore, the measurement of binding and unbinding state of NADH from FLIM allows the inference of the ETC flux, irrespective of the detailed behaviors of the oxidative reactions.
Remarkably, all the binding and unbinding rates of the NADH redox model are coarsegrained into two effective parameters: $\alpha$ and ${\beta}_{\mathrm{e}\mathrm{q}}$, which can be experimentally measured. We determined the value of $\alpha$ from an OCR measurement (Appendix 5, Equation S48), and we determined the value of ${\beta}_{\mathrm{e}\mathrm{q}}$ from FLIM of NADH at low oxygen levels or from rotenone perturbation (Appendix 5, Figure 5—figure supplement 1h). In MII mouse oocytes, $\alpha$ do not significantly vary in response to oxygen or drug and nutrient perturbations. This is demonstrated by the agreement between the predicted ETC flux and the measured OCR with a constant $\alpha$ of 5.4±0.2 s^{1} across a variety of conditions (Figure 5). $\alpha$ is predicted to depend only on the coarsegrained unbinding rates of NADH from the enzymes (Equation S43), so the observed constancy of $\alpha$ implies that the perturbations in this study primarily impacted the reduction/oxidation reaction rates (and not the unbinding rates). In other scenarios, such as when the concentrations of enzymes change, the coarsegrained unbinding rates might change, so $\alpha$ might not be a constant. In contrast, ${\beta}_{\mathrm{e}\mathrm{q}}$ does vary with drug and nutrient perturbations, but not with oxygen level, allowing ${\beta}_{\mathrm{e}\mathrm{q}}$ to be obtained at the lowest oxygen level for different drug and nutrient conditions (Figure 5—figure supplement 1h and Figure 8—figure supplement 1). Using these two parameters, we inferred the effective turnover rate of free NADH, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, from FLIM measurements of NADH. By multiplying this turnover rate with the concentration of free NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ (also obtained from FLIM measurements using Equation 2a), we inferred the ETC flux from Equation 5a. Thus, all the complex behaviors of the binding and unbinding and reaction rates are captured by the variations in FLIM parameters of NADH, and our coarsegrained model provides a generic way to interpret these variations.
While we found that ${\beta}_{\mathrm{e}\mathrm{q}}$ is smaller than $\beta $ in mouse oocytes, this does not generically have to be true. Thus, if a perturbation is observed to decrease the NADH bound ratio $\beta ,$ it does not necessarily imply a decrease of the ETC flux. Similarly, a decrease of NADH long lifetime is not necessarily associated with an increase of the ETC flux. Therefore, measurements of $\alpha $ and ${\beta}_{\mathrm{e}\mathrm{q}}$ are required to use Equations 5ac to infer ETC flux from FLIM measurements of NADH.
The underlying assumptions and limitations of the NADH redox model
In this section, we clarify the underlying assumptions and limitations of the model to facilitate the accurate interpretation of FLIM measurements of NADH in different biological contexts.
To use the coarsegrained NADH model, segmentation needs to be performed to separate the mitochondrial NADH signal from the cytoplasmic NADH signal, because they encode different metabolic fluxes. In mouse oocytes, the segmentation can be reliably performed based on NADH images due to the higher NADH intensity in mitochondria than cytoplasm. Mitochondrial movements are also slow in MII oocytes (Video 1); hence, long exposure times can be used to obtain high contrast NADH images. For cells where NADH contrast is low, such as in yeast cells (Papagiannakis et al., 2017; Shaw and Nunnari, 2002), MitoTracker dye (Appendix 1 Figure 1—figure supplement 1) or mitochondrial associated fluorescent proteins (Westermann and Neupert, 2000) will likely be needed for reliable segmentation of mitochondria.
One of the most important assumptions that enables the coarsegrained model to be used to predict fluxes is that the NADH redox cycle can be well approximated as being at steadystate, that is, the rate of change of NADH concentrations is much slower than the kinetic rates, including the binding/unbinding rates and the reaction rates. This is true for mouse oocytes, where the NADH intensity does not significantly change over the course of hours. This assumption also holds for slow processes such as the cell cycle (Papagiannakis et al., 2017), which occurs on the timescale of hours compared to timescales of seconds for the kinetic rates. This claim is supported by the success of the model on human tissue culture cells. The steadystate approximation could fail for rapid dynamics of NADH, such as the transient overshoot of NADH in neurons induced by acute external stimulus (DíazGarcía et al., 2021), but this needs to be tested experimentally.
While NADH and NADPH share the same fluorescence spectrum, NADH concentration is 40 times greater than the concentration of NADPH for the whole mouse oocytes and presumably even higher for mitochondria (Bustamante et al., 2017). NADPH concentration can be comparable to that of NADH for other cell types such as tissue culture cells (Park et al., 2016). However, we have shown that the presence of NADPH signal and other background fluorescence signals only affect the equilibrium bound ratio ${\beta}_{\mathrm{e}\mathrm{q}}$ or the prefactor $\alpha$, and hence does not affect the flux inference procedure if ${\beta}_{\mathrm{e}\mathrm{q}}$ can be reliably determined and $\alpha$ remains a constant (Appendix 5). This was validated in tissue culture cells by comparing predicted ETC flux (Figure 7) with previous OCR measurements (MacVicar and Lane, 2014).
Finally, when relating NADH FLIM measurements to the ETC flux, we did not explicitly consider the contribution to the flux through FADH_{2}. This is a valid approximation when the FADH_{2} oxidative flux is much smaller than the NADH oxidative flux, as is often the case since pyruvate dehydrogenase plus the TCA cycle yields four NADH molecules but only one FADH_{2} molecule per cycle. Alternatively, if the FADH_{2} flux is proportional to the NADH flux, then a rescaled value of $\alpha$ can be used in Equation 5b to effectively account for both fluxes. The proportionality of FADH_{2} flux and NADH flux is expected when NADH and FADH_{2} are produced from the same redox cycle with fixed stoichiometry, such as the pyruvate dehydrogenase and TCA cycle. This proportionality will break down if significant amounts of NADH and FADH_{2} are produced in independent cycles where the stoichiometry varies, for example, when the glycerol phosphate shuttle acts as a reductase in mitochondria for FADH_{2} but not for NADH.
Given these underlying assumptions, the model needs to be tested before being applied to other biological systems. The present study provides an example for such tests in mouse oocytes and human tissue culture cells by comparing the predicted ETC flux from FLIM with direct measurements of OCR across a wide range of perturbations.
Towards spatiotemporal regulations of metabolic fluxes in cells
Cells transduce energy from nutrients to power various cellular processes. The ETC flux represents the total rate of energy conversion by mitochondria. Despite the detailed knowledge of the biochemistry of mitochondrial metabolism, it is still unclear what cellular processes determine ETC flux or how cells partition energetic fluxes to different cellular processes, including biosynthesis, ion pumping, and cytoskeleton assemblies. Energetic costs of specific cellular processes have been estimated from theoretical calculations (Stouthamer, 1973) or through inhibition experiments (Mookerjee et al., 2017). The latter typically involves measurements of the change of metabolic fluxes, such as OCR, upon inhibition of specific cellular processes, and interpreting this change as the energetic cost of the inhibited process. This interpretation is valid if metabolic flux is determined by the energy demand of different cellular processes in an additive manner. This assumption has not been thoroughly tested. Using the NADH redox model, we discovered a homeostasis of ETC flux in mouse oocytes where perturbing energy demand and supply do not impact ETC flux despite significantly changing NADH metabolic state. On the other hand, perturbing ATP synthesis and proton leak greatly impacted the ETC flux. From these results, we concluded that it is the intrinsic rates of mitochondrial respiration, rather than energy supply or demand, that controls the ETC flux in mouse oocytes. While NADH metabolic state significantly changed in response to perturbing energy demand and supply, indicating cell metabolism was indeed impacted, it is unclear if these perturbations also influenced ATP, ADP, or AMP levels. Future work, including direct measurements of ATP, ADP, and AMP levels, will be required to uncover the mechanism of flux homeostasis. More broadly, our work demonstrates that it is a prerequisite to understand the regulation of ETC fluxes in order to correctly interpret the changes of ETC flux upon inhibiting subcellular processes.
The mechanism of the homeostasis of ETC flux is unclear. One possibility is the presence of flux buffering pathways, where the change of ATP fluxes induced by process inhibition is offset by the opposing change of fluxes through the buffering pathways. Enzymes such as adenylate kinase are known to buffer concentrations of adenine nucleotide (De la Fuente et al., 2014), but it is unclear if they also buffer fluxes. Another possibility is a global coupling of cellular processes, where the change of ATP consumption by one process is offset by the change of others. Changes in proton leak could also compensate for changes in ATP production. Additional work will be required to distinguish between these (and other) possibilities.
FLIM data is obtained with optical resolution, enabling subcellular measurements of NADH metabolic state. Interpreting these measurements using the NADH redox model enables inference of metabolic fluxes with subcellular resolution. Using this method, we discovered a subcellular spatial gradient of ETC flux in mouse oocytes, where the ETC flux is higher in mitochondria closer to the cell periphery. We found that this flux gradient is primarily a result of a spatially heterogeneous mitochondrial proton leak. It will be an exciting aim for future research to uncover the causes and consequences of the subcellular spatial variation in mitochondrial activity.
Materials and methods
Culturing of mouse oocytes
Request a detailed protocolFrozen MII mouse oocytes (Strain B6C3F1) were purchased from EmbryoTech. Oocytes were thawed and cultured in droplets of AKSOM media purchased from MilliporeSigma in plastic petri dish. Mineral oil from VitroLife was applied to cover the droplets to prevent evaporation of the media. Oocytes were then equilibrated in an incubator at 37°C, with 5% CO_{2} and air saturated oxygen before imaging. For imaging, oocytes were transferred to a 2μl media droplet in a 35mm glass bottom FluoroDish from WPI covered with 400–500 μl of oil. The glass bottom dish was placed in an Ibidi chamber with temperature and gas control during imaging. Temperature was maintained at 37°C via heated chamber and objective heater. CO_{2} was maintained at 5% using gas tanks from Airgas.
Cell lines
Request a detailed protocolThe hTERTRPE1 cell line is an established wildtype cell line received from the Cheeseman lab that has been validated based on behavior and properties. The hTERTRPE1 cell line was maintained and tested for mycoplasma contamination in the Needleman lab on a regular basis (Southern Biotech).
Culturing of hTERTRPE1 cells
Request a detailed protocolCell lines were maintained at 37°C and 5% CO_{2}. Cells were grown in Dulbecco’s modified Eagle’s medium (DMEM) (11966025, Gibco) supplemented with 10% fetal bovine serum (FBS), 0.5 mM sodium pyruvate, 5 mM HEPES, 1% penicillin and streptomycin, and either 10 mM glucose or 10 mM galactose. Cells were passaged in glucose or galactose at least three times before imaging. Cells were plated on 35 mm glass bottom FluoroDishes from WPI for imaging. Right before imaging, the media was replaced with 1 ml of phenol redfree DMEM (A1443001, Gibco) supplemented with 0.5 mM sodium pyruvate, 4 mM Lglutamine, 10 mM HEPES, and either 10 mM glucose or 10 mM galactose.
FLIM measurements
Request a detailed protocolOur FLIM system consists of a twophoton confocal microscope with a 40× 1.25 NA water immersion Nikon objective, Becker and Hickle Time Correlated Single Photon Counting (TCSPC) acquisition system and a pulsed MaiTai DeepSee Ti:Sapphire laser from SpectraPhysics. NADH autofluorescence was obtained at 750 nm excitation wavelength with a 460/50 nm emission filter. Laser power at the objective was maintained at 3 mW. The scanning area was 512 by 512 pixels with a pixel size of 420 nm. Acquisition time was 30 s per frame. Oocytes were imaged with optical sectioning across their equators. A histogram of NADH fluorescence decay times was obtained at each pixel of the image.
Oxygen measurements
Request a detailed protocolOxygen level was measured in the Ibidi chamber with an electrodebased oxygen sensor (GasLab). Since the oil layer covering the media droplet was very thin, the oxygen level in the droplet was assumed to be in instant equilibration with the chamber.
Image and FLIM data analysis
Request a detailed protocolTo separate mitochondrial NADH signal from cytoplasmic signal, we performed machine learningbased segmentation algorithms on NADH intensity images. We used the freeware Ilastik (Berg et al., 2019), which implements a supervised learning algorithm for pixel classification. The classifiers were trained to separate mitochondrial pixels from cytoplasmic pixels with a greater than 80% accuracy, as tested by MitoTracker Red FM (Appendix 1, Figure 1—figure supplement 1). We grouped photons from all mitochondrial pixels to obtain a histogram of NADH decay times for each oocyte and for each image of tissue culture cells. To extract the FLIM parameters of NADH bound fraction $f,$ long lifetime ${\tau}_{\mathrm{l}}$ and short lifetime ${\tau}_{\mathrm{s}}$, we fitted the histogram with $G=\mathrm{IRF*}{(C}_{1}F+{C}_{2})$, where $\mathrm{*}$ indicates a convolution, and $\mathrm{I}\mathrm{R}\mathrm{F}$ is the instrument response function of the FLIM system, measured using a urea crystal. $F\left(\tau \right)=f\cdot \mathrm{exp}\left(\frac{\tau}{{\tau}_{\mathrm{l}}}\right)+\left(1f\right)\cdot \mathrm{exp}\left(\frac{\tau}{{\tau}_{\mathrm{s}}}\right)$ is the twoexponential model for the NADH fluorescence decay. ${C}_{1}$ is the amplitude of the decay and ${C}_{2}$ is the background noise. The fitting was performed with a custom MATLAB code using a LevenbergMarquardt algorithm (Yoo, 2018). To obtain the intensity, $I$, of mitochondrial NADH, we first measured the average number of photons per mitochondrial pixel, and divided it by the pixel area, 0.185 μm^{2}, and pixel scanning time 4.09 μs. The flux of ETC is inferred using Equations 5ac for each oocyte and for tissue culture cells in a single image. Heatmaps of mean NADH fluorescence decay times were obtained by computing NADH fluorescence decay time of each mitochondrial pixel and averaging over neighboring mitochondrial pixels weighted by a Gaussian kernel with a standard deviation of 20 pixels. All FLIM measurements were taken from distinct individual oocytes and distinct images of tissue culture cells. Error bars in all figures of FLIM represent standard error of the mean across different individual oocytes or across different images for tissue culture cells. Number of oocytes is reported with n. Number of images for tissue culture cells is reported with N.
Error analysis
Request a detailed protocolFLIM curves were independently fit for each individual oocyte. The reported error bars in this manuscript are standard errors of the mean (SEMs) across these measurements, which depends on the level of variation (the standard deviation) between the oocytes. Two sources of variation in FLIM measurements across the oocytes are: (1) true biological variations between oocytes and (2) fitting errors in the FLIM analysis. To estimate the error of fitting, we performed bootstrapping with randomly drawn points with substitution from each fluorescence decay curve for 53 oocytes. There are ~66,000 photons per oocyte, from which we generated 10 bootstrapped decay curves per oocyte to estimate the fitting error. The fitting error is computed as the variance and covariance of the fitted parameters across bootstrapped decay curves and averaged over 53 oocytes.
At high oxygen level in the AKSOM condition, the bootstrapping yields a variance of 2.2×10^{4}, 4.6×10^{3} ns^{2}, and 6.0×10^{4} ns^{2} for bound fraction, long lifetime, and short lifetime, respectively. The celltocell variances obtained from a single fit per oocyte are 4.4×10^{4}, 9.5×10^{3} ns^{2}, and 1.6×10^{3} ns^{2} for bound fraction, long lifetime, and short lifetime, respectively. Hence the bootstrapping error accounts for 50%, 49%, and 40% of the celltocell variance in bound fraction, long lifetime, and short lifetime, respectively. The bootstrapping yields a covariance of 1.0×10^{3} ns between bound fraction and long lifetime, which only accounts for ~20% of the covariance between these two variables during oxygen drop experiment. The inferred mean flux for oocytes at high oxygen levels in AKSOM is $\u27e8{J}_{\mathrm{o}\mathrm{x}}\u27e9=56.6\text{}\mu \mathrm{M}\cdot {\mathrm{s}}^{1}$. Propagating the error of fitting in all parameters from the bootstrapping analysis to the inferred flux gives a standard error of the mean in ${J}_{\mathrm{o}\mathrm{x}}$ of 1.1 μM·s^{1}. The standard error of the mean in ${J}_{\mathrm{o}\mathrm{x}}$ obtained from a single fit per oocyte was 2.0 μM·s^{1}. Thus, fitting errors account for ∼50% of the standard error of the mean in ${J}_{\mathrm{o}\mathrm{x}}$ .
Metabolic and demand perturbations
Request a detailed protocolOxygen drop experiments for oocytes were performed by mixing nitrogenbalanced 5% O_{2} gas with 0% O_{2} gas at different ratios to create a continuous oxygen drop profile. CO_{2} was maintained at 5%. Oocytes were imaged for 10 min at 5% O_{2}, 30 min during the continuous drop from 5% O_{2} to approximately 0% O_{2}, and 20 min after quickly returning to 5% O_{2}. Oxygen levels were simultaneously monitored with an electrodebased oxygen sensor in the Ibidi chamber. 5% O_{2} corresponds to ~50 μM of oxygen concentration in the culturing media. All the drug perturbations for oocytes were performed by equilibrating oocytes in the AKSOM media containing the corresponding drug for 15–30 min before the oxygen drop experiments. Pyruvate and potassium perturbations were performed by making KSOM media following Cold Spring Harbor Laboratory protocols with varying concentrations of sodium pyruvate and potassium, respectively. For oligomycin, FCCP, rotenone and pyruvate perturbations, 9 mM of sodium oxamate was also added to the media to suppress cytoplasmic NADH signal for better mitochondrial segmentation. The addition of the oxamate does not change the ETC flux of the mitochondria (Figure 5b).
For hTERTRPE1 cells, drug perturbations were performed by replacing the media with drugcontaining media through pipetting. Cells were imaged for 20–30 min immediately after drug perturbations.
All drugs were purchased from SigmaAldrich. Temperature was maintained at 37°C. CO_{2} was maintained at 5%.
Oxygen consumption rate measurement
Request a detailed protocolThe OCR of the oocytes was measured using the nanorespirometer from Unisense (Lopes et al., 2005). A batch of 10–15 oocytes was placed at the bottom of a glass capillary with a diameter of 0.68 mm and a height of 3 mm. The capillary well is filled with AKSOM media or drugcontaining media for metabolic perturbations. After an equilibration time of ~2 hr, a steadystate linear oxygen gradient is established in the capillary well due to the balance of oocyte respiration and oxygen diffusion. A motorcontrolled electrodebased oxygen sensor (Unisense) is used to measure the oxygen gradient. The OCR is calculated as the product of the oxygen gradient, diffusivity of oxygen in the media, taken to be 3.37×10^{5} cm^{2}/s, and the crosssectional area of the capillary well, which was 0.36 mm^{2}. The entire system was enclosed in a custombuilt chamber with temperature and gas control. Temperature was maintained at 37°C. Oxygen level was continuously varied during oxygen drop experiments by slowly mixing 20% O_{2} with 0% O_{2} from gas tanks, and maintained at the air saturation level for drug and pyruvate perturbations. OCR was measured on a group of 10–15 oocytes at a time. Singleoocyte OCR was obtained by dividing the measured OCR by the number of oocytes in the group. Error bars in all figures of OCR represent standard error of the mean across different groups of oocytes normalized by the number of oocytes in each group. Number of oocytes is reported with n. Number of groups is reported with N.
Statistical analysis
Request a detailed protocolFor the comparison between inferred ETC flux and measured ETC flux of the oocytes, twosample ttest was performed on the vectors of inferred singlecell ETC flux (with n elements, where n is the number of oocytes) and the batch OCR measurements (with N elements, where N is the number of batch groups). For the comparison between inferred ETC flux and measured ETC flux of the tissue culture cells, twosample ttest was performed on the vectors of inferred relative change of ETC flux (with n elements, where n is the number of images) and the relative change of OCR estimated from Figure 1A of MacVicar and Lane, 2014 (with N elements, where N is the estimated number of OCR data points).
Mitochondrial membrane potential measurement
Request a detailed protocolThe spatial distribution of mitochondrial membrane potential within oocytes was measured with a potentialsensitive dye TMRM (SigmaAldrich). Oocytes were cultured in AKSOM with 100 nM TMRM for 30 min before imaging. TMRM signal was obtained at 830 nm excitation wavelength with 560/40 nm emission filter. Mitochondrial TMRM intensity in different regions of the oocyte was computed by dividing the total number of photons from that region by the number of pixels in the same region. Heatmaps of mitochondrial TMRM intensity were obtained by computing photon counts for each mitochondrial pixel and averaging over neighboring mitochondrial pixels weighted by a Gaussian kernel with a standard deviation of 20 pixels. To normalize TMRM intensity by mitochondrial mass, we cultured oocytes in AKSOM with 100 nM MitoTracker Red FM and 25 nM TMRM for 30 min before imaging. We also cultured oocytes in AKSOM with 1 μg/ml JC1 dye for 3 hr before imaging.
Mitochondrial membrane potential of hTERTRPE1 cells was measured with TMRM. The cells were cultured in DMEM with 100 nM TMRM for 15–30 min before imaging. To measure membrane potential under drug perturbations, the original media was pipetted out and replaced with media containing both 100 nM TMRM and the drug. The cells were imaged for 20–30 min immediately after drug perturbations. TMRM intensity ratio was obtained by normalizing the mitochondrial TMRM intensity by the cytoplasmic TMRM intensity.
Appendix 1
Segmentation of mitochondria and calculation of NADH concentrations
Segmentation of mitochondria
We used Ilastik, a machine learningbased software for image analysis, to classify pixels in the NADH intensity images containing mitochondria (Berg et al., 2019). For each experiment, we generated a timelapse movie of NADH (Video 1). We used a few images in the movie as the training data set to train the software to classify mitochondrial pixels by manually selecting clustered high brightness pixels. Other pixels are classified as either cytoplasm or background. We then applied the trained pixel classifier to generate a mitochondrial probability map for each image in the entire movie, with each pixel assigned a probability between 0 and 1 to be mitochondrial pixel. Pixels with a probability higher than 0.7 were considered to be mitochondrial pixels.
To test the accuracy of this segmentation algorithm, we immersed the oocytes in AKSOM media containing MitoTracker Red FM, a dye that specifically labels mitochondria. Pixels with intensity above 60 percentile in the MitoTracker image were considered to be mitochondrial pixels. We imaged NADH and MitoTracker for the same oocyte and compared the resulting distribution of mitochondria (Figure 1—figure supplement 1). We defined the accuracy of the NADHbased segmentation as the fraction of photons originating from true mitochondrial pixels. The accuracy of the segmentation is 80.6±1% (SEM) for the control condition as averaged over seven oocytes. We repeated the analysis for oxamate, oligomycin, FCCP and rotenone perturbations, and obtained an accuracy of segmentation of 78.6±1.4%, 84.1±1.6%, 83.7±0.5%, and 81.7±2%, respectively, similar to the control condition.
Converting NADH intensity to NADH concentrations
Since the molecular brightness of NADH depends on the fluorescence lifetime of NADH, which changes drastically upon binding enzymes, the NADH concentration is not linearly proportional to NADH intensity. FLIM provides an accurate way of measuring NADH concentrations by simultaneously measuring fluorescence intensity and lifetime. We now derive the NADH intensityconcentration relation from the FLIM measurements. Assuming molecular brightness is proportional to the fluorescence lifetime, and therefore that free and bound NADH have different contributions to the measured intensity, we have
where $I$ is the intensity of NADH and ${c}_{\mathrm{s}}$ is a calibration factor that depends on the laser power. From Equation S1, we obtained the concentrations of free and bound NADH:
where $f$ is the fraction of bound NADH.
To get the calibration factor ${c}_{\mathrm{s}}$, we titrated NADH in AKSOM solutions and fitted the calibration curve using:
where ${\tau}_{\mathrm{s}\mathrm{o}\mathrm{l}}$ is the lifetime of NADH in solution. ${\tau}_{\mathrm{s}\mathrm{o}\mathrm{l}}$ was directly measured by FLIM (Figure 1—figure supplement 2), allowing us to obtain ${c}_{\mathrm{s}}$ from the fit (Figure 1—figure supplement 2).
FLIM can be used to accurately measure concentrations of bound and free NADH in vitro
To test if absolute concentrations of free and bound NADH can be accurately measured from FLIM of NADH, we prepared solutions with known total concentration of NADH, and titrated the concentration of LDH, an enzyme to which NADH can bind. We prepared the solutions with 50 mM TRIS buffer, 150 mM NaCl at pH 7.6 and 37°C. We added a total concentration of 50 μM NADH to the solution and titrated LDH concentrations at 0, 1.4, and 3.5 μM. We first performed single exponential fitting of the NADH decay curve at 0 μM LDH, where all NADH are free, to obtain the NADH short lifetime (Figure 1—figure supplement 3d). From the NADH intensity (Figure 1—figure supplement 3a), we obtained the calibration factor ${c}_{\mathrm{s}}$ using Equation S4 with $\left[\mathrm{N}\mathrm{A}\mathrm{D}{\mathrm{H}}_{\mathrm{s}\mathrm{o}\mathrm{l}}\right]=50\mu M$. We then fixed the short lifetime and performed twoexponential fitting of the NADH decay curve at LDH concentrations of 1.4 μM and 3.5 μM to obtain the bound ratio (Figure 1—figure supplement 3b) and long lifetime (Figure 1—figure supplement 3c). As expected, NADH bound ratio increases with LDH concentrations, as there is more enzyme for NADH to bind. Finally, we calculated free NADH concentration $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and bound NADH concentration $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$ using Equations S2 and S3 from the FLIM parameters. Remarkably, the free and bound concentrations of NADH both change with LDH concentrations but the total concentration remains at 50 μM (Figure 1—figure supplement 3e). This result shows that Equations S2S3 can be used to accurately measure the concentrations of free and bound NADH from FLIM measurements of NADH.
Appendix 2
Reversible MichaelisMenten kinetics, full and reduced notations
The kinetic equations of the reversible MichaelisMenten kinetics (Figure 3a, left) for the ith oxidase are
where $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ is the concentration of free NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}}_{\mathrm{f}}^{+}\right]$ is the concentration of free NAD^{+}, $\left[{\mathrm{O}\mathrm{x}}_{i}\right]$ is the concentration of free oxidase, and $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet {\mathrm{o}\mathrm{x}}_{i}\right]$ is the concentration of the NADHoxidase complex. ${k}_{1}$, ${k}_{1}$, ${k}_{2}$, and ${k}_{2}$ are the forward and reverse reaction rates.
In the reduced notation as introduced in Figure 3a (right), the same enzyme kinetics is described by
where ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{b}}$ and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{u}}$ are the effective binding and unbinding rates of NADH to the oxidase, and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{b}}$ and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{u}}$ are the effective binding and unbinding rates of NAD^{+} to the oxidase. In this reduced notation, the concentration of the free oxidase $\left[\mathrm{O}{\mathrm{x}}_{i}\right]$ is absorbed into the effective binding rates: that is, ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{u}}={k}_{1}$, ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{b}}={k}_{1}\left[{\mathrm{O}\mathrm{x}}_{i}\right]$, ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{u}}={k}_{2}$, and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{b}}={k}_{2}\left[{\mathrm{O}\mathrm{x}}_{i}\right]$. Hence $\left[\mathrm{O}{\mathrm{x}}_{i}\right]$ becomes an implicit variable whose behavior is not evident from the reduced notation diagram (Figure 3a, right). Modeling the full dynamics of a reversible MichalisMenten enzyme requires specifying the equation for $\left[\mathrm{O}{\mathrm{x}}_{i}\right]$:
The reduced notation (Figure 3a, right; Equations S9–S12) and the full notation (Figure 3a, left; Equations S5S8) are mathematically identical and describe the exact same kinetics.
Generalized enzyme kinetics, reduced notation
The reduced notation for the ith oxidase displaying generalized enzyme kinetics (Figure 3b) refers to the following class of mathematical models:
where $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ is the concentration of free NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}}_{\mathrm{f}}^{+}\right]$ is the concentration of free NAD^{+}, $\left[{\mathrm{O}\mathrm{x}}_{i}\right]$ is the concentration of free oxidase, $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet {\mathrm{O}\mathrm{x}}_{i}\right]$ is the concentration of the NADHoxidase complex, and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}}^{+}\bullet {\mathrm{O}\mathrm{x}}_{i}\right]$ is the concentration of the NAD^{+}oxidase complex. ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{b}}$ and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{u}}$ are the effective binding and unbinding rates of NADH to the oxidase, ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{b}}$ and ${k}_{{\mathrm{o}\mathrm{x}}_{i}}^{\mathrm{`}\mathrm{u}}$ are the effective binding and unbinding rates of NAD^{+} to the oxidase, and ${r}_{{\mathrm{o}\mathrm{x}}_{i}}^{+}$ and ${r}_{{\mathrm{o}\mathrm{x}}_{i}}^{}$ are the forward and reverse oxidation rates. These rates can be arbitrary functions of implicit variables, such as the concentration of free oxidase, $\left[{\mathrm{O}\mathrm{x}}_{i}\right]$, the mitochondrial membrane potential, $\mathrm{\Delta}{G}_{\mathrm{H}}$, $\mathrm{p}\mathrm{H}$, and other factors:
Thus, while Equations S13S16 superficially appear to be linear and first order, they can actually refer to nonlinear reactions of any order because the rate can depend on $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and other variables (Equations S17a–f).
The implicit variables that these rates depend on can each be governed by their own dynamics that are arbitrary functions of other variables:
Describing the dynamics of the enzyme requires specifying the implicit variables that the rates depend on, the functional form of these dependencies, and the additional equations for the dynamics of the implicit variables (Equations S18). However, we will show that the predicted relationship between FLIM measurements of NADH and fluxes does not depend on these modeling choices. Thus, the reduced notation is convenient for deriving these relations for a broad class of models.
Appendix 3
Coarsegraining the detailed NADH redox model
We consider an NADH redox loop consisting of M reductases and N oxidases (Figure 4—figure supplement 1), each of which is described by the generalized enzyme kinetics (Figure 3b; Equations S13S16). We coarsegrain this detailed NADH redox model by coarsegraining all oxidases into a single effective oxidase and all reductases into a single effective reductase (Figure 4b). We relate the kinetic rates of the coarsegrained model to those of the detailed model by keeping the global binding and unbinding fluxes and the global reaction fluxes through the oxidases and reductases the same as the detailed model.
We first coarsegrain the oxidases and reductases:
We require the global binding and unbinding fluxes of NADH to the effective oxidase and reductase to be equal to the sum of their binding and unbinding fluxes to all of the individual oxidases and reductases:
which leads to
We require the global forward and reverse reaction flux through the effective oxidase and reductase to be equal to the sum of the reaction fluxes through all of the individual oxidases and reductases:
which leads to
By applying the same procedure to ${\mathrm{N}\mathrm{A}\mathrm{D}}^{+}$, we can obtain the effective reduction rates ${r}_{\mathrm{r}\mathrm{e}}^{+}$, ${r}_{\mathrm{o}\mathrm{x}}^{}$ and the effective binding and unbinding rates of ${\mathrm{N}\mathrm{A}\mathrm{D}}^{+}$ : ${k}_{\mathrm{o}\mathrm{x}}^{\mathrm{`}\mathrm{b}}$, ${k}_{\mathrm{r}\mathrm{e}}^{\mathrm{`}\mathrm{b}}$, ${k}_{\mathrm{o}\mathrm{x}}^{\mathrm{`}\mathrm{u}}$, ${k}_{\mathrm{r}\mathrm{e}}^{\mathrm{`}\mathrm{u}}$. We omit the derivation here because these rates are not needed to infer ETC flux. We hence explicitly related the kinetic rates of the coarsegrained model (Figure 4b) to those of the detailed model (Figure 4—figure supplement 1). We note that under the generalized enzyme kinetics (Figure 3b; Equations S13S16), all kinetic rates are considered to be general functions of enzyme concentrations, metabolite concentrations, and other factors, and thus, all of the rates in the coarsegrained model can also depend on all of those factors. These implicit variables can obey their own dynamical equations (Equation S18). The coarsegraining presented here is mathematically exact and independent of both the functional forms of these rates and the functional form of the dynamic equations of the implicit variables.
Appendix 4
Predicting the ETC flux using the coarsegrained NADH redox model
The coarsegrained NADH redox model
We start with the equations characterizing the dynamics of the coarsegrained NADH redox model as described in Figure 4b:
where $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}}_{\mathrm{f}}^{+}\right]$ are the concentrations of free NADH and free NAD^{+}; $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{R}\mathrm{e}\right]$ and $[{\mathrm{N}\mathrm{A}\mathrm{D}}^{+}\bullet \mathrm{R}\mathrm{e}]$ are concentrations of reductasebound NADH and NAD^{+}; $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}\right]$ and $[{\mathrm{N}\mathrm{A}\mathrm{D}}^{+}\bullet \mathrm{O}\mathrm{x}]$ are concentrations of oxidasebound NADH and NAD^{+}; $k$ denotes binding (b) and unbinding (u) rates, with subscript re and ox denoting reductase and oxidase, respectively; ${r}_{\mathrm{r}\mathrm{e}}^{+}$ and ${r}_{\mathrm{r}\mathrm{e}}^{}$ are the forward and reverse reaction rates of the reductase; ${r}_{\mathrm{o}\mathrm{x}}^{+}$ and ${r}_{\mathrm{o}\mathrm{x}}^{}$ are the forward and reverse reaction rates of the oxidase. The reaction rates, and binding and unbinding rates, can be arbitrary functions of metabolite concentrations, enzyme concentrations, and other variables (such as membrane potential, oxygen concentration, etc, each of which can obey their own dynamical equations).
Predicting the ETC flux
The flux through the ETC is
At steadystate (or in the quasistatic limit), all the time derivatives are zero. Setting $\mathrm{d}[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\cdot \mathrm{O}\mathrm{x}]/\mathrm{d}\mathrm{t}$ (Equation S32) to zero, we obtained the steadystate flux through the ETC:
Setting $\mathrm{d}[\mathrm{N}\mathrm{A}\mathrm{D}{\mathrm{H}}_{\mathrm{f}}]/\mathrm{d}\mathrm{t}$ (Equation S31) to zero gives:
and using:
from which we solved for $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}\right]:$
Substituting $\left[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}\right]$ in Equation S37 with Equation S40, we obtained our central result:
From Equation S41, we see that the flux through the ETC is a product of the turnover rate of free NADH, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, and the concentration of free NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, where
and
where we defined the NADH bound ratio and its equilibrium counterpart as:
Appendix 5
Flux inference procedures using the coarsegrained NADH redox model
Equations S41S42, or equivalently Equations 5ac from the main text, can be used to infer the flux through the ETC, ${J}_{\mathrm{o}\mathrm{x}}$, from FLIM measurements of NADH across a wide range of metabolic perturbations (Figure 5—figure supplement 1). To do so, we infer the turnover rate of free NADH, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, and the concentration of free NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$. The product of $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ gives ${J}_{\mathrm{o}\mathrm{x}}$. $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ can be obtained using Equation S2 (Figure 5—figure supplement 1). In this section, we describe two procedures to obtain $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$: one from the measurement of NADH bound ratio $\beta $, and the other from the measurement of NADH long fluorescence lifetime ${\tau}_{\mathrm{l}}$.
Inferring $\stackrel{\mathbf{~}}{\mathit{r}}}_{\mathbf{o}\mathbf{x}$ from NADH bound ratio $\mathit{\beta}$
Equation S42, ${\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}}=\text{}\alpha \left(\beta {\beta}_{\mathrm{e}\mathrm{q}}\right)$, provides a method to obtain $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$. We measure the NADH bound ratio, $\beta $, using $\beta =f\phantom{\rule{thinmathspace}{0ex}}/(1f)$, where $f$ is the NADH bound fraction obtained by fitting the fluorescence decay curve of NADH (see Materials and methods). We obtain the equilibrium bound ratio, ${\beta}_{\mathrm{e}\mathrm{q}}$, by dropping the oxygen level to the lowest achievable value with our setup: $\left[{\mathrm{O}}_{2}\right]=0.26\pm 0.04\mathrm{\mu}\mathrm{M}$, and assuming ${\beta}_{\mathrm{e}\mathrm{q}}$ does not change with oxygen levels. Note that ${\beta}_{\mathrm{e}\mathrm{q}}$ does change with drug perturbations, and therefore needs to be separately determined for each condition (Figure 5—figure supplement 1h). We obtain $\alpha $ using direct measurement of ${J}_{\mathrm{o}\mathrm{x}}$ from OCR measurements:
where $V}_{\mathrm{m}}=9.5\times {10}^{4}\text{}{\mu \mathrm{m}}^{3$ is the average volume of mitochondria per oocyte approximated from the area fraction of mitochondria based on the segmentation, where the mitochondrial area fraction is estimated at 46% and oocyte volume at 2×10^{5} µm^{3}. Using $\mathrm{O}\mathrm{C}\mathrm{R}=2.68\pm 0.06\phantom{\rule{thinmathspace}{0ex}}\mathrm{f}\mathrm{m}\mathrm{o}\mathrm{l}/\mathrm{s}$ per oocyte in the control condition (AKSOM media at 50 μM oxygen level), we get $\alpha =5.4\pm 0.2\text{}{\mathrm{s}}^{1}$. $\alpha $ is approximated as a constant that does not vary with perturbations, hence $\alpha $ calibrated at one condition can be used for all other conditions (as confirmed by the agreement between FLIM based inference and OCR measurements in Figures 5 and 8).
Once $\alpha $ is calibrated at the control condition using Equation S48, and ${\beta}_{\mathrm{e}\mathrm{q}}$ is determined from an oxygen drop experiment, then subsequent FLIM measurements of $\beta $ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ can be used with Equation S42 and S41 to determine the absolute value of $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ and ${J}_{\mathrm{o}\mathrm{x}}$ for all conditions (Figure 5).
Inferring $\stackrel{\mathbf{~}}{\mathit{r}}}_{\mathbf{o}\mathbf{x}$ from NADH long fluorescence lifetime $\mathit{\tau}}_{\mathbf{l}$
In this section, we derive an alternative procedure for determining the turnover rate of free NADH $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$, and hence ${J}_{\mathrm{o}\mathrm{x}}$, using changes in the NADH long fluorescence lifetime. The NADH long fluorescence lifetime, ${\tau}_{\mathrm{l}}$, is associated with enzymebound NADH (Sharick et al., 2018). In the coarsegrained NADH redox model described above, and in Figure 4b, the enzymebound NADH consists of reductasebound NADH ($[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{R}\mathrm{e}]$) and oxidasebound NADH ($[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}]$). We therefore assume that the experimentally measured NADH long lifetime, ${\tau}_{\mathrm{l}}$, is a linear combination of the lifetimes of $[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}]$ and $[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{R}\mathrm{e}]$:
where ${\tau}_{\mathrm{o}\mathrm{x}}$ and ${\tau}_{\mathrm{r}\mathrm{e}}$ are the fluorescence lifetimes corresponding to the oxidasebound NADH and reductasebound NADH, respectively. Solving for $[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{O}\mathrm{x}]$ and $[\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}\bullet \mathrm{R}\mathrm{e}]$ as a function of $\beta $ using Equations S39S40 and S44, and substituting into Equation S49, we predict that the NADH long fluorescence lifetime ${\tau}_{\mathrm{l}}$ is linearly related to the inverse of the NADH bound ratio $1/\beta $:
with
The predicted linear relationship between ${\tau}_{\mathrm{l}}$ and $1/\beta $ is empirically observed during oxygen drop experiments, as shown in Figure 6a and Figure 6—figure supplement 1. This is a selfconsistency check that argues for the validity of the assumption in Equation S49.
At equilibrium, when there is no flux through the ETC (i.e., ${J}_{\mathrm{o}\mathrm{x}}=0$), Equation S50 gives:
where ${\tau}_{\mathrm{e}\mathrm{q}}$ is the NADH long lifetime at equilibrium. Solving for $\beta $ and ${\beta}_{\mathrm{e}\mathrm{q}}$ as a function of ${\tau}_{\mathrm{l}}$ and ${\tau}_{\mathrm{e}\mathrm{q}}$ from Equations S50 and S53 and substituting into Equation S42, we obtain $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ in terms of ${\tau}_{\mathrm{l}}$ :
where A and B are the slope and offset of the linear relation between ${\tau}_{\mathrm{l}}$ and $1/\beta $ in Equation S50.
We experimentally measured A and B for each oocyte from the slope and offset of a linear fit between ${\tau}_{\mathrm{l}}$ and $1/\beta $ during oxygen drop experiments across all drug perturbations (Figure 6a; Figure 6—figure supplement 1). We obtained the equilibrium long lifetime, ${\tau}_{\mathrm{e}\mathrm{q}}$, by FLIM measurements at the lowest achievable oxygen level in our set up: $\left[{\mathrm{O}}_{2}\right]=0.26\pm 0.04\mathrm{\mu}\mathrm{M}$. Once A, B, and ${\tau}_{\mathrm{e}\mathrm{q}}$ are measured, Equation S54 can be used to determine $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ from FLIM measurements of ${\tau}_{\mathrm{l}}$. If $\alpha $ is not known, this procedure can only be used to obtain $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ up to a constant of proportionality. If $\alpha $ is independently measured from Equation S48 at one condition, then Equation S54 can be used to determine the absolute value of $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ for all conditions (Figure 6b).
As described in the main text, $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ inferred from ${\tau}_{\mathrm{l}}$ using Equation S54 produces the same results as $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ inferred from $\beta $ using Equation S42 (Figure 6b). The agreement between these two methods is a strong selfconsistency check of the NADH redox model.
Accounting for NADPH and other background fluorescence
Equation S41 provides a method to infer ETC flux since all factors in it, except for the constant of proportionality, $\alpha $, depend only on $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$, which can be measured from FLIM of NADH. One potential complication with this procedure is that NADPH, another autofluorescent electron carrier, shares a similar fluorescence spectrum with NADH, resulting in a mixed NAD(P)H signal from the autofluorescence measurement. While NADH concentration is 40 times greater than the concentration of NADPH for the whole mouse oocytes (Bustamante et al., 2017), and presumably even higher for mitochondria, NADPH concentration can be comparable to that of NADH for other cell types such as tissue culture cells (Park et al., 2016). In this section, we generalize Equation S41 to predict the ETC flux by explicitly considering the potential contributions of other fluorescence species, such as NADPH, to the measured autofluorescence signal.
We start from the fact that concentrations of bound and free fluorescent species measured from FLIM using Equations S2S3, $\left[{\mathrm{N}}_{\mathrm{f}}\right]$ and $\left[{\mathrm{N}}_{\mathrm{b}}\right]$, could be different from the actual concentrations of free and bound NADH, $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$. If the signal from NADPH and other additional fluorescence species is additive, then:
where ${C}_{\mathrm{f}}$ and ${C}_{\mathrm{b}}$ are the nonNADH contributions to the measured concentrations of free and bound fluorescent species. We substitute Equations S55S56 to the predicted ETC flux in Equation S41 and obtain
Rearranging, we obtain
where
Comparing Equation S58 with Equation S41, we notice that the background fluorescence does not change the form of the equation of the predicted ETC flux because the concentrations of the background fluorescent species are incorporated into the equilibrium bound ratio ${\beta}_{\mathrm{N},\mathrm{e}\mathrm{q}}$ . If ${\beta}_{\mathrm{N},\mathrm{e}\mathrm{q}}$ can be reliably measured, the background fluorescence will not affect the flux inference procedures. In other words, $\left[{\mathrm{N}}_{\mathrm{f}}\right]$ and $\left[{\mathrm{N}}_{\mathrm{b}}\right]$ can be used for flux inference in place of $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{f}}\right]$ and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$ in Equations S41S42. Therefore, an additive offset to the measured concentrations of free and bound species will not affect the flux inference procedure, whether that additive offset comes from NADPH or from other sources of fluorescent background.
Alternatively, if the signal from background fluorescence changes proportionally with NADH, then:
we have
where
Comparing Equation S63 with Equation S41, we again obtained the same form for ${J}_{\mathrm{o}\mathrm{x}}$, but with a rescaled ${\alpha}_{\mathrm{N}}$. Therefore, a background fluorescence signal that changes proportionally with NADH will not affect the flux inference procedure, whether that background comes from NADPH or from other sources.
Appendix 6
Spatial gradient of mitochondrial metabolism in mouse oocytes
${\mathit{\beta}}_{\mathbf{e}\mathbf{q}}$ is uniform within the oocyte
To obtain subcellular ETC flux as a function of distance to the oocyte’s center using Equations 5ac in the main text, we need to know the spatial variation of ${\beta}_{\mathrm{e}\mathrm{q}}$. While the NADH bound ratio at the lowest oxygen level gives a good approximation for the average ${\beta}_{\mathrm{e}\mathrm{q}}$ of the cell (Figure 5—figure supplement 1h), subpopulations of mitochondria closer to the cell periphery are exposed to slightly higher oxygen level than those away from the cell periphery, obscuring the determination of the spatial variation of ${\beta}_{\mathrm{e}\mathrm{q}}$ from oxygen drop experiment. Hence to obtain the spatial variation of ${\beta}_{\mathrm{e}\mathrm{q}}$ throughout the oocyte, we inhibited the ETC completely using 15 μM of rotenone, an inhibitor of complex I in the ETC, for an extended period of time until the NADH bound ratio reaches the lowest level. We then fitted the NADH decay curves from mitochondrial pixels within equaldistanced concentric rings (Figure 9—figure supplement 1a) to obtain ${\beta}_{\mathrm{e}\mathrm{q}}$ as a function of distance from the oocyte’s center (Figure 9—figure supplement 1b). A linear fit yielded a slope of 0.001±0.0012 (SEM), which is statistically indistinguishable from 0 (p=0.42). Therefore, the resulting ${\beta}_{\mathrm{e}\mathrm{q}}$ is uniform throughout the oocyte and is equal to the average ${\beta}_{\mathrm{e}\mathrm{q}}$ obtained by fitting the decay curve from all mitochondrial pixels in the oocyte at the lowest oxygen level (Figure 5—figure supplement 1h). Hence, we used a constant ${\beta}_{\mathrm{e}\mathrm{q}}$ throughout the oocyte to compute the subcellular ETC flux (Figure 9d).
Subcellular spatial gradient of mitochondrial membrane potential
As shown in the main text, we observed a strong spatial gradient of the intensity of TMRM in mitochondria in oocytes. TMRM is a potentialsensitive dye that preferentially accumulates in mitochondria with higher membrane potential (Figure 9g and h). To test whether this spatial gradient is due to the subcellular variation of mitochondrial membrane potential or the variation in mitochondrial mass, we labelled mitochondria with a potentialinsensitive dye MitoTracker Red FM to quantify mitochondrial mass, together with TMRM. We did not observe a strong gradient of MitoTracker intensity (Figure 9—figure supplement 2b, f) as compared to TMRM intensity (Figure 9—figure supplement 2a, e) within the same oocyte, indicating the mitochondrial mass is uniformly distributed. We further normalized the TMRM intensity by the MitoTracker intensity, and observed a strong spatial gradient of the ratio (Figure 9—figure supplement 2c, g). These results suggest that the spatial gradient of TMRM is due to the variation of mitochondrial membrane potential, rather than the variation of mitochondrial mass. Finally, to test the robustness of the result, we used an alternative potentialsensitive dye JC1, and observed a similar spatial gradient of mitochondrial membrane potential (Figure 9—figure supplement 2d, h). Taken together, these results show that the subcellular spatial gradient of mitochondrial membrane potential is a robust observation that does not depend on the variation of mitochondrial mass or the type of dye used.
Appendix 7
Flux prediction for a NADH redox model with each enzyme described by the reversible MichaelisMenten kinetics
In this section, we derive the ETC flux for a NADH redox model where each of the N oxidase and M reductase obeys reversible MichaelisMenten kinetics (Equations S9S11). We achieve this by reducing the flux prediction of the generalized enzyme kinetics to that of the reversible MichaelisMenten kinetics.
We first consider an NADH redox model with a single oxidase and a single reductase, each of which obeys reversible MichaelisMenten kinetics (i.e., N=M=1). The flux through the oxidase is:
We show that Equations S41S43 that characterize the flux of the generalized NADH redox model can be reduced to Equation S67 that characterizes the flux of the reversible MichaelisMenten model in the limit:
In this limit, we have from Equation S43:
where ‘MM’ stands for ‘MichaelisMenten.’ Similarly, from Equation S45S47, we have
From Equation S40, we have in this limit:
Substituting the expressions for ${\alpha}_{\mathrm{M}\mathrm{M}}$, ${\beta}_{\mathrm{e}\mathrm{q}}^{\mathrm{M}\mathrm{M}}$, and $\left[{\mathrm{N}\mathrm{A}\mathrm{D}\mathrm{H}}_{\mathrm{b}}\right]$ from Equations S69S71 to the predicted ETC flux in Equation S41, we obtain
Thus, we have shown that the flux of the generalized NADH redox model reduces to the flux of the reversible MichaelisMenten model in the limit where the unbinding rate of NADH from the reductase is much faster than any other rates in the model.
We note that the predicted fluxconcentration relation for the reversible MichaelisMenten model (Equation S72) remains exactly the same as the generalized model (Equation S41), but with different expressions for $\alpha $ and ${\beta}_{\mathrm{e}\mathrm{q}}$ as expressed in Equations S69S70.
Next, we generalize the results in Equation S72 to a detailed model with N oxidase and M reductase, each of which is described by the reversible MichaelisMenten kinetics. Unpacking the coarsegrained binding and unbinding rates from Equations S24S25, we obtain
where ${k}_{1,i}$ and ${k}_{1,i}$ denote the binding and unbinding rates of NADH to the ith oxidase.
Connecting the NADH redox model to detailed biophysical models of mitochondrial metabolism
In this section, we show that the coarsegrained NADH redox model described above, and in Figure 4b of the main text, can be directly related to detailed biophysical models of mitochondrial metabolism, including previously published models (Beard, 2005; Korzeniewski and Zoladz, 2001; Hill, 1977; Jin and Bethke, 2002; Chang et al., 2011).
In mitochondria, NADH oxidation is catalyzed by complex I of the ETC, which has the overall reaction:
where two electrons are transferred from NADH to ubiquinone Q, and four protons are pumped out of the mitochondrial matrix. To connect our model with detailed model of complex I, we rewrite the flux through the ETC:
using
where
where
The last equality in Equation S84 is obtained by assuming that the total concentration of NADH plus NAD^{+} is constant:
Equation S81 allows us to connect our coarsegrained model to previously published detailed models of complex I. By equating the flux through complex I, ${J}_{\mathrm{C}1}$, in previous models to the flux through the ETC in our NADH redox model, ${J}_{ox}$, we can determine $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}$ (and $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}}^{+$ and $\stackrel{~}{r}}_{\mathrm{o}\mathrm{x}}^{$) in terms of variables defined in those more detailed models. In Appendix 7—table 1, we summarize the relationship between the NADH redox model and several previously published models of complex I.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting file; Source Data files have been provided for Figures 2, Figure 5, Figure 5—figure supplement 1, Figure 6, Figure 6—figure supplement 1, Figure 7, Figure 8, Figure 8—figure supplement 1.
References

The spatiotemporal dynamics of mitochondrial membrane potential during oocyte maturationMolecular Human Reproduction 25:695–705.https://doi.org/10.1093/molehr/gaz055

Fluorescence lifetime imagingtechniques and applicationsJournal of Microscopy 247:119–136.https://doi.org/10.1111/j.13652818.2012.03618.x

A genetically encoded fluorescent reporter of ATP:ADP ratioNature Methods 6:161–166.https://doi.org/10.1038/nmeth.1288

ilastik: interactive machine learning for (bio)image analysisNature Methods 16:1226–1232.https://doi.org/10.1038/s4159201905829

The development of fertilized human ova to the blastocyst stage in KSOM(AA) medium: is a twostep protocol necessary?Reproductive Biomedicine Online 5:133–140.https://doi.org/10.1016/s14726483(10)61615x

Separating NADH and NADPH fluorescence in live cells and tissues using FLIMNature Communications 5:3936.https://doi.org/10.1038/ncomms4936

Assessing mitochondrial dysfunction in cellsThe Biochemical Journal 435:297–312.https://doi.org/10.1042/BJ20110162

The role of mitochondria in agingThe Journal of Clinical Investigation 123:951–957.https://doi.org/10.1172/JCI64125

RESPIRATORY ENZYMES IN OXIDATIVE PHOSPHORYLATIONJournal of Biological Chemistry 217:383–393.https://doi.org/10.1016/S00219258(19)571897

The role of mitochondrial function in the oocyte and embryoCurrent Topics in Developmental Biology 77:21–49.https://doi.org/10.1016/S00702153(06)770028

Advances in measuring cellular bioenergetics using extracellular fluxDrug Discovery Today 13:268–274.https://doi.org/10.1016/j.drudis.2007.12.008

BookNatural Biomarkers for Cellular Metabolism Biology, Techniques, and ApplicationsCRC Press Taylor & Francis Group.https://doi.org/10.1201/b17427

Mitochondrial respiration in the low oxygen environment of the cell. Effect of ADP on oxygen kineticsBiochimica et Biophysica Acta 1365:249–254.https://doi.org/10.1016/s00052728(98)000760

Oxygen consumption and energy metabolism of the early mouse embryoMolecular Reproduction and Development 44:476–485.https://doi.org/10.1002/(SICI)10982795(199608)44:4<476::AIDMRD7>3.0.CO;2I

Kinetics of electron transfer through the respiratory chainBiophysical Journal 83:1797–1808.https://doi.org/10.1016/S00063495(02)739453

A model of oxidative phosphorylation in mammalian skeletal muscleBiophysical Chemistry 92:17–34.https://doi.org/10.1016/s03014622(01)001843

Lactate regulates pyruvate uptake and metabolism in the preimplantation mouse embryoBiology of Reproduction 62:16–22.https://doi.org/10.1095/biolreprod62.1.16

Extraction and Quantitation of Nicotinamide Adenine Dinucleotide Redox CofactorsAntioxidants & Redox Signaling 28:167–179.https://doi.org/10.1089/ars.2017.7014

Kinetics of the Reversible MichaelisMenten Mechanism and the Applicability of the Steadystate Approximation 1Journal of the American Chemical Society 80:5146–5151.https://doi.org/10.1021/ja01552a034

Quantifying intracellular rates of glycolytic and oxidative ATP production and consumption using extracellular flux measurementsThe Journal of Biological Chemistry 292:7189–7207.https://doi.org/10.1074/jbc.M116.774471

Metabolite concentrations, fluxes and free energies imply efficient enzyme usageNature Chemical Biology 12:482–489.https://doi.org/10.1038/nchembio.2077

Permeabilized cell and skinned fiber techniques in studies of mitochondrial function in vivoMolecular and Cellular Biochemistry 184:81–100.https://doi.org/10.1023/A:1006834912257

Mitochondrial dynamics and division in budding yeastTrends in Cell Biology 12:178–184.https://doi.org/10.1016/s09628924(01)022462

In vivo kinetics and the reversible MichaelisMenten modelJournal of Chemical Education 69:981.https://doi.org/10.1021/ed069p981

Metabolic fluxes and metabolic engineeringMetabolic Engineering 1:1–11.https://doi.org/10.1006/mben.1998.0101

A theoretical study on the amount of ATP required for synthesis of microbial cell materialAntonie van Leeuwenhoek 39:545–565.https://doi.org/10.1007/BF02578899

A brief history of the development of the KSOM family of mediaJournal of Assisted Reproduction and Genetics 30:995–999.https://doi.org/10.1007/s1081501300978

Metabolic heterogeneity in clonal microbial populationsCurrent Opinion in Microbiology 45:30–38.https://doi.org/10.1016/j.mib.2018.02.004

Increase of mitochondria surrounding spindle causes mouse oocytes arrested at metaphase I stageBiochemical and Biophysical Research Communications 527:1043–1049.https://doi.org/10.1016/j.bbrc.2020.05.037

Fueling thought: Management of glycolysis and oxidative phosphorylation in neuronal metabolismThe Journal of Cell Biology 217:2235–2246.https://doi.org/10.1083/jcb.201803152

Redistribution of mitochondria leads to bursts of ATP production during spontaneous mouse oocyte maturationJournal of Cellular Physiology 224:672–680.https://doi.org/10.1002/jcp.22171
Decision letter

Lydia WS FinleyReviewing Editor; Memorial Sloan Kettering Cancer Center, United States

Naama BarkaiSenior Editor; Weizmann Institute of Science, Israel

Jason W LocasaleReviewer; Duke University School of Medicine, United States

Denis V TitovReviewer; University of California, Berkeley, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
[Editors’ note: the authors submitted for reconsideration following the decision after peer review. What follows is the decision letter after the first round of review.]
Thank you for submitting your work entitled "Coarsegrained model of mitochondrial metabolism enables subcellular flux inference from fluorescence lifetime imaging of NADH" for consideration by eLife. Your article has been reviewed by 4 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by a Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Jason W. Locasale (Reviewer #2); Edmund J. Crampin (Reviewer #3).
We are sorry to say that, after consultation with the reviewers, we have decided that your work will not be considered further for publication by eLife at this time. The reviewers raised three major issues regarding the clarity of the manuscript, the validity of the model, and the potential of the methodology to be used to generate biological insights in oocytes or additional cell systems. There were several concerns about the validity of the assumptions of the model, most notably the concern that as the coarsegrained model describes a substrateenzyme binding reaction with a rate that depends only on the substrate concentration, this model would be valid only under conditions in which enzyme concentrations were orders of magnitude higher than substrate concentration. More broadly, the reviewers raised questions as to whether the assumptions in the model are valid in other systems and whether the methodology will be broadly applicable for the generation of novel biological insights. Given these significant concerns, we are returning the manuscript reviews to you so that they may guide you as you seek publication elsewhere. Should you wish to resubmit at a future date, we would be happy to reconsider a revised version that (1) more clearly justifies the validity of the approach and its assumptions and (2) provides an application of the new method to generate new biological insight.
Reviewer #1:
In this manuscript, the authors present a model to relate FLIM measurements to mitochondrial metabolic fluxes. Using mouse oocytes, which have little NADPH, the authors develop a coarsegrained model to infer mitochondrial NADH oxidation by exploiting NAD(P)H FLIM. Using this approach, the authors uncover regional variation in mitochondrial fluxes in mouse oocytes. The modeled mitochondrial flux shows a strong negative correlation with mitochondrial membrane potential and no correlation with mitochondrial content. While this is not the first paper to use NAD(P)H FLIM to show subcellular metabolic variability, this manuscript does present a model to connect NAD(P)H FLIM to mitochondrial redox cycles. Therefore, the major utility of the model lies in its ability to provide subcellular information about mitochondrial NAD(P)H oxidation. The authors provide a comprehensive and accessible discussion of the assumptions, caveats, and conclusions enabled by their modeling. At present, however, it is not clear to this reviewer how generalizable this method will prove beyond mouse oocytes. This concern stems from the potential difficulty in establishing key parameters of the model in other cell types in which assumptions safely made in mouse oocytes may not be appropriate.
To demonstrate the utility of their model, the authors should test key parameters in at least one additional cell type. In particular, the following issues should be addressed:
1. A key requirement of the model is the ability to determine the equilibrium NADH bound ratio. Here, the authors use low oxygen (or rotenone) to establish this parameter. Will this be feasible in other cell types, for example those with active NNT?
2. The authors note that the confounding signal from NADPH can be ignored in mouse oocytes, which have 40fold higher NADH than NADPH. How generalizable is this? Will other mammalian cell types be amenable to this method?
3. The authors test the assumption that NADH signal originates in the mitochondria by comparing with mitotracker signal under control conditions. This should also be repeated for key conditions (e.g. low oxygen or oxamate treatment).
Reviewer #2:
In the manuscript "Coarsegrained model of mitochondrial metabolism enables subcellular flux inference from fluorescence lifetime imaging of NADH", the authors use fluorescence imaging to estimate NADH/NAD turnover flux and electron transfer rate in the mitochondria of mouse oocytes. Because of high spatial resolution of microscopy, the authors could also observe significant subcellular spatial gradient of oxidative flux in oocytes.
The fluorescence imaging and quantification of flux are generally solid and convincing, but there are issues that need to be addressed.
– In figure 1c, the author estimated two parameters τl and τs under different oxygen levels. They should be constant in all oxygen levels if this model is valid, but they vary a lot when oxygen level is below 10µM. However, the NADH concentration and bound fraction only vary a lot in this oxygen range. This should be addressed.
– The author used a mixture of LDH and NADH to prove the FLIM works in vitro. However, there are lots of different types of enzyme and complex in mitochondria that can bind NADH, and the author's model combines them together to do the calculations. Some justification of this is needed.
– In the spatial model, does variation of thickness of oocyte from center to periphery affect fluorescence levels? If yes, have author corrected this effect and how to correct?
– The nucleus will also cause heterogeneous distribution of mitochondria, which might also need to be considered in modeling spatial distribution.
Reviewer #3:
This paper describes an analysis of fluorescence lifetime imaging (FLIM) of NADH in mitochondria in intact mouse oocytes, using a mathematical model to interpret the fluorescence data to infer mitochondrial NADH redox fluxes. The authors measure FLIM data for varying oxygen concentrations and using several other perturbations to mitochondrial respiration, in order to infer consequential changes to key mitochondrial metabolic fluxes. One striking observation is of subcellular spatial gradients in the inferred metabolic flux across the oocytes.
The authors tackle an important issue in measurement and understanding mitochondrial function in intact cells. The analysis of the FLIM data is dependent on a mathematical model that the authors develop. The correctness and suitability of this model is not clear to me from the way it is described in the manuscript.
The analysis is based on a model, presented in Appendix 2, which considers 'course graining' of a 'detailed' NADH redox model. The latter considers N oxidases and M reductases acting on NADH and NAD+. The aim of the course graining approach is to reduce this model to an equivalent model with one effective oxidase and one effective reductase, and to calculate the effective binding and unbinding coefficients for this reduced model as functions of the binding and unbinding coefficients of the full model. This reduced model is then used to analyse the fluorescence lifetime imaging data, to infer mitochondrial redox fluxes, and to draw some conclusions on spatial gradients of mitochondrial function within the oocytes.
As the analysis and interpretation of the fluorescence lifetime imaging data depends closely on this model, it is necessary to (a) be convinced that the full model is an appropriate representation of the underlying system, and that (b) the course graining methodology is valid.
In relation to (a), the full model (Figure 2 in Appendix 2, but seemingly never presented as a system of mathematical equations), it is not clear why the authors have chosen this particular kinetic scheme, and there is seemingly no specific justification given for it. A simpler scheme would be (for example, for the ith oxidase)
NADH_{f} + Ox_{i} ⇌ Complex ⇌ NAD^{+}_{f} + Ox_{i}
i.e. a standard reversible MichaelisMenten scheme. It isn't clear why the authors have chosen to represent these reactions with two complexes, when a simpler scheme might suffice.
A more serious concern is in relation to (b), the validity of the course graining procedure that is subsequently outlined. The schema that is presented in Appendix 2 has the binding rates to be independent of free enzyme concentration. Thus, in Equation S6, for example, the binding rate for NADH_{f} for the ith oxidase is k_{oi},_{b}∙[NADH_{f}]. But as per the standard analysis of enzyme reactions (Michaelis Menten, etc), following from basic mass action principles, this should depend on the free enzyme concentration, and should instead be given by k_{oi},_{b}∙[NADH_{f}]∙[Ox_{i}] for the ith enzyme. No justification appears to be given as to why there is no dependence on the free enzyme concentration.
Unfortunately, a consequence of this is that the simple factorisation that allows equations S6S9 to simplify into S10 and S11 is no longer possible. Essentially, by omitting the free enzyme concentration the authors end up considering a linear system which can be factored in the way they have presented, whereas the nonlinear system that is obtained when the free enzyme concentration is included does not allow this simple factorisation. It is this nonlinearity that generates the saturating behaviour in standard enzyme kinetics, for example.
Can the authors justify their linear model and omission of the free enzyme concentration? This does not appear to have been justified in the manuscript. One possibility that I can think of may be that the enzymes are all present in very high concentration, such that enzyme concentration is effectively constant. It is not clear that this is an appropriate regime. Alternatively, the authors should justify why the free enzyme concentration has been omitted in these equations. One way to demonstrate the validity of their approach would be through simulation, for example by selecting arbitrary parameters for a model with N and M oxidases and reductases and comparing full simulations of the nonlinear ODE system generated for this model with simulations of the reduced model derived from the course graining approach. As far as I could see there was no demonstration of the validity of the course grained model, however.
Otherwise, if the model is indeed incorrect, it is not clear what is the consequence of this will be on the subsequent data analysis (or indeed if the data can be analysed in this manner), given that Equation 1 in the main text and all subsequent analysis appear to follow directly from this assumption.
I hope that I am wrong, but if I am then I would strongly encourage the authors to provide further justification as to why their model and their course graining approach is correct and valid.
Reviewer #4:
Not every single cell is the same in terms of its metabolism. To study the causes of such celltocell differences, we need microscopic tools to assess metabolic properties, such as metabolite levels and metabolic fluxes, on the single cell level or even beyond. While sensors exist to visualize certain metabolite levels, we still largely lack methods to assess metabolic fluxes in single cells. The work of Yang and Needleman presents a method that can assess under certain assumptions the flux through electron transport chain (ETC) in mitochondria of single mouse oocytes at quasi steadystates with subcellular resolution.
For their method, the authors use FLIM (fluorescence lifetime imaging microscopy) to determine the concentration of free and bound NADH in mitochondria, and these measurements are then used in a simple coarsegrained model to infer the flux through the ETC. This coarsegrained steadystate model describes the oxidation of NADH with one oxidase (resembling the ETC) and one NADH reductase (resembling all the 3 TCA cycle NADH dehydrogenases plus pyruvate dehydrogenase, but neglecting the FADH2dependent succinate dehydrogenase) with only two free model parameters.
Strikingly, when fed with the FLIM data, this coarsegrained model could describe the outcomes of a number of perturbations, where the oxygen uptake rate (i.e. a proxy for the flux through the ETC) was independently measured with a different method. Applying the method, the authors also suggest that the ETC flux is higher in mitochondria that are rather located at the outside of the oocyte.
While FLIM measurements of bound and unbound NADH have been done before, the main strength of the paper is that it presents a method to infer metabolic activity in an oocyte, where the novelty resides on the development of the simple coarsegrained model and on showing that the modelbased analysis of the FLIM data can allow to obtain quasisteadystate ETC fluxes.
The main weakness of the paper is the following: Unfortunately, the work falls short on the application side. One would have wished that for a novel method like this, if it is indeed relevant, it should have been easy for the authors to add exciting application cases that would indeed generate novel biological insight.
While the main strength is the paper is the method (i.e. inference of ETC flux of modelbased analysis of FLIM data), I feel that the description of the method, its assumptions etc falls short, which made assessment of the method and its potential limitations challenging. I feel that this is due to the fact that the writing of the manuscript is suboptimal. While the biochemistry is described/introduced on a very detailed textbook level, the methods, the measurements, the analyses of the measurement data in the result section and in the method section are described in a very short, condensed, and sometime convoluted, manner. As this is primarily meant to be a method paper, the authors need to do a better job in describing what they have done (i.e. model development, model assumptions, inference procedure, etc) in a clearer manner.
I felt that a strong point was that the two different versions of how the experimental data is used in the model, i.e. lifetime (tau) and bound ratio (β), leads to similarly inferred r_{ox}. However, due to the above criticized too short explanations, I could not tell whether this would be trivial or not. Also, the whole method boils down to this equation J_{ox} = α * (β – β_{eq}) * [NADH_{f}], describing the full complexity of mitochondrial metabolism (TCA cycle, the electron transport chain, metabolite exchange between mitochondria and cytoplasm) with a single equation with only two free parameters (α, β_{eq}). For this reviewer, also this part still remains somewhat elusive.
For publication of a new method in a journal such as eLife, I would expect that the manuscript also shows an application of the new method, which generates new biological insight, thereby demonstrating the power and value of the new method. With the measurements of the ETC flux in the spatially located mitochondria the authors start doing this, but unfortunately the manuscript does not develop this into anything interesting (i.e. unraveling of what could cause this gradient, where for instance one could do a dynamic addition of an inhibitor that would diffuse from the outside of the cell, etc). I think it is important to show a compelling case here.
[Editors’ note: further revisions were suggested prior to acceptance, as described below.]
Thank you for submitting your article "A coarsegrained NADH redox model enables inference of subcellular metabolic fluxes from fluorescence lifetime imaging" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Naama Barkai as the Senior Editor. Overall, all 3 reviewers were very positive about your revised manuscript. We would be pleased to accept a revised version that includes the textual modifications suggested by Reviewer #3.
Because we could not reach two of the initial reviewers, we added an additional reviewer at this stage. The following individuals involved in review of your submission have agreed to reveal their identity: Jason W. Locasale (Reviewer #2); Denis V Titov (Reviewer #3).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Please see the suggestions made by Reviewer #3 for clarifying points in the text. No additional experiments  only text revisions  are required.
Reviewer #1:
The authors have significantly improved the manuscript, which is much clearer than before. I think it is very interesting, although I cannot fully evaluate the strength of the assumptions. The authors take great care to explicitly address the major concerns raised by reviewers with regards to the assumptions of the coarsegraining model. Whether these clarified assumptions are valid is not obvious to me, as such modeling is outside of my expertise, but the concern appears to be addressed as far as I can tell. At least, it is clear in the text what assumptions the authors are making.
The finding that ETC flux is not regulated by substrate availability but rather by intrinsic ATP synthesis and proton flux is extremely interesting. While not a criticism of the manuscript, it is worth noting, though, that this result could have been achieved by OCR alone  as the authors show. Therefore, a major strength of this method is the ability to quantify ETC flux with subcellular resolution.
Reviewer #2:
Accept
Reviewer #3:
I was not one of the original reviewers of this manuscript. I have carefully read the manuscript, previous comments by 4 reviewers and authors response. I feel that authors addressed all of the comments raised during original review.
Below I provide comments based on my review of this manuscript, but I want to highlight that I feel that the manuscript is strong in current form and that model is interesting and well validated and that this model will provide interesting insight into activity of ETC in single cells or even subcellular compartment that is difficult to measure currently.
1) Interpretation of Figure 8 data/ homeostasis of ETC flux in MII oocytes. I agree with the authors that Figure 8 clearly demonstrates that various perturbations can lead to changes in redox status of cells without affect Jox/OCR and it is impressive validation of the model that not of these changes affect Jox validated by no effect on OCR. It feels to this reviewer that it is still entirely possible that textbook view that Jox/OCR is regulated by ATP demand is still correct here and it just happened that none of the perturbation that authors used were enough to significantly change ATP consumption in this system and instead all perturbations somehow affected substrate supply without changing demand. Since no measurements of ATP, ADP, AMP have been reported by authors, it is difficult to be sure that any of the perturbation actually affected ATP demand in a significant way although I agree that would have been my expectation. Perhaps, authors should consider adding this caveat to the text that it is possible that none of the perturbations changes ATP demand or if authors strongly believe this hypothesis is correct then showing that large changes in ATP, ADP, AMP are observed without change in Jox/OCR would strongly support this hypothesis.
2) Why does lifetime of free NADH change with low oxygen and with other treatments in other figures)? I can imagine how protein bound lifetime might change as the fractional contributions of proteins to which NADH binds might change with changes in [NADH] but I would have expected free NADH lifetime to stay the same. Perhaps authors should comment on this in the text to clarify this for readers and maybe provide examples of factors that can change free NADH lifetime.
3) Figure 1c. What are the confidence intervals for fittingbased estimates of lifetime and fraction bound of each sample? I could be wrong, but it seems to me that fitting of double exponential equation to decay data might produce redundant values for estimates of lifetimes and fraction bound. I think it would be useful if authors could add confidence intervals (e.g., using bootstrapping with randomly drawn points with substitution from each lifetime curve) that estimate the uncertainty of the fittingbased estimates to show that the fitting procedure used by authors produces nonredundant values of lifetimes and fraction bound. These estimates of uncertainty could also be propagated to Jox to provide a better estimate of Jox uncertainty compared to using one set of values from each FLIM trace that authors use currently if I understood correctly.
https://doi.org/10.7554/eLife.73808.sa1Author response
[Editors’ note: the authors resubmitted a revised version of the paper for consideration. What follows is the authors’ response to the first round of review.]
We are sorry to say that, after consultation with the reviewers, we have decided that your work will not be considered further for publication by eLife at this time. The reviewers raised three major issues regarding the clarity of the manuscript, the validity of the model, and the potential of the methodology to be used to generate biological insights in oocytes or additional cell systems. There were several concerns about the validity of the assumptions of the model, most notably the concern that as the coarsegrained model describes a substrateenzyme binding reaction with a rate that depends only on the substrate concentration, this model would be valid only under conditions in which enzyme concentrations were orders of magnitude higher than substrate concentration. More broadly, the reviewers raised questions as to whether the assumptions in the model are valid in other systems and whether the methodology will be broadly applicable for the generation of novel biological insights. Given these significant concerns, we are returning the manuscript reviews to you so that they may guide you as you seek publication elsewhere. Should you wish to resubmit at a future date, we would be happy to reconsider a revised version that (1) more clearly justifies the validity of the approach and its assumptions and (2) provides an application of the new method to generate new biological insight.
We thank the reviewers for their detailed and constructive comments on our manuscript. We have taken great effort to thoroughly address the issues raised by the reviewers. Our revised manuscript now clearly explains the assumptions of the model, validates the model on an additional cell type and provides two applications of the method to generate new biological insights. We are hereby resubmitting the manuscript for reconsideration for publication at eLife.
Regarding the binding reaction issue raised by Reviewer #3, we clarified that we did not assume a linear reaction model, but kept the model entirely general with all nonlinearities (including in the free enzyme concentrations) incorporated into the kinetic rates. The predicted relationship between FLIM measurements and ETC flux generally holds, independent of these nonlinearities. We introduced Figure 2 and Appendix 2 to illustrate this point and built the model systematically from this generalized enzyme kinetics. To demonstrate the general applicability of the model, we validated the model in an additional cell type: human tissue culture cells (Figure 7). This addresses the concerns from Reviewer #1. In response to the suggestions from Reviewer #4, we provided two applications of the method to gain biological insight: we discovered homeostasis (Figure 8) and subcellular heterogeneities (Figure 9) of ETC flux in mouse oocytes. We concluded from these observations that ETC fluxes in mouse oocytes are not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration. Thus, our work provides new insights into the spatiotemporal regulation of metabolic fluxes in cells.
We addressed all reviewers’ comments itembyitem below and made corresponding changes in the manuscript.
With these significant improvements of the manuscript, we hope our work can be reconsidered for publication at eLife.
Reviewer #1:
In this manuscript, the authors present a model to relate FLIM measurements to mitochondrial metabolic fluxes. Using mouse oocytes, which have little NADPH, the authors develop a coarsegrained model to infer mitochondrial NADH oxidation by exploiting NAD(P)H FLIM. Using this approach, the authors uncover regional variation in mitochondrial fluxes in mouse oocytes. The modeled mitochondrial flux shows a strong negative correlation with mitochondrial membrane potential and no correlation with mitochondrial content. While this is not the first paper to use NAD(P)H FLIM to show subcellular metabolic variability, this manuscript does present a model to connect NAD(P)H FLIM to mitochondrial redox cycles. Therefore, the major utility of the model lies in its ability to provide subcellular information about mitochondrial NAD(P)H oxidation. The authors provide a comprehensive and accessible discussion of the assumptions, caveats, and conclusions enabled by their modeling. At present, however, it is not clear to this reviewer how generalizable this method will prove beyond mouse oocytes. This concern stems from the potential difficulty in establishing key parameters of the model in other cell types in which assumptions safely made in mouse oocytes may not be appropriate.
To demonstrate the utility of their model, the authors should test key parameters in at least one additional cell type. In particular, the following issues should be addressed:
We appreciate the reviewers’ constructive comments. We have now validated our model in an additional cell type: human tissue culture cells. The results are summarized in the new section “The NADH redox model enables accurate prediction of ETC flux in human tissue culture cells” and in Figure 7. These results demonstrated the generality of the method. We address the reviewers’ specific comments item by item below.
1. A key requirement of the model is the ability to determine the equilibrium NADH bound ratio. Here, the authors use low oxygen (or rotenone) to establish this parameter. Will this be feasible in other cell types, for example those with active NNT?
We have demonstrated experimentally that we can also use rotenone to obtain the equilibrium NAD(P)H bound ratio in human tissue culture cells, which leads to accurate prediction of the ETC flux under various mitochondrial inhibitors and nutrient perturbations. The results are summarized in the section “NADH redox model enables accurate prediction of ETC flux in tissue culture cells” and in Figure 7 and Figure 7—figure supplement 1.
2. The authors note that the confounding signal from NADPH can be ignored in mouse oocytes, which have 40fold higher NADH than NADPH. How generalizable is this? Will other mammalian cell types be amenable to this method?
Concentration of NADPH can be of the same order as NADH in other cell types, including mammalian tissue culture cells. Therefore, we have now explicitly accounted for the impact of NADPH signal and other fluorescence background signal on the flux inference procedure. We showed that if the background fluorescence is either constant under perturbations or changes proportionally with NADH, then it does not impact the inference procedure. This derivation is presented in Appendix 5 section “Accounting for NADPH and other background fluorescence”. We further tested this method experimentally in human tissue culture cells, in which NADPH is substantial, and accurately predicted the change of ETC flux under a variety of metabolic perturbations. Our results demonstrated that the presented flux inference method can work on other cell types, despite of the potential enrichment of NADPH.
3. The authors test the assumption that NADH signal originates in the mitochondria by comparing with mitotracker signal under control conditions. This should also be repeated for key conditions (e.g. low oxygen or oxamate treatment).
We repeated the test for oxamate, oligomycin, fccp and rotenone conditions, which yielded a NADH signal originating from mitochondria of 78.6% ± 1.4%, 84.1% ± 1.6%, 83.7% ± 0.5%, 81.7% ± 2.0%, respectively, similar to the control condition. We added these data to the manuscript.
Reviewer #2:
In the manuscript "Coarsegrained model of mitochondrial metabolism enables subcellular flux inference from fluorescence lifetime imaging of NADH", the authors use fluorescence imaging to estimate NADH/NAD turnover flux and electron transfer rate in the mitochondria of mouse oocytes. Because of high spatial resolution of microscopy, the authors could also observe significant subcellular spatial gradient of oxidative flux in oocytes.
The fluorescence imaging and quantification of flux are generally solid and convincing, but there are issues that need to be addressed.
We appreciate the positive comments of the reviewer. We address the issues below.
– In figure 1c, the author estimated two parameters τl and τs under different oxygen levels. They should be constant in all oxygen levels if this model is valid, but they vary a lot when oxygen level is below 10µM. However, the NADH concentration and bound fraction only vary a lot in this oxygen range. This should be addressed.
The fluorescence lifetimes $\tau}_{\mathrm{l}$ and $\tau}_{\mathrm{s}$ do not have to be constant for the model to be valid. In fact, we have shown that the variation of $\tau}_{\mathrm{l}$ provides an additional method, in a selfconsistent way with the bound ratio method, to infer ETC flux. Under this framework, the variations of $\tau}_{\mathrm{l}$ could be attributed to the change of the relative abundance of bound NADH to the oxidases and reductases in the cell (Equation 6). We have now included the derivations and experimental validations in the section “Accurately predicting ETC flux from FLIM of NADH using the NADH redox model” and in Appendix 5.
– The author used a mixture of LDH and NADH to prove the FLIM works in vitro. However, there are lots of different types of enzyme and complex in mitochondria that can bind NADH, and the author's model combines them together to do the calculations. Some justification of this is needed.
FLIM is well established to be able to distinguish bound and free NADH in vivo based on the large change of fluorescence lifetime when NADH binds to enzymes (Bird et al., 2005; Skala et al., 2007; Heikal, 2010; Sharick et al., 2018; Sanchez et al., 2019). Even though how much lifetime changes depends on the specific enzyme NADH binds to, the enzymebound NADH always has a much longer fluorescence lifetime than free NADH, so a two exponential fitting of the fluorescence decay curve will yield an average lifetime and fraction of NADH bound to all enzymes in vivo. Therefore, the method to calculate free and bound concentrations of NADH from FLIM measurements is expected to hold in vivo. We have now clarified this point in the manuscript.
– In the spatial model, does variation of thickness of oocyte from center to periphery affect fluorescence levels? If yes, have author corrected this effect and how to correct?
We have shown that the intensity gradient is present for NADH and TMRM, but not for
Mitotracker (Figure 9—figure supplement 2), hence ruling out the impact of thickness variations.
– The nucleus will also cause heterogeneous distribution of mitochondria, which might also need to be considered in modeling spatial distribution.
We are using oocytes arrested at Meiosis II, which do not have a nucleus. The presence of nucleus in other cell types will indeed need to be considered in studying spatial distribution of metabolism.
Reviewer #3:
This paper describes an analysis of fluorescence lifetime imaging (FLIM) of NADH in mitochondria in intact mouse oocytes, using a mathematical model to interpret the fluorescence data to infer mitochondrial NADH redox fluxes. The authors measure FLIM data for varying oxygen concentrations and using several other perturbations to mitochondrial respiration, in order to infer consequential changes to key mitochondrial metabolic fluxes. One striking observation is of subcellular spatial gradients in the inferred metabolic flux across the oocytes.
The authors tackle an important issue in measurement and understanding mitochondrial function in intact cells. The analysis of the FLIM data is dependent on a mathematical model that the authors develop. The correctness and suitability of this model is not clear to me from the way it is described in the manuscript.
The analysis is based on a model, presented in Appendix 2, which considers 'course graining' of a 'detailed' NADH redox model. The latter considers N oxidases and M reductases acting on NADH and NAD+. The aim of the course graining approach is to reduce this model to an equivalent model with one effective oxidase and one effective reductase, and to calculate the effective binding and unbinding coefficients for this reduced model as functions of the binding and unbinding coefficients of the full model. This reduced model is then used to analyse the fluorescence lifetime imaging data, to infer mitochondrial redox fluxes, and to draw some conclusions on spatial gradients of mitochondrial function within the oocytes.
As the analysis and interpretation of the fluorescence lifetime imaging data depends closely on this model, it is necessary to (a) be convinced that the full model is an appropriate representation of the underlying system, and that (b) the course graining methodology is valid.
In relation to (a), the full model (Figure 2 in Appendix 2, but seemingly never presented as a system of mathematical equations), it is not clear why the authors have chosen this particular kinetic scheme, and there is seemingly no specific justification given for it. A simpler scheme would be (for example, for the ith oxidase)
NADH_{f} + Ox_{i} ⇌ Complex ⇌ NAD^{+}_{f} + Ox_{i}
i.e. a standard reversible MichaelisMenten scheme. It isn't clear why the authors have chosen to represent these reactions with two complexes, when a simpler scheme might suffice.
A more serious concern is in relation to (b), the validity of the course graining procedure that is subsequently outlined. The schema that is presented in Appendix 2 has the binding rates to be independent of free enzyme concentration. Thus, in Equation S6, for example, the binding rate for NADH_{f} for the ith oxidase is k_{oi},_{b}∙[NADH_{f}]. But as per the standard analysis of enzyme reactions (Michaelis Menten, etc), following from basic mass action principles, this should depend on the free enzyme concentration, and should instead be given by k_{oi},_{b}∙[NADH_{f}]∙[Ox_{i}] for the ith enzyme. No justification appears to be given as to why there is no dependence on the free enzyme concentration.
Unfortunately, a consequence of this is that the simple factorisation that allows equations S6S9 to simplify into S10 and S11 is no longer possible. Essentially, by omitting the free enzyme concentration the authors end up considering a linear system which can be factored in the way they have presented, whereas the nonlinear system that is obtained when the free enzyme concentration is included does not allow this simple factorisation. It is this nonlinearity that generates the saturating behaviour in standard enzyme kinetics, for example.
Can the authors justify their linear model and omission of the free enzyme concentration? This does not appear to have been justified in the manuscript. One possibility that I can think of may be that the enzymes are all present in very high concentration, such that enzyme concentration is effectively constant. It is not clear that this is an appropriate regime. Alternatively, the authors should justify why the free enzyme concentration has been omitted in these equations. One way to demonstrate the validity of their approach would be through simulation, for example by selecting arbitrary parameters for a model with N and M oxidases and reductases and comparing full simulations of the nonlinear ODE system generated for this model with simulations of the reduced model derived from the course graining approach. As far as I could see there was no demonstration of the validity of the course grained model, however.
Otherwise, if the model is indeed incorrect, it is not clear what is the consequence of this will be on the subsequent data analysis (or indeed if the data can be analysed in this manner), given that Equation 1 in the main text and all subsequent analysis appear to follow directly from this assumption.
I hope that I am wrong, but if I am then I would strongly encourage the authors to provide further justification as to why their model and their course graining approach is correct and valid.
We thank the reviewer for the detailed and indepth comments. We agree that we did not do a good job of describing the model in the original manuscript, which caused confusion regarding the notation and assumptions underlying the model. One of the sources of confusion is our usage of a reduced notation which was not adequately explained. We have now introduced Figure 2 together with explicit mathematical descriptions (Appendix 2) to explain the reduced notation.
Because of the use of our reduced notation, it might seem that we have assumed that the binding rate is linear in NADH concentration and does not depend on the free enzyme concentration and other factors. However, this is not the case. In our reduced notation, these dependencies are embedded in the rates themselves. As a result, all of the rates can depend on many variables including the free enzyme concentration (such as [Ox_{i}]), NADH concentration and other factors. These additional variables will obey their own dynamical equations. We used this reduced notation because we want to capture a broad class of models, including models where the rates have nonlinear dependencies on these variables. Modeling the dynamics of the redox pathways would require specify these dependencies: i.e, writing down the explicit functional form through which the rates depend on these variables and dynamical equations for the temporal variation of these variables. Remarkably, the relationship between FLIM parameters and the predicted ETC flux (Equations 5ac) does not depend on these modeling details. Thus, using the reduced notation allows us to derive these general results without specifying the detailed enzyme kinetics. For clarity, we now also explicitly show that our general model includes the reversible MichaelisMenten model as a special case (Appendix 7).
The coarsegraining procedure is mathematically exact and independent of the functional forms of the kinetic rates. In our reduced notation, all the coarsegrained parameters can also have complex dependencies on free enzyme concentrations, NADH concentrations, etc. To completely describe the dynamics of the coarsegrained model would also require specifying these dependencies and the associated dynamical equations. However, they do not affect the predicted relationship between FLIM parameters and ETC flux. As an example, we have also explicitly demonstrated the coarsegraining for a redox cycle in which each of the enzymes obey reversible MichaelisMenten kinetics (Appendix 7).
The reason our generalized model works is because the predicted flux only depends on two coarsegrained parameters, $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$, which are experimentally determined. The detailed kinetics of a specific enzyme model only modifies the functional forms of $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$. For example, we have explicitly shown the functional dependence of $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$ on the parameters of the reversible MichaelisMenten kinetics (Appendix 7, Equations S73S74). Since $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$ can be measured experimentally with FLIM of NADH, they can be used for the flux inference without knowing their functional dependencies.
For our inference method to work, we assumed that the NADH redox loop is at steady state. This assumption is sufficient to predict the relationship between FLIM parameters and the ETC flux (Equations 5ac). For this relationship to be useful for flux inference, we need to be able to determine $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$ reliably from experiments. To determine $\beta}_{\mathrm{e}\mathrm{q}$, we need to perform an experiment to inhibit ETC flux without changing the value of $\beta}_{\mathrm{e}\mathrm{q}$. In our work, we chose oxygen drop experiment or rotenone inhibition for this purpose. We also assumed that $\alpha$ is a constant across perturbations, and hence we determined $\alpha$ in a control condition and used the same value across all other perturbations. Thus, the three key assumptions are: (1) that the NADH redox loop can be well approximated as being at steadystate; (2) that $\beta}_{\mathrm{e}\mathrm{q}$ does not vary due to the perturbation used to inhibit ETC flux and measure $\beta}_{\mathrm{e}\mathrm{q}$; (3) that $\alpha$ is constant across perturbations. It is an empirical question if these assumptions hold. We tested our model predictions by comparing the predicted ETC flux with direct measurements of oxygen consumption rate (OCR) across a wide range of perturbations for both mouse oocytes (Figure 5, Figure 8) and human tissue culture cells (Figure 7). The quantitative agreement between the predicted flux and OCR robustly validated our model assumptions.
With the improved model presentation and additional experimental validation, we hope we have convinced the reviewer of the correctness and suitability of our method. We thank the reviewer again for giving us the opportunity to significantly improve our manuscript.
Reviewer #4:
Not every single cell is the same in terms of its metabolism. To study the causes of such celltocell differences, we need microscopic tools to assess metabolic properties, such as metabolite levels and metabolic fluxes, on the single cell level or even beyond. While sensors exist to visualize certain metabolite levels, we still largely lack methods to assess metabolic fluxes in single cells. The work of Yang and Needleman presents a method that can assess under certain assumptions the flux through electron transport chain (ETC) in mitochondria of single mouse oocytes at quasi steadystates with subcellular resolution.
For their method, the authors use FLIM (fluorescence lifetime imaging microscopy) to determine the concentration of free and bound NADH in mitochondria, and these measurements are then used in a simple coarsegrained model to infer the flux through the ETC. This coarsegrained steadystate model describes the oxidation of NADH with one oxidase (resembling the ETC) and one NADH reductase (resembling all the 3 TCA cycle NADH dehydrogenases plus pyruvate dehydrogenase, but neglecting the FADH2dependent succinate dehydrogenase) with only two free model parameters.
Strikingly, when fed with the FLIM data, this coarsegrained model could describe the outcomes of a number of perturbations, where the oxygen uptake rate (i.e. a proxy for the flux through the ETC) was independently measured with a different method. Applying the method, the authors also suggest that the ETC flux is higher in mitochondria that are rather located at the outside of the oocyte.
While FLIM measurements of bound and unbound NADH have been done before, the main strength of the paper is that it presents a method to infer metabolic activity in an oocyte, where the novelty resides on the development of the simple coarsegrained model and on showing that the modelbased analysis of the FLIM data can allow to obtain quasisteadystate ETC fluxes.
The main weakness of the paper is the following: Unfortunately, the work falls short on the application side. One would have wished that for a novel method like this, if it is indeed relevant, it should have been easy for the authors to add exciting application cases that would indeed generate novel biological insight.
While the main strength is the paper is the method (i.e. inference of ETC flux of modelbased analysis of FLIM data), I feel that the description of the method, its assumptions etc falls short, which made assessment of the method and its potential limitations challenging. I feel that this is due to the fact that the writing of the manuscript is suboptimal. While the biochemistry is described/introduced on a very detailed textbook level, the methods, the measurements, the analyses of the measurement data in the result section and in the method section are described in a very short, condensed, and sometime convoluted, manner. As this is primarily meant to be a method paper, the authors need to do a better job in describing what they have done (i.e. model development, model assumptions, inference procedure, etc) in a clearer manner.
I felt that a strong point was that the two different versions of how the experimental data is used in the model, i.e. lifetime (tau) and bound ratio (β), leads to similarly inferred r_{ox}. However, due to the above criticized too short explanations, I could not tell whether this would be trivial or not. Also, the whole method boils down to this equation J_{ox} = α * (β – β_{eq}) * [NADH_{f}], describing the full complexity of mitochondrial metabolism (TCA cycle, the electron transport chain, metabolite exchange between mitochondria and cytoplasm) with a single equation with only two free parameters (α, β_{eq}). For this reviewer, also this part still remains somewhat elusive.
We thank the reviewer for the detailed and indepth review of our manuscript. We appreciate the reviewer’s suggestion to add application cases to demonstrate the usefulness of our method. We now added two application of our flux inference procedure to the revised manuscript. The first case is the discovery of homeostasis of ETC flux in mouse oocytes: perturbations of nutrient supply and energy demand do not change ETC flux despite significantly impacting NADH metabolic state (Figure 8). The second case is the discovery of the intracellular spatial gradient of ETC flux in mouse oocytes. As suggested by the reviewer, we have used metabolic inhibitors to help reveal the cause of this gradient and found that this gradient is primarily a result of a spatially heterogeneous mitochondrial proton leak (Figure 9). We concluded from these observations that ETC flux in mouse oocytes is not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration.
We thank the reviewer for their suggestions to improve the presentation of this work. We have significantly rewritten the paper to clearly describe the model development, model assumptions, data analysis procedures and results. Regarding the comparison of the two inference methods, we presented details of the assumptions and derivations in the Results section and demonstrated that the agreement between these two methods is not trivial, and is a robust selfconsistency check of the method. We also now explained the coarsegraining procedure in detail in the main text and in Appendix 2 and 3 to demonstrate how all the model complexities are coarsegrained into only two free parameters $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$. In a nutshell, $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$ would be functions of the kinetic rates of the model, with the kinetic rates depending on the details of mitochondrial metabolism. However, using the model to infer ETC flux does not require knowing the functional forms of $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$, because $\alpha$ and $\beta}_{\mathrm{e}\mathrm{q}$ can be experimentally measured with FLIM. The only assumptions required are that $\alpha$ remains a constant under perturbations and $\beta}_{\mathrm{e}\mathrm{q}$ can be determined from ETC inhibitions. These assumptions are validated experimentally in mouse oocytes and human tissue culture cells from the agreement between predicted ETC flux and direct measurements of OCR.
We also validated our model in an additional cell type of human tissue culture cells, demonstrating the generality of our method (Figure 7).
For publication of a new method in a journal such as eLife, I would expect that the manuscript also shows an application of the new method, which generates new biological insight, thereby demonstrating the power and value of the new method. With the measurements of the ETC flux in the spatially located mitochondria the authors start doing this, but unfortunately the manuscript does not develop this into anything interesting (i.e. unraveling of what could cause this gradient, where for instance one could do a dynamic addition of an inhibitor that would diffuse from the outside of the cell, etc). I think it is important to show a compelling case here.
We have now strengthened the paper with two compelling cases of application: homeostasis and heterogeneity of ETC flux. We have shown that perturbations of nutrient supply and energy demand impact NADH metabolic state but do not impact ETC flux in mouse oocytes (Figure 8). We have also unraveled the cause of the ETC flux gradient in the oocyte using oligomycin and discovered that the flux gradient is caused by a higher rate of proton leak of mitochondria closer to the cell periphery (Figure 9). We concluded from these observations that ETC flux in mouse oocytes is not controlled by energy demand or supply, but by the intrinsic rates of mitochondrial respiration. Our work provided new insights into the spatiotemporal regulations of metabolic fluxes in cells.
[Editors’ note: what follows is the authors’ response to the second round of review.]
Please see the suggestions made by Reviewer #3 for clarifying points in the text. No additional experiments  only text revisions  are required.
Reviewer #3:
I was not one of the original reviewers of this manuscript. I have carefully read the manuscript, previous comments by 4 reviewers and authors response. I feel that authors addressed all of the comments raised during original review.
Below I provide comments based on my review of this manuscript, but I want to highlight that I feel that the manuscript is strong in current form and that model is interesting and well validated and that this model will provide interesting insight into activity of ETC in single cells or even subcellular compartment that is difficult to measure currently.
We thank the reviewer for the positive comments and detailed suggestions on our manuscript.
1) Interpretation of Figure 8 data/ homeostasis of ETC flux in MII oocytes. I agree with the authors that Figure 8 clearly demonstrates that various perturbations can lead to changes in redox status of cells without affect Jox/OCR and it is impressive validation of the model that not of these changes affect Jox validated by no effect on OCR. It feels to this reviewer that it is still entirely possible that textbook view that Jox/OCR is regulated by ATP demand is still correct here and it just happened that none of the perturbation that authors used were enough to significantly change ATP consumption in this system and instead all perturbations somehow affected substrate supply without changing demand. Since no measurements of ATP, ADP, AMP have been reported by authors, it is difficult to be sure that any of the perturbation actually affected ATP demand in a significant way although I agree that would have been my expectation. Perhaps, authors should consider adding this caveat to the text that it is possible that none of the perturbations changes ATP demand or if authors strongly believe this hypothesis is correct then showing that large changes in ATP, ADP, AMP are observed without change in Jox/OCR would strongly support this hypothesis.
We thank the reviewer for the insightful comments. We agree that while NADH metabolic state significantly changed in response to perturbing energy demand and supply, indicating cell metabolism was indeed impacted, it is unclear if these perturbations also influenced ATP, ADP or AMP levels. Future work, including direct measurements of ATP, ADP and AMP levels, will be required to uncover the mechanism of flux homeostasis. We added texts in the discussion (pg. 26) to clarify this point.
2) Why does lifetime of free NADH change with low oxygen and with other treatments in other figures)? I can imagine how protein bound lifetime might change as the fractional contributions of proteins to which NADH binds might change with changes in [NADH] but I would have expected free NADH lifetime to stay the same. Perhaps authors should comment on this in the text to clarify this for readers and maybe provide examples of factors that can change free NADH lifetime.
We have clarified in the text (pg. 4) that in addition to enzyme binding (Sharick et al., 2018, Ghukasyan and Heikal, 2015), NADH fluorescence lifetimes are also impacted by other factors including viscosity and pH (Ghukasyan and Heikal, 2015). We expect that these latter effects might be responsible for the observed changes in the lifetime of free NADH.
3) Figure 1c. What are the confidence intervals for fittingbased estimates of lifetime and fraction bound of each sample? I could be wrong, but it seems to me that fitting of double exponential equation to decay data might produce redundant values for estimates of lifetimes and fraction bound. I think it would be useful if authors could add confidence intervals (e.g., using bootstrapping with randomly drawn points with substitution from each lifetime curve) that estimate the uncertainty of the fittingbased estimates to show that the fitting procedure used by authors produces nonredundant values of lifetimes and fraction bound. These estimates of uncertainty could also be propagated to Jox to provide a better estimate of Jox uncertainty compared to using one set of values from each FLIM trace that authors use currently if I understood correctly.
We thank the reviewer for the suggestions. In this manuscript, we independently fit FLIM curves for each individual oocyte. The reported error bars in this manuscript are standard errors of the mean (SEMs) across these measurements, which depends on the level of variation (the standard deviation) between the oocytes. Two sources of variation in FLIM measurements across the oocytes are: (1) true biological variations between oocytes and (2) fitting errors in the FLIM analysis. To estimate the error of fitting, we performed bootstrapping with randomly drawn points with substitution from each fluorescence decay curve for 53 oocytes. There are ~66000 photons per oocyte, from which we generated 10 bootstrapped decay curves per oocyte to estimate the fitting error. The fitting error is computed as the variance and covariance of the fitted parameters across bootstrapped decay curves and averaged over 53 oocytes.
For oocytes at high oxygen levels in AKSOM, the bootstrapping yields a variance of 2.2×10^{4}, 4.6×10^{3} ns^{2}, 6.0×10^{4} ns^{2} for bound fraction, long lifetime and short lifetime, respectively. The celltocell variances obtained from a single fit per oocyte are 4.4×10^{4}, 9.5×10^{3} ns^{2}, 1.6×10^{3} ns^{2} for bound fraction, long lifetime, short lifetime, respectively. Hence fitting errors account for 50%, 49% and 40% of the celltocell variance in bound fraction, long lifetime and short lifetime, respectively.
The inferred mean flux for oocytes at high oxygen levels in AKSOM is $\u27e8{J}_{\text{ox}}\u27e9=56.6\text{}\mu M\cdot {s}^{1}$. Propagating the error of fitting in all parameters from the bootstrapping analysis to the inferred flux gives a standard error of the mean in $J}_{\text{ox}$ of 1.1 µM·s^{1}. The standard error of the mean in $J}_{\text{ox}$ obtained from a single fit per oocyte was 2.0 µM·s^{1}. Thus, fitting errors account for ∼50% of the standard error of the mean in $J}_{\text{ox}$.
Reviewer 3 is correct that fitting errors can induce covariations between long lifetime and fraction bound: our bootstrapping analysis indicates that with ~66000 photons per oocyte this covariance is of 1.0×10^{3} ns. We believe that such fitting errors do not significantly impact our results for two reasons: (1) The level of covariation between these two parameters in this study is much greater than can be accounted for by such errors. For example, the covariance between long lifetime and bound fraction during oxygen drop experiments is 4.6×10^{3} ns, so fitting errors can only account for ~20% of this covariation; (2) Under other conditions, such as during mouse preimplantation embryo development, we observe that long lifetime and bound fraction are positively correlated during certain developmental stages (Sanchez et al. 2019).
We have added a section of error analysis in the methods section of the manuscript (pg. 30).
https://doi.org/10.7554/eLife.73808.sa2Article and author information
Author details
Funding
National Institutes of Health (R01HD09255001)
 Dan Needleman
National Science Foundation (PFITT1827309)
 Dan Needleman
National Science Foundation (PHY2013874)
 Dan Needleman
National Science Foundation (MCB2052305)
 Dan Needleman
National Science Foundation (PHY1748958)
 Xingbo Yang
 Dan Needleman
National Institutes of Health (R25GM067110)
 Xingbo Yang
 Dan Needleman
Gordon and Betty Moore Foundation (2919.02)
 Xingbo Yang
 Dan Needleman
National Science Foundation (1764269)
 Gloria Ha
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank Easun Arunachalam, YuChen Chao, Will Conway, Carlos Manlio DíazGarcía, Peter Foster, Bill Ireland, Denis Titov, and Gary Yellen for suggestions, advice, and comments on the manuscript. This work is supported by the National Institutes of Health (R01HD09255001) and the National Science Foundation (PFITT1827309, PHY2013874, and MCB2052305). G.H. acknowledges support from the NSFSimons Center for Mathematical and Statistical Analysis of Biology at Harvard (Award number #1764269) and the Harvard Quantitative Biology Initiative. X.Y. and D.J.N. acknowledge discussions with participants of the ‘Cellenergy19’ and the ‘Active20’ KITP programs, supported in part by the National Science Foundation Grant no. NSF PHY1748958, NIH Grant no. R25GM067110, and the Gordon and Betty Moore Foundation Grant no. 2919.02.
Senior Editor
 Naama Barkai, Weizmann Institute of Science, Israel
Reviewing Editor
 Lydia WS Finley, Memorial Sloan Kettering Cancer Center, United States
Reviewers
 Jason W Locasale, Duke University School of Medicine, United States
 Denis V Titov, University of California, Berkeley, United States
Publication history
 Preprint posted: November 21, 2020 (view preprint)
 Received: September 11, 2021
 Accepted: November 19, 2021
 Accepted Manuscript published: November 22, 2021 (version 1)
 Version of Record published: March 17, 2022 (version 2)
Copyright
© 2021, Yang 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

 1,409
 Page views

 263
 Downloads

 2
 Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
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)
Further reading

 Biochemistry and Chemical Biology
 Structural Biology and Molecular Biophysics
Selfcleaving ribozymes are RNA molecules that catalyze the cleavage of their own phosphodiester backbones. These ribozymes are found in all domains of life and are also a tool for biotechnical and synthetic biology applications. Selfcleaving ribozymes are also an important model of sequencetofunction relationships for RNA because their small size simplifies synthesis of genetic variants and selfcleaving activity is an accessible readout of the functional consequence of the mutation. Here, we used a highthroughput experimental approach to determine the relative activity for every possible single and double mutant of five selfcleaving ribozymes. From this data, we comprehensively identified nonadditive effects between pairs of mutations (epistasis) for all five ribozymes. We analyzed how changes in activity and trends in epistasis map to the ribozyme structures. The variety of structures studied provided opportunities to observe several examples of common structural elements, and the data was collected under identical experimental conditions to enable direct comparison. Heatmapbased visualization of the data revealed patterns indicating structural features of the ribozymes including paired regions, unpaired loops, noncanonical structures, and tertiary structural contacts. The data also revealed signatures of functionally critical nucleotides involved in catalysis. The results demonstrate that the data sets provide structural information similar to chemical or enzymatic probing experiments, but with additional quantitative functional information. The largescale data sets can be used for models predicting structure and function and for efforts to engineer selfcleaving ribozymes.

 Biochemistry and Chemical Biology
 Cell Biology
Dystroglycan (DG) requires extensive posttranslational processing and Oglycosylation to function as a receptor for extracellular matrix (ECM) proteins containing lamininGlike (LG) domains. Matriglycan is an elongated polysaccharide of alternating xylose (Xyl) and glucuronic acid (GlcA) that binds with highaffinity to ECM proteins with LGdomains and is uniquely synthesized on αdystroglycan (αDG) by likeacetylglucosaminyltransferase1 (LARGE1). Defects in the posttranslational processing or Oglycosylation of αDG that result in a shorter form of matriglycan reduce the size of αDG and decrease laminin binding, leading to various forms of muscular dystrophy. Previously, we demonstrated that Protein OMannose Kinase (POMK) is required for LARGE1 to generate fulllength matriglycan on αDG (~150250 kDa) (Walimbe et al., 2020). Here, we show that LARGE1 can only synthesize a short, nonelongated form of matriglycan in mouse skeletal muscle that lacks the DG Nterminus (αDGN), resulting in a ~100125 kDa αDG. This smaller form of αDG binds laminin and maintains specific force but does not prevent muscle pathophysiology, including reduced force production after eccentric contractions or abnormalities in the neuromuscular junctions. Collectively, our study demonstrates that αDGN, like POMK, is required for LARGE1 to extend matriglycan to its full mature length on αDG and thus prevent muscle pathophysiology.