Locating Regions of Maximum Chiral Discrimination: A Computational Study of Enantioselection on a Popular Chiral Stationary Phase Used in Chromatography

Molecular Simulation of the Separation of Some Amino Acid Enantiomers by β-Cyclodextrin in Gas-Phase

The complexes formed by β-cyclodextrin and some amino acids (alanine, valine, leucine, and isoleucine) in vacuo are studied by molecular mechanics and dynamics simulations. These methods have been improved with respect to our previous studies with amino acids, regarding the determination of molecular structures or initial enantiomer dispositions in the molecular dynamics trajectories. The greatest contribution to the interaction energy is from the van der Waals term, although the discrimination between enantiomers is due mainly to the electrostatic contribution. The lowest energy structures of the complexes obtained from molecular mechanics are inclusion complexes in which the carboxylic end of amino acids is pointing toward the narrow (D-) or wide rim (L-) of β-cyclodextrin. The position probability density provided by molecular dynamics also confirms inclusion complex formation, because the guests spend most time inside the cavity of β-cyclodextrin along its axis, with the carboxylic end pointing toward the narrow rim. The L-amino acids are the first eluted enantiomers in all cases and chiral discrimination increases with the size of guests, except leucine, which has the lowest capacity to discriminate. During the simulation, Ala and Val remain in weakly enantioselective regions, while Leu and Ile stay in zones with great chiral selectivity.

Noncovalent Complexes of Cyclodextrin with Small Organic Molecules: Applications and Insights into Host-Guest Interactions in the Gas Phase and Condensed Phase

Cyclodextrins (CDs) have drawn a lot of attention from the scientific communities as a model system for host–guest chemistry and also due to its variety of applications in the pharmaceutical, cosmetic, food, textile, separation science, and essential oil industries. The formation of the inclusion complexes enables these applications in the condensed phases, which have been confirmed by nuclear magnetic resonance (NMR) spectroscopy, X-ray crystallography, and other methodologies. The advent of soft ionization techniques that can transfer the solution-phase noncovalent complexes to the gas phase has allowed for extensive examination of these complexes and provides valuable insight into the principles governing the formation of gaseous noncovalent complexes. As for the CDs’ host–guest chemistry in the gas phase, there has been a controversial issue as to whether noncovalent complexes are inclusion conformers reflecting the solution-phase structure of the complex or not. In this review, the basic principles governing CD’s host–guest complex formation will be described. Applications and structures of CDs in the condensed phases will also be presented. More importantly, the experimental and theoretical evidence supporting the two opposing views for the CD–guest structures in the gas phase will be intensively reviewed. These include data obtained via mass spectrometry, ion mobility measurements, infrared multiphoton dissociation (IRMPD) spectroscopy, and density functional theory (DFT) calculations.

Molecular Simulation of the Separation of Isoleucine Enantiomers by β-Cyclodextrin

Molecular mechanics and dynamics simulations were carried out to study the capacity of isoleucine enantiomers to form inclusion complexes with β⁻cyclodextrin, and to be discriminated by this chiral compound, in vacuo and with different solvents. Solvents were characterized not only by the value of dielectric constant ε in the Coulombic interaction energy, but also by the neutral and zwitterion configurations of isoleucine. Whereas the discrimination between the enantiomers for ε ≤ 2 is due to the electrostatic contribution, these differences are mainly due to the Lennard-Jones potential for ε > 2. The most enantioselective regions are located near the cavity walls, independently of the solvent. D-Ile is more stable than L-Ile in broader regions in vacuo, but L-Ile presents more stable locations with water. Isoleucine can form inclusion complexes with β⁻cyclodextrin in vacuo and with different solvents. Two probable configurations are deduced from the molecular dynamics simulation, in which the guest is always inside the cavity and with the carboxylic end of the amino acid oriented towards either rim of β⁻CD. In the simulation, the enantiomers preferentially occupy regions with greater chiral discrimination. The first eluted enantiomer in vacuo and with different solvents is L-Ile, independently of the solvent polarity.

Theoretical Study of the β-Cyclodextrin Inclusion Complex Formation of Eugenol in Water

The interaction between eugenol and β-cyclodextrin in the presence of water is studied by molecular mechanics and dynamics simulations. A force field model is used in molecular mechanics to determine the interaction energy and the complex configuration at the absolute minimum. The van der Waals term is the main contribution to the total energy, and so directly determines the configuration of the inclusion complex. The formation of inclusion complexes is simulated by molecular dynamics, in which their configurations are deduced from the position probability density that represents the preferred location and orientation of the guest in the simulation. When eugenol approaches from the rims of β-cyclodextrin, it tends to enter the cavity, remain inside for a short period and then exit from it. The guest tends to include the phenyl ring inside the cavity in the most probable configurations. Two inclusion complex configurations are proposed, each with the hydroxyl and methoxyl groups pointing towards one different rim of β-cyclodextrin. The initial guest orientation is the main factor determining these configurations. The model presented in this study reproduces the experimental findings on inclusion complex formation and proposes two possible complex configurations, one previously suggested by different authors.

Molecular Dynamics Study of the Factors Influencing the β-Cyclodextrin Inclusion Complex Formation of the Isomers of Linear Molecules

The influence of the size, composition, and atomic distribution of linear guests on β-cyclodextrin inclusion complex formation is clarified by means of a molecular dynamics simulation at constant temperature. The intermolecular energy is modelled by a Lennard-Jones potential, where the molecular composition is represented by various parameters and by a continuum description of the guest and cavity walls. It is concluded that the parameters related to the atomic size require minimum values for the confinement of linear molecules inside the cavity. The isomer with optimal affinity for β-cyclodextrin as predicted by the free energy presents an asymmetrical molecular structure, and the position probability density shows that the isomer tends to insert the portion with largest atoms into the cavity, although the preferential binding site of the guest is not always located in regions of the host with maximum discriminatory power.

Self-assembly of perfunctionalized beta-cyclodextrins on a quartz crystal microbalance for real-time chiral recognition

This study reports the development of novel chiral sensors based on the self-assembly of perfunctionalized beta-cyclodextrins (beta-CD) on a quartz crystal microbalance transducer for real time chiral recognition. Ten chiral sensors immobilized with mercaptyl perfunctionalized beta-cyclodextrins were found to exhibit promising enantioselectivity in the gas phase. Well-defined sizes of molecular cavities of the modified beta-CDs associated with cooperative weak interactions with the host molecules afforded enhanced chiral discrimination. This study contributes to the realization of novel chiral sensors applicable for real-time recognition and analysis of enantiomeric alcohols and lactates.

Nitroxide Radicals as Probes for Exploring the Binding Properties of the Cucurbit[7]uril Host

EPR spectroscopy was used for the first time to explore the binding properties of cucurbit[7]uril (CB7), a representative member of the cucurbituril family. Evidence for the formation of a complex between nitroxide radicals and the host system in an aqueous solution was provided by large changes in the nitrogen hyperfine splitting, attributed to the different polar environments experienced by the included radical. In the presence of alkali cations, the EPR spectra of benzyl tert-butyl nitroxide were characterised by new signals attributed to the radical hosted in the CB7 cavity in which one metal cation is in close contact with the nitroxidic oxygen. The formation of the coordination complex results in a substantial increase in the electron spin density on the nitrogen in inverse order with respect to the size of the cation owing to increased localisation of negative charge on the oxygen atom from bonding to the alkali cation. The EPR spectra showed selective line-broadening effects as a result of metal exchange between bulk water and the coordination complex. Analysis of the EPR linewidth variations allowed us to measure the corresponding kinetic rate constants for the first time. NMR spectroscopy showed that this behaviour is not peculiar to nitroxides but is also exhibited by the related carbonyl compounds. These data allowed us to quantify the template effect and to reach the conclusion that, in the presence of a guest having a coordinating lone pair, the formation of ternary metal-guest-CB complexes must be taken into account when discussing the complexation behaviour of cucurbituril derivatives in the presence of salts.

Adsorption behavior of the (+/-)-Tröger's base enantiomers in the phase system of a silica-based packing coated with amylose tri(3,5-dimethyl carbamate) and 2-propanol and molecular modeling interpretation

The binary adsorption isotherms of the enantiomers of Tröger's base in the phase system made of Chiral Technologies ChiralPak AD [a silica-based packing coated with amylose tri(3,5-dimethyl carbamate)] as the chiral stationary phase (CSP) and 2-propanol as the mobile phase were measured by the perturbation method. The more retained enantiomer exhibits a S-shaped adsorption isotherm with a clear inflection point, the concentration of the less retained enantiomer having practically no competitive influence on this isotherm: In the entire range of concentrations studied, dq2/dC1 approximately 0. By contrast, the less retained enantiomer has a Langmuir adsorption isotherm when pure. At constant mobile phase concentrations, however, its equilibrium concentration in the adsorbed phase increases with increasing concentration of the more retained enantiomer and dq1/dC2 > 0. This cooperative adsorption behavior, opposed to the classical competitive behavior, is exceedingly rare but was clearly demonstrated in this case. Two adsorption isotherm equations that account for these physical observations were derived. They are based on the formation of an adsorbed multi-layer, as suggested by the isotherm data. The excellent agreement between the experimental overloaded elution profiles of binary mixtures and the profiles calculated with the equilibrium-dispersive model validates this binary isotherm model. The adsorption energies calculated by molecular mechanics (MM) and by molecular dynamics (MD) indicate that the chiral recognition arising from the different interactions between the functional groups of the CSP and the molecules of the Tröger's base enantiomers are mainly driven by their Van der Waals interactions. The MD data suggest that the interactions of the (-)-Tröger's base with the CSP are more favored by 8+/-(5) kJ/mol than those of (+)-Tröger's base. This difference seems to be a contributing factor to the increased retention of the - enantiomer on this chromatographic system. The modeling of the data also indicates that both enantiomers can form high stoichiometry complexes while binding onto the stationary phase, in agreement with the results of the equilibrium isotherm studies.

Synthesis and application of mono-2A-azido-2A-deoxyperphenylcarbamoylated beta-cyclodextrin and mono-2A-azido-2A-deoxyperacetylated beta-cyclodextrin as chiral stationary phases for high-performance liquid chromatography

Two novel chiral stationary phases (CSPs) were prepared based upon the regioselective immobilizations of beta-cyclodextrin (beta-CD) at its C2 position to the silica support. The mono-2A-azido-2A-deoxyperphenylcarbamoylated beta-cyclodextrin and mono-2A-azido-2A-deoxyperacetylated beta-cyclodextrin were synthesized by selective tosylation and azidolysis followed by perfunctionalisation. The derivatised cyclodextrins were then immobilized onto the aminised silica gel via the Staudinger reaction to provide new chiral stationary phases. Their application to high-performance liquid chromatography (HPLC) enantioseparation of racemic compounds was demonstrated using beta-adrenergic blockers, flavonone compounds, benzodiazepinones, antihistamines and weakly protolytic compounds, of which good separations were achieved for some racemic compounds, for instance, bendroflumethiazide (Rs 6.26), oxazepam (Rs 5.99), temazepam (Rs 2.85) and althiazide (Rs 1.13) when compared with the corresponding CSPs where the beta-CD molecule was regioselectively immobilized at the C6 position. The enantiodiscriminatory properties of these CSPs were found to be affected by the orientation of the CD cavity under reversed-phase conditions, and also by the derivitising groups of the CD. The HPLC results inferred that the mono-6A-azido-6A-deoxyperphenylcarbamoylated CD CSP (CD bonded at C6 position to silica) exhibited slightly better chiral recognition ability than mono-2A-azido-2A-deoxyperphenylcarbamoylated CD CSP under the normal-phase and reversed-phase modes on the separation of 31 different racemic compounds and drugs. On the contrary, higher chiral recognition abilities were observed on the mono-2(A)-azido-2A-deoxyperacetylated CD CSP compared to mono-6A-azido-6A-deoxyperacetylated CD CSP.

Theoretical reassessment of Whelk-O1 as an enantioselective receptor for 1-(4-halogeno-phenyl)-1-ethylamine derivatives

A combination of molecular mechanics and first principles calculations was used to explore the enantioselectivity of receptors, taking into account experimental data from the CHIRBASE database. Interactions between the Whelk-O1 HPLC chiral stationary phase with the complete series of 1-(4-halogeno-phenyl)-1-ethylamine derivative racemates were studied. The objective was to extract information from the interactions between the chiral Whelk-O1 stationary phase and the enantiomers, hence probing the origin of the enantioselective behavior. Calculations correctly reproduce the elution orders and reasonably describe the experimental enantioselectivities and retention factors. Different binding modes were observed for the first eluted enantiomer complexes, whereas the second eluted show only one prevalent diastereomeric binding fashion. Natural bond orbital (NBO) analysis was used on the global minima bound-complexes to quantify donor-acceptor interactions among chiral stationary phase and ligand moieties. Intermolecular hydrogen bonding was found to be the essential energetic interaction for all systems studied. CH-pi, aromatic stacking and various charge transfer interactions were found to be smaller in magnitude but still important for the global enantioselective behavior. The three-point interaction model is discussed, pointing out the difficulty of its application for the qualitative prediction of elution orders (absolute configurations).

Molecular dynamics (MD) simulations for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (β-cyclodextrin, β-CD) based on the MM–PBSA (molecular mechanics–Poisson–Boltzmann surface area) approach

Molecular dynamics (MD) simulations were performed for the prediction of chiral discrimination of N-acetylphenylalanine enantiomers by cyclomaltoheptaose (beta-cyclodextrin, beta-CD). Binding free energies and various conformational properties were obtained using by the MM-PBSA (molecular mechanics Poisson-Boltzmann/surface area) approach. The calculated relative difference (DeltaDeltabinding) of binding free energy was in fine agreement with the experimentally determined value. The difference of rotameric distributions of guest N-acetylphenylalanine enantiomers complexed with the host, beta-CD, was observed after the conformational analyses, suggesting that the conformational changes of guest captured within host cavity would be a decisive factor for enantiodifferentiation at a molecular level.

Probability rule for chiral recognition

Molecular Chirality is of central interest in biological studies because enantiomeric compounds, while indistinguishable by most inanimate systems, show profoundly different properties in biochemical environments. Enantioselective separation methods, based on the differential recognition of two optical isomers by a chiral selector, have been amply documented. Also, great effort has been directed towards a theoretical understanding of the fundamental mechanisms underlying the chiral recognition process. Here we report a comprehensive data examination of enantio separation measurements for over 72000 chiral selector-select and pairs from the chiral selection compendium CHIRBASE. The distribution of alpha = k'(D)/k'(L) values was found to follow a power law, equivalent to an exponential decay for chiral differential free energies. This observation is experimentally relevant in terms of the number of different individual or combinatorial selectors that need to be screened in order to observe alpha values higher than a preset minimum. A string model for enantiorecognition (SMED) formalism is proposed to account for this observation on the basis of an extended Ogston three-point interaction model. Partially overlapping molecular interaction domains are analyzed in terms of a string complementarity model for ligand-receptor complementarity. The results suggest that chiral selection statistics may be interpreted in terms of more general concepts related to biomolecular recognition.

Chiral recognition: design and preparation of chiral stationary phases using selectors derived from ugi multicomponent condensation reactions and a combinatorial approach

Combinatorial approaches together with high-throughput screening have been used to develop highly selective stationary phases for chiral recognition. Libraries of potential chiral selectors have been prepared by the Ugi multicomponent condensation reactions and screened for their enantioselectivity using the reciprocal approach involving a chiral stationary phase with immobilized model target compound N-(3,5-dinitrobenzoyl)-alpha-l-leucine. The best candidates were identified from the library of phenyl amides of 2-oxo-azetidineacetic acid derivatives. This screening also enabled specification of the functionalities of the selector desired to achieve the highest level of chiral recognition. The substituents of the phenyl ring adjacent to the chiral center of the selector candidates exhibited the most profound effect on the chiral recognition. The best candidate was then synthesized on a larger scale, resolved into single enantiomers using preparative enantioselective HPLC, and attached to porous poly(2-hydroxyethyl methacrylate-co-ethylene dimethacrylate) beads via an ester linkage to afford the desired stationary phase. Selectivities alpha as high as 3.2 were found for the separation of a variety of amino acid derivatives.

Macrocycle conformation and self-inclusion phenomena in octakis(3-O-butanoyl-2,6-di-O-pentyl)-gamma-cyclodextrin (Lipodex E) by NMR spectroscopy and molecular dynamics

Octakis(3-O-butanoyl-2,6-di-O-pentyl)-gamma-cyclodextrin (Lipodex E) is a lipophilic chiral selector successfully used for the enantioselective gas chromatographic separation of a multitude of racemic analytes. NMR data (13C chemical shifts, 3J(HH), rotating frame NOEs (ROEs)) and molecular dynamics (MD) simulations point out that the macrocycle is distorted with respect to the canonical truncated-cone shape of native cyclodextrins, although C(8) symmetry is retained on the NMR timescale. ROE data and MD trajectories provide evidence for self-inclusion of one 6-O-pentyl pendant chain within the cavity of Lipodex E. The interpretation of long-range and low-intensity ROEs is supported by the calculation of average internuclear distances by using the radial distribution function (RDF) calculated from MD trajectories. MD simulations are eventually used to compare the flexibility of the macrocycle of Lipodex E with that of native gammaCD.

NMR studies of chiral discrimination relevant to the enantioseparation of N-acylarylalkylamines by an (R)-phenylglycinol-derived chiral selector

Recently, it was reported that the chiral recognition ability of (R)-N-3,5-dinitrobenzoyl phenylglycinol derivative was examined as a new HPLC chiral stationary phase (CSP 1) for the resolution of racemic N-acylnaphthylalkylamines. However, the mechanism of chiral discrimination on the CSP remained elusive until now. In this study, a spectroscopic investigation of the chiral discrimination mechanism of CSP 1 was undertaken using mixtures of (R)-N-3,5-dinitrobenzoyl phenylglycinol-derived chiral selector (2) and each of the enantiomers of N-acylnaphthylalkylamines (3) by NMR study. First, the differences in free energy changes (DeltaDeltaG) upon diastereomeric complexation in solution between the complex of each isomer with chiral selector 2 by NMR titration were calculated. The values were then compared with those estimated by chiral HPLC. The chemical shift changes of each proton on the chiral selector and analytes were also checked and it was found that the chemical shift changes decreased continuously as the acyl group on analytes increased in length. This observation was consistent with the HPLC data. From these experimental results, the interaction mechanism of chiral discrimination between the chiral selector and the analytes is more precisely explained.

Atomistic modeling of enantioselection in chromatography

A review of atomistic molecular modeling studies related to chromatographic separations of enantiomers is presented. Only those types of calculations where direct interactions between a selector and a selectand are involved are described in this review; omitted are regression models. An emphasis is placed on comparing methods used for sampling potential energy surfaces implementing different methodologies like quantum and molecular mechanics for energy calculations, and molecular dynamics and Monte Carlo sampling strategies for simulations. Type I-V chiral stationary phases and additives for capillary electrophoresis and ion-pair chromatography are covered in this review.

Separation of enantiomers by gas chromatography

The separation of enantiomers by gas chromatography is performed on chiral stationary phases (CSPs) via hydrogen bonding, coordination and inclusion. Thus, typical chiral selectors are amino acid derivatives, terpene-derived metal coordination compounds and modified cyclodextrins. In Chirasil-type stationary phases the chiral selector is anchored to a polysiloxane backbone improving gas chromatographic performance. The present review article describes the state-of-the-art, scope and limitations, applications and mechanistic considerations at the advent of the millennium incorporating 16 figures and 168 references.

Molecular-dynamics simulation for liquid chromatographic interactions: effect of mobile phase composition

A molecular-dynamics simulation method has been applied to investigate the influence of the mobile-phase composition on the retention of solutes in HPLC. The distribution profiles of the distance between two atoms in ODS ligands were constructed to characterize the conformation of ODS ligand molecules. The distinct difference of ODS conformation is observed by comparing molecular models consisting of solvent molecules at each solvent composition. The distribution profiles of the distance between the mobile-phase solvent molecules and ODS ligand molecules were also constructed to characterize the distribution of the solvent molecules at each composition. In all distribution profiles, the difference in the distribution due to a change in the solvent compositions was very clearly found, and the facts seem to be very reasonable. The distribution profiles of the distance between the solute, n-propylbenzene, and the terminal carbon atom in the ODS ligand, and between the solute and the silicon atom in the ODS ligand have been also constructed to see the distribution of the solutes in the separation system. The calculated solute distribution in the ODS-methanol/water system is very consistent with the actual chromatographic retention behaviors.

Enantiodiscrimination by a quinine-based chiral stationary phase: a computational study

A detailed computational study of a derivatized quinine chiral stationary phase (CSP) interacting with enantiomeric 3, 5-dinitrobenzoyl derivatives of leucine was carried out to understand where and how chiral discrimination takes place. The most stable structure of the CSP derived from a conformer search gave a structure whose geometry agrees with an X-ray structure (rmsd 0.6 A). The computed retention order and enantiodiscriminating free energy differences also agree with chromatographic data. The location and characteristics of the analyte binding site were assessed. An evaluation of total energies and intermolecular energies responsible for complex formation and for chiral discrimination was performed. Molecular dynamics trajectories of those intermolecular forces as well as distributions of the stabilizing and destabilizing forces are presented. A partitioning of the CSP into molecular fragments and the role each fragment plays in complexation and chiral recognition is also described.