Jun 10, 2009 - Ernster, L., Hoberman, H. D., Howard, R. L., King, T. E., Lee, C. P., Mack- ler, B., and Sottocasa, ..... Ben-Tal, N. (2005) Nucl. Acids Res.
THE JOURNAL OF BIOLOGICAL CHEMISTRY VOL. 284, NO. 43, pp. 29773–29783, October 23, 2009 © 2009 by The American Society for Biochemistry and Molecular Biology, Inc. Printed in the U.S.A.
Structural Basis for the Mechanism of Respiratory Complex I*□ S
Received for publication, June 10, 2009, and in revised form, July 15, 2009 Published, JBC Papers in Press, July 27, 2009, DOI 10.1074/jbc.M109.032144
John M. Berrisford and Leonid A. Sazanov1 From the Medical Research Council Mitochondrial Biology Unit, Cambridge CB2 0XY, United Kingdom
Complex I (NADH:ubiquinone oxidoreductase, EC 1.6.5.3) is the first enzyme of the mitochondrial and bacterial respiratory chains. It catalyzes the transfer of two electrons from NADH to quinone, coupled to the translocation of approximately four protons across the membrane, contributing to the proton-motive force required for the synthesis of ATP (1, 2). The mitochondrial enzyme consists of 45 subunits (3) with a combined mass of ⬃980 kDa. The prokaryotic enzyme is simpler, consisting of ⬃14 subunits conserved from bacteria to humans, and has a total mass of ⬃550 kDa (2). The mitochon-
* This work was supported by the Medical Research Council. □ S
The on-line version of this article (available at http://www.jbc.org) contains supplemental text, Table S1, and Figs. S1–S5. The atomic coordinates and structure factors (codes 3IAS, 3I9V, and 3IAM) have been deposited in the Protein Data Bank, Research Collaboratory for Structural Bioinformatics, Rutgers University, New Brunswick, NJ (http://www.rcsb.org/). 1 To whom correspondence should be addressed: Medical Research Council Mitochondrial Biology Unit, Wellcome Trust/MRC Bldg., Hills Road, Cambridge CB2 0XY, UK. Fax: 44-1223-252915; E-mail: sazanov@mrc-mbu. cam.ac.uk.
OCTOBER 23, 2009 • VOLUME 284 • NUMBER 43
drial and bacterial enzymes contain equivalent redox components and have a similar L-shaped structure, with the hydrophobic arm embedded in the membrane and the hydrophilic peripheral arm protruding into the mitochondrial matrix or the bacterial cytoplasm (2, 4). Thus, the bacterial enzyme represents a “minimal” model of complex I. Because of the central role of complex I in respiration, mutations in individual subunits can lead to many human neurodegenerative diseases (5). Complex I, along with complex III (bc1), has been suggested to be a major source of reactive oxygen species (ROS)2 in mitochondria, which can damage mitochondrial DNA and may be one of the causes of aging (6). Parkinson disease, at least in its sporadic form (which represents ⬃95% of cases), may be caused by increased ROS production from malfunctioning complex I (7). We have previously determined the crystal structure of the hydrophilic domain (eight different subunits of 280 kDa total mass) of complex I from Thermus thermophilus, establishing the electron transfer pathway from the primary electron acceptor flavin mononucleotide (FMN) through seven conserved iron-sulfur clusters to the quinone-binding site (Q-site) at the interface with the membrane domain (8). Two additional ironsulfur clusters, which are not part of the main redox chain, may represent an evolutionary remnant (cluster N7) and a possible anti-oxidant (cluster N1a; cluster names are assigned to structural motifs as in Ref. 8). The membrane-spanning part of the enzyme lacks covalently bound prosthetic groups (9) but must contain essential components of the proton translocating machinery. Its atomic structure is currently unknown. The mechanism of the highly efficient coupling between electron transfer and proton pumping, conserving nearly 100% of the available energy, remains a mystery. Two models are being discussed: direct (redox-driven through chemical intermediates, usually employing modifications of the Q cycle, with quinol as a mobile proton/electron carrier) (10) and indirect or conformation-driven coupling (2, 4, 11, 12). Sequence comparisons indicate that the three largest hydrophobic subunits of complex I, Nqo12, 13, and 14 (Thermus nomenclature), are homologous to each other and to the antiporter family (Mrp) (13, 14) and so are likely to participate in proton translocation. Two of these subunits, Nqo12 and Nqo13, are located ⬃100 Å away from the Q-site (15), which implies the need for conformational coupling as at least a part of the mechanism. We have 2
The abbreviations used are: ROS, reactive oxygen species; FMN, flavin mononucleotide; Q-site, quinone-binding site; bis-tris, 2-[bis(2-hydroxyethyl)amino]-2-(hydroxymethyl)propane-1,3-diol; ASU, asymmetric unit; NCS, noncrystallographic symmetry; structures of hydrophilic domain: RND, reduced by NADH and dithionite; RN, reduced by NADH; RD, reduced by dithionite; RA1/2, reduced by NADH and exposed to air; O2, oxidized, 2 mol/ASU; O4, oxidized, 4 mol/ASU.
JOURNAL OF BIOLOGICAL CHEMISTRY
29773
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
Complex I plays a central role in cellular energy production, coupling electron transfer between NADH and quinone to proton translocation. The mechanism of this highly efficient enzyme is currently unknown. Mitochondrial complex I is a major source of reactive oxygen species, which may be one of the causes of aging. Dysfunction of complex I is implicated in many human neurodegenerative diseases. We have determined several x-ray structures of the oxidized and reduced hydrophilic ˚ domain of complex I from Thermus thermophilus at up to 3.1 A resolution. The structures reveal the mode of interaction of complex I with NADH, explaining known kinetic data and providing implications for the mechanism of reactive oxygen species production at the flavin site of complex I. Bound metals were identified in the channel at the interface with the frataxinlike subunit Nqo15, indicating possible iron-binding sites. Conformational changes upon reduction of the complex involve adjustments in the nucleotide-binding pocket, as well as small but significant shifts of several ␣-helices at the interface with the membrane domain. These shifts are likely to be driven by the reduction of nearby iron-sulfur clusters N2 and N6a/b. Cluster N2 is the electron donor to quinone and is coordinated by unique motif involving two consecutive (tandem) cysteines. An unprecedented “on/off switch” (disconnection) of coordinating bonds between the tandem cysteines and this cluster was observed upon reduction. Comparison of the structures suggests a novel mechanism of coupling between electron transfer and proton translocation, combining conformational changes and protonation/deprotonation of tandem cysteines.
Mechanism of Respiratory Complex I now determined several structures of the oxidized and reduced hydrophilic domain of complex I from T. thermophilus, which show how NADH interacts with the complex and provide novel insights into the coupling mechanism.
29774 JOURNAL OF BIOLOGICAL CHEMISTRY
VOLUME 284 • NUMBER 43 • OCTOBER 23, 2009
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
EXPERIMENTAL PROCEDURES Crystallization—The hydrophilic domain of complex I from T. thermophilus was purified and crystallized as described previously (8, 16). To improve the diffraction properties of the crystals, a wide range of additives were tested. The most useful was the replacement of CaCl2 for either MgCl2 or MnCl2, as well as the addition of low amounts of polyethylene glycol 400 in some cases. To obtain structures of the reduced domain, extensive crystallization set-ups were performed in the presence of NADH and/or sodium dithionite. To avoid oxidative damage and enzymatic turnover with oxygen, this was performed anaerobically using a Belle Technologies (Dorset, UK) glove box with oxygen levels below 20 ppm. Crystallization in the presence of NAD⫹ was performed aerobically, with the enzyme oxidized as purified. The crystals were grown at 22 °C by the vapor diffusion method, by mixing 1 l of protein solution (12–15 mg/ml, in 20 mM bis-tris propane pH 7.3, 150 mM NaCl, 0.002% phenylmethylsulfonyl fluoride, and 1% octyl--D-glucopyranoside) with 1 l of reservoir solution (supplemental Table S1). The crystals were fully formed as large thin plates after ⬃5 days. They were cryo-protected by adding harvest solution containing reservoir solution with polyethylene glycol concentration increased by 2% and with the addition of 1% octyl--D-glucopyranoside and 25% ethylene glycol. The crystals were frozen by plunging them into liquid nitrogen. For reduced protein it was important to minimize the exposure of crystals to air during harvesting and freezing, with exposures of longer than 5 min leading to a gradual decline in the diffraction resolution and exposures of more than 20 –30 min resulting in the gradual loss of electron density for the bound NADH. Crystals frozen with air exposures of less than 5 min did not suffer any loss in resolution and are referred to as being frozen anaerobically (“anaerobic” state). We have also shown by spectroscopic measurements that sodium dithionite at concentrations used for crystal harvesting (such as for RND and RD structures) remains nearly fully reduced for more than 5 min after exposure to air, thus preserving the reducing environment. Crystals frozen after exposure to air for ⬃30 min are referred to as being frozen aerobically (“aerobic” state). Data Processing and Structure Determination—Data were collected at 100 K with ADSC Q315 CCD detectors at Beamlines ID29 and ID23-1 at the European Synchrotron Radiation Facility (Grenoble, France) or at I02 Beamline at the Diamond Light Source (Didcot, UK). The image data were processed with MOSFLM and SCALA from the CCP4 suite (17). The crystals belonged to space group P21, with one exception (data set RD is in P1). For a range of data sets at intermediate resolutions (between 4.5 and 3.5 Å), electron density for NADH was usually present in molecular replacement solutions (performed with PHASER (18) using the known structure as a search model (Protein Data Bank code 2FUG (8)), whereas NAD⫹ was not seen despite being added to mother liquor at up to 30 mM concentration. This is consistent with low affinity of NAD⫹ to the
oxidized enzyme (⬃1–2 mM (19)), which most likely results from overall shift in the nucleotide charge toward positive (compared with NADH), in the vicinity of the positively charged binding site surface (4). In contrast, NADH has high affinity for both the oxidized (⬃10 M) and the reduced enzyme (⬃100 M) (4, 19). The native data used for the previously published structure (protein oxidized as purified, 4 molecules/ASU (8)) had significant anisotropy and had subsequently been subjected to anisotropic scaling and ellipsoidal truncation using the diffraction anisotropy server (20), to the resolution limits of 3.15, 3.15, and 3.5 Å, where the F/ ratio drops to ⬃2.6 along a*, b*, and c*, respectively (data set O4). Experimental electron density was then calculated with this data and previous iron peak/derivative data (8) using SHARP (21) and solvent flattening in DM from the CCP4 suite, as described previously (8). This lead to improved experimental and modeled electron density. Two of the new data sets, RND (supplemental Table S1), from the crystal of the reduced domain grown with NADH/dithionite, and data set O2, from the crystal of oxidized domain grown in the presence of NAD⫹ (but not containing any bound ordered nucleotide), were both in a novel packing arrangement containing 2 molecules/ASU, and both extended to 3.1 Å resolution. The O2 data set was collected from two different areas on the same crystal because of radiation damage, and so R factors for this structure are somewhat higher than for the others. This structure is especially useful for comparisons because the crystal packing is the same as in most structures of the reduced enzyme. For data set RND, data from the same crystal were also collected at the iron peak wavelength. The initial positions of the iron-sulfur clusters were determined by molecular replacement with PHASER using the known structure as a search model (Protein Data Bank code 2FUG (8)). MAD phasing in SHARP, using iron peak and native data sets, was then performed essentially as previously described (8). Subsequent density modification using solvent flipping in CNS (22) produced significantly better maps than solvent flattening in DM, probably because only two molecules were available for NCS averaging, as compared with our original case of four molecules (when there was no significant difference between the two approaches). CNS maps were of very good quality (supplemental Fig. S1), showing the power of two-wavelength MAD phasing using iron-sulfur clusters. This allowed us to avoid model bias during the solution of the structure of the reduced domain. Improved resolution and the availability of experimental density in two different crystal forms allowed us to increase the completeness of the model and to improve its accuracy in some of surface-exposed, not very well ordered areas. In particular, residues 226 –243 from subunit Nqo4 and residues 714 –724 and 748 –758 from Nqo3 were retraced. The register was adjusted in some residues neighboring these areas, as well as in residues 2–26 from Nqo1, residues 173–180 and 511–527 from Nqo3, residues 348 –350 from Nqo4, and residues 2–13 and 107–140 from Nqo5. These changes do not involve any of the functionally important residues discussed previously (4, 8). Overall, these adjustments lead to improved geometry and lower R factors for the oxidized domain structure and for good quality models in new crystal forms (supplemental Table S1).
Mechanism of Respiratory Complex I
3
The prefixes before the residue names indicate the subunit numbers.
OCTOBER 23, 2009 • VOLUME 284 • NUMBER 43
(data set RD). Furthermore, two data sets at 4.2 Å resolution (RA1 and RA2 and supplemental Table S1) were collected from aerobically frozen crystals of NADH-reduced domain. Thus, we compared two structures of the oxidized enzyme, O4 and O2, with five structures of the reduced enzyme in somewhat different states. All low resolution structures (data sets RN, RD, RA1, and RA2) were refined by the same procedure in CNS, optimized by following Rfree. Molecular replacement solutions (obtained using the O2 structure of the oxidized domain) were subjected to a round of torsion angle dynamics-simulated annealing with starting temperature of 2500 K, followed by slow cooling and two macro cycles of 20 coordinate minimization steps and five group B factor refinement steps. NCS restraints with 100 kcal/mol/Å2 weight were used throughout. To keep good model geometry it was important to use phi/psi restraints to the starting structure angles, as implemented in CNS. These restraints also lowered the Rfree, and overall such a procedure resulted in a large drop of Rfree (between 3 and 8%) from the initial values. For structures with bound NADH, the refinement resulted in a movement of Nqo1 loop 202–207 and in shifts of neighboring ␣-helices similar to those seen in the higher resolution RND structure, validating such an approach. If the RND structure was used as a starting model instead of O2 (as a control), the final structures after refinement with any particular data set were very similar, again validating the procedure used. Overall, the refined structures fitted well to the calculated electron density. These low resolution structures were not rebuilt manually, because the main purpose was to check for gross shifts of ␣-helices, which were indeed observed. For comparisons, the structures were superposed using CCP4 program LSQKAB (17). The figures were prepared with PyMOL (DeLano Scientific).
RESULTS AND DISCUSSION Structures—To determine the structure of the reduced domain, the crystals were grown anaerobically in the presence of NADH and/or dithionite. The structure with bound substrate NADH has been determined to 3.1 Å resolution in a novel packing arrangement of 2 mol/ASU in space group P21 (NADH and dithionite, data set/structure RND; supplemental Table S1). Additional structures were determined to 3.7 Å resolution: RN (NADH, 2 mol/ASU in P21) and RD (dithionite, 4 mol/ASU in P1) (supplemental Table S1). These crystals were frozen with minimal exposure to oxygen (⬍5 min) during cryo protection to ensure that a fully reduced (anaerobic) state is represented. Some crystals were also grown anaerobically in the presence of NADH but were harvested with longer exposures to oxygen (⬃30 min). The resulting structures (NADH aerobic, 4.2 Å resolution, 4 mol/ASU in P21, data sets RA1 and RA2; supplemental Table S1) were included in this analysis because they appear to represent a redox state similar to the RD structure. One of the crystals was grown aerobically in the presence of NAD⫹, and the resulting structure of the oxidized domain represents, usefully, a similar packing arrangement to the reduced domain structures (2 mol/ASU in P21, 3.1 Å resolution, data set O2; supplemental Table S1). It does not contain any ordered NAD⫹ because of its low affinity to the oxidized enzyme. The resolution of the original structure of the oxidized domain (4 mol/ JOURNAL OF BIOLOGICAL CHEMISTRY
29775
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
The improved resolution also allowed us to observe a partially disordered loop between two N-terminal -strands in Nqo4, which led us to adjust register adjustment in these strands. Relatively poor electron density of this loop (residues 32–38) did not allow us to build it unambiguously, and so it is not included in the final models. This loop is facing the Q cavity in the vicinity of functionally important 4Tyr87.3 Although the loop is partly disordered, the remaining electron density appears to be stronger in structures of the oxidized enzyme compared with the reduced state. This flexible loop contains two residues that are highly conserved and essential for activity (23), histidines 34 and 38, which may serve as intermediates in quinone protonation and could interact with 4Tyr87. As noted previously, in all structures some of the residues with flexible side chains, exposed to the solvent at the surface of the molecule, have clear electron density for backbone atoms, but not for side chains. These have not been mutated to alanine but can be identified by high B factors (more than ⬃110 Å2). In the O2 and RND structures, unexplained electron density is present near the entrance to the putative iron-binding channel, in the vicinity of 2Glu94 and 2Asp98. It might represent several manganese or magnesium atoms (or possibly polyethylene glycol 400 in the O2 structure) and was not modeled because of the irregular shape of the density. In the RND structure, 6Cys46 was not modeled as disconnected from cluster N2, despite the lack of electron density, because the nature of the replacement iron ligand is unknown. However, the length of the iron-sulfur bond was increased in refinement to place the cysteine backbone atoms closer to their electron density. Model building was carried out in O (24), using the PHASER molecular replacement solutions as starting points. The higher resolution structures (O4, O2, and RND) were refined in CNS version 1.2 (22) using NCS restraints with 100 kcal/mol/Å2 weight and restrained individual B factor refinement. Experimental phase information was used for refinement of the structures O4 and RND. The final round of refinement was carried out in Phenix (25), which resulted in slightly improved R factors and geometry. According to PROCHECK (26) and MOLPROBITY (27) analyses, these structures are of better than average quality of similar resolution structures. The models of the oxidized domain (O4 and O2) consist of the following residues: subunit Nqo1, 2– 438; Nqo2, 3–180; Nqo3, 1–55, 73–143, and 150 –777; Nqo4, 26 –31 and 39 – 409; Nqo5, 1–196; Nqo6, 15–57 and 74 –175; Nqo9, 26 –179; and Nqo15, 3–129. Each model also contains seven [4Fe-4S] clusters, two [2Fe-2S] clusters, and one FMN molecule. The structure of the reduced domain (RND) contains an NADH molecule. It also contains residues 2Gly2 and 4Val25, whereas residue 6Arg57 is not included. To verify whether the observed differences between structures of the oxidized and reduced domain are reproducible, we also determined structures from crystals grown with NADH only (without dithionite) at 3.7 Å resolution (data set RN) and from crystals grown with dithionite only also at 3.7 Å resolution
Mechanism of Respiratory Complex I
29776 JOURNAL OF BIOLOGICAL CHEMISTRY
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
ASU in P21; Protein Data Bank code 2FUG) was improved to 3.15 Å (from 3.3 Å) by anisotropic scaling (supplemental text) (data set O4; supplemental Table S1). The availability of experimental electron density in two different crystal forms and the increased resolution allowed us to improve the completeness (from 2333 to 2369 residues) and the accuracy of the original model in some of surface-exposed areas (see “Experimental Procedures”). Altogether seven structures have been used in the analysis, and any features discussed here have been reproduced at least twice. Mode of NADH Binding—Electron density for a bound nucleotide, adjacent to the FMN, was identified in the RND structure (supplemental Table S1 and Fig. S1) and has been modeled as NADH (rather than NAD⫹) in accordance with crystal growth conditions (supplemental text). The FMN is also likely to be reduced in these conditions (supplemental text). Hydrophobic stacking interactions are formed between the nicotinamide ring of NADH and the exposed face of the isoalloxazine ring of the bound FMN, similarly to many other nucleotide-binding flavoenzymes (29, 30). The adenine ring of NADH is stabilized by stacking interactions with the aromatic rings of three conserved phenylalanines (Fig. 1A and supplemental Figs. S2 and S3). NADH binding induces movement of the loop formed by Nqo1 residues 202–207, which puts 1Phe205 ⬃1.7 Å closer to the adenine ring, promoting a stacking interaction (Fig. 1B). The glycine-rich loop formed by Nqo1 residues 65–71 (as well as the preceding N-terminal ␣-helical domain) moves ⬃0.7 Å away from NADH, allowing it to occupy the binding pocket. In contrast to most enzymes interacting with nucleotides (29), this exposed glycine-rich loop (which is not part of classic Rossmann fold (8)) is not involved in hydrogen bonding with the pyrophosphate moiety, despite being in close proximity to it. Instead, upon NADH binding, conserved 1Lys202 loses an H-bond to 1Glu184 and makes an H-bond to O2A of the pyrophosphate. The carbonyls of 1Gly66 and 1Gly67 from the glycine-rich loop interact with N7N of the carboxamide group and O3D of the nicotinamide ribose. Additionally, conserved 1Lys75 loses H-bonds to FMN O5⬘ and O2P and instead makes an H-bond to O3B of the adenosine ribose. The carboxyl group of conserved 1Glu185 also interacts with the adenosine ribose, in a similar manner to the aspartate conserved in most dehydrogenases (29). Altogether, ⬃10 putative hydrogen bonds are formed between NADH and the protein (Fig. 1A and supplemental Fig. S2). These hydrogen bonds and the stacking interactions extend over the whole length of nucleotide, which fits tightly in its extended form into the binding pocket, explaining the high affinity of complex I for its substrate. The O2B atom of adenosine ribose points into the protein interior, but there is some limited unoccupied space in this area, which may accommodate the extra phosphate of NADPH. This explains the observable (but low) activity of complex I with NADPH as a substrate (31, 32). The N6A atom of the adenine moiety points into solvent, without any interactions with the protein, which explains why complex I can effectively use deamino-NADH (in which this atom is replaced by oxygen), unlike many other dehydrogenases. Mechanism of Hydride Transfer—The B-face of the nicotinamide ring is stacked (slightly at an angle) against the re-face of
FIGURE 1. The NADH-binding site. A, the site viewed from the solvent-exposed side. FMN and residues involved in NADH binding are shown as sticks with carbon in yellow and NADH with carbon in salmon. Hydrogen bonds are shown as dotted lines in green, stacking interactions are in gray, and the hydride transfer path is in red. The van der Waal’s contact between 1Glu97 and C4N of NADH is shown in gray. Prefixes before residue names indicate the subunit number. B, conformational changes in the nucleotide-binding site upon reduction by NADH. Superposed structures of the oxidized (O2, green) and NADH-reduced (RND, yellow) domain are shown. C, surface representation without NADH; FMN is shown as sticks with carbon in cyan. D, surface representation with bound NADH, colored salmon.
VOLUME 284 • NUMBER 43 • OCTOBER 23, 2009
Mechanism of Respiratory Complex I
OCTOBER 23, 2009 • VOLUME 284 • NUMBER 43
a drug directed at minimizing ROS production by complex I would be a compound with high affinity for the reduced flavin site but low affinity toward oxidized enzyme, to avoid interference with the main physiological activity. Conformational Changes upon Reduction—From a mechanistic point of view it is of particular interest to know whether any long range conformational changes occur upon reduction of the enzyme. Comparison of the structure of NADH-reduced hydrophilic domain (RND) to the oxidized form (O2 or O4) does not reveal major rearrangements (root mean square deviation ⫽ ⬃0.6 Å). However, some reproducible differences are observed, mostly at the extremities of the structure (Fig. 2, A and B), near the nucleotide-binding site (tops of the panels) and near the interface with the membrane domain (bottoms of the panels). Changes near the nucleotide-binding site appear to be directly related to NADH binding, as discussed above. Near to the interface with the membrane domain, reduction by NADH leads to a shift of ⬃1 Å in the positions of helices H1 and H2 (residues 18 –35 and 143–160 of Nqo6) toward the tip of the membrane domain. There is also a concerted shift (of ⬃1 Å) of a four-helix bundle in Nqo4 (residues 110 –297) toward the membrane interior (Fig. 2C). Even though these shifts are small, they are determined unambiguously at our resolution and are reproducible (supplemental text). The presence or absence of dithionite did not lead to any significant differences between the two structures of the NADH-reduced domain (RND and RN). Cross-linking studies have suggested that, upon reduction, contacts between the subunits of complex I in the bovine, Escherichia coli and T. thermophilus enzymes are diminished (32, 42, 43). Consistently, the number of intersubunit contacts (calculated using the CCP4 program CONTACT with 3.5 Å cutoff) is reduced by ⬃7– 8% in our structures of the reduced enzyme (RND and RN) as compared with the oxidized form (O4 and O2). This likely reflects small adjustments occurring upon the reduction of iron-sulfur clusters, some of which build up to the shifts of helices described above. It cannot be excluded that the extent of these movements might be different in the intact complex. However, our cross-linking studies have shown a similar degree of change in cross-linking patterns upon reduction both for the intact complex and for the isolated hydrophilic domain (32). Surprisingly, movements near the interface with the membrane domain are not accompanied by any significant changes in the central part of the structure, and so they do not appear to be communicated directly by protein interactions. Therefore, it is likely that these shifts of the helices are driven by a change in the redox state of nearby iron-sulfur clusters (N2, N6a, and N6b (8)). The structure of the dithionite-reduced domain (RD) shows similar movements of helices H1 and H2, but interestingly, no shifts in the four-helix bundle of Nqo4 (supplemental Fig. S4B). Dithionite reduces complex I clusters slowly, on a time scale of minutes (44), and so it is likely that in our RD case, after a short exposure of crystal to air during harvesting, cluster N2 remains reduced because of its high potential, whereas other clusters are oxidized. On the other hand, NADH reduces most of complex I clusters, including N6b, on a millisecond time scale (34), and so they will remain reduced in RND and RN JOURNAL OF BIOLOGICAL CHEMISTRY
29777
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
the isoalloxazine ring, in a position that is typical for direct hydride transfer in flavoproteins. This arrangement is consistent with the known stereospecificity of the complex I reaction for the 4B hydrogen of NADH (31, 33). Upon NADH binding, the side chain of conserved 1Glu97 moves away from the flavin, allowing the nicotinamide ring to occupy its position. The carboxyl group of 1Glu97 is involved in several hydrogen bonds to the protein backbone and the side chain of the conserved 180 , which may stabilize the observed conformation. The 1Tyr  C atom of this glutamate is in van der Waal’s contact (3.2 Å) with the C4N atom of the nicotinamide, likely pushing it closer to the N5 atom of FMN and thus promoting hydride transfer from C4N to N5. The distance between these two atoms (3.2 Å) is slightly shorter than average for flavoenzymes (30), consistent with fast rates of hydride transfer in complex I (34). The particular location of C4N, as a site of oxidative attack, relative to the isoalloxazine ring (its projected position falls between the two lines defined by the C4A-C10A and N5-N10 atoms of FMN), is very similar to that seen in other flavoenzymes, consistent with a notion of highly conserved, stereochemically optimized geometry of the active site (30). Upon NADH binding, the isoalloxazine ring shifts ⬃0.5 Å away from neighboring ironsulfur cluster N3. This is consistent with the observed line width narrowing of the EPR signal of flavosemiquinone radical in the presence of NADH, suggested to occur because of diminished dipole-dipole interaction with cluster N3 (35). This movement of the flavin is likely due to an overall shift of the structure around the glycine-rich loop away from the NADH. Thus, this structure shows how an unusual Rossmann fold domain of complex I, which evolved to bind both flavin and NADH by the addition of an extra glycine-rich loop (8), incorporates familiar elements of nucleotide binding in a different structural environment. It explains the kinetic properties of NADH oxidation by complex I, such as high affinity, stereospecificity, substrate specificity, and fast hydride transfer rates. All of the residues involved in interactions with the cofactors and most residues lining the binding cavity are very well conserved (supplemental Fig. S3), suggesting that the mechanism of dehydrogenation by complex I is conserved throughout the species from bacteria to humans. During turnover in vivo (kcat ⫽ ⬃200 s⫺1), most of the iron-sulfur clusters are reduced (36), and the likely rate-limiting steps are NAD⫹ release (34) and/or quinone/quinol binding/release. Implications for Generation of ROS by Complex I—Flavin is now generally regarded as a main source of ROS in complex I (4, 37–39), although cluster N2 (40) and (semi)quinone (41) are also being discussed. In the absence of NADH, the isoalloxazine ring is exposed to the solvent (Fig. 1C), and so if the flavin is reduced, ROS formation is inevitable. However, when a nucleotide is bound, the flavin is shielded from the solvent (Fig. 1D). Therefore, during normal turnover, ROS formation would be minimal (⬃1% of electrons escape to oxygen (37, 39)), because electron transfer from flavin down the redox chain is faster than NAD⫹ release (34). However, under conditions when reoxidation of flavin is slow (the presence of complex I inhibitors, a highly reduced quinone pool, or reverse electron transport), the rate of ROS production increases (39), consistent with our structural interpretation. Therefore, an attractive candidate for
Mechanism of Respiratory Complex I
29778 JOURNAL OF BIOLOGICAL CHEMISTRY
VOLUME 284 • NUMBER 43 • OCTOBER 23, 2009
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
and neighboring clusters in the NADH aerobic state (RA1 and RA2) is similar to that of the dithionite state (RD). Consistently, in the NADH aerobic state helices H1/H2 undergo similar shifts as in the RD (and RND) structures, whereas the four-helix bundle of Nqo4 shows limited or no movement (as in the RD structure). This suggests that the H1/H2 movement is mainly driven by the reduction of cluster N2, to which these helices are connected either directly (H2) or through a short -strand (H1). EPR studies of the intact complex suggest that although in the NADHreduced enzyme cluster N6b is likely to be reduced in addition to N2, N6a may not be (34, 46, 47). Therefore it is likely that the shift of the four-helix bundle, which occurs in this case (RND structure), is driven by the reduction of cluster N6b or N6b/N2 as a pair. Contribution from N6a (or even more distant clusters) cannot be excluded. However, the large distance from cluster N6a to the preceding N5 (8), combined with the possibility of a low potential of N6a (47), may allow only very transient reduction of N6a under steady state conditions, so that the N6b/N2 cluster pair works with electrons delivered effectively in pairs from NADH (and donated in pairs to quinone). Cluster N6b is hydrogen-bonded to semi-conserved 9His41, whereas N6a does not FIGURE 2. Conformational changes upon reduction by NADH. A and B, the structure of the NADH-reduced have such a partner and so may have domain (RND) was superposed with that of oxidized domain (O2) and colored according to C␣ deviations between the structures - blue to green to yellow as deviation increases. Higher deviations near the nucleotide- a lower potential. The conclusion binding site (top) and at the interface with the membrane domain (bottom) are evident. NADH is shown as that the reduction of iron-sulfur spheres in magenta. C, close-up of the interface with the membrane domain (structures RND and O4 were clusters drives conformational superposed). The approximate (expected) border of the membrane domain is shaded in blue. Helices H1 and H2 from Nqo6 subunit and the four-helix bundle from Nqo4 are indicated. Clusters N2 and N6b are shown as changes, rather than nucleotide spheres. The arrows indicate movements of helix H1 and the four-helix bundle upon reduction. binding per se, is consistent with our cross-linking data (32). In the intact structures, despite a similar (short) exposure of crystals to air. complex, helix H1 is likely to extend toward the antiporter-like After a prolonged exposure, as in the RA1 and RA2 structures, subunits Nqo12–14 (4). The long and rigid four-helix bundle, the density for bound nucleotide becomes weaker than in the with its apparent concerted movement into the membrane anaerobic state (RND and RN). It is likely that in the aerobic domain, and the amphipathic helix H1 moving like a “lever” state, as NADH is consumed to produce ROS (37) and NADH/ along the surface of the membrane domain (Fig. 2C), together NAD⫹ ratio in the solution drops, most of the redox chain in may be able to drive rearrangements of transmembrane helices, the complex is oxidized, whereas cluster N2 remains reduced. leading to proton translocation. Environment of Cluster N2—The terminal cluster N2 is the Indeed, EPR experiments on the Neurospora crassa enzyme have shown that after reduction of the complex by low amounts electron donor to quinone (2, 8, 46). It is coordinated, uniquely, of NADH, cluster N2 remained reduced after exposure to air, by two consecutive cysteines of Nqo6 subunit (“tandem cyswhereas other clusters were reoxidized (45). Similar N2 coor- teine motif”). This leads to an unfavorable, strained geometry of dination patterns (see below) suggest that the redox state of N2 the side chains of the two cysteines. Recently, a similar cluster
Mechanism of Respiratory Complex I
OCTOBER 23, 2009 • VOLUME 284 • NUMBER 43
JOURNAL OF BIOLOGICAL CHEMISTRY
29779
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
reduced domain (structures RND and RN), the density connecting 46 to the cluster disappears 6Cys completely (Fig. 3, C and D). It cannot be excluded that 6Cys46 remains a coordinating ligand but with much increased flexibility. However, on the basis of the electron density (the backbone density for the cysteine is too far from the potential S␥ position; Fig. 3, C and D), it is more likely that the iron atom of the reduced cluster is now coordinated by another ligand, such as an hydroxide or a water molecule, as seen in aconitase (49). In principle, solvent-coordinated iron atom might be protonated directly. However, importantly, such a disconnection is likely to lead to protonation of the free 6Cys46. Strikingly, in the structure of dithionite-reduced domain (RD), 45 becomes disconnected in6Cys stead, whereas 6Cys46 retains good connecting density (Fig. 3E). In both structures of aerobically harvested NADH-reduced enzyme (RA1 and RA2), the density connecting 6Cys46 to the cluster also comes back fully (i.e. is similar to that in oxidized enzyme), whereas density connecting 6Cys45 to the cluster also disappears (Fig. 3F). Therefore, RA1, RA2, and RD structures are likely to represent a similar redox state, consistent with crystal harvesting conditions, as discussed above. This suggests that when cluster N2 is reduced, one of the tandem cysFIGURE 3. Environment of cluster N2. A and B, oxidized domain (structure O2). Two copies of the molecule in the ASU are shown. C and D, NADH-reduced domain (structure RND). Two copies of the molecule in the ASU are teines disconnects, the choice being shown. E, dithionite-reduced domain (structure RD). F, NADH-reduced domain exposed to air for ⬃30 min dependent on the redox state of (structure RA1). Cluster ligands and residues discussed in the text are shown as sticks. The backbones of subunits are colored green for Nqo4 and red for Nqo6. Electron density is from a A-weighted 2Fo ⫺ Fc maps neighboring clusters. When only N2 contoured at 1. The cluster is shown as spheres of 0.3 van der Waal’s radius. The prefixes before the residue is reduced (RD, RA1, and RA2 strucnames indicate the subunit number. tures), 6Cys45 is disconnected, but when N6b is also reduced (RND and coordination was discovered in adenosine 5⬘-phosphosulfate RN structures), 6Cys46 is disconnected instead. The redox state reductase (48). The motif is fully conserved in complex I and so of N6b may be communicated to N2 via conformational is likely to be important for the function. The strained geometry changes that we observe (shift of four-helix bundle in RND is likely to lead to flexibility in this area, which appears to be structure, as compared with RD). Furthermore, the interaction important for the conformational changes that we observe, par- may also be mediated via charged residues. For example, ticularly around helices H1/H2. Another possibility is that such 4Arg84, which is conserved and essential for activity (50), is coordination may allow the cluster to be protonated upon hydrogen-bonded to cluster N2 but is also only 6 Å away from reduction, because its redox potential is pH-dependent (redox- cluster N6b ligand 9Cys101 and may be hydrogen-bonded to the Bohr effect) (46). main chain carbonyl of 9Ile99, located between the two Cys In both of our structures of the oxidized domain, the electron ligands of N6b. density connecting both tandem cysteines to the cluster is very Similarly to 6Cys46 in the NADH anaerobic state, cluster clear (Fig. 3, A and B, and Ref. 8). Remarkably, in the NADH- coordination by 6Cys45 in the dithionite anaerobic (or NADH
Mechanism of Respiratory Complex I
29780 JOURNAL OF BIOLOGICAL CHEMISTRY
FIGURE 4. Bound cations. A, an overview. The three structures (O4, O2, and RND) containing three different types of bound cations were superposed, and the cations are shown as spheres. Oxidized domain structure (O4) is shown in cartoon representation with subunit Nqo1 is in yellow, Nqo2 is in blue, Nqo3 is in salmon, Nqo4 is in green, Nqo5 is in violet, Nqo6 is in red, Nqo9 is in cyan, and Nqo15 is in orange. AG indicates the acidic groove. Ca2⫹ ions are shown in green except the tightly bound Ca2⫹ present in all structures, which is in cyan. Mg2⫹ is in blue, and Mn2⫹ is in magenta except the Mn2⫹ atom bound in the putative ironbinding channel, which is in orange. B, close-up of the tightly bound Ca2⫹ ion with interacting residues indicated. C, close-up of the Mn2⫹ ion bound in the putative iron-binding channel at the interface with frataxin-like subunit Nqo15. Residues interacting with the metal are indicated, including possible weak interaction with carbonyl of 2Lys121. The prefixes before the residue names indicate the subunit number.
activity (4). The fourth Ca2⫹ atom is involved in crystal contacts (as well as some of Mg2⫹ and Mn2⫹ atoms in the other two structures), explaining why high concentrations of divalent cations are necessary for crystallization. In the RND structure (crystals grown in 0.1 M MgCl2), we could resolve 5–7 bound cations for each molecule, three of which were also found in the acidic groove (Fig. 4A), confirming that this area is involved in interactions with divalent cations (4). Most cations refined well as Mg2⫹ atoms, except for one tightly bound (almost buried under the surface) at the interface VOLUME 284 • NUMBER 43 • OCTOBER 23, 2009
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
aerobic) state is likely to be replaced by a solvent ligand. We have found previously that locally bound nonchelatable metal may be available for the Fenton reaction near cluster N2 (32), consistent with the possibility of cluster ligation by solvent. Under no conditions did we observe a structure with both tandem cysteines disconnected. The similarity of dithionite anaerobic and NADH aerobic states suggests that the latter is not the result of oxidative damage. The disconnection of tandem cysteines is also not likely to result from radiation damage during data collection (supplemental text). Thus, it is likely that during the catalytic cycle 6Cys45 and 46 would become disconnected from the cluster (and so 6Cys protonated) in sequence when cluster N2 is reduced. 6Cys45 is next to invariant and essential for activity 4Tyr87 (50), which points into the quinone-binding cavity, and so one can envisage a 6Cys45 3 4Tyr87 3 (Q-site) proton delivery pathway. Although 4Tyr87 is one of likely secondary proton acceptors, there are other possibilities, such as nearby conserved 6Arg83, which is very well positioned for this role, making apparent hydrogen bonds to the sulfurs of the cluster and 6Cys45. It will be of interest to study the role of this arginine by mutagenesis. Protons used during quinone reduction in complex I are delivered from the cytosol side in bacteria (or matrix in mitochondria). However, the quinone-binding site resembles a funnel pointing into the membrane (8) and may not be directly accessible from the solvent in the intact complex. Interestingly, a densely packed path of conserved ionizable residues connects the cytosol-exposed part of the Nqo4 subunit to the 4Tyr8745 area (supplemental Fig. S5), possibly representing a pro6Cys ton delivery pathway. A similar pathway, with at least four residues conserved between complex I and hydrogenases (4Arg350, 401 , 4His129, and 4Arg279), was suggested for proton trans4Asp fer in the homologous [Ni-Fe]-hydrogenases (51). Other proton delivery pathways to either one of the tandem cysteines cannot be excluded without further mutagenesis data. In the large subunit of [Ni-Fe]-hydrogenases, invariant Glu18 is essential for proton transfer (51) and superposes by structural alignment to 87 in complex I, consistent with the possible role of this 4Tyr tyrosine in proton delivery. Thus, proton transfer toward the Q-site may be controlled, through tandem cysteines, by the redox state of clusters N2/N6b. Such a control is not necessary for simple “chemical” protonation of quinone but would be needed as a part of coupling mechanism (see below). Bound Divalent Cations—Each of the three higher resolution structures described here has been solved in the presence of different divalent cations, Ca2⫹, Mg2⫹, or Mn2⫹, mostly in an effort to improve the resolution. This has also allowed us to identify cations bound to the protein. In the O4 structure we could resolve at least four ions bound to each molecule (Fig. 4A). The presence of 0.1 M CaCl2 in the mother liquor, B factor refinement, the geometry and nature of coordination ligands (carboxyls of Glu and Asp, main chain carbonyl and N⑀ of His) allowed us to assign the bound metal as Ca2⫹. Interestingly, three of four cations are bound in the previously identified acidic groove (Fig. 4A), a negatively charged protein surface between the C-terminal domain of Nqo3 and the rest of the complex (4). We have suggested it as a possible area of interaction with cations, which are known to modulate complex I
Mechanism of Respiratory Complex I
OCTOBER 23, 2009 • VOLUME 284 • NUMBER 43
JOURNAL OF BIOLOGICAL CHEMISTRY
29781
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
in a similar position, which represents one of the previously predicted iron-binding sites (4). Additionally, in the O2 and RND structures, an unexplained electron density might represent a cation coordinated by 3His208, 15His92, and 94 15His , in the vicinity of another previously suggested iron-binding site (4). However, this putative manganese/magnesium atom was not included in the models because of apparently less than full occupancy. Subunit Nqo15 is structurally homologous to frataxins and iron chaperones and represents the first example of frataxin-like protein permanently bound to its partner (4, 16). Frataxin deficiency in humans leads to Friedreich’s ataxia, a severe neurodegenerative disorder (52). However, exactly how frataxin delivers iron (Fe2⫹) to its protein partners is not known, and T. thermophilus complex I may provide such a model. Iron and manganese are neighboring transition metals with similar properties. Often, they can be exchanged in biological reactions, and so Mn2⫹ can be used to FIGURE 5. A model of complex I mechanism, based on structures described here. Iron-sulfur clusters N2 mimic easily oxidized Fe2⫹. The fact and N6b are depicted as oxidized or reduced by empty or filled circles, respectively. C46 and C45 indicate the 2⫹ tandem cysteines from Nqo6 subunit, whereas Q/QH2 indicate quinone/quinol. One of possible proton accep- that Mn was identified as a bound 87 49 tors from C45 is 4Tyr (Y-O) and from C46-6Glu (E-O). H-path indicates proton delivery pathway from the metal in the putative iron-binding cytosol to tandem cysteines, such as that shown in supplemental Fig. 5, which may be the same for both channel suggests that iron is also cysteines (although other pathways are also possible). H1 and 4HB indicate helix H1 from Nqo6 and the fourhelix bundle from Nqo4, respectively, and are shown on dark backgrounds when helices are shifted relative to likely to bind in the same location. oxidized state. The proposed catalytic cycle is described in detail in the text. The metal ion shown in Fig. 4C is not directly coordinated by residues of subunits Nqo3 and Nqo15 (Fig. 4B). B factor refinement sug- of Nqo15 but interacts with 1His350, 2Glu123, and 2Ser68. This gested that it is a heavier atom, with a mass similar to Ca2⫹. This puts the metal within ⬃10 Å of cluster N1a. This cluster might atom is present in all three structures. It appears to stabilize the have unique role as an antioxidant (4, 8), and if iron binds in this Nqo3-Nqo15 interaction, and it is known that calcium stabi- location, it could be used for regeneration of the cluster in case lizes bacterial complex I (32). On the basis of these observa- of iron loss. Overall, the structure is consistent with iron being tions, we tentatively assigned this atom as Ca2⫹ in all structures, delivered by frataxin multimers, because monomers do not on the assumption that it remains tightly bound to the complex bind it with high affinity (52), so that upon delivery iron would throughout purification and crystallization. This remains to be be coordinated by a protein partner. Complex I Mechanism—The two main models of coupling confirmed by further investigation. The only other cation found in the same position in all three structures is coordinated by the between electron transfer and proton translocation in complex carboxyl group of 3Asp302 and the main chain carbonyl of I, conformational and direct coupling, are not mutually exclu274 , also in the acidic groove. It cannot be excluded that this sive. To account for the high proton/electron stoichiometry, it 3Leu is a tightly bound Ca2⫹ as well, although B factor refinement was proposed that a combination of the two mechanisms is possible (8, 11, 53). As shown by our structures, upon reduction was less conclusive in this case. In the O2 structure (crystals grown in 0.1 M MnCl2), six to some rather small but significant and concerted shifts of several eight cations bound per molecule were resolved. Again, three ␣-helices occur at the interface with the membrane domain. cations were found in the acidic groove, whereas one Mn2⫹ ion These movements may be sufficient to drive proton transloca(only in the first copy of the molecule in ASU) was found in the tion by inducing conformation changes in transmembrane helpreviously identified putative iron-binding channel (8) at the ices. Indeed, our cross-linking studies indicated that the interinterface of subunit Nqo15 with the rest of the complex (Fig. face between two hydrophobic subunits is changed upon 4C). A bound Mg2⫹ atom was identified in the RND structure reduction of E. coli complex I (32). Because there are three sim-
Mechanism of Respiratory Complex I
29782 JOURNAL OF BIOLOGICAL CHEMISTRY
coordination of cluster N2 by tandem cysteines is central both for conformational coupling (allowing movements of helices H1/H2) and for direct coupling via cysteine protonation. One of the consequences of our mechanism would be pH dependence of the midpoint redox potential of cluster N2, because its reduction is accompanied by the protonation of 6Cys45 (which is thus predicted to act as a redox-Bohr group). Such pH dependence of N2 potential is indeed observed (46), consistent with our proposal. An alternative candidate for the redox-Bohr group has been discussed, but the evidence is ambiguous (supplemental text). Previously, the redox-Bohr effect has been repeatedly proposed to play some role in the proton translocation mechanism of complex I, and our mechanism shows the likely details of the process. Mutagenesis of the residues involved (tandem cysteines (28), 4Tyr87(50), and 6Glu49(54)) leads to the loss of oxidoreductase and proton-translocation activity, consistent with our model. Elements of the mechanism suggested here may be applicable to other enzymes harboring a similar cluster coordinating motif (48). It appears that in complex I, unfavorable and unstable coordination of the cluster was evolutionarily selected because of the associated high efficiency of proton pumping. The full details of the coupling mechanism can be established by solving the structure of the intact enzyme and by further mutagenesis studies. Acknowledgments—We thank the staffs of Beamlines ID23-1 and ID29 (European Synchrotron Radiation Facility) and I02 (Diamond Light Source) for assistance. REFERENCES 1. Walker, J. E. (1992) Q. Rev. Biophys. 25, 253–324 2. Yagi, T., and Matsuno-Yagi, A. (2003) Biochemistry 42, 2266 –2274 3. Carroll, J., Fearnley, I. M., Skehel, J. M., Shannon, R. J., Hirst, J., and Walker, J. E. (2006) J. Biol. Chem. 281, 32724 –32727 4. Sazanov, L. A. (2007) Biochemistry 46, 2275–2288 5. Schapira, A. H. (1998) Biochim. Biophys. Acta 1364, 261–270 6. Balaban, R. S., Nemoto, S., and Finkel, T. (2005) Cell 120, 483– 495 7. Dawson, T. M., and Dawson, V. L. (2003) Science 302, 819 – 822 8. Sazanov, L. A., and Hinchliffe, P. (2006) Science 311, 1430 –1436 9. Carroll, J., Fearnley, I. M., and Walker, J. E. (2006) Proc. Natl. Acad. Sci. U.S.A. 103, 16170 –16175 10. Dutton, P. L., Moser, C. C., Sled, V. D., Daldal, F., and Ohnishi, T. (1998) Biochim. Biophys. Acta 1364, 245–257 11. Friedrich, T. (2001) J. Bioenerg. Biomembr. 33, 169 –177 12. Brandt, U. (2006) Annu. Rev. Biochem. 75, 69 –92 13. Fearnley, I. M., and Walker, J. E. (1992) Biochim. Biophys. Acta 1140, 105–134 14. Mathiesen, C., and Ha¨gerha¨ll, C. (2003) FEBS Lett. 549, 7–13 15. Baranova, E. A., Morgan, D. J., and Sazanov, L. A. (2007) J. Struct. Biol. 159, 238 –242 16. Hinchliffe, P., Carroll, J., and Sazanov, L. A. (2006) Biochemistry 45, 4413– 4420 17. CCP4 (1994) Acta Crystallogr. Sect. D 50, 760 –763 18. McCoy, A. J., Grosse-Kunstleve, R. W., Adams, P. D., Winn, M. D., Storoni, L. C., and Read, R. J. (2007) J. Appl. Crystallogr. 40, 658 – 674 19. Vinogradov, A. D. (2008) Biochim. Biophys. Acta 1777, 729 –734 20. Strong, M., Sawaya, M. R., Wang, S., Phillips, M., Cascio, D., and Eisenberg, D. (2006) Proc. Natl. Acad. Sci. U.S.A. 103, 8060 – 8065 21. La Fortelle, E. D., and Bricogne, G. (1997) Methods Enzymol. 276, 472– 494 22. Brunger, A. T. (2007) Nat. Prot. 2, 2728 –2733
VOLUME 284 • NUMBER 43 • OCTOBER 23, 2009
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
ilar antiporter-like subunits in the membrane domain, it is possible that such conformational coupling would lead to the translocation of three protons (one/each subunit) in one catalytic cycle. The remaining fourth proton (accounting for the 4H⫹/2e⫺ stoichiometry) could be translocated near the Q-site with the participation of cluster N2/tandem cysteines. A putative model of such a mechanism, based on our structures, is shown in Fig. 5. As discussed above, the main determinant appears to be the redox state of clusters N2 and N6b (although contribution from more distant clusters cannot be excluded). Starting from the oxidized state 1 (O4 and O2 structures), the arrival of the first electron from NADH via the redox chain to cluster N2 will result in disconnection and protonation of 6Cys45. In this state, state 2 (dithionite anaerobic/NADH aerobic structures, RD, RA1, and RA2), helices H1/H2 of Nqo6 are shifted, and the proton is delivered to 6Cys45. Although we illustrate the figure with two proton delivery pathways for two tandem cysteines, there may be only one pathway, such as that suggested in supplemental Fig. S5, bifurcating near the cysteines. The reduction of cluster N6b by the second electron will result in state 3, where 6Cys45 is reconnected to the cluster N2 (and so deprotonated), whereas 6Cys46 is disconnected and so protonated (NADH anaerobic structures, RND and RN). The four-helix bundle of Nqo4 is shifted; the proton leaving 6Cys45 is transferred to one of the residues nearby, such as 4Tyr87 or 83 (Y-OH). Then, after cluster N2 donates the first electron 6Arg to quinone, creating a semiquinone radical, it is quickly reduced again by cluster N6b (because of the high potential of N2), and state 4 is created. The redox state of the clusters here is similar to state 2, and so 6Cys46 is reconnected to the cluster, whereas 45 is disconnected and so protonated again. The proton 6Cys leaving 6Cys46 is transferred to one of ionizable groups nearby, such as 6Glu49 (E-OH). This glutamate is essential for activity of E. coli complex I, and as discussed previously (54), the rotation of the side chain can bring it to within hydrogen-bonding distance of 6Cys46. In the next step, N2 donates an electron to the semiquinone, which is protonated by the Y-OH and E-OH groups, so that quinol is produced. Helices H1/H2 return to their starting positions (because N2 is oxidized), and this step is linked to the translocation of three protons through antiporterlike subunits Nqo12–14 via conformational changes. Conformational coupling probably involves the four-helix bundle as well and so is initiated at earlier steps. Its depiction at this step reflects its final stage and the fact that semiquinone radical is strongly stabilized by the proton-motive force (55). Because cluster N2 is likely to be oxidized first, before Y-OH donates its proton to the quinone, the proton released upon 6Cys45 reconnection to the cluster can take a different pathway from that shown in step 3, and so it can be ejected into the periplasm. The exit pathway for this “pumped” proton is likely to be located in the adjacent hydrophobic subunits. Finally, quinol is released from the binding site, quinone binds instead, and the cycle restarts. The sum of the reactions depicted in Fig. 5 will be the transfer of two electrons from NADH and two protons from the cytoplasm to the quinone, coupled to the translocation of four protons across the membrane, consistent with known stoichiometry. In our “combined” complex I mechanism, the unique
Mechanism of Respiratory Complex I
OCTOBER 23, 2009 • VOLUME 284 • NUMBER 43
1757, 553–561 40. Fato, R., Bergamini, C., Leoni, S., Strocchi, P., and Lenaz, G. (2008) Neurochem. Res. 33, 2487–2501 41. Lambert, A. J., and Brand, M. D. (2004) J. Biol. Chem. 279, 39414 –39420 42. Belogrudov, G., and Hatefi, Y. (1994) Biochemistry 33, 4571– 4576 43. Mamedova, A. A., Holt, P. J., Carroll, J., and Sazanov, L. A. (2004) J. Biol. Chem. 279, 23830 –23836 44. Kowal, A. T., Morningstar, J. E., Johnson, M. K., Ramsay, R. R., and Singer, T. P. (1986) J. Biol. Chem. 261, 9239 –9245 45. Rasmussen, T., Scheide, D., Brors, B., Kintscher, L., Weiss, H., and Friedrich, T. (2001) Biochemistry 40, 6124 – 6131 46. Ohnishi, T. (1998) Biochim. Biophys. Acta 1364, 186 –206 47. Euro, L., Bloch, D. A., Wikstro¨m, M., Verkhovsky, M. I., and Verkhovskaya, M. (2008) Biochemistry 47, 3185–3193 48. Chartron, J., Carroll, K. S., Shiau, C., Gao, H., Leary, J. A., Bertozzi, C. R., and Stout, C. D. (2006) J. Mol. Biol. 364, 152–169 49. Beinert, H., Kennedy, M. C., and Stout, C. D. (1996) Chem. Rev. 96, 2335–2374 50. Kashani-Poor, N., Zwicker, K., Kerscher, S., and Brandt, U. (2001) J. Biol. Chem. 276, 24082–24087 51. Fontecilla-Camps, J. C., Volbeda, A., Cavazza, C., and Nicolet, Y. (2007) Chem. Rev. 107, 4273– 4303 52. Karlberg, T., Schagerlo¨f, U., Gakh, O., Park, S., Ryde, U., Lindahl, M., Leath, K., Garman, E., Isaya, G., and Al-Karadaghi, S. (2006) Structure 14, 1535–1546 53. Sazanov, L. A., and Walker, J. E. (2000) J. Mol. Biol. 302, 455– 464 54. Flemming, D., Hellwig, P., Lepper, S., Kloer, D. P., and Friedrich, T. (2006) J. Biol. Chem. 281, 24781–24789 55. Yano, T., Dunham, W. R., and Ohnishi, T. (2005) Biochemistry 44, 1744 –1754
JOURNAL OF BIOLOGICAL CHEMISTRY
29783
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
23. Grgic, L., Zwicker, K., Kashani-Poor, N., Kerscher, S., and Brandt, U. (2004) J. Biol. Chem. 279, 21193–21199 24. Jones, T. A., and Kjeldgaard, M. (1997) Methods Enzymol. 277, 173–208 25. Adams, P. D., Grosse-Kunstleve, R. W., Hung, L. W., Ioerger, T. R., McCoy, A. J., Moriarty, N. W., Read, R. J., Sacchettini, J. C., Sauter, N. K., and Terwilliger, T. C. (2002) Acta Crystallogr. Sect. D 58, 1948 –1954 26. Morris, A. L., MacArthur, M. W., Hutchinson, E. G., and Thornton, J. M. (1992) Proteins 12, 345–364 27. Davis, I. W., Leaver-Fay, A., Chen, V. B., Block, J. N., Kapral, G. J., Wang, X., Murray, L. W., Arendall, W. B., 3rd, Snoeyink, J., Richardson, J. S., and Richardson, D. C. (2007) Nucleic Acids Res. 35, W375–383 28. Flemming, D., Schlitt, A., Spehr, V., Bischof, T., and Friedrich, T. (2003) J. Biol. Chem. 278, 47602– 47609 29. Lesk, A. M. (1995) Curr. Opin. Struct. Biol. 5, 775–783 30. Fraaije, M. W., and Mattevi, A. (2000) Trends Biochem. Sci. 25, 126 –132 31. Vinogradov, A. D. (1998) Biochim. Biophys. Acta 1364, 169 –185 32. Berrisford, J. M., Thompson, C. J., and Sazanov, L. A. (2008) Biochemistry 47, 10262–10270 33. Ernster, L., Hoberman, H. D., Howard, R. L., King, T. E., Lee, C. P., Mackler, B., and Sottocasa, G. (1965) Nature 207, 940 –941 34. Verkhovskaya, M. L., Belevich, N., Euro, L., Wikstro¨m, M., and Verkhovsky, M. I. (2008) Proc. Natl. Acad. Sci. U.S.A. 105, 3763–3767 35. Sled, V. D., Rudnitzky, N. I., Hatefi, Y., and Ohnishi, T. (1994) Biochemistry 33, 10069 –10075 36. Krishnamoorthy, G., and Hinkle, P. C. (1988) J. Biol. Chem. 263, 17566 –17575 37. Kussmaul, L., and Hirst, J. (2006) Proc. Natl. Acad. Sci. U.S.A. 103, 7607–7612 38. Kudin, A. P., Bimpong-Buta, N. Y., Vielhaber, S., Elger, C. E., and Kunz, W. S. (2004) J. Biol. Chem. 279, 4127– 4135 39. Grivennikova, V. G., and Vinogradov, A. D. (2006) Biochim. Biophys. Acta
SUPPLEMENTARY TEXT The nature of bound nucleotide. Crystals with bound nucleotide (datasets RND and RN) were grown anaerobically in the presence of excess NADH (3 mM). It is likely that after initial NADH binding and hydride transfer to FMN, the electrons will be passed on to the Fe-S clusters, including the high-potential terminal cluster N2, and the flavin will be re-oxidized. Due to the low affinity of NAD+ to the oxidized enzyme, NAD+ will leave the binding site and another NADH molecule will bind, donating further two electrons into the Fe-S chain. The process will be repeated until the redox chain cannot accept any more electrons. The flavin will then remain reduced and the nucleotide species found in the crystals is likely to be NADH due to its large excess over NAD+. Even if the redox chain is fully reduced using ~5 NADH molecules (10 electrons for 8 Fe-S clusters and FMN) per enzyme (which is present at about 10 μM concentration during crystallization), this will result in about 60-fold excess of NADH over NAD+ at equilibrium. Therefore, the bound nucleotide is modeled in the structure as NADH. Often (but not always) the nicotinamide ring is puckered in NADH bound to the protein (1); however, this cannot be discerned at our current resolution. Flavin environment. Instead of a positively charged entity as in other flavoenzymes (2), the N1C2=O2 locus of FMN is H-bonded to Nδ atom of 1Asn220, consistent with a low redox potential of flavin in complex I (~ -350 mV vs. ~ -220 mV of free flavin (3,4)). Despite an overall positive charge in the NADH binding cavity (5), there are several conserved negatively charged residues in the vicinity of FMN, some of which may also be responsible for its low potential. There is no ionizable residue (such as His) in the vicinity of N1 atom of FMN, therefore only bulk solvent appears to be a likely source of protons when and if FMNH- to FMNH2 protonation happens. Are the shifts of helices upon reduction significant? Relative to 3.1 Å resolution of O2/RND structures, 1 Å shifts may appear to be small. However, they are far larger than cross-validated coordinate error, which is 0.55 Å for RND structure and 0.57 Å for O2 structure, as determined by the refinement program (Phenix). Furthermore, although such errors could result in random XYZ shifts of individual atoms, they cannot lead to a concerted shift of an entire α-helix, consisting of hundreds of atoms. Therefore, even for a single structure such a shift can be determined unambiguously at our resolution. However, the shifts of these helices were reproducible – the two structures of NADH-reduced enzyme (RND and RN) show similar differences from the two oxidized forms (O4 and O2) (Supplementary Fig. 4A). It could also be argued that crystal contacts, in which helix H1 is involved, might play a role in observed changes. However, the two oxidized states do not show such shifts despite being in different crystal forms with
1
different crystal contacts near helix H1. On the other hand, the two structures of the reduced (RND) and oxidized (O2) enzyme in the same crystal form with the same crystal contacts do show such shifts (Supplementary Fig. 4A). Thus, we have observed similar helical shifts in several independently determined structures (and for several copies of the molecule in each asymmetric unit), further confirming the reliability of the results. The disconnection of tandem cysteines is not caused by radiation damage during data collection. It may be argued that radiation damage and reduction of the Fe-S cluster N2 by X-rays might have led to disconnection. However, most of our datasets were collected with similar strategy and similar exposure to X-rays, and systematic differences in tandem Cys coordination were observed only between oxidized and reduced domain states. Importantly, the crystal used for the highest resolution dataset for the reduced domain (RND) was used subsequently for data collection at Fe-peak, resulting in good quality experimental phases (Supplementary Fig. 1A). Also, scaling did not indicate any significant radiation damage in either dataset. Therefore, the first native RND dataset was clearly well below damage or “overreduction” limits. However, in the case of the O2 crystal there was some radiation damage, as mentioned above, and so the O2 crystal clearly was subject to a higher radiation dose than RND. Nevertheless, Cys disconnection was observed only in the RND structure, suggesting that disconnection is not caused by Xray “over-reduction”. Furthermore, X-ray damage may lead to the disappearance of some side chains, however it cannot result in the concerted movement of the protein backbone, since data collection is performed with frozen crystals, where viscosity is similar to that of glass. In our case, the backbone around tandem cysteines is shifted away from cluster N2 in the reduced state (see electron density in Fig. 3C and D), which must have happened before crystals were frozen. Since such a shift is allowed only by Cys disconnection, it further confirms that it is a genuine observation and is not caused by X-ray damage. Redox-Bohr group in complex I. Our mechanism suggests that 6Cys45 is a likely redox-Bohr group, conferring pH dependence on midpoint redox potential of cluster N2. One publication (6) might argue for another residue (His226 of the 49 kDa subunit in Yarrowia lipolytica) as a redox-Bohr group. However, there are some technical problems and alternative explanations, which make that conclusion questionable. In particular, redox titrations were performed not with purified enzyme, but with the mitochondrial membranes. However, in the mutant involving the redox-Bohr group proposed by authors (H226M), N2 signal is hardly visible in the EPR spectra of membranes (as authors themselves acknowledge in the previous publication (7)). Furthermore, even if the titration of such a small signal was accurate, the pH dependence was studied only in the pH range of 6 to 8, due to limited stability of the mutant. The histidine in question (4His169 in Thermus) is within hydrogen-bonding distance not only to cluster N2 (which would
2
explain observed negative shift of N2 redox potential upon its mutation to methionine), but also to 6Cys45 (and thus a H226M mutation could change the pK of this cysteine to the value outside the studied pH range). So, the pH dependence could in fact have been preserved in H226M mutant (just shifted outside the studied range). This would also be consistent with the preserved catalytic activity in the mutant. Therefore, there is no clear evidence to suggest that the redox-Bohr group is 4His169, rather than 6Cys45, which we propose.
3
Dataset / Structure
O4 (oxidised)
O2 (oxidised)
RND (reduced)
Crystallisation conditions
Ca
NAD+, Mn
Mol/ASU, Space group
4, P21
Beamline, wavelength (Å) Resolution range (Å)
RD (reduced)
RA1 (reduced)
Dithionite, Mg
2, P21
NADH, dithionite, NADH, Mn Mg 2, P21 2, P21
4, P1
NADH, Ca NADH, Mg Frozen aerobically Frozen aerobically 4, P21 4, P21
ID29 0.97607 30-3.1 (3.27-3.1)
ID29 ID29 ID29 0.97654 1.73372 0.97650 30-3.1 30-4.0 30-3.7 (3.9-3.7)
I02 0.97960 30-3.7 (3.9-3.7)
ID23-1 0.97950 30-4.2 (4.43-4.2)
ID29 1.00000 30-4.2 (4.43-4.2)
(3.27-3.1)
(4.22-4.0)
198977 (29429)
134538 (19991) 53359 (8035) 90.1 (93.0) 2.5 (2.5)
176390 (26241)
220779 (32522)
236174 (344920)
239105 (34009)
61699 (9273)
129721 (19268)
79492 (116340
84562 (12543)
81.5 (83.7)
84.5 (85.8)
80.2 (80.5)
84.3 (85.8)
1.9 (1.9)
293589 (43457) 111667 (16621) 88.2 (90.2) 2.6 (2.6)
2.9 (2.8)
1.7 (1.7)
3.0 (3.0)
2.8 (2.7)
5.9 (1.2)
7.4 (1.4) 5.2 (1.2) 4.1 (1.2)
4.2 (1.2)
5.0 (1.6)
4.5 (1.5)
13.8 (53.8)
18.6 (61.0)
30.1 (72.5)
33.8 (77.0)
23.5 (33.9)
15.1 20.5 (58.6) (66.5) 23.1 (35.7)
27.5 (60.5)
Rcryst (%)
ID29 0.97564 30-3.15, anisosc. (3.32-3.15) 690795 (93222) - before anisosc. 195404 (14803) -after anisosc. 82.5 (44.0) -after anisosc. 3.2 (3.0) - before anisosc. 6.0 (1.2) - before anisosc. 18.0 (75.9) - before anisosc. 23.7 (36.3)
29.6 (41.2)
27.3 (44.4)
27.0 (38.0)
27.5 (35.7)
Rfree (%)
27.0 (37.8)
28.3 (35.7)
26.2 (39.7)
33.7 (41.9)
30.2 (44.2)
31.8 (40.8)
29.4 (37.4)
Average B factor (Å2)
72.8
64.9
62.2
87.6
101.5
118.1
102.7
Ramachandran plot (%) core/allow/gener/disall PDB accession no.
80.4/18.9/0.6/0.0
80.4/19.0/0.6/0.0
81.6/17.7/0.6/0.1
76.1/23.0/1.0/0.0
77.7/21.7/0.6/0.0
78.9/20.4/0.6/0.1
78.2/20.9/0.9/0.0
3IAS
3I9V
3IAM
Observed reflections Unique reflections Completeness (%) Redundancy Rmerge (%)
102583 (15411) 82.0 (84.4)
RN (reduced)
Supplementary Table 1. Data collection and refinement statistics. The legend is on the following page.
4
RA2 (reduced)
Supplementary Table 1 legend. Values in parentheses are for the highest-resolution shell. Rfree sets consist of 2% of reflections. For the O4 dataset, the statistics are shown either before or after anisotropic scaling and truncation. Ramachandran plot values are from the PROCHECK analysis (8) and correspond to the most-favored / additional allowed / generously allowed / disallowed dihedral angles. Crystallisation conditions and unit cell dimensions for the datasets are as follows. O4 - 0.1 M HEPES pH 7.5, 0.4 M NaCl, 0.1 M CaCl2, 8% PEG4000. a = 135.1 Å, b = 266.1 Å, c = 201.7 Å and β = 104.70 O2 – 20 mM NAD+, 0.1 M HEPES pH 7.5, 0.4 M NaCl, 0.1 M MnCl2, 6% PEG4000, 2.5% PEG400. a = 109.0 Å, b = 150.5 Å, c = 216.7 Å and β = 92.20 RND – 3 mM NADH, 5 mM sodium dithionite, 0.1 M HEPES pH 7.5, 0.05 M NaCl, 0.1 M MgCl2, 7% PEG4000. a = 109.0 Å, b = 151.0 Å, c = 216.6 Å and β = 93.10 RN - 3 mM NADH, 0.1M HEPES pH 7.5, 0.1 M NaCl, 0.15 M MnCl2, 4% PEG4000. a = 110.1 Å, b = 151.6 Å, c = 219.0 Å and β = 90.60 RD - 10 mM sodium dithionite, 0.1 M HEPES pH 7.5, 0.1 M MgCl2, 5.5% PEG4000, 2.5% PEG400. a = 123.5 Å, b = 148.8 Å, c = 210.7 Å; α = 98.60, β = 90.80 and γ = 104.40. RA1 - 3 mM NADH, 0.1M HEPES pH 7.5, 0.1 M NaCl, 0.1M CaCl2, 8% PEG4000. a = 135.6 Å, b = 265.5 Å, c = 201.1 Å and β = 104.40 RA2 - 2 mM NADH, 0.1M HEPES pH 7.5, 0.1 M MgCl2, 7% PEG4000. a = 135.4 Å, b = 266.8 Å, c = 201.1 Å and β = 104.80
5
Supplementary Figure 1. Experimental (A) and calculated (B) electron density around bound nucleotide. The RND structure is viewed from the solvent towards nucleotide binding site, in a similar orientation to Fig.1. NADH is shown with carbon in salmon. In A, the map was obtained by phasing in SHARP using Fe-peak and native data, followed by density modification with solvent flipping in CNS. In B, a final σA-weighted 2Fo-Fc map is shown. Both maps are contoured at 1σ level.
6
Supplementary Figure 2. Interactions between NADH and the protein. The figure was plotted using structure RND (Supplementary Table 1) and program LIGPLOT v 4.4.2 (9). Some weak hydrogen bonds were added on the basis of MOLPROBITY analysis (10). NADH is residue NAI with carbon in purple. “(1)” indicates that residues are from subunit Nqo1.
7
Supplementary Figure 3. Conserved residues in the NADH-binding cavity. Subunit Nqo1 is shown with residues colored blue to red according to the degree of conservation (high to low). High conservation in the vicinity of NADH (shown as sticks with carbon in salmon) and FMN (carbon in yellow) is evident. The figure was produced using structure RND (Supplementary Table 1) and the ConSurf server (11) with default parameters.
8
Supplementary Figure 4. Comparison of the structures of the oxidized and reduced hydrophilic domain. The interface with the membrane domain is shown. (A) Superposition of the structures of the oxidized (O4, in green and O2, in cyan) and NADH-reduced (RND, in orange and RN, in yellow) domain. Similar shifts of helices H1/H2 and the 4-helix bundle (4HB) are evident in the NADH-reduced states as compared to oxidized. (B) Superposition of the structures of the oxidized (O4, in green) and dithionite-reduced (RD, in yellow) domain. In this case, only helices H1/H2 are shifted upon reduction.
9
Supplementary Figure 5. Putative proton delivery pathway from the cytosol to the cluster N2 area. Nqo4 residues, proposed to be involved, are shown as sticks with carbon in green. They are, roughly starting from cluster N2 area, H169, H170, D86, R350, D401, H129, R279, H89, R125, E122, R249, Y257, Y254, Y260, R296 (conserved residues are in bold). 4Tyr87 and tandem cysteines are shown with carbon in yellow. Subunits are colored as in Fig. 4. Clusters N2 and N6b are shown as spheres. Interface with the membrane domain is at the bottom of the figure.
10
References
1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11.
Meijers, R., Morris, R. J., Adolph, H. W., Merli, A., Lamzin, V. S., and Cedergren-Zeppezauer, E. S. (2001) J. Biol. Chem. 276, 9316-9321 Fraaije, M. W., and Mattevi, A. (2000) Trends Biochem. Sci. 25, 126-132 Euro, L., Bloch, D. A., Wikstrom, M., Verkhovsky, M. I., and Verkhovskaya, M. (2008) Biochemistry 47, 3185-3193 Sled, V. D., Rudnitzky, N. I., Hatefi, Y., and Ohnishi, T. (1994) Biochemistry 33, 10069-10075 Sazanov, L. A. (2007) Biochemistry 46, 2275-2288 Zwicker, K., Galkin, A., Drose, S., Grgic, L., Kerscher, S., and Brandt, U. (2006) J. Biol. Chem. 281, 23013-23017 Grgic, L., Zwicker, K., Kashani-Poor, N., Kerscher, S., and Brandt, U. (2004) J. Biol. Chem. 279, 21193-21199 Morris, A. L., MacArthur, M. W., Hutchinson, E. G., and Thornton, J. M. (1992) Proteins 12, 345-364 Wallace, A. C., Laskowski, R. A., and Thornton, J. M. (1995) Protein Eng. 8, 127-134 Davis, I. W., Leaver-Fay, A., Chen, V. B., Block, J. N., Kapral, G. J., Wang, X., Murray, L. W., Arendall, W. B., 3rd, Snoeyink, J., Richardson, J. S., and Richardson, D. C. (2007) Nucleic Acids Res. 35, W375-383 Landau, M., Mayrose, I., Rosenberg, Y., Glaser, F., Martz, E., Pupko, T., and Ben-Tal, N. (2005) Nucl. Acids Res. 33, W299-302
11
Metabolism and Bioenergetics: Structural Basis for the Mechanism of Respiratory Complex I John M. Berrisford and Leonid A. Sazanov J. Biol. Chem. 2009, 284:29773-29783. doi: 10.1074/jbc.M109.032144 originally published online July 27, 2009
Find articles, minireviews, Reflections and Classics on similar topics on the JBC Affinity Sites. Alerts: • When this article is cited • When a correction for this article is posted Click here to choose from all of JBC's e-mail alerts
Supplemental material: http://www.jbc.org/content/suppl/2009/07/27/M109.032144.DC1.html This article cites 55 references, 17 of which can be accessed free at http://www.jbc.org/content/284/43/29773.full.html#ref-list-1
Downloaded from http://www.jbc.org/ by guest on October 29, 2015
Access the most updated version of this article at doi: 10.1074/jbc.M109.032144