In Silico Evaluation of Enzymatic Tunnels in the Biotransformation of α-Tocopherol Esters

Motivation: α-Tocopherol is a molecule obtained primarily from plant sources that are important for the pharmaceutical and cosmetics industry. However, this component has some limitations such as sensitivity to oxygen, presence of light, and high temperatures. For this molecule to become more widely used, it is important to carry out a structural modification so that there is better stability and thus it can carry out its activities. To carry out this structural modification, some modifications are carried out, including the application of biotransformation using enzymes as biocatalysts. Thus, the application of a computational tool that helps in understanding the transport mechanisms of molecules in the tunnels present in the enzymatic structures is of fundamental importance because it promotes a computational screening facilitating bench applications. Objective: The aim of this work was to perform a computational analysis of the biotransformation of α-tocopherol into tocopherol esters, observing the tunnels present in the enzymatic structures as well as the energies which correspond to the transport of molecules. Method: To carry out this work, 9 lipases from different organisms were selected; their structures were analyzed by identifying the tunnels (quantity, conformation, and possibility of transport) and later the calculations of substrate transport for the biotransformation reaction in the identified tunnels were carried out. Additionally, the transport of the product obtained in the reaction through the tunnels was also carried out. Results: In this work, the quantity of existing tunnels in the morphological conformational characteristics in the lipases was verified. Thus, the enzymes with fewer tunnels were RML (3 tunnels), LBC and RNL (4 tunnels), PBLL (5 tunnels), CALB (6 tunnels), HLG (7 tunnels), and LCR and LTL (8 tunnels) and followed by the enzyme LPP with the largest number of tunnels (39 tunnels). However, the enzyme that was most likely to transport substrates in terms of α-tocopherol biotransformation (in relation to the Emax and Ea energies of ligands and products) was CALB, as it obtains conformational and transport characteristics of molecules with a particularity. The most conditions of transport analysis were α-tocopherol tunnel 3 (Emax: −4.6 kcal/mol; Ea: 1.1 kcal/mol), vinyl acetate tunnel 1 (Emax: −2.4 kcal/mol; Ea: 0.1 kcal/mol), and tocopherol acetate tunnel 2 (Emax: −3.7 kcal/mol; Ea: 2 kcal/mol).


INTRODUCTION
Better known as vitamin E, α-tocopherol was discovered in 1922; it is a molecule with a fat-soluble characteristic that can be extracted from different plant sources and has great importance for the body's physiological functions acting on oxidative stress, age, inflammatory diseases, and cardiovascular diseases as well as neurodegenerative diseases (Bartolini et al., 2021). Although there are other forms of tocopherol (β, γ, and δ), which perform biological activities, only α-tocopherol has shown to be essential in the prevention of diseases and can be called a vitamin (Azzi, 2018). The most well-known characteristic of α-tocopherol is its antioxidant activity as it has the function of eliminating free radicals, donating its phenolic hydrogen present in its structure to free oxygen radicals. When applied in the topic form, α-tocopherol realizes a primordial compound function in care with pollution, UV-rays, and antiaging. However, as α-tocopherol is unstable under light, oxygen, and high temperatures, an alternative for inserting this component in formulations in a different way with better properties is to perform the structural modification (Freitas et al., 2008;Ben-Shabat et al., 2013;Abla and Banga, 2014). Further route of α-tocopherol is the structural modification to tocopherol esters since these esters have advantages over the unmodified molecule. Some esters such as acetate, succinate, ferulate, and α-tocopherol nicotinate have better stability than α-tocopherol and are more widely used in different products of the pharmaceutical industry for internal use (nutraceuticals and supplements) and topical products (products of hygiene, cosmetics, and semi-solid formulations) (Yang and McClements, 2013). The topical use, in addition to carrying out the antioxidant action, it slows down the degradation process of collagen existing in the skin, moreover to protect against ultraviolet rays, which suggests a wide use in antiaging and sun protection products (Abla and Banga, 2014). To obtain these esters, the industry uses chemical routes to obtain the final product.
The most commonly used catalysts in reactions are Lewis acid, BrÔnsted acid, or a combination of acids (ZnCl 2 , ZnCl 2 /HCl, BF 3 , AlCl 3 , and FeCl 2 /Fe/HCl). These catalysts have a low cost, but their large-scale production can generate numerous disadvantages such as wasted water, raw material, low recovery rate and catalyst use, and generating more by-products (Shi and Zhang, 2017;Zou et al., 2021).
Thus, the use of enzymatic catalysis has been a useful alternative in the reaction process with conditions of greater specificity and speed. Among the various types of enzymes that can be applied in industrial processes, lipases are biocatalysts that have different applications in industry and are considered to be quite useful for being present in industrial processes such as the food, beverage, pharmaceutical, and cosmetics industries (Sarmah et al., 2018). Therefore, obtaining a product that has a process with a more environmentally favorable route, such as the use of biocatalysts, is of great interest.
The first report of obtaining tocopherol esters by the enzymatic route was obtained by Torres et al. (2008) in which the acylation of the phenolic group of α-tocopherol was carried out through the transesterification reaction using vinyl acetate in 2-methyl-2-butanol to obtain α-tocopherol acetate using different lipases as biocatalysts. Lipase (triacylglycerol acyl hydrolase EC 3.1.1.3) has several advantages. They have the ability to carry out acetylation, transesterification, hydrolysis, and esterification reactions in order to obtain different products for a wide industrial range. They are safe for human consumption and an abundant source of essential amino acids, antioxidants, and minerals and have found application in health foods, supplements, and cosmetics. These enzymes have been extensively studied due to their different biochemical and physiological properties (Kumar et al., 2019;Kublicki et al., 2021). Since, the catalytic residues present in its triad occur in a different order in the protein sequence for lipases (Ser-His-Asp/ Glu). This group of enzymes has a molecular weight ranging between 20 and 75 kDa, optimal pH for activity 4-9, and performing activity in the temperature range between 25 and 70°C (Uppenberg et al., 1994;Rios et al., 2018). The screening was performed with different lipases [lipase B from Candida antarctica (Novozyme 435), lipase from Thermomyces lanuginosus (TL IM), lipase from Rhizomucor miehei (RM IM), lipase from Pseudomonas cepacia (Lipase 1 PS-C), lipase from Alcaligenes sp (Lipase PLG), lipase from Alcaligenes sp (Lipase QLG), cholesterol esterase from Pseudomonas sp., cholesterol esterase from Ophiostoma piceae (immobilized in Dilbeads TA), lipase from porcine pancreas, lipase from Candida rugosa, protease from Bacillus licheniformis, esterases from bovine rumen microflora, and esterase from the Urania hypersaline basin]. In 2011, Chunhua et al. (2011) also synthesized enzymatically, using Candida antarctica lipase B (Novozym 435), to obtain tocopherol succinate, modifying the biocatalyst to improve its catalytic performance using acetic anhydride, propionic anhydride, and succinic anhydride separately, and with the best modifying agent succinic anhydride (1: 5) yielded 94.4%. Jiaojiao et al. (2021) used a nanogel-modified Candida rugosa lipase to obtain vitamin E succinate and obtained a yield of approximately 62% with a reaction time of 15 h. Xin et al. (2011), on the other hand, verified the ability of lipases to catalyze through the transesterification reaction of α-tocopherol and ethyl ferulate to obtain tocopherol ferulate, and it can be observed that only Candida antarctica lipase B and Candida rugosa lipase catalyzed the transesterification reaction. A yield of 25.2% tocopherol ferulate was obtained using a 5:1 molar ratio of α-tocopherol to solvent-free ethyl ferulate in 72 h of reaction. To predict the behavior of substrates in the face of the processes that biocatalysis uses, applying molecular dynamics tools has been a useful strategy.
Computational screening is a well-established technique that is used to observe the structure of components used in the pharmaceutical industry and in the academic environment (Wolf et al., 2010). It is a strategy that assists in the first stages of the drug preparation process to generate a large database of compounds in one given receptor and then selecting the binders obtained from commercial libraries Marques et al., 2020). The tracking and analysis of the behavior of the ligand along the path of biomolecular systems in molecular dynamics represent a strategy that enhances the protein design process, to highlight important regions for the transport of ligands, that is, molecular tunnels, channels, and gates, which establish the association of the ligand with dissociation mechanisms (Surpeta et al., 2020).
Enzyme tunnels and channels are the fundamental transport pathways that allow ligands to move from the outer to the inner medium of the protein. Consequently, there are web servers that can be applied to study the transport of molecules in tunnels such as MolAxis, MOLE, ChExVis, 3V, and CAVER. Thus, the use of the Caver web server (a combination of Caver and CaverDock online) is an innovative tool to investigate the passage of the ligand through the channels present in biomolecules Marques et al., 2020;Singh et al., 2021). Enzymatic tunnels are characterized by having a single opening that allows the transport of substrates, products, and other molecules inside and outside the active site. These enzymatic structures are essential for a good catalytic behavior of the enzyme, as they directly affect its specificity, stability, and substrate activity. The conformation and physicochemical properties of the tunnels also function as a way of protecting the hydrophobic nucleus of the enzymes to restrict the access of solvent molecules and inhibitors .
Recently, identification studies on the catalytic behavior of enzymes as fatty acids developed by our group were essential for understanding the reactions at the molecular level. However, the study of the transport of substances through enzymatic tunnels has not yet been explored. Consequently, the study of molecular dynamics in relation to biocatalysis can promote an optimization of experimental time in the bench, as these studies can predict the behavior of the biocatalyst against the substrate of interest (Barbosa et al., 2019;Brandão et al., 2020;Santana et al., 2020;Rodrigues et al., 2021). The study of tunnels to access and transport of ligands using experimental techniques in silico is not too investigated, not so well understood, and not yet widely used, as these analyses present a broad challenge, since most methods used to verify tunnels are based on molecular dynamics simulation practices . Therefore, the objective of this work was to identify the tunnels present in the biocatalysts and their characteristics, as well as to select the tunnels with the highest probability of substrate transport (in the acetylation reaction) in silico, carry out the transport of products in the tunnels, identify the most viable route of obtaining, and, therefore, identify the best biocatalyst to be used in the reaction to obtain α-tocopherol esters applying molecular dynamics tools.

MATERIALS AND METHODS
To carry out this study, a screening of lipases was carried out with potential applications for the biotransformation of α-tocopherol.
Altogether, 9 lipases were selected ( Table 1), and with the aid of the Discovery Studio program, the amino acids that make up the active site of the right lipases were selected. For each lipase, the PDBs (Protein Data Banks) were able to verify the crystallographic structures. The SMILES (Simplified Molecular Input Line Entry System) codes available at PubChem for each molecule subjected to analysis were used ( Table 2). Subsequently, the PDBs of each lipase were identified using the computer web server CaverWeb, together with the SMILES codes of the substances to analyze the possible tunnels contained in the structures as well as the possible transport of each molecule in the biocatalyst. The CaverWeb 1.0 server is available for free; it is simple to realize the analysis and allows the online execution of the steps of screening. First, the enzyme specifications were carried out followed by specifications; after which the calculation of tunnels according to energy following the screening of trajectories was carried out, and the data were analyzed followed by the results .
After that, the structure of each lipase, reaction analysis, transport of each molecule in the biocatalyst, behavior of ligands in the tunnels, and E max and E a as well as bottleneck radius (Å), length (Å), curvature, and throughput were analyzed according to the web server CaverWeb. Thus, to understand in a computational way the behavior of ligands in the tunnels present in the lipase structures is of great interest 2 | SMILES codes for binders and products, molecular formula, and structure of the substances applied in this study.

Substance
Canonical SMILES Molecular formula Frontiers in Bioengineering and Biotechnology | www.frontiersin.org January 2022 | Volume 9 | Article 805059 4 for biocatalysis because it allows the understanding of the reaction possibilities between the substrate and enzyme. The simulation was started with the identification of the number of tunnels and their properties. Then, the SMILES codes of the respective ligands were added, and the server calculated the energies [E bound , E max , E surface , E a , and ΔEBS (E surface ) (Kcal/ mol)] for the screening of trajectories in the tunnels. Based on this observation, data analysis was performed, and thus we predicted which tunnel was most likely to transport the substrates and products .
The energy expenditure of each path is transformed into a new measure of the importance of a direction which is known as throughput. The values obtained by the transfer rate can vary between 0 and 1, and the closer the rate to 1, the greater the probability of transporting binders in the tunnel. However, tunnels that have a lower throughput have lower possibilities of transporting (Chovancova et al., 2012;Vavra et al., 2019). The possibilities of transporting binders in tunnels can be analyzed using E max (highest binding energy in the path) and E a (activation energy of the E max association-E surface for binders) and that the lower the values of E max and E a , closer to 0, the greater the probability that a ligand has to go through the tunnel to reach the active site. Thus, to identify which tunnel is most likely to be transported, the morphological characteristics of the tunnels obtained in CaverWeb were observed along with the E a and E max energy values Vavra et al., 2019;Singh et al., 2021). A scheme of steps about this computational analysis is visualized in Scheme 1. SCHEME 1 | Representative scheme of the steps that were performed in the computational analysis.
Step 1: Obtaining the PDB's of the analyzed enzymes in the Protein Data Bank.
Step 2: Structuring and identifying the enzymatic structure as well as identifying the active site of each enzyme.
Step 3: Search and selection of the binders and products (SMILES codes) of this analysis.
Step 4: Inserting the SMILES codes of the substances into the analyzed enzymes using the CaverWeb server.
FIGURE 1 | Structural representation of the amino acids present in the CALB access tunnels. Tunnel 1 in blue, tunnel 2 in green and tunnel 3 in red taken from CaverWeb (https://loschmidt.chemi.muni.cz/caverweb/) with PDB 5A6V.

RESULTS
In in silico analysis of biomolecules, mainly lipases can be applied to verify enzymatic behavior under optimal conditions not including interference from temperature changes (Table 1), pressure, or immobilization techniques that can promote the best enzymatic performance in addition to being a low tool cost and easy access for the scientific community. Molecular docking is a tool used to estimate the chemical and biological activities of molecules, employing computational calculations prior to experimental processes and that demonstrates the energy parameters of the ligand molecules. The numerical values of these parameters are applied to compare the biological activities of mottraslecules against interactions between molecules and enzymes (Taslimi et al., 2020). Among the most different applications of molecular dynamics, the study of enzymatic tunnels and their analysis of substance transport can be highlighted for a better understanding of biocatalytic processes.

Tunnel Analysis
The structural information that could be observed through CaverWeb was the binder inlet diameter (Å), tunnel length (Å), curvature, and the possibility of transporting through each tunnel present. This information is important to molecular dynamics, as it demonstrates the possibility of observing which pathway is most likely to enter the ligand to pass through the lipase structure and reach the active site to carry out the reaction of interest as well as demonstrate along the way that the binder can go through each tunnel. Each enzymatic tunnel has its own characteristic (Figure 1), not being the same tunnel with the same morphological or functional structure in the most different existing enzymatic structures. Thus, the conformational study of tunnels present in enzymes has great relevance to biocatalytic reactions, as it promotes the optimization of several enzymatic properties .
To analyze the transport possibilities of the tunnels, it is necessary to understand the information about the energies that are demonstrated in CaverWeb ( Figure 2). The energy values [E bound , E max , E surface , E a , and ΔEBS (kcal/mol)] that were obtained through the in silico analysis of the binders (α-tocopherol, vinyl acetate, vinyl succinate, vinyl nicotinate, and vinyl ferulate) and donated products (α-tocopherol acetate, α-tocopherol succinate, α-tocopherol ferulate, and α-tocopherol nicotinate) were verified about the tunnels present in the enzymatic structures individually ( Figure 3). The lower absolute affinity energy (E a ) value indicates that the binder has greater ease of substrate entry. ΔEBS corresponds to the energy difference between the active site present in the enzyme and its surface. Minimum bonding values, corresponding to the ligand that is connected to the active site and the surface, are related to the ligand that is connected to the surface Stourac et al., 2019). Substance transport analyses were performed with the 9 selected lipases. Consequently, it can be verified that Thermomyces lanuginosus lipase has 8 tunnels, among them only 2 were with high ligand transport probabilities that α-tocopherol presented a higher probability of entry through tunnel 1 with FIGURE 3 | Graphical representation of the CALB crystallographic structure (PDB: 5A6V) in gray, in yellow the location of the active site, and the visualization of the 6 access tunnels in different colors that compose it, obtained in Caver Web and visualized in Discovery Studio, and the characteristics of tunnel size (Å), Length (Å), curvature, and the throughtput (possibility of transport).
Frontiers in Bioengineering and Biotechnology | www.frontiersin.org January 2022 | Volume 9 | Article 805059 (kcal/mol)]. Human gastric lipase has 7 tunnels in its structure; however, only 2 were with greater possibilities for transporting substances. α-Tocopherol demonstrated entry into tunnels 1 and 2 with E max 11.4 (kcal/mol)and E a 8. Thus, understanding the behavior of a ligand regarding the enzymatic structure is important to analyze the path that these molecules take to perform the biocatalysis. Figure 4 represents the conformational structure of CALB where the 3 tunnels (analyzed in the CaverWeb) were with the highest probability of transporting substances, confirming the experimental results obtained by Torres et al. (2008). CALB is a biocatalyst that realizes great catalytic reactions such as esterification (Santin et al., 2014), transesterification (Ferrer et al., 2005), acetylation (Torres et al., 2008), and others can be observed due to their morphological characteristics. To understand how the transport of molecules through enzymatic tunnels works, it is important to study molecules and substances through enzymatic tunnels.

Substrate Transport
To catalyze a substrate, it is necessary to understand the following points: the substrate must fit in the binding tunnel present in the enzyme, and the substrate must be able to pass through the tunnel. In enzymatic structures, the active site is located in the internal cavity, while the tunnels are connections between the active site and the substrate (Wang et al., 2021). Lipases are a family of enzymes that contain a tunnel region that connects the acyl end of their substrates. They have a length of 8-22 Å and a geometry that can vary greatly between different lipases. By sketching the tunnel wrap residues, specificity for the length of the acyl chain can be achieved .
To consist how substrate transport happens, the ligand is transported toward or away from the active site within the receptor through a tunnel. However, each receiver can have multiple tunnels that lead to the same active location. In CaverWeb, the active site is considered as the beginning of the tunnel and the surface of the receiver as the end of the tunnel. The calculation of the tunnels starts by selecting the starting point at the active site through energies, where the minimum size of the probe defines the minimum radius that the tunnel must have to be found. The depth and radius of the tunnel define the surface of the substrate in the receiver (Hozzová et al., 2021).
For the energetic understanding regarding the substrates submitted to CALB, Table 3 shows the input energies of the substrates in the main access tunnels. For α-tocopherol, the best entry tunnel was tunnel 3 with E max −4.6 and E a 1.1. Vinyl acetate showed very close entry possibilities in the 3 tunnels (tunnel 1: E max −2.4 and E a 0.1; tunnel 2: E max −2.6 and E a 0; and tunnel 3: E max −2.5 and E a 0.1) which indicates that this ligand has different input possibilities.
To comprehend analysis of the transport of substances or ligands in a specific tunnel, as we can see in Figure 5, the representation of the energy trajectory of α-tocopherol through the tunnel with the highest transport probability was represented by tunnel 3. In Figure 6, the vinyl acetate energy trajectory was through the tunnel with the highest probability of transport (tunnel 1). In both images, the direction of the protein surface as well as the direction of the binding site can be seen. α-Tocopherol, through its structural conformation in the best transport probability tunnel, as well as the interactions with the amino acids that are disposed there can be visualized in Figures 7A,B.
To get the interactions of vinyl acetate with the protein (Table 4), it is observed that the existing interactions are van  Frontiers in Bioengineering and Biotechnology | www.frontiersin.org January 2022 | Volume 9 | Article 805059 9 der Waals interactions, which demonstrate an interaction, even if weak, with the amino acids Leu 140, Ile 189, and Glu 188. It is important to highlight that van der Waals interactions play an important role in structure formation, energetic stability, and reaction mechanisms for a wide variety of molecules (Kleshchonok and Tkatchenko, 2018).
Therefore, obtaining knowledge of the transport of ligands in the access tunnels of a protein structure as well as understanding the types of bonds that the ligands have with the respective amino acids that are disposed of there is of fundamental importance.

Product Transport
Through the energetic information of the substrates that was obtained, the characteristics of the 3 tunnels present the greatest possibilities of the output of the product (α-tocopherol esters). The information presented regarding the conformational study of the tunnels present in CALB is considered of primordial importance, as it is through these characteristics that we can understand the energetic behavior of binders and products.
However, for the products analyzed in this research, the 4 products showed transport potential through the first 3 tunnels located in CALB but with a greater possibility of exit through tunnel 2 ( Figure 8) which were α-tocopherol acetate with E max −3.7 kcal/mol and E a 2 kcal/mol, α-ferulate tocopherol with E max −4.7 kcal/mol and E a 1.9 kcal/mol, and α-tocopherol nicotinate with E max −4.7 kcal/mol and E a 1.9 kcal/mol and through tunnel 3 was α-tocopherol succinate with E max −4.8 kcal/mol and E a 1.8 kcal/mol. It is important to note that although α-tocopherol acetate is the most commonly used ester in the pharmaceutical and cosmetic industry, α-tocopherol succinate, α-tocopherol nicotinate, and α-tocopherol ferulate are of fundamental importance for biological studies and biotransformation, as they are also addressed in the literature (Xin et al., 2011;Jiaojiao et al., 2021).
Molecular dynamics tools can perform an orientation that obtains a higher probability of one or more molecules in the active sites of the lipase. As molecular docking will generate a series of possible complexes, the best orientation will be selected through scoring functions and used to predict binding affinity. However, it is clear that there are different tunnels present in the lipase, and Table 5 shows the summary of data for the tunnels in the target lipase. Based on information about the structure of biomolecules and prior knowledge of ligands, the application of molecular dynamics tools to predict their behavior and optimize future reaction processes is of fundamental importance.
To analyze the output of the product of interest, α-tocopherol acetate passes through an energy trajectory that can be seen in Figure 8. This representation demonstrates the energy profile through the molecular chain in tunnel 2, which is the tunnel with great possibility of passage of this substance. Figures 9A,B show the position of α-tocopherol acetate about the tunnel and the receptor (CALB) and the disposition of this substance about the exit tunnel of the product as well as the interactions of this molecule with the amino acids present in the enzymatic structure of CALB.

Substance Amino acid Interaction
Vinyl acetate Leu 140 van der Waals Ile 189 Glu 188 FIGURE 8 | Representation of the trajectory binding energies and ray size for α-tocopherol acetate in tunnel 2.
Frontiers in Bioengineering and Biotechnology | www.frontiersin.org January 2022 | Volume 9 | Article 805059 Table 6 shows the amino acids that interact with α-tocopherol acetate as well as the types of interaction that are related to this molecule. Amino acids Ala 279, Ala 282, Glu 188, Leu 140, Thr 138, and Thr 40 have van der Waals-type interactions. It can be verified that the interactions that exist are van der Waals interactions, which express an interaction of the substance, even if weak, with the amino acids present in the enzyme. These interactions play an important role in the formation of structures, energy stability, and reaction mechanisms for a wide variety of molecules (KLESHCHONOK; TKATCHENKO., 2018). The amino acid Leu 278, on the other hand, has a pi-alkyl-type interaction, and the amino acids Val 154, Ile 285, and Leu 144 perform alkyl-type interactions. Since van der Waals-type interactions are weaker interactions, the pi-alkyl and alkyl interactions are characterized by hydrophobic bonds. Figure 10 represents how the entry of ligands (α-tocopherol and vinyl acetate) occurs to obtain α-tocopherol acetate. The characteristics of CALB tunnels demonstrated in CaverWeb represent a particularity, which allows this lipase to have a significant reaction potential compared to other analyzed lipases.

Potential of Computational Analysis
Tunnel 1 has shorter lengths (Å) (4) and curvature (Å) (1); however, it has a higher efficiency (0.86). Tunnels 2 and 3, on the other hand, have the same throughput (0.79) and the same curvature (1.7), but tunnel 3 has a greater length (Å) (11.1) than that of tunnel 2 (10.6). The neck radius Å for the 3 tunnels was similar. What stands out in this analysis was the verification of the condescending behavior of the binders in the tunnels. Note that α-tocopherol enters through tunnel 3, vinyl acetate through tunnel 1, and the product (α-tocopherol acetate) exits through tunnel 2. This implies that demonstrating the passage of these substances through different tunnels denotes ease of entry of binders and exit from the product of interest.

Tunnel Analysis
Access tunnels provide information about different important functions of enzymes, such as the distinction between substrates with higher transport probability, the control of  the entry of co-substrates or solvent, and the observation of possible enzymatic damage caused by transition metals. as well as the prevention of cell damage by highly reactive intermediates and metabolites, the simultaneous order of individual catalytic steps, synchronization of reactions that require a large number of substrates, intermediates and/or cofactors and, finally, the spatial safety concentration of the reactant species . This knowledge provides the obtention of transport interest (substrates and products) on tunnels. According to Park et al. (2013), the substrate entrance is delimited by hydrophobic residues, but unlike other lipases, CALB does not have a cap and does not undergo interfacial activation. Four hydrophobic residues, Ile 189, Leu 278, Ala 282, and Ile 285, form a narrow tunnel that must be traversed by the substrate before the formation of the second tetrahedral intermediate. These amino acids can be observed in tunnels 1, 2, and 3 that make up the CALB represented in Figure 1.
Most tunnels function as extensions of the active site of an enzyme, which are characterized in their structural conformation by an entrance and exit, and also in the functioning of a natural connection between substances. The conformational characteristic and size of a tunnel can not only predict the maximum size that a ligand can pass but also inform the number of interactions that can be performed or the barriers that a ligand can face when being transported to the tunnel active location in an energy process (Kingsley and Lill, 2015).
Molecular dynamics analyses are used to calculate the binding affinities of the biomolecule to the ligand. These methods use an energetic approach to calculate the relative free energy binding between a receptor and two ligands based on the thermodynamic cycle between the molecules . Thus, analyzing individually and particularly the structure of each tunnel present in the enzymatic structure (Figures 3, 4) is a fundamental importance to understand the transport possibilities that a ligand is subjected to energetically.
The active site is a specific space that is usually hidden in the protein, since the chemical reaction between the protein and the ligand can take place in this area. Consequently, the importance of access tunnels for these ligands is directly linked to the amino acids that are located around the tunnel, which substantially enable the structural and kinetic stability of the studied protein (Vonásek et al., 2019). Thus, most of the non-active sites reported so far are located in the flexible loop regions. According to their functionalities, the positions of non-active sites can be divided as follows: protein surface, region of turns, and the substrate tunnel, and the latter being considered an alternative approach to understand catalyst mechanisms and enzymatic engineering (Wang et al., 2021). Therefore, studying the transport of ligands in enzyme tunnels is of fundamental importance for biocatalysis.
FIGURE 10 | Representative diagram of the 3 most promising tunnels for transport in CALB, in pink (tunnel 1), blue (tunnel 2), and green (tunnel 3) and the demonstration of binder entry and exit of the product.  (Kleshchonok and Tkatchenko, 2018;Ranjbar et al., 2020;El Aissouq et al., 2021). According to Kumaresan et al. (2011), the activity and enantioselectivity of CALB are directly influenced by the interaction of the substrate with three binding pockets that are located close to the active site, being differentiated into 3 types with different constituent amino acids. Among them are medium pocket (G39, T42, S47, W104, and A225), large pocket (A141,  L144, V154, I285, P289, and K290), and acyl (D134, T138, Q157,  I189, and V190). From this understanding, it is possible to analyze some amino acids that α-tocopherol interacts with corresponding to the large pocket (Ala 141, Leu 144, Val 154, and Ile 285) as well as in the acyl pocket with the amino acids Ile 189 and Thr 138. Ferrari et al. (2014) carried out a study based on the backbone protein and the other based on the oleoyl chain and the backbone residues (CALB) that constitute the catalytic cavity. Among them can be observed the amino acids: Thr 138, Ile 189, Val 190, Asp 134, Thr 40, Gln 157, Ile 285, Gly 44, and Val 149 as well as the residues that constitute the oxyanion hole: Thr 40 and Gln 106. Kim et al. (2014) identified that of all residues present in CALB, two isoleucines, ILE-189 and ILE-285, play a key role in the openclosed dynamics of the cavity and can hinder the passage of the substrate to the catalytic site. Consequently, it is believed that the positions of the amino acids that form the biocatalyst in its cavity influence the enzymatic activity, since the CALB catalytic cavity is divided into an acyl and alcohol group next to the ILE close to the cavity, and this conformation can affect the enzyme activity.
Therefore, performing the analysis of the ligands in the access tunnels of the substrates in the cavities present in the biocatalyst in relation to the amino acids present is important for the understanding of biocatalytic reactions. In Figures 11A,B, it is possible to observe the transport of vinyl acetate corresponding to tunnel 1 of CALB. Vinyl acetate is an acylating agent which undergoes nucleophilic attack by the amino group of an amino acid or peptide fragment that can result in the formation of a peptide bond between them (Machado et al., 2004;Barbosa et al., 2011).

Product Transport
To analyze the functioning of the transport of substances about the interaction with the active site, Chaput et al. (2012)   For the output of CALB products, the best tunnel to exit α-tocopherol acetate (E max −3.7 kcal/mol and E a 2 kcal/mol) (Figure 6), α-tocopherol succinate (E max −4.8 kcal/mol and E a 1.8 kcal/mol), and α-tocopherol nicotinate (E max −4 .5 kcal/mol and E a 1.3 kcal/mol) is tunnel 3. Through the binders and products on the tunnels present in the CALB, it can be seen that the substrates' entrance and exit tunnels are distinct, which denotes an optimized possibility of obtaining reaction products.

Potential of Computational Analysis
For better knowledge of this mechanism, Kyriakou et al. (2012) carried out in their work a in silico study of how quercetin (a type of flavonoid) behaves in CALB to perform an acylation of this component. It can then be seen that this flavonoid can be favorably accommodated for transformation into the substratebinding site of CALB.
Consequently, product formation was observed in situ with a yield of 45% only in the presence of the enzyme and optimal reaction conditions (solvent and acyl donor). The great interest also for biocatalytic processes is to consider the importance and high interest of biochemistry in given industrial applications. Molecular dynamics tools can be applied to shorten times in search of new ecologically correct routes; some studies using CALB are already known for their high efficiency and high selectivity applied in different industrial processes or suitable catalytic activities to enable biotransformation (Idris and Bukhari, 2012;Duncan and Suzuki, 2017). Thus, the biotransformation using CALB as a biocatalyst to obtain α-tocopherol acetate is already known, and the molecular dynamics tool showed that it was possible to understand the behavior of this biocatalyst through the substrate.

CONCLUSION
Through the computational screening performed in this work, it can be verified that the number of existing tunnels in the enzymatic structure does not imply a better biocatalytic result, since the lipase with the largest number of tunnels submitted to in silico analysis in this work was the LPP with 39 tunnels, and no significant results were obtained. On the other hand, the lipase that best obtained the results was CALB (with 6 tunnels) because the characteristics of the tunnels present in its structure, compared to the binders submitted for analysis, were more likely to transport binders because the net tunnel for substrates was distinct from the product exit tunnel (α-tocopherol acetate). However, compared to results described in the literature using CALB as a biocatalyst to obtain α-tocopherol acetate, it was reinforced as well as other biocatalysts that were used (LBC and LCR) did not obtain promising results. Therefore, the in silico analysis of biocatalytic reactions has shown significant results in analyses before in situ reactions. It is important to point out that when performing the in silico analysis, it demonstrates optimization of the execution time in bench experiments, as it performs a prediction of the submitted biocatalytic behavior.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
TA and LS worked on the methodology used and investigation carried out for the study and wrote the original draft of the manuscript. EF carried out the investigation and validation for the study. MP worked on the conceptualization, methodology, and investigation for the study and helped with the writing, reviewing, and editing of the manuscript. CF contributed to the conceptualization, project administration, funding acquisition, supervision, and writing, reviewing, and editing of the manuscript.

FUNDING
The authors are grateful for the financial support from the following Brazilian research funding agencies: Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) with finance code 001 and process number 88882.365586/2019-01; Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq); and Fundação de Apoio à Pesquisa e à Inovação Tecnológica do Estado de Sergipe (FAPITEC/SE).