One of the goals of the present manuscript was the identification and characterization of the synthetized La(2b’)3 coordination complex through a detailed study of the spectroscopic results. In addition, it would be interesting to observe the effect of the La(III) ion on the molecular structure and atomic charge distribution of 2b ligand, particularly in the triazole ring. Different basis sets were tested and the best one was selected.
3.2. Molecular structure of the complex
Since lanthanide metals coordinate more preferably to oxygen atoms than to nitrogen atoms [
14], the starting structure optimized was that with the lanthanum(III) ion coordinated with three 2b ligands through the carboxylate group (COO
−),
Figure 1.
Figure 1.
Labelling of the atoms and plot of the optimized La(2b’)3 structure with 2b: sodium 2-(4-chlorophenyl)-5-(pyrrolidin-1-yl)-2H-1,2,3-triazole-4-carboxylate. Front and lateral view forms were included in the optimized structures at the: (a) M06-2X/Lanl2dz level, and (b) M06-2X/Lanl2mb level. Several bond length values of interest calculated are included in the figure, together with the energy values of the system. 1 AU = 2625.5 kJ/mol.
Figure 1.
Labelling of the atoms and plot of the optimized La(2b’)3 structure with 2b: sodium 2-(4-chlorophenyl)-5-(pyrrolidin-1-yl)-2H-1,2,3-triazole-4-carboxylate. Front and lateral view forms were included in the optimized structures at the: (a) M06-2X/Lanl2dz level, and (b) M06-2X/Lanl2mb level. Several bond length values of interest calculated are included in the figure, together with the energy values of the system. 1 AU = 2625.5 kJ/mol.
Two views of the optimized structure at the M06-2X/lanl2dz level are shown in
Figure 1a, with an almost symmetrical arrangement. By rotation of the carboxylate group around the C
9-C
11 bond, another conformer can be obtained, but is slightly less stable than that plotted in this figure. The notation used for labelling the atoms is in accordance with the one previously reported for this very same 2b ligand [
12]. Few bond length values are also included in this figure. In the bottom of it, the total energy (
E) value is shown, which includes the ZPE (zero-point vibrational energy) correction, and the Gibbs energy (
G). Several selected optimized geometrical parameters, namely bond lengths, bond angles and torsional angles of interest in one of the ligands (labelled as
I) are collected in
Table 4. For comparison purposes, the MP2 values obtained in the isolated 2b structure [
22] are listed in the last column of the table.
The intramolecular H-bond between the carbonyl oxygen O
12 and the pyrrolidine hydrogen H
18 is very weak by M06-2X (2.443 Å) in the La-(2b’)
3 complex due to La-O
12 bond formation, but it is present in the isolated ligand by MP2 (2.034 Å) [
22].
The average calculated value of the La-O bond length (2.5 Å) agrees well with the experimental value of the odd diacids and it is specific for the monometallic complexes of lanthanides [
14,
23]. An increase in this value has been related to a closing of the OCO angle [
24]. The calculated bond angle C=O
12-La (94.0º) is slightly different to the C=O
13-La (95.7º) angle,
Table 4, which indicates that the C-O-La bond is not fully linear due to the slightly higher negative charge of O
13 (−1.266
e), compared to O
12 (−1.250
e),
Table 5. A higher negative charge on O
13 leads to a shortened of the O
13-La bond length (2.491 Å) vs. O
12-La (2.523 Å), and therefore an opening of the C=O
13-La angle. This very slight asymmetry is also observed in the O
12-La-O′
13 angle (103.0º), which is slightly more open than O
12-La-O′
13 (102.0º). These features also lead to rotation of the ligands, and thus the torsional angle C
11-O
12···O′
12-C′
11 (13.2º) is quite different from the C
9-C
11···C′
11-C′
9 angle (−3.2º), as well as the (C
11···La·· C′
11) angle (121.1º), compared to (C′
11···La···C′′
11) (118.3º).
In this optimized structure, at the M06-2X/lanl2dz level, the aryl ring is fully planar and almost coplanar with the triazole ring, with a C5-C4-N4-N10 torsional angle of 0.8º, although it appears slightly rotated (−2.3º) by MP2. The triazole ring is also planar with values of the torsional angles lower than 1º. The pyrrolidine ring, however, is noticeably out-of-plane with a value of the torsional angle N14-C15-C16-C17 of −32.2º, in accordance to that found in a cyclopyrrol ring. This pyrrolidine ring is also non-coplanar with the triazole ring plane with a value of the torsional angle C9-C8-N14-C18 of −18.3º, although it is 32.0º by MP2 indicating a failure in this angle calculation by M06-2X method with lower values.
The carbon C11 atom of the carboxylate group appears slightly titled out of the triazole ring plane and remarkably rotated, but this rotation is larger with O13, than with O12, with a value of the torsional angle N10-C9-C=O13 of −11.2º vs. C8-C9-C=O12 of −7.5º. For the binding to the lanthanide ion, the OCO angle closes up to 117.4º vs. 131.3º by MP2 in the free ligand, and the CO bond lengths are noticeably lengthened. This large flexibility of the carboxylic oxygens facilitates bonding to the lanthanum(III) ion. The lengthening of the CO bonds leads to a shortening of the C9-C11 bond length, which noticeably reduces the double bond character of N7=C8 and C9=N10 bonds of the triazole ring and, consequently, the N4-N10 and C8-N14 bonds are shortened. This feature slightly modifies the bond angles of the triazole ring with the lanthanide binding.
The optimized structure at the M06-2X/lanl2mb level appears deformed,
Figure 1b, perhaps due to its smaller basis set. However, this more planar arrangement can be more probably expected to be present for the solid state sample due to the crystal packing forces compressing the structure. The selected optimized geometric parameters at this level are also included in
Table 4. It is noted that the intramolecular O
12···H
18 H-bond only appears in ligand
I, and it is very weak in ligands
II and
III. The C-N and N-N bond lengths are noticeably larger and the La-O bonds - shorter than by M06-2X/lanl2dz. Because of the structural deformation, the C=O-La and O
12-La-O′
13 angles are closer and the O
12-La-O′
13 angle larger than with the lanl2dz basis set. The pyrrolidine and carboxylate groups are also more rotated related to triazole ring, and the arrangement of the ligands differ remarkably, with a (C′
11···La···C′′
11) angle of 84.7º vs. 118.3º with lanl2mb basis set, and torsional angle (C
11-O
12···O′
12-C′
11) of −169.5º vs. 13.2º with lanl2mb basis set. However, the value of the OCO angle remains almost the same in both levels.
Due to the optimized structure at the B3LYP/Cep-4g appearing largely deformed, it was included in Figure S1 (Supplementary Material). To this feature contributes the low basis set used, as well as that the B3LYP method does not consider long range interactions that stabilize the system to a symmetric arrangement as in
Figure 1a. The optimized values at this level appear with noticeably longer bond lengths than by M06-2X and MP2. With longer bond lengths, the triazole substituents are more twisted, especially the pyrrolidine ring, which can interact with other ligands.
3.3. Atomic charges and molecular properties calculation
The lanthanide ion has the highest positive charge, which is 3.164
e by M06-2X/Lanl2dz and 2.884
e by B3LYP/Cep-4G,
Table 5. Uncommonly, it is rather small (2.066
e) by M06-2X/Lanl2mb. As expected, the highest negative charge corresponds to the oxygen atoms, with the charge on O
13 slightly higher in negative value than on O
12, and in contrast to that calculated in the free 2b molecule by MP2. These oxygens appear as the most reactive and therefore they are bonded to the lanthanide ion. They are expected to participate in the biological activity of these molecules. The negative charge on the nitrogen atoms N
4 and N
14 is lower but it is expected that they can also contribute to the biological activity of this compound. By contrast, a high positive charge appears in the carbon atoms C
8 and C
11 because they are bonded to large negative atoms.
Table 5.
APT (Atomic Polar Tensor) charges calculated at different DFT levels in La-(2b’)3 complex. The values of ligand I are only shown.
Table 5.
APT (Atomic Polar Tensor) charges calculated at different DFT levels in La-(2b’)3 complex. The values of ligand I are only shown.
|
M06-2X |
B3LYP |
MP2 |
atom |
Lanl2dz |
Lanl2MB |
Cep-4G |
6-31G(d,p)a
|
La CL C1 C4 N4 N7 C8 C9 N10 C11 O12 O13 N14
|
3.164 −0.423 0.452 0.518 −0.566 −0.293 0.714 −0.835 0.340 1.821 −1.250 −1.266 −0.924 |
2.066 −0.496 0.476 0.723 −0.832 −0.069 0.465 −0.716 0.398 1.506 −0.925 −1.033 −0.771 |
2.884 −0.254 0.267 0.620 −0.742 −0.066 0.373 −0.844 0.401 1.381 −1.014 −1.089 −0.713 |
- −0.054 −0.080 0.210 −0.073 −0.402 0.419 0.031 −0.240 0.955 −0.887 −0.839 −0.557 |
Several thermodynamic parameters, rotational constants and dipole moments were also calculated and included in
Table 6. Because of the symmetry of the complex at M06-2X/Lanl2dz level, the rotational constant values in the three directions (A, B, C) have almost the same value. They differ, however, by the other DFT methods. One value is larger (A axis) - corresponding to the ligand noted as
I. Two values are lower (B and C axes) - corresponding to ligands
II and
III. The C
v and entropy (S) values were close to M06-2X and lower than those calculated by B3LYP. This feature may be due to the lower symmetry of the complex obtained by B3LYP. The calculated dipole moment value indicates that this complex presents a slight water solubility, which could facilitate its biomedical use.
With the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) values the global chemical reactivity descriptors were calculated [
25,
26] to better understand the reactivity and stability of the La(2b’)
3 complex under study. The energy gap (E
g) between HOMO and LUMO frontier orbitals is a significant characteristic of molecules and helps characterize their chemical reactivity and kinetic stability. A high Eg in a molecule indicates that it is less polarizable. It is generally associated with low chemical reactivity and high kinetic stability. Therefore, our low calculated Eg values in the La(2b’)
3 complex by all three DFT methods indicates large chemical reactivity and small excitation energies to the manifold of excited states.
The following formulae were used to calculate the global chemical reactivity descriptors:
The calculated value for the ionization potential (IP) is somewhat low in accordance with the large reactivity of the complex. The electron affinity (EA) is remarkably lower than IP and positive, with the exception of the M06-2X/Lanl2mb calculation. The electronegativity (χ) is low in accordance with a neutral system. Chemical hardness (η) and global softness (S) express the resistance of a system to a change in its number of electrons. When η is weak the molecule is called soft and it has a small HOMO–LUMO gap and when it is high, the molecule is called hard. Our low values correspond to a soft molecule, with a small gap, and with an electron density that can change easily.
3.4. Vibrational Analysis
For the identification and characterization of the synthesized La(2b’)
3 coordination complex, a detailed analysis of the calculated and experimental spectra was carried out. The wavenumbers with high IR or Raman intensity only were included in
Table 7, while in Table S1 to Table S3 (Supplementary Material) were collected all calculated values corresponding to M06-2X/Lanl2dz, M06-2X/Lanl2mb and B3LYP/Cep-4g levels, respectively. All data correspond to the most stable conformer, with the ligands orientation such as it is plotted in
Figure 1.
The first column lists the calculated wavenumbers at the M06-2X/Lanl2mb or M06-2X/Lanl2dz levels. Three values appear for each vibrational mode corresponding to the three ligands in the complex. Of these three values, the one with higher calculated IR intensity is shown in bold type, while that with higher Raman intensity is shown in italic type. This notation is omitted if these three values are the same, or with almost null IR or Raman intensity. The second and third columns collect the relative IR and Raman intensities (%) calculated in the wavenumber shown in bold type and in italic type, respectively. The relative intensities are obtained by normalizing each calculated value to the intensity of the strongest one in the spectrum. The scaled values by the LSE or PSE procedures are included in another column. The notation used for the values was the same as in the first column. In the next two columns the experimental IR and Raman values observed in the spectra, are listed. The last column corresponds to the main characterization of the vibrations at the M06-2X/Lanl2mb or M06-2X/Lanl2dz levels. In few cases, the % contribution of the different modes to a computed value (PEDs) is included. The calculated values at the B3LYP/Cep-4g level were only included in Table S3 due to their large differences.
The whole IR spectra in the 3750-400 cm
-1 range obtained with the scaled values at the three DFT levels used is shown in
Figure S2. It can be easily compared to the experimental one plotted at the same scale. In a general comparison of these IR spectra the following is observed:
A very broad band at ca. 3500 cm-1 is observed in the experimental spectrum, which can only correspond to the O-H stretching ν(O-H) mode of hydration water strongly H-bonded to the ligands in the La-(2b’)3 complex system. Because the spatial arrangement of the ligands in the complex, many holes appear in the structure that can be occupied by water molecules associated with the synthesis of the complex. Due to the large negative charge around the three carboxylate groups, water molecules can be H-bonded through their hydrogen atoms. Moreover, this band does not appear in the Raman spectrum, as it is expected.
A broad and very strong band centered at 1577.8 cm-1 in the experimental spectrum, in which the in-plane bending δ(O-H) mode of these hydrated water molecules contributes to its broadness.
A large similarity between the scaled spectra at the M06-2X/Lanl2dz and M06-2X/Lanl2mb levels with the experimental one appears, while the B3LYP/Cep-4g level differs remarkably. This feature is in accordance with a more symmetric and better optimized structure by M06-2X method, compared to B3LYP. The best accordance appears at the M06-2X/Lanl2mb level, although the characterization obtained at the M06-2X/Lanl2dz level was also used for the assignment of experimental spectrum.
The coordination of the 2b ligands to the lanthanide ion noticeably changes the IR and Raman spectra. They seem different of those obtained with the 2b ligand molecule alone [
22].
The same comparison but in the whole Raman spectra is included in
Figure S3. Unfortunately, the experimental spectrum shows large background noise that impedes the detection of all weak bands. In this case, the three theoretical scaled spectra appear to have a large similarity to the experimental one. The main difference appears in the experimental Raman line observed at 72 cm
-1, which was not reproduced in the theoretical scaled ones.
For a detailed and better analysis of the different experimental and scaled vibrational wavenumbers of these figures, the spectra are divided in three ranges: from 3750-2600 cm
-1 (
Figure 2), from 1800-1000 cm
-1 (
Figure 3) and from 1000-400 cm
-1 (
Figure 4). Due to a lower quality of the experimental Raman spectrum,
Figure 5 shows only the spectra comparison in the 1800-800 cm
-1 range. The theoretical comparison in other spectral ranges is included as Supplementary Material. Therefore, the region from 3400-2600 cm
-1 is shown in Figure S4, from 1800-1000 cm
-1 in Figure S5, and from 1000-50 cm
-1 in
Figure S6. The assignment of the strong and characteristic Raman lines is included in these figures.
In a general comparison of the scaled spectra, it is noted that most of the calculated normal modes appear in their expected ranges. This feature, together with the fact that most of the scaled strongest vibrations by M06-2X method appear close in its wavenumber to the experimental ones, confirm both the scaling carried out on the calculated wavenumbers, as well as the M06-2X method used. Therefore, the assignments in general could be considered correct.
Because the main objective was the identification and characterization of the synthetized La(2b’)
3 coordination complex, we focus the main attention in the most characteristics and stronger IR and Raman bands to confirm the structure proposed and optimized in
Figure 1. Because the spectra at the B3LYP/Cep-4g level were the worst, their values were omitted in the discussion. For this purpose, the assignment and further discussed was done under the following sections: (i) The COO group modes, (ii) the triazole ring modes, and (iii) the aryl ring modes.
3.4.1. The carboxylate COO group modes
In the isolated state of 2b ligand, the ν
as(COO) stretching mode was predicted (scaled) with very strong IR intensity at 1710 cm
-1 [
22]. However, in the La(2b’)
3 complex it is expected to be significantly red shifted to lower wavenumbers because of the noticeable lengthening of the CO bonds to form the six O-La bonds. At the M06-2X/Lanl2mb level, this stretching mode was clearly identified in the calculated wavenumber with the strongest IR intensity at 1595 cm
-1 (scaled at 1472 cm
-1) in very good accordance to the experimental very strong IR band at 1484.5 cm
-1. However, at the M06-2X/Lanl2dz level, the symmetric ν
s(COO) mode was clearly characterized first at higher wavenumbers, compared to the asymmetric mode. It was at 1596 cm
-1 (scaled at 1516 cm
-1), with strong IR intensity, and at 1607 cm
-1 (scaled at 1527 cm
-1) with very strong Raman intensity. Both scaled values were in good concordance with the experimental strong IR band at 1500.2 cm
-1 and the strong Raman line at 1504 cm
-1. However, this assignment differs from that found at the M06-2X/Lanl2mb level.The scaled spectrum at M06-2X/Lanl2mb level appears to be a slightly better fit with the experimental ones (
Figure 3). The asymmetric stretching vibrations for the carboxylate group (COO
−) are reported to appear as strong IR absorption near 1600-1560 cm
-1 for solid state samples [
28,
29] and at higher wavenumbers than the symmetric ones. For these reasons we assign the experimental IR spectrum mainly to that calculated at the M06-2X/Lanl2mb level.
The symmetric ν
s(COO) stretching mode appears scaled at the M06-2X/Lanl2mb level at 1368 cm
-1 (IR) and at 1372 cm
-1 (Raman) with strong intensity, in excellent agreement with the experimental, very strong IR band at 1372.3 cm
-1 and the Raman line at 1375 cm
-1. These values are also in accordance the ones reported in solid state samples of related compounds [
28,
29] - near the 1420-1400 cm
-1 range for this symmetric stretching mode. At the M06-2X/Lanl2dz level, the asymmetric (instead of the symmetric) mode was the one characterized in calculated wavenumber at 1440 cm
-1 (scaled at 1375 cm
-1), also in good accordance with the experimental spectra at 1372.3 cm
-1 (IR) and 1375 cm
-1 (Raman).
3.4.2. The triazole ring modes
The characteristic normal modes of the 1,2,3-triazole ring have been reported [
30,
31,
32], which are in accordance with our calculations. For simplicity, the discussion was only focused on the assignment of the strongest bands.
NNN modes: The ν
s(NNN) stretching appears strongly coupled with the ν
s(COO) mode, as well as with other ring modes. At the M06-2X/Lanl2mb level, the highest contribution of this mode was identified in the scaled wavenumber at 1375 cm
-1 with strong Raman intensity, in excellent fit with the very strong Raman line at the same value, 1375 cm
-1. A large contribution of this mode was also observed in the scaled wavenumber at 1368 cm
-1, whose major contribution corresponds to the ν
s(COO) stretching mode. At the M06-2X/Lanl2dz level this mode was identified at the same scaled wavenumber, but strongly coupled with the ν
as(COO) mode. In the isolated state of 2b ligand, it was scaled at 1360 cm
-1 and related to the experimental IR band at 1340.5 cm
-1 [
22]. Another triazole ring stretching mode with symmetric character and strongly coupled with an aryl ring mode appears identified in the scaled wavenumber at 1097 cm
-1 and in good relation to the very strong experimental IR band at 1091.2 cm
-1 and to the Raman line at 1091 cm
-1.
The νas(NNN) stretching mode appears characterized at the M06-2X/Lanl2mb level in the scaled wavenumber with very strong IR intensity at 1304 cm-1, also in excellent agreement with the very strong experimental IR band at 1302.0 cm-1.
C8-N14 modes: The stretching mode is predicted at the M06-2X/Lanl2mb level with very high IR intensity and medium Raman activity at 1574 cm
-1 in excellent agreement with the very strong and broad band observed in the IR spectrum and centered at 1577.8 cm
-1, in which this mode contributes in addition to δ(O-H) of the hydrated water molecules. A similar result was obtained at the M06-2X/Lanl2dz level, but it was predicted at 1581 cm
-1 with the highest IR intensity of the spectrum. In the 2b molecule alone [
22], it was scaled at 1556 cm
-1 in accordance to the experimental IR band with medium intensity at 1543.9 cm
-1 and to the Raman line at 1550.6 cm
-1.
3.4.3. The aryl ring modes
The assignments for several aryl ring modes are obvious and require no further discussion, therefore the attention was focused only on the strongest vibrations to confirm the structure of the synthesized complex. For the assignments of the ring modes was followed the Varsanyi notation [
27] for a 1,4-disubstituted benzene.
The aromatic C-H stretching vibrations are generally observed in the 3200-2950 cm-1 range and they are predicted theoretically as almost pure modes (100% PED) with weak and very weak IR intensity and weak Raman activity. Therefore, the only 2b mode with the scaled wavenumber at 2956 cm-1 was related to the experimental IR band at 2968.1 cm-1.
The aromatic C-C stretching vibration, modes 8a and 8b, appears characterized as almost pure modes with %PED higher than 90%, and they are observed in the experimental spectra at similar wavenumbers. However, mode 8a is predicted with almost null IR intensity and with the highest Raman intensity. It is in good agreement with the very strong Raman line observed at 1597 cm-1. Mode 8b is predicted with almost null IR and Raman intensity and thus it was not observed in the experimental spectra.
Mode 19a is predicted at the M06-2X/Lanl2mb level with strong IR and Raman intensity at 1511 cm-1 in good accordance to the experimental strong IR band at 1500.2 cm-1 and to the Raman line 1504 cm-1. The C-C stretching, mode 14, appears predicted at 1234 cm-1 with medium IR intensity and null Raman intensity in accordance with the experimental IR band with medium intensity observed at 1218.0 cm-1 and to its non-detection in the Raman spectrum.
3.4.4. Low-frequency vibrations
Low-lying molecular vibrations have been studied by different authors [
33,
34,
35,
36] and their presence indicates a high flexibility in the molecular structure that may be caused by various factors. Moreover, in our La(2b’)
3 complex the normal vibrational modes that include the lanthanide ion mainly appear below 200 cm
-1. Thus, this kind of vibrations may play an important role in some biological functions, such as transcription and replication of the double DNA helix, specific interactions with proteins and drugs, and gene expression [
36]. Our ligand molecule
2b possesses high conformational flexibility, since it was characterized by ten vibrational modes below 200 cm
-1. These modes may be also important for local recognition in the protein cavity. Moreover, the pattern of individual vibrations for every ligand also reflects the details of its conformation. Taking into account the high deformability of our molecule 2b, these low frequencies should be very sensitive to changes in the local environment, such as hydration and interactions with other biomolecules.
Additionally, the spectra in the frequency region under 200 cm
−1 are quite interesting, since they offer information about metal–ligand vibrations. Thus, in the present work a theoretical analysis of the low-lying vibrations of our La(2b’)
3 complex with frequencies below 200 cm
-1 was carried out. Because the large number of vibrational modes below 200 cm
-1 in our complex, 42 in total, only the six most characteristic ones were considered and plotted in schematic form in
Figure 6. For simplicity and clarity, only the carboxylic and triazole ring vibrations of two of the three ligands were included in this figure. The larger number of these vibrations and its shape, compared to the 2b ligand, agree well with a noticeable increase in the La(2b’)
3 complex flexibility.
The first detail observed in these vibrations is that the lanthanide ion motions do not follow the bond length lines with the carboxylic oxygen atoms. It appears that its motion depends on the modulus and direction of the displacement vectors of the surrounding oxygen atoms. The second detail observed is that most of the vibrational modes follow the displacement vectors calculated in the dimer form of benzoic acid (BA) [
37] and in agreement with those found by far-infrared [
38]. The lanthanide ion does not appear to affect the displacement vectors of the carboxylic group.
As it is expected in the low frequency range, in some atoms the in-plane and out-of-plane character of the displacement vectors is mixed, which makes it difficult to assign it.
Figure 6 shows only the most characteristic modes that can be well related to those corresponding to BA molecule and in which the lanthanide ion has a noticeable large displacement vector. In the remaining vibrational modes, the motions of other groups remarkably prevail over those of the carboxylic group and the lanthanide ion and therefore, they were not included in this
Figure 6. There are also modes that are difficult to be characterized. An analysis of the shape of the selected vibrations reveals that most of the modes involves a strong deformation of complex, and the motion of the whole structure seems, for example as a γ (butterfly) and as a γ (tilting “slanting”).
3.5. Radical-Scavenging assays
It has been reported[
39] that the sodium salt of the ligand (2b), together with its non-ionic, conjugate acid (1b) suppress
in vitro 2-deoxyribose degradation at a concentration of 3·10
-5 M, or higher. Scavenging of DPPH
● was not observed with 2b, mild effect was observed with 1·10
-4 M 1b (RSA=10±2%), decreasing in a concentration-dependent manner to practically zero at 1·10
-5 M. Both compounds scavenge ABTS
●+ at 1·10
-4 M (RSA= 17±1% with 2b and RSA= 13±2% with 1b), again – the effect drops to zero at 1·10
-5 M. Those results will be presented together with the new data, derived from the equivalent experiments on the complex La(2b’)
3.
3.5.1. Impact of 1b, 2b and La(2b’)3 on 2-Deoxyribose degradation
The effect of 2b, and La(2b’)
3 compounds on the UV-induced degradation of 2-deoxyribose are presented in
Figure 7. They were investigated at concentrations between 1·10
-6 M and 1·10
-4 M. A statistically significant concentration-dependent effect with respect to scavenging of OH
● is observed with these three compounds. Unlike 2b and 1b, the complex manifests significant scavenging activity at a lower concentration (SPh-SI= 72±2% at 1·10
-5 M). At concentrations above 3·10
-6 М the scavenging activity of the complex is not different than the scavenging activity of 2b at three times higher concentration. Each La(III) ion coordinates three 2b ligands. This feature may be a clue that all three ligands in the complex participate equally in the scavenging process.
3.5.2. Impact of 1b, 2b, and La(2b’)3 on a model system containing the stable radical DPPH●
The ability of the tested compounds to exchange hydrogen with is presented in
Figure 8. Previous research showed that 2b practically does not interact with DPPH
●, while its conjugate acid 1b interacts only weakly at 1·10
-4 M (RSA= 10±2%), the effect decreasing in a concentration-dependent manner to zero at 1·10
-5 M. For this reason, lower concentrations were not studied. The La(III) complex was investigated at three times lower concentrations, manifesting very mild activity within the tested concentration range. An additional concentration of 1·10
-6 M yielded a result of RSA= 0,5±0,4%. The low activity toward DPPH
● means all three compounds tend to scavenge it via HAT only weakly. One has to consider that DPPH
● has a large molecule and steric hindrances may play a part in the observed low activity. The model system with 2-deoxyribose produces OH
● that are known to be scavenged via HAT reactions [
40]. Unlike DPPH
●, OH
● is small in size and very mobile. La(2b’)
3 and 2b demonstrate clear scavenging effect toward OH
●, generated by water radiolysis at concentrations higher than 1·10
-6 M. The mild, concentration-dependent DPPH-scavenging activity of 1b may be due to the presence of an active carboxyl hydrogen atom in its molecule. La(2b’)
3 is more active than concentrations three times higher of 2b. That may be due to the impact of the La(III) ion on electron density distribution within the coordinated ligands, thus yielding active hydrogen atoms, increasing interaction with DPPH
●.
3.5.3. Impact of 2b, 1b and La(2b’)3 on a model system containing the stable radical ABTS●+
The ability of the tested compounds to participate in SET reactions with ABTS
●+ is presented in
Figure 9. At the highest tested concentration of 1·10
-4 M 2b and 1b have a mild effect on this model system, decreasing to practically statistically zero at 1·10
-5 M. For that reason, lower concentrations were not tested. La(2b’)
3 was investigated at three times lower concentrations. Its activity at 3·10
-5 M (RSA= 9±1%) was greater than that of 2b at the same concentration (RSA= 5±1%). On the other hand, 1·10
-4 M 2b (three times greater concentration) is more active (RSA= 17.3±0,8%) than 3·10
-5 M La(2b’)
3. A similar, statistically significant observation is noted when comparing 1·10
-5 M La(2b’)
3 (RSA= 2.1±1.6%) with 3·10
-5 M 2b (RSA= 5±1%). At concentrations of 1·10
-5 M the ligand and the complex do not demonstrate significant electron exchange with ABTS
●+. La(2b’)
3 tends to be less active than three times more concentrated 2b. That may be due to the fact that the active sites for electron exchange in the ligand might be sterically hindered after complexation with La(III) occurs. Such complexation may also cause electron density redistribution due to the presence of the metal coordination center.
3.5.4. Impact of 2b and La(2b’)3 on MTT-formazan transformation via Fenton reaction derived hydroxyl radicals
The Fenton reaction is a well-known chemical process with clear clinical significance [
41]. “Free” transition metal ions within living organisms create the conditions for its fast progress that yields the highly reactive hydroxyl radicals. OH
● tend to attack molecular sites bearing conjugated double bonds, causing a free radical chain reaction, molecular fragmentation, lipid peroxidation and MDA formation. Antioxidants can stop the aforementioned processes, thus preventing damage to biomolecules and possible pathogenesis.
The water radiolysis model system can serve as a good demonstration of the interaction between OH
● and a potential antioxidant. On the other hand, the Fenton reaction is a significant biochemical process that takes place at physiological conditions in living organisms. Therefore, what would be expected is similarity between the results of these model systems. That is what makes the results presented in
Figure 10 even more interesting.
The complex was tested at concentrations between 3·10-6 M and 3·10-5 M. The ligand 2b was tested at three times greater molarities. In presence of both compounds, the behavior of this model system was extremely unstable. Experiments were repeated twice, on different days to lower the probability of human error or equipment failure as much as possible. Based on observations made, both compounds increase formazan production, i.e. increase RS formation, acting as prooxidants. At concentrations up to 3·10-5 M the ligand 2b manifests no statistically significant impact on absorbance at 578 nm, hence on RSA. At 1·10-4 M the effect is strong and statistically significant (RSA= 115±72%). At three times lower molarity, La(2b’)3 also seems to acts as a prooxidant, though to a smaller degree (RSA= 68±16%). At lower concentrations the effect statistically is zero. These results are unexpected for the authors – in the water radiolysis system both ligand and complex act as potent hydroxyl radical scavengers while at the same time MTT-formazan transformation, caused by Fenton-generated OH● seems to be potentiated by these same substances.
3.5.5. Impact of 2b and La(2b’)3 on LDCL in presence of KO2
The superoxide radical-ion is a type of RS that normally occurs in the human body, produced by one-electron reduction of oxygen. A by-product of oxygen metabolism, and various dedicated enzymes it also plays a role in defence against pathogens and in a variety of cell signalling pathways [
42,
43]. The ability of the ligand 2b and Ln(2b’)
3 to scavenge KO
2 derived superoxide is presented in
Figure 11.
At concentrations 3·10-3 M and below both compounds are inactive. At 1·10-4 M the ligand increases LDCL slightly (CL-SI= 108±7%). The effect becomes more evident at 3·10-4 M (CL-SI= 121±2.%) – evidence for slight prooxidant action. At 1·10-4 M La(2b’)3 has CL-SI= 118±6.%. The limited solubility of both compounds does not allow testing higher molarities. In terms of in vitro KO2 superoxide scavenging both ligand and its complex could be viewed as relatively inert.
3.5.6. Impact of 2b and La(2b’)3 on LDCL in presence of NaClO
Hypochlorous acid is naturally produced in the human body as an important component of immune defence [
44]. It is produced by neutrophils during the respiratory burst - the enzyme myeloperoxidase uses H
2O
2 and chloride ions as substrates to synthesize this highly toxic specie, a strong oxidant that damages a variety of molecular targets, acting as a bactericide. Its high chemical reactivity and non-specific action are the reason HClO is associated with variety of human pathologies [
45]. The ability of the ligand 2b and Ln(2b’)
3 to scavenge
in vitro NaClO derived superoxide is presented in
Figure 12.
When it comes to
in vitro interaction with hypochlorite, the divergent trends of behaviour are quite stark. At the lowest tested concentration, the ligand 2b increases LDCL (CL-SI= 139±12%) quite significantly, compared to the control. As molarity rises, CL-SI decreases in a concentration-dependent manner down to 64±40% at 1·10
-4 M. The opposite is observed with Ln(2b’)
3. At 1·10
-6 M it does not significantly impact LDCL (CL-SI=104±4%). At the next, higher concentration (3·10
-6 M) CL-SI increases and is statistically the same as that of 2b. That increase continues right up to 201±10% at 1·10
-4 M. A previous study [
46], using the same model system with another La(III) complex has demonstrated clear correlation between the behavior of both ligand and complex. In this case, the coordination of 2b with La(III) seems to show some kind of drastic change in behavior.