Hirshfeld Surface Analysis and Energy Framework for Crystals of Quinazoline Methylidene Bridged Compounds

: The crystal structures of 4-(3,4-dimethoxyphenylethylamino)-methylidene-2,3,4,10-tetrahydro-1H-pyrido[2,1-b]-quinazolin-10-one ( 1 ) and 4-(3,4-methylene-dioxyphenylethylamino)-methylidene-2,3,4,10-tetrahydro-1H-pyrido[2,1-b]-quina-zolin-10-one hydrochloride ( 2 ) were studied by single crystal X-ray diffraction. Their molecular and crystal structures are described in the context of intra- and inter-molecular interactions and of stereoisomerism. The crystallographic study established mixed E , Z configuration about the C4=C12 bond for ( 1 ) and E configuration about the C4=C12 bond for ( 2 ). For molecular crystals, Hirshfeld surface analyses may provide insight into intermolecular interactions, and energy framework analyses allow one to quantify different contributions to the overall energy. These analyses were performed to pinpoint intermolecular interactions in ( 1 ) and ( 2 ). According to our results, the molecules are associated with intra- and intermolecular hydrogen bonds, C-H··· π and π -stacking interactions. The three-dimensional Hirshfeld surface analyses and two-dimensional fingerprint plots revealed that the structures are dominated by H···H, H···C/C···H and H···O/O···H contacts. The intermolecular energy analysis confirmed a significant contribution of dispersion to the stabilization of molecular packings in ( 1 ) and ( 2 ).


Crystallographic Details
Crystal diffraction data were collected: for (1) on a Bruker D8 diffractometer equipped with an APEX CCD detector and multilayer optics microfocus tube with MoKα radiation and for (2) on the EH1 Kappa diffractometer equipped with a Dectris CdTe area detector using synchrotron radiation at beamline P24, DESY, Hamburg. The crystals of (1) were extremely small; because of this, it was used for synchrotron radiation. Absorption correction was not applied for (1). Crystal structures of (1) and (2) were solved by direct methods with the program SHELXS [17] and refined by full-matrix least-squares on F 2 with the SHELXL [18] package.
In the crystal structure (1), H atoms attached to C were positioned geometrically and treated as riding on their parent atoms, with C-H = 0.95 (aromatic), 0.98 (methyl), 0.99 (methylene) Å and were refined with Uiso(H) = 1.5Ueq(C) for methyl H atoms and 1.2Ueq(C) otherwise. The H atom in the enamine group was refined with a distance restraint (target distance N-H = 0.95 Å) and with Uiso(H) = 1.2Ueq(N). The residual density is high (0.85 e Å -3 ) due to the extensive disorder.
In the crystal structure of (2), H atoms attached to C were positioned geometrically and treated as riding on their parent atoms, with C-H = 0.95 (aromatic), 0.99 (methylene) Å and were refined with 1.2Ueq(C) otherwise. The H atom bonded to nitrogen was constrained with a distance 0.88 A˚ and with Uiso(H) = 1.2Ueq(N).
Details of the crystallographic data collection, structural determination and refinement for (1) and (2) are summarized in Table S1 (see Supplementary Materials). The full dataset is available in form of CIF files, deposited with the CCDC (2039173, 2039174), and may be obtained free of charge via https://www.ccdc.cam.ac.uk/structures.

Hirshfeld Surface Calculations
The Hirshfeld surface (HS) is the area around a molecule in crystal space which separates two regions, that of the inner reference molecule from the outer neighboring molecules. The separation of space by the HS allows the analysis of intermolecular interactions of fingerprints in a crystalline medium [19,20]. HS analysis can be employed to visualize and quantify various non-covalent interactions that stabilize the crystal packing. HS can be mapped with different properties, namely, dnorm, electrostatic potential, shape index and curvature. The dnorm property is a symmetric function of distances to the surface between nuclei inside and outside the Hirshfeld surface (di and de, respectively), relative to their respective van der Waals radii. The regions with red and blue color on the dnorm represent the shorter and longer inter contacts while the white color indicates the contacts around the van der Waals radii. Two-dimensional fingerprint plots provide relevant information of intermolecular contacts in the crystal. HS analysis has become a very useful tool for explaining the nature of intermolecular interactions that affect the packing of molecules in crystals [21]. Hirshfeld surfaces analyses are widely used to study the phenomena of polymorphism [22,23], cocrystallization [24], the inclusion of small molecules in the cavities of macromolecules [25,26], and the search for correlations between the strength of interactions and the melting point [27].
The HS analyses and their associated 2D fingerprint plots (full and decomposed) [28] were carried out employing the CrystalExplorer 17.5 program [29]. The dnorm surface was mapped with the color scale in the range −0.050 au (red) to 0.600 au (blue). The 2D fingerprint plots (di vs. de) were displayed using the expanded 0.6-2.8 Å range.

Structural Description
The compounds (1) and (2) crystallize in the monoclinic system in space group P21/c. Displacement ellipsoid plots and numbering schemes for both molecules are provided in Figure 1. For the structure of (1), X-ray diffraction revealed that the binding enamine group is disordered, leading to the formation of E and Z isomers about the C4=C12 double bond. In case of the E configuration, an O···H-N intermolecular hydrogen bond links the N atom of the enamine group to the keto group of the quinazoline moiety and gives rise to a 8 chain motif connecting adjacent molecules to zigzag like chains propagating along the b-axis (Figure 2a). In the case of the Z configuration, an N9B'-H9B'···N5 intramolecular hydrogen bond (2.812(8) Å) is formed, which generates an 6 ring motif (Figure 2b). In both cases, an additional C6-H6···O2 hydrogen bond leads to a two-dimensional network parallel to (102) plane (Figure 2a,b and Figure S1). The hydrogenbond geometry for (1) is listed in Table S2.  Figure 3a shows the interaction between adjacent molecules along the crystallographic a-axis. Cg4 denotes the centroid of the 3,4-dimethoxyphenylethylamino ring (C1′-C6′). In the Zconfiguration of (1), both contacts C7B'-H7BB···Cg4 and C11'-H11B···Cg4 participate in the interaction between adjacent molecules, whereas in the case of the E-configuration, only the latter contact contributes to C-H···π interactions.
Structure (2) consists of Cl − counter anions and the organic cations. The latter are composed of tricyclic quinazoline and 3,4-methylenedioxyphenylethylamino units linked by an enamine group (Figure 1b). The crystallographic study of structure (2) revealed the E configuration about the C4=C12 double bond, similar to one of the alternative configurations in structure (1) (Figure 1b). N5-H5 Cl1 (3.088(5)Å) and N10'-10' Cl1(3.187(6)Å) intermolecular hydrogen bonds links the Clanion to N atoms in the enamine group and the pyrimidine ring and give rise to a centrosymmetric dimer ( Figure S2 and Table S3)
The most prominent type of contacts in (1) structure corresponds to H···H contacts; they contribute 56.5% to the overall surface contacts (Figure 5b). The presence of C-H···π interactions gives rise to a pair of characteristic wings in the fingerprint plot decomposed into C···H/H···C contacts (Figure 5c) contributing 21.4% to the HS. Short O···H/H···O contacts (Figure 5d) contribute 14.5% to the HS and form a pair of spikes. N···H/H···N contacts, contributing 2.6% to the overall surface, are depicted in Figure 5e as widely scattered wings.
For the structure (2), the fingerprint profile is dominated by H···H (44.4%) surface contacts (Figure 5g). C···H/H···C interactions contribute 17.0% of the total HS (Figure 5h (2) (i); N···H/H···C for (1) (e) and Cl···H/H···Cl for (2) (j). The outline of the full fingerprint plot is shown in grey. Surfaces to the right highlight the relevant dnorm surface patches associated with the specific contacts. The percentage of contribution is specified for each contact.

Interaction Energies
Interaction energies in (1) and (2) (1) and (2). The cylindrical radius is proportional to the relative strength of the corresponding energy and was adjusted to the scale factor of 80 kJ mol -1 with a cut-off value of 10 kJ mol −1 .

Conclusions
A detailed structural analysis of 4-(3,4-dimethoxyphenylethylamino)-methylidene-2,3,4,10tetrahydro-1H-pyrido[2,1-b]-quinazolin-10-one revealed disorder of the enamine group, leading to the concomitant presence of E and Z isomers in the crystal. In the E isomer, an intramolecular N-H···O hydrogen bond results in the formation of a C(8) ring motif. In the Z isomer, an intramolecular N-H···N hydrogen bond results in the formation of an S(6) pattern. The large number of H···H, H···C/C···H and H···O/O···H contacts revealed by a HS analysis indicate that van der Waals interactions and hydrogen bonding mostly contribute to crystal packing. Initial results on energy frameworks suggest that dispersion forces are relevant to the topology of the overall interaction energies in these crystals of quinazoline derivatives.
Author Contributions: A.T., S.Z., U.E. and K.T. conceived and designed the experiments, conceptualized the work, and prepared the manuscript for publication; U.E. conducted the X-ray analysis, reviewed, and edited the manuscript; U.E., R.W., B.T. and K.T. validated and formal analyzed; I.K. crystallization experiments; S.Z. and V.V. provided chemical synthesis; U.E. reviewed, and edited the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This work was supported by the Istedod Foundation of the Republic of Uzbekistan.