Human NQO1 [NAD(H):quinone oxidoreductase 1] is a multi-functional and stress-inducible dimeric protein involved in the antioxidant defense, the activation of cancer prodrugs and the stabilization of oncosuppressors. Despite its roles in human diseases, such as cancer and neurological disorders, a detailed characterization of its enzymatic cycle is still lacking. In this work, we provide a comprehensive analysis of the NQO1 catalytic cycle using rapid mixing techniques, including multiwavelength and spectral deconvolution studies, kinetic modeling and temperature-dependent kinetic isotope effects (KIEs). Our results systematically support the existence of two pathways for hydride transfer throughout the NQO1 catalytic cycle, likely reflecting that the two active sites in the dimer catalyze two-electron reduction with different rates, consistent with the cooperative binding of inhibitors such as dicoumarol. This negative cooperativity in NQO1 redox activity represents a sort of half-of-sites activity. Analysis of KIEs and their temperature dependence also show significantly different contributions from quantum tunneling, structural dynamics and reorganizations to catalysis at the two active sites. Our work will improve our understanding of the effects of cancer-associated single amino acid variants and post-translational modifications in this protein of high relevance in cancer progression and treatment.

The multifunctional nature of human flavoproteins is critically linked to their ability to populate multiple conformational states. Ligand binding, post-translational modifications and disease-associated mutations can reshape this functional landscape, although the structure-function relationships of these effects are not well understood. Herein, we characterized the structural and functional consequences of two mutations (the cancer-associated P187S and the phosphomimetic S82D) on different ligation states which are relevant to flavin binding, intracellular stability and catalysis of the disease-associated NQO1 flavoprotein. We found that these mutations affected the stability locally and their effects propagated differently through the protein structure depending both on the nature of the mutation and the ligand bound, showing directional preference from the mutated site and leading to specific phenotypic manifestations in different functional traits (FAD binding, catalysis and inhibition, intracellular stability and pharmacological response to ligands). Our study thus supports that pleitropic effects of disease-causing mutations and phosphorylation events on human flavoproteins may be caused by long-range structural propagation of stability effects to different functional sites that depend on the ligation-state and site-specific perturbations. Our approach can be of general application to investigate these pleiotropic effects at the flavoproteome scale in the absence of high-resolution structural models.

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.

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.

The biosynthesis of the flavin mononucleotide (FMN) and flavin adenine dinucleotide (FAD), cofactors used by 2% of proteins, occurs through the sequential action of two ubiquitous activities: a riboflavinkinase (RFK) that phosphorylates the riboflavin (RF) precursor to FMN, and a FMN:adenylyltransferase (FMNAT) that transforms FMN into FAD. In most mammals two different monofunctional enzymes have each of these activities, but in prokaryotes a single bifunctional enzyme, FAD synthase (FADS), holds them. Differential structural and functional traits for RFK and FMNAT catalysis between bacteria and mammals, as well as within the few bacterial FADSs so far characterized, has envisaged the potentiality of FADSs from pathogens as targets for the development of species‐specific inhibitors. Here, we particularly characterize the FADS from the ovine pathogen Brucella ovis (BoFADS), causative agent of brucellosis. We show that BoFADS has RFK activity independently of the media redox status, but its FMNAT activity (in both forward and reverse senses) only occurs under strong reducing conditions. Moreover, kinetics for flavin and adenine nucleotides binding to the RFK site show that BoFADS binds preferentially the substrates of the RFK reaction over the products and that the adenine nucleotide must bind prior to flavin entrapment. These results, together with multiple sequence alignments and phylogenetic analysis, point to variability in the less conserved regions as contributing to the species‐specific features in prokaryotic FADSs, including those from pathogens, that allow them to adopt alternative strategies in FMN and FAD biosynthesis and overall flavin homeostasis.

New treatments for diseases caused by antimicrobial-resistant microorganisms can be developed by identifying unexplored therapeutic targets and by designing efficient drug screening protocols. In this study, we have screened a library of compounds to find ligands for the flavin-adenine dinucleotide synthase (FADS) -a potential target for drug design against tuberculosis and pneumonia- by implementing a new and efficient virtual screening protocol. The protocol has been developed for the in silico search of ligands of unexplored therapeutic targets, for which limited information about ligands or ligand-receptor structures is available. It implements an integrative funnel-like strategy with filtering layers that increase in computational accuracy. The protocol starts with a pharmacophore-based virtual screening strategy that uses ligand-free receptor conformations from molecular dynamics (MD) simulations. Then, it performs a molecular docking stage using several docking programs and an exponential consensus ranking strategy. The last filter, samples the conformations of compounds bound to the target using MD simulations. The MD conformations are scored using several traditional scoring functions in combination with a newly-proposed score that takes into account the fluctuations of the molecule with a Morse-based potential. The protocol was optimized and validated using a compound library with known ligands of the Corynebacterium ammoniagenes FADS. Then, it was used to find new FADS ligands from a compound library of 14,000 molecules. A small set of 17 in silico filtered molecules were tested experimentally. We identified five inhibitors of the activity of the flavin adenylyl transferase module of the FADS, and some of them were able to inhibit growth of three bacterial species: C. ammoniagenes, Mycobacterium tuberculosis, and Streptococcus pneumoniae, where the last two are human pathogens. Overall, the results show that the integrative VS protocol is a cost-effective solution for the discovery of ligands of unexplored therapeutic targets. Developing cures for antimicrobial-resistant microorganisms is a pressing necessity. Addressing this problem requires the discovery of novel therapeutic targets -for example, bacterial proteins with no human homologues- and the development of cost-effective drug screening protocols. In this work, we tackled the problem on both sides. We developed an efficient and successful integrative computational protocol for screening inhibitory-molecules for unexplored targets. We used it to discover five novel inhibitors of flavin-adenine dinucleotide synthase (FADS), a promising protein target of pathogens causing tuberculosis and pneumonia.