Theoretical studies on the tautomerism and intramolecular hydrogen shifts of 5-Amino-tetrazole in the gas phase

The tautomerism and intramolecular hydrogen shifts of 5-amino-tetrazole in the gas phase were studied in the present work. The minimum energy path (MEP) information of 5-amino-tetrazole was obtained at the CCSD(T)/6–311G**//MP2/6–311G** level of theory. The six possible tautomers of 1H, 4H-5-imino-tetrazole (a), 1H-5-amino-tetrazole (b), 2H-5-amino-tetrazole (c), 1H, 2H-5-imino-tetrazole (d), the mesoionic form (e) and 2H, 4H-5-imino-tetrazole (f) were investigated. Among these tautomers, there are 2 amino- forms, 3 imino- forms, and 1 mesoionic structure form. In all the tautomers, 2-H form (c) is the energetically preferred one in the gas phase. In the imino- tautomers, the energy value of the compound d is similar as that of the compound f but it is higher than the energy value of the compound a. The potential energetic surface (PES) and kinetics for five reactions have been investigated. Reaction 2 (b→c) was hydrogen shifts only in which the 1-H and 2-H rearrangement. This means that the reaction 2 (b→c) is energetically favorable having an activation barrier of 45.66 kcal·mol⁻¹ and the reaction energies (ΔE) is only 2.67 kcal·mol⁻¹. However, the reaction energy barrier for tautomerism of reaction 1 (b→e) is 54.90 kcal·mol⁻¹. Reaction 1 (b→a), reaction 3 (c→d), and reaction 5 (c→f) were amino- →imino- tautomerism reactions. The energy barriers of amino- →imino- tautomerism reactions required are 59.39, 65.57, 73.61 kcal·mol⁻¹ respectively in the gas phase. The calculated values of rate constants using TST, TST/Eckart, CVT, CVT/SCT and CVT/ZCT methods using the optimized geometries obtained at the MP2/6–311G** level of theory show the variational effects are small over the whole temperature range, while tunneling effects are big in the lower temperature range for all tautomerism reactions. Graphical Abstract Figure (DOC 45.0 KB)     

Initial step of B12-dependent enzymatic catalysis: energetic implications regarding involvement of the one-electron-reduced form of adenosylcobalamin cofactor

Density functional theory analysis was performed to elucidate the impact of one-electron reduction upon the initial step of adenosylcobalamin-dependent enzymatic catalysis. The transition state (TS) corresponding to the Co–C bond cleavage and subsequent hydrogen abstraction from the substrate was located. The intrinsic reaction coordinate calculations predicted that the reaction consisting of Co–C5′ bond cleavage in [Co^III(corrin^•)]–Rib (where Rib is ribosyl) and hydrogen-atom abstraction from the CH₃–CH₂–CHO substrate occurs in a concerted fashion. The computed activation energy barrier of the reaction (15.0 kcal/mol) was lowered by approximately 54.5% in comparison with the reaction involving the positively charged cofactor model (Im–[Co^III(corrin)]–Rib⁺, where Im is imidazole; energy barrier = 33.0 kcal/mol). The Im base was detached during the TS search in the reaction involving the one-electron-reduced analogue. Thus, to compare the energetics of the two reactions, the axial Im ligand detachment energy for the Im–[Co^III(corrin^•)]–Rib model was computed [7.6 kcal/mol (gas phase); 4.6 kcal/mol (water)]. Consequently, the effective activation energy barrier for the reaction mediated by the Im-off [Co^III(corrin^•)]–Rib was estimated to be 22.6 kcal/mol, which implied an overall 31.5% reduction in the energetic demands of the reaction. Considering that the lengthened Co–N_axial bond has been observed in X-ray crystal structure studies of B₁₂-dependent mutases, the catalytic impact induced by one-electron reduction of the cofactor is expected to be higher in the presence of the enzymatic environment.     

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.     

Conformational stabilities, infrared, and vibrational dichroism spectroscopy studies of tris(ethylenediamine) zinc(II) chloride

The conformational stabilities of the transition metal complex of Zn (en)₃Cl₂ were studied using density functional theory (DFT). Deformational potential energy profiles (PEPs), and pathways between the different isomeric conformational energies were calculated using DFT/B3LYP/6–31G. The relative conformational energies of Δ(λλλ), Δ(λλδ), Δ(λδδ) and Δ(δδδ) are 10.48, 7.08, 3.56, and 0.0 kcal/mol, respectively, which are small compared to the barrier heights for reversible phase transitions (49.56, 49.55, 49.52 kcal/mol, respectively). Frequency assignment was carried out by decomposing Fourier transform infrared (FTIR) spectra using Gaussian and Gaussview. The theoretical IR and vibrational dichroism spectroscopy (VCD) absorption spectra are presented for all conformations within the range of 400–3,500 cm^-1.     

Effect of carboxylate and Na ions on the tautomeric status of 9-methylguanine

Interactions of 9-methylguanine (m9Gua) with carboxylate ion of acetic acid (CH3COO-) and Na+ were studied by 1H NMR spectroscopy and ab initio quantum chemical calculations of the B3LYP/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels of theory. Changes in the m9Gua 1H NMR spectrum in the presence of the equimolar amount of sodium acetate (NaAc), which in anhydrous DMSO dissociates to CH3COO- and Na+, were interpreted as a consequence of a complex formation of m9Gua in the amino-keto-N1H tautomeric form (m9GuaN1H) with carboxylate ions via two H-bonds involving amino and N1H-imino protons. Quantum chemical calculations of interactions of the m9GuaN1H ground-state tautomer and the m9GuaN3H high energy one with relative energy 20.01 kcal/mol show that the ground state tautomer forms the ground-state complex CH3COO-:mgGuaN1H, by 5.57 kcal/mol more stable than the CH3COO-:m9GuaN3H complex, and coordination of Na+ with the O6 and N7 atoms reduces this energy difference to 2.57 kcal/mol. Such a coordination of Na+ with tautomer m9GuaN3H therewith decreases its relative energy only to 13.31 kcal/mol. Non-additivity of the two ligands contributions to the 8-times reduction of the relative energy of the high energy tautomer in the CH3COO-:m9GuaN3H:Na+(O6,N7) triple complex was concluded, the role of CH3COO- being dominant. Besides, coordination with Na+ resulted in an iminoproton transfer from the base to CH3COO- in the triple complexes of both tautomers, according to calculations in vacuum. Biological significance of the results is noticed.     

De(side chain) model of epothilone: bioconformer interconversions DFT study

Using ab initio methods, we have studied conformations of the de(sidechain)de(dioxy)difluoroepothilone model to quantify the effect of stability change between the exo and endo conformers of the epoxy ring. The DFT minimization of the macrolactone ring reveals four low energy conformers, although MP2 predicted five stable structures. The model tested with DFT hybride functional (B3LYP/6–31+G(d,p)) exhibits the global minimum for one of the exo forms (C), experimentally observed in the solid state, but unexpectedly with the MP2 electron correlation method for the virtual endo form (W). Using the QST3 technique, several pathways were found for the conversion of the low energy conformers to the other low energy exo representatives, as well as within the endo analog subset. The potential energy relationships obtained for several exo forms suggest a high conformational mobility between three, experimentally observed, conformers. The high rotational barrier, however, excludes direct equilibrium with experimental EC-derived endo form S. The highest calculated transition state for the conversion of the most stable exo M interligand to the endo S form is approximately a 28 kcal/mol above the energy of the former. The two-step interconversion of the exo H conformer to the endo S requires at least 28 kcal/mol. Surprisingly, we found that the transition state energy of the H form to the virtual endo W has the acceptable value of about 9 kcal/mol and the next energy barrier for free interconversion of endo W to endo S is 13 kcal/mol. Figure DFT Model of Epothilone Interconversions Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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     

Quantum chemical investigation of intramolecular thione-thiol tautomerism of 1,2,4-triazole-3-thione and its disubstituted derivatives

The one step intramolecular thione-thiol tautomerism of 1,2,4-triazole-3-thione and its disubstituted derivatives has been studied through the use of electronic structure methods. Due to the absence of experimental data for the parent molecule of 1,2,4-triazole-3-thione the structure and energetics of aforementioned tautomers were derived using various basis sets and levels including HF, B3LYP, and MP2 methods. The gas phase results show that in all different levels of theory the most stable tautomer is the thione form. It has also been revealed that B3LYP/6-31G(d,p) level is quite well suited and reliable to investigate these kinds of tautomerism. To account the influence of substituents on the mentioned tautomerization, the tautomerism and conformational properties as well as vibrational analysis of 20 halophenyl and isopyridyl derivatives were investigated using B3LYP/6-31G(d,p) calculations. In all cases the calculations indicate that substituents have no considerable effects on relative stabilities and energy barriers for the thione-thiol proton transfer and the thione forms are the predominant species in the gas phase. In order to figure out the relative stabilities of the species involved in the tautomerism, geometrical and natural bond orbital (NBO) analyses have been employed. It has also been shown that the computed vibrational frequencies of tautomers with different scaling factors could be used to interpret the vibrational frequencies in IR spectrum of similar species.     

Mechanism of peptide hydrolysis by co-catalytic metal centers containing leucine aminopeptidase enzyme: a DFT approach

In this density functional theory study, reaction mechanisms of a co-catalytic binuclear metal center (Zn1–Zn2) containing enzyme leucine aminopeptidase for two different metal bridging nucleophiles (H₂O and –OH) have been investigated. In addition, the effects of the substrate (l-leucine-p-nitroanilide → l-leucyl-p-anisidine) and metal (Zn1 → Mg and Zn2 → Co, i.e., Mg1–Zn2 and Mg1–Co2 variants) substitutions on the energetics of the mechanism have been investigated. The general acid/base mechanism utilizing a bicarbonate ion followed by this enzyme is divided into two steps: (1) the formation of the gem-diolate intermediate, and (2) the cleavage of the peptide bond. With the computed barrier of 17.8 kcal/mol, the mechanism utilizing a hydroxyl nucleophile was found to be in excellent agreement with the experimentally measured barrier of 18.7 kcal/mol. The rate-limiting step for reaction with l-leucine-p-nitroanilide is the cleavage of the peptide bond with a barrier of 17.8 kcal/mol. However, for l-leucyl-p-anisidine all steps of the mechanism were found to occur with similar barriers (18.0–19.0 kcal/mol). For the metallovariants, cleavage of the peptide bond occurs in the rate-limiting step with barriers of 17.8, 18.0, and 24.2 kcal/mol for the Zn1–Zn2, Mg1–Zn2, and Mg1–Co2 enzymes, respectively. The nature of the metal ion was found to affect only the creation of the gem-diolate intermediate, and after that all three enzymes follow essentially the same energetics. The results reported in this study have elucidated specific roles of both metal centers, the nucleophile, indirect ligands, and substrates in the catalytic functioning of this important class of binuclear metallopeptidases.     

Consequences of one-electron oxidation and one-electron reduction for 4-aminopyrimidine—DFT studies

The consequences of one-electron oxidation and one-electron reduction were studied for 4-aminopyrimidine (4APM), which displays prototropic tautomerism. Since experimental techniques are incapable of detecting less than 0.1% of minor tautomers, quantum-chemical calculations [DFT(B3LYP)/6-311+G(d,p)] were carried out for all possible tautomers of neutral 4AMP and its redox forms, 4APM (+ ?) and 4APM (- ?). Four tautomers were considered: one amine and three imine tautomers (two NH and one CH form). Geometric isomerism of the exo?=?NH group was also taken into account. One-electron oxidation (4APM - e → 4APM (+??)) has no significant effect on the tautomeric preferences; it influences solely the composition of the tautomeric mixture. The amine tautomer is favored for both 4APM (+??) and 4APM. An interesting change in the tautomeric preference occurs for 4APM (- ?). One-electron reduction (4APM?+?e → 4APM (- ?)) favors the C5 atom for the labile proton. The preference of the imine CH tautomer in the tautomeric mixture of 4APM (- ?) may partially explain the origin of CH tautomers in nucleobases.     

Accurate micro- and macro- gas phase basicities of hydroxyl-radical-modified pyrimidines estimated by advanced quantum chemistry methods

The gas phase basicities (GPB) of 16 hydroxyl-radical-derived analogues of model pyrimidine nucleosides were analysed at B3LYP/aug-cc-pVDZ and G3MP2B3 levels of theory. All possible tautomeric equilibriums of neutral and protonated forms of pyrimidine analogues were taken into account. The oxidation of pyrimidines usually reduces the GPB values with respect to canonical values. The only exception was observed in the case of 5,6-dihydroxycytidine, which is more basic than model cytidine if both macro- and micro-scopic measures of basicities are used for comparison. In addition, 6-hydroxycytidine is characterised more by the basic character of the O₂ atom compared to the more basic centre of cytidine, despite the fact that the GPB values of both these compounds are almost identical. Although the B3LYP/aug-cc-pVDZ approach seems to be an accurate and robust method for GPB estimation, the microscopic protonation properties are much more sensitive to this method since the difference in energy between some tautomers is often less than 1 kcal/mol with method-dependent succession. In such cases, G3MP2B3 is recommended.     

Mechanisms on inhibition of polyethylene electrical tree aging: a theoretical study

A theoretical investigation is completed on inhibition mechanism of polyethylene electrical tree aging. Foremost it elucidates that it is one of the important factors for inhibiting initiation and propagation of polyethylene electrical tree through keto-enol tautomerism of acetophenone and its analogues. Geometries of the keto tautomer and the enol tautomer of acetophenone and its analogues, and its transition states are optimized at the B3LYP/6?311+G(d,p) level, the harmonic vibration frequencies of the equilibrium geometries are calculated at the same level. The minimum energy path (MEP) is obtained by the intrinsic reaction coordinate (IRC) theory at the same level. The calculated results show that the energy barriers of keto-enol tautomerism of acetophenone and its analogues are smaller than the average C—C bond energy of polyethylene, the acetophenone and its analogues adulterated in polyethylene composites can improve the strength of alternate current puncture that PE can endure as well as inhibit polyethylene electrical tree from initiation and propagation.     

Assessing the reactivation efficacy of hydroxylamine anion towards VX-inhibited AChE: a computational study

Oximate anions are used as potential reactivating agents for OP-inhibited AChE because of they possess enhanced nucleophilic reactivity due to the α-effect. We have demonstrated the process of reactivating the VX–AChE adduct with formoximate and hydroxylamine anions by applying the DFT approach at the B3LYP/6-311 G(d,p) level of theory. The calculated results suggest that the hydroxylamine anion is more efficient than the formoximate anion at reactivating VX-inhibited AChE. The reaction of formoximate anion and the VX–AChE adduct is a three-step process, while the reaction of hydroxylamine anion with the VX–AChE adduct seems to be a two-step process. The rate-determining step in the process is the initial attack on the VX of the VX–AChE adduct by the nucleophile. The subsequent steps are exergonic in nature. The potential energy surface (PES) for the reaction of the VX–AChE adduct with hydroxylamine anion reveals that the reactivation process is facilitated by the lower free energy of activation (by a factor of 1.7 kcal mol⁻¹) than that of the formoximate anion at the B3LYP/6-311 G(d,p) level of theory. The higher free energy of activation for the reverse reactivation reaction between hydroxylamine anion and the VX–serine adduct further suggests that the hydroxylamine anion is a very good antidote agent for the reactivation process. The activation barriers calculated in solvent using the polarizable continuum model (PCM) for the reactivation of the VX–AChE adduct with hydroxylamine anion were also found to be low. The calculated results suggest that V-series compounds can be more toxic than G-series compounds, which is in accord with earlier experimental observations.     

Elucidation of insulin degrading enzyme catalyzed site specific hydrolytic cleavage of amyloid β peptide: a comparative density functional theory study

In this B3LYP study, the catalytic mechanisms for the hydrolysis of the three different peptide bonds (Lys28-Gly29, Phe19-Phe20, and His14-Gln15) of Alzheimer amyloid beta (Aβ) peptide by insulin-degrading enzyme (IDE) have been elucidated. For all these peptides, the nature of the substrate was found to influence the structure of the active enzyme–substrate complex. The catalytic mechanism is proposed to proceed through the following three steps: (1) activation of the metal-bound water molecule, (2) formation of the gem-diol intermediate, and (3) cleavage of the peptide bond. With the computed barrier of 14.3, 18.8, and 22.3 kcal/mol for the Lys28-Gly29, Phe19-Phe20, and His14-Gln15 substrates, respectively, the process of water activation was found to be the rate-determining step for all three substrates. The computed energetics show that IDE is the most efficient in hydrolyzing the Lys28-Gly29 (basic polar–neutral nonpolar) peptide bond followed by the Phe19-Phe20 (neutral nonpolar–neutral nonpolar) and His14-Gln15 (basic polar–neutral polar) bonds of the Aβ substrate.     

Understanding structural/functional properties of amidase from Rhodococcus erythropolis by computational approaches

The 3D structure of the amidase from Rhodococcus erythropolis (EC 3.5.1.4) built by homology-based modeling is presented. Propionamide and acetamide are docked to the amidase. The reaction models were used to characterize the explicit enzymatic reaction. The calculated free energy barrier at B3LYP/6-31G* level of Model A (Ser194 + propionamide) is 19.72 kcal mol⁻¹ in gas (6.47 kcal mol⁻¹ in solution), and of Model B (Ser194 + Gly193 + propionamide) is 18.71 kcal mol⁻¹ in gas (4.57 kcal mol⁻¹ in solution). The docking results reveal that propionamide binds more strongly than acetamide due to the ethyl moiety of propionamide, which makes the carboxyl oxygen center of the substrate slightly more negative, making formation of the positively charged tetrahedral intermediate slightly easier. The quantum mechanics results demonstrate that Ser194 is essential for the acyl-intermediate, and Gly193 plays a secondary role in stabilizing acyl-intermediate formation as the NH groups of Ser194 and Gly193 form hydrogen bonds with the carbonyl oxygen of propionamide. The new structural and mechanistic insights gained from this computational study should be useful in elucidating the detailed structures and mechanisms of amidase and other homologous members of the amidase signature family.     

The interactions of square platinum(II) complexes with guanine and adenine: a quantum-chemical ab initio study of metalated tautomeric forms

The influence of binding of square planar platinum complexes on tautomeric equilibria of the DNA bases guanine and adenine was investigated using the density functional B3LYP method. Neutral trans-dichloro(amine)-, +1 charged chloro(diamine)-, and +2 charged triamine-platinum(II) species were chosen for coordination to bases. Only the N₇ interaction site of the bases was considered. The calculations demonstrate that the neutral platinum adduct does not change the tautomeric equilibria of the bases. Furthermore, N₇ binding of the neutral Pt adduct moderately reduces the probability of protonation of the N₁ position of adenine. Larger effects can be observed for +1 and mainly +2 adducts, but these can be rationalized by electrostatic effects. Since the electrostatic effects are expected to be efficiently compensated for by a charged backbone of DNA and counterions in a polar solvent, no dramatic increase in mispair formation is predicted for Pt(II) adducts, which is in agreement with experiment. The interaction energies between Pt adducts and the nucleobases were also evaluated. These interaction energies range from ca. 210 kJ/mol for neutral adducts, interacting with both bases and their tautomers, up to 500 kJ/mol for the +2 charged adducts, interacting with the keto-guanine tautomer and the anti-imino-adenine tautomer. The surprisingly large interaction energy for the latter structure is due to the strong H-bond between the NH₃ ligand group of the metal adduct and the N₆ nitrogen atom of the base. Received: 6 July 1999 / Accepted: 7 December 1999     

Quantum chemical studies of the catalytic mechanism of N-terminal nucleophile hydrolase

Modeling of the catalytic mechanism of penicillin acylase, a member of the N-terminal nucleophile hydrolase superfamily, is for the first time conducted at ab initio quantum chemistry level. The uniqueness of this family of enzymes is that their active site lacks His and Asp (Glu) residues, comprising together with a serine residue the classical catalytic triad. The current investigation confirms that the amino group of the N-terminal serine residue in N-terminal hydrolases is capable of activating its own hydroxyl group. Using the MP2/RHF method with the 6−31+G** basis set, stationary points on the potential energy surface of the considered molecular system were located, corresponding to local minima (complexes of reagents, products, intermediate) and to saddle points (transition states). It turned out that the stage of acyl-serine formation proceeds via two transition states; the first one, which separates reagents from the so-called tetrahedral intermediate, has the highest relative energy (30 kcal/mol). In contrast to recently proposed empiric suggestions, we have found that participation of a bridging water molecule in proton shuttling is not necessary for the catalysis. The quantum chemical calculations showed a crucial role of a specific solvation in decreasing the activation barrier of the reaction by approximately 10 kcal/mol. Published in Russion in Biokhimiya, 2007, Vol. 72, No. 5, pp. 615–621.     

A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances

Abstract

In this study for the first time we have revealed and investigated in details 123 different prototropic tautomers of the most stable conformer of the quercetin molecule using quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of QM theory. We have found that in the most energetically favorable prototropic tautomer mobile hydrogen atoms are localized at the О3, О3′, О4′, О5, and О7 exocyclic oxygen atoms. Molecular tautomers are in the range of the Gibbs free energies from 0.0 to 69.8?kcal·mol^?1, while zwitterionic ones – from 30.1 до 172.8?kcal·mol^?1 at normal conditions. It was also reliably established that the weakest point causing the decyclization of the molecule is its C ring – this reaction is launched by the transition of the proton from the C8H group to the endocyclic O1 oxygen atom. All prototropic tautomers, except two cases, are joined by the intramolecular cooperative specific interactions (from 1 to 5) – H-bonds and attractive van der Waals contacts, which have been revealed and characterized by QTAIM analysis.

Communicated by Ramaswamy H. Sarma     

The tri–μ–hydrido–bis[(η5–C5Me5)aluminum(III)] theoretical study, the assets of sandwiched M2H3 (M of 13th group elements) stability

The stability of the tri–μ–hydrido–bis[(η⁵–C₅Me₅)aluminum], Cp*₂Al₂H₃, 1 is studied at B3LYP/6–311+G(d,p), CCSD(T)//B3LYP/6–311+G(d,p) and MP4//B3LYP/6–311+G(d,p) levels. The coordination between Al₂H₃ entity and both C₅(CH₃)₅ groups is ensured by strong electrostatic and orbital interactions. The orbital analysis of the interacting fragments shows that Al₂H₃ acceptor, which keeps its tribridged structure, implies the vacant ( \texta₁￠ ) \left( {{\text{a}}_1^\prime } \right) and five antibonding (a₂^￠￠ a_2^{\prime \prime } , e′ and e″) molecular orbitals to interact with two orbitals mixtures, b₁ and e" of the donors (C₅Me₅). When we take into account the solvent effect, the computation shows that 1 seems to be stable in condensed phase with a tribridged bond between the Al atoms [Cp*Al(μ-H)₃AlCp*], whereas in the gas phase, the monobridged Cp*AlH(μ-H)AlHCp* 4 is slightly favored (4 kcal mol⁻¹). We propose that 1 could be prepared thanks to Cp*Al (2) and Cp*AlH₂ (3) reaction in acidic medium. The experimental treatment of this type of metallocenes would contribute to the development of the organometallic chemistry of 13th group elements.      

Conformational analysis of perezone and dihydroperezone using vibrational circular dichroism

The vibrational circular dichroism (VCD) spectra of perezone and dihydroperezone measured from CDCl₃ solutions were quite similar, suggesting analogous conformations for both molecules. Their absolute configurations were confirmed by comparison of the experimental VCD spectrum of each compound with curves generated from theoretical calculations using density functional theory (DFT) at the B3LYP/DGDZVP level of theory taking into account their conformational mobility. Conformational analysis of the 8-(R) enantiomer showed 19 low energy conformers in a 2.4 kcal/mol energy range, while for 8-(R), with the saturated side alkyl chain, 34 conformers were considered in the first 2 kcal/mol. Initial analyses were carried out using a Monte Carlo searching with the MMFF94 molecular mechanics force field, all MMFF94 conformers were geometrically optimized using DFT at the B3LYP/6-31G(d) level of theory, followed by reoptimization and calculations of their vibrational frequencies at the B3LYP/DGDZVP level. Good agreement between the theoretical 8-(R) enantiomers and experimental VCD curves were observed for both.