Sequence-based bioprospecting of myo-inositol oxygenase (Miox) reveals new homologues that increase glucaric acid production in Saccharomyces cerevisiae.

myo-Inositol oxygenase (Miox) is a rate-limiting enzyme for glucaric acid production via microbial fermentation. The enzyme converts myo-inositol to glucuronate, which is further converted to glucaric acid, a natural compound with industrial uses that range from detergents to pharmaceutical synthesis to polymeric materials. More than 2,000 Miox sequences are available in the Uniprot database but only thirteen are classified as reviewed in Swiss-Prot (August 2019). In this study, sequence similarity networks were used to identify new homologues to be expressed in Saccharomyces cerevisiae for glucaric acid production. The expression of four homologues did not lead to product formation. Some of these enzymes may have a defective "dynamic lid" - a structural feature important to close the reaction site - which might explain the lack of activity. Thirty-one selected Miox sequences did allow for product formation, of which twenty-five were characterized for the first time. Expression of Talaromyces marneffei Miox led to the accumulation of 1.76 ±â€¯0.33 g glucaric acid/L from 20 g glucose/L and 10 g/L myo-inositol. Specific glucaric acid titer with TmMiox increased 44 % compared to the often-used Arabidopsis thaliana variant AtMiox4 (0.258 vs. 0.179 g glucaric acid/g biomass). AtMiox4 activity decreased from 12.47 to 0.40 nmol/min/mg protein when cells exited exponential phase during growth on glucose, highlighting the importance of future research on Miox stability in order to further improve microbial production of glucaric acid.

Target discovery using biobanks and human genetics.

Large-scale biobanks can yield unprecedented insights into our health and provide discoveries of new and potentially targetable biomarkers. Several protective loss-of-function alleles have been identified, including variants that protect against cardiovascular disease, obesity, type 2 diabetes, and asthma and allergic diseases. These alleles serve as indicators of efficacy, mimicking the effects of drugs and suggesting that inhibiting these genes could provide therapeutic benefit, as has been observed for PCSK9. We provide a context for these findings through a multifaceted review covering the use of genetics in drug discovery efforts through genome-wide and phenome-wide association studies, linking deep mutation scanning data to molecular function and highlighting some additional tools that might help in the interpretation of newly discovered variants.

Rational design of thiolase substrate specificity for metabolic engineering applications.

Metabolic engineering efforts require enzymes that are both highly active and specific toward the synthesis of a desired output product to be commercially feasible. The 3-hydroxyacid (3HA) pathway, also known as the reverse ß-oxidation or coenzyme-A-dependent chain-elongation pathway, can allow for the synthesis of dozens of useful compounds of various chain lengths and functionalities. However, this pathway suffers from byproduct formation, which lowers the yields of the desired longer chain products, as well as increases downstream separation costs. The thiolase enzyme catalyzes the first reaction in this pathway, and its substrate specificity at each of its two catalytic steps sets the chain length and composition of the chemical scaffold upon which the other downstream enzymes act. However, there have been few attempts reported in the literature to rationally engineer thiolase substrate specificity. In this study, we present a model-guided, rational design study of ordered substrate binding applied to two biosynthetic thiolases, with the goal of increasing the ratio of C6/C4 products formed by the 3HA pathway, 3-hydroxy-hexanoic acid and 3-hydroxybutyric acid. We identify thiolase mutants that result in nearly 10-fold increases in C6/C4 selectivity. Our findings can extend to other pathways that employ the thiolase for chain elongation, as well as expand our knowledge of sequence-structure-function relationship for this important class of enzymes.

The human noncoding genome defined by genetic diversity.

Understanding the significance of genetic variants in the noncoding genome is emerging as the next challenge in human genomics. We used the power of 11,257 whole-genome sequences and 16,384 heptamers (7-nt motifs) to build a map of sequence constraint for the human species. This build differed substantially from traditional maps of interspecies conservation and identified regulatory elements among the most constrained regions of the genome. Using new Hi-C experimental data, we describe a strong pattern of coordination over 2 Mb where the most constrained regulatory elements associate with the most essential genes. Constrained regions of the noncoding genome are up to 52-fold enriched for known pathogenic variants as compared to unconstrained regions (21-fold when compared to the genome average). This map of sequence constraint across thousands of individuals is an asset to help interpret noncoding elements in the human genome, prioritize variants and reconsider gene units at a larger scale.

Biocuration in the structure-function linkage database: the anatomy of a superfamily.

With ever-increasing amounts of sequence data available in both the primary literature and sequence repositories, there is a bottleneck in annotating molecular function to a sequence. This article describes the biocuration process and methods used in the structure-function linkage database (SFLD) to help address some of the challenges. We discuss how the hierarchy within the SFLD allows us to infer detailed functional properties for functionally diverse enzyme superfamilies in which all members are homologous, conserve an aspect of their chemical function and have associated conserved structural features that enable the chemistry. Also presented is the Enzyme Structure-Function Ontology (ESFO), which has been designed to capture the relationships between enzyme sequence, structure and function that underlie the SFLD and is used to guide the biocuration processes within the SFLD. Database URL: http://sfld.rbvi.ucsf.edu/.

Porting the synthetic D-glucaric acid pathway from Escherichia coli to Saccharomyces cerevisiae.

D-Glucaric acid can be produced as a value-added chemical from biomass through a de novo pathway in Escherichia coli. However, previous studies have identified pH-mediated toxicity at product concentrations of 5 g/L and have also found the eukaryotic myo-inositol oxygenase (MIOX) enzyme to be rate-limiting. We ported this pathway to Saccaromyces cerevisiae, which is naturally acid-tolerant and evaluate a codon-optimized MIOX homologue. We constructed two engineered yeast strains that were distinguished solely by their MIOX gene - either the previous version from Mus musculus or a homologue from Arabidopsis thaliana codon-optimized for expression in S. cerevisiae - in order to identify the rate-limiting steps for D-glucaric acid production both from a fermentative and non-fermentative carbon source. myo-Inositol availability was found to be rate-limiting from glucose in both strains and demonstrated to be dependent on growth rate, whereas the previously used M. musculus MIOX activity was found to be rate-limiting from glycerol. Maximum titers were 0.56 g/L from glucose in batch mode, 0.98 g/L from glucose in fed-batch mode, and 1.6 g/L from glucose supplemented with myo-inositol. Future work focusing on the MIOX enzyme, the interplay between growth and production modes, and promoting aerobic respiration should further improve this pathway.

Mechanistic and bioinformatic investigation of a conserved active site helix in α-isopropylmalate synthase from Mycobacterium tuberculosis, a member of the DRE-TIM metallolyase superfamily.

The characterization of functionally diverse enzyme superfamilies provides the opportunity to identify evolutionarily conserved catalytic strategies, as well as amino acid substitutions responsible for the evolution of new functions or specificities. Isopropylmalate synthase (IPMS) belongs to the DRE-TIM metallolyase superfamily. Members of this superfamily share common active site elements, including a conserved active site helix and an HXH divalent metal binding motif, associated with stabilization of a common enolate anion intermediate. These common elements are overlaid by variations in active site architecture resulting in the evolution of a diverse set of reactions that include condensation, lyase/aldolase, and carboxyl transfer activities. Here, using IPMS, an integrated biochemical and bioinformatics approach has been utilized to investigate the catalytic role of residues on an active site helix that is conserved across the superfamily. The construction of a sequence similarity network for the DRE-TIM metallolyase superfamily allows for the biochemical results obtained with IPMS variants to be compared across superfamily members and within other condensation-catalyzing enzymes related to IPMS. A comparison of our results with previous biochemical data indicates an active site arginine residue (R80 in IPMS) is strictly required for activity across the superfamily, suggesting that it plays a key role in catalysis, most likely through enolate stabilization. In contrast, differential results obtained from substitution of the C-terminal residue of the helix (Q84 in IPMS) suggest that this residue plays a role in reaction specificity within the superfamily.

Bioprospecting in the genomic age.

The genomic revolution promises great advances in the search for useful biocatalysts. Function-based metagenomic approaches have identified several enzymes with properties that make them useful candidates for a variety of bioprocesses. As DNA sequencing costs continue to decline, the volume of genomic data, along with their corresponding predicted protein sequences, will continue to increase dramatically, necessitating new approaches to leverage this information for gene-based bioprospecting efforts. Additionally, as new functions are discovered and correlated with this sequence information, the knowledge of the often complex relationship between a protein's sequence and function will improve. This in turn will lead to better gene-based bioprospecting approaches and facilitate the tailoring of desired properties through protein engineering projects. In this chapter, we discuss a number of recent advances in bioprospecting within the context of the genomic age.

The Structure-Function Linkage Database.

The Structure-Function Linkage Database (SFLD, http://sfld.rbvi.ucsf.edu/) is a manually curated classification resource describing structure-function relationships for functionally diverse enzyme superfamilies. Members of such superfamilies are diverse in their overall reactions yet share a common ancestor and some conserved active site features associated with conserved functional attributes such as a partial reaction. Thus, despite their different functions, members of these superfamilies 'look alike', making them easy to misannotate. To address this complexity and enable rational transfer of functional features to unknowns only for those members for which we have sufficient functional information, we subdivide superfamily members into subgroups using sequence information, and lastly into families, sets of enzymes known to catalyze the same reaction using the same mechanistic strategy. Browsing and searching options in the SFLD provide access to all of these levels. The SFLD offers manually curated as well as automatically classified superfamily sets, both accompanied by search and download options for all hierarchical levels. Additional information includes multiple sequence alignments, tab-separated files of functional and other attributes, and sequence similarity networks. The latter provide a new and intuitively powerful way to visualize functional trends mapped to the context of sequence similarity.

Residues required for activity in Escherichia coli o-succinylbenzoate synthase (OSBS) are not conserved in all OSBS enzymes.

Understanding how enzyme specificity evolves will provide guiding principles for protein engineering and function prediction. The o-succinylbenzoate synthase (OSBS) family is an excellent model system for elucidating these principles because it has many highly divergent amino acid sequences that are <20% identical, and some members have evolved a second function. The OSBS family belongs to the enolase superfamily, members of which use a set of conserved residues to catalyze a wide variety of reactions. These residues are the only conserved residues in the OSBS family, so they are not sufficient to determine reaction specificity. Some enzymes in the OSBS family catalyze another reaction, N-succinylamino acid racemization (NSAR). NSARs cannot be segregated into a separate family because their sequences are highly similar to those of known OSBSs, and many of them have both OSBS and NSAR activities. To determine how such divergent enzymes can catalyze the same reaction and how NSAR activity evolved, we divided the OSBS family into subfamilies and compared the divergence of their active site residues. Correlating sequence conservation with the effects of mutations in Escherichia coli OSBS identified two nonconserved residues (R159 and G288) at which mutations decrease efficiency ≥200-fold. These residues are not conserved in the subfamily that includes NSAR enzymes. The OSBS/NSAR subfamily binds the substrate in a different orientation, eliminating selective pressure to retain arginine and glycine at these positions. This supports the hypothesis that specificity-determining residues have diverged in the OSBS family and provides insight into the sequence changes required for the evolution of NSAR activity.

The evolution of function in strictosidine synthase-like proteins.

The exponential growth of sequence data provides abundant information for the discovery of new enzyme reactions. Correctly annotating the functions of highly diverse proteins can be difficult, however, hindering use of this information. Global analysis of large superfamilies of related proteins is a powerful strategy for understanding the evolution of reactions by identifying catalytic commonalities and differences in reaction and substrate specificity, even when only a few members have been biochemically or structurally characterized. A comparison of >2500 sequences sharing the six-bladed ß-propeller fold establishes sequence, structural, and functional links among the three subgroups of the functionally diverse N6P superfamily: the arylesterase-like and senescence marker protein-30/gluconolactonase/luciferin-regenerating enzyme-like (SGL) subgroups, representing enzymes that catalyze lactonase and related hydrolytic reactions, and the so-called strictosidine synthase-like (SSL) subgroup. Metal-coordinating residues were identified as broadly conserved in the active sites of all three subgroups except for a few proteins from the SSL subgroup, which have been experimentally determined to catalyze the quite different strictosidine synthase (SS) reaction, a metal-independent condensation reaction. Despite these differences, comparison of conserved catalytic features of the arylesterase-like and SGL enzymes with the SSs identified similar structural and mechanistic attributes between the hydrolytic reactions catalyzed by the former and the condensation reaction catalyzed by SS. The results also suggest that despite their annotations, the great majority of these >500 SSL sequences do not catalyze the SS reaction; rather, they likely catalyze hydrolytic reactions typical of the other two subgroups instead. This prediction was confirmed experimentally for one of these proteins.