A three-dimensional topology of complex I inferred from evolutionary correlations
© Kensche et al.; licensee BioMed Central Ltd. 2012
Received: 13 February 2012
Accepted: 28 June 2012
Published: 3 August 2012
The quaternary structure of eukaryotic NADH:ubiquinone oxidoreductase (complex I), the largest complex of the oxidative phosphorylation, is still mostly unresolved. Furthermore, it is unknown where transiently bound assembly factors interact with complex I. We therefore asked whether the evolution of complex I contains information about its 3D topology and the binding positions of its assembly factors. We approached these questions by correlating the evolutionary rates of eukaryotic complex I subunits using the mirror-tree method and mapping the results into a 3D representation by multidimensional scaling.
More than 60% of the evolutionary correlation among the conserved seven subunits of the complex I matrix arm can be explained by the physical distance between the subunits. The three-dimensional evolutionary model of the eukaryotic conserved matrix arm has a striking similarity to the matrix arm quaternary structure in the bacterium Thermus thermophilus (rmsd=19 Å) and supports the previous finding that in eukaryotes the N-module is turned relative to the Q-module when compared to bacteria. By contrast, the evolutionary rates contained little information about the structure of the membrane arm. A large evolutionary model of 45 subunits and assembly factors allows to predict subunit positions and interactions (rmsd = 52.6 Å). The model supports an interaction of NDUFAF3, C8orf38 and C2orf56 during the assembly of the proximal matrix arm and the membrane arm. The model further suggests a tight relationship between the assembly factor NUBPL and NDUFA2, which both have been linked to iron-sulfur cluster assembly, as well as between NDUFA12 and its paralog, the assembly factor NDUFAF2.
The physical distance between subunits of complex I is a major correlate of the rate of protein evolution in the complex I matrix arm and is sufficient to infer parts of the complex’s structure with high accuracy. The resulting evolutionary model predicts the positions of a number of subunits and assembly factors.
KeywordsEukaryotic complex I Quaternary topology Assembly Mirror-tree method Co-evolution
NADH:ubiquinone oxidoreductase (complex I) is with about 1000 kDa [1, 2] the largest of the five complexes of the oxidative phosphorylation (OXPHOS) and a major contributor to the proton motive force that drives the ATP production by ATP-synthase . Complex I has an L-shape with a hydrophilic matrix arm that protrudes into the cytoplasm in bacteria or the mitochondrial matrix in eukaryotes and a hydrophobic membrane arm. The canonical “core” of complex I consists of 14 subunits that originate from three pre-existing evolutionary modules . The N-module at the distal end of the matrix arm contains flavin-mononucleotide (FMN) that accepts electrons from a donor, usually NADH. The electrons are transported through a chain of iron-sulfur (FeS) clusters along the matrix arm towards the joint of the two arms at the membrane. This membrane-proximal part of the matrix arm represents the Q-module in which the electrons are transferred to ubiquinone (Q). The energy freed by the electron-transfer is transmitted along the P-module (NADH1-6/4L) that uses the energy to pump protons across the membrane [5–7].
In diverse taxa, the canonical core of complex I has been extended by further subunits. For instance, complex I in Thermus thermophilus contains an additional subunit located at the interface of the N- and Q-modules  and a recent analysis of complex I in the α-proteobacterium Paracoccus denitrificans identified three additional subunits . Eukaryotes obtained complex I with the endosymbiotic uptake of an α-proteobacterium that gave rise to the mitochondria. Following the endosymbiosis, the mitochondrial genome was reduced and the genes encoding matrix arm subunits of complex I were transferred to the nucleus. Additionally, complex I was extended to up to 45 subunits by so-called “accessory” or “supernumerary” subunits [1, 10]. This set of permanent subunits is further extended by a number of assembly factors absent from the mature complex [11–19].
Up to now, the structures of the complete complex in the eubacterium Thermus thermophilus and the eukaryote Yarrowia lipolytica have been published. However, the latter structure is of a too low resolution to allow identification of the positions of the supernumerary subunits. Approximate subunit positions within the eukaryotic complex are hinted at by various types of experiments, mostly from sub-complexes observed by fractionation or as assembly intermediates (e.g. [1, 20]). For instance, the application of chaotropic detergents to the bovine complex produces the three sub-complexes Iα, Iλ, and Iγ. Because these sub-complexes are large, they provide only rough information about subunit positions. For instance, Iα represents an extended Iλ sub-complex and the additional subunits could in principle be located anywhere on the surface of the Iλ sub-complex. Only limited data are available from yeast-two-hybrid [12, 21], co-immunoprecipitation , or cross-linking  experiments. The identification of the positions of the assembly factors is hampered by the temporariness of the assembly intermediates and our incomplete understanding of the assembly process.
The increasing number of genome sequences allows making predictions of physical interactions by evolutionary correlation methods, including the co-occurrence of genes or phylogenetic profiling [23–26], the mirror-tree approach , and residue correlation [28–31], which have successfully identified new complex I subunits and assembly factors [21, 25] (reviewed in ) and predicted relations between the five OXPHOS complexes . Of these methods, residue correlation is based on the most direct evidence of physical interaction, namely the compensatory mutations at sites of interacting residues to maintain the structure of a protein or complex. By contrast, the mirror-tree method detects co-evolution more indirectly by correlating sequence similarity matrices between orthologous groups . The similarities between protein sequences depend both on species divergence times and on rates of evolution. By removing the similarity due to the species divergence times [34, 35], one obtains similarities that are more related to evolutionary rates. A high correlation in evolutionary rates between protein families can be evidence of a direct physical relation between proteins. For instance, if there is selection to maintain the interaction of two proteins then disrupting mutations have to be compensated for at the rate that they occur. Therefore, to maintain the interaction, an increased rate of change in one protein needs to be compensated for at a similar rate in the other protein. Note that the pairwise correlations in evolutionary rates between proteins can also be due to indirect interactions .
Here we ask whether we can use evolutionary rate correlations to predict the three-dimensional (3D) conformation of complex I. To this aim, we analyzed the evolutionary rate correlation between 38 subunits and 7 assembly factors of human complex I using the mirror-tree method  and find that subunits that are known to be physically close in complex I tend to show a higher degree of correlation in evolutionary rates than those that are physically distant. In the conserved core of the matrix arm, this correlation is strong enough to construct a 3D model with striking similarity to the bacterial reference structure. In a second evolutionary model that includes the 14 canonical core subunits, the membrane and matrix arms appear as clearly separated groups. Finally, we calculated a third evolutionary model including 38 subunits and seven assembly factors. This last model retains some features of the physical structure, including the separation of the matrix and membrane arms and a proximodistal axis in the matrix arm. We discuss the positions of the seven assembly factors in this model and make specific predictions about the association of some assembly factors with each other and with the permanent subunits.
Correlation of physical and co-evolutionary distances in the conserved core
Because our study aims at complex I of human we selected 38 subunits and seven assembly factors of the human complex that have a sufficient number of orthologs for the application of the mirror-tree method (Additional file 1: Table S1). Shortly, we collected homologs of the 45 proteins by querying the nr database. All orthologous sequence sets were aligned and highly variable alignment columns were filtered out using BMGE . Note that NADH3, 4L and 6 could not unambiguously be located in our reference – the structure of the complete complex I of T. thermophilus (PDB:3M9S). However, these subunits are known to be direct neighbors and we decided to treat them as a single unit, termed NADH34L6. We calculated maximum-likelihood trees from the alignments  and obtained a distance matrix for each protein family from which we removed the common signal of the phylogeny . The phylogeny-corrected matrices were correlated and the resulting correlation matrix was transformed into the distance matrix (see Methods). For the purpose of this article we call these distances “co-evolutionary” but note that the signal measured by the mirror-tree method is also determined by expression and general functional relatedness [39, 40].
After the prediction of the topology of the seven matrix arm core subunits, we predicted the topology of the complete core of 14 matrix- and membrane arm subunits that are conserved among bacteria and eukaryotes. Again, the co-evolutionary distances of the subunits can be well embedded in three dimensions (P 3 =0.84) and result in a configuration with separate membrane and matrix arms (Figure 2c). The rmsd of this evolutionary model and bacterial structure is 47.6 Å, which corresponds to about 25% of the length of the membrane arm . Also in the complete-core model, the proximodistal axis of the matrix arm is recovered (r2=0.96, p=9.3×10-5) although the rmsd of the matrix arm core subunits to the bacterial structure is lower (34.0 Å) than that for the configuration that only contained matrix arm subunits (18.7 Å). As expected from the correlation between the pairwise physical distances and their evolutionary correlation, the accuracy of the positioning in the membrane arm is poor (rmsd=61.8 Å). A closer inspection of the evolutionary model suggests that the membrane arm subunits are located as a cluster more sideways of the matrix arm than in the physical structure. In the superimposition of the two structures, this manifests in a tendency of the Q-module away from the membrane and a slight tendency of the distal N-module towards the membrane.
Eukaryotic complex I and assembly factors
Our results show that the evolutionary rates of complex I subunits contain a significant amount of information about the complex’s quaternary structure. For the matrix arm we found that about 61% of the correlation in evolutionary rate could be explained by the distances of the subunit centers. This is even more striking if we consider that the evolutionary model was derived from eukaryotic sequences and thus should reflect the matrix arm structure in eukaryotes, while our reference structure is from a bacterium. Indeed, the evolutionary 3D model revealed a twist between the N-module and the Q-module when compared to the bacterial structure, a finding that is supported by experimental data (, personal communication).
In the two models that include the membrane arm, mitochondria-encoded subunits were predicted to be separate from nucleus-encoded subunits, which is in line with previous results. The independent variation of evolutionary rates in the nuclear and mitochondrial genomes  may have contributed to the isolation of the mitochondria encoded subunits in our models. Nevertheless, we stress that the position of the membrane arm core subunits, specifically at the proximal end of the matrix arm in both models, indicates a signal of the physical structure in the evolutionary correlation data. Furthermore, the strict separation of the nucleus-encoded Iβ subunits and their mitochondria-encoded counterparts NADH4 and 5 may be explained by other factors, as these two groups also behave differently in experiments [1, 52, 53]. Interestingly, despite its nuclear encoding, the membrane-integral subunit NDUFA1 of the Iα sub-complex [1, 2, 54] is positioned close to the membrane arm core, in particular close to NADH2 (Figures 3a and 3b). In T. thermophilus, NADH2 is located between the two subunits NADH1 and NADH4  both of which are known physical interactors of NDUFA1 . A direct physical interaction of NDUFA1 with NADH2 is therefore likely.
The evolutionary correlation failed to identify the correct topology of the membrane arm core. A number of biological reasons could explain such a lack of signal. First, long-range structural constraints [5, 6] may interfere with the distance-dependent structural constraints that are necessary for a distance-dependent strength of evolutionary rate correlation. Second, the formation of OXPHOS super-complexes with complex I dimers may result in correlations between distant subunits. Indeed, despite their positions at opposite ends of the membrane arm, NADH1 and 5 show a high correlation in evolutionary rates with each other and with subunit CYTb of complex III  consistent with their proximity in OXPHOS complexes organized into respiratory strings . Third, the lack of correlation with physical distance may result from non-adaptive variation in the mitochondria-encoded genes caused by variable and, at least in some eukaryotic taxa, heterogeneous mutation-pressure . Indeed, in a number of animal taxa changes in gene order or mutation-pressure led to non-adaptive changes in mitochondrial genes [57–59]. The mitochondrial genomes of some taxa in our study, such as plants, are clearly different from those in animals (reviewed in ) and their genes are likely under different mutation-pressures . Fourth, the embedding of the membrane proteins in two dimensions might reduce the evolutionary constraints to maintain interactions in comparison to proteins that are embedded in three dimensions.
The integration of multiple proteins in a single model assumes that the interactions are permanent and non-competitive. This is clearly not the case for the model of 45 proteins because it includes assembly factors. This model can therefore not exactly represent a physical structure. According to current models, complex I assembles from independent subcomplexes . Of the assembly factors required for this process and included in our study, only NDUFAF1 (AF1) is required for the assembly of the distal membrane arm sub-complex [13, 63, 64]. In our model, AF1 is located close to the matrix arm, which supports an indirect rather than a direct involvement of AF1 in membrane arm assembly . The distal membrane arm further combines with a pre-formed membrane-anchored proximal matrix/membrane arm that contains the subunits NDUFS2 (S2) and NADH1 (1) and possibly NDUFS3 (S3) and NDUFS7 (S7) [62, 64] and whose assembly involves NDUFAF3 (AF3) and possibly C8orf38 (O38) [17, 66]. Although the membrane-association of AF3 and O38 is not reflected in our data, they form a tightly co-evolving triple with C2orf56 (O56), which is known to bind the proximal matrix arm subunit S2 . The high correlation in evolutionary rates between AF3, O38, and O56 suggest strong selective constraints on their cooperation during the assembly of the proximal matrix/membrane arm sub-complex. The fourth assembly factor that has been experimentally linked to the proximal membrane arm, C20orf7 (O7) [18, 64], is indeed placed close to the proximal matrix arm subunits S2, S3 (Figure 3b, right bottom), and the proximal membrane arm subunits A3, A6, and A8 [2, 47, 48].
After the joining of the two membrane arm intermediates, the proximal matrix arm is further extended. This step involves the NUBPL-mediated assembly of at least one FeS-cluster into the distal matrix arm [11, 67]. In the evolutionary configuration the assembly factor NUBPL is positioned side by side with the permanent subunit NDUFA2 (A2; Figure 3b, right top). Like NUBPL, A2 is associated with the distal matrix arm . The highly conserved A2 subunit is structurally similar to thioredoxin-like proteins with a loop-region of probably variable conformation that contains two cysteines in human (C24 and C58) . These cysteines can form a revertible disulfide bridge with an in-vitro redox-potential in the range of the large majority of isopotential FeS-clusters of complex I [69, 70]. Although the cysteines are not fully conserved, occasionally FeS-clusters are bound by serine, histidine, or aspartate . Indeed, the human serine 30 in NDUFA2 is a good candidate for FeS-cluster binding because it is perfectly conserved in all species, with the notable exceptions of Trypanosoma and Leishmania, in which it is substituted by cysteine. Together these observations and the very strong evolutionary rate correlation of A2 and NUBPL support an involvement of A2 in complex I associated FeS-cluster assembly or maintenance. The peripheral position of A2 and NUBPL in the model could be a consequence of other strong evolutionary constraints not directly related to complex I.
Also NDUFAF2 (AF2, B17.2L) has been linked to the assembly of the distal matrix arm [14, 64]. Interestingly, the evolutionary data position AF2 directly besides its paralog NDUFA12 (A12, B17.2) [10, 14]. Like AF2, A12 is known to be associated to the distal-matrix arm to which it is directly recruited from the mitochondrial matrix . The correlation in evolutionary rates and the independent co-loss in multiple complex I lacking taxa  support an evolutionarily conserved functional relationship of AF2 and A12. It is tempting to speculate that AF2 temporarily binds at the binding site of A12, e.g. to stabilize the local structural context, and is later substituted by its paralog. Such close positioning and physical interaction of homologous proteins within the same protein complex is one of the prevailing trends in the “fate” of duplicated proteins in complexes . Complex I appears to add another twist to this pattern in the sense that the predicted interaction is only temporary.
The rate of protein evolution is influenced by diverse factors , in particular expression and general functional relatedness [39, 40, 74]. It is therefore even more remarkable that we found physical distance to be the major determinant of the evolutionary rate correlation for the complex I matrix arm. However, this result does not apply to the whole complex. Thus, to establish whether the mirror-tree/MDS combination is a good general method to predict quaternary structures, other complexes need to be analyzed. Furthermore, instead of using the mirror-tree method one could use residue correlation to measure the co-evolution of subunits more directly. Residue correlation has been used to predict contact interfaces for protein pairs [30, 75] and to investigate a rotation-symmetric homo-multimeric complex . A simple implementation would be to integrate pairwise residue correlations  or correlations that account for indirect correlations [30, 31, 77] or phylogenetic dependency [76, 78] by in-silico two-hybrid  into subunit distances and map these into three dimensions by multidimensional scaling.
The correlations of evolutionary rates between subunits of the eukaryotic complex I contain detailed information about the structural arrangement of the matrix arm subunits. This allowed us to make specific predictions about the positions of supernumerary subunits and assembly factors of the matrix arm, which may guide further experimental investigations. Multidimensional scaling could not reconstruct the structure of the membrane arm core. A future analysis will have to investigate what may cause this lack a spatial signal along the membrane arm and thus clarify in particular the relevance of conformational dynamics and super-complex arrangement into a respiratory string for the sequence evolution of complex I.
We included 38 permanent subunits and seven assembly factors of human complex I for which a sufficient number of sequences were available. We collected homologous sequences from the nr database  using PSI-BLAST (default parameters). Multiple queries from different species were used whenever PSI-BLAST failed to find known homologs [see Additional file 1. For A6, B9, A12, and AF2, orthologous groups were manually identified in neighbor-joining trees constructed with identity matrices and correcting for multiple substitutions. Species overlap between the partitions was used to divide the trees into separate orthologous groups. The remaining subunits were treated by a different protocol. First, to ensure a separation of the paralogs NADH2, 4, and 5, we built a set of trusted orthologs of NADH2, 4, and 5 from those sequences that had the best bi-directional hit with the human query using PSI-BLAST. From these seed sequence sets we computed three HMM profiles and sorted the remaining homologs into the orthologous group to which they showed the best profile-alignment . For all sequence sets we selected as single ortholog per species the sequence with the highest NEEDLE score in a pairwise alignment to the human query  (default parameters) and/or manual selection based on multiple alignments (MAFFT , CLUSTALW [84, 85], HMMER , HHSEARCH ). The kinetoplastida were excluded from our analysis due to their high level of sequence divergence. To gain high quality alignments, we aligned all sequence sets with CLUSTALW and manually fixed misalignments. The manually curated alignments are provided in Additional file 3. Next, we filtered alignment columns with BMGE  (−m BLOSUM30 -g 0.50 -b 4), removed sequences that had more than 33% gaps, and restricted the alignments to those species for which we found at least eight subunits of the complex. Of the 43 alignments, 39 had more than 75 sequences and there was no alignment with less than 44 sequences. Finally, we calculated phylogenetic trees using RAXML  (Version 7.2.6, PROTGAMMAMTREV for NADH1/2/3/4/4L/5/6, otherwise PROTGAMMAJTT; 4 rate categories) [see Additional file 4. A single tree was calculated for the concatenated alignment of NADH3, 4L, and 6.
with the row vector p T . We derived the reference vector p directly from the subunits’ distance vectors, as suggested by Kann et al.. Specifically, the reference distance between a pair of species was calculated as the average of the distances between these species in the trees of the complex I subunits. Note that this choice of reference as an average of the analyzed vectors also removes the specific pattern of co-variation in evolutionary rates that reflects selective pressure on the complex as a whole. It thus focuses the results on distances between the subunits rather than their distances to unrelated proteins. Finally, the corrected distance vectors were correlated with each other by Spearman rank correlation to yield the subunits' co-evolutionary similarity . We required that the species pairs were present in at least five of our 43 distance vectors. Species pairs occurring in fewer than five vectors were ignored in the correlations. Note that our choice of the set of subunits included the requirement that there are at least 15 species in all pairs of alignments. Only for 17 out of 903 subunit pairs, the correlation values were based on less than 30 species. The mirror-tree method has the advantages of being easily implemented and it requires low computational resources, even with correction for the basic correlation due to the shared phylogeny.
Multidimensional scaling (cMDS)
where λ i is the i-th largest eigenvalue of B Δ . Note that the relation between the co-evolutionary dissimilarity and the distance in the cMDS configuration (Shepard diagram, Additional file 1) indicates that the 3D configuration reflects the raw co-evolutionary distances over its whole range.
Superimposition of configurations
where |t i - c i | is the distance between the bacterial and predicted center of the i- th subunit.
We thank Richard Notebaart and Fiona Nielsen for valuable discussions. PRK was funded by the BioRange program of the Netherlands Bioinformatics Centre (NBIC) supported by the Netherlands Genomics Initiative (NGI) and the European Union’s Sixth Framework Program EPISTEM (CT-2005-019067). ID was funded by the Portuguese Foundation for Science and Technology (SFRH/32966/2006) and by Bolsas Rui Tavares 2010.
- Carroll J, Fearnley IM, Shannon RJ, Hirst J, Walker JE: Analysis of the subunit composition of complex I from bovine heart mitochondria. Mol Cell Proteomics 2003, 2(2):117–126. 10.1074/mcp.M300014-MCP200PubMedGoogle Scholar
- Hirst J, Carroll J, Fearnley IM, Shannon RJ, Walker JE: The nuclear encoded subunits of complex I from bovine heart mitochondria. Biochim Biophys Acta 2003, 1604(3):135–150. 10.1016/S0005-2728(03)00059-8PubMedGoogle Scholar
- Garcia-Vallve S: Contribution of each complex of the mitochondrial respiratory chain in the generation of the proton-motive force. Biochem Mol Biol Edu 2004, 32(1):17–19. 10.1002/bmb.2004.494032010308Google Scholar
- Friedrich T, Abelmann A, Brors B, Guénebaut V, Kintscher L, Leonard K, Rasmussen T, Scheide D, Schlitt A, Schulte U, et al.: Redox components and structure of the respiratory NADH:ubiquinone oxidoreductase (complex I). Biochim Biophys Acta 1998, 1365(1–2):215–219.PubMedGoogle Scholar
- Efremov RG, Baradaran R, Sazanov LA: The architecture of respiratory complex I. Nature 2010, 465(7297):441–445. 10.1038/nature09066PubMedGoogle Scholar
- Hunte C, Zickermann V, Brandt U: Functional modules and structural basis of conformational coupling in mitochondrial complex I. Science 2010, 329(5990):448–451. 10.1126/science.1191046PubMedGoogle Scholar
- Efremov RG, Sazanov LA: Structure of the membrane domain of respiratory complex I. Nature 2011, 476(7361):414–420. 10.1038/nature10330PubMedGoogle Scholar
- Sazanov LA, Hinchliffe P: Structure of the hydrophilic domain of respiratory complex I from Thermus thermophilus. Science 2006, 311(5766):1430–1436. 10.1126/science.1123809PubMedGoogle Scholar
- Yip C-, Harbour ME, Jayawardena K, Fearnley IM, Sazanov LA: Evolution of respiratory complex I: 'supernumerary' subunits are present in the α-proteobacterial enzyme. J Biol Chem 2010, 286(7):5023–5033.PubMed CentralPubMedGoogle Scholar
- Gabaldón T, Rainey D, Huynen MA: Tracing the evolution of a large protein complex in the eukaryotes, NADH:ubiquinone oxidoreductase (Complex I). J Mol Biol 2005, 348(4):857–870. 10.1016/j.jmb.2005.02.067PubMedGoogle Scholar
- Bych K, Kerscher S, Netz DJA, Pierik AJ, Zwicker K, Huynen MA, Lill R, Brandt U, Balk J: The iron-sulphur protein Ind1 is required for effective complex I assembly. EMBO J 2008, 27(12):1736–1746. 10.1038/emboj.2008.98PubMed CentralPubMedGoogle Scholar
- Carilla-Latorre S, Gallardo ME, Annesley SJ, Calvo-Garrido J, Grana O, Accari SL, Smith PK, Valencia A, Garesse R, Fisher PR, et al.: MidA is a putative methyltransferase that is required for mitochondrial complex I function. J Cell Sci 2010, 123(10):1674–1683. 10.1242/jcs.066076PubMedGoogle Scholar
- Küffner R, Rohr A, Schmiede A, Krüll C, Schulte U: Involvement of two novel chaperones in the assembly of mitochondrial NADH:Ubiquinone oxidoreductase (complex I). J Mol Biol 1998, 283(2):409–417. 10.1006/jmbi.1998.2114PubMedGoogle Scholar
- Ogilvie I, Kennaway NG, Shoubridge EA: A molecular chaperone for mitochondrial complex I assembly is mutated in a progressive encephalopathy. J Clin Invest 2005, 115(10):2784–2792. 10.1172/JCI26020PubMed CentralPubMedGoogle Scholar
- Pagliarini DJ, Calvo SE, Chang B, Sheth SA, Vafai SB, Ong S-E, Walford GA, Sugiana C, Boneh A, Chen WK, et al.: A mitochondrial protein compendium elucidates complex I disease biology. Cell 2008, 134(1):112–123. 10.1016/j.cell.2008.06.016PubMed CentralPubMedGoogle Scholar
- Saada A, Edvardson S, Rapoport M, Shaag A, Amry K, Miller C, Lorberboum-Galski H, Elpeleg O: C6ORF66 is an assembly factor of mitochondrial complex I. Am J Hum Genet 2008, 82(1):32–38. 10.1016/j.ajhg.2007.08.003PubMed CentralPubMedGoogle Scholar
- Saada A, Vogel RO, Hoefs SJ, van den Brand MA, Wessels HJ, Willems PH, Venselaar H, Shaag A, Barghuti F, Reish O, et al.: Mutations in NDUFAF3 (C3ORF60), Encoding an NDUFAF4 (C6ORF66)-Interacting Complex I Assembly Protein, Cause Fatal Neonatal Mitochondrial Disease. Am J Hum Genet 2009, 84: 718–727. 10.1016/j.ajhg.2009.04.020PubMed CentralPubMedGoogle Scholar
- Sugiana C, Pagliarini DJ, McKenzie M, Kirby DM, Salemi R, Abu-Amero KK, Dahl H-HM, Hutchison WM, Vascotto KA, Smith SM, et al.: Mutation of C20orf7 disrupts complex I assembly and causes lethal neonatal mitochondrial disease. Am J Hum Genet 2008, 83(4):468–478. 10.1016/j.ajhg.2008.09.009PubMed CentralPubMedGoogle Scholar
- Vogel RO, Janssen RJRJ, van den Brand MAM, Dieteren CEJ, Verkaart S, Koopman WJH, Willems PHGM, Pluk W, van den Heuvel LPWJ, Smeitink JAM, et al.: Cytosolic signaling protein Ecsit also localizes to mitochondria where it interacts with chaperone NDUFAF1 and functions in complex I assembly. Genes Dev 2007, 21(5):615–624. 10.1101/gad.408407PubMed CentralPubMedGoogle Scholar
- Ugalde C, Vogel R, Huijbens R, Heuvel BVD, Smeitink J, Nijtmans L: Human mitochondrial complex I assembles through the combination of evolutionary conserved modules: a framework to interpret complex I deficiencies. Hum Mol Genet 2004, 13(20):2461–2472. 10.1093/hmg/ddh262PubMedGoogle Scholar
- Gershoni M, Fuchs A, Shani N, Fridman Y, Corral-Debrinski M, Aharoni A, Frishman D, Mishmar D: Co-evolution predicts direct interactions between mtDNA and nuclear DNA-encoded subunits of oxidative phosphorylation complex I. J Mol Biol 2010, 404(1):158–171. 10.1016/j.jmb.2010.09.029PubMedGoogle Scholar
- Yamaguchi M, Hatefi Y: Mitochondrial NADH:ubiquinone oxidoreductase (complex I): proximity of the subunits of the flavoprotein and the iron-sulfur protein subcomplexes. Biochemistry 1993, 32(8):1935–1939. 10.1021/bi00059a008PubMedGoogle Scholar
- Huynen MA, Bork P: Measuring genome evolution. Proc Natl Acad Sci U S A 1998, 95(11):5849–5856. 10.1073/pnas.95.11.5849PubMed CentralPubMedGoogle Scholar
- Pellegrini M, Marcotte EM, Thompson MJ, Eisenberg D, Yeates TO: Assigning protein functions by comparative genome analysis: protein phylogenetic profiles. Proc Natl Acad Sci U S A 1999, 96(8):4285–4288. 10.1073/pnas.96.8.4285PubMed CentralPubMedGoogle Scholar
- Gabaldón T, Huynen MA: Lineage-specific gene loss following mitochondrial endosymbiosis and its potential for function prediction in eukaryotes. Bioinformatics 2005, 21(Suppl 2):ii144-ii150. 10.1093/bioinformatics/bti1124PubMedGoogle Scholar
- Kensche PR, van Noort V, Dutilh BE, Huynen MA: Practical and theoretical advances in predicting the function of a protein by its phylogenetic distribution. J R Soc Interface 2008, 5(19):151–170. 10.1098/rsif.2007.1047PubMed CentralPubMedGoogle Scholar
- Pazos F, Valencia A: Similarity of phylogenetic trees as indicator of protein-protein interaction. Protein Eng 2001, 14(9):609–614. 10.1093/protein/14.9.609PubMedGoogle Scholar
- Göbel U, Sander C, Schneider R, Valencia A: Correlated mutations and residue contacts in proteins. Proteins 1994, 18(4):309–317. 10.1002/prot.340180402PubMedGoogle Scholar
- Neher E: How frequent are correlated changes in families of protein sequences? Proc Natl Acad Sci U S A 1994, 91(1):98–102. 10.1073/pnas.91.1.98PubMed CentralPubMedGoogle Scholar
- Burger L, van Nimwegen E: Accurate prediction of protein-protein interactions from sequence alignments using a Bayesian method. Mol Syst Biol 2008, 4: 165.PubMed CentralPubMedGoogle Scholar
- Weigt M, White RA, Szurmant H, Hoch JA, Hwa T: Identification of direct residue contacts in protein-protein interaction by message passing. Proc Natl Acad Sci U S A 2009, 106(1):67–72. 10.1073/pnas.0805923106PubMed CentralPubMedGoogle Scholar
- Huynen MA, de Hollander M, Szklarczyk R: Mitochondrial proteome evolution and genetic disease. Biochim Biophys Acta 2009, 1792(12):1122–1129. 10.1016/j.bbadis.2009.03.005PubMedGoogle Scholar
- Yang M, Ge Y, Wu J, Xiao J, Yu J: Coevolution study of mitochondria respiratory chain proteins: Toward the understanding of protein-protein interaction. J Genet Genomics 2011, 38(5):201–207. 10.1016/j.jgg.2011.04.003PubMedGoogle Scholar
- Pazos F, Ranea JA, Juan D, Sternberg MJ: Assessing protein co-evolution in the context of the tree of life assists in the prediction of the interactome. J Mol Biol 2005, 352(4):1002–1015. 10.1016/j.jmb.2005.07.005PubMedGoogle Scholar
- Sato T, Yamanishi Y, Kanehisa M, Toh H: The inference of protein-protein interactions by co-evolutionary analysis is improved by excluding the information about the phylogenetic relationships. Bioinformatics 2005, 21(17):3482–3489. 10.1093/bioinformatics/bti564PubMedGoogle Scholar
- Juan D, Pazos F, Valencia A: High-confidence prediction of global interactomes based on genome-wide coevolutionary networks. Proc Natl Acad Sci 2008, 105(3):934–939. 10.1073/pnas.0709671105PubMed CentralPubMedGoogle Scholar
- Criscuolo A, Gribaldo S: BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol Biol 2010, 10(1):210. 10.1186/1471-2148-10-210PubMed CentralPubMedGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22(21):2688–2690. 10.1093/bioinformatics/btl446PubMedGoogle Scholar
- Hakes L, Lovell SC, Oliver SG, Robertson DL: Specificity in protein interactions and its relationship with sequence diversity and coevolution. Proc Natl Acad Sci U S A 2007, 104(19):7999–8004. 10.1073/pnas.0609962104PubMed CentralPubMedGoogle Scholar
- Clark NL, Alani E, Aquadro CF: Evolutionary rate covariation reveals shared functionality and coexpression of genes. Genome Res 2012, 22(4):714–720. 10.1101/gr.132647.111PubMed CentralPubMedGoogle Scholar
- Everitt BS, Dunn G: Applied multivariate data analysis /Brian S. Everitt and Graham Dunn. Edward Arnold, London ; Melbourne; 1991.Google Scholar
- Marques I, Dencher NA, Videira A, Krause F: Supramolecular Organization of the Respiratory Chain in Neurospora crassa Mitochondria. Eukaryot Cell 2007, 6(12):2391–2405. 10.1128/EC.00149-07PubMed CentralPubMedGoogle Scholar
- Marques I, Ushakova AV, Duarte M, Videira A: Role of the Conserved Cysteine Residues of the 11.5 kDa Subunit in Complex I Catalytic Properties. J Biochem 2007, 141(4):489–493. 10.1093/jb/mvm049PubMedGoogle Scholar
- Fearnley IM, Finel M, Skehel JM, Walker JE: NADH:ubiquinone oxidoreductase from bovine heart mitochondria. cDNA sequences of the import precursors of the nuclear-encoded 39 kDa and 42 kDa subunits. Biochem J 1991, 278: 821–829.PubMed CentralPubMedGoogle Scholar
- Abdrakhmanova A, Zwicker K, Kerscher S, Zickermann V, Brandt U: Tight binding of NADPH to the 39-kDa subunit of complex I is not required for catalytic activity but stabilizes the multiprotein complex. Biochimica et Biophysica Acta (BBA). Bioenergetics 2006, 1757(12):1676–1682. 10.1016/j.bbabio.2006.09.003Google Scholar
- Antonicka H, Ogilvie I, Taivassalo T, Anitori RP, Haller RG, Vissing J, Kennaway NG, Shoubridge EA: Identification and characterization of a common set of complex I assembly intermediates in mitochondria from patients with complex I deficiency. J Biol Chem 2003, 278(44):43081–43088. 10.1074/jbc.M304998200PubMedGoogle Scholar
- Heinrich H, Azevedo JE, Werner S: Characterization of the 9.5-kDa ubiquinone-binding protein of NADH:ubiquinone oxidoreductase (complex I) from Neurospora crassa. Biochemistry 1992, 31(46):11420–11424. 10.1021/bi00161a021PubMedGoogle Scholar
- Angerer H, Zwicker K, Wumaier Z, Sokolova L, Heide H, Steger M, Kaiser S, Nübel E, Brutschy B, Radermacher M, et al.: A scaffold of accessory subunits links the peripheral arm and the distal proton pumping module of mitochondrial complex I. Biochem J 2011, 437: 279–288. 10.1042/BJ20110359PubMed CentralPubMedGoogle Scholar
- da Silva MV, Alves PC, Duarte M, Mota N, da Cunha AL, Harkness TA, Nargang FE, Videira A: Disruption of the nuclear gene encoding the 20.8-kDa subunit of NADH: ubiquinone reductase of Neurospora mitochondria. Mol Gen Genet 1996, 252(1–2):177–183.PubMedGoogle Scholar
- Nehls U, Friedrich T, Schmiede A, Ohnishi T, Weiss H: Characterization of assembly intermediates of NADH:ubiquinone oxidoreductase (complex I) accumulated in Neurospora mitochondria by gene disruption. J Mol Biol 1992, 227(4):1032–1042. 10.1016/0022-2836(92)90519-PPubMedGoogle Scholar
- Lynch M, Koskella B, Schaack S: Mutation Pressure and the Evolution of Organelle Genomic Architecture. Science 2006, 311(5768):1727–1730. 10.1126/science.1118884PubMedGoogle Scholar
- Sazanov LA, Peak-Chew SY, Fearnley IM, Walker JE: Resolution of the membrane domain of bovine complex I into subcomplexes: implications for the structural organization of the enzyme. Biochemistry 2000, 39(24):7229–7235. 10.1021/bi000335tPubMedGoogle Scholar
- Lemma-Gray P, Valusová E, Carroll CA, Weintraub ST, Musatov A, Robinson NC: Subunit analysis of bovine heart complex I by reversed-phase high-performance liquid chromatography, electrospray ionization-tandem mass spectrometry, and matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry. Anal Biochem 2008, 382(2):116–121. 10.1016/j.ab.2008.07.029PubMed CentralPubMedGoogle Scholar
- Marques I, Duarte M, Videira A: The 9.8 kDa Subunit of Complex I, Related to Bacterial Na+−translocating NADH Dehydrogenases, is Required for Enzyme Assembly and Function in Neurospora crassa. J Mol Biol 2003, 329(2):283–290. 10.1016/S0022-2836(03)00443-1PubMedGoogle Scholar
- Dudkina NV, Kouřil R, Peters K, Braun H-P, Boekema EJ: Structure and function of mitochondrial supercomplexes. Biochim Biophys Acta 2010, 1797(6–7):664–670.PubMedGoogle Scholar
- Reyes A, Gissi C, Pesole G, Saccone C: Asymmetrical directional mutation pressure in the mitochondrial genome of mammals. Mol Biol Evol 1998, 15(8):957–966. 10.1093/oxfordjournals.molbev.a026011PubMedGoogle Scholar
- Delorme MO, Hénaut A: Codon usage is imposed by the gene location in the transcription unit. Curr Genet 1991, 20(5):353–358. 10.1007/BF00317061PubMedGoogle Scholar
- Asakawa S, Kumazawa Y, Araki T, Himeno H, Miura K, Watanabe K: Strand-specific nucleotide composition bias in echinoderm and vertebrate mitochondrial genomes. J Mol Evol 1991, 32(6):511–520. 10.1007/BF02102653PubMedGoogle Scholar
- Wei S-J, Shi M, Chen X-X, Sharkey MJ, van Achterberg C, Ye G-Y, He J-H: New views on strand asymmetry in insect mitochondrial genomes. PLoS One 2010, 5(9):e12708. 10.1371/journal.pone.0012708PubMed CentralPubMedGoogle Scholar
- Nosek J, Tomaska L: Mitochondrial genome diversity: evolution of the molecular architecture and replication strategy. Curr Genet 2003, 44(2):73–84. 10.1007/s00294-003-0426-zPubMedGoogle Scholar
- Nikolaou C, Almirantis Y: Deviations from Chargaff's second parity rule in organellar DNA: Insights into the evolution of organellar genomes. Gene 2006, 381: 34–41.PubMedGoogle Scholar
- Lazarou M, Thorburn DR, Ryan MT, McKenzie M: Assembly of mitochondrial complex I and defects in disease. Biochim Biophys Acta, Mol Cell Res 2009, 1793(1):78–88. 10.1016/j.bbamcr.2008.04.015PubMedGoogle Scholar
- Dunning CJR, McKenzie M, Sugiana C, Lazarou M, Silke J, Connelly A, Fletcher JM, Kirby DM, Thorburn DR, Ryan MT: Human CIA30 is involved in the early assembly of mitochondrial complex I and mutations in its gene cause disease. EMBO J 2007, 26(13):3227–3237. 10.1038/sj.emboj.7601748PubMed CentralPubMedGoogle Scholar
- Lazarou M, McKenzie M, Ohtake A, Thorburn DR, Ryan MT: Analysis of the assembly profiles for mitochondrial- and nuclear-DNA-encoded subunits into complex I. Mol Cell Biol 2007, 27(12):4228–4237. 10.1128/MCB.00074-07PubMed CentralPubMedGoogle Scholar
- Vogel RO, van den Brand MAM, Rodenburg RJ, van den Heuvel LPWJ, Tsuneoka M, Smeitink JAM, Nijtmans LGJ: Investigation of the complex I assembly chaperones B17.2L and NDUFAF1 in a cohort of CI deficient patients. Mol Genet Metab 2007, 91(2):176–182. 10.1016/j.ymgme.2007.02.007PubMedGoogle Scholar
- McKenzie M, Tucker EJ, Compton AG, Lazarou M, George C, Thorburn DR, Ryan MT: Mutations in the Gene Encoding C8orf38 Block Complex I Assembly by Inhibiting Production of the Mitochondria-Encoded Subunit ND1. J Mol Biol 2011, 414(3):413–426. 10.1016/j.jmb.2011.10.012PubMedGoogle Scholar
- Sheftel AD, Stehling O, Pierik AJ, Netz DJA, Kerscher S, Elsasser H-P, Wittig I, Balk J, Brandt U, Lill R: Human Ind1, an Iron-Sulfur Cluster Assembly Factor for Respiratory Complex I. Mol Cell Biol 2009, 29(22):6059–6073. 10.1128/MCB.00817-09PubMed CentralPubMedGoogle Scholar
- Dieteren CEJ, Willems PHGM, Vogel RO, Swarts HG, Fransen J, Roepman R, Crienen G, Smeitink JAM, Nijtmans LGJ, Koopman WJH: Subunits of Mitochondrial Complex I Exist as Part of Matrix- and Membrane-associated Subcomplexes in Living Cells. J Biol Chem 2008, 283(50):34753–34761. 10.1074/jbc.M807323200PubMed CentralPubMedGoogle Scholar
- Brockmann C, Diehl A, Rehbein K, Strauss H, Schmieder P, Korn B, Kühne R, Oschkinat H: The oxidized subunit B8 from human complex I adopts a thioredoxin fold. Structure 2004, 12(9):1645–1654. 10.1016/j.str.2004.06.021PubMedGoogle Scholar
- Brandt U: Energy converting NADH:quinone oxidoreductase (complex I). Annu Rev Biochem 2006, 75: 69–92. 10.1146/annurev.biochem.75.103004.142539PubMedGoogle Scholar
- Johnson DC, Dean DR, Smith AD, Johnson MK: Structure, function, and formation of biological iron-sulfur clusters. Annu Rev Biochem 2005, 74: 247–281. 10.1146/annurev.biochem.74.082803.133518PubMedGoogle Scholar
- Szklarczyk R, Huynen M, Snel B: Complex fate of paralogs. BMC Evol Biol 2008, 8(1):337. 10.1186/1471-2148-8-337PubMed CentralPubMedGoogle Scholar
- Pál C, Papp B, Lercher MJ: An integrated view of protein evolution. Nat Rev Genet 2006, 7(5):337–348. 10.1038/nrg1838PubMedGoogle Scholar
- Drummond DA, Raval A, Wilke CO: A single determinant dominates the rate of yeast protein evolution. Mol Biol Evol 2006, 23(2):327–337.PubMedGoogle Scholar
- Schug A, Weigt M, Onuchic JN, Hwa T, Szurmant H: High-resolution protein complexes from integrating genomic information with molecular simulation. Proc Natl Acad Sci U S A 2009, 106(52):22124–22129. 10.1073/pnas.0912100106PubMed CentralPubMedGoogle Scholar
- Morcos F, Pagnani A, Lunt B, Bertolino A, Marks DS, Sander C, Zecchina R, Onuchic JN, Hwa T, Weigt M: Direct-coupling analysis of residue coevolution captures native contacts across many protein families. Proc Natl Acad Sci U S A 2011, 108(49):E1293-E1301. 10.1073/pnas.1111471108PubMed CentralPubMedGoogle Scholar
- Sreekumar J, ter Braak CJ, van Ham RC, van Dijk AD: Correlated mutations via regularized multinomial regression. BMC Bioinforma 2011, 12: 444. 10.1186/1471-2105-12-444Google Scholar
- Burger L, van Nimwegen E: Disentangling direct from indirect co-evolution of residues in protein alignments. PLoS Comput Biol 2010, 6(1):e1000633. 10.1371/journal.pcbi.1000633PubMed CentralPubMedGoogle Scholar
- Pazos F, Valencia A: In silico two-hybrid system for the selection of physically interacting protein pairs. Proteins 2002, 47(2):219–227. 10.1002/prot.10074PubMedGoogle Scholar
- Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, Federhen S, et al.: Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 2006, 35(suppl 1):D5-D12.PubMed CentralPubMedGoogle Scholar
- Söding J: Protein homology detection by HMM–HMM comparison. Bioinformatics 2005, 21(7):951–960. 10.1093/bioinformatics/bti125PubMedGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet 2000, 16(6):276–277. 10.1016/S0168-9525(00)02024-2PubMedGoogle Scholar
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinform 2008, 9(4):286–298. 10.1093/bib/bbn013PubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22(22):4673–4680. 10.1093/nar/22.22.4673PubMed CentralPubMedGoogle Scholar
- Li K-B: ClustalW-MPI: ClustalW analysis using distributed and parallel computing. Bioinformatics 2003, 19(12):1585–1586. 10.1093/bioinformatics/btg192PubMedGoogle Scholar
- Eddy SR: Profile hidden Markov models. Bioinformatics 1998, 14(9):755–763. 10.1093/bioinformatics/14.9.755PubMedGoogle Scholar
- Strang G: Introduction to Linear Algebra. Wellesley-Cambridge Press, In.; 1998.Google Scholar
- Kann MG, Shoemaker BA, Panchenko AR, Przytycka TM: Correlated evolution of interacting proteins: looking behind the mirrortree. J Mol Biol 2009, 385(1):91–98. 10.1016/j.jmb.2008.09.078PubMed CentralPubMedGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria; 2011.Google Scholar
- Borg I, Groenen PJF: Modern Multidimensional Scaling. Theory and Applications, Springer; 2005.Google Scholar
- Singer MS, Vriend G, Bywater RP: Prediction of protein residue contacts with a PDB-derived likelihood matrix. Protein Eng 2002, 15(9):721–725. 10.1093/protein/15.9.721PubMedGoogle Scholar
- Lê S, Josse J, Husson F: FactoMineR: An R Package for Multivariate Analysis. J Stat Softw 2008, 25(1):1–18.Google Scholar
- Vogel RO, Smeitink JAM, Nijtmans LGJ: Human mitochondrial complex I assembly: a dynamic and versatile process. Biochim Biophys Acta 2007, 1767(10):1215–1227. 10.1016/j.bbabio.2007.07.008PubMedGoogle Scholar
- Carroll J, Shannon RJ, Fearnley IM, Walker JE, Hirst J: Definition of the nuclear encoded protein composition of bovine heart mitochondrial complex I. Identification of two new subunits. J Biol Chem 2002, 277(52):50311–50317. 10.1074/jbc.M209166200PubMedGoogle Scholar
- Fearnley IM, Carroll J, Shannon RJ, Runswick MJ, Walker JE, Hirst J: GRIM-19, a cell death regulatory gene product, is a subunit of bovine mitochondrial NADH:ubiquinone oxidoreductase (complex I). J Biol Chem 2001, 276(42):38345–38348. 10.1074/jbc.C100444200PubMedGoogle Scholar
- Rual J-F, Venkatesan K, Hao T, Hirozane-Kishikawa T, Dricot A, Li N, Berriz GF, Gibbons FD, Dreze M, Ayivi-Guedehoussou N, et al.: Towards a proteome-scale map of the human protein-protein interaction network. Nature 2005, 437(7062):1173–1178. 10.1038/nature04209PubMedGoogle Scholar
- Murray J, Taylor SW, Zhang B, Ghosh SS, Capaldi RA: Oxidative Damage to Mitochondrial Complex I Due to Peroxynitrite. J Biol Chem 2003, 278(39):37223–37230. 10.1074/jbc.M305694200PubMedGoogle Scholar
- Mishmar D, Ruiz-Pesini E, Mondragon-Palomino M, Procaccio V, Gaut B, Wallace DC: Adaptive selection of mitochondrial complex I subunits during primate radiation. Gene 2006, 378: 11–18.PubMedGoogle Scholar
- Nouws J, Nijtmans L, Houten SM, van den Brand M, Huynen M, Venselaar H, Hoefs S, Gloerich J, Kronick J, Hutchin T, et al.: Acyl-CoA dehydrogenase 9 is required for the biogenesis of oxidative phosphorylation complex I. Cell Metab 2010, 12(3):283–294. 10.1016/j.cmet.2010.08.002PubMedGoogle Scholar
- Szklarczyk D, Franceschini A, Kuhn M, Simonovic M, Roth A, Minguez P, Doerks T, Stark M, Muller J, Bork P: The STRING database in 2011: functional interaction networks of proteins, globally integrated and scored. Nucleic Acids Res 2011, 39(Database issue):D561–568.PubMed CentralPubMedGoogle Scholar