A Novel Fungal Lipase With Methanol Tolerance and Preference for Macaw Palm Oil

Macaw palm is a highly oil-producing plant, which presents high contents of free fatty acids, being a promising feedstock for biofuel production. The current chemical routes are costly and complex, involving highly harsh industrial conditions. Enzymatic processing is a potential alternative; however, it is hampered by the scarce knowledge on biocatalysts adapted to this acidic feedstock. This work describes a novel lipase isolated from the thermophilic fungus Rasamsonia emersonii (ReLip), which tolerates extreme conditions such as the presence of methanol, high temperatures, and acidic medium. Among the tested feedstocks, the enzyme showed the highest preference for macaw palm oil, producing a hydrolyzate with a final free fatty acid content of 92%. Crystallographic studies revealed a closed conformation of the helical amphipathic lid that typically undergoes conformational changes in a mechanism of interfacial activation. Such conformation of the lid is stabilized by a salt bridge, not observed in other structurally characterized homologs, which is likely involved in the tolerance to organic solvents. Moreover, the lack of conservation of the aromatic cluster IxxWxxxxxF in the lid of ReLip with the natural mutation of the phenylalanine by an alanine might be correlated with the preference of short acyl chains, although preserving catalytic activity on insoluble substrates. In addition, the presence of five acidic amino acids in the lid of ReLip, a rare property reported in other lipases, may have contributed to its ability to tolerate and be effective in acidic environments. Therefore, our work describes a new fungal biocatalyst capable of efficiently hydrolyzing macaw oil, an attractive feedstock for the production of “drop-in” biofuels, with high desirable feature for industrial conditions such as thermal and methanol tolerance, and optimum acidic pH. Moreover, the crystallographic structure was elucidated, providing a structural basis for the enzyme substrate preference and tolerance to organic solvents.


INTRODUCTION
The increasing demand for energy that complies with the current restrictions on emissions of polluting gases is promoting the development and production of alternative fuels from sustainable and renewable sources (Yan et al., 2015;Wang et al., 2018). The biohydrocarbons, in particular, medium-and long-chain fatty alkanes and alkenes, represent very promising alternatives for replacing the fossil fuels used in the transportation sector, once they have similar chemical composition and physical characteristics to those of conventional fuels derived from petroleum (Xu et al., 2017). Therefore, biohydrocarbons can share the infrastructure for storage and transportation already used for the distribution of gasoline, diesel, and jet fuel (Yan et al., 2015;Sousa et al., 2018) and can be combined at any proportion with those fuels without any modification in the engines and other mechanical components of the vehicles (Sousa et al., 2018).
Macaw palm (Acrocomia aculeate) is a native species distributed in the tropical and subtropical Americas, and it occurs naturally from southern Mexico to northern Argentina (Ciconini et al., 2013;Pires et al., 2013;Silva et al., 2014;da Conceição et al., 2017). It is predicted that a commercial plantation under appropriate agronomic conditions can produce from 16,000 to 25,000 kg of fruit per hectare (Pires et al., 2013). Its fruits are oil rich, accumulating up to 70% of oil (dry weight) and yielding approximately 6,200 kg of oil per hectare (Pires et al., 2013;Silva et al., 2014). The macaw palm has interesting agronomic and ecological features because it can occupy degraded areas or agroforestry systems given that it has high plasticity to grow in different ecosystems (Moura et al., 2009;Pires et al., 2013), avoiding conflict with areas used for food production. Macaw pulp oil presents a high content of free fatty acids (FFA) and, thus, elevated acidity (Lopes et al., 2013), with oleic acid (C18:1) being the predominant fatty acid (Bora and Rocha, 2004). Thus, it can generate a high-quality biofuel, with high amounts of monounsaturated compounds (Aguieiras et al., 2014).
Lipases (triacylglycerol acylhydrolases, E.C. 3.1.1.3) are watersoluble enzymes (Reis et al., 2009) that can catalyze oil hydrolysis, esterification, and transesterification (Goswami et al., 2013). Most of them are interfacial enzymes, which means that, at the oilwater interface, they undergo structural changes that involve a movement of a short helical segment (lid) that covers the active site. Thus, during the activation process, the lid changes its conformation, exposing the catalytic triad and increasing the aliphatic surface surrounding the active site, required for substrate binding (Derewenda et al., 1992;Maruyama et al., 2000;Gruber and Pleiss, 2012). In addition, due to their transesterification activity, during the last decade, lipases have been extensively studied as a promising industrial biocatalyst for biodiesel production (Bajaj et al., 2010;Hwang et al., 2014;Norjannah et al., 2016;Amini et al., 2017). It has recently been reported that few vegetable lipases, especially those from castor bean seeds, are more effective on the hydrolysis of macaw oil because they usually act at acidic pHs compared to fungal lipases that perform under more basic conditions (Ory et al., 1962;Muto and Beevers, 1974;Avelar et al., 2013;Bressani et al., 2015).
Therefore, considering the interest to make macaw palm an alternative platform to produce renewable biofuel and other bioproducts, and the very limited knowledge of biocatalysts available that can convert this specific acidic oil, the aim of this work was to develop a fungal lipase that could be compatible and efficient for this feedstock. Rasamsonia emersonii is an important microbial platform for cellulose degradation since its genome possesses many effective glycoside hydrolases (Gudmundsson et al., 2016;Martínez et al., 2016); nevertheless, it remains underexplored regarding enzymes active on triglycerides. For the first time, we have cloned and heterologously expressed the lipase from Rasamsonia emersonii (ReLip). The functional analysis revealed that ReLip has an optimum acidic pH profile and was active on a wide range of substrates from short to long carbon length chains. The enzyme was remarkably active on macaw oil, being able to produce a total of 92% FFA content. Furthermore, ReLip is highly thermo-tolerant under extreme conditions such as in the presence of methanol. Crystallographic analysis revealed the closed conformation of this lipase indicating structural properties driving substrate selectivity and tolerance to organic solvents. In this sense, this work, besides providing a novel effective biocatalyst suitable for typical harsh industrial conditions, reveals new molecular aspects associated with substrate preference and tolerance to extreme conditions.

MATERIALS AND METHODS
The lipase gene from Rasamsonia emersonii cloned into the pET28a(+) vector with a histidine-tag at the N-terminus was purchased from Genscript (Piscataway, NJ, United States) (GenBank accession number: 915129218). Gum Arabic and ethanol (minimum 99.5%) were purchased from Synth and acetonitrile (minimum 99.9%) from Tedia. Commercial soybean (Liza), canola (Liza), olive (Gallo), sunflower (Liza), corn (Liza), and coconut oils were obtained at local markets. Palm and macaw oils were kindly supplied by Agropalma and Dr. Carlos Colombo [Agronomic Institute of Campinas (IAC), respectively]. FFA composition and percentage of FFA (acidity) of the vegetable oils used in this work are described in Supplementary Table S1. The FFA content was measured according to AOCS standard number Ca 5a-40 (Walker, 1990). The fatty acid composition of each oil was determined using a Pegasus HT (Leco) connected to a gas chromatographer (7890A, Agilent), according AOCS standard Ce 1-62 method (AOCS, 2005).

Protein Expression and Purification
The plasmid ReLip1-pET28a was produced in E. coli BL21(DE3)pLysS strain with pRARE2 plasmid. The cells were cultured in Luria-Bertani (LB) agar containing 25 µg/ml of kanamycin and chloramphenicol. For growth, one colony was picked and cultured in 5 ml of LB broth containing 25 µg/ml of kanamycin and chloramphenicol, with overnight shaking at 37 • C and 250 rpm. After this period, the material was diluted 100-fold with fresh selective terrific broth (TB) and incubated at 37 • C until an optical density (OD 600nm ) of 1.2. Thus, the protein expression was subsequently induced by 0.5 mM isopropyl β-D-1-thiogalactopyroanoside (IPTG) and cultivated for 16 h at 20 • C and 250 rpm. The cells were collected by centrifugation at 8,000 × g for 20 min at 4 • C and resuspended in 20 ml of lysis buffer (50 mM phosphate pH 7.4, 300 mM NaCl, 30 mM imidazole, 1 mM phenylmethylsulfonyl fluoride (PMSF), and 0.5 mg ml −1 of lysozyme). The cells were, then, lysed by sonication, and the soluble protein was extracted by centrifugation at 12,000 × g for 30 min at 4 • C. Thereafter, the enzyme was purified from the supernatant by two chromatographic steps: (i) metal-affinity chromatography, using a non-linear imidazole gradient (20-500 mM) in a 5-ml HisTrap HP column (GE Healthcare), previously equilibrated with 50 mM phosphate buffer, pH 7.4, 20 mM imidazole, and 300 mM NaCl and (ii) size-exclusion chromatography using a Superdex 200 HiLoad 16/60 column (GE Healthcare), previously equilibrated with 50 mM phosphate buffer, pH 7.4 and 150 mM NaCl coupled to an AKTÄ Purifier FPLC system (GE Healthcare Life Sciences). All steps of heterologous expression and enzyme purification were analyzed by sodium dodecyl sulfate-polyacrylamide gel (12% polyacrylamide) electrophoresis (SDS-PAGE). The protein concentration was spectroscopically determined using the molar extinction coefficient calculated from the amino acid composition 1 , which is 38,765 M −1 cm −1 .
The optimum temperature and pH of the enzyme was determined spectrophotometrically following the hydrolysis of p-nitrophenylbutyrate (pNPB; C4:0; Sigma-Aldrich Co., St. Louis, United States) at 348 nm. The reactions were performed as described above. The thermostability and the optimal temperature of the ReLip were evaluated in a temperature range of 35 • C to 75 • C, at pH 4.0. The optimum pH for ReLip activity was determined in the pH range of 2.0 to 6.0 at 65 • C, using the following buffers: HCl-glycine (pH 2.0 and 2.5) and citrate (pH 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, and 6.0), at a final concentration of 40 mM.
Each additive was investigated in individual 100-µl reactions containing 5 mM substrate pNPB, 100 mM citrate-phosphate buffer (pH 3.5), and 0.5 g/L of enzyme. The samples were incubated at 65 • C for 20 min, and then, 100 µl of pure acetonitrile was added to stop the reaction. The activity was determined by absorbance measurement at 348 nm and compared to the control, which was incubated without any compound. The activity measured without additives was defined as 100%. Then, the evaluation was expressed as relative activity (%).

Circular Dichroism (CD) Analysis
Circular dichroism (CD) measurements were carried out in a Jasco J-815 spectropolarimeter (Jasco, United States), using a Peltier-type temperature control system for temperature maintenance inside the cell. For the analysis of the secondary structure, the purified and homogeneous enzyme was diluted to multiple buffers in different conditions. Data were collected from 260 to 190 nm at 20 • C, using a N 2 flow rate of 10 ml/min, quartz cuvette with 0.1 cm of path length, scan speed of 50 nm/min, response time of 1.0 s, spectral bandwidth of 1.0 nm and spectral resolution of 0.1 nm. The final CD spectrum obtained was an average of 20 accumulations.
The CD spectroscopy was also applied to determine the protein stability at denaturing conditions, by measuring thermalinduced unfolding and melting temperature. For this purpose, the sample was heated from 20 • C to 100 • C at a rate of 1.0 • C/min. The reversibility of the temperature effect was evaluated by cooling the denatured sample from 100 • C to 20 • C, using the same parameters described above. The values obtained in CD measurements (mDeg) were normalized to residual molar ellipticity (MRE), and the temperature at the midpoint transition (T m ) was obtained by fitting the CD data to a sigmoidal Boltzmann function. Data were treated using the software Origin 8.1 (OriginLAB Corporation).

Determination of Hydrolytic Activity Using Triacylglycerol Substrates
The hydrolytic activity of the ReLip was determined on the hydrolysis of emulsified vegetable oils (Soares et al., 1999), aiming at verifying the optimum temperature and pH using triacylglycerol substrates and assessing the substrate specificity of the lipase. For the substrate specificity tests, vegetable oils such as soybean, canola, olive, sunflower, palm, macaw, corn, and coconut were used. The substrates were prepared by the emulsion of 7.5 g of vegetable oil, 67.5 g of water, 6.57 g of gum arabic (10% of water weight), and 60 ml of citrate buffer (pH 4.0, 50 mM), resulting in a final oil concentration of 5%. For the experiments, protein (at a final concentration of 2 g/L) and substrate were added (final reaction volume of 10 ml) into 125-ml Erlenmeyer flasks and incubated at 60 • C and 200 rpm for 5 min in an incubator shaker (Innova 44, New Brunswick Scientific, United States). After the incubation period, the reactions were ended by the addition of 10 ml of ethanol P.A. (Synth).
For the analysis of optimum temperature and pH and the design of experiments, the hydrolysis reactions were performed using pulp macaw oil. The emulsion was prepared as described above. Thus, substrate and protein, at pre-defined proportions and final volume of 10 ml, were added into 125-ml Erlenmeyer flasks and incubated at 35 • C and 200 rpm. After the desired time of reaction, the experiments were stopped by the addition of 10 ml of acetonitrile. The fatty acids produced in the hydrolysis were titrated with KOH solution, 0.2 mol L −1 , using phenolphthalein as indicator. The hydrolysis degree was calculated according to Eq. (1) (Rooney and Weatherley, 2001): (1) where V is the volume of potassium hydroxide solution required in the titration process (ml), M KOH is the KOH solution concentration (mol/L), MM FFA is the average molecular mass of fatty acids of the macaw oil (276.01 g/mol), w is the weight of the sample titrated (g), and f is the fraction of oil in the begging of the reaction.

Design of Experiments and Optimization
The experiments were performed at 35 • C and pH 4.0. A matrix of 12 experiments were performed in duplicate, using alpha for orthogonality (α) of 1.21 and four replications at the center point. The independent parameters were defined as follows: catalyst concentration (C) and time of reaction (t). Supplementary  Table S2 shows the variables used in the central composite design (CCD), with their five coded and uncoded levels. The response variable was defined as the hydrolysis degree. The hydrolysis degree obtained were fitted in a quadratic model using regression analysis, and for selecting the significant terms of the model, parameters with values of p < 0.05 were considered statistically significant. The range of each parameter was selected according to preliminary tests. For the optimization process, canonical analysis technique was employed (Box and Tiao, 1977). Thus, aiming at validating the optimization, the enzymatic hydrolysis was performed under the optimal experimental conditions suggested by the analysis. In this work, the CCD and RSM were developed by Statistica software version 7.0 (Statsoft, United States), and the canonical analysis was implemented using the software Maple 17.

Dynamic Light Scattering (DLS)
Dynamic light scattering (DLS) technique was used to determine the hydrodynamic behavior of ReLip in buffer purification at different concentrations (1-10 mg/ml). The experiments were carried out at room temperature using a Malvern Zetasizer Nano ZS90 (Malvern Instruments, Worcestershire, United Kingdom) with a 633-nm laser, in a quartz cell with a scattering angle of 90 • . The diffusion coefficient (DT) was determined from the analysis of measured time-dependent fluctuations in the scattering intensity and used to calculate the hydrodynamic radius (Rh) of the protein according to the Stokes-Einstein equation.

Small-Angle X-Ray Scattering (SAXS)
Small-angle X-ray scattering measurements were acquired using a monochromatic X-ray beam (λ = 1.488 Å) from the D01A-SAXS2 beamline at the Brazilian Synchrotron Light Laboratory (LNLS, Brazil). ReLip (2-4 mg/ml) samples were prepared in buffer purification. Prior to conducting the SAXS experiments, all samples were centrifuged for 15 min at 20,000 × g and 4 • C to remove any potential residual aggregates. The sample-todetector distance was set as 1,000 mm, resulting in a scattering vector (q) range of 0.02 Å −1 < q < 0.50 Å −1 , where the q-vector magnitude is defined as q = 4πsinθ/λ in which 2θ is the scattering angle. Samples were analyzed at 20 • C and placed in 1-mm-path-length mica cells, and the scattering profiles were recorded in 10 successive frames (30 s each) to monitor radiation damage. The buffer contribution in each SAXS profile was subtracted taking into account the attenuation and integrated in the sample using the FIT2D software (Hammersley, 2016). Sample mono dispersity was checked by means of the Guinier's law, and all studied systems presented here were found to be monodisperse, and no aggregation took place over the SAXS curves (data not shown). The Gnom program (Svergun et al., 1988) was also used to generate the pair distance distribution function [p(r)] and the protein maximum dimension (D max ) from scattering profiles. Employing the p(r) function, the DAMMIN software (Svergun, 1999) was applied to obtain ab initio models for ReLip (dummy atom model) by a simulated annealing optimization routine that yields a best fit to the experimental scattering data. Shapes were reconstructed by averaging a minimum of 10 different ab initio models using the DAMAVER package (Volkov and Svergun, 1995). The experimentally derived low-resolution envelopes were superimposed on structures obtained by molecular modeling using the SUPCOMB software (Kozin and Svergun, 2001).

Crystallization, Data Collection, and Processing
Multiple amino acid sequence alignment of fungal lipase sequences was performed using Clustal Omega and ESPript 3.0. The following lipases were aligned: 1DT3:  Figure  S1). Protein sample at a concentration of 16.5 mg ml −1 in 50 mM sodium phosphate, pH 7.4, and 150 mM NaCl buffer was used in the crystallization experiments. The sample was under sitting-drops prepared at 18 • C using a Cartesian HoneyBee 963 system (Genomic Solutions). Conditions (544) were tested, based on commercially available crystallization kits from Hamptom Research (SaltRx, Crystal Screens I and II), Emerald BioSystems (Precipitant Synergy and Wizard I and II), and Qiagen/Nextal (PACT and JCSG+). The drop, composed of 0.7 µl of the protein solution and 0.3 µl of the condition, was equilibrated over the reservoir containing 80 µl of the respective solution. Crystal optimization was performed using a systematic grid in which PEG 8000 concentration (from 20% to 8%) was varied in function of PEG 400 (from 20% to 8%) in a solution of 100 mmol/L MgCl 2 and 100 mmol/L Tris buffer pH 8.5. Two crystals with approximate dimensions of 30 × 30 µm were obtained. The solution consisted of 100 mmol/L MgCl 2 , 20% PEG 8000, 20% PEG 400, 100 mmol/L Tris buffer pH 8.5. The crystallization conditions are summarized in Supplementary Table S3.
Diffraction data were collected at the BL12-2 beamline from the Stanford Synchrotron Radiation Lightsource (Menlo Park, CA, United States). Crystals were cooled in liquid nitrogen and kept under a nitrogen gas stream during the collection at 100 K. A total of 1,800 images were collected, and data were indexed, integrated, and scaled using XDS package (Kabsch, 2010). Pointless (Evans, 2006(Evans, , 2011 and Zanuda (Lebedev and Isupov, 2014) were used to select and validate the chosen space group. Data analysis was performed with Xtriage (Zwart et al., 2005;Adams et al., 2010), and detection of anisotropy was made using the diffraction anisotropy server (Strong et al., 2006). Supplementary Table S4 presents the parameters used in the data processing.

ReLip Presents an Uncommon Acidic pH-Activity Profile and Its Catalytic Activity Is Not Metal Dependent
ReLip was overexpressed in the soluble and stable form in E. coli BL21(DE3)plysS yielding 120 mg/L of TB medium. After steps of purification, fractioned samples were analyzed by SDS-PAGE (Supplementary Figure S2). ReLip showed a single band between 31 and 45 kDa after size exclusion chromatography, corresponding to the expected molecular weight of the construct with an N-terminal His-tag (31.75 kDa).
The effects of pH and temperature on ReLip activity were investigated using pNPB as substrate. Figure 1A shows the effect of temperature variation on the enzymatic activity, from 35 • C to 75 • C at pH 4.0. The optimum temperature profile revealed that ReLip possesses a high thermal stability, with a relative activity of 60% at lower temperatures and a maximum activity around 65 • C when analyzed using a synthetic substrate. In Figure 1B, it is possible to observe the effect of pH, varying from 2 to 6. According to the pH profile, ReLip works better at acidic conditions, with an optimum pH around 3.5. Based on these results, ReLip shows an interesting feature when compared to most of the industrial lipases from fungi. Usually, lipases, mainly from microbial origin, act in an alkaline spectrum (Papaparaskevas et al., 1992;Soares et al., 1999;Hiol et al., 2000;Saxena et al., 2003;Yu et al., 2007;Utsugi et al., 2009;Belhaj et al., 2011). Only a few lipases, like the lipase from Aspergillus niger, show such acidic profile (Romero et al., 2007).
Knowing that additives can improve the enzymatic activity of lipases (Priyanka et al., 2019), ReLip was tested in the presence of a wide range of additives, such as detergents, chelating and modifying agents, cations, anions, and gum arabic. The results are shown in Table 1. Among the salts tested, only NaCl increased the relative activity. Considering this result, an investigation using different concentrations of NaCl was carried out, aiming at better understanding the effect of ionic strength on enzymatic activity. In the presence of 50 mM NaCl, ReLip had its relative activity increased twofold with small increments in the activity up to 0.5 M. Inhibition effects were observed only at concentrations higher than 1 M. Calcium was also evaluated in the range from 0.01 to 2 M, since this divalent cation is usually an important cofactor for lipase activity (Shibata et al., 1998;Simons et al., 1999;Choi et al., 2005;Tayyab et al., 2011;Lan et al., 2015). Different from typical lipases, no significant improvement in the activity was observed with calcium, indicating that ReLip is not a metal-dependent enzyme. This was further supported by the fact that the chelating agent, EDTA, did not substantially affect the catalytic activity.
The presence of gum arabic and detergents was also tested ( Table 1), once those compounds are typically used in substrate emulsions to stabilize and improve emulsion quality (Glogauer et al., 2011). Regarding gum arabic, its presence did not considerably increase the activity or promote inhibition at any of the concentrations tested, indicating that this gum can be used in the preparation of substrate emulsions for the application of ReLip. Among the detergents, SDS, Tween 80% and Tween 20% decreased the activity of the enzyme, especially SDS, with a remaining relative activity of only 11%. However, Triton X-100 showed a great and significant effect as additive for ReLip activity. An increase in the concentration of this additive promoted an activation on the activity, achieving the highest relative activity of 177% at a concentration of 0.5%. The inhibition process only started at concentrations above 10%, possibly due to a denaturation process (Glogauer et al., 2011). Considering that Tween and Triton X-100 are non-ionic detergents (with low hydrophilic/lipophilic balance value), they probably did not FIGURE 1 | Effect of temperature (A) and pH (B) on the activity of the purified lipase from Rasamsonia emersonii (ReLip). The temperature effect was evaluated from 35 • C to 75 • C, at pH 4.0, while pH influence was assessed at 65 • C using different buffers (pH 2.0 to 6.0). (C) ReLip specificity against synthetic substrates (pNPs) with different chain lengths (65 • C, pH 3.5, incubated for 20 min); 100% of activity corresponds to 59.34 U/g. (D) Hydrolytic activity of ReLip using different vegetable oils (5 wt%, 60 • C, pH 4.0, incubated for 5 min); 100% of activity corresponds to 2,000.36 U/g. All measurements were conducted in triplicate. The relative activity (%) was calculated considering the maximum catalytic activity observed for the biological unit of the enzyme, and the specific activity was defined as the number of units per g of protein.
interact extensively with the surface of proteins. On the other hand, ionic detergents, such as SDS, may have non-specific interactions with the surface of the proteins, facilitating structure denaturation (Priyanka et al., 2019).

ReLip Is a Lipase With Preference for Macaw Oil as Substrate
In order to identify the specificity of ReLip, the enzymatic activity was evaluated as a function of different chain lengths of p-nitrophenyl substrates ( Figure 1C). The results showed that ReLip was able to hydrolyze all substrates tested. However, ReLip preferentially hydrolyzed short acyl chains, showing higher activities against pNP butyrate (C4:0), pNP octanoate (C8:0), and pNP acetate (C2:0), with relative activities of 100%, 90.9%, and 80.1%, respectively. The lowest relative activity was observed against pNP palmitate (C16:0). It is known that one of the main differences between esterase and lipase is the fact that while esterase acts only on soluble substrates with smaller chains, lipases can hydrolyze insoluble substrates with larger chains above 10 carbons. The assays in the presence of pNPs with different sizes revealed that ReLip has the ability to act on both short and long chains and consequently, ReLip can be classified as an enzyme with esterase and lipase activity (Sharma et al., 2001(Sharma et al., , 2013Jaeger et al., 2002;Silva et al., 2009). Thus, in order to prove the activity on insoluble substrates and confirm that ReLip is a true lipase, the hydrolytic activity using different emulsified vegetable oils, such as soybean, canola, olive, sunflower, palm, macaw, corn, and coconut, was evaluated ( Figure 1D). It is possible to observe that ReLip was active against all oils tested, presenting a sequential preference for macaw oil (2,000 U/g), followed by olive oil (811 U/g), and sunflower oil (788 U/g). These results showed that the hydrolytic activity depends on a wide range of factors and not only on the length of the acyl chain and the amount of unsaturation. Probably, the high amounts of FFA in the macaw oil led to an acidification of the medium that benefited the hydrolytic activity of ReLip.
In view of the preference for hydrolyzing macaw oil, the optimum temperature and pH profiles using this oil as substrate were also evaluated (Supplementary Figure S3). The enzyme showed a distinct temperature dependence using the macaw oil as substrate compared to the synthetic substrate (pNPB), with an activity plateau in the temperature range from 25 • C to 50 • C (Supplementary Figure S3A). Despite that, the enzyme retained up to 73.3% of its activity after 24 h of incubation in this temperature range, indicating a high stability under these conditions (Supplementary Figure S3C). Regarding the pH profile, both synthetic and natural substrates exhibited the optimum pH at 4.0 (Supplementary Figure S3B). Thus, to obtain the kinetic profile of macaw oil hydrolysis, the experiments were set at 35 • C and pH 4.0. The data revealed a hydrolysis degree of 50% with 3 h of reaction, and 63% was obtained with 8 h (Figure 2A). Consequently, a central composite design (CCD) coupled with response surface methodology (RSM) was performed in order to evaluate the effect of the experimental conditions "reaction time" and "enzyme concentration" on the enzymatic hydrolysis of macaw pulp oil. The enzyme concentration ranged from 0.5 to 1.8 mg/ml, according to preliminary tests. The range of reaction time was defined from 1 to 15 h, in which the center point was set at 8 h.
where X C and X t correspond to coded values for lipase concentration and reaction time, respectively. Eq. (2) had a correlation coefficient R 2 of 0.9393 and adjusted R 2 of 0.9166, which means that the statistical model accurately describes the correlation between the actual and predicted responses (Jun et al., 2019). The analysis of variance (ANOVA) for the model is presented in the Supplementary Material Section (Supplementary Table S5).
The regression coefficients of Eq.
(2) show that reaction time (X t ) was the most significant parameter on the hydrolytic degree, with a coefficient of 11.97, followed by the quadratic effect of enzyme concentration (X 2 C ), with a coefficient of 5.86, and, finally, the linear effect of enzyme concentration (X C ), with 5.60. The interaction between the parameters was not significant. Furthermore, the signals of the regression coefficient show that all linear parameters positively affect the response of the process. In contrast, the quadratic term of enzyme concentration negatively affects hydrolytic activity.
Comparing the experiments in which ReLip concentration was fixed and the time of reaction was changed (runs 1-2, 3-4, 7-8-center points), it is possible to observe the positive effect of reaction time on hydrolysis degree. Raising the variable from level −1 to level +1 promoted an increase of 1.35fold in the hydrolysis degree. In addition, comparing runs where the time of reaction was fixed and enzyme concentration was changed (runs 1-3, 2-4, 5-6-center points), it is possible to observe that an increase in enzyme concentration also promoted an increase in the hydrolysis degree. However, this effect was less pronounced. From level −1 to level +1, an increase of 1.18-fold in the hydrolysis degree occurred. The effect of the independent variables on hydrolysis degree can also be observed at the 3D response surface plot ( Figure 2B). It is possible to note that the region with the highest hydrolysis degree is obtained at the highest values of enzyme concentration and time of reaction, in accordance with the positive effect obtained in the statistic equation [Eq. (2)].
In order to maximize the hydrolysis degree, a canonical analysis for the complete regression model was applied. The canonical model is presented in Eq. (3). All characteristic roots (λ 1 and λ 2 ) presented negative signals, which means that a maximum stationary point was obtained. The optimal conditions suggested by the optimization were 3.0 mg/ml of enzyme (coded value 3.43) and 14 h 06 min of reaction (coded value 1.05), with a predicted hydrolysis degree of 84%. An experimental validation was conducted. The experiments at the optimal conditions were performed in duplicate. A hydrolysis degree of 81.47% was obtained, with a final FFA content of 92%, confirming the result predicted by the canonical analysis.
The results obtained in the present work, which were hydrolysis degree of 81.5% and hydrolyzate with final FFA of 92.0%, are very promising, since ReLip presented similar results on the hydrolysis of macaw pulp oil to those conducted with commercial and vegetable lipases ( Table 3). As previously mentioned, macaw palm is a very productive oleaginous tree, adapted to various types of soil and to semiarid ecosystems (Moura et al., 2009), which avoids conflict with areas used for food production. Additionally, its acidic oil cannot be used as food (Aguieiras et al., 2014), making it a great alternative as feedstock for biofuel production.

ReLip Is Highly Thermo-Tolerant Under Acidic Conditions and in the Presence of Organic Solvents
To evaluate the conformational stability of ReLip under harsh conditions such as low pH and the presence of organic solvents, the secondary structure profile was monitored by circular dichroism (CD) in a broad temperature range. CD analyses showed that ReLip presents a typical α-helical profile, with minimum molar ellipticities at 208 and 222 nm. Furthermore, at  different pH conditions (from 3 to 6.5), the secondary structure of the enzyme remained practically unchanged ( Figure 3A). The analysis of the melting temperature (T m ), by the measurement of thermal denaturation profiles at 220 nm ( Figure 3B), revealed that the T m varied from 73 • C to 84 • C for different pHs, with the maximum values obtained at pH 4.5 and 5.0, with a T m of 84 • C and 82 • C, respectively. These results corroborate the data obtained previously, in which higher activities were observed under similar conditions. Thermo-tolerant lipases such as Bacillus thermoamylovorans, Thermomyces lanuginosus, and Geobacillus thermoleovorans have been reported in the literature, with thermal stability up to approximately 60 • C (Sun et al., 2016;Mehta et al., 2017). ReLip shows T m values higher than other lipases considered thermostable, with an advantage of featuring high stability under acidic conditions, which is extremely desirable for many biotechnological applications, especially for biofuel production from crude or residual vegetable oils, which present high FFA contents. The conformational changes in the secondary structures were also evaluated with ReLip in the presence of organic solvents such as methanol and ethanol, compounds that are widely used in the biofuel industries. Surprisingly, ReLip had its secondary structure unaffected in the presence of up to 60% ethanol and methanol at room temperature (Supplementary  Figures S4A,B), which confirms its ability to tolerate organic solvents in high quantities, similar to that used in the industry. Figures 3C,D show the thermal denaturation profiles of ReLip in the presence of these organic solvents. Changes in the enzyme thermal stability were observed; T m values decreased as the amount of solvent was increased. In the presence of ethanol, from 30% to 60%, T m varied from 64 • C to 56 • C.
In the presence of methanol, it changed from 68 • C to 54 • C. Despite the decrease detected for both solvents, T m values are significantly high, showing the good thermal stability of the enzyme and its high potential for application in reactions that involve those solvents. The phenomenon of T m reduction in the presence of solvents may occur due to an alteration in the water-solvating layer (caused by organic solvents soluble in water) that surrounds the protein in aqueous solution, which may compromise the structural integrity of the enzyme. In addition, hydrophobic interactions of the structure can be affected by the presence of methanol and ethanol (Kamal et al., 2013;Sharma et al., 2016).

ReLip Has Unique Structural Properties That May Have Contributed to Its Acid and Methanol Tolerance
Despite a wealth of functional data on lipases, structural and mechanistic information are rather limited, especially for those from fungal origin. This is due to the difficulty of heterologous expression in a soluble and stable form, and the inherent recalcitrance to crystallization processes. The quaternary arrangement required for functional and structural stability is also partially understood with monomeric (Rhizomucor miehei) and dimeric (Thermomyces lanuginosus) lipolytic enzymes reported.
Consequently, the oligomeric arrangement of ReLip was also assessed by DLS (Supplementary Figure S5A) and SAXS (Supplementary Figures S6A,B). DLS analysis revealed that ReLip behaved as a monodisperse population (Pd < 15%) of monomers with an average hydrodynamic radius (Rh) of 3.05 nm. SAXS data showed that ReLip has a gyration radius of 2.4 nm, which is in agreement with DLS data. In addition, the obtained envelope from the SAXS curve showed that the low-resolution structure of ReLip consists of a monomeric globular form that corroborates with the three-dimensional structure obtained by X-ray crystallography (Supplementary Figure S6C).
Aiming to get structural information regarding the interesting functional properties of ReLip, the enzyme was submitted to crystallization tests. The obtained crystal (Supplementary Figure S5B) diffracted to 3.0 Å resolution, and the structure was solved using molecular replacement methods with the lipase from Penicillium cyclopium (PDB: 5CH8, 48.09% identity) as template. The ReLip structure (PDB: 6UNV) presents an α/β hydrolase fold and consists of a major eight-stranded mixed β-sheet, two minor two-stranded β-sheet arrangements and five α-helices (Figure 4). Unlike its orthologs, ReLip lacks the glycosylation site at position 33, which has a natural mutation with alanine instead of an asparagine. The literature reports that this site is important for the binding and recognition of micelles by TLL but apparently does not interfere in the catalytic activity (Peters et al., 2002). The lack of a glycosylation site in ReLip structure may have contributed to its expression in a folded and stable form in a bacterial organism.
As expected, the catalytic triad of ReLip consists of the residues Ser177, His290, and Asp232, which are fully conserved in other lipases and esterases ( Figure 4A). ReLip possesses three disulfide bonds (C67-C72, C53-C300, C135-C138) that are also conserved in TLL (Figure 4). One of them (C135-C138) is located close to the lid subdomain and is present in both lipases and esterases. The C53-C300 disulfide bridge is connecting the N-and C-termini, and is adjacent to the other C67-C72, conferring high structural stability to ReLip.
Usually, the lipolytic activity of some lipases increases significantly at the critical micelle concentration of the substrate. This phenomenon of interfacial activation is related to the presence of a hydrophobic patch of the lid subdomain that is protected from the solvent adopting a closed conformation (Skjold-Jørgensen et al., 2014;Khan et al., 2017). With the recognition of the substrate, the lid undergoes conformational changes to expose the hydrophobic patch to interact with the substrate. However, in some lipases that have minilids or do not even possess them, the interfacial activation does not occur, but these enzymes still present lipolytic activities. The three-dimensional structure of ReLip ( Figure 4C) confirmed the presence of a helical amphipathic lid, indicating that interfacial activation might be involved in ReLip mechanism of substrate binding such as its ortholog TLL. Indeed, in the crystallographic closed conformation of ReLip, the substrate could not be accommodated in the active-site cleft (Supplementary Figure S7), supporting that a conformational change is required for substrate binding. However, to better comprehend the molecular basis of the putative interfacial activation of ReLip, an in-depth investigation is needed.
Lipases are also known by the occurrence of cluster IxxWxxxxxF, which confers a high hydrophobic character to the lid subdomain (Skjold-Jørgensen et al., 2014). The literature shows that the presence of the tryptophan residue is essential for the interfacial activation and substrate binding. ReLip conserves the residues isoleucine and tryptophan, but the residue phenylalanine is replaced by an alanine. FAE also did not conserve the phenylalanine; instead, it has a tyrosine residue, which contributes to make its lid less hydrophobic (Hermoso et al., 2004). The lack of this aromatic residue in ReLip might be associated with the preference for shorter acyl chains compared to TLL, while preserving the ability to cope with water-insoluble substrates. , and esterase from Aspergillus niger (PDB code: 1USW). The diffraction experiment was performed at SSRL-SLAC beamline BL12-2 with a wavelength of 0.98 Å. Crystal diffracted and collected at 2.4-Å resolution. The phases were obtained by molecular replacement using the three-dimensional structure of MDLA lipase from Penicillium cyclopium (PDB: 5CH8). The refinement was performed to the resolution of 3.0 Å (due to crystallographic anisotropy) using Phenix. Another singular feature of ReLip is the presence of a salt bridge formed between Arg291 and Glu118 from the lid subdomain ( Figure 4B). This salt bridge is not conserved in other structurally characterized lipases with the lid subdomain such as TLL. In TLL, the glutamate residue is conserved in the lid, but the arginine is substituted by a leucine. The presence of this salt bridge in ReLip implicates in a higher protection of the hydrophobic pocket to bulky solvent, which can be correlated with the high tolerance of ReLip to organic solvents. On the other hand, this salt bridge can represent an additional energetic cost for ReLip to undergo the conformational changes required for substrate binding. The lid belonging to ReLip also has a rare property when compared to other lipases and even esterases, which is the presence of five amino acids with acidic nature, which may have contributed to its robust performance under acidic conditions, thus favoring the hydrolysis of macaw oil, a substrate with high acidity.
Although esterases and lipases have a fully conserved catalytic site, charge distribution on the catalytic interface is reported to be quite unlike (Hermoso et al., 2004). While esterases exhibit a more neutral surface charge distribution, lipases have a negatively charged surface. Surface charge distribution has been associated with the better performance of FAE at lower pHs. However, the charge distribution of ReLip largely differs from FAE, resembling that of TLL, despite having a better catalytic activity at lower pHs like FAE. This suggests that other factors aside from surface charge distribution at the catalytic interface are involved in the preference for acidic pHs for hydrolytic activity of these enzymes (Figure 5).

CONCLUSION
Several lipases have been studied for the hydrolysis of vegetable oils; nevertheless, only a few, especially those from castor bean seeds, were capable of hydrolyzing macaw oil, an extremely acidic substrate coming from a native palm that is cultivated in a wide variety of soil types. In addition, the enzyme showed other desirable properties for vegetable oil processing, such as tolerance to organic solvents and high temperatures. The crystallographic structure was elucidated indicating the presence of a helical amphipathic lid, which might be involved in a mechanism of interfacial activation with the exposure of the hydrophobic patch for substrate binding. The typical aromatic cluster in the lid IxxWxxxxxF is not fully conserved in ReLip with the natural mutation of the phenylalanine by an alanine. It decreases the hydrophobicity of the lid subdomain and is associated with the preference of short acyl chains, although preserving catalytic activity on insoluble substrates. The presence of a salt bridge between the lid and the vicinity of the active-site cleft confers a higher protection of the hydrophobic active site to bulky solvent, which might explain the tolerance of ReLip to organic solvents.

DATA AVAILABILITY STATEMENT
The atomic coordinates and structure factors of ReLip were deposited in the Protein Data Bank (http://ww.pdb.org/) under accession ID 6UNV. The original contributions presented in the study are included in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
LR, MS, NM, and CS carried out the expression and purification of ReLip and the oil hydrolysis reactions. RM and NM performed the functional experiments using synthetic substrates. NM performed the DLS analyses. MS carried out the CD experiments. PV, MM, and LZ designed and performed the crystallization trials, structure solution and refinement, and structural and SAXS analyses. LR, BK, HC, and AB contributed with the design of the experiments and with the analysis, treatment and discussion of the data obtained. LR wrote the article. LZ directed the overall study, analyzed the data, and wrote the manuscript. All authors read and approved the final manuscript.