Acessibilidade / Reportar erro

Modeling gallic acid production rate by empirical and statistical analysis

Abstracts

For predicting the rate of enzymatic reaction empirical correlation based on the experimental results obtained under various operating conditions have been developed. Models represent both the activation as well as deactivation conditions of enzymatic hydrolysis and the results have been analyzed by analysis of variance (ANOVA). The tannase activity was found maximum at incubation time 5 min, reaction temperature 40ºC, pH 4.0, initial enzyme concentration 0.12 v/v, initial substrate concentration 0.42 mg/ml, ionic strength 0.2 M and under these optimal conditions, the maximum rate of gallic acid production was 33.49 mumoles/ml/min.

Gallic acid; modeling; tannase; Caesalpinia digyna


Para predizer a taxa das reações enzimaticas uma correlação empírica baseada nos resultados experimentais foi desenvolvida. Os modelos representam a ativação e a desativativação da hydrolise enzimatica. Os resultados foram avaliados pela análise de variança (ANOVA). A atividade máxima da tannase foi obtida após 5 minutos de incubação, temperatura 40ºC, pH 4,0, concentração inicial da enzima de 0,12 v/v, concentração inicial do substrato 0,42 mg/ml, força iônica 0,2 M. Sob essas condições a taxa máxima de produção ácido galico foi de 33,49 µmoles/ml/min.


Modeling Gallic Acid Production Rate by Empirical and Statistical Analysis

Bratati Kar1, Rintu Banerjee1* * Author for correspondence and Bimal Chandra Bhattacharyya2

1Microbial Biotechnology & Downstream Processing Laboratory,Agricultural & Food Engineering Department, Indian Institute of Technology, Kharagpur - 721 302, India;

2Visiting Professor, College of Leather Technology, Calcutta, India

ABSTRACT

For predicting the rate of enzymatic reaction empirical correlation based on the experimental results obtained under various operating conditions have been developed. Models represent both the activation as well as deactivation conditions of enzymatic hydrolysis and the results have been analyzed by analysis of variance (ANOVA). The tannase activity was found maximum at incubation time 5 min, reaction temperature 40oC, pH 4.0, initial enzyme concentration 0.12 v/v, initial substrate concentration 0.42 mg/ml, ionic strength 0.2 M and under these optimal conditions, the maximum rate of gallic acid production was 33.49 mmoles/ml/min.

Key words: Gallic acid, modeling, tannase, Caesalpinia digyna

INTRODUCTION

Bioconversion of chemical feed stock is an economically viable option to chemical process especially in view of its reduced scale of pollution and energy consumption. Caesalpinia digyna (Teri pod) cover was used as the raw material for the conversion of tannin present in it, to gallic acid under modified solid-state fermentation, with the aid of the enzyme, tannase (tannin acyl hydrolase, EC 3.1.1.20) produced by the microorganism Rhizopus oryzae (RO IIT RB-13 NRRL 21498), locally isolated from the soil of IIT campus. Gallic acid is a phenolic compound (3,4,5-trihydroxy benzoic acid) and is used mainly in pharmaceutical industries for manufacturing trimethoprim (TMP). Trimethoprim is an antibacterial agent, which is usually given along with sulfonamide, and together it has a broad spectrum of action. As microbial processes are less predictable and often more sensitive to environmental variations than chemical reactions (Namdev et al., 1991; Patnaik, 1995) because of the complexity of microbial physiology and their interactions with transport processes, there have been few attempts of theoretical analysis. Mathematical modeling is the application of mathematics to obtain answers to real world problems. It may be used in a wide range of applications, some of which do not appear initially to be mathematics in nature. Mathematical model together with carefully designed experiments makes it possible to evaluate the behavior of the systems more rapidly than with laboratory experiments only. Models often allow quick and cheap evaluation of alternatives, leading to optimal solutions, which are not otherwise obvious. Moreover, by using models we can better understand or explain natural phenomena and under some conditions we can make predictions in a deterministic or probabilistic sense.

A mathematical model for a microbial process can be expressed using two different mechanisms: structured and unstructured models(Sinclair & Kristiansen, 1987). The building of kinetic model consists of comparing the predicted values of assumed models with experimental data in order to obtain more relevant correlation. Such model enables the bioengineer to design and control microbial process(Moser and Steiner, 1975) by predicting the behaviour of the same.

Knowledge and an understanding of the kinetics of gallic acid production are of great economical importance in view of the fact that it has got immense application in the pharmaceutical industry. Industrial fermenta-tion is gradually moving away from the traditional and largely empirical approach towards simpler and better controlled process (Trilli 1990).

MATERIALS AND METHODS

Raw Material:Teri pod (Caesalpinia digyna) shrub is found in the forests of the Assam, West Bengal, Andhra Pradesh, Madhya Pradesh and Orissa. The Teri pod cover was readily available from Orissa during the processing of the pod for recovery of oil. It was dried at a temperature of 60oC in an oven to remove the moisture. After drying, the pod cover was finely powdered in a grinder and screened to 74-422 mm. Oversized particles were recycled. The powder was stored in a cool and dry place at room temperature.

Media: Modified Czapek-Dox medium (Hadi et al., 1994) was used to prepare induced inoculum (Chen, 1969). In the present study, the medium constituted of sodium nitrate (0.25% (w/v)), potassium dihydrogen orthophosphate (0.1% (w/v)), magnesium sulphate (0.05% (w/v)), potassium chloride (0.05% (w/v)) and Caesalpinia digyna seed cover powder (2.0 % (w/v)).

Micro-organism: A fungal strain, Rhizopus oryzae (RO IIT RB-13 NRRL 21498), which was isolated from the soil of Indian Institute of Technology campus, Kharagpur, was used for the production of tannase. The culture was maintained on 2% (w/w) malt extract agar slant (Banerjee, 1992).

Fermentation: Biosynthesis experiments were carried out in 250-ml Erlenmeyer flasks containing 100 ml of the modified Czapek-Dox medium without the addition of Caesalpinia digyna pod cover powder, to which 1.0 ml of induced inoculum was added and incubated at 37oC for 72 h in an incubator shaker (Sambros and Instruments India Ltd., India).

Study of tannase kinetics: After fermentation the insoluble material was sedimented by centrifugation in a cold centrifuge (Remi C-24, India). The supernatant was used for study of tannase kinetics under various operating conditions following the assay procedure (Kar et al., 1999).

Modelling: An empirical correlation for predicting the optimum rate of gallic acid production vis-à-vis enzymatic reaction influenced by various physicochemical parameters, namely, incubation time (t, min), reaction temperature (T, oC), pH (H), initial enzyme concentration (E, v/v), initial substrate concentration (S, g/ml), ionic strength of the medium (I, M) has been proposed in the following form:

P = Ø (t) a (T) b (H) c (E) d (S) e (I) f .(1)

where P (mmoles/ml/min) is the gallic acid production rate, Ø is the dimensional proportionality constant and a, b, c, d, e, f are the exponential coefficients. The exponential coefficients and the proportionality constant were evaluated from the experimental results using Excel spreadsheet (Microsoft Office 97, Microsoft Corp.). All the experiments were conducted in triplicate and the average values were taken. Finally the analysis of the results was carried out by the analysis of variance (ANOVA).

RESULTS

Results of the enzyme kinetic experiments for the rate of gallic acid production were obtained for the following conditions: incubation time (min), 1, 2, 3, 4, 5, 7, 10, 15, 18, 20, 25 and 30; reaction temperature (oC), 8, 22.5, 30, 37, 40, 45, 50, 55 and 60; pH 3, 3.5, 4, 4.5, 5, 5.5, 6; initial enzyme concentration (v/v), 0.01, 0.02, 0.03, 0.04, 0.08, 0.12, 0.16, 0.20; initial substrate concentration (mg/ml), 0.028, 0.035, 0.07, 0.14, 0.21, 0.28, 0.35, 0.42, 0.49, 0.56, 0.70, 1.4, 2.8; ionic strength (M), 0.025, 0.05, 0.075, 0.1, 0.15, 0.175, 0.2, 0.25. Two trends, activation and deactivation of the rate of enzymatic reaction were observed from the results under the above conditions. Therefore, for each of the above two trends, a correlation between dependent and independent variables was obtained using the Excel spreadsheet. The equations obtained were:

P = 2.9551x10-6 (t) 0.9614 (T) 1.7138 (H) 4.3848

Statistical analysis of model 1 and model 2 were performed using Excel spreadsheet and the results are shown in Tables 1 and 2 respectively.

(E)1.6664 (S)0.7094 (I)1.1922

(for activation) (2)

P = 1.9832 x 108 (t) -1.3451 (T) -2.7373 (P) - 4.8017

(for deactivation) (3)

Model 1 (Figure 1) and model 2 (Figure 2) were developed using equations 2 & 3 respectively, which depict the results of P (predicted) versus P (experimental).

From tables 1 and 2 it can be seen that F-critical values are greater than the F-values, which indicate that model 1 and model 2 show good agreement with the experimental data. The coefficient of determination (R2) and the coefficient of correlation (R) for the rate of gallic acid production in case of model 1 and model 2 were calculated by regression analysis of the experimental data using Excel spreadsheet. The values were as follows: R2 = 0.979 (model 1) and 0.991 (model 2); R = 0.989 (model 1) and 0.995 (model 2).

DISCUSSION

For the development of the model, the system was divided into two parts: model 1 representing the activation trend (rate of gallic acid production from minimum to maximum level) and model 2 representing the deactivation trend (rate of gallic acid production from maximum to minimum level). The model of deactivation trend is the control model for activation trend as the latter system represents only increasing trend.

Further, the models are valid only for the following ranges of the parameters studied:

The objective of the present investigation was to find a suitable correlation for selecting the best combination of microbiological parameters in order to obtain maximum rate of gallic acid production by the enzyme, tannin acyl hydrolase. The coefficient of determination (R2) is 0.979 and 0.991 for model 1 and model 2 respectively. The value of R2 is the measure of the total variation of the observed value of gallic acid production about the mean explained by the fitted model, which is often expressed in percentage (Rosander, 1951). This implies that 97.9% and 99.1% of the total variation in the gallic acid production which have been explained by the fitted models can be attributed to the independent variables, viz., t, T, H, E, S and I.

The values of R2 obtained, indicated a good agreement between the experimental and the predicted values of rate of gallic acid production (model 1 and model 2). The regression coefficients are valid for any of the variables studied, placed in the analyzed ranges of the corresponding variables. The coefficient of correlation (R) is 0.989 and 0.995 for model 1 and model 2 respectively. It explains the correlation between the experimental and predicted values from model 1 and model 2 respectively and was considered satisfactory.

The F-values corresponding to tannin acyl hydrolase was 0.005 and 0.0009 for model 1 and model 2 respectively (Tables 1 and 2). The F-value is the ratio of the mean square due to regression to the mean square due to error (Rosander, 1951). In general, the critical F-value should be several times greater than the computed F-values, the null hypothesis is rejected when P-value is greater than the F-value and infers that the variation accounted for by the models was significantly greater than the unexplained variations (Rosander, 1951). In the present study, the null hypothesis was rejected at 0.94 and 0.97 (for model 1 and 2 respectively) for the enzyme, tannase.

The optimal levels of parameters thus obtained from equations (2) & (3) for maximum tannase hydrolysis were: incubation time, 5 min; reaction temperature, 40oC; pH, 4.0; initial enzyme concentration, 0.12 (v/v); initial substrate concentration, 0.42 mg/ml; ionic strength of medium, 0.2 M at which gallic acid production was 33.49 (mmoles/ml/min).

RESUMO

Para predizer a taxa das reações enzimaticas uma correlação empírica baseada nos resultados experimentais foi desenvolvida. Os modelos representam a ativação e a desativativação da hydrolise enzimatica. Os resultados foram avaliados pela análise de variança (ANOVA). A atividade máxima da tannase foi obtida após 5 minutos de incubação, temperatura 40oC, pH 4,0, concentração inicial da enzima de 0,12 v/v, concentração inicial do substrato 0,42 mg/ml, força iônica 0,2 M. Sob essas condições a taxa máxima de produção ácido galico foi de 33,49 µmoles/ml/min..

Received March 18, 2000;

Revised: April 23, 2000;

Accepted September 09, 2000.

  • Banerjee, R. (1992), Some studies on production, purification and characterization of extra-cellular protease by Rhizopus oryzae (RO, IITKGP). Ph.D. thesis, IIT, Kharagpur, India.
  • Chen, T. (1969), Ph.D. Thesis of Biochemistry. University of California, Berkeley, U.S.A.
  • Hadi, T. A., Banerjee, R and Bhattacharyya, B. C. (1994), Optimization of tannase synthesis by newly isolated Rhizopus oryzae. Bioproc. Engg, 11, 239-243.
  • Kar, B., Banerjee, R. and Bhattacharyya, B. C. (1999), Microbial production of gallic acid by modified solid sate fermentation. J. Ind. Microbiol. Biotechnol, 23, 173-177.
  • Moser, A. and Steiner, W. (1975), The influence of the term kd for endogenous metabolism on the evaluation of Monod kinetics for biotechnological processes. Eur. J. Appl. Microbiol., 1, 281-289.
  • Namdev, P. K., Yegneswaran, P. K., Thompson, B. G. and Gray, M. R. (1991), Experimental simulation of large scale bioreactor environments using Monte-Carlo method. Can. J. Chem. Eng., 69, 513-519.
  • Patnaik, P. R. (1995), Sensitivity of recombinant fermentation with run-away plasmids: a structural analysis of the effect of dilution rate. Chem. Eng. Commun, 131, 125-140.
  • Rosander, A. C. (1951), Elementary Principles of Statistics, D. Van Nostrand Company, Inc., New York, pp. 514-539.
  • Sinclair C. G. and Kristiansen, B. (1987), Fermentation kinetics and modeling. Milton Keynes, Open Univ Press.
  • Trilli A. (1990), Kinetics of secondary metabolites production. In- Microbial growth dynamics, eds. R. K. Poole, M. J. Bazin and C. W. Keevil, Oxford: IRL Press, Oxford Univ Press, pp. 102-126.
  • *
    Author for correspondence
  • Publication Dates

    • Publication in this collection
      25 Oct 2001
    • Date of issue
      2000

    History

    • Accepted
      09 Sept 2000
    • Reviewed
      23 Apr 2000
    • Received
      18 Mar 2000
    Instituto de Tecnologia do Paraná - Tecpar Rua Prof. Algacyr Munhoz Mader, 3775 - CIC, 81350-010 Curitiba PR Brazil, Tel.: +55 41 3316-3052/3054, Fax: +55 41 3346-2872 - Curitiba - PR - Brazil
    E-mail: babt@tecpar.br