Research Article | Open Access
Palladium(II) Complexes Containing Mixed Nitrogen-Sulphur Donor Ligands: Interaction of [Pd(Methionine Methyl Ester)(H2O)2]2+ with Biorelevant Ligands
Pd(MME)Cl2 complex (MME = methionine methyl ester) was synthesised and characterized by physicochemical measurements. The reaction of [Pd(MME)(H2O)2]2+ with amino acids, peptides, or dicarboxylic acids was investigated at 25°C and 0.1 M ionic strength. Amino acids and dicarboxylic acids form 1 : 1 complexes. Peptides form both 1 : 1 complexes and the corresponding deprotonated amide species. The stability of the complexes formed was determined and the binding centres of the ligands were assigned. Effect of solvent on the stability constant of Pd(MME)-CBDCA complex, taken as a representative example, shows that the complex is more favoured in a medium of low dielectric constant. The concentration distribution diagrams of the complexes were evaluated.
cis-platin [cis-diamminedichloroplatinum(II)] is one of the most active antitumor agents in clinical use . However, it has a narrow spectrum of activity, and its clinical use is limited by undesirable side effects, including nephrotoxicity, ototoxicity, nausea, vomiting, and myelosuppression. Significant efforts were made to produce novel platinum-containing compounds in order to improve structural and physiological disadvantages, such as several side effects, drug resistance, and limitation in the field of application. Most of the developed and studied platinum complexes include nitrogen-containing monodentate, bidentate, or tridentate ligands [2–8]. The ultimate aim of the modifications of the parent drug is to make related analogues that produce a different spectrum of DNA lesions and so circumvent the problem of resistance to cisplatin [2, 9] Therefore, nonclassical platinum derivatives were also investigated, which may violate the classical structure-activity relationship . Complexes include sulphur-containing ligands and show antitumor activity . An example is [Pt(CH3SCH2CH2SCH3)Cl2], Pt(dt), a sulphur analogue of the well-studied [Pt(H2NCH2CH2NH2)Cl2], Pt(en)Cl2 complex . Baltic et al. examined the antitumor activity of Pt(dt) against the human breast cancer cell line and found the complex to inhibit the growth of MCF-7 cells in a dose and time dependent manner. Another example is the complex dichloro(2-methylthiomethylpyridine)platinum(II)Pt(mtp)Cl2 a complex with nitrogen as well as sulphur donor ligands. In the literature different synthetic pathways are described [11–13] and the bidentate N, S-complex was found to be promising cytostatic agent .
Pd(II) and Pt(II) amine complexes have the same structure, with five orders of magnitude higher reactivity in the case of Pd(II) complexes, but similar thermodynamic parameters. Pd(II) complexes are good models for the analogous Pt(II) complexes in solution. Recent work in our laboratories focused on the kinetics and equilibria of complex-formation reactions of cis-(diamine)palladium(II) complexes with DNA, the major target in chemotherapy of tumours, and biorelevant ligands as amino acids, peptides, dicarboxylic acids, and esters [15–18]. The present investigation describes the equilibria associated with the interaction of [Pd(MME)(H2O)2]2+ with biorelevant ligands as amino acids, peptides, and dicarboxylic acids. The results may give information regarding the distribution of the antitumor agent in biological fluid.
2.1. Materials and Reagents
All reagents were of analytical grade. K2PdCl4 and Methionine methyl ester hydrochloride were obtained from Aldrich. The amino acids and related compounds (glycine, alanine, -phenylalanine, proline, methionine imidazole, lysine, histidine, and histamine) were provided by Sigma Chemical Co. The peptides and amides used (glycinamide, glycylglycine, asparagine, and glutamine) and the dibasic acids used (cyclobutane 1,1-dicarboxylic acid, oxalic acid, malonic acid, succinic acid, and adipic acid) were all provided by BDH-Biochemicals Ltd., Poole, England. For equilibrium studies, [Pd(MME)Cl2] was converted into the diaqua complex by treating it with two equivalents of AgNO3 as described before . The ligands in the form of hydrochlorides were converted into the corresponding hydronitrates  All solutions were prepared in deionized water.
Pd(MME)Cl2 was prepared by dissolving K2PdCl4 (2 mmol) in 10 mL water with stirring. The clear solution of [PdCl4]2- was filtered and Methionine methyl ester hydrochloride (2 mmol), dissolved in 10 mL H2O, was added dropwise to the stirred solution. A yellowish-brown precipitate of [Pd(MME)Cl2] was formed and stirred for a further 30 minutes at 50°C. After filtering of the precipitate. A yellow powder was obtained. An orange crystalline precipitate was obtained; yield 98.9%. Anal. Calcd., for PdC6H13SNO2Cl2 (F. Wt = 340.38): C, 21.16; H, 3.81; N, 4.11; S, 9.40. Found: C, 21.2; H, 4.04; N, 3.7; S, 9.1.
Potentiometric titrations were performed with a Metrohm 686 titroprocessor equipped with a 665 Dosimat. The titroprocessor and electrode were calibrated with standard buffer solutions, prepared according to NBS specification . The pH meter readings were converted to hydrogen ion concentration by titrating a standard HNO3 solution (0.01 M), the ionic strength of which was adjusted to 0.1 M with NaNO3, with standard NaOH (0.05 M). The pH was plotted against p[H]. The relationship pH − p[H] = 0.05 was observed. All titrations were carried out at °C in purified nitrogen atmosphere using a titration vessel described previously . The in dioxane-water solution was determined as described previously [16, 17]. For this purpose, various amounts of standard NaOH solution were added to a solution containing 0.1 M NaNO3. The [OH−] was calculated from the amount of base added. The [H+] was calculated from the pH value. The product of [OH−] and [H+] was taken. IR spectra were measured on an 8001-PC FT-IR Shimadzu spectrophotometer using KBr pellets.
2.4. Procedure and Measuring Technique
The acid dissociation constants of the ligands were determined by titrating 1 mmole samples of each with standard NaOH solutions. Ligands were converted into their protonated form with standard HNO3 solutions. The acid dissociation constants of the coordinated water molecules in [Pd(MME)(H2O)2]2+ were determined by titrating 1 mmole of complex with standard 0.05 M NaOH solution. The formation constants of the complexes were determined by titrating solution mixtures of [Pd(MME)(H2O)2]2+ (1 mmole) and the ligand in the concentration ratio of 1 : 1 for amino acids, peptides, and dicarboxylic acids. The titrated solution mixtures each had a volume of 40 mL and the titrations were carried out at 25°C and 0.1 M ionic strength (adjusted with NaNO3), such a condition is similar to what exists in biological system. A standard 0.05 M NaOH solution was used as titrant.
The species formed were characterized by the general equilibrium (1) for which the formation constants are given by where M, L, and H stand for [Pd(MME)(H2O)2]2+ ion, ligand, and proton, respectively. The calculations were performed using the computer program  MINIQUAD-75. The stoichiometry and stability constants of the complexes formed were determined by trying various possible composition models for the systems studied. The model selected was that which gave the best statistical fit and was chemically consistent with the magnitudes of various residuals, as described elsewhere . The accepted model gave smallest standard deviation and sum of square of residuals. The reliability of the proposed model is tested by comparing the experimental titration data and the theoretically simulated data obtained by MINIQUAD-75 program.
Tables 1, 2, 3, and 4 list the stability constants together with their standard deviations and the sum of the squares of the residuals derived from the MINIQUAD output. The concentration distribution diagrams were obtained with the program SPECIES  under the experimental condition used.
aM, L, and H are the stoichiometric coefficients corresponding to Pd(MME), amino acid, and H+, respectively; the coefficient −1 refers to a proton loss; b of Pd(MME)-amino acids complexes, standard deviations are given in parentheses; sum of square of residuals is less than ; cthe p of the ligands and the protonated complexes. |
aM, L, and H are the stoichiometric coefficients corresponding to Pd(MME), peptides, and H+, respectively; the coefficient −1 refers to a proton loss; b of Pd(MME)-peptide complexes, standard deviations are given in parentheses and sum of square of residuals is less than ; cthe p of the peptides or of coordinated peptides.|
aM, L, and H are the stoichiometric coefficients corresponding to Pd(MME), dibasic acids, and H+, respectively. b of Pd(MME)-dibasic acid complexes; standard deviations are given in parentheses; sum of square of residuals is less than . cThe p of ligands or the protonated complexes.|
3. Results and Discussion
The analytical data indicated that the complex is of 1 : 1 stoichiometry and of formula Pd(MMA)Cl2. From comparing the IR spectrum of MME and Pd(MME)Cl2 complex, the stretching vibration of NH group of methionine methyl ester located at 2954 cm−1 is broadened in the complex spectrum due to coordination by the amino group. The ester carbonyl group of the ligand exhibits a sharp band at 1745 cm−1. This band was not affected upon complex formation. This indicates that the ester group is not participated in complex formation and methionine methyl ester is bound to Pd(II) by the amino and sulphur atom.
Before investigating the complex formation equilibria involving [Pd(MME)(H2O)2]2+ the hydrolysis of the ester group of methionine methyl ester was tested by pH-stat technique . The results show that the extent of ester hydrolysis in the pH range of the complex formation study is very small and can be neglected.
3.1. Acid-Base Equilibria of [Pd(MME)(H2O)2]2+
The acid-base equilibria of the [Pd(MME)(H2O)2]2+ complex is characterized by fitting the potentiometric data to various models. The best fit model was found to be consistent with species 10-1, 10-2, and 20-1. The first two species, 10-1 and 10-2, are due to deprotonation of the two coordinated water molecules as given by (3) and (4).(3)(4)The third species, 20-1, is the dimeric species formed as a result of the combination of the monoaqua hydroxo species (10-1) with the diaqua species (100) [Pd(MME)(H2O)2]2+ as given by:(5)
The formation of the dihydroxo-bridged dimer (20-2), found for most Pd-diimine complexes, is not favoured in the case of the Pd-MME complex. This may be accounted for on the basis that the strong labilization effect of the S-donor atom will cause the dimeric form (20-2) to be strained and consequently energetically not favoured .
The concentration distribution diagram for [Pd(MME)(H2O)2]2+ and its hydrolysed species is shown in Figure 1. The concentration of monohydroxo species, 10-1, and the dimeric species increase with increasing pH. The dimeric species, 20-1, attains maximum concentration of ~24% at pH ~5.2. The monohydroxo species (10-1) is the main species in the pH rang ~5.5 to 7.5; that is, it is the main species present in the physiological conditions. A further increase in pH is accompanied by an increase in the dihydroxo species (10-2) concentration, which is the main species above pH ca.8.
3.2. Complex Formation Equilibria Involving Amino Acids
Fitting of pH titration data for Pd(MME)-amino acid equilibria indicated the formation of 1 : 1 complexes with high stability constant values. This indicates that amino acids bind through the amino and carboxylate groups. Histidine is a tridentate ligand having amino, imidazole, and carboxylate groups as binding sites. With [Pd(MMA)(H2O)2]2+, only two of the three binding sites are involved in the complex formation; hence histidine coordinates in either a glycine-like or a histamine-like mode. The stability constant of the histidine complexes is significantly higher than that of histamine and amino acids. This indicates that histidine interacts with the complex by the amino and imidazole groups. The proline complex has the highest value among the simple amino acids. This may be due to the highest basicity of the proline amino group as reflected by the highest value. The stability constant of lysine complex is not so higher, as compared with the stability constants of the other amino acids. This may be explained on the premise that lysine coordinates by the two amino groups forming a less stable eight-membered chelate ring. The stability of methionine complex (logβ = 9.05) is lower than that of most simple amino acids. This may be explained on the premise that the amino group of methionine is less basic than that of other amino acids as reflected by values.
3.3. Complex Formation Equilibria Involving Peptides
The potentiometric data of the peptide (HL) complexes were fitted assuming the formation of the species [Pd(MME)(L)]+ (110 species) and [Pd(MME)(LH−1)] (11-1 species). The former species is formed by coordination through the amino group and carbonyl oxygen atom. On increasing the pH, the coordination site should switch from the carbonyl oxygen to the amide nitrogen with release of the amide hydrogen, forming the complex [Pd(MME)(LH−1)]. Such changes in coordination centers are now well documented [23, 24]. The of the coordinated amide group was calculated using (7) and is given in Table 2:
The for the glycinamide complex is lower than that of the other peptides. This may be explained on the premise that the more bulky substituent group on the peptide may serve to hinder the structural change in going from protonated to deprotonated complexes. Asparagine complex has the highest stability constant value, most probably due to the presence of α-amino group that can coordinate firstly as glycine does. The -amino group of asparagine is more basic than that of other peptides, which result in an increase of stability constant of its complex.
The concentration distribution diagrams of peptide complexes indicate that all peptides form the complex species (110) at low pH and thus prevent the hydrolysis of Pd(II) ion; that is, the hydrolysed species (10-1) and (10-2) are either not formed or formed in very low concentration. The induced ionization of the peptide hydrogen of glycinamide and glycylglycine starts at pH ~3. However asparagine ionization starts above pH ~8. Therefore, under normal physiological condition (pH 6-7), the peptides would coordinate to [Pd(MME)(H2O)2]2+ in entirely different ways. Glycinamide and glycylglycine are present entirely in the deprotonated form (11-1), whereas asparagine exists solely in the protonated form.
3.4. Complex Formation Equilibria Involving Dicarboxylic Acids
In the case of dicarboxylic acids, the potentiometric data were fitted on the basis of the formation of 110 complex and its monoprotonated form (111). The formation constants of the (110) complexes formed with oxalic, CBDCA and malonic acids, where five- and six-membered chelate rings are higher than those involving seven membered, as in succinic, and nine-membered chelate rings as in adipic acid. This may be explained on the premise that the five- and six-membered rings are more favored energetically than the seven- and the nine-membered rings. It is interesting to note that 1,1-CBDCA has a higher stability constant than that of malonic acid, although both form 6-membered chelate rings. This may be due to the higher values of the former than the latter dicarboxylic acid. The value of the protonated complex species of [Pd(MME)(CBDCA)] is 1.82. This value is lower than that of free CBDCAH-, which indicates acidification of the second carboxylic group upon coordination of [Pd(MME)(H2O)2]2+ with the first carboxylate group. The value of this protonated species in case of [Pd(en)HCBDCA]+ was estimated before from u.v./vis measurements to be ca.2.5 at 25°C and 0.1 M ionic strength . This species was documented to be the active form in the case of carboplatin . The concentration distribution diagram of CBDCA complex, taken as an example, is given in Figure 2 and shows that the protonated species (111) is stable only at low pH and the species (110) is predominating at the physiological pH. The hydrolysed species for all dibasic acid complexes predominates only at high pH.
3.5. Effect of Solvent on Complex Formation Equilibria
It was reported that a lower polarity has been detected in some biochemical microenvironments, such as active sites of enzymes and side chains in proteins [27–30]. These properties approximately correspond to those (or can be simulated by those) existing in the water/dioxane mixtures. Consequently, a study of the Pd(MME)-CBDCA complex formation, taken as a typical example, in dioxane-water solutions of different compositions, could be of biological significance. In order to characterize the formation equilibria of the Pd(MME)-CBDCA complex in dioxane-water solutions, all other equilibria involved, namely, acid-base equilibria of CBDCA and [Pd(MME)(H2O)2]2+, have to be studied in the same solvent. The equilibrium constants are reported in Table 4.
The hydrolysis of [Pd(MME)(H2O)2]2+ complex in dioxane-water solution leads to the formation of mono- and dihydroxy species. The dihydroxo bridged dimer was not detected. The values of CBDCA and those of the coordinated water molecules in [Pd(MME)(H2O)2]2+ increase linearly with increasing of dioxane concentration. This may be correlated with the ability of a solvent of relatively low dielectric constant to increase the electrostatic attraction between the proton and ligand anion in case of CBDCA and that between a proton and the hydrolyzed form of Pd(II) species. The variation in stability constant of the [Pd(MME)(H2O)2]2+ complex with CBDCA as a function of solvent composition of range from 12.5 to 62.5% is investigated. The results are displayed in Figures 3 and 4. The stability constant for the Pd(MME)-CBDCA complex increases with increasing dioxane concentration. This is explained in terms of complex formation involving oppositely charged ions as in the Pd(MME)-CBDCA complex, which is favored by the low dielectric constant of the medium, that is, with increasing dioxane concentration. The results show that the CBDCA complex with [Pd(MME)(H2O)2]2+ will be more favored in biological environments of lower dielectric constant. It is emphasized that the trend of stability constants as a function of solvent dielectric constant is particularly evident when charged reactants are involved in complex formation.
Pd(MMA)Cl2 complex was synthesised and characterized. In comparison of stability constants of [Pd(MME)(H2O)2]2+ complexes with biorelevant ligands, it would be possible to evaluate the speciation of Pd(II) complexes in biological fluid. This would form a clear basis for understanding the mode of action of such metal species under physiological condition.
Pd(MME)(CBDCA) is taken as a model for carboplatin drug. The study of its stability in dioxane-water solutions of different compositions could be of biological significance. The results show that the formation of the complex is more favoured in biological environments of lower dielectric constant.
Conflict of Interests
The authors declare that they have no conflict of interests regarding the publication of this paper.
- B. Lippert, Ed., Chemistry and Biochemistry of Leading Anticancer Drugs, Weinheim, Germany, Wiley-VCH, 1999.
- E. Wong and C. M. Giandornenico, “Current status of platinum-based antitumor drugs,” Chemical Reviews, vol. 99, no. 9, pp. 2451–2466, 1999.
- L. A. S. Costa, W. R. Rocha, W. B. de Almeida, and H. F. Dos Santos, “The solvent effect on the aquation processes of the cis-dichloro(ethylenediammine)platinum(II) using continuum solvation models,” Chemical Physics Letters, vol. 387, no. 1–3, pp. 182–187, 2004.
- N. Farrell, “DNA-binding and chemistry of dinuclear platinum complexes,” Comments on Inorganic Chemistry, vol. 16, no. 6, pp. 373–389, 1995.
- V. Brabec, J. Kašpárková, O. Vrána et al., “DNA modifications by a novel bifunctional trinuclear platinum phase I anticancer agent,” Biochemistry, vol. 38, no. 21, pp. 6781–6790, 1999.
- Y. Qu and N. Farrell, “Interaction of bis(platinum) complexes with the mononucleotide 5'-guanosine monophosphate. Effect of diamine linker and the nature of the bis(platinum) complex on product formation,” Journal of the American Chemical Society, vol. 113, no. 13, pp. 4851–4857, 1991.
- N. Farrell and Y. Qu, “Chemistry of bis(platinum) complexes. Formation of trans derivatives from tetraamine complexes,” Inorganic Chemistry, vol. 28, no. 18, pp. 3416–3420, 1989.
- N. Farrell, Y. Qu, L. Feng, and B. Van Houten, “Comparison of chemical reactivity, cytotoxicity, interstrand cross-linking and DNA sequence specificity of bis(platinum) complexes containing monodentate or bidentate coordination spheres with their monomeric analogues,” Biochemistry, vol. 29, no. 41, pp. 9522–9531, 1990.
- L. Kelland, “The resurgence of platinum-based cancer chemotherapy,” Nature Reviews Cancer, vol. 7, no. 8, pp. 573–584, 2007.
- G. Bogdanović, V. Kojić, T. Srdić et al., “Growth effects of some platinum(Ii) complexes with sulfur-containing carrier ligands on MCF7 human breast cancer cell line upon simultaneous administration with taxol,” Metal-Based Drugs, vol. 9, no. 1-2, pp. 33–43, 2002.
- A. Baldo, G. Chessa, G. Marangoni, and B. Pitteri, “The reactions between tetrachloroplatinate(II) anion and bidentate sulphur-nitrogen chelating agents: a direct route for the synthesis of the neutral dichloro species,” Polyhedron, vol. 4, no. 8, pp. 1428–1431, 1985.
- R. C. Jones, R. L. Madden, B. W. Skelton et al., “Solution, structural and catalytic studies of neutral MCl2 (M = Pd, Pt) complexes of the N/E mixed-donor ligands 2-(RECH2)C5H4N (RE = MeS, PhS, MeSe),” European Journal of Inorganic Chemistry, vol. 2005, no. 6, pp. 1048–1055, 2005.
- R. C. Jones, B. W. Skelton, V. A. Tolhurst, A. H. White, A. J. Wilson, and A. J. Canty, “Synthesis and solid-state structural characterization of Pt(II, IV) bromide complexes containing palladium(II) and biologically relevant nucleophiles,” European Journal of Inorganic Chemistry, pp. 2261–2270, 2009.
- G. H. Du Preez, “Platinum(II) complexes, preparation and use, document number,” CA Patent 2547275, 2005.
- T. Soldatović, M. Shoukry, R. Puchta, Ž. D. Bugarčić, and R. Ven Eldik, “Equilibrium and kinetic studies of the reactions between aqua[1-(2-aminoethyl)piperazinejpalladium(II) and biologically relevant nucleophiles,” European Journal of Inorganic Chemistry, no. 15, pp. 2261–2270, 2009.
- M. R. Shehata, M. M. Shoukry, F. H. Abdel-Shakour, and R. V. Eldik, “Equilibrium studies on complex-formation reactions of pd[(2-(2-aminoethyl) pyridine)(H2O)2]2+ with ligands of biological significance and displacement reactions of DNA constituents,” European Journal of Inorganic Chemistry, no. 26, pp. 3912–3920, 2009.
- M. R. Shehata, M. M. Shoukry, F. M. H. Nasr, and R. van Eldik, “Complex-formation reactions of dicholoro (S-methyl-L-cysteine) palladium (II) with bio-relevant ligands. Labilization induced by S-donor chelates,” Dalton Transactions, no. 6, pp. 779–786, 2008.
- A. A. El-Sherif, M. M. Shoukry, and R. Van Eldik, “Complex-formation reactions and stability constants for mixed-ligand complexes of diaqua(2-picolylamine)palladium(II) with some bio-relevant ligands,” Dalton Transactions, no. 7, pp. 1425–1432, 2003.
- R. G. Bates, Determination of pH: Theory and Practice, Wiley Interscience, New York, NY, USA, 2nd edition, 1975.
- Ž. D. Bugarčič, M. M. Shoukry, and R. van Eldik, “Equilibrium and kinetic data for the interaction of diaqua-(S-methyl-L-cysteine)palladium(II) with biologically relevant ligands,” Journal of the Chemical Society, Dalton Transactions, no. 21, pp. 3945–3951, 2002.
- P. Gans, A. Sabatini, and A. Vacca, “An improved computer program for the computation of formation constants from potentiometric data,” Inorganica Chimica Acta, vol. 18, no. C, pp. 237–239, 1976.
- L. Pettit, Personal Communication, University of Leeds, Leeds, UK, 1993.
- M. Lim, “Mixed-ligand complexes of palladium. 5. Diaqua(ethylenediamine)palladium(II) complexes of ethanolamine, L-serine, L-threonine, L-homoserine, and L-hydroxyproline,” Inorganic Chemistry, vol. 20, no. 5, pp. 1377–1379, 1981.
- M. C. Lim, “Mixed-ligand complexes of palladium(II). Part 1. Diaqua(ethylenediamine) palladium(II) complexes of glycylglycine and glycinamide,” Journal of the Chemical Society, Dalton Transactions, no. 1, pp. 15–17, 1977.
- A. Shoukry, T. Rau, M. Shoukry, and R. van Eldik, “Kinetics and mechanisms of the ligand substitution reactions of bis(amine)(cyclobutane-1,1-dicarboxylato)palladium(II),” Journal of the Chemical Society—Dalton Transactions, no. 18, pp. 3105–3112, 1998.
- U. Frey, J. D. Ranford, and P. J. Sadler, “Ring-opening reactions of the anticancer drug carboplatin: NMR characterization of cis-[Pt(NH3)2(CBDCA-O)(5′-GMP-N7)] in solution,” Inorganic Chemistry, vol. 32, no. 8, pp. 1333–1340, 1993.
- U. Kramer-Schnabel and P. W. Linder, “Substituent effects in the protonation and complexation with copper(II) ions of organic monophosphate esters. A potentiometric and calorimetric study,” Inorganic Chemistry, vol. 30, no. 6, pp. 1248–1254, 1991.
- D. C. Rees, “Experimental evaluation of the effective dielectric constant of proteins.,” Journal of Molecular Biology, vol. 141, no. 3, pp. 323–326, 1980.
- N. K. Rogers, G. R. Mooreb, and M. J. E. Strenberg, “Electronic interactions in globular proteins: calculation of the pH dependence of the redox potential of cytochrome,” Journal of Molecular Biology, vol. 182, pp. 613–616, 1985.
- G. Akerlof and O. A. Short, “Chemistry of the aliphatic esters of thiophosphoric acids,” Journal of the American Chemical Society, vol. 75, no. 24, p. 6357, 1953.
Copyright © 2014 Mohamed M. Shoukry and Sameya M. T. Ezzat. 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.