Evolutionary history of a specialized p450 propane monooxygenase

Strategies found not to be suitable for stabilizing high steroid hydroxylation activities of CYP450 BM3-based whole-cell biocatalysts

The implementation of biocatalytic steroid hydroxylation processes plays a crucial role in the pharmaceutical industry due to a plethora of medicative effects of hydroxylated steroid derivatives and their crucial role in drug approval processes. Cytochrome P450 monooxygenases (CYP450s) typically constitute the key enzymes catalyzing these reactions, but commonly entail drawbacks such as poor catalytic rates and the dependency on additional redox proteins for electron transfer from NAD(P)H to the active site. Recently, these bottlenecks were overcome by equipping Escherichia coli cells with highly active variants of the self-sufficient single-component CYP450 BM3 together with hydrophobic outer membrane proteins facilitating cellular steroid uptake. The combination of the BM3 variant KSA14m and the outer membrane pore AlkL enabled exceptionally high testosterone hydroxylation rates of up to 45 U gCDW-1 for resting (i.e., living but non-growing) cells. However, a rapid loss of specific activity heavily compromised final product titers and overall space-time yields. In this study, several stabilization strategies were evaluated on enzyme-, cell-, and reaction level. However, neither changes in biocatalyst configuration nor variation of cultivation media, expression systems, or inducer concentrations led to considerable improvement. This qualified the so-far used genetic construct pETM11-ksa14m-alkL, M9 medium, and the resting-cell state as the best options enabling comparatively efficient activity along with fast growth prior to biotransformation. In summary, we report several approaches not enabling a stabilization of the high testosterone hydroxylation rates, providing vital guidance for researchers tackling similar CYP450 stability issues. A comparison with more stable natively steroid-hydroxylating CYP106A2 and CYP154C5 in equivalent setups further highlighted the high potential of the investigated CYP450 BM3-based whole-cell biocatalysts. The immense and continuously developing repertoire of enzyme engineering strategies provides promising options to stabilize the highly active biocatalysts.

Utilization of low-stability variants in protein evolutionary engineering

Evolutionary engineering involves repeated mutations and screening and is widely used to modify protein functions. However, it is important to diversify evolutionary pathways to eliminate the bias and limitations of the variants by using traditionally unselected variants. In this study, we focused on low-stability variants that are commonly excluded from evolutionary processes and tested a method that included an additional restabilization step. The esterase from the thermophilic bacterium Alicyclobacillus acidocaldarius was used as a model protein, and its activity at its optimum temperature of 65 °C was improved by evolutionary experiments using random mutations by error-prone PCR. After restabilization using low-stability variants with low-temperature (37 °C) activity, several re-stabilizing variants were obtained from a large number of variant libraries. Some of the restabilized variants achieved by removing the destabilizing mutations showed higher activity than that of the wild-type protein. This implies that low-stability variants with low-temperature activity can be re-evolved for future use. This method will enable further diversification of evolutionary pathways.

Robust genetic codes enhance protein evolvability

The standard genetic code defines the rules of translation for nearly every life form on Earth. It also determines the amino acid changes accessible via single-nucleotide mutations, thus influencing protein evolvability-the ability of mutation to bring forth adaptive variation in protein function. One of the most striking features of the standard genetic code is its robustness to mutation, yet it remains an open question whether such robustness facilitates or frustrates protein evolvability. To answer this question, we use data from massively parallel sequence-to-function assays to construct and analyze 6 empirical adaptive landscapes under hundreds of thousands of rewired genetic codes, including those of codon compression schemes relevant to protein engineering and synthetic biology. We find that robust genetic codes tend to enhance protein evolvability by rendering smooth adaptive landscapes with few peaks, which are readily accessible from throughout sequence space. However, the standard genetic code is rarely exceptional in this regard, because many alternative codes render smoother landscapes than the standard code. By constructing low-dimensional visualizations of these landscapes, which each comprise more than 16 million mRNA sequences, we show that such alternative codes radically alter the topological features of the network of high-fitness genotypes. Whereas the genetic codes that optimize evolvability depend to some extent on the detailed relationship between amino acid sequence and protein function, we also uncover general design principles for engineering nonstandard genetic codes for enhanced and diminished evolvability, which may facilitate directed protein evolution experiments and the bio-containment of synthetic organisms, respectively.

Choose Your Own Adventure: A Comprehensive Database of Reactions Catalyzed by Cytochrome P450 BM3 Variants

Cytochrome P450 BM3 monooxygenase is the topic of extensive research as many researchers have evolved this enzyme to generate a variety of products. However, the abundance of information on increasingly diversified variants of P450 BM3 that catalyze a broad array of chemistry is not in a format that enables easy extraction and interpretation. We present a database that categorizes variants by their catalyzed reactions and includes details about substrates to provide reaction context. This database of >1500 P450 BM3 variants is downloadable and machine-readable and includes instructions to maximize ease of gathering information. The database allows rapid identification of commonly reported substitutions, aiding researchers who are unfamiliar with the enzyme in identifying starting points for enzyme engineering. For those actively engaged in engineering P450 BM3, the database, along with this review, provides a powerful and user-friendly platform to understand, predict, and identify the attributes of P450 BM3 variants, encouraging the further engineering of this enzyme.

Dimer Stabilization by SpyTag/SpyCatcher Coupling of the Reductase Domains of a Chimeric P450 BM3 Monooxygenase from Bacillus spp. Improves its Stability, Activity, and Purification

The vast majority of known enzymes exist as oligomers, which often gives them high catalytic performance but at the same time imposes constraints on structural conformations and environmental conditions. An example of an enzyme with a complex architecture is the P450 BM3 monooxygenase CYP102A1 from Bacillus megaterium. Only active as a dimer, it is highly sensitive to dilution or common immobilization techniques. In this study, we engineered a thermostable P450BM3 chimera consisting of the heme domain of a CYP102A1 variant and the reductase domain of the homologous CYP102A3. The dimerization of the hybrid was even weaker compared to the corresponding CYP102A1 variant. To create a stable dimer, we covalently coupled the C-termini of two monomers of the chimera via SpyTag003/SpyCatcher003 interaction. As a result, purification, thermostability, pH stability, and catalytic activity were improved. Via a bioorthogonal two-step affinity purification, we obtained high purity (94 %) of the dimer-stabilized variant being robust against heme depletion. Long-term stability was increased with a half-life of over 2 months at 20 °C and 80-90 % residual activity after 2 months at 5 °C. Most catalytic features were retained with even an enhancement of the overall activity by ~2-fold compared to the P450BM3 chimera without SpyTag003/SpyCatcher003.

Efficiency aspects of regioselective testosterone hydroxylation with highly active CYP450-based whole-cell biocatalysts

Steroid hydroxylations belong to the industrially most relevant reactions catalysed by cytochrome P450 monooxygenases (CYP450s) due to the pharmacological relevance of hydroxylated derivatives. The implementation of respective bioprocesses at an industrial scale still suffers from several limitations commonly found in CYP450 catalysis, that is low turnover rates, enzyme instability, inhibition and toxicity related to the substrate(s) and/or product(s). Recently, we achieved a new level of steroid hydroxylation rates by introducing highly active testosterone-hydroxylating CYP450 BM3 variants together with the hydrophobic outer membrane protein AlkL into Escherichia coli-based whole-cell biocatalysts. However, the activity tended to decrease, which possibly impedes overall productivities and final product titres. In this study, a considerable instability was confirmed and subject to a systematic investigation regarding possible causes. In-depth evaluation of whole-cell biocatalyst kinetics and stability revealed a limitation in substrate availability due to poor testosterone solubility as well as inhibition by the main product 15β-hydroxytestosterone. Instability of CYP450 BM3 variants was disclosed as another critical factor, which is of general significance for CYP450-based biocatalysis. Presented results reveal biocatalyst, reaction and process engineering strategies auguring well for industrial implementation of the developed steroid hydroxylation platform.

Computational redesign of cytochrome P450 CYP102A1 for highly stereoselective omeprazole hydroxylation by UniDesign

Cytochrome P450 CYP102A1 is a prototypic biocatalyst that has great potential in chemical synthesis, drug discovery, and biotechnology. CYP102A1 variants engineered by directed evolution and/or rational design are capable of catalyzing the oxidation of a wide range of organic compounds. However, it is difficult to foresee the outcome of engineering CYP102A1 for a compound of interest. Here, we introduce UniDesign as a computational framework for enzyme design and engineering. We tested UniDesign by redesigning CYP102A1 for stereoselective metabolism of omeprazole (OMP), a proton pump inhibitor, starting from an active but nonstereoselective triple mutant (TM: A82F/F87V/L188Q). To shift stereoselectivity toward (R)-OMP, we computationally scanned three active site positions (75, 264, and 328) for mutations that would stabilize the binding of the transition state of (R)-OMP while destabilizing that of (S)-OMP and picked three variants, namely UD1 (TM/L75I), UD2 (TM/A264G), and UD3 (TM/A328V), for experimentation, based on computed energy scores and models. UD1, UD2, and UD3 exhibit high turnover rates of 55 ± 4.7, 84 ± 4.8, and 79 ± 5.7 min-1, respectively, for (R)-OMP hydroxylation, whereas the corresponding rates for (S)-OMP are only 2.2 ± 0.19, 6.0 ± 0.68, and 14 ± 2.8 min-1, yielding an enantiomeric excess value of 92, 87, and 70%, respectively. These results suggest the critical roles of L75I, A264G, and A328V in steering OMP in the optimal orientation for stereoselective oxidation and demonstrate the utility of UniDesign for engineering CYP102A1 to produce drug metabolites of interest. The results are discussed in the context of protein structures.

Evolvability-enhancing mutations in the fitness landscapes of an RNA and a protein

Can evolvability-the ability to produce adaptive heritable variation-itself evolve through adaptive Darwinian evolution? If so, then Darwinian evolution may help create the conditions that enable Darwinian evolution. Here I propose a framework that is suitable to address this question with available experimental data on adaptive landscapes. I introduce the notion of an evolvability-enhancing mutation, which increases the likelihood that subsequent mutations in an evolving organism, protein, or RNA molecule are adaptive. I search for such mutations in the experimentally characterized and combinatorially complete fitness landscapes of a protein and an RNA molecule. I find that such evolvability-enhancing mutations indeed exist. They constitute a small fraction of all mutations, which shift the distribution of fitness effects of subsequent mutations towards less deleterious mutations, and increase the incidence of beneficial mutations. Evolving populations which experience such mutations can evolve significantly higher fitness. The study of evolvability-enhancing mutations opens many avenues of investigation into the evolution of evolvability.

Engineering Catalytically Self-Sufficient P450s

The P450 superfamily comprises some of the most powerful and versatile enzymes for the site-selective oxidation of small molecules. One of the main drawbacks for the applications of the P450s in biotechnology is that the majority of these enzymes is multicomponent in nature and requires the presence of suitable redox partners to support their functions. Nevertheless, the discovery of several self-sufficient P450s, namely those from Classes VII and VIII, has served as an inspiration for fusion approaches to generate chimeric P450 systems that are self-sufficient. In this Perspective, we highlight the domain organizations of the Class VII and Class VIII P450 systems, summarize recent case studies in the engineering of catalytically self-sufficient P450s based on these systems, and outline outstanding challenges in the field, along with several emerging technologies as potential solutions.

Enhanced and specific epoxidation activity of P450 BM3 mutants for the production of high value terpene derivatives

Terpenes are natural molecules of valuable interest for different industrial applications. Cytochromes P450 enzymes can functionalize terpenoids to form high value oxidized derivatives in a green and sustainable manner, representing a valid alternative to chemical catalysis. In this work, an enhanced and specific epoxidation activity of cytochrome P450 BM3 mutants was found for the terpenes geraniol and linalool. This is the first report showing the epoxidation of linalool by P450 BM3 and its mutant A2 (Asp251Gly/Gln307His) with the formation of valuable oxide derivatives, highlighting the relevance of this enzymes for industrial applications.

Adaptive evolvability through direct selection instead of indirect, second-order selection

Can evolvability itself be the product of adaptive evolution? To answer this question is challenging, because any DNA mutation that alters only evolvability is subject to indirect, "second order" selection on the future effects of this mutation. Such indirect selection is weaker than "first-order" selection on mutations that alter fitness, in the sense that it can operate only under restrictive conditions. Here I discuss a route to adaptive evolvability that overcomes this challenge. Specifically, a recent evolution experiment showed that some mutations can enhance both fitness and evolvability through a combination of direct and indirect selection. Unrelated evidence from gene duplication and the evolution of gene regulation suggests that mutations with such dual effects may not be rare. Through such mutations, evolvability may increase at least in part because it provides an adaptive advantage. These observations suggest a research program on the adaptive evolution of evolvability, which aims to identify such mutations and to disentangle their direct fitness effects from their indirect effects on evolvability. If evolvability is itself adaptive, Darwinian evolution may have created more than life's diversity. It may also have helped create the very conditions that made the success of Darwinian evolution possible.

Valorization of Small Alkanes by Biocatalytic Oxyfunctionalization

The oxidation of alkanes into valuable chemical products is a vital reaction in organic synthesis. This reaction, however, is challenging, owing to the inertness of C-H bonds. Transition metal catalysts for C-H functionalization are frequently explored. Despite chemical alternatives, nature has also evolved powerful oxidative enzymes (e. g., methane monooxygenases, cytochrome P450 oxygenases, peroxygenases) that are capable of transforming C-H bonds under very mild conditions, with only the use of molecular oxygen or hydrogen peroxide as electron acceptors. Although progress in alkane oxidation has been reviewed extensively, little attention has been paid to small alkane oxidation. The latter holds great potential for the manufacture of chemicals. This Minireview provides a concise overview of the most relevant enzyme classes capable of small alkanes (C_<6 ) oxyfunctionalization, describes the essentials of the catalytic mechanisms, and critically outlines the current state-of-the-art in preparative applications.

Viral protein instability enhances host-range evolvability

Viruses are highly evolvable, but what traits endow this property? The high mutation rates of viruses certainly play a role, but factors that act above the genetic code, like protein thermostability, are also expected to contribute. We studied how the thermostability of a model virus, bacteriophage λ, affects its ability to evolve to use a new receptor, a key evolutionary transition that can cause host-range evolution. Using directed evolution and synthetic biology techniques we generated a library of host-recognition protein variants with altered stabilities and then tested their capacity to evolve to use a new receptor. Variants fell within three stability classes: stable, unstable, and catastrophically unstable. The most evolvable were the two unstable variants, whereas seven of eight stable variants were significantly less evolvable, and the two catastrophically unstable variants could not grow. The slowly evolving stable variants were delayed because they required an additional destabilizing mutation. These results are particularly noteworthy because they contradict a widely supported contention that thermostabilizing mutations enhance evolvability of proteins by increasing mutational robustness. Our work suggests that the relationship between thermostability and evolvability is more complex than previously thought, provides evidence for a new molecular model of host-range expansion evolution, and identifies instability as a potential predictor of viral host-range evolution.

Highly stereoselective and enantiodivergent synthesis of cyclopropylphosphonates with engineered carbene transferases

Organophosphonate compounds have represented a rich source of biologically active compounds, including enzyme inhibitors, antibiotics, and antimalarial agents. Here, we report the development of a highly stereoselective strategy for olefin cyclopropanation in the presence of a phosphonyl diazo reagent as carbene precursor. In combination with a ‘substrate walking’ protein engineering strategy, two sets of efficient and enantiodivergent myoglobin-based biocatalysts were developed for the synthesis of both (1R,2S) and (1S,2R) enantiomeric forms of the desired cyclopropylphosphonate ester products. This methodology enables the efficient transformation of a broad range of vinylarene substrates at a preparative scale (i.e. gram scale) with up to 99% de and ee. Mechanistic studies provide insights into factors that contribute to make this reaction inherently more challenging than hemoprotein-catalyzed olefin cyclopropanation with ethyl diazoacetate investigated previously. This work expands the range of synthetically useful, enzyme-catalyzed transformations and paves the way to the development of metalloprotein catalysts for abiological carbene transfer reactions involving non-canonical carbene donor reagents.

Two enantiocomplementary myoglobin-based carbene transfer biocatalysts were developed for the synthesis of cyclopropylphosphonate esters with high diastereo- and enantioselectivity and in high yields.

Learning Strategies in Protein Directed Evolution

Synthetic biology is a fast-evolving research field that combines biology and engineering principles to develop new biological systems for medical, pharmacological, and industrial applications. Synthetic biologists use iterative "design, build, test, and learn" cycles to efficiently engineer genetic systems that are reliable, reproducible, and predictable. Protein engineering by directed evolution can benefit from such a systematic engineering approach for various reasons. Learning can be carried out before starting, throughout or after finalizing a directed evolution project. Computational tools, bioinformatics, and scanning mutagenesis methods can be excellent starting points, while molecular dynamics simulations and other strategies can guide engineering efforts. Similarly, studying protein intermediates along evolutionary pathways offers fascinating insights into the molecular mechanisms shaped by evolution. The learning step of the cycle is not only crucial for proteins or enzymes that are not suitable for high-throughput screening or selection systems, but it is also valuable for any platform that can generate a large amount of data that can be aided by machine learning algorithms. The main challenge in protein engineering is to predict the effect of a single mutation on one functional parameter-to say nothing of several mutations on multiple parameters. This is largely due to nonadditive mutational interactions, known as epistatic effects-beneficial mutations present in a genetic background may not be beneficial in another genetic background. In this work, we provide an overview of experimental and computational strategies that can guide the user to learn protein function at different stages in a directed evolution project. We also discuss how epistatic effects can influence the success of directed evolution projects. Since machine learning is gaining momentum in protein engineering and the field is becoming more interdisciplinary thanks to collaboration between mathematicians, computational scientists, engineers, molecular biologists, and chemists, we provide a general workflow that familiarizes nonexperts with the basic concepts, dataset requirements, learning approaches, model capabilities and performance metrics of this intriguing area. Finally, we also provide some practical recommendations on how machine learning can harness epistatic effects for engineering proteins in an "outside-the-box" way.

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.

Comprehensive Structure-Activity Profiling of Micheliolide and its Targeted Proteome in Leukemia Cells via Probe-Guided Late-Stage C-H Functionalization

The plant-derived sesquiterpene lactone micheliolide was recently found to possess promising antileukemic activity, including the ability to target and kill leukemia stem cells. Efforts toward improving the biological activity of micheliolide and investigating its mechanism of action have been hindered by the paucity of preexisting functional groups amenable for late-stage derivatization of this molecule. Here, we report the implementation of a probe-based P450 fingerprinting strategy to rapidly evolve engineered P450 catalysts useful for the regio- and stereoselective hydroxylation of micheliolide at two previously inaccessible aliphatic positions in this complex natural product. Via P450-mediated chemoenzymatic synthesis, a broad panel of novel micheliolide analogs could thus be obtained to gain structure-activity insights into the effect of C2, C4, and C14 substitutions on the antileukemic activity of micheliolide, ultimately leading to the discovery of "micheliologs" with improved potency against acute myelogenic leukemia cells. These late-stage C-H functionalization routes could be further leveraged to generate a panel of affinity probes for conducting a comprehensive analysis of the protein targeting profile of micheliolide in leukemia cells via chemical proteomics analyses. These studies introduce new micheliolide-based antileukemic agents and shed new light onto the biomolecular targets and mechanism of action of micheliolide in leukemia cells. More broadly, this work showcases the value of the present P450-mediated C-H functionalization strategy for streamlining the late-stage diversification and elucidation of the biomolecular targets of a complex bioactive molecule.

Optimisation of Cytochrome P450 BM3 Assisted by Consensus-Guided Evolution

Cytochrome P450 enzymes have attracted much interest over the years given their ability to insert oxygen into saturated carbon-hydrogen bonds, a difficult feat to accomplish by traditional chemistry. Much of the activity in this field has centered on the bacterial enzyme CYP102A1, or BM3, from Bacillus megaterium, as it has shown itself capable of hydroxylating/acting upon a wide range of substrates, thereby producing industrially relevant pharmaceuticals, fine chemicals, and hormones. In addition, unlike most cytochromes, BM3 is both soluble and fused to its natural redox partner, thus facilitating its use. The industrial use of BM3 is however stifled by its instability and its requirement for the expensive NADPH cofactor. In this work, we added several mutations to the BM3 mutant R966D/W1046S that enhanced the turnover number achievable with the inexpensive cofactors NADH and NBAH. These new mutations, A769S, S847G, S850R, E852P, and V978L, are localized on the reductase domain of BM3 thus leaving the oxidase domain intact. For NBAH-driven reactions by new mutant NTD5, this led to a 5.24-fold increase in total product output when compared to the BM3 mutant R966D/W1046S. For reactions driven by NADH by new mutant NTD6, this enhanced total product output by as much as 2.3-fold when compared to the BM3 mutant R966D/W1046S. We also demonstrated that reactions driven by NADH with the NTD6 mutant not only surpassed total product output achievable by wild-type BM3 with NADPH but also retained the ability to use this latter cofactor with greater total product output as well.

A Series of Crystallographically Characterized Linear and Branched σ-Alkane Complexes of Rhodium: From Propane to 3-Methylpentane

Using solid-state molecular organometallic (SMOM) techniques, in particular solid/gas single-crystal to single-crystal reactivity, a series of σ-alkane complexes of the general formula [Rh(Cy2PCH2CH2PCy2)(ηn:ηm-alkane)][BArF4] have been prepared (alkane = propane, 2-methylbutane, hexane, 3-methylpentane; ArF = 3,5-(CF3)2C6H3). These new complexes have been characterized using single crystal X-ray diffraction, solid-state NMR spectroscopy and DFT computational techniques and present a variety of Rh(I)···H-C binding motifs at the metal coordination site: 1,2-η2:η2 (2-methylbutane), 1,3-η2:η2 (propane), 2,4-η2:η2 (hexane), and 1,4-η1:η2 (3-methylpentane). For the linear alkanes propane and hexane, some additional Rh(I)···H-C interactions with the geminal C-H bonds are also evident. The stability of these complexes with respect to alkane loss in the solid state varies with the identity of the alkane: from propane that decomposes rapidly at 295 K to 2-methylbutane that is stable and instead undergoes an acceptorless dehydrogenation to form a bound alkene complex. In each case the alkane sits in a binding pocket defined by the {Rh(Cy2PCH2CH2PCy2)}+ fragment and the surrounding array of [BArF4]- anions. For the propane complex, a small alkane binding energy, driven in part by a lack of stabilizing short contacts with the surrounding anions, correlates with the fleeting stability of this species. 2-Methylbutane forms more short contacts within the binding pocket, and as a result the complex is considerably more stable. However, the complex of the larger 3-methylpentane ligand shows lower stability. Empirically, there therefore appears to be an optimal fit between the size and shape of the alkane and overall stability. Such observations are related to guest/host interactions in solution supramolecular chemistry and the holistic role of 1°, 2°, and 3° environments in metalloenzymes.

Pervasive cooperative mutational effects on multiple catalytic enzyme traits emerge via long-range conformational dynamics

Multidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.

Connecting conformational dynamics and epistasis has so far been limited to a few proteins and a single fitness trait. Here, the authors provide evidence of positive epistasis on multiple catalytic traits in the evolution and dynamics of engineered cytochrome P450 monooxygenase, offering insights for in silico protein design.

Selectivity of Rh⋅⋅⋅H-C Binding in a σ-Alkane Complex Controlled by the Secondary Microenvironment in the Solid State

Single-crystal to single-crystal solid-state molecular organometallic (SMOM) techniques are used for the synthesis and structural characterization of the σ-alkane complex [Rh(tBu2 PCH2 CH2 CH2 PtBu2 )(η2 ,η2 -C7 H12 )][BArF 4 ] (ArF =3,5-(CF3 )2 C6 H3 ), in which the alkane (norbornane) binds through two exo-C-H⋅⋅⋅Rh interactions. In contrast, the bis-cyclohexyl phosphine analogue shows endo-alkane binding. A comparison of the two systems, supported by periodic DFT calculations, NCI plots and Hirshfeld surface analyses, traces this different regioselectivity to subtle changes in the local microenvironment surrounding the alkane ligand. A tertiary periodic structure supporting a secondary microenvironment that controls binding at the metal site has parallels with enzymes. The new σ-alkane complex is also a catalyst for solid/gas 1-butene isomerization, and catalyst resting states are identified for this.

Computational design of enzymes for biotechnological applications

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

Kinetic Evidence for an Induced Fit Mechanism in the Binding of the Substrate Camphor by Cytochrome P450cam

Bacterial cytochrome P450 (P450) 101A1 (P450cam) has served as a prototype among the P450 enzymes and has high catalytic activity towards its cognate substrate, camphor. X-ray crystallography and NMR and IR spectroscopy have demonstrated the existence of multiple conformations of many P450s, including P450cam. Kinetic studies have indicated that substrate binding to several P450s is dominated by a conformational selection process, in which the substrate binds an individual conformer(s) of the unliganded enzyme. P450cam was found to differ in that binding of the substrate camphor is dominated by an induced fit mechanism, in which the enzyme binds camphor and then changes conformation, as evidenced by the equivalence of binding eigenvalues observed when varying both camphor and P450cam concentrations. The accessory protein putidaredoxin had no effect on substrate binding. Estimation of the rate of dissociation of the P450cam·camphor complex (15 s-1) and fitting of the data yield a minimal kinetic mechanism in which camphor binds (1.5 × 107 M-1 s-1) and the initial P450cam•camphor complex undergoes a reversible equilibrium (k forward 112 s-1, k reverse 28 s-1) to a final complex. This induced fit mechanism differs from those reported for several mammalian P450s and bacterial P450BM-3, indicative of the diversity of how P450s recognize multiple substrates. However, similar behavior was not observed with the alternate substrates (+)-α-pinene and 2-adamantanone, which probably utilize a conformational selection process.

Engineered dual selection for directed evolution of SpCas9 PAM specificity

The widely used Streptococcus pyogenes Cas9 (SpCas9) nuclease derives its DNA targeting specificity from protein-DNA contacts with protospacer adjacent motif (PAM) sequences, in addition to base-pairing interactions between its guide RNA and target DNA. Previous reports have established that the PAM specificity of SpCas9 can be altered via positive selection procedures for directed evolution or other protein engineering strategies. Here we exploit in vivo directed evolution systems that incorporate simultaneous positive and negative selection to evolve SpCas9 variants with commensurate or improved activity on NAG PAMs relative to wild type and reduced activity on NGG PAMs, particularly YGG PAMs. We also show that the PAM preferences of available evolutionary intermediates effectively determine whether similar counterselection PAMs elicit different selection stringencies, and demonstrate that negative selection can be specifically increased in a yeast selection system through the fusion of compensatory zinc fingers to SpCas9.

Engineered SAM Synthetases for Enzymatic Generation of AdoMet Analogs with Photocaging Groups and Reversible DNA Modification in Cascade Reactions

Methylation and demethylation of DNA, RNA and proteins has emerged as a major regulatory mechanism. Studying the function of these modifications would benefit from tools for their site‐specific inhibition and timed removal. S‐Adenosyl‐L‐methionine (AdoMet) analogs in combination with methyltransferases (MTases) have proven useful to map or block and release MTase target sites, however their enzymatic generation has been limited to aliphatic groups at the sulfur atom. We engineered a SAM synthetase from Cryptosporidium hominis (PC‐ChMAT) for efficient generation of AdoMet analogs with photocaging groups that are not accepted by any WT MAT reported to date. The crystal structure of PC‐ChMAT at 1.87 Å revealed how the photocaged AdoMet analog is accommodated and guided engineering of a thermostable MAT from Methanocaldococcus jannaschii. PC‐MATs were compatible with DNA‐ and RNA‐MTases, enabling sequence‐specific modification (“writing”) of plasmid DNA and light‐triggered removal (“erasing”).

Expanding the applicability of cytochrome P450s and other haemoproteins

Cytochrome P450BM3 has long been regarded as a promising candidate for use as a biocatalyst, owing to its excellent efficiency for the hydroxylation of unactivated C-H bonds. However, because of its high substrate specificity, its possible applications have been severely limited. Consequently, various approaches have been proposed to overcome the enzyme's natural limitations, thereby expanding its substrate scope to encompass non-native substrates, evoking chemoselectivity, regioselectivity and stereoselectivity and enabling previously inaccessible chemical conversions. Herein, these approaches will be classified into three categories: (1) mutagenesis including directed evolution, (2) haem substitution with artificial cofactors and (3) use of substrate mimics, 'decoy molecules'. Herein, we highlight the representative work that has been conducted in above three categories for discussion of the future outlook of P450BM3 in green chemistry.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.

Automated Continuous Evolution of Proteins in Vivo

We present automated continuous evolution (ACE), a platform for the hands-free directed evolution of biomolecules. ACE pairs OrthoRep, a genetic system for continuous targeted mutagenesis of user-selected genes in vivo, with eVOLVER, a scalable and automated continuous culture device for precise, multiparameter regulation of growth conditions. By implementing real-time feedback-controlled tuning of selection stringency with eVOLVER, genes of interest encoded on OrthoRep autonomously traversed multimutation adaptive pathways to reach desired functions, including drug resistance and improved enzyme activity. The durability, scalability, and speed of biomolecular evolution with ACE should be broadly applicable to protein engineering as well as prospective studies on how selection parameters and schedules shape adaptation.

Crystals in Minutes: Instant On-Site Microcrystallisation of Various Flavours of the CYP102A1 (P450BM3) Haem Domain

Despite CYP102A1 (P450BM3) representing one of the most extensively researched metalloenzymes, crystallisation of its haem domain upon modification can be a challenge. Crystal structures are indispensable for the efficient structure-based design of P450BM3 as a biocatalyst. The abietane diterpenoid derivative N-abietoyl-l-tryptophan (AbiATrp) is an outstanding crystallisation accelerator for the wild-type P450BM3 haem domain, with visible crystals forming within 2 hours and diffracting to a near-atomic resolution of 1.22 Å. Using these crystals as seeds in a cross-microseeding approach, an assortment of P450BM3 haem domain crystal structures, containing previously uncrystallisable decoy molecules and diverse artificial metalloporphyrins binding various ligand molecules, as well as heavily tagged haem-domain variants, could be determined. Some of the structures reported herein could be used as models of different stages of the P450BM3 catalytic cycle.

Directed evolution methods for overcoming trade-offs between protein activity and stability

Engineered proteins are being widely developed and employed in applications ranging from enzyme catalysts to therapeutic antibodies. Directed evolution, an iterative experimental process composed of mutagenesis and library screening, is a powerful technique for enhancing existing protein activities and generating entirely new ones not observed in nature. However, the process of accumulating mutations for enhanced protein activity requires chemical and structural changes that are often destabilizing, and low protein stability is a significant barrier to achieving large enhancements in activity during multiple rounds of directed evolution. Here we highlight advances in understanding the origins of protein activity/stability trade-offs for two important classes of proteins (enzymes and antibodies) as well as innovative experimental and computational methods for overcoming such trade-offs. These advances hold great potential for improving the generation of highly active and stable proteins that are needed to address key challenges related to human health, energy and the environment.

Promiscuous Ribozymes and Their Proposed Role in Prebiotic Evolution

The ability of enzymes, including ribozymes, to catalyze side reactions is believed to be essential to the evolution of novel biochemical activities. It has been speculated that the earliest ribozymes, whose emergence marked the origin of life, were low in activity but high in promiscuity, and that these early ribozymes gave rise to specialized descendants with higher activity and specificity. Here, we review the concepts related to promiscuity and examine several cases of highly promiscuous ribozymes. We consider the evidence bearing on the question of whether de novo ribozymes would be quantitatively more promiscuous than later evolved ribozymes or protein enzymes. We suggest that while de novo ribozymes appear to be promiscuous in general, they are not obviously more promiscuous than more highly evolved or active sequences. Promiscuity is a trait whose value would depend on selective pressures, even during prebiotic evolution.

A Continuing Career in Biocatalysis: Frances H. Arnold

On the occasion of Professor Frances H. Arnold's recent acceptance of the 2018 Nobel Prize in Chemistry, we honor her numerous contributions to the fields of directed evolution and biocatalysis. Arnold pioneered the development of directed evolution methods for engineering enzymes as biocatalysts. Her highly interdisciplinary research has provided a ground not only for understanding the mechanisms of enzyme evolution but also for developing commercially viable enzyme biocatalysts and biocatalytic processes. In this Account, we highlight some of her notable contributions in the past three decades in the development of foundational directed evolution methods and their applications in the design and engineering of enzymes with desired functions for biocatalysis. Her work has created a paradigm shift in the broad catalysis field.

A 96-multiplex capillary electrophoresis screening platform for product based evolution of P450 BM3

The main challenge that prevents a broader application of directed enzyme evolution is the lack of high-throughput screening systems with universal product analytics. Most directed evolution campaigns employ screening systems based on colorimetric or fluorogenic surrogate substrates or universal quantification methods such as nuclear magnetic resonance spectroscopy or mass spectrometry, which have not been advanced to achieve a high-throughput. Capillary electrophoresis with a universal UV-based product detection is a promising analytical tool to quantify product formation. Usage of a multiplex system allows the simultaneous measurement with 96 capillaries. A 96-multiplexed capillary electrophoresis (MP-CE) enables a throughput that is comparable to traditional direct evolution campaigns employing 96-well microtiter plates. Here, we report for the first time the usage of a MP-CE system for directed P450 BM3 evolution towards increased product formation (oxidation of alpha-isophorone to 4-hydroxy-isophorone; highest reached total turnover number after evolution campaign: 7120 mol_4-OH mol_P450⁻¹). The MP-CE platform was 3.5-fold more efficient in identification of beneficial variants than the standard cofactor (NADPH) screening system.

Hoodwinking Cytochrome P450BM3 into Hydroxylating Non-Native Substrates by Exploiting Its Substrate Misrecognition

Bacterial cytochrome P450s (P450s) are at the focus of attention as potential biocatalysts for applications in green synthetic chemistry, as they possess high activity for the hydroxylation of inert substrate C-H bonds. The high activity of bacterial P450s, such as P450BM3, is chiefly due to their high substrate specificity, and consequently, the catalytic activity of P450BM3 toward non-native substrates is very low, limiting the utility of bacterial P450s as biocatalysts. To enable oxidation of non-native substrates by P450BM3 without any mutagenesis, we have developed a series of "decoy molecules", inert dummy substrates, with structures that resemble those of the native substrates. Decoy molecules fool P450BM3 into generating the active species, so-called Compound I, enabling the catalytic oxidation of non-native substrates other than fatty acids. Perfluorinated carboxylic acids (PFCs) serve as decoy molecules to initiate the activation of molecular oxygen in the same manner as long-alkyl-chain fatty acids, due to their structural similarity, and induce the generation of Compound I, but, unlike the native substrates, PFCs are not oxidizable by Compound I, allowing the hydroxylation of non-native substrates, such as gaseous alkanes and benzene. The catalytic activity for non-native substrate hydroxylation was significantly enhanced by employing second generation decoy molecules, PFCs modified with amino acids (PFC-amino acids). Cocrystals of P450BM3 with PFC9-Trp revealed clear electron density in the fatty-acid-binding channel that was readily assigned to PFC9-Trp. The alkyl chain terminus of PFC9-Trp does not reach the active site owing to multiple hydrogen bonding interactions between the carboxyl and carbonyl groups of PFC9-Trp and amino acids located at the entrance of the substrate binding channel of P450BM3 that fix it in place. The remaining space above the heme after binding of PFC9-Trp can be utilized to accommodate non-native substrates. Further developments revealed that third generation decoy molecules, N-acyl amino acids, such as pelargonoyl-l-phenylalanine (C9-Phe), can serve as decoy molecules, indicating that the rationale "fluorination is required for decoy molecule function" can be safely discarded. Diverse carboxylic acids including dipeptides could now be exploited as building blocks, and a library of decoy molecules possessing diverse structures was prepared. Among the third-generation decoy molecules examined N-enanthyl-l-proline modified with l-phenylalanine (C7-Pro-Phe) afforded the maximum turnover rate for benzene hydroxylation. The structural diversity of third-generation decoy molecules was also utilized to control the stereoselectivity of hydroxylation for the benzylic hydroxylation of Indane, showing that decoy molecules can alter stereoselectivity. As both the catalytic activity and enantioselectivity are dependent upon the structure of the decoy molecules, their design allows us to regulate reactions catalyzed by wild-type enzymes. Furthermore, decoy molecules can also activate intracellular P450BM3, allowing the use of E. coli expressing wild-type P450BM3 as an efficient whole-cell bioreactor. It should be noted that Mn-substituted full-length P450BM3 (Mn-P450BM3) is also active for the hydroxylation of propane in which the regioselectivity diverged from that of Fe-P450BM3. The results summarized in this Account represent good examples of how the reactive properties of P450BM3 can be controlled for the monooxygenation of non-native substrates in vitro as well as in vivo to expand the potential of P450BM3.

Mechanisms of Cytochrome P450-Catalyzed Oxidations

Enzymes are complex biological catalysts and are critical to life. Most oxidations of chemicals are catalyzed by cytochrome P450 (P450, CYP) enzymes, which generally utilize mixed-function oxidase stoichiometry, utilizing pyridine nucleotides as electron donors: NAD(P)H + O₂ + R → NAD(P)⁺ + RO + H₂O (where R is a carbon substrate and RO is an oxidized product). The catalysis of oxidations is largely understood in the context of the heme iron-oxygen complex generally referred to as Compound I, formally FeO³⁺, whose basis was in peroxidase chemistry. Many X-ray crystal structures of P450s are now available (≥ 822 structures from ≥146 different P450s) and have helped in understanding catalytic specificity. In addition to hydroxylations, P450s catalyze more complex oxidations, including C-C bond formation and cleavage. Enzymes derived from P450s by directed evolution can even catalyze more unusual reactions, e.g. cyclopropanation. Current P450 questions under investigation include the potential role of the intermediate Compound 0 (formally Fe^III-O₂ ^-) in catalysis of some reactions, the roles of high- and low-spin forms of Compound I, the mechanism of desaturation, the roles of open and closed structures of P450s in catalysis, the extent of processivity in multi-step oxidations, and the role of the accessory protein cytochrome b ₅. More global questions include exactly how structure drives function, prediction of catalysis, and roles of multiple protein conformations.

Evolutionary repurposing of a sulfatase: A new Michaelis complex leads to efficient transition state charge offset

The recruitment and evolutionary optimization of promiscuous enzymes is key to the rapid adaptation of organisms to changing environments. Our understanding of the precise mechanisms underlying enzyme repurposing is, however, limited: What are the active-site features that enable the molecular recognition of multiple substrates with contrasting catalytic requirements? To gain insights into the molecular determinants of adaptation in promiscuous enzymes, we performed the laboratory evolution of an arylsulfatase to improve its initially weak phenylphosphonate hydrolase activity. The evolutionary trajectory led to a 100,000-fold enhancement of phenylphosphonate hydrolysis, while the native sulfate and promiscuous phosphate mono- and diester hydrolyses were only marginally affected (≤50-fold). Structural, kinetic, and in silico characterizations of the evolutionary intermediates revealed that two key mutations, T50A and M72V, locally reshaped the active site, improving access to the catalytic machinery for the phosphonate. Measured transition state (TS) charge changes along the trajectory suggest the creation of a new Michaelis complex (E•S, enzyme-substrate), with enhanced leaving group stabilization in the TS for the promiscuous phosphonate (β_leavinggroup from -1.08 to -0.42). Rather than altering the catalytic machinery, evolutionary repurposing was achieved by fine-tuning the molecular recognition of the phosphonate in the Michaelis complex, and by extension, also in the TS. This molecular scenario constitutes a mechanistic alternative to adaptation solely based on enzyme flexibility and conformational selection. Instead, rapid functional transitions between distinct chemical reactions rely on the high reactivity of permissive active-site architectures that allow multiple substrate binding modes.

Enhancement of Metallosphaera sedula Bioleaching by Targeted Recombination and Adaptive Laboratory Evolution

Thermophilic and lithoautotrophic archaea such as Metallosphaera sedula occupy acidic, metal-rich environments and are used in biomining processes. Biotechnological approaches could accelerate these processes and improve metal recovery by biomining organisms, but systems for genetic manipulation in these organisms are currently lacking. To gain a better understanding of the interplay between metal resistance, autotrophy, and lithotrophic metabolism, a genetic system was developed for M. sedula and used to evaluate parameters governing the efficiency of copper bioleaching. Additionally, adaptive laboratory evolution was used to select for naturally evolved M. sedula cell lines with desirable phenotypes for biomining, and these adapted cell lines were shown to have increased bioleaching capacity and efficiency. Genomic methods were used to analyze mutations that led to resistance in the experimentally evolved cell lines, while transcriptomics was used to examine changes in stress-inducible gene expression specific to the environmental conditions.

Enzyme engineering: reaching the maximal catalytic efficiency peak

The practical need for highly efficient enzymes presents new challenges in enzyme engineering, in particular, the need to improve catalytic turnover (k_cat) or efficiency (k_cat/K_M) by several orders of magnitude. However, optimizing catalysis demands navigation through complex and rugged fitness landscapes, with optimization trajectories often leading to strong diminishing returns and dead-ends. When no further improvements are observed in library screens or selections, it remains unclear whether the maximal catalytic efficiency of the enzyme (the catalytic 'fitness peak') has been reached; or perhaps, an alternative combination of mutations exists that could yield additional improvements. Here, we discuss fundamental aspects of the process of catalytic optimization, and offer practical solutions with respect to overcoming optimization plateaus.