Skip to main content
  • Original article
  • Open access
  • Published:

Construction of a biodynamic model for Cry protein production studies

Abstract

Mathematical models have been used from growth kinetic simulation to gen regulatory networks prediction for B. thuringiensis culture. However, this culture is a time dependent dynamic process where cells physiology suffers several changes depending on the changes in the cell environment. Therefore, through its culture, B. thuringiensis presents three phases related with the predominance of three major metabolic pathways: vegetative growth (Embded-Meyerhof-Parnas pathway), transition (γ-aminobutiric cycle) and sporulation (tricarboxylic acid cycle). There is not available a mathematical model that relates the different stages of cultivation with the metabolic pathway active on each one of them. Therefore, in the present study, and based on published data, a biodynamic model was generated to describe the dynamic of the three different phases based on their major metabolic pathways. The biodynamic model is used to study the interrelation between the different culture phases and their relationship with the Cry protein production. The model consists of three interconnected modules where each module represents one culture phase and its principal metabolic pathway. For model validation four new fermentations were done showing that the model constructed describes reasonably well the dynamic of the three phases. The main results of this model imply that poly-β-hydroxybutyrate is crucial for endospore and Cry protein production. According to the yields of dipicolinic acid and Cry from poly-β-hydroxybutyrate, calculated with the model, the endospore and Cry protein production are not just simultaneous and parallel processes they are also competitive processes.

Introduction

Bacillus thuringiensis (B. thuringiensis) is a Gram-positive endospore-forming bacterium isolated from soil that synthesizes a crystalline δ-endotoxin protein, named Cry protein (Ito et al. [2004]). The main application of the Cry protein is the biological control of certain insects of several orders (Ito et al. [2004]). It has been also reported that some non-toxic Cry protein for insects are highly cytotoxic to a wide range of mammalian cells, particularly human cancer cells (Ito et al. [2004]). Cry protein production by B. thuringie nsis has been largely studied. Many different culture media formulations, carbon-nitrogen ratios, operating conditions and production systems have been used in order to improve the cost-productivity relation (Farrera et al. [1998]; Navarro et al. [2006]; Ozcan et al. [2010]; Zhuang et al. [2011]). Also, several mathematical models describing the growth and endospore formation kinetics of B. thuringiensis are available (Kraemer-Schafhalter and Moser [1996]; Liu and Tzeng [2000]; Popovic et al. [2001]; Rivera et al. [1999]; Starzak and Bajpai [1991]).

B. thuringiensis culture is a time dependent dynamic process where the physiology of cells suffers several changes. During the culture the number of cells increase causing changes in the environment and subsequently the cells respond to those changes modifying their activity. This dynamics continues until the death of the cells when active entomopathogenic products (endospores and Cry protein) are obtained. The response of the cell to the environmental changes is a complicated process that involves cell metabolism, cell cycle regulation, signaling and gene regulatory networks and gene expression. As result of the response of the cell to the environmental changes, the culture typically shows three different phases related with the predominance of one of three major metabolic pathways (Figure 1a): vegetative growth-Embded-Meyerhof-Parnas pathway (EMP), transition-γ-aminobutiric cycle (GABA) and sporulation-tricarboxylic acid cycle (TCA) (Anderson [1990]; Rowe [1990]).

Figure 1
figure 1

Schematic representation of (a) B. thuringiensis culture with their key compounds and (b) model modules.

Mathematical models are powerful tools to describe processes or natural systems that are reproducible and predictable, which also simplify reality. Models are designed to focus on certain aspects of the system of study while other aspects are abstracted away. The mathematical modeling serves as aids to biological investigation in a number of ways. Different models have been used to describe the different phases present in B. thuringiensis culture, specifically: the vegetative growth phase (Kraemer-Schafhalter and Moser [1996]; Popovic et al. [2001]; Rivera et al. [1999]) and the sporulation phase (Liu and Tzeng [2000]). These models consider just two phases of the culture and represent only one or another. Starzak and Bajpai et al. ([1991]) model considers the growth and the sporulation phases but it assumes that during sporulation there are no energy requirements so no link between both phases is considered. Thus, there is no mathematical model which links the different culture phases with its metabolic pathway.

In the present study, a biodynamic model was constructed to describe the complete evolution of the B. thuringiensis culture. Also, it was used to study the interrelation between the different culture phases and its relationship with the Cry protein production. The considerations for the construction of the model were: a) each phase of B. thuringiensis culture can be represented by its principal metabolic pathways (EMP for vegetative growth, GABA for transition and TCA for sporulation); and b) each pathway can be follow with its key compound (cell for EMP, PHB for GABA and DPA-Cry protein for TCA).

Materials and methods

Model development

The previously reported kinetics of the key compounds were analyzed to identify the main variables affecting the three phases of B. thuringiensis culture. It is not straightforward to compare the cells/biomass and PHB kinetics data reported in the literature because different conditions and analytic techniques were used. Thus, in order to be able to compare these data a normalization was done (Additional file 1). The normalization is the mathematical process used to reduce the data to absolute values. The normalized cells/biomass and PHB kinetics show the sigmoidal curve typically well represented by the Gompertz model. Also, the two sets of endospores’ kinetics data published (Farrera et al. [1998]; Navarro et al. [2008]) show a sigmoidal curve. Therefore, the Gompertz model was used to simulate and estimate the kinetic parameters of the key compound data.

Model formulation

The assumptions for the biodynamic model construction were: i) PHB affects positively the Cry production; ii) Yield of DPA, Cry protein and PHB from biomass are independent from initial substrate concentration; iii) Embded-Meyerhof-Parnas, γ-aminobutyric acid cycle and tricarboxylic acid cycle are the main pathways for vegetative growth, transition and sporulation phases, respectively; iv) Every metabolic change has a key compound associated and can be represented by: cells, poly-β-hydroxybutyrate (PHB) and DPA-Cry protein.

As show in Figure 1b the biodynamic model has three principal modules interconnected, one for each phase. Module 1 represents the vegetative growth, module 2 represents the transition and module 3 represents the sporulation phase. In general terms, during vegetative growth the substrate is used to generate new cells and carbon dioxide. The CO2 is released due to the respiration process. During transition the substrate is accumulated as PHB. While during sporulation: there is no substrate consumption; PHB is consumed; the vegetative cell is transformed into endospore; and in the same period of time the Cry protein is produced from protein turnover (Anderson et al. [2002]). The PHB is the energy source for endospore and Cry protein formation processes. Thus, module 1 and module 2 are related with module 3.

Module 1

Phase: Vegetative growth; metabolic pathway: Embded-Meyerhof-Parnas pathway; key compound: cells. The vegetative cell evolution requires a two-term equation: growth and dead (Equation 1). The first term is simulated using the Gompertz model while the second term takes into account that endospores (DPA) are formed from biomass. The meanings of the symbols were defined in Table 1.

X= X max e - e - u max t - t c - X max - DPA Y DPA / X
(1)
Table 1 Nomenclature

Module 2

Phase: Transition; metabolic pathway: γ-aminobutyric acid cycle; key compound: Poly-β-hydroxybutyrate. The PHB also requires a two-term equation: production and consumption (Equation 2). In the first term, PHB production was modeled using the Gompertz model. It has been reported that the PHB serves as an endogenous reserve of carbon and energy during the sporulation and it is necessary for the δ-endotoxin synthesis (Benoit et al. [1990]; Kraemer-Schafhalter and Moser [1996]; Navarro et al. [2006]; Rowe and Margaritis [1987]). Therefore, the second term takes into account that the production of DPA and Cry is based on the consumption of PHB.

PHB=PH B max e - e μ maxp t - t c p - P H B max - DPA Y DPA / PHB - Cry Y Cry / P H B
(2)

Module 3

Phase: Sporulation; metabolic pathway: tricarboxylic acid cycle; key compound: Cry protein and DPA. DPA and Cry protein production are assumed as simultaneous and parallel processes both represented by the Gompertz model, Cry protein (Equation 3) and DPA (endospore) (Equation 4):

Cry = C r y max e - e ( - μ maxc t - t c c
(3)
DPA=DP A max e - e ( - μ maxd t - t c d
(4)

The biodynamic model constructed (Equations 1-4) was solved with ModelMaker 3.0.3. The model parameters were estimated by nonlinear minimum sum of squares analysis (i.e. the differences between the predicted and measured values) using the Marquardt algorithm.

The biodynamic model was constructed using the kinetic data sets available in the literature and the model validation was done by the comparison of the predicted values and those obtained in four new experiments.

New data sets for model validation

The biodynamic model validation was made using the data of the key compounds kinetics obtained in four new fermentations (F1, F2, F3, F4). These fermentations were done using the same culture media, carbon:nitrogen relation, operating condition and bioreactor. In order to establish the relationship between the carbon source and the biomass and Cry protein production, different concentration of the substrates in the culture media were used in each run.

Organism

Bacillus thuringiensis var. kurstaki HD-73 (ATCC-35866), which produces a 133.3 kDa Cry1A(c) insecticidal crystal protein, was obtained from the USDA Insect Pathology Research Unit, Brownsville, Texas. The microorganism conservation was made in endospore filter-paper disk and the inoculum for fermentation was prepared in two stages, according to Navarro et al. ([2006]). The culture of the second stage was used as inoculum for the bioreactor experiments.

Fermentation procedures

All the fermentations were done in a 7 L glass reactor. The reactor was fully automated and equipped with two Rushton turbines. The operating and the sterilization conditions were done according to Navarro et al. ([2006]). The end-of-fermentation criterion was set to 85-90% of endospores released determined microscopically. Samples were collected every 2 h during the fermentation.

Media

The media was adjusted for all the fermentations to a carbon:nitrogen ratio of 7:1; to achieve this ratio several glucose-soybean meal concentrations were used (Table 2). Additionally, all media contains: 5.8 g/L yeast extract, 9.2 g/L corn steep liquor, 3 g/L KCl, 0.2 mg/L MgSO4 • 7H2O, 40 mg/L MnSO4, 30 mg/L CoCl2, 7.5 mg/L CuSO4 • 5H2O, 7.5 g/L ZnSO4 • 7H2O, 1.35 mg/L FeSO4 and 7 ml/L H3PO4 85% (w/v).

Table 2 Parameters estimated with structural model

Analytical methods

Bacillus and endospore count

The counting of cells and endospores was done in a Neubauer chamber with a Nikon Eclipse 55i microscope (Nikon, Tokyo, Japan) under dark field illumination, with a 40 X objective and 10 X ocular lenses, as previously described (Navarro et al. [2008]). In order to express the number of cells as biomass concentration it was consider that the weight of one B. thuringiensis cell is 2.3 pg (Rodríguez and de la Torre [1996]).

Poly-β-hydroxybutyrate

The PHB concentration was determined in a FID gas chromatograph (SRI Instrument, SRI 8610C, USA) equipped with an AT-1000 column (Navarro et al. [2006]). Each value reported is the average of four replicates. The PHB content, reported as percentage of dry cell weight was calculated with the maximum PHB and biomass concentration obtained in each culture.

Cry protein

The concentration of Cry protein was determined by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE). Bovine serum albumin (BSA) was used as protein standard. The samples preparation was done according to Navarro et al. ([2006]). The solubilization process was done using the method described by Farrera et al. ([1998]). Each value reported is the average of four replicates. The Cry percentage of dry cell weight was calculated with the maximum Cry and biomass concentration obtained in each culture.

Dipicolinic acid

Pyridine-2,6-dicarboxylic acid (dipicolinic acid-DPA) is a unique constituent of all endospores of Bacillus and Clostridium genus. In order to obtain more accurate result, endospore kinetics was follow by DPA photoluminescence technique (Navarro et al. [2008]). The photoluminescence of all samples were measured by duplicate in a spectrofluorometer (SLM 48000C Instrument, Urban IL).

Results

The normalized data from the literature show that the cells/biomass and PHB kinetics (Additional file 1: Figures S1 and S2) always have the same behavior and that they are independent to: culture systems, operating conditions, culture media, bioreactors used or the techniques used to measure. The Gompertz model was statistically sufficient to describe the growth and PHB data for all previously reported kinetics (Additional file 1).

Experimental key compounds kinetics analysis for model validation

Vegetative growth phase

The cell kinetics for B. thuringiensis var. kurstaki HD-73 cultivations are shown in Figure 2a. The cell kinetics show the different phases of the culture (Figure 2a): lag phase (first 3 ± 1 h of culture), vegetative growth (from 3 ± 1 h to 7 ± 1 h of culture), transition (from 7 ± 1 h to 9 ± 1 h of culture) and sporulation (after 9 ± 1 h of culture). The productivity from the experimental data and from the literature in batch culture is shown in Table 3.

Figure 2
figure 2

Experimental and simulated kinetics of key compounds of B. thuringiensis var. kurstaki HD-73. Key compound kinetics: (a) biomass; (b) Poly-β-hidroxybutyrate; (c) Endospores and dipicolinic acid (DPA); and (d) Cry protein. Experimental kinetics of the key compounds (dots) obtained at different glucose and soybean meal concentration. Simulated kinetics of key compounds calculated with the biodynamic model (dotted lines).

Table 3 Parameters culture of Bacillus thuringiensis reported in the literature

Transition phase

The PHB kinetics for B. thuringiensis var. kurstaki HD-73 cultivations are shown in Figure 2b. During the first 6 ± 1 h of culture no PHB was detected, then its accumulation begins and reaches a maximum at 11 ± 1 h. After this time, PHB concentration gradually diminishes until the end of the culture, this means that the PHB consumption outstrips its production (Figure 2b). The maximum production of PHB in the batch culture of B. thuringiensis is far away from the highest production values but agrees with those reported in similar experiments (Table 3).

Sporulation phase

The endospore and DPA kinetics for B. thuringiensis var. kurstaki HD-73 cultivations are shown in Figure 2c. From the first 10 to 14 hours of culture no endospore or DPA formation were detected. Then, the endospore and DPA are accumulated until the end of the culture (Figure 2c).

The kinetics of Cry protein during B. thuringiensis var. kurstaki HD-73 cultivations are shown in Figure 2d. From the first 10 to 12 hours of culture there is no Cry formation. Then, the Cry protein is accumulated until the end of the culture (Figure 2d). At the same time, cells are dying, so active biomass diminish, there is no carbon source consumption and the PHB is consumed (Figure 2).

Biodynamic model validation

The model results fit well the experimental data (Figure 2). In general, biodynamic model predicts the kinetics of all selected key compounds (X, PHB, Cry, DPA). The comparison between the experimental and the simulated data is shown in Figure 3a,b,c,d the deviation between them is within the experimental deviations values and they can be fitted with a straight line. The coefficient of determination for all kinetics was always greater than 90%. The parameters obtained with the biodynamic model constructed are listed in Table 2. For biodynamic model validation it is very important that the difference between the values of the parameters obtained theoretically and experimentally were near to zero.

Figure 3
figure 3

Comparison of simulated and experimental key compound concentrations of B. thuringiensis var. kurstaki HD-73. Batch fermentations at different glucose and soybean meal concentration: F1 at 25.1 & 4.4 g/L; F2 at 34.74 & 14.73 g/L; F3 at 44.42 & 25.05 g/L; and F4 at 54.1 & 35.37 g/L, respectively. Comparison between simulated and experimental data of (a) biomass; (b) Poly-B-hydroxybutyrate; (c) dipicolinic acid (DPA); and (d) Cry protein.

Discussion

Key compounds productivity

In general, the experimental data of the biomass and PHB production agree with those previously reported. Avignone-Rossa and Mignone ([1995]) reported the highest cells productivity for B. thuringiensis (1.6 × 1010 CFU/ml or 36.8 g/L, Table 3) up to date. However, this concentration was obtained in a Fed-Batch system. Pal et al. ([2009]) and Rohini et al. ([2006]) reported two of the highest percentages of PHB per gram dry weight of B. thuringiensis cells, 60.3% and 64.1% respectively (Table 3). But, it is worth to mention that these experiments were conducted mainly to improve PHB production.

The sequence of endospore development in B. thuringiensis var. kurstaki is a seven stages process (Bechtel and Bulla [1976]).

Stage I starts at the beginning of the sporulation process (t0); stage II, the first hour; stage III, 1 h after t0; stage IV, 2 h after t0; stage V, 4 h after t0; stage VI, 5 h after t0; and stage VII, more than 5 h after t0. Figure 2c shows that, in general, the DPA is detected two hours earlier than endospore. This would be expected because during stage V, DPA synthetase is encoded by spoVF (Paidhungat et al. [2000]), and after stage VII, mature endospores are released (Bechtel and Bulla [1976]). During sporulation process, cells are dying, viable cells decrease, the PHB is consumed, but there is no carbon source consumption (Figure 2). Therefore, the PHB becomes the carbon and energy source for Cry protein production as it have been suggested by Benoit et al. ([1990]), Kraemer-Schafhalter and Moser ([1996]), Navarro et al. ([2006]) and Rowe and Margaritis ([1987]). Figures 2c and d show that, in general, the Cry protein is detected before than endospore. This would be expected because the Cry protein is first observed during stage III and the crystal protein is almost full-sized during stage IV of sporulation (Bechtel and Bulla [1976]).

According to the experimental results, the Cry protein represents approximately 6% of the dry cell weight. This result is below some previously reported values (Table 3). However, there is not enough reports using SDS-PAGE for Cry protein quantification and therefore it is difficult to establish a comparison. The Cry protein concentrations reported in literature are scattered and confused. For example, while Rowe and Margaritis ([1987]) reported that after cell lysis the endospore and Cry protein mass were 15% of the biomass dry weight, Agaisse and Lereclus ([1995]) reported that Cry protein can account for up to 25% of the dry weight of the sporulated cells. Even some values from the literature of Cry protein concentration are much higher, for example: Anderson and Jayaraman ([2005]) reported 98% of the biomass dry weight (calculated from the reported values of 15.7 g/L of Cry protein from 16 g/L of biomass, Table 3). Considering these data, the yields coefficient of Cry from biomass are 0.98 g Cry/g biomass (Table 3), which means that almost every cell is transformed into Cry protein. In the same way, the yield coefficient of Cry from biomass obtained by Chang et al. ([2008]) and Dang Vu et al. ([2009]) seems excessively high (Table 3).

Biodynamic model parameters

The theoretical values of the parameters like the maximum concentrations of the key compounds (Xmax, PHBmax and Crymax) agree with the experimental and the reported ones (Tables 2 and 3, Additional file 1: Figures S3 and S5). Different maximum specific growth rate (μmax) and growth yield coefficients (YX/S) have been reported for different culture conditions (Table 3). For example, Popovic et al. ([2001]) reported a growth yield coefficient of 2.25 g active biomass/g limiting substrate (assuming glucose as the limiting substrate). According to the stoichiometry for anaerobic growth, the maximum YX/S expected would be 0.863 g biomass/g glucose base on the general chemical composition of biomass (CH1.8O0.5N0.2), 1 C-mol of glucose (CH2O) and no respiration process (no formation of CO2). Then, the growth yield coefficient for aerobic growth must be always below 0.863 g of biomass/g glucose and values greater than 0.863 imply that glucose is not the limiting substrate. In the present study, the μmax and the YX/S values obtained agree with previously reported (Tables 2 and 3). The Xmax, PHBmax, DPAmax, Crymax, tc and tcp relationship with the initial glucose concentration were analyzed from the four new fermentations. The Xmax, PHBmax and Crymax have a linear relationship with the glucose concentration in the media, but no relationship with DPAmax was found (Table 2). The lag period of each key compound (indirectly represented by tc) shows a linear relationship with the glucose concentration in the media (Table 2). The same behavior was observed with literature data. The tcx y tcp obtained with the model for fermentation in literature with high glucose concentration in formulation is greater than those from fermentation with low concentration (Additional file 1: Figure S4 and S6). There is no relationship between the specific rates of the key compounds with their maximum concentrations (Table 2).

Dynamics of B. thuringiensis culture

The biodynamic model represents the evolution of the B. thuringiensis fermentation and it was used to study the dynamic of the process. The three phases of the fermentation are related between them and the Cry protein production is the result of the interaction between those phases. Therefore, the endospore formation is related with the vegetative growth and the endospore-Cry formation is related with the transition phase. This means that the bacilli disappearance is related with DPA (endospore) formation through YDPA/X. In the same way, the PHB consumption is related with DPA-Cry production through YDPA/PHB and YCry/PHB. These yields were calculated with the model (Table 2) and they show a different relationship with the glucose concentration in the media. The YCry/PHB increases and the YDPA/PHB decreases when the initial nutrient concentration increases. These results suggest that endospore formation consumes preferably the PHB and excess (if any) is used in the Cry protein formation. There is not data available in the literature to compare these results and further experiments have to be done. Also, the biodynamic model simulation results indicate that the Cry protein production is a function of development of the three fermentation phases, and that endospore and Cry protein production are competitive processes. Therefore, the optimization of each phase cannot be done independently. The reports in the literature usually do not take into account the transition phase neither the PHB production as important variables for Cry production. However, the biodynamic model implies that PHB is crucial for endospore and Cry protein production.

In summary, the new biodynamic model reported here represents well the dynamic of the B. thuringiensis fermentation and it can be used as a tool to understand the complex processes related to the physiology of the bacilli. Besides, it can be applied for fermentation with different production systems, operating conditions and culture media. Also, it can be used to study other dynamic processes that occur during fermentation as the cell cycle regulation, the signaling and gen regulatory networks and the changes in gene expression, between others. Although, the biodynamic model is the first model that represents the entire fermentation of B. thuringiensis and it is also the first one that links each phase of the culture with its predominantly active metabolic pathway. This model is a powerful tool to study the dynamics of the B. thuringiensis culture and the interrelation between their phases. Also, the model can be used to know the physiological phase during the fermentation. From the simulation results it can be concluded that the PHB is crucial for endospore and Cry protein production; and that those processes are simultaneous, parallel and competitive. In terms of the DPA and Cry yield from PHB, the endospore formation is preferred over the Cry protein synthesis.

Additional file

References

  • Agaisse H, Lereclus D: How does Bacillus thuringiensis produce so much insecticidal crustal protein? J Bacteriol 1995, 177: 6027–6032.

    PubMed Central  CAS  PubMed  Google Scholar 

  • Anderson TB: Effects of carbon:nitrogen ratio and oxygen on the growth kinetics of Bacillus thuringiensis and yield of bioinsecticidal crystal protein. University of Western Ontario, Dissertation; 1990.

    Google Scholar 

  • Anderson RKI, Jayaraman K: Impact of balanced substrate flux on the metabolic process employing fuzzy logic during the cultivation of Bacillus thuringiensis var. Galleriae . World J Microbiol Biotechnol 2005, 21: 127–133. 10.1007/s11274-004-3043-1

    Article  CAS  Google Scholar 

  • Anderson RKI, Jayaraman K, Voisard D, Marison IW, von Stockar U: Heat flux as an on-line indicator of metabolic activity in pilot scale bioreactor during the production of Bacillus thuringiensis var. galleriae -based biopesticides. Thermochim Acta 2002, 386: 127–138. 10.1016/S0040-6031(01)00709-2

    Article  CAS  Google Scholar 

  • Aslim B, Yüksekdağ ZN, Beyatli Y: Determination of PHB growth quantities of certain Bacillus species isolated from soil. Turk Electron J Biotechnol 2002, (Special Issue):24–30. http://www.biyotekder.hacettepe.edu.tr/5.pdf

  • Avignone-Rossa C, Mignone CF: Bacillus thuringiensis growth and toxicity. Basic and applied considerations. Mol Biotechnol 1995, 4: 55–71. 10.1007/BF02907471

    Article  CAS  PubMed  Google Scholar 

  • Bechtel DB, Bulla LA: Electron microscope study of sporulation and parasporal crystal formation in Bacillus thuringiensis . J Bacteriol 1976, 127: 1472–1481.

    PubMed Central  CAS  PubMed  Google Scholar 

  • Benoit DG, Wilson GR, Baugh CL: Fermentation during growth and sporulation of Bacillus thuringiensis HD-1. Lett Appl Microbiol 1990, 10: 15–18. 10.1111/j.1472-765X.1990.tb00084.x

    Article  CAS  Google Scholar 

  • Berbert-Molina MA, Prata AMR, Pessanha LG, Silveira MM: Kinetics of Bacillus thuringiensis var. Israelensis growth on high glucose concentration. J Ind Microbiol Biotechnol 2008, 35: 1397–1404. 10.1007/s10295-008-0439-1

    Article  CAS  PubMed  Google Scholar 

  • Chang M, Zhou SG, Lu N, Ni JR: Starch processing wastewater as a new medium for production of Bacillus thuringiensis . World J Microbiol Biotechnol 2008, 24: 441–447. 10.1007/s11274-007-9491-7

    Article  CAS  Google Scholar 

  • Dang Vu K, Tyagi RD, Valero JR, Surampalli RY: Impact of different pH control agents on biopesticidal activity of Bacillus thuringiensis during the fermentation of starch industry wastewater. Bioproc Biosyst Eng 2009, 32: 511–519. 10.1007/s00449-008-0271-z

    Article  Google Scholar 

  • Farrera RR, Pérez-Guevara F, de la Torre M: Carbon:nitrogen ratio interacts with initial concentration of total solids on insecticidal crystal protein and spore production in Bacillus thuringiensis HD-73. Appl Microbiol Biotechnol 1998, 49: 758–765. 10.1007/s002530051243

    Article  CAS  Google Scholar 

  • Holmberg A, Sievänen R: Fermentation of Bacillus thuringiensis for exotoxin production: process analysis study. Biotechnnol Bioeng 1980, 22: 1707–1724. 10.1002/bit.260220812

    Article  CAS  Google Scholar 

  • Ito A, Sasaguri Y, Kitada S, Kusaka Y, Kuwano K, Masutomi K, Mizuki E, Akao T, Ohba M: A Bacillus thuringiensis crystal protein with selective cytocidal action to human cells. J Biol Chem 2004, 279: 21282–21286. 10.1074/jbc.M401881200

    Article  CAS  PubMed  Google Scholar 

  • Kraemer-Schafhalter A, Moser A: Kinetic study of Bacillus thuringiensis var. israelensis in lab-scale process. Bioproc Eng 1996, 14: 139–144. 10.1007/BF00369431

    Article  Google Scholar 

  • Liu B, Tzeng Y: Characterization study of the sporulation kinetics of Bacillus thuringiensis . Biotechnol Bioeng 2000, 68: 11–17. 10.1002/(SICI)1097-0290(20000405)68:1<11::AID-BIT2>3.0.CO;2-T

    Article  CAS  PubMed  Google Scholar 

  • Navarro AK, Farrera RR, Lopez R, Perez-Guevara F: Relationship between Poly-β-hydroxybutyrate production and δ-endotoxin for Bacillus thuringiensis var. kurstaki . Biotechnol Lett 2006, 28: 641–644. 10.1007/s10529-006-0029-0

    Article  CAS  PubMed  Google Scholar 

  • Navarro AK, Peña A, Pérez-Guevara F: Endospore-Dipicolinic acid detection during Bacillus thuringiensis culture. Lett Appl Microbiol 2008, 46: 166–170. 10.1111/j.1472-765X.2007.02277.x

    Article  CAS  PubMed  Google Scholar 

  • Ozcan O, Icgenb B, Ozcengiz G: Pretreatment of poultry litter improves Bacillus thuringiensis -based biopesticides production. Bioresour Technol 2010, 101: 2401–2404. 10.1016/j.biortech.2009.11.048

    Article  CAS  PubMed  Google Scholar 

  • Paidhungat M, Setlow B, Driks A, Setlow P: Characterization of spores of Bacillus subtilis which lack dipicolinic acid. J Bacteriol 2000, 182: 5505–5512. 10.1128/JB.182.19.5505-5512.2000

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  • Pal A, Prabhu A, Kumar AA, Rajagopal B, Dadhe K, Ponnamma V, Shivakumar S: Optimization of process parameters for maximum poly(-β-)hydroxybutyrate (PHB) production by Bacillus thuringiensis IAM 12077. Pol J Microbiol 2009, 58: 149–154.

    CAS  PubMed  Google Scholar 

  • Popovic M, Liu W, Iannotti EL, Bajpai RK: A mathematical model for vegetative growth of Bacillus thuriniensis . Eng Life Sci 2001, 2: 85–90. 10.1002/1618-2863(200108)1:2<85::AID-ELSC85>3.0.CO;2-M

    Article  Google Scholar 

  • Rivera D, Margaritis A, De Lasa H: A sporulation kinetic model for bacth growth of B. thuringiensis . Can J Chem Eng 1999, 77: 903–910. 10.1002/cjce.5450770516

    Article  CAS  Google Scholar 

  • Rodríguez MM, de la Torre M: Effect of the dilution rate on the biomass yield of Bacillus thuringiensis and determination of its rate coefficients under steady-state conditions. Appl Microbiol Biotechnol 1996, 45: 546–550. 10.1007/BF00578469

    Article  Google Scholar 

  • Rohini D, Phadnis S, Rawal SK: Synthesis and characterization of Poly-β-hydroxybutyrate from Bacillus thuringiensis R1. Indian J Biotechnol 2006, 5: 276–283.

    CAS  Google Scholar 

  • Rowe GE: Central metabolism of Bacillus thuringiensis during growth and sporulation. University of Western Ontario, Dissertation; 1990.

    Google Scholar 

  • Rowe GE, Margaritis A: Bioprocess development in the production of bioinsecticides by Bacillus thuringiensis . CRC Crit Rev Biotechnol 1987, 6: 87–127. 10.3109/07388558709086986

    Article  CAS  Google Scholar 

  • Starzak M, Bajpai RK: A structured model for vegetative growth and sporulation in Bacillus thuringiensis . Appl Biochem Biotechnol 1991, 28(29):699–718. 10.1007/BF02922643

    Article  PubMed  Google Scholar 

  • Zhuang L, Zhou S, Wang Y, Chang M: Mosquito biolarvicide production by sequential fermentation with dual strains of Bacillus thuringiensis subsp. Israelensis and Bacillus sphaericus using sewage sludge. Bioresour Technol 2011, 102: 1574–1580. 10.1016/j.biortech.2010.08.090

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

The authors acknowledge the funding for this research provided by Consejo Nacional de Ciencia y Tecnología (2003-43127 and 2008-105057) and by Programa de Mejoramiento del Profesorado (103.5/10/0246). A.K. Navarro-Mtz also acknowledges the fellowship provided by Consejo Nacional de Ciencia y Tecnología.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ana Karin Navarro-Mtz.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AKN-M had done substantial contributions to the acquisition, analysis and interpretation of data and she had been involved in drafting the manuscript. FP-G had done substantial contributions to the conception, experiment design and data analysis. Both authors read and approved the final manuscript.

Electronic supplementary material

13568_2014_79_MOESM1_ESM.pdf

Additional file 1: Normalization of cells/biomass and PHB kinetics from B. thuringiensis culture. Figure S1. Normalization of the biomass kinetics from reported data. The data correspond to different culture systems, different operating conditions, different culture media, different bioreactors, etc., and biomass concentrations were measured using different techniques. The line represents an average biomass kinetic for B. thuringiensis var. kurstaki HD-73. Figure S2. Normalization of the Poly-β-hidroxybutyrate kinetics from reported data. The data correspond to different culture systems, different operating conditions, different culture media, different bioreactors, etc. and Poly-β-hidroxybutyrate concentrations were determined using different techniques. The line represents an average PHB kinetic for B. thuringiensis var. kurstaki HD-73. Figure S3. Effect of the glucose concentration on the Xmax for B. thuringiensis culture (bibliographic data). Figure S4. Effect of the glucose concentration on the tcx for B. thuringiensis culture (bibliographic data). Figure S5. Effect of the glucose concentration on the PHBmax and for B. thuringiensis culture (bibliographic data). Figure S6. Effect of the glucose concentration on the tcp for B. thuringiensis culture (bibliographic data). (PDF 9 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0), which permits use, duplication, adaptation, distribution, and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Navarro-Mtz, A.K., Pérez-Guevara, F. Construction of a biodynamic model for Cry protein production studies. AMB Expr 4, 79 (2014). https://doi.org/10.1186/s13568-014-0079-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13568-014-0079-y

Keywords