N-(Sulfonio)Sulfilimine Reagents: Non-Oxidizing Sources of Electrophilic Nitrogen Atom for Skeletal Editing.

The one-pot synthesis of λ4-dibenzothiophen-5-imino-N-dibenzothiophenium triflate (1) in multigram scale is reported. This compound reacts with Rh2(esp)2 (esp=α,α,α',α'-tetramethyl-1,3-benzenedipropionic acid) generating a Rh-coordinated sulfonitrene species, which is able to transfer the electrophilic nitrene moiety to olefins. When indenes are used as substrates, isoquinolines are obtained in good yields. We assumed that after formation of the corresponding N-sulfonio aziridine, a ring expansion occurs via selective C-C bond cleavage and concomitant elimination of dibenzothiophene. Unexpectedly, a similar protocol transforms 1-arylcyclobutenes into 1-cyano-1-arylcyclopropanes. Our calculations indicate that aziridination is not favored in this case; instead, sulfilimine-substituted cyclobutyl carbocations are initially formed, and these evolve to the isolated cyclopropanes via ring contraction. Both procedures are operationally simple, tolerate a range of functional groups, including oxidation-sensitive alcohols and aldehydes, and enable the convenient preparation of valuable 15N-labelled products. These results demonstrate the potential of 1 to provide alternative pathways for the selective transfer of N-atoms in organic molecules.

Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach.

The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.

Dispersion Control over Molecule Cohesion: Exploiting and Dissecting the Tipping Power of Aromatic Rings.

ConspectusWe have learned over the past years how London dispersion forces can be effectively used to influence or even qualitatively tip the structure of aggregates and the conformation of single molecules. This happens despite the fact that single dispersion contacts are much weaker than competing polar forces. It is a classical case of strength by numbers, with the importance of London dispersion forces scaling with the system size. Knowledge about the tipping points, however difficult to attain, is necessary for a rational design of intermolecular forces. One requires a careful assessment of the competing interactions, either by sensitive spectroscopic techniques for the study of the isolated molecules and aggregates or by theoretical approaches. Of particular interest are the systems close to the tipping point, when dispersion interactions barely outweigh or approach the strength of the other interactions. Such subtle cases are important milestones for a scale-up to realistic multi-interaction situations encountered in the fields of life and materials science. In searching for examples that provide ideal competing interactions in complexes and small clusters, aromatic systems can offer a diverse set of molecules with a variation of dispersion and electrostatic forces that control the dominant and peripheral interactions. Our combined spectroscopic and theoretical investigations provide valuable insights into the balance of intermolecular forces because they typically allow us to switch the aromatic substituent on and off. High-resolution rotational spectroscopy serves as a benchmark for molecular structures, as correct calculations should be based on correct geometries. When discussing the competition with other noncovalent interactions, obvious competitors are directional hydrogen bonds. As a second counterweight to aryl interactions, we will discuss aurophilic/metallophilic interactions, which also have a strong stabilization with a small number of atoms involved. Vibrational spectroscopy is most sensitive to interactions of light atoms, and the competition of OH hydrogen bonds with dispersion forces in a molecular aggregate can be judged well by the OH stretching frequency. Experiments in the gas phase are ideal for gauging the accuracy of quantum chemical predictions free of solvent forces. A tight collaboration utilizing these three methods allows experiment vs experiment vs theory benchmarking of the overall influence of dispersion in molecular structures and energetics.

Multiple redox switches of the SARS-CoV-2 main protease in vitro provide opportunities for drug design.

Besides vaccines, the development of antiviral drugs targeting SARS-CoV-2 is critical for preventing future COVID outbreaks. The SARS-CoV-2 main protease (Mpro), a cysteine protease with essential functions in viral replication, has been validated as an effective drug target. Here, we show that Mpro is subject to redox regulation in vitro and reversibly switches between the enzymatically active dimer and the functionally dormant monomer through redox modifications of cysteine residues. These include a disulfide-dithiol switch between the catalytic cysteine C145 and cysteine C117, and generation of an allosteric cysteine-lysine-cysteine SONOS bridge that is required for structural stability under oxidative stress conditions, such as those exerted by the innate immune system. We identify homo- and heterobifunctional reagents that mimic the redox switching and inhibit Mpro activity. The discovered redox switches are conserved in main proteases from other coronaviruses, e.g. MERS-CoV and SARS-CoV, indicating their potential as common druggable sites.

Nuclear Quantum Effects Made Accessible: Local Density Fitting in Multicomponent Methods.

The simulation of nuclear quantum effects (NQEs) is crucial for an accurate description of systems and processes involving light nuclei, such as hydrogen atoms. Within the last years, the importance of those effects has been highlighted for a vast range of systems with tremendous implications in chemistry, biology, physics, and materials sciences. However, while electronic structure theory methods have become routine tools for quantum chemical investigations, there is still a lack of approaches to address NQEs that are computationally accessible and straightforward to use. To address this, we present the first combination of the nuclear-electronic orbital Hartree-Fock approach with both local and density fitting approximations (LDF-NEO-HF). This results in a low-order scaling approach that enables the inclusion of NQEs for large systems within a fraction of a day and for small to medium size systems in minutes. Moreover, we demonstrate the qualitative accuracy and robustness of our approach to retrieve NQEs for three real-use cases motivated by chemical, biological, and materials science applications.

Dynamic Protonation States Underlie Carbene Formation in ThDP-Dependent Enzymes: A Theoretical Study.

The activation mechanism of thiamine diphosphate (ThDP) in enzymes has long been the subject of intense research and controversial discussion. Particularly contentious is the formation of a carbene intermediate, the first one observed in an enzyme. For the formation of the carbene to take place, both intramolecular and intermolecular proton transfer pathways have been proposed. However, the physiologically relevant pH of ThDP-dependent enzymes around neutrality does not seem to be suitable for the formation of such reactive chemical species. Herein, we investigate the general mechanism of activation of the ThDP cofactor in human transketolase (TKT), by means of electronic structure methods. We show that in the case of the human TKT, the carbene species is accessible through a pKa shift induced by the electrostatics of a neighboring histidine residue (H110), whose protonation state change modulates the pKa of ThDP and suppresses the latter by more than 6 pH units. Our findings highlight that ThDP enzymes activate the cofactor beyond simple geometric constraints and the canonical glutamate. Such observations in nature can pave the way for the design of biomimetic carbene catalysts and the engineering of tailored enzymatic carbenes.

The first HyDRA challenge for computational vibrational spectroscopy.

Vibrational spectroscopy in supersonic jet expansions is a powerful tool to assess molecular aggregates in close to ideal conditions for the benchmarking of quantum chemical approaches. The low temperatures achieved as well as the absence of environment effects allow for a direct comparison between computed and experimental spectra. This provides potential benchmarking data which can be revisited to hone different computational techniques, and it allows for the critical analysis of procedures under the setting of a blind challenge. In the latter case, the final result is unknown to modellers, providing an unbiased testing opportunity for quantum chemical models. In this work, we present the spectroscopic and computational results for the first HyDRA blind challenge. The latter deals with the prediction of water donor stretching vibrations in monohydrates of organic molecules. This edition features a test set of 10 systems. Experimental water donor OH vibrational wavenumbers for the vacuum-isolated monohydrates of formaldehyde, tetrahydrofuran, pyridine, tetrahydrothiophene, trifluoroethanol, methyl lactate, dimethylimidazolidinone, cyclooctanone, trifluoroacetophenone and 1-phenylcyclohexane-cis-1,2-diol are provided. The results of the challenge show promising predictive properties in both purely quantum mechanical approaches as well as regression and other machine learning strategies.

Pentafluorocyclopropanation of (Hetero)arenes Using Sulfonium Salts: Applications in Late-Stage Functionalization.

The evaluation of the pentafluorocyclopropyl group as a chemotype in crop protection and medicinal chemistry has been hampered in the past by the lack of suitable methodologies that enable the practical incorporation of this moiety into advanced synthetic intermediates. Herein, we report the gram-scale synthesis of an unprecedented sulfonium salt, 5-(pentafluorocyclopropyl)dibenzothiophenium triflate, and its use as a versatile reagent for the photoinduced C-H pentafluorocyclopropylation of a broad series of non-previously functionalized (hetero)arenes through a radical mediated mechanism. The scope and potential benefits of the protocol developed are further demonstrated by the late-stage introduction of the pentafluorocyclopropyl unit into biologically relevant molecules and widely used pharmaceuticals.

Revisiting the origin of the bending in group 2 metallocenes AeCp2 (Ae = Be-Ba).

Metallocenes are well-established compounds in organometallic chemistry, and can exhibit either a coplanar structure or a bent structure according to the nature of the metal center (E) and the cyclopentadienyl ligands (Cp). Herein, we re-examine the chemical bonding to underline the origins of the geometry and stability observed experimentally. To this end, we have analysed a series of group 2 metallocenes [Ae(C5R5)2] (Ae = Be-Ba and R = H, Me, F, Cl, Br, and I) with a combination of computational methods, namely energy decomposition analysis (EDA), polarizability model (PM), and dispersion interaction densities (DIDs). Although the metal-ligand bonding nature is mainly an electrostatic interaction (65-78%), the covalent character is not negligible (33-22%). Notably, the heavier the metal center, the stronger the d-orbital interaction with a 50% contribution to the total covalent interaction. The dispersion interaction between the Cp ligands counts only for 1% of the interaction. Despite that orbital contributions become stronger for heavier metals, they never represent the energy main term. Instead, given the electrostatic nature of the metallocene bonds, we propose a model based on polarizability, which faithfully predicts the bending angle. Although dispersion interactions have a fair contribution to strengthen the bending angle, the polarizability plays a major role.

Mechanisms of Cysteine-Lysine Covalent Linkage-The Role of Reactive Oxygen Species and Competition with Disulfide Bonds.

Recently, a new naturally occurring covalent linkage was characterised, involving a cysteine and a lysine, bridged through an oxygen atom. The latter was dubbed as the NOS bond, reflecting the individual atoms involved in this uncommon bond which finds little parallel in lab chemistry. It is found to form under oxidising conditions and is reversible upon addition of reducing agents. Further studies have identified the bond in crystal structures across a variety of systems and organisms, potentially playing an important role in regulation, cellular defense and replication. Not only that, double NOS bonds have been identified and even found to be competitive in relation to the formation of disulfide bonds. This raises several questions about how this exotic bond comes to be, what are the intermediates involved in its formation and how it competes with other pathways of sulfide oxidation. With this objective in mind, we revisited our first proposed mechanism for the reaction with model electronic structure calculations, adding information about the reactivity with alternative reactive oxygen species and other potential competing products of oxidation. We present a network with more than 30 reactions which provides one of the most encompassing pictures for cysteine oxidation pathways to date.

The Catalytic Mechanism of Acetoacetate Decarboxylase: A Detailed Study of Schiff Base Formation, Protonation States, and Their Impact on Catalysis.

The enzyme acetoacetate decarboxylase (AAD) has a crucial function in the process of decarboxylating the substrate acetoacetate (AA). It has been extensively studied over the years, but its exact catalytic mechanism has remained partly unsolved due to the difficulty in assessing reaction intermediates. In this study, we combine molecular dynamics and electronic structure calculations to rediscover its catalytic mechanism. Our results show that the presence of the substrate, the acetoacetate, significantly influences the electrostatic potential of the active site. Furthermore, our simulations show that the decarboxylation reaction can take place by means of a direct proton transfer instead of via an enamine intermediate, which is thought to be strictly necessary. This work provides new insights into the role of the electrostatic interactions on the catalytic activity of AAD and for the first time connects it to the catalytic mechanism of other decarboxylases.

Modulating Secondary Structure Motifs Through Photo-Labile Peptide Staples.

Peptide-protein interactions (PPIs) are facilitated by the well-defined three-dimensional structure of bioactive peptides, interesting compounds for the development of new therapeutic agents. Their secondary structure and thus their propensity to engage in PPIs can be influenced by the introduction of peptide staples on the side chains. In particular, light-controlled staples based on azobenzene photoswitches and their structural influence on helical peptides have been studied extensively. In contrast, photolabile staples bearing photocages as a structural key motif, have mainly been used to block supramolecular interactions. Their influence on the secondary structure of the target peptide is under-investigated. Thus, in this study we use a combination of spectroscopic techniques and in silico simulations to systematically study a series of helical peptides with varying length of the photo-labile staple to obtain a detailed insight into the structure-property relationship in such photoresponsive biomolecules.

Maximized axial helicity in a Pd2L4 cage: inverse guest size-dependent compression and mesocate isomerism.

Helicity is an archetypal structural motif of many biological systems and provides a basis for molecular recognition in DNA. Whilst artificial supramolecular hosts are often helical, the relationship between helicity and guest encapsulation is not well understood. We report a detailed study on a significantly coiled-up Pd2L4 metallohelicate with an unusually wide azimuthal angle (∼176°). Through a combination of NMR spectroscopy, single-crystal X-ray diffraction, trapped ion mobility mass spectrometry and isothermal titration calorimetry we show that the coiled-up cage exhibits extremely tight anion binding (K of up to 106 M-1) by virtue of a pronounced oblate/prolate cavity expansion, whereby the Pd-Pd separation decreases for mono-anionic guests of increasing size. Electronic structure calculations point toward strong dispersion forces contributing to these host-guest interactions. In the absence of a suitable guest, the helical cage exists in equilibrium with a well-defined mesocate isomer that possesses a distinct cavity environment afforded by a doubled Pd-Pd separation distance.

Automated handling of complex chemical structures in Z-matrix coordinates-The chemcoord library.

In this work, we present a fully automated method for the construction of chemically meaningful sets of hierarchical nonredundant internal coordinates (ICs; also commonly denoted as Z-matrices) from the Cartesian coordinates of a molecular system. Particular focus is placed on avoiding ill-definitions of angles and dihedrals due to linear arrangements of atoms, to consistently guarantee a well-defined transformation to Cartesian coordinates, even after structural changes. The representations thus obtained are particularly well suited for pathway construction in double-ended methods for transition state search and optimizations with nonlinear constraints. Analytical gradients for the transformation between the coordinate systems were derived for analytical geometry optimizations purely in Z-matrix coordinates. The geometry optimization was coupled with a Symbolic Algebra package to support arbitrary nonlinear constraints in Z-matrix coordinates, while retaining analytical energy gradient conversion. The difference to the commonly used nonhierarchical IC transformations is discussed. Sample applications are provided for a number of common chemical reactions and illustrative examples.

5-(Trifluorovinyl)dibenzothiophenium Triflate: Introducing the 1,1,2-Trifluoroethylene Tether in Drug-Like Structures.

This manuscript reports the synthesis and structure of an unprecedented sulfonium salt, 5-(trifluorovinyl)dibenzothiophenium triflate, and its use as a versatile reagent for the introduction of the bioisosteric 1,1,2-trifluoroethylene linker in drug-like structures. The protocol developed consists of the reaction of this compound with alcohols and phenols to deliver a complete set of 1,2,2-trifluoro-2-(alkoxy-/aryloxy)ethyl sulfonium salts, which have been purified by column chromatography and fully characterized. Subsequent single electron reduction under mild photochemical conditions efficiently affords the corresponding fluoroalkyl radicals that are trapped either intra- or intermolecularly through their reaction with (hetero)arenes. Theoretical calculations are used to evaluate the conformational consequences derived from the presence of the CF2 -CHF tether.

Ligand-Enabled Disproportionation of 1,2-Diphenylhydrazine at a PV -Center.

We present herein the synthesis of a nearly square-pyramidal chlorophosphorane supported by the tetradentate bis(amidophenolate) ligand, N,N'-bis(3,5-di-tert-butyl-2-phenoxy)-1,2-phenylenediamide. After chloride abstraction the resulting phosphonium cation efficiently promotes the disproportionation of 1,2-diphenylhydrazine to aniline and azobenzene. Mechanistic studies, spectroscopic analyses and theoretical calculations suggest that this unprecedented reactivity mode for PV -centres is induced by the high electrophilicity at the cationic PV -center, which originates from the geometry constraints imposed by the rigid pincer ligand, combined with the ability of the o-amidophenolate moieties to act as electron reservoir. This study illustrates the promising role of cooperativity between redox-active ligands and phosphorus for the design of organocatalysts able to promote redox processes.

Setting up the HyDRA blind challenge for the microhydration of organic molecules.

The procedure leading to the first HyDRA blind challenge for the prediction of water donor stretching vibrations in monohydrates of organic molecules is described. A training set of 10 monohydrates with experimentally known and published water donor vibrations is presented and a test set of 10 monohydrates with unknown or unpublished water donor vibrational wavenumbers is described together with relevant background literature. The rules for data submissions from computational chemistry groups are outlined and the planned publication procedure after the end of the blind challenge is discussed.

Trapping X-ray Radiation Damage from Homolytic Se-C Bond Cleavage in BnSeSeBn Crystals (Bn=benzyl, CH2 C6 H5 ).

Irradiation of dibenzyl diselenide BnSeSeBn with X-ray or UV-light cleaves the Se-C and the Se-Se bonds, inducing stable and metastable radical states. They are inevitably important to all natural and life sciences. Structural changes due to X-ray-induced Se-C bond-cleavage could be pin-pointed in various high-resolution X-ray diffraction experiments for the first time. Extended DFT methods were applied to characterize the solid-state structure and support the refinement of the observed residuals as contributions from the BnSeSe⋅ radical species. The X-ray or UV-irradiated crystalline samples of BnSeSeBn were characterized by solid-state EPR. This paper provides insight that in the course of X-ray structure analysis of selenium compounds not only organo-selenide radicals like RSe⋅ may occur, but also organo diselenide BnSeSe⋅ radicals and organic radicals R⋅ are generated, particularly important to know in structural biology.