Learning from directed evolution: Further lessons from theoretical investigations into cooperative mutations in lipase enantioselectivity

Learning from Protein Engineering by Deconvolution of Multi-Mutational Variants

This review analyzes a development in biochemistry, enzymology and biotechnology that originally came as a surprise. Following the establishment of directed evolution of stereoselective enzymes in organic chemistry, the concept of partial or complete deconvolution of selective multi-mutational variants was introduced. Early deconvolution experiments of stereoselective variants led to the finding that mutations can interact cooperatively or antagonistically with one another, not just additively. During the past decade, this phenomenon was shown to be general. In some studies, molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) computations were performed in order to shed light on the origin of non-additivity at all stages of an evolutionary upward climb. Data of complete deconvolution can be used to construct unique multi-dimensional rugged fitness pathway landscapes, which provide mechanistic insights different from traditional fitness landscapes. Along a related line, biochemists have long tested the result of introducing two point mutations in an enzyme for mechanistic reasons, followed by a comparison of the respective double mutant in so-called double mutant cycles, which originally showed only additive effects, but more recently also uncovered cooperative and antagonistic non-additive effects. We conclude with suggestions for future work, and call for a unified overall picture of non-additivity and epistasis.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

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.

Artificial cysteine-lipases with high activity and altered catalytic mechanism created by laboratory evolution

Engineering artificial enzymes with high activity and catalytic mechanism different from naturally occurring enzymes is a challenge in protein design. For example, many attempts have been made to obtain active hydrolases by introducing a Ser → Cys exchange at the respective catalytic triads, but this generally induced a breakdown of activity. We now report that this long-standing dogma no longer pertains, provided additional mutations are introduced by directed evolution. By employing Candida antarctica lipase B (CALB) as the model enzyme with the Ser-His-Asp catalytic triad, a highly active cysteine-lipase having a Cys-His-Asp catalytic triad and additional mutations W104V/A281Y/A282Y/V149G can be evolved, showing a 40-fold higher catalytic efficiency than wild-type CALB in the hydrolysis of 4-nitrophenyl benzoate, and tolerating bulky substrates. Crystal structures, kinetics, MD simulations and QM/MM calculations reveal dynamic features and explain all results, including the preference of a two-step mechanism involving the zwitterionic pair Cys105⁻/His224⁺ rather than a concerted process.

Candida antarctica lipase B (CALB) is a serine lipase. Here, the authors use directed evolution to exchange serine with cysteine in the catalytic triad of the enzyme, thereby obtaining a highly active CALB variant that — unlike the wild type — accommodates bulky substrates.

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.

QM/MM free energy simulations: recent progress and challenges

Due to the higher computational cost relative to pure molecular mechanical (MM) simulations, hybrid quantum mechanical/molecular mechanical (QM/MM) free energy simulations particularly require a careful consideration of balancing computational cost and accuracy. Here we review several recent developments in free energy methods most relevant to QM/MM simulations and discuss several topics motivated by these developments using simple but informative examples that involve processes in water. For chemical reactions, we highlight the value of invoking enhanced sampling technique (e.g., replica-exchange) in umbrella sampling calculations and the value of including collective environmental variables (e.g., hydration level) in metadynamics simulations; we also illustrate the sensitivity of string calculations, especially free energy along the path, to various parameters in the computation. Alchemical free energy simulations with a specific thermodynamic cycle are used to probe the effect of including the first solvation shell into the QM region when computing solvation free energies. For cases where high-level QM/MM potential functions are needed, we analyze two different approaches: the QM/MM-MFEP method of Yang and co-workers and perturbative correction to low-level QM/MM free energy results. For the examples analyzed here, both approaches seem productive although care needs to be exercised when analyzing the perturbative corrections.

Using experimental evolution to probe molecular mechanisms of protein function

Directed evolution is a powerful tool for engineering protein function. The process of directed evolution involves iterative rounds of sequence diversification followed by assaying activity of variants and selection. The range of sequence variants and linked activities generated in the course of an evolution are a rich information source for investigating relationships between sequence and function. Key residue positions determining protein function, combinatorial contributors to activity and even potential functional mechanisms have been revealed in directed evolutions. The recent application of high throughput sequencing substantially increases the information that can be retrieved from directed evolution experiments. Combined with computational analysis this additional sequence information has allowed high-resolution analysis of individual residue contributions to activity. These developments promise to significantly enhance the depth of insight that experimental evolution provides into mechanisms of protein function.

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.

Determination of lipolytic enzyme activities

Pseudomonas aeruginosa is a versatile human opportunistic pathogen that produces and secretes an arsenal of enzymes, proteins and small molecules many of which serve as virulence factors. Notably, about 40 % of P. aeruginosa genes code for proteins of unknown function, among them more than 80 encoding putative, but still unknown lipolytic enzymes. This group of hydrolases (EC 3.1.1) is known already for decades, but only recently, several of these enzymes have attracted attention as potential virulence factors. Reliable and reproducible enzymatic activity assays are crucial to determine their physiological function and particularly assess their contribution to pathogenicity. As a consequence of the unique biochemical properties of lipids resulting in the formation of micellar structures in water, the reproducible preparation of substrate emulsions is strongly dependent on the method used. Furthermore, the physicochemical properties of the respective substrate emulsion may drastically affect the activities of the tested lipolytic enzymes. Here, we describe common methods for the activity determination of lipase, esterase, phospholipase, and lysophospholipase. These methods cover lipolytic activity assays carried out in vitro, with cell extracts or separated subcellular compartments and with purified enzymes. We have attempted to describe standardized protocols, allowing the determination and comparison of enzymatic activities of lipolytic enzymes from different sources. These methods should also encourage the Pseudomonas community to address the wealth of still unexplored lipolytic enzymes encoded and produced by P. aeruginosa.

Biocatalysis in organic chemistry and biotechnology: past, present, and future

Enzymes as catalysts in synthetic organic chemistry gained importance in the latter half of the 20th century, but nevertheless suffered from two major limitations. First, many enzymes were not accessible in large enough quantities for practical applications. The advent of recombinant DNA technology changed this dramatically in the late 1970s. Second, many enzymes showed a narrow substrate scope, often poor stereo- and/or regioselectivity and/or insufficient stability under operating conditions. With the development of directed evolution beginning in the 1990s and continuing to the present day, all of these problems can be addressed and generally solved. The present Perspective focuses on these and other developments which have popularized enzymes as part of the toolkit of synthetic organic chemists and biotechnologists. Included is a discussion of the scope and limitation of cascade reactions using enzyme mixtures in vitro and of metabolic engineering of pathways in cells as factories for the production of simple compounds such as biofuels and complex natural products. Future trends and problems are also highlighted, as is the discussion concerning biocatalysis versus nonbiological catalysis in synthetic organic chemistry. This Perspective does not constitute a comprehensive review, and therefore the author apologizes to those researchers whose work is not specifically treated here.

Biocatalytic and structural properties of a highly engineered halohydrin dehalogenase

Two highly engineered halohydrin dehalogenase variants were characterized in terms of their performance in dehalogenation and epoxide cyanolysis reactions. Both enzyme variants outperformed the wild-type enzyme in the cyanolysis of ethyl (S)-3,4-epoxybutyrate, a conversion yielding ethyl (R)-4-cyano-3-hydroxybutyrate, an important chiral building block for statin synthesis. One of the enzyme variants, HheC2360, displayed catalytic rates for this cyanolysis reaction enhanced up to tenfold. Furthermore, the enantioselectivity of this variant was the opposite of that of the wild-type enzyme, both for dehalogenation and for cyanolysis reactions. The 37-fold mutant HheC2360 showed an increase in thermal stability of 8 °C relative to the wild-type enzyme. Crystal structures of this enzyme were elucidated with chloride and ethyl (S)-3,4-epoxybutyrate or with ethyl (R)-4-cyano-3-hydroxybutyrate bound in the active site. The observed increase in temperature stability was explained in terms of a substantial increase in buried surface area relative to the wild-type HheC, together with enhanced interfacial interactions between the subunits that form the tetramer. The structures also revealed that the substrate binding pocket was modified both by substitutions and by backbone movements in loops surrounding the active site. The observed changes in the mutant structures are partly governed by coupled mutations, some of which are necessary to remove steric clashes or to allow backbone movements to occur. The importance of interactions between substitutions suggests that efficient directed evolution strategies should allow for compensating and synergistic mutations during library design.

Computational enzyme design

Recent developments in computational chemistry and biology have come together in the "inside-out" approach to enzyme engineering. Proteins have been designed to catalyze reactions not previously accelerated in nature. Some of these proteins fold and act as catalysts, but the success rate is still low. The achievements and limitations of the current technology are highlighted and contrasted to other protein engineering techniques. On its own, computational "inside-out" design can lead to the production of catalytically active and selective proteins, but their kinetic performances fall short of natural enzymes. When combined with directed evolution, molecular dynamics simulations, and crowd-sourced structure-prediction approaches, however, computational designs can be significantly improved in terms of binding, turnover, and thermal stability.

Iterative combinatorial mutagenesis as an effective strategy for generation of deacetoxycephalosporin C synthase with improved activity toward penicillin G

An iterative combinatorial mutagenesis (ICM) strategy was used to engineer deacetoxycephalosporin C synthase of Streptomyces clavuligerus (scDAOCS) for improved activity toward penicillin G. Seven mutational sites were repeatedly combined onto a starter mutant (C155Y Y184H V275I C281Y) of scDAOCS. Eleven improved combinatorial mutants were identified from 24 mutants in four rounds of ICM.

Why is chemical synthesis and property optimization easier than expected?

Identifying optimal conditions for chemical and material synthesis as well as optimizing the properties of the products is often much easier than simple reasoning would predict. The potential search space is infinite in principle and enormous in practice, yet optimal molecules, materials, and synthesis conditions for many objectives can often be found by performing a reasonable number of distinct experiments. Considering the goal of chemical synthesis or property identification as optimal control problems provides insight into this good fortune. Both of these goals may be described by a fitness function J that depends on a suitable set of variables (e.g., reactant concentrations, components of a material, processing conditions, etc.). The relationship between J and the variables specifies the fitness landscape for the target objective. Upon making simple physical assumptions, this work demonstrates that the fitness landscape for chemical optimization contains no local sub-optimal maxima that may hinder attainment of the absolute best value of J. This feature provides a basis to explain the many reported efficient optimizations of synthesis conditions and molecular or material properties. We refer to this development as OptiChem theory. The predicted characteristics of chemical fitness landscapes are assessed through a broad examination of the recent literature, which shows ample evidence of trap-free landscapes for many objectives. The fundamental and practical implications of OptiChem theory for chemistry are discussed.

Iterative saturation mutagenesis accelerates laboratory evolution of enzyme stereoselectivity: rigorous comparison with traditional methods

Efficacy in laboratory evolution of enzymes is currently a pressing issue, making comparative studies of different methods and strategies mandatory. Recent reports indicate that iterative saturation mutagenesis (ISM) provides a means to accelerate directed evolution of stereoselectivity and thermostability, but statistically meaningful comparisons with other methods have not been documented to date. In the present study, the efficacy of ISM has been rigorously tested by applying it to the previously most systematically studied enzyme in directed evolution, the lipase from Pseudomonas aeruginosa as a catalyst in the stereoselective hydrolytic kinetic resolution of a chiral ester. Upon screening only 10,000 transformants, unprecedented enantioselectivity was achieved (E = 594). ISM proves to be considerably more efficient than all previous systematic efforts utilizing error-prone polymerase chain reaction at different mutation rates, saturation mutagenesis at hot spots, and/or DNA shuffling, pronounced positive epistatic effects being the underlying reason.