- Research article
- Open Access
Optimization of physical conditions for the production of thermostable T1 lipase in Pichia guilliermondii strain SO using response surface methodology
BMC Biotechnologyvolume 17, Article number: 78 (2017)
Pichia guilliermondii was found capable of expressing the recombinant thermostable lipase without methanol under the control of methanol dependent alcohol oxidase 1 promoter (AOXp 1). In this study, statistical approaches were employed for the screening and optimisation of physical conditions for T1 lipase production in P. guilliermondii.
The screening of six physical conditions by Plackett-Burman Design has identified pH, inoculum size and incubation time as exerting significant effects on lipase production. These three conditions were further optimised using, Box-Behnken Design of Response Surface Methodology, which predicted an optimum medium comprising pH 6, 24 h incubation time and 2% inoculum size. T1 lipase activity of 2.0 U/mL was produced with a biomass of OD600 23.0.
The process of using RSM for optimisation yielded a 3-fold increase of T1 lipase over medium before optimisation. Therefore, this result has proven that T1 lipase can be produced at a higher yield in P. guilliermondii.
Lipases (triacylglycerol acylhydrolases) are natural catalysts, catalysing the hydrolysis of long chain triacylglyerides into their simpler di- and monoacylglycerol, as well as glycerol forms [1, 2]. Lipases are abundantly found in nature originated from plants, animals and microbes . Meanwhile, microbial thermostable lipases are known to be originated from bacterial, fungi and yeast with many applications in industrial and biotechnological processes [4,5,6].
Some lipase producing yeasts belonged to Pichia, Kluyveromyces, Candida and Torulaspora genre are methylotrophic in nature [7,8,9]. They metabolise methanol as their carbon and energy source [10,11,12]. Recombinant protein expression in methylotrophic yeasts are driven by the strong regulatory AOXp 1. Furthermore, activation of AOXp 1 for gene expression and protein production in some methylotrophs is independent of induction by methanol, even in the presence of AOXp 1, which is a mechanism referred to as de-repression and methanol induced expression [10, 13, 14] . Conventionally, biomass generation and gene expression phases are distinct to each other The carbon source used is either glycerol or glucose, which is mostly repressive to gene expression as it vary from each yeast and is a determining factor on the level of biomass generated . Yurimoto and co-authors illustrated that the level of AOXp 1 activities in different carbon sources is related to their gene expression . In some yeasts it is known as methanol induced expression while in some, it is a de-repressed activation of gene expression depending on the carbon source, level of de-repression and choice of host [12, 16].
However, a new yeast P. guilliermondii isolated from spoilt orange was used as a recombinant host for the expression of thermostable T1 lipase gene isolated from Geobacillus zalihae [17, 18]. Similar to methylotrophic yeasts, gene expression and protein production in P. guilliermondii have been found to be driven by methanol induced AOXp 1 . In this study, P. guilliermondii was found capable of generating biomass using glycerol as a carbon source and relatively undergoing gene expression without any carbon source as observed in the study on Hansenula polymorpha  depending on the level of de-repression.
Previously, Escherichia coli has been used as a recombinant host for T1 lipase production, but was characterised by low expression . It is important to note that T1 lipase gene has been also cloned and expressed in P. pastoris but recently, it was observed that P. guilliermondii was observed as a better producer of T1 lipase under the control of AOXp 1 . Therefore, it is appropriate to screen and optimise the enzyme production, initially by media optimisation, followed by optimisation of physical conditions significant for the enhancement of T1 lipase production in P. guilliermondii.
Varying one parameter at a time (Conventional method) in a bid to optimise production process is labour intensive, time consuming and ignores the co-interaction between parameters, leading to errors in identifying the optimum production conditions [20, 21]. As such, recent studies have shown the use of statistical approaches, which eliminates the limitation associated with the use of the conventional method of optimisation. Plackett-Burman Design (PBD) has been used in several biotechnology applications for screening multiple parameters and identifying the best parameters with significant effects . Box-Behnken Design (BBD) optimisation method has been also employed for effective optimisation of PBD experiments . However, there is no known study made on the optimisation of physical conditions during the production of recombinant T1 lipase. This gap in knowledge could affect the accurate and effective production of T1 lipase in P. guilliermondii. Thus, this study aimed at optimising buffered and non-buffered media as well as physical conditions for the production of recombinant T1 lipase using PBD and BBD Response Surface Methodology.
Materials and methods
Strain and inoculum preparation
The recombinant yeast P. guilliermondii from previous study  was maintained on YPD- Zeocin agar plate (1% yeast, 2% peptone, 2% Bacteriological agar, 20% dextrose and 100 μg/mg Zeocin). A single colony was inoculated into 10 mL YPD broth at 30 °C with an agitation of 250 rpm overnight.
Two different media were used for T1 lipase production, which are buffered and non-buffered media (BMY and YPT, respectively). BMGY medium (Buffered complex glycerol medium - 1% yeast extract, 2% peptone, 1% 100 mM potassium phosphate pH 6.0, 1.34% yeast nitrogen broth, 0.4 mg/ L biotin and 1% glycerol), and YPTG (Yeast extract, Peptone, Tryptic Soy Broth and Glycerol) were used for biomass generation. Meanwhile, BMY (Buffered complex medium) with composition similar to BMGY without glycerol and YPT with composition similar to YPTG without glycerol were used as lipase production media.
Determination of lipase activity.
Lipase activity was assayed colorimetrically  with slight modifications . 2.5 mL olive oil, 50 mM phosphate buffer (pH 9) emulsion (1:1, v/v) and 20 μl of 0.02 M CaCl2 were added to 1 mL of culture filtrate. The reaction mixture was incubated for 30 min at 70 °C. Enzyme reaction emulsion was terminated by adding 1 mL of 6 N HCl and 5 mL of isooctane. The reaction mixture was vortexed for 30 s and two layers were formed. 4 mL of the upper isooctane layer containing the fatty acid was transferred into a new test tube for analysis after that 1 mL copper (II) pyridine reagent (pH 6.1) was added and vortexed for 30 s. The absorbance of the upper layer was read at 715 nm. One unit of lipase activity was defined as the rate of free fatty acid formed in μmole per minute.
Plackett-Burman (PB) screening design
Factorial designs were evaluated at levels −1 and +1 for low and upper level respectively. The parameters screened are inoculum size, incubation time, temperature, pH, culture volume and agitation speed. Twelve sets of trial experiments were generated from Design Expert (version 7.0.) software package. PBD was patterned in line with the first-order model (Y = Σβ0 + ΣβpXp), where, Y is the response (T1 lipase), β0 is the intercept, β is the coefficient of the variable and Xp is the independent variable. Screened parameters were represented on Pareto chart of standardised effects.
Box-Behnken design (BBD)
Physical parameters were evaluated at three different interactive levels, which are low (−), middle (0) and high levels (+). Statistical design was used for optimisation involving 17 sets of trial experiment, to determine the co-efficient of the model, as well as to predict the response and its experimental validation. The response was fitted into a second order polynomial equation (Eq. 1) and optimum levels were represented in response surface plots.
where, Y is the response, β0 the intercept, β the coefficient of the variable, Xp, Xq, Kz are equal to variable (parameter) 1, 2 and 3, Xpqz represent the linear effect, X(p2q2z2) are the quadratic form, X(pq,pz,qz) are the interactive effects and ϵ is the residual relative to the experiment .
Responses observed, were analysed using two way analysis of variance (ANOVA). This analysis was generated having the definition for independent variables experimented. The linear, quadratic and interaction regression coefficient of each term in the model were determined. Using the F-value at a probability (P) of 0.05 and confidence level of 95%, the significance of all terms in the polynomial was statistically analysed and all coefficients were computed using Expert Software version 7.0.
Results and discussion
Effects of different media on T1 lipase production
Alterations in the pH of a cultivation medium can lead to a change in the charge of a medium, which affects cell growth and recombinant protein production. pH is mentioned as a crucial fermentation factor affecting cell growth and recombinant proteins production if not controlled [26,27,28]. In this study, recombinant T1 lipase was produced in a pH controlled medium (BMY) and non-pH controlled medium (YPT). Figure 1 illustrates that T1 lipase production was higher in BMY with T1 lipase yield of about 0.7 U/mL compared to YPT with a yield of 0.5 U/mL. Both media gave optimum production after 28 h of cultivation and declined afterwards. On the other hand, both media provided vitamins and trace elements, as growth enhancers for biomass and energy generation during enzyme production. Furthermore, BMY provided a constant pH condition for T1 lipase production.
Six parameters were screened, and their effects on T1 lipase production were observed (Table 1). These parameters are pH, incubation time (IT), inoculum size (IS), culture volume (CV), temperature (T) and agitation speed (AS). Using Pareto chart (Fig. 2), significant and non-significant physical parameters were placed above and below the horizontal line (t-value limit) respectively, in descending order. The vertical line shown in the chart indicted the statistical significance (P > 0.05), at 95% confidence level. The results obtained for each parameter showed that pH > incubation time > inoculum size, exerted significant effects on T1 lipase production.
Optimisation of screened parameters using BBD
Parameters identified by PBD experiment (pH, incubation time and inoculum size) having significant effects on T1 lipase production, were further examined and optimised through BBD methodology. This method allows the interaction of pH, incubation time and inoculum size in a 17 trial randomised experiments (Table 2) with their optimum levels observed as maximum amount of lipase was produced.
Statistical analysis of variance (Table 3) was used to investigate the effectiveness of BBD model. The F-value (Fisher’s statistical analysis) and the p-value (<0.0001) were used as tools for evaluating the significance of the model. The model terms of Prob > F (less than 0.05) were deemed significant while, those with prob. > F greater than 0.10 were considered insignificant [29, 30]. The analysis of BBD model had a prob. > F of 0.0004, which was considered significant to T1 lipase production. The ‘Lack of Fit’ (Table 3), which was insignificant to the model, was considered good.
Generally, the coefficient for determination for model variable (R2), and its prediction of the response is significant when, its value is closer to 1, which is 100% [31, 32]. Thus, the model was considered significant with over 98% of model variability (R2–0.9835) explained. The adj R2 with value of 0.9605 was in reasonable agreement with the predictable R2–0.8814. The signal to noise ratio was measured by a term called adequate precision. Any ratio greater than 4 is acceptable [32, 33]. In this model, the ratio was19.512, which is greater than 4 thus, it was considered acceptable.
To further explain the dependence of T1 lipase production on physical parameters, a second order polynomial equation (Eq. 2) was generated from the statistical analysis of BBD model.
Where A = pH, B = induction time and C = inoculum size.
Diagnostic case tools
Normality percentage plot measures standard deviations, adequate prediction and interpretation of the model. Normality percentage plot was used as a tool to check the efficiency of the model and to infer if the internally studentised residuals follows a normal distribution where data points fall along the straight line of regression. This is necessary prior response surface optimisation of the model to avoid misleading results from the model [34, 35]. Values observed (Fig. 3a) aligned close to the ideal line of regression indicated the accuracy and significance of the model, which was in line with the normal distribution assumption. Besides, perturbation plot (Fig. 3b), was used as a comparative tool to determine the sensitivity of individual parameter to the response (T1 lipase). A deviation from the reference point indicated the magnitude of effect exerted by a parameter on the response, which confirmed the efficiency of the model. Parameter ‘A’ (pH) was observed to have the most deviation (curvature) from the reference point. This indicated the sensitivity of the response to pH and the significance of pH to T1 lipase (response) production. Deviation of parameter ‘C’ (inoculum size) from the reference point was reported though the deviation was less intense compared to Parameter ‘A’. The deviation of parameter ‘B’ (induction time) to the reference point was observed to be insignificant to T1 lipase production. Studentised residuals (Fig. 3c and d) are often used to explain the systematic deviation from the hypothesis that errors examined are often independent of each other and evenly distributed [36, 37]. It also measures the values of standard deviation that distinguishes experimental values (actual values) and values generated from BB model design (predicted) experiments . In this study, residual errors were within the allowable limits having an observed random error distribution without confirming any pattern but maintaining a constant variation. The adequacies of these diagnostic plots indicated the applicability of the models produced by RSM for optimally producing T1 lipase in P. guilliermondii.
3D response surface and contour plots
The interaction of two parameters and their effects on lipase production were evaluated while, maintaining the third parameter at a constant level. 3D responses aided in the visual determination of optimum levels of each parameter as they interact. Basically, the 3D responses have either elliptical representation, which is considered significance or the circular representation which is considered to be negligible [39, 40]. As the level of these parameters changes, the interactive effects they exert on the response (recombinant T1 lipase) are also changed. Interaction between pH and inoculum size was studied while, maintaining incubation time at a constant value of 36 h (Fig. 4a). At a very high pH of 10 and inoculum size of 2% lipase production was optimum with a yield of 1.9 U/mL. Decrease in lipase production became trendy as a decrease in pH was observed with subsequent increase in inoculum size. In Fig. 4b, lipase production was found to be optimum in 24 h incubation time, whereas 2% of the inoculum size, pH was maintained at a constant range of eight (8). The level of lipase produced at these optimum parameters levels is 1.7 U/mL.
Validation of the model
The accuracy of the model was determined (Fig. 5). A point prediction of the model was generated by the statistical analysis software. pH (6), incubation time (24 h), inoculum size (2%) lipase yield of 2.1 U/mL and biomass of OD600 24.6 were predicted. These predicted parameters were experimented, demonstrating the experimental yield for lipase production of 2.0 U/ mL and a biomass of OD600 23.0. This observed result showed a high degree of correlation with the predicted one. Hence, the model was considered accurate and reliable.
This study is the first attempt of optimising T1 lipase yield in P. guilliermondii under the regulatory control of AOXp 1, without methanol. Statistical methods of PBD and BBD were used to screen and optimise physical conditions (parameters). This methods allowed the elimination of physical conditions that were insignificant to T1 lipase production while, providing a platform to study the individual and combined effects of a reduced number of physical conditions essential to lipase production. It was observed that pH, incubation time and inoculum size were physical conditions significance to T1 lipase production in P. guilliermondii. These conditions were optimised at various randomised fermentation levels while the predicted conditions for lipase yield were experimentally validated. The observed relationship between the predicted and experimented values as they aligned in straight line to the line of regression has indicated the suitability of the model in predicting conditions and optimised levels at which lipase production was optimum. A 3-fold increase in lipase yield after optimisation was observed, with 300% increase in T1 lipase production. RSM has been proven to be effective in enhancing T1 lipase production in P. guilliermondii without induction by methanol, thus reducing the cost of production and hazards associated with methanol usage.
Buffered complex – glycerol medium
Buffered complex - medium
- CaCl2 :
Coefficient of variation
Degree of freedom
- OD600 :
Optical Density at 600 nm
Plackett- Burman Design
- R2 :
Rotation per minute
Sum of Square
- Std Dev:
Volume per volume
Weight per volume
Yeast extract Peptone Dextrose
Yeast extract Peptone Tryptic soy
Yeast extract Peptone Tryptic soy Glycerol
Li H, Zhang X. Characterization of thermostable lipase from thermophilic Geobacillus sp. TW1. Protein Expr Purif. 2005;42(1):153–9.
Balan A, et al. Purification and characterization of a thermostable lipase from Geobacillus thermodenitrificans IBRL-nra. Enzyme research. 2012:2012.
Abdel-Fattah YR, Gaballa AA. Identification and over-expression of a thermostable lipase from Geobacillus thermoleovorans Toshki in Escherichia Coli. Microbiol Res. 2008;163(1):13–20.
Salihu A, Alam MZ. Thermostable lipases: an overview of production, purification and characterization. Biosciences Biotechnology Research Asia. 2014;11(3):1095–107.
Masomian M, et al. A new thermostable and organic solvent-tolerant lipase from Aneurinibacillus thermoaerophilus strain HZ. Process Biochem. 2013;48(1):169–75.
Gandhi NN. Applications of lipase. J Am Oil Chem Soc. 1997;74(6):621–34.
Singh AK, Mukhopadhyay M. Overview of fungal lipase: a review. Appl Biochem Biotechnol. 2012;166(2):486–520.
Romo-Sánchez S, et al. Yeast biodiversity from oleic ecosystems: study of their biotechnological properties. Food Microbiol. 2010;27(4):487–92.
Gupta R, Gupta N, Rathi P. Bacterial lipases: an overview of production, purification and biochemical properties. Appl Microbiol Biotechnol. 2004;64(6):763–81.
Yurimoto H, Oku M, Sakai Y. Yeast methylotrophy: metabolism, gene regulation and peroxisome homeostasis. International journal of microbiology. 2011;2011
Ahmad M, et al. Protein expression in Pichia Pastoris: recent achievements and perspectives for heterologous protein production. Appl Microbiol Biotechnol. 2014;98(12):5301–17.
Hartner FS, Glieder A. Regulation of methanol utilisation pathway genes in yeasts. Microb Cell Factories. 2006;5(1):39.
Yurimoto H, Sakai Y. Methanol-inducible gene expression and heterologous protein production in the methylotrophic yeast Candida Boidinii. Biotechnol Appl Biochem. 2009;53(2):85–92.
Yurimoto H. Molecular basis of methanol-inducible gene expression and its application in the methylotrophic yeast Candida Boidinii. Biosci Biotechnol Biochem. 2009;73(4):793–800.
Cox H, et al. Constitutive expression of recombinant proteins in the methylotrophic yeast Hansenula Polymorpha using the PMA1 promoter. Yeast. 2000;16(13):1191–203.
Gellissen G, et al. New yeast expression platforms based on methylotrophic Hansenula Polymorpha and Pichia Pastoris and on dimorphic Arxula adeninivorans and Yarrowia lipolytica–a comparison. FEMS Yeast Res. 2005;5(11):1079–96.
Oslan SN, et al. Locally isolated yeasts from Malaysia: identification, phylogenetic study and characterization. Acta Biochim Pol. 2012;59(2):225–9.
Oslan SN, et al. A newly isolated yeast as an expression host for recombinant lipase. Cell Mol Biol Lett. 2015;20(2):279–93.
Oslan SN, et al. Pichia Pastoris as a host to overexpress the thermostable T1 lipase from Geobacillus zalihae. GSTF International Journal on Bioformatics & Biotechnology (JBio). 2014;3(1):7.
Morowvat MH, et al. Optimization of fermentation conditions for recombinant human interferon beta production by Escherichia Coli using the response surface methodology. Jundishapur journal of microbiology. 2015;8(4)
Weuster-Botz D. Experimental design for fermentation media development: statistical design or global random search? J Biosci Bioeng. 2000;90(5):473–83.
Plackett RL, Burman JP. The design of optimum multifactorial experiments. Biometrika. 1946;33(4):305–25.
Box, G.E., W.G. Hunter, and J.S. Hunter, Statistics for experimenters: an introduction to design, data analysis, and model building. Vol. 1. 1978: JSTOR.
Kwon DY, Rhee JS. A simple and rapid colorimetric method for determination of free fatty acids for lipase assay. J Am Oil Chem Soc. 1986;63(1):89–92.
Leow TC, et al. High level expression of thermostable lipase from Geobacillus sp. strain T1. Biosci Biotechnol Biochem. 2004;68(1):96–103.
Lan D, et al. Enhancing production of lipase MAS1 from marine Streptomyces sp. strain in Pichia Pastoris by chaperones co-expression. Electron J Biotechnol. 2016;22:62–7.
de Almeida, A.F., S.M. Tauk-Tornisielo, and E.C. Carmona, Acid lipase from Candida viswanathii: production, biochemical properties, and potential application. Biomed Res Int, 2013. 2013.
Zouaoui B, Bouziane A. Isolation, optimisation and purification of lipase production by Pseudomonas Aeruginosa. J Biotechnol Biomater. 2011;1(7):1–4.
Li X, Liu Z, Chi Z. Production of phytase by a marine yeast Kodamaea ohmeri BG3 in an oats medium: optimization by response surface methodology. Bioresour Technol. 2008;99(14):6386–90.
Salihu A, et al. Optimization of lipase production by Candida Cylindracea in palm oil mill effluent based medium using statistical experimental design. J Mol Catal B Enzym. 2011;69(1):66–73.
Haaland, P.D., Experimental design in biotechnology. Vol. 105. 1989: CRC press.
Kaushik R, et al. Statistical optimization of medium components and growth conditions by response surface methodology to enhance lipase production by Aspergillus carneus. J Mol Catal B Enzym. 2006;40(3):121–6.
Chauhan B, Gupta R. Application of statistical experimental design for optimization of alkaline protease production from bacillus sp. RGR-14. Process Biochem. 2004;39(12):2115–22.
Kar S, Ray RC, Mohapatra U. amylase production by Streptomyces Erumpens MTCC 7317 in solid state fermentation using response surface methodology (RSM). Pol J Microbiol. 2008;57(4):289–96.
Khusro A, et al. Optimization of thermo-alkali stable amylase production and biomass yield from bacillus sp. under submerged cultivation. Fermentation. 2017;3(1):7.
Manderson D, Dempster R, Chisti Y. Production of an active recombinant Aspin antigen in Escherichia Coli for identifying animals resistant to nematode infection. Enzym Microb Technol. 2006;38(5):591–8.
Mohajeri A, et al. Expression and secretion of Endostar protein by Escherichia Coli: optimization of culture conditions using the response surface methodology. Mol Biotechnol. 2016;58(10):634–47.
Akbari V, et al. Optimization of a single-chain antibody fragment overexpression in Escherichia Coli using response surface methodology. Research in pharmaceutical sciences. 2015;10(1):75.
Fu J, Zhao Y, Wu Q. Optimising photoelectrocatalytic oxidation of fulvic acid using response surface methodology. J Hazard Mater. 2007;144(1):499–505.
Qiu P, et al. Application of box-Behnken design with response surface methodology for modeling and optimizing ultrasonic oxidation of arsenite with H2O2. Cent Eur J Chem. 2014;12(2):164–72.
This research work was funded by the Fundamental Research Grant Scheme Projects, (Code: 02-02-14-1499). I wanted to thank the Ministry of Higher Education Malaysia, and Tertiary Education Trust Fund (Tetfund) Nigeria for granting the academic scholarship for A. M. L.
This research work was funded by the Fundamental Research Grant Scheme Projects, (Code: 02–02–14-1499).
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that there is no conflict of interests regarding the publication of this paper.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.