Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

High Cell Density Cultivation Combined with high Specific Enzyme Activity: Cultivation Protocol for the Production of an Amine Transaminase from Bacillus megaterium in E. coli

High cell density cultivation is an established method for the production of various industrially important products such as recombinant proteins. However, these protocols are not always suitable for biocatalytic processes as the focus often lies on biomass production rather than high specific activities of the enzyme inside the cells. In contrast, a range of shake flask protocols are well known with high specific activities but rather low cell densities. To overcome this gap, we established a tailor-made fed-batch protocol combining both aspects: high cell density and high specific activities of heterologously produced enzyme. Using the example of an industrially relevant amine transaminase from Bacillus megaterium, we describe a strategy to optimize the cultivation yield based on the feed rate, IPTG concentration, and post-induction temperature. By adjusting these key parameters, we were able to increase the specific activity by 2.6-fold and the wet cell weight by even 17-fold compared to shake flasks. Finally, we were able to verify our established protocol by transferring it to another experimenter. With that, our optimization strategy can serve as a template for the production of high titers of heterologously produced, active enzymes and might enable the availability of these catalysts for upscaling biocatalytic processes.

Rational design of transaminases based on comparative analysis of catalytically active and distance-free modes of the high-energy intermediate state

Bioproduction of chiral amines is limited by low transaminase (TA) activity on nonnatural substrates, leading to the need for protein engineering. To address the challenge of quickly and precisely identifying the engineering targets, a strategy was proposed based on analyzing the mode changes in the high-energy intermediate state (H-state) of the substrate-enzyme complex during catalysis. By substituting the residues with minimal structural changes in catalytically active mode (A-mode) and distance-free mode (F-mode) of the H-state complex with more conserved ones to stabilize it, a TA mutant M5(T295C/L387A/V436A) with 121.9-fold higher activity for synthesizing the (S)-Rivastigmine precursor (S)-1-(3-methoxyphenyl)ethylamine was created. The applicability of this strategy was also validated by engineering another TA for 1.52-fold higher activity and >99% selectivity toward (R)-3-amino-1-butanol biopreparation. The much higher stereoselectivity of the mutant compared with the wild type (28.3%) demonstrated that this strategy is not only advantageous in engineering enzyme activity but also applicable for modulating stereoselectivity.

1,3-Proton Transfer of Pyridoxal 5'-Phosphate Schiff Base in the Branched-Chain Aminotransferase: Concerted or Stepwise Mechanism?

The branched-chain aminotransferase from Mycobacterium tuberculosis (MtIlvE) is a pyridoxal 5'-phosphate (PLP) dependent enzyme, and it is essential for the synthesis of the branched-chain amino acids. Ketimine is an important intermediate in the catalytic process. We have investigated the mechanism of ketimine formation and the energy landscape using the multiple computational methods. It is found that the 1,3-proton transfer involved in ketimine formation occurs through a stepwise process rather than a one-step process. Lys204 is identified as a key residue for ligand binding and as a base that abstracts the Cα proton from the PLP-Glu Schiff base, yielding a carbanionic intermediate. The first proton transfer is the rate-limiting step with an energy barrier of 17.8 kcal mol-1. Our study disclosed the detailed pathway of the proton transfer from external aldimine to ketimine, providing novel insights into the catalytic mechanism of other PLP-dependent enzymes.

Reprogramming biocatalytic futile cycles through computational engineering of stereochemical promiscuity to create an amine racemase

Repurposing the intrinsic properties of natural enzymes can offer a viable solution to current synthetic challenges through the development of novel biocatalytic processes. Although amino acid racemases are ubiquitous in living organisms, an amine racemase (AR) has not yet been discovered despite its synthetic potential for producing chiral amines. Here, we report the creation of an AR based on the serendipitous discovery that amine transaminases (ATAs) can perform stereoinversion of 2-aminobutane. Kinetic modeling revealed that the unexpected off-pathway activity results from stereochemically promiscuous futile cycles due to incomplete stereoselectivity for 2-aminobutane. This finding motivated us to engineer an S-selective ATA through in silico alanine scanning and empirical combinatorial mutations, creating an AR with broad substrate specificity. The resulting AR, carrying double point mutations, enables the racemization of both enantiomers of diverse chiral amines in the presence of a cognate ketone. This strategy may be generally applicable to a wide range of transaminases, paving the way for the development of new-to-nature racemases.

The Quantum Chemical Cluster Approach in Biocatalysis

ConspectusThe quantum chemical cluster approach has been used for modeling enzyme active sites and reaction mechanisms for more than two decades. In this methodology, a relatively small part of the enzyme around the active site is selected as a model, and quantum chemical methods, typically density functional theory, are used to calculate energies and other properties. The surrounding enzyme is modeled using implicit solvation and atom fixing techniques. Over the years, a large number of enzyme mechanisms have been solved using this method. The models have gradually become larger as a result of the faster computers, and new kinds of questions have been addressed. In this Account, we review how the cluster approach can be utilized in the field of biocatalysis. Examples from our recent work are chosen to illustrate various aspects of the methodology. The use of the cluster model to explore substrate binding is discussed first. It is emphasized that a comprehensive search is necessary in order to identify the lowest-energy binding mode(s). It is also argued that the best binding mode might not be the productive one, and the full reactions for a number of enzyme-substrate complexes have therefore to be considered to find the lowest-energy reaction pathway. Next, examples are given of how the cluster approach can help in the elucidation of detailed reaction mechanisms of biocatalytically interesting enzymes, and how this knowledge can be exploited to develop enzymes with new functions or to understand the reasons for lack of activity toward non-natural substrates. The enzymes discussed in this context are phenolic acid decarboxylase and metal-dependent decarboxylases from the amidohydrolase superfamily. Next, the application of the cluster approach in the investigation of enzymatic enantioselectivity is discussed. The reaction of strictosidine synthase is selected as a case study, where the cluster calculations could reproduce and rationalize the selectivities of both the natural and non-natural substrates. Finally, we discuss how the cluster approach can be used to guide the rational design of enzyme variants with improved activity and selectivity. Acyl transferase from Mycobacterium smegmatis serves as an instructive example here, for which the calculations could pinpoint the factors controlling the reaction specificity and enantioselectivity. The cases discussed in this Account highlight thus the value of the cluster approach as a tool in biocatalysis. It complements experiments and other computational techniques in this field and provides insights that can be used to understand existing enzymes and to develop new variants with tailored properties.

A computational study of the reaction mechanism and stereospecificity of dihydropyrimidinase

Dihydropyrimidinase (DHPase) is a key enzyme in the pyrimidine pathway, the catabolic route for synthesis of β-amino acids. It catalyses the reversible conversion of 5,6-dihydrouracil (DHU) or 5,6-dihydrothymine (DHT) to the corresponding N-carbamoyl-β-amino acids. This enzyme has the potential to be used as a tool in the production of β-amino acids. Here, the reaction mechanism and origin of stereospecificity of DHPases from Saccharomyces kluyveri and Sinorhizobium meliloti CECT4114 were investigated and compared using a quantum mechanical cluster approach based on density functional theory. Two models of the enzyme active site were designed from the X-ray crystal structure of the native enzyme: a small cluster to characterize the mechanism and the stationary points and a large model to probe the stereospecificity and the role of stereo-gate-loop (SGL) residues. It is shown that a hydroxide ion first performs a nucleophilic attack on the substrate, followed by the abstraction of a proton by Asp358, which occurs concertedly with protonation of the ring nitrogen by the same residue. For the DHT substrate, the enzyme displays a preference for the L-configuration, in good agreement with experimental observation. Comparison of the reaction energetics of the two models reveals the importance of SGL residues in the stereospecificity of catalysis. The role of the conserved Tyr172 residue in transition-state stabilization is confirmed as the Tyr172Phe mutation increases the activation barrier of the reaction by ∼8 kcal mol^-1. A detailed understanding of the catalytic mechanism of the enzyme could offer insight for engineering in order to enhance its activity and substrate scope.

Catalytic Mechanism of Pyridoxal 5'-Phosphate-Dependent Aminodeoxychorismate Lyase: A Computational QM/MM Study

Aminodeoxychorismate lyase (ADCL) is a kind of pyridoxal-5'-phosphate (PLP)-dependent enzyme that catalyzes the conversion of 4-amino-4-deoxychorismate (ADC) to p-aminobenzoate (PABA), which is a key step for the biosynthesis of folate. To illuminate the reaction details at the atomistic level, an enzyme-substrate reactant model has been constructed, and QM/MM calculations have been performed. Our calculation results reveal that the overall catalytic cycle contains 11 elementary steps, which can be described by three stages, including the transamination reaction of PLP, the release of pyruvate and aromatization of ADC, and the recovery to the initial aldimine. During the reaction, a series of intramolecular proton transfer are involved, which are the key for the C-N bond formation and cleavage as well as the aromatization of the ADC ring. In addition to forming the Schiff base with the pocket residue Lys251 and substrate in the internal aldimine and the external aldimine, respectively, the coenzyme PLP also plays a critical role in the intramolecular proton transfer by employing its hydroxyl oxygen anion and phosphate group. These findings may provide useful information for further understanding the catalytic mechanism of other PLP-dependent enzymes.

Mechanism-Guided Computational Design of ω-Transaminase by Reprograming of High-Energy-Barrier Steps

ω-Transaminases (ω-TAs) show considerable potential for the synthesis of chiral amines. However, their low catalytic efficiency towards bulky substrates limits their application, and complicated catalytic mechanisms prevent precise enzyme design. Herein, we address this challenge using a mechanism-guided computational enzyme design strategy by reprograming the transition and ground states in key reaction steps. The common features among the three high-energy-barrier steps responsible for the low catalytic efficiency were revealed using quantum mechanics (QM). Five key residues were simultaneously tailored to stabilize the rate-limiting transition state with the aid of the Rosetta design. The 14 top-ranked variants showed 16.9-143-fold improved catalytic activity. The catalytic efficiency of the best variant, M9 (Q25F/M60W/W64F/I266A), was significantly increased, with a 1660-fold increase in k_cat /K_m and a 1.5-26.8-fold increase in turnover number (TON) towards various indanone derivatives.

Construction and yield optimization of a cinnamylamine biosynthesis route in Escherichia coli

Background

With the development of metabolic engineering and synthetic biology, the biosynthesis of aromatic compounds has attracted much attention. Cinnamylamine is an aromatic compound derived from l-phenylalanine, which is used in the synthesis of biologically active molecules, including drugs, and energetic materials. Cinnamylamine has been mainly synthesized by chemical methods to date, and few reports have focused on the biosynthesis of cinnamylamine. Therefore, it is desirable to establish an efficient biosynthesis method for cinnamylamine.

Results

The ω-aminotransferase Cv-ωTA from Chromobacterium violaceum has been demonstrated to have high enzyme activity in the conversion of cinnamaldehyde to cinnamylamine. To prevent the preferable conversion of cinnamaldehyde to cinnamyl alcohol in wild-type Escherichia coli, the E. coli MG1655 strain with reduced aromatic aldehyde reduction (RARE) in which six aldehyde ketone reductase and alcohol dehydrogenase genes have been knocked out was employed. Then, the carboxylic acid reductase from Neurospora crassa (NcCAR) and phosphopantetheinyl transferase (PPTase) from E. coli were screened for a high conversion rate of cinnamic acid to cinnamaldehyde. To shift the equilibrium of the reaction toward cinnamylamine, saturation mutagenesis of Cv-ωTA at key amino acid residues was performed, and Cv-ωTA Y168G had the highest conversion rate with 88.56 mg/L cinnamylamine obtained after 4 h of fermentation. Finally, by optimizing the substrates and the supply of the cofactors, PLP and NADPH, in the fermentation, the yield of cinnamylamine in engineered E. coli reached 523.15 mg/L.

Conclusion

We achieved the first biosynthesis of cinnamylamine using cinnamic acid as the precursor in E. coli using a combinatorial metabolic engineering strategy. This study provides a reference for the biosynthesis of other amine compounds and lays a foundation for the de novo synthesis of cinnamylamine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02199-7.

A DFT study of the active role of the phosphate group of an internal aldimine in a transamination reaction

A transamination reaction from an internal aldimine ([PLP]) and (S)-alanine to pyridoxamine phosphate (PMP) and pyruvic acid was investigated by DFT calculations. As [PLP], a model where the lysine (-Lys) part was approximated by -CH[-NH-C(O)-CH₃]-C(O)-NH-CH₃ was adopted. (H₂O)₄ was also included to trace reaction paths involving proton transfers. 13 elementary processes were obtained. For (the external aldimine → quinoid), (quinoid → ketimine) and (ketimine → carbinol amine) processes, the water dimer was found to connect a phosphate-group oxygen with the moving proton. The connection promoted the Grotthuss-type proton transfer in transition states. It was revealed that the phosphate group is not a mere substituent but has the central role in the transfer.

Predicting Effects of Site-Directed Mutagenesis on Enzyme Kinetics by QM/MM and QM Calculations: A Case of Glutamate Carboxypeptidase II

Quantum and molecular mechanics (QM/MM) and QM-only (cluster model) modeling techniques represent the two workhorses in mechanistic understanding of enzyme catalysis. One of the stringent tests for QM/MM and/or QM approaches is to provide quantitative answers to real-world biochemical questions, such as the effect of single-point mutations on enzyme kinetics. This translates into predicting the relative activation energies to 1-2 kcal·mol^-1 accuracy; such predictions can be used for the rational design of novel enzyme variants with desired/improved characteristics. Herein, we employ glutamate carboxypeptidase II (GCPII), a dizinc metallopeptidase, also known as the prostate specific membrane antigen, as a model system. The structure and activity of this major cancer antigen have been thoroughly studied, both experimentally and computationally, which makes it an ideal model system for method development. Its reaction mechanism is quite well understood: the reaction coordinate comprises a "tetrahedral intermediate" and two transition states and experimental activation Gibbs free energy of ∼17.5 kcal·mol^-1 can be inferred for the known k_cat ≈ 1 s^-1. We correlate experimental kinetic data (including the E424H variant, newly characterized in this work) for various GCPII mutants (k_cat = 8.6 × 10^-5 s^-1 to 2.7 s^-1) with the energy profiles calculated by QM/MM and QM-only (cluster model) approaches. We show that the near-quantitative agreement between the experimental values and the calculated activation energies (ΔH^⧧) can be obtained and recommend the combination of the two protocols: QM/MM optimized structures and cluster model (QM) energetics. The trend in relative activation energies is mostly independent of the QM method (DFT functional) used. Last but not least, a satisfactory correlation between experimental and theoretical data allows us to provide qualitative and fairly simple explanations of the observed kinetic effects which are thus based on a rigorous footing.

Computational Prediction of ω-Transaminase Specificity by a Combination of Docking and Molecular Dynamics Simulations

ω-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 Unexplored Importance of Fleeting Chiral Intermediates in Enzyme-Catalyzed Reactions

Decades of extensive research efforts by biochemists, organic chemists, and protein engineers have led to an understanding of the basic mechanisms of essentially all known types of enzymes, but in a formidable number of cases an essential aspect has been overlooked. The occurrence of short-lived chiral intermediates formed by symmetry-breaking of prochiral precursors in enzyme catalyzed reactions has been systematically neglected. We designate these elusive species as fleeting chiral intermediates and analyze such crucial questions as "Do such intermediates occur in homochiral form?" If so, what is the absolute configuration, and why did Nature choose that particular stereoisomeric form, even when the isolable final product may be achiral? Does the absolute configuration of a chiral product depend in any way on the absolute configuration of the fleeting chiral precursor? How does this affect the catalytic proficiency of the enzyme? If these issues continue to be unexplored, then an understanding of the mechanisms of many enzyme types remains incomplete. We have systematized the occurrence of these chiral intermediates according to their structures and enzyme types. This is followed by critical analyses of selected case studies and by final conclusions and perspectives. We hope that the fascinating concept of fleeting chiral intermediates will attract the attention of scientists, thereby opening an exciting new research field.

Computational Redesign of an ω-Transaminase from Pseudomonas jessenii for Asymmetric Synthesis of Enantiopure Bulky Amines

ω-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.

Structural Evidence of Active Site Adaptability towards Different Sized Substrates of Aromatic Amino Acid Aminotransferase from Psychrobacter Sp. B6

Aromatic amino acid aminotransferases present a special potential in the production of drugs and synthons, thanks to their ability to accommodate a wider range of substrates in their active site, in contrast to aliphatic amino acid aminotransferases. The mechanism of active site adjustment toward substrates of psychrophilic aromatic amino acid aminotransferase (PsyArAT) from Psychrobacter sp. B6 is discussed based on crystal structures of complexes with four hydroxy-analogs of substrates: phenylalanine, tyrosine, tryptophan and aspartic acid. These competitive inhibitors are bound in the active center of PsyArAT but do not undergo transamination reaction, which makes them an outstanding tool for examination of the enzyme catalytic center. The use of hydroxy-acids enabled insight into substrate binding by native PsyArAT, without mutating the catalytic lysine and modifying cofactor interactions. Thus, the binding mode of substrates and the resulting analysis of the volume of the catalytic site is close to a native condition. Observation of these inhibitors' binding allows for explanation of the enzyme's adaptability to process various sizes of substrates and to gain knowledge about its potential biotechnological application. Depending on the character and size of the used inhibitors, the enzyme crystallized in different space groups and showed conformational changes of the active site upon ligand binding.

Theoretical Study on the Mechanism of the Acylate Reaction of β-Lactamase

Using density functional theory and a cluster approach, we study the reaction potential surface and compute Gibbs free energies for the acylate reaction of β-lactamase with penicillin G, where the solvent effect is important and taken into consideration. Two reaction paths are investigated: one is a multi-step process with a rate-limit energy barrier of 19.1 kcal/mol, which is relatively small, and the reaction can easily occur; the other is a one-step process with a barrier of 45.0 kcal/mol, which is large and thus makes the reaction hard to occur. The reason why the two paths have different barriers is explained.

Efficient Amino Donor Recycling in Amination Reactions: Development of a New Alanine Dehydrogenase in Continuous Flow and Dialysis Membrane Reactors

Transaminases have arisen as one of the main biocatalysts for amine production but despite their many advantages, their stability is still a concern for widespread application. One of the reasons for their instability is the need to use an excess of the amino donor when trying to synthesise amines with unfavourable equilibria. To circumvent this, recycling systems for the amino donor, such as amino acid dehydrogenases or aldolases, have proved useful to push the equilibria while avoiding high amino donor concentrations. In this work, we report the use of a new alanine dehydrogenase from the halotolerant bacteria Halomonas elongata which exhibits excellent stability to different cosolvents, combined with the well characterised CbFDH as a recycling system of L-alanine for the amination of three model substrates with unfavourable equilibria. In a step forward, the amino donor recycling system has been co-immobilised and used in flow with success as well as re-used as a dialysis enclosed system for the amination of an aromatic aldehyde.

Engineering the large pocket of an (S)-selective transaminase for asymmetric synthesis of (S)-1-amino-1-phenylpropane

Amine transaminases offer an environmentally benign chiral amine asymmetric synthesis route.

A QM/MM simulation study of transamination reaction at the active site of aspartate aminotransferase: Free energy landscape and proton transfer pathways

Transaminase is a key enzyme for amino acid metabolism, which reversibly catalyzes the transamination reaction with the help of PLP (pyridoxal 5^' -phosphate) as its cofactor. Here we have investigated the mechanism and free energy landscape of the transamination reaction involving the aspartate transaminase (AspTase) enzyme and aspartate-PLP (Asp-PLP) complex using QM/MM simulation and metadynamics methods. The reaction is found to follow a stepwise mechanism where the active site residue Lys258 acts as a base to shuttle a proton from α-carbon (CA) to imine carbon (C4A) of the PLP-Asp Schiff base. In the first step, the Lys258 abstracts the CA proton of the substrate leading to the formation of a carbanionic intermediate which is followed by the reprotonation of the Asp-PLP Schiff base at C4A atom by Lys258. It is found that the free energy barrier for the proton abstraction by Lys258 and that for the reprotonation are 17.85 and 3.57 kcal/mol, respectively. The carbanionic intermediate is 7.14 kcal/mol higher in energy than the reactant. Hence, the first step acts as the rate limiting step. The present calculations also show that the Lys258 residue undergoes a conformational change after the first step of transamination reaction and becomes proximal to C4A atom of the Asp-PLP Schiff base to favor the second step. The active site residues Tyr70* and Gly38 anchor the Lys258 in proper position and orientation during the first step of the reaction and stabilize the positive charge over Lys258 generated at the intermediate step.

Asymmetric Synthesis of Optically Pure Aliphatic Amines with an Engineered Robust ω-Transaminase

The production of chiral amines by transaminase-catalyzed amination of ketones is an important application of biocatalysis in synthetic chemistry. It requires transaminases that show high enantioselectivity in asymmetric conversion of the ketone precursors. A robust derivative of ω-transaminase from Pseudomonasjessenii (PjTA-R6) that naturally acts on aliphatic substrates was constructed previously by our group. Here, we explore the catalytic potential of this thermostable enzyme for the synthesis of optically pure aliphatic amines and compare it to the well-studied transaminases from Vibrio fluvialis (VfTA) and Chromobacterium violaceum (CvTA). The product yields indicated improved performance of PjTA-R6 over the other transaminases, and in most cases, the optical purity of the produced amine was above 99% enantiomeric excess (e.e.). Structural analysis revealed that the substrate binding poses were influenced and restricted by the switching arginine and that this accounted for differences in substrate specificities. Rosetta docking calculations with external aldimine structures showed a correlation between docking scores and synthetic yields. The results show that PjTA-R6 is a promising biocatalyst for the asymmetric synthesis of aliphatic amines with a product spectrum that can be explained by its structural features.

Structural studies reveal flexible roof of active site responsible for ω-transaminase CrmG overcoming by-product inhibition

Amine compounds biosynthesis using ω-transaminases has received considerable attention in the pharmaceutical industry. However, the application of ω-transaminases was hampered by the fundamental challenge of severe by-product inhibition. Here, we report that ω-transaminase CrmG from Actinoalloteichus cyanogriseus WH1-2216-6 is insensitive to inhibition from by-product α-ketoglutarate or pyruvate. Combined with structural and QM/MM studies, we establish the detailed catalytic mechanism for CrmG. Our structural and biochemical studies reveal that the roof of the active site in PMP-bound CrmG is flexible, which will facilitate the PMP or by-product to dissociate from PMP-bound CrmG. Our results also show that amino acceptor caerulomycin M (CRM M), but not α-ketoglutarate or pyruvate, can form strong interactions with the roof of the active site in PMP-bound CrmG. Based on our results, we propose that the flexible roof of the active site in PMP-bound CrmG may facilitate CrmG to overcome inhibition from the by-product.

A study on the structure, mechanism, and biochemistry of kanamycin B dioxygenase (KanJ)-an enzyme with a broad range of substrates

Kanamycin A is an aminoglycoside antibiotic isolated from Streptomyces kanamyceticus and used against a wide spectrum of bacteria, including Mycobacterium tuberculosis. Biosynthesis of kanamycin involves an oxidative deamination step catalyzed by kanamycin B dioxygenase (KanJ), thereby the C2' position of kanamycin B is transformed into a keto group upon release of ammonia. Here, we present for the first time, structural models of KanJ with several ligands, which along with the results of ITC binding assays and HPLC activity tests explain substrate specificity of the enzyme. The large size of the binding pocket suggests that KanJ can accept a broad range of substrates, which was confirmed by activity tests. Specificity of the enzyme with respect to its substrate is determined by the hydrogen bond interactions between the methylamino group of the antibiotic and highly conserved Asp134 and Cys150 as well as between hydroxyl groups of the substrate and Asn120 and Gln80. Upon antibiotic binding, the C terminus loop is significantly rearranged and Gln80 and Asn120, which are directly involved in substrate recognition, change their conformations. Based on reaction energy profiles obtained by density functional theory (DFT) simulations, we propose a mechanism of ketone formation involving the reactive Fe^IV  = O and proceeding either via OH rebound, which yields a hemiaminal intermediate or by abstraction of two hydrogen atoms, which leads to an imine species. At acidic pH, the latter involves a lower barrier than the OH rebound, whereas at basic pH, the barrier leading to an imine vanishes completely. DATABASES: Structural data are available in PDB database under the accession numbers: 6S0R, 6S0T, 6S0U, 6S0W, 6S0V, 6S0S. Diffraction images are available at the Integrated Resource for Reproducibility in Macromolecular Crystallography at http://proteindiffraction.org under DOIs: 10.18430/m36s0t, 10.18430/m36s0u, 10.18430/m36s0r, 10.18430/m36s0s, 10.18430/m36s0v, 10.18430/m36s0w. A data set collection of computational results is available in the Mendeley Data database under DOI: 10.17632/sbyzssjmp3.1 and in the ioChem-BD database under DOI: 10.19061/iochem-bd-4-18.

Protein-ligand binding with the coarse-grained Martini model

The detailed understanding of the binding of small molecules to proteins is the key for the development of novel drugs or to increase the acceptance of substrates by enzymes. Nowadays, computer-aided design of protein–ligand binding is an important tool to accomplish this task. Current approaches typically rely on high-throughput docking essays or computationally expensive atomistic molecular dynamics simulations. Here, we present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein–ligand interactions of small drug-like molecules. Remarkably, we achieve high accuracy without the need of any a priori knowledge of binding pockets or pathways. Our approach is applied to a range of systems from the well-characterized T4 lysozyme over members of the GPCR family and nuclear receptors to a variety of enzymes. The presented results open the way to high-throughput screening of ligand libraries or protein mutations using the coarse-grained Martini model.

Computer-aided design of protein-ligand binding is important for the development of novel drugs. Here authors present an approach to use the recently re-parametrized coarse-grained Martini model to perform unbiased millisecond sampling of protein-ligand binding interactions of small drug-like molecules.

Using Atomic Confining Potentials for Geometry Optimization and Vibrational Frequency Calculations in Quantum-Chemical Models of Enzyme Active Sites

Quantum-chemical studies of enzymatic reaction mechanisms sometimes use truncated active-site models as simplified alternatives to mixed quantum mechanics molecular mechanics (QM/MM) procedures. Eliminating the MM degrees of freedom reduces the complexity of the sampling problem, but the trade-off is the need to introduce geometric constraints in order to prevent structural collapse of the model system during geometry optimizations that do not contain a full protein backbone. These constraints may impair the efficiency of the optimization, and care must be taken to avoid artifacts such as imaginary vibrational frequencies. We introduce a simple alternative in which terminal atoms of the model system are placed in soft harmonic confining potentials rather than being rigidly constrained. This modification is simple to implement and straightforward to use in vibrational frequency calculations, unlike iterative constraint-satisfaction algorithms, and allows the optimization to proceed without constraint even though the practical result is to fix the anchor atoms in space. The new approach is more efficient for optimizing minima and transition states, as compared to the use of fixed-atom constraints, and also more robust against unwanted imaginary frequencies. We illustrate the method by application to several enzymatic reaction pathways where entropy makes a significant contribution to the relevant reaction barriers. The use of confining potentials correctly describes reaction paths and facilitates calculation of both vibrational zero-point and finite-temperature entropic corrections to barrier heights.

On the reaction mechanism of an endoperoxide ring formation by fumitremorgin B endoperoxidase. The right arrangement makes a difference

Fumitremorgin B endoperoxidase (FtmF) belongs to 2-oxoglutarate dependent dioxygenases and catalyzes an unusual oxidative reaction of endoperoxide formation at the final stage of biosynthesis of verruculogen - a mycotoxin produced by Aspergillus and Penicillinum strains. The published crystal structure of FtmF (PDB: ), which is of overall good quality, contains a model of the substrate bound in the active site, which, however, has very low occupancy and its conformation does not comply with the small molecule crystal structure. Moreover, a previous computational study that employed a model based on this crystal structure revealed a substantial reaction barrier, which might indicate that the model of FtmF/substrate complex can have serious errors. The purpose of this work was to model with computational methods the structure of the enzyme-substrate complex and to investigate the mechanisms of the enzymatic reaction. Docking, molecular dynamics simulation and DFT results, all indicate the substrate most likely binds in the active site in a configuration very different from that originally suggested. Moreover, for this newly proposed structure of the enzyme-substrate complex, the reaction energy profile is characterised exclusively by low barriers and it successfully explains the observed regiospecificity of the enzymatic process. Finally, a plausible binding site for ascorbate was found and it is suggested that ascorbate is involved in the final step of the FtmF reaction.

PLP-dependent enzymes as important biocatalysts for the pharmaceutical, chemical and food industries: a structural and mechanistic perspective

PLP-dependent enzymes described on this review are attractive targets for enzyme engineering towards their application in an industrial biotechnology framework.

Computer-Guided Surface Engineering for Enzyme Improvement

Protein engineering strategies are often guided by our understanding of how the structure of a protein determines its function. However, our understanding is generally restricted to small regions of a protein, namely the active site and its immediate vicinity, while the remainder of the protein is something of an enigma. Studying highly homologous transaminases with strictly conserved active sites, but different substrate preferences and activities, we predict and experimentally validate that the surface of the protein far from the active site carries out a decisive role in substrate selectivity and catalytic efficiency. Using a unique molecular dynamics approach and novel trajectory analysis, we demonstrate the phenomenon of surface-directed ligand diffusion in this well-known protein family for the first time. Further, we identify the residues involved in directing substrate, design surface channel variants endowed for improved kinetic properties and establish a broadly applicable new approach for protein engineering.

A Theoretical Study of the Benzoylformate Decarboxylase Reaction Mechanism

Density functional theory calculations are used to investigate the detailed reaction mechanism of benzoylformate decarboxylase, a thiamin diphosphate (ThDP)-dependent enzyme that catalyzes the nonoxidative decarboxylation of benzoylformate yielding benzaldehyde and carbon dioxide. A large model of the active site is constructed on the basis of the X-ray structure, and it is used to characterize the involved intermediates and transition states and evaluate their energies. There is generally good agreement between the calculations and available experimental data. The roles of the various active site residues are discussed and the results are compared to mutagenesis experiments. Importantly, the calculations identify off-cycle intermediate species of the ThDP cofactor that can have implications on the kinetics of the reaction.

Stereoelectronic effects in the reaction of aromatic substrates catalysed by Halomonas elongata transaminase and its mutants

A transaminase from Halomonas elongata and four mutants generated by an in silico-based design were recombinantly produced in E. coli, purified and applied to the amination of mono-substituted aromatic carbonyl-derivatives. While benzaldehyde derivatives were excellent substrates, only NO₂-acetophenones were transformed into the (S)-amine with a high enantioselectivity. The different behaviour of wild-type and mutated transaminases was assessed by in silico substrate binding mode studies.

Application of ω-Transaminases in the Pharmaceutical Industry

Chiral amines are valuable building blocks for the pharmaceutical industry. ω-TAms have emerged as an exciting option for their synthesis, offering a potential "green alternative" to overcome the drawbacks associated with conventional chemical methods. In this review, we explore the application of ω-TAms for pharmaceutical production. We discuss the diverse array of reactions available involving ω-TAms and process considerations of their use in both kinetic resolution and asymmetric synthesis. With the aid of specific drug intermediates and APIs, we chart the development of ω-TAms using protein engineering and their contribution to elegant one-pot cascades with other enzymes, including carbonyl reductases (CREDs), hydrolases and monoamine oxidases (MAOs), providing a comprehensive overview of their uses, beginning with initial applications through to the present day.

Recent Trends in Quantum Chemical Modeling of Enzymatic Reactions

The quantum chemical cluster approach is a powerful method for investigating enzymatic reactions. Over the past two decades, a large number of highly diverse systems have been studied and a great wealth of mechanistic insight has been developed using this technique. This Perspective reviews the current status of the methodology. The latest technical developments are highlighted, and challenges are discussed. Some recent applications are presented to illustrate the capabilities and progress of this approach, and likely future directions are outlined.

Quantum Chemical Study of Dual-Substrate Recognition in ω-Transaminase

ω-Transaminases are attractive biocatalysts for the production of chiral amines. These enzymes usually have a broad substrate range. Their substrates include hydrophobic amines as well as amino acids, a feature referred to as dual-substrate recognition. In the present study, the reaction mechanism for the half-transamination of l-alanine to pyruvate in (S)-selective Chromobacterium violaceum ω-transaminase is investigated using density functional theory calculations. The role of a flexible arginine residue, Arg416, in the dual-substrate recognition is investigated by employing two active-site models, one including this residue and one lacking it. The results of this study are compared to those of the mechanism of the conversion of (S)-1-phenylethylamine to acetophenone. The calculations suggest that the deaminations of amino acids and hydrophobic amines follow essentially the same mechanism, but the energetics of the reactions differ significantly. It is shown that the amine is kinetically favored in the half-transamination of l-alanine/pyruvate, whereas the ketone is kinetically favored in the half-transamination of (S)-1-phenylethylamine/acetophenone. The calculations further support the proposal that the arginine residue facilitates the dual-substrate recognition by functioning as an arginine switch, where the side chain is positioned inside or outside of the active site depending on the substrate. Arg416 participates in the binding of l-alanine by forming a salt bridge to the carboxylate moiety, whereas the conversion of (S)-1-phenylethylamine is feasible in the absence of Arg416, which here represents the case in which the side chain of Arg416 is positioned outside of the active site.

Development of in situ product removal strategies in biocatalysis applying scaled-down unit operations

An experimental platform based on scaled-down unit operations combined in a plug-and-play manner enables easy and highly flexible testing of advanced biocatalytic process options such as in situ product removal (ISPR) process strategies. In such a platform, it is possible to compartmentalize different process steps while operating it as a combined system, giving the possibility to test and characterize the performance of novel process concepts and biocatalysts with minimal influence of inhibitory products. Here the capabilities of performing process development by applying scaled-down unit operations are highlighted through a case study investigating the asymmetric synthesis of 1-methyl-3-phenylpropylamine (MPPA) using ω-transaminase, an enzyme in the sub-family of amino transferases (ATAs). An on-line HPLC system was applied to avoid manual sample handling and to semi-automatically characterize ω-transaminases in a scaled-down packed-bed reactor (PBR) module, showing MPPA as a strong inhibitor. To overcome the inhibition, a two-step liquid-liquid extraction (LLE) ISPR concept was tested using scaled-down unit operations combined in a plug-and-play manner. Through the tested ISPR concept, it was possible to continuously feed the main substrate benzylacetone (BA) and extract the main product MPPA throughout the reaction, thereby overcoming the challenges of low substrate solubility and product inhibition. The tested ISPR concept achieved a product concentration of 26.5 g_MPPA  · L^-1 , a purity up to 70% g_MPPA  · g_tot^-1 and a recovery in the range of 80% mol · mol^-1 of MPPA in 20 h, with the possibility to increase the concentration, purity, and recovery further. Biotechnol. Bioeng. 2017;114: 600-609. © 2016 Wiley Periodicals, Inc.

Identification of (S)-selective transaminases for the asymmetric synthesis of bulky chiral amines

The use of transaminases to access pharmaceutically relevant chiral amines is an attractive alternative to transition-metal-catalysed asymmetric chemical synthesis. However, one major challenge is their limited substrate scope. Here we report the creation of highly active and stereoselective transaminases starting from fold class I. The transaminases were developed by extensive protein engineering followed by optimization of the identified motif. The resulting enzymes exhibited up to 8,900-fold higher activity than the starting scaffold and are highly stereoselective (up to >99.9% enantiomeric excess) in the asymmetric synthesis of a set of chiral amines bearing bulky substituents. These enzymes should therefore be suitable for use in the synthesis of a wide array of potential intermediates for pharmaceuticals. We also show that the motif can be engineered into other protein scaffolds with sequence identities as low as 70%, and as such should have a broad impact in the field of biocatalytic synthesis and enzyme engineering.

The Industrial Age of Biocatalytic Transamination

During the last decade the use of ω-transaminases has been identified as a very powerful method for the preparation of optically pure amines from the corresponding ketones. Their immense potential for the preparation of chiral amines, together with their ease of use in combination with existing biocatalytic methods, have made these biocatalysts a competitor to any chemical methodology for (asymmetric) amination. An increasing number of examples, especially from industry, shows that this biocatalytic technology outmaneuvers existing chemical processes by its simple and flexible nature. In the last few years numerous publications and patents on synthetic routes, mainly to pharmaceuticals, involving ω-transaminases have been published. The review gives an overview of the application of ω-transaminases in organic synthesis with a focus on active pharmaceutical ingredients (APIs) and the developments during the last few years.