Molecular simulation study of the bonded-phase structure in reversed-phase liquid chromatography with neat aqueous solvent

Characterization of ribostamycin and its impurities using a nano-quantity analyte detector: Systematic comparison of performance among three different aerosol detectors

Characterization of aminoglycoside antibiotics like ribostamycin is important due to the complex composition and common toxic impurities. Aerosol detectors are often employed for determination of these non-absorbent analytes. In this work, a robust and cost-effective method was developed for simultaneous detection of ribostamycin and its related substances using high-performance liquid chromatography (HPLC) with a relative new aerosol detector named nano-quantity analyte detector (NQAD). With the introduction of less toxic but more compatible ion-pairs pentafluoropropionic acid (PFPA) and trifluoroacetic acid (TFA) in the eluent, an optimized separation effect was achieved. Compared with the other two aerosol detectors namely ELSD (evaporative light scattering detector) and CAD (charged aerosol detector), method verification and quantitative detection results revealed that NQAD had higher sensitivity than ELSD with a 0.8 μg/mL limit of detection, as well as wider linear range (from 2 μg/mL to 1000 μg/mL) than both CAD (from 2 μg/mL to 200 μg/mL) and ELSD (from 8 μg/mL to 200 μg/mL) detector. The performance of NQAD helped to realize detection of ribostamycin and its impurities with significant concentration differences in a single run. With a cation suppressor to eliminate the ion-suppression caused by the ion-pairs in the eluent, the structure of nine impurities in ribostamycin sample was characterized by liquid chromatography-mass spectrum (LC-MS). Both external standard and area normalization calculation were investigated, and NQAD obtained more accurate results due to its full-range linear response-to-concentration relationship, providing an alternative for routine quality control of multi analyte systems.

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.

The combination of partition, size exclusion, and hydrodynamic models in chromatography, and application to bonded phases on porous supports

Liquid-liquid partition chromatography has been used for many years as a model and teaching introduction to column chromatography. However, the partition model does not describe separations on bonded phases with porous supports particularly well, especially regarding the thermodynamics controlling solute distribution. Further difficulties arise when more than one mechanism is involved in solute retention. Nomenclature is not perfectly aligned with the underlying thermodynamic descriptors and is inconsistently applied over various chromatographic techniques. Presented here is a general description of retention that spans partition, size exclusion, and hydrodynamic separation processes, and is then applied to bonded-phase separations on porous supports. The model provides a general description applicable to adsorption, reversed-phase, hydrophilic interaction, size-exclusion, hydrodynamic chromatography, and any combination of these techniques including liquid chromatography at the critical condition. Further expansion to include retention by ion-exchange and field-flow fractionation appears to be possible. Recommendations on retention factor definition and evaluation are given.

Simulation and theory of open-tube dispersion in short and long capillaries with slip boundaries and retention

Using random walk techniques, high resolution simulations of zone shape are conducted in open capillary tubes for short and long tube conditions. Finite size solutes are used as tracers in this treatment. Slip flow boundary conditions and wall retention are utilized as needed. These simulations are able to reproduce previous work in short and long tubes. For the short tube case where dispersion does not asymptotically approach the classic Taylor-Aris and Golay solutions, the effect of slip flow boundaries in the transient region shows zone shapes with abbreviated tails where the larger slip flow values cause zone compression. The use of slip flow to lower dispersion in capillary-based, wall-coated separations is shown to favor long tube behavior. This is because slip flow is relevant for cases where slip lengths are fractions of small capillary tube diameters. Incorporating slip flow into transport in capillaries favors a very small capillary radius where the cross-sectional diffusion length is very small and sampling times are fast. The purely convective zone shape with slip flow boundaries is derived analytically. Applications for this type of separation, guided by both analytical theory and simulation, show the potential for nano-sized capillary tubes less than 1 μm in diameter and favor very fast isocratic separations. Using long tube retention theory with slip boundaries shows that the dispersion-reducing region is most important in the range 0 ≤ k' ≤ 1, a relatively small retention window. Further discussion of the gradient elution technique and dispersion in packed beds suggests that the general usage of slip flow boundaries is restricted in liquid phase separation systems.

Computer simulation for the study of the liquid chromatographic separation of explosive molecules

The application of high performance liquid chromatography (HPLC) to separate explosive chemicals was investigated by molecular dynamics (MD) simulations. The explosive ingredients including NG, RDX, HMX and TNT were assigned as solutes, while methanol (CH₃OH) and acetonitrile (CH₃CN) were assigned as solvents in the solution system. The polymeric-molecular siloxanes (SiC8) and poly-1,2-methylenedioxy-4-propenyl benzene (PISAF) compounds were treated as stationary phase in the simulation. The simulation results showed that the different species of explosive ingredients were separated successfully in the solutions by each of the constructed stationary phase of SiC8 and PISAF after a total simulation time of 12.0 ps approximately, which were consistent with the experimental analysis of HPLC spectra. The origin for the separation was found due to the electrostatic interactions between polymer and explosives.

Investigations into the separation behaviour of perfluorinated C8 and undecanoic acid modified silica hydride stationary phases

In this study, the surface charge properties of perfluorinated C8 (PerfluoroC8) and undecanoic acid (UDA) modified silica hydride stationary phases have been investigated. The zeta potential values of these stationary phases were measured in aqueous/acetonitrile mobile phases of different pH, buffer concentrations and acetonitrile contents. The retention behaviour of several basic, acidic and neutral compounds were then examined with these two stationary phases, with U-shaped retention dependencies evident with regard to the organic solvent content of the mobile phase. Plots of the logarithmic retention factor versus buffer concentration revealed slopes ≥ -0.41 for both stationary phases, indicating the involvement of mixed mode retention mechanisms with contributions from both ionic and non-ionic interactions. Using a linear solvation energy relationship approach, the origins of these interactions under different mobile phase conditions were differentiated and quantified. The PerfluoroC8 stationary phase exhibited stronger retention for basic compounds under high acetonitrile content mobile phase conditions, whilst stronger retention was observed for all compounds with the UDA stationary phase under high aqueous content mobile phase conditions. The more negative zeta potentials of the UDA stationary phase correlated with higher total charge density, surface charge density and charge density at the beta plane (the outer plane of the double layer) compared to the PerfluoroC8 stationary phase. With mobile phases of low buffer concentrations, more negative zeta potential values were unexpectedly observed for the PerfluoroC8 stationary phase with slight increases in the C descriptor value, reflecting also the greater accessibility of the analytes to the stationary phase surface. Comparison of the retention behaviours on these phases with other types of silica hydride stationary phases has revealed different patterns of selectivity.

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.

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.

Computational study of enantioseparation by amylose tris(3,5-dimethylphenylcarbamate)-based chiral stationary phase

The mechanism of chiral separation on amylose tris(3,5-dimethylphenylcarbamate) is studied with docking simulations of enantiomers by molecular dynamics. All-atom models of amylose tris(3,5-dimethylphenylcarbamate) on the modified silica gel surface were constructed for the docking simulations of metalaxyl and benalaxyl. The elution orders and energetic differences were also predicted based on the intermolecular interactions, which were in agreement with the experimental results. The radial distribution function was employed to analyze the structural features of the enantiomer-chiral stationary phase complex and used to elucidate the mechanism of chiral separation. The separation of metalaxyl and benalaxyl is mainly controlled by the hydrogen bond. And the binding sites had slight differences for the pair of enantiomers, but obvious differences between different chemicals.

Effects of nanopore size on properties of modified inner surfaces

By analyzing sorption isotherm curves of surface treated MCM-41 samples, it is noticed that if the nanopore size is relatively small, the end group dominates the solid-liquid interaction and the influence of the side group is relatively weak, which can be attributed to the confinement effect of nanopore walls.

Adsorption of Acridine Orange at a C8,18/Water/Acetonitrile Interface

Fully atomistic simulations are used to characterize the molecular dynamics (MD) of acridine orange (3,6-dimethylaminoacridine) at a chromatographic interface. Multiple 1 ns MD simulations were performed for acridine orange at the interface between three different acetonitrile/water mixtures (0/100, 20/80, and 50/50) with C8 and C18 alkyl chains. The diffusion coefficient, D, of acridine orange in pure solvent was found to be 4 times smaller at the water/C18 interface (D = 0.022 x 10(-4) cm2/s) than in bulk water (D = 0.087 x 10(-4) cm2/s), in qualitative agreement with experiment. Rotational reorientation times were 20 and 700 ps, which also agree favorably with the measured time scales of 130 and 740 ps. Contrary to experiment, the simulations found that for increasing surface coverage, the diffusion coefficient for acridine decreased. Detailed analysis of the solvent structure showed that the transport properties of acridine were primarily governed by the solvent distribution above the functionalized surface. The solvent structure, in turn, was largely determined by the surface consisting of the silica layer, the alkyl chains, and their functionalization.

Retention Mechanism in Reversed-Phase Liquid Chromatography: A Molecular Perspective

A detailed, molecular-level understanding of the retention mechanism in reversed-phase liquid chromatography (RPLC) has eluded analytical chemists for decades. Through validated, particle-based Monte Carlo simulations of a model RPLC system consisting of dimethyloctadecylsilanes at a coverage of 2.9 micro mol/m2 on an explicit silica substrate with unprotected residual silanols in contact with a water/methanol mobile phase, we show that the molecular-level retention processes for nonpolar and polar analytes, such as alkanes and alcohols, are much more complex than what has been previously deduced from thermodynamic and theoretical arguments. In contrast to some previous assumptions, the simulations indicate that both partitioning and adsorption play a key role in the separation process and that the stationary phase in RPLC behaves substantially different from a bulk hydrocarbon phase. The retention of nonpolar methylene segments is dominated by lipophilic interactions with the retentive phase, while solvophilic interactions are more important for the retention of the polar hydroxyl group.

Effect of temperature on gradient reequilibration in reversed-phase liquid chromatography

The effect of mobile phase modifier and temperature on gradient reequilibration is examined using three different stationary phases. The stationary phases studied are a traditional C18 phase, a polar endcapped C18 phase, and an alkyl phase with a polar embedded group. It was observed that both temperature and choice of mobile phase organic modifier had an effect on gradient reequilibration volume on both the traditional C18 stationary phase and the polar endcapped phase. On both these phases, at any given temperature, the reequilibration volume was generally smaller when methanol was used as the mobile phase modifier as compared to acetonitrile. As the temperature is increased from 10 to 50 degrees C, significant reductions in reequilibration volume were observed with both mobile phase modifiers. In contrast, neither temperature nor choice of modifier appeared to have much effect on reequilibration volume when the polar embedded group stationary phase was considered.

Selectivity in reversed-phase separations: General influence of solvent type and mobile phase pH

The influence of the mobile phase on retention is studied in this paper for a group of over 70 compounds with a broad range of multiple functional groups. We varied the pH of the mobile phase (pH 3, 7, and 10) and the organic modifier (methanol, acetonitrile (ACN), and tetrahydrofuran (THF)), using 15 different stationary phases. In this paper, we describe the overall retention and selectivity changes observed with these variables. We focus on the primary effects of solvent choice and pH. For example, transfer rules for solvent composition resulting in equivalent retention depend on the packing as well as on the type of analyte. Based on the retention patterns, one can calculate selectivity difference values for different variables. The selectivity difference is a measure of the importance of the different variables involved in method development. Selectivity changes specific to the type of analyte are described. The largest selectivity differences are obtained with pH changes.

Solvation of the Whelk-O1 chiral stationary phase: A molecular dynamics study

Density functional theory calculations and molecular dynamics simulations are employed to explore the solvation of the Whelk-O1 chiral stationary phase. First, a semi-flexible representation of the Whelk-O1 selective molecule is extracted from an extensive series of B3LYP/6-311+ G(2d,p) calculations. The resulting model is used to build a chiral surface, including end-caps, for molecular dynamics study of the interface between solvent and Whelk-O1. Three solvent environments in common use for Whelk-O1 HPLC have been examined: a normal-phase solvent of n-hexane/2-propanol; a reversed-phase solvent of water/methanol; and a supercritical solvent of CO(2) and methanol. In each case, we analyze the interface with an emphasis on solvent composition and solvent hydrogen bonding to the Whelk-O1 selector.

Chain conformation and solvent partitioning in reversed-phase liquid chromatography: Monte Carlo simulations for various water/methanol concentrations

Many structural models for the stationary phase in reversed-phase liquid chromatography (RPLC) systems have been suggested from thermodynamic and spectroscopic measurements and theoretical considerations. To provide a molecular picture of chain conformation and solvent partitioning in a typical RPLC system, a particle-based Monte Carlo simulation study is undertaken for a dimethyl octadecyl (C(18)) bonded stationary phase on a model siliceous substrate in contact with mobile phases having different methanol/water concentrations. Following upon previous simulations for gas-liquid chromatography and liquid-liquid phase equilibria, the simulations are conducted using the configurational-bias Monte Carlo method in the Gibbs ensemble and the transferable potentials for phase equilibria force field. The simulations are performed for a chain surface density of 2.9 micromol/m(2), which is a typical bonded-phase coverage for mono-functional alkyl silanes. The solvent concentrations used here are pure water, approximately 33 and 67% mole fraction of methanol and pure methanol. The simulations show that the chain conformation depends only weakly on the solvent composition. Most chains are conformationally disordered and tilt away from the substrate normal. The interfacial width increases with increasing methanol content and, for mixtures, the solvent shows an enhancement of the methanol concentration in a 10 Angstrom region outside the Gibbs dividing surface. Residual surface silanol groups are found to provide hydrogen bonding sites that lead to the formation of substrate bound water and methanol clusters, including bridging clusters that penetrate from the solvent/chain interfacial region all the way to the silica surface.

Identification of isolated cavity features within molecular dynamics simulated chromatographic surfaces

Highly ordered morphological features were characterized for molecular dynamics simulated alkyl-modified silica models that represent chromatographic materials with enhanced shape recognition capability. Deep cavities (8-10A wide) within the alkyl chains were identified for C18 polymeric models corresponding to shape-selective RPLC stationary phases. The all-trans conformational distal-end segments of these isolated cavities averaged over a 100 ps simulation time interval were observed to increase (up to 15 A) in models with an increase in both surface coverage and corresponding shape selectivity. Similar-structure cavities with significant alkyl chain ordered regions (>11A) were isolated from two independent C18 models (differing in bonding chemistry, density and temperature) that represent highly shape-selective materials. The size and depth of these ordered regions increased (up to 28 A) for the extended-length C30 alkyl phase models. These initial results offer a physical representation of alkyl-modified surfaces that may facilitate the identification of potential molecular features that may be involved in the shape-selective retentive processes, as well as illustrating the potential for such computational techniques to predict the molecular recognition capabilities of novel analyte-specific sorbents.

Adsorption mechanism in reversed-phase liquid chromatography

The effect of the bonding density of the octadecyl chains onto the same silica on the adsorption and retention properties of low molecular weight compounds (phenol, caffeine, and sodium 2-naphthalene sulfonate) was investigated. The same mobile phase (methanol:water, 20:80, v/v) and temperature (T = 298 K) were applied and two duplicate columns (A and B) from each batch of packing material (neat silica, simply endcapped or C1 phase, 0.42, 1.01, 2.03, and 3.15 micromol/m2 of C18 alkyl chains) were tested. Adsorption data of the three compounds were acquired by frontal analysis (FA) and the adsorption energy distributions (AEDs) were calculated using the expectation-maximization method. Results confirmed earlier findings in linear chromatography of a retention maximum at an intermediate bonding density. From a general point of view, the saturation capacity of the adsorbent tends to decrease with increasing bonding density, due to the vanishing space intercalated between the C18 bonded chains and to the decrease of the specific surface area of the stationary phase. The equilibrium constants are maximum for an intermediary bonding density (approximately 2 micromol/m2). An enthalpy-entropy compensation was found for the thermodynamic parameters of the isotherm data. Weak equilibrium constants (small deltaH) and high saturation capacities (large deltaS) were observed at low bonding densities, higher equilibrium constants and lower saturation capacities at high bonding densities, the combinations leading to similar apparent retention in RPLC. The use of a low surface coverage column is recommended for preparative purposes.

Influence of solvents on the conformational order of C18 alkyl modified silica gels

The influence of solvents on the conformational order of C(18) alkyl modified silica gels is studied by means of variable temperature FTIR spectroscopy. Symmetric and anti-symmetric CH(2) stretching modes were utilized for getting qualitative information about the changes in alkyl chain conformational order as a function of both protonated and perdeuterated solvents. It was found that interaction between the C(18) alkyl modified silica gels and the mobile phase results in pronounced changes of the alkyl chain conformational order. Furthermore, it was observed that some perdeuterated solvents exhibit isotope effects, which again is reflected by a different alkyl chain conformational order as compared to the corresponding stationary phases with protonated solvents.

Molecular Dynamics Simulations of Alkylsilane Stationary-Phase Order and Disorder. 1. Effects of Surface Coverage and Bonding Chemistry

"Shape-selective" polymeric alkylsilane stationary phases are routinely employed over the more common monomeric phases in reversed-phase liquid chromatography (RPLC) to improve the separation of geometric isomers of shape-constrained solutes. We have investigated the molecular dynamics of chromatographic models that represent both monomeric and polymeric stationary phases with alkylsilane surface coverages and bonding chemistries typical of actual materials in an effort to elucidate the molecular-level structural features that control shape-selective separations. The structural characterization of these models is consistent with previous experimental observations of alkyl chain order and disorder: (1) alkyl chain order increases with increased surface coverage; and (2) monomeric and polymeric phases with similar surface coverages yield similar alkyl chain order (although subtle differences exist). In addition, a significant portion of the alkyl chain proximal to the silica surface is disordered (primarily gauche conformations) and the distal end is most ordered. Models that represent shape-selective RPLC phases possess a significant region of distal end chain order with primarily trans dihedral angle conformations. This is consistent with the view that the alkyl chains comprising polymeric stationary phases contain a series of well-defined and rigid voids in which shape-constrained solutes can penetrate and hence be selectively retained.

Temperature Dependence of Hydrogen Bonding: An Investigation of the Retention of Primary and Secondary Alcohols in Gas−Liquid Chromatography

Configurational-bias Monte Carlo simulations in the Gibbs Ensemble were carried out to investigate the analyte partitioning of n-pentane, n-hexane, n-heptane, 1-propanol, and 2-propanol into a dioctyl ether retentive (stationary) phase used in gas-liquid chromatography. The united-atom version of the TraPPE (transferable potentials for phase equilibria) force field was used to model all analytes and the solvent. The analyte partition coefficients, Gibbs free energies of transfer, and Kovats retention indexes were calculated at four different temperatures ranging from 303.15 to 348.15 K. Although hydrogen bonding is a major contributor to the retention of the alcohol analytes over the entire temperature range, its importance for the separation factor between the primary and secondary alcohol decreases substantially with increasing temperature. The enthalpies and entropies for hydrogen bond formation were also estimated from the temperature dependence of the corresponding equilibrium constants. In agreement with experimental measurements, it is observed that the hydrogen bond involving 1-propanol is enthalpically favored, but entropically disfavored compared to 2-propanol.