Determining the Factors accounting for Reaction Selectivity A Relative Energy Gradient - Interacting Quantum Atoms and Natural Bonding Orbitals Study.

Identifying the main physicochemical properties accounting for the course of a reaction is of utmost importance to rationalize chemical syntheses. To this aim, the relative energy gradient (REG) method is an appealing approach because it is an unbiased and automatic process to extract the most relevant pieces of energy information. Initially formulated within the interacting quantum atoms (IQA) framework for a single reaction, here we extend the REG method to natural bond orbitals (NBO) analysis and to the case of two competitive processes. This development enables the determination of the driving forces of any chemical selectivity. We illustrate the extended REG method on the case study of ring opening in cyclobutenes, which is an important instance of the so-called torquoselectivity.

Expanding horizons in conceptual density functional theory: Novel ensembles and descriptors to decipher reactivity patterns.

Conceptual density functional theory (CDFT) and the quantum reactivity descriptors stemming from it have proven to be valuable tools for understanding the chemical behavior of molecules. This article is presented as being intrinsically of dual character. In a first part, it briefly reviews, in a deliberately didactical way, the main ensembles in CDFT, while the second half presents two additional ensembles, where the chemical hardness acts as a natural variable, and their respective reactivity descriptors. The evaluation of these reactivity descriptors on common organic chemical reagents are presented and discussed.

Aza-Michael Addition in Explicit Solvent: A Relative Energy Gradient-Interacting Quantum Atoms Study.

Aza-Michael additions are key reactions in organic synthesis. We investigate, from a theoretical and computational point of view, several examples ranging from weak to strong electrophiles in dimethylsulfoxide treated as explicit solvent. We use the REG-IQA method, which is a quantum topological energy decomposition (Interacting Quantum Atoms, IQA) coupled to a chemical-interpretation calculator (Relative Energy Gradient, REG). We focus on the rate-limiting addition step in order to unravel the different events taking place in this step, and understand the influence of solvent on the reaction, with an eye on predicting the Mayr electrophilicity. For the first time, a link is established between an REG-IQA analysis and experimental values.

Metalloenzyme-Mediated Thiol-Yne Addition Towards Photoisomerizable Fluorescent Dyes.

Proteins are able to irreversibly assemble biologically active ligands from building blocks bearing complementary reactive functions due their spatial proximity, through a kinetic target-guided synthetic process (also named in situ click chemistry). Although linkages thus formed are mostly passive, some of them have shown to significantly contribute to the protein binding through for instance hydrogen bonding and stacking interactions. Biocompatible reactions and click chemistry are a formidable source of inspiration for developing such new protein-directed ligations. This study reports a proximity-induced thiol-yne synthesis of carbonic anhydrase inhibitors. Not only this example widens the arsenal of kinetic target-guided synthesis (KTGS) eligible reactions, but the obtained product displayed unsuspected photophysical properties. The corresponding vinyl sulfide linkage conjugated to a coumarin core proved to be engaged in a monodirectional Z to E photoisomerization process. Further investigations guided by theoretical calculations showed that fine-tuning of the nature of the substituents on the coumarin moiety allows to obtain a bidirectional photochemical process, thus discovering a new photoswitching moiety, displaying moreover fluorescence properties. Due to the spectral tunability of coumarin derivatives, this work should open new opportunities for the design of vinyl sulfide-based photoswitch systems with modular photophysical properties.

Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations.

We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.

Polarisation of Electron Density and Electronic Effects: Revisiting the Carbon-Halogen Bonds.

Electronic effects (inductive and mesomeric) are of fundamental importance to understand the reactivity and selectivity of a molecule. In this article, polarisation temperature is used as a principal index to describe how electronic effects propagate in halogeno-alkanes and halogeno-alkenes. It is found that as chain length increases, polarisation temperature decreases. As expected, polarisation is much larger for alkenes than for alkanes. Finally, the polarisation mode of the carbon-fluorine bond is found to be quite different and might explain the unusual reactivity of fluoride compounds.

Electrophilicity Indices and Halogen Bonds: Some New Alternatives to the Molecular Electrostatic Potential.

We assess the ability of various atomic and molecular electrophilicity descriptors to predict the strength of halogen bonds. To this aim, several physicochemical quantities rooted within the framework of conceptual density functional theory were derived using second and third order Taylor expansions of the electronic energy, and their correlation to binding energies were compared with those obtained for more usual electronic descriptors. This benchmark was performed for a large and representative database of noncovalent complexes involving fluorine, chlorine, and bromine atoms, and showed that some of these new quantities, in particular the atomic cubic electrophilicity index, exhibited more transferability and broadness of application than did more common descriptors such as the molecular electrostatic potential.

Characterization of ß-turns by electronic circular dichroism spectroscopy: a coupled molecular dynamics and time-dependent density functional theory computational study.

Electronic circular dichroism is one of the most used spectroscopic techniques for peptide and protein structural characterization. However, while valuable experimental spectra exist for α-helix, ß-sheet and random coil secondary structures, previous studies showed important discrepancies for ß-turns, limiting their use as a reference for structural studies. In this paper, we simulated circular dichroism spectra for the best-characterized ß-turns in peptides, namely types I, II, I' and II'. In particular, by combining classical molecular dynamics simulations and state-of-the-art quantum time-dependent density functional theory (with the polarizable embedding multiscale model) computations, two common electronic circular dichroism patterns were found for couples of ß-turn types (namely, type I/type II' and type II/type I'), at first for a minimal di-peptide model (Ace-Ala-Ala-NHMe), but also for all sequences tested with non-aromatic residues in the central positions. On the other hand, as expected, aromatic substitution causes important perturbations to the previously found ECD patterns. Finally, by applying suitable approximations, these patterns were subsequently rationalized based on the exciton chirality rule. All these results provide useful predictions and pave the way for a possible experimental characterization of ß-turns based on circular dichroism spectroscopy.

Rational design of novel fluorescent analogues of cholesterol: a "step-by-step" computational study.

In this paper we show a theoretical rational design approach on a series of intrinsically fluorescent analogues of cholesterol (FLACs), called polyene-sterols (P-sterols), followed by a step-by-step selection of potential candidates, employing, sequentially, state-of-the-art quantum mechanical (QM) computations of the optical properties (single- and multiphoton absorption electronic spectroscopies and emission), molecular dynamics (MD) simulations in model membranes, and multiscale approaches (polarizable embedding). This selection converged to a promising candidate that shows simultaneously interesting single- and multiphoton absorption properties as well as emitting properties and good abilities to mimic cholesterol order effects in model membranes.

Nine questions on energy decomposition analysis.

The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.

Can molecular and atomic descriptors predict the electrophilicity of Michael acceptors?

In this paper, we assess the ability of various intrinsic molecular and atomic descriptors, grounded in the conceptual density functional theory and the quantum theory of atoms-in-molecules frameworks, to predict the electrophilicity of Michael acceptors, which are fundamental organic reagents involved in the formation of carbon-carbon bonds. To this aim, linear and multilinear regressions between these theoretical properties and the experimental values gathered in Mayr-Patz' scale were performed. The relevance of quantum chemical descriptors are then discussed.

Decomposition of Møller-Plesset Energies within the Quantum Theory of Atoms-in-Molecules.

We discuss two main approaches to decompose the Møller-Plesset perturbation theory molecular energies into atomic contributions within the interacting quantum atoms (IQA) formalism, as implemented in the programs Morphy and AIMAll. For this purpose, the so-called intraatomic energies (also known as self-energies) of a representative set of 55 small molecules are compared with each other. The origin of the possible discrepancies between both approaches is analyzed, and linear models linking the two approaches are proposed for each atom type.

Revisiting absorption and electronic circular dichroism spectra of cholesterol in solution: a joint experimental and theoretical study.

Cholesterol is doubtless one of the most studied bio-molecules, which unfortunately features low emitting properties, precluding its in vivo study by fluorescence experiments. The design of fluorescent analogues of cholesterol is thus an appealing challenge in biochemistry, which simultaneously requires minor changes in its chemical structure (to retain main biological properties) and considerable enhancement of light emission. To this aim, the photochemical behaviour of the native molecule has to be deeply understood. In this work, we focused our attention on the electronic absorption of cholesterol in several common organic solutions, combining experimental (through ultraviolet-visible and electronic circular dichroism spectroscopy) and theoretical approaches (at the time-dependent density functional theory level) in order to solve the important discrepancies previously reported in the literature on the maximum absorption wavelengths and on the nature (Rydberg and/or π → π*) of the associated electronic transition.

A Unique (3+2) Annulation Reaction between Meldrum's Acid and Nitrones: Mechanistic Insight by ESI-IMS-MS and DFT Studies.

The fragile intermediates of the domino process leading to an isoxazolidin-5-one, triggered by unique reactivity between Meldrum's acid and an N-benzyl nitrone in the presence of a Brønsted base, were determined thanks to the softness and accuracy of electrospray ionization mass spectrometry coupled to ion mobility spectrometry (ESI-IMS-MS). The combined DFT study shed light on the overall organocatalytic sequence that starts with a stepwise (3+2) annulation reaction that is followed by a decarboxylative protonation sequence encompassing a stereoselective pathway issue.

Anthracenyl polar embedded stationary phases with enhanced aromatic selectivity. Part II: A density functional theory study.

New polar embedded aromatic stationary phases (mono- and trifunctional versions) that contain an amide-embedded group coupled with a tricyclic aromatic moiety were developed for chromatographic applications and described in the first paper of this series. These phases offered better separation performance for PAHs than for alkylbenzene homologues, and an enhanced ability to differentiate aromatic planarity to aromatic tridimensional conformation, especially for the trifunctional version and when using methanol instead of acetonitrile. In this second paper, a density functional theory study of the retention process is reported. In particular, it was shown that the selection of the suitable computational protocol allowed for describing rigorously the interactions that could take place, the solvent effects, and the structural changes for the monofunctional and the trifunctional versions. For the first time, the experimental data coupled with these DFT results provided a better understanding of the interaction mechanisms and highlighted the importance of the multimodal character of the designed stationary phases: alkyl spacers for interactions with hydrophobic solutes, amide embedded groups for dipole-dipole and hydrogen-bond interactions, and aromatic terminal groups for π-π interactions.

On Atoms-in-Molecules Energies from Kohn-Sham Calculations.

Herein, we discuss three methods to partition the total molecular energy into additive atomic contributions within the framework of Bader's atoms-in-molecules theory and in the particular context of Kohn-Sham density functional theory. The first method is derived from the virial theorem, whereas the two other schemes, termed "standard" and "model", are based on Pendás' interacting-quantum-atoms decomposition. The methods are then compared for a dataset of molecules of interest for direct application in organic chemistry and biochemistry. Finally, the relevance of the three methods for the prediction of intrinsic reactivity properties (e.g., electrophilicity) or for unravelling the nature of chemical bonding (e.g., in halogen bonds, beyond the pure electrostatic point of view), is examined and paves the way for their more systematic use for the in silico design of new reactants.

Palladium-Catalyzed Synthesis of 3-Trifluoromethyl-Substituted 1,3-Butadienes by Means of Directed C-H Bond Functionalization.

A palladium-catalyzed C-H bond functionalization of acrylamides was developed to build up stereoselectively trifluoromethylated 1,3-butadienes. Using a tertiary amide as a directing group, olefins were selectively functionalized with 2-bromo-3,3,3-trifluoropropene to access these important fluorinated compounds. The methodology was extended to the construction of pentafluoroethyl-substituted 1,3-dienes. Mechanistic studies supported by density functional theory calculations suggested a redox neutral mechanism for this transformation.

On the importance of intramolecular hydrogen bond cooperativity in d-glucose - an NMR and QTAIM approach.

The idea that hydrogen bond cooperativity is responsible for the structure and reactivity of carbohydrates is examined. Density functional theory and gauge-including atomic orbital calculations on the known conformers of the α and ß anomers of d-glucopyranose in the gas phase are used to compute proton NMR chemical shifts and interatomic distances, which are taken as criteria for probing intramolecular interactions. Atom-atom interaction energies are calculated by the interacting quantum atoms approach in the framework of the quantum theory of atoms in molecules. Association of OH1 in the counterclockwise conformers with a strong acceptor, pyridine, is accompanied by cooperative participation from OH2, but there is no significant change in the bonding of the two following 1,2-diol motifs. The OH6... O5 (G-g+/cc/t and G+g-/cc/t conformers) or OH6... O4 (Tg+/cc/t conformer) distance is reduced, and the OH6 proton is slightly deshielded. In the latter case, this shortening and the associated increase in the OH6-O4 interaction energy may be interpreted as a small cooperative effect, but intermolecular interaction energies are practically the same for all three conformers. In most of the pyridine complexes, one ortho proton interacts with the endocyclic oxygen O5. Analogous results are obtained when the clockwise conformer, G-g+/cl/g-, detected for the α anomer, and a hypothetical conformer, Tt/cl/g-, are complexed with pyridine through OH6. Generally, the cooperative effect does not go beyond the first two OH groups of a chain. Copyright © 2017 John Wiley & Sons, Ltd.

A theoretical study of the diastereoselective allylation of aldehydes with new chiral allylsilanes.

In this paper, we investigate the reaction mechanism for the synthesis of new functionalized chiral hydroxyamides theoretically, using a concerted density functional theory (DFT)-conceptual DFT-quantum theory of atoms-in-molecules (QTAIM) approach, which casts light on the main physicochemical properties responsible for the observed selectivity. We use a particular nucleophilic addition step to illustrate the strengths of this general strategy, which could, in principle, be applied to unravel the main features of any chemical synthesis.

Molecular Relaxations in Supercooled Liquid and Glassy States of Amorphous Quinidine: Dielectric Spectroscopy and Density Functional Theory Approaches.

In this article, we conduct a comprehensive molecular relaxation study of amorphous Quinidine above and below the glass-transition temperature (Tg) through broadband dielectric relaxation spectroscopy (BDS) experiments and theoretical density functional theory (DFT) calculations, as one major issue with the amorphous state of pharmaceuticals is life expectancy. These techniques enabled us to determine what kind of molecular motions are responsible, or not, for the devitrification of Quinidine. Parameters describing the complex molecular dynamics of amorphous Quinidine, such as Tg, the width of the α relaxation (ßKWW), the temperature dependence of α-relaxation times (τα), the fragility index (m), and the apparent activation energy of secondary γ relaxation (Ea-γ), were characterized. Above Tg (> 60 °C), a medium degree of nonexponentiality (ßKWW = 0.5) was evidenced. An intermediate value of the fragility index (m = 86) enabled us to consider Quinidine as a glass former of medium fragility. Below Tg (< 60 °C), one well-defined secondary γ relaxation, with an apparent activation energy of Ea-γ = 53.8 kJ/mol, was reported. From theoretical DFT calculations, we identified the most reactive part of Quinidine moieties through exploration of the potential energy surface. We evidenced that the clearly visible γ process has an intramolecular origin coming from the rotation of the CH(OH)C9H14N end group. An excess wing observed in amorphous Quinidine was found to be an unresolved Johari-Goldstein relaxation. These studies were supplemented by sub-Tg experimental evaluations of the life expectancy of amorphous Quinidine by X-ray powder diffraction and differential scanning calorimetry. We show that the difference between Tg and the onset temperature for crystallization, Tc, which is 30 K, is sufficiently large to avoid recrystallization of amorphous Quinidine during 16 months of storage under ambient conditions.