Enantioresolution and stereochemical characterization of two chiral sulfoxides endowed with COX-2 inhibitory activity

A journey in unraveling the enantiorecognition mechanism of 3,5-dinitrobenzoyl-amino acids with two Cinchona alkaloid-based chiral stationary phases: The power of molecular dynamic simulations

BACKGROUND

Innovations in computer hardware and software capabilities have paved the way for advances in molecular modelling techniques and methods, leading to an unprecedented expansion of their potential applications. In contrast to the docking technique, which usually identifies the most stable selector-selectand (SO-SA) complex for each enantiomer, the molecular dynamics (MD) technique enables the consideration of a distribution of the SO-SA complexes based on their energy profile. This approach provides a more truthful representation of the processes occurring within the column. However, benchmark procedures and focused guidelines for computational treatment of enantioselectivity at the molecular level are still missing.

RESULTS

Twenty-eight molecular dynamics simulations were performed to study the enantiorecognition mechanisms of seven N-3,5-dinitrobenzoylated α- and β-amino acids (DNB-AAs), occurring with the two quinine- and quinidine-based (QN-AX and QD-AX) chiral stationary phases (CSPs), under polar-ionic conditions. The MD protocol was optimized in terms of box size, simulation run time, and frame recording frequency. Subsequently, all the trajectories were analyzed by calculating both the type and amount of the interactions engaged by the selectands (SAs) with the two chiral selectors (SOs), as well as the conformational and interaction energy profiles of the formed SA-SO associates. All the MDs were in strict agreement with the experimental enantiomeric elution order and allowed to establish (i) that salt-bridge and H-bond interactions play a pivotal role in the enantiorecognition mechanisms, and (ii) that the π-cation and π-π interactions are the discriminant chemical features between the two SOs in ruling the chiral recognition mechanism.

SIGNIFICANCE

The results of this work clearly demonstrate the high contribution given by MD simulations in the comprehension of the enantiorecognition mechanism with Cinchona alkaloid-based CSPs. However, from this research endeavor it clearly emerged that the MD protocol optimization is crucial for the quality of the produced results.

Our research cooperation with Professor Yoshio Okamoto

This article summarizes our cooperation with the research group of Prof. Yoshio Okamoto at Nagoya University during the period of time between 1992 and 2005. Although the text deals entirely with enantioseparations in high-performance liquid chromatography, capillary electrophoresis, and capillary electrochromatography, this is not a detailed review in any of these areas. The text highlights selected aspects of these techniques, which have been the subject of our joint research and in part their reflection in follow-up research by our and other research groups. Together with more systematically studied topics, aspects such as ultrafast separation of enantiomers, uncommonly high separation factor of enantiomers and other related issues are also addressed.

Lipase Assisted (S)-Ketoprofen Resolution from Commercially Available Racemic Mixture

Ketoprofen is a commercially available drug sold as a racemic mixture that belongs to the family of non-steroidal anti-inflammatory drugs known as profens. It has been demonstrated (in vitro) that (S)-ketoprofen is around 160 times more potent than its enantiomer (R)-ketoprofen, while accumulation of (R)-ketoprofen can cause serious side effects, such as dyspepsia, gastrointestinal ulceration/bleeding, pain, salt and fluid retention, and hypertension. In this work, four commercially available lipases were systematically assessed. Parameters such as conversion, enantiomeric excess, and enantioselectivity were considered. Among them, and by evaluating lipase load, temperature, solvent, and alcohol, Candida rugosa lipase exhibited the best results in terms of enantioselectivity E = 185 ((S)-enantiopreference) with esterification conversions of c = 47% (out of 50%) and enantiomeric excess of 99%. The unreacted (R)-enantiomer was recovered by liquid-liquid extraction and racemized under basic media, which was recycled as starting material. Finally, the (S)-alkyl ketoprofen ester was successfully enzymatically hydrolyzed to the desired (S)-ketoprofen with c = 98.5% and 99% ee. This work demonstrated the benefit and efficiency of using Candida rugosa lipase to kinetically resolve racemic ketoprofen by an environmentally friendly protocol and with the recycling of the undesired (R)-ketoprofen.

Synthesis and Configurational Assignment of Vinyl Sulfoximines and Sulfonimidamides

Vinyl sulfones and sulfonamides are valued for their use as electrophilic warheads in covalent protein inhibitors. Conversely, the S(VI) aza-isosteres thereof, vinyl sulfoximines and sulfonimidamides, are far less studied and have yet to be applied to the field of protein bioconjugation. Herein, we report a range of different synthetic methodologies for constructing vinyl sulfoximine and vinyl sulfonimidamide architectures that allows access to new areas of electrophilic chemical space. We demonstrate how late-stage functionalization can be applied to these motifs to incorporate alkyne tags, generating fully functionalized probes for future chemical biology applications. Finally, we establish a workflow for determining the absolute configuration of enantioenriched vinyl sulfoximines and sulfonimidamides by comparing experimentally and computationally determined electronic circular dichroism spectra, enabling access to configurationally assigned enantiomeric pairs by separation.

Enantioseparation, Stereochemical Assignment and Chiral Recognition Mechanism of Sulfoxide-Containing Drugs

The distinct pharmacodynamic and pharmacokinetic properties of enantiopure sulfoxide drugs have stimulated us to systematically investigate their chiral separation, stereochemical assignment, and chiral recognition mechanism. Herein, four clinically widely-used sulfoxide drugs were chosen and optically resolved on various chiral stationary phases (CSPs). Theoretical simulations including electronic circular dichroism (ECD) calculation and molecular docking were adopted to assign the stereochemistry and reveal the underlying chiral recognition mechanism. Our results showed that the sequence of calculated mean binding energies between each pair of enantiomers and CSP matched exactly with experimentally observed enantiomeric elution order (EEO). It was also found that the length of hydrogen bond might contribute dominantly the interaction between two enantiomers and CSP. We hope our study could provide a fresh perspective to explore the stereochemistry and chiral recognition mechanism of chiral drugs.

Determination of absolute configuration and binding efficacy of benzimidazole-based FabI inhibitors through the support of electronic circular dichroism and MM-GBSA techniques

We have previously reported benzimidazole-based compounds to be potent inhibitors of FabI for Francisella tularensis (FtFabI), making them promising antimicrobial hits. Optically active enantiomers exhibit markedly differing affinities toward FtFabI. The IC₅₀ of benzimidazole (-)-1 is ∼100× lower than the (+)-enantiomer, with similar results for the 2 enantiomers. Determining the absolute configuration for these optical compounds and elucidating their binding modes is important for further design. Electronic circular dichroism (ECD) quantum calculations have become important in determining absolute configurations of optical compounds. We determined the absolute configuration of (-)/(+)-1 and (-)/(+)-2 by comparing experimental spectra and theoretical density functional theory (DFT) simulations of ECD spectra at the B3LYP/6-311+G(2d, p) level using Gaussian09. Comparison of experimental and calculated ECD spectra indicates that the S configuration corresponds to the (-)-rotation for both compounds 1 and 2, while the R configuration corresponds to the (+)-rotation. Further, molecular dynamics simulations and MM-GBSA binding energy calculations for these two pairs of enantiomers with FtFabI show much tighter binding MM-GBSA free energies for S-1 and S-2 than for their enantiomers, R-1 and R-2, consistent with the S configuration being the more active one, and with the ECD determination of the S configuration corresponding to (-) and the R configuration corresponding to (+). Thus, our computational studies allow us to assign (-) to (S)- and (+) to (R)- for compounds 1 and 2, and to further evaluate structural changes to improve efficacy.