Hypoxia and intra-complex genetic suppressors rescue complex I mutants by a shared mechanism.

The electron transport chain (ETC) of mitochondria, bacteria, and archaea couples electron flow to proton pumping and is adapted to diverse oxygen environments. Remarkably, in mice, neurological disease due to ETC complex I dysfunction is rescued by hypoxia through unknown mechanisms. Here, we show that hypoxia rescue and hyperoxia sensitivity of complex I deficiency are evolutionarily conserved to C. elegans and are specific to mutants that compromise the electron-conducting matrix arm. We show that hypoxia rescue does not involve the hypoxia-inducible factor pathway or attenuation of reactive oxygen species. To discover the mechanism, we use C. elegans genetic screens to identify suppressor mutations in the complex I accessory subunit NDUFA6/nuo-3 that phenocopy hypoxia rescue. We show that NDUFA6/nuo-3(G60D) or hypoxia directly restores complex I forward activity, with downstream rescue of ETC flux and, in some cases, complex I levels. Additional screens identify residues within the ubiquinone binding pocket as being required for the rescue by NDUFA6/nuo-3(G60D) or hypoxia. This reveals oxygen-sensitive coupling between an accessory subunit and the quinone binding pocket of complex I that can restore forward activity in the same manner as hypoxia.

Genetic Screen for Cell Fitness in High or Low Oxygen Highlights Mitochondrial and Lipid Metabolism.

Human cells are able to sense and adapt to variations in oxygen levels. Historically, much research in this field has focused on hypoxia-inducible factor (HIF) signaling and reactive oxygen species (ROS). Here, we perform genome-wide CRISPR growth screens at 21%, 5%, and 1% oxygen to systematically identify gene knockouts with relative fitness defects in high oxygen (213 genes) or low oxygen (109 genes), most without known connection to HIF or ROS. Knockouts of many mitochondrial pathways thought to be essential, including complex I and enzymes in Fe-S biosynthesis, grow relatively well at low oxygen and thus are buffered by hypoxia. In contrast, in certain cell types, knockout of lipid biosynthetic and peroxisomal genes causes fitness defects only in low oxygen. Our resource nominates genetic diseases whose severity may be modulated by oxygen and links hundreds of genes to oxygen homeostasis.

Metabolic regulation of species-specific developmental rates.

Animals display substantial inter-species variation in the rate of embryonic development despite a broad conservation of the overall sequence of developmental events. Differences in biochemical reaction rates, including the rates of protein production and degradation, are thought to be responsible for species-specific rates of development1-3. However, the cause of differential biochemical reaction rates between species remains unknown. Here, using pluripotent stem cells, we have established an in vitro system that recapitulates the twofold difference in developmental rate between mouse and human embryos. This system provides a quantitative measure of developmental speed as revealed by the period of the segmentation clock, a molecular oscillator associated with the rhythmic production of vertebral precursors. Using this system, we show that mass-specific metabolic rates scale with the developmental rate and are therefore higher in mouse cells than in human cells. Reducing these metabolic rates by inhibiting the electron transport chain slowed down the segmentation clock by impairing the cellular NAD+/NADH redox balance and, further downstream, lowering the global rate of protein synthesis. Conversely, increasing the NAD+/NADH ratio in human cells by overexpression of the Lactobacillus brevis NADH oxidase LbNOX increased the translation rate and accelerated the segmentation clock. These findings represent a starting point for the manipulation of developmental rate, with multiple translational applications including accelerating the differentiation of human pluripotent stem cells for disease modelling and cell-based therapies.

Loss of LUC7L2 and U1 snRNP subunits shifts energy metabolism from glycolysis to OXPHOS.

Oxidative phosphorylation (OXPHOS) and glycolysis are the two major pathways for ATP production. The reliance on each varies across tissues and cell states, and can influence susceptibility to disease. At present, the full set of molecular mechanisms governing the relative expression and balance of these two pathways is unknown. Here, we focus on genes whose loss leads to an increase in OXPHOS activity. Unexpectedly, this class of genes is enriched for components of the pre-mRNA splicing machinery, in particular for subunits of the U1 snRNP. Among them, we show that LUC7L2 represses OXPHOS and promotes glycolysis by multiple mechanisms, including (1) splicing of the glycolytic enzyme PFKM to suppress glycogen synthesis, (2) splicing of the cystine/glutamate antiporter SLC7A11 (xCT) to suppress glutamate oxidation, and (3) secondary repression of mitochondrial respiratory supercomplex formation. Our results connect LUC7L2 expression and, more generally, the U1 snRNP to cellular energy metabolism.

Hepatic NADH reductive stress underlies common variation in metabolic traits.

The cellular NADH/NAD+ ratio is fundamental to biochemistry, but the extent to which it reflects versus drives metabolic physiology in vivo is poorly understood. Here we report the in vivo application of Lactobacillus brevis (Lb)NOX1, a bacterial water-forming NADH oxidase, to assess the metabolic consequences of directly lowering the hepatic cytosolic NADH/NAD+ ratio in mice. By combining this genetic tool with metabolomics, we identify circulating α-hydroxybutyrate levels as a robust marker of an elevated hepatic cytosolic NADH/NAD+ ratio, also known as reductive stress. In humans, elevations in circulating α-hydroxybutyrate levels have previously been associated with impaired glucose tolerance2, insulin resistance3 and mitochondrial disease4, and are associated with a common genetic variant in GCKR5, which has previously been associated with many seemingly disparate metabolic traits. Using LbNOX, we demonstrate that NADH reductive stress mediates the effects of GCKR variation on many metabolic traits, including circulating triglyceride levels, glucose tolerance and FGF21 levels. Our work identifies an elevated hepatic NADH/NAD+ ratio as a latent metabolic parameter that is shaped by human genetic variation and contributes causally to key metabolic traits and diseases. Moreover, it underscores the utility of genetic tools such as LbNOX to empower studies of 'causal metabolism'.

Circulating N-lactoyl-amino acids and N-formyl-methionine reflect mitochondrial dysfunction and predict mortality in septic shock.

INTRODUCTION: Sepsis is a highly morbid condition characterized by multi-organ dysfunction resulting from dysregulated inflammation in response to acute infection. Mitochondrial dysfunction may contribute to sepsis pathogenesis, but quantifying mitochondrial dysfunction remains challenging. OBJECTIVE: To assess the extent to which circulating markers of mitochondrial dysfunction are increased in septic shock, and their relationship to severity and mortality. METHODS: We performed both full-scan and targeted (known markers of genetic mitochondrial disease) metabolomics on plasma to determine markers of mitochondrial dysfunction which distinguish subjects with septic shock (n = 42) from cardiogenic shock without infection (n = 19), bacteremia without sepsis (n = 18), and ambulatory controls (n = 19) - the latter three being conditions in which mitochondrial function, proxied by peripheral oxygen consumption, is presumed intact. RESULTS: Nine metabolites were significantly increased in septic shock compared to all three comparator groups. This list includes N-formyl-L-methionine (f-Met), a marker of dysregulated mitochondrial protein translation, and N-lactoyl-phenylalanine (lac-Phe), representative of the N-lactoyl-amino acids (lac-AAs), which are elevated in plasma of patients with monogenic mitochondrial disease. Compared to lactate, the clinical biomarker used to define septic shock, there was greater separation between survivors and non-survivors of septic shock for both f-Met and the lac-AAs measured within 24 h of ICU admission. Additionally, tryptophan was the one metabolite significantly decreased in septic shock compared to all other groups, while its breakdown product kynurenate was one of the 9 significantly increased. CONCLUSION: Future studies which validate the measurement of lac-AAs and f-Met in conjunction with lactate could define a sepsis subtype characterized by mitochondrial dysfunction.

A natural fusion of flavodiiron, rubredoxin, and rubredoxin oxidoreductase domains is a self-sufficient water-forming oxidase of Trichomonas vaginalis.

Microaerophilic pathogens such as Giardia lamblia, Entamoeba histolytica, and Trichomonas vaginalis have robust oxygen consumption systems to detoxify oxygen and maintain intracellular redox balance. This oxygen consumption results from H2O-forming NADH oxidase (NOX) activity of two distinct flavin-containing systems: H2O-forming NOXes and multicomponent flavodiiron proteins (FDPs). Neither system is membrane bound, and both recycle NADH into oxidized NAD+ while simultaneously removing O2 from the local environment. However, little is known about the specific contributions of these systems in T. vaginalis. In this study, we use bioinformatics and biochemical analyses to show that T. vaginalis lacks a NOX-like enzyme and instead harbors three paralogous genes (FDPF1-3), each encoding a natural fusion product between the N-terminal FDP, central rubredoxin (Rb), and C-terminal NADH:Rb oxidoreductase domains. Unlike a "stand-alone" FDP that lacks Rb and oxidoreductase domains, this natural fusion protein with fully populated flavin redox centers directly accepts reducing equivalents of NADH to catalyze the four-electron reduction of oxygen to water within a single polypeptide with an extremely high turnover. Furthermore, using single-particle cryo-EM, we present structural insights into the spatial organization of the FDP core within this multidomain fusion protein. Together, these results contribute to our understanding of systems that allow protozoan parasites to maintain optimal redox balance and survive transient exposure to oxic conditions.

MitoCarta3.0: an updated mitochondrial proteome now with sub-organelle localization and pathway annotations.

The mammalian mitochondrial proteome is under dual genomic control, with 99% of proteins encoded by the nuclear genome and 13 originating from the mitochondrial DNA (mtDNA). We previously developed MitoCarta, a catalogue of over 1000 genes encoding the mammalian mitochondrial proteome. This catalogue was compiled using a Bayesian integration of multiple sequence features and experimental datasets, notably protein mass spectrometry of mitochondria isolated from fourteen murine tissues. Here, we introduce MitoCarta3.0. Beginning with the MitoCarta2.0 inventory, we performed manual review to remove 100 genes and introduce 78 additional genes, arriving at an updated inventory of 1136 human genes. We now include manually curated annotations of sub-mitochondrial localization (matrix, inner membrane, intermembrane space, outer membrane) as well as assignment to 149 hierarchical 'MitoPathways' spanning seven broad functional categories relevant to mitochondria. MitoCarta3.0, including sub-mitochondrial localization and MitoPathway annotations, is freely available at http://www.broadinstitute.org/mitocarta and should serve as a continued community resource for mitochondrial biology and medicine.

Top-down characterization of endogenous protein complexes with native proteomics.

Protein complexes exhibit great diversity in protein membership, post-translational modifications and noncovalent cofactors, enabling them to function as the actuators of many important biological processes. The exposition of these molecular features using current methods lacks either throughput or molecular specificity, ultimately limiting the use of protein complexes as direct analytical targets in a wide range of applications. Here, we apply native proteomics, enabled by a multistage tandem MS approach, to characterize 125 intact endogenous complexes and 217 distinct proteoforms derived from mouse heart and human cancer cell lines in discovery mode. The native conditions preserved soluble protein-protein interactions, high-stoichiometry noncovalent cofactors, covalent modifications to cysteines, and, remarkably, superoxide ligands bound to the metal cofactor of superoxide dismutase 2. These data enable precise compositional analysis of protein complexes as they exist in the cell and demonstrate a new approach that uses MS as a bridge to structural biology.

An informatic framework for decoding protein complexes by top-down mass spectrometry.

Efforts to map the human protein interactome have resulted in information about thousands of multi-protein assemblies housed in public repositories, but the molecular characterization and stoichiometry of their protein subunits remains largely unknown. Here, we report a computational search strategy that supports hierarchical top-down analysis for precise identification and scoring of multi-proteoform complexes by native mass spectrometry.

ß-Glucocerebrosidase Modulators Promote Dimerization of ß-Glucocerebrosidase and Reveal an Allosteric Binding Site.

ß-Glucocerebrosidase (GCase) mutations cause Gaucher's disease and are a high risk factor in Parkinson's disease. The implementation of a small molecule modulator is a strategy to restore proper folding and lysosome delivery of degradation-prone mutant GCase. Here, we present a potent quinazoline modulator, JZ-4109, which stabilizes wild-type and N370S mutant GCase and increases GCase abundance in patient-derived fibroblast cells. We then developed a covalent modification strategy using a lysine targeted inactivator (JZ-5029) for in vitro mechanistic studies. By using native top-down mass spectrometry, we located two potentially covalently modified lysines. We obtained the first crystal structure, at 2.2 Å resolution, of a GCase with a noniminosugar modulator covalently bound, and were able to identify the exact lysine residue modified (Lys346) and reveal an allosteric binding site. GCase dimerization was induced by our modulator binding, which was observed by native mass spectrometry, its crystal structure, and size exclusion chromatography with a multiangle light scattering detector. Finally, the dimer form was confirmed by negative staining transmission electron microscopy studies. Our newly discovered allosteric site and observed GCase dimerization provide a new mechanistic insight into GCase and its noniminosugar modulators and facilitate the rational design of novel GCase modulators for Gaucher's disease and Parkinson's disease.

Mapping Proteoforms and Protein Complexes From King Cobra Venom Using Both Denaturing and Native Top-down Proteomics.

Characterizing whole proteins by top-down proteomics avoids a step of inference encountered in the dominant bottom-up methodology when peptides are assembled computationally into proteins for identification. The direct interrogation of whole proteins and protein complexes from the venom of Ophiophagus hannah (king cobra) provides a sharply clarified view of toxin sequence variation, transit peptide cleavage sites and post-translational modifications (PTMs) likely critical for venom lethality. A tube-gel format for electrophoresis (called GELFrEE) and solution isoelectric focusing were used for protein fractionation prior to LC-MS/MS analysis resulting in 131 protein identifications (18 more than bottom-up) and a total of 184 proteoforms characterized from 14 protein toxin families. Operating both GELFrEE and mass spectrometry to preserve non-covalent interactions generated detailed information about two of the largest venom glycoprotein complexes: the homodimeric l-amino acid oxidase (∼130 kDa) and the multichain toxin cobra venom factor (∼147 kDa). The l-amino acid oxidase complex exhibited two clusters of multiproteoform complexes corresponding to the presence of 5 or 6 N-glycans moieties, each consistent with a distribution of N-acetyl hexosamines. Employing top-down proteomics in both native and denaturing modes provides unprecedented characterization of venom proteoforms and their complexes. A precise molecular inventory of venom proteins will propel the study of snake toxin variation and the targeted development of new antivenoms or other biotherapeutics.

Characterizing the Structure and Oligomerization of Major Royal Jelly Protein 1 (MRJP1) by Mass Spectrometry and Complementary Biophysical Tools.

Royal jelly (RJ) triggers the development of female honeybee larvae into queens. This effect has been attributed to the presence of major royal jelly protein 1 (MRJP1) in RJ. MRJP1 isolated from royal jelly is tightly associated with apisimin, a 54-residue α-helical peptide that promotes the noncovalent assembly of MRJP1 into multimers. No high-resolution structural data are available for these complexes, and their binding stoichiometry remains uncertain. We examined MRJP1/apisimin using a range of biophysical techniques. We also investigated the behavior of deglycosylated samples, as well as samples with reduced apisimin content. Our mass spectrometry (MS) data demonstrate that the native complexes predominantly exist in a (MRJP14 apisimin4) stoichiometry. Hydrogen/deuterium exchange MS reveals that MRJP1 within these complexes is extensively disordered in the range of residues 20-265. Marginally stable secondary structure (likely antiparallel ß-sheet) exists around residues 266-432. These weakly structured regions interchange with conformers that are extensively unfolded, giving rise to bimodal (EX1) isotope distributions. We propose that the native complexes have a "dimer of dimers" quaternary structure in which MRJP1 chains are bridged by apisimin. Specifically, our data suggest that apisimin acts as a linker that forms hydrophobic contacts involving the MRJP1 segment 316VLFFGLV322. Deglycosylation produces large soluble aggregates, highlighting the role of glycans as aggregation inhibitors. Samples with reduced apisimin content form dimeric complexes with a (MRJP12 apisimin1) stoichiometry. The information uncovered in this work will help pave the way toward a better understanding of the unique physiological role played by MRJP1 during queen differentiation.

Native Electron Capture Dissociation Maps to Iron-Binding Channels in Horse Spleen Ferritin.

Native electron capture dissociation (NECD) is a process during which proteins undergo fragmentation similar to that from radical dissociation methods, but without the addition of exogenous electrons. However, after three initial reports of NECD from the cytochrome c dimer complex, no further evidence of the effect has been published. Here, we report NECD behavior from horse spleen ferritin, a ∼490 kDa protein complex ∼20-fold larger than the previously studied cytochrome c dimer. Application of front-end infrared excitation (FIRE) in conjunction with low- and high-m/z quadrupole isolation and collisionally activated dissociation (CAD) provides new insights into the NECD mechanism. Additionally, activation of the intact complex in either the electrospray droplet or the gas phase produced c-type fragment ions. Similar to the previously reported results on cytochrome c, these fragment ions form near residues known to interact with iron atoms in solution. By mapping the location of backbone cleavages associated with c-type ions onto the crystal structure, we are able to characterize two distinct iron binding channels that facilitate iron ion transport into the core of the complex. The resulting pathways are in good agreement with previously reported results for iron binding sites in mammalian ferritin.

Native GELFrEE: a new separation technique for biomolecular assemblies.

The cadre of protein complexes in cells performs an array of functions necessary for life. Their varied structures are foundational to their ability to perform biological functions, lending great import to the elucidation of complex composition and dynamics. Native separation techniques that are operative on low sample amounts and provide high resolution are necessary to gain valuable data on endogenous complexes. Here, we detail and optimize the use of tube gel separations to produce samples proven compatible with native, multistage mass spectrometry (nMS/MS). We find that a continuous system (i.e., no stacking gel) with a gradient in its extent of cross-linking and use of the clear native buffer system performs well for both fractionation and native mass spectrometry of heart extracts and a fungal secretome. This integrated advance in separations and nMS/MS offers the prospect of untargeted proteomics at the next hierarchical level of protein organization in biology.

Probing Asymmetric Charge Partitioning of Protein Oligomers during Tandem Mass Spectrometry.

Dissociation of gaseous protein complexes produced by native electrospray often induces an asymmetric partitioning of charge between ejected subunits. We present a simple asymmetric charge partitioning factor (ACPF) to quantify the magnitude of asymmetry in this effect. When applied to monomer ejection from the cytochrome c dimer and ß-amylase tetramer, we found that the ~60-70% of precursor charge ending up in the ejected monomers corresponds to ACPFs of 1.38 and 2.51, respectively. Further, we used site-specific fragmentation from electron transfer dissociation (ETD) to identify differences in fragmentation and characterize domains of secondary-structure present in the dimer, ejected monomers, and monomers obtained directly from electrospray ionization (ESI). We found evidence of structural changes between the dimer and ejected monomer, but also that the ejected monomer had a nearly identical set of fragment ions produced by ETD as the ESI monomer with the same charge state. Surprisingly, APCF values for ETD fragment ions generated directly from the dimer revealed that the fragments undergo asymmetric charge partitioning at over twice the magnitude of that observed for ejection of the monomer.

Fragmentation of integral membrane proteins in the gas phase.

Integral membrane proteins (IMPs) are of great biophysical and clinical interest because of the key role they play in many cellular processes. Here, a comprehensive top down study of 152 IMPs and 277 soluble proteins from human H1299 cells including 11 087 fragments obtained from collisionally activated dissociation (CAD), 6452 from higher-energy collisional dissociation (HCD), and 2981 from electron transfer dissociation (ETD) shows their great utility and complementarity for the identification and characterization of IMPs. A central finding is that ETD is ∼2-fold more likely to cleave in soluble regions than threshold fragmentation methods, whereas the reverse is observed in transmembrane domains with an observed ∼4-fold bias toward CAD and HCD. The location of charges just prior to dissociation is consistent with this directed fragmentation: protons remain localized on basic residues during ETD but easily mobilize along the backbone during collisional activation. The fragmentation driven by these protons, which is most often observed in transmembrane domains, both is of higher yield and occurs over a greater number of backbone cleavage sites. Further, while threshold dissociation events in transmembrane domains are on average 10.1 (CAD) and 9.2 (HCD) residues distant from the nearest charge site (R, K, H, N-terminus), fragmentation is strongly influenced by the N- or C-terminal position relative to that site: the ratio of observed b- to y-fragments is ∼1:3 if the cleavage occurs >7 residues N-terminal and ∼3:1 if it occurs >7 residues C-terminal to the nearest basic site. Threshold dissociation products driven by a mobilized proton appear to be strongly dependent on not only relative position of a charge site but also N- or C-terminal directionality of proton movement.

Applying label-free quantitation to top down proteomics.

With the prospect of resolving whole protein molecules into their myriad proteoforms on a proteomic scale, the question of their quantitative analysis in discovery mode comes to the fore. Here, we demonstrate a robust pipeline for the identification and stringent scoring of abundance changes of whole protein forms <30 kDa in a complex system. The input is ~100-400 µg of total protein for each biological replicate, and the outputs are graphical displays depicting statistical confidence metrics for each proteoform (i.e., a volcano plot and representations of the technical and biological variation). A key part of the pipeline is the hierarchical linear model that is tailored to the original design of the study. Here, we apply this new pipeline to measure the proteoform-level effects of deleting a histone deacetylase (rpd3) in S. cerevisiae. Over 100 proteoform changes were detected above a 5% false positive threshold in WT vs the Δrpd3 mutant, including the validating observation of hyperacetylation of histone H4 and both H2B isoforms. Ultimately, this approach to label-free top down proteomics in discovery mode is a critical technical advance for testing the hypothesis that whole proteoforms can link more tightly to complex phenotypes in cell and disease biology than do peptides created in shotgun proteomics.

Top Down proteomics: facts and perspectives.

The rise of the "Top Down" method in the field of mass spectrometry-based proteomics has ushered in a new age of promise and challenge for the characterization and identification of proteins. Injecting intact proteins into the mass spectrometer allows for better characterization of post-translational modifications and avoids several of the serious "inference" problems associated with peptide-based proteomics. However, successful implementation of a Top Down approach to endogenous or other biologically relevant samples often requires the use of one or more forms of separation prior to mass spectrometric analysis, which have only begun to mature for whole protein MS. Recent advances in instrumentation have been used in conjunction with new ion fragmentation using photons and electrons that allow for better (and often complete) protein characterization on cases simply not tractable even just a few years ago. Finally, the use of native electrospray mass spectrometry has shown great promise for the identification and characterization of whole protein complexes in the 100 kDa to 1 MDa regime, with prospects for complete compositional analysis for endogenous protein assemblies a viable goal over the coming few years.