Materials used in this study
All chemicals used in this study were of analytical grade and purchased from Sigma Aldrich (USA) and Fischer scientific unless stated otherwise. All restriction enzymes, DNA ligase and Q5 Taq DNA polymerase used for the PCR and cloning were purchased from New England Biolabs (NEB) (USA).
Microorganism, culture maintenance and inoculum preparation
The current study made use of strain originated from adaptive evolution of engineered Y. lipolytica PSA02004 with deletion of Ylsdh5 gene encoding a sub-unit of succinate dehydrogenase [26]. The recombinant Y. lipolytica strain was preserved in 20 % glycerol (v/v) at -80 °C and maintained on a petri dish containing YPD agar medium (1% yeast extract, 2% peptone, 2% dextrose and 2% agar) at pH 7.0 and 30 °C. The seed culture was grown in a 250 mL Erlenmeyer flask containing 50 mL minimal medium (see Section Submerged cultivations in shake flask). The flasks for seed culture were inoculated by transferring a loopful of 48 h culture grown on YPD plate. The final pH of the medium prior to sterilization was adjusted to 6.8. Cultivation was carried out for 24 h at 30 °C on a rotary shaker at an agitation speed of 250 rpm.
Cloning and expression of heterologous xylose assimilation gene in Y. lipolytica strain
Escherichia coli (DH5α) strain used for cloning and plasmid propagation. The strain was cultivated in the Lysogeny broth (LB) liquid medium at 37 °C. The gene encoding xylose reductase (XR) (YALI0D07634), xylitol dehydrogenase (XDH) (YALI0E12463) and xylulokinase (XK) (YALIF10923) were extracted from the genome of Po1d using the appropriate primers (Table S2), the golden gateway (GG) assembly was constructed according to former studies [51,52]. The GG was constructed with Scaffold of three genes comprising three transcription units and selection marker, flanked with integration targeting sequences, constructed on a destination vector backbone. Each gene was flanked with 396 nt of TEF promoter and 122 nt of Lip2 terminator sequences, both native to Y. lipolytica. URA3 (1289 nt) gene was used as selection marker in this assembly. Random integrations in Y. lipolytica PSA02004 were driven through zeta sequences (305 nt and 395 nt for UP and DOWN respectively). The expression vector is linearized using NotI enzyme and gel purified before transformation in Y. lipolytica. The overexpression cassette was transformed in the genome of Y. lipolytica using lithium acetate method described by Le Dall et al. [54]. The transformants were selected on YNBUra plates, the genomic DNA was isolated using the protocol developed by Lõoke et al. [55], and the positive transformants were identified with PCR. All the plasmids and strains used in this study are listed in Table S1 (supplementary information).
Submerged cultivations in shake flask
The minimal medium used for fermentation had the following composition: xylose, 20 g/L; yeast nitrogen base (YNB), 1.7 g/L; NH4Cl, 1.5 g/L. The medium was prepared in 50 mM phosphate buffer. In case of co-fermentation with two carbon sources, each one was used at a level of 20 g/L. The initial pH was adjusted to 6.8 before inoculation by using 5 N NaOH. The submerged cultivations were carried out in 500 mL shake flasks containing 100 mL working volume. The flasks were inoculated with fresh inoculum at OD600 of 0.1 and kept at 30 0C under constant shaking at 250 rpm on a rotary shaker. The xylose-rich lignocellulosic hydrolysate (with an initial xylose concentration of 42.8 g/L) from sugarcane bagasse was obtained from Nova Pangea Technologies, UK.
Measurement of xylose reductase (XR) and xylitol dehydrogenase (XDH) activities
For measuring the enzymatic activities, cell free extract was prepared. Initially, the cells were harvested by centrifugation at 8,000 x g and 4°C for 10 minutes. The cell pellet was then washed twice with 50 mM phosphate buffer (pH 7.2) and resuspended in the buffer. The cell disruption was performed in homogenizer by mixing the above-mentioned cells with 0.5 g (0.3 mm) glass beads and vortexed for 10 minutes. The homogenized mixture was centrifuged at 8,000 x g and 4°C for 10 minutes, the supernatant was collected and used for quantifying enzyme activities. The protein concentration was determined by the Bradford method [55].
The activities of xylose reductase (XR) and xylitol dehydrogenase (XDH) were measured using a UV spectrophotometer (Jenway 6310, UK). The molar extinction coefficient of NADPH and NAD+ used for calculation is 6,220 m-1cm-1. The XR activity was measured by the reduction of the coenzyme NADPH at 30 °C in a reaction medium consisting of 0.17 mM NADPH, 0.17 M xylose, 0.25 mg cell extract, and the final volume was made up to 0.5 mL using 0.1 M phosphate buffer. One unit of XR enzyme activity was defined as the amount of enzyme that catalyzed the oxidation of 1 μmol of NADPH per minute at 30°C. The quantification of XDH activity was based on reduction of the coenzyme NAD+ at 30°C. For XDH measurement, the reaction mixture consists of 1.5 mM NAD+, 0.15 M xylitol, 0.25 mg cell extract with a total volume of 0.5 mL made up by 0.1 M Tris buffer. One unit of XDH was defined as the amount of enzyme catalyzing the oxidation of 1 μmol of NAD+ per minute at 30°C [56].
Substrate inhibition kinetic studies for recombinant Y. lipolytica
The substrate inhibition studies were performed in shake flask using a minimal medium as mentioned in Section 2.4 at different concentration of xylose (20, 40, 60, 80, 100, 120 g/L). Specific growth rate (µ) were calculated using Equation (1):
Where X= DCW (g/L) and t= fermentation time.
Elevated levels of substrates may hamper the microbial growth and productivity of Y. lipolytica as the multiple substrate can bind to the same site. Several inhibition models derived from the Monod’s kinetics are used in this study to predict the kinetic parameters from the experimental data, and these are presented in Table 1. The parameters of different models were estimated from the experimental results using MATLAB (Mathwork R, 7.1). Since the models had non-linear coefficients, the parameters were quantified iteratively with the aid of non-linear least square algorithm.
Comparison of models for acceptability
High regression coefficient or low Mean Square error value for any kinetic model resulted in the best fitted model, which is achieved at a cost of complexity in the model. Since models with different complexity (i.e. degree of freedom) were chosen, it was a key to test the model which was consistent with the experimental data[57]. In this study best fit 3 and 4 parameter models were compared using Akaike information criterion (AIC).
The Akaike information criterion (AIC) is defined by the following equation:
where “b” = prm + 1. When there are few data points, the corrected AIC (AICc) is used.
The model with lower AICc value is more likely to be correct and the probability (pAIC) that the more complex model is correct is given by
Bioreactor studies
The batch experiments were performed in a 2.5 L bench-top bioreactor (Electrolab Bioreactors, UK) with 1.0 L working volume. The minimal medium with 60 g/L xylose was used for running bioreactor experiments. In case of lignocellulosic hydrolysate, the xylose concentration was 40 g/L. The temperature, agitation speed and aeration rate were controlled at 30 ⁰C, 600 rpm and 2.0 L/min, respectively. The starting pH was 6.8, and it remained uncontrolled during the fermentation. For fed-batch fermentations, the residual xylose concentration was maintained at or above 10 g/L with concentrated feed containing 500 g/L xylose and 5 g/L yeast extract.
Analytical methods
The samples were withdrawn periodically and analyzed for OD600, pH, residual glucose, xylose, succinic and acetic acid. Cell growth was quantified by measuring the optical density at 600 nm wavelength in a 1 mm-path-length cuvette using a double beam spectrophotometer (Jenway 6310, UK). One unit of absorbance at 600 nm corresponded to a cell dry weight (CDW) of 0.21 g/L. The concentrations of glucose, xylose, succinic and acetic acid were measured by high performance liquid chromatography (Agilent Technologies 1200 series, USA). The supernatants obtained by centrifugation of the culture samples at 10,000 × g for 10 min, were filtered through a 0.22 µm PVDFmembrane (Sartorious, Germany)) and eluted using Rezex ROA-Organic Acid H+ (Phenomenex, USA) column at 60 °C attached with refractive index detector (RID) and Diode Array Detector (DAD). The mobile phase and flow rate were 0.5 mM H2SO4 and 0.4 mL/min, respectively. All measurements were conducted in triplicates and the values were averaged. The standard deviation was no more than 10 %.