Preparation of heterocycle-masked beta-enamino acids

The preparation of masked N-substituted beta-enamino acid derivatives 3a-n by reaction of 2-alkyl-oxa(thia)zolines 2a,b with imidoylbenzotriazoles 1a-g in the presence of LDA was described.     

A theoretical study on the structure and properties of phenothiazine derivatives and their radical cations

Semiempirical CNDO, AM1, PM3 and ab initio HF/STO-3G, HF/3-21G(d), and HF/6-31(d) methods were employed in the geometry optimization of the phenothiazine and the corresponding radical cation. The results obtained from the PM3 performances were as good as those from the ab initio calculations in the structure optimization of both phenothiazine and phenothiazine radical cation. The PM3 method was used to optimize the structures of a series of N-substituted phenothiazine derivatives and their radical cations. The PM3-optimized results were then analyzed with the ab initio calculation at the 6-311G(d,p) level, which yielded the total energy, frontier molecular orbitals, dipole moments, and charge and spin density distributions of the phenothiazine derivatives and their radical cations.     

New strategy for the stereoselective synthesis of fluorinated beta-amino acids

Racemic and chiral nonracemic alpha-substituted and alpha-unsubstituted beta-fluoroalkyl beta-amino acid derivatives 6 and 9 have been synthesized in two steps starting from fluorinated imidoyl chlorides 1 and ester enolates. This approach is based on the chemical reduction of previously obtained gamma-fluorinated beta-enamino esters 4 by using ZnI(2)/NaBH(4) in a nonchelated aprotic medium (dry CH(2)Cl(2)) as the reducing agent. A metal-chelated six-membered model has been suggested to explain the stereochemical outcome of the reduction reaction. The process takes place with high yields and with moderate to good diastereoselectivity. The best results related to diastereoselective reduction of chiral beta-enamino esters 4 were provided by the use of (-)-8-phenylmenthol as a chiral auxiliary.     

Comparative molecular field analysis of artemisinin derivatives: Ab initio versus semiempirical optimized structures

Based on the belief that structural optimization methods, producing structures more closely to the experimental ones, should give better, i.e. more relevant, steric fields and hence more predictive CoMFA models, comparative molecular field analyses of artemisinin derivatives were performed based on semiempirical AM1 and HF/3-21G optimized geometries. Using these optimized geometries, the CoMFA results derived from the HF/3-21G method are found to be usually but not drastically better than those from AM1. Additional calculations were performed to investigate the electrostatic field difference using the Gasteiger and Marsili charges, the electrostatic potential fit charges at the AM1 level, and the natural population analysis charges at the HF/3-21G level of theory. For the HF/3-21G optimized structures no difference in predictability was observed, whereas for AM1 optimized structures such differences were found. Interestingly, if ionic compounds are omitted, differences between the various HF/3-21G optimized structure models using these electrostatic fields were found.     

A comparative molecular field analysis (CoMFA) study using semiempirical, density functional, ab initio methods and pharmacophore derivation using DISCOtech on sigma 1 ligands

The Comparative Molecular Field Analysis (CoMFA) was developed to investigate a three-dimensional quantitative structure activity relationship (3D-QSAR) model of ligands for the sigma 1 receptor. The starting geometry of sigma-1 receptor ligands was obtained from the Tripos force field minimizations and conformations were decided from DISCOtech using the SYBYL 6.8. program. The structures of 48 molecules were fully optimized at the ab initio HF/3-21G* and semiempirical AM1 calculations using GAUSSIAN 98. The electrostatic charges were calculated using several methods such as semiempirical AM1, density functional B3LYP/3-21G*, and ab initio HF/3-21G*, MP2/3-21G* calculations within GAUSSIAN 98. Using the optimized geometries, the CoMFA results derived from the HF/3-21G method were better than those from AM1. The best CoMFA was obtained from HF/3-21G* optimized geometry and charges (R2 = 0.977). Using the optimized geometries, the CoMFA results derived from the HF/3-21G methods were better than those from AM1 calculations. The training set of 43 molecules gave higher R2 (0.989-0.977) from HF/3-21G* optimized geometries than R2 (0.966-0.911) values from AM1 optimized geometries. The test set of five molecules also suggested that HF/3-21G* optimized geometries produced good CoMFA models to predict bioactivity of sigma 1 receptor ligands but AM1 optimized geometries failed to predict reasonable bioactivity of sigma 1 receptor ligands using different calculations for atomic charges.     

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.     

3-酰基-2,4-吡咯烷二酮及其N-酰基衍生物相对稳定性的理论计算

对 3 酰基 2 ,4 吡咯烷二酮 (1)及其N 酰基衍生物 (2 )的各种异构体 (a -e)分别进行了HF/6 31G和HF/6 31G 水平上的优化计算 ,并考虑电子相关效应 (MP2 /6 31G )对其能量进行了校正 ,探讨了计算水平对体系相对稳定性的影响。HF/6 31G 及MP2 /6 31G 上的计算结果表明 1d和 2a分别为化合物 1和 2的最稳定结构形式 ,与实验观察到的结果一致。同时对N 位取代酰基引起的优势结构的变化给出了较合理的解释。     

Synthesis and Antinociceptive Activity of N -Substituted 4-Aryl-4-oxo-2-[(3-thiophen-2-yl)amino]but-2-enamides

2-{[5-Aryl-2-oxofuran-3(2H)-ylidene]amino}thiophene-3-carboxylic acid derivatives reacted with substituted amines to give new N-substituted 4-aryl-4-oxo-2-[(3-thiophen-2-yl)amino]but-2-enamides. Antinociceptive activity of the synthesized compounds was studied.     

Synthesis of exo- and endocyclic enamino derivatives of 2-(3-arylprop-2-enoyl)cyclohexane-1,3-diones

Condensation of 2-acetylcyclohexane-1,3-diones with aromatic aldehydes in the presence of piperidine, pyrrolidine, hexamethyleneimine, or morpholine gave the corresponding 2-(3-arylprop-2-enoyl)cyclohexane-1,3-diones as the major products and 2-[3-aryl-1-cycloalkylaminoprop-2-en-1-ylidene]cyclohexane-1,3-diones as minor ones. Endocyclic enamino derivatives were synthesized in two steps through enol methyl ethers which reacted with amines. Endocyclic enamino derivatives of 2-(3-arylprop-2-enoyl)cyclohexane-1,3-diones containing an aryl group on the nitrogen atom readily underwent cyclization to 1,2,3,4,5,6,7,8-octahydroquinoline-4,5-dione derivatives.     

The structure, Raman spectra, and force fields of methylbromosilanes CH3SiClnBr3−n (n = 0, 1, 2)

The laser-excited Raman spectra of liquid CH₃SiCl_nBr_3?n (n = 0, 1, 2) were studied. Quantumchemical calculations of these substances with geometry optimization were performed to determine their harmonic force fields and vibrational frequencies. The calculations were made using the HF/6-31G* and HF/6-311++G** approximations and density functional theory at the B3LYP/6-31G* and B3LYP/6-311++G** levels. An interpretation of the spectra was suggested and the calculated force fields were discussed in comparison with the data on related compounds.     

An Ab Initio Study and NBO Analysis of the Stability and Conformational Properties of Hexakis(trimethylelementhyl)benzene (Element = C,Si, Ge,and Sn)

Ab initio molecular orbital and density functional theory were used to investigate energetic and structural properties of the various conformations of hexa-tertbutylbenzene (1), hexakis(trimethylsilyl)benzene (2), hexakis (trimethylgermyl)benzene (3), and hexakis(trimethylstannyl)benzene (4). HF/3-21G//HF/3-21G and B3LYP/3-21G//HF/3-21G results revealed that the Twist-Boat (TB) conformer of compound 1 is more stable than the 1-Chair (C), 1-Boat (B), and 1-Planar (P) conformers. B3LYP/3-21G//HF/3-21G results show that the 1- TB conformer is more stable than 1- C, 1- B, and 1- P conformers of about 1.13, 4.34, and 99.94 kcal mol^?1 , respectively. Contrary to the stability order of compound 1 conformers, the C conformer of compounds 2–4 is more stable than TB, B, and P conformations, as calculated by B3LYP/3-21G//HF/3-21G and HF/3-21G//HF/3-21G levels of theory. The energy gap between the C and P conformers in compounds 1–4 is decreased in the following order: ΔE(4: C, P) < ΔE (3: C, P) < ΔE(2: C, P) < ΔE (1: C, P). This fact can be explained in terms of the increase of C _aromatic -M (M═C, Si, Ge, and Sn) bond lengths and the decrease of steric (van der Waals) repulsions in the previously discussed compounds. For compounds 1–3, the calculations were also performed at the B3LYP/ 6-31G*//HF/3-21G level of theory. However, the comparison showed that the results at B3LYP/3-21G//HF/3-21G methods correlated well with those obtained at the B3LYP/6-31G*// HF/6-31G method. Further, NBO analysis revealed that in compounds 1–4, the resonance energy associated with the σ_M-C1 to σ*_C2-C3 delocalization is 5.20, 9.68, 11.15, and 12.27 kcal mol^?1, respectively. These resonance energy values could explain the easiness of the ring flipping processes of C, B, and TB conformers of compounds 4 to 1. Also, the NBO results showed that by an increase of the σ_M-C1 → σ *_C2-C3 resonance energies in compounds 1–4, the σ_M-C1 bonding orbital occupancies decrease. This fact could fairly explain the increase of the C_aryl-M bond length from compound 1 to 4. The NBO results are also in good agreement with the calculated energy barriers for the ring flipping of the chair conformations in compounds 1–4, as calculated by B3LYP and HF methods.     

Theoretical studies on the gas-phase pyrolysis of 2-phenoxycarboxylic acids: an ONIOM approach

Various ONIOM combinations-ONIOM(HF/6-31G*: PM3), ONIOM(B3LYP/6-31G*: PM3), ONIOM(MP2/6-31G*: PM3), and ONIOM(MP2/6-31G*: HF/3-21G)--were applied to investigate thermal decomposition mechanisms of four 2-phenoxycarboxylic acids (2-phenoxyacetic acid, 2-phenoxypropionic acid, 2-phenoxybutyric acid, and 2-phenoxyisobutyric acid) in the gas phase. All the transition states and intermediates of the reaction paths were optimized. The reaction pathway of four reactants yielding the phenol, CO, and the corresponding carbonyl compound was characterized on the potential energy surface and found to proceed stepwise. The first step corresponds to the elimination of phenol and the formation of alpha-lactone intermediate through a five-membered ring transition state, and the second step is the cycloreversion process of alpha-lactone intermediate to form CO and the corresponding carbonyl compound. The reaction pathway of latter three compounds to produce the carboxylic acid and phenol via a four-membered cyclic transition structure was also examined theoretically. Comparison with experiment indicates that the activation parameters for the fist reaction channel are accurately predicted at the ONIOM(MP2/6-31G*: HF/3-21G) level of theory.     

Computational structural determination and energy landscape analysis of the hepatic carcinogen 2-(acetylamino)fluorene

 2-(Acetylamino)fluorene (AAF), a potent mutagen and a prototypical example of the mutagenic aromatic amines, forms covalent adducts to DNA after metabolic activation in the liver. A benchmark study of AAF is presented using a number of the most widely used molecular mechanics and semiempirical computational methods and models. The results are compared to higher-level quantum calculations and to experimentally obtained crystal structures. Hydrogen bonding between AAF molecules in the crystal phase complicates the direct comparison of gas-phase theoretical calculations with experiment, so Hartree–Fock (HF) and Becke–Perdew (BP) density functional theory (DFT) calculations are used as benchmarks for the semiempirical and molecular mechanics results. Systematic conformer searches and dihedral energy landscapes were carried out for AAF using the SYBYL and MMFF94 molecular mechanics force fields; the AM1, PM3 and MNDO semiempirical quantum mechanics methods; HF using the 3-21G*and 6-31G* basis sets; and DFT using the nonlocal BP functional and double numerical polarization basis sets. MMFF94, AM1, HF and DFT calculations all predict the same planar structures, whereas SYBYL, MNDO and PM3 all predict various nonplanar geometries. The AM1 energy landscape is in substantial agreement with HF and DFT predictions; MMFF94 is qualitatively similar to HF and DFT; and the MNDO, PM3 and SYBYL results are qualitatively different from the HF and DFT results and from each other. These results are attributed to deficiencies in MNDO, PM3 and SYBYL. The MNDO, PM3 and SYBYL models may be unreliable for compounds in which an amide group is immediately adjacent to an aromatic ring. Received: 26 May 2002 / Accepted: 12 December 2002 / Published online: 14 February 2003     

Conformational structures and vibrational spectra of isolated pyrimidine nucleosides: Fourier transform infrared matrix isolation study of 2-deoxyuridine

The conformational structures of 2-deoxyuridine (dU) were investigated using Fourier transform infrared (FTIR) matrix isolation spectroscopy. For the first time the FTIR spectra of dU in Ar matrices were obtained in the range 4000-200 cm(-1). The stabilities of conformers were estimated by the methods HF/3-21G (p), HF/6-31G (d,p) and MP2/6-31G (d,p). Ab initio calculations of the infrared spectra were performed by the methods HF/3-21G (p) and HF/6-31G (d,p). The actual occupancy of conformational isomers in matrix samples was determined. It was shown that anti-conformers of dU are dominant. The ribose rings of the main anti-conformers dU _a0, dU _a1 are in the C2'-endo conformation, but the ribose rings of minor anti-conformers dU_a2, dU_a3 have the C3'-endo conformation, stabilized by intramolecular hydrogen bonds O3'H...O5' and O5'H...O3', accordingly. Syn-conformers of dU are stabilized by the intramolecular hydrogen bond O5'H...O2 and the dominant conformation of the ribose ring is C2'-endo.     

Quantum chemical calculations based on ONIOM and the DFT methods in the inclusion complex: doxycycline/2-O-Me-β-cyclodextrin

A computational study of inclusion complexes of 2-methyl-βCD with Doxycycline tautomeric (enol and keto form) has been performed with several combinations of ONIOM hybrid calculations. The reliability of the ONIOM2 calculations at the integrated level, ONIOM2 (M05-2X/6-31G(d): M05-2X/3-21G*), ONIOM2 (M05-2X/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):HF/3-21G*), ONIOM2 (B3LYP/6-31G(d):B3LYP/3-21G*) and ONIOM2 (B3PW91/6-31G(d):B3PW91/3-21G*) was examined. Their complexation, binding, deformation and stabilization energies, and geometrical data were compared with those of the target geometry structure optimized at the M05-2X/6-31G(d) level of theory. Mixed combinations ONIOM2 (M05-2X 6-31G(d):HF 3-21G*) and ONIOM2 (B3LYP 6-31G(d):HF 3-21G*) reproduces nearly the target geometry structure and provides realistic energetic results at a relatively low computational cost.     

Stable conformations of CH3CH2OCH2CH2OH: A comparison of theoretical methods

This article presents the results of an extensive examination of the stable conformations of CH₃CH₂OCH₂CH₂OH at various levels of theory. In particular, 41 initial conformations are optimized using the MM2 force field in BIGSTRN-3; the MINDO/3, MNDO, and AM1 Hamiltonians in AMPAC 2.2; the PM3 Hamiltonian in MOPAC 7.0; and at the HF/STO-3G and HF/3-21G levels using Gaussian 92. The optimized HF/3-21G structures are reoptimized at the HF/6-31G(d) level, and the unique structures are optimized again at the MP2 = FULL/6-31G(d) level. In addition, single-point MP2/6-31G(d) calculations are performed using the HF/6-31G(d) geometries. The goal is to determine the relative accuracy of each method and discuss their strengths and weaknesses. © 1994 by John Wiley & Sons, Inc.     

Quantum chemical study of the reactivity of C60HR and C60(CHR) derivatives

In the present work a quantum chemical study of a series of substituted hydrofullerenes, C(60)HR, and a series of methanofullerenes, C(60)(CHR), is presented. Their reactivity and geometrical, energetic, electronic, and magnetic properties, as well as the influence of the substituent, are discussed. As a probe of the reactivity, the acidic properties of these fullerene derivatives were predicted, based on the calculated deprotonation energies, with a previously set up scheme. The electronic delocalization upon deprotonation was described, and the global (magnetizabilities) and local aromaticity (nucleus-independent chemical shifts) was analyzed and compared with respect to the group properties for the series of functional groups. The geometries of both acidic and basic forms were fully optimized at the AM1 level, and all property calculations were performed at the HF/3-21G and the B3LYP/6-31G* level of theory.     

Novel synthesis of 2-chloroquinolines from 2-vinylanilines in nitrile solvent

2-Vinyl- or heteroaryl-substituted anilines were reacted with diphosgene in acetonitrile solution via a reactive imidoyl moiety to afford the corresponding 2-chloroquinolines. Facile syntheses of nine 2-chloroquinoline derivatives from several anilines and their postulate mechanism is described. The postulate mechanism of 2-chloroquinoline formation via imidoyl moiety as a good leaving group shows that the reaction consists of the following three steps: (1) generation of phenylisocyanate, (2) quinoline ring formation, and (3) chlorination on C2 position of quinoline.     

Molecular structures of benzoic acid and 2-hydroxybenzoic acid, obtained by gas-phase electron diffraction and theoretical calculations

The structures of benzoic acid (C6H5COOH) and 2-hydroxybenzoic acid (C6H4OHCOOH) have been determined in the gas phase by electron diffraction using results from quantum chemical calculations to inform restraints used on the structural parameters. Theoretical methods (HF and MP2/6-311+G(d,p)) predict two conformers for benzoic acid, one which is 25.0 kJ mol(-1) (MP2) lower in energy than the other. In the low-energy form, the carboxyl group is coplanar with the phenyl ring and the O-H group eclipses the C=O bond. Theoretical calculations (HF and MP2/6-311+G(d,p)) carried out for 2-hydroxybenzoic acid gave evidence for seven stable conformers but one low-energy form (11.7 kJ mol(-1) lower in energy (MP2)) which again has the carboxyl group coplanar with the phenyl ring, the O-H of the carboxyl group eclipsing the C=O bond and the C=O of the carboxyl group oriented toward the O-H group of the phenyl ring. The effects of internal hydrogen bonding in 2-hydroxybenzoic acid can be clearly observed by comparison of pertinent structural parameters between the two compounds. These differences for 2-hydroxybenzoic acid include a shorter exocyclic C-C bond, a lengthening of the ring C-C bond between the substituents, and a shortening of the carboxylic single C-O bond.     

Intermolecular Interaction of HMX： an Application of ONIOM Methodology

IntroductionExtensive studies on intermolecular interac-tions have been made in the past decades due totheir importance in a wide range of physical,chem-ical and biological fields.Researches on the weakintermolecular interactions began with hydrogenbonds.Scheiner summarized the ab initio investiga-tions on hydrogen bonding in detail[1] .With the abinitio method supermolecular structures and bind-ing energies can be predicted notonly for H- bondedsystems but also for other systems[2 _ 5] when …