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

Classification and functional origins of stereocomplementary alcohol dehydrogenases for asymmetric synthesis of chiral secondary alcohols: A review

Alcohol dehydrogenases (ADHs) mediated biocatalytic asymmetric reduction of ketones have been widely applied in the synthesis of optically active secondary alcohols with highly reactive hydroxyl groups ligated to the stereogenic carbon and divided into (R)- and (S)-configurations. Stereocomplementary ADHs could be applied in the synthesis of both enantiomers and are increasingly accepted as the "first of choice" in green chemistry due to the high atomic economy, low environmental factor, 100 % theoretical yield, and high environmentally friendliness. Due to the equal importance of complementary alcohols, development of stereocomplementary ADHs draws increasing attention. This review is committed to summarize recent advance in discovery of naturally evolved and tailor-made stereocomplementary ADHs, unveil the molecular mechanism of stereoselective catalysis in views of classification and functional basis, and provide guidance for further engineering the stereoselectivity of ADHs for the industrial biosynthesis of chiral secondary alcohol of industrial relevance.

Engineering diaryl alcohol dehydrogenase KpADH reveals importance of retaining hydration shell in organic solvent tolerance

Alcohol dehydrogenases (ADHs) are synthetically important biocatalysts for the asymmetric synthesis of chiral alcohols. The catalytic performance of ADHs in the presence of organic solvents is often important since most prochiral ketones are highly hydrophobic. Here, the organic solvent tolerance of KpADH from Kluyveromyces polyspora was semi-rationally evolved. Using tolerant variants obtained, meticulous experiments and computational studies were conducted to explore properties including stability, activity and kinetics in the presence of various organic solvents. Compared with WT, variant V231D exhibited 1.9-fold improvement in ethanol tolerance, while S237G showed a 6-fold increase in catalytic efficiency, a higherT 50 15 $$ {\mathrm{T}}_{50}^{15} $$ , as well as 15% higher tolerance in 7.5% (v/v) ethanol. Based on 3 × 100 ns MD simulations, the increased tolerance of V231D and S237G against ethanol may be ascribed to their enhanced ability in retaining water molecules and repelling ethanol molecules. Moreover, 6.3-fold decreased KM value of V231D toward hydrophilic ketone substrate confirmed its capability of retaining hydration shell. Our results suggest that retaining hydration shell surrounding KpADH is critical for its tolerance to organic solvents, as well as catalytic performance. This study provides useful guidance for engineering organic solvent tolerance of KpADH and other ADHs.

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.

Mechanistic Insights into the Halophilic Xylosidase Xylo-1 and Its Role in Xylose Production

The Xylo-1 xylosidase, which belongs to the GH43 family, exhibits a high salt tolerance. The present study demonstrated that the catalytic activity of Xylo-1 increased by 195% in the presence of 5 M NaCl. Additionally, the half-life of Xylo-1 increased 25.9-fold in the presence of 1 M NaCl. Through comprehensive analysis including circular dichroism, fluorescence spectroscopy, and molecular dynamics simulations, we elucidated that the presence of Na+ ions increased the contact frequency between the surface acidic amino acids and the surrounding water molecules. This resulted in the stabilization of the surrounding hydration layer of Xylo-1. Additionally, Na+ ions also stabilized the substrate-binding conformation and the fluctuation of water molecules within the active site, which enhanced the catalytic activity of Xylo-1 by increasing the nucleophilic attack by the water molecules. Ultimately, the optimal reaction conditions for the production of xylose by synergistic catalysis with Xylo-1 and xylanase were determined. The results demonstrated that the conversion yield of the method was high for various sources of xylan, indicating the method could have potential industrial applications. This study explored the structure-activity relationship of catalysis in Xylo-1 under high-salt conditions, provides novel insights into the mechanism of halophilic enzymes, and serves as a reference for the industrial application of Xylo-1.

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.

Engineering All-Round Cellulase for Bioethanol Production

One strategy to decrease both the consumption of crude oil and environmental damage is through the production of bioethanol from biomass. Cellulolytic enzyme stability and enzymatic hydrolysis play important roles in the bioethanol process. However, the gradually increased ethanol concentration often reduces enzyme activity and leads to inactivation, thereby limiting the final ethanol yield. Herein, we employed an optimized Two-Gene Recombination Process (2GenReP) approach to evolve the exemplary cellulase CBHI for practical bioethanol fermentation. Two all-round CBHI variants (named as R2 and R4) were obtained with simultaneously improved ethanol resistance, organic solvent inhibitor tolerance, and enzymolysis stability in simultaneous saccharification and fermentation (SSF). Notably, CBHI R4 had a 7.0- to 34.5-fold enhanced catalytic efficiency (kcat/KM) in the presence/absence of ethanol. Employing the evolved CBHI R2 and R4 in the 1G bioethanol process resulted in up to 10.27% (6.7 g/L) improved ethanol yield (ethanol concentration) than non-cellulase, which was far more beyond than other optimization strategies. Besides bioenergy fields, this transferable protein engineering routine holds the potential to generate all-round enzymes that meet the requirement in biotransformation and bioenergy fields.

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.

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.

Deletion of the Loop Linking Two Domains of Exo-Inulinase InuAMN8 Diminished the Enzymatic Thermo-Halo-Alcohol Tolerance

Inulin is the rich water-soluble storage polysaccharide after starch in nature, and utilization of inulin through hydrolysis of exo-inulinases has attracted much attention. Thermo-halo-alcohol tolerance is essential for exo-inulinase applications, while no report reveals the molecular basis involved in halo-alcohol tolerance of exo-inulinases via experimental data. In this study, two loops of exo-inulinase InuAMN8, including the loop built with ³⁶⁰GHVRLGPQP³⁶⁸ linking domains of Glyco_hydro_32N and Glyco_hydro_32C and another loop built with ¹⁶⁹GGAG¹⁷² in the catalytic domain, were deleted to generate mutants MutG360Δ9 and MutG169Δ4, respectively. After heterologous expression, purification, and dialysis, InuAMN8, MutG169Δ4, and MutG360Δ9 showed half-lives of 144, 151, and 7 min at 50°C, respectively. InuAMN8 and MutG169Δ4 were very stable, while MutG360Δ9 showed a half-life of approximately 60 min in 5.0% (w/v) NaCl, and they showed half-lives of approximately 60 min in 25.0, 25.0, and 5.0% (w/v) ethanol, respectively. Structural analysis indicated that two cation-π bonds, which contributed to thermal properties of InuAMN8 at high temperatures, broke in MutG360Δ9. Four basic amino acid residues were exposed to the structural surface of MutG360Δ9 and formed positive and neutral electrostatic potential that caused detrimental effects on halo-alcohol tolerance. The study may provide a better understanding of the loop-function relationships that are involved in thermo-halo-alcohol adaptation of enzymes in extreme environment.

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.

Immobilization of Eversa Lipases on Hydrophobic Supports for Ethanolysis of Sunflower Oil Solvent-Free

Lipases are an important group of biocatalysts for many industrial applications. Two new commercial low-cost lipases Eversa® Transform and Eversa® Transform 2.0 was immobilized on four different hydrophobic supports: Lewatit-DVB, Purolite-DVB, Sepabeads-C18, and Purolite-C18. The performance of immobilized lipases was investigated in the transesterification of sunflower oil solvent-free in an anhydrous medium. Interesting results were obtained for both lipases and the four supports, but with Sepabeads support the lipases Eversa showed high catalytic activity. However, the more stable and efficient derivative was Eversa® Transform immobilized on Sepabeads C-18. A 98 wt% of ethyl ester of fatty acid (FAEE) was obtained, in 3 h at 40ºC, ethanol/sunflower oil molar ratio of 3:1 and a 10 wt% of the immobilized biocatalyst. After 6 reaction cycles, the immobilized biocatalyst preserved 70 wt% of activity. Both lipases immobilized in Sepabeads C-18 were highly active and stable in the presence of ethanol. The immobilization of Eversa Transform and Eversa Transform 2.0 in hydrophobic supports described in this study appears to be a promising alternative to the immobilization and application of these news lipases still unexplored.

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.

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.

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.