Fundamental investigation of impact of water and TFA additions in peptide sub/supercritical fluid separations.

The retention of three peptides was studied under analytical and overloaded conditions at different concentrations of trifluoroacetic acid (TFA) and water added to the co-solvent methanol (MeOH). Four columns with different stationary phase properties, i.e., silica, diol, 2-ethylpyridine and cyanopropyl (CN) columns, were evaluated in this investigation. The overall aim was to get a deeper understanding on how column chemistry as well as water and TFA in the co-solvent affect the analytical and overloaded elution profiles using multivariate design of experiments and adsorption measurements of co-solvent components. Multivariate experimental design modeling indicated that water had on average around five times higher effect on the retention than the addition of TFA. The results also showed that the retention increases with the addition of TFA and water to the co-solvent on all columns except the CN column, on which the retention decreased. When examining the effect of adding water to the co-solvent, evidence of a hydrophilic interaction liquid chromatography (HILIC)-like retention mechanism was found on the three other columns with more polar stationary phases. However, on the CN column water acted as an additive, decreasing the retention due to competition with the peptide for available adsorption surface. Adsorption isotherm measurements of the polar solvent MeOH showed that MeOH adsorbs much weaker on the CN column than on the other columns. Addition of TFA and water to the co-solvent substantially sharpened the elution profiles under both overloaded and analytical conditions. Adding a small amount of TFA (from 0 % to 0.05 %) to the co-solvent substantially improved the peak shape of the elution profiles, while further addition (from 0.05 % to 0.15 %) had only a minor effect on the elution profile shape. The reduced retention on the CN column could not be explained by TFA adsorption, which was very weak on all studied columns (retention factor, 0.05-0.15). One could therefore speculate that the ion-pairing complex formed between the peptide and TFA in the mobile phase, reduce the retention due to its reduced polarity. On the other columns displaying HILIC-like properties, the TFA probably just decreased the pH of the mobile phase, thereby promoting the partitioning of the peptide into the water-rich layer. Finally, peak deformation due to diluent-eluent mismatch was observed under overloaded conditions. This was most severe in the cases where MeOH adsorption to the stationary phase was strong and the peptides were only mildly retained. Adding 1,4-dioxan to the diluent resolved this issue.

Prediction of overloaded concentration profiles under ultra-high-pressure liquid chromatographic conditions.

In this study, overloaded elution profiles under ultra-high-pressure liquid chromatographic (UHPLC) conditions and accounting for the severe pressure and temperature gradients generated, are compared with experimental data. The model system consisted of an C18 column packed with 1.7-µm particles (i.e., a UHPLC column) and the solute was 1,3,5-tri­tert-butylbenzene eluted with a mobile phase composed of 85/15 (v/v) acetonitrile/water. Two thermal modes were considered, and the solute was eluted at the very high inlet pressures necessary to achieve a highly efficient and rapid chromatographic process, as provided by using columns packed with small particles. However, the high inlet pressure and high linear velocity of the mobile phase caused the production of a significant amount of heat, and consequently, the formation of axial and radial temperature gradients. Due to these gradients, the retention and the mobile phase velocity were no longer constant. Thus, simple mathematical models consisting only of the mass balance equations are unsuitable to properly model the elution profiles. Here, the elution concentration profiles were predicted using a combined two-dimensional heat and mass transfer model, also including the calculation of the mobile phase velocity distribution. The isotherm adsorption model was the bi-Langmuir isotherm model with Henry constants that depended on the local temperature and pressure in the column. These adjustments allowed us to precisely account for changes in the shape and retention of the overloaded concentration profiles in the mobile phase. The proposed model provided accurate predictions of the overloaded concentration profiles, demonstrating good agreement with experimental profiles eluted under severe pressure and temperature gradients in the column even in the most extreme cases where the pressure drops reached 846 bar and the temperature gradients equaled 0.15 K mm-1 and 0.95 K mm-1 in the axial and the radial directions, respectively. In such cases 36 % decrease of the retention factor was observed along the column and 2 % increase in radial direction. These changes, combined with the velocity distribution, shifted the overloaded elution profile's shock towards the center of the column, advancing approximately 3 mm from its initial position close to the column wall. Ultimately, this resulted in the broadening of the elution band.

Strategies for predictive modeling of overloaded oligonucleotide elution profiles in ion-pair chromatography.

Due to their potential for gene regulation, oligonucleotides have moved into focus as one of the preferred modalities modulating currently undruggable disease-associated targets. In the course of synthesis and storage of oligonucleotides a significant number of compound-related impurities can be generated. Purification protocols and analytical methods have become crucial for the therapeutic application of any oligonucleotides, be they antisense oligonucleotides (ASOs), small interfering ribonucleic acids (siRNAs) or conjugates. Ion-pair chromatography is currently the standard method for separating and analyzing therapeutic oligonucleotides. Although mathematical modeling can improve the accuracy and efficiency of ion-pair chromatography, its application remains challenging. Simple models may not be suitable to treat advanced single molecules, while complex models are still inefficient for industrial oligonucleotide optimization processes. Therefore, fundamental research to improve the accuracy and simplicity of mathematical models in ion-pair chromatography is still a necessity. In this study, we predict overloaded concentration profiles of oligonucleotides in ion-pair chromatography and compare relatively simple and more advanced predictive models. The experimental system consists of a traditional C18 column using (dibutyl)amine as the ion-pair reagent and acetonitrile as organic modifier. The models were built and tested based on three crude 16-mer oligonucleotides with varying degrees of phosphorothioation, as well as their respective n - 1 and (P = O)1 impurities. In short, the proposed models were suitable to predict the overloaded concentration profiles for different slopes of the organic modifier gradient and column load.

Capillary electrokinetic chromatography for studying interactions between ß-blockers and Intralipid emulsion.

Toxicity of ß-blockers is one of the most common causes of poison-induced cardiogenic shock throughout the world. Therefore, methodologies for in vivo removal of the drugs from the body have been under investigation. Intralipid emulsion (ILE) is a common commercial lipid emulsion used for parenteral nutrition, but it has also been administered to patients suffering from drug toxicities. In this work, a set of ß-blockers of different hydrophobicity's (log KD values ranging from 0.16 to 3.8) were investigated. The relative strength of the interactions between these compounds and the ILE was quantitatively assessed by means of binding constants and adsorption constants of the formed ß-blocker-ILE complexes. The binding constants were determined by capillary electrokinetic chromatography and the adsorption constants were calculated based on different adsorption isotherms. Expectedly, the binding constants were strongly related to the log KD values of the ß-blockers. The binding and adsorption constants also show that less hydrophobic ß-blockers interact with ILE, suggesting that this emulsion could be useful for capturing such compounds in cases of their overdoses. Thus, the use of ILE for treatment of toxicities caused by a larger range of ß-blockers is worth further investigation.

Assessment of physicochemical properties of sorbent materials in passive and active sampling systems towards gaseous nitrogen-containing compounds.

The adsorption and desorption behavior of volatile nitrogen-containing compounds in vapor phase by solid-phase microextraction Arrow (SPME-Arrow) and in-tube extraction (ITEX) sampling systems, were investigated experimentally using gas chromatography-mass spectrometry. Three different SPME-Arrow coating materials, DVB/PDMS, MCM-41, and MCM-41-TP and two ITEX adsorbents, TENAX-GR and MCM-41-TP were compared to clarify the selectivity of the sorbents towards nitrogen-containing compounds. In addition, saturated vapor pressures for these compounds were estimated, both experimentally and theoretically. In this study, the adsorption of nitrogen-containing compounds on various adsorbents followed the Elovich model well, while a pseudo-first-order kinetics model best described the desorption kinetics. Pore volume and pore sizes of the coating sorbents were essential parameters for the determination of the adsorption performance for the SPME-Arrow sampling system. MCM-41-TP coating with the smallest pore size gave the slowest adsorption rate compared to that of DVB/PDMS and MCM-41 in the SPME-Arrow sampling system. Both adsorbent and adsorbate properties, such as hydrophobicity and basicity, affected the adsorption and desorption kinetics in SPME-Arrow system. The adsorption and desorption rates of studied C6H15N isomers in the MCM-41 and MCM-41-TP sorbent materials of SPME-Arrow system were higher for dipropylamine and triethylamine (branched amines) than for hexylamine (linear chain amines). DVB/PDMS-SPME-Arrow gave fast adsorption rates for the aromatic-ringed pyridine and o-toluidine. All studied nitrogen-containing compounds demonstrated high desorption rates with DVB/PDMS-SPME-Arrow. Chemisorption and physisorption were the sorption mechanisms in MCM-41- and MCM-41-TP- SPME-Arrow, but additional experiments are needed to confirm this. An active sampling technique ITEX gave comparable adsorption and desorption rates on the selective MCM-41-TP and universal TENAX-GR sorbent materials for all the compounds studied. Vapor pressures of nitrogen-containing compounds were experimentally estimated by using retention index approach and these values were compared with the theoretical ones, calculated using the COnductor-like Screening MOdel for Real Solvent (COSMO-RS) model. Both values agreed well with those found in the literature proving that these methods can be successfully used in predicting VOC's vapor pressures, e.g. for the formation of secondary organic aerosols.

Development of a unified gradient theory for ion-pair chromatography using oligonucleotide separations as a model case.

Ion-pair chromatography is the de facto standard for separating oligonucleotides and related impurities, particularly for analysis but also often for small-scale purification. Currently, there is limited understanding of the quantitative modeling of both analytical and overloaded elution profiles obtained during gradient elution in ion-pair chromatography. Here we will investigate a recently introduced gradient mode, the so-called ion-pairing reagent gradient mode, for both analytical and overloaded separations of oligonucleotides. The first part of the study demonstrates how the electrostatic theory of ion-pair chromatography can be applied for modeling gradient elution of oligonucleotides. When the ion-pair gradient mode is used in a region where the electrostatic surface potential can be linearized, a closed-form expression of retention time can be derived. A unified retention model was then derived, applicable for both ion-pair reagent gradient mode as well as co-solvent gradient mode. The model was verified for two different experimental systems and homo- and heteromeric oligonucleotides of different lengths. Quantitative modeling of overloaded chromatography using the ion-pairing reagent gradient mode was also investigated. Firstly, a unified adsorption isotherm model was developed for both gradient modes. Then, adsorption isotherms parameter of a model oligonucleotide and two major synthetic impurities were estimated using the inverse method. Secondly, the parameters of the adsorption isotherm were then used to investigate how the productivity of oligonucleotide varies with injection volume, gradient slope, and initial retention factor. Here, the productivity increased when using a shallow gradient slope combined with a low initial retention factor. Finally, experiments were conducted to confirming some of the model predictions. Comparison with the conventional co-solvent gradient mode showed that the ion-pairing reagent gradient leads to both higher yield and productivity while consuming less co-solvent.

A closer study of overloaded elution bands and their perturbation peaks in ion-pair chromatography.

There is strong renewed interest in ion-pair chromatography (IPC) because of its great importance for separating new-generation biosimilar pharmaceuticals such as oligonucleotides. Due to the complexity of the IPC process, its mathematical modeling is challenging, especially in preparative mode. In a recent study, Lesko et al. (2021) developed a mathematical model for predicting, with good accuracy, overloaded concentration profiles for sodium benzenesulfonate, describing how the overloaded solute concentration profiles change from Langmuirian to complicated U-shaped, and then back again to Langmuirian profiles, with increasing concentration of the ion-pair reagent in the mobile phase. This study identifies and explains the underlying mechanism generating these complex peak shapes and band-shape transformations; this was only possible by visualizing and modeling the underlying equilibrium perturbations that occur upon injection in preparative IPC. In the 2021 study, the model was derived based on the concentration profiles obtained using a conventional UV detector principle, so the concentration gradients and perturbation zones of the mobile-phase components were not visualized. In this study, the necessary mechanistic information was obtained via complementary experiments combining two detection principles, i.e., refractive index detection and UV detection, with modeling efforts. The models correctly described the invisible equilibrium perturbations and how these formed internal gradients of the mobile-phase components. The models also explained the complex overloaded solute-band deformations reported in the recent study. In addition, a rule of thumb was developed for predicting experimental conditions that could result in deformed solute elution profiles and/or for avoiding these deformations. The latter is crucial for the practical chromatographer, since such U-shaped solute-band profiles are undesirable in preparative separation due to the broader elution zones, resulting in lower productivity than that of normal band shapes.

Building machine-learning-based models for retention time and resolution predictions in ion pair chromatography of oligonucleotides.

Support vector regression models are created and used to predict the retention times of oligonucleotides separated using gradient ion-pair chromatography with high accuracy. The experimental dataset consisted of fully phosphorothioated oligonucleotides. Two models were trained and validated using two pseudo-orthogonal gradient modes and three gradient slopes. The results show that the spread in retention time differs between the two gradient modes, which indicated varying degree of sequence dependent separation. Peak widths from the experimental dataset were calculated and correlated with the guanine-cytosine content and retention time of the sequence for each gradient slope. This data was used to predict the resolution of the n - 1 impurity among 250 000 random 12- and 16-mer sequences; showing one of the investigated gradient modes has a much higher probability of exceeding a resolution of 1.5, particularly for the 16-mer sequences. Sequences having a high guanine-cytosine content and a terminal C are more likely to not reach critical resolution. The trained SVR models can both be used to identify characteristics of different separation methods and to assist in the choice of method conditions, i.e. to optimize resolution for arbitrary sequences. The methodology presented in this study can be expected to be applicable to predict retention times of other oligonucleotide synthesis and degradation impurities if provided enough training data.

Kinetics and interaction studies of anti-tetraspanin antibodies and ICAM-1 with extracellular vesicle subpopulations using continuous flow quartz crystal microbalance biosensor.

Continuous flow quartz crystal microbalance (QCM) was utilized to study binding kinetics between EV subpopulations (exomere- and exosome-sized EVs) and four affinity ligands: monoclonal antibodies against tetraspanins (anti-CD9, anti-CD63, and anti-CD81) and recombinant intercellular adhesion molecule-1 (ICAM-1) or CD54 protein). High purity CD9+, CD63+, and CD81+ EV subpopulations of <50 nm exomeres and 50-80 nm exosomes were isolated and fractionated using our recently developed on-line coupled immunoaffinity chromatography - asymmetric flow field-flow fractionation system. Adaptive Interaction Distribution Algorithm (AIDA), specifically designed for the analysis of complex biological interactions, was used with a four-step procedure for reliable estimation of the degree of heterogeneity in rate constant distributions. Interactions between exomere-sized EVs and anti-tetraspanin antibodies demonstrated two interaction sites with comparable binding kinetics and estimated dissociation constants Kd ranging from nM to fM. Exomeres exhibited slightly higher affinity compared to exosomes. The highest affinity with anti-tetraspanin antibodies was achieved with CD63+ EVs. The interaction of EV subpopulations with ICAM-1 involved in cell internalization of EVs was also investigated. EV - ICAM-1 interaction was also of high affinity (nM to pM range) with overall lower affinity compared to the interactions of anti-tetraspanin antibodies and EVs. Our findings proved that QCM is a valuable label-free tool for kinetic studies with limited sample concentration, and that advanced algorithms, such as AIDA, are crucial for proper determination of kinetic heterogeneity. To the best of our knowledge, this is the first kinetic study on the interaction between plasma-derived EV subpopulations and anti-tetraspanin antibodies and ICAM-1.

Biosensor-Enabled Deconvolution of the Avidity-Induced Affinity Enhancement for the SARS-CoV-2 Spike Protein and ACE2 Interaction.

Avidity is an effective and frequent phenomenon employed by nature to achieve extremely high-affinity interactions. As more drug discovery efforts aim to disrupt protein-protein interactions, it is becoming increasingly common to encounter systems that utilize avidity effects and to study these systems using surface-based technologies, such as surface plasmon resonance (SPR) or biolayer interferometry. However, heterogeneity introduced from multivalent binding interactions complicates the analysis of the resulting sensorgram. A frequently applied practice is to fit the data based on a 1:1 binding model, and if the fit does not describe the data adequately, then the experimental setup is changed to favor a 1:1 binding interaction. This reductionistic approach is informative but not always biologically relevant. Therefore, we aimed to develop an SPR-based assay that would reduce the heterogeneity to enable the determination of the kinetic rate constants for multivalent binding interactions using the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein and the human receptor angiotensin-converting enzyme 2 (ACE2) as a model system. We employed a combinatorial approach to generate a sensor surface that could distinguish between monovalent and multivalent interactions. Using advanced data analysis algorithms to analyze the resulting sensorgrams, we found that controlling the surface heterogeneity enabled the deconvolution of the avidity-induced affinity enhancement for the SARS-CoV-2 spike protein and ACE2 interaction.

Experimental and theoretical investigation of high- concentration elution bands in ion-pair chromatography.

The effective separation of many solutes, including pharmaceuticals, can be performed using an ion-pair reagent (IPR) in the mobile phase. However, chromatographic separation and mathematical modelling are a challenge in ionpair chromatography (IPC), especially in preparative mode, due to the complicated chromatographic process. In this study, we present a retention mechanism and a mathematical model that predict overloaded concentration profiles in IPC using a system with X-Bridge C18 as stationary phase and tetrabutylammonium bromide in the 0 - 15 mM concentration range as the IPR. Two different mobile phases were used: (i) 15/85 [v/v] acetonitrile/water, (ii) 25/75 methanol/water. The model compounds were sodium salts of organic compounds with sulfonic acid functions. The analytical and preparative elution profiles were obtained for specified conditions. The analytical data were utilized to calculate the difference in electrical potential between the surface and bulk solution using firm electrostatic theory. In the preparative mode in a certain range of IPR concentrations, complicated U-shaped overloaded profiles were observed. In the other considered cases, Langmuir overloaded elution profiles were recorded. A multilayer adsorption model was derived, which is consistent with the dynamic ion exchange models. The model assumes that lipophilic IPR adsorbs on the stationary phase, creating charged active sites that serve as exchange sites for the solutes. The molecules of the solute can adsorb on the already formed IPR layer. It was also assumed that a subsequent layer of solute can form on the formed layer of complexes due to interactions between the solute molecules. The model takes into account the electrostatic attraction and repulsion of the molecules, depending on the considered situation. The proposed model allowed prediction of the overloaded concentration profiles with very good agreement for the model solute and followed the progression from Langmuirian, through U-shaped, to again Langmuirian profiles.

Selectivity limits of and opportunities for ion pair chromatographic separation of oligonucleotides.

Here it was investigated how oligonucleotide retention and selectivity factors are affected by electrostatic and non-electrostatic interactions in ion pair chromatography. A framework was derived describing how selectivity depends on the electrostatic potential generated by the ion-pair reagent concentration, co-solvent volume fraction, charge difference between the analytes, and temperature. Isocratic experiments verified that, in separation problems concerning oligonucleotides of different charges, selectivity increases with increasing surface potential and analyte charge difference and with decreasing co-solvent volume fraction and temperature. For analytes of the same charge, for example, diastereomers of phosphorothioated oligonucleotides, selectivity can be increased by decreasing the co-solvent volume fraction or the temperature and has only a minor dependency on the ion-pairing reagent concentration. An important observation is that oligonucleotide retention is driven predominantly by electrostatic interaction generated by the adsorption of the ion-pairing reagent. We therefore compared classical gradient elution in which the co-solvent volume fraction increases over time versus gradient elution with a constant co-solvent volume fraction but with decreasing ion-pair reagent concentration over time. Both modes decrease the electrostatic potential. Oligonucleotide selectivity was found to increase with decreasing ion-pairing reagent concentration. The two elution modes were finally applied to two different model antisense oligonucleotide separation problems, and it was shown that the ion-pair reagent gradient increases the selectivity of non-charge-based separation problems while maintaining charge-difference-based selectivity.

A Retention-Matching Strategy for Method Transfer in Supercritical Fluid Chromatography: Introducing the Isomolar Plot Approach.

A strategy to match any retention shifts due to increased or decreased pressure drop during supercritical fluid chromatography (SFC) method transfer is presented. The strategy relies on adjusting the co-solvent molarity without the need to adjust the back-pressure regulator. Exact matching can be obtained with minimal changes in separation selectivity. To accomplish this, we introduce the isomolar plot approach, which shows the variation in molar co-solvent concentration depending on the mass fraction of co-solvent, pressure, and temperature, here exemplified by CO2-methanol. This plot allowed us to unify the effects of the co-solvent mass fraction and density on retention in SFC. The approach, which was verified on 12 known empirical retention models for each enantiomer of six basic pharmaceuticals, allowed us to numerically calculate the apparent retention factor for any column pressure drop. The strategy can be implemented either using a mechanistic approach if retention models are known or empirically by iteratively adjusting the co-solvent mass fraction. As a rule of thumb for the empirical approach, we found that the relative mass fraction adjustment needed is proportional to the relative change in the retention factor caused by a change in the pressure drop. Different proportionality constants were required to match retention in the case of increasing or decreasing pressure drops.

Predictions of overloaded concentration profiles in supercritical fluid chromatography.

Here, overloaded concentration profiles were predicted in supercritical fluid chromatography using a combined two-dimensional heat and mass transfer model. The heat balance equation provided the temperature and pressure profiles inside the column. From this the density, viscosity, and mobile phase velocity profiles in the column were calculated. The adsorption model is here expressed as a function of the density and temperature of the mobile phase. The model system consisted of a Kromasil Diol column packed with 2.2-µm particles (i.e., a UHPSFC column) and the solute was phenol eluted with neat carbon dioxide at three different outlet pressures and five different mobile phase flow rates. The proposed model successfully predicted the eluted concentration profiles in all experimental runs with good agreement even with high-density drops along the column. It could be concluded that the radial temperature and density gradients did not significantly influence the overloaded concentration elution profiles.

Impact of stationary-phase pore size on chromatographic performance using oligonucleotide separation as a model.

A combined experimental and theoretical study was performed to understand how the pore size of packing materials with pores 60-300 Å in size affects the separation of 5-50-mer oligonucleotides. For this purpose, we developed a model in which the solutes were described as thin rods to estimate the accessible surface area of the solute as a function of the pore size and solute size. First, an analytical investigation was conducted in which we found that the selectivity increased by a factor of 2.5 when separating 5- and 15-mer oligonucleotides using packing with 300 Å rather than 100 Å pores. We complemented the analytical investigation by theoretically demonstrating how the selectivity is dependent on the column's accessible surface area as a function of solute size. In the preparative investigation, we determined adsorption isotherms for oligonucleotides using the inverse method for separations of a 9- and a 10-mer. We found that preparative columns with a 60 Å-pore-size packing material provided a 10% increase in productivity as compared with a 300 Å packing material, although the surface area of the 60 Å packing is as much as five time larger.

Impact of Methanol Adsorption on the Robustness of Analytical Supercritical Fluid Chromatography in Transfer from SFC to UHPSFC.

In supercritical fluid chromatography (SFC), the retention of a solute depends on the temperature, density, pressure, and cosolvent fraction. Here, we investigate how the adsorption of the cosolvent MeOH changes with pressure and temperature and how this affects the retention of several solutes. The lower the pressure, the stronger the MeOH adsorption to the stationary phase; in addition, at low pressure, perturbing the pressure results in significant changes in the amounts of MeOH adsorbed to the stationary phase. The robustness of the solute retention was lowest when operating the systems at low pressures, high temperatures, and low cosolvent fractions in the eluent. Here, we found a clear relationship between the sensitivity of MeOH adsorption to the stationary phase and the robustness of the separation system. Finally, we show that going from classical SFC to ultrahigh-performance SFC (UHPSFC), that is, separations conducted with much smaller packing diameters, results in retention factors that are more sensitive to fluctuations in the flow rate than with traditional SFC. The calculated density profiles indicate only a slight density drop over the traditional SFC column (3%, visualized as lighter → darker blue in the TOC), whereas the drop for the UHPSFC one was considerably larger (20%, visualized as dark red → light green in the TOC). The corresponding temperature drops were calculated to be 0.8 and 6.5 °C for the SFC and UHPSFC systems, respectively. These increased density and temperature drops are the underlying reasons for the decreased robustness of UHPSFC.

Advanced Analysis of Biosensor Data for SARS-CoV-2 RBD and ACE2 Interactions.

The traditional approach for analyzing interaction data from biosensors instruments is based on the simplified assumption that also larger biomolecules interactions are homogeneous. It was recently reported that the human receptor angiotensin-converting enzyme 2 (ACE2) plays a key role for capturing SARS-CoV-2 into the human target body, and binding studies were performed using biosensors techniques based on surface plasmon resonance and bio-layer interferometry. The published affinity constants for the interactions, derived using the traditional approach, described a single interaction between ACE2 and the SARS-CoV-2 receptor binding domain (RBD). We reanalyzed these data sets using our advanced four-step approach based on an adaptive interaction distribution algorithm (AIDA) that accounts for the great complexity of larger biomolecules and gives a two-dimensional distribution of association and dissociation rate constants. Our results showed that in both cases the standard assumption about a single interaction was erroneous, and in one of the cases, the value of the affinity constant KD differed more than 300% between the reported value and our calculation. This information can prove very useful in providing mechanistic information and insights about the mechanism of interactions between ACE2 and SARS-CoV-2 RBD or similar systems.

Evaluating the advantages of higher heat conductivity in a recently developed type of core-shell diamond stationary phase particle in UHPLC.

In recent studies, the nature and magnitude of the temperature gradients developed in ultra-high pressure liquid chromatography (UHPLC), were found to be dependent on the heat conductivity properties of the column matrices, but also, on the principle used for controlling the temperature over the column. Here, we investigated the potential of using highly heat conductive diamond-based stationary phases (85 times higher than silica), for reducing the temperature gradients. The stationary phases investigated were a (i) Diamond Analytics FLARE column, based on particles comprised of a graphite core surrounded by a very thin diamond shell, and two silica hybrid columns: (ii) a core-shell silica Kromasil Eternity Shell column and (iii) a fully porous silica Kromasil Eternity XT column. Models were developed based on two-dimensional heat transfer theory and mass transfer theory, which were used to model the temperature profiles and the migration of an analyte band accounting for column efficiencies at different flow rates. For the silica-based columns, using water-controlled temperature mode, the temperature gradients along the column axes are suppressed whereas temperature gradients in the radial direction prevails resulting in decreased column efficiencies. Using these columns with air-controlled temperature mode, the radial temperature gradients are reduced whereas temperature gradients along the column prevails resulting in decreased retention times. With the Diamond FLARE column, there was no loss in column efficiency using the water-controlled temperature mode and the van Deemter curves are almost identical using both temperature control modes. Thus, for the Diamond FLARE column, in contrast to the silica-based columns, there are almost no losses of column efficiencies due to reduced radial temperature gradients independent on how the column temperature was controlled.

Systematic investigations of peak distortions due to additives in supercritical fluid chromatography.

The impact of eluent components added to improve separation performance in supercritical fluid chromatography was systematically, and fundamentally, investigated. The model system comprised basic pharmaceuticals as solutes and eluents containing an amine (i.e., triethylamine, diethylamine, or isopropylamine) as additive with MeOH as the co-solvent. First, an analytical-scale study was performed, systematically investigating the impact of the additives/co-solvent on solute peak shapes and retentions, using a design of experiments approach; here, the total additive concentration in the eluent ranged between 0.021 and 0.105 % (v/v) and the MeOH fraction in the eluent between 16 and 26 % (v/v). The co-solvent fraction was found to be the most efficient tool for adjusting retentions, whereas the additive fraction was the prime tool for improving column efficiency and peak analytical performance. Next, the impacts of the amine additives on the shapes of the so-called overloaded solute elution profiles were investigated. Two principal types of preparative peak deformations appeared and were investigated in depth, analyzed using computer simulation with mechanistic modeling. The first type of deformation was due to the solute eluting too close to the additive perturbation peak, resulting in severe peak deformation caused by co-elution. The second type of deformation was also due to additive-solute interactions, but here the amine additives acted as kosmotropic agents, promoting the multilayer adsorption to the stationary phase of solutes with bulkier aryl groups.

Analytical and preparative separation of phosphorothioated oligonucleotides: columns and ion-pair reagents.

Oligonucleotide drugs represent an emerging area in the pharmaceutical industry. Solid-phase synthesis generates many structurally closely related impurities, making efficient separation systems for purification and analysis a key challenge during pharmaceutical drug development. To increase the fundamental understanding of the important preparative separation step, mass-overloaded injections of a fully phosphorothioated 16mer, i.e., deoxythymidine oligonucleotide, were performed on a C18 and a phenyl column. The narrowest elution profiles were obtained using the phenyl column, and the 16mer could be collected with high purity and yield on both columns. The most likely contribution to the successful purification was the quantifiable displacement of the early-eluting shortmers on both columns. In addition, the phenyl column displayed better separation of later-eluting impurities, such as the 17mer impurity. The mass-overloaded injections resulted in classical Langmuirian elution profiles on all columns, provided the concentration of the ion-pairing reagent in the eluent was sufficiently high. Two additional column chemistries, C4 and C8, were also investigated in terms of their selectivity and elution profile characteristics for the separation of 5-20mers fully phosphorothioated deoxythymidine oligonucleotides. When using triethylamine as ion-pairing reagent to separate phosphorothioated oligonucleotides, we observed peak broadening caused by the partial separation of diastereomers, predominantly seen on the C4 and C18 columns. When using the ion-pair reagent tributylamine, to suppress diastereomer separation, the greatest selectivity was found using the phenyl column followed by C18. The present results will be useful when designing and optimizing efficient preparative separations of synthetic oligonucleotides.