Automated Quantum Chemistry-Based Calculation of Optical Rotation for Large Flexible Molecules

Emerging Conformational-Analysis Protocols from the RTCONF55-16K Reaction Thermochemistry Conformational Benchmark Set

RTCONF55-16K is a new, reactive conformational data set based on cost-efficient methods to assess different conformational analysis protocols. Our reference calculations underpinned the accuracy of the CENSO (Grimme et al. J. Phys. Chem. A, 2021, 125, 4039) procedure and resulted in alternative recipes with different cost-accuracy compromises. Our general-purpose and economical protocols (CENSO-light and zero, respectively) were found to be 10-30 times faster than the original algorithm, adding only 0.4-0.7 kcal/mol absolute error to the relative free energy estimates.

Study of optical rotation based on the molecular structure in fused oligomers of macrocycles

We designed a unique oligomer form in which several helically twisted macrocycles (M- or P-helicity) are arranged through fusion. We investigated the optical rotation of a series of fused oligomers of macrocycles with a difference in the number and arrangement of elements associated with point-chiral auxiliary. Some oligomers cooperatively attained a situation where an identical sense of twisting was preferred throughout the entire molecule. On the basis of these results, we estimated diastereomeric excess induced in each oligomer. We revealed that the molar optical rotation per element was modulated with a rotational angle between elements: an increase via 0° rotational arrangement, a decrease via 180° rotational arrangement, or a decrease via cyclic arrangement. Alternatively, for other oligomers in which several diastereomeric conformers coexist, we uniquely attempted to consider the optical rotation based on the molecular structure through the assessment of a change ratio of the absorption dissymmetry factor before and after complexation with an achiral guest. We found that molar optical rotation could be different based on the arrangement, although actual measured values were similar.

Reciprocal polarization imaging of optical activity in reflection

We present reciprocal polarization imaging for the optical activity of chiral media in reflection geometry. The method is based on the reciprocal polar decomposition of backscattering Mueller matrices accounting for the reciprocity of light waves in forward and backward scattering paths. Anisotropic depolarization is introduced to gain sensitivity to optical activity in backscattering. Experiments with glucose solutions show that while the Lu-Chipman decomposition of the backscattering Mueller matrices produces erroneous results, reciprocal polarization imaging correctly retrieves the optical activity of chiral media. The recovered optical rotation agrees with that obtained in the forward geometry and increases linearly with the concentration and thickness of the chiral media. The potential for in vivo glucose monitoring based on optical activity sensing using reciprocal polarization imaging is then discussed.

Supramolecular Chirogenesis in a Sterically Hindered Porphyrin: A Critical Theoretical Analysis

The determination of molecular stereochemistry and absolute configuration is an important part of modern chemistry, pharmacology, and biology. Electronic circular dichroism (ECD) spectroscopy is a widely used tool for chirality assignment, especially with porphyrin macrocycles employed as reporter chromophores. However, the mechanisms of induced ECD in porphyrin complexes are yet to be comprehensively rationalized. In this work, the ECD spectra of a sterically hindered hexa-cationic porphyrin with two camphorsulfonic acids in dichloromethane and chloroform were experimentally measured and computationally analyzed. The influence of geometric factors such as the position of chiral guest molecules, distortion of the porphyrin macrocycle, and orientation of aromatic and non-aromatic peripheral substituents on the ECD spectra was theoretically studied. Various potential pitfalls, such as a lack of significant conformations and accidental agreement of experimental and simulated spectra, are considered and discussed.

Pseudo-Enantiomeric Paraphaeolactones and Their Biosynthetic Derivatives from Paraphaeosphaeria sp. KT4192: A Proposition of a Favorskii Rearrangement for Their Biosynthesis

Paraphaeolactones A1, A2, B1, and B2 (1-4, respectively), known arthropsadiol D (5), massariphenone (6) and its positional isomer 7, and massarilactones E (8) and G (9) were isolated from the culture broth of Paraphaeosphaeria sp. KT4192. Although the structural resemblance between 1 and 2 implies that these comprised a diastereomeric pair at the C-2 stereogenic center, electronic circular dichroism (ECD) spectral analyses revealed that they were pseudo-enantiomers possessing the common (2R)-configuration. Paraphaeolactones B1 and B2 (3 and 4) were the derivatives of 2, which equipped the 3-(1-hydroxy-2-oxopropyl)-4-methylcatechol moiety via an acetal bond at C-10. The relative configurations of their acetal carbons were elucidated by NOE experiments, and those of C-8' were deduced independently by ECD spectral analysis. The present study disclosed that 1-5, 8, and 9 contain a methylcyclohexene substructure with the same absolute configuration. This prompted us to reinvestigate the absolute configurations of known structurally related fungal metabolites, allowing us to conclude that the methylcyclohexene moieties of these natural products have the same absolute configuration despite the variety of configurations of other stereogenic centers. The plausible biosynthetic routes for 1-9 are discussed on the basis of the above conclusion. We propose a Favorskii rearrangement as the key transformation for biosyntheses of 1-4.

Theoretical-Computational Modeling of CD Spectra of Aqueous Monosaccharides by Means of Molecular Dynamics Simulations and Perturbed Matrix Method

The electronic circular dichroism (ECD) spectra of aqueous d-glucose and d-galactose were modeled using a theoretical-computational approach combining molecular dynamics (MD) simulations and perturbed matrix method (PMM) calculations, hereafter termed MD-PMM. The experimental spectra were reproduced with a satisfactory accuracy, confirming the good performances of MD-PMM in modeling different spectral features in complex atomic-molecular systems, as already reported in previous studies. The underlying strategy of the method was to perform a preliminary long timescale MD simulation of the chromophore followed by the extraction of the relevant conformations through essential dynamics analysis. On this (limited) number of relevant conformations, the ECD spectrum was calculated via the PMM approach. This study showed that MD-PMM was able to reproduce the essential features of the ECD spectrum (i.e., the position, the intensity, and the shape of the bands) of d-glucose and d-galactose while avoiding the rather computationally expensive aspects, which were demonstrated to be important for the final outcome, such as (i) the use of a large number of chromophore conformations; (ii) the inclusion of quantum vibronic coupling; and (iii) the inclusion of explicit solvent molecules interacting with the chromophore atoms within the chromophore itself (e.g., via hydrogen bonds).

Epimeric Mixture Analysis and Absolute Configuration Determination Using an Integrated Spectroscopic and Computational Approach-A Case Study of Two Epimers of 6-Hydroxyhippeastidine

Structural elucidation has always been challenging, and misassignment remains a stringent issue in the field of natural products. The growing interest in discovering unknown, complex natural structures accompanies the increasing awareness concerning misassignments in the community. The combination of various spectroscopic methods with molecular modeling has gained popularity in recent years. In this work, we demonstrated, for the first time, its power to fully elucidate the 2-dimensional and 3-dimensional structures of two epimers in an epimeric mixture of 6-hydroxyhippeastidine. DFT calculation of chemical shifts was first performed to assist the assignment of planar structures. Furthermore, relative and absolute configurations were established by three different ways of computer-assisted structure elucidation (CASE) coupled with ORD/ECD/VCD spectroscopies. In addition, the significant added value of OR/ORD computations to relative and absolute configuration determination was also revealed. Remarkably, the differentiation of two enantiomeric scaffolds (crinine and haemanthamine) was accomplished via OR/ORD calculations with cross-validation by ECD and VCD.

Quantum chemical calculation of the vapor pressure of volatile and semi volatile organic compounds

The vapor pressure is a specific and temperature-dependent parameter that describes the volatility of a substance and thus its driving force for evaporation or sublimation into the gas phase. Depending on the magnitude of the vapor pressure, there are different methods for experimental determination. However, these are usually associated with a corresponding amount of effort and become less accurate as the vapor pressure decreases. For purposes of vapor pressure prediction, algorithms were developed that are usually based on quantitative structure-activity relationships (QSAR). The quantum mechanical (QM) approach followed here applies an alternative, much less empirical strategy, where the change in Gibbs free energy for the transition from the condensed to the gas phase is obtained from conformer ensembles computed for each phase separately. The results of this automatic, so-called CRENSO workflow are compared with experimentally determined vapor pressures for a large set of environmentally relevant compounds. In addition, comparisons are made with the single structure-based COSMO-RS QM approach, linear-free-energy relationships (LFER) as well as results from the SPARC program. We show that our CRENSO workflow is superior to conventional prediction models and provides reliable vapor pressures for liquids and sub-cooled liquids over a wide pressure range.

A computational approach for modeling electronic circular dichroism of solvated chromophores

The present study consists in a novel computational protocol to model the UV-circular dichroism spectra of solvated species. It makes use of quantum-chemical calculations on a series of conformations of a flexible chromophore or on a series of chromophore/solvent clusters extracted from molecular dynamic simulations. The protocol is described and applied to the aqueous cationic tripeptide GAG⁺ and to the aqueous neutral decapeptide (GVGVP)₂ . The protocol has proven able to: (i) properly consider the conformational motion of solute in the given environment; (ii) give the actual statistical weight of each conformational state; (iii) provide a reliable quantum mechanical method able to reproduce the spectral features. Temperature effects on conformations and spectral properties are properly taken into account. The role of explicit solvent on the conformational analysis and the spectra calculation is discussed. The comparison of the calculated circular dichroism spectra with experimental ones recorded at different temperatures represents a strict validation test of the method.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

A Computational Protocol for Vibrational Circular Dichroism Spectra of Cyclic Oligopeptides

Cyclic peptides are a promising class of compounds for next-generation antibiotics as they may provide new ways of limiting antibiotic resistance development. Although their cyclic structure will introduce some rigidity, their conformational space is large and they usually have multiple chiral centers that give rise to a wide range of possible stereoisomers. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) are used to assign stereochemistry and discriminate enantiomers of chiral molecules, often in combination with electronic structure methods. The reliable determination of the absolute configuration of cyclic peptides will require robust computational methods than can identify all significant conformers and their relative population and reliably assign their stereochemistry from their chiroptical spectra by comparison with ab initio calculated spectra. We here present a computational protocol for the accurate calculation of the VCD spectra of a series of flexible cyclic oligopeptides. The protocol builds on the Conformer-Rotamer Ensemble Sampling Tool (CREST) developed by Grimme and co-workers ( Phys. Chem. Chem. Phys. 2020, 22, 7169-7192 and J. Chem. Theory. Comput. 2019, 15, 2847-2862) in combination with postoptimizations using B3LYP and moderately sized basis sets. Our recommended computational protocol for the computation of VCD spectra of cyclic oligopeptides consists of three steps: (1) conformational sampling with CREST and tight-binding density functional theory (xTB); (2) energy ranking based on single-point energy calculations as well as geometry optimization and VCD calculations of conformers that are within 2.5 kcal/mol of the most stable conformer using B3LYP/6-31+G*/CPCM; and (3) VCD spectra generation based on Boltzmann weighting with Gibbs free energies. Our protocol provides a feasible basis for generating VCD spectra also for larger cyclic peptides of biological/pharmaceutical interest and can thus be used to investigate promising compounds for next-generation antibiotics.

Computer-aided drug design, quantum-mechanical methods for biological problems

Quantum chemistry enables to study systems with chemical accuracy (<1 kcal/mol from experiment) but is restricted to a handful of atoms due to its computational expense. This has led to ongoing interest to optimize and simplify these methods while retaining accuracy. Implementing quantum mechanical (QM) methods on modern hardware such as multiple-GPUs is one example of how the field is optimizing performance. Multiscale approaches like the so-called QM/molecular mechanical method are gaining popularity in drug discovery because they focus the application of QM methods on the region of choice (e.g., the binding site), while using efficient MM models to represent less relevant areas. The creation of simplified QM methods is another example, including the use of machine learning to create ultra-fast and accurate QM models. Herein, we summarize recent advancements in the development of optimized QM methods that enhance our ability to use these methods in computer aided drug discovery.

Optimization of the r2SCAN-3c Composite Electronic-Structure Method for Use with Slater-Type Orbital Basis Sets

The "Swiss army knife" composite density functional electronic-structure method r²SCAN-3c (J. Chem. Phys. 2021, 154, 064103) is extended and optimized for the use with Slater-type orbital basis sets. The meta generalized-gradient approximation (meta-GGA) functional r²SCAN by Furness et al. is combined with a tailor-made polarized triple-ζ Slater-type atomic orbital (STO) basis set (mTZ2P), the semiclassical London dispersion correction (D4), and a geometrical counterpoise (gCP) correction. Relativistic effects are treated explicitly with the scalar-relativistic zeroth-order regular approximation (SR-ZORA). The performance of the new implementation is assessed on eight geometry and 74 energy benchmark sets, including the extensive GMTKN55 database as well as recent sets such as ROST61 and IONPI19. In geometry optimizations, the STO-based r²SCAN-3c is either on par with or more accurate than the hybrid density functional approximation M06-2X-D3(0)/TZP. In energy calculations, the overall accuracy is similar to the original implementation of r²SCAN-3c with Gaussian-type atomic orbitals (GTO), but basic properties, intermolecular noncovalent interactions, and barrier heights are better described with the STO approach, resulting in a lower weighted mean absolute deviation (WTMAD-2(STO) = 7.15 vs 7.50 kcal mol^-1 with the original method) for the GMTKN55 database. The STO-optimized r²SCAN-3c outperforms many conventional hybrid/QZ approaches in most common applications at a fraction of their cost. The reliable, robust, and accurate r²SCAN-3c implementation with STOs is a promising alternative to the original implementation with GTOs and can be generally used for a broad field of quantum chemical problems.

Hydrocarbon Macrocycle Conformer Ensembles and 13 C-NMR Spectra

NMR as a routine analytical method provides important three-dimensional structure information of compounds in solution. Here we apply the recently released CRENSO computational workflow for the automated generation of conformer ensembles to the quantum mechanical calculation of ¹³ C-NMR spectra of a series of flexible cycloalkanes up to C₂₀ H₄₀ . We evaluate the computed chemical shifts in comparison with corresponding experimental data in chloroform. It is shown that accurate and properly averaged theoretical NMR data can be obtained in about a day of computation time on a standard workstation computer. The excellent agreement between theory and experiment enables one to deduce preferred conformations of large, non-rigid macrocycles under ambient conditions from our automated procedure.

Quantum Chemical Calculation and Evaluation of Partition Coefficients for Classical and Emerging Environmentally Relevant Organic Compounds

Octanol/water (K_OW), octanol/air (K_OA), and hexadecane/air (K_HdA) partition coefficients are calculated for 67 organic compounds of environmental concern using computational chemistry. The extended CRENSO workflow applied here includes the calculation of extensive conformer ensembles with semiempirical methods and refinement through density functional theory, taking into account solvation models, especially COSMO-RS, and thermostatistical contributions. This approach is particularly advantageous for describing large and nonrigid molecules. With regard to K_OW and K_HdA, one can refer to many experimental data from direct and indirect measurement methods, and very good matches with results from our quantum chemical workflow are evident. In the case of the K_OA values, however, good matches are only obtained for the experimentally determined values. Larger systematic deviations between data computed here and available, nonexperimental quantitative structure-activity relationship literature data occur in particular for phthalic acid esters and organophosphate esters. From a critical analysis of the coefficients calculated in this work and comparison with available literature data, we conclude that the presented quantum chemical composite approach is the most powerful so far for calculating reliable partition coefficients because all physical contributions to the conformational free energy are considered and the structure ensembles for the two phases are generated independently and consistently.