Search

The synthesis of aryl amines from 3-alkynyl-2-pyrones and various amines is described. Mechanistically, the aryl amines are proposed to arise from the 3-alkynyl-2-pyrone substrates through their selective opening in a 1,6-fashion by secondary amines followed by decarboxylation and an unexpected rearrangement. The proposed mechanism is supported by quantum chemical transition-state calculations, which are consistent with the regiochemical outcome. The scope of this transformation spans a variety of 3-alkynyl-2-pyrones and a range of secondary amines. The influence of the secondary amine coupling partners on reaction efficiency was elucidated through data-driven modeling as well as scope exploration. These latter studies revealed that the steric bulk of the secondary amine coupling partner under the reaction conditions serves as a strong indicator of overall reaction efficiency., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Designing solvent systems is key to achieving the facile synthesis and separation of desired products from chemical processes, so many machine learning models have been developed to predict solubilities. However, breakthroughs are needed to address deficiencies in the model's predictive accuracy and generalizability; this can be addressed by expanding and integrating experimental and computational solubility databases. To maximize predictive accuracy, these two databases should not be trained separately, and they should not be simply combined without reconciling the discrepancies from different magnitudes of errors and uncertainties. Here, we introduce self-evolving solubility databases and graph neural networks developed through semi-supervised self-training approaches. Solubilities from quantum-mechanical calculations are referred to during semi-supervised learning, but they are not directly added to the experimental database. Dataset augmentation is performed from 11 637 experimental solubilities to >900 000 data points in the integrated database, while correcting for the discrepancies between experiment and computation. Our model was successfully applied to study solvent selection in organic reactions and separation processes. The accuracy (mean absolute error around 0.2 kcal mol -1 for the test set) is quantitatively useful in exploring Linear Free Energy Relationships between reaction rates and solvation free energies for 11 organic reactions. Our model also accurately predicted the partition coefficients of lignin-derived monomers and drug-like molecules. While there is room for expanding solubility predictions to transition states, radicals, charged species, and organometallic complexes, this approach will be attractive to predictive chemistry areas where experimental, computational, and other heterogeneous data should be combined., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Long-lived organic radicals are promising candidates for the development of high-performance energy solutions such as organic redox batteries, transistors, and light-emitting diodes. However, "stable" organic radicals that remain unreactive for an extended time and that can be stored and handled under ambient conditions are rare. A necessary but not sufficient condition for organic radical stability is the presence of thermodynamic stabilization, such as conjugation with an adjacent π-bond or lone-pair, or hyperconjugation with a σ-bond. However, thermodynamic factors alone do not result in radicals with extended lifetimes: many resonance-stabilized radicals are transient species that exist for less than a millisecond. Kinetic stabilization is also necessary for persistence, such as steric effects that inhibit radical dimerization or reaction with solvent molecules. We describe a quantitative approach to map organic radical stability, using molecular descriptors intended to capture thermodynamic and kinetic considerations. The comparison of an extensive dataset of quantum chemical calculations of organic radicals with experimentally-known stable radical species reveals a region of this feature space where long-lived radicals are located. These descriptors, based upon maximum spin density and buried volume, are combined into a single metric, the radical stability score, that outperforms thermodynamic scales based on bond dissociation enthalpies in identifying remarkably long-lived radicals. This provides an objective and accessible metric for use in future molecular design and optimization campaigns. We demonstrate this approach in identifying Pareto-optimal candidates for stable organic radicals., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Nuclear magnetic resonance (NMR) is one of the primary techniques used to elucidate the chemical structure, bonding, stereochemistry, and conformation of organic compounds. The distinct chemical shifts in an NMR spectrum depend upon each atom's local chemical environment and are influenced by both through-bond and through-space interactions with other atoms and functional groups. The in silico prediction of NMR chemical shifts using quantum mechanical (QM) calculations is now commonplace in aiding organic structural assignment since spectra can be computed for several candidate structures and then compared with experimental values to find the best possible match. However, the computational demands of calculating multiple structural- and stereo-isomers, each of which may typically exist as an ensemble of rapidly-interconverting conformations, are expensive. Additionally, the QM predictions themselves may lack sufficient accuracy to identify a correct structure. In this work, we address both of these shortcomings by developing a rapid machine learning (ML) protocol to predict 1 H and 13 C chemical shifts through an efficient graph neural network (GNN) using 3D structures as input. Transfer learning with experimental data is used to improve the final prediction accuracy of a model trained using QM calculations. When tested on the CHESHIRE dataset, the proposed model predicts observed 13 C chemical shifts with comparable accuracy to the best-performing DFT functionals (1.5 ppm) in around 1/6000 of the CPU time. An automated prediction webserver and graphical interface are accessible online at http://nova.chem.colostate.edu/cascade/. We further demonstrate the model in three applications: first, we use the model to decide the correct organic structure from candidates through experimental spectra, including complex stereoisomers; second, we automatically detect and revise incorrect chemical shift assignments in a popular NMR database, the NMRShiftDB; and third, we use NMR chemical shifts as descriptors for determination of the sites of electrophilic aromatic substitution., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Radical cations generated from the oxidation of C[double bond, length as m-dash]C π-bonds are synthetically useful reactive intermediates for C-C and C-X bond formation. Radical cation formation, induced by sub-stoichiometric amounts of external oxidant, are important intermediates in the Woodward-Hoffmann thermally disallowed [2 + 2] cycloaddition of electron-rich alkenes. Using density functional theory (DFT), we report the detailed mechanisms underlying the intermolecular heterodimerisation of anethole and β-methylstyrene to give unsymmetrical, tetra-substituted cyclobutanes. Reactions between trans -alkenes favour the all-trans adduct, resulting from a kinetic preference for anti -addition reinforced by reversibility at ambient temperatures since this is also the thermodynamic product; on the other hand, reactions between a trans -alkene and a cis -alkene favour syn -addition, while exocyclic rotation in the acyclic radical cation intermediate is also possible since C-C forming barriers are higher. Computations are consistent with the experimental observation that hexafluoroisopropanol ( HFIP ) is a better solvent than acetonitrile, in part due to its ability to stabilise the reduced form of the hypervalent iodine initiator by hydrogen bonding, but also through the stabilisation of radical cationic intermediates along the reaction coordinate., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Atropisomeric biaryl systems are privileged architectures used in asymmetric synthesis and pharmaceutical structures. We report that by simply removing a single-electron, the resistance of biaryls towards racemization is reduced dramatically. Even though the steric properties are unaltered, biaryl oxidation changes atropisomerization into a two step mechanism with considerably smaller activation barriers than closed-shell biaryls. The effect is general for a series of biaryls and helicenes studied and results from the dependence of frontier molecular orbital energies on biaryl conformation.

The first enantioselective sulfa-Michael addition of alkyl thiols to alkenyl benzimidazoles, enabled by a bifunctional iminophosphorane (BIMP) organocatalyst, is described. The iminophosphorane moiety of the catalyst provides the required basicity to deprotonate the thiol nucleophile while the chiral scaffold and H-bond donor control facial selectivity. The reaction is broad in scope with respect to the thiol and benzimidazole reaction partners with the reaction proceeding in up to 98% yield and 96 : 4 er.

A new asymmetric conjugate addition method was developed for β-substituted cyclopentenones to form quaternary centres using alkylzirconocene nucleophiles giving up to 97% yield and 92% ee. Key to the reaction's success was the design of suitable phosphoramidite ligands which was aided by a linear quantitative structure-selectivity relationship (QSSR). QSSR models were created from the ligand screening data (a total of 36 ligands) which revealed important electronic and steric requirements and led to the synthesis of more enantioselective ligands. DFT calculations of competing transition structures enable the interpretation of the electronic and steric terms present in the QSSR models.

We have studied the cation-π interactions of neutral aromatic ligands with the cationic amino acid residues arginine, histidine and lysine using ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of all available Protein Data Bank (PDB) X-ray structures. Quantum chemical potential energy surfaces (PES) for these interactions were obtained at the DLPNO-CCSD(T) level of theory and compared against the empirical distribution of 2012 unique protein-ligand cation-π interactions found in X-ray crystal structures. We created a workflow to extract these structures from the PDB, filtering by interaction type and residue p K a . The gas phase cation-π interaction of lysine is the strongest by more than 10 kcal mol -1 , but the empirical distribution of 582 X-ray structures lies away from the minimum on the interaction PES. In contrast, 1381 structures involving arginine match the underlying calculated PES with good agreement. SAPT analysis revealed that underlying differences in the balance of electrostatic and dispersion contributions are responsible for this behavior in the context of the protein environment. The lysine-arene interaction, dominated by electrostatics, is greatly weakened by a surrounding dielectric medium and causes it to become essentially negligible in strength and without a well-defined equilibrium separation. The arginine-arene interaction involves a near equal mix of dispersion and electrostatic attraction, which is weakened to a much smaller degree by the surrounding medium. Our results account for the paucity of cation-π interactions involving lysine, even though this is a more common residue than arginine. Aromatic ligands are most likely to interact with cationic arginine residues as this interaction is stronger than for lysine in higher polarity surroundings.

Sulfones feature prominently in biologically active molecules and are key functional groups for organic synthesis. We report a mild, photoredox-catalyzed reaction for sulfonylation of aniline derivatives with sulfinate salts, and demonstrate the utility of the method by the late-stage functionalization of drugs. Key features of the method are the straightforward generation of sulfonyl radicals from bench-stable sulfinate salts and the use of simple aniline derivatives as convenient readily available coupling partners.

The diastereoselective coupling of O -substituted arylhydroxamates and cyclopropenes mediated by Rh(iii) catalysis was successfully developed. Through ligand development, the diastereoselectivity of this reaction was improved using a heptamethylindenyl (Ind*) ligand, which has been rationalized using quantum chemical calculations. In addition, the nature of the O -substituted ester of benzhydroxamic acid proved important for high diastereoselectivity. This transformation tolerates a variety of benzamides and cyclopropenes that furnish cyclopropa[ c ]dihydroisoquinolones with high diastereocontrol, which could then be easily transformed into synthetically useful building blocks for pharmaceuticals and bio-active molecules.

We report a 12-step catalytic enantioselective formal synthesis of malhamensilipin A ( 3 ) and diastereoisomeric analogues from ( E )-2-undecenal. The convergent synthesis relied upon iterative epoxidation and phosphorus(v)-mediated deoxydichlorination reactions as well a titanium-mediated epoxide-opening to construct the C11-C16 stereohexad. The latter transformation occurred with very high levels of stereoretention regardless of the C13 configuration of the parent epoxide, implicating anchimeric assistance of either the γ- or δ-chlorine atoms, and the formation of chloretanium or chlorolanium ions, respectively. A computational analysis of the chloronium ion intermediates provided support for the involvement of chlorolanium ions, whereas the potential chloretanium ions were found to be less likely intermediates on the basis of their greater carbocationic character.

Understanding the regioselectivity of C-H activation in the absence of directing groups is an important step towards the design of site-selective C-H functionalizations. The Pd(ii)-catalyzed direct arylation of chromones and enaminones provides an intriguing example where a simple substitution leads to a divergence in substrate-controlled site-selectivity. We describe computational and experimental studies which reveal this results from a switch in mechanism and therefore the selectivity-determining step. We present computational results and experimentally measured kinetic isotope effects and labelling studies consistent with this proposal. The C-H activation of these substrates proceeds via a CMD mechanism, which favors more electron rich positions and therefore displays a pronounced kinetic selectivity for the C3-position. However, C2-selective carbopalladation is also a competitive pathway for chromones so that the overall regiochemical outcome depends on which substrate undergoes activation first. Our studies provide insight into the site-selectivity based on the favorability of two competing CMD and carbopalladation processes of the substrates undergoing coupling. This model can be utilized to predict the regioselectivity of coumarins which are proficient substrates for carbopalladation. Furthermore, our model is able to account for the opposite selectivities observed for enaminone and chromone, and explains how a less reactive coupling partner leads to a switch in selectivity.

The nucleophilic reactivity of fluoride ion is altered in the presence of hydrogen-bond donors, including alcohols. Relatively little is known about the coordination involved; to rectify this, the X-ray structures of fourteen novel fluoride-alcohol complexes with tetrabutylammonium as the counterion have been determined. The coordination number varies from two to four depending on the steric bulk of the alcohol and is closely linked to trends in reactivity. This diversity in coordination stoichiometry is unprecedented but significant, as it implies differences in the ability of the fluoride-alcohol complexes to dissociate in solution with release of a more active and/or selective fluoride source.