Research Article  Open Access
Adebayo A. Adeniyi, Peter A. Ajibade, "Insights into the Intramolecular Properties of η6AreneRuBased Anticancer Complexes Using Quantum Calculations", Journal of Chemistry, vol. 2013, Article ID 892052, 14 pages, 2013. https://doi.org/10.1155/2013/892052
Insights into the Intramolecular Properties of η6AreneRuBased Anticancer Complexes Using Quantum Calculations
Abstract
The factors that determine the stability and the effects of noncovalent interaction on the 6arene ruthenium anticancer complexes are determined using DFT method. The intramolecular and intraatomic properties were computed for two models of these halfsandwich ruthenium anticancer complexes and their respective hydrated forms. The results showed that the stability of these complexes depends largely on the network of hydrogen bonds (HB), strong nature of charge transfer, polarizability, and electrostatic energies that exist within the complexes. The hydrogen bonds strength was found to be related to the reported anticancer activities and the activation of the complexes by hydration. The metal–ligand bonds were found to be closed shell systems that are characterised by high positive Laplacian values of electron density. Two of the complexes are found to be predominantly characterised by LMCT while the other two are predominately characterised by MLCT.
1. Introduction
There have been several research efforts to synthesize Rubased anticancer complexes as alternative to cisplatin in cancer therapy [1–3]. Among the most studied compounds are the halfsandwich complexes of ruthenium due to their unique properties [4–6]. Among the most studied complexes are the halfsandwich complexes of ruthenium. Several of these halfsandwich ruthenium complexes have found numerous applications as catalysts for organic transformations, in the supramolecular field and in medicinal chemistry [7]. The applications of these complexes as anticancer agent have been reported [3, 8–12].
Some of the properties of interest are the existing noncovalent interactions and the effect of hydration on the interatomic interactions in the complexes. The noncovalent interactions such as hydrogen bonding, anionπ, cationπ, and ππ interactions and other weak forces are important in chemical reactions, molecular recognition, and regulation of biochemical processes [13, 14]. Deep understanding of these interactions has been pointed out to be of great importance in rationalizing their effects [14]. Using Bader’s quantum theory of atoms in molecules (QTAIM) [15], the atomic properties such as electronic population, energies, and (de)localization are evaluated over the atomic basins. Computer simulation is known to be helpful in giving detailed atomic structural properties and in interpreting experimental data at atomic level of interaction to show the mechanisms of biomolecular function [16]. Also, the quantum calculation plays significant roles in determination of force fields [17, 18] necessary for the in silico drug designs which is known to be pivotal in discovering new drugs and designing more efficient ones [19, 20].
In this research work, we have selected two of the models compounds (Figure 1) which are [Ru(η6pbenzene)Cl_{2}(pta)] (RAPTAH named as complexes 1 and 2) and [Ru(η6pcymene)Cl_{2}(pta)] (RAPTAC named as complexes 3 and 4) reported by Chatterjee et al. as anticancer agents [21]. The hydrated form of these complexes which are [Ru(η6pbenzene)Cl(H_{2}O)(pta)] named complex 2 and [Ru(η6pcymene)Cl(H_{2}O)(pta)] named complex 4 (Figure 1) is considered since activation of Ru complexes is known to occur through hydration [10, 22–24]. Complexes with PTA ligand have in recent years received attention because of their water solubility and applications as catalyst [25]. The only difference between complexes 1 and 3 is the use of cymene as arene unit in 3 while benzene is used in complex 1. This little change in ligand has been reported to enhance the anticancer activities of complex 3 compare to 1 [21]. The choice of ligand is important because a too strongly bound ligand could render the drug inactive, while a labile ligand could be easily hydrolysed or replaced [26]. Many of these rutheniumarene complexes are known to have complicated and unstable ligand exchange [2]. In order to improve their anticancer activities and obtain a better lead compounds, their stability must be improved [2]. In this study, we have used the quantum theory of atoms in molecules (QTAIM) to understand the effects of noncovalent interactions on the stability and hydration of these complexes.
(a)
(b)
(c)
(d)
2. Computational Method
In this work, the geometries of the complexes were first optimized with PBE0 [27] functional and mixed basis sets SBKJC VDZ with effective core potential (ESP) [28] for Ru, P, and Cl while basis set 631G* was applied on other atoms in each of the complexes (this will subsequently be referred to as ECP (Ru,P,Cl)631G*). In the second optimization, the SBKJC VDZ is limited to only the ruthenium atom while the scaledup basis set 631+G(d,p) was applied on other atoms in the complexes which shall subsequently be referred to as ECP(Ru)631+G(d,p). The external basis sets were obtained from EMSL basis set library [29, 30] and were incorporated into the input files in a format that each FIREFLY and Gaussian 09 (G09) can read. SBKJC VDZ ECP basis set with PBE0 functional has been shown to be effective in treating complexes with large number of electrons and has been applied in computing properties of many metal clusters [31, 32]. Other properties of the complexes are computed at B3LYP hybrid functional level of theories [33] using basis set DGDVZP applied on Ru atom while others are treated with 631+G(d,p) which will be referred to as DGDVZP(Ru)631+G(d,p) subsequently. Also lower basis set 321G [34] was applied on all atoms of the complexes in order to compare its values with DGDVZP(Ru)631+G(d,p) systems. The Bader quantum theory of atoms in molecules (QTAIM) analysis was done mainly using the wavefunction obtained from both DGDVZP(Ru)631+G(d,p) and 321G basis sets treated systems. A topological analysis was performed in order to calculate the charge density and its second Laplacian derivative of charge density for the bond critical points (BCP).
The major computation tools used are FIREFLY 7.1.G [35] which is partially based on the GAMESS (US) [36] source code, Gaussian 09 [37], and AIMAll 12.06.03 [38]. NBO 5.0G program [39] as implemented in FIREFLY 7.1.G was used for the natural bond orbital (NBO) analysis [40] and natural energy decomposition analysis (NEDA) [41].
3. Result and Discussion
Stable stationary geometries were obtained for the four complexes which are characterized with zero imaginary number from the thermodynamic calculation. When the geometries of the complexes were further optimized from PBE0/ECP(Ru,P,Cl)631G* to PBE0/ECP(Ru)631+G(d,p), the highest contraction of bond lengths up to 0.05 Å was observed for Ru–Cl bonds (Table 1). There is also a slight change in bond orders up to 0.1 when comparing the two steps of optimization. Applying uniform basis set 321G on all atoms of the complexes, there is no significant change in the total energies from the first optimization (−6307.84, −5925.61, −6464.66, and −6082.04 A.U., resp.) to the second optimization (−6307.84, −5925.61, −6464.26, and −6082.04 A.U.). The bond orders of the Ru–Cl bonds increase in the hydrated complexes 2 and 4 which suggest that the possibility of displacing the second Cl atom with aqua ligand will be difficult [42]. But on the other hand, the Ru–P bond order decreases upon hydration of the complexes. The properties computed using combination of higher basis sets DGDVZP(Ru)631+G(d,p) and those computed using lower basis sets 321G follow the same trend which is an indication that lower basis sets can give insight into the trend though it either underestimates or overestimates the computed properties.

3.1. The Natural Bond Orbitals (NBO) Analysis
The total electron densities are characterised mostly as Lewis electrons which are higher than 97% in all the complexes indicating that they are stable. Also, the negligible nature of extravalence orbitals (i.e., Rydberg) which is a slight departure from a localized Lewis structure model (0.14% in all the complexes) suggests the stability of the complexes. However, the relatively high electron delocalization which is the nonLewis orbitals (little above 2.00%) indicates that they play significant roles in the stability of the complexes. The 4dorbitals of the ruthenium atom are lower in energy than 5sorbital and are used preferentially in bonding with other atoms in the complexes.
The polarization features of some visible Ru–L bonding orbitals are shown in Table 2, while Table 3 shows the nature of the electron delocalization orbitals with their secondorder perturbation energy (). The number of the electron transfer into each of the acceptor orbitals in Table 3 are shown in Table 4 with the values of their energy level. The bonding polarizations of the Ru–L are directed towards the ligand atoms (Table 2). The reverse is the case in all the antibonding orbitals interactions as a result of backbonding of electrons into the antibonding lone pair of ruthenium atoms. The features of the charge transfer are clearly seen from the delocalized orbitals that have notable interaction with the Ru atom or bonds and have significant stabilization energy equal or greater than 10 kcal/mol as presented in Table 3. There are many bond to bond, atom to bond, and atom to atom delocalization orbitals of electrons in the complexes. The most significant feature is the presence of metal to ligand charge transfer (MLCT) from the lone pair of Ru to the antibonding lone pair of arene C atoms in all the complexes. There is also an observed MLCT from the lone pair of ruthenium to the bonding lone pair of the same arene C atoms in complexes 2 and 3. In all the complexes there is also ligand to metal charge transfer from the arene C atoms into the antibonding lone pair of ruthenium atoms or antibonding orbitals of the rutheniumligand (σRuL) bonds. All the backbonding of electronic charges is responsible for the significant effect of the nonLewis backbonding contribution to the stability of the complexes. From the NEDA analysis, if the total charge transfer from ligand to the lone pair of ruthenium metals (LMCT) is subtracted from that of metal to ligand (MLCT), the overall features show that complexes 1 and 4 are predominantly characterised with LMCT with respective total charges of −0.56209 and −0.47667 transfer to the lone pair of ruthenium atoms from ligand atoms while complexes 2 and 3 are predominately characterised with MLCT of total charges of 0.37489 and 0.45986 transfer from the Ru lone pair to the ligand atoms on bonds. These features are typical nature of π ligand metal complexes with intense intramolecular CT between metal and ligand (MLCT or LMCT) transitions and associated with π backbonding [43, 44].
 
Polarization coefficient is the values with starred superscript, and the square of it is percentage of the NBO ( squared) on each hybrid orbital (in parentheses). 
 
The orbital analyses in the table are defined in terms of the second perturbation energy or stability energy ( in kcal/mol). 

The features of the electronic orbitals in Figure 2 give an insight into observed MLCT and LMCT in the complexes. The ruthenium atom is characterised as part of both HOMO and LUMO as a result of bonding electron contribution of metal and the backbonding contribution of the ligands. The PTA ligand dominates the HOMO orbitals of the hydrated complexes 2 and 4 compared to the Ru and Cl atoms which characterise the HOMO of the unhydrated complexes 1 and 3. The arene ligand is predominantly the LUMO which confirms the existence of MLCT between the Ru atom and the π orbital of the arene moiety of the ligand. The HOMO nature of the Ru and the Cl atom in complexes 1 and 3, significantly reduced upon being hydrolysed to complexes 2 and 4, respectively.
(a)
(b)
3.2. Natural Energy Decomposition Analysis
Results obtained from the NEDA analysis are shown in Table 5. The complexes 1, 2, 3, and 4 are fragmented into four, three, two, and six units, respectively. The respective fragmentations of complexes 1, 2, 3, and 4 are shown in Table 6 except the last two fragments of complex 4 which are the two hydrogen atoms of the aqua ligand that have negligible zero contribution. The number of the fragmentations in the complexes significantly correlates with their total interaction energy (Table 5). The strength of the interaction energy does not correlate with the number of hydrogen bond in the complexes. There are two, one, four, and five hydrogen bonds in complexes 1, 2, 3, and 4, respectively, as shown in Figure 3. The result shows that either polarization (POL) (complexes 1, 2, and 3) or charge transfer (CT) (complex 4) has the greatest contribution to the stability of the complexes. The significant values of the POL and ES result in the electrical energy having the highest contribution to the stability of the complexes (Table 5). The hydrated complex 4 has the highest interaction energy as a result of its significant high CT, ES, and POL. On the contrary, the hydrated complex 2 has lower interaction energy compared to its unhydrated complex 1. The high stability of interaction energy of RAPTAC (i.e., hydrated complex 4) should significantly contribute to its experimentally reported higher anticancer activity than the other RAPTA complexes [21]. In all the complexes, the high negative values of CT, ES, and POL balance up the strong core repulsion energy that is made up of the sum of exchange interaction (XC) and deformation energy (DEF).


(a)
(b)
(c)
(d)
The first fragment in each of the complexes comprises the Ru and the PTA atoms while the second fragment is made up of their respective arene ligands. The arene fragment stability is greatly enhanced in the presence of other fragments in the complexes resulting in high induced energy (Table 6). The significant high induced dipole moment of the remaining chloride atom of complex 4 will definitely enhance the second hydration by exposing better to electrophilic attack. The high induced dipole moment and energy of the fragments in each complex are an indication that the fragments are synergistic. The improvement in the stability energy of one unit on the other is due to the synergistic effect of the ligands that are coordinated to the metal centre which is a unique property of metal complexes that is suggested to be responsible for the observed cytotoxic effect of many ligands [45]. The significant values of the induced dipoles in the complex especially from the acceptor unit imply that there is no strict orthogonality required from wavefunctions of the fragments which is responsible for the significant polarization, thereby exposing them as better acceptor and donor to nucleophilic and electrophilic attacks.
3.3. The Bond Critical Points Analysis
The bonds properties from both the B3LYP/DGDVZP(Ru)631+G(d,p) treated systems (Table 7) and B3LYP/321G (see Supplementary Table S1 available online at http://dx.doi.org/10.1155/2013/892052) methods of computation were obtained from the quantum theory of atoms in molecules (QTAIM) analysis of their respective electron density wavefunctions. The features of the interatomic interaction were explained in terms of the critical points density to obtain information about the existence of bonds and the Laplacian of the electron density at that point which gives the kind of interactions and their strength based on the sign and magnitude, respectively. The and of all the Ru–L bonds are shown in Table 7 and S1. All the Ru–L bonds of the complexes are characterized with positive but higher (Table 7 and S1 and Figure 3) than hydrogen bonds confirming that they are closed shell interactions like dative, hydrogen, ionic, and van der Waals bonds [46, 47]. The features of the bonds are shown using the contour plot of the along the plane of the P, Ru, and Cl or O atoms where the hydrated complexes are considered (Figure 3). Many of the atoms of the complexes are found to be on the chosen plane and all the ligand–ligand bonds are characterised with negative (solid lines) while the metalligand bonds are characterised with positive (dash lines). All the covalent bonds in the complexes are characterised with high negative values of , lower , and higher values of the Laplacian of the electronnuclear attractive contribution to virial field (Ven). All the strong noncovalent bonds like Ru–L in the complexes are also characterised with the same properties of the covalent bond except that their values are lower and the has positive values.
 
Ven: Laplacian of the electronnuclear attractive contribution to virial field, : Laplacian of the electron density. 
Complexes 1, 2, 3, and 4 have two, one, four, and five hydrogen bonds, respectively (Figure 3). In the literature there have been reports of unusual alkyl and halogen Hbonds [48–54] which are also observed in these complexes. There are unusual Hbonds between the alkyl C atoms of the cymene unit of complexes 3 and 4 and also H bonds between the chloride atoms and H atom of aqua, PTA, and arene ligands which contribute significantly to the H bond networks of complexes 1, 3, and 4. Hydrogen bonding is known to be a stabilizing interaction that influences the structure, function, and dynamics of many chemical and biological systems in the gaseous, liquid and solid states [49, 50]. Hydrogen bonds are significantly important in defining the crystal packing of many organic and organometallic molecules and the 3D structure of biological macromolecules as well as modulating the reactivity of different groups within a molecule [49]. Peculiar to complexes 3 and 4, one chloride atom forms two hydrogen bonds with the nearby hydrogen atoms from PTA, aqua, and arene ligands. The introduction of cymene as source of π bond in the arene moiety of complexes 3 and 4 other than the benzene of complexes 1 and 2 brings about significant differences in their HB network. This network of HB further gives insight into the nature of electron communication that results in the synergistic effects of the coordinated ligands. The higher HB network in complexes 3 and 4 indicates that many of the atoms in these complexes will be sensitive to macromolecular interactions and should be responsible for their experimental reports of high cytotoxic activities than complexes 1 and 2 [21]. The hydrated forms (complexes 2 and 4) are characterised with extra and stronger HB suggesting that the activation mechanism of these complexes by hydration [10, 22–24] can possibly be the result of increase in the sensitivity of the atoms as a result of increased networks of HB and electronic interactions. There is a unique HB that exists between Cl and H atoms of the water molecules in the hydrated forms that is stronger than any other existing HB based on a higher value of and (Figure 3). The values of the ratio of the magnitude of potential energy (P.E) and the kinetic energy (K.E) of the electrons in Table 7 and S1 are relatively low indicating a balance between stronger bond interactions and high density in the regions [15]. A high value of corresponds to higher density and stronger bonds except where the is small (Table 7 and S1 and correlation in Table S2). The unique HB in the hydrated complexes that are also recognised in the bond order analysis is characterised with higher values of which is closed to that of the metal–ligand bonds than other HB. The very low average kinetic energy which is responsible for higher has been pointed out through the uncertainty principle to be the nature of loosely bound density [15].
The strength of the two Ru–Cl bonds in terms of the and is not equal in complex 3 compared to complex 1 (Table 7 and S1) which is also observed in their bond order (Table 1). The imbalance features of the two Ru–Cl in complex 3 can enhance the aqua substitution of the weaker Ru–Cl to form Ru–OH_{2} which may also play significant role in the reported higher anticancer activity of RAPTAC since easy hydration has direct correlation with the higher anticancer activity [24]. The of Ru–P bond in all the complexes is found to be the lowest in all the existing metal–ligand bonds which would have suggested it as the most possible leaving unit except for its higher , lower , and relatively high Ven values (Table 7 and S1). The of Ru–P bond slightly reduces in value in the hydrated complexes 2 and 4 compared to complexes 1 and 3, respectively, while the of the RuCl bond of the retained Cl atom slightly increased just as it was observed for the bond order analysis.
The relationships between the computed factors of the bonds are constructed over all the existing bonds in the complexes (Supplementary Table S2). A very high negative value of is an indication of strong covalent bond while a high positive value corresponds to a strong noncovalent bond. The general relation is that high negative value of is directly proportional to high , high Ven, lower bond stretch (BPLGBL_I), lower kinetic energy , high negative values of potential energy , lower ellipticity , and relatively high value of . The higher electronic kinetic energy is an indication of lower values of , , and potential energy but a higher value of (Supplementary Table S2). These observed correlations further confirmed the reported nature of very flat electron density region that is characterised with a very low average values for , which is usually found to have a relatively high [55].
The QTAIM analysis using the wavefunction file obtained from the mixed basis set where the ECP is applied on Ru, P, and Cl atoms (PBE0/ECP(Ru,P,Cl)631G* treated systems) gives spurious topological features of the complexes. We observed common nonnuclear attractor (NNA) critical points (CP), that is, (3, −3), between the Ru–P and Ru–Cl bonds which were completely absent when all electron basis sets DGDVZP(Ru)631+G(d,p) and 321G were used. The existence of NNA in only ECP models suggests it is just a computational artefact of ECP basis set. However, the existence of NNA has been reported to be not just an artefact of computational methods but a genuine feature of the electron density distribution of some complexes [56]. An instance of the first unambiguous experimental evidence of such a feature in a stable molecule was reported by Platts et al. [15]. Our observation shows that the pseudopotentials that are used to represent the core electrons cannot conveniently account for all the electrons distribution around these complexes and it is responsible for the introduced NNA. This agrees with report which shows that NNA originate from the shape of valence molecular orbitals that might occur in bonds of low polarity in which core contributions are negligible and the radial form of the valence orbitals dominates the total density [15, 57].
3.4. The IntraAtomic Properties
The properties of some selected atoms in each of the four complexes are presented in Table 8 and supplementary Table S3. The integrated Lagrangian values of all the atomic basins are approximately equal to zero, which is an indication of satisfactory numerical integration [15]. The distribution of the atomic charges shows that the arene C atoms have gain charges confirming the existence of MLCT. For all the atoms in the complexes other than H atoms, the values of the percentage localization are higher than the percentage delocalization . The reverse nature of the and of the H atoms indicates that they can be easily perturbed by an external electric field [15]. The atomic volume (Vol) of Ru and P atoms increases in the respective hydrated complexes 2 and 4 while that of the Cl atom decreases. The of the Ru and arene C atoms increases in the respective hydrated complexes 2 and 4 while that of the P and Cl atoms decreases which consequentially leads to a reverse trend in their . Both bonding and total dipole moments of the Ru and Cl atoms are higher than the intraatomic dipole moment in all the complexes, but the reverse is the case of the P atom. The magnitude of the three dipole moments of the Ru atom decreases in the respective hydrated complexes 2 and 4. The bonding dipole of the P atom increased appreciably in the hydrated complexes 2 and 4 but still lower than its intraatomic dipole moment. The contribution of Ru atom to electronic energy of molecule slightly decreases in the hydrated complexes and also decreases sequentially from complexes 1 to 4 while the values of the increases. The Ru atom of complex 4 has a significantly higher and different value of its out of plane magnetizability contribution compared to its values in other complexes as a result of its higher intraatomic magnetizability than the bonding magnetizability. This further distinguishes hydrated RAPTAC complex which is experimentally reported as having better anticancer activities than others. The N atoms from the PTA which have been reported to usually take part in HB interaction of the complexes with macromolecule [58] are observed to have virtually the same charges in all the complexes and a very little reduction in the in the hydrated forms. The H atoms from the water molecule of the hydrated forms which are involved in HB are characterised with higher , lower , and lower compared to other H atoms from the PTA and arene units that are also involved in HB.
 
is net charge of atom , is Lagrangian of atom , is average number of electrons in atom , is electronic kinetic energy of atom (Hamiltonian form), % is percentage of average number of electrons localized in atom , % is the percentage of electron delocalization index of atom , is approximation to virialbased total energy of atom , is magnitude of intraatomic dipole moment of atom , is contribution of atom to electronic energy of molecule, Vol() is the volume bounded by an isosurface of the electron density distribution (0.001) and by interatomic surfaces of atom , and is out of plane magnetizability contribution of atom . 
The features of the relation within the computed intramolecular properties can be seen from the constructed correlation in Table S4. Higher is found to be as a result of higher bonding dipole moment contribution of the atoms , high number of electrons and relatively low that are poorly correlated. The high value of is an indication of low , high , and high of atoms. The atoms that are characterised with bigger volume will be associated with higher negative values of the . The atomic energy contribution to the virial energy of the system depends significantly on the number of the electrons and the number of localized electrons. The out of plane magnetizability is proportional to the number of the atomic localized electrons and inversely proportional to the atomic delocalized number of electrons.
The total sum of the intraatomic properties computed over all existing atoms shows the changes in the computed atomic properties from one complex to another (Supplementary Table S5). The decrease in the values of the atomic energy contributions like and of the respective hydrated complexes shows the significant contribution of the Cl atom to the total energy of the systems. The average values of the percentage electron delocalization over all the atoms increase in the respective hydrated complexes while the average electron localization decreases which correspond with the theory of activation of this complexes by hydration. Also, the average electron delocalization and the three types of dipole moments in complexes 3 and 4 are higher than complexes 1 and 2 which further distinguishes RAPTAC complex of reported better anticancer activities.
4. Conclusion
The quantum properties of complexes 1, 2, 3, and 4 were computed based on the NBO, NPA, NEDA, and QTAIM analyses using two mixed basis sets and single minimum basis set. Stable stationary geometries were obtained for the four complexes which are characterized with zero imaginary number. The bond orders, , and of the Ru–Cl bonds increase in the hydrated complexes 2 and 4 which suggests that the possibility of displacing the last Cl atom with aqua ligand will be difficult [42]. The results from the NBO and NEDA analyses with the features HOMO and LUMO show that complexes 1 and 4 are predominantly characterised with LMCT while complexes 2 and 3 are predominately characterised with MLCT. The result shows that either polarization (POL) (complexes 1, 2, and 3) or charge transfer (CT) (complex 4) has the greatest contribution to the stability of the complexes. The hydrated complex 4 has the highest interaction energy as a result of its significant high CT, ES, and POL which gives a possible distinguishing features that can contribute to the reported higher anticancer activities of RAPTAC [21]. In all the complexes, the high negative values of CT, ES and POL balance up the strong core repulsion energy that is made up of the sum of exchange interaction (XC) and deformation energy (DEF). The strength of the charge transfer within the complexes depends mostly on the fragmentation of each complex. The NEDA analysis of complexes shows that the coordinated ligands have high synergistic effects which may be responsible for their impacted anticancer activities. The differences in HB networks suggest a difference in sensitivity of the atoms in each complex and are found to correlate with the reported anticancer activities of these complexes. The networks of charge transfer and HB significantly contribute to the observed synergistic effects of the coordinated ligands confirming the possibility of these effects resulting in the high cytotoxicity as reported experimentally [45]. Several bond properties and intraatomic properties have been used to explain the chemistry of these complexes, and the correlation of the computed properties was considered. The QTAIM analysis shows that all the metal–ligand bonds are closed shell interactions characterised with dative bonds. The HB interactions in the complexes increased when the benzene is substituted with cymene as the arene moiety of the ligand, and the hydrated forms are characterised with stronger HB interactions indicating that HB may have significantly affected their experimentally reported activation by hydration. The intraatomic properties show that only the H atoms are characterised with high and low indicating that they can be easily perturbed than any other atoms in the complexes. There is a little increase in of Ru and P atoms in the hydrated forms.
Acknowledgments
The authors gratefully acknowledge the financial support of Govan Mbeki Research and Development Centre, University of Fort Hare, South Africa. The CHPC and CSIR, Department of Science and Technology, South Africa, is gratefully acknowledged for providing the computing facilities and some of the software that is used for the computation.
Supplementary Materials
Supplementary Tables: Table S1: The QTAIM properties of selected bonds in complex 1, 2, 3 and 4 that are associated with metal using 321G. Table S2: The correlation in the properties of the electron density critical point analysis of complexes 1, 2, 3 and 4. Table S3: Selected atomic properties derived through QTAIM analysis for complex 1, 2, 3 and 4 using 321G. Table S4: The correlation in the atomic properties of complexes 1, 2, 3 and 4. Table S5: The sum total of the computed atomic properties over all the atoms in complexes 1, 2, 3 and 4 from DGDVZP(Ru)631+G(d,p) systems.
References
 E. Meggers, “Targeting proteins with metal complexes,” Chemical Communications, no. 9, pp. 1001–1010, 2009. View at: Publisher Site  Google Scholar
 C. S. Allardyce, A. Dorcier, C. Scolaro, and P. J. Dyson, “Development of organometallic (organotransition metal) pharmaceuticals,” Applied Organometallic Chemistry, vol. 19, no. 1, pp. 1–10, 2005. View at: Publisher Site  Google Scholar
 P. J. Dyson and G. Sava, “Metalbased antitumour drugs in the post genomic era,” Dalton Transactions, no. 16, pp. 1929–1933, 2006. View at: Publisher Site  Google Scholar
 Y. N. V. Gopal, D. Jayaraju, and A. K. Kondapi, “Inhibition of topoisomerase II catalytic activity by two ruthenium compounds: a liganddependent mode of action,” Biochemistry, vol. 38, no. 14, pp. 4382–4388, 1999. View at: Publisher Site  Google Scholar
 C. S. K. Rajapakse, A. Martínez, B. Naoulou et al., “Synthesis, characterization, and in vitro antimalarial and antitumor activity of new ruthenium(II) complexes of chloroquine,” Inorganic Chemistry, vol. 48, no. 3, pp. 1122–1131, 2009. View at: Publisher Site  Google Scholar
 M. Groessl, Y. O. Tsybin, C. G. Hartinger, B. K. Keppler, and P. J. Dyson, “Ruthenium versus platinum: interactions of anticancer metallodrugs with duplex oligonucleotides characterised by electrospray ionisation mass spectrometry,” Journal of Biological Inorganic Chemistry, vol. 15, no. 5, pp. 677–688, 2010. View at: Publisher Site  Google Scholar
 B. Dutta, R. Scopelliti, and K. Severin, “Synthesis, structure, and reactivity of the methoxybridged dimer [Cp∧Ru(μOMe)]2 (Cp∧ = η51methoxy2,4ditertbutyl3neopentylcyclopentadienyl),” Organometallics, vol. 27, no. 3, pp. 423–429, 2008. View at: Publisher Site  Google Scholar
 B. Dutta, C. Scolaro, R. Scopelliti, P. J. Dyson, and K. Severin, “Importance of the πligand: remarkable effect of the cyclopentadienyl ring on the cytotoxicity of ruthenium PTA compounds,” Organometallics, vol. 27, no. 7, pp. 1355–1357, 2008. View at: Publisher Site  Google Scholar
 W. H. Ang and P. J. Dyson, “Classical and nonclassical rutheniumbased anticancer drugs: towards targeted chemotherapy,” European Journal of Inorganic Chemistry, no. 20, pp. 4003–4018, 2006. View at: Publisher Site  Google Scholar
 G. Gasser, I. Ott, and N. MetzlerNolte, “Organometallic anticancer compounds,” Journal of Medicinal Chemistry, vol. 54, no. 1, pp. 3–25, 2011. View at: Publisher Site  Google Scholar
 S. BetanzosLara, L. Salassa, A. Habtemariam et al., “Photoactivatable organometallic pyridyl ruthenium(II) arene complexes,” Organometallics, vol. 31, no. 9, pp. 3466–3479, 2012. View at: Google Scholar
 E. A. Hillard and G. Jaouen, “Bioorganometallics: future trends in drug discovery, analytical chemistry, and catalysis,” Organometallics, vol. 30, no. 1, pp. 20–27, 2011. View at: Publisher Site  Google Scholar
 E. A. Meyer, R. K. Castellano, and F. Diederich, “Interactions with aromatic rings in chemical and biological recognition,” Angewandte Chemie, vol. 42, no. 11, pp. 1210–1250, 2003. View at: Publisher Site  Google Scholar
 I. Alkorta, F. Blanco, and J. Elguero, “Simultaneous interaction of tetrafluoroethene with anions and hydrogenbond donors: a cooperativity study,” Journal of Chemical Theory and Computation, vol. 5, no. 4, pp. 1186–1194, 2009. View at: Publisher Site  Google Scholar
 J. A. Platts, J. Overgaard, C. Jones, B. B. Iversen, and A. Stasch, “First experimental characterization of a nonnuclear attractor in a dimeric magnesium(I) compound,” Journal of Physical Chemistry A, vol. 115, no. 2, pp. 194–200, 2011. View at: Publisher Site  Google Scholar
 I. M. Kapetanovic, “Computeraided drug discovery and development (CADDD): in silicochemicobiological approach,” ChemicoBiological Interactions, vol. 171, no. 2, pp. 165–176, 2008. View at: Publisher Site  Google Scholar
 Y. Duan, C. Wu, S. Chowdhury et al., “A pointcharge force field for molecular mechanics simulations of proteins based on condensedphase quantum mechanical calculations,” Journal of Computational Chemistry, vol. 24, no. 16, pp. 1999–2012, 2003. View at: Publisher Site  Google Scholar
 N. Rai, S. P. Tiwari, and E. J. Maginn, “Force field development for actinyl ions via quantum mechanical calculations: an approach to account for many body solvation effects,” Journal of Physical Chemistry B, vol. 116, pp. 10885–10897, 2012. View at: Google Scholar
 K. Cobb, “Dock this: drug design feeds drug development,” Biomedical Computation Review, pp. 20–30, 2007. View at: Google Scholar
 A. L. Harvey, “Natural products in drug discovery,” Drug Discovery Today, vol. 13, no. 1920, pp. 894–901, 2008. View at: Publisher Site  Google Scholar
 S. Chatterjee, S. Kundu, A. Bhattacharyya, C. G. Hartinger, and P. J. Dyson, “The ruthenium(II)arene compound RAPTAC induces apoptosis in EAC cells through mitochondrial and p53JNK pathways,” Journal of Biological Inorganic Chemistry, vol. 13, no. 7, pp. 1149–1155, 2008. View at: Publisher Site  Google Scholar
 A. E. Egger, C. G. Hartinger, A. K. Renfrew, and P. J. Dyson, “Metabolization of [Ru(η6C6H5CF 3)(pta)Cl2]: a cytotoxic RAPTAtype complex with a strongly electron withdrawing arene ligand,” Journal of Biological Inorganic Chemistry, vol. 15, no. 6, pp. 919–927, 2010. View at: Publisher Site  Google Scholar
 M. Hanif, H. Henke, S. M. Meier et al., “Is the reactivity of M(II)arene complexes of 3hydroxy2(1 H)pyridones to biomolecules the anticancer activity determining parameteŕ,” Inorganic Chemistry, vol. 49, no. 17, pp. 7953–7963, 2010. View at: Publisher Site  Google Scholar
 F. Wang, A. Habtemariam, E. P. L. van der Geer et al., “Controlling ligand substitution reactions of organometallic complexes: tuning cancer cell cytotoxicity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 51, pp. 18269–18274, 2005. View at: Publisher Site  Google Scholar
 A. GarcíaFernández, J. Díez, M. P. Gamasa, and E. Lastra, “Facile modification of 1,3,5triaza7phosphatricyclo[3.3.1.1 3,7]decane phosphanes coordinated to ruthenium(II),” Inorganic Chemistry, vol. 48, no. 6, pp. 2471–2481, 2009. View at: Publisher Site  Google Scholar
 W. H. Ang, E. Daldini, C. Scolaro, R. Scopelliti, L. JuilleratJeannerat, and P. J. Dyson, “Development of organometallic rutheniumarene anticancer drugs that resist hydrolysis,” Inorganic Chemistry, vol. 45, no. 22, pp. 9006–9013, 2006. View at: Publisher Site  Google Scholar
 C. Adamo and V. Barone, “Toward reliable density functional methods without adjustable parameters: the PBE0 model,” Journal of Chemical Physics, vol. 110, no. 13, pp. 6158–6170, 1999. View at: Google Scholar
 W. J. Stevens, M. Krauss, H. Basch, and P. G. Jasien, “Relativistic compact effective potentials and efficient, sharedexponent basis sets for the third, fourth, and fifthrow atoms,” Canadian Journal of Chemistry, vol. 70, pp. 612–630, 1992. View at: Google Scholar
 D. Feller, “The role of databases in support of computational chemistry calculations,” Journal of Computational Chemistry, vol. 17, no. 13, pp. 1571–1586, 1996. View at: Google Scholar
 K. L. Schuchardt, B. T. Didier, T. Elsethagen et al., “Basis set exchange: a community database for computational sciences,” Journal of Chemical Information and Modeling, vol. 47, no. 3, pp. 1045–1052, 2007. View at: Publisher Site  Google Scholar
 R. Marchal, P. Carbonnière, D. Begue, and C. Pouchan, “Structural and vibrational determination of small galliumarsenide clusters from CCSD(T) and DFT calculations,” Chemical Physics Letters, vol. 453, no. 1–3, pp. 49–54, 2008. View at: Publisher Site  Google Scholar
 R. Marchal, P. Carbonnière, and C. Pouchan, “Structural and vibrational properties prediction of SnnTen clusters (n = 28) using the GSAM approach,” Computational and Theoretical Chemistry, vol. 990, pp. 100–105, 2012. View at: Google Scholar
 A. D. Becke, “Densityfunctional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993. View at: Google Scholar
 K. D. Dobbs and W. J. Hehre, “Molecular orbital theory of the properties of inorganic and organometallic compounds 5. Extended basis sets for firstrow transition metals,” Journal of Computational Chemistry, vol. 6, pp. 861–879, 1987. View at: Google Scholar
 A. A. Granovsky, “Firefly version 7. 1.G,” http://classic.chem.msu.su/gran/firefly/index.html. View at: Google Scholar
 M. W. Schmidt, K. K. Baldridge, J. A. Boatz et al., “General atomic and molecular electronic structure system,” Journal of Computational Chemistry, vol. 14, no. 11, pp. 1347–1363, 1993. View at: Google Scholar
 M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., G03a: GAUSSIAN03, Revision D. 01, Gaussian, Wallingford, UK, 2009.
 T. A. Keith, “AIMAll (Version 12. 06. 03),” Overland Park, Kan, USA, TK Gristmill Software, 2012, http://aim.tkgristmill.com. View at: Google Scholar
 A. E. Reed, L. A. Curtiss, and F. Weinhold, “Intermolecular interactions from a natural bond orbital, donoracceptor viewpoint,” Chemical Reviews, vol. 88, no. 6, pp. 899–926, 1988. View at: Google Scholar
 J. P. Foster and F. Weinhold, “Natural hybrid orbitals,” Journal of the American Chemical Society, vol. 102, no. 24, pp. 7211–7218, 1980. View at: Google Scholar
 E. D. Glendening and A. Streitwieser, “Natural energy decomposition analysis: an energy partitioning procedure for molecular interactions with application to weak hydrogen bonding, strong ionic, and moderate donoracceptor interactions,” The Journal of Chemical Physics, vol. 100, no. 4, pp. 2900–2909, 1994. View at: Google Scholar
 Z. Chval, Z. Futera, and J. V. Burda, “Comparison of hydration reactions for “pianostool” RAPTAB and [Ru(η6 Arene)(en)Cl]+ complexes: density functional theory computational study,” Journal of Chemical Physics, vol. 134, no. 2, Article ID 024520, 2011. View at: Publisher Site  Google Scholar
 I. C. de Silva, R. M. de Silva, and K. M. N. de Silva, “Investigations of nonlinear optical (NLO) properties of Fe, Ru and Os organometallic complexes using high accuracy density functional theory (DFT) calculations,” Journal of Molecular Structure, vol. 728, no. 1–3, pp. 141–145, 2005. View at: Publisher Site  Google Scholar
 P. J. Mendes, J. P. Prates Ramalho, A. J. E. Candeias, M. P. Robalo, and M. H. Garcia, “Density functional theory calculations on η5 monocyclopentadienylnitrilecobalt complexes concerning their secondorder nonlinear optical properties,” Journal of Molecular Structure, vol. 729, no. 12, pp. 109–113, 2005. View at: Publisher Site  Google Scholar
 I. N. Stepanenko, A. Casini, F. Edafe et al., “Conjugation of organoruthenium(ii) 3(1hbenzimidazol2 yl)pyrazolo[3,4b]pyridines and indolo[3,2d]benzazepines to recombinant human serum albumin: a strategy to enhance cytotoxicity in cancer cells,” Inorganic Chemistry, vol. 50, no. 24, pp. 12669–12679, 2011. View at: Publisher Site  Google Scholar
 M. J. Calhorda and P. E. M. Lopes, “An “atoms in molecules” (AIM) analysis of the dihydrogen bond in organometallic compounds,” Journal of Organometallic Chemistry, vol. 609, no. 12, pp. 53–59, 2000. View at: Google Scholar
 D. Tiana, Organometallic Chemistry from the Interacting Quantum Atoms Approach [Ph.D. thesis], Università degli Studi di Milano, Milan, Italy, 2009/2010.
 I. Rozas, I. Alkorta, and J. Elguero, “Behavior of ylides containing N, O, and C atoms as hydrogen bond acceptors,” Journal of the American Chemical Society, vol. 122, no. 45, pp. 11154–11161, 2000. View at: Publisher Site  Google Scholar
 I. Alkorta, I. Rozas, and J. Elguero, “Nonconventional hydrogen bonds,” Chemical Society Reviews, vol. 27, no. 2, pp. 163–170, 1998. View at: Google Scholar
 S. Hammerum, “Alkyl radicals as hydrogen bond acceptors: computational evidence,” Journal of the American Chemical Society, vol. 131, no. 24, pp. 8627–8635, 2009. View at: Publisher Site  Google Scholar
 J. E. D. Bene, I. Alkorta, G. SanchezSanz, and J. Elguero, “Variations in the structures and binding energies of binary complexes with HBO,” Chemical Physics Letters, vol. 538, pp. 5–9, 2012. View at: Google Scholar
 K. Tang and F. Q. Shi, “Comparative analysis of blueshifted hydrogen bond versus conventional hydrogen bond in methyl radical complexes,” International Journal of Quantum Chemistry, vol. 107, no. 3, pp. 665–669, 2007. View at: Publisher Site  Google Scholar
 E. S. Shubina, N. V. Belkova, and L. M. Epstein, “Novel types of hydrogen bonding with transition metal πcomplexes and hydrides,” Journal of Organometallic Chemistry, vol. 536537, pp. 17–29, 1997. View at: Google Scholar
 M. Solimannejad, M. Malekani, and I. Alkorta, “Theoretical study of the halogenhydride complexes between XeH2 and carbon halogenated derivatives,” Journal of Molecular Structure, vol. 955, no. 1–3, pp. 140–144, 2010. View at: Publisher Site  Google Scholar
 L. J. Farrugia and H. M. Senn, “On the unusual weak intramolecular C⋯C interactions in Ru3(CO)12: a case of bond path artifacts introduced by the multipole model?” Journal of Physical Chemistry A, vol. 116, no. 1, pp. 738–746, 2012. View at: Publisher Site  Google Scholar
 Q. K. Timerghazin, I. Rizvi, and G. H. Peslherbe, “Can a dipolebound electron form a pseudoatom? An atomsinmolecules study of the hydrated electron,” Journal of Physical Chemistry A, vol. 115, no. 45, pp. 13201–13209, 2011. View at: Publisher Site  Google Scholar
 K. E. Edgecombe, R. O. Esquivel, V. H. Smith Jr., and F. MüllerPlathe, “Pseudoatoms of the electron density,” The Journal of Chemical Physics, vol. 97, no. 4, pp. 2593–2599, 1992. View at: Google Scholar
 A. GarcíaFernández, J. Díez, Á. Manteca, J. Sánchez, M. P. Gamasa, and E. Lastra, “Novel hydridotris(pyrazolyl)borate ruthenium(II) complexes containing the watersoluble phosphane 1,3,5triaza7phosphaadamantane: synthesis and evaluation of DNA binding properties,” Polyhedron, vol. 27, no. 4, pp. 1214–1228, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Adebayo A. Adeniyi and Peter A. Ajibade. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.