Determination of the absolute configuration of flexible molecules by ab initio ORD calculations: a case study with cytoxazones and isocytoxazones

Ab initio calculations of the optical rotatory power of the natural cytokine modulator cytoxazone 1 and its trans-diastereomer 2, as well as the structural isomers cis-3 and trans-4 isocytoxazones, have been performed at four different wavelengths (589, 546, 435, and 405 nm) by Density Functional Theory. The calculation of ORD curves provides a reliable method for the assignment of absolute configuration of these conformationally flexible molecules. The absolute configurations of isocytoxazones has been established as (+)-(4R,5S)-cis-3 and (+)-(4S,5S)-trans-4.     

Determining absolute configuration in flexible molecules: a case study

Assigning absolute configuration of molecules continues to be a major problem. Determining absolute configuration in conformationally flexible systems is challenging, even for experts. Here, we present a case study in which we use a combination of molecular modeling, solution NMR, and X-ray crystallography to illustrate why it is difficult to use solution methods alone for configuration assignment. For the case examined, a comparison of calculated and experimental optical rotatory dispersion (ORD) data provides the most straightforward way to assign the absolute configuration.     

Absolute configuration of C76 from optical rotatory dispersion.

The absolute configuration of C76 has been determined as (+)589-(fC)-C76 , for the first time, by comparing the experimental and predicted optical rotatory dispersion (ORD) patterns. The experimental ORD pattern was derived from the experimental electronic circular dichroism (ECD) spectrum using the Kramers-Kronig (KK) transform. The theoretical ORD spectra were calculated in the resonant region using linear response theory, and also using the KK transform of the theoretical ECD spectrum, at different theoretical levels, namely BHLYP/6-31G*, B3LYP/6-31G*, BLYP/6-31G*, and HF/6-31G*. Good agreement noted between experimental and predicted spectra allows for an unambiguous determination of the absolute configuration.     

Comparison of time-dependent density-functional theory and coupled cluster theory for the calculation of the optical rotations of chiral molecules

A comparison of the abilities of time-dependent density-functional theory (TDDFT) and coupled cluster (CC) theory to reproduce experimental sodium D-line specific rotations for 13 conformationally rigid organic molecules is reported. The test set includes alkanes, alkenes, and ketones with known absolute configurations. TDDFT calculations make use of gauge-including atomic orbitals and give origin-independent specific rotations. CC rotations are computed using both the origin-independent dipole-velocity and origin-dependent dipole-length representations. The mean absolute deviations of calculated and experimental rotations are of comparable magnitudes for all three methods. The origin-independent DFT and CC methods give the same sign of [alpha]D for every molecule except norbornanone. For every large-rotation ketone and alkene for which DFT and CC yield the incorrect sign as compared to liquid-phase experimental data, the corresponding optical rotatory dispersion (ORD) curve is bisignate, suggesting that the two models cannot reliably reproduce the relative excitation energies and antagonistic rotational strengths of multiple competing electronic states that contribute to the total long-wavelength rotation. Several potential sources of error in the theoretical treatments are considered, including basis set incompleteness, vibrational and temperature effects, electron correlation, and solvent effects.     

Time dependent density functional theory modeling of specific rotation and optical rotatory dispersion of the aromatic amino acids in solution

Time Dependent Density Functional Theory (TDDFT) along with the COnductor-like Screening MOdel (COSMO) has been applied to model the specific rotation at 589.3 nm and the optical rotatory dispersion (ORD) of the aromatic amino acids phenylalanine, tyrosine, histidine, and tryptophan. Solution structures at low, neutral, and high pH were determined. Both the anomalous dispersion absorbing (resonance) region and the lower energy (transparent) region of the ORD of the compounds were modeled. Linear response calculation of the specific rotation and ORD as well as Kramers-Kronig transformations of calculated circular dichroism spectra to model resonant ORD were compared with experimental data from the literature.     

Bijvoet in solution reveals unexpected stereoselectivity in a Michael addition

The absolute configuration of small crystallizable molecules can be determined with anomalous X-ray diffraction as shown by Bijvoet in 1951. For the majority of compounds that can neither be crystallized nor easily be converted into crystallizable derivatives, stereocontrolled organic synthesis is still required to establish their absolute configuration. In this contribution, a new fundamental methodology for resolving the absolute configuration will be presented that does not require crystallization. With residual dipolar coupling enhanced NMR spectroscopy, ensembles of a limited number of structures are created reflecting the correct conformations and relative configuration. Subsequently, from these ensembles, optical rotation dispersion (ORD) spectra are predicted by DFT calculations and compared to experimental results. The combination of these two steps reveals the absolute configuration of a flexible molecule in solution, which is a big challenge to chiroptical methods and DFT in the absence of NMR spectroscopy. Here the absolute stereochemistry of the product of a new Michael addition, synthesized via a niobium(V) chiral enolate, will be elucidated by using the new methodology.     

Absolute configuration of C2-symmetric spiroselenurane: 3,3,3',3'-tetramethyl-1,1'-spirobi[3 H,2,1]benzoxaselenole

The enantiomers of 3,3,3',3'-tetramethyl-1,1'-spirobi[3 H,2,1]benzoxaselenole have been separated on a chiral preparative chromatographic column. The experimental vibrational circular dichroism (VCD) spectra have been obtained for both enantiomers in CH(2)Cl(2). The theoretical VCD spectra have been obtained by means of density functional theoretical calculations with the B3 LYP density functional. From a comparison of experimental and theoretical VCD spectra, the absolute configuration of an enantiomer with positive specific rotation in CH(2)Cl(2) at 589 nm is determined to be R. This conclusion has been verified by comparing results of experimental optical rotatory dispersion (ORD) and electronic circular dichroism (ECD) to predictions of the same properties using the B3 LYP functional for the title compound.     

The optical activity of carvone: a theoretical and experimental investigation

The optical rotatory dispersion (ORD) and circular dichroism of the conformationally flexible carvone molecule has been investigated in 17 solvents and compared with results from calculations for the "free" (gas phase) molecule. The G3 method was used to determine the relative energies of the six conformers. The optical rotation of (R)-(-)-carvone at 589 nm was calculated using coupled cluster and density functional methods, including temperature-dependent vibrational corrections. Vibrational corrections are significant and are primarily associated with normal modes involving the stereogenic carbon atom and the carbonyl group, whose n → π? excitation plays a significant role in the chiroptical response of carvone. Without the inclusion of vibrational corrections the optical rotation calculated with CCSD and DFT has the opposite sign of experimental data. Calculations of optical rotation performed in solution using the polarizable continuum model were also opposite in sign when compared to that of the experiment.     

Optical rotatory dispersion from liquid crystalline solutions and films of hydroxypropylcellulose

The chiroptical properties of (hydroxypropyl)cellulose (HPC) in methanol are measured by optical rotatory dispersion (ORD), both in isotropic solution and in the cholesteric liquid crystalline phase. The ORD spectra of HPC films cast from lyotropic solution are also examined. Isotropic solutions of HPC in methanol provide no ORD evidence for the presence of a helical conformation, or for concentration-induced changes in conformation. The ORD curve of HPC in methanol, when expressed as specific rotation, is independent of polymer concentration for isotropic solutions containing between 1 and 43 per cent polymer by weight. From the slope of Drude plots, the Cotton effect responsible for the observed ORD curve occurs at 175-180 nm. However a 45 per cent liquid crystal-line solution exhibits plain positive dispersion, and the magnitude of the specific rotation is also much greater than that found for solutions containing less than 43 per cent polymer. Results for more concentrated solutions confirm that ORD spectra of lyotropic liquid crystals of HPC in methanol contain a significant positive component that is not accounted for by the de Vries equation for cholesteric reflection. ORD measurements in other solvents and on dry films show that this contribution to the optical activity vanishes when the solvent is removed or when hydrogen bonding is disrupted. The effect is tentatively ascribed to a hydrogen bonded structure that contributes strongly to the optical rotation.     

Determination of the absolute configurations of natural products via density functional theory calculations of vibrational circular dichroism, electronic circular dichroism, and optical rotation: the iridoids plumericin and isoplumericin

The absolute configurations (ACs) of the iridoid natural products, plumericin (1) and isoplumericin (2), have been re-investigated using vibrational circular dichroism (VCD) spectroscopy, electronic circular dichroism (ECD) spectroscopy, and optical rotatory dispersion (ORD). Comparison of DFT calculations of the VCD spectra of 1 and 2 to the experimental VCD spectra of the natural products, (+)-1 and (+)-2, leads unambiguously to the AC (1R,5S,8S,9S,10S)-(+) for both 1 and 2. In contrast, comparison of time-dependent DFT (TDDFT) calculations of the ECD spectra of 1 and 2 to the experimental spectra of (+)-1 and (+)-2 does not permit definitive assignment of their ACs. On the other hand, TDDFT calculations of the ORD of (1R,5S,8S,9S,10S)-1 and -2 over the range of 365-589 nm are in excellent agreement with the experimental data of (+)-1 and (+)-2, confirming the ACs derived from the VCD spectra. Thus, the ACs initially proposed by Albers-Sch?nberg and Schmid are shown to be correct, and the opposite ACs recently derived from the ECD spectra of 1 and 2 by Els?sser et al. are shown to be incorrect. As a result, the ACs of other iridoid natural products obtained by chemical correlation with 1 and 2 are not in need of revision.     

Non-empirical assignment of the absolute configuration of (-)-naringenin, by coupling the exciton analysis of the circular dichroism spectrum and the ab initio calculation of the optical rotatory power

The non-empirical assignment of the absolute configuration of (-)-naringenin, the aglycone of (-)-naringin, a flavanone glycoside abundant in the albedo of immature grapefruits and showing several interesting biological properties, has been approached by two different methods: (a) the exciton analysis of the circular dichroism (CD) spectrum and (b) the ab initio calculation of the optical rotatory power. Both the methods indicate the configurational correlation (-)/(S), as empirically suggested by Gaffield. A comparison of advantages and limitations of the two methods of analysis is also presented.     

A single chiroptical spectroscopic method may not be able to establish the absolute configurations of diastereomers: dimethylesters of hibiscus and garcinia acids

Electronic circular dichroism (ECD), optical rotatory dispersion (ORD), and vibrational circular dichroism (VCD) spectra of hibiscus acid dimethyl ester have been measured and analyzed in combination with quantum chemical calculations of corresponding spectra. These results, along with those reported previously for garcinia acid dimethyl ester, reveal that none of these three (ECD, ORD, or VCD) spectroscopic methods, in isolation, can unequivocally establish the absolute configurations of diastereomers. This deficiency is eliminated when a combined spectral analysis of either ECD and VCD or ORD and VCD methods is used. It is also found that the ambiguities in the assignment of absolute configurations of diastereomers may also be overcome when unpolarized vibrational absorption is included in the spectral analysis.     

Determining the absolute configuration of (+)-mefloquine HCl, the side-effect-reducing enantiomer of the antimalaria drug Lariam

Even though the important antimalaria drug rac-erythro-mefloquine HCl has been on the market as Lariam for decades, the absolute configurations of its enantiomers have not been determined conclusively. This is needed, since the (-) enantiomer is believed to cause adverse side effects in malaria treatment resulting from binding to the adenosine receptor in the human brain. Since there are conflicting assignments based on enantioselective synthesis and anomalous X-ray diffraction, we determined the absolute configuration using a combination of NMR, optical rotatory dispersion (ORD), and circular dichroism (CD) spectroscopy together with density functional theory calculations. First, structural models of erythro-mefloquine HCl compatible with NMR-derived (3)J(HH) scalar couplings, (15)N chemical shifts, rotational Overhauser effects, and residual dipolar couplings were constructed. Second, we calculated ORD and CD spectra of the structural models and compared the calculated data with the experimental values. The experimental results for (-)-erythro-mefloquine HCl matched our calculated chiroptical data for the 11R,12S model. Accordingly, we conclude that the assignment of 11R,12S to (-)-erythro-mefloquine HCl is correct.     

Microscale determination of the absolute configuration of alpha-aryl-substituted alcohols by the CD exciton chirality method

The absolute configurations of a broad spectrum of aryl alcohols 1 have been determined for the first time by the CD exciton chirality method. The configurational assignment is additionally verified by computer modeling and lipase-catalyzed acetylation of the racemic alcohols. The CD-spectroscopic data have revealed that the S enantiomers of the benzoate derivatives 2 display a positive first Cotton effect and the R enantiomers a negative one at around 228 nm. Thus, the sense of the first Cotton effect of the benzoate derivative 2 allows a reliable assignment of the absolute configuration of the corresponding alcohol 1.     

Synthesis of the Natural Product Building Block 5‐(3‐Bromophenyl)‐4‐hydroxy‐5‐methylhexan‐2‐one and its Chiral Characterization by Using Chiroptical Spectroscopy

The absolute configuration of 5‐(3‐bromophenyl)‐4‐hydroxy‐5‐methylhexan‐2‐one, an intermediate in the synthesis of various natural products, is assigned by using vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and optical rotatory dispersion (ORD). Experimental spectra were compared to density functional theory (DFT) calculations of the molecule with known configuration. These three techniques independently confirm that the absolute configuration is (S)‐5‐(3‐bromophenyl)‐4‐hydroxy‐5‐methylhexan‐2‐one, thus enabling us to assign the absolute configuration with high reliability. The reliability of the VCD analysis was assessed quantitatively by using the CompareVOA program. We found that, in cases in which the agreement between theory and experiment was very good, a value of 10 cm^?1 for the triangular weighting function gave a more‐realistic discriminative power between enantiomers than the default value of 20 cm^?1.     

Ab initio calculation of optical rotation in (P)-(+)-[4]triangulane

Optical rotation, the angle through which plane-polarized light rotates when passed through an enantiomerically pure medium, plays a vital role in the determination of the absolute configurations of chiral molecules such as natural products. We describe new quantum mechanical methodology designed to assist in this endeavor by providing high-accuracy computational optical rotatory dispersion data for matching to experimental results. Comparison between theory and experiment for the rigid, helical molecule trispiro[2.0.0.2.1.1]nonane [also known as (P)-(+)-[4]triangulane], recently synthesized with enantiomeric purity, shows that the coupled cluster quantum chemical model provides superb agreement for optical rotation across a wide range of wavelengths (589-365 nm), with errors averaging only 1%.     

Electron donor-acceptor complexes in polypeptides. 3. Optical activity of intra-and inter-molecular charge-transfer associated with absorption bands

Optical activity, i. e. circular dichroism (CD) and optical rotatory dispersion (ORD) associated with charge-transfer (CT) absorption bands were found for intra- and intermolecular electron-donor-acceptor (EDA)-complexes containing centers of configurational asymmetry.     

Spectroscopic investigation of the structures of dialkyl tartrates and their cyclodextrin complexes

Structures of three dialkyl tartrates, namely, dimethyl tartrate, diethyl tartrate, and diisopropyl tartrate, in CCl4, dimethyl sulfoxide (DMSO)/DMSO-d6, and H2O/D2O solvents have been investigated using vibrational absorption (VA), vibrational circular dichroism (VCD), and optical rotatory dispersion (ORD). VA, VCD, and ORD spectra are found to be dependent on the solvent used. Density functional theory (DFT) calculations are used to interpret the experimental data in CCl4 and DMSO. The trans-COOR conformer with hydrogen bonding between the OH group and the C=O group attached to the same chiral carbon is dominant for dialkyl tartrates both in vacuum and in CCl4. The experimental VA, VCD, and ORD data of dialkyl-D-tartrates in CCl4 correlated well with those predicted for dimethyl-(S,S)-tartrate molecule as both isolated and solvated in CCl4. In DMSO solvent, dialkyl tartrate molecules favor formation of intermolecular hydrogen bonding with DMSO molecules. Clusters of dimethyl-(S,S)-tartrate, with one molecule of dimethyl-(S,S)-tartrate hydrogen bonded to two DMSO molecules, are used for the DFT calculations. A trans-COOR cluster and a trans-H cluster are needed to obtain a reasonable agreement between the predicted and experimental data of dimethyl tartrate in DMSO solvent. VA, VCD, and optical rotations are also measured for dialkyl tartrate-cyclodextrin complexes. It is noted that these properties are barely affected by complexation of dialkyl tartrates with cyclodextrins, indicating weak interaction between tartrates and cyclodextrin. Binding constants of alpha-CD and beta-CD with diethyl L-tartrate in both H2O and DMSO have been determined using isothermal titration calorimetry technique. The smaller binding constants (less than 100) confirmed the weak interaction between tartrates and cyclodextrin in the solution state.     

Calculation of circular dichroism spectra from optical rotatory dispersion, and vice versa, as complementary tools for theoretical studies of optical activity using time-dependent density functional theory

A comparison of two theoretical methods based on time-dependent density functional theory for the calculation of the linear dispersive and absorptive properties of chiral molecules has been made. For this purpose, a recently proposed computational method for the calculation of circular dichroism (CD) spectra from the imaginary part of the optical rotation parameter has been applied to six rigid organic molecules. The results have been compared to the CD spectra obtained from the rotatory strengths and from the Kramers-Kronig transformation of optical rotatory dispersion (ORD) curves. We have also investigated a criterion based on the Kramers-Kronig integration formula to determine a number of excitations in truncated CD spectra which may yield a reasonable low frequency resonant ORD. It has been tested by calculating the ORD from the sum-over-states formula both in the nonresonant and resonant regions. Finally, we have applied these methods to model the resonant optical activity of proline at low pH.     

Chiroptical spectroscopic determination of molecular structures of chiral sulfinamides: t-butanesulfinamide

The absolute configuration of t-butanesulfinamide has been determined as (-)-(S) using three different chiroptical spectroscopic methods, namely, vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and optical rotatory dispersion (ORD). Furthermore, the predominant conformation of this molecule is determined to have S=O and NH2 groups staggered with respect to the three methyl groups and to have amine hydrogens in gauche orientation with respect to S=O. The quality of predictions obtained for vibrational properties, namely, vibrational absorption and VCD, is found to be satisfactory with the B3LYP functional and 6-31G* basis set. However, this basis set is found to be inadequate for obtaining reliable predictions of electronic properties, namely, electronic absorption and ECD, but a larger aug-cc-pVDZ basis set is found to provide satisfactory prediction of electronic properties. t-Butanesulfinamide serves as an example which invalidates the recommendation of using the 6-31G* basis set for molecules that exhibit the same sign for the long-wavelength ECD band and ORD. This molecule also emphasizes the importance of simultaneous investigation of ECD and ORD, and the use of multiple chiroptical spectroscopic methods, for reliable determination of molecular stereochemistry.