A Retention-Matching Strategy for Method Transfer in Supercritical Fluid Chromatography: Introducing the Isomolar Plot Approach.

A strategy to match any retention shifts due to increased or decreased pressure drop during supercritical fluid chromatography (SFC) method transfer is presented. The strategy relies on adjusting the co-solvent molarity without the need to adjust the back-pressure regulator. Exact matching can be obtained with minimal changes in separation selectivity. To accomplish this, we introduce the isomolar plot approach, which shows the variation in molar co-solvent concentration depending on the mass fraction of co-solvent, pressure, and temperature, here exemplified by CO2-methanol. This plot allowed us to unify the effects of the co-solvent mass fraction and density on retention in SFC. The approach, which was verified on 12 known empirical retention models for each enantiomer of six basic pharmaceuticals, allowed us to numerically calculate the apparent retention factor for any column pressure drop. The strategy can be implemented either using a mechanistic approach if retention models are known or empirically by iteratively adjusting the co-solvent mass fraction. As a rule of thumb for the empirical approach, we found that the relative mass fraction adjustment needed is proportional to the relative change in the retention factor caused by a change in the pressure drop. Different proportionality constants were required to match retention in the case of increasing or decreasing pressure drops.

Impact of Methanol Adsorption on the Robustness of Analytical Supercritical Fluid Chromatography in Transfer from SFC to UHPSFC.

In supercritical fluid chromatography (SFC), the retention of a solute depends on the temperature, density, pressure, and cosolvent fraction. Here, we investigate how the adsorption of the cosolvent MeOH changes with pressure and temperature and how this affects the retention of several solutes. The lower the pressure, the stronger the MeOH adsorption to the stationary phase; in addition, at low pressure, perturbing the pressure results in significant changes in the amounts of MeOH adsorbed to the stationary phase. The robustness of the solute retention was lowest when operating the systems at low pressures, high temperatures, and low cosolvent fractions in the eluent. Here, we found a clear relationship between the sensitivity of MeOH adsorption to the stationary phase and the robustness of the separation system. Finally, we show that going from classical SFC to ultrahigh-performance SFC (UHPSFC), that is, separations conducted with much smaller packing diameters, results in retention factors that are more sensitive to fluctuations in the flow rate than with traditional SFC. The calculated density profiles indicate only a slight density drop over the traditional SFC column (3%, visualized as lighter → darker blue in the TOC), whereas the drop for the UHPSFC one was considerably larger (20%, visualized as dark red → light green in the TOC). The corresponding temperature drops were calculated to be 0.8 and 6.5 °C for the SFC and UHPSFC systems, respectively. These increased density and temperature drops are the underlying reasons for the decreased robustness of UHPSFC.

Advanced Analysis of Biosensor Data for SARS-CoV-2 RBD and ACE2 Interactions.

The traditional approach for analyzing interaction data from biosensors instruments is based on the simplified assumption that also larger biomolecules interactions are homogeneous. It was recently reported that the human receptor angiotensin-converting enzyme 2 (ACE2) plays a key role for capturing SARS-CoV-2 into the human target body, and binding studies were performed using biosensors techniques based on surface plasmon resonance and bio-layer interferometry. The published affinity constants for the interactions, derived using the traditional approach, described a single interaction between ACE2 and the SARS-CoV-2 receptor binding domain (RBD). We reanalyzed these data sets using our advanced four-step approach based on an adaptive interaction distribution algorithm (AIDA) that accounts for the great complexity of larger biomolecules and gives a two-dimensional distribution of association and dissociation rate constants. Our results showed that in both cases the standard assumption about a single interaction was erroneous, and in one of the cases, the value of the affinity constant KD differed more than 300% between the reported value and our calculation. This information can prove very useful in providing mechanistic information and insights about the mechanism of interactions between ACE2 and SARS-CoV-2 RBD or similar systems.

Analytical and preparative separation of phosphorothioated oligonucleotides: columns and ion-pair reagents.

Oligonucleotide drugs represent an emerging area in the pharmaceutical industry. Solid-phase synthesis generates many structurally closely related impurities, making efficient separation systems for purification and analysis a key challenge during pharmaceutical drug development. To increase the fundamental understanding of the important preparative separation step, mass-overloaded injections of a fully phosphorothioated 16mer, i.e., deoxythymidine oligonucleotide, were performed on a C18 and a phenyl column. The narrowest elution profiles were obtained using the phenyl column, and the 16mer could be collected with high purity and yield on both columns. The most likely contribution to the successful purification was the quantifiable displacement of the early-eluting shortmers on both columns. In addition, the phenyl column displayed better separation of later-eluting impurities, such as the 17mer impurity. The mass-overloaded injections resulted in classical Langmuirian elution profiles on all columns, provided the concentration of the ion-pairing reagent in the eluent was sufficiently high. Two additional column chemistries, C4 and C8, were also investigated in terms of their selectivity and elution profile characteristics for the separation of 5-20mers fully phosphorothioated deoxythymidine oligonucleotides. When using triethylamine as ion-pairing reagent to separate phosphorothioated oligonucleotides, we observed peak broadening caused by the partial separation of diastereomers, predominantly seen on the C4 and C18 columns. When using the ion-pair reagent tributylamine, to suppress diastereomer separation, the greatest selectivity was found using the phenyl column followed by C18. The present results will be useful when designing and optimizing efficient preparative separations of synthetic oligonucleotides.

Investigation of factors influencing the separation of diastereomers of phosphorothioated oligonucleotides.

This study presents a systematic investigation of factors influencing the chromatographic separation of diastereomers of phosphorothioated pentameric oligonucleotides as model solutes. Separation was carried out under ion-pairing conditions using an XBridge C18 column. For oligonucleotides with a single sulfur substitution, the diastereomer selectivity was found to increase with decreasing carbon chain length of the tertiary alkylamine used as an ion-pair reagent. Using an ion-pair reagent with high selectivity for diastereomers, triethylammonium, it was found the selectivity increased with decreased ion-pair concentration and shallower gradient slope. Selectivity was also demonstrated to be dependent on the position of the modified linkage. Substitutions at the center of the pentamer resulted in higher diastereomer selectivity compared to substitutions at either end. For mono-substituted oligonucleotides, the retention order and stereo configuration were consistently found to be correlated, with Rp followed by Sp, regardless of which linkage was modified. The type of nucleobase greatly affects the observed selectivity. A pentamer of cytosine has about twice the diastereomer selectivity of that of thymine. When investigating the retention of various oligonucleotides eluted using tributylammonium as the ion-pairing reagent, no diastereomer selectivity could be observed. However, retention was found to be dependent on both the degree and position of sulfur substitution as well as on the nucleobase. When analyzing fractions collected in the front and tail of overloaded injections, a significant difference was found in the ratio between Rp and Sp diastereomers, indicating that the peak broadening observed when using tributylammonium could be explained by partial diastereomer separation.

Reliable Strategy for Analysis of Complex Biosensor Data.

When using biosensors, analyte biomolecules of several different concentrations are percolated over a chip with immobilized ligand molecules that form complexes with analytes. However, in many cases of biological interest, e.g., in antibody interactions, complex formation steady-state is not reached. The data measured are so-called sensorgram, one for each analyte concentration, with total complex concentration vs time. Here we present a new four-step strategy for more reliable processing of this complex kinetic binding data and compare it with the standard global fitting procedure. In our strategy, we first calculate a dissociation graph to reveal if there are any heterogeneous interactions. Thereafter, a new numerical algorithm, AIDA, is used to get the number of different complex formation reactions for each analyte concentration level. This information is then used to estimate the corresponding complex formation rate constants by fitting to the measured sensorgram one by one. Finally, all estimated rate constants are plotted and clustered, where each cluster represents a complex formation. Synthetic and experimental data obtained from three different QCM biosensor experimental systems having fast (close to steady-state), moderate, and slow kinetics (far from steady-state) were evaluated using the four-step strategy and standard global fitting. The new strategy allowed us to more reliably estimate the number of different complex formations, especially for cases of complex and slow dissociation kinetics. Moreover, the new strategy proved to be more robust as it enables one to handle system drift, i.e., data from biosensor chips that deteriorate over time.

Optimizing Column Length and Particle Size in Preparative Batch Chromatography Using Enantiomeric Separations of Omeprazole and Etiracetam as Models: Feasibility of Taguchi Empirical Optimization.

The overreaching purpose of this study is to evaluate new approaches for determining the optimal operational and column conditions in chromatography laboratories, i.e., how best to select a packing material of proper particle size and how to determine the proper length of the column bed after selecting particle size. As model compounds, we chose two chiral drugs for preparative separation: omeprazole and etiracetam. In each case, two maximum allowed pressure drops were assumed: 80 and 200 bar. The processes were numerically optimized (mechanistic modeling) with a general rate model using a global optimization method. The numerical predictions were experimentally verified at both analytical and pilot scales. The lower allowed pressure drop represents the use of standard equipment, while the higher allowed drop represents more modern equipment. For both compounds, maximum productivity was achieved using short columns packed with small-particle size packing materials. Increasing the allowed backpressure in the separation leads to an increased productivity and reduced solvent consumption. As advanced numerical calculations might not be available in the laboratory, we also investigated a statistically based approach, i.e., the Taguchi method (empirical modeling), for finding the optimal decision variables and compared it with advanced mechanistic modeling. The Taguchi method predicted that shorter columns packed with smaller particles would be preferred over longer columns packed with larger particles. We conclude that the simpler optimization tool, i.e., the Taguchi method, can be used to obtain "good enough" preparative separations, though for accurate processes, optimization, and to determine optimal operational conditions, classical numerical optimization is still necessary.

Thermodynamic and kinetic approaches for evaluation of monoclonal antibody - Lipoprotein interactions.

Two complementary instrumental techniques were used, and the data generated was processed with advanced numerical tools to investigate the interactions between anti-human apoB-100 monoclonal antibody (anti-apoB-100 Mab) and apoB-100 containing lipoproteins. Partial Filling Affinity Capillary Electrophoresis (PF-ACE) combined with Adsorption Energy Distribution (AED) calculations provided information on the heterogeneity of the interactions without any a priori model assumptions. The AED calculations evidenced a homogenous binding site distribution for the interactions. Quartz Crystal Microbalance (QCM) studies were used to evaluate thermodynamics and kinetics of the Low-Density Lipoprotein (LDL) and anti-apoB-100 Mab interactions. High affinity and selectivity were observed, and the emerging data sets were analysed with so called Interaction Maps. In thermodynamic studies, the interaction between LDL and anti-apoB-100 Mab was found to be predominantly enthalpy driven. Both techniques were also used to study antibody interactions with Intermediate-Density (IDL) and Very Low-Density (VLDL) Lipoproteins. By screening affinity constants for IDL-VLDL sample in a single injection we were able to distinguish affinity constants for both subpopulations using the numerical Interaction Map tool.

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.

Combining Chemometric Models with Adsorption Isotherm Measurements to Study Omeprazole in RP-LC.

The adsorption of the proton-pump inhibitor omeprazole was investigated using RP-LC with chemometric models combined with adsorption isotherm modelling to study the effect of pH and type of organic modifier (i.e., acetonitrile or methanol). The chemometric approach revealed that omeprazole was tailing with methanol and fronting with acetonitrile along with increased fronting at higher pH. The increased fronting with higher pH for acetonitrile was explored using a pH-dependent adsorption isotherm model that was determined using the inverse method and it agreed well with the experimental data. The model indicated that the peaks exhibit more fronting at high pH due to a larger fraction of charged omeprazole molecules. This model could accurately predict the shape of elution profiles at arbitrary pH levels in the studied interval. Using a two-layer adsorption isotherm model, the difference between acetonitrile and methanol was studied at the lowest pH at which almost all omeprazole molecules are neutral. Omeprazole had adsorbate-adsorbate interactions that were similar in strength for the acetonitrile and methanol mobile phases, while the solute-adsorbent interactions were almost twice as strong with methanol. The difference in the relative strengths of these two interactions likely explains the different peak asymmetries (i.e., tailing/fronting) in methanol and acetonitrile. In conclusion, thermodynamic modelling can complement chemometric modeling in HPLC method development and increase the understanding of the separation.

Partial filling affinity capillary electrophoresis including adsorption energy distribution calculations--towards reliable and feasible biomolecular interaction studies.

In this work, a method to study and analyze the interaction data in free solution by exploiting partial filling affinity capillary electrophoresis (PF-ACE) followed by adsorption energy distribution calculations (AED) prior model fit to adsorption isotherms will be demonstrated. PF-ACE-AED approach allowed the possibility to distinguish weak and strong interactions of the binding processes between the most common apolipoprotein E protein isoforms (apoE2, apoE3, apoE4) of high density lipoprotein (HDL) and apoE-containing HDL2 with major glycosaminoglycan (GAG) chain of proteoglycans (PGs), chondroitin-6-sulfate (C6S). The AED analysis clearly revealed the heterogeneity of the binding processes. The major difference was that they were heterogeneous with two different adsorption sites for apoE2 and apoE4 isoforms, whereas interestingly for apoE3 and apoE-containing HDL2, the binding was homogeneous (one site) adsorption process. Moreover, our results allowed the evaluation of differences in the binding process strengths giving the following order with C6S: apoE-containing HDL2 > apoE2 > apoE4 > apoE3. In addition, the affinity constant values determined could be compared with those obtained in our previous studies for the interactions between apoE isoforms and another important GAG chain of PGs - dermatan sulfate (DS). The success of the combination of AED calculations prior to non-linear adsorption isotherm model fit with PF-ACE when the concentration range was extended, confirmed the power of the system in the clarification of the heterogeneity of biological processes studied.

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.

Partial-filling affinity capillary electrophoresis and quartz crystal microbalance with adsorption energy distribution calculations in the study of biomolecular interactions with apolipoprotein E as interaction partner.

Adsorption energy distribution (AED) calculations were successfully applied to partial-filling affinity capillary electrophoresis (PF-ACE) to facilitate more detailed studies of biomolecular interactions. PF-ACE with AED calculations was employed to study the interactions between two isoforms of apolipoprotein E (apoE) and dermatan sulfate (DS), and a quartz crystal microbalance (QCM) was used in combination with AED calculations to examine the interactions of the 15-amino-acid peptide fragment of apoE with DS. The heterogeneity of the interactions was elucidated. Microscale thermophoresis was used to validate the results. The interactions studied are of interest because, in vivo, apolipoprotein E localizes on DS-containing regions in the extracellular matrix of human vascular subendothelium. Two-site binding was demonstrated for the isoform apoE3 and DS, but only one-site binding for apoE2-DS. Comparable affinity constants were obtained for the apoE2-DS, apoE3-D3, and 15-amino-acid peptide of apoE-DS using the three techniques. The results show that combining AED calculations with modern biosensing techniques can open up another dimension in studies on the heterogeneity and affinity constants of biological molecules.

Regeneration of a silica monolithic rod column using harsh methods followed by firm thermodynamic and kinetic validation.

In this study, a numerical tool is introduced--based on thermodynamic and kinetic separation theory--for validating the regeneration of monolithic rod columns after cutting their inlet sections. A long-used RP-18e monolithic column was deemed to be unfit for further coffee analysis because of poor separation performance. The columns brownish inlet section was physically removed with a lathe, leaving a clean white inlet section. The original and regenerated columns were extensively analyzed and compared using numerical tools for processing adsorption data. The perturbation peak method was used to measure the adsorption isotherm of phenol on the original and regenerated monolith and the adsorption energy distributions were calculated for identifying any change in the degree of heterogeneity. Although peak shapes improved considerably after regeneration, no significant differences were found in the detailed characterization of the processed adsorption data between the original column and the regenerated one. This indicates that the removal of a section of the monolithic bed can be undertaken without damaging the column and columns in which their inlet head sections are contaminated may still function with normal adsorption behavior. In addition, the combined thermodynamic and kinetic methodology could accurately be used to evaluate any regeneration method of columns.

Modelling the pH dependent retention and competitive adsorption of charged and ionizable solutes in mixed-mode and reversed-phase liquid chromatography.

This study investigated the influence of pH on the retention of solutes using a mixed-mode column with carboxyl (-COOH) groups acting as weak cation exchanger bonded to the terminal of C18 ligands (C18-WCX column) and a traditional reversed-phase C18 column. First, a model based on electrostatic theory was derived and successfully used to predict the retention of charged solutes (charged, and ionizable) as a function of mobile phase pH on a C18-WCX column. While the Horváth model predicts the pH-dependent retention of ionizable solutes in reversed-phase liquid chromatography (RPLC) solely based on solute ionization, the developed model incorporates the concept of surface potential generated on the surface of the stationary phase and its variation with pH. To comprehensively understand the adsorption process, adsorption isotherms for these solutes were individually acquired on the C18-WCX and reversed-phase C18 columns. The adsorption isotherms followed the Langmuir model for the uncharged solute and the electrostatically modified Langmuir model for charged solutes. The elution profiles for the single components were calculated from these isotherms using the equilibrium dispersion column model and were found to be in close agreement with the experimental elution profiles. To enable modelling of two-component cases involving charged solute(s), a competitive adsorption isotherm model based on electrostatic theory was derived. This model was later successfully used to calculate the elution profiles of two components for scenarios involving (a) a C18 Column: two charged solutes, (b) a C18 Column: one charged and one uncharged solute, and (c) a C18-WCX Column: two charged solutes. The strong alignment between the experimental and calculated elution profiles in all three scenarios validated the developed competitive adsorption model.

Prediction of overloaded concentration profiles under ultra-high-pressure liquid chromatographic conditions.

In this study, overloaded elution profiles under ultra-high-pressure liquid chromatographic (UHPLC) conditions and accounting for the severe pressure and temperature gradients generated, are compared with experimental data. The model system consisted of an C18 column packed with 1.7-µm particles (i.e., a UHPLC column) and the solute was 1,3,5-tri­tert-butylbenzene eluted with a mobile phase composed of 85/15 (v/v) acetonitrile/water. Two thermal modes were considered, and the solute was eluted at the very high inlet pressures necessary to achieve a highly efficient and rapid chromatographic process, as provided by using columns packed with small particles. However, the high inlet pressure and high linear velocity of the mobile phase caused the production of a significant amount of heat, and consequently, the formation of axial and radial temperature gradients. Due to these gradients, the retention and the mobile phase velocity were no longer constant. Thus, simple mathematical models consisting only of the mass balance equations are unsuitable to properly model the elution profiles. Here, the elution concentration profiles were predicted using a combined two-dimensional heat and mass transfer model, also including the calculation of the mobile phase velocity distribution. The isotherm adsorption model was the bi-Langmuir isotherm model with Henry constants that depended on the local temperature and pressure in the column. These adjustments allowed us to precisely account for changes in the shape and retention of the overloaded concentration profiles in the mobile phase. The proposed model provided accurate predictions of the overloaded concentration profiles, demonstrating good agreement with experimental profiles eluted under severe pressure and temperature gradients in the column even in the most extreme cases where the pressure drops reached 846 bar and the temperature gradients equaled 0.15 K mm-1 and 0.95 K mm-1 in the axial and the radial directions, respectively. In such cases 36 % decrease of the retention factor was observed along the column and 2 % increase in radial direction. These changes, combined with the velocity distribution, shifted the overloaded elution profile's shock towards the center of the column, advancing approximately 3 mm from its initial position close to the column wall. Ultimately, this resulted in the broadening of the elution band.

Fundamental investigation of impact of water and TFA additions in peptide sub/supercritical fluid separations.

The retention of three peptides was studied under analytical and overloaded conditions at different concentrations of trifluoroacetic acid (TFA) and water added to the co-solvent methanol (MeOH). Four columns with different stationary phase properties, i.e., silica, diol, 2-ethylpyridine and cyanopropyl (CN) columns, were evaluated in this investigation. The overall aim was to get a deeper understanding on how column chemistry as well as water and TFA in the co-solvent affect the analytical and overloaded elution profiles using multivariate design of experiments and adsorption measurements of co-solvent components. Multivariate experimental design modeling indicated that water had on average around five times higher effect on the retention than the addition of TFA. The results also showed that the retention increases with the addition of TFA and water to the co-solvent on all columns except the CN column, on which the retention decreased. When examining the effect of adding water to the co-solvent, evidence of a hydrophilic interaction liquid chromatography (HILIC)-like retention mechanism was found on the three other columns with more polar stationary phases. However, on the CN column water acted as an additive, decreasing the retention due to competition with the peptide for available adsorption surface. Adsorption isotherm measurements of the polar solvent MeOH showed that MeOH adsorbs much weaker on the CN column than on the other columns. Addition of TFA and water to the co-solvent substantially sharpened the elution profiles under both overloaded and analytical conditions. Adding a small amount of TFA (from 0 % to 0.05 %) to the co-solvent substantially improved the peak shape of the elution profiles, while further addition (from 0.05 % to 0.15 %) had only a minor effect on the elution profile shape. The reduced retention on the CN column could not be explained by TFA adsorption, which was very weak on all studied columns (retention factor, 0.05-0.15). One could therefore speculate that the ion-pairing complex formed between the peptide and TFA in the mobile phase, reduce the retention due to its reduced polarity. On the other columns displaying HILIC-like properties, the TFA probably just decreased the pH of the mobile phase, thereby promoting the partitioning of the peptide into the water-rich layer. Finally, peak deformation due to diluent-eluent mismatch was observed under overloaded conditions. This was most severe in the cases where MeOH adsorption to the stationary phase was strong and the peptides were only mildly retained. Adding 1,4-dioxan to the diluent resolved this issue.

Three complementary techniques for the clarification of temperature effect on low-density lipoprotein-chondroitin-6-sulfate interaction.

A rigorous processing of adsorption data from quartz crystal microbalance technology was successfully combined with the data obtained by partial filling affinity capillary electrophoresis and molecular dynamics for the clarification of the temperature effect on the interaction of a major glycosaminoglycan chain chondroitin-6-sulfate (C6S) of proteoglycans with low-density lipoprotein (LDL) and with a peptide fragment of apolipoprotein B-100 (residues 3359-3377 of LDL, PPBS). Two experimental techniques and computational atomistic methods demonstrated a nonlinear pattern of the affinity of C6S at temperatures above 38.0 °C to both LDL and PPBS. The temperature affects the interaction of C6S with LDL and PPBS by influencing the structural behavior of glycosaminoglycan C6S and/or that of LDL.

Sample conditions to avoid pH distortion in RP-LC.

Band deformations might take place for acids and bases in preparative applications and adsorption studies where it is necessary to use overloaded injections. In this study, we focus on how deformations can be prevented without using highly concentrated buffers that may precipitate in the eluent. We have systematically investigated how the elution zones depend on which protolytic form the analyte has when it is dissolved. Basic and acidic model compounds are investigated using eluents with different pH values and the resulting elution profiles are compared when the analytes are dissolved in their protonated and deprotonated form, i.e., in uncharged form or as different kinds of salts. Depending on the analyte's protolytic form, a sample zone is created at the column inlet whose pH deviates, more or less, from the bulk eluent's. If the local adsorption strength in this sample zone is greater than the bulk eluent's, the elution profiles are compressed. Under opposite conditions, the eluted bands are more or less deformed and may even be split; completely different deformations can even take place for different kinds of salt combinations. Explanations of these, and other, effects, together with detailed guidelines for proper sample preparation to avoid peak deformations, are given.

Evaluation of a combined linear-nonlinear approach for column characterization using modern alkaline-stable columns as model.

This study investigates if deeper understanding is achieved when combining nonlinear and linear chromatographic column characterization methods. As test systems, two hybrid columns (Phenomenex Gemini-NX C18 and Kromasil Eternity C18) and one classic one (Kromasil-C18) were selected. The nonlinear methods were based on firm adsorption theory and involved determination of adsorption isotherms followed by calculations with a new numerical tool, adsorption energy distribution, on probe components at different pH values. The linear methods involved the hydrophobic subtraction model and selected probe components retention factors as a function of pH. The combined analysis indicated that both complementary and confirmative information can be achieved regarding the actual model systems.