Retention of ionizable compounds in high-performance liquid chromatography. 14. Acid-base pK values in acetonitrile-water mobile phases

Using competitive frontal analysis, the binary adsorption isotherms of the enantiomers of 1-phenyl-l-propanol were measured on a microbore column packed with a chiral stationary phase based on cellulose tribenzoate. These measurements were carried out using only the racemic mixture. The experimental data were fitted to four different isotherm models: Langmuir, BiLangmuir, Langmuir-Freundlich and Tóth. The BiLangmuir and the Langmuir-Freundlich models accounted best for the competitive adsorption data. An excellent agreement between the experimental and the calculated overloaded band profiles for various samples of racemic mixture was obtained when the equilibrium dispersive model of chromatography was used together with the BiLangmuir competitive isotherm. The isotherm parameters measured under competitive conditions were used to calculate the overloaded band profiles of large samples of the pure S- and R-enantiomers, too. A satisfactory agreement between the experimental and calculated band profiles was observed when using in the computation the corresponding single component BiLangmuir isotherm derived from the binary isotherm previously determined. Thus oniy data derived from the racemic mixture are required for computer optimization of the preparative chromatography separation of the enantiomers.     

Adsorption behavior and prediction of the band profiles of the enantiomers of 3-chloro-1-phenyl-1-propanol. Influence of the mass transfer kinetics

The single-component and competitive adsorption isotherms of the enantiomers of 3-chloro-1-phenyl-1-propanol were measured by frontal analysis. The stationary phase was a cellulose tribenzoate coated on silica, the mobile phase an n-hexane-ethyl acetate (95:5) solution. The adsorption data measured fitted well to the Langmuir isotherm model. The band profiles of single components and of their mixtures were calculated using the equilibrium-dispersive model. These profiles were found to match quite satisfactorily the experimental band profiles. However, the agreement between calculated and experimental band profiles was significantly improved when a more complex model taking into account the mass transfer kinetics was used. The mass transfer rate coefficients, k(f), for both single components were determined by using the transport-dispersive model of chromatography. The coefficients obtained were used to predict the band profiles of mixtures of the two enantiomers to good agreement.     

Study of the adsorption behavior of the enantiomers of 1-phenyl-1-propanol on a cellulose-based chiral stationary phase

Using single-step frontal analysis, we measured single-component and competitive adsorption isotherm data for the two enantiomers of 1-phenyl-1-propanol (PP). These experimental data were fitted to several competitive bi-Langmuir models (with 8, 6, 5 and 4 parameters) and to the competitive Langmuir model. The latter model accounted well for the behavior of both PP enantiomers on Chiracel OB (cellulose tribenzoate coated on silica gel). The parameters obtained were used in numerical calculations to predict the band profiles of the two single components and of their mixtures under overloaded conditions. The equilibrium-dispersive model provides satisfactory results, with minor differences between the calculated and the experimental profiles. These differences became negligible when a more complex kinetic model was used, with a concentration-dependent rate coefficient.     

Application of the general rate model and the generalized Maxwell-Stefan equation to the study of the mass transfer kinetics of a pair of enantiomers

The general rate model of chromatography can be coupled with the generalized Maxwell-Stefan equation that describes the surface diffusion flux. The resulting model is useful to describe the behavior of two enantiomers during their separation on chiral phases, cases in which the mass transfer kinetics is known to be sluggish. A case in point is the modeling of the elution profiles of the racemic mixture of the two enantiomers of 1-phenyl-1-propanol on cellulose tribenzoate coated on silica, a popular chiral stationary phase. The competitive equilibrium isotherm behavior of the two enantiomers on the chiral stationary phase was described using the competitive Tóth isotherm model. An excellent agreement between the experimental and the calculated profiles was observed in the whole range of experimental conditions investigated, at low and high column loadings.     

Separation mechanism of nortriptyline and amytriptyline in RPLC

The single and the competitive equilibrium isotherms of nortriptyline and amytriptyline were acquired by frontal analysis (FA) on the C18- bonded discovery column, using a 28/72 (v/v) mixture of acetonitrile and water buffered with phosphate (20 mM, pH 2.70). The adsorption energy distributions (AED) of each compound were calculated from the raw adsorption data. Both the fitting of the adsorption data using multi-linear regression analysis and the AEDs are consistent with a trimodal isotherm model. The single-component isotherm data fit well to the tri-Langmuir isotherm model. The extension to a competitive two-component tri-Langmuir isotherm model based on the best parameters of the single-component isotherms does not account well for the breakthrough curves nor for the overloaded band profiles measured for mixtures of nortriptyline and amytriptyline. However, it was possible to derive adjusted parameters of a competitive tri-Langmuir model based on the fitting of the adsorption data obtained for these mixtures. A very good agreement was then found between the calculated and the experimental overloaded band profiles of all the mixtures injected.     

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.     

Comparison of the binary equilibrium isotherms of the 1-indanol enantiomers on three high-performance liquid chromatography columns of different sizes

The competitive isotherm data for the enantiomers of 1-indanol were measured on three columns, a microbore column (15 cm x 0.1 cm), a conventional analytical column (15 cm x 0.46 cm), and a semi-preparative column (20 cm x 1.0 cm), packed with Chiralcel OB. The sets of isotherm data measured on each one of these three columns could be fitted well by a bi-Langmuir isotherm model. The experimental elution band profiles of mixtures of the 1-indanol isomers were recorded on the three columns. The isotherm model, combined with the equilibrium dispersive model of chromatography, gave calculated profiles that are in excellent agreement with the experimental profiles in all cases investigated. It was found that the value of the inner diameter of the column is an important parameter in the calculation of the isotherm parameters from the measured isotherm data. In order to use isotherm data obtained on one column to account for the phase equilibrium on another one, the inner diameters of these columns must be measured accurately. The diameters of the three columns were all slightly off their nominal value. Without correction, an important systematic error was made on the isotherm data obtained with the microbore column while only negligible errors were made on the data obtained with the other two columns. After due correction for this effect, the relative difference between the isotherm data for the microbore and the semi-preparative column is still, on the average, about 10%, a difference that might be explained by the limited precision of the measurement of the microbore column diameter. The relative difference between the isotherm data for the analytical and the semi-preparative columns was about 1%, a reasonable value, since the two columns came from different batches of the same packing material.     

Adsorption of the enantiomers of 3-chloro-1-phenyl-propanol on silica-bonded chiral quinidine carbamate

The interactions of 3-chloro-1-phenyl-propanol with a quinidine carbamate-bonded chiral stationary phase under NPLC conditions were studied by measuring the adsorption isotherm data of its enantiomers by frontal analysis, modeling these data with a suitable isotherm model, and comparing the experimental overloaded elution band profiles with those calculated with this isotherm and the equilibrium dispersive model of liquid chromatography. The affinity energy distribution was calculated from the adsorption isotherm data. The results show that the surface of the adsorbent is heterogeneous and exhibits a bimodal adsorption energy distribution. This fact is interpreted in terms of the presence of two different types of adsorption sites on the stationary phase, nonselective and enantioselective sites. Albeit the bi-Langmuir isotherm model successfully accounts for the single-component data corresponding to both enantiomers, the competitive bi-Langmuir isotherm model does not allow an accurate prediction of the overloaded band profiles of the racemic mixture. Thermodynamic data are drawn for explanation. Some aspects of the retention mechanism are discussed in the light of the data obtained.     

Numerical determination of competitive adsorption isotherm of mandelic acid enantiomers on cellulose-based chiral stationary phase

The use of inverse method for the determination of competitive adsorption isotherm of mandelic acid enantiomers on cellulose tris(3,5-diethylphenyl carbamate) stationary phase is proposed in this work. Non-dominated sorting genetic algorithm with jumping genes (NSGA-II-JG) was applied to acquire the isotherm parameters by minimizing the sum of square deviations of the model predictions from the measured elution profiles. Three different competitive isotherm models, i.e., Langmuir, biLangmuir and Tóth, combined with transport-dispersive chromatographic model were used in predicting the elution profiles. Orthogonal collocation on finite element (OCFE) method was applied to obtain the calculated elution profiles. Results indicate that biLangmuir isotherm and Tóth isotherm give remarkably similar equilibrium isotherms within the investigated liquid concentration range. Band profiles calculated from both isotherm models are in good agreement with the experimental data. The validity of the determined parameters was verified by comparing the model predictions with experimental elution profiles at various experimental conditions.     

Band splitting in overloaded isocratic elution chromatography II. New competitive adsorption isotherms

The equations of two new binary competitive isotherms models are derived. The first of these models assumes that the isotherms of the two pure, single compounds have distinct monolayer capacities. Its derivation is based on kinetic arguments. The ideal adsorbed solution (IAS) framework was applied to derive the second model that is a thermodynamically consistent competitive isotherm. This second model predicts the competitive adsorption isotherm behavior of a mixture of two compounds that have single-component adsorption behavior following a BET and/or a Langmuir isotherms. Both models apply well to the binary adsorption of ethylbenzoate and 4-tert.-butylphenol on a Kromasil-C18 column (with methanol-water, 62:38, v/v, as the mobile phase). The best single-solute adsorption isotherms of these two compounds are the liquid-solid extended multilayer BET and the Langmuir isotherms, respectively. The kinetic and thermodynamic new competitive models were compared, regarding the accuracy of their prediction of the elution band profiles of mixtures of these two compounds. A better agreement between experimental and calculated profiles was observed with the kinetic model. The IAS model failed because the behavior of the ethylbenzoate/4-tert.-butylphenol adsorbed phase mixture is probably non-ideal. The most striking result is the qualitative prediction by these models of the peak splitting of 4-tert.-butylphenol during its elution in presence of ethylbenzoate.     

Study of the competitive isotherm model and the mass transfer kinetics for a BET binary system

The competitive adsorption behavior of the binary mixture of phenetole (ethoxy-benzene) and propyl benzoate in a reversed-phase system was investigated. The adsorption equilibrium data of the single-component systems were acquired by frontal analysis. The same data for binary mixtures were acquired by the perturbation method. For both compounds, the single-component isotherm data fit best to the multilayer BET model. The experimental overloaded band profiles are in excellent agreement with the profiles calculated with either the general rate model or the modified transport-dispersive models. The competitive adsorption data were modeled using the ideal adsorbed solution (IAS) theory. The numerical values of the coefficients were derived by fitting the retention times of the perturbation pulses to those calculated using the IAS theory compiled with the coherence conditions. Finally, the elution profiles of binary mixtures were recorded. They compared very well with those calculated. As a characteristic feature of this case, an unusual retainment effect of the chromatographic band of the more retained component by the less retained one was observed. The combination of the General Rate Model and the adsorption isotherm model allowed an accurate prediction of the band profiles.     

Determination of the single component and competitive adsorption isotherms of the 1-indanol enantiomers by the inverse method

The inverse method of isotherm determination consists in calculating the numerical values of the coefficients of an isotherm model that give a set of chromatographic profiles in best possible agreement with the set of experimental profiles available. This method was applied to determine the adsorption isotherms of the 1-indanol enantiomers on a cellulose tribenzoate chiral stationary phase. Both single-component and competitive isotherms were determined by using no more than one or two overloaded band profiles. The isotherms determined from the overloaded band profiles agreed extremely well with the isotherms determined by frontal analysis. Several isotherm models were used and tested. The best-fit isotherm was selected by means of statistical evaluation of the results. The results show that the adsorption is best characterized with a model describing heterogeneous adsorption with bimodal adsorption energy distribution.     

Numerical determination of the competitive isotherm of enantiomers

A numerical method was developed and used to determine adsorption isotherms in chromatography. The numerical parameters of an isotherm model were derived from the recorded band profiles of the racemic mixture of the 1-phenyl-1-propanol enantiomers, by means of a nonlinear least-squares method. We used the equilibrium-dispersive model of chromatography with several isotherm models. The numerical constants of the isotherm models were tuned so that the calculated and the measured band profiles match as much as possible. We show that this numerical inverse method can be applied even without the knowledge of the individual band profile of the pure enantiomers. The isotherms determined from the--usually unresolved--overloaded band profiles matched extremely well the isotherms determined by frontal analysis. Several isotherm models were used and tested--such as Langmuir, biLangmuir, Tóth, Langmuir-Freundlich. The best-fit isotherm was selected by means of statistical evaluation of the results.     

Repeatability and reproducibility of high concentration data in reversed-phase liquid chromatography. I. Overloaded band profiles on Kromasil-C18

Single-component adsorption-isotherm data were acquired by frontal analysis (FA) for six low-molecular-mass compounds (phenol, aniline, caffeine, theophylline, ethylbenzene and propranolol) on one Kromasil-C18 column, using water-methanol solutions (between 70:30 and 20:80, v/v) as the mobile phase. Propranolol data were also acquired using an acetate buffer (0.2 M) instead of water. The data were modeled for best agreement between calculated and experimental overloaded band profiles. The adsorption energy distribution was also derived and used for the selection of the best isotherm model. Widely different isotherm models were found to model best the data obtained for these compounds, convex upward (i.e. Langmuirian), convex downward (i.e. anti-Langmuirian), and S-shaped isotherms. Using the same sample size for all columns (loading factor, Lf approximately 10%), overloaded band profiles were recorded on four different columns packed with the same batch of Kromasil-C18 and five other columns packed with different batches of Kromasil-C18. These experimental band profiles were compared to the profile calculated from the isotherm measured by FA on the first column. The repeatability as well as the column-to-column and the batch-to-batch reproducibilities of the band profiles are better than 4%.     

Determination of competitive isotherms of enantiomers by a hybrid inverse method using overloaded band profiles and the periodic state of the simulated moving-bed process

A procedure for determination of adsorption isotherms in simulated moving-bed (SMB) chromatography is presented. The parameters of a prescribed adsorption isotherm model and rate constants are derived using a hybrid inverse method, which incorporates overloaded band profiles of the racemic mixture and breakthrough data from a single frontal experiment. The latter are included to reduce the uncertainty on the estimated saturation capacity, due to the dilution of the chromatograms with respect to the injected concentrations. The adsorption isotherm model is coupled with an axially dispersed flow model with finite mass-transfer rate to describe the experimental band profiles. The numerical constants of the isotherm model are tuned so that the calculated and measured band profiles match as much as possible. The accuracy of the isotherm model is then checked against the cyclic steady state (CSS) of the target SMB process, which is readily and cheaply obtained experimentally on a single-column set-up. This experiment is as expensive and time consuming as just a few breakthrough experiments. If necessary, the isotherm parameters are adjusted by applying the inverse method to the experimental CSS concentration profile. The method is successfully applied to determine the adsorption isotherms of Tr?gers base enantiomers on Chiralpak AD/methanol system. The results indicate that the proposed inverse method offers a reliable and quick approach to determine the competitive adsorption isotherms for a specific SMB separation.     

Modeling of the adsorption behavior and the chromatographic band profiles of enantiomers : Behavior of methyl mandelate on immobilized cellulose

The adsorption isotherms of (−)- and (+)-methyl mandelate from a hexane-isopropanol (90:10) solution were measured on a chromatographic column packed with 4-methylcellulose tribenzoate coated on silica. These isotherms are accounted for by a bi-Langmuir isotherm model, the two Langmuir terms having widely different initial slopes and saturation capacities, but each term having the same saturation capacity for the two enantiomers. The competitive isotherms were also measured. They are in excellent agreement with the prediction of a competitive bi-Langmuir model based on the single-component isotherms. The individual band profiles are in agreement with the profiles calculated from these isotherms. Thus, a simplified competitive isotherm can be used to model a separation on a chiral stationary phase the recognition mechanism of which is not well identified and the adsorption behavior of which is certainly not ideal.     

New thermodynamically consistent competitive adsorption isotherm in RPLC

A new equation of competitive isotherms was derived in the framework of the ideal adsorbed solution (IAS) that predicts multisolute adsorption isotherms from single-solute isotherms. The IAS theory makes this new isotherm thermodynamically consistent, whatever the saturation capacities of these single-component isotherms. On a Kromasil-C(18) column, with methanol-water (80/20 v/v) as the mobile phase, the best single-solute adsorption isotherm of both toluene and ethylbenzene is the liquid-solid extended multilayer BET isotherm. Despite a significant difference between the monolayer capacities of toluene (370 g/l) and ethylbenzene (170 g/l), the experimental adsorption data fit very well to single-component isotherms exhibiting the same capacities (200 g/l). The new competitive model was used for the modeling of the elution band profiles of mixtures of the two compounds. Excellent agreement between experimental and calculated profiles was observed, suggesting that the behavior of the toluene-ethylbenzene adsorbed phase on the stationary phase is close to ideal. For example, the concentrations measured for the intermediate plateau obtained in frontal analysis differ by less than 2% from those predicted by the IAS model.     

Enantioselective capillary electrophoresis for the assessment of CYP3A4-mediated ketamine demethylation and inhibition in vitro

Enantioselective CE with sulfated cyclodextrins as chiral selectors was used to determine the CYP3A4-catalyzed N-demethylation kinetics of ketamine to norketamine and its inhibition in the presence of ketoconazole in vitro. Ketamine, a chiral phencyclidine derivative, was incubated with recombinant human CYP3A4 from a baculovirus expression system as racemic mixture and as single enantiomer. Alkaline liquid/liquid extracts of the samples were analyzed with a pH 2.5 buffer comprising 50 mM Tris and phosphoric acid together with either multiple isomer sulfated β-cyclodextrin (10 mg/mL) or highly sulfated γ-cyclodextrin (2%, w/v). Data obtained in the absence of ketoconazole revealed that the N-demethylation occurred stereoselectively with Michaelis-Menten (incubation of racemic ketamine) and Hill (separate incubation of single enantiomers) kinetics. Data generated in the presence of ketoconazole as the inhibitor could best be fitted to a one-site competitive model and inhibition constants were calculated using the equation of Cheng and Prusoff. No stereoselective difference was observed, but inhibition constants for the incubation of racemic ketamine were found to be larger compared with those obtained with the incubation of single ketamine enantiomers.     

Chromatographic resolution and isotherm determination of (R,S)-mandelic acid on Chiralcel-OD column

Resolution of racemic mandelic acid ((R,S)-MA) and numerical determination of binary competitive isotherm of (R,S)-MA on Chiralcel-OD column have been investigated in this study. The effects of the alcohol modifier and acidic additive in the mobile phase on the retention and enantioseparation of (R,S)-MA were studied at first. The inverse method was then used to determine the competitive isotherm parameters of (R,S)-MA by minimizing the sum of square deviations of the model predictions from the measured elution profiles. The results indicate that the mobile phase with 85% hexane/15% isopropanol/0.3% trifluoroacetic acid mixture gives the best resolution of (R,S)-MA and competitive-modified Langmuir isotherm provides the more accurate sorption mechanism of (R,S)-MA on the cellulose tris(3,5-dimethylphenylcarbamate) chiral stationary phase.