Vollautomatisierte quantenchemische Berechnung von Spin-Spin- gekoppelten magnetischen Kernspinresonanzspektren

Predicting 195 Pt NMR Chemical Shifts in Water-Soluble Inorganic/Organometallic Complexes with a Fast and Simple Protocol Combining Semiempirical Modeling and Machine Learning

Water-soluble Pt complexes are the key components in medicinal chemistry and catalysis. The well-known cisplatin family of anticancer drugs and industrial hydrosylilation catalysts are two leading examples. On the molecular level, the activity mechanisms of such complexes mostly involve changes in the Pt coordination sphere. Using ¹⁹⁵ Pt NMR spectroscopy for operando monitoring would be a valuable tool for uncovering the activity mechanisms; however, reliable approaches for the rapid correlation of Pt complex structure with ¹⁹⁵ Pt chemical shifts are very challenging and not available for everyday research practice. While NMR shielding is a response property, molecular 3D structure determines NMR spectra, as widely known, which allows us to build up 3D structure to ¹⁹⁵ Pt chemical shift correlations. Accordingly, we present a new workflow for the determination of lowest-energy configurational/conformational isomers based on the GFN2-xTB semiempirical method and prediction of corresponding chemical shifts with a Machine Learning (ML) model tuned for Pt complexes. The workflow was designed for the prediction of ¹⁹⁵ Pt chemical shifts of water-soluble Pt(II) and Pt(IV) anionic, neutral, and cationic complexes with halide, NO₂ ^- , (di)amino, and (di)carboxylate ligands with chemical shift values ranging from -6293 to 7090 ppm. The model offered an accuracy (normalized root-mean-square deviation/RMSD) of 1.08 %/145.02 ppm on the held-out test set.

Unusual Enantiodivergence in Chiral Brønsted Acid‐Catalyzed Asymmetric Allylation with β‐Alkenyl Allylic Boronates

We report herein a rare example of enantiodivergent aldehyde addition with β-alkenyl allylic boronates via chiral Brønsted acid catalysis. 2,6-Di-9-anthracenyl-substituted chiral phosphoric acid-catalyzed asymmetric allylation using β-vinyl substituted allylic boronate gave alcohols with R absolute configuration. The sense of asymmetric induction of the catalyst in these reactions is opposite to those in prior reports. Moreover, in the presence of the same acid catalyst, the reactions with β-2-propenyl substituted allylic boronate generated homoallylic alcohol products with S absolute configuration. Unusual substrate-catalyst C-H⋅⋅⋅π interactions in the favoured reaction transition state were identified as the origins of observed enantiodivergence through DFT computational studies.

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.

Titanocene-Catalyzed [2+2] Cycloaddition of Bisenones and Comparison with Photoredox Catalysis and Established Methods

Cp2 Ti(TFA) is a broadly applicable catalyst for the [2+2] cycloaddition of bisenones by inner-sphere electron transfer. The attractiveness of this mechanism is shown by comparison with outer-sphere ET methods. DFT calculations show that the reaction proceeds through a unique unfavorable 5-exo (the rate-determining step) and a favorable 4-exo cyclization.

Matrix Isolation-Vibrational Circular Dichroism Spectroscopic Study of Conformations and Non-Covalent Interactions of Tetrahydro-2-Furoic Acid

The conformational landscape and aggregation behaviour of tetrahydro-2-furoic acid (THFA) were investigated by using matrix isolation-vibrational circular dichroism (MI-VCD). The well-resolved experimental MI-IR and MI-VCD features in an argon matrix at 10 K allow one to identify two dominant monomeric conformations as trans-THFA where the hydroxyl and carbonyl groups of COOH are at opposite sides, as well as one cis-conformer. At 24 K and 30 K deposition temperatures, the experimental IR and VCD spectral features reveal further growth of the binary THFA aggregates. Systematic conformational searches identified three vastly different binary binding topologies, resulting in a few hundred stable (THFA)₂ conformers. Interestingly, the main binary structures observed correspond to an unusual type of structure which is made of two trans-THFA subunits, in contrast to the usual double H-bonded ring binary structures, identified in a previous solution study. The present work showcases the power of MI-VCD spectroscopy in revealing unusual structures formed in a cold rare gas matrix.

Surprising Homolytic Gas Phase Co-C Bond Dissociation Energies of Organometallic Aryl-Cobinamides Reveal Notable Non-Bonded Intramolecular Interactions

Aryl-cobalamins are a new class of organometallic structural mimics of vitamin B12 designed as potential 'antivitamins B12 '. Here, the first cationic aryl-cobinamides are described, which were synthesized using the newly developed diaryl-iodonium method. The aryl-cobinamides were obtained as pairs of organometallic coordination isomers, the stereo-structure of which was unambiguously assigned based on homo- and heteronuclear NMR spectra. The availability of isomers with axial attachment of the aryl group, either at the 'beta' or at the 'alpha' face of the cobalt-center allowed for an unprecedented comparison of the organometallic reactivity of such pairs. The homolytic gas-phase bond dissociation energies (BDEs) of the coordination-isomeric phenyl- and 4-ethylphenyl-cobinamides were determined by ESI-MS threshold CID experiments, furnishing (Co-Csp 2 )-BDEs of 38.4 and 40.6 kcal mol-1 , respectively, for the two β-isomers, and the larger BDEs of 46.6 and 43.8 kcal mol-1 for the corresponding α-isomers. Surprisingly, the observed (Co-Csp 2 )-BDEs of the Coβ -aryl-cobinamides were smaller than the (Co-Csp 3 )-BDE of Coβ -methyl-cobinamide. DFT studies and the magnitudes of the experimental (Co-Csp 2 )-BDEs revealed relevant contributions of non-bonded interactions in aryl-cobinamides, notably steric strain between the aryl and the cobalt-corrin moieties and non-bonded interactions with and among the peripheral sidechains.

Design Platform for Sustainable Catalysis with Radicals: Electrochemical Activation of Cp2 TiCl2 for Catalysis Unveiled

The combination of synthesis, rotating ring-disk electrode (RRDE) and cyclic voltammetry (CV) measurements, and computational investigations with the aid of DFT methods shows how a thiourea, a squaramide, and a bissulfonamide as additives affect the Eq Cr equilibrium of Cp2 TiCl2 . We have, for the first time, provided quantitative data for the Eq Cr equilibrium and have determined the stoichiometry of adduct formation of [Cp2 Ti(III)Cl2 ]- , [Cp2 Ti(III)Cl] and [Cp2 Ti(IV)Cl2 ] and the additives. By studying the structures of the complexes formed by DFT methods, we have established the Gibbs energies and enthalpies of complex formation as well as the adduct structures. The results not only demonstrate the correctness of our use of the Eq Cr equilibrium as predictor for sustainable catalysis. They are also a design platform for the development of novel additives in particular for enantioselective catalysis.

Oxidation Under Reductive Conditions: From Benzylic Ethers to Acetals with Perfect Atom-Economy by Titanocene(III) Catalysis

Described here is a titanocene-catalyzed reaction for the synthesis of acetals and hemiaminals from benzylic ethers and benzylic amines, respectively, with pendant epoxides. The reaction proceeds by catalysis in single-electron steps. The oxidative addition comprises an epoxide opening. An H-atom transfer, to generate a benzylic radical, serves as a radical translocation step, and an organometallic oxygen rebound as a reductive elimination. The reaction mechanism was studied by high-level dispersion corrected hybrid functional DFT with implicit solvation. The low-energy conformational space was searched by the efficient CREST program. The stereoselectivity was deduced from the lowest lying benzylic radical structures and their conformations are controlled by hyperconjugative interactions and steric interactions between the titanocene catalyst and the aryl groups of the substrate. An interesting mechanistic aspect is that the oxidation of the benzylic center occurs under reducing conditions.

Unlocking the Water Trimer Loop: Isotopic Study of Benzophenone-(H2 O)1-3 Clusters with Rotational Spectroscopy

Examined here are the structures of complexes of benzophenone microsolvated with up to three water molecules by using broadband rotational spectroscopy and the cold conditions of a molecular jet. The analysis shows that the water molecules dock sideways on benzophenone for the water monomer and dimer moieties, and they move above one of the aromatic rings when the water cluster grows to the trimer. The rotational spectra shows that the water trimer moiety in the complex adopts an open‐loop arrangement. Ab initio calculations face a dilemma of identifying the global minimum between the open loop and the closed loop, which is only solved when zero‐point vibrational energy correction is applied. An OH⋅⋅⋅π bond and a Bürgi‐Dunitz interaction between benzophenone and the water trimer are present in the cluster. This work shows the subtle balance between water–water and water–solute interactions when the solute molecule offers several different anchor sites for water molecules.

Quantification of Noncovalent Interactions in Azide-Pnictogen, -Chalcogen, and -Halogen Contacts

The noncovalent interactions between azides and oxygen‐containing moieties are investigated through a computational study based on experimental findings. The targeted synthesis of organic compounds with close intramolecular azide–oxygen contacts yielded six new representatives, for which X‐ray structures were determined. Two of those compounds were investigated with respect to their potential conformations in the gas phase and a possible significantly shorter azide–oxygen contact. Furthermore, a set of 44 high‐quality, gas‐phase computational model systems with intermolecular azide–pnictogen (N, P, As, Sb), –chalcogen (O, S, Se, Te), and –halogen (F, Cl, Br, I) contacts are compiled and investigated through semiempirical quantum mechanical methods, density functional approximations, and wave function theory. A local energy decomposition (LED) analysis is applied to study the nature of the noncovalent interaction. The special role of electrostatic and London dispersion interactions is discussed in detail. London dispersion is identified as a dominant factor of the azide–donor interaction with mean London dispersion energy‐interaction energy ratios of 1.3. Electrostatic contributions enhance the azide–donor coordination motif. The association energies range from −1.00 to −5.5 kcal mol⁻¹.

Discovering the Elusive Global Minimum in a Ternary Chiral Cluster: Rotational Spectra of Propylene Oxide Trimer

The chirality controlled conformational landscape of the trimer of propylene oxide (PO), a prototypical chiral molecule, was investigated using rotational spectroscopy and a range of theoretical tools for conformational searches and for evaluating vibrational contributions to effective structures. Two sets of homochiral (PO)3 rotational transitions were assigned and the associated conformers identified with theoretical support. One set of heterochiral (PO)3 transitions was assigned, but no structures generated by one of the latest, advanced conformational search codes could account for them. With the aid of a Python program, the carbon atom backbone and then the heterochiral (PO)3 structure were generated using 13 C isotopic data. Excellent agreement between theoretical and experimental rotational constants and relative dipole moment components of all three conformers was achieved, especially after applying vibrational corrections to the rotational constants.

Predicting 19 F NMR Chemical Shifts: A Combined Computational and Experimental Study of a Trypanosomal Oxidoreductase-Inhibitor Complex

The absence of fluorine from most biomolecules renders it an excellent probe for NMR spectroscopy to monitor inhibitor-protein interactions. However, predicting the binding mode of a fluorinated ligand from a chemical shift (or vice versa) has been challenging due to the high electron density of the fluorine atom. Nonetheless, reliable 19 F chemical-shift predictions to deduce ligand-binding modes hold great potential for in silico drug design. Herein, we present a systematic QM/MM study to predict the 19 F NMR chemical shifts of a covalently bound fluorinated inhibitor to the essential oxidoreductase tryparedoxin (Tpx) from African trypanosomes, the causative agent of African sleeping sickness. We include many protein-inhibitor conformations as well as monomeric and dimeric inhibitor-protein complexes, thus rendering it the largest computational study on chemical shifts of 19 F nuclei in a biological context to date. Our predicted shifts agree well with those obtained experimentally and pave the way for future work in this area.

Structure Optimisation of Large Transition-Metal Complexes with Extended Tight-Binding Methods

Large transition-metal complexes are used in numerous areas of chemistry. Computer-aided theoretical investigations of such complexes are limited by the sheer size of real systems often consisting of hundreds to thousands of atoms. Accordingly, the development and thorough evaluation of fast semi-empirical quantum chemistry methods that are universally applicable to a large part of the periodic table is indispensable. Herein, we report on the capability of the recently developed GFNn-xTB method family for full quantum-mechanical geometry optimisation of medium to very large transition-metal complexes and organometallic supramolecular structures. The results for a specially compiled benchmark set of 145 diverse closed-shell transition-metal complex structures for all metals up to Hg are presented. Further the GFNn-xTB methods are tested on three established benchmark sets regarding reaction energies and barrier heights of organometallic reactions.

Computational Chemistry: The Fate of Current Methods and Future Challenges

"Where do we go from here?" is the underlying question regarding the future (perhaps foreseeable) developments in computational chemistry. Although this young discipline has already permeated practically all of chemistry, it is likely to become even more powerful with the rapid development of computational hard- and software.