Substrate Tunnel Engineering Aided by X-ray Crystallography and Functional Dynamics Swaps the Function of MIO-Enzymes

Navigating the landscape of enzyme design: from molecular simulations to machine learning

Global environmental issues and sustainable development call for new technologies for fine chemical synthesis and waste valorization. Biocatalysis has attracted great attention as the alternative to the traditional organic synthesis. However, it is challenging to navigate the vast sequence space to identify those proteins with admirable biocatalytic functions. The recent development of deep-learning based structure prediction methods such as AlphaFold2 reinforced by different computational simulations or multiscale calculations has largely expanded the 3D structure databases and enabled structure-based design. While structure-based approaches shed light on site-specific enzyme engineering, they are not suitable for large-scale screening of potential biocatalysts. Effective utilization of big data using machine learning techniques opens up a new era for accelerated predictions. Here, we review the approaches and applications of structure-based and machine-learning guided enzyme design. We also provide our view on the challenges and perspectives on effectively employing enzyme design approaches integrating traditional molecular simulations and machine learning, and the importance of database construction and algorithm development in attaining predictive ML models to explore the sequence fitness landscape for the design of admirable biocatalysts.

Engineered Phenylalanine Ammonia-Lyases for the Enantioselective Synthesis of Aspartic Acid Derivatives

Biocatalytic hydroamination of alkenes is an efficient and selective method to synthesize natural and unnatural amino acids. Phenylalanine ammonia-lyases (PALs) have been previously engineered to access a range of substituted phenylalanines and heteroarylalanines, but their substrate scope remains limited, typically including only arylacrylic acids. Moreover, the enantioselectivity in the hydroamination of electron-deficient substrates is often poor. Here, we report the structure-based engineering of PAL from Planctomyces brasiliensis (PbPAL), enabling preparative-scale enantioselective hydroaminations of previously inaccessible yet synthetically useful substrates, such as amide- and ester-containing fumaric acid derivatives. Through the elucidation of cryo-electron microscopy (cryo-EM) PbPAL structure and screening of the structure-based mutagenesis library, we identified the key active site residue L205 as pivotal for dramatically enhancing the enantioselectivity of hydroamination reactions involving electron-deficient substrates. Our engineered PALs demonstrated exclusive α-regioselectivity, high enantioselectivity, and broad substrate scope. The potential utility of the developed biocatalysts was further demonstrated by a preparative-scale hydroamination yielding tert-butyl protected l-aspartic acid, widely used as intermediate in peptide solid-phase synthesis.

Phenylalanine ammonia-lyases and 4-coumaric acid coenzyme A ligases in Chara braunii, Marchantia polymorpha, and Physcomitrium patens as extant model organisms for plant terrestrialization

The conquest of land posed severe problems to plants which they had to cope with by adapting biosynthetic capacities. Adaptations to respond to UV irradiation, water loss, pathogen and herbivore defense, and the earth's pull were essential. Chemical compounds alleviating these problems can be synthesized by the phenylpropanoid pathway, the core of which are three enzymes: phenylalanine ammonia-lyase (PAL), cinnamic acid 4-hydroxylase, and 4-coumaric acid coenzyme A-ligase (4CL). The genomes of model organisms, Chara braunii as aquatic alga and the two bryophytes Physcomitrium patens and Marchantia polymorpha, were searched for sequences encoding PAL and 4CL and selected sequences heterologously expressed in Escherichia coli for biochemical characterization. Several possible isoforms were identified for both enzymes in Marchantia polymorpha and Physcomitrium patens, while only one or two isoforms could be retrieved for Chara braunii. Active forms of both enzymes were found in all three organisms, although the catalytic efficiencies varied in a wide range. l-Phenylalanine was accepted as best substrate by all PAL-like enzymes, despite annotations in some cases suggesting different activities. The substrate spectrum of 4CLs was more diverse, but caffeic and/or 4-coumaric acids generally were the best-accepted substrates. Our investigations show that PAL and 4CL, important enzymes for the formation of phenolic compounds, are present and active in extant charophytes and bryophytes as model organisms for the conquest of land.

Efficient and easible biocatalysts: Strategies for enzyme improvement. A review

Owing to the environmental friendliness and vast advantages that enzymes offer in the biotechnology and industry fields, biocatalysts are a prolific investigation field. However, the low catalytic activity, stability, and specific selectivity of the enzyme limit the range of the reaction enzymes involved in. A comprehensive understanding of the protein structure and dynamics in terms of molecular details enables us to tackle these limitations effectively and enhance the catalytic activity by enzyme engineering or modifying the supports and solvents. Along with different strategies including computational, enzyme engineering based on DNA recombination, enzyme immobilization, additives, chemical modification, and physicochemical modification approaches can be promising for the wide spread of industrial enzyme usage. This is attributed to the successful application of biocatalysts in industrial and synthetic processes requires a system that exhibits stability, activity, and reusability in a continuous flow process, thereby reducing the production cost. The main goal of this review is to display relevant approaches for improving enzyme characteristics to overcome their industrial application.

Engineered Zea mays phenylalanine ammonia-lyase for improve the catalytic efficiency of biosynthesis trans-cinnamic acid and p-coumaric acid

Phenylalanine ammonia-lyase (PAL) plays a pivotal role in the biosynthesis of phenylalanine. PAL from Zea mays (ZmPAL2) exhibits a bi-function of direct deamination of L-phenylalanine (L-Phe) or L-tyrosine(-L-Tyr) to form trans-cinnamic acid or p-coumaric acid. trans-Cinnamic acid and p-coumaric acid are mainly used in flavors and fragrances, food additives, pharmaceutical and other fields. Here, the Activity of ZmPAL2 toward L-Phe or L-Tyr was improved by using semi-rational and rational designs. The catalytic efficiency (kcat/Km) of mutant PT10 (V258I/I459V/Q484N) against L-Phe was 30.8 μM-1 s-1, a 4.5-fold increase compared to the parent, and the catalytic efficiency of mutant PA1 (F135H/I459L) to L-tyrosine exhibited 8.6 μM-1 s-1, which was 1.6-fold of the parent. The yield of trans-cinnamic acid in PT10 reached 30.75 g/L with a conversion rate of 98%. Meanwhile, PA1 converted L-Tyr to yield 3.12 g/L of p-coumaric acid with a conversion rate of 95%. Suggesting these two engineered ZmPAL2 to be valuable biocatalysts for the synthesis of trans-cinnamic acid and p-coumaric acid. In addition, MD simulations revealed that the underlying mechanisms of the increased catalytic efficiency of both mutant PT10 and PA1 are attributed to the substrate remaining stable within the pocket and closer to the catalytically active site. This also provides a new perspective on engineered PAL.

Roles of the N-terminal motif in improving the activity and soluble expression of phenylalanine ammonia lyases in Escherichia coli

Phenylalanine ammonia-lyase (PAL) has various applications in fine chemical manufacturing and the pharmaceutical industry. In particular, PAL derived from Anabaena variabilis (AvPAL) is used as a therapeutic agent to the treat phenylketonuria in clinical settings. In this study, we aligned the amino acid sequences of AvPAL and PAL derived from Nostoc punctiforme (NpPAL) to obtain several mutants with enhanced activity, expression yield, and thermal stability via amino acid substitution and saturation mutagenesis at the N-terminal position. Enzyme kinetic experiments revealed that the kcat values of NpPAL-N2K, NpPAL-I3T, and NpPAL-T4L mutants were increased to 3.2-, 2.8-, and 3.3-fold that of the wild-type, respectively. Saturation mutagenesis of the fourth amino acid in AvPAL revealed that the kcat values of AvPAL-L4N, AvPAL-L4P, AvPAL-L4Q and AvPAL-L4S increased to 4.0-, 3.7-, 3.6-, and 3.2-fold, respectively. Additionally, the soluble protein yield of AvPAL-L4K increased to approximately 14 mg/L, which is approximately 3.5-fold that of AvPAL. Molecular dynamics studies further revealed that maintaining the attacking state of the reaction and N-terminal structure increased the rate of catalytic reaction and improved the solubility of proteins. These findings provide new insights for the rational design of PAL in the future.

Comparative biochemical characterization of mammalian-derived CYP11A1s with cholesterol side-chain cleavage activities

Steroid drugs, the second largest class of pharmaceuticals after antibiotics, have shown significant anti-inflammatory, anti-allergic, and endocrine-regulating effects. A group of cytochrome P450 enzymes, namely, CYP11A1 isoenzymes from different organisms are capable of converting cholesterol into pregnenolone, which is a pivotal reaction in both steroid metabolism and (bio)synthetic network of steroid products. However, the low activity of CYP11A1s greatly restricts the industrial application of these cholesterol side-chain cleavage enzymes. Herein, we investigate ten CYP11A1 enzymes of different origins and in vitro characterize two CYP11A1s with a relatively higher expression level from Capra hircus and Sus scrofa, together with the CYP11A1s from Homo sapiens and Bos taurus as references. Towards five selected sterol substrates with different side chain structures, S. scrofa CYP11A1 displays relatively higher activities. Through redox partners combination screening, we reveal the optimal redox partner pair of S. scrofa adrenodoxin and C. hircus adrenodoxin reductase. Moreover, the semi-rational mutagenesis for the active sites and substrate entrance channels of human and bovine CYP11A1s is performed based on comparative analysis of their crystal structures. The mutant mBtCYP11A1-Q377A derived from mature B. taurus CYP11A1 shows a 1.46 times higher activity than the wild type enzyme. These results not only demonstrate the tunability of the highly conserved CYP11A1 isoenzymes, but also lay a foundation for the following engineering efforts on these industrially relevant P450 enzymes.

Enhancing Enzyme Activity and Thermostability of Bacillus amyloliquefaciens Chitosanase BaCsn46A Through Saturation Mutagenesis at Ser196

Chitosanase plays an important role in chitooligosaccharides (COS) production. We found that the chitosanase (BaCsn46A) of Bacillus amyloliquefacien was a good candidate for chitosan hydrolysis of COS. In order to further improve the enzyme properties of BaCsn46A, the S196 located near the active center was found to be a critical site impacts on enzyme properties by sequence alignment analysis. Herein, saturation mutation was carried out to study role of 196 site on BaCsn46A catalytic function. Compared with WT, the specific enzyme activity of S196A increased by 118.79%, and the thermostability of S196A was much higher than WT. In addition, we found that the enzyme activity of S196P was 2.41% of that of WT, indicating that the type of amino acid in 196 site could significant affect the catalytic activity and thermostability of BaCsn46A. After molecular docking analysis we found that the increase in hydrogen bonds and decrease in unfavorable bonds interacting with the substrate were the main reason for the change of enzyme properties which is valuable for future studies on Bacillus species chitosanase.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

Inside Out Computational Redesign of Cavities for Improving Thermostability and Catalytic Activity of Rhizomucor Miehei Lipase

Cavities are created by hydrophobic interactions between residue side chain atoms during the folding of enzymes. Redesigning cavities can improve the thermostability and catalytic activity of the enzyme; however, the synergistic effect of cavities remains unclear. In this study, Rhizomucor miehei lipase (RML) was used as a model to explore volume fluctuation and spatial distribution changes of the internal cavities, which could reveal the roles of internal cavities in the thermostability and catalytic activity. We present an inside out cavity engineering (CE) strategy based on computational techniques to explore how changes in the volumes and spatial distribution of cavities affect the thermostability and catalytic activity of the enzyme. We obtained 12 single-point mutants, among which the melting temperatures (Tm) of 8 mutants showed an increase of more than 2°C. Sixteen multipoint mutations were further designed by spatial distribution rearrangement of internal cavities. The Tm of the most stable triple variant, with mutations including T21V (a change of T to V at position 21), S27A, and T198L (T21V/S27A/T198L), was elevated by 11.0°C, together with a 28.7-fold increase in the half-life at 65°C and a specific activity increase of 9.9-fold (up to 5,828 U mg-1), one of the highest lipase activities reported. The possible mechanism of decreased volumes and spatial rearrangement of the internal cavities improved the stability of the enzyme, optimizing the outer substrate tunnel to improve the catalytic efficiency. Overall, the inside out computational redesign of cavities method could help to deeply understand the effect of cavities on enzymatic stability and activity, which would be beneficial for protein engineering efforts to optimize natural enzymes. IMPORTANCE In the present study, R. miehei lipase, which is widely used in various industries, provides an opportunity to explore the effects of internal cavities on the thermostability and catalytic activity of enzymes. Here, we execute high hydrostatic pressure molecular dynamics (HP-MD) simulations to screen the critical internal cavity and reshape the internal cavities through site-directed mutation. We show that as the global internal cavity volume decreases, cavity rearrangement can improve the stability of the protein while optimizing the substrate channel to improve the catalytic efficiency. Our results provide significant insights into understanding the mechanism of action of the internal cavity. Our strategy is expected to be applied to other enzymes to promote increases in thermostability and catalytic activity.

Phenylalanine ammonia-lyases: combining protein engineering and natural diversity

In this study, rational design and saturation mutagenesis efforts for engineering phenylalanine ammonia-lyase from Petroselinum crispum (PcPAL) provided tailored PALs active towards challenging, highly valuable di-substituted substrates, such as the L-DOPA precursor 3,4-dimethoxy-L-phenylalanine or the 3-bromo-4-methoxy-phenylalanine. The rational design approach and saturation mutagenesis strategy unveiled identical PcPAL variants of improved activity, highlighting the limited mutational variety of the substrate specificity-modulator residues, L134, F137, I460 of PcPAL. Due to the restricted catalytic efficiency of the best performing L134A/I460V and F137V/I460V PcPAL variants, we imprinted these beneficial mutations to PALs of different origins. The variants of PALs from Arabidopsis thaliana (AtPAL) and Anabaena variabilis (AvPAL) showed higher catalytic efficiency than their PcPAL homologues. Further, the engineered PALs were also compared in terms of catalytic efficiency with a novel aromatic ammonia-lyase from Loktanella atrilutea (LaAAL), close relative of the metagenome-derived aromatic ammonia-lyase AL-11, reported recently to possess atypically high activity towards substrates with electron-donor aromatic substituents. Indeed, LaAAL outperformed the engineered Pc/At/AvPALs in the production of 3,4-dimethoxy-L-phenylalanine; however, in case of 3-bromo-4-methoxy derivatives it showed no activity, with computational results supporting the occurrence of steric hindrance. Transferring the unique array of selectivity modulator residues from LaAAL to the well-characterized PALs did not enhance their activity towards the targeted substrates. Moreover, applying the rational design strategy valid for these well-characterized PALs to LaAAL decreased its activity. These results suggest that distinct tailoring rationale is required for LaAAL/AL-11-like aromatic ammonia-lyases, which might represent a distinct PAL subclass, with natural reaction and substrate scope modified through evolutionary processes. KEY POINTS: • PAL-activity for challenging substrates generated by protein engineering • Rational/semi-rational protein engineering reveals constrained mutational variability • Engineered PALs are outperformed by novel ALs of distinct catalytic site signature.

Mutability-Landscape-Guided Engineering of l-Threonine Aldolase Revealing the Prelog Rule in Mediating Diastereoselectivity of C-C Bond Formation

l-threonine aldolase (LTA) catalyzes C-C bond synthesis with moderate diastereoselectivity. In this study, with LTA from Cellulosilyticum sp (CpLTA) as an object, a mutability landscape was first constructed by performing saturation mutagenesis at substrate access tunnel amino acids. The combinatorial active-site saturation test/iterative saturation mutation (CAST/ISM) strategy was then used to tune diastereoselectivity. As a result, the diastereoselectivity of mutant H305L/Y8H/V143R was improved from 37.2 %_syn to 99.4 %_syn . Furthermore, the diastereoselectivity of mutant H305Y/Y8I/W307E was inverted to 97.2 %_anti . Based on insight provided by molecular dynamics simulations and coevolution analysis, the Prelog rule was employed to illustrate the diastereoselectivity regulation mechanism of LTA, holding that the asymmetric formation of the C-C bond was caused by electrons attacking the carbonyl carbon atom of the substrate aldehyde from the re or si face. The study would be useful to expand LTA applications and guide engineering of other C-C bond-forming enzymes.

A growth selection system for the directed evolution of amine-forming or converting enzymes

Fast screening of enzyme variants is crucial for tailoring biocatalysts for the asymmetric synthesis of non-natural chiral chemicals, such as amines. However, most existing screening methods either are limited by the throughput or require specialized equipment. Herein, we report a simple, high-throughput, low-equipment dependent, and generally applicable growth selection system for engineering amine-forming or converting enzymes and apply it to improve biocatalysts belonging to three different enzyme classes. This results in (i) an amine transaminase variant with 110-fold increased specific activity for the asymmetric synthesis of the chiral amine intermediate of Linagliptin; (ii) a 270-fold improved monoamine oxidase to prepare the chiral amine intermediate of Cinacalcet by deracemization; and (iii) an ammonia lyase variant with a 26-fold increased activity in the asymmetric synthesis of a non-natural amino acid. Our growth selection system is adaptable to different enzyme classes, varying levels of enzyme activities, and thus a flexible tool for various stages of an engineering campaign.

Towards a general approach for tailoring the hydrophobic binding site of phenylalanine ammonia-lyases

Unnatural substituted amino acids play an important role as chiral building blocks, especially for pharmaceutical industry, where the synthesis of chiral biologically active molecules still represents an open challenge. Recently, modification of the hydrophobic binding pocket of phenylalanine ammonia-lyase from Petroselinum crispum (PcPAL) resulted in specifically tailored PcPAL variants, contributing to a rational design template for PAL-activity enhancements towards the differently substituted substrate analogues. Within this study we tested the general applicability of this rational design model in case of PALs, of different sources, such as from Arabidopsis thaliana (AtPAL) and Rhodosporidium toruloides (RtPAL). With some exceptions, the results support that the positions of substrate specificity modulating residues are conserved among PALs, thus the mutation with beneficial effect for PAL-activity enhancement can be predicted using the established rational design model. Accordingly, the study supports that tailoring PALs of different origins and different substrate scope, can be performed through a general method. Moreover, the fact that AtPAL variants I461V, L133A and L257V, all outperformed in terms of catalytic efficiency the corresponding, previously reported, highly efficient PcPAL variants, of identical catalytic site, suggests that not only catalytic site differences influence the PAL-activity, thus for the selection of the optimal PAL-biocatalysts for a targeted process, screening of PALs from different origins, should be included.

Integrating dynamics into enzyme engineering

Enzyme engineering has become a widely adopted practice in research labs and industry. In parallel, the past decades have seen tremendous strides in characterizing the dynamics of proteins, using a growing array of methodologies. Importantly, links have been established between the dynamics of proteins and their function. Characterizing the dynamics of an enzyme prior to, and following, its engineering is beginning to inform on the potential of 'dynamic engineering', i.e. the rational modification of protein dynamics to alter enzyme function. Here we examine the state of knowledge at the intersection of enzyme engineering and protein dynamics, describe current challenges and highlight pioneering work in the nascent area of dynamic engineering.

Robust, site-specifically immobilized phenylalanine ammonia-lyases for the enantioselective ammonia addition of cinnamic acids

Phenylalanine ammonia-lyases (PALs) catalyse the non-oxidative deamination of l-phenylalanine to trans-cinnamic acid, while in the presence of high ammonia concentration, the synthetically attractive reverse reaction occurs. Although they have been intensively studied, the wider application of PALs for the large scale synthesis of non-natural amino acids is still rather limited, mainly due to the decreased operational stability of PALs under the high ammonia concentration conditions of ammonia addition. Herein, we describe the development of a highly stable and active immobilized PAL-biocatalyst obtained through site-specific covalent immobilization onto single-walled carbon nanotubes (SWCNTs), employing maleimide/thiol coupling of engineered enzymes containing surficial Cys residues. The immobilization method afforded robust biocatalysts (by strong covalent attachment to the support) and allowed modulation of enzymatic activity (by proper selection of binding site, controlling the orientation of the enzyme attached to the support). The novel biocatalysts were investigated in PAL-catalyzed reactions, focusing on the synthetically challenging ammonia addition reaction. The optimization of the immobilization (enzyme load) and reaction conditions (substrate : biocatalyst ratio, ammonia source, reaction temperature) involving the best performing biocatalyst SWCNTNH2 -SS-PcPAL was performed. The biocatalyst, under the optimal reaction conditions, showed high catalytic efficiency, providing excellent conversion (c ∼90% in 10 h) of cinnamic acid into l-Phe, and more importantly, possesses high operational stability, maintaining its high efficiency over >7 reaction cycles. Moreover, the site-specifically immobilized PcPAL L134A/S614C and PcPAL I460V/S614C variants were successfully applied in the synthesis of several l-phenylalanine analogues of high synthetic value, providing perspectives for the efficient replacement of classical synthetic methods for l-phenylalanines with a mild, selective and eco-friendly enzymatic alternative.