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

Pnictogen bonding in imide derivatives for chiral folding and self-assembly

Pnictogen bonding (PnB) is an attraction interaction that originates from the anisotropic distribution of electron density of pnictogen elements, which however has been rarely found in nitrogen atoms. In this work, for the first time, we unveil the general presence of N-involved PnB in aromatic or aliphatic imide groups and reveal its implications in chiral self-assembly of folding. This long-neglected interaction was consolidated by Cambridge structural database (CSD) searching as well as subsequent computational studies. Though the presence of PnB has limited effects on spectroscopic properties in the solution phase, conformation locking effects are sufficiently expressed in the chiral folding and self-assembly behavior. PnB anchors the chiral conformation to control the emergence and inversion of chiroptical signals, while intramolecular PnB induces the formation of supramolecular tilt chirality. It also enables the chiral folding of imide-containing amino acid or peptide derivatives, which induces the formation of unique secondary structural sequences such as β-sheets. Finally, the effects of PnB in directing folded helical structures were revealed. Examples of cysteine and cystine derivatives containing multiple N⋯O and N⋯S PnBs constitute an α-helix like secondary structure with characteristic circular dichroism. This work discloses the comprehensive existence of imide-involved PnB, illustrates its important role in folding and self-assembly, and sheds light on the rational fabrication of conformation-locked compounds and polymers with controllable chiroptical activities.

A Deeper Insight into the Supramolecular Activation of Oxidative Addition Reactions Involving Pincer-Rhodium(I) Complexes

The factors governing the acceleration of the oxidative addition of methyl iodide to pincer rhodium(I)-complexes induced by coronene have been computationally explored in detail using quantum chemical methods. Both the parent reaction and the coronene-mediated process proceed via a stepwise SN2-type mechanism. It is found that the acceleration of the process derives from the formation of an initial supramolecular complex, mainly stabilized by electrostatic and π-π interactions, which significantly increases the electron richness of the complex. The impact of this effect on the reaction barrier has been quantitatively analyzed by applying the activation strain model in combination with the energy decomposition analysis method. In addition, the influence of other polycyclic aromatic hydrocarbons on the oxidative reaction has been also considered.

Computational study of ground-state properties of μ2 -bridged group 14 porphyrinic sandwich complexes

The structural properties of μ₂ -bridged porphyrinic double-decker complexes are investigated and the influence of various ligands, metals, substituents, and bridging atoms on the dominant structural motif is elucidated. A variety of quantum chemical methods including semiempirical (SQM) methods and density functional theory (DFT) is assessed for the calculation of ecliptic and staggered conformational energies. Local coupled cluster (DLPNO-CCSD(T1)) data are generated for reference. The r² SCAN-3c composite scheme as well as the B2PLYP-D4/def2-QZVPP approach are identified as reliable methods. Energy decomposition analyses (EDA) and localized molecular orbital analyses (LMO) are used to investigate the bonding situation and the nature of the inter-ligand interaction energy underlining the crucial role of attractive London dispersion interactions. Targeted modification of the bridging atom, e.g., by replacing O^2- by S^2- is shown to drastically change the major structural features of the investigated complexes. Further, the influence of different substituents of varying size at the phthalocyanine ligand regarding the dominant conformation is described.

Involvement of Arsenic Atom of AsF3 in Five Pnicogen Bonds: Differences between X-ray Structure and Theoretical Models

Bonding within the AsF₃ crystal is analyzed via quantum chemical methods so as to identify and quantify the pnicogen bonds that are present. The structure of a finite crystal segment containing nine molecules is compared with that of a fully optimized cluster of the same size. The geometries are qualitatively different, with a much larger binding energy within the optimized nonamer. Although the total interaction energy of a central unit with the remaining peripheral molecules is comparable for the two structures, the binding of the peripherals with one another is far larger in the optimized cluster. This distinction of much stronger total binding within the optimized cluster is not limited to the nonamer but repeats itself for smaller aggregates as well. The average binding energy of the cluster rises quickly with size, asymptotically approaching a value nearly triple that of the dimer.

Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods

The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S₀/T₁ minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S₀/S₁ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.

Towards understanding solvation effects on the conformational entropy of non-rigid molecules

The absolute molecular entropy is a fundamental quantity for the accurate description of thermodynamic properties. For non-rigid molecules, a substantial part of the entropy can be attributed to a conformational contribution. Systems and properties where this is relevant, e.g., protein-ligand binding affinities or pK_a values refer usually to the liquid phase. In this work, the influence of solvation on the conformational entropy is investigated. A recently introduced state-of-the-art and automated computational protocol for the computation of conformational entropies [Pracht et al., Chem. Sci., 2021, 12, 6551-6568.] is applied in combination with fast and accurate semiempirical quantum-chemical methods and implicit solvation models for a set of 25 commercially available drug molecules and five transition metal compounds. Computed gas-phase conformational entropies are compared with values obtained in implicit n-hexane and water. It is found that implicit solvation can have a substantial effect of several cal mol^-1 K^-1 on the entropy as a result of large conformational changes in the different phases. We conclude that for flexible molecules chemical accuracy for free energies in solution can only be achieved if solvation effects on the conformational ensemble are considered.

Surprising Chemistry of 6-Azidotetrazolo[5,1-a]phthalazine: What a Purported Natural Product Reveals about the Polymorphism of Explosives

6-Azidotetrazolo[5,1-a]phthalazine (ATPH) is a nitrogen-rich compound of surprisingly broad interest. It is purported to be a natural product, yet it is closely related to substances developed as explosives and is highly polymorphic despite having a nearly planar structure with little flexibility. Seven solid forms of ATPH have been characterized by single-crystal X-ray diffraction. The structures show diverse patterns of molecular organization, including both stacked sheets and herringbone packing. In all cases, N···N and C-H···N interactions play key roles in ensuring molecular cohesion. The high polymorphism of ATPH appears to arise in part from the ability of virtually every atom of nitrogen and hydrogen in the molecule to take part in close N···N and C-H···N contacts. As a result, adjacent molecules can adopt many different relative orientations that are energetically similar, thereby generating a polymorphic landscape with an unusually high density of potential structures. This landscape has been explored in detail by the computational prediction of crystal structures. Studying ATPH has provided insights into the field of energetic materials, where access to multiple polymorphs can be used to improve performance and clarify how it depends on molecular packing. In addition, our work with ATPH shows how valuable insights into molecular crystallization, often gleaned from statistical analyses of structural databases, can also come from in-depth empirical and theoretical studies of single compounds that show distinctive behavior.

Synthesis, Experimental and Theoretical Study of Azidochromones

A series of 2-(haloalkyl)-3-azidomethyl and 6-azido chromones has been synthetized, characterized and studied by theoretical (DFT calculations) and spectroscopic methods (UV-Vis, NMR). The crystal structure of 3-azidomethyl-2-difluoromethyl chromone, determined by X-ray diffraction methods, shows a planar framework due to extended π-bond delocalization. Its molecular packing is stabilized by F···H, N···H and O···H hydrogen bonds, π···π stacking and C–O···π intermolecular interactions. Moreover, AIM, NCI and Hirshfeld analysis evidenced that azido moiety has a significant role in the stabilization of crystal packing through weak intermolecular interactions, where analysis of electronic density suggested closed-shell (CS) interatomic interactions.

DeepNCI: DFT Noncovalent Interaction Correction with Transferable Multimodal Three-Dimensional Convolutional Neural Networks

A multimodal deep learning model, DeepNCI, is proposed for improving noncovalent interactions (NCIs) calculated via density functional theory (DFT). DeepNCI is composed of a three-dimensional convolutional neural network (3D CNN) for abstracting critical and comprehensive features from 3D electron density, and a neural network for modeling one-dimensional quantum chemical properties. By merging features from two networks, DeepNCI is able to reduce the root-mean-square error of DFT-calculated NCI from 1.19 kcal/mol to ∼0.2 kcal/mol for a NCI molecular database (>1000 molecules). The representativeness of the joint features can be visualized by t-distributed stochastic neighbor embedding (t-SNE), where they can distinguish categorized NCI systems quite well. Therefore, the fused model performs better than its component networks. In addition, the 3D CNN takes electron density as inputs that are in the same range, despite the size of molecular systems, so it can promote model applicability and transferability. To clarify the applicability of DeepNCI, an application domain (AD) has been defined with merged features using the K-nearest-neighbor method. The calculations for external test sets are shown that AD can properly monitor the reliability for a prediction. The model transferability is tested with a small database of homolysis bond dissociation energy including only dozens of samples. With NCI database pretrained parameters, the same or better performance than the reported results is achieved by transfer learning. This suggests that the DeepNCI model is transferable and it may transfer to other relative tasks, which possibly can resolve some small sampling problems. The source code of DeepNCI can be freely accessed at https://github.com/wenzelee/DeepNCI.

London Dispersion Helps Refine Steric A-Values: Dispersion Energy Donor Scales

We suggest a scale of dispersion energy donors (DEDs) that allows for direct comparisons with steric effects. This scale is based on the classic A-values and allows groups to reorient to minimize strain, thereby providing an advantage over raw group polarizabilities. The A-value can no longer be considered purely a steric factor. Even for groups that do not participate in charge transfer or electrostatic interactions, the A-value includes Pauli repulsion (steric hindrance) and attractive London dispersion (LD) interactions. Although the common assumption is that, at the distances found in monosubstituted cyclohexanes, steric demands are the key factors influencing conformer preferences, we show in this computational study that there is a non-negligible LD part. We use this system to build a DED scale and a complementary steric scale. These scales are quantitatively comparable, as they are based on the same system, and allow for comparison of the two competing interactions in experimentally relevant settings. In addition, we show that LD interactions can be used to explain puzzling data regarding relative group sizes.

cis-3-Azido-2-methoxyindolines as safe and stable precursors to overcome the instability of fleeting 3-azidoindoles

Use of 3-azidoindoles in organic synthesis remains a difficult task owing to their instabilities. Herein, we report a general and concise approach for tackling this problem by using 3-azidoindole surrogates. The surrogates are bench-stable, presumably due to the observed intramolecular O-N_β bonding. The resultant fleeting intermediates undergo capturing in situ to afford 3-substitued indoles through formal ipso-substitution of the azide group by nucleophiles. In these investigations, we found that the fleeting 3-azidoindoles show a C3-electrophilic character for the first time.

Preparation of Nitrogen Analogues of Ceramide and Studies of Their Aggregation in Sphingomyelin Bilayers

Ceramides can regulate biological processes probably through the formation of laterally segregated and highly packed ceramide-rich domains in lipid bilayers. In the course of preparation of its analogues, we found that a hydrogen-bond-competent functional group in the C1 position is necessary to form ceramide-rich domains in lipid bilayers [Matsufuji; Langmuir 2018]. Hence, in the present study, we newly synthesized three ceramide analogues: CerN₃, CerNH₂, and CerNHAc, in which the 1-OH group of ceramide is substituted with a nitrogen functionality. CerNH₂ and CerNHAc are capable of forming hydrogen bonds in their headgroups, whereas CerN₃ is not. Fluorescent microscopy observation and differential scanning calorimetry analysis disclosed that these ceramide analogues formed ceramide-rich phases in sphingomyelin bilayers, although their thermal stability was slightly inferior to that of normal ceramides. Moreover, wide-angle X-ray diffraction analysis showed that the chain packing structure of ceramide-rich phases of CerNHAc and CerN₃ was similar to that of normal ceramide, while the CerNH₂-rich phase showed a slightly looser chain packing due to the formation of CerNH₃⁺. Although the domain formation of CerN₃ was unexpected because of the lack of hydrogen-bond capability in the headgroup, it may become a promising tool for investigating the mechanistic link between the ceramide-rich phase and the ceramide-related biological functions owing to its Raman activity and applicability to click chemistry.

Azide···Oxygen Interaction: A Crystal Engineering Tool for Conformational Locking

We have designed, synthesized and crystallized 36 compounds, each containing an azide group and an oxygen atom separated by three bonds. Crystal structure analysis revealed that each of these molecules adopts a conformation in which the azide and oxygen groups orient syn to each other with a short O ··· N b  contact. Geometry-optimized structures [using  M06-2X/6-311G(d,p) level of theory ] also showed the syn conformation in all 36 of these cases, suggesting that this not merely a crystal packing effect. Quantum topological analysis using Bader's Atoms in Molecules (AIM) theory revealed bond paths and bond critical points (BCP) in these structures suggesting its nature and energetics to be similar to weak hydrogen bonding. The NCI-RDG plot clearly revealed the attractive interaction consisting of electrostatic or dispersive components in all the 36 systems. NBO analysis suggested a weak orbital-relaxation (charge-transfer) contribution of energy for a few  (sp2) O-donor systems. Natural population analysis (NPA) and molecular electrostatic potential mapping (MESP) of these crystal structures further revealed the existence of favorable azide-oxygen interaction.   A CSD search indicated the frequent and consistent occurrence of this interaction and its role dictating the syn conformation of azide and oxygen in molecules where these groups are separated by 2-4 bonds.

From QCEIMS to QCxMS: A Tool to Routinely Calculate CID Mass Spectra Using Molecular Dynamics

Mass spectrometry (MS) is a powerful tool in chemical research and substance identification. For the computational modeling of electron ionization MS, we have developed the quantum-chemical electron ionization mass spectra (QCEIMS) program. Here, we present an extension of QCEIMS to calculate collision-induced dissociation (CID) spectra. The more general applicability is accounted for by the new name QCxMS, where "x" refers to EI or CID. To this end, fragmentation and rearrangement reactions are computed "on-the-fly" in Born-Oppenheimer molecular dynamics (MD) simulations with the semiempirical GFN2-xTB Hamiltonian, which provides an efficient quantum mechanical description of all elements up to Z = 86 (Rn). Through the explicit modeling of multicollision processes between precursor ions and neutral gas atoms as well as temperature-induced decomposition reactions, QCxMS provides detailed insight into the collision kinetics and fragmentation pathways. In combination with the CREST program to determine the preferential protonation sites, QCxMS becomes the first standalone MD-based program that can predict mass spectra based solely on molecular structures as input. We demonstrate this for six organic molecules with masses ranging from 159 to 296 Da, for which QCxMS yields CID spectra in reasonable agreement with experiments.

Comment on "The Nature of Chalcogen-Bonding-Type Tellurium-Nitrogen Interactions": Fixing the Description of Finite-Temperature Effects Restores the Agreement Between Experiment and Theory

Mitzel and co-workers recently presented an intriguing molecule displaying a tellurium-nitrogen interaction. Structural data obtained in the solid and in gas phase indicated a large increase of the Te-N equilibrium distance re from 2.64 to 2.92 Å, respectively. Although some DFT calculations appear to support the large re in gas phase, we argue that the lions share of the increase is due to an incomplete description of finite-temperature effects in the back-corrected experimental data. This hypothesis is based on high-level coupled-cluster (CC) and periodic DFT calculations, which consistently point towards a much smaller re in the isolated molecule. Further support comes through MD simulations with a tuned GFN2-xTB Hamiltonian: Calibrated against a CC reference, these show a six-times larger influence of temperature than with the originally used GFN1-xTB. Taking this into account, the back-corrected re in gas phase becomes 2.67±0.08 Å, in good agreement with high-level CC theory and most DFT methods.

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⁻¹.