First Principles Simulations of Optical Rotation of Chiral Molecular Crystals

In this work, we present simulations of the optical rotation (OR) for five molecular crystals at density functional theory level with periodic boundary conditions (DFT-PBC). Calculations are compared with experimental measurements and show semi-quantitative agreement with experimental data for three of the crystals: tartatic acid, benzil, and pentaerythritol. For the other two crystals, aspartic acid and glutamic acid, the calculated data are in qualitative agreement with, but two orders of magnitude smaller than, the experimental data. We provide some arguments that support the theoretical predictions and suggest that the experiments should be revisited. We also find that the position of H centers provided in experimental X-ray data is not sufficiently reliable for simulating OR, and better results are obtained when H atoms are allowed to relax while keeping heavier elements fixed at the experimental positions. Comparison with molecular cluster calculations with a better functional and a larger basis set indicate that the role of intermolecular interactions (reproduced with the PBC technique) is as or more important than the choice of model chemistry. Despite the current limitations in the level of theory that can be employed, these simulations provide a promising avenue to investigate the effect of intermolecular interactions on this sensitive electronic property of molecules and materials.

Conformational analysis of lipid membrane mimetics modified with Aβ42 peptide by Raman spectroscopy and computer simulations

Peptide-lipid interactions play an important role in maintaining the integrity and function of the cell membrane. Even slight changes in these interactions can induce the development of various diseases. Specifically, peptide misfolding and aggregation in the membrane is considered to be one of the triggers of Alzheimer's disease (AD), however its exact mechanism is still unclear. To this end, an increase of amyloid-beta (Aβ) peptide concentration in the human brain is widely accepted to gradually produce cytotoxic Aβ aggregates (plaques). These plaques initiate a sequence of pathogenic events ending up in observable symptoms of dementia. Understanding the mechanism of the Aβ interaction with cells is crucial for early detection and prevention of Alzheimer's disease. Hence, in this work, a comprehensive Raman analysis of the Aβ42 conformational dynamics in water and in liposomes and lipodiscs that mimic the membrane system is presented. The obtained results show that the secondary structure of Aβ42 in liposomes is dominated by the α-helix conformation, which remains stable over time. However, it comes as a surprise to reveal that the lipodisc environment induces the transformation of the Aβ42 secondary structure to a β-turn/random coil. Our Raman spectroscopy findings are supported with molecular dynamics (MD) and density functional theory (DFT) simulations, showing their good agreement.Communicated by Ramaswamy H. Sarma.

Effective fully polarizable QM/MM approaches to compute Raman and Raman Optical Activity spectra in aqueous solution

Raman and Raman Optical Activity (ROA) signals are amply affected by solvent effects, especially in the presence of strongly solute-solvent interactions such as Hydrogen Bonding (HB). In this work, we extend the fully atomistic polarizable Quantum Mechanics/Molecular Mechanics approach, based on the Fluctuating Charges and Fluctuating Dipoles force field to the calculation of Raman and ROA spectra. Such an approach is able to accurately describe specific HB interactions, by also accounting for anisotropic contributions due to the inclusion of fluctuating dipoles. To highlight the potentiality of the novel approach, Raman and ROA spectra of L-Serine and L-Cysteine dissolved in aqueous solution are computed and compared both with alternative theoretical approaches and experimental measurements.

Matrix effects in MI-VCD spectra of two chiral oxiranes and their potential microscopic origin

Combining vibrational circular dichroism (VCD) spectroscopy with the matrix isolation (MI) technique opens up interesting possibilities to study chiral molecules. MI involves the isolation of guest species in inert solid matrices at cryogenic temperatures. Hence, MI-VCD measures are solid-state VCD measurements, and as such, can suffer from mostly birefringance-related artefacts in the same way as common solid-state VCD measurements. In this contribution, we demonstrate that the sample preparation condition have tremendous impact on the quality and reliability of the recorded MI-VCD spectra. While MI-IR spectra are basically blind to these artefacts, the variation of deposition temperatures and host flow rates seem to control whether high quality MI-VCD spectra are obtained or if depolarization effects lead to completely obscured spectra. For two selected examples, styrene oxide (SO) and 1-phenyl propylene oxide (PPO), we discuss how the various experimental conditions may lead to the aforementioned effects and give a microscopic description of their origin.

Synthesis and VCD Spectroscopic Characterization of a Series of Azacryptands from a Chiral Valine-Based Derivative of Tris(2-aminoethyl)amine (TREN)

Utilizing experimental and computational vibrational circular dichroism (VCD) spectroscopy, we explored the conformational preferences of a series of chiral C3 -symmetric octaazacryptands with tris(2-aminoethyl)-amine head groups derived from valine. While the spectra of the smallest azacryptand with p-phenyl linkers and its elongated derivative with p-biphenyls linker were found to match well with the computed spectra, the computed conformational preferences of the m-biphenyl-based azacryptand did not seem to reflect the conformations dominating in chloroform solution. A detailed analysis revealed that structural changes resulting in a collapsed cage structure gave a notably better match with the experiment. It could subsequently be concluded from the VCD analysis, that the octaazacryptands prefer a collapsed structure, which is not predicted by density functional theory (DFT) calculations as the global minimum structures. These findings are expected to have consequences also for future studies on inclusion complexes of such azacryptands.

Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach

Vibrational circular dichroism (VCD) spectroscopy has become an important part of the (stereo-)chemists' toolbox as a reliable method for the determination of absolute configurations. Being the chiroptical version of infrared spectroscopy, it has also been recognized as being very sensitive to conformational changes and intermolecular interactions. This sensitivity originates from the fact that the VCD spectra of individual conformers are often more different than their IR spectra, so that changes in conformational distributions or band positions and intensities become more pronounced. What is an advantage for studies focussing on intermolecular interactions can, however, quickly turn into a major obstacle during AC determinations: solute-solvent interactions can have a strong influence on spectral signatures and they must be accurately treated when simulating VCD and IR spectra. In this perspective, we showcase selected examples which exhibit particularly pronounced solvent effects. It is demonstrated that it is typically sufficient to model solute-solvent interactions by placing single solvent molecules near hydrogen bonding sites of the solute and subsequently use the optimized structures for spectra simulations. This micro-solvation approach works reasonably well for medium-sized, not too conformationally flexible molecules. We thus also discuss its limitations and outline the next steps that method development needs to take in order to further improve the workflows for VCD spectra predictions.

Isolated and solvated chiroptical behavior in conformationally flexible butanamines

The nonresonant optical activity of two highly flexible aliphatic amines, (2R)-3-methyl-2-butanamine (R-MBA) and (2R)-(3,3)-dimethyl-2-butanamine (R-DMBA), has been probed under isolated and solvated conditions to examine the roles of conformational isomerism and to explore the influence of extrinsic perturbations. The optical rotatory dispersion (ORD) measured in six solvents presented uniformly negative rotatory powers over the 320-590 nm region, with the long-wavelength magnitude of chiroptical response growing nearly monotonically as the dielectric constant of the surroundings diminished. The intrinsic specific optical rotation,α λ T (in deg dm-1 [g/mL]-1 ), extracted for ambient vapor-phase samples of R-MBA [-11.031(98) and -2.29 (11)] and R-DMBA [-9.434 (72) and -1.350 (48)] at 355 and 633 nm were best reproduced by counterintuitive solvents of high polarity (yet low polarizability) like acetonitrile and methanol. Attempts to interpret observed spectral signatures quantitatively relied on the linear-response frameworks of density-functional theory (B3LYP, cam-B3LYP, and dispersion-corrected analogs) and coupled-cluster theory (CCSD), with variants of the polarizable continuum model (PCM) deployed to account for the effects of implicit solvation. Building on the identification of several low-lying equilibrium geometries (nine for R-MBA and three for R-DMBA), ensemble-averaged ORD profiles were calculated at T = 300 K by means of the independent-conformer ansatz, which enabled response properties predicted for the optimized structure of each isomer to be combined through Boltzmann-weighted population fractions derived from corresponding relative internal-energy or free-energy values, the latter of which stemmed from composite CBS-APNO and G4 analyses. Although reasonable accord between theory and experiment was realized for the isolated (vapor-phase) species, the solution-phase results were less satisfactory and tended to degrade progressively as the solvent polarity increased. These trends were attributed to solvent-mediated changes in structural parameters and energy metrics for the transition states that separate and putatively isolate the equilibrium conformations supported by the ground electronic potential-energy surface, with the resulting displacement of barrier locations and/or decrease of barrier heights compromising the underlying premise of the independent-conformer ansatz.

Diaryl-Pyrano-Chromenes Atropisomers: Stereodynamics and Conformational Studies

The dynamic scenario of di-aryls-pyrano-chromenes was investigated using DFT calculations. The symmetry of the chromene scaffold and the presence of two ortho-substituted aryls substituents can generate two syn/anti diastereoisomers and conformational enantiomers with different rotational barriers. The relative conformations and configurations were derived using NOESY-1D experiments. Depending on the energies related to the conformational exchange, the experimental energy barriers were determined through Dynamic NMR, Dynamic HPLC or kinetic studies. The atropisomeric pairs were resolved in the latter scenario, and their absolute configuration was assigned using the ECD/TD-DFT method.

Can the absolute configuration of cyclic peptides be determined with vibrational circular dichroism?

Cyclic peptides show a wide range of biological activities, among others as antibacterial agents. These peptides are often large and flexible with multiple chiral centers. The determination of the stereochemistry of molecules with multiple chiral centers is a challenging and important task in drug development. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) can distinguish between different stereoisomers. The absolute configuration (AC) of a stereoisomer can be determined by comparing its experimental spectra to computed spectra of stereoisomers with known AC. In this way, the AC of rigid molecules with up to seven chiral centers has been assigned (Bogaerts et al., Phys. Chem. Chem. Phys., 2020, 22, 18014). The question arises whether this is possible with more conformationally flexible molecules such as cyclic peptides. We here investigate to what extent the AC of cyclic peptides can be determined with VCD. More specifically, we investigate the maximum number of chiral centers a cyclic peptide can have in order to be able to unambiguously assign the AC with VCD. We present experimental and computed IR and VCD spectra for a series of eight tetrapeptides and hexapeptides with two, three and four chiral centers. We use our recently developed computational protocol with a conformational search based on sampling with meta-dynamics. We use visual inspection to compare the computed spectra of different stereoisomers with an experimental spectrum of the corresponding cyclic peptide with known AC. We find that the AC of the investigated cyclic peptides with two chiral centers can be unambiguously assigned with VCD. This is however not possible for all of the cyclic peptides with three chiral centers and for none of those with four chiral centers. At best, one can limit the number of possible stereoisomers in those cases. Our work shows that other techniques are needed to assign the AC of cyclic peptides with three or more chiral centers. Our study also constitutes a warning that the spectra of all stereoisomers should be computed before attempting to match to an experimental spectrum, to avoid an accidental erroneous match.

Atropostatin: Design and Total Synthesis of an Atropisomeric Lactone–Atorvastatin Prodrug

Atorvastatins play an important role in the inhibition of HMG-CoA reductase, an enzyme present in the liver that takes part in the biosynthesis of cholesterol. In this article, we report the total synthesis of a lactone–atorvastatin prodrug with additional atropisomeric features. Conformational and experimental studies of model compounds were designed to test the stability of the chiral axis. Docking calculations were performed to evaluate the constant inhibition of a library of atorvastatins. Full synthesis of the best candidate was achieved and thermally stable atropisomeric lactone–atorvastatin was obtained. The absolute configuration of the chiral axis of the atropisomers was assigned by means of chiroptical ECD spectroscopy coupled with TD-DFT calculations.

Matrix-isolation and cryosolution-VCD spectra of α-pinene as benchmark for anharmonic vibrational spectra calculations

The inclusion of anharmonicity in vibrational spectral analyis remains associated to small molecular systems with up to a dozen of atoms, with half a dozen of non-hydrogen atoms, typically thesize of propylene oxide. One may see two reasons for this: first of all, larger systems are often thought to be computationally too demanding (high computational costs) for a full anharmonic vibrational analysis. Second, the identification of resonances and their correction is often considered something only expert theoreticians could address because of the lack of unequivocal criteria. In this contribution, we illustrate that resonances can indeed become a complex problem, which can be handled almost transparently thanks to recent advances in vibrational perturbation theory (VPT2). The study also emphasizes the importance and the central role played by experiment in benchmarking novel theoretical approaches. In fact, we herein provide the currently highest resolution VCD spectra available for α- and β-pinene obtained under matrix-isolation conditions and in liquid Xenon as solvent. They are interpreted by VPT2 with novel tests for the identification of resonances. Hence, the study demonstrates the mutual stimulation of advances in both experimental techniques and computational models.

Dynamic Stereochemistry of a Biphenyl-Bisprolineamide Model Catalyst and its Imidazolidinone Intermediates

In this study, we characterize the dynamic stereochemistry of a biphenyl-2,2'-bis(proline amide) catalyst in chloroform and DMSO as representative weakly and strongly hydrogen bonding solvents. Using vibrational circular dichroism (VCD) spectroscopy and density functional theory (DFT) based spectra calculations, we show that the preferred axial stereochemistry of the catalyst is determined by solute-solvent interactions. Explicitly considering solvation with DMSO molecules is found to be essential to correctly predict the conformational preferences of the catalyst. Furthermore, we investigate the stereochemistry of the corresponding enamines and imidazolidinones that are formed upon reaction with isovaleraldehyde. The enamines are found to rapidly convert to endo-imidazolidinones and the thermodynamically favored exo-imidazolidinones are formed only slowly. The present study demonstrates that the stereochemistry of these imidazolidinones can be deduced directly from the VCD spectra analysis without any further detailed analysis of NMR spectra. Hence, we herein exemplify the use of VCD spectroscopy for an in situ characterization of intermediates relevant in asymmetric catalysts.

Solvation of the Boc-Val-Phe-nPr peptide characterized by VCD spectroscopy and DFT calculations

The conformational preferences of peptides are strongly determined by hydrogen bonding interactions. Intermolecular solute-solvent interactions compete with intramolecular interactions, which typically stabilize the secondary structure of the peptide. The analysis of vibrational circular dichroism (VCD) spectra can give insights into solvation-induced changes in the conformational distribution of small peptides. Here we describe the VCD spectroscopic characterization of the model peptide Boc-Val-Phe-nPr in chloroform as representative for a weakly interacting solvent and dimethyl sulfoxide (DMSO-d₆) as a strongly hydrogen bonding solvent. We show that the conformational preferences of the peptide in chloroform are well-described by the computationally predicted distribution of the isolated molecule assuming only implicit solvation effects through a continuum solvation model. In order to simulate the spectra recorded in DMSO-d₆, solvation was accounted for explicitly by computed microsolvated structures containing one to three solvent molecules. A good match of the computed spectra with the experimental data is obtained by this method. Comparing the conformational distributions in deuterated chloroform-d₁ and DMSO-d₆, structures with intramolecular hydrogen bonds such as the (δ,δ)-conformer family contribute to the conformational distribution only when there is no strong interaction with the solvent. This is in contrast to the results for the related Boc-Pro-Phe-nPr studied before, for which the intramolecular interaction was found to persist in DMSO-d₆. Furthermore, we discuss the influence of hydrogen bonding to different numbers of solvent molecules on the spectral signatures and show that the structure of the peptide in DMSO-d₆ is best described as a mixture of twofold-solvated (δ,β)- and threefold-solvated (β,β)-conformers.

Proton-coupled electron transfer of catechin in tea wine: the enhanced mechanism of anti-oxidative capacity

Tea wine is a Chinese traditional alcoholic drink made by cereal and tea leaves. It is rich in tea polyphenols, caffeine, amino acids, and protons and possesses various healthcare functions. In this work, electrochemical methods, as well as density functional theory (DFT) calculations, were adopted to reveal the proton-coupled electron-transfer process of catechin in tea wine. The electrochemical results showed that the catechin preferred hydrogen-bonding with ethanol and formed molecular clusters. Thus, the direct electron-transfer process of catechin changed to proton-coupled electron transfer. This procedure reduced the energy barrier of the redox reaction and enhanced the anti-oxidative capacity. Subsequently, DFT calculations were employed to explore the bond length, bond energy, and HOMO-LUMO energy gap of catechin, which confirmed the above-mentioned mechanism. Our work offers some positive value for the scientific promotion of traditional food and a greater understanding of the health mechanisms in terms of chemistry.

The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution

We present a computational study of vibrational circular dichroism (VCD) in solutions of (S)-lactic acid, relying on ab initio molecular dynamics (AIMD) and full solvation with bulk water. We discuss the effect of the hydrogen bond network on the aggregation behaviour of the acid: while aggregates of the solute represent conditions encountered in a weakly interacting solvent, the presence of water drastically interferes with the clusters - more strongly than originally anticipated. For both scenarios we computed the VCD spectra by means of nuclear velocity perturbation theory (NVPT). The comparison with experimental data allows us to establish a VCD-structure relationship that includes the solvent network around the chiral solute. We suggest that fundamental modes with strong polarisation such as the carbonyl stretching vibration can borrow VCD from the chirally restructured solvent cage, which extends the common explanatory models of VCD generation in aqueous solution.

Optical Activity in Saturated Cyclic Amines: Untangling the Roles of Nitrogen-Inversion and Ring-Puckering Dynamics

The dispersive optical activity of two saturated cyclic amines, (R)-2-methylpyrrolidine (R-2MPY) and (S)-2-methylpiperidine (S-2MPI), has been interrogated under isolated and solvated conditions to elucidate the roles of large-amplitude motion associated with nitrogen-center inversion and ring-puckering dynamics. Experimental optical rotatory dispersion profiles were almost mirror images of one another and displayed parallel solvent dependencies. Quantum-chemical analyses built on density-functional and coupled-cluster methods revealed four low-lying conformers for each molecule, which are distinguished by axial/equatorial orientations of their amino hydrogens and methyl substituents. Chiroptical signatures predicted for these species were combined through an independent-conformer ansatz to simulate the ensemble-averaged response, with a polarizable continuum model (PCM) being used to treat implicit solute-solvent interactions. The intrinsic behavior observed for isolated (gaseous) R-2MPY and S-2MPI was reproduced best by merging coupled-cluster (CCSD) estimates of rotatory powers with thermal population fractions deduced from complete basis set (CBS-APNO) free-energy calculations. Although prior claims of sizable chiroptical contributions arising from helically twisted (chiral) heterocyclic frameworks could be discounted, less satisfactory agreement between experiment and theory was realized for solution phases. Response properties sustained modest isomer-dependent changes in the presence of PCM solvation, but the corresponding energy metrics showed systematic trends, whereby structures having larger electric-dipole moments were stabilized preferentially in media of high polarity. Despite the fact that R-2MPY conformations were predicted to undergo a progressive reordering of their relative energies across the six solvents of interest, S-2MPI was found to exhibit more pronounced solvent-induced perturbations at long wavelengths (viz., in regions far removed from electronic resonances). Experimental results are discussed in terms of the distinct ring-puckering mechanisms for R-2MPY and S-2MPI, which are expected to be dominated by hindered pseudorotation among envelope/twist motifs and semi-inversion between chairlike antipodes, respectively.

Modeling Complex Solvent Effects on the Optical Rotation of Chiral Molecules: A Combined Molecular Dynamics and Density Functional Theory Study

The challenge of assigning the absolute stereochemical configuration to a chiral compound can be overcome via accurate ab initio predictions of optical rotation, a sensitive molecular property that is further complicated by solvent effects. The solvent's "chiral imprint"-the transfer of the chirality from the solute to the surrounding achiral solvent-is explored here using conformational averaging and time-dependent density-functional theory. These complex solvent effects are taken into account via simple averaging over a molecular dynamics trajectory together with the explicit quantum mechanical consideration of the solvent molecules within the solute's cybotactic region and implicit modeling of the bulk solvent. We consider several axes along which the system's optical rotation varies, including the sampling of the dynamical trajectory, the quality of the one-electron basis set, and the use of continuum solvent models to account for bulk effects.

Natural and magnetic circular dichroism spectra of nucleosides: effect of the dynamics and environment

Chiroptical spectroscopic methods are excellent tools to study structure and interactions of biomolecules. However, their sensitivity to different structural aspects varies. To understand the dependence of absorption, electronic and magnetic circular dichroism (ECD, MCD) intensities on the structure, dynamics and environment, we measured and simulated spectra of nucleosides and other nucleic acid model components. The conformation space was explored by molecular dynamics (MD), the electronic spectra were generated using time dependent density functional theory (TDDFT). The sum over state (SOS) method was employed for MCD. The results show that accounting for the dynamics is crucial for reproduction of the experiment. While unpolarized absorption spectroscopy is relatively indifferent, ECD reflects the conformation and geometry dispersion more. MCD spectra provide variable response dependent on the wavelength and structural change. In general, MCD samples the structure more locally than ECD. Simple computational tests suggest that the optical spectroscopies coupled with the computational tools provide useful information about nucleic acid components, including base pairing and stacking.

Basis set dependence of S[double bond, length as m-dash]O stretching frequencies and its consequences for IR and VCD spectra predictions

Benchmarking functionals and basis sets for the computational prediction of molecular properties is usually done on very small model systems. Larger organic molecules containing heavier second row atoms are not the typical model structures. We herein present the first survey of basis sets and functionals for the prediction of the IR and VCD spectra of chiral tosylates and sulfinates as we noted drastic deviations between computed harmonic frequencies obtained at B3LYP/6-311++G(2d,p) level of theory and those observed in experimental solution phase IR and VCD spectra. We show that the harmonic frequencies of the asymmetric and symmetric S[double bond, length as m-dash]O stretching modes of tosylates are predicted at significantly too low vibrational frequencies if the employed basis set does not provide higher order polarization functions. The results of our benchmarks show that at least the 6-311G(3df,2dp) basis (or equivalent Dunning and Ahlrichs variants) should be used.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Solvation and the secondary structure of a proline-containing dipeptide: insights from VCD spectroscopy

In this study we investigate the IR and VCD spectra of the diastereomeric dipeptide Boc-Pro-Phe-(n-propyl) 1 in chloroform-d₁ (CDCl₃) and the strongly hydrogen bonding solvent dimethylsulfoxide-d₆ (DMSO-d₆). From comparison of the experimental spectra, the amide II spectral region is identified as marker signature for the stereochemistry of the dipeptide: the homochiral LL-1 features a (+/-)-pattern in the amide II region of the VCD spectrum, while the amide II signature of the diastereomer LD-1 is inverted. Computational analysis of the IR and VCD spectra of LL-1 reveals that the experimentally observed amide II signature is characteristic for a β_I-turn structure of the peptide. Likewise, the inverted pattern found for LD-1 arises from a β_II-turn structure of the dipeptide. Following a micro-solvation approach, the experimental spectra recorded in DMSO-d₆ are computationally well reproduced by considering only a single solvent molecule in a hydrogen bond with N-H groups. Considering a second solvent molecule, which would lead to a cleavage of intramolecular hydrogen bonds in 1, is found to give a significantly worse match with the experiment. Hence, the detailed computational analysis of the spectra of LL- and LD-1 recorded in DMSO-d₆ confirms that the intramolecular hydrogen bonding pattern, that stabilizes the β-turns and other conformations of LL- and LD-1 in apolar solvents, remains intact. Our findings also show that it is essential to consider solvation explicitly in the analysis of the IR and VCD spectra of dipeptides in strongly hydrogen bonding solvents. As the solute-solvent interactions affect both conformational preferences and spectral signatures, it is also demonstrated that this inclusion of solvent molecules cannot be circumvented by applying fitting procedures to non-solvated structures.

Approximate versus Exact Embedding for Chiroptical Properties: Reconsidering Failures in Potential and Response

We investigate the suitability of subsystem time-dependent density-functional theory (sTDDFT) for describing chiroptical properties with a focus on optical rotation parameters. Our starting point is a new implementation of the recently proposed projection-based, coupled frozen-density embedding (FDEc) framework. We adapt the generalized, non-Hermitian formulation of TDDFT and derive corresponding expressions for regular and damped response properties from subsystem TDDFT. We verify that our implementation of this "exact" formulation allows to reproduce supermolecular results of electronic circular dichroism (ECD) spectra, of optical rotatory dispersion, and of polarizabilities. We present a systematic test of the main approximations typically introduced in practical frozen-density embedding (FDE) calculations of response properties: (i) the use of approximate nonadditive kinetic-energy (NAKE) functionals, which can be avoided through projection techniques, (ii) the use of monomer (subsystem) basis sets rather than supersystem basis sets, and (iii) the neglect of intersubsystem response coupling within the so-called uncoupled FDE (or FDEu) approximation. While approximation (i) is known to generally lead to large errors for covalently bound subsystems, we present cases in which either the basis set or the coupling step are similarly or even (much) more important. In particular, we explicitly demonstrate by comparison to a fully coupled calculation that missing intersubsystem response couplings are responsible for the failure of FDE reported in a previous study [ J. Chem. Theory Comput. 2015, 11, 5305-5315]. We show that good agreement with reference results can be obtained in this case even with standard NAKE approximations for the FDE potentials and efficient monomer basis sets, making calculations for larger systems well accessible.

Chiroptical Sensing: A Conceptual Introduction

Chiroptical responses have been an essential tool over the last decades for chemical structural elucidation due to their exceptional sensitivity to geometry and intermolecular interactions. In recent times, there has been an increasing interest in the search for more efficient sensing by the rational design of tailored chiroptical systems. In this review article, advances made in chiroptical systems towards their implementation in sensing applications are summarized. Strategies to generate chiroptical responses are illustrated. Theoretical approaches to assist in the design of these systems are discussed. The development of efficient chiroptical reporters in different states of matter, essential for the implementation in sensing devises, is reviewed. In the last part, remarkable examples of chiroptical sensing applications are highlighted.

More complex, less complicated? Explicit solvation of hydroxyl groups for the analysis of VCD spectra

With increasing size of the molecules, hydrogen bonding induced solvent effects on the IR and VCD spectra become more negligible.

How many solvent molecules are required to solvate chiral 1,2-diols with hydrogen bonding solvents? A VCD spectroscopic study

A VCD spectroscopic analysis of selected model systems for solute–solvent interactions of chiral diols with hydrogen bonding solvents DMSO and ACN.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.

Absolute Configurations of Synthetic Molecular Scaffolds from Vibrational CD Spectroscopy

Vibrational circular dichroism (VCD) spectroscopy is one of the most powerful techniques for the determination of absolute configurations (AC), as it does not require any specific UV/vis chromophores, no chemical derivatization, and no growth of suitable crystals. In the past decade, it has become increasingly recognized by chemists from various fields of synthetic chemistry such as total synthesis and drug discovery as well as from developers of asymmetric catalysts. This perspective article gives an overview about the most important experimental aspects of a VCD-based AC determination and explains the theoretical analysis. The comparison of experimental and computational spectra that leads to the final conclusion about the AC of the target molecules is described. In addition, the review summarizes unique VCD studies carried out in the period 2008-2018 that focus on the determination of unknown ACs of new compounds, which were obtained in its enantiopure form either through direct asymmetric synthesis or chiral chromatography.

A combined experimental and theoretical study of optical rotatory dispersion for (R)-glycidyl methyl ether in aqueous solution

The dispersive optical activity for aqueous solutions of non-rigid (R)-glycidyl methyl ether (R-GME) has been explored synergistically from experimental and theoretical perspectives. Density functional theory analyses performed with the polarizable continuum model for implicit solvation identified nine low-lying stable conformers that are interconverted by rotation about two large-amplitude torsional coordinates. The antagonistic chiroptical signatures predicted for these structural isomers were averaged under a Boltzmann-weighting ansatz to estimate the behavior expected for a thermally equilibrated ensemble. This led to optical rotatory dispersion profiles that reproduced the overall shape of observations but failed to achieve uniform agreement with measured specific-rotation values even when anharmonic vibrational corrections were applied. A mixed QM/FQ paradigm, whereby quantum-mechanical (QM) calculations of optical activity were combined with classical molecular dynamics simulations of explicit solvation that included mutual-polarization effects by means of fluctuating charges (FQ), was enlisted to elucidate the microsolvation environment and gauge its impact upon conformer distributions and response properties. Although quantitative accord with experiments remained elusive, this approach revealed strong variations in the magnitude and sign of rotatory powers for R-GME as the configuration of surrounding water molecules evolved, thereby highlighting the inherently dynamical nature of the solvated chiroptical response, calling into question the validity of "static" descriptions based on the presumption of distinct energy minima, and giving insight into the inherent complexity posed by the modeling of such properties for solvated systems.

On the Absolute Stereochemistry of Tolterodine: A Circular Dichroism Study

Tolterodine (1) is a potent muscarinic receptor antagonist used in the treatment of overactive urinary bladder (OAB) syndrome. Tolterodine is chiral and it was patented, and is currently marketed, as the l-tartrate salt of the (R)-enantiomer. However, the existing literature does not offer an ultimate proof of a stereoselective mode of action of 1. A second open stereochemical issue concerns the absolute configuration (AC) of 1. Neither the original patents nor subsequent studies have established the AC of 1 in an unambiguous way, although the AC of the l-tartrate salt of 1 was assigned by X-ray diffractometry. Finally, neither electronic nor vibrational circular dichroism (ECD and VCD) spectra of 1 are reported so far. We performed a thorough ECD/VCD study of 1 in different solvents and at variable temperatures. Solvent and temperature dependence highlighted the existence of moderate flexibility which was confirmed by molecular modelling. ECD calculations with time-dependent density functional theory (TDDFT) accurately reproduced the experimental spectra and allowed us to confirm the AC of 1 in an independent way.

Insight into vibrational circular dichroism of proteins by density functional modeling

Vibrational circular dichroism (VCD) spectroscopy is an excellent method to determine the secondary structure of proteins in solution. Comparison of experimental spectra with quantum-chemical simulations represents a convenient and objective way to extract information on the structure. This has been difficult for such large molecules where approximate theoretical models have to be used. In the present study we applied the Cartesian-coordinate based tensor transfer (CCT) making it possible to extend the density functional theory (DFT) and model spectral intensities of large globular proteins nearly at quantum-chemical precision. Indeed, comparison with experiment provided a better understanding of the dependence of VCD spectral shapes on the geometry, their sensitivity to fine structural details and interactions with the environment. On a model set of globular proteins the simulated spectra correlated well with experimental data and revealed which structural information can (and cannot) be obtained from this kind of spectroscopy. Although the VCD technique has been regarded as being rather insensitive to side-chain variations, we found that the spectra of human and hen lysozyme differing by a few amino acids only are quite distinct. This has been explained by long-distance coupling of the amide vibrations. Likewise, the modeling reproduced some spectral changes caused by protein deuteration even when the protein structure was conserved.

The determination of enantiomer composition of 1-((3-chlorophenyl)-(phenyl)methyl) amine and 1-((3-chlorophenyl)(phenyl)-methyl) urea (Galodif) by NMR spectroscopy, chiral HPLC, and polarimetry

For the first time, a method for enantiomer resolution of the anticonvulsant Galodif (1-((3-chlorophenyl)(phenyl)methyl) urea) by chiral HPLC was developed, whereas the enantiomeric composition of 1-((3-chlorophenyl)(phenyl)methyl) amine-precursor in Galodif synthesis-cannot be resolved by this method. However, starting 1-((3-chlorophenyl)(phenyl)methyl) amine quantitatively forms diastereomeric N-((3-chlorophenyl)(phenyl)methyl)-1-camphorsulfonamides in reaction with chiral (1R)-(+)- or (1S)-(-)-camphor-10-sulfonyl chlorides. The diastereomeric ratio of obtained camphorsulfonamides can be easily determined by NMR ¹ H and ¹³ C spectroscopy. The DFT calculations of specific rotation of Galodif enantiomers showed good agreement with experimental data. The absolute configuration of enantiomers was proposed for the first time.

Polyoxygenated cyathane diterpenoids from the mushroom Cyathus africanus, and their neurotrophic and anti-neuroinflammatory activities

In a previous study, we reported ten new polyoxygenated cyathane diterpenoids, neocyathins A-J, and their anti-neuroinflammatory effects from the liquid culture of the medicinal Basidiomycete Cyathus africanus. In the present study, eight new highly polyoxygenated cyathane diterpenoids, named neocyathins K-R (1-8), were isolated from the solid culture of C. africanus cultivated on cooked rice, together with three known congeners (9-11). The structures and the absolute configurations of the new compounds were elucidated through comprehensive NMR and HRESIMS spectroscopic data, electronic circular dichroism (ECD) data, and chemical conversion. Compounds 1 and 2 represent the first reported naturally occurring compounds with 4,9-seco-cyathane carbon skeleton incorporating an unprecedented medium-sized 9/7 fused ring system, while the 3,4-seco-cyathane derivative (3) was isolated from Cyathus species for the first time. All compounds were evaluated for their neurotrophic and anti-neuroinflammatory activity. All the isolates at 1-25 μM displayed differential nerve growth factor (NGF)-induced neurite outgrowth-promoting activity in PC-12 cells, while one of the compounds, allocyathin B₂ (11), inhibited NO production in lipopolysaccharide (LPS)-stimulated microglia BV-2 cells. In addition, molecular docking studies showed that compound 11 generated interactions with the inducible nitric oxide synthase (iNOS) protein.

Solvent-sensitive signs and magnitudes of circularly polarised luminescence and circular dichroism spectra: probing two phenanthrenes as emitters endowed with BINOL derivatives

A C₂-symmeric binaphthyl framework bearing phenanthrene as the emitter exhibited circularly polarised luminescence (CPL) in dilute solutions.

Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution

A polarizable QM/MM approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems is proposed and applied to aqueous solutions of (l)-methyl lactate and (S)-glycidol.

Specific optical rotation is a versatile tool for the identification of critical micelle concentration and micellar growth of tartaric acid-based diastereomeric amphiphiles

Four novel tartaric acid-based diastereomeric chiral amphiphiles, two being enantiomers of the other two, have been synthesized and investigated using chiroptical spectroscopic methods, along with tensiometry and dynamic light scattering experiments. We found that an inflection point in specific optical rotation (SOR) values at ~0.32 mM corresponds to the critical micelle concentration (CMC). The increase in magnitude of SOR values beyond CMC corresponds to the growth of aggregates. For enantiomers, oppositely signed SOR values were observed, ruling out the possibility for the presence of aggregation size mediated artefacts. SOR values did not exhibit concentration dependence for a chiral tartaric acid based non-aggregating analogue further establishing the absence of artefacts or anomalous interaction of tartaric acid based head group with solvent. Electronic circular dichroism spectra showed no significant changes in band positions or intensities with concentration. Due to the requirement for higher concentrations (~200 mM) needed to obtain vibrational circular dichroism spectra, these measurements are not found to be useful for studying concentration dependent properties of chiral amphiphiles.

Chirality sensing of bioactive compounds with amino alcohol unit via circular dichroism

The aim of the present work was to test various chiroptical techniques, including in particular the in situ dirhodium methodology, to assign the absolute configuration of 1,2- and 1,3-amino alcohols. As models, we selected mainly compounds that have both an additional strongly absorbing and interfering chromophoric system and application in medicinal chemistry. Determination of the absolute configuration (AC) of the tested molecules such as cinchona alkaloids, Tamiflu, and others was carried out using a combination of electronic and vibrational circular dichroism (ECD, VCD) spectroscopy. It has been demonstrated that both 1,2- and 1,3-aminol moieties are subject to the same sector rule correlating stereostructure of formed Rh₂ -complex with chiroptical properties, and that the changes in the position of the stereogenic center do not affect its proper use.

The shape of the electronic circular dichroism spectrum of (2,6-dimethylphenyl)(phenyl)methanol: interplay between conformational equilibria and vibronic effects

Analysis of the interplay between conformational equilibria, solvent effects and vibronic contributions in the ECD spectra.

Bimodal concentration-dependent reactivity pattern of a glycosyl donor: Is the solution structure involved?

Changes in concentration (0.001-0.1 M) of an arabinofuranosyl donor (1) have been shown to modulate the temperature T at which activation of 1 occurs (from -23 °C to +7 °C), the reaction time (from 1.5 h to 3 days) and the yield of the disaccharide formed (from 14% to 82%). At concentrations exceeding 0.01 M, these parameters, as well as the specific optical rotation of the solution of 1, virtually do not depend on concentration suggesting formation of reacting species (supramers) of glycosyl donor with similar structures, hence reactivities, but considerably different from those formed in more dilute solutions. The found critical concentration (0.01 M) separates two concentration ranges of reaction solutions corresponding to two types of solution structure that are featured by the presence of fundamentally different supramers of glycosyl donor, which have distinct chemical properties. These results allow a fresh look at the problems of reactivity of chemical compounds and selectivity of the reactions in which they participate.

Optical Rotation Calculations for a Set of Pyrrole Compounds

Optical rotation of 14 molecules containing the pyrrole group is calculated by employing both time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional and the second-order approximate coupled-cluster singles and doubles (CC2) method. All optical rotations have been provided using the aug-cc-pVDZ basis set at λ = 589 nm. The two methods predict similar results for both sign and magnitude for the optical rotation of all molecules. The obtained signs are consistent with experiments as well, although several conformers for four molecules needed to be studied to reproduce the experimental sign. We have also calculated excitation energies and rotatory strengths for the six lowest lying electronic transitions for several conformers of the two smallest molecules and found that each rotatory strength has various contributions for each conformer which can cause different optical rotations for different conformers of a molecule. Our results illustrate that both methods are able to reproduce the experimental optical rotations, and that the CAM-B3LYP functional, the least computationally expensive method used here, is an applicable and reliable method to predict the optical rotation for these molecules in line with previous studies.