Is Promiscuous CALB a Good Scaffold for Designing New Epoxidases?

Mechanistic studies of a lipase unveil effect of pH on hydrolysis products of small PET modules

Biocatalysis is a key technology enabling plastic recycling. However, despite advances done in the development of plastic-degrading enzymes, the molecular mechanisms that govern their catalytic performance are poorly understood, hampering the engineering of more efficient enzyme-based technologies. In this work, we study the hydrolysis of PET-derived diesters and PET trimers catalyzed by the highly promiscuous lipase B from Candida antarctica (CALB) through QM/MM molecular dynamics simulations supported by experimental Michaelis-Menten kinetics. The computational studies reveal the role of the pH on the CALB regioselectivity toward the hydrolysis of bis-(hydroxyethyl) terephthalate (BHET). We exploit this insight to perform a pH-controlled biotransformation that selectively hydrolyzes BHET to either its corresponding diacid or monoesters using both soluble and immobilized CALB. The discoveries presented here can be exploited for the valorization of BHET resulting from the organocatalytic depolymerization of PET.

Computational Studies Suggest Promiscuous Candida antarctica Lipase B as an Environmentally Friendly Alternative for the Production of Epoxides

Environmentally friendly processes are nowadays a trending topic to get highly desired chemical compounds and, in this sense, the use of enzyme-catalyzed routes is becoming a promising alternative to traditional synthetic methods. In the present paper, a hybrid quantum mechanics/molecular mechanics (QM/MM) computational study on the epoxidation of alkenes catalyzed by the Ser105Ala variant of the promiscuous Candida antarctica lipase B (CALB) is presented in an attempt to search for alternative paths to get useful intermediates in industries. The catalyzed reaction, described at the atomistic level with a model of the full solvated in a box of water molecules, is compared with the alternative epoxidation of alkenes by peroxy acids in chloroform. Free-energy profiles obtained at the density functional theory (DFT)/MM level show how Ser105Ala CALB is capable of epoxide short alkenes in a two-step process with free-energy barriers, in agreement with available experimental data, that are significantly lower than those of the single-step reaction in solution. The possible (R)-enantioselectivity dictated by the binding step, explored by means of alchemical QM/MM free-energy perturbation (FEP) methods, and the preference for the (S)-enantiomer derived from the free-energy landscape of the chemical steps would cancel out, thus predicting the lack of enantioselectivity experimentally observed. In general, our results provide general information on the molecular mechanism employed by a highly promiscuous enzyme, with potential applications in biotechnology.

Combined Theoretical and Experimental Study to Unravel the Differences in Promiscuous Amidase Activity of Two Nonhomologous Enzymes

Convergent evolution has resulted in nonhomologous enzymes that contain similar active sites that catalyze the same primary and secondary reactions. Comparing how these enzymes achieve their reaction promiscuity can yield valuable insights to develop functions from the optimization of latent activities. In this work, we have focused on the promiscuous amidase activity in the esterase from Bacillus subtilis (Bs2) and compared with the same activity in the promiscuous lipase B from Candida antarctica (CALB). The study, combining multiscale quantum mechanics/molecular mechanics (QM/MM) simulations, deep machine learning approaches, and experimental characterization of Bs2 kinetics, confirms the amidase activity of Bs2 and CALB. The computational results indicate that both enzymes offer a slightly different reaction environment reflected by electrostatic effects within the active site, thus resulting in a different reaction mechanism during the acylation step. A convolutional neural network (CNN) has been used to understand the conserved amino acids among the evolved protein family and suggest that Bs2 provides a more robust protein scaffold to perform future mutagenesis studies. Results derived from this work will help reveal the origin of enzyme promiscuity, which will find applications in enzyme (re)design, particularly in creating a highly active amidase.

Stability comparison of four lipases and catalytic mechanism during the synthesis of 1,3-di-oleic-2-medium chain triacylglycerols in a trace water-in-oil system: Experimental analyses and computational simulations

In the present study, a kind of structured lipids, namely 1,3-di-oleic-2-medium chain (OMO) triacylglycerols, were synthesized through lipase-catalyzed reactions using coconut oil and rapeseed acid as materials in a trace water-in-oil system. Experimental analysis and computational simulations were undertaken to compare the stability of four lipases including Lipozyme RMIM, Lipozyme TLIM, Novozym 435, and Aspergillus oryzae immobilized lipase (AOIM), and illustrate catalytic mechanism of Novozym 435 during the synthesis of OMO. Fourier transform infrared and molecular dynamics simulation results demonstrated that a decrease in ordered structure (α-helix and β-sheet) led to a reduction in enzyme activity. Compared with Lipozyme RMIM and Novozym 435, Lipozyme TLIM and AOIM exhibited better stability due to a short-chain lid in TLIM, which covers activity sites, and hydrogen bonds formed between activity center of AOIM and water. Among four lipases, AOIM exhibited best catalytic performance: a OMO yield of 30.7% at 3 hr and a good stability of long term (48 hr). Density functional theory results demonstrated that specifically, during the synthesis of OMO triacylglycerol, the addition of Novozym 435 (derived from Candida antarctica lipase B, CALB) substantially lowered reaction barriers (64.4 KJ/mol with CALB vs. 332.7 KJ/mol with no lipase), aiding in the generation of OMO because of the formations of transitional tetrahedral intermediates. A trace water-in-oil system was a green and efficient alternative for lipase-catalyzed production of OMO, and this study provided crucial insights into the stability/instability and catalytic mechanisms of lipase in the synthesis of structured lipids. PRACTICAL APPLICATIONS: We compared the stability of Lipozyme RMIM, Lipozyme 435, Lipozyme TLIM, and AOIM during the synthesis of the OMO triacylglycerols in a trace water-in-oil system, where exhibited a high catalytic activity of lipase in water-oil interface. AOIM had the highest stability and showed the best catalytic performance due to the formation of hydrogen bonds. Besides, for the first time, the transition tetrahedral structure was revealed in the enzymatic synthesis of medium- and long-chain triacylglycerols. This study provides a rational approach to compare lipase stability and a possible hint to choose appropriate enzyme in a specific catalytic condition.

Computational Studies of Candida Antarctica Lipase B to Test Its Capability as a Starting Point To Redesign New Diels-Alderases

The design of new biocatalysts is a target that is receiving increasing attention. One of the most popular reactions in this regard is the Diels-Alder cycloaddition because of its applications in organic synthesis and the absence of efficient natural enzymes that catalyze it. In this paper, the possibilities of using the highly promiscuous Candida Antarctica lipase B as a protein scaffold to redesign a Diels-Alderase has been explored by means of theoretical quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. Free energy surfaces have been computed for two reactions in the wild-type and in several mutants with hybrid AM1/MM potentials with corrections at M06-2X/MM level. The study of the counterpart reactions in solution has allowed performing comparative analysis that render interesting conclusions. Since the dienophile anchors very well in the oxyanion hole of all tested protein variants, the slight electronic changes from reactant complex to the transition state suggest that mutations should be focused in favoring the formation of reactive conformations of a reactant complex that, in turn, would reduce the energy barrier.