PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of springeropenLink to Publisher's site
Journal of Molecular Modeling
 
J Mol Model. 2017; 23(4): 107.
Published online 2017 March 8. doi:  10.1007/s00894-017-3284-1
PMCID: PMC5343086

Esters with imidazo [1,5-c] quinazoline-3,5-dione ring spectral characterization and quantum-mechanical modeling

Abstract

1-phenyl-2H,6H-imidazo[1,5-c]quinazoline-3,5-dione reacts with ethyl bromoacetate under mild conditions to give 2-(ethoxycarbonylmethyl)-1-phenyl-6H-imidazo[1,5-c]quinazoline-3,5-dione (MEPIQ) and next 2,6-bis(ethoxycarbonylmethyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione (BEPIQ). The products were isolated at high yield and identified on the basis of IR, 1H- and 13C-NMR, UV spectroscopy, and X-ray crystallography. Diester (BEPIQ) can be presented by 16 possible pair of enantiomers. Only one pair of them is the most stable and crystallizes which is shown crystallographic research. Based on quantum-mechanical modeling, with the use of DFT method, which conformers of mono- and diester and why they were formed was explained. It was calculated that 99.93% of the monoester (MEPIQ) is formed at position No. 2 and one pair of the monoester conformers, from six possible, has the largest share (51.63%). These results afforded to limit the number of diester conformers to eight. Unfortunately, the quantum-mechanical calculations performed that their shares are similar. Further quantum-mechanical modeling showed that conformers are able to undergo mutual transformations. As a result only one pair of diester conformers forms crystals. These conformers have substituents in trans position and these substituents are located parallel to imidazoquinazoline ring. This allows for the denser packing of the molecules in the unit cell.

Electronic supplementary material

The online version of this article (doi:10.1007/s00894-017-3284-1) contains supplementary material, which is available to authorized users.

Keywords: Ester, Imidazo[1,5-c]quinazoline ring, Quantum-mechanical modeling, Spectral characterization

Introduction

Compounds containing imidazoquinoline and imidazoquinazoline moieties have been reported to possess interesting properties. Until now, more than 100 imidazo[4,5-c]quinolin-2-one derivatives have been described in the literature, of which at least half show biological activity in various aspects. More recently, due to their biological properties, there has also been great interest in 3,3-disubstituted quinoline-2,4-diones such as, for example, 3,3-diazidoquinoline-2,4-dione which behaes as an inhibitor of the platelet aggregation and 3-hydroxy-3-alkylquinoline-2,4-dion which occurs in the bacteria organisms and it works as an antibiotic [1].

Imidazoquinazolines have been widely applied in pharmacy and medicine due to their versatile biological activities, such as antitumor, antiviral, antibacterial, and anticonvulsant activity. There are drugs containing imidazoquinazolines currently available on the market [2]. One of them is anagrelide, i.e., 6,7-dichloro-1H,5H-imidazo[2,1-b]quinazoline-2(3H)-one, used in the treatment of thrombosis, essential thrombocythemia, or the chronic myelogenous leukemia. Another example is quazinone, i.e., (3R)-6-chloro-3-methyl-5H,10H-imidazo[2,1-b]quinazoline-2(3H)-one, which is widely distributed as a component of heart disease drugs. These types of substances are also derivatives of imidazo[1,5-a]quinazoline such as NNC 14-0185 and 14-0189 NNC, which exhibit an anticonvulsant activity in rats and mice. In the future, they may be used for the production of antiepileptic drugs [2, 3].

Research has shown that compounds containing imidazoquinazoline moiety can be used as the neutralizing agents’ free radicals and thus prevent lipid peroxidation and cell damage. Tested imidazo[1,2-c]quinazoline derivatives were found to be good inhibitors of the lipid peroxidation [4].

Imidazo[1,5-c]quinazoline-3,5-dione derivatives, obtained during reaction 3-aminoquinolinediones with urea as a result of rearrangement [5, 6], can also have biological activity. Unfortunately, they are sparingly soluble and because of that fact their application is difficult. Therefore, there is a need to find more soluble derivatives of these compounds. One of the ways is to obtain their ester derivatives.

This paper illustrates the example of 1-phenyl-2H,6H-imidazo[1,5-c]quinazoline-3,5-dione (PIQ). With the use of this compound, mono- and diester were obtained and their spectral and physical properties characterization was made. In addition, the spatial structure of the obtained esters was explained based on the quantum-mechanical calculations by means of DFT method.

Methods

Materials

PIQ was prepared according to literature procedure [7]. The rest of the reagents were purchased and used as received: potassium carbonate, pure for analysis, POCH, Poland; N,N-dimethylformamide (DMF), pure for analysis, Chempur, Poland; ethyl bromoacetate, 98%, Sigma-Aldrich, US; chloroform, pure for analysis, Chempur, Poland; sodium sulfate, anhydrous, ≥ 99.0%, ACS Reagent, US; benzene, pure for analysis, Chempur, Poland.

Synthetic procedures

An external file that holds a picture, illustration, etc.
Object name is 894_2017_3284_Figa_HTML.jpg

A mixture of PIQ (2.77 g, 10 mmol) and potassium carbonate (3.04 g, 22 mmol) in DMF (40 cm3) was stirred at room temperature. After 10 min, ethyl bromoacetate (1.33 cm3, 12 mmol or 2.66 cm3, 24 mmol) was added. Reaction mixture was protected from moisture with a tube filled with potassium hydroxide and was stirred for 4 h at room temperature and 1 h at temperature 40 °C. After cooling, the reaction mixture was extracted with chloroform. The extract was dried sodium sulfate. After evaporation of chloroform, the crude product was recrystallized from benzene. The course of reaction and also the purity of substances were monitored by TLC (elution systems chloroform-ethanol, 9:1) on Alugram SIL G/UV254 foils (Macherey-Nagel).

2-(ethoxycarbonylmethyl)-1-phenyl-6H-imidazo[1,5-c]quinazoline-3,5-dione (MEPIQ); Yield 70%, light-yellow crystals, mp. 174 °C (benzene), IR: 3570-3320 (s, N−H valence), 3066 (w, CH, deformation of phenyl ring), 2982 (w, CH3- and - CH2 -, asymmetric, valence), 2941 (w, CH3- and - CH2 -, symmetric, valence), 1762 (s, C = O, valence), 1683, 1613, 1506, 1443 (s, skeletal of phenyl ring), 1639, 1587, 1481 (s, skeletal of quinazoline ring), 1201 (s, C−O(O), symmetric, valence), 1029 (m, C−O(O), asymmetric, valence), 752, 700 (w, C-H in Ph ring, nonplanar deformation), [cm−1]; APCI-MS: m/z 365 [M + H]+ (100%). APCIMS/MS of precursor ion m/z 365: m/z 337 [M + H-CO]+, 319 [M + H-CO-H2O]+, 294 [M + H–NHCO–CO] + (100%), 263, 217; EA: Anal. calcd (found) for C20H17N3O4: C 69.31 (69.16); H 4.72 (4.69); N 11.56 (11.51); 1H-NMR (500 MHz, d6-DMSO), δ = 1.10 (3H, t, CH3, J22,23 = 6.26 Hz), 4.05 (2 H, q, CH2 , J 22,23 = 6.26 Hz), 4.22 (2 H, s, CH2), 6.72 (1 H, d, J10,11 = 8.05 Hz), 6.78 (1 H, t, J 10,9 = 7.56 Hz), 7.02 (1 H, d, J 9,8 = 8.05 Hz), 7.14 (1 H, t, J 11,10 = 7.56 Hz), 7.46 (2 H, m), 7.62 (3 H, m), 10.80 (1 H, s); 13C-NMR (d6-DMSO), δ = 167.67 (C21), 147.85 (C3), 144.99 (C5), 134.54 (C13), 130.76 (C15), 130.39 (C19), 129.67 (C17), 128.48 (C10), 129.33 (C14), 127.53 (C16 and C18), 123.18 (C9), 121.53 (C11), 118.05 (C1), 114.86 (C8), 113.32 (C7), 111.02 (C12), 62.33 (C22), 43.40 (C20), 13.92 (C23), [ppm]; UV: 206, 261, 332, [nm].

2,6-bis(ethoxycarbonylmethyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione (BEPIQ); Yield 85%, colorless crystals, mp. 182 °C (benzene), IR: 3057 (w, CH, deformation of phenyl ring), 2983 (w, CH3- and - CH2 -, asymmetric, valence), 2941 (w, CH3- and - CH2 -, symmetric, valence), 1745 (s, C = O, valence), 1681, 1604, 1505, 1444 (s, skeletal of phenyl ring), 1641, 1583, 1486 (s, skeletal of quinazoline ring), 1211 (s, C−O(O), symmetric, valence), 1026 (m, C−O(O), asymmetric, valence), 749, 664 (w, C-H in Ph ring, nonplanar deformation), [cm−1]; APCI-MS: m/z 450 [M + H]+ (100%); APCIMS/MS of precursor ion m/z 450: m/z 422 [M + H–CO]+, 404 [M + H–CO-H2O]+, 376 [M + H–NHCO–CO] + (100%), 348, 302; EA: Anal. calcd (found) for C24H23N3O6: C 64.11 (64.09); H 5.16 (5.23); N 9.34 (9.14); 1H-NMR (500 MHz, d6-DMSO), δ = 1.09 (3H, t, CH3, J22,23 = 6.99 Hz), 1.22 (3H, t, CH3, J26,27 = 6.99 Hz), 4.05 (2 H, q, CH2 , J 22,23 = 6.99 Hz), 4.19 (2H, q, CH2, J 26,27 = 6.99 Hz), 4.24 (2 H, s, CH2), 4.89 (2 H, s, CH2), 6.82 (1 H, d, J10,11 = 7.59 Hz), 6.88 (1 H, t, J 10,9 = 7.59 Hz), 7.12 (1 H, d, J 9,8 = 8.39 Hz), 7.24 (1 H, t, J 11,10 = 7.79 Hz), 7.48 (2 H, m), 7.61 (3 H, m), [ppm]; 13C-NMR (d6-DMSO), δ = 168.20 (C25), 167.67 (C21), 147.80 (C3), 144.99 (C5), 134.54 (C13), 130.68 (C15), 130.34 (C19), 129.63 (C17), 128.59 (C10), 128.22 (C14), 127.18 (C16 and C18), 123.16 (C9), 121.48 (C11), 117.63 (C1), 114.90 (C8), 113.43 (C7), 112.08 (C12), 61.21 (C26), 61.15 (C22), 43.97 (C24), 42.26 (C20), 13.96 (C23), 13.80 (C27), [ppm]; UV: 207, 263, 332, [nm].

Methods

IR spectra

Infrared spectra (4000−400 cm−1), as obtained from KBr disks, were recorded on a Bruker ALPHA FT-IR instrument, with a resolution of 0.01 cm−1.

NMR spectra

NMR spectra were recorded using Bruker 500 MHz spectrometer in deuterated dimethyl sulfoxide (DMSO-d6). 1H and 13C chemical shifts were given on the δ scale (ppm) and were referenced to internal tetramethylsilane (TMS). All two-dimensional (2D) experiments for correlation spectroscopy (1H,1H-COSY), heteronuclear single-quantum correlation spectroscopy (HSQC), heteronuclear multiple-bond correlation spectroscopy (HMBC) were performed using manufacturer’s software. Proton spectra were assigned using COSY. Protonated carbons were assigned by HSQC. Quaternary carbons were assigned by HMBC.

Mass spectra

The positive-ion APCI mass spectra were measured on an ion trap analyzer Esquire 3000 (Bruker Daltonics, Bremen, Germany) within the mass range m/z = 50-1000. Samples were dissolved in chloroform and analyzed by direct infusion at the flow rate of 40 mL/min. The ion source temperature was 300 °C, the APCI probe temperature was 350 °C, the flow rate and the pressure of nitrogen were 3 L/min and 25 psi, respectively. For MS/MS measurements, the collision amplitude was 0.9 V and the isolation width of precursor ions was 4 m/z.

UV spectra

Ultraviolet-visible spectroscopy (UV-VIS) measurements were carried out at 25 °C in the range of 200-700 nm with a spectrometer Hewlett Packard 8943 in the methanol solution of a measuring cell with a thickness of 1 cm.

Elemental analysis

Elemental analyses (C, H, N) of products were carried out on an elemental analyzer Vario ELIII instrument. It should be noted that all kinds of other atoms were not notified of the test samples even in the residual contents.

X-ray crystallography

Single-crystal X-ray diffraction data were collected on a Agilent Technologies SuperNova single source diffractometer with MoKα radiation (λ = 0.71073 Å) at 100 (2) K using CrysAlis RED software [8]. The multi-scan empirical absorption correction using spherical harmonics was applied as implemented in SCALE3 ABSPACK scaling algorithms [8]. The structural determination procedure was carried out using the SHELX package [9]. The structure was solved with direct methods and then successive least-square refinements were carried out based on the full-matrix least-squares on F 2 using the XLMP program [9]. All H-atoms bound to C-atoms were positioned geometrically, with C–H equal to 0.93 Å, 0.96 Å, and 0.97 Å for the aromatic, methyl, and methylene H-atoms, respectively, and constrained to ride on their parent atoms with U iso (H) = xU eq (C), where x = 1.2 for the aromatic and methylene H-atoms, and x = 1.5 for the methyl H-atoms. The figures for this publication were prepared using ORTEP-3 and Olex2 programs [10, 11].

Quantum mechanical calculations

The quantum-mechanical calculations were performed using a set of computational methods based on the electron density of the tested system in the stationary state, i.e., density functional theory (DFT) [12]. The main objective of DFT is to find a value of the functional, which requires an application of successive approximations. Among the proven and widely used functionals, the B3LYP functional was selected (Becke 3-term correlation functional; Lee, Yang, and Parr exchange functional) [13, 14] as the most suitable for the tasks associated with organic compounds. Every functional requires a functional sets, which allows the replacement of differential equations-integral by a system of algebraic equations. In the preliminary studies, the following function bases were taken into account: 6-31G(d,p) [15] 6-311++G(d,p) [16] and aug -CC-pVDZ [17, 18]. They are the most often used function bases in an optimization of the spatial structure of compounds which are structurally similar. Quantum-mechanical calculations were performed by using Gaussian program (version 9) [19].

It should also be mentioned that the computing cluster called HYDRA (model Blade System Actina Solar Hewlett-Packard) working under the control of the operating system Scientific Linux 6 was used for the calculations. The calculated values of Gibbs free energy and other necessary parameters were read using Notepad++ and EDA-Reader applications [20, 21]. Visualizations of the esters structures were performed with the use of Gauss View and Mercury software [22, 23].

Results and discussion

Synthesis and spectral characterization of esters with imidazo[1,5-c]quinazoline ring

In the course of our work, we have focused on the subject of the synthesis and characterization of the esters with imidazo[1,5-c]quinazoline ring, i.e., 2-(ethoxycarbonylmethyl)-1-phenyl-6H-imidazo[1,5-c]quinazoline-3,5-dione (MEPIQ) and 2,6-bis(ethoxycarbonylmethyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione (BEPIQ).

The products were obtained by simple reaction of PIQ with an equimolar amount or 2-molar excess of ethyl bromoacetate in the presence of a potassium carbonate catalyst. They have been isolated in pure form in good yield 70% (MEPIQ) and 85% (BEPIQ), respectively.

The solubility of both esters is much better than PIQ. These compounds have good solubility (besides DMSO and DMF as PIQ) in benzene, chloroform, acetone, ethyl acetate, and ethyl alcohol.

MEPIQ is formed as the only reaction product of equimolar amounts of the PIQ with ethyl bromoacetate. The reasons for the reaction chemoselectivity during the obtaining of monosubstituted derivatives of PIQ were explained in detail earlier in paper [24].

The new esters have been identified by elemental analysis, the 1H- and 13C-NMR (Fig. 1S, 2S, 4S and 5S), and IR (3S and 8S) and APCIMS/MS measurements (see Experimental part).

All the signals in 1H- and 13C-NMR spectra of BEPIQ were assigned on the basis of COSY, HSQC, and HMBC experiments to the corresponding atoms. Signals in 1H- and 13C-NMR spectra of MEPIQ were assigned to the proper atoms based on 1H- and 13C-NMR spectra of BEPIQ. The signal positions are in accord with the proposed structure MEPIQ and BEPIQ.

A comparison of the elemental analysis results of MEPIQ or BEPIQ and the calculated amount of corresponding elements confirmed the esters compositions.

In the IR spectra of the new esters (Fig. 3S and 8S), the significant changes are observed in the range of 2000-3500 cm−1. The band of the valence vibrations of N-H bonds disappeared partially or totally in the case of MEPIQ and BEPIQ, respectively. The band of the valence vibrations C-H bonds of the aromatic ring is revealed at 3030 cm−1. Furthermore, there are visible bands of the symmetrical and asymmetrical valence vibrations of C-H bonds of the methylene and methyl groups at 2890 and 2960 cm−1. In the IR spectra (Fig. 3S and 8S), the ν(C=O) stretching vibration bands of MEPIQ (1762 cm−1) and BEPIQ (1745 cm−1) shift slightly upon substitution of PIQ (1738 cm−1) similar to all skeletal bands of imidazoquinazoline and phenyl rings. The ν(C-O) stretching asymmetrical and symmetrical vibration bands of ester (CO)-O bond appear at 1201 (1211) and 1029 (1026) cm−1 in the spectrum of MEPIQ (BEPIQ), indicating that substitution at nitrogen atoms takes place.

In the 13C-NMR spectrum of MEPIQ (Fig. 2S), the three signals of the sp3 hybridized carbon atoms are found at 62.33 (C22), 43.40 (C20), and 13.92 (C23) ppm. In the 13C-NMR spectrum of BEPIQ (Fig. 5S), the three additional signals of the sp3 carbon atoms are present. Thus, we observe six signals at 61.21 (C26), 61.15 (C22), 43.97 (C24), 42.26 (C20), 13.96 (C23), and 13.80 (C27) ppm. The occurrence in the 13C NMR spectra, shift of the carbonyl group at 167.67 (C21) (MEPIQ) or 168.20 (C25) and 167.67 (C21) (BEPIQ), is typical for the ester groups.

Confirmation of the ester structures was obtained by a complete interpretation of NMR spectra. NH protons of PIQ resonating at ca. 11.15 ppm disappear in the 1H-NMR spectrum of MEPIQ (Fig. 1S) compared to PIQ spectrum. It confirms the substitution at nitrogen atom No. 2. Proton signals of ethyl group appear as a triplet at 1.10 ppm and a quartet 4.05 ppm. Additionally, a singlet of methylene protons linked to nitrogen atom No. 2. lies at 4.22 ppm.

Proton signals at carbon atoms No. 8, 9, 10, and 11 in the quinazoline ring are clearly separated compared to 1H-NMR spectrum of PIQ and appear at 7.12, 6.88, 7.25, and 6.83 ppm, respectively. Similarly, the phenyl ring protons give two distinct signals at 7.46 and 7.62 ppm (see Experimental part).

The 1H-NMR spectrum of BEPIQ (Fig. 4S) is quite similar to that of MEPIQ. The signals of both esters differ only slightly, both in the proton and carbon NMR spectra, and have a relative integral ratio of proton signals corresponding to calculated ones from the esters structures. Nevertheless, in 1H-NMR spectrum of BEPIQ, proton signal NH at 10.80 ppm is not observed. Moreover, three new signals appear: a triplet at 1.22 ppm (CH3), a quartet at 4.19 ppm (O-CH2), and a singlet at 4.89 ppm (N-CH2), from the ester group substituted at nitrogen atom No. 6.

The positive-ion atmospheric pressure chemical ionization (APCI) mass spectra of both esters studied show the protonated molecules [M + H]+ as the only ions in the spectra, which confirms unambiguously the expected molecular weights. The fragment ions observed in MS/MS spectra are in accordance with the suggested structures. The primary cleavage usually leading to the base peak of MS/MS spectrum is the neutral loss of the side chain, as the losses of CO, H2O, and others.

Crystallographic characterization of 2,6-bis(ethoxycarbonylmethyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione

Single crystal X-ray diffraction data collection was performed for the BEPIQ compound, crystals of which were obtained by recrystallization of a crude product from benzene. It appears that BEPIQ crystallizes in the orthorhombic P212121 space group in the form of two enantiomers despite of the lack of stereogenic center. The enantiomers are formed in statistical 50%:50% proportions and they are not distinguishable by NMR technique.

Asymmetric unit of the crystal lattice except of two molecules of the investigated compound contains also a pair of solvent molecules (benzene) (Fig. 1).

Fig. 1
Molecular structure of benzene-solvated BEPIQ with atomic labels. Displacement ellipsoids are drawn at 25% probability level and H-atoms are shown as small spheres of arbitrary radius. The C–H · · · O ...

Table Table11 below contains crystallographic data. Selected bond lengths are presented in the caption of Fig. 1 (full list of geometric parameters, including bond lengths, valence and torsion angles can be found in the Supporting information—see Table 1S).

Table 1
Crystallographic data and structural refinement details of benzene-solvated BEPIQ

The analysis of molecular overlay (Fig. 2) revealed that the geometry of benzene molecules denoted as A and B is essentially the same. Both BEPIQ molecules have substituents in the trans position and the ethoxycarbonylmethylene groups located parallel to the imidazoquinazoline ring. However, some visible differences in their geometry can be noticed. The respective values of the average deviations from planarity of imidazo[1,5-c]quinazoline-3,5-dione ring system labeled as A and B are 0.089 Å and 0.029 Å, respectively. The phenyl substituent in BEPIQ molecules A and B is inclined to the aforementioned ring system at a dihedral angle of 64.7° and 65.8°. More evident differences occur in orientations of the ethoxycarbonylmethylene substituents. Those changes can be illustrated by specifying the values of dihedral angles between the mean-square quadratic planes defined by the atoms of carboxylic groups and the imidazo [1,5-c] quinazoline-3,5-dione moiety. Thus, the carboxylic groups C21A/O22A/O25A and C21B/O22B/25B are inclined to the specified ring system, respectively, by the angle of 77.2 and 85.7°. The carboxylic groups C29A/O30A/O33A and C29B/O30B/O33B are oriented in such a way that the dihedral angle relative to the imidazo[1,5-c]quinazoline-3,5-dione fragment equals 79.3 and 72.9°.

Fig. 2
The alignment of benzene-solvated BEPIQ molecules labeled as A and B. Alignment RMSD {N2B, C1B, C13B, …} to {N2A, C1A, C13A, …} with inversion 0.382 Å, while the alignment RMSD {C39B, C38B, C37B, …} to {C39B, C38B, ...

As mentioned previously, in the case of both BEPIQ conformers, the ester substituents are located parallel to the imidazoquinazoline ring. It allows for a denser packing of the adjacent molecules in the crystal lattice. The search for intermolecular interactions using PLATON [25] has shown that adjacent molecules of BEPIQ are linked by the network of the C–H · · · O, C = O · · · π and π-π intermolecular interactions, which causes the creation of supramolecular framework surrounding 1D-channels filled with the solvent (Fig. 3). The benzene molecules denoted as A are interacting by some weak C–H · · · O intermolecular interactions engaging the O26B carbonyl oxygen atoms of BEPIQ molecules denoted as B (Fig. 1).

Fig. 3
The arrangement of molecules in the crystal lattice of benzene-solvated BEPIQ viewed along a a and b-direction b. Hydrogen atoms not involved in H-bond type interactions have been omitted for clarity. 1D-channels filled with solvent molecules have been ...

Quantum-mechanical modeling of esters with imidazo[1,5-c]quinazoline-3,5-dione ring

The BEPIQ monocrystalline study has revealed that only one pair of its conformers (Fig. 1) exists, while theoretically there are 32 possible enantiomers (16 pairs) (Fig. 10S). In order to explain this phenomenon, conformers of monoester were analyzed, at the beginning.

Spatial structures of possible conformers of PIQ monoesters substituted at nitrogen atom No. 2 or No. 6 are shown schematically in Fig. 4. The imidazoquinazoline ring is a flat molecule, whereas the phenyl ring is set at a determined angle toward the plane of the imidazoquinazoline ring. The phenyl ring can take two boundary positions. In turn, the ester substituent (ethoxycarbonylmethylene group) for a given substitution reaction can be located above or below the imidazoquinazoline ring plane. In each of these positions, it is located either perpendicular or parallel to the mentioned ring. Therefore, each monoesters substituted at the nitrogen atom No. 2 or No. 6 can theoretically exists in the form of eight conformers (Fig. 11S).

Fig. 4
Possible conformers of monoesters PIQ

Quantum-mechanical calculations revealed that there are eight stable conformers of monoester substituted at the nitrogen atom No. 6 and only six stable conformers substituted at the nitrogen atom No. 2. Quantum-mechanical calculations were carried out based on the total energy of the conformer as an optimization criterion and on the assumption that there is no interaction between the conformers.

Next, the percentage share of all 14 conformers of MEPIQ was determined taking into account their value of Gibbs free energy (enthalpy) (Table 5S). The results, shown in Table Table2,2, indicate that conformers substituted at nitrogen atom No. 6 almost do not exist. It has been shown that in 99.93% of cases, the formation of the monoester occurs at position No. 2, which was confirmed experimentally.

Table 2
The percentage of N-substituted monoesters of PIQ

Taking into account the most probable conformers of the MEPIQ obtained by substituting at nitrogen atom No. 2, the number of diester conformers decreases to 24 (Fig. 11S).

Among the six conformers of 1-phenyl-2-(ethoxycarbonylmethyl)-2H-imidazo[1,5-c]quinazoline-3,5-dione the one pair of enantiomers (E’ and F’) has the largest share (51.63%). These enantiomers have parallel situated ester substituent to the imidazoquinazoline ring. The rest are the enantiomers with the ethoxycarbonylmethylene substituent located perpendicular to the imidazoquinazoline ring. A parallel arrangement of the substituent is responsible for the spatial discharge (as in the case of both boat and chair conformations of substituted cyclohexane derivatives [26]) and lower energy of conformers, and thus a larger participation in the total population.

Figure 5 shows the four proposed conformers for monoester substituted at the nitrogen atom No. 2 (G’, F’, H’, E’). These prevailing conformers (E’ and F’) are shown in full colors.

Fig. 5
Suggested conformers monoester MEPIQ substituted at the nitrogen atom No. 2 (G’, F’, H’, E’)

In both conformer E’ and F’, the plane of the phenyl ring is parallel to the plane of the C-C bond (a methylene group–a carbonyl carbon). Therefore, the distance between the carbonyl group and the phenyl ring is the shortest, and the distance between the oxygen atom of the carbonyl group and the nearest hydrogen atom of the phenyl ring is observed to be 2.46 Å [27, 28]. This distance enables an electrostatic interaction between the phenyl ring and the carbonyl group — a high electron density on the oxygen atom of the carbonyl group and a low electron density at the edge of the phenyl ring (Fig. 6). Thus, the conformer becomes more stable. In turn, both longer distances (~4 Å) and weaker interactions are observed for the G’ and H’ conformers.

Fig. 6
Electrostatic potential of conformer F’

Therefore, due to the electrostatic interactions, the ring is “pulled” toward the carbonyl group by the charge existing on the oxygen atom, and the more stable conformer is formed.

This electrostatic interaction is responsible for the conformational transition of the H’ conformer into the E’ and the G’ conformer into the F’ as a result of changing the angle of inclination of the phenyl ring plane (Fig. 7). This explains a reduction in the number of stable conformers of 1-phenyl-2-(ethoxycarbonylmethyl)-2H-imidazo [1,5-c] quinazoline-3,5-dione from eight to six, mentioned earlier.

Fig. 7
The comp image scheme of transformation conformer H’ into conformer E’

Taking into account the predominant conformers E’ and F’ of MEPIQ, four pairs of BEPIQ conformers were determined (Fig. 8). Therefore, eight conformers of diesters have been subjected to quantum-mechanical modeling in order to clarify the reason for the formation of only one conformer in the form of two enantiomers G-E’ and H-F’ (Fig. 1).

Fig. 8
Possible conformers of 1-phenyl-2,6-bis(ethoxycarbonylmehtyl)imidazo[1,5-c]quinazoline-3,5-dione, drawn schematically

As a result of the calculation, the values of Gibbs free energy were provided and then used to determine the shares of the various conformers, i.e., four pairs of enantiomers in the total population. The shares are as shown in Table Table2.2. The content of the enantiomer pair (H-F’ and G-E’) found by crystallography method equals 21.09 mol% (Table (Table33).

Table 3
The percentage of 2,6-bis(ethoxycarbonylmehtyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione conformers

In further investigations, the possibility of mutual transformations of diester conformers was considered. Therefore, the simulation of the total energy value changes of the selected conformers was performed. At first, an analysis of the transformation of D-F’ into H-F’ (one ester substituent is below and the second one is above the imidazoquinazoline ring plane, i.e., they are in trans position) by rotation around the C-C bond (bond between a carbonyl group–a methylene group) was performed (Fig. 9). As a result of the rotation, the substituent which was perpendicular to imidazoquinazoline ring plane becomes parallel and now both substituents of H-F’ take parallel position. It was noted that the energy barrier of the free rotation depends on the rotation direction. A rotation out of the plane of the imidazoquinazoline ring involves energy about 10 kJ mol−1. In turn, a rotation above the plane of the imidazoquinazoline ring requires energy about 18 kJ mol−1. These values indicate that free rotation is slightly limited but it is possible. For comparison, the difference in total energy values of the terminal ethane conformers’ amounts to 12 kJ mol−1 and a rate of the conformation change is 106/s at room temperature [29].

Fig. 9
Energy torsional strain in function of the rotation angle about the C-C bond (carbonyl group–methylene group) in the molecule in 2,6-bis(ethoxycarbonyomehtyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione

As shown in Fig. 10, an analogous situation refers to transformations of the A-F’ isomer into the F-F’ (both substituent lie at the same side of the imidazoquinazoline ring plane; they are in cis position). The ester substituent also changes its position from the perpendicular to the parallel.

Fig. 10
Energy torsional strain in function of the rotation angle about the C-C bond (carbonyl group - methylene group) in the molecule in 2,6-bis(ethoxycarbonyomehtyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione

Furthermore, the transformation of conformers B-E’ (cis conformer) into G-E’ (trans conformer), when the configuration change takes palce was considered. Then, the ester substituent changes its location form one side of the imidazoquinazoline ring plane to the second one, and its position changes from the parallel to the perpendicular. It was calculated that the energetic barrier of configuration change (cis into trans) B-E’ in G-E’, amounts to 29 kJ mol−1 (or 53 kJ mol−1 in the case of conversion from the steric hindrance). It should be noted that the energetic barrier of cyclohexane ring inversion is 46.2 kJ mol−1 at the rate of 105/s at room temperature [30].

The obtained resuts indicate the possibility of spontaneous reconfiguration of all diester conformers.

Therefore, it can be thought that during the process of crystal forming, a free transformation occurs according to the scheme shown in Fig. 11.

Fig. 11
Diagram of the possible conformer transformations

Finally, only the H-F’ conformer and its mirror image (G-E’) form crystals, as is demonstrated by crystallographic studies (Fig. 1). These conformers have substituents in trans positions and the substituents lie parallel to imidazoquinazoline ring plane. It allows for the denser packing of the molecules in the unit cell. As X-ray diffraction analysis shows (Fig. 12), there are four pairs of enantiomers in the unit cell and, moreover, their imidazoquinazoline rings are parallel to each other. The detailed explanation of the diester conformer packing will be the subject of the next paper.

Fig. 12
Unit cell of crystallographic network of 1-phenyl-2,6-bis(ethoxycarbonylmehtyl)imidazo[1,5-c]quinazoline-3,5-dione

Crystallographic structure versus quantum-mechanical model

Comparing the bond lengths from the energetically optimized model from Gaussian with the values measured in crystal (Table 24S), one can see that differences are mostly not bigger than 2%. Merely exceptions are lengths of all bonds between hydrogen and carbon atoms (values of differences amounts from 12 to 17%) and lengths of two bonds hydrogen-oxygen atoms in one of the enantiomers, where deviations reach 11%.

So significant values of differences could be caused by the kind of basis set we used: 6-31++G(d,p). Changing the accuracy of the used data base (e.g., 6-31G(d,p) or 6-311++G(d,p)), did not improve the values.

Differences in angles between bonds (Table 25S) do not exceed 3%, however, maximum value of difference is 3° and it is 111° for O22A-C23A-C24A angle. It is related directly to differences for dihedral angles (Table 26S) of ethoxycarbonylmethyl groups, where relatively large deviations are observed in: C24A—C23A—O22A—C21A, C32A—C31A—O30A—C29A, C32B—C31B—O30B—C29B, and C24B—C23B—O22B—C21B. As we can see, those deviations are related to positions of ethoxycarbonylmethyl groups in crystals relative to simulated model, and probably caused by getting stuck in local minimum by optimization algorithm.

Conclusions

This work provides the synthesis and spectral characterization of two new compounds — mono- and diesters with the imidazoquinazoline ring.

Diester structure was investigated by single-crystal X-ray diffraction. Crystallographic studies have shown that only one pair of 2,6-bis(ethoxycarbonylmehtyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione conformers exists among 16 theoretically possible pairs.

Quantum-mechanical modeling of mono- and diesters explains why only one pair of diester enantiomers is formed in a crystalline structure. The modeling results clearly indicate that the diester conformers are capable of mutual transformations. The energy barriers of the conformer changes are not very high and afford the conformer transformations. Therefore, finally, only one pair of 2,6-bis(ethoxycarbonylmehtyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione enantiomers forms crystals. The ester substituents of these conformers have beneficial parallel position toward the imidazoquinazoline ring plane. The imidazoquinazoline rings of these conformers are located parallel to each other in the cell unit for the improved packaging of molecules.

Electronic supplementary material

Below is the link to the electronic supplementary material.

ESM 1(4.9M, doc)

(DOC 5084 kb)

Acknowledgements

NMR spectra was made in the Laboratory of Spectrometry, Faculty of Chemistry, Rzeszow University of Technology and was financed from DS budget.

The authors would like to thank Interdisciplinary Centre for Mathematical and Computational Modelling in Warsaw for providing computer facilities for possibility of performing of the quantum-mechanical calculations under grant G49-12.

Antonín Klásek thanks the financial support from the internal grant of TBU in Zlín (No. IGA/FT/20017/005), funded from the resources of specific university research.

Authors thank Hana Gerzova (Department of Chemistry, Faculty of Technology, Tomas Bata University in Zlin, CZ-762 72 Zlin, Czech Republic) for help in synthesis of 2,6-bis(ethoxycarbonylmethyl)-1-phenylimidazo[1,5-c]quinazoline-3,5-dione.

Footnotes

This paper belongs to Topical Collection 7th Conference on Modeling & Design of Molecular Materials in Trzebnica (MDMM 2016)

CCDC-1499889 entries contain the supplementary crystallographic data (cif file) for this paper. These data can be obtained free of charge from The Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/data_request/cif.

References

1. Stadlbauer W, Laschober R, Lutschounig H, Schindler G, Kappe T. Halogenation reactions in position -3 of quinoline-2,4-dione systems by electrophilic substitution and halogen exchange. Monatsh Chem. 1992;123:617–636. doi: 10.1007/BF00816857. [Cross Ref]
2. Li G, Kakarla R, Gerritz SW, Pendri A, Mac B. A facile one-step synthesis of 5-chloro-imidazo[1,5-a]quinazoline by microwave irradiation. Tetrahedron Lett. 2009;50:6048–6052. doi: 10.1016/j.tetlet.2009.08.054. [Cross Ref]
3. Chen Z, Huang X, Yang HY, Ding WB, Gao LY, Yea ZQ, Zhang YD, Yu YP, Lou YJ. Anti-tumor effects of B-2, a novel 2,3-disubstituted 8-arylamino-3H-imidazo[4,5-g]quinazoline derivative, on the human lung adenocarcinoma A549 cell line in vitro and in vivo. Chem Biol Interact. 2011;189:90–99. doi: 10.1016/j.cbi.2010.11.004. [PubMed] [Cross Ref]
4. Domány G, Gizur T, Gere A, Takács-Novák K, Farsang G, Ferenczy GG, Tárkányi G, Demeter M. Imidazo[1,2-c]quinazolines with lipid peroxidation inhibitory effect. Eur J Med Chem. 1998;33:181–187. doi: 10.1016/S0223-5234(98)80007-X. [Cross Ref]
5. Mrkvicka V, Lycka A, Rudolf O, Klásek A. Reaction of 3-aminoquinoline-2,4-diones with isothiocyanic acid - an easy pathway to thioxo derivatives of imidazo[1,5-c]quinazolin-5-ones and imidazo[4,5-c]quinolin-4-ones. Tetrahedron. 2010;66:8441–8445. doi: 10.1016/j.tet.2010.08.056. [Cross Ref]
6. Mrkvicka V, Rudolf O, Lycka A, Klasek A. Reaction of 1-substituted 3-aminoquinolinediones with isocyanic and isothiocyanic acid. Tetrahedron. 2011;67:2407–2413. doi: 10.1016/j.tet.2011.02.002. [Cross Ref]
7. Klasek A, Koristek K, Lycka A, Holcapek M. Unprecedented reactivity of 3-amino-1H,3H-quinoline-2,4-diones with urea: an efficient synthesis of 2,6-dihydroimidazo[1,5-c]quinazoline-3,5-diones. Tetrahedron. 2003;59:1283–1288. doi: 10.1016/S0040-4020(03)00028-0. [Cross Ref]
8. CrysAlis CCD, CrysAlis RED. Oxford diffraction. Yarnton: Oxford Diffraction Ltd; 2008.
9. Sheldrick GM. A short history of SHELX. Acta Cryst Sect A. 2008;64:112–122. doi: 10.1107/S0108767307043930. [PubMed] [Cross Ref]
10. Farrugia LJ. ORTEP-3 for windows - a version of ORTEP-III with a Graphical User Interface (GUI) J. Appl. Crystallogr. 1997;30:565–572. doi: 10.1107/S0021889897003117. [Cross Ref]
11. Dolomanov OV, Bourhis LJ, Gildea RJ, Howard JAK, Puschmann H. OLEX2: a complete structure solution, refinement and analysis program. J. Appl. Crystallogr. 2009;42:339–341. doi: 10.1107/S0021889808042726. [Cross Ref]
12. Becke AD. Density‐functional thermochemistry. III. The role of exact exchange. J Chem Phys. 1998;98:5648–5658. doi: 10.1063/1.464913. [Cross Ref]
13. Becke AD. Density‐functional thermochemistry. IV. A new dynamical correlation functional and implications for exact‐exchange mixing. J Chem Phys. 1996;104:1040–1048. doi: 10.1063/1.470829. [Cross Ref]
14. Lee C, Yang W, Parr RG. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B. 1988;37:785–791. doi: 10.1103/PhysRevB.37.785. [PubMed] [Cross Ref]
15. Ditchfield R, Hehre WJ, Pople JA. Self‐consistent molecular‐orbital methods IX. An extended gaussian‐type basis for molecular‐orbital studies of organic molecules. J Chem Phys. 1971;54:724–730. doi: 10.1063/1.1674902. [Cross Ref]
16. Krishnan R, Binkley JS, Seeger R, Pople JA. Self-consistent molecular orbital methods XX. A basis set for correlated wave functions. J Chem Phys. 1980;72:650–658. doi: 10.1063/1.438955. [Cross Ref]
17. Dunning TH. Gaussian basis sets for use in correlated molecular calculations I. The atoms boron through neon and hydrogen. J Chem Phys. 1989;90:1007–1015. doi: 10.1063/1.456153. [Cross Ref]
18. Kendall RA, Dunning TH, Harrison RJ. Electron affinities of the first‐row atoms revisited systematic basis sets and wave functions. J Chem Phys. 1992;96:6796–6802. doi: 10.1063/1.462569. [Cross Ref]
19. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA, Nakatsuji H, Caricato M, Li X, Hratchian HP, Izmaylov AF, Bloino J, Zheng G, Sonnenberg JL, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Vreven T, Montgomery JAJr, Peralta JE, Ogliaro F, Bearpark M, Heyd JJ, Brothers E, Kudin KN, Staroverov VN, Kobayashi R, Normand J, Raghavachari K, Rendell A, Burant JC, Iyengar SS, Tomasi J, Cossi M, Rega N, Millam MJ, Klene M, Knox JE, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Martin RL, Morokuma K, Zakrzewski VG, Voth GA, Salvador P, Dannenberg JJ, Dapprich S, Daniels AD, Farkas Ö, Foresman JB, Ortiz JV, Cioslowski J, Fox DJ (2009) Gaussian 09, revision D.01. Gaussian Inc., Wallingford
21. Hęclik K, Dębska B, Dobrowolski JCZ. RSC Adv. 2014;4:17337–17346. doi: 10.1039/c4ra02294g. [Cross Ref]
22. Dennington R, Keith T, Millam J (2009) Gauss View, version 5. Semichem Inc., Shawnee Mission, KS
23. Macrae CF, Edgington PR, McCabe P, Pidcock E, Shields GP, Taylor R, Towler M, van de Streek J. J. Appl. Crystallogr. 2006;39:453–457. doi: 10.1107/S002188980600731X. [Cross Ref]
24. Szyszkowska A, Hęclik K, Lubczak J, Trzybiński D, Woźniak K, Klasek A, Zarzyka I. Chemoselectivity of 1-phenyl-2H,6H-imidazo[1,5-c]quinazoline-3,5-dione reaction with oxiranes. J Mol Struct. 2017;1127:708–715. doi: 10.1016/j.molstruc.2016.08.025. [Cross Ref]
25. Spek AL. Structure validation in chemical crystallography. Acta Crystallogr. Sect. D. 2009;65:148–155. doi: 10.1107/S090744490804362X. [PubMed] [Cross Ref]
26. Wiberg KB, Hammer JD, Castejon H, Bailey WF, DeLeon EL, Jarret RM. Conformational studies in the cyclohexane series. 1. Experimental and computational investigation of methyl, ethyl, isopropyl, and tert-Butylcyclohexanes. J Org Chem. 1999;64:2085–2095. doi: 10.1021/jo990056f. [PubMed] [Cross Ref]
27. Steiner T. The hydrogen bond in the solid state. Angew Chem Int Ed. 2002;41:48–76. doi: 10.1002/1521-3773(20020104)41:1<48::AID-ANIE48>3.0.CO;2-U. [PubMed] [Cross Ref]
28. Perrin CL, Nielson JB. Strong hydrogen bonds in chemistry and biology. Annu Rev Phys Chem. 1997;48:511–544. doi: 10.1146/annurev.physchem.48.1.511. [PubMed] [Cross Ref]
29. Bickelhaupt FM, Baerends EJ. The case for steric repulsion causing the staggered conformation of ethane. Angew Chem Int Ed. 2003;42:4183–4188. doi: 10.1002/anie.200350947. [PubMed] [Cross Ref]
30. Leventis N, Hanna SB, Sotiriou-Leventis C. A three-dimensional energy surface for the conformational inversion of cyclohexane. J Chem Educ. 1997;74:813–814. doi: 10.1021/ed074p813. [Cross Ref]

Articles from Springer Open Choice are provided here courtesy of Springer