Optimization of Lactoperoxidase and Lactoferrin Separation on an Ion-Exchange Chromatography Step

Lactoperoxidase (LP), which is a high-value minor whey protein, has recently drawn extensive attention from research scientists and industry due to its multi-function and potential therapeutic applications. In this study, the separation and optimization of two similar-sized proteins, LP and lactoferrin (LF) were investigated using strong cation exchange column chromatography. Optimization was started with central composite design based experiments to characterize the importance of different decision variables. The three variables used in the optimization were flow rate, length of gradient and final salt concentration in the linear elution gradient step. The obtained empiric functional model represented the effect of the significant factors on the yield as the objective function. Afterwards, the calibrated mechanistic model was employed to predict accurate optimal set of variables. The optimal operating points were found and the results were compared with validation experiments. Predictions respecting yield confirmed a very good agreement with experimental results while keeping purity, a product quality characteristic, equal or above to a predefined value.


Introduction
Nowadays, the extraction and purification of high-value proteins from waste food streams, which can be re-used as nutritionally and pharmacologically functional ingredients has attracted great attention [1].Whey, the liquid residue from dairy processing, is one of the biggest reservoirs of food protein available today.World whey output was approximately 180 million tonnes in 2013, containing about 1.5 million tonnes of high-value proteins and 8.6 million tonnes of lactose [2].Whey proteins are not only nutritious supplements, but also valuable bioproducts demonstrating a wide range of biological features, such as antivirus properties, anticarcinogenic activities and immune system modulation [3].It is therefore of utmost importance to develop cost-effective techniques for large-scale isolation of pure whey protein components.
Whey proteins have been adequately separated into different fractions, the major and minor proteins [4][5][6][7][8].Isolation and purification of major proteins, e.g., β-lactoglobulin, α-Lactalbumin and serum albumin, through various chromatographic and membrane-separation techniques have been extensively studied [7,9,10].However, the efficient purification of high-value minor proteins such as lactoperoxidase (LP) and lactoferrin (LF) still remains a challenge due to their extremely low concentrations.LP, a peroxidase enzyme with molecular weight of approximately 77.5 kDa [11] and an isoelectric point approaching 9.5, contributes to the humoral immune defense against pathogens by oxidation of thiocyanate and iodide [12].LF, an iron-binding glycoprotein with a molecular weight of about 80.0 kDa and an isoelectric point around 8.7, shows many biological functions, including: antibacterial, antitumor and anti-inflammatory functions, as well as anabolic functions for bone [13].Both LP and LF have great therapeutic potential and are expected to find many important applications.Over the past two decades, research efforts have been invested in isolating LP and LF by different chromatographic methods [8,14,15].However, the average recovery rates and final purity for LP and LF were quite low [16], and technical development to improve the yield and purity of LP and LF deserves further study.
Ion-exchange chromatography (IEC) is one of the most powerful and widely-used techniques for isolation and purification of LP and LF [8,17].Nevertheless, the similar molecular weights and isoelectric points of lactoperoxidase and lactoferrin pose difficulties in the isolation of LP and/or LF in a pure form by ICE.To improve the resolution, either salt-gradient or pH-gradient elution can be employed [18], which complicates the design and the operation of an IEC process.Optimization of the operating parameters, such as protein loading, elution flow rate, pH, salt concentration and gradient length, necessitates a large number of experiments, which are expensive and time-consuming.Optimization of LP and LF separation by IEC based on experimentally-verified mechanistic chromatographic model is a better choice because of its high efficiency and accuracy.
Successful approaches to the optimization of an IEC always include a detailed consideration of the physicochemical properties of involved components as well as interactions between the proteins and the adsorbent phase.To the best of our knowledge, no model-integrated approach has been applied to optimizing minor milk protein separation in IEC, making it a topic worthy of investigation.The present work aims to develop a two-step optimization strategy to maximize yield of lactoperoxidase while respecting constraints on purity.A statistical model based on response surface methodology was first employed to examine the influence of the operation variables (concentration of salt at the end of gradient, length of gradient, and velocity).Subsequently, productive information was used to assist the mechanistic model to gain more insight into the optimization of IEC.Methods developed in this study provide a reasonable way to solve complex high-dimensional optimization problems, which can be applied for other chromatographic separations of proteins.

Response Surface Modeling and Design of Experiments
Response surface methodology (RSM), a collection of mathematical-statistical techniques based on design of experiments (DoE), has been successfully used for optimization studies of different bio-separation processes [19].RSM has been widely adopted to investigate the effects of several design factors influencing a response by varying them simultaneously in a limited set of experiments.The concept of DoE-techniques outlines various statistical approaches to maximize specific information in an experimental planning and after all determine the most favorable direction to move in order to find a true optimum.Central composite design (CCD) is an ideal choice as a symmetrical experimental design for sequential experimentation, and allows reasonable information to calculate model lack of fit while reducing the number of design points.In general, CCD is the better-known class of quadratic design that consists of: (1) a factorial (or cubic) design; (2) an additional design, often a star design, with all points set equidistantly from the center; and (3) at least one center point [20].Therefore, full uniformly-routable central composite designs involve the total number of N = k 2 + 2k + c p points, where k is the number of factors, and c p is the number of replicate runs performed at the center point.In this design, the distance α from the center point depends on the number of the factors and can be calculated by α = 2 (k−p)/4 .All factors have to be adjusted at five levels (−α, −1, 0, +1, +α) [20].The quadratic polynomial model for the measured values of the results from experiments variable, Y with k factors is given by: where x i and x j are the design factors in coded values, β 0 is the constant parameter, β i , β ij , and β ii are the coefficient of the linear, interaction, and quadratic terms of the model, respectively.The coefficients of Equation (1) are estimated using statistical software packages (e.g., Minitab, Design Expert, SPSS).

Mechanistic Modeling of Chromatography
A mechanistic model is used to describe the physical phenomena based on a set of mathematical equations.Two types of physical phenomena dominate chromatography: movement of solutes through the packed bed of porous particles via mass transfer mechanisms, and adsorption based on the fundamental thermodynamic interactions between migrating solutes and the stationary phase.The general system of equations used to describe the mass transfer phenomena consist of two sets of partial differential mass conservation equations.The general rate model for a chromatographic process includes convective and diffusive flows through porous particles on the column level and imitates mass transfer resistances and surface interactions on particle level.In IEC, an external film surrounding adsorbent particles is commonly presumed to model the movement of components from column to particle level; the sorption of protein on the particle surface can be described by the steric mass-action (SMA) model, developed by Iyer et al. [21] and generally used for the modeling of salt gradient elution in IEC, for example in [22].
On column level, concentration change for the ith component with respect to the time and position, is described by: The first term on the right-hand side represents the convective transport through the column while the second and third terms represent respectively the dispersive transport and the mass transfer to the particle surface.The symbol u int indicates the interstitial velocity, ε c the column voidage, r p the particle radius, D ax the axial dispersion representing combined effect of dispersion and diffusive processes, and k eff,i epitomize combined effect of both the internal and external mass transfer resistances in one lumped film diffusion coefficient.Analogously, on a particle level, concentration change for the ith component is expressed by: where q i denotes the concentration of component i within the particle and ε p the particle voidage.
The first term on right-hand side describes adsorption and desorption processes on particle level, i.e., the interaction between mobile and particle-bound phase.
For the description of sorption kinetics, the steric mass-action (SMA) isotherm developed by Brooks and Cramer [23] was embedded into the mechanistic model.It has been successfully used to describe protein adsorption in IEC.Based upon the stoichiometric exchange of charges and steric hindrance of binding sites, the isotherm can be described by the reaction: The parameter ν i is the characteristic charge of protein, which represents the average number of ligands, or binding sites, interacting during adsorption.q salt is the concentration of adsorbed salt counter-ions that are available for exchange, c salt is the salt concentration of the bulk.
where the parameter Λ is the ionic capacity of the adsorbent and σ i is the steric factor, which represents the average number of counter-ions shielded per adsorbed protein molecule.
SMA isotherm in its kinetic form can be expressed as: k ads,i and k des,i respectively denote the adsorption and desorption rates.

Optimization for the ICE Separation of LP and LF
Optimization maximizes/minimizes an objective function by varying one or more variables, called decision variables, to obtain an optimal purification step.A chromatographic purification step involves many decision variables, such as loading-, washing-and elution times, salt concentration, flow rate for different steps, column length, column diameter, and gradient for the elution step.The decision variables can be narrowed down by including lower and upper boundaries.To limit the optimization further, constraints can also be added as equality or inequality functions.The most typical objective functions used in a preparative chromatographic purification step are production rate and yield.In this work, the objective function studied was the yield, and the inequality constraints set on the optimization was purity.
The yield is calculated as the fraction of the target component captured: where c captured,i is the concentration of the target component i leaving the column during the elution.The purity is defined as: where m captured,i is the amount of substance captured of component i.

Materials, Column and Software
The study aims at an optimal separation of lactoperoxidase on the adsorbent SP Sepharose FF by linear gradients.Lactoperoxidase from bovine milk, ≥90% pure, and Lactoferrin from bovine milk, approximately 95% pure, were used.Sodium monobasic phosphate, sodium dibasic phosphate and sodium chloride for buffers preparation were purchased from Sigma-Aldrich (Oakville, ON, Canada).NaOH was used for pH-adjustment.The running buffer in all experiments was a 0.02 M sodium phosphate buffer at pH 6.7.The buffer for elution purposes contained additional NaCl.Proteins were diluted into phosphate buffers by slowly stirring to prevent any foam formation.Afterward, the protein solution was filtered through a 0.2 µm hydrophilic polypropylene membrane filter to remove any fine particles.The chromatographic setup consisted of a prepacked HiScreen™ SP Sepharose FF 4.7 mL column (10.0 cm length, 0.77 cm ID) and an ÄKTA purifier 100 system, both purchased from GE Healthcare (Mississauga, ON, Canada).The software Minitab ® 17 (State College, PA, USA) was used as a statistical tool for handling response surface methodology.The software MATLAB R2014a (Natik, MA, USA) was used to execute the mechanistic model.

Experimental Methods
In all experimental setups the column was at first equilibrated with running buffer for 5 column volumes (CV).This step was followed by an automated sample load of 2 mL protein mixture.Then the column was flushed for another 2 CV to remove unbound proteins, before initiating a linear elution gradient.The elution gradient was applied from 0% to 100% high-salt elution buffer, followed by a 5 CV high-salt wash step, and regenerated and re-equilibrated with 1 M NaOH and running buffer respectively.Conductivity and UV-absorbance were measured online at column outlet.The data collected from these measurements was further analyzed and taken into account to estimate SMA parameter by the inverse method.Three design factors were employed to describe the gradient profile: The experimental ranges used for these factors included (0.35-1.35 M) for final concentration of salt in elution step, (1.11-5.98 cm/min) for flow velocity, and (9.88-35.1 CV) for length of gradient.All the buffers were prepared using ultra-pure water (Barnstead easy-pure RODI, Fisher Scientific), filtered with a 0.45 µm membrane and degassed prior to use.

Protein Quantification
Lactoperoxidase concentration was determined by Lowry protein assay.The collected fractions of lactoperoxidase were analyzed using measurement of absorbance at 280 nm (extinction coefficient of lactoperoxidase ε 1% = 14.9 [24]) and 412 nm.Lactoperoxidase absorbs radiation at 280 nm as well as 412 nm; it has maximum absorbance at 412 nm [25] and its purity is estimated as a ratio of A 412 /A 280 .The BCA (bicinchoninic acid) protein assay kit from Sigma (Oakville, ON, Canada) with BSA protein standard was used for further analysis of protein fractions from ÄKTA.

Screening Experiments to Determine Importance of Design Factors
Concentration of salt was studied over a broad range, as for the elution step a lower salt concentration can result in a lower purity for lactoperoxidase, and less lactoperoxidase will elute in the elution step.On the other hand, a higher salt concentration in the elution step is less critical, as lactoferrin is far from eluting at lower conductivities.Flow rate is remarked to be easy to control and was investigated over a narrower range, with an upper bound recommended by the manufacturer.Initial operating conditions to maximize chromatography performance were determined by means of a functional relationship between the experimental designs combined with response surface modeling.To ensure the consistency of prediction error, the value (α) was adjusted as 1.68 ( 4 √ 8), and it was assumed that all points with an equal distance from the design center have a constant prediction variance.Table 1 shows design factors in their coded and un-coded (experimental) forms in an experimental design performed in a random order to avoid systematic error.As seen in Table 1, the design experiments had 20 runs in total.The analysis of variance (ANOVA) was used to test the significance of the fit of the empiric model.The significance level of 0.05 was chosen to establish the statistical significance in all results.The inverse method [26] was used to calibrate the model parameters to experimental data.SMA parameters were estimated based on 20 chromatograms seeking to achieve a best fit between model response and measured chromatogram data.The optimal set of the parameter values, {ν, k ads , k des , σ}, can be generated by minimizing the error function F(p), defined as follows: where p is the total number of the experimental data.The minimization of Equation ( 9) was performed with the MATLAB function lsqnonlin.Method of lines (MoL) [27] was used to discretize the column in spatial dimensions.The boundary conditions of the column were Danckwert's boundary conditions [28].The discretized model was solved with a stiff ordinary differential equation solver (ode15s) using variable-order method in MATLAB [29].
For estimating the dispersion coefficient's dependency on flow rate, the mean particle Péclet number (Pe), estimated to be 0.5 [30], the dispersion coefficient was calculated from: The film mass transfer coefficient k f can be estimated from the correlation [31]: where Sh = k eff •d p /D m is the Sherwood number, Re = u 0 ρd p /η is the Reynolds number and Sc = η/ρD m is the Schmidt number.
To ascertain that the model is correct, the confidence interval of the parameters calculated from the Jacobian and the residual.The relative residual error is plotted against the cumulative distribution to show how well the model fits the data.The model is then compared to the validation experiments, and if the simulation fits the experiments, the model is valid.

Optimization Method
The calibrated mechanistic model was then employed for optimization with respect to an objective function.In this study, yield as the fraction of target protein eluted was defined the objective function.The purity of lactoperoxidase was used as the nonlinear inequality constraint with a requirement of 85%.Parameters such as concentration of feed and pH are determined by the composition of natural milk.Other parameters, such as buffer, eluting salt and stationary phase material, were kept the same during the optimization step.The main decision variables were determined according to the results of the DoE-RSM; length of gradient has no significant effect on yield, thus the response to variation of the final salt concentration in the elution and mobile phase flow velocity were studied by computer simulation.The optimization problem with respect to the objective was solved using fmincon in MATLAB.fmincon finds the optimum of an objective function Y with defined lower and upper boundaries on the decision variables.The optimization problem can be defined as: Decision variables, x = [velocity, salt concentration at the end of gradient] LB ≤ x ≤ UB subject to purity ≥ purity requirement.

Results of Response Surface Modeling
In this study, the best-fitting response surface provided a functional relationship between the objective function yield and the gradient explanatory factors salt concentration, length of gradient, and flow velocity.The normal probability plot of the residuals from the analysis was normally distributed, indicating no evidence of non-normality, skewedness, outliers, or the existence of unidentified variables.An R 2 of 0.89 is probably due to the variations in the experiments at the six center points.
According to the analysis of variance (ANOVA), the calculated probability (p-value) of a test statistic was less than 0.005 for the term final concentration of salt which translates its significance within a 95% confidence interval; whereas the p-values of gradient length and flow velocity justifiy that the coefficients are zero.The mixed effects/interaction terms of each two factors were also quantified and revealed no considerable effect in the response model.The coefficient plot (Figure 1) shows the scaled and centered coefficients for the most important factors and their influence on the objective function.The height and direction of the bars illustrate the corresponding significance of each factor.In addition, the coefficient plot displays relatively large confidence intervals for the coefficients.The broad confidence intervals are probably due to variations in the experimental equipment as the six center points show fairly large variance.Another explanation could be the fact that the investigated system reveals non-linear performance that the statistical model cannot describe.quantified and revealed no considerable effect in the response model.The coefficient plot (Figure 1) shows the scaled and centered coefficients for the most important factors and their influence on the objective function.The height and direction of the bars illustrate the corresponding significance of each factor.In addition, the coefficient plot displays relatively large confidence intervals for the coefficients.The broad confidence intervals are probably due to variations in the experimental equipment as the six center points show fairly large variance.Another explanation could be the fact that the investigated system reveals non-linear performance that the statistical model cannot describe.Analysis of the model response shows that gradient final salt concentration is an important factor for yield.This is to be expected as salt will displace protein as elution proceeds, and as a result Analysis of the model response shows that gradient final salt concentration is an important factor for yield.This is to be expected as salt will displace protein as elution proceeds, and as a result more protein will be present in the collected fractions.It suggests that gradients with a high salt concentration at gradient end were most successful with respect to the separation problem.
The variation in the two factors gradient length and flow velocity were not high enough to cause an effect on the separation of lactoperoxidase in the process.It assures that at higher flow rates, there is no broadening effect or leakage in the loading step.The significance of the coefficient for velocity × velocity term with regard to the p-value predicted by the model indicates that the hypersurface exhibits a curvature meaning that there is a maximum or minimum somewhere in the direction of flow velocity (Figure 2).
Separations 2017, 4, 10 8 of 14 more protein will be present in the collected fractions.It suggests that gradients with a high salt concentration at gradient end were most successful with respect to the separation problem.The variation in the two factors gradient length and flow velocity were not high enough to cause an effect on the separation of lactoperoxidase in the process.It assures that at higher flow rates, there is no broadening effect or leakage in the loading step.The significance of the coefficient for velocity × velocity term with regard to the p-value predicted by the model indicates that the hypersurface exhibits a curvature meaning that there is a maximum or minimum somewhere in the direction of flow velocity (Figure 2). Figure 3 illustrates the projection of the response surface as a two-dimensional plane for three levels of the factor length of gradient: 15, 22.5 and 30 CV.As can be seen in Figure 3, the maximum velocity of around 3.8 cm/min with highest concentration of salt in elution buffer is predicted to achieve maximum protein yield to for all gradient lengths.The optimal region is quite large and the concentration of salt within the contour plot is steep.This indicates that small changes in the concentration of salt in elution buffer will have significant effects on the product quality of the investigated process.
Based on this empiric model function and the modeling surfaces, the factor setups for maximum qualities of separation with respect to yield were predicted (see Table 2) and experiments to evaluate the results predicted by the model were performed.The comparison between experimental results Figure 3 illustrates the projection of the response surface as a two-dimensional plane for three levels of the factor length of gradient: 15, 22.5 and 30 CV.As can be seen in Figure 3, the maximum velocity of around 3.8 cm/min with highest concentration of salt in elution buffer is predicted to Separations 2017, 4, 10 9 of 15 achieve maximum protein yield to for all gradient lengths.The optimal region is quite large and the concentration of salt within the contour plot is steep.This indicates that small changes in the concentration of salt in elution buffer will have significant effects on the product quality of the investigated process.Based on this empiric model function and the modeling surfaces, the factor setups for maximum qualities of separation with respect to yield were predicted (see Table 2) and experiments to evaluate the results predicted by the model were performed.The comparison between experimental results and the empiric model predictions shows low predictability of the RSM model.As can be seen in Figure 4, the predicted coefficient of determination, pred-R 2 of about 0.66, indicates a narrow predictability.In other words, quadratic RSM can only describe 66% of the variety in the experimental data.In this section, the bed parameters of a 4.7 mL prepacked HiScreen™ SP Sepharose FF column characterized in [32] were set into the mechanistic model.The adsorption parameters with their confidence intervals are presented in Table 3.The parameters were estimated based on DoE experiments in section 4.1 to keep the predictions over the similar design space.The steric factor (σ) had negligible influence on the fitting result; thus it was fixed during the optimization of Equation ( 8) to previously calibrated values in [32].Broad statistical range in the determination of σ has been reported before in [21,22].The successful predictions based on the calibrated and verified mechanistic model is a key associated to a model-integrated process development.The validation experiments were carried out as a gradient elution experiment with higher protein concentrations that were not used in the model calibration.The model fitted the validation experiment with relatively good accuracy; in Figure 5 an example chromatogram for the gradient elution experiments with 3.55 cm/min flow velocity, an elution gradient volume of 22.5 CV at 0.26, and 1.47 mg/mL for LP and LF concentrations is shown.The elution gradient was started after 50 min and continued for around 50 min, at which point 100% high-salt buffer with a salt concentration of 0.725 M NaCl was reached.The first breakthrough is LP and the second is LF.In the washing step, some of the LP and LF were washed out; the model did not capture this amount of protein loss, causing the lack of fit in the elution step at the chromatogram.
The cumulative distribution of the error residual, including both the calibration and the validation experiments, is plotted in Figure 6.The 90% figure for the error residual is below a relative error of 4.8% and 99% of the error residual is below 14%, which may seem high; at high protein loading (validation experiments) a small deviation in the retention volume will result in a large residual error.However, these results confirm that the model is capable of describing the problem.If a more accurate solution is sought, the model should be recalibrated based on a new set of DoE experiments in this region.
not capture this amount of protein loss, causing the lack of fit in the elution step at the chromatogram.
The cumulative distribution of the error residual, including both the calibration and the validation experiments, is plotted in Figure 6.The 90% figure for the error residual is below a relative error of 4.8% and 99% of the error residual is below 14%, which may seem high; at high protein loading (validation experiments) a small deviation in the retention volume will result in a large residual error.However, these results confirm that the model is capable of describing the problem.If a more accurate solution is sought, the model should be recalibrated based on a new set of DoE experiments in this region.Cumulative distribution (%)

Optimization Predictions Based on the Mechanistic Model
The elution gradient was numerically optimized with respect to the objective function.On Table 4, the prediction for the optimal gradient based on SMA parameters derived from chromatograms of the 4.7 mL column on the Ä KTA system is given in numbers.To quantitatively evaluate the optimizing performance of the calibrated model, the experimental validation was performed with the corresponding data; chromatogram was transformed beforehand from the UV absorbance measurements on the Ä KTA system to mg/ml of protein using the extinction coefficient.However, the experimental yield is close to the predictions with good precision, and describes the mechanistic model to be successful and predictive.The predictions for the yield of lactoperoxidase with respect to optimal operating conditions deviate slightly; one explanation for this deviation can be the dynamic effects that the model cannot handle, such as pH-variations caused by increasing high-salt buffer in the system as suggested in [23].However, as the constraint on purity (≥85%) was acceptably satisfied in the optimization step, deviation in the yield of protein can be ignored.The lowest expected purity can be determined by means of the confidence intervals of the process parameters.

Optimization Predictions Based on the Mechanistic Model
The elution gradient was numerically optimized with respect to the objective function.On Table 4, the prediction for the optimal gradient based on SMA parameters derived from chromatograms of the 4.7 mL column on the ÄKTA system is given in numbers.To quantitatively evaluate the optimizing performance of the calibrated model, the experimental validation was performed with the corresponding data; chromatogram was transformed beforehand from the UV absorbance measurements on the ÄKTA system to mg/ml of protein using the extinction coefficient.However, the experimental yield is close to the predictions with good precision, and describes the mechanistic model to be successful and predictive.The predictions for the yield of lactoperoxidase with respect to optimal operating conditions deviate slightly; one explanation for this deviation can be the dynamic effects that the model cannot handle, such as pH-variations caused by increasing high-salt buffer in the system as suggested in [23].However, as the constraint on purity (≥85%) was acceptably satisfied in the optimization step, deviation in the yield of protein can be ignored.The lowest expected purity can be determined by means of the confidence intervals of the process parameters.

Conclusions
The method presented in this manuscript constitutes data from factorial design of an experimental approach to quickly identifying the significance of process parameters on the objective function, as well as complexity, of the system.The empirical multi-variate model revealed a low capability to predict the behavior in the IEC column with respect to optimal operating settings for separation as a result of lack of fit and limitations in robustness of the linear gradient elution process.Nevertheless, it was reasonably accurate at proving the significant factors.In the next step, a mechanistic model was calibrated based on the DoE-planned experiments by the inverse method.The calibrated model with a reduced number of variables was then applied to find the optimal operating conditions and provide insight into the knowledge of process performance with respect to yield of lactoperoxidase and more accurate predictions with respect to process variations.The optimal operation was successfully predicted, yet the prediction of yield deviated slightly from the experimental results.This may be related to some effects such as random configurations of protein interacting with binding sites of the adsorbent, or slight changes in pH during the salt gradient elution step that the model was unable to take into account.
In summary, model-integrated process development proved to be efficient with regard to the objective of the optimization of process and led to finding the true optimum process parameters of flow velocity and concentration of salt at the end of gradient.

Figure 1 .
Figure 1.The coefficient plot resulting from the response surface regression of the screening experiments.The values on y-axis represents corresponding values of variables depicted on x-axis.

Figure 1 .
Figure 1.The coefficient plot resulting from the response surface regression of the screening experiments.The values on y-axis represents corresponding values of variables depicted on x-axis.

Figure 2 .
Figure 2. 3D surface plot based on the central composite design (CCD) method.

Figure 2 .
Figure 2. 3D surface plot based on the central composite design (CCD) method.

Figure 5 .
Figure 5.The validation experiment with an elution gradient of 22.5 column values (CV).

Figure 5 .
Figure 5.The validation experiment with an elution gradient of 22.5 column values (CV).

Figure 6 .
Figure 6.The cumulative distribution for the calibration and validation experiments plotted against the relative residual error.

Figure 6 .
Figure 6.The cumulative distribution for the calibration and validation experiments plotted against the relative residual error.

Table 1 .
Coded and un-coded values of the process factors in screening experiments of central composite.

Table 2 .
Optimum factor set for maximum yield of lactoperoxidase based on the design of experiments-response surface methodology (DoE-RSM) approach.

Table 3 .
SMA parameters for lactoperoxidase and lactoferrin on HiScreen™ SP FF column.

Table 4 .
Optimal operating conditions for the separation of LP and LF.

Table 4 .
Optimal operating conditions for the separation of LP and LF.
c p,i concentration of protein i in the pores of the adsorbent (M) c salt salt concentration in the pores of the adsorbent (M) c i protein concentration i in the mobile phase (M) D ax axial dispersion coefficient (mm 2 /s) D m molecular diffusivity in mobile phase (mm 2 /s) k eff,i effective mass transfer coefficient of protein i (mm/s) k ads,i adsorption coefficient of protein i in the SMA isotherm k des,i desorption coefficient of protein i in the SMA isotherm k i steric factor of protein i in the SMA isotherm