How Ternary Mobile Phases Allow Tuning of Analyte Retention in Hydrophilic Interaction Liquid Chromatography

Inception and development of a LC-MS/MS assay for the multiplexed quantitation of nine human drug transporter biomarkers

Background: It has become common practice to assess solute carrier transporter (SLC)-mediated drug-drug interactions (DDIs) by quantitating various individual endogenous compounds as biomarkers in human plasma and urine. The goal of this work was to develop biomarker multiplex assays that could be utilized during first in human studies to support the simultaneous assessment of clinical DDI risk across various SLCs. Methodology: Hydrophilic interaction chromatography-MS/MS methods were developed, and validations were performed. Results: The multiplex assays were applied to a first in human study. Placebo/reference subject biomarker data were consistent with single assay in-house and published data. Conclusion: This work demonstrates the utility of these multiplex methods to support the concurrent evaluation of clinical DDI risk across various SLCs.

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.

Retention and mass transfer properties of the series of unbonded, amide-bonded, and alkylsulfobetaine-bonded ethylene bridged hybrid hydrophilic interaction liquid chromatography columns

This work investigates the link between the retentivity and the stationary phase to mobile phase mass transfer resistance of hydrophilic interaction liquid chromatography (HILIC) columns packed with the same base ethylene-bridged hybrid particles (BEH). The retention volumes, the plate heights, and the volume of the adsorbed water layer were measured for the ACQUITY^TM UPLC^TM BEH^TM 130 Å HILIC Column (unbonded BEH), ACQUITY UPLC BEH 130 Å Amide Column (amide group attached), and Atlantis^TM Premier BEH 95 Å Z-HILIC (zwitterionic group attached) Column. The method of Guo (toluene retention volumes in pure acetonitrile and in the HILIC eluent) was validated from the UNIFAC group-contribution method and applied to measure accurately the water layer volumes in these columns. A strong correlation was found between the retention volumes of most neutral polar analytes and the volume of the water layer adsorbed in the HILIC column. The fraction of the pore volume occupied by the water layer increases significantly from the BEH HILIC Column to the BEH Amide Column, and to the BEH Z-HILIC Column. This is explained by the water solvation of the attached ligands in the pore volume of the BEH Particles and to the smaller average mesopore size of the BEH Z-HILIC Particles. A second and strong correlation is also observed between the water content in the HILIC particle and the stationary phase to mobile phase mass transfer resistance of the HILIC columns at high mobile phase linear velocities. The measured intra-particle diffusivity normalized to the bulk diffusion coefficient decreased from 0.33 (BEH HILIC Column) to 0.10 (BEH Amide Column) and to only 0.03 (BEH Z-HILIC Column) for comparable retention of cytosine. These results are fully consistent with the higher viscosity of the internal eluent (higher water content) and higher internal obstruction for diffusion (smaller mesopores and internal porosity) in the BEH Z-HILIC Particles. Still, in gradient elution mode, the peak capacity was found to be 18% higher for the BEH Z-HILIC Column than that on the BEH Amide Column because the retention factors at elution were smaller when maintaining the same analysis time and starting eluent composition.

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.

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.

Collective Solvation and Transport at Tetrahydrofuran-Silica Interfaces for Separation of Aromatic Compounds: Insight from Molecular Dynamics Simulations

We have performed umbrella sampling molecular dynamics simulations to study the separation mechanism of aromatic compounds at the tetrahydrofuran (THF)-methanol-silica interface by liquid chromatography. Solute molecules with different polarities (naphthol and naphthalene) are selected as representative aromatic compounds. For the polar solute (naphthol), the free energy profile shows a deep minimum near the THF-silica interface, suggesting strong interactions with the polar surface. When methanol is added to the interface, there is a sharp increase in naphthol's free energy minimum, and the corresponding diffusion dynamics also undergoes a dramatic change. These findings explain the fast separation mechanism in recent experiments of separating fused ring compounds in asphaltenes with liquid chromatography. Further solvation structure and orientation analysis suggest that apolar and polar solutes may find their own comfort zones several angstroms away from the interface, and their phenyl ring's orientations would undergo a parallel-to-perpendicular transition as the solute molecule moves away from the surface. Extending our simulation studies to systems with different solute concentrations reveals that there is a decrease in the adsorption free energy accompanied by enhanced surface diffusion as the solute concentration increases, which is related to the crowding in the interfacial layers. Our simulation analysis gives a detailed microscopic description of solute solvation and transport at the THF-silica chromatography interface and will be helpful for improving separation protocols in future applications.

Insights from molecular simulations about dead time markers in reversed-phase liquid chromatography

Among the most popular compounds to estimate the hold-up time in reversed-phase liquid chromatography (RPLC) are acetone and uracil, which are considered as too small and too polar, respectively, for retention by the hydrophobic stationary phase, although their observed elution behavior does not fully support this assumption. We investigate how acetone and uracil as solutes interact with the chromatographic interface through molecular dynamics simulations in an RPLC mesopore model of a silica-supported, endcapped, C₁₈ phase equilibrated with a water (W)‒acetonitrile (ACN) mobile phase. The simulation results provide a molecular-level explanation for the observed elution behavior of acetone and uracil, but also question whether true dead time markers for RPLC exist. Both solutes have a density maximum in the interfacial region in addition to a low presence in the bonded-phase region, but these density peaks clearly differ from the adsorption and partitioning peaks of true analytes. Acetone partially behaves like a co-solvent of ACN and partially like the analyte acetophenone. Like ACN, acetone can be found in the first and second layer of solvent molecules at the silica surface; like acetophenone, acetone adsorbs to the bonded-phase chains by orienting its polar group to the bulk region to sustain hydrogen bonds with W molecules. Uracil behavior is governed by a need for extensive hydrogen-bond coordination by W molecules. Uracil adsorbs to the very edge of the bonded-phase chains, on the bulk-region side of the ACN density maximum in the interfacial region. Further penetration into the chains is prevented by the absence of W molecules, which are not found deeper in the bonded phase, except at the silica surface. Contrary to true analytes, accumulation of uracil and acetone in the interfacial region ceases at an equimolar presence of W and ACN in the mobile phase (at 70‒80% ACN volume fraction). Uracil achieves a closer approximation of the stationary-phase limit than acetone, but carries the risk of HILIC retention at high ACN fraction in the mobile phase.

Recent progress in the fundamental understanding of hydrophilic interaction chromatography (HILIC)

With the exponential growth in the application of the HILIC technique, there has been a significant progress in understanding the fundamental aspects of hydrophilic interaction chromatography.