Protein adsorption isotherm behavior in hydrophobic interaction chromatography

From protein structure to an optimized chromatographic capture step using multiscale modeling

Optimizing a biopharmaceutical chromatographic purification process is currently the greatest challenge during process development. A lack of process understanding calls for extensive experimental efforts in pursuit of an optimal process. In silico techniques, such as mechanistic or data driven modeling, enhance the understanding, allowing more cost-effective and time efficient process optimization. This work presents a modeling strategy integrating quantitative structure property relationship (QSPR) models and chromatographic mechanistic models (MM) to optimize a cation exchange (CEX) capture step, limiting experiments. In QSPR, structural characteristics obtained from the protein structure are used to describe physicochemical behavior. This QSPR information can be applied in MM to predict the chromatogram and optimize the entire process. To validate this approach, retention profiles of six proteins were determined experimentally from mixtures, at different pH (3.5, 4.3, 5.0, and 7.0). Four proteins at different pH's were used to train QSPR models predicting the retention volumes and characteristic charge, subsequently the equilibrium constant was determined. For an unseen protein knowing only the protein structure, the retention peak difference between the modeled and experimental peaks was 0.2% relative to the gradient length (60 column volume). Next, the CEX capture step was optimized, demonstrating a consistent result in both the experimental and QSPR-based methods. The impact of model parameter confidence on the final optimization revealed two viable process conditions, one of which is similar to the optimization achieved using experimentally obtained parameters. The multiscale modeling approach reduces the required experimental effort by identification of initial process conditions, which can be optimized.

Laminin Adsorption and Adhesion of Neurons and Glial Cells on Carbon Implanted Titania Nanotube Scaffolds for Neural Implant Applications

Interfacing neurons persistently to conductive matter constitutes one of the key challenges when designing brain-machine interfaces such as neuroelectrodes or retinal implants. Novel materials approaches that prevent occurrence of loss of long-term adhesion, rejection reactions, and glial scarring are highly desirable. Ion doped titania nanotube scaffolds are a promising material to fulfill all these requirements while revealing sufficient electrical conductivity, and are scrutinized in the present study regarding their neuron-material interface. Adsorption of laminin, an essential extracellular matrix protein of the brain, is comprehensively analyzed. The implantation-dependent decline in laminin adsorption is revealed by employing surface characteristics such as nanotube diameter, ζ-potential, and surface free energy. Moreover, the viability of U87-MG glial cells and SH-SY5Y neurons after one and four days are investigated, as well as the material's cytotoxicity. The higher conductivity related to carbon implantation does not affect the viability of neurons, although it impedes glial cell proliferation. This gives rise to novel titania nanotube based implant materials with long-term stability, and could reduce undesirable glial scarring.

An Experimental and Modeling Combined Approach in Preparative Hydrophobic Interaction Chromatography

Chromatography is a technique widely used in the purification of biopharmaceuticals, and generally consists of several chromatographic steps. In this work, Hydrophobic Interaction Chromatography (HIC) is investigated as a polishing step for the purification of therapeutic proteins. Adsorption mechanisms in hydrophobic interaction chromatography are still not completely clear and a limited amount of published data is available. In addition to new data on adsorption isotherms for some proteins (obtained both by high-throughput and frontal analysis method), and a comparison of different models proposed in the literature, two different approaches are compared in this work to investigate HIC. The predictive approach exploits an in-house code that simulates the behavior of the component in the column using the model parameters found from the fitting of experimental data. The estimation approach, on the other hand, exploits commercial software in which the model parameters are found by the fitting of a few experimental chromatograms. The two approaches are validated on some bind-elute runs: the predictive approach is very informative, but the experimental effort needed is high; the estimation approach is more effective, but the knowledge gained is lower. The second approach is also applied to an in-development industrial purification process and successfully resulted in predicting the behavior of the system, allowing for optimization with a reduction in the time and amount of sample needed.

Use of Microfluidic Capillary Electrophoresis for the Determination of Multi-Component Protein Adsorption Isotherms: Application to High-Throughput Analysis for Hydrophobic Interaction Chromatography

Chromatography is a widely used separation process for purification of biopharmaceuticals that is able to obtain high purities and concentrations. The phenomena that occur during separation, mass transfer and adsorption are quite complex. To better understand these phenomena and their mechanisms, multi-component adsorption isotherms must be investigated. High-throughput methodologies are a very powerful tool to determine adsorption isotherms and they waste very small amounts of sample and chemicals, but the quantification of component concentrations is a real bottleneck in multi-component isotherm determination. The behavior of bovine serum albumin, Corynebacterium diphtheriae CRM₁₉₇ protein and lysozyme, selected as model proteins in binary mixtures with hydrophobic resin, is investigated here. In this work we propose a new method for determining multi-component adsorption isotherms using high-throughput experiments with filter plates, by exploiting microfluidic capillary electrophoresis. The precision and accuracy of the microfluidic capillary electrophoresis platform were evaluated in order to assess the procedure; they were both found to be high and the procedure is thus reliable in determining adsorption isotherms for binary mixtures. Multi-component adsorption isotherms were determined with a totally high-throughput procedure that turned out to be a very fast and powerful tool. The same procedure can be applied to every kind of high-throughput screening.

Protein-protein interactions and reduced excluded volume increase dynamic binding capacity of dual salt systems in hydrophobic interaction chromatography

Deploying two salts in hydrophobic interaction chromatography can significantly increase dynamic binding capacities. Nevertheless, the mechanistic understanding of this phenomenon is lacking. Here, we investigate whether surface tension or ionic strength govern dynamic binding capacities of the chromatographic resin Toyopearl Butyl-650 M in dual salt systems. Small-angle X-ray scattering was employed to analyze the model proteins and the protein-resin adduct in the respective dual salt systems. The dual salt systems incorporate sodium citrate and a secondary sodium salt (acetate, sulfate, or phosphate). As model proteins, we used lysozyme, GFP, and a monoclonal antibody (adalimumab). Moreover, for the protein-resin adduct, we determined the model parameters of a self-avoiding random walk model fitted into the pair density distribution function of the SAXS data. Ionic strength is more predictive for dynamic binding capacities in HIC dual salt systems than surface tension. However, dynamic binding capacities still differ by up to 30 % between the investigated dual salt systems. The proteins exhibit extensive protein-protein interactions in the studied dual salt HIC buffers. We found a correlation of protein-protein interactions with the well-known Hofmeister series. For systems with elevated protein-protein interactions, adsorption isotherms deviate from Langmuirian behavior. This highlights the importance of lateral protein-protein interactions in protein adsorption, where monomolecular protein layers are usually assumed. SAXS analysis of the protein-resin adduct indicates an inverse correlation of the binding capacity and the excluded volume parameter. This is indicative of the deposition of proteins in the cavities of the stationary phase. We hypothesize that increasing protein-protein interactions allow the formation of attractive clusters and multilayers in the cavities, respectively.

Water on hydrophobic surfaces: mechanistic modeling of polyethylene glycol-induced protein precipitation

For the purification of biopharmaceutical proteins, liquid chromatography is still the gold standard. Especially with increasing product titers, drawbacks like slow volumetric throughput and high resin costs lead to an intensifying need for alternative technologies. Selective preparative protein precipitation is one promising alternative technique. Although the capability has been proven, there has been no precipitation process realized for large-scale monoclonal antibody (mAb) production yet. One reason might be that the mechanism behind protein phase behavior is not completely understood and the precipitation process development is still empirical. Mechanistic modeling can be a means for faster, material-saving process development and a better process understanding at the same time. In preparative chromatography, mechanistic modeling was successfully shown for a variety of applications. Lately, a new isotherm for hydrophobic interaction chromatography (HIC) under consideration of water molecules as participants was proposed, enabling an accurate description of HIC. In this work, based on similarities between protein precipitation and HIC, a new precipitation model was derived. In the proposed model, the formation of protein-protein interfaces is thought to be driven by hydrophobic effects, involving a reorganization of the well-ordered water structure on the hydrophobic surfaces of the protein-protein complex. To demonstrate model capability, high-throughput precipitation experiments with pure or prior to the experiments purified proteins lysozyme, myoglobin, bovine serum albumin, and one mAb were conducted at various pH values. Polyethylene glycol (PEG) 6000 was used as precipitant. The precipitant concentration as well as the initial protein concentration was varied systematically. For all investigated proteins, the initial protein concentrations were varied between 1.5 mg/mL and 12 mg/mL. The calibrated models were successfully validated with experimental data. This mechanistic description of protein precipitation process offers mathematical explanation of the precipitation behavior of proteins at PEG concentration, protein concentration, protein size, and pH.

Water on hydrophobic surfaces: Mechanistic modeling of hydrophobic interaction chromatography

Mechanistic models are successfully used for protein purification process development as shown for ion-exchange column chromatography (IEX). Modeling and simulation of hydrophobic interaction chromatography (HIC) in the column mode has been seldom reported. As a combination of these two techniques is often encountered in biopharmaceutical purification steps, accurate modeling of protein adsorption in HIC is a core issue for applying holistic model-based process development, especially in the light of the Quality by Design (QbD) approach. In this work, a new mechanistic isotherm model for HIC is derived by consideration of an equilibrium between well-ordered water molecules and bulk-like ordered water molecules on the hydrophobic surfaces of protein and ligand. The model's capability of describing column chromatography experiments is demonstrated with glucose oxidase, bovine serum albumin (BSA), and lysozyme on Capto™ Phenyl (high sub) as model system. After model calibration from chromatograms of bind-and-elute experiments, results were validated with batch isotherms and prediction of further gradient elution chromatograms.

The influence of mixed salts on the capacity of HIC adsorbers: A predictive correlation to the surface tension and the aggregation temperature

Hydrophobic interaction chromatography (HIC) is one of the most frequently used purification methods in downstream processing of biopharmaceuticals. During HIC, salts are the governing additives contributing to binding strength, binding capacity, and protein solubility in the liquid phase. A relatively recent approach to increase the dynamic binding capacity (DBC) of HIC adsorbers is the use of salt mixtures. By mixing chaotropic with kosmotropic salts, the DBC can strongly be influenced. For salt mixtures with a higher proportion of chaotropic than kosmotropic salt, higher DBCs were achieved compared with single salt approaches. By measuring the surface tensions of the protein salt solutions, the cavity theory-proposed by Melander and Horváth-that higher surface tensions lead to higher DBCs, was found to be invalid for salt mixtures. Aggregation temperatures of lysozyme in the salt mixtures, as a degree of hydrophobic forces, were correlated to the DBCs. Measuring the aggregation temperatures has proven to be a fast analytical methodology to estimate the hydrophobic interactions and thus can be used as a measure for an increase or decrease in the DBCs. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:346-354, 2016.

Prediction of protein retention times in hydrophobic interaction chromatography by robust statistical characterization of their atomic-level surface properties

The correlation between the dimensionless retention times (DRT) of proteins in hydrophobic interaction chromatography (HIC) and their surface properties were investigated. A ternary atomic-level hydrophobicity scale was used to calculate the distribution of local average hydrophobicity across the proteins surfaces. These distributions were characterized by robust descriptive statistics to reduce their sensitivity to small changes in the three-dimensional structure. The applicability of these statistics for the prediction of protein retention behaviour was looked into. A linear combination of robust statistics describing the central tendency, heterogeneity and frequency of highly hydrophobic clusters was found to have a good predictive capability (R2  = 0.78), when combined a factor to account for protein size differences. The achieved error of prediction was 35% lower than for a similar model based on a description of the protein surface on an amino acid level. This indicates that a robust and mathematically simple model based on an atomic description of the protein surface can be used for the prediction of the retention behaviour of conformationally stable globular proteins with a well determined 3D structure in HIC. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:372-381, 2016.

Comparison of polymer induced and solvent induced trypsin denaturation: the role of hydrophobicity

Trypsin adsorption from aqueous buffer by various copolymers of allyl glycidyl ether-ethylene glycol dimethacrylate (AGE-EGDM) copolymer with varying crosslink density increases with increasing crosslink density and the effect slowly wears off after reaching a plateau at 50% crosslink density. The copolymer with 25% crosslink density was reacted with different amines with alkyl/aryl side chains to obtain a series of copolymers with 1,2-amino alcohol functional groups and varying hydrophobicity. Trypsin binding capacity again increases with hydrophobicity of the reacting amine and a good correlation between logPoctanol of the amine and protein binding is observed. The bound trypsin is denatured to the extent of 90% in spite of the presence of hydrophilic hydroxyl and amino groups. The behavior was comparable to that in mixtures of aqueous buffer and water-miscible organic co-solvents where the solvent concentration required to deactivate 50% of the enzyme (C50) is dependent on logPoctanol of the co-solvent.

Suitability of commercial hydrophobic interaction sorbents for temperature-controlled protein liquid chromatography under low salt conditions

The effect of temperature in the range from 10°C to 40°C and comparatively low ammonium sulfate (AS) concentrations of up to 0.5M on the adsorption of bovine serum albumin (BSA) on four different commercially available sepharose-based stationary phases was investigated. The determined isotherms were fitted by the Langmuir equation, and thermodynamic values were calculated by van't Hoff analysis. The adsorption of BSA onto the chromatographic resin Butyl Sepharose 4FF showed the strongest temperature influence; however, protein unfolding effects occurred when characterizing this system by dynamic column experiments, with an unfolded BSA fraction strongly attached to the sorbent. The percentage of the unfolding fraction was determined for different operating conditions and found to increase with the concentration of the cosmotropic salt, but even stronger with increasing temperature. Temperature-induced cyclic adsorption and desorption experiments were carried out to investigate the long-term performance of Butyl Sepharose 4FF by applying purely temperature-controlled regeneration. Over a period of five cycles, the working capacity remained stable, but BSA also started to accumulate on the column due to incomplete regeneration. Finally, the possibility to fractionate different proteins with a single temperature shift was shown by the complete separation of lysozyme and BSA. The results presented indicate that temperature-induced binding and elution may offer a possibility to shift the operation conditions of HIC resins toward reduced salt concentrations, thus saving chemicals and facilitating salt removal in further downstream processing stages.

Microcalorimetric study of the adsorption of PEGylated lysozyme and PEG on a mildly hydrophobic resin: influence of ammonium sulfate

Adsorption of native as well as mono-, di-, and tri-PEGylated lysozyme on Toyopearl PPG-600M, a mildly hydrophobic resin is studied by isothermal titration calorimetry and by independent adsorption equilibrium measurements in sodium phosphate buffer at pH 7.0 and 25 °C. For PEGylation two different PEG sizes are used (5 and 10 kDa) which leads to six different forms of PEGylated lysozyme all of which are systematically studied. Additionally, the adsorption of five pure PEGs is explored. The ammonium sulfate concentration is varied from 600 to 1200 mM. The molar enthalpy of adsorption Δh(p)(ads) is determined from the calorimetric and the adsorption equilibrium data. It is found to be endothermic in all experiments. The comparison of the adsorption of different PEGylated forms shows that the adsorption of PEGylated lysozyme is driven by the adsorption of the PEG chain. The results provide insight into the adsorption mechanisms of polymer-modified proteins on hydrophobic chromatographic resins.

A review of protein adsorption on bioceramics

Bioceramics, because of its excellent biocompatible and mechanical properties, has always been considered as the most promising materials for hard tissue repair. It is well know that an appropriate cellular response to bioceramics surfaces is essential for tissue regeneration and integration. As the in vivo implants, the implanted bioceramics are immediately coated with proteins from blood and body fluids, and it is through this coated layer that cells sense and respond to foreign implants. Hence, the adsorption of proteins is critical within the sequence of biological activities. However, the biological mechanisms of the interactions of bioceramics and proteins are still not well understood. In this review, we will recapitulate the recent studies on the bioceramic-protein interactions.

Protein stability and structure in HIC: hydrogen exchange experiments and COREX calculations

Hydrogen exchange mass spectrometry (HXMS) coupled to proteolytic digestion has been used to probe the conformation of bovine β-lactoglobulin (BLG), bovine α-lactalbumin (BLA), and human serum albumin (HSA) in solution and while adsorbed to the hydrophobic interaction chromatography media Phenyl Sepharose 6FF. All three proteins show evidence of EX1 exchange kinetics, indicating a loss of stability on the surface. HX protection patterns for all three proteins also indicate that the unfolded form is only partially solvent exposed. The hydrogen-deuterium exchange patterns of BLG and BLA on the surface suggest a structure that resembles each protein's respective solution phase molten globule state. The low stability of Domain II of HSA observed on Phenyl Sepharose 6FF also suggests a link to solution stability because Domain II is frequently cited as the least stable domain in solution unfolding pathways. COREX, an algorithm used to compute protein folding stabilities, correctly predicts solution hydrogen-deuterium exchange patterns for BLG and offers insight into its adsorbed phase stabilities but is unreliable for BLA predictions. The results of this work demonstrate a link between solution-phase local stability patterns and the nature of partially unfolded states that proteins can adopt on HIC surfaces.

Changes in solvent exposure reveal the kinetics and equilibria of adsorbed protein unfolding in hydrophobic interaction chromatography

Hydrogen exchange has been a useful technique for studying the conformational state of proteins, both in bulk solution and at interfaces, for several decades. Here, we propose a physically based model of simultaneous protein adsorption, unfolding and hydrogen exchange in HIC. An accompanying experimental protocol, utilizing mass spectrometry to quantify deuterium labeling, enables the determination of both the equilibrium partitioning between conformational states and pseudo-first order rate constants for folding and unfolding of adsorbed protein. Unlike chromatographic techniques, which rely on the interpretation of bulk phase behavior, this methodology utilizes the measurement of a molecular property (solvent exposure) and provides insight into the nature of the unfolded conformation in the adsorbed phase. Three model proteins of varying conformational stability, alpha-chymotrypsinogen A, beta-lactoglobulin B, and holo alpha-lactalbumin, are studied on Sepharose HIC resins possessing assorted ligand chemistries and densities. alpha-Chymotrypsinogen, conformationally the most stable protein in the set, exhibits no change in solvent exposure at all the conditions studied, even when isocratic pulse-response chromatography suggests nearly irreversible adsorption. Apparent unfolding energies of adsorbed beta-lactoglobulin B and holo alpha-lactalbumin range from -4 to 3 kJ/mol and are dependent on resin properties and salt concentration. Characteristic pseudo-first order rate constants for surface-induced unfolding are 0.2-0.9 min(-1). While poor protein recovery in HIC is often associated with irreversible unfolding, this study documents that non-eluting behavior can occur when surface unfolding is reversible or does not occur at all. Further, this hydrogen exchange technique can be used to assess the conformation of adsorbed protein under conditions where the protein is non-eluting and chromatographic methods are not applicable.

Studies of lysozyme binding to histamine as a ligand for hydrophobic charge induction chromatography

Histamine was immobilized on Sepharose CL-6B (Sepharose) for use as a ligand of hydrophobic charge induction chromatography (HCIC) of proteins. Lysozyme adsorption onto Histamine-Sepharose (HA-S) was studied by adsorption equilibrium and calorimetry to uncover the thermodynamic mechanism of the protein binding. In both the experiments, the influence of salt (ammonium sulfate and sodium sulfate) was examined. Adsorption isotherms showed that HA-S exhibited a high salt tolerance in lysozyme adsorption. This property was well explained by the combined contributions of hydrophobic interaction and aromatic stacking. The isotherms were well fitted to the Langmuir equation, and the equilibrium parameters for lysozyme adsorption were obtained. In addition, thermodynamic parameters (DeltaH(ads), DeltaS(ads), and DeltaG(ads)) for the adsorption were obtained by isothermal titration calorimetry by titrating lysozyme solutions into the adsorbent suspension. Furthermore, free histamine was titrated into lysozyme solution in the same salt-buffers. Compared with the binding of lysozyme to free histamine, lysozyme adsorption onto HA-S was characterized by a less favorable DeltaG(ads) and an unfavorable DeltaS(ads) because histamine was covalently attached to Sepharose via a three-carbon-chain spacer. Consequently, the immobilized histamine could only associate with the residues on the protein surface rather than those in the hydrophobic pocket, causing a less favorable orientation between histamine and lysozyme. Further comparison of thermodynamic parameters indicated that the unfavorable DeltaS(ads) was offset by a favorable DeltaH(ads), thus exhibiting typical enthalpy-entropy compensation. Moreover, thermodynamic analyses indicated the importance of the dehydration of lysozyme molecule and HA-S during the adsorption and a substantial conformational change of the protein during adsorption. The results have provided clear insights into the adsorption mechanisms of lysozyme onto the new HCIC material.

Adsorption induced enzyme denaturation: the role of polymer hydrophobicity in adsorption and denaturation of alpha-chymotrypsin on allyl glycidyl ether (AGE)-ethylene glycol dimethacrylate (EGDM) copolymers

Effects of changes in hydrophobicity of polymeric support on structure and activity of alpha-chymotrypsin (E.C. 3.4.21.1) have been studied with copolymers of allyl glycidyl ether (AGE) and ethylene glycol dimethacrylate (EGDM) with increasing molar ratio of EGDM to AGE (cross-link density 0.05 to 1.5). The enzyme is readily adsorbed from aqueous buffer at room temperature following Langmuir adsorption isotherms in unexpectedly large amounts (25% w/w). Relative hydrophobicity of the copolymers has been assessed by studying adsorption of naphthalene and Fmoc-methionine by the series of copolymers from aqueous solutions. Polymer hydrophobicity appears to increase linearly on increasing cross-link density from 0.05 to 0.25. Further increase in cross-link density causes a decrease in naphthalene binding but has little effect on binding of Fmoc-Met. Binding of alpha-chymotrypsin to these copolymers follow the trend for Fmoc-methionine binding, rather than naphthalene binding, indicating involvement of polar interactions along with hydrophobic interactions during binding of protein to the polymer. The adsorbed enzyme undergoes extensive denaturation (ca. 80%) with loss of both tertiary and secondary structure on contact with the copolymers as revealed by fluorescence, CD and Raman spectra of the adsorbed protein. Comparison of enzyme adsorption behavior with Eupergit C, macroporous Amberlite XAD-2, and XAD-7 suggests that polar interactions of the EGDM ester functional groups with the protein play a significant role in enzyme denaturation.

A new thermodynamic model describes the effects of ligand density and type, salt concentration and protein species in hydrophobic interaction chromatography

A new thermodynamic model is derived that describes both loading and pulse-response behavior of proteins in hydrophobic interaction chromatography (HIC). The model describes adsorption in terms of protein and solvent activities, and water displacement from hydrophobic interfaces, and distinguishes contributions from ligand density, ligand type and protein species. Experimental isocratic response and loading data for a set of globular proteins on Sepharose resins of various ligand types and densities are described by the model with a limited number of parameters. The model is explicit in ligand density and may provide insight into the sensitivity of protein retention to ligand density in HIC as well as the limited reproducibility of HIC data.