Induced allostery in the directed evolution of an enantioselective Baeyer-Villiger monooxygenase

Engineering of a Baeyer-Villiger monooxygenase to Improve Substrate Scope, Stereoselectivity and Regioselectivity

Baeyer-Villiger monooxygenases belong to a family of flavin-binding proteins that catalyze the Baeyer-Villiger (BV) oxidation of ketones to produce lactones or esters, which are important intermediates in pharmaceuticals or sustainable materials. Phenylacetone monooxygenase (PAMO) from Thermobifida fusca with moderate thermostability catalyzes the oxidation of aryl ketone substrates, but is limited by high specificity and narrow substrate scope. In the present study, we applied loop optimization by loop swapping followed by focused saturation mutagenesis in order to evolve PAMO mutants capable of catalyzing the regioselective BV oxidation of cyclohexanone and cyclobutanone derivatives with formation of either normal or abnormal esters or lactones. We further modulated PAMO to increase enantioselectivity. Crystal structure studies indicate that rotation occurs in the NADP-binding domain and that the high B-factor region is predominantly distributed in the catalytic pocket residues. Computational analyses further revealed dynamic character in the catalytic pocket and reshaped hydrogen bond interaction networks, which is more favorable for substrate binding. Our study provides useful insights for studying enzyme-substrate adaptations.

Protein Allostery Study in Cells Using NMR Spectroscopy

Protein allostery is commonly observed in vitro. But how protein allostery behaves in cells is unknown. In this work, a protein monomer-dimer equilibrium system was built with the allosteric effect on the binding characterized using NMR spectroscopy through mutations away from the dimer interface. A chemical shift linear fitting method was developed that enabled us to accurately determine the dissociation constant. A total of 28 allosteric mutations were prepared and grouped to negative allosteric, nonallosteric, and positive allosteric modulators. ∼ 50% of mutations displayed the allosteric-state changes when moving from a buffered solution into cells. For example, there were no positive allosteric modulators in the buffered solution but eight in cells. The change in protein allostery is correlated with the interactions between the protein and the cellular environment. These interactions presumably drive the surrounding macromolecules in cells to transiently bind to the monomer and dimer mutational sites and change the free energies of the two species differently which generate new allosteric effects. These surrounding macromolecules create a new protein allostery pathway that is only present in cells.

Engineering Rieske oxygenase activity one piece at a time

Enzyme engineering plays a central role in the development of biocatalysts for biotechnology, chemical and pharmaceutical manufacturing, and environmental remediation. Rational design of proteins has historically relied on targeting active site residues to confer a protein with desirable catalytic properties. However, additional "hotspots" are also known to exist beyond the active site. Structural elements such as subunit-subunit interactions, entrance tunnels, and flexible loops influence enzyme catalysis and serve as potential "hotspots" for engineering. For the Rieske oxygenases, which use a Rieske cluster and mononuclear iron center to catalyze a challenging set of reactions, these outside of the active site regions are increasingly being shown to drive catalytic outcomes. Therefore, here, we highlight recent work on structurally characterized Rieske oxygenases that implicates architectural pieces inside and outside of the active site as key dictators of catalysis, and we suggest that these features may warrant attention in efforts aimed at Rieske oxygenase engineering.

Discovery of New Phenylacetone Monooxygenase Variants for the Development of Substituted Indigoids through Biocatalysis

Indigoids are natural pigments obtained from plants by ancient cultures. Romans used them mainly as dyes, whereas Asian cultures applied these compounds as treatment agents for several diseases. In the modern era, the chemical industry has made it possible to identify and develop synthetic routes to obtain them from petroleum derivatives. However, these processes require high temperatures and pressures and large amounts of solvents, acids, and alkali agents. Thus, enzyme engineering and the development of bacteria as whole-cell biocatalysts emerges as a promising green alternative to avoid the use of these hazardous materials and consequently prevent toxic waste generation. In this research, we obtained two novel variants of phenylacetone monooxygenase (PAMO) by iterative saturation mutagenesis. Heterologous expression of these two enzymes, called PAMOHPCD and PAMOHPED, in E. coli was serendipitously found to produce indigoids. These interesting results encourage us to characterize the thermal stability and enzyme kinetics of these new variants and to evaluate indigo and indirubin production in a whole-cell system by HPLC. The highest yields were obtained with PAMOHPCD supplemented with L-tryptophan, producing ~3000 mg/L indigo and ~130.0 mg/L indirubin. Additionally, both enzymes could oxidize and produce several indigo derivatives from substituted indoles, with PAMOHPCD being able to produce the well-known Tyrian purple. Our results indicate that the PAMO variants described herein have potential application in the textile, pharmaceutics, and semiconductors industries, prompting the use of environmentally friendly strategies to obtain a diverse variety of indigoids.

Redesigning Enzymes for Biocatalysis: Exploiting Structural Understanding for Improved Selectivity

The discovery of new enzymes, alongside the push to make chemical processes more sustainable, has resulted in increased industrial interest in the use of biocatalytic processes to produce high-value and chiral precursor chemicals. Huge strides in protein engineering methodology and in silico tools have facilitated significant progress in the discovery and production of enzymes for biocatalytic processes. However, there are significant gaps in our knowledge of the relationship between enzyme structure and function. This has demonstrated the need for improved computational methods to model mechanisms and understand structure dynamics. Here, we explore efforts to rationally modify enzymes toward changing aspects of their catalyzed chemistry. We highlight examples of enzymes where links between enzyme function and structure have been made, thus enabling rational changes to the enzyme structure to give predictable chemical outcomes. We look at future directions the field could take and the technologies that will enable it.

Leveraging intrinsic flexibility to engineer enhanced enzyme catalytic activity

Dynamic motions of enzymes occurring on a broad range of timescales play a pivotal role in all steps of the reaction pathway, including substrate binding, catalysis, and product release. However, it is unknown whether structural information related to conformational flexibility can be exploited for the directed evolution of enzymes with higher catalytic activity. Here, we show that mutagenesis of residues exclusively located at flexible regions distal to the active site of Homo sapiens kynureninase (HsKYNase) resulted in the isolation of a variant (BF-HsKYNase) in which the rate of the chemical step toward kynurenine was increased by 45-fold. Mechanistic pre–steady-state kinetic analysis of the wild type and the evolved enzyme shed light on the underlying effects of distal mutations (>10 Å from the active site) on the rate-limiting step of the catalytic cycle. Hydrogen-deuterium exchange coupled to mass spectrometry and molecular dynamics simulations revealed that the amino acid substitutions in BF-HsKYNase allosterically affect the flexibility of the pyridoxal-5′-phosphate (PLP) binding pocket, thereby impacting the rate of chemistry, presumably by altering the conformational ensemble and sampling states more favorable to the catalyzed reaction.

Unlocking the Stereoselectivity and Substrate Acceptance of Enzymes: Proline-Induced Loop Engineering Test

Protein stability and evolvability influence each other. Although protein dynamics play essential roles in various catalytically important properties, their high flexibility and diversity makes it difficult to incorporate such properties into rational engineering. Therefore, how to unlock the potential evolvability in a user-friendly rational design process remains a challenge. In this endeavor, we describe a method for engineering an enantioselective alcohol dehydrogenase. It enables synthetically important substrate acceptance for 4-chlorophenyl pyridine-2-yl ketone, and perfect stereocontrol of both (S)- and (R)-configured products. Thermodynamic analysis unveiled the subtle interaction between enzyme stability and evolvability, while computational studies provided insights into the origin of selectivity and substrate recognition. Preparative-scale synthesis of the (S)-product (73 % yield; >99 % ee) was performed on a gram-scale. This proof-of-principle study demonstrates that interfaced proline residues can be rationally engineered to unlock evolvability and thus provide access to new biocatalysts with highly improved catalytic performance.

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.

Enzymatic Kinetic Resolution by Addition of Oxygen

Kinetic resolution using biocatalysis has proven to be an excellent complementary technique to traditional asymmetric catalysis for the production of enantioenriched compounds. Resolution using oxidative enzymes produces valuable oxygenated structures for use in synthetic route development. This Minireview focuses on enzymes which catalyse the insertion of an oxygen atom into the substrate and, in so doing, can achieve oxidative kinetic resolution. The Baeyer-Villiger rearrangement, epoxidation, and hydroxylation are included, and biological advancements in enzyme development, and applications of these key enantioenriched intermediates in natural product synthesis are discussed.

Symmetry-related residues as promising hotspots for the evolution of de novo oligomeric enzymes

Directed evolution has provided us with great opportunities and prospects in the synthesis of tailor-made proteins. It, however, often requires at least mid to high throughput screening, necessitating more effective strategies for laboratory evolution. We herein demonstrate that protein symmetry can be a versatile criterion for searching for promising hotspots for the directed evolution of de novo oligomeric enzymes. The randomization of symmetry-related residues located at the rotational axes of artificial metallo-β-lactamase yields drastic effects on catalytic activities, whereas that of non-symmetry-related, yet, proximal residues to the active site results in negligible perturbations. Structural and biochemical analysis of the positive hits indicates that seemingly trivial mutations at symmetry-related spots yield significant alterations in overall structures, metal-coordination geometry, and chemical environments of active sites. Our work implicates that numerous artificially designed and natural oligomeric proteins might have evolutionary advantages of propagating beneficial mutations using their global symmetry.

Symmetry-related residues located at the rotational axes can be promising hotspots for the evolution of de novo oligomeric enzymes even though they are distantly located from the active site pocket.

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.

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.

The mutagenesis of a single site for enhancing or reversing the enantio- or regiopreference of cyclohexanone monooxygenases

The mutagenesis of a "second sphere" switch residue of CHMO_Acineto could control its enantio- and regiopreference. Replacing phenylalanine (F) at position 277 of CHMO_Acineto into larger tryptophan (W) enabled a significant enhancement of enantio- or regioselectivity toward structurally diverse substrates, moreover, a complete reversal of enantio- or regiopreference was realized by mutating F277 into a range of smaller amino acids (A/C/D/E/G/H/I/K/L/M/N/P/Q/R/S/T/V).

Natural Variation in the ‘Control Loop’ of BVMOAFL210 and Its Influence on Regioselectivity and Sulfoxidation

Baeyer-Villiger monooxygenases (BVMOs) are flavin-dependent enzymes that primarily convert ketones to esters, but can also catalyze heteroatom oxidation. Several structural studies have highlighted the importance of the ‘control loop’ in BVMOs, which adopts different conformations during catalysis. Central to the ‘control loop’ is a conserved tryptophan that has been implicated in NADP(H) binding. BVMOAFL210 from Aspergillus flavus, however, contains a threonine in the equivalent position. Here, we report the structure of BVMOAFL210 in complex with NADP+ in both the ‘open’ and ‘closed’ conformations. In neither conformation does Thr513 contact the NADP+. Although mutagenesis of Thr513 did not significantly alter the substrate scope, changes in peroxyflavin stability and reaction rates were observed. Mutation of this position also brought about changes in the regio- and enantioselectivity of the enzyme. Moreover, lower rates of overoxidation during sulfoxidation of thioanisole were also observed.

Functional tunability from a distance: Rheostat positions influence allosteric coupling between two distant binding sites

For protein mutagenesis, a common expectation is that important positions will behave like on/off “toggle” switches (i.e., a few substitutions act like wildtype, most abolish function). However, there exists another class of important positions that manifests a wide range of functional outcomes upon substitution: “rheostat” positions. Previously, we evaluated rheostat positions located near the allosteric binding sites for inhibitor alanine (Ala) and activator fructose-1,6-bisphosphate (Fru-1,6-BP) in human liver pyruvate kinase. When substituted with multiple amino acids, many positions demonstrated moderate rheostatic effects on allosteric coupling between effector binding and phosphoenolpyruvate (PEP) binding in the active site. Nonetheless, the combined outcomes of all positions sampled the full range of possible allosteric coupling (full tunability). However, that study only evaluated allosteric tunability of “local” positions, i.e., positions were located near the binding sites of the allosteric ligand being assessed. Here, we evaluated tunability of allosteric coupling when mutated sites were distant from the allosterically-coupled binding sites. Positions near the Ala binding site had rheostatic outcomes on allosteric coupling between Fru-1,6-BP and PEP binding. In contrast, positions in the Fru-1,6-BP site exhibited modest effects on coupling between Ala and PEP binding. Analyzed in aggregate, both PEP/Ala and PEP/Fru-1,6-BP coupling were again fully tunable by amino acid substitutions at this limited set of distant positions. Furthermore, some positions exhibited rheostatic control over multiple parameters and others exhibited rheostatic effects on one parameter and toggle control over a second. These findings highlight challenges in efforts to both predict/interpret mutational outcomes and engineer functions into proteins.

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.

Alternative Splicing of the Aflatoxin-Associated Baeyer⁻Villiger Monooxygenase from Aspergillus flavus: Characterisation of MoxY Isoforms

Aflatoxins are carcinogenic mycotoxins that are produced by the filamentous fungus Aspergillus flavus, a contaminant of numerous food crops. Aflatoxins are synthesised via the aflatoxin biosynthesis pathway, with the enzymes involved encoded by the aflatoxin biosynthesis gene cluster. MoxY is a type I Baeyer–Villiger monooxygenase (BVMO), responsible for the conversion of hydroxyversicolorone (HVN) and versicolorone (VN) to versiconal hemiacetal acetate (VHA) and versiconol acetate (VOAc), respectively. Using mRNA data, an intron near the C-terminus was identified that is alternatively spliced, creating two possible MoxY isoforms which exist in vivo, while analysis of the genomic DNA suggests an alternative start codon leading to possible elongation of the N-terminus. These four variants of the moxY gene were recombinantly expressed in Escherichia coli, and their activity evaluated with respect to their natural substrates HVN and VN, as well as surrogate ketone substrates. Activity of the enzyme is absolutely dependent on the additional 22 amino acid residues at the N-terminus. Two MoxY isoforms with alternative C-termini, MoxYAltN and MoxYAltNC, converted HVN and VN, in addition to a range of ketone substrates. Stability and flavin-binding data suggest that MoxYAltN is, most likely, the dominant isoform. MoxYAltNC is generated by intron splicing, in contrast to intron retention, which is the most prevalent type of alternative splicing in ascomycetes. The alternative C-termini did not alter the substrate acceptance profile, or regio- or enantioselectivity of the enzyme, but did significantly affect the solubility and stability.

Directed Evolution Mimics Allosteric Activation by Stepwise Tuning of the Conformational Ensemble

Allosteric enzymes contain a wealth of catalytic diversity that remains distinctly underutilized for biocatalysis. Tryptophan synthase is a model allosteric system and a valuable enzyme for the synthesis of noncanonical amino acids (ncAA). Previously, we evolved the β-subunit from Pyrococcus furiosus, PfTrpB, for ncAA synthase activity in the absence of its native partner protein PfTrpA. However, the precise mechanism by which mutation activated TrpB to afford a stand-alone catalyst remained enigmatic. Here, we show that directed evolution caused a gradual change in the rate-limiting step of the catalytic cycle. Concomitantly, the steady-state distribution of the intermediates shifts to favor covalently bound Trp adducts, which have increased thermodynamic stability. The biochemical properties of these evolved, stand-alone TrpBs converge on those induced in the native system by allosteric activation. High-resolution crystal structures of the wild-type enzyme, an intermediate in the lineage, and the final variant, encompassing five distinct chemical states, show that activating mutations have only minor structural effects on their immediate environment. Instead, mutation stabilizes the large-scale motion of a subdomain to favor an otherwise transiently populated closed conformational state. This increase in stability enabled the first structural description of Trp covalently bound in a catalytically active TrpB, confirming key features of catalysis. These data combine to show that sophisticated models of allostery are not a prerequisite to recapitulating its complex effects via directed evolution, opening the way to engineering stand-alone versions of diverse allosteric enzymes.

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.

Directed Evolution of Proteins Based on Mutational Scanning

Directed evolution has emerged as one of the most effective protein engineering methods in basic research as well as in applications in synthetic organic chemistry and biotechnology. The successful engineering of protein activity, allostery, binding affinity, expression, folding, fluorescence, solubility, substrate scope, selectivity (enantio-, stereo-, and regioselectivity), and/or stability (temperature, organic solvents, pH) is just limited by the throughput of the genetic selection, display, or screening system that is available for a given protein. Sometimes it is possible to analyze millions of protein variants from combinatorial libraries per day. In other cases, however, only a few hundred variants can be screened in a single day, and thus the creation of smaller yet smarter libraries is needed. Different strategies have been developed to create these libraries. One approach is to perform mutational scanning or to construct "mutability landscapes" in order to understand sequence-function relationships that can guide the actual directed evolution process. Herein we provide a protocol for economically constructing scanning mutagenesis libraries using a cytochrome P450 enzyme in a high-throughput manner. The goal is to engineer activity, regioselectivity, and stereoselectivity in the oxidative hydroxylation of a steroid, a challenging reaction in synthetic organic chemistry. Libraries based on mutability landscapes can be used to engineer any fitness trait of interest. The protocol is also useful for constructing gene libraries for deep mutational scanning experiments.

Reply to the Comment by Littlechild and Isupov

I thank Drs. Littlechild and Isupov for their recent comments, which are considered below. Before addressing these specifically, their correspondence raises two more general issues which require initial clarification.[...].

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.

One-step combined focused epPCR and saturation mutagenesis for thermostability evolution of a new cold-active xylanase

Enzymes active at low temperature are of great interest for industrial bioprocesses due to their high efficiency at a low energy cost. One of the particularities of naturally evolved cold-active enzymes is their increased enzymatic activity at low temperature, however the low thermostability presented in this type of enzymes is still a major drawback for their application in biocatalysis. Directed evolution of cold-adapted enzymes to a more thermostable version, appears as an attractive strategy to fulfill the stability and activity requirements for the industry. This paper describes the recombinant expression and characterization of a new and highly active cold-adapted xylanase from the GH-family 10 (Xyl-L), and the use of a novel one step combined directed evolution technique that comprises saturation mutagenesis and focused epPCR as a feasible semi-rational strategy to improve the thermostability. The Xyl-L enzyme was cloned from a marine-Antarctic bacterium, Psychrobacter sp. strain 2-17, recombinantly expressed in E. coli strain BL21(DE3) and characterized enzymatically. Molecular dynamic simulations using a homology model of the catalytic domain of Xyl-L were performed to detect flexible regions and residues, which are considered to be the possible structural elements that define the thermolability of this enzyme. Mutagenic libraries were designed in order to stabilize the protein introducing mutations in some of the flexible regions and residues identified. Twelve positive mutant clones were found to improve the T₅₀¹⁵ value of the enzyme, in some cases without affecting the activity at 25°C. The best mutant showed a 4.3°C increase in its T₅₀¹⁵. The efficiency of the directed evolution approach can also be expected to work in the protein engineering of stereoselectivity.

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.

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.

The role of side chain entropy and mutual information for improving the de novo design of Kemp eliminases KE07 and KE70

Side chain entropy and mutual entropy information between residue pairs have been calculated for two de novo designed Kemp eliminase enzymes, KE07 and KE70, and for their most improved versions at the end of laboratory directed evolution (LDE).

Protein engineering: a new frontier for biological therapeutics

Protein engineering holds the potential to transform the metabolic drug landscape through the development of smart, stimulusresponsive drug systems. Protein therapeutics are a rapidly expanding segment of Food and Drug Administration approved drugs that will improve clinical outcomes over the long run. Engineering of protein therapeutics is still in its infancy, but recent general advances in protein engineering capabilities are being leveraged to yield improved control over both pharmacokinetics and pharmacodynamics. Stimulus- responsive protein therapeutics are drugs which have been designed to be metabolized under targeted conditions. Protein engineering is being utilized to develop tailored smart therapeutics with biochemical logic. This review focuses on applications of targeted drug neutralization, stimulus-responsive engineered protein prodrugs, and emerging multicomponent smart drug systems (e.g., antibody-drug conjugates, responsive engineered zymogens, prospective biochemical logic smart drug systems, drug buffers, and network medicine applications).

(Dis)similarity Index To Compare Correlated Motions in Molecular Simulations

Molecular dynamics (MD) simulations are widely used to complement or guide experimental studies in the characterization of protein dynamics, thanks to improvements in force-field accuracy, along with in the software and hardware to sample the conformational landscape of proteins. Among the different applications of MD simulations, the study of correlated motions is largely employed for different purposes. Several metrics have been developed to describe correlated motions in the MD ensemble, such as methods based on Pearson Correlation or Mutual Information. Cross-correlation analysis of MD trajectories is indeed appealing not only to identify residues characterized by coupled fluctuations in protein structures but also since it can be used to extrapolate motions along directions in which major conformational changes should occur, for example on longer time scales than the ones that are actually simulated. Nevertheless, most of the MD studies employ average correlation maps and mostly in a qualitative way, even when different systems or different replicates of the same system are compared. The broad application of correlation metrics in the analysis of MD simulations, especially for comparative purposes, requires a step forward toward more quantitative and accurate comparisons. We thus here employed a simple but effective index, which is based on a normalized Frobenius norm of the differences between protein correlation maps, to compare correlated motions. We applied this index for a quantitative comparison of correlated motions from MD simulations of seven proteins of different size and fold. We also employed the index to assess the robustness of correlation description when multi-replicate MD simulations of a same system are used, and we compared our index to metrics for comparison of structural ensembles such as Root Mean Square Inner Product and the Bhattacharyya Coefficient.

The taming of oxygen: biocatalytic oxyfunctionalisations

The scope and limitations of oxygenases as catalysts for preparative organic synthesis is discussed.

Directed evolution of phenylacetone monooxygenase as an active catalyst for the Baeyer-Villiger conversion of cyclohexanone to caprolactone

Phenylacetone monooxygenase (PAMO) is an exceptionally robust Baeyer-Villiger monooxygenase, which makes it ideal for potential industrial applications. However, its substrate scope is limited, unreactive cyclohexanone being a prominent example. Such a limitation is unfortunate, because this particular transformation in an ecologically viable manner would be highly desirable, the lactone and the respective lactam being of considerable interest as monomers in polymer science. We have applied directed evolution in search of an active mutant for this valuable C-C activating reaction. Using iterative saturation mutagenesis (ISM), several active mutants were evolved, with only a minimal trade-off in terms of stability. The best mutants allow for quantitative conversion of 2 mM cyclohexanone within 1 h reaction time. In order to circumvent the NADP(+) regeneration problem, whole E. coli resting cells were successfully applied. Molecular dynamics simulations and induced fit docking throw light on the origin of enhanced PAMO activity. The PAMO mutants constitute ideal starting points for future directed evolution optimization necessary for an industrial process.

Design of allosterically regulated protein catalysts

Activity of allosteric protein catalysts is regulated by an external stimulus, such as protein or small molecule binding, light activation, pH change, etc., at a location away from the active site of the enzyme. Since its original introduction in 1961, the concept of allosteric regulation has undergone substantial expansion, and many, if not most, enzymes have been shown to possess some degree of allosteric regulation. The ability to create new catalysts that can be turned on and off using allosteric interactions would greatly expand the chemical biology toolbox and will allow for detection of environmental pollutants and disease biomarkers and facilitate studies of cellular processes and metal homeostasis. Thus, design of allosterically regulated protein catalysts represents an actively growing area of research. In this paper, we describe various approaches to achieving regulation of catalysis.