Molecular-Level Comparison of Alkylsilane and Polar-Embedded Reversed-Phase Liquid Chromatography Systems

Organic-solvent ditch overlap in reversed-phase liquid chromatography: A molecular dynamics simulation study in cylindrical 6-12 nm-diameter pores

Mass transport through the mesopore space of a reversed-phase liquid chromatography (RPLC) column depends on the properties of the chromatographic interface, particularly on the extent of the organic-solvent ditch that favors the analyte surface diffusivity. Through molecular dynamics simulations in cylindrical RPLC mesopore models with pore diameters between 6 and 12 nm we systematically trace the evolution of organic-solvent ditch overlap due to spatial confinement in the mesopore space of RPLC columns for small-molecule separations. Each pore model of a silica-based, endcapped, C18-stationary phase is equilibrated with two mobile phases of comparable elution strength, namely 70/30 (v/v) water/acetonitrile and 60/40 (v/v) water/methanol, to consider the influence of the mobile-phase composition on the onset of organic-solvent ditch overlap. The simulations show that, as the pore diameter decreases from 9 to 6 nm, the bonded-phase density extends and compacts towards the pore center, which leads to increased accumulation of organic-solvent excess and thus enhanced organic-solvent diffusivity in the ditch. Because the acetonitrile ditch is more pronounced than the methanol ditch, acetonitrile ditch overlap sets in at less severe spatial confinement than methanol ditch overlap. The pore-averaged methanol and acetonitrile diffusivities are considerably raised by ditch overlap in the 6 nm-diameter pore, but also benefit from the ditch (without overlap) in the 7 to 12 nm-diameter pores, whereby local and pore-averaged effects are generally larger for acetonitrile than methanol.

The Solvation Shell of Small Solutes in Aqueous-Organic Solvent Mixtures and Its Implications for Reversed-Phase Liquid Chromatography

Reversed-phase liquid chromatography (RPLC) operates with water-organic solvent (W-OS) mobile phases where preferential solvation (PS) of solutes is likely. To investigate the relevance of the solute solvation shell in the mobile phase for RPLC retention, we combine data from molecular dynamics simulations of small, neutral solutes (six analytes and two dead time markers) in W-methanol (MeOH) and W-acetonitrile (ACN) mixtures with corresponding retention data obtained on an RPLC column over a wide range of W/OS ratios. Data derived from Kirkwood-Buff integrals show PS by the OS for analytes vs low or negative PS for dead time markers. W-ACN mixtures generate a higher amount of PS than W-MeOH mixtures, which contributes to the higher eluent strength of ACN in RPLC. Difference spatial distribution functions reveal anisotropic solvation shells with OS excess at hydrocarbon elements and W excess at functional groups, predicting that retention by the hydrophobic stationary phase is favored by hydrocarbon elements and limited by functional groups. Analysis of solute-solvent hydrogen bonds pinpoints the hydrogen-bond requirements toward W as the retention-limiting factor. The relation between the solute solvation shell and retention confirms the importance of W-OS and solute-W hydrogen bonding for RPLC retention.

Atomistic Details of Peptide Reversed-Phase Liquid Chromatography from Molecular Dynamics Simulations

Peptide separations by reversed-phase liquid chromatography (RPLC) are an integral part of bottom-up proteomics. These separations typically employ C18 columns with water/acetonitrile gradient elution in the presence of formic acid. Despite the widespread use of such workflows, the exact nature of peptide interactions with the stationary and mobile phases is poorly understood. Here, we employ microsecond molecular dynamics (MD) simulations to uncover details of peptide RPLC. We examined two tryptic peptides, a hydrophobic and a hydrophilic species, in a slit pore lined with C18 chains that were grafted onto SiO₂ support. Our simulations explored peptide trapping, followed by desorption and elution. Trapping in an aqueous mobile phase was initiated by C18 contacts with Lys butyl moieties. This was followed by extensive anchoring of nonpolar side chains (Leu/Ile/Val) in the C18 layer. Exposure to water/acetonitrile triggered peptide desorption in a stepwise fashion; charged sites close to the termini were the first to lift off, followed by the other residues. During water/acetonitrile elution, both peptides preferentially resided close to the pore center. The hydrophilic peptide exhibited no contacts with the stationary phase under these conditions. In contrast, the hydrophobic species underwent multiple transient Leu/Ile/Val binding interactions with C18 chains. These nonpolar interactions represent the foundation of differential peptide retention, in agreement with the experimental elution behavior of the two peptides. Extensive peptide/formate ion pairing was observed in water/acetonitrile, particularly at N-terminal sites. Overall, this work uncovers an unprecedented level of RPLC molecular details, paving the way for MD simulations as a future tool for improving retention prediction algorithms and for the design of novel column materials.

Mobile-Phase Contributions to Organic-Solvent Excess Adsorption and Surface Diffusion in Reversed-Phase Liquid Chromatography

Fast transport of retained analytes in reversed-phase liquid chromatography occurs through surface diffusion in the organic-solvent (OS)-enriched interfacial "ditch" region between the hydrophobic stationary phase and the water (W)-OS mobile phase. Through molecular dynamics simulations that recover the OS excess adsorption isotherms of a typical C₁₈-stationary phase for methanol and acetonitrile, we explore the relation between OS properties, OS excess adsorption, and surface diffusion. The emerging molecular-level picture attributes the mobile-phase contribution to surface diffusion to the hydrogen-bond capability and the eluting power of the OS. The higher affinity of methanol for the formation of W-OS hydrogen bonds at the soft, hydrophobic surface presented by the bonded-phase (C₁₈) chains reduces the OS excess and the related viscosity drop in the ditch. The lower eluting power of methanol, however, translates to increased bonded-phase contacts for analytes, which can increase their mobility gain from surface diffusion above the gain observed with acetonitrile.

Volume and composition of semi-adsorbed stationary phases in hydrophilic interaction liquid chromatography. Comparison of water adsorption in common stationary phases and eluents

Pycnometric and homologous series retention methods are used to determine the volume and mean composition of the water-rich layers partially adsorbed on the surface of several hydrophilic interaction liquid chromatography (HILIC) column fillings with acetonitrile-water and methanol-water as eluents. The findings obtained in this work confirm earlier studies using direct methods for measuring the stationary phase water content performed by Jandera's and Irgum's research groups. Water is preferentially adsorbed on the surface of the HILIC bonded phase in hydroorganic eluents containing more than 40% acetonitrile or 70% methanol, and a gradient of several water-rich transition layers between the polar bonded phase and the poorly polar bulk mobile phase is formed. These layers of reduced mobility act as HILIC stationary phases, retaining polar solutes. The volume of these layers and concentration of adsorbed water is much larger for acetonitrile-water than for methanol-water mobile phases. In hydroorganic eluents with less than 20-30% acetonitrile or 40% methanol the amount of preferentially adsorbed water is very small, and the observed retention behavior is close to the one in reversed-phase liquid chromatography (RPLC). In eluents with intermediate acetonitrile-water or methanol-water compositions a mixed HILIC-RPLC behavior is presented. Comparison of several HILIC columns shows that the highest water enrichment in the HILIC retention region for acetonitrile-water mobile phases is observed for zwitterionic and aminopropyl bonded phases, followed in minor grade for diol and polyvinyl alcohol functionalizations. Pentafluorophenyl bonded phase, usually considered a HILIC column, does not show significant water adsorption, nor HILIC retention.

Consequences of Cylindrical Pore Geometry for Interfacial Phenomena in Reversed-Phase Liquid Chromatography

The interfacial phenomena behind analyte separation in a reversed-phase liquid chromatography column take place nearly exclusively inside the silica mesopores. Their cylindrical geometry can be expected to shape the properties of the chromatographic interface with consequences for the analyte density distribution and diffusivity. To investigate this topic through molecular dynamics simulations, we introduce a cylindrical pore inside a slit pore configuration, where the inner curved and outer planar silica surface bear the same bonded phase. The present model replicates an average-sized (9 nm) mesopore in an endcapped C₁₈ column equilibrated with a mobile phase of 70/30 (v/v) water/acetonitrile. Simulations performed for ethylbenzene and acetophenone show that the surface curvature shifts the bonded phase and analyte density toward the pore center, decreases the solvent density in the bonded-phase region, increases the acetonitrile excess in the interfacial region, and considerably enhances the surface diffusivity of both analytes. Overall, the cylindrical pore provides a more hydrophobic environment than the slit pore. Ethylbenzene density is decidedly increased in the cylindrical pore, whereas acetophenone density is nearly equally distributed between the cylindrical and slit pore. The cylindrical pore geometry thus sharpens the discrimination between the apolar and moderately polar analytes while enhancing the mass transport of both.

Emerging or Underestimated Silica-Based Stationary Phases in Liquid Chromatography

Several newly synthesized or forgotten silica-based stationary phases proposed for liquid chromatography are described, including non-endcapped, short-chain alkyl phases; hydrophilic and polar-endcapped stationary phases; polar-embedded alkyl phases; long-chain alkyl phases. Stationary phases with aromatic, cyanopropyl, diol and aminopropyl functionalities are also reviewed. Stationary phases of particular interest are biomolecular materials - based on immobilized cholesterol, aminoacids, peptides, proteins or lipoproteins. Packing materials involving macrocyclic chemistry (crown ethers; calixarenes; aza-macrocycles; oligo-and polysaccharides including these of marine origin - chitin- or chitosan-based; macrocyclic antibiotics) are discussed. Since many stationary phases developed for one type of applications (e.g. chiral separation) have been found useful in solving other analytical problems (e.g. drug's plasma protein binding ability), it seemed reasonable to discuss particular chemistries behind the stationary phases presented in this review rather than specific types of interactions or chromatographic modes.

Column Characterization and Selection Systems in Reversed-Phase High-Performance Liquid Chromatography

Reversed-phase high-performance liquid chromatography (RP-HPLC) is the most popular chromatographic mode, accounting for more than 90% of all separations. HPLC itself owes its immense popularity to it being relatively simple and inexpensive, with the equipment being reliable and easy to operate. Due to extensive automation, it can be run virtually unattended with multiple samples at various separation conditions, even by relatively low-skilled personnel. Currently, there are >600 RP-HPLC columns available to end users for purchase, some of which exhibit very large differences in selectivity and production quality. Often, two similar RP-HPLC columns are not equally suitable for the requisite separation, and to date, there is no universal RP-HPLC column covering a variety of analytes. This forces analytical laboratories to keep a multitude of diverse columns. Therefore, column selection is a crucial segment of RP-HPLC method development, especially since sample complexity is constantly increasing. Rationally choosing an appropriate column is complicated. In addition to the differences in the primary intermolecular interactions with analytes of the dispersive (London) type, individual columns can also exhibit a unique character owing to specific polar, hydrogen bond, and electron pair donor-acceptor interactions. They can also vary depending on the type of packing, amount and type of residual silanols, "end-capping", bonding density of ligands, and pore size, among others. Consequently, the chromatographic performance of RP-HPLC systems is often considerably altered depending on the selected column. Although a wide spectrum of knowledge is available on this important subject, there is still a lack of a comprehensive review for an objective comparison and/or selection of chromatographic columns. We aim for this review to be a comprehensive, authoritative, critical, and easily readable monograph of the most relevant publications regarding column selection and characterization in RP-HPLC covering the past four decades. Future perspectives, which involve the integration of state-of-the-art molecular simulations (molecular dynamics or Monte Carlo) with minimal experiments, aimed at nearly "experiment-free" column selection methodology, are proposed.

Recent development in liquid chromatography stationary phases for separation of Traditional Chinese Medicine components

Traditional Chinese Medicine (TCM) is an ancient medical practice which has been used to prevent and cure diseases for thousands of years. TCMs are frequently multi-component systems with mainly unidentified constituents. The study of the chemical compositions of TCMs remains a hotspot of research. Different strategies have been developed to manage the significant complexity of TCMs, in an attempt to determine their constituents. Reversed-phase liquid chromatography (RPLC) is still the method of choice for the separation of TCMs, but has many problems related to limited selectivity. Recently, enormous efforts have been concentrated on the development of efficient liquid chromatography (LC) methods for TCMs, based on selective stationary phases. This can improve the resolution and peak capacity considerably. In addition, high-efficiency stationary phases have been applied in the analysis of TCMs since the invention of ultra high-performance liquid chromatography (UHPLC). This review describes the advances in LC methods in TCM research from 2010 to date, and focuses on novel stationary phases. Their potential in the separation of TCMs using relevant applications is also demonstrated.

A remarkable enhancement of selectivity towards versatile analytes by a strategically integrated H-bonding site containing phase

A double β-alanylated L-glutamide-derived organic phase has been newly designed and synthesized in such a way that integrated H-bonding (interaction) sites make it very suitable for the separation of versatile analytes, including shape-constrained isomers, and nonpolar, polar and basic compounds. The β-alanine residues introduced into two long-chain alkyl group moieties provide ordered polar groups through H-bonding among the amide groups.

Silica, hybrid silica, hydride silica and non-silica stationary phases for liquid chromatography

Free silanols on the surface of silica are the "villains", which are responsible for detrimental interactions of those compounds and the stationary phase (i.e., bad peak shape, low efficiency) as well as low thermal and chemical stability. For these reasons, we began this review describing new silica and hybrid silica stationary phases, which have reduced and/or shielded silanols. At present, in liquid chromatography for the majority of analyses, reversed-phase liquid chromatography is the separation mode of choice. However, the needs for increased selectivity and increased retention of hydrophilic bases have substantially increased the interest in hydrophilic interaction chromatography (HILIC). Therefore, stationary phases and this mode of separation are discussed. Then, non-silica stationary phases (i.e., zirconium oxide, titanium oxide, alumina and porous graphitized carbon), which afford increased thermal and chemical stability and also selectivity different from those obtained with silica and hybrid silica, are discussed. In addition, the use of these materials in HILIC is also reviewed.

Adsorption-induced changes of the structure of the tethered chain layers in a simple fluid

We use density functional theory to study the influence of fluid adsorption on the structure of grafted chain layer. The chains are modeled as freely jointed spheres. The chain segments and spherical molecules of the fluid interact via the Lennard-Jones potential. The fluid molecules are attracted by the substrate. We calculate the excess adsorption isotherms, the average height of tethered chains, and the force acting on selected segments of the chains. The parameters that were varied include the length of grafted chains, the grafting density, the parameters characterizing fluid-chain and fluid-surface interactions, the bulk fluid density, and temperature. We show that depending on the density of the bulk fluid the height of the bonded layer increases, remains constant, or decreases with increasing temperature.

Molecular simulation studies of reversed-phase liquid chromatography

Over the past 20 years, molecular simulation methods have been applied to the modeling of reversed-phase liquid chromatography (RPLC). The purpose of these simulations was to provide a molecular-level understanding of: (i) the structure and dynamics of the bonded phase and its interface with the mobile phase, (ii) the interactions of analytes with the bonded phase, and (iii) the retention mechanism for different analytes. However, the investigation of chromatographic systems poses significant challenges for simulations with respect to the accuracy of the molecular mechanics force fields and the efficiency of the sampling algorithms. This review discusses a number of aspects concerning molecular simulation studies of RPLC systems including the historical development of the subject, the background needed to understand the two prevalent techniques, molecular dynamics (MD) and Monte Carlo (MC) methods, and the wealth of insight provided by these simulations. Examples from the literature employing MD approaches and from the authors' laboratory using MC methods are discussed. The former can provide information on chain dynamics and transport properties, whereas the latter techniques are uniquely suited for the investigation of phase and sorption equilibria that underly RPLC retention, and both can be used to elucidate the bonded-chain conformations and solvent distributions.

The use of methyl-β-cyclodextrin to solubilize cholesterol prior to coating onto a C18 stationary phase

The use of methyl-β-cyclodextrin (MBCD) as a mobile phase additive in reversed-phase liquid chromatography is explored, with the primary goal of using MBCD to solubilize cholesterol in reversed-phase mobile phases for cholesterol-coating of C18 stationary phases. MBCD is shown to increase the solubility of cholesterol in typical reversed-phase mobile phases, especially when the stoichiometric ratio of MBCD to cholesterol exceeds 2:1. Additional equivalents of MBCD further increase solubility, or allow for weaker solvents to be used. The use of weaker solvents allows for larger coating levels of cholesterol onto a C18 stationary phase than are possible without the use of MBCD. Stationary phases coated with cholesterol using MBCD as a co-additive have different selectivity than uncoated phases, especially with regards to phenyl and shape selectivity. Further, the use of MBCD as a mobile phase additive for the elution of cholesterol is examined. It is seen via van't Hoff analysis that the reduction in retention of cholesterol when MBCD is added to the mobile phase is enthalpically driven.

Multicomponent adsorption of alcohols onto silicalite-1 from aqueous solution: isotherms, structural analysis, and assessment of ideal adsorbed solution theory

Configurational-bias Monte Carlo (CBMC) simulations in the isobaric-isothermal version of the Gibbs ensemble (GE) were carried out to probe the adsorption from aqueous solutions of methanol and/or ethanol onto silicalite-1. This methodology does require neither specification of the chemical potential nor any reference to activity models based on experimental data. The CBMC-GE methodology can be applied to the complete range of mixture compositions from pure water to pure alcohols and can also be used when multiple solute types are present at high concentration. The simulations demonstrate high selectivities for the alcohols (α(ethanol) > α(methanol)) almost over the entire composition range. The ideal adsorbed solution theory is found to substantially underpredict the amount of sorbed water and leads to very large errors for low alcohol solution concentrations. The simulations indicate that, at lower loadings, the adsorbed alcohol molecules can serve as seeds for water adsorption but, at higher loadings, alcohols displace water molecules from their preferred region.

Molecular simulations of retention in chromatographic systems: use of biased Monte Carlo techniques to access multiple time and length scales

The use of configurational-bias Monte Carlo simulations in the Gibbs ensemble allows for the sampling of phenomena that occur on vastly different time and length scales. In this review, applications of this simulation approach to probe retention in gas and reversed-phase liquid chromatographic systems are discussed. These simulations provide an unprecedented view of the retention processes at the molecular-level and show excellent agreement with experimental retention data.