About Underappreciated Yet Active Conformations of Thiourea Organocatalysts

Formation of chalcogen-bonding interactions and their role in the trans-trans conformation of thiourea

The chalcogen-bonding interactions formed at both sides of the thiocarbonyl group in thiourea were investigated. In particular, the role of these chalcogen-bonding interactions in the trans-trans conformation of thiourea was evaluated via single-crystal X-ray diffraction analysis and DFT calculations. The obtained results indicate that the Se⋯S⋯Se dual chalcogen-bonding interactions play a stronger role in controlling the planar structure than the S⋯S⋯S interactions.

Metallocenium incorporated charge-enhanced thiourea catalysts

Several metallocenium derivatives were prepared including N-ferrocenium-N'-phenylthiourea tetrakis-[3,5-bis(trifluoromethyl)phenyl]borate (BArF4, 4a), N-methylferrocenium-N'-phenylthiourea BArF4 (4b) and N-cobaltocenium-N'-phenylthiourea BArF4 (5). These compounds are competent catalysts for the Friedel-Crafts alkylation of indoles with trans-β-nitrostyrenes, and are much more active than their reduced non-charged forms. The iron derivatives are less stable than the cobalt analog and were generated and used in situ whereas the cobalt-containing thiourea was isolated and characterized by X-ray crystallography. All three of these metallocenium salts have Lewis and Brønsted acidic sites which were exploited in tandem to afford charge-activated catalysts.

Elucidating Conformation and Hydrogen-Bonding Motifs of Reactive Thiourea Intermediates

Substituted diphenylthioureas (DPTUs) are efficient hydrogen-bonding organo-catalysts, and substitution of DPTUs has been shown to greatly affect catalytic activity. Yet, both the conformation of DPTUs in solution and the conformation and hydrogen-bonded motifs within catalytically active intermediates, pertinent to their mode of activation, have remained elusive. By combining linear and ultrafast vibrational spectroscopy with spectroscopic simulations and calculations, we show that different conformational states of thioureas give rise to distinctively different N-H stretching bands in the infrared spectra. In the absence of hydrogen-bond-accepting substrates, we show that vibrational structure and dynamics are highly sensitive to the substitution of DPTUs with CF3 groups and to the interaction with the solvent environment, allowing for disentangling the different conformational states. In contrast to bare diphenylthiourea (0CF-DPTU), we find the catalytically superior CF3-substituted DPTU (4CF-DPTU) to favor the trans-trans conformation in solution, allowing for donating two hydrogen bonds to the reactive substrate. In the presence of a prototypical substrate, DPTUs in trans-trans conformation hydrogen bond to the substrate's C=O group, as evidenced by a red-shift of the N-H vibration. Yet, our time-resolved infrared experiments indicate that only one N-H group forms a strong hydrogen bond to the carbonyl moiety, while thiourea's second N-H group only weakly interacts with the substrate. Our data indicate that hydrogen-bond exchange between these N-H groups occurs on the timescale of a few picoseconds for 0CF-DPTU and is significantly accelerated upon CF3 substitution. Our results highlight the subtle interplay between conformational equilibria, bonding states, and bonding lifetimes in reactive intermediates in thiourea catalysis, which help rationalize their catalytic activity.

Bis-Thiourea Chiral Sensor for the NMR Enantiodiscrimination of N-Acetyl and N-Trifluoroacetyl Amino Acid Derivatives

A C2-symmetrical bis-thiourea chiral solvating agent (CSA), TFTDA, for NMR spectroscopy has been obtained by reacting (1R,2R)-1,2-bis(2-hydroxyphenyl)ethylenediamine and 3,5-bis(trifluoromethyl)phenyl isothiocyanate. TFTDA shows remarkable propensity to enantiodiscriminate N-trifluoroacetyl (N-TFA) and N-acetyl (N-Ac) derivatives of amino acids with free carboxyl functions, with the co-presence of 1,4-diazabicyclo[2.2.2]octane (DABCO) as the third achiral additive, which is needed for substrate solubilization. TFTDA shows enhanced enantiodiscriminating efficiency in comparison with the corresponding monomeric counterpart, TFTMA, pointing out cooperativity between its two symmetrical entities. A wide range of amino acid derivatives have been efficiently enantiodiscriminated in CDCl₃, with high enantioresolution quotients, which guarantee high quality in applications devoted to the quantification of enantiomers. High enantiodiscriminating efficiency is maintained also in diluted 5 mM conditions or in the presence of sub-stoichiometric amounts of CSA (0.3 equiv). The role of phenolic hydroxyls in the DABCO-mediated interaction mechanism between TFTDA and the two enantiomeric substrates has been pointed out by means of diffusion-ordered spectroscopy (DOSY) and rotating frame Overhauser effect spectroscopy (ROESY) experiments. A conformational model for both the CSA and its diastereomeric solvates formed with the two enantiomers of N-acetyl leucine has also been conceived on the basis of ROE data in order to give a chiral discrimination rationale.

Anion-binding of a chiral tris(2-aminoethyl)amine-based tripodal thiourea: A spectroscopic and computational study

Thioureas are well-known structural motifs in supramolecular anion recognition. Their conformational preferences are typically characterized by detailed NMR spectroscopy and crystallography, which are often complemented with computational results from geometry...

Non-Covalent Interactions in Enantioselective Organocatalysis: Theoretical and Mechanistic Studies of Reactions Mediated by Dual H-Bond Donors, Bifunctional Squaramides, Thioureas and Related Catalysts

Chiral bifunctional dual H-bond donor catalysts have become one of the pillars of organocatalysis. They include squaramide, thiosquaramide, thiourea, urea, and even selenourea-based catalysts combined with chiral amines, cinchona alkaloids, sulfides, phosphines and more. They can promote several types of reactions affording products in very high yields and excellent stereoselectivities in many cases: conjugate additions, cycloadditions, the aldol and Henry reactions, the Morita–Baylis–Hilman reaction, even cascade reactions, among others. The desire to understand mechanisms and the quest for the origins of stereoselectivity, in attempts to find guidelines for developing more efficient catalysts for new transformations, has promoted many mechanistic and theoretical studies. In this review, we survey the literature published in this area since 2015.

Accuracy of DLPNO-CCSD(T): Effect of Basis Set and System Size

The DLPNO-CCSD(T) method is designed to study large molecular systems at significantly reduced cost relative to its canonical counterpart. However, the error in this approach is also size-extensive and relies on cancellation of errors for the calculation of relative energies. This work provides a direct comparison of canonical CCSD(T) and TightPNO DLPNO-CCSD(T) calculations of reaction energies and barriers of a broad range of chemical reactions. The dataset includes acidities, anion binding affinities, enolization, Diels-Alder, nucleophilic substitution, and atom transfer reactions and complements existing theoretical datasets in terms of system size as well as new reaction types (e.g., anion binding affinities and chlorine atom transfer reactions). The performance of DLPNO-CCSD(T) was further examined with respect to systematic variation of basis set and system size and amounts of nonbonded interaction present in the system. The errors in the DLPNO-CCSD(T) were found to be relatively insensitive to the choice of basis set for small systems but increase monotonically with system size. Additionally, calculations of barriers appear to be more challenging than reaction energies with errors exceeding 5 kJ mol^-1 for many Diels-Alder reactions. Further tests on three realistic organic reactions reveal the impact of the DLPNO approximation in calculating absolute and relative barriers that are important for predictions such as stereoselectivity.

Anion Binding Affinity: Acidity versus Conformational Effects

High-level quantum chemical calculations were used to elucidate the gas- and solution-phase conformational equilibria for a series of symmetrically substituted (thio)ureas, (thio)squaramides, and croconamides. Gas-phase calculations predict that the thermodynamic conformer of many of these anion receptors is not the dual-hydrogen-bond-facilitating anti-anti conformer as is commonly assumed. For N,N'-diaryl thiosquaramides and croconamides, the syn-syn conformer is typically the predominant conformer. Solution-phase calculations show that the anti-anti conformer is increasingly stabilized as the polarity of the solvent increases. However, the syn-syn conformer remains the lowest energy conformation for croconamides. These predictions are used to explain the acidity versus chloride binding affinity correlations recently reported for some of these compounds. The chloride binding constants for thioureas and croconamides are significantly lower than expected on the basis of their pK_a values, and this may be due in part to the need for these receptors to reorganize into the anti-anti conformer. Experimental NMR nuclear Overhauser effect (NOE) measurements of an asymmetrically substituted squaramide and its thio analogue are consistent with the syn-syn conformation being predominant at 298 K. The conformational equilibria should therefore be an important consideration for the design and development of future anion receptors and organocatalysts.

Theoretical free energy profile and benchmarking of functionals for amino-thiourea organocatalyzed nitro-Michael addition reaction

Amino-thiourea organocatalysis is an important catalytic process for enantioselective conjugate addition reactions. The interaction of the reactants with the catalyst has a substantial effect of dispersion forces and is a challenge for a reliable description when applying density functional theory. In this report, the classical addition of acetylacetone to β-nitro-styrene catalyzed by Takemoto's catalyst in toluene was studied using the PBE functional for geometry optimization and the DLPNO-CCSD(T) benchmark method for single point energy. The complete free energy profile calculated for the reaction was able to explain all experimental observations, including the fact that the carbon-carbon bond formation step is rate-determining. The overall barrier was calculated to be 22.8 kcal mol^-1 (experimental value approximately 20 kcal mol^-1), and the enantiomeric excess was calculated to be 88% (experimental value in the range of 84 to 92%). Some functionals were tested for single point energy. The hybrid B3LYP presented a high mean absolute deviation (MAD) from the DLPNO-CCSD(T) benchmark method by approximately 20 kcal mol^-1. The inclusion of empirical dispersion correction in the B3LYP method decreased the MAD to 6 kcal mol^-1. Even the double-hybrid mPW2-PLYP and B2GP-PLYP methods had MAD values of approximately 5 kcal mol^-1. The inclusion of the dispersion correction decreased the MAD to 3.6 kcal mol^-1. M06-2X and ωB97X-D3 were the most accurate among the tested functionals, with MADs of 2.5 kcal mol^-1 and 1.8 kcal mol^-1, respectively. Additivity approximation of the correlation energy was also tested and presented a MAD of only 0.6 kcal mol^-1.

Mechanistic Insight toward Understanding the Role of Charge in Thiourea Organocatalysis

Pyranylation and glycosylation are pivotal for accessing a myriad of natural products, pharmaceuticals, and drug candidates. Catalytic approaches for enabling these transformations are of utmost importance and integral to advancing this area of synthesis. In exploring this chemical space, a combined experimental and computational mechanistic study of pyranylation and 2-deoxygalactosylation catalyzed by a cationic thiourea organocatalyst is reported. To this end, a thiourea-cyclopropenium organocatalyst was employed as a model system in combination with an arsenal of mechanistic techniques, including ¹³C kinetic isotope effect experiments, deuterated labeling studies, variable-temperature ¹H NMR spectroscopy, and density functional theory calculations. From these studies, two distinct reaction pathways were identified for this transformation corresponding to either dual hydrogen bond (H-bond) activation or Brønsted acid catalysis. The former involving thiourea orchestrated bifurcated hydrogen bonding proceeded in an asynchronous concerted fashion. In contrast, the latter stepwise mechanism involving Brønsted acid catalysis hinged upon the formation of an oxocarbenium intermediate accompanied by subsequent alcohol addition.

Unravelling the mechanism of the organocatalyzed aminolytic kinetic resolution of α-nitroepoxides: a theoretical study

The mechanism of the aminolytic kinetic resolution of 2-methyl-2-nitro-3-phenyl oxirane with aniline in the presence of Takemoto's catalyst is disclosed by way of DFT calculations.

Dynamic Refolding of Ion-Pair Catalysts in Response to Different Anions

Four distinct folding patterns are identified in two foldamer-type urea-thiourea catalysts bearing a basic dimethylamino unit by a combination of X-ray crystallography, solution NMR studies, and computational studies (DFT). These patterns are characterized by different intramolecular hydrogen bonding schemes that arise largely from different thiourea conformers. The free base forms of the catalysts are characterized by folds where the intramolecular hydrogen bonds between the urea and the thiourea units remain intact. In contrast, the catalytically relevant salt forms of the catalyst, where the catalyst forms an ion pair with the substrate or substrate analogues, appear in two entirely different folding patterns. With larger anions that mimic the dialkyl malonate substrates, the catalysts maintain their native fold both in the solid state and in solution, but with smaller halide anions (fluoride, chloride, and bromide), the catalysts fold around the halide anion (anion receptor fold), and the intramolecular hydrogen bonds are disrupted. Titration of catalyst hexafluoroacetylacetonate salt with tetra-n-butylammonium chloride results in dynamic refolding of the catalyst from the native fold to the anion receptor fold.

Dynamic optical activity induction in the N-alkyl-N'-trityl ureas and thioureas

Considered to be rigid, the urea and thiourea functionalities, often used in material chemistry and in asymmetric organocatalysis, are able to transmit information regarding 3D structure from a permanently chiral inducer part to a dynamically chiral (reporter) part of the molecule. Despite a considerable distance between the inducer and the reporter parts of the molecule, the chirality transfer phenomenon has been demonstrated for a series of secondary N-alkyl-N'-trityl ureas and thioureas. The induction of helicity in a stereodynamic trityl propeller is revealed by rising non-zero Cotton effects in the area of trityl absorption. The information regarding the 3D structure of the inducer is transferred to the reporter part of the system through a set of weak but complementary electrostatic interactions. The presence of two supramolecular motifs in the same molecule, characterized by opposite properties, significantly affected the molecular solid state structure of the thioureas and their abilities to assemble. In the crystalline phase, the model, a chiral N-tert-butyl-thiourea derivative that retains the extended Z,Z conformation of the linker, is prone to form a supramolecular network typical of secondary ureas and thioureas. In contrast, the presence of the hydrophobic trityl group suppresses the thioamide NHS[double bond, length as m-dash]C hydrogen bonds. Therefore, trityl acts as a supramolecular protecting group for thioamide functionality, hampering the formation of hydrogen bonded networks in the solid state.

Transition state analysis of an enantioselective Michael addition by a bifunctional thiourea organocatalyst

The mechanism of the enantioselective Michael addition of diethyl malonate to trans-β-nitrostyrene catalyzed by a tertiary amine thiourea organocatalyst is explored using experimental 13C kinetic isotope effects and density functional theory calculations. Large primary 13C KIEs on the bond-forming carbon atoms of both reactants suggest that carbon-carbon bond formation is the rate-determining step in the catalytic cycle. This work resolves conflicting mechanistic pictures that have emerged from prior experimental and computational studies.

Stereoselective organocatalyzed glycosylations - thiouracil, thioureas and monothiophthalimide act as Brønsted acid catalysts at low loadings

Thiouracil catalyzes stereoselective glycosylations with galactals in loadings as low as 0.1 mol%. It is proposed that in these glycosylations thiouracil, monothiophthalimide, and the previously reported catalyst, Schreiner's thiourea, do not operate via a double H-bonding mechanism but rather by Brønsted acid/base catalysis. In addition to the synthesis of 2-deoxyglycosides and glycoconjugates, we report the first organocatalytic synthesis of 1,1'-linked trehalose-type sugars.

Why electrostatically enhanced thiourea is better than Schreiner's thiourea in both catalytic activity and regioselectivity?

Electrostatically enhanced thiourea is more active and efficient than Schreiner's thiourea in the ring-opening aminolysis of styrene oxide with aniline, and the underlying reasons were explored by DFT calculations.

Cyclopropenium Enhanced Thiourea Catalysis

An integral part of modern organocatalysis is the development and application of thiourea catalysts. Here, as part of our program aimed at developing cyclopropenium catalysts, the synthesis of a thiourea-cyclopropenium organocatalyst with both cationic hydrogen-bond donor and electrostatic character is reported. The utility of the this thiourea organocatalyst is showcased in pyranylation reactions employing phenols, primary, secondary, and tertiary alcohols under operationally simple and mild reaction conditions for a broad substrate scope. The addition of benzoic acid as a co-catalyst facilitating cooperative Brønsted acid catalysis was found to be valuable for reactions involving phenols and higher substituted alcohols. Mechanistic investigations, including kinetic and ¹H NMR binding studies in conjunction with density function theory calculations, are described that collectively support a Brønsted acid mode of catalysis.

Unraveling the Origin of Solvent Induced Enantioselectivity in the Henry Reaction with Cinchona Thiourea as Catalyst

In this work, we report an energy decomposition and electronic structure analysis using DFT calculations for the C-C coupling step in the Henry reaction with cinchona thiourea as catalyst and DMF solvent to unravel the origin of enantioselectivity. We found that the conformation of flexible thiourea moiety is affected by the solvent, and in the preferred conformation of thiourea in strong Lewis basic DMF solvent, the N-H sites are in the opposite direction, i.e., in trans conformation. Hence, the thiourea moiety acts via single hydrogen bonding with substrates. The conformation of the substrates with respect to the forming C-C bond plays critical role to increase orbital interaction between two substrates and enhances hydrogen bond strength between substrates and catalyst, which in turn stabilizes the positive charge developing on the catalyst at the transition state for one of the enantiomers ( S). Thus, the enantioselectivity has electronic structure origin. The stronger H-bond formation in the S enantiomer has been confirmed by the calculated IR spectra and is in agreement with thus far experimental and computational results.

Quantification of Catalytic Activity for Electrostatically Enhanced Thioureas via Reaction Kinetics and UV-vis Spectroscopic Measurement

Charged thiourea derivatives containing one and two methylated or octylated pyridinium ion centers and a tetraarylborate or triflate counteranion are reported. These novel catalysts are much more active in the Friedel-Crafts reactions of trans-β-nitroalkenes with N-methylindoles than the privileged N, N'-bis(3,5-bis(trifluoromethyl)phenyl)thiourea (i.e., Schreiner's thiourea) by up to 2-3 orders of magnitude. A previously reported UV-vis spectroscopic method by Kozlowski et al. was exploited to rationalize their reactivity order along with noncharged analogues. These results offer a new design strategy for organocatalysts by introducing positively charged centers without adding additional N-H, O-H, or S-H hydrogen bond donor sites.