Synthesis of Two Novel Copper (II) Complexes as Potential Inhibitors of HIV-1 Protease Enzyme: Experimental and Theoretical Investigations

: In this study, we report the synthesis of two new copper complexes: [Cu(C 11 H 7 O 2 )(SCN)(C 10 H 8 N 2 )], denoted as (C-1), and [Cu(C 11 H 7 O 2 ) (C 12 H 8 N 2 ) Cl] · H 2 O, denoted as (C-2). They are based on 2,2 (cid:48) -bipyridine or 1,10-phenanthroline and 2-hydroxy-1-naphtaldehyde ligands. The obtained complexes were characterized by FT-IR, UV-visible spectroscopy, and single-crystal X-ray diffraction analysis. Molecular docking was employed to predict the binding mode involved in the interaction between the two synthetic copper (II) complexes and HIV-1 protease enzyme. The X-ray structural analysis revealed that the crystal structures of both complexes are mainly stabilized by several intra-and intermolecular hydrogen bonds. The ﬁngerprint plots associated with the Hirshfeld surfaces of both complexes clearly show that H ··· H interactions provide the largest contributions. According to the docking results, the synthesized complexes exhibit promising features which enable them to be bound to the HIV-protease enzyme.


Introduction
In recent years, there have been several reports dedicated to the development of transition metal complexes that are based on the 2-hydroxy-1-naphthaldehyde C 11 H 8 O 2 ligand [1].Indeed, 2-hydroxy-1-naphthaldehyde ligands can be used as a starting block for the synthesis of new ligands [2,3], or as a potential chelating agent for ions; however, only a few studies were focused on the synthesis of copper (II) complexes of bidentate ligands [4,5].
Moreover, the 2,2 -bipyridine C 10 H 8 N 2 and 1,10-phenanthroline C 12 H 8 N 2 bidentate chelating ligands are widely used as complexing moieties to form stable coordination complexes with miscellaneous transition metal ions [6].Transition metal complexes that have these ligands are recognized in the literature as potent antibacterial and antimicrobial agents [7][8][9].It was also reported that copper complexes are good candidates for the design of novel types of protease enzyme inhibitors; here, the main goal is to treat the acquired immunodeficiency syndrome (AIDS) infection induced by the human immunodeficiency virus (HIV).Since the protease of the human immunodeficiency virus is key to the replication of the virus, the enzyme is therefore considered an attractive target for antiviral therapy; thus, the treatment strategy for AIDS involves the inhibition of HIV-1 protease [10].To date, real progress in terms of therapy for AIDS has been achieved by using peptide-based therapeutics; however, the toxicity of peptides reduces the effectiveness of this treatment [11,12].The copper complexes, which contain non-peptide ligands that have therapeutic potential, could be an alternative methodology for the treatment of AIDS-induced diseases.
In this contribution, we describe the synthesis of two new complexes, namely, the [(2,2 -bipyridine-K 2 N,N ) (1-formyl-2-naphtholato-K 2 O,O ) (Thiocyanate-KN) copper(II)], denoted as C-1, and [(1-formyl-2-naphtholato-K 2 O,O ) (1,10-phenanthroline-K 2 N,N ) chloride copper], denoted as C-2.FT-IR, UV-visible spectroscopy, and single-crystal X-ray diffraction analysis were used to characterize and analyze the structure of the two synthesized complexes.DFT calculations using the B3LYP [13][14][15][16] functional were conducted on the molecular structures of C-1 and C-2, and the findings were compared with the X-ray crystal data.The docking of C-1 and C-2 into the HIV-1 protease active site was assessed through the calculation of binding energy with the Autodock vina program [17], and the obtained results were compared with those of the (dimethanol)bis[N-(4-methyl-2-pyridyl)-2,3,6-trimethoxybenzamide] copper (II) dimethanol diperchlorate complex as a reference (complex-ref), which was found to be a non-peptide HIV-1 protease inhibitor with an IC50 of 1.5 µM [18].It is well known that the most effective inhibitor is the one that forms a bond deep in the hydrophobic cavity of the HIV-1 protease active site, and is positioned in the axes of an octahedron with the four subsites S1/S2 and S 1/S 2 forming the equatorial plane.The binding pocket of this enzyme includes the Asp25 (25 )-Thr26 (26 )-Gly27 (27) catalytic triad and flap regions, which are most likely involved in the binding process of substrates or inhibitors [19,20].

Materials and Instruments
All reagents and solvents were purchased from commercial sources (Sigma-Aldrich, St. Louis, MI, USA), and they were used in the form in which they were received.
Infrared spectra were recorded in the range of 4000-400 cm −1 on a FT-IR Bruker ATR Vertex 70 Spectrometer.UV-visible spectra were recorded on a UV-visible Spectrophotometer Optizen 1220.The absorbance of ligands and their metals complexes were measured in DMSO in a concentration of 10 −5 mol•L −1 , in a wavelength range of 280-440 nm.
Crystals suitable for single crystal X-ray diffraction were selected, and the lattice parameters were determined using an APEXII Bruker diffractometer.

General
Procedure for the Synthesis of C-1 and C-2 For the synthesis of C-1, Cu(CH 3 COO) 2 .H 2 O (0.25 mmol) and 2-hydroxy-1-naphthaldehyde (0.125 mmol) were added and stirred into a DMSO solution containing 0.25 mmol of 2,2 -bipyridine and 0.125 mmol of KSCN.The mixture was stirred for 1 h at T = 25 • C (Scheme 1a).After several days of slow evaporation, light green crystals appeared.This complex is stable in the air and soluble only in DMSO.Yield: 61%.IR (cm

X-ray Diffraction
Data collection was performed using MoKα radiation with a BRUKER APEX2 diffractometer [21].Cell refinement and data reduction were performed through the SAINT program.SHELXS and SHELXL [22] were used for structure solution and refinement.Structure refinement and crystal data details are reported in Table 1.All hydrogen atoms were refined in geometrically idealized positions, with C-H = 0.93 Å for aromatic rings (2,2′-bipyridine, 1,10-phenanthroline), and they were allowed to ride on their parent atoms, with Uiso(H) = 1.2 Ueq(C).The molecular structures of the two complexes are illustrated in Figure 1.

X-ray Diffraction
Data collection was performed using MoKα radiation with a BRUKER APEX2 diffractometer [21].Cell refinement and data reduction were performed through the SAINT program.SHELXS and SHELXL [22] were used for structure solution and refinement.Structure refinement and crystal data details are reported in Table 1.All hydrogen atoms were refined in geometrically idealized positions, with C-H = 0.93 Å for aromatic rings (2,2 -bipyridine, 1,10-phenanthroline), and they were allowed to ride on their parent atoms, with Uiso(H) = 1.2 Ueq(C).The molecular structures of the two complexes are illustrated in Figure 1.

Computational Methodology
The structural optimization of C-1 and C-2 was performed using the framework of the density functional theory (DFT).The B3LYP hybrid functional was used with a 6-31 + G(d) basis set for C, H, N, Cl, O, S, and SDD for Cu.The DFT computations were conducted using the Gaussian 09W software package [23].The PDB file of the ligand was generated using OpenBabel, version 2.4.1 [24,25].The HIV-1 protease-ligand interaction properties were assessed by performing docking studies using AutoDock vina software.The threedimensional structure of HIV-1 protease (PDB ID: 3TLH) was retrieved from the Protein Data Bank (http://www.rcsb.org/pdb(accessed on 8 March 2022)), and it was used to prepare the target site for docking calculation.The water molecules and ligand were removed from the protein, and the missing hydrogen atoms and charge were corrected.The search grid of the HIV-1 protease was identified as center_x: −6.118, center_y: 23.382, center_z: 36.63 with dimensions size_x: 62, size_y: 82, and size_z: 76.The best-scoring conformation was visualized and analyzed using Autodock vina tools.

X-ray Diffraction
Data collection was performed using MoKα radiation with a BRUKER APEX2 diffractometer [21].Cell refinement and data reduction were performed through the SAINT program.SHELXS and SHELXL [22] were used for structure solution and refinement.Structure refinement and crystal data details are reported in Table 1.All hydrogen atoms were refined in geometrically idealized positions, with C-H = 0.93 Å for aromatic rings (2,2′-bipyridine, 1,10-phenanthroline), and they were allowed to ride on their parent atoms, with Uiso(H) = 1.2 Ueq(C).The molecular structures of the two complexes are illustrated in Figure 1.

Results and Discussions
The single-crystal X-ray diffraction shows that the [Cu (C 11 H 7 O 2 ) (SCN) (C 10 H 8 N 2 )] complex crystallizes in the monoclinic space group P2 1 /n with Z = 4.The asymmetric unit of this structure is formed by a Cu(II) that coordinates with a thiocyanate ligand in an apical position, together with a nitrogen atom (N3).The coordination sphere of Cu(II) is completed by one 2,2 -bipyridine molecule acting as a bidentate ligand (through the two nitrogen atoms (N1, N2)), and one 2-hydroxy-1-naphthaldehyde (through the two oxygen atoms (O1, O2)).The Cu(II) atom has five coordinates, and it displays a square-pyramidal coordination geometry, as quantified by the value of the structural index [26]  of this structure is formed by a Cu(II) that coordinates with a thiocyanate ligand in an apical position, together with a nitrogen atom (N3).The coordination sphere of Cu(II) is completed by one 2,2′-bipyridine molecule acting as a bidentate ligand (through the two nitrogen atoms (N1, N2)), and one 2-hydroxy-1-naphthaldehyde (through the two oxygen atoms (O1, O2)).The Cu(II) atom has five coordinates, and it displays a square-pyramidal coordination geometry, as quantified by the value of the structural index [26] S1); this is in accordance with the data in the literature [27].These parameters were also calculated in the gas phase using the B3LYP functional, and the results indicate that the calculated bond lengths in the copper coordination sphere agree with the corresponding experimental values within a maximum deviation of 2.3% (Table S1).The bond angles exhibit a maximum deviation of 10%.This discrepancy between the experimental and theoretical values may be due to environmental effects; indeed, the theoretical values were obtained from gas-phase DFT calculations, whereas the experimental values were taken from crystallographic data.
The crystalline structure of C-1 consists of four hydrogen bond networks (Table 2).Two intermolecular hydrogen bonds are observed between H7 of C11 with S1 (½ − x, ½ + y, ½ − z), and H14 of C20 with S1(x + 1, y, z).Moreover, two intramolecular hydrogen bonds of the C-H•••O type are established between 2,2′-bipyridine and 2-hydroxy-1-naphthaldehyde, subsequently generating S(5) ring motifs [28,29]; thus, the presence of hydrogen-bond networks points to the significant role of intra and intermolecular hydrogen bonding in terms of the stabilization of complex structures.Indeed, the networks facilitate  S1); this is in accordance with the data in the literature [27].These parameters were also calculated in the gas phase using the B3LYP functional, and the results indicate that the calculated bond lengths in the copper coordination sphere agree with the corresponding experimental values within a maximum deviation of 2.3% (Table S1).The bond angles exhibit a maximum deviation of 10%.This discrepancy between the experimental and theoretical values may be due to environmental effects; indeed, the theoretical values were obtained from gas-phase DFT calculations, whereas the experimental values were taken from crystallographic data.
The crystalline structure of C-1 consists of four hydrogen bond networks (Table 2).Two intermolecular hydrogen bonds are observed between H7 of C11 with S1 ( of this structure is formed by a Cu(II) that coordinates with a thiocyanate ligand in an apical position, together with a nitrogen atom (N3).The coordination sphere of Cu(II) is completed by one 2,2′-bipyridine molecule acting as a bidentate ligand (through the two nitrogen atoms (N1, N2)), and one 2-hydroxy-1-naphthaldehyde (through the two oxygen atoms (O1, O2)).The Cu(II) atom has five coordinates, and it displays a square-pyramidal coordination geometry, as quantified by the value of the structural index [26] S1); this is in accordance with the data in the literature [27].These parameters were also calculated in the gas phase using the B3LYP functional, and the results indicate that the calculated bond lengths in the copper coordination sphere agree with the corresponding experimental values within a maximum deviation of 2.3% (Table S1).The bond angles exhibit a maximum deviation of 10%.This discrepancy between the experimental and theoretical values may be due to environmental effects; indeed, the theoretical values were obtained from gas-phase DFT calculations, whereas the experimental values were taken from crystallographic data.
The crystalline structure of C-1 consists of four hydrogen bond networks (Table 2).Two intermolecular hydrogen bonds are observed between H7 of C11 with S1 (½ − x, ½ + y, ½ − z), and H14 of C20 with S1(x + 1, y, z).Moreover, two intramolecular hydrogen bonds of the C-H•••O type are established between 2,2′-bipyridine and 2-hydroxy-1-naphthaldehyde, subsequently generating S(5) ring motifs [28,29]; thus, the presence of hydrogen-bond networks points to the significant role of intra and intermolecular hydrogen bonding in terms of the stabilization of complex structures.Indeed, the networks facilitate  ance with the data in the literature [27].These s phase using the B3LYP functional, and the re-ngths in the copper coordination sphere agree es within a maximum deviation of 2.3% (Table deviation of 10%.This discrepancy between the e due to environmental effects; indeed, the the-ase DFT calculations, whereas the experimental ata.

½ + y, z).
Moreover, two intramolecular hydrogen between 2,2′-bipyridine and 2-hydroxy-1-naph-ring motifs [28,29]; thus, the presence of hydro-ant role of intra and intermolecular hydrogen plex structures.Indeed, the networks facilitate  he apical Cu1-N3 bond 2.188 Å is longer than those of accordance with the data in the literature [27].These the gas phase using the B3LYP functional, and the re-ond lengths in the copper coordination sphere agree tal values within a maximum deviation of 2.3% (Table imum deviation of 10%.This discrepancy between the s may be due to environmental effects; indeed, the the-gas-phase DFT calculations, whereas the experimental aphic data.
This complex crystallizes in the triclinic system, space group P-1, with Z = 2.The asymmetric unit in this structure is composed of a copper (II) cation that is coordinated with the nitrogen atoms of the 1,10-phenanthroline, two oxygen atoms of the 2-hydroxy-1naphtaldehyde in the equatorial plane, and a chlorine anion in the axial position.We also note the presence of a water molecule.

FT-IR and UV-Visible Absorption Spectra
Figure 5a shows the FT-IR spectra of the two complexes.Overall, both complexes show similar IR bands.The strong absorptions at 1604 cm −1 and 1119 cm −1 , correspond to the C=O and C-O bonds of the aldehyde, respectively.The C-1 complex exhibits a broad band in the region of 3100-2800 cm −1 , which is attributed to the υ(C-H) of the aromatic ring.The bands appearing at 1562 cm −1 and 1014 cm −1 are due to the υ (C=N) and υ(C-N) of the 2-2 -bipyridine, respectively.The υ (C=N) band of the thieno ligand is observed at 2063 cm −1 .In the FT-IR spectrum of C-2, the bands appearing at 3100 cm −1 and 2855 cm −1 correspond to the aromatic υ(C-H).The spectrum also shows bands at 1434, 1274, and 730 cm −1 , which are assigned to the coordinated phenanthroline molecule [31,32].
The two complexes exhibit a d-d transition band at 624-640 cm −1 , and a band at 821-833 cm −1 , which is typical for square pyramidal geometry [30].
Figure 5b shows the UV-visible spectra in the range of 280-440 nm for C-1 and C-2.The absorption peaks of these two complexes are very similar.Nevertheless, those of C-2 are very weak.In the UV region, the two complexes show peaks near 290-300 nm due to the π → π* aromatic ring.Moreover, we noticed two intense peaks, or a shoulder, in the (305 → 315 nm) (I)-(310 → 325 nm) (II) region, which can be assigned to the charge transfer between the ligand and the metal.

FT-IR and UV-Visible Absorption Spectra
Figure 5a shows the FT-IR spectra of the two complexes.Overall, both complexes show similar IR bands.The strong absorptions at 1604 cm −1 and 1119 cm −1 , correspond to the C=O and C-O bonds of the aldehyde, respectively.The C-1 complex exhibits a broad band in the region of 3100-2800 cm −1 , which is attributed to the υ(C-H) of the aromatic ring.The bands appearing at 1562 cm −1 and 1014 cm −1 are due to the υ (C=N) and υ(C-N) of the 2-2′-bipyridine, respectively.The υ (C=N) band of the thieno ligand is observed at 2063 cm −1 .In the FT-IR spectrum of C-2, the bands appearing at 3100 cm −1 and 2855 cm −1 correspond to the aromatic υ(C-H).The spectrum also shows bands at 1434, 1274, and 730 cm −1 , which are assigned to the coordinated phenanthroline molecule [31,32].
The two complexes exhibit a d-d transition band at 624-640 cm −1 , and a band at 821-833 cm −1 , which is typical for square pyramidal geometry [30].Figure 5b shows the UV-visible spectra in the range of 280-440 nm for C-1 and C-2.The absorption peaks of these two complexes are very similar.Nevertheless, those of C-2 are very weak.In the UV region, the two complexes show peaks near 290-300 nm due to the π → π* aromatic ring.Moreover, we noticed two intense peaks, or a shoulder, in the (305→315 nm) (I)-(310→325 nm) (II) region, which can be assigned to the charge transfer between the ligand and the metal.No d-d bands were observed in the spectra of the two complexes.These bands should be observed in the region of 500-600 nm.Their absence is due to the low concentration (10 −5 mol•L −1 ) of the solutions of the complexes.

Hirshfeld Surface Analysis
The intermolecular interactions and the packing modes that occur in crystalline structures can be studied quantitatively using Hirshfeld surface analysis [33].Hirshfeld surfaces (HS) and fingerprint plots [34] were carried out using Crystal Explorer software [35].
The HS of the C-1 complex was mapped over d norm in the range of −0.1832 to +1.3833 a.u., as shown in Figure 6a.Color coding is used for each specific region to identify the nature of the occurring intermolecular interactions, wherein contacts that are shorter, equal to, and longer than the sum of the van der Waals radii, correspond to the red, white and blue areas, respectively.The major interactions are caused by H 6a).
The shape-index surface (Figure 6b) indicates the presence of blue and red triangles that are characteristic of π•••π interactions between neighboring molecules, whereas the flat surface patches in the curvedness plots (Figure 6c) reveal planar stacking.
The 2-D fingerprint plots for the C-1 complex show the relative contributions of H

Inhibitory Activity of the C-1 and C-2 Complexes against the HIV-1 Protease Enzyme
To provide insight into the binding mode of C-1 and C-2, a molecular docking investigation was carried out.As shown in Figure 8 (upper panel), the tested complexes occupied a common cavity in the receptor, suggesting that C-1 and C-2 block access to the active site of the HIV-1 protease.They are located in the hydrophobic pocket of HIV-1 protease, and they are surrounded by the Val32, Pro81, and Pro79 residues of the active site, and Gly48, Gly49, Ile50, Ile84, Ile47, and Leu84 of the subsites S1/S1 , S2/S2 , and S3/S3 , wherein a strong hydrophobic interaction was observed, as illustrated in Figure 8 (lower panel).The ligand-receptor complex indicates that most of the ligand is positioned in the flap region of the enzyme, as represented by the Ile50, Gly51, Gly52, and Gly49 residues.No hydrogen bonds were found in the docked structures.

Figure 2 .
Figure 2. The copper square-pyramidal environment (a), and the dihedral angle between the two ligands (b).The dihedral angle between the two planes that contain ligands is 3.65°.The Cu-N and Cu-O distances in the basal plane are [(Cu1-N1 = 2.005 Å, Cu1-N2 = 2.001 Å, Cu1-O1 = 1.908Å, Cu1-O2 = 1.943Å], and the apical Cu1-N3 bond 2.188 Å is longer than those of the basal plane (TableS1); this is in accordance with the data in the literature[27].These parameters were also calculated in the gas phase using the B3LYP functional, and the results indicate that the calculated bond lengths in the copper coordination sphere agree with the corresponding experimental values within a maximum deviation of 2.3% (TableS1).The bond angles exhibit a maximum deviation of 10%.This discrepancy between the experimental and theoretical values may be due to environmental effects; indeed, the theoretical values were obtained from gas-phase DFT calculations, whereas the experimental values were taken from crystallographic data.The crystalline structure of C-1 consists of four hydrogen bond networks (Table2).Two intermolecular hydrogen bonds are observed between H7 of C11 with S1 (½ − x, ½ + y, ½ − z), and H14 of C20 with S1(x + 1, y, z).Moreover, two intramolecular hydrogen bonds of the C-H•••O type are established between 2,2′-bipyridine and 2-hydroxy-1-naphthaldehyde, subsequently generating S(5) ring motifs[28,29]; thus, the presence of hydrogen-bond networks points to the significant role of intra and intermolecular hydrogen bonding in terms of the stabilization of complex structures.Indeed, the networks facilitate

Figure 2 .
Figure 2. The copper square-pyramidal environment (a), and the dihedral angle between the two ligands (b).The dihedral angle between the two planes that contain ligands is 3.65 • .The Cu-N and Cu-O distances in the basal plane are [(Cu1-N1 = 2.005 Å, Cu1-N2 = 2.001 Å, Cu1-O1 = 1.908Å, Cu1-O2 = 1.943Å], and the apical Cu1-N3 bond 2.188 Å is longer than those of the basal plane (TableS1); this is in accordance with the data in the literature[27].These parameters were also calculated in the gas phase using the B3LYP functional, and the results indicate that the calculated bond lengths in the copper coordination sphere agree with the corresponding experimental values within a maximum deviation of 2.3% (TableS1).The bond angles exhibit a maximum deviation of 10%.This discrepancy between the experimental and theoretical values may be due to environmental effects; indeed, the theoretical values were obtained from gas-phase DFT calculations, whereas the experimental values were taken from crystallographic data.The crystalline structure of C-1 consists of four hydrogen bond networks (Table2).Two intermolecular hydrogen bonds are observed between H7 of C11 with S1 ( τ5 = 0.019 (τ5 = (β − α)/60, where β = N1-Cu1-O2 = 164.11°and α = N2-Cu1-O1 = 162.97°.(Perfect coordination geometries of square pyramidal and trigonal bipyramidal arrangements are associated with the τ5 values of zero and unity, respectively).The basal plane is formed by (O1, O2, N1, N2) and N3 atoms in the axial site (Figure2a).

Figure 2 .
Figure 2. The copper square-pyramidal environment (a), and the dihedral angle between the two ligands (b).The dihedral angle between the two planes that contain ligands is 3.65°.The Cu-N and Cu-O distances in the basal plane are [(Cu1-N1 = 2.005 Å, Cu1-N2 = 2.001 Å, Cu1-O1 = 1.908Å, Cu1-O2 = 1.943Å], and the apical Cu1-N3 bond 2.188 Å is longer than those of the basal plane (TableS1); this is in accordance with the data in the literature[27].These parameters were also calculated in the gas phase using the B3LYP functional, and the results indicate that the calculated bond lengths in the copper coordination sphere agree with the corresponding experimental values within a maximum deviation of 2.3% (TableS1).The bond angles exhibit a maximum deviation of 10%.This discrepancy between the experimental and theoretical values may be due to environmental effects; indeed, the theoretical values were obtained from gas-phase DFT calculations, whereas the experimental values were taken from crystallographic data.The crystalline structure of C-1 consists of four hydrogen bond networks (Table2).Two intermolecular hydrogen bonds are observed between H7 of C11 with S1 (½ − x, ½ + y, ½ − z), and H14 of C20 with S1(x + 1, y, z).Moreover, two intramolecular hydrogen bonds of the C-H•••O type are established between 2,2′-bipyridine and 2-hydroxy-1-naphthaldehyde, subsequently generating S(5) ring motifs[28,29]; thus, the presence of hydrogen-bond networks points to the significant role of intra and intermolecular hydrogen bonding in terms of the stabilization of complex structures.Indeed, the networks facilitate at coordinates with a thiocyanate ligand in an tom (N3).The coordination sphere of Cu(II) is e acting as a bidentate ligand (through the two xy-1-naphthaldehyde (through the two oxygen coordinates, and it displays a square-pyramidal e value of the structural index [26] τ5 = 0.019 (τ5 1° and α = N2-Cu1-O1 = 162.97°.(Perfect coor-and trigonal bipyramidal arrangements are as-ity, respectively).The basal plane is formed by l site (Figure 2a).

Figure 5 .
Figure 5.The FT-IR Spectrum (a) and the overlapping UV-visible spectrum in the DMSO (b) of complexes C-1 and C-2.
to the Hirshfeld surface; these are illustrated in Figure6d.Contact between H•••H comprise the predominant form of interactions (33.7%).The HS of the C-2 complex is depicted in Figure7a.The dominant interactions between (water) O-H, C-H, and Cl atoms are represented by the bright-red areas, which are marked as 1, 2, and 3.The light-red spots denoted as 4, 5, 6 and 7 are due to C-H•••O and C-H•••Cl interactions.The π•••π interactions in the C-2 complex are revealed by the characteristic blue and red triangles on the shape-index surface (Figure7b), whereas planar stacking is identified by the curvedness plots (Figure7c) that exhibit characteristic flat surface patches.The relative contributions of the H•••H, H•••C/C•••H, H•••O/O•••H, H•••Cl/Cl•••H and H•••N/N•••H contactsto the Hirshfeld surface in the C-2 complex are determined using the 2-D fingerprint plots, as illustrated in Figure7d.As in the C-1complex, the fingerprint plots of the C-2 complex show that the H•••H contacts provide the most significant contribution to the Hirshfeld surface (45.1%).

Figure 6 .
Figure 6.A view of the Hirshfeld surface for (I), which is mapped over (a) d norm in the range of −0.1832 to +1.3833 arbitrary units.(b) Shape-index, (c) curvedness, and (d) 2-D fingerprint plots for the C-1 complex, which reveal the contributions of all contacts.

Figure 6 .
Figure 6.A view of the Hirshfeld surface for (I), which is mapped over (a) dnorm in the range −0.1832 to +1.3833 arbitrary units.(b) Shape-index, (c) curvedness, and (d) 2-D fingerprint plots f the C-1 complex, which reveal the contributions of all contacts.

Figure 7 .
Figure 7.A view of the Hirshfeld surface for (II), which is mapped over (a) d norm in the range of −0.4175 to +1.3253 arbitrary units.(b) Shape-index, (c) curvedness, and (d) 2-D fingerprint plots for the C-2 complex, which reveal the contributions of all contacts.