Solution structure, mutagenesis, and NH exchange studies of the MutT enzyme–Mg-8-oxo-dGMP complex

The MutT pyrophosphohydrolase from E. coli (129 residues) catalyzes the hydrolysis of nucleoside triphosphates (NTP), including 8-oxo-dGTP, by substitution at Pβ, to yield NMP and pyrophosphate. The product, 8-oxo-dGMP is an unusually tight binding, slowly exchanging inhibitor with a K_D=52 nM, (ΔG°=−9.8 kcal/mol) which is 6.1 kcal/mol tighter than the binding of dGMP (ΔG°=−3.7 kcal/mol). The higher affinity for 8-oxo-dGMP results from a more favorable ΔH_binding (−32 kcal/mol) despite an unfavorable −TΔS°_binding (+22 kcal/mol). The solution structure of the MutT–Mg²⁺-8-oxo-dGMP complex shows a narrowed, hydrophobic nucleotide-binding cleft with Asn-119 and Arg-78 among the few polar residues. The N119A, N119D, R78K and R78A single mutations, and the R78K+N119A double mutant all showed largely intact active sites, on the basis of small changes in the kinetic parameters of dGTP hydrolysis and in ¹H–¹⁵N HSQC spectra. However, the N119A mutation profoundly weakened the active site binding of 8-oxo-dGMP by 4.3 kcal/mol (1650-fold). The N119D mutation also weakened 8-oxo-dGMP binding but only by 2.1 kcal/mol (37-fold), suggesting that Asn-119 functioned both as a hydrogen bond donor to C8=O, and a hydrogen bond acceptor from N7H of 8-oxo-dGMP, while aspartate at position −119 functioned as an acceptor of a single hydrogen bond. Much smaller weakening effects (0.3–0.4 kcal/mol) on the binding of dGMP and dAMP were found, indicating specific hydrogen bonding of Asn-119 to 8-oxo-dGMP. While formation of the wild type MutT–Mg²⁺-8-oxo-dGMP complex slowed the backbone NH exchange rates of 45 residues distributed throughout the protein, the same complex of the N119A mutant slowed the exchange rates of only 11 residues at or near the active site, indicating an increase in conformational flexibility of the N119A mutant. The R78K and R78A mutations weakened the binding of 8-oxo-dGMP by 1.7 and 1.1 kcal/mol, respectively, indicating a lesser role of Arg-78 than of Asn-119 in the selective binding of 8-oxo-dGMP, likely donating a single hydrogen bond to its C6=O. The R78K+N119A double mutant weakened the binding of 8-oxo-dGMP (K_I^slope=3.1 mM) by 6.5±0.2 kcal/mol which overlaps, within error with the sum of the effects of the two single mutants (6.0±0.3 kcal/mol). Such additive effects of the two single mutants in the double mutant are most simply explained by the independent functioning of Asn-119 and Arg-78 in the binding of 8-oxo-dGMP. Independent functioning of these two residues in nucleotide binding is consistent with their locations in the MutT–Mg²⁺-8-oxo-dGMP complex, on opposite sides of the active site cleft, with a distance of 8.4±0.5 Å between their side chain nitrogens.     

DFT and CASPT2 study of two thermal reactions of nitromethane: C–N bond cleavage and nitro-to-nitrite isomerization. An example of the inverse symmetry breaking deficiency in density functional calculations of an homolytic dissociation

The reaction paths of nitromethane leading to the dissociation products or isomerization to methyl nitrite have been computationally investigated at the CAS-SCF and DFT levels of theory. Additionally, the CAS-SCF wave functions were used as reference in a second-order perturbation treatment, CASPT2, in order to obtain a good estimate for the activation energy of each reaction path. Both methods predict the isomerization as a concerted reaction. However, the behavior of the two approximations with respect to dissociation is rather different; while CASPT2 predicts a barrier height of (≈59 kcal/mol) in good accordance with the experimental activation energy (59.0 kcal/mol), B3-LYP/6-31G^* calculations overestimate the barrier for more than 30 kcal/mol. The DFT prediction of the dissociation channel exhibits inverse symmetry breaking, dissociating to the unphysical absurd CH₃^δ+ plus NO₂^δ−.     

Solvent effect on K to Na ion mutation: a Monte Carlo simulation study

Monte Carlo simulation studies of statistical perturbation theory (SPT) have been carried out to investigate the solvent effects on the relative free energies of solvation and the difference in partition coefficients (log P) for K⁺ to Na⁺ ion mutation in the several solvents. We compared the relative free energies for interconversion of K⁺ to Na⁺, in H₂O (TIP4P) in this study with those published works, that in H₂O (TIP4P) is −16.55 kcal/mol in this study, those of the published works are −17.6, −17.3 and −17.31 kcal/mol and that of the experiment is −17.6 kcal/mol, respectively. Comparing the relative free energies for interconversion of K⁺ to Na⁺, in CH₃OH in this study with those published works, that in CH₃OH is −18.08±0.28 kcal/mol in this study, that of molecular dynamic simulation is −19.6±0.4 kcal/mol and that of the experimental work is −17.3 kcal/mol, respectively. There is good agreement among the several studies if we consider both methods of obtaining the solvation (or hydration) free energies and the standard deviations. For the present K⁺ and Na⁺ ions, the relative free energies of solvation vs Born's function of solvents are decreased with increasing Born's function of solvent except for CH₃OH, THF and MEOME. There is also good agreement between the calculated structural properties in this study and the computer simulation, ab initio and experimental works.     

An ab initio study of the Cl(P)+C2H6→C2H5+HCl abstraction reaction

Electronic energies, geometries, and harmonic vibration frequencies for the reactants, products, and transition state for the Cl(³P)+C₂H₆→C₂H₅+HCl abstraction reaction were evaluated at the HF and MP2 levels using several correlation consistent polarized-valence basis sets. Single-point calculations at PMP2, MP4, QCISD(T), and CCSD(T) levels were also carried out. The values of the forward activation energies obtained at the MP4/cc-pVTZ, QCISD(T)/cc-pVTZ, and CCSD(T)/cc-pVTZ levels using the MP2/cc-pVTZ structures are equal to −0.1, −0.4, and −0.3 kcal/mol, respectively. The experimental value is equal to 0.3±0.2 kcal/mol. We found that the MP2/aug-cc-pVTZ adiabatic vibration energy for the reaction (−2.4 kcal/mol) agrees well with the experimental value −(2.2–2.6) kcal/mol. Rate constants calculated with the zeroth-order interpolated variational transition state (IVTST-0) method are in good agreement with experiment. In general, the theoretical rate constants differ from experiment by, at most, a factor of 2.6.     

Dual binding modes of Congo red to amyloid protofibril surface observed in molecular dynamics simulations

Congo red has been used to identify amyloid fibrils in tissues for more than 80 years and is also a weak inhibitor to both amyloid-beta fibril formation and toxicity. However, the specificity of the binding and its inhibition mechanism remain unclear. Using all-atom molecular dynamics simulations with the explicit solvent model, we have identified and characterized two specific binding modes of Congo red molecules to a protofibril formed by an amyloidogenic fragment (GNNQQNY) of the yeast prion protein Sup35. The observation of dual-mode was consistent with the experimentally observed dual-mode binding to Abeta fibrils by a series of compounds similar to Congo red. In the primary mode, Congo red bound to a regular groove formed by the first three residues (GNN) of the beta-strands along the beta-sheet extension direction. Comparative simulations demonstrated that Thioflavin T also bound to the grooves on KLVFFAE protofibril surface. Because of the ubiquitous long grooves on the amyloid fibril surface, we propose that this binding interaction could be a general recognition mode of amyloid fibrils by Congo red, Thioflavin T, and other long flat molecules. In the secondary mode, Congo red bound parallel to the beta-strands on the edge or in the middle of a beta-sheet. The primary binding mode of Congo red and GNNQQNY protofibril was more stable than the secondary mode by -5.7 kcal/mol as estimated by the MM-GBSA method. Detailed analysis suggests that the hydrophobic interactions play important roles for burial of the hydrophobic part of the Congo red molecules. Two potential inhibition mechanisms of disrupting beta-sheet stacking were inferred from the primary mode, which could be exploited for the development of non-peptidic amyloid-specific inhibitors.     

Structure and stability of thioureate anions with water, DFT calculations

Structure and stability of thioureate anion with the water clusters CSNH₂NH⁻-(H₂O)_n = 1–7, are studied, using density functional theory. Molecular structures and the stability of the clusters are discussed based on the calculation results of the stable conformers and their relative energies. All the clusters are stable thermodynamically in gas phase with respect to separate monomers. The clusters are stabilized progressively with an increasing number of water molecules, as indicated by the increasing of the binding energies. The binding energies of CSNH₂NH⁻ and a water molecule are 14.34 and 16.36 kcal/mol for cis CSNH₂NH⁻ and trans CSNH₂NH⁻, respectively. As the reaction in aqueous solution progresses, the CS bond distance increases monotonically, indicating that the CS bond of the thioureate anion unit in the clusters is de-stabilized with an increasing number of water molecules.     

Ab initio calculations of minimum-energy pathways of the nucleophilic addition of the H anion, LiH molecule and Li/H ion pair to acetylene and methylacetylene

Ab initio calculations were performed for special points of the minimal energy pathways (MEP) of the nucleophilic addition reactions of the isolated H⁻ anion, LiH molecule and Li⁺/H⁻ ion pair to acetylene (A) and methylacetylene (MA) molecules, proceeding in accordance (M) and against (aM) the Markovnikov's rule. All structural parameters were optimized using the restricted Hartree–Fock (RHF) method. For the addition of H⁻, the 6-31++G^* basis set was used and for the reactions of LiH and Li⁺/H⁻ the 6-31G^* basis set with the subsequent recalculation of single point energies, taking into account of electron correlation energy by means of the second-order Möller–Plesset perturbation theory at the MP2/6-31++G^** level. The results of calculations demonstrate, that the energy characteristics of both M- and aM-additions with H⁻ do not differ sufficiently (0.1–1.2 kcal/mol for the activation energies (ΔE_a) and the reaction heats (ΔQ)). The substitution of the H atom by the CH₃ group in A molecule results in practically the same values of ΔQ and ΔE_a. On the contrary, for the LiH molecule and Li⁺/H⁻ ionic pair, the M-addition is favorable (charge control). It is found that the presence of electrophile decreases the activation energy by 3–5 kcal/mol as compared with the addition of the isolated hydride ion H⁻.     

Kinetic modeling of promotion and inhibition of temperature on photocatalytic degradation of benzene vapor

This study investigated the effects of temperature, humidity, and benzene concentration on the photocatalytic oxidation of benzene vapor over titanium dioxide. An annular packed-bed photocatalytic reactor was employed to determine the intrinsic oxidation rates for the photocatalysis of benzene. Degussa P-25 TiO₂ was used as the photocatalyst and a 15 W near-UV lamp (350 nm) was used as the light source. The experiments were conducted at influent benzene concentrations of 250–450 ppmv, water vapor concentrations of 13,500–27,500 ppmv, and reaction temperatures ranging from 100 to 200 °C. Benzene oxidation rates increased with temperature below 160–180 °C, but decreased with temperature above 160–180 °C. Raising the reaction temperature increased the chemical reaction rates but reduced the reactant adsorption rate on TiO₂ surfaces. The overall reaction rate increased with temperature, indicating that the reduction of reactant adsorption rate did not affect the overall reaction, and thus the chemical reaction was the rate-limiting step. As the chemical reaction rate gradually exceeds the reactant adsorption rate with temperature, the rate-limiting step was shifted from the chemical reaction to the reactant adsorption. Additionally, the competitive adsorption between benzene and water for the active sites on TiO₂ resulted in the promotion and inhibition of reaction rate by humidity. This study developed a modified bimolecular Langmuir–Hinshelwood kinetic model to simulate the temperature and humidity related promotion and inhibition of the photocatalysis of benzene. The correlation developed here was used as a basis for determining the apparent activation energy of 0.76 kcal/mol and adsorption enthalpies of benzene and water of −20.1 and −13.7 kcal/mol.     

Far infrared spectra, conformational equilibria, vibrational assignments, ab initio calculations and structural parameters for 2-bromoethanol

The far infrared spectrum from 370 to 50 cm⁻¹ of gaseous 2-bromoethanol, BrCH₂CH₂OH, was recorded at a resolution of 0.10 cm⁻¹. The fundamental O–H torsion of the more stable gauche (Gg′) conformer, where the capital G refers to internal rotation around the C–C bond and the lower case g to the internal rotation around the C–O bond, was observed as a series of Q-branch transitions beginning at 340 cm⁻¹. The corresponding O–H torsional modes were observed for two of the other high energy conformers, Tg (285 cm⁻¹) and Tt (234 cm⁻¹). The heavy atom asymmetric torsion (rotation around C–C bond) for the Gg′ conformer has been observed at 140 cm⁻¹. Variable temperature (−63 to −100°C) studies of the infrared spectra (4000–400 cm⁻¹) of the sample dissolved in liquid xenon have been recorded. From these data the enthalpy differences have been determined to be 411±40 cm⁻¹ (4.92±0.48 kJ/mol) for the Gg′/Tt and 315±40 cm⁻¹ (3.76±0.48 kJ/mol) for the Gg′/Tg, with the Gg′ conformer the most stable form. Additionally, the infrared spectrum of the gas, and Raman spectrum of the liquid phase are reported. The structural parameters, conformational stabilities, barriers to internal rotation and fundamental frequencies have been obtained from ab initio calculations utilizing different basis sets at the restricted Hartree–Fock or with full electron correlation by the perturbation method to second order. The theoretical results are compared to the experimental results when appropriate. Combining the ab initio calculations with the microwave rotational constants, r₀ adjusted parameters have been obtained for the three 2-haloethanols (F, Cl and Br) for the Gg′ conformers.     

Flow injection hydride generation electrothermal atomic absorption spectrometric determination of toxicologically relevant arsenic in urine

Analytical procedure for the determination of toxicologically relevant arsenic (the sum of arsenite, arsenate, monomethylarsonate and dimethylarsinate) in urine by flow injection hydride generation and collection of generated inorganic and methylated hydrides on an integrated platform of a transverse-heated graphite atomizer for electrothermal atomic absorption spectrometric determination (ETAAS) is elaborated. Platforms are pre-treated with 2.7 μmol of zirconium and then with 0.10 μmol of iridium which serve both as an efficient hydride sequestration medium and permanent chemical modifier. Arsine, monomethylarsine and dimethylarsine are generated from diluted urine samples (10–25-fold) in the presence of 50 mmol L⁻¹ hydrochloric acid and 70 mmol L⁻¹ l-cysteine. Collection, pyrolysis and atomization temperatures are 450, 500, 2100 and 2150 °C, respectively. The characteristic mass, characteristic concentration and limit of detection (3σ) are 39 pg, 0.078 μg L⁻¹ and 0.038 μg L⁻¹ As, respectively. The limits of detection in urine are ca. 0.4 and 1 μg L⁻¹ with 10- and 25-fold dilutions. The sample throughput rate is 25 h⁻¹. Applications to several urine CRMs are given.     

Water clusters in a confined nonpolar environment

The stability and structure of water clusters in a confined nonpolar environment is investigated theoretically by examining the encapsulation of water molecules inside a fullerene (C₆₀) cage. While the Hartree–Fock (6-31G) calculations suggest H₂O@C₆₀ to be marginally more stable (−0.5 kcal/mol) than the isolated water and C₆₀ molecules, second order Møller–Plesset perturbation theory suggests it to be much more stable (−9.9 kcal/mol). It is shown that encapsulation results in the breaking of hydrogen bonds and rearrangement of water clusters. The tetramer inside the cage, for example, is tetrahedral in arrangement, in contrast to a square planar geometry observed in the gas phase.     

The molecular structure and conformation of trichloronitromethane as determined by gas-phase electron diffraction and theoretical calculations

The molecular structure of trichloronitromethane has been studied in the gas phase using electron diffraction data. The molecules are found to undergo low barrier rotation about the CN bond with a planar CNO₂ moiety in agreement with HF/MP2/B3LYP/6-311G(d,p) calculations. The experimental data are consistent with a dynamic model using a potential function for the torsion of V = (V₆/2)(1 − cos 6τ). The major geometrical parameters (r_g and ) for the eclipsed form, obtained from least squares analysis of the data are as follows: r(NO₃) = r(NO₄) = 1.213(2) Å, r(CN) = 1.592(6) Å, r(CCl)_av = 1.749(1) Å, Cl₅CN/Cl₆CN = 109. 6°/106.3°(2), O₃NC/O₄NC = 117. 6°/114.1°(4), τCl₅C₁N₂O₃ = 0.0°, and V₆ = 0.20(25) kcal/mol.     

Co-aggregation of fullerene C60 and thiophene in the non-aromatic solvent cyclohexene

Co-aggregation of fullerene C₆₀ and thiophene has been studied calorimetrically in cyclohexene at 25 °C. The total aggregation heat is found to depend on initial concentration of thiophene and fall between −1.9 and −5.8 kJ mol⁻¹. The corresponding thiophene/fullerene molar ratio (“co-aggregation number”) ranges from 7 to 12. The data are rationalized by formation of heteromolecular nanoaggregates with intermolecular contacts of both fullerene–thiophene and fullerene–fullerene types. A physical model describing interaction between fullerenes and π-donors in solution is substantiated and used to explain heterogeneity of composites containing fullerenes.     

Infrared and Raman spectra, conformational stability, ab initio calculations of structure, and vibrational assignment of 2-hexyne

The infrared spectra (3500–50 cm⁻¹) of the gas and solid and the Raman spectra (3500–50 cm⁻¹) of the liquid and solid have been recorded for 2-hexyne, CH₃–CC–CH₂CH₂CH₃. Variable temperature studies of the infrared spectrum (3500–400 cm⁻¹) of 2-hexyne dissolved in liquid krypton have also been recorded. Utilizing four anti/gauche conformer pairs, the anti(trans) conformer is found to be the lower energy form with an enthalpy difference of 74±8 cm⁻¹ (0.88±0.10 kJ/mol) determined from krypton solutions over the temperature range −105 to −150 °C. At room temperature it is estimated that there is 42% of the anti conformer present. Equilibrium geometries and energies of the two conformers have been determined by ab initio (HF and MP2) and hybrid DFT (B3LYP) methods using a number of basis sets. Only the HF and DFT methods predict the anti conformer as the more stable form as found experimentally. A vibrational assignment is proposed based on the force constants, relative intensities, depolarization ratios from the ab initio and DFT calculations and on rotational band contours obtained using the calculated equilibrium geometries. From calculated energies it is shown that the CH₃ group exhibits almost completely free rotation which is in agreement with the observation of sub-band structure for the degenerate methyl vibrations from which values of the Coriolis coupling constants, ζ, have been determined. The results are compared to similar properties of some corresponding molecules.     

Kinetic study of trifluoromethylation with s-(trifluoromethyl) dibenzothiophenium salts

The kinetic parameters were determined for C-trifluoromethylation of aniline with S-(trifluoromethyl)dibenzothiophenium triflate (1), its 3,7-dinitro derivative (2) and S-(trifluoromethyl)diphenylsulfonium triflate (3) in DMF-d₇. The higher reactivity of heterocyclic 1 compared with non-heterocyclic 3 could be explained on the basis of its greatly enhanced activation entropy (ΔS^≠: −11.2 cal mol ⁻¹ K⁻¹ for 1; −47.1 for 3), but not its enhanced activation enthalpy (ΔH^≠: 21.2 kcal mol⁻¹ for 1; 12.1 for 3). The aromatic delocalization of the heterocyclic ring may thus be only slightly disturbed by the S-trifluoromethyl substituent. The high reactivity of 2 was attributed to the great electron deficiency caused by two nitro groups in addition to the heterocyclic salt system (ΔH^≠ 17.0 kcal mol⁻¹, ΔS^≠ −9.1 cal mol⁻¹ K⁻¹ for 2). The reaction mechanism is discussed; the conventional S_N2 attack mechanism was ruled out and a mechanism for a side-on attack to the CF₃-S bond may possibly be applicable.     

Infrared and Raman spectra, conformational stability, ab initio calculations, and vibrational assignments for methylaminothiophosphoryl difluoride

Infrared spectra (4000–50 cm⁻¹) of the vapor, amorphous and crystalline solids and Raman spectra (3600–10 cm⁻¹) of the liquid with qualitative depolarization data as well as the amorphous and crystalline solids of methylaminothiophosphoryl difluoride, CH₃N(H)P(=S)F₂, and three deuterated species, CD₃N(H)P(=S)F₂, CH₃N(D)P(=S)F₂, and CD₃N(D)P(=S)F₂, have been recorded. The spectra indicate that in the vapor, liquid and amorphous solid a small amount of a second conformer is present, whereas only one conformer remains in the low temperature crystalline phase. The near-infrared spectra of the vapor confirms the existence of two conformers in the gas phase. Asymmetric top contour simulation of the vapor shows that the trans conformer is the predominant vapor phase conformer. From a temperature study of the Raman spectrum of the liquid the enthalpy difference between the trans and near-cis conformers was determined to be 368±15 cm⁻¹ (4.41±0.2 kJ/mol), with the trans conformer being thermodynamically preferred. Ab Initio calculations with structure optimization using the 6-31G(d) and 6-311+G(d,p) basis sets at the restricted Hartree–Fock (RHF) and/or with full electron correlation by the perturbation method to second order (MP2) support the occurrence of near-trans (5° from trans) and near-cis (20° from cis) conformers. From the RHF/6-31G(d) calculation the near-trans conformer is predicted to be the more stable form by 451 cm⁻¹ (5.35 kJ/mol) and from the MP2/6-311+G(d,p) calculation by 387 cm⁻¹ (4.63 kJ/mol). All of the normal modes of the near-trans rotamer have been assigned based on infrared band contours, depolarization values and group frequencies and the assignment is supported by the normal coordinate calculation utilizing harmonic force constants from the MP2/6-31G(d) ab initio calculations.     

Theoretical calculations of energetics, structures, and rate constants for the H + CH3OH hydrogen abstraction reactions

Barrier heights, structures, reaction energies, and rate constants are calculated with the DFT, MP2, and CCSD(T) methods for the first two channels of the H + CH₃OH reactions. The best estimate results based on CCSD(T) calculations give reaction enthalpies () for the first (−8.4 kcal/mol) and second (0.8 kcal/mol) reactive channels, which are comparable to the experimental values, −8.8 ± 0.9 and −0.3 ± 0.9 kcal/mol, respectively. Rate constants and activation energies calculated with the IVTST-0 method using CCSD(T)/cc-pVTZ geometries and frequencies are also in excellent agreement with experiment data.     

Role of axial ligand on the electronic structures of active intermediates in cytochrome P-450, peroxidase and catalase

DFT method (B3LYP) with 6-31G^* basis set was utilized in the computation of a fully optimized structure, net atomic charges and spin densities of the intermediate of cytochrome P-450-oxoiron(IV) porphyrin cation radical, compound I – in the presence of axial ligand such as thiolate (SMe⁻) imidazole (IM), phenoxide (OPh⁻), methoxide (OMe⁻) and chloride (Cl⁻). The results show doublet states in compound I are about 2–4 kcal/mol more stable than quartet states for all aforementioned ligands, and the doublet state is the ground state in all cases. However, electron donor ability of the ligands are in the order of SMe⁻ > IM > OMe⁻ > OPh⁻ > Cl⁻. Also the active oxidant intermediate of cytochrome P-450 between different mesomeric structures select sulfur oxygen radical type structure and can be viewed as (RS^↓)Fe^↑(IV)(O^↑)(Por). In horseraddish peroxidase (HRP) and peroxidase with histidine axial ligand π cation radical character of porphyrin ring is preferred (Im)Fe^↑(IV)(O^↑)(Por^↓). For the ligands such as OMe⁻, OPh⁻ and Cl⁻ oxidation mainly took place on the iron and the active intermediate can be viewed as (L)Fe^↑(V)(O)(Por) with one unpaired electron localized on the iron.     

Isothermal kinetic analysis of the thermal decomposition of magnesium hydroxide using thermogravimetric data

In the present study, the kinetics of the thermal decomposition of magnesium hydroxide is investigated, using isothermal methods of kinetic analysis. For this purpose, experiments in thermogravimetric analyser were carried out in standard values of temperature (350°, 400°, 450° and 500°C) which resulted in weight loss percent as a function of time. The data were further modified to give fraction reacted ‘' versus time to be tested in various forms of ‘' functions. In order to determine the mechanism of the magnesium hydroxide decomposition and the form of the conversion function which governs the dehydroxylation of Mg(OH)₂, four different methods of isothermal kinetic analysis were used. Applying each of these methods to the data, it was concluded that the nucleation mechanism predominates the Mg(OH)₂, decomposition for all values of temperature tested; at 350°C the kinetic model which represents the experimental data is that of reaction at phase boundaries (random nucleation), F₁: ln(1−)=kt) while for the higher temperatures 400°, 450° and 500°C the kinetic equation of nucleation and development in two dimensions, A₂: [−ln (1−)]^1/2=kt was found to fit better the experimental results. The activation energy was evaluated applying two alternative methods; the Arrhenius plot, using maximum rates of reaction, from which the activation energy was evaluated to be 20.54 kcal/mol. An alternative method based on plots of ln t versus 1/T corresponding to the same value of ‘' gave values of 10.72, 13.82 and 16.31 kcal/mol for ‘' values of 0.25, 0.50 and 0.75, respectively.     

Ab initio studies of the conformations of ethynyl formate, cyano formate and related ethynyl and cyano compounds

The optimized geometries, relative free energies and related thermodynamic properties, harmonic frequencies, and dipole moments have been calculated at the HF and MP2 levels for ethynyl formate (1a), ethynyl acetate (1b), cyano formate, HCO₂CN (1c), cyano acetate (1d), S-ethynyl thioformate (2a), S-ethynyl thioacetate (2b), S-cyano thioformate (2c), S-cyano thioacetate (2d), N-ethynylformamide (3a), N-ethynylacetamide (3b), N-cyanoformamide (3c), and N-cyanoacetamide (3d) with the gaussian 98 program. For ethynyl formate, the calculation for 25 °C at the MP2/6-311++G(df,pd) level predicts that the Z isomer is more stable by 1.23 kcal/mol. For S-ethynyl thioformate, calculations at the MP2/6-311++G(2d,2p) level predict that the E isomer is favored by 0.71 kcal/mol at 25 °C. The E isomers of N-ethynylformamide and N-ethynylacetamide were found at all levels to be more stable than the Z isomers at 25 °C. For cyano formate and cyano acetate, calculations at the MP2/6-311++G(df,pd) level predict that the Z isomers are more stable at 25 °C by 1.50 and 2.72 kcal/mol, respectively. At this level and temperature, the Z isomers of 2c, 2d, 3c, and 3d are predicted to have free energies of 0.46, −0.07, 1.22, and 2.28 kcal/mol, respectively, relative to the E conformations. Z to E free-energy barriers at 25 °C of 8.63, 10.64, 17.63, 7.39, and 14.03 kcal/mol were calculated for 1a, 2a, 3a, 1c, and 3c at the HF/6-311G(d,p) level, and at the HF/6-311+G(d,p) level, the free-energy barrier for 2c was 7.08 kcal/mol.