Assignment of absolute configuration of sulfinyl dilactones: Optical rotations and 1H NMR experiment and DFT calculations

https://doi.org/10.1016/j.molstruc.2010.11.076Get rights and content

Abstract

In previous work, two sulfinyl dilactone compounds were synthesized as key intermediates in the development of new proteasome inhibitors. Two stereogenic centers were introduced during synthesis, including one chiral sulfoxide in a five-membered ring. The absolute configurations (AC) of these two compounds were hard to determine just based on experimental techniques. Therefore, theoretical calculations of chiro-optical and spectroscopic features were used to determine stereochemistry. Ab initio calculations at DFT/B3LYP/6-31G (d, p) level of theory obtained satisfactory agreement with the experimental results. Furthermore, there was good agreement between the calculations of optical rotations and hydrogen chemical shifts in 1H NMR spectra for the final AC assignment as (1′S, 2R, 5R)- and (1′S, 2S, 5S)-5-[1′-(t-Boc-amino)-3′-methylbutyl]-1,3,2-dioxathiolane-2-oxide. This is a successful application of ab initio calculations for the AC assignment for small flexible molecules.

Introduction

The proteasome functionality is indispensable in cellular regulation and protein degradation of eukaryotic cells, and intensive investigations have been focused on the development of new compounds which could interact with the catalytic threonine to inhibit proteasomal proteolysis [1], [2]. Because of its sensitivity to nucleophilic attack, the sulfinyl dilactone functional group was rationally designed as the C-terminus of a potential proteasome inhibitor which could form an irreversible covalent adduct with the catalytic threonine to terminate proteasome function. To obtain a proteasome inhibitor containing a sulfinyl dilactone moiety, we first synthesized a protected aminoalkyl sulfinyl dilactone as the key intermediate [3]. Scheme 1 outlines the general synthetic route. During synthesis, two stereogenic centers were introduced in addition to the S carbon (Cα) already present in the starting material. A stereocarbon was introduced by dihydroxylation of olefin 5 to form diol 6, and a stereo sulfur was introduced by insertion of a sulfoxide moiety.

In a sulfoxide, the oxygen contributes unpaired p electrons to a d orbital of the sulfur, called dπpπ bonding (or dπ) bonding [4]. Two resonance structures can be used to describe the sulfoxide bond: R1(R2)Sdouble bondO  R1(R2)S+single bondO. For sulfinyl dilactones with a sulfoxide in a five-membered ring, the pyramidal inversion barrier has been calculated to be >130 kcal mol−1 [5], [6], [7], [8], and therefore stable, separable stereoisomers should be obtained under ambient conditions. N-protected aminoalkyl sulfinyl dilactones containing two undetermined asymmetric centers could form four possible diastereomers. To eliminate confusion, the S or R stereochemical descriptors can be labeled using subscript C for carbon or S for sulfur, so that the four diastereomers can be represented as SCRCRS, SCRCSS, SCSCRS or SCSCSS (see Fig. 1), with the middle descriptor referring to the ring carbon, and the initial descriptor referring to the acyclic carbon, which is always S based on the starting materials and synthetic route.

Experimentally, we obtained two stereometrically purified sulfinyl dilactones, compounds 1 and 2, from HPLC separation [3]. Our goal in this paper was to determine the absolute configuration (AC) of these two compounds. There were two difficulties: first, the needle-like crystals of the sulfinyl dilactones were too slim to be used for X-ray diffraction; second, there were no literature precedents to be used as reference for assigning the AC for a pyramidal sulfoxide in a five-membered ring. Therefore, ab initio calculations of chiro-optical and spectroscopic features were used to assist in the AC assignment. Density functional theory (DFT) has been extensively exploited to calculate chiro-optical properties such as optical rotation (OR) [9], [10], [11], [12], [13] and circular dichroism (CD) [14], and increasingly, consistency has been obtained between theoretical and experimental results [10]. So we used DFT to calculate OR at the sodium D line wavelength and then compared with experimental results. A medium-sized basis set was used to balance accuracy and computational cost. Before the ab initio calculations, the conformational space was investigated for these flexible diastereomers, in order to be able to include as many as possible of the stable conformations which contribute to the overall chiro-optical and spectroscopic features.

For spectroscopic features, the nuclear magnetic shielding constant at or close to the center of an anisotropic group can be employed to determine quantitatively the stereochemistry of potentially proximal nuclei [15]. A number of such applications have been published [16], [17], [18]. In the sulfinyl dilactone molecule, the introduction of the sulfoxide and amide groups into the molecular framework can cause anisotropy on nearby protons [19], and such anisotropic effects can be seen through comparison of the 1H NMR spectra of compounds 1 and 2. Therefore, we further calculated the proton chemical shifts of the four diastereomers and compared them with the experimental 1H NMR spectra to help specify the absolute configurations of 1 and 2.

Section snippets

Experimental methods

In order to obtain a direct resolution of diastereomers, the sulfinyl dilactones were first re-crystallized in methanol, and then the white needle-like crystalline material was analyzed by high performance liquid chromatography. The analysis was performed using an Agilent 1100 HPLC system (Agilent, Waldbronn, Germany). The mixture of the diastereomers was resuspended in a mobile phase consisting of acetonitrile (A) and H2O (B) with a gradient range from 30% to 60% A for t = 0–20 min and from 60%

Computational methods

Upon inspection of the flexible molecular structure, we expected many low-energy conformations to be relevant for the N-protected aminoalkyl sulfinyl dilactones. The three chiral centers are in close proximity to each other, so there would be no benefit to studying molecular fragments. Hence, stable conformations were generated for the whole molecules. Initial conformations were generated using the random search module implemented in Sybyl8.0 [20], which performs a conformational search

Results and discussion

After the conformational search and geometry optimizations, topological analysis showed that among the four diastereomers, only SCRCSS possessed several folded conformational minima containing internal hydrogen bonding between the polar hydrogen attached to the secondary amine and the sulfoxide oxygen. The other three diastereomers possessed only extended conformations in which the sulfoxide pointed away from the amide group. Fig. 2 depicts the lowest energy geometries for each of the four

Conclusion

In the present work, the ab initio calculations of optical rotations (OR) and 1H chemical shifts using DFT/B3LYP methodology were employed to assign ACs of previously reported novel sulfinyl dilactones. The two different calculated properties produced consistent results, and hence we were able confidently to assign the AC of compound 1 as SCRCRS and compound 2 as SCSCSS. From the OR calculations, we can see that the energetic population analysis was very important for flexible molecules which

Acknowledgments

Partial funding from the University of Mississippi and computer time from the Mississippi Center for Supercomputer Research are greatly appreciated. This investigation was conducted in a facility constructed with support from Research Facilities Improvements Program (C06 RR-14503-01) from the NIH National Center for Research Resources. The synthetic work was supported by the National Natural Science Foundation of China (Grants 30772650 and 20772008).

References (36)

  • A.F. Kisselev et al.

    Chem. Biol.

    (2001)
  • F. Yuste et al.

    Tetrahedron

    (2002)
  • C. Bertucci et al.

    J. Pharm. Biomed. Anal.

    (2010)
  • J.G. Delcros et al.

    Curr. Med. Chem.

    (2003)
  • G. Fu et al.

    J. Chin. Pharm. Sci.

    (2007)
  • R. Bentley

    Chem. Soc. Rev.

    (2005)
  • S. Allenmark et al.

    Chirality

    (1997)
  • H. Marom et al.

    Chirality

    (2007)
  • E.N. Suslova et al.

    Russ. J. General Chem.

    (2007)
  • T.D. Crawford

    Theor. Chem. Acc.

    (2006)
  • T.D. Crawford et al.

    J. Phys. Chem. A

    (2007)
  • G. Mazzeo et al.

    J. Org. Chem.

    (2010)
  • R. Mohammed et al.

    Aust. J. Chem.

    (2010)
  • P.L. Polavarapu

    Chem. Rec.

    (2007)
  • E. Kleinpeter et al.

    J. Phys. Chem. A

    (2008)
  • G. Toth et al.

    Magn. Reson. Chem.

    (2001)
  • E. Kleinpeter et al.

    J. Org. Chem.

    (2004)
  • E. Kleinpeter et al.

    J. Org. Chem.

    (2005)
  • Cited by (7)

    • Isolation and identification of renoprotective substances from the mushroom Ganoderma lucidum

      2015, Tetrahedron
      Citation Excerpt :

      For this purpose, DFT and TD-DFT calculations were carried out at 298 K in the gas phase with Gaussian 09.11 Conformational searches were performed at the molecular mechanics level of theory employing the MMFF force fields.12–18 The conformers with relative energies within 10 kcal/mol of the lowest energy conformer were selected and further geometry optimized at the B3LYP/6-311++G(2d,p) level.

    • Highly unsaturated pyranone derivatives from the basidiomycete Junghuhnianitida

      2014, Tetrahedron Letters
      Citation Excerpt :

      All calculations were carried out at 298 K in the gas phase with Gaussian 09.14 Conformational searches were carried out via employing MMFF94 force fields.20–27 The conformers with relative energy below 6 kcal/mol were selected for optimization at the B3LYP/6-311++G(2d,p) level.23

    • Prenyleudesmanes, rare natural diterpenoids from Dysoxylum densiflorum

      2014, Tetrahedron
      Citation Excerpt :

      All the lowest-energy conformers, which correspond to 99% of the total Boltzmann distribution, were selected for ECD spectra calculation. The Boltzmann factor for each conformer was calculated based on Gibbs free energy.37 Vibrational analysis at the B3LYP/6-311++G(2d,p) level of theory resulted in no imaginary frequencies, confirming the considered conformers as real minima.

    • Organic stereochemistry: Experimental and computational methods

      2015, Organic Stereochemistry: Experimental and Computational Methods
    View all citing articles on Scopus
    View full text