Search

Quantum chemical calculations were carried out to quantitatively understand the origin of the Felkin-Anh(-Eisenstein) model, widely used to rationalize the π-facial stereoselectivity in the nucleophilic addition reaction to carbonyl groups directly attached to a stereogenic center. To this end, the possible approaches of cyanide to both ( S )-2-phenylpropanal and ( S )-3-phenylbutan-2-one have been explored in detail. With the help of the activation strain model of reactivity and the energy decomposition analysis method, it is found that the preference for the Felkin-Anh addition is mainly dictated by steric factors which manifest in a less destabilizing strain-energy rather than, as traditionally considered, in a lower Pauli repulsion. In addition, other factors such as the more favorable electrostatic interactions also contribute to the preferred approach of the nucleophile. Our work, therefore, provides a different, more complete rationalization, based on quantitative analyses, of the origin of this seminal and highly useful concept in organic chemistry., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

A common goal in chemistry is to study the relationship between a measured signal and the variability of certain factors. To this end, researchers often use Design of Experiment to decide which experiments to conduct and (Multiple) Linear Regression, and/or Analysis of Variance to analyze the collected data. Among the assumptions to the very foundation of this strategy, all the experiments are independent, conditional on the settings of the factors. Unfortunately, due to the presence of uncontrollable factors, real-life experiments often deviate from this assumption, making the data analysis results unreliable. In these cases, Mixed-Effects modeling, despite not being widely used in chemometrics, represents a solid data analysis framework to obtain reliable results. Here we provide a tutorial for Linear Mixed-Effects models. We gently introduce the reader to these models by showing some motivating examples. Then, we discuss the theory behind Linear Mixed-Effect models, and we show how to fit these models by making use of real-life data obtained from an exposome study. Throughout the paper we provide R code so that each researcher is able to implement these useful model themselves., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

We have quantum chemically investigated how microsolvation affects the various E2 and S N 2 pathways, their mutual competition, and the α-effect of the model reaction system HOO - (H 2 O) n + CH 3 CH 2 Cl, at the CCSD(T) level. Interestingly, we identify the dual nature of the α-nucleophile HOO - which, upon solvation, is in equilibrium with HO - . This solvent-induced dual appearance gives rise to a rich network of competing reaction channels. Among both nucleophiles, S N 2 is always favored over E2, and this preference increases upon increasing microsolvation. Furthermore, we found a pronounced α-effect, not only for S N 2 substitution but also for E2 elimination, i.e. , HOO - is more reactive than HO - in both cases. Our activation strain and quantitative molecular orbital analyses reveal the physical mechanisms behind the various computed trends. In particular, we demonstrate that two recently proposed criteria, required for solvent-free nucleophiles to display the α-effect, must also be satisfied by microsolvated HOO - (H 2 O) n nucleophiles.

Although in vivo extracellular microenvironments are dynamic, most in vitro studies are conducted under static conditions. Here, we exposed diffuse large B-cell lymphoma (DLBCL) cells to gradient increases in the concentration of hydrogen peroxide (H 2 O 2 ), thereby capturing some of the dynamics of the tumour microenvironment. Subsequently, we measured the phosphorylation response of B-cell receptor (BCR) signalling proteins CD79a, SYK and PLCγ2 at a high temporal resolution via single-cell phospho-specific flow cytometry. We demonstrated that the cells respond bimodally to static extracellular H 2 O 2 , where the percentage of cells that respond is mainly determined by the concentration. Computational analysis revealed that the bimodality results from a combination of a steep dose-response relationship and cell-to-cell variability in the response threshold. Dynamic gradient inputs of varying durations indicated that the H 2 O 2 concentration is not the only determinant of the signalling response, as cells exposed to more shallow gradients respond at lower H 2 O 2 levels. A minimal model of the proximal BCR network qualitatively reproduced the experimental findings and uncovered a rate-dependent sensitivity to H 2 O 2 , where a lower rate of increase correlates to a higher sensitivity. These findings will bring us closer to understanding how cells process information from their complex and dynamic in vivo environments., (© 2024. The Author(s).)

The origin of the electrophilicity of a series of cyclohexanones and benzaldehydes is investigated using the activation strain model and quantitative Kohn-Sham molecular orbital (MO) theory. We find that this electrophilicity is mainly determined by the electrostatic attractions between the carbonyl compound and the nucleophile (cyanide) along the entire reaction coordinate. Donor-acceptor frontier molecular orbital interactions, on which the current rationale behind electrophilicity trends is based, appear to have little or no significant influence on the reactivity of these carbonyl compounds., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

The stereoselective introduction of glycosidic bonds (glycosylation) is one of the main challenges in the chemical synthesis of carbohydrates. Glycosylation reaction mechanisms are difficult to control because, in many cases, the exact reactive species driving product formation cannot be detected and the product outcome cannot be explained by the primary reaction intermediate observed. In these cases, reactions are expected to take place via other low-abundance reaction intermediates that are in rapid equilibrium with the primary reaction intermediate via a Curtin-Hammett scenario. Despite this principle being well-known in organic synthesis, mechanistic studies investigating this model in glycosylation reactions are complicated by the challenge of detecting the extremely short-lived reactive species responsible for product formation. Herein, we report the utilization of the chemical equilibrium between low-abundance reaction intermediates and the stable, readily observed α-glycosyl triflate intermediate in order to infer the structure of the former species by employing exchange NMR. Using this technique, we enabled the detection of reaction intermediates such as β-glycosyl triflates and glycosyl dioxanium ions. This demonstrates the power of exchange NMR to unravel reaction mechanisms as we aim to build a catalog of kinetic parameters, allowing for the understanding and eventual prediction of glycosylation reactions.

Understanding the geometrical preferences in chemical reactions is crucial for advancing the field of organic chemistry and improving synthetic strategies. One such preference, the Bürgi-Dunitz angle, is central to nucleophilic addition reactions involving carbonyl groups. This study successfully employs a novel two-dimensional Distortion-Interaction/Activation-Strain Model in combination with a two-dimensional Energy Decomposition Analysis to investigate the origins of the Bürgi-Dunitz angle in the addition reaction of CN - to (CH 3 ) 2 C═O. We constructed a 2D potential energy surface defined by the distance between the nucleophile and carbonylic carbon atom and by the attack angle, followed by an in-depth exploration of energy components, including strain and interaction energy. Our analysis reveals that the Bürgi-Dunitz angle emerges from a delicate balance between two key factors: strain energy and interaction energy. High strain energy, as a result of the carbonyl compound distorting to avoid Pauli repulsion, is encountered at high angles, thus setting the upper bound. On the other hand, interaction energy is shaped by a dominant Pauli repulsion when the angles are lower. This work emphasizes the value of the 2D Energy Decomposition Analysis as a refined tool, offering both quantitative and qualitative insights into chemical reactivity and selectivity.

We have quantum chemically studied the influence of ring strain on the competition between the two mechanistically different S N 2 and E2 pathways using a series of archetypal ethers as substrate in combination with a diverse set of Lewis bases (F - , Cl - , Br - , HO - , H 3 CO - , HS - , H 3 CS - ), using relativistic density functional theory at ZORA-OLYP/QZ4P. The ring strain in the substrate is systematically increased on going from a model acyclic ether to a 6- to 5- to 4- to 3-membered ether ring. We have found that the activation energy of the S N 2 pathway sharply decreases when the ring strain of the system is increased, thus on going from large to small cyclic ethers, the S N 2 reactivity increases. In contrast, the activation energy of the E2 pathway generally rises along this same series, that is, from large to small cyclic ethers. The opposing reactivity trends induce a mechanistic switch in the preferred reaction pathway for strong Lewis bases from E2, for large cyclic substrates, to S N 2, for small cyclic substrates. Weak Lewis bases are unable to overcome the higher intrinsic distortivity of the E2 pathway and, therefore, always favor the less distortive S N 2 reaction., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

The Bürgi-Dunitz (BD) angle plays a pivotal role in organic chemistry to rationalize the nucleophilic addition to carbonyl groups. Yet, the origin of the obtuse trajectory of the nucleophile remains incompletely understood. Herein, we quantify the importance of the underlying physical factors quantum chemically. The obtuse BD angle appears to originate from the concerted action of a reduced Pauli repulsion between the nucleophile HOMO and carbonyl π bond, a more stabilizing HOMO-π*-LUMO(C=O) interaction, as well as a more favorable electrostatic attraction., (© 2023 The Authors. ChemPhysChem published by Wiley-VCH GmbH.)

Interactions between RNA-binding proteins and RNA molecules are at the center of multiple biological processes. Therefore, accurate characterization of the composition of ribonucleoprotein complexes (RNPs) is crucial. Ribonuclease (RNase) for mitochondrial RNA processing (MRP) and RNase P are highly similar RNPs that play distinct roles at the cellular level; as a consequence, the specific isolation of either of these complexes is essential to study their biochemical function. Since their protein components are nearly identical, purification of these endoribonucleases using protein-centric methods is not feasible. Here, we describe a procedure employing an optimized high-affinity streptavidin-binding RNA aptamer, termed S1m, to purify RNase MRP free of RNase P. This report details all steps from the RNA tagging to the characterization of the purified material. We show that using the S1m tag allows efficient isolation of active RNase MRP., Competing Interests: Competing interests The authors declare no conflict of interest., (Copyright © 2023 The Authors; exclusive licensee Bio-protocol LLC.)

We have quantum chemically studied the iron-mediated CX bond activation (X = H, Cl, CH 3 ) by d 8 -FeL 4 complexes using relativistic density functional theory at ZORA-OPBE/TZ2P. We find that by either modulating the electronic effects of a generic iron-catalyst by a set of ligands, that is, CO, BF, PH 3 , BN(CH 3 ) 2 , or by manipulating structural effects through the introduction of bidentate ligands, that is, PH 2 (CH 2 ) n PH 2 with n = 6-1, one can significantly decrease the reaction barrier for the CX bond activation. The combination of both tuning handles causes a decrease of the CH activation barrier from 10.4 to 4.6 kcal mol -1 . Our activation strain and Kohn-Sham molecular orbital analyses reveal that the electronic tuning works via optimizing the catalyst-substrate interaction by introducing a strong second backdonation interaction (i.e., "ligand-assisted" interaction), while the mechanism for structural tuning is mainly caused by the reduction of the required activation strain because of the pre-distortion of the catalyst. In all, we present design principles for iron-based catalysts that mimic the favorable behavior of their well-known palladium analogs in the bond-activation step of cross-coupling reactions., (© 2022 The Authors. Journal of Computational Chemistry published by Wiley Periodicals LLC.)

We have quantum chemically studied the base-catalyzed Diels-Alder (DA) reaction between 3-hydroxy-2-pyrone and N-methylmaleimide using dispersion-corrected density functional theory. The uncatalyzed reaction is slow and is preceded by the extrusion of CO 2 via a retro-DA reaction. Base catalysis, for example, by triethylamine, lowers the reaction barrier up to 10 kcal mol -1 , causing the reaction to proceed smoothly at low temperature, which quenches the expulsion of CO 2 , yielding efficient access to polyoxygenated natural compounds. Our activation strain analyses reveal that the base accelerates the DA reaction via two distinct electronic mechanisms: i) by the HOMO-raising effect, which enhances the normal electron demand orbital interaction; and ii) by donating charge into 3-hydroxy-2-pyrone which accumulates in its reactive region and promotes strongly stabilizing secondary electrostatic interactions with N-methylmaleimide., (© 2022 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Cardiac valves ensure unidirectional blood flow through the heart, and altering their function can result in heart failure. Flow sensing via wall shear stress and wall stretching through the action of mechanosensors can modulate cardiac valve formation. However, the identity and precise role of the key mechanosensors and their effectors remain mostly unknown. Here, we genetically dissect the role of Pkd1a and other mechanosensors in atrioventricular (AV) valve formation in zebrafish and identify a role for several pkd and piezo gene family members in this process. We show that Pkd1a, together with Pkd2, Pkd1l1, and Piezo2a, promotes AV valve elongation and cardiac morphogenesis. Mechanistically, Pkd1a, Pkd2, and Pkd1l1 all repress the expression of klf2a and klf2b, transcription factor genes implicated in AV valve development. Furthermore, we find that the calcium-dependent protein kinase Camk2g is required downstream of Pkd function to repress klf2a expression. Altogether, these data identify, and dissect the role of, several mechanosensors required for AV valve formation, thereby broadening our understanding of cardiac valvulogenesis., (© 2023. The Author(s).)

The traditional and widespread rationale behind the stability trend of alkyl-substituted carbocations is incomplete. Through state-of-the-art quantum chemical analyses, we quantitatively established a generally overlooked driving force behind the stability of carbocations, namely, that the parent substrates are substantially destabilized by the introduction of substituents, often playing a dominant role in solution. This stems from the repulsion between the substituents and the C-X bond.

The gauche conformer in 1-X,2-Y-disubstituted ethanes, that is, the staggered orientation in which X and Y are in closer contact, is only favored for relatively small substituents that do not give rise to large X···Y steric repulsion. For more diffuse substituents, weakly attractive orbital interactions between antiperiplanar bonds (i.e., hyperconjugation) cannot overrule the repulsive forces between X and Y. Our quantum chemical analyses of the rotational isomerism of XCH 2 CH 2 Y (X = F, OH; Y = I) at ZORA-BP86-D3(BJ)/QZ4P reveal that indeed the anti conformer is generally favored due to a less destabilizing I···F and I···O-H steric repulsion. The only case when the gauche conformer is preferred is when the hydroxyl hydrogen is oriented toward the iodine atom in the 2-iodoethanol. This is because of the significantly stabilizing covalent component of the I···H-O intramolecular hydrogen bond. Therefore, we show that strong intramolecular interactions can overcome the steric repulsion between bulky substituents in 1,2-disubstituted ethanes and cause the gauche effect. Our quantum chemical computations have guided nuclear magnetic resonance experiments that confirm the increase in the gauche population as X goes from F to OH.

The C-X bond activation (X = H, C) of a series of substituted C(n°)-H and C(n°)-C(m°) bonds with C(n°) and C(m°) = H 3 C- (methyl, 0°), CH 3 H 2 C- (primary, 1°), (CH 3 ) 2 HC- (secondary, 2°), (CH 3 ) 3 C- (tertiary, 3°) by palladium were investigated using relativistic dispersion-corrected density functional theory at ZORA-BLYP-D3(BJ)/TZ2P. The effect of the stepwise introduction of substituents was pinpointed at the C-X bond on the bond activation process. The C(n°)-X bonds become substantially weaker going from C(0°)-X, to C(1°)-X, to C(2°)-X, to C(3°)-X because of the increasing steric repulsion between the C(n°)- and X-group. Interestingly, this often does not lead to a lower barrier for the C(n°)-X bond activation. The C-H activation barrier, for example, decreases from C(0°)-X, to C(1°)-X, to C(2°)-X and then increases again for the very crowded C(3°)-X bond. For the more congested C-C bond, in contrast, the activation barrier always increases as the degree of substitution is increased. Our activation strain and matching energy decomposition analyses reveal that these differences in C-H and C-C bond activation can be traced back to the opposing interplay between steric repulsion across the C-X bond versus that between the catalyst and substrate., (© 2022 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Hierarchical, convergent ab initio benchmark computations were performed followed by a systematic analysis of DFT performance for five pericyclic reactions comprising Diels-Alder, 1,3-dipolar cycloaddition, electrocyclic rearrangement, sigmatropic rearrangement, and double group transfer prototypes. Focal point analyses (FPA) extrapolating to the ab initio limit were executed via explicit quantum chemical computations with electron correlation treatments through CCSDT(Q) and correlation-consistent Gaussian basis sets up to aug'-cc-pV5Z. Optimized geometric structures and vibrational frequencies of all stationary points were obtained at the CCSD(T)/cc-pVTZ level of theory. The FPA reaction barriers and energies exhibit convergence to within a few tenths of a kcal mol -1 . The FPA benchmarks were used to evaluate the performance of 60 density functionals (eight dispersion-corrected), covering the local-density approximation (LDA), generalized gradient approximations (GGAs), meta-GGAs, hybrids, meta-hybrids, double-hybrids, and range-separated hybrids. The meta-hybrid M06-2X functional provided the best overall performance [mean absolute error (MAE) of 1.1 kcal mol -1 ] followed closely by the double-hybrids B2K-PLYP, mPW2K-PLYP, and revDSD-PBEP86 [MAE of 1.4-1.5 kcal mol -1 ]. The regularly used GGA functional BP86 gave a higher MAE of 5.8 kcal mol -1 , but it qualitatively described the trends in reaction barriers and energies. Importantly, we established that accurate yet efficient meta-hybrid or double-hybrid DFT potential energy surfaces can be acquired based on geometries from the computationally efficient and robust BP86/DZP level.

We have quantum chemically explored the competition between the S N 2 and S N 2' pathways for X - + H 2 C═CHCH 2 Y (X, Y = F, Cl, Br, I) using a combined relativistic density functional theory and coupled-cluster theory approach. Bimolecular nucleophilic substitution reactions at allylic systems, i.e., C γ ═C β -C α -Y, bearing a leaving-group at the α-position, proceed either via a direct attack at the α-carbon (S N 2) or via an attack at the γ-carbon, involving a concerted allylic rearrangement (S N 2'), in both cases leading to the expulsion of the leaving-group. Herein, we provide a physically sound model to rationalize under which circumstances a nucleophile will follow either the aliphatic S N 2 or allylic S N 2' pathway. Our activation strain analyses expose the underlying physical factors that steer the S N 2/S N 2' competition and, again, demonstrate that the concepts of a reaction's "characteristic distortivity" and "transition state acidity" provide explanations and design tools for understanding and predicting reactivity trends in organic synthesis.

Magnetic skyrmions are quasiparticles with nontrivial topology, envisioned to play a key role in next-generation data technology while simultaneously attracting fundamental research interest due to their emerging topological charge. In chiral magnetic multilayers, current-generated spin-orbit torques or ultrafast laser excitation can be used to nucleate isolated skyrmions on a picosecond time scale. Both methods, however, produce randomly arranged skyrmions, which inherently limits the precision on the location at which the skyrmions are nucleated. Here, we show that nanopatterning of the anisotropy landscape with a He + -ion beam creates well-defined skyrmion nucleation sites, thereby transforming the skyrmion localization into a deterministic process. This approach allows control of individual skyrmion nucleation as well as guided skyrmion motion with nanometer-scale precision, which is pivotal for both future fundamental studies of skyrmion dynamics and applications.

We have studied the palladium-mediated activation of C(sp n )-X bonds (n = 1-3 and X = H, CH 3 , Cl) in archetypal model substrates H 3 C-CH 2 -X, H 2 C=CH-X and HC≡C-X by catalysts PdL n with L n = no ligand, Cl - , and (PH 3 ) 2 , using relativistic density functional theory at ZORA-BLYP/TZ2P. The oxidative addition barrier decreases along this series, even though the strength of the bonds increases going from C(sp 3 )-X, to C(sp 2 )-X, to C(sp)-X. Activation strain and matching energy decomposition analyses reveal that the decreased oxidative addition barrier going from sp 3 , to sp 2 , to sp, originates from a reduction in the destabilizing steric (Pauli) repulsion between catalyst and substrate. This is the direct consequence of the decreasing coordination number of the carbon atom in C(sp n )-X, which goes from four, to three, to two along this series. The associated net stabilization of the catalyst-substrate interaction dominates the trend in strain energy which indeed becomes more destabilizing along this same series as the bond becomes stronger from C(sp 3 )-X to C(sp)-X., (© 2021 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

It has been experimentally observed that water-ice-embedded polycyclic aromatic hydrocarbons (PAHs) form radical cations when exposed to vacuum UV irradiation, whereas ammonia-embedded PAHs lead to the formation of radical anions. In this study, we explain this phenomenon by investigating the fundamental electronic differences between water and ammonia, the implications of these differences on the PAH-water and PAH-ammonia interaction, and the possible ionization pathways in these complexes using density functional theory (DFT) computations. In the framework of the Kohn-Sham molecular orbital (MO) theory, we show that the ionic state of the PAH photoproducts results from the degree of occupied-occupied MO mixing between the PAHs and the matrix molecules. When interacting with the PAH, the lone pair-type highest occupied molecular orbital (HOMO) of water has poor orbital overlap and is too low in energy to mix with the filled π-orbitals of the PAH. As the lone-pair HOMO of ammonia is significantly higher in energy and has better overlap with filled π-orbitals of the PAH, the subsequent Pauli repulsion leads to mixed MOs with both PAH and ammonia character. By time-dependent DFT calculations, we demonstrate that the formation of mixed PAH-ammonia MOs opens alternative charge-transfer excitation pathways as now electronic density from ammonia can be transferred to unoccupied PAH levels, yielding anionic PAHs. As this pathway is much less available for water-embedded PAHs, charge transfer mainly occurs from localized PAH MOs to mixed PAH-water virtual levels, leading to cationic PAHs., Competing Interests: The authors declare no competing financial interest., (© 2022 The Authors. Published by American Chemical Society.)

We have quantum chemically investigated how solvation influences the competition between the S N 2 and E2 pathways of the model F - + C 2 H 5 Cl reaction. The system is solvated in a stepwise manner by going from the gas phase, then via microsolvation of one to three explicit solvent molecules, then last to bulk solvation using relativistic density functional theory at (COSMO)-ZORA-OLYP/QZ4P. We explain how and why the mechanistic pathway of the system shifts from E2 in the gas phase to S N 2 upon strong solvation of the Lewis base (i.e., nucleophile/protophile). The E2 pathway is preferred under weak solvation of the system by dichloromethane, whereas a switch in reactivity from E2 to S N 2 is observed under strong solvation by water. Our activation strain and Kohn-Sham molecular orbital analyses reveal that solvation of the Lewis base has a significant impact on the strength of the Lewis base. We show how strong solvation furnishes a weaker Lewis base that is unable to overcome the high characteristic distortivity associated with the E2 pathway, and thus the S N 2 pathway becomes viable.

RNase MRP is a ribonucleoprotein complex involved in the endoribonucleolytic cleavage of different RNAs. Mutations in the RNA component of the RNP are the cause of cartilage hair hypoplasia. Patients with cartilage hair hypoplasia are characterized by skeletal dysplasia. Biochemical purification of RNase MRP is desired to be able to study its biochemical function, composition and activity in both healthy and disease situations. Due to the high similarity with RNase P, a method to specifically isolate the RNase MRP complex is currently lacking. By fusing a streptavidin-binding RNA aptamer, the S1m-aptamer, to the RNase MRP RNA we have been able to compare the relative expression levels of wildtype and mutant MRP RNAs. Moreover, we were able to isolate active RNase MRP complexes. We observed that mutant MRP RNAs are expressed at lower levels and have lower catalytic activity compared to the wildtype RNA. The observation that a single nucleotide substitution at position 40 in the P3 domain but not in other domains of RNase MRP RNA severely reduced the binding of the Rpp25 protein subunit confirmed that the P3 region harbours the main binding site for this protein. Altogether, this study shows that the RNA aptamer tagging approach can be used to identify RNase MRP substrates, but also to study the effect of mutations on MRP RNA expression levels and RNase MRP composition and endoribonuclease activity.

The reactivity of phenothiazine (PS), phenoselenazine (PSE), and phenotellurazine (PTE) with different reactive oxygen species (ROS) has been studied using density functional theory (DFT) in combination with the QM-ORSA (Quantum Mechanics-based Test for Overall Free Radical Scavenging Activity) protocol for an accurate kinetic rate calculation. Four radical scavenging mechanisms have been screened, namely hydrogen atom transfer (HAT), radical adduct formation (RAF), single electron transfer (SET), and the direct oxidation of the chalcogen atom. The chosen ROS are HO . , HOO . , and CH 3 OO . . PS, PSE, and PTE exhibit an excellent antioxidant activity in water regardless of the ROS due to their characteristic diffusion-controlled regime processes. For the HO . radical, the primary active reaction mechanism is, for all antioxidants, RAF. But, for HOO . and CH 3 OO . , the dominant mechanism strongly depends on the antioxidant: HAT for PS and PSE, and SET for PTE. The scavenging efficiency decreases dramatically in lipid environment and remains only significant (via RAF) for the most reactive radical (HO . ). Therefore, PS, PSE, and PTE are excellent antioxidant molecules, especially in aqueous, physiological environments where they are active against a broad spectrum of harmful radicals. There is no advantage or significant difference in the scavenging efficiency when changing the chalcogen since the reactivity mainly derives from the amino hydrogen and the aromatic sites., (© 2021 The Authors. ChemMedChem published by Wiley-VCH GmbH.)

The concept of dipolar repulsion has been widely used to explain several phenomena in organic chemistry, including the conformational preferences of carbonyl compounds. This model, in which atoms and bonds are viewed as point charges and dipole moment vectors, respectively, is however oversimplified. To provide a causal model rooted in quantitative molecular orbital theory, we have analyzed the rotational isomerism of haloacetaldehydes OHC-CH 2 X (X = F, Cl, Br, I), using relativistic density functional theory. We have found that the overall trend in the rotational energy profiles is set by the combined effects of Pauli repulsion (introducing a barrier around gauche that separates minima at syn and anti ), orbital interactions (which can pull the anti minimum towards anticlinal to maximize hyperconjugation), and electrostatic interactions. Only for X = F, not for X = Cl-I, electrostatic interactions push the preference from syn to anti . Our bonding analyses show how this trend is related to the compact nature of F versus the more diffuse nature of the heavier halogens.

We investigate the propagation of magnons after ultrashort perturbations of the exchange interaction in the prototype two-dimensional Heisenberg antiferromagnet. Using the recently proposed neural quantum states, we predict highly anisotropic spreading in space constrained by the symmetry of the perturbation. Interestingly, the propagation speed at the shortest length scale and timescale is up to 40% higher than the highest magnon velocity. We argue that the enhancement stems from extraordinary strong magnon-magnon interactions, suggesting new avenues for manipulating information transfer on ultrashort length scales and timescales.

Brønsted acid-catalyzed inverse-electron demand (IED) aza-Diels-Alder reactions between 2-aza-dienes and ethylene were studied using quantum chemical calculations. The computed activation energy systematically decreases as the basic sites of the diene progressively become protonated. Our activation strain and Kohn-Sham molecular orbital analyses traced the origin of this enhanced reactivity to i) "Pauli-lowering catalysis" for mono-protonated 2-aza-dienes due to the induction of an asynchronous, but still concerted, reaction pathway that reduces the Pauli repulsion between the reactants; and ii) "LUMO-lowering catalysis" for multi-protonated 2-aza-dienes due to their highly stabilized LUMO(s) and more concerted synchronous reaction path that facilitates more efficient orbital overlaps in IED interactions. In all, we illustrate how the novel concept of "Pauli-lowering catalysis" can be overruled by the traditional concept of "LUMO-lowering catalysis" when the degree of LUMO stabilization is extreme as in the case of multi-protonated 2-aza-dienes., (© 2021 The Authors. Published by Wiley-VCH GmbH.)

Chemical reactions are ubiquitous in the universe, they are at the core of life, and they are essential for industrial processes. The drive for a deep understanding of how something occurs, in this case, the mechanism of a chemical reaction and the factors controlling its reactivity, is intrinsically valuable and an innate quality of humans. The level of insight and degree of understanding afforded by computational chemistry cannot be understated. The activation strain model is one of the most powerful tools in our arsenal to obtain unparalleled insight into reactivity. The relative energy of interacting reactants is evaluated along a reaction energy profile and related to the rigidity of the reactants' molecular structure and the strength of the stabilizing interactions between the deformed reactants: ΔE(ζ) = ΔEstrain(ζ) + ΔEint(ζ). Owing to the connectedness between the activation strain model and Kohn-Sham molecular orbital theory, one is able to obtain a causal relationship between both the sterics and electronics of the reactants and their mutual reactivity. Only when this is accomplished one can eclipse the phenomenological explanations that are commonplace in the literature and textbooks and begin to rationally tune and optimize chemical transformations. We showcase how the activation strain model is the ideal tool to elucidate fundamental organic reactions, the activation of small molecules by metallylenes, and the cycloaddition reactivity of cyclic diene- and dipolarophiles.

A fundamental and ubiquitous phenomenon in chemistry is the contraction of both C-H and C-C bonds as the carbon atoms involved vary, in s-p hybridization, along sp 3 to sp 2 to sp. Our quantum chemical bonding analyses based on Kohn-Sham molecular orbital theory show that the generally accepted rationale behind this trend is incorrect. Inspection of the molecular orbitals and their corresponding orbital overlaps reveals that the above-mentioned shortening in C-H and C-C bonds is not determined by an increasing amount of s-character at the carbon atom in these bonds. Instead, we establish that this structural trend is caused by a diminishing steric (Pauli) repulsion between substituents around the pertinent carbon atom, as the coordination number decreases along sp 3 to sp 2 to sp., (© 2021 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Organic chemistry has undoubtedly had a profound impact on humanity. Day in and day out, we find ourselves constantly surrounded by organic compounds. Pharmaceuticals, plastics, fuels, cosmetics, detergents, and agrochemicals, to name a few, are all synthesized by organic reactions. Very often, these reactions require a catalyst in order to proceed in a timely and selective manner. Lewis acids and organocatalysts are commonly employed to catalyze organic reactions and are considered to enhance the frontier molecular orbital (FMO) interactions. A vast number of textbooks and primary literature sources suggest that the binding of a Lewis acid or an iminium catalyst to a reactant (R1) stabilizes its LUMO and leads to a smaller HOMO(R2)-LUMO(R1) energy gap with the other reactant (R2), thus resulting in a faster reaction. This forms the basis for the so-called LUMO-lowering catalysis concept. Despite the simplicity and popularity of FMO theory, a number of deficiencies have emerged over the years, as a consequence of these FMOs not being the operative factor in the catalysis. LUMO-lowering catalysis is ultimately incomplete and is not always operative in catalyzed organic reactions. Our groups have recently undertaken a concerted effort to generate a unified framework to rationalize and predict chemical reactivity using a causal model that is rooted in quantum mechanics. In this Account, we propose the concept of Pauli repulsion-lowering catalysis to understand the catalysis in fundamental processes in organic chemistry. Our findings emerge from state-of-the-art computational methods, namely, the activation strain model (ASM) of reactivity in conjunction with quantitative Kohn-Sham molecular orbital theory (KS-MO) and a matching energy decomposition analysis (EDA). The binding of the catalyst to the substrate not only leads to a stabilization of its LUMO but also induces a significant reduction of the two-orbital, four-electron Pauli repulsion involving the key molecular orbitals of both reactants. This repulsion-lowering originates, for the textbook Lewis acid-catalyzed Diels-Alder reaction, from the catalyst polarizing the occupied π orbital of the dienophile away from the carbon atoms that form new bonds with the diene. This polarization of the occupied dienophile π orbital reduces the occupied orbital overlap with the diene and constitutes the ultimate physical factor responsible for the acceleration of the catalyzed process as compared to the analogous uncatalyzed reaction. We show that this physical mechanism is generally applicable regardless of the type of reaction (Diels-Alder and Michael addition reactions) and the way the catalyst is bonded to the reactants (i.e., from pure covalent or dative bonds to weaker hydrogen or halogen bonds). We envisage that the insights emerging from our analysis will guide future experimental developments toward the design of more efficient catalytic transformations.

We have quantum chemically investigated the rotational isomerism of 1,2-dihaloethanes XCH 2 CH 2 X (X = F, Cl, Br, I) at ZORA-BP86-D3(BJ)/QZ4P. Our Kohn-Sham molecular orbital (KS-MO) analyses reveal that hyperconjugative orbital interactions favor the gauche conformation in all cases (X = F-I), not only for X = F as in the current model of this so-called gauche effect. We show that, instead, it is the interplay of hyperconjugation with Pauli repulsion between lone-pair-type orbitals on the halogen substituents that constitutes the causal mechanism for the gauche effect. Thus, only in the case of the relatively small fluorine atoms, steric Pauli repulsion is too weak to overrule the gauche preference of the hyperconjugative orbital interactions. For the larger halogens, X⋅⋅⋅X steric Pauli repulsion becomes sufficiently destabilizing to shift the energetic preference from gauche to anti, despite the opposite preference of hyperconjugation., (© 2021 The Authors. ChemPhysChem published by Wiley-VCH GmbH.)

The so-called peroxidatic cysteines and selenocysteines in proteins reduce hydroperoxides through a dual attack to the peroxide bond in a two-step mechanism. First, a proton dislocation from the thiol/selenol to a close residue of the enzymatic pocket occurs. Then, a nucleophilic attack of the anionic cysteine/selenocysteine to one O atom takes place, while the proton is shuttled back to the second O atom, promoting the formation of a water molecule. In this computational study, we use a molecular model of GPx to demonstrate that the enzymatic environment significantly lowers the barrier of the latter S N 2 step. Particularly, in our Se-based model the energy barriers for the two steps are 29.82 and 2.83 kcal mol -1 , both higher than the corresponding barriers computed in the enzymatic cluster, i. e., 21.60 and null, respectively. Our results, obtained at SMD-B3LYP-D3(BJ)/6-311+G(d,p), cc-pVTZ//B3LYP-D3(BJ)/6-311G(d,p), cc-pVTZ level of theory, show that the mechanistic details can be well reproduced using an oversimplified model, but the energetics is definitively more favorable in the GPx active site. In addition, we pinpoint the role of the chalcogen in the peroxide reduction process, rooting the advantages of the presence of selenium in its acidic and nucleophilic properties., (© 2020 Wiley-VCH GmbH.)

A judiciously oriented external electric field (OEEF) can catalyze a wide range of reactions and can even induce endo/exo stereoselectivity of cycloaddition reactions. The Diels-Alder reaction between cyclopentadiene and maleic anhydride is studied by using quantitative activation strain and Kohn-Sham molecular orbital theory to pinpoint the origin of these catalytic and stereoselective effects. Our quantitative model reveals that an OEEF along the reaction axis induces an enhanced electrostatic and orbital interaction between the reactants, which in turn lowers the reaction barrier. The stronger electrostatic interaction originates from an increased electron density difference between the reactants at the reactive center, and the enhanced orbital interaction arises from the promoted normal electron demand donor-acceptor interaction driven by the OEEF. An OEEF perpendicular to the plane of the reaction axis solely stabilizes the exo pathway of this reaction, whereas the endo pathway remains unaltered and efficiently steers the endo/exo stereoselectivity. The influence of the OEEF on the inverse electron demand Diels-Alder reaction is also investigated; unexpectedly, it inhibits the reaction, as the electric field now suppresses the critical inverse electron demand donor-acceptor interaction., (© 2020 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

The selectivity and rate enhancement of bifunctional hydrogen bond donor-catalyzed Diels-Alder reactions between cyclopentadiene and acrolein were quantum chemically studied using density functional theory in combination with coupled-cluster theory. (Thio)ureas render the studied Diels-Alder cycloaddition reactions exo selective and induce a significant acceleration of this process by lowering the reaction barrier by up to 7 kcal mol -1 . Our activation strain and Kohn-Sham molecular orbital analyses uncover that these organocatalysts enhance the Diels-Alder reactivity by reducing the Pauli repulsion between the closed-shell filled π-orbitals of the diene and dienophile, by polarizing the π-orbitals away from the reactive center and not by making the orbital interactions between the reactants stronger. In addition, we establish that the unprecedented exo selectivity of the hydrogen bond donor-catalyzed Diels-Alder reactions is directly related to the larger degree of asynchronicity along this reaction pathway, which is manifested in a relief of destabilizing activation strain and Pauli repulsion., (© 2020 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

We have studied the activation of dihydrogen by metallylenes using relativistic density functional theory (DFT). Our detailed activation strain and Kohn-Sham molecular orbital analyses have quantified the physical factors behind the decreased reactivity of the metallylene on going down Group 14, from carbenes to stannylenes. Along this series, the reactivity decreases due to a worsening of the back-donation interaction between the filled lone-pair orbital of the metallylene and the σ*-orbital of H 2 , which, therefore, reduces the metallylene-substrate interaction and increases the reaction barrier. As the metallylene ligand is varied from nitrogen to phosphorus to arsenic a significant rate enhancement is observed for the activation of H 2 due to (i) a reduced steric (Pauli) repulsion between the metallylene and the substrate; and (ii) less activation strain, as the metallylene becomes increasingly more predistorted. Using a rationally designed metallylene with an optimal Group 14 atom and ligand combination, we show that a number of small molecules ( i.e. HCN, CO 2 , H 2 , NH 3 ) may also be readily activated. For the first time, we show the ability of our H 2 activated designer metallylenes to hydrogenate unsaturated hydrocarbons. The results presented herein will serve as a guide for the rational design of metallylenes toward the activation of small molecules and subsequent reactions., Competing Interests: The authors declare no conflict of interest., (This journal is © The Royal Society of Chemistry.)

Analyses of multifactorial experimental designs are used as an explorative technique describing hypothesized multifactorial effects based on their variation. The procedure of analyzing multifactorial designs is well established for univariate data, and it is known as analysis of variance (ANOVA) tests, whereas only a few methods have been developed for multivariate data. In this work, we present the weighted-effect ASCA, named WE-ASCA, as an enhanced version of ANOVA-simultaneous component analysis (ASCA) to deal with multivariate data in unbalanced multifactorial designs. The core of our work is to use general linear models (GLMs) in decomposing the response matrix into a design matrix and a parameter matrix, while the main improvement in WE-ASCA is to implement the weighted-effect (WE) coding in the design matrix. This WE-coding introduces a unique solution to solve GLMs and satisfies a constrain in which the sum of all level effects of a categorical variable equal to zero. To assess the WE-ASCA performance, two applications were demonstrated using a biomedical Raman spectral data set consisting of mice colorectal tissue. The results revealed that WE-ASCA is ideally suitable for analyzing unbalanced designs. Furthermore, if WE-ASCA is applied as a preprocessing tool, the classification performance and its reproducibility can significantly improve.

The tuberculosis vaccine bacillus Calmette-Guérin (BCG) protects against some heterologous infections, probably via induction of non-specific innate immune memory in monocytes and natural killer (NK) cells, a process known as trained immunity. Recent studies have revealed that the induction of trained immunity is associated with a bias toward granulopoiesis in bone marrow hematopoietic progenitor cells, but it is unknown whether BCG vaccination also leads to functional reprogramming of mature neutrophils. Here, we show that BCG vaccination of healthy humans induces long-lasting changes in neutrophil phenotype, characterized by increased expression of activation markers and antimicrobial function. The enhanced function of human neutrophils persists for at least 3 months after vaccination and is associated with genome-wide epigenetic modifications in trimethylation at histone 3 lysine 4. Functional reprogramming of neutrophils by the induction of trained immunity might offer novel therapeutic strategies in clinical conditions that could benefit from modulation of neutrophil effector function., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

We have quantum chemically analyzed the competition between the bimolecular nucleophilic substitution (S N 2) and base-induced elimination (E2) pathways for F - + CH 3 CH 2 Cl and PH 2 - + CH 3 Cl using the activation strain model and Kohn-Sham molecular orbital theory at ZORA-OLYP/QZ4P. Herein, we correct an earlier study that intuitively attributed the mechanistic preferences of F 2 and PH - , i.e., E2 and S 2 - , i.e., E2 and S N 2, respectively, to a supposedly unfavorable shift in the polarity of the abstracted β-proton along the PH 2 - -induced E2 pathway while claiming that ″ ...no correlation between the thermodynamic basicity and E2 rate should be expected. ″ Our analyses, however, unequivocally show that it is simply the 6 kcal mol -1 Cl and hence to preferentially react via the E2 pathway, despite the higher characteristic distortivity (more destabilizing activation strain) associated with this pathway. On the other hand, the less basic PH - that enables this base to engage in a more stabilizing orbital interaction with CH 3 CH 2 Cl and is, therefore, unable to overcome the characteristic distortivity of the E2 pathway. Therefore, the mechanistic preference of PH 2 - has a weaker stabilizing interaction with CH 3 CH 2 Cl and is, therefore, unable to overcome the characteristic distortivity of the E2 pathway. Therefore, the mechanistic preference of PH 2 - is steered to the S N 2 reaction channel (less-destabilizing activation strain).

Chains of magnetic atoms with either strong spin-orbit coupling or spiral magnetic order which are proximity-coupled to superconducting substrates can host topologically non-trivial Majorana bound states. The experimental signature of these states consists of spectral weight at the Fermi energy which is spatially localized near the ends of the chain. However, topologically trivial Yu-Shiba-Rusinov in-gap states localized near the ends of the chain can lead to similar spectra. Here, we explore a protocol to disentangle these contributions by artificially augmenting a candidate Majorana spin chain with orbitally-compatible nonmagnetic atoms. Combining scanning tunneling spectroscopy with ab-initio and tight-binding calculations, we realize a sharp spatial transition between the proximity-coupled spiral magnetic order and the non-magnetic superconducting wire termination, with persistent zero-energy spectral weight localized at either end of the magnetic spiral. Our findings open a new path towards the control of the spatial position of in-gap end states, trivial or Majorana, via different chain terminations, and the realization of designer Majorana chain networks for demonstrating topological quantum computation.

We have quantum chemically studied the reactivity, site-, and regioselectivity of the 1,3-dipolar cycloaddition between methyl azide and various allenes, including the archetypal allene propadiene, heteroallenes, and cyclic allenes, by using density functional theory (DFT). The 1,3-dipolar cycloaddition reactivity of linear (hetero)allenes decreases as the number of heteroatoms in the allene increases, and formation of the 1,5-adduct is, in all cases, favored over the 1,4-adduct. Both effects find their origin in the strength of the primary orbital interactions. The cycloaddition reactivity of cyclic allenes was also investigated, and the increased predistortion of allenes, that results upon cyclization, leads to systematically lower activation barriers not due to the expected variations in the strain energy, but instead from the differences in the interaction energy. The geometric predistortion of cyclic allenes enhances the reactivity compared to linear allenes through a unique mechanism that involves a smaller HOMO-LUMO gap, which manifests as more stabilizing orbital interactions., (© 2020 The Authors. Published by Wiley-VCH GmbH.)

A study on the reactivity of the N-heterocyclic silylene Dipp 2 NHSi (1,3-bis(diisopropylphenyl)-1,3-diaza-2-silacyclopent-4-en-2-yliden) with the transition metal complexes [Ni(CO) 4 ], [M(CO) 6 ] (M=Cr, Mo, W), [Mn(CO) 5 (Br)] and [(η 5 -C 5 H 5 )Fe(CO) 2 (I)] is reported. We demonstrate that N-heterocyclic silylenes, the higher homologues of the now ubiquitous NHC ligands, show a remarkably different behavior in coordination chemistry compared to NHC ligands. Calculations on the electronic features of these ligands revealed significant differences in the frontier orbital region which lead to some peculiarities of the coordination chemistry of silylenes, as demonstrated by the synthesis of the dinuclear, NHSi-bridged complex [{Ni(CO) 2 (μ-Dipp 2 NHSi)} 2 ] (2), complexes [M(CO) 5 (Dipp 2 NHSi)] (M=Cr 3, Mo 4, W 5), [Mn(CO) 3 (Dipp 2 NHSi) 2 (Br)] (9) and [(η 5 -C 5 H 5 )Fe(CO) 2 (Dipp 2 NHSi-I)] (10). DFT calculations on several model systems [Ni(L)], [Ni(CO) 3 (L)], and [W(CO) 5 (L)] (L=NHC, NHSi) reveal that carbenes are typically the much better donor ligands with a larger intrinsic strength of the metal-ligand bond. The decrease going from the carbene to the silylene ligand is mainly caused by favorable electrostatic contributions for the NHC ligand to the total bond strength, whereas the orbital interactions were often found to be higher for the silylene complexes. Furthermore, we have demonstrated that the contribution of σ- and π-interaction depends significantly on the system under investigation. The σ-interaction is often much weaker for the NHSi ligand compared to NHC but, interestingly, the π-interaction prevails for many NHSi complexes. For the carbonyl complexes, the NHSi ligand is the better σ-donor ligand, and contributions of π-symmetry play only a minor role for the NHC and NHSi co-ligands., (© 2020 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

In a liquid-STM setup environment, the redox behavior of manganese porphyrins was studied at various solid-liquid interfaces. In the presence of a solution of Mn(III)Cl porphyrins in 1-phenyloctane, which was placed at a conductive surface, large and constant additional currents relative to a set tunneling current were observed, which varied with the magnitude of the applied bias voltage. These currents occurred regardless of the type of surface (HOPG or Au(111)) or tip material (PtIr, Au or W). The additional currents were ascribed to the occurrence of redox reactions in which chloride is oxidized to chlorine and the Mn(III) center of the porphyrin moiety is reduced to Mn(II). The resulting Mn(II) porphyrin products were identified by UV-vis analysis of the liquid phase. For solutions of Mn(III) porphyrins with non-redox active acetate instead of chloride axial ligands, the currents remained absent., (Copyright © 2020, Habets et al.; licensee Beilstein-Institut.)

In this combined experimental and theoretical study, a computational protocol is reported to predict the excited states in D-π-A compounds containing the B( F Xyl) 2 ( F Xyl = 2,6-bis(trifluoromethyl)phenyl) acceptor group for the design of new thermally activated delayed fluorescence (TADF) emitters. To this end, the effect of different donor and π-bridge moieties on the energy gaps between local and charge-transfer singlet and triplet states is examined. To prove this computationally aided design concept, the D-π-B( F Xyl) 2 compounds 1 - 5 were synthesized and fully characterized. The photophysical properties of these compounds in various solvents, polymeric film, and in a frozen matrix were investigated in detail and show excellent agreement with the computationally obtained data. Furthermore, a simple structure-property relationship is presented on the basis of the molecular fragment orbitals of the donor and the π-bridge, which minimize the relevant singlet-triplet gaps to achieve efficient TADF emitters., Competing Interests: The authors declare no conflict of interest., (© 2020 The Authors. Published by WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim.)

αB-crystallin is highly expressed in the heart and slow skeletal muscle; however, the roles of αB-crystallin in the muscle are obscure. Previously, we showed that αB-crystallin localizes at the sarcomere Z-bands, corresponding to the focal adhesions of cultured cells. In myoblast cells, αB-crystallin completely colocalizes with microtubules and maintains cell shape and adhesion. In this study, we show that in beating cardiomyocytes α-tubulin and αB-crystallin colocalize at the I- and Z-bands of the myocardium, where it may function as a molecular chaperone for tubulin/microtubules. Fluorescence recovery after photobleaching (FRAP) analysis revealed that the striated patterns of GFP-αB-crystallin fluorescence recovered quickly at 37°C. FRAP mobility assay also showed αB-crystallin to be associated with nocodazole-treated free tubulin dimers but not with taxol-treated microtubules. The interaction of αB-crystallin and free tubulin was further confirmed by immunoprecipitation and microtubule sedimentation assay in the presence of 1-100 μM calcium, which destabilizes microtubules. Förster resonance energy transfer analysis showed that αB-crystallin and tubulin were at 1-10 nm apart from each other in the presence of colchicine. These results suggested that αB-crystallin may play an essential role in microtubule dynamics by maintaining free tubulin in striated muscles, such as the soleus or cardiac muscles., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Japanese Biochemical Society. All rights reserved.)

Ultraviolet-visible diffuse reflectance spectroscopy (UV-Vis DRS) combined with chemometrics was used for the first time to differentiate Angelicae Sinensis Radix (ASR) from four other similar herbs (either from the same genus or of similar appearance). A total of 191 samples, including 40 ASR, 39 Angelicae Pubescentis Radix (APR), 38 Chuanxiong Rhizoma (CR), 35 Atractylodis Macrocephalae Rhizoma (AMR) and 39 Angelicae Dahuricae Radix (ADR), were collected and divided into the training and prediction sets. Principal component analysis (PCA) was used for observing the sample cluster tendency of the calibration set. Different preprocessing methods were investigated and the optimal preprocessing combination was selected according to spectral signal characteristics and three-dimensional PCA (3D PCA) clustering results. The final discriminant model was built using extreme learning machine (ELM). The exploratory studies on the raw spectra and their 3D PCA scores indicate that the classification of the five herbs cannot be achieved by PCA of the raw spectra. Autoscaling, continuous wavelet transform (CWT) and Savitzky-Golay (SG) smoothing can improve the clustering results to different degrees. Furthermore, their combination in the order of CWT + autoscaling + SG smoothing can enhance the spectral resolution and obtain the best clustering result. These results are also validated using ELM models of raw and different preprocessing methods. By using CWT + autoscaling + SG smoothing + ELM, 100% classification accuracy can be achieved in both the calibration set and the prediction set. Therefore, the developed method could be used as a rapid, economic and effective method for discriminating the five herbs used in this study.

The Diels-Alder reactions between cyclopentadiene and various α,β-unsaturated aldehyde, imine, and iminium dienophiles were quantum chemically studied using a combined density functional theory and coupled-cluster theory approach. Simple iminium catalysts accelerate the Diels-Alder reactions by lowering the reaction barrier up to 20 kcal mol -1 compared to the parent aldehyde and imine reactions. Our detailed activation strain and Kohn-Sham molecular orbital analyses reveal that the iminium catalysts enhance the reactivity by reducing the steric (Pauli) repulsion between the diene and dienophile, which originates from both a more asynchronous reaction mode and a more significant polarization of the π-system away from the incoming diene compared to aldehyde and imine analogs. Notably, we establish that the driving force behind the asynchronicity of the herein studied Diels-Alder reactions is the relief of destabilizing steric (Pauli) repulsion and not the orbital interaction between the terminal carbon of the dienophile and the diene, which is the widely accepted rationale., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)