Abstract

The changes in the interatomic distances and the corresponding spin-spin coupling as a result of the hydrolysis of the ruthenium complexes and the effects of different derivatives of the pyrazole ligands and the substituents methyl, carboxylic, and phenyl on the pyrazole rings were studied. A good agreement was obtained between the experimental and the theoretical proton NMR. Significant changes are observed in the isotropic and anisotropic shielding tensor of the atoms and related spin-spin coupling of their bonds due to hydrolysis of the complexes. This observation gives more insight into the known mechanism of activation of the ruthenium complexes by hydrolysis. There are no direct effects of interatomic distances on many of the computed spin-spin couplings with the exception of 1J(Ru-N) which shows significant changes especially within the pair of 1J(Ru-N) in the complexes with two nitrogen atoms of the bis-pyrazole moiety. The magnitude of interatomic spin-spin coupling of the Ru-X follows the order of Ru-Cl > Ru-N > Ru-C > Ru-O. The Ramsey term Fermi contact (FC) has the most significant contribution in most of the computed spin-spin interactions except in 1J(Ru-Cl) and 1J() which are predominantly defined by the contribution from the paramagnetic spin orbit (PSO).

1. Introduction

Molecular properties like the spin-spin coupling constant between a pair of atoms in a molecule can give very useful information for analyses of their structure and onset of intramolecular bonding interactions [1, 2]. The spin-spin coupling is an interaction between the magnetic moments of the coupling nuclei. Unlike direct dipolar coupling which is a through-space interaction, the spin-spin coupling is sensitive to bonding interactions and is mediated by the polarization of the spins of the intervening bonding electrons [3]. The classical NMR experiments [4] of the molecules dissolved in liquid crystals [5, 6] are a well-known experimental approach for the determination of the spin-spin interactions but interpretation of the experimental data is often challenging [4]. Through the computational method, the interpretation of the spin-spin coupling is much easier by partition of the spin-spin property into contributing Ramsey terms [7] which are made up of the Fermi contact (FC), spin dipole (SD), diamagnetic spin orbit (DSO), and paramagnetic spin orbit (PSO) [8].

The general expectation is that there should be diminishing in the values of the coupling interactions as the number of bonds between two nuclei increases [1, 2]. However, it is possible to obtain large coupling constants even when coupling atoms are many bonds apart as far as there is close proximity between the two nuclei. In this research six η6-cymene ruthenium complexes of pyrazole derivatives (Figure 1) which have methyl, carboxylic, and phenyl group on C3 and C5 were synthesised and their proton NMR spectra were related to the theoretical modelling. The effects of hydrolysis on the NMR and intramolecular spin-spin interactions in the mono- and bidentate derivatives of pyrazole complexes were considered. This makes the total number of theoretical molecules added up to eleven complexes. Hydrolysis is known to be an important process which defines the activation of metal complexes in their biological application as anticancer drugs [912].

Some of the research interests are to know the extent of the effects of pyrazole unit substituents and the types of the pyrazole ligand (mono-, bi-, and tridentate) on the ruthenium-ligand spin-spin coupling and also on the 3J(H⋯H) spin-spin coupling of the four hydrogen atoms (4H) of the benzene ring of the cymene ligand in the complexes. Also of interest is the change in the intramolecular spin-spin interactions in the complexes as a result of hydrolysis which is well defined mechanism for the activation of the metal complexes in anticancer application.

2. The Computational and Experimental Methods

2.1. Computational Details

All the calculations were accomplished using the Gaussian 09 (G09) version D01 [13]. The optimized geometries of the eleven derivatives of η6-cymene Ru(II) complexes were obtained using PBE exchange functional [14, 15] and no negative harmonic force constants were obtained in all the complexes as a clear indication that local minimums were obtained. During the optimization, all atoms in each of the complex were treated with 6-31+G(d,p) [1618] except Ru atom which was treated with basis set SBKJC VDZ [19] (obtained from EMSL Basis Set Library [20, 21]) and a total of 28 core electrons were removed from Ru (1s, 2s, 2p, 3s, 3p, and 3d) atom. The basis set SBKJC VDZ comes with effective core potential and has been applied in computing properties of many metal clusters [22, 23]. DFT PBE exchange functional has been found to be good for the optimization of metal complexes and give similar results with hybrid functional PBE0 and mPW1PW91 [24]. It is also good for studying weak-interacting systems [25] but just like any other DFT methods, it is limited in the study of noncovalent intermolecular interactions because of its long-range deficiency which makes hybrid DFT preferable [2628]. Even though all the properties of interest are purely intramolecular interaction, the NMR, spin-spin coupling constants J(A, B) [29, 30], and other stationary properties were computed using B3LYP hybrid DFT method which is Becke’s three-parameter exchange [31] and Lee-Yang-Parr’s correlation nonlocal functional. In the stationary phase computation of NMR and their spin-spin interactions, all electron basis set DGDZVP [32, 33] was used for Ru atom while other atoms were still treated with 6-31+G(d,p) as in optimization step. The Gauge-Independent Atomic Orbital (GIAO) method was used to compute the NMR properties. The direct and fitting methods were used for the calculation of the theoretical 1H-NMR chemical shift. In direct method, the isotropic shielding of each proton in the complexes was directly subtracted from that of the reference tetramethylsilane (TMS) while in fitting method the reported equation = 31.0 − 0.97 in the literature [34] was applied.

2.2. Experimental Synthesis of the Complexes

The first step is the synthesis of the precursor [(cym)RuCl2]2 by reacting RuCl3·3H2O with 8.0 equivalent mmol of p-phellandrene in ethanol. A similar approach as was reported in a literature was adopted [35] only that 8.0 mmol equivalent of p-phellandrene was used instead of approximately 10 mmol equivalent. All the complexes are then obtained from the precursor [η6-(cym)RuCl2]2 interaction with 2.2 equivalent of each ligand and reflux for 6–12 h. The synthesis of the complexes follows a similar method used for the synthesis of RAPTA complexes [36, 37].

3. Result and Discussion

3.1. The Geometries Properties

Six η6-cymene ruthenium complexes dmpzRu(II)Cym (1a), dcpzRu(II)Cym (2a), bdmpzmRu(II)Cym (3a), bdmpzaRu(II)Cym (4a), bphpzaRu(II)Cym (5a), and bdmpzpyRu(II)Cym (6) (Figure 1) were synthesised and their proton NMR property was determined. All the six synthesised complexes and their respective hydrated form dmpzRu(II)WCym (1b), dcpzRu(II)WCym (2b), bdmpzmRu(II)WCym (3b), bdmpzaRu(II)WCym (4b), and bphpzaRu(II)WCym (5b) with tridentate derivatives 6 were studied computationally. The peculiar features which are of interest in this research are the effect of the methyl, carboxylic, phenyl group, and hydrolysis on their NMR and spin-spin coupling.

The bond distances which define all the spin-spin coupling of interest in the complexes are shown in Table 1 and selected bond angles are shown in Table 2. The spin-spin couplings which are of interest are those which involved the one-bond distance coupling of ruthenium-ligand bonds (1J(Ru-X) where X stands for C, Cl, N, and O), three-bond distance (3J(Hc⋯Hc)) of the 4H atoms of the benzene ring of cymene ligand, the four-bond distance coupling of the two coordinating nitrogen atoms of the bidentate ligands 4J(N⋯N), three- and six-bond distance coupling of three coordinating nitrogen atoms of the tridentate ((3J(N⋯Np)) and 6J(N⋯N) where Np stands for the coordinating nitrogen atom of the pyridine unit), one-bond 1J(C-H) coupling of the C4-H in pyrazole unit, and the three-bond distance hydrogen atoms of the pyrazole unit in the ligand 5a (3J(H⋯H)). There is gradual decrease in the N⋯N interatomic distance of complexes with bidentate as a result of hydrolysis which should be direct effect of their decrease in N-Ru-N bond angle (Table 2). The bond distances of N-H decrease also upon hydrolysis while C-H bond distances are relatively the same in all the complexes. All the metal-ligand bonds (Ru-X) follow the order of Ru-Cl > Ru-C, Ru-O > Ru-N.

It is interesting to point out the relationship in the properties of the pair of Ru-N bonds in the bis-pyrazole derivatives of the complexes. The bidentate and tridentate derivatives show a very little difference between bond distances of the pair of Ru-N ranges from 0.01 Å in 3a and its hydrated form 3b to 0.03 Å in complex 4b (Table 1). The differences in the pair of Ru-N bonds’ distances are not the direct effect of their monodentate ligand Cl or H2O position since the Cl atom in 4a forms almost the same angle with the two coordinating N atoms (having N-Ru-Cl bond angles of 84.52° and 84.13° as shown in Table 2) and yet has the highest variation in the pair of its Ru-N bonds compared to other complexes that have higher bond differences with their monodentate ligand Cl or H2O.

3.2. The 1H-NMR, Isotropic, and Anisotropic Shielding Tensor

The experimental 1H-NMR spectra of the complexes 3a, 5a, and 6 derivatives are shown in Figure 2 and the comparison of all the synthesised experimental shifts with the theoretical direct and fitting methods is shown in Table 3. The unique feature in complex 3a which is absent in the 5a and 6e is the proton shift of the pyrazole linking CH2 which shows a long singlet peak around −3.8948 (Figure 2). The pyridine unit in 6 makes the span of its hydrogen shift higher than 3a and the two phenyl units in the complex 5a result in a much more wider shift of its proton atoms. There is a good agreement between the experimental and the theoretical 1H-NMR especially for the CH3 substituent (Table 3). The methyl protons which are pointing towards the Ru atom are very low in shift while those pointing towards the chloride ion are higher in chemical shifts. It is only in complex 6 that proton shift of C4-H of the pyrazole unit spans a wider range theoretically compared to their experimental range.

There is a significant increase in the magnitude of both isotropic and anisotropic shielding tensor of ruthenium atom due to hydrolysis. This can serve as one of the indicators for the activation process of the metal complexes in their anticancer application. Likewise, the shielding of the coordinating nitrogen atoms and their anisotropic tensor increases significantly upon hydrolysis. In all of the selected atoms in Table 4, the magnitudes of the isotropic shielding tensor of the atoms are lower than their corresponding anisotropic shielding tensors except for the Cl atoms. The isotropic shielding tensor of the coordinating nitrogen atoms (N) is lower than that of the other nitrogen atoms () which are not involved in coordination while their anisotropic tensors are reverse of the isotropic tensor. It is only in the complex 2a that the coordinating nitrogen is found to be deshielded (negative isotropic value) due to the presence of two carboxylic groups on the pyrazole unit. The isotropic shielding and anisotropic tensors of Hc atoms of cymene are little higher than Hz atoms of the pyrazole. The isotropic shielding of Hc atoms is within a very close range in all the complexes but their anisotropic tensors are found higher in the complexes with monodentate pyrazole units. In the complexes with monodentate pyrazole unit, the anisotropic properties of the Hn atoms significantly reduce upon hydrolysis while there is a slight increase in their isotropic shielding tensors.

3.3. The Spin-Spin Coupling Interaction

The graphical representation of selected spin-spin interactions in the complexes is shown in Figure 3 and the details of all their spin-spin are shown in Table 5. In many of the experimental characterisations of derivatives of η6-cymene Ru(II) complexes, the values of the spin-spin 3J(Hc⋯Hc) of benzene ring are often presented as part of their characterisation and have been reported in the literature to have experimental ranges of 6–8 [38, 39], 6 [40, 41], 5.7–6.2 [42], 6.4 [43], and 6.02–6.56 [44] which agree with our computed ranges of 6.124 in 1a to 7.151 in 6. The changes in the spin-spin coupling of the complexes as corresponding effect of the type of the monodentate, bidentate, and tridentate ligands and hydrolysis of the complexes are examined. All the spin-spin couplings of 1J(Ru-C), 1J(Ru-Cl), 1J(Ru-N), 1(N-N), 3J(N⋯Np), and are negative in values while 1J(Ru-O), 1J(C-H), 3J(Hc⋯Hc), and 3J(H⋯H) are positive. The conventional signs of the coupling do not directly follow the proposed convection of positive for one-bond coupling constants, negative for two-bond coupling constants, and positive for three-bond coupling constants based on the generalized Dirac vector mode [4]. It is expected that as the number of bonds between the two nuclei increases their coupling should diminish [1, 2]. The magnitude of interatomic spin-spin coupling of the Ru-X follows the order of Ru-Cl > Ru-N > Ru-C > Ru-O which does not correspond to the expected reverse of their corresponding interatomic bond distances. The possible explanation for the higher spin-spin coupling of Ru-Cl above other Ru bonds in the compounds is the predominant HOMO feature of the Cl atom in the complexes as shown in Figure 4. It would also be expected that values of the coupling interaction of the 1J(N-N) are higher than all other spin-spin couplings because of its lower bond distance except for the 1J(C-H) and 1J(N-H) (Table 1), but its values are lower than many of the other interatomic couplings with longer bond distances. In addition, since there is a gradual decrease in N⋯N interatomic distances upon hydrolysis of complexes with bidentate nitrogen chelating ligands (Table 1), one would expect the spin-spin interaction of 3J(N⋯N) in hydrated complexes to be higher but the reverse is the case (Table 5). Also, the bond distances of N⋯Np (Table 1) are shorter than many of the N⋯N bond distances but its spin-spin coupling values are lower than many of the spin-spin values of the N⋯N (Table 5). The coupling and FC of some intramolecular atoms of molecules have been found to be independent of the bond types and distances [4] unlike the coupling of intermolecular hydrogen bonding. The reason has been traced to the consequence of the absence of an explicit distance-dependent term in the FC operator [4]. It is possible to obtain large coupling constants even when lying many bonds apart as far as there is close proximity between two nuclei.

The little differences in the Ru-N bond distances affect their pair of spin-spin couplings 1J(Ru-N) as the complexes 4a and 5a which have a bond differences of 0.03 and 0.02 Å in their pair of Ru-N also have the highest differences between their pair of spin-spin couplings, 1.188 and 0.944, respectively, compared to other complexes (Table 5). Out of all the 1J(Ru-O) of the hydrated complexes, 5b has the highest magnitude while 2b has the lowest (Figure 3 and Table 5) which is inversely proportional to their bond distances (Table 1). All the 1J(Ru-C) span the range of −7.033 in 6 to −14.610 in 3b. The magnitude of spin-spin interaction 1J(Ru-Cl) is lower in monodentate compared to bidentate complexes and even much lower after hydrolysis (Figure 3). There is significant increase in the magnitude of 1J(C-H) and 1J(N-H) upon hydrolysis. The coupling of 1J(N-N) decreases slightly upon hydrolysis (except for its increased hydrolysed complex 2b compared to its nonhydrolysed form 2a) and its values range from −1.293 in the complex 2a to −2.043 in complex 5a.

3.4. Ramsey Terms and Their Correlation with Spin-Spin Interatomic Interaction

The Ramsey terms SD, FC, PSO, and DSO which defined each of the spin-spin coupling interactions for the selected interatomic distances of interest in the complexes are presented in Supplementary Table S1 (see Supplementary Material available online at http://dx.doi.org/10.1155/2015/832045). The most significant Ramsey term that contributes mostly to the spin-spin interactions is FC except for the Ru-Cl and which are predominantly defined by the contribution from the PSO. Also, PSO contribution to the Ru-O spin-spin interaction is higher than FC in complexes with monodentate pyrazole derivatives. The properties of the Ramsey terms for the Ru-Cl are very unique; after the most significant contribution from the PSO, there is also a significant contribution from SD which in complexes like 3a and 4a is higher than the contribution from the FC. Also, PSO is the most significant Ramsey term in all the spin-spin followed by the contribution of FC except for the complexes 4a, 5a, and 6 where FC contribution is higher than the PSO. The entire Ru-C and Ru-N spin-spin coupling are defined mostly by the contribution from their FC followed by small contribution of the PSO but the PSO contribution improves slightly upon hydrolysis. Also FC is the most significant Ramsey term which defines all H⋯H, C-H, and N-H spin-spin coupling but next contributing Ramsey terms is their DSO.

The correlation () and the linear equation of the Ramsey terms, the interatomic distances, and isotropic and anisotropic tensors of the two atoms involved in the spin-spin coupling and their respective sum are shown in Figure 5 while the correlations within all the terms are shown in Table 6. Just as the FC is found to be the most significant Ramsey terms that define many of the interatomic spin-spin interactions in the complexes, the trend of changes in its values is also found to be highly correlated () with the trend of changes in all the spin-spin interactions of the complexes (Figure 5). Even though in many of the interactions of interest the PSO contribution is significant yet it has a very poor correlation () with the trend of changes in the spin-spin interaction that is poorer than the SD term () which has a very negligible contribution to most of the spin-spin interactions. The correlation of DSO () with the trend of changes in the spin-spin coupling is relatively better followed by the correlation of their respective interatomic distances (). As shown in Table 6, there is significant and inverse correlation between the interatomic distances and their FC and spin-spin interaction. There is significant correlation between the SD and PSO and also between the FC and DSO. The Iso1 and Aniso1 which are predominantly from the ruthenium atoms are significantly correlated with the DSO while the Iso2 and Aniso2 are significantly correlated with the SD and PSO. The averages of the isotropic and anisotropic of the two atoms involved in the spin-spin coupling are significantly correlated mostly with the DSO.

4. Conclusion

The experimental and the theoretical proton 1H-NMR of η6-cymene Ru(II) complexes with monodentate, bidentate, and tridentate of pyrazole derivatives were studied and a very good agreement was obtained between the experimental and theoretical methods. The features and the relations between the geometrical interatomic distances, angles, and the spin-spin coupling which involve all the Ru-ligand bonds and one bond or many bonds’ interaction of interest as results of the nature of the substituents and the ligands were presented. Through the theoretical modelling, the very low proton shifts in the experimental results are found to be the methyl protons which are pointing towards the Ru atom while those pointing towards the chloride ion have higher chemical shifts. A gradual decrease in the N⋯N interatomic distance of complexes was observed in the bidentate derivatives as a result of hydrolysis which is a direct effect of the decrease in the N-Ru-N bond angle of the hydrolysed complexes compared to the nonhydrolysed ones. There is a significant increase in the magnitude of both isotropic and anisotropic shielding tensor of ruthenium atom due to hydrolysis. The anisotropic properties of the Hn atoms of complexes with monodentate pyrazole unit significantly reduce upon hydrolysis while their isotropic shielding tensors slightly increased. Likewise, the shielding of the coordinating nitrogen atoms and their anisotropic tensor increases significantly upon hydrolysis. The changes in the anisotropic and isotropic tensors of the complexes upon hydrolysis can serve as a good parameter to define the mechanism of their activation by hydrolysis which is important for their application as anticancer. The computed spin-spin coupling of the benzene ring of the cymene unit of the complexes is within the experimental ranges that have been reported for similar complexes. All the spin-spin couplings of 1J(Ru-C), 1J(Ru-Cl), 1J(Ru-N), 1(N-N), 3J(N⋯Np), and are negative in values while 1J(Ru-O), 1J(C-H), 3J(Hc⋯Hc), and 3J(H⋯H) are positive. In many of the computed spin-spin couplings, there is no direct relation between their interatomic distances and their spin-spin couplings. The magnitude of interatomic spin-spin coupling of the Ru-X follows the order of Ru-Cl > Ru-N > Ru-C > Ru-O which does not correspond to the expected reverse of their corresponding interatomic bond distances. The little difference between bonds distances of the pair of Ru-N in the complexes with two pyrazole coordinating nitrogen results also in corresponding differences in their spin-spin coupling. The magnitude of spin-spin interaction 1J(Ru-Cl) is lower in monodentate compared to bidentate complexes and even much lower upon their hydrolysis. Also, there is significant increase in the magnitude of 1J(C4-H) and 1J(N-H) upon hydrolysis. While the coupling of 1J(N-N) decreases slightly upon hydrolysis (except for its increased hydrolysed complex 2b compared to its nonhydrolysed form 2a). The most significant Ramsey term that contributes significantly in most of the spin-spin interactions is FC except for the Ru-Cl and which are predominantly defined by the contribution from the PSO. Also, it is only in 1J(Ru-Cl) that the contribution from SD was found higher than the FC contribution. Even though in many of the interactions of interest the PSO contribution is significant than DSO and SD values yet it has poorer correlation with the trend of changes in the spin-spin interaction than DSO and SD.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors gracefully acknowledge the financial support of Govan Mbeki Research and Development Centre, University of Fort Hare, South Africa. The CHPC in Republic of South Africa is gracefully acknowledged for providing the computing facilities and some of the software used for the computation.

Supplementary Materials

Table S1: The Ramsey terms and the spin-spin J-coupling interaction of selected bonds in complexes of pyrazole derivatives.

  1. Supplementary Material