The Conversion of UDP-Glc to UDP-Man: In Silico and Biochemical Exploration To Improve the Catalytic Efficiency of CDP-Tyvelose C2-Epimerases.

A promiscuous CDP-tyvelose 2-epimerase (TyvE) from Thermodesulfatator atlanticus (TaTyvE) belonging to the nucleotide sugar active short-chain dehydrogenase/reductase superfamily (NS-SDRs) was recently discovered. TaTyvE performs the slow conversion of NDP-glucose (NDP-Glc) to NDP-mannose (NDP-Man). Here, we present the sequence fingerprints that are indicative of the conversion of UDP-Glc to UDP-Man in TyvE-like enzymes based on the heptagonal box motifs. Our data-mining approach led to the identification of 11 additional TyvE-like enzymes for the conversion of UDP-Glc to UDP-Man. We characterized the top two wild-type candidates, which show a 15- and 20-fold improved catalytic efficiency, respectively, on UDP-Glc compared to TaTyvE. In addition, we present a quadruple variant of one of the identified enzymes with a 70-fold improved catalytic efficiency on UDP-Glc compared to TaTyvE. These findings could help the design of new nucleotide production pathways starting from a cheap sugar substrate like glucose or sucrose.

Complete Switch of Reaction Specificity of an Aldolase by Directed Evolution In Vitro: Synthesis of Generic Aliphatic Aldol Products.

A structure-guided engineering of fructose-6-phosphate aldolase was performed to expand its substrate promiscuity toward aliphatic nucleophiles, that is, unsubstituted alkanones and alkanals. A "smart" combinatorial library was created targeting residues D6, T26, and N28, which form a binding pocket around the nucleophilic carbon atom. Double-selectivity screening was executed by high-performance TLC that allowed simultaneous determination of total activity as well as a preference for acetone versus propanal as competing nucleophiles. D6 turned out to be the key residue that enabled activity with non-hydroxylated nucleophiles. Altogether 25 single- and double-site variants (D6X and D6X/T26X) were discovered that show useful synthetic activity and a varying preference for ketone or aldehyde as the aldol nucleophiles. Remarkably, all of the novel variants had completely lost their native activity for cleavage of fructose 6-phosphate.

Bioinformatic analysis of fold-type III PLP-dependent enzymes discovers multimeric racemases.

Pyridoxal-5'-phosphate (PLP)-dependent enzymes are ubiquitous in nature and catalyze a variety of important metabolic reactions. The fold-type III PLP-dependent enzyme family is primarily comprised of decarboxylases and alanine racemases. In the development of a multiple structural alignment database (3DM) for the enzyme family, a large subset of 5666 uncharacterized proteins with high structural, but low sequence similarity to alanine racemase and decarboxylases was found. Compared to these two classes of enzymes, the protein sequences being the object of this study completely lack the C-terminal domain, which has been reported important for the formation of the dimer interface in other fold-type III enzymes. The 5666 sequences cluster around four protein templates, which also share little sequence identity to each other. In this work, these four template proteins were solubly expressed in Escherichia coli, purified, and their substrate profiles were evaluated by HPLC analysis for racemase activity using a broader range of amino acids. They were found active only against alanine or serine, where they exhibited Michaelis constants within the range of typical bacterial alanine racemases, but with significantly lower turnover numbers. As the already described racemases were proposed to be active and appeared to be monomers as judged from their crystal structures, we also investigated this aspect for the four new enzymes. Here, size exclusion chromatography indicated the presence of oligomeric states of the enzymes and a native-PAGE in-gel assay showed that the racemase activity was present only in an oligomeric state but not as monomer. This suggests the likelihood of a different behavior of these enzymes in solution compared to the one observed in crystalline form.

Switch in Cofactor Specificity of a Baeyer-Villiger Monooxygenase.

Baeyer-Villiger monooxygenases (BVMOs) catalyze the oxidation of ketones to esters or lactones by using molecular oxygen and a cofactor. Type I BVMOs display a strong preference for NADPH. However, for industrial purposes NADH is the preferred cofactor, as it is ten times cheaper and more stable. Thus, we created a variant of the cyclohexanone monooxygenase from Acinetobacter sp. NCIMB 9871 (CHMOAcineto ); this used NADH 4200-fold better than NADPH. By combining structure analysis, sequence alignment, and literature data, 21 residues in proximity of the cofactor were identified and targeted for mutagenesis. Two combinatorial variants bearing three or four mutations showed higher conversions of cyclohexanone with NADH (79 %) compared to NADPH (58 %) as well as specificity. The structural reasons for this switch in cofactor specificity of a type I BVMO are especially a hydrogen-bond network coordinating the two hydroxy groups of NADH through direct interactions and bridging water molecules.

Two subtle amino Acid changes in a transaminase substantially enhance or invert enantiopreference in cascade syntheses.

Amine transaminases (ATAs) are powerful enzymes for the stereospecific production of chiral amines. However, the synthesis of amines incorporating more than one stereocenter is still a challenge. We developed a cascade synthesis to access optically active 3-alkyl-substituted chiral amines by combining two asymmetric synthesis steps catalyzed by an enoate reductase and ATAs. The ATA wild type from Vibrio fluvialis showed only modest enantioselectivity (14 % de) in the amination of (S)-3-methylcyclohexanone, the product of the enoate-reductase-catalyzed reaction step. However, by protein engineering we created two variants with substantially improved diastereoselectivities: variant Leu56Val exhibited a higher R selectivity (66 % de) whereas the Leu56Ile substitution caused a switch in enantiopreference to furnish the S-configured diastereomer (70 % de). Addition of 30 % DMSO further improved the selectivity and facilitated the synthesis of (1R,3S)-1-amino-3-methylcyclohexane with 89 % de at 87 % conversion.

Simultaneous use of in silico design and a correlated mutation network as a tool to efficiently guide enzyme engineering.

In order to improve the efficiency of directed evolution experiments, in silico multiple-substrate clustering was combined with an analysis of the variability of natural enzymes within a protein superfamily. This was applied to a Pseudomonas fluorescens esterase (PFE I) targeting the enantioselective hydrolysis of 3-phenylbutyric acid esters. Data reported in the literature for nine substrates were used for the clustering meta-analysis of the docking conformations in wild-type PFE I, and this highlighted a tryptophan residue (W28) as an interesting target. Exploration of the most frequently, naturally occurring amino acids at this position suggested that the reduced flexibility observed in the case of the W28F variant leads to enhancement of the enantioselectivity. This mutant was subsequently combined with mutations identified in a library based on analysis of a correlated mutation network. By interrogation of <80 variants a mutant with 15-fold improved enantioselectivity was found.

Alteration of the Donor/Acceptor Spectrum of the (S)-Amine Transaminase from Vibrio fluvialis.

To alter the amine donor/acceptor spectrum of an (S)-selective amine transaminase (ATA), a library based on the Vibrio fluvialis ATA targeting four residues close to the active site (L56, W57, R415 and L417) was created. A 3DM-derived alignment comprising fold class I pyridoxal-5'-phosphate (PLP)-dependent enzymes allowed identification of positions, which were assumed to determine substrate specificity. These positions were targeted for mutagenesis with a focused alphabet of hydrophobic amino acids to convert an amine:α-keto acid transferase into an amine:aldehyde transferase. Screening of 1200 variants revealed three hits, which showed a shifted amine donor/acceptor spectrum towards aliphatic aldehydes (mainly pentanal), as well as an altered pH profile. Interestingly, all three hits, although found independently, contained the same mutation R415L and additional W57F and L417V substitutions.

NucleaRDB: information system for nuclear receptors.

The NucleaRDB is a Molecular Class-Specific Information System that collects, combines, validates and disseminates large amounts of heterogeneous data on nuclear hormone receptors. It contains both experimental and computationally derived data. The data and knowledge present in the NucleaRDB can be accessed using a number of different interactive and programmatic methods and query systems. A nuclear hormone receptor-specific PDF reader interface is available that can integrate the contents of the NucleaRDB with full-text scientific articles. The NucleaRDB is freely available at http://www.receptors.org/nucleardb.

Consensus engineering of sucrose phosphorylase: the outcome reflects the sequence input.

Consensus engineering, which is replacing amino acids by the most frequently occurring one at their positions in a multiple sequence alignment (MSA), is a known strategy to increase the stability of a protein. The application of this concept to the entire sequence of an enzyme, however, has been tried only a few times mainly because of the problems determining the consensus in highly variable regions. We show that this problem can be solved by replacing such problematic regions by the corresponding sequence of the natural homologue closest to the consensus. When one or a few sub-families are overrepresented in the MSA the consensus sequence is a biased representation of the sequence space. We examine the influence of this bias by constructing three consensus sequences using different MSAs of sucrose phosphorylase (SP). Each consensus enzyme contained about 70 mutations compared to its closest natural homologue and folded correctly and displayed activity on sucrose. Correlation analysis revealed that the family's co-evolution network was kept intact, which is one of the main advantages of full-length consensus design. The consensus enzymes displayed an "average" thermostability, that is, one that is higher than some but not all known representatives. We cautiously present practical rules for the design of consensus sequences, but warn that the measure of success depends on which natural enzyme is used as point of comparison.

Structure-function relationships in NDP-sugar active SDR enzymes: Fingerprints for functional annotation and enzyme engineering.

Short-chain Dehydrogenase/Reductase enzymes that are active on nucleotide sugars (abbreviated as NS-SDR) are of paramount importance in the biosynthesis of rare sugars and glycosides. Some family members have already been extensively characterized due to their direct implication in metabolic disorders or in the biosynthesis of virulence factors. In this review, we combine the knowledge gathered from studies that typically focused only on one NS-SDR activity with an in-depth analysis and overview of all of the different NS-SDR families (169,076 enzyme sequences). Through this structure-based multiple sequence alignment of NS-SDRs retrieved from public databases, we could identify clear patterns in conservation and correlation of crucial residues. Supported by this analysis, we suggest updating and extending the UDP-galactose 4-epimerase "hexagonal box model" to an "heptagonal box model" for all NS-SDR enzymes. This specificity model consists of seven conserved regions surrounding the NDP-sugar substrate that serve as fingerprint for each specificity. The specificity fingerprints highlighted in this review will be beneficial for functional annotation of the large group of NS-SDR enzymes and form a guide for future enzyme engineering efforts focused on the biosynthesis of rare and specialty carbohydrates.

Novel tools for extraction and validation of disease-related mutations applied to Fabry disease.

Genetic disorders are often caused by nonsynonymous nucleotide changes in one or more genes associated with the disease. Specific amino acid changes, however, can lead to large variability of phenotypic expression. For many genetic disorders this results in an increasing amount of publications describing phenotype-associated mutations in disorder-related genes. Keeping up with this stream of publications is essential for molecular diagnostics and translational research purposes but often impossible due to time constraints: there are simply too many articles to read. To help solve this problem, we have created Mutator, an automated method to extract mutations from full-text articles. Extracted mutations are crossreferenced to sequence data and a scoring method is applied to distinguish false-positives. To analyze stored and new mutation data for their (potential) effect we have developed Validator, a Web-based tool specifically designed for DNA diagnostics. Fabry disease, a monogenetic gene disorder of the GLA gene, was used as a test case. A structure-based sequence alignment of the alpha-amylase superfamily was used to validate results. We have compared our data with existing Fabry mutation data sets obtained from the HGMD and Swiss-Prot databases. Compared to these data sets, Mutator extracted 30% additional mutations from the literature.

3DM: systematic analysis of heterogeneous superfamily data to discover protein functionalities.

Ten years of experience with molecular class-specific information systems (MCSIS) such as with the hand-curated G protein-coupled receptor database (GPCRDB) or the semiautomatically generated nuclear receptor database has made clear that a wide variety of questions can be answered when protein-related data from many different origins can be flexibly combined. MCSISes revolve around a multiple sequence alignment (MSA) that includes "all" available sequences from the entire superfamily, and it has been shown at many occasions that the quality of these alignments is the most crucial aspect of the MCSIS approach. We describe here a system called 3DM that can automatically build an entire MCSIS. 3DM bases the MSA on a multiple structure alignment, which implies that the availability of a large number of superfamily members with a known three-dimensional structure is a requirement for 3DM to succeed well. Thirteen MCSISes were constructed and placed on the Internet for examination. These systems have been instrumental in a large series of research projects related to enzyme activity or the understanding and engineering of specificity, protein stability engineering, DNA-diagnostics, drug design, and so forth.

Determinants of the Nucleotide Specificity in the Carbohydrate Epimerase Family 1.

In recent years, carbohydrate epimerases have attracted increasing attention as promising biocatalysts for the production of specialty sugars and derivatives. The vast majority of these enzymes are active on nucleotide-activated sugars, rather than on their free counterparts. Although such epimerases are known to have a clear preference for a particular nucleotide (UDP, GDP, CDP, or ADP), very little is known about the determinants of the respective specificities. In this work, sequence motifs are identified that correlate with the different nucleotide specificities in one of the main epimerase superfamilies, carbohydrate epimerase 1 (CEP1). To confirm their relevance, GDP- and CDP-specific residues are introduced into the UDP-glucose 4-epimerase from Thermus thermophilus, resulting in a 3-fold and 13-fold reduction in KM for GDP-Glc and CDP-Glc, respectively. Moreover, several variants are severely crippled in UDP-Glc activity, which further underlines the crucial role of the identified positions. Hence, the analysis should prove to be valuable for the further exploration and application of epimerases involved in carbohydrate synthesis.

Correlated mutation analyses on super-family alignments reveal functionally important residues.

Correlated mutation analyses (CMA) on multiple sequence alignments are widely used for the prediction of the function of amino acids. The accuracy of CMA-based predictions is mainly determined by the number of sequences, by their evolutionary distances, and by the quality of the alignments. These criteria are best met in structure-based sequence alignments of large super-families. So far, CMA-techniques have mainly been employed to study the receptor interactions. The present work shows how a novel CMA tool, called Comulator, can be used to determine networks of functionally related residues in enzymes. These analyses provide leads for protein engineering studies that are directed towards modification of enzyme specificity or activity. As proof of concept, Comulator has been applied to four enzyme super-families: the isocitrate lyase/phoshoenol-pyruvate mutase super-family, the hexokinase super-family, the RmlC-like cupin super-family, and the FAD-linked oxidases super-family. In each of those cases networks of functionally related residue positions were discovered that upon mutation influenced enzyme specificity and/or activity as predicted. We conclude that CMA is a powerful tool for redesigning enzyme activity and selectivity.

nanoDSF as screening tool for enzyme libraries and biotechnology development.

Enzymes are attractive tools for synthetic applications. To be viable for industrial use, enzymes need sufficient stability towards the desired reaction conditions such as high substrate and cosolvent concentration, non-neutral pH and elevated temperatures. Thermal stability is an attractive feature not only because it allows for protein purification by thermal treatment and higher process temperatures but also due to the associated higher stability against other destabilising factors. Therefore, high-throughput screening (HTS) methods are desirable for the identification of thermostable biocatalysts by discovery from nature or by protein engineering but current methods have low throughput and require time-demanding purification of protein samples. We found that nanoscale differential scanning fluorimetry (nanoDSF) is a valuable tool to rapidly and reliably determine melting points of native proteins. To avoid intrinsic problems posed by crude protein extracts, hypotonic extraction of overexpressed protein from bacterial host cells resulted in higher sample quality and accurate manual determination of several hundred melting temperatures per day. We have probed the use of nanoDSF for HTS of a phylogenetically diverse aldolase library to identify novel thermostable enzymes from metagenomic sources and for the rapid measurements of variants from saturation mutagenesis. The feasibility of nanoDSF for the screening of synthetic reaction conditions was proved by studies of cosolvent tolerance, which showed protein melting temperature to decrease linearly with increasing cosolvent concentration for all combinations of six enzymes and eight water-miscible cosolvents investigated, and of substrate affinity, which showed stabilisation of hexokinase by sugars in the absence of ATP cofactor. ENZYMES: Alcohol dehydrogenase (NADP+ ) (EC 1.1.1.2), transketolase (EC 2.2.1.1), hexokinase (EC 2.7.1.1), 2-deoxyribose-5-phosphate aldolase (EC 4.1.2.4), fructose-6-phosphate aldolase (EC 4.1.2.n).

Distant Non-Obvious Mutations Influence the Activity of a Hyperthermophilic Pyrococcusfuriosus Phosphoglucose Isomerase.

The cupin-type phosphoglucose isomerase (PfPGI) from the hyperthermophilic archaeon Pyrococcus furiosus catalyzes the reversible isomerization of glucose-6-phosphate to fructose-6-phosphate. We investigated PfPGI using protein-engineering bioinformatics tools to select functionally-important residues based on correlated mutation analyses. A pair of amino acids in the periphery of PfPGI was found to be the dominant co-evolving mutation. The position of these selected residues was found to be non-obvious to conventional protein engineering methods. We designed a small smart library of variants by substituting the co-evolved pair and screened their biochemical activity, which revealed their functional relevance. Four mutants were further selected from the library for purification, measurement of their specific activity, crystal structure determination, and metal cofactor coordination analysis. Though the mutant structures and metal cofactor coordination were strikingly similar, variations in their activity correlated with their fine-tuned dynamics and solvent access regulation. Alternative, small smart libraries for enzyme optimization are suggested by our approach, which is able to identify non-obvious yet beneficial mutations.

Identification of fungal oxaloacetate hydrolyase within the isocitrate lyase/PEP mutase enzyme superfamily using a sequence marker-based method.

Aspergillus niger produces oxalic acid through the hydrolysis of oxaloacetate, catalyzed by the cytoplasmic enzyme oxaloacetate acetylhydrolase (OAH). The A. niger genome encodes four additional open reading frames with strong sequence similarity to OAH yet only the oahA gene encodes OAH activity. OAH and OAH-like proteins form subclass of the isocitrate lyase/PEP mutase enzyme superfamily, which is ubiquitous present filamentous fungi. Analysis of function-specific residues using a superfamily-based approach revealed an active site serine as a possible sequence marker for OAH activity. We propose that presence of this serine in family members correlates with presence of OAH activity whereas its absence correlates with absence of OAH. This hypothesis was tested by carrying out a serine mutagenesis study with the OAH from the fungal oxalic acid producer Botrytis cinerea and the OAH active plant petal death protein as test systems.

A quantitative indicator diagram for lytic polysaccharide monooxygenases reveals the role of aromatic surface residues in HjLPMO9A regioselectivity.

Lytic polysaccharide monooxygenases (LPMOs) have changed our understanding of lignocellulosic degradation dramatically over the last years. These metalloproteins catalyze oxidative cleavage of recalcitrant polysaccharides and can act on the C1 and/or C4 position of glycosidic bonds. Structural data have led to several hypotheses, but we are still a long way from reaching complete understanding of the factors that determine their divergent regioselectivity. Site-directed mutagenesis enables the investigation of structure-function relationship in enzymes and will be of major importance in unraveling this intriguing matter. In this context, it is crucial to have an enzyme assay or screening approach with a direct correlation with the desired functionality. LPMOs render this search extra challenging due to their insoluble substrates, complex pattern of reaction products and lack of synthetic standards of most oxidized products. Here, we describe a regioselectivity indicator diagram based on the time-course of only 2 HPAEC-PAD signals. The diagram was successfully used to confirm the hypothesis that aromatic surface residues influence the C1/C4 oxidation ratio in Hypocrea jecorina LPMO9A. Consequently, the diagram should become a valuable tool in the search towards better understanding and engineering of regioselectivity in LPMOs.

NewProt - a protein engineering portal.

The NewProt protein engineering portal is a one-stop-shop for in silico protein engineering. It gives access to a large number of servers that compute a wide variety of protein structure characteristics supporting work on the modification of proteins through the introduction of (multiple) point mutations. The results can be inspected through multiple visualizers. The HOPE software is included to indicate mutations with possible undesired side effects. The Hotspot Wizard software is embedded for the design of mutations that modify a proteins' activity, specificity, or stability. The NewProt portal is freely accessible at http://newprot.cmbi.umcn.nl/ and http://newprot.fluidops.net/.

Erratum: Enzyme discovery beyond homology: a unique hydroxynitrile lyase in the Bet v1 superfamily.

This corrects the article DOI: 10.1038/srep46738.