Predicting the behaviour of proteins in hydrophobic interaction chromatography. 1: Using the hydrophobic imbalance (HI) to describe their surface amino acid distribution

This paper focuses on the prediction of the dimensionless retention time of proteins (DRT) in hydrophobic interaction chromatography (HIC) by means of mathematical models based on characteristics of the surface hydrophobicity distribution. We introduce a new parameter, called hydrophobic imbalance (HI), obtained from the three-dimensional structure of proteins. This parameter quantifies the displacement of the superficial geometric centre of the protein when the effect of the hydrophobicity of each amino acid is considered. This parameter is simpler and less expensive than those reported previously. We use HI as a way to incorporate information about the surface hydrophobicity distribution in order to improve the prediction of DRT. We tested the performance of our DRT predictive models in a set of 15 proteins. This set includes four proteins whose DRTs are known as very difficult to predict. By means of the variable HI, it was possible to improve the predictive characteristics obtained by models based on the average surface hydrophobicity (ASH) by 9.1%. Also, we studied linear multivariable models based on characteristics determined from the HI. By using this multivariable model, a correlation coefficient of 0.899 was obtained. With this model, we managed to improve the predictive characteristics shown by previous models based on ASH by 31.8%.     

Prediction of retention times of proteins in hydrophobic interaction chromatography using only their amino acid composition

This paper focuses on the prediction of the dimensionless retention time of proteins (DRT) in hydrophobic interaction chromatography (HIC) by means of mathematical models based, essentially, only on aminoacidic composition. The results show that such prediction is indeed possible. Our main contribution was the design of models that predict the DRT using the minimal information concerning a protein: its aminoacidic composition. The performance is similar to that observed in models that use much more sophisticated information such as the three-dimensional structure of proteins. Three models that, in addition to the amino acid composition, use different assumptions about the amino acids tendency to be exposed to the solvent, were evaluated in 12 proteins with known experimental DRT. In all the cases analyzed, the model that obtained the best results was the one based on a linear estimation of the aminoacidic surface composition. The models were adjusted using a collection of 74 vectors of aminoacidic properties plus a set of 6388 vectors derived from these using two mathematical tools: k-means and self-organizing maps (SOM) algorithms. The best vector was generated by the SOM algorithm and was interpreted as a hydrophobicity scale based partly on the tendency of the amino acids to be hidden in proteins. The prediction error (MSE(JK)) obtained by this model was almost 35% smaller than that obtained by the model that supposes that all the amino acids are completely exposed and 40% smaller than that obtained by the model that uses a simple correction factor considering the general tendency of each amino acid to be exposed to the solvent. In fact, the performance of the best model based on the aminoacidic composition was 5% better than that observed in the model based on the three-dimensional structure of proteins.     

Thermodynamics of hydrophobic interaction chromatography of acetyl amino acid methyl esters

Thermodynamic analysis of hydrophobic interaction chromatography of amino acid methyl esters showed entropy-driven adsorption, consistent with solvophobic theory, except for phenyl ester on the Toyopearl resins. All esters adsorbed more strongly to the Toyopearl resins, including the polymethacrylate base matrix, than to Butyl Sepharose. Enthalpy changes were more favorable with the former, explaining the retention difference between Toyopearl Butyl and Butyl Sepharose. An enthalpy change versus heat capacity change plot showed Van der Waals interactions predominantly with the resin matrix. Literature data revealed the same effect for dansylamino acids, shown by isothermodynamic temperature analysis to adsorb more entropically than the esters.     

Influence of nonionic surfactants on the chromatographic behaviour of proteins in hydrophobic interaction chromatography

The chromatographic characteristics of proteins in the presence of additives of nonionic surfactants Brij-35 and Tween-80 in the conditions of descending gradient of ammonium sulfate and phenyl-coated polymeric stationary phase were investigated. It was revealed that retention factors of proteins may be regulated by use of mentioned additives. The improvement of resolution is achieved for some hardly separated pairs of proteins, viz. albumin egg/albumin bovine, aldolase/tripsin. A reversion of the elution order is observed for tripsin/chymotrypsinogen A.     

Effect of surface hydrophobicity distribution on retention of ribonucleases in hydrophobic interaction chromatography

The effect of surface hydrophobicity distribution of proteins on retention in hydrophobic interaction chromatography (HIC) was investigated. Average surface hydrophobicity as well as hydrophobic contact area between protein and matrix were estimated using a classical thermodynamic model. The applicability of the model to predict protein retention in HIC was investigated on ribonucleases with similar average surface hydrophobicity but different surface hydrophobicity distribution. It was shown experimentally that surface hydrophobicity distribution could have an important effect on protein retention in HIC. The parameter "hydrophobic contact area," which comes from the thermodynamic model, was able to represent well the protein retention in HIC with salt gradient elution. Location and size of the hydrophobic patches can therefore have an important effect on protein retention in HIC, and the hydrophobic contact area adequately describes this.     

Specific salt effects in hydrophobic interaction chromatography of proteins

Summary It has been noted in the literature that certain salts enter into specific interaction with proteins. As a result of this, they may act as salting-in agents. We have investigated the effect of magnesium chloride which is known to possess such unusual properties on the retention of proteins in hydrophobic-interaction chromatography. First the retention behaviour of amino acids and small peptides having a wide polarity range was studied on reversed-phase columns using eluents containing (NH₄)₂SO₄, MgSO₄ or MgCl₂, in wide the concentration ranges. For less polar eluites plots of the logarithmic retention factors against the salt concentration were found to be linear, whereas the more polar species showed irregular behavior. The retention of a wide range of proteins was measured on a TSK Phenyl-5-PW column using eluents containing (NH₄)₂SO₄, MgSO₄ or MgCl₂ at different concentrations.The salt-mediated retention was regular with (NH₄)₂ SO₄ and MgSO₄ although MgSO₄ showed a lesser effect than that predicted by the surface-tension increment. The effect of MgCl₂ was quite irregular: the retention factors either increased or decreased or remained unchanged depending on the protein. These results corroborate earlier observations regarding the particular effect of MgCl₂ and suggest the modulation of selectivity in hyrophobic-interaction chromatography by the addition of MgCl₂ to the eluent.     

Multidimensional analysis of denatured milk proteins by hydrophobic interaction chromatography coupled to a dynamic surface tension detector

Multidimensional analysis of denatured milk proteins is reported using high-performance liquid chromatography (HPLC) combined with dynamic surface tension detection (DSTD). A hydrophobic interaction chromatography (HIC) column (a TSK-Gel Phenyl-5PW column, TosoBiosep), in the presence of 3.0 M guanidine hydrochloride (GdmHCl) as denaturing agent is employed as the mobile phase. Dynamic surface tension is measured through the differential pressure across the liquid-air interface of repeatedly growing and detaching drops. Continuous surface tension measurement throughout the entire drop growth (50 ms to 4 s) is achieved, for each eluting drop of 4 s length, providing insight into both the kinetic and thermodynamic behavior of molecular orientation processes at the liquid-air interface. An automated calibration procedure and data analysis method is applied with the DSTD system, which allows two unique solvents to be used, the HIC mobile phase for the sample and a second solvent (water for example) for the standard, permitting real-time dynamic surface tension data to be obtained. Three-dimensional data is obtained, with surface tension as a function of drop time first converted to surface pressure, which is plotted as a function of the chromatographic elution time axis. Experiments were initially performed using flow injection analysis (FIA) with the DSTD system for investigating commercial single standard milk proteins (alpha-lactalbumin, beta-lactoglobulin, alpha-, beta-, kappa-casein and a casein mixture) denatured by GdmHCl. These FIA-DSTD experiments allowed the separation and detection conditions to be optimized for the HIC-DSTD experiments. Thus, the HIC-DSTD system has been optimized and successfully applied to the selective analysis of surface-active casein fractions (alpha s1- and beta-casein) in a commercial casein mixture, raw milk samples (cow's, ewe's and goat's milk) and other diary products (yogurt, stracchino, mozzarella, parmesan cheese and chocolate cream). The different samples were readily distinguished based upon the selectivity provided by the HIC-DSTD method. The selectivity advantage of using DSTD relative to absorbance detection is also demonstrated.     

Hydrophobic contribution of amino acids in peptides measured by hydrophobic interaction chromatography

The adsorption behaviors of amino acids in short chain peptides were examined. Each amino acid, aliphatic or charged, was inserted between the two tryptophans of a peptide, GWWG. The capacity factors of these peptides on an Ocytl-Sepharose column were measured. The adsorption enthalpies, entropies, and the number of repelled water molecules after adsorption were estimated to analyze the contribution of each different amino acid to its hydrophobic adsorption. The peptides inserted with aliphatic amino acids owned the highest capacity factors but released the least amount of adsorption heat among all the peptides under examination. It was found that the hydrophobic contribution of aliphatic amino acids was derived from the entropy gain by repelling the ordered water surrounding them. The insertion of negatively charged amino acids greatly reduced the capacity factors but still repelled a significant number of water molecules after adsorption. This indicated that the water molecules surrounding ionic amino acids were not orderly aligned. The dehydration cost energy but the water repelling did not offer enough entropy to drive the adsorption. Subsequently, lower retention was obtained from the peptides inserted with negatively charged ionic amino acids. The insertion of lysine increased the adsorption enthalpy but repelled no water molecules after adsorption. It was speculated that the inserted lysine still interacted with hydrophobic ligands but disturbed the interaction between ligands and adjacent tryptophans. Therefore, the adsorption enthalpy increased and the capacity factors decreased. Different amino acids contributed to hydrophobic interaction in different ways. The simultaneous analysis of capacity factor, adsorption enthalpy, adsorption entropy, and the number of repelled water molecules facilitated the understanding of the adsorption processes.     

Use of ion-exchange chromatography and hydrophobic interaction chromatography in the preparation and recovery of polyethylene glycol-linked proteins

Cation- and anion-exchange chromatography can be used to purify a polyethylene glycol-linked protein dimer (PEG dimer) made with M, 20 000 PEG bis-vinylsulfone, even when there are no net charge differences between the components that are being separated. The retention time on ion-exchange generally is inversely proportional to the PEG:protein ratio (on a mass basis). One of the biggest challenges in developing the process for making this PEG dimer was the quality of the PEG linker. Reversed-phase HPLC can be used to determine both size heterogeneity and the degree of end-group activation of Mr 20 000 PEG bis-vinylsulfone. In addition, we have found that hydrophobic interaction chromatography can be used make more size homogeneous preparations of Mr 20000 PEG bis-vinylsulfone, which significantly increased the recovery of the PEG dimer.     

Prediction of the partitioning behaviour of proteins in aqueous two-phase systems using only their amino acid composition

The prediction of the partition behaviour of proteins in aqueous two-phase systems (ATPS) using mathematical models based on their amino acid composition was investigated. The predictive models are based on the average surface hydrophobicity (ASH). The ASH was estimated by means of models that use the three-dimensional structure of proteins and by models that use only the amino acid composition of proteins. These models were evaluated for a set of 11 proteins with known experimental partition coefficient in four-phase systems: polyethylene glycol (PEG) 4000/phosphate, sulfate, citrate and dextran and considering three levels of NaCl concentration (0.0% w/w, 0.6% w/w and 8.8% w/w). The results indicate that such prediction is feasible even though the quality of the prediction depends strongly on the ATPS and its operational conditions such as the NaCl concentration. The ATPS 0 model which use the three-dimensional structure obtains similar results to those given by previous models based on variables measured in the laboratory. In addition it maintains the main characteristics of the hydrophobic resolution and intrinsic hydrophobicity reported before. Three mathematical models, ATPS I-III, based only on the amino acid composition were evaluated. The best results were obtained by the ATPS I model which assumes that all of the amino acids are completely exposed. The performance of the ATPS I model follows the behaviour reported previously, i.e. its correlation coefficients improve as the NaCl concentration increases in the system and, therefore, the effect of the protein hydrophobicity prevails over other effects such as charge or size. Its best predictive performance was obtained for the PEG/dextran system at high NaCl concentration. An increase in the predictive capacity of at least 54.4% with respect to the models which use the three-dimensional structure of the protein was obtained for that system. In addition, the ATPS I model exhibits high correlation coefficients in that system being higher than 0.88 on average. The ATPS I model exhibited correlation coefficients higher than 0.67 for the rest of the ATPS at high NaCl concentration. Finally, we tested our best model, the ATPS I model, on the prediction of the partition coefficient of the protein invertase. We found that the predictive capacities of the ATPS I model are better in PEG/dextran systems, where the relative error of the prediction with respect to the experimental value is 15.6%.     

Application of a pH responsive multimodal hydrophobic interaction chromatography medium for the analysis of glycosylated proteins

Protein glycosylation has significant effects on the structure and function of proteins. The efficient separation and enrichment of glycoproteins from complex biological samples is one key aspect and represents a major bottleneck of glycoproteome research. In this paper, we have explored pH multimodal hydrophobic interaction chromatography to separate glycosylated from non-glycosylated forms of proteins. Three different proteins, ribonuclease, invertase and IgG, have been examined and different glycoforms have been identified. The media itself shows strong responsiveness to small variations in pH, which makes it possible to fine-tune the chromatographic conditions according to the properties of the protein isolated. Optimal glycoprotein separation has been obtained at pH 4. The pH responsive multimodal HIC medium in contrast to conventional HIC media is able to resolve contaminating DNA.     

High-performance liquid chromatography of amino acids, peptides and proteins. CXV. Thermodynamic behaviour of peptides in reversed-phase chromatography.

The thermodynamic behaviour of three peptides, bombesin, beta-endorphin and glucagon, was studied under reversed-phase high-performance liquid chromatographic conditions. Experimental data related to the interactive surface contact area (S values) and solute affinity (log k0) were derived over a range of temperatures between 5 and 85 degrees C. These experimental conditions allowed changes in the secondary structure of the solute to be monitored. The influence of the nature of the stationary phase ligand on the relative conformational stability of the three peptides was analysed by acquiring data with n-octadecyl silica (C18) and n-butyl silica (C4) sorbents. Values for the relative changes in entropy and enthalpy associated with the interactive process were also determined. The results provide further insight into the factors involved with the stabilization of secondary structure and the mechanism of the interaction of peptides with hydrophobic surfaces.     

Purification of metallothionein proteins from the crab, Portunus pelagicus-selectivity of hydrophobic interaction chromatography

Both hydrophobic interaction chromatography (HIC) and ion-exchange chromatography (IEC) are commonly used for the purification of proteins as their mild elution conditions preclude the dangers of denaturation. Of the two methods, IEC which fractionates proteins on the basis of charge differences, exhibits a higher selectivity. In the case of the metallothionein (MT) proteins, the charge states of the proteins are very susceptible to change due to oxidation of the abundant thiol groups. This complicates fractionation with IEC. Separation using HIC, on the other hand, is based on the exposed hydrophobic groups, which remain relatively intact if the tertiary structure of the proteins is not disrupted. In this work, the successful isolation of two MT isoforms from the tropical crab species, Portunus pelagicus, using HIC serves to demonstrate the high selectivity of this technique, in addition to the fact that it is indifferent to the state of oxidation of the MT proteins during the purification procedure.     

Separation of proteins from human plasma by sample displacement chromatography in hydrophobic interaction mode

Sample displacement chromatography (SDC) in reversed-phase and ion-exchange modes was introduced approximately 20 years ago. This method was first used for the preparative purification of peptides and proteins. Recently, SDC in ion-exchange mode was also successfully used for enrichment of low-abundance proteins from human plasma. In this paper, the use of SDC for the separation of plasma proteins in hydrophobic interaction mode is demonstrated. By use of two or more columns coupled in series during sample application, and subsequent elution of detached columns in parallel, additional separation of bound proteins was achieved. Further low-abundance, physiologically active proteins could be highly enriched and detected by ESI-MS/MS.     

Methodology for predicting the separation of proteins by hydrophobic interaction chromatography and its application to a cell extract

Ligand-exchange micellar electrokinetic capillary chromatography was used for the chiral resolution of underivatized and dansyl amino acid enantiomers simultaneously. The separation was achieved by chiral copper(II)-L-valine complexes incorporated in micelles of sodium dodecyl sulfate (SDS). The enantioresolution was strongly affected by SDS and a concentration of 20 mM SDS was shown to be necessary for the separation. Other impacting factors were investigated including pH, the molar ratio of copper(II) to L-valine and the total concentration of complex. Using the proposed method, 11 different dansyl amino acids and two underivatized amino acids were separated successfully with a running electrolyte of 20 mM NH4OAc, 4 mM CuSO4, 8 mM L-valine and 20 mM SDS at pH 9.0 in less than 25 min. Experiments were also performed with other amino acid ligands in order to vary the stability and the sterical arrangement of the copper(II) complexes and the possible chiral recognition mechanism was also discussed briefly.     

Modeling of adsorption in hydrophobic interaction chromatography systems using a preferential interaction quadratic isotherm

A preferential interaction quadratic isotherm model for hydrophobic interaction chromatographic systems is presented in this paper. In this isotherm, the nonlinear effect of salt on the capacity factor is described using the preferential interaction model developed by Perkins et al. [J. Chromatogr. A, 766 (1997) 1]. This is then coupled with a quadratic nonlinear isotherm to describe nonlinear adsorption behavior at high solute concentrations. The resulting preferential interaction quadratic isotherm is examined for its ability to describe solute adsorption behavior under both linear and nonlinear conditions over a wide range of salt concentrations in HIC systems. The results indicate that this isotherm is well suited for predicting nonlinear adsorption behavior in HIC systems for both proteins and low-molecular mass HIC displacers.