Feedback regulation in preparative elution chromatography.

Input-output models utilizing discrete variables are developed for elution chromatography and used together with control theory to investigate behavior patterns between output variables (e.g., yield, purity, and production rate) and input variables (e.g., cut point locations and feed slug size). Variable interactions and controllability characteristics for isocratic and stepwise gradient elution in linear and nonlinear chromatography are investigated by using the relative gain array, singular value decomposition, and direct simulation. Inverse-based regulatory systems are evaluated for updating process inputs by using information on process outputs such that product quality and system performance parameters are efficiently maintained for the case where the column is overloaded and the solute peaks overlap with each other. Strategies for feedback optimizing control are also discussed.     

Optimizing the sample size and the retention parameters to achieve maximum production rates for enantiomers in chiral chromatography

The optimum experimental conditions (sample size and mobile phase composition) are calculated for maximum production rate of either one of two enantiomers contained in feeds of different compositions (1/1, 1/10, and 10/1). The products are obtained at 99% purity. The calculations use the equilibrium-dispersive model of chromatography and the equilibrium isotherms determined experimentally from the rear, diffuse boundary of overloaded elution profiles. The production rate measured experimentally under the optimum conditions calculated agree with 4% of the calculated values. There is an optimum value for the retention factor which is higher than predicted by a model assuming constant separation factor, because both separation factor and retention decrease with increasing organic solvent concentration in the mobile phase. (c) 1992 John Wiley & Sons, Inc.     

Purification of recombinant brain derived neurotrophic factor using reversed phase displacement chromatography

This work investigates the utility of RPLC displacement chromatography for the purification of recombinant brain derived neurotrophic factor (rHu-BDNF) from its variants and E. coli. protein (ECP) impurities. The closely associated variants (six in total) differ by one amino acid from the native BDNF and thus pose a challenging separation problem. Several operational parameters were investigated to study their effects on the yield of the displacement process. The results indicated that the concentration of trifluoroacetic acid (TFA) in the buffer was a key factor in achieving the desired purification. Displacement chromatography on an analytical scale column resulted in extremely high purity and yield in a single chromatographic step. The process was successfully scaled-up with respect to particle and column diameter. The production rate of a pilot scale RPLC displacement process was shown to be 23 times higher than the combined production rates of the current preparative ion exchange and hydrophobic interaction gradient elution steps that are used to remove variant and ECP impurities, respectively.     

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.     

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.     

Rapid isolation of monosialogangliosides from bovine brain gangliosides by selective-overload chromatography.

A procedure for rapid isolation of monosialogangliosides from purified bovine brain gangliosides has been developed. It utilizes the selective difference in association between monosialogangliosides and polysialogangliosides for the ion-exchange resin Q-Sepharose. When the ion-exchange column is overloaded with a bovine brain ganglioside mixture in the proper ganglioside to column bed-volume ratio, the polysialogangliosides are selectively retained by the column while the monosialogangliosides emerge with the void volume without the use of salt for elution. With the critical ganglioside to bed-volume ratio (1 g:8.32 ml), and an appropriate column bed-height to column radius ratio of 6.9, monosialogangliosides are reproducibly obtained in high purity with greater than 90% yield. The method has been used at both the analytical and preparative scale. We call this separation technique selective-overload chromatography.     

Comparison of batch and continuous multi‐column protein A capture processes by optimal design

Multi‐column capture processes show several advantages compared to batch capture. It is however not evident how many columns one should use exactly. To investigate this issue, twin‐column CaptureSMB, 3‐ and 4‐column periodic counter‐current chromatography (PCC) and single column batch capture are numerically optimized and compared in terms of process performance for capturing a monoclonal antibody using protein A chromatography. Optimization is carried out with respect to productivity and capacity utilization (amount of product loaded per cycle compared to the maximum amount possible), while keeping yield and purity constant. For a wide range of process parameters, all three multi‐column processes show similar maximum capacity utilization and performed significantly better than batch. When maximizing productivity, the CaptureSMB process shows optimal performance, except at high feed titers, where batch chromatography can reach higher productivity values than the multi‐column processes due to the complete decoupling of the loading and elution steps, albeit at a large cost in terms of capacity utilization. In terms of trade‐off, i.e. how much the capacity utilization decreases with increasing productivity, CaptureSMB is optimal for low and high feed titers, whereas the 3‐column process is optimal in an intermediate region. Using these findings, the most suitable process can be chosen for different production scenarios.     

Continuous separation of multicomponent protein mixtures by annular displacement chromatography

Displacement chromatography is a predominantly nonlinear mode of chromatography, which has certain advantages over the elution mode for preparative bioseparations. Whereas continuous production (and separation) processes have their theoretical benefits in this context, protein displacement chromatography has up to now only been performed in the batch mode. In this contribution, we demonstrate that the principle of continuous annular chromatography can be adapted to displacement chromatography. Separations of up to three standard proteins (two whey proteins, soybean trypsin inhibitor) were developed and optimized using a small (4 x 250 mm) batch column. These separations were subsequently transferred directly to the continuous system (500-mL column). Separations of similar quality in terms of final product purity and recovery yield were obtained using the continuous system.     

Predicting the elution behavior of proteins in affinity chromatography on non-porous particles.

Affinity chromatography on non-porous particles of microsize is particularly useful for the rapid analysis and micropreparative separation of proteins. The elution behavior of proteins in an affinity column packed with non-porous copolymerized particles of styrene, methyl methacrylate and glycidyl methacrylate was investigated both theoretically and experimentally, using the lysozyme-Cibacron Blue 3G-A affinity system. Equations used to predict the elution profiles, resulting from the elution by increasing the ionic strength (NaCl concentration) in the mobile phase, were obtained. The maximum adsorbate concentration, desorption rate constant and equilibrium constant under elution conditions were determined by matching experimental data with predicted elution profiles. Based on the parameters determined at a flow-rate of 0.5 ml/min and with 1 M NaCl in the elution buffer, the model equations could predict the elution profiles for other experimental runs, where different flow-rates and sodium chloride concentrations were used. Both the experimental and predicted results revealed that the affinity interaction kinetics are not significantly influenced by the flow-rate and, hence, the film mass transfer. To elute bound lysozyme from immobilized dye ligand, a higher value of the ionic strength leads to a faster elution and a sharper elution peak. The influence of elution conditions on the kinetic and thermodynamic parameters and, consequently, on the elution peak profiles was evaluated. The model equations can also predict the behavior of protein elution from an affinity column by changing the pH of the mobile phase, according to a previous study.     

Screening of protein-ligand interactions by affinity chromatography

This paper examines affinity chromatography (AC) as an alternative tool for the determination of protein-ligand interactions for the particular case in which the ligand is the same protein. The methodology is less labor-intensive and more sample-efficient than traditional methods used to measure the second virial coefficient (B(22)), a parameter commonly used to evaluate protein-protein interactions. The chromatographic capacity factor (k') was studied for lysozyme and equine serum albumin for a wide range of experimental solution conditions such as crystallizing agent concentration, protein concentration and pH. Parallel experiments using AC to determine k' and static light scattering (SLS) to determine B(22) showed that the two parameters were highly correlated. Two different column volumes ( approximately 1 and approximately 0.1 mL) were tested and gave essentially the same values for k', showing the feasibility of miniaturization.     

铜藻岩藻黄素提取及纯化工艺研究

为了对岩藻黄素的提取、纯化进行系统研究,进而为高纯度岩藻黄素的工业化生产提供研究基础,筛选了适用于提取铜藻（Sargassum horneri）鲜藻中岩藻黄素的有机溶剂,并通过单因素实验和正交实验确定了最佳的提取溶剂浓度、提取温度、提取时间、料液比等工艺参数。随后采用硅胶柱层析法进行纯化,并通过单因素实验确定了最佳的硅胶柱床高度、上样量和洗脱流速。最后采用制备液相法对经层析纯化的岩藻黄素进一步纯化。结果表明,有机溶剂萃取的最佳工艺条件为：甲醇浓度90%,提取温度50 ℃,提取时间1 h,料液比1∶10,此条件下岩藻黄素提取率达到(0.258 9±0.003 6) mg·g-1鲜重（FW）［(1.078 8±0.015 0) mg·g-1干重（DW）］。硅胶柱层析的最佳工艺条件为：硅胶柱床高度10 cm,上样量6 g,洗脱流速10 mL·min-1,此条件下岩藻黄素得率为0.176 5 mg·g-1FW（0.735 3 mg·g-1 DW）,纯度为87.01%±0.88%。经制备液相进一步纯化后,岩藻黄素得率为0.127 1 mg·g-1 FW（0.529 4 mg·g-1 DW）,纯度为99.27%±0.22%。研究所用工艺简单,岩藻黄素得率高,为高纯度岩藻黄素的制备提供了实验基础。     

True and apparent temperature dependence of protein adsorption equilibrium in reversed-phase HPLC

The adsorption behavior of bovine insulin on a C(8)-bonded silica stationary phase was investigated at different column pressures and temperatures in isocratic reversed-phase HPLC. Changes in the molar volume of insulin (deltaV(m)) upon adsorption were derived from the pressure dependence of the isothermal retention factor (k'). The values of deltaV(m) were found to be practically independent of the temperature between 25 and 50 degrees C at -96 mL/mol and to increase with increasing temperature, up to -108 mL/mol reached at 50 degrees C. This trend was confirmed by two separate series of measurements of the thermal dependence of ln(k'). In the first series the average column pressure was kept constant. The second series involved measurements of ln(k') under constant mobile-phase flow rate, the average column pressure varying with the temperature. In both cases, a parabolic shape relationship was observed between ln(k') and the temperature, but the values obtained for ln k' were higher in the first than in the second case. The relative difference in ln(k'), caused by the change in pressure drop induced by the temperature, is equivalent to a systematic error in the estimate of the Gibbs free energy of 12%. Thus, a substantial error is made in the estimates of the enthalpy and entropy of adsorption when neglecting the pressure effects associated with the change in the molar volume of insulin. This work proves that the average column pressure must be kept constant during thermodynamic measurements of protein adsorption constants, especially in RPLC and HIC. Our results show also that there is a critical temperature, T(c) approximately equals 53 degrees C, at which ln(k') is maximum and the insulin adsorption process changes from an exothermic to an endothermic one. This temperature determines also the transition point in the molecular mechanism of insulin adsorption that involves successive unfolding of the protein chain.     

A Stereoselective Synthesis of (±)-endo-Brevicomin,a Pheromone Inhibitor Produced by Dendroctonus Bark Beetles

Prediction of elution curves in gel chromatography was attempted on the basis of a mass balance model which gave consideration to gel phase diffusion and longitudinal dispersion in a column. The basic differential equations for the model were solved by means of Laplace transformation, and then the solution in Laplace domain was inverted into time domain numerically. The calculated elution curves were in good agreement with the experimental ones of NaCl and myoglobin with various Sephadex gel columns. This indicated the validity of the calculation method and the model employed in this study.

Furthermore, the elution curves were calculated tentatively for various combinations of the parameters appearing in the mass balance model. Then, the magnitude of peak asymmetry, the shift of peak position and the maximum peak height of the elution curves were correlated with various parameters and operational variables. These correlations might permit prediction of suitable operational conditions for gel chromatography, especially for molecular weight determination.     

Cation exchange chromatography performed in overloaded mode is effective in removing viruses during the manufacturing of monoclonal antibodies

Viral safety is a critical concern with regard to monoclonal antibody (mAb) products produced in mammalian cells such as Chinese hamster ovary cells. Manufacturers are required to ensure the safety of such products by validating the clearance of viruses in downstream purification steps. Cation exchange (CEX) chromatography is widely used in bind/elute mode as a polishing step in mAb purification. However, bind/elute modes require a large volume of expensive resin. To reduce the production cost, the use of CEX chromatography in overloaded mode has recently been investigated. The viral clearance ability in overloaded mode was evaluated using murine leukemia virus (MLV). Even under high-load conditions such as 2,000 g mAb/L resin, MLV was removed from mAb solutions. This viral clearance ability was not significantly affected by resin type or mAb type. The overloaded mode can also remove other types of viruses such as pseudorabies virus and reovirus Type 3 from mAb solutions. Based on these results, this cost-effective overloaded mode is comparable to the bind-elute mode in terms of viral removal.     

Chemical dynamic approach to synergetics and instability in biological systems

Based on Haken's theory, self-organization or synergetics is discussed using chemical dynamics to represent an autocatalytic reaction. In a simple case the changes in a self-organizing system are given by a set of two rate equations for a primary and a partial system. When these systems mutually form a feedback loop under the adiabatic condition, the rate equation of self-organization is described by a generalized Gibbs' free energy change delta U (delta x) followed by the reaction. The sign of the parameter k'3 (k0-kex; k0, kex: rate constants with or without an external stimulus) determines the instability of the coupled system in quasi-equilibrium (k'3 approximately greater than 0; k0 greater than kex). When the product exceeds the threshold (k'3 less than 0; k0 less than kex), the system transfers into a new state, or a phase transition appears. Considering the Boltzmann distribution, the transition parameter k'3 is evaluated by an average distribution of the states and the instability is discussed using the reaction velocities vqe and vqe in the quasi-equilibrium state. As an example of this model membrane excitation is discussed briefly.     

Reaction mechanism of the magnesium ion-dependent adenosine triphosphatase of frog muscle myosin and subfragment 1.

The Mg2+-dependent ATPase (adenosine 5'-triphosphatase) mechanism of myosin and subfragment 1 prepared from frog leg muscle was investigated by transient kinetic technique. The results show that in general terms the mechanism is similar to that of the rabbit skeletal-muscle myosin ATPase. During subfragment-1 ATPase activity at 0-5 degrees C pH 7.0 and I0.15, the predominant component of the steady-state intermediate is a subfragment-1-products complex (E.ADP.Pi). Binary subfragment-1-ATP (E.ATP) and subfragment-1-ADP (E.ADP) complexes are the other main components of the steady-state intermediate, the relative concentrations of the three components E.ATP, E.ADP.Pi and E.ADP being 5.5:92.5:2.0 respectively. The frog myosin ATPase mechanism is distinguished from that of the rabbit at 0-5 degrees C by the low steady-state concentrations of E.ATP and E.ADP relative to that of E.ADP.Pi and can be described by: E + ATP k' + 1 in equilibrium k' - 1 E.ATP k' + 2 in equilibrium k' - 2 E.ADP.Pi k' + 3 in equilibrium k' - 3 E.ADP + Pi k' + 4 in equilibrium k' - 4 E + ADP. In the above conditions successive forward rate constants have values: k' + 1, 1.1 X 10(5)M-1.S-1; k' + 2 greater than 5s-1; k' + 3, 0.011 s-1; k' + 4, 0.5 s-1; k'-1 is probably less than 0.006s-1. The observed second-order rate constants of the association of actin to subfragment 1 and of ATP-induced dissociation of the actin-subfragment-1 complex are 5.5 X 10(4) M-1.S-1 and 7.4 X 10(5) M-1.S-1 respectively at 2-5 degrees C and pH 7.0. The physiological implications of these results are discussed.     

两种水溶性抗菌活性物质的分离提取

对水溶性、不解离的极性物质分离时,一般采用吸附层析和凝胶层析等途径。实验通过硅胶柱层析、葡聚糖凝胶柱层析以及硅胶ＧＦ２５４制备型薄板层析,从发酵液样品中分离出两种有抗菌活性的纯物质。薄层层析的展开剂为二氯甲烷－四氢呋喃－甲醇－水（２５：３０：２）,分离出的组分中Ｒｆ＝０．７和Ｒｆ＝０．８两种物质有抗菌活性。硅胶柱层析洗脱过程为梯度洗脱,先用１５０ｍｌ上述展开剂洗脱,再用二氯甲烷－甲醇（２０：８０）     

Polymer displacement/shielding in protein chromatography

An overview of different applications of polymer interactions with ion-exchange and dye-affinity chromatographic matrices is presented here. The strength of interaction between the ligand and the polymer plays a crucial role in deciding the mode of chromatographic application. Charged, non-ionic and thermosensitive polymers such as poly(ethylene imine), poly(N-vinyl pyrrolidone) and poly(vinyl caprolactam) respectively, show different degrees of interaction with the dye molecules in dye ligand chromatography. Polymers, with their ability of multipoint and hence strong attachment to the chromatographic matrices, were used as efficient displacers in displacement chromatography. The polymer displacement resulted in better recoveries and sharper elution profiles than traditional salt elutions. The globule–coil transition of the thermosensitive reversible soluble–insoluble polymer, poly(vinyl caprolactam), can be exploited in dye-affinity columns for the temperature induced displacement of the bound protein. In another situation, prior to the column chromatography of crude protein extract, polymers formed complexes with the dye matrix and “shielded” the column. The polymer shielding decreased the nonspecific interactions without affecting the specific interactions of the target protein to the dye matrix.     

Zonal rate model for stacked membrane chromatography part II: characterizing ion-exchange membrane chromatography under protein retention conditions

The Zonal Rate Model (ZRM) has previously been shown to accurately account for contributions to elution band broadening, including external flow nonidealities and radial concentration gradients, in ion-exchange membrane (IEXM) chromatography systems operated under nonbinding conditions. Here, we extend the ZRM to analyze and model the behavior of retained proteins by introducing terms for intra-column mass transfer resistances and intrinsic binding kinetics. Breakthrough curve (BTC) data from a scaled-down anion-exchange membrane chromatography module using ovalbumin as a model protein were collected at flow rates ranging from 1.5 to 20 mL min(-1). Through its careful accounting of transport nonidealities within and external to the membrane stack, the ZRM is shown to provide a useful framework for characterizing putative protein binding mechanisms and models, for predicting BTCs and complex elution behavior, including the common observation that the dynamic binding capacity can increase with linear velocity in IEXM systems, and for simulating and scaling separations using IEXM chromatography. Global fitting of model parameters is used to evaluate the performance of the Langmuir, bi-Langmuir, steric mass action (SMA), and spreading-type protein binding models in either correlating or fundamentally describing BTC data. When combined with the ZRM, the bi-Langmuir, and SMA models match the chromatography data, but require physically unrealistic regressed model parameters to do so. In contrast, for this system a spreading-type model is shown to accurately predict column performance while also providing a realistic fundamental explanation for observed trends, including an observed increase in dynamic binding capacity with flow rate.