Quantum Chemical Study on the Circular Dichroism Spectra and Specific Rotation of Donor−Acceptor Cyclophanes

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Chiroptical properties of dithia[3.3]cyclophanes composed of anthracene and pyridine/pyridinium moieties: A combined experimental and theoretical study

Circular dichroisms (CDs) of neutral and protonated [3.3]anthracenopyridinophane (1 and 1-H⁺ ) were investigated experimentally and theoretically. Introducing an anthracene moiety with extended conjugation affected the cyclophane structure with the bent angles being appreciably  reduced from those of parent [3.3]pyridinophane. The Cotton effects (CEs) observed at the ¹ B_b band for both 1 and 1-H⁺ were fairly strong and apparently bisignate, which, however, turned out not to be a simple exciton couplet but to be composed of multiple transitions. In contrast, the CEs were much weaker in the ¹ L_a band region. The spectral changes upon protonation were less significant compared with the parent pyridinophane, being dominated by the local transitions of anthracene. Nevertheless, the CD spectra of 1 and 1-H⁺ were well reproduced by theoretical calculations to allow us an unambiguous absolute configuration determination of the first high-performance liquid chromatography (HPLC) elute (from Chiralcel IB column) as S_p . The transannular interactions between the anthracene and pyridine/pyridinium units were examined by UV-vis and fluorescence spectroscopy to reveal a charge-transfer (CT) band in the low-energy region, particularly for 1-H⁺ . Despite the comparable CT interactions, the CE at the CT band was much stronger for the anthracenopyridinophane than for the parent pyridinophane, affording an anisotropy (g) factor as large as 4 × 10^-3 .

Absolute configuration determination through the unique intramolecular excitonic coupling in the circular dichroisms of o,p′-DDT and o,p′-DDD. A combined experimental and theoretical study

Circular dichroisms (CDs) of theo,p′-isomers of 1,1,1-trichloro- and 1,1-dichloro-2,2-bis(chlorophenyl)ethanes (DDTandDDD) were investigated experimentally and theoretically.

Vibronic structure in the far-UV electronic circular dichroism spectra of proteins

The Franck–Condon effect is considered and the vibrational structure of the π_nbπ* transition of the peptide backbone is incorporated into matrix method calculations of the electronic circular dichroism (CD) spectra of proteins in the far-ultraviolet. We employ the state-averaged CASPT2 method to calculate the ground and π_nbπ* excited state geometries and frequencies of N-methylacetamide (NMA), which represents the peptide chromophore. The results of these calculations are used to incorporate vibronic levels of the excited states into the matrix method calculation. The CD spectra of a set of 49 proteins, comprising a range of structural types, are calculated to assess the influence of the vibrational structure. The calculated spectra of α-helical proteins are better resolved using the vibronic parameters and correlation between the experimental and the calculated intensity of less regular β structure proteins improves over most wavelengths in the far-UV. No obvious improvement is observed in the calculated spectra of regular β-sheet proteins. Our high-level ab initio calculations of the vibronic structure of the π_nbπ* transition in NMA have provided some further insight into the physical origins of the nature of protein CD spectra in the far-UV.

Investigation on CD inversion at visible region caused by a tilt of the π-conjugated substituent: theoretical and experimental approaches by using an asymmetric framework of diarylethene annulated isomer

A substituent effect of asymmetric diarylethene annulated isomer on their chiroptical properties was investigated by means of theoretical and experimental approaches. The absolute configuration of the annulated isomer was determined by X-ray structural analysis and DFT calculation. The TD-DFT calculation successfully reproduced not only the sign but also the shape and magnitude of experimental CD spectrum by considering the Boltzmann-weighted average of four atrop-isomers. A fragment decomposition (FD) analysis of rotatory strength clearly revealed a noteworthy effect; the tilting motion concomitant with the rotating motion of the substituent affects the sign and magnitude of CD signals. It was found that even when the absolute structure of the chiral core moiety does not change, the slight motion of the substituent can trigger the inversion of the CD signal.

Absolute configuration of atropisomeric polychlorinated biphenyl 183 enantiomerically enriched in human samples

Polychlorinated biphenyls (PCBs) are still of serious concern as a potential health hazard due to their persistency and bioacumulation. Of 209 possible PCB congeners, with varying number and position of chlorine atom(s), 19 are chiral. These are mostly highly chlorinated and tend to remain longer against the biological decompositions, suffering biological deracemization in the environment. In this work, we have unequivocally determined the absolute configurations of important chiral PCBs 183 and 171, as well as 132, through the combined theoretical and experimental investigations of the chiroptical properties (circular dichroism and optical rotation), which will be valuable in elucidating the mechanism of biological enantiomer enrichment of PCBs in the environment.

Determination of Absolute Configuration of Natural Products: Theoretical Calculation of Electronic Circular Dichroism as a Tool

Determination of absolute configuration (AC) is one of the most challenging features in the structure elucidation of chiral natural products, especially those with complex structures. With revolutionary advancements in the area of quantum chemical calculations of chiroptical spectroscopy over the past decade, the time dependent density functional theory (TDDFT) calculation of electronic circular dichroism (ECD) spectra has emerged as a very promising tool. The principle is simply based on the comparison of the calculated and experimental ECD spectra: the more closely they match, the more reliable conclusion for the AC assignment can be drawn. This review attempts to use several examples representing monomeric flavonoids, rotationally restricted biflavonoids, complex hexahydroxydiphenoyl-containing flavonoids, conformationally flexible and restrained sesquiterpenoids, cembrane-africanene terpenoids, dihydropyranocoumarins, alkaloids, and dihydroxanthones to illustrate the applicability of this approach in determining the AC of structurally diverse natural products. The findings clearly indicate that the TDDFT calculation of ECD spectra can quantify the contribution of individual conformers and the interaction of multiple chromophores, making it possible to determine the AC of complex chiral molecules. The calculated electronic transitions and molecular orbitals provide new insight into the interpretation of ECD spectra at the molecular level.

Absolute configurations of DNA lesions determined by comparisons of experimental ECD and ORD spectra with DFT calculations

The usefulness of modern density functional theory (DFT) methods is considered for establishing the absolute configurations of DNA lesions by comparisons of computed and experimentally measured optical rotatory dispersion (ORD) and electronic circular dichroism (ECD) spectra. Two rigid, structurally different DNA lesions (two spiroiminodihydantoin stereoisomers and four equine estrogen 4-hydoxyequilenin-DNA stereoisomeric adducts) have been investigated. In all cases, the signs and shapes of the computed ORD spectra reproduced the experimentally measured ORD spectra, although the magnitudes of the computed and experimental ORD values do not coincide exactly. The computed ECD spectra also reproduced the shapes of the experimental ECD spectra rather well, but are blue-shifted by 10-20 nm. Since the assignments of the absolute configurations of the DNA lesions studied based on computed and experimental ORD and ECD spectra are fully consistent with one another, the computational DFT method shows significant promise for determining the absolute configurations of DNA lesions. Establishing the stereochemistry of DNA lesions is highly useful for understanding their biological impact, especially when sufficient amounts of material are not available for other methods of structural characterization.

A Combined Experimental and Theoretical Study on the Conformation of Multiarmed Chiral Aryl Ethers

Four series of multiarmed chiral aryl ethers carrying two, three, five, or eight side-chains on a variety of aromatic core molecules (2-5) were prepared. The structure and conformation of 2 and 3 (in the solid state) were determined by the X-ray crystallographic analyses. While a pair of alternated (anti) conformers (i.e, up-down and down-up) were found in the crystal of 2, three side-arms in 3 were aligned in the same direction to give a C(3)-symmetric syn-conformation. Examinations by dispersion-corrected density functional (DFT-D) calculations revealed that two out of six anti- and two out of four syn-conformers of 2 are energetically most important. Two calculated structures of anti-conformers are in good agreement with those found in the solid state by X-ray analysis. Similarly, relevant conformations of syn-3, fully alternated 4, and C(5)-symmetric 5 were optimized at the DFT-D-B-LYP/TZVP level. The structure and conformation of the side-arms in 2-5 in solution were further studied by temperature dependent (1)H NMR and UV-vis spectroscopy. In addition, comparative experimental and theoretical CD spectral studies were carried out in order to elucidate the contribution of the thermodynamically less-stable minor isomers in solution. The CD spectral changes observed for 2 and 3 at varying temperatures were quite different, while the parent chiral arene 1, as well as 4 and 5, only showed an increased intensity of the negative Cotton effect for the (1)L(b) band. The latter behavior is readily accounted for in terms of the conformational freezing of the chiral groups at low temperatures. The unusual CD spectral behavior observed for 2 and 3 was rationalized by the conformational alteration of the side-arms. Because of attractive van der Waals interactions between the aromatic units of the arms in nonpolar solvents, the syn-conformations become gradually more important for 2 at low temperatures, which eventually results in a weak positive Cotton effect for the (1)L(b) band. This was also supported by the SCS-MP2/TZVPP single-point energy calculations for the relevant conformers of 2. For 3, the contribution of the C(3)-symmetrical conformer becomes more important than the less-symmetrical isomers at low temperatures. The conformations of 2 and 3 in their excited states as well as in the oxidized states were also examined.