Density Functional Theory (DFT) Study of Edaravone Derivatives as Antioxidants

Quantum chemical calculations at the B3LYP/6-31G* level of theory were employed for the structure-activity relationship and prediction of the antioxidant activity of edaravone and structurally related derivatives using energy (E), ionization potential (IP), bond dissociation energy (BDE), and stabilization energies(Delta E-iso). Spin density calculations were also performed for the proposed antioxidant activity mechanism. The electron abstraction is related to electron-donating groups (EDG) at position 3, decreasing the IP when compared to substitution at position 4. The hydrogen abstraction is related to electron-withdrawing groups (EDG) at position 4, decreasing the BDECH when compared to other substitutions, resulting in a better antioxidant activity. The unpaired electron formed by the hydrogen abstraction from the C-H group of the pyrazole ring is localized at 2, 4, and 6 positions. The highest scavenging activity prediction is related to the lowest contribution at the carbon atom. The likely mechanism is related to hydrogen transfer. It was found that antioxidant activity depends on the presence of EDG at the C-2 and C-4 positions and there is a correlation between IP and BDE. Our results identified three different classes of new derivatives more potent than edaravone.

A theoretical study of SnF2+, SnCl2+, and SnO2+ and their experimental search

We present a detailed theoretical study of the stability of the gas-phase diatomic dications SnF2+, SnCl2+, and SnO2+ using ab initio computer calculations. The ground states of SnF2+, SnCl2+, and SnO2+ are thermodynamically stable, respectively, with dissociation energies of 0.45, 0.30, and 0.42 eV. Whereas SnF2+ dissociates into Sn2+ + F, the long range behaviour of the potential energy curves of SnCl2+ and SnO2+ is repulsive and wide barrier heights due to avoided crossing act as a kind of effective dissociation energy. Their equilibrium internuclear distances are 4.855, 5.201, and 4.852 a(0), respectively. The double ionisation energies (T-e) to form SnF2+, SnCl2+, and SnO2+ from their respective neutral parents are 25.87, 23.71, and 25.97 eV. We combine our theoretical work with the experimental results of a search for these doubly positively charged diatomic molecules in the gas phase. SnO2+ and SnF2+ have been observed for prolonged oxygen (O-16(-)) ion beam sputtering of a tin metal foil and of tin (II) fluoride (SnF2) powder, respectively, for ion flight times of about 10(-5) s through a magnetic-sector mass spectrometer. In addition, SnCl2+ has been detected for O-16(-) ion surface bombardment of stannous (tin (II)) chloride (SnCl2) powder. To our knowledge, SnF2+ is a novel gas-phase molecule, whereas SnCl2+ had been detected previously by electron-impact ionization mass spectrometry, and SnO2+ had been observed before by spark source mass spectrometry as well as by atom probe mass spectrometry. We are not aware of any previous theoretical studies of these molecular systems. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4758475]

On-surface chemical reactions on an insulating substrate

Diese Arbeit beschreibt zum ersten Mal die kovalente Verknüpfung organischer Moleküle auf einer Isolatoroberfläche, motiviert im Hinblick auf die Nutzung der Synthesemethode für die molekulare Elektronik und verwandte Anwendungen. Durch die Verwendung der Nichtkontakt-Rasterkraftmikroskopie und der Kelvinprobe-Mikroskopie bei Raumtemperatur wurden grundlegende molekulare Prozesse der Wechselwirkungen zwischen Molekülen und der Calcit(10.4) Oberfläche sowie die chemische Reaktivität der Moleküle auf der Oberfläche analysiert. Das Zusammenspiel zwischen intermolekularen und Molekül-Oberfläche Wechselwirkungen zeigt sich für Biphenyl-4,4'-dicarbonsäure (BPDCA) durch die Koexistenz zweier unterschiedlicher molekularer Strukturen, die einen Einblick in die treibenden Kräfte der molekularen Selbstorganisation bieten. Die sehr ausgeprägte Reihenstruktur basiert auf der optimalen geometrischen Struktur der BPDCA Moleküle zu den Abmessungen des Substrats, während die zweite Struktur durch Wasserstoffbrücken zwischen den Molekülen gekennzeichnet ist. Der Deprotonierungsvorgang von 2,5-Dihydroxybenzoesäure (DHBA)-Molekülen auf Calcit wird bei Zimmertemperatur gezeigt. Zwei Phasen werden beobachtet, die nach Aufbringen der Moleküle koexistieren. Mit der Zeit geht eine bulk-ähnliche Phase in eine stabile, dicht gepackte Phase über. Der Übergang wird durch Betrachtung des Protonierungszustands der Moleküle erklärt. Die bulk-ähnliche Phase benötigt Wasserstoffbrückbindungen zur Strukturbildung. Werden die Moleküle deprotoniert, so wird die resultierende dicht gepackte Phase durch die elektrostatische Wechselwirkung der deprotonierten Carboxylatgruppen mit den Oberflächen-Calciumkationen stabilisiert. 4-Iodbenzoesäure (IBA)-Moleküle bilden auf Calcit nur Inseln an Stufenkanten, was auf die schwache Molekül-Oberflächen-Wechselwirkung zurückzuführen ist. Für einen stärkeren Einfluss des Substrats durchlaufen die Moleküle einen kontrollierten Übergangsschritt vom protonierten zum deprotonierten Zustand. Im deprotonierten Zustand nehmen die Moleküle eine wohldefinierte Adsorptionsposition auf dem Substrat ein. Die deprotonierte Säuregruppe wird ausgenutzt, um die Desorption der halogensubstituierten Benzoesäure-Moleküle bei der thermischer Aktivierung für die Vernetzungsreaktion zu vermeiden. Darüber hinaus wird die Carboxylatgruppe als starker Elektronendonor verwendet um die Phenyl-Halogen-Bindung zu schwächen und somit die homolytische Spaltung dieser Bindung auch bei moderaten Temperaturen zu ermöglichen. Diesem Konzept folgend ist die erste erfolgreiche kovalente Verknüpfung von 2,5-Diiod-benzoesäure, 2,5-Dichlorbenzoesäure, 3,5-Diiod Salicylsäure und 4-Iod-benzoesäure zu durchkonjugierten molekularen Drähten, Zick-Zack-Strukturen sowie Dimere gezeigt durch Ausnutzen von unterschiedlichen Substitutionsposition sowie Ändern der Anzahl der substituierten Halogenatome. Aufbauend auf diesem Erfolg, wird eine zweistufige Vernetzungsreaktion vorgestellt. Zum Induzieren der ortsspezifischen und sequentiellen kovalenten Verknüpfung wird ein Ausgangsmolekül gewählt, das sowohl eine Bromphenyl als auch eine Chlorphenyl Gruppe mit unterschiedlichen Dissoziationsenergien für die homolytische Spaltung besitzt. Die Reaktionsstellen und sequentielle Reihenfolge für die Reaktion sind somit in der molekularen Struktur einkodiert und bisher unerreichte Reaktionspfade können mithilfe der kovalente Verknüpfung organischer Moleküle auf einer Isolatoroberfläche beschritten werden.

Watson–Crick and Sugar-Edge Base Pairing of Cytosine in the Gas Phase: UV and Infrared Spectra of Cytosine·2-Pyridone

While keto-amino cytosine is the dominant species in aqueous solution, spectroscopic studies in molecular beams and in noble gas matrices show that other cytosine tautomers prevail in apolar environments. Each of these offers two or three H-bonding sites (Watson–Crick, wobble, sugar-edge). The mass- and isomer-specific S1 ← S0 vibronic spectra of cytosine·2-pyridone (Cyt·2PY) and 1-methylcytosine·2PY are measured using UV laser resonant two-photon ionization (R2PI), UV/UV depletion, and IR depletion spectroscopy. The UV spectra of the Watson–Crick and sugar-edge isomers of Cyt·2PY are separated using UV/UV spectral hole-burning. Five different isomers of Cyt·2PY are observed in a supersonic beam. We show that the Watson–Crick and sugar-edge dimers of keto-amino cytosine with 2PY are the most abundant in the beam, although keto-amino-cytosine is only the third most abundant tautomer in the gas phase. We identify the different isomers by combining three different diagnostic tools: (1) methylation of the cytosine N1–H group prevents formation of both the sugar-edge and wobble isomers and gives the Watson–Crick isomer exclusively. (2) The calculated ground state binding and dissociation energies, relative gas-phase abundances, excitation and the ionization energies are in agreement with the assignment of the dominant Cyt·2PY isomers to the Watson–Crick and sugar-edge complexes of keto-amino cytosine. (3) The comparison of calculated ground state vibrational frequencies to the experimental IR spectra in the carbonyl stretch and NH/OH/CH stretch ranges strengthen this identification.

Gibbs energies of formation of rare earth oxysulfides

The standard Gibbs energy change accompanying the conversion of rare earth oxides to oxysulfides by reaction of rare earth oxides with diatomic sulfur gas has been measured in the temperature range 870 to 1300 K using the solid state cell: Pt/Cu+Cu2S/R2O2S+R2O3‖(CaO)ZrO2‖Ni+NiO, Pt where R=La, Nd, Sm, Gd, Tb, and Dy. The partial pressure of diatomic sulfur over a mixture of rare earth oxide (R2O3) and oxysulfide (R2O2S) is fixed by the dissociation of Cu2S to Cu in a closed system. The buffer mixture of Cu+Cu2S is physically separated from the rare earth oxide and oxysulfide to avoid complications arising from interaction between them. The corresponding equilibrium oxygen partial pressure is measured with an oxide solid electrolyte cell. Gibbs energy change for the conversion of oxide to the corresponding oxysulfide increases monotonically with atomic number of the rare earth element. Second law enthalpy of formation also shows a similar trend. Based on this empirical trend Gibbs energies of formation of oxysulfides of Pr, Eu, Ho, and Er are estimated as a function of temperature.

Gibbs’ free energies of formation of Cu2Ln2O5 (Ln = Tb, Dy, Er, Yb) compounds

The standard Gibbs' free energies of formation of compounds of type Cu2L%05 (Ln = Tb,Dy,Er,Yb) were measured using the solid state cell in the temperature range of 970 to 1323 K For formation of Cu2L?O5 compounds from their binary component oxides according to the reaction 2 CUO (s) + L%03 (s) -, Cu,,L%05 (s),the Gibbs' free energy changes can be represented by the following equations:AGO = 13 080 - 13.70 'I" (+80) J mol-' (Ln = Tb)AGq = 11 480 - 13.51 T (260) J mol-I (Ln = Dy)AGO = 10 750 - 13.99 T (260) J mol-I (Ln = Er)AGO = 9 920 - 13.90 T (260) J mol-' (Ln = Yb) Since formation of the compounds is endothermic, the compounds become thermodynamically unstable with respect to their component oxides below 955 K for Cu2Tb205, 850 K for Cu2Dy205, 768 K for Cu2Er205 and 714 K for Cu2Yb2OS When the oxygen partial pressure over Cu2L%05 is lowered, they decompose according to the scheme, 2 CU,L%O, (s) -r 2 L%03 (s) +2 cu20 (s) + 02(g)The equilibrium chemical potentials of oxygen corresponding to the dissociation reactions are computed from the emf data and auxiliary information on Cu20 and CuO. The computed decomposition temperatures at an oxygen partial pressure of 5.0 x ld Pa are compared with those obtained directly from combined thermogravimetric (TGA) and differential thermal analyses (DTA).The free energy, enthalpy and entropy of formation of Cu2Ln205 compounds show systematic variation with the ionic radius of the trivalent lanthanide ion. The trends obtained in this study are compared with information available in the literature. The staZbility of Cu2Ln205 compounds increases with the decrease in ionic radii of the ~ n ion~. +

Intense-field dissociation dynamics of D2+ molecular ions using ultrafast laser pulses

The dynamics of dissociation of pre-ionized D2+ molecules using intense (10^12–10^15 W cm-2), ultrashort (50 fs), infrared (? = 790 nm) laser pulses are examined. Use of an intensity selective scan technique has allowed the deuterium energy spectrum to be measured over a broad range of intensity. It is found that the dominant emission shifts to lower energies as intensity is increased, in good agreement with corresponding wavepacket simulations. The results are consistent with an interpretation in terms of bond softening, which at high intensity (approximately >3 × 10^14 W cm-2) becomes dominated by dissociative ionization. Angular distribution measurements reveal the presence of slow molecular dissociation, an indication that vibrational trapping mechanisms occur in this molecule.

Insight into the adsorption competition and the relationship between dissociation and association reactions in ammonia synthesis

Ammonia synthesis on three metal surfaces (Zr, Ru, and Pd) is investigated using density functional theory calculations. In addition to N-2 dissociation, all the transition states of the hydrogenation reactions from N to NH3 are located and the reaction energy profiles at both low and high surface coverages are compared and analyzed. The following are found: (i) Surface coverage effect on dissociation reactions is more significant than that on association reactions. (ii) The difference between N and H chemisorption energies, the so-called chemisorption energy gap which is a measure of adsorption competition, is vital to the reactivity of the catalysts. (iii) The hydrogenation barriers can considerably affect the overall rate of ammonia synthesis. A simple model to describe the relationship between dissociation and association reactions is proposed. (c) 2007 American Institute of Physics.

General trends in CO dissociation on transition metal surfaces

Dissociative adsorption is one of the most important reactions in catalysis. In this communication we propose a model aiming to generalize the important factors that affect dissociation reactions. Specifically, for a dissociation reaction, say AB -->A + B, the model connects the dissociation barrier with the association barrier, the chemisorption energies of A and B at the final state and the bonding energy of AB in the gas phase. To apply this model, we have calculated CO dissociation on Ru(0001), Rh(111), Pd(111) (4d transition metals), Os(0001), Ir(111), and Pt(111) (5d transition metals) using density function theory (DFT). All the barriers are determined. We find that the DFT results can be rationalized within the model. The model can also be used to explain many experimental observations. (C) 2001 American Institute of Physics.

Can mass dissociation patterns of transition-metal complexes be predicted from electrochemical data?

The Cooks kinetic method has been very convenient to correlate the relative dissociation rates obtained by collision-induced fragmentation experiments with the energies of two related bonds in molecules and complexes in the gas phase. Reliable bond energy data are, however, not always available, particularly for polynuclear transition-metal complexes, such as the triruthenium acetate clusters of the general formula [Ru(3) (mu(3)-O)(mu-CH(3)COO)(6)(py)(2)(L)](+), where L = ring substituted N-heterocyclic ligands. Accordingly, their gas-phase collision-induced tandem mass spectrometry (CID MS/MS) dissociation patterns have been analyzed pursuing a relationship with the more easily accessible redox potentials (E(1/2)) and Lever`s E(L) parameters. In fact, excellent linear correlations of In(1/2A(L)/A(py)), where A(py) and A(L) are the abundance of the fragments retaining the pyridine (py) and L ligand, respectively, with E(1/2) and E(L) were found. This result shows that those electrochemical parameters are correlated with bond energies and can be used in the analysis of the dissociation data. Such modified Cooks method can be used, for example, to determine the electronic effects of substituents on the metal-ligand bonds for a series of transition-metal complexes. Copyright (C) 2008 John Wiley & Sons, Ltd.

Role of the range of the dipole function in the classical dynamics of molecular dissociation

The dissociation dynamics of heteronuclear diatomic molecules induced by infrared laser pulses is investigated within the framework of the classical driven Morse oscillator. The interaction between the molecule and the laser field described in the dipole formulation is given by the product of a time-dependent external field with a position-dependent permanent dipole function. The effects of changing the spatial range of the dipole function in the classical dissociation dynamics of large ensembles of trajectories are studied. Numerical calculations have been performed for distinct amplitudes and carrier frequencies of the external pulses and also for ensembles with different initial energies. It is found that there exist a set of values of the dipole range for which the dissociation probability can be completely suppressed. The dependence of the dissociation on the dipole range is explained through the examination of the Fourier series coefficients of the dipole function in the angle variable of the free system. In particular, the suppression of dissociation corresponds to dipole ranges for which the Fourier coefficients associated with nonlinear resonances are null and the chaotic region in the phase space is reduced to thin layers. In this context, it is shown that the suppression of dissociation of heteronuclear molecules for certain frequencies of the external field is a consequence of the finite range of the corresponding permanent dipole. © 2013 American Physical Society.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin

NMR investigations have been carried out of complexes between bovine chymotrypsin Aα and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-l-Phe-CF3, N-Acetyl-Gly-l-Phe-CF3, N-Acetyl-l-Val-l-Phe-CF3, and N-Acetyl-l-Leu-l-Phe-CF3. The D/H fractionation factors (φ) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hδ1) in these four complexes at 5°C were in the range φ = 0.32–0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hɛ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 8.97–9.18), the exchange rate of the His 57-Hδ1 proton with bulk water protons (284–12.4 s−1), and the activation enthalpies for this hydrogen exchange (14.7–19.4 kcal⋅mol−1). It was found that the previously noted correlations between the inhibition constants (Ki 170–1.2 μM) and the chemical shifts of His 57-Hδ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 18.61–18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940–11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hδ1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

Unbinding forces of single antibody-antigen complexes correlate with their thermal dissociation rates

Point mutants of three unrelated antifluorescein antibodies were constructed to obtain nine different single-chain Fv fragments, whose on-rates, off-rates, and equilibrium binding affinities were determined in solution. Additionally, activation energies for unbinding were estimated from the temperature dependence of the off-rate in solution. Loading rate-dependent unbinding forces were determined for single molecules by atomic force microscopy, which extrapolated at zero force to a value close to the off-rate measured in solution, without any indication for multiple transition states. The measured unbinding forces of all nine mutants correlated well with the off-rate in solution, but not with the temperature dependence of the reaction, indicating that the same transition state must be crossed in spontaneous and forced unbinding and that the unbinding path under load cannot be too different from the one at zero force. The distance of the transition state from the ground state along the unbinding pathway is directly proportional to the barrier height, regardless of the details of the binding site, which most likely reflects the elasticity of the protein in the unbinding process. Atomic force microscopy thus can be a valuable tool for the characterization of solution properties of protein-ligand systems at the single molecule level, predicting relative off-rates, potentially of great value for combinatorial chemistry and biology.