Similarities and Differences in Benzene Reduction with Ca, Sr, Yb and Sm: Strong Evidence for Tetra-Anionic Benzene

Inverse sandwich complexes of Yb and Sm stabilized by a bulky β-diketiminate (BDI) ligand have been prepared: (BDI)Ln(η6,η6-C6H6)Ln(BDI); Ln=lanthanide. Coordinated benzene ligands can be neutral, di-anionic or, often controversially discussed, even tetra-anionic. The formal charge on benzene is correlated to assignment of the metal oxidation state which generally poses a problem. Herein, we take advantage of the structural similarities found when comparing CaII with YbII, and SrII with SmII complexes. In this work, we found an excellent overlap of the Ca/Yb inverse sandwich structures but a striking difference for the Sr/Sm pair. The much shorter Sm-N and Sm-C6H6 distances are strong evidence for a SmIII-benzene-4-SmIII assignment. This was further supported by NMR spectroscopy, magnetic susceptibility, reactivity and comprehensive computational investigation.

New Perspectives on Delocalization Pathways in Aromatic Molecular Chameleons

This study comprehensively analyzes the magnetically induced current density of polycyclic compounds labeled as "aromatic chameleons" since they can arrange their π-electrons to exhibit aromaticity in both the ground and the lowest triplet state. These compounds comprise benzenoid moieties fused to a central skeleton with 4n π-electrons and traditional magnetic descriptors are biased due to the superposition of local magnetic responses. In the S0 state, our analysis reveals that the molecular constituent fragments preserve their (anti)aromatic features in agreement with two types of resonant structures: one associated with aromatic benzenoids and the other with a central antiaromatic ring. Regarding the T1 state, a global and diatropic ring current is revealed. Our aromaticity study is complemented with advanced electronic and geometric descriptors to consider different aspects of aromaticity, particularly important in the evaluation of excited state aromaticity. Remarkably, these descriptors consistently align with the general features on the main delocalization pathways in polycyclic hydrocarbons consisting of fused 4n π-electron rings. Moreover, our study demonstrates an inverse correlation between the singlet-triplet energy difference and the antiaromatic character of the central ring in S0.

A combined analysis of two prospective randomised studies exploring the impact of extended post-radiation temozolomide on survival outcomes in newly diagnosed glioblastoma

PURPOSE

The optimal duration of post-radiation temozolomide in newly diagnosed glioblastoma remains unclear, with no published phase III randomised trials. Standard-of-care stipulates 6 months. However, in routine care, it is often extended to 12 months, despite lacking robust supporting data.

METHODS

GEINO14-01 (Spain) and EX-TEM (Australia) studies enrolled glioblastoma patients without progression at the end of 6 months post-radiation temozolomide. Participants were randomised 1:1 to six additional months of temozolomide or observation. Primary endpoint was 6-month progression free survival from date of randomisation (6mPFS). Secondary endpoints included overall survival (OS) and toxicity. 204 patients were required to detect an improvement in 6mPFS from 50 to 60% (80% power). Neither study recruited sufficient patients. We performed a combined analysis of individual patient data.

RESULTS

205 patients were recruited: 159 in GEINO14-01 (2014-2018) and 46 in EX-TEM (2019-2022). Median follow-up was 20.0 and 14.5 months. Baseline characteristics were balanced. There was no significant improvement in 6mPFS (57.2% vs 64.0%, OR0.75, p = 0.4), nor across any subgroups, including MGMT methylated; PFS (HR0.92, p = 0.59, median 7.8 vs 9.7 months); or OS (HR1.03, p = 0.87, median 20.1 vs 19.4 months). During treatment extension, 64% experienced any grade adverse event, mainly fatigue and gastrointestinal (both 54%). Only a minority required treatment changes: 4.5% dose delay, 7.5% dose reduction, 1.5% temozolomide discontinuation.

CONCLUSION

For glioblastoma patients, extending post-radiation temozolomide from 6 to 12 months is well tolerated but does not improve 6mPFS. We could not identify any subset that benefitted from extended treatment. Six months should remain standard-of-care.

T-Shaped Palladium and Platinum {MNO}10 Nitrosyl Complexes

The synthesis and characterization of a homologous series of T-shaped {MNO}10 nitrosyl complexes of the form [M(PR3)2(NO)]+ (M = Pd, Pt; R = tBu, Ad) are reported. These diamagnetic nitrosyls are obtained from monovalent or zerovalent precursors by treatment with NO and NO+, respectively, and are notable for distinctly bent M-NO angles of ∼123° in the solid state. Adoption of this coordination mode in solution is also supported by the analysis of isotopically enriched samples by 15N NMR spectroscopy. Effective oxidation states of M0/NO+ are calculated, and metal-nitrosyl bonding has been interrogated using DFT-based energy decomposition analysis techniques. While a linear nitrosyl coordination mode was found to be electronically preferred, the M-NO and P-M-P angles are inversely correlated to the extent that binding in this manner is prevented by steric repulsion between the bulky ancillary phosphine ligands.

Aromaticity and Magnetic Behavior in Benzenoids: Unraveling Ring Current Combinations

Nowadays, an active research topic is the connection between Clar's rule, aromaticity, and magnetic properties of polycyclic benzenoid hydrocarbons. In the present work, we employ a meticulous magnetically induced current density analysis to define the net current flowing through any cyclic circuit, connecting it to aromaticity based on the ring current concept. Our investigation reveals that the analyzed polycyclic systems display a prominent global ring current, contrasting with subdued semi-local and local ring currents. These patterns align with Clar's aromatic π-sextets only in cases where migrating π-sextet structures are invoked. The results of this study will enrich our comprehension of aromaticity and magnetic behavior in such systems, offering significant insights into coexisting ring current circuits in these systems.

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.

Isomerism tunes the diradical character of difluorenopyrroles at constant Hückel-level anti-aromaticity

A new diradical based on diindenocarbazole or difluorenopyrrole was synthesized and experimentally characterized by optical, electrochemical, and magnetic techniques, as well as quantum chemical calculations. The isomerism of these structures tunes the diradical character and the associated properties, representing a unique case of such important modulation. A full study of the electronic structure was carried out considering the perturbative interactions between different canonical forms as well as the anti-aromatic character of the molecular cores. Such a study reveals how we can tune diradical character simply by reorganizing the bonding patterns at constant chemical costs (composition).

Application of Inverse Design Approaches to the Discovery of Nonlinear Optical Switches

Molecular switches, in which a stimulus induces a large and reversible change in molecular properties, are of significant interest in the domain of photonics. Due to their commutable redox states with distinct nonlinear optical (NLO) properties, hexaphyrins have emerged as a novel platform for multistate switches in nanoelectronics. In this study, we employ an inverse design algorithm to find functionalized 26R→28R redox switches with maximal βHRS contrast. We focus on the role of core modifications, since a synergistic effect with meso-substitutions was recently found for the 30R-based switch. In contrast to these findings, the inverse design optima and subsequent database analysis of 26R-based switches confirm that core modifications are generally not favored when high NLO contrasts are targeted. Moreover, while push-pull combinations enhance the NLO contrast for both redox switches, they prefer a different arrangement in terms of electron-donating and electron-withdrawing functional groups. Finally, we aim at designing a three-state 26R→28R→ 30R switch with a similar NLO response for both ON states. Even though our best-performing three-state switch follows the design rules of the 30R-based component, our chemical compound space plots show that well-performing three-state switches can be found in regions shared by high-responsive 26R and 30R structures.

Designing hexaphyrins for high-potential NLO switches: the synergy of core-modifications and meso-substitutions

Due to the enormous size of the chemical compound space, usually only small regions are traversed with traditional direct molecular design approaches making the discovery for novel functionalized molecules for nonlinear optical applications challenging. By applying inverse molecular design algorithms, we aim to efficiently explore larger regions of the compound space in search of promising hexaphyrin-based molecular switches as measured by their first-hyperpolarizability (β_HRS) contrast. We focus on the 28R → 30R switch with a functionalization pattern allowing for centrosymmetric OFF states yielding zero β_HRS response. This switch is particularly challenging as full meso-substitution with a single type of functional group or core-modifications result in almost no contrast enhancement. We carried out four inverse design procedures during which two sets of core-modifications and three sets of meso-substitutions sites were systematically optimized. All 4 optimal switches are characterized by a mix of meso-substitutions and core-modifications, of which the best performing switch yields a 10-fold improvement over the parent macrocycle. Throughout the inverse design procedures, we collected and analyzed a database biased towards high NLO contrasts that contains 277 different patterns for hexaphyrin-based switches. We derived three design rules to obtain highly functional 28R → 30R NLO switches: (I) a combination of 2 strong EWG and 1 EDG group is the ideal recipe for increasing the NLO contrast, though their position also plays an important role. (II) The type of core-modification is less important when only the diagonal positions are core-modified. Switches with 4 core-modifications show a clear preference for oxygen. (III) Keeping centrosymmetry in the OFF state remains highly beneficial given the investigated functionalization pattern. Finally, we have demonstrated that combining meso-substitutions with core-modifications can synergistically improve the NLO contrast.

Understanding the Reactivity of Supported Late Transition Metals on a Bare Anatase (101) Surface: A Periodic Conceptual DFT Investigation

The rapidly growing interest for new heterogeneous catalytic systems providing high atomic efficiency along with high stability and reactivity triggered an impressive progress in the field of single-atom catalysis. Nevertheless, unravelling the factors governing the interaction strength between the support and the adsorbed metal atoms remains a major challenge. Based on periodic density functional theory (DFT) calculations, this paper provides insight into the adsorption of single late transition metals on a defect-free anatase surface. The obtained adsorption energies fluctuate, with the exception of Pd, between -3.11 and -3.80 eV and are indicative of a strong interaction. Depending on the considered transition metal, we could attribute the strength of this interaction with the support to i) an electron transfer towards anatase (Ru, Rh, Ni), ii) s-d orbital hybridisation effects (Pt), or iii) a synergistic effect between both factors (Fe, Co, Os, Ir). The driving forces behind the adsorption were also found to be strongly related to Klechkowsky's rule for orbital filling. In contrast, the deviating behaviour of Pd is most likely associated with the lower dissociation enthalpy of the Pd-O bond. Additionally, the reactivity of these systems was evaluated using the Fermi weighted density of states approach. The resulting softness values can be clearly related to the electron configuration of the catalytic systems as well as with the net charge on the transition metal. Finally, these indices were used to construct a model that predicts the adsorption strength of CO on these anatase-supported d-metal atoms. The values obtained from this regression model show, within a 95 % probability interval, a correlation of 84 % with the explicitly calculated CO adsorption energies.

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.

Ultrasound-Guided Percutaneous Cryoanalgesia for Pectus Excavatum: When Should It be Applied?

INTRODUCTION

 The addition of ultrasound-guided percutaneous cryoanalgesia (PCr) for pain management after pectus excavatum (PE) surgery offers a new and advantageous approach. Our aim is to describe our experience with PCr applied on the same day, 24 hours, and 48 hours prior to PE surgery.

MATERIAL AND METHODS

 Prospective pilot study in patients undergoing ultrasound-guided PCr (2019-2022) was divided into three groups: PCr on the same day of surgery (PCrSD), PCr 24 hours before (PCr24), and PCr 48 hours before (PCr48). We describe the application of technique and data obtained by comparing the three groups.

RESULTS

 We present 42 patients (25 PCrSD, 11 PCr24, 6 PCr48). PCr24 had a shorter procedure duration than PCrSD (65.8 vs. 91.2 minute; p = 0.048). Related to analgesia, PCr24 and PCr48 showed lower opioid consumption than PCrSD in PCA volume (48.5 and 49.6 vs. 75.1 mL; p = 0.015) and PCA time (23.3 and 23.8 vs. 34.3 hours; p = 0.01). Degree of pain (VAS scale) on the day of surgery and on the second postoperative day was lower in PCr24 and PCr48 than in PCrSD (4 and 2 vs. 5; p = 0.012; 0 and 1 vs. 2; p = 0.01, respectively) as well as shorter hospital stay (3 and 3.5 vs. 5 days; p = 0.021). In addition, PCr24 showed lower opioid consumption and hospital stay than PCr48 (p > 0.05). The greatest savings in hospital costs were obtained in the PCr24 group.

CONCLUSION

 PCr48 and PCr24 prior to PE surgery offers lower opioid consumption, less pain and shorter hospital stay than PCrSD. PCr24 is comparable to PCr48, but seems to show advantages and simpler logistics for the patient and the hospital.

Quest for the Most Aromatic Pathway in Charged Expanded Porphyrins

Despite the central role of aromaticity in the chemistry of expanded porphyrins, the evaluation of aromaticity remains difficult for these extended macrocycles. The presence of multiple conjugation pathways and different planar and nonplanar π-conjugation topologies makes the quantification of global and local aromaticity even more challenging. In neutral expanded porphyrins, the predominance of the aromatic conjugation pathway passing through the imine-type nitrogens and circumventing the amino NH groups is established. However, for charged macrocycles, the question about the main conjugation circuit remains open. Accordingly, different conjugation pathways in a set of neutral, anionic, and cationic expanded porphyrins were investigated by means of several aromaticity indices rooted in the structural, magnetic, and electronic criteria. Overall, our results reveal the predominance of the conjugation pathway that passes through all nitrogen atoms to describe the aromaticity of deprotonated expanded porphyrins, while the outer pathway through the perimeter carbon atoms becomes the most aromatic in protonated macrocycles. In nonplanar and charged macrocycles, a discrepancy between electronic and magnetic descriptors is observed. Nevertheless, our work demonstrates AV_min remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Searching for a prognostic index in lupus nephritis

BACKGROUND

Currently we do not have an ideal biomarker in lupus nephritis (LN) that should help us to identify those patients with SLE at risk of developing LN or to determine those patients at risk of renal progression. We aimed to evaluate the development of a prognostic index for LN, through the evaluation of clinical, analytical and histological factors used in a cohort of lupus. We have proposed to determine which factors, 6 months after the diagnosis of LN, could help us to define which patients will have a worse evolution of the disease and may be, more aggressive treatment and closer follow-up.

METHODS

A retrospective study to identify prognostic factors was carried out. We have included patients over 18 years of age with a clinical diagnosis of systemic lupus erythematosus (SLE) and kidney involvement confirmed by biopsy, who are followed up in our centre during the last 20 years. A multi-step statistical approach will be used in order to obtain a limited set of parameters, optimally selected and weighted, that show a satisfactory ability to discriminate between patients with different levels of prognosis.

RESULTS

We analysed 92 patients with LN, although only 73 have been able to be classified according to whether or not they have presented poor renal evolution. The age of onset (44 vs. 32; p = 0.024), the value of serum creatinine (1.41 vs. 1.04; p = 0.041), greater frequency of thrombocytopenia (30 vs. 7%; p = 0.038), higher score in the renal chronicity index (2.47 vs. 1.04; p = 0.015), proliferative histological type (100%) and higher frequency of interstitial fibrosis (67 vs. 32%; p = 0.017) and tubular atrophy (67 vs. 32%; p = 0.018) was observed between two groups. The multivariate analysis allowed us to select the best predictive model for poor outcome at 6 months based on different adjustment and discrimination parameters.

CONCLUSION

We have developed a prognostic index of poor renal evolution in patients with LN that combines demographic, clinical, analytical and histopathological factors, easy to use in routine clinical practice and that could be an effective tool in the early detection and management.

Does Supramolecular Gelation Require an External Trigger?

The supramolecular gelation of small molecules is typically preceded by an external stimulus to trigger the self-assembly. The need for this trigger stems from the metastable nature of most supramolecular gels and can limit their applicability. Herein, we present a small urea-based molecule that spontaneously forms a stable hydrogel by simple mixing without the addition of an external trigger. Single particle tracking experiments and observations made from scanning electron microscopy indicated that triggerless gelation occurred in a similar fashion as the archetypical heat-triggered gelation. These results could stimulate the search for other supramolecular hydrogels that can be obtained by simple mixing. Furthermore, the mechanism of the heat-triggered supramolecular gelation was elucidated by a combination of molecular dynamics simulations and quantitative NMR experiments. Surprisingly, hydrogelation seemingly occurs via a stepwise self-assembly in which spherical nanoparticles mature into an entangled fibrillary network.

Modulation of Metal Carbonyl Stretching Frequencies in the Second Coordination Sphere through the Internal Stark Effect

Spectroscopic and computational examination of a homologous series of rhodium(I) pybox carbonyl complexes has revealed a correlation between the conformation of the flanking aryl-substituted oxazoline donors and the carbonyl stretching frequency. This relationship is also observed experimentally for octahedral rhodium(III) and ruthenium(II) variants and cannot be explained through the classical, Dewar-Chatt-Duncanson, interpretation of metal-carbonyl bonding. Instead, these findings are reconciled by local changes in the magnitude of the electric field that is projected along the metal-carbonyl vector: the internal Stark effect.

Defining patients needs and expectations for eHealth-based cardiac rehabilitation in Germany and Spain: living lab data from the TIMELY study

Background

eHealth innovations have set the stage to optimize personalized care and provide assistance for disease control to patients with CAD. TIMELY is the first artificial intelligence (AI)-driven eHealth approach that employs internet of things devices and is based on cardiac rehabilitation (CR) components. To guarantee acceptance and usability of eHealth solutions, patients are actively participating in development through a Living Lab approach.

Purpose

To define patients' needs for an eHealth-based lifestyle intervention and self-care support.

Methods

The Living Lab approach included a guided survey conducted among CAD patients at CR centers in Germany and Spain during inpatient or outpatient CR, respectively. Questions referred to current use of technology and patients' opinion on the usefulness of suggested features of a future eHealth application. Ratings for usefulness/importance were recorded on a 5-point Likert scale reported as median score.

Results

79 patients (20% female) were interviewed (DE, n=49; ES, n=30). Patients' mean age was 57 years (range 37–79), educational level was 87% ≤ high school and 13% > high school. All patients owned a smartphone that they also used for information (76%) and documentation (43%). Patients rated the importance of all CR components (regular exercise, healthy diet, stress management, smoking cessation, risk factor reduction) along the continuum of care as “very high” (5/5). Individual need for regular exercise support after structured CR was rated “high” (4/5). Exercise reminders, suggestions on activities, update of recommended training heart rate, evaluation of training progress, and achieved goals were rated “useful” (4/5) to “very useful” (5/5). The importance of support for diet, stress management, overall risk factor management, and medication was rated 3/5. Usefulness of a learning/education tool was rated 4/5, while motivational messages scored 3/5 and individual feedback of a person or virtual agent scored 4/5 and 3/5, respectively. The availability of electronic health records was evaluated as “very useful” (5/5) remote ECG monitoring and blood pressure management were rated as “useful” (4/5). The range of all items accessed was 1–5, indicating large interindividual differences. No significant differences existed between female and male patients or German and Spanish patients, though the importance of sharing training progress with family/ friends was rated higher by men (3/5) compared to women (1/5) and higher by Spanish (4/5) compared to German patients (2/5).

Conclusions

CR patients expressed a strong need for individual support of regular physical exercise and updated training recommendations. The observed inter-individual differences regarding usefulness and acceptance indicate the need for a highly adaptive system to prevent exclusion from eHealth access. CR patients from Germany and Spain showed equally high acceptance of eHealth components for the management of CAD.

Funding Acknowledgement

Type of funding sources: Public grant(s) – EU funding. Main funding source(s): European Commission Horizon 2020

Correction: SuFEx-enabled, chemoselective synthesis of triflates, triflamides and triflimidates

Correction for ‘SuFEx-enabled, chemoselective synthesis of triflates, triflamides and triflimidates’ by Bing-Yu Li et al., Chem. Sci., 2022, 13, 2270–2279, DOI: 10.1039/D1SC06267K.

Percutaneous Cryoanalgesia: A New Strategy for Pain Management in Pectus Excavatum Surgery

INTRODUCTION

 In recent years, pain protocols for pectus excavatum (PE) have incorporated cryoanalgesia through thoracoscopic approach. Since 2019, ultrasound-guided percutaneous cryoanalgesia (PCr) has been applied at our institution, either on the same day as the Nuss procedure or 48 hours before surgery. We carried out a preliminary retrospective review of patients with PE in whom PCr prior to surgery was performed at our institution between 2019 and 2021.

MATERIALS AND METHODS

 Two groups were evaluated: PCr on the same day (PCrSD) and PCr 48 hours before surgery (PCr48). Despite PCr, patients were treated with "patient-controlled analgesia" (PCA) with opioids for at least 24 hours, switching to conventional intravenous analgesia and oral analgesia in the following days. Demographic, clinical-radiological variables, PCA opioid use, pain grade according to the visual analog scale (VAS), and length of stay (LOS) were compared between the groups. A total of 20 patients were included (12 with PCrSD and 8 with PCr48), without significant differences in demographics or clinical-radiological variables. The overall median time of PCr was 65 minutes (55-127), with no differences between the groups.

RESULTS

 PCr48 group presented with significantly lower median number of hours of continuous PCA (24 vs. 32 hours; p = 0.031), lower median number of rescue boluses (11 vs. 18; p = 0.042), lower median VAS in the early postoperative hours (2 vs. 5.5; p = 0.043), and lower median LOS (3.5 vs. 5 days).

CONCLUSION

 PCr performed 48 hours prior to surgery is more effective in terms of PCA requirements, VAS, and LOS when compared with cryoanalgesia on the same day.

Conceptual Density Functional Theory under Pressure: Part I. XP-PCM Method Applied to Atoms

High pressure chemistry offers the chemical community a range of possibilities to control chemical reactivity, develop new materials and fine-tune chemical properties. Despite the large changes that extreme pressure brings to the table, the field has mainly been restricted to the effects of volume changes and thermodynamics with less attention devoted to electronic effects at the molecular scale. This paper combines the conceptual DFT framework for analyzing chemical reactivity with the XP-PCM method for simulating pressures in the GPa range. Starting from the new derivatives of the energy with respect to external pressure, an electronic atomic volume and an atomic compressibility are found, comparable to their enthalpy analogues, respectively. The corresponding radii correlate well with major known sets of this quantity. The ionization potential and electron affinity are both found to decrease with pressure using two different methods. For the electronegativity and chemical hardness, a decreasing and increasing trend is obtained, respectively, and an electronic volume-based argument is proposed to rationalize the observed periodic trends. The cube of the softness is found to correlate well with the polarizability, both decreasing under pressure, while the interpretation of the electrophilicity becomes ambiguous at extreme pressures. Regarding the electron density, the radial distribution function shows a clear concentration of the electron density towards the inner region of the atom and periodic trends can be found in the density using the Carbó quantum similarity index and the Kullback–Leibler information deficiency. Overall, the extension of the CDFT framework with pressure yields clear periodic patterns.

Conceptual DFT has provided a framework in which to study chemical reactivity. Since high pressure is more and more a tool to control reactions and fine-tune chemical properties, this variable is introduced into the CDFT framework.

SuFEx-enabled, chemoselective synthesis of triflates, triflamides and triflimidates

Sulfur(vi) Fluoride Exchange (SuFEx) chemistry has emerged as a next-generation click reaction, designed to assemble functional molecules quickly and modularly. Here, we report the ex situ generation of trifluoromethanesulfonyl fluoride (CF₃SO₂F) gas in a two chamber system, and its use as a new SuFEx handle to efficiently synthesize triflates and triflamides. This broadly tolerated protocol lends itself to peptide modification or to telescoping into coupling reactions. Moreover, redesigning the S^VI–F connector with a S Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O → SNR replacement furnished the analogous triflimidoyl fluorides as SuFEx electrophiles, which were engaged in the synthesis of rarely reported triflimidate esters. Notably, experiments showed H₂O to be the key towards achieving chemoselective trifluoromethanesulfonation of phenols vs. amine groups, a phenomenon best explained—using ab initio metadynamics simulations—by a hydrogen bonded termolecular transition state for the CF₃SO₂F triflylation of amines.

Triflyl fluoride gas (CF₃SO₂F) and its aza analogues are reported as new SuFEx activators. These S^VI–F reagents react efficiently with a variety of nucleophiles, yet the presence of water grants complete chemoselectivity to phenols.

Fine-Tuning of Nonlinear Optical Contrasts of Hexaphyrin-Based Molecular Switches Using Inverse Design

In the search for new nonlinear optical (NLO) switching devices, expanded porphyrins have emerged as ideal candidates thanks to their tunable chemical and photophysical properties. Introducing meso-substituents to these macrocycles is a successful strategy to enhance the NLO contrasts. Despite its potential, the influence of meso-substitution on their structural and geometrical properties has been scarcely investigated. In this work, we pursue to grasp the underlying pivotal concepts for the fine-tuning of the NLO contrasts of hexaphyrin-based molecular switches, with a particular focus on the first hyperpolarizability related to the hyper-Rayleigh scattering (β_HRS). Building further on these concepts, we also aim to develop a rational design protocol. Starting from the (un)substituted hexaphyrins with various π-conjugation topologies and redox states, structure-property relationships are established linking aromaticity, photophysical properties and β_HRS responses. Ultimately, inverse molecular design using the best-first search algorithm is applied on the most favorable switches with the aim to further explore the combinatorial chemical compound space of meso-substituted hexaphyrins in search of high-contrast NLO switches. Two definitions of the figure-of-merit of the switch performance were used as target objectives in the optimization problem. Several meso-substitution patterns and their underlying characteristics are identified, uncovering molecular symmetry and the electronic nature of the substituents as the key players for fine-tuning the β_HRS values and NLO contrasts of hexaphyrin-based switches.

Uncovering the Activity of Alkaline Earth Metal Hydrogenation Catalysis Through Molecular Volcano Plots

Recent advances in alkaline earth (Ae) metal hydrogenation catalysis have broadened the spectrum of potential catalysts to include candidates from the main group, providing a sustainable alternative to the commonly used transition metals. Although Ae-amides have already been demonstrated to catalyze hydrogenation of imines and alkenes, a lucid understanding of how different metal/ligand combinations influence the catalytic activity is yet to be established. In this article, we use linear scaling relationships and molecular volcano plots to assess the potential of the Ae metal-based catalysts for the hydrogenation of alkenes. By analyzing combinations of eight metals (mono-, bi-, tri-, and tetravalent) and seven ligands, we delineate the impact of metal-ligand interplay on the hydrogenation activity. Our findings highlight that the catalytic activity is majorly determined by the charge and the size of the metal ions. While bivalent Ae metal cations delicately regulate the binding and the release of the reactants and the products, respectively, providing the right balance for this reaction, ligands play only a minor role in determining their catalytic activity. We show how volcano plots can be utilized for the rapid screening of prospective Ae catalysts to establish a guideline to achieve maximum activity in facilitating the hydrogenation process.

Computational Tools to Rationalize and Predict the Self-Assembly Behavior of Supramolecular Gels

Supramolecular gels form a class of soft materials that has been heavily explored by the chemical community in the past 20 years. While a multitude of experimental techniques has demonstrated its usefulness when characterizing these materials, the potential value of computational techniques has received much less attention. This review aims to provide a complete overview of studies that employ computational tools to obtain a better fundamental understanding of the self-assembly behavior of supramolecular gels or to accelerate their development by means of prediction. As such, we hope to stimulate researchers to consider using computational tools when investigating these intriguing materials. In the concluding remarks, we address future challenges faced by the field and formulate our vision on how computational methods could help overcoming them.

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.

Reactivity of Single Transition Metal Atoms on a Hydroxylated Amorphous Silica Surface: A Periodic Conceptual DFT Investigation

The drive to develop maximal atom-efficient catalysts coupled to the continuous striving for more sustainable reactions has led to an ever-increasing interest in single-atom catalysis. Based on a periodic conceptual density functional theory (cDFT) approach, fundamental insights into the reactivity and adsorption of single late transition metal atoms supported on a fully hydroxylated amorphous silica surface have been acquired. In particular, this investigation revealed that the influence of van der Waals dispersion forces is especially significant for a silver (98 %) or gold (78 %) atom, whereas the oxophilicity of the Group 8-10 transition metals plays a major role in the interaction strength of these atoms on the irreducible SiO₂ support. The adsorption energies for the less-electronegative row 4 elements (Fe, Co, Ni) ranged from -1.40 to -1.92 eV, whereas for the heavier row 5 and 6 metals, with the exception of Pd, these values are between -2.20 and -2.92 eV. The deviating behavior of Pd can be attributed to a fully filled d-shell and, hence, the absence of the hybridization effects. Through a systematic analysis of cDFT descriptors determined by using three different theoretical schemes, the Fermi weighted density of states approach was identified as the most suitable for describing the reactivity of the studied systems. The main advantage of this scheme is the fact that it is not influenced by fictitious Coulomb interactions between successive, charged reciprocal cells. Moreover, the contribution of the energy levels to the reactivity is simultaneously scaled based on their position relative to the Fermi level. Finally, the obtained Fermi weighted density of states reactivity trends show a good agreement with the chemical characteristics of the investigated metal atoms as well as the experimental data.

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.

Single-molecule conductance in a unique cross-conjugated tetra(aminoaryl)ethene

A 1,1,2,2-tetrakis(4-aminophenyl)ethene with three paths of π-conjugation, linear-cis, linear-trans and a cross-conjugation, has been prepared. The molecule is able to bind to gold electrodes forming molecular junctions for single-molecule conductance measurements. Only two regimes of conduction are found experimentally. The modelling of the conductance allows to assign them to through-bond transmission in the linear case, while the cross-conjugated channel is further assisted by through-space transmission, partially alleviating the destructive quantum interference.

Synthesis and Reactivity of Spirocarbocycles as Scaffolds for Nucleoside Analogues

A novel class of substituted spiro[3.4]octanes can be accessed via a [2 + 2]-cycloaddition of dichloroketene on a readily prepared exo-methylene cyclopentane building block. This reaction sequence was found to be robust on a multigram scale and afforded a central spirocyclobutanone scaffold for carbocyclic nucleosides. The reactivity of this constrained building block was evaluated and compared to the corresponding 4'-spirocyclic furanose analogues. Density functional theory calculations were performed to support the observed selectivity in the carbonyl reduction of spirocyclobutanone building blocks. Starting from novel spirocyclic intermediates, we exemplified the preparation of an undescribed class of carbocyclic nucleoside analogues and provided a proof of concept for application as inhibitors for the protein methyltransferase target PRMT5.

Mutations in the voltage-gated sodium channel gene associated with deltamethrin resistance in commercially sourced Phytoseiulus persimilis

The implementation of Integrated Pest Management in current agricultural practice is a convenient and very effective strategy to keep pest populations under control. The use of biological control agents, such as Phytoseiulus persimilis, is key for the success of such an approach. This predatory mite is widely used as it is very effective for controlling Tetranychus urticae, one of the most devastating crop pests. Here, we identify several mutations located in the voltage-gated sodium channel (VGSC) of commercially sourced P. persimilis that correlate with a reduced susceptibility to the pyrethroid deltamethrin. We found that the mites sourced from two different biocontrol product companies have intrinsic genotypic differences that correlate with their phenotype when tested with different concentrations of deltamethrin. Mites from Syngenta Bioline, carrying the mutations M918L and A1536T, were able to survive deltamethrin concentrations of up to 10 ppm, while the mites from Koppert Biological Systems, with the combination M918L, L925V and S1539T, survived treatment with 40 ppm. All of the point mutations identified in the predatory mite samples are located in a particular region of the VGSC, previously proposed as the binding site for this family of pesticides and identified as a 'hot spot' for resistance.

Molecular epidemiology of an enterovirus A71 outbreak associated with severe neurological disease, Spain, 2016

Introduction

Enterovirus A71 (EV-A71) is an emerging pathogen that causes a wide range of disorders including severe neurological manifestations. In the past 20 years, this virus has been associated with large outbreaks of hand, foot and mouth disease with neurological complications in the Asia-Pacific region, while in Europe mainly sporadic cases have been reported. In spring 2016, however, an EV-A71 outbreak associated with severe neurological cases was reported in Catalonia and spread further to other Spanish regions.

Aim

Our objective was to investigate the epidemiology and clinical characteristics of the outbreak.

Methods

We carried out a retrospective study which included 233 EV-A71-positive samples collected during 2016 from hospitalised patients. We analysed the clinical manifestations associated with EV-A71 infections and performed phylogenetic analyses of the 3’-VP1 and 3Dpol regions from all Spanish strains and a set of EV-A71 from other countries.

Results

Most EV-A71 infections were reported in children (mean age: 2.6 years) and the highest incidence was between May and July 2016 (83%). Most isolates (218/233) were classified as subgenogroup C1 and 217 of them were grouped in one cluster phylogenetically related to a new recombinant variant strain associated with severe neurological diseases in Germany and France in 2015 and 2016. Moreover, we found a clear association of EV-A71-C1 infection with severe neurological disorders, brainstem encephalitis being the most commonly reported.

Conclusion

An emerging recombinant variant of EV-A71-C1 was responsible for the large outbreak in 2016 in Spain that was associated with many severe neurological cases.

Stereoselective Reductions of 3-Substituted Cyclobutanones: A Comparison between Experiment and Theory

The stereoselective reduction of carbonyls is of key importance in the total synthesis of natural products and in medicinal chemistry. Nevertheless, models for rationalizing the stereoselectivity of the hydride reductions of cyclobutanones toward cyclobutanols are largely lacking, unlike cyclohexanone reductions. In order to elucidate the factors that control the stereoselectivity of these reductions, we have investigated the effect of the reaction temperature, solvent, substituent, and type of reducing agent using a synergistic experimental-computational approach. On the experimental side, the hydride reduction of 3-substituted cyclobutanones was proven to be highly selective for the formation of cis alcohol (>90%), irrespective of the size of the hydride reagent. The pronounced selectivity can be further enhanced by lowering the reaction temperature or decreasing the solvent polarity. On the computational side, density functional theory and noncovalent interaction analysis reveal that torsional strain plays a major role in the preference for the antifacial hydride approach, consistent with the Felkin-Anh model. In the presence of the benzyloxy substituent, the high selectivity for the cis isomer is also driven by repulsive electrostatic interactions in the case of a syn-facial hydride attack. The computed cis/trans ratios are in good agreement with the experimental ones and thus show the potential of computational chemistry for predicting and rationalizing the stereoselectivity of hydride reductions of cyclobutanones.

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Hückel-Möbius Interconversions in Expanded Porphyrins

Localized orbital coupled cluster theory has recently emerged as a nonempirical alternative to DFT for large systems. Intuitively, one might expect such methods to perform less well for highly delocalized systems. In the present work, we apply both canonical CCSD(T) approximations and a variety of localized approximations to a set of flexible expanded porphyrins—macrocycles that can switch between Hückel, figure-eight, and Möbius topologies under external stimuli. Both minima and isomerization transition states are considered. We find that Möbius(-like) structures have much stronger static correlation character than the remaining structures, and that this causes significant errors in DLPNO-CCSD(T) and even DLPNO-CCSD(T₁) approaches, unless TightPNO cutoffs are employed. If sub-kcal mol^–1 accuracy with respect to canonical relative energies is required even for Möbius-type systems (or other systems plagued by strong static correlation), then Nagy and Kallay’s LNO-CCSD(T) method with “tight” settings is the suitable localized approach. We propose the present POLYPYR21 data set as a benchmark for localized orbital methods or, more broadly, for the ability of lower-level methods to handle energetics with strongly varying degrees of static correlation.

Hand perfusion following radial or ulnar forearm free flap harvest for oral cavity reconstruction: A prospective study

The radial forearm free flap (RFFF) and ulnar forearm free flap (UFFF) are used in head and neck reconstruction because they provide a thin and pliable skin paddle as well as a long vascular pedicle. However, in spite of several studies showing the safety of the UFFF, the RFFF is more popular among reconstructive surgeons based on concerns about hand ischaemia. A prospective study was designed in which 10 UFFF and 11 RFFF surgeries were performed in 20 patients undergoing oral cavity reconstruction between January 2017 and July 2018. Hand vascular parameters were evaluated preoperatively and postoperatively using Doppler ultrasound and plethysmography. The preoperative and postoperative diameters of the radial and ulnar arteries, and the flow velocities through the remainder of the forearm artery were measured preoperatively and at 3 months postoperative. Additionally, a comparison was performed between the preoperative and postoperative fingertip perfusion values according to impedance plethysmography. The preoperative mean diameter of the radial artery (2.89±0.47mm) was significantly greater than that of the ulnar artery (2.35±0.48mm) at the level of the wrist; however, 3 months after the surgery, the mean diameters of the two arteries did not differ significantly. There were no differences in digital perfusion when a UFFF was used compared with an RFFF.

Performance of Electronic Structure Methods for the Description of Hückel-Möbius Interconversions in Extended π-Systems

Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. DFT remains the workhorse for modeling such extended macrocycles, when taking into account their size and huge conformational flexibility. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, π···π stacking, steric effects, ring strain, and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wave function methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a stronger degree of static correlation than the Hückel and twisted-Hückel structures, and as a result the relative energies of singly twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between the extrapolated MP2/CBS energies estimated from cc-pV{T,Q}Z, diffuse-augmented heavy-aug-cc-pV{T,Q}Z and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n-particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol-1 with the CCSD(T) relative energies. Regarding DFT methods, the range-separated double hybrids, such as ωB97M(2) and B2GP-PLYP, outperform other functionals with RMSDs of 0.6 and 0.8 kcal mol-1, respectively. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended π-systems, the errors drop down to 1.9 kcal mol-1 for the revised revDOD-PBEP86-NL, which eliminates the same-spin correlation energy. Minnesota meta-GGA functionals with high fractions of exact exchange (M06-2X and M08-HX) also perform reasonably well, outperforming more robust and significantly less empirically parametrized functionals like SCAN0-D3.

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.

Factors associated with sepsis development in 606 Spanish adult patients with cellulitis

Background

Cellulitis, a frequent cause of admission of adult patients to medical wards, occasionally evolves to sepsis. In this study we analyze the factors related to sepsis development.

Methods

Prospective and observational study of 606 adult patients with cellulitis admitted to several Spanish hospitals. Comorbidities, microbiological, clinical, lab, diagnostic, and treatment data were analyzed. Sepsis was diagnosed according to the criteria of the 2016 International Sepsis Definitions Conference. Multiple logistic regression modelling was performed to determine the variables independently associated with sepsis development.

Results

Mean age was 63.4 years and 51.8% were men. Overall 65 (10.7%) patients developed sepsis, 7 (10.8%) of whom died, but only 4 (6.2%) due to cellulitis. Drawing of blood (P < 0.0001) or any (P < 0.0001) culture, and identification of the agent (P = 0.005) were more likely among patients with sepsis. These patients had also a longer duration of symptoms (P = 0.04), higher temperature (P = 0.03), more extensive cellulitis (P = 0.02), higher leukocyte (P < 0.0001) and neutrophil (P < 0.0001) counts, serum creatinine (P = 0.001), and CRP (P = 0.008) than patients without sepsis. Regarding therapy, patients with sepsis were more likely to undergo changes in the initial antimicrobial regimen (P < 0.0001), received more antimicrobials (P < 0.0001), received longer intravenous treatment (P = 0.03), and underwent surgery more commonly (P = 0.01) than patients without sepsis. Leukocyte counts (P = 0.002), serum creatinine (P = 0.003), drawing of blood cultures (P = 0.004), change of the initial antimicrobial regimen (P = 0.007) and length of cellulitis (P = 0.009) were independently associated with sepsis development in the multivariate analysis.

Conclusions

Increased blood leukocytes and serum creatinine, blood culture drawn, modification of the initial antimicrobial regimen, and maximum length of cellulitis were associated with sepsis in these patients.

Switching between Hückel and Möbius aromaticity: a density functional theory and information-theoretic approach study

Benziporphyrins are versatile macrocycles exhibiting aromaticity switching behaviors. The existence of both Hückel and Möbius (anti)aromaticity has been reported in these systems, whose validity is respectively governed by the [4n + 2] and [4n] π-electron rule on the macrocyclic pathway. Despite the experimental evidence on the floppiness of benziporphyrins, the switching mechanism between Hückel and Möbius structures is still not clear, as well as the factors influencing the stability of the different π-conjugation topologies. For these reasons, we performed a systematic study on A,D-di-p-benzihexaphyrins(1.1.1.1.1.1) with two redox states corresponding to [28] and [30] π-electron conjugation pathways. Whereas benzi[28]hexaphyrin obeys Möbius aromaticity, benzi[30]hexaphyrin follows Hückel aromaticity. The dynamic interconversion between Möbius and Hückel aromaticity is investigated through the rotation of a phenylene ring, which acts as the topology selector. Further analyses of the energy profiles using energy decomposition and information-theoretic approaches provide new insights into conformational stability, aromaticity and antiaromaticity for these species. Strong and opposite cross correlations between aromaticity indexes and information-theoretic quantities were found for the two macrocyclic systems with opposite global aromaticity and antiaromaticity behaviors. These results indicate that Hückel and Möbius aromaticity and antiaromaticity, though qualitatively different, are closely related and can be interchanged, and information-theoretic quantities provide a novel understanding about their relevance. Our present results should provide in-depth insights to appreciate the nature and origin about Möbius (anti)aromaticity and its close relationship with Hückel (anti)aromaticity.

Rationalising Supramolecular Hydrogelation of Bis-Urea-Based Gelators through a Multiscale Approach

Invited for this month's cover is the Quantum Chemistry group (ALGC) of the Free University of Brussels (VUB) in collaboration with the Molecular Design and Synthesis (MolDesignS) group of the KU Leuven. The image shows how a variety of computational techniques, such as density functional theory, molecular dynamics, and the noncovalent interaction index, can be used to virtually zoom in on the noncovalent interactions that are key for the supramolecular hydrogelation of bis-urea-based gelators. Read the full text of the article at 10.1002/cplu.201900551.

Molecular dynamics based descriptors for predicting supramolecular gelation

Four molecular dynamics-based descriptors were derived able to classify gelator–solvent combinations as a gel, precipitate or clear solution.

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.