Computational Studies of Reactions of 1,2,4,5-Tetrazines with Enamines in MeOH and HFIP.

The reaction between 1,2,4,5-tetrazines and alkenes in polar solvents proceeds through a Diels-Alder cycloaddition along the C-C axis (C3/C6 cycloaddition) of the tetrazine, followed by dinitrogen loss. By contrast, the reactions of 1,2,4,5-tetrazines with enamines in hexafluoroisopropanol (HFIP) give 1,2,4-triazine products stemming from a formal Diels-Alder addition across the N-N axis (N1/N4 cycloaddition). We explored the mechanism of this interesting solvent effect through DFT calculations in detail and revealed a novel reaction pathway characterized by C-N bond formation, deprotonation, and a 3,3-sigmatropic rearrangement. The participation of an HFIP molecule was found to be crucial to the N1/N4 selectivity over C3/C6 due to the more favored initial C-N bond formation than C-C bond formation.

Cycloaddition Reactivities Analyzed by Energy Decomposition Analyses and the Frontier Molecular Orbital Model.

This Account describes our quest to understand and predict organic reactivity, a principal goal of physical and theoretical organic chemistry. The focus is on the development and testing of models for the prediction of cycloaddition reactivities and selectivities. We describe the involvement of the Houk group, and other groups, in the evolution of theoretical models that can achieve ever greater accuracy as well as provide practical heuristic models for understanding and prediction.Is the venerable frontier molecular orbital (FMO) model, the basis of Kenichi Fukui's 1981 Nobel Prize, still useful, or must it be replaced with more advanced models? In particular, models such as Conceptual Density Functional, the Pauli Exclusion Model, and the recent popularity of Electrostatic Potential Plots and Dispersion Energies have not only added to our understanding, but they have also created uncertainty about whether the simple FMO heuristic model has a place in 21st century discussions. This Account addresses this issue and asserts the value of the FMO model.Beginning with brief descriptions of selected models for cycloaddition reactivity starting with early donor-acceptor (nucleophile-electrophile) charge-transfer concepts, this Account reviews Fukui's frontier molecular orbital model, Salem and Klopman's orbital, electrostatic and Pauli repulsion model, the conceptual DFT model by Parr and later by Domingo and others, the recent Houk and Bickelhaupt Distortion/Interaction Activation Strain model, and the Bickelhaupt-Hamlin's Pauli-repulsion lowering model.Computations and analyses of four well-studied Diels-Alder cycloadditions, both normal and inverse electron-demand types, are presented. Most were studied earlier in our published work but are presented here with new insights from calculations with modern methods. Depending on the types of substrates (cycloaddends), the dominant factors controlling reactivity can be orbital interactions, electrostatics and polarization, or Pauli repulsion and dispersion effects, or a combination of all of these.By comparing orbital interactions, especially the frontier molecular orbital interactions, with the other factors that influence reactivity, we show why the FMO model is such a powerful─and theoretically meaningful─heuristic for understanding and predicting reactivity. We also present a method to understand Pauli repulsion effects on activation barriers, ρ(1.1). The use of a new reaction coordinate, the extent of Pauli repulsion along the reaction path, is advocated to emphasize the role of repulsive occupied orbital interactions on reactivity.Fukui's frontier molecular orbital model is effective because FMO interactions parallel all the quantities that influence reactivity. The FMO model continues to provide a practical model to understand and guide experiments.

"Holographic" Autostereoscopic Displays: A Perspective on Their Technology and Potential Impact in Chemistry.

The three-dimensional arrangement of atoms in molecules is essential for understanding their properties and behavior. Traditional 2D representations and digital 3D models presented on 2D media often fall short in conveying the complexity of molecular structures. Autostereoscopic displays, often marketed as "holographic" displays, pose a potential solution to this challenge. These displays, with their multi-view and single-view configurations, promise to advance chemistry education and research by offering accurate 3D representations with depth and parallax. In this perspective, I delve into the possibilities and limitations of autostereoscopic displays in chemistry, discussing the underlying technology and potential applications, from research to teaching and science communication. Multi-view autostereoscopic displays excel in facilitating collaborative work by enabling multiple viewers to simultaneously perceive the same 3D structure from different angles. However, they currently suffer from low resolution and high cost, which could limit their immediate widespread adoption. Conversely, single-view autostereoscopic displays with eye-tracking, while limited to one viewer at a time, provide higher resolution at a lower cost, thus suggesting that they might become the technology of the future given the balance of price to performance. Despite current limitations, autostereoscopic displays possess undeniable potential for shaping the future of chemistry education and research.

Substituent Effects in Bioorthogonal Diels-Alder Reactions of 1,2,4,5-Tetrazines.

1,2,4,5-Tetrazines are increasingly used as reactants in bioorthogonal chemistry due to their high reactivity in Diels-Alder reactions with various dienophiles. Substituents in the 3- and 6-positions of the tetrazine scaffold are known to have a significant impact on the rate of cycloadditions; this is commonly explained on the basis of frontier molecular orbital theory. In contrast, we show that reactivity differences between commonly used classes of tetrazines are not controlled by frontier molecular orbital interactions. In particular, we demonstrate that mono-substituted tetrazines show high reactivity due to decreased Pauli repulsion, which leads to a more asynchronous approach associated with reduced distortion energy. This follows the recent Vermeeren-Hamlin-Bickelhaupt model of reactivity increase in asymmetric Diels-Alder reactions. In addition, we reveal that ethylene is not a good model compound for other alkenes in Diels-Alder reactions.

Oxidative Desymmetrization Enables the Concise Synthesis of a trans-Cyclooctene Linker for Bioorthogonal Bond Cleavage.

Modified trans-cyclooctenes (TCO) are capable of highly efficient molecular manipulations in biological environments, driven by the bioorthogonal reaction with tetrazines (Tz). The development of click-cleavable TCO has fueled the field of in vivo chemistry and enabled the design of therapeutic strategies that have already started to enter the clinic. A key element for most of these approaches is the implementation of a cleavable TCO linker. So far, only one member of this class has been developed, a compound that requires a high synthetic effort, mainly to fulfill the multilayered demands on its chemical structure. To tackle this limitation, we developed a dioxolane-fused cleavable TCO linker (dcTCO) that can be prepared in only five steps by applying an oxidative desymmetrization to achieve diastereoselective introduction of the required functionalities. Based on investigation of the structure, reaction kinetics, stability, and hydrophilicity of dcTCO, we demonstrate its bioorthogonal application in the design of a caged prodrug that can be activated by in-situ Tz-triggered cleavage to achieve a remarkable >1000-fold increase in cytotoxicity.

To Bond or Not to Bond: Metal-Metal Interaction in Heterobimetallic Rare-Earth Metal-Silver Complexes.

The reaction of 2,4-tBu2-6-(PPh2)PhOH (HOArP) with silver(I) triflate in a 3:1 molar ratio gave the mononuclear coinage metal complex (HOArP-κP)3AgIOTf (1). Treatment of HOArP with LnIII[N(SiMe3)2]3 (Ln = La, Sm, Y, Yb) in a 3:1 molar ratio yielded the mononuclear rare-earth metal complexes LnIII(OArP-κ2O,P)3 (2-Ln). The heterobimetallic rare-earth metal-silver complexes LnIII(OTf)(µ-OArP-1κ1O,2κ1P)3AgI (3-Ln) were prepared from monometallic precursors by reactions of equimolar amounts of 1 with LnIII[N(SiMe3)2]3 or 2-Ln with silver(I) triflate, respectively. The compounds were characterized by NMR, ultraviolet-visible (UV-vis), and infrared (IR) spectroscopy, single-crystal X-ray diffraction, elemental analysis, and the effective magnetic moments of the paramagnetic complexes were determined via the Evans NMR method. Computational studies were conducted on 3-La and 3-Y.

Computational generation of an annotated gigalibrary of synthesizable, composite peptidic macrocycles.

Peptidomimetic macrocycles have the potential to regulate challenging therapeutic targets. Structures of this type having precise shapes and drug-like character are particularly coveted, but are relatively difficult to synthesize. Our laboratory has developed robust methods that integrate small-peptide units into designed scaffolds. These methods create macrocycles and embed condensed heterocycles to diversify outcomes and improve pharmacological properties. The hypothetical scope of the methodology is vast and far outpaces the capacity of our experimental format. We now describe a computational rendering of our methodology that creates an in silico three-dimensional library of composite peptidic macrocycles. Our open-source platform, CPMG (Composite Peptide Macrocycle Generator), has algorithmically generated a library of 2,020,794,198 macrocycles that can result from the multistep reaction sequences we have developed. Structures are generated based on predicted site reactivity and filtered on the basis of physical and three-dimensional properties to identify maximally diverse compounds for prioritization. For conformational analyses, we also introduce ConfBuster++, an RDKit port of the open-source software ConfBuster, which allows facile integration with CPMG and ready parallelization for better scalability. Our approach deeply probes ligand space accessible via our synthetic methodology and provides a resource for large-scale virtual screening.

Mechanistic Insights into the Reaction of Amidines with 1,2,3-Triazines and 1,2,3,5-Tetrazines.

1,2,3-Triazines and 1,2,3,5-tetrazines react rapidly, efficiently, and selectively with amidines to form pyrimidines/1,3,5-triazines, exhibiting an orthogonal reactivity with 1,2,4,5-tetrazine-based conjugation chemistry. Whereas the mechanism of the reaction of the isomeric 1,2,4-triazines and 1,2,4,5-tetrazines with alkenes is well understood, the mechanism of the 1,2,3-triazine/1,2,3,5-tetrazine-amidine reaction as well as its intrinsic reactivity remains underexplored. By using 15N-labeling, kinetic investigations, and kinetic isotope effect studies, complemented by extensive computational investigations, we show that this reaction proceeds through an addition/N2 elimination/cyclization pathway, rather than the generally expected concerted or stepwise Diels-Alder/retro Diels-Alder sequence. The rate-limiting step in this transformation is the initial nucleophilic attack of an amidine on azine C4, with a subsequent energetically favored N2 elimination step compared with a disfavored stepwise formation of a Diels-Alder cycloadduct. The proposed reaction mechanism is in agreement with experimental and computational results, which explains the observed reactivity of 1,2,3-triazines and 1,2,3,5-tetrazines with amidines.

Uncovering the Key Role of Distortion in Bioorthogonal Tetrazine Tools That Defy the Reactivity/Stability Trade-Off.

The tetrazine/trans-cyclooctene ligation stands out from the bioorthogonal toolbox due to its exceptional reaction kinetics, enabling multiple molecular technologies in vitro and in living systems. Highly reactive 2-pyridyl-substituted tetrazines have become state of the art for time-critical processes and selective reactions at very low concentrations. It is widely accepted that the enhanced reactivity of these chemical tools is attributed to the electron-withdrawing effect of the heteroaryl substituent. In contrast, we show that the observed reaction rates are way too high to be explained on this basis. Computational investigation of this phenomenon revealed that distortion of the tetrazine caused by intramolecular repulsive N-N interaction plays a key role in accelerating the cycloaddition step. We show that the limited stability of tetrazines in biological media strongly correlates with the electron-withdrawing effect of the substituent, while intramolecular repulsion increases the reactivity without reducing the stability. These fundamental insights reveal thus far overlooked mechanistic aspects that govern the reactivity/stability trade-off for tetrazines in physiologically relevant environments, thereby providing a new strategy that may facilitate the rational design of these bioorthogonal tools.

Synergistic Experimental and Computational Investigation of the Bioorthogonal Reactivity of Substituted Aryltetrazines.

Tetrazines (Tz) have been applied as bioorthogonal agents for various biomedical applications, including pretargeted imaging approaches. In radioimmunoimaging, pretargeting increases the target-to-background ratio while simultaneously reducing the radiation burden. We have recently reported a strategy to directly 18F-label highly reactive tetrazines based on a 3-(3-fluorophenyl)-Tz core structure. Herein, we report a kinetic study on this versatile scaffold. A library of 40 different tetrazines was prepared, fully characterized, and investigated with an emphasis on second-order rate constants for the reaction with trans-cyclooctene (TCO). Our results reveal the effects of various substitution patterns and moreover demonstrate the importance of measuring reactivities in the solvent of interest, as click rates in different solvents do not necessarily correlate well. In particular, we report that tetrazines modified in the 2-position of the phenyl substituent show high intrinsic reactivity toward TCO, which is diminished in aqueous systems by unfavorable solvent effects. The obtained results enable the prediction of the bioorthogonal reactivity and thereby facilitate the development of the next generation of substituted aryltetrazines for in vivo applications.

Origin of Increased Reactivity in Rhenium-Mediated Cycloadditions of Tetrazines.

Pyridyl tetrazines coordinated to metals like rhenium have been shown to be more reactive in [4 + 2] cycloadditions than their uncomplexed counterparts. Using density functional theory calculations, we found a more favorable interaction energy caused by stronger orbital interactions as the origin of this increased reactivity. Additionally, the high regioselectivity is due to a greater degree of charge stabilization in the transition state, leading to the major product.

How the Lewis Base F- Catalyzes the 1,3-Dipolar Cycloaddition between Carbon Dioxide and Nitrilimines.

The mechanism of the Lewis base F- catalyzed 1,3-dipolar cycloaddition between CO2 and nitrilimines is interrogated using DFT calculations. F- activates the nitrilimine, not CO2 as proposed in the literature, and imparts a significant rate enhancement for the cycloaddition. The origin of this catalysis is in the strength of the primary orbital interactions between the reactants. The Lewis base activated nitrilimine-F- has high-lying filled FMOs. The smaller FMO-LUMO gap promotes a rapid nucleophilic attack and overall cycloaddition with CO2.

Computational Exploration of Ambiphilic Reactivity of Azides and Sustmann's Paradigmatic Parabola.

We examine the theoretical underpinnings of the seminal discoveries by Reiner Sustmann about the ambiphilic nature of Huisgen's phenyl azide cycloadditions. Density functional calculations with ωB97X-D and B2PLYP-D3 reproduce the experimental data and provide insights into ambiphilic control of reactivity. Distortion/interaction-activation strain and energy decomposition analyses show why Sustmann's use of dipolarophile ionization potential is such a powerful predictor of reactivity. We add to Sustmann's data set several modern distortion-accelerated dipolarophiles used in bioorthogonal chemistry to show how these fit into the orbital energy criteria that are often used to understand cycloaddition reactivity. We show why such a simple indicator of reactivity is a powerful predictor of reaction rates that are actually controlled by a combination of distortion energies, charge transfer, closed-shell repulsion, polarization, and electrostatic effects.

Origins of Endo Selectivity in Diels-Alder Reactions of Cyclic Allene Dienophiles.

Strained cyclic allenes, first discovered in 1966 by Wittig and co-workers, have recently emerged as valuable synthetic building blocks. Previous experimental investigations, and computations reported here, demonstrate that the Diels-Alder reactions of furans and pyrroles with 1,2-cyclohexadiene and oxa- and azaheterocyclic analogs proceed with endo selectivity. This endo selectivity gives the adduct with the allylic saturated carbon of the cyclic allene endo to the diene carbons. The selectivity is very general and useful in synthetic applications. Our computational study establishes the origins of this endo selectivity. We analyze the helical frontier molecular orbitals of strained cyclic allenes and show how secondary orbital and electrostatic effects influence stereoselectivity. The LUMO of carbon-3 of the allene (C-3 is not involved in primary orbital interactions) interacts in a stabilizing fashion with the HOMO of the diene in such a way that the carbon of the cyclic allene attached to C-1 favors the endo position in the transition state. The furan LUMO, allene HOMO interaction reinforces this preference. These mechanistic studies are expected to prompt the further use of long-avoided strained cyclic allenes in chemical synthesis.

Live Monitoring of Strain-Promoted Azide Alkyne Cycloadditions in Complex Reaction Environments by Inline ATR-IR Spectroscopy.

The strain-promoted azide alkyne cycloaddition (SPAAC) is a powerful tool for forming covalent bonds between molecules even under physiological conditions, and therefore found broad application in fields ranging from biological chemistry and biomedical research to materials sciences. For many applications, knowledge about reaction kinetics of these ligations is of utmost importance. Kinetics are commonly assessed and studied by NMR measurements. However, these experiments are limited in terms of temperature and restricted to deuterated solvents. By using an inline ATR-IR probe we show that the cycloaddition of azides and alkynes can be monitored in aqueous and even complex biological fluids enabling the investigation of reaction kinetics in various solvents and even human blood plasma under controlled conditions in low reaction volumes.

Concerted [4 + 2] and Stepwise (2 + 2) Cycloadditions of Tetrafluoroethylene with Butadiene: DFT and DLPNO-UCCSD(T) Explorations.

Tetrafluoroethylene and butadiene form the 2 + 2 cycloadduct under kinetic control, but the Diels-Alder cycloadduct is formed under thermodynamic control. Borden and Getty showed that the preference for 2 + 2 cycloaddition is due to the necessity for syn-pyramidalization of the two CF2 groups in the 4 + 2 transition state. We have explored the full potential energy surface for the concerted and stepwise reactions of tetrafluoroethylene and butadiene with density functional theory, DFT (B3LYP and M06-2X), DLPNO-UCCSD(T), and CASSCF-NEVPT2 methods and with the distortion/interaction-activation strain model to explain the energetics of different pathways. The 2 + 2 cycloadduct is formed by an anti-transition state followed by two rotations and a final bond formation transition state. Energetics are compared to the reaction of maleic anhydride and ethylene.

Secondary Orbital Interactions Enhance the Reactivity of Alkynes in Diels-Alder Cycloadditions.

We have investigated the inverse electron-demand Diels-Alder reactions of trans-cyclooctene (TCO) and endo-bicyclo[6.1.0]nonyne (BCN) with a 1,2,4,5-tetrazine, a cyclopentadienone, and an ortho-benzoquinone. Tetrazines react significantly faster with TCO compared to BCN because the highest occupied molecular orbital (HOMO) of TCO is significantly higher in energy than the HOMO of BCN and there is less distortion of the tetrazine. Despite the different HOMO energies, TCO and BCN have similar reactivities toward cyclopentadienones, while BCN is significantly more reactive than TCO in the cycloaddition with ortho-benzoquinone. We find that the higher reactivity of BCN compared to TCO with ortho-benzoquinone is due to secondary orbital interactions of the BCN HOMO-1 with the diene LUMO.

autoDIAS: a python tool for an automated distortion/interaction activation strain analysis.

The distortion/interaction activation strain (DIAS) analysis is a powerful tool for the investigation of energy barriers. However, setup and data analysis of such a calculation can be cumbersome and requires lengthy intervention of the user. We present autoDIAS, a python tool for the automated setup, performance, and data extraction of the DIAS analysis, including automated detection of fragments and relevant geometric parameters. © 2019 Wiley Periodicals, Inc.

Chemoselectivity of Tertiary Azides in Strain-Promoted Alkyne-Azide Cycloadditions.

The strain-promoted alkyne-azide cycloaddition (SPAAC) is the most commonly employed bioorthogonal reaction with applications in a broad range of fields. Over the years, several different cyclooctyne derivatives have been developed and investigated in regard to their reactivity in SPAAC reactions with azides. However, only a few studies examined the influence of structurally diverse azides on reaction kinetics. Herein, we report our investigations of the reactivity of primary, secondary, and tertiary azides with the cyclooctynes BCN and ADIBO applying experimental and computational methods. All azides show similar reaction rates with the sterically non-demanding cyclooctyne BCN. However, due to the increased steric demand of the dibenzocyclooctyne ADIBO, the reactivity of tertiary azides drops by several orders of magnitude in comparison to primary and secondary azides. We show that this chemoselective behavior of tertiary azides can be exploited to achieve semiorthogonal dual-labeling without the need for any catalyst using SPAAC exclusively.

HPMA-Based Nanoparticles for Fast, Bioorthogonal iEDDA Ligation.

Fast and bioorthogonally reacting nanoparticles are attractive tools for biomedical applications such as tumor pretargeting. In this study, we designed an amphiphilic block copolymer system based on HPMA using different strategies to introduce the highly reactive click units 1,2,4,5-tetrazines (Tz) either at the chain end (Tz-CTA) or statistical into the hydrophobic block. This reactive group undergoes a rapid, bioorthogonal inverse electron-demand Diels-Alder reaction (iEDDA) with trans-cyclooctenes (TCO). Subsequently, this polymer platform was used for the preparation of different Tz-covered nanoparticles, such as micelles and colloids. Thereby it was found that the reactivity of the polymeric micelles is comparable to that of the low molar mass tetrazines. The core-cross-linked micelles can be successfully conjugated at rather low concentrations to large biomacromolecules like antibodies, not only in physiological buffer, but also in human blood plasma, which was confirmed by fluorescence correlation spectroscopy (FCS).