Determination of the absolute configurations of natural products via density functional theory calculations of vibrational circular dichroism, electronic circular dichroism and optical rotation: the schizozygane alkaloid schizozygine

The development of density functional theory (DFT) methods for the calculation of vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and transparent spectral region optical rotation (OR) has revolutionized the determination of the absolute configurations (ACs) of chiral molecules using these chiroptical properties. We report the first concerted application of DFT calculations of VCD, ECD, and OR to the determination of the AC of a natural product whose AC was previously undetermined. The natural product is the alkaloid schizozygine, isolated from Schizozygia caffaeoides. Comparison of DFT calculations of the VCD, ECD, and OR of schizozygine to experimental data leads, for each chiroptical technique, to the AC 2R,7S,20S,21S for the naturally occurring (+)-schizozygine. Three other alkaloids, schizogaline, schizogamine, and 6,7-dehydro-19beta-hydroxyschizozygine, have also been isolated from S. caffaeoides and shown to have structures closely related to schizozygine. Assuming a common biosynthetic pathway, their ACs are defined by that of schizozygine.     

Determination of the absolute configurations of pharmacological natural products via density functional theory calculations of vibrational circular dichroism: the new cytotoxic iridoid prismatomerin

A new highly cytotoxic iridoid has very recently been isolated from Prismatomeris tetrandra and shown to have the structure 3, similar to that of the iridoid oruwacin, 2. We report the determination of the absolute configuration (AC) of the new iridoid, prismatomerin, using vibrational circular dichroism (VCD) spectroscopy. The VCD spectrum of the acetate derivative of 3, 4, is analyzed using the Stephens theory of VCD and density functional theory (DFT). The AC of the naturally occurring 3 is shown to be 1R,5S,8S,9S,10S, identical to that of the naturally occurring iridoid plumericin, 1, also determined using VCD spectroscopy. The [alpha]D values of the natural products 3 and 1 are negative and positive, respectively. Since the ACs of 3 and 1 are identical, it follows that the AC of 3 cannot be correctly determined by empirical comparison of the signs of the [alpha]D values of 3 and 1.     

Determination of absolute configuration using density functional theory calculations of optical rotation and electronic circular dichroism: chiral alkenes

In principle, the absolute configuration (AC) of a chiral molecule can be deduced from its optical rotation (OR) and/or its electronic circular dichroism (ECD). In practice, this requires reliable methodologies for predicting OR and ECD. The recent application of ab initio time-dependent density functional theory (TDDFT) to the calculation of transparent spectral region OR and ECD has greatly enhanced the reliability with which these phenomena can be predicted. TDDFT calculations of OR and ECD are being increasingly utilized in determining ACs. Nevertheless, such calculations are not perfect, and as a result, ACs determined are not 100% reliable. In this paper, we examine the reliability of the TDDFT methods in the case of chiral alkenes. Sodium d line specific rotations, [alpha]D, are predicted for 26 conformationally rigid alkenes of known AC, ranging in size from 5 to 20 C atoms, and with [alpha]D values in the range of 0-500. The mean absolute deviation of predicted [alpha]D values from experimental values is 28.7. With one exception, beta-pinene, the signs of [alpha]D are correctly predicted. Errors in calculated [alpha]D values are approximately random. Our results define a "zone of indeterminacy" within which calculated [alpha]D values cannot be used to determine ACs with >95% confidence. TDDFT ECD spectra are predicted for eight of the alkenes and compared to experimental spectra. Agreement ranges from modestly good to poor, leading to the conclusion that TDDFT calculations of ECD spectra are not yet of sufficient accuracy to routinely provide highly reliable ACs. TDDFT OR calculations for two conformationally flexible alkenes, 3-tert-butylcyclohexene and trans-4-carene, are also reported. For the former, predicted rotations are incorrect in sign over the range 589-365 nm. It is possible that the AC of this molecule has been incorrectly assigned.     

Determination of the absolute configuration of chiral molecules via density functional theory calculations of vibrational circular dichroism and optical rotation: The chiral alkane D3-anti-trans-anti-trans-anti-trans-perhydrotriphenylene

The Absolute configuration (AC) of the chiral alkane D₃-anti-trans-anti-trans-anti-trans-perhydrotriphenylene (PHTP), 1, is determined by comparison of density functional theory (DFT) calculations of its vibrational circular dichroism (VCD) and optical rotation (OR) to the experimental VCD and OR of (+)−1, obtained in high enantiomeric excess using chiral gas chromatography. Conformational analysis of 1 demonstrates that the all-chair (CCCC) conformation is the lowest in energy and that other conformations are too high in energy to be significantly populated at room temperature. The B3PW91/TZ2P calculated IR spectrum of the CCCC conformation of 1 is in excellent agreement with the experimental IR spectrum, confirming the conformational analysis and demonstrating the excellent accuracy of the B3PW91 functional and the TZ2P basis set. The B3PW91/TZ2P calculated VCD spectrum of the CCCC conformation of S-1 is in excellent agreement with the experimental VCD spectrum of (+)−1, unambiguously defining the AC of 1 to be S(+)/R(−). The B3LYP/aug-cc-pVDZ calculated OR of S-1 over the range 589–365 nm has the same sign and dispersion as the experimental OR of (+)−1, further supporting the AC S(+)/R(−). Our results confirm the AC proposed earlier by Farina and Audisio. This study provides a further demonstration of the excellent accuracy of VCD spectra predicted using Stephens’ equation for vibrational rotational strengths together with the ab initio DFT methodology, and further documents the utility of VCD spectroscopy in determining the ACs of chiral molecules.     

Determination of the absolute configurations using electronic and vibrational circular dichroism measurements and quantum chemical calculations

The stereochemistry of products obtained via a chemical reaction may not always be obvious from the reaction scheme utilized. The identification of convenient methods to determine the stereochemistry in such cases is highly desirable. To identify these methods, we considered a substituted 4-vinyl-1-azabicyclo[3.2.0]hept-3-en-7-one that undergoes spontaneous oxidation in the atmosphere at room temperature, yielding an epoxide with unknown absolute configuration. To determine the absolute configuration of the resulting epoxide, three different approaches have been utilized: (a) experimental NOE measurements; (b) experimental electronic circular dichroism (ECD) spectroscopic measurements and their analysis using corresponding quantum chemical predictions at the B3LYP/aug-cc-pVDZ level; (c) experimental vibrational circular dichroism (VCD) spectroscopic measurements and their analysis using corresponding quantum chemical predictions at the B3LYP/aug-cc-pVDZ level. It was found that the NOE data could not provide enough proof for assigning the absolute configuration, while ECD data could not provide enough discrimination to distinguish between the two possible stereoisomers. On the other hand, VCD spectroscopic analysis provided enough discrimination to distinguish between the two possible stereoisomers, and the absolute configuration could be assigned with confidence.     

Determination of the absolute configuration of a chiral oxadiazol-3-one calcium channel blocker, resolved using chiral chromatography, via concerted density functional theory calculations of its vibrational circular dichroism, electronic circular dichroism, and optical rotation

The chiral oxadiazol-3-one 2 has recently been shown to exhibit myocardial calcium entry channel blocking activity, substantially higher than that of diltiazem. To determine the enantioselectivity of this activity, the enantiomers of 2 have been resolved using chiral chromatography. The absolute configuration (AC) of 2 has been determined by comparison of density functional theory (DFT) calculations of its vibrational circular dichroism (VCD) spectrum, electronic circular dichroism (ECD) spectrum, and optical rotation (OR) to experimental VCD, ECD, and OR data. All three chiroptical properties yield identical ACs; the AC of 2 is unambiguously determined to be S(+)/R(-).     

Determination of the absolute configuration of [3(2)](1,4)barrelenophanedicarbonitrile using concerted time-dependent density functional theory calculations of optical rotation and electronic circular dichroism

The technique of time-dependent density functional theory (TDDFT) has very recently been applied to the calculation of both transparent spectral region optical rotations and electronic circular dichroism (CD). Here, we report the concerted application of the new methodologies to the determination of the absolute configuration (AC) of [3(2)](1,4)barrelenophanedicarbonitrile, 1, the first optically active barrelenophane. 1 is conformationally flexible: the two three-carbon bridges of 1 can each exhibit two conformations, leading to three inequivalent conformations of 1: a, b, and c. Conformational structures and energies are predicted using DFT at the B3LYP/6-31G level. Comparison of the calculated structures to structures obtained via X-ray crystallography of (+)-1 shows that (remarkably) all three conformations a-c are simultaneously present in crystalline (+)-1. The sodium D line specific rotations, [alpha](D), and CD spectra of a-c are calculated using TDDFT at the B3LYP/aug-cc-pVDZ level. Comparison of the conformationally averaged specific rotation and CD spectrum to the experimental data of Matsuda-Sentou and Shinmyozu leads to the AC 9S,12S(+)/9R,12R(-). The same AC is obtained both from [alpha](D) and from the CD, strongly supporting its reliability.     

Density functional theory calculations and vibrational circular dichroism of aromatic foldamers

Ab initio calculations together with vibrational circular dichroism (VCD) have been used for studying the conformations of a quinoline-derived oligoamide bearing a terminal chiral residue. Three helically folded conformers of the dimer, trimer, and tetramer forms of the oligomer were optimized at the density functional theory (DFT) level using the B3LYP functional and the 6-31G* basis set. For each form, the three conformers differ in their helical handedness and in the conformation of the chiral end group. The calculated structures of the tetramer and also the proportions predicted between them based on their calculated Gibbs free energies differences match remarkably well with experimental data collected on an octamer. Specifically, a R-phenethyl terminal group gives rise to a 91:9 ratio between left handed and right handed helices. The predicted VCD spectrum calculated from the Boltzmann population of the individual conformer reproduces very well the experimental VCD spectrum of the tetramer in CDCl3 solution. The DFT calculations performed for the trimer also allow one to assess the preferred handedness of the helix and the conformation of the chiral end group, but the calculated relative populations differ slightly from experimental data. Finally, this study shows that the dimer fragment is not sufficient to obtain valuable information on the conformation of this aromatic oligoamide foldamer.     

Determination of absolute configuration using optical rotation calculated using density functional theory

[structure: see text] We report the first determinations of the absolute configurations (ACs) of chiral molecules using discrete frequency, transparent spectral region optical rotations calculated using density functional theory (DFT). The ACs of 2H-naphtho[1,8-bc]thiophene 1-oxide (3), naphtho[1,8-cd]-1,2-dithiole 1-oxide (4), and 9-phenanthryl methyl sulfoxide (5) are determined by comparison of their specific rotations to values calculated via the time-dependent DFT/gauge-invariant atomic orbital (TDDFT/GIAO) methodology using the B3LYP functional and the aug-cc-pVDZ basis set.     

Time dependent density functional theory calculations for electronic circular dichroism spectra and optical rotations of conformationally flexible chiral donor-acceptor dyad

Twelve conformations of a chiral donor-acceptor (charge-transfer) dyad and six conformations of its dimer complex were structurally optimized by using the Kohn-Sham density functional theory (BLYP/TZV2P) incorporating a recently developed empirical correction scheme that uses C6/R6 potentials for van der Waals interactions (DFT-D). Subsequent time-dependent DFT calculations with BH-LYP and B3-LYP functionals (with triple-zeta basis set) were performed to obtain theoretical circular dichroism (CD) spectra. The experimental CD spectra obtained independently were properly reproduced by averaging the calculated spectra of individual conformers according to a Boltzmann population derived from single-point SCS-MP2 energies. The optical rotations of the monomer were also calculated by using the same functionals with an aug-cc-pVDZ basis set. Dielectric continuum solvation models (COSMO) applied to correct the relative energies from the isolated molecule calculations resulted in conformer distributions that piled the same or even poorer level of agreement with the experimental CD spectrum. Our results clearly show the advantage of the DFT-D method for the geometry optimization of large systems with donor-acceptor interactions and the TD-DFT/BH-LYP calculations for reproducing the experimental CD spectra. As compared with the calculated optical rotations, the wealthy information embedded in the experimental/calculated CD spectra is requisite for the configurational and/or conformational analyses of relatively large and flexible chiral organic molecules in solution.     

Determination of absolute configuration using concerted ab Initio DFT calculations of electronic circular dichroism and optical rotation: bicyclo[3.3.1]nonane diones

The concerted use of ab initio time-dependent density functional theory (TDDFT) calculations of transparent spectral region optical rotation and of circular dichroism has recently become practicable, permitting the concerted use of transparent spectral region optical rotation and circular dichroism in determining the absolute configurations of chiral molecules. Here, we report concerted TDDFT calculations of the transparent spectral region specific rotations and of the circular dichroism spectra originating in n --> pi C=O group excitations of four bicyclo[3.3.1]nonane diones, 1-4. Comparison to experiment yields absolute configurations for 1-4. For each dione, specific rotations and circular dichroism spectra give identical absolute configurations. Our results are consistent with previous work, with the exception of the Octant Rule-derived absolute configuration of the 2,9-dione.     

Investigating cyclic sotolon,maple furanone and their dimers in solution using optical rotation,electronic circular dichroism and vibrational circular dichroism

The experimental optical rotation (OR), electronic circular dichroism (ECD) and vibrational circular dichroism (VCD) spectra of (R)-3-hydroxy-4,5-dimethylfuran-2(5H)-one (sotolon, 1) and (R)-5-ethyl-3-hydroxy-4-methylfuran-2(5H)-one (maple furanone, 2) taken in chloroform were compared to their spectra calculated with time-dependent density functional theory (TDDFT). Sotolon was shown to exist as a dimer in chloroform while maple furanone remains a monomer. Transition state barriers for the enol/keto tautomerization of sotolon were calculated and found to be high. The VCD method offers promise to ultimately distinguish between the presence of monomers or dimers.     

Determination of absolute configuration using density functional theory calculation of optical rotation: chiral alkanes

The recently developed Gauge-Invariant (Including) Atomic Orbital (GIAO) based Time-Dependent Density Functional Theory (TDDFT) methodology for the calculation of transparent spectral region optical rotations of chiral molecules provides a new approach to the determination of absolute configurations. Here, we discuss the application of the TDDFT/GIAO methodology to chiral alkanes. We report B3LYP/aug-cc-pVDZ calculations of the specific rotations of the 22 chiral alkanes, 2-23, of well-established Absolute Configuration. The average absolute deviation of calculated and experimental [alpha](D) values for molecules 2-22 is 24.8. In two of the molecules 2-23, trans-pinane, 10, and endo-isocamphane, 13, the sign of [alpha](D) is incorrectly predicted. Our results demonstrate that absolute configurations of alkanes can be reliably assigned by using B3LYP/aug-cc-pVDZ TDDFT/GIAO calculations if, but only if, [alpha](D) is significantly greater than 25. In the case of (-)-anti-trans-anti-trans-anti-trans-perhydrotriphenylene, 1, [alpha](D) is -93 and TDDFT/GIAO calculations reliably lead to the absolute configuration R(-).     

Enantiomeric discorhabdin alkaloids and establishment of their absolute configurations using theoretical calculations of electronic circular dichroism spectra

Enantiomeric pairs of the cytotoxic pyrroloiminoquinone marine alkaloids discorhabdins B (2), G*/I (3), L (4), and W (5) have been isolated from Latrunculia species sponges collected at different locations around the coast of New Zealand. The absolute configuration of all compounds was secured by comparison of observed data with the results of time-dependent density functional theory (TDDFT) calculations of electronic circular dichroism (ECD) spectra. Enantiomeric discorhabdins exhibit equipotent antiproliferative biological activity.     

The current state of ab initio calculations of optical rotation and electronic circular dichroism spectra

The current ability of ab initio models to compute chiroptical properties such as optical rotatory dispersion and electronic circular dichroism spectra is reviewed. Comparison between coupled cluster linear response theory and experimental data (both gas and liquid phase) yields encouraging results for small to medium-sized chiral molecules including rigid species such as (S)-2-chloropropionitrile and (P)-[4]triangulane, as well as conformationally flexible molecules such as (R)-epichlorohydrin. More problematic comparisons are offered by (S)-methyloxirane, (S)-methylthiirane, and (1S,4S)-norbornenone, for which the comparison between theory and experiment is much poorer. The impact of basis-set incompleteness, electron correlation, zero-point vibration, and temperature are discussed. In addition, future prospects and obstacles for the development of efficient and reliable quantum chemical models of optical activity are discussed, including the problem of gauge invariance, scaling of the coupled cluster approach with system size, and solvation.     

Determination of absolute configurations from vibrational circular dichroism: (+)-2-methylthiirane-3,3-d(2)

Experimental vibrational circular dichroism (VCD) spectra for the dextrorotatory enantiomer and theoretical VCD spectra obtained with localized molecular orbital theory using 6-31G* basis set for the (R) configuration of 2-methylthiirane-3,3-d(2) in the 700-1500 cm(-1) region are presented. The observed and predicted VCD signs are in very good agreement suggesting that the dextrorotatory enantiomer has the (R) configuration. This conclusion is also supported by the optical rotational data.     

Determination of the absolute configuration of three as-hydrindacene compounds by vibrational circular dichroism

[Structure: see text]. The absolute configurations of three compounds with a rigid 1,8-disubstituted as-hydrindacene skeleton have been determined using vibrational circular dichroism spectroscopy and quantum chemical calculations. Experimental spectra were compared to B3LYP/6-31G and B3LYP/cc-pVTZ level predicted spectra. Based on the agreement between the predicted and experimental spectra, the stereochemistry could be assigned with high confidence. The results were found to be in agreement with ECD determinations and/or predictions based on the applied asymmetric methods in the synthetic route.     

Solvation of propylene oxide in water: vibrational circular dichroism, optical rotation, and computer simulation studies

The solvation of propylene oxide (PO) in water has been studied using vibrational circular dichroism (VCD) spectroscopy, optical rotation dispersion (ORD) spectroscopy, molecular dynamics simulations, and ab initio calculations. VCD and ORD measurements were carried out for PO as neat liquid, in CCl4, and in water solutions. The classical molecular dynamics simulations were carried out for the PO + water binary mixtures at different concentrations, and the solvation information was derived from the radial distribution functions obtained in the simulations. The total number of water molecules within the closest vicinity of PO was predicted to be about 3. The geometry optimizations, vibrational frequencies, and VCD intensities were evaluated for the PO monomer and the PO-(H2O)n clusters with n = 1-3 , using density functional theory calculations at the B3LYP/aug-cc-pVTZ level of theory. The chirality transfer VCD feature, which is a direct result of the explicit H-bonding between water and the chiral PO solute, was detected experimentally at the water bending band region. This feature exhibits high sensitivity to the solvation structure around PO. Comparison of the calculated and experimental chirality transfer features leads to the conclusion that the PO-water binary complex is the dominating species in aqueous solution at room temperature and the anti conformation, where water is on the opposite side of the oxirane ring of the PO methyl group, is preferred over the syn one. This conclusion is also supported by the complementary ORD studies. Possible contributions from the ternary and quaternary PO-water complexes are also discussed.     

Palmarumycins BG1-BG7 and preussomerin BG1: establishment of their absolute configurations using theoretical calculations of electronic circular dichroism spectra

Palmarumycins BG1-BG7 (1-7), seven new compounds related to palmarumycins, were isolated from the aerial parts of Bruguiera gymnorrhiza as well as a new preussomerin derivative BG1 (8). The structures of these compounds were determined mainly by the analysis of their NMR and MS data, and their relative configurations were assigned on the basis of their (3)J(H,H) coupling constants. Compounds 4 and 7 have a sulfate group that is unprecedented among members of spirodioxynaphthalene-type natural products. The absolute configurations of 1-8 were determined by TDDFT CD calculations of the solution conformers. Compound 5 displayed inhibitory activity against HL 60 and MCF-7 cell lines.     

Vacuum-ultraviolet electronic circular dichroism of L-alanine in aqueous solution investigated by time-dependent density functional theory

The electronic circular dichoism (ECD) of L-alanine in the vacuum-ultraviolet region was calculated for various optimized structures using time-dependent density functional theory (TDDFT) to assign the CD spectrum observed experimentally in aqueous solution down to 140 nm [Matsuo, et al. Chem. Lett. 2002, 826]. The structure of L-alanine in vacuo was optimized using density functional theory (DFT) at the B3LYP/6-31G* level. Its hydrated structure was optimized with nine water molecules (six and three around carboxyl and amino groups, respectively) using DFT and a continuum model (Onsager model). The dihedral angles of carboxyl and amino groups in the optimized hydrated structure differed greatly from those in the crystal and in nonhydrated structures optimized using a continuum model only. The ECD spectrum calculated for the hydrated structure had two successive positive peaks with molar ellipticities of about 2000 deg cm2 dmol(-1) at around 205 and 185 nm, which were close to those observed experimentally. These positive peaks were attributable to n pi* transitions of the carboxyl group, with the latter peak also influenced by the pi pi* transition of the carboxyl group that originates below 175 nm. A small negative peak observed at around 252 nm was also predicted from the hydrated structure. These results demonstrate that the hydrated water molecules around the zwitterions play a crucial role in stabilizing the conformation of L-alanine in aqueous solution and that TDDFT is useful for the ab initio assignment of ECD spectra down to the vacuum-ultraviolet region.