Combination of computational prescreening and experimental library construction can accelerate enzyme optimization by directed evolution

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.

Combinatorial Approach for Exploring Conformational Space and Activation Barriers in Computer-Aided Enzyme Design

Computer-aided enzyme design is a field of great potential importance for biotechnological applications, medical advances, and a fundamental understanding of enzyme action. However, reaching a predictive ability in this direction is extremely challenging. It requires both the ability to predict quantitatively the activation barriers in cases where the structure and sequence are known and the ability to predict the effect of different mutations. In this work, we propose a protocol for predicting reasonable starting structures of mutants of proteins with known structures and for calculating the activation barriers of the generated mutants. Our approach also allows us to use the predicted structures of the generated mutant to predict structures and activation barriers for subsequent set of mutations. This protocol is used to examine the reliability of the in silico directed evolution of Kemp eliminase and haloalkane dehalogenase. We also used the results of single and double mutations as a base for predicting the effect of transition-state stabilization by multiple concurrent mutations. This strategy seems to be useful in creating an activity funnel that provides a qualitative ranking of the catalytic power of different mutants.

Computational Identification of Amino-Acid Mutations that Further Improve the Activity of a Chalcone-Flavonone Isomerase from Glycine max

Protein design for improving enzymatic activity remains a challenge in biochemistry, especially to identify target amino-acid sites for mutagenesis and to design beneficial mutations for those sites. Here, we employ a computational approach that combines multiple sequence alignment, positive selection detection, and molecular docking to identify and design beneficial amino-acid mutations that further improve the intramolecular-cyclization activity of a chalcone-flavonone isomerase from Glycine max (GmCHI). By this approach, two GmCHI mutants with higher activities were predicted and verified. The results demonstrate that this approach could determine the beneficial amino-acid mutations for improving the enzymatic activity, and may find more applications in engineering of enzymes.

Valence bond and enzyme catalysis: a time to break down and a time to build up

Understanding enzyme catalysis and developing ability to control of it are two great challenges in biochemistry. A few successful examples of computational-based enzyme design have proved the fantastic potential of computational approaches in this field, however, relatively modest rate enhancements have been reported and the further development of complementary methods is still required. Herein we propose a conceptually simple scheme to identify the specific role that each residue plays in catalysis. The scheme is based on a breakdown of the total catalytic effect into contributions of individual protein residues, which are further decomposed into chemically interpretable components by using valence bond theory. The scheme is shown to shed light on the origin of catalysis in wild-type haloalkane dehalogenase (wt-DhlA) and its mutants. Furthermore, the understanding gained through our scheme is shown to have great potential in facilitating the selection of non-optimal sites for catalysis and suggesting effective mutations to enhance the enzymatic rate.

Enhancing the efficiency of directed evolution in focused enzyme libraries by the adaptive substituent reordering algorithm

Directed evolution is a broadly successful strategy for protein engineering in the quest to enhance the stereoselectivity, activity, and thermostability of enzymes. To increase the efficiency of directed evolution based on iterative saturation mutagenesis, the adaptive substituent reordering algorithm (ASRA) is introduced here as an alternative to traditional quantitative structure-activity relationship (QSAR) methods for identifying potential protein mutants with desired properties from minimal sampling of focused libraries. The operation of ASRA depends on identifying the underlying regularity of the protein property landscape, allowing it to make predictions without explicit knowledge of the structure-property relationships. In a proof-of-principle study, ASRA identified all or most of the best enantioselective mutants among the synthesized epoxide hydrolase from Aspergillus niger, in the absence of peptide seeds with high E-values. ASRA even revealed a laboratory error from irregularities of the reordered E-value landscape alone.

Rapid sequence scanning mutagenesis using in silico oligo design and the Megaprimer PCR of whole plasmid method (MegaWHOP)

A wide variety of random- and site-directed mutagenesis techniques have been developed to investigate the structure-function relationship in proteins and intergenic regions like promoter sequences. Similar techniques can be employed to optimize protein properties like enantioselectivity, substrate specificity, and stability in a directed evolution approach. Due to the tremendous genetic diversity that is created by common random-mutagenesis methods, directed evolution techniques usually require the time-consuming and cumbersome screening of large numbers of variants. A gene-scanning saturation-mutagenesis approach represents one efficient way to limit the screening effort by reducing the created genetic diversity. In structure/function studies often a similar method, e.g., alanine- or arginine-scanning mutagenesis, is used to probe the role of specific amino acids in a protein. Here, we present a standardized mutagenesis strategy that can speed up the process of scanning whole proteins for structure/function studies and, furthermore, allows for the fast and efficient generation of gene-scanning saturation-mutagenesis libraries to be used in the directed evolution of enzyme functions and properties. The described method uses automated computer-assisted oligonucleotide design, and a two-step PCR-mutagenesis protocol to amplify site-specifically mutated circular plasmids that can be directly transformed in Escherichia coli expression strains.

Combined use of experimental and computational screens to characterize protein stability

One of the primary goals of protein design is to engineer proteins with improved stability. Protein stability is a key issue for chemical, biotechnology and pharmaceutical industries. The development of robust proteins/enzymes with the ability to withstand the potentially harsh conditions of industrial operations is of high importance. A number of strategies are currently being employed to achieve this goal. Two particular approaches, (i) directed evolution and (ii) computational protein design, are quite powerful yet have only recently been combined or applied and analyzed in parallel. In directed evolution, libraries of variants are searched experimentally for clones possessing the desired properties. With computational methods, protein design algorithms are utilized to perform in silico screening for stable protein sequences. Here, we used gene libraries of an unstable variant of streptococcal protein G (Gbeta1) and an in vivo screening method to identify stabilized variants. Many variants with notably increased thermal stabilities were isolated and characterized. Concomitantly, computational techniques and protein design algorithms were used to perform in silico screening of the same destabilized variant of Gbeta1. The combined use, and critical analysis, of these methods promises to advance the field of protein design.

Inter-conversion of catalytic abilities in a bifunctional carboxyl/feruloyl-esterase from earthworm gut metagenome

Carboxyl esterases (CE) exhibit various reaction specificities despite of their overall structural similarity. In present study we have exploited functional metagenomics, saturation mutagenesis and experimental protein evolution to explore residues that have a significant role in substrate discrimination. We used an enzyme, designated 3A6, derived from the earthworm gut metagenome that exhibits CE and feruloyl esterase (FAE) activities with p‐nitrophenyl and cinnamate esters, respectively, with a [(k_cat/K_m)]_CE/[(k_cat/K_m)]_FAE factor of 17. Modelling‐guided saturation mutagenesis at specific hotspots (Lys²⁸¹, Asp²⁸², Asn³¹⁶ and Lys³¹⁷) situated close to the catalytic core (Ser¹⁴³/Asp²⁷³/His³⁰⁵) and a deletion of a 34‐AA–long peptide fragment yielded mutants with the highest CE activity, while cinnamate ester bond hydrolysis was effectively abolished. Although, single to triple mutants with both improved activities (up to 180‐fold in k_cat/K_m values) and enzymes with inverted specificity ((k_cat/K_m)_CE/(k_cat/K_m)_FAE ratio of ∼0.4) were identified, no CE inactive variant was found. Screening of a large error‐prone PCR‐generated library yielded by far less mutants for substrate discrimination. We also found that no significant changes in CE activation energy occurs after any mutation (7.3 to −5.6 J mol⁻¹), whereas a direct correlation between loss/gain of FAE function and activation energies (from 33.05 to −13.7 J mol⁻¹) was found. Results suggest that the FAE activity in 3A6 may have evolved via introduction of a limited number of ‘hot spot’ mutations in a common CE ancestor, which may retain the original hydrolytic activity due to lower restrictive energy barriers but conveys a dynamic energetically favourable switch of a second hydrolytic reaction.

Force-field parameters for the simulation of tetrahedral intermediates of serine hydrolases

CHARMM force-field parameters are reported for the tetrahedral intermediate of serine hydrolases. The fitting follows the standard protocol proposed for CHARMM22. The reference data include ab initio (RHF/6-31G*) interaction energies of complexes between water and the model compound 1,1-dimethoxyethoxide, torsional profiles of related model compounds from correlated ab initio (MP2/6-311+G*//B3LYP/6-31+G*) calculations, as well as molecular geometries and vibrational frequencies from density functional theory (B3LYP/6-31+G*). The optimized parameters reproduce the target data well. Their utility is demonstrated by a QM/MM study of the tetrahedral intermediate in Bacillus subtilis lipase A, and by classical molecular modeling of enantioselectivity in Pseudomonas aeruginosa lipase and its mutants.

Influence of delta-functional groups on the enantiorecognition of secondary alcohols by Candida antarctica lipase B

The selectivity of acetylation of delta-functionalized secondary alcohols catalyzed by Candida antarctica lipase B has been examined by molecular dynamics. The results from the simulation show that a delta-alcohol functionality forms a hydrogen bond with the carbonyl group of Thr 40. This interaction stabilizes the tetrahedral intermediate and thus leads to selective acetylation of the R enantiomer. A stabilizing interaction of the delta-(R)-acetoxy group with the peptide NH of alanine 282 was also observed. No stabilizing interaction could be found for the delta-keto functionality, and it is proposed that this is the reason for the experimentally observed decrease in enantioselectivity. From these results, it was hypothesized that the enantioselectivity could be restored by mutating the alanine in position 281 for serine. The mutation was made experimentally, and the results show that the E value increased from 9 to 120.

Better library design: data-driven protein engineering

Data-driven protein engineering is increasingly used as an alternative to rational design and combinatorial engineering because it uses available knowledge to limit library size, while still allowing for the identification of unpredictable substitutions that lead to large effects. Recent advances in computational modeling and bioinformatics, as well as an increasing databank of experiments on functional variants, have led to new strategies to choose particular amino acid residues to vary in order to increase the chances of obtaining a variant protein with the desired property. Strategies for limiting diversity at each position, design of small sub-libraries, and the performance of scouting experiments, have also been developed or even automated, further reducing the library size.

Recent progress in engineering alpha/beta hydrolase-fold family members

The members of the alpha/beta hydrolase-fold family represent a functionally versatile group of enzymes with many important applications in biocatalysis. Given the technical significance of alpha/beta hydrolases in processes ranging from the kinetic resolution of enantiomeric precursors for pharmaceutical compounds to bulk products such as laundry detergent, optimizing and tailoring enzymes for these applications presents an ongoing challenge to chemists, biochemists, and engineers alike. A review of the recent literature on alpha/beta hydrolase engineering suggests that the early successes of "random processes" such as directed evolution are now being slowly replaced by more hypothesis-driven, focused library approaches. These developments reflect a better understanding of the enzymes' structure-function relationship and improved computational resources, which allow for more sophisticated search and prediction algorithms, as well as, in a very practical sense, the realization that bigger is not always better.

Biocatalytic production of enantiopure cyclohexane-trans-1,2-diol using extracellular lipases from Bacillus subtilis

Two extracellular lipases from Bacillus subtilis, B. subtilis lipase A and lipase B, have been expressed in the heterologous host Escherichia coli, biochemically characterized and used for the kinetic resolution of (rac)-trans-1,2-diacetoxycyclohexane. Both enzymes were selectively acting on the (R,R)-enantiomer of the racemic substrate, highly specifically hydrolyzing only one of the two ester groups present, thus allowing the preparation of enantiopure (R,R)- and (S,S)-cyclohexane-trans-1,2-diol. The reaction conditions for the use of purified enzyme and crude cell lyophilizate were optimized and reactions in batch and repetitive batch modes were carried out on a preparative scale to yield enantiopure product (>99% enantiomeric excess).