Computer-Aided Protein Directed Evolution: a Review of Web Servers, Databases and other Computational Tools for Protein Engineering

Biotechnological Plastic Degradation and Valorization Using Systems Metabolic Engineering

Various kinds of plastics have been developed over the past century, vastly improving the quality of life. However, the indiscriminate production and irresponsible management of plastics have led to the accumulation of plastic waste, emerging as a pressing environmental concern. To establish a clean and sustainable plastic economy, plastic recycling becomes imperative to mitigate resource depletion and replace non-eco-friendly processes, such as incineration. Although chemical and mechanical recycling technologies exist, the prevalence of composite plastics in product manufacturing complicates recycling efforts. In recent years, the biodegradation of plastics using enzymes and microorganisms has been reported, opening a new possibility for biotechnological plastic degradation and bio-upcycling. This review provides an overview of microbial strains capable of degrading various plastics, highlighting key enzymes and their role. In addition, recent advances in plastic waste valorization technology based on systems metabolic engineering are explored in detail. Finally, future perspectives on systems metabolic engineering strategies to develop a circular plastic bioeconomy are discussed.

In silico prediction of mutation sites for anthranilate synthase from Serratia marcesens to deregulate tryptophan feedback inhibition

Allosteric feedback inhibition of the committed step in amino acid biosynthetic pathways is a major concern for production of amino acids at industrial scale. Anthranilate synthase (AS) catalyzes the first reaction of tryptophan biosynthetic pathway found in microorganisms and is feedback inhibited by its own product i.e. tryptophan. Here, we identified new mutant sites in AS using computational mutagenesis approach. MD simulations (20 ns) followed by MMPBSA and per residue decomposition energy analysis identified seven amino acid residues with best binding affinity for tryptophan. All 19 mutant structures were generated for each identified amino acid residue followed by simulation to evaluate effect of mutation on protein stability. Later, molecular docking studies were employed to generate mutant-tryptophan complex and structures with binding energies (kcal/mol) much higher than wild-type AS were selected. Finally, two mutants i.e., S37W and S37H were identified on the basis of positive binding scores and loss of tryptophan binding inside pocket. Further, MD simulations run for 200 ns were performed over these mutant-tryptophan complexes followed by RMSD, RMSF, radius of gyration , solvent accessible surface area , intra-protein hydrogen bond numbers, principal component analysis, free energy landscape (FEL) and secondary structure analysis to rationale effect of mutations on stability of protein. Cross correlation analysis of mutant site amino acids (S37W) with key residues of catalytic site (G325, T326, H395 and G482) was done to evaluate the effect of mutations on catalytic site conformation. Current computational mutagenesis approach predicted two mutants S37W and S37H with proposed deregulated feedback inhibition by tryptophan and retained catalytic activity.Communicated by Ramaswamy H. Sarma.

NMR-guided directed evolution

Directed evolution is a powerful tool for improving existing properties and imparting completely new functionalities to proteins1-4. Nonetheless, its potential in even small proteins is inherently limited by the astronomical number of possible amino acid sequences. Sampling the complete sequence space of a 100-residue protein would require testing of 20100 combinations, which is beyond any existing experimental approach. In practice, selective modification of relatively few residues is sufficient for efficient improvement, functional enhancement and repurposing of existing proteins5. Moreover, computational methods have been developed to predict the locations and, in certain cases, identities of potentially productive mutations6-9. Importantly, all current approaches for prediction of hot spots and productive mutations rely heavily on structural information and/or bioinformatics, which is not always available for proteins of interest. Moreover, they offer a limited ability to identify beneficial mutations far from the active site, even though such changes may markedly improve the catalytic properties of an enzyme10. Machine learning methods have recently showed promise in predicting productive mutations11, but they frequently require large, high-quality training datasets, which are difficult to obtain in directed evolution experiments. Here we show that mutagenic hot spots in enzymes can be identified using NMR spectroscopy. In a proof-of-concept study, we converted myoglobin, a non-enzymatic oxygen storage protein, into a highly efficient Kemp eliminase using only three mutations. The observed levels of catalytic efficiency exceed those of proteins designed using current approaches and are similar with those of natural enzymes for the reactions that they are evolved to catalyse. Given the simplicity of this experimental approach, which requires no a priori structural or bioinformatic knowledge, we expect it to be widely applicable and to enable the full potential of directed enzyme evolution.

Protein Design: From the Aspect of Water Solubility and Stability

Water solubility and structural stability are key merits for proteins defined by the primary sequence and 3D-conformation. Their manipulation represents important aspects of the protein design field that relies on the accurate placement of amino acids and molecular interactions, guided by underlying physiochemical principles. Emulated designer proteins with well-defined properties both fuel the knowledge-base for more precise computational design models and are used in various biomedical and nanotechnological applications. The continuous developments in protein science, increasing computing power, new algorithms, and characterization techniques provide sophisticated toolkits for solubility design beyond guess work. In this review, we summarize recent advances in the protein design field with respect to water solubility and structural stability. After introducing fundamental design rules, we discuss the transmembrane protein solubilization and de novo transmembrane protein design. Traditional strategies to enhance protein solubility and structural stability are introduced. The designs of stable protein complexes and high-order assemblies are covered. Computational methodologies behind these endeavors, including structure prediction programs, machine learning algorithms, and specialty software dedicated to the evaluation of protein solubility and aggregation, are discussed. The findings and opportunities for Cryo-EM are presented. This review provides an overview of significant progress and prospects in accurate protein design for solubility and stability.

Critical assessment of structure-based approaches to improve protein resistance in aqueous ionic liquids by enzyme-wide saturation mutagenesis

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for green industrial processes involving biocatalysts, but often reduce enzyme activity. Experimental and computational methods are applied to predict favorable substitution sites and, most often, subsequent site-directed surface charge modifications are introduced to enhance enzyme resistance towards aIL. However, almost no studies evaluate the prediction precision with random mutagenesis or the application of simple data-driven filtering processes. Here, we systematically and rigorously evaluated the performance of 22 previously described structure-based approaches to increase enzyme resistance to aIL based on an experimental complete site-saturation mutagenesis library of Bacillus subtilis Lipase A (BsLipA) screened against four aIL. We show that, surprisingly, most of the approaches yield low gain-in-precision (GiP) values, particularly for predicting relevant positions: 14 approaches perform worse than random mutagenesis. Encouragingly, exploiting experimental information on the thermostability of BsLipA or structural weak spots of BsLipA predicted by rigidity theory yields GiP = 3.03 and 2.39 for relevant variants and GiP = 1.61 and 1.41 for relevant positions. Combining five simple-to-compute physicochemical and evolutionary properties substantially increases the precision of predicting relevant variants and positions, yielding GiP = 3.35 and 1.29. Finally, combining these properties with predictions of structural weak spots identified by rigidity theory additionally improves GiP for relevant variants up to 4-fold to ∼10 and sustains or increases GiP for relevant positions, resulting in a prediction precision of ∼90% compared to ∼9% in random mutagenesis. This combination should be applicable to other enzyme systems for guiding protein engineering approaches towards improved aIL resistance.

Improving the catalytic efficiency and dimeric stability of Cu,Zn superoxide dismutase by combining structure-guided consensus approach with site-directed mutagenesis

Superoxide dismutase (SOD) leads the front line of defense against injuries mediated by the reactive oxygen species (ROS). The SOD from a high-altitude plant Potentilla atrosanguinea is a unique thermostable enzyme. In this study, we applied a structure-guided consensus approach on Cu,Zn SOD from Potentilla atrosanguinea plant, to improve its enzymatic properties. The polar uncharged amino acid (threonine) at position 97 of wild-type (WT) SOD was selected as a target residue for substitution by aspartate (T97D) through site-directed mutagenesis. The WT and T97D were examined by a combinative approach consisting of robust computational and experimental tools. The in-silico analysis indicated improved dimeric stability in T97D as compared to the WT. The strong interactions between the monomers were related to improved dimerization and enhanced catalytic efficiency of T97D. These results were validated by in-vitro assays showing improved dimer stability and catalytic efficiency in T97D than WT. Moreover, the mutation also improved the thermostability of the enzyme. The combined structural and functional data described the basis for improved specific activity and thermostability. This study could expand the scope of interface residue to be explored as targets for designing of SODs with improved kinetics.

Illustrative Tutorials for ProThermDB: Thermodynamic Database for Proteins and Mutants

ProThermDB (https://web.iitm.ac.in/bioinfo2/prothermdb/index.html) is a primary resource for protein stability, which contains experimentally determined thermodynamic data for proteins and their mutants. The most recent version of ProThermDB accumulates the data obtained from both high- and low-throughput experimental biophysical methods. It includes comprehensive information at four different levels, i.e.: (i) protein sequence and structure; (ii) experimental conditions; (iii) thermodynamic parameters such as Gibbs free energy, melting temperature, enthalpy, etc.; and (iv) literature. In the following protocols, we present detailed tutorials for retrieving data using different search, display and sorting options, interpretation of search results, description of each entry-level information category, data upload and download, cross-links with other databases, and visualization options. This protocol consists of six pictorial exercises, which are useful for biologists/users to understand the contents and organization of data in ProThermDB. Further, potential applications of ProThermDB in protein engineering are discussed. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Retrieval of experimental thermodynamic data for wild-type and mutants of a specific protein using a simple query Basic Protocol 2: Retrieval of stabilizing point mutations, which are located at the interior of α-helical regions, and obtaining data by thermal denaturation methods Basic Protocol 3: Retrieval of destabilizing point mutations, which are in β-sheets of exposed regions, and obtaining data by chemical denaturation methods (urea and GdnHCl) Basic Protocol 4: Retrieval of stabilizing and destabilizing point mutations in a range of physiological conditions (pH: 6-9 and T: 20°C-25°C) and publication years (2010-2020) Support Protocol: Downloading the entire data of the database for academic research purposes and submission of new data in ProThermDB.

Microbial proteases: ubiquitous enzymes with innumerable uses

Proteases are ubiquitous enzymes, having significant physiological roles in both synthesis and degradation. The use of microbial proteases in food fermentation is an age-old process, which is today being successfully employed in other industries with the advent of ‘omics’ era and innovations in genetic and protein engineering approaches. Proteases have found application in industries besides food, like leather, textiles, detergent, waste management, agriculture, animal husbandry, cosmetics, and pharmaceutics. With the rising demands and applications, researchers are exploring various approaches to discover, redesign, or artificially synthesize enzymes with better applicability in the industrial processes. These enzymes offer a sustainable and environmentally safer option, besides possessing economic and commercial value. Various bacterial and fungal proteases are already holding a commercially pivotal role in the industry. The current review summarizes the characteristics and types of proteases, microbial source, their current and prospective applications in various industries, and future challenges. Promoting these biocatalysts will prove significant in betterment of the modern world.

Information System for Selection of Conditions and Equipment for Mammalian Cell Cultivation

Over the past few decades, animal cell culture technology has advanced significantly. It is now considered a reliable, functional, and relatively well-developed technology. At present, biotherapeutic drugs are synthesized using cell culture techniques by large manufacturing enterprises that produce products for commercial use and clinical research. The reliable implementation of mammalian cell culture technology requires the optimization of a number of variables, including the culture environment and bioreactor conditions, suitable cell lines, operating costs, efficient process management and, most importantly, quality. Successful implementation also requires an appropriate process development strategy, industrial scale, and characteristics, as well as the certification of sustainable procedures that meet the requirements of current regulations. All of this has led to a trend of increasing research in the field of biotechnology and, as a result, to a great accumulation of scientific information which, however, remains fragmentary and non-systematic. The development of information and network technologies allow us to solve this problem. Information system creation allows for implementation of the modern concept of integrating various structured and unstructured data, as well as the collection of information from internal and external sources. We propose and develop an information system which contains the conditions and various parameters of cultivation processes. The associated ranking system is the result of the set of recommendations—both from technological and hardware solutions—which allow for choosing the optimal conditions for the cultivation of mammalian cells at the stage of scientific research, thereby significantly reducing the time and cost of work. The proposed information system allows for the accumulation of experience regarding existing technologies for the cultivation of mammalian cells, along with application to the development of new technologies. The main goal of the present work is to discuss information systems, the organizational support of scientific research in the field of mammalian cell cultivation, and to provide a detailed description of the developed system and its main modules, including the conceptual and logical scheme of the database.

Computational design of enzymes for biotechnological applications

Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.

Protein Engineering for Improving and Diversifying Natural Product Biosynthesis

Proteins found in nature have traditionally been the most frequently used biocatalysts to produce numerous natural products ranging from commodity chemicals to pharmaceuticals. Protein engineering has emerged as a powerful biotechnological toolbox in the development of metabolic engineering, particularly for the biosynthesis of natural products. Recently, protein engineering has become a favored method to improve enzymatic activity, increase enzyme stability, and expand product spectra in natural product biosynthesis. This review summarizes recent advances and typical strategies in protein engineering, highlighting the paramount role of protein engineering in improving and diversifying the biosynthesis of natural products. Future prospects and research directions are also discussed.

Rozhkova A, N. Zorov I, Korotkova O, P. Sinitsyn A, Schwaneberg U, D. Davari M. Engineering Robust Cellulases for Tailored Lignocellulosic Degradation Cocktails

Lignocellulosic biomass is a most promising feedstock in the production of second-generation biofuels. Efficient degradation of lignocellulosic biomass requires a synergistic action of several cellulases and hemicellulases. Cellulases depolymerize cellulose, the main polymer of the lignocellulosic biomass, to its building blocks. The production of cellulase cocktails has been widely explored, however, there are still some main challenges that enzymes need to overcome in order to develop a sustainable production of bioethanol. The main challenges include low activity, product inhibition, and the need to perform fine-tuning of a cellulase cocktail for each type of biomass. Protein engineering and directed evolution are powerful technologies to improve enzyme properties such as increased activity, decreased product inhibition, increased thermal stability, improved performance in non-conventional media, and pH stability, which will lead to a production of more efficient cocktails. In this review, we focus on recent advances in cellulase cocktail production, its current challenges, protein engineering as an efficient strategy to engineer cellulases, and our view on future prospects in the generation of tailored cellulases for biofuel production.

Robust Prediction of Single and Multiple Point Protein Mutations Stability Changes

Accurate prediction of protein stability changes resulting from amino acid substitutions is of utmost importance in medicine to better understand which mutations are deleterious, leading to diseases, and which are neutral. Since conducting wet lab experiments to get a better understanding of protein mutations is costly and time consuming, and because of huge number of possible mutations the need of computational methods that could accurately predict effects of amino acid mutations is of greatest importance. In this research, we present a robust methodology to predict the energy changes of a proteins upon mutations. The proposed prediction scheme is based on two step algorithm that is a Holdout Random Sampler followed by a neural network model for regression. The Holdout Random Sampler is utilized to analysis the energy change, the corresponding uncertainty, and to obtain a set of admissible energy changes, expressed as a cumulative distribution function. These values are further utilized to train a simple neural network model that can predict the energy changes. Results were blindly tested (validated) against experimental energy changes, giving Pearson correlation coefficients of 0.66 for Single Point Mutations and 0.77 for Multiple Point Mutations. These results confirm the successfulness of our method, since it outperforms majority of previous studies in this field.

Current Trends in Protein Engineering: Updates and Progress

Proteins are one of the most important and resourceful biomolecules that find applications in health, industry, medicine, research, and biotechnology. Given its tremendous relevance, protein engineering has emerged as significant biotechnological intervention in this area. Strategic utilization of protein engineering methods and approaches has enabled better enzymatic properties, better stability, increased catalytic activity and most importantly, interesting and wide range applicability of proteins. In fact, the commercialization of engineered proteins have manifested in economically beneficial and viable solutions for industry and healthcare sector. Protein engineering has also evolved to become a powerful tool contributing significantly to the developments in both synthetic biology and metabolic engineering. The present review revisits the current trends in protein engineering approaches such as rational design, directed evolution, de novo design, computational approaches etc. and encompasses the recent progresses made in this field over the last few years. The review also throws light on advanced or futuristic protein engineering aspects, which are being explored for design and development of novel proteins with improved properties or advanced applications.

Model-guided mechanism discovery and parameter selection for directed evolution

Directed evolution is frequently applied to identify genetic variants with improvements in a single or multiple properties. When used to improve multiple properties simultaneously, a common strategy is to apply alternating rounds of selection criteria to enrich for variants with each desirable trait. In particular, counterselection, or selection against undesired traits rather than for desired ones, has been successfully employed in many studies. Although the sequence and stringency of alternating selective pressures for different traits are known to be highly consequential for the outcome of the screen, the effects of these parameters have not been systematically evaluated. We developed a method for producing a statistical modeling framework to elucidate these effects. The model uses single-cell fluorescence intensity distributions to estimate the proportions of phenotypic populations within a library and then predicts the changes in these proportions depending on specified positive selective or counterselective pressures. We validated the approach using recently described systems for metabolite-responsive bacterial transcription factors and yeast G-protein-coupled receptors. Finally, we applied the model to identify biological sources that exert undesirable selective pressure on libraries during sorting. Notably, these pressures produce substantial artifacts that, if unaddressed, can lead to failure of the screen. This method for model generation can be applied to FACS-based directed evolution experiments to create a quantitative framework that identifies subtle population effects. Such models can guide the choice of experimental design parameters to better enrich for true positive genetic variants and improve the chance of successful directed evolution.

Computer-aided re-engineering of nonribosomal peptide and polyketide biosynthetic assembly lines

Covering: 2014 to 2019Nonribosomal peptide synthetases (NRPSs) and polyketide synthases (PKSs) have been the subject of engineering efforts for multiple decades. Their modular assembly line architecture potentially allows unlocking vast chemical space for biosynthesis. However, attempts thus far are often met with mixed success, due to limited molecular compatibility of the parts used for engineering. Now, new engineering strategies, increases in genomic data, and improved computational tools provide more opportunities for major progress. In this review we highlight some of the challenges and progressive strategies for the re-design of NRPSs & type I PKSs and survey useful computational tools and approaches to attain the ultimate goal of semi-automated and design-based engineering of novel peptide and polyketide products.

a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).

Structure-based design and application of an engineered glutathione transferase for the development of an optical biosensor for pesticides determination

In the present work, a structure-based design approach was used for the generation of a novel variant of synthetic glutathione transferase (PvGmGSTU) with higher sensitivity towards pesticides. Molecular modelling studies revealed Phe117 as a key residue that contributes to the formation of the hydrophobic binding site (H-site) and modulates the affinity of the enzyme towards xenobiotic compounds. Site-saturation mutagenesis of position Phe117 created a library of PvGmGSTU variants with altered kinetic and binding properties. Screening of the library against twenty-five different pesticides, showed that the mutant enzyme Phe117Ile displays 3-fold higher catalytic efficiency and exhibits increased affinity towards α-endosulfan, compared to the wild-type enzyme. Based on these catalytic features the mutant enzyme Phe117Ile was explored for the development of an optical biosensor for α-endosulfan. The enzyme was entrapped in alkosixylane sol-gel system in the presence of two pH indicators (bromocresol purple and phenol red). The sensing signal was based on the inhibition of the sol-gel entrapped GST, with subsequent decrease of released [H+] by the catalytic reaction, measured by sol-gel entrapped indicators. The assay response at 562 nm was linear in the range pH = 4-7. Linear calibration curves were obtained for α-endosulfan in the range of 0-30 μΜ. The reproducibility of the assay response, expressed by relative standard deviation, was in the order of 4.1% (N = 28). The method was successfully applied to the determination of α-endosulfan in real water samples without sample preparation steps.

Capturing variation impact on molecular interactions in the IMEx Consortium mutations data set

The current wealth of genomic variation data identified at nucleotide level presents the challenge of understanding by which mechanisms amino acid variation affects cellular processes. These effects may manifest as distinct phenotypic differences between individuals or result in the development of disease. Physical interactions between molecules are the linking steps underlying most, if not all, cellular processes. Understanding the effects that sequence variation has on a molecule’s interactions is a key step towards connecting mechanistic characterization of nonsynonymous variation to phenotype. We present an open access resource created over 14 years by IMEx database curators, featuring 28,000 annotations describing the effect of small sequence changes on physical protein interactions. We describe how this resource was built, the formats in which the data is provided and offer a descriptive analysis of the data set. The data set is publicly available through the IntAct website and is enhanced with every monthly release.

Genetic variants might exert their functional effects via influencing molecular interaction. Here, the authors present a resource featuring almost 28,000 annotations describing the effect of small sequence changes on physical protein interactions, curated by IMEx Consortium curators.

Histidine substitution in the most flexible fragments of firefly luciferase modifies its thermal stability

Molecular dynamics (MD) at two temperatures of 300 and 340 K identified two histidine residues, His461 and His489, in the most flexible regions of firefly luciferase, a light emitting enzyme. We therefore designed four protein mutants H461D, H489K, H489D and H489M to investigate their enzyme kinetic and thermodynamic stability changes. Substitution of His461 by aspartate (H461D) decreased ATP binding affinity, reduced the melting temperature of protein by around 25 °C and shifted its optimum temperature of activity to 10 °C. In line with the common feature of psychrophilic enzymes, the MD data showed that the overall flexibility of H461D was relatively high at low temperature, probably due to a decrease in the number of salt bridges around the mutation site. On the other hand, substitution of His489 by aspartate (H489D) introduced a new salt bridge between the C-terminal and N-terminal domains and increased protein rigidity but only slightly improved its thermal stability. Similar changes were observed for H489K and, to a lesser degree, H489M mutations. Based on our results we conclude that the MD simulation-based rational substitution of histidines by salt-bridge forming residues can modulate conformational dynamics in luciferase and shift its optimal temperature activity.

Modification of the peroxygenative:peroxidative activity ratio in the unspecific peroxygenase from Agrocybe aegerita by structure-guided evolution

Unspecific peroxygenase (UPO) is a heme-thiolate peroxidase capable of performing with high-selectivity C-H oxyfunctionalizations of great interest in organic synthesis through its peroxygenative activity. However, the convergence of such activity with an unwanted peroxidative activity encumbers practical applications. In this study, we have modified the peroxygenative:peroxidative activity ratio (P:p ratio) of UPO from Agrocybe aegerita by structure-guided evolution. Several flexible loops (Glu1-Pro35, Gly103-Asp131, Ser226-Gly243, Gln254-Thr276 and Ty293-Arg327) were selected on the basis on their B-factors and ΔΔG values. The full ensemble of segments (43% of UPO sequence) was subjected to focused evolution by the Mutagenic Organized Recombination Process by Homologous IN vivo Grouping (MORPHING) method in Saccharomyces cerevisiae. Five independent mutant libraries were screened in terms of P:p ratio and thermostability. We identified several variants that harbored substitutions at positions 120 and 320 with a strong enhancement in the P:p ratio albeit at the cost of stability. The most thermostable mutant of this process (S226G with an increased T50 of 2°C) was subjected to further combinatorial saturation mutagenesis on Thr120 and Thr320 yielding a collection of variants with modified P:p ratio and recovered stability. Our results seem to indicate the coexistence of several oxidation sites for peroxidative and peroxygenative activities in UPO.

Computational tools for enzyme improvement: why everyone can - and should - use them

This review presents computational methods that experimentalists can readily use to create smart libraries for enzyme engineering and to obtain insights into protein-substrate complexes. Computational tools have the reputation of being hard to use and inaccurate compared to experimental methods in enzyme engineering, yet they are essential to probe datasets of ever-increasing size and complexity. In recent years, bioinformatics groups have made a huge leap forward in providing user-friendly interfaces and accurate algorithms for experimentalists. These methods guide efficient experimental planning and allow the enzyme engineer to rationalize time and resources. Computational tools nevertheless face challenges in the realm of transient modern technology.

CADEE: Computer-Aided Directed Evolution of Enzymes

The tremendous interest in enzymes as biocatalysts has led to extensive work in enzyme engineering, as well as associated methodology development. Here, a new framework for computer-aided directed evolution of enzymes (CADEE) is presented which allows a drastic reduction in the time necessary to prepare and analyze in silico semi-automated directed evolution of enzymes. A pedagogical example of the application of CADEE to a real biological system is also presented in order to illustrate the CADEE workflow.

The Need for Integrated Approaches in Metabolic Engineering

This review highlights state-of-the-art procedures for heterologous small-molecule biosynthesis, the associated bottlenecks, and new strategies that have the potential to accelerate future accomplishments in metabolic engineering. We emphasize that a combination of different approaches over multiple time and size scales must be considered for successful pathway engineering in a heterologous host. We have classified these optimization procedures based on the "system" that is being manipulated: transcriptome, translatome, proteome, or reactome. By bridging multiple disciplines, including molecular biology, biochemistry, biophysics, and computational sciences, we can create an integral framework for the discovery and implementation of novel biosynthetic production routes.

Design and engineering of artificial oxygen-activating metalloenzymes

Many efforts are being made in the design and engineering of metalloenzymes with catalytic properties fulfilling the needs of practical applications. Progress in this field has recently been accelerated by advances in computational, molecular and structural biology. This review article focuses on the recent examples of oxygen-activating metalloenzymes, developed through the strategies of de novo design, miniaturization processes and protein redesign. Considerable progress in these diverse design approaches has produced many metal-containing biocatalysts able to adopt the functions of native enzymes or even novel functions beyond those found in Nature.

Beyond the outer limits of nature by directed evolution

For more than thirty years, biotechnology has borne witness to the power of directed evolution in designing molecules of industrial relevance. While scientists all over the world discuss the future of molecular evolution, dozens of laboratory-designed products are being released with improved characteristics in terms of turnover rates, substrate scope, catalytic promiscuity or stability. In this review we aim to present the most recent advances in this fascinating research field that are allowing us to surpass the limits of nature and apply newly gained attributes to a range of applications, from gene therapy to novel green processes. The use of directed evolution in non-natural environments, the generation of catalytic promiscuity for non-natural reactions, the insertion of unnatural amino acids into proteins or the creation of unnatural DNA, is described comprehensively, together with the potential applications in bioremediation, biomedicine and in the generation of new bionanomaterials. These successful case studies show us that the limits of directed evolution will be defined by our own imagination, and in some cases, stretching beyond that.

Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering

Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.

Synthetic biology for the directed evolution of protein biocatalysts: navigating sequence space intelligently

The amino acid sequence of a protein affects both its structure and its function. Thus, the ability to modify the sequence, and hence the structure and activity, of individual proteins in a systematic way, opens up many opportunities, both scientifically and (as we focus on here) for exploitation in biocatalysis. Modern methods of synthetic biology, whereby increasingly large sequences of DNA can be synthesised de novo, allow an unprecedented ability to engineer proteins with novel functions. However, the number of possible proteins is far too large to test individually, so we need means for navigating the 'search space' of possible protein sequences efficiently and reliably in order to find desirable activities and other properties. Enzymologists distinguish binding (Kd) and catalytic (kcat) steps. In a similar way, judicious strategies have blended design (for binding, specificity and active site modelling) with the more empirical methods of classical directed evolution (DE) for improving kcat (where natural evolution rarely seeks the highest values), especially with regard to residues distant from the active site and where the functional linkages underpinning enzyme dynamics are both unknown and hard to predict. Epistasis (where the 'best' amino acid at one site depends on that or those at others) is a notable feature of directed evolution. The aim of this review is to highlight some of the approaches that are being developed to allow us to use directed evolution to improve enzyme properties, often dramatically. We note that directed evolution differs in a number of ways from natural evolution, including in particular the available mechanisms and the likely selection pressures. Thus, we stress the opportunities afforded by techniques that enable one to map sequence to (structure and) activity in silico, as an effective means of modelling and exploring protein landscapes. Because known landscapes may be assessed and reasoned about as a whole, simultaneously, this offers opportunities for protein improvement not readily available to natural evolution on rapid timescales. Intelligent landscape navigation, informed by sequence-activity relationships and coupled to the emerging methods of synthetic biology, offers scope for the development of novel biocatalysts that are both highly active and robust.

Application of synthetic biology for production of chemicals in yeast Saccharomyces cerevisiae

Synthetic biology and metabolic engineering enable generation of novel cell factories that efficiently convert renewable feedstocks into biofuels, bulk, and fine chemicals, thus creating the basis for biosustainable economy independent on fossil resources. While over a hundred proof-of-concept chemicals have been made in yeast, only a very small fraction of those has reached commercial-scale production so far. The limiting factor is the high research cost associated with the development of a robust cell factory that can produce the desired chemical at high titer, rate, and yield. Synthetic biology has the potential to bring down this cost by improving our ability to predictably engineer biological systems. This review highlights synthetic biology applications for design, assembly, and optimization of non-native biochemical pathways in baker's yeast Saccharomyces cerevisiae We describe computational tools for the prediction of biochemical pathways, molecular biology methods for assembly of DNA parts into pathways, and for introducing the pathways into the host, and finally approaches for optimizing performance of the introduced pathways.

Directed evolution 2.0: improving and deciphering enzyme properties

A KnowVolution: knowledge gaining directed evolution including four phases is proposed in this feature article, which generates improved enzyme variants and molecular understanding.

SDS, a structural disruption score for assessment of missense variant deleteriousness

We have developed a novel structure-based evaluation for missense variants that explicitly models protein structure and amino acid properties to predict the likelihood that a variant disrupts protein function. A structural disruption score (SDS) is introduced as a measure to depict the likelihood that a case variant is functional. The score is constructed using characteristics that distinguish between causal and neutral variants within a group of proteins. The SDS score is correlated with standard sequence-based deleteriousness, but shows promise for improving discrimination between neutral and causal variants at less conserved sites. The prediction was performed on 3-dimentional structures of 57 gene products whose homozygous SNPs were identified as case-exclusive variants in an exome sequencing study of epilepsy disorders. We contrasted the candidate epilepsy variants with scores for likely benign variants found in the EVS database, and for positive control variants in the same genes that are suspected to promote a range of diseases. To derive a characteristic profile of damaging SNPs, we transformed continuous scores into categorical variables based on the score distribution of each measurement, collected from all possible SNPs in this protein set, where extreme measures were assumed to be deleterious. A second epilepsy dataset was used to replicate the findings. Causal variants tend to receive higher sequence-based deleterious scores, induce larger physico-chemical changes between amino acid pairs, locate in protein domains, buried sites or on conserved protein surface clusters, and cause protein destabilization, relative to negative controls. These measures were agglomerated for each variant. A list of nine high-priority putative functional variants for epilepsy was generated. Our newly developed SDS protocol facilitates SNP prioritization for experimental validation.

Computational tools for designing smart libraries

Traditional directed evolution experiments are often time-, labor- and cost-intensive because they involve repeated rounds of random mutagenesis and the selection or screening of large mutant libraries. The efficiency of directed evolution experiments can be significantly improved by targeting mutagenesis to a limited number of hot-spot positions and/or selecting a limited set of substitutions. The design of such "smart" libraries can be greatly facilitated by in silico analyses and predictions. Here we provide an overview of computational tools applicable for (a) the identification of hot-spots for engineering enzyme properties, and (b) the evaluation of predicted hot-spots and selection of suitable amino acids for substitutions. The selected tools do not require any specific expertise and can easily be implemented by the wider scientific community.

Engineering proteins for thermostability through rigidifying flexible sites

Engineering proteins for thermostability is an exciting and challenging field since it is critical for broadening the industrial use of recombinant proteins. Thermostability of proteins arises from the simultaneous effect of several forces such as hydrophobic interactions, disulfide bonds, salt bridges and hydrogen bonds. All of these interactions lead to decreased flexibility of polypeptide chain. Structural studies of mesophilic and thermophilic proteins showed that the latter need more rigid structures to compensate for increased thermal fluctuations. Hence flexibility can be an indicator to pinpoint weak spots for enhancing thermostability of enzymes. A strategy has been proven effective in enhancing proteins' thermostability with two steps: predict flexible sites of proteins firstly and then rigidify these sites. We refer to this approach as rigidify flexible sites (RFS) and give an overview of such a method through summarizing the methods to predict flexibility of a protein, the methods to rigidify residues with high flexibility and successful cases regarding enhancing thermostability of proteins using RFS.

A frequency-based linguistic approach to protein decoding and design: Simple concepts, diverse applications, and the SCS Package

Protein structure and function information is coded in amino acid sequences. However, the relationship between primary sequences and three-dimensional structures and functions remains enigmatic. Our approach to this fundamental biochemistry problem is based on the frequencies of short constituent sequences (SCSs) or words. A protein amino acid sequence is considered analogous to an English sentence, where SCSs are equivalent to words. Availability scores, which are defined as real SCS frequencies in the non-redundant amino acid database relative to their probabilistically expected frequencies, demonstrate the biological usage bias of SCSs. As a result, this frequency-based linguistic approach is expected to have diverse applications, such as secondary structure specifications by structure-specific SCSs and immunological adjuvants with rare or non-existent SCSs. Linguistic similarities (e.g., wide ranges of scale-free distributions) and dissimilarities (e.g., behaviors of low-rank samples) between proteins and the natural English language have been revealed in the rank-frequency relationships of SCSs or words. We have developed a web server, the SCS Package, which contains five applications for analyzing protein sequences based on the linguistic concept. These tools have the potential to assist researchers in deciphering structurally and functionally important protein sites, species-specific sequences, and functional relationships between SCSs. The SCS Package also provides researchers with a tool to construct amino acid sequences de novo based on the idiomatic usage of SCSs.