Acessibilidade / Reportar erro

Fermentation optimization of rennet-producing Bacillus amyloliquefaciens GSBa-1 for high-density culture and its kinetic model

Abstract

The milk-clotting enzyme (MCE) produced by Bacillus amyloliquefaciens GSBa-1 achieved the scale-up criterion from 0.25-L erlenmeyer flask to 5-L bioreactor sucessfully and exhibited remarkable milk-clotting activity (MCA) (626.3 SU/mL) that was 1.45-fold higher than the control by optimizing both medium composition and fermentation conditions.The growth of GSBa-1 was optimized through response surface methodology employing Plackett-Burman design and Box-Behnken design. The final optimized fermentation conditions were as follows: maltose (7.258g/L), corn steep liquor (8.15g/L), and Tween 80 (1.97g/L); constant fermentation pH (7.0), temperature (45 °C) and stirrer speed (200 rpm). Kinetic models of microbial growth (X), product formation (P), and substrate consumption (S) were constructed, respectively. The derived model was the first reported model for GSBa-1 rennet production, contributing visual description for the rennet fermentation process, as well as the improvement of its productivity and efficiency.

Keywords:
Bacillus amyloliquefaciens GSBa-1; rennet; fermentation optimization; high-density culture; kinetic model

1 Introduction

Bacillus amyloliquefaciens belongs to the super kingdom of bacteria with a Bacillaceae family and genus of Bacillus (Teng et al., 2017Teng, J., Zhao, X., Yang, Y., Zhang, J., Zhao, A., Jiang, Y., Li, L., Zheng, Z., & Yang, Z. (2017). Isolation and identification of microbial strains producing rennet from Jiuqu, a traditional Chinese fermentation starter. Shipin Kexue, 38(16), 23-28.; WoldemariamYohannes et al., 2020WoldemariamYohannes, K., Wan, Z., Yu, Q., Li, H., Wei, X., Liu, Y., Wang, J., & Sun, B. (2020). Prebiotic, probiotic, antimicrobial, and functional food applications of Bacillus amyloliquefaciens. Journal of Agricultural and Food Chemistry, 68(50), 14709-14727. http://dx.doi.org/10.1021/acs.jafc.0c06396. PMid:33280382.
http://dx.doi.org/10.1021/acs.jafc.0c063...
). Even though B. amyloliquefaciens has different strains with different properties and applications, it is considered as a safe and nontoxic producing microbe (Chen et al., 2020Chen, Y.-C., Huang, S.-D., Tu, J.-H., Yu, J.-S., Nurlatifah, A. O., Chiu, W.-C., Su, Y.-H., Chang, H.-L., Putri, D. A., & Cheng, H.-L. (2020). Exopolysaccharides of Bacillus amyloliquefaciens modulate glycemic level in mice and promote glucose uptake of cells through the activation of Akt. International Journal of Biological Macromolecules, 146, 202-211. http://dx.doi.org/10.1016/j.ijbiomac.2019.12.217. PMid:31887379.
http://dx.doi.org/10.1016/j.ijbiomac.201...
; Shahzad et al., 2020Shahzad, R., Shehzad, A., Bilal, S., & Lee, I.-J. (2020). Bacillus amyloliquefaciens RWL-1 as a new potential strain for augmenting biochemical and nutritional composition of fermented soybean. Molecules, 25(10), 2346. http://dx.doi.org/10.3390/molecules25102346. PMid:32443519.
http://dx.doi.org/10.3390/molecules25102...
; Tran et al., 2020Tran, T. D., Cid, C., Hnasko, R., Gorski, L., & McGarvey, J. A. (2020). Bacillus amyloliquefaciens ALB65 inhibits the growth of Listeria monocytogenes on cantaloupe melons. Applied and Environmental Microbiology, 87(1), e01920-e01926. http://dx.doi.org/10.1128/AEM.01926-20. PMid:33097500.
http://dx.doi.org/10.1128/AEM.01926-20...
). The rennet produced by B. amyloliquefacies GSBa-1 is a promising alternative to calf rennet in cheese processing (Teng et al., 2019Teng, J.-W., Liu, Z.-M., Su, M.-Y., Gao, H.-Y., Zheng, Y.-R., Li, L., Zheng, Z., Zhang, J., & Yang, Z.-N. (2019). Application of rennet from Bacillus amyloliquefaciens GSBa-1 in mozzarella cheese. Shipin Kexue, 40, 85-92.; Zhao et al., 2020aZhao, D., Liu, L., Jiang, J., Guo, S., Ping, W., & Ge, J. (2020a). The response surface optimization of exopolysaccharide produced by Weissella confusa XG-3 and its rheological property. Preparative Biochemistry & Biotechnology, 50(10), 1014-1022. http://dx.doi.org/10.1080/10826068.2020.1780609. PMid:32589090.
http://dx.doi.org/10.1080/10826068.2020....
). In order to apply this rennet in cheese production properly, it is necessary to optimize the conditions for the high milk-clotting activity rennet produced by B. amyloliquefaciens GSBa-1(Mota et al., 2018Mota, M. J., Lopes, R. P., Koubaa, M., Roohinejad, S., Barba, F. J., Delgadillo, I., & Saraiva, J. A. (2018). Fermentation at non-conventional conditions in food- and bio-sciences by the application of advanced processing technologies. Critical Reviews in Biotechnology, 38(1), 122-140. http://dx.doi.org/10.1080/07388551.2017.1312272. PMid:28423948.
http://dx.doi.org/10.1080/07388551.2017....
). High-density fermentation technology has been widely used in the fermentation optimization of eukaryotic microorganisms, prokaryotic microorganisms and algae, which can significantly increase the yield of microorganisms (Freudenberg et al., 2021Freudenberg, R. A., Baier, T., Einhaus, A., Wobbe, L., & Kruse, O. (2021). High cell density cultivation enables efficient and sustainable recombinant polyamine production in the microalga Chlamydomonas reinhardtii. Bioresource Technology, 323, 124542. http://dx.doi.org/10.1016/j.biortech.2020.124542. PMid:33385626.
http://dx.doi.org/10.1016/j.biortech.202...
; Kleman & Strohl, 1994Kleman, G. L., & Strohl, W. R. (1994). Developments in high cell density and high productivity microbial fermentation. Current Opinion in Biotechnology, 5(2), 180-186. http://dx.doi.org/10.1016/S0958-1669(05)80033-3. PMid:7764797.
http://dx.doi.org/10.1016/S0958-1669(05)...
; Liu et al., 2019Liu, W.-C., Inwood, S., Gong, T., Sharma, A., Yu, L.-Y., & Zhu, P. (2019). Fed-batch high-cell-density fermentation strategies for Pichia pastoris growth and production. Critical Reviews in Biotechnology, 39(2), 258-271. http://dx.doi.org/10.1080/07388551.2018.1554620. PMid:30599783.
http://dx.doi.org/10.1080/07388551.2018....
; Wang et al., 2020Wang, T., Lu, Y., Yan, H., Li, X., Wang, X., Shan, Y., Yi, Y., Liu, B., Zhou, Y., & Lü, X. (2020). Fermentation optimization and kinetic model for high cell density culture of a probiotic microorganism: Lactobacillus rhamnosus LS-8. Bioprocess and Biosystems Engineering, 43(3), 515-528. http://dx.doi.org/10.1007/s00449-019-02246-y. PMid:31712884.
http://dx.doi.org/10.1007/s00449-019-022...
). And the fermentation of microorganism is highly dependent on the composition of medium as well as their fermentation conditions (Ju et al., 2018Ju, Y., Son, K.-H., Jin, C., Hwang, B. S., Park, D.-J., & Kim, C.-J. (2018). Statistical optimization of culture medium for improved production of antimicrobial compound by Streptomyces rimosus AG-P1441. Food Science and Biotechnology, 27(2), 581-590. PMid:30263783.; Martinez-Burgos et al., 2021Martinez-Burgos, W. J., Sydney, E. B., Paula, D. R., Medeiros, A. B. P., Carvalho, J. C., Molina, D., & Soccol, C. R. (2021). Hydrogen production by dark fermentation using a new low-cost culture medium composed of corn steep liquor and cassava processing water: process optimization and scale-up. Bioresource Technology, 320(Pt B), 124370. http://dx.doi.org/10.1016/j.biortech.2020.124370. PMid:33220544.
http://dx.doi.org/10.1016/j.biortech.202...
; Wei et al., 2009Wei, G., Yang, X., Gan, T., Zhou, W., Lin, J., & Wei, D. (2009). High cell density fermentation of Gluconobacter oxydans DSM 2003 for glycolic acid production. Journal of Industrial Microbiology & Biotechnology, 36(8), 1029-1034. http://dx.doi.org/10.1007/s10295-009-0584-1. PMid:19434434.
http://dx.doi.org/10.1007/s10295-009-058...
).

Medium optimization is a very practical and valuable tool for fermentation industry to improve product yield and minimize by-products as well as reduce overall manufacturing costs. The metabolites produced by microorganisms mainly depend on the strain, the composition of the culture medium and culture conditions. Medium composition has a significant impact on the formation of target products (Jung & Lee, 2020Jung, S., & Lee, J.-H. (2020). Characterization of transcriptional response of Lactobacillus plantarum under acidic conditions provides insight into bacterial adaptation in fermentative environments. Scientific Reports, 10(1), 19203. http://dx.doi.org/10.1038/s41598-020-76171-6. PMid:33154427.
http://dx.doi.org/10.1038/s41598-020-761...
; Rajeswari et al., 2015Rajeswari, P., Jose, P. A., Amiya, R. A., & Jebakumar, S. R. D. (2015). Characterization of saltern based Streptomyces sp. and statistical media optimization for its improved antibacterial activity. Frontiers in Microbiology, 5, 753. PMid:25653640.). Statistical optimization is a good mean screening significant variables in a large experimental design (Abdel-Fattah et al., 2005Abdel-Fattah, Y. R., Saeed, H. M., Gohar, Y. M., & El-Baz, M. A. (2005). Improved production of Pseudomonas aeruginosa uricase by optimization of process parameters through statistical experimental designs. Process Biochemistry, 40(5), 1707-1714. http://dx.doi.org/10.1016/j.procbio.2004.06.048.
http://dx.doi.org/10.1016/j.procbio.2004...
; Saavedra et al., 2021Saavedra, S., Alejandro-Paredes, L., Flores-Santos, J. C., Flores-Fernández, C. N., Arellano-García, H., & Zavaleta, A. I. (2021). Optimization of lactic acid production by Lactobacillus plantarum strain Hui1 in a medium containing sugar cane molasses. Food Science and Technology, 39(1), 98-107.). Variables with statistically significant effects on milk-clotting activity were identified by screening using the Plackett-Burman design (PBD) (Elsayed & Abdelwahed, 2020Elsayed, E. A., & Abdelwahed, N. A. (2020). Medium optimization by response surface methodology for improved cholesterol oxidase production by a newly isolated streptomyces rochei NAM-19 strain. BioMed Research International, 2020, 1870807. http://dx.doi.org/10.1155/2020/1870807. PMid:32337225.
http://dx.doi.org/10.1155/2020/1870807...
; Jiang et al., 2020Jiang, J., Zu, Y., Li, X., Meng, Q., & Long, X. (2020). Recent progress towards industrial rhamnolipids fermentation: process optimization and foam control. Bioresource Technology, 298, 122394. http://dx.doi.org/10.1016/j.biortech.2019.122394. PMid:31757615.
http://dx.doi.org/10.1016/j.biortech.201...
; Plackett & Burman, 1946Plackett, R. L., & Burman, J. P. (1946). The design of optimum multifactorial experiments. Biometrika, 33(4), 305-325. http://dx.doi.org/10.1093/biomet/33.4.305.
http://dx.doi.org/10.1093/biomet/33.4.30...
). These were further optimized by Box-Behnken Design (BBD) experiments (Yang et al., 2020Yang, X., Guo, P., Li, M., Li, H., Hu, Z., Liu, X., & Zhang, Q. (2020). Optimization of culture conditions for amoxicillin degrading bacteria screened from pig manure. International Journal of Environmental Research and Public Health, 17(6), 1973. http://dx.doi.org/10.3390/ijerph17061973. PMid:32192171.
http://dx.doi.org/10.3390/ijerph17061973...
) and Response Surface Methodology (RSM) (Ibrahim et al., 2019Ibrahim, A. A., El-Housseiny, G. S., Aboshanab, K. M., Yassien, M. A., & Hassouna, N. A. (2019). Paromomycin production from Streptomyces rimosus NRRL 2455: statistical optimization and new synergistic antibiotic combinations against multidrug resistant pathogens. BMC Microbiology, 19(1), 18. http://dx.doi.org/10.1186/s12866-019-1390-1. PMid:30658584.
http://dx.doi.org/10.1186/s12866-019-139...
; Zhao et al., 2020bZhao, X., Cai, M., Yang, Z., Luo, T., Chen, C., Cao, Y., & Yang, Z. (2020b). Proteolysis and bioactivity of cheddar cheese as influenced by milk-clotting enzyme from Bacillus amyloliquefaciens GSBa-1. Shipin Kexue, 41, 103-111.).

Kinetic modeling has long been used to provide crucial information about the metabolic capabilities of microorganisms during their cultivation (Olivares-Marin et al., 2018Olivares-Marin, I. K., González-Hernández, J. C., Regalado-Gonzalez, C., & Madrigal-Perez, L. A. (2018). Saccharomyces cerevisiae exponential growth kinetics in batch culture to analyze respiratory and fermentative metabolism. Journal of Visualized Experiments, 139, 58192. PMid:30320748.; Zhang et al., 2019Zhang, Y., Xu, S., Cui, M., & Wong, J. W. C. (2019). Effects of different thermal pretreatments on the biodegradability and bioaccessibility of sewage sludge. Waste Management, 94, 68-76. http://dx.doi.org/10.1016/j.wasman.2019.05.047. PMid:31279397.
http://dx.doi.org/10.1016/j.wasman.2019....
). In order to further study the internal laws and dynamic changes during the fermentation of B. amyloliquefaciens GSBa-1, mathematical models were used to quantitatively analyze the changes in bacterial growth, milk-clotting activity and the consumption of total sugars during the fermentation process.

In our previous research, the rennet from B. amyloliquefaciens GSBa-1 has been proved to have better metal ion chelating ability and antioxidant activity compared with commercial rennet. It improved biological activity of cheese while meeting general cheese requirements (Zhao et al., 2020aZhao, D., Liu, L., Jiang, J., Guo, S., Ping, W., & Ge, J. (2020a). The response surface optimization of exopolysaccharide produced by Weissella confusa XG-3 and its rheological property. Preparative Biochemistry & Biotechnology, 50(10), 1014-1022. http://dx.doi.org/10.1080/10826068.2020.1780609. PMid:32589090.
http://dx.doi.org/10.1080/10826068.2020....
). In addition, GSBa-1 rennet cheese has slightly stronger proteolytic activity, which leads to more types flavors than those made with commercial rennet (Teng et al., 2019Teng, J.-W., Liu, Z.-M., Su, M.-Y., Gao, H.-Y., Zheng, Y.-R., Li, L., Zheng, Z., Zhang, J., & Yang, Z.-N. (2019). Application of rennet from Bacillus amyloliquefaciens GSBa-1 in mozzarella cheese. Shipin Kexue, 40, 85-92.). Therefore, GSBa-1 rennet has great application potential in cheese production. In this study, the high-density culture strategy was approached by optimizing the composition of the medium and fermentation conditions, taking the milk-clotting activity as an indicator. At the same time, monitoring and understanding the fermentation characteristics of B. amyloliquefaciens GSBa-1 and establishing its fermentation kinetic model laid the foundation for the application of B. amyloliquefaciens GSBa-1 rennet in food industry.

2 Materials and method

2.1 Microorganisms and medium

B. amyloliquefaciens GSBa-1 was isolated from wine koji. The strain was stored in 50% glycerol at -80 °C.

The seed culture medium contained the following per liter: tryptone (10 g/L), yeast extract (5 g/L), and sodium chloride (10 g/L). The pH was adjusted to 7.0 and then autoclaved at 121 °C for 20 min.

The fermentation medium contained maltose (6.74 g/L), beef extract (8 g/L), corn syrup (7.5 g/L), Tween 80 (1.9 g/L). The pH was adjusted to 7.0 and then autoclaved at 121 °C for 20 min.

2.2 Erlenmeyer flask culture research

In order to develop a medium suitable for the production of high-activity rennet from B. amyloliquefaciens GSBa-1, the medium(carbon source, nitrogen source and inducer) was optimized in an Erlenmeyer flask.

2.3 Identification of the most significant nutrients by PBD

The PBD plays an important role in the rapid determination of key components in the medium (Singh et al., 2009Singh, A. K., Mehta, G., & Chhatpar, H. S. (2009). Optimization of medium constituents for improved chitinase production by Paenibacillus sp. D1 using statistical approach. Letters in Applied Microbiology, 49(6), 708-714. http://dx.doi.org/10.1111/j.1472-765X.2009.02731.x. PMid:19780958.
http://dx.doi.org/10.1111/j.1472-765X.20...
). The four main variables of the current research, namely corn steep liquor, maltose, beef extract, and Tween 80, were represented by X1, X2, X3, and X4, respectively (Table S1) . Variables were evaluated at high (+1) and low (-1) levels. According to the model, 12 experimental runs were carried out. The PBD was based on the first-order polynomial equation (Equation 1):

Y = β 0 + β i X i (1)

where Y is the activity of GSBa-1 rennet, β0 is the coefficient of the model, βi is the linear coefficient, and Xi is the level of each independent factor.

From the regression analysis, variables that were significant at or above 95% level and at probability value of p < 0.05 were considered to have a great impact on the milk-clotting activity of GSBa-1 and were further statistically optimized by BBD.

2.4 Box-Behnken Design (BBD)

BBD is one of the RSM design with a three-level factorial design, which has been used as the experimental design model to optimize the influential parameters for improving milk-clotting activity (Wang et al., 2019Wang, W., Wu, X., Tan, J., Zhu, L., Mou, Y., Zhang, D., & Gao, J. (2019). Using response surface methodology optimize culture conditions for human lactoferrin production in desert Chlorella. Protein Expression and Purification, 155, 130-135. http://dx.doi.org/10.1016/j.pep.2018.11.004. PMid:30508587.
http://dx.doi.org/10.1016/j.pep.2018.11....
). According to the results of PBD experiment, the levels of the variables and the experimental design were shown in Table S2. BBD was used to investigate the combined effect of three variables such as maltose, corn syrup, and Tween 80 concentration on rennet production. The corresponding fermentation was performed in 500 mL shake flasks with 50 mL medium and incubated at 30 °C in a shaker (160 rpm) for about 24 h. The broth was refrigerated at 4 °C and centrifuged at 10,000 rpm for 10 min, and the supernatant was obtained as crude enzyme solution that was used for milk-clotting activity analysis.

The accuracy of the fitted model was proved by analysis of variance (ANOVA) and the coefficient of R2. The meaning of all terms in the polynomial model is calculated by calculating the statistical judgment F in the probability-value (value) is 0.05. Minitab software 17.0 is used to regress experimental data and generate 3D contour plots. These plots are generated by keeping two variable constants at 0 level and changing other variables within the scope of the experiment. The experiment was performed in triplicate, and the average response was used for analysis. Minitab was used for experimental design and analysis. Design-Expert 11.0 was used to design the experiment and analysis data.

2.5 Fermentation parameters optimization in bioreactor

The culture was used to inoculate the pre-prepared seed culture (3%, v/v) to optimize fermentation of natural pH medium cultivating in the fermenter (BioFlo 115, 5 L vessel volume, China) containing 3 L of minimal medium. After 24 h of fermentation, the broth was refrigerated at 4 °C and centrifuged at 10,000 rpm for 10 min, and the supernatant was obtained as crude enzyme solution that was used for milk-clotting activity analysis. The effects of fermentation conditions (pH, temperature and stirrer speed on milk-clotting activity of B. amyloliquefaciens GSBa-1 were studied to further improve the milk-clotting activity in the bioreactor.

2.6 Kinetic model of B. amyloliquefaciens GSBa-1

At present, various structured and unstructured kinetic models have been applied in microbial fermentation (Phukoetphim et al., 2017Phukoetphim, N., Salakkam, A., Laopaiboon, P., & Laopaiboon, L. (2017). Kinetic models for batch ethanol production from sweet sorghum juice under normal and high gravity fermentations: logistic and modified Gompertz models. Journal of Biotechnology, 243, 69-75. http://dx.doi.org/10.1016/j.jbiotec.2016.12.012. PMid:27988216.
http://dx.doi.org/10.1016/j.jbiotec.2016...
; Procentese et al., 2015Procentese, A., Raganati, F., Olivieri, G., Russo, M. E., Salatino, P., & Marzocchella, A. (2015). Continuous lactose fermentation by Clostridium acetobutylicum-assessment of solventogenic kinetics. Bioresource Technology, 180, 330-337. http://dx.doi.org/10.1016/j.biortech.2015.01.008. PMid:25621726.
http://dx.doi.org/10.1016/j.biortech.201...
; Zhou et al., 2020Zhou, Q., Liu, Y., & Yuan, W. (2020). Kinetic modeling of butyric acid effects on butanol fermentation by Clostridium saccharoperbutylacetonicum. New Biotechnology, 55, 118-126. http://dx.doi.org/10.1016/j.nbt.2019.10.004. PMid:31626983.
http://dx.doi.org/10.1016/j.nbt.2019.10....
). The kinetic modeling of the fermentation process aims to describe the time evolution of different metabolite concentrations by explaining the formation and consumption rates of different metabolites in the metabolic pathway (Buehler & Mesbah, 2016Buehler, E. A., & Mesbah, A. (2016). Kinetic study of acetone-butanol-ethanol fermentation in continuous culture. PLoS One, 11(8), e0158243. http://dx.doi.org/10.1371/journal.pone.0158243. PMid:27486663.
http://dx.doi.org/10.1371/journal.pone.0...
; Thilakavathi et al., 2007Thilakavathi, M., Basak, T., & Panda, T. (2007). Modeling of enzyme production kinetics. Applied Microbiology and Biotechnology, 73(5), 991-1007. http://dx.doi.org/10.1007/s00253-006-0667-0. PMid:17338073.
http://dx.doi.org/10.1007/s00253-006-066...
). Unstructured models make it easier to monitor and predict batch fermentation processes than structured kinetic models (Phukoetphim et al., 2017Phukoetphim, N., Salakkam, A., Laopaiboon, P., & Laopaiboon, L. (2017). Kinetic models for batch ethanol production from sweet sorghum juice under normal and high gravity fermentations: logistic and modified Gompertz models. Journal of Biotechnology, 243, 69-75. http://dx.doi.org/10.1016/j.jbiotec.2016.12.012. PMid:27988216.
http://dx.doi.org/10.1016/j.jbiotec.2016...
).

Logistic equation is usually used to express the growth kinetics of microorganisms, which describes the typical S-shaped curve of strains. The Logistic model is not affected by the substrate, and can better show the influence of the change of the bacterial content on its growth during the fermentation process (Tashiro & Yoshimura, 2019Tashiro, T., & Yoshimura, F. (2019). A neo-logistic model for the growth of bacteria. Physica A, 525(03), 199-215. http://dx.doi.org/10.1016/j.physa.2019.03.049.
http://dx.doi.org/10.1016/j.physa.2019.0...
). The GSBa-1 has an S-type growth curve, so the Logistic equation is used to express the growth law of B. amyloliquefaciens GSBa-1. The logistic equation is (Equation 2):

d X d t = μ m 1 X X m X (2)

where 𝑑𝑋/𝑑𝑡 is cell growth rate (cfu·mL-1·h-1), t is the fermentation time (h), X is the cell concentration (× 107 cfu·mL-1), Xm is the maximum cell concentration (× 107 cfu·mL-1), μ is the specific growth rate (h-1), μm is the maximum specific growth rate (h-1). When t = 0, X0 is the initial biomass concentration (× 107 cfu·mL-1), the integrated form of Equation 2 can use X = X0 gives a sigmoidal variation of X as a function of t (Equation 3), which can represent both exponential- and stationary-phase.

X = X 0 e b t 1 1 e b t X 0 X m (3)

The Luedeking-Piret equation is widely used to characterize the relationship between product formation and microbial growth (Dhagat & Jujjavarapu, 2015Dhagat, S., & Jujjavarapu, S. E. (2015). Isolation of a novel thermophilic bacterium capable of producing highyield bioemulsifier and its kinetic modelling aspects along with proposed metabolic pathway. Brazilian Journal of Microbiology, 51, 135-143.). GSBa-1 rennet also uses the Luedeking-Piret equation to describe its synthesis kinetics. The Luedeking-Piret equation is (Equation 4):

d P d t = α d X d t + β X (4)

where P is rennet activity (SU·mL-1), 𝛼 and 𝛽 are rennet synthesis kinetic parameters, 𝛼 is the product synthesis coefficient coupled with bacterial growth rate, 𝛽 is coupled with bacterial concentration Product synthesis constant. When t = 0, P0 is the initial rennet activity (SU·mL-1), the integrated form of Equation 3 can be changed to Equation 5.

P = P 0 α X + α X 0 X m e μ m t X m X 0 + X 0 e μ m t + B X m μ m ln X m X 0 + X 0 e μ m t X m (5)

B. amyloliquefaciens GSBa-1 consumes the substrate, which is used for bacterial growth, maintenance of bacterial life activities and rennet formation. The Luedeking-Piret modified equation was used to describe the substrate consumption process of the rennet produced by GSBa-1. The Luedeking-Piret modified equation is (Equation 6):

d S d t = dX ω dt + γ X (6)

where S is the substrate concentration (g·L-1), 𝜔 is the yield constant of the substrate for bacterial growth, γ is the yield constant of the formation of rennet.

When t = 0, S0 is the initial substrate concentration (g/L), the integrated form of Equation 6 can be changed to Equation 7.

S = S 0 + X 0 ω X 0 X m e μ m t ω X m X 0 + X 0 e μ m t + γ X m μ m ln X m X 0 + X 0 e μ m t X m (7)

These model parameters were estimated by Origin software (2018 version) and used for nonlinear fitting.

2.7 Milk-clotting activity assay

Milk-clotting activity in the fermentation broths was analyzed according to Arima's (Yu et al., 2014Yu, J., Tamura, G., & Arima, K. (2014). Milk-clotting enzyme from microorganisms: part VIII. Active center of milk-clotting enzyme (Mucor-rennin) isolated from Mucor pusillus Lindt. Agricultural and Biological Chemistry, 35(8), 1194-1199.) method. Two milliliters of 0.01 mol/L CaCl2 solution containing 10% skimmed milk was accurately measured and kept it at 35 °C for 5 min, and two hundred microliters of fermented crude enzyme solution was added and mixed quickly. The time between addition and appearance of clots was recorded. The formula for calculating milk-clotting activity is as follows (Equation 8):

MCA = 2400 T × 5 0.5 × D (8)

where MCA is the rennet activity, SU·mL-1; T is the coagulation time, s; D is the dilution multiple of the crude enzyme solution.

3 Results and discussion

3.1 Evaluation of the most significant factors affecting rennet activity

The optimization of key nutrients in fermentation medium for rennet production by B.amyloliquefaciens GSBa-1 using PBD was carried out in submerged cultures. Plackett-Burman design was carried out using MINITAB 17.0 software (Minitab, LLC, State College PA, USA) . The number of tests selected was N = 12, with four parameters X1, X2, X3, and X4 representing the maltose, corn steep liquor, beef extract, and Tween 80, respectively (Table S1). It can be seen that the obtained milk-clotting activities of the rennet ranged from 375 SU/mL to 437.125 SU/mL, which reflected the importance of medium optimization to enhancing milk-clotting activity.

The results of statistical analysis on the main effect were given in Table S3. The significance of each factor was assessed using the P-value at the significance level of 0.05. Pareto chart displayed that the significance of the five variables in terms of their effects on the milk-clotting activity. It was ranked as: B (Maltose) > A (Corn steep liquor) > D (Tween 80) > C (Beef extract) (Figure 1). Bars extending beyond this line correspond to statistically significant effects at a confidence level of 95% (Mayerhoff et al., 2006Mayerhoff, Z. D. V. L., Roberto, I. C., & Franco, T. T. (2006). Response surface methodology as an approach to determine the optimal activities of xylose reductase and xylitol dehydrogenase enzymes from Candida Mogii. Applied Microbiology and Biotechnology, 70(6), 761-767. http://dx.doi.org/10.1007/s00253-005-0304-3. PMid:16505992.
http://dx.doi.org/10.1007/s00253-005-030...
). Accordingly, the three variables were further considered for final optimization using BBD.

Figure 1
Pareto chart for the effect of medium nutrient components according to PBD.

3.2 RSM optimization for enhancing B. amyloliquefaciens GSBa-1 rennet activity

Based on the PBD and single factor test results, the BBD was used for optimization with maltose and corn syrup, maltose and Tween 80, corn syrup and Tween 80 as variables, and the following regression equation prediction model is obtained with quadratic multiple regression simulation after integration (Equation 9):

Y = 602.33 + 45.72 A + 26.16 B + 31.75 C 30.58 A B 24.87 A C + 10.19 B C 59.04 A 2 40.6 B 2 49.72 C 2 (9)

where Y is milk-clotting activity (SU/mL), A is maltose concentration (g/L), B is corn steep liquor concentration (g/L), and C is Tween 80 concentration(g/L), respectively.

The results were analyzed by variance analysis, and statistical significance was tested by F value. It has been shown from the obtained result that the model was highly significant, as suggested by the calculated F value (7.08) and low probability value (< 0.0001). The lack of fit F value of 0.3112 indicated that this value is not significant to the pure error. The coefficient of determination (R2) was calculated as 0.9758 for enhancing milk-clotting activity, representing that the statistical model can explain the adequate variability in the response (Table S4).

Furthermore, three-dimensional response surface plots showed the interactive effects of the three significant fermentation factors on the milk-clotting activity (Figure 2). It indicated that the maximum milk-clotting activity should occur with medium levels of maltose, corn syrup and tween 80. We found that the interactions between maltose and corn syrup (AB), and maltose and tween 80 (AC) were elliptical, indicating that a substantial correlation between each two variables and their influence on improving milk-clotting activity. The result was same as the p-value in Table S4. But The corn syrup and tween 80 (BC) contour lines were round, indicating a lack of significance. With the help of numerical optimization, predicted maximum milk-clotting activity (616.657 SU/mL) could be obtained with 7.258 g/L maltose, 8.15 g/L corn syrup and 1.97 g/L Tween 80. The milk-clotting activity by GSBa-1 is higher when compared with the thermophilic fungus N31 rennet (60.5 U/mL) and the fungal Quambalaria cyanescens rennet (117 SU/mL), but lower than a novel BL312 milk-clotting enzyme (MCE) (865 ± 20 SU/mL) exhibited high-level expression by optimizing induction conditions in recombinant Escherichia. coli harboring pET24a-proMCE (Almeida et al., 2015Almeida, C. M., Gomes, D., Faro, C., & Simões, I. (2015). Engineering a cardosin B-derived rennet for sheep and goat cheese manufacture. Applied Microbiology and Biotechnology, 99(1), 269-281. http://dx.doi.org/10.1007/s00253-014-5902-5. PMid:24986551.
http://dx.doi.org/10.1007/s00253-014-590...
; Wang et al., 2014Wang, C., Rao, H., Zhou, Z., & Zhang, Y. (2014). Mutagenic breeding of Quambalaria cyanescens strain for the production of milk-clotting enzymes. Pakistan Journal of Pharmaceutical Sciences, 27(Suppl. 3), 649-657. PMid:24816696.; Zhang et al., 2021Zhang, X., Tong, Y., Wang, J., Lyu, X., & Yang, R. (2021). Screening of a Bacillus subtilis strain producing both nattokinase and milk-clotting enzyme and its application in fermented milk with thrombolytic activity. Journal of Dairy Science, 104(9), 9437-9449. http://dx.doi.org/10.3168/jds.2020-19756. PMid:34218912.
http://dx.doi.org/10.3168/jds.2020-19756...
). What is more, bacterial MCEs have more potential than other MCEs due to their greater biochemical diversity and easier genetic modification. Therefore, the attempts should be made to increase the MCA/PA ratios of GSBa-1 MCEs by mutation or to study the changes in cleavage sites and functional peptides after genetic modification in the future.

Figure 2
Three-dimensional response surface plots showing the influence of variable pairs on rennet activity: A-maltose and corn syrup, B-maltose and Tween 80, C-corn steep liquor and Tween 80.

3.3 Optimization of fermentation condition in bioreactor

Bioreactors are practical tools that are used for economical, time-conserving and large-scale production of biomass from cell cultivation (Manjarres-Pinzón et al., 2022Manjarres-Pinzón, K., Mendoza-Meza, D., Arias-Zabala, M., Correa-Londoño, G., & Rodriguez-Sandoval, E. (2022). Effects of agitation rate and dissolved oxygen on xylose reductase activity during xylitol production at bioreactor scale. Food Science and Technology, 42, e04221. http://dx.doi.org/10.1590/fst.04221.
http://dx.doi.org/10.1590/fst.04221...
). Different constant pH, stirrer speed and fermentation temperature on the activity of GSBa-1 rennet were investigated in a fermenter with 3 L effective volume (Figure 3).

Figure 3
Optimization of fermentation parameters in bioreactor: A-pH, B-rotating speed, C-temperature. The bars represented by different letters are significantly different (p < 0.05).

The pH affects product growth and yield in microbial deep fermentation (Anggela et al., 2022Anggela, Harmayani, E., Setyaningsih, W., & Wichienchot, S. (2022). Prebiotic effect of porang oligo-glucomannan using fecal batch culture fermentation. Food Science and Technology, 42, e06321. http://dx.doi.org/10.1590/fst.06321.
http://dx.doi.org/10.1590/fst.06321...
; Kumar et al., 2010Kumar, A., Grover, S., Sharma, J., & Batish, V. K. (2010). Chymosin and other milk coagulants: sources and biotechnological interventions. Critical Reviews in Biotechnology, 30(4), 243-258. http://dx.doi.org/10.3109/07388551.2010.483459. PMid:20524840.
http://dx.doi.org/10.3109/07388551.2010....
). Since the acidity of the fermentation broth was constantly changing during the fermentation process, the method of automatic feeding was adopted in this experiment, and sodium hydroxide and hydrochloric acid were selected to control the pH to maintain a fixed value. Maximum milk-clotting activity of 525.3 SU/mL was obtained with the optimal pH around 7 (Figure 3A). In the acidic environment with a pH of 5.5, the milk-clotting activity of this broth decreased significantly, which indicated that the pH could significantly affect the rennet activity.

The shaker speed was related to the dissolved oxygen content in the culture solution (Ducros et al., 2009Ducros, E., Ferrari, M., Pellegrino, M., Raspanti, C., & Bogni, C. (2009). Effect of aeration and agitation on the protease production by Staphylococcus aureus mutant RC128 in a stirred tank bioreactor. Bioprocess and Biosystems Engineering, 32(1), 143-148. http://dx.doi.org/10.1007/s00449-008-0233-5. PMid:18491147.
http://dx.doi.org/10.1007/s00449-008-023...
). Figure 3B showed that the maximum milk-clotting activity of B. amyloliquefaciens GSBa-1 was enhanced when a faster stirring rate (200 rpm) was applied as compared to the cell density under the stirring rate of 160 rpm. But by further increasing the stirring speed, the milk-clotting activity was decreased.

Milk-clotting activity was also significantly regulated by different temperature. Maximum milk-clotting activity of 540.98 SU/mL was obtained with the optimal temperature around 45 °C (Figure 3C), which is often the case that most strains share the same trend, despite that the optimal culture temperatures differ. Rennet bioprocess based on B. amyloliquefaciens GSBa-1 was successfully scaled up based on rennet activity as the scale-up criterion from 0.25-L erlenmeyer flask to 5-L bioreactor.

3.4 Fermentation characteristics of high-density culture of B. amyloliquefaciens GSBa-1 to produce rennet

The changes of OD value and MCA during the fermentation of GSBa-1 were observed (Figure 4). The growth rate of GSBa-1 increased rapidly after 2 h, and the logarithmic growth phase began; and after 12 h, it entered the stable phase. At this time, the OD value of the crude enzyme solution was 1.6024, and the milk-clotting activity was 501.33 SU/mL, which increased by 12% than before. At 24 h, the accumulation of metabolites of the strain reach the highest, and the activity of rennet reached the highest value, which was 626.3 SU/mL.

Figure 4
Changes of MCA and OD600 during the process of culturing GSBa-1 to produce rennet in a fermenter.

3.5 Fermentation kinetic parameter fitting

To be of value for industrial biotechnology, mathematical models should be able to assist in the rational design of cell factory properties or in the production processes in which they are utilized. Kinetic models are particularly suitable towards this end because they are capable of representing the complex biochemistry of cells in a more complete way compared to most other types of models. They can, at least in principle, be used to in detail understand, predict, and evaluate the effects of adding, removing, or modifying molecular components of a cell factory and for supporting the design of the bioreactor or fermentation process. The initial values of biomass concentration, milk-clotting activity and substrate concentration were derived from the experimental results with mean values of 130.54, 0 and 30.608, respectively. And other unknown parameters in Equations 4, 6 and 8 were estimated by fitting the experimental data obtained from the batch fermentation.

Table S5 summarized the fitting curves of cell growth, milk-clotting activity and substrate consumption kinetics can be observed, and the R2 values of the model equations were 0.9927, 0.9168 and 0.9846, respectively (Figure 5). The results showed that the test value and the predicted value of the model have a high degree of fit, and the selected model could well represent the changes of various substances during the batch fermentation of B. amyloliquefaciens GSBa-1.

Figure 5
B. amyloliquefaciens GSBa-1 kinetic model fitting curve: A-GSBa-1 growth kinetic model, B-milk-clotting activity synthesis kinetic model, C-substrate consumption kinetic model.

4 Conclusion

Under the optimized medium and culture conditions, the GSBa-1 milk-clotting enzyme achieved the scale-up criterion from 0.25-L erlenmeyer flask to 5-L bioreactor successfully, and the milk-clotting activity of GSBa-1 rennet was increased to 626.3 SU/mL, which was an increase of 45.13% compared to that before optimization. What is more, a kinetic simulation model has been developed to predict accurately the dynamic behavior of metabolites in rennet fermentation by B. amyloliquefaciens GSBa-1 in the optimized medium. Therefore, these findings were valuable for the development of effective fermentation methods for the production of GSBa-1 rennet, and laid the foundation for the application of GSBa-1 rennet in the industrial production of cheese.

Acknowledgements

The National Key Research and Development Program of China (2018YFC1604302), National Natural Science Foundation of China (Grant no. 31871823), and Beijing Talent Cultivation Quality Construction—First-class Professional Construction (Municipal Level)—Food Science and Engineering (PXM2019_014213_000010) are kindly acknowledged.

  • Practical Application: The rennet produced by B. amyloliquefaciens GSBa-1 is a promising alternative to calf rennet in processing of cheese.

References

  • Abdel-Fattah, Y. R., Saeed, H. M., Gohar, Y. M., & El-Baz, M. A. (2005). Improved production of Pseudomonas aeruginosa uricase by optimization of process parameters through statistical experimental designs. Process Biochemistry, 40(5), 1707-1714. http://dx.doi.org/10.1016/j.procbio.2004.06.048
    » http://dx.doi.org/10.1016/j.procbio.2004.06.048
  • Almeida, C. M., Gomes, D., Faro, C., & Simões, I. (2015). Engineering a cardosin B-derived rennet for sheep and goat cheese manufacture. Applied Microbiology and Biotechnology, 99(1), 269-281. http://dx.doi.org/10.1007/s00253-014-5902-5 PMid:24986551.
    » http://dx.doi.org/10.1007/s00253-014-5902-5
  • Anggela, Harmayani, E., Setyaningsih, W., & Wichienchot, S. (2022). Prebiotic effect of porang oligo-glucomannan using fecal batch culture fermentation. Food Science and Technology, 42, e06321. http://dx.doi.org/10.1590/fst.06321
    » http://dx.doi.org/10.1590/fst.06321
  • Buehler, E. A., & Mesbah, A. (2016). Kinetic study of acetone-butanol-ethanol fermentation in continuous culture. PLoS One, 11(8), e0158243. http://dx.doi.org/10.1371/journal.pone.0158243 PMid:27486663.
    » http://dx.doi.org/10.1371/journal.pone.0158243
  • Chen, Y.-C., Huang, S.-D., Tu, J.-H., Yu, J.-S., Nurlatifah, A. O., Chiu, W.-C., Su, Y.-H., Chang, H.-L., Putri, D. A., & Cheng, H.-L. (2020). Exopolysaccharides of Bacillus amyloliquefaciens modulate glycemic level in mice and promote glucose uptake of cells through the activation of Akt. International Journal of Biological Macromolecules, 146, 202-211. http://dx.doi.org/10.1016/j.ijbiomac.2019.12.217 PMid:31887379.
    » http://dx.doi.org/10.1016/j.ijbiomac.2019.12.217
  • Dhagat, S., & Jujjavarapu, S. E. (2015). Isolation of a novel thermophilic bacterium capable of producing highyield bioemulsifier and its kinetic modelling aspects along with proposed metabolic pathway. Brazilian Journal of Microbiology, 51, 135-143.
  • Ducros, E., Ferrari, M., Pellegrino, M., Raspanti, C., & Bogni, C. (2009). Effect of aeration and agitation on the protease production by Staphylococcus aureus mutant RC128 in a stirred tank bioreactor. Bioprocess and Biosystems Engineering, 32(1), 143-148. http://dx.doi.org/10.1007/s00449-008-0233-5 PMid:18491147.
    » http://dx.doi.org/10.1007/s00449-008-0233-5
  • Elsayed, E. A., & Abdelwahed, N. A. (2020). Medium optimization by response surface methodology for improved cholesterol oxidase production by a newly isolated streptomyces rochei NAM-19 strain. BioMed Research International, 2020, 1870807. http://dx.doi.org/10.1155/2020/1870807 PMid:32337225.
    » http://dx.doi.org/10.1155/2020/1870807
  • Freudenberg, R. A., Baier, T., Einhaus, A., Wobbe, L., & Kruse, O. (2021). High cell density cultivation enables efficient and sustainable recombinant polyamine production in the microalga Chlamydomonas reinhardtii. Bioresource Technology, 323, 124542. http://dx.doi.org/10.1016/j.biortech.2020.124542 PMid:33385626.
    » http://dx.doi.org/10.1016/j.biortech.2020.124542
  • Ibrahim, A. A., El-Housseiny, G. S., Aboshanab, K. M., Yassien, M. A., & Hassouna, N. A. (2019). Paromomycin production from Streptomyces rimosus NRRL 2455: statistical optimization and new synergistic antibiotic combinations against multidrug resistant pathogens. BMC Microbiology, 19(1), 18. http://dx.doi.org/10.1186/s12866-019-1390-1 PMid:30658584.
    » http://dx.doi.org/10.1186/s12866-019-1390-1
  • Jiang, J., Zu, Y., Li, X., Meng, Q., & Long, X. (2020). Recent progress towards industrial rhamnolipids fermentation: process optimization and foam control. Bioresource Technology, 298, 122394. http://dx.doi.org/10.1016/j.biortech.2019.122394 PMid:31757615.
    » http://dx.doi.org/10.1016/j.biortech.2019.122394
  • Ju, Y., Son, K.-H., Jin, C., Hwang, B. S., Park, D.-J., & Kim, C.-J. (2018). Statistical optimization of culture medium for improved production of antimicrobial compound by Streptomyces rimosus AG-P1441. Food Science and Biotechnology, 27(2), 581-590. PMid:30263783.
  • Jung, S., & Lee, J.-H. (2020). Characterization of transcriptional response of Lactobacillus plantarum under acidic conditions provides insight into bacterial adaptation in fermentative environments. Scientific Reports, 10(1), 19203. http://dx.doi.org/10.1038/s41598-020-76171-6 PMid:33154427.
    » http://dx.doi.org/10.1038/s41598-020-76171-6
  • Kleman, G. L., & Strohl, W. R. (1994). Developments in high cell density and high productivity microbial fermentation. Current Opinion in Biotechnology, 5(2), 180-186. http://dx.doi.org/10.1016/S0958-1669(05)80033-3 PMid:7764797.
    » http://dx.doi.org/10.1016/S0958-1669(05)80033-3
  • Kumar, A., Grover, S., Sharma, J., & Batish, V. K. (2010). Chymosin and other milk coagulants: sources and biotechnological interventions. Critical Reviews in Biotechnology, 30(4), 243-258. http://dx.doi.org/10.3109/07388551.2010.483459 PMid:20524840.
    » http://dx.doi.org/10.3109/07388551.2010.483459
  • Liu, W.-C., Inwood, S., Gong, T., Sharma, A., Yu, L.-Y., & Zhu, P. (2019). Fed-batch high-cell-density fermentation strategies for Pichia pastoris growth and production. Critical Reviews in Biotechnology, 39(2), 258-271. http://dx.doi.org/10.1080/07388551.2018.1554620 PMid:30599783.
    » http://dx.doi.org/10.1080/07388551.2018.1554620
  • Manjarres-Pinzón, K., Mendoza-Meza, D., Arias-Zabala, M., Correa-Londoño, G., & Rodriguez-Sandoval, E. (2022). Effects of agitation rate and dissolved oxygen on xylose reductase activity during xylitol production at bioreactor scale. Food Science and Technology, 42, e04221. http://dx.doi.org/10.1590/fst.04221
    » http://dx.doi.org/10.1590/fst.04221
  • Martinez-Burgos, W. J., Sydney, E. B., Paula, D. R., Medeiros, A. B. P., Carvalho, J. C., Molina, D., & Soccol, C. R. (2021). Hydrogen production by dark fermentation using a new low-cost culture medium composed of corn steep liquor and cassava processing water: process optimization and scale-up. Bioresource Technology, 320(Pt B), 124370. http://dx.doi.org/10.1016/j.biortech.2020.124370 PMid:33220544.
    » http://dx.doi.org/10.1016/j.biortech.2020.124370
  • Mayerhoff, Z. D. V. L., Roberto, I. C., & Franco, T. T. (2006). Response surface methodology as an approach to determine the optimal activities of xylose reductase and xylitol dehydrogenase enzymes from Candida Mogii. Applied Microbiology and Biotechnology, 70(6), 761-767. http://dx.doi.org/10.1007/s00253-005-0304-3 PMid:16505992.
    » http://dx.doi.org/10.1007/s00253-005-0304-3
  • Mota, M. J., Lopes, R. P., Koubaa, M., Roohinejad, S., Barba, F. J., Delgadillo, I., & Saraiva, J. A. (2018). Fermentation at non-conventional conditions in food- and bio-sciences by the application of advanced processing technologies. Critical Reviews in Biotechnology, 38(1), 122-140. http://dx.doi.org/10.1080/07388551.2017.1312272 PMid:28423948.
    » http://dx.doi.org/10.1080/07388551.2017.1312272
  • Olivares-Marin, I. K., González-Hernández, J. C., Regalado-Gonzalez, C., & Madrigal-Perez, L. A. (2018). Saccharomyces cerevisiae exponential growth kinetics in batch culture to analyze respiratory and fermentative metabolism. Journal of Visualized Experiments, 139, 58192. PMid:30320748.
  • Phukoetphim, N., Salakkam, A., Laopaiboon, P., & Laopaiboon, L. (2017). Kinetic models for batch ethanol production from sweet sorghum juice under normal and high gravity fermentations: logistic and modified Gompertz models. Journal of Biotechnology, 243, 69-75. http://dx.doi.org/10.1016/j.jbiotec.2016.12.012 PMid:27988216.
    » http://dx.doi.org/10.1016/j.jbiotec.2016.12.012
  • Plackett, R. L., & Burman, J. P. (1946). The design of optimum multifactorial experiments. Biometrika, 33(4), 305-325. http://dx.doi.org/10.1093/biomet/33.4.305
    » http://dx.doi.org/10.1093/biomet/33.4.305
  • Procentese, A., Raganati, F., Olivieri, G., Russo, M. E., Salatino, P., & Marzocchella, A. (2015). Continuous lactose fermentation by Clostridium acetobutylicum-assessment of solventogenic kinetics. Bioresource Technology, 180, 330-337. http://dx.doi.org/10.1016/j.biortech.2015.01.008 PMid:25621726.
    » http://dx.doi.org/10.1016/j.biortech.2015.01.008
  • Rajeswari, P., Jose, P. A., Amiya, R. A., & Jebakumar, S. R. D. (2015). Characterization of saltern based Streptomyces sp. and statistical media optimization for its improved antibacterial activity. Frontiers in Microbiology, 5, 753. PMid:25653640.
  • Saavedra, S., Alejandro-Paredes, L., Flores-Santos, J. C., Flores-Fernández, C. N., Arellano-García, H., & Zavaleta, A. I. (2021). Optimization of lactic acid production by Lactobacillus plantarum strain Hui1 in a medium containing sugar cane molasses. Food Science and Technology, 39(1), 98-107.
  • Shahzad, R., Shehzad, A., Bilal, S., & Lee, I.-J. (2020). Bacillus amyloliquefaciens RWL-1 as a new potential strain for augmenting biochemical and nutritional composition of fermented soybean. Molecules, 25(10), 2346. http://dx.doi.org/10.3390/molecules25102346 PMid:32443519.
    » http://dx.doi.org/10.3390/molecules25102346
  • Singh, A. K., Mehta, G., & Chhatpar, H. S. (2009). Optimization of medium constituents for improved chitinase production by Paenibacillus sp. D1 using statistical approach. Letters in Applied Microbiology, 49(6), 708-714. http://dx.doi.org/10.1111/j.1472-765X.2009.02731.x PMid:19780958.
    » http://dx.doi.org/10.1111/j.1472-765X.2009.02731.x
  • Tashiro, T., & Yoshimura, F. (2019). A neo-logistic model for the growth of bacteria. Physica A, 525(03), 199-215. http://dx.doi.org/10.1016/j.physa.2019.03.049
    » http://dx.doi.org/10.1016/j.physa.2019.03.049
  • Teng, J., Zhao, X., Yang, Y., Zhang, J., Zhao, A., Jiang, Y., Li, L., Zheng, Z., & Yang, Z. (2017). Isolation and identification of microbial strains producing rennet from Jiuqu, a traditional Chinese fermentation starter. Shipin Kexue, 38(16), 23-28.
  • Teng, J.-W., Liu, Z.-M., Su, M.-Y., Gao, H.-Y., Zheng, Y.-R., Li, L., Zheng, Z., Zhang, J., & Yang, Z.-N. (2019). Application of rennet from Bacillus amyloliquefaciens GSBa-1 in mozzarella cheese. Shipin Kexue, 40, 85-92.
  • Thilakavathi, M., Basak, T., & Panda, T. (2007). Modeling of enzyme production kinetics. Applied Microbiology and Biotechnology, 73(5), 991-1007. http://dx.doi.org/10.1007/s00253-006-0667-0 PMid:17338073.
    » http://dx.doi.org/10.1007/s00253-006-0667-0
  • Tran, T. D., Cid, C., Hnasko, R., Gorski, L., & McGarvey, J. A. (2020). Bacillus amyloliquefaciens ALB65 inhibits the growth of Listeria monocytogenes on cantaloupe melons. Applied and Environmental Microbiology, 87(1), e01920-e01926. http://dx.doi.org/10.1128/AEM.01926-20 PMid:33097500.
    » http://dx.doi.org/10.1128/AEM.01926-20
  • Wang, C., Rao, H., Zhou, Z., & Zhang, Y. (2014). Mutagenic breeding of Quambalaria cyanescens strain for the production of milk-clotting enzymes. Pakistan Journal of Pharmaceutical Sciences, 27(Suppl. 3), 649-657. PMid:24816696.
  • Wang, T., Lu, Y., Yan, H., Li, X., Wang, X., Shan, Y., Yi, Y., Liu, B., Zhou, Y., & Lü, X. (2020). Fermentation optimization and kinetic model for high cell density culture of a probiotic microorganism: Lactobacillus rhamnosus LS-8. Bioprocess and Biosystems Engineering, 43(3), 515-528. http://dx.doi.org/10.1007/s00449-019-02246-y PMid:31712884.
    » http://dx.doi.org/10.1007/s00449-019-02246-y
  • Wang, W., Wu, X., Tan, J., Zhu, L., Mou, Y., Zhang, D., & Gao, J. (2019). Using response surface methodology optimize culture conditions for human lactoferrin production in desert Chlorella. Protein Expression and Purification, 155, 130-135. http://dx.doi.org/10.1016/j.pep.2018.11.004 PMid:30508587.
    » http://dx.doi.org/10.1016/j.pep.2018.11.004
  • Wei, G., Yang, X., Gan, T., Zhou, W., Lin, J., & Wei, D. (2009). High cell density fermentation of Gluconobacter oxydans DSM 2003 for glycolic acid production. Journal of Industrial Microbiology & Biotechnology, 36(8), 1029-1034. http://dx.doi.org/10.1007/s10295-009-0584-1 PMid:19434434.
    » http://dx.doi.org/10.1007/s10295-009-0584-1
  • WoldemariamYohannes, K., Wan, Z., Yu, Q., Li, H., Wei, X., Liu, Y., Wang, J., & Sun, B. (2020). Prebiotic, probiotic, antimicrobial, and functional food applications of Bacillus amyloliquefaciens. Journal of Agricultural and Food Chemistry, 68(50), 14709-14727. http://dx.doi.org/10.1021/acs.jafc.0c06396 PMid:33280382.
    » http://dx.doi.org/10.1021/acs.jafc.0c06396
  • Yang, X., Guo, P., Li, M., Li, H., Hu, Z., Liu, X., & Zhang, Q. (2020). Optimization of culture conditions for amoxicillin degrading bacteria screened from pig manure. International Journal of Environmental Research and Public Health, 17(6), 1973. http://dx.doi.org/10.3390/ijerph17061973 PMid:32192171.
    » http://dx.doi.org/10.3390/ijerph17061973
  • Yu, J., Tamura, G., & Arima, K. (2014). Milk-clotting enzyme from microorganisms: part VIII. Active center of milk-clotting enzyme (Mucor-rennin) isolated from Mucor pusillus Lindt. Agricultural and Biological Chemistry, 35(8), 1194-1199.
  • Zhang, X., Tong, Y., Wang, J., Lyu, X., & Yang, R. (2021). Screening of a Bacillus subtilis strain producing both nattokinase and milk-clotting enzyme and its application in fermented milk with thrombolytic activity. Journal of Dairy Science, 104(9), 9437-9449. http://dx.doi.org/10.3168/jds.2020-19756 PMid:34218912.
    » http://dx.doi.org/10.3168/jds.2020-19756
  • Zhang, Y., Xu, S., Cui, M., & Wong, J. W. C. (2019). Effects of different thermal pretreatments on the biodegradability and bioaccessibility of sewage sludge. Waste Management, 94, 68-76. http://dx.doi.org/10.1016/j.wasman.2019.05.047 PMid:31279397.
    » http://dx.doi.org/10.1016/j.wasman.2019.05.047
  • Zhao, D., Liu, L., Jiang, J., Guo, S., Ping, W., & Ge, J. (2020a). The response surface optimization of exopolysaccharide produced by Weissella confusa XG-3 and its rheological property. Preparative Biochemistry & Biotechnology, 50(10), 1014-1022. http://dx.doi.org/10.1080/10826068.2020.1780609 PMid:32589090.
    » http://dx.doi.org/10.1080/10826068.2020.1780609
  • Zhao, X., Cai, M., Yang, Z., Luo, T., Chen, C., Cao, Y., & Yang, Z. (2020b). Proteolysis and bioactivity of cheddar cheese as influenced by milk-clotting enzyme from Bacillus amyloliquefaciens GSBa-1. Shipin Kexue, 41, 103-111.
  • Zhou, Q., Liu, Y., & Yuan, W. (2020). Kinetic modeling of butyric acid effects on butanol fermentation by Clostridium saccharoperbutylacetonicum. New Biotechnology, 55, 118-126. http://dx.doi.org/10.1016/j.nbt.2019.10.004 PMid:31626983.
    » http://dx.doi.org/10.1016/j.nbt.2019.10.004

Publication Dates

  • Publication in this collection
    03 June 2022
  • Date of issue
    2022

History

  • Received
    22 Mar 2022
  • Accepted
    06 May 2022
Sociedade Brasileira de Ciência e Tecnologia de Alimentos Av. Brasil, 2880, Caixa Postal 271, 13001-970 Campinas SP - Brazil, Tel.: +55 19 3241.5793, Tel./Fax.: +55 19 3241.0527 - Campinas - SP - Brazil
E-mail: revista@sbcta.org.br