Navigating the landscape of enzyme design: from molecular simulations to machine learning

Global environmental issues and sustainable development call for new technologies for fine chemical synthesis and waste valorization. Biocatalysis has attracted great attention as the alternative to the traditional organic synthesis. However, it is challenging to navigate the vast sequence space to identify those proteins with admirable biocatalytic functions. The recent development of deep-learning based structure prediction methods such as AlphaFold2 reinforced by different computational simulations or multiscale calculations has largely expanded the 3D structure databases and enabled structure-based design. While structure-based approaches shed light on site-specific enzyme engineering, they are not suitable for large-scale screening of potential biocatalysts. Effective utilization of big data using machine learning techniques opens up a new era for accelerated predictions. Here, we review the approaches and applications of structure-based and machine-learning guided enzyme design. We also provide our view on the challenges and perspectives on effectively employing enzyme design approaches integrating traditional molecular simulations and machine learning, and the importance of database construction and algorithm development in attaining predictive ML models to explore the sequence fitness landscape for the design of admirable biocatalysts.

Evolution of the catalytic mechanism at the dawn of the Baeyer-Villiger monooxygenases

Enzymes are crucial for the emergence and sustenance of life on earth. How they became catalytically active during their evolution is still an open question. Two opposite explanations are plausible: acquiring a mechanism in a series of discrete steps or all at once in a single evolutionary event. Here, we use molecular phylogeny, ancestral sequence reconstruction, and biochemical characterization to follow the evolution of a specialized group of flavoprotein monooxygenases, the bacterial Baeyer-Villiger monooxygenases (BVMOs). These enzymes catalyze an intricate chemical reaction relying on three different elements: a reduced nicotinamide cofactor, dioxygen, and a substrate. Characterization of ancestral BVMOs shows that the catalytic mechanism evolved in a series of steps starting from a FAD-binding protein and further acquiring reactivity and specificity toward each of the elements participating in the reaction. Together, the results of our work portray how an intrinsically complex catalytic mechanism emerged during evolution.

LibGENiE - A bioinformatic pipeline for the design of information-enriched enzyme libraries

Enzymes are potent catalysts with high specificity and selectivity. To leverage nature's synthetic potential for industrial applications, various protein engineering techniques have emerged which allow to tailor the catalytic, biophysical, and molecular recognition properties of enzymes. However, the many possible ways a protein can be altered forces researchers to carefully balance between the exhaustiveness of an enzyme screening campaign and the required resources. Consequently, the optimal engineering strategy is often defined on a case-by-case basis. Strikingly, while predicting mutations that lead to an improved target function is challenging, here we show that the prediction and exclusion of deleterious mutations is a much more straightforward task as analyzed for an engineered carbonic acid anhydrase, a transaminase, a squalene-hopene cyclase and a Kemp eliminase. Combining such a pre-selection of allowed residues with advanced gene synthesis methods opens a path toward an efficient and generalizable library construction approach for protein engineering. To give researchers easy access to this methodology, we provide the website LibGENiE containing the bioinformatic tools for the library design workflow.

Biomanufacturing by In Vitro Biotransformation (ivBT) Using Purified Cascade Multi-enzymes

In vitro biotransformation (ivBT) refers to the use of an artificial biological reaction system that employs purified enzymes for the one-pot conversion of low-cost materials into biocommodities such as ethanol, organic acids, and amino acids. Unshackled from cell growth and metabolism, ivBT exhibits distinct advantages compared with metabolic engineering, including but not limited to high engineering flexibility, ease of operation, fast reaction rate, high product yields, and good scalability. These characteristics position ivBT as a promising next-generation biomanufacturing platform. Nevertheless, challenges persist in the enhancement of bulk enzyme preparation methods, the acquisition of enzymes with superior catalytic properties, and the development of sophisticated approaches for pathway design and system optimization. In alignment with the workflow of ivBT development, this chapter presents a systematic introduction to pathway design, enzyme mining and engineering, system construction, and system optimization. The chapter also proffers perspectives on ivBT development.

Engineering of Cyclodextrin Glycosyltransferase through a Size/Polarity Guided Triple-Code Strategy with Enhanced α-Glycosyl Hesperidin Synthesis Ability

Hesperidin, a flavonoid enriched in citrus peel, can be enzymatically glycosylated using CGTase with significantly improved water solubility. However, the reaction catalyzed by wild-type CGTase is rather inefficient, reflected in the poor production rate and yield. By focusing on the aglycon attacking step, seven residues were selected for mutagenesis in order to improve the transglycosylation efficiency. Due to the lack of high-throughput screening technology regarding to the studied reaction, we developed a size/polarity guided triple-code strategy in order to reduce the library size. The selected residues were replaced by three rationally chosen amino acids with either changed size or polarity, leading to an extremely condensed library with only 32 mutants to be screened. Twenty-five percent of the constructed mutants were proved to be positive, suggesting the high quality of the constructed library. Specific transglycosylation activity of the best mutant Y217F was assayed to be 935.7 U/g, and its k_cat/K_mA is 6.43 times greater than that of the wild type. Homology modeling and docking computation suggest the source of notably enhanced catalytic efficiency is resulted from the combination of ligand transfer and binding effect. IMPORTANCE Size/polarity guided triple-code strategy, a novel semirational mutagenesis strategy, was developed in this study and employed to engineer the aglycon attacking site of CGTase. Screening pressure was set as improved hesperidin glucoside synthesis ability, and eight positive mutants were obtained by screening only 32 mutants. The high quality of the designed library confirms the effectiveness of the developed strategy is potentially valuable to future mutagenesis studies. Mechanisms of positive effect were explained. The best mutant exhibits 6.43 times enhanced k_cat/K_mA value and confirmed to be a superior whole-cell catalyst with potential application value in synthesizing hesperidin glucosides.

Key aspects of the past 30 years of protein design

Proteins are the workhorse of life. They are the building infrastructure of living systems; they are the most efficient molecular machines known, and their enzymatic activity is still unmatched in versatility by any artificial system. Perhaps proteins' most remarkable feature is their modularity. The large amount of information required to specify each protein's function is analogically encoded with an alphabet of just ∼20 letters. The protein folding problem is how to encode all such information in a sequence of 20 letters. In this review, we go through the last 30 years of research to summarize the state of the art and highlight some applications related to fundamental problems of protein evolution.

Insertions and deletions in protein evolution and engineering

Protein evolution or engineering studies are traditionally focused on amino acid substitutions and the way these contribute to fitness. Meanwhile, the insertion and deletion of amino acids is often overlooked, despite being one of the most common sources of genetic variation. Recent methodological advances and successful engineering stories have demonstrated that the time is ripe for greater emphasis on these mutations and their understudied effects. This review highlights the evolutionary importance and biotechnological relevance of insertions and deletions (indels). We provide a comprehensive overview of approaches that can be employed to include indels in random, (semi)-rational or computational protein engineering pipelines. Furthermore, we discuss the tolerance to indels at the structural level, address how domain indels can link the function of unrelated proteins, and feature studies that illustrate the surprising and intriguing potential of frameshift mutations.

Making Enzymes Suitable for Organic Chemistry by Rational Protein Design

This review outlines recent developments in protein engineering of stereo‐ and regioselective enzymes, which are of prime interest in organic and pharmaceutical chemistry as well as biotechnology. The widespread application of enzymes was hampered for decades due to limited enantio‐, diastereo‐ and regioselectivity, which was the reason why most organic chemists were not interested in biocatalysis. This attitude began to change with the advent of semi‐rational directed evolution methods based on focused saturation mutagenesis at sites lining the binding pocket. Screening constitutes the labor‐intensive step (bottleneck), which is the reason why various research groups are continuing to develop techniques for the generation of small and smart mutant libraries. Rational enzyme design, traditionally an alternative to directed evolution, provides small collections of mutants which require minimal screening. This approach first focused on thermostabilization, and did not enter the field of stereoselectivity until later. Computational guides such as the Rosetta algorithms, HotSpot Wizard metric, and machine learning (ML) contribute significantly to decision making. The newest advancements show that semi‐rational directed evolution such as CAST/ISM and rational enzyme design no longer develop on separate tracks, instead, they have started to merge. Indeed, researchers utilizing the two approaches have learned from each other. Today, the toolbox of organic chemists includes enzymes, primarily because the possibility of controlling stereoselectivity by protein engineering has ensured reliability when facing synthetic challenges. This review was also written with the hope that undergraduate and graduate education will include enzymes more so than in the past.

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.

Directed Evolution: Methodologies and Applications

Directed evolution aims to expedite the natural evolution process of biological molecules and systems in a test tube through iterative rounds of gene diversifications and library screening/selection. It has become one of the most powerful and widespread tools for engineering improved or novel functions in proteins, metabolic pathways, and even whole genomes. This review describes the commonly used gene diversification strategies, screening/selection methods, and recently developed continuous evolution strategies for directed evolution. Moreover, we highlight some representative applications of directed evolution in engineering nucleic acids, proteins, pathways, genetic circuits, viruses, and whole cells. Finally, we discuss the challenges and future perspectives in directed evolution.

Flavoprotein monooxygenases: Versatile biocatalysts

Flavoprotein monooxygenases (FPMOs) are single- or two-component enzymes that catalyze a diverse set of chemo-, regio- and enantioselective oxyfunctionalization reactions. In this review, we describe how FPMOs have evolved from model enzymes in mechanistic flavoprotein research to biotechnologically relevant catalysts that can be applied for the sustainable production of valuable chemicals. After a historical account of the development of the FPMO field, we explain the FPMO classification system, which is primarily based on protein structural properties and electron donor specificities. We then summarize the most appealing reactions catalyzed by each group with a focus on the different types of oxygenation chemistries. Wherever relevant, we report engineering strategies that have been used to improve the robustness and applicability of FPMOs.

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.

Computational design of enzymes for biotechnological applications

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

Synthesis of libraries and multi-site mutagenesis using a PCR-derived, dU-containing template

Directed DNA libraries are useful because they focus genetic diversity in the most important regions within a sequence. Ideally, all sequences in such libraries should appear with the same frequency and there should be no significant background from the starting sequence. These properties maximize the number of different sequences that can be screened. Described herein is a method termed SLUPT (Synthesis of Libraries via a dU-containing PCR-derived Template) for generating highly targeted DNA libraries and/or multi-site mutations wherein the altered bases may be widely distributed within a target sequence. This method is highly efficient and modular. Moreover, multiple distinct sites, each with one or more base changes, can be altered in a single reaction. There is very low background from the starting sequence, and SLUPT libraries have similar representation of each base at the positions selected for variation. The SLUPT method utilizes a single-stranded dU-containing DNA template that is made by polymerase chain reaction (PCR). Synthesis of the template in this way is significantly easier than has been described earlier. A series of oligonucleotide primers that are homologous to the template and encode the desired genetic diversity are extended and ligated in a single reaction to form the mutated product sequence or library. After selective inactivation of the template, only the product library is amplified. There are no restrictions on the spacing of the mutagenic primers except that they cannot overlap.

Structure-guided evolution of a ketoreductase for efficient and stereoselective bioreduction of bulky α-amino β-keto esters

Chiral vicinal amino alcohols were generated with excellent stereoselectivity and high conversion from bulky α-amino β-keto esters by an engineered ketoreductase called M30.

Identification two key residues at the intersection of domains of a thioether monooxygenase for improving its sulfoxidation performance

AcCHMO, a cyclohexanone monooxygenase from Acinetobacter calcoaceticus, is a typical Type I Baeyer-Villiger monooxygenase (BVMO). We previously obtained the AcCHMO_M6 mutant, which oxidizes omeprazole sulfide (OPS) to the chiral sulfoxide drug esomeprazole. To further improve the catalytic efficiency of the AcCHMO_M6 mutant, a focused mutagenesis strategy was adopted at the intersections of the FAD-binding domain, NADPH-binding domain, and α-helical domain based on structural characteristics of AcCHMO. By using focused mutagenesis and subsequent global evolution two key residues (L55 and P497) at the intersections of the domains were identified. Mutant of L55Y improved catalytic efficiency significantly, whereas the P497S mutant alleviated substrate inhibition remarkably. AcCHMO_M7 (L55Y/P497S) was obtained by combining the two mutations, which increased the specific activity from 18.5 (M6) to 108 U/g, and an increase in the K_i of the substrate OPS from 34 to 265 μM. The results indicate that catalytic performance can be elevated by modification of the sensitive sites at the intersection of the domains of AcCHMO. The results also provided some insights for the engineering of other Type I BVMOs or other multidomain proteins.

Enhancing promiscuous chemistries of a Schiff-base forming enzyme by divergent evolution

Directed evolution has emerged as a powerful technique for the rapid tailoring of enzymes toward particular synthetic demands, spawning a number of enzymes capable of complex chemical transformations. During random mutagenesis of a protein, changes in fitness must be assayed in order to quantify and understand the relative effect a given mutation has, and the assay employed must be carefully chosen to report directly on the transformation of interest. Here, we describe a series of medium-throughput screening techniques that have been utilized for the evolution and engineering of an artificial carboligase, RA95.5-8, resulting in improvement of catalytic efficiency of a number of promiscuous chemistries. The methods make use of common analytical chemistry equipment and low-cost materials, and may help inspire development of novel screening workflows for related transformations.

Protein design under competing conditions for the availability of amino acids

Isolating the properties of proteins that allow them to convert sequence into the structure is a long-lasting biophysical problem. In particular, studies focused extensively on the effect of a reduced alphabet size on the folding properties. However, the natural alphabet is a compromise between versatility and optimisation of the available resources. Here, for the first time, we include the impact of the relative availability of the amino acids to extract from the 20 letters the core necessary for protein stability. We present a computational protein design scheme that involves the competition for resources between a protein and a potential interaction partner that, additionally, gives us the chance to investigate the effect of the reduced alphabet on protein-protein interactions. We devise a scheme that automatically identifies the optimal reduced set of letters for the design of the protein, and we observe that even alphabets reduced down to 4 letters allow for single protein folding. However, it is only with 6 letters that we achieve optimal folding, thus recovering experimental observations. Additionally, we notice that the binding between the protein and a potential interaction partner could not be avoided with the investigated reduced alphabets. Therefore, we suggest that aggregation could have been a driving force in the evolution of the large protein alphabet.

Efficient Synthesis of Methyl 3-Acetoxypropionate by a Newly Identified Baeyer-Villiger Monooxygenase

Baeyer-Villiger monooxygenases (BVMOs) are an emerging class of promising biocatalysts for the oxidation of ketones to prepare corresponding esters or lactones. Although many BVMOs have been reported, the development of highly efficient enzymes for use in industrial applications is desirable. In this work, we identified a BVMO from Rhodococcus pyridinivorans (BVMO_Rp) with a high affinity toward aliphatic methyl ketones (K_m < 3.0 μM). The enzyme was highly soluble and relatively stable, with a half-life of 23 h at 30°C and pH 7.5. The most effective substrate discovered so far is 2-hexanone (k _cat = 2.1 s^-1; K_m = 1.5 μM). Furthermore, BVMO_Rp exhibited excellent regioselectivity toward most aliphatic ketones, preferentially forming typical (i.e., normal) products. Using the newly identified BVMO_Rp as the catalyst, a high concentration (26.0 g/liter; 200 mM) of methyl levulinate was completely converted to methyl 3-acetoxypropionate after 4 h, with a space-time yield of 5.4 g liter^-1 h^-1 Thus, BVMO_Rp is a promising biocatalyst for the synthesis of 3-hydroxypropionate from readily available biobased levulinate to replace the conventional fermentation.IMPORTANCE BVMOs are emerging as a green alternative to traditional oxidants in the BV oxidation of ketones. Although many BVMOs are discovered and used in organic synthesis, few are really applied in industry, especially in the case of aliphatic ketones. Herein, a highly soluble and relatively stable monooxygenase from Rhodococcus pyridinivorans (BVMO_Rp) was identified with high activity and excellent regioselectivity toward most aliphatic ketones. BVMO_Rp possesses unusually high substrate loading during the catalysis of the oxidation of biobased methyl levulinate to 3-hydroxypropionic acid derivatives. This study indicates that the synthesis of 3-hydroxypropionate from readily available biobased levulinate by BVMO_Rp-catalyzed oxidation holds great promise to replace traditional fermentation.

Genetic Code Evolution Investigated through the Synthesis and Characterisation of Proteins from Reduced-Alphabet Libraries

The universal genetic code of 20 amino acids is the product of evolution. It is believed that earlier versions of the code had fewer residues. Many theories for the order in which amino acids were integrated into the code have been proposed, considering factors ranging from prebiotic chemistry to codon capture. Several meta-analyses combined these theories to yield a feasible consensus chronology of the genetic code's evolution, but there is a dearth of experimental data to test the hypothesised order. We used combinatorial chemistry to synthesise libraries of random polypeptides that were based on different subsets of the 20 standard amino acids, thus representing different stages of a plausible history of the alphabet. Four libraries were comprised of the five, nine, and 16 most ancient amino acids, and all 20 extant residues for a direct side-by-side comparison. We characterised numerous variants from each library for their solubility and propensity to form secondary, tertiary or quaternary structures. Proteins from the two most ancient libraries were more likely to be soluble than those from the extant library. Several individual protein variants exhibited inducible protein folding and other traits typical of intrinsically disordered proteins. From these libraries, we can infer how primordial protein structure and function might have evolved with the genetic code.

In Silico-Directed Evolution Using CADEE

Recent years have seen an explosion of interest in both sequence- and structure-based approaches toward in silico-directed evolution. We recently developed a novel computational toolkit, CADEE, which facilitates the computer-aided directed evolution of enzymes. Our initial work (Amrein et al., IUCrJ 4:50-64, 2017) presented a pedagogical example of the application of CADEE to triosephosphate isomerase, to illustrate the CADEE workflow. In this contribution, we describe this workflow in detail, including code input/output snippets, in order to allow users to set up and execute CADEE simulations on any system of interest.

Improvement of the Activity and Stability of Starch-Debranching Pullulanase from Bacillus naganoensis via Tailoring of the Active Sites Lining the Catalytic Pocket

Pullulanases are well-known debranching enzymes that hydrolyze α-1,6-glycosidic linkages in starch and oligosaccharides. However, most of the pullulanases exhibit limited activity for practical applications. Here, two sites (787 and 621) lining the catalytic pocket of Bacillus naganoensis pullulanase were identified as being critical for enzymatic activity by triple-code saturation mutagenesis. Subsequently, both sites were subjected to NNK-based saturation mutagenesis to obtain positive variants. Among the variants showing enhanced activity, the enzymatic activity and specific activity of D787C were 1.5-fold higher than those of the wild-type (WT). D787C also showed a 1.8-fold increase in k_cat and a 1.7-fold increase in k_cat/ K_m. In addition, D787C maintained higher activity compared with that of WT at temperatures over 60 °C. All the positive variants showed higher acid resistance, with D787C maintaining 90% residual activity at pH 4.0. Thus, enzymes with improved properties were obtained by saturation mutagenesis at the active site.

Side-Chain Pruning Has Limited Impact on Substrate Preference in a Promiscuous Enzyme

Detoxifying enzymes such as flavin-containing monooxygenases deal with a huge array of highly diverse xenobiotics and toxic compounds. In addition to being of high physiological relevance, these drug-metabolizing enzymes are useful catalysts for synthetic chemistry. Despite the wealth of studies, the molecular basis of their relaxed substrate selectivity remains an open question. Here, we addressed this issue by applying a cumulative alanine mutagenesis approach to cyclohexanone monooxygenase from Thermocrispum municipale, a flavin-dependent Baeyer–Villiger monooxygenase which we chose as a model system because of its pronounced thermostability and substrate promiscuity. Simultaneous removal of up to eight noncatalytic active-site side chains including four phenylalanines had no effect on protein folding, thermostability, and cofactor loading. We observed a linear decrease in activity, rather than a selectivity switch, and attributed this to a less efficient catalytic environment in the enlarged active-site space. Time-resolved kinetic studies confirmed this interpretation. We also determined the crystal structure of the enzyme in complex with a mimic of the reaction intermediate that shows an unaltered overall protein conformation. These findings led us to propose that this cyclohexanone monooxygenase may lack a distinct substrate selection mechanism altogether. We speculate that the main or exclusive function of the protein shell in promiscuous enzymes might be the stabilization and accessibility of their very reactive catalytic intermediates.

Kinetic characterization of acetone monooxygenase from Gordonia sp. strain TY-5

Acetone monooxygenase (ACMO) is a unique member of the Baeyer–Villiger monooxygenase (BVMO) family based on its ability to act on small ketones, such as acetone. Herein, we performed a kinetic analysis of ACMO from the propane-utilizing bacterium Gordonia sp. strain TY-5 to assess its preference for smaller ketone substrates. Steady state kinetic analysis of ACMO with a range of linear (C3–C7) and cyclic ketones (C4–C6) reveals that the enzyme elicits the highest catalytic efficiency towards butanone and cyclobutanone. Stopped-flow and inhibition studies further revealed that ACMO has a relatively weak binding affinity for the coenzyme with a dissociation constant of 120 μM. We show through mutagenesis that sequence variation in the residue that coordinates to the 2′-phosphate of NADP(H) partially accounts for the weaker binding affinity observed. As for shown for related BVMOs, NADP⁺ stabilizes the C4a-peroxyflavin intermediate in ACMO; however, the rate of its formation is considerably slower in ACMO. The observed rate constant for NADPH-dependent flavin reduction was 60 s⁻¹ at 25 °C, and experiments performed with 4(R)-[4-²H]NADPH confirm that the C4-pro-R-hydride from the nicotinamide ring is transferred to the FAD. The latter experimental result suggests that the nicotinamide ring rotates within the active site to carry out its two functional roles: reduction of the FAD cofactor and stabilization of the C4a-peroxyflavin adduct.

Resolution of Atropisomeric Cyclic Catechol Monoether O-Sulfate Esters by a Molluscan Sulfatase

Atropisomeric cyclic catechol ethers are notoriously difficult to resolve by classical chiral phase high-performance liquid chromatography. Here, we show the first application of sulfatase enzymes for the kinetic resolution of O-sulfato-catechol ethers with enantioselectivities ranging from 30 to 65% ee, as determined by preparation of their Marfey's ether derivatives. Substrate-structure dependence was briefly explored.

Inducing high activity of a thermophilic enzyme at ambient temperatures by directed evolution

The long-standing problem of achieving high activity of a thermophilic enzyme at low temperatures and short reaction times with little tradeoff in thermostability has been solved by directed evolution, an alcohol dehydrogenase found in hot springs serving as the catalyst in enantioselective ketone reductions.

Spatial requirement for PAMO for transformation of non-native linear substrates

Phenylacetone monooxygenase is the most stable and thermo-tolerant member of the Baeyer-Villiger monooxygenases family, and therefore it is an ideal candidate for the synthesis of industrially relevant ester or lactone compounds. However, its limited substrate scope has largely limited its industrial applications. Linear substrates are interesting from an industrial point of view, it is thus necessary to identify the essential spatial requirement for achieving high conversions for non-native linear substrates. Here using molecular dynamics simulations, we compared the conversion of a non-native linear substrate 2-octanone and the native substrate phenylacetone, catalyzed by the WT enzyme and a quadruple variant P253F/G254A/R258M/L443F that exhibits significantly improved activity towards 2-octanone. We uncovered that a remarkable movement of L289 is crucial for a reshaping of the active site of the quadruple variant so as to prevent the aliphatic substrate from moving away from the C4a-peroxyflavin, thus enabling it to keep a catalytically relevant pose during the oxygenation process. By performing steady-state kinetic analysis of two single-mutation variants at position 258, we further validated that the L289 reposition is attributed to the combined effect of quadruple mutations. In order to further explore the substrate scope of PAMO we also studied the binding of cyclopentanone and 2-phenylcyclohexanone, which are the typical substrates of CPMO in group I and CHMO in group III, respectively. Our study provides fundamental atomic-level insights in rational engineering of PAMO for wide applications in industrial biocatalysis, in particular, in the biotransformation of long-chain aliphatic oils into potential biodiesels.

Catalytic mechanism of phenylacetone monooxygenases for non-native linear substrates

Phenylacetone monooxygenase (PAMO) is the most stable and thermo-tolerant member of the Baeyer-Villiger monooxygenase family, and therefore it is an ideal candidate for the synthesis of industrially relevant compounds. However, its limited substrate scope has largely limited its industrial applications. In the present work, we provide, for the first time, the catalytic mechanism of PAMO for the native substrate phenylacetone as well as for a linear non-native substrate 2-octanone, using molecular dynamics simulations, quantum mechanics and quantum mechanics/molecular mechanics calculations. We provide a theoretical basis for the preference of the enzyme for the native aromatic substrate over non-native linear substrates. Our study provides fundamental atomic-level insights that can be employed in the rational engineering of PAMO for wide applications in industrial biocatalysis, in particular, in the biotransformation of long-chain aliphatic oils into potential biodiesels.

Reconstruction of cysteine biosynthesis using engineered cysteine-free enzymes

Amino acid biosynthesis pathways observed in nature typically require enzymes that are made with the amino acids they produce. For example, Escherichia coli produces cysteine from serine via two enzymes that contain cysteine: serine acetyltransferase (CysE) and O-acetylserine sulfhydrylase (CysK/CysM). To solve this chicken-and-egg problem, we substituted alternate amino acids in CysE, CysK and CysM for cysteine and methionine, which are the only two sulfur-containing proteinogenic amino acids. Using a cysteine-dependent auxotrophic E. coli strain, CysE function was rescued by cysteine-free and methionine-deficient enzymes, and CysM function was rescued by cysteine-free enzymes. CysK function, however, was not rescued in either case. Enzymatic assays showed that the enzymes responsible for rescuing the function in CysE and CysM also retained their activities in vitro. Additionally, substitution of the two highly conserved methionines in CysM decreased but did not eliminate overall activity. Engineering amino acid biosynthetic enzymes to lack the so-produced amino acids can provide insights into, and perhaps eventually fully recapitulate via a synthetic approach, the biogenesis of biotic amino acids.

Biocatalysts for the pharmaceutical industry created by structure-guided directed evolution of stereoselective enzymes

Enzymes have been used for a long time as catalysts in the asymmetric synthesis of chiral intermediates needed in the production of therapeutic drugs. However, this alternative to man-made catalysts has suffered traditionally from distinct limitations, namely the often observed wrong or insufficient enantio- and/or regioselectivity, low activity, narrow substrate range, and insufficient thermostability. With the advent of directed evolution, these problems can be generally solved. The challenge is to develop and apply the most efficient mutagenesis methods which lead to highest-quality mutant libraries requiring minimal screening. Structure-guided saturation mutagenesis and its iterative form have emerged as the method of choice for evolving stereo- and regioselective mutant enzymes needed in the asymmetric synthesis of chiral intermediates. The number of (industrial) applications in the preparation of chiral pharmaceuticals is rapidly increasing. This review features and analyzes typical case studies.

Development of Halogenase Enzymes for Use in Synthesis

Nature has evolved halogenase enzymes to regioselectively halogenate a diverse range of biosynthetic precursors, with the halogens introduced often having a profound effect on the biological activity of the resulting natural products. Synthetic endeavors to create non-natural bioactive small molecules for pharmaceutical and agrochemical applications have also arrived at a similar conclusion: halogens can dramatically improve the properties of organic molecules for selective modulation of biological targets in vivo. Consequently, a high proportion of pharmaceuticals and agrochemicals on the market today possess halogens. Halogenated organic compounds are also common intermediates in synthesis and are particularly valuable in metal-catalyzed cross-coupling reactions. Despite the potential utility of organohalogens, traditional nonenzymatic halogenation chemistry utilizes deleterious reagents and often lacks regiocontrol. Reliable, facile, and cleaner methods for the regioselective halogenation of organic compounds are therefore essential in the development of economical and environmentally friendly industrial processes. A potential avenue toward such methods is the use of halogenase enzymes, responsible for the biosynthesis of halogenated natural products, as biocatalysts. This Review will discuss advances in developing halogenases for biocatalysis, potential untapped sources of such biocatalysts and how further optimization of these enzymes is required to achieve the goal of industrial scale biohalogenation.

A General Tool for Engineering the NAD/NADP Cofactor Preference of Oxidoreductases

The ability to control enzymatic nicotinamide cofactor utilization is critical for engineering efficient metabolic pathways. However, the complex interactions that determine cofactor-binding preference render this engineering particularly challenging. Physics-based models have been insufficiently accurate and blind directed evolution methods too inefficient to be widely adopted. Building on a comprehensive survey of previous studies and our own prior engineering successes, we present a structure-guided, semirational strategy for reversing enzymatic nicotinamide cofactor specificity. This heuristic-based approach leverages the diversity and sensitivity of catalytically productive cofactor binding geometries to limit the problem to an experimentally tractable scale. We demonstrate the efficacy of this strategy by inverting the cofactor specificity of four structurally diverse NADP-dependent enzymes: glyoxylate reductase, cinnamyl alcohol dehydrogenase, xylose reductase, and iron-containing alcohol dehydrogenase. The analytical components of this approach have been fully automated and are available in the form of an easy-to-use web tool: Cofactor Specificity Reversal-Structural Analysis and Library Design (CSR-SALAD).

Manipulating the stereoselectivity of the thermostable Baeyer–Villiger monooxygenase TmCHMO by directed evolution

The thermostable Baeyer–Villiger monooxygenase TmCHMO and evolved mutants are viable catalysts in stereoselective reactions of structurally different ketones.

Enzymatic site-selectivity enabled by structure-guided directed evolution

This review covers recent advances in the directed evolution of enzymes for controlling site-selectivity of hydroxylation, amination and chlorination.

Substitution of a Single Amino Acid Reverses the Regiospecificity of the Baeyer-Villiger Monooxygenase PntE in the Biosynthesis of the Antibiotic Pentalenolactone

In the biosynthesis of pentalenolactone (1), PenE and PntE, orthologous proteins from Streptomyces exfoliatus and S. arenae, respectively, catalyze the flavin-dependent Baeyer-Villiger oxidation of 1-deoxy-11-oxopentalenic acid (4) to the lactone pentalenolactone D (5), in which the less-substituted methylene carbon has migrated. By contrast, the paralogous PtlE enzyme from S. avermitilis catalyzes the oxidation of 4 to neopentalenolactone D (6), in which the more substituted methane substitution has undergone migration. We report the design and analysis of 13 single and multiple mutants of PntE mutants to identify the key amino acids that contribute to the regiospecificity of these two classes of Baeyer-Villiger monooxygenases. The L185S mutation in PntE reversed the observed regiospecificity of PntE such that all recombinant PntE mutants harboring this L185S mutation acquired the characteristic regiospecificity of PtlE, catalyzing the conversion of 4 to 6 as the major product. The recombinant PntE mutant harboring R484L exhibited reduced regiospecificity, generating a mixture of lactones containing more than 17% of 6. These in vitro results were corroborated by analysis of the complementation of the S. avermitilis ΔptlED double deletion mutant with pntE mutants, such that pntE mutants harboring L185S produced 6 as the major product, whereas complemention of the ΔptlED deletion mutant with pntE mutants carrying the R484L mutation gave 6 as more than 33% of the total lactone product mixture.

Facilitating the Evolution of Esterase Activity from a Promiscuous Enzyme (Mhg) with Catalytic Functions of Amide Hydrolysis and Carboxylic Acid Perhydrolysis by Engineering the Substrate Entrance Tunnel

Promiscuous enzymes are generally considered to be starting points in the evolution of offspring enzymes with more specific or even novel catalytic activities, which is the molecular basis of producing new biological functions. Mhg, a typical α/β fold hydrolase, was previously reported to have both γ-lactamase and perhydrolase activities. However, despite having high structural similarity to and sharing an identical catalytic triad with an extensively studied esterase from Pseudomonas fluorescens, this enzyme did not show any esterase activity. Molecular docking and sequence analysis suggested a possible role for the entry of the binding pocket in blocking the entrance tunnel, preventing the ester compounds from entering into the pocket. By engineering the entrance tunnel with only one or two amino acid substitutions, we successfully obtained five esterase variants of Mhg. The variants exhibited a very broad substrate acceptance, hydrolyzing not only the classical p-nitrophenol esters but also various types of chiral esters, which are widely used as drug intermediates. Site 233 at the entrance tunnel of Mhg was found to play a pivotal role in modulating the three catalytic activities by adjusting the size and shape of the tunnel, with different amino acid substitutions at this site facilitating different activities. Remarkably, the variant with the L233G mutation was a very specific esterase without any γ-lactamase and perhydrolase activities. Considering the amino acid conservation and differentiation, this site could be a key target for future protein engineering. In addition, we demonstrate that engineering the entrance tunnel is an efficient strategy to regulate enzyme catalytic capabilities.

IMPORTANCE

Promiscuous enzymes can act as starting points in the evolution of novel catalytic activities, thus providing a molecular basis for the production of new biological functions. In this study, we identified a critical amino acid residue (Leu233) at the entry of the substrate tunnel of a promiscuous enzyme, Mhg. We found that substitution of this residue with smaller amino acids such as Gly, Ala, Ser, or Pro endowed the enzyme with novel esterase activity. Different amino acids at this site can facilitate different catalytic activities. These findings exhibited universal significance in this subset of α/β fold hydrolases, including Mhg. Furthermore, we demonstrate that engineering the entrance tunnel is an efficient strategy to evolve new enzyme catalytic capabilities. Our study has important implications for the regulation of enzyme catalytic promiscuity and development of protein engineering methodologies.

Structural Analysis of an Evolved Transketolase Reveals Divergent Binding Modes

The S385Y/D469T/R520Q variant of E. coli transketolase was evolved previously with three successive smart libraries, each guided by different structural, bioinformatical or computational methods. Substrate-walking progressively shifted the target acceptor substrate from phosphorylated aldehydes, towards a non-phosphorylated polar aldehyde, a non-polar aliphatic aldehyde, and finally a non-polar aromatic aldehyde. Kinetic evaluations on three benzaldehyde derivatives, suggested that their active-site binding was differentially sensitive to the S385Y mutation. Docking into mutants generated in silico from the wild-type crystal structure was not wholly satisfactory, as errors accumulated with successive mutations, and hampered further smart-library designs. Here we report the crystal structure of the S385Y/D469T/R520Q variant, and molecular docking of three substrates. This now supports our original hypothesis that directed-evolution had generated an evolutionary intermediate with divergent binding modes for the three aromatic aldehydes tested. The new active site contained two binding pockets supporting π-π stacking interactions, sterically separated by the D469T mutation. While 3-formylbenzoic acid (3-FBA) preferred one pocket, and 4-FBA the other, the less well-accepted substrate 3-hydroxybenzaldehyde (3-HBA) was caught in limbo with equal preference for the two pockets. This work highlights the value of obtaining crystal structures of evolved enzyme variants, for continued and reliable use of smart library strategies.

A part toolbox to tune genetic expression in Bacillus subtilis

Libraries of well-characterised components regulating gene expression levels are essential to many synthetic biology applications. While widely available for the Gram-negative model bacterium Escherichia coli, such libraries are lacking for the Gram-positive model Bacillus subtilis, a key organism for basic research and biotechnological applications. Here, we engineered a genetic toolbox comprising libraries of promoters, Ribosome Binding Sites (RBS), and protein degradation tags to precisely tune gene expression in B. subtilis. We first designed a modular Expression Operating Unit (EOU) facilitating parts assembly and modifications and providing a standard genetic context for gene circuits implementation. We then selected native, constitutive promoters of B. subtilis and efficient RBS sequences from which we engineered three promoters and three RBS sequence libraries exhibiting ∼14 000-fold dynamic range in gene expression levels. We also designed a collection of SsrA proteolysis tags of variable strength. Finally, by using fluorescence fluctuation methods coupled with two-photon microscopy, we quantified the absolute concentration of GFP in a subset of strains from the library. Our complete promoters and RBS sequences library comprising over 135 constructs enables tuning of GFP concentration over five orders of magnitude, from 0.05 to 700 μM. This toolbox of regulatory components will support many research and engineering applications in B. subtilis.