Abstract
Due to great commercial application of protease, it is necessary to study kinetic characterization of this enzyme in order to improve design of enzymatic reactors. In this study, mathematical modeling of protease enzyme production kinetics which is derived from Bacillus licheniformis BBRC 100053 was studied (at 37°C, pH 10 after 73 h in stationary phase, and 150 rpm). The aim of the present paper was to determine the best kinetic model and kinetic parameters for production of protease and calculating (inhibition constant) of different inhibitors to find the most effective one. The kinetic parameters (Michaelis-Menten constant) and (maximum rate) were calculated 0.626 mM and 0.0523 mM/min. According to the experimental results, using DFP (diisopropyl fluorophosphate) and PMSF (phenylmethanesulfonyl fluoride) as inhibitors almost 50% of the enzyme activity could be inhibited when their concentrations were 0.525 and 0.541 mM, respectively. for DFP and PMSF were 0.46 and 0.56 mM, respectively. Kinetic analysis showed that the Lineweaver-Burk model was the best fitting model for protease production kinetics DFP was more effective than PMSF and both of them should be covered in the group of noncompetitive inhibitors.
1. Introduction
Microbial proteases dominate the worldwide enzyme market, accounting for a two-thirds share of the detergent industry [1]. Microbial protease is one of the most important groups of industrial enzymes that caters to the requirement of nearly 60% of the worldwide enzyme market [2]. Of these, the alkaline proteases are of interest from a biotechnological perspective and find numerous applications in the food industry and meat tenderization process, peptide synthesis, infant formula preparations, baking, and brewing. Furthermore, they are used in pharmaceutical and medical diagnosis, in the detergent industry as additives, and in textile industry in the process of dehairing and leather processing [3, 4].
Currently, many economically important industrial enzymes are produced by cultivation of bacteria, such as Bacillus sp. [5, 6]. Over 300 tons of enzymes, mainly proteases, are being annually produced from Bacillus sp.
Proteases with high activity and stability in high alkaline range and high temperatures are interesting for bioengineering and biotechnological applications [7].
Kinetic study of enzyme constitutes the information pertaining to rates of activation and inactivation of enzymes and actually gives the rate at which a process occurs. The kinetic parameters controlling the rates of enzyme catalyzed reaction are incubation time, initial substrate concentration, initial enzyme concentration, pH, and reaction temperature [8].
Regarding previous studies, the optimized conditions for protease production, enzyme characterization, and enzyme immobilization were determined [9, 10]. Although there have been a number of studies on protease production by Bacillus species, little information on kinetic analysis of the protease production process is available in the literature. This research followed the study kinetic model for the production of protease enzyme based on experimental data. In this paper, the experiments were conducted to study the effects of various inhibitors on enzymatic activity of protease by determination of kinetic parameters and inhibition constants. Kinetic parameters of protease enzyme have been determined based on some Bacillus genus in which Bacillus licheniformis is less mentioned.
2. Materials and Methods
2.1. Bacterial Strain
The organism used was Bacillus licheniformis BBRC 100053 obtained from Biochemical and Bioenvironmental Research Center, a local culture collection in Sharif University Of Technology, Iran. The culture was maintained on nutrient agar medium at 30°C for 2 days and stored at 4°C [9]. Two inhibitors, DFP (diisopropyl fluorophosphate) and PMSF (phenylmethanesulfonyl fluoride), were obtained from Merck Co., Germany. All other reagents used were also of analytical reagent (AR) grade.
2.2. Inoculum Preparation and Protease Production
Inocula were prepared by adding a loop full of pure culture into 25 mL of sterile Luria-Bertani (LB) broth medium (composed of (g/L) peptone 10, yeast extract 5, and NaCl 5) [9]. Broth cultures were diluted with the appropriate broth to obtain 107 cfu/mL as estimated by absorbance at 600 nm. Target absorbance to obtain these populations was 0.2 for Gram-positive bacilli, based upon preliminary experiments [9–11].
About 8–10% inoculums from this culture were added to 100 mL Erlenmeyer flasks containing 20 mL production media containing lactose 1%, yeast extract 0.5%, peptone 0.5%, KH2PO4 0.1%, MgSO4·7H2O 0.02%, and pH 8 [9]. Media were autoclaved at 121°C for 20 min in 15 psi [3].
After incubation at 37°C for 73 h under shaking (150 rpm) the cultures were centrifuged and the supernatants were used for estimation of proteolytic activity. The growth of the microorganism was monitored by measuring absorbance at 660 nm. All experiments were carried out in duplicate.
2.3. Assay for Proteolytic Activity
Protease activity was determined by Anson’s modified method using casein as substrate.
0.5 mL of enzyme solution was added to 4.5 mL of substrate solution (1% w/v, casein with 50 mM Tris-HCl buffer, and pH 8) and incubated at 30°C for 10 min. The reaction was stopped by adding 5 mL of 5% TCA (trichloroacetic acid) mixture (5% TCA, 9% Na-acetate, and 9% acetic acid) followed by 30 min holding at room temperature and centrifugation at 8000 rpm [12].
The precipitates were removed by filtration through Whatman-42 filter paper and absorbance of the filtrate was measured at 280 nm. One unit of protease activity was defined as the amount of enzyme liberating 1 μgr of tyrosine per min under assay conditions. Enzyme units were measured using tyrosine solutions (0–100 mg) as standard [13].
2.4. Determination of Kinetic Parameters (, )
Due to great commercial application of protease, it is necessary to study kinetic characterization of this enzyme in order to improve design of enzymatic reactors [14].
Enzymes are natural catalysts that speed up the chemical reactions. However, the speed of any fastidious reaction being catalysed by a particular enzyme can only reach a certain maximum value. This rate is known as while can be defined as the concentration of substrate at which half of the maximal velocity was obtained [15].
Michaelis-Menten equation:
Taking the reciprocal gives
In order to investigate and , specific concentration of enzyme with different concentration of casein as substrate (0.25%, 0.45%, 0.65%, 0.8%, and 1%) (w/v), was prepared for identifying approximate values. The enzyme (0.5 mL) was added to 4.5 mL casein (1% w/v in 50 mM Tris-HCl buffer, pH 8) and the reaction mixture was incubated at 45°C for 10 min before the addition of 5 mL of TCA mixture (1.8% trichloroacetic acid, 1.8% sodium acetate, and 1.98% acetic acid). In different times (1, 2, 3, 4, and 5 minutes) the absorbance of the filtrate was measured at 280 nm. The slope of the plot of absorbance versus time was 0.1 to 0.4 using the best concentrations. The slope of the plot of absorbance (280 nm) versus time could determine ().
2.5. Effect of Inhibitors on Protease Activity
Different inhibitor concentrations (volume ratio of inhibitor to enzyme) were used and increased up to inhibit more than 50% of the enzyme activity. For PMSF 0.025, 0.0075, 0.002, and 0.001 and for DFP 0.01, 0.005, 0.001, and 0.0004 volume ratios were selected and added at the beginning of exam. The reaction was started by adding 50 μL of enzyme solution to 450 μL of various substrates solutions. The reaction was stopped, by adding 50 μL of 5% TCA to this mixture, followed by filtration and measuring the absorbance at 280 nm. The slope of the plot of versus time determined for different volume ratio (or inhibitor concentration) using various substrate concentrations (0.75%, 1%, 1.25%, and 1.5%).
The plot of versus identified the kind of group in which the inhibitors should be covered [16].
3. Results and Discussions
3.1. Kinetic Constant Calculation of Protease Enzyme Production
According to different for various initial substrate concentrations, the Lineweaver-Burk model was selected as the best kinetic model and and were calculated by the plot of versus . The obtained final ranges of initial substrate concentration were examined to identify amounts of (Table 1). Through the slope and the interception of the plot of versus , the exact values of and were determined 0.626 mM and 0.0523 mM/min, respectively (Figure 1).
3.2. Inhibition Kinetics and Constant
The absorbance of the samples with different inhibitor concentrations and in the presence of several substrate concentrations was measured (280 nm) for DFP and PMSF, as shown in Tables 2 and 3.
Similar experiments were studied in the absence of inhibitors with respect to control (Table 4).
Afterwards, different values were calculated from the slope of the plot of absorbance versus time. According to the plots of versus as shown in Figures 2 and 3 and various types of inhibitors, apparent kinetic parameters ( and ) were calculated (based on Lineweaver-Burk equation) for different inhibitor concentrations. Results showed that both inhibitors (PMSF and DFP) should be classified as noncompetitive inhibitors:
Secondary plot was used to determine inhibition constant (). The slop of the plot and the intersection of the line with X axis showed an increase in ( and were described as the maximum velocity in the presence and absence of inhibitors). As indicated in Figures 4 and 5, obtained 0.56 and 0.46 mM for PMSF and DFP, respectively.
3.3. The Plot of IC50
The specific concentration of inhibitors which was able to inhibit almost 50% of the enzyme activity was identified through Figures 6 and 7. In addition via the plot of enzyme activity () versus inhibitor concentration ([]) and were identified as initial velocity in the presence and absence of inhibitor. The inhibitor concentration met the required value when was 0.5.
50% of the enzyme activity could be inhibited in the presence of 0.525 and 0.541 mM of DFP and PMSF, respectively (Figures 6 and 7). The identical enzyme activity should be intended for all steps.
4. Conclusions
In this paper, mathematical modeling of protease enzyme production kinetics derived from B. licheniformis, under optimum condition (37°C, pH 10), was studied.
Based on the results, the Lineweaver-Burk model was the best fitting model for protease production kinetics The kinetic parameters , were 0.626 mM and 0.0523 mM/min, respectively.
The results showed that DFP and PMSF were noncompetitive inhibitors. In order to reduce the enzyme activity up to 50% of its initial activity, the required concentrations of DFP and PMSF were 0.525 and 0.541 mM. In addition inhibition constant () is obtained 0.46 and 0.56 mM for DFP and PMSF, respectively.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.