Theoretical Determination of Influence of the Metallic State of Oxidation toward Cytotoxic Activity: Case of Ruthenium Complexes ()
1. Introduction
Cancer is known as a widespread and mortal disease. Thus, investigations are currently made to eradicate it. Besides, since discovered, many forms of cancer are known today [1] [2]. They require different types of treatment such as chemotherapy, radiotherapy, surgery, and even their combination through photodynamic therapy (PDT) [3]. Although the tremendous success of researchers improves the technic to eradicate the disease, the patient is still exposed to the side effects. Namely, the patient can be confronted with the loss of healthy cells. He can also undergo nausea thereby losing appetite or losing absolutely an organ. He can also undergo vomiting, hair loss, and damage to the kidney [4].
In chemotherapy, many types of drugs are used. The first one has hitherto been the cisplatin molecule and its derivatives [5] [6] [7] [8] [9]. Since these molecules were assumed to present more side effects although their efficiency, other molecules are explored. We can number natural molecules like coumarins [10] [11] and metallic complexes [12]. Regarding these metallic complexes, ruthenium complexes are of great interest. Their particularity is that they have many oxidation states and can easily make difference between ill and healthy cells thereby reducing drastically the toxic effects [13]. They are also of great interest because of their capacity to mimic the transferrin when they are bonded to the biological molecules in the blood [14]. The first complexes discovered are ruthenium complexes of bipyridine, azopyridine, and others like NAMI.
The most common fundamental oxidation states of ruthenium known within the complex molecules are states II and III. Both are assumed to react on the DNA by binding covalently to the N7 position of guanin to form 1,2- or 1,3-intrastrand crosslinks, and interstrand crosslinks. Sometimes, molecules that display halogen atoms in their structures are assumed to undergo hydrolysis and the loss of those atoms that favor the linking to the DNA of the infected cells [15]. Yet, certain complexes do not present any halogen atom and they are known for their efficiency and stability as drugs against Cancer diseases [16] [17] [18] [19]. Therefore, how can we explain this discrepancy between those ruthenium complexes; however, they all are discovered active molecules at fundamental state?
Through this work, we intend to find out the key element that renders the ruthenium complexes so active cytotoxically regardless their states of oxidation and classify them according to their reactivity. We have focused our research on a series of ruthenium complexes that display states Ru(II) and Ru(III) as fundamental oxidation states. We have complexes α-RuCl2(Azpy)2 [20] and RuCl3(H2O)2Gly [21], RuCl3Terpyridine [22], [RuCl4(Him)2]− (NAMI), [RuCl4(Him)(DMSO)]−(NAMI-A) [23],
,
and
[24] [25]. Respectively, Gly stands for glycoside, Him for imidazoline, Azpy for 2-phenylazopyridine and Bipy for bipyridine. Their first difference that can be noticed is through their ligands that are different from one complex to another. Besides, before pursuing the investigation of the discrepancy that regards their state, we will summarize the description of each complex.
1.1. RuCl3(H2O)2Gly
This molecule was studied by De Sousa et al. [21]. it is considered as a +III state molecule. This molecule is assumed to possess antileishmanial activity [26]. As it presents three isomers, the fac-cis-diaquo-trischloro-glycinatoruthenium III isomer presented in Figure 1 displays the medical activity.
1.2. α-RuCl2(Azpy)2
The structure of this molecule is presented in Figure 2. It has much been studied since the discovery of its biological activity [15] [27] [28] [29] [30]. Although the complex shows five isomers due to the antisymetrical shape of the ligand, the cis-RuCl isomer was tested to display cytotoxic activities against renal cancer
(A498), lung cancer (H226), ovarian cancer (IGROV), breast cancer (MCF-7) and colon cancer (WIDR). The authors assumed that this complex presents the ruthenium atom in its +II state.
1.3. RuCl3 Terpyridine
Figure 3 shows the structure of this molecule. It was discovered by Paul M. Van Vliet et al. [22], it presents also two isomers where the mer isomer is admitted as biologically active molecule against leukemia cells. Presenting the Ru atom in +III state, the molecule is particularly assumed to bind to DNA via the guanine [31].
1.4. NAMI and NAMI-A
Both molecules are anti-cancer drugs. Their activities are noticed by the substitution of the chlorine atoms through hydrolysis. They were both discovered after noticing the side effects provoked by cis-platin and its derivatives [32] [33]. Even if they display four chloride atoms in their structure, both complexes are considered with the metal Ru in its +III state [34] [35]. Figure 4 shows both structures to highlight their difference.
Figure 4. Structure of both NAMI and NAMI-A.
1.5.
and
Both complexes are isomers. They show the Ru atom in the state of +II. Among them, the merisomer is assumed to be cytotoxic agent against human tumors [36]. What is more, this isomer was especially assumed to be very active against M19 cells more than the cisplatin [24]. However, the Fac isomer that is C3 symetrical is reported to appear unlikely [37] inactive owing to the hindrance created by the orientation of the three azpy ligands.
1.6.
This complex is assumed to fix the ruthenium atom at +II state of oxidation. It is a D3 symmetrical and its cytotoxicity was assumed to depend on the metal, the ligand and the shape of the molecule [38]. It is admitted hampering the acetylcholinesterase, which is an enzyme mainly involved in neurotransmittermetabolism [39]. Figure 5 displays the both isomers
and
molecules.
Therefore, we will theoretically deal with the difference between both types of ruthenium complexes. First and foremost, the Ru(II) are assumed to be more active than R(III) that are biologically inert [40].
2. Methods
We performed All calculations by using the Gaussian program package 09 [41]. All the structures were optimized with density functional theory methods at Wb97xd/Lanl2dz and B3lyp/lanl2dz levels. The pseudo potential Lanl2dz was used for Ru atom which is a heavy atom and displays many inner electrons. These inner electrons are at the origin of relativity that cannot be handled by standard basis set. The calculation of the frequencies shows a lack of imaginaries values thereby confirming the effective stability of the molecules. The Time Dependent Density Functional Theory (TD-DFT) and the Natural Bond Orbital (NBO) calculations were performed at the Wb97xd/Lanl2dz and B3lyp/lanl2dz levels.
The NBO was preformed to find out the charge and the electronic structure of central metal ruthenium and TD-DFT consisted in performing the spectrum of each complex.
Figure 5. Structure of
,
and
.
3. Results and Discussion
3.1. Charge of the Ruthenium in the Complex
All the studied complexes were assumed to be synthesized through RuCl3·3H2O. yet, they present two types of state according to literature. For the state +II of the ruthenium named Ru(II), we have the complexes α-RuCl2(Azpy)2,
,
and
. Regarding the +III state, we have the complexes NAMI, NAMI-A, RuCl3Terpy and RuCl3(H2O)2Gly including RuCl3·3H2O. They are Ru(III) complexes. Table 1 displays the charge calculated through NBO calculation and the electronic valence structure of Ru. We can notice that the values showed by Ru in each complex cannot make difference between the state of the complexes. Here, we can see that in the first Ru(II) group calculated with wb97xd method, the charge of Ru is comprised between 0 to 0.26. For instance, in both Mer and Fac
the charge is the same indicating that the stereoisomer has no effect on the ruthenium charge. Thus, there is not any difference between both isomers regarding the charge of Ru. Whereas α-RuCl2(Azpy)2, shows a negative charge for Ru certainly due to both chloride atoms and since both Mer and Fac
derives from it. It means that the substitution of both chloride atoms by a bidentate Azpy ligand increases the ruthenium’s charge. Regarding B3lyp method, we can notice that both Mer and Fac
have still the same charge confirming the lack of difference between them. Whereas the α-RuCl2(Azpy)2 at B3lyp/Lanl2dz level, it shows that the presence of Cl atoms permits to increase the charge of Ru.
Concerning the second Ru(III) group, we can see in NAMI that Ru has a very small charge of 0.02. However, NAMI-A presents a negative charge of −0.21 and RuCl3(Terpy) shows 0.13 as charge of Ru. Besides, RuCl3·3H2O presents 0.31 as Ru charge that is smaller than the charge displayed in RuCl3(H2O)2Gly (0.34). Yet, the RuCl3·3H2O is also known to present a state +III. Therefore, we can assume that the ruthenium charge does not predict the oxidation state of the
Table 1. Charge and natural valence structure of ruthenium calculated at Wb97xd/Lanl2dz and B3lyp/Lanl2dz levels.
complex for Wb97xd method. At B3lyp/Lanl2dz level for the Ru(III) complexes, we can see that Ru charge is high for NAMI, RuCl3Terpy and RuCl3(H2O)3 with respectively 0.70, 0.73 and 0.90. However, NAMI-A and RuCl3(H2O)2Gly present low charges for Ru that are 0.54 and 0.30 respectively. These charges are lower than those of certain Ru(II) complexes such as
with 0.62.
Regarding the electronic structure of Ru, we recall that with the NBO program, there are the Lewis orbitals and the Rydberg orbitals. The occupied orbitals are known as Lewis’ orbitals and the unoccupied ones are called Rydberg’s orbitals. All the complexes show the same structure concerning both groups of molecules except for α-RuCl2(Azpy)2. Thus, this fact cannot also explain any difference between both groups of state. Therefore, let us investigate the state of the atomic orbital d of the ruthenium in each complex. Table 2 displays the structure of the valence layer of Ru in the complexes calculated at Wb97xd/Lanl2dz level. The structure here regards only the Lewis atomic orbitals.
In octahedral metallic complexes, the structure of valence orbitals in the metal is usually of the type (n − 1)dxnsynpz where x, y and z stand for the numbers of electrons respectively for Atomic Orbitals (AO)d, s and p. The five atomic d orbitals are degenerated in the isolated metal Ru. However, in a complex, according to the crystallin field theory, the five atomic orbitals are divided into two groups: t2gspans AO dxy, dxzand dyz with equal energy and eg comprises dx²-y²and
Table 2. Electronic structure of valence orbitals of ruthenium in the complex calculated at Wb97xd/Lanl2de level.
dz²also degenerated as d orbitals. One of these groups is assumed to theoretically contribute to the complex formation. In the case of octahedral complexes, t2g is assumed to participate in the molecular orbitals [42]. In Table 2, according to NBO calculation, the structure of the valence orbital of Ru is 4dx5sy5pz and the five atomic orbitals in the subshell d of Ru are dxy, dxz, dyz, dx²-y², dz² do not display the same energy nor are they divided in two groups. Therefore, they have not the equal number of electrons. Yet, in this study we can notice that the complexes are different regardless of their octahedral structures. Hens, each complex presents its own structure. Therefore, we can conclude that the state of Ru does not depend on the electronic structure of the d orbitals. For instance, regarding the group Ru(II), while 4dx²-y² is the most stable in α-RuCl2(Azpy)2 it is replaced by 4dz² in
and in
. The particularity of these complexes is that all the five dorbitals have the lowest energy. However in the Ru(III) group, we can see that they have high energy compared to the Ru(II) complexes. Also, the AOs d are randomly involved in the formation of the molecular orbital with the ligands inasmuch as there are no longer degenerated. Besides, Table 3 displays the structure of the valence orbitals of Ru calculated at B3lyp/Lanl2dz level.
Table 3. Electronic structure of valence orbitals of ruthenium in the complex calculated at B3lyp/Lanl2dz level.
Table 3 shows that the structure of ruthenium’s d orbitals does not depend on the method used for calculation. In fact at B3lyp/Lanl2dz level, regarding the Ru(II) group, we can notice that the five d orbitals have not the same energy and they don’t also respect the crystallin field theory. Moreover, with the Ru(III) group, the same remark is made and the atomic orbitals of Ru are randomly involved in the complexes. Furthermore, no particular difference can be observed between both groups apart from the complexes of Ru(II) group seem to display lower energy. Table 4 displays the atomic orbitals d of Ru classified after the remove of the degeneracy.
Table 4 shows that both functional describe the complexes almost the same way. This result shows the Jahn-Teller effect where α-RuCl2(Azpy)2 undergo a compression as 4dx²-y² is more stable than 4dz². However, all the other complexes irrespective of their state display a structure with elongation with the inversion of the stability between both orbitals [42]. Also, this statement cannot explain any difference between both groups Ru(II) and Ru(III). Moreover, Table 5 presents the reactivity parameters relative to frontier orbitals and the properties of the wavelength recorded at both B3lyp/Lanl2dz and Wb97xd/Lanl2dz levels.
Table 5 shows that the gap energy ΔEH-Ldisplays no difference between the
Table 4. Classification of d Atomic Orbitals of Ru after removing the degeneracy.
Table 5. Reactivity parameters with frontier orbitals and absorption properties of Ru(II) and Ru(III) complexes calculated at both Wb97xd/Lanl2dz and B3lyp/Lanl2dz levels, the properties of the wavelength and its corresponding energy and the frequency.
Ru(II) and Ru(III) complexes both with B3lyp/Lan2dz and wb97xd/Lanl2dz levels. At wb97xd/Lanl2dz level we can see that RuCl3Terpy and α-Rucl2(azpy)2 have the smallest values regarding the gap energy. Moreover, both isomers
and
have almost the same gap energy regardless the difference in their geometrical structure. Furthermore, the same remark is observed at B3lyp/Lan2dz level. Among all the complexes, the Ru(III) compound RuCl3Terpy is the most reactive with 1.72 eV as gap energy followed by Ru(II) complex α-Rucl2(azpy)2 with 2.22 eV. Thus, we can assume that the energy gap cannot make any real difference between Ru(II) and Ru(III) complexes regarding their reactivities. Besides both molecules, the third most reactive complex at both levels is the originate RuCl3(H2O)3. This fact shows that the reactivity of Ru is the most effective when the ligands are Azpy and Terpy in presence of Cl atoms.
Regarding the absorption properties, we can see again that only RuCl3Terpy and α-Rucl2(azpy)2 present wavelengths in the visible domain up to 500 nm. This data is characteristic of MLCT transfer [43]. They present respectively 521.35 nm and 624.06 nm at Wb97xd/Lanl2dz level and 536.41 nm and 640.44 nm at B3lyp/Lanl2dz level. Moreover, they require the least energies to perform these metals to ligand charge transfer transitions. Basically, B3lyp offers all the complexes to be active in visible domain regardless the state of the metal.
Figure 6 displays the energetic diagram of the Ru complexes recorded at Wb97xd/Lanl2dz level. We can notice that the energy of the frontier orbitals increases with the number of chloride atom of the structure. And comparatively to the energy of the HOMO of the DNA (−2.08) [44] that is assumed to provide electron to the complex sensitizer, Ru(III) complexes are found effectively more capable of binding to the target. Regarding NAMI and NAMI-A, Even if they
Figure 6. Energetic diagramm of Ru complexes calculated at wb97xd/Lanl2dz level. Comparatively to the energy HOMO of the DNA we can effectively see that only Ru(III) complexes are likely available to link to the DNA.
have their LUMO orbitalsrelatively higher than the HOMO of the DNA, they nevertheless are less efficient comparatively to the other Ru(III) complexes’ LUMO. We can assume that this is certainly due to the presence of four chloride atoms. Concerning the Ru(II) complexes, α-RuCl2(Azpy)2 presents the most active LUMO that can react with the cancer molecule thereby confirming its activity as cytotoxic molecule [45]. Regarding the cation molecules such as
,
and
, they have the lowest energy concerning both their HOMO and LUMO orbitals. However, they are far from reacting with DNA as receivers of electron. Their cytotoxicity must certainly be explained by other mechanisms.
When we consider B3lyp/Lanl2dz level, we can notice the same order of energy of the frontier orbitals of the complexes. Regarding the energy of the LUMO orbitals, they are classified according to the following order: LUMONAMI > LUMONAMI-A >
>
>
>
>
>
>
. Herein, Figure 3 displays the energy diagram of Ru complexes performed à B3lyp/Lanl2dz.
Figure 7 shows that
and
have almost the same energy regarding both the HOMO and the LUMO. Yet,
is assumed to be active cytotoxically and
is however inactive. The difference between both complexes may certainly roots in relation of the complexes. Besides, Table 6 displays the structure of the frontier orbitals.
Figure 7. Energy diagram of Ru complexes calculated at B3lyp/Lanl2dz level.
Table 6. Frontier orbitals and their compositions calculated at Wb97xd/Lanl2dz and B3lyp/Lanl2dz levels.
When we consider the groups of molecules involved in the structure of frontier orbitals through Table 6, we can notice that the HOMO orbitals are mainly represented by orbitals of Ru in all the complexes. However, regarding the LUMO orbitals for Ru(II) complexes, we notice in this case that Ru does not intervene in the orbital both with Wb97xd and B3lyp functionals. Here again we cannot see any difference that can explain the inability for
to react as therapy molecular agent comparatively with
. Therefore, let us study the photophysical properties of the complexes.
3.2. TDDFT Analysis of the Complexes
3.2.1. Singlet Excited States
Table 7 displays the lowest singlets of Ru(II) and Ru(III) calculated at B3lyp/Lanl2dz level.
Table 7. Lowest singlets and triplets calculated at B3lyp/Lanl2dz of Ru(II) and Ru(III) complexes. It displays their energy (eV) the wavelength (nm) and the oscillator strength (f).
We can see in Table 7 that almost all the Ru(II) complexes present their first wavelengths higher than 600 nm thereby indication the capacity to behave as photo therapy agent. In fact, they wavelengths are situated in the therapeutic window [46]. However, only
presents 489.43 nm suitable for shallow surface diseases [47].
Regarding the Ru(III) molecules, NAMI and RuCl3(H2O)2Gly are active only for superficial diseases with respectively 481.72 nm and 478.97 nm. However, NAMI-A, RuCl3Terpy with RuCl3(H2O)3 show their ability to act against deeply localized cancers.
3.2.2. Triplets Excited States
This is the most important state for the complexes since its lifetime is very longer than the singlet’s lifetime [48]. And it is responsible for the photosensitized reactions [49]. Table 7 displays the six lowest excited triplet states. We can notice that the lowest energy T1 (ET1) shows a value of a frequency different from zero apart from the case of α-RuCl2(Azpy)2,
and
. These molecules belong all to the family of Ru(II). It results from this statement that those three complexes cannot generate a T state. Therefore, they cannot be used for photodynamic therapy. However,
that is reported to not display any cytotoxic activity appears the lone to show an excited state with 0.085 eV as ET1 and its wavelength at 14,542.30 nm. Regarding the Ru(III) complexes, we can see that all of them display T energies with high frequency. Therefore, they can well be used as PDT agents. Then, our investigation will concern only these Ru(III) complexes and
. Thus, let’s determine the vertical electron affinity and vertical potential ionization both at ground and excited states.
3.2.3. Vertical Electron Affinity (VEA) and Vertical Ionization Potential (VIP)
Table 8 and Table 9 display the vertical electron affinity (VEA) and vertical ionization potential (VIP) both at ground state S0 and excited T1 state. There are assumed to be respectively at the origin of electron donation and electron acceptance in reactions. We can see for both levels of calculus that VEAS0 is negative for all the complexes except for NAMI-A. Therefore, the coming work will be performed only with B3lyp/Lanl2dz. Here, NAMI appears to be the most ionizable complex with the lowest VIPS0 (0.15 hartree) followed by NAMI-A. the complexes without chloride atoms are the least ionizable. It may indicate that the strength of ionization depends on the numbers of chloride atoms in the structure. Concerning VEAS0, it shows that the molecules most able to accept electrons are the cations. And the anions with four chloride atoms namely NAMI and NAMI-A are the least acceptors. Besides, both
and
don’t show any difference.
Moreover, Table 8 displays the vertical electronic affinity and vertical ionization potential in the excited state T1 for the Ru complexes. As we said before that α-RuCl2(Azpy)2,
and
do not show any triplet
Table 8. Vertical electronic affinity and vertical ionization potential recorded at ground state S0 both at B3lyp/Lanl2dz and Wb97xd/Lanl2dz levels. Ep, Ea and Ec are calculated in hartree.
aVEAS0 = Ea − Ep. bVIPS0 = Ec − Ep.
Table 9. Excitation energies of the lowest triplets (ET1 in eV), vertical electron affinities (VEA in eV) and vertical ionization potentials (VIP in eV) of Ru complexes in vacuum.
aVEAT1 = VEAS0 − ET1; bVIPT1 = VIPS0 − ET1.
state, they therefore are assumed not to be active photo dynamically. So, we will analyze the Ru(III) molecules and
in so far that they show energy values at triplet state. VEAT1 and VIPT1 designate respectively the ability to receive and to donate an electron at the excited triplet state. In Table 8, we can see that
is the most ionizable molecule with VIPT1 = 0.382 a.u. conversely, RuCl3(H2O)2Gly is the most able to receive electron at excited state as it shows the smallest value of VEAT1.
3.3. Study of Photosensitization Capacity of Ru Complexes
The objective here is to find out the mechanism of the reaction of the photosensitized complex with the DNA or the RNA bases which are responsible for the damage of the cell. Two mechanisms known as I and II are hitherto admitted in literature [50]. Through mechanism I, the complex is said to act directly on the DNA or RNA base that is assumed to donate the electron. Also, this reaction s assumed to take place in an area without oxygen molecule. In this case, the complex must be either in its excited state or as a cation. In our previous article, we calculated the VIPS0 and VEAS0 of the DNA and RNA base and they are displayed in Table 10 [49]. Regarding the mechanism II, it is performed in presence of a triplet 3O2 that will be excited into singlet 1O2, responsible for DNA or RNA damage.
3.3.1. Mechanism I
It was admitted that Guanine and Uracil are respectively the most donor and acceptor of electron. It means that those bases must be the targets of the complex in the cell-free of oxygen. Therefore, the prediction of the mechanism must be performed only with both molecules. The reactions concerned by the mechanism I are expressed by the following equations. These equations resume two pathways.
(1)
(2)
(3)
(4)
Equation (1) indicates that Ru complex (Ru-K) being in its excited state T1 will extract electron from the DNA base. This reaction will be governed by VEAT1 of Ru-K and VIP of the DNA base and the sum of both data must be negative to consider the reaction as thermodynamically feasible.
Secondly, a radical cation of Ru-K performed through Equations (2) and (3) can also receive electron from the DNA base. That reaction is expressed by Equation (4). Equation (2) is also expressed by sum either VEAT1 + VIPS0 or VIPT1 + VEAS0 of Ru-K. whereas Equation (3) is governed by the sum VEAT1 + VIPT1.
Table 11 shows that Equation (1) is positive for
and all the Ru(III) complexes. It indicates that the DNA or RNA cells cannot react directly
Table 10. Vertical ionization potential and vertical electronic affinity of DNA and RNA bases performed at B3lyp/Lanl2dz level.
Table 11. Sum of the VEAT1 parameters of the complex and VIP parameters of the Guanine reflecting the interactions according to Equation (1), Sum VEAT1 + VIPS0 or VIPT1 + VEAS0 both reflecting the auto ionization reaction and VIPT1 + VEAT1, all performed at B3lyp/Lanl2dz level.
on the complexes since the sum of vertical electronic affinity of the complexes and the potential ionized vertical of the cells are positive. Whereas the auto ionization Equations (2) and (3), the
still show positive value indicating its impossibility to react without oxygen as photo dynamic therapy agent directly over DNA or RNA molecules. Even the cationic formation of
shows the positive value. We can conclude that the Ru(II) complex cannot react on DNA or RNA by the reaction type I.
Regarding all Ru(III) complexes, we can notice that they all are active through the auto-ionization process both through Equation (2) and Equation (3). We can therefore assume that the reaction type I is allowed for the Ru(III) complexes only if they are radical cations
.
3.3.2. Mechanism II
In this case, the Ru complex must be in excited triplet state. And we can have two reactions: The first reaction regards the Ru complex in an excited state where, it will transfer energy to the triplet oxygen (in its fundamental state) allowing it to become a singlet molecule according to Equation (5). The singlet oxygen owing to its strength is therefore responsible for the damage of all anormal molecules in its environment [51].
(5)
Besides, Equation (5) requires that the energy necessary to excite the complex must be higher than the energy admitted for oxygen to pass from ground state to its excited state. According to literature, the energy necessary to excite 3O2 is 1.06 eV [52] [53]. Furthermore, all the Ru(III) complexes display an ET1 according to Table 11 higher than this energy thereby concluding their ability to excite the oxygen triplet. However,
shows 0.08 eV as ET1 energy. Thus, we can also assume that the Ru(II) complex cannot generate 1O2 molecule.
Table 12. Sum of VIPT1 of Ru complex and AEA(O2) of oxygen in eV.
The second reaction regards a transfer of electron from Ru complex to 3O2 to generate a superoxide radical anion
through Equation (6).
(6)
The condition to perform this reaction requires that the sum of the vertical ionization potential at the triplet state VIPT1 of the complex and the adiabatic affinity of oxygen AEA(O2) must be negative. And the AEA(O2) in vacuum is assumed to be −0.59 eV [46] [54]. Table 12 shows the sum of VIPT1 of Ru complex and AEA(O2).
All the values in Table 12 regarding the sum VIPT1 + AEA(O2) are negative indicating that Ru complexes can generate the radical superoxide anion. Even here,
allows this reaction but with difficulty comparatively to Ru(III) complexes. We can conclude in the mechanism II that
is active only when it concerns the transfer of electron to 3O2 to generate
. Whereas Ru(III) complexes in general, they can be used as photodynamic agent regardless the nature of the mechanism.
4. Conclusions
We have been comparing the reactivity of Ru(II) and Ru(III) complexes through their oxidation state in this paper. We used the DFT method and both functional B3lyp and Wb97xd with the pseudopotential Lanl2dz for this investigation. Both groups of Ru molecules comprise α-RuCl2(Azpy)2,
,
and
for the group Ru(II) and NAMI, NAMI-A, RuCl3Terpy, RuCl3(H2O)2Gly and RuCl3(H2O)3 for the group Ru(III). All these molecules except RuCl3(H2O)3 are assumed to be therapeutic agents at their fundamental state. Moreover, they don’t seem to show any difference regarding the electronic structure of the ruthenium atom irrespective of their oxidation state. In addition, the analysis of the Ru five d orbitals of the complexes shows that they undergo the Jahn-Teller effect that doesn’t affect their therapeutic activity. Only the structure of α-RuCl2(Azpy)2 is discovered different from the others. Besides, regarding their frontier orbitals, they show that the cations display the lowest energy and the energy is increasing with the number of chloride atoms. Also, RuCl3terpy and α-RuCl2(Azpy)2, display the lowest LUMO energy suitable for DNA and RNA binding without excitation however they belong to different groups. Hence, there is no real difference between both groups.
When it comes to the PDT activity, the TD-DFT method was performed to determine vertical electron affinity and vertical ionization potential both at the ground and excited states. This calculus shows that Ru(II) complexes are inactive due certainly to their singlet state where valence electrons are coupled. However,
that is assumed not to be active over DNA or RNA cells at fundamental state is the lone Ru(II) complex able to display the exciting T1 state. Maybe due to the ligand Azpy or its C3 symmetry. Yet, this complex does not form a cation by auto-ionization. And it hardly reacts only over oxygen to generate superoxide radical anion.
Whereas Ru(III) complexes are active both at fundamental and excited states in the presence or lack of oxygen. In the latter state, they can react through auto-ionization to generate a cationic complex without oxygen, or they can generate the radical superoxide
. All these forms can destroy the affected DNA cells. We can retain that the ability for Ru(III) complexes to react better than Ru(II) complexes is certainly due to the presence of the non-coupled electron that does not require much energy to move up to a high shell within Ru(III) complexes. In consequence, for the coming paper, we intend to check whether Ru(II) complex can be active under the effect of light if it is in a triplet charge state. This will certainly be a tremendous breakthrough for therapeutic research.