Mechanochemical Felkin-Anh Model: Achieving Forbidden Reaction Outcomes with Mechanical Force.

Anti-Felkin-Anh diastereoselectivity can be achieved for nucleophilic additions to α-chiral ketones upon stretching the ketone with a mechanical pulling force. Herein, a mechanochemical Felkin-Anh model is proposed for predicting the outcome of a nucleophilic addition to an α-chiral ketone. Essentially, the fully stretched chiral ketone has one substituent shielding each side of the carbonyl, in contrast to the Felkin-Anh model, in which free rotation around a bond is required to achieve the two rotamers of the ketone. Depending on the pulling scenario, either Felkin-Anh or anti-Felkin-Anh diastereoselectivity is obtained. The model is entirely based on the distance between the pulling points, which is maximized in the anti-periplanar arrangement. The major diastereomer is associated with the approach with the least steric interactions. The intuitive model is validated by means of mechanochemical density functional theory calculations. Importantly, the ketone is fully stretched in the sub 1 nN force regime, thus minimizing the risk of undesired homolytic bond rupture. Moreover, the mechanical force is not used for lowering the reaction barriers associated with the nucleophilic addition; instead, it is solely applied for locking the conformation of a molecule and provoking otherwise inaccessible reaction pathways on the force-modified potential energy surface.

Wandering through quantum-mechanochemistry: from concepts to reactivity and switches.

Mechanochemistry has experienced a renaissance in recent years witnessing, at the molecular level, a remarkable interplay between theory and experiment. Molecular mechanochemistry has welcomed a broad spectrum of quantum-chemical methods to evaluate the influence of an external mechanical force on molecular properties. In this contribution, an overview is given on recent work on quantum mechanochemistry in the Brussels Quantum Chemistry group (ALGC). The effect of an external force was scrutinized both in fundamental topics, like reactivity descriptors in Conceptual DFT, and in applied topics, such as designing molecular force probes and tuning the stereoselectivity of certain types of reactions. In the conceptual part, a brief overview of the techniques introducing mechanical forces into a quantum-mechanical description of a molecule is followed by an introduction to conceptual DFT. The evolution of the electronic chemical potential (or electronegativity), chemical hardness and electrophilicity are investigated when a chemical bond in a series of diatomics is put under mechanical stress. Its counterpart, the influence of mechanical stress on bond angles, is analyzed by varying the strain present in alkyne triple bonds by applying a bending force, taking the strain promoted alkyne-azide coupling cycloaddition as an example. The increase of reactivity of the alkyne upon bending is probed by Fukui functions and the local softness. In the applied part, a new molecular force probe is presented based on an intramolecular 6π-electrocyclization in constrained polyenes operating under thermal conditions. A cyclic process is conceived where ring opening and closure are triggered by applying or removing an external pulling force. The efficiency of mechanical activation strongly depends on the magnitude of the applied force and the distance between the pulling points. The idea of pulling point distances as a tool to identify new mechanochemical processes is then tested in [28]hexaphyrins with an intricate equilibrium between Möbius aromatic and Hückel antiaromatic topologies. A mechanical force is shown to trigger the interconversion between the two topologies, using the distance matrix as a guide to select appropriate pulling points. In a final application, the Felkin-Anh model for the addition of nucleophiles to chiral carbonyls under the presence of an external mechanical force is scrutinized. By applying a force for restricting the conformational freedom of the chiral ketone, otherwise inaccessible reaction pathways are promoted on the force-modified potential energy surfaces resulting in a diastereoselectivity different from the force-free reaction.

Mechanochemically Triggered Topology Changes in Expanded Porphyrins.

A hitherto unexplored class of molecules for molecular force probe applications are expanded porphyrins. This work proves that mechanical force is an effective stimulus to trigger the interconversion between Hückel and Möbius topologies in [28]hexaphyrin, making these expanded porphyrins suitable to act as conformational mechanophores operating at mild (sub-1 nN) force conditions. A straightforward approach based on distance matrices is proposed for the selection of pulling scenarios that promote either the planar Hückel topology or the three lowest lying Möbius topologies. This approach is supported by quantum mechanochemical calculations. Force distribution analyses reveal that [28]hexaphyrin selectively allocates the external mechanical energy to molecular regions that trigger Hückel-Möbius interconversions, explaining why certain pulling scenarios favor the Hückel two-sided topology and others favor Möbius single-sided topologies. The meso-substitution pattern on [28]hexaphyrin determines whether the energy difference between the different topologies can be overcome by mechanical activation.

Designing Force Probes Based on Reversible 6π-Electrocyclizations in Polyenes Using Quantum Chemical Calculations.

The conjugated π-system in polyenes can be interrupted by electrocyclic ring-closure reactions. In this work, this 6π-electrocylization is shown by means of density functional calculations to be reversible by the application of an external mechanical pulling force at the terminal ends of the interrupted polyene chain. The test systems were constrained in a fused ring system, thus locking the orientation of three π-bonds and generally promoting 6π-electrocyclic ring-closure reactions. For several systems, the forward reaction is exergonic and the corresponding reaction barrier is comparable to those reported in the literature. The reverse reaction is triggered by an external pulling force of 2 nN (nano-Newton) or less and also becomes exergonic in all investigated polyenes under these force conditions. Moreover, it proceeds via a low reaction barrier when a pulling force of 2 nN is active, indicating that the mechanical force is an efficient stimulus for triggering ring-opening reactions. Analysis of the strain energy induced by this mechanical force confirms an optimal activation of the corresponding C-C σ-bond that breaks upon ring opening when the pulling positions are located on the polyene chain.

Alkaline Earth Metals Activate N2 and CO in Cubic Complexes Just Like Transition Metals: A Conceptual Density Functional Theory and Energy Decomposition Analysis Study.

Following the recent discovery of stable octa-coordinated alkaline earth metals with N2 and CO, the role of group II metals in the catalytic reduction of these ligands by means of density functional theory (DFT) calculations and conceptual DFT-based reactivity indices is investigated. Cubic group IV and octahedral group VI transition metal complexes as well as the free ligands are computed for reference. The outer and most accessible atoms of N2 and CO become much more nucleophilic and electrophilic in all complexes, relevant for N2 fixation, as probed by the Fukui function and local softness. Within one row of the periodic table, the alkaline earth complexes often show the strongest activation. On the contrary, the electrostatic character is found to be virtually unaffected by complexation. Trends in the soft frontier orbital and hard electrostatic character are in agreement with calculated proton affinities and energy decomposition analyses of the protonated structures, demonstrating the dominance of the soft (HOMO-LUMO) orbital interactions.

Ambident Nucleophilic Substitution: Understanding Non-HSAB Behavior through Activation Strain and Conceptual DFT Analyses.

The ability to understand and predict ambident reactivity is key to the rational design of organic syntheses. An approach to understand trends in ambident reactivity is the hard and soft acids and bases (HSAB) principle. The recent controversy over the general validity of this principle prompted us to investigate the competing gas-phase SN 2 reaction channels of archetypal ambident nucleophiles CN- , OCN- , and SCN- with CH3 Cl (SN 2@C) and SiH3 Cl (SN 2@Si), using DFT calculations. Our combined analyses highlight the inability of the HSAB principle to correctly predict the reactivity trends of these simple, model reactions. Instead, we have successfully traced reactivity trends to the canonical orbital-interaction mechanism and the resulting nucleophile-substrate interaction energy. The HOMO-LUMO orbital interactions set the trend in both SN 2@C and SN 2@Si reactions. We provide simple rules for predicting the ambident reactivity of nucleophiles based on our Kohn-Sham molecular orbital analysis.

Implementing the mechanical force into the conceptual DFT framework: understanding and predicting molecular mechanochemical properties.

In molecular mechanochemistry, the chemical properties of a reactant system are modified through the absorption of mechanical energy at the single molecular level, as opposed to thermochemistry (heat), electrochemistry (electrical current) and photochemistry (light). In this proof-of-concept paper, we present a theoretical approach to rationalize and predict the change in chemical properties and concepts when an external stretching force, Fext, is applied to a chemical bond. By considering perturbations to the total molecular energy due to Fext, a series of mechanochemical response indices are obtained within the framework of conceptual density functional theory (DFT). Trends in these indices monitoring the change in hardness and electrophilicity among others, were rationalized for a number of diatomic molecules in terms of the ground-state geometry and the frontier molecular orbitals of the molecule. Finally, we present a set of rules that can be easily adopted for polyatomic molecules to predict the changing chemical reactivity of single molecules when subjected to an external force.

A Quantum Chemical Deep-Dive into the π-π Interactions of 3-Methylindole and Its Halogenated Derivatives-Towards an Improved Ligand Design and Tryptophan Stacking.

Non-covalent π-π stacking interactions often play a key role in the stability of the secondary and tertiary structures of peptides and proteins, respectively, and can be a means of ensuring the binding of ligands within protein and enzyme binding sites. It is generally accepted that minor structural changes to the aromatic ring, such as substitution, can have a large influence on these interactions. Nevertheless, a thorough understanding of underpinning phenomena guiding these key interactions is still limited. This is especially true for larger aromatic structures. To expand upon this knowledge, elaborate ab initio calculations were performed to investigate the effect of halogenation on the stability of 3-methylindole stacking. 3-Methylindole served as a representation of the tryptophan side chain, and is a frequently used motif in drug design and development. Moreover, an expression is derived that is able to accurately predict the interaction stability of stacked halogenated 3-methylindole dimers as well as halogenated toluene dimers, based on monomer level calculated DFT descriptors. We aim for this expression to provide the field with a straightforward and reliable method to assess the effect of halogenation on the π-π stacking interactions between aromatic scaffolds.

The hunt for reactive alkynes in bio-orthogonal click reactions: insights from mechanochemical and conceptual DFT calculations.

In our effort to implement the mechanical force used to activate single molecules in mechanochemistry in the context of conceptual density functional theory, we present a theoretical investigation of strained alkynes for rationalizing structural trends as well as the reactivity of cyclic alkynes that are of great importance in in vivo click reactions. The strain on the triple bond in cyclic alkynes is modeled by angular constraints in a 2-butyne fragment and the corresponding bending force is calculated by means of an extended COGEF (constrained geometries simulate external forces) model. In general, the force required to bend the triple bond is smaller with electron-withdrawing groups on the propargylic C-atom, which elegantly results in smaller angles around the triple bond in cyclic alkynes with such substitution pattern. By means of conceptual DFT descriptors, the electrophilic and nucleophilic character of bent triple bonds was investigated revealing moderate activation for small distortions from the linear geometry (0° to 15°) and a drastically more reactive π-space if the triple bond is bent further. This analysis of the intrinsic reactivity of the triple bond is in line with experimental observations, explaining the reactive nature of cyclooctynes and cycloheptynes, whereas larger cyclic systems do not drastically activate the triple bond.