Structural basis of the pleiotropic and specific phenotypic consequences of missense mutations in the multifunctional NAD(P)H:quinone oxidoreductase 1 and their pharmacological rescue

Precision proteogenomics reveals pan-cancer impact of germline variants

We investigate the impact of germline variants on cancer patients' proteomes, encompassing 1,064 individuals across 10 cancer types. We introduced an approach, "precision peptidomics," mapping 337,469 coding germline variants onto peptides from patients' mass spectrometry data, revealing their potential impact on post-translational modifications, protein stability, allele-specific expression, and protein structure by leveraging the relevant protein databases. We identified rare pathogenic and common germline variants in cancer genes potentially affecting proteomic features, including variants altering protein abundance and structure and variants in kinases (ERBB2 and MAP2K2) impacting phosphorylation. Precision peptidome analysis predicted destabilizing events in signal-regulatory protein alpha (SIRPA) and glial fibrillary acid protein (GFAP), relevant to immunomodulation and glioblastoma diagnostics, respectively. Genome-wide association studies identified quantitative trait loci for gene expression and protein levels, spanning millions of SNPs and thousands of proteins. Polygenic risk scores correlated with distal effects from risk variants. Our findings emphasize the contribution of germline genetics to cancer heterogeneity and high-throughput precision peptidomics.

The lack of trade-off between conformational stability and binding affinity in a nanobody with therapeutic potential for a misfolding disease

To improve protein pharmaceuticals, we need to balance protein stability and binding affinity with in vivo efficiency. We have recently developed a nanobody (NB-AGT-2) against the alanine:glyoxylate aminotransferase with high stability (Tm ∼ 86 °C) that may be useful to treat a misfolding disease called primary hyperoxaluria type 1. In this work, we characterize the relationships between protein stability and binding affinity in NB-AGT-2 by generating single and double cavity-creating mutants in its hydrophobic core. These mutations decrease thermal stability by 10-20 °C, reflecting changes in thermodynamic stability of up to 8 kcal·mol-1, hardly affecting their binding affinity for its target. Statistical mechanical analysis support long-range propagation of stability effects due to mutations. Our results thus show that NB stability can be largely challenged without an effect on its binding.

Phosphorylation of cytosolic hPGK1 affects protein stability and ligand binding: implications for its subcellular targeting in cancer

Human phosphoglycerate kinase 1(hPGK1) is a key glycolytic enzyme that regulates the balance between ADP and ATP concentrations inside the cell. Phosphorylation of hPGK1 at S203 and S256 has been associated with enzyme import from the cytosol to the mitochondria and the nucleus respectively. These changes in subcellular locations drive tumorigenesis and are likely associated with site-specific changes in protein stability. In this work, we investigate the effects of site-specific phosphorylation on thermal and kinetic stability and protein structural dynamics by hydrogen-deuterium exchange (HDX) and molecular dynamics (MD) simulations. We also investigate the binding of 3-phosphoglycerate and Mg-ADP using these approaches. We show that the phosphomimetic mutation S256D reduces hPGK1 kinetic stability by 50-fold, with no effect of the mutation S203D. Calorimetric studies of ligand binding show a large decrease in affinity for Mg-ADP in the S256D variant, whereas Mg-ADP binding to the WT and S203D can be accurately investigated using protein kinetic stability and binding thermodynamic models. HDX and MD simulations confirmed the destabilization caused by the mutation S256D (with some long-range effects on stability) and its reduced affinity for Mg-ADP due to the strong destabilization of its binding site (particularly in the apo-state). Our research provides evidence suggesting that modifications in protein stability could potentially enhance the translocation of hPGK1 to the nucleus in cancer. While the structural and energetic basis of its mitochondrial import remain unknown.

Structural dynamics and functional cooperativity of human NQO1 by ambient temperature serial crystallography and simulations

The human NQO1 (hNQO1) is a flavin adenine nucleotide (FAD)-dependent oxidoreductase that catalyzes the two-electron reduction of quinones to hydroquinones, being essential for the antioxidant defense system, stabilization of tumor suppressors, and activation of quinone-based chemotherapeutics. Moreover, it is overexpressed in several tumors, which makes it an attractive cancer drug target. To decipher new structural insights into the flavin reductive half-reaction of the catalytic mechanism of hNQO1, we have carried serial crystallography experiments at new ID29 beamline of the ESRF to determine, to the best of our knowledge, the first structure of the hNQO1 in complex with NADH. We have also performed molecular dynamics simulations of free hNQO1 and in complex with NADH. This is the first structural evidence that the hNQO1 functional cooperativity is driven by structural communication between the active sites through long-range propagation of cooperative effects across the hNQO1 structure. Both structural results and MD simulations have supported that the binding of NADH significantly decreases protein dynamics and stabilizes hNQO1 especially at the dimer core and interface. Altogether, these results pave the way for future time-resolved studies, both at x-ray free-electron lasers and synchrotrons, of the dynamics of hNQO1 upon binding to NADH as well as during the FAD cofactor reductive half-reaction. This knowledge will allow us to reveal unprecedented structural information of the relevance of the dynamics during the catalytic function of hNQO1.

Insights into the pathogenesis of primary hyperoxaluria type I from the structural dynamics of alanine:glyoxylate aminotransferase variants

Primary hyperoxaluria type I (PH1) is caused by deficient alanine:glyoxylate aminotransferase (AGT) activity. PH1-causing mutations in AGT lead to protein mistargeting and aggregation. Here, we use hydrogen-deuterium exchange (HDX) to characterize the wild-type (WT), the LM (a polymorphism frequent in PH1 patients) and the LM G170R (the most common mutation in PH1) variants of AGT. We provide the first experimental analysis of AGT structural dynamics, showing that stability is heterogeneous in the native state and providing a blueprint for frustrated regions with potentially functional relevance. The LM and LM G170R variants only show local destabilization. Enzymatic transamination of the pyridoxal 5-phosphate cofactor bound to AGT hardly affects stability. Our study, thus, supports that AGT misfolding is not caused by dramatic effects on structural dynamics.

Gentisic acid ameliorates type 2 diabetes induced by Nicotinamide-Streptozotocin in male mice by attenuating pancreatic oxidative stress and inflammation through modulation of Nrf2 and NF-кB pathways

AIMS

There is a close link between oxidative stress, inflammation, and type 2 diabetes mellitus (T2DM). Gentisic acid (GA) is a di-phenolic compound and an active metabolite of aspirin that possesses antioxidant and anti-inflammatory properties, but its potential anti-diabetic effects have not been evaluated so far. Therefore, this study aimed to evaluate GA's potential antidiabetic effects through the Nuclear Factor Erythroid 2-Related Factor (Nrf2) and Nuclear Factor Kappa Beta (NF-кB) signaling pathways.

MATERIAL AND METHODS

In this study, T2DM induced by a single intraperitoneal injection of STZ (65 mg/kg B.W) after 15 min nicotinamide (120 mg/kg B.W) injection. After seven days of injections, fasting blood glucose (FBS) was measured. Seven days after FBS monitoring treatments started. Grouping and treatments were as follows: 1) Normal control group; NC, 2) Diabetic control group; DC, 3) Metformin group; MT (150 mg/kg B.W, daily), 4) Test group; GA (100 mg/kg B.W, daily). Treatments continued for 14 consecutive days.

KEY FINDINGS

Diabetic mice treatment with GA significantly decreased FBS, improved plasma lipid profiles and pancreatic antioxidant status. GA modulated Nrf2 pathway by upregulation of Nrf2 protein, NAD(P)H: quinone oxidoreductase 1 (Nqo1), and p21, and downregulation of miR-200a, Kelch-like ECH-associated protein 1 (Keap1), and nicotinamide adenine dinucleotide phosphate oxidase-2 (NOX2). Also, GA attenuated inflammation by upregulation of metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) and interleukin-10 (IL-10) and downregulation of miR-125b, NF-кB, tumor necrosis factor-alpha (TNF-α), and interleukin-1 beta (IL-1ß).

SIGNIFICANCE

GA attenuates T2DM, possibly by improving antioxidant status through the Nrf2 pathway and attenuation of inflammation.

Modular droplet injector for sample conservation providing new structural insight for the conformational heterogeneity in the disease-associated NQO1 enzyme

Droplet injection strategies are a promising tool to reduce the large amount of sample consumed in serial femtosecond crystallography (SFX) measurements at X-ray free electron lasers (XFELs) with continuous injection approaches. Here, we demonstrate a new modular microfluidic droplet injector (MDI) design that was successfully applied to deliver microcrystals of the human NAD(P)H:quinone oxidoreductase 1 (NQO1) and phycocyanin. We investigated droplet generation conditions through electrical stimulation for both protein samples and implemented hardware and software components for optimized crystal injection at the Macromolecular Femtosecond Crystallography (MFX) instrument at the Stanford Linac Coherent Light Source (LCLS). Under optimized droplet injection conditions, we demonstrate that up to 4-fold sample consumption savings can be achieved with the droplet injector. In addition, we collected a full data set with droplet injection for NQO1 protein crystals with a resolution up to 2.7 Å, leading to the first room-temperature structure of NQO1 at an XFEL. NQO1 is a flavoenzyme associated with cancer, Alzheimer's and Parkinson's disease, making it an attractive target for drug discovery. Our results reveal for the first time that residues Tyr128 and Phe232, which play key roles in the function of the protein, show an unexpected conformational heterogeneity at room temperature within the crystals. These results suggest that different substates exist in the conformational ensemble of NQO1 with functional and mechanistic implications for the enzyme's negative cooperativity through a conformational selection mechanism. Our study thus demonstrates that microfluidic droplet injection constitutes a robust sample-conserving injection method for SFX studies on protein crystals that are difficult to obtain in amounts necessary for continuous injection, including the large sample quantities required for time-resolved mix-and-inject studies.

Riboflavin kinase and pyridoxine 5′-phosphate oxidase complex formation envisages transient interactions for FMN cofactor delivery

Enzymes catalysing sequential reactions have developed different mechanisms to control the transport and flux of reactants and intermediates along metabolic pathways, which usually involve direct transfer of metabolites from an enzyme to the next one in a cascade reaction. Despite the fact that metabolite or substrate channelling has been widely studied for reactant molecules, such information is seldom available for cofactors in general, and for flavins in particular. Flavin adenine dinucleotide (FAD) and flavin mononucleotide (FMN) act as cofactors in flavoproteins and flavoenzymes involved in a wide range of physiologically relevant processes in all type of organisms. Homo sapiens riboflavin kinase (RFK) catalyses the biosynthesis of the flavin mononucleotide cofactor, and might directly interplay with its flavin client apo-proteins prior to the cofactor transfer. Non-etheless, none of such complexes has been characterized at molecular or atomic level so far. Here, we particularly evaluate the interaction of riboflavin kinase with one of its potential FMN clients, pyridoxine-5′-phosphate oxidase (PNPOx). The interaction capacity of both proteins is assessed by using isothermal titration calorimetry, a methodology that allows to determine dissociation constants for interaction in the micromolar range (in agreement with the expected transient nature of the interaction). Moreover, we show that; i) both proteins become thermally stabilized upon mutual interaction, ii) the tightly bound FMN product can be transferred from RFK to the apo-form of PNPOx producing an efficient enzyme, and iii) the presence of the apo-form of PNPOx slightly enhances RFK catalytic efficiency. Finally, we also show a computational study to predict likely RFK-PNPOx binding modes that can envisage coupling between the FMN binding cavities of both proteins for the potential transfer of FMN.

Phenotypic Modulation of Cancer-Associated Antioxidant NQO1 Activity by Post-Translational Modifications and the Natural Diversity of the Human Genome

Human NAD(P)H:quinone oxidoreductase 1 (hNQO1) is a multifunctional and antioxidant stress protein whose expression is controlled by the Nrf2 signaling pathway. hNQO1 dysregulation is associated with cancer and neurological disorders. Recent works have shown that its activity is also modulated by different post-translational modifications (PTMs), such as phosphorylation, acetylation and ubiquitination, and these may synergize with naturally-occurring and inactivating polymorphisms and mutations. Herein, I describe recent advances in the study of the effect of PTMs and genetic variations on the structure and function of hNQO1 and their relationship with disease development in different genetic backgrounds, as well as the physiological roles of these modifications. I pay particular attention to the long-range allosteric effects exerted by PTMs and natural variation on the multiple functions of hNQO1.

Molecular Mechanisms, Genotype-Phenotype Correlations and Patient-Specific Treatments in Inherited Metabolic Diseases

Advances in DNA sequencing technologies are revealing a vast genetic heterogeneity in human population, which may predispose to metabolic alterations if the activity of metabolic enzymes is affected [...].

Phosphorylation of Thr9 Affects the Folding Landscape of the N-Terminal Segment of Human AGT Enhancing Protein Aggregation of Disease-Causing Mutants

The mutations G170R and I244T are the most common disease cause in primary hyperoxaluria type I (PH1). These mutations cause the misfolding of the AGT protein in the minor allele AGT-LM that contains the P11L polymorphism, which may affect the folding of the N-terminal segment (NTT-AGT). The NTT-AGT is phosphorylated at T9, although the role of this event in PH1 is unknown. In this work, phosphorylation of T9 was mimicked by introducing the T9E mutation in the NTT-AGT peptide and the full-length protein. The NTT-AGT conformational landscape was studied by circular dichroism, NMR, and statistical mechanical methods. Functional and stability effects on the full-length AGT protein were characterized by spectroscopic methods. The T9E and P11L mutations together reshaped the conformational landscape of the isolated NTT-AGT peptide by stabilizing ordered conformations. In the context of the full-length AGT protein, the T9E mutation had no effect on the overall AGT function or conformation, but enhanced aggregation of the minor allele (LM) protein and synergized with the mutations G170R and I244T. Our findings indicate that phosphorylation of T9 may affect the conformation of the NTT-AGT and synergize with PH1-causing mutations to promote aggregation in a genotype-specific manner. Phosphorylation should be considered a novel regulatory mechanism in PH1 pathogenesis.

FAD-deficient P187S mutation of NAD(P)H:quinone oxidoreductase 1 (NQO1*2) binds and accelerates β-amyloid aggregation

Alzheimer's disease (AD) is one of the most prominent neurodegenerative diseases. Results from animal and cellular models suggest that FAD-deficient forms of NAD(P)H quinone oxidoreductase 1 (NQO1) may accelerate the aggregation of Alzheimer's amyloid-β peptide (Aβ1-42). Here, we examined in vitro whether NQO1 and its FAD-deficient P187S mutation (NQO1*2) directly interact with Aβ1-42 and modify its rate of aggregation. When monitored using the fluorescence of either noncovalent thioflavin T (ThT) or HiLyte Fluor 647 (HF647) dye covalently attached to the Aβ1-42 peptide, the aggregation kinetics of Aβ1-42 were markedly more rapid in the presence of NQO1*2 than the wild-type (WT) NQO1. Experiments using apo-NQO1 indicate that this increase is linked to the inability of NQO1*2 to bind to FAD. Furthermore, dicoumarol, an NQO1 inhibitor that binds near the FAD-binding site and stabilizes NQO1*2, markedly decreased the aggregation kinetics of Aβ1-42. Imaging flow cytometry confirmed in-vitro coaggregation of NQO1 isoforms and Aβ1-42. Aβ1-42 alone forms rod-shaped fibril structures while in the presence of NQO1 isoforms, Aβ1-42 is incorporated in the middle of larger globular protein aggregates surrounded by NQO1 molecules. Isothermal titration calorimetry (ITC) analysis indicates that Aβ1-42 interacts with NQO1 isoforms with a specific stoichiometry through a hydrophobic interaction with positive enthalpy and entropy changes. These data define the kinetics, mechanism, and shape of coaggregates of Aβ1-42 and NQO1 isoforms and the potential relevance of FAD-deficient forms of NQO1 for amyloid aggregation diseases.

Pey A. Effect of naturally-occurring mutations on the stability and function of cancer-associated NQO1: Comparison of experiments and computation

Recent advances in DNA sequencing technologies are revealing a large individual variability of the human genome. Our capacity to establish genotype-phenotype correlations in such large-scale is, however, limited. This task is particularly challenging due to the multifunctional nature of many proteins. Here we describe an extensive analysis of the stability and function of naturally-occurring variants (found in the COSMIC and gnomAD databases) of the cancer-associated human NAD(P)H:quinone oxidoreductase 1 (NQO1). First, we performed in silico saturation mutagenesis studies (>5,000 substitutions) aimed to identify regions in NQO1 important for stability and function. We then experimentally characterized twenty-two naturally-occurring variants in terms of protein levels during bacterial expression, solubility, thermal stability, and coenzyme binding. These studies showed a good overall correlation between experimental analysis and computational predictions; also the magnitude of the effects of the substitutions are similarly distributed in variants from the COSMIC and gnomAD databases. Outliers in these experimental-computational genotype-phenotype correlations remain, and we discuss these on the grounds and limitations of our approaches. Our work represents a further step to characterize the mutational landscape of NQO1 in the human genome and may help to improve high-throughput in silico tools for genotype-phenotype correlations in this multifunctional protein associated with disease.

Different phenotypic outcome due to site-specific phosphorylation in the cancer-associated NQO1 enzyme studied by phosphomimetic mutations

Protein phosphorylation is a common phenomenon in human flavoproteins although the functional consequences of this site-specific modification are largely unknown. Here, we evaluated the effects of site-specific phosphorylation (using phosphomimetic mutations at sites S40, S82 and T128) on multiple functional aspects as well as in the structural stability of the antioxidant and disease-associated human flavoprotein NQO1 using biophysical and biochemical methods. In vitro biophysical studies revealed effects of phosphorylation at different sites such as decreased binding affinity for FAD and structural stability of its binding site (S82), conformational stability (S40 and S82) and reduced catalytic efficiency and functional cooperativity (T128). Local stability measurements by H/D exchange in different ligation states provided structural insight into these effects. Transfection of eukaryotic cells showed that phosphorylation at sites S40 and S82 may reduce steady-levels of NQO1 protein by enhanced proteasome-induced degradation. We show that site-specific phosphorylation of human NQO1 may cause pleiotropic and counterintuitive effects on this multifunctional protein with potential implications for its relationships with human disease. Our approach allows to establish relationships between site-specific phosphorylation, functional and structural stability effects in vitro and inside cells paving the way for more detailed analyses of phosphorylation at the flavoproteome scale.

Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype-phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Allosteric Communication in the Multifunctional and Redox NQO1 Protein Studied by Cavity-Making Mutations

Allosterism is a common phenomenon in protein biochemistry that allows rapid regulation of protein stability; dynamics and function. However, the mechanisms by which allosterism occurs (by mutations or post-translational modifications (PTMs)) may be complex, particularly due to long-range propagation of the perturbation across protein structures. In this work, we have investigated allosteric communication in the multifunctional, cancer-related and antioxidant protein NQO1 by mutating several fully buried leucine residues (L7, L10 and L30) to smaller residues (V, A and G) at sites in the N-terminal domain. In almost all cases, mutated residues were not close to the FAD or the active site. Mutations L→G strongly compromised conformational stability and solubility, and L30A and L30V also notably decreased solubility. The mutation L10A, closer to the FAD binding site, severely decreased FAD binding affinity (≈20 fold vs. WT) through long-range and context-dependent effects. Using a combination of experimental and computational analyses, we show that most of the effects are found in the apo state of the protein, in contrast to other common polymorphisms and PTMs previously characterized in NQO1. The integrated study presented here is a first step towards a detailed structural–functional mapping of the mutational landscape of NQO1, a multifunctional and redox signaling protein of high biomedical relevance.

Targeting HIF-1α Function in Cancer through the Chaperone Action of NQO1: Implications of Genetic Diversity of NQO1

HIF-1α is a master regulator of oxygen homeostasis involved in different stages of cancer development. Thus, HIF-1α inhibition represents an interesting target for anti-cancer therapy. It was recently shown that the HIF-1α interaction with NQO1 inhibits proteasomal degradation of the former, thus suggesting that targeting the stability and/or function of NQO1 could lead to the destabilization of HIF-1α as a therapeutic approach. Since the molecular interactions of NQO1 with HIF-1α are beginning to be unraveled, in this review we discuss: (1) Structure–function relationships of HIF-1α; (2) our current knowledge on the intracellular functions and stability of NQO1; (3) the pharmacological modulation of NQO1 by small ligands regarding function and stability; (4) the potential effects of genetic variability of NQO1 in HIF-1α levels and function; (5) the molecular determinants of NQO1 as a chaperone of many different proteins including cancer-associated factors such as HIF-1α, p53 and p73α. This knowledge is then further discussed in the context of potentially targeting the intracellular stability of HIF-1α by acting on its chaperone, NQO1. This could result in novel anti-cancer therapies, always considering that the substantial genetic variability in NQO1 would likely result in different phenotypic responses among individuals.

A single evolutionarily divergent mutation determines the different FAD-binding affinities of human and rat NQO1 due to site-specific phosphorylation

The phosphomimetic mutation S82D in the cancer-associated, FAD-dependent human NADP(H):quinone oxidoreductase 1 (hNQO1) causes a decrease in flavin-adenine dinucleotide-binding affinity and intracellular stability. We test in this work whether the evolutionarily recent neutral mutation R80H in the vicinity of S82 may alter the strong functional effects of S82 phosphorylation through electrostatic interactions. We show using biophysical and bioinformatic analyses that the reverse mutation H80R prevents the effects of S82D phosphorylation on hNQO1 by modulating the local stability. Consistently, in rat NQO1 (rNQO1) which contains R80, the effects of phosphorylation were milder, resembling the behaviour found in hNQO1 when this residue was humanized in rNQO1 (by the R80H mutation). Thus, apparently neutral and evolutionarily divergent mutations may determine the functional response of mammalian orthologues towards phosphorylation.