Standardized cloning vectors for protein production and generation of large gene libraries in Escherichia coli

Biocatalytic Conversion of Carrageenans for the Production of 3,6-Anhydro-D-galactose

Marine biomass stands out as a sustainable resource for generating value-added chemicals. In particular, anhydrosugars derived from carrageenans exhibit a variety of biological functions, rendering them highly promising for utilization and cascading in food, cosmetic, and biotechnological applications. However, the limitation of available sulfatases to break down the complex sulfation patterns of carrageenans poses a significant limitation for the sustainable production of valuable bioproducts from red algae. In this study, we screened several carrageenolytic polysaccharide utilization loci for novel sulfatase activities to assist the efficient conversion of a variety of sulfated galactans into the target product 3,6-anhydro-D-galactose. Inspired by the carrageenolytic pathways in marine heterotrophic bacteria, we systematically combined these novel sulfatases with other carrageenolytic enzymes, facilitating the development of the first enzymatic one-pot biotransformation of ι- and κ-carrageenan to 3,6-anhdyro-D-galactose. We further showed the applicability of this enzymatic bioconversion to a broad series of hybrid carrageenans, rendering this process a promising and sustainable approach for the production of value-added biomolecules from red-algal feedstocks.

Improving enzyme functional annotation by integrating in vitro and in silico approaches: The example of histidinol phosphate phosphatases

Advances in sequencing technologies have led to a rapid growth of public protein sequence databases, whereby the fraction of proteins with experimentally verified function continuously decreases. This problem is currently addressed by automated functional annotations with computational tools, which however lack the accuracy of experimental approaches and are susceptible to error propagation. Here, we present an approach that combines the efficiency of functional annotation by in silico methods with the rigor of enzyme characterization in vitro. First, a thorough experimental analysis of a representative enzyme of a group of homologues is performed which includes a focused alanine scan of the active site to determine a fingerprint of function-determining residues. In a second step, this fingerprint is used in combination with a sequence similarity network to identify putative isofunctional enzymes among the homologues. Using this approach in a proof-of-principle study, homologues of the histidinol phosphate phosphatase (HolPase) from Pseudomonas aeruginosa, many of which were annotated as phosphoserine phosphatases, were predicted to be HolPases. This functional annotation of the homologues was verified by in vitro testing of several representatives and an analysis of the occurrence of annotated HolPases in the corresponding phylogenetic groups. Moreover, the application of the same approach to the homologues of the HolPase from the archaeon Nitrosopumilus maritimus, which is not related to the HolPase from P. aeruginosa and was newly discovered in the course of this work, led to the annotation of the putative HolPase from various archaeal species.

Retracing the Rapid Evolution of an Herbicide-Degrading Enzyme by Protein Engineering

The mechanisms underlying the rapid evolution of novel enzymatic activities from promiscuous side activities are poorly understood. Recently emerged enzymes catalyzing the catabolic degradation of xenobiotic substances that have been spread out into the environment during the last decades provide an exquisite opportunity to study these mechanisms. A prominent example is the herbicide atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine), which is degraded through a number of enzymatic reactions constituting the Atz pathway. Here, we analyzed the evolution of the hydroxyatrazine ethylaminohydrolase AtzB, a Zn(II)-dependent metalloenzyme that adopts the amidohydrolase fold and catalyzes the second step of the Atz pathway. We searched for promiscuous side activities of AtzB, which might point to the identity of its progenitor. These investigations revealed that AtzB has low promiscuous guanine deaminase activity. Furthermore, we found that the two closest AtzB homologues, which have not been functionally annotated up to now, are guanine deaminases with modest promiscuous hydroxyatrazine hydrolase activity. Based on sequence comparisons with the closest AtzB homologues, the guanine deaminase activity of AtzB could be increased by three orders of magnitude through the introduction of only four active site mutations. Interestingly, introducing the inverse four mutations into the AtzB homologues significantly enhanced their hydroxyatrazine hydrolase activity, and in one case is even equivalent to that of wild-type AtzB. Molecular dynamics simulations elucidated the structural and molecular basis for the mutation-induced activity changes. The example of AtzB highlights how novel enzymes with high catalytic proficiency can evolve from low promiscuous side activities by only few mutational events within a short period of time.

Engineering of brick and staple components for ordered assembly of synthetic repeat proteins

Synthetic ɑRep repeat proteins are engineered as Brick and Staple protein pairs that together self-assemble into helical filaments. In most cases, the filaments spontaneously form supercrystals. Here, we describe an expanded series of ɑRep Bricks designed to stabilize the interaction between consecutive Bricks, to control the length of the assembled multimers, or to alter the spatial distribution of the Staple on the filaments. The effects of these Brick modifications on the assembly, on the final filament structure and on the crystal symmetry are analyzed by biochemical methods, electron microscopy and small angle X-ray scattering. We further extend the concept of Brick/Staple protein origami by designing a new type of "Janus"-like Brick protein that is equally assembled by orthogonal staples binding its inner or outer surfaces and thus ending inside or outside the filaments. The relative roles of longitudinal and lateral associations in the assembly process are discussed. This set of results demonstrates important proofs-of-principle for engineering these remarkably versatile proteins toward nanometer-to-micron scale constructions.

Investigating the Prevalence of RNA-Binding Metabolic Enzymes in E. coli

An open research field in cellular regulation is the assumed crosstalk between RNAs, metabolic enzymes, and metabolites, also known as the REM hypothesis. High-throughput assays have produced extensive interactome data with metabolic enzymes frequently found as hits, but only a few examples have been biochemically validated, with deficits especially in prokaryotes. Therefore, we rationally selected nineteen Escherichia coli enzymes from such datasets and examined their ability to bind RNAs using two complementary methods, iCLIP and SELEX. Found interactions were validated by EMSA and other methods. For most of the candidates, we observed no RNA binding (12/19) or a rather unspecific binding (5/19). Two of the candidates, namely glutamate-5-kinase (ProB) and quinone oxidoreductase (QorA), displayed specific and previously unknown binding to distinct RNAs. We concentrated on the interaction of QorA to the mRNA of yffO, a grounded prophage gene, which could be validated by EMSA and MST. Because the physiological function of both partners is not known, the biological relevance of this interaction remains elusive. Furthermore, we found novel RNA targets for the MS2 phage coat protein that served us as control. Our results indicate that RNA binding of metabolic enzymes in procaryotes is less frequent than suggested by the results of high-throughput studies, but does occur.

Structural and Functional Characterization of the Ureidoacrylate Amidohydrolase RutB from Escherichia coli

We present a detailed structure-function analysis of the ureidoacrylate amidohydrolase RutB from Eschericha coli, which is an essential enzyme of the Rut pathway for pyrimidine utilization. Crystals of selenomethionine-labeled RutB were produced, which allowed us to determine the first structure of the enzyme at a resolution of 1.9 Å and to identify it as a new member of the isochorismatase-like hydrolase family. RutB was co-crystallized with the substrate analogue ureidopropionate, revealing the mode of substrate binding. Mutation of residues constituting the catalytic triad (D24A, D24N, K133A, C166A, C166S, C166T, C166Y) resulted in complete inactivation of RutB, whereas mutation of other residues close to the active site (Y29F, Y35F, N72A, W74A, W74F, E80A, E80D, S92A, S92T, S92Y, Q105A, Y136A, Y136F) leads to distinct changes of the turnover number (k_cat) and/or the Michaelis constant (K_M). The results of our structural and mutational studies allowed us to assign specific functions to individual residues and to formulate a plausible reaction mechanism for RutB.

Experimental and computational analysis of the ancestry of an evolutionary young enzyme from histidine biosynthesis

The conservation of fold and chemistry of the enzymes associated with histidine biosynthesis suggests that this pathway evolved prior to the diversification of Bacteria, Archaea, and Eukaryotes. The only exception is the histidinol phosphate phosphatase (HolPase). So far, non-homologous HolPases that possess distinct folds and belong to three different protein superfamilies have been identified in various phylogenetic clades. However, their evolution has remained unknown to date. Here, we analyzed the evolutionary history of the HolPase from γ-Proteobacteria (HisB-N). It has been argued that HisB-N and its closest homologue d-glycero-d-manno-heptose-1,7-bisphosphate 7-phosphatase (GmhB) have emerged from the same promiscuous ancestral phosphatase. GmhB variants catalyze the hydrolysis of the anomeric d-glycero-d-manno-heptose-1,7-bisphosphate (αHBP or βHBP) with a strong preference for one anomer (αGmhB or βGmhB). We found that HisB-N from Escherichia coli shows promiscuous activity for βHBP but not αHBP, while βGmhB from Crassaminicella sp. shows promiscuous activity for HolP. Accordingly, a combined phylogenetic tree of αGmhBs, βGmhBs, and HisB-N sequences revealed that HisB-Ns form a compact subcluster derived from βGmhBs. Ancestral sequence reconstruction and in vitro analysis revealed a promiscuous HolPase activity in the resurrected enzymes prior to functional divergence of the successors. The following increase in catalytic efficiency of the HolP turnover is reflected in the shape and electrostatics of the active site predicted by AlphaFold. An analysis of the phylogenetic tree led to a revised evolutionary model that proposes the horizontal gene transfer of a promiscuous βGmhB from δ- to γ-Proteobacteria where it evolved to the modern HisB-N.

Hexamerization and thermostability emerged very early during geranylgeranylglyceryl phosphate synthase evolution

A large number of archaea live in hyperthermophilic environments. In consequence, their proteins need to adopt to these harsh conditions, including the enzymes that catalyze the synthesis of their membrane ether lipids. The enzyme that catalyzes the formation of the first ether bond in these lipids, geranylgeranylglyceryl phosphate synthase (GGGPS), exists as a hexamer in many hyperthermophilic archaea, and a recent study suggested that hexamerization serves for a fine‐tuning of the flexibility – stability trade‐off under hyperthermophilic conditions. We have recently reconstructed the sequences of ancestral group II GGGPS enzymes and now present a detailed biochemical characterization of nine of these predecessors, which allowed us to trace back the evolution of hexameric GGGPS and to draw conclusions about the properties of extant GGGPS branches that were not accessible to experiments up to now. Almost all ancestral GGGPS proteins formed hexamers, which demonstrates that hexamerization is even more widespread among the GGGPS family than previously assumed. Furthermore, all experimentally studied ancestral proteins showed high thermostability. Our results indicate that the hexameric oligomerization state and thermostability were present very early during the evolution of group II GGGPS, while the fine tuning of the flexibility – stability trade‐off developed very late, independent of the emergence of hexamerization.

Reprogramming the Specificity of a Protein Interface by Computational and Data-Driven Design

The formation of specific protein complexes in a cell is a non-trivial problem given the co-existence of thousands of different polypeptide chains. A particularly difficult case are two glutamine amidotransferase complexes (anthranilate synthase [AS] and aminodeoxychorismate synthase [ADCS]), which are composed of homologous pairs of synthase and glutaminase subunits. We have attempted to identify discriminating interface residues of the glutaminase subunit TrpG from AS, which are responsible for its specific interaction with the synthase subunit TrpEx and prevent binding to the closely related synthase subunit PabB from ADCS. For this purpose, TrpG-specific interface residues were grafted into the glutaminase subunit PabA from ADCS by two different approaches, namely a computational and a data-driven one. Both approaches resulted in PabA variants that bound TrpEx with higher affinity than PabB. Hence, we have accomplished a reprogramming of protein-protein interaction specificity that provides insights into the evolutionary adaptation of protein interfaces.

Analysis of allosteric communication in a multienzyme complex by ancestral sequence reconstruction

Tryptophan synthase (TS) is a heterotetrameric αββα complex. It is characterized by the channeling of the reaction intermediate indole and the mutual activation of the α-subunit TrpA and the β-subunit TrpB via a complex allosteric network. We have analyzed this allosteric network by means of ancestral sequence reconstruction (ASR), which is an in silico method to resurrect extinct ancestors of modern proteins. Previously, the sequences of TrpA and TrpB from the last bacterial common ancestor (LBCA) have been computed by means of ASR and characterized. LBCA-TS is similar to modern TS by forming a αββα complex with indole channeling taking place. However, LBCA-TrpA allosterically decreases the activity of LBCA-TrpB, whereas, for example, the modern ncTrpA from Neptuniibacter caesariensis allosterically increases the activity of ncTrpB. To identify amino acid residues that are responsible for this inversion of the allosteric effect, all 6 evolutionary TrpA and TrpB intermediates that stepwise link LBCA-TS with ncTS were characterized. Remarkably, the switching from TrpB inhibition to TrpB activation by TrpA occurred between 2 successive TS intermediates. Sequence comparison of these 2 intermediates and iterative rounds of site-directed mutagenesis allowed us to identify 4 of 413 residues from TrpB that are crucial for its allosteric activation by TrpA. The effect of our mutational studies was rationalized by a community analysis based on molecular dynamics simulations. Our findings demonstrate that ancestral sequence reconstruction can efficiently identify residues contributing to allosteric signal propagation in multienzyme complexes.

Library Selection with a Randomized Repertoire of (βα)8-Barrel Enzymes Results in Unexpected Induction of Gene Expression

The potential of the frequently encountered (βα)₈-barrel fold to acquire new functions was tested by an approach combining random mutagenesis and selection in vivo. For this purpose, the genes encoding 52 different phosphate-binding (βα)₈-barrel proteins were subjected to error-prone PCR and cloned into an expression plasmid. The resulting mixed repertoire was used to transform different auxotrophic Escherichia coli strains, each lacking an enzyme with a phosphate-containing substrate. After plating of the different transformants on minimal medium, growth was observed only for two strains, lacking either the gene for the serine phosphatase SerB or the phosphoserine aminotransferase SerC. The same mutants of the E. coli genes nanE (encoding a putative N-acetylmannosamine-6-phosphate 2-epimerase) and pdxJ (encoding the pyridoxine 5'-phosphate synthase) were responsible for rescuing both ΔserB and ΔserC. Unexpectedly, the complementing NanE and PdxJ variants did not catalyze the SerB or SerC reactions in vitro. Instead, RT-qPCR, RNAseq, and transcriptome analysis showed that they rescue the deletions by enlisting the help of endogenous E. coli enzymes HisB and HisC through exclusive up-regulation of histidine operon transcription. While the promiscuous SerB activity of HisB is well-established, our data indicate that HisC is promiscuous for the SerC reaction, as well. The successful rescue of ΔserB and ΔserC through point mutations and recruitment of additional amino acids in NanE and PdxJ provides another example for the adaptability of the (βα)₈-barrel fold.

Light-Regulation of Tryptophan Synthase by Combining Protein Design and Enzymology

The spatiotemporal control of enzymes by light is of growing importance for industrial biocatalysis. Within this context, the photo-control of allosteric interactions in enzyme complexes, common to practically all metabolic pathways, is particularly relevant. A prominent example of a metabolic complex with a high application potential is tryptophan synthase from Salmonella typhimurium (TS), in which the constituting TrpA and TrpB subunits mutually stimulate each other via a sophisticated allosteric network. To control TS allostery with light, we incorporated the unnatural amino acid o-nitrobenzyl-O-tyrosine (ONBY) at seven strategic positions of TrpA and TrpB. Initial screening experiments showed that ONBY in position 58 of TrpA (aL58ONBY) inhibits TS activity most effectively. Upon UV irradiation, ONBY decages to tyrosine, largely restoring the capacity of TS. Biochemical characterization, extensive steady-state enzyme kinetics, and titration studies uncovered the impact of aL58ONBY on the activities of TrpA and TrpB and identified reaction conditions under which the influence of ONBY decaging on allostery reaches its full potential. By applying those optimal conditions, we succeeded to directly light-activate TS(aL58ONBY) by a factor of ~100. Our findings show that rational protein design with a photo-sensitive unnatural amino acid combined with extensive enzymology is a powerful tool to fine-tune allosteric light-activation of a central metabolic enzyme complex.

Light Regulation of Enzyme Allostery through Photo-responsive Unnatural Amino Acids

Imidazole glycerol phosphate synthase (ImGPS) is an allosteric bienzyme complex in which substrate binding to the synthase subunit HisF stimulates the glutaminase subunit HisH. To control this stimulation with light, we have incorporated the photo-responsive unnatural amino acids phenylalanine-4'-azobenzene (AzoF), o-nitropiperonyl-O-tyrosine (NPY), and methyl-o-nitropiperonyllysine (mNPK) at strategic positions of HisF. The light-mediated isomerization of AzoF at position 55 (fS55AzoF^E ↔ fS55AzoF^Z) resulted in a reversible 10-fold regulation of HisH activity. The light-mediated decaging of NPY at position 39 (fY39NPY → fY39) and of mNPK at position 99 (fK99mNPK → fK99) led to a 4- to 6-fold increase of HisH activity. Molecular dynamics simulations explained how the unnatural amino acids interfere with the allosteric machinery of ImGPS and revealed additional aspects of HisH stimulation in wild-type ImGPS. Our findings show that unnatural amino acids can be used as a powerful tool for the spatiotemporal control of a central metabolic enzyme complex by light.

Functional characterisation of two Δ12-desaturases demonstrates targeted production of linoleic acid as pheromone precursor in Nasonia

Insect pheromones are often derived from fatty acid metabolism. Fatty acid desaturases, enzymes introducing double bonds into fatty acids, are crucial for the biosynthesis of these chemical signals. Δ12-desaturases catalyse the biosynthesis of linoleic acid by introducing a second double bond into oleic acid, but have been identified in only a few animal species. Here, we report the functional characterisation of two Δ12-desaturases, Nvit_D12a and Nvit_D12b, from the parasitic wasp Nasonia vitripennis. We demonstrate that Nvit_D12a is expressed in the rectal vesicle of males where they produce a linoleic acid-derived sex pheromone to attract virgin females. ¹³C-labelling experiments with Urolepis rufipes, a closely related species belonging to the 'Nasonia group', revealed that females, but not males, are able to synthesise linoleic acid. U. rufipes males produce an isoprenoid sex pheromone in the same gland and do not depend on linoleic acid for pheromone production. This suggests that Δ12-desaturases are common in the 'Nasonia group', but acquired a specialised function in chemical communication of those species that use linoleic acid as a pheromone precursor. Phylogenetic analysis suggests that insect Δ12-desaturases have evolved repeatedly from Δ9-desaturases in different insect taxa. Hence, insects have developed a way to produce linoleic acid independent of the omega desaturase subfamily which harbours all of the eukaryotic Δ12-desaturases known so far.

Sequence selection by FitSS4ASR alleviates ancestral sequence reconstruction as exemplified for geranylgeranylglyceryl phosphate synthase

For evolutionary studies, but also for protein engineering, ancestral sequence reconstruction (ASR) has become an indispensable tool. The first step of every ASR protocol is the preparation of a representative sequence set containing at most a few hundred recent homologs whose composition determines decisively the outcome of a reconstruction. A common approach for sequence selection consists of several rounds of manual recompilation that is driven by embedded phylogenetic analyses of the varied sequence sets. For ASR of a geranylgeranylglyceryl phosphate synthase, we additionally utilized FitSS4ASR, which replaces this time-consuming protocol with an efficient and more rational approach. FitSS4ASR applies orthogonal filters to a set of homologs to eliminate outlier sequences and those bearing only a weak phylogenetic signal. To demonstrate the usefulness of FitSS4ASR, we determined experimentally the oligomerization state of eight predecessors, which is a delicate and taxon-specific property. Corresponding ancestors deduced in a manual approach and by means of FitSS4ASR had the same dimeric or hexameric conformation; this concordance testifies to the efficiency of FitSS4ASR for sequence selection. FitSS4ASR-based results of two other ASR experiments were added to the Supporting Information. Program and documentation are available at https://gitlab.bioinf.ur.de/hek61586/FitSS4ASR.

Mapping key amino acid residues for the epimerase efficiency and stereospecificity of the sex pheromone biosynthetic short-chain dehydrogenases/reductases of Nasonia

Males of the parasitic wasp genus Nasonia use blends of chiral hydroxylactones as sex pheromones to attract conspecific females. Whereas all Nasonia species use a mixture of (4R,5S)-5-hydroxy-4-decanolide (RS) and 4-methylquinazoline (MQ) as sex pheromones, Nasonia vitripennis evolved (4R,5R)-5-hydroxy-4-decanolide (RR) as an extra sex pheromone component. We recently identified and functionally characterized three short-chain dehydrogenases/reductases (SDRs) NV10127, NV10128, and NV10129 that are capable of catalyzing the epimerization of RS to RR via (4R)-5-oxo-4-decanolide (ODL) as intermediate. Despite their very high sequence identities of 88-98%, these proteins differ drastically in their ability to epimerize RS to RR and in their stereoselectivity when reducing ODL to RR/RS. Here, in order to unravel the sequence differences underlying these varying functional properties of NV1027, NV10128 and NV10129, we created chimeras of the three enzymes and monitored their catalytic activities in vitro. The results show that a few amino acid changes at the C-termini and active sites of Nasonia vitripennis SDRs lead to substantially altered RS to RR epimerization and ODL-reduction activities. Thus, our study adds to the understanding of pheromone evolution by showing that subtle mutations in key biosynthetic enzymes can result in drastic effects on the composition of chemical signals.