The CoA ligase from Metallosphaera sedula (MsACL) can be used for the chemoenzymatic synthesis of amides from carboxylic acids. In this CoA-independent conversion, the enzyme catalyzes the adenylation of a carboxylic acid with the help of ATP, followed by the uncatalyzed cleavage of acyl-AMP by a nucleophilic amine to yield an amide. With ω-amino acids as substrates this reaction may result in formation of lactams, but unfortunately the substrate preference of the wild-type enzyme is rather limited. To allow structure-based protein engineering and expand the substrate scope of the enzyme, crystal structures of MsACL were solved in the thioesterification conformational state with AMP, CoA and with the reaction intermediate acetyl-AMP bound in the active site. Using substrate docking and by comparing the crystals structures and sequence of MsACL to those of related CoA ligases, mutations were predicted which increase the affinity in the carboxylic acid binding pocket for ω-amino acids. The resulting mutations transformed a non-active enzyme into an active enzyme for ε-caprolactam synthesis, highlighting the potential of the thermophilic CoA ligase for this synthetic and biotechnologically relevant reaction.

Aspartate ammonia lyases catalyze the reversible amination of fumarate to l-aspartate. Recent studies demonstrate that the thermostable enzyme from Bacillus sp. YM55-1 (AspB) can be engineered for the enantioselective production of substituted beta-amino acids. This reaction would be attractive for the conversion of acrylic acid to beta-alanine, which is an important building block for the preparation of bioactive compounds. Here we describe a bioinformatics and computational approach aimed at introducing the beta-alanine synthesis activity. Three strategies were used: First, we redesigned the alpha-carboxylate binding pocket of AspB to introduce activity with the acrylic acid. Next, different template enzymes were identified by genome mining, equipped with a redesigned alpha-carboxylate pocket, and investigated for beta-alanine synthesis, which yielded variants with better activity. Third, interactions of the SS-loop that covers the active site and harbors a catalytic serine were computationally redesigned using energy calculations to stabilize reactive conformations and thereby further increase the desired beta-alanine synthesis activity. Different improved enzymes were obtained and the best variants showed k cat values with acrylic acid of at least 0.6-1.5 s-1 with K M values in the high mM range. Since the beta-alanine production of wild-type enzyme was below the detection limit, this suggests that the k cat/K m was improved by at least 1000-fold. Crystal structures of the 6-fold mutant of redesigned AspB and the similarly engineered aspartase from Caenibacillus caldisaponilyticus revealed that their ligand-free structures have the SS-loop in a closed (reactive) conformation, which for wild-type AspB is only observed in the substrate-bound enzyme. AlphaFold-generated models suggest that other aspartase variants redesigned for acrylic acid hydroamination also prefer a 3D structure with the loop in a closed conformation. The combination of binding pocket redesign, genome mining, and enhanced active-site loop closure thus created effective beta-alanine synthesizing variants of aspartase.

The P450 monooxygenase CYP109A2 from Bacillus megaterium DSM319 was previously found to convert vitamin D3 (VD3) to 25‐hydroxyvitamin D3. Here, we show that this enzyme is also able to convert testosterone in a highly regio‐ and stereoselective manner to 16β‐hydroxytestosterone. To reveal the structural determinants governing the regio‐ and stereoselective steroid hydroxylation reactions catalyzed by CYP109A2, two crystal structures of CYP109A2 were solved in similar closed conformations, one revealing a bound testosterone in the active site pocket, albeit at a nonproductive site away from the heme‐iron. To examine whether the closed crystal structures nevertheless correspond to a reactive conformation of CYP109A2, docking and molecular dynamics (MD) simulations were performed with testosterone and vitamin D3 (VD3) present in the active site. These MD simulations were analyzed for catalytically productive conformations, the relative occurrences of which were in agreement with the experimentally determined stereoselectivities if the predicted stability of each carbon‐hydrogen bond was taken into account. Overall, the first‐time determination and analysis of the catalytically relevant 3D conformation of CYP109A2 will allow for future small molecule ligand screening in silico, as well as enabling site‐directed mutagenesis toward improved enzymatic properties of this enzyme.

Whereas directed evolution and rational design by structural inspection are established tools for enzyme redesign, computational methods are less mature but have the potential to predict small sets of mutants with desired properties without laboratory screening of large libraries. We have explored the use of computational enzyme redesign to change the enantioselectivity of a highly thermostable alcohol dehydrogenase from Thermus thermophilus in the asymmetric reduction of ketones. The enzyme reduces acetophenone to (S)-1-phenylethanol. To invert the enantioselectivity, we used an adapted CASCO workflow which included Rosetta for enzyme design and molecular dynamics simulations for ranking. To correct for unrealistic binding modes, we used Boltzmann weighing of binding energies computed by a linear interaction energy approach. This computationally cheap method predicted four variants with inverted enantioselectivity, each with 6-8 mutations around the substrate-binding site, causing only modest reduction (2- to 7-fold) of k(cat)/K-M values. Laboratory testing showed that three variants indeed had inverted enantioselectivity, producing (R)-alcohols with up to 99 % enantiomeric excess. The broad substrate range allowed reduction of acetophenone derivatives with full conversion to highly enantioenriched alcohols. The results demonstrate the use of computational methods to control ketoreductase stereoselectivity in asymmetric transformations with minimal experimental screening.

Protein engineering is a powerful and widely applied tool for tailoring enzyme properties to meet application-specific requirements. An attractive group of biocatalysts are PLP-dependent amine transaminases which are capable of converting prochiral ketones to the corresponding chiral amines by asymmetric catalysis. The enzymes often display high enantioselectivity and accept various amine donors. Practical applications of these amine transaminases can be hampered by enzyme instability and by their limited substrate scope. Various strategies to improve robustness of amine transaminases and to redirect their substrate specificity have been explored, including directed evolution, rational design and computation-supported engineering. The approaches used and results obtained are reviewed in this paper, showing that different strategies can be used in a complementary manner and can expand the applicability of amine transaminases in biocatalysis.

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

Caprolactamase is the first enzyme in the caprolactam degradation pathway of Pseudomonas jessenii. It is composed of two subunits (CapA and CapB) and sequence-related to other ATP-dependent enzymes involved in lactam hydrolysis, like 5-oxoprolinases and hydantoinases. Low sequence similarity also exists with ATP-dependent acetone- and acetophenone carboxylases. The caprolactamase was produced in Escherichia coli, isolated by His-tag affinity chromatography, and subjected to functional and structural studies. Activity toward caprolactam required ATP and was dependent on the presence of bicarbonate in the assay buffer. The hydrolysis product was identified as 6-aminocaproic acid. Quantum mechanical modeling indicated that the hydrolysis of caprolactam was highly disfavored (ΔG0 '= 23 kJ/mol), which explained the ATP dependence. A crystal structure showed that the enzyme exists as an (αβ)2 tetramer and revealed an ATP-binding site in CapA and a Zn-coordinating site in CapB. Mutations in the ATP-binding site of CapA (D11A and D295A) significantly reduced product formation. Mutants with substitutions in the metal binding site of CapB (D41A, H99A, D101A, and H124A) were inactive and less thermostable than the wild-type enzyme. These residues proved to be essential for activity and on basis of the experimental findings we propose possible mechanisms for ATP-dependent lactam hydrolysis.

Diaminopimelate decarboxylases (DAPDCs) are highly selective enzymes that catalyze the common final step in different lysine biosynthetic pathways, i.e. the conversion of meso-diaminopimelate (DAP) to L-lysine. We examined the modification of the substrate specificity of the thermostable decarboxylase from Thermotoga maritima with the aim to introduce activity with 2-aminopimelic acid (2-APA) since its decarboxylation leads to 6-aminocaproic acid (6-ACA), a building block for the synthesis of nylon-6. Structure-based mutagenesis of the distal carboxylate binding site resulted in a set of enzyme variants with new activities toward different D-amino acids. One of the mutants (E315T) had lost most of its activity toward DAP and primarily acted as a 2-APA decarboxylase. We next used computational modeling to explain the observed shift in catalytic activities of the mutants. The results suggest that predictive computational protocols can support the redesign of the catalytic properties of this class of decarboxylating PLP-dependent enzymes.

ω-Transaminases (ω-TA) are attractive biocatalysts for the production of chiral amines from prochiral ketones via asymmetric synthesis. However, the substrate scope of ω-TAs is usually limited due to steric hindrance at the active site pockets. We explored a protein engineering strategy using computational design to expand the substrate scope of an (S)-selective ω-TA from Pseudomonas jessenii (PjTA-R6) toward the production of bulky amines. PjTA-R6 is attractive for use in applied biocatalysis due to its thermostability, tolerance to organic solvents, and acceptance of high concentrations of isopropylamine as amino donor. PjTA-R6 showed no detectable activity for the synthesis of six bicyclic or bulky amines targeted in this study. Six small libraries composed of 7–18 variants each were separately designed via computational methods and tested in the laboratory for ketone to amine conversion. In each library, the vast majority of the variants displayed the desired activity, and of the 40 different designs, 38 produced the target amine in good yield with >99% enantiomeric excess. This shows that the substrate scope and enantioselectivity of PjTA mutants could be predicted in silico with high accuracy. The single mutant W58G showed the best performance in the synthesis of five structurally similar bulky amines containing the indan and tetralin moieties. The best variant for the other bulky amine, 1-phenylbutylamine, was the triple mutant W58M + F86L + R417L, indicating that Trp58 is a key residue in the large binding pocket for PjTA-R6 redesign. Crystal structures of the two best variants confirmed the computationally predicted structures. The results show that computational design can be an efficient approach to rapidly expand the substrate scope of ω-TAs to produce enantiopure bulky amines.

Omniligase-1 is a broadly applicable enzyme for peptide bond formation between an activated acyl donor peptide and a non-protected acyl acceptor peptide. The enzyme is derived from an earlier subtilisin variant called peptiligase by several rounds of protein engineering aimed at increasing synthetic yields and substrate range. To examine the contribution of individual mutations on S/H ratio and substrate scope in peptide synthesis, we selected peptiligase variant M222P/L217H as a starting enzyme and introduced successive mutations. Mutation A225N in the S1' pocket and F189W of the S2' pocket increased the synthesis to hydrolysis (S/H) ratio and overall coupling efficiency, whereas the I107V mutation was added to S4 pocket to increase the reaction rate. The final omniligase variants appeared to have a very broad substrate range, coupling more than 250 peptides in a 400-member library of acyl acceptors, as indicated by a high-throughput FRET assay. Crystal structures and computational modelling could rationalize the exceptional properties of omniligase-1 in peptide synthesis.

Protein engineering has shown widespread use in improving the industrial application of enzymes and broadening the conditions they are able to operate under by increasing their thermostability and solvent tolerance. Here, we show that protein engineering can be used to increase the thermostability of an artificial metalloenzyme. Thermostable variants of the human steroid carrier protein 2L, modified to bind a metal catalyst, were created by rational design using structural data and a 3DM database. These variants were tested to identify mutations that enhanced the stability of the protein scaffold, and a significant increase in melting temperature was observed with a number of modified metalloenzymes. The ability to withstand higher reaction temperatures resulted in an increased activity in the hydroformylation of 1-octene, with more than fivefold improvement in turnover number, whereas the selectivity for linear aldehyde remained high up to 80%.

ω-Transaminases (ω-TAs) catalyze the conversion of ketones to chiral amines, often with high enantioselectivity and specificity, which makes them attractive for industrial production of chiral amines. Tailoring ω-TAs to accept non-natural substrates is necessary because of their limited substrate range. We present a computational protocol for predicting the enantioselectivity and catalytic selectivity of an ω-TA from Vibrio fluvialis with different substrates and benchmark it against 62 compounds gathered from the literature. Rosetta-generated complexes containing an external aldimine intermediate of the transamination reaction are used as starting conformations for multiple short independent molecular dynamics (MD) simulations. The combination of molecular docking and MD simulations ensures sufficient and accurate sampling of the relevant conformational space. Based on the frequency of near-attack conformations observed during the MD trajectories, enantioselectivities can be quantitatively predicted. The predicted enantioselectivities are in agreement with a benchmark dataset of experimentally determined ee% values. The substrate-range predictions can be based on the docking score of the external aldimine intermediate. The low computational cost required to run the presented framework makes it feasible for use in enzyme design to screen thousands of enzyme variants.

The use of enzymes in preparative biocatalysis often requires tailoring enzyme selectivity by protein engineering. Herein we explore the use of computational library design and molecular dynamics simulations to create variants of limonene epoxide hydrolase that produce enantiomeric diols from meso-epoxides. Three substrates of different sizes were targeted: cis-2,3-butene oxide, cyclopentene oxide, and cis-stilbene oxide. Most of the 28 designs tested were active and showed the predicted enantioselectivity. Excellent enantioselectivities were obtained for the bulky substrate cis-stilbene oxide, and enantiocomplementary mutants produced (S,S)- and (R,R)-stilbene diol with >97 % enantiomeric excess. An (R,R)-selective mutant was used to prepare (R,R)-stilbene diol with high enantiopurity (98 % conversion into diol, >99 % ee). Some variants displayed higher catalytic rates (kcat ) than the original enzyme, but in most cases KM values increased as well. The results demonstrate the feasibility of computational design and screening to engineer enantioselective epoxide hydrolase variants with very limited laboratory screening.

CYP154C5 from Nocardia farcinica is a P450 monooxygenase able to hydroxylate a range of steroids with high regio- and stereoselectivity at the 16α-position. Using protein and substrate engineering based on the crystal structure of CYP154C5, an altered regioselectivity of the enzyme in steroid hydroxylation could be achieved. Thus, conversion of progesterone by mutant CYP154C5 F92A resulted in formation of the corresponding 21-hydroxylated product 11-deoxycorticosterone in addition to 16α-hydroxylation. Using MD simulation, this altered regioselectivity appeared to result from an alternate binding mode of the steroid in the active site of mutant F92A. MD simulation further suggested that water entrance to the active site caused higher uncoupling in this mutant. Moreover, exclusive 15α-hydroxylation was observed for wild-type CYP154C5 in the conversion of 5α-androstan-3-one, lacking an oxy-functional group at C17. Overall, our data give valuable insight into the structure-function relationship of this cytochrome P450 monooxygenase for steroid hydroxylation.