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

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.

From sound check to encore: A journey through high-resolution mass spectrometry-based food analyses and metabolomics

This manuscript presents a comprehensive review of high-resolution mass spectrometry in the field of food analysis and metabolomics. We have followed the historical evolution of metabolomics, its associated techniques and technologies, and its increasing role in food science and research. The review provides a critical comparison and synthesis of tentative identification guidelines proposed for over 15 years, offering a condensed resource for researchers in the field. We have also examined a wide range of recent metabolomics studies, showcasing various methodologies and highlighting key findings as a testimony of the versatility of the field and the possibilities it offers. In doing so, we have also carefully provided a compilation of the software tools that may be employed in this type of studies. The manuscript also explores the prospects of high-resolution mass spectrometry and metabolomics in food science. By covering the history, guidelines, applications, and tools of metabolomics, this review attempts to become a comprehensive guide for researchers in a rapidly evolving field.

Exploring the application limits of different hold-up time markers in supercritical fluid chromatography

The study focuses on the application range of nitrous oxide as a hold-up time marker in supercritical fluid chromatography (SFC). This compound has been suggested a decade ago to be used as unretained marker, something that the field of SFC was missing for a long time, since its beneficial properties make it an ideal candidate as hold-up time marker. Determination of the hold-up volume and actual volumetric flow rates have always been problematic in SFC due to the compressibility of carbon dioxide and one part of this is the difficulty of hold-up time measurements. Depending on the mobile phase, different methods have been used to measure the hold-up time with varying results. Nitrous oxide and other molecules have been compared in different conditions, mobile phases and stationary phases. In all cases, nitrous oxide gave the lowest elution times. However, detection was difficult in mobile phases containing 10% or more of organic modifier, because most solvents mask the signal of nitrous oxide. Interestingly, the choice of stationary phase also had a slight effect on detection, while different pressure and temperature settings affected each compound in a different manner.

BART: A transferable liquid chromatography retention time library for bile acids

Identification of unknown bile acids, especially the distinguishment between isomers, requires retention times of a large number of reference standards, which are often not commercially available. Meanwhile, published retention information cannot be directly transferred across labs due to the differences between liquid chromatography (LC) systems, such as different extra column volume and dwell volume. To improve this situation, a transferrable retention time library for bile acids named BART was developed. BART was composed of isocratic retention models of 272 bile acids and a software tool to predict their gradient retention times on various LC systems. The isocratic retention times of bile acids were acquired on a Waters BEH C18 column with mobile phases of acidic ammonium acetate buffer and acetonitrile, and fit to the quadratic solvent strength model (QSSM). Segmented linear gradient retention times were calculated with holdup time (t0), dwell time (tD) and actual gradient profile corrected using 21 bile acid calibration standards. In addition to the reference system where the isocratic retention times were acquired, this approach has been validated on four other LC-MS systems in four labs with two gradient methods. Average root mean square errors (RMSE) between predicted and experimental retention times were 0.052 and 0.054 min for the two gradients tested, which were 9-fold more accurate than referring to a static retention time library. The library is freely available at https://bafinder.github.io/.

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.

Towards an accurate method for column void volume determination using liquid chromatography-mass spectrometry

Prediction of analyte retention times requires prior knowledge of the column void volume, the measurement of which is still highly contested within the literature and therefore experimental based prediction is often used. In this study, we investigated deuterated acetonitrile as an isotopically labelled mobile phase component to observe its elution behaviour in a binary mixture with water at 25 different mobile phase compositions (from 5 to 95 vol.% of acetonitrile), on two stationary phases (C8 and C18), and at two temperatures (30 and 40 °C) using LC-MS. The same experimental design was additionally used for three commonly used neutral void volume markers: uracil, phloroglucinol and N,N-dimethylformamide. Temperature was observed to influence the elution of acetonitrile in an inversely proportional manner with higher temperatures coinciding with lower elution times. By utilizing a three-way ANOVA, the composition of the mobile phase has been shown to have a significant effect on deuterated acetonitrile and other investigated void volume markers, demonstrating the fact that both void volume markers and acetonitrile itself exhibit retention-like behaviour. Excess adsorption isotherms for acetonitrile were calculated using deuterated acetonitrile elution data. The comparison of void volumes, obtained with conventional neutral void volume markers, revealed the former to be 24-36% lower than the void volume obtained using deuterated acetonitrile, as an isotopically labelled mobile phase component. For a water:acetonitrile mobile phase, the minor disturbance method using deuterated acetonitrile to obtain an integral average void volume (2.08 and 2.05 mL for C18 at 30 and 40 °C, respectively and 2.16 and 2.13 mL for C8 at 30 and 40 °C, respectively) was found to be the most appropriate method for determining the elusive column void volume.

Harmonization of Experimental Methods Used to Measure the True Hold-Up Volume of Liquid Chromatography Columns

There are as many measurement values of the true column hold-up volume, V0, as techniques applied to evaluate this most important property in liquid chromatography (LC). The relative errors made on V0 measurements using conventional “non-retained” markers—such as acetone, uracil, or thiourea in reversed-phase liquid chromatography (RPLC), or benzene or acenaphthene in hydrophilic interaction chromatography (HILIC)—can be as large as ±30%. This situation is extremely confusing for LC users who wish to classify and predict the retention behavior of LC columns. In this work, along with advances in mass spectrometry (MS) instrumentation, the hold-up volume of any LC column—including, but not limited to, RPLC, HILIC, ion exchange chromatography (IEX), and mixed-mode columns—is accurately measured by injecting labeled deuterated acetonitrile (CD3CN) molecules and detecting them selectively by MS-single ion reaction (m/z = 45) using non-labeled and pure acetonitrile (CH3CN) as the eluent. This proposed harmonization of all conventional V0 measurement methods is illustrated and successfully applied to RPLC, HILIC, anion exchange (AEX), and RP-AEX mixed-mode chromatography, irrespective of the mobile phase composition selected.

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.

Confinement Effects on Distribution and Transport of Neutral Solutes in a Small Hydrophobic Nanopore

Molecular dynamics simulations are used to study confinement effects in small cylindrical silica pores with extended hydrophobic surface functionalization as realized, for example, in reversed-phase liquid chromatography (RPLC) columns. In particular, we use a 6 nm cylindrical and a 10 nm slit pore bearing the same C₁₈ stationary phase to compare the conditions inside the smaller-than-average pores within an RPLC column to column-averaged properties. Two small, neutral, apolar to moderately polar solutes are used to assess the consequences of spatial confinement for typical RPLC analytes with water (W)-acetonitrile (ACN) mobile phases at W/ACN ratios between 70/30 and 10/90 (v/v). The simulated data show that true bulk liquid behavior, as observed over an extended center region in the 10 nm slit pore, is not recovered within the 6 nm cylindrical pore. Instead, the ACN-enriched solvent layer around the C₁₈ chain ends (the ACN ditch), a general feature of hydrophobic interfaces equilibrated with aqueous-organic liquids, extends over the entire pore lumen of the small cylindrical pore. This renders the entire pore a highly hydrophobic environment, where, contrary to column-averaged behavior, neither the local nor the pore-averaged sorption and diffusion of analytes scales directly with the W/ACN ratio of the mobile phase. Additionally, the solute polarity-related discrimination between analytes is enhanced. The consequences of local ACN ditch overlap in RPLC columns are reminiscent of ion transport in porous media with charged surfaces, where electrical double-layer overlap occurring locally in smaller pores leads to discrimination between co- and counterionic species.

Retention modeling and adsorption mechanisms in reversed-phase liquid chromatography

To interpret the dependence of solute retention behavior on modifier content in reversed-phase liquid chromatography, a theoretical framework, based on the concentration dependence of solvophobic forces imposed on solutes and the competitive adsorptions of solutes and solvent modifiers, was proposed. The generality of the developed model was demonstrated by comparing the model with conventional retention models. The linear dependence of the Gibbs energy change of solute adsorption with respect to the modifier concentration was assumed, and the model was fitted to the experimental results, with good agreement demonstrated between the experimental data and the model. Retention behaviors were inferred to be determined by two key dimensionless groups that represented the reductions in the retention factors resulting from a weakened solvophobic interaction and modifier competitive adsorption. The retention behaviors were successfully deconvoluted for each contribution as a function of the modifier concentration by using the fitted parameters. The effects of both contributions on the retention behaviors were enhanced for the solutes with aromatic groups. The standard Gibbs energy change _SL^o of benzene adsorption was found to depend linearly on the number of modifier molecules present but independent of modifier identity. For the solutes associated with hydrogen-bonding groups, the degree of reduction in the solvophobic interactions was considerably reduced. Hence, the relative contributions of both mechanisms to solute retention depend greatly on the solute structure. Perturbation method was performed to investigate the modifier adsorption mechanisms. The results show that the standard Gibbs energy change _SL^o for the first-layer adsorption of modifiers changed linearly with the carbon number of modifier molecule. These results demonstrated that the proposed model can offer a physically consistent quantitative description of retention when solvent composition is varied.

Brownian sieving enhancement of microcapillary hydrodynamic chromatography. Analysis of the separation performance based on Brenner's macro-transport theory

In a recent article [Analytical Chemistry, 93(17), 6808-6816 (2021)], an unconventional device configuration enforcing a Brownian sieving mechanism was proposed as proof of concept for the efficient implementation of microcapillary hydrodynamic chromatography (MHDC). In this article, we perform a thorough analysis of the device geometry and of operating conditions, in order to single out the optimal configuration maximizing separation resolution. Brenner's macro-transport theory provides the technical picklock to perform the search for the optimum over a wide choice of device geometries and for a range of values of the particle Péclet number covering most conditions encountered in practical implementations of MHDC. Specifically, effective transport coefficients defining the dynamics of the suspended phase are obtained by the solution of a two-dimensional steady-state advection-diffusion equation defined onto the channel cross-section. The eigenvalue/eigenfunction structure of the associated transient problem is exploited in order to quantify the timescale for reaching the macro-transport regime conditions. Based on this timescale and on the effective transport parameters, an estimate of the column length necessary to achieve a prescribed level of separation resolution is obtained. We identify device geometry and operating conditions where the capillary length is shrunk down by a factor above ten compared to the standard MHDC configuration. Lagrangian stochastic statistics of particle ensembles are used to validate the results obtained through Brenner's macro-transport approach. The method proposed can be readily generalized to other classes of device geometries enforcing the same Brownian sieving mechanism.

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.

Perspective on the Future Approaches to Predict Retention in Liquid Chromatography

The demand for rapid column screening, computer-assisted method development and method transfer, and unambiguous compound identification by LC/MS analyses has pushed analysts to adopt experimental protocols and software for the accurate prediction of the retention time in liquid chromatography (LC). This Perspective discusses the classical approaches used to predict retention times in LC over the last three decades and proposes future requirements to increase their accuracy. First, inverse methods for retention prediction are essentially applied during screening and gradient method optimization: a minimum number of experiments or design of experiments (DoE) is run to train and calibrate a model (either purely statistical or based on the principles and fundamentals of liquid chromatography) by a mere fitting process. They do not require the accurate knowledge of the true column hold-up volume V₀, system dwell volume V_dwell (in gradient elution), and the retention behavior (k versus the content of strong solvent φ, temperature T, pH, and ionic strength I) of the analytes. Their relative accuracy is often excellent below a few percent. Statistical methods are expected to be the most attractive to handle very complex retention behavior such as in mixed-mode chromatography (MMC). Fundamentally correct retention models accounting for the simultaneous impact of φ, I, pH, and T in MMC are needed for method development based on chromatography principles. Second, direct methods for retention prediction are ideally suited for accurate method transfer from one column/system configuration to another: these quality by design (QbD) methods are based on the fundamentals and principles of solid-liquid adsorption and gradient chromatography. No model calibration is necessary; however, they require universal conventions for the accurate determination of true retention factors (for 1 < k < 30) as a function of the experimental variables (φ, T, pH, and I) and of the true column/system parameters (V₀, V_dwell, dispersion volume, σ, and relaxation volume, τ, of the programmed gradient profile at the column inlet and gradient distortion at the column outlet). Finally, when the molecular structure of the analytes is either known or assumed, retention prediction has essentially been made on the basis of statistical approaches such as the linear solvation energy relationships (LSERs) and the quantitative structure retention relationships (QSRRs): their ability to accurately predict the retention remains limited within 10-30%. They have been combined with molecular similarity approaches (where the retention model is calibrated with compounds having structures similar to that of the targeted analytes) and artificial intelligence algorithms to further improve their accuracy below 10%. In this Perspective, it is proposed to adopt a more rigorous and fundamental approach by considering the very details of the solid-liquid adsorption process: Monte Carlo (MC) or molecular dynamics (MD) simulations are promising tools to explain and interpret retention data that are too complex to be described by either empirical or statistical retention models.