versão On-line ISSN 1678-4324
Braz. arch. biol. technol. vol.54 no.6 Curitiba nov./dez. 2011
FOOD/FEED SCIENCE AND TECNOLOGY
Ganga S. Meena; Suneel Gupta; Gautam C. Majumdar; Rintu Banerjee*
Microbial Biotechnology and Downstream Processing Laboratory; Department of Agricultural and Food Engineering; Indian Institute of Technology; Kharagpur-721 302 - India
The aim of this work was to optimize the biomass production by Bifidobacterium bifidum 255 using the response surface methodology (RSM) and artificial neural network (ANN) both coupled with GA. To develop the empirical model for the yield of probiotic bacteria, additional carbon and nitrogen content, inoculum size, age, temperature and pH were selected as the parameters. Models were developed using ¼ fractional factorial design (FFD) of the experiments with the selected parameters. The normalized percentage mean squared error obtained from the ANN and RSM models were 0.05 and 0.1%, respectively. Regression coefficient (R2) of the ANN model showed higher prediction accuracy compared to that of the RSM model. The empirical yield model (for both ANN and RSM) obtained were utilized as the objective functions to be maximized with the help of genetic algorithm. The optimal conditions for the maximal biomass yield were 37.4 °C, pH 7.09, inoculum volume 1.97 ml, inoculum age 58.58 h, carbon content 41.74% (w/v), and nitrogen content 46.23% (w/v). The work reported is a novel concept of combining the statistical modeling and evolutionary optimization for an improved yield of cell mass of B. bifidum 255.
Key words: Probiotics, response surface methodology (RSM), FFD, artificial neural network (ANN), genetic algorithms (GA)
Bifidobacterium is the most prominent member of plethora class of bacterial species with probiotic properties. The popularity of this group of bacteria is based on the millennia of use in the food and feed that are used in the probiotic dairy drinks and yoghurts since long (Sanders, 1999). At present, in India, the production of probiotics is reported to grow annually about 22.6 % until 2015 and the market of the probiotics is ~20.6 million rupees (€320,000). The market demand indicates that it is economically viable product. The probiotics have immense application in the food/healthcare sector. There are plenty of industries venturing into the production and selling of the probiotics sachets to meet the increasing demand. Most common bacteria targeted by the industries for the probiotic sachet preparation includes Bifidiobacterium. Microbial colonization of the human intestine starts immediately after the birth (Gibson and Roberfroid, 1995). The predominant bacteria at the infancy stage are Bifidobacteria which colonize within the first 4-7 days of birth with the numbers ranging from 109-1010 CFU/g of faeces in breast-fed infants (Gismondo et al., 1999). Bifidobacterium sp. is one of the major microorganisms in the gastrointestinal tract flora of the children and adults. These bacteria have a strong stimulatory effect for the normal development of microbiota and maturation of gut associated lymphoid tissue (Schezenmeir and De Vrese, 2001). Probiotic bacteria such as Bifidobacterium and Lactobacillus sp. in the gastrointestinal tract can play an important role in promoting the human health (Savage, 1977; Mitsuoka, 1990). These microorganisms can contribute to digestion, immune stimulation and inhibition of the pathogens such as Bacteroides, Escherichia, Clostriduim and Proteus which are potentially harmful bacteria found in the gastrointestinal tract (Ziemer and Gibson, 1998).
The primary mechanism for probiotic action is known as competitive colonization or competitive suppression. It is best described as the proliferation of the probiotic bacteria in the human intestine, leaving little space for the growth of any pathogens (Ballongue, 1992; Biavati et al., 2000). To develop the growth model of probiotic bacteria through the traditional method, i.e. one variable at-a-time is time consuming and interactions of different variables can also affect the yield. Unlike the conventional optimization, the statistical optimization methods can take into account the interactions of the variables in generating the process response. Process optimization through the statistical method is a technique in which changes or adjustments are made in a process to get better results (Myers and Montgomery, 2002). There are several techniques for process optimization, i.e., Response Surface Methodology (RSM), Artificial Neural Networks (ANN), Genetic Algorithms (GA), etc. In these engineering applications, a response of interest is usually influenced by several variables and the objective of the engineering applications is to find the variables that can optimize the response. RSM is a tool on that basis we find the optimal process parameters that produce a maximum or minimum value of the response and represent the direct and interactive effects of the process parameters through two and three-dimensional plots (Gangadharan et al., 2008). Artificial neural networks are computational models of nervous systems. Natural organisms, however, do not possess only nervous systems but also genetic information stored in the nucleus of their cells (genotype). The nervous system is part of the phenotype which is derived from this genotype through the process of development (Rajasekaran and Vijaylakshmi, 2004). Using the method of neural networks (NN), the relationship between a set of independent variables X and the dependent variables Y can be obtained. From the given pairs of input X and output Y data, neural network directly learns, and then develops a relationship between them but does not yield any mathematical equation relating the variables. After the learning, this network is able to predict the correct output from an input data set that has not been previously used during the learning. Genetic algorithms (GA) are a tool by which the optimization problems can be accurately solved within a limited use of computer time (Das, 2005). The objective of this work was to optimize and improve the yield of probiotic bacteria, Bifidobacterium bifidum by optimizing the growth parameters such as temperature, pH, inoculum volume, inoculum age and additional effect of different carbon and nitrogen sources with the help of Response Surface Methodology, Artificial Neural Network and Genetic Algorithms.
MATERIALS AND METHODS
Organism and growth condition
Pure culture of Bifidobacterium bifidum 255 was obtained from the National Collection of Dairy Cultures (NCDC) Karnal, Haryana (India). The culture was grown in a modified MRS media containing 1% (w/v) sodium thiosulphate at 30ºC under anaerobic condition. Biomass growth was determined by measuring the optical density (OD) at 600 nm.
Selection of initial parameters
For the selection of initial parameters, 'one variable at a time method' was used. The different variables viz. temperature, pH, volume of inoculum, age of inoculum and additional carbon and nitrogen sources were selected for growth of B. bifidum.
Empirical model development
To find out the effect of different growth parameters on the predicted value of the bacterial growth, Yp was obtained by conducting the experiments on different combination of independent variables (growth parameters), which was obtained from a standard experimental design. During the experiments, the 'response' or values of 'dependent variables' obtained from each of the combinations of independent variables was measured. A mathematical relationship between the independent and dependent variables was developed. This relationship was called 'model'. Using this model, the predicted values of responses were found out within the domain of limiting values of independent variables. For the different growth parameters, a polynomial model was developed between the growth and growth parameters to find out the following relationship between the coded values x1, x2, x3, x4, x5 and x6 of independent variables and dependent variable Yp as shown below
Where bo, b1, b2.........etc. are the regression constants.
Fractional factorial design
Using two levels (+1 and -1) factorial design, two values of l and s for two sacrificing interactions were l1, s1, l2 and s2. With the help of factorial design, s values were identified as (s1= 0, s2 =0), (s1= 0, s2 =1), (s1= 1, s2 =0), and (s1= 1, s2 =1). In this study, all the experiments were conducted according to s1= 0 and s2 =0 design.
Neural Network modeling
ANN chosen was a radial basis function network with supervised learning. The model was based on feed forward back propagation training method. In this process, the network computed the error between the desired output (predicted) and the actual (experimental) output. It trained the network to make adjustments to minimize the error and back propagate the same.
In this optimization study, GA was applied to the developed ANN based model as shown in the Fig 2. The prime objective of this study was to maximize the biomass yield of gap B. bifidum by monitoring the growth parameters such as temperature, pH, inoculum volume, inoculum age, carbon % and nitrogen %. It was posed as the minimization of problem associated with the optimization studies. Genetic optimization continued till the termination condition i.e. maximum biomass yield was obtained.
For proper execution of ANN and GA, MATLAB 7.0 was used to develop the empirical model.
RESULTS AND DISCUSSION
Selection of initial parameters
Fig 3 (A-F) shows the effect of temperature, initial pH, initial inoculum volume, initial incubation period, supplementation of additional carbon and nitrogen sources on the growth of the bacterial culture. All these parameters, their variation and optimum values are given in Table 1.
Empirical model development
From the above results, the maximum and minimum values of six independent parameters for B. bifidum were fixed as shown in Table 2. For developing the model between coded values x1, x2, x3, x4, x5, x6 of independent variables and dependent variable Yp, the experiments were conducted according to the fractional factorial design. All these combinations have been given in Table 3 with their corresponding l and s values.Various combination of process variable found at s1=0, s2=0 is shown in the Table 4 with their experimental value Ye for the growth of B. bifidum.
The experimental data were fitted to the full quadratic equation. The design matrix and the fitness of each term were analyzed by means of the ANOVA (Kumari et al., 2008). Figure 4 shows the corresponding model coefficients (R2 0.840) together with the regression coefficient of determination, which is a measure of how well the regression model can be made to fit the raw data.
A self-organizing feature map network was used to predict the growth condition parameters. Different factors, viz. temperature, pH, inoculum volume, Inoculum age, additional carbon and nitrogen sources were used as each unit of input layer. The output layer was composed of one response variable, the growth of B. bifidum. A set of factors was used for training and fed into the computer. Several iterations were conducted with different numbers of neurons of hidden layer in order to determine the optimal ANN structure.
The optimum number of neurons in the hidden layer was iteratively determined by changing the number of neurons. This was started with two neurons and the number of neurons was increased up to six. The least MSE value and a good prediction of the outputs of both training and validation sets were obtained with four neurons in the hidden layer (Dutta et al., 2004). The R2 value between the actual and estimated responses was determined as 0.930 (Fig. 5). In ANN modeling, the replicates at center point did not improve the rediction capability of the network because of the similar inputs.
Using MATLAB 7.0, the constants of regression equation and predicted value of dependent variable (OD) were found out. The 'model' which was obtained for B. bifidum 255 is given below.
The predicted value of independent variable and corresponding experimental value for B. bifidum 255 is shown in Table 5. Genetic algorithms were applied on the data obtained from the neural network using MATLAB 7.0 The optimum values or the combination of different process parameters on which the bacterial growth measured by the optical density (OD) was maximum for B. bifidum which is given in the Table 6.
There are several reports on the optimization of growth of the probiotic bacteria which are very close to the present result. Kiviharju et al. (2005) reported maximum production of B.longum at 40ºC. Ram and Chander, (2003) reported maximum growth of Bifidobacteria at 37 ºC and pH 7.0. Laxmi et al. (2011) reported the addition of carbon and nitrogen sources for enhanced growth of Bifidobacterium sp. In the present study, the RSM/ANN coupled with GA methodology resulted in an enhanced biomass yield. This is a new approach not reported earlier. However, optimization studies based on the ANN-GA for improved performance of biological systems have been reported earlier by Haider et al. (2008) and Sivapathasekaran et al. (2010).
In the present study, MATLAB 7.0 was used to fit the experimental values into a regression equation which predicted the yield of B. bifidum 255. The RSM and ANN methodologies coupled with GA were used for optimizing the input parameters. Both the models provided similar quality predictions for the above independent variables in terms of the growth conditions with ANN with more accuracy in estimation. The regression coefficients (R2) of ANN and RSM were 0.9368 and 0.8838, respectively, which clearly reflected that the ANN was better than RSM. The optimum values obtained after the GA study were 37.4°C, pH 7.09, inoculum volume 1.97 ml, inoculum age 58.58 h, carbon content 41.74% (w/v), nitrogen content 46.23% (w/v), resulting the maximum yield of probiotic bacteria. It was further noticed that ANN coupled with GA was the best combination for model development of B. bifidum.
Ballongue, J. (1992), Bifidobacteria and probiotic action probiotics, the scientific basis. J. Dairy Sci., 7, 357-413. [ Links ]
Biavati, B.; Vescovo, M.; Torriabi, S. and Bottazzi, V. (2000), Bifidobacteria: history, ecology, physiology and applications. Analyt. Microbiol., 50, 117-131. [ Links ]
Das, H. (2005), Hand book of food processing operations analysis. Asian book private limited, New Delhi. [ Links ]
Dutta, J.R.; Dutta, P.K. and Banerjee, R. (2004), Optimization of culture parameters for extracellular protease production from a newly isolated Pseudomonas sp. using response surface and artificial neural network models. Process Biochem., 39, 2193-98. [ Links ]
Gangadharan, D.; Sivaramakrishnan, S.; Nampoothiri, K.M.; Sukumaran, R.K. and Pandey, A. (2008), Response surface methodology for the optimization of alpha amylase production by Bacillus amyloliquefaciens. Bioresource Technol., 99, 4597-02. [ Links ]
Gibson, G.R. and Roberfroid, M. (1995), Dietary modulation of the human colonic microbiota: introducing the concept of prebiotics. J. Nutr., 125, 1401-1412. [ Links ]
Gismondo, M.R.; Drago, L. and Lombardi, A. (1999), Review of probiotics available to modify gastrointestinal flora. Int. J. Antimicrob. Agents, 12, 287-292. [ Links ]
Goldberg, D. (1989), Genetic Algorithms in search, Optimization and machine learning. Pearson Education, Asia. [ Links ]
Gulati, T.; Chakrabarti, M.; Singh, A.; Duvuuri, M. and Banerjee, R. (2010), Comparative study of response surface methodology, artificial neural network and genetic algorithms for optimization of Soybean hydration. Food Technol. Biotech., 48, 11-18. [ Links ]
Haider, M.A.; Pakshirajan, K.; Singh, A. and Chaudhury, S. (2008), Artificial neural network-genetic algorithm approach to optimize media constituents for enhancing lipase production by a soil microorganism. Appl. Biochem. Biotechnol., 144, 225-235. [ Links ]
Kiviharju, K.; Leisola, M. and Eerikäinen, T. (2005), Optimization of a Bifidobacterium longum production process. J. Biotechnol., 117, 299-308. [ Links ]
Kumari, K.S.; Babu, I.S. and Rao, G.H. (2008), Process optimization for citric acid production from raw glycerol using response surface methodology, Indian J. Biotech., 7, 496-501. [ Links ]
Laxmi, N.P.; Mutamed, M.A. and Nagendra, P.S. (2011), Effect of carbon and nitrogen sources on growth of Bifidobacterium animalis Bb12 and Lactobacillus delbrueckii ssp. bulgaricus ATCC 11842 and production of β-galactosidase under different culture conditions. Int. Food Res. J., 18, 373-380. [ Links ]Mitsuoka, T. (1990), Bifidobacteria and their role in human health. J. Ind. Microbiol., 6, 263-268. [ Links ]
Myers, R.M. and Montgomery, D.C. (2002), Response surface methodology. John Wiley and Sons, Inc., New York. [ Links ]
Rajasekaran, S. and Vijaylakshmi, P.G.S. (2004), Neural networks fuzzy Logic and genetic algorithms. Prentice Hall of India, New Delhi. [ Links ]
Ram, C. and Chander, H. (2003), Optimization of culture conditions of probiotic bifidobacteria for maximal adhesion to hexadecane. World J. Microbiol. Biotechnol., 19, 407-410. [ Links ]
Sanders, M.E. (1999), Probiotics. Food Biotechnol.,53, 67-77. [ Links ]
Savage, D. (1977), Microbial ecology of the gastrointestinal tract. J. Microbiol., 31, 107-133. [ Links ]
Schezenmeir, J. and De Vrese, M. (2001), Probiotics, prebiotics, and synbiotics -approaching a definition. Am. J. Clin. Nutr., 73, 361-364. [ Links ]
Sen, R. and Babu K.S. (2005), Modeling and optimization of the process conditions for biomass production and sporulation of a probiotic culture. Process Biochem., 40, 2531-38. [ Links ]
Sivapathasekaran, C.; Mukherjee, S.; Ray, A.; Gupta, A. and Sen, R. (2010). Artificial neural network modeling and genetic algorithm based medium optimization for the improved production of marine biosurfactant. Bioresource Technol., 101, 1884-87. [ Links ]
Ziemer, C.J. and Gibson, G.R. (1998), An Overview of probiotics, prebiotics and synbiotics in the functional food concept: perspectives and future strategies. Int. Dairy. J., 8, 473-479. [ Links ]
Received: June 09, 2010
Revised: December 28, 2010
Accepted: September 12, 2011.