SNAr or Sulfonylation? Chemoselective Amination of Halo(het)arene Sulfonyl Halides for Synthetic Applications and Ultralarge Compound Library Design.

The chemoselectivity of halo(het)arene sulfonyl halide aminations is studied thoroughly under parallel synthesis conditions, and the scope and limitations of the method are established. It is shown that SNAr-reactive sulfonyl halides typically undergo sulfonamide synthesis during the first step; the second amination is also possible provided that the SNAr-active center is sufficiently reactive. On the contrary, sulfonyl fluorides bearing an arylating moiety undergo selective transformation at the latter reactive center under proper control. Further sulfur-fluoride exchange (SuFEx) is also possible, which can be especially valuable for some sulfonyl halide classes. The developed two-step parallel double amination protocol provides access to a 6.67-billion compound synthetically tractable REAL-type chemical space (76% expected synthesis success rate).

When Yield Prediction Does Not Yield Prediction: An Overview of the Current Challenges.

Machine Learning (ML) techniques face significant challenges when predicting advanced chemical properties, such as yield, feasibility of chemical synthesis, and optimal reaction conditions. These challenges stem from the high-dimensional nature of the prediction task and the myriad essential variables involved, ranging from reactants and reagents to catalysts, temperature, and purification processes. Successfully developing a reliable predictive model not only holds the potential for optimizing high-throughput experiments but can also elevate existing retrosynthetic predictive approaches and bolster a plethora of applications within the field. In this review, we systematically evaluate the efficacy of current ML methodologies in chemoinformatics, shedding light on their milestones and inherent limitations. Additionally, a detailed examination of a representative case study provides insights into the prevailing issues related to data availability and transferability in the discipline.

2-Oxabicyclo[2.2.2]octane as a new bioisostere of the phenyl ring.

The phenyl ring is a basic structural element in chemistry. Here, we show the design, synthesis, and validation of its new saturated bioisostere with improved physicochemical properties - 2-oxabicyclo[2.2.2]octane. The design of the structure is based on the analysis of the advantages and disadvantages of the previously used bioisosteres: bicyclo[1.1.1]pentane, bicyclo[2.2.2]octane, and cubane. The key synthesis step is the iodocyclization of cyclohexane-containing alkenyl alcohols with molecular iodine in acetonitrile. 2-Oxabicyclo[2.2.2]octane core is incorporated into the structure of Imatinib and Vorinostat (SAHA) drugs instead of the phenyl ring. In Imatinib, such replacement leads to improvement of physicochemical properties: increased water solubility, enhanced metabolic stability, and reduced lipophilicity. In Vorinostat, such replacement results in a new bioactive analog of the drug. This study enhances the repertoire of available saturated bioisosteres of (hetero)aromatic rings for the use in drug discovery projects.

C-C Coupling through Nitrogen Deletion: Application to Library Synthesis.

Invited for the cover of this issue are Oleksandr Grygorenko and his Ukrainian colleagues at Enamine Ltd., Taras Shevchenko National University of Kyiv, National Academy of Sciences of Ukraine, and ChemSpace, as well as Mark Levin at the University of Chicago. The image depicts application of a nitrogen-deleting anomeric amide to parallel C(sp3 )-C(sp3 ) coupling. Read the full text of the article at 10.1002/chem.202203470.

C-C Coupling through Nitrogen Deletion: Application to Library Synthesis.

A protocol for parallel C(sp3 )-C(sp3 ) coupling of (hetero)aromatic aldehydes and (hetero)arylmethyl amines based on a reductive amination - "nitrogen deletion" reaction sequence has been developed. After preliminary validation experiments, an illustrative compound library of 25 members was prepared with 76 % synthetic efficiency. The estimated chemical space accessible by the proposed approach covers almost 600 000 representatives that are scarcely represented in current compound databases.

Creation of targeted compound libraries based on 3D shape recognition.

In the emerging field of drug discovery, rapid virtual screening methods become extremely valuable, especially when dealing with ultra-large databases of organic small bioactive molecules. In this work, we present a fast, computationally resource-efficient, and simple workflow for screening targeted compound libraries generated from ultra-large virtual chemical space. This workflow aims to find compounds with similar molecular 3D shapes with reference ones, and at the same time to expand chemical diversity and to identify new and potentially active scaffolds. This pipeline ensures the enrichment of the generated libraries with novel chemotypes. Also, it was shown that delicate tailoring of the physicochemical parameters of the search set ensures that all library compounds will possess desired property distributions. A visual inspection has shown that found structures bind to the receptor in the same way as the reference ones. Using our screening workflow, we have created a number of conventional protein-targeted libraries: the GPCRs Targeted Library (531 K compounds) and the Protein Kinases Targeted Library (113 K compounds). The described pipeline and scripts are freely accessible at: https://github.com/ChemSpace-LLC/usrcat_sim .

Virtual Screening in Search for a Chemical Probe for Angiotensin-Converting Enzyme 2 (ACE2).

We elaborate new models for ACE and ACE2 receptors with an excellent prediction power compared to previous models. We propose promising workflows for working with huge compound collections, thereby enabling us to discover optimized protocols for virtual screening management. The efficacy of elaborated roadmaps is demonstrated through the cost-effective molecular docking of 1.4 billion compounds. Savings of up to 10-fold in CPU time are demonstrated. These developments allowed us to evaluate ACE2/ACE selectivity in silico, which is a crucial checkpoint for developing chemical probes for ACE2.

Analysis of External and Internal Disorder to Understand Band-Like Transport in n-Type Organic Semiconductors.

Charge transport in organic semiconductors is notoriously extremely sensitive to the presence of disorder, both internal and external (i.e., related to interactions with the dielectric layer), especially for n-type materials. Internal dynamic disorder stems from large thermal fluctuations both in intermolecular transfer integrals and (molecular) site energies in weakly interacting van der Waals solids and sources transient localization of the charge carriers. The molecular vibrations that drive transient localization typically operate at low-frequency (

Rotator side chains trigger cooperative transition for shape and function memory effect in organic semiconductors.

Martensitic transition is a solid-state phase transition involving cooperative movement of atoms, mostly studied in metallurgy. The main characteristics are low transition barrier, ultrafast kinetics, and structural reversibility. They are rarely observed in molecular crystals, and hence the origin and mechanism are largely unexplored. Here we report the discovery of martensitic transition in single crystals of two different organic semiconductors. In situ microscopy, single-crystal X-ray diffraction, Raman and nuclear magnetic resonance spectroscopy, and molecular simulations combined indicate that the rotating bulky side chains trigger cooperative transition. Cooperativity enables shape memory effect in single crystals and function memory effect in thin film transistors. We establish a molecular design rule to trigger martensitic transition in organic semiconductors, showing promise for designing next-generation smart multifunctional materials.

Ab initio random structure searching of organic molecular solids: assessment and validation against experimental data.

This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated.

Donor-acceptor stacking arrangements in bulk and thin-film high-mobility conjugated polymers characterized using molecular modelling and MAS and surface-enhanced solid-state NMR spectroscopy.

Conjugated polymers show promising properties as cheap, sustainable and solution-processable semiconductors. A key challenge in the development of these materials is to determine the polymer chain structure, conformation and packing in both the bulk polymer and in thin films typically used in devices. However, many characterisation techniques are unable to provide atomic-level structural information owing to the presence of disorder. Here, we use molecular modelling, magic-angle spinning (MAS) and dynamic nuclear polarisation surface-enhanced NMR spectroscopy (DNP SENS) to characterise the polymer backbone group conformations and packing arrangement in the high-mobility donor-acceptor copolymer diketopyrrolo-pyrrole-dithienylthieno[3,2-b]thiophene (DPP-DTT). Using conventional 1H and 13C solid-state MAS NMR coupled with density functional theory calculations and molecular dynamics simulations, we find that the bulk polymer adopts a highly planar backbone conformation with a laterally-shifted donor-on-acceptor stacking arrangement. DNP SENS enables acquisition of 13C NMR data for polymer films, where sensitivity is limiting owing to small sample volumes. The DNP signal enhancement enables a two-dimensional 1H-13C HETCOR spectrum to be recorded for a drop-cast polymer film, and a 13C CPMAS NMR spectrum to be recorded for a spin-coated thin-film with a thickness of only 400 nm. The results show that the same planar backbone structure and intermolecular stacking arrangement is preserved in the films following solution processing and annealing, thereby rationalizing the favourable device properties of DPP-DTT, and providing a protocol for the study of other thin film materials.

Emulating Natural Product Conformation by Cooperative, Non-Covalent Fluorine Interactions.

Pervasive in Nature, the propane unit is an essential component of numerous bioactive molecules. These range from acyclic systems, such as the neurotransmitter γ-aminobutyric acid, through to the bicyclic nuclei of various chromanes and dihydrobenzofurans. In the latter case, cyclisation via cyclic ether formation ensures a highly pre-organised structure, whilst linear scaffolds display more dynamic conformational behaviour resulting from rotation about the two internal C(sp3 )-C(sp3 ) bonds. In this study, the replacement of -[CH2 ]- units by -[CHF]- centres is evaluated as a strategy to achieve acyclic conformational control by hindering these internal rotations. Reinforcing, non-covalent fluorine interactions are validated as powerful design features that result in programmable conformational behaviours: These are encoded by the relative configuration of each centre. By exploiting cooperative neighbouring stereoelectronic effects in a multi-vicinal fluoroalkane it is possible to emulate the overall conformation of the dihydrobenzofuran scaffold found in a variety of natural products with an acyclic mimic. This is described as a function of two bond vectors at the chain termini and validated by combined theoretical, crystallographic and spectroscopic analyses. In view of the favourable physicochemical properties associated with fluorine introduction, this approach to bioactive scaffold design may prove to be expansive.

Unraveling Unprecedented Charge Carrier Mobility through Structure Property Relationship of Four Isomers of Didodecyl[1]benzothieno[3,2-b][1]benzothiophene.

The structural and electronic properties of four isomers of didodecyl[1]-benzothieno[3,2-b][1]benzothiophene (C12-BTBT) have been investigated. Results show the strong impact of the molecular packing on charge carrier transport and electronic polarization properties. Field-induced time-resolved microwave conductivity measurements unravel an unprecedented high average interfacial mobility of 170 cm(2) V(-1) s(-1) for the 2,7-isomer, holding great promise for the field of organic electronics.

Ammonia as a case study for the spontaneous ionization of a simple hydrogen-bonded compound.

Modern ab initio calculations predict ionic and superionic states in highly compressed water and ammonia. The prediction apparently contradicts state-of-the-art experimentally established phase diagrams overwhelmingly dominated by molecular phases. Here we present experimental evidence that the threshold pressure of ~120 GPa induces in molecular ammonia the process of autoionization to yet experimentally unknown ionic compound--ammonium amide. Our supplementary theoretical simulations provide valuable insight into the mechanism of autoionization showing no hydrogen bond symmetrization along the transformation path, a remarkably small energy barrier between competing phases and the impact of structural rearrangement contribution on the overall conversion rate. This discovery is bridging theory and experiment thus opening new possibilities for studying molecular interactions in hydrogen-bonded systems. Experimental knowledge on this novel ionic phase of ammonia also provides strong motivation for reconsideration of the theory of molecular ice layers formation and dynamics in giant gas planets.

Coexistence of helical morphologies in columnar stacks of star-shaped discotic hydrazones.

Discotic hydrazone molecules are of particular interest as they form discotic phases where the discs are rigidified by intramolecular hydrogen bonds. Here, we investigate the thermotropic behavior and solid-state organizations of three discotic hydrazone derivatives with dendritic groups attached to their outer peripheries, containing six, eight, and ten carbons of linear alkoxy chains. On the basis of two-dimensional wide angle X-ray scattering (2DWAXS), the elevated temperature liquid crystalline (LC) phases were assigned to a hexagonal columnar (Colh) organization with nontilted hydrazone discs for all three compounds. With WAXS, advanced solid-state nuclear magnetic resonance (SSNMR) techniques, and ab initio computations, the compounds with six and ten carbons of achiral alkoxy side chains were further subjected to studies at 25 °C, revealing complex crystalline phases with rigid columns and flexible side chains. This combined approach led to models of coexisting helical columnar stacking morphologies for both systems with two different tilt/pitch angles between successive hydrazone molecules. The differences in tilt/pitch angles between the two compounds illustrate that the columns with short alkoxy chains (six carbons) are more influenced by the presence of other stacks in their vicinity, while those with long side chains are less tilted due to a larger alkoxy (ten carbons) buffer zone. The formation of different packing morphologies in the crystalline phase of a columnar LC has rarely been reported so far, which suggests the possibility of complex stacking structures of similar organic LC systems, utilizing small molecules as potential materials for applications in organic electronics.

Exploiting the Synergy of Powder X-ray Diffraction and Solid-State NMR Spectroscopy in Structure Determination of Organic Molecular Solids.

We report a strategy for structure determination of organic materials in which complete solid-state nuclear magnetic resonance (NMR) spectral data is utilized within the context of structure determination from powder X-ray diffraction (XRD) data. Following determination of the crystal structure from powder XRD data, first-principles density functional theory-based techniques within the GIPAW approach are exploited to calculate the solid-state NMR data for the structure, followed by careful scrutiny of the agreement with experimental solid-state NMR data. The successful application of this approach is demonstrated by structure determination of the 1:1 cocrystal of indomethacin and nicotinamide. The 1H and 13C chemical shifts calculated for the crystal structure determined from the powder XRD data are in excellent agreement with those measured experimentally, notably including the two-dimensional correlation of 1H and 13C chemical shifts for directly bonded 13C-1H moieties. The key feature of this combined approach is that the quality of the structure determined is assessed both against experimental powder XRD data and against experimental solid-state NMR data, thus providing a very robust validation of the veracity of the structure.

Solid-State Organization and Ambipolar Field-Effect Transistors of Benzothiadiazole-Cyclopentadithiophene Copolymer with Long Branched Alkyl Side Chains.

The solid-state organization of a benzothiadiazole-cyclopentadithiophene copolymer with long, branched decyl-tetradecyl side chains (CDT-BTZ-C14,10) is investigated. The C14,10 substituents are sterically demanding and increase the π-stacking distance to 0.40 nm from 0.37 nm for the same polymer with linear hexadecyls (C16). Despite the bulkiness, the C14,10 side chains tend to crystallize, leading to a small chain-to-chain distance between lamellae stacks and to a crystal-like microstructure in the thin film. Interestingly, field-effect transistors based on solution processed layers of CDT-BTZ-C14,10 show ambipolar behavior in contrast to CDT-BTZ-C16 with linear side chains, for which hole transport was previously observed. Due to the increased π-stacking distance, the mobilities are only 6 × 10-4 cm²/Vs for electrons and 6 × 10-5 cm²/Vs for holes, while CDT-BTZ-C16 leads to values up to 5.5 cm²/Vs. The ambipolarity is attributed to a lateral shift between stacked backbones provoked by the bulky C14,10 side chains. This reorganization is supposed to change the transfer integrals between the C16 and C14,10 substituted polymers. This work shows that the electronic behavior in devices of one single conjugated polymer (in this case CDT-BTZ) can be controlled by the right choice of the substituents to place the backbones in the desired packing.

Probing the relation between charge transport and supramolecular organization down to ångström resolution in a benzothiadiazole-cyclopentadithiophene copolymer.

Molecular modeling shows that longitudinal displacement of the backbones by a couple of ångströms has a profound impact on the electronic coupling mediating charge transport in a conjugated copolymer. These changes can be probed by monitoring the calculated X-ray scattering patterns and NMR chemical shifts as a function of sliding of the polymer chains and comparing them to experiment.

A strategy for revealing the packing in semicrystalline π-conjugated polymers: crystal structure of bulk poly-3-hexyl-thiophene (P3HT).

To tilt or not to tilt: The crystal structure for bulk P3HT (phase I) was determined by "multi-technique crystallography", which combines X-ray diffraction, solid-state NMR spectroscopy, and DFT calculations. The results showed that this semiconducting polymer crystallizes in the monoclinic space group P2(1)/c with nontilted π-stacks at a distance of 3.9 Š(see picture).

Use of X-ray diffraction, molecular simulations, and spectroscopy to determine the molecular packing in a polymer-fullerene bimolecular crystal.

The molecular packing in a polymer: fullerene bimolecular crystal is determined using X-ray diffraction (XRD), molecular mechanics (MM) and molecular dynamics (MD) simulations, 2D solid-state NMR spectroscopy, and IR absorption spectroscopy. The conformation of the electron-donating polymer is significantly disrupted by the incorporation of the electron-accepting fullerene molecules, which introduce twists and bends along the polymer backbone and 1D electron-conducting fullerene channels.