Abstract

The interaction between nickel (Ni2+), copper (Cu2+), and zinc (Zn2+) ions and 1-methylimidazole has been studied by exploring the geometries of eleven crystal structures in the Cambridge Structural Database (CSD). The coordination behavior of the respective ions was further investigated by means of density functional theory (DFT) methods. The gas-phase complexes were fully optimized using B3LYP/GENECP functionals with 6-31G and LANL2DZ basis sets. The Ni2+ and Cu2+ complexes show distorted tetrahedral geometries around the central ions, with Zn2+ being a perfect tetrahedron. Natural bond orbital (NBO) analysis and natural population analysis (NPA) show substantial reduction in the formal charge on the respective ions. The interaction between metal d-orbitals (donor) and ligand orbitals (acceptor) was also explored using second-order perturbation of the Fock matrix. These interactions followed the order Ni2+ > Cu2+ > Zn2+ with Zn2+ having the least interaction with the ligand orbitals. Examination of the frontier orbitals shows the stability of the complexes in the order Ni2+ > Cu2+ < Zn2+ which is consistent with the Irving–Williams series.

1. Introduction

Metalloproteins play important roles in the structure and physiology of cells. They account for nearly half of all proteins in nature [1]. Some of their cellular functions include water oxidation, photosynthesis, respiration, molecular oxygen reduction, and zinc fingers [1, 2]. The binding stability of divalent transition metals in metalloproteins has been largely studied in terms of semiempirical and qualitative theories such as the hard and soft acid and base principles of Parr and Pearson and the Irving–Williams series of stability constants [3]. Yet the fundamental role of metal ions in the structure and function of metalloproteins is still a matter of ongoing research [4, 5]. The amino acid histidine is a common ligand in metalloproteins [6]. The imidazole ring of histidine is an important five-membered heterocycle that is widely present in natural products and most synthetic molecules. Because of its unique electronic structure, imidazole-based compounds are widely used as anticancer, antifungal, antibacterial, anti-inflammatory, antihistaminic, and other medicinal agents [710]. It readily binds with a variety of enzymes and receptors in biological systems through hydrogen bonds, coordinate covalent bonds, ion-dipole interactions, π-π interactions, and van der Waals forces, thereby exhibiting broad bioactivities [7].

Imidazole is largely considered an N-donor in most of its coordination complexes [11]. This interaction largely affects the formal charge of the coordinating electrophile. According to Hasegawa et al., the total charge of 4-methylimidazole increased upon binding to Zn2+, indicating that some negative charge is transferred to the Zn2+ in the complex. This charge delocalization has a tendency of reducing the formal charge on Zn2+ to 1.70. Similar atomic charges were observed in the nitrogen atoms of the 4-methylimidazole ligand [12]. This charge delocalization could greatly affect the reactivity and catalytic activity of imidazole-based transition metal complexes.

A study by Rulı́šekr and Vondrášek on the coordination geometries of selected metal ions in metalloproteins with data from the Protein Data Bank suggests that Ni2+, Cu2+, and Zn2+ tend to bind in tetrahedral coordination modes although there are other geometries such as octahedral and square planar [3]. A search in the Cambridge Structural Database version 5.38 (November 2016) plus one update revealed one Ni complex bearing 1-methylimidazole ligand [13]. In this complex, the nickel ion is octahedrally coordinated with the trien ligand bonded through its four nitrogen atoms with the remaining two coordination sites occupied by 1-methylimidazole ligands. Ten copper complexes bearing 1-methylimidazole ligands were also found. These include molecules with the following reference codes: BEJGUS [14], CAHJAW [15], and CUSHON [16] with square-planar coordination around the copper ion through two carboxylic oxygen atoms and two 1-methylimidazole nitrogen atoms; MACCUA10 [17] of square-pyramidal geometry with the four nitrogen atoms of the macrocycle forming the basal plane and 1-methylimidazole nitrogen at the apex; GALLAG [18] of square-planar coordination through four 1-methylimidazole ligands; GALLOU [18], a tetragonally distorted octahedral complex with two water molecules trans to each other and four 1-methylimidazole ligands in the square plane; KAYPEF [19] (chloro(glycinato)(1-methylimidazole)copper(II) complex) with a square-pyramidal structure having four close ligating atoms (N2OCl) and an axial chlorine ligand; CEZLOI [20], a Cu+ complex with S4 site symmetry in which the Cu ion is tetrahedrally coordinated with four 1-methylimidazole ligands; and finally, the complexes GALLEK and GALLIO with tetrahedral geometries in which the copper ions are linked to a central oxygen atom with μ2-Cl atoms above each edge of the tetrahedron and terminal 1-methylimidazole ligand. However, there was no zinc complex bearing 1-methylimidazole ligand at the time of preparing this manuscript. Considering the importance of nickel, copper, and zinc ions in the structure and function of metalloproteins, this article investigates the effect of bonding on the energies of the d-orbitals of Ni2+, Cu2+, and Zn2+ in tetrahedral ligand fields. It also reports the effect on the formal charges of the ions upon coordination to the 1-methylimidazole (1-MeIm) ligand.

2. Methodology

2.1. CSD Analysis

The ligand and transition metal complexes were analyzed with version 5.38 of the CSD (November 2016) plus one update. The CSD program ConQuest Version 1.19 was used to perform substructure searches of Ni, Cu, and Zn complexes bearing 1-methylimidazole ligands. The accepted entries met the following criteria: 3D coordinates determined; crystallographic R factor ≤ 0.05; no disorder in the crystal structure; no errors in the structure; no polymeric bonding; no ions; and no powder structures and only organometallic structures (according to standard CSD definitions) [21]. The search revealed eleven crystal structures: one nickel and ten copper complexes. No zinc complex was found in the database. The geometry of N-methylimidazole in a novel 4- and 5-coordinated silicon complex (reference code: GAGXER) [22] was used as the reference. The 3D structures of all the complexes were visualized using version 3.9 of the CSD program Mercury [23].

2.2. Computational Studies

All input files were prepared using the GaussView 5.0.8 molecular structure viewer [24]. The Gaussian 09 program was used to perform all the computations. All structures were optimized with the density functional theory (DFT) using B3LYP three-parameter hybrid functionals with no constraints in the respective geometries. Effective core potentials (ECPs) were used to represent the valence electrons of Ni2+, Cu2+, and Zn2+ and the basis set of double-ζ quality associated with pseudopotentials known as LANL2DZ [25]. Because of the size of the complexes and the accompanying computational cost, the 6-31G basis set was used for all other atoms. Similar basis sets have been used to study the geometries and molecular orbitals of transition metal complexes [26]. Frequency analyses were employed to confirm that the optimized geometries were at stationary points corresponding to local minima with no imaginary frequencies [21]. The optimized geometries with the accompanying numbering scheme and the Cartesian atomic coordinates are shown in Figure 1.

3. Results and Discussion

3.1. Geometry Optimization

Figure 1 shows the gas-phase optimized geometries and the selected atomic numbering scheme of 1-methylimidazole and the respective complexes. It also indicates the Cartesian atomic coordinates. Selected bond lengths of the optimized structures are shown in Table 1. A search in the version 5.38 of the CSD (November 2016) showed that there was no crystal structure of 1-methylimidazole in the database. However, the geometry of N-methylimidazole in a novel 4- and 5-coordinated silicon complex (reference code: GAGXER) [22] was used as the reference. The bond lengths did not differ significantly from that of the optimized structure. The slight differences could be attributed to the fact that the optimization was carried out on a single molecule in the gas phase with no intermolecular interactions, whereas in the crystal structure, there are lattice interactions which affect the bond parameters [21].

Generally, the intramolecular bond lengths in 1-MeIm were not affected by the coordination. Notable exceptions include the C1-N4, C2-N8, and C3-N8 bonds which recorded slight decreases in bond lengths as a result of coordination to the metal ion. The metal-ligand bond lengths (M-N) ranged between 2.018 and 2.036 Å in the Ni-1MeIm complex and between 2.017 and 2.038 Å in the Cu-1MeIm analogue but remained the same (2.085 Å) in the zinc complex. The differences in the Ni-N and Cu-N bond lengths could be attributed to Jahn–Teller distortions for the (e)4(t2)4 and (e)4(t2)5 electronic configurations in the respective Ni-1MeIm and Cu-1MeIm complexes [27]. Zn2+ has fully occupied d-orbitals and interacts weakly with ligand orbitals and maintains a spherical isotropy. The C3-H7 bond length did not show much variation as a result of the bonding interaction of the ligand with the metal ions. The observed differences in bond lengths were further investigated by a series of density functional theory (DFT) calculations.

3.2. Natural Bond Orbital (NBO) Analysis and Charge Distribution

The effect of coordination and electron distribution on the ligands and metal d-orbitals was assessed by the NBO and second-order perturbation of the Fock matrix. This provides details about the type of hybridization, nature of bonding, and strength of interaction between the metal ion and ligand [28]. Table 2 shows the occupancy and electron density (ED) of molecular orbitals in 1-MeIm. The C1-N4 σ-bond in the uncoordinated ligand has an electron density of 1.980 and an occupancy of 41.4% C1 and 58.6% N4. Upon coordination, the electron density reduced to 0.990 in both Ni-1MeIm and Cu-1MeIm with the occupancies also changing to 38.1% and 61.9% for C1 and N4, respectively.

The reduced electron density in the C1-N4 σ-bond could have been delocalized onto the metal ions [29] to reduce their formal charges from +2 to −0.2985 and +0.2132, respectively, as shown in Table 3. Despite the fact that occupancy changed to 38.8% C1 and 61.7% N4 in Zn-1MeIm, the ED was not affected. Similar effects were observed in the C3-N4 σ-bond. The C3-N4 π-bond observed in the free ligand was, however, nonexistent in the Ni2+ and Cu2+ complexes. These π-electrons could have contributed to reducing the strength of the formal charges on the respective metal ions as observed in Table 3. Generally, although there were substantial changes in the occupancy of the bonds in the 1-MeIm upon coordination to the Zn2+ ion, the ED on the ligand remained unchanged as compared to the free ligand.

From Table 3, the calculated formal charges on the central metal ions have been reduced to −0.2985, +0.2132, and +1.3850, respectively. This reduction in the ionic charge is higher in Ni-1MeIm and Cu-1MeIm than in Zn-1MeIm. As shown in Table 2, there is higher delocalization of ligand electrons onto Ni2+ and Cu2+ than in the case of Zn2+. The natural charge on the N4 is −0.4842 in the free ligand but −0.03501 and −0.3524 in the Ni2+ and Cu2+ complexes, respectively. It however increased to −0.6791 in the Zn2+ complex.

Similarly, natural charges on C1, C2, C3, N8, and C9 which decreased in the Ni-1MeIm and Cu-1MeIm complexes, however, increased in the Zn2+ analogue. These support the fact that although there is substantial delocalization of ligand-bonding electrons onto Ni2+ and Cu2+ d-orbitals [30], these electrons are, however, localized on the ligand ring of the Zn-1MeIm complex, making it more nucleophilic.

The 3d-orbital occupancy and energies of the complexes are assessed in Table 4. It is shown that the t2 orbitals (, , and ) of Ni-1MeIm have higher energies (−0.611 eV), with the e orbitals ( and ) having lower energy (−0.626 eV). This is typical of metal d-orbitals in a tetrahedral field [27].

There is, however, a reduced occupancy (≈0.9950) in all the orbitals. These electrons could be delocalized onto antibonding or non-Lewis type (Rydberg) molecular orbitals [31]. Similar occupancies were observed in the Cu2+ ion. With regard to the energy levels, the and orbitals were degenerated (−0.612 eV), but the orbital sunk low in energy (−0.660 eV) below the e energy levels. The orbital was also lower in energy than the orbital. These differences in the energy support the Jahn–Teller distortions observed in Table 1 as a result of unequal occupancy of the t2 orbitals in a tetrahedral field. Similarly, the Cu-1MeIm bond lengths also showed variations as observed in Table 1. In the case of Zn-1MeIm, the d-orbital occupancy was the same (≈2), showing that the d-electrons were localized within the Zn2+ d-orbitals [12]. Furthermore, the d-orbitals were degenerated with no d-d splitting. This explains the uniformity in the Zn-1MeIm bond lengths observed in Table 1.

To further probe the low occupancy of the Ni2+ and Cu2+ d-orbitals, second-order perturbation analysis of the Fock matrix was performed to investigate the interaction between metal d-orbitals (donor) and non-Lewis orbitals (acceptor) on the coordinating atom of the ligands. Tables 57 summarize the interaction between the d-electrons on the metal ions and the acceptor sites on the respective ligands. The highest energies, E(2), were recorded for all the complexes.

On the contrary, weaker interactions were observed in the case of Zn-1MeIm. Notable interactions include donation of LP(5) on Zn49 to RY ∗ (2)N1, RY ∗ (2)N8, RY ∗ (2)N15, and RY ∗ (2)N22 of energies 0.11, 0.10, 0.10, and 0.10 kcal/mol, respectively.

These interaction energies are, however, lower than those observed in the Ni-1MeIm and Cu-1MeIm. This explains the higher occupancies of the Zn2+ d-orbitals compared with those of the Ni2+ and Cu2+ analogues as shown in Table 4.

The energies of the frontier orbitals are important in describing the chemical properties of molecules [32]. The energy of the highest occupied molecular orbital (HOMO) gives an indication of the electron-donating ability of the complex [33].

The higher the energy of the HOMO, the easier it is to donate electrons into the lowest unoccupied molecular orbital (LUMO). Also, the larger the energy difference between these frontier orbitals, the more stable the complex. Figure 2 shows the energy gap between the frontier orbitals in the respective complexes. The gap is in the order Ni-1MeIm > Cu-1MeIm < Zn-1MeIm. This order is in agreement with the Irving–Williams series [27].

4. Conclusion

This article assessed the coordination behavior of Ni2+, Cu2+, and Zn2+ ions in tetrahedral 1-methylimidazole ligand fields. Measurement of the bond lengths of the optimized structures revealed Jahn–Teller distortions in the Ni-1MeIm and Cu-1MeIm complexes. These distortions were, however, absent in Zn-1MeIm. A natural bond orbital (NBO) analysis of the molecular orbitals shows substantial decrease in the formal charges on the metal ions upon coordination to the ligands. This was in the order Ni2+ ˂ Cu2+ ˂ Zn2+ with Zn2+ being the most electropositive one. Second-order perturbation of the Fock matrix shows higher interaction energies between Ni2+ 3d-orbitals (donor) and ligand orbitals (acceptor) than in the case of Zn2+ 3d-orbitals (donor). Finally, the frontier orbitals were assessed as a measure of the stability of the complexes. The energy of the gap was in the order Ni-1MeIm > Cu-1MeIm < Zn-1MeIm which is consistent with the Irving–Williams series.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the preparation and publication of this manuscript.

Acknowledgments

The author is grateful to the Cambridge Crystallographic Data Centre (CCDC) for the opportunity to use the Cambridge Structural Database (CSD) for the substructure molecular searches and visualizations.