Biomolecular Chemical Simulations on Enantioselectivity and Reactivity of Lipase Enzymes to Azulene Derivatives

Abstract

Biomolecular chemical simulations have recently become a useful research method in the fields of organic chemistry and bioscience. In the last few years, we have been focusing on the biomolecular computational simulation on lipase enzyme and ligand complexes to predict the enantioselectivity and reactivity of lipases toward non-natural organic compounds. In this paper, we describe the molecular simulations including molecular dynamics (MD) and fragment molecular orbital (FMO) calculations for the complexes of Candida antarctica lipase type A (CALA) and trifluoromethylazulene alcohol derivatives. From the MD calculations, we found that the fast-reacting enantiomer of esters with high enantioselectivity stays in the vicinity of the active site of CALA, while the slow-reacting enantiomer leaves the active site of CALA. On the other hand, both (R)- and (S)-enantiomers of ester with low ensntioselectivity were found to keep near to near the active site of CALA. Further, for the esters that do not react with lipase enzyme, we found that both (R)- and (S)-enantiomers move away from the active site of lipase enzyme. From the FMO calculations, we found that each fast-reacting enantiomer of esters with high enantioselectivity strongly interacts with certain particular amino acid residues in CALA containing Asp95, while both (R)- and (S)-enantiomers of ester with low enantioselectivity interact with same amino acid residues in CALA including Asp95. These results suggest that it is possible to predict not only the enantioselectivity but also the reactivity of CALA and to identify the amino acid residues important to the enzymatic reaction. Therefore, we consider that our computational simulations would be a useful method for predicting and understanding the reactivity and the enantioselectivity of lipase-catalyzed biotransformations.

Keywords

bioactivity azulene derivatives lipase enantioselectivity biomolecular chemical simulation organic synthesis

Introduction

Today, the dramatic developments in computer technology and simulation software have made it possible to perform biomolecular chemical calculations not only for small organic compounds but also for biopolymers such as proteins. In Particular, simulation methods such as molecular dynamics (MD) calculation and fragment molecular orbital (FMO) calculation^1–4 have been utilized as a research method in the fields of organic chemistry and bioscience. FMO calculation is a quantum chemical calculation method in which biopolymers such as proteins are divided into fragments (amino acid residues), and enable a detailed evaluation of the interaction energy between compounds and each amino acid residue in protein. From this fact, FMO calculation is expected to be an effective method for molecular design of drug discovery and functional new compounds and for modification of enzyme function. However, due to problems with the accuracy of simulations and the analytical environment, it seems that the manufacturing which connected computational simulations and experiments well is still insufficient. In other words, the methodology, such as the analytical validity and the creativity, to make good use of the computational simulations has not been established. Therefore, the computational simulation is not necessarily widespread enough in the industry for the actual manufacturing. As a first step toward solving these problems, we have aimed to reproduce and even to predict the experimental results by the computational simulations. In recent years, we are focusing on biomolecular chemical simulations on the complexes of lipases and organic compounds to predict the enantioselectivity and reactivity of lipase in organic synthesis, and have selected lipases such as Burkholderia cepacia lipase (BCL) and Candida antarctica lipase typeB (CALB) as target proteins.^5,6

Lipase enzymes are widely used as environmentally friendly and efficient biocatalysts for chemoselective, regioselective, and enantioselective reactions under extremely mild conditions due to their ease of use. Lipases are very special enzymes having a peculiar mechanism of reaction.^7–11 Lipases are capable of recognizing many different substrates, and can catalyze hydrolysis, acidolysis, esterification, transesterification, and amination reactions.^8,12–22 Enzymes with excellent catalytic properties have been well studied as catalysts because of their unique physicochemical behavior and have been recently used as a potential biocatalyst^15,23–25 in a large number of biotechnological sciences, such as dairy products, detergents, pharmaceuticals, chemicals,²⁶ agriculture products, and oil chemistry.²⁷ From these facts, the production and utilization of enzymes may enable enzymes to be a better alternative of chemical catalysts.^8,9,26

Until a dozen years ago, our research group has investigated the enantioselectivity of such enzyme lipases as BCL, CALB, and Candida antarctica lipase typeA (CALA) through their biocatalytic synthetic reactions with lipases.^28–31 The degree of enantioselectivity is given by an experimentally obtained E value that is an enantiomeric ratio for each lipase-catalyzed reaction.³² As the E value increases, the enantioselectivity becomes higher. Although a large amount of experimental data is available for BCL, CALB, and CALA lipases, including our results, the mechanism responsible for their enantioselectivity is still not completely elucidated.

In the last few years, we have performed biomolecular chemical simulations, including MD calculation and FMO calculation, on the complexes of BCL and CALB with various primary and secondary alcohol esters.^5,6 We have shown that our simulation is a useful tool for predicting the enantioselectivity of BCL and CALB for various alcohol esters.^5,6

In this paper, we performed the biomolecular chemical simulations toward CALA complexes with the 7 different trifluoromethylazulene secondary alcohol esters (Figure 1). CALA is made from a single polypeptide chain bearing 431 amino acid residues. The active site amino acid residue of CALA is Ser184 and the catalytic triad comprises Ser184, Asp334, and His366.³³ CALA has 2 disulfide bonds, Cys101-Cys273 and Cys350-Cys394. The substrate esters, azulene and its derivatives, comprise “a characteristic fused 5-7 bicyclic aromatic ring system” and are typical “nonbenzenoid aromatic hydrocarbons.”³¹ Azulenes are interesting substances due to their expected reaction specificity and physicochemical and pharmacological properties.³¹ Especially, it is well known that the azulene derivative, sodium guaiazulene sulfonate, has been used as a blue dye in gargles. Over several decades, our research group Konan Chemical Industry Co., Ltd has been developing and producing the water-soluble azulene, sodium guaiazulene sulfonate. They have a lot of experimental data on the azulene derivatives and have reported a part of them in 2004.³¹ In this work, we used their previous experimental data³¹ for comparison with the present computational results. In Figure 1, we show the 7 different trifluoromethylazulene secondary alcohol esters and experimentally obtained E values of CALA toward each azulene alcohol ester.³¹

Figure 1.

Seven different trifluoromethylazulene secondary alcohol esters and experimentally obtained E values of CALA toward each azulene alcohol ester.³¹ CALA shows high enantioselectivity toward esters 1-4 with a large E value, whereas for esters 5 with a very small E value, CALA gives low enantioselectivity. Esters 6 and 7 do not react with CALA.

As shown in Figures 1, CALA shows high enantioselectivity toward esters 1-4 with a large E value, whereas for esters 5 with a very small E value, CALA gives low enantioselectivity. It is observed that for the esters 1-4 with high enantioselectivity, (S)-enantiomer react faster than (R)-enantiomer.³¹ Thus, for esters 1-4, (S)-enantiomer is “the fast-reacting enantiomer,” and (R)-enantiomer is “the slow-reacting enantiomer.” For ester 5, due to its low enantioselectivity, it is difficult to determine whether (R)- or (S)-enantiomer is “the fast-reacting enantiomer.” Here, the azulene alcohol esters 6 and 7 are not converted by CALA. In addition, azulene alcohol ester 1 has been found to react with CALA, but not with BCL and CALB.

Through a series of our work combining of experimental chemistry and computational simulations, we aim to obtain knowledge on predicting lipase enantioselectivity and the extent of the enantioselectivity, and on elucidating the mechanism for enantiomeric recognition and reactivity of enzyme lipases.

Materials and Methods

Construction of CALA Complexes with Azulene Derivatives and Structural Optimization

The structures of (R)- and (S)-enantiomers of 7 different azulene alcohol esters were modeled on a computer using GaussView. GaussView is a graphical user interface used with Gaussian, a quantum chemistry calculation program. Then, the structures of esters were optimized by 2 steps quantum chemical calculations at HF/6-31G* level and at B3LYP/6-31G* level using Gaussian03.³⁴ Next, we downloaded the x-ray crystallographic structure of CALA from Protein Data Bank³⁵ (PDB code: 3GUU). This pdb structure was obtained by x-ray structural analysis with resolution of 2.10 Å and does not include the structural information of hydrogen atoms. In this study, we assigned the protonation state of histidine as delta protonated state (HID). And all crystallization water was removed. We used AMBER11 software³⁶ both of construction of the CALA-ester complexes and the structural optimization of the complexes by the energy minimization. Here, N-and C-terminus of CALA were capped by NH3⁺ and COO⁻, respectively.⁴ Each azulene alcohol ester is manually placed near the active site Ser184 of CALA (Figure 2(a)). We assigned the charge of esters by AM1-BCC of antechamber³⁷ in AMBER11. As a first step of our study, we focused on the enzymatic reactions in water. We solvated the structures of the CALA-ester complexes with TIP3P model water molecules within 8.0 Å of each complex. The number of TIP3P model water molecules placed around each complex is over 10,000. Counterions were placed to neutralize the systems. Subsequently, we carried out the energy minimizations for all atoms in the presence of TIP3P water molecules by a molecular mechanics (MM) calculation using AMBER11 with the ff99 and general AMBER force field for CALA and azulene alcohol esters, respectively.^38,39 The obtained structure takes an atomic configuration that minimizes the energy of the system. Here, we took the disulfide bonds in CALA into account in our all calculations: MM, MD, and FMO calculations.

Figure 2.

(a) Complex of CALA and (S)-enantiomer of ester 1 [(S)-1]. (b) Structure near the active site of CALA and (S)-1 complex. C-O interatomic distance, R_C-O, between the carbonyl carbon of (S)-1 and the oxygen of the active site Ser184 side chain OH in CALA-ester complex is shown with a black double-headed arrow. Color code: green C, red O, blue N, white H.

Molecular Dynamics Calculations

After energy minimization, CALA-ester complexes with TIP3P model water molecules were subsequently subjected to a MD calculation over a period of 2000 ps at 300 K using AMBER11. We employed periodic boundary condition by controlling the pressure. The temperature was kept constant using the weak-coupling algorithm⁴⁰ with a coupling time of 1.0 ps. We constrained only bond lengths involving hydrogen atoms by the SHAKE method.⁴¹ We calculated non-bonding interactions using the Particle Mesh Ewald method with 10.0 Å cutoff distance. The time step for the MD calculation was 1 fs. The whole system temperature was gradually increased by heating to 300 K for 50 ps and kept at 300 K for the next 2000 ps. Then, we estimated the C-O interatomic distance, R_C-O, between the carbonyl carbon of azulene alcohol esters and the oxygen of the active site Ser184 side chain OH in each CALA-ester complex (Figure 2(b)) using VMD (Visual Molecular Dynamics) Version 1.8.6⁴² and examined the time dependence of R_C-O during the MD trajectory. Here, the C-O interatomic distance, so called R_C-O, is also defined in Refs^5,6. Moreover, to examine the reactivity of lipase enzymes in detail, we carried out the similar MD simulations on the complexes of BCL with azulene alcohol ester 1 and of CALB and azulene alcohol ester 1 which are prepared by above method. These complexes are unreactive ones.

Fragment Molecular Orbital Calculations

We performed the FMO calculations for the complexes in which the ligand remains in the active site. Since the depth of the active pocket of CALA, the distance from the entrance of the active pocket to Ser184 is about 9.0 Å, we concluded that the ligand is separated from the active site for complexes with R_CO > 9.0 Å at 2000 ps of MD calculations. Therefore, we selected the complexes of CALA and (S)-enantiomer of esters 1-4 with high enantioselectivity, and the complexes of CALA and both (R)- and (S)-enantiomers of ester 5 with low enantioselectivity, for the FMO calculation. After MD calculations, for the structures at 2000 ps, we carried out the structural optimizations by MM calculations using AMBER11. Then, we removed the surrounding TIP3P model water molecules and counterions from the system. Subsequently, we performed an FMO calculation on the obtained complexes using ABINIT-MP/BioStation program^43–45 at MP2/6-31G level. The Cholesky decomposition with adaptive metric method⁴⁶ was applied to the FMO calculations. In our FMO calculations, each amino acid residue except cysteines with disulfide bond of CALA and the enantiomer of azulene alcohol esters were treated as a single fragment, respectively. We treated 2 cysteines bonding with disulfide bond as one fragment. The computations were performed under vacuum conditions. We obtained the interaction energy (inter-fragment interaction energy: IFIE) between the azulene alcohol esters and the amino acid residues in CALA.

Results and Discussion

Molecular Dynamics Calculations

Figures 3 and 4 show the time dependence of R_C-O in the complexes of CALA with the (R)- and (S)-enantiomer of azulene alcohol esters 1-4 and azulene alcohol ester 5, respectively. As shown in Figure 3, for the esters 1-4 with a large E value (high enantioselectivity), we found that R_C-O of the fast-reacting (S)-enantiomer did not change almost around 5 Å, but R_C-O of the slow-reacting (R)-enantiomer increased with time. This result indicates that the fast-reacting enantiomer of esters with high enantioselectivity stays in the vicinity of the active site of CALA, while the slow-reacting enantiomer leaves the active site of CALA. On the other hand, as shown in Figure 4, for ester 5 with a very small E value (low enantioselectivity), both (R)- and (S)-enantiomers were found to keep near to the active site of CALA. Similar MD calculations were performed for some initial structures, and we obtained the same results mentioned above.

Figure 3.

Time dependence of the C–O interatomic distance, R_C-O, obtained by MD calculations on the complexes of CALA with (R)- and (S)-enantiomers of esters 1-4 with a large E value (high enantioselectivity). For all complexes, R_C-O of (R)-enantiomer is larger than that of (S)-enantiomer.

Figure 4.

Time dependence of the C-O interatomic distance, R_C-O, obtained by MD calculations on the complexes of CALA with (R)- and (S)-enantiomers of ester 5 with a very small E value (low enantioselectivity). R_C-O of both (R)- and (S)-enantiomers are about 5 Å.

Then, we estimated the difference in R_C-O between the (R)- and (S)-enantiomer complexes, ΔR_C-O, at 2000 ps of the MD trajectory. Namely, ΔR_C-O is defined by the following equation:

Δ R_{C - O} = R_{C - O, R} – R_{C - O, S}

Here, R_{C-O, R} is R_C-O for the (R)-enantiomer complex, and R_{C-O, S} is R_C-O for the (S)-enantiomer complex. Figure 5 shows the relationship between the experimentally obtained E value and ΔR_C-O for alcohol esters 1-5. For the esters 1-4 with a large E value, ΔR_C-O was more than 7.0 Å, while that for the ester 5 with a very small E value, ΔR_C-O was 0.2 Å. The latter indicates that (R)-enantiomer of ester 5 is slightly distant from the active site of CALA, but both (R)- and (S)-enantiomers are located at similar position in the active site of CALA. This suggests that both (R)- and (S)-enantiomers of ester 5 react with CALA. It is seemed that the values of ΔR_C-O obtained by MD calculations are related to the experimentally obtained E values, as shown in previous works.^5,6 In the present simulation, we are not yet able to calculate the E values directly, but we consider that we may predict E values by examining in detail the correlations between the calculated values, such as R_C-O and ΔR_C-O, and many accurate experimental E values.

Figure 5.

Relationship between the difference in the C-O interatomic distance of (R)- and (S)-enantiomer complexes, ΔR_C-O, obtained by MD calculations and experimentally obtained E values for esters 1-5. For the esters 1-4 with a large E value, ΔR_C-O is more than 7.0 Å, while that for the ester 5 with a very small E value, ΔR_C-O is 0.2 Å.

In addition, we show the MD computations on the unreactive complexes such as CALA with ester 6, CALA with ester 7, BCL with ester 1, and CALB with ester 1 in Figure 6. For the esters that do not react with lipase enzyme, we found that both (R)- and (S)-enantiomers move away from the active site of lipase enzyme.

Figure 6.

Time dependence of the C-O interatomic distance, R_C-O, obtained by MD calculations on the unreactive complexes such as CALA with ester 6, CALA with ester 7, BCL with ester 1, and CALB with ester 1. For both (R)- and (S)-enantiomers, R_C-O increases with time.

As a supplementary addition, we show the time dependence of root mean square deviation (RMSD) during the MD simulations on the complexes of lipases and trifluoromethylazulene alcohol esters in Supplemental Figures S1, S2, and S3.

Fragment Molecular Orbital Calculations

As shown in our previous works,^5,6 the IFIE analysis is a useful method for gaining detailed information on protein-ligand binding. Figure 7 shows the IFIE values between the fast-reacting enantiomer of esters 1-4 with high enantioselectivity and the selected amino acid residues in CALA. We found that each fast-reacting enantiomer of esters with high enantioselectivity showed a large interaction energy with Asp95 in CALA. In addition to the interaction with Asp95, it was also found that each fast-reacting (S)-enantiomer interacted with certain particular amino acid residues in CALA such as Leu225, Phe230, Phe233, Glu336, and Ile336. These amino acid residues are located near the active site of CALA. In contrast, the interactions of each slow-reacting (R)-enantiomer with these amino acid residues were very weak. The reason for this result is considered to be that the slow-reacting enantiomer departs from the active site of CALA. In this paper, we omitted showing the IFIEs for the slow-reacting enantiomer. As an example, the structures of the fast-reacting enantiomers of esters 1, 2 and the interacting amino acid residues are illustrated in Figure 8. The IFIE between Asp95 and (S)-2 is larger than that between Asp95 and (S)-1. The distance between the oxygen atom in side chain of Asp95 and (S)-2 is about 2.4 Å, whereas that between the oxygen atom in side chain of Asp95 and (S)-1 is about 3.4 Å. The IFIE between Glu335 and (S)-1 is larger than that between Glu335 and (S)-2. The side chain of Glu335 orients toward (S)-1, but not toward (S)-2.

Figure 7.

Results of FMO calculations: IFIE between the fast-reacting (S)-enantiomer of esters 1-4 with high enantioselectivity and the selected amino acid residues in CALA. (S)-enantiomer of esters show a large interaction energy with Asp95 in CALA. In addition to Asp95, each (S)-enantiomer interacts with certain particular amino acid residues in CALA such as Leu225, Phe230, Phe233, Glu335, and Ile336. These amino acid residues are located near the active site of CALA.

Figure 8.

Structures of (S)-enantiomers of esters 1(a), 2(b), and the interacting amino acid residues. The distance between the oxygen atom in side chain of Asp95 and (S)-1 is about 3.4 Å (a), whereas that between the oxygen atom in side chain of Asp95 and (S)-2 is about 2.4 Å (b). The side chain of Glu335 orients toward (S)-1, but not toward (S)-2.

On the other hand, as shown in Figure 9, for ester 5 with low enantioselectivity, both (R)- and (S)-enantiomers were found to interact with the same amino acid residues containing Asp95. We also found that the interactions of both enantiomers with amino acid residues in CALA were similar to each other in their intensities and patterns. Figure 9 shows the IFIE values between (R)- and (S)-enantiomers of esters 5 and the selected amino acid residues in CALA.

Figure 9.

Results of FMO calculations: IFIE between the ester 5 with low enantioselectivity and the selected amino acid residues in CALA. Both (R)- and (S)-enantiomers interact with the same amino acid residues containing Asp95.

Interestingly, regardless of the extent of enantioselectivity, the interaction between each enantiomer and the active site amino acid residue Ser184 was very weak or hardly observed. These FMO computational results suggest that the IFIEs revealed in the present calculations can also be related to the enantioselectivity of CALA toward the azulene alcohol esters observed by synthetic investigations.

Conclusions

To predict the lipase enantioselectivity and the extent of enantioselectivity toward non-natural organic compounds, for the complexes of CALA and 7 different trifluoromethylazulene secondary alcohol esters, we performed the MD calculation over a period of 2000 ps and estimated the C-O interatomic distance, R_C-O, between the carbonyl carbon of azulene alcohol esters and the oxygen of the active site Ser184 side chain OH for each CALA-ester complex. The FMO simulations on the CALA-ester complex structures obtained by MD calculations were subsequently carried out and the IFIE between esters and amino acid residues in CALA were calculated.

From the MD calculations, we found that the fast-reacting (S)-enantiomer of esters 1-4 with high enantioselectivity stays in the vicinity of the active site of CALA, while the slow-reacting (R)-enantiomer leaves the active site of CALA. On the other hand, both (R)- and (S)-enantiomers of ester 5 with low enantioselectivity were found to keep near to the active site of CALA. In addition, for the esters that do not react with lipase enzyme not only CALA but also BCL and CALB, both (R)- and (S)-enantiomers were found to leave the active site of lipase enzymes. The results strongly suggest that not only the enantioselectivity but also the reactivity of lipases such as CALA, BCL, and CALB can be predicted by the MD simulations.

Another aim of our study is to make clear the mechanism of enantiomeric recognition of enzyme lipases. For the esters showing high enantioselectivity, it is also found that each fast-reacting enantiomer shows strong interaction with the particular amino acid residues including Asp95 in CALA. On the contrary, each slow-reacting enantiomer hardly interacted with these amino acid residues. The reason for this result is considered to be that the slow-reacting enantiomer of esters with high enantioselectivity leaves the active site of CALA. In contrast, for the esters bearing low enantioselectivity, we observed that both (R)- and (S)-enantiomers interact with the same amino acid residues containing Asp95. The IFIE for the amino acid residue Asp95 is particularly noteworthy. It is predictable that Asp95 may play an important role in the chiral recognition of enantiomers by CALA. Alternatively, Asp95 might play the role of anchoring the esters to the active site of CALA during the reaction. This result is similar to that of His286 in BCL and Thr40 in CALB.^5,6

The very weak interaction between esters and the active site amino acid residue Ser184 in CALA may be due to the fact that Ser184 is not yet activated. The additional more accurate FMO calculations that take into account the charge state of Ser184 might give more detailed information on the interactions between esters and amino acid residues in CALA. Otherwise, we may need to examine the structure changes of complexes and the time change of IFIE between esters and Ser184. Moreover, we made CALA-ester complexes manually using AMBER in the present work, but we are planning to carry out a docking simulation for constructing CALA-ester complexes. Moreover, in order to make our computational method more accurate and realistic, we have some issues to be solved in the future: (1) protonation states or charge states of the particular amino acid residues, such as histidine and the active site amino acids including Ser184, (2) identification of the amino acid residues involved in the molecular recognition mechanism or reaction mechanism/pathway of CALA and esters, (3) simulations under various situations, such as in organic solvents and exposure to the medium, (4) roles of surrounding water molecules or organic solvents (solvent effects) in CALA-ester binding, (5) comparison with new experiments with higher accuracy. Even though there are various issues at present, we believe that our computational approach including the MD calculation and FMO calculations may enable us to elucidate the mechanism responsible for enantiomeric recognition of enzyme lipases. Previously, we have performed similar computational simulations on the complexes of HIV-1 protease and its inhibitors, and found the following results: (1) a strong correlation between the experimentally obtained activity values and the calculated values such as binding energy, and (2) identification of the amino acid residues that strongly interact with the inhibitor.⁴⁷ Further, we are planning the computational simulations toward the complexes of a human protease and its inhibitor. We consider that our computational method is sufficiently applicable not only to the lipase-ester complexes in this work but also to other systems, such as the protease-inhibitor complexes and the enzyme-surrogate complexes. Through our continuous works, we have reached a conclusion that biomolecular chemical simulations are useful tools for predicting and understanding the reactivity and the enantioselectivity of lipase-catalyzed biotransformations.

Supplemental Material

sj-docx-1-npx-10.1177_1934578X221108572 - Supplemental material for Biomolecular Chemical Simulations on Enantioselectivity and Reactivity of Lipase Enzymes to Azulene Derivatives

Supplemental material, sj-docx-1-npx-10.1177_1934578X221108572 for Biomolecular Chemical Simulations on Enantioselectivity and Reactivity of Lipase Enzymes to Azulene Derivatives by Yoichiro Yagi, Takatomo Kimura and Makoto Kamezawa in Natural Product Communications

Footnotes

Acknowledgments

We would like to express our gratitude to Konan Chemical Industry Co., Ltd for their great cooperation from the early stages of this research. This research was performed as activities of the FMO drug design consortium (FMODD).

Declaration of Conflicting Interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The authors received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Yoichiro Yagi

Supplemental material

Supplemental material for this article is available online.

References

Kitaura

Ikeo

Asada

Nakano

Uebayasi

. Fragment molecular orbital method: an approximate computational method for large molecules. Chem Phys Lett. 1999;313(3-4):701-706. doi: 10.1016/S0009-2614(99)00874-X

Fedorov

Kitaura

. The Fragment Molecular Orbital Method: PRACTICAL APPLICATION TO LARGE MOLECULAR SYSTEM. CRC Press; 2009.

Tanaka

Mochizuki

Komeiji

Okiyama

Fukuzawa

. Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems. Phys Chem Chem Phys. 2014;16(22):10310-10344. doi:10.1039/C4CP00316K

Mochizuki

Tanaka

Fukuzawa

. Recent Advances of the Fragment Molecular Orbital Method Enhanced Performance and Applicability. Springer; 2021.

Yagi

Tanaka

Imagawa

, et al. Large-Scale biomolecular chemical computations toward the prediction of Burkholderia cepacia lipase enantioselectivity. J Adv Simul Sci Eng. 2014;1(1):141-160. doi:10.15748/jasse.1.141

Yagi

Kimura

Kamezawa

Naoshima

. Biomolecular chemical simulations toward elucidation of the enantioselectivity and reactivity of lipases in organic synthesis. CBI J. 2018;18:21-31. doi:10.1273/cbij.18.21

dos Santos

JCS

Garcia-Galan

Rodrigues

de Sant’Ana

Gonçalves

LRB

Fernandez-Lafuente

. Stabilizing hyperactivated lecitase structures through physical treatment with ionic polymers. Process Biochem. 2014;49(9):1511-1515. doi:10.1016/j.procbio.2014.05.009

Rios

Neto

DMA

dos Santos

JCS

Fechine

PBA

Fernández-Lafuente

Gonçalves

LRB

. Comparison of the immobilization of lipase from Pseudomonas fluorescens on divinylsulfone or p-benzoquinone activated support. Int J Biol Macromol. 2019;134:936-945. doi:10.1016/j.ijbiomac.2019.05.106

Fernandez-Lopez

Bartolome-Cabrero

Rodriguez

dos Santos

Rueda

Fernandez-Lafuente

. Stabilizing effects of cations on lipases depend on the immobilization protocol. RSC Adv 2015;5(102):83868–83875. doi:10.1039/c5ra18344h

10.

da S. Souza

de Oliveira

Alexandre

JYNH

, et al. A comprehensive review on the use of metal–organic frameworks (MOFs) coupled with enzymes as biosensors. Electrochem. 2022;3(1):89-113. doi:10.3390/electrochem3010006

11.

Cavalcante

FTT

Cavalcante

ALG

de Sousa

Neto

dos Santos

JCS

. Current Status and future perspectives of supports and protocols for enzyme immobilization. Catalysts. 2021;11(10):1222. doi:10.3390/catal11101222

12.

Bonazza

Manzo

dos Santos

JCS

Mammarella

. Operational and thermal stability analysis of Thermomyces lanuginosus lipase covalently immobilized onto modified chitosan supports. Appl Biochem Biotechnol. 2018;184(1):182-196. doi:10.1007/s12010-017-2546-9

13.

Virgen-Ortíz

dos Santos

JCS

Ortiz

, et al. Lecitase ultra: a phospholipase with great potential in biocatalysis. Mol Catal. 2019;473:110405. doi:10.1016/j.mcat.2019.110405

14.

Bezerra

Monteiro

RRC

Neto

DMA

, et al. A new heterofunctional support for enzyme immobilization: PEI functionalized Fe3O4 MNPs activated with divinyl sulfone. Application in the immobilization of lipase from Thermomyces lanuginosus. Enzyme Microb Technol. 2020;138:109560. doi:10.1016/j.enzmictec.2020.109560

15.

Moreira

de Oliveira

ALB

, et al. Lipase from Rhizomucor miehei immobilized on magnetic nanoparticles: performance in fatty acid ethyl ester (FAEE) optimized production by the Taguchi method. Front Bioeng Biotechnol. 2020;8:693. doi:10.3389/fbioe.2020.00693

16.

Melo

ADQ

Silva

FFM

dos Santos

JCS

Fernández-Lafuente

Lemos

TLG

Filho

FAD

. Synthesis of benzyl acetate catalyzed by lipase immobilized in nontoxic chitosan-polyphosphate beads. Molecules. 2017;22(12):2165. doi:10.3390/molecules22122165

17.

de Oliveira

UMF

de Matos

LJBL

de Souza

MCM

Pinheiro

dos Santos

JCS

Gonçalves

LRB

. Efficient biotechnological synthesis of flavor esters using a low-cost biocatalyst with immobilized Rhizomucor miehei lipase. Mol Biol Rep. 2019;46(1):597-608. doi:10.1007/s11033-018-4514-z

18.

Monteiro

RRC

Lima

PJM

Pinheiro

, et al. Immobilization of lipase A from Candida Antarctica onto chitosan-coated magnetic nanoparticles. Int J Mol Sci. 2019;20(16):4018. doi:10.3390/ijms20164018

19.

Manoel

Pinto

dos Santos

JCS

, et al. Design of a core–shell support to improve lipase features by immobilization. RSC Adv. 2016;6(67):62814–62824. doi:10.1039/C6RA13350A

20.

da Fonseca

de Freitas Ícaro

BÍB

Soares

, et al. Synthesis, biological activity, and in silico study of bioesters derived from bixin by the CALB enzyme. Biointerface Res Appl Chem. 2022;12(5):5901-5917. doi:10.33263/BRIAC125.59015917

21.

Mota

de Sousa

de Oliveira

ALB

, et al. Biodiesel production from microalgae using lipase-based catalysts: current challenges and prospects. Algal Res. 2022;62:102616. doi:10.1016/j.algal.2021.102616

22.

Verdasco-Martín

Villalba

dos Santos

JCS

Tobajas

Fernandez-Lafuente

Oteroa

. Effect of chemical modification of novozym 435 on its performance in the alcoholysis of camelina oil. Biochem Eng J. 2016;111:75-86. doi:10.1016/j.bej.2016.03.004

23.

Lima

da Silva

de Sousa Fonseca

, et al. Chemoenzymatic synthesis of (S)-pindolol using lipases. Appl Catal A-Gen. 2017;546:7-14. doi:10.1016/j.apcata.2017.08.003

24.

Moreira

de Oliveira

ALB

de Moura

Jr. , et al. Taguchi design-assisted co-immobilization of lipase A and B from Candida Antarctica onto chitosan: characterization, kinetic resolution application, and docking studies. Chem Eng Res Des. 2022;177:223-244. doi:10.1016/j.cherd.2021.10.033

25.

Valério

RBR

Cavalcante

ALG

Mota

, et al. Understanding the biocatalytic potential of lipase from Rhizopus chinensis. Biointerface Res Appl Chem. 2022;12(3):4230-4260. doi:10.33263/BRIAC123.42304260

26.

Lima

PJM

da Silva

Neto

CACG

, et al. An overview on the conversion of glycerol to value-added industrial products via chemical and biochemical routes. Appl Biochem Biotechnol. 2021. doi:10.1002/bab.2098. Early View. (Online Version of Record before inclusion in an issue).

27.

Cavalcante

FTT

Neto

de Aguiar Falcão

, et al. Opportunities for improving biodiesel production via lipase catalysis. Fuel. 2021;288:119577. doi:10.1016/j.fuel.2020.119577

28.

Naoshima

Kamezawa

Tachibana

, et al. Enzymatic preparation of enantiomerically pure alkan-2- and -3-ols by lipase-catalysed hydrolysis with Pseudomonas cepacia in the presence of organic media. J Chem Soc, Perkin Trans 1. 1993;(5):557-561. doi:10.1039/P19930000557

29.

Ohtani

Nakatsukasa

Kamezawa

Tachibana

Naoshima

. Enantioselectivity of Candida Antarctica lipase for some synthetic substrates including aliphatic secondary alcohols. J Mol Catal B: Enzym. 1998;4(1-2):53-60. doi:10.1016/S1381-1177(97)00021-0

30.

Naoshima

Kamezawa

Kimura

, et al. Enantioselectivity and reactivity of Candida and Pseudomonas sp. Lipases for synthetic substrates: biocatalytic synthesis of chiral building blocks with alkanol skeletons. Recent Res Devel Org Bioorg Chem. 2001;4:1-16.

31.

Naoshima

Kimura

Mori

, et al. The first synthesis of chiral azulene alcohols possessing a trifluoromethyl group by lipase-mediated biotransformation. Recent Res Devel Org Bioorg Chem. 2004;6:1-9.

32.

Chen

C-S

Fujimoto

Girdaukas

Sih

. Quantitative analyses of biochemical kinetic resolutions of enantiomers. J Am Chem Soc. 1982;104(25):7294-7299. doi:10.1021/ja00389a064

33.

Ericsson

Kasrayan

Johansson

, et al. x-ray structure of Candida Antarctica lipase A shows a novel lid structure and a likely mode of interfacial activation. J Mol Biol. 2008;376(1):109-119. doi:10.1016/j.jmb.2007.10.079

34.

Frisch

Trucks

Schlegel

, et al. Gaussian03, Revision C.02. Gaussian, Inc.; 2004. https://gaussian.com/g03citation/

35.

wwPDB consortium. Protein data bank: the single global archive for 3D macromolecular structure data. Nucleic Acids Res. 2019;47(D1):D520–D528. doi:10.1093/nar/gky949

36.

Case

Darden

Cheatham

et al. AMBER11. University of California; 2010.

37.

Jakalian

Jack

Bayly

. Fast, efficient generation of high-quality atomic charges. AM1–BCC model: II. Parameterization and validation. J Comput Chem 2002;23(16):1623-1641. doi:10.1002/jcc.10128

38.

Wang

Cieplak

Kollman

. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J Comput Chem. 2000;21(12):1049-1074. doi:10.1002/1096-987X(200009)21:12<1049::AID-JCC3>3.0.CO;2-F

39.

Wang

Wolf

Caldwell

Kollman

Case

. Development and testing of a general amber force field. J Comput Chem. 2004;25(9):1157-1174. doi:10.1002/jcc.20035

40.

Berendsen

HJC

Postma

JPM

van Gunsteren

DiNola

Haak

. Molecular dynamics with coupling to external bath. J Chem Phys. 1984;81(8):3684-3690. doi:10.1063/1.448118

41.

Ryckaert

J-P

Ciccotti

Berendsen

HJC

. Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977;23(3):327-341. doi:10.1016/0021-9991(77)90098-5

42.

Humphrey

Dalke

Schulten

. VMD: visual molecular dynamics. J Mol Graph. 1996;14(1):33-38. doi:10.1016/0263-7855(96)00018-5

43.

Nakano

Kaminuma

Sato

, et al. Fragment molecular orbital method: use of approximate electrostatic potential. Chem Phys Lett. 2002;351(5–6):475-480. doi:10.1016/S0009-2614(01)01416-6

44.

Nakano

Mochizuki

Fukuzawa

Amari

Tanaka

. Developments and applications of ABINIT-MP software based on the fragment molecular orbital method. In: Starikov

Lewis

Tanaka

, eds. Modern Methods for Theoretical Physical Chemistry of Biopolymers. ELSEVIER; 2006:39-52.

45.

Nakano

Mochizuki

Kato

, et al. Developments of FMO methodology and graphical user interface in ABINIT-MP. In: Fedorov

Kitaura

, eds. the Fragment Molecular Orbital Method: Practical Application to Large Molecular Systems. CRC Press; 2009:37-62.

46.

Okiyama

Nakano

Yamashita

Mochizuki

Taguchi

Tanaka

. Acceleration of fragment molecular orbital calculations with Cholesky decomposition approach. Chem Phys Lett. 2010;490(1-3):84-89. doi:10.1016/j.cplett.2010.03.001

47.

Yagi

Hattori

Naoshima

. Large-Scale quantum chemical calculation on the complexes of HIV-1 protease and inhibitors: a relationship of calculation energies with drug effects of inhibitors. Trans Jpn Soc Simul Tech. 2012;4(2):41-50. doi:10.11308/tjsst.4.41

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.03 MB