Sugar-Improved Enzymatic Synthesis of Biodiesel with Yarrowia

May 2, 2017 - In this research, it was found that in Yarrowia lipolytica Lipase 2 (YLLIP2)-catalyzed biodiesel production, the content of fatty acid m...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of South Dakota

Article

Sugar-Improved Enzymatic Synthesis of Biodiesel with Yarrowia lipolytica Lipase 2 Hao Cao, Meng Wang, Li Deng, Luo Liu, Ulrich Schwaneberg, Tianwei Tan, Fang Wang, and Kaili Nie Energy Fuels, Just Accepted Manuscript • Publication Date (Web): 02 May 2017 Downloaded from http://pubs.acs.org on May 8, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Energy & Fuels is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Sugar-Improved Enzymatic Synthesis of Biodiesel with Yarrowia lipolytica Lipase 2 Hao Cao,1,2 Meng Wang, *,1 Li Deng,1 Luo Liu,1 Ulrich Schwaneberg, 2,3

Tianwei Tan,1 Fang Wang,1 Kaili Nie*,1

1 Beijing Bioprocess Key Laboratory, College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, PR China

2 Lehrstuhl für Biotechnologie, RWTH Aachen University, Worringerweg 3, 52074 Aachen, Germany

3 DWI-Leibniz Institut für Interaktive Materialien, Forckenbeckstraße 50, 52056 Aachen, Germany

1

ACS Paragon Plus Environment

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 32

KEYWORDS Yarrowia lipolytica Lipase 2; Biodiesel; Glucose; Sugars; Methanol resistance; Molecular dynamics.

ABSTRACT In this research, it was found that in Yarrowia lipolytica Lipase 2 (YLLIP2) catalyzed biodiesel production, the content of fatty acid methyl esters (FAMEs) was increased by about 10% with the addition of

D(+)-glucose.

The result indicated that the

D(+)-glucose

could be used as an

effective additive in YLLIP2 catalyzed biodiesel production. According to the results above, the single factor experiments of key parameters in the process were first carried out. Based on the single factor experiment results, a five-factor, three-level response surface method was adopted to obtain the optimal reaction conditions: lipase dosage 40 IU/g oil, D(+)-glucose to lipase 1:1.05 (w/w), water content 1.95%, and reaction temperature 39.4 °C. A stoichiometric amount of methanol (11.5%, methanol/waste oil) was added in 6 times (every 4 hours for each addition). A subsequent pilot scale production in a 5-ton reactor was carried out to check the performance of this method, and a good biodiesel content of 91.4% was obtained. Finally, molecular dynamics (MD) simulation was adopted to help to explain the possible functions of D(+)-glucose acting on the lipase in the process. The simulation results indicated that one of the functions maybe attributed to

D(+)-glucose

preventing methanol to diffuse into YLLIP2, thus resulting in the

prevention of the methanol disruption to protein.

ACS Paragon Plus Environment

2

Page 3 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Introduction Lipase (glycerol ester hydrolase EC 3.1.1.3) has a broad range of applications in industry through hydrolysis, esterification and transesterification reactions1. In the last decade, chemical methods were frequently used for biodiesel production, but these methods have drawbacks, such as the need for removal of catalysts, high energy consumption and negative environmental issues. For those reasons, lipase has been widely investigated in biodiesel production to avoid the drawbacks of chemical methods2. However, the sensitivity of lipase to polar solvents, such as methanol, limits the widely utilization of lipases in the industrial applications. Methanol is a commonly used alcohol in biodiesel synthesis due to its low cost3 and availability4. Moreover, the use of methanol is advantageous as methanol could be more reactive than other short chain alcohols for biodiesel production5 and be easily and simultaneously separated with glycerol in the process4. In principle, the theoretically optimal methanol: triacylglycerols molar ratio is 3:1 in biodiesel production. An increase in methanol concentration was suggested for to a higher content of biodiesel. However, when the concentration of methanol is excess to a certain amount or even at the stoichiometric condition, methanol still exhibits the inactivation, inhibition, or denaturation effect on lipase. The explanations of methanol caused lipase catalytic activity decrease were previously reported in the view of disruption of enzyme-water interaction6, contribution of substrate-solvent interaction7, disruption of intra-protein interactions8 and competitive inhibition9. Consequently, several strategies were employed to overcome the deleterious effect of the presence of methanol on lipase (Table 1). Host screening strategy indicated that lipases could be

ACS Paragon Plus Environment

3

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 32

classified as sensitivity or robustness to methanol. Lipase from Burkholderia and Stenotrophomonas genus appeared to be highly tolerant to methanol, while Candida genus lipase and several others were methanol sensitive ones. In the case of protein engineering, modifications of surface amino acid residues of lipases could change the surface hydrophobicity in order to influence the contact between lipase and solvent, resulting in improving the lipase stability in methanol solvent. Towards process engineering, a stepwise addition strategy requires a critical selection of the most suitable bioreactors when applied to the industrial scale; both cosolvent and additive strategies have the problem of raw materials’ cost and the requirement for separation from the products. Table 1. Improved strategies for to overcome the negative effect of the presence of methanol on lipase Strategy

Microorganism engineering

Protein engineering

Process engineering

Strain

Method

Result

Ref.

Burkholderia glumae

Host screening

Highest activity in the presence of 50–70 % methanol

10

Stenotrophomonas maltophilia

Host screening

11

Proteus mirabilis

Directed evolution

Geobacillus stearothermophilus

Semi-directed evolution

Candida antarctica

Rational design

Thermomyces lanuginosus Pseudomonas fluorescens Yarrowia lipolytica

Stepwise addition Cosolvent

82.2 % residual activity in 20% methanol solvent for 24 h 50 fold longer halfinactivation time in 50% methanol 23 fold longer halfinactivation time in 70% methanol 75 % residual activity in 80% methanol solvent for 24 h Optimal yield reached 97 % from 85% Optimal yield reached 87.6 % from 61.6% Optimal yield reached 89 % from 68 %

Additive

12

13

14

15

16

17

ACS Paragon Plus Environment

4

Page 5 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Candida antarctica

Replacement of Acyl acceptor

Optimal yield reached 92 % with stoichiometric methyl acetate

18

In our previous work, β-cyclodextrin (β-CD) was used as additive to improve the biodiesel synthesis with Yarrowia lipolytica Lipase 2 (YLLIP2)

17

. The results of molecular dynamics

(MD) simulation indicated that the mechanism of the action of this circular saccharide on the protein was the interaction between the hydroxyl groups of β-CD and the polar residues on the surface of lipase, thus enhancing the methanol resistance19. β-CD is constructed by 7 glucose units. It has been known that sugars may increase the thermal stability of proteins, and can be used as protective agent in protein lyophilization 20, 21. However, literature scarcely reported that sugar could protect lipase against being denatured by polar solvents through direct addition into the reaction system. The aim of this research was to investigate whether the glucose could be an additive to improve the efficiency of the biodiesel synthesis process catalyzed by YLLIP2, and to find out the optimal reaction conditions. Furthermore, MD simulations were utilized to understand and explain the positive effect of glucose on the YLLIP2 in a methanol solvent in order to develop a rational strategy to improve the lipase-catalyzed reaction using additives. Materials and methods Experimental Materials Waste cooking oil (WCOs) was obtained from Lvming Co. Ltd (Shanghai, China). It contains 83.9% of free fatty acids (FFAs), 0.5% of monoacylglycerols (MAGs), 6.9% of diacylglycerols (DAGs), and 8.7% of triacylglycerols (TAGs). Since Shanghai is a coastal city in South China, the temperature and air humidity are relatively high year-round. For this reason, the waste cooking oil generally had a serious rancidification during collection, storage, and transportation processes. Most of the WCOs used by Lvming had a free fatty acid content around

ACS Paragon Plus Environment

5

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 32

70%-80% by weight. The average fatty acid molecular mass was 275.9 g/mol according to the fatty acid compositions (Table 2). According to the composition of WCOs and average fatty acid molecular mass, the stoichiometric amount of methanol was 1.15g with 10 g waste cooking oil by calculation. The WCOs with high free fatty acid content were difficult to transform using traditional chemical method. Table 2. Fatty acid compositions of the waste cooking oil Fatty acid

wt%

Molecular formula

Name of fatty acid

C14:0

Myristic acid

0.8±0.1

C16:0

Palmitic acid

21.2±1.1

C16:1

Palmitoleic acid

1.2±0.1

C18:0

Stearic acid

6.6±0.6

C18:1

Oleic acid

34.6±1.2

C18:2

Linoleic acid

30.6±0.8

C18:3

Linolenic acid

3.4±0.2

C20:0

Erucic acid

1.0±0.2

C22:1

Arachidic acid

0.6±0.1

Lipase from Candida sp. 99-125 belonging to Yarrowia lipolytica with enzyme activity about 20,000 IU/g was used as catalyst (Kaitai Biochemical Technology Company, Beijing, China).

D(+)-glucose

(using glucose as abbreviation in the following descriptions) and other

chemicals were of analytical grade and purchased from Beijing Chemical Factory (Beijing, China). Lipase activity assay

ACS Paragon Plus Environment

6

Page 7 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

In this research, the hydrolytic activity of lipase was quantified as the initial and residual enzyme activity. The hydrolytic activity of lipase was measured by titrimetric assay according to a modified olive oil emulsion method19, 20. Olive oil [5%, (v/v)] was emulsified in distilled water containing 2% (w/v) of polyvinyl alcohol (PVA) in a homogenizer for 6 min at maximum speed. The assay mixture consisted of emulsion (5 mL), phosphate buffer (4 mL, 100 mmol·L-1, pH 8.0) and lipase (1 mL, concentrated or diluted, depending on the quantity of lipase). The oil hydrolysis was incubated at 40 oC for 10 min with agitation. The reaction was stopped by adding 20 mL ethanol. The enzyme activity was determined by titration of the fatty acids using 50 mmol·L-1 NaOH. One activity unit of lipase was defined as the amount of enzyme that released 1 mmol of fatty acids per minute under the assay conditions. General process for methanolysis of waste oil The enzymatic synthesis was carried out in a 100 mL conical flask at a 200 rpm in an air bath rotating bed. In single factor experiment, 10 g waste cooking oil was used, and stoichiometric amounts of methanol were added six times. The reaction time and temperature were set at 24 h and under 40oC following the previous research22. The amount of lipase, glucose and water content were individually optimized, meaning that when the effect of a given parameter was investigated, the other parameters were kept constant. In response surface optimization process, 10 g waste cooking oil and stoichiometric amounts of methanol were used as substrate. Other parameters, such as lipase dosage, water content, glucose/lipase, reaction temperature and methanol adding frequency were optimized with response surface methodology. A five-factor, three-level response surface was designed by Design Expert 7.0.0 software using Box-Behnken design methodology and shown in SI Table 1 and SI Table 2. The reaction results were also analyzed by Design Expert 7.0.0 software. During the reaction process, 10µL of

ACS Paragon Plus Environment

7

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 32

reaction mixture was taken every 2 h and dissolved in 1 mL n-hexane. After centrifugation at 4800 ×g, the n-hexane solution was used for gas chromatography analysis. The data are shown as averages with their standard deviation (n = 3). Methanol tolerance assay Parallel reactions were performed under the following conditions. 10 g oil, 2.5 g methanol, and 2.5 g water were added into a 200 mL conical flask. 0.5 g lipase with 0.5 g glucose or other sugars were mixed, and then added into the system. The flasks were incubated in an air bath rotating bed with a 200 rpm stirring at 40 °C, 45 °C and 50 °C, respectively. Samples were taken every 2 h. The lipase was collected after centrifugation (4800×g). The residual activity of lipase was detected. The control experiments without glucose were performed and the lipase activity was also assayed. The data are shown as averages with their standard deviation (n = 3). Gas chromatography analysis The free fatty acid and esters, such as free acid methyl ester (FAMEs), contents in the mixture were quantified by a GC-2010 gas chromatography (GC, Shimadzu 2010 Japan). The GC analytical method was the same as reported23, 24. Heptadecanoic acid and its methyl ester purchased from Sigma were used as an internal standard. The specific compound content was defined as the weight ratio of a specific compound divided by the total weight of FFAs, FAMEs, MAGs, DAGs and TAGs, as shown in Equation (1) Specific compound content/% =

Weight of a specific compound × 100% Total weight of FFAs, FAMEs, MAGs, DAGs and TAGs …………………………………Equation (1)

ACS Paragon Plus Environment

8

Page 9 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Circular dichroism measurement The protein secondary structure of YLLIP2 was measured by circular dichroism spectrum (CD) between 180 nm and 250 nm at 25°C. The system without enzyme was used as a blank control. The secondary structure element content was estimated using the DICHROWEB application package based on the SELCON3 algorithm described by Sreerama et al.25, 26. Computational procedures Protein structure used The initial opened conformation of YLLIP2, based on the X-ray structure of YLLIP2 (PDB ID: 4JEI), was obtained in our previous study

27

.It was a semi-open conformation and was

similar with the report of Bordes et al.28. Crystal waters around the surface were preserved in the molecular dynamics simulations. Molecular Dynamics Simulations All molecular simulations were performed with the YASARA software package (version 16.3.8) using the AMBER03 force field29. The lipase was embedded in a cubic box with a 1.0 nm space left around the lipase for adding 50 % (v/v) methanol-water mixture solvent system, in which the TIP3P model was used for water molecules and methanol, α-D-glucopyranose and β-Dglucopyranose molecules were optimized by semi-empirical quantum mechanics. The rough ratio (α:β=1:1,mol/mol) of D-glucopyranose was employed due to the variable conformation of glucose in the solvent (SI Table 4). 108 α-D-glucopyranose molecules and 108 β-Dglucopyranose molecules were added into the solvent system (glucose:lipase = 1:1, w/w). The production simulations of the protein-solvent systems were performed at 313 K and at atmospheric

pressure

for

10

ns

using

the

standard

md_run

macro

(http://www.yasara.org/macros.htm) on YASARA.

ACS Paragon Plus Environment

9

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 32

In order to reduce the bias of the initial atom velocities and position of glucose molecules, three independent MD simulations were carried out. Control experiments without glucose were also performed according to the same methods. Simulation procedure The root-mean-square displacement (RMSD) of backbone atoms (the carboxylic carbon atom, the α carbon atom and the amino nitrogen atom of each amino acid residue) and the root mean square distance of all atoms of each residue from the mass center of protein (radius of gyration) were calculated to evaluate the overall motion of YLLIP2. The radial distribution function (RDF), g(r) in a system of particles, describes how the density varies as a function of distance from a reference particle. In the present work, RDF was calculated to detect the distribution of solvent molecules around YLLIP230. The dynamic cross-correlation matrix (DCCM) was calculated to analyze the movements of residues31. The values in the DCCM range from -1 (perfectly anti-correlated) to +1 (perfectly correlated). The values along the diagonal are always +1 (because the motion of an atom is perfectly correlated to itself). The 10 ns trajectory was used to calculate the above values. YASARA's definition of a hydrogen bond is that the hydrogen bond energy is better than 6.25 kJ/mol. The structural snapshot of YLLIP2 at 10 ns was obtained to view and analyze by Pymol package 32. Results and discussion Effects of glucose on biodiesel production Previously, β-CD was successfully used as additive in enzymatic synthesis of biodiesel

17

.

Recently, glucose, as an additive, was used to improve the biodiesel production. The comparative experiments for YLLIP2 lipase were undertaken both with and without glucose. The

ACS Paragon Plus Environment

10

Page 11 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

content of fatty acid methyl esters (FAMEs) was monitored to evaluate the effect of glucose in the biodiesel production. As shown in Figure 1, the content of FAMEs was 77.4 % with a lipase dosage of 1000 IU/g WCOs without glucose. However, the content reached 85.6 % under the same condition with glucose present. The results indicated that glucose could improve the methanolysis reaction with YLLIP2 for biodiesel production.

Figure 1. YLLIP2 catalyzed methanolysis both with and without glucose. (Conditions: waste oil (10 g), water content (5 wt %), lipase dosage (1000 IU/g WCOs), glucose:lipase = 1:1(w/w), 40°C, intermittent methanol addition of 6 times, during 40 h reaction.) To investigate the positive function of glucose, the methanol tolerance assay was designed. The effect of methanol and temperature on the lipase with or without glucose was checked. More methanol amount (25%wt, methanol/oil) than the methanol resistance of YLLIP233 was added to make the comparative results more obvious. As shown in Figure 2, the residual activity of

ACS Paragon Plus Environment

11

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 32

YLLIP2 without glucose was only 40 % of initial activity at 40 °C in 10 h. Whereas, YLLIP2 with glucose present could maintain 65 % of its initial activity. With increasing of incubation temperature, the residual activity of YLLIP2 when using glucose was 2.9 times and 4.6 times than that of YLLIP2 without using glucose at 45 °C and 50 °C in 10 h, respectively. Several researches indicated that methanol, especially methanol-water mixtures, inhibited the activity of lipases. The denaturation of lipase by methanol was explained that methanol molecule can penetrate into the enzyme molecules (or strip the critical water from the enzyme surface), then resulted in solvent-swollen or unfolding of the tertiary structure34, 35. On the other hand, lyophilized enzymes with sugars exhibited the desired performance in organic solvents in which the anhydrous environment presumably locks the enzyme molecule kinetically in its prior conformation. However, the lyophilization also could significantly cause protein denature because of the dehydration, and this detrimental effect can be greatly minimized or even prevented by lyophilizing enzymes in the presence of sugars36. Interestingly, in the present work, the simple addition of glucose to the reaction system could do well to the enzymatic synthesis of biodiesel with YLLIP2, indicating that glucose could help to improve the properties of YLLIP2. The results of CD analysis further proved that glucose slightly impacts the secondary structure of YLLIP2. At all tested temperatures, the α-helix of lipase with glucose was always higher than that without glucose, which indicated that the structure of protein might be more stable in the presence of glucose (SI Figure 1). In the following parts, MD simulation was employed to further investigate the effect of glucose on the methanol resistance of YLLIP2.

ACS Paragon Plus Environment

12

Page 13 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 2. Residual activity of YLLIP2 with glucose and without glucose after incubation in oleic triacylglycerol with 50% (v/v) methanol-water solvent at 40 °C, 45 °C and 50 °C for 10 h. Optimization of sugar-assisted lipase catalyzed biodiesel production To make this sugar-assisted lipase-catalyzed reaction more economically feasible for further industrial operation, the reaction conditions were optimized using single factor experiments and response surface method. Lipase amount, additive dosage, and water content had been investigated as the important parameters for biodiesel synthesis in our previous work 17. Normally, the feasibility of an enzymatic process depends on the dosage of biocatalyst. Therefore, the lipase dosage was first optimized in this research. The results indicated that the ester content increased significantly when lipase dosage increased from 10 IU/g WCOs to 40 IU/g WCOs. The content leveled off when the lipase dosage was excess to 40 IU/g oil (Figure 3 A). Based on the minimization of the cost of the process, an enzyme dosage of 40 IU/g WCOs was chosen as optimal condition.

ACS Paragon Plus Environment

13

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

The effect of lipase and additive weight ratio was shown in Figure 3 B. The content of FAMEs increased with lipase: glucose ratio (w/w) from 1:0 to 1:1, and then decreased when the glucose was excess to lipase. The results were in accordance with our work with β-CD as additive for biodiesel production17. Lower amounts of glucose than lipase were insufficient to the lipase, whereas overload of glucose might block the mass transfer between the substrate and lipase. As a result, the optimal amount of glucose dosage was chosen as lipase: glucose ratio (w/w) of 1:1. The optimal water content for enzymatic methanolysis is different for various lipases. For example, the synthesis of alkylesters by Candida antarctica lipase B decreases with increasing water content, whereas water improves the yields of reactions catalysed by Rhizopus oryzae, Candida sp. 99-125, and Pseudomonas cepacia37. In our previous work, the optimal water content for Yarrowia lipolytica Lipase 2 without any additives was 5-10 wt% (water/oil, w/w) in the methanolysis reaction22. In this work, the optimal water content was reduced to 2 wt% (water/oil, w/w) in the presence of glucose (Figure 3 C). A specific amount of water is clearly essential for maintaining the active three-dimensional structure of lipase in the reaction system38. It is assumed that glucose may replace the essential water to stabilize the lipase structure, which leads to a significant reduction of the optimal water content.

ACS Paragon Plus Environment

14

Page 15 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 3. Effects of enzyme dosage (A), glucose dosage (B) and water content (C) on methanolysis with (■) and without glucose (●). Conditions: 10 g of waste oil, 40 ° C, intermittent methanol addition of 6 times, during 24 h reaction. Based on these optimal parameters from single factor experiments, further investigations were made by response surface methodology to figure out the interaction that may exist among the individual variables. A five-factor (lipase dosage, water content, glucose/lipase weight ratio, temperature, and methanol adding frequency), three-level response surface was designed. Design Expert 7.0.0 program was used to evaluate the effect of each factor and their interactions on FAMEs content.

ACS Paragon Plus Environment

15

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 32

The experimental data obtained from Box-Behnken design were analyzed by response surface methodology. A mathematical model, following a second-order polynomial Eq. (2) which includes all interaction terms was used to calculate the predicted response: 1

#$%&'(&' = )* + , )- .- + -34

1

, )-- .-/ -34

1

+ , )-0 .- .0 -20

…………………………………Equation (2) The experimental results were illustrated in SI Table 1 and SI Table 2. The content of FAMEs varied from 50.75% to 91.69% (shown in SI Table 2), and the experiments of 2 and 3 gave the minimum and maximum contents, respectively. According to the results in SI Table 2, the model expressed can be by Eq. (3). The statistical significance of Eq. (3) was controlled by F-test and the analysis of variance (ANOVA) for response surface quadratic model and is given in SI Table 3. Values of probability (P)>F less than 0.005 indicate that model terms are significant. The relationship between predicted and experimental FAMEs content was shown in SI Figure 2.

Content/% = 85.64 + 1.57A − 0.72B − 0.72C − 1.86D + 11.87E − 4.60AB + 0.06AC + 4.03AD − 3.18AE − 0.42BC + 4.07BD − 0.46BE + 0.57CD − 0.32CE − 2.4DE − 0.90A/ − 2.32B/ + 0.42C/ − 2.31D/ − 8.75E / …………………………………Equation (3) With A (Lipase dosage), B (Water content), C (Glucose/lipase (w/w)), D (Temperature) and E (Methanol adding frequency), the optimum waste cooking oil methyl ester content was

ACS Paragon Plus Environment

16

Page 17 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

determined from the model (Eq. (3)). The optimum conditions of lipase catalyzed biodiesel synthesis from WCOs were: lipase dosage 40 IU/g oil, water content 1.95%, glucose/lipase (w/w) 1.05:1, temperature 39.4 oC and methanol adding frequency 6 times. The predicted WCOs methyl ester content was 92.3% from the model. To confirm the accuracy of the model, lipase catalyzed biodiesel synthesis from WCOs was carried out under the optimal conditions. The experimental WCOs methyl ester content was 92.1%. According to the results, verification experiments confirmed the validity of the predicted model. The reusability of the lipase in the glucose system was also tested. After each batch, the reaction mixture stood for half an hour to allow the lipase and glucose to settle down to the bottom of the reactor. Then the upper oils layer (biodiesel and oil) could be separated. Fresh WCOs was added to the reactor for the next batch reaction. The results were shown in Table 3. Following the procedures, lipase could be reused for 4 batches with the content of FAMEs over 80%. Table 3. The reuse of lipase under optimal conditions Content (%) Batches FAMEs

FFAs

MAGs

DAGs

TAGs

1

90.2±0.9

4.79±0.36

2.97±0.15

1.14±0.12

0.11±0.05

2

88.1±1.2

5.82±0.27

3.20±0.23

1.31±0.09

0.62±0.11

3

86.2±1.5

7.39±0.52

3.22±0.12

1.37±0.11

0.72±0.18

4

83.0±1.4

9.87±0.29

3.05±0.18

1.77±0.08

0.95±0.07

5

74.1±1.1

14.37±0.48

3.51±0.21

2.98±0.15

2.96±0.13

ACS Paragon Plus Environment

17

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 32

Furthermore, the process has been scaled up to a 5-ton reactor (Figure 4). Due to the insufficient mass transfer in the pilot scale reactor, the reaction time was prolonged to 40 hours to obtain a better FAMEs content and the result was shown in Figure 5. The reuse of lipase was also tested in this scaled-up production. The content of FAMEs in five batches were 91.4%, 89.5%, 83.4%, 82.5%, and 76.2% respectively. The content of FAMEs when using glucose was similar with that by using β-CD in the previous work. The cost of additive was significantly reduced since the price of glucose is much cheaper than that of β-CD, which would make the process more feasible for industrial application. After the separation process, the biodiesel properties in scaled-up experiments meet the EU biodiesel standard (EN14214), and the data were shown in SI Table 5.

Figure 4. The structure and photographic illustration of the 5 ton reactor used in the research.

ACS Paragon Plus Environment

18

Page 19 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 5. Biodiesel production with glucose as additive under optimal conditions in a 5 ton reactor. Conditions: 3.5 ton waste oil, water content 2 wt%, lipase dosage (40000 IU/kg waste oil, 0.2 wt% to the weight of oil), glucose:lipase = 1:1(w/w), 40 °C, intermittent methanol addition of 6 times, during the 40 hr reaction. Components of fatty acid methyl esters (FAMEs), free

fatty

acids

(FFAs),

monoacylglycerols

(MAGs),

diacylglycerols

(DAGs)

and

triacylglycerols (TAGs) were illustrated. In order to compare the difference of sugars on this process, a series of sugars were tested (Table 4). The results indicated that all kinds of common sugars could be used as additives in the YLLIP2 catalyzed biodiesel production. Glucose and lactose gave the highest FAMEs content improvement, following by fructose, galactose, maltose, and sucrose. Furthermore, this result validated the MD simulation results, which indicated that the protein-sugar interaction could help to improve the stability of YLLIP2. Table 4. Effect of different sugars on FAMEs content catalyzed by YLLIP2

ACS Paragon Plus Environment

19

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 32

Additives

FAMEs Content (%)

Blank Control (without additives)

( ( /% 74.3±1.4

D(+)-Glucose

90.2±0.8

D(-)-Fructose

89.9±1.1

D(+)-Galactose

89.7±1.5

α-D-Glucosido-(1→2)β-D-fructofuranoside

85.4±0.9

α-D-Glucopyranosido-(1→4)α-glucopyranose

88.2±1.2

β-D-galactopyranosyl-(1→4)α-glucopyranose

90.5±0.7

Reaction Conditions: 10g waste cooking oil was used as substrate, lipase dosage 40 IU/g oil, sugars:lipase=1:1 (w/w), water content 2%, temperature 40 oC, stoichiometric amount of methanol was added 6 times in 40 h.

Mechanism investigation by MD simulation for the glucose case RMSD is the most widely used measurement for comparing protein structures39. For each simulation system with and without glucose, three independent trajectories were presented in the supplementary information (SI Figure 3 and SI Figure 4). These trajectories (RMSD of backbone atoms vs. simulation time) reached the similar plateau before 10000 ps, indicating the convergence of the simulation calculation. Moreover, The RMSD of backbone atoms, Rg of YLLIP2, RDF of solvent system and DCCM of per-residues were calculated from the average of the three trajectories in the system with and without glucose, respectively, which could provide more reliable information for protein motion in solution than single trajectory does. In

ACS Paragon Plus Environment

20

Page 21 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

comparison to the relatively larger fluctuation of the RMSD value in the control group, the change of RMSD value in the presence of glucose was less pronounced, indicating that glucose stabilized the structure of YLLIP2. Especially, regarding the system in absence of glucose, the trajectory of RMSD exhibited the obvious upward fluctuation at 4000 ps, which was consistent with the increment of Rg at 4000 ps (Figure 6). In contrast, the Rg values remain fairly stable, comparing to the initial conformation in the system with glucose (Figure 6). The above results meant that glucose improved the overall structure stability of YLLIP2 in methanol.

Figure 6. RMSD of backbone atoms (A) Rg (B)of YLLIP2 in 50 % (v/v) methanol-water mixture solvent system with and without glucose, starting from the initial structure vs. simulation time. To understand the molecular mechanism underlying glucose stabilizing the protein structure, the distribution of the solvent on the protein surface was calculated using RDF between the solvent molecules and the protein. As shown in Figure 7 A, the high density of water molecules is distributed around the protein from about 0.1 nm to 0.5 nm, while the distribution of methanol and glucose molecules is 0.1 nm distance away from water molecules. In contrast, Figure 7 B showed the distribution of methanol in absence of glucose. The peak of methanol from about 0.1 to 0.2 nm was higher than that of water. This indicated that many

ACS Paragon Plus Environment

21

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 32

methanol molecules diffused into the water shell of protein, and some water molecules were expelled from the protein surface. These results were in agreement with other methanol-induced protein denaturation mechanisms, in which methanol molecules penetrated into the enzyme molecules or stripped the critical water from enzyme surface, thus resulting in solvent-swollen or unfolding of the tertiary structure34, 35. In the system with glucose, glucose molecules could form hydrogen bonds with methanol and then overcome the tendency of methanol molecule to accumulate in the critical water layer and diffuse to the protein chains. The number of hydrogen bonds of methanol-residue and residue-residue were calculated (Figure 8). The fewer number of hydrogen bonds between methanol-residue and more hydrogen presence between residue-residue with glucose demonstrated that glucose could prevent methanol molecules denaturing YLLIP2.

Figiure 7. RDF for solvent molecules around the surface of YLLIP2 in 50 % (v/v) methanolwater mixture solvent system, with (A) and without glucose (B).

ACS Paragon Plus Environment

22

Page 23 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

Figure 8. Average hydrogen bond number between residue-residue of YLLIP2 (A) and methanol molecule-residue of YLLIP2 (B) in 50 % (v/v) methanol-water mixture solvent system with and without glucose. In order to characterize the protection of the structure stability of YLLIP2 from glucose, four regions with obvious difference between the system both with and without glucose were pointed in the DCCM of pre-residues. The region Ⅲ and Ⅳ including the residues of the catalytic triad, and the independent stabilization of catalytic triad is important for lipase to keep its catalytic activity. As shown in Figure 9, the movement of buried residues (Glu54-Gly141, Lys20-Pro152 and Ile156-Gly240) in the regionsⅠ, Ⅲ and Ⅳ exhibited the anti-correlated to the surface regions (Lys39-Phe50, Asn261-Leu290 and Asp145-Pro152) in the presence of glucose. The result indicated that the protein internal was not disturbed by the solvent environment. However, without the protection of glucose, the movement of the core is more correlated with the movement of surface residues caused by the diffusion of the methanol (Figure 9).

ACS Paragon Plus Environment

23

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

Figure 9. DCCM of per-residues of YLLIP2 in 50 % (v/v) methanol-water mixture solvent system with (A) and without glucose (B). The values in the DCCM ranged from -1 (perfectly anti-correlated) to +1 (perfectly correlated). The values along the diagonal are always +1 (because the motion of an atom is perfectly correlated to itself). The critical regions in DCCM are pointed and presented in the structural model of YLLIP2.

ACS Paragon Plus Environment

24

Page 25 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

With respect to the abnormal regionⅡ, both two correlated regions (Arg 99-Thr 106 and Asn 107-Gln 150) are located on the surface of the protein. Hence, the interaction between glucose and regionⅡ was analyzed, and a snapshot at 10 ns was show in Figure 10. As shown in Figure 10, the glucose molecules bound the secondary structures of the protein together through hydrogen and Van der Waals bonds, and could lead to the coupled movement between the Arg 99-Thr 106 region and the Asn 107-Gln 150 region. Emil Fischer indicated the specific interaction between enzyme and glucoside in 198440. Later, Researchers have further investigated the protein-sugar interaction. Quiocho et.al presented a review of the biochemical characteristics of carbohydrate-binding sites and identified the polar residues (Asn, Asp, Gln, Glu, Arg) as the most frequently involved residues in hydrogen bonding41.

Figure 10. Graphical representation of interaction among surface residues (Asn134 and Ala111) of YLLIP2, glucose, methanol and water molecules. Speculative hydrogen bond was pointed by black dashed. Conclusion

ACS Paragon Plus Environment

25

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 32

The process using glucose as an additive in lipase catalyzed biodiesel synthesis was developed. This efficient and economic method provides a feasible enzymatic catalysis at an industrial scale. Moreover, the using of glucose provided a way to improve the methanol tolerance of YLLIP2. The MD simulation results showed that the glucose molecules not only self-assemble at the surface of the protein through hydrogen bond and Van der Waals attraction, but could also become a "member" of surface residues to defense the impaction of the environment on the structure of proteins. The present work is a successful example of combining experiments and MD simulations in order to explain and understand the protection mechanism of additives in the enzymatic process. In future work, an integrated semi-rational strategy with MD simulation could be developed to predict and guide the application of additives to improve the enzymatic reactions. ASSOCIATED CONTENT Supporting Information. Experimental setup of response surface (SI Table 1), Experimental data of response surface (SI Table 2), Analysis of variance (ANOVA) (SI Table 3), Optical rotation of

D-glucose

in

methanol-water solution(SI Table 4), Fuel properties of the biodiesel product Secondary (SI Table 5), structure measurement (SI Figure 1.), Predicted fatty acid methyl ester content versus experimental fatty acid methyl ester content (SI Figure 2.), RMSD values (SI Figure 3.) and Rg values (SI Figure 4.) (.docx)

AUTHOR INFORMATION Corresponding Author

ACS Paragon Plus Environment

26

Page 27 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

* Corresponding author: Kaili Nie Tel: +86 18611510432 Email: [email protected] Beijing Bioprocess Key Laboratory, College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, PR China Meng Wang Tel: +86 159011005745 Email: [email protected] Beijing Bioprocess Key Laboratory, College of Life Science and Technology, Beijing University of Chemical Technology, Beijing, 100029, PR China

Notes The authors declare no competing financial interest. Acknowledgements This work was supported by the National 973 Basic Research Program of China (2014CB745100), National Natural Science Foundation of China (21406011, 21676016).

ACS Paragon Plus Environment

27

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 32

REFERENCES (1) Hasan, F.; Shah, A. A.; Hameed, A. Enzyme. Microb. Technol. 2006, 39, 235-251. (2) Shimada, Y.; Watanabe, Y.; Samukawa, T.; Sugihara, A.; Noda, H.; Fukuda, H.; Tominaga, Y. In Conversion of vegetable oil to biodiesel using immobilized Candida antarctica lipase, 1999; pp 789-793. (3) Knothe, G.; Van Gerpen, J. H.; Krahl, J., The biodiesel handbook. AOCS press Champaign, 2005; pp 83-114. (4) Demirbas, A. Prog. Energ. Combust. Sci. 2005, 31, 466-487. (5) Ma, F.; Hanna, M. A. Bioresource. Technol. 1999, 70, 1-15. (6) Branco, R. J.; Graber, M.; Denis, V.; Pleiss, J. ChemBioChem 2009, 10, 2913-2919. (7) Sandoval, G.; Condoret, J.; Monsan, P.; Marty, A. Biotechnol.Bioeng. 2002, 78, 313-320. (8) Lousa, D.; Baptista, A. N. M.; Soares, C. U. M. J. Chem. Inf. Model. 2012, 52, 465-473. (9) Kulschewski, T.; Sasso, F.; Secundo, F.; Lotti, M.; Pleiss, J. J. Biotechnol. 2013, 168, 462469. (10) Santambrogio, C.; Sasso, F.; Natalello, A.; Brocca, S.; Grandori, R.; Doglia, S. M.; Lotti, M. Appl. Microbiol. Biotechnol. 2013, 97, 8609-8618. (11)

Li, M.; Yang, L. R.; Xu, G.; Wu, J. P. Bioresource. Technol. 2013, 148, 114-120.

(12) Korman, T. P.; Sahachartsiri, B.; Charbonneau, D. M.; Huang, G. L.; Beauregard, M.; Bowie, J. U. Biotechnol. Biofuels. 2013, 6, 1.

ACS Paragon Plus Environment

28

Page 29 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

(13) Dror, A.; Shemesh, E.; Dayan, N.; Fishman, A. Appl. Environ. Microbiol. 2014, 80, 15151527. (14)

Park, H. J.; Joo, J. C.; Park, K.; Yoo, Y. J. Biotechnol. Bioprocess. Eng. 2012, 17, 722-

728. (15) Dizge, N.; Keskinler, B.; Tanriseven, A. Biochem. Eng. J. 2009, 44, 220-225. (16) Dors, G.; Freitas, L.; Mendes, A. A.; Furigo Jr, A.; de Castro, H. F. Energ. Fuel. 2012, 26, 5977-5982. (17) Nie, K.; Wang, M.; Zhang, X.; Hu, W.; Liu, L.; Wang, F.; Deng, L.; Tan, T. Fuel. 2015, 146, 13-19. (18) Du, W.; Xu, Y.; Liu, D.; Zeng, J. J. Mol. Catal. B Enzym. 2004, 30, 125-129. (19)

Cao, H.; Jiang, Y.; Zhang, H.; Nie, K.; Lei, M.; Deng, L.; Wang, F.; Tan, T. Enzyme.

Microb. Technol. 2017, 96, 157-162. (20)

Lamare, S.; Sanchez, M. J.; Legoy, M. D. Ann. N. Y. Acad. Sci. 1992, 672, 171-177.

(21)

Shariat, S. S.; Jafari, N.; Tavakoli, N.; Najafi, R. B. Res. Pharm. Sci. 2015, 10, 152.

(22) Nie, K.; Xie, F.; Wang, F.; Tan, T. J. Mol. Catal. B Enzym. 2006, 43, 142-147. (23)

Munari, F.; Cavagnino, D.; Cadoppi, A.; Scientific, T. F. 2009.

(24) Liu, S.; Nie, K.; Zhang, X.; Wang, M.; Deng, L.; Ye, X.; Wang, F.; Tan, T. J. Mol. Catal. B Enzym. 2014, 99, 43-50. (25)

Whitmore, L.; Wallace, B. A. Biopolymers. 2008, 89, 392-400.

ACS Paragon Plus Environment

29

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 32

(26)

Sreerama, N.; Woody, R. W. Anal. Biochem. 2000, 287, 252-260.

(27)

Cao, H.; Deng, L.; Lei, M.; Wang, F.; Tan, T. J. Mol. Catal. B. Enzym. 2014, 109, 101-

108. (28) Bordes, F.; Barbe, S.; Escalier, P.; Mourey, L.; André, I.; Marty, A.; Tranier, S. Biophys. J. 2010, 99, 2225-2234. (29) Hess, B.; Kutzner, C.; Van Der Spoel, D.; Lindahl, E. J. Chem. Theory. Comput. 2008, 4, 435-447. (30) Berendsen, H. J. C.; Postma, J. P. M.; Gunsteren, W. F. V.; Hermans, J., Interaction Models for Water in Relation to Protein Hydration. Springer Netherlands: 1981; pp 331-342. (31) Hünenberger, P. H.; Mark, A. E.; van Gunsteren, W. F. J. Mol. Biol. 1995, 252, 492-503. (32) DeLano, W. L. The PyMOL molecular graphics system. 2002. (33) Yu, M.; Qin, S.; Tan, T. Process. Biochem. 2007, 42, 384-391. (34) Li, L.; Jiang, Y.; Zhang, H.; Feng, W.; Chen, B.; Tan, T. J. Phys. Chem. B 2014, 118, 19761983. (35) Klibanov, A. M. Nature. 2001, 409, 241-246. (36) Carpenter, J. F.; Chang, B. S.; Garzon-Rodriguez, W.; Randolph, T. W. In Rational design of stable protein formulations, Springer: US, 2002; pp 109-133. (37) Lotti, M.; Pleiss, J.; Valero, F.; Ferrer, P. Biotechnol. J. 2015, 10, 22-30. (38) Deng, L.; Tan, T.; Wang, F.; Xu, X. Eur. J. Lipid. Sci. Technol. 2003, 105, 727-734.

ACS Paragon Plus Environment

30

Page 31 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Energy & Fuels

(39) Betancourt, M. R.; Skolnick, J. Biopolymers. 2001, 59, 305-309. (40) Fischer, E. Berichte der deutschen chemischen Gesellschaft 1894, 27, 2985-2993. (41) Quiocho, F. A.; Vyas, N. K.; Spurlino, J. C. Protein-carbohydrate interactions. American Crystallographic Association, Buffalo; NY, 1991; pp 23-35.

ACS Paragon Plus Environment

31

Energy & Fuels 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 32

Table of Contents

ACS Paragon Plus Environment

32