The geometric, electronic structure, effect of the substitution, and structure physical-chemistry relationship for oxazoles derivatives have been studied by ab initio and DFT theory. In the present work, the calculated values, namely, net charges, bond lengths, dipole moments, electron affinities, heats of formation, and QSAR properties are reported and discussed in terms of the reactivity of oxazole derivatives.

1. Introduction

The emerging resistance to antimicrobial drugs demands the synthesis of new remedies for microbial infections, which are effective against organisms resistant to currently available drugs [1].

Oxazoles and isoxazoles have an important role in the synthesis of potential medicines including antitumor, antimicrobial, anti-infective, cardiovascular, and nervous system agents [2].

There are various experimental methods developed for the synthesis of oxazole and isoxazole derivatives [3, 4]. But there are few relatively theoretical works on oxazoles and isoxazoles [57].

Quantum chemistry methods play an important role in obtaining molecular geometries and predicting various properties. To obtain highly accurate geometries and physical properties for molecules that are built from electronegative elements, expensive ab initio/HF electron correlation methods are required [8]. Density functional theory methods offer an alternative use of inexpensive computational methods which could handle relatively large molecules [9].

Quantitative structure-activity relationships (QSAR) are attempts to correlate molecular structure, or properties derived from molecular structure, with a particular kind of chemical or biochemical activity. The kind of activity is a function of the interest of the user: QSAR is widely used in pharmaceutical, environmental, and agricultural chemistry in the search for particular properties. The molecular properties used in the correlations relate as directly as possible to the key physical or chemical processes taking place in the target activity [10].

In present case, we have studied the structure of oxazole Figure 1 and some of oxazole derivatives, by using ab initio method. For a complete and comparative study, we have also used DFT method.

2. Computational Methods

All calculations were performed by using HyperChem 8.03 software [11]. The geometries of oxazole and its methyl, cyanide derivatives and the series of oxazole derivatives were first fully optimized by molecular mechanics, with MM force field (rms = 0.001 Kcal/Å).

We also used the molecular dynamics for the conformational research, with the following options: 1000 K, in vacuo, steep size: 0.001 ps, relaxation time: 0.1 ps. Further, geometries were fully reoptimized by using ab initio/HF (6-31G**) [12].

A parallel study has been made by DFT/B3LYP exchange-correlation potential with 6-31G** basis [13].

The calculation of QSAR properties is performed by the module (QSAR Properties, version 8.0). QSAR Properties is a module that, together with HyperChem, allows several properties commonly used in QSAR studies to be calculated. The calculations are empirical and so, generally, are fast. The calculated results have been reported in the present work.

3. Results and Discussion

3.1. Geometric and Electronic Structure of Oxazole and Oxazole Systems

The efficiency of DFT/B3LYP method may be scrutinized by comparison with the results obtained by more elaborate calculation such as ab initio/HF. Present results concern bond length values for oxazole (Table 1) and charge densities (Table 2). We can note a good correlation between the ab initio and DFT for bond lengths. Charge densities calculated by the ab initio are similar than DFT method. The geometry of the oxazole is almost planar. The Dihedral angles of this molecule vary between 0.001 and 8.426 degree (Table 3). The lower deviations of registered dihedral angles compared with normal values were imposed essentially by a cyclic chain [7].

The calculated values of methyl-substituted oxazoles and cyanide-substituted oxazoles (Figure 2) are given in Tables 47. In Tables 4 and 5, heat of formation, dipole moment, HOMO (highest occupied molecular orbital), LUMO (lowest unoccupied molecular orbital), and their difference (ΔE) are reported for oxazole and its methyl and cyanide derivatives. In Tables 6 and 7, net atomic charges are reported also.

Series 1:(1)R1=R2= R3=H(2)R1=CH3, R2=R3=H(3)R1=R3= H, R2=CH3(4)R1=R2=H, R3=CH3(5)R1=R2=CH3, R3=H(6)R1=R3=CH3, R2=H(7)R1=H, R2=R3=CH3(8)R1=R2=R3=CH3

Series 2:(1)R1=R2=R3=H(2)R1=CN, R2=R3=H(3)R1=R3=H, R2=CN(4)R1=R2=H, R3=CN(5)R1=R2=CN, R3=H(6)R1=R2=CN, R3=H(7)R1=H, R2=R3=CN(8)R1=R2=R3=CN

It can be seen from the heat of formation data that approximately 9 kcal/mol is increased at each addition of methyl group, in the base compound oxazole irrespective of the number of substitutions.

The ionization potential values in compounds 18 show a decreasing trend which depicts increasing trend in the easy flow of charges in higher energy states of these compounds. Oxygen and nitrogen contribute eight and seven electron density of oxazoles, respectively.

The negative atomic charge on oxygen is increased considerably for methyl derivatives, but on nitrogen is enhanced except for compounds 4 (Table 6). In the monosubstituted methyl group category, the 5-methyl oxazole (compound 4) shows maximum charge on 2th position carbon (0.390) which leads to nucleophilic substitution (Table 6). This is further supported by the least HOMO-LUMO energy gap (13.76) (Table 4) which depicts the chemical reactivity of the compound; the higher is the HOMO-LUMO energy gap, the lesser is the flow of electrons to the higher energy state, making the molecular hard and less reactive.

On the other hand in lesser HOMO-LUMO gap, there is easy flow of electrons to the higher energy state making it softer and more reactive (HSAB principle: hard and soft acids and bases). Hard bases have highest occupied molecular orbitals (HOMO) of low energy, and hard acids have lowest unoccupied molecular orbitals (LUMO) of high energy [14].

Compound 4 also shows maximum dipole moment value. These results are in close agreement with the experiment [15].

In the case of dimethyl-substituted, oxazole the C-2-position (compound 6) shows maximum charge (0.564), least HOMO-LUMO energy gap (13.58), and high dipole moment value (Tables 46) which leads to preferential site of nucleophilic attack.

This conclusion finds support from experimental evidence. In search of basicity, N atom is predicted to be the main basic centre of the oxazole systems in accordance with the electron densities (Table 4).

The C–H hyperconjugation is the principal mode of electron release by the methyl group (pseudohetero atom) and stabilizes excited states more than ground state [16]. In the order of increasing number of conjugated methyl groups, ionization potentials (IPs) decreases in the case of compounds 1–8 as expected from those listed in Table 6.

The 2, 4, 5-trimethyl oxazole (compound 8) is predicted to be the most reactive with least HOMO-LUMO energy gap of all the oxazole systems and, respectively, C2, C4, and C5 are the most preferential sites for nucleophilic attack (Table 4).

These results are in close agreement with the experiment [15], and we found in our literature that the majority of trisubstituted oxazoles have biological activity [3, 4, 10, 15].

We note also that the methyl substituent (donor effect) has the effect of increasing the energy of the HOMO, with little change of the LUMO (Table 4).

In the present work, we have studied cyanide-substituted oxazoles along the same line of methyl-substituted oxazoles whish is for a comparative study. It is interesting to note that the heat of formation approximately 36 Kcal/mol is increased, for each addition of cyanide group irrespective of oxygen in the ring. For nitrogen atom, the negative atomic charge is considerably weakened for these cyanide derivatives (Table 7).

In monosubstituted cyanide derivatives, 2-cyano oxazole (compound 2) is predicted to be more chemically reactive than 5-cyano oxazole and 4-cyano oxazole on the basis of least HOMO-LUMO energy gap (Table 5).

The carbon C-2 in 2-cyano oxazole shows maximum positive charge (0.6101) leading to favored site for nucleophilic attack (Table 6). In disubstituted cyanide derivatives, 2,5-dicyano oxazole (compound 6) seems to be more reactive than the other two composites 7 and 5, due to least HOMO-LUMO energy gap (11.97) (Table 5).

The compound 8 is predicted to be the most reactive with least HOMO-LUMO energy gap of all oxazole systems. We note also that the cyanide substituent (attractor effect) lowers the energies of HOMO and LUMO. His influence on the energy of the LUMO is more important (Table 5).

3.2. Study of Structure Physical-Chemistry Property Relationship for Oxazole Derivatives

We have studied six physical-chemical properties of oxazole derivatives (ten compounds) by HyperChem software, for example, in Figure 3, the favored conformation in 3D of the compound 4. We will continue this work in the future by a quantitative calculation. QSAR properties are van der Waals-surface-bounded molecular volume, the log of the octanol-water partition coefficient ( ), polarizability, and solvent-accessible surface-bounded molecular volume and molecular mass (M). Calculation of is carried out using atomic parameters derived by Viswanadhan and coworkers [17].

Computation of molar refractivity was made via the same method as . Ghose and Crippen presented atomic contributions to the refractivity [18]. The solvent-accessible surface bounded molecular volume and van der Waals-surface-bounded molecular volume calculations are based on a grid method derived by Bodor et al. [19], using the atomic radii of Gavezzotti [20]. The polarizability was estimated from an additivity scheme given by Miller [21] with a precision on the calculation of 3%, where different increments are associated with different atom types. The hydration energy is a key factor determining the stability of different molecular conformations [22].

The calculation is based on exposed surface area and employs the surface area as computed by the approximate method (above), weighted by atom type.

3.2.1. Structural Comparison of the Oxazole Derivatives

Based on our conclusions on the effect of substitution on the oxazole molecule. We chose a series of oxazole derivatives, some having a biological activity. This series of oxazole derivatives has been synthesized and characterized by Shablykin et al. [4, 23] and Rybakov et al. [24]. Initially, we performed a structural comparison of this series (Figure 4). These molecules have a weak conformational flexibility, with regard to the other macrocycles of macrolide type [2531]. In a window of 2 kcal/mol (MM), only one favored conformations are found, for each structure.

3.2.2. Structure Physical-Chemistry Property Relationship

The values of polarizability are generally proportional to the values of surfaces and of volumes; the decreasing order of polarizability for these studied oxazoles is 3, 10, 9, 8, 5, 2, 6, 7, 4, and 1 (Table 8). The order of polarizability is the same one for volumes, but it is a bit little different for surfaces; that is due primarily to the folding up of surfaces of some structures compared to the majority of the extended structures. This is also explained by the relation between polarizability and volume, for the relatively nonpolar molecules.

The polarizability of the molecule depends only on its volume; the thermal agitation of the nonpolar molecules does not have any influence on the appearance of dipole moments in these molecules. On the other hand for the polar molecules, the polarizability of the molecule does not depend solely on volume but also depends on other factors such as the temperature, because of the presence of the permanent dipole [32].

The surface and the volume of distribution of these molecules are definitely higher than that of more polar molecules like the lipopeptides or beta-lactams. For example, Deleu et al. used TAMMO software [33] on the surfactins C13, C14, and C15 having cores similar to the macrolides. They found that their surfaces vary from 129 to 157 Å2 [34]; contrarily for these oxazoles, surfaces vary from 409 to 629 Å2. These oxazoles have a great variation of distribution volume, in particular compound 3 and compound 10 which have respective volumes: 1080 and 960 Å3 (Table 8).

The hydration energy in the absolute value, most important is that of the compound 6 (17.7 kcal/mol) and the weakest is that of compound 3 (5.6 kcal/mol) (Table 8). Indeed, in the biological environments the polar molecules are surrounded by water molecules. They are established hydrogen bonds between a water molecule and these molecules. The donor sites of proton interact with the oxygen atom of water, and the acceptor sites of proton interact with the hydrogen atom.

The first corresponds to the complex with the strongest hydrogen bond. These hydrated molecules are dehydrated at least partially before and at the time of their interaction. These interactions of weak energy, which we observe in particular between messengers and receivers, are generally reversible [35]. Compound 6 has a site donor of proton (1NH on the alkyl group/R3) and seven acceptor sites of proton (5N on the alkyl groups and 1N, 1O on the principal cycle). On the other hand, the compound 3 does not possess any donor site but possess five acceptor sites of proton (3N, on the alkyl groups, 1N and 1O on the principal cycle,) With the first having higher value, it has two more acceptor sites of protons and one donor site. This property supports the first compound, not only by fixing the receiver, but also activates it. It is thus about an agonist. The compound 9 presents the low coefficient of division (−3.41) and comes after the compound 10 (−3.01). These molecules are the most absorbent products. When the coefficient of division is rather low, it has as a consequence a better gastric tolerance. The compounds 6 and 8 which have, respectively, higher values 2.51 and 1.85 have capacities to be dependent on plasmatic proteins.

4. Conclusion

The present work on the oxazoles reveals that the substitution of methyl and cyanide group does not affect the heat of formation but the electronic parameters due to charge disturbance in the ring. In the substituted methyl group, 2, 4, 5-trimethyl oxazole is predicted to be the most reactive with least HOMO-LUMO energy gap of all oxazole systems.

The DFT and ab initio method can be used quite satisfactorily in predicting the chemical reactivity of the molecules and the effect of substitution of either donor electron or acceptor electron.

The compound 9 in the series of oxazole derivatives presents the lower coefficient of division ( ).

This molecule is the most absorbent product. Compound 6 has important hydration energy. It has as a consequence a better distribution in fabrics.