Ancestral Haloalkane Dehalogenases Show Robustness and Unique Substrate Specificity

SUPREM: an engineered non-site-specific m6A RNA methyltransferase with highly improved efficiency

N 6-Methyladenine (m6A) RNA methylation plays a key role in RNA processing and translational regulation, influencing both normal physiological and pathological processes. Yet, current techniques for studying RNA methylation struggle to isolate the effects of individual m6A modifications. Engineering of RNA methyltransferases (RNA MTases) could enable development of improved synthetic biology tools to manipulate RNA methylation, but it is challenging due to limited understanding of structure-function relationships in RNA MTases. Herein, using ancestral sequence reconstruction, we explore the sequence space of the bacterial DNA methyltransferase EcoGII (M.EcoGII), a promising target for protein engineering due to its lack of sequence specificity and its residual activity on RNA. We thereby created an efficient non-specific RNA MTase termed SUPer RNA EcoGII Methyltransferase (SUPREM), which exhibits 8-fold higher expression levels, 7°C higher thermostability and 12-fold greater m6A RNA methylation activity compared with M.EcoGII. Immunofluorescent staining and quantitative liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis confirmed SUPREM's higher RNA methylation activity compared with M.EcoGII in mammalian cells. Additionally, Nanopore direct RNA sequencing highlighted that SUPREM is capable of methylating a larger number of RNA methylation sites than M.EcoGII. Through phylogenetic and mutational analysis, we identified a critical residue for the enhanced RNA methylation activity of SUPREM. Collectively, our findings indicate that SUPREM holds promise as a versatile tool for in vivo RNA methylation and labeling.

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

Engineering an Ancestral Glycosyltransferase for Biosynthesis of 2-Phenylethyl-β-d-Glucopyranoside and Salidroside

Phenylethanoid glycosides (PhGs) are naturally occurring glycosides derived from plants with various biological activities. Glycosyltransferases catalyze the production of PhGs from phenylethanols via a transglycosylation reaction. The low activity and stability of glycosyltransferase limit its industrial application. An ancestral glycosyltransferase, UGTAn85, with heat resistance, alkali resistance, and high stability was resurrected using ancestral sequence reconstruction technology. This enzyme can efficiently convert phenylethanols to PhGs. The optimal reaction temperature and pH for UGTAn85 were found to be 70 °C and pH 10.0, respectively. This study employed a combination of structure-guided rational design and co-evolution analysis to enhance its catalytic activity. Potential mutation sites were identified through computer-aided design, including homology modeling, molecular docking, Rosetta dock design, molecular dynamics simulation, and co-evolution analysis. By targeted mutagenesis, the UGTAn85 mutant Q23E/N65D exhibited a 2.2-fold increase in enzyme activity (11.85 U/mg) and elevated affinity (Km = 0.11 mM) for 2-phenylethanol compared to UGTAn85. Following a fed-batch reaction, 36.16 g/L 2-phenylethyl-β-d-glucopyranoside and 51.49 g/L salidroside could be produced within 24 h, respectively. The findings in this study provide a new perspective on enhancing the stability and activity of glycosyltransferases, as well as a potential biocatalyst for the industrial production of PhGs.

Engineering Biocatalysts for the C-H Activation of Fatty Acids by Ancestral Sequence Reconstruction

Selective, one-step C-H activation of fatty acids from biomass is an attractive concept in sustainable chemistry. Biocatalysis has shown promise for generating high-value hydroxy acids, but to date enzyme discovery has relied on laborious screening and produced limited hits, which predominantly oxidise the subterminal positions of fatty acids. Herein we show that ancestral sequence reconstruction (ASR) is an effective tool to explore the sequence-activity landscape of a family of multidomain, self-sufficient P450 monooxygenases. We resurrected 11 catalytically active CYP116B ancestors, each with a unique regioselectivity fingerprint that varied from subterminal in the older ancestors to mid-chain in the lineage leading to the extant, P450-TT. In lineages leading to extant enzymes in thermophiles, thermostability increased from ancestral to extant forms, as expected if thermophily had arisen de novo. Our studies show that ASR can be applied to multidomain enzymes to develop active, self-sufficient monooxygenases as regioselective biocatalysts for fatty acid hydroxylation.

Evolution and separation of actinobacterial pyranose and C-glycoside-3-oxidases

FAD-dependent pyranose oxidase (POx) and C-glycoside-3-oxidase (CGOx) are both members of the glucose-methanol-choline superfamily of oxidoreductases and belong to the same sequence space. Pyranose oxidases had been studied for their oxidation of monosaccharides such as D-glucose, but recently, a bacterial C-glycoside-3-oxidase that is phylogenetically related to POx and that reacts with C-glycosides such as carminic acid, mangiferin or puerarin has been described. Since these actinobacterial CGOx enzymes belong to the same sequence space as bacterial POx, they must have evolved from the same ancestor. Here, we performed a phylogenetic analysis of actinobacterial sequences and resurrected seven ancestral enzymes of the POx/CGOx sequence space to study the evolutionary trajectory of substrate preferences for monosaccharides and C-glycosides. Clade I, with its dimeric member POx from Kitasatospora aureofaciens, shows strict preference for monosaccharides (D-glucose and D-xylose) and does not react with any of the glycosides tested. No extant member of clade II has been studied to date. The two extant members of clades III and IV, monomeric POx/CGOx from Pseudoarthrobacter siccitolerans and Streptomyces canus, oxidized both monosaccharides as well as various C-glycosides (homoorientin, isovitexin, mangiferin, and puerarin). Steady-state kinetic parameters of several clades III and IV ancestral enzymes indicate that the generalist ancestor N35 slowly evolved to present-day enzymes with a much higher preference for C-glycosides than monosaccharides. Based on structural predictions of ancestors, we hypothesize that the strict specificity of bacterial clade I POx (and also fungal POx) is the result of oligomerization, which in turn results from the evolution of protein segments that were shown to be important for oligomerization, the arm, and the head domain.IMPORTANCEC-Glycosides often form active compounds in various plants. Breakage of the C-C bond in these glycosides to release the aglycone is challenging and proceeds via a two-step reaction, the oxidation of the sugar and subsequent cleavage of the C-C bond. Recently, an enzyme from a soil bacterium, FAD-dependent C-glycoside-3-oxidase (CGOx), was shown to catalyze the initial oxidation reaction. Here, we show that CGOx belongs to the same sequence space as pyranose oxidase (POx), and that an actinobacterial ancestor of the POx/CGOx family evolved into four clades, two of which show a high preference for C-glycosides.

Machine Learning-Guided Protein Engineering

Recent progress in engineering highly promising biocatalysts has increasingly involved machine learning methods. These methods leverage existing experimental and simulation data to aid in the discovery and annotation of promising enzymes, as well as in suggesting beneficial mutations for improving known targets. The field of machine learning for protein engineering is gathering steam, driven by recent success stories and notable progress in other areas. It already encompasses ambitious tasks such as understanding and predicting protein structure and function, catalytic efficiency, enantioselectivity, protein dynamics, stability, solubility, aggregation, and more. Nonetheless, the field is still evolving, with many challenges to overcome and questions to address. In this Perspective, we provide an overview of ongoing trends in this domain, highlight recent case studies, and examine the current limitations of machine learning-based methods. We emphasize the crucial importance of thorough experimental validation of emerging models before their use for rational protein design. We present our opinions on the fundamental problems and outline the potential directions for future research.

Multimeric structure of a subfamily III haloalkane dehalogenase-like enzyme solved by combination of cryo-EM and x-ray crystallography

Haloalkane dehalogenase (HLD) enzymes employ an SN 2 nucleophilic substitution mechanism to erase halogen substituents in diverse organohalogen compounds. Subfamily I and II HLDs are well-characterized enzymes, but the mode and purpose of multimerization of subfamily III HLDs are unknown. Here we probe the structural organization of DhmeA, a subfamily III HLD-like enzyme from the archaeon Haloferax mediterranei, by combining cryo-electron microscopy (cryo-EM) and x-ray crystallography. We show that full-length wild-type DhmeA forms diverse quaternary structures, ranging from small oligomers to large supramolecular ring-like assemblies of various sizes and symmetries. We optimized sample preparation steps, enabling three-dimensional reconstructions of an oligomeric species by single-particle cryo-EM. Moreover, we engineered a crystallizable mutant (DhmeAΔGG ) that provided diffraction-quality crystals. The 3.3 Å crystal structure reveals that DhmeAΔGG forms a ring-like 20-mer structure with outer and inner diameter of ~200 and ~80 Å, respectively. An enzyme homodimer represents a basic repeating building unit of the crystallographic ring. Three assembly interfaces (dimerization, tetramerization, and multimerization) were identified to form the supramolecular ring that displays a negatively charged exterior, while its interior part harboring catalytic sites is positively charged. Localization and exposure of catalytic machineries suggest a possible processing of large negatively charged macromolecular substrates.

Ion-Powered Rotary Motors: Where Did They Come from and Where They Are Going?

Molecular motors are found in many living organisms. One such molecular machine, the ion-powered rotary motor (IRM), requires the movement of ions across a membrane against a concentration gradient to drive rotational movement. The bacterial flagellar motor (BFM) is an example of an IRM which relies on ion movement through the stator proteins to generate the rotation of the flagella. There are many ions which can be used by the BFM stators to power motility and different ions can be used by a single bacterium expressing multiple stator variants. The use of ancestral sequence reconstruction (ASR) and functional analysis of reconstructed stators shows promise for understanding how these proteins evolved and when the divergence in ion use may have occurred. In this review, we discuss extant BFM stators and the ions that power them as well as recent examples of the use of ASR to study ion-channel selectivity and how this might be applied to further study of the BFM stator complex.

Efficient Exploration of Sequence Space by Sequence-Guided Protein Engineering and Design

The rapid growth of sequence databases over the past two decades means that protein engineers faced with optimizing a protein for any given task will often have immediate access to a vast number of related protein sequences. These sequences encode information about the evolutionary history of the protein and the underlying sequence requirements to produce folded, stable, and functional protein variants. Methods that can take advantage of this information are an increasingly important part of the protein engineering tool kit. In this Perspective, we discuss the utility of sequence data in protein engineering and design, focusing on recent advances in three main areas: the use of ancestral sequence reconstruction as an engineering tool to generate thermostable and multifunctional proteins, the use of sequence data to guide engineering of multipoint mutants by structure-based computational protein design, and the use of unlabeled sequence data for unsupervised and semisupervised machine learning, allowing the generation of diverse and functional protein sequences in unexplored regions of sequence space. Altogether, these methods enable the rapid exploration of sequence space within regions enriched with functional proteins and therefore have great potential for accelerating the engineering of stable, functional, and diverse proteins for industrial and biomedical applications.

Engineering the Thermostability of a d-Carbamoylase Based on Ancestral Sequence Reconstruction for the Efficient Synthesis of d-Tryptophan

Employing ancestral sequence reconstruction and consensus sequence analysis, the thermostability of a novel d-carbamoylase derived from Nitratireductor indicus (NiHyuC) was engineered through greedy-oriented iterative combinatorial mutagenesis. A mutant S202P/E208D/R277L (M4Th3) was obtained with significantly elevated thermostability. M4Th3 has a half-life of 36.5 h at 40 °C, about 28.5 times of 1.3 h of its parent M4. For the reaction at 40 °C, M4Th3 can catalyze 10 mM N-carbamoyl-d-tryptophan to produce d-tryptophan with a conversion ratio of 96.4% after 12 h, which is significantly higher than 64.1% of M4. MD simulation reveals that new hydrogen bonds emerging from E208D on the surface can increase the hydrophobicity of the protein, leading to improved stability. More importantly, R277L could contribute to enhanced interface stability of homodimeric M4. This study provides a thermostable d-carbamoylase for the "hydantoinase process", which has potential in the industrial synthesis of optically pure natural and non-natural amino acids.

Ancestral sequence reconstruction as a tool to study the evolution of wood decaying fungi

The study of evolution is limited by the techniques available to do so. Aside from the use of the fossil record, molecular phylogenetics can provide a detailed characterization of evolutionary histories using genes, genomes and proteins. However, these tools provide scarce biochemical information of the organisms and systems of interest and are therefore very limited when they come to explain protein evolution. In the past decade, this limitation has been overcome by the development of ancestral sequence reconstruction (ASR) methods. ASR allows the subsequent resurrection in the laboratory of inferred proteins from now extinct organisms, becoming an outstanding tool to study enzyme evolution. Here we review the recent advances in ASR methods and their application to study fungal evolution, with special focus on wood-decay fungi as essential organisms in the global carbon cycling.

Engineering functional thermostable proteins using ancestral sequence reconstruction

Natural proteins are often only slightly more stable in the native state than the denatured state, and an increase in environmental temperature can easily shift the balance toward unfolding. Therefore, the engineering of proteins to improve protein stability is an area of intensive research. Thermostable proteins are required to withstand industrial process conditions, for increased shelf-life of protein therapeutics, for developing robust 'biobricks' for synthetic biology applications, and for research purposes (e.g., structure determination). In addition, thermostability buffers the often destabilizing effects of mutations introduced to improve other properties. Rational design approaches to engineering thermostability require structural information, but even with advanced computational methods, it is challenging to predict or parameterize all the relevant structural factors with sufficient precision to anticipate the results of a given mutation. Directed evolution is an alternative when structures are unavailable but requires extensive screening of mutant libraries. Recently, however, bioinspired approaches based on phylogenetic analyses have shown great promise. Leveraging the rapid expansion in sequence data and bioinformatic tools, ancestral sequence reconstruction can generate highly stable folds for novel applications in industrial chemistry, medicine, and synthetic biology. This review provides an overview of the factors important for successful inference of thermostable proteins by ancestral sequence reconstruction and what it can reveal about the determinants of stability in proteins.

Tuning Enzyme Thermostability via Computationally Guided Covalent Stapling and Structural Basis of Enhanced Stabilization

Enhancing the thermostability of enzymes without impacting their catalytic function represents an important yet challenging goal in protein engineering and biocatalysis. We recently introduced a novel method for enzyme thermostabilization that relies on the computationally guided installation of genetically encoded thioether "staples" into a protein via cysteine alkylation with the noncanonical amino acid O-2-bromoethyl tyrosine (O2beY). Here, we demonstrate the functionality of an expanded set of electrophilic amino acids featuring chloroacetamido, acrylamido, and vinylsulfonamido side-chain groups for protein stapling using this strategy. Using a myoglobin-based cyclopropanase as a model enzyme, our studies show that covalent stapling with p-chloroacetamido-phenylalanine (pCaaF) provides higher stapling efficiency and enhanced stability (thermodynamic and kinetic) compared to the other stapled variants and the parent protein. Interestingly, molecular simulations of conformational flexibility of the cross-links show that the pCaaF staple allows fewer energetically feasible conformers than the other staples, and this property may be a broader indicator of stability enhancement. Using this strategy, pCaaF-stapled variants with significantly enhanced stability against thermal denaturation (ΔTm' = +27 °C) and temperature-induced heme loss (ΔT50 = +30 °C) were obtained while maintaining high levels of catalytic activity and stereoselectivity. Crystallographic analyses of singly and doubly stapled variants provide key insights into the structural basis for stabilization, which includes both direct interactions of the staples with protein residues and indirect interactions through adjacent residues involved in heme binding. This work expands the toolbox of protein stapling strategies available for protein stabilization.

Structure-function relationship between soluble epoxide hydrolases structure and their tunnel network

Enzymes with buried active sites maintain their catalytic function via a single tunnel or tunnel network. In this study we analyzed the functionality of soluble epoxide hydrolases (sEHs) tunnel network, by comparing the overall enzyme structure with the tunnel's shape and size. sEHs were divided into three groups based on their structure and the tunnel usage. The obtained results were compared with known substrate preferences of the studied enzymes, as well as reported in our other work evolutionary analyses data. The tunnel network architecture corresponded well with the evolutionary lineage of the source organism and large differences between enzymes were observed from long fragments insertions. This strategy can be used during protein re-engineering process for large changes introduction, whereas tunnel modification can be applied for fine-tuning of enzyme.

Combining Ancestral Reconstruction with Folding-Landscape Simulations to Engineer Heterologous Protein Expression

Obligate symbionts typically exhibit high evolutionary rates. Consequently, their proteins may differ considerably from their modern and ancestral homologs in terms of both sequence and properties, thus providing excellent models to study protein evolution. Also, obligate symbionts are challenging to culture in the lab and proteins from uncultured organisms must be produced in heterologous hosts using recombinant DNA technology. Obligate symbionts thus replicate a fundamental scenario of metagenomics studies aimed at the functional characterization and biotechnological exploitation of proteins from the bacteria in soil. Here, we use the thioredoxin from Candidatus Photodesmus katoptron, an uncultured symbiont of flashlight fish, to explore evolutionary and engineering aspects of protein folding in heterologous hosts. The symbiont protein is a standard thioredoxin in terms of 3D-structure, stability and redox activity. However, its folding outside the original host is severely impaired, as shown by a very slow refolding in vitro and an inefficient expression in E. coli that leads mostly to insoluble protein. By contrast, resurrected Precambrian thioredoxins express efficiently in E. coli, plausibly reflecting an ancient adaptation to unassisted folding. We have used a statistical-mechanical model of the folding landscape to guide back-to-ancestor engineering of the symbiont protein. Remarkably, we find that the efficiency of heterologous expression correlates with the in vitro (i.e., unassisted) folding rate and that the ancestral expression efficiency can be achieved with only 1-2 back-to-ancestor replacements. These results demonstrate a minimal-perturbation, sequence-engineering approach to rescue inefficient heterologous expression which may potentially be useful in metagenomics efforts targeting recent adaptations.

Structural Analysis of the Ancestral Haloalkane Dehalogenase AncLinB-DmbA

Haloalkane dehalogenases (EC 3.8.1.5) play an important role in hydrolytic degradation of halogenated compounds, resulting in a halide ion, a proton, and an alcohol. They are used in biocatalysis, bioremediation, and biosensing of environmental pollutants and also for molecular tagging in cell biology. The method of ancestral sequence reconstruction leads to prediction of sequences of ancestral enzymes allowing their experimental characterization. Based on the sequences of modern haloalkane dehalogenases from the subfamily II, the most common ancestor of thoroughly characterized enzymes LinB from Sphingobium japonicum UT26 and DmbA from Mycobacterium bovis 5033/66 was in silico predicted, recombinantly produced and structurally characterized. The ancestral enzyme AncLinB-DmbA was crystallized using the sitting-drop vapor-diffusion method, yielding rod-like crystals that diffracted X-rays to 1.5 Å resolution. Structural comparison of AncLinB-DmbA with their closely related descendants LinB and DmbA revealed some differences in overall structure and tunnel architecture. Newly prepared AncLinB-DmbA has the highest active site cavity volume and the biggest entrance radius on the main tunnel in comparison to descendant enzymes. Ancestral sequence reconstruction is a powerful technique to study molecular evolution and design robust proteins for enzyme technologies.

Highly thermostable and promiscuous β-1,3-xylanasen designed by optimized ancestral sequence reconstruction

The ancestor of β-1,3-xylanases (AncXyl09) were reconstructed by the optimized ancestral sequences reconstruction strategy to solve the poor catalytic performances of existing β-1,3-xylanases. The results showed that the half-life at 50 °C was 65.08 h, indicating good thermostability. The large number of hydrogen bonds and the disulfide bonds were the major attributes related with the thermal stability of Anxyl09. Interestingly, AncXyl09 could hydrolyze lichen besides the original substrate of β-1, 3-xylan, which is the first reported β-1,3-xylanase with substrate promiscuity. Moreover, the hydrolytic products are mainly disaccharides, the content of β-1,3-xylobiose and lichoridiose more than 70% as determined by high performance liquid chromatography (HPLC), which could significantly facilitate the separation and purification of oligosaccharides. The successful design of AncXyl09 was the representative of the semi-rationally engineered β-1, 3-xylanase, which will shield a new light on the β-1,3-xylanase engineering, active oligosaccharide preparation and marine algae resource utilization.

Fully Automated Ancestral Sequence Reconstruction using FireProtASR

Protein evolution and protein engineering techniques are of great interest in basic science and industrial applications such as pharmacology, medicine, or biotechnology. Ancestral sequence reconstruction (ASR) is a powerful technique for probing evolutionary relationships and engineering robust proteins with good thermostability and broad substrate specificity. The following protocol describes the setting up and execution of an automated FireProt^ASR workflow using a dedicated web site. The service allows for inference of ancestral proteins automatically, from a single protein sequence. Once a protein sequence is submitted, the server will build a dataset of homology sequences, perform a multiple sequence alignment (MSA), build a phylogenetic tree, and reconstruct ancestral nodes. The protocol is also highly flexible and allows for multiple forms of input, advanced settings, and the ability to start jobs from: (i) a single sequence, (ii) a set of homologous sequences, (iii) an MSA, and (iv) a phylogenetic tree. This approach automates all necessary steps and offers a way for novices with limited exposure to ASR techniques to improve the properties of a protein of interest. The technique can even be used to introduce catalytic promiscuity into an enzyme. A web server for accessing the fully automated workflow is freely accessible at https://loschmidt.chemi.muni.cz/fireprotasr/. © 2021 Wiley Periodicals LLC. Basic Protocol: ASR using the Web Server FireProt^ASR.

Ancestral sequence reconstruction for protein engineers

In addition to its value in the study of molecular evolution, ancestral sequence reconstruction (ASR) has emerged as a useful methodology for engineering proteins with enhanced properties. Proteins generated by ASR often exhibit unique or improved activity, stability, and/or promiscuity, all of which are properties that are valued by protein engineers. Comparison between extant proteins and evolutionary intermediates generated by ASR also allows protein engineers to identify substitutions that have contributed to functional innovation or diversification within protein families. As ASR becomes more widely adopted as a protein engineering approach, it is important to understand the applications, limitations, and recent developments of this technique. This review highlights recent exemplifications of ASR, as well as technical aspects of the reconstruction process that are relevant to protein engineering.

A Reconstructed Common Ancestor of the Fatty Acid Photo-decarboxylase Clade Shows Photo-decarboxylation Activity and Increased Thermostability

Light-dependent enzymes are a rare type of biocatalyst with high potential for research and biotechnology. A recently discovered fatty acid photo-decarboxylase from Chlorella variabilis NC64A (CvFAP) converts fatty acids to the corresponding hydrocarbons only when irradiated with blue light (400 to 520 nm). To expand the available catalytic diversity for fatty acid decarboxylation, we reconstructed possible ancestral decarboxylases from a set of 12 extant sequences that were classified under the fatty acid decarboxylases clade within the glucose-methanol choline (GMC) oxidoreductase family. One of the resurrected enzymes (ANC1) showed activity in the decarboxylation of fatty acids, showing that the clade indeed contains several photo-decarboxylases. ANC1 has a 15 °C higher melting temperature (Tm ) than the extant CvFAP. Its production yielded 12-fold more protein than this wild type decarboxylase, which offers practical advantages for the biochemical investigation of this photoenzyme. Homology modelling revealed amino acid substitutions to more hydrophilic residues at the surface and shorter flexible loops compared to the wild type. Using ancestral sequence reconstruction, we have expanded the existing pool of confirmed fatty acid photo-decarboxylases, providing access to a more robust catalyst for further development via directed evolution.

Engineering of Ancestors as a Tool to Elucidate Structure, Mechanism, and Specificity of Extant Terpene Cyclase

Structural information is crucial for understanding catalytic mechanisms and to guide enzyme engineering efforts of biocatalysts, such as terpene cyclases. However, low sequence similarity can impede homology modeling, and inherent protein instability presents challenges for structural studies. We hypothesized that X-ray crystallography of engineered thermostable ancestral enzymes can enable access to reliable homology models of extant biocatalysts. We have applied this concept in concert with molecular modeling and enzymatic assays to understand the structure activity relationship of spiroviolene synthase, a class I terpene cyclase, aiming to engineer its specificity. Engineering a surface patch in the reconstructed ancestor afforded a template structure for generation of a high-confidence homology model of the extant enzyme. On the basis of structural considerations, we designed and crystallized ancestral variants with single residue exchanges that exhibited tailored substrate specificity and preserved thermostability. We show how the two single amino acid alterations identified in the ancestral scaffold can be transferred to the extant enzyme, conferring a specificity switch that impacts the extant enzyme's specificity for formation of the diterpene spiroviolene over formation of sesquiterpenes hedycaryol and farnesol by up to 25-fold. This study emphasizes the value of ancestral sequence reconstruction combined with enzyme engineering as a versatile tool in chemical biology.

FireProtASR: A Web Server for Fully Automated Ancestral Sequence Reconstruction

There is a great interest in increasing proteins’ stability to widen their usability in numerous biomedical and biotechnological applications. However, native proteins cannot usually withstand the harsh industrial environment, since they are evolved to function under mild conditions. Ancestral sequence reconstruction is a well-established method for deducing the evolutionary history of genes. Besides its applicability to discover the most probable evolutionary ancestors of the modern proteins, ancestral sequence reconstruction has proven to be a useful approach for the design of highly stable proteins. Recently, several computational tools were developed, which make the ancestral reconstruction algorithms accessible to the community, while leaving the most crucial steps of the preparation of the input data on users’ side. FireProt^ASR aims to overcome this obstacle by constructing a fully automated workflow, allowing even the unexperienced users to obtain ancestral sequences based on a sequence query as the only input. FireProt^ASR is complemented with an interactive, easy-to-use web interface and is freely available at https://loschmidt.chemi.muni.cz/fireprotasr/.

Resurrection of ancestral effector caspases identifies novel networks for evolution of substrate specificity

Apoptotic caspases evolved with metazoans more than 950 million years ago (MYA), and a series of gene duplications resulted in two subfamilies consisting of initiator and effector caspases. The effector caspase genes (caspases-3, -6, and -7) were subsequently fixed into the Chordata phylum more than 650 MYA when the gene for a common ancestor (CA) duplicated, and the three effector caspases have persisted throughout mammalian evolution. All caspases prefer an aspartate residue at the P1 position of substrates, so each caspase evolved discrete cellular roles through changes in substrate recognition at the P4 position combined with allosteric regulation. We examined the evolution of substrate specificity in caspase-6, which prefers valine at the P4 residue, compared with caspases-3 and -7, which prefer aspartate, by reconstructing the CA of effector caspases (AncCP-Ef1) and the CA of caspase-6 (AncCP-6An). We show that AncCP-Ef1 is a promiscuous enzyme with little distinction between Asp, Val, or Leu at P4. The specificity of caspase-6 was defined early in its evolution, where AncCP-6An demonstrates a preference for Val over Asp at P4. Structures of AncCP-Ef1 and of AncCP-6An show a network of charged amino acids near the S4 pocket that, when combined with repositioning a flexible active site loop, resulted in a more hydrophobic binding pocket in AncCP-6An. The ancestral protein reconstructions show that the caspase-hemoglobinase fold has been conserved for over 650 million years and that only three substitutions in the scaffold are necessary to shift substrate selection toward Val over Asp.

Ancestral Resurrection and Directed Evolution of Fungal Mesozoic Laccases

Ancestral sequence reconstruction and resurrection provides useful information for protein engineering, yet its alliance with directed evolution has been little explored. In this study, we have resurrected several ancestral nodes of fungal laccases dating back ∼500 to 250 million years. Unlike modern laccases, the resurrected Mesozoic laccases were readily secreted by yeast, with similar kinetic parameters, a broader stability, and distinct pH activity profiles. The resurrected Agaricomycetes laccase carried 136 ancestral mutations, a molecular testimony to its origin, and it was subjected to directed evolution in order to improve the rate of 1,3-cyclopentanedione oxidation, a β-diketone initiator commonly used in vinyl polymerization reactions.IMPORTANCE The broad variety of biotechnological uses of fungal laccases is beyond doubt (food, textiles, pulp and paper, pharma, biofuels, cosmetics, and bioremediation), and protein engineering (in particular, directed evolution) has become the key driver for adaptation of these enzymes to harsh industrial conditions. Usually, the first requirement for directed laccase evolution is heterologous expression, which presents an important hurdle and often a time-consuming process. In this work, we resurrected a fungal Mesozoic laccase node which showed strikingly high heterologous expression and pH stability. As a proof of concept that the ancestral laccase is a suitable blueprint for engineering, we performed a quick directed evolution campaign geared to the oxidation of the β-diketone 1,3-cyclopentanedione, a poor laccase substrate that is used in the polymerization of vinyl monomers.

Structures of hyperstable ancestral haloalkane dehalogenases show restricted conformational dynamics

Ancestral sequence reconstruction is a powerful method for inferring ancestors of modern enzymes and for studying structure-function relationships of enzymes. We have previously applied this approach to haloalkane dehalogenases (HLDs) from the subfamily HLD-II and obtained thermodynamically highly stabilized enzymes (ΔT m up to 24 °C), showing improved catalytic properties. Here we combined crystallographic structural analysis and computational molecular dynamics simulations to gain insight into the mechanisms by which ancestral HLDs became more robust enzymes with novel catalytic properties. Reconstructed ancestors exhibited similar structure topology as their descendants with the exception of a few loop deviations. Strikingly, molecular dynamics simulations revealed restricted conformational dynamics of ancestral enzymes, which prefer a single state, in contrast to modern enzymes adopting two different conformational states. The restricted dynamics can potentially be linked to their exceptional stabilization. The study provides molecular insights into protein stabilization due to ancestral sequence reconstruction, which is becoming a widely used approach for obtaining robust protein catalysts.

Exploration of enzyme diversity: High-throughput techniques for protein production and microscale biochemical characterization

Enzymes are being increasingly utilized for acceleration of industrially and pharmaceutically critical chemical reactions. The strong demand for finding robust and efficient biocatalysts for these applications can be satisfied via the exploration of enzyme diversity. The first strategy is to mine the natural diversity, represented by millions of sequences available in the public genomic databases, by using computational approaches. Alternatively, metagenomic libraries can be targeted experimentally or computationally to explore the natural diversity of a specific environment. The second strategy, known as directed evolution, is to generate man-made diversity in the laboratory using gene mutagenesis and screen the constructed library of mutants. The selected hits must be experimentally characterized in both strategies, which currently represent the rate-limiting step in the process of diversity exploration. The traditional techniques used for biochemical characterization are time-demanding, cost, and sample volume ineffective, and low-throughput. Therefore, the development and implementation of high-throughput experimental methods are essential for discovering novel enzymes. This chapter describes the experimental protocols employing the combination of robust production and high-throughput microscale biochemical characterization of enzyme variants. We validated its applicability against the model enzyme family of haloalkane dehalogenases. These protocols can be adapted to other enzyme families, paving the way towards the functional characterization and quick identification of novel biocatalysts.

Practically useful protein-design methods combining phylogenetic and atomistic calculations

Our ability to design new or improved biomolecular activities depends on understanding the sequence-function relationships in proteins. The large size and fold complexity of most proteins, however, obscure these relationships, and protein-optimization methods continue to rely on laborious experimental iterations. Recently, a deeper understanding of the roles of stability-threshold effects and biomolecular epistasis in proteins has led to the development of hybrid methods that combine phylogenetic analysis with atomistic design calculations. These methods enable reliable and even single-step optimization of protein stability, expressibility, and activity in proteins that were considered outside the scope of computational design. Furthermore, ancestral-sequence reconstruction produces insights on missing links in the evolution of enzymes and binders that may be used in protein design. Through the combination of phylogenetic and atomistic calculations, the long-standing goal of general computational methods that can be universally applied to study and optimize proteins finally seems within reach.

Highly thermostable carboxylic acid reductases generated by ancestral sequence reconstruction

Carboxylic acid reductases (CARs) are biocatalysts of industrial importance. Their properties, especially their poor stability, render them sub-optimal for use in a bioindustrial pipeline. Here, we employed ancestral sequence reconstruction (ASR) - a burgeoning engineering tool that can identify stabilizing but enzymatically neutral mutations throughout a protein. We used a three-algorithm approach to reconstruct functional ancestors of the Mycobacterial and Nocardial CAR1 orthologues. Ancestral CARs (AncCARs) were confirmed to be CAR enzymes with a preference for aromatic carboxylic acids. Ancestors also showed varied tolerances to solvents, pH and in vivo-like salt concentrations. Compared to well-studied extant CARs, AncCARs had a Tm up to 35 °C higher, with half-lives up to nine times longer than the greatest previously observed. Using ancestral reconstruction we have expanded the existing CAR toolbox with three new thermostable CAR enzymes, providing access to the high temperature biosynthesis of aldehydes to drive new applications in biocatalysis.

Bacillus subtilis Spore Surface Display of Haloalkane Dehalogenase DhaA

The haloalkane dehalogenase DhaA can degrade sulfur mustard (2,2'-dichlorethyl sulfide; also known by its military designation HD) in a rapid and environmentally safe manner. However, DhaA is sensitive to temperature and pH, which limits its applications in natural or harsh environments. Spore surface display technology using resistant spores as a carrier to ensure enzymatic activity can reduce production costs and extend the range of applications of DhaA. To this end, we cloned recombinant Bacillus subtilis spores pHY300PLK-cotg-dhaa-6his/DB104(FH01) for the delivery of DhaA from Rhodococcus rhodochrous NCIMB 13064. A dot blotting showed that the fusion protein CotG-linker-DhaA accounted for 0.41% ± 0.03% (P < 0.01) of total spore coat proteins. Immunofluorescence analyses confirmed that DhaA was displayed on the spore surface. The hydrolyzing activity of DhaA displayed on spores towards the HD analog 2-chloroethyl ethylsulfide was 1.74 ± 0.06 U/mL (P < 0.01), with a specific activity was 0.34 ± 0.04 U/mg (P < 0.01). This is the first demonstration that DhaA displayed on the surface of B. subtilis spores retains enzymatic activity, which suggests that it can be used effectively in real-world applications including bioremediation of contaminated environments.

Protein engineers turned evolutionists-the quest for the optimal starting point

The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.

Biotechnological and protein-engineering implications of ancestral protein resurrection

Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.