Research Article  Open Access
Ameel M. R. AlMayah, "Simulation of Enzyme Catalysis in Calcium Alginate Beads", Enzyme Research, vol. 2012, Article ID 459190, 13 pages, 2012. https://doi.org/10.1155/2012/459190
Simulation of Enzyme Catalysis in Calcium Alginate Beads
Abstract
A general mathematical model for a fixed bed immobilized enzyme reactor was developed to simulate the process of diffusion and reaction inside the biocatalyst particle. The modeling and simulation of starch hydrolysis using immobilized αamylase were used as a model for this study. Corn starch hydrolysis was carried out at a constant pH of 5.5 and temperature of . The substrate flow rate was ranging from 0.2 to 5.0 mL/min, substrate initial concentrations 1 to 100 g/L. αamylase was immobilized on to calcium alginate hydrogel beads of 2 mm average diameter. In this work MichaelisMenten kinetics have been considered. The effect of substrate flow rate (i.e., residence time) and initial concentration on intraparticle diffusion have been taken into consideration. The performance of the system is found to be affected by the substrate flow rate and initial concentrations. The reaction is controlled by the reaction rate. The model equation was a nonlinear second order differential equation simulated based on the experimental data for steady state condition. The simulation was achieved numerically using FINITE ELEMENTS in MATLAB software package. The simulated results give satisfactory results for substrate and product concentration profiles within the biocatalyst bead.
1. Introduction
Enzyme immobilization on to supports (or carriers) and their applications as catalysts have grown considerably during the last three decades, and during the last few years have become the most exciting aspects of biotechnology [1–3]. Several methods of enzyme immobilization exist and can be classified into three main categories: carrier binding, cross linking, and entrapment [1]. A number of advantages of enzyme immobilization on to support and several major reasons are the ability to stop the reaction rapidly by removing the enzyme from the reaction solution (or vice versa), products being free of enzyme (especially useful in the food and pharmaceutical industries), reduced effluent disposal problems, suitability for continuous reactor operation, and multiple or respective use of a single batch of enzymes, especially if the enzymes are scarce or expensive, their applicability to continuous processes, and the minimization of pH and substrateinhibition effects. This has an obvious economic impact and allows the utilization of reactors with high enzyme loads [4].
Enzyme entrapment within a gel matrix is one of the enzyme immobilization ways. In this way, the enzyme is surrounded by a semipermeable membrane. Enzyme support of a specific structure permits the contact between the substrate and the biocatalyst in an appropriate way [5, 6].
Enzyme entrapment in calcium alginate beads has been shown to be a relatively simple and safe technique [7–9]. Calcium alginate beads made with 2% (w/v) solution have an average pore diameter of 80 to 100 Å [6, 7, 10]. Starch molecules are very large, often reaching a molecular weight of 80 million Daltons [9, 11]. It is expected that starch hydrolysis reaction occurs more effectively if enzyme is bound to the surface. Several methodshave been compared and reported for the immobilization of αamylase on different supports which provide useful information on the efficiency of the hydrolysis of starch into smaller sugars [12–15].
The catalytic activity of the immobilized enzyme is affected mainly by the limitations of internal and external mass transfer. External mass transfer limitations can be reduced by changing the reactor hydraulic conditions (e.g., the level of agitation), while the internal mass transfer limitations are severe and much more difficult to solve [4]. Enzyme biochemical properties and reaction type and kinetics as well as support chemical and mechanical properties all affect the internal mass transfer [4].
To overcome these limitations, small carrier particle size is proposed to be used. Use of small particle size complicates the reactor operation due to the increasing the pressure drop (i.e., in case of packed bed reactor) or increasing catalyst washout (i.e., fluidized bed reactors). However the idea cannot be generalized for every enzyme. For example, in the case of enzymes that are inhibited by substrate concentration, the operation can be improved by the reduction of the local substrate concentration.
Many mathematical models have been developed by a large number of researchers supported by experimental data, for different types of reactors and mode of operation containing immobilized enzymes [16–19]. Simulation of these models can contribute to improve the understanding of the immobilized systems as well as in the prediction of substrate consumption and product formation rates.
Some of theses models consider the mass transfer of substrate from the bulk to the active sites where the enzyme is immobilized inside the carrier where the reaction takes place. Fick's law is used to model the mass transfer inside the biocatalyst particle [17, 20]. Enzyme catalysis is nonlinear which makes mathematical models more complex, even for the simple MichaelisMenten kinetics model [1, 20–22].
The present study aimed to simulate the process of diffusion and reaction inside the biocatalyst particle. A set of differential equations obtained for spherical immobilized biocatalyst particles allowings the determination of the concentration profile of substrate and product within the biocatalyst in terms of particle geometry (radius) and the concentration of substrate in the bulk liquid phase.
The hydrolysis reaction of corn starch catalyzed by immobilized αamylase within hydrogel matrix was chosen and a mathematical model including the effects of diffusion rate and reaction rate parameters in steady state conditions with scaling analysis were used in this study.
2. Kinetic Parameters
Starch is a major component in many crops such as wheat, maize, tapioca, corn, and potato [2]. Starch consists of a mixture of amylose (15–30%) and amylopectin (70–85%). Starch is hydrolyzed to low molecular weight hydrocarbons by the action of either acids or enzymes [23].
αamylase enzyme hydrolyzes the internal α1,4glycosidic links that exist in amylase and amylopectin to produce low molecular weight products [23–25]. Starch hydrolysis is influenced by several factors such as crystal structure, particle size, amylase and amylopectin content, and the presence of enzyme inhibitors [2, 24]. The kinetics of starch hydrolysis follows the MichaelisMenten model as confirmed by results of many authors [2, 24] and as explained by (1). where is the enzyme concentration, is the substrate concentration, is the rate constant for the forward reaction between enzyme and substrate, is the rate constant for the backward reaction, and denote the rate constant for the ES complex dissociation.
The rate of reactant consumption or product formation can be expressed as [26–28]: where denoted the initial enzyme concentration.
For plug flow reactor, (2) can be rewritten as follows [28–30]:
where represent the rate of disappearance of substrate.
Rearranging (5) and integration for the boundary conditions [][] at and [][] at time , yields the integrated form of the MichaelisMenten model:
For fixed bed reactor, the reaction time or the residence time the reactant spends in the reactor is equal to , where , is the reactor volume and is the substrate volumetric flow rate. Reactor voidage can be expressed as (where denoted the volume of the immobilized enzyme beads), then the residence time can be expressed as .
Rearranging (7) yields Equation (8) can be further arranged in terms of conversion as: or where .
A plot of versus gives a straight line with slope equal to and intercept of .
3. Mathematical Model
In the present study, the enzyme was immobilized on alginate beads; the intraparticle mass transfer resistance can be affecting the rate of enzyme reaction. In order to derive an equation that shows the effect of masstransfer on the effectiveness of immobilized enzyme. The following assumptions were made in developing and solving the mathematical model.(1)The reaction takes place at a constant temperature.(2)Enzyme concentration is constant and uniformly distributed within the beads.(3)Enzyme activity is uniform within the beads.(4)The pressure drop across the reactor and other mechanical effects are negligible.(5)All physical and transport properties are constant except rate constant.(6)Steady state conditions.(7)The substrate concentration is constant within the bulk liquid.(8)Enzyme kinetics is well described by the MichaelisMenten model.(9)External mass transfer limitations are negligible. (10)The immobilized enzyme bead is spherical axissymmetric geometry.(11)Constant diffusivities of substrate and products within the beads.(12)Mass transfer through the immobilized enzyme occurs via molecular diffusion.
In this model cylindrical coordinates were used to locate the domain where diffusion and enzyme activity take place. Based on the above assumptions, the following unsteady state diffusionreaction partial differential equations of mass balances for the product and substrate concentrations within the immobilized enzyme bead (and as shown in Figure 1) can be written as: where and are denoted to the substrate and product concentration, and are the cylindrical coordinates, and and denoted the effective diffusivity of substrate and product.
The axis is the axis of symmetry of the enzyme bead that measures the distance from the planar bead. The variable measures the radial distance from the axis.
The initial and boundary conditions to be assumed for the present problem are as follows:
In order to simplify the solution of (12) the following dimensionless parameters were introduced here as [29]: Equation (12) can be rewritten in the normalized form with the dimensionless parameters as:
The dimensionless parameter known as Thiele modulus is very important because it relates the reaction rate with diffusion rate. In order to simplify the calculations by relating the product concentration to substrate concentration (15) and (16) can be reduced to one linear partial differential as follows:
Let , then, for a steady state conditions, the change in substrate concentration, , and are equal to zero, then (15) will become:
To simulate the reaction inside the biocatalyst, (17) and are used in the present study with the following initial boundary conditions:
4. Experimental Work
4.1. Materials
Sodium alginate salt (type Spectrum Chemical Mfg. Corp.), calcium chloride (CaCl_{2}·2H_{2}O, of BDH type), fungal amylase powder (EC 3.2.1.1; 1,4αDglucose glucanohydrolase, of MP Biomedicals type), corn starch (BDH), and HCl (BDH) were used in this study. All other chemicals were analytical grade reagents.
4.2. Buffer Solution
8 g of KH_{2}PO_{4} (monobasic phosphate) and 0.2 g of K_{2}HPO_{4} (dibasic phosphate) were dissolved in 0.5 liter of water to make a pH buffer solution in order to keep the pH value near to 5.5 during the enzymatic hydrolysis of the starch.
4.3. Iodine Solution
100 mL of iodine solution (0.5% KI and 0.15% I_{2}) was diluted with distilled water till the final volume reached 300 mL. It was used to make a complex compound by reaction with residual starch in the collected samples. The final complex compound has a deep blue color in order to measure the absorbency of this compound.
4.4. Preparation of Substrate Solution
40 g of corn starch powder was mixed with 50 mL of water in a beaker. The slurry was added to 900 mL of warm water in a large beaker. During this period, the slurry was mixed well using a magnetic stirrer and then cooled to room temperature to get the final gelatinized starch solution. Additional volume of water was added in order to bring the total volume to 1 liter.
A few drops of iodine were added to the solution, and the solution color changed to blue which indicates the presence of starch in the solution.
An equal volume of the above buffer and starch solutions were mixed. The resulting solution contains 20 g/L of starch in a buffered environment.
4.5. Enzyme Immobilization
3% (w/v) solution of sodium alginate was prepared by dissolving 15 g of it in 500 mL of water. During the preparation, sodium alginate powder was added slowly (to prevent clomping) to the beaker of water while stirring on a magnetic stirrer. After that, 0.01 g of αamylase enzyme powder was mixed with 150 mL of sodium alginate solution. The final mixture was dropped using a syringe into 500 mL solution of 0.2 M CaCl_{2}. Finally, the beads were left for 2 hours in the calcium solution to get the final hardened form of 2 mm average bead diameter. The final beads were removed from the calcium solution and washed five times with distilled water to remove the excess calcium chloride.
4.6. Immobilized Enzyme Assay
20 mL of 3% starch solution and 2 g of calcium alginate beads were assayed by using Bernfeld method [31]; also this method is recommended by SigmaAldrich Company [32]. Details of this test are stated in this reference. One unit was defined as the amount of amylase that produced 1 mmole of reducing sugar under assay condition per gram of bead.
4.7. Fixed Bed Reactor
A 2 cm 20 cm glass column was used in the present study as a fixed bed reactor. The reactor is surrounded with a jacket of water in order to keep the reaction in an isothermal condition. Two layers of cotton 2 mm thick were placed at the two ends of the column in order to supports the beads and distribute the substrate solution uniformly. The substrate flow rate was adjusted using a dosing pump (B. Braun Melsungen AG, model: 870602) and in a range of 0.2–5.0 mL/min. Samples of product effluents were collected at specific time intervals when the conditions were in steady state. In this study the hydrolysis of corn starch was carried out at constant temperature of 50°C, atmospheric pressure, and at constant pH of 5.5. Figure 2 shows the experimental fixed bed reactor setup.
4.8. Packed Bed Void Fraction
The void fraction () was determined experimentally using liquid impregnation method which can be illustrated as follows.
A 60 cm^{3} reactor volume was used in this work. First, the reactor was filled with the immobilized enzyme beads, and then distilled water was poured in it till it covers all the beads. After that the reactor was drained to measure the volume of water which is equal to the volume of the void. The void fraction can be calculated using the following equation:
4.9. Enzyme Beads
Enzyme beads were placed in the reactor between two layers of cotton. The beads have an average diameter of 2 mm and density of 1.2 g/cm^{3}. The packed bed has an average void fraction of 0.42.
4.10. Stopping Solution
0.1 N HCl solution was used as a stopping solution to stop the hydrolysis reaction of starch for the collected samples in order to analyze the starch content.
4.11. Reactor Operational Efficiency
The entrapment efficiency was used to express the reactor operational efficiency. The entrapment efficiency was calculated at the same hydrolysis conditions. A buffer solution (pH 5.5) was used to wash the beads and the effect of feed flow rate (i.e., residence time) on the entrapment efficiency was determined in the range of 0.2–5 mL/min. The efficiency of entrapment was evaluated during 20 h continuous operation in a fixed bed reactor by collecting the samples from the reactor outlet solution at the end of 1, 2, 4, 6, 8, 10, 15, and 20 h, and calculating the unentrapped enzyme using the same procedure of enzyme assay. The entrapment efficiency of the enzyme was calculated using the following equation: where is the enzyme initial amount and is the enzyme amount in the outlet solution.
4.12. Starch and Products Analysis
The collected samples at timed intervals were analysed using Cintra 5 Double Beam UVSpectrophotometer for residual starch content. Sample absorbency was measured at 620 nm.
5. Results and Discussion
5.1. Starch Hydrolysis
The effect of substrate (starch) flow rate on substrate concentration is explained in Figure 3. It can be observed that as the substrate flow rate increases substrate conversion decreases and at different values of substrate initial concentration, because the residence time is inversely proportional to substrate flow rate. Figure 3 shows the relation between residence time and substrate concentration. In this figure it can also be observed that the relation is linear at the first quarter of the time domain, and then the rate of change decreases rapidly. This is due to the fact that at the beginning of the reaction there is not a lot of substrate present near the enzymes, and the rate increases as the substrate increases because this will give the enzymes more substrate to work on; the rate of change continues till at a certain point, the rate of change decreases, because most of the active sites of enzyme within the bead will be saturated with the substrate to act on. After this point, the substrate concentration becomes too much for the enzyme to work on and the rate of change does not increase further so the rate of reaction becomes nearly constant.
5.2. Reactor Performance
Volumetric activity is another important parameter for bioreactor; it allows decreasing reactor volume and reduces production costs. In this section, the performance as stability of the entrapped αamylase enzyme in the beads was evaluated in a continuous mode. Figure 4 shows the relationship between the entrapment efficiency and feed flow rate (i.e., residence time). A low decrease in reactor stability with time of operation can be observed leading to a small change in reactor stability. Less than 5% reduction in the entrapment efficiency (95%) was observed at the end of 20 h continuous operation and at 0.21 mL/min feed flow rate. As the flow rate increases entrapment efficiency decreases and reached 70% at 5.04 mL/min feed flow rate and for the same operation. This can be attributed to the enzymatic leakage into the buffer solution from the alginate beads. As the reactor volumetric activity increases (i.e., low residence time) enzyme leakage increases which in turn decrease entrapment efficiency.
5.3. Kinetic Parameters
Kinetic parameters, MichaelisMenten constant maximum activity , , and were determined at different substrate flow rates (i.e., residence times) and initial concentrations using (4) and (10). These values are estimated from the slope and intercept of the straight lines shown in Figure 5. The values of , , and at different residence times were listed in Table 1. These values were drawn against residence time and the relation between them can be shown in Figures 6, 7, and 8. As the residence time increases (i.e., substrate flow rate decreases), value decreases and it also increases with increasing initial substrate concentration. The (which is named molecular activity, or turnover number of an enzyme, and is the number of substrate molecules converted to product by an enzyme molecule per unit time when the enzyme is fully saturated with substrate) value also decreases with increasing contact time and increases the initial substrate concentration because it is direct proportional to value according to (4). These observations can be attributed to the reason that with increasing residence time there is an increasing amount of substrate supplied to the enzyme and this phenomenon is not expected to continue for long where saturation has to be reached at some sites in the bead.

It can also be observed that value increases with increasing residence time and initial substrate concentration. Thus, the lower the value of , the higher the affinity of enzyme for substrate. The velocity of the enzymecatalyzed reaction is limited by the rate of breakdown of the ES complex. Then with increasing value more products are formed.
5.4. Simulation of Hydrolysis Reaction within the Hydrogel Bead
The solution for (17) was achieved using FINITE ELEMENTS in MATLAB V. 2008A software package.
Figure 9 represent the algorithm for the computer simulation which is used to simulate substrate and product concentrations.
The effective diffusivity of substrate was calculated at the proposed conditions of the present work and according to the references [33, 34]. The details of the method and equations are illustrated in these references. This value is equal to cm^{2}/min.
The simulated dimensionless substrate and product concentration profiles are shown in Figures 10, 11, 12, 13, 14, 15, and 16. The colorscale map was used to study the substrate and product concentration profiles within the bead. It was assumed that two phases exist: solid bead phase and bulk liquid phase. At low substrate residence time (i.e., high substrate flow rate), the substrate drops rapidly only near the interface between the bead phase and the bulk liquid phase, and in this case the substrate reacts in a fast manner and it will never diffuse into the internal part of the bead and this is can be shown very well in Figure 10(a); on the other hand, the product is formed near the interface and diffuses at a very slow rate (as shown in Figure 10(b)). As the residence time increases (as shown in Figure 10(a) and Figure 16(a)) the reaction region increases and the substrate diffuses into the internal part of the bead. On the other hand, the product formed in a very slow manner and diffused at a slow rate so its concentration remains nearly low (as shown in Figure 10(b) to Figure 16(b)). According to assumption IV listed above, the system is at a steady state. Thus the composition and mass must be unchanged; substrate cannot accumulate in the shell.
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
(a) Dimensionless substrate concentration
(b) Dimensionless product concentration
The relation between Thiele modulus () and the residence time can be well implemented by Figure 17. The is a measure of whether the process is reaction rate controlled at low or diffusion rate controlled at high. It can be observed that increases with decreasing residence time. As increases, this means that mass transfer is much slower than the reaction; it is possible that all substrate entering the particle will be consumed before reaching the center of the bead. In this case, the concentration drops rapidly within the solid as illustrated by Figure 18. The active sites occupied by the immobilized enzyme near the center are starved of substrate and the core of the bead becomes inactive. At low value the concentration of substrate within the bead is naturally higher or lower than in the liquid phase. It can be stated from Figure 17 that starch hydrolysis at the proposed conditions is reaction rate controlled. Its reaction rate controllability decreases as the value of decreases.
6. Conclusions
In the present study, the performance of starch hydrolysis using αamylase immobilized enzyme in a fixed bed reactor at steady state conditions was presented. The numerical solution of mass balances differential equations expressed in terms of a set of dimensionless parameters with the assumed initial and boundary conditions gave the following conclusions.(i)The system performance is strongly affected by the substrate flow rate (i.e., residence time) and initial concentration.(ii)The results of reactor performance experiments demonstrated that the immobilized amylase retained 95% of activity after 20 h of continuous operation at 0.21 ml/min feed flow rate while at 5.04 mL/min the retained activity is 70%. Hence, it could be concluded that the immobilization of amylase in calcium alginate is very useful for continuous starch hydrolysis. (iii)The determined Thiele modulus () values indicated that the reaction rate was controlled by reaction rate within the calcium alginate hydrogel beads. (iv)A decrease in values determines an improvement of substrate conversion. (v)The simulation, which has been performed with experimental data, gave satisfactory results for the substrate and product concentration profiles.(vi)The simulation results comprise a valuable tool for immobilized enzyme reactor design by providing a quantitative relation of enzyme performance with operational variables like substrate flow rate, initial concentration, conversion, and particle size.
Nomenclature
Dimensionless concentration (—)  
:  Effective diffusivity (/min) 
:  Enzyme concentration (g/L) 
:  Initial enzyme concentration (g/L) 
:  Enzymesubstrate complex (—) 
:  Rate constant () 
:  Michaelis constant (g/L) 
:  Product concentration (g/L) 
:  Substrate flowrate, (/min) 
:  Radius of immobilized enzyme bead (cm) 
:  Radius of immobilized enzyme bead (cm) 
:  Dimensionless radius (—) 
:  Substrate concentration (g/L) 
:  time (min) 
:  Rate of reactant decomposition () 
:  Rate of reaction in (2) () 
:  Volume () 
:  Dimensionless concentration (—) 
:  Conversion (%) 
:  Diminsionless hight (—). 
:  Entrapped enzyme amount in the beads 
:  Unentrapped enzyme amount in the reactor outlet solution 
1:  Forward reaction between enzyme and substrate 
1:  Backward reaction 
2:  ES complex dissociation 
:  Bead in equation 
:  Bulk solution 
:  Maximum value 
0:  Initial value (—) 
:  Product 
:  Reactor 
:  Substrate 
:  Concentration at time . 
:  Thiel’s modulus, (—) 
:  Dimensionless time (—) 
:  Dimensionless ratio (—) 
:  Voidage (—). 
References
 P. Valencia, L. Wilson, C. Aguirre, and A. Illanes, “Evaluation of the incidence of diffusional restrictions on the enzymatic reactions of hydrolysis of penicillin G and synthesis of cephalexin,” Enzyme and Microbial Technology, vol. 47, no. 6, pp. 268–276, 2010. View at: Publisher Site  Google Scholar
 D. Djabali, N. Belhaneche, B. Nadjemi, V. Dulong, and L. Picton, “Relationship between potato starch isolation methods and kinetic parameters of hydrolysis by free and immobilised αamylase on alginate (from Laminaria digitata algae),” Journal of Food Composition and Analysis, vol. 22, no. 6, pp. 563–570, 2009. View at: Publisher Site  Google Scholar
 D. Gangadharan, K. M. Nampoothiri, S. Sivaramakrishnan, and A. Pandey, “Immobilized bacterial αamylase for effective hydrolysis of raw and soluble starch,” Food Research International, vol. 42, no. 4, pp. 436–442, 2009. View at: Publisher Site  Google Scholar
 A. Illanes, “Immobilized biocatalysts,” in Comprehensive Biotechnology, vol. 1, pp. 25–39, Elsevier, 2nd edition, 2011. View at: Google Scholar
 P. M. Doran, “Heterogeneous reactions,” in Bioprocess Engineering Principles, pp. 297–332, Academic Press, 1995. View at: Google Scholar
 S. Gülay, Immobilization of thermophilic recombinant esterase enzyme by entrapment in coated Caalginate beads [M.S. thesis], The İzmir Institute of Technology, 2007.
 Z. Konsoula and M. LiakopoulouKyriakides, “Thermostable αamylase production by Bacillus subtilis entrapped in calcium alginate gel capsules,” Enzyme and Microbial Technology, vol. 39, no. 4, pp. 690–696, 2006. View at: Publisher Site  Google Scholar
 S. Talekar and S. Chavare, “Optimization of immobilization of αamylase in alginate gel and its comparative biochemical studies with free αamylase,” Recent Research in Science and Technology, vol. 4, no. 2, pp. 01–05, 2012. View at: Google Scholar
 L. N. Meyer, Effect of immobilization method on activity of alphaamylase [thesis], The Ohio State University, 2007.
 E. Demirkan, S. Dincbas, N. Sevinc, and F. Ertan, “Immobilization of B. amyloliquefaciensαamylase and comparison of some of its enzymatic properties with the free form,” Romanian Biotechnological Letters, vol. 16, no. 6, pp. 6690–6701, 2011. View at: Google Scholar
 A. Riaz, S. A. U. Qader, A. Anwar, and S. Iqbal, “Immobilization of a thermostable Aamylase on calcium alginate beads from Bacillus subtilis KIBGEHAR,” Australian Journal of Basic and Applied Sciences, vol. 3, no. 3, pp. 2883–2887, 2009. View at: Google Scholar
 G. Dey, B. Singh, and R. Banerjee, “Immobilization of αamylase produced by Bacillus circulans GRS 313,” Brazilian Archives of Biology and Technology, vol. 46, no. 2, pp. 167–176, 2003. View at: Google Scholar
 S. Dincbas and E. Demirkan, “Comparison of hydrolysis abilities onto soluble and commercial raw starches of immobilized and free B. amyloliquefaciensαamylase,” Journal of Biological & Environmental Sciences, vol. 4, no. 11, pp. 87–95, 2010. View at: Google Scholar
 A. Anwar, S. A. U. Qader, A. Raiz, S. Iqbal, and A. Azhar, “Calcium Alginate: a support material for immobilization of proteases from newly isolated strain of Bacillus subtilis KIBGEHAS,” World Applied Sciences Journal, vol. 7, no. 10, pp. 1281–1286, 2009. View at: Google Scholar
 D. Park, S. Haam, K. Jang, I. S. Ahn, and W. S. Kim, “Immobilization of starchconverting enzymes on surfacemodified carriers using single and coimmobilized systems: properties and application to starch hydrolysis,” Process Biochemistry, vol. 40, no. 1, pp. 53–61, 2005. View at: Publisher Site  Google Scholar
 A. Horta, J. R. Álvarez, and S. Luque, “Analysis of the transient response of a CSTR containing immobilized enzyme particles. Part I. Model development and analysis of the influence of operating conditions and process parameters,” Biochemical Engineering Journal, vol. 33, no. 1, pp. 72–87, 2007. View at: Publisher Site  Google Scholar
 D. Jeison, G. Ruiz, F. Acevedo, and A. Illanes, “Simulation of the effect of intrinsic reaction kinetics and particle size on the behaviour of immobilised enzymes under internal diffusional restrictions and steady state operation,” Process Biochemistry, vol. 39, no. 3, pp. 393–399, 2003. View at: Publisher Site  Google Scholar
 M. Dadvar and M. Sahimi, “Pore network model of deactivation of immobilized glucose isomerase in packedbed reactors II: Threedimensional simulation at the particle level,” Chemical Engineering Science, vol. 57, no. 6, pp. 939–952, 2002. View at: Publisher Site  Google Scholar
 M. Dadvar, M. Sohrabi, and M. Sahimi, “Pore network model of deactivation of immobolized glucose isomerase in packedbed reactors I: Twodimensional simulations at the particle level,” Chemical Engineering Science, vol. 56, no. 8, pp. 2803–2819, 2001. View at: Publisher Site  Google Scholar
 S. Loghambal and L. Rajendran, “Mathematical modeling in amperometric oxidase enzymemembrane electrodes,” Journal of Membrane Science, vol. 373, no. 12, pp. 20–28, 2011. View at: Publisher Site  Google Scholar
 A. E. AlMuftah and I. M. AbuReesh, “Effects of internal mass transfer and product inhibition on a simulated immobilized enzymecatalyzed reactor for lactose hydrolysis,” Biochemical Engineering Journal, vol. 23, no. 2, pp. 139–153, 2005. View at: Publisher Site  Google Scholar
 S. Bhatia, W. S. Long, and A. H. Kamaruddin, “Enzymatic membrane reactor for the kinetic resolution of racemic ibuprofen ester: modeling and experimental studies,” Chemical Engineering Science, vol. 59, no. 2223, pp. 5061–5068, 2004. View at: Publisher Site  Google Scholar
 H. Zhang and Z. Jin, “Preparation of products rich in resistant starch from maize starch by an enzymatic method,” Carbohydrate Polymers, vol. 86, no. 4, pp. 1610–1614, 2011. View at: Publisher Site  Google Scholar
 V. Varatharajan, R. Hoover, J. Li et al., “Impact of structural changes due to heatmoisture treatment at different temperatures on the susceptibility of normal and waxy potato starches towards hydrolysis by porcine pancreatic alpha amylase,” Food Research International, vol. 44, no. 9, pp. 2594–2606, 2011. View at: Publisher Site  Google Scholar
 A. Mukherjee, A. K. Ghosh, and S. Sengupta, “Purification and characterization of a thiol amylase over produced by a noncereal nonleguminous plant, Tinospora cordifolia,” Carbohydrate Research, vol. 345, no. 18, pp. 2731–2735, 2010. View at: Publisher Site  Google Scholar
 R. S. S. Kumar, K. S. Vishwanath, S. A. Singh, and A. G. A. Rao, “Entrapment of αamylase in alginate beads: single step protocol for purification and thermal stabilization,” Process Biochemistry, vol. 41, no. 11, pp. 2282–2288, 2006. View at: Publisher Site  Google Scholar
 Z. Konsoula and M. LiakopoulouKyriakides, “Starch hydrolysis by the action of an entrapped in alginate capsules αamylase from Bacillus subtilis,” Process Biochemistry, vol. 41, no. 2, pp. 343–349, 2006. View at: Publisher Site  Google Scholar
 A. G. Marangoni, Enzyme Kinetics a Modern Approach, John Wiley & Sons, 2003.
 J. M. Lee, Biochemical Engineering, PrenticeHall, 2011.
 M. Y. Arica, V. Hasirci, and N. G. Alaeddinoǧlu, “Covalent immobilization of αamylase onto pHEMA microspheres: preparation and application to fixed bed reactor,” Biomaterials, vol. 16, no. 10, pp. 761–768, 1995. View at: Publisher Site  Google Scholar
 P. Bernfeld, S. P. Colowick, and N. O. Kaplan, “Amylases α and β,” in Methods in Enzymology, vol. 1, pp. 149–158, Academic Press, New York, NY, USA, 1955. View at: Google Scholar
 http://www.sigmaaldrich.com/sigmaaldrich/technicaldocuments/protocols/biology/enzymaticassayofaamylase.html.
 P. Grunwald, K. Hansen, and W. Gunber, “The determination of effective diffusion coefficients in a polysaccharide matrix used for the immobilization of biocatalysts,” Solid State Ionics, vol. 101–103, part 2, pp. 863–867, 1997. View at: Google Scholar
 J. Coulson, M. J. F. Richardson, J. R. Backhurst, and J. H. Marker, Chemical Engineering, vol. 1, ButterworthHeinemann, 6th edition, 1999.
Copyright
Copyright © 2012 Ameel M. R. AlMayah. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.