Modeling of overloaded gradient elution of nociceptin/orphanin FQ in reversed-phase liquid chromatography

Impact of changes in physicochemical parameters of the mobile phase along the column on the retention time in gradient liquid chromatography. Part A - temperature gradient

Gradient chromatography has been widely applied in analytical and preparative chromatography since it provides better peak shapes and the ability to elute analytes in a shorter time frame. Apart from changes in the composition of a mobile phase also alteration of process temperature can be applied during separation procedures to improve efficiency. However, proper mathematical modeling of the gradient chromatography and further correct prediction of solutes' retention behavior have become a serious challenge as it involves the need to develop computational procedures that accurately account for the time and spatial gradients of crucial parameters. In this work, a computational procedure including the equilibrium-dispersive two-dimensional mass transfer model, the two-dimensional (2D) heat transfer model together with Darcy's law and the continuity equation have been proposed. Additionally, the calculation procedure was simplified by replacing the 2D model with the one-dimensional (1D) mass transfer model in order to speed up the computations. Both proposed solutions were validated employing external experimental data of temperature gradient HPLC [1] as well as with predictions based on the linear elution strength (LES) model available therein. The proposed procedures made it possible to efficiently predict the concentration profiles with average relative errors of calculated retention times not exceeding 3.22%. Moreover, the effect of the axial dispersion coefficient determination method on the obtained peak shapes was examined involving the Gunn, the Wen-Fan, and the Chung-Wen correlations, indicating that the latter produces the most accurate results. Finally, the proposed mathematical procedures were tested under UHPLC conditions, and due to a significant difference in retention times found the 2D model is strongly advised.

Downstream Processing of Therapeutic Peptides by Means of Preparative Liquid Chromatography

The market of biomolecules with therapeutic scopes, including peptides, is continuously expanding. The interest towards this class of pharmaceuticals is stimulated by the broad range of bioactivities that peptides can trigger in the human body. The main production methods to obtain peptides are enzymatic hydrolysis, microbial fermentation, recombinant approach and, especially, chemical synthesis. None of these methods, however, produce exclusively the target product. Other species represent impurities that, for safety and pharmaceutical quality reasons, must be removed. The remarkable production volumes of peptide mixtures have generated a strong interest towards the purification procedures, particularly due to their relevant impact on the manufacturing costs. The purification method of choice is mainly preparative liquid chromatography, because of its flexibility, which allows one to choose case-by-case the experimental conditions that most suitably fit that particular purification problem. Different modes of chromatography that can cover almost every separation case are reviewed in this article. Additionally, an outlook to a very recent continuous chromatographic process (namely Multicolumn Countercurrent Solvent Gradient Purification, MCSGP) and future perspectives regarding purification strategies will be considered at the end of this review.

Modern trends in downstream processing of biotherapeutics through continuous chromatography: The potential of Multicolumn Countercurrent Solvent Gradient Purification

Single-column (batch) preparative chromatography is the technique of choice for purification of biotherapeutics but it is often characterized by an intrinsic limitation in terms of yield-purity trade-off, especially for separations containing a larger number of product-related impurities. This drawback can be alleviated by employing multicolumn continuous chromatography. Among the different methods working in continuous mode, in this paper we will focus in particular on Multicolumn Countercurrent Solvent Gradient Purification (MCSGP) which has been specifically designed for challenging separations of target biomolecules from their product-related impurities. The improvements come from the automatic internal recycling of the impure fractions inside the chromatographic system, which results in an increased yield without compromising the purity of the pool. In this article, steps of the manufacturing process of biopharmaceuticals will be described, as well as the advantages of continuous chromatography over batch processes, by particularly focusing on MCSGP.

From batch to continuous chromatographic purification of a therapeutic peptide through multicolumn countercurrent solvent gradient purification

A twin-column Multicolumn Countercurrent Solvent Gradient Purification (MCSGP) process has been developed for the purification of a therapeutic peptide, glucagon, from a crude synthetic mixture. This semi-continuous process uses two identical columns operating either in interconnected or in batch mode, thus enabling the internal recycle of the portions of the eluting stream which do not comply with purity specifications. Because of this feature, which actually results in the simulated countercurrent movement of the stationary phase with respect to the mobile one, the yield-purity trade-off typical of traditional batch preparative chromatography can be alleviated. Moreover, the purification process can be completely automatized. Aim of this work is to present a simple procedure for the development of the MCSGP process based on a single batch experiment, in the case of a therapeutic peptide of industrial relevance. This allowed to recover roughly 90% of the injected glucagon in a purified pool with a purity of about 90%. A comparison between the performance of the MCSGP process and the classical single column batch process indicates that percentage increase in the recovery of target product is +23% when transferring the method from batch conditions to MCSGP, with an unchanged purity of around 89%. This improvement comes at the expenses of a reduction of about 38% in productivity.

Modeling the nonlinear behavior of a bioactive peptide in reversed-phase gradient elution chromatography

The thermodynamic behavior of octreotide, a cyclic octapeptide with important pharmaceutical functions, has been simulated under reversed-phase gradient elution conditions. To this end, adsorption behavior was firstly investigated in isocratic conditions, under a variety of water/acetonitrile + 0.02% (v/v) trifluoroacetic acid (TFA) mixtures as mobile phase by using a Langmuir isotherm. Organic modifier was varied in the range between 23 and 28% (v/v). Adsorption isotherms were determined by means of the so-called Inverse Method (IM) with a minimum amount of peptide. The linear solvent strength (LSS) model was used to find the correlation between isotherm parameters and mobile phase composition. This study contributes to enlarge our knowledge on the chromatographic behavior under nonlinear gradient conditions of peptides. In particular, it focuses on a cyclic octapeptide.

Estimation of Nonlinear Adsorption Isotherms in Gradient Elution RP-LC of Peptides in the Presence of an Adsorbing Additive

ABSTRACT

In electrostatic repulsive interaction chromatography, using a charged surface hybrid sorbent carrying positive charges can improve the peak shape of peptides in reversed-phase liquid chromatography (RP-LC), especially in overloaded conditions, compared with standard C18 sorbents. However, the positive surface charges can interact with anionic additives commonly used in peptide separations, e.g., trifluoroacetic acid (TFA), complicating adsorption isotherm estimation. We investigated how the competition for available adsorption sites between TFA and two peptides influenced the adsorption isotherm in gradient elution. A model accounting for the competition with TFA was compared with a model neglecting TFA adsorption. We found that the two models predicted elution profiles with the same accuracy. We also found that the adsorption isotherms were extremely similar in shape, leading to the conclusion that neglecting the competition with TFA is a valid approximation enabling faster and more robust adsorption isotherm estimation for the studied type of sorbent.

GRAPHICAL ABSTRACT

The importance of ion-pairing in peptide purification by reversed-phase liquid chromatography

The adsorption mechanism for three peptides was studied under overloaded conditions through adsorption isotherm measurements in the presence of an ion-pairing reagent, trifluoroacetic acid (TFA), on an end-capped C₁₈-bonded stationary phase. The overall aim of the study was to obtain a better understanding of how the acetonitrile and the TFA fractions in the eluent affected the overloaded elution profiles and the selectivity between peptides using mechanistic modelling and multivariate design of experiments. When studying the effect of TFA, direct evidence for ion pair formation between a peptide and TFA in acetonitrile-water solutions was provided by fluorine-proton nuclear Overhauser NMR enhancement experiments and the adsorption of TFA on the stationary phase was measured by frontal analysis. The adsorption isotherms for each peptide were then determined by the inverse method at eight TFA concentrations ranging from 2.6mM to 37.3mM (0.02-0.29vol-%) in isocratic elution. The equilibrium between the peptide ion and the peptide-TFA complex was modelled by coupling the mass-balance to reaction kinetics and determining separate adsorption isotherms for the two species. We found that a Langmuir isotherm described the elution profile of peptide-TFA complex well while the peptide ion was described by a bi-Langmuir adsorption isotherm since it exhibited strong secondary interactions. The elution profiles had an unfavorable shape at low TFA concentrations consisting of a spike in their front and a long tailing rear due to the secondary interactions for the peptide ion having very low saturation capacity. The acetonitrile dependence on the adsorption isotherms was studied by determination of adsorption isotherms directly from elution profiles obtained in gradient elution which enabled a broad acetonitrile interval to be studied. Here, it was found that the column saturation capacity was quickly reached at very low acetonitrile fractions and that there were significant variations in adsorption with the molecular weight. Finally, practical implications for method development are discussed based on an experimental design where gradient slope and TFA concentrations are used as factors.

Choice of Model for Estimation of Adsorption Isotherm Parameters in Gradient Elution Preparative Liquid Chromatography

The inverse method is a numerical method for fast estimation of adsorption isotherm parameters directly from a few overloaded elution profiles and it was recently extended to adsorption isotherm acquisition in gradient elution conditions. However, the inverse method in gradient elution is cumbersome due to the complex adsorption isotherm models found in gradient elution. In this case, physicochemically correct adsorption models have very long calculation times. The aim of this study is to investigate the possibility of using a less complex adsorption isotherm model, with fewer adjustable parameters, but with preserved/acceptable predictive abilities. We found that equal or better agreement between experimental and predicted elution profiles could be achieved with less complex models. By being able to select a model with fewer adjustable parameters, the calculation times can be reduced by at least a factor of 10.

Estimation of adsorption isotherm parameters with inverse method—Possible problems

In recent years the inverse method (IM) has been frequently applied to estimate of isotherm parameters. The IM has been used for adsorption process modeling for one, two and even three component chromatography. This method requires only a few injections with various sample concentrations, so the solute consumption and time requirements are very modest. The successful estimation of isotherm parameters with IM depends on applied chromatography column model and a numerical method used to solve the model. For HPLC column the classical equilibrium-dispersive (ED) model can be used. This model is solved frequently with very fast Rouchon finite difference method. However, the accuracy of computations with Rouchon method is decreasing with increase of the number of analyzed components. The aim of this work is the comparison of the results obtained with inverse method when ED model was solved with Rouchon or orthogonal collocation on finite element (OCFE) scheme. Assuming that solution of ED model with OCFE method can be regarded as real a solution, it was found that the Rouchon scheme may not give satisfactory results even for column with 10,000 theoretical plates for three component chromatography. Moreover, the optimal conditions for separation, calculated with Rouchon method, can be remarkably different from that obtained with the OCFE method. The next aim of this work is the presentation of Craig method application to estimation of model parameters with IM.

Modeling and optimization of preparative reversed-phase liquid chromatography for insulin purification

This paper presents a model for reversed-phase purification of insulin from desamido insulin. The system is described by a reaction dispersive model with a competitive Langmuir isotherm. A model building and calibration method is presented and the model's region of validity is defined. The model is calibrated using only two-component experiments on the raw mixture by the inverse method and then experimentally validated. The model is then used to optimize the system's production rate with both purity and yield requirements. The yield requirement is varied between 80 and 95% to study the effect on the production rate and the operating point. The operating points found with the optimization were found outside the model's region of validity, but the experimental validation of the operating points shows that the model can be extrapolated to the interesting operating points.

Analysis of the Ideal Model for a Single Component in Preparative Gradient Elution Chromatography

The analytical solutions of the ideal model of chromatography for a single component were studied by using the method of characteristics for solving partial differential equations. By assuming that the mobile-phase gradient entering the column is not affected by any sorption of mobile-phase solvents on the column, we demonstrated that the stationary-phase concentration of the component is kept constant along the characteristic curve irrespective of the isotherm and gradient profile used. By using this property and assuming a Langmiur isotherm, we also calculated the band profiles for single linear, stepwise, and ladderlike gradients, respectively. These band profiles were compared with those computed by using the finite-difference method. The results were found to be in good agreement.

Determination of competitive adsorption isotherm parameters of pindolol enantiomers on alpha1-acid glycoprotein chiral stationary phase

In this paper, inverse method (IM) was used to determine the binary competitive adsorption isotherm of pindolol enantiomers by a least-square fitting of the proposed model to the experimentally measured elution curves of racemic pindolol. The isotherm parameters were determined by minimizing the least-square error using an adaptation of genetic algorithm, non-dominated sorting genetic algorithm with jumping genes (NSGA-II-JG). An equilibrium dispersive (ED) model combined with bi-Langmuir isotherm was used in predicting the elution profiles. The determined parameters show good agreement with the experimental profiles at various experimental conditions such as sample volume, concentration and flow rates of the racemic mixture. Robustness and validity of the isotherm parameters were also verified by frontal analyses at various step inputs. Results from both the pulse tests and the frontal analysis indicate that adsorption isotherm derived from the inverse method is quite reliable. This method requires relatively less number of experiments to be performed and therefore, lower experimental costs confirming that inverse method is an attractive alternative approach of experimental technique in determining the competitive adsorption isotherm for binary systems.