Study on four polymorphs of bifendate based on X-ray crystallography.

Bifendate, a synthetic anti-hepatitis drug, exhibits polycrystalline mode phenomena with 2 polymorphs reported (forms A and B). Single crystals of the known crystalline form B and 3 new crystallosolvates involving bifendate solvated with tetrahydrofuran (C), dioxane (D), and pyridine (E) in a stoichiometric ratio of 1:1 were obtained and characterized by X-ray crystallography, thermal analysis, and Fourier transform infrared (FT-IR) spectroscopy. The differences in molecular conformation, intermolecular interaction and crystal packing arrangement for the four polymorphs were determined and the basis for the polymorphisms was investigated. The rotation of single bonds resulted in different orientations for the biphenyl, methyl ester and methoxyl groups. All guest solvent molecules interacted with the host molecule via an interesting intercalative mode along the [1 0 0] direction in the channel formed by the host molecules through weak aromatic stacking interactions or non-classical hydrogen bonds, of which the volume and planarity played an important role in the intercalation of the host with the guest. The incorporation of solvent-augmented rotation of the C-C bond of the biphenyl group had a striking effect on the host molecular conformation and contributed to the formation of bifendate polymorphs. Moreover, the simulated powder X-ray diffraction (PXRD) patterns for each form were calculated on the basis of the single-crystal data and proved to be unique. The single-crystal structures of the four crystalline forms are reported in this paper.

Introduction

Polycrystalline phenomena with solvates are not uncommon for many known drugs. The incorporation of solvent molecules results in the formation of new crystalline forms. This has attracted the attention of drug researchers. It is well known that the polymorphs of a substance are chemically identical but can differ significantly in their pharmaceutically relevant properties, such as solubility1, 2, dissolution rate, and density, which can affect their bioavailability, membrane permeability and efficacy in drug formulations. The term “polymorphism” describes a phenomenon in which crystal structures of a substance are defined by different unit cells, and yet each of the forms has the same elemental composition. Along this same line, “solvatomorphism” encompasses crystals of a substance defined by different unit cells, but these unit cells differ in their elemental composition through the inclusion of solvent molecules. For greater insight into the regularities of solvents in crystal lattices, single-crystal X-ray diffraction (SCXRD) is suitable for detecting solvent structure and molecular interactions with the drug substance, and is particularly useful in the research of polymorphs. In drug development, the formation of solvates has some relationship with the drug׳s preparation and production, and guest solvent molecules contribute significantly to the stabilization of crystalline samples and their molecular conformation. For instance, solvated forms of paclitaxel, an anticancer drug, were found to be useful intermediates. Solvates can also be used as the medicinal forms in clinical treatment. For example, cabazitaxel, another anticancer drug, is a 1:1 acetone solvate (form A) as a commercially available raw material.

This paper focuses on the well-known drug bifendate (dimethyl-4,4′-dimethoxy-5,6,5′,6′-dimethylenedioxybiphenyl-2,2′-dicarboxylate, Fig. 1A), a synthetic hepatoprotective agent derived from schisandrin C12, 13. Bifendate displays a remarkable hepatoprotective ability against chemical liver injury induced by tetrachloride or ethanol as well as other solvents by lowering alanine transaminase (ALT) activity and protecting the liver in patients15, 16, and is widely used for the management of chronic hepatitis in China17, 18, 19. Bifendate exhibits the phenomenon of polymorphism. However, to date, only two polymorphs of bifendate with melting point values of 178–180 °C (form A) and 159–161 °C (form B), respectively, are currently established in the literature. No studies on the crystal structures of bifendate have yet been reported. Consequently, in order to understand both the molecular conformational states in crystal lattices and their packing architectures, to investigate the reasons for polymorphism or solvatomorphism at the atomic scale, and to obtain the standard PXRD patterns for different forms of bifendate, herein a study of the preparation and structural characterization of the bifendate polymorphs and solvatomorphs using SCXRD, thermal analysis and Fourier transform infrared spectrometry (FT-IR) was carried out. Following numerous crystallization attempts the known polymorph, form B, was produced, plus three new solvates, referred to as forms C, D, and E denoting bifendate with tetrahydrofuran, dioxane, and pyridine, respectively, in 1:1 molar ratios. SCXRD confirmed that the molecular conformations were different among distinct forms of bifendate and the simulated PXRD patterns of four forms were calculated followed by the crystallographic data of forms B–E. This provided some scientific basis for the formation of bifendate polymorphs.

Figure 1

Chemical structure and structural overlay for forms B–E of bifendate. Eight molecules of form B are in different colors, solvate C is in red, solvate D is in blue, solvate E is in green. (A) Chemical structure of bifendate; (B) overlay of eight molecules of form B; (C) overlay of solvates C–E.

Materials and methods

Materials

Bifendate (white powder, C20H18O10, MW 418.36, ≥99.6%) was purchased from Hubei Jianyuan Chemical Engineering Inc. (Batch No. 20100707, Hubei, China) and used without further purification. All reagents used for this study were of analytical grade and purchased from Sinopharm Chemical Reagent Company Ltd. (Shanghai, China). Demineralized distilled water was used to prepare buffers.

Preparation of samples

Bifendate crystalline form B and crystallosolvates with tetrahydrofuran (C), dioxane (D) and pyridine (E) were obtained from the filtered acetonitrile–water mixture (3.5:1, v/v), tetrahydrofuran, dioxane–water mixture (3:1, v/v) and pyridine solutions, respectively, by slow evaporation at ambient conditions after 10, 4, 6 and 20 days, respectively.

General techniques

Single-crystal X-ray diffraction (SCXRD)

The X-ray diffraction experiments were carried out on a Rigaku MicroMax-002+ CCD diffractometer with CuKα radiation (λ=1.54178 Å) at 293 K (Rigaku, Americas, the Woodlands, Texas, USA). Absorption correction and integration of the collected data were handled using the CrystalClear software package (Rigaku Americas). The crystal structures were solved by the direct method followed by Fourier syntheses with SHELXS–97 and then refined by full-matrix least-squares procedures using SHELXL–97 on F2 with anisotropic displacement parameters (ADPs) for non-hydrogen atoms. The H atoms were placed in calculated positions and included in the final cycles of refinement in the riding modes. To help data convergence, some restraints were introduced to the refinement of disorders. The PLATON program with the probe radius of 1.2 Å was also employed to calculate the volume of the channels where the solvents resided in an asymmetric unit of solvates.

Calculation of powder X-ray diffractometric patterns

The MERCURY software (version 3.3, Cambridge crystallographic data center, Cambridge, UK), with the conditions of start angle, 3° stop angle, 40° step size, 0.02°, was adopted to calculate the simulated powder X-ray pattern for a CuKα radiation.

Thermal analysis

To investigate the stoichiometry between the guests and the hosts for solvates, thermal gravimetric analysis (TGA) was performed on a Mettler Toledo DSC/TGA 1 calorimeter (Mettler Toledo, Greifensee, Switzerland). Approximately 10 mg of the sample was added to an aluminum oxide cell and heated over the temperature range of 30–500 °C at a constant heating rate of 10 °C/min and purged with a stream of flowing dry nitrogen throughout the experiment at 50 mL/min.

To better understand the thermal events, differential scanning calorimetry (DSC) was also introduced to monitor the desolvation process and determine the melting temperature of corresponding solvates using a Mettler Toledo DSC 1 calorimeter (Mettler Toledo) with the STARe software package (Mettler Toledo) for Microsoft Windows. The samples were weighed into perforated aluminum sample pans (40 μL) and heated at the rate of 10 °C/min from 30 to 200 °C.

FT-IR spectroscopy

FT-IR spectra were recorded with a PerkinElmer Spectrum 400 FT-IR spectrophotometer (PerkinElmer, Waltham, Massachusetts, USA) using an attenuated total reflectance sampling accessory in the scan range of 650–4000 cm−1 with a resolution of 4 cm−1, and then processed via the SPECTRUM software suite (PerkinElmer).

Results

Crystal structure

As internationally recognized, SCXRD can provide precise structural information for polymorphic compositions about molecular conformations, spatial arrangements of molecules and the way that molecules interact. To acquire the relevant information for different forms of bifendate, SCXRD was used to characterize forms B–E. Crystal data and refinement conditions for form B to E are summarized in Table 1. In all four forms, the molecular conformations were comparable, comprising two essentially planar units with the planes linked by the C–C bond of biphenyl group. Additionally, the molecular arrangement mode was also considered together with intermolecular interactive mode, and the simulated PXRD patterns of the four forms were calculated.

Table 1

Crystallographic data and experimental details for four crystalline forms of bifendate.

Parameter	Form B	Form C	Form D	Form E
Empirical formula	C20H18O10	C20H18O10·C4H8O	C20H18O10·C4H8O2	C20H18O10·C5H5N
Molecular weight (Da)	418.36	490.45	506.45	497.44
Crystal system	Triclinic	Orthorhombic	Orthorhombic	Orthorhombic
Space group	P 1¯	P na21	P 212121	P na21
Color/shape	Colorless/Prism	Colorless/Prism	Colorless/Prism	Colorless/Prism
Crystal size (mm3)	0.45×0.60×0.79	0.19×0.27×0.61	0.39×0.78×0.83	0.59×0.76×0.99
Temperature (K)	293(2)	293(2)	293(2)	293(2)
a (Å)	16.305(3)	8.516(2)	8.473(3)	8.555(2)
b (Å)	16.312(3)	21.957(4)	12.529(3)	22.243(4)
c (Å)	29.042(6)	12.429(2)	22.728(5)	12.207(2)
α (°)	97.68(3)	90.00	90.00	90.00
β (°)	97.64(3)	90.00	90.00	90.00
γ (°)	90.03(3)	90.00	90.00	90.00
Volume (Å3)	7,586(3)	2,324.0(8)	2,412.8(12)	2,322.8(8)
Z	16	4	4	4
dcalcd (g/cm3)	1.465	1.402	1.394	1.422
F (0 0 0)	3,488	1,032	1,064	1,040
Absorption coefficient (mm−1)	1.023	0.950	0.965	0.943
Theta range for data collection (°)	1.55–72.65	4.03–72.21	3.89–71.72	3.97–71.97
Reflections collected	107,349	12,118	14,611	14,776
Independent reflections	28,817	3,871	4,521	3,291
Rint	0.0933	0.0487	0.0623	0.0636
R1 [I>2σ(I)]	0.0654	0.0530	0.0454	0.0522
wR2 [I>2σ(I)]	0.1571	0.1445	0.1149	0.1279
R1 (all)	0.0904	0.0579	0.0482	0.0535
wR2 (all)	0.1765	0.1533	0.1190	0.1299
Goodness-of-fit on F2	0.998	1.045	1.059	1.064
Completeness	0.956	0.990	0.986	0.978
CCDC deposition number	1,432,083	1,432,084	1,432,085	1,432,086

SCXRD analysis revealed that form B of bifendate belonged to P space group, triclinic crystal system while solvates C and E crystallized in the orthorhombic space group P na21, P 212121, and P na21, respectively. Form B was a non-solvated polymorph, with the asymmetric unit composed of 8 molecules, and 16 molecules located in the unit cell. The basic molecular skeleton for bifendate was a biphenyl structure, consisting of 2 symmetrical fragments, in which both of the 5- and 6-membered rings were almost planar and coplanar. The guest solvent molecules of solvate C–E were tetrahydrofuran, dioxane and pyridine, respectively, in stoichiometry 1:1 with the host molecule. For the three solvates, the basic unit contained one bifendate molecule and one solvent molecule, which had a higher temperature factor because of the disordered thermal motion of the solvent molecule. The dihedral angles between two planar fragments linked by the C–C bond of the biphenyl group of forms B–E are shown in Table 2, which displayed obvious differences among the four forms.

Table 2

Dihedral angles between the ring planes linked by C–C bond of the biphenyl group.

Form	∠R1–R2 (°)a
B (A)	71.3
C	125.3
D	125.7
E	126.7

R1 denotes the C1/C2/C3/C4/C5/O4/C16/O5/C6 ring; R2 denotes the C7/C12/C11/C10/C9/O7/C17/O6/C8 ring.

The crystal density value of form B was higher than those of forms C to E, illustrating that the bifendate molecules were arranged tightly in three-dimensional space, and the calculated volume (103.7 Å3) of solvent-accessible voids per the unit-cell volume in form B was 1.4%. Comparing that with form B, it was easy to see that all the molecules in the structures of forms C to E were arranged loosely with the corresponding calculated volume of 537.2 Å3 (C), 644.6 Å3 (D), and 534.9 Å3 (E) accounting for approximately 23.1%, 26.7%, and 23.0% of unit-cell volume.

Conformational analysis

Fig. 1B shows a structure overlay for form B of bifendate, with two molecules in different conformational states in an asymmetry unit, and in which the two symmetrical fragments linked by the C–C bond from the biphenyl group were not coplanar, with the corresponding torsion angles being 69.3° and −69.3°, respectively. Apart from that, and despite two methyl ester side chains in each molecule directed to the identical side owing to the single-bonding rotations, two different orientations were shown for the molecules, including the methoxyl fragments. To better describe the molecular conformational state the relative torsion angles in the basic unit for the 4 forms were calculated and are listed in Table 3. Structure overlay for 3 solvates of bifendate are shown in Fig. 1C; it is clear that the host molecules also displayed 2 different molecular conformations for the solvates. Among these, the host molecules in forms C and E adopted the same conformation of the C–C bond of biphenyl group rotating outwards from the paper, with the torsion angle of 120.7(3)° and 122.1(3)°, respectively, while that of solvate D rotated in an inward direction and the torsion angle was −121.3(2)°. In addition, with regard to the single bond rotations, significant differences in the orientation of methyl ester and methoxyl were observed among the three solvates; several of the methoxyl fragments of form D were turned to the identical side in comparison with those in forms C and E, which were pointed to two different sides. Secondly, although a pair of methyl ester groups for each solvate were turned to the identical side, those in solvate D were discovered to adopt to the opposite orientation with respect to those of solvates C and E.

Table 3

Torsion angle values describing conformational states of bifendate molecules in four forms studied.

Form	∠C6–C1–C7–C8	∠C15–O3–C4–C3	∠C1–C2–C13–O1	∠C7–C12–C19–O9	∠C18–O8–C10–C11
B (A/C)	69.3(3)/−69.3(3)	−6.4(3)/5.5(3)	21.8(4)/−25.1(4)	24.4(4)/−23.5(3)	−3.4(3)/7.9(3)
C	120.7(3)	9.5(6)	−25.6(5)	−28.7(5)	−5.8(5)
D	−121.3(2)	−11.8(4)	27.8(3)	27.2(3)	−1.6(4)
E	122.1(3)	−7.2(5)	−30.2(4)	−25.9(4)	9.6(5)

Intermolecular interaction

As depicted in Fig. 2, in the crystal structure of form B, non-classical hydrogen bonds and weak π-π stacking interactions played an important role in forming a supramolecular structure. The host molecules were assembled into a one-dimensional chain along the [0 0 1] direction via weak π-π stacking interactions between the phenyl rings and the oxygenic heterocyclic rings. Together with the non-classical hydrogen bonds, a two-dimensional layer parallel to the (0 0 1) plane can be constructed. When the non-classical hydrogen bonds and van der Waals forces were considered together, a three-dimensional supramolecular structure was completed in the crystal structure. Meanwhile, layered structures of solvates C–E with solvent molecules located in the channels between the layers for solvates C to E are shown in Fig. 3A–C. It was disclosed that the guest solvent molecules were, in turn, arranged in the channel formed by the host molecules, which resulted in a “Z” chain structure parallel to the a axis. The host molecules interacted with the guest molecules through an intercalation mode. Fig. 3D–F display a space-filling representation of bifendate showing the continuous channels along [1 0 0] viewed down crystallographic a axis for the three solvates. As seen in the figure, among these solvates, channels of the forms C and E were similar but different from that of form D.

Figure 2

Stereochemical structure of form B: (A) a view of the one-dimensional catenoid structure of form B along [0 0 1] direction; (B) a view of the two-dimensional layer structure of form B along (0 0 1) plane. Molecules are shown as wireframe models and hydrogen atoms are omitted for clarity.

Figure 3

Upper row: layered structures of solvates C–E with solvent molecules located in the channels between the layers for (A) solvate C: tetrahydrofuran, view down c axis; (B) solvate D: dioxane, view down b axis; (C) solvate E: pyridine, view down c axis. Solvent molecules are represented by spheres corresponding with the van der Waals radii of the respective atoms. Hydrogen atoms are omitted for clarity. Lower row: space-filling representation of bifendate showing the continuous channels along [1 0 0] viewed down crystallographic a axis for (D) solvate C: tetrahydrofuran; (E) solvate D: dioxane and (F) solvate E: pyridine. Bifendate molecules are represented by spheres corresponding with the van der Waals radii of the respective atoms. Solvents are omitted for clarity.

Simulated powder X-ray diffraction analysis

The simulated powder X-ray diffraction patterns for each form of bifendate and solvates C–E after deducting solvents based on the crystallographic data of forms B–E are given in Fig. 4. As shown, the simulated PXRD pattern for forms B–E was unique, and the topologies of PXRD patterns of solvates C to E were similar, especially for forms C and E. Three strong characteristic peaks appeared at the 2θ range of 5–10° and 22–28° for forms C–E, and the values of d were 10.98, 3.67, 3.53 Å (form C); 11.36, 3.77, 3.56 Å (form D); 11.12, 3.63, 3.49 Å (form E) corresponding to the indices of lattice plane (0 2 0), (1 1 3), (1 2 3); (0 0 2), (1 2 4), (1 3 2) and (0 2 0), (1 1 3), (1 2 3), for solvates C, D, and E, respectively. The presence of characteristic peaks coincided well with the results of SCXRD. Moreover, Fig. 4B–D clearly indicates that the varying peaks for solvates C–E after deducting solvents, which can be assigned to the corresponding diffraction peaks of tetrahydrofuran, dioxane, and pyridine, respectively. Thus, the varying peaks for solvates C–E, with the d-spacings of 7.94, 6.69, 5.41, 5.07, 4.89, 4.57 , 4.32, 4.14, 4.07, 3.96, 3.90, 3.78, 3.67, 3.53 Å (form C), 7.94 , 5.65, 5.49, 5.15, 5.04, 4.92, 4.60, 4.42 , 4.20, 4.11, 4.00, 3.95, 3.77, 3.56 Å (form D), and 7.98, 6.68, 6.33, 5.60, 5.35, 4.84, 4.35, 4.12, 3.97, 3.79, 3.70, 3.62, 3.49 Å (form E) correlated to tetrahydrofuran, dioxane, and pyridine, respectively.

Figure 4

Simulated PXRD pattern for form B and comparison of the simulated PXRD patterns for solvates C–E and those after deducting solvents of bifendate. (A) form B; (B) solvate C: tetrahydrofuran; 0.5 solvated; unsolvated; (C) solvate D: dioxane; 0.5 solvated; unsolvated; (D) solvate E: pyridine; 0.5 solvated and unsolvated.

Thermal analysis

To elucidate the thermodynamic behaviors of forms B, C, D, and E for bifendate, differential scanning calorimetry (DSC) and thermal gravimetric analysis (TGA) were employed to assess the thermal events. According to the results of DSC and TGA (Fig. 5), it further confirmed that, except for form B, forms C–E were all solvates. For form B, the DSC heating curve showed a single endothermic peak with a value of 162.82 °C while those of solvates C–E presented double endothermic peaks. The first endothermic peak values in forms C–E were 77.21 °C, 128.27 °C, and 81.42 °C, respectively, and the second were 182.63 °C, 181.65 °C, and 182.11 °C, respectively. The DSC data for forms B–E of bifendate are listed in Table 4. The weight loss of solvents for solvates C, D, and E of bifendate by TGA were 14.57%, 17.13%, and 15.13%, respectively, corresponding to the number of the solvent in each asymmetry unit of solvates of 0.99, 0.98, and 0.96, respectively, which was in agreement with the results of SCXRD.

Figure 5

Thermal analysis profiles of forms B–E. (A) DSC profiles of forms B–E; (B) TGA profiles of forms B–E. Profiles of forms B–E are indicated with black, red, blue, and green, respectively. Top to bottom: form B; solvate C: tetrahydrofuran; solvate D: dioxane; and solvate E: pyridine.

Table 4

DSC data for forms B–E of bifendate.

Form	Tdesolv (°C, peak)	Tfus (°C, peak)
B	–	162.82
C	77.21	182.63
D	128.37	181.65
E	81.42	182.11

IR spectra

With the purpose of further clarifying the distinctions between the different polymorphs, the FT-IR spectra for forms B–E of bifendate were studied and assigned on the basis of a careful comparison between crystalline form B and solvates. As illustrated in Fig. 6, it was noteworthy that the band appearing at 737 cm−1 was the specialized absorption band for form B used to distinguish from form A but not seen in other solvates, which was in accord with the study of DSC. Comparing with form B, it can be found that the intensity of the bending vibration of the C–H bond from the methylene group (1464 cm−1), together with the stretching vibration of the C–O bond from methoxyl group (1262 cm−1) decreased sharply in solvates C and E, while that in solvate D almost disappeared. In addition, the band associated with the δC–C presenting at 1136 cm−1 in form B was shifted considerably toward lower frequency at 1119 cm−1 in form D and was hardly evident in solvates C and E. Lastly, an apparent distinction in the fingerprint region with the range of 1000–850 cm−1 for all the four forms was also observed.

Figure 6

Infrared spectra for forms B to E of bifendate with scanning region of 4000–650 cm−1. Spectra of forms B to E are indicated with solid lines in black, red, blue, and green, respectively. Top to bottom: form B; solvate C: tetrahydrofuran; solvate D: dioxane; and solvate E: pyridine.

Discussion

Based on the results of SCXRD, form B was a non-solvated polymorph and forms C–E were solvates. The different space groups illustrated that the molecular symmetry was different between the non-solvated form and the three solvates. As seen in Table 1, the densities of solvates C–E were lower than that of form B, which indicated the molecules in the structures of forms C–E are arranged loosely, and also demonstrated indirectly the greater thermodynamic stability of form B relative to forms C–E. Furthermore, the larger calculated volumes of solvent-accessible voids per the unit-cell volume in solvates C–E implied that big molecular channels may exist in the stereochemical structure to hold the solvents molecules for the three solvates.

In the structure of form B, the benzene rings of the biphenyl group connect with each other in a skew-skew fashion; this led to two symmetrical fragments linked by the C–C bond of the biphenyl group and resulted in two different conformational molecules presented in an asymmetry unit. For the diversity of the molecular conformation, eight bifendate molecules were assembled into an asymmetry unit for form B. In comparison with form B, two different molecular conformations were also found for solvates C–E; however, because of the adducted solvents, the C–C bond further rotated and the torsion angle increased, resulting in two molecular conformations that were different from that in form B. These adducted solvents had an important influence on the conformations of host molecules.

The interactions of the guest with the host evoked our interest. Theoretically speaking, neither the host molecules nor the guest molecules contained potentially strong hydrogen bond donors or acceptors capable of interacting by classical hydrogen bonds in solvates C–E. In spite of some rigid and planar phenyl rings present in the host molecules, neither the saturated 5-membered ring tetrahydrofuran molecules nor the 6-membered ring dioxane molecules could interact with host molecules by π-π stacking interactions, while the minimum centroid–centroid separation of 4.340(3) Å was too far for the pyridine rings to generate aromatic π-π stacking interactions with the phenyl rings of the host. In fact, as the structure of form B showed, a dihedral angle of approximately 71° between the two planar fragments existed in the host molecule because of the rotation of the C–C bond in the biphenyl skeleton. All the guest solvent molecules for forms C, D, and E involving tetrahydrofuran, dioxane, and pyridine, respectively, were planar with suitable volumes, and could be available to insert into the gap formed by two bifendate molecules. The four host molecules formed a large cavity in which two guest molecules could reside through intercalative interactions, with the host molecules stabilized by weak π-π stacking interactions and non-classical hydrogen bonds to form a molecular channel, where the guest molecules arranged by turns forming a “Z” chain structure as described above. The intercalation mode between guest and host molecules gave rise to further rotation of the C–C bond and resulted in the dihedral angle increasing, with the corresponding value of 125.3° (form C), 125.7° (form D), and 126.7° (form E).

On these grounds, we deduced that to incorporate into the channel formed by the host molecules, guest molecules should meet some structural requirements: on one hand, the guest molecules should be planar so that they could insert into the gap assembled by the host molecules; on the other hand, the expected volumes of the guests should be 100–300 Å3 to maintain the stability of the molecular structures. This explains why we could not generate additional solvates with small molecules such as methanol or ethanol.

The simulated powder X-ray diffraction analysis showed that the simulated PXRD patterns of forms B to E were unique and the patterns of all the solvates after deducting the solvents were also very different. The stated analysis showed that the main differences between the four forms were the solvents, the molecular arrangements and the modes of interaction, which led to the differences in the simulated PXRD patterns of forms B–E. Furthermore, the simulated PXRD patterns of all the solvates after deducting solvents may allow prediction of other new polymorphs for bifendate. In the DSC profiles, the Tfus value of form B was 162.82 °C and was in agreement with its reported melting point, which confirmed that form B was a pure non-solvated polymorph. Among the double endothermic peaks for solvates C–E, the first melting peak can be attributed to that of the solvent while the second was the endothermic peak of polymorph, with the Tfus value of 182.63 °C (form C), 181.65 °C (form D), and 182.11 °C (form E), respectively, consistent with the reported melting point of form A. This indicated that the three solvates could be transformed into form A under suitable temperature conditions because of its greater thermodynamic stability, which was confirmed by desolvation. The differences of the FT-IR spectra for forms B–E of bifendate indicated that the incorporation of solvents could affect the molecular arrangement and intermolecular interaction mode, and resulted in the changes of the infrared absorption spectra.

Conclusions

Our study demonstrates that bifendate has the characteristic of solvatomorphism. Crystalline form B and three new solvates of bifendate with structurally similar solvent molecules were obtained and their crystal structures were characterized. By comparing the crystal structures of form B with the three solvates C–E, it can be concluded that, firstly, weak π-π stacking interactions and non-classical hydrogen bonds made an important contribution to the three-dimensional supramolecular structures for the four forms; secondly, the rotation of the C–C bond of the biphenyl group together with methyl ester and methoxyl fragments led to distinct molecular conformations in crystal lattices, which may result in 8 bifendate molecules per asymmetry unit for form B; thirdly, all the guest solvent molecules interacted with the host molecules through an intercalation pathway arranged by turns in the channel formed by the host molecules; lastly, the incorporation of the guest solvent molecules had an important influence on the host molecular conformation. Consequently, we conclude that the guest molecules should be planar with the expected volumes of approximately 100–300 Å3 to interact with the bifendate molecules. Moreover, all the three solvates were unstable under high temperature conditions. SCXRD analysis illustrated that adducted solvents, the differences in the molecular conformation and the spatial arrangement caused the formation of polymorphs of bifendate. In addition, the unique simulated powder X-ray diffraction patterns were calculated from single crystal experiments of forms B–E and provided a standard for identifying the purity of polymorphs and solvates. Further investigations on this and similar systems are in progress in our laboratory.