Reaction Coordinate Research Articles

Restriction on molecular fluxionality by substitution: A case study for the 1,10-dicyanobullvalene.

We show herein that 1,10-dicyano substitution restricts the paragon fluxionality of bullvalene to just 14 isomers which isomerize along a single cycle. The restricted fluxionality of 1,10-dicyanobullvalene (DCB) is investigated by means of: (i) Bonding analyses of the isomer structures using the adaptive natural density partitioning (AdNDP). (ii) Quantum dynamical simulations of the isomerizations along the cyclic intrinsic reaction coordinate of the potential energy surface (PES). The PES possesses 14 equivalent potential wells supporting 14 isomers which are separated by 14 equivalent potential barriers supporting 14 transition states. Accordingly, at low temperatures, DCB appears as a hindered molecular rotor, without any delocalization of the wavefunction in the 14 potential wells, without any nuclear spin isomers, and with completely negligible tunneling. These results are compared and found to differ from those for molecular boron rotors. (iii) Born-Oppenheimer molecular dynamics (BOMD) simulations of thermally activated isomerizations. (iv) Calculations of the rate constants in the frame of transition state theory (TST) with reasonable agreement achieved with the BOMD results. (v) Simulations of the equilibration dynamics using rate equations for the isomerizations with TST rate coefficients. Accordingly, in the long-time limit, isomerizations of the 14 isomers, each with Cs symmetry, approach the "14 Cs → C7v" thermally averaged structure. This is a superposition of the 14 equally populated isomer structures with an overall C7v symmetry. By extrapolation, the results for DCB yield working hypotheses for so far un-explored properties e.g. for the equilibration dynamics of C10H10.

ℏ 4 quantum corrections to semiclassical transmission probabilities

The combination of vibrational perturbation theory with the replacement of the harmonic oscillator quantization condition along the reaction coordinate with an imaginary action to be used in the uniform semiclassical approximation for the transmission probability has been shown in recent years to be a practical method for obtaining thermal reaction rates. To date, this theory has been developed systematically only up to second order in perturbation theory. Although it gives the correct leading order term in an ℏ2 expansion, its accuracy at lower temperatures, where tunneling becomes important, is not clear. In this paper, we develop the theory to fourth order in the action. This demands developing the quantum perturbation theory up to sixth order. Remarkably, we find that the fourth order theory gives the correct ℏ4 term in the expansion of the exact thermal rate. The relative magnitude of the fourth order correction as compared to the second order term objectively indicates the accuracy of the second order theory. We also extend the previous modified second order theory to the fourth order case, creating an ℏ2 modified potential for this purpose. The resulting theory is tested on the standard examples—symmetric and asymmetric Eckart potentials and a Gaussian potential. The modified fourth order theory is remarkably accurate for the asymmetric Eckart potential.

QM/MM Free Energy Calculations of Long-Range Biological Protonation Dynamics by Adaptive and Focused Sampling.

Water-mediated proton transfer reactions are central for catalytic processes in a wide range of biochemical systems, ranging from biological energy conversion to chemical transformations in the metabolism. Yet, the accurate computational treatment of such complex biochemical reactions is highly challenging and requires the application of multiscale methods, in particular hybrid quantum/classical (QM/MM) approaches combined with free energy simulations. Here, we combine the unique exploration power of new advanced sampling methods with density functional theory (DFT)-based QM/MM free energy methods for multiscale simulations of long-range protonation dynamics in biological systems. In this regard, we show that combining multiple walkers/well-tempered metadynamics with an extended system adaptive biasing force method (MWE) provides a powerful approach for exploration of water-mediated proton transfer reactions in complex biochemical systems. We compare and combine the MWE method also with QM/MM umbrella sampling and explore the sampling of the free energy landscape with both geometric (linear combination of proton transfer distances) and physical (center of excess charge) reaction coordinates and show how these affect the convergence of the potential of mean force (PMF) and the activation free energy. We find that the QM/MM-MWE method can efficiently explore both direct and water-mediated proton transfer pathways together with forward and reverse hole transfer mechanisms in the highly complex proton channel of respiratory Complex I, while the QM/MM-US approach shows a systematic convergence of selected long-range proton transfer pathways. In this regard, we show that the PMF along multiple proton transfer pathways is recovered by combining the strengths of both approaches in a QM/MM-MWE/focused US (FUS) scheme and reveals new mechanistic insight into the proton transfer principles of Complex I. Our findings provide a promising basis for the quantitative multiscale simulations of long-range proton transfer reactions in biological systems.

Physical Prior Mean Function-Driven Gaussian Processes Search for Minimum-Energy Reaction Paths with a Climbing-Image Nudged Elastic Band: A General Method for Gas-Phase, Interfacial, and Bulk-Phase Reactions.

The climbing-image nudged elastic band (CI-NEB) method serves as an indispensable tool for computational chemists, offering insight into minimum-energy reaction paths (MEPs) by delineating both transition states (TSs) and intermediate nonstationary structures along reaction coordinates. However, executing CI-NEB calculations for reactions with extensive reaction coordinate spans necessitates a large number of images to ensure a reliable convergence of the MEPs and TS structures, presenting a computationally demanding optimization challenge, even with mildly costly electronic-structure methods. In this study, we advocate for the utilization of physically inspired prior mean function-based Gaussian processes (GPs) to expedite MEP exploration and TS optimization via the CI-NEB method. By incorporating reliable prior physical approximations into potential energy surface (PES) modeling, we demonstrate enhanced efficiency in multidimensional CI-NEB optimization with surrogate-based optimizers. Our physically informed GP approach not only outperforms traditional nonsurrogate-based optimizers in optimization efficiency but also on-the-fly learns the reaction path valley during optimization, culminating in significant advancements. The surrogate PES derived from our optimization exhibits high accuracy compared to true PES references, aligning with our emphasis on leveraging reliable physical priors for robust and efficient posterior mean learning in GPs. Through a systematic benchmark study encompassing various reaction pathways, including gas-phase, bulk-phase, and interfacial/surface reactions, our physical GPs consistently demonstrate superior efficiency and reliability. For instance, they outperform the popular fast inertial relaxation engine optimizer by approximately a factor of 10, showcasing their versatility and efficacy in exploring reaction mechanisms and surface reaction PESs.

Understanding the Energy Landscape of Intrinsically Disordered Protein Ensembles.

A substantial portion of various organisms' proteomes comprises intrinsically disordered proteins (IDPs) that lack a defined three-dimensional structure. These IDPs exhibit a diverse array of conformations, displaying remarkable spatiotemporal heterogeneity and exceptional conformational flexibility. Characterizing the structure or structural ensemble of IDPs presents significant conceptual and methodological challenges owing to the absence of a well-defined native structure. While databases such as the Protein Ensemble Database (PED) provide IDP ensembles obtained through a combination of experimental data and molecular modeling, the absence of reaction coordinates poses challenges in comprehensively understanding pertinent aspects of the system. In this study, we leverage the energy landscape visualization method (JCTC, 6482, 2019) to scrutinize four IDP ensembles sourced from PED. ELViM, a methodology that circumvents the need for a priori reaction coordinates, aids in analyzing the ensembles. The specific IDP ensembles investigated are as follows: two fragments of nucleoporin (NUL: 884-993 and NUS: 1313-1390), yeast sic 1 N-terminal (1-90), and the N-terminal SH3 domain of Drk (1-59). Utilizing ELViM enables the comprehensive validation of ensembles, facilitating the detection of potential inconsistencies in the sampling process. Additionally, it allows for identifying and characterizing the most prevalent conformations within an ensemble. Moreover, ELViM facilitates the comparative analysis of ensembles obtained under diverse conditions, thereby providing a powerful tool for investigating the functional mechanisms of IDPs.

Concerning for the solvent-polarity-dependent conformational equilibrium and ESIPT mechanism in Pz3HC system: A novel insight

In this report, we propose a new insight into the interaction between the solvent-polarity-dependent conformational equilibrium and excited state intramolecular proton transfer (ESIPT) behavior of Pz3HC system in four different polar solvents (polarity order: ACN > THF > TOL > CYC). Using quantum chemistry method, we first announce a coexistence mechanism between Pz3HC-1 and Pz3HC-3 in the ground state in four solvents based on the Boltzmann distribution. In particular, Pz3HC-1 is the principal configuration in non-polar solvent, but Pz3HC-3 is the principal configuration in polar solvent. In addition, the simulated fluorescence spectra interprets the negative solvatochromism effect of Pz3HC-1 and Pz3HC-3 in four solvents. The evidence from intramolecular hydrogen bonding (IHB) parameters and electronic perspective collectively confirms the light-induced IHB enhancement and intramolecular charge transfer (ICT) properties in Pz3HC-1 and Pz3HC-3, which raises the likelihood of the ESIPT process. Combining the calculation of potential energy curve (PEC) and intrinsic reaction coordinate (IRC), we demonstrate that the ESIPT ease of Pz3HC-1 in different polar solvents obeys the order of CYC > TOL > THF > ACN, while the order of ESIPT ease in Pz3HC-3 is opposite. Notably, the ESIPT process of Pz3HC-3 in CYC solvent is accompanied by the twisted intramolecular charge transfer (TICT) process. In addition, we also reveal that the enol* and keto* fluorescence peaks of Pz3HC-3 in CYC solvent are quenched by ISC and TICT process, respectively. Our work not only provides a satisfactory explanation of the novel dynamics mechanism for Pz3HC system, but also brings light to the design and application of new sensing molecules in the future.

Deep learning path-like collective variable for enhanced sampling molecular dynamics.

Several enhanced sampling techniques rely on the definition of collective variables to effectively explore free energy landscapes. The existing variables that describe the progression along a reactive pathway offer an elegant solution but face a number of limitations. In this paper, we address these challenges by introducing a new path-like collective variable called the "deep-locally non-linear-embedding," which is inspired by principles of the locally linear embedding technique and is trained on a reactive trajectory. The variable mimics the ideal reaction coordinate by automatically generating a non-linear combination of features through a differentiable generalized autoencoder that combines a neural network with a continuous k-nearest neighbor selection. Among the key advantages of this method is its capability to automatically choose the metric for searching neighbors and to learn the path from state A to state B without the need to handpick landmarks a priori. We demonstrate the effectiveness of DeepLNE by showing that the progression along the path variable closely approximates the ideal reaction coordinate in toy models, such as the Müller-Brown potential and alanine dipeptide. Then, we use it in the molecular dynamics simulations of an RNA tetraloop, where we highlight its capability to accelerate transitions and estimate the free energy of folding.

Theoretical investigation on the promotional role of Ag for reverse water-gas shift reaction on Ni catalysts during CO2 reduction

Catalytic reduction of CO2 to CO via reverse water-gas shift reaction on Ni-based catalysts is challenging due to the competing methanation reaction. In this work, the reduction of CO2 toward CO and CH4 has been studied by density functional theory calculation and microkinetic simulations on Ni(211) and Ag@Ni(211) surfaces. On both surfaces, direct CO2 dissociation is more favorable over formate and carboxyl pathways toward CO, and direct CO dissociation is the major pathway toward CH4. The presence of Ag shifts the d-band center of adjacent Ni atoms away from the Fermi level, resulting in reduced affinity to CO2/CO/intermediates and slightly increased barriers for most elementary steps. In particular, the apparent barrier for direct CO2 dissociation is slightly reduced by ∼0.1 eV while direct CO dissociation is strongly inhibited by ∼0.3 eV along the reaction coordinate, leading to CO as the predominant product on Ag@Ni(211), which agrees well with the experimental results.

Preparation and performance of carboxyl modified NdFeB and derived magnetic self-healing PU composites

With the rapid development of today’s society, magnetic self-healing materials are widely used in many fields due to their excellent magnetism and self-healing ability. This article presents a magnetic self-healing composite with good mechanical, self-healing and magnetic properties and magnetic drive performance. Tetraethyl orthosilicate (TEOS), silane coupling agent and succinic anhydride (SA) were used for surface modification of NdFeB magnetic particles (MPs) and NdFeB@SiO2-COOH MPs was synthesized. Iron(III) ions were added to the self-polymerization process of the natural small molecule thioctic acid (TA) to form the PTA-Fe prepolymer. PTA-Fe prepolymer, isophorone diisocyanate (IPDI) and polyether amine (PEA) reacted to synthesis PTA-Fe-PU matrix with good self-healing performance at room temperature. The prepared NdFeB@SiO2-COOH MPs was mixed with the self-healing matrix and cured to obtain the magnetic room-temperature self-healing composite MPTA-Fe-PU. The best tensile strength of the composite was 7.31 MPa and the ultimate elongation was 833 %. After the prepared MPTA-Fe-PU sample strip was completely cut, the incision on the surface almost disapeared and the tensile strength recovered to 6.35 MPa after 12 h self-healing at room temperature. The self-healing efficiency was up to 87 %. The highly efficient self-healing properties were mainly derived from the disulfide bonds exchange reaction and iron(III) ion coordination in the matrix.

Investigating the catalytic activity of Mgn (n = 4–8) clusters for the hydrogen evolution reaction using density functional theory

AbstractTo efficiently desorb H2, pure Mgn (n = 4–8) clusters were chosen for the hydrogen evolution reaction with H2O. At the PBE0/def2‐TZVP level and the PBE0‐D3/def2‐TZVP level, the lowest energy structures of Mgn (n = 4–8) clusters and the most stable structures of Mgn@H2O (n = 4–8) complexes were searched in the local region. The transition state was predicted, and then the hydrogen evolution reaction channel was obtained by using the intrinsic reaction coordinate (IRC) to confirm the transition state. To better analyze the hydrogen reaction mechanism, the character of Mgn@H2O (n = 4–8) complexes and MgnO (n = 4–8) clusters, as well as the atomic charge change trend, were investigated using interaction region indicator function analysis (IRI) and natural population analysis (NPA). The reaction effect of Mg4 cluster and H2O is the worst. The energy barrier does, however, progressively lower as the cluster atom count rises, improving the reaction effect.

Effects of Ionic Liquids on the Nucleofugality of Bromide.

The nucleofugality of bromide was measured in solvent mixtures containing ionic liquids. The solvolysis rate constants of the bromides of well-defined electrofuges were determined in mixtures containing different proportions of 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide in ethanol. Temperature-dependent kinetic studies allowed an explanation of the observed solvent effects in different mixtures in terms of interactions in solution. Using the solvolysis data, the nucleofugality of bromide in these systems was determined. Likewise, nucleofugality data for bromide were determined in mixtures containing high proportions of seven further ionic liquids. These data allowed quantification of the effects of both varying the amount of ionic liquid and the nature of ionic liquid components on the nucleofugality of bromide. Importantly, ionic liquid mixtures were shown to affect the nucleofugality in a manner similar to chloride, providing a method for predicting the effects of ionic liquids on other electrofuges. Further, the ionic liquids were shown to move the transition state earlier along the reaction coordinate, meaning that there is less charge development in the transition state.

C-H Bond Oxidation by MnIV-Oxo Complexes: Hydrogen-Atom Tunneling and Multistate Reactivity.

The reactivity of six MnIV-oxo complexes in C-H bond oxidation has been examined using a combination of kinetic experiments and computational methods. Variable-temperature studies of the oxidation of 9,10-dihydroanthracene (DHA) and ethylbenzene by these MnIV-oxo complexes yielded activation parameters suitable for evaluating electronic structure computations. Complementary kinetic experiments of the oxidation of deuterated DHA provided evidence for hydrogen-atom tunneling in C-H bond oxidation for all MnIV-oxo complexes. These results are in accordance with the Bell model, where tunneling occurs near the top of the transition-state barrier. Density functional theory (DFT) and DLPNO-CCSD(T1) computations were performed for three of the six MnIV-oxo complexes to probe a previously predicted multistate reactivity model. The DFT computations predicted a thermal crossing from the 4B1 ground state to a 4E state along the C-H bond oxidation reaction coordinate. DLPNO-CCSD(T1) calculations further confirm that the 4E transition state offers a lower energy barrier, reinforcing the multistate reactivity model for these complexes. We discuss how this multistate model can be reconciled with recent computations that revealed that the kinetics of C-H bond oxidation by this set of MnIV-oxo complexes can be well-predicted on the basis of the thermodynamic driving force for these reactions.

A review on dynamics of permeability-glycoprotein in efflux of chemotherapeutic drugs

Permeability-glycoprotein (P-gp) belongs to the ABS transporter protein family, with a high expression rate in cancerous cells. The substrate/inhibitors of the protein are structurally diverse, with no lucid mechanism of inhibition. There are two schools of thought on the inhibition mechanism: (i) P-gp inhibitors bind to the huge hydrophobic cavity between two Trans-Membrane Domains (TMDs), supported by ample literary proof and (ii) P-gp inhibitors bind to the vicinity of Nucleotide-Binding Sites (NBSs). Structural biologists have presented several experimental and theoretical structures of P-gp with bound nucleotides and inhibitors to explain the same. However, the available experimental P-gp structures are insufficient to address the catalytic transition path of mammalian P-gp in detail, thus the dynamics and mechanism by which drugs are effluxed is still unknown. Targeted Molecular Dynamics (targeted MD) could be used to minutely analyse and explore the catalytic transition inward open (IO) to outward open (OO) and relaxation path (OO to IO). Finally, analysis of targeted MD trajectory may help to explore different conformational states of Pg-p (reaction coordinate of catalytic transition/relaxation), efflux of compounds aided by the dynamics of Nucleotide Binding Domains/NBDs (ATP coupled process) and TMDs (peristalsis-like movement pushes the bound molecule). This review presents an understanding of the catalytic transition and dynamics of protein which provides insights at the efflux of chemotherapeutic drug using in cancer treatment.

Computational efficiency study of a micro-macro Markov chain Monte Carlo method for molecular dynamics

We study the numerical properties of a recently introduced micro-macro Markov chain Monte Carlo (mM-MCMC) scheme that accelerates sampling of Gibbs distributions when there is a time-scale separation between the complete (microscopic) molecular dynamics and the slow dynamics of a low dimensional reaction coordinate. The micro-macro Markov chain works in three steps: 1) compute the reaction coordinate value associated to the current molecular state; 2) generate a new macroscopic proposal using a proposal distribution; 3) reconstruct a molecular configuration that is consistent with the newly sampled macroscopic value. Here, we systematically study the impact of method parameters on the efficiency of the micro-macro Markov chain Monte Carlo method, namely the chosen macroscopic dynamics and the microscopic reconstruction distribution. We specifically investigate the impact of the macroscopic and reconstruction proposal distributions and their parameters on this efficiency on two prototypical molecular systems, a three-atom molecule and butane. We find, through detailed computational experiments, that the efficiency is largest when these proposal distributions are close to their time-invariant counterparts and that it is robust with respect to parameter changes.

Molecular Simulations Reveal the Free Energy Landscape and Transition State of Membrane Electroporation.

The formation of pores over lipid membranes by the application of electric fields, termed membrane electroporation, is widely used in biotechnology and medicine to deliver drugs, vaccines, or genes into living cells. Continuum models for describing the free energy landscape of membrane electroporation were proposed decades ago, but they have never been tested against spatially detailed atomistic models. Using molecular dynamics (MD) simulations with a recently proposed reaction coordinate, we computed potentials of mean force of pore nucleation and pore expansion in lipid membranes at various transmembrane potentials. Whereas the free energies of pore expansion are compatible with previous continuum models, the experimentally important free energy barrier of pore nucleation is at variance with established models. The discrepancy originates from different geometries of the transition state; previous continuum models assumed the presence of a membrane-spanning defect throughout the process, whereas, according to the MD simulations, the transition state of pore nucleation is typically passed before a transmembrane defect has formed. A modified continuum model is presented that qualitatively agrees with the MD simulations. Using kinetics of pore opening together with transition state theory, our free energies of pore nucleation are in excellent agreement with previous experimental data.

Tandem hybridization chain reaction and selective coordination enable fluorescence detection of exosomes in lung cancer

Liquid biopsy based on tumor-derived exosomes (TDEs) is potentially effective in assisting early cancer diagnosis. However, accurate quantification of the approach remains challenging due to the low concentration and heterogeneity of TDEs. This work involved designing a homogeneous enzyme-free fluorescent sensing strategy by targeting programmed cell death ligand 1 (PD-L1) that was overexpressed on the surface of TDEs, coupled with the tandem hybridization chain reaction (HCR) and the selective coordination of Cu2+. The aptamer was used as a specific recognition probe for binding with the target protein, thereby regulating the tandem HCR reactions and achieving signal amplification. The coordination between Cu2+ and double-stranded DNA (dsDNA) can further regulate the fluorescence intensity of the signal molecule, calcein, for precise quantification of TDEs. The results showed that the detection limit (LOD) for TDEs can be as low as 100 particles/mL. Moreover, the proposed strategy can accurately differentiate between lung cancer patients (n = 25) and healthy donors (n = 12), with a sensitivity of 92%, specificity of 100%, and an area under the receiver operating characteristic curve (AUC) of 0.963. The fluorescence analysis also demonstrated a strong correlation with the computed tomography scans and pathological analysis. Therefore, the proposed method in this study can provide more insights into assisting in the early diagnosis of lung cancer through TDEs-based liquid biopsy.

Intermolecular interaction of azide, cyano and alkyne-N-phenethylacetamide dimers: Experimental and quantum chemical approach

In this study, we present the synthesis and structure characterization by single-crystal X-ray Diffraction (SC-XRD), Ultraviolet Spectroscopy (UV), Fourier Transform Infrared Spectroscopy (FTIR), Nuclear Magnetic Resonance (NMR) and High-resolution Mass Spectrometry (HRMS) of 2-azido-N-phenethylacetamide (4) and 2-cyano-N-phenethylacetamide (6). Further, Density Functional Theory (DFT) was employed to get a better insight regarding the properties exhibited by the compounds. The computational chemistry has been used to explain and investigate the bimolecular nucleophilic substitution (SN2) and concerted acyl substitution mechanisms for the synthesis of compounds 4 and 6, respectively. Additionally, intrinsic reaction coordinate (IRC) calculations confirm the transition state for these reactions. The X-ray diffraction analysis showed that the compounds 4 and 6 crystallized as monoclinic systems, in space group P21/c and P21/n, respectively. The Hirshfeld Surface Analysis revealed that the main contributions in the crystal packing have H···H, H···N/N···H, H···C/C···H and H···O/O···H interactions. The frontier molecular orbital energies were calculated and the HOMO-LUMO energy gap of 4 and 6 were 10.011 and 10.278 eV, respectively. Theoretical values of electronic transitions, vibrational frequencies and NMR chemical shifts were calculated and compared with experimental results, they all showed good results.

Derivation of the Invariant Free-Energy Landscape Based on Langevin Dynamics.

A conventionally defined free-energy landscape (FEL) exhibits unphysical dependence on the choice of reaction coordinates and hence lacks universal predictive ability. We here show that three physically plausible requirements uniquely determine the FEL formula for a given reaction coordinate. Our FEL is expressed solely in terms of quantities obtained through time-series data analysis, namely, the probability distribution and the diffusion matrix. It is free from any unphysical coordinate dependence and coincides with the conventional FEL in special cases. The uniqueness and robustness of the formula strongly suggest that our FEL has universal predictive power.

The reaction mechanism of the Ideonella sakaiensis PETase enzyme

Polyethylene terephthalate (PET), the most abundantly produced polyester plastic, can be depolymerized by the Ideonella sakaiensis PETase enzyme. Based on multiple PETase crystal structures, the reaction has been proposed to proceed via a two-step serine hydrolase mechanism mediated by a serine-histidine-aspartate catalytic triad. To elucidate the multi-step PETase catalytic mechanism, we use transition path sampling and likelihood maximization to identify optimal reaction coordinates for the PETase enzyme. We predict that deacylation is likely rate-limiting, and the reaction coordinates for both steps include elements describing nucleophilic attack, ester bond cleavage, and the “moving-histidine” mechanism. We find that the flexibility of Trp185 promotes the reaction, providing an explanation for decreased activity observed in mutations that restrict Trp185 motion. Overall, this study uses unbiased computational approaches to reveal the detailed reaction mechanism necessary for further engineering of an important class of enzymes for plastics bioconversion.

Role of the Vibrational and Translational Energies in the CN(v)+C2H6(ν1, ν2, ν5 and ν9) Reactions. A Theoretical QCT Study.

Quasi-classical trajectory (QCT) calculations were conducted on the newly developed full-dimensional potential energy surface, PES-2023, to analyse two critical aspects: the influence of vibrational versus translational energy in promoting reactivity, and the impact of vibrational excitation within similar vibrational modes. The former relates to Polanyi's rules, while the latter concerns mode selectivity. Initially, the investigation revealed that independent vibrational excitation by a single quantum of ethane's symmetric and asymmetric stretching modes (differing by only 15 cm-1) yielded comparable dynamics, reaction cross-sections, HCN(v) vibrational product distributions, and scattering distributions. This observation dismisses any significant mode selectivity. Moreover, an equivalent amount of energy provided as translational energy (at total energies of 9.6 and 20.0 kcal mol-1) gave rise to slightly lower reactivity compared to the same amount of energy provided as vibrational energy. This effect is more evident at low energies, presenting a counterintuitive scenario in an 'early transition state' reaction. These findings challenge the straightforward application of Polanyi's rules in polyatomic systems. Regarding CN(v) vibrational excitation, our calculations reveal that the reaction cross-section remains practically unaffected by this vibrational excitation, suggesting that the CN stretching mode is a spectator mode. The results were rationalized by considering several factors: the strong coupling between different vibrational modes, and between vibrational modes and the reaction coordinate; and a significant vibrational energy redistribution within the ethane reactant before collision. This redistribution creates an unphysical energy flow, resulting in loss of adiabaticity and vibrational memory before the reactants' collision. These theoretical findings require future confirmation through experimental or theoretical quantum mechanical studies, which are currently unavailable.