Systematic interpolation method predicts protein chromatographic elution with salt gradients,pH gradients and combined salt/pH gradients

A methodology is presented to predict protein elution behavior from an ion exchange column using both individual or combined pH and salt gradients based on high‐throughput batch isotherm data. The buffer compositions are first optimized to generate linear pH gradients from pH 5.5 to 7 with defined concentrations of sodium chloride. Next, high‐throughput batch isotherm data are collected for a monoclonal antibody on the cation exchange resin POROS XS over a range of protein concentrations, salt concentrations, and solution pH. Finally, a previously developed empirical interpolation (EI) method is extended to describe protein binding as a function of the protein and salt concentration and solution pH without using an explicit isotherm model. The interpolated isotherm data are then used with a lumped kinetic model to predict the protein elution behavior. Experimental results obtained for laboratory scale columns show excellent agreement with the predicted elution curves for both individual or combined pH and salt gradients at protein loads up to 45 mg/mL of column. Numerical studies show that the model predictions are robust as long as the isotherm data cover the range of mobile phase compositions where the protein actually elutes from the column.     

Systematic interpolation method predicts protein chromatographic elution from batch isotherm data without a detailed mechanistic isotherm model

Predicting protein elution for overloaded ion exchange columns requires models capable of describing protein binding over broad ranges of protein and salt concentrations. Although approximate mechanistic models are available, they do not always have the accuracy needed for precise predictions. The aim of this work is to develop a method to predict protein chromatographic behavior from batch isotherm data without relying on a mechanistic model. The method uses a systematic empirical interpolation (EI) scheme coupled with a lumped kinetic model with rate parameters determined from HETP measurements for non‐binding conditions, to numerically predict the column behavior. For two experimental systems considered in this work, predictions based on the EI scheme are in excellent agreement with experimental elution profiles under highly overloaded conditions without using any adjustable parameters. A qualitative study of the sensitivity of predicting protein elution profiles to the precision, granularity, and extent of the batch adsorption data shows that the EI scheme is relatively insensitive to the properties of the dataset used, requiring only that the experimental ranges of protein and salt concentrations overlap those under which the protein actually elutes from the column and possess a ±10% measurement precision.     

Evaluation of differences between dual salt‐pH gradient elution and mono gradient elution using a thermodynamic model: Simultaneous separation of six monoclonal antibody charge and size variants on preparative‐scale ion exchange chromatographic resin

The efficiencies of mono gradient elution and dual salt‐pH gradient elution for separation of six mAb charge and size variants on a preparative‐scale ion exchange chromatographic resin are compared in this study. Results showed that opposite dual salt‐pH gradient elution with increasing pH gradient and simultaneously decreasing salt gradient is best suited for the separation of these mAb charge and size variants on Eshmuno^® CPX. Besides giving high binding capacity, this type of opposite dual salt‐pH gradient also provides better resolved mAb variant peaks and lower conductivity in the elution pools compared to single pH or salt gradients. To have a mechanistic understanding of the differences in mAb variants retention behaviors of mono pH gradient, parallel dual salt‐pH gradient, and opposite dual salt‐pH gradient, a linear gradient elution model was used. After determining the model parameters using the linear gradient elution model, 2D plots were used to show the pH and salt dependencies of the reciprocals of distribution coefficient, equilibrium constant, and effective ionic capacity of the mAb variants in these gradient elution systems. Comparison of the 2D plots indicated that the advantage of opposite dual salt‐pH gradient system with increasing pH gradient and simultaneously decreasing salt gradient is the noncontinuous increased acceleration of protein migration. Furthermore, the fitted model parameters can be used for the prediction and optimization of mAb variants separation in dual salt‐pH gradient and step elution. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 34:973–986, 2018     

Apolipoprotein A-I(Milano) anion exchange chromatography: Self association and adsorption equilibrium

The self-associative properties of apolipoprotein A-I(Milano) (apoA-I(M)) were investigated in relationship to its anion exchange behavior on Q-Sepharose-HP with and without the addition of urea as a denaturant. Self-association was dependent on protein and urea concentration and both influenced interactions of the protein with the chromatographic surface. In the absence of urea, apoA-I(M) was highly associated and existed primarily as a mixture of homodimer, tetramer and hexamer forms. Under these conditions, since the binding strength was greater for the oligomer forms, broad, asymmetrical peaks were obtained in both isocratic and gradient elution. Adding urea depressed self-association and caused unfolding. This resulted in sharper peaks but also decreased the binding strength. Thus, under these conditions chromatographic elution occurred at lower salt concentrations. The adsorption isotherms obtained at high protein loadings were also influenced by self-association and by the varying binding strength of the differently associated and unfolded forms. The isotherms were thus dependent on protein, urea, and salt concentration. Maximum binding capacity was obtained in the absence of urea, where adsorption of oligomers was shown to be dominant. Adding urea reduced the apparent binding capacity and weakened the apparent binding strength. A steric mass action model accounting for competitive binding of the multiple associated forms was used to successfully describe the equilibrium binding behavior using parameters determined from isocratic elution and isotherm experiments.     

Mathematical modeling of elution curves for a protein mixture in ion exchange chromatography applied to high protein concentration

Protein elution curves in ion exchange chromatography (IEC) were simulated with a rate model. Three pure proteins and their mixture were used (α‐lactalbumin, BSA, and conalbumin) under different operational conditions. The anionic matrix Q‐Sepharose FF was used packed in a 1 mL column. A high protein concentration (37.5 mg/mL of total protein injected into the column) was used in order to extend the utility of the model. Mass transfer parameters were calculated using empiric correlations, where the axial dispersion was negligible (Pe > 300) and the mass transfer was controlled by the intraparticle diffusion (Bi > 10). The model assumes a modulator–eluite relationship were the equilibrium constant of the Langmuir isotherm was a function of salt concentration. Adsorption kinetic parameters were estimated from experimental data. The parameters for pure proteins were determined, and elution curves for changes in flow rate, ionic strength gradient, concentration, and sample size were predicted by the model. Then the kinetic parameters of the mixture were determined under the same operational conditions and some of the parameters had to be modified to take into account effects such as protein–protein interactions, competition, and displacement. Experimental elution curves obtained for changes in operational conditions such as flow rate and ionic strength gradient were simulated by the rate model for the protein mixture with a relative error in retention time of visible peaks <5%. IEC operational conditions and the peak fraction collection can be selected using a cost function of the production process which considers yield, purity, concentration, and process time that are obtained from simulations. Operational conditions that gave the minimum cost were selected. Simulations allows to diminish experimental time and cost. Biotechnol. Bioeng. 2009; 104: 572–581 © 2009 Wiley Periodicals, Inc.     

Influence of pH and ionic strength on the steric mass-action model parameters around the isoelectric point of protein

The ion-exchange equilibrium and the dependence of the parameters in the steric mass-action (SMA) model on salt concentration and buffer pH around the isoelectric point of protein were studied. Bovine serum albumin (BSA, isoelectric point = 5.4) was used as a model protein and DEAE Sepharose FF as an ion exchanger. Finite batch adsorption experiments and isocratic elution chromatography were performed for the determination of the model parameters (i.e., characteristic charge, equilibrium constant, and steric factor). The results showed that pH had significant effects on the parameters. With an increase of pH from 4.5 to 6.5, the characteristic charge increased from 0.9 to 3.0 and leveled off as a plateau at pH above 5.5. The charge groups in the contact region of protein surface were considered to play a crucial role on the characteristic charge. The decrease of pH and increase of salt concentration lowered the absolute value of the zeta potential of the protein surface and led to a decrease of the equilibrium constant. The steric factor remained unchanged at about 31 at pH 5.5 and 6.0 and increased to 44.5 at pH 5.0 and 96.8 at pH 4.5, mainly as a result of the lower adsorption capacity of BSA at pH <5.5. Furthermore, the increase of the molecular volume of BSA at pH 4.5 would be an additional reason for the increase of the steric factor. Taking into account the effect of the pH and salt concentration on these parameters, the SMA model described the ion exchange equilibrium of protein more accurately.     

Adsorption-desorption of BSA to highly substituted dye-ligand adsorbent: quantitative study of the effect of ionic strength

Cibacron Blue 3GA was immobilized on Sepharose CL-6B to obtain a highly substituted dye-ligand adsorbent which dye concentration was 17.4?μmol dye per gram wet gel. This adsorbent had a highly binding capacity for bovine serum albumin (BSA). The effects of ionic strength on the adsorption and desorption of BSA to the adsorbent were studied. Adsorption isotherms were expressed by the Langmuir model. The quantitative relationships between the model parameters and the ionic strength were obtained. The desorptions were performed by adding salt to the BSA solutions in which adsorption equilibria had been reached. Adding salt to the solution resulted in the desorption of the bound protein. It was found that the isotherm obtained from the desorption experiments agreed well to the isotherm obtained from the adsorption experiments at the same ionic strength. The result demonstrated that the adsorption of BSA to the highly substituted adsorbent was reversible.     

Cation-exchange displacement chromatography of proteins with protamine displacers: effect of induced salt gradients.

Protamine was investigated for its utility as a protein displacer in cation-exchange systems. Although the protamine solution contained several variants of the molecule, the high affinity of all of the components in this heterogeneous biopolymer enabled it to act as an efficient protein displacer. To facilitate parameter estimation of the protamine, a preliminary purification was carried out by preparative elution chromatography. Chromatographic parameters of both the feed proteins and protamine displacer were obtained for use in a multicomponent steric mass action ion-exchange displacement model. Model simulations were compared to displacement results under both moderate and intense induced salt gradient conditions. In both cases, excellent agreement was obtained between the displacement experiments and theoretical predictions. In addition, these studies serve to dramatize the importance of induced salt gradients in ion-exchange displacement systems.     

A model that links growth and secondary metabolite production in plant cell suspension cultures

Plant cell suspensions of grape cells (Vitis vinifera L. cv. Gamay Fréaux) were grown in shake flasks operated both in the batch and semicontinuous mode. A mathematical model was developed to describe grape cell growth, sucrose uptake, and secondary metabolite (anthocyanin) production. Parameters were estimated from batch studies data. The model was able to predict results for semicontinuous experiments by only modifying the value of four of these parameters. The modified parameters (maximum specific rate of biomass production, maximum specific rate of substrate consumption for maintenance, maximum specific rate of anthocyanin production, and degradation constant of anthocyanins) were related to the kinetics rather than to the yield of the process. The model introduces the concept of primary and secondary metabolism substrate concentration-dependent competition for precursors. Further, the model was able to predict the evolution of the cell system when substrate is scarce, as the value of the different kinetic constants determines the portion of substrate that is used for biomass production, secondary metabolite production, and cell maintenance. (c) 1995 John Wiley & Sons, Inc.     

Protein adsorption on Ion exchange resin: Estimation of equilibrium isotherm parameters from batch kinetic data

The simple Langmuir isotherm is frequently employed to describe the equilibrium behavior of protein adsorption on a wide variety of adsorbents. The two adjustable parameters of the Langmuir isotherm—the saturation capacity, orq _m, and the dissociation constant,K _d—are usually estimated by fitting the isotherm equation to the equilibrium data acquired from batch equilibration experiments. In this study, we have evaluated the possibility of estimatingq _m andK _d for the adsorption of bovine serum albumin to a cation exchanger using batch kinetic data. A rate model predicated on the kinetic form of the Langmuir isotherm, with three adjustable parameters (q _m,K _d, and a rate constant), was fitted to a single kinetic profile. The value ofq _m determined as the result of this approach was quantitatively consistent with theq _m value derived from the traditional batch equilibrium data. However, theK _d value could not be retrieved from the kinetic profile, as the model fit proved insensitive to this parameter. Sensitivity analysis provided significant insight into the identifiability of the three model parameters.     

Hydrophobic interaction adsorption of whey proteins: effect of temperature and salt concentration and thermodynamic analysis

The adsorptive behavior of bovine serum albumin (BSA) and beta-lactoglobulin (beta-lg) on hydrophobic adsorbent was studied at four temperatures and different salt concentrations. The Langmuir model was fitted by experimental equilibrium data showing that an increase in temperature and salt concentration results in an increase on the capacity factor of both proteins. A thermodynamic analysis coupled with isotherm measurements showed that salt concentration and temperature affected the enthalpic and entropic behavior of the adsorption process of both proteins, mainly to the beta-lg. The fast variation in the Z value for temperature over than 303.1K suggest a great conformational change occurring in the beta-lg structure, which almost duplicated the maximum adsorption capacity of this protein.     

Enrichment and characterization of histones by two-dimensional hydroxyapatite/reversed-phase liquid chromatography-mass spectrometry

Here we report a novel two-dimensional liquid chromatography-mass spectrometry (2D LC-MS) method that combines offline hydroxyapatite (HA) chromatography with online reversed-phase liquid chromatography-mass spectrometry (HA/RP LC-MS). The 2D LC-MS method was used to enrich and characterize histones and their posttranslational modifications. The 2D HA/RP LC-MS approach separates histones based on their relative binding affinity to DNA and relative hydrophobicity. HA/RP separations showed improvement in the number of histone isoforms detected as compared with one-dimensional RP LC-MS of acid-extracted histones. The improved histone fractionation resulted in better detection of lower abundant histone variants as well as their posttranslationally modified isoforms. Histones eluted from the HA/RP in the following order: H1, H2A/H2B heterodimers followed by H3/H4 heterotetramers, as predicted from their spatial organization in nucleosomes for binding affinity to DNA. Comparison between HA-purified and acid-extracted histones revealed similar histone profiles with the exception that the HA fractions showed a greater number of H1 isoforms. Two elution conditions were also examined: batch elution and salt gradient elution. Although both elution techniques were able to fractionate the histones sufficiently, the salt gradient approach has the most potential for purification of selected histone isoforms.     

Hydrophobic interaction adsorption of hen egg white proteins albumin, conalbumin, and lysozyme

Hydrophobic adsorption equilibrium data of the hen egg white proteins albumin, conalbumin, and lysozyme were obtained in batch systems, at 25 degrees C, using the Streamline Phenyl resin as adsorbent. The influence of three types of salt, NaCl, Na(2)SO(4), or (NH(4))(2)SO(4), and their concentration on the equilibrium data were evaluated. The salt Na(2)SO(4) showed the higher interaction with the studied proteins, thus favoring the adsorption of proteins by the adsorbent, even though each type of salt interacted in a distinct manner with each protein. The isotherm models of Langmuir, Langmuir exponential, and Chen and Sun were well fitted to the equilibrium data, with no significant difference being observed at the 5% level of significance. The mass transfer model applied simulated correctly adsorption kinetics of the proteins under the studied conditions.     

Purification of recombinant α-amylase with immuno-affinity chromatography using monoclonal antibody

A monoclonal antibody against recombinant thermostable α-amylase produced by Escherichia coli was isolated from serum-free medium and immobilized on Sepharose 4B. The adsorption equilibrium between α-amylase and the immobilized immuno-adsorbent showed a Langmuir type isotherm. The breakthrough curve calculated numerically using the averaged volumetric coefficient coincided well with the experimental data. More than 90% of the activity of bound α-amylase could be recovered by eluting with glycine-HCl buffer (pH 2.5). The elution profile at pH 2.5 became sharper with increasing temperature. By using an immuno-affinity column, the recombinant α-amylase produced by E. coli could be purified homogeneously from crude extract enzyme solution with two-step elution.     

Role of protein conformation and aggregation in pumping water in and out of a cell

Dialysis cassettes containing BSA solutions were used to simulate passive in vivo conditions to assess the effect of protein conformation and aggregation on cell water content. The cassettes were suspended in dextran solutions to provide a range of fixed osmotic stress values simulating blood plasma. The system was placed on a shaker for 24 h to attain equilibrium. Four manipulation methods; pH, cosolute salt concentration, e.g. NaCl, temperature annealing and urea concentration denaturant were varied to produce well-known manipulations of BSA conformation. It was observed that the cell water content varied from +14% to about -13% with changes in protein conformation and aggregation. The findings demonstrate that a change in protein conformation and aggregation, pumps water in and out of a cell to maintain equilibrium % water content matching the protein conformational hydration parameter. This concept supplements existing theories on cell volume regulation.     

Optimum Sorption Isotherm by Linear and Nonlinear Methods for Safranin onto Alkali-Treated Rice Husk

Rice husk, a lignocellulosic by-product of the agroindustry, was treated with alkali and used as a low-cost adsorbent for the removal of safranin from aqueous solution in batch adsorption procedure. In order to estimate the equilibrium parameters, the equilibrium adsorption data were analyzed using the following two-parameter isotherms: Freundlich, Langmuir, and Temkin. A comparison of linear and nonlinear regression methods in selecting the optimum adsorption isotherm was applied on the experimental data. Six linearized isotherm models (including four linearized Langmuir models) and three nonlinear isotherm models are thus discussed in this paper. In order to determine the best-fit isotherm predicted by each method, seven error functions namely, coefficient of determination (r ²), the sum of the squares of the errors (SSE), sum of the absolute errors (SAE), average relative error (ARE), hybrid fractional error-function (HYBRID), Marquardt's percent standard deviation (MPSD), and the chi-square test (χ²) were used. It was concluded that the nonlinear method is a better way to obtain the isotherm parameters and the data were in good agreement with the Langmuir isotherm model.     

Role of biotic interactions and physical factors in determining the distribution of marsh species along an estuarine salinity gradient

Understanding the mechanisms that shape plant distribution patterns is a major goal in ecology. We investigated the role of biotic interactions (competition and facilitation) and abiotic factors in creating horizontal plant zonation along salinity gradients in the Elbe estuary. We conducted reciprocal transplant experiments with four dominant species from salt and tidal freshwater marshes at two tidal elevations. Ten individuals of each species were transplanted as sods to the opposing marsh type and within their native marsh (two sites each). Transplants were placed at the centre of 9‐m² plots along a line parallel to the river bank. In order to disentangle abiotic and biotic influences, we set up plots with and without neighbouring vegetation, resulting in five replicates per site. Freshwater species (Bolboschoenus maritimus and Phragmites australis) transplanted to salt marshes performed poorly regardless of whether neighbouring vegetation was present or not, although 50–70% of the transplants did survive. Growth of Phragmites transplants was impaired also by competition in freshwater marshes. Salt marsh species (Spartina anglica and Puccinellia maritima) had extremely low biomass when transplanted to freshwater marshes and 80–100% died in the presence of neighbours. Without neighbours, biomass of salt marsh species in freshwater marshes was similar to or higher than that in salt marshes. Our results indicate that salt marsh species are precluded from freshwater marshes by competition, whereas freshwater species are excluded from salt marshes by physical stress. Thus, our study provides the first experimental evidence from a European estuary for the general theory that species boundaries along environmental gradients are determined by physical factors towards the harsh end and by competitive ability towards the benign end of the gradient. We generally found no significant impact of competition in salt marshes, indicating a shift in the importance of competition along the estuarine gradient.     

Ion exchange chromatography of proteins-predictions of elution curves and operating conditions. II. Experimental verification

The applicability and validity of the model developed in Part I were confirmed experimentally. In this article, various proteins were eluted both by stepwise and linear gradient elution on DEAE ion exchangers under a variety of experimental conditions. Adsorption isotherms were measured as a function of ionic strength in batch experiments. The moment method was empolyed for the determination of various parameteres such as the gel-phase diffusion coefficient and the longitudinal dispersion coefficient. By use of these parameters and the experimentally measured ionic strength of the peak position, the number opf plates was determined according to the method described in Part I. Theoretical elution curves were calculated with the experimentally measured adsorption eqluilibria and the number of plates. Good agreement was observed between theory an experiments. Various factors affecting the separation were investigated. It was found that the effect of the number of plates for salts, N'(p), was negligible except the case of stepwise elution of high ionic strength buffer. When elution curves were symmetrical, the widths of the elution curves were inversely proportional to the square root of the number of plates of proteins, N(p), as in other chromatographic techniques. A simple graphical method for prediction of the peak position in linear gradient elution described in Part I was found applicable when the elution curves were symmetrical. A useful correlation of prediction of the peak width in a linear gradient elution was proposed on the basis of the approximate solution derived in Part I of this study. This graphical method and correlation permit easy prediction of the peak position and peak width in linear gradient elution in the case of symmetrical elution curves.