Computational Design of Stable and Soluble Biocatalysts

Partial Consensus Design and Enhancement of Protein Function by Secondary-Structure-Guided Consensus Mutations

Consensus design (CD) is a representative sequence-based protein design method that enables the design of highly functional proteins by analyzing vast amounts of protein sequence data. This study proposes a partial consensus design (PCD) of a protein as a derivative approach of CD. The method replaces the target protein sequence with a consensus sequence in a secondary-structure-dependent manner (i.e., regionally dependent and divided into α-helix, β-sheet, and loop regions). In this study, we generated several artificial partial consensus l-threonine 3-dehydrogenases (PcTDHs) by PCD using the TDH from Cupriavidus necator (CnTDH) as a target protein. Structural and functional analysis of PcTDHs suggested that thermostability would be independently improved when consensus mutations are introduced into the loop region of TDHs. On the other hand, enzyme kinetic parameters (k_cat/K_m) and average productivity would be synergistically enhanced by changing the combination of the mutations-replacement of one region of CnTDH with a consensus sequence provided only negative effects, but the negative effects were nullified when the two regions were replaced simultaneously. Taken together, we propose the hypothesis that there are protein regions that encode individual protein properties, such as thermostability and activity, and that the introduction of consensus mutations into these regions could additively or synergistically modify their functions.

Design and construction of novel biocatalyst for bioprocessing: Recent advances and future outlook

Bioprocess, a biocatalysis-based technology, is becoming popular in many research fields and widely applied in industrial manufacturing. However, low bioconversion, low productivity, and high costs during industrial processes are usually the limitation in bioprocess. Therefore, many biocatalyst strategies have been developed to meet these challenges in recent years. In this review, we firstly discuss protein engineering strategies, which are emerged for improving the biocatalysis activity of biocatalysts. Then, we summarize metabolic engineering strategies that are promoting the development of microbial cell factories. Next, we illustrate the necessity of using the combining strategy of protein engineering and metabolic engineering for efficient biocatalysts. Lastly, future perspectives about the development and application of novel biocatalyst strategies are discussed. This review provides theoretical guidance for the development of efficient, sustainable, and economical bioprocesses mediated by novel biocatalysts.

Reviewing Challenges of Predicting Protein Melting Temperature Change Upon Mutation Through the Full Analysis of a Highly Detailed Dataset with High-Resolution Structures

Predicting the effects of mutations on protein stability is a key problem in fundamental and applied biology, still unsolved even for the relatively simple case of small, soluble, globular, monomeric, two-state-folder proteins. Many articles discuss the limitations of prediction methods and of the datasets used to train them, which result in low reliability for actual applications despite globally capturing trends. Here, we review these and other issues by analyzing one of the most detailed, carefully curated datasets of melting temperature change (ΔTm) upon mutation for proteins with high-resolution structures. After examining the composition of this dataset to discuss imbalances and biases, we inspect several of its entries assisted by an online app for data navigation and structure display and aided by a neural network that predicts ΔTm with accuracy close to that of programs available to this end. We pose that the ΔTm predictions of our network, and also likely those of other programs, account only for a baseline-like general effect of each type of amino acid substitution which then requires substantial corrections to reproduce the actual stability changes. The corrections are very different for each specific case and arise from fine structural details which are not well represented in the dataset and which, despite appearing reasonable upon visual inspection of the structures, are hard to encode and parametrize. Based on these observations, additional analyses, and a review of recent literature, we propose recommendations for developers of stability prediction methods and for efforts aimed at improving the datasets used for training. We leave our interactive interface for analysis available online at http://lucianoabriata.altervista.org/papersdata/proteinstability2021/s1626navigation.html so that users can further explore the dataset and baseline predictions, possibly serving as a tool useful in the context of structural biology and protein biotechnology research and as material for education in protein biophysics.

Recent advances in user-friendly computational tools to engineer protein function

Progress in technology and algorithms throughout the past decade has transformed the field of protein design and engineering. Computational approaches have become well-engrained in the processes of tailoring proteins for various biotechnological applications. Many tools and methods are developed and upgraded each year to satisfy the increasing demands and challenges of protein engineering. To help protein engineers and bioinformaticians navigate this emerging wave of dedicated software, we have critically evaluated recent additions to the toolbox regarding their application for semi-rational and rational protein engineering. These newly developed tools identify and prioritize hotspots and analyze the effects of mutations for a variety of properties, comprising ligand binding, protein-protein and protein-nucleic acid interactions, and electrostatic potential. We also discuss notable progress to target elusive protein dynamics and associated properties like ligand-transport processes and allosteric communication. Finally, we discuss several challenges these tools face and provide our perspectives on the further development of readily applicable methods to guide protein engineering efforts.

Implicit Solvents for the Polarizable Atomic Multipole AMOEBA Force Field

Computational protein design, ab initio protein/RNA folding, and protein-ligand screening can be too computationally demanding for explicit treatment of solvent. For these applications, implicit solvent offers a compelling alternative, which we describe here for the polarizable atomic multipole AMOEBA force field based on three treatments of continuum electrostatics: numerical solutions to the nonlinear and linearized versions of the Poisson-Boltzmann equation (PBE), the domain-decomposition conductor-like screening model (ddCOSMO) approximation to the PBE, and the analytic generalized Kirkwood (GK) approximation. The continuum electrostatics models are combined with a nonpolar estimator based on novel cavitation and dispersion terms. Electrostatic model parameters are numerically optimized using a least-squares style target function based on a library of 103 small-molecule solvation free energy differences. Mean signed errors for the adaptive Poisson-Boltzmann solver (APBS), ddCOSMO, and GK models are 0.05, 0.00, and 0.00 kcal/mol, respectively, while the mean unsigned errors are 0.70, 0.63, and 0.58 kcal/mol, respectively. Validation of the electrostatic response of the resulting implicit solvents, which are available in the Tinker (or Tinker-HP), OpenMM, and Force Field X software packages, is based on comparisons to explicit solvent simulations for a series of proteins and nucleic acids. Overall, the emergence of performative implicit solvent models for polarizable force fields opens the door to their use for folding and design applications.

Web-based tools for computational enzyme design

Enzymes are in high demand for very diverse biotechnological applications. However, natural biocatalysts often need to be engineered for fine-tuning their properties towards the end applications, such as the activity, selectivity, stability to temperature or co-solvents, and solubility. Computational methods are increasingly used in this task, providing predictions that narrow down the space of possible mutations significantly and can enormously reduce the experimental burden. Many computational tools are available as web-based platforms, making them accessible to non-expert users. These platforms are typically user-friendly, contain walk-throughs, and do not require deep expertise and installations. Here we describe some of the most recent outstanding web-tools for enzyme engineering and formulate future perspectives in this field.

Computational design of highly stable and soluble alcohol dehydrogenase for NADPH regeneration

Nicotinamide adenine dinucleotide phosphate (NADPH), as a well-known cofactor, is widely used in the most of enzymatic redox reactions, playing an important role in industrial catalysis. However, the absence of a comparable method for efficient NADP+ to NADPH cofactor regeneration radically impairs efficient green chemical synthesis. Alcohol dehydrogenase (ADH) enzymes, allowing the in situ regeneration of the redox cofactor NADPH with high specific activity and easy by-product separation process, are provided with great industrial application potential and research attention. Accordingly, herein a NADP+-specific ADH from Clostridium beijerinckii was selected to be engineered for cofactor recycle, using an automated algorithm named Protein Repair One-stop Shop (PROSS). The mutant CbADH-6M (S24P/G182A/G196A/H222D/S250E/S254R) exhibited a favorable soluble and highly active expression with an activity of 46.3 U/mL, which was 16 times higher than the wild type (2.9 U/mL), and a more stable protein conformation with an enhanced thermal stability: Δ T 1 / 2 60 min =  + 3.6 °C (temperature of 50% inactivation after incubation for 60 min). Furthermore, the activity of CbADH-6M was up-graded to 2401.8 U/mL by high cell density fermentation strategy using recombinant Escherichia coli, demonstrating its industrial potential. Finally, the superb efficiency for NADPH regeneration of the mutant enzyme was testified in the synthesis of some fine chiral aromatic alcohols coupling with another ADH from Lactobacillus kefir (LkADH).

Power of Biocatalysis for Organic Synthesis

Biocatalysis, using defined enzymes for organic transformations, has become a common tool in organic synthesis, which is also frequently applied in industry. The generally high activity and outstanding stereo-, regio-, and chemoselectivity observed in many biotransformations are the result of a precise control of the reaction in the active site of the biocatalyst. This control is achieved by exact positioning of the reagents relative to each other in a fine-tuned 3D environment, by specific activating interactions between reagents and the protein, and by subtle movements of the catalyst. Enzyme engineering enables one to adapt the catalyst to the desired reaction and process. A well-filled biocatalytic toolbox is ready to be used for various reactions. Providing nonnatural reagents and conditions and evolving biocatalysts enables one to play with the myriad of options for creating novel transformations and thereby opening new, short pathways to desired target molecules. Combining several biocatalysts in one pot to perform several reactions concurrently increases the efficiency of biocatalysis even further.

FireProtDB: database of manually curated protein stability data

The majority of naturally occurring proteins have evolved to function under mild conditions inside the living organisms. One of the critical obstacles for the use of proteins in biotechnological applications is their insufficient stability at elevated temperatures or in the presence of salts. Since experimental screening for stabilizing mutations is typically laborious and expensive, in silico predictors are often used for narrowing down the mutational landscape. The recent advances in machine learning and artificial intelligence further facilitate the development of such computational tools. However, the accuracy of these predictors strongly depends on the quality and amount of data used for training and testing, which have often been reported as the current bottleneck of the approach. To address this problem, we present a novel database of experimental thermostability data for single-point mutants FireProt^DB. The database combines the published datasets, data extracted manually from the recent literature, and the data collected in our laboratory. Its user interface is designed to facilitate both types of the expected use: (i) the interactive explorations of individual entries on the level of a protein or mutation and (ii) the construction of highly customized and machine learning-friendly datasets using advanced searching and filtering. The database is freely available at https://loschmidt.chemi.muni.cz/fireprotdb.

SoluProt: Prediction of Soluble Protein Expression in Escherichia coli

Motivation

Poor protein solubility hinders the production of many therapeutic and industrially useful proteins. Experimental efforts to increase solubility are plagued by low success rates and often reduce biological activity. Computational prediction of protein expressibility and solubility in Escherichia coli using only sequence information could reduce the cost of experimental studies by enabling prioritization of highly soluble proteins.

Results

A new tool for sequence-based prediction of soluble protein expression in E.coli, SoluProt, was created using the gradient boosting machine technique with the TargetTrack database as a training set. When evaluated against a balanced independent test set derived from the NESG database, SoluProt’s accuracy of 58.5% and AUC of 0.62 exceeded those of a suite of alternative solubility prediction tools. There is also evidence that it could significantly increase the success rate of experimental protein studies. SoluProt is freely available as a standalone program and a user-friendly webserver at https://loschmidt.chemi.muni.cz/soluprot/.

Availability and implementation

https://loschmidt.chemi.muni.cz/soluprot/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Computational design of noncanonical amino acid-based thioether staples at N/C-terminal domains of multi-modular pullulanase for thermostabilization in enzyme catalysis

Enzyme thermostabilization is considered a critical and often obligatory step in biosynthesis, because thermostability is a significant property of enzymes that can be used to evaluate their feasibility for industrial applications. However, conventional strategies for thermostabilizing enzymes generally introduce non-covalent interactions and/or natural covalent bonds caused by natural amino acid substitutions, and the trade-off between the activity and stability of enzymes remains a challenge. Here, we developed a computationally guided strategy for constructing thioether staples by incorporating noncanonical amino acid (ncAA) into the more flexible N/C-terminal domains of the multi-modular pullulanase from Bacillus thermoleovorans (BtPul) to enhance its thermostability. First, potential thioether staples located in the N/C-terminal domains of BtPul were predicted using RosettaMatch. Next, eight variants involving stable thioether staples were precisely predicted using FoldX and Rosetta ddg_monomer. Six positive variants were obtained, of which T73(O2beY)-171C had a 157% longer half-life at 70 °C and an increase of 7.0 °C in T m, when compared with the wild-type (WT). T73(O2beY)-171C/T126F/A72R exhibited an even more improved thermostability, with a 211% increase in half-life at 70 °C and a 44% enhancement in enzyme activity compared with the WT, which was attributed to further optimization of the local interaction network. This work introduces and validates an efficient strategy for enhancing the thermostability and activity of multi-modular enzymes.

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/.

EnzymeMiner: automated mining of soluble enzymes with diverse structures, catalytic properties and stabilities

Millions of protein sequences are being discovered at an incredible pace, representing an inexhaustible source of biocatalysts. Despite genomic databases growing exponentially, classical biochemical characterization techniques are time-demanding, cost-ineffective and low-throughput. Therefore, computational methods are being developed to explore the unmapped sequence space efficiently. Selection of putative enzymes for biochemical characterization based on rational and robust analysis of all available sequences remains an unsolved problem. To address this challenge, we have developed EnzymeMiner—a web server for automated screening and annotation of diverse family members that enables selection of hits for wet-lab experiments. EnzymeMiner prioritizes sequences that are more likely to preserve the catalytic activity and are heterologously expressible in a soluble form in Escherichia coli. The solubility prediction employs the in-house SoluProt predictor developed using machine learning. EnzymeMiner reduces the time devoted to data gathering, multi-step analysis, sequence prioritization and selection from days to hours. The successful use case for the haloalkane dehalogenase family is described in a comprehensive tutorial available on the EnzymeMiner web page. EnzymeMiner is a universal tool applicable to any enzyme family that provides an interactive and easy-to-use web interface freely available at https://loschmidt.chemi.muni.cz/enzymeminer/.

Decoding the intricate network of molecular interactions of a hyperstable engineered biocatalyst

Computational design of protein catalysts with enhanced stabilities for use in research and enzyme technologies is a challenging task. Using force-field calculations and phylogenetic analysis, we previously designed the haloalkane dehalogenase DhaA115 which contains 11 mutations that confer upon it outstanding thermostability (T m = 73.5 °C; ΔT m > 23 °C). An understanding of the structural basis of this hyperstabilization is required in order to develop computer algorithms and predictive tools. Here, we report X-ray structures of DhaA115 at 1.55 Å and 1.6 Å resolutions and their molecular dynamics trajectories, which unravel the intricate network of interactions that reinforce the αβα-sandwich architecture. Unexpectedly, mutations toward bulky aromatic amino acids at the protein surface triggered long-distance (∼27 Å) backbone changes due to cooperative effects. These cooperative interactions produced an unprecedented double-lock system that: (i) induced backbone changes, (ii) closed the molecular gates to the active site, (iii) reduced the volumes of the main and slot access tunnels, and (iv) occluded the active site. Despite these spatial restrictions, experimental tracing of the access tunnels using krypton derivative crystals demonstrates that transport of ligands is still effective. Our findings highlight key thermostabilization effects and provide a structural basis for designing new thermostable protein catalysts.

Reengineering biocatalysts: Computational redesign of chondroitinase ABC improves efficacy and stability

Maintaining biocatalyst stability and activity is a critical challenge. Chondroitinase ABC (ChABC) has shown promise in central nervous system (CNS) regeneration, yet its therapeutic utility is severely limited by instability. We computationally reengineered ChABC by introducing 37, 55, and 92 amino acid changes using consensus design and forcefield-based optimization. All mutants were more stable than wild-type ChABC with increased aggregation temperatures between 4° and 8°C. Only ChABC with 37 mutations (ChABC-37) was more active and had a 6.5 times greater half-life than wild-type ChABC, increasing to 106 hours (4.4 days) from only 16.8 hours. ChABC-37, expressed as a fusion protein with Src homology 3 (ChABC-37-SH3), was active for 7 days when released from a hydrogel modified with SH3-binding peptides. This study demonstrates the broad opportunity to improve biocatalysts through computational engineering and sets the stage for future testing of this substantially improved protein in the treatment of debilitating CNS injuries.

One Pot Use of Combilipases for Full Modification of Oils and Fats: Multifunctional and Heterogeneous Substrates

Lipases are among the most utilized enzymes in biocatalysis. In many instances, the main reason for their use is their high specificity or selectivity. However, when full modification of a multifunctional and heterogeneous substrate is pursued, enzyme selectivity and specificity become a problem. This is the case of hydrolysis of oils and fats to produce free fatty acids or their alcoholysis to produce biodiesel, which can be considered cascade reactions. In these cases, to the original heterogeneity of the substrate, the presence of intermediate products, such as diglycerides or monoglycerides, can be an additional drawback. Using these heterogeneous substrates, enzyme specificity can promote that some substrates (initial substrates or intermediate products) may not be recognized as such (in the worst case scenario they may be acting as inhibitors) by the enzyme, causing yields and reaction rates to drop. To solve this situation, a mixture of lipases with different specificity, selectivity and differently affected by the reaction conditions can offer much better results than the use of a single lipase exhibiting a very high initial activity or even the best global reaction course. This mixture of lipases from different sources has been called “combilipases” and is becoming increasingly popular. They include the use of liquid lipase formulations or immobilized lipases. In some instances, the lipases have been coimmobilized. Some discussion is offered regarding the problems that this coimmobilization may give rise to, and some strategies to solve some of these problems are proposed. The use of combilipases in the future may be extended to other processes and enzymes.

Evaluating Protein Engineering Thermostability Prediction Tools Using an Independently Generated Dataset

Engineering proteins to enhance thermal stability is a widely utilized approach for creating industrially relevant biocatalysts. The development of new experimental datasets and computational tools to guide these engineering efforts remains an active area of research. Thus, to complement the previously reported measures of T 50 and kinetic constants, we are reporting an expansion of our previously published dataset of mutants for β-glucosidase to include both measures of T M and ΔΔG. For a set of 51 mutants, we found that T 50 and T M are moderately correlated, with a Pearson correlation coefficient and Spearman's rank coefficient of 0.58 and 0.47, respectively, indicating that the two methods capture different physical features. The performance of predicted stability using nine computational tools was also evaluated on the dataset of 51 mutants, none of which are found to be strong predictors of the observed changes in T 50, T M, or ΔΔG. Furthermore, the ability of the nine algorithms to predict the production of isolatable soluble protein was examined, which revealed that Rosetta ΔΔG, FoldX, DeepDDG, PoPMuSiC, and SDM were capable of predicting if a mutant could be produced and isolated as a soluble protein. These results further highlight the need for new algorithms for predicting modest, yet important, changes in thermal stability as well as a new utility for current algorithms for prescreening designs for the production of mutants that maintain fold and soluble production properties.

Robust ω-Transaminases by Computational Stabilization of the Subunit Interface

Transaminases are attractive catalysts for the production of enantiopure amines. However, the poor stability of these enzymes often limits their application in biocatalysis. Here, we used a framework for enzyme stability engineering by computational library design (FRESCO) to stabilize the homodimeric PLP fold type I ω-transaminase from Pseudomonas jessenii. A large number of surface-located point mutations and mutations predicted to stabilize the subunit interface were examined. Experimental screening revealed that 10 surface mutations out of 172 tested were indeed stabilizing (6% success), whereas testing 34 interface mutations gave 19 hits (56% success). Both the extent of stabilization and the spatial distribution of stabilizing mutations showed that the subunit interface was critical for stability. After mutations were combined, 2 very stable variants with 4 and 6 mutations were obtained, which in comparison to wild type (T _m ^app = 62 °C) displayed T _m ^app values of 80 and 85 °C, respectively. These two variants were also 5-fold more active at their optimum temperatures and tolerated high concentrations of isopropylamine and cosolvents. This allowed conversion of 100 mM acetophenone to (S)-1-phenylethylamine (>99% enantiomeric excess) with high yield (92%, in comparison to 24% with the wild-type transaminase). Crystal structures mostly confirmed the expected structural changes and revealed that the most stabilizing mutation, I154V, featured a rarely described stabilization mechanism: namely, removal of steric strain. The results show that computational interface redesign can be a rapid and powerful strategy for transaminase stabilization.