Computer-Assisted Recombination (CompassR) Teaches us How to Recombine Beneficial Substitutions from Directed Evolution Campaigns

Engineering of Methionine Adenosyltransferase toward Mitigated Product Inhibition for Efficient Production of S-Adenosylmethionine

S-Adenosylmethionine (SAM) is a crucial metabolic intermediate playing irreplaceable roles in organismal activities. However, the synthesis of SAM by methionine adenosyltransferase (MAT) is hindered by low conversion due to severe product inhibition. Herein structure-guided semirational engineering was conducted on MAT from Escherichia coli (EcMAT) to mitigate the product inhibitory effect. Compared with the wild-type EcMAT, the best variant E56Q/Q105R exhibited an 8.13-fold increase in half maximal inhibitory concentration and a 4.46-fold increase in conversion (150 mM ATP and l-methionine), leading to a SAM titer of 47.02 g/L. Another variant, E56N/Q105R, showed superior thermostability with an impressive 85.30-fold increase in half-life (50 °C) value. Furthermore, molecular dynamics (MD) simulation results demonstrate that the alleviation in product inhibitory effect could be attributed to facilitated product release. This study offers molecular insights into the mitigated product inhibition, and provides valuable guidance for engineering MAT toward enhanced catalytic performance.

Multi-strategy orthogonal enhancement and analysis of aldo-keto reductase thermal stability

Given their outstanding efficiency and selectivity, enzymes are integral in various domains such as drug synthesis, the food industry, and environmental management. However, the inherent instability of natural enzymes limits their widespread industrial application. In this study, we underscore the efficacy of enhancing protein thermal stability through comprehensive protein design strategies, encompassing elements such as the free energy of protein folding, internal forces within proteins, and the overall structural design. We also demonstrate the efficiency and precision of combinatorial screening in the thermal stability design of aldo-keto reductase (AKR7-2-1). In our research, three single-point mutations and five combinatorial mutations were strategically introduced into AKR7-2-1, using multiple computational techniques. Notably, the E12I/S235I mutant showed a significant increase of 25.4 °C in its melting temperature (Tm). Furthermore, the optimal mutant, E12V/S235I, maintained 80 % of its activity while realizing a 16.8 °C elevation in Tm. Remarkably, its half-life at 50 °C was increased to twenty times that of the wild type. Structural analysis indicates that this enhanced thermal stability primarily arises from reduced oscillation in the loop region and increased internal hydrogen bonding. The promising results achieved with AKR7-2-1 demonstrate that our strategy could serve as a valuable reference for enhancing the thermal stability of other industrial enzymes.

Empowering Protein Engineering through Recombination of Beneficial Substitutions

Directed evolution stands as a seminal technology for generating novel protein functionalities, a cornerstone in biocatalysis, metabolic engineering, and synthetic biology. Today, with the development of various mutagenesis methods and advanced analytical machines, the challenge of diversity generation and high-throughput screening platforms is largely solved, and one of the remaining challenges is: how to empower the potential of single beneficial substitutions with recombination to achieve the epistatic effect. This review overviews experimental and computer-assisted recombination methods in protein engineering campaigns. In addition, integrated and machine learning-guided strategies were highlighted to discuss how these recombination approaches contribute to generating the screening library with better diversity, coverage, and size. A decision tree was finally summarized to guide the further selection of proper recombination strategies in practice, which was beneficial for accelerating protein engineering.

Corner Engineering: Tailoring Enzymes for Enhanced Resistance and Thermostability in Deep Eutectic Solvents

Deep eutectic solvents (DESs), heralded for their synthesis simplicity, economic viability, and reduced volatility and flammability, have found increasing application in biocatalysis. However, challenges persist due to a frequent diminution in enzyme activity and stability. Herein, we developed a general protein engineering strategy, termed corner engineering, to acquire DES-resistant and thermostable enzymes via precise tailoring of the transition region in enzyme structure. Employing Bacillus subtilis lipase A (BSLA) as a model, we delineated the engineering process, yielding five multi-DESs resistant variants with highly improved thermostability, such as K88E/N89 K exhibited up to a 10.0-fold catalytic efficiency (kcat /KM ) increase in 30 % (v/v) choline chloride (ChCl): acetamide and 4.1-fold in 95 % (v/v) ChCl: ethylene glycol accompanying 6.7-fold thermal resistance improvement than wild type at ≈50 °C. The generality of the optimized approach was validated by two extra industrial enzymes, endo-β-1,4-glucanase PvCel5A (used for biofuel production) and esterase Bs2Est (used for plastics degradation). The molecular investigations revealed that increased water molecules at substrate binding cleft and finetuned helix formation at the corner region are two dominant determinants governing elevated resistance and thermostability. This study, coupling corner engineering with obtained molecular insights, illuminates enzyme-DES interaction patterns and fosters the rational design of more DES-resistant and thermostable enzymes in biocatalysis and biotransformation.

Directed Evolution of Material Binding Peptide for Polylactic Acid-specific Degradation in Mixed Plastic Wastes

In order to preserve our livelihood for future generations, responsible use of plastics in a climate-neutral and circular economy has to be developed so that plastics can be used in an environmentally friendly way by future generations. The prerequisite is that bioplastic polymers such as polylactic acid (PLA) can be efficiently recycled from petrochemical based plastic. Here, a concept in which accelerated PLA degradation in the mixed suspension of PLA and polystyrene (PS) nanoparticles has been achieved through an engineered material binding peptide. After comparison of twenty material binding peptides, Cg-Def is selected due to its PLA binding specificity. Finally, a suitable high-throughput screening system is developed for enhancing material-specific binding toward PLA in presence of PS. Through KnowVolution campaign, a variant Cg-Def YH (L9Y/S19H) with 2.0-fold improved PLA binding specificity compared to PS is generated. Contact angle and surface plasmon resonance measurements validated higher surface coverage of Cg-Def YH on PLA surface and the fusion of Cg-Def YH with PLA degrading enzyme confirmed the accelerated PLA depolymerization (two times higher than only enzyme) in mixed PLA/PS plastics.

A Competitive High-Throughput Screening Platform for Designing Polylactic Acid-Specific Binding Peptides

Among biobased polymers, polylactic acid (PLA) is recognized as one of the most promising bioplastics to replace petrochemical-based polymers. PLA is typically blended with other polymers such as polypropylene (PP) for improved melt processability, thermal stability, and stiffness. A technical challenge in recycling of PLA/PP blends is the sorting/separation of PLA from PP. Material binding peptides (MBPs) can bind to various materials. Engineered MBPs that can bind in a material-specific manner have a high potential for material-specific detection or enhanced degradation of PLA in mixed PLA/PP plastics. To obtain a material-specific MBP for PLA binding (termed PLAbodies ), protein engineering of MBP Cg-Def for improved PLA binding specificity is reported in this work. In detail, a 96-well microtiter plate based high-throughput screening system for PLA specific binding (PLABS) was developed and validated in a protein engineering (KnowVolution) campaign. Finally, the Cg-Def variant V2 (Cg-Def S19K/K10L/N13H) with a 2.3-fold improved PLA binding specificity compared to PP was obtained. Contact angle and surface plasmon resonance measurements confirmed improved material-specific binding of V2 to PLA (1.30-fold improved PLA surface coverage). The established PLABS screening platform represents a general methodology for designing PLAbodies for applications in detection, sorting, and material-specific degradation of PLA in mixed plastics.

A step forward to the optimized HlyA type 1 secretion system through directed evolution

Secretion of proteins into the extracellular space has great advantages for the production of recombinant proteins. Type 1 secretion systems (T1SS) are attractive candidates to be optimized for biotechnological applications, as they have a relatively simple architecture compared to other classes of secretion systems. A paradigm of T1SS is the hemolysin A type 1 secretion system (HlyA T1SS) from Escherichia coli harboring only three membrane proteins, which makes the plasmid-based expression of the system easy. Although for decades the HlyA T1SS has been successfully applied for secretion of a long list of heterologous proteins from different origins as well as peptides, but its utility at commercial scales is still limited mainly due to low secretion titers of the system. To address this drawback, we engineered the inner membrane complex of the system, consisting of HlyB and HlyD proteins, following KnowVolution strategy. The applied KnowVolution campaign in this study provided a novel HlyB variant containing four substitutions (T36L/F216W/S290C/V421I) with up to 2.5-fold improved secretion for two hydrolases, a lipase and a cutinase. KEY POINTS: • An improvement in protein secretion via the use of T1SS • Reaching almost 400 mg/L of soluble lipase into the supernatant • A step forward to making E. coli cells more competitive for applying as a secretion host.

Improving thermostability of a PL 5 family alginate lyase with combination of rational design strategies

Alginate lyases with strict substrate specificity possess potential in directed production of alginate oligosaccharides with specific composition. However, their poor thermostability hampered their applications in industry. In this study, an efficient comprehensive strategy including sequence-based analysis, structure-based analysis, and computer-aid ΔΔGfold value calculation was proposed. It was successfully performed on alginate lyase (PMD) with strict poly-β-D-mannuronic acid substrate specificity. Four single-point variants A74V, G75V, A240V, and D250G with increased T_m of 3.94 °C, 5.21 °C, 2.56 °C, and 4.80 °C, respectively, were selected out. After ordered combined mutations, a four-point mutant (M4) was finally generated which displayed remarkable increase on thermostability. The T_m of M4 increased from 42.25 °C to 51.59 °C and its half-life at 50 °C was about 58.9-fold of PMD. Meanwhile, there was no obvious loss of enzyme activity (more than 90% retained). Molecular dynamics simulation analysis insisted that the improvement of thermostability might be attribute to the rigidified region A which might be caused by the newly formed hydrogen bonds and salt bridges introduced by mutations, the lower distance of original hydrogen bonds, and the more compact overall structures.

Evolving Robust and Interpretable Enzymes for the Bioethanol Industry

Obtaining a robust and applicable enzyme for bioethanol production is a dream for biorefinery engineers. Herein, we describe a general method to evolve an all-round and interpretable enzyme that can be directly employed in the bioethanol industry. By integrating the transferable protein evolution strategy InSiReP 2.0 (In Silico guided Recombination Process), enzymatic characterization for actual production, and computational molecular understanding, the model cellulase PvCel5A (endoglucanase II Cel5A from Penicillium verruculosum) was successfully evolved to overcome the remaining challenges of low ethanol and temperature tolerance, which primarily limited biomass transformation and bioethanol yield. Remarkably, application of the PvCel5A variants in both first- and second-generation bioethanol production processes (i. Conventional corn ethanol fermentation combined with the in situ pretreatment process; ii. cellulosic ethanol fermentation process) resulted in a 5.7-10.1 % increase in the ethanol yield, which was unlikely to be achieved by other optimization techniques.

Problem of Domain/Building Block Preservation in the Evolution of Biological Macromolecules and Evolutionary Computation

Structurally and functionally isolated domains in biological macromolecular evolution, both natural and artificial, are largely similar to "schemata", building blocks (BBs), in evolutionary computation (EC). The problem of preserving in subsequent evolutionary searches the already found domains / BBs is well known and quite relevant in biology as well as in EC. Both biology and EC are seeing parallel and independent development of several approaches to identifying and preserving previously identified domains / BBs. First, we notice the similarity of DNA shuffling methods in synthetic biology and multi-parent recombination algorithms in EC. Furthermore, approaches to computer identification of domains in proteins that are being developed in biology can be aligned with BB identification methods in EC. Finally, approaches to chimeric protein libraries optimization in biology can be compared to evolutionary search methods based on probabilistic models in EC. We propose to validate the prospects of mutual exchange of ideas and transfer of algorithms and approaches between evolutionary systems biology and EC in these three principal directions. A crucial aim of this transfer is the design of new advanced experimental techniques capable of solving more complex problems of in vitro evolution.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Protein-Glutaminase Engineering Based on Isothermal Compressibility Perturbation for Enhanced Modification of Soy Protein Isolate

Protein-glutaminase plays a significant role in future food (e.g., plant-based meat) processing as a result of its ability to improve the solubility, foaming, emulsifying, and gel properties of plant-based proteins. However, poor stability, activity, high pressure, and high shear processing environments hinder its application. Therefore, we developed an application-oriented method isothermal compressibility perturbation engineering strategy to improve enzyme performance by simulating the high-pressure environment. The best variant with remarkable improvement in specific activity and half-time, N16M/Q21H/T113E, exhibited a 4.28-fold increase compared to the wild type in specific activity (117.18 units/mg) and a 1.23-fold increase in half-time (472 min), as one of the highest comprehensive performances ever reported. The solubility of the soy protein isolate deaminated by the N16M/Q21H/T113E mutant was 55.74% higher than that deaminated by the wild type, with a tinier particle size and coarser texture. Overall, this strategy has the potential to improve the functional performance of enzymes under complex food processing conditions.

Identification and Mutation Analysis of Nonconserved Residues on the TIM-Barrel Surface of GH5_5 Cellulases for Catalytic Efficiency and Stability Improvement

Exploring the potential functions of nonconserved residues on the outer side of α-helices and systematically optimizing them are pivotal for their application in protein engineering. Based on the evolutionary structural conservation analysis of GH5_5 cellulases, a practical molecular improvement strategy was developed. Highly variable sites on the outer side of the α-helices of the GH5_5 cellulase from Aspergillus niger (AnCel5A) were screened, and 14 out of the 34 highly variable sites were confirmed to exert a positive effect on the activity. After the modular combination of the positive mutations, the catalytic efficiency of the mutants was further improved. By using CMC-Na as the substrate, the catalytic efficiency and specific activity of variant AnCel5A_N193A/T300P/D307P were approximately 2.0-fold that of AnCel5A (227 ± 21 versus 451 ± 43 ml/s/mg and 1,726 ± 19 versus 3,472 ± 42 U/mg, respectively). The half-life (t_1/2) of variant AnCel5A_N193A/T300P/D307P at 75°C was 2.36 times that of AnCel5A. The role of these sites was successfully validated in other GH5_5 cellulases. Computational analyses revealed that the flexibility of the loop 6-loop 7-loop 8 region was responsible for the increased catalytic performance. This work not only illustrated the important role of rapidly evolving positions on the outer side of the α-helices of GH5_5 cellulases but also revealed new insights into engineering the proteins that nature left as clues for us to find. IMPORTANCE A comprehensive understanding of the residues on the α-helices of the GH5_5 cellulases is important for catalytic efficiency and stability improvement. The main objective of this study was to use the evolutionary conservation and plasticity of the TIM-barrel fold to probe the relationship between nonconserved residues on the outer side of the α-helices and the catalytic efficiency of GH5_5 cellulases by conducting structure-guided protein engineering. By using a four-step nonconserved residue screening strategy, the functional role of nonconserved residues on the outer side of the α-helices was effectively identified, and a variant with superior performance and capability was constructed. Hence, this study proved the effectiveness of this strategy in engineering GH5_5 cellulases and provided a potential competitor for industrial applications. Furthermore, this study sheds new light on engineering TIM-barrel proteins.

A Plurizyme with Transaminase and Hydrolase Activity Catalyzes Cascade Reactions

Engineering dual‐function single polypeptide catalysts with two abiotic or biotic catalytic entities (or combinations of both) supporting cascade reactions is becoming an important area of enzyme engineering and catalysis. Herein we present the development of a PluriZyme, TR₂E₂, with efficient native transaminase (k_cat: 69.49±1.77 min⁻¹) and artificial esterase (k_cat: 3908–0.41 min⁻¹) activities integrated into a single scaffold, and evaluate its utility in a cascade reaction. TR₂E₂ (pH_opt: 8.0–9.5; T_opt: 60–65 °C) efficiently converts methyl 3‐oxo‐4‐(2,4,5‐trifluorophenyl)butanoate into 3‐(R)‐amino‐4‐(2,4,5‐trifluorophenyl)butanoic acid, a crucial intermediate for the synthesis of antidiabetic drugs. The reaction proceeds through the conversion of the β‐keto ester into the β‐keto acid at the hydrolytic site and subsequently into the β‐amino acid (e.e. >99 %) at the transaminase site. The catalytic power of the TR₂E₂PluriZyme was proven with a set of β‐keto esters, demonstrating the potential of such designs to address bioinspired cascade reactions.

Sustainable isomaltulose production in Corynebacterium glutamicum by engineering the thermostability of sucrose isomerase coupled with one-step simplified cell immobilization

Sucrose isomerase (SI), catalyzing sucrose to isomaltulose, has been widely used in isomaltulose production, but its poor thermostability is still resisted in sustainable batches production. Here, protein engineering and one-step immobilized cell strategy were simultaneously coupled to maintain steady state for long-term operational stabilities. First, rational design of Pantoea dispersa SI (PdSI) for improving its thermostability by predicting and substituting the unstable amino acid residues was investigated using computational analysis. After screening mutagenesis library, two single mutants (PdSIV280L and PdSIS499F) displayed favorable characteristics on thermostability, and further study found that the double mutant PdSIV280L/S499F could stabilize PdSIWT better. Compared with PdSIWT, PdSIV280L/S499F displayed a 3.2°C-higher T m , and showed a ninefold prolonged half-life at 45°C. Subsequently, a one-step simplified immobilization method was developed for encapsulation of PdSIV280L/S499F in food-grade Corynebacterium glutamicum cells to further enhance the recyclability of isomaltulose production. Recombinant cells expressing combinatorial mutant (RCSI2) were successfully immobilized in 2.5% sodium alginate without prior permeabilization. The immobilized RCSI2 showed that the maximum yield of isomaltulose by batch conversion reached to 453.0 g/L isomaltulose with a productivity of 41.2 g/l/h from 500.0 g/L sucrose solution, and the conversion rate remained 83.2% after 26 repeated batches.

Modulating the Coupling Efficiency of P450 BM3 by Controlling Water Diffusion through Access Tunnel Engineering

Cytochromes P450 have gained much interest for their broad substrate scope in the catalysis of oxidation reactions for pharmaceuticals, plastics, and hormones. However, achieving high coupling efficiency by the engineering of P450s is still a big challenge. The presence of extra water around the active site is deemed to be related to uncoupling. In this study, the access tunnels of P450 BM3 from Bacillus megaterium are engineered to control water access from bulk solvent to the active site. Nine residues located in tunnels are investigated by site-saturation mutagenesis to reduce water diffusion, thereby improving the coupling efficiency. The recombined variant N319L/T411V/T436A shows improved coupling efficiency (from 31.2 % to 52.6 %). Tunnel polarity analysis and molecular dynamics simulation further indicate that reduced water molecules around the active site lead to higher coupling efficiency. Overall, this study provides valuable insight on improving coupling efficiency by controlling water diffusion through tunnel engineering.

Polar Substitutions on the Surface of a Lipase Substantially Improve Tolerance in Organic Solvents

Biocatalysis in organic solvents (OSs) enables more efficient routes to the synthesis of various valuable chemicals. However, OSs often reduce enzymatic activity, which limits the use of enzymes in OSs. Herein, we report a comprehensive understanding of interactions between surface polar substitutions and DMSO by integrating molecular dynamics (MD) simulations of 45 variants from Bacillus subtilis lipase A (BSLA) and substitution landscape into a "BSLA-SSM" library. By systematically analyzing 39 structural-, solvation-, and interaction energy-based observables, we discovered that hydration shell maintenance, DMSO reduction, and decreased local flexibility simultaneously govern the stability of polar variants in OS. Moreover, the fingerprints of 1631 polar-related variants in three OSs demonstrated that substituting aromatic to polar amino acid(s) hold great potential to highly improve OSs resistance. Hence, surface polar engineering is a powerful strategy to generate OS-tolerant lipases and other enzymes, thereby adapting the catalyst to the desired reaction and process with OSs.

Computer-Aided Targeted Mutagenesis of Thermoclostridium caenicola d-Allulose 3-Epimerase for Improved Thermostability

d-Allulose 3-epimerase (DAEase) is a key enzyme in d-allulose bioproduction. DAEase from Thermoclostridium caenicola suffers from poor thermostability, hampering its large-scale applications in industry. In this study, mutants A70P, G107P, F155Y, and D162T with increased melting point temperature (T_m) were obtained by targeted mutagenesis based on the calculation of the free energy of folding. The optimal single-point mutant G107P showed 11.08 h, 5, and 5.70 °C increases in the values of half-life (t_1/2) at 60 °C, the optimum temperature (T_opt), and T_m, respectively. Beneficial mutations were combined by ordered recombination mutagenesis, and the combinational mutant Var3 (G107P/F155Y/D162T/A70P) was generated with ΔT_opt of 10 °C and ΔT_m of 12.25 °C. Its t_1/2 value at 65 °C was more than 140 times higher than that of the wild-type enzyme. Molecular dynamics simulations and homology modeling analysis indicated that the enhanced overall rigidity, increased hydrogen bonds between subunits, and redistributed surface electrostatic charges might be responsible for the improved thermostability of the mutant Var3.

Using Molecular Simulation to Guide Protein Engineering for Biocatalysis in Organic Solvents

Biocatalysis in organic solvents (OSs) is very appealing for the industry in producing bulk and/or fine chemicals, such as pharmaceuticals, biodiesel, and fragrances. The poor performance of enzymes in OSs (e.g., reduced activity, insufficient stability, and deactivation) negates OSs' excellent solvent properties. Molecular dynamics (MD) simulations provide a complementary method to study the relationship between enzymes dynamics and the stability in OSs. Here we describe computational procedure for MD simulation of enzymes in OSs with an example of Bacillus subtilis lipase A (BSLA) in dimethyl sulfoxide (DMSO) cosolvent with software GROMACS. We discuss main essential practical issues considered (such as choice of force field, parameterization, simulation setup, and trajectory analysis). The core part of this protocol (enzyme-OS system setup, analysis of structural-based and solvation-based observables) is transferable to other enzymes and any OS systems. Combining with experimental studies, the obtained molecular knowledge is most likely to guide researchers to access rational protein engineering approaches to tailor OS resistant enzymes and expand the scope of biocatalysis in OS media. Finally, we discuss potential solutions to overcome the remaining challenges of computational biocatalysis in OSs and briefly draw future directions for further improvement in this field.

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.

Recombination of Compatible Substitutions by 2GenReP and InSiReP

The CompassR rule enables to identify the beneficial substitutions, which can be recombined in directed evolution with gradually improving the enzymatic properties. However, the question of how to efficiently explore the protein sequence space when ten or more beneficial substitutions are identified has not yet been addressed. Two recombination strategies 2GenReP and InSiReP employing CompassR are systematically investigated to minimize experimental efforts and maximize possible improvements. Here we describe the details of the 2GenReP and InSiReP procedure with an example of recombining 15 substitutions and discuss some important practical issues that should be considered for the application of 2GenReP and InSiReP, such as placing the substitutions into subsets. The core part of the protocol (Step1 to Step5) is transferable to other enzymes and any recombination of potential substitutions.

Recombination of Single Beneficial Substitutions Obtained from Protein Engineering by Computer-Assisted Recombination (CompassR)

A large number of beneficial substitutions can be obtained from a successful directed enzyme evolution campaign and/or (semi)rational design. It is expected that the recombination of some beneficial substitutions leads to a much higher degree of performance through synergistic effect. However, systematic recombination studies show that poorly performing variants are often obtained after recombination of three to four individual beneficial substitutions and this limits protein engineers to exploit nature's potential in generating better performing enzymes. Computer-assisted Recombination (CompassR) strategy allows the recombination of identified beneficial substitutions in an effective and efficient manner in order to generate active enzymes with improved performance. Here, we describe in detail the CompassR procedure with an example of recombining four substitutions and discuss some important practical issues that should be considered (such as the selection of protein structures, number of FoldX runs, evaluation of calculations) for application of the CompassR rule. The core part of this protocol (system setup, ΔΔG_fold calculation, and CompassR application) is transferable to other enzymes and any recombination of single beneficial substitutions.

PyPEF-An Integrated Framework for Data-Driven Protein Engineering

Data-driven strategies are gaining increased attention in protein engineering due to recent advances in access to large experimental databanks of proteins, next-generation sequencing (NGS), high-throughput screening (HTS) methods, and the development of artificial intelligence algorithms. However, the reliable prediction of beneficial amino acid substitutions, their combination, and the effect on functional properties remain the most significant challenges in protein engineering, which is applied to develop proteins and enzymes for biocatalysis, biomedicine, and life sciences. Here, we present a general-purpose framework (PyPEF: pythonic protein engineering framework) for performing data-driven protein engineering using machine learning methods combined with techniques from signal processing and statistical physics. PyPEF guides the identification and selection of beneficial proteins of a defined sequence space by systematically or randomly exploring the fitness of variants and by sampling random evolution pathways. The performance of PyPEF was evaluated concerning its predictive accuracy and throughput on four public protein and enzyme data sets using common regression models. It was proved that the program could efficiently predict the fitness of protein sequences for different target properties (predictive models with coefficient of determination values ranging from 0.58 to 0.92). By combining machine learning and protein evolution, PyPEF enabled the screening of proteins with various functions, reaching a screening capacity of more than 500,000 protein sequence variants in the timeframe of only a few minutes on a personal computer. PyPEF displayed significant accuracies on four public data sets (different proteins and properties) and underlined the potential of integrating data-driven technologies for covering different philosophies by either predicting the fitness of the variants to the highest accuracy accounting for epistatic effects or capturing the general trend of introduced mutations on the fitness in directed protein evolution campaigns. In essence, PyPEF can provide a powerful solution to current sequence exploration and combinatorial problems faced in protein engineering through exhaustive in silico screening of the sequence space.

Aqueous ionic liquids redistribute local enzyme stability via long-range perturbation pathways

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for biocatalysis due to their unique properties. On the other hand, the incubation of enzymes in IL or aIL often reduces enzyme activity. Recent studies proposed various aIL-induced effects to explain the reduction, classified as direct effects, e.g., local dehydration or competitive inhibition, and indirect effects, e.g., structural perturbations or disturbed catalytic site integrity. However, the molecular origin of indirect effects has largely remained elusive. Here we show by multi-μs long molecular dynamics simulations, free energy computations, and rigidity analyses that aIL favorably interact with specific residues of Bacillus subtilis Lipase A (BsLipA) and modify the local structural stability of this model enzyme by inducing long-range perturbations of noncovalent interactions. The perturbations percolate over neighboring residues and eventually affect the catalytic site and the buried protein core. Validation against a complete experimental site saturation mutagenesis library of BsLipA (3620 variants) reveals that the residues of the perturbation pathways are distinguished sequence positions where substitutions highly likely yield significantly improved residual activity. Our results demonstrate that identifying these perturbation pathways and specific IL ion-residue interactions there effectively predicts focused variant libraries with improved aIL tolerance.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

An engineered cellobiohydrolase I for sustainable degradation of lignocellulosic biomass

This study provides computational-assisted engineering of the cellobiohydrolase I (CBH-I) from Penicillium verruculosum with simultaneous enhanced thermostability and tolerance in ionic liquids, deep eutectic solvent, and concentrated seawater without affecting its wild-type activity. Engineered triple variant CBH-I R1 (A65R-G415R-S181F) showed 2.48-fold higher thermostability in terms of relative activity at 65°C after 1 h of incubation when compared with CBH-I wild type. CBH-I R1 exhibited 1.87-fold, 1.36-fold, and 1.57-fold higher specific activities compared with CBH-I wild type in [Bmim]Cl (50 g/L), [Ch]Cl (50 g/L), and two-fold concentrated seawater, respectively. In the multicellulases mixture, CBH-I R1 showed higher hydrolytic efficiency to hydrolyze aspen wood compared with CBH-I wild type in the buffer, [Bmim]Cl (50 g/L), and two-fold concentrated seawater, respectively. Structural analysis revealed a molecular basis for the higher stability of the CBH-I structure in which A65R and G415R substitutions form salt bridges (D64 … R65, E411 … R415) and S181F forms π-π interaction (Y155 … F181), leading to stabilize surface-exposed flexible α-helixes and loop in the multidomain β-jelly roll fold structure, respectively. In conclusion, the variant CBH-I R1 could enable efficient lignocellulosic biomass degradation as a cost-effective alternative for the sustainable production of biofuels and value-added chemicals.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Tunnel engineering for modulating the substrate preference in cytochrome P450BsβHI

An active site is normally located inside enzymes, hence substrates should go through a tunnel to access the active site. Tunnel engineering is a powerful strategy for refining the catalytic properties of enzymes. Here, P450BsβHI (Q85H/V170I) derived from hydroxylase P450Bsβ from Bacillus subtilis was chosen as the study model, which is reported as a potential decarboxylase. However, this enzyme showed low decarboxylase activity towards long-chain fatty acids. Here, a tunnel engineering campaign was performed for modulating the substrate preference and improving the decarboxylation activity of P450BsβHI. The finally obtained BsβHI-F79A variant had a 15.2-fold improved conversion for palmitic acid; BsβHI-F173V variant had a 3.9-fold improved conversion for pentadecanoic acid. The study demonstrates how the substrate preference can be modulated by tunnel engineering strategy.

A Computational Method to Predict Effects of Residue Mutations on the Catalytic Efficiency of Hydrolases

With scientific and technological advances, growing research has focused on engineering enzymes that acquire enhanced efficiency and activity. Thereinto, computer-based enzyme modification makes up for the time-consuming and labor-intensive experimental methods and plays a significant role. In this study, for the first time, we collected and manually curated a data set for hydrolases mutation, including structural information of enzyme-substrate complexes, mutated sites and Kcat/Km obtained from vitro assay. We further constructed a classification model using the random forest algorithm to predict the effects of residue mutations on catalytic efficiency (increase or decrease) of hydrolases. This method has achieved impressive performance on a blind test set with the area under the receiver operating characteristic curve of 0.86 and the Matthews Correlation Coefficient of 0.659. Our results demonstrate that computational mutagenesis has an instructive effect on enzyme modification, which may expedite the design of engineering hydrolases.

CompassR Yields Highly Organic-Solvent-Tolerant Enzymes through Recombination of Compatible Substitutions

The CompassR (computer-assisted recombination) rule enables, among beneficial substitutions, the identification of those that can be recombined in directed evolution. Herein, a recombination strategy is systematically investigated to minimize experimental efforts and maximize possible improvements. In total, 15 beneficial substitutions from Bacillus subtilis lipase A (BSLA), which improves resistance to the organic cosolvent 1,4-dioxane (DOX), were studied to compare two recombination strategies, the two-gene recombination process (2GenReP) and the in silico guided recombination process (InSiReP), employing CompassR. Remarkably, both strategies yielded a highly DOX-resistant variant, M4 (I12R/Y49R/E65H/N98R/K122E/L124K), with up to 14.6-fold improvement after screening of about 270 clones. M4 has a remarkably enhanced resistance in 60 % (v/v) acetone (6.0-fold), 30 % (v/v) ethanol (2.1-fold), and 60 % (v/v) methanol (2.4-fold) compared with wild-type BSLA. Molecular dynamics simulations revealed that attracting water molecules by charged surface substitutions is the main driver for increasing the DOX resistance of BSLA M4. Both strategies and obtained molecular knowledge can likely be used to improve the properties of other enzymes with a similar α/β-hydrolase fold.

KnowVolution of prodigiosin ligase PigC towards condensation of short-chain prodiginines

One round of KnowVolution enhanced the catalytic activity of prodigiosin ligase PigC with short-chain monopyrroles, opening access to anticancer prodiginines.

Understanding substrate binding and the role of gatekeeping residues in PigC access tunnels

Prodigiosin ligase PigC has been engineered by semi-rational design to accept short chain-pyrroles, providing molecular understanding of access tunnels and the substrate-binding pocket.

Can constraint network analysis guide the identification phase of KnowVolution? A case study on improved thermostability of an endo-β-glucanase

Cellulases are industrially important enzymes, e.g., in the production of bioethanol, in pulp and paper industry, feedstock, and textile. Thermostability is often a prerequisite for high process stability and improving thermostability without affecting specific activities at lower temperatures is challenging and often time-consuming. Protein engineering strategies that combine experimental and computational are emerging in order to reduce experimental screening efforts and speed up enzyme engineering campaigns. Constraint Network Analysis (CNA) is a promising computational method that identifies beneficial positions in enzymes to improve thermostability. In this study, we compare CNA and directed evolution in the identification of beneficial positions in order to evaluate the potential of CNA in protein engineering campaigns (e.g., in the identification phase of KnowVolution). We engineered the industrially relevant endoglucanase EGLII from Penicillium verruculosum towards increased thermostability. From the CNA approach, six variants were obtained with an up to 2-fold improvement in thermostability. The overall experimental burden was reduced to 40% utilizing the CNA method in comparison to directed evolution. On a variant level, the success rate was similar for both strategies, with 0.27% and 0.18% improved variants in the epPCR and CNA-guided library, respectively. In essence, CNA is an effective method for identification of positions that improve thermostability.

MicroGelzymes: pH-Independent Immobilization of Cytochrome P450 BM3 in Microgels

Microgels are an emerging class of "ideal" enzyme carriers because of their chemical and process stability, biocompatibility, and high enzyme loading capability. In this work, we synthesized a new type of permanently positively charged poly(N-vinylcaprolactam) (PVCL) microgel with 1-vinyl-3-methylimidazolium (quaternization of nitrogen by methylation of N-vinylimidazole moieties) as a comonomer (PVCL/VimQ) through precipitation polymerization. The PVCL/VimQ microgels were characterized with respect to their size, charge, swelling degree, and temperature responsiveness in aqueous solutions. P450 monooxygenases are usually challenging to immobilize, and often, high activity losses occur after the immobilization (in the case of P450 BM3 from Bacillus megaterium up to 100% loss of activity). The electrostatic immobilization of P450 BM3 in permanently positively charged PVCL/VimQ microgels was achieved without the loss of catalytic activity at the pH optimum of P450 BM3 (pH 8; ∼9.4 nmol 7-hydroxy-3-carboxy coumarin ethyl ester/min for free and immobilized P450 BM3); the resulting P450-microgel systems were termed P450 MicroGelzymes (P450 μ-Gelzymes). In addition, P450 μ-Gelzymes offer the possibility of reversible ionic strength-triggered release and re-entrapment of the biocatalyst in processes (e.g., for catalyst reuse). Finally, a characterization of the potential of P450 μ-Gelzymes to provide resistance against cosolvents (acetonitrile, dimethyl sulfoxide, and 2-propanol) was performed to evaluate the biocatalytic application potential.

Computer-Assisted Recombination (CompassR) Teaches us How to Recombine Beneficial Substitutions from Directed Evolution Campaigns

A main remaining challenge in protein engineering is how to recombine beneficial substitutions. Systematic recombination studies show that poorly performing variants are usually obtained after recombination of 3 to 4 beneficial substitutions. This limits researchers in exploiting nature's potential in generating better enzymes. The Computer-assisted Recombination (CompassR) strategy provides a selection guide for beneficial substitutions that can be recombined to gradually improve enzyme performance by analysis of the relative free energy of folding (ΔΔGfold ). The performance of CompassR was evaluated by analysis of 84 recombinants located on 13 positions of Bacillus subtilis lipase A. The finally obtained variant F17S/V54K/D64N/D91E had a 2.7-fold improved specific activity in 18.3 % (v/v) 1-butyl-3-methylimidazolium chloride ([BMIM][Cl]). In essence, the deducted CompassR rule allows recombination of beneficial substitutions in an iterative manner and empowers researchers to generate better enzymes in a time-efficient manner.