Protein engineering using mutability landscapes: Controlling site-selectivity of P450-catalyzed steroid hydroxylation

Directed evolution and rational design have been used widely in engineering enzymes for their application in synthetic organic chemistry and biotechnology. With stereoselectivity playing a crucial role in catalysis for the synthesis of valuable chemical and pharmaceutical compounds, rational design has not achieved such wide success in this specific area compared to directed evolution. Nevertheless, one bottleneck of directed evolution is the laborious screening efforts and the observed trade-offs in catalytic profiles. This has motivated researchers to develop more efficient protein engineering methods. As a prime approach, mutability landscaping avoids such trade-offs by providing more information of sequence-function relationships. Here, we describe an application of this efficient protein engineering method to improve the regio-/stereoselectivity and activity of P450BM3 for steroid hydroxylation, while keeping the mutagenesis libraries small so that they will require only minimal screening.

Rational enzyme design for enabling biocatalytic Baldwin cyclization and asymmetric synthesis of chiral heterocycles

Chiral heterocyclic compounds are needed for important medicinal applications. We report an in silico strategy for the biocatalytic synthesis of chiral N- and O-heterocycles via Baldwin cyclization modes of hydroxy- and amino-substituted epoxides and oxetanes using the limonene epoxide hydrolase from Rhodococcus erythropolis. This enzyme normally catalyzes hydrolysis with formation of vicinal diols. Firstly, the required shutdown of the undesired natural water-mediated ring-opening is achieved by rational mutagenesis of the active site. In silico enzyme design is then continued with generation of the improved mutants. These variants prove to be versatile catalysts for preparing chiral N- and O-heterocycles with up to 99% conversion, and enantiomeric ratios up to 99:1. Crystal structural data and computational modeling reveal that Baldwin-type cyclizations, catalyzed by the reprogrammed enzyme, are enabled by reshaping the active-site environment that directs the distal RHN and HO-substituents to be intramolecular nucleophiles.

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.

The Unexplored Importance of Fleeting Chiral Intermediates in Enzyme-Catalyzed Reactions

Decades of extensive research efforts by biochemists, organic chemists, and protein engineers have led to an understanding of the basic mechanisms of essentially all known types of enzymes, but in a formidable number of cases an essential aspect has been overlooked. The occurrence of short-lived chiral intermediates formed by symmetry-breaking of prochiral precursors in enzyme catalyzed reactions has been systematically neglected. We designate these elusive species as fleeting chiral intermediates and analyze such crucial questions as "Do such intermediates occur in homochiral form?" If so, what is the absolute configuration, and why did Nature choose that particular stereoisomeric form, even when the isolable final product may be achiral? Does the absolute configuration of a chiral product depend in any way on the absolute configuration of the fleeting chiral precursor? How does this affect the catalytic proficiency of the enzyme? If these issues continue to be unexplored, then an understanding of the mechanisms of many enzyme types remains incomplete. We have systematized the occurrence of these chiral intermediates according to their structures and enzyme types. This is followed by critical analyses of selected case studies and by final conclusions and perspectives. We hope that the fascinating concept of fleeting chiral intermediates will attract the attention of scientists, thereby opening an exciting new research field.

Machine Learning Enables Selection of Epistatic Enzyme Mutants for Stability Against Unfolding and Detrimental Aggregation

Machine learning (ML) has pervaded most areas of protein engineering, including stability and stereoselectivity. Using limonene epoxide hydrolase as the model enzyme and innov'SAR as the ML platform, comprising a digital signal process, we achieved high protein robustness that can resist unfolding with concomitant detrimental aggregation. Fourier transform (FT) allows us to take into account the order of the protein sequence and the nonlinear interactions between positions, and thus to grasp epistatic phenomena. The innov'SAR approach is interpolative, extrapolative and makes outside-the-box, predictions not found in other state-of-the-art ML or deep learning approaches. Equally significant is the finding that our approach to ML in the present context, flanked by advanced molecular dynamics simulations, uncovers the connection between epistatic mutational interactions and protein robustness.

Towards self-doping multimetal porphyrin systems

An approach for the possible production of novel bimetallic self-doped porphyrin-based compounds of potential interest in material science is reported. Heating Cu(II)tetraphenylporphyrin (TPPCu) with chromocene at 120°C in benzonitrile affords the crystalline multimetal porphyrin system TPPCu/TPPCr in good yield. The X-ray single crystal structural analysis reveals a random distribution of TPPCu and TPPCr, with a Cu:Cr ratio of 71(2):29(2)%. Exploratory DFT calculations of TPPCu/TPPCr indicate little if any electron transfer. In contrast, calculations of a hypothetical cationic TPPCu/TPPRu system indicates the possibility of self-doping.

Regio- and Stereoselective Steroid Hydroxylation at C7 by Cytochrome P450 Monooxygenase Mutants

Steroidal C7β alcohols and their respective esters have shown significant promise as neuroprotective and anti-inflammatory agents to treat chronic neuronal damage like stroke, brain trauma, and cerebral ischemia. Since C7 is spatially far away from any functional groups that could direct C-H activation, these transformations are not readily accessible using modern synthetic organic techniques. Reported here are P450-BM3 mutants that catalyze the oxidative hydroxylation of six different steroids with pronounced C7 regioselectivities and β stereoselectivities, as well as high activities. These challenging transformations were achieved by a focused mutagenesis strategy and application of a novel technology for protein library construction based on DNA assembly and USER (Uracil-Specific Excision Reagent) cloning. Upscaling reactions enabled the purification of the respective steroidal alcohols in moderate to excellent yields. The high-resolution X-ray structure and molecular dynamics simulations of the best mutant unveil the origin of regio- and stereoselectivity.

P450-BM3-Catalyzed Sulfoxidation versus Hydroxylation: A Common or Two Different Catalytically Active Species?

While the mechanism of the P450-catalyzed oxidative hydroxylation of organic compounds has been studied in detail for many years, less is known about sulfoxidation. Depending upon the structure of the respective substrate, heme-Fe=O (Cpd I), heme–Fe(III)–OOH (Cpd 0), and heme–Fe(III)–H₂O₂ (protonated Cpd 0) have been proposed as reactive intermediates. In the present study, we consider the transformation of isosteric substrates via sulfoxidation and oxidative hydroxylation, respectively, catalyzed by regio- and enantioselective mutants of P450-BM3 which were constructed by directed evolution. 1-Thiochromanone and 1-tetralone were used as the isosteric substrates because, unlike previous studies involving fully flexible compounds such as thia-fatty acids and fatty acids, respectively, these compounds are rigid and cannot occur in a multitude of different conformations and binding modes in the large P450-BM3 binding pocket. The experimental results comprising activity and regio- and enantioselectivity, flanked by molecular dynamics computations within a time scale of 300 ns and QM/MM calculations of transition-state energies, unequivocally show that heme-Fe=O (Cpd I) is the common catalytically active intermediate in both sulfoxidation and oxidative hydroxylation.

Exploiting Designed Oxidase-Peroxygenase Mutual Benefit System for Asymmetric Cascade Reactions

A unique P450 monooxygenase–peroxygenase mutual benefit system was designed as the core element in the construction of a biocatalytic cascade reaction sequence leading from 3-phenyl propionic acid to (R)-phenyl glycol. In this system, P450 monooxygenase (P450-BM3) and P450 peroxygenase (OleT_JE) not only function as catalysts for the crucial initial reactions, they also ensure an internal in situ H₂O₂ recycle mechanism that avoids its accumulation and thus prevents possible toxic effects. By directed evolution of P450-BM3 as the catalyst in the enantioselective epoxidation of the styrene-intermediate, formed from 3-phenyl propionic acid, and the epoxide hydrolase ANEH for final hydrolytic ring opening, (R)-phenyl glycol and 9 derivatives thereof were synthesized from the respective carboxylic acids in one-pot processes with high enantioselectivity.

Directed Evolution of Artificial Metalloenzymes: A Universal Means to Tune the Selectivity of Transition Metal Catalysts?

Transition metal catalysts mediate a wide variety of chemo-, stereo-, and regioselective transformations, and therefore play a pivotal role in modern synthetic organic chemistry. Steric and electronic effects of ligands provide organic chemists with an exceedingly useful tool. More than four decades ago, chemists began to think about a different approach, namely, embedding achiral ligand/metal moieties covalently or noncovalently in protein hosts with formation of artificial metalloenzymes. While structurally fascinating, this approach led in each case only to a single (bio)catalyst, with its selectivity and activity being a matter of chance. In order to solve this fundamental problem, my group proposed in 2000-2002 the idea of directed evolution of artificial metalloenzymes. In earlier studies, we had already demonstrated that directed evolution of enzymes constitutes a viable method for enhancing and inverting the stereoselectivity of enzymes as catalysts in organic chemistry. We speculated that it should also be possible to manipulate selectivity and activity of artificial metalloenzymes, which would provide organic chemists with a tool for optimizing essentially any transition metal catalyzed reaction type. In order to put this vision into practice, we first turned to the Whitesides system for artificial metalloenzyme formation, comprising a biotinylated diphosphine/Rh moiety, which is anchored noncovalently to avidin or streptavidin. Following intensive optimization, proof of principle was finally demonstrated in 2006, which opened the door to a new research area. This personal Account critically assesses these early studies as well as subsequent efforts from my group focusing on different protein scaffolds, and includes briefly some of the most important current contributions of other groups. Two primary messages emerge: First, since organic chemists continue to be extremely good at designing and implementing man-made transition metal catalysts, often on a large scale, those scientists that are active in the equally intriguing field of directed evolution of artificial metalloenzymes should be moderate when generalizing claims. All factors required for a truly viable catalytic system need to be considered, especially activity and ease of upscaling. Second, the most exciting and thus far very rare cases of directed evolution of artificial metalloenzymes are those that focus on selective transformations that are not readily possible using state of the art transition metal catalysts.

Chemo- and Regioselective Dihydroxylation of Benzene to Hydroquinone Enabled by Engineered Cytochrome P450 Monooxygenase

Hydroquinone (HQ) is produced commercially from benzene by multi-step Hock-type processes with equivalent amounts of acetone as side-product. We describe an efficient biocatalytic alternative using the cytochrome P450-BM3 monooxygenase. Since the wildtype enzyme does not accept benzene, a semi-rational protein engineering strategy was developed. Highly active mutants were obtained which transform benzene in a one-pot sequence first into phenol and then regioselectively into HQ without any overoxidation. A computational study shows that the chemoselective oxidation of phenol by the P450-BM3 variant A82F/A328F leads to the regioselective formation of an epoxide intermediate at the C3=C4 double bond, which departs from the binding pocket and then undergoes fragmentation in aqueous medium with exclusive formation of HQ. As a practical application, an E. coli designer cell system was constructed, which enables the cascade transformation of benzene into the natural product arbutin, which has anti-inflammatory and anti-bacterial activities.

A machine learning approach for reliable prediction of amino acid interactions and its application in the directed evolution of enantioselective enzymes

Directed evolution is an important research activity in synthetic biology and biotechnology. Numerous reports describe the application of tedious mutation/screening cycles for the improvement of proteins. Recently, knowledge-based approaches have facilitated the prediction of protein properties and the identification of improved mutants. However, epistatic phenomena constitute an obstacle which can impair the predictions in protein engineering. We present an innovative sequence-activity relationship (innov'SAR) methodology based on digital signal processing combining wet-lab experimentation and computational protein design. In our machine learning approach, a predictive model is developed to find the resulting property of the protein when the n single point mutations are permuted (2n combinations). The originality of our approach is that only sequence information and the fitness of mutants measured in the wet-lab are needed to build models. We illustrate the application of the approach in the case of improving the enantioselectivity of an epoxide hydrolase from Aspergillus niger. n = 9 single point mutants of the enzyme were experimentally assessed for their enantioselectivity and used as a learning dataset to build a model. Based on combinations of the 9 single point mutations (29), the enantioselectivity of these 512 variants were predicted, and candidates were experimentally checked: better mutants with higher enantioselectivity were indeed found.

Manipulating the stereoselectivity of a thermostable alcohol dehydrogenase by directed evolution for efficient asymmetric synthesis of arylpropanols

Chiral arylpropanols are valuable components in important pharmaceuticals and fragrances, which is the motivation for previous attempts to prepare these building blocks enantioselectively in asymmetric processes using either enzymes or transition metal catalysts. Thus far, enzymes used in kinetic resolution proved to be best, but several problems prevented ecologically and economically viable processes from being developed. In the present study, directed evolution was applied to the thermostable alcohol dehydrogenase TbSADH in the successful quest to obtain mutants that are effective in the dynamic reductive kinetic resolution (DYRKR) of racemic arylpropanals. Using rac-2-phenyl-1-propanal in a model reaction, (S)- and (R)-selective mutants were evolved which catalyzed DYRKR of this racemic substrate with formation of the respective (S)- and (R)-alcohols in essentially enantiomerically pure form. This was achieved on the basis of an unconventional form of iterative saturation mutagenesis (ISM) at randomization sites lining the binding pocket using a reduced amino acid alphabet. The best mutants were also effective in the DYRKR of several other structurally related racemic aldehydes.

Rapid and Error-Free Site-Directed Mutagenesis by a PCR-Free In Vitro CRISPR/Cas9-Mediated Mutagenic System

The quality and efficiency of any PCR-based mutagenesis technique may not be optimal due to, among other things, amino acid bias, which means that the development of efficient PCR-free methods is desirable. Here, we present a highly efficient in vitro CRISPR/Cas9-mediated mutagenic (ICM) system that allows rapid construction of designed mutants in a PCR-free manner. First, it involves plasmid digestion by utilizing a complex of Cas9 with specific single guide RNA (sgRNA) followed by degradation with T5 exonuclease to generate a 15 nt homologous region. Second, primers containing the desired mutations are annealed to form the double-stranded DNA fragments, which are then ligated into the linearized plasmid. In theory, neither the size of the target plasmid nor the unavailable restriction enzyme site poses any problems that may arise in traditional techniques. In this study, single and multiple site-directed mutagenesis were successfully performed even for a large size plasmid (up to 9.0 kb). Moreover, a PCR-free site-saturation mutagenesis library on single site and two adjacent sites of a green fluorescent protein was also generated with promising results. This demonstrates the great potential of the ICM system for creating high-quality mutant libraries in directed evolution as an alternative to PCR-based saturation mutagenesis, thus facilitating research on synthetic biology.

Overriding Traditional Electronic Effects in Biocatalytic Baeyer-Villiger Reactions by Directed Evolution

Controlling the regioselectivity of Baeyer-Villiger (BV) reactions remains an ongoing issue in organic chemistry, be it by synthetic catalysts or enzymes of the type Baeyer-Villiger monooxygenases (BVMOs). Herein, we address the challenging problem of switching normal to abnormal BVMO regioselectivity by directed evolution using three linear ketones as substrates, which are not structurally biased toward abnormal reactivity. Upon applying iterative saturation mutagenesis at sites lining the binding pocket of the thermostable BVMO from Thermocrispum municipale DSM 44069 (TmCHMO) and using 4-phenyl-2-butanone as substrate, the regioselectivity was reversed from 99:1 (wild-type enzyme in favor of the normal product undergoing 2-phenylethyl migration) to 2:98 in favor of methyl migration when applying the best mutant. This also stands in stark contrast to the respective reaction using the synthetic reagent m-CPBA, which provides solely the normal product. Reversal of regioselectivity was also achieved in the BV reaction of two other linear ketones. Kinetic parameters and melting temperatures revealed that most of the evolved mutants retained catalytic activity, as well as thermostability. In order to shed light on the origin of switched regioselectivity in reactions of 4-phenyl-2-butanone and phenylacetone, extensive QM/MM and MD simulations were performed. It was found that the mutations introduced by directed evolution induce crucial changes in the conformation of the respective Criegee intermediates and transition states in the binding pocket of the enzyme. In mutants that destabilize the normally preferred migration transition state, a reversal of regioselectivity is observed. This conformational control of regioselectivity overrides electronic control, which normally causes preferential migration of the group that is best able to stabilize positive charge. The results can be expected to aid future protein engineering of BVMOs.

1-Butanol as a Solvent for Efficient Extraction of Polar Compounds from Aqueous Medium: Theoretical and Practical Aspects

The extraction of polar molecules from aqueous solution is a challenging task in organic synthesis. 1-Butanol has been used sporadically as an eluent for polar molecules, but it is unclear which molecular features drive its efficiency. Here, we employ free energy simulations to study the partitioning of 15 solutes between water and 1-butanol. The simulations demonstrate that the high affinity of polar molecules to the wet 1-butanol phase is associated with its nanostructure. Small inverse micelles of water are able to accommodate polar solutes and locally mimic an aqueous environment. We verify the simulations based on partition coefficients between water and 1-octanol, and include a blind prediction of the water/1-butanol partition coefficient of cyclohexane-1,2-diol. The calculations are in excellent agreement with experiment, reaching root-mean-square deviations below 0.7 kcal/mol. Actual extractions of cyclohexane-1,2-diol from buffer solutions that mimic cell lysates and suspensions in biocatalytic reactions further exemplify our findings. The yields highlight that extractions with 1-butanol can be significantly more efficient than the conventional protocol based on ethyl acetate.

Boosting the efficiency of site-saturation mutagenesis for a difficult-to-randomize gene by a two-step PCR strategy

Site-saturation mutagenesis (SSM) has been used in directed evolution of proteins for a long time. As a special form of saturation mutagenesis, it involves individual randomization at a given residue with formation of all 19 amino acids. To date, the most efficient embodiment of SSM is a one-step PCR-based approach using NNK codon degeneracy. However, in the case of difficult-to-randomize genes, SSM may not deliver all of the expected 19 mutants, which compels the user to invest further efforts by applying site-directed mutagenesis for the construction of the missing mutants. To solve this problem, we developed a two-step PCR-based technique in which a mutagenic primer and a non-mutagenic (silent) primer are used to generate a short DNA fragment, which is recovered and then employed as a megaprimer to amplify the whole plasmid. The present two-step and older one-step (partially overlapped primer approach) procedures were compared by utilizing cytochrome P450-BM3, which is a "difficult-to-randomize" gene. The results document the distinct superiority of the new method by checking the library quality on DNA level based on massive sequence data, but also at amino acid level. Various future applications in biotechnology can be expected, including the utilization when constructing mutability landscapes, which provide semi-rational information for identifying hot spots for protein engineering and directed evolution.

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.

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.

Correction: Simultaneous engineering of an enzyme's entrance tunnel and active site: the case of monoamine oxidase MAO-N

Correction for ‘Simultaneous engineering of an enzyme’s entrance tunnel and active site: the case of monoamine oxidase MAO-N’ by Guangyue Li et al., Chem. Sci., 2017, 8, 4093–4099.

Simultaneous engineering of an enzyme's entrance tunnel and active site: the case of monoamine oxidase MAO-N

A new directed evolution approach is presented to enhance the activity of an enzyme and to manipulate stereoselectivity by focusing iterative saturation mutagenesis (ISM) simultaneously on residues lining the entrance tunnel and the binding pocket. This combined mutagenesis strategy was applied successfully to the monoamine oxidase from Aspergillus niger (MAO-N) in the reaction of sterically demanding substrates which are of interest in the synthesis of chiral pharmaceuticals based on the benzo-piperidine scaffold. Reversal of enantioselectivity of Turner-type deracemization was achieved in the synthesis of (S)-1,2,3,4-tetrahydro-1-methyl-isoquinoline, (S)-1,2,3,4-tetrahydro-1-ethylisoquinoline and (S)-1,2,3,4-tetrahydro-1-isopropylisoquinoline. Extensive molecular dynamics simulations indicate that the altered catalytic profile is due to increased hydrophobicity of the entrance tunnel acting in concert with the altered shape of the binding pocket.