Crystallographic and Computational Study of Purine : Caffeine Derivative

The crystal structure of substituted purine derivative, 8-(3-butyl-4-phenyl-2,3-dihydrothiazol-2-ylidene)hydrazino-3,7-dihydro1,3,7-trimethyl-1H-purine-2,6-diones, caffeine derivative, has been determined. It crystallized in monoclinic system and space group P2 1 /c with unit cell parameters a = 15.2634 (9), b = 13.4692 (9), c = 11.9761 (7) Å, and β = 108.825 (3). Although each constituting moiety of the structure individually is planar, nonplanar configuration for the whole molecule was noticed. Molecular mechanics computations indicated the same nonplanar feature of thewholemolecule. A network of intermolecular hydrogen bonds contacts and π interactions stabilized the structure.

Purine ring system is one of the most heterocyclic ring systems in nature possessing the potential to impact several areas, such as a better understanding of the biological effect of DNA damaging agents, enzymes/substrate interactions, and the development of more potent medicines, such as antineoplastic (anti abnormal tissue growth, tumor), antileukemic (blood cancer), anti-HIV (anti-human immunodeficiency virus), and antimicrobial.Moreover, caffeine has been found to enhance anticancer activity of some chemotherapeutic agents and ionizing radiation.Previously, it was reported that methylxanthines may protect cells against the cytotoxic (poisonous to cells) effects and significantly decrease the mutagenicity of the anticancer aromatic drugs [2].
Recently, the activity of the present compound has been studied and it was found to have a great pharmaceutical and pharmacological interest [2].Neither the crystal structure nor the molecular modeling has been reported yet.Therefore, in the present work, we introduce the X-ray crystallographic analysis of the molecular structure and discuss the results with the molecular mechanics computations; such results should be valuable information in the fields of pharmaceutics and pharmacology.
In computational terms, molecular mechanics is the least expensive and fastest method.It is a method to calculate the structure and energy of molecules for providing excellent structural parameters in terms of bond distances, angles, and so forth, for the most stable conformation of the molecules [3,4].Also, many literatures have mentioned successful applications for combining X-ray single crystal structure analysis with computational studies, such as molecular mechanics in vacuo calculations (M.M.) [5][6][7].

Materials and Methods
2.1.Synthesis.The target compound was prepared as the reported procedure [2] (Scheme 1), and the IR spectra were recorded using KBr discs on a Perkin-Elmer 1430 spectrophotometer.Melting point was determined in open-glass capillaries on a Gallenkamp melting point apparatus and is uncorrected.The starting 8-hydrazinocaffeine (1) was prepared by the treatment of 8-chlorocaffeine with hydrazine hydrate.A mixture of 8-(N-butylthiocarbamoylhydrazino)-3,7-dihydro-1,3,7-trimethyl-1H-purine-2,6-dione (2 mmole) and phenacyl bromide (2 mmole) in absolute ethanol (20 mL) was heated under reflux for 30 minutes.The reaction mixture was then concentrated and left to cool to room temperature.The separated crystalline product was filtered, dried, and recrystallized from ethanol.

X-Ray Single Crystal
Measurements.Crystal was selected and checked for imperfections such as cracks, bubbles, twining, or voids and mounted onto thin glass fibers and glued with epoxy glue.X-ray diffraction data was collected at room temperature on an Enraf-Nonius 590 Kappa CCD single crystal diffractometer with graphite monochromated Mo-K ( = 0.71073 Å) radiation, at the National Research Center of Egypt [8,9].
Cell refinement and data reduction were carried using Denzo and Scalepak programs [10].The crystal structure was solved by direct method using SHELXS-97 program [11,12] which revealed the positions of all nonhydrogen atoms and refined by the full matrix least squares refinement based on F 2 using maXus package [13].The anisotropic displacement parameters of all nonhydrogen atoms were refined, and then the hydrogen atoms were introduced as a riding model with C-H = 0.96 Å and refined isotropically.The molecular graphics were prepared using ORTEP [14], Diamond [15], and Qmol [16] programs.
The crystal data is listed in  by Steepest Descents Method.The conformational energy of the molecule was calculated.

Results and Discussions
3.1.Crystal Structure Description.An ORTEP diagram with 50% probability displacement ellipsoids of the molecular structure is shown in Figure 1.The structure consists mainly of caffeine group linked with thiazole ring through hydrazine at C8, thiazole ring attached with butane, and phenyl ring.Each constituting moiety of the compound individually showed planar configuration; which can be noticed from the best plane calculations; where the maximum deviations were corresponding to C11, −0.0251 (5) in thiazole, C16, 0.0099 (6) in phenyl, and C7, 0.0967 (5) Å in caffeine group.However, nonplanar configuration for the whole molecule is noted, whereas the dihedral angle between phenyl ring and thiazole ring is 66.22(3) ∘ .Also, the tilting angle between the thiazole ring and caffeine group is 68.94(4) ∘ and butane plane is 84.98(4) ∘ .The nonplanarity feature that has been observed in this molecule may be attributed to the effect of the steric hindrance interaction between the different moieties composing it and the heavy substitution effect of butane moiety at N1.
The average values of bond lengths and angles are almost within the expected range and in consent with similar structures [20].The structure packing was stabilized by a network of intermolecular hydrogen bonds contacts, N-H⋅ ⋅ ⋅ Br, and  interactions between C13-H13 and the center of gravity (cg) of C12-C17 ring, as shown in Table 2 and Figure 2.

Molecular Mechanics.
Figure 3 represents the obtained molecular structure of the present compound using molecular mechanics calculations; comparison with that obtained crystallographically is given in Figure 4.The minimum energy structure obtained by molecular mechanics of the investigated compound to some extent matches the crystal structures obtained experimentally.The global minimum energy conformation of the molecule has the values 30.3 and 31.46 kcal/mol, for the molecular structures obtained experimentally (Exp.) and using molecular mechanics (M.M.) respectively.Unexpectedly, the crystal structure has the lower value, although it was reported that the crystal structure causes the molecules to adopt higherenergy conformations, which corresponds to local minima in the molecular potential energy surface.This is may be due to that the total potential energy in crystal structure of the present compound is affected by the steric hindrance.
It is noticeable from Table 3 and Figure 3 that the dimensions of the caffeine group, thiazole ring, phenyl ring, and butane obtained theoretically agree to some extent with those obtained experimentally with X-ray diffraction.However,   the structure obtained theoretically as a whole is not in a good agreement with the X-ray crystal structure.Moreover, the values of C10-C11-C12-C13 and C10-C11-N1-C18 torsion angles have a considerable difference between theoretical and experimental results, as shown in Figure 4.This variation may be because of the steric hindrance effect within the whole molecule, which has been noticed in the crystal structure.Also, the effect of the above-mentioned intermolecular contacts (Figure 2) is in agreement with what have been found in similar studies [5][6][7].

Scheme 1 :
Scheme 1: Chemical diagram of synthesis of the target compound.

Figure 1 :
Figure 1: The 50% probability displacement ellipsoids representation of the present compound.

Figure 2 :Figure 3 :
Figure 2: The molecular packing of the title compound with the intermolecular interactions.The N-H⋅ ⋅ ⋅ Br and C-H⋅ ⋅ ⋅  (cg) interactions are shown as green and brown dashed lines, respectively.

Figure 4 :
Figure 4: The X-ray crystal structure (green) and the structure of the one obtained by M.M. (red) of the target compound.

Table 1 .
The crystallographic supplementary data of C 21 H 26 BrN 7 O 2 S can be obtained

Table 1 :
Crystal data of the title compound.
[4,18,19]ular Mechanics Computations.Molecular mechanics in vacuo computations were carried out using HyperChem package[17].The Molecular Mechanics (MM+) force field was used as it is developed principally for organic molecules[4,18,19].The process of energy minimization was carried out

Table 3 :
Selected geometrical values of the experimentally and molecular mechanics obtained structures of the compound.