Novel enaminone derived from thieno [2,3-b] thiene: Synthesis, x-ray crystal structure, HOMO, LUMO, NBO analyses and biological activity

Background Due to their structural and therapeutic diversity, thienothiophene derivatives have attracted much synthetic interest because of their reactivity and biological activity. The thieno [2,3-b] thiophene moiety has been used in the design of a novel pharmaceutical therapies. Additionally, its enaminones derivatives are versatile synthons and have a lot of synthetic applications such as N-heterocycles, wide variety of naturally occurring alkaloids and pharmaceutical drugs. Results Synthesis of (2E,2′E)-1,1′-(3,4-diphenylthieno [2,3-b] thiophene-2,5-diyl) bis (3-(dimethylamino) prop-2-en-1-one) 5 was reported. The structure of compound 5 was deduced by spectroscopic techniques. The compound was crystallizes in the monoclinic system with space group P-1 with cell coordinates a=9.9685 (8) Å, b=10.1382 (8) Å, c=13.3220 (11) Å, α=101.018 (2) °, β=94.480 (2) °, γ=107.207 (1) °, V=1249.3 (1) Å3, and Z=2. In the crystal molecules are packed in chains formed via weak intermolecular C21–H21A… O1, C22–H22A…O2 and C27–H27A…O2 hydrogen bondings. Theoretical quantum chemical calculations have been performed on the studied compound using the DFT B3LYP/6-311G (d, p) method. The geometric parameters of the optimized structure are in good agreement with the experimental data obtained from our reported X-ray structure. The two benzene rings and the two side chains are not coplanar with the fused thiophene rings. The electronic spectra of the studied compound have been calculated using the TD-DFT method at the same level of theory. The transition bands at 352.9 nm (f=0.5549) and 332.1 nm (f=0.2190) are due to the H-1 → L (72%) and H → L + 1 (82%) excitations respectively. The NBO calculations were performed to predict the natural atomic charges at the different atomic sites and to study the different intramolecular charge transfer (ICT) interactions occurring in the studied system. It is found that the O and N-atoms have the highest negative charge densities while the S-atoms are the most electropositive. These results give idea about how our molecule could react with the receptor active sites. Compound 5 was evaluated against ant-microbial activity. Conclusions Synthesis, molecular structure and spectroscopic invesitgation of (2E,2′E)-1,1′-(3,4-diphenylthieno [2,3-b] thiophene-2,5-diyl) bis (3- (dimethylamino) prop-2-en-1-one) 5 was studied. Graphical Abstract Molecular structure investigation of novel enaminone derived from thieno [2,3-b] thiene. Electronic supplementary material The online version of this article (doi:10.1186/s13065-015-0100-9) contains supplementary material, which is available to authorized users.


Background
Enaminones are versatile precursors and have a lot of synthetic applications in organic chemistry. Such compounds are key synthons for the preparation of a wide variety of naturally occurring alkaloids [1,2] and nitrogen-containing heterocycles [3][4][5][6]. They have also been employed as important synthesitic intermediate for the synthesis of pharmaceutical drugs with antiviral, and larvicidal [7] anticonvulsant [8][9][10], properties. Due to their rich applications, many efficient approaches to these compounds have been developed. The structural range and biological importance of functional thienothiophenes has made them attractive research targets over the past decades.
Thienothiophenes skeletons are important in pharmaceutical research because of their versatile biological activities, such as antitumor, antiviral antibacterial, anticancer, antioxidant and α-glucuronidase and α-glucosidase inhibition, antiglaucoma activity, and inhibitors of platelet aggregation properties [11][12][13][14][15]. The thieno [2,3-b] thiophenes have been the focus of active research in recent years. For this reason we have focused to prepare certain bis-heterocycles containing thieno [2,3-b] thiophene derivatives. The molecule 5 that was prepared was found to be potent bacteria and fungus inhibitor.

Results and discussion
Chemistry Compound 5 was synthesized as depicted in Scheme 1, in 75-80% Yield. The structure was deduced by combined use of IR, 1 H-NMR, 13 C-NMR, and mass spectral data [23][24][25]. Accordingly, the assigned structure was unambiguously established via single-crystal X-ray diffraction.
Crystal structure of compound 5 Slow evaporation of glacial acetic acid solution of pure compound 5 yielded colorless crystals. A crystal of dimensions 0.54 × 0.53 × 0.30 mm was selected for X-ray diffraction analysis. Data were collected on a Bruker APEX-II diffractometer equipped with CCD detector and graphite monochromatic Mo Kα radiation (τ=71073 A) at 293 (2)°K. Cell refinement and data reduction were carried out by Bruker SAINT. SHELXS-97 [26,27] was used to solve structure ( Table 1). The final refinement was carried out by full-matrix least-squares techniques with anisotropic thermal data for nonhydrogen atoms on F2. All the hydrogen atoms were placed in calculated positions. The crystal structure 5 ( Figure 1) was finally refined with R factor of 4.43% for 4390 unique reflections. Molecules were found to be packed in crystal lattice through intermolecular hydrogen bonding (Tables 2 and 3).
H27A… O2 interaction to form chain arranged as observed in Figure 2.

Computational details
All the quantum chemical calculations of the studied compound were performed by applying DFT method with the B3LYP functional and 6-311G (d, p) basis set using Gaussian 03 software [28]. The input file was taken from the CIF obtained from our reported X-ray single crystal measurement. The geometry was optimized by minimizing the energies with respect to all the geometrical parameters without imposing any molecular symmetry constraints. GaussView4.1 [29] and Chemcraft [30] programs have been used to draw the structures of the optimized geometries. Frequency calculations at the optimized geometry confirmed that the optimized structure is an energy minimum where no imaginary frequency modes were obtained. The electronic spectra of the studied compound were calculated by the TD-DFT method. The natural bond orbital analyses were performed using the NBO calculations as implemented in the Gaussian 03 package [31] at the DFT/B3LYP level.

Optimized molecular geometry
The optimized bond lengths and bond angles obtained for the studied compound using the B3LYP method with 6-311G (d, p) basis set are given in Table 4; while the atom numbering of the optimized structure is given in Figure 3. The studied compound possesses C 1 point group. The optimized geometry of the studied compound is compared with the structural parameters obtained from the CIF data. Generally, the bond lengths and bond angles are predicted very well. Most of bond lengths are overestimated except the C=O bonds. The maximum deviations of the calculated bond length and bond angle values from the experimental data are 0.04 Å (C48-C49) and 4.3°( C48-C49-C51) respectively. These deviations are attributed to the phase difference between the calculations and Figure 1 The ORTEP diagram of the final X-ray model of compound 5 with displacement ellipsoids drawn at 30% probability level. the experiment. The calculations refer to an isolated molecule in the gas phase, while the experimental data are those for the molecule in the solid phase. Most of the deviations occur for the geometric parameters around the side chain moieties. These side chains have more freedom in the gas phase for rotation around the C-C and C-N bonds causing bigger variations from the solid state molecular structure. The calculated C-C-C bond angle values of the benzene rings are in the range of 118.8-120.6°(exp. 118.4-120.6) [32]. The two benzene rings and the two fused thiophene rings having nearly planar structure where the C-C-C-C and C-C-C-S dihedral angles ring do not exceed 0.47°and 1.16°respectively. The two benzene rings and two side chain moieties are not coplanar with the thiophene ring plane (Figure 4b). The angle between the plane passing through each benzene ring and the fused thiophene ring plane is calculated to be 70. 9.
In agreement with the X-ray structure, the DFT calculations predicted the S--O intramolecular distances to be 2.794 Å (exp. 2.799-2.804 Å). This value is smaller than the sum of the Batsanov's van der Waal radii of the two elements [33] which indicate the presence of some intramolecular S--O interactions. These intramolecular noncovalent interactions play crucial role in modulating the conformational preferences of the molecule as well as in applications of drug design [34].

Natural atomic charge
Distribution of positive and negative charges has vital role in the application of quantum chemical calculations to molecular system because of atomic charges affect dipole moment, molecular polarizability, electronic structure, acidity-basicity behavior and more lot of properties of molecular system [35]. These electronic properties have strong relations to the biological activity of compound. In this regards, the natural atomic charges (NAC) were calculated using the DFT B3LYP/6-311G (d, p) and the results were given in Table 5. The studied molecule has O, N and S heteroatoms able to interact with the receptor reactive sites. The O and N-atoms have negative natural charge densities where the O-atoms are the most electronegative atomic sites in the molecule. The natural atomic charges at the O and N-atoms are calculated to be-0.6275 and-0.4109 respectively. In contrast, the S-atoms have electropositive nature with NAC value around +0.4768. The negatively charged atomic sites (O and N) are more likely to interact with positive part of the receptor. On the contrary, the most positively charged part (S-atom) will interact quite easily with negatively charged part of the receptor. These interactions can play crucial role in bioactivity of the studied compound.
Moreover, all the C-atoms have moderate negative natural charge values except C35, C38, C48 and C51 atoms which are electropositive. The most electropositive C-sites are C48 and C51 as these atoms bonded to the highest electronegative atoms (O) in the molecule. All the H-atoms are electropositive with NAC values in the range of 0.1793-0.2135.

Frontier molecular orbitals
The energies and electron densities of the frontier molecular orbitals, HOMO and LUMO, are important electronic parameters. The latter were used to determine the most reactive sites in the unsaturated system [36]. Also, the HOMO-LUMO energy difference (ΔE) is a measure of the intramolecular charge transfer and was used in biological activity studies [37,38]. The HOMO and LUMO pictures calculated using the B3LYP/6-311G (d, p) are shown in Figure 4. The E HOMO and E LUMO values are-5.5811 and-1.5962 eV respectively. The ΔE value represents the smallest amount of energy needed for electronic excitation which belongs mainly to π-π* transition. In the studied compound, the HOMO-LUMO energy gap (ΔE) of the studied compound is 3.9849 eV. The more accurate electronic transitions have been calculated using the timedependant density functional theory (TD-DFT). The calculated electronic spectra using the TD-DFT method is shown in Figure 5. The results of the TD calculations which represent the calculated λ max values with their major contributions of molecular orbitals to the formation of bands are given in Additional file 1: Table S2.

Nonlinear optical properties
Nonlinear optical materials were used as key materials for photonic communications which use light instead of electron for data transmission. With the development of laser technology, nonlinear optical materials have been extensively applied to industry, national defense, medicine and research [39,40]. Several organic materials were used for such applications. These organic compounds were characterized by their high polarizability (α 0 ) and low HOMO-LUMO gap (ΔE). The α 0 value of the studied compound is calculated to be 409.93 Bohr 3 which is 15 times higher than the urea, the reference used for comparison of the NLO activity [32]. Moreover, the small HOMO-LUMO energy gap is an important requirement for NLO materials. The small ΔE value indicates low transition energy and hence more shift for electronic transition band to the visible region. The studied compound has lower (3.9849 eV) energy gap (ΔE) than urea (7.6644 eV). Based on these calculations, the studied compound could be predicted as better NLO material than urea [41].

Natural bond orbital analysis
The stabilization energies E (2) due to intramolecular charge transfer (ICT) interactions have been calculated using the second-order perturbation theory [42]. The significant ICT interactions are given in Table 6. The E (2) values indicate the intensity of the electron delocalization i.e. extent of conjugation of the whole system [43]. The ICT interactions formed by the orbital overlap between π → π*, n → π* and n → σ* causing stabilization of the system upto 51.63 kcal/mol. The most significant interactions in this molecule, is electron donation from LP (1) N to the antibonding π*-NBO of the adjacent C=C bonds. For example, the ICT interactions from the LP (1) N5/ LP (1) N6 to the antibonding π*(C36-C38)/C49-C51 having E (2) of 51.63 kcal/mol. Also, the π (C36-C38) → π*(O3-C35) and π (C49-C51) → π*(O4-C48) interactions resulting in stabilization of the system by 26.57 kcal/mol. These results indicate the presence of strong electron delocalizations between the C=C, C=O and the lone pair of the Natom of the side chain moiety. Interestingly, the studied compound has electron deficient bivalent sulphur atoms (~+ 0.4769) have relatively low-lying σ*-orbitals (0.209-0.235 A.U.) related to the C-S bonds that are available for interaction with the electron donor O-atoms neighboring to them. The NBO calculations predicted the energies of the σ*(C-S) NBO  The calculated dihedral angles given in Additional file 1: Table S1.
Antimicrobial activity of compound 5 Compound 5 was tested against its antimicrobial activity representing Gram-negative bacteria (Escherichia coli and Pseudomonas aeruginosa), Gram-positive bacteria (Staphylococcus pneumoniae and Bacillus subtilis), and fungi (Candida albicans and Aspergillus fumigatus), and the activities were compared with standard antimicrobial Figure 3 The optimized molecular structure and atom numbering scheme of the studied compound (a); the benzene rings and the side chains are not coplanar with the fused thiophene rings (b).

Figure 4
The ground state isodensity surface plots for the frontier molecular orbitals.
drug, specified in US pharmacopeia at 25 μg/mL. Compound 5 showed a relatively moderate inhibitory effect against Gram positive bacteria (Bacillus subtilis and Staphylococcus pneumoniae) as compared to standard drug ampicillin (Table 7).

Experimental section General
General All the chemicals were purchased from various suppliers, including Sigma-Aldrich and Fluka, and were used without further purification, unless otherwise stated. All melting points were measured on a Gallenkamp melting point apparatus in open glass capillaries and are uncorrected. IR Spectrum was recorded as KBr pellets on a Nicolet 6700 FT-IR spectrophotometer. The NMR spectra were recorded on a Varian Mercury Jeol-400 NMR spectrometer. 1 HNMR (400 MHz) and 13 C-NMR (100 MHz) were run in deuterated dimethyl sulphoxide (DMSO-d 6 ). Chemical shifts (δ) are referred to in ppm while J-coupling constants were represented in Hz. Mass spectra were recorded on a Jeol of JMS-600 H. Elemental analysis was carried out on Elmer 2400 Elemental Analyzer, CHN mode. The single-crystal X-ray diffraction measurements were performed using Bruker SMART APEX II CCD diffractometer.   Figure 5 The calculated electronic spectra of the studied compound using TD-DFT method.
(m, 5H, C 6 H 5 ); 13  Antifungal activity of compound (5) Tested sample was screened in vitro for its antifungal activity against various fungi, namely, Aspergillus fumigatus (RCMB 002568) and Candida albicans (RCMB 05036). The antifungal activity was performed by agar well diffusion method. Fungal strains were grown in 5 mL sabouraud dextrose broth (glucose/peptone; 40/10) for 3-4 days to obtain 105 CFU/mL cells. The fungal culture (0.1 mL) was spread uniformly on the sabouraud dextrose agar plates by sterilized triangular folded glass rod. Plates were left for 5-10 min so the culture is properly adsorbed on the surface. Now small wells of size 4 mm × 2 mm were cut into the plates with the help of well cutter and bottom of the wells was sealed with 0.8% soft agar to prevent the flow of test sample at the bottom of the well. 100 μL of the tested samples (10 mg/mL) was loaded into the wells of the plates. Compound 5 dissolved in DMSO, while pure DMSO was also used as control. The plates were kept for incubation at 30°C for 3-4 days and then examined for the formation of zones of inhibition. The test was performed three times for each fungus. Amphotericin B was used as standard antifungal drug.

Antibacterial activity of compound (5)
Antibacterial activities were investigated by using agar well diffusion method, against the Staphylococcus pneumonia (RCMB 010010) and Bacillus subtilis (RCMB 010067) {as Gram-positive bacteria} and Pseudomonas aeruginosa (RCMB 010043) and Escherichia coli (RCMB 0100052) {as Gram-negative bacteria}. The solution of 5 mg/mL of compound in DMSO was prepared for testing against bacteria. Centrifuged pellets of bacteria from 24 h old culture containing approximately 104-106 CFU (colony forming unit) per mL were spread on the surface of nutrient agar (type tone 1%, yeast extract 0.5%, NaCl 0.5%, agar, and 1000 mL of distilled water, pH 7.0) which was autoclaved under 121°C for at least 20 min. Wellswere created in   BD (2)  medium with the help of sterile metallic bores and then cooled down to 45°C. The activity was determined by measuring the diameter of the inhibition zone (in mm). A volume of 100 μL of the tested samples (10 mg/mL) was loaded into the wells of the plates. Solution of compound was prepared in DMSO while DMSO was also loaded as control. The plates were kept for incubation at 37°C for 24 h and then the plates were examined for the formation of zone of inhibition. Each inhibition zone was measured three times by caliper to get an average value. The test was performed three times for each bacterium. Ampicillin and was used as antibacterial standard drug.

Conclusion
The synthesis and characterization of a new 1,1′-(3,4diphenylthieno [2,3-b] thiene-2,5-diyl) bis (3-dimethylaminoprop-2-en-1-one) (5) were successfully achieved in high yield. The DFT B3LYP/6-311G (d, p) method were used to calculate the optimized molecular structure of the studied compound. The optimized molecular structure showed good agreement with our reported X-ray crystal structure. The natural atomic charges were calculated using the NBO method. The O and S atoms are the most probable sites to react with the active sites of the receptor molecule. The α 0 and HOMO-LUMO energy gap (ΔE) values indicated that the studied molecule is considered as a better NLO material than urea. The calculated electronic spectra using TD-DFT method showed two intense transition bands at 352.9 nm (f = 0.5549) and 332.1 nm (f = 0.2190) due to the H-1 → L (72%) and H → L + 1 (82%) transitions respectively. There is strong electron delocalization from the Lp (1) N to the adjacent C=C and extended to the C=O group.

Additional file
Additional file 1: Additional information.