Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.     

Ab initio study of the conformational equilibrium of ethylene glycol

Summary The conformational equilibrium of ethylene glycol (CH₂OHCH₂OH) has been examined by performing geometry optimizations at the 6-31G*, MP2/6-31G* and 6-31G** levels. Final energies have been calculated at the MP3 level with the optimized geometries. The two most stable conformers are atGg andgGg but it is verified that the inclusion of electronic correlations reduces their energy difference of 0.6 kcal/mol at the HF level to less than 0.2 kcal/mol. The possible coexistence of the two most stable conformers is in agreement with some previous studies of Frei et al. For thetXg conformer a detailed analysis of the intramolecular potential as a function of rotation around the C-C bond is also reported.     

An Ab Initio Molecular Orbital Theory and Density Functional Theory (DFT) Study of Conformers and Rotamers of 4-Substituted (Methyl, Hydroxymethyl, Methanoyl, Ethanoyl, Cyano, Fluoro, Chloro, Bromo, Acetoxy) Tetrahydro-2H-thiopyran-1,1-dioxides

The structures and energies of axial and equatorial conformers and rotamers of 4-substituted tetrahydro-2H-thiopyran-1,1-dioxides (tetrahydrothiopyran-1,1-dioxides, thiacyclohexane-1,1-dioxides, thiane-1,1-dioxides, and 1,1-dioxothianes; CH₃, CH₂OH, CHO, COCH₃, CN, F, Cl, Br, and OCOCH₃) were calculated using the hybrid density functionals B3LYP, B3P86, and B3PW91, as well as MP2 and the 6-31G(d), 6-31G(2d), 6-31G(3d), 6-31G(d,p), and 6-31+G(d) basis sets. MP2/6-31+G(d)/ /HF/6-31+G(d) [–G° = 1.73 kcal/mol], B3P86/6-31G(d) [–G° = 1.75 kcal/mol], and B3PW91/6-31G(d) [–G° = 1.85 kcal/mol] gave conformational free energy (G°) values at 180 K for 4-methyltetrahydro-2H-thiopyran-1,1-dioxide which were similar to the reported experimental values for methylcyclohexane (–G° = 1.80 kcal/mol), 4-methyltetrahydro-2H-thiopyran (–G° = 1.80 kcal/mol), and other 4-methyl-substituted heterocycles. All levels of theory showed that the conformational preferences of the 4-methanoyl (4-formyl), 4-ethanoyl (4-acetyl), and 4-cyano substituents were small. The HF calculations gave conformational free energy (G°) values for 4-chlorotetrahydro-2H-thiopyran-1,1dioxide which were closer to the experimental value than the MP2 and density functional methods. The best agreement with available experimental data for 4-bromotetrahydro-2H-thiopyran-1,1-dioxide was obtained from the HF/6-31G(2d), HF/6-31G(3d), and B3LYP/6-31G(2d) calculations, and, for 4-acetoxytetrahydro-2H-thiopyran-1,1-dioxide, from the HF/6–31G(3d) calculations. The conformational free energies (G°) and relative energies (E) of the conformers and rotamers have been compared with the correspondingly substituted cyclohexanes and tetrahydro-2H-thiopyrans and are discussed in terms of dipole–dipole (electrostatic) interactions and repulsive nonbonded interactions (steric) in the most stable axial and equatorial conformers. The axial S=O bond lengths are shorter than the equatorial S=O bond lengths and the C2–C3 bond lengths in the substituents with carbon-bonded to the ring are shorter than the C3–C4 and C4–C-5 bond lengths. In contrast, the C2–C3 bond lengths in the 4-halogen and 4-acetoxy substituents are longer than the C3–C4 and C4–C-5 bond lengths.     

Estimation of the relative stability of chair forms of hexahydropyrimidine

Nonempirical approximations HF/6-31G(d), MP2/6-31G(d)//HF/6-31G(d), PBE/3z, and RI-MP2/λ2 were applied to the calculation of the relative stability and parameters of the conformational equilibrium of three chair forms of the unsubstituted hexahydropyrimidine with different orientation of the N-H bonds. It was established that in the framework of the used methods the most stable are the conformers of diaxial and axial-equatorial chairs. The occupancy of the diequatorial chair form at the room temperature does not exceed 1.7%.     

Comparative conformational analysis of 1,3-dioxane and 1,3-dithiane

It was shown by ab initio quantum-chemical approximations HF/6-31G(d) and MP2/6-31G(d)//HF/6-31G(d) that the conformational isomerization of 1,3-dioxane and 1,3-dithiane proceeded along common routes. The potential energy surface of both compounds contains six minima including the chair invertomers and enantiomeric flexible forms. They are separated by several potential barriers. It was established by molecular dynamics method that the flexible conformers at heating and keeping at 295–300 K transformed into each other and in the chair conformer.     

A QM/MM-based computational investigation on the catalytic mechanism of saccharopine reductase

Saccharopine reductase from Magnaporthe grisea, an NADPH-containing enzyme in the α-aminoadipate pathway, catalyses the formation of saccharopine, a precursor to L-lysine, from the substrates glutamate and α-aminoadipate-δ-semialdehyde. Its catalytic mechanism has been investigated using quantum mechanics/molecular mechanics (QM/MM) ONIOM-based approaches. In particular, the overall catalytic pathway has been elucidated and the effects of electron correlation and the anisotropic polar protein environment have been examined via the use of the ONIOM(HF/6-31G(d):AMBER94) and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) methods within the mechanical embedding formulism and ONIOM(MP2/6-31G(d)//HF/6-31G(d):AMBER94) and ONIOM(MP2/6-311G(d,p)//HF/6-31G(d):AMBER94) within the electronic embedding formulism. The results of the present study suggest that saccharopine reductase utilises a substrate-assisted catalytic pathway in which acid/base groups within the cosubstrates themselves facilitate the mechanistically required proton transfers. Thus, the enzyme appears to act most likely by binding the three required reactant molecules glutamate, α-aminoadipate-δ-semialdehyde and NADPH in a manner and polar environment conducive to reaction.     

Theoretical study of the gas-phase thermolysis reaction of alkyl (ethyl,isopropyl, and tert-butyl) N,N-dimethylcarbamates and N,N-diethylcarbamates

Theoretical studies on the thermolysis in the gas phase of alkyl N,N-dialkylcarbamates were carried out using ab initio theoretical methods, at the MP2/6-31G(d), MP2/6-31++G(d,p) and MP2/6-311++G(2d,p)//MP2/6-31G(d) levels. The reactions have two steps: the first one corresponds to the formation of an alkene and a neutral dialkylcarbamic acid intermediate via a six-membered cyclic transition state; the second one is the decarboxylation of this intermediate via a four-membered cyclic transition state, leading to carbon dioxide and the corresponding dialkylamine. The progress of the reactions was followed by means of the Wiberg bond indices. The results indicate that the transition states have character intermediate between reactants and products, and the calculated synchronicities show that the reactions are slightly asynchronous. The bond-breaking processes are more advanced than the bond-forming ones, indicating a bond deficiency in the transition states. The rate constants calculated for all the reactions agree very well with the available experimental data.From the Proceedings of the 28th Congreso de Químicos Teóricos de Expresión Latina (QUITEL 2002)     

Conformational analysis of 2-methyl-5-nitro-1,3,2-dioxaborinane

Quantum-chemical methods HF/6-31G(d), HF/6-31+G(d), MP2/6-31G(d)//HF/6-31G(d), and MP2/6-31+G(d)//HF/6-31+G(d) were used to investigate the conformational isomerization of 2-methyl-5-nitro-1,3,2-dioxaborinane. It has been shown that a potential energy surface of this compound includes two minima: an axial form of semi-chair and equatorial sofa together with a transition state belonging to the conformation of 2,5-twist-form. A comparison between experimental NMR ¹H and theoretical vicinal coupling constants was used to determine the quantitative conformational composition of cyclic boric acid ester and a value of ΔG ⁰ for nitro group at the ring carbon atom C⁵ in CCl₄ and C₆D₅NO₂ solutions.     

Oxo-hydroxy tautomerism of 5-fluorouracil: water-assisted proton transfer

Post-Hartree-Fock ab initio quantum chemical calculations were performed for 5-fluorouracil in the gas phase and in a three-water cluster. Full geometry optimizations of the 5-fluorouracil-water complexes were carried out at the MP2/6-31+G(d,p) level of theory. MP4/6-31+G(d,p)//MP2/6-31+G(d,p) and MP4/6-31++G(d,p)//MP2/6-31+G(d,p) single-point calculations were performed to obtain more accurate energies. In water solution, 5-fluorouracil exists mainly in the 2,4-dioxo form (A). We propose that the populations of the 2-hydroxy-4-oxo (B) and 4-hydroxy-2-oxo (D) tautomers are 1 x 10(-4)% and 3.9 x 10(-8)%, respectively, on the basis of the relative stabilities of the tautomers calculated at the MP4/6-31++G(d,p)//MP2/6-31+G(d,p) level of theory. A profound difference between isolated and hydrated 5-fluorouracil is noted for the height of the tautomerization barrier. In the absence of water, the process of proton transfer is very slow. The addition of water molecules decreases the barrier by 2.3 times, making the process much faster. The minimum energy path (MP2/6-31+G(d,p)) for water-assisted proton transfer in trihydrated 5-fluorouracil was followed. CNDO/S-CI calculations predict singlet pi-pi(*) electron transitions at 312 nm for B and at 318 nm for D. The fluorescence spectrum of 5-fluorouracil in water confirms the presence of the hydroxy tautomer.     

Theoretical study of the addition and abstraction reactions of hydroxyl radical with uracil

The addition as well as abstraction reactions of hydroxyl radical (OH) with the nucleic acid base, uracil (U), in the gas phase has been explored at the B3LYP/6-31+G(d,p) level of density functional theory (DFT). The energy barrier of the OH addition to both the C5 and C6 positions of the uracil is less than 1 kcal/mol while the hydrogen abstractions (H-abstractions) from either the N1 or the N3 positions are ∼9.5 kcal/mol. Further the energetics of these reactions are assessed by applying the effect of aqueous medium through the polarizable continuum model (PCM). Both the gas and the solution phase data established that the thermodynamic and kinetic factors are more favorable for the OH addition to either C5 or C6 positions of the uracil than the H-abstraction reactions. Moreover, calculations at the MPW1K/6-31+G(d,p), CCSD(T)/6-31+G(d,p)//B3LYP/6-31+G(d,p) and CCSD(T)/6-31+G(d,p)//MPW1K/6-31+G(d,p) levels of theoretical methods qualitatively supported the B3LYP/6-31+G(d,p) results.     

Modern valence-bond description of chemical reaction mechanisms: the 1,3-dipolar addition of diazomethane to ethene

The electronic mechanism for the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane (CH2N2) to ethene (C2H4) is described through spin-coupled (SC) calculations at a sequence of geometries along the intrinsic reaction coordinate obtained at the MP2/6-31G(d) level of theory. It is shown that the bonding rearrangements occurring during the course of this reaction follow a heterolytic pattern, characterized by the movement of three well-identifiable orbital pairs, which are initially responsible for the pi bond in ethene and the C-N pi bond and one of the N-N pi bonds in diazomethane and are retained throughout the entire reaction path from reactants to product. Taken together with our previous SC study of the electronic mechanism of the 1,3-dipolar cycloaddition of fulminic acid (HCNO) to ethyne (C2H2) (Theor. Chim. Acc. 1998, 100, 222), the results of the present work suggest strongly that most gas-phase concerted 1,3-dipolar cycloaddition reactions can be expected to follow a heterolytic mechanism of this type, which does not involve an aromatic transition state. The more conventional aspects of the gas-phase concerted 1,3-dipolar cycloaddition of diazomethane to ethene, including optimized transition structure geometry, electronic activation energy, activation barrier corrected for zero-point energies, standard enthalpy, entropy and Gibbs free energy of activation, have been calculated at the HF/6-31G(d), B3LYP/6-31G(d), MP2/6-31G(d), MP2/6-31G(d,p), QCISD/6-31G(d) and CCD/6-31G(d) levels of theory. We also report the CCD/6-311++G(2d, 2p)//CCD/6-31G(d), MP4(SDTQ)/6-311++G(2d,2p)//CCD/6-31G(d) and CCSD(T)/6-311++G(2d, 2p)//CCD/6-31G(d) electronic activation energies.     

A theoretical study of 1-amino-3-butene and 3-butene-1-thiol

Conformational analysis of 1-amino-3-butene and 3-butene-1-thiol was carried out using the 4-21G basis set. The conformers obtained were subjected to 6-31G^* single-point analysis for the calculation of energies, charge distributions, and dipole moments. The geometries and stabilities obtained are in good agreement with available experimental data. The results are interpreted in terms of intramolecular hydrogen bonding and anomeric interactions: Some of the most stable conformers of both molecules have intramolecular hydrogen bonds between the hydrogens of the amino or thiol groups and the electrons of the double bond. The 4-21G geometries were refined to obtain rotational constants closer to the experimental values.     

Intramolecular hydrogen bond in the push–pull CF3-aminoenones: DFT and FTIR study,NBO analysis

Postulated conformers of trifluoromethylated β-aminoenones stabilized by intramolecular NH?O and N?HO bonds were studied by IR and NMR spectroscopy and evaluated with quantum chemical calculations (B3LYP/6-311+G(d,p), MP2/6-311+G(d,p)//B3LYP/6-311+G(d,p) and MP2/6-31G(d,p)) and NBO analysis. The influence of the nature of EWG, substituents at the nitrogen atom and double bond, and of orbital interactions of heteroatoms and double bonds in these structures on the proton affinity of basic and acid centers, strength of hydrogen bonds, and the energy of tautomeric transfers is discussed. The theoretical results agree satisfactorily with the experimental observations.     

Steric and Electronic Structure of Selenoanisole: A Quantum-Chemical Study

The potential functions of internal rotation around the C_sp2-Se bond in selenoanisole were ob- tained by quantum-chemical calculations in the approximations HF/3-21G(d), HF/6-31G(d), MP2(f)/6-31G(d), and B3LYP/6-31G(d). The calculations were performed in the range of variation of the torsion angle (between the planes of the benzene ring and C_sp2-Se-C_sp3 bonds) from 0° to 90° with 15° step. The energy minimum is in the region of the orthogonal conformation ( 90°), and the energy maximum, in the region of the planar form ( 0°). The rotation barriers (kJ mol^{- 1}) are as follows: HF/3-21G(d), 9.20; HF/6-31G(d), 13.13; MP2(f)/6-31G(d), 10.25; and B3LYP/6-31G(d), 6.41. The geometric parameters, Koopmans ionization potentials, and dipole moments are given. The energies, degrees of hybridization, populations of the lone electron pairs of Se, energies of their interaction with the antibonding * orbitals of the benzene ring, and electron density distributions were determined in terms of the natural bond orbital approach.     

Molecular structure of 6-(N,N-dimethylamino)fulvene from a joint gas-phase electron diffraction and quantum-chemical study

The molecular structure of 6-(N,N-dimethylamino)fulvene was studied by gas-phase electron diffraction and quantum-chemical methods (HF/6-31G(d), MP2/3Z, MP2/4Z, density functional theory with the B3LYP/6-31G(d) and PBE/3Z functionals). Pronounced flattening of the nitrogen atom and equalization of the intracyclic C—C bonds were found to be a consequence of the electron delocalization in the molecule.     

Quantum-chemical investigation of conformational characteristics of the complex of 2-methyl-1,3,2-dioxaborinane with dimethyl ether

With the help of the quantum-chemical approximation MP2/6-31G(d)//HF/6-31G(d) conformational characteristics were investigated of donor-acceptor associate of 2-methyl-1,3,2-dioxaborinane with dimethyl ether, 1:1. The character of the conformational isomerization of the adduct differs from that established for the proper cyclic boric ester and involves an interconversion between two conformers nondegenerate by the energy. The stability of the associate under study is comparable with those of the molecular complexes of 2-methyl-1,3,2-dioxaborinane with water and methanol.     

Conformational Study of Naringenin in the Isolated and Solvated States by Semiempirical and Ab Initio Methods

Naringenin is a natural widespread flavanone occurring in different foodstuffs that presents several important biological activities. Although its properties are well documented, its mechanisms of action are still controversial. The present article reports a conformational analysis of naringenin, using the semiempirical AM1 and ab initio methods, at the Hartree–Fock level of theory. The 3-21G, 3-21G*, 6-31G, and 6-31G** basis sets were used. The electron correlation effects were included through the Møller–Plesset second-order perturbation theory. The solvation of naringenin has been investigated through the standard SCRF, the supermolecule (SM), and the combined SM/SCRF models. The results have shown that there are two degenerate forms of naringenin, differing mainly by the orientation of a hydroxyl group (C4—OH). The energy barrier for the interconversion between them is ca. 6 kcal.mol^–1, suggesting some conjugation between the -system of the aromatic B ring and the hydroxyl group (C4—OH).     

On the dimerization process of nitroso compounds

Summary Many organic C-nitroso compounds R-NO form stable dimers with a covalent NN bond. To gain insight into the dimerization reaction 2 R-NO (R-NO)₂ a theoretical study of the dimerization to atrans-form was performed using HNO as a model compound. Complete geometry optimizations were carried out at the HF, MP2 and QCISD levels using a 6–31G* basis. In the stationary points energies were calculated at the MP4(SDTQ) and QCISD(T) levels. For the equilibrium structure of the monomer and dimers stable RHF solutions were found, whereas for the TS UHF and UMPn calculations were applied. Extensive spin contamination was found in the UHF wavefunction, and projections up tos+4 were invoked. Relative energies were corrected for differences in ZPE. Calculations were made (a) for the least-motion path (C _2h symmetry) and (b) for a path with complete relaxation of all internal coordinates. Along the latter path a TS having virtuallyC _i symmetry was found. Along path (a) an activation energy of around 150 kcal/mol was predicted, in conformity with a symmetry forbidden reaction. On the relaxed path (b) the barrier to dimerization was estimated to be 10.7 kcal/mol at the MP4(SDTQ)//MP2 level, and 10.9 kcal/mol at the QCISD(T)//QCISD level. Unscaled ZPE corrections, calculated at the SCF level, changed these values to 12.7 and 12.9 kcal/mol, respectively. The reaction energy for the dimerization process is predicted to be – 17.2 kcal/mol at the MP4(SDTQ)//MP2 level corrected for ZPE. Calculations at the G1 level gave a corresponding value of – 16.4 kcal/mol. The equilibrium constant for the association to thetrans dimer is estimated to beK _p =259 atm, indicating that the dimer should be an observable species in the gas phase.     

Absence of p,π-conjugation in the Cl2C=CHOMe molecule from data ofab initio calculations

ab initio Calculation of the Cl₂C=CHOMe molecule has been carried out with geometrical optimization in the split valence basis set in the RHF/6-31G^*//RHF/6-31G^* approximation. An analysis of the populations of the p-orbitals of the Cl atoms implies the absence of p,-conjugation between the unshared electron pairs of these atoms and the -bond. The asymmetry parameters of the electric field gradient at the³⁵Cl nuclei, calculated from the 3p-orbital populations of the Cl atoms, practically coincide with the experimental values.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2369–2372, December, 1995.The authors are grateful to G. S. Beloglazov for assistance in accomplishingab initio calculations.