Enantioseparation and molecular docking study of selected chiral pharmaceuticals on a commercialized phenylcarbamate-β-cyclodextrin column using polar organic mode

The chiral separation capability of Chiral-CD-Ph column, containing phenylcarbamate-β-cyclodextrin as the chiral selector in polar organic mode was investigated. A total of twenty-five compounds with different structures and acid–base properties were evaluated, and twenty of them were separated using acetonitrile or methanol as eluent. The effects of various chromatographic parameters, such as the type and proportion of organic modifier, flow rate, and column temperature were analyzed in detail in relation to chromatographic performance. A U-shape retention curve was observed when a mixture of acetonitrile and methanol was used as the eluent, indicating different types of interactions in different solvent mixtures. Van 't Hoff analysis was used for calculation of thermodynamic parameters which revealed that the enantioseparation is mainly enthalpy controlled; however, entropic control was also observed. The enantiomer recognition ability at the atomic level was also investigated through a molecular docking study, which revealed surface binding in polar organic mode instead of inclusion complexation. Our work proves that the phenylcarbamate-β-cyclodextrin-based chiral stationary phase can be effectively used in polar organic mode for the chiral separation of structurally diverse compounds. Furthermore, it is important to note that our study demonstrated that surface binding is responsible for the formation of supramolecular complexes in certain cyclodextrin derivatives.

www.nature.com/scientificreports/Cyclodextrins (CDs) are among one of the most widely used chiral selectors.CDs are cyclic oligosaccharides formed by d-glucopyranosyl units linked through 1,4-linkages.These macrocycles feature a toroidal shape with a hydrophobic interior cavity and hydrophilic exterior.Due to their chiral nature, they are capable of enantiodiscrimination through transient diastereomeric complex formation either through inclusion complexation or surface binding 6 .In liquid chromatographic enantioseparations, the choice of mobile phase can impact the binding mode to the chiral selector.In an aqueous system, inclusion complexation is more likely.The hydrophobic cavity of the CD provides a confined environment for the lipophilic part of the guest molecule.This inclusion effect shields the guest molecule from the surrounding aqueous environment, leading to a decrease in its entropy.When CD-based CSPs are used in polar organic or normal-phase modes, the inner cavity is obstructed by solvent molecules, preventing inclusion complexation.Nevertheless, in such media, hydrophilic interactions may be enhanced, where solutes with hydrophilic groups bind to the polar surface of the CD.Nevertheless, without the application of additional techniques such as molecular modeling or spectroscopy, determining the binding mode is not possible 5 .Diverse types of CDs are available, including native-and derivatized CDs.The latter were developed to offer additional intermolecular interactions, such as π-π interactions, hydrogen bonding, dipole-dipole interactions, and ion-pairing, resulting in improved enantioseparation capabilities 7,8 .Although there is a continuous interest in the development and implementation of new CD derivatives as silica bound CSPs, it is worth noting that only a limited number of these are commercially available 9 .Armstrong et al. made significant contributions to the field by developing both native and derivatized α-, β-, and γ-CD CSPs, which enabled the efficient enantioseparation of dansylated amino acids, barbiturates, as well as aromatic and aliphatic sulfoxides through high-performance liquid chromatography (HPLC) [10][11][12] .After Armstrong's research group demonstrated that CD-based stationary phases are suitable for the direct separation of enantiomers, numerous new developments were initiated.Wang et al. successfully developed native β-CD, permethylated β-CD, and perphenylcarbamated β-CD CSPs for the enantioseparation of several dansyl amino acids, flavonoids, β-blockers, and α-ionone derivatives 13 .Additionally, Tang's group has studied the enantiomeric separation of various analytes including flavonoids, aromatic alcohols, β-blockers, and amino acids using three β-CD CSPs prepared by clicking different per-chloro-methyl-phenylcarbamated moiety onto a silica support 14,15 .In a recent article, Li et al. prepared a novel 3,5-dichloro-phenylcarbamate-β-CD chemically bonded CSP and used it for the enantioseparation of different proton-pump inhibitors using normal-phase LC 9 .With the help of the 3,5-dimethylphenylcarbamate-β-CD column, Rezenka et al. were able to determine the presence of two diacylglycerol enantiomers in bacteria, suggesting the existence of heterochiral membranes 16 .At present native CD-based columns as well as acetylated β-, permethylated β-, dimethylated β-hydroxypropylated β-, napthylethylcarbamated β-, 3,5-dimethylphenylcarbamate-β-, and perphenylcarbamate-β-CD are commercially available on the market as CSPs.
The Chiral CD-Ph column features a phenylcarbamated β-CD-modified silica CSP.It is produced by Osaka Soda (Shiseido) and distributed by Phenomenex.The column has already been used in several cases for the separation of different enantiomers for a specific purpose, such as besifloxacin hydrochloride 17 , carvone 18 , citalopram 19 , fexofenadine 20 , omeprazole and 5-hydroxyomeprazole 21 .However, a comprehensive study investigating the utilization of this commercialized column for diverse components has not yet been published.Further clarification is needed regarding the performance characteristics of the column.
One of the key advantages of using CD-based CSPs in chiral separations is the ability to utilize multimodal elution, which enables the resolution of a wide range of stereoisomeric compounds through the combination of different elution modes, including reversed-phase, normal-phase, and polar organic-mode 22 .Several studies recommend using reversed-phase or normal-phase mode with CD-based CSPs 23,24 .However, phenylcarbamated CDs, similar to polysaccharide selectors with phenylcarbamate groups, can also exhibit good performance in polar organic mode [25][26][27] .Despite this, inferior performance was reported with in-house prepared phenylcarbamated-β-CD column in polar organic mode according to previous literature 28 .In polar organic mode, the mobile phase consists solely of polar organic solvents, such as pure alcohols (methanol (MeOH), ethanol (EtOH), and 2-propanol (IPA)) or pure acetonitrile (ACN), individually or in combination.To prevent the formation of ionic species, an acid (e.g., acetic acid (Acac)) or a base (e.g., diethylamine (DEA)) can be introduced to the mobile phase, depending on the analyte's acid-base properties.Polar organic mode offers several advantages, including shorter run times, high efficiency, easier mass spectrometric coupling, and enhanced solubility of the analytes in the mobile phase.Additionally, polar organic mode is suitable for both analytical and preparative applications [29][30][31][32][33] .The present study aimed to investigate the enantioseparation performance of the Chiral CD-Ph column in polar organic mode for a range of structurally diverse compounds and to clarify the mechanism of chiral recognition through a molecular docking study.

Materials.
Racemic thalidomide, pomalidomide, (S)-pomalidomide, (R)-rabeprazole, lenalidomide, (R)and (S)-apremilast were purchased from Beijing Mesochem Technology Co. Ltd. (Beijing, China).(R)-thalidomide, (S)-enantiomer-enriched lenalidomide 34 , guaifenesin, warfarin, ibuprofen, 1-aminoindane, amphetamine sulfate, methamphetamine hydrochloride, bamethan hemisulfate, terbutaline sulfate, propranolol hydrochloride, metoprolol tartrate, norepinephrine hydrochloride, and omeprazole were ordered from Sigma-Aldrich, Hungary (Budapest, Hungary).Stiripentol, naproxen, naringenin, hesperetin, bisoprolol hemifumarate, ofloxacin was from TCI Chemicals (Tokyo, Japan).Rabeprazole sodium was United States Pharmacopoeia (USP) reference standard (Rockville, USA).Iminoflavan derivatives were synthesized as previously described 35,36 .The structure of investigated compounds is depicted in Fig. 1, while certain physico-chemical parameters are summarized in Table 1.Gradient grade MeOH and ACN as well as Acac and DEA were purchased from Thomasker Finechemicals Ltd. (Budapest, Hungary).Chiral CD-Ph column (250 × 4.6 mm, particle size 5 μm), Nucleodex  The optimized phenylcarbamate-β-CD structure has been used as a host in docking calculations.The 3D structures of thalidomide, pomalidomide, lenalidomide, stiripentol, hesperetin, metoprolol, norepinephrine, and rabeprazole enantiomers were prepared and the corresponding PDBQT files were created which includes information on the torsional degree of freedom, atomic partial charges and atom types, polar hydrogens, but not hydrogens bonded to carbon atoms.These files were used as inputs docked to the prepared phenylcarbamateβ-CD by using AutoDock Vina 38 and the binding affinity of the complexes was computed.A maximum of 9 binding modes were considered for each enantiomer and the exhaustiveness of the calculations was set to 8. The best binding modes (lowest in terms of energy) were selected, and the corresponding binding affinities (E A ) were compared.The grid box was 30 × 30 × 30 Å 3 in each case.

Results and discussion
Screening in polar organic mode.To evaluate the performance of the Chiral CD Ph column in polar organic mode, 25 compounds were selected with diverse chemical structures and acid-base properties.ACN or MeOH was used as eluent, supplemented with 0.1% Acac (v/v) for acidic compounds and 0.1% DEA (v/v) for basic compounds.For amphoteric compounds, both modifiers (Acac and DEA) were used in equal proportions (0.1-0.1% (v/v)).The obtained results are summarized in Table 2, and representative chromatograms are presented in Fig. 3.
The results in Table 2 show that the Chiral CD-Ph column exhibits enantiorecognition towards 20 out of 25 racemic compounds tested.Notably, twelve compounds were separated with baseline resolution (R s ≥ 1.5) www.nature.com/scientificreports/without requiring any method development.For five compounds, including iminoflavan derivatives, hesperetin, guaifenesin, and stiripentol, the resolution exceeded three.However, for five compounds, including thalidomide, warfarin, naringenin, 1-aminoindane, and omeprazole, the obtained enantioresolution values were lower than 0.7.
To ensure comparability with other commercially available columns, we repeated the same measurements on three additional chiral columns.Among these, two are CD-based CSPs, commercially available in the market (Astec Cyclobond I 2000, containing β-CD as chiral selector, and Nucleodex β-PM, containing permethylated β-CD as chiral selector).Lux Amylose-1 was the third column, a polysaccharide-type CSP, containing amylose tris(3,5-dimethylphenyl)-carbamate as chiral selector.The latter being generally recognized as one of the most effective chiral columns available on the market.Enantioseparation results obtained on these three columns are summarized in Supplementary Tables 1, 2, and 3, respectively.The results clearly demonstrate that in polar organic mode, the performance of the Chiral CD-Ph column is significantly superior compared to other CDbased CSPs.Using the Astec Cyclobond I 2000 and the Nucleodex β-PM columns, only a total of nine enantioselective separations were observed, but these were generally characterized by poor resolution values.The highest enantioresolution value achieved on these columns in polar organic mode was a mere 1.2 (warfarin on Astec Cyclobond I 2000, using ACN-based mobile phase).The obtained results indicate that CD-based columns, where the structure of the chiral selector contains no aromatic groups, exhibit subpar performance in polar organic mode.Among these columns, ACN tends to yield better results than MeOH due to the absence of aromatic rings that could potentially form π-π interactions in alcohol-based mobile phases.In contrast, not surprisingly, out of all four columns, the Lux Amylose-1 demonstrated the most remarkable performance, frequently yielding very high resolution values.Nonetheless, the performance of the Lux Amylose-1 column and the Chiral CD-pH column is comparable.A comparison of enantioseparation systems using three different characterization approaches is illustrated in Supplementary Fig. 1.Our results clearly indicate that the Chiral CD-Ph column is a promising choice for enantiomer separation in polar organic mode.
It is also noteworthy that three out of the five compounds where enantioseparation was not observed, all present a carboxylic group in their structure (naproxen, ibuprofen, and ofloxacin).This observation is interesting and may suggest a correlation between the presence of a carboxylic group and the lack of enantiorecognition on the Chiral CD-Ph column in polar organic mode.Upon comparing the performances of the separation systems Table 2. Chromatographic data during the screening method on Chiral CD-Ph column (k 1 -retention factor of the first eluting enantiomer, k 2 -retention factor of the second eluting enantiomer, α-selectivity, R s -resolution).*For neutral compounds (1-7) pure eluent, for acidic compounds (8-12) the eluent was modified with 0.1% Acac (v/v), for basic compounds (13-21) the eluent was modified with 0.1% DEA (v/v), for amphoteric compounds (22-25) the eluent was modified with the mixture of DEA-Acac 0.1-0.1% (v/v).www.nature.com/scientificreports/using ACN and MeOH-based eluents, it can be easily observed that MeOH appears to be more effective than ACN for these analytes, particularly for acidic compounds.However, there are some analytes (thalidomide, pomalidomide, omeprazole) that show better results with ACN-based mobile phases.Thus, it is recommended to investigate both eluents for screening.Based on the results from the sample set employed, no clear pattern can be discerned regarding retention factors or enantioselectivity.From the data, in most cases, the retention factor is not lower in ACN than in MeOH.Additionally, we can find examples where a low retention factor is paired with high enantioresolution value (e.g., stiripentol), while in several cases, despite the high retention factor, enantioseparation was not observed (e.g., apremilast, omeprazole).The compounds selected exhibit structural diversity, but certain groups of compounds have similar structures, such as thalidomide analogs (thalidomide, pomalidomide, lenalidomide, apremilast), flavanones (naringenin and hesperetin), β-blockers (propranolol, bisoprolol, metoprolol), proton-pump inhibitors (omeprazole and rabeprazole), iminoflavan derivatives, and the amphetamine-methamphetamine pair.The chromatographic behavior was compared within these groups.It was observed that the chromatographic behavior of amphetamine and methamphetamine, as well as the iminoflavan derivatives, is similar despite the differences in acid-base properties of iminoflavans.Thalidomide and pomalidomide exhibited related results, both compounds being separated only in ACN.However, lenalidomide, which has a minor structural difference of an oxo group compared to pomalidomide, also displayed enantioseparation in MeOH.In contrast, apremilast, which has significant structural variations from its parent compound, was not successfully separated.Naringenin and hesperetin exhibited similar chromatographic behavior and could only be separated using a MeOH-based eluent.However, while the two compounds are structurally similar, naringenin had a relatively low resolution of 0.6 compared to the notably higher resolution of 3.0 for hesperetin.Among the proton pump inhibitors, rabeprazole exhibits a higher resolution when separated in methanolic eluent compared to omeprazole.It is also worth noting that while omeprazole can be separated in ACN, rabeprazole cannot.

MeOH
In conclusion, it can be inferred that predicting enantiomer separation on the Chiral-CD-Ph column is challenging as even slight structural modifications in the analytes can significantly affect the chromatographic outcome.However, this is also generally true for other chiral columns.Exploring the interactions between the selector and the selectand can aid in more accurately predicting the success of chiral separations 39,40 .
The influence of various chromatographic parameters, including eluent composition, eluent modifier, temperature, and flow rate on the enantioseparation performance was also investigated on the Chiral CD-Ph column, using a restricted set of analytes.www.nature.com/scientificreports/ The impact of various chromatographic parameters on chiral separation.For a more in-depth analysis, we have chosen to investigate a range of compounds including thalidomide derivatives and stiripentol as representative of neutral compounds, naringenin and hesperetin from acidic compounds, bisoprolol, metoprolol, and basic iminoflavan form basic compounds, and norepinephrine and rabeprazole from amphoteric compounds.
The effects of various eluent additives were thoroughly investigated, including the influence of MeOH-ACN mixtures on enantioseparation.Additionally, we evaluated the impact of changing flow rate and temperature on chromatographic performance.To study the impact of acidic additive, different percentages of Acac (0.05%, 0.1%, 0.15% (v/v)) were added to MeOH, and for the basic additive different percentages of DEA (0.05%, 0.1%, 0.15% (v/v)) were utilized.For neutral and amphoteric compounds, the effects of both modifiers were examined separately and in combination.In the case of neutral compounds, the addition of modifiers did not significantly influence the separation performance.However, for basic compounds such as bisoprolol and metoprolol, enantioseparation was only achieved in the presence of a basic modifier with the best resolution observed using 0.15% DEA (v/v).Furthermore, the basic iminoflavan compound was found to be baseline separated even in neat MeOH, which can possibly be attributed to its highly lipophilic nature and the large distance between the basic nitrogen and the chiral center (Supplementary Table 4).The addition of Acac in the case of naringenin and hesperetin improved resolution up to 0.1% Acac (v/v) (Supplementary Table 5).Our investigation of the amphoteric compound rabeprazole, revealed that enantioseparation is possible in pure MeOH, while this was not the case for norepinephrine.When using DEA as a modifier in MeOH, norepinephrine enantiomers were baseline separated, however, upon increasing the DEA concentration of the mobile phase, enantioresolution significantly worsened.Using Acac as modifier alone resulted in the lack of noticeable enantioseparation for amphoteric compounds.The best resolution was obtained using a combination of both 0.1% Acac (v/v) and 0.1% DEA (v/v) (Supplementary Table 6).
The impact of binary mobile phases containing varying percentages of MeOH and ACN on resolution and retention factors was also evaluated.Figure 4 shows the effect of MeOH content in ACN on the retention and resolution of stiripentol as an example (similar U-shaped retention curves were also observed for the other set of investigated analytes).It is interesting to observe that the lowest retention and resolution values were observed at 50% ACN content in MeOH (v/v).From this local minimum increasing the MeOH or ACN content led to an increase in enantioresolution.The highest resolution can be observed by applying 100% ACN or 100% MeOH.This U-shape curve indicates that the mechanism of separation is different in MeOH compared to ACN.This is explained by the fact that the probability of H-bond formation between the stiripentol enantiomers and the chiral selector is higher in ACN than in MeOH, where most probably π-π interactions are dominant.Similar observations have already been made on polysaccharide-and cyclofructan-based columns in polar organic mode 26,41,42 .Additionally, it should be noted that the polarity and the dielectric constant of the two solvents also play a role in determining the strength of the interactions between the enantiomers and the chiral selector, which can further influence separation performance.
Variations in flow rate did not elicit any unexpected results.As the flow rate increased, retention decreased while the selectivity remained unchanged.In general, a higher flow rate resulted in a slight reduction in resolution, but the retention times can be reduced, and the peak shapes can be improved by applying higher flow rates (Supplementary Fig. 2).www.nature.com/scientificreports/Thermodynamic analyses.The separation temperature is a critical factor that can significantly impact the selectivity, resolution, retention and EEO of enantiomers during HPLC enantioseparation.Thermodynamic calculations are highly significant and commonly employed methods for establishing a chiral recognition mechanism using direct HPLC.The aim of the thermodynamic study in our research was to provide insights into the interaction energies and driving forces behind chiral recognition on the Chiral CD-Ph column in polar organic mode.By employing thermodynamic calculations, we aimed to elucidate the underlying mechanisms that contribute to the enantioseparation efficiency observed in our experimental results 5,43,44 .
To investigate the energetic interactions between the analytes and the stationary phase a classical van 't Hoff approach was used, assuming only enantioselective interaction sites on the stationary phase.Van 't Hoff plots were constructed by plotting the natural logarithm of the retention factor as a function of the inverse of the absolute temperature over a temperature range of 10-40 °C (applying 10 °C, 20 °C, 25 °C, 30 °C, and 40 °C).
where R stands for the universal gas constant, T is the temperature in Kelvin, k is the retention factor of the enantiomers.ΔH° denotes the standard enthalpy, while ΔS° is the standard entropy change of transfer of the solute from the mobile phase to the stationary phase, and is the phase ratio of the Chiral CD-Ph column.If ΔH° is constant in the selected temperature range, a linear relationship is obtained between lnk and 1/T, with a slope of − ΔH°/R and an intercept of ΔS°/R + ln .Since the value of the phase ratio is seldom known, ΔS°* values (ΔS°* = ΔS° + Rln ) are often used, to compensate the uncertainty in 5,45 .
Similarly, the difference in change of standard enthalpy Δ(ΔH°) and standard entropy Δ(ΔS°) for the two enantiomers moving from the mobile phase to the stationary phase were also calculated according to modified van 't Hoff equation: Isoenantioselective temperatures (T iso ) were also calculated as the ratio between ΔΔH° and ΔΔS°.T iso represents the temperature at which the enthalpy contribution is compensated by the entropic term, indicating that Gibbs free energy (ΔΔG°) equals zero.At T iso , the two enantiomers co-elute, and no separation is achieved.ΔΔG° provides information on the strength of binding between selector and selectant, more negative values indicating a more efficient binding.Consequently, a value of Δ(ΔG°) = 0 means that there is no difference between the binding strength of the enantiomers, therefore, at T iso the two enantiomers co-elute, and no separation is achieved.
The obtained thermodynamic data are summarized in Table 3.By analyzing the obtained data, it is possible to gain valuable insights into the factors influencing enantioselectivity on the Chiral CD-Ph column and to optimize the separation conditions for enhanced selectivity and resolution.
Based on the data obtained at various temperatures, it is evident that the retention factor consistently decreases as separation temperature increases, irrespective of the analyte under examination.As expected, selectivity decreased with increasing temperature except for naringenin in MeOH-based mobile phase where selectivity slightly increased with increasing temperature.Representative chromatograms for the thermodynamic study are depicted in Fig. 5.
Linear regression analyses were performed for all cases, comparing lnk vs. 1/T and lnα vs. 1/T, indicating that the binding mode and conformation of the selector remained unchanged within the temperature range investigated.Supplementary Fig. 3 displays fitted van 't Hoff plots, illustrating the relationship between lnα and 1/T.The Δ(ΔH°) values provide insights into the relative ease of analyte transfer from the mobile phase to the stationary phase.These values ranged from − 22.49 kJ/mol (stiripentol in ACN) to 0.99 kJ/mol (naringenin in MeOH).The significantly low Δ(ΔH°) value observed for stiripentol in ACN suggests that the second eluted enantiomer forms an additional strong secondary interaction, such as a hydrogen bond with the selector, indicating relatively strong levels of chiral recognition.Negative Δ(ΔH°) values for enantiomeric pairs are accompanied by negative Δ(ΔS°) values.A negative Δ(ΔS°) reflects an increase in order and/or a reduction in degrees of freedom during enantioselective interactions between the selector and selectand.Just as the highest negative Δ(ΔH°) value is found for stiripentol in ACN, the highest Δ(ΔS°) value is also found for stiripentol in ACN, which is − 66.10 J/ molK.The relative contributions of enthalpic and entropic terms to the free energy of adsorption can be elucidated by examining the enthalpy/entropy ratio, denoted as Q.By comparing the Q values, it was observed that the enantioseparation was predominantly controlled by enthalpy in most cases, as indicated by Q values larger than 1.Interestingly, for naringenin, enantioseparation controlled by entropy (Q < 1) were observed.This finding suggests that the binding of naringenin to the chiral selector may exhibit distinct characteristics compared to other analytes, potentially involving different interactions or mechanisms.Temperature-dependent EEO reversal was not observed.When the column temperature was T iso , the isomers would be eluted simultaneously.However, if the column temperature was higher than T iso , the EEO would be reversed.In all investigated cases the T iso value is out of the range of use.
The results indicate that the calculated binding energies are within the range of − 4.5 to − 7.2 kcal/mol.The negative values suggest that the formation of the host-guest complexes is a spontaneous process driven by enthalpy.It is important to emphasize that in all cases where the enantiomer elution order was determined, the molecular docking results and experimental findings are in good agreement with each other (Table 4).By analyzing the optimized structure of phenylcarbamate-β-CD it was found that the cavity of the CD host is blocked by the phenylcarbamate groups.Thus, surface binding occurred, and the analytes established interactions with the surface of the phenylcarbamate-β-CD (Fig. 6).For the clear presentation of the accessibility of the cavity in native CDs and phenylcarbamate-β-CD, the molecular surface of the 3D structure was depicted in Supplementary Fig. 4. Based on our this, attempting to "force" the analyte into the cyclodextrin cavity in case of the phenylcarbamoylated-CDs is an erroneous approach.
The difference in binding energy (ΔE A ) between a pair of enantiomers can be considered a measure of chiral recognition, and a non-zero ΔE A is a criterion for enantiomer separation.Although ΔE A in some cases is only 0.1 kcal/mol (e.g.rabeprazole), the two enantiomers are still baseline separated.In theory, the ΔE A values should correspond with the selectivity magnitudes obtained from the enantioseparation conditions, but this is not always the case.These discrepancies arise from the limitations of docking which prevents taking into account all experimental factors such as the impact of additives, solvation effects, unreacted silica particles in the columns, and others.Despite the listed limitations, docking calculations were able to predict the observed enantiomer elution order in four cases.Additionally, four additional analytes were also studied, but for those, experimental data regarding the enantiomer elution order were not available due to the absence of enantiomerically pure standards.Various types of interactions were established between the analytes and the host including stronger hydrogen bonds and hydrophobic contacts, and π-π interactions (Fig. 6).To validate the docking results, singlepoint energy calculations were conducted on specific complexes (such as pomalidomide-Chiral CD-Ph).The calculations utilized the APFD hybrid density functional, and the influence of solvent effects was considered by applying the conductor-like polarizable continuum model (CPCM) The outcomes revealed that the enantiomer energy hierarchy remained unchanged; however, there was an observed increase in relative energy (see Supplementary Table 7).

Conclusions
It was found that the Chiral CD-Ph column demonstrated good performance in separating a diverse range of 25 racemic compounds with varying chemical structures and acid-base properties in polar organic mode.The high success rate in separating the tested compounds highlights the effectiveness of the Chiral CD-Ph column in polar organic mode, especially when using MeOH as the mobile phase.We investigated the influence of mobile phase compositions, additives, flow rate, and temperature on enantioseparation in detail and found that the appropriate mobile phase is crucial for optimal results.Based on our results, MeOH-based mobile phases lead to better outcomes than ACN-based ones.Thermodynamic analysis revealed mainly enthalpy-controlled enantioseparations.We also investigated enantiorecognition at the molecular level using molecular docking and found that molecules bind to the surface of the CD.Although there were some limitations to our study, such as the lack of enantiomerically pure standards for some analytes, docking calculations were still able to predict the observed EEO in several cases.Overall, our work demonstrates that the Chiral CD-Ph column is a reliable and versatile tool for chiral analysis in polar organic mode.

Figure 1 .
Figure 1.Chemical structures of the studied test compounds (chiral centers marked by an asterisk).

Figure 2 .
Figure 2. 2D (left) and 3D (middle) representations of phenylcarbamate-β-cyclodextrin.R = phenylcarbamate.The molecular surface of the 3D structure (right) is also presented, where red indicates exposed regions (to the solvent), green signifies hydrophobic groups, and magenta represents polar groups.

Figure 4 .
Figure 4. Plots of the retention and resolution factors as a function of the MeOH content on Chiral CD-Ph column.(Chromatographic conditions: flow rate, 0.5 mL/min; column temperature, 20 °C).

Figure 6 .
Figure 6.Interactions between the best binding modes of thalidomide enantiomers and phenylcarbamateβ-CD.The molecular surface of the 3D structure (right) is also presented, where red indicates exposed regions (to the solvent), green signifies hydrophobic groups, and magenta represents polar groups.

Table 4 .
Binding affinities (E A ) of the R and S enantiomers of the studied analytes towards phenylcarbamateβ-CD.The types of main interactions are also listed along with the observed experimental enantiomer elution order where it is available.n.d.no data.

Table 1 .
Calculated physico-chemical parameters of the studied compounds.