*Article* **Use of a Sequential Fermentation Method for the Production of** *Aspergillus tamarii* **URM4634 Protease and a Kinetic/ Thermodynamic Study of the Enzyme**

**Rodrigo Lira de Oliveira <sup>1</sup> , Emiliana de Souza Claudino <sup>1</sup> , Attilio Converti 2,\* and Tatiana Souza Porto <sup>3</sup>**


**Abstract:** Microbial proteases are commonly produced by submerged (SmF) or solid-state fermentation (SSF), whose combination results in an unconventional method, called sequential fermentation (SF), which has already been used only to produce cellulolytic enzymes. In this context, the aim of the present study was the development of a novel SF method for protease production using wheat bran as a substrate. Moreover, the kinetic and thermodynamic parameters of azocasein hydrolysis were estimated, thus providing a greater understanding of the catalytic reaction. In SF, an approximately 9-fold increase in protease activity was observed compared to the conventional SmF method. Optimization of glucose concentration and medium volume by statistical means allowed us to achieve a maximum protease activity of 180.17 U mL−<sup>1</sup> . The obtained enzyme had an optimum pH and temperature of 7.0 and 50 ◦C, respectively. Kinetic and thermodynamic parameters highlighted that such a neutral protease is satisfactorily thermostable at 50 ◦C, a temperature commonly used in many applications in the food industry. The results obtained suggested not only that SF could be a promising alternative to produce proteases, but also that it could be adapted to produce several other enzymes.

**Keywords:** *Aspergillus*; protease; sequential fermentation; kinetics; thermodynamics

## **1. Introduction**

Proteases are enzymes that catalyze the hydrolysis of proteins resulting in their breakdown into peptides and amino acids. They represent an important group of industrial enzymes, accounting for almost 60% of the entire enzyme market and 40% of the global enzyme sales [1]. Proteases are in fact used in various industrial sectors, including those of detergents, leather goods, food, fine chemicals, pharmaceuticals, and cosmetics, and can be obtained from different sources such as plants, animal, and microorganisms, with the latter source contributing to nearly two-third of industrially used proteases. This is mainly due to the possibility of producing them quickly in large amounts by established fermentation processes that also allow for their regular supply [2,3]. Filamentous fungi have an extraordinary capability to secrete different types of proteases, with the main fungal strains for their industrial production belonging to the genera *Trichoderma* and *Aspergillus* [2].

Commonly, microbial proteases are produced by submerged fermentation (SmF), a process in which microorganisms grow in a liquid medium with high free water content, or by solid state fermentation (SSF), where they grow in an environment without free water or with a very low content of free water [4]. Both processes have their own advantages

**Citation:** de Oliveira, R.L.; de Souza Claudino, E.; Converti, A.; Porto, T.S. Use of a Sequential Fermentation Method for the Production of *Aspergillus tamarii* URM4634 Protease and a Kinetic/Thermodynamic Study of the Enzyme. *Catalysts* **2021**, *11*, 963. https://doi.org/10.3390/ catal11080963

Academic Editors: Chun-Yung Huang, Chwen-Jen Shieh, Cheng-Di Dong and Chia-Hung Kuo

Received: 27 June 2021 Accepted: 10 August 2021 Published: 11 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

and disadvantages depending on the operating conditions. SmF has advantages related to the presence of instrumentation and process control, and consequently it is the preferred process for industrial enzyme production. On the other hand, SSF can be particularly advantageous to cultivate filamentous fungi, because it simulates their natural habitat leading to higher enzyme productivity compared to SmF [5,6].

A combination of these two cultivation methods, called sequential fermentation (SF or SeqF), was proposed initially by Cunha et al. [7] and applied effectively to produce cellulolytic enzymes. It is characterized by a preculture preparation with an initial step of fungal growth under solid state, followed by a submerged step [8]. This fermentation technique, so far, has been poorly documented and used exclusively to produce lignocellulolytic enzymes using filamentous fungi such as *Aspergillus niger* and *Trichoderma reesei* [9,10] or the white rot fungus *Pleurotus ostreatus* [11]. Different configurations were validated for this approach, using shake flasks, stirred tanks and bubble column bioreactors, which allowed an increase in enzyme productivity compared to the conventional SmF [8]. Most studies have reported the use of sugarcane bagasse as a substrate, nonetheless other substrates such as wheat bran and soybean bran have also been used successfully [12].

Although SF is used for a specific category of enzymes [8], it could be adapted to produce different enzymes of industrial interest. Based on that background, sequential fermentation is used for the first time in this study to produce *Aspergillus tamarii* protease. Particularly, the aims of this study were: (a) to test the sequential fermentation method for protease production by an *Aspergillus tamarii* strain using wheat bran as a substrate, (b) to biochemically characterize the crude extract, and (c) to determine and investigate the kinetic and thermodynamic parameters related to both substrate hydrolysis and enzyme thermal denaturation.

#### **2. Results and Discussion**

#### *2.1. Protease Production by Sequential Fermentation*

Protease production by *Aspergillus tamarii* URM4634 in submerged fermentation (SmF) and sequential fermentation (SF) was evaluated after 48 h of cultivation for comparison purposes. Even though other fermentation times were evaluated (data not shown), 48 h was shown to be the most suitable for enzyme production in both fermentation processes. SF allowed an almost double enzyme production (39.08 U mL−<sup>1</sup> ) compared to SmF (20.0 U mL−<sup>1</sup> ), confirming the positive trend observed by several authors who used this approach to produce cellulolytic enzymes [8,10,12]. A possible explanation for this increase is related to the morphology of filamentous fungi, since in SF there was a predominance of dispersed filamentous mycelium, while in the conventional SmF early formation of fungal pellets took place [9]. However, the dispersed hyphal morphology results in a more viscous medium that can make mixing and aeration difficult [8]; furthermore, other factors should be considered for a complete understanding of SF, especially if intended for the production of proteolytic enzymes.

Glucose concentration, medium volume and inoculum size were selected as the independent variables for protease production by SF, because they are among those that mostly influence the conventional fermentation processes. Their effects on such a response were investigated through cultivations carried out according to a 2<sup>3</sup> -full factorial design whose results are listed in Table 1. The highest protease production (127.67 U mL−<sup>1</sup> ) was obtained in the run 6 performed using a glucose concentration of 50 g L−<sup>1</sup> , a medium volume of 15 mL g−<sup>1</sup> of substrate and an inoculum size of 10<sup>8</sup> spores g−<sup>1</sup> of substrate. According to the statistical analysis (Table 2), all the individual variables and their interactions were statistically significant. In particular, glucose concentration showed a positive effect probably because, as known, in fermentation processes the carbon source not only acts as the major constituent for the construction of cell material but is also used in the synthesis of polysaccharides and as an energy source [13]. However, it should be remembered that glucose can have a repressive effect on the synthesis of other enzymes such as pectinases [14]. A similar positive effect was exerted by the inoculum size, while the medium volume

had a negative influence on protease production. A larger inoculum size reduces the time that cells take to consume all the substrate and produce the enzyme of interest. However, beyond a certain limit, enzyme production can decrease due to the depletion of nutrients needed for biomass growth, which results in a decrease in metabolic activity [15,16]. To the best of our knowledge, this is the first report on the evaluation of the influence of medium volume on SF, and the trend observed at lower medium volumes may be related to the prevention of enzyme dilution in the culture medium itself.

**Table 1.** Experimental conditions and results of protease production by *Aspergillus tamarii* URM4634 in sequential fermentations using wheat bran as a substrate. Runs were carried out according to a 2 3 -full factorial design plus three central points.


**Table 2.** Effects of the independent variables and their interactions on protease production by *Aspergillus tamarii* URM4634 in sequential fermentations carried out according to the 2<sup>3</sup> -full factorial design shown in Table 1 using wheat bran as a substrate.


\* Statistically significant estimates at 95% confidence level (*p* < 0.05).

Interaction effects are interpreted differently from what is usually done with the individual variables. A positive interaction between two variables means that they have a synergistic effect, i.e., the levels of both must increase or decrease together to enhance the response. A negative interaction, on the other hand, means that the variables have an antagonistic effect, i.e., the level of one variable must increase and that of the other decrease to enhance the response. In the case of interaction among three variables, the negative sign means that one of them has an antagonistic effect with the others. More specifically, as shown by the results of Table 2, the medium volume had a negative effect, while the other variables a positive effect. The interaction had a negative sign, but the behavior of each variable involved in the interaction must be evaluated. The interaction among the three factors can be better seen in the cubic graph depicted in Figure 1.

factors can be better seen in the cubic graph depicted in Figure 1.

variable involved in the interaction must be evaluated. The interaction among the three

**Figure 1.** Cubic plot of the effects of glucose concentration, medium volume, and inoculum size on the production of *Aspergillus tamarii* URM4634 protease by sequential fermentation. **Figure 1.** Cubic plot of the effects of glucose concentration, medium volume, and inoculum size on the production of *Aspergillus tamarii* URM4634 protease by sequential fermentation.

The antagonist ternary interaction among the independent variables confirms the highest levels of glucose and inoculum size and the lowest one of medium volume as the best fermentation conditions. For optimization purposes and due to the difficulty of carrying out industrial production at higher spore concentrations, a 2<sup>2</sup> -CCRD was performed using only glucose concentration and medium volume as the independent variables, while keeping the inoculum size constant at the highest level (10<sup>8</sup> spores g−1), whose experimental conditions and results are listed in Table 3. The highest protease activity (180.17 U mL−1) was observed in the run 4 carried out at a glucose concentration of 55 g The antagonist ternary interaction among the independent variables confirms the highest levels of glucose and inoculum size and the lowest one of medium volume as the best fermentation conditions. For optimization purposes and due to the difficulty of carrying out industrial production at higher spore concentrations, a 2<sup>2</sup> -CCRD was performed using only glucose concentration and medium volume as the independent variables, while keeping the inoculum size constant at the highest level (10<sup>8</sup> spores g−<sup>1</sup> ), whose experimental conditions and results are listed in Table 3. The highest protease activity (180.17 U mL−<sup>1</sup> ) was observed in the run 4 carried out at a glucose concentration of 55 g L−<sup>1</sup> and a medium volume of 17.5 mL per g of substrate.

The analysis of variance (ANOVA) (Table 4) indicates that the quadratic model showed a significant lack of fit and an unsatisfactory adjustment to the experimental data (R<sup>2</sup> = 0.44), therefore, the response surface was not shown. Nonetheless, the experiments **Table 3.** Experimental conditions and results of protease production by *Aspergillus tamarii* URM4634 in sequential fermentations using wheat bran as a substrate. Optimization tests were carried out according to a central composite rotational design.


L−1 and a medium volume of 17.5 mL per g of substrate.

The analysis of variance (ANOVA) (Table 4) indicates that the quadratic model showed a significant lack of fit and an unsatisfactory adjustment to the experimental data (R<sup>2</sup> = 0.44), therefore, the response surface was not shown. Nonetheless, the experiments performed were effective in increasing protease production by as much as 4.6 times compared to the initial SF runs. Only the medium volume showed a significant linear negative effect (−8.05), and its interaction with glucose concentration, contrary to the statistical analysis of the 2<sup>3</sup> -full factorial design, was significant and synergistic (6.55) (results not shown), likely due to the different experimental regions investigated.

**Table 4.** Results of the analysis of variance applied to the results of the central composite rotational design used for protease production by *Aspergillus tamarii* URM4634 in sequential fermentations using wheat bran as a substrate.


GC = Glucose concentration; MV = Medium volume; L = Linear effect; Q = Quadratic effect.

#### *2.2. Effect of Temperature and pH on Protease Activity*

The rate of enzyme-catalyzed reactions can be influenced by various environmental factors through reversible or irreversible changes in protein structure, of which temperature and pH are the most important [17]. The effect of a temperature rise on enzyme activity is the result of two phenomena acting simultaneously but in opposite directions, i.e., an increase in the reaction rate and a progressive enzyme inactivation due to its denaturation. On the other hand, the enzyme activity is influenced by pH because the ionization state of amino acids constituting the enzyme structure depends on the concentration of hydrogen ions [17,18]. In the present study, the optimal temperature and pH of *A. tamarii* URM4634 protease were shown to be 50 ◦C (Figure 2A) and 7.0 (Figure 2B), respectively.

In a previous study by our research group, the protease obtained by SSF using the same fungal strain and substrate (wheat bran) showed different optimal conditions (40 ◦C and pH 8.0) [19], suggesting that the fermentation method may directly influence the characteristics of the biocatalyst. A similar optimum temperature was observed for proteases obtained from different *Aspergillus* species, such as *A. heteromorphus* [20], *A. niger* [21] and *A. oryzae* [22], in conventional fermentation processes (SSF and SmF).

Based on these and previous results, the *A. tamarii* URM4634 protease under investigation was characterized as a neutral protease, being active under neutral, weakly acidic or weakly alkaline conditions. Proteases belonging to this family are considered especially important for the food industry, mainly for baking and brewing processes, as they lead to less bitter food protein hydrolysates thanks to an intermediate reaction rate and are poorly sensitive to plant proteinase inhibitors [23,24]. Fungal neutral proteases, of which *A. oryzae* is the most important producer [25], are constituents of various enzyme preparations, and the search for new sources of proteases has aroused considerable interest. Therefore, the neutral protease from *A. tamarii* URM4634 obtained in this study by the SF method may be an interesting alternative for further studies involving applications in the food industry. A neutral protease from another *A. tamarii* strain was already produced by our research group by extractive fermentation [26], which provides confirmation of the potential of this species as a producer of neutral proteases.

**Figure 2.** Effects of temperature (**A**) and pH (**B**) on the activity of *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substrate. Different letters (a–f) indicate statistically significant differences (*p* < 0.05). **Figure 2.** Effects of temperature (**A**) and pH (**B**) on the activity of *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substrate. Different letters (a–f) indicate statistically significant differences (*p* < 0.05).

#### *2.3. Kinetic and Thermodynamic Parameters of Azocasein Hydrolysis*

Based on these and previous results, the *A. tamarii* URM4634 protease under investigation was characterized as a neutral protease, being active under neutral, weakly acidic or weakly alkaline conditions. Proteases belonging to this family are considered especially important for the food industry, mainly for baking and brewing processes, as they lead to less bitter food protein hydrolysates thanks to an intermediate reaction rate and are poorly sensitive to plant proteinase inhibitors [23,24]. Fungal neutral proteases, of which *A. oryzae* is the most important producer [25], are constituents of various enzyme preparations, and the search for new sources of proteases has aroused considerable interest. Therefore, the neutral protease from *A. tamarii* URM4634 obtained in this study by the SF method may be an interesting alternative for further studies involving applications in the food industry. A neutral protease from another *A. tamarii* strain was already produced by our research group by extractive fermentation [26], which provides confirmation of the potential of this species as a producer of neutral proteases. The crude extract obtained by SF showed a typical Michaelis-Menten kinetic profile, and the kinetic parameters of azocasein hydrolysis were estimated from the slope of Lineweaver-Burk plot (Figure 3A) with satisfactory correlation (R<sup>2</sup> = 0.993). As shown in Table 5, the enzyme exhibited values of the Michaelis constant (Km), maximum rate (Vmax) and turnover number (kcat) of 16.26 mg mL−<sup>1</sup> , 147.06 mg mL−<sup>1</sup> min−<sup>1</sup> and 195.37 s−<sup>1</sup> , respectively. Such a K<sup>m</sup> value is lower than that obtained in SSF (18.7 mg mL−<sup>1</sup> ) [19] using the same fungal strain and substrate, while that of Vmax is 5.15-fold higher (28.57 mg mL−<sup>1</sup> min−<sup>1</sup> ) [19]. These results indicate an increase in the affinity for substrate and the catalytic potential of protease produced by SF compared to that produced by SSF. The value of kcat (195.37 s−<sup>1</sup> ), which corresponds to the maximum number of substrate molecules converted into product per active site per unit of time, is almost 20% higher than that reported by Lee et al. [27] for purified *A. oryzae* protease obtained by SmF using casein as a substrate (163.5 s−<sup>1</sup> ).

*2.3. Kinetic and Thermodynamic Parameters of Azocasein Hydrolysis* The crude extract obtained by SF showed a typical Michaelis-Menten kinetic profile, and the kinetic parameters of azocasein hydrolysis were estimated from the slope of Lineweaver-Burk plot (Figure 3A) with satisfactory correlation (R<sup>2</sup> = 0.993). As shown in Table 5, the enzyme exhibited values of the Michaelis constant (Km), maximum rate (Vmax) and turnover number (kcat) of 16.26 mg mL−1, 147.06 mg mL−1 min−1 and 195.37 s−1, respectively. Such a K<sup>m</sup> value is lower than that obtained in SSF (18.7 mg mL−1) [19] using the same fungal strain and substrate, while that of Vmax is 5.15-fold higher (28.57 mg mL−1 min−1) From the slope of the straight line to the right of the Arrhenius-type plot (Figure 3B), an activation energy (E\*a) of 40.38 kJ mol−<sup>1</sup> was estimated for azocasein hydrolysis (Table 5). This value is lower than that reported for casein hydrolysis by *A. fumigatus* protease (62 kJ mol−<sup>1</sup> ) [28] and slightly higher than that of bread protein hydrolysis by *A. awamori* protease (36.8 kJ mol−<sup>1</sup> ) [29]. It is worth remembering, for comparative purposes, that a lower E\* <sup>a</sup> value indicates that less energy is required to form the activated complex of azocasein hydrolysis, so the reaction is favored. Another important parameter to evaluate the temperature influence on a catalyzed reaction is the temperature quotient (Q10), which indicates the rate increase resulting from a 10 ◦C temperature rise. At temperatures

s

s

between 30 and 50 ◦C, the mean value of Q<sup>10</sup> (1.17) was very close to those reported by de Castro et al. [30] for *A. niger* protease (1.20 to 1.28) in the temperature range of 35–55 ◦C. It has been stressed that deviations of Q<sup>10</sup> values from the typical range of enzyme reactions (1.0 to 2.0) are indicative of the involvement of factors other than temperature in the control of the reaction rate [31]. Therefore, it can be inferred that azocasein hydrolysis by *A. tamarii* protease was kinetically controlled by temperature within the entire temperature range investigated. [19]. These results indicate an increase in the affinity for substrate and the catalytic potential of protease produced by SF compared to that produced by SSF. The value of kcat (195.37 −1), which corresponds to the maximum number of substrate molecules converted into product per active site per unit of time, is almost 20% higher than that reported by Lee et al. [27] for purified *A. oryzae* protease obtained by SmF using casein as a substrate (163.5 −1).

**Figure 3.** (**A**) Lineweaver-Burk plot to calculate K<sup>m</sup> and Vmax of *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using azocasein as a substrate (R<sup>2</sup> = 0.993), and (**B**) Arrheniustype plot (on the right) to estimate the activation energy of the azocasein hydrolysis reaction (R<sup>2</sup> = 0.998). **Figure 3.** (**A**) Lineweaver-Burk plot to calculate K<sup>m</sup> and Vmax of *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using azocasein as a substrate (R<sup>2</sup> = 0.993), and (**B**) Arrhenius-type plot (on the right) to estimate the activation energy of the azocasein hydrolysis reaction (R<sup>2</sup> = 0.998).

**Table 5.** Kinetic and thermodynamic parameters of the azocasein hydrolysis reaction catalyzed by the *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substrate. 25 °C was selected as a reference temperature. **Table 5.** Kinetic and thermodynamic parameters of the azocasein hydrolysis reaction catalyzed by the *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substrate. 25 ◦C was selected as a reference temperature.


<sup>g</sup>ΔS\* (J K−1 mol−1) −73.94 <sup>a</sup>Michaelis constant; <sup>b</sup>Maximum rate; <sup>c</sup> Turnover number; <sup>d</sup>Activation energy; <sup>e</sup> Activation Gibbs <sup>a</sup> Michaelis constant; <sup>b</sup> Maximum rate; <sup>c</sup> Turnover number; <sup>d</sup> Activation energy; <sup>e</sup> Activation Gibbs free energy; <sup>f</sup> Activation enthalpy; <sup>g</sup> Activation entropy.

free energy; <sup>f</sup> Activation enthalpy; <sup>g</sup>Activation entropy.

**T (°C)**

**<sup>a</sup>k<sup>d</sup> (min−1)**

**R<sup>2</sup>**

50 0.0030 0.979 231.05 767.53

denaturation, and <sup>h</sup>Activation entropy of thermal denaturation.

The other thermodynamic parameters of azocasein hydrolysis, namely activation enthalpy (∆H\*), entropy (∆S\*) and Gibbs free energy (∆G\*), were also calculated by Equations (2) to (4) and listed in Table 5.

The relatively low ∆H\* value (37.90 kJ mol−<sup>1</sup> ) indicates that the formation of the transition state or activated enzyme-substrate complex occurred effectively; however, this value is higher than that reported by Fernandes et al. [20] for azocasein hydrolysis by *A. heteromorphus* protease (21.8 kJ mol−<sup>1</sup> ). As known, ∆S\* is correlated to the order degree of a reaction system; therefore, in enzyme-catalyzed reactions negative values such as that estimated in this study (−73.94 J K−<sup>1</sup> mol−<sup>1</sup> ) suggest that the structure of enzyme-substrate at transition state is more ordered than that of the enzyme-substrate complex. Similar qualitative behavior was observed for proteases [32] and other hydrolases such as β-fructofuranosidases [31], pectinases [33] and levansucrases [34]. As known, ∆G\* is the most suitable thermodynamic parameter to measure the feasibility and extent of a reaction, as the lower the ∆G\* value the more feasible the reaction, i.e., the conversion of the reactant into the product is more spontaneous [35]. The value of ∆G\* estimated for azocasein hydrolysis catalyzed by *A. tamarii* URM4634 protease (59.94 kJ mol−<sup>1</sup> ) is lower than those of casein and azocasein hydrolysis catalyzed by other fungal (*A. heteromorphus*: ∆G\* = 65.8 kJ mol−<sup>1</sup> [20]) and bacterial (*Bacillus stearothermophilus*: ∆G\* = 91.71 kJ mol−<sup>1</sup> [32]; *Bacillus cereus*: ∆G\* = 68.68–69.42 kJ mol−<sup>1</sup> [36]) proteases, which indicates that the conversion of the transition state of enzyme-substrate complex into the product is more spontaneous.

#### *2.4. Kinetic and Thermodynamic Parameters of Protease Thermal Denaturation*

The time-dependent thermal inactivation of biocatalysts (denaturation) can be a major problem in enzyme-catalyzed processes, and knowledge of the impact of this phenomenon is particularly useful before running any industrial-scale bioprocess. Then, to evaluate the protease thermostability, residual activity tests were performed in the temperature range of 50 to 80 ◦C, the results of which are shown in semi-log plots (Figure 4). *Catalysts* **2021**, *11*, x FOR PEER REVIEW 9 of 16

**Figure 4.** Semi-log plots of the coefficient of residual activity (A/A0) of *Aspergillus tamarii* URM4634 protease versus time used to estimate the first-order inactivation rate constant (kd) at different temperatures (°C): () 50; () 60; () 70; () 80. **Figure 4.** Semi-log plots of the coefficient of residual activity (A/A<sup>0</sup> ) of *Aspergillus tamarii* URM4634 protease versus time used to estimate the first-order inactivation rate constant (kd) at different temperatures (◦C): (•) 50; () 60; ( ) 70; () 80.

**Table 6.** Kinetic and thermodynamic parameters of thermal denaturation of crude *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substate. **<sup>b</sup> t1/2 (min) <sup>c</sup> D-Value (min) <sup>d</sup>Z-Value (°C) <sup>e</sup>E\*d (kJ mol−1) <sup>f</sup>ΔG\*<sup>d</sup> (kJ mol−1) <sup>g</sup>ΔH\*<sup>d</sup> (kJ mol−1) <sup>h</sup>ΔS\*<sup>d</sup> (J mol−1 K−1)** The slopes of the straight lines obtained at different temperatures were used to estimate the first-order inactivation rate constant (kd), whose values were estimated with a satisfactory correlation (0.969 ≤ R <sup>2</sup> <sup>≤</sup> 0.998) and are listed in Table <sup>6</sup> together with those of the other kinetic and thermodynamic parameters. One can see that k<sup>d</sup> progressively increased with temperature from 0.0030 to 0.0227 min−<sup>1</sup> , which means that denaturation

> At 50 °C, a temperature commonly used in several industrial processes, the crude protease obtained by SF exhibited much longer half-life (t1/2 = 231.05 min) and decimal reduction time (D-value = 767.53 min) than that produced by the same fungal strain and substrate by SSF, either in crude (t1/2: 96.3 min; D-value: 319.9 min) or purified (t1/2: 70 min; D-value: 232.6 min) form [37]. These results indicate that the protease produced by SF is more thermostable than that obtained by SSF; therefore, one of the next efforts will be an enzyme structural analysis that can provide some additional information for a better understanding of the thermostabilization mechanism under SF conditions. Moreover, the protease produced in the present study appears to be more thermostable than other *Aspergillus* proteases at same temperature (50 °C), such as those produced by *A. fumigatus* (t1/2: 96.3 min) [28], *A. niger* (t1/2: 135.91 min; D-value: 451.49 min) [30] and *A. heteromorphus*

> As proof of protease thermostability, a thermal resistance constant (Z-value) of 33.89 °C was estimated with good correlation (R<sup>2</sup> = 0.945) from the slope of the semi-log plot of D-value versus temperature (data not shown). In general, high Z-values, such as that obtained in the present study, indicate a greater sensitivity to the duration of the heat treatment than to the temperature increase [35]. The Z-value obtained is lower than that detected for the protease produced by SSF either in crude (44.2 °C) or purified (75.8 °C) form [37], which suggests a noticeable reduction in sensitivity to duration of thermal treatment

105.97 62.09 −135.76

60 0.0083 0.969 83.51 277.42 106.51 62.01 −133.57

80 0.0227 0.998 30.53 101.43 110.12 61.84 −136.71 <sup>a</sup> First-order inactivation rate constant; <sup>b</sup> Half-life; <sup>c</sup> Decimal reduction time; <sup>d</sup> Thermal resistance constant; <sup>e</sup> Activation energy of thermal denaturation; <sup>f</sup> Activation Gibbs free energy of thermal denaturation; <sup>g</sup> Activation enthalpy of thermal

(t1/2: 54.6 min; D-value: 181.3 min) [20].

when the enzyme is produced by SF.

became progressively more significant, likely due to the breaking of an increasing number of strong electrostatic bonds.

**Table 6.** Kinetic and thermodynamic parameters of thermal denaturation of crude *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substate.


<sup>a</sup> First-order inactivation rate constant; <sup>b</sup> Half-life; <sup>c</sup> Decimal reduction time; <sup>d</sup> Thermal resistance constant; <sup>e</sup> Activation energy of thermal denaturation; <sup>f</sup> Activation Gibbs free energy of thermal denaturation; <sup>g</sup> Activation enthalpy of thermal denaturation, and <sup>h</sup> Activation entropy of thermal denaturation.

> At 50 ◦C, a temperature commonly used in several industrial processes, the crude protease obtained by SF exhibited much longer half-life (t1/2 = 231.05 min) and decimal reduction time (D-value = 767.53 min) than that produced by the same fungal strain and substrate by SSF, either in crude (t1/2: 96.3 min; D-value: 319.9 min) or purified (t1/2: 70 min; D-value: 232.6 min) form [37]. These results indicate that the protease produced by SF is more thermostable than that obtained by SSF; therefore, one of the next efforts will be an enzyme structural analysis that can provide some additional information for a better understanding of the thermostabilization mechanism under SF conditions. Moreover, the protease produced in the present study appears to be more thermostable than other *Aspergillus* proteases at same temperature (50 ◦C), such as those produced by *A. fumigatus* (t1/2: 96.3 min) [28], *A. niger* (t1/2: 135.91 min; D-value: 451.49 min) [30] and *A. heteromorphus* (t1/2: 54.6 min; D-value: 181.3 min) [20].

> As proof of protease thermostability, a thermal resistance constant (Z-value) of 33.89 ◦C was estimated with good correlation (R<sup>2</sup> = 0.945) from the slope of the semi-log plot of Dvalue versus temperature (data not shown). In general, high Z-values, such as that obtained in the present study, indicate a greater sensitivity to the duration of the heat treatment than to the temperature increase [35]. The Z-value obtained is lower than that detected for the protease produced by SSF either in crude (44.2 ◦C) or purified (75.8 ◦C) form [37], which suggests a noticeable reduction in sensitivity to duration of thermal treatment when the enzyme is produced by SF.

> As for the thermodynamic parameters of protease thermal denaturation, an activation energy (E\*d) of 64.78 kJ mol−<sup>1</sup> was estimated from the Arrhenius-type plot illustrated in Figure 5 with satisfactory correlation (R<sup>2</sup> = 0.955). As known, E\* <sup>d</sup> is related to the activation enthalpy of thermal denaturation (∆H\*d), i.e., the total amount of thermal energy required to denature the enzyme through disruption of non-covalent bonds, including hydrophobic interactions [31], which varied between 61.84 and 62.09 kJ mol−<sup>1</sup> . Comparison of the above E\*<sup>d</sup> and ∆H\* <sup>d</sup> values with those of protease obtained by SSF either in crude (E\*<sup>d</sup> = 49.7 kJ mol−<sup>1</sup> ; ∆H\*<sup>d</sup> = 46.7–47.0 kJ mol−<sup>1</sup> ) or purified (E\*<sup>d</sup> = 28.8 kJ mol−<sup>1</sup> ; ∆H\*<sup>d</sup> = 25.9–26.1 kJ mol−<sup>1</sup> ) form [37] confirms the higher thermostability of the enzyme produced by SF. To the best of our knowledge, studies comparing the thermostability of enzymes produced by different fermentation processes are very scarce. Nonetheless, the work of Saqib et al. [38] can be cited, in which the authors reported, on the basis of kinetic and thermodynamic parameters, a greater thermostability of an *A. fumigatus* endoglucanase produced by SSF rather than by SmF. Therefore, further studies are needed to investigate possible structural differences in enzymes produced in different fermentation processes, including unconventional methods such as SF.

including unconventional methods such as SF.

drophobic interactions [31], which varied between 61.84 and 62.09 kJ mol−1

**Figure 5.** Arrhenius-type plot used to estimate the activation energy of thermal denaturation of *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substrate. **Figure 5.** Arrhenius-type plot used to estimate the activation energy of thermal denaturation of *Aspergillus tamarii* URM4634 protease produced by sequential fermentation using wheat bran as a substrate.

As for the thermodynamic parameters of protease thermal denaturation, an activation energy (E\*d) of 64.78 kJ mol−1 was estimated from the Arrhenius-type plot illustrated in Figure 5 with satisfactory correlation (R<sup>2</sup> = 0.955). As known, E\*d is related to the activation enthalpy of thermal denaturation (ΔH\*d), i.e., the total amount of thermal energy required to denature the enzyme through disruption of non-covalent bonds, including hy-

of the above E\*<sup>d</sup> and ΔH\*d values with those of protease obtained by SSF either in crude (E\*<sup>d</sup> = 49.7 kJ mol−1; ΔH\*<sup>d</sup> = 46.7–47.0 kJ mol−1) or purified (E\*<sup>d</sup> = 28.8 kJ mol−1; ΔH\*<sup>d</sup> = 25.9– 26.1 kJ mol−1) form [37] confirms the higher thermostability of the enzyme produced by SF. To the best of our knowledge, studies comparing the thermostability of enzymes produced by different fermentation processes are very scarce. Nonetheless, the work of Saqib et al. [38] can be cited, in which the authors reported, on the basis of kinetic and thermodynamic parameters, a greater thermostability of an *A. fumigatus* endoglucanase produced by SSF rather than by SmF. Therefore, further studies are needed to investigate possible structural differences in enzymes produced in different fermentation processes,

. Comparison

Usually, a breakdown of the enzyme structure occurs during thermal denaturation, which is accompanied by an increase in the degree of disorder or randomness, i.e., a positive value of the activation entropy (ΔS\*d) [28]. However, negative values of ΔS\*<sup>d</sup> of thermal denaturation ranging from −133.52 to −136.71 J mol−1 K−1 were estimated for *A. tamarii* URM4634 protease produced by SF, suggesting an even higher order degree of the transition state compared to the separate reactants and then, high thermostability in the investigated range of temperature [39]. When ΔS is negative and ΔH positive, a process is not spontaneous at any temperature, but the reverse process is [29], which means that protease denaturation was reversible and not spontaneous at all temperatures investigated. As for the Gibbs free energy of thermal denaturation (ΔG\*d), small or negative values of this parameter are related to a more spontaneous process; contrariwise, large or positive values like those estimated in this study (105.97 ≤ ΔG\*<sup>d</sup> ≤ 110.12 kJ mol−1) indicate special resistance to denaturation, confirming the thermostability indicated by the other Usually, a breakdown of the enzyme structure occurs during thermal denaturation, which is accompanied by an increase in the degree of disorder or randomness, i.e., a positive value of the activation entropy (∆S\*d) [28]. However, negative values of ∆S\*<sup>d</sup> of thermal denaturation ranging from <sup>−</sup>133.52 to <sup>−</sup>136.71 J mol−<sup>1</sup> <sup>K</sup>−<sup>1</sup> were estimated for *A. tamarii* URM4634 protease produced by SF, suggesting an even higher order degree of the transition state compared to the separate reactants and then, high thermostability in the investigated range of temperature [39]. When ∆S is negative and ∆H positive, a process is not spontaneous at any temperature, but the reverse process is [29], which means that protease denaturation was reversible and not spontaneous at all temperatures investigated. As for the Gibbs free energy of thermal denaturation (∆G\*d), small or negative values of this parameter are related to a more spontaneous process; contrariwise, large or positive values like those estimated in this study (105.97 <sup>≤</sup> <sup>∆</sup>G\*<sup>d</sup> <sup>≤</sup> 110.12 kJ mol−<sup>1</sup> ) indicate special resistance to denaturation, confirming the thermostability indicated by the other kinetic and thermodynamic parameters.

#### kinetic and thermodynamic parameters. **3. Materials and Methods**

#### *3.1. Materials*

Azocasein was purchased from Sigma-Aldrich (St. Louis, MO, USA). Potato dextrose agar (PDA) medium, glucose, yeast extract and all other reagents used were analytical grade. Wheat bran used as a substrate in SmF and SF was acquired in a local market in the city of Garanhuns, PE, Brazil.

#### *3.2. Microorganism and Inoculum Preparation*

The *Aspergillus tamarii* URM4634 strain was provided by the "Micoteca-URM" of Mycology Department, Centre of Biosciences of Federal University of Pernambuco (UFPE), Recife, PE, Brazil, preserved in mineral oil and maintained at 25 ± 1 ◦C in Czapek Dox agar medium. For sporulation, the strain was inoculated in PDA medium for 7 days at 30 ◦C. The spores were suspended in NaCl solution (0.9%) containing Tween 80 (0.01%), and the inoculum was adjusted to a final concentration of 10<sup>5</sup> spores mL−<sup>1</sup> and 10<sup>7</sup> spores g−<sup>1</sup> of substrate for the initial SmF and SF runs, respectively.

#### *3.3. Protease Production by Sequential Fermentation*

The sequential fermentation (SF) was performed at 30 ◦C in two steps. First the solidstate fermentation (SSF) was conducted for 24 h, then the nutrient medium was added, and the submerged fermentation (SmF) was carried out for the same time. The fermentation

was performed in 250 mL-Erlenmeyer flasks containing 5 g of substrate (wheat bran standardized to a particle size between 0.5 and 2.0 mm), nutrition solution (1.0% glucose and 0.5% yeast extract) and a spore suspension (10<sup>7</sup> spores g−<sup>1</sup> ), corresponding to a 60% moisture content. After 24 h of growth, we added, in the proportion of 20 mL g−<sup>1</sup> of substrate, a nutrient medium adapted from Cunha et al. [7] having the following composition: 0.14% (NH4)2SO4, 0.20% KH2PO4, 0.03% CaCl2, 0.02% MgSO4·7H2O, 0.50% peptone, 0.20% yeast extract, 0.03% urea, 0.10% Tween 80, 0.10% mineral solution (5 mg L−<sup>1</sup> FeSO4·7H2O, 1.6 mg L−<sup>1</sup> MnSO4·H2O and 1.4 mg L−<sup>1</sup> ZnSO4·7H2O) and supplemented with glucose (40.0 g L−<sup>1</sup> ). After 24 h of SmF carried out under stirring (150 rpm), the fermented medium was filtered with a vacuum pump, and the crude enzyme extract was stored at −20 ◦C for further analysis.

Protease production by SF was investigated using a 2<sup>3</sup> -full factorial design plus four central points, in which glucose concentration (30, 40 and 50 g L−<sup>1</sup> ), medium volume (15, 20 and 25 mL g−<sup>1</sup> of substrate) and inoculum size (10<sup>6</sup> , 10<sup>7</sup> and 10<sup>8</sup> spores g−<sup>1</sup> of substrate) were selected as the independent variables (Table 1). Based on the results obtained with this statistical design (Table 2), a central composite rotational design (CCRD) was performed to better identify the experimental conditions to maximize protease production, in which the independent variables were glucose concentration and medium volume (Table 3).

The statistical analysis of both experimental designs was performed using the software package Statistica 7.0 (Statsoft Inc., Tulsa, OK, USA). The best fitting quadratic model equation was determined by the Fisher's test. The fitting ability of the model was assessed by the coefficient of determination (R<sup>2</sup> ) and the analysis of variance (ANOVA).

#### *3.4. Protease Production by Submerged Fermentation*

For comparison purposes, submerged fermentations were also performed at 130 rpm and 30 ◦C for 48 h, according to Porto et al. [40] with some modifications, using wheat flour as a substrate in 250 mL-Erlenmeyer flasks containing 50 mL of culture medium with the following composition (%): filtered wheat flour, 4.0; NH4Cl, 0.1; MgSO4·7H2O, 0.06; yeast extract, 0.1; K2HPO4, 0.435; glucose, 0.5; and 1.0% (v v−<sup>1</sup> ) mineral solution (100 mg FeSO4·7H2O; 100 mg MnCl2·4H2O; and 100 mg ZnSO4·H2O in 100 mL of distilled water). The medium pH was adjusted to 7.0 before fermentation.

#### *3.5. Analytical Determinations*

Protease activity was determined as described by Ginther [41] with adaptations using azocasein as a substrate. For this purpose, a 1.0% azocasein solution was prepared in Tris-HCl buffer (0.2 M, pH 7.2) containing 1.0 mM CaCl2. The reaction was performed using 0.15 mL of enzyme extract and 0.25 mL of azocasein solution at room temperature (25 ± 1 ◦C). After 1 h, the reaction was stopped by addition of 1.0 mL of 10% trichloroacetic acid (TCA). Then, the samples were centrifuged at 8000× *g* for 20 min at 4 ◦C, 0.8 mL of supernatant were added to 0.2 mL of 1.8 M NaOH, and the absorbance of samples was measured at 420 nm with a spectrophotometer (Biochrom Libra S6, Cambridge, UK). Protease activity was defined as the enzyme amount responsible for a 0.1 increase in absorbance at 420 nm within 1 h. The total protein concentration was determined by the Bradford method [42] using bovine serum albumin (BSA) as a standard.

#### *3.6. Effect of pH and Temperature on Protease Activity*

The optimum pH of enzyme activity was determined through activity tests carried out at 25 ◦C as described in the previous section using substrate solutions at pH values ranging from 5.0 to 12.0. For this purpose, the following buffers were used at 0.1 M concentration: citrate-phosphate (5.0–7.0), Tris-HCl (7.0–9.0) and glycine-NaOH (9.0–12.0). On the other hand, the influence of temperature on protease activity was investigated through activity tests carried out at different temperatures (30–80 ◦C) and constant pH (7.2). The results were evaluated by the Tukey's test to verify statistically significant differences (*p* < 0.05) among different samples at different values of pH or temperature.

#### *3.7. Kinetic and Thermodynamic Parameters of Azocasein Hydrolysis*

The Michaelis constant (Km) and maximum rate (Vmax) of azocasein hydrolysis were estimated by the Lineweaver-Burk double reciprocal plot. For this purpose, protease activity assays were performed using different substrate concentrations (2.0–100 mg mL−<sup>1</sup> ) at pH 7.2 and 25 ± 1 ◦C. The catalytic constant or turnover number (kcat) was calculated as the ratio of Vmax to the protein concentration. The activation energy (E\* <sup>a</sup>) of substrate hydrolysis was estimated by the Arrhenius-type plot of lnA<sup>0</sup> versus 1/T, being A<sup>0</sup> the starting activity and T the absolute temperature, using the results of enzyme activity tests. For a better understanding of temperature influence on the hydrolysis reaction, the temperature quotient (Q10), which is the factor by which the enzyme activity increases due to a 10 ◦C temperature increase, was calculated by the following equation [43]:

$$\mathbf{Q}\_{10} = \text{anti } \log\left(\frac{\mathbf{E}\_{\text{a}}^{\*} \times 10}{\mathbf{RT}^{2}}\right) \tag{1}$$

where R is the gas constant (8.314 J mol−<sup>1</sup> K −1 ).

The thermodynamic parameters, namely the activation enthalpy (∆H\*), Gibbs free energy (∆G\*) and entropy (∆S\*) of azocasein hydrolysis, were calculated using the following equations:

$$
\Delta \mathbf{H}^\* = \mathbf{E\_a^\*} - \mathbf{R} \mathbf{T} \tag{2}
$$

$$
\Delta \mathbf{G}^\* = -\mathbf{R} \mathbf{T} \ln \left(\frac{\mathbf{k}\_{\text{cat}} \mathbf{h}}{\mathbf{k}\_{\text{b}} \mathbf{T}}\right) \tag{3}
$$

$$
\Delta \mathbf{S}^\* = \frac{\Delta \mathbf{H}^\* - \Delta \mathbf{G}^\*}{\mathbf{T}} \tag{4}
$$

where h (6.626 <sup>×</sup> <sup>10</sup>−<sup>34</sup> J s) is the Planck constant, k<sup>b</sup> (1.381 <sup>×</sup> <sup>10</sup>−<sup>23</sup> J K−<sup>1</sup> ) the Boltzmann constant and *T* is the reference temperature (298.15 K).

#### *3.8. Kinetic and Thermodynamic Parameters of Protease Thermal Denaturation*

As known, the increase in temperature not only accelerates the enzyme reaction but also promotes the reversible enzyme unfolding and the subsequent denaturation of its structure [31]. The enzyme thermal denaturation can then be kinetically described by the equation:

$$\frac{d\mathbf{A}}{dt} = -\mathbf{k}\_{\mathbf{d}} \cdot \mathbf{A} \tag{5}$$

where k<sup>d</sup> is the rate constant of denaturation, A the enzyme activity and t the time.

Accordingly, the k<sup>d</sup> values were estimated through the linearized version of this equation:

$$\ln\left(\frac{\mathbf{A}}{\mathbf{A}\_0}\right) = -\mathbf{k}\_\mathbf{d}\mathbf{t} \tag{6}$$

from the slopes of the straight lines obtained by plotting the experimental data of ln(A/A0) versus time in the temperature range of 50 to 80 ◦C, where A<sup>0</sup> is the starting activity at t = 0.

From the k<sup>d</sup> values, we also calculated the half-life (t1/2), which is defined as the time after which the enzyme activity is reduced to one-half the initial value:

$$\mathbf{t}\_{1/2} = \frac{\ln 2}{\mathbf{k}\_d} \tag{7}$$

and the decimal reduction time (D-value), i.e., the time required for a 10-fold activity reduction:

$$\text{D}-\text{value} = \frac{\ln 10}{\text{k}\_{\text{d}}} \tag{8}$$

The D-value is often accompanied by the thermal resistance constant (Z-value) that corresponds to the temperature increase required to achieve a 10-fold reduction of the D-value. This parameter was estimated from the slope of the thermal death time plot of logD versus T (◦C).

The thermodynamic parameters related to protease thermal denaturation were also estimated from an Arrhenius-type plot. Particularly, the activation energy (E\*d) was estimated from the slope of the straight line of lnk<sup>d</sup> vs. 1/T, while the activation enthalpy (∆H\*d), Gibbs free energy (∆G\*d) and entropy (∆S\*d) were calculated similarly to those of enzyme-catalyzed reaction using k<sup>d</sup> instead of kcat in Equations (2)–(4).

#### **4. Conclusions**

The sequential fermentation (SF) process proved to be an interesting approach for protease production, showing higher outcomes (180.17 U mL−<sup>1</sup> ) compared to the conventional submerged fermentation (20.0 U mL−<sup>1</sup> ). Several differences in the enzyme characteristics, especially an improved thermostability, were observed compared to the protease obtained in a previous study by solid-state fermentation using the same fungal strain and substrate, thereby indicating the need to carry out further studies on the protein structure for a better understanding of SF features. The enzyme produced has been classified as a neutral protease (optimum pH of 7.0), which would offer several advantages in typical food industry applications such as baking and brewing processes. Moreover, the kinetic and thermodynamic study pointed out a satisfactory thermostability of the protease at 50 ◦C, a temperature commonly used in many industrial processes, evidenced especially by a half-life of 231.05 min and a decimal reduction time of 767.53 min. The results obtained in this study can contribute to the development of the SF technique to produce proteases and various other enzymes of industrial interest.

**Author Contributions:** Conceptualization, R.L.d.O., A.C. and T.S.P.; methodology, R.L.d.O., A.C. and T.S.P.; software, R.L.d.O. and T.S.P.; validation, A.C. and T.S.P.; formal analysis, R.L.d.O., E.d.S.C. and T.S.P.; investigation, R.L.d.O. and E.d.S.C.; resources, T.S.P.; data curation, R.L.d.O. and A.C.; writing—original draft preparation, R.L.d.O. and E.d.S.C.; writing—review and editing, R.L.d.O., A.C. and T.S.P.; visualization, A.C. and T.S.P.; supervision, R.L.d.O.; project administration, T.S.P.; funding acquisition, T.S.P. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Coordenacão de Aperfeicoamento de Pessoal de Nível Superior (CAPES, Brazil) grant number 23038.003634/2013–15 and the National Council for Scientific and Technological Development (CNPq, Brazil) grant number 471773/2013–1.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Acknowledgments:** Rodrigo Lira de Oliveira is grateful to FACEPE (Foundation for Science and Technology of the State of Pernambuco, Brazil) for his post-doctoral scholarship (grant BFP-0105- 5.07/20). The authors are grateful to the Coordenacão de Aperfeicoamento de Pessoal de Nível Superior (CAPES, Brazil) and the National Council for Scientific and Technological Development (CNPq, Brazil) for financial funding, and the Federal Rural University of Pernambuco for the laboratory infrastructure.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

