Next Article in Journal
Applying Nutrient Profiling Systems to Packaged Foods and Drinks Sold in Jamaica
Next Article in Special Issue
Quantification of Cheese Yield Reduction in Manufacturing Parmigiano Reggiano from Milk with Non-Compliant Somatic Cells Count
Previous Article in Journal
Consumers’ Perceptions and Sensory Properties of Beef Patty Analogues
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whey Protein Concentrate as a Novel Source of Bifunctional Peptides with Angiotensin-I Converting Enzyme Inhibitory and Antioxidant Properties: RSM Study

by
Fatima Abdelhameed Hussein
1,2,
Shyan Yea Chay
1,
Mohammad Zarei
3,
Shehu Muhammad Auwal
4,
Azizah Abdul Hamid
1,
Wan Zunairah Wan Ibadullah
1 and
Nazamid Saari
1,*
1
Department of Food Science, Faculty of Food Science and Technology, University Putra Malaysia, Serdang, Selangor 43400, Malaysia
2
Department of Dairy Production, Faculty of Animal Production, University of Khartoum, PO Box 32, Khartoum North 13314, Sudan
3
Department of Food Science and Technology, Faculty of Applied Sciences, Universiti Teknologi MARA, Shah Alam, Selangor 40450, Malaysia
4
Department of Biochemistry, Faculty of Basic Medical Sciences, Bayero University, Kano 700231, Nigeria
*
Author to whom correspondence should be addressed.
Foods 2020, 9(1), 64; https://doi.org/10.3390/foods9010064
Submission received: 18 September 2019 / Revised: 7 October 2019 / Accepted: 9 October 2019 / Published: 8 January 2020
(This article belongs to the Special Issue Cheese and Whey)

Abstract

:
Whey protein concentrate (WPC) is a unique source of protein with numerous nutritional and functional values due to the high content of branched-chain amino acid. This study was designed to establish the optimum conditions for Alcalase-hydrolysis of WPC to produce protein hydrolysates with dual biofunctionalities of angiotensin-I converting enzyme (ACE) inhibitory and antioxidant activities via response surface methodology (RSM). The results showed that the optimum conditions were achieved at temperature = 58.2 °C, E/S ratio = 2.5%, pH = 7.5 and hydrolysis time = 361.8 min in order to obtain the maximum DH (89.2%), ACE-inhibition (98.4%), DPPH• radical scavenging activity (50.1%) and ferrous ion chelation (73.1%). The well-fitted experimental data to predicted data further validates the regression model adequacy. Current study demonstrates the potential of WPC to generate bifunctional hydrolysates with ACE inhibition and antioxidant activity. This finding fosters the use of WPC hydrolysate as a novel, natural ingredient for the development of functional food products.

1. Introduction

Hypertension is considered a key risk factor in cardiovascular-related diseases including stroke, peripheral and coronary heart diseases. The total number of adults who had elevated blood pressure was 594 million in 1975 and it increased tremendously to 1.13 billion in 2015 [1], attributed to unhealthy diets and poor lifestyle such as alcohol intake and smoking. By 2025, it is projected that hypertension will affect more than 1.5 billion people globally [2]. Angiotensin-I converting enzyme (ACE) has been identified as the key enzyme which increases blood pressure in human body. It cleaves the inactive decapeptide, angiotensin-I, to an octapeptide, angiotensin-II, which acts as a potent vasoconstrictor, as well as inactivates the vasodilating nonapeptide, bradykinin. The dual action of ACE results in a rise in the blood pressure. Inhibition of ACE by reducing angiotensin II generation and increasing bradykinin generation is necessary to control high blood pressure. In recent years, the increasing health awareness drives consumers to look for natural alternatives that are safer with minimal side effects compared to synthetic blood pressure-controlling drugs. Interestingly, positive results have been garnered with the discovery of bioactive peptides derived from various food proteins.
Whey proteins concentrate (WPC), being a natural valuable peptide source, provides considerable nutritional and health benefits for humans due to its high level of branched-chain amino acids, high protein content (30–80%) [3,4,5] and presence of promising functional molecules including immunoglobulin, α-lactalbumin, lactoperoxidase, albumin, lactoferrin and caseinomacropeptide [6,7]. However, WPC is less studied as compared to whey protein (the crude, non-purified form) and whey protein isolate (higher protein content than WPC). Hydrolysis of these whey products would generate smaller peptide fragments with different bioactivities, such as antioxidative [8], antihypertensive [9,10,11], immunomodulatory [12], antithrombotic [13], antimicrobial [14] and opiate properties [15]. Of all these bioactivities, ACE inhibition and antioxidant activity are extensively studied, but mainly on a separate basis. For instance, Guo, Pan and Tanokura [9] and O’keeffe, Conesa and FitzGerald [16] reported ACE inhibition between 40.0–84.4% for different whey products while Ajlouni and Pan [17] as well as Zhang et al. [18] reported antioxidant activities between 0.0–72.1%, respectively.
Hydrolysis parameters rely heavily on the purpose of hydrolysis. Different purposes would yield different sets of optimum conditions and enzyme selection. WPC has been widely used to evaluate the optimum hydrolysis conditions using different enzymes. For instance, Guo, Pan and Tanokura [9] reported the optimum temperature = 38.9 °C, pH = 9.2 and E/S ratio = 0.60 for L. helveticus LB13 crude proteinase hydrolysis while Naik et al. [19] reported the optimum temperature = 35.5 °C, time = 8 h and pH = 7.3 when Flavourzyme was used. Also, WPC has been hydrolysed using Alcalase, the same enzyme in current work, in several other studies. According to Athira et al. [20], the optimum conditions for Alcalase-assisted WPC hydrolysis were 55 °C, 8 h and pH 9 while Fenoglio et al. [21] suggested to use 50 °C, pH 8–9 and 25 min to achieve optimum hydrolysis. While these studies used the same substrate and same enzyme as in current work, they produced different sets of optimum temperatures, pH and hydrolysis time, depending on the different objectives (i.e., the measured responses) targeted by different works. For instance, Athira et al. [20] focused on single bioactivity of WPC, Fenoglio et al. [21] measured the production of bitter peptides for dessert making while current work measured multiple bioactivities from WPC.
Unlike most works that focused on single bioactivity of WPC, current study evaluates the dual functionality of WPC hydrolysate in terms of ACE inhibition and antioxidant activities by employing 4-factors design using response surface methodology (RSM). RSM is a collection of statistical methods for building a model by setting a combination of factor levels to get an optimum response, by employing a polynomial equation to evaluate the relationship between factors as well as their key and combined effect to induce a desired response [22]. This study serves to determine the optimized conditions to generate WPC hydrolysates that possess maximum ACE-inhibition and antioxidant activities by employing response surface methodology according to Box-Behnken design.

2. Materials and Methods

2.1. Materials

Whey protein concentrate (WPC) was purchased from Purefit Company (Selangor, Malaysia). ACE from rabbit lung, N-Hippuryl-His-Leu hydrate powder and O-phthaldialdehyde (OPA) were obtained from Sigma Chemical Co., (St Louis, MO, USA). Alcalase was purchased from Novozymes (Copenhagen, Denmark). All other chemicals employed were of analytical grade and obtained from either Acros Organics (Geel, Belgium), Merck KGaA (Darmstadt, Germany), Fisher Scientific (Waltham, MA, USA) or J.T. Baker (Phillipsburg, NJ, USA).

2.2. Proximate Analysis of WPC

The proximate analysis was determined according to the method of Association of Official Analytical Chemists [23]. Oven method was employed to measure the moisture content using a forced draft oven (model: Quimis Q-314M242, serial 020, São Paulo, Brazil) at 105 °C until reaching a constant weight. Ash was determined via incineration with a muffle furnace at 550 °C. Soxhlet extraction method with ethylic ether was employed to determine the fat content. The micro Kjeldahl method was used to determine the protein content, whereby a conversion of 6.38 was multiplied. All analyses were performed in triplicates.

2.3. Sodium Dodecyl Sulphate Polyacrylamide Gel Electrophoresis (SDS-PAGE)

SDS-PAGE was performed on WPC hydrolysate as described by Laemmli [24], using 40% acrylamide for protein separation. The hydrolysates were adjusted to a 10 mg protein/mL concentration. The sample solutions and the sample buffer (0.5 M Tris-HCl, pH 6.8, containing 4% SDS, 20% glycerol) were mixed at a 3:1 (v/v) ratio and heated for 10 min at 95 °C, prior to loading. The concentrations of stacking and resolving gels were 4% and 12%, respectively. The gel was run at 150 V. Aliquots of 10 μL were loaded into individual wells and a constant current was passed through the gel for 1 h to complete the peptide separation. The gel was stained using coomassie blue (5.0 g coomassie brilliant blue R-250 mixed with 500 mL ethanol, 100 mL acetic acid and 400 mL deionized water), distained with distaining solution containing 200 mL of 40% ethanol and 80 mL of 6% acetic acid with 720 mL deionized water, then preserved in 5% acetic acid. The molecular weights were determined by approximation using low range molecular weight standards (Bio-Rad, Hercules, CA, USA).

2.4. Response Surface Methodology (RSM) Design for WPC Hydrolysis

While optimum hydrolysis conditions are recommended by the enzyme manufacturer, these parameters are given in ranges without pinpointing exact values. Optimum conditions are subjective and multifactor-dependent (purpose of study, sample nature and equipment availability). Thus, optimum conditions should be specifically considered on a case-by-case basis, justifying the need to perform optimization study on WPC hydrolysis. Box-Behnken design (BBD) was employed for optimisation of WPC hydrolysis parameters, whereby each parameter had three levels coded as −1, 0, +1, based on four independent variables: temperature (°C, x1), enzyme to substrate ratio (E/S, x2), pH (x3) and time (min, x4). BBD was chosen over central composite design because it is more economic and suitable for current study which manipulates four independent variables and four dependent variables (responses). A total of twenty-seven runs with three replicates of the central point was performed. For each treatment, the mean pertaining to triplicate measurements was considered as responses for degree of hydrolysis (%, y1), ACE-inhibition (%, y2), DPPH• (%, y3) and ferrous ion chelating (%, y4). Randomisation of experimental runs was performed to mitigate the impacts of unexpected variability to actual responses.
Briefly, Alcalase was mixed with WPC in 100 mL borate buffer at different ranges of temperatures (40–70 °C), enzyme/substrate (E/S) ratios (0.5%–2.5%), pHs (pH 6.5–10) and times (30–480 min), then incubated in a water bath shaker at an agitation rate of 150 rpm. At the end of hydrolysis, the mixture was placed in boiling water (100 °C) for 10 min to inactivate Alcalase. Upon cooling, centrifugation at 10,000 g was performed for 15 min. The supernatant was stored at −20 °C prior to analysis of degree of hydrolysis (DH), biological activities including ACE-inhibition, DPPH•, and ferrous ions chelating activity.
For data analysis, analysis of variance (ANOVA) and regression equation were employed to determine the regression coefficients, statistical significance of the model terms and to fit the predicted mathematical models with that of the experimental data. The least-squares technique was employed to determine the multiple regression coefficients for predicting quadratic and linear polynomial models pertaining to the calculated response variables [25,26]. The general polynomial model for determining the response is given as follows (Equation (1)):
y = β0 + β1 x1 + β2 x2 + β3 x3 + β4 x4 + β11 x12 + β22 x22 + β33 x32 + β44 x4212 x1 x213 x1 x3 + β14 x1 x4 + β23 x2 x3 + β24 x2 x4 + β34 x3 x4
whereby, y represents the response; β0 signifies the offset term for the design; β1, β2, β3 and β4 denote the regression coefficients defining the linear effect terms; β11, β22, β33 and β44 signify the quadratic effects; β12, β13, β14, β23, β24 and β34 represent the interaction effects; and x1, x2, x3, and x4 represent independent factors in this model.
All linear effects, whether significant or not, are included in the original regression equations (non-reduced model). On the other hand, for square and interaction effects, only significant terms were taken into consideration for reduced polynomial regression pertaining to degree of hydrolysis (y1, Equation (2)), ACE-inhibition (y2, Equation (3)), DPPH• radical scavenging (y3, Equation (4)) and ferrous ion chelating activity (y4, Equation (5)).
y1 = −116.761 + 3.421x1 + 55.453x2 − 6.253x3 + 0.219x4 − 0.056x1 x1 + 0.36x1 x3 − 5.603x2 x3
y2 = −57.8367 − 0.943x1 + 33.1444x2 + 34.294x3 + 0.0542x4 − 0.0286x12 − 3.347x32 − 0.0002x42 + 0.4761x1 x3 + 0.0014x1 x4 − 4.0186x2 x3
y3 = −190.337 +6.167x1 − 0.313x2 + 5.799x3 + 0.158x4 − 0.050x1 x1 − 0.018x3 x4
y4 = −31.8992 + 1.9866x1 − 5.3641x2 + 13.9083x3 + 0.0609x4 − 0.0174x1 x1 − 0.9596x3 x3 − 0.0002x4 x4 + 0.0157x2 x4
For each response (y-value), four important elements, i.e., R2, R2-adj, Fisher values (f-values) and p-values for regression were taken into account to determine the adequacy of the models in predicting experimental values. Based on reduced models, all significant (p < 0.05) interaction effects from two independent factors were visualized as three dimensional (3D) surface plots, with the responses (DH, ACE-inhibitory activity, DPPH• radical scavenging and ferrous ion chelating) plotted on y-axis. Drawing of the 3D plots was obtained by maintaining two constant variables at the centre point while modifying the other two variables within the experimental range. The Minitab statistical package (Minitab Inc., State College, PA, USA) was employed for data analysis, optimisation procedure and experimental design conditions. Specifically, for mathematical optimisation, response optimiser was applied to calculate specific optimum levels of concurrent and individual multiple response optimisations, resulting in appropriate response objectives. The experimental data was compared with the predicted values obtained from the equations to test the adequacy of regression equations. Significant differences were identified at p < 0.05.

2.5. Degree of Hydrolysis (DH)

OPA method was employed following the procedure of Church et al. [27] and Salami et al. [28] with minor modifications to measure DH at various hydrolysis time. To prepare OPA reagent, three solutions, A, B and C, were prepared separately. Solution A comprised of 7.62 g of sodium tetrahydroborate and 200 mg of sodium dodecyl sulphate (SDS) dissolved in 150 mL deionised water; Solution B consisted of 160 mg of OPA dissolved in 4 mL of ethanol 96%; and Solution C consisted of 400 μL of β-mercaptoethanol, in 50 mL deionised water. Mixing the 3 solutions would yield fresh OPA reagent. Then, 36 μL of sample (protein concentration = 20 mg/mL) was mixed with 270 μL of OPA reagent in a 96-well plate, incubated for 2 min at room temperature, and the absorbance was measured at 340 nm by employing a microplate reader system (model: PowerWave X340, Biotek Instruments Inc., Winooski, VT, USA). The following equation was employed to calculate the degree of hydrolysis according to Adler-Nissen [29]:
DH (%) = {(Asample − Aprotein)/Atotal} × 100
in which, Asample represents the sample absorbance after hydrolysis, Aprotein signifies the sample absorbance before hydrolysis (negative control) and Atotal denotes the absorbance of total protein after 24 h complete hydrolysis using 6M HCl.

2.6. Angiotensin Converting Enzyme (Ace) Inhibitory Activity Assay

ACE–inhibitory was determined according to Cushman and Cheung [30] and Ferreira et al. [31] with certain modifications. The assay was performed by first mixing 10 μL of sample (protein concentration = 20 mg/mL) and 10 μL of ACE (100 mU/mL) then incubated for 10 min, followed by addition of 50 μL 0.1 M sodium borate buffer (pH 8.3) comprising of 5 mM HHL and 0.3 M NaCl into the incubated sample. For control and blank, 10 μL and 20 μL of distilled water were respectively added. After incubation for 60 min at 37 °C, the reaction was terminated by introducing 75 μL of 1 M HCl. After the reaction was stopped, addition of 150 μL of pyridine followed by 75 uL of benzene sulphonyl chloride (BSC) was done. Then, the solution was vortexed for 1 min and cooled in ice bath. After cooling, 200 μL of the solution was transferred to a 96-well plate, and the absorbance was measured at 410 nm. ACE inhibition was calculated as follow:
ACE inhibition (%) = {(CS)/(CB)} × 100
where, C is the absorbance of control (ACE + substrate), S is the absorbance of sample (peptide + ACE + substrate) and B is the absorbance of blank (only substrate).

2.7. Determination of Antioxidant Activity

2.7.1. 1, 1-diphenyl-2-picrylhydrazyl Free Radical Scavenging Assay

The percentage for 1, 1-diphenyl-2-picrylhydrazyl (DPPH•) free radical scavenging activity was obtained according to the method of Hwang et al. [32] with some modification. Briefly, 50 μL of WPC hydrolysate (WPCH, protein concentration = 20 mg/mL) was added to100 μL of DPPH• (0.15 mM, 100% methanol) and 50 μL deionised water. At room temperature, the mixture was placed in a dark room for 45 min and then the absorbance was measured at 517 nm. The scavenging ability of WPCH was calculated as follow:
DPPH• radical scavenging activity (%) = {(Acontrol − Asample)/Acontrol} × 100
whereby, Acontrol represents the absorbance of control and Asample denotes the absorbance of samples at 517 nm.

2.7.2. Ferrous Ion Chelating Activity Assay

The ferrous ion chelating activity of WPCH was measured based on a previously defined method put forward by Dinis, Madeira and Almeida [33] with some changes. WPCH (100 μL) of concentrations between 0.005 mg/mL to 10 mg/mL, was mixed with 185 μL of double distilled water and 5 μL of 2 mmol/L iron dichloride solution, followed by addition of 10 μL of 5 mmol/L ferrozine solution and mixed in a vigorous manner. At room temperature, the mixture was incubated for 10 min, prior to absorbance measurement at 562 nm.
The following equation was employed to determine the chelating effect:
Metal chelating activity (%) = {(Acontrol − Asample)/Acontrol} × 100
where, Acontrol represents the absorbance of control and Asample denotes the absorbance of samples at 562 nm.

3. Results and Discussion

3.1. Proximate Analysis

Table 1 shows the proximate composition of WPC with the values of 5.64%, 76.13%, 3.91% and 2.96% for moisture, protein, fat and ash, respectively. Except for moisture and carbohydrate, all values remained comparable to that reported by Adamson [3], Suthar, Jana and Balakrishnan [4] and Whetstine, Croissant and Drake [5]. In particular, the protein content of WPC is high. Therefore, WPC is considered a valuable low-cost protein source to generate hydrolysate (a mixture of bioactive peptides) with dual functionalities upon optimized hydrolysis.

3.2. Molecular Weight of WPC Hydrolysates

The size of the generated peptides is essential for ACE inhibitory, antioxidative and other biological activities. SDS-PAGE was performed to investigate the molecular weight distribution of peptides from WPC hydrolysate at different hydrolysis times. Figure 1 shows that Alcalase efficiently degrades major WPC proteins, i.e., β-lactoglobulin and α-lactalbumin (Mw: 15–20 kDa), into peptides of smaller fragments (below 10 kDa) after 30 min of hydrolysis. Also, serum albumin with heavier molecular weight (around 75 kDa), produces no visible band on the gel after 30 min. These observations are in line with the work performed by Pena-Ramos and Xiong [34], who reported a rapid degradation of WPC proteins into smaller peptides by protease A (an endopeptidase obtained from Bacillus licheniformis, now widely known as Alcalase).

3.3. Model Performance Appraisal on the Regression Elements

In terms of statistics, five parameters are crucial, which include estimated regression coefficient (R2), adjusted coefficient of determination (R2-adj), lack of fit values, Fisher test value (f-value) and regression p-values. The combination of high R2, a non-significant lack of fit, high f-value and low p-values signify a good fitness of experimental values to predicted values obtained from regression equations, and that the regression model is statistically sufficient to explain the experimental data. The R2-adj further improves the reliability of regression model by taking into consideration only the factors (x-values) that significantly affect the response (y-values).
Table 2 presents the predicted and experimental responses of DH (y1), ACE-inhibitory activity (y2), DPPH• (y3) and ferrous ion chelating (y4) under various hydrolysis conditions, at different combinations of temperature (x1), E/S ratio (x2), pH (x3) and time (x4). It can be seen that most experimental values were similar to predicted values. From Table 3, the values for R2 were 91.3% (y1), 90.0% (y2), 86.8% (y3) and 90.8% (y4) while the result for R2-adj were 87.4%, 83.7%, 71.4% and 80.1% for y1, y2, y3 and y4. The lack of fit values for y1, y2, y3 and y4 were 0.665, 0.060, 0.013 and 0.027 while the f-test values of y1, y2, y3 and y4 were 23.54, 14.33, 9.80 and 12.33, respectively. These results confirmed the validity of selected regression model to predict the responses (y-values).

3.4. Effects of Temperature, E/S Ratio, pH and Time on Different Responses

3.4.1. Degree of Hydrolysis (DH)

From Table 2, DH between 17.9–93.3% was recorded for hydrolysis under different combinations of temperature, E/S ratio, pH and time. High DH (above 80%) was observed for 7 runs, out of a total of 27 runs, and this is supported by the work of See, Hoo and Babji [35] who reported a DH of 77.03% for fish skin hydrolysis that used a single enzyme of Alcalase. High DH indicates that the selected combinations of temperature, E/S ratio, pH and time were optimum to perform hydrolysis on large, bulky whey protein molecules to effectively cut them into smaller peptides.
The analysis of variance (p < 0.05) and regression coefficients for DH are depicted in Table 3. As shown, the regression model for DH, when expressed in full quadratic equation, is significant at p = 0.000. This confirms the adequacy of model to fit the DH experimental data. The squares effects were identified as temp*temp and time*time while the interaction effects were found to be temp*pH and E/S*pH. The high R2, non-significant lack-of-fit and large f-value implies a good fit of experimental values to predicted model. Figure 2 depicts the interaction effects from independent factors, comprised of temperature, E/S ratio and pH, which contributed to the changes in DH. From Figure 2c, DH increased with hydrolysis temperature, reached a maximum point then decreased following higher temperature, forming a quadratic surface plot. This is due to enzyme denaturation of Alcalase at high temperatures which reduces its efficiency to hydrolyse WPC molecules into smaller peptide fragments, and thus results in lower DH values.

3.4.2. ACE-inhibitory Activity

The analysis of variance (p < 0.05) and regression coefficients for ACE-inhibitory activity were generated from reduced full quadratic model and were compiled in Table 3. Six terms were found to be significant to explain the changes in ACE-inhibitory activity, namely, temp*temp, pH*pH, time*time, temp*pH, temp*time and E/S*pH. The regression model for ACE-inhibitory activity with high R2 of 89.96%, nonsignificant lack-of-fit of 0.06 and large f-value of 14.33 denotes the goodness of fit of experimental results to the predicted ones. The surface plots for ACE-inhibitory activity were constructed based on the interaction effects of temperature, E/S ratio, pH and time (Figure 3).
ACE-inhibitory activity of WPCHs was highly affected by hydrolysis parameters as well as the resulting DH. At higher DH, a higher ACE inhibitory bioactivity was achieved. This result is in line with the work reported by Silvestre et al. [36] and Ghanbari et al. [37], in which the ACE-inhibitory activity of protein hydrolysates was affected by the type of enzymes as well as the degree of hydrolysis. Also, the increased level of E/S ratio and alkali pH resulted in high values of DH to indicate the release of potent peptides which further contributed towards the strong ACE-inhibitory activity. These findings are comparable to another work reported by van der Ven et al. [38] which used the same sample of whey protein hydrolysate.

3.4.3. Antioxidant Activities

Unlike typical works which perform only single antioxidant assay in RSM study, current work reported two antioxidant activities, i.e., DPPH• radical scavenging and ferrous chelating activities. The reduced quadratic model of DPPH• radical scavenging is shown in Table 3. Even though this model contains only two significant terms, contributed by temp*temp and pH*time, the overall regression equation is significant (p = 0.000) due to the strong effect from linear, squares and interaction. Figure 4 presents the 3D plots regarding DPPH• radical scavenging activity. Similar to the effect of temperature on DH, as previously discussed in Section 3.4, the effect of temperature on DPPH• radical scavenging activity is also prominent. As the hydrolysis temperature increased, the scavenging activity increased until a maximum point, then decreased with higher temperature. This is, again, due to the application of temperature which allow/deter optimum Alcalase activity to generate peptides with DPPH• radical scavenging activity. The ability for scavenging free radicals could be attributed to the peptides’ ability to donate hydrogen, neutralise or stabilise free radicals and then cease their propagation. This potency is strongly affected by the amino acid composition in samples, especially aromatic amino acids like tryptophan, phenylalanine, histidine and tyrosine as well as hydrophobic amino acids like alanine, leucine and valine, with methionine playing the most critical role on affecting the scavenging activity [39,40,41].
Ferrous chelating activity, which helps to minimise lipid oxidation, was exhibited by WPC hydrolysate. The quadratic model for this bioactivity is depicted in Table 3. The model is significant (p = 0.000) and encompasses all linear, squares, interaction effects with p values of 0.000, 0.000 and 0.024, respectively. Four significant terms were reported for ferrous chelating activity, which consisted of mainly squares effect, i.e., temp*temp, pH*pH, time*time, along with an interaction effect of E/S*time. The 3D plots for ferrous chelating activity are shown in Figure 5. The potency of peptides are thought to be affected by sources of protein, type of proteases, amino acid composition and time of hydrolysis [42,43]. The higher ferrous ion chelating activity demonstrated by WPCH proves the suitability of Alcalase as the choice of enzyme to produce biopeptide molecules.

3.5. Optimized Hydrolysis Condition For WPC and Model Validation

Optimization procedure was performed using Box-Behnken design. Figure 6 presents the optimum predicted hydrolysis conditions for WPC that could theoretically produces the maximum DH, ACE-inhibitory, DPPH• radical scavenging and ferrous ion chelating activities. The predicted optimal conditions were as follow: temperature = 58.2 °C, E/S ratio = 2.5%, pH = 7.5, and hydrolysis time = 361.8 min. The predicted parameters were employed in actual experiments for model validation purpose. The model was confirmed to be validate via two observations: (i) The overall desirability (D-value) of 0.88 indicated a high level of confidence for the model to produce the responses as predicted; and (ii) The experimental values for DH (89.2% vs. theoretical 89.6%), ACE-inhibitory activity (98.4% vs. theoretical 98.8%), DPPH• radical scavenging (50.1% vs. theoretical 50.6%) and ferrous ion chelating (73.1% vs. theoretical 74.0%) were identical to their respective predicted values. On the other hand, the goodness-of-fit between experimental and predicted data is depicted in Figure 7a–d. In all figures, the data points fall closely to the regression line. Also, high R2 of 91.3%, 90.0%, 74.6% and 84.6% for DH, ACE inhibitory, DPPH• scavenging and ferrous ion chelating, respectively, indicated that most of the variabilities of these responses were explained by their linear models.

4. Conclusions

This study represents a pioneer work on the optimization of whey protein concentrate hydrolysis using Alcalase to generate biopeptides with dual functionalities (ACE inhibition and antioxidant activities). The optimized hydrolysis parameters were determined as follow: Temperature = 58.2 °C, E/S ratio = 2.5%, pH = 7.5 and hydrolysis time = 361.8 min. Under such conditions, the experimental values did not show any statistical differences from the predicted values but fitted closely to them. Results showed that the selected hydrolysis parameters were able to produce whey protein hydrolysate with dual bioactivities of ACE inhibition and antioxidant at desirable levels. This study proved the reliability of the selected regression models to sufficiently explain the factor-response relationship during WPC hydrolysis and that the predicted optimum hydrolysis conditions are valid to generate, from WPC, bioactive peptides with ACE inhibitory and antioxidative activities. Further work on the in vivo efficacy and subsequent clinical trials of WPC hydrolysate is deemed appropriate.

Author Contributions

Conceptualization: A.A.H., N.S.; data curation: S.Y.C., M.Z.; investigation: F.A.H.; methodology: S.M.A., W.Z.W.I.; supervision and resources: N.S., writing—original draft: F.A.H., S.Y.C., M.Z., S.M.A. writing—review and editing: A.A.H., W.Z.W.I., N.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by Universiti Putra Malaysia under project GP-IPS/2018/9643400 (vot no: 9643400) and the article processing charges is funded by Universiti Putra Malaysia.

Conflicts of Interest

All authors declare no conflict of interest.

References

  1. Zhou, B.; Bentham, J.; Di Cesare, M.; Bixby, H.; Danaei, G.; Cowan, M.J.; Paciorek, C.J.; Singh, G.; Hajifathalian, K.; Bennett, J.E. Worldwide trends in blood pressure from 1975 to 2015: A pooled analysis of 1479 population-based measurement studies with 19 1 million participants. Lancet 2017, 389, 37–55. [Google Scholar] [CrossRef] [Green Version]
  2. Kearney, P.M.; Whelton, M.; Reynolds, K.; Muntner, P.; Whelton, P.K.; He, J. Global burden of hypertension: Analysis of worldwide data. Lancet 2005, 365, 217–223. [Google Scholar] [CrossRef]
  3. Adamson, N. Whey processing. In Dairy Processing Handbook; Bylund, G., Ed.; Tetra Pak Processing Systems: Lund, Sweden, 2015; pp. 331–352. [Google Scholar]
  4. Suthar, J.; Jana, A.; Balakrishnan, S. High protein milk ingredients-A tool for value-addition to dairy and food products. J. Dairy Vet. Anim. Res. 2017, 6, 00171. [Google Scholar]
  5. Whetstine, M.C.; Croissant, A.; Drake, M. Characterization of dried whey protein concentrate and isolate flavor. J. Dairy Sci. 2005, 88, 3826–3839. [Google Scholar] [CrossRef]
  6. del Mar Contreras, M.; Hernández-Ledesma, B.; Amigo, L.; Martín-Álvarez, P.J.; Recio, I. Production of antioxidant hydrolyzates from a whey protein concentrate with thermolysin: Optimization by response surface methodology. LWT Food Sci. Technol. 2011, 44, 9–15. [Google Scholar] [CrossRef]
  7. Muro Urista, C.; Álvarez Fernández, R.; Riera Rodriguez, F.; Arana Cuenca, A.; Tellez Jurado, A. Production and functionality of active peptides from milk. Food Sci. Technol. Int. 2011, 17, 293–317. [Google Scholar] [CrossRef] [PubMed]
  8. Lin, S.; Tian, W.; Li, H.; Cao, J.; Jiang, W. Improving antioxidant activities of whey protein hydrolysates obtained by thermal preheat treatment of pepsin, trypsin, alcalase and flavourzyme. Int. J. Food Sci. Technol. 2012, 47, 2045–2051. [Google Scholar] [CrossRef]
  9. Guo, Y.; Pan, D.; Tanokura, M. Optimisation of hydrolysis conditions for the production of the angiotensin-I converting enzyme (ACE) inhibitory peptides from whey protein using response surface methodology. Food Chem. 2009, 114, 328–333. [Google Scholar] [CrossRef]
  10. Morais, H.A.; Silvestre, M.P.C.; Silva, M.R.; Silva, V.D.M.; Batista, M.A.; e Silva, A.C.S.; Silveira, J.N. Enzymatic hydrolysis of whey protein concentrate: Effect of enzyme type and enzyme: Substrate ratio on peptide profile. J. Food Sci. Technol. 2015, 52, 201–210. [Google Scholar] [CrossRef]
  11. Tavares, T.; Malcata, F. Whey proteins as source of bioactive peptides against hypertension. In Bioactive Food Peptides in Health and Disease; Hernández-Ledesma, B., Ed.; InTech Publisher: London, UK, 2013; pp. 75–114. [Google Scholar]
  12. Wu, Q.; Zhang, X.; Jia, J.; Kuang, C.; Yang, H. Effect of ultrasonic pretreatment on whey protein hydrolysis by alcalase: Thermodynamic parameters, physicochemical properties and bioactivities. Process Biochem. 2018, 67, 46–54. [Google Scholar] [CrossRef]
  13. Silveira, S.T.; Martínez Maqueda, D.; Recio, I.; Hernández Ledesma, B. Dipeptidyl peptidase-IV inhibitory peptides generated by tryptic hydrolysis of a whey protein concentrate rich in β-lactoglobulin. Food Chem. 2013, 141, 1072–1077. [Google Scholar] [CrossRef] [PubMed]
  14. Boyacı, D.; Korel, F.; Yemenicioğlu, A. Development of activate-at-home-type edible antimicrobial films: An example pH-triggering mechanism formed for smoked salmon slices using lysozyme in whey protein films. Food Hydrocoll. 2016, 60, 170–178. [Google Scholar] [CrossRef] [Green Version]
  15. Madureira, A.; Tavares, T.; Gomes, A.M.P.; Pintado, M.; Malcata, F.X. Invited review: Physiological properties of bioactive peptides obtained from whey proteins. J. Dairy Sci. 2010, 93, 437–455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. O’keeffe, M.B.; Conesa, C.; FitzGerald, R.J. Identification of angiotensin converting enzyme inhibitory and antioxidant peptides in a whey protein concentrate hydrolysate produced at semi-pilot scale. Int. J. Food Sci. Technol. 2017, 52, 1751–1759. [Google Scholar] [CrossRef]
  17. Ajlouni, S.; Pan, Y. Effect of pH and whey protein isolate to glucose ratios on the formation of Maillard reaction products as antioxidants. Kasetsart J. (Nat. Sci.) 2014, 48, 759–768. [Google Scholar]
  18. Zhang, Q.X.; Wu, H.; Ling, Y.F.; Lu, R.R. Isolation and identification of antioxidant peptides derived from whey protein enzymatic hydrolysate by consecutive chromatography and Q-TOF MS. J. Dairy Res. 2013, 80, 367–373. [Google Scholar] [CrossRef]
  19. Naik, L.; Mann, B.; Bajaj, R.; Sangwan, R.; Sharma, R. Process optimization for the production of bio-functional whey protein hydrolysates: Adopting response surface methodology. Int. J. Pept. Res. 2013, 19, 231–237. [Google Scholar] [CrossRef]
  20. Athira, S.; Mann, B.; Saini, P.; Sharma, R.; Kumar, R.; Singh, A.K. Production and characterisation of whey protein hydrolysate having antioxidant activity from cheese whey. J. Sci. Food Agric. 2015, 95, 2908–2915. [Google Scholar] [CrossRef]
  21. Fenoglio, C.; Vierling, N.; Manzo, R.; Ceruti, R.; Sihufe, G.; Mammarella, E. Whey protein hydrolysis with free and immobilized alcalase®: Effects of operating parameters on the modulation of peptide profiles obtained. Am. J. Food Technol. 2016, 11, 152–158. [Google Scholar]
  22. Bezerra, M.; Santelli, R.; Oliveira, E.; Villar, L.; Escaleira, L. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta 2008, 76, 965–977. [Google Scholar] [CrossRef]
  23. AOAC International. Official Methods of Analysis of AOAC International, 18th ed.; Association of Official Analytical Chemists (AOAC) International: Rockville, MD, USA, 2005. [Google Scholar]
  24. Laemmli, U.K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature 1970, 227, 680. [Google Scholar] [CrossRef] [PubMed]
  25. Alexopoulos, E. Introduction to multivariate regression analysis. Hippokratia 2010, 14, 23–28. [Google Scholar] [PubMed]
  26. Mayer, R.; Montgomery, D. Response Surface Methodology: Process and Product Optimization Using Designed Experiments; John Wiley and Sons, Inc.: New York, NY, USA, 1995. [Google Scholar]
  27. Church, F.; Swaisgood, H.; Porter, D.; Catignani, G. Spectrophotometric assay using o-phthaldialdehyde for determination of proteolysis in milk and isolated milk proteins. J. Dairy Sci. 1983, 66, 1219–1227. [Google Scholar] [CrossRef]
  28. Salami, M.; Yousefi, R.; Ehsani, M.R.; Dalgalarrondo, M.; Chobert, J.M.; Haertlé, T.; Razavi, S.H.; Saboury, A.A.; Niasari-Naslaji, A.; Moosavi-Movahedi, A.A. Kinetic characterization of hydrolysis of camel and bovine milk proteins by pancreatic enzymes. Int. Dairy J. 2008, 18, 1097–1102. [Google Scholar] [CrossRef]
  29. Adler-Nissen, J. Enzymic Hydrolysis of Food Proteins; Elsevier Applied Science Publishers: New York, NY, USA, 1986. [Google Scholar]
  30. Cushman, D.; Cheung, H. Spectrophotometric assay and properties of the angiotensin-converting enzyme of rabbit lung. Biochem. Pharm. 1971, 20, 1637–1648. [Google Scholar] [CrossRef]
  31. Ferreira, I.; Pinho, O.; Mota, M.; Tavares, P.; Pereira, A.; Goncalves, M.; Torres, D.; Rocha, C.; Teixeira, J. Preparation of ingredients containing an ACE-inhibitory peptide by tryptic hydrolysis of whey protein concentrates. Int. Dairy J. 2007, 17, 481–487. [Google Scholar] [CrossRef] [Green Version]
  32. Hwang, J.Y.; Shyu, Y.S.; Wang, Y.T.; Hsu, C.K. Antioxidative properties of protein hydrolysate from defatted peanut kernels treated with esperase. LWT Food Sci. Technol. 2010, 43, 285–290. [Google Scholar] [CrossRef]
  33. Dinis, T.C.; Madeira, V.M.; Almeida, L.M. Action of phenolic derivatives (acetaminophen, salicylate, and 5-aminosalicylate) as inhibitors of membrane lipid peroxidation and as peroxyl radical scavengers. Arch. Biochem. Biophys. 1994, 315, 161–169. [Google Scholar] [CrossRef]
  34. Pena-Ramos, E.; Xiong, Y. Antioxidative activity of whey protein hydrolysates in a liposomal system. J. Dairy Sci. 2001, 84, 2577–2583. [Google Scholar] [CrossRef]
  35. See, S.F.; Hoo, L.L.; Babji, A.S. Optimization of enzymatic hydrolysis of Salmon (Salmo salar) skin by Alcalase. Int. Food Res. J. 2011, 18, 1359–1365. [Google Scholar]
  36. Silvestre, M.P.C.; Silva, M.R.; Silva, V.D.M.; Souza, M.W.S.; Junior, L.; Oliveira, C.; Afonso, W.O. Analysis of whey protein hydrolysates: Peptide profile and ACE inhibitory activity. Braz. J. Pharm. Sci. 2012, 48, 747–757. [Google Scholar] [CrossRef] [Green Version]
  37. Ghanbari, R.; Zarei, M.; Ebrahimpour, A.; Abdul-Hamid, A.; Ismail, A.; Saari, N. Angiotensin-I converting enzyme (ACE) inhibitory and anti-oxidant activities of sea cucumber (Actinopyga lecanora) hydrolysates. Int. J. Mol. Sci. 2015, 16, 28870–28885. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. van der Ven, C.; Gruppen, H.; de Bont, D.; Voragen, A. Optimisation of the angiotensin converting enzyme inhibition by whey protein hydrolysates using response surface methodology. Int. Dairy J. 2002, 12, 813–820. [Google Scholar] [CrossRef]
  39. Chi, C.F.; Hu, F.Y.; Wang, B.; Li, Z.R.; Luo, H.Y. Influence of amino acid compositions and peptide profiles on antioxidant capacities of two protein hydrolysates from skipjack tuna (Katsuwonus pelamis) dark muscle. Mar. Drugs 2015, 13, 2580–2601. [Google Scholar] [CrossRef]
  40. Luo, S.; Levine, R.L. Methionine in proteins defends against oxidative stress. FASEB J. 2009, 23, 464–472. [Google Scholar] [CrossRef] [Green Version]
  41. Zou, T.B.; He, T.P.; Li, H.B.; Tang, H.W.; Xia, E.Q. The structure-activity relationship of the antioxidant peptides from natural proteins. Molecules 2016, 21, 72. [Google Scholar] [CrossRef]
  42. Karamać, M.; Kosińska-Cagnazzo, A.; Kulczyk, A. Use of different proteases to obtain flaxseed protein hydrolysates with antioxidant activity. Int. J. Mol. Sci. 2016, 17, 1027. [Google Scholar] [CrossRef] [Green Version]
  43. O’Loughlin, I.B.; Kelly, P.M.; Murray, B.A.; FitzGerald, R.J.; Brodkorb, A. Molecular characterization of whey protein hydrolysate fractions with ferrous chelating and enhanced iron solubility capabilities. J. Agric. Food Chem. 2015, 63, 2708–2714. [Google Scholar] [CrossRef]
Figure 1. SDS-PAGE of native whey protein concentrate (WPC) and Alcalase-treated WPC at different hydrolysis hours.
Figure 1. SDS-PAGE of native whey protein concentrate (WPC) and Alcalase-treated WPC at different hydrolysis hours.
Foods 09 00064 g001
Figure 2. 3D surface plots for degree of hydrolysis (%) as affected by (a) E/S ratio and pH; (b) Time and temperature and (c) Temperature and E/S ratio.
Figure 2. 3D surface plots for degree of hydrolysis (%) as affected by (a) E/S ratio and pH; (b) Time and temperature and (c) Temperature and E/S ratio.
Foods 09 00064 g002
Figure 3. 3D surface plots for ACE-inhibitory activity (%) as affected by (a) Temperature and pH; (b) Time and pH and (c) Temperature and E/S ratio.
Figure 3. 3D surface plots for ACE-inhibitory activity (%) as affected by (a) Temperature and pH; (b) Time and pH and (c) Temperature and E/S ratio.
Foods 09 00064 g003
Figure 4. 3D surface plots for DPPH• radical scavenging (%) as affected by (a) Time and pH; (b) Time and temperature and (c) Temperature and E/S ratio.
Figure 4. 3D surface plots for DPPH• radical scavenging (%) as affected by (a) Time and pH; (b) Time and temperature and (c) Temperature and E/S ratio.
Foods 09 00064 g004
Figure 5. 3D surface plots for ferrous chelating activity as affected by (a) Time and pH; (b) pH and temperature and (c) Temperature and E/S ratio.
Figure 5. 3D surface plots for ferrous chelating activity as affected by (a) Time and pH; (b) pH and temperature and (c) Temperature and E/S ratio.
Foods 09 00064 g005
Figure 6. The predicted maximum responses (DH, ACE-inhibitory, DPPH• radical scavenging, ferrous ion chelating activity) based on predicted optimum hydrolysis parameters (temperature, E/S ratio, pH and hydrolysis time). Composite desirability = overall desirability encompassing all responses and d = individual desirability. Optimum hydrolysis parameters are denoted as red lines while the maximum responses are denoted as blue dotted lines.
Figure 6. The predicted maximum responses (DH, ACE-inhibitory, DPPH• radical scavenging, ferrous ion chelating activity) based on predicted optimum hydrolysis parameters (temperature, E/S ratio, pH and hydrolysis time). Composite desirability = overall desirability encompassing all responses and d = individual desirability. Optimum hydrolysis parameters are denoted as red lines while the maximum responses are denoted as blue dotted lines.
Foods 09 00064 g006
Figure 7. Fitted line plots for predicted and experimental values for (a) degree of hydrolysis, (b) ACE-inhibition, (c) DPPH• radical scavenging activity and (d) ferrous ion chelating activity. All values are reported in %.
Figure 7. Fitted line plots for predicted and experimental values for (a) degree of hydrolysis, (b) ACE-inhibition, (c) DPPH• radical scavenging activity and (d) ferrous ion chelating activity. All values are reported in %.
Foods 09 00064 g007aFoods 09 00064 g007b
Table 1. Comparison of proximate composition of whey protein concentrate with literature. All values are reported in dry basis (%).
Table 1. Comparison of proximate composition of whey protein concentrate with literature. All values are reported in dry basis (%).
ParametersPresent Study (Mean ± sd)Literature Data *
Moisture (%)5.64 ± 0.033.8–4.5
Protein (%)76.13 ± 0.2374.8–80.0
Crude fat (%)3.91 ± 0.751.8–7.2
Ash (%)2.96 ± 0.012.7–7.5
Carbohydrates (%)11.93 ± 0.003.5–8.9
* Reported data were compiled from three literatures: Adamson [3] Suthar et al. [4] and Whetstine et al. [5]. sd = standard deviation from 3 readings.
Table 2. Behnken design, predicted and experimental values of four response variables (y1, y2, y3, y4) as affected by temperature (°C), E/S ratio (%), pH and hydrolysis time (min).
Table 2. Behnken design, predicted and experimental values of four response variables (y1, y2, y3, y4) as affected by temperature (°C), E/S ratio (%), pH and hydrolysis time (min).
Run Orderx1x2x3x4% DH (y1)%ACE Inhibition (y2)% DPPH• (y3)% Ferrous Ion Chelating (y4)
ExperimentalPredictedExperimentalPredictedExperimentalPredictedExperimentalPredicted
1400.58.2525530.7145.2227.5326.8527.5326.8573.7872.09
2700.58.2525544.6352.3651.1046.2851.1046.2876.3074.33
3402.58.2525567.4263.6720.9726.2320.9726.2365.0469.37
4702.58.2525562.0970.8147.1345.6547.1345.6574.7071.61
5551.56.53022.6725.2231.8735.5931.8735.5973.4868.97
6551.5103042.6443.2951.4754.0351.4754.0358.0062.23
7551.56.548087.3577.5049.1055.0349.1055.0365.6364.82
8551.51048091.0795.5740.7745.5340.7745.5358.3758.08
9401.58.253017.9018.0934.2023.8134.2023.8163.4463.51
10701.58.253034.4625.2341.3343.2341.3343.2364.0465.75
11401.58.2548072.3770.3721.6729.2721.6729.2760.9359.36
12701.58.2548066.5277.5251.7748.751.7748.759.0061.60
13550.56.525540.5842.5442.8045.6242.845.6277.0077.55
14552.56.525579.6580.6057.7345.0057.7345.0074.9674.83
15550.51025582.0680.2248.3350.0948.3350.0969.9670.81
16552.51025581.9179.0642.9349.4642.9349.4667.5968.09
17401.56.525558.1554.8724.524.3124.524.3166.5671.16
18701.56.525551.0143.0932.843.7332.843.7373.1173.40
19401.51025559.7254.0230.428.7730.428.7770.1964.42
20701.51025590.4380.0851.7048.2051.7048.2066.2266.66
21550.58.253027.4425.0342.3745.1242.3745.1273.2273.43
22552.58.253035.2543.4844.5744.5044.5744.5065.3763.64
23550.58.2548083.4877.3160.1350.5960.1350.5961.2262.21
24552.58.2548093.2895.7655.1749.9655.1749.9667.5166.56
25551.58.2525560.2170.6148.1747.5448.1747.5475.9375.76
26551.58.2525576.9570.6149.5347.5449.5347.5474.9375.76
27551.58.2525572.9470.6148.3047.5448.3047.5475.4175.76
Remark: Independent variables: x1 = Temperature; x2 = E/S ratio; x3 = pH; x4 = hydrolysis time.
Table 3. Of variance at p < 0.05 and regression coefficients for DH (y1), ACE inhibition (y2), DPPH• (y3) and ferrous ion chelating activity (y4), expressed in full quadratic models.
Table 3. Of variance at p < 0.05 and regression coefficients for DH (y1), ACE inhibition (y2), DPPH• (y3) and ferrous ion chelating activity (y4), expressed in full quadratic models.
Factorsy1p-Valuey2p-Valuey3p-Valuey4p-Value
Regression 0.000 0.000 0.000 0.000
Linear 0.000 0.004 0.000 0.000
Squares 0.001 0.000 0.000 0.000
Interaction 0.013 0.000 0.038 0.024
Constant−116.7610.203−57.83670.408−190.3370.000−31.89920.355
x13.4210.110−0.94300.4286.1670.0001.98660.003
x255.4530.01033.14440.007−0.3130.865−5.36410.009
x3−6.2530.51134.29400.0075.7990.02013.90830.041
x40.2190.0000.05420.2100.1580.0270.06090.001
x1*x1−0.0560.001−0.02860.003−0.0500.000−0.01740.004
x2*x2 NS NS NS NS
x3*x3 NS−3.34720.000 NS−0.95960.022
x4*x40.0000.005−0.00020.000 NS−0.00020.000
x1*x2 NS NS NS NS
x1*x30.3600.0320.47610.000 NS NS
x1*x4 NS0.00140.050 NS NS
x2*x3−5.6030.027−4.08160.006 NS NS
x2*x4 NS NS NS0.01570.024
x3*x4 NS NS−0.0180.038 NS
R291.27% 89.96% 74.62% 84.56%
R2-adj87.40% 83.68% 67.01% 77.70%
Lack of fit 0.665 0.060 0.013 0.027
f-value23.54 14.33 9.80 12.33
NS means not significant (p > 0.05) and thus the respective regression coefficient is not reported.

Share and Cite

MDPI and ACS Style

Hussein, F.A.; Chay, S.Y.; Zarei, M.; Auwal, S.M.; Hamid, A.A.; Wan Ibadullah, W.Z.; Saari, N. Whey Protein Concentrate as a Novel Source of Bifunctional Peptides with Angiotensin-I Converting Enzyme Inhibitory and Antioxidant Properties: RSM Study. Foods 2020, 9, 64. https://doi.org/10.3390/foods9010064

AMA Style

Hussein FA, Chay SY, Zarei M, Auwal SM, Hamid AA, Wan Ibadullah WZ, Saari N. Whey Protein Concentrate as a Novel Source of Bifunctional Peptides with Angiotensin-I Converting Enzyme Inhibitory and Antioxidant Properties: RSM Study. Foods. 2020; 9(1):64. https://doi.org/10.3390/foods9010064

Chicago/Turabian Style

Hussein, Fatima Abdelhameed, Shyan Yea Chay, Mohammad Zarei, Shehu Muhammad Auwal, Azizah Abdul Hamid, Wan Zunairah Wan Ibadullah, and Nazamid Saari. 2020. "Whey Protein Concentrate as a Novel Source of Bifunctional Peptides with Angiotensin-I Converting Enzyme Inhibitory and Antioxidant Properties: RSM Study" Foods 9, no. 1: 64. https://doi.org/10.3390/foods9010064

APA Style

Hussein, F. A., Chay, S. Y., Zarei, M., Auwal, S. M., Hamid, A. A., Wan Ibadullah, W. Z., & Saari, N. (2020). Whey Protein Concentrate as a Novel Source of Bifunctional Peptides with Angiotensin-I Converting Enzyme Inhibitory and Antioxidant Properties: RSM Study. Foods, 9(1), 64. https://doi.org/10.3390/foods9010064

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop