Absolute configuration of (−)-myrtenal by vibrational circular dichroism

The VCD spectrum of the monoterpene (−)-myrtenal (1) was compared with theoretical spectra using ab initio density functional theory (DFT) calculations at the B3LYP/6-31G(d,p), B3LYP/6-31G+(d,p), B3LYP/6-311G+(d,p), B3LYP/DGDZVP, and B3PW91/DGTZVP levels of theory. Conformational analysis of 1 indicated that the lowest energy conformer was s-trans-C2-C10, which contributes more than 98.5% to the total conformational population regardless of the employed level of theory. The use of a recently developed confidence level algorithm demonstrated that VCD spectra calculated for the main conformer, using the indicated hybrid functionals and basis set, gave no significant changes, from where it follows that B3LYP/DGDZVP calculations provide a superior balance between computer cost and VCD spectral accuracy. The DGDZVP basis set demanded around a quarter the time than the 6-311G+(d,p) basis set while providing similar results. The spectral comparison also provided evidence that the levorotatory enantiomer of myrtenal has the 1R absolute configuration.     

Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.     

Racemization of the gastrointestinal antisecretory chiral drug esomeprazole magnesium via the pyramidal inversion mechanism: A theoretical study

The pyramidal inversion mechanisms of the 6‐methoxy and the 5‐methoxy tautomers of (S)‐omeprazole were studied, employing ab initio and DFT methods. The conformational space of the model molecule (S)‐2‐[(3‐methyl‐2‐pyridinyl)methyl]sulfinyl‐1H‐benzimidazole was calculated, with respect to rotations around single bonds, at the B3LYP/6‐311G(d,p) level. All of the resulting conformations were used as starting points for full optimizations of (S)‐omeprazole, at B3LYP/6‐31G(d), B3LYP/6‐311G(d,p), B3LYP/6‐311++G(d,p), B3LYP/6‐311G(2df,2pd), MP2/6‐31G(d), and MP2/6‐311G(d,p) levels. Four distinct pathways were found for enantiomerization via the pyramidal inversion mechanism for each of the tautomers of (S)‐omeprazole. Each transition state, in which the sulfur, the oxygen and the two carbon atoms connected directly to the sulfur are in one plane, connects two diastereomeric minima. The enantiomerization is completed by free rotation around the sulfur–methylene bond, and around the methylene–pyridine ring bond. The effective Gibbs' free energy barrier for racemization ΔG of the two tautomers of (S)‐omeprazole are 39.8 kcal/mol (5‐methoxy tautomer) and 40.0 kcal/mol (6‐methoxy tautomer), indicating that the enantiomers of omeprazole are stable at room temperature (in the gas phase). The 5‐methoxy tautomer of (S)‐omeprazole was found to be slightly more stable than the 6‐methoxy tautomer, in the gas phase. The energy barrier (ΔG^?) for the(S,M) (S,P) diastereomerization of (S)‐omeprazole due to the rotation around the pyridine chiral axis was very low, 5.8 kcal/mole at B3LYP/6‐311G(d,p). Chirality 2010. © 2010 Wiley‐Liss, Inc.     

Ionic interactions for substituted MCH1R inhibitors studied by pKa values

MCH1R inhibitors with the quinoline moiety having the aromatic amine and aliphatic amine chain were selected, and then the effect of substituents of the quinoline ring on the ionic interaction were studied by calculating pK_a values for these amines at the B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d) level in the gas phase and in water. For substituent with C, N, O, and S atoms next to the quinoline ring, respectively, the pK_a values of aromatic amines are estimated to be 8.98, 12.19, 4.64, and 4.33 and those of the aliphatic amines are 12.65, 10.82, 9.94, and 11.55, respectively.     

An assessment to evaluate the validity of different methods for the description of some corrosion inhibitors

New research and development efforts using computational chemistry in studying an assessment of the validity of different quantum chemical methods to describe the molecular and electronic structures of some corrosion inhibitors were introduced. The standard and the highly accurate CCSD method with 6-311++G(d,p), ab initio calculations using the HF/6-31G++(d,p) and MP2 with 6-311G(d,p), 6-31++G(d,p), and 6-311++G(2df,p) methods as well as DFT method at the B3LYP, BP86, B3LYP*, M06L, and M062x/6-31G++(d,p) basis set level were performed on some triazole derivatives and sulfur containing compounds used as corrosion inhibitors. Quantum chemical parameters, such as the energy of the highest occupied molecular orbital energy (E_HOMO), the energy of the lowest unoccupied molecular orbital energy (E_LUMO), energy gap (ΔE), dipole moment (μ), sum of total negative charges (TNC), chemical potential (Pi), electronegativity (χ), hardness (η), softness (σ), local softness (s), Fukui functions (f ⁺,f ^?), electrophilicity (ω), the total energy change (?E_T) and the solvation energy (S.E), were calculated. Furthermore, the accuracy and the applicability of these methods were estimated relative to the highest accuracy and standard CCSD with 6-311++G(d,p) method. Good correlations between the quantum chemical parameters and the corresponding inhibition efficiency (IE%) were found.     

Studies on molecular structure and tautomerism of a vitamin B<Subscript>6</Subscript> analog with density functional theory

This work presents a computational study on the molecular structure and tautomeric equilibria of a novel Schiff base L derived from pyridoxal (PL) and o-phenylenediamine by using the density functional method B3LYP with basis sets 6-31 G(d,p), 6-31++G(d,p), 6-311 G(d,p) and 6-311++G(d,p). The optimized geometrical parameters obtained by B3LYP/6-31 G(d,p) method showed the best agreement with the experimental values. Tautomeric stability study of L inferred that the enolimine form is more stable than its ketoenamine form in both gas phase and solution. However, protonation of the pyridoxal nitrogen atom (LH) have accelerated the formation of ketoenamine form, and therefore, both ketoenamine and enolimine forms could be present in acidic media.     

Computational study of decomposition mechanisms and thermodynamic properties of molecular-type cracking patterns for the highly energetic molecule GZT

This study uses the Gaussian 03 program and density functional theory B3LYP with three basis set methods—[B3LYP/6-311+G(d,p), B3LYP/6-31+G(2d,p), and B3LYP/6-31G(d,p)]—to model the highly energetic ionic compound diguanidinium 5,5′-azotetrazolate (GZT) to research its decomposition mechanisms and thermodynamic properties. Molecular-type cracking patterns are proposed, which were initiated by heterocyclic ring opening, sequential cracking of the two five-membered rings of GZT, and simultaneous release of N₂ molecules; whereas proton transfer, bond-breaking, and atomic rearrangements were performed subsequently. Finally, 15 reaction paths and five transition states were obtained. All possible decomposition species and transition states, including intermediates and products, were identified, and their corresponding enthalpy and Gibbs free energy values were obtained. The results revealed that (1) the maximum activation energy required is 187.8 kJ mol^–1, and the enthalpy change (ΔH) and Gibbs free-energy change (ΔG) of the net reaction are ?525.1 kJ mol^–1 and ?935.6 kJ mol^–1, respectively; (2) GZT can release large amounts of energy, the main contribution being from the disintegration of the 5,5'-azotetrazolate anion (ZT^2?) skeleton (ΔH?=??598.3 kJ mol^–1); and (3) the final products contained major amounts of N₂ gas, but remaining gas molecules such as HCN and NH₃ were obtained, which are in agreement with experimental results. The detailed decomposition simulation results demonstrated the feasibility of this method to calculate the energies of the thermodynamic reactions for the highly energetic GZT and predict the most feasible pathways and the final products.     

Ab initio and DFT study of the inner mechanism and dynamic stereochemistry of electrophilic addition reaction of bromine to bisbenzotetracyclo[6.2.2.23,6.02,7]tetradeca-4,9,11,13-tetraene

The inner mechanism and dynamic stereochemistry of electrophilic addition of bromine to bisbenzotetracyclo[6.2.2.2(3,6).0(2,7)]tetradeca-4,9,11,13-tetraene(BBTT) molecule have been investigated by the methods of quantum chemistry. The structure of the BBTT molecule has been studied by ab initio and DFT/B3LYP methods using the 6-31G(d) and 6-311G(d) basis sets. The double bonds of BBTT molecule are endo-pyramidalized. The structure and stability of the cationic intermediates and products of the addition reaction have been investigated by HF/6-311G(d), HF/6-311G(d,p), B3LYP/6-311G(d) and B3LYP/6-311++G(2d,p)//B3LYP/6-311G(d) methods. The bridged bromonium cation isomerized into the more stable nonclassical delocalized N- and U-type cations and the difference between the stability of these cations is small. For the determination of the direction of addition reaction and the stereochemistry of the products, the stability of nonclassical delocalized N- and U-type ions and the structure of their cationic centres play a vital role. Since the cationic centre of the N-type ion is in interaction with the benzene ring from the exo face, the nucleofilic attackof the bromide anion to this centre occurs from the endo face and the exo,endo-isomer of the N-type product is obtained. The attack of bromide anion, towards the cationic centre of U-type ion from the endo face is sterically hindered by the hydrogen atom therefore the attack occurs from the exo face, which interacts with the benzene ring and the more stable exo,exo-isomer of U-type product is formed. Although, the U-type cation was 2.232 kcal mol(-1) more stable than the N-type cation, the U-type product was 0.587 kcal mol(-1) less stable than the N-type product.     

The (α-1,6) glycosidic bond of isomaltose: a tricky system for theoretical conformational studies

Stable conformations of β-isomaltose (α-d-glucopyranosyl-(1→6)-β-d-glucose) in gas-phase and aqueous solution are investigated in this study using quantum mechanical calculations. Conformational maps are calculated at HF/6-31G(d,p) level and lower energy structures are sampled in the most stable regions. Entropic and thermal corrections are considered and the Boltzmann population is obtained for conformers that are representative of the 18 most stable regions found on the potential energy surface. B3LYP/6-31+G(d,p) and B3LYP/6-311+G(2d,2p) calculations are used in conformational samplings. Solvation effects are considered through the polarizable continuum model approach. Hydroxymethyl group orientations are investigated for the most stable conformers. The influence of electronic correlation and solvation on the glycosidic linkage preference (TG, GT, and GG) and hydroxymethyl group orientation (tg, gt, and gg) are discussed. Heteronuclear spin coupling constants (³J_C,H) along the glycosidic linkage are calculated and comparison with other theoretical results and experiments is used to validate the obtained structures.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

Molecular structure,vibrational spectra and HOMO,LUMO analysis of 4-piperidone by density functional theory and ab initio Hartree–Fock calculations

Vibrational frequencies and geometrical parameters of 4-piperidone (4-PID) in the ground state have been calculated by using the Hartree–Fock (HF) and density functional methods (B3LYP) with 6-311++G(d,p) and 6-311+G(3df,2p) basis sets. These methods are proposed as a tool to be applied in the structural characterisation of 4-PID (C₅H₉NO). The title molecule has C _s point group symmetry, thus providing useful support in the interpretation of experimental IR and Raman data. The DFT-B3LYP/6-311+G(3df,2p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 4-PID. The calculated highest occupied molecular orbital and lowest unoccupied molecular orbital energies show that charge transfer occurs within the molecule. The theoretical spectrograms for FT-IR and FT-Raman spectra of the title molecule have been constructed.     

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

Structure and electronic properties of azadirachtin

We performed a combined DFT and Monte Carlo ¹³C NMR chemical-shift study of azadirachtin A, a triterpenoid that acts as a natural insect antifeedant. A conformational search using a Monte Carlo technique based on the RM1 semiempirical method was carried out in order to establish its preferred structure. The B3LYP/6-311++G(d,p), wB97XD/6-311++G(d,p), M06/6-311++G(d,p), M06-2X/6-311++G(d,p), and CAM-B3LYP/6-311++G(d,p) levels of theory were used to predict NMR chemical shifts. A Monte Carlo population-weighted average spectrum was produced based on the predicted Boltzmann contributions. In general, good agreement between experimental and theoretical data was obtained using both methods, and the ¹³C NMR chemical shifts were predicted highly accurately. The geometry was optimized at the semiempirical level and used to calculate the NMR chemical shifts at the DFT level, and these shifts showed only minor deviations from those obtained following structural optimization at the DFT level, and incurred a much lower computational cost. The theoretical ultraviolet spectrum showed a maximum absorption peak that was mainly contributed by the tiglate group.     

σ*-Aromaticity in Three-Membered Rings

The structures and energies of unsaturated three-membered rings of the general formula (CH)₂XY_n^m, with charge = m and n substituents Y at X (Y = fluorine, chlorine, bromine, hydrogen, X = phosphorus, silicon, nitrogen, carbon) are compared to their saturated analogs. The structures were optimized with B3LYP/6-311+G(2d,p) and at MP2/6-31+G(d), with single point energy calculations on the latter geometries at MP4SDTQ/6-31+G(d). The geometrical changes in bond lengths and angles, which correlate with substituent electronegativities, are discussed for the different ring systems. The relative stabilities of unsaturated and saturated rings are compared using isodesmic ring-opening reactions and homodesmic substituent-exchange reactions. C*-Aromaticity, a hyperconjugative effect found in the disubstituted rings, causes lowering of ring strain energies for the unsaturated rings and preference of unsaturated rings over saturated ones for the more electronegative substituents. For the mono-substituted ?-aromatic silacyclopropenes and cyclopropenes, a destabilization by more electronegative ligands is found. For the neutral rings monosubstituted at main group V atoms like the 1H-phosphirenes and also the isoelectronic negatively charged rings with main group IV atoms like the silacyclopropenium anions, no correlation of stabilization energies or geometrical changes with ligand electronegativity is found.     

Density functional computations of enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes

The enantioselective alkynylation of aldehyde catalyzed by chiral zinc(II)-complexes was studied by means of the density functional theory (DFT). All the structures were optimized completely at the B3LYP/6-31G(d,p) level. To obtain more exact energies, single-point energy calculations at B3LYP/6-31+G(d,p) level were carried out on the B3LYP/6-31G(d,p) geometries. As shown, this enantioselective alkynylation was endothermic. The chirality-determining step for the alkynylation was the formation of the catalyst–ethanol complexes and the transition states for this step involved a six-membered ring. The dominant products predicted theoretically were of (R)-chirality, in good agreement with experiment.      

Bis-dibenzo[a.i]fluorenylidene, does it exist as stable 1,2-diradical?

The geometries, energies, and electronic properties of the two possible configurations of bis-[dibenzo[a.i]fluorenylidene] were investigated theoretically by density functional theory DFT B3LYP at the UB3LYP/6-311?+?G(2d,p) // UB3LYP/6-31?+?G(d,p) level of theory. According to the performed calculations, it was found that the singlet is 3.4?kcal?mol^-1 lower in energy compared to triplet state at room temperature. This gap is compared with those of other alkenes like ethylene, (61.9?kcal?mol^-1) tetra-tert-butyethylene, (6.4?kcal?mol^-1) and bis-fluorenylidene (19.5?kcal?mol^-1). These results confirm the experimental findings of the paramagnetic properties determined by Franzen and Joschek. The low singlet-triplet gap in the case of bis-[dibenzo[a.i]fluorenylidene] is the result of a steric destabilization of the singlet due to strain and stabilization of the triplet electronic state by delocalization of each free electron within each aromatic moiety. This correlates with the special electronic structure of the triplet state of this compound, where facial interaction of two hydrogen atoms lying close to the lobes of each p-orbital occupied with a single electron at the distorted double bond in the triplet electronic state.

Conformational preferences of α,α-trehalose in gas phase and aqueous solution

This work presents an investigation on the conformational preferences of α,α-trehalose in gas phase and aqueous solution. Eighty-one systematically selected structures were studied at the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d) level, giving rise to 40 unique conformers. The 19 lower energy structures and some selected other were further re-optimized at the B3LYP/6-311++G(d,p) level. The main factors accounting for the conformer’s stability were pointed out and discussed. NBO and QTAIM analyses were performed in some selected conformers in order to address the anomeric and exo-anomeric effects as well as intramolecular hydrogen bonding. The effect of solvent water on the relative stability of the conformers was accounted for by applying the conductor-like polarizable continuum model, CPCM.     

Structure,energetics and properties of some molecules with potent anti-HIV activity: a theoretical study

Density functional theory (DFT; B3LYP) and Hartree–Fock (HF; 3-21G, 6-31G(d) and 6-311G(d,p)) calculations with complete geometry optimisations are carried out in the ground state on five 6-aminoquinolone derivatives, which have been proved to be highly effective in inhibiting HIV replication, to study their structures, energetics and HOMO–LUMO correlation with physiological action. The gas-phase calculations and single-point polarisable continuum model water-phase calculations show that the molecules are highly effective in inhibiting HIV replication, which is in excellent agreement with the experiment. Structural features, energies, charge densities and HOMO–LUMO correlation have been found to substantiate the experimental findings. Compound 4 (pyrazine) shows some special features in DFT calculations which are not found in HF calculations. In the present series, HF results are more reliable as expected.     

ONIOM DFT/PM3 calculation on the interaction between STI-571 and abelson tyrosine kinase

The ONIOM2 (B3LYP/6–31G (d, p): PM3) and B3LYP/6–31G (d, p) methods were applied to investigate the interaction between STI-571 and abelson tyrosine kinase binding site. The complex of N-[4-methyl-3-(4-pyridin-3-yl-pyrimidin-2-ylamino)- phenyl]-benzamide (part of STI-571) and related 16 amino acid residues were found at B3LYP/6–31G (d, p) level to have hydrogen bonds and π....π stacking interaction, their binding energy via HAF optimization was −20.4 kcal mol⁻¹. The results derived from this study agreed well with the reported observation. Figure Optimized structure of STI-571 and Thr315 in abelson tyrosine kinase based on ONIOM2 method