Retention mechanism in combined hydrodynamic and slalom chromatography for analyzing large nucleic acid biopolymers relevant to cell and gene therapies.

Slalom chromatography (SC) was discovered in 1988 for analyzing double-stranded (ds) DNA. However, its progress was impeded by practical issues such as low-purity particles, sample loss, and lack of a clear retention mechanism. With the rise of cell and gene therapies and the availability today of bio-inert ultra-high-pressure liquid chromatography (UHPLC) columns and systems, SC has regained interest. In SC, the elution order is opposite to that observed in hydrodynamic chromatography (HDC): larger DNA molecules are more retained than small ones. Yet, the underlying SC retention mechanism remains elusive. We provide the physicochemical background necessary to explain, at a microscopic scale, the full transition from a HDC to a SC retention mechanism. This includes the persistence length of the DNA macromolecule (representing DNA stiffness), their relaxation time (τR) from the non-equilibrium contour length to the equilibrium entropic configuration, and the relationship between the mobile phase shear rate (〈γ̇〉) in packed columns and the DNA extended length. We propose a relevant retention model to account for the simultaneous impact of hydrodynamic chromatography (HDC) and SC on the retention factors of a series of large and linear dsDNAs (ranging from 2 to 48 kbp). SC data were acquired using bio-inert MaxPeakTM Columns packed with 1.7µm BEHTM 45 Å, 1.8µm BEH 125 Å, 2.4µm BEH 125 Å, 5.3µm BEH 125 Å, and 11.3µm BEH 125 Å Particles, an ACQUITYTM UPLCTM I-class PLUS System, and either 1 × PBS (pH 7.4) or 100 mM phosphate buffer (pH 8) as the mobile phase. SC is a non-equilibrium retention mode that is dominant when the Weissenberg number (Wi=〈γ̇〉τR) is much larger than 10 and the average extended length of DNA exceeds the particle diameter. HDC, on the other hand, is an equilibrium retention mode that dominates when Wi<1 (DNA chains remaining in their non-extended configuration). Maximum dsDNA resolution is observed in a mixed HDC-SC retention mode when the extended length of the DNA is approximately half the particle diameter. This work facilitates the development of methods for characterizing various plasmid DNA mixtures, containing linear, supercoiled, and relaxed circular dsDNAs which all have different degree of molecular stiffness.

Preparation and investigation of two-component silica-modified hydrophobic layer for minimizing retention loss of reversed-phase chromatographic columns using fully aqueous mobile phases.

Chromatographers often employ fully aqueous mobile phases to retain highly polar compounds in reversed-phase liquid chromatography (RPLC). However, when the flow rate is interrupted, either accidentally or intentionally, a substantial loss in retention occurs due to the spontaneous dewetting of water from the hydrophobic surface of conventional RPLC-C18 particles. Previous studies have shown that maintaining a low C18 surface coverage (approximately 1.5 µmol/m2) can mitigate water dewetting by increasing chain disorder, facilitating the intercalation of water clusters between the C18-bonded chains, and keeping the mesopores wetted. In this research, we explore the potential and additional benefits of using two-component surface bonding materials (C8/C18 and PhenylHexyl (PhHx)/C18) at a constant and low total surface coverage of 1.51 ± 0.15 µmol/m2. We synthesized seven one- and two-component modified silica particles with a volume average particle size of 5.22 µm and an average mesopore size of 104 Å. The surface coverage was increased from 0 to 0.54, 1.00, and to 1.66 µmol2 for C8 chains and from 0 to 0.52, 0.70, and to 1.65 µmol2 for PhHx ligands. To prevent interactions between water and any unreacted silanols, all seven derivatized particles were heavily endcapped with trimethylsilane (TMS) reagent. The fraction of the surface area remaining in contact with water was determined by measuring the retention times of weakly (thiourea) and strongly (thymine) retained compounds at intervals of 1, 2, 4, 8, 16, 32, and 64 minutes following the cessation of flow. Two distinct column temperatures, 24°C and 60°C, were employed in the experiments. Retention losses were found to be minimized in the presence of a small quantity of C8 chains (less than 40% of the total surface coverage). Additionally, it is essential to consider substantial fractions of PhHx chains, as long as the presence of the PhHx ligand does not significantly impact retention and selectivity. Combining mixed RPLC bondings with a low total surface coverage of approximately 1.5 µmol/m2 emerges as a viable solution for further minimizing retention loss in standard C18-bonded RPLC columns, particularly within the surface coverage range of 2.5-3.0 µmol/m2.

Physical origin of the peak tailing of monoclonal antibodies in size-exclusion chromatography using bio-compatible systems and columns.

The analysis of mixtures containing monoclonal antibody (mAb) (approximately 150 kDa molecular weight) and sub-unit impurities (approximately 100 kDa) is challenging, even when adopting the latest ultra-high-pressure liquid chromatography (UHPLC) columns (4.6 mm [Formula: see text] 150 mm coated hardware, 1.7 [Formula: see text]m 250 BEH[Formula: see text] Surface-modified Particles) and systems (ACQUITY[Formula: see text] UPLC[Formula: see text] I-class Bio Plus). The main issue still encountered is a persistent tail of the mAb peak. Here, the physical origin(s) of such peak tailing in size-exclusion chromatography (SEC) are investigated from both fundamental and practical approaches. Up to five relevant physical origins are analyzed: sample heterogeneity (glycoforms), UHPLC system dispersion, strong residual binding of the mAb to the SEC particles (via hydrophobic and/or electrostatic interactions) and to the stainless steel column/system hardware, slow escape kinetics of the mAb from the SEC particles, and flow heterogeneity caused by the non-ideal slurry packing of SEC columns. Experiments (testing sample heterogeneity, system dispersion, and strong residual interactions) and calculations (predicting the transient absorption/escape kinetics in a single SEC particle and the two-dimensional peak concentration profiles) altogether unambiguously demonstrate that the observed mAb peak tailing is caused primarily by the long-range velocity biases across the SEC column combined with the slow transverse dispersion of mAbs. Therefore, improvement in the resolution between mAb and sub-unit fragment impurities can only be achieved by increasing the column length, e.g., by applying recycling chromatography at acceptable pressures.

Resolution limits of size exclusion chromatography columns identified from flow reversal and overcome by recycling liquid chromatography to improve the characterization of manufactured monoclonal antibodies.

The flow reversal (FR) technique consists of reversing the flow direction along a chromatographic column. It is used to reveal the origin (such as poor column packing, active sites, or slow absorption/escape kinetics) for the resolution limit of 4.6 mm × 150 mm long columns packed with 1.7 µm 200 Å Bridge-Ethylene-Hybrid (BEHTM) Particles. These columns are used to separate manufactured monoclonal antibodies (mAb, ∼ 150 kDa) from their close impurities (or IdeS fragments, ∼ 100 kDa) by size exclusion chromatography (SEC). FR unambiguously demonstrates that the resolution limit of these SEC columns is primarily due to long-range flow velocity biases covering distances of at least 500 µm across the column diameter. This confirms the existence of center-to-wall flow heterogeneities which cause undesirable tailing for the mAb peak. Because the transverse dispersion coefficient (Dt=1.1 × 10-6 cm2/s) of mAbs across the column diameter is intrinsically low, the bandspreading of the mAb in a single flow direction is in part reversible upon reversing the flow direction. For the very same residence time in the column, the column efficiency is found to increase by +85% relative to that observed under conventional elution mode. The observed peak tailing of the mAb and its sub-units is not caused by active surface sites or by slow absorption/escape from the BEH Particles. Therefore, the most critical mAb impurities (hydrolytic degradation Fab/c and IdeS [Formula: see text] fragments) can only be successfully separated and quantified with acceptable accuracy by adopting alternate pumping recycling liquid chromatography (APRLC). APRLC enables the full baseline separation of the mAb and 100 kDa mAb fragments and partial separation of Fab/c and [Formula: see text] fragments after increasing the number of cycles to ten. It was made possible to accurately measure the relative abundances of the mAb (99.0 ± 0.1%), [Formula: see text] fragment (0.88 ± 0.03%), and Fab/c immunogenic fragment (0.13 ± 0.02%) in less than 45 min for a total mAb sample load of only 5 µg. Still, further improvements are needed to increase the sensitivity of the APRLC method and to reduce the solvent consumption by adopting narrow-bore 2.1 mm i.d. SEC columns.

Absorption and escape kinetics of spherical biomolecules from fully porous particles utilized in size exclusion chromatography.

The increasing demand for the characterization of large biomolecules such as monoclonal antibodies, double-stranded deoxyribonucleic acid (dsDNA), and virus-like particles (VLPs) is raising fundamental questions pertaining to their absorption (ingress) and escape (egress) kinetics from fully porous particles. The exact expression of their concentration profiles is derived as a function of time and radial position across a single sub-3 µm Bridge-Ethylene-Hybrid (BEHTM) Particle present in size exclusion chromatography (SEC) columns. The boundary condition at the external surface area of the particle is a rectangular concentration profile mimicking the passage of the chromatographic zone. Four different BEH Particles were considered in the calculations depending on the molecular size of the analyte: 2.0 µm 100 Å BEH Particles for small molecules, 2.0 µm 200 Å BEH Particles for monoclonal antibodies, 2.0 µm 300 Å BEH Particles for dsDNA (100 base pairs), and 2.5 µm 900 Å BEH Particles for virus-like particles (VLPs). The calculated concentration profiles of small molecules and monoclonal antibodies confirm that all BEH Particles present in the column reach quasi-instantaneously thermodynamic equilibrium with the bulk mobile phase during the passage of the chromatographic band. This is no longer the case for larger biomolecules such as dsDNA or VLPs, especially when the SEC particle is located near the column inlet and for high velocities. The kinetics of biomolecule egress is slower than its kinetics of ingress leading to pronounced peak tailing. The mean concentration of the largest biomolecules in the SEC particles remains always smaller than the maximum bulk concentration. This persistent and transient intra-particle diffusion regime has direct implications on the theoretical expressions of the observed retention factors and plate heights. Classical theories of chromatography assume uniform spatial distribution of the analyte in the particle volume: this hypothesis is not verified for the largest biomolecules. These results imply that non-porous particles or monolithic structures are the most promising stationary phases for the separation and purification of the largest biomolecules in life science.

Multiple-open-tubular column enabling transverse diffusion. Part 3: Simulation of solute dispersion along a real three dimensional-printed column with quadratic channels.

Multiple-open-tubular columns enabling transverse diffusion (MOTTD) consist of straight and parallel flow-through channels separated by a mesoporous stationary phase. In Part 1, a stochastic model of band broadening along MOTTD columns accounting for longitudinal diffusion, trans-channel velocity bias, and mass transfer resistance in the stationary phase was derived to demonstrate the intrinsic advantage of MOTTD columns over classical particulate columns. In Part 2, the model was refined for the critical contribution of the channel-to-channel polydispersity and applied to address the best trade-off between analysis speed and performance. In this Part 3, a MOTTD column with a square array of quadratic channels is fabricated by 3D-printing (combining polymer stereolithography with photolithography using photomasks) to deliver unprecedently small apparent channel diameters of 117.6 ± 5.0 µm. The colors in the microscopy photographs of the actual 3D-printed channels are binarized to delimitate the mobile phase volume from the stationary phase volume. The same numerical simulations as those in Part 2 are then performed for two MOTTD columns (external porosity ϵe=31.7%, same apparent channel diameter 117.6 µm): one containing 16 virtual perfect quadratic channels and the other 16 real 3D-printed channels. The reduced velocities (or Peclet numbers) are varied over a wide range from 0.2 to 5000 and the zone retention factors were fixed at k1=1.04, 5, and 25. The results demonstrate that smoothing the edges of the targeted quadratic channels by the 3D-printed technique is advantageous in terms of solute dispersion. It outperforms the negative effect of the channel-to-channel polydispersity which is mitigated by transverse diffusion of the analyte in the stationary phase. For Peclet numbers larger than 50, the HETP of the 3D-printed MOTTD column is found 7%, 15%, and 16% smaller than that of the MOTTD column consisting of a square array of perfect quadratic channels. This confirms the known effect of channel geometry on solute dispersion in microfluidic systems. Flow channels in fabricated MOTTD columns are preferred to be circular so that the distribution of transverse diffusion lengths across the open channels remains as tight as possible. Finally, the general theory of nonuniform columns of Giddings reveals that the polydispersity of the cross-sectional area (RSD 8.4%) along a single 3D-printed channel has no negative impact on solute dispersion in MOTTD columns. Overall, MOTTD columns could become a serious alternative technology to conventional particulate columns. This implies a novel fabrication process that delivers circular channel diameters smaller than 10 µm, cross-sectional area polydispersity no larger than 25%, external porosities in a range from 15% (high speed separations) to 75% (high performance separations), and conventional mesoporous silica as the stationary phase. It adresses new synthesis routes based on either organic fibers or tubular micelle templating agents in suspension with silica gel solutions.

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 ACQUITYTM UPLCTM BEHTM 130 Å HILIC Column (unbonded BEH), ACQUITY UPLC BEH 130 Å Amide Column (amide group attached), and AtlantisTM 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.

Prediction of surface excess adsorption and retention factors in reversed-phase liquid chromatography from molecular dynamics simulations.

An alternative method to the classical fit of semi-empirical, statistical, or artificial intelligence-based models to retention data is proposed to predict surface excess adsorption and retention factors in liquid chromatography. The approach is based on a fundamental, microscopic description of the liquid-to-solid adsorption of analytes taking place at the interface between a bulk liquid phase and a solid surface. Molecular dynamics (MD) simulations are performed at T=300 K in a 100 Å wide slit-pore model (ß-cristobalite-C18 surface in contact with an acetonitrile/water mobile phase) to quantify a priori the retention factors of small molecules expected in reversed phase liquid chromatography (RPLC). Uracil is chosen as the reference "non-retained" marker, whereas benzyl alcohol, acetophenone, benzene, and ethylbenzene are four selected retained, neutral compounds. The MD simulations allow to determine the pore-level density profiles of these five compounds, i.e., the variation of the analyte concentration as a function of distance from the silica surface. The retention factors of the retained analytes are expressed using their respective calculated surface excess adsorption relative to uracil. By definition, the retention factors are proportional to the surface excess adsorbed and the proportionality constant is directly scaled to the retention time of the "non-retained" marker. Experimentally, a 4.6 mm × 150 mm RPLC-C18 column packed with 5 µm 100 Å High Strength Silica (HSS)-C18 particles is used and the retention times of these five compounds are measured. The volume fraction of acetonitrile in water increases from 20 to 90% generating a wide range of retention factors from 0.15 to 183 at T=300 K. The results demonstrate very good agreement between the MD-predicted surface excess adsorption data and measured retention factors (R2> 0.985). A systematic error is observed as the proportionality constant is not exactly scaled to the retention time of uracil. This is most likely caused by the differences between the chemical and morphological features of the slit-pore model adopted in the MD simulations and those of the actual HSS-C18 particles: the average surface coverage with C18 chains, the geometry of the mesopores, and the pore size distribution. Specifically, the impact on RPLC retention of slight, local variations in surface chemistry (e.g., functional group density and uniformity) and how this aspect is affected by the pore space morphology (e.g., pore curvature and size) is worth investigating by future MD simulations.

Modeling of the transient diffusion regime in fully porous particles-Application to the analysis of large biomolecules by ultra-high pressure liquid chromatography.

The transient diffusion regime of large biomolecules such as monoclonal antibodies, double-stranded (ds) DNA (base pair number ∼ 100), and virus-like particles is modeled in a single fully porous particle utilized as size exclusion chromatography (SEC) packing materials in ultra-high pressure liquid chromatography (UHPLC). The expression of the time and space dependent concentration profiles is derived for a step concentration change. Four different UHPLC particles were considered in the calculations depending on the size of the analyte : 2.0 µm 125 Å Bridge-Ethylene-Hybrid (BEH) XBridgeTM particles for small molecules (size: 4 Å), 2.0 µm 125 Å and 250 Å BEH particles for monoclonal antibodies (size: 55 Å), 2.0 µm 250 Å and 2.5 µm 450 Å BEH particles for dsDNA (size: 160 Å), and 2.5 µm 900 Å BEH particles for virus-like particles (size: 500 Å). The accessible porosities and the hindrance diffusion factors of these analytes were determined from the physical reconstruction of the internal structure of a 2.0 µm 130 Å BEH particle and from the simulation of the analyte mobility in the reconstructed mesopores. The analysis of the transient diffusion profiles reveals that there is insufficient time for the largest analytes to fully equilibrate the BEH particles in UHPLC. The dynamic capacity of the BEH particles is estimated to decrease from more than 99% for small molecules to 96% for monoclonal antibodies, 74% for 100 base pair dsDNA, and to less than 2% for virus-like particles. A full but reduced column capacity relative to fully porous particles can be recovered for monoclonal antibodies by considering superficially porous particles with a core-to-particle diameter ratio of about 0.8. In contrast, this is not possible for either large dsDNA (base pair number ≥ 100) or virus-like particles because the shell thickness would become smaller than the mesopore size required. The presented results open up two research avenues in order to analyze such large biomolecules by UHPLC: prepare ultra-high pressure-resistant sub-3 µm silica-based particles with pore sizes much larger than 1000 Å and design proper nonporous macrostructures to load, trap, separate, and elute rapidly by gradient UHPLC.

On the road toward highly efficient and large volume three-dimensional-printed liquid chromatography columns?

The current performance of commercially packed liquid chromatography columns is limited by the random structure of the packed bed and by the wall-to-center heterogeneity of its structure. The minimum reduced plate heights observed are not smaller than 1.4, whereas they could theoretically be as low as 0.1 for dense and perfectly ordered packings of spheres. To bridge this gap, a wide inner diameter column with an ordered macroporous structure is printed in three dimensions by stereolithography of poly(ethylene glycol diacrylate) resin. Feature sizes below 100 µm are achieved by combining conventional polymer stereolithography with photolithography using photomasks. A layer-by-layer polymerization is performed by alternating two distinct photomasks having horizontally and vertically oriented patterns. Despite the inevitable printing imperfections, minimum reduced plate heights around unity are measured for nonretained analytes. The next challenges for the successful printing of highly efficient and large volume liquid chromatography columns are threefold: reducing the feature size down to below 10 µm, keeping minimum the unevenness of the flow channel dimensions, and tackling additive manufacturing of silica aerogels at such small feature sizes for higher mechanical stability and broader range of retention/selectivity than those delivered by polymer materials.

Model of retention time and density of gradient peak capacity for improved LC-MS method optimization: Application to metabolomics.

A general and deterministic model is derived from the fundamentals of liquid chromatography to calculate retention time, peak width, peak capacity, and density of peak capacity in gradient liquid chromatography. The calculation of these chromatographic properties accounts for 1) the presence of initial (separation of the earliest eluters) and final (column wash) isocratic steps before and after the linear gradient, respectively, 2) the pre- (flow through needle and preheater tubes) and post-column (outlet and emitter tubes before MS detection) dispersion, 3) the compression of the chromatographic band, and 4) the retention of the organic modifier onto the RPLC column. The multiple and variable method parameters may include the column dimensions, particle size, flow rate, temperature, initial and final isocratic hold times, gradient time, gradient steepness, column conditioning/sample load time, and the pre- and post-column tube dimensions. The model enables the users to perform robust multi-dimensional optimization of UHPLC-MS methods and offers the possibility to predict the expected MS feature density for increased method performance. Method optimization can be further improved by matching the observed MS feature density (number of metabolites detected as function of time) to the predicted density of peak capacity. It is directly applied to the optimization of high-throughput RPLC separation methods specifically designed for large-scale urinary metabolic phenotyping.

Extraction of intrinsic column peak profiles of narrow-bore and microbore columns by peak deconvolution methods.

The intrinsic peak profiles (free from the delay and dispersion caused by state-of-the art UHPLC systems) generated by narrow-bore and microbore chromatographic columns used in liquid chromatography-mass spectrometry (LC-MS) proteomic analyses are extracted from two different deconvolution methods. The first method is based on the classical discrete Fourier transform (DFT) while the second method refers to the Taylor expansion of the continuous Fourier transform (FT). The two numerical methods are compared regarding the accurate determination of the intrinsic peak profiles of the non-retained compound (toluene) expected on a narrow-bore 2.1 mm × 100 mm column packed with 1.6 µm CORTECS-C18 superficially porous particles and installed on three different LC systems (ACQUITY i-class UPLC, ACQUITY H-class UPLC, and Arc LC systems). The DFT-based method is most relevant when the low-frequency band of the chromatographic peak does not overlap with the high-frequency bands related to the experimental baseline noise (pump/detector). The Taylor expansion-based method is successful for the extraction of the intrinsic peak profiles of narrow-bore 2.1 mm i.d. columns packed with sub-2 µm particles installed on standard UHPLC systems. When the LC system dispersion significantly exceeds that of the column, the DFT-based method is preferred over the Taylor expansion-based method and is successfully applied to extract the intrinsic peak profiles generated by a microbore 1.0 mm × 100 mm column packed with 1.8 µm HSS-C18 fully porous particles (volume variance ∼ 0.15 µL2 for the non-retained compound toluene) run on the low-dispersion ACQUITY i-class UPLC system (∼ 1 µL2 volume variance). This result opens up promising avenues for the development, quality control, and LC-MS analyses of microbore 1 mm i.d. columns using the state-of-the-art UHPLC instruments at flow rates larger than 0.1 mL/min.

Rebirth of recycling liquid chromatography with modern chromatographic columns : Extension to gradient elution.

Twin column recycling semi-preparative liquid chromatography (TCRLC) is revived to prepare small amount (∼ 1 mg) of a pure targeted compound, which cannot be isolated by conventional preparative liquid chromatography. In this work, TCRLC is extended to gradient elution. The first step of this modified process consists of a gradient step, which eliminates both early and late impurities. If not discarded, some late impurities could echo during the second isocratic recycling step of the process and compromise the purity level required for the targeted compound. Additionally, the entire gradient TCRLC (GTCRLC) process is automated regarding the eluent composition programmed and the actuation times of two valves: one two-position four-port divert valve enables to shave the targeted compound from early and late impurities during the initial gradient step. The second two-position six-port recycling valve ensures the complete baseline resolution between the band of the targeted compound and those of the closest impurities, which are not fully eliminated after the initial gradient step. The automation of the whole GTCRLC process is achieved by running four preliminary scouting gradient runs (at four different relative gradient times, tgt0= 2, 6, 18, and 54, where t0 is the hold-up column time) for the accurate determination of the thermodynamics (lnk versus φ plots of the retention factor as a function of the mobile phase composition) of the first impurity, the targeted compound(s), and of the last impurity. The automated GTCRLC process was successfully applied for the isolation of a polycyclic aromatic hydrocarbon (PAH), chrysene, from a complex mixture of PAHs containing two nearly co-eluting impurities (benzo[a]anthracene and triphenylene) and nine other early/late impurities (sample volume injected: 1 mL, 7.8 mm × 150 mm Sunfire-C18 column, acetonitrile/water eluent mixtures, T= 55 ∘C, 20 cycles, baseline separation in less than two hours). Additionally, the GTCRLC process is advantageously used to isolate and baseline separate the vitamins D2 and D3 initially present in a milk extract mixture (0.3 mL sample injection volume, 7.8 mm × 150 mm Sunfire-C18 column, methanol/water eluent mixtures, T= 65 ∘C, 14 cycles needed in 1.5 hours). These results open promising avenues toward an effective preparation of unknown targeted compounds before further physico-chemical characterization and unambiguous identification.

Theoretical framework for mixer design for noise reduction and gradient fidelity.

The mixing of two or more solvent streams to deliver a stable and accurate solvent composition is crucial to the performance, repeatability and reproducibility of a liquid chromatographic separation. We provide a theoretical treatment of axial mixing of a sequence of solvent packets with the framework of continuous stirred tank reactors (CSTRs) in series and investigate the tradeoffs presented between the primary goal of mixers (noise reduction) and it's necessary side-effects of gradient deformation and asymmetry. An experimental setup to mimic CSTR conditions was created using a stop-flow setup where the fluid flow was periodically paused and sonicated within pods of a certain volume. The effects of mixer volume relative to the volume of pump strokes and gradient volumes were investigated and discussed. A total mixer volume that is six-fold the pump stroke volume was found to be necessary to achieve sufficient (95%) noise reduction necessary for certain applications. A series of two or more CSTR elements was found to outperform a single CSTR element for larger mixer-to-pump stroke volume ratio in dampening baseline noise. For linear gradients, a gradient volume that is ten times larger than the mixer volume was found to sufficiently maintain gradient fidelity. For very small gradient volumes relative to the mixer volume, deformation of linear gradients was found to be significantly greater than predicted by the analytical solution. Furthermore, the nature of the solvent gradient deformation was asymmetric, with the latter half of the solvent gradient deforming significantly more than the first half. Combining analytically and numerically derived solutions for multiple CSTRs connected in series with experimental data, several suggestions can be made on mixer dimensions and design for a certain pump system and method transfer, given a pump stroke volume and gradient time.

Mitigation of analyte loss on metal surfaces in liquid chromatography.

The adsorptive loss of acidic analytes in liquid chromatography was investigated using metal frits. Repetitive injections of acidic small molecules or an oligonucleotide were made on individual 2.1 or 4.6 mm i.d. column frits. Losses were observed for adenosine 5'-(α,ß-methylene) diphosphate, 2-pyridinol 1-oxide and the 25-mer phosphorothioate oligonucleotide Trecovirsen (GEM91) on stainless steel and titanium frits. Analyte adsorption was greatest at acidic pH due to the positive charge on the metal oxide surface. Analyte recovery increased when a series of injections was performed; this effect is known as sample conditioning. Nearly complete recovery was achieved when the metal adsorptive sites were saturated with the analyte. A similar effect was achieved by conditioning the frits with phosphoric, citric or etidronic acids, or their buffered solutions. These procedures can be utilized to mitigate analyte loss. However, the effect is temporary, as the conditioning agent is gradually removed by the running mobile phase. Metal frits modified with hybrid organic/inorganic surface technology were shown to mitigate analyte-to-metal surface interactions and improve recovery of acidic analytes. Quantitative recovery of a 15-35 mer oligodeoxythymidine mixture was achieved using column hardware modified with hybrid surface technology, without a need for column conditioning prior to analysis.

Theoretical study of the efficiency of liquid chromatography columns with particle size gradient.

Modern analytical applications of liquid chromatography require more and more efficient columns. In this work, the possibility of utilizing particle size gradient in the chromatographic column was studied by a theoretical approach. In the course of our work three different scenarios of particle size gradients were considered with different shapes (linear, convex and concave). The evolution of bandwidth inside the column was plotted for each scenario. As a reference point, the bandwidth of the uniform column was used, which had the same pressure drop as the non-uniform column. According to our calculations, in isocratic elution mode, the non-uniform column does not offer any advantage compared to the uniform column, regardless the type of the particle size gradient. In gradient elution mode, however, extra band compression occurs was found. For negative particle size gradients, the final physical bandwidth was found to be approximately 1-4 % smaller than for uniform columns. This slight gain in efficiency in terms of bandwidth compression can be expanded to 5-8 % by the optimization of the limiting particle sizes. These optimized results are obtained when the final particle size is approximately 40% of the initial particle diameter.

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 V0, system dwell volume Vdwell (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 (V0, Vdwell, 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.

Multiple-open-tubular column enabling transverse diffusion. Part 2: Channel size distribution and structure optimization.

Multiple-open-tubular columns enabling transverse diffusion (MOTTD) are made of straight, parallel, and cylindrical flow channels separated by a mesoporous stationary phase. In Part 1, a model of band broadening along MOTTD columns accounting for longitudinal diffusion, the trans-channel velocity bias, and mass transfer resistance in the stationary phase was proposed and validated. In this Part 2, the model is completed by considering the impact of short-range inter-channel velocity biases on the MOTTD plate number. These velocity biases are caused by the wide distribution of the channel diameters. Different ratios, ρ, of the average inner diameter, 2, of the flow channels to their closest center-to-center distance d (d= 5 µm, ρ= 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, and 0.9) with a relative standard deviation (RSD) increasing from 0 to 50% are considered. The zone retention factor k1 was increased from 1 to 25. The complete model of band broadening is validated after adjustment to dispersion data obtained by 1) the lattice-Boltzmann method for modeling fluid flow, 2) a random-walk particle-tracking (RWPT) technique to address advective-diffusive transport, and 3) by considering two distinct populations of flow channels (inner radii rc,1=(1-RSD) and rc,2=(1+RSD)) arranged at the nodes of a hexagonal compact array. The completed model of band broadening in MOTTD columns reveals that the RSD of the channel diameters has only a moderate impact on the optimum plate number of MOTTD columns: the relative increase of the minimum plate height do not exceed 30% even for the largest RSDs. However, when the mass transfer of the analyte is governed by its slow rate of transverse diffusion across the MOTTD column, the plate height can be increased by up to 100% at high average velocities. Regarding the best trade-off between analysis speed and column performance at a fixed pressure drop of 400 bar, irrespective of the zone retention factor and RSD of the distribution of the channel diameters, the fastest analyses are recommended for MOTTD columns having a small structural parameter ρ. In contrast, for the longest analysis times, the largest values of ρ are required to maximize the performance of MOTTD columns.

Theoretical performance of multiple size-exclusion chromatography columns connected in series.

The fundamental relationships are derived for the retention, peak width, and peak capacity of non-retained polymers eluting from multiple standard size-exclusion chromatography (SEC) columns connected in series. The standard SEC columns may have different dimensions and are packed with particles having distinct average particle diameters (APDs) and average mesopore sizes (AMSs). The performances (peak capacity, local resolution power, and sensitivity) of three standard SEC columns connected in series (called a tri-SEC column) packed with bridged-ethylene-hybrid (BEH) fully porous particles (FPPs) having three different APDs (1.7, 2.5, and 3.5 µm) and AMSs (200, 450, and 900 Å, respectively) are calculated as a function of the applied flow rate and size of polystyrene standards. Irrespective of the APD and AMS, the present investigation assumes isomorphological materials relative to the mesopore space of the three different BEH particles. The advantage of a 15 cm long tri-SEC column over a single reference SEC column (APD=3.5 µm, AMS=900 Å), which generates the same back pressure and separation window as those of the tri-SEC column, is expected at flow rates larger than the optimum flow rate generating the maximum peak capacity. The calculations predict a significant relative increase of the peak capacity (from +25% to +85%), resolution of small molecules (from +75% to +225%), and of the detection limit of intermediate size (from +15% to +70%) and largest polymers (from +25 to +110%). This is explained by 1) the exclusion of the largest polymers from the internal volume of the particles having the smallest mesopores (restricted access media) and 2) the minimum dispersion along the columns packed with the smallest particle sizes in the tri-SEC column. The main benefit of multi-SEC columns is to easily adjust the desired pore size distribution by properly selecting the lengths of each individual SEC column. The user can then control the pore size distribution for any specific separation problem. A potential application is theoretically demonstrated for the fast purification of monoclonal antibodies from metabolites, host cell proteins, aggregated forms, and from virus-like particles.