Anti-Markovnikov Intermolecular Hydroamination of Alkenes and Alkynes: A Mechanistic View.

Hydroamination, the addition of an N-H bond across a C-C multiple bond, is a reaction with a great synthetic potential. Important advances have been made in the last decades concerning catalysis of these reactions. However, controlling the regioselectivity in the amine addition toward the formation of anti-Markovnikov products (addition to the less substituted carbon) still remains a challenge, particularly in intermolecular hydroaminations of alkenes and alkynes. The goal of this review is to collect the systems in which intermolecular hydroamination of terminal alkynes and alkenes with anti-Markovnikov regioselectivity has been achieved. The focus will be placed on the mechanistic aspects of such reactions, to discern the step at which regioselectivity is decided and to unravel the factors that favor the anti-Markovnikov regioselectivity. In addition to the processes entailing direct addition of the amine to the C-C multiple bond, alternative pathways, involving several reactions to accomplish anti-Markovnikov regioselectivity (formal hydroamination processes), will also be discussed in this review. The catalysts gathered embrace most of the metal groups of the Periodic Table. Finally, a section discussing radical-mediated and metal-free approaches, as well as heterogeneous catalyzed processes, is also included.

Chemoselective Ru-Catalyzed Oxidative Lactamization vs Hydroamination of Alkynylamines: Insights from Experimental and Density Functional Theory Studies.

The Ru-catalyzed intramolecular oxidative amidation (lactamization) of aromatic alkynylamines with 4-picoline N-oxide as an external oxidant has been developed. This chemoselective process is very efficient to achieve medium-sized ε- and ζ-lactams (seven- and eight-membered rings) but not for the formation of common δ-lactams (six-membered rings). DFT studies unveiled the capital role of the chain length between the amine and the alkyne functionalities: the longer the connector, the more favored the lactamization process vs hydroamination.

Catalysis by [Ga4 L6 ]12- Metallocage on the Nazarov Cyclization: The Basicity of Complexed Alcohol is Key.

The Nazarov cyclization is investigated in solution and within K12 [Ga4 L6 ] supramolecular organometallic cage by means of computational methods. The reaction needs acidic condition in solution but works at neutral pH in the presence of the metallocage. The reaction steps for the process are analogous in both media: (a) protonation of the alcohol group, (b) water loss and (c) cyclization. The relative Gibbs energies of all the steps are affected by changing the environment from solvent to the metallocage. The first step in the mechanism, the alcohol protonation, turns out to be the most critical one for the acceleration of the reaction inside the metallocage. In order to calculate the relative stability of protonated alcohol inside the cavity, we propose a computational scheme for the calculation of basicity for species inside cavities and can be of general use. These results are in excellent agreement with the experiments, identifying key steps of catalysis and providing an in-depth understanding of the impact of the metallocage on all the reaction steps.

Origin of the Rate Acceleration in the C-C Reductive Elimination from Pt(IV)-complex in a [Ga4 L6 ]12- Supramolecular Metallocage.

The reductive elimination on [(Me3 P)2 Pt(MeOH)(CH3 )3 ]+ , 2P, complex performed in MeOH solution and inside a [Ga4 L6 ]12- metallocage are computationally analysed by mean of QM and MD simulations and compared with the mechanism of gold parent systems previously reported [Et3 PAu(MeOH)(CH3 )2 ]+ , 2Au. The comparative analysis between the encapsulated Au(III) and Pt(IV)-counterparts shows that there are no additional solvent MeOH molecules inside the cavity of the metallocage for both systems. The Gibbs energy barriers for the 2P reductive elimination calculated at DFT level are in good agreement with the experimental values for both environments. The effect of microsolvation and encapsulation on the rate acceleration are evaluated and shows that the latter is far more relevant, conversely to 2Au. Energy decomposition analysis indicates that the encapsulation is the main responsible for most of the energy barrier reduction. Microsolvation and encapsulation effects are not equally contributing for both metal systems and consequently, the reasons of the rate acceleration are not the same for both metallic systems despite the similarity between them.

Modeling Kinetics and Thermodynamics of Guest Encapsulation into the [M4L6]12- Supramolecular Organometallic Cage.

The encapsulation of molecular guests into supramolecular hosts is a complex molecular recognition process in which the guest displaces the solvent from the host cavity, while the host deforms to let the guest in. An atomistic description of the association would provide valuable insights on the physicochemical properties that guide it. This understanding may be used to design novel host assemblies with improved properties (i.e., affinities) toward a given class of guests. Molecular simulations may be conveniently used to model the association processes. It is thus of interest to establish efficient protocols to trace the encapsulation process and to predict the associated magnitudes ΔGbind and ΔGbind⧧. Here, we report the calculation of the Gibbs energy barrier and Gibbs binding energy by means of explicit solvent molecular simulations for the [Ga4L6]12- metallocage encapsulating a series of cationic molecules. The ΔGbind⧧ for encapsulation was estimated by means of umbrella sampling simulations. The steps involved were identified, including ion-pair formation and naphthalene rotation (from L ligands of the metallocage) during the guest's entrance. The ΔGbind values were computed using the attach-pull-release method. The results reveal the sensitivity of the estimates on the force field parameters, in particular on atomic charges, showing that higher accuracy is obtained when charges are derived from implicit solvent quantum chemical calculations. Correlation analysis identified some indicators for the binding affinity trends. All computed magnitudes are in very good agreement with experimental observations. This work provides, on one side, a benchmarked way to computationally model a highly charged metallocage encapsulation process. This includes a nonstandard parameterization and charge derivation procedure. On the other hand, it gives specific mechanistic information on the binding processes of [Ga4L6]12- at the molecular level where key motions are depicted. Taken together, the study provides an interesting option for the future design of metal-organic cages.

Computational Analysis on the Pd-Catalyzed C-N Coupling of Ammonia with Aryl Bromides Using a Chelate Phosphine Ligand.

The Buchwald-Hartwig amination of arylhalides with the Pd-Josiphos complex is a very useful process for the generation of primary amines using ammonia as a reactant. Density-functional theory (DFT) calculations are carried out to examine the reaction mechanism for this process. Although the general mechanism for the C-N cross-coupling reaction is known, there are still some open questions regarding the effect of a chelate phosphine ligand and the role of the base in the process. Reaction pathways involving the release of one of the arms of the phosphine ligand are compared with those where the chelate phosphine remains fully coordinated. Conformational analysis for the complex with the open chelate phosphine is required to properly evaluate the proposed pathways. The role played by the added base (t-BuO-) as a possible ligand or just as a base was also evaluated. The understanding of all of these aspects allowed us to propose a complete reaction mechanism for the Pd-catalyzed C-N coupling of arylhalides with ammonia using the chelate Josiphos ligand.

Reaction Rate Inside the Cavity of [Ga4 L6 ]12- Supramolecular Metallocage is Regulated by the Encapsulated Solvent.

In the present study the dependence of the reaction rate of carbon-carbon reductive elimination from R3 PAu(MeOH)(CH3 )2 (R=Me, Et) complexes inside [Ga4 L6 ]12- metallocage on the nature of the phosphine ligand is investigated by computational means. The reductive elimination mechanism is analyzed in methanol solution and inside the metallocage. Classical molecular dynamics simulations reveal that the smaller the gold complex (which depends on the phosphine ligand size) the larger the number of solvent molecules encapsulated. The size of the phosphine ligands defines the space that is left available inside the cavity that can be occupied by solvent molecules. The Gibbs energy barriers calculated at DFT level, in excellent agreement with experiment both in solution and in the metallocage, show that the presence/absence of explicit solvent molecules inside the cavity significantly modifies the reaction rate.

A resorcinarene-based tetrabenzoimidazolylidene complex of rhodium.

A tetrabenzoimidazolium-resorcinarene cavitand was used for the preparation of a tetra-benzoimidazolylidene of rhodium, which is unprecedented in the field of poly-NHC metal complexes. Both the experimental and computational analyses of the molecule reveal a distorted vase conformation as the most stable one, although several non-interconverting conformational isomers due to the restricted rotation about the Rh-C(carbene) bond coexist in the product. There is a fluxional behaviour involving the vase-kite interconversion. The main interactions between the arms of the cavitand are mostly concentrated on the terminal organometallic fragments attached to NHC, along with those between -CH3 and a N-heterocyclic carbene ring from benzoimidazoles.

Catalytic Regioselective Isomerization of 2,2-Disubstituted Oxetanes to Homoallylic Alcohols.

The selective isomerization of strained heterocyclic compounds is an important tool in organic synthesis. An unprecedented regioselective isomerization of 2,2-disubstituted oxetanes into homoallylic alcohols is described. The use of tris(pentafluorophenyl)borane (B(C6 F5 )3 ), a commercially available Lewis acid was key to obtaining good yields and selectivities since other Lewis acids afforded mixtures of isomers and substantial polymerization. The reaction took place under exceptionally mild reaction conditions and very low catalyst loading (0.5 mol %). DFT calculations disclose the mechanistic features of the isomerization and account for the high selectivity displayed by the B(C6 F5 )3 catalyst. The synthetic applicability of the new reaction is demonstrated by the preparation of γ-chiral alcohols using iridium-catalyzed asymmetric hydrogenation.

Skeletal diversity in Pt- and Au-catalyzed annulations of allenedienes: dissecting unconventional mechanistic pathways.

We describe the discovery of unprecedented annulation processes of 1,7-allenedienes, promoted by Pt or Au catalysts. These transformations revealed mechanistic pathways that had not been previously observed in reactions involving carbophilic catalysis. In particular, we have found that allenedienes bearing a silyl ether in the carbon tether connecting the diene and the allene divergently afford cyclopropane-embedded tricyclic derivatives, 6,6-fused bicarbocyclic products or 5,6-fused bicarbocyclic systems, depending on the type of Au or Pt catalyst used. We have carried out experimental and computational studies that shed light on the mechanistic reasons behind this rich and unusual skeletal divergence, and provide new lessons on the drastic influence of platinum ancillary ligands on the reaction outcome.

A Reversible Phase Transition of 2D Coordination Layers by B-H∙∙∙Cu(II) Interactions in a Coordination Polymer.

Materials that combine flexibility and open metal sites are crucial for myriad applications. In this article, we report a 2D coordination polymer (CP) assembled from CuII ions and a flexible meta-carborane-based linker [Cu2(L1)2(Solv)2]•xSolv (1-DMA, 1-DMF, and 1-MeOH; L1: 1,7-di(4-carboxyphenyl)-1,7-dicarba-closo-dodecaborane). 1-DMF undergoes an unusual example of reversible phase transition on solvent treatment (i.e., MeOH and CH2Cl2). Solvent exchange, followed by thermal activation provided a new porous phase that exhibits an estimated Brunauer-Emmett-Teller (BET) surface area of 301 m2 g-1 and is capable of a CO2 uptake of 41 cm3 g-1. The transformation is reversible and 1-DMF is reformed on addition of DMF to the porous phase. We provide evidence for the reversible process being the result of the formation/cleavage of weak but attractive B-H∙∙∙Cu interactions by a combination of single-crystal (SCXRD), powder (PXRD) X-ray diffraction, Raman spectroscopy, and DFT calculations.

Microsolvation and Encapsulation Effects on Supramolecular Catalysis: C-C Reductive Elimination inside [Ga4L6]12- Metallocage.

The host effect of the supramolecular [Ga4L6]12- tetrahedral metallocage on reductive elimination of substrate by encapsulated Au(III) complexes is investigated by means of computational methods. The behavior of the reactants in solution and within the metallocage is initially evaluated by means of classical molecular dynamics simulations. These results guided the selection of proper computational models to describe the reaction in solution and inside the metallocage at the DFT level. The calculated Gibbs energy barriers are in very good agreement with experiment both in solution and inside the metallocage. The analysis in solution revealed that microsolvation around the Au(III) complex increases the Gibbs energy barrier. The analysis within the metallocage shows that its encapsulation favors the reaction. The process can be formally described as removing explicit microsolvation around the gold complex and encapsulating the metal complex inside the metallocage. Both processes are important for the reaction, but the removal of the solvent molecules surrounding the Au(III) metal complex is fundamental for the reduction of the reaction barrier. The energy decomposition analysis of the barrier among strain, interaction, and thermal terms shows that strain term is very low whereas the contribution of thermal (entropic) effects is moderate. Interestingly, the key term responsible for reducing the Gibbs energy barrier is the interaction. This term can be mainly associated with electrostatic interactions in agreement with previous examples in the literature.

GARLEEK: Adding an extra flavor to ONIOM.

The ONIOM method, developed in the group of Keiji Morokuma, is one of the most successful examples of quantum mechanics/molecular mechanics (QM/MM) treatments, and of multilayer methods in general. Its implementation in the Gaussian program package is in particular widely used. This implementation has access to the wide variety of QM methods available in Gaussian, but is limited to only three specific force fields. The current article presents the GARLEEK interface, which expands the availability of molecular mechanics methods to the wide variety of force fields available in MM packages. The focus is in the simple installation and use. Two examples of the performance of the interface with selected systems are provided. GARLEEK is MIT-licensed and freely available at https://github.com/insilichem/garleek. © 2018 Wiley Periodicals, Inc.

Direct Asymmetric Hydrogenation of N-Methyl and N-Alkyl Imines with an Ir(III)H Catalyst.

A novel cationic [IrH(THF)(P,N)(imine)] [BArF] catalyst containing a P-stereogenic MaxPHOX ligand is described for the direct asymmetric hydrogenation of N-methyl and N-alkyl imines. This is the first catalytic system to attain high enantioselectivity (up to 94% ee) in this type of transformation. The labile tetrahydrofuran ligand allows for effective activation and reactivity, even at low temperatures. Density functional theory calculations allowed the rationalization of the stereochemical course of the reaction.

Iridium-Catalyzed Isomerization of N-Sulfonyl Aziridines to Allyl Amines.

The Crabtree's reagent catalyzes the isomerization of N-sulfonyl 2,2-disubstituted aziridines to allyl amines. The selectivity of allyl amine vs imine is very high (up to 99/1). The unprecedented isomerization takes place in mild conditions without activation of the catalyst by hydrogen. The mechanism has been studied computationally by DFT calculations; instead of the usual hydrogenation of COD, the catalytic species is formed by a loss of the pyridine ligand. Approaching of aziridine to this unsaturated species leads to a carbocation intermediate through a low energy barrier. A metal-mediated tautomerization involving sequentially γ-H elimination and N-H reductive elimination affords selectively the allyl amine. The readiness of the CγH bond to participate in the H elimination step accounts for the selectivity toward the allyl amine product.

Calculation of Reaction Free Energies in Solution: A Comparison of Current Approaches.

The result of the application of different approaches based on the ideal gas/rigid rotor/harmonic oscillator (IGRRHO) model, commonly used in popular software packages, for the calculation of free energies in solution is compared with that of ab initio molecular dynamics for a process involving ligand exchange in palladium complexes. The IGRRHO-based approaches considered differ in most cases in the extent to which the rotational and translational contributions are included in the correction. Our study supports the use the free energy values directly obtained from dispersion-corrected DFT functionals without any correction or with minor corrections at most.

Mechanistic implications of the enantioselective addition of alkylzinc reagents to aldehydes catalyzed by nickel complexes with α-amino amide ligands.

The enantioselective alkylation of aldehydes catalysed by nickel(ii)-complexes derived from α-amino amides was studied by means of density functional theory (DFT) and ONIOM (B3LYP:UFF) calculations. A mechanism was proposed in order to investigate the origin of enantioselectivity. The chirality-determining step for the alkylation was the formation of the intermediate complexes with the involvement of a 5/4/4-fused tricyclic transition state. The predominant products predicted theoretically were of (S)-configuration, in good agreement with experimental observations. The scope of the reaction was examined and high yields and enantioselectivities were observed for the enantioselective addition of Et2Zn and Me2Zn to aromatic and aliphatic aldehydes.

Rim, Side Arms, and Cavity: Three Sites for the Recognition of Anions by Tetraazolium Resorcinarene Cavitands.

Two tetrabenzoimidazolium-resorcinarene cavitands were prepared and used for the recognition of chloride, bromide, iodide, cyanide, nitrate, perchlorate, hexanoate, benzenesulfonate, and p-toluenesulfonate. Binding affinities of the two cavitands were determined by 1 H NMR titration and computational analysis. The observed spectral changes were related to specific interaction sites, which were supported by the computational studies. In the case of the C2-H tetrabenzoimidazolium-resorcinarene, the recognition region of the inorganic anions and hexanoate was located at the rim of the cavitand, although chloride and bromide also interacted with the aromatic C-H bonds located between adjacent arms of the cavitand. By contrast, the recognition of the two anions with an aromatic ring (benzenesulfonate and p-toluenesulfonate) results from encapsulation of the aromatic part of the anions inside the hydrophobic cavity of the host. In the case of the C2-Me tetrabenzoimizazolium-resorcinarene receptor, the ability of the molecule to bind all inorganic anions and hexanoate was suppressed, but the receptor maintained its ability to strongly bind benzenesulfonate and p-toluenesulfonate. This is interpreted in terms of suppression of the ability of the cavitand to form hydrogen bonds at the rim of the molecule due to replacement of the C2-H proton by a methyl group, while the hydrophobic pocket of the molecule maintains its binding abilities.

Acid Activation in Phenyliodine Dicarboxylates: Direct Observation, Structures, and Implications.

The use of the hypervalent iodine reagents in oxidative processes has become a staple in modern organic synthesis. Frequently, the reactivity of λ3 iodanes is further enhanced by acids (Lewis or Brønsted). The origin of such activation, however, has remained elusive. Here, we use the common combination of PhI(OAc)2 with BF3·Et2O as a model to fully explore this activation phenomenon. In addition to the spectroscopic assessment of the dynamic acid-base interaction, for the first time the putative PIDA·BF3 complex has been isolated and its structure determined by X-ray diffraction. Consequences of such activation are discussed from a structural and electronic (DFT) points of views, including the origins of the enhanced reactivity.

First-Principles Molecular Dynamics Studies of Organometallic Complexes and Homogeneous Catalytic Processes.

Computational chemistry is a valuable aid to complement experimental studies of organometallic systems and their reactivity. It allows probing mechanistic hypotheses and investigating molecular structures, shedding light on the behavior and properties of molecular assemblies at the atomic scale. When approaching a chemical problem, the computational chemist has to decide on the theoretical approach needed to describe electron/nuclear interactions and the composition of the model used to approximate the actual system. Both factors determine the reliability of the modeling study. The community dedicated much effort to developing and improving the performance and accuracy of theoretical approaches for electronic structure calculations, on which the description of (inter)atomic interactions rely. Here, the importance of the model system used in computational studies is highlighted through examples from our recent research focused on organometallic systems and homogeneous catalytic processes. We show how the inclusion of explicit solvent allows the characterization of molecular events that would otherwise not be accessible in reduced model systems (clusters). These include the stabilization of nascent charged fragments via microscopic solvation (notably, hydrogen bonding), transfer of charge (protons) between distant fragments mediated by solvent molecules, and solvent coordination to unsaturated metal centers. Furthermore, when weak interactions are involved, we show how conformational and solvation properties of organometallic complexes are also affected by the explicit inclusion of solvent molecules. Such extended model systems may be treated under periodic boundary conditions, thus removing the cluster/continuum (or vacuum) boundary, and require a statistical mechanics simulation technique to sample the accessible configurational space. First-principles molecular dynamics, in which atomic forces are computed from electronic structure calculations (namely, density functional theory), is certainly the technique of choice to investigate chemical events in solution. This methodology is well established and thanks to advances in both algorithms and computational resources simulation times required for the modeling of chemical events are nowadays accessible, though the computational requirements use to be high. Specific applications reviewed here include mechanistic studies of the Shilov and Wacker processes, speciation in Pd chemistry, hydrogen bonding to metal centers, and the dynamics of agostic interactions.