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).
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).
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).
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).
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).
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).
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).
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
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).
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).
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
Article 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.
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,657
- views
-
- 389
- downloads
-
- 43
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
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
- Cell Biology
Activation of the Wnt/β-catenin pathway crucially depends on the polymerization of dishevelled 2 (DVL2) into biomolecular condensates. However, given the low affinity of known DVL2 self-interaction sites and its low cellular concentration, it is unclear how polymers can form. Here, we detect oligomeric DVL2 complexes at endogenous protein levels in human cell lines, using a biochemical ultracentrifugation assay. We identify a low-complexity region (LCR4) in the C-terminus whose deletion and fusion decreased and increased the complexes, respectively. Notably, LCR4-induced complexes correlated with the formation of microscopically visible multimeric condensates. Adjacent to LCR4, we mapped a conserved domain (CD2) promoting condensates only. Molecularly, LCR4 and CD2 mediated DVL2 self-interaction via aggregating residues and phenylalanine stickers, respectively. Point mutations inactivating these interaction sites impaired Wnt pathway activation by DVL2. Our study discovers DVL2 complexes with functional importance for Wnt/β-catenin signaling. Moreover, we provide evidence that DVL2 condensates form in two steps by pre-oligomerization via high-affinity interaction sites, such as LCR4, and subsequent condensation via low-affinity interaction sites, such as CD2.
-
- Biochemistry and Chemical Biology
The development of proteolysis targeting chimeras (PROTACs), which induce the degradation of target proteins by bringing them into proximity with cellular E3 ubiquitin ligases, has revolutionized drug development. While the human genome encodes more than 600 different E3 ligases, current PROTACs use only a handful of them, drastically limiting their full potential. Furthermore, many PROTAC development campaigns fail because the selected E3 ligase candidates are unable to induce degradation of the particular target of interest. As more and more ligands for novel E3 ligases are discovered, the chemical effort to identify the best E3 ligase for a given target is exploding. Therefore, a genetic system to identify degradation-causing E3 ligases and suitable target/E3 ligase pairs is urgently needed. Here, we used the well-established dimerization of the FKBP12 protein and FRB domain by rapamycin to bring the target protein WDR5 into proximity with candidate E3 ligases. Strikingly, this rapamycin-induced proximity assay (RiPA) revealed that VHL, but not Cereblon, is able to induce WDR5 degradation - a finding previously made by PROTACs, demonstrating its predictive power. By optimizing the steric arrangement of all components and fusing the target protein with a minimal luciferase, RiPA can identify the ideal E3 for any target protein of interest in living cells, significantly reducing and focusing the chemical effort in the early stages of PROTAC development.