Exogenous Glutathione Protects IPEC-J2 Cells against Oxidative Stress through a Mitochondrial Mechanism

The accumulation of reactive oxygen species (ROS) triggers oxidative stress in cells by oxidizing and modifying various cellular components, preventing them from performing their inherent functions, ultimately leading to apoptosis and autophagy. Glutathione (GSH) is a ubiquitous intracellular peptide with multiple functions. In this study, a hydrogen peroxide (H2O2)-induced oxidative damage model in IPEC-J2 cells was used to investigate the cellular protection mechanism of exogenous GSH against oxidative stress. The results showed that GSH supplement improved the cell viability reduced by H2O2-induced oxidative damage model in IPEC-J2 cells in a dose-dependent manner. Moreover, supplement with GSH also attenuated the H2O2-induced MMP loss, and effectively decreased the H2O2-induced mitochondrial dysfunction by increasing the content of mtDNA and upregulating the expression TFAM. Exogenous GSH treatment significantly decreased the ROS and MDA levels, improved SOD activity in H2O2-treated cells and reduced H2O2-induced early apoptosis in IPEC-J2 cells. This study showed that exogenous GSH can protect IPEC-J2 cells against apoptosis induced by oxidative stress through mitochondrial mechanisms.     

Anharmonic Vibrational Frequencies of Water Borane and Associated Molecules

Water borane (BH3OH2) and borinic acid (BH2OH) have been proposed as intermediates along the pathway of hydrogen generation from simple reactants: water and borane. However, the vibrational spectra for neither water borane nor borinic acid has been investigaged experimentally due to the difficulty of isolating them in the gas phase, making accurate quantum chemical predictions for such properties the most viable means of their determination. This work presents theoretical predictions of the full rotational and fundamental vibrational spectra of these two potentially application-rich molecules using quartic force fields at the CCSD(T)-F12b/cc-pCVTZ-F12 level with additional corrections included for the effects of scalar relativity. This computational scheme is further benchmarked against the available gas-phase experimental data for the related borane and HBO molecules. The differences are found to be within 3 cm−1 for the fundamental vibrational frequencies and as close as 15 MHz in the B0 and C0 principal rotational constants. Both BH2OH and BH3OH2 have multiple vibrational modes with intensities greater than 100 km mol−1, namely ν2 and ν4 in BH2OH, and ν1, ν3, ν4, ν9, and ν13 in BH3OH2. Finally, BH3OH2 has a large dipole moment of 4.24 D, which should enable it to be observable by rotational spectroscopy, as well.     

Photochemistry of Thymine in Protic Polar Nanomeric Droplets Using Electrostatic Embeding TD-DFT/MM

Thymine photochemistry is important for understanding DNA photodamage. In the gas phase, thymine undergoes a fast non-radiative decay from S 2 to S 1 . In the S 1 state, it gets trapped for several picoseconds until returning to the ground-state S 0 . Here, we explore the electrostatic effects of nanomeric droplets of methanol and water on the excited states of thymine. For this purpose, we develop and implement an electrostatic embedding TD-DFT/MM method based on a QM/MM coupling defined through electrostatic potential fitting charges. We show that both in methanol and water, the mechanism is similar to the gas phase. The solvent molecules participate in defining the branching plane of S 0 /S 1 intersection and have a negligible effect on the S 1 /S 2 intersection. Despite the wrong topology of the ground/excited state intersections, electrostatic embedding TD-DFT/MM allows for a fast exploration of the potential energy surfaces and a qualitative picture of the photophysics of thymine in solvent droplets.     

Numerical Interchain Mean-Field Theory for the Specific Heat of the Bimetallic Ferromagnetically Coupled Chain Compound MnNi(NO2)4(en)2 (en = Ethylenediamine)

We present a detailed study of the field-dependent specific heat of the bimetallic ferromagnetically coupled chain compound MnNi(NO2)4(en)2, en = ethylenediamine. For this material, which in zero field orders antiferromagnetically below TN=2.45 K, small fields suppress magnetic order. Instead, in such fields, a double-peak-like structure in the temperature dependence of the specific heat is observed. We attribute this behavior to the existence of an acoustic and an optical mode in the spin-wave dispersion as a result of the existence of two different spins per unit cell. We compare our experimental data to numerical results for the specific heat obtained by exact diagonalization and Quantum Monte Carlo simulations for the alternating spin-chain model, using parameters that have been derived from the high-temperature behavior of the magnetic susceptibility. The interchain coupling is included in the numerical treatment at the mean-field level. We observe remarkable agreement between experiment and theory, including the ordering transition, using previously determined parameters. Furthermore, the observed strong effect of an applied magnetic field on the ordered state of MnNi(NO2)4(en)2 promises interesting magnetocaloric properties.     

On the Coexistence of the Carbene⋯H-D Hydrogen Bond and Other Accompanying Interactions in Forty Dimers of N-Heterocyclic-Carbenes (I,IMe2, IiPr2, ItBu2, IMes2, IDipp2, IAd2; I = imidazol-2-ylidene) and Some Fundamental Proton Donors (HF,HCN, H2O,MeOH, NH3)

The subject of research is forty dimers formed by imidazol-2-ylidene (I) or its derivative (IR2) obtained by replacing the hydrogen atoms in both N-H bonds with larger important and popular substituents of increasing complexity (methyl = Me, iso-propyl = iPr, tert-butyl = tBu, phenyl = Ph, mesityl = Mes, 2,6-diisopropylphenyl = Dipp, 1-adamantyl = Ad) and fundamental proton donor (HD) molecules (HF, HCN, H2O, MeOH, NH3). While the main goal is to characterize the generally dominant C⋯H-D hydrogen bond engaging a carbene carbon atom, an equally important issue is the often omitted analysis of the role of accompanying secondary interactions. Despite the often completely different binding possibilities of the considered carbenes, and especially HD molecules, several general trends are found. Namely, for a given carbene, the dissociation energy values of the IR2⋯HD dimers increase in the following order: NH3< H2O < HCN ≤ MeOH ≪ HF. Importantly, it is found that, for a given HD molecule, IDipp2 forms the strongest dimers. This is attributed to the multiplicity of various interactions accompanying the dominant C⋯H-D hydrogen bond. It is shown that substitution of hydrogen atoms in both N-H bonds of the imidazol-2-ylidene molecule by the investigated groups leads to stronger dimers with HF, HCN, H2O or MeOH. The presented results should contribute to increasing the knowledge about the carbene chemistry and the role of intermolecular interactions, including secondary ones.     

Spectral Signatures of Hydrogen Thioperoxide (HOSH) and Hydrogen Persulfide (HSSH): Possible Molecular Sulfur Sinks in the Dense ISM

For decades, sulfur has remained underdetected in molecular form within the dense interstellar medium (ISM), and somewhere a molecular sulfur sink exists where it may be hiding. With the discovery of hydrogen peroxide (HOOH) in the ISM in 2011, a natural starting point may be found in sulfur-bearing analogs that are chemically similar to HOOH: hydrogen thioperoxide (HOSH) and hydrogen persulfide (HSSH). The present theoretical study couples the accuracy in the anharmonic fundamental vibrational frequencies from the explicitly correlated coupled cluster theory with the accurate rotational constants provided by canonical high-level coupled cluster theory to produce rovibrational spectra for use in the potential observation of HOSH and HSSH. The ν6 mode for HSSH at 886.1 cm−1 is within 0.2 cm−1 of the gas-phase experiment, and the B0 rotational constant for HSSH of 6979.5 MHz is within 9.0 MHz of the experimental benchmarks, implying that the unknown spectral features (such as the first overtones and combination bands) provided herein are similarly accurate. Notably, a previous experimentally-attributed 2ν1 mode, at 7041.8 cm−1, has been reassigned to the ν1+ν5 combination band based on the present work’s ν1+ν5 value at 7034.3 cm−1. The most intense vibrational transitions for each molecule are the torsions, with HOSH having a more intense transition of 72 km/mol compared to HSSH’s intensity of 14 km/mol. Furthermore, HOSH has a larger net dipole moment of 1.60 D compared to HSSH’s 1.15 D. While HOSH may be the more likely candidate of the two for possible astronomical observation via vibrational spectroscopy due to the notable difference in their intensities, both HSSH and HOSH have large enough net dipole moments to be detectable by rotational spectroscopy to discover the role these molecules may have as possible molecular sulfur sinks in the dense ISM.     

Thermodynamic Properties of Crystalline Cellulose Allomorphs Studied with Dispersion-Corrected Density Functional Methods

The phonon properties and thermodynamics of four crystalline cellulose allomorphs, Iα, Iβ, II, and III1, have been investigated using dispersion-corrected density functional theory (DFT). In line with experimental findings, the free energy differences between the studied cellulose allomorphs are small, less than 1 kJ/mol per atom. The calculated specific heat at constant volume (Cv) has been compared with the available experimental data in the temperature range 10–300 K. Quasiharmonic approximation has been employed to study thermodynamics and specific heat at constant pressure (Cp). For the studied temperature range of 10–400 K, the specific heat of all cellulose allomorphs shows very similar behavior. The calculated and experimental specific heat agree well at low temperatures below 100 K, but the deviation between theory and experiment increases with temperature. This may be due to increasing phonon anharmonicity as the temperature increases.     

How Water Interacts with the NOH Group: The Rotational Spectrum of the 1:1 N,N-diethylhydroxylamine·Water Complex

The rotational spectrum of the 1:1 N,N-diethylhydroxylamine-water complex has been investigated using pulsed jet Fourier transform microwave spectroscopy in the 6.5–18.5 GHz frequency region. The most stable conformer has been detected as well as the 13C monosubstituted isotopologues in natural abundance and the 18O enriched water species, allowing to determine the nitrogen nuclear quadrupole coupling constants and the molecular structure in the vibrational ground state. The molecule has a Cs symmetry and the water lies in the bc symmetry plane forming two hydrogen bonds with the NOH frame with length: dHOH·NOH = 1.974 Å and dH2O·HON = 2.096 Å. From symmetry-adapted perturbation theory calculations coupled to atoms in molecule approach, the corresponding interaction energy values are estimated to be 24 and 13 kJ·mol−1, respectively. The great strength of the intermolecular interaction involving the nitrogen atom is in agreement with the high reactivity of hydroxylamine compounds at the nitrogen site.     

Optical Dynamic Nuclear Polarization of 13C Spins in Diamond at a Low Field with Multi-Tone Microwave Irradiation

Majority of dynamic nuclear polarization (DNP) experiments have been requiring helium cryogenics and strong magnetic fields for a high degree of nuclear polarization. In this work, we instead demonstrate an optical hyperpolarization of naturally abundant 13C nuclei in a diamond crystal at a low magnetic field and the room temperature. It exploits continuous laser irradiation for polarizing electronic spins of nitrogen vacancy centers and microwave irradiation for transferring the electronic polarization to 13C nuclear spins. We have studied the dependence of 13C polarization on laser and microwave powers. For the first time, a triplet structure corresponding to the 14N hyperfine splitting has been observed in the 13C polarization spectrum. By simultaneously exciting three microwave frequencies at the peaks of the triplet, we have achieved 13C bulk polarization of 0.113 %, leading to an enhancement of 90,000 over the thermal polarization at 17.6 mT. We believe that the multi-tone irradiation can be extended to further enhance the 13C polarization at a low magnetic field.     

Millimeter-Wave Spectroscopy of Methylfuran Isomers: Local vs. Global Treatments of the Internal Rotation

Methylfurans are methylated aromatic heterocyclic volatile organic compounds and primary or secondary pollutants in the atmosphere due to their capability to form secondary organic aerosols in presence of atmospheric oxidants. There is therefore a significant interest to monitor these molecules in the gas phase. High resolution spectroscopic studies of methylated furan compounds are generally limited to pure rotational spectroscopy in the vibrational ground state. This lack of results might be explained by the difficulties arisen from the internal rotation of the methyl group inducing non-trivial patterns in the rotational spectra. In this study, we discuss the benefits to assign the mm-wave rotational-torsional spectra of methylfuran with the global approach of the BELGI-Cs code compared to local approaches such as XIAM and ERHAM. The global approach reproduces the observed rotational lines of 2-methylfuran and 3-methylfuran in the mm-wave region at the experimental accuracy for the ground vt=0 and the first torsional vt=1 states with a unique set of molecular parameters. In addition, the V3 and V6 parameters describing the internal rotation potential barrier may be determined with a high degree of accuracy with the global approach. Finally, a discussion with other heterocyclic compounds enables the study of the influence of the electronic environment on the hindered rotation of the methyl group.     

The Relationship between the IC50 Values and the Apparent Inhibition Constant in the Study of Inhibitors of Tyrosinase Diphenolase Activity Helps Confirm the Mechanism of Inhibition

Tyrosinase is the enzyme involved in melanization and is also responsible for the browning of fruits and vegetables. Control of its activity can be carried out using inhibitors, which is interesting in terms of quantitatively understanding the action of these regulators. In the study of the inhibition of the diphenolase activity of tyrosinase, it is intriguing to know the strength and type of inhibition. The strength is indicated by the value of the inhibition constant(s), and the type can be, in a first approximation: competitive, non-competitive, uncompetitive and mixed. In this work, it is proposed to calculate the degree of inhibition (iD), varying the concentration of inhibitor to a fixed concentration of substrate, L-dopa (D). The non-linear regression adjustment of iD with respect to the initial inhibitor concentration [I]0 allows for the calculation of the inhibitor concentration necessary to inhibit the activity by 50%, at a given substrate concentration (IC₅₀), thus avoiding making interpolations between different values of iD. The analytical expression of the IC₅₀, for the different types of inhibition, are related to the apparent inhibition constant (KIapp). Therefore, this parameter can be used: (a) To classify a series of inhibitors of an enzyme by their power. Determining these values at a fixed substrate concentration, the lower IC₅₀, the more potent the inhibitor. (b) Checking an inhibitor for which the type and the inhibition constant have been determined (using the usual methods), must confirm the IC₅₀ value according to the corresponding analytical expression. (c) The type and strength of an inhibitor can be analysed from the study of the variation in iD and IC₅₀ with substrate concentration. The dependence of IC₅₀ on the substrate concentration allows us to distinguish between non-competitive inhibition (iD does not depend on [D]0) and the rest. In the case of competitive inhibition, this dependence of iD on [D]0 leads to an ambiguity between competitive inhibition and type 1 mixed inhibition. This is solved by adjusting the data to the possible equations; in the case of a competitive inhibitor, the calculation of KI1app is carried out from the IC₅₀ expression. The same occurs with uncompetitive inhibition and type 2 mixed inhibition. The representation of iD vs. n, with n=[D]0/KmD, allows us to distinguish between them. A hyperbolic iD vs. n representation that passes through the origin of coordinates is a characteristic of uncompetitive inhibition; the calculation of KI2app is immediate from the IC₅₀ value. In the case of mixed inhibitors, the values of the apparent inhibition constant of meta-tyrosinase (Em) and oxy-tyrosinase (Eox), KI1app and the apparent inhibition constant of metatyrosinase/Dopa complexes (EmD) and oxytyrosinase/Dopa (EoxD), KI2app are obtained from the dependence of iD vs. n, and the results obtained must comply with the IC₅₀ value.     

Cross-Correlated Motions in Azidolysozyme

The changes in the local and global dynamics of azide-labelled lysozyme compared with that of the wild type protein are quantitatively assessed for all alanine residues along the polypeptide chain. Although attaching -N3 to alanine residues has been considered to be a minimally invasive change in the protein it is found that depending on the location of the alanine residue, the local and global changes in the dynamics differ. For Ala92, the change in the cross-correlated motions are minimal, whereas attaching -N3 to Ala90 leads to pronounced differences in the local and global correlations as quantified by the cross-correlation coefficients of the Cα atoms. We also demonstrate that the spectral region of the asymmetric azide stretch distinguishes between alanine attachment sites, whereas changes in the low frequency, far-infrared region are less characteristic.     

On Some Topological Indices Defined via the Modified Sombor Matrix

Let G be a simple graph with the vertex set V={v1,…,vn} and denote by dvi the degree of the vertex vi. The modified Sombor index of G is the addition of the numbers (dvi2+dvj2)−1/2 over all of the edges vivj of G. The modified Sombor matrix AMS(G) of G is the n by n matrix such that its (i,j)-entry is equal to (dvi2+dvj2)−1/2 when vi and vj are adjacent and 0 otherwise. The modified Sombor spectral radius of G is the largest number among all of the eigenvalues of AMS(G). The sum of the absolute eigenvalues of AMS(G) is known as the modified Sombor energy of G. Two graphs with the same modified Sombor energy are referred to as modified Sombor equienergetic graphs. In this article, several bounds for the modified Sombor index, the modified Sombor spectral radius, and the modified Sombor energy are found, and the corresponding extremal graphs are characterized. By using computer programs (Mathematica and AutographiX), it is found that there exists only one pair of the modified Sombor equienergetic chemical graphs of an order of at most seven. It is proven that the modified Sombor energy of every regular, complete multipartite graph is 2; this result gives a large class of the modified Sombor equienergetic graphs. The (linear, logarithmic, and quadratic) regression analyses of the modified Sombor index and the modified Sombor energy together with their classical versions are also performed for the boiling points of the chemical graphs of an order of at most seven.     

Fluorescent Liquid Tetrazines

Tetrazines with branched alkoxy substituents are liquids at ambient temperature that despite the high chromophore density retain the bright orange fluorescence that is characteristic of this exceptional fluorophore. Here, we study the photophysical properties of a series of alkoxy-tetrazines in solution and as neat liquids. We also correlate the size of the alkoxy substituents with the viscosity of the liquids. We show using time-resolved spectroscopy that intersystem crossing is an important decay pathway competing with fluorescence, and that its rate is higher for 3,6-dialkoxy derivatives than for 3-chloro-6-alkoxytetrazines, explaining the higher fluorescence quantum yields for the latter. Quantum chemical calculations suggest that the difference in rate is due to the activation energy required to distort the tetrazine core such that the nπ*S1 and the higher-lying ππ*T2 states cross, at which point the spin-orbit coupling exceeding 10 cm−1 allows for efficient intersystem crossing to occur. Femtosecond time-resolved anisotropy studies in solution allow us to measure a positive relationship between the alkoxy chain lengths and their rotational correlation times, and studies in the neat liquids show a fast decay of the anisotropy consistent with fast exciton migration in the neat liquid films.     

Data-Driven and Multiscale Modeling of DNA-Templated Dye Aggregates

Dye aggregates are of interest for excitonic applications, including biomedical imaging, organic photovoltaics, and quantum information systems. Dyes with large transition dipole moments (μ) are necessary to optimize coupling within dye aggregates. Extinction coefficients (ε) can be used to determine the μ of dyes, and so dyes with a large ε (>150,000 M⁻¹cm⁻¹) should be engineered or identified. However, dye properties leading to a large ε are not fully understood, and low-throughput methods of dye screening, such as experimental measurements or density functional theory (DFT) calculations, can be time-consuming. In order to screen large datasets of molecules for desirable properties (i.e., large ε and μ), a computational workflow was established using machine learning (ML), DFT, time-dependent (TD-) DFT, and molecular dynamics (MD). ML models were developed through training and validation on a dataset of 8802 dyes using structural features. A Classifier was developed with an accuracy of 97% and a Regressor was constructed with an R2 of above 0.9, comparing between experiment and ML prediction. Using the Regressor, the ε values of over 18,000 dyes were predicted. The top 100 dyes were further screened using DFT and TD-DFT to identify 15 dyes with a μ relative to a reference dye, pentamethine indocyanine dye Cy5. Two benchmark MD simulations were performed on Cy5 and Cy5.5 dimers, and it was found that MD could accurately capture experimental results. The results of this study exhibit that our computational workflow for identifying dyes with a large μ for excitonic applications is effective and can be used as a tool to develop new dyes for excitonic applications.     

An Exchange Mechanism for the Magnetic Behavior of Er3+ Complexes

We study the magnetic properties of the erbium based compounds, Na9[Er(W5O18)2] and [(Pc)Er{Pc{N(C4H9)2}8}]·/−, in the framework of an effective spin exchange model involving delocalized electrons occupying molecular orbitals. The calculations successfully reproduce the experimental data available in the literature for the magnetic spectrum, magnetization and molar susceptibility in dc and ac fields. Owing to their similar molecular geometry, the compounds’ magnetic behaviors are interpreted in terms of the same set of active orbitals and thus the same effective spin coupling scheme. For all three complexes, the model predicts a prompt change in the ground state from a Kramer’s doublet at zero fields to a fully polarized quartet one brought about by the action of an external magnetic field without Zeeman splitting. This alteration is attributed to the enhancement of the effect of orbital interactions over the spin exchange as the magnitude of the external magnetic field increases.     

Naphthazarin Derivatives in the Light of Intra- and Intermolecular Forces

Our long-term investigations have been devoted the characterization of intramolecular hydrogen bonds in cyclic compounds. Our previous work covers naphthazarin, the parent compound of two systems discussed in the current work: 2,3-dimethylnaphthazarin (1) and 2,3-dimethoxy-6-methylnaphthazarin (2). Intramolecular hydrogen bonds and substituent effects in these compounds were analyzed on the basis of Density Functional Theory (DFT), Møller–Plesset second-order perturbation theory (MP2), Coupled Clusters with Singles and Doubles (CCSD) and Car-Parrinello Molecular Dynamics (CPMD). The simulations were carried out in the gas and crystalline phases. The nuclear quantum effects were incorporated a posteriori using the snapshots taken from ab initio trajectories. Further, they were used to solve a vibrational Schrödinger equation. The proton reaction path was studied using B3LYP, ωB97XD and PBE functionals with a 6-311++G(2d,2p) basis set. Two energy minima (deep and shallow) were found, indicating that the proton transfer phenomena could occur in the electronic ground state. Next, the electronic structure and topology were examined in the molecular and proton transferred (PT) forms. The Atoms In Molecules (AIM) theory was employed for this purpose. It was found that the hydrogen bond is stronger in the proton transferred (PT) forms. In order to estimate the dimers’ stabilization and forces responsible for it, the Symmetry-Adapted Perturbation Theory (SAPT) was applied. The energy decomposition revealed that dispersion is the primary factor stabilizing the dimeric forms and crystal structure of both compounds. The CPMD results showed that the proton transfer phenomena occurred in both studied compounds, as well as in both phases. In the case of compound 2, the proton transfer events are more frequent in the solid state, indicating an influence of the environmental effects on the bridged proton dynamics. Finally, the vibrational signatures were computed for both compounds using the CPMD trajectories. The Fourier transformation of the autocorrelation function of atomic velocity was applied to obtain the power spectra. The IR spectra show very broad absorption regions between 700 cm−1–1700 cm−1 and 2300 cm−1–3400 cm−1 in the gas phase and 600 cm−1–1800 cm−1 and 2200 cm−1–3400 cm−1 in the solid state for compound 1. The absorption regions for compound 2 were found as follows: 700 cm−1–1700 cm−1 and 2300 cm−1–3300 cm−1 for the gas phase and one broad absorption region in the solid state between 700 cm−1 and 3100 cm−1. The obtained spectroscopic features confirmed a strong mobility of the bridged protons. The inclusion of nuclear quantum effects showed a stronger delocalization of the bridged protons.     

Galvano-Fenton Engineering Solution with Spontaneous Catalyst’s Generation from Waste: Experimental Efficiency,Parametric Analysis and Modeling Interpretation Applied to a Clean Technology for Dyes Degradation in Water

In this paper, the degradation of the diazo dye naphthol blue black (NBB) using the Galvano-Fenton process is studied experimentally and numerically. The simulations are carried out based on the anodic, cathodic, and 34 elementary reactions evolving in the electrolyte, in addition to the oxidative attack of NBB by HO• at a constant rate of 3.35×107 mol−1·m3·s−1 during the initiation stage of the chain reactions. The selection of the operating conditions including the pH of the electrolyte, the stirring speed, and the electrodes disposition is performed by assessing the kinetics of NBB degradation; these parameters are set to 3, 350 rpm and a parallel disposition with a 3 cm inter-electrode distance, respectively. The kinetics of Fe(III) in the electrolyte were monitored using the principles of Fricke dosimetry and simulated numerically. The model showed more than a 96% correlation with the experimental results in both the blank test and the presence of the dye. The effects of H2O2 and NBB concentrations on the degradation of the dye were examined jointly with the evolution of the simulated H2O2, Fe2+, and HO• concentrations in the electrolyte. The model demonstrated a good correlation with the experimental results in terms of the initial degradation rates, with correlation coefficients exceeding 98%.     

Semiempirical Potential in Kinetics Calculations on the HC3N + CN Reaction

The reaction between the cyano radical CN and cyanoacetylene molecule HC3N is of great interest in different astronomical fields, from star-forming regions to planetary atmospheres. In this work, we present a new synergistic theoretical approach for the derivation of the rate coefficient for gas phase neutral-neutral reactions. Statistic RRKM calculations on the Potential Energy Surface are coupled with a semiempirical analysis of the initial bimolecular interaction. The value of the rate coefficient for the HC3N + CN → H + NCCCCN reaction obtained with this method is compared with previous theoretical and experimental investigations, showing strengths and weaknesses of the new presented approach.