Mechanochemical Synthesis of the Catechol-Theophylline Cocrystal: Spectroscopic Characterization and Molecular Structure

Pharmaceutical cocrystallization offers the possibility to modify the physicochemical and biopharmaceutical properties of active pharmaceutical ingredients. The mechanochemical synthesis and spectroscopic characterization of the catechol-theophylline (CAT-TEO) cocrystal is reported. The cocrystal was prepared by the solvent-assisted grinding method. The ATR-IR spectroscopy study allowed to determine the formation of the cocrystal because the O-H and C=O stretching bands in the CAT-TEO cocrystal were shifted with respect to the starting materials, suggesting the formation of the C=O···H-O hydrogen bond interaction. Infrared spectroscopy also allowed to discard hydration of the cocrystal, and polymorphic transitions of the starting products as a consequence of the mechanochemical grinding. The X-ray powder diffraction and thermal studies confirmed the formation of a new solid phase. In the solid state 13C NMR spectra of the cocrystal, the signals were shifted with respect to the starting products. The 13C NMR chemical shifts of the CAT-TEO cocrystal were simulated by using the gauge including the atomic orbital (GIAO) method. These results showed a good correlation between the experimental and calculated 13C NMR results. Theoretical calculations and natural bonding orbital analysis (NBO) at a B3LYP/6-31G(d,p) level of theory were performed to obtain structural information of the cocrystal.


Introduction
Theophylline is a methyl-xanthine drug that occurs naturally in products such as coffee, tea and cocoa nuts, and is used for the treatment of asthma and chronic obstructive pulmonary disease [1]. Molecular recognition of theophylline has been exploited in solution and a solid state. In solution, it has led to the development of molecular receptors [2,3]; meanwhile, in solid state, it has been applied in the design and synthesis of pharmaceutical cocrystals. Pharmaceutical cocrystals are multi-component crystalline solids formed by a pharmaceutical active ingredient (API) and a cocrystal former, stabilized by noncovalent interactions, principally hydrogen bonds [4].
Theophylline is a model molecule [5] for cocrystallization studies because bears in their chemical structure carbonyl groups are able to act as hydrogen bond acceptors (theophylline can also act as a donor because it has a N-H group), allowing the formation of cocrystals with different cocrystal coformers, such as: organic acids, diclofenac, nicotinamide, acesulfame, to mention some [6][7][8][9].
On the other hand, catechol or 1,2-dihydroxybenzene contain hydroxyl groups, which can act as hydrogen bond donors for the synthesis of cocrystals due to their ability to form O-H···X (X = O, N) hydrogen bonds, with the advantage that they form neutral hydrogen bonds, discarding the formation of salts. Some examples of API's that have formed cocrystals with catechol are nalidixic acid, tenoxicam, brexpiprazole, and isoniazid [10][11][12][13].
From the reported cocrystallization methods, the most used are the solvent evaporation method and the solvent-assisted grinding method [14]. In the solvent evaporation method, the API and the cocrystal former are dissolved in a suitable solvent, which is left to evaporate to obtain the cocrystal. This method is used to obtain suitable monocrystals for single crystal X-ray diffraction; however, it depends on the solubility of the compounds, it is necessary to wait a few days for the evaporation of the solvent, and the presence of the solvent can lead to the formation of solvates. On the other hand, in the solventassisted method, the API and the cocrystal former are ground together with some solvent drops. The mechanochemical energy favors the rearrangement of the molecules in solid state to form the cocrystal. This method offers the advantage of reproducibility, rapid screening and independence of the solubility of the starting products [14].
Solid-state monitoring techniques such as powder and single X-ray diffraction, solidstate nuclear magnetic resonance (NMR) and infrared spectroscopy (IR) are employed in the characterization of pharmaceutical cocrystals [15]. IR spectroscopy has been used to determine the formation of cocrystals, polymorphs, salts and solvates [16][17][18]. Attenuated total reflection (ATR) IR spectroscopy offers the advantage of rapid determination (the sample does not need previous preparation for the analysis) of the formation of the new solid between the API and the cocrystal former by the shift of the IR absorption bands of the functional groups in the cocrystal with respect to the starting components [19,20]. Solid state nuclear magnetic resonance is used for the characterization of pharmaceutical cocrystals and polymorphs because it gives structural and crystallographic information and is sensitive to hydrogen bonding [21,22].
In order to exploit the ability of the dihydroxybenzenes to form O-H···O=C hydrogen bond interactions with theophylline, here we report the mechanochemical synthesis of the catechol-theophylline cocrystal (CAT-TEO) ( Figure 1) obtained by the solvent-assisted grinding method. The ATR-IR spectroscopy has played an important role because it allowed to rapidly determine the formation of the powder cocrystal, discarding the hydrate formation. The characterization of the CAT-TEO cocrystal was completed using solidstate 13 C NMR, powder X-ray diffraction and DSC calorimetry. The molecular structure of the CAT-TEO cocrystal was obtained by theoretical calculations.

Mechanochemical Synthesis
Theophylline anhydrous, catechol and dichloromethane were purchased from Aldrich. All the compounds were used as received.
The CAT-TEO cocrystal was prepared by grinding 0.486 g of anhydrous theophylline (2.72 mmol) with 0.300 g of catechol (2.7 mmol) in a pestle with a mortar for 15 min. Before the start of grinding, 0.5 mL of dichloromethane was added; then, the mixture was ground for 3 min. After grinding, the powder was collected in the center of the mortar. After that, 0.5 mL of dichloromethane was added to the powder mixture and ground for an additional 3 min. The cycle of adding 0.5 mL of dichloromethane and grinding for 3 min was repeated five times until 15 min of grinding was completed. The powder obtained was stored in a glass vial.
In order to discard polymorphic transitions as a consequence of the grinding, the starting products were ground individually (2.72 mmol) with dichloromethane under the same conditions of the CAT-TEO mixture.
To evaluate the water absorption effect over the catechol and theophylline, 2.72 mmol of each starting product were ground with 0.5 mL of water for 15 min. Water was added before the start of grinding. The 1:1 mixture of CAT-TEO was treated in the same way.
All the attempts to obtain single crystals for diffraction were unsuccessful.

Instruments
Infrared spectra of solid samples of the API s, the crystal formers and the ground products were obtained neat in a Bruker Tensor-27 spectrophotometer (16 scans, spectral range 600 to 4000 cm −1 , resolution 4 cm −1 ) equipped with an attenuated total reflection (ATR) system. IR spectra were acquired and processed using the OPUS software version 7.2 (Bruker Optik GmbH).
Solid-state 13C NMR spectra were recorded on a Bruker 400 Avance III ( 13 C, 100 MHz) instrument with a Cross Polarization Magic/Angle Spinning (CP/MAS) unit at 25 • C. The spinning rate and acquisition time were 8 kHz and 32 msec. A 4-mm bullet-type Kel-F zirconia rotor contained about 100 mg of sample. The recycle time of the pulse was 3 s. The signal of the adamantane was used as external reference (δ = 38.48 ppm).
Differential scanning calorimetry was performed on a TA Instruments Q2000 DSC. Temperature calibration was performed using an indium standard. Three milligrams of the sample were placed in an aluminum pan with a gradient of 5.00 • C /min from 50 • C to 200 • C. Thermal gravimetric analysis was performed using a TA instruments Q5000IR analyzer. The samples were loaded in an aluminum pan. The samples were heated at the rate of 5.00 • C/min, under nitrogen.

Theoretical Calculations
Two optimized geometries models as CAT-TEO 1 and CAT-TEO 2 were performed using the density functional theory (DFT) calculations, with the B3LYP hybrid functional and the 6-31G(d,p) basis set, using the Gaussian 09 software [23]. Vibrational frequencies were calculated to ensure that the optimized geometries correspond to local minima on the potential energy surface. Natural bond orbital (NBO) Fock matrix analysis at the same level of theory was performed to elucidate the electron density distribution within the molecules of the complexes, while 13 C NMR chemical shifts and magnetic shieldings were simulated using the gauge including atomic orbital (GIAO) method at B3LYP/6-31G(d,p) and referenced with the TMS HF/6-31G(d) option of the GaussView software [23].

IR Spectroscopy
ATR-IR spectroscopy study was performed to determine the formation of the new crystalline solid-phase solid by monitoring the shifts of the hydroxyl (O-H), and carbonyl (C=O) stretching frequencies in the CAT-TEO ( Figure 2) ground solid with respect to the starting products (CAT and TEO, Figure 2). The IR spectra were acquired directly after the CAT-TEO ground process, and from the bulk CAT and TEO. TEO and CAT were ground in cycles of 3 min, until a total of 15 min of grinding time. IR spectra of CAT and TEO were assigned according to previous reports [16,18].
(C=O) stretching frequencies in the CAT-TEO ( Figure 2) ground solid with respect to the starting products (CAT and TEO, Figure 2). The IR spectra were acquired directly after the CAT-TEO ground process, and from the bulk CAT and TEO. TEO and CAT were ground in cycles of 3 min, until a total of 15 min of grinding time. IR spectra of CAT and TEO were assigned according to previous reports [16,18].
In the IR spectrum of CAT-TEO obtained after 15 min of grinding ( Figure 2), the O-H stretching bands (Table 1) were shifted with respect to the starting products and the physical mixture (CAT-TEO-MIX, Figure 2). This suggests the formation of a C=O···H-O hydrogen bond between CAT and TEO. The free form of CAT showed two O-H stretching bands ( Figure 2)-one in the 3500 to 3400 cm −1 range and the other in the 3300 cm −1 range of the phenolic O-H stretching frequency. As a consequence of the rearrangement in the intermolecular hydrogen bonding pattern from O-H···O-H in the free form to the O-H···O=C in the cocrystal, the band in the 3500 to 3400 cm −1 range was shifted to lower frequencies with Δ(νO-H) = −293 cm −1 . Meanwhile, the band in the 3300 cm −1 range shifted to higher frequencies with Δ(νO-H) = 105 cm −1 . This shift is in agreement with previous reports about molecular complexes formed by O-H···O=C hydrogen bond interactions between catechol with oxalamates and carbamates [24,25].  In the IR spectrum of CAT-TEO obtained after 15 min of grinding (Figure 2), the O-H stretching bands (Table 1)    The carbonyl C=O stretching bands of TEO ( Figure 2) were shifted in the IR spectrum of the CAT-TEO ground mixture as a consequence of the formation of the cocrystal (Table 1), with absolute values of ∆(νC=O) less than 20 cm −1 . This shift is in accordance with the formation of the neutral C=O···H-O hydrogen bond interaction, discarding the proton transfer [16].
The solid form of TEO presents polymorphism, and each polymorph gives a different vibrational spectrum [26]. On the other hand, there are not any reports about polymorphs of CAT. To discard polymorphic transitions as a consequence of the grinding, CAT and TEO were ground individually with dichloromethane. After 15 min of grinding, the IR spectra were similar to those obtained before grinding, which discards polymorphic transformations ( Figure S1 of supplementary material) as a consequence of the grinding.
Possible effects of water absorption from the air during mechanochemical grinding over the starting products and the mixtures were evaluated by grinding CAT and TEO with water (CAT-H and TEO-H, Figure 2); the CAT-TEO mixture was ground with water (CAT-TEO-H, Figure 2). The IR spectrum of TEO ground with water, recorded after 15 min of grinding (TEO-H, Figure 2), was different with respect to the anhydrous form (TEO, Figure 2), showing shifts in the N-H, C=O and C-H(alkyl) stretching frequencies, and the appearance of the broad O-H band in the range of the 3350 to 3360 cm −1 . The spectrum (TEO-H, Figure 2) is similar to that reported for the TEO hydrate [18]. Two characteristic bands in the IR spectrum of TEO-H are the N-H stretching frequency at 3139 cm −1 and the C-H alkyl stretching frequency at 3106 cm −1 .
Contrarily to TEO, the IR frequencies of CAT ground with water (CAT-H, Figure 2) did not show any significant shift with respect to the anhydrous form. The O-H band of the water appeared to overlap with the O-H bands of CAT. It is noteworthy that after the grinding of CAT-H, a liquid product (catechol solution) was obtained, so the IR spectrum of CAT-H was obtained in liquid form. The fact that water acts as an interference in the IR spectrum of CAT-H means that it is not part of the crystalline lattice of the solid, discarding the hydration. This is in agreement with the absence of reports of hydrated solid forms of CAT.
The IR spectrum of the CAT-TEO 1:1 mixture ground with water (CAT-TEO-H, Figure 2) showed a different pattern with respect to the anhydrous ground mixture (CAT-TEO, Figure 2). The IR spectrum of CAT-TEO-H showed two O-H bands at 3517 and 3454 cm −1 ; meanwhile, CAT-TEO showed a single O-H band at 3151 cm −1 . In the same way, the IR spectrum of CAT-TEO-H showed a different pattern of C=O bands with respect to the CAT-TEO spectrum. This allows determining that the CAT-TEO-H ground solid is the cocrystal hydrate.

Powder X-ray Diffraction
In order to confirm the formation of a new solid phase belonging to the CAT-TEO cocrystal, powder X-ray diffraction studies (PXRD) were performed. The PXRD pattern of CAT-TEO ground with dichloromethane after 15 min of grinding time was different to the starting solids CAT and TEO, and also to the hydrated CAT-TEO-H ground mixture (Figure 3). The PXRD pattern of TEO-H is similar to that reported for the theophylline hydrate [18]. The absence of the signals of CAT at low 2θ in the XRPD diffractograms of the anhydrous CAT-TEO ground mixture suggests the complete conversion of CAT and TEO into the CAT-TEO cocrystal. The PXRD patterns of CAT-TEO showed characteristic peaks at 2θ = 9.0 • 10.5 • , 13 CAT-TEO ground with dichloromethane after 15 min of grinding time was different to the starting solids CAT and TEO, and also to the hydrated CAT-TEO-H ground mixture ( Figure 3). The PXRD pattern of TEO-H is similar to that reported for the theophylline hydrate [18].

Thermal Analysis
Differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) experiments were carried out to study the thermal behavior of the CAT-TEO cocrystal. DSC results revealed that the melting point of the CAT-TEO cocrystal (134.16 °C) was different from CAT (100 to 103 °C) [27] and TEO (270 to 274 °C) [28], indicating the formation of a new crystalline solid phase. Meanwhile, the TGA curve of the CAT-TEO cocrystal indicates the absence of hydration (DSC and TGA curves are shown in Figure S2 of supplementary material).

Solid-State 13 C NMR
Solid-state 13 C nuclear magnetic resonance spectra of CAT, TEO and the CAT-TEO cocrystal ( Figure 4) were obtained to acquire information about hydrogen bonding interactions and crystallography. The 13 C NMR chemical shifts of the CAT and TEO starting products were assigned according to previous reports [8,29]. The 13 C NMR spectrum of free CAT showed the half of total signals due to the C2 symmetry.
Most of the 13 C NMR signals in the spectrum of the CAT-TEO cocrystal were shifted with respect to the starting CAT and TEO (Table 2); this is due to the change of the chemical environment as a consequence of the cocrystal formation. Carbonyl carbons of TEO are sensitive to hydrogen bonding formation. The carbonyl C2 signal of TEO was shifted in the CAT-TEO cocrystal with respect to the free TEO (Δδ 13 C = 3.6 ppm), due to the formation of the C=O···H-O hydrogen bond interaction. The C5 signal was sensitive to the change of chemical environment with Δδ 13 C = −3 ppm. C6 and C7 also showed significant shifts C6 Δδ 13 C = −2.4 ppm, and C7 Δδ 13 C = 1.9 ppm.

Thermal Analysis
Differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) experiments were carried out to study the thermal behavior of the CAT-TEO cocrystal. DSC results revealed that the melting point of the CAT-TEO cocrystal (134.16 • C) was different from CAT (100 to 103 • C) [27] and TEO (270 to 274 • C) [28], indicating the formation of a new crystalline solid phase. Meanwhile, the TGA curve of the CAT-TEO cocrystal indicates the absence of hydration (DSC and TGA curves are shown in Figure S2 of supplementary material).

Solid-State 13 C NMR
Solid-state 13 C nuclear magnetic resonance spectra of CAT, TEO and the CAT-TEO cocrystal ( Figure 4) were obtained to acquire information about hydrogen bonding interactions and crystallography. The 13 C NMR chemical shifts of the CAT and TEO starting products were assigned according to previous reports [8,29]. The 13 C NMR spectrum of free CAT showed the half of total signals due to the C2 symmetry.
Most of the 13 C NMR signals in the spectrum of the CAT-TEO cocrystal were shifted with respect to the starting CAT and TEO (Table 2); this is due to the change of the chemical environment as a consequence of the cocrystal formation. Carbonyl carbons of TEO are sensitive to hydrogen bonding formation. The carbonyl C2 signal of TEO was shifted in the CAT-TEO cocrystal with respect to the free TEO (∆δ 13 C = 3.6 ppm), due to the formation of the C=O···H-O hydrogen bond interaction. The C5 signal was sensitive to the change of chemical environment with ∆δ 13 C = −3 ppm. C6 and C7 also showed significant shifts C6 ∆δ 13 C = −2.4 ppm, and C7 ∆δ 13 C = 1.9 ppm. Appl. Sci. 2021, 11, x FOR PEER REVIEW 7 of 11  All the 13 C NMR signals of CAT in the CAT-TEO spectrum were shifted with respect to free CAT as a consequence of the formation of the cocrystal ( Table 2). The C13 and C16 signals appeared to be split, indicating an asymmetric chemical environment in these carbons.

Theoretical Calculations
To validate the information obtained from the 13 C CP/MAS NMR study (δexp), the gauge including the atomic orbital (GIAO) method at a B3LYP/6-31G(d,p) level of theory was used to calculate the isotropic magnetic shielding (σiso) and predict the 13 C NMR chemical shifts (δcalc) of the CAT-TEO1 and CAT-TEO2 models (Table 3, Figure 5). The calculated chemical shifts are in agreement with the experimental results because the correlation coefficient R 2 ( Figure S3 supplementary material) is 0.997 and 0.996, respectively.   All the 13 C NMR signals of CAT in the CAT-TEO spectrum were shifted with respect to free CAT as a consequence of the formation of the cocrystal ( Table 2). The C13 and C16 signals appeared to be split, indicating an asymmetric chemical environment in these carbons.

Theoretical Calculations
To validate the information obtained from the 13 C CP/MAS NMR study (δexp), the gauge including the atomic orbital (GIAO) method at a B3LYP/6-31G(d,p) level of theory was used to calculate the isotropic magnetic shielding (σiso) and predict the 13 C NMR chemical shifts (δcalc) of the CAT-TEO 1 and CAT-TEO 2 models (Table 3, Figure 5). The calculated chemical shifts are in agreement with the experimental results because the correlation coefficient R 2 ( Figure S3 supplementary material) is 0.997 and 0.996, respectively.
Density functional theory methods (DFT) allow the calculation of molecular properties of the molecular complexes. One of these properties is the optimization of the geometry in the gas phase, which gives an idea about the molecular structure of a complex [30,31]. In the same way, natural bond orbital (NBO) analysis gives information about the electron density distribution involved in the formation of hydrogen bonds among the molecules of the complexes [32,33]. Table 3. Calculated (δcalc) and experimental (δexp) CP/MAS solid state 13 C NMR chemical shifts (ppm) of TEO and CAT in the 1:1 anhydrous ground mixture, and calculated (root mean square error less than 4.0) GIAO isotropic magnetic shielding tensors (σiso). Density functional theory methods (DFT) allow the calculation of molecular properties of the molecular complexes. One of these properties is the optimization of the geometry in the gas phase, which gives an idea about the molecular structure of a complex [30,31]. In the same way, natural bond orbital (NBO) analysis gives information about the electron density distribution involved in the formation of hydrogen bonds among the molecules of the complexes [32,33].

C1
In order to obtain structural information about the intermolecular hydrogen bond patterns involved in the formation of the CAT-TEO cocrystal, DFT calculations at B3LYP/6-31G(d,p) level of theory were performed in CAT-TEO1 and CAT-TEO2 models ( Figure 5) considering the two possible C=O···H-O hydrogen bond patterns in the CAT-TEO cocrystal. The CAT-TEO2 model is the most stable structure (by only 2.404 kcal/mol) involving the C=O(2)···H-O (11) and the N-H···O(11) hydrogen bond interactions ( Figure  5); meanwhile, the CAT-TEO1 model is in accordance with the 13 C NMR shifts because the carbonyl C2 signal showed a Δδ 13 C = 3.6 ppm in the CAT-TEO cocrystal with respect to the free TEO. This deshielding suggests the formation of the C=O···H-O hydrogen bonding in this position, contrary to the carbonyl C1, which showed a Δδ 13 C = −0.4 ppm. This is also in accordance with previous reports about synthon hierarchies in cocrystals of TEO with hydroxybenzoic acids with respect to the hydrogen bonds involving phenolic O-H groups, in which the O-H···O=C(urea) interaction is preferred over the O-H···O=C(amide) interaction because the O=C(amide) is involved in the dimerization of TEO by the N-H···O=C(amide) hydrogen bond synthon [5].  Table S1 of supplementary material). The superficial analysis at DFT-B3LYP/6-31G(d,p) confirmed the preference for the C=O(2)···H-O in the CAT-TEO2 model cocrystal (Figure 6), as suggested by the 13 C NMR In order to obtain structural information about the intermolecular hydrogen bond patterns involved in the formation of the CAT-TEO cocrystal, DFT calculations at B3LYP/6-31G(d,p) level of theory were performed in CAT-TEO 1 and CAT-TEO 2 models ( Figure 5) considering the two possible C=O···H-O hydrogen bond patterns in the CAT-TEO cocrystal. The CAT-TEO 2 model is the most stable structure (by only 2.404 kcal/mol) involving the C=O(2)···H-O (11) and the N-H···O(11) hydrogen bond interactions ( Figure 5); meanwhile, the CAT-TEO 1 model is in accordance with the 13 C NMR shifts because the carbonyl C2 signal showed a ∆δ 13 C = 3.6 ppm in the CAT-TEO cocrystal with respect to the free TEO. This deshielding suggests the formation of the C=O···H-O hydrogen bonding in this position, contrary to the carbonyl C1, which showed a ∆δ 13 C = −0.4 ppm. This is also in accordance with previous reports about synthon hierarchies in cocrystals of TEO with hydroxybenzoic acids with respect to the hydrogen bonds involving phenolic O-H groups, in which the O-H···O=C(urea) interaction is preferred over the O-H···O=C(amide) interaction because the O=C(amide) is involved in the dimerization of TEO by the N-H···O=C(amide) hydrogen bond synthon [5].
The C=O···H-O hydrogen bond interaction in the CAT-TEO 2 model is formed by the overlap of the lone pair electrons of the carbonylic oxygen (nO) with the σ*(H-O) anti-bonding orbital (the energetic parameters are shown in Table S1 of supplementary material). The superficial analysis at DFT-B3LYP/6-31G(d,p) confirmed the preference for the C=O(2)···H-O in the CAT-TEO 2 model cocrystal (Figure 6), as suggested by the 13 C NMR results because in the spectra the C2 carbon showed a higher shift with respect to the C1. The hydrogen bond distances are similar to the reported for cocrystals of dihydroxybenzenes [15,24].
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 11 results because in the spectra the C2 carbon showed a higher shift with respect to the C1. The hydrogen bond distances are similar to the reported for cocrystals of dihydroxybenzenes [15,24].

Conclusions
The CAT-TEO cocrystal was prepared by the solvent-assisted grinding method and characterized by infrared spectroscopy, powder X-ray diffraction and solid-state CP/MAS 13 C NMR. The infrared spectroscopy study allowed to determine that the formation of the cocrystal, directed by C=O···H-O hydrogen bond interactions and confirmed by 13 C solidstate NMR and theoretical calculations. Hydration of the CAT-TEO cocrystal was discarded because the IR spectrum of the anhydrous CAT-TEO mixture was different with respect to the hydrated ground mixture (CAT-TEO-H), and the TGA curve does not show hydration. IR spectroscopy allowed to determine the hydration of theophylline (TEO-H) and the CAT-TEO cocrystal (CAT-TEO-H) as a consequence of the grinding with water. The formation of the solid phase belonging to the CAT-TEO cocrystal was confirmed by powder X-ray diffraction. Thermal analysis revealed that the melting point of the CAT-TEO cocrystal is different from the starting compounds. The GIAO 13 C NMR chemical shift values were correlated with the experimental data.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Figure S1: Figure S1. IR spectra of (a) bulk catechol; (b) catechol ground with dichloromethane after 15 min of grinding; c) bulk theophylline; d) theophylline ground with dichloromethane after 15 min of grinding; Figure S2: DSC and TGA curves of the CAT-TEO cocrystal; Figure S3: Correlations of the experimental and theoretical shielding parameters for CAT-TEO1 and CAT-TEO2 models. Table S1: Theoretical optimization energy and second order perturbation theory analysis of Fock matrix in NBO basis of the optimized CAT-TEO models.

Conclusions
The CAT-TEO cocrystal was prepared by the solvent-assisted grinding method and characterized by infrared spectroscopy, powder X-ray diffraction and solid-state CP/MAS 13 C NMR. The infrared spectroscopy study allowed to determine that the formation of the cocrystal, directed by C=O···H-O hydrogen bond interactions and confirmed by 13 C solid-state NMR and theoretical calculations. Hydration of the CAT-TEO cocrystal was discarded because the IR spectrum of the anhydrous CAT-TEO mixture was different with respect to the hydrated ground mixture (CAT-TEO-H), and the TGA curve does not show hydration. IR spectroscopy allowed to determine the hydration of theophylline (TEO-H) and the CAT-TEO cocrystal (CAT-TEO-H) as a consequence of the grinding with water. The formation of the solid phase belonging to the CAT-TEO cocrystal was confirmed by powder X-ray diffraction. Thermal analysis revealed that the melting point of the CAT-TEO cocrystal is different from the starting compounds. The GIAO 13 C NMR chemical shift values were correlated with the experimental data.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/app11093810/s1, Figure S1: Figure S1. IR spectra of (a) bulk catechol; (b) catechol ground with dichloromethane after 15 min of grinding; c) bulk theophylline; d) theophylline ground with dichloromethane after 15 min of grinding; Figure S2: DSC and TGA curves of the CAT-TEO cocrystal; Figure S3: Correlations of the experimental and theoretical shielding parameters for CAT-TEO 1 and CAT-TEO 2 models. Table S1: Theoretical optimization energy and second order perturbation theory analysis of Fock matrix in NBO basis of the optimized CAT-TEO models.

Conflicts of Interest:
The authors declare no conflict of interest.