Facial selectivity in epoxidation of 2-cyclohexen-1-ol with peroxy acids. A computational DFT study

We addressed the mechanism of epoxidation of 2-cyclohexen-1-ol by locating all the transition structures (TSs) for the reaction of peroxyformic acid (PFA) with both pseudoequatorial and pseudoaxial cyclohexenol conformers (five TSs for each conformer) and, for purpose of comparison, also those for the PFA epoxidation of cyclohexene. Geometry optimizations were performed at the B3LYP/6-31G level, energies refined with single point B3LYP/6-311+G// B3LYP/6-31G calculations and solvent effects introduced with the CPCM method. Our results can be summarized as follows: (i) all TSs exhibit a spiro-like structure, that is, the dihedral angle between the peroxy acid plane and the forming oxirane plane is closer to 90 degrees than to 0 degrees (or 180 degrees ); (ii) there is a stabilizing hydrogen bonding interaction in syn TSs that, however, is partly counteracted by unfavorable entropic effects; (iii) syn,exo TSs with hydrogen bonding at the PFA peroxy oxygens are definitely more stable than syn,endo TSs hydrogen bonded at the PFA carbonyl oxygen; (iv) facial selectivity of epoxidation of both cyclohexenol conformers is mostly the result of competition between only two TSs, namely, an anti,exo TS and its syn,exo counterpart. The latter TS is more stable than the former one, as stabilization by hydrogen bonding overrides the unfavorable entropic and solvent effects; (v) calculations correctly predict both the experimental dominance of attack leading to syn epoxide for both cyclohexenol conformers and the higher syn selectivity observed for the pseudoequatorial as compared to the pseudoaxial derivative. Moreover, also the experimental relative and absolute epoxidation rates for cyclohexene and cyclohexenol as well as for pseudoaxial and pseudoequatorial cyclohexenol derivatives are fairly well reproduced by computational data.     

Bisaziridine tetracarboxylates as building blocks in the stereoselective synthesis of C60-fullerene diads and dumbbell-like bis-C60-fullerene triads

The synthesis of alkoxycarbonyl-substituted bisaziridines with the two aziridine units connected by conjugated p-phenylene, partly conjugated 1,1'-biphenyl-4,4'-diyl, and nonconjugated 4,4'-methylenediphenyl linkers was developed. The reaction of azomethine ylides derived from the bisaziridines with fullerene C(60) was optimized and used for the stereoselective preparation of both the monoadducts (C(60)-linker-aziridine dicarboxylate), and the dumbbell bisadducts (C(60)-linker-C(60)). The reasons for the observed selectivity of the azomethine ylide formation and cycloaddition were theoretically studied at the DFT B3LYP/6-31G(d) level or at the ONIOM B3LYP/6-31G(d):B3LYP/STO-3G level for fullerene-containing molecules.     

Cyclic vinyl p-tolyl sulfilimines as chiral dienophiles: theoretical study on the stereoselectivity, Lewis acid catalysis, and solvent effects in their Diels-Alder reactions

The theoretical study reported in the present work deals with chiral cyclic vinyl sulfilimines and their reactivity as dienophiles in [4 + 2] cycloaddition reactions, using B3LYP/6-31G(d)//AM1 and B3LYP/6-31G(d)//B3LYP/6-31G(d) model chemistries. Consideration of Lewis acid catalysis, illustrated by BF(3), decreases the activation energies of the cycloaddition process while the charge transfer from the diene to the sulfilimine is augmented. The [4 + 2] cycloaddition reactions of sulfilimines with both furan and cyclopentadiene occur in the gas phase with endo stereoselectivity, which is more pronounced with the latter diene. Endo-exo energy differences in the gas phase with the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31G(d)//AM1 model chemistries are almost the same. Solvent effects are responsible for the inversion of the stereoselectivity in the reactions of sulfilimines with furan because of the great difference in the dipole moments in endo and exo approaches.     

Endo-exo and facial stereoselectivity in the Diels-Alder reactions of 3-substituted cyclopropenes with butadiene

A computational examination of the four modes of addition in the Diels-Alder reactions of 3-substituted cyclopropene derivatives (substituents: BH(2), CH(3), SiH(3), NH(2), PH(2), OH, SH, F, and Cl) with butadiene have been carried out at the B3LYP/6-31++G(d)//HF/6-31++G(d) level. The degree of stabilization of these derivatives at the ground state correlates with the electronegativity of the substituent. This attenuation of reactivity and differences in steric interactions are the only factors needed to explain both the high facial selectivity and the differences in the endo-exo selectivity seen in these reactions. Furthermore, evidence is presented that indicates that stabilization by an interaction involving the syn C-3 hydrogen of cyclopropene and butadiene is small or irrelevant in controlling the endo-exo selectivity of the Diels-Alder reaction.     

Septanose carbohydrates: synthesis and conformational studies of methyl alpha-D-glycero-D-idoseptanoside and methyl beta-D-glycero-D-guloseptanoside

We report the synthesis of methyl alpha-D-glycero-D-idoseptanoside (1) and methyl beta-D-glycero-D-guloseptanoside (2) and the characterization of their preferred solution conformations by computational chemistry and (1)H NMR (3)J(H,H) coupling constant analysis. Central to the synthetic approach was the epoxidation of glucose-derived oxepine 3 using DMDO. Nucleophilic attack on the resulting 1,2-anhydroseptanose using NaOCH(3) in CH(3)OH followed by deprotection provided the 1,2-trans diastereomers 1 and 2. The computational approach for determining the preferred low energy septanose conformations began with a pseudo Monte Carlo search for each isomer using minimization with the AMBER force field. Single-point energy calculations (HF/6-31G *and B3LYP/6-31+G**) as well as full geometry optimizations in a model for aqueous solvent were then conducted using the conformers within 5 kcal/mol of the AMBER global minimum. Calculated (3)J(H,H) values, based on a Boltzmann distribution of the computed low energy conformers, were compared to experimental (3)J(H,H) values from (1)H NMR coupling constant analyses. The correlation between calculated and observed values suggest that septanose carbohydrates are not so flexible and should generally prefer one twist-chair (TC) conformation.     

A theoretical rationalization of the asymmetric induction in sulfinyl-directed [5C + 2C] intramolecular cycloadditions

A computational DFT examination (B3LYP/6-31G and B3LYP/6-311+G//B3LYP/6-31G) of the thermal [5C + 2C] cycloaddditions of 6-acetoxy-3-pyranones and 3-silyloxy-4-pyrones with tethered alkenyl sulfoxides confirms that the high level of diastereofacial selectivity obtained is ultimately due to the preference of the alkenesulfinyl group to adopt a well-defined conformation in the transition state of the reaction. This conformation, which is different from that found in the ground state, is most probably dictated by dipolar interaction effects between the sulfoxide and the oxidopyrylium ylide intermediate.     

Development of predictive models of π-facial selectivity; a critical study of nucleophilic addition to sterically unbiased ketones

Quantum chemical calculations at B3LYP/6-31G* and semiempirical levels have been performed on a series of sterically unbiased ketones, where facial differentiation during nucleophilic additions is electronically induced through distal functional groups. The face selectivity data for fifty-four substrates representing nine different skeleta were computed and compared with the available experimental data on thirty-eight of them. The predictive abilities of various computational methods such as, charge model, hydride model, LiH transition state model, Cieplak hyperconjugation effect estimated by NBO analysis and the cation complexation model have been evaluated. A comparison of the computed and experimental face-selectivity data indicates that the hydride model and the LiH transition state model at the semiempirical levels are the best choices to predict diastereoselectivity. Unexpectedly, the performance of charge, hydride and LiH transition state models are inferior at the B3LYP level compared to the semiempirical methods in predicting the facial selectivities. On the other hand, the Cieplak type hyperconjugation evaluated using the NBO analysis, and cation complexation model are less reliable despite the fact that these two involve higher (B3LYP/6-31G*) level calculations. The inadequate performance of the charge model, NBO and the cation complexation models were traced to their emphasis on only one or two factors which are responsible for stereodifferentiation and undermining of the other subtle aspects involving a combination of orbital and electrostatic effects. On the other hand, the hydride and LiH transition state models, at semiempirical levels, provide reliable results to model the face-selectivities.     

Quantum-Chemical Study on Positional Selectivity in the Trimethylsilylation and Sulfonation of Pyrrole and <Emphasis Type="Italic">N</Emphasis>-Alkylpyrroles

Positional selectivity (α:β ratio) of electrophilic substitution in pyrrole, N-methylpyrrole, and N-tert-butylpyrrole was analyzed by ab initio [RHF/6-31G(d), MP2/6-31G(d)//RHF/6-31G(d)] and DFT [B3LYP/6-31G(d)] calculations of some σ-complexes derived from the substrates. The results of calculations with the use as model electrophilic species of trimethylsilyl cation [MP2/6-31G(d)//RHF/6-31G(d) and B3LYP/6-31G(d)] and SO₃ molecule [B3LYP/6-31G(d)] instead of proton are fairly consistent with the experimental data, according to which trimethylsilylation of pyrrole and its N-substituted derivatives with trimethylsilyl trifluoromethanesulfonate, as well as sulfonation with pyridine-sulfur trioxide complex, gives the corresponding β-substituted products.     

丁二酰亚胺的结构、振动频率和热力学性质计算

用密度泛函理论（ＤＦＴ）和从头计算（ａｂ ｉｎｉｔｉｏ）方法,在Ｂ３ＬＹＰ／６－３１Ｇ、Ｂ３ＬＹＰ／６－３１１Ｇ和ＭＰ２／６－３１１Ｇ水平上,全优化计算了了二酰亚胺的分子几何构型和电子结构．进行了简正振动频率分析并用校正后的频率计算了２００－６００Ｋ温度范围的标准热力学函数,对计算结果进行了比较和讨论．     

Comparison of Some Computational Methods for Geometry Optimisation of Phosphorus Acid Derivatives

Several economical methods for geometry optimization, that should be applicable to larger molecules, have been evaluated for 19 phosphorus acid derivatives. MP2/cc-pVDZ geometry optimizations are used as reference points and the geometries obtained from the other methods are evaluated with respect to deviations in bond lengths and angles, from the reference geometries. The geometry optimization methods are also compared to the much used B3LYP/6-31G(d) method. Single point energies obtained by subsequent EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) calculations on the respective equilibrium geometries are also reported relative to the energies obtained from the reference geometries. The geometries from HF/MIDI! optimizations were closer to those of the references than the geometries of the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! optimizations. The EDF1/6-31+G(d) or B3LYP/6-31+G(d,p) single point energies obtained from the HF/3-21G(d), HF/6-31G(d), and B3LYP/MIDI! geometries gave a mean absolute deviation (MAD) from that of the reference geometries of 1.4-3.9 kcal mol m 1 . The HF/MIDI! geometries, however, gave EDF1/6-31+G(d) and B3LYP/6-31+G(d,p) energies with a MAD of only about 0.5 and 0.55 kcal mol m 1 respectively from the energies obtained with the reference geometries. Thus, use of HF/MIDI! for geometry optimization of phosphorus acids is a method that gives geometries of near-MP2 quality, resulting in a fair accuracy of energies in subsequent single point calculations, at a much lower computational cost other methods that give similar accuracies.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

2,3,7,8-四氯苯并二英分子结构与热力学性质的理论研究

用密度泛函理论(DFT)和从头算(ab initio)方法,在B3LYP/6-31G、 B3LYP/6-31G*、 B3LYP/6-311G*和MP2/6-31G*水平上全优化计算了2,3,7,8-四氯苯并二英(2,3,7,8-TCDD)的几何构型、电子结构和振动频率,并用校正后的频率计算了298～1500 K的标准热力学函数,同时用半经验的PM3 SCF-MO进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

电场作用下分子导线的理论研究

摘要利用从头计算法分别在HF/6-31G, HF/6-31G*, HF/6-31G**, HF/6-31+G, HF/6-31++G, HF/6-31+G*, HF/6-31+G**, HF/D95+*, B3LYP/6-31G*和B3LYP/6-31+G*水平上计算了5个单体的聚乙炔分子导线, 从几何构型、 SCF能量和分子轨道能级三个方面讨论了外电场对分子导线的影响, 给出了聚乙炔分子导线性质与外电场变化的定量关系.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Computational modeling of a stereoselective epoxidation: reaction of carene with peroxyformic acid

Electronic structure theory was used to model the epoxidation of 3-carene by peroxyformic acid. Reactants, products, and transition states were optimized at the B3LYP/6-31G* level of theory, followed by B3LYP/6-311+G** and MP2/6-311+G** single point calculations. The reaction pathway yielding the trans-epoxide product was found to have a significantly lower reaction barrier (7.8 kcal/mol) than that leading to the cis-epoxide product (9.4 kcal/mol), in agreement with expectations. Magnetic shieldings of the two isomeric carene epoxides were also calculated, using the GIAO method, and compared to experimental (1)H and (13)C NMR spectra. Although the calculated carbon spectra proved inconclusive, the proton shieldings calculated for the trans-epoxide correlated much more closely to the experimental values for the major epoxidation product than did the shieldings calculated for the cis-epoxide, serving to verify the identity of the major product.     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

Regioselective enzymatic acylations of polyhydroxylated eudesmanes: semisynthesis, theoretical calculations, and biotransformation of cyclic sulfites

Different lipase enzymes have been tested in order to perform regioselective acetylations on the eudesmane tetrol from vulgarin. High yields (95%) of 1,12-diacetoxy derivative (4) were achieved in 1 h with Candida antarctica lipase (CAL). However, only the 12-acetyl derivative (6) was obtained in similar yield with Mucor miehei (MML) or Candida cylindracea (CCL) lipases. The enzymatic protection at C-1 and C-12 has been used to form eudesmane cyclic-sulfites between C-6 and C-4 atoms. The R/S-sulfur configuration has been assigned by means of the experimental and theoretical (13)C and (1)H NMR chemical shifts. The theoretical shifts were calculated using the GIAO method, with a MM+ geometry optimization followed by a single-point calculation at the B3LYP/6-31G(*) level (B3LYP/6-31G(*)//MM+). Moreover, B3LYP/6-31G(*) geometry optimizations were carried out to test the B3LYP/6-31G(*)//MM+ results, for the deacetylated sulfites (12 and 15). In addition to the delta(C) and delta(H) shifts, the (3)J(HH) coupling constants were also calculated and compared with the experimental values when available. Finally, different reactivities have been checked in both sulfites by biotransformation with Rhizopus nigricans. While the R-sulfite gave 2 alpha- and 11 beta-hydroxylated metabolites, the S-sulfite yielded only regioselective deacetylations. Furthermore, both sulfites showed different reactivities in redox processes.     

Non-aqueous solvation of n-octanol and ethanol: spectroscopic and computational studies

Raman spectroscopy was used to examine the interactions of the free O-H bonds in n-octanol and ethanol with the organic solvents carbon tetrachloride (CCl(4)), cyclohexane, and benzene. These spectra reveal that the solvents CCl(4) and cyclohexane have a small effect on the free O-H peak of alcohols, whereas benzene as a solvent significantly red-shifts the free O-H band. Calculated spectra were generated via MP2/6-31G* calculations and the B3LYP/6-31+G**//MP2/6-31G*-derived Boltzmann populations of each ethanol complex and are consistent with the experimental results. Additional spectra were calculated using Boltzmann populations derived from single-point energies at the polarizable continuum model (PCM) level with the B3LYP/6-31+G** level of theory to take overall solvent effects into account, and these simulated spectra are also largely consistent with the experimental results. Analysis of the computational results reveals a lengthening of the O-H bond from the O-H interaction with the delocalized electronic structure of benzene as well as a bimodal distribution of the free O-H peak of the alcohol/benzene mixtures due to two distinctly different types of alcohol/benzene complexes.     

An efficient method for calculating atomic charges of peptides and proteins from electronic populations

An efficient method is presented to calculate atomic charges of peptides and proteins derived from Mulliken electronic populations for terminally blocked amino acids (Ac-X-NHMe, X = any neutral or charged residue) calculated at the B3LYP/6-31G(d,p)//HF/6-31G(d,p) level of theory. This electronic population-derived atomic charges (EPAC) method is based on the geometry-dependent atomic charge (GDAC) method proposed by Cho et al. (J. Phys. Chem. B 2001, 105, 3624), in which atomic charges are calculated by using the partial equalization of atomic electronegativities with electronegativity parameters and damping factors given by interatomic distances between covalently bonded atoms in a molecule. The overall mean absolute difference (mad) and root-mean-square deviation (rmsd) between dipole moments micro(EPAC) and micro(B3LYP), obtained from EPAC charges and from B3LYP/6-31G(d,p) level calculations, respectively, for Ac-X-NHMe are estimated to be 0.38 and 0.59 D, respectively, for 26 representative conformations in the training set, and 0.54 and 0.79 D, respectively, for 172 representative conformations not used for parametrization. For Ac-(Ala)(n)-NHMe (n = 2-6), the EPAC method reasonably predicts the increase of the dipole moment with increase of the chain length, although the deviations from the micro(B3LYP) values are somewhat larger. For Ac-Ala-NHMe and Ac-(Ala)(3)-NHMe, the EPAC charge for a specific type of atom does not depend on its position in the sequence or on the length of the sequence. In addition, charge neutrality holds for any Ala residue of these two peptides. Thus, these results suggest that the EPAC charges derived from B3LYP/6-31G(d,p) Mulliken populations can be used reliably for conformational analysis of peptides and proteins.     

Si~2Cl~6分子的内旋转能垒和振动基频的理论研究

用从头算方法HF/6-31G^*^*和密度函方法B3LYP/6-31G^*^*,对Si~2Cl~6分子的平衡几何构型进行优化,优化的结果与实验结果吻合得较好.并用上述两种不同的方法计算Si~2Cl~6分子的内旋转能垒,结果分别为8.786和6.694kJ/mol,其中DFT方法的计算结果与实验结果4.18kJ/mol吻合得较好.对Si~2Cl~6分子的振动基频进行计算.用HF/6-31G^*^*SQM力场所计算的频率理论值与实验值的平均误差为7.3cm^-^1,用B3LYP/6-31G^*^*未标度的力场所计算的频率理论值与实验值的平均误差为6.0cm^-^1.该密度泛函方法(B3LYP/~6-31G^*^*)的理论计算值比用HF/6-31G^*^*标度后的SQM力场计算的频率与实验值(除Si--Si键扭转振动基频之外的11条振动基频)吻合得更好.并给出了Si--Si键扭转振动基频的预测值。