|
|
||||||||
1 Department of Food Bioengineering, Faculty of Food Science and Engineering, University Dunarea de Jos, 111 Domneasca Street, 800201 Galati, Romania
2 Laboratory of Food Technology, Department of Food and Microbial Technology, Katholieke Universiteit Leuven, Kasteelpark Arenberg 22, B-3001 Heverlee, Belgium
Corresponding author: D. Borda; e-mail: daniela.borda{at}ugal.ro.
| ABSTRACT |
|---|
|
|
|---|
Key Words: plasmin high pressure kinetic mathematical model
Abbreviation key: HP = high pressure
| INTRODUCTION |
|---|
|
|
|---|
Generally, HP treatment at room temperature allows inactivation of vegetative bacterial cells in the range 300 to 600 MPa, whereas microbial spores are more pressure-stable, and inactivation requires HP (typically above 700 MPa) in combination with elevated temperature (typically above 70°C) (Hendrickx et al., 1998; Heinz and Knorr, 2001). Hence, considering the specific properties of milk constituents and the problems associated with dairy product inocuity, a combined HP mild-temperature treatment could be a more appropriate approach. Some of these applications include shelf-life extension of milk, the acceleration of cheese ripening, or the design of novel dairy products with new or improved functional properties. In this context, understanding the complex phenomena that take place during combined HP/thermal treatments requires further study. Detailed analysis of pressure-temperature inactivation kinetics for several milk components in model systems can help to apply such treatments to real systems.
Plasmin (EC 3.4.21.7) is a native milk proteinase that can play an important role in coagulation and ripening of many cheese varieties and in influencing the stability of UHT milk. This enzyme is very important in regulating casein hydrolysis in dairy products and readily hydrolyzes ß-casein,
s2-casein and, more slowly,
s1-casein. Plasmin is an alkaline serine proteinase with a pH optimum of 7.5 (Fox and McSweeney, 1996). As part of a complex system comprising plasminogen (representing the zymogen), plasminogen activators, plasminogen activator inhibitors, and plasmin inhibitors, plasmin exhibits interesting behavior during HP thermal inactivation. Plasminogen, the inactive form of plasmin, can be converted into plasmin either in the presence of urokinase or tissue type activators (Nielsen, 2003).
In a previous study (Borda et al., 2004), we investigated combined thermal and HP inactivation of plasmin in a model system containing plasminogen and plasmin, and the kinetics in terms of inactivation rate constants and activation energies was reported.
The aim of the current research was to quantify the combined thermal and HP inactivation kinetics of the plasmin system, after conversion of plasminogen into plasmin with commercially available urokinase, and to compare the results obtained with those for the system containing both plasminogen and plasmin.
At the same time, the second objective was to develop mathematical models for the combined pressure-temperature inactivation of plasmin in the 2 systems, namely a plasmin system with both plasminogen and plasmin present, and a system in which all plasminogen was converted into plasmin with urokinase. Furthermore, we tried to relate differences in the proposed models to possible differences in structural changes that occur during HP thermal inactivation in these 2 systems.
| MATERIALS AND METHODS |
|---|
|
|
|---|
To prepare the system with plasminogen converted into plasmin, 100 µL of plasmin system (containing both plasmin and plasminogen) was mixed with 100 µL of urokinase solution (1000 Ploug U/mL) and 1000 µL of 0.05 M phosphate buffer, pH 6.6, and then incubated for 2 h at 37°C.
Isobaric-isothermal (pressure/temperature) treatments were conducted in multivessel (8 vessels of 8 mL) HP equipment (Resato, Roden, The Netherlands) with a maximum pressure of 1000 MPa, as described previously (Borda et al., 2004).
After the thermal or combined HP/thermal treatment, the remaining plasmin activity was measured by fluorescence spectrophotometry using 1 mM N-succi-nyl-L-alanyl-L-phenylalanyl-L-lysyl-7-amido-4-methyl coumarin as a substrate, based on a modified procedure described in the literature (Richardson, 1983). The procedure detects the fluorescent intensity of the reaction product between plasmin and coumarin peptide, 7-amino-4-methyl coumarin.
As a basis to formulate a plasmin inactivation model for the system containing both plasminogen and plasmin, 58 different pressure/temperature combinations in the 0.1 to 800 MPa pressure range and 30 to 90°C temperature range were considered. Formulation of an appropriate mathematical model for plasmin inactivation, in the system where plasminogen was converted into plasmin with urokinase, was based on 43 pressure-temperature combinations in the same range.
Data Analysis
The thermal HP inactivation of plasmin in both systems could be described by first-order kinetics. In case of constant extrinsic/intrinsic factors the first order model can be integrated to obtain:
![]() | ([1]) |
where A0 is the initial activity (at t = 0) and A the residual activity at time t, after the treatment.
The temperature dependence at constant pressure of the plasmin inactivation rate constant (kobs) could be described by the Arrhenius equation:
![]() | ([2]) |
where k0 is the inactivation rate constant at reference temperature T0, Ea is the activation energy, and R is the universal gas constant (R = 8.314 J/mol K). The activation energy at constant pressure was estimated using a linear regression analysis.
To express the pressure-temperature dependence of the inactivation rate constant of enzymes, a thermodynamically based equation is frequently used (Smeller and Heremans, 1997; Heremans, 2002; Ludikhuyze et al., 2002).
Using the Eyring transition state theory, this thermodynamic model can be converted into a kinetic model to obtain equation 3, which is frequently used to describe the HP thermal denaturation of proteins in general and enzymes in particular (Heremans and Smeller, 1998; Van Loey et al., 2002):
![]() | ([3]) |
where
Cp
is the heat capacity change (T 
S/
T)P, 

is the thermal expansibility factor (T 
V/
T)P = (
S/
P)T, and
k
is the compressibility factor (
V/
P)T.
Recently, Smeller (2002) suggested a second-degree polynomial approximation based on the following substitution:
![]() | ([4]) |
When applying the resulting model it is important to recognize both the power and the limitations of the model (Heremans, 2002). According to Morlid (1981), this equation, in principle, describes the pressure-temperature behavior of most phenomena. One limitation of this equation is that only the difference in thermal expansion, compressibility factor, and heat capacity between denatured and native state are obtained. Also, it should be taken into account that this model (and thus the elliptical shape of the iso-rate P, T contour) results from the fact that the series expansion is truncated after the second-degree terms (Smeller, 2002). Leaving out higher degree terms implies that
Cp
, 

, and 

are independent of temperature and pressure, and if any of these show pressure or temperature dependence, higher degree terms appear not to be negligible and as a consequence the ellipse can be distorted (Smeller and Heremans, 1997).
In a recent study on pectin methylesterase inactivation (Ly-Nguyen et al., 2003), 3 third-degree terms were required in the model to accurately describe its pressure-temperature inactivation kinetics.
Equations 3 and 4 were used as a starting point to obtain mathematical models for plasmin inactivation in the 2 systems under study, using multilinear regression analysis (SAS, 2001).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
For this system, HP treatment at temperatures exceeding 30°C clearly resulted in plasmin inactivation. Linear curves were obtained when the residual plasmin activity was plotted as a function of time on a log linear scale, indicating first-order kinetics (as observed for thermal inactivation). Using a linear regression approach, the corresponding inactivation rate constants values were estimated (Table 1
).
|
|
Comparison Between the Combined Pressure and Temperature Plasmin Inactivation in the System With and Without Urokinase
The inactivation rate constants of the system containing plasminogen and plasmin in the pressure range 300 to 800 MPa and temperatures between 30 and 65°C are summarized in Table 3
(results from Borda et al., 2004).
|
|
The absence of plasminogen in the system converted with urokinase and the fewer disulfide bonds stabilizing the plasmin structure, as compared with the system with plasminogen and plasmin, resulted in less dramatic changes above 600 MPa, so no antagonistic effect was noticed and only the stabilization phenomenon was observed. Only further in-depth structural studies will allow confirmation of this hypothesis.
Mathematical Models for Inactivation-Rate Constant vs. Pressure and Temperature
In a first attempt, the second-degree polynomial (combining equations 3 and 4) was used to describe the pressure-temperature dependence of the plasmin inactivation-rate constants. For both systems, the second-degree polynomial was not suitable to describe the inactivation-rate constant as function of pressure and temperature (results not shown).
To model the data presented in Figure 1A and B
, and to formulate a general inactivation model, the thermodynamic model (equation 3) was expanded including third-degree terms using a Taylor series development for a function of 2 variables (P,T) resulting in equation 5:
![]() | ([5]) |
Equation 5 can be transformed into a linearized equation (equation 6), so that multilinear regression analysis can be applied.
If we carry out the following substitutions:
![]() |
we obtain:
![]() | ([6]) |
By using multiple linear regression analysis we can perform multivariate tests across the multiple dependent variables and select the number of significant variables in the model. The "forward" selection procedure in the REG multiple linear regression analysis (SAS, 2001) allows one to start with no variables in the model, and for each independent variable, it will calculate the F-statistic that evaluates whether the contribution of the variable to the model is significant. Based on this selection procedure, for the plasmin system containing both plasminogen and plasmin, equation 6 was reduced to:
![]() | ([7]) |
Hence, one extra term (P P0)3 was finally introduced in the second-degree polynomial approach. The simulated iso-rate plot has now a similar shape as the experimental iso-rate plot and showed the same bending (Figure 2A
) at pressures above 600 MPa.
|
![]() | ([8]) |
As compared with equation 7 (presenting the plasmin-plasminogen mixture), equation 8 (presenting the behavior of plasmin only) has an extra term in (P P0)2(T T0) and one term less, representing the (P P0)2 dependence of the inactivation rate constant. The predicted iso-rate contour (Figure 2B
) does not have an elliptical shape because a third-degree term is included. The predicted shape is obviously closer to the experimental iso-rate contour and allows taking into account the stabilization phenomenon for pressures above 600 MPa.
In addition to the corrected correlation coefficient for the model, R2, the Mallows Mw statistic was computed for each model selected; Mw is a measure for the total sum of squared errors and is defined as:
![]() | ([9]) |
where s2 is the mean error sum of squares for the full model, SSEp is the error sum of squares for a model with p parameters including the intercept, N is the number of observation, and p is the number of parameters including the intercept.
The best situation occurs when Mw is very close to p. For equations 7 and 8 the Mw value was equal to p.
The plots of the residuals (differences between experimental and predicted inactivation rate constants) as a function of temperature and pressure showed an unbiased scattering for both cases. The random distribution of the residuals indicates good agreement between the proposed model and the experimental data.
Based on equations 7 and 8, parity charts representing the difference between the natural logarithm of predicted inactivation rate constants and experimental values were plotted, providing a good fit of the experimental data compared with the predicted data (Figure 3A, B
). The deviation from the bisector can be considered an indicator for the accuracy of the proposed model. The estimated parameters and associated standard errors, as well as some statistical information on the quality of the fit, are summarized in Table 4
.
|
|
| CONCLUSIONS |
|---|
|
|
|---|
This effect is of interest for further studies as well as for the formulation of models for HP thermal inactivation of plasmin in systems containing ß-casein and ß-lactoglobulin.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication May 9, 2004. Accepted for publication August 19, 2004.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |