Atomic structure of a mitochondrial complex I intermediate from vascular plants
Abstract
Respiration, an essential metabolic process, provides cells with chemical energy. In eukaryotes, respiration occurs via the mitochondrial electron transport chain (mETC) composed of several large membrane-protein complexes. Complex I (CI) is the main entry point for electrons into the mETC. For plants, limited availability of mitochondrial material has curbed detailed biochemical and structural studies of their mETC. Here, we present the cryoEM structure of the known CI assembly intermediate CI* from Vigna radiata at 3.9 Å resolution. CI* contains CI’s NADH-binding and CoQ-binding modules, the proximal-pumping module and the plant-specific γ-carbonic-anhydrase domain (γCA). Our structure reveals significant differences in core and accessory subunits of the plant complex compared to yeast, mammals and bacteria, as well as the details of the γCA domain subunit composition and membrane anchoring. The structure sheds light on differences in CI assembly across lineages and suggests potential physiological roles for CI* beyond assembly.
eLife digest
Respiration is the process used by all forms of life to turn organic matter from food into energy that cells can use to live and grow. The final stage of this process relies on an intricate chain of protein complexes which produce the molecule that cells use for energy. Complexes in the chain are made up of specific proteins that are carefully assembled, often into discrete modules or intermediate complexes, before coming together to form the full protein complex. Understanding how these complexes are assembled provides important insights into how respiration works.
The precise three-dimensional structure of these complexes has been identified for bacteria, yeast and mammals. However, less is known about how these respiration complexes form in plants. For this reason, Maldonado et al. studied the structure of an intermediate complex that is only found in plants, called Cl*. This intermediate structure goes on to form complex I – the largest complex in the respiration chain.
A technique called cryo-electron microscopy was used to obtain a structure of Cl* at a near-atomic level of detail. This structure revealed how the proteins that make up Cl* fit together, highlighting differences and similarities in how plants assemble complex I compared to bacteria, yeast and mammals. Maldonado et al. also studied the activity of Cl*, leading to the suggestion that this complex may be more than just a stepping stone towards building the full complex I and could have its own role in the cell.
The structure of this complex provides new insights into the respiration mechanism of plants and could help scientists improve crop production. For instance, new compounds may be able to block respiration in pests, while leaving the crop unharmed; or genetic modifications could create plants that respire more efficiently in different environments.
Introduction
Respiration is an essential metabolic process that provides the energy and intermediate metabolites needed for growth and maintenance of all eukaryotes. In plants, respiratory pathways are not only involved in energy conversion but also play crucial roles in the procurement of biosynthetic precursors and in the balancing of the cellular redox state (O'Leary et al., 2019). Plant respiratory processes are also closely intertwined with photosynthetic pathways. Despite the importance of respiratory processes to plants’ biomass accumulation, carbon flux and acclimation (O'Leary et al., 2019; Amthor et al., 2019; Heskel et al., 2016), the fundamental mechanisms by which the plant mitochondrial electron transport chain (mETC) produces proton (H+) gradients that are converted into chemical energy remain poorly understood. Molecular knowledge of the structures and mechanisms of the plant mETC components, which differ significantly in their assembly and composition from better-studied mammalian systems, is essential to understand how plants efficiently convert energy and balance respiration with photosynthesis.
Plant mitochondria possess a ‘canonical’ mETC shared with most eukaryotes that is composed of four large membrane protein complexes (complexes I-IV, CI-IV) and an associated ATP synthase in the inner mitochondrial membrane (IMM). Complexes I-IV couple oxidoreduction reactions to H+ pumping against the concentration gradient across the IMM to produce a large H+ electrochemical potential (‘proton motive force’) that is then dissipated through ATP synthase’s rotary mechanism to produce ATP in the mitochondrial matrix. Additionally, plants also possess an ‘alternative’ mETC that dissipates reduction equivalents in a non-H+-pumping, non-energy-conserving fashion (Millar et al., 2011; Schertl and Braun, 2014).
Complex I (CI) is the main energy-conserving entry point for electrons into the mETC. In plants, as in most eukaryotes so far studied, CI is the largest (~1 MDa) and mechanistically least understood component of the mETC (Sazanov, 2015; Hirst, 2013). CI oxidizes NADH and reduces coenzyme Q (CoQ, ubiquinone), pumping four H+ per two electrons from NADH (Jones et al., 2017). CI is an L-shaped multiprotein complex, with a membrane arm and a peripheral arm. In eukaryotes, the peripheral arm of CI extends into the mitochondrial matrix, while the membrane arm is buried within the IMM. Both arms are composed of ‘modules’ with specific functions and distinct evolutionary origins (Efremov and Sazanov, 2012). The peripheral arm contains the NADH dehydrogenase N-module and the CoQ-reducing Q-module, which provide the binding sites for NADH and quinone, respectively, as well as the chain of FeS clusters needed for electron transfer (Figure 1A). The membrane arm contains four proton pumps, two of which are located in the proximal-pumping module (PP), with the remaining two pumps in the distal-pumping module (PD; Figure 1A; Dröse et al., 2011). Through a still poorly understood mechanism, the energy released from NADH-CoQ oxidoreduction in the peripheral arm (N- and Q-modules) is coupled to conformational changes along the membrane arm (PP and PD), resulting in proton pumping from the mitochondrial matrix into the mitochondrial intermembrane space (IMS).

The structure of CI* from Vigna radiata.
(A) An overview of the conserved modular structure of CI using the Thermus thermophilus bacterial core subunits as a simple model (PDB: 4HEA) (Baradaran et al., 2013). (B) CryoEM density map of CI* from V. radiata highlighting its modular architecture. N, NADH-binding module; Q, quinone-binding module; PP, proximal-pump module; PD, distal-pump module; γCA, carbonic anhydrase domain, see also Video 1). (C) Atomic model of V. radiata CI* with all 30 assigned subunits labeled. The additional N-terminal helix of NDUS8 is indicated with an asterisk (*).
Across the studied eukaryotes, mitochondrial CI is composed of 14 highly conserved ‘core’ subunits that are responsible for electron transport and H+ pumping, and 30–35 ‘accessory’ subunits that are involved in CI’s assembly, stability and regulation (Millar et al., 2011; Meyer, 2012). The exact number of subunits in plant mitochondrial CI is still unclear, with several mass spectrometry measurements revealing differing compositions (Meyer, 2012). Nonetheless, it is known that several plant CI accessory subunits are not found in fungi and metazoans (opisthokonts). Most notably, five gamma-type carbonic anhydrase (γCA) proteins (CA1, CA2, CA3, CAL1, and CAL2) have been shown to be associated with CI in plants (Sunderhaus et al., 2006; Perales et al., 2004). These proteins are located on the matrix side of CI’s membrane arm, likely as a heterotrimer of CAL1 or CAL2 monomer plus a CA1/CA2 hetero- or homodimer (Fromm et al., 2016). Hence, only a subset of the five γCA proteins are expected to be simultaneously associated with CI. Although the exact γCA protein combinations are likely tissue- and development-stage-dependent (Cï Rdoba et al., 2019), the role of the γCA domain in plant CI’s function is unknown (Martin et al., 2009).
Another major difference between plants and metazoans occurs in the CI assembly pathway. In metazoans, the N-module (which is responsible for NADH oxidation) is assembled onto the rest of the complex (Q-, PP- and PD-modules) as the final step of assembly (Formosa et al., 2018; Guerrero-Castillo et al., 2017; Garcia et al., 2017; Stroud et al., 2016; Figure 1—figure supplement 1A). In plants, more similar to what occurs in bacterial CI assembly (Friedrich et al., 2016), the final assembly step is the attachment of the PD-module onto an intermediate (termed CI*) that already contains the N-, Q- and PP-modules (Ligas et al., 2019; Figure 1—figure supplement 1B). This difference in the order of assembly of CI in plants vs. metazoans is significant: in metazoans, adding the NADH dehydrogenase N-module last ensures that no assembly intermediate is capable of transferring electrons from NADH to CoQ. This is believed to have protective roles, to prevent the formation of reactive oxygen species during the CI assembly process (Parey et al., 2019). In contrast, the plant CI* intermediate contains all the subunits and co-factors needed to carry out NADH:CoQ oxidoreduction.
In contrast to the large number of recent high-resolution structures of mammalian and yeast respiratory complexes and supercomplexes, the most detailed plant CI structures known were obtained by negative-stain electron microscopy (EM) two-dimensional (2D) classifications from Solanum tuberosum (potato) and Arabidopsis thaliana or sub-tomogram averaged reconstructions that lack secondary structure details (Bultema et al., 2009; Dudkina et al., 2005; Davies et al., 2018). The paucity of functional and structural data for plant mETC complexes stems in large part from the limited availability of sufficient protein sample needed for structural analysis (Dudkina et al., 2015). Indeed, it has been difficult to obtain intact plant mitochondria in sufficient amounts for preparative biochemical fractionation. A typical reported yield of mitochondria is ~0.2–0.5 mg mitochondria/g fresh weight of starting plant material (Luster and Fites, 1987), which contrasts with a yield of ~30 mg mitochondria/g fresh weight from mammalian sources. In light of these challenges, most of the biochemical data on plant mETC have used intact mitochondria (e.g. oxygen-consumption experiments) or complexes that have been electro-eluted from electrophoretic gels (Bultema et al., 2009; Dudkina et al., 2005; Dudkina et al., 2006; Eubel et al., 2005). Although such electro-eluted protein samples have yielded the low-resolution structures described above and have proven suitable for proteomic studies, the low yields and low activities of these protein samples have so far thwarted detailed functional or structural analyses of the plant mETC complexes. A detailed understanding of the energy-converting mechanisms of plant respiratory mETC complexes and supercomplexes requires improved protocols for their extraction from plant mitochondrial membranes, and their purification in sufficient amounts while maintaining them in a functionally active state.
Here, we present a cryoEM structure of an ~800 kDa assembly intermediate of plant mitochondrial CI from etiolated Vigna radiata (mung bean) hypocotyls at 3.9 Å resolution. This assembly intermediate, CI* (Ligas et al., 2019), contains the intact peripheral arm (N- and Q-modules) as well as the PP-module and γCA domain, but lacks the PD-module. Our structure allowed us to build the first atomic model for any mitochondrial CI species from the plant kingdom and revealed important differences in the CI core and accessory subunits between plants, mammals, yeast and bacteria. Such subunit differences shed light on the known differences in CI assembly in plants versus opisthokonts. The structure also allowed us to define the interface between the γCA domain and the membrane arm of CI and revealed a key role for lipids in this interaction. We also discuss the implications of our findings on the possibility that CI* may provide additional flexibility to plants’ mETC.
Results
Structure of a mitochondrial Complex I assembly intermediate from etiolated V. radiata (mung bean)
In order to investigate the plant mitochondrial electron transport chain, we identified V. radiata (mung bean) as an optimal model system. V. radiata offers several advantages for plant mitochondrial research: i) it can be easily sprouted and harvested within six days, ii) it can be grown in the dark (etiolated) to minimize development of chloroplasts, which would otherwise contaminate the mitochondrial preparations, iii) its age and growth conditions can be controlled experimentally, iv) its genome has been sequenced and v) its mitochondrial content has been reported to be higher than other plant sources previously used for plant mitochondrial research (Luster and Fites, 1987). Moreover, we have optimized standard plant mitochondria isolation protocols (Millar et al., 2007) to routinely obtain ~1 g of wet weight mitochondria per 1 kg of etiolated V. radiata hypocotyls, approximately 3–4 times what has been previously reported (Luster and Fites, 1987).
Isolation of the mitochondrial electron transport complexes of V. radiata was performed by extraction from washed mitochondrial membranes using the gentle detergent digitonin, followed by exchange into the amphipathic polymer A8-35 to further stabilize the complexes. The presence of complex I (CI)-containing bands was analyzed using a standard in-gel NADH-dehydrogenase activity assay for CI on a blue-native gel (BN-PAGE) (Schertl and Braun, 2015). As expected from previously reported plant mitochondrial extractions (Bultema et al., 2009; Dudkina et al., 2005; Eubel et al., 2004a; Eubel et al., 2004b; Eubel et al., 2003; Krause et al., 2004), we observed a number of bands with NADH-dehydrogenase activity, representing CI in different assembly states, such as in mitochondrial supercomplexes (Bultema et al., 2009; Dudkina et al., 2005; Eubel et al., 2004a; Eubel et al., 2004b; Eubel et al., 2003; Krause et al., 2004; Dudkina et al., 2010; Figure 1—figure supplement 2A). The amphipol-stabilized complexes and supercomplexes were separated on a linear sucrose gradient (Figure 1—figure supplement 2B–C). Two peaks displaying NADH-dehydrogenase activity were of sufficient amount to be further purified by size-exclusion chromatography (Figure 1—figure supplement 2D). These purified fractions retained their NADH-dehydrogenase activity by in-gel activity assays (Figure 1—figure supplement 2E). Moreover, these fractions also showed NADH-decylubiquinone oxidoreductase activity using a standard CI spectroscopic activity assay (Huang et al., 2015; Figure 1—figure supplement 2F). These fractions were investigated by single-particle cryoEM. Here, we present results from the lower molecular weight fraction (‘peak 2’) (Figure 1—figure supplement 2G–H).
Structural analysis revealed that this fraction contained an ~800 kDa CI subcomplex, previously identified as a plant mitochondrial CI assembly intermediate termed complex I* (CI*, Figure 1B), which we were able to resolve to a nominal resolution of 3.9 Å (Figure 1C, Tables 1–2, Video 1). The existence of this assembly intermediate has been determined by genetic and mitochondrial proteomics experiments of CI's assembly pathway in etiolated seedlings (Heazlewood et al., 2003) and non-etiolated seedlings and leaves of Arabidopsis thaliana (Ligas et al., 2019; Meyer et al., 2011; Schertl et al., 2012; Schimmeyer et al., 2016; Senkler et al., 2017), as well in non-etiolated leaves of Nicotiana sylvestris (Pineau et al., 2008). Moreover, the A. thaliana and N. sylvestris CI* intermediate shows NADH-dehydrogenase activity by the same in-gel activity assay used in our preparation (Meyer et al., 2011; Pineau et al., 2008; Haïli et al., 2013). CI* contains CI’s intact peripheral arm (N- and Q-modules), PP-module and γCA domain. However, it is missing the two membrane arm core subunits NU4M and NU5M and their associated accessory subunits that form the PD-module (Figure 1B). As expected from complexome profiling analyses (Ligas et al., 2019; Senkler et al., 2017), our structure of CI* is composed of over 30 subunits of the N-module, Q-module, PP-module and the γCA domain. Throughout this manuscript, we use the plant nomenclature for the subunits (see Table 3 for subunit name conversions).
CryoEM density for the CI* composite map.
Cryo-EM data collection, refinement and validation statistics.
Data Collection and processing | |||
---|---|---|---|
Microscope | Titan krios, (UCSF) | ||
Camera | K3 | ||
Magnification | 60,010 | ||
Voltage (kV) | 300 kV | ||
Electron exposure (e-/Å2) | 86.4 | ||
Defocus range (µm) | −0.5 to −2.0 | ||
Pixel size (Å) | 0.8332 | ||
Software | SerialEM | ||
Reconstruction | CI* Peripheral Arm | CI* Membrane Arm | CI* Composite Map |
Number of particles | 34,407 | 34,407 | The CI* Peripheral Arm and Membrane Arm Maps were combined in Phenix to generate this composite map |
Accuracy of rotations (°) | 0.68 | 1.489 | |
Accuracy of translations (pixels) | 0.655 | 0.881 | |
Box size (pixels) | 512 | 512 | |
Final resolution (Å) | 3.8 | 3.9 | |
Map sharpening B factor (Å2) | −90 | −96 | |
EMDB ID | 22093 | 22092 | 22090 |
Refinement | |||
Software | Phenix | ||
Initial model (PDB code) | 6Q9D | 6Q9B and 1QRG | 6Q9D, 6Q9B and 1QRG |
Map/model correlation | |||
Model resolution (Å) | 3.9 | 4.0 | 3.9 |
d99 (Å) | 3.9 | 4.0 | 4.0 |
FSC model 0.5 (Å) | 3.9 | 3.9 | 3.9 |
Map CC (around atoms) | 0.82 | 0.86 | 0.87 |
Model composition | |||
Non-hydrogen atoms | 26,001 | 19,052 | 45047 |
Protein residues | 3284 | 2453 | 5736 |
Number of chains | 17 | 18 | 34 |
Number of ligands and cofactors | 11 | 1 | 12 |
Number of lipids | 0 | 6 | 6 |
Atomic Displacement Parameters (ADP) | |||
Protein average (Å2) | 68.78 | 58.40 | 64.39 |
Ligand average (Å2) | 48.59 | 48.59 | 48.59 |
R.m.s. deviations | |||
Bond lengths (Å) | 0.007 | 0.007 | 0.007 |
Bond angles (°) | 1.187 | 1.122 | 0.845 |
Ramachandran Plot | |||
Favored (%) | 82.90 | 88.03 | 84.98 |
Allowed (%) | 16.76 | 11.88 | 14.79 |
Disallowed (%) | 0.34 | 0.08 | 0.23 |
Validation | |||
MolProbity score | 2.41 | 2.31 | 2.38 |
Clash score | 16.79 | 16.21 | 16.42 |
Rotamer outliers (%) | 0.25 | 0.20 | 0.23 |
EMRinger score | 1.47 | 2.09 | 2.17 |
PDB ID | --- | --- | 6X89 |
Model building statistics by subunit.
Subunit name | Uniprot ID | Chain ID | Total residues | Atomic residues | Poly-Ala | Un-modeled residues | % atomic | TMH | Identified RNA editing sites* | Ligands, lipids |
---|---|---|---|---|---|---|---|---|---|---|
Peripheral arm core subunits | ||||||||||
NDUS1 | A0A1S3TQ85 | S1 | 746 | 57–744 | 57 | 1–56, 745–746 | 92.1% | 4Fe4S×2, 2Fe2S | ||
NDUS2 | E9KZN6 | S2 | 394 | 9–17,21-394 | 1–8, 18–20 | 98.0% | S26L, 246L, S67F, H82Y, S84L, R106C, S112L, S193L, S233L, H242Y, S245L, P247F, R257C, R353C, S360F, S363L, S368F, P375L | |||
NDUS3 | E9KZM7 | S3 | 190 | 1–184 | 185–190 | 96.8% | S31F, S56L, P100S, R110W, S133L, L147F | |||
NDS7 | A0A1S3U8J5 | S7 | 213 | 56–213 | 1–55 | 74.2% | 4Fe4S, PC | |||
NDS8 | A0A1S3VGS8 | S8 | 222 | 42–222 | 1–41 | 81.5% | 4Fe4S×2 | |||
NDUV1 | A0A1S3V7V2 | V1 | 491 | 59–491 | 1–58 | 88.2% | 4Fe4S, FMN | |||
NDUV2 | A0A1S3U769 | V2 | 251 | 28–243 | 1–27, 244–251 | 86.1% | 2Fe2S | |||
Peripheral arm accessory subunits | ||||||||||
NDUA2 | A0A1S3TVC7 | A2 | 98 | 2–93 | 1, 94–98 | 93.9% | ||||
NDUA5 | A0A1S3U023 | A5 | 169 | 12–137 | 1–11, 138–169 | 74.6% | ||||
NDUA6 | A0A1S3W1K8 | A6 | 132 | 118–131 | 1–117, 132 | 11.4% | ||||
NDUA7 | A0A1S3UVC7 | A7 | 127 | 19–127 | 1–18 | 85.8% | ||||
NDUA9 | A0A1S3V8W7 | A9 | 396 | 47–381 | 1–46, 382–396 | 84.6% | NADPH | |||
NDUA12 | A0A1S3VNK7 | AL | 156 | 21–155 | 1–20, 156 | 86.5% | ||||
NDUS4 | A0A1S3UIW7 | S4 | 146 | 42–142 | 142 | 1–41, 143–146 | 69.2% | |||
NDUS6 | A0A1S3VYF3 | S6 | 103 | 31–102 | 1–30, 103 | 69.9% | Zn2+ | |||
Membrane arm core subunits | ||||||||||
NU1M | A0A1S4ETV6/E9KZL0 | 1M | 325 | 2–213, 220–325 | 1, 214–219 | 97.8% | 8 | R89W, P164S, R165C, S167L, S179F, R225C, P242L, P248L, P252L, R300W, R310W | ||
NU2M | E9KZK9 | 2M | 488 | 1–487 | 488 | 99.8% | 14 | S19F, S103F, S104F, P119L, P121S, R123C, H132Y, P143L, S166LL, S221F, P307L, H310Y, R320C, S376L, S467L, S468F, S486L | PC×2 | |
NU3M | Q9XPB4 | 3M | 118 | 1–28, 56–118 | 29–55 | 77.1% | 3 | P70F, P83S, P84L, S115L, R117W | ||
NU4LM | A0A1S4ETY3/E9KZN8 | 4L | 100 | 1–86 | 87–100 | 86.0% | 3 | S14F, P29L, S32L, P34S, S37L, S53L, S63L, S66L | ||
NU6M | E9KZM5 | 6M | 205 | 1–72, 111–172 | 73–110 | 173–205 | 65.4% | 5 | P9L, A18V, P30F, P32L, R35C, P54L, H57Y | |
Membrane arm accessory subunits | ||||||||||
CA1 | A0A1S3VT00 | G1 | 270 | 3–222 | 223–233 | 1–2, 234–270 | 81.5% | |||
CA2 | A0A1S3U544 | G2 | 273 | 2–237 | 1, 238–273 | 86.4% | ||||
CAL2 | A0A1S3UI49 | L2 | 256 | 49–129, 134–254 | 1–48, 130–133, 255–256 | 80.5% | ||||
NDUX1 | A0A1S3VI15 | X1 | 101 | 1–100 | 101 | 99.0% | 2 | |||
NDUC2 | A0A1S3UPL8 | C2 | 81 | 5–68 | 1–4, 69–81 | 79.0% | 2 | |||
NDUA8 | A0A1S3VVN6 | A8 | 106 | 2–106 | 1 | 99.1% | ||||
NDUA13 | A0A1S3UYW0 | AM | 143 | 2–143 | 1 | 99.3% | 1 | |||
NDUA1 | A0A1S3TU57 | A1 | 65 | 2–63 | 1, 64–65 | 95.4% | 1 | PC | ||
NDS5 | A0A1S3TQ33 | S5 | 399 | 2–70 | 1, 71–399 | 17.3% | ||||
NDUA3 | A0A1S3TCK0 | A3 | 63 | 2–45 | 1, 46–63 | 69.8% | 1 | |||
P2 | A0A1S3TGE7 | P2 | 115 | 83–106 | 77–82 | 1–76, 107–115 | 20.9% | |||
Unassigned density | ||||||||||
A | 1–18 | |||||||||
B | 1–24 | |||||||||
C | 1–43 | 1 |
-
*RNA editing of mitochondrially encoded subunits: amino acids were changed at the listed positions as detailed. The changes were based on the reported equivalent A. thaliana RNA edits (Giegé and Brennicke, 1999; Bentolila et al., 2008) and were only made when density was unambiguously correct for the edited V. radiata amino acid in the cryoEM map.
Complex I subunit homologues in plants, mammals, yeast and bacteria.
V. radiata homologues were obtained by performing BLASTp searches of the Arabidopsis thaliana genes (Meyer et al., 2019; Braun et al., 2014). Mammalian, yeast and bacterial homologues were obtained from Letts and Sazanov, 2015. Additional BLASTp searches were performed wherever necessary. Given the high sequence similarity between the carbonic anhydrase (CA) paralogues, the names of the V. radiata CA proteins appear to have been mis-assigned in the genetic databases relative to their A. thaliana homologues. The CA1, CA2, CA2-like nomenclature used in the table is the one that, based on our sequence alignments, best represents homology to the A. thaliana CA proteins. N, NADH-binding module; Q, quinone-binding module; PP, proximal-pumps module; PD, distal-pumps module; CA, carbonic anhydrase domain.
Module | Vigna radiata protein name | Vigna radiata gene | Vigna radiata uniprot identifier | Arabidopsis thaliana protein name | Arabidopsis thaliana gene | Homo sapiens name | Ovis aries name | Mus musculus name | Yarrowia lipolytica name | Thermus thermophilus name |
---|---|---|---|---|---|---|---|---|---|---|
CORE peripheral arm | ||||||||||
N | NDUS1 | LOC106757688 | A0A1S3TQ85 | 75 kDa | At5g37510 | NDUFS1 | NDUFS1 | NDUFS1 | NUAM | Nqo3 |
N | NDUV1 | LOC106772405 | A0A1S3V7V2 | 51 kDa | At5g08530 | NDUFV1 | NDUFV1 | NDUFV1 | NUBM | Nqo1 |
N | NDUV2 | LOC106762461 | A0A1S3U769 | 24 kDa | At4g02580 | NDUFV2 | NDUFV2 | NDUFV2 | NUHM | Nqo2 |
Q | NDUS2 | nad7 | E9KZN6 | Nad7 | AtMg00510 | NDUFS2 | NDUFS2 | NDUFS2 | NUCM | Nqo4 |
Q | NDUS3 | nad9 | E9KZM7 | Nad9 | AtMg00070 | NDUFS3 | NDUFS3 | NDUFS3 | NUGM | Nqo5 |
Q | NDS7 | LOC106762764 | A0A1S3U8J5 | PSST | At5g11770 | NDUFS7 | NDUFS7 | NDUFS7 | NUKM | Nqo6 |
Q | NDS8 | LOC106775047 | A0A1S3VGS8 | TYKY | At1g79010, At1g16700 | NDUFS8 | NDUFS8 | NDUFS8 | NUIM | Nqo9 |
CORE membrane arm | ||||||||||
PP | NU1M | nad1 | A0A1S4ETV6 | Nad1 | AtMg00516, AtMg01120, AtMg01275 | MT-ND1 | MT-ND1 | MT-ND1 | NU1M | Nqo8 |
PP | NU2M | nad2 | E9KZK9 | Nad2 | AtMg00285, AtMg01320 | MT-ND2 | MT-ND2 | MT-ND2 | NU2M | Nqo14 |
PP | NU3M | nad3 | Q9XPB4 | Nad3 | AtMg00990 | MT-ND3 | MT-ND3 | MT-ND3 | NU3M | Nqo7 |
PP | NU4LM | nad4L | A0A1S4ETY3 | Nad4L | AtMg00650 | MT-ND4L | MT-ND4L | MT-ND4L | NULM | Nqo11 |
PP | NU6M | nad6 | E9KZM5 | Nad6 | AtMg00270 | MT-ND6 | MT-ND6 | MT-ND6 | NU6M | Nqo10 |
PD | NU4M* | nad4 | E9KZL8 | Nad4 | AtMg00580 | MT-ND4 | MT-ND4 | MT-ND4 | NU4M | Nqo13 |
PD | NU5M* | nad5 | E9KZL1 | Nad5 | AtMg00060, AtMg00513, AtMg00665 | MT-ND5 | MT-ND5 | MT-ND5 | NU5M | Nqo12 |
ACCESSORY membrane arm | ||||||||||
PP | NDUA1 | LOC106758834 | A0A1S3TU57 | MWFE | At3g08610 | NDUFA1 | NDUFA1 | NDUFA1 | NIMM | - |
PP | NDUA3 | LOC106754061 | A0A1S3TCK0 | B9 | At2g46540 | NDUFA3 | NDUFA3 | NDUFA3 | NI9M | - |
PP | NDUA8-B | LOC106778955 | A0A1S3VVN6 | PGIV | At3g06310, At5g18800 | NDUFA8 | NDUFA8 | NDUFA8 | NUPM | - |
PP | NDUA13-A | LOC106769964 | A0A1S3UYW0 | B16.6 | At2g33220, At1g04630 | NDUFA13 | NDUFA13 | NDUFA13 | NB6M | - |
PP | NDUC1 | LOC106771273 | A0A1S3V2Z3 | - | - | NDUFC1 | NDUFC1 | NDUFC1 | - | - |
PP | NDUS5 | LOC106757655 | A0A1S3TQ33 | 15 kDa | At3g62790, At2g47690 | NDUFS5 | NDUFS5 | NDUFS5 | NIPM | - |
PP | NDUB8 | LOC106765859 | A0A1S3UJ95 | ASHI | At5g47570 | NDUFB8 | NDUFB8 | NDUFB8 | NIAM | - |
PP | NDUB10-B | LOC106774903 | A0A1S3VGT1 | PDSW | At1g49140, At3g18410 | NDUFB10 | NDUFB10 | NDUFB10 | NIDM | - |
PP | NDUA11* | LOC106756741 | A0A1S3TLY8 | B14.7 | At2g42210 | NDUFA11 | NDUFA11 | NDUFA11 | NUJM | - |
Module | Vigna radiata protein name | Vigna radiata gene | Vigna radiata Uniprot identifier | Arabidopsis thaliana protein name | Arabidopsis thaliana gene | Homo sapiens name | Ovis aries name | Mus musculus name | Yarrowia lipolytica name | Thermus thermophilus name |
ACCESSORY membrane arm | ||||||||||
PD | NDUB1* | LOC106775330 | A0A1S3VI15 | MNLL | At4g16450 | NDUFB1 | NDUFB1 | NDUFB1 | - | - |
PD | NDUC2* | LOC106767534 | A0A1S3UPL8 | B14.5b | At4g20150 | NDUFC2 | NDUFC2 | NDUFC2 | NEBM | - |
PD | NDUB2* | LOC106754955 | A0A1S3TFG6 | AGGG | At1g76200 | NDUFB2 | NDUFB2 | NDUFB2 | - | - |
PD | NDUB3* | LOC106769121 | A0A1S3UVV0 | B12 | At2g02510, At1g14450 | NDUFB3 | NDUFB3 | NDUFB3 | NB2M | - |
PD | NDUB4* | LOC106766640 | A0A1S3ULL3 | B15 | At2g31490 | NDUFB4 | NDUFB4 | NDUFB4 | NB5M | - |
PD | NDUB5* | LOC106767179 | A0A1S3UND4 | - | - | NDUFB5 | NDUFB5 | NDUFB5 | NUNM | - |
PD | NDUB7* | LOC106770979 | A0A1S3V2B8 | B18 | At2g02050 | NDUFB7 | NDUFB7 | NDUFB7 | NB8M | - |
PD | NDUB9* | LOC106760947 | A0A1S3U1J6 | B22 | At4g34700 | NDUFB9 | NDUFB9 | NDUFB9 | NI2M | - |
PD | NDUB11* | LOC106771273 | A0A1S3V2Z3 | ESSS | At2g42310, At3g57785 | NDUFB11 | NDUFB11 | NDUFB11 | NESM | - |
ACCESSORY peripheral arm | ||||||||||
N | NDUA2 | LOC106759195 | A0A1S3TVC7 | B8 | At5g47890 | NDUFA2 | NDUFA2 | NDUFA2 | NI8M | - |
N | NDUA12 | LOC106776991 | A0A1S3VNK7 | B17.2 | At3g03100 | NDUFA12 | NDUFA12 | NDUFA12 | N7BM | - |
N | NDUS4 | LOC106765762 | A0A1S3UIW7 | 18 kDa | At5g67590 | NDUFS4 | NDUFS4 | NDUFS4 | NUYM | - |
N | NDUS6 | LOC106779709 | A0A1S3VYF3 | 13 kDa | At3g03070 | NDUFS6 | NDUFS6 | NDUFS6 | NUMM | - |
Q | NDUA5 | LOC106760411 | A0A1S3U023 | B13 | At5g52840 | NDUFA5 | NDUFA5 | NDUFA5 | NUFM | - |
Q | NDUA6† | LOC106780789 | A0A1S3W1K8 | B14 | At3g12260 | NDUFA6 | NDUFA6 | NDUFA6 | NB4N | - |
Q | NDUA7 | LOC106768957 | A0A1S3UVC7 | B14.5a | At5g08060 | NDUFA7 | NDUFA7 | NDUFA7 | NUZM | - |
Q | NDUA9 | LOC106772694 | A0A1S3V8W7 | 39 kDa | At2g20360 | NDUFA9 | NDUFA9 | NDUFA9 | NUEM | - |
Plant-specific accessory | ||||||||||
CA | CA1‡ | LOC106778103 | A0A1S3VT00 | Gamma-CA 1 | At1g19580 | - | - | - | - | - |
CA | CA2§ | LOC106761992, LOC106761993 | A0A1S3U566, A0A1S3U544 | Gamma-CA 2 | At1g47260 | - | - | - | - | - |
CA | CA2-L¶ | LOC106765552 | A0A1S3UI49 | Gamma CA-like 2 | At3g48680 | - | - | - | - | - |
CA | CA3* | n.a.** | n.a.** | Gamma-CA 3 | At5g66510 | - | - | - | - | - |
CA | CA1-L* | n.a.** | n.a.** | Gamma-CA-like 1 | At5g63510 | - | - | - | - | - |
PP | NDUX1†† | LOC106775330 | A0A1S3VI15 | 20.9 kDa | At4g16450 | - | - | - | NUXM | - |
PP | P2/16 kDA | LOC106755236 | A0A1S3TGE7 | P2 | At2g27730 | - | - | - | - | - |
Plant-specific accessory | ||||||||||
Unconfirmed plant CI subunits (not seen in CI*) | MICOS (DUF543) | LOC106779628 | A0A1S3VY06 | MICOS subunit Mic10 | At1g72165 | - | - | - | - | - |
Uncharacterized protein LOC106758628 | LOC106758628 | A0A1S3TTD7 | NDU10 | At4g00585 | - | - | - | - | ||
P1/11 kDA | LOC106761134 | A0A1S3U2B9 | P1 | At1g67350 | - | - | - | - | - | |
P3 | LOC106755586 | A0A1S3THM0 | P3 | At5g14105 | - | - | - | - | - | |
P4 | LOC106767179 | A0A1S3UND4 | P4 | At1g67785 | - | - | - | - | - | |
TIM22−4 × 1 | LOC106779665 | A0A3Q0EN44 | TIM22-4 | At1g18320 | - | - | - | - | - | |
TIM22−4 × 2 | LOC106779665 | A0A1S3VZ08 | TIM22-1 | At3g10110 | - | - | - | - | - | |
TIM23-2 | LOC106761237 | A0A1S3U2K1 | TIM23-2 | At1g72750 | - | - | - | - | - | |
Uncharacterized protein LOC106768488 isoform X4 | LOC106768488 | A0A1S3UST2 | SH3/FCH domain protein | At1g68680 | - | - | - | - | - | |
UDP-galactose transporter 1 | LOC106762681 | A0A1S3U838 | TPT domain-containing protein | At1g72180 | - | - | - | - | - | |
Gravitropic in the light 1 | LOC106779790 | A0A1S3VYR1 | DUF641 domain-containing protein | At2g28430 | - | - | - | - | - |
-
*Not seen in CI*.
†Only the C-terminus seen in CI* (see Main body and Discussion).
-
‡Called gamma carbonic anhydrase one in Uniprot.
§Called gamma carbonic anhydrase 1, mitochondrial in Uniprot (mis-assigned in the database).
-
¶Called gamma carbonic anhydrase-like 2, mitochondrial in Uniprot (mis-assigned in the database).
**Homologue not found using BLASTp.
-
††New identified subunit.
Key differences in observed core subunits
The peripheral and membrane arm core subunits present in the structure of CI* are structurally homologous to the bacterial, yeast and mammalian CI core subunits, with a few notable differences.
The N-terminus of core Q-module subunit NDUS2 is shortened in V. radiata compared to NDUS2 from Y. lipolytica and mammals, in which the N-terminus of NDUS2 extends from the interface of the peripheral and membrane arms of the complex along the matrix side of the membrane arm. Whereas in Y. lipolytica the N-terminus of NDUS2 binds to the matrix surface of core H+-pumping subunit NU2M, in mammals the N-terminus of NDUS2 extends further along the membrane arm and binds to the matrix surface of core H+-pumping subunit NU4M, bridging across the PP- and PD-modules. In contrast, V. radiata NDUS2 is ~40 amino acid residues shorter on the N-terminus compared to mammals and does not extend along the membrane arm. Moreover, the equivalent path for the Y. lipolytica or mammalian NDUS2 N-terminus in V. radiata is blocked by the γCA domain to the plant PP-module on the membrane arm.
The N-terminus of core peripheral arm subunit NDUS8 is also divergent between plants, fungi and mammals. In V. radiata, the N-terminus possesses an additional α-helix that binds between the Q-module accessory subunit NDUA5 and the PP-module core membrane subunit NU2M, enlarging the interaction interface between the peripheral and membrane arms (Figure 1C). In Y. lipolytica, the N-terminus of NDUS8 forms an extended coil that reaches up along the peripheral arm between the Q-module accessory subunits NDUA5 and NDUA7, making contact with the core Q-module subunit NDUS3. In contrast, the N-terminus of mammalian NDUS8 folds back along the surface of the membrane arm and tucks underneath the Q-module accessory subunit NDUA7. In Y. lipolytica, this binding site, underneath the NDUA7 homologue (NUZM), is occupied by NUZM’s C-terminus, which folds back under itself. However, in V. radiata the binding site underneath NDUA7 is occupied by an unidentified subunit that extends from this pocket under NDUA7 toward the core transmembrane subunits adjacent to the NU3M transmembrane helix (TMH) 1–2 loop and the NU6M TMH3, which undergo conformational changes during CI’s enzymatic turnover in the fungal structures (Agip et al., 2018; Letts et al., 2019; Parey et al., 2018). Although the identity of this sequence in the V. radiata structure remains unclear, it appears to be unique to plant CI.
Core subunit NU2M in V. radiata CI* contains three N-terminal transmembrane helices that are present in yeast and bacterial complexes, but lost in the metazoan lineage (Birrell and Hirst, 2010). Moreover, V. radiata CI* contains a homologue of Y. lipolytica’s accessory subunit NUXM (absent in metazoans), which binds to the NU2M N-terminal transmembrane helices. Based on the Y. lipolytica subunit name, we coined this subunit of V. radiata CI NDUX1. The presence of this subunit in both plants and fungi suggests that this subunit was present in the ancestral eukaryotic CI before the unikont/bikont lineage divergence but was lost in metazoans when NU2M became N-terminally truncated. The first transmembrane helix of NU2M in Y. lipolytica is notably short (only 15 amino acids), enters only to the midplane of the membrane and is bound by a membrane-penetrating loop of the accessory subunit NUXM. In contrast, in bacteria (T. thermophilus and E. coli) and V. radiata, the first transmembrane helix of NU2M spans the full length of the membrane. Furthermore, the loop connecting V. radiata’s NU2M TMH1-2 in the mitochondrial matrix is longer than in any of the other CI structures and extends into the matrix, where it contacts the N-terminal helix of NDUS8 discussed above. Given the universality of the hinging motion between CI’s peripheral and membrane arms, seen in the structures of several organisms (Agip et al., 2018; Letts et al., 2019; Parey et al., 2018), the additional interaction surface formed by NDUS8 and NU2M in V. radiata CI is likely functionally relevant.
Key differences in observed accessory subunits
Although the majority of the accessory subunits present in CI* have homologues in fungi and mammals (opisthokonts), there are a number of notable differences.
In the plant complex, the peripheral arm accessory subunit NDUS6 lacks an N-terminal domain that is seen in both the Y. lipolytica and mammalian structures (Figure 2A). In Y. lipolytica, mammals and V. radiata, the C-terminal, Zn2+-containing domain of NDUS6 binds mainly to the core subunits NDUS1, NDUS8 and NDUS2 at the interface of the N- and Q-modules. However, in opisthokonts, the N-terminal domain of NDUS6 binds to the Q-module at an additional site through contacts with the membrane-anchored NDUA9 accessory subunit (Figure 2A). In order to bind across these two locations, NDUS6 in opisthokonts extends above the C-terminus of accessory subunit NDUA12. This arrangement determines the order of assembly of these subunits in opisthokonts, as NDUA12 must be bound to the peripheral arm before the N-terminal domain of NDUS6 binds. However, due to the lack of the N-terminal domain in V. radiata’s NDUS6, there is no interaction with NDUA9 nor traversing of the NDUA12 C-terminus. This difference has important implications for the assembly of CI in plants versus opisthokonts. In opisthokonts, the interaction between NDUS6, NDUA12 and the NDUA12-homologous assembly factor NDUFAF2 establishes an important checkpoint for assembly of the peripheral arm. Thus, the lack of the NDUS6 N-terminus may in part explain observed differences between the assembly pathways of plant and opisthokont CI (see Discussion).

Key differences in CI accessory subunits between V. radiata and opisthokonts.
Accessory subunits NDUS6 and NDUA12, NDUA8 and NDUC2 of V. radiata (this study), Y. lipolytica (PDB:6RFR) and O. aries (PDB: 6QA9) are shown as surface for comparison. (A) NDUS6 (green) and NDUA12 (orange), with an additional label for NDUA9. (B) NDUA8 (maroon), with additional label for the V. radiata’s carbonic anhydrase domain (CA). (C) NDUC2 (blue), with additional label for the V. radiata’s CA.
Other key differences can be seen on the intermembrane space side of the membrane arm in accessory subunits NDUA8 and NDUC2. Compared to both Y. lipolytica and mammals, the double-CHCH domain of the PP-module NDUA8 subunit, which binds to the ‘heel’ of the complex on the intermembrane space (Figure 2B), is C-terminally truncated in V. radiata. In the Y. lipolytica structure, the C-terminus of NDUA8 folds back onto itself with an additional α-helix, forming a bulkier subunit and a further interaction interface with the core transmembrane subunit NU1M. More interestingly, in mammals, the C-terminus of NDUA8 extends as a long coil halfway along the membrane arm and binds in a pocket between NU2M and NU4M at the interface of the PP-module and PD-module. The PP-module accessory subunit NDUC2 is also C-terminally truncated in V. radiata and Y. lipolytica relative to NDUC2 in mammals (Figure 2C). In all mitochondrial CI structures to date, this subunit binds to the final transmembrane helix of the core NU2M subunit. However, in mammals, the NDUC2 C-terminus forms an extended coil on the intermembrane space side of the complex that extends along the membrane arm to interact with NDUB10 and NDUB11, bridging the PP- and PD- modules. This bridging interaction is also present in Y. lipolytica via an extended loop on the PD-module core subunit NU4M.
This pattern of truncated core and accessory subunits or missing interactions (e.g. NDUS2, NDUA8 and NDUC2; Table 4) in V. radiata relative to those in opisthokonts likely diminishes the stability of the attachment of PP-module to the PD-module, which may have consequences for CI’s function and assembly (see Discussion).
PP- and PD-module bridging subunits in mammalian, Y. lipolytica and V. radiata CI.
Subunits discussed in the manuscript are marked with two asterisks (**). Bridging interactions are shaded in green. Lack of interactions by existing subunits or lack of homologues are shaded in orange. Lack of the PD subunits in V. radiata CI* is shaded in yellow. PP, proximal pumping domain; PD, distal pumping domain.
Location | Subunit | Mammals | Y. lipolytica | V. radiata |
---|---|---|---|---|
Inter-membrane space (IMS) | NDUA8** | Extends along membrane arm, bridges NU2M (PP) and NU4M (PD) | Does not extend to the PP/PD-module interface but has an additional helix interacts with NU1M | C-terminally truncated (does not bridge) |
NDUC2** | C-terminus bridges NDUB10 (PP) and NDUB11 (PD) | C-terminally truncated, but bridging interaction replaced by extended loop on NU4M | C-terminally truncated (does not bridge) | |
NDUB5 | Bridging interactions | Bridging interactions | N- and C-terminally truncated (subunit not present in CI*) | |
NDUA11 | Does not bridge in the IMS | C-terminal extension binding to NU4M | Subunit not present in CI* | |
Membrane | NDUA11 | Binds to the lateral helix of NU5M, connecting NU5M and NU2M | Binds to the lateral helix of NU5M, connecting NU5M and NU2M | Subunit not present in CI* |
Matrix | NDUS2** | Bridging interactions | Does not bridge | N-terminally truncated (does not bridge) |
NU5M | Lateral helix extends into PP | Lateral helix extends into PP | Subunit not present in CI* | |
NDUA10 | Bridging interactions | No homologue present | No homologue present | |
NDUB11 | Bridging interactions | Does not bridge | Subunit not present in CI* | |
NDUB4 | Does not bridge | N-terminus extends along matrix arm and binds to NU2M | N-terminally truncated (subunit not present in CI*) |
Known Q-module accessory subunits not present in CI*
Compared to the mammalian and Y. lipolytica structures, two accessory subunits are absent from the Q-module in the CI* structure, namely the LYR-protein subunit NDUA6 and its accompanying acyl-carrier protein (ACPM1). The absence of the NDUA6 and ACPM1 subunits in CI* is notable given that, when the Y. lipolytica NDUA6 homologue is knocked out or mutated, this severely impacts the activity of the complex (Angerer et al., 2014). Therefore, although it is not completely understood how NDUA6 modulates the activity of CI, the lack of NDUA6 in CI* may be a way to regulate the activity of the assembly intermediate.
Although densities for NDUA6 and ACPM1 are absent in our CI* structure, density can be seen for a short α-helix bound under NDUS1, where the C-terminus of NDUA6 binds in both the Y. lipolytica and mammalian structures. This suggests that NDUA6 may be bound to CI* via its C-terminus, without fully engaging with the complex. Although this would be surprising, the density for the amino acid sidechains in this region is consistent with the sequence of the NDUA6 C-terminus; thus, this density was modelled as such. If correct, this suggests that NDUA6 may be attached to the Q module but unable to fully bind to its main site on NDUS2.
Plant-specific accessory subunits
V. radiata CI* does not have any plant-specific accessory subunits on the peripheral arm. Notwithstanding the unique features of NDUS6 and the absence of NDUA6 and ACPM1 discussed above, all of the V. radiata CI* N- and Q-module subunits have homologues in fungi and metazoans. However, this is not the case for the PP-module. Most notably, a large (~90 kDa) hetero-trimeric γCA domain lies on top of the core membrane arm subunit NU2M (Figure 1C).
The identity of the components of the plant γCA has remained elusive, with different three-way combinations of the five plant γCA proteins proposed based on different genetic and biochemical studies (Sunderhaus et al., 2006; Perales et al., 2004; Fromm et al., 2016; Cï Rdoba et al., 2019). Our structure allowed us to unambiguously assign the identity of the subunits of the γCA domain despite high sequence identity between the five carbonic anhydrase proteins in plants. Based on unambiguous density for key non-conserved residues, we were able to definitively assign the three different subunits of V. radiata CI* as CA1, CA2 and CA2L (Figure 3A).

V. radiata γ-carbonic-anhydrase (γCA) domain, Zn2+ coordination and associated lipid cavity.
(A) Top view of the carbonic anhydrase domain with its CA1 (green), CA2 (lime) and CAL2 (lime green). Key residues at subunit interfaces for Zn2+ coordination shown as sticks; Zn2+ shown as grey sphere. Only the CA1-CA2 interface has all three key Zn2+ coordinating histidines in place. (B) Zoom-in of Zn2+ coordination site in (A), with map density for the three histidines and Zn2+ shown as blue meshes. (C) Two phosphatidylcholines (spheres) are placed in the lipid cavity between the γCA and the PP-module. Asterisk indicates the N-terminal amphipathic helices of CA1 and CA2. (D Zoom-in of the lipid cavity in C), with lipid density shown as blue mesh and key interacting residues shown as sticks.
The interaction surface between the γCA domain and the PP-module (subunits NU2M, NDUC2, P2 and NDUX1) is large, covering an approximate surface of 3,740 Å2. As expected (Sunderhaus et al., 2006), the γCA interacts with the PP-module tightly, with an approximate gain of solvation free energy of −210 kcal/mol, which is almost twice as large as the solvation energy gain of association of the γCA hetero-trimer itself (Figure 3A, Table 5).
Quantification of interfaces within the γ-carbonic-anhydrase (γCA) domain and between γCA and the proximal pumping domain (PP) of CI*.
Interface residues, surface areas, solvation free energies and P-values were determined by uploading the molecular model of CI* into the the PDBePISA tool for the exploration of macromolecular interfaces (Krissinel and Henrick, 2007). The table with the full list of interaction surfaces for CI* was filtered for the interfaces involving CA1, CA2 or CAL2. Total values were obtained by adding the relevant two-way interactions, as per PDBePISA guidelines.
Subunit 1 | Subunit 2 | Inter-subunit interface | ||||||
---|---|---|---|---|---|---|---|---|
Subunit | # Interfacing residues | Interfacing surface area (Å2) | Subunit | # Interfacing residues | Interfacing surface area (Å2) | Interface surface area (Å2) | Interface solvation free energy (kcal/mol) | Solvation free energy gain P-value |
Within γCA domain | ||||||||
CA1 | 88 | 20,989 | CA2 | 80 | 22,088 | 2581.4 | −45.6 | 0.000 |
CA1 | 79 | 20,989 | CAL2 | 72 | 20,073 | 2637.8 | −38.2 | 0.001 |
CA2 | 73 | 22,088 | CAL2 | 66 | 20,073 | 2148.1 | −33.2 | 0.000 |
Total | 7367.3 | −117 | ||||||
Between γCA domain and membrane arm (PP) | ||||||||
NU2M | 16 | 7610 | CA1 | 16 | 20,989 | 259.2 | −17.1 | 0.001 |
NU2M | 16 | 7610 | CA2 | 16 | 22,088 | 261.1 | −17.1 | 0.001 |
NU2M | 16 | 7610 | CAL2 | 16 | 20,073 | 263.6 | −18.0 | 0.003 |
NDUC2 | 21 | 13,708 | CA1 | 21 | 20,989 | 427.8 | −20.2 | 0.000 |
NDUC2 | 16 | 13,708 | CA2 | 16 | 22,088 | 242.3 | −15.8 | 0.000 |
NDUC2 | 21 | 13,708 | CAL2 | 22 | 20,073 | 416.2 | −18.7 | 0.000 |
P2 | 21 | 10,433 | CA1 | 21 | 20,989 | 362.6 | −16.5 | 0.000 |
P2 | 26 | 10,433 | CA2 | 22 | 22,088 | 512.7 | −21.5 | 0.001 |
P2 | 17 | 10,433 | CAL2 | 17 | 20,073 | 244.6 | −16.8 | 0.001 |
NDUX1 | 16 | 13,955 | CA1 | 16 | 20,989 | 236.5 | −15.9 | 0.000 |
NDUX1 | 17 | 13,955 | CA2 | 17 | 22,088 | 268.1 | −15.0 | 0.000 |
NDUX1 | 16 | 13,955 | CAL2 | 16 | 20,073 | 240.9 | −16.7 | 0.000 |
Total | 3735.6 | −209.3 |
As has been previously demonstrated by proteomic analysis, the N-terminal mitochondrial signal pre-sequences for CA1 and CA2 remain uncleaved (Klodmann et al., 2010). We show here that these two N-terminal sequences together form a short α-helical coiled-coil-like structure (Figure 3C). This coiled coil is amphipathic and binds on the matrix surface of the inner mitochondrial membrane, contacting the NDUC2 and P2 subunits (see below) adjacent to the NU2M core subunit. In contrast, no density was observed for the N-terminal pre-sequence of CA2L, consistent with it being post-translationally cleaved (Huang et al., 2009).
The physiological role of the γCA domain on plant CI is unknown. Although recombinant mitochondrial γCA from plants has been shown to bind bicarbonate (HCO3-), it remains unclear whether it exhibits enzymatic activity (Martin et al., 2009). The canonical γCA trimer possesses three active sites, one at each interface between two protomers. Each active site is formed by three essential Zn2+-coordinating histidine residues. At each active site, two histidine residues are provided by one subunit and the third is provided by the adjacent subunit. However, in the plant CI γCA heterotrimer, the CA2L subunit is lacking two of the three essential histidine residues (Ala-147 and Arg-152 in V. radiata) that would be necessary to form active sites at the interfaces with the CA1 and CA2 subunits. This renders two of the possible three catalytic sites non-functional (Figure 3A, Figure 3—figure supplement 1). Furthermore, the V. radiata CA1 subunit is also missing one of the three Zn2+-coordinating histidine residues (Gln-135). Therefore, only one potentially catalytically active interface with all three Zn2+ coordinating residues remains in V. radiata’s γCA—namely, the site between CA1 and CA2 at the "top" (most matrix-exposed periphery) of the domain. Clear density for a Zn2+ can only be seen at this site (Figure 3B). In contrast, no Zn2+ is seen at either of the two other sites, whose mutated residues are chemically incompatible with ion coordination. It is also important to note that the plant CA1, CA2 and CAL2 proteins belong to the CamH subclass of γCAs, which lack the acidic loop containing the catalytically important ‘proton shuttle’ glutamate residue (Glu89 in the canonical γCA from Methanosarcina thermophila) (Zimmerman et al., 2010). While some members of the CamH subclass are catalytically active, some are not (Soto et al., 2006; Jeyakanthan et al., 2008). Therefore, carbonic anhydrase activity of the γCA domain of CI must be confirmed experimentally (Ferry, 2010).
The other plant-specific subunit we were able to assign in CI* was the single-transmembrane subunit P2. This subunit binds on top of NDUX1, adjacent to NU2M and directly underneath the γCA domain. The N-terminus of P2 interacts directly with the γCA domain in the matrix. Together, P2, NDUX1, NU2M and NDUC2 form a lipid-filled cavity positioned directly below the γCA domain (Figure 3C and D). Several positively charged residues from the γCA domain subunits can be seen interacting with these lipids, demonstrating that this lipid pocket also forms an important part of the γCA domain/membrane arm interface.
Unassigned density
We were unable to assign four small regions of density in the CI* structure. One is the region near the N-terminus of NDUS8 discussed above (Figure 4A). Another is the likely C-terminal helix of NDUA6 also discussed above (Figure 4B). The third is on the intermembrane space side of the membrane arm (Figure 4C). In both Y. lipolytica and mammalian CI, this binding site is occupied by the C-terminus of the PP- and PD-module-spanning subunit NDUB5. In Y. lipolytica and mammals, NDUB5 spans nearly the entire length of the membrane arm. In V. radiata CI*, the density for this subunit follows the equivalent path of NDUB5 in Y. lipolytica and mammals but becomes disordered by the PP-module's core subunit NU2M, which is adjacent to the C-terminus of accessory subunit NDUC2. The final stretch of unassigned density is for a single-transmembrane accessory subunit bound above NU6M TMH1 that contacts NU6M and NDUS5 on the intermembrane space side of the membrane arm (Figure 4D). This unassigned subunit protrudes away from CI* toward the location where CIII2 binds in the mammalian supercomplex I+III2 (Letts and Sazanov, 2015), suggesting a possible role for this subunit in supercomplex formation. No equivalent subunit is seen in either Y. lipolytica or mammalian CI, suggesting that this is a plant-specific subunit. However, due to local disorder, the density was too poor to assign the sequence from the reconstruction alone.
Catalytic sites
All the cofactors necessary for the transfer of electrons between NADH and CoQ are present in the CI* intermediate. This includes the flavin mononucleotide (FMN) in NDUV1, all seven FeS-clusters of the main electron transport pathway (N3[V1], N1b[S1], N4[S1], N5[S1], N6a[S8], N6b[S8], N2[S7]), and the off-pathway FeS cluster N1a[V2] (Figure 5A). Moreover, density can be seen in the cryoEM map in the region of the Q-tunnel, in an equivalent position to that of CoQ in the Y. lipolytica structure (Parey et al., 2019; Figure 5B). This likely represents a CoQ molecule bound at the entry of the CI* Q-tunnel. However, this density is indistinct and thus we have not modeled a CoQ at this position. Analogously to the Y. lipolytica structure, no density for CoQ can be seen deeper in the Q-tunnel where CoQ would need to bind to accept electrons from the terminal FeS cluster (Figure 5B).

Structure of the redox centers, Q cavity and the hydrophilic axis of V. radiata CI*.
(A) V. radiata’s FMN (stick) and iron-sulfur clusters (spheres) are labeled by nearest-atom center-to-center distances, overlaid with those from T. thermophilus (transparent grey). (B) Key residues (stick) delineating the Q cavity and the nearby N2 iron-sulfur cluster (spheres). Unassigned density in the Q cavity, potentially corresponding to quinone, shown as blue mesh. (C) Key CI* residues constituting the hydrophilic axis within the membrane domain shown as sticks.
The loops that cap the Q-tunnel at the interface of the peripheral and membrane arms of the complex, namely the NU3M TMH1-2 and NU1M TMH5-6 loops, are disordered. This is analogous to what is observed in the open or deactive structures of the mammalian and Y. lipolytica complexes (Agip et al., 2018; Letts et al., 2019; Parey et al., 2018). Conformational changes in these loops are thought to play an important role in CI’s coupling mechanism, which transduces the energy of NADH-quinone oxidoreduction in the Q module to proton pumping along the membrane arm (Parey et al., 2018; Cabrera-Orefice et al., 2018). In particular, a π-bulge in NU6M’s TMH3 in mammals has been seen to undergo a major conformational change, refolding into an α-helix during complex I’s open-to-closed transition (Agip et al., 2018; Letts et al., 2019). This π-bulge in NU6M’s TMH3 is also present in V. radiata CI*.
The ‘E-channel’ (Baradaran et al., 2013) and the hydrophilic axis of polar amino acid residues that are involved in proton translocation and span the membrane arm of CI are also evident in V. radiata CI* (Figure 5C). Given the lack of additional accessory subunits or assembly factors to cap the end of CI*’s shortened membrane arm, hydrophilic-axis residue Lys399 on NU2M’s TMH12 is exposed to the midplane of the membrane. In all other structures of CI, the final transmembrane core subunit NU5M contains a transmembrane helix (TMH15) that caps the hydrophilic axis at the end of the transmembrane arm of full-length CI. The lack of such a cap on V. radiata NU2M in CI* suggests that, although Lys399 of NU2M is mostly surrounded by protein, the core hydrophilic axis may be in contact with lipid.
Discussion
Protein sample
The structure of V. radiata CI* presented here is the first atomic resolution structure of any plant mitochondrial electron transport chain complex and reveals several key features of mitochondrial CI from vascular plants.
CI* is an established assembly intermediate of plant CI, previously identified with genetic and proteomic studies in non-etiolated seedlings and mature leaves of A. thaliana and N. sylvestris (Ligas et al., 2019; Meyer et al., 2011; Schertl et al., 2012; Schimmeyer et al., 2016; Senkler et al., 2017; Pineau et al., 2008). Furthermore, CI* exhibits NADH-dehydrogenase activity in in-gel activity assays (Meyer et al., 2011; Pineau et al., 2008; Haïli et al., 2013). Thus, it is unlikely that CI* in our mitochondrial preparations is a peculiarity of our etiolating growth conditions or our choice of model organism. Nevertheless, it may be the case that etiolating conditions promote the accumulation of CI* in V. radiata hypocotyls compared to seedlings grown in the light (see Appendix).
Moreover, it is also unlikely that CI* is a degradation product of CI rather than the assembly intermediate. Firstly, our membrane-extraction conditions (1% w:v digitonin, 4:1 g:g detergent:protein; see Materials and methods) are very gentle and were chosen after optimization to preserve protein:protein interactions in protein complexes and supercomplexes. Furthermore, immediately after extraction, we stabilize the detergent-extracted complexes with amphipathic polymers, which wrap around the complexes and further protect them from degradation/dissociation (Breibeck and Rompel, 2019). A large section of membrane stabilized and co-purified by our gentle digitonin/amphipol treatment is clearly seen around the perimeter of CI* at low contour (Figure 1—figure supplement 4E). Secondly, using digitonin at a higher concentration (5% w:v), an A. thaliana complexome profiling study (Senkler et al., 2017) obtained not only full-length CI and CI*, but also full-length CI in a higher order assembly with complex III (supercomplex SC I+III2) [Bultema et al., 2009; Dudkina et al., 2005; Eubel et al., 2004a; Eubel et al., 2004b]. Protein:protein interactions between complexes in supercomplexes are known to be more labile than intra-complex protein:protein interactions. Given that the more fragile CI:CIII2 interactions are maintained in 5% digitonin (Senkler et al., 2017), this argues that the presence of CI* —both in Senkler et al., 2017 and in this study— is not due to a digitonin-induced dissociation of the PD domain, but rather that it is the true assembly intermediate. Thirdly, controlled-degradation experiments of plant CI in the presence of harsh detergents have shown that, analogous to mammalian CI, plant CI’s detergent-induced dissociation occurs via detachment of the full peripheral arm (PP-PD) from the full matrix arm (N-Q) (Klodmann et al., 2010), not by dissociation between the PP and PD modules. Fourthly, we have reproducibly obtained the CI* fraction, which retains its in-gel and spectroscopic NADH-oxidase activity and chromatographic peak for several days, even after freeze/thaw cycles. For these reasons, it is evident that our structure corresponds to the CI* assembly intermediate, rather than to a degradation product of V. radiata CI.
Carbonic anhydrase domain of plant CI
A major unique feature of plant CI compared to the other known structures is the large γCA domain located on the mitochondrial matrix side of the membrane arm of the complex (Sunderhaus et al., 2006).
Here, we were able to define the interface and anchoring interactions between the γCA domain and the rest of the complex at high resolution (Figure 3). In line with expectations from the early biochemical experiments on the plant γCA domain (Sunderhaus et al., 2006), the structure clearly shows that the interface between the γCA domain and the PP-module is extensive and strong (Table 5). Additionally, we established that the γCA domain is membrane-targeted via two amphipathic helices that contact the CI membrane arm and through specific interactions with lipids in a lipid-filled pocket formed by core subunit NU2M, accessory subunits NDUX1, NDUC2 and plant-specific accessory subunit P2. Furthermore, our structure unambiguously resolves the identities of the hetero-trimeric components of the γCA domain of etiolated V. radiata as CA1, CA2 and CA2L. Unexpectedly, our structure also reveals that, due to this composition, only one out of the three potential active sites formed at the interfaces between CA1, CA2 and CA2L is capable of coordinating the Zn2+ ion required for carbonic anhydrase catalysis. Nevertheless, whether the combination of γCA subunits and, consequently, the active site arrangements are different in different species, tissues or developmental stages (Sunderhaus et al., 2006; Perales et al., 2004; Fromm et al., 2016; Cï Rdoba et al., 2019) remains to be confirmed.
Structure alone is not sufficient to demonstrate catalytic ability of the plant CI γCA domain. Indeed, only bicarbonate binding to the plant mitochondrial γCAs has been shown (Martin et al., 2009) and, despite extensive attempts, no catalytic activity has been measured to date (Fromm et al., 2016; Martin et al., 2009). Further functional and structural studies with purified CI or CI* samples are necessary to determine whether the γCA domain possesses enzymatic activity.
Structural insights on plant CI assembly
Less is known about CI assembly in plants than in fungi or metazoans (opisthokonts). In metazoans, detailed models of CI assembly have been generated and over a dozen CI assembly factors have been identified (Formosa et al., 2018; Guerrero-Castillo et al., 2017; Garcia et al., 2017). In plants, only three assembly factors have been thus far identified: L-galactono-1,4-lactone dehydrogenase (GLDH) (Senkler et al., 2017), the FeS protein INDH (Wydro et al., 2013) and an LYR protein termed CIAF1 (Ivanova et al., 2019). One possibility is that some of the unassigned densities observed in our reconstruction correspond to assembly factors that are bound to CI*. Current models of plant CI biogenesis predict that, of these three, only GLDH should be bound to the CI* intermediate (Ligas et al., 2019). However, GLDH is a large (~60 kDa) globular enzyme (Leferink et al., 2008), for which we do not see any consistent density in our structure. Nonetheless, it is possible that GLDH is bound via a flexible loop and thus averaged out in our reconstructions. Further assembly factors have been predicted to bind and cap NU2M in the membrane (Ligas et al., 2019). However, as noted above, we do not observe any additional transmembrane subunits capping the end of the shortened transmembrane arm.
There are major differences in CI assembly between plants and metazoans (Figure 1—figure supplement 1). In metazoans, the N-module (responsible for NADH oxidation) is assembled onto the Q-, PP- and PD-modules last (Formosa et al., 2018; Guerrero-Castillo et al., 2017; Garcia et al., 2017). This ensures that no assembly intermediate is capable of transferring electrons from NADH to CoQ. In contrast, in plants the final assembly step is the attachment of the PD-module onto the CI* intermediate (Ligas et al., 2019). As noted above, the V. radiata CI* intermediate contains all of the subunits and co-factors needed to carry out NADH:CoQ oxidoreduction: CI* is, in principle, catalytically competent. Indeed, we were able to measure NADH-DQ oxidoreductase activity in the isolated CI* fraction (Figure 1—figure supplement 2).
The V. radiata CI* structure presented here reveals that this difference in assembly may in part stem from a significant difference in the structure of the peripheral-arm accessory subunit NDUS6. The plant NDUS6 subunit lacks an N-terminal domain relative to the NDUS6 homologues of opisthokonts. In opisthokonts, the N-terminal domain of NDUS6 binds over top of NDUA12 to interact with the Q-module accessory subunit NDUA9 (Figure 2A). Moreover, the assembly factor NDUFAF2 –a paralogue of NDUA12 that occupies the same binding site—sterically prevents the binding of NDUS6 (Parey et al., 2019). Thus, in opisthokonts, NDUFAF2 must be removed and replaced with NDUA12 before NDUS6 can bind on the peripheral arm to complete the assembly of CI. In plants, a NDUFAF2 homologue on CI has yet to be observed experimentally (Meyer et al., 2019). Additionally, due to the lack of the N-terminal domain on NDUS6, plant NDUS6 does not cross over NDUA12 but binds next to it on the surface of the peripheral arm. Thus, in plants, NDUS6 may assemble on CI independent of the status of NDUFAF2/NDUA12. Furthermore, attaching the N-module before the PD-module in plants may provide additional flexibility to their mitochondrial ETC (see discussion below and Appendix).
It is clear from the currently available structures that the interface between the PP-module and PD-module is more extensively stabilized by accessory subunit interactions in mammals than in Y. lipolytica or V. radiata (Table 4). Although we currently only have the structure of the CI* intermediate for V. radiata (which only contains the PP-module), key truncations in core subunit NDUS2 and accessory subunits NDUA8 and NDUC2, discussed above (Figure 2B and C), already make this distinction clear. The lack of the NDUA8 and NDUC2 bridging interactions suggest that the interface between the PP- and PD-modules in plants may be weaker, which may also help explain the differences in the CI assembly pathway in plants versus opisthokonts. Identification of other possible bridging interactions across the PP- and PD-modules in plants will have to await the structure of full-length plant CI.
Potential roles for CI* beyond CI assembly
The bioenergetic regulation of plants, which generate their energy through respiration and photosynthesis, is more intricate and dynamic than that of heterotrophs, whose main bioenergetic process is respiration. Mitochondrial respiration is the major source of ATP in plants’ non-photosynthetic tissues such as roots. In photosynthetic tissue in the light, the role of mitochondrial respiration in ATP production is debated (Shameer et al., 2019; Gardeström and Igamberdiev, 2016) (see Appendix). Moreover, in photosynthetic tissue, conditions of intense light may lead to an over-production of reducing equivalents (NAD(P)H), which could be detrimental to the cells via the production of reactive oxygen species (ROS). To mitigate this, the plant mitochondrial electron transport chain (mETC) contains several ‘alternative’ oxidoreductases and oxidases that shunt electrons to molecular oxygen without pumping H+, thus preventing the over-reduction of the NADH pool (Millar et al., 2011; Schertl and Braun, 2014). However, given that alternative complexes do not pump any H+, energy is instead dissipated as heat.
Based on the fact that CI* is missing two of its four standard H+ pumps (those in the PD module), and on our finding that CI* shows NADH-DQ oxidoreduction activity (Figure 1—figure supplement 2), we hypothesize that CI* may be an NADH-CoQ oxidoreductase with a lower H+-pumping-to-electron-transfer ratio than full-length CI. Namely, we hypothesize that CI* could pump protons at a 2H+:2e- ratio rather than the 4H+:2e- of full-length CI (Jones et al., 2017).
Decreased H+:e- ratios have previously been reported in functional yeast and bacterial CI mutants (Dröse et al., 2011; Steimle et al., 2011). A mutant of Y. lipolytica CI in which the PD-module accessory subunit NB8M (homologue of plant NDUB7) is deleted (nb8mΔ) fails to assemble the PD-module (Dröse et al., 2011). The resulting CI subcomplex is analogous to CI*, as it lacks only the PD-module. The nb8mΔ mutant CI is a functional H+-pumping NADH-CoQ oxidoreductase. However, its H+:e- ratio, which is normally 4H+:2e- in fully assembled CI, is reduced to 2H+:2e- (Dröse et al., 2011). This is consistent with two of the four H+-pumping subunits (NU4M and NU5M) being absent in the nb8mΔ mutant subcomplex. Similar results are seen in E. coli mutants with mutations in its distal H+-pumping subunit NuoL (homologue of plant NU5M). Deletion of NuoL or truncation of its transmembrane helices 15–16, which bridge the PP and PD modules, result in a functional CI mutant whose H+:e- coupling is 2H+:2e- (Steimle et al., 2011).
We hypothesize that a lower-H+-pumping CI* could provide additional flexibility to plants’ bioenergetic regulation, beyond the interplay between the canonical and alternative pathways of the mETC. For instance, having a 2H+:2e- ratio would allow CI* to contribute to ATP generation in situations where the mitochondrial [NAD+]/[NADH] ratio would not support H+ pumping by CI (see Appendix for an in-depth discussion). Thus, CI* may provide additional energy-converting flexibility to plants’ electron flow and energy conservation. This would be analogous to the flexibility seen for the electron transport chain of chloroplasts, which employ several dynamic mechanisms at different levels of regulation to adjust the H+:e- coupling and the energetic and redox outputs to changing environmental conditions (Heber and Kirk, 1975; Scheibe et al., 2005; Rochaix, 2011; Murchie and Ruban, 2020).
Conclusion
Here, we present the structure of a mitochondria CI assembly intermediate, CI*, isolated from etiolated hypocotyls of V. radiata. CI* showed NADH-dehydrogenase activity in native in-gel and spectroscopic activity assays. Although we did not introduce experimental manipulations to prevent the assembly of mitochondrial CI, we were nonetheless able to isolate sufficient amounts of the CI* assembly intermediate for structure determination. This suggests that there are significant steady-state amounts of CI*in V. radiata mitochondria under these etiolating conditions and that CI* may be playing an independent physiological function beyond its role in CI assembly. The structure of V. radiata CI* presented here provides a wealth of information on mitochondrial CI composition, assembly and evolution and raises several questions on the dynamics and regulation of plant respiration. In order to address these questions, further research is needed into the structures of the fully assembled plant mitochondrial CI, as well as of its supercomplex with CIII2. In addition, biochemical, cell biological and genetic approaches are paramount to test hypotheses on the potential functions of CI*.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Biological sample (Vigna radiata) | V. radiata seeds | Todd’s Tactical Group | TS-229 | Lot SMU2-8HR; DOB 2/25/2019 |
Commercial assay or kit | Pierce BCA assay kit | Thermo Fisher | 23225 | |
Commercial assay or kit | 3–12% NativePAGE gels and buffers | Invitrogen | BN1001BOX; BN2001; BN2002 | |
Chemical compound, drug | Digitonin, high purity | EMD Millipore | 300410 | |
Chemical compound, drug | A8-35 | Anatrace | A835 | |
Chemical compound, drug | Gamma-cyclodextrin | EMD Millipore | C4892 | |
Chemical compound, drug | NADH | VWR Life Sciences | 97061–536 | |
Chemical compound, drug | Nitrotetrazoleum | EMD Millipore | 74032 | |
Software, algorithm | SerialEM | University of Colorado, Schorb et al., 2019 | RRID:SCR_017293 | |
Software, algorithm | RELION 3.0 | Zivanov et al., 2018 | RRID:SCR_016274 | |
Software, algorithm | Motioncor2 | Zheng et al., 2017 | ||
Software, algorithm | Ctffind4 | Rohou and Grigorieff, 2015 | RRID:SCR_016732 | |
Software, algorithm | crYOLO | Wagner et al., 2019; Wagner and Raunser, 2020 | RRID:SCR_016732 | |
Software, algorithm | Phyre2 | Kelley et al., 2015 | ||
Software, algorithm | Coot | Emsley and Cowtan, 2004 | RRID:SCR_014222 | |
Software, algorithm | PHENIX | Liebschner et al., 2019; Goddard et al., 2018; Pettersen et al., 2004 | RRID:SCR_014224 | |
Software, algorithm | UCSF Chimera | Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, Pettersen et al., 2004 | RRID:SCR_004097 | |
Software, algorithm | PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC. | Schrödinger, LLC | RRID:SCR_000305 | Version 2.0 |
Other | Holey carbon grids | Quantifoil | Q310CR1.3 | 1.2/1.3 300 mesh |
Vigna radiata mitochondria purification
Request a detailed protocolV. radiata seeds were purchased from Todd’s Tactical Group (Las Vegas, NV). Seeds were incubated in 1% (v:v) bleach for 20 min and rinsed until the water achieved neutral pH. Seeds were subsequently imbibed in a 6 mM CaCl2 solution for 20 hr in the dark. The following day, the imbibed seeds were sown in plastic trays on damp cheesecloth layers, at a density of 0.1 g/cm2 and incubated in the dark at 20°C for 6 days. The resulting etiolated mung beans were manually picked, and the hypocotyls were separated from the roots and cotyledons. The hypocotyls were further processed for mitochondria purification based on established protocols (Millar et al., 2007). Briefly, hypocotyls were homogenized in a Waring blender with homogenization buffer (0.4 M sucrose, 1 mM EDTA, 25 mM MOPS-KOH, 10 mM tricine, 1% w:v PVP-40, freshly added 8 mM cysteine and 0.1% w:v BSA, pH 7.8) before a centrifugation of 10 min at 1000 x g (4°C). The supernatant was collected and centrifuged for 30 min at 12,000 x g (4°C). The resulting pellet was resuspended with wash buffer (0.4 M sucrose, 1 mM EDTA, 25 mM MOPS-KOH, freshly added 0.1% w:v BSA, pH 7.2) and gently centrifuged at 1000 x g for 5 min (4°C). This supernatant was then centrifuged for 45 min at 12,000 x g. The resulting pellet was resuspended in wash buffer, loaded on to sucrose step gradients (35% w:v, 55% w:v, 75% w:v) and centrifuged for 60 min at 52,900 x g. The sucrose gradients were fractionated with a BioComp Piston Gradient Fractionator (Fredericton, Canada) connected to a Gilson F203B fraction collector, following absorbance at 280 nm. The fractions containing mitochondria were pooled, diluted 1:5 in 10 mM MOPS-KOH, 1 mM EDTA, pH 7.2 and centrifuged for 20 min at 12,000 x g (4°C). The pellet was resuspended in final resuspension buffer (20 mM HEPES, 50 mM NaCl, 1 mM EDTA, 10% glycerol, pH 7.5) and centrifuged for 20 min at 16,000 x g (4°C). The supernatant was removed, and the pellets were frozen and stored in a −80°C freezer. The yield of these mitochondrial pellets was 0.8–1 mg per gram of hypocotyl.
Vigna radiata mitochondrial membrane wash
Request a detailed protocolFrozen V. radiata mitochondrial pellets were thawed at 4°C, resuspended in 10 ml of chilled (4°C) double-distilled water per gram of pellet and homogenized with a cold Dounce glass homogenizer. Chilled KCl was added to the homogenate to a final concentration of 0.15 M and further homogenized. The homogenate was centrifuged for 45 min at 32,000 x g (4°C). The pellets were resuspended in cold Buffer M (20 mM Tris, 50 mM NaCl, 1 mM EDTA, 2 mM DTT, 0.002% PMSF, 10% glycerol, pH 7.4) and further homogenized before centrifugation at 32,000 x g for 45 min (4°C). The pellets were resuspended in 3 ml of Buffer M per gram of starting material and further homogenized. The protein concentration of the homogenate was determined using a Pierce BCA assay kit (Thermo Fisher, Waltham, MA), and the concentration was adjusted to a final concentration of 10 mg/ml and 30% glycerol.
Extraction and purification of mitochondrial complexes
Request a detailed protocolWashed membranes were thawed at 4°C. Digitonin (EMD Millipore, Burlington, MA) was added to the membranes at a final concentration of 1% (w:v) and a digitonin:protein ratio of 4:1. Membranes complexes were extracted by tumbling this mixture for 60 min at 4°C. The extract was centrifuged at 16,000 x g for 45 min (4°C). Amphipol A8-35 (Anatrace, Maumee, OH) was added to the supernatant at a final concentration of 0.2% w:v and tumbled for 30 min at 4°C, after which gamma-cyclodextrin (EMD Millipore, Burlington, MA) was added to a final amount of 1.2x gamma-cyclodextrain:digitonin (mole:mole). The mixture was centrifuged at 137,000 x g for 60 min (4°C). The supernatant was concentrated with centrifugal protein concentrators (Pall Corporation, NY, NY) of 100,000 MW cut-off, loaded onto 10–45% (w:v) or 15–45% (w:v) linear sucrose gradients in 15 mM HEPES, 20 mM KCl, pH 7.8 produced using factory settings of a BioComp Instruments (Fredericton, Canada) gradient maker and centrifuged for 16 hr at 37,000 x g (4°C). The gradients were subsequently fractionated with BioComp Piston Fractionatr connected to a Gilson F203B fraction collector, following absorbance at 280 nm. Select fractions were pooled, concentrated with protein concentrators (Pall Corporation, NY, NY) of 100,000 MW cut-off and purified on a Superose6 10/300 chromatography column (GE Healthcare, Chicago, IL) using an NGC 10 Medium-Pressure chromatography system (Biorad, Hercules, CA). For grid preparation, the relevant fractions were buffer-exchanged into 20 mM HEPES, 150 mM NaCl, 1 mM EDTA, pH 7.8 (no sucrose) and concentrated to a final protein concentration of 6 mg/ml and mixed one-to-one with the same buffer containing 0.2% digitonin (w:v),for a final concentration of 0.1% digitonin (w:v).
BN-PAGE
Request a detailed protocolMitochondrial membrane extractions were diluted in 2X BN-loading buffer (250 mM aminocaproic acid, 100 mM Tris-HCl, pH 7.4, 50% glycerol, 2.5% (w:v) Coomassie G-250), loaded on pre-cast 3–12% NativePAGE Bis-Tris gels (Invitrogen, Carlsbad, CA) and run at 4°C. The cathode buffer was 50 mM Tricine, 50 mM BisTris-HCl, pH 6.8 plus 1X NativePAGE Cathode Buffer Additive (0.02% Coomassie G-250) (Invitrogen, Carlsbad, CA) and the anode buffer was 50 mM Tricine, 50 mM BisTris-HCl, pH 6.8. Gels were run at 200 V constant voltage for ∼30 min, after which the cathode buffer was switched for a ‘light blue’ cathode buffer containing 50 mM Tricine, 50 mM BisTris-HCl, pH 6.8 plus 0.1X NativePAGE Cathode Buffer Additive (0.002% Coomassie G-250) (Invitrogen, Carlsbad, CA). The settings were changed to 7 mA constant amperage and run for another ∼90 min.
Activity assays
Request a detailed protocolThe CI in-gel NADH dehydrogenase activity assay was performed based on Schertl and Braun, 2015. The BN-PAGE gel was incubated in 10 ml of freshly prepared reaction buffer (1 mg/ml nitrotetrazoleum blue in 10 mM Tris-HCl pH 7.4). Freshly thawed NADH was added to the container with the gel, to a final concentration of 150 μM. The gel with the complete reaction buffer was rocked at room temperature for ∼10 min. Once purple bands indicating NADH-dehydrogenase activity appeared, the reaction was quenched with a solution of 50% methanol (v:v) and 10% acetic acid (v:v).
The spectroscopic NADH dehydrogenase activity assay was performed based on Huang et al., 2015; Letts et al., 2019. CI NADH:decylubiquinone (DQ) activity was measured by spectroscopic observation of NADH oxidation at 340 nm wavelength at 30°C using a Molecular Devices (San Jose, CA) Spectramax M2 spectrophotometer. Reactions were carried out in 96-well plates. Protein samples were added to 190 μL of reaction buffer (100 mM HEPES, pH 7.4, 50 mM NaCl, 10% glycerol, 4 μM KCN, 1 mg/ml BSA, 10 μM cyt c, with or without 100 μM DQ as required) and mixed by pipetting. The reaction was initiated by addition of NADH to a final concentration of 150 μM and briefly mixed by pipetting and plate stirring for 10 s before recording. Measurements were done in triplicate, averaged and background-corrected. The known extinction co-efficient of NADH (6.22 mM−1 cm−1) was used in the calculations. Statistical significance was determined using a two-tailed t-test.
CryoEM data acquisition
Request a detailed protocolThe CI* sample (6 mg/ml protein in 20 mM HEPES, 150 mM NaCl, 1 mM EDTA, 0.1% digitonin, pH 7.8) was applied onto glow-discharged holey carbon grids (Quantifoil, 1.2/1.3 300 mesh) followed by a 60 s incubation and blotting for 9 s at 15°C with 100% humidity and flash-freezing in liquid ethane using a FEI Vitrobot Mach III.
CryoEM data acquisition was performed on a 300 kV Titan Krios electron microscope equipped with an energy filter and a K3 detector at the UCSF W.M. Keck Foundation Advanced Microscopy Laboratory, accessed through the Bay Area Cryo-EM Consortium. Automated data collection was performed with the SerialEM package (Schorb et al., 2019). Micrographs were recorded at a nominal magnification of 60,010 X, resulting in a pixel size of 0.8332 Å2. Defocus values varied from 1.5 to 3.0 µm. The dose rate was 20 electrons per pixel per second. Exposures of 3 s were dose-fractionated into 118 frames, leading to a dose of 0.72 electrons per Å2 per frame and a total accumulated dose of 86.4 electrons per Å2. A total of 9816 micrographs were collected, 8541 of which were used for further analysis.
Data processing
Request a detailed protocolSoftware used in the project was installed and configured by SBGrid (Morin et al., 2013). All processing steps were done using RELION 3.0 (Zivanov et al., 2018) unless otherwise stated. Motioncor2 (Zheng et al., 2017) was used for whole-image drift correction of each micrograph. Contrast transfer function (CTF) parameters of the corrected micrographs were estimated using Ctffind4 (Rohou and Grigorieff, 2015) and refined locally for each particle in RELION. Automated particle picking using crYOLO (Wagner et al., 2019; Wagner and Raunser, 2020) resulted in ~1.5 million particles. The particles were extracted using 4002 pixel box binned two-fold and sorted by reference-free 2D classification followed by re-extraction at 5122 pixel box. Reference-free 2D classification resulted in the identification of 190,951 CI* particles. An ab initio model was generated in RELION from these particles (Punjani et al., 2017). This model, lowpass-filtered at 30 Å, was used for initial 3D classification with a regularization parameter T of 4. This initial processing resulted in ~34,000 particles of good quality, which separated into a single class (Figure 1—figure supplement 3C). The best class was refined to a nominal resolution of 3.9 Å according to the gold standard FSC criteria (Scheres and Chen, 2012). It was clear that the local resolution of this refinement was impacted by hinge-like motions between the membrane and peripheral arms of the complex. Therefore, sub-region refinements were also performed masking around the membrane arm and peripheral arm, respectively (Figure 1—figure supplement 3C). This resulted in significantly, improved map quality, especially for the γCA domain on the membrane arm (Figure 1—figure supplement 3C). These improved maps were used for model building and refinement. The two focused refined maps were then combined into a composite map using Phenix.
Model building and refinement
Request a detailed protocolStarting models for isolated ovine CI (Letts and Sazanov, 2015) and bacterial γCA (Iverson et al., 2000), corrected for the V. radiata sequence, were used as templates. Additionally, starting models were generated using the Phyre2 web portal (Kelley et al., 2015). These models were split and fit into the highest-resolution focused refinement maps for separate atomic model building of the CI* peripheral arm and CI* membrane arm in Coot (Emsley and Cowtan, 2004). Real-space refinement of the model was done in PHENIX (Liebschner et al., 2019; Goddard et al., 2018; Pettersen et al., 2004) and group atomic displacement parameters (ADPs) were refined in reciprocal space. The single cycle of group ADP refinement was followed by three cycles of global minimization, followed by an additional cycle of group ADP refinement and finally three cycles of global minimization (Letts et al., 2019).
Model interpretation and figure preparation
Request a detailed protocolMolecular graphics and analyses were performed with UCSF Chimera (Pettersen et al., 2004), developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311, as well as the PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC.
Appendix 1
Bioenergetic Considerations of CI* H+ pumping
1.Background
The purpose of this Appendix is to discuss the bioenergetic implications of CI*’s potential function as a 2e-:2H+ redox-coupled proton pump and what, if any, possible roles it could play in plant bioenergetic physiology.
For plant cells, there is much debate about the roles of the various pathways and organelles in supplying energy to the cytoplasm and reduction equivalents to the peroxisome (O'Leary et al., 2019; Gardeström and Igamberdiev, 2016). Energy is provided to the cytoplasm and peroxisome in the form of ATP and/or reduction equivalents (NAD(P)H), from the chloroplasts and mitochondria through many pathways. Pathways that export ATP or reduction equivalents from the chloroplast to the cytoplasm include: 1) the chloroplast malate valve; 2) the triose phosphate-3-phosphglycerate (TP-2PGA) shuttle and 3) the phosphoenolpyruvate (PEP)-pyruvate shuttle. ATP can also be imported into the chloroplast from the cytosol via the plastidial nucleotide translocase to energize the chloroplast at night (Flügge et al., 2011). Pathways that export ATP or reduction equivalents directly from the mitochondria to the cytoplasm include: 1) the mitochondrial malate valve; 2) the mitochondrial adenylate nucleotide translocase.
Furthermore, photosynthesis and mitochondrial respiration are linked via the photorespiratory C2 cycle and the peroxisome (Husic et al., 1987). In this pathway, 2-phosphoglycerate generated by the oxygenase activity of RuBisCO is converted to glycolate and transported into the peroxisome, where it is converted into glycine. Glycine is then transported into the mitochondria, where it is converted by the glycine decarboxylase complex (GDC) into serine, CO2 and NH3+, reducing NAD+ to NADH in the mitochondrial matrix. These reduction equivalents can then be exported to the cytoplasm by the mitochondrial malate valve or be fed into oxidative phosphorylation via complex I (CI). In order to complete the photorespiratory pathway, serine must be transported back to the peroxisome, where is it reduced to glycerate, which is then transported back into the chloroplast. The reduction equivalents for the conversion of serine to glycerate in the peroxisome are provided from the cytoplasm by the malate-oxaloacetate shuttle of the peroxisome.
Recent detailed modelling of the energetic coupling between organelles in plant cells revealed that, although chloroplasts can theoretically generate sufficient ATP to satisfy the energy requirements of the entire plant cell, this would require unrealistic light-use efficiency and higher-than-available levels of ATP export from the chloroplast (Shameer et al., 2019). Although still controversial, these modeling results in conjunction with some experimental results (reviewed in Gardeström and Igamberdiev, 2016) suggest that during photosynthesis the bulk of cytosolic ATP is provided by the mitochondria. Furthermore, these studies suggest that rather than being exported for use by the peroxisome in photorespiration, the reduction equivalents (NADH) generated by GDC in the mitochondrial matrix are used directly for ATP production by the mitochondrial electron transport chain (mETC) and the needs of the peroxisome are mainly met by chloroplast-derived reduction equivalents (Shameer et al., 2019). Moreover, it is important to note that the photorespiratory C2 cycle is a high-flux pathway in C3 plants such as V. radiata, with a flux approximately equal on a molar basis to the flux through the photosynthetic C3 cycle (Oliver, 1994). The total CO2 released by GDC is ~25% of the moles of C fixed by RuBisCO (Husic et al., 1987; Oliver, 1994).
Given the above considerations, it is reasonable to conclude that a major source of NADH for mitochondrial respiration in actively photosynthesizing cells is generated by GDC and that these levels may fluctuate as a function of photosynthetic output, driven by light availability.
In this Appendix, we consider the effect of the fluctuations of the mitochondrial-matrix [NAD+]/[NADH] ratio on the reaction catalyzed by CI. We then discuss how alternative metabolic routes in the mETC, such as alternative NADH dehydrogenases (NDs), the alternative oxidase and possibly CI*, may provide plants with respiratory flexibility to manage a range of different and changing conditions (see also the review by O'Leary et al., 2019).
2.Complex I reaction
Mitochondrial CI catalyzes the reaction:
Where:
NADH is the reduced form of nicotinamide adenine dinucleotide
CoQ is the oxidized form of coenzyme Q (ubiquinone)
represents a proton on the negative (N) side of the membrane (mitochondrial matrix)
is the number of H+ pumped across the inner mitochondrial membrane by CI
NAD+ is the oxidized from of nicotinamide adenine dinucleotide
CoQH2 is the reduced from of coenzyme Q (ubiquinol)
represents a proton on the positive (P) side of the membrane (inter-membrane space)
The Gibbs energy change (ΔGCI) of the CI reaction can be determined by splitting the reaction into its separate electron transfer and H+-pumping parts. For completeness, we will briefly derive the expression for these two parts here and then combine them into the final expression for ΔGCI.
2.1 Electron Transfer
The above oxidoreduction reactions for NADH and CoQ can be represented by two half reactions:
The midpoint potential at which the concentrations of the reduced and oxidized forms are equal at pH 7.0 () for these half reactions are known to be -320 mV for Reaction 2 and 4 mV for Reaction 3 (see Appendix 1—table 1 for references).
Values used in the calculations.
Variable | Value | Source |
---|---|---|
160 mV | Values of 140–190 mV have been reported from respiring cells (Ripple et al., 2013); a value of 200 mV was reported for isolated etiolated V. radiata mitochondria after addition of 1 mM NADH, which defines an upper limit for steady-state respiration (Moore and Bonner, 1981) | |
R | 8.314 kJ K−1 mol−1 | Physical Constant |
T | 300 K | Approximately 27 °C |
F | 96,485 C mol−1 | Physical Constant |
4 mV | This value varies as a function of pH so should only be considered an estimate (Nicholls, 2013) | |
10 | Kim et al., 2012 | |
34 mV | Calculated from and using Equation 1 | |
−320 mV | This value varies as a function of pH, so should only be considered an estimate (Nicholls, 2013) |
The redox potential of the half reactions at pH 7 can be calculated using the following equation:
Where:
is the redox potential
is the midpoint potential
R is the gas constant (8.314 kJ K−1 mol−1)
The factor of 2.3 originates from converting the natural logarithm to log10
T is the absolute temperature (K)
is the number of electrons transferred in the half reaction
F is Faraday’s constant (96,485 C mol−1)
[oxidized] is the actual concentration of the oxidized form
[reduced] is the actual concentration of the reduced form
The redox potential difference between the NADH and CoQ pools is defined as the difference in their redox potential:
Where:
is the redox potential difference
is the redox potential for CoQ
is the redox potential for NADH
as presented in Equation 2 is also known as the redox span of CI (). The redox span of CI is related to the Gibbs energy change accompanying the electron transfer (ΔGET) between the couples by:
Where:
ΔGET is the Gibbs energy change of the electron transfer
2 is the number of electrons transferred
F is Faraday’s constant (96,485 C mol−1)
2.2 Proton pumping
In the general case for the Gibbs energy change (ΔG) accompanying the transport of an ion across a membrane, the ion will be affected by both concentrative and electrical gradients:
Where:
is the charge of the ion
F is Faraday’s constant (96,485 C mol−1)
ΔΨ is the membrane potential
R is the gas constant (8.314 kJ K−1 mol−1)
T is the absolute temperature (K)
[Xm+]P is the concentration of ions on the P side of the membrane
[Xm+]N is the concentration of ions on the N side of the membrane
This is often expressed as the ion electrochemical gradient with units of kJ mol−1. For a proton electrochemical gradient , Equation 4 can be simplified as pH is a logarithmic function of [H+]:
Where:
ΔpH is defined as the pH on the P side of the membrane minus the pH on the N side (pHP-pHN)
The factor of 2.3 comes from converting the natural logarithm to log10
The proton motive force (PMF or Δp) was defined by Peter Mitchell to convert into units of voltage to facilitate comparison with redox potential differences:
Where:
F is Faraday’s constant (96,485 C mol−1)
2.3 Combining terms for overall ΔG expression
Given Equations 3 and 5, the overall ΔG of Reaction 1 catalyzed by CI can be given as:
Using Equation 6 we obtain:
Using Equation 2 we obtain:
Finally, using Equation 1 we obtain:
Equation 10 allows us to express the Gibbs energy change for CI at a given proton motive force ( and redox poise of the CoQ pool as a function of the number of H+ pumped and the ratio of NAD+ to NADH in the mitochondrial matrix (Appendix 1—table 1, Appendix 1—figure 1).

The Gibbs energy change of the CI reaction (ΔGCI) as a function of the redox poise of the mitochondrial NADH pool.
The Gibbs energy change was calculated using equation 10 and the values presented in Table A1, for reactions in which CI pumps 4 H+ (blue; representative of the standard, full-length CI pumping with a 4H+:2e- ratio) or 2 H+ (red; representative of a putative CI* pumping with a 2H+/2e- ratio). The horizontal dashed line indicates equilibrium state (ΔGCI = 0) for the different [NAD+]/[NADH] ratios. The vertical dashed line indicates the [NAD+]/[NADH] ratio at which full-length CI (blue) attains equilibrium ( = 1). The highlighted orange region corresponds to conditions in which thermodynamics would favor reverse electron transport (RET) by full-length CI ( > 1).
At negative values of ΔGCI, the CI reaction occurs in the forward direction, that is, NADH oxidation and proton pumping into the intermembrane space. At positive values of ΔGCI, the CI reaction occurs in the reverse direction, usually associated with large amounts of reactive oxygen species (ROS) production, which are known to be detrimental to the cell (Robb et al., 2018). This reverse reaction is also called ‘reverse electron transport’ (RET).
2.4. Thermodynamic efficiency
The thermodynamic efficiency () of the reaction is defined as the fraction of the energy released on electron transfer that is transduced into the proton motive force:
By definition, < 1 if the CI reaction is in the forward direction (i.e., oxidation of NADH), = 1 at equilibrium (where ) and > 1 during RET by CI (i.e. ROS production).
3.Discussion
Although the energetic analysis presented here does not take into account the dynamics of the electron transport system (i.e., fluctuations in and caused by fluctuations in [NAD+]/[NADH]), several important conclusions can be drawn. (For clarity, we call a CI pumping with a 4H+:2e-ratio simply ‘CI’. For the sake of argument, we assume that CI* is a CI entity that pumps protons with a 2H+:2e- ratio.):
When CI is operating near equilibrium, CI* runs irreversibly in the forward direction. CI* would be incapable of operating in ROS-generating RET mode in a range of [NAD+]/[NADH] up to values of ~105.
Under conditions that favor RET from CI (e.g. a drop in [NADH], leading to an increased [NAD+]/[NADH] ratio), CI* continues to work in the forward direction. This could help maintain the inner mitochondrial membrane proton motive force () at high [NAD+]/[NADH].
CI has a higher thermodynamic efficiency than CI* under all conditions that favor the forward direction for CI (ΔGCI < 0). Nevertheless, under these conditions, CI* would still transduce energy from NADH into . Due to the factor of 2 difference in H+-pumping, CI* would always transduce half the amount of energy compared to CI (see Equation 11). For example, under conditions of near equilibrium for CI (), CI* would still transduce energy at 50% efficiency ().
From this analysis, we conclude that one of the key advantages of having a CI* pumping at 2H+:2e- would be that it could still work in the forward direction, maintaining the inner mitochondrial membrane’s proton motive force, in situations where full-length CI would operate in the reverse direction and generate ROS. However, in order for plants to take full advantage of this potential bioenergetic benefit of partial energy transduction by CI*, plant CI should display a strong rectification that prevents RET. In other words, there should be a mechanism that strongly inhibits CI turnover in conditions where it would otherwise run in reverse and lead to oxidative damage.
This type of rectification of CI has been proposed to exist in mammals via CI’s active-to-deactive transition as a way to prevent oxidative damage upon ischemic reperfusion (Chouchani et al., 2013). The active-to-deactive transition of mitochondrial CI has been observed in several but not all studied fungi and metazoans; moreover, it is absent in all prokaryotic CI thus far examined (reviewed in Babot et al., 2014). It is currently unknown whether plant CI displays an active-to-deactive transition.
The alternative NDs irreversibly operate in the forward directly over the large range of [NAD+]/[NADH] in which CI operates in reverse. As discussed above, this would also be the case for a CI* that pumped protons at 2H+:2e-. The simultaneous activity of alternative NDs and CI would continuously push the [NAD+]/[NADH] ratio towards the CI RET regime, due to the irreversible oxidation of NADH and reduction of CoQ by the NDs. Thus, the potential existence of a 2H+:2e--pumping CI* does not generate additional bioenergetic problems beyond those already created by the existence of the alternative NDs (which do not pump any protons at all). The plant cell must already have regulatory mechanisms to deal with the threat of RET by CI imposed by the NDs. The degree to which these alternative NDs are employed and regulated in vivo remains poorly understood (O'Leary et al., 2019). We predict that some type of rectification operates on plant CI as a mechanism to prevent ROS production under any conditions that favors RET by CI.
This analysis also proposes a possible answer to why our preparations of etiolated V. radiata contain such a significant amount of CI*, compared to the previously reported lower abundance of CI* in non-etiolated tissues (Ligas et al., 2019; Senkler et al., 2017). To the best of our knowledge, the mitochondrial [NAD+]/[NADH] ratio of etiolated hypocotyls has not been investigated. However, given the lack of input of reducing equivalents by the C2 cycle via GDH in the dark (an otherwise high-flux pathway), it is conceivable that the [NAD+]/[NADH] ratio in etiolated hypocotyls is higher than in photosynthesizing cells. A high [NAD+]/[NADH] ratio may favor the use of CI* over CI in order to ensure maintenance of the proton motive force (), at the expense of thermodynamic efficiency. It is conceivable that, as hypocotyls develop under etiolating conditions and their only source of energy (i.e. the seed oils) diminishes, the ratio of CI* to CI present in the mitochondrial membranes may be dynamically regulated to increase CI* levels.
Although CI*’s proton-pumping ratio remains to be characterized, the theoretical analysis above suggests that that a 2H+:2e--pumping entity may be beneficial for plants’ bioenergetic flexibility if a rectification mechanism for CI exists in plants. Further studies are needed to test these hypotheses.
Data availability
Deposition of structural models and maps to the PDB and EMDB are: EMD-22090, PDB ID 6X89, EMD-22091, EMD-22092, EMD-22093.
-
RCSB Protein Data BankID 6X89. Vigna radiata mitochondrial complex I*.
-
Electron Microscopy Data BankID EMD-22090. Vigna radiata mitochondrial complex I*.
-
Electron Microscopy Data BankID EMD-22091. Vigna radiata mitochondrial complex I*.
-
Electron Microscopy Data BankID EMD-22092. Vigna radiata complex I* membrane arm.
-
Electron Microscopy Data BankID EMD-22093. Vigna radiata complex I* peripheral arm.
References
-
Cryo-EM structures of complex I from mouse heart mitochondria in two biochemically defined statesNature Structural & Molecular Biology 25:548–556.https://doi.org/10.1038/s41594-018-0073-1
-
Characterisation of the active/de-active transition of mitochondrial complex IBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1837:1083–1092.https://doi.org/10.1016/j.bbabio.2014.02.018
-
The life of plant mitochondrial complex IMitochondrion 19:295–313.https://doi.org/10.1016/j.mito.2014.02.006
-
Successful amphiphiles as the key to crystallization of membrane proteins: bridging theory and practiceBiochimica Et Biophysica Acta (BBA) - General Subjects 1863:437–455.https://doi.org/10.1016/j.bbagen.2018.11.004
-
Megacomplex organization of the oxidative phosphorylation system by structural analysis of respiratory supercomplexes from potatoBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1787:60–67.https://doi.org/10.1016/j.bbabio.2008.10.010
-
Respiratory chain supercomplexes in the plant mitochondrial membraneTrends in Plant Science 11:232–240.https://doi.org/10.1016/j.tplants.2006.03.007
-
Structure and function of mitochondrial supercomplexesBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1797:664–670.https://doi.org/10.1016/j.bbabio.2009.12.013
-
The coupling mechanism of respiratory complex I — A structural and evolutionary perspectiveBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1817:1785–1795.https://doi.org/10.1016/j.bbabio.2012.02.015
-
Coot: model-building tools for molecular graphicsActa Crystallographica. Section D, Biological Crystallography 60:2126–2132.https://doi.org/10.1107/S0907444904019158
-
Identification and characterization of respirasomes in potato mitochondriaPlant Physiology 134:1450–1459.https://doi.org/10.1104/pp.103.038018
-
Respiratory supercomplexes in plant mitochondriaActa Physiologiae Plantarum 26:86–87.https://doi.org/10.1134/S1021443714060168
-
The gamma class of carbonic anhydrasesBiochimica Et Biophysica Acta (BBA) - Proteins and Proteomics 1804:374–381.https://doi.org/10.1016/j.bbapap.2009.08.026
-
The role of transporters in supplying energy to plant plastidsJournal of Experimental Botany 62:2381–2392.https://doi.org/10.1093/jxb/erq361
-
Building a complex complex: assembly of mitochondrial respiratory chain complex ISeminars in Cell & Developmental Biology 76:154–162.https://doi.org/10.1016/j.semcdb.2017.08.011
-
Assembly of the Escherichia coli NADH:ubiquinone oxidoreductase (respiratory complex I)Biochimica Et Biophysica Acta (BBA) - Bioenergetics 1857:214–223.https://doi.org/10.1016/j.bbabio.2015.12.004
-
The carbonic anhydrase domain of plant mitochondrial complex IPhysiologia Plantarum 157:289–296.https://doi.org/10.1111/ppl.12424
-
The origin of cytosolic ATP in photosynthetic cellsPhysiologia Plantarum 157:367–379.https://doi.org/10.1111/ppl.12455
-
The pentatricopeptide repeat MTSF1 protein stabilizes the nad4 mRNA in Arabidopsis mitochondriaNucleic Acids Research 41:6650–6663.https://doi.org/10.1093/nar/gkt337
-
Mitochondrial complex I from Arabidopsis and rice: orthologs of mammalian and fungal components coupled with plant-specific subunitsBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1604:159–169.https://doi.org/10.1016/S0005-2728(03)00045-8
-
Flexibility of coupling and stoichiometry of ATP formation in intact chloroplastsBiochimica Et Biophysica Acta (BBA) - Bioenergetics 376:136–150.https://doi.org/10.1016/0005-2728(75)90212-1
-
Mitochondrial complex IAnnual Review of Biochemistry 82:551–575.https://doi.org/10.1146/annurev-biochem-070511-103700
-
Activity assay for plant mitochondrial enzymesMethods in Molecular Biology 1305:139–149.https://doi.org/10.1007/978-1-4939-2639-8_10
-
The oxidative photosynthetic Carbon-Cycle or C2 cycleCrc Critical Reviews in Plant Sciences 5:45–100.https://doi.org/10.1080/07352688709382234
-
A mitochondrial LYR protein is required for complex I assemblyPlant Physiology 181:1632–1650.https://doi.org/10.1104/pp.19.00822
-
Observation of a calcium-binding site in the gamma-class carbonic anhydrase from Pyrococcus horikoshiiActa Crystallographica Section D Biological Crystallography 64:1012–1019.https://doi.org/10.1107/S0907444908024323
-
Respiratory complex I in Bos taurus and Paracoccus denitrificans Pumps Four Protons across the Membrane for Every NADH OxidizedJournal of Biological Chemistry 292:4987–4995.https://doi.org/10.1074/jbc.M116.771899
-
The Phyre2 web portal for protein modeling, prediction and analysisNature Protocols 10:845–858.https://doi.org/10.1038/nprot.2015.053
-
"Respirasome"-like supercomplexes in green leaf mitochondria of spinachJournal of Biological Chemistry 279:48369–48375.https://doi.org/10.1074/jbc.M406085200
-
Inference of macromolecular assemblies from crystalline stateJournal of Molecular Biology 372:774–797.https://doi.org/10.1016/j.jmb.2007.05.022
-
Gaining mass: the structure of respiratory complex I-from bacterial towards mitochondrial versionsCurrent Opinion in Structural Biology 33:135–145.https://doi.org/10.1016/j.sbi.2015.08.008
-
Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in phenixActa Crystallographica Section D Structural Biology 75:861–877.https://doi.org/10.1107/S2059798319011471
-
The assembly pathway of complex I in Arabidopsis thalianaThe Plant Journal : For Cell and Molecular Biology 97:447–459.https://doi.org/10.1111/tpj.14133
-
Insights into the composition and assembly of the membrane arm of plant complex I through analysis of subcomplexes in Arabidopsis mutant linesJournal of Biological Chemistry 286:26081–26092.https://doi.org/10.1074/jbc.M110.209601
-
Proteomic investigations of complex I composition: how to define a subunit?Frontiers in Plant Science 3:106.https://doi.org/10.3389/fpls.2012.00106
-
Assembly of the complexes of the oxidative phosphorylation system in land plant mitochondriaAnnual Review of Plant Biology 70:23–50.https://doi.org/10.1146/annurev-arplant-050718-100412
-
Isolation and subfractionation of mitochondria from plantsMethods in Cell Biology 80:65–90.https://doi.org/10.1016/S0091-679X(06)80003-8
-
Organization and regulation of mitochondrial respiration in plantsAnnual Review of Plant Biology 62:79–104.https://doi.org/10.1146/annurev-arplant-042110-103857
-
A comparison of the phosphorylation potential and electrochemical proton gradient in mung bean mitochondria and phosphorylating sub-mitochondrial particlesBiochimica Et Biophysica Acta (BBA) - Bioenergetics 634:117–128.https://doi.org/10.1016/0005-2728(81)90132-8
-
The glycine decarboxylase complex from plant mitochondriaAnnual Review of Plant Physiology and Plant Molecular Biology 45:323–337.https://doi.org/10.1146/annurev.pp.45.060194.001543
-
Gamma carbonic anhydrase like complex interact with plant mitochondrial complex IPlant Molecular Biology 56:947–957.https://doi.org/10.1007/s11103-004-6324-z
-
UCSF chimera--a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
L-galactono-1,4-lactone dehydrogenase is required for the accumulation of plant respiratory complex IJournal of Biological Chemistry 283:32500–32505.https://doi.org/10.1074/jbc.M805320200
-
Mammalian complex I pumps 4 protons per 2 electrons at high and physiological proton motive force in living cellsJournal of Biological Chemistry 288:5374–5380.https://doi.org/10.1074/jbc.M112.438945
-
Control of mitochondrial superoxide production by reverse electron transport at complex IJournal of Biological Chemistry 293:9869–9879.https://doi.org/10.1074/jbc.RA118.003647
-
Regulation of photosynthetic electron transportBiochimica Et Biophysica Acta (BBA) - Bioenergetics 1807:375–383.https://doi.org/10.1016/j.bbabio.2010.11.010
-
CTFFIND4: fast and accurate defocus estimation from electron micrographsJournal of Structural Biology 192:216–221.https://doi.org/10.1016/j.jsb.2015.08.008
-
A giant molecular proton pump: structure and mechanism of respiratory complex INature Reviews Molecular Cell Biology 16:375–388.https://doi.org/10.1038/nrm3997
-
Strategies to maintain redox homeostasis during photosynthesis under changing conditionsJournal of Experimental Botany 56:1481–1489.https://doi.org/10.1093/jxb/eri181
-
Prevention of overfitting in cryo-EM structure determinationNature Methods 9:853–854.https://doi.org/10.1038/nmeth.2115
-
L-galactono-1,4-lactone dehydrogenase (GLDH) forms part of three subcomplexes of mitochondrial complex I in Arabidopsis thalianaJournal of Biological Chemistry 287:14412–14419.https://doi.org/10.1074/jbc.M111.305144
-
Respiratory electron transfer pathways in plant mitochondriaFrontiers in Plant Science 5:163.https://doi.org/10.3389/fpls.2014.00163
-
Activity measurements of mitochondrial enzymes in native gelsMethods in Molecular Biology 1305:131–138.https://doi.org/10.1007/978-1-4939-2639-8_9
-
The mitochondrial complexome of Arabidopsis thalianaThe Plant Journal : For Cell and Molecular Biology 89:1079–1092.https://doi.org/10.1111/tpj.13448
-
Identification and preliminary characterization of two cDNAs encoding unique carbonic anhydrases from the marine alga Emiliania huxleyiApplied and Environmental Microbiology 72:5500–5511.https://doi.org/10.1128/AEM.00237-06
-
Carbonic anhydrase subunits form a matrix-exposed domain attached to the membrane arm of mitochondrial complex I in plantsJournal of Biological Chemistry 281:6482–6488.https://doi.org/10.1074/jbc.M511542200
Decision letter
-
David M KramerReviewing Editor; Michigan State University, United States
-
Cynthia WolbergerSenior Editor; Johns Hopkins University School of Medicine, United States
-
Toshiharu ShikanaiReviewer
-
Nicholas FisherReviewer; Michigan State University, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
This structural work represents a significant advance in our understanding of the mechanism and regulation of mitochondrial Complex I, particularly in plants, which have been relatively little studied. Further, the system that was developed will serves the field more generally by providing a useful platform for exploring a broad range of open questions about the function, regulation, assembly and Complex I in plants.
Decision letter after peer review:
Thank you for submitting your article "Atomic Structure of a Mitochondrial Complex I Intermediate from Vascular Plants" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Cynthia Wolberger as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Toshiharu Shikanai (Reviewer #2); Nicholas Fisher (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission. Your manuscript is of interest, but as described below that additional experiments are required before it is published. We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is “in revision at eLife”. Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)
Summary:
The structural work is of high quality and has interesting, possibly strong, implications for the mechanism of Complex I (CI) and its roles in diverse systems. The work focuses on a subcomplex, CI*, that is thought to be an assembly intermediate, and is missing the distal subcomplex of the proton-pumping component of the complex. An interesting hypothesis is presented in which CI* is functional, but has lower capacity to pump protons, allowing the mitochondrion to modulate energy coupling.
Essential revisions:
All reviewers had concerns about the relevance of the CI* subcomplex, and felt that it is critical to assess, at least to some extent, whether it is functional. Indeed, the title of the paper and part of the discussion assumes that CI* is an assembly intermediate, while much of the Discussion centers around the possibility that CI* is functional, but with a lower H+ stoichiometry. The manuscript does not provide any support for either of these roles.
While the reviewers recognize the difficulty in obtaining new data in the current COVID-19 crisis, knowing if the subcomplex is a functional component, a true assembly intermediate, both or an artifact of isolation is essential for the interpretation of much of the results.
Reviewers #2 and 3 provide straightforward approaches to test for these roles. First, it was thought important to address the question of whether the CI* subcomplex is an artifact of isolation, and if more intact complexes were isolated in other tissues. Second, it is essential to determine if the CI* complex has enzymatic activity. Although not considered to be essential, the paper would have greater impact if it was possible to demonstrate a change in H+/e- pumping stoichiometry. The title and Discussion should be modified to in response to these additional data.
Reviewer #1:
This manuscript describes a cryo-EM structure of a portion of plant mitochondrial Complex I (1). This subcomplex, termed CI*, is thought to represent an intermediate in assembly.
The work shows key details about some of the structural distinctions between plant and non-plant CI, including the additional subunit related to carbonic anhydrase (CA). The analysis of the structure leads to some interesting discussion, but it lacks firm hypothesis testing that would allow for stronger conclusions.
1) "Whereas the exact γCA protein combinations are likely tissue- and development-stage-dependent (Cordoba et al., 2019), the role of the γCA domain in plant CI's function is unknown (Martin et al., 2009)."
How does this relate to the CA thought to be involved in the CCM in cyanobacteria?
2) How do we know if this is really the CI* intermediate and not some sort of breakdown product?
3) I am confused by this:
“At each active site, two histidine residues are provided by one subunit and the third is provided by the adjacent subunit. However, in the plant CI γCA heterotrimer, the γCA2L subunit is lacking two of the three essential histidine residues (Ala-147 and Arg-152 in V. radiata) that would be necessary to form active sites at the interfaces with the CA1 and CA2 subunits. This renders two of the possible three catalytic sites non-functional (Figure 3A, Figure 3—figure supplement 2) …Our subunit assignments are consistent with this nonfunctional residue on CA1 being at the interface with CA2L, which would already be inactive due to the lack of histidine on CA2L at position."
It is unclear is the structure provides independent support that there is one of three active sites, or is their assignment based on this assumption?
In other words, there are hypotheses that the proteins represent one-of-three active CA or are completely inactive for CA activity. To what extent can the new structure actually test these possibilities? The test states that the Zn2+ was observed only one site. How strong is this evidence and how strongly does this support the one-in-three CA model? In this regard is it far more important what is inconsistent with the structure than what can be made to be consistent with it. In other words, if the structure can eliminate all but one of these hypotheses, it would be a big advance.
4) "Therefore, carbonic anhydrase activity of the plant γCA domain on CI must be confirmed experimentally." Why not test for activity in the amazing new prep?
5) "The lack of such a cap on V. radiata NU2M in CI* suggests that, although the Lys399 of NU2M is mostly surrounded by protein, the core hydrophilic axis may be in contact with lipid." Alternatively, some subunits may be missing or difficult to resolve in the prep.
6) "This leaves open the possibility that the γCA domain does not interconvert CO2 and bicarbonate but acts only as a sensor of CO2 or bicarbonate concentration. In this scenario, conformational changes induced by bicarbonate binding could propagate into the membrane arm of CI and potentially regulate the catalytic turnover of CI. Such sensory and regulatory roles have also been proposed for other CI accessory subunits in other organisms."
Is there any evidence from the structure in support of are inconsistent with the proposed conformation propagation mechanism for CO2 sensing? Perhaps I am missing something, but the CA subunit seems to be rather loosely tethered, and it is difficult to see from the images provided how small conformational changes at the metal site would strongly affect the structure.
Also, given that this is a major point of the discussion, elaborate on why such signaling is interesting and what does the structure say about it. In particular, see points above on whether the structure confirms or not the proposed lack of catalytic activity.
7) The remaining differences between this and other structures are framed as important for assembly. Yet, how is it known that this is a real assembly intermediate (see point 2)?
8) I am totally not convinced that there is any evidence from the structure one way or the other for a variable H+/e- transmission. Please convince me.
9) "In heterotrophs that such as metazoans that obtain all their reduction equivalents from the breakdown of sugars, fats and proteins and that depend almost entirely on oxidative phosphorylation for energy production, there would be no advantage to being able to adjust CI's H+-pumping-to-electron-transfer ratio." This statement seems to be contradicted by the text that follows, which argues that adjusting H+/e- stoichiometry might be important for maintaining PMF (though I am not sure I agree with those).
10) Regarding the paragraph from "In autotrophs such as plants, which generate their energy through a combination of respiration and photosynthesis, the bioenergetic processes necessarily have to be more dynamic to balance production of energy (ATP) and of reducing equivalents (NADH) throughout the cell and the organism."
This is highly speculative, which in itself is OK, but I'm not sure what, if anything in the structure, supports this view.
11) The authors are correct that the energy balance of the cell will be strongly impacted by activating photosynthesis. However, there are several factors that need to be considered: 1) It is known for some time that respiration is down-regulated during photosynthesis, so that little if any ATP is produced by the mitochondria; 2) The critical factor is often the balance of energy in the forms of ATP/ADP+Pi and NAD(P)H/NADP+ rather than the total energy availability. In this case, a variable proton pumping stoichiometry might be useful, but; 3) There are already alternative mechanisms of shunting electrons to O2 without production of ATP, e.g. though the alternative oxidases or Water-water cycles.
12) "We posit that, in addition to bypassing H+-pumping altogether through the alternative complexes, it is conceivable that plants may find additional advantages in regulating the H+-pumping-to-electron-transfer ratio in the "canonical" electron transport chain. If the NADH/NAD+ ratio dropped in photosynthetic cells due to decreased light availability, this would lead to a decrease in free energy available to CI for H+ pumping.
13) "Full-length CI must pump 4 protons per 2 NADH electrons, it cannot operate in a graduated manner. Consequently, at a decreased NADH/NAD+ ratio, with lower free energy available, full-length CI would initially not be able to pump any protons at all. This would result in a reduction in the mitochondrial membrane potential, decreasing the proton-motive force and lowering ATP generation.
First, I think this point, if a full-length CI can pump lower Hstoichiometry of H+/e- is debatable (and being debated), so one cannot conclude that it cannot.
Second, the logic connecting the two sentences above seems to be incorrect. The smaller the coupling stoichiometry, the higher (not lower) the PMF that can be maintained at static head for a given NADH/NAD+ ratio.
CI is a reversible proton pumping-NADH:UQ oxidoreductase, so it will reach quasi-equilibrium (static head) when the free energy in the reaction
NADH + UQ + nH(n)+ ⇓ ◊ NAD+ UQH2 + nH+(p)
equals that in PMF.
The text then states: "however, if the plant ETC could find a way to keep pumping protons even at a decreased NADH/NAD+ ratio…" This would seem to be self-contradictory. One possibility is that the authors intend to include the total proton translocation of the system, and not just "pumping". In this case, PMF includes contributions from the chemical protons that are taken up during UQ reduction, which are likely to be deposited on the p-side of the membrane during UQH2 oxidation etc, e.g. at the Q-cycle.
The text continues, "however, if the plant ETC could find a way to keep pumping protons even at a decreased NADH/NAD+ ratio, this would allow for better maintenance of the proton motive force." The authors are apparently conflating near equilibrium and steady-state systems. With a smaller coupling stoichiometry, the number of protons pumped per NADH oxidized decreases, but the extent of PMF that can be generated at near equilibrium conditions increases. In a situation where the PMF is far out of equilibrium from the free energy in oxidation of NADH by UQ, however, the PMF is determined by the competition between proton pumping by CI and that by the ATP synthase. If this is the case, and CI rates are limiting, one might expect to see a decrease in PMF as pumping stoichiometry decreases.
14) "Although the enzymatic activity of CI* and its H+-pumping-to electron-transfer ratio still remain to be confirmed, we venture that, by potentially pumping protons at a 2H+ :2e- ratio, CI* may play novel roles in the regulation of electron flow and energy conservation, adding flexibility to the mitochondrial electron transport chain of plants."
First, I object to the term "novel" in this context. How can one claim novelty when one does not know how all the others work?
15) "These observations also suggest that, throughout eukaryotic evolution, different bioenergetic strategies have demanded a trade-off between respiratory efficiency and dynamism." What observations are being referred to here?
It seems to be inferred, but never clearly stated, that CI* is a true, active form of CI, but with half of the proton pumping capacity.
The only place this is really clearly stated are in statements like "this suggests that there are significant steady-state amounts of this assembly intermediate in V. radiata mitochondria under these conditions and that CI* may be playing an independent physiological function."
Yet, the text also seems to imply that this form is an assembly intermediate. In any case, if this is an active intermediate, they should demonstrate this with at least a rudimentary activity assay.
16) The text starts to speculate about "dynamism", which is an interesting topic worthy of a review/hypothesis paper, but the data does not really provide very strong arguments in favor of such flexibility. If the authors want to make the case that their structure does indicate such flexibility, they will need to bolster their analysis.
On the other hand, there are some interesting observations in the data, especially the CA subunit, that could provide the bases of a more powerful discussion.
Reviewer #2:
An assembly intermediate of the mitochondrial complex I (CI*) was isolated from etiolated seedlings of mung bean and its structure was determined at 3.9 A resolution. The structure revealed some diversity from the complexes of bacteria, yeast and mammals including the presence of γCA subunits. Based on the structure, they propose that the CI* may be functional and provides the flexibility on respiratory electron transport.
This is the first report of the structure of plant complex I (an assembly intermediate) and provides important scientific message. The quality of science is high, and text is written clearly in general. I have several opinions to improve the manuscript.
Specific comments
1) The authors did not provide any evidence on the CI* to be functional. The idea of the regulation of the H+/e- ratio is attractive but the discussion sounds rather speculative and redundant (the last four paragraphs).
2) Plant mitochondria function in consuming excessive reducing power generated in chloroplasts and the proposed idea is interesting. But the CI* complexes were isolated from etiolated seedlings. Is it possible to detect the CI* complexes in green tissues?
3) They did not show the activity of the CI* complex. Is there any genetic evidence supporting that the complex I is functional in the absence of PD?
4) I am unsure why they focused on the assembly intermediate before clarifying the structure of the fully assembled complex I.
5) P5, L18. NUM2?
Reviewer #3:
The authors present a cryoEM structure at a nominal 3.9 Å resolution for an intermediary form of mitochondrial respiratory Complex I isolated from etiolated mung bean tissue. This is an exciting and timely manuscript and represents the first report of a (pseudo)atomic structure of Complex I (CI) from plants.
I am, in general, enthusiastic about this manuscript and consider it within the remit for publication within eLife. It is well written, and the description of the structure especially seems very sound. I do, however, have questions for the authors concerning the methodology, and in particular, the discussion.
i) The enzyme, as isolated from etiolated mung bean by detergent treatment, is lacking two major (proton-translocating) subunits from the membrane domain (the “PD” domain) and is described throughout the manuscript as CI*, as an intermediate, assembly form of Complex I. This is, in principle, acceptable, as this CI* intermediate has been documented by others elsewhere. However, it does raise the question whether this is an artefact of isolation (i.e. JBC 278: 43114 , 2003). Did the authors observe any intact CI in their preparation, or evidence for a solubilised NuoL/M subcomplex? It also remains to be seen if CI* is a peculiarity to etiolated tissue, and perhaps the authors can expand on this.
ii) The enzymatic activity of the CI* preparation was not determined. The authors suggest, that it functions as a 2H+/2e pump, by analogy with mutant forms of the fungal enzyme. This in itself is a reasonable suggestion, as this stoichiometry has also been observed with similar mutants of E. coli Complex I (Biochemistry 50: 3386, 2011), however, it is a pity that this was not investigated here as a discussion of the bioenergetics of CI* form a central part of the manuscript.
iii) The authors write that "Moreover, in photosynthetic tissue, conditions of intense light may lead to an over-production of NADH". I think this statement needs qualification for the benefit of non-expert readers. Obviously, the primary reductant produced by photosynthetic linear electron flow is NADPH and plants have many mechanisms to deal with sudden increases in light intensity (induction of NPQ, photosynthetic control at the level of cytochrome b6f etc.). Now, it is true that “excess” reducing equivalents are transferred from the chloroplast via the malate shunt, but this is a relatively low flux pathway, and so I assume that the authors are referring to NADH produced by the photorespiratory pathway via glycine decarboxylase activity. It is reasonable to assume that this may cause relatively large fluctuations in intramitochondrial [NADH]. (This also, obviously, will not be a factor in non-photosynthetic/etiolated tissue).
Related to point ii) above, the authors write "If the NADH/NAD+ ratio dropped in photosynthetic cells due to decreased light availability, this would lead to a decrease in free energy available to CI for H+ pumping. Full-length CI must pump 4 protons per 2 NADH electrons, it cannot operate in a graduated manner. Consequently, at a decreased NADH/NAD+ ratio, with lower free energy available, full-length CI would initially not be able to pump any protons at all. This would result in a reduction in the mitochondrial membrane potential, decreasing the proton-motive force and lowering ATP generation."
The thinking seems muddled here. Firstly, Complex I can be assumed to be operating in equilibrium (or near equilibrium) with the PMF, and so a “variable” H+/e pumping stoichiometry afforded by CI*/CI provides no benefit. The authors can demonstrate this for themselves through elementary thermodynamic principles if they seek reassurance (i.e. the equilibrium relationships (H+/2e)*PMF = 2ΔEh and also K = 10^((ΔEh-(PMF*(H+/2e))/30) etc) – in fact, I would encourage this as a supplementary figure. I doubt that there will be circumstances under which a supply of reducing equivalents to the plant mitochondrial respiratory chain will become limiting, particularly given the menagerie of electron donors that supply it. I can think of one instance in which 2H+/2e translocating CI* would provide a benefit over a 4H+/2e counterpart; it would be effectively irreversible at high (>200 mV) PMF, which may limit ROS production due to reverse ET (clearly the non-protonmotive NDA enzymes can be considered irreversible under all circumstances).
iv) The suggestion that the γCA subunit may be acting as a bicarbonate sensor by conformational change transmitted to the membrane domain seems very speculative.
v) Typo – Figure 1—figure supplement 3 – “accesory” in many places.
https://doi.org/10.7554/eLife.56664.sa1Author response
Reviewer #1:
[…]
1) "Whereas the exact γCA protein combinations are likely tissue- and development-stage-dependent, the role of the γCA domain in plant CI's function is unknown."
How does this relate to the CA thought to be involved in the CCM in cyanobacteria?
The carbonic anhydrase family is a highly diverse family of enzymes, composed of seven genetically distinct classes. While the seven classes possess carbonic anhydrase/CO2 hydratase activity, the structure and domain composition of the seven classes is diverse, and the sequence similarity across families is low. In α- and γ-CA classes, the metal ion in the active site is coordinated by three histidine residues. In contrast, the metal is coordinated by one histidine and two cysteine residues in the β-CA class. Moreover, α-CAs are active as monomers or dimers and β-CAs (which have a distinct structure compared to α-CAs) are active as tetramers. Finally, γ-CAs have a characteristic left-handed parallel β-sheet fold and are only active as trimers with active sites typically between subunits.
Genomic and phylogenetic studies have shown that bacteria may encode for α, β or γ CAs. In contrast, cyanobacteria have only been proven to encode for α- or β-CA in their various CO2-concentrating mechanisms (CCM). Moreover, searches of DNA databases of vascular plants and algae have not found clear homologues of cyanobacterial CCM genes. It is therefore believed that CCM genes, including those for cyanobacterial CAs, are almost exclusively restricted to cyanobacteria (and some proteobacteria). Nevertheless, cyanobacterial CA proteins involved in CCM do belong to the large CA protein family.
Therefore, the association seen between CA proteins and CI in cyanobacteria and plants is a case of convergent evolution, the physiology of which is incompletely understood. Given the recency of the first structures of cyanobacterial CI (Schuller et al., 2020) and plant CI (presented here), it still remains unclear whether the association of CI and CA in the two domains of life was driven to fulfil the same physiological function.
2) How do we know if this is really the CI* intermediate and not some sort of breakdown product?
The reviewer brings up an important point, which we should have emphasized in our manuscript. We present the arguments that support that CI* is an assembly intermediate and not a degradation product below. We have included these points in a new section in the Discussion of our revised manuscript (Discussion – Protein sample) and have added a new supplementary figure showing details of CI*’s purification and enzymatic activity (Figure 1—figure supplement 2).
We believe our sample is the CI* and not a degradation product for the following reasons:
– The 800-kDa CI subcomplex missing the peripheral-pumps domain (i.e. CI*) is a known subcomplex of plant CI, as determined by genetic and complexome profiling experiments of CI’s assembly pathway in Arabidopsis thaliana (Meyer et al., 2011, Schertk et al., 2012, Schimmeyer et al., 2016, Senkler et al., 2017, Ligas, et al., 2019). This intermediate is seen in non-etiolated seedlings and mature leaves in A. thaliana and N. sylvestris (Pineau et al., 2008). Moreover, the A. thaliana and N. sylvestris CI* showed NADH-dehydrogenase activity by the same in-gel activity assay we use in our preparation (Meyer et al., 2011, Pineau et al., 2008).
– The membrane-extraction conditions used in our prep (1% w:v digitonin, 4:1 g:g detergent:protein) are very gentle conditions for membrane protein extraction. These gentle conditions were chosen after optimization of different percentages, amounts and types of detergent to preserve protein:protein interactions in protein complexes.
– Immediately after extraction, we trap the detergent-extracted complexes in amphipatic polymers, which stabilize the complexes and protect them from degradation/dissociation (Breibeck and Rompel, 2019).
– These conditions are sufficiently gentle that they do not fully dissociate a large portion of lipid membrane that co-purifies with CI* and can be seen at low contour in our reconstructions (Figure 1—figure supplement 4E).
– Using the same detergent, digitonin, an A. thaliana complexome profiling study (Senkler, 2017) not only obtained full-length CI and CI*, but also full-length CI in a higher-order assembly with complex III, i.e. the “supercomplex” SC I+III2. The protein:protein interactions between complexes in a supercomplex more labile than intra-complex protein:protein interactions. Given that the more fragile CI:CIII2 interactions are maintained in 5% digitonin (Senkler et al., 2017), this argues that the presence of CI* – both in the Senkler and in our study – is not due to digitonin-induced dissociation of the distal domain, but rather that it is the true assembly intermediate.
– Published controlled-degradation experiments of plant CI in the presence of increasingly harsh detergents (increasingly strong disruption of protein:protein interactions within the complex) have shown that, similar to mammalian CI, plant CI detergent-induced dissociation occurs by dissociation of the full peripheral arm from the full matrix arm (Klodmann et al., 2010). In other words, CI detergent dissociation occurs at the “elbow” between CI’s arms and not by dislodging the distal pumps from the matrix-arm-associated-proximal pumps, which would need to occur if CI* were a degradation product.
– We have reproducibly obtained the CI* fraction, which retains in-gel and spectroscopic NADH-oxidase activity and chromatographic peak even after several days of manipulation. These data have been included in the new Figure 1—figure supplement 2.
3) I am confused by this:
“At each active site, two histidine residues are provided by one subunit and the third is provided by the adjacent subunit. However, in the plant CI γCA heterotrimer, the γCA2L subunit is lacking two of the three essential histidine residues (Ala-147 and Arg-152 in V. radiata) that would be necessary to form active sites at the interfaces with the CA1 and CA2 subunits. This renders two of the possible three catalytic sites non-functional (Figure 3A, Figure 3—figure supplement 2) …Our subunit assignments are consistent with this nonfunctional residue on CA1 being at the interface with CA2L, which would already be inactive due to the lack of histidine on CA2L at position."
It is unclear is the structure provides independent support that there are one of three active sites, or is their assignment based on this assumption?
In other words, there are hypotheses that the proteins represent one-of-three active CA or are completely inactive for CA activity. To what extent can the new structure actually test these possibilities? The test states that the Zn2+ was observed only one site. How strong is this evidence and how strongly does this support the one-in-three CA model? In this regard is it far more important what is inconsistent with the structure than what can be made to be consistent with it. In other words, if the structure can eliminate all but one of these hypotheses, it would be a big advance.
We thank the reviewer for pointing out these confusing statements. We have re-worded the revised document to clarify this issue.
The structure provides independent support that only one of the putative active sites could be active, as there is only one site where all three Zn-coordinating residues are present. The EM map shows clear density for the three histidine sidechains and the Zn2+ ion at the active site 1. Likewise, the density for the sidechains in active sites 2 and 3 is also clear. It is unambiguous that no histidine sidechain density is visible at active sites 2 and 3. Chemically, no Zn2+ ion could be coordinated by the residues in active sites 2 and 3. Consequently, no density for the Zn2+ ion is seen in those locations. This structural data provides independent confirmation of our subunit assignment from our sequence alignments. The structure eliminates all but one possible subunit arrangement to fit in the density seen and shows that only one of the three potential active sites allows for Zn2+ coordination and could potentially be actually enzymatically active.
4) "Therefore, carbonic anhydrase activity of the plant γCA domain on CI must be confirmed experimentally." Why not test for activity in the amazing new prep?
Once we are able to return to the lab, we are indeed planning on testing the potential carbonic anhydrase activity of CI* and CI, as well as the potential interplay between CO2 sensing and CI turnover. However, this is beyond the scope of this current manuscript.
5) "The lack of such a cap on V. radiata NU2M in CI* suggests that, although the Lys399 of NU2M is mostly surrounded by protein, the core hydrophilic axis may be in contact with lipid." Alternatively, some subunits may be missing or difficult to resolve in the prep.
Although it is theoretically possible that a subunit at the end of the proximal domain may have been lost, given our gentle extraction conditions and our exchange into amphipol, this is unlikely, as discussed in response to reviewer’s point 2. We agree with the reviewer that it is surprising to observe a lack of capping subunits on the end of the membrane arm, which is why we qualify the statement carefully with “suggests that” and “may be.” Nonetheless, the EM density does not show the presence of any subunits in this location – not even weak density that could be interpreted as a hard-to-resolve subunit – and we must report what we observe.
6) "This leaves open the possibility that the γCA domain does not interconvert CO2 and bicarbonate but acts only as a sensor of CO2 or bicarbonate concentration. In this scenario, conformational changes induced by bicarbonate binding could propagate into the membrane arm of CI and potentially regulate the catalytic turnover of CI. Such sensory and regulatory roles have also been proposed for other CI accessory subunits in other organisms."
Is there any evidence from the structure in support of are inconsistent with the proposed conformation propagation mechanism for CO2 sensing? Perhaps I am missing something, but the CA subunit seems to be rather loosely tethered, and it is difficult to see from the images provided how small conformational changes at the metal site would strongly affect the structure.
Also, given that this is a major point of the discussion, elaborate on why such signaling is interesting and what does the structure say about it. In particular, see points above on whether the structure confirms or not the proposed lack of catalytic activity.
Since the initial biochemical experiments on plant’s CA domain (Sunderhaus et al., 2006), it has been recognized that the CA domain is tightly linked to the membrane arm. Our structure shows that the interface between the CA domain and the membrane arm is extensive, with protein:protein interactions of the CA domain with several of the membrane-arm subunits (NU2M, NDUC2, P2, NDUX1). We were able to quantify the size and strength of the CA-membrane arm interface using the standard PDBePISA protein interface tool. We estimate that the interface between the CA domain and the membrane arm is ~3740 A2, with a solvation energy gain of ~210 kcal/mol. This gain in free energy is almost twice as large as the gain from the association of the CA domain itself, demonstrating that it is a tight interaction.
We have added these points to our description of the CA domain in the main text and have included a new supplementary table with the details of the interface size and strength (Table 5).
Given this large surface and the number of protein:protein interactions, it is feasible that a conformational change in the CA domain could propagate to the N or Q modules to affect NADH-ubiquinone oxidoreduction. However, our structure does not provide evidence for the catalytic activity of the CA domain (a long-standing hypothesis in the field) or for any potential influence of CA conformational changes on CI’s NADH-ubiquinone oxidoreduction. We have therefore removed the speculation about the potential carbonate sensing from our discussion.
7) The remaining differences between this and other structures are framed as important for assembly. Yet, how is it known that this is a real assembly intermediate (see point 2)?
CI* is a known assembly intermediate of plant CI, as discussed in response to point 2).
8) I am totally not convinced that there is any evidence from the structure one way or the other for a variable H+/e- transmission. Please convince me.
The nature of CI* is that it lacks two proton pumps (those forming the distal domain) compared to full-length CI, which contains four proton pumps and is known to pump protons at a 4:2 proton:electron ratio. Our structure shows the molecular arrangement of CI*, which contains all of the subunits and cofactors needed for the electron transfer reactions between NADH and CoQ but lacks two out of four proton-pumps. It is correct that our structure does not show that CI* pumps protons. However, the structure does strongly imply that if CI* were capable of pumping protons, it would be at a lower proton-to-electron ratio than full-length CI due to the lack of two proton pumping domains. This is a hypothesis at this stage, and its testing requires functional experiments. Nonetheless, genetic manipulations of CI in Y. lipolytica and E. coli have resulted in equivalent CI subcomplexes lacking the two distal H+-pumps. These CI subcomplexes are capable of NADH-CoQ oxidoreduction and proton-pumping at a 2:2 proton:electron ratio, rather than the 4:2 proton:electron ratio seen for full length CI (Drose et al., 2011, Steimle et al., 2011). Given these analogous studies, we conclude that given our structure it is reasonable to hypothesize that CI* could also be a functional subcomplex that pumps protons at a 2:2 ratio.
Moreover, in this revised manuscript, we add a new figure with biochemical data showing that CI* is capable of NADH-ubiquinone oxidoreduction (Figure 1—figure supplement 2). We have also re-written our discussion on CI*s potential lower proton:electron ratio to shorten the section and emphasize its speculative nature. As suggested by reviewer 3, we have moved the bioenergetic discussion into an Appendix, where we provide a more detailed and quantitative explanation of our rational for the possible role of CI* in plants, to address the reviewers’ bioenergetic concerns.
We intend to more thoroughly characterize the function of CI* and carefully quantify its proton:eletron ratio using a reconstituted system. However, given the competitive situation surrounding this manuscript and the pandemic-related limitations on wet-lab research, we believe it is better to develop this evidence in a follow-up manuscript focusing on a deep functional and enzymatic characterization of CI*.
9) "In heterotrophs that such as metazoans that obtain all their reduction equivalents from the breakdown of sugars, fats and proteins and that depend almost entirely on oxidative phosphorylation for energy production, there would be no advantage to being able to adjust CI's H+-pumping-to-electron-transfer ratio." This statement seems to be contradicted by the text that follows, which argues that adjusting H+/e- stoichiometry might be important for maintaining PMF (though I am not sure I agree with those).
The statement above presented the situation in heterotrophs in contrast with that in autotrophs. Nonetheless, for clarity and to focus more on our structural findings, we have removed these paragraphs from this section of the Discussion.
10) Regarding the paragraph from "In autotrophs such as plants, which generate their energy through a combination of respiration and photosynthesis, the bioenergetic processes necessarily have to be more dynamic to balance production of energy (ATP) and of reducing equivalents (NADH) throughout the cell and the organism. "
This is highly speculative, which in itself is OK, but I'm not sure what, if anything in the structure, supports this view.
We have re-written our Discussion to better align with the results from our structure and to remove unnecessary speculation. We have lengthened and made a separate section on the CA domain and shortened the speculative section on CI*’s potential roles. Moreover, we have included an Appendix to further explain our bioenergetic arguments for the possible roles of CI*.
11) The authors are correct that the energy balance of the cell will be strongly impacted by activating photosynthesis. However, there are several factors that need to be considered: 1) It is known for some time that respiration is down-regulated during photosynthesis, so that little if any ATP is produced by the mitochondria; 2) The critical factor is often the balance of energy in the forms of ATP/ADP+Pi and NAD(P)H/NADP+ rather than the total energy availability. In this case, a variable proton pumping stoichiometry might be useful, but; 3) There are already alternative mechanisms of shunting electrons to O2 without production of ATP, e.g. though the alternative oxidases or Water-water cycles.
We agree with the reviewer, although our understanding is that point 1 remains controversial, as recent experimental results and metabolic simulations suggest an important role for mitochondria in the supply of ATP to the cytosol under conditions in which the ATP consumption in photosynthetic CO2 fixation is sufficiently high (see references below). We have clarified these points in the Discussion and clarified the discussion in an Appendix.
12) "We posit that, in addition to bypassing H+-pumping altogether through the alternative complexes, it is conceivable that plants may find additional advantages in regulating the H+-pumping-to-electron-transfer ratio in the "canonical" electron transport chain. If the NADH/NAD+ ratio dropped in photosynthetic cells due to decreased light availability, this would lead to a decrease in free energy available to CI for H+ pumping.
13) "Full-length CI must pump 4 protons per 2 NADH electrons, it cannot operate in a graduated manner. Consequently, at a decreased NADH/NAD+ ratio, with lower free energy available, full-length CI would initially not be able to pump any protons at all. This would result in a reduction in the mitochondrial membrane potential, decreasing the proton-motive force and lowering ATP generation.
First, I think this point, if a full-length CI can pump lower Hstoichiometry of H+/e- is debatable (and being debated), so one cannot conclude that it cannot.
We have removed these paragraphs and clarified our bioenergetic arguments in an Appendix.
Second, the logic connecting the two sentences above seems to be incorrect. The smaller the coupling stoichiometry, the higher (not lower) the PMF that can be maintained at static head for a given NADH/NAD+ ratio.
CI is a reversible proton pumping-NADH:UQ oxidoreductase, so it will reach quasi-equilibrium (static head) when the free energy in the reaction
NADH + UQ + nH(n)+ ⇓ ◊ NAD+ UQH2 + nH+(p)
equals that in PMF.
The text then states: "however, if the plant ETC could find a way to keep pumping protons even at a decreased NADH/NAD+ ratio…" This would seem to be self-contradictory. One possibility is that the authors intend to include the total proton translocation of the system, and not just "pumping". In this case, PMF includes contributions from the chemical protons that are taken up during UQ reduction, which are likely to be deposited on the p-side of the membrane during UQH2 oxidation etc, e.g. at the Q-cycle.
The text continues, "however, if the plant ETC could find a way to keep pumping protons even at a decreased NADH/NAD+ ratio, this would allow for better maintenance of the proton motive force." The authors are apparently conflating near equilibrium and steady-state systems. With a smaller coupling stoichiometry, the number of protons pumped per NADH oxidized decreases, but the extent of PMF that can be generated at near equilibrium conditions increases. In a situation where the PMF is far out of equilibrium from the free energy in oxidation of NADH by UQ, however, the PMF is determined by the competition between proton pumping by CI and that by the ATP synthase. If this is the case, and CI rates are limiting, one might expect to see a decrease in PMF as pumping stoichiometry decreases.
We agree with the reviewer that this discussion was unclear. In this revised manuscript, we have made this discussion more succinct in the main text. In order to better express our arguments, we have included a more quantitative discussion in an Appendix.
We fully appreciate the distinctions pointed out by the reviewer between near-equilibrium and steady-state limitations. We address this in more detail in the Appendix to better support our claims about the possible roles of CI*.
The steady-state limit of the proton motive force (PMF) was determined in etiolated V. radiata mitochondria to be 200 mV (Moore and Bonner in 1981). Using reasonable assumptions (see Appendix for details), it is possible to determine that pumping 4 protons against a 200 mV PMF would require an [NAD+]/[NADH] ratio of 0.3 or less (i.e. the pool would have to be >70% reduced). Conversely, CI* pumping 2 protons would work in the forward direction against a 200 mV PMF up to a [NAD+]/[NADH] of >105. As the reviewer states, if CI is rate-limiting under these conditions, a decrease in pumping stoichiometry may result in a decrease in PMF, as ATP synthase outpaces the supply of protons. Nonetheless, the presence of a functional CI* under conditions of high [NAD+]/[NADH], even if rate-limiting, would still maintain the PMF to a greater extent than if all proton-pumping by CI ceased due to the loss of free energy of NADH oxidation. The degree to which this would be the case requires a steady-state model of the plant electron transport chain and would need good estimates of stoichiometries of the plant mitochondrial electron transport complexes (including CI*) and ATP synthase which currently (to the best of our knowledge) are not available and go beyond the scope of this paper.
Of course, 200 mV is an upper limit of steady-state proton pumping and likely the PMF in the plant cell is somewhat lower. Nonetheless, under any conditions in which full CI is operating near equilibrium, a two-proton-pumping CI* would be irreversible. This is more so the case for the non-proton pumping alternate NADH dehydrogenases (NDs). The point we failed to make clearly in the original manuscript was that if CI is operating near equilibrium and the [NAD+]/[NADH] ratio increases due to metabolic fluctuations (perhaps even due to changes in light, due to the input of the glycine dehydrogenase complex into mitochondrial NADH levels via the C2 photorespiration cycle), the PMF will decrease. The presence of CI* could be used to maintain the PMF at or closer to steady-state levels (with a 50% loss in thermodynamic efficiency) under a large range of [NAD+]/[NADH] ratios that favor reverse electron transport by full CI. Of course, many questions still remain, and this could lead to other issues of reactive oxygen species by CI via reverse electron transport. However, this problem already exists in plant mitochondria owing to the activity of NDs, which have an even stronger thermodynamic drive to push the energetics towards CI reverse electron transport. Many unknowns remain with respect to the co-expression of CI, the NDs – and potentially CI* – in mitochondrial membranes.
Moreover, the levels of CI* in our etiolated V. radiata sprouts appear higher than in reports from non-etiolated samples. We suspect this may be due to the seedlings’ major energy source (seed oils) being consumed, and cellular NADH levels dropping without energetic input from photosynthesis. This may favor a push towards less efficient but forward proton-pumping through the use of CI*. This will need to be tested experimentally and is beyond the scope of this manuscript.
14) "Although the enzymatic activity of CI* and its H+-pumping-to electron-transfer ratio still remain to be confirmed, we venture that, by potentially pumping protons at a 2H+ :2e- ratio, CI* may play novel roles in the regulation of electron flow and energy conservation, adding flexibility to the mitochondrial electron transport chain of plants."
First, I object to the term "novel" in this context. How can one claim novelty when one does not know how all the others work?
We take the reviewer’s point. We removed this sentence.
15) "These observations also suggest that, throughout eukaryotic evolution, different bioenergetic strategies have demanded a trade-off between respiratory efficiency and dynamism." What observations are being referred to here?
In order to improve the clarity of the manuscript and reduce speculation, we have removed this discussion from the manuscript.
It seems to be inferred, but never clearly stated, that CI* is a true, active form of CI, but with half of the proton pumping capacity.
The only place this is really clearly stated are in statements like "this suggests that there are significant steady-state amounts of this assembly intermediate in V. radiata mitochondria under these conditions and that CI* may be playing an independent physiological function."
Yet, the text also seems to imply that this form is an assembly intermediate. In any case, if this is an active intermediate, they should demonstrate this with at least a rudimentary activity assay.
The reviewer is correct that this is our hypothesis. We have added a section to the Discussion (Discussion – Protein sample) to explain the reasons why we believe our sample is CI* rather than a degradation product (see above). Our revised manuscript includes a new figure (Figure 1—figure supplement 2) showing CI*’s NADH oxidoreductase capabilities with in-gel activity and spectroscopic assays. This is in line with previous reports of CI* in-gel activity, as discussed above.
As pointed out by the reviewer, and as discussed above, we still have not tested the proton:electron ratio of CI*. We intend to carry out this quantification in the next stage of the project, where we will carry out a detailed biochemical and enzymatic characterization of CI* in solution and in reconstituted lipid vesicles.
16) The text starts to speculate about "dynamism", which is an interesting topic worthy of a review/hypothesis paper, but the data does not really provide very strong arguments in favor of such flexibility. If the authors want to make the case that their structure does indicate such flexibility, they will need to bolster their analysis.
On the other hand, there are some interesting observations in the data, especially the CA subunit, that could provide the bases of a more powerful discussion.
We take the reviewer’s point. We have re-written and shortened our main discussion to focus on the structure, minimizing the speculation on CI*’s potential roles, and expanding on the findings for the CA domain, as described above.
Reviewer #2:
[…]
Specific comments:
1) The authors did not provide any evidence on the CI* to be functional. The idea of the regulation of the H+/e- ratio is attractive but the discussion sounds rather speculative and redundant (the last four paragraphs).
Our revised manuscript includes a new figure (Figure 1—figure supplement 2) showing CI*’s NADH oxidoreductase capabilities with standard in-gel activity and spectroscopic assays. While we intend to carry out a deeper biochemical and enzymatic characterization of CI* both in solution and in reconstituted vesicles, we believe this is outside the scope of the current manuscript, especially given the practical limitations imposed by the pandemic at this time. We have re-written and shortened our discussion to minimize the speculation on CI*’s potentially altered proton pumping. We also provide supplementary discussion in Appendix with further details of the bioenergetic arguments.
2) Plant mitochondria function in consuming excessive reducing power generated in chloroplasts and the proposed idea is interesting. But the CI* complexes were isolated from etiolated seedlings. Is it possible to detect the CI* complexes in green tissues?
Although obtaining CI* from green tissues is beyond the scope of our manuscript, several other groups have detected CI* in green and non-etiolated tissues (both seedlings and mature leaves) of A. thaliana and N. sylvestris (Pineau et al., 2008, Meyer et al., 2011, Schertl et al., 2012, Schimmeyer et al., 2016, Senkler et al., 2017, Ligas, et al., 2019). We have clarified this in our revised main text and Discussion.
A possible reason for the apparently higher levels of CI* in our preparations of etiolated V. radiata compared to the previous reports with non-etiolated samples may be due to the seedlings’ consumption of their oil seeds during development and lack of other energy sources. We have expanded our discussion of these issues in an Appendix to ensure clarity and brevity in the main text.
3) They did not show the activity of the CI* complex. Is there any genetic evidence supporting that the complex I is functional in the absence of PD?
There is genetic evidence for this in other organisms. Knockout of Yarrowia lipolytica (yeast) mitochondrial CI subunit nb8m (nb8m∆) has been reported to prevent assembly of its CI distal-pump domain. The nb8m∆ CI subcomplex, which is analogous to CI*, is functional and pumps protons at half the H+:e- ratio of fully assembled CI (i.e. 2H+:2e- as opposed to the standard 4H+:2e- of intact CI). Additionally, deletion or truncation of E. coli’s NuoL (proximal-pump domain homologue of NU5M, which bridges the proximal and distal pumps) results in a functional CI that pumps protons at a 2H+:2e- ratio (Steimle et al., 2011). We have clarified this in the revised discussion and added the E. coli reference, as per reviewer 3’s recommendation. As described above, we have also included an additional figure showing CI* NADH-oxidoreductase activity (Figure 1—figure supplement 2).
The study of nad4 or nad5 (PD pumps) deletion mutants is challenging given that these genes are mitochondrially encoded. In plants, to our knowledge there is no direct genetic evidence supporting or refuting the hypothesis that CI* is a functional NADH-quinone dehydrogenase. However, it is known that complete absence of CI is not lethal to plants (Kuhn et al., 2015, Colas des Francs-Small and Small, 2014). Moreover, some “surrogate mutants” in nad4 or nad5 expression do not show macroscopic phenotypes (Colas des Francs-Small and Small, 2014 and references therein).
A paper studying the differences between Arabidopsis plants completely lacking CI due to mutation of the core flavoprotein subunit NDUFV1 (matrix arm) and those with reduced CI activity due to the mutation of the accessory subunit NDUFS4 (matrix arm) claims that CI subcomplexes are non-functional (Kuhn et al., 2015). This is based on their indirect evidence that ndufv1 and ndufs4 mutants both accumulate the same CI assembly subcomplexes. The argument is that, given that the CI subcomplexes in the severe ndufv1 mutants are not able to rescue at least some CI activity, this implies that CI subcomplexes in general (Kuhn et al., 2015), and CI* in particular (Schimmeyer et al., 2016), are non-functional. We find this unconvincing. Firstly, there is no characterization of the CI subcomplexes of either mutant beyond one unclear Western blot. Secondly, more importantly, the subcomplexes of CI (including CI*) of the ndufv1 mutants would also be lacking NDUFV1. Thus, these subcomplexes are not the same as wild-type subcomplexes and, given the lack of NDUFV1, would also not be expected to be functional in this mutant.
On the other hand, several studies have found that plant surrogate mutants with impaired expression of nad4 and nad5 (the PD pumps) accumulate CI* (Karpova and Newton, 1999, Pineau et al., 2008, Haili et al., 2013). A study of maize (Karpova and Newton, 1999) found that heteroplasmic nad4 mutants accumulate a “smaller, faster migrating, partially assembled” CI subcomplex that has NADH-dehydrogenase activity in the standard in-gel activity assay. Interestingly, the more severe the mutant phenotype, the higher the levels of the CI subcomplex. Given the timing of this study, however, it was not determined whether this subcomplex actually is CI*, as the CI* designation had yet to be defined. A study of Arabidopsis mutants of mtsf1 – an mRNA-binding protein essential for the stabilization and 3′-end processing of mitochondrial nad4 mRNA – found that these plants accumulate CI*, which is active in in-gel assays (Haili et al., 2013). Moreover, these mutants showed significantly higher oxygen consumption and lower carbon dioxide assimilation than wild-type. Although not understood, these observations were deemed “likely to stem from both a higher commitment of electrons to the AOX pathway and inefficient electron donation, thereby leading to an uncoupling effect”. Although speculative, these observations are also consistent with the activity of CI* pumping with a 2H+:2e- ratio.
4) I am unsure why they focused on the assembly intermediate before clarifying the structure of the fully assembled complex I.
As discussed in the manuscript, we believe the structure of CI* already offers a wealth of information about the differences between plant/mammalian/yeast/bacteria CI, its assembly pathway and additional potential biological implications.
5) P5, L18. NUM2?
Thank you. Corrected to NU2M.
Reviewer #3:
[…]
I am, in general, enthusiastic about this manuscript and consider it within the remit for publication within eLife. It is well written, and the description of the structure especially seems very sound. I do, however, have questions for the authors concerning the methodology, and in particular, the Discussion.
i) The enzyme, as isolated from etiolated mung bean by detergent treatment, is lacking two major (proton-translocating) subunits from the membrane domain (the “PD” domain) and is described throughout the manuscript as CI*, as an intermediate, assembly form of Complex I. This is, in principle, acceptable, as this CI* intermediate has been documented by others elsewhere. However, it does raise the question whether this is an artefact of isolation (i.e. JBC 278: 43114 , 2003). Did the authors observe any intact CI in their preparation, or evidence for a solubilised NuoL/M subcomplex? It also remains to be seen if CI* is a peculiarity to etiolated tissue, and perhaps the authors can expand on this.
We agree that this is an important point, which we should have emphasized in our manuscript. The arguments that support that CI* is an assembly intermediate and not a degradation product or a peculiarity of etiolated tissue are below. We have included these points in a new section in the Discussion of our revised manuscript (Discussion – Protein sample) and have added a new supplementary figure showing details of CI*’s purification and enzymatic activity (Figure 1—figure supplement 2).
Please also see response to Reviewer #1 point 2.
ii) The enzymatic activity of the CI* preparation was not determined. The authors suggest, that it functions as a 2H+/2e pump, by analogy with mutant forms of the fungal enzyme. This in itself is a reasonable suggestion, as this stoichiometry has also been observed with similar mutants of E. coli Complex I (Biochemistry 50: 3386, 2011), however, it is a pity that this was not investigated here as a discussion of the bioenergetics of CI* form a central part of the manuscript.
The revised manuscript includes a new figure (Figure 1—figure supplement 2) showing CI* that is enzymatically active in standard in-gel and spectroscopic activity assays. This is consistent with CI*’s previously reported in-gel activity, as described above. We thank the reviewer for the additional reference, which we now include in our revised Discussion. We have also added an Appendix with details of our bioenergetic arguments for the potential roles of CI*.
iii) The authors write that "Moreover, in photosynthetic tissue, conditions of intense light may lead to an over-production of NADH". I think this statement needs qualification for the benefit of non-expert readers. Obviously, the primary reductant produced by photosynthetic linear electron flow is NADPH and plants have many mechanisms to deal with sudden increases in light intensity (induction of NPQ, photosynthetic control at the level of cytochrome b6f etc.). Now, it is true that “excess” reducing equivalents are transferred from the chloroplast via the malate shunt, but this is a relatively low flux pathway, and so I assume that the authors are referring to NADH produced by the photorespiratory pathway via glycine decarboxylase activity. It is reasonable to assume that this may cause relatively large fluctuations in intramitochondrial [NADH]. (This also, obviously, will not be a factor in non-photosynthetic/etiolated tissue).
Related to point ii) above, the authors write "If the NADH/NAD+ ratio dropped in photosynthetic cells due to decreased light availability, this would lead to a decrease in free energy available to CI for H+ pumping. Full-length CI must pump 4 protons per 2 NADH electrons, it cannot operate in a graduated manner. Consequently, at a decreased NADH/NAD+ ratio, with lower free energy available, full-length CI would initially not be able to pump any protons at all. This would result in a reduction in the mitochondrial membrane potential, decreasing the proton-motive force and lowering ATP generation."
The thinking seems muddled here. Firstly, Complex I can be assumed to be operating in equilibrium (or near equilibrium) with the PMF, and so a “variable” H+/e pumping stoichiometry afforded by CI*/CI provides no benefit. The authors can demonstrate this for themselves through elementary thermodynamic principles if they seek reassurance (i.e. the equilibrium relationships (H+/2e)*PMF = 2ΔEh and also K = 10^((ΔEh-(PMF*(H+/2e))/30) etc) – in fact, I would encourage this as a supplementary figure. I doubt that there will be circumstances under which a supply of reducing equivalents to the plant mitochondrial respiratory chain will become limiting, particularly given the menagerie of electron donors that supply it. I can think of one instance in which 2H+/2e translocating CI* would provide a benefit over a 4H+/2e counterpart; it would be effectively irreversible at high (>200 mV) PMF, which may limit ROS production due to reverse ET (clearly the non-protonmotive NDA enzymes can be considered irreversible under all circumstances).
We thank the reviewer for this comment, which drove the creation of our Appendix, where we examine the bioenergetic arguments in more detail. Please see this Appendix and our response to reviewer #1’s points 12 and 13.
iv) The suggestion that the γCA subunit may be acting as a bicarbonate sensor by conformational change transmitted to the membrane domain seems very speculative.
Indeed, our structure does not provide evidence for any potential influence of CA conformational changes on CI’s NADH-ubiquinone oxidoreduction, and this currently remains speculative. However, our structure shows that the interface between the CA domain and the membrane arm is extensive, with protein:protein interactions of the CA domain with several of the membrane-arm subunits. Using the PDBePISA protein interface tool, we estimate that the interface between the CA domain and the membrane arm is ~3,740 A2. We have added the estimated size of the interface to the main text, together with an additional supplementary table (Table 5). Given this large surface and the number of protein:protein interactions, it is feasible that a conformational change in the CA domain could propagate to the N or Q modules to affect NADH-ubiquinone oxidoreduction.
However, given the speculative nature of the hypothesis and the reviewers’ comments, we have removed it from the Discussion.
v) Typo – Figure 1—figure supplement 3 – “accesory” in many places.
Thank you. Corrected.
References:
Colas des Francs-Small, C and Small, I. Surrogate mutants for studying mitochondrially encoded functions, Biochimie, Volume 100, 2014, Pages 234-242, https://doi.org/10.1016/j.biochi.2013.08.019.
Heddy Soufari, Camila Parrot, Lauriane Kuhn, Florent Waltz, Yaser Hashem
bioRxiv 2020.02.21.959148; doi: https://doi.org/10.1101/2020.02.21.959148
Iñiguez C, Capó-Bauçà S, Niinemets Ü, Stoll H, Aguiló-Nicolau P, Galmés J. Evolutionary trends in RuBisCO kinetics and their co-evolution with CO2 concentrating mechanisms. Plant J. 2020;101(4):897‐918. doi:10.1111/tpj.14643
Klanchui et al. (2017). Exploring Components of the CO2-Concentrating Mechanism in Alkaliphilic Cyanobacteria Through Genome-Based Analysis. Comput Struct Biotechnol J. 15:340-350. doi: 10.1016/j.csbj.2017.05.001
Kühn K, Obata T, Feher K, Bock R, Fernie AR, Meyer EH. Complete Mitochondrial Complex I Deficiency Induces an Up-Regulation of Respiratory Fluxes That Is Abolished by Traces of Functional Complex I. Plant Physiol. 2015;168(4):1537‐1549. doi:10.1104/pp.15.00589
Long et al. (2016). Cyanobacterial CO2-concentrating mechanism components: function and prospects for plant metabolic engineering. Curr Opin Plant Biol.31:1-8. doi: 10.1016/j.pbi.2016.03.002.
Supuran and Capasso (2017). An Overview of the Bacterial Carbonic Anhydrases. Metabolites 7. doi: 10.3390/metabo7040056
Schuller, J. M. et al. (2020). Redox-coupled proton pumping drives carbon concentration in the photosynthetic complex I. Nat. Commun. 11, 494–7.
https://doi.org/10.7554/eLife.56664.sa2Article and author information
Author details
Funding
No external funding was received for this work.
Acknowledgements
We are grateful to K Abe, MG Zaragoza, C Goodwin, R Murguia and C Bower for help with V. radiata growth and mitochondrial isolations. Data was collected at the UCSF WM Keck Foundation Advanced Microscopy Laboratory, accessed through the Bay Area Cryo-EM Consortium BACEM, with the assistance of D Bulkley and Z Yu. We are grateful to J Al-Bassam, E Baldwin, JC Lagarias and J Callis for helpful input on the manuscript and to W Broadly of the UC Davis High-Performance Cluster for technical assistance. MM acknowledges funding from the UC Davis POP Program.
Senior Editor
- Cynthia Wolberger, Johns Hopkins University School of Medicine, United States
Reviewing Editor
- David M Kramer, Michigan State University, United States
Reviewers
- Toshiharu Shikanai
- Nicholas Fisher, Michigan State University, United States
Publication history
- Received: March 5, 2020
- Accepted: July 22, 2020
- Version of Record published: August 25, 2020 (version 1)
Copyright
© 2020, Maldonado 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
-
- 3,358
- Page views
-
- 344
- Downloads
-
- 33
- Citations
Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.
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
A complex interplay between structure, conformational dynamics and pharmacology defines distant regulation of G protein-coupled receptors.
-
- Biochemistry and Chemical Biology
The glucagon-like peptide-1 receptor (GLP1R) is a broadly expressed target of peptide hormones with essential roles in energy and glucose homeostasis, as well as of the blockbuster weight-loss drugs semaglutide and liraglutide. Despite its large clinical relevance, tools to investigate the precise activation dynamics of this receptor with high spatiotemporal resolution are limited. Here, we introduce a novel genetically encoded sensor based on the engineering of a circularly permuted green fluorescent protein into the human GLP1R, named GLPLight1. We demonstrate that fluorescence signal from GLPLight1 accurately reports the expected receptor conformational activation in response to pharmacological ligands with high sensitivity (max ΔF/F0=528%) and temporal resolution (τON = 4.7 s). We further demonstrated that GLPLight1 shows comparable responses to glucagon-like peptide-1 (GLP-1) derivatives as observed for the native receptor. Using GLPLight1, we established an all-optical assay to characterize a novel photocaged GLP-1 derivative (photo-GLP1) and to demonstrate optical control of GLP1R activation. Thus, the new all-optical toolkit introduced here enhances our ability to study GLP1R activation with high spatiotemporal resolution.