Food Engineering Progress
Korean Society for Food Engineering
Article

Optimization of Extrusion Variables for Improving the Qualities of Textured Vegetable Protein with Green Tea Using Response Surface Methodology

Xuelian Ma, Bonyeob Gu, Gihyung Ryu*
Department of Food Science and Technology, Kongju National University
*Corresponding author: Department of Food Science and Technology, Kongju National University, Chungnam, 324-39, South Korea +82-41-330-1484; +82-41-335-5944ghryu@kongju.ac.kr

© Korean Society for Food Engineering. All rights reserved.

Received: Jan 16, 2018; Revised: Feb 4, 2018; Accepted: Feb 17, 2018

Abstract

The aim of this study was to evaluate the optimization extrusion variables on quality of textured vegetable protein by using response surface methodology. In this study, 50% soy protein isolate, 40% wheat gluten, and 10% corn starch were blended and 15% of the mixture was substituted with green tea. The moisture content (45, 50, and 55%), barrel temperature (130, 140, and 150°C), and screw speed (100, 150, and 200 rpm) were varied. A Box-Behnken design was used in this experiment. Second order polynomial regression equations were developed to relate the response to extrusion variables as well as to obtain a response surface plot. The independent variables had significant effects on the quality of the products and moisture content was the most significant. The lower moisture content led to the higher integrity index, lower nitrogen solubility index, lower water absorption capacity, higher texture, and higher cutting strength. The optimum conditions were identified as moisture content 47.78%, barrel temperature 150.00°C, and screw speed 196.05 rpm. Incorporation of green tea into protein materials could effectively improve the nutritional value of the product. Understanding these optimized extrusion variables on the product quality was useful for producing textured vegetable protein in the future.

Keywords: textured vegetable protein; response surface methodology; green tea

Introduction

Extrusion is a high temperature short time (HTST) process to produce commercially shelf stable products, because of low cost, high productivity, energy efficiency and versatility (Angela, 2007). Extrusion has been widely used in the food industry like ready-to-eat cereal breakfast, pet food, pasta, candies (Lin et al., 2007) and texturized vegetable protein (TVP).

TVP is a kind of food products made from edible protein sources and characterized by having a structural integrity and identifiable structure such that each unit will withstand hydration and cooking, and other procedures used in preparing the food for consumption. TVP is considerably similar to meat in taste, color, and structure (Asgar et al., 2010). It could partially or completely replace meat to reduce obesity, hypertension, hyperlipaemia, cardiovascular disease, diabetes, and other modern “civilized disease” (Chen, 2010). Also, it can supply protein for vegetarians or Muslims (Asgar et al., 2010). Therefore, it has drawn great attention.

Soy protein, containing eight essential amino acids, is considered the most popular food ingredients in 21st century (Yu, 2011). Soy protein isolate (SPI) is obtained by removing oil at low temperatures, which contains more than 90% protein (Yu, 2011). SPI is the ideal material for extrusion. Gluten is a wheat protein obtained when separating non-protein components from wheat flour (Sarkki& Tehtaat, 1979). Wheat Gluten (WG) is also one of the main raw materials for producing extruded tissue protein products, which can significantly improve the degree of texturization and fibrosis of texturized vegetable protein (Day et al., 2006). Green tea contains considerable amounts of catechins, which contribute to antioxidant properties (Lu et al., 2010; Shen et al., 2014). Several studies reported that green tea has anti-arteriosclerotic, anti-carcinogenic, and anti-microbial (Wang & Zhou, 2004; Sharma & Zhou, 2011).

It is reported that there are many factors affecting the texturization and formability of the product like feed rate, material moisture, barrel temperature, material pH, screw speed, or protein content (Wei et al., 2009; Maurya & Said, 2014). In this study, Box-Behken Design (BBD) in response surface methodology (RSM) was used to investigate the optimization extrusion variables (moisture content, barrel temperature, and screw speed) on integrity index, nitrogen solubility index, water absorption capacity, hardness, springiness, cohesive- ness, chewiness, and cutting strength of TVP with 15% green tea addition.

Materials and Method

Materials

Soy protein isolate (Wachsen Industry Co., Qingdao, China), wheat gluten (Roquette Freres, Lestrem, France), corn starch (Qone Co., Incheon, Korea), and green tea powder (Green tea institute, Hadong, Korea) were purchased. The moisture contents of SPI, WG, corn starch, and green tea were 6.20, 6.20, 10.87, and 6.10%, respectively. SPI, WG, and corn starch were mixed in ratio of 5:4:1 that was determined from preliminary studies. Green tea partially replaced 15% of basic formulation.

Extrusion process

The mixture was extruded with a co-rotating intermeshing twin-screw extruder (THK 31T, Incheon Machinery, Incheon, Korea) with a screw length of 690 mm and a screw diameter of 30 mm (L/D = 23:1). The screw configuration is shown in Fig. 1. Feed rate was fixed to 100 g/min during extrusion process. The extrusion variables were set up as follows: M: moisture content (45, 50, and 55%); B: barrel temperature (130, 140, and 150°C); S: screw speed (100, 150, and 200 rpm). The TVP was dried at 50°C for 12 h and was cut to 2±0.2 cm pellet to measure the physical properties. The ground sample (50-70 meshes) was used for nitrogen solubility index.

fep-22-1-1-g1
Fig. 1. Screw configuration of twin-screw extruder.
Download Original Figure
Integrity index

Integrity index was determined by slightly modified method of Park et al. (2017). Around 5 g of dried samples were soaked in 100 mL distilled water for 30 min at 80°C. The soaked samples were heated for 15 min at 121°C in autoclave (DS-PAC 60, Lab house, Anyang, Korea), then was washed with tap water for 15 s. The washed samples with 100 mL distilled water were homogenized at 14450 rpm for 1 min. The samples were filtered by 20 meshes and dried at 105°C until equilibration.

Integrity index  ( % ) = Dry residue weight initial weight × 100
(1)

Nitrogen solubility index

The soluble nitrogen was extracted from the ground sample (0.2 g) with 10 mL of 0.5% KOH. The mixture was stirred at 120 rpm for 20 min. A subsample was centrifuged at 2000 rpm for 10 min. After centrifugation, 0.1 mL of the supernatant was collected. Total nitrogen was hydrated from the 0.2 g ground sample with 1 mL of 6 N HCl solutions for 24 h at 100°C. After hydration, the mixture was re-dissolved into 10 mL of distilled water and then centrifuged at 3000 rpm for 30 min, then 0.1 mL of supernatant was collected. The soluble and total nitrogen content was analyzed by ninhydrin method. Albumin was used as a standard. nitrogen solubility index (NSI) was presented as soluble nitrogen content divided by total nitrogen content, and multiplied by 100.

Water absorption capacity

Water absorption capability (WAC) was determined by slightly modified method of Lin et al. (2002). The dried samples (5 g) were weighed and rehydrated in 100 mL distilled water for 2 h at 80°C. The rehydrated samples were drained for 5 min by using 20 meshes screen in room temperature and then weighed. The WAC was calculated according to the following equation:

WAC  ( % ) = [ ( W 2 W 1 ) / W 1 ] × 100
(2)

where W1 and W2 are weights before and after rehydration, respectively.

Texture profile analysis

The texture profile analysis (TPA) of TVP was determined by a Texture Analyzer (Compac-100II, Sun Sci. Co., Tokyo, Japan). The TVP samples (2.0×1.0×1.0 cm) were rehydrated for 30 min at 80°C and drained for 10 min in room temperature by using 20 meshes sieve. Hardness, springiness, cohesiveness, and chewiness of TPA were recorded at 60 mm/min travel speed and 10 kg load by using cylindrical probe with 25 mm diameter (Park et al., 2017). Ten replications were measured.

Cutting Strength

Cutting strength was determined by a Texture Analyzer (Compac-100II, Sun Sci. Co., Tokyo, Japan). The TVP samples (2.0×1.0×1.0 cm) were rehydrated for 30 min at 80°C and drained for 10 min in room temperature. For cutting strength, a knife probe with 2 kg load and 60 mm/min were used to cut the sample. Ten samples were measured.

Experimental design and statistical analysis

The optimization experiment was carried out using RSM with a BBD experimental design of three levels and three factors with three replicates at the center point. Different extrusion processes and levels are shown in Table 1. The experimental data for all of responses were analyzed using Design-expert version 8.06 (Stat-Ease, Inc., Minneapolis, Minnesota, USA) software. Significance of results was determined at the p<0.05 level for all data analysis. Experimental model can be analyzed by stepwise and expressed with the following equation:

Table 1. Factors and levels of experiment
Factors Symbol Values
-1 0 1
Moisture content (%) M 45 50 55
Barrel temperature (°C) B 130 140 150
Screw speed (rpm) S 100 150 200
Download Excel Table

Y=β 0 + β 1 M+β 2 B+β 3 S+β 12 MB+β 013 MS+β 23 BS+β 11 M 2 + β 22 B 2 + β 33 S 2 + ε
(3)

with Y representing the estimated response, β0 representing the equation parameters for the constant term, βi representing the linear terms, βij (i,j = 1-3) representing the cross terms, βii representing the quadratic terms for a single variable, and ε representing the random error.

Results and Discussion

The results on the quality of integrity index, nitrogen solubility index, water absorption capacity, springiness, cohesiveness, chewiness, and cutting strength of TVP are shown in Table 2. The coefficient and analysis of variance (ANOVA) at p<0.05 of TVP is shown in Table 3. In order to understand the combined effects of two individual treatment factors on all of responses, the response surface plot is created as the function of two independent variables, while keeping the third one which is the most insignificant factor at the center.

Table 2. Experimental results of texturization of TVP
Run M (%) B (° C) S (rpm) II (%) NSI (%) WAC (%) Hardness (g) Springiness (%) Cohesiveness (%) Chewiness (g) Cutting strength (g/cm 2)
1 1(55) 0(140) 1(200) 1.70 40.37 268.84 1307.21 82.20 72.84 996.95 361.23
2 0(50) 0(140) 0(150) 38.82 10.16 298.56 2088.34 81.10 66.22 1363.09 529.45
3 -1(45) 1(150) 0(150) 55.95 7.19 284.02 2533.49 82.51 69.33 1736.61 678.46
4 0(50) -1(130) 1(200) 66.33 4.91 215.50 3564.38 79.51 61.47 2415.98 745.21
5 -1(45) 0(140) 1(200) 51.31 4.62 266.69 2359.24 76.88 62.18 1622.83 829.52
6 1(55) 0(140) -1(100) 1.64 42.12 288.88 2041.11 79.54 64.83 1229.36 261.99
7 0(50) 0(140) 0(150) 63.89 4.12 250.44 2540.00 74.61 60.87 1384.18 508.16
8 0(50) 1(150) 1(200) 61.35 5.80 259.06 2362.13 79.89 67.88 1570.82 738.21
9 1(55) -1(130) 0(150) 1.12 43.03 245.71 2891.74 81.36 62.44 1831.32 296.69
10 1(55) 1(150) 0(150) 5.78 36.94 290.74 1606.68 75.03 69.80 1126.71 403.44
11 -1(45) -1(130) 0(150) 65.44 4.95 206.41 4453.74 78.44 57.73 2985.91 378.46
12 0(50) -1(130) -1(100) 55.06 7.55 229.25 2926.45 77.00 55.79 1572.24 504.45
13 0(50) 1(150) -1(100) 56.86 7.58 265.38 2453.61 72.19 67.10 1532.76 805.91
14 -1(45) 0(140) -1(100) 14.08 27.02 237.22 2585.70 76.39 59.97 1803.75 429.69
15 0(50) 0(140) 0(150) 36.16 11.65 270.89 1761.00 77.94 66.25 1248.90 530.40

M=Moisture content; B=Barrel temperature; S=Screw speed; II=Integrity index; NSI=Nitrogen solubility index; WAC=Water absorption capacity.

Download Excel Table
Integrity index

The values of integrity index of TVP are presented in Table 2. The highest value of integrity index was 66.33% when extrusion variables were 50% moisture content, 130°C barrel temperature, and 200 rpm screw speed. The coefficient of predicted model and ANOVA evaluation for integrity index of TVP are shown in Table 3. The independent variables moisture content and quadratic terms of moisture content and barrel temperature significantly (p<0.05) affected integrity index. The regression equation could be expressed as:

Integrity index  ( % ) = 1301.24 + 94.86 M-50.21B-0 .99M 2 + 0.18 B 2 R 2 = 0.8630
(4)

According to Fig. 2, the value of integrity index significantly (p<0.05) decreased with increasing moisture content. This result was in agreement with the description of Lin et al. (2002) and Park et al. (2017), who reported that texturization of TVP increased as the moisture content decreased. This might be because increase in friction and shear due to lower moisture content. Therefore, it led to higher degree of texturization and fiber formation (Lin et al., 2002).

fep-22-1-1-g2
Fig. 2. Response surface plot of integrity index of TVP at 150 rpm screw speed.
Download Original Figure
Nitrogen solubility index

The values of NSI of TVP ranged from 4.12 to 43.03% (Table 2). As shown in Table 3, moisture content and screw speed presented significant effect on NSI. The interaction between moisture content with screw speed also significantly (p<0.05) affected NSI. In addition, moisture content presented quadratic effect on NSI with the regression equation of:

NSI ( % ) = 1862.93 7366 M-1.10S-+0.02M × S+0 .74M 2 R 2 = 0.9399
(5)

Table 3. Coefficient of prediction models and ANOVA evaluation for each response variable
Source II (%) NSI (%) WAC (%) Hardness (g) Springiness (%) Cohesiveness (%) Chewiness (g) Cutting strength (g/cm 2)
Coefficient p-value Coefficient p-value Coefficient p-value Coefficient p-value p-value Coefficient p-value Coefficient p-value Coefficient p-value
estimate Prob> F estimate Prob>F estimate Prob>F estimate Prob>F Prob > F estimate Prob>F estimate Prob>F estimate Prob>F
Model 1301.25 0.0006 1862.93 <0.0001 -223.23 0.0011 16300.15 0.0149 0.2105 -32.01 0.0003 101800.74 0.0003 -14906.73 0.0066
M 94.86 0.0004 -73.66 <0.0001 2.5 0.0427 -102.14 0.0247 0.6445 0.52 0.0106 -74.12 0.0018 619.89 0.0088
B -50.21 0.8125 2.53 0.0007 -61.00 0.0100 0.2543 0.46 0.0002 -1346.18 0.0024 8.77 0.0431
S -1.1 0.0496 -1.03 0.7972 0.1785 0.04 0.0308 -10.08 0.0495
MB 0.2282
MS 0.02 0.0465 0.3067
BS 0.7021
M 2 -0.99 0.0036 0.74 <0.0001 0.2338 -6.45 0.0082
B 2 0.18 0.0201 0.0858 4.68 0.0047
S 2 0.0951 0.04 0.0372
Lack of Fit 0.118 0.4875 0.6273 0.3462 0.8786 0.8083 0.0643 0.519

M=Moisture content; B=Barrel temperature; S=Screw speed; II=Integrity index; NSI=Nitrogen solubility index; WAC=Water absorption capacity; Nonsignificant factors (p>0.05) were removed.

Download Excel Table

NSI values significantly (p<0.05) decreased with the increase in moisture content and screw speed (Fig. 3). It was consistent with the reported of Park et al. (2017), who claimed NSI value increased as moisture increased from 50 to 55% when SPI and WG were 4:1. This may be because of the high moisture content which decreased protein denaturation and increased in protein solubility (Wei et al., 2009).

fep-22-1-1-g3
Fig. 3. Response surface plot of nitrogen solubility index of TVP at 50% moisture content.
Download Original Figure
Water absorption capacity

The WAC of TVP is shown in Table 2. The highest value of WAC appeared at 50% moisture content, 140°C barrel temperature, and 150 rpm screw speed. A linear coordination of WAC with moisture content and barrel temperature was observed (Table 3). The equation is expressed as follows:

WAC ( % ) = 223.23 + 2.49 M+2 .53B R 2 = 0.7090
(6)

WAC significantly (p<0.05) increased with the increase in barrel temperature and moisture content (Fig. 4). Sun & Muthukumarappan (2002) and Yu (2011) also demonstrated that higher barrel temperature resulted in higher rehydration rate. Sun& Muthukumarappan (2002) and Seker (2005) didn’t find any significant effect of screw speed on rehydration rate of TVP. Yu (2011) reported that high barrel temperature probably creates more air cells in product. This may cause increased absorption of water during rehydration of TVP, subsequently, the WAC of TVP increased with the increasing barrel temperature.

fep-22-1-1-g4
Fig. 4. Response surface plot of water absorption capacity of TVP at 150 rpm screw speed.
Download Original Figure
Hardness

Hardness is important quality attributes of un-expansion extruded products and depend on structural characteristics arising from extrusion condition. As shown in Table 3, analysis of variance indicated that the extrusion moisture content and barrel temperature were significant to hardness of product. The regression model for process variables and product hardness is shown:

Hardness  ( g ) = 16300.15 102.14 M-61.00B-1 .03S R 2 = 0.5996
(7)

As barrel temperature and moisture content increased, the hardness of products decreased sharply (Fig. 5a). Kang (2007) reported that increasing moisture content within the range of 35-55% could reduce the hardness. We got similar trend with Kang (2007). The samples extruded at 55% moisture had the lowest hardness, which were likely due to more water contained within the samples (Lin et al., 2000).

fep-22-1-1-g5
Fig. 5. Response surface plot of hardness (a) and chewiness (d) of TVP at 150 rpm screw speed; response surface plot of cohesiveness (c) and springiness (b) of TVP at 50% moisture content.
Download Original Figure
Springiness

The values of springiness of TVP are presented in Table 2. The springiness ranged between 72.19 and 82.51%. The coefficient of predicted model and ANOVA evaluation of springiness are shown in Table 3. The ANOVA study results showed that model of springiness was not significant meanwhile none of experimental variables exhibited a significant effect on springiness (Table 3). These extrusion variables might not be high enough to significantly change springiness. The springiness response surface plot of TVP is shown in Fig. 5(B). Springiness increased with increasing barrel temperature from 130 to 140°C while it decreased with the increase in barrel temperature from 140 to 150°C. That might be because materials did not reach a molten state under lower barrel temperature. However, the moisture of product was evaporated rapidly under higher barrel temperature, which led to lower springiness.

Cohesiveness

The experimental values of cohesiveness of TVP products under different extrusion conditions are shown in Table 2. Cohesiveness ranged from 55.79 to 72.84%. Minimum cohesiveness (55.79%) occurred at 50% moisture content, 130°C barrel temperature, and 100 rpm screw speed, while the maximum (72.84%) was obtained at 55% moisture content, 140°C barrel temperature, and 200 rpm screw speed (Table 2). Table 3 shows the coefficient of predicted model and ANOVA evaluation for cohesiveness of TVP. After removing all of insignificant terms, the regression coefficients generated a linear equation:

Cohesiveness ( % ) = 32.01 + 0.52 M+0.46B+0 .04S R 2 = 0.8042
(8)

The response surface plot between cohesiveness and moisture content with barrel temperature at 150 rpm screw speed is presented in Fig. 5(C). The results showed that the cohesiveness increased with an increase in moisture content and barrel temperature. Wei et al. (2006), who investigated the effects of extrusion variables on physical properties of defatted soy meal, reported that cohesiveness increased with the increase in moisture content. It was consistent with our results.

Chewiness

The results of chewiness for different experimental conditions are shown in Table 2. The chewiness of TVP ranged between 996.95 g (55% moisture content, 140°C barrel temperature, and 200 rpm screw speed) and 2,985.91 g (45% moisture content with 130°C barrel temperature and 150 rpm screw speed). The regression model was significant (Table 3). Moisture content, barrel temperature, and the quadratic term of barrel temperature significantly (p<0.05) affected the chewiness (Table 3). The fitted model was created as follow:

Chewness  ( g ) = 101800.74 74.12 M-1346.18B+4 .68B R 2 = 0.8893
(9)

Increasing moisture content and barrel temperature, the chewiness significantly decreased (Fig. 5d). Sun (2009), who studied textured soy protein produced by high moisture extrusion technology, claimed that chewiness decreased with increasing moisture content. This was consistent with our result.

Cutting strength

The experimental values of cutting strength of TVP products are shown in Table 2. The highest cutting strength value was 829.52 g/cm2 when the extrusion condition was 45% moisture content, 140°C barrel temperature, and 200 rpm screw speed. As shown in Table 3, the moisture content, barrel temperature, and screw speed were the significant factor. In addition, moisture content, screw speed was also presented the quadratic effect with the regression equation as:

Cutting strength  ( g/cm 2 ) = 14906.73 + 619.89 M+8.77B-10.08S-6 .45M 2 + 0.04 S 2  R 2 = 0.7929
(10)

According to Fig. 6, the cutting strength increased with increasing moisture content from 45 to 50% while it decreased with increasing moisture content from 50 to 55%. The higher barrel temperature resulted in higher cutting strength. The sample was cut in cross-linking direction. Higher cutting strength stands for higher integrity index of products, which indicates that it has better texturization degree. Barrel temperature and moisture content could influence the materials viscosity in the barrel and affect the friction between materials and screw. Thereby they seemed to influence the maturation of materials in barrel, ultimately affecting the texture characteristics of the product.

fep-22-1-1-g6
Fig. 6. Response surface plot of cutting strength of TVP at 150 rpm screw speed.
Download Original Figure

Conclusion

The optimization of three extrusion variables on the quality of SPI, WG, corn flour, and green tea blends were studied. Three important extrusion parameters (moisture content, barrel temperature, screw speed) had different effects. None of the variables presented significant effect on springiness, but hardness, cohesiveness, cutting strength was significant affected by all the three independent variables. Moisture content and barrel temperature exhibited significant quadratic effect on integrity index.

Base on integrity index results, the optimum conditions of TVP given by the model as follows: moisture content 47.78%, barrel temperature 150.00°C, and screw speed 196.05 rpm. Under this condition, the predicted values of integrity index, nitrogen solubility index, water absorption capacity, hardness, cohesiveness, chewiness, and cutting strength by the model were 64.30%, 4.51%, 275.24%, 2067.97 g, 69.68%, 1632.29 g, 863.47 g/cm2, respectively.

Acknowledgment

This research was supported by Basic Science Research Program through the National Research Foundation (NRF) funded by the Ministry of Education, Science and Technology (2016H1C1A1035974).

References

1.

Angela C. 2007 Protein-protein interaction of soy protein isolate from extrusion processing. MS. thesis. Missouri-Columbia Univ. .

2.

Asgar M.A., Fazilah A., Huda N., Bhat R., Karim A.A. Nonmeat protein alternatives as meat extenders and meat analogs.. Compr. Rev. Food Sci. Food Saf. 2010; 9:513-529.

3.

Chen F.L. 2010 Role of water in extrusion texturization of vegetable protein. Ph.D. thesis. CAAS. .

4.

Day L., Augustin M.A., Batey I.L. Wheat-gluten uses and industry needs.. Trends Food Sci. Technol. 2006; 17:82-90.

5.

Kang L.N. 2007 Texturization technology and mechanism of soybean protein by high moisture extrusion. Doctor Dissertation of Northwest A&F University in China. .

6.

Lin S., Huff H.E., Hsieh F. Extrusion process parameters, sensory characteristics, and structural properties of a high moisture soy protein meat analog.. J. Food Sci. 2002; 67:1066-1072.

7.

Lu T.M., Lee C.C., Mau J.L., Lin S.D. Quality and antioxidant property of green tea sponge cake.. Food Chem. 2010; 125:186-192.

8.

Maurya A.K., Said P.P. Extrusion processing on physical and chemical properties of protein rich products-an overview.. J. Bioresour. Eng. Technol. 2014; 1:67-73.

9.

Park J.H., Chatpaisarn A., Ryu G.H. Effects of gluten and moisture contents on texturization of extruded soy protein isolate.. J. Korean Soc. Food Sci. Nutr. 2017; 46:473-480.

10.

Sarkki M.L., Tehtaat R. Food uses of wheat gluten.. J. Am. Oil Chem. Soc. 1979; 56:443-446.

11.

Seker M. Selected properties of native or modified maize starch/soy protein mixtures extruded at varying screw speed.. J. Sci. Food Agric. 2005; 85:1161-1165.

12.

Sharma A., Zhou W.B. A stability study of green tea catechins during the biscuit making process.. Food Chem. 2011; 126:568-573.

13.

Shen X.J., Han J.Y., Ryu G.H. Effects of the addition of green tea powder on the quality and antioxidant properties of vacuum-puffed and deep-fried Yukwa.. Food Sci. Technol. 2014; 55:362-367.

14.

Sun Y., Muthukumarappan K. Changes in functionality of soy-based extrudates during single-screw extrusion processing.. Int. J. Food Prop. 2002; 5:379-389.

15.

Sun Z.X. 2009 Study on textured soy protein produced by high moisture extrusion technology. MS. thesis. Heilongjiang. .

16.

Wang H.W., Ma L.Q., Zhou J.G., Lin B.J. Effect of soy protein raw materials system on extrusion texturization.. J. Chin. Inst. Food Sci. Technol. 2002; 2:33-38.

17.

Wang R., Zhou W. Stability of tea catechins in the bread making process.. J. Agric. Food Chem. 2004; 52:8224-8229.

18.

Wei Y.M., Kang L.N., Zhang B. Processing and mechanism of high moisture textured soy protein.. Agric. Eng. 2006; 22:193-197.

19.

Wei Y.M., Kang L.N., Zhang T. Food extrusion theory and technology. Bejing, China: Light Industry Press. 2009; p. 49-87.

20.

Yu L. 2011 Extrusion processing of protein rich food formulations. Ph.D. thesis. McGill Univ. .