1. Introduction
Photodynamic therapy involves the marking of diseased tissue with a photosensitizer followed by exposure to light (the wavelength is generally in visible spectrum) to induce a phototoxic reaction which leads to its destruction. The main field of application of photodynamic therapy is oncology. Photosensitizers preferentially accumulate in neoplastic tissues, allowing a priori selective destruction of tumors [
1,
2].
Phthalocyanines are promising photosensitizers for photodynamic therapy (PDT) due to their strong absorption of red light that penetrates tissues, easy derivatization, and high singlet oxygen generation efficiency [
3]. Three classes of metallophthalocyanines have been used in clinical trials, namely, silicon (IV) phthalocyanine [
4,
5], zinc (II) phthalocyanine [
6], and aluminum (III) phthalocyanine [
7]. However, the main problems with phthalocyanine derivatives are their low solubility in biological environments and poor retention in malignant tissues [
8]. The study of their photosensitization properties in a heterogeneous environment and the evaluation of their interactions with biological molecules and cellular components have thus received increased interest [
8]. The application of computational chemistry methods can provide necessary information on the stability and reactivity of these compounds.
Zhang et al. [
9] used DFT methods to investigate the structure and spectroscopic properties of phthalocyanines substituted by methoxy and sulfoxyl groups. In addition, theoretical chemistry methods have been used to study zinc phthalocyanines with grafted naphthoxy and phenylazonaftoxy moieties [
10] as well as phthalocyanines substituted with four imidazolyl rings [
11]. The efficiency of certain symmetrical (substituted with four tert-butyl residues) and asymmetrical (substituted with three tert-butyl residues and an iodine group) phthalocyanines as photosensitizers was investigated in [
12]. Symmetrically substituted phthalocyanines have been shown to be effective in PDT applied to the lung, intestine, ovarian, and pancreatic tumors, while asymmetrically substituted phthalocyanines are promising imaging agents for colon and pancreatic tumors [
12]. In [
13], the authors outlined the high biological activity of fluorinated phthalocyanines, with their results showing complete regression of tumors in 20% of treated mice.
Another area of interest is represented by the study of delivery systems based on colloidal structures due to their potential for transporting different drugs while protecting them against degradation. Drug delivery systems with adjustable properties have the ability to host a variety of molecules. For example, vesicularly cross-linked hydrogels are able to provide different environments for pharmaceuticals comprising bulk water, sites on polymeric chains, the interior of vesicles, and bilayers [
14]. Most phthalocyanines are hydrophobic compounds that require association with other drug delivery systems. A number of phthalocyanines derivatized with chlorine and aluminum have shown promising results in PDT, producing high concentrations of reactive oxygen species. These results have been confirmed by a series of in vitro tests on A549 cancer lines, showing high phototoxicity [
15]. Pyridine derivatives are a type of important pharmaceutical intermediates which have excellent biocompatibility; they can form hydrogen bonds with biological molecules thanks to the presence of a nitrogen atom. Thus, [
16] investigated the in vitro photodynamic activity of Zn (II) phthalocyanines substituted with pyridine groups. Other results have shown that a tetrasubstituted zinc phthalocyanine has a low tendency to aggregate in the biological environment and a high level of reactive oxygen species, giving it high photodynamic activity [
17].
The present study deals with the investigation of variously substituted metal-phthalocyanines (with amino and methyl groups) as possible photosensitizers, the central metal being Fe, Ni, and Zn, respectively. In the following, amino-substituted phthalocyanines are denoted by I, methyl derivatives are denoted by II, and unsubstituted phthalocyanines by III. According to data from the literature, a compound with applications in photodynamic therapy is characterized by the following properties [
18]: low HOMO orbital energy (which is an advantage for electronic transfer); a small HOMO–LUMO gap; low dipole moment; and UV–Vis absorption in the range of 400–650 nm.
The general structure of the investigated metal-phthalocyanines is depicted in the
Figure 1.
The second part of this study deals with the evaluation of four fullerenes as nanocarriers for these phthalocyanines. We employed an all-carbon fullerene (C
52) and three fullerenes doped with nitrogen and phosphorus atoms (C
46N
3P
3) with different substitution patterns. The heteroatomic substitution pattern is depicted below for both sumanene, a circulene-type precursor of fullerene C
52, and the corresponding fullerenes. Our previous studies have dealt with the evaluation of the properties and aromatic characteristics of various all-carbon and nitrogen- and phosphorus-doped nanostructures, including both fullerenes and their polycyclic aromatic hydrocarbon (mostly coronene and sumanene) precursors [
19,
20].
The heteroatom substitution pattern is illustrated in
Figure 2.
The present paper aims to (i) evaluate the properties of nine metal-phthalocyanines in order to establish whether they are suitable photosensitizers and (ii) investigate the interactions between the aforementioned compounds and several heteroatom-doped fullerenes in order to use the latter compounds as delivery systems. Both the influence of the central metal and the role of the substituents on the properties of the title compounds are evaluated.
2. Materials and Methods
The phthalocyanines were first optimized at the semi-empirical level (PM6 method implemented in Gaussian 09W), and the obtained structures were refined at the HF/LanL2DZ level of theory. A similar approach was employed for the four fullerenes. These computations were performed by means of G09W software [
21]. In addition, the energies of the HOMO and LUMO orbitals and the dipole moment were computed by means of G09 W at the HF/LanL2DZ level of theory.
Chemical potential (µ), chemical hardness (η), and electrophilicity (ω) were calculated using the following equations:
The steric parameters [
22], namely, ovality, Connolly accessible area, and Connolly solvent excluded volume, as well as the partition coefficient logP, were obtained with Chem3D software. Chem3D was employed for the computations of the number of hydrogen bond donors and acceptors as well. AutoDock Vina software [
23] was employed for the molecular docking studies.
The fullerenes were considered as receptors and a grid box of 40 × 40 × 40 Å was used, with the center of the grid box considered the center of the fullerenes. The optimized structures of the investigated phthalocyanines were the ligands, and the torsions were assigned along the rotatable bonds. The binding constant K
B was calculated using the following equation:
where ΔG
B is the binding affinity (J·mol
−1), R is the gas constant (J·mol
−1·K
−1), and T is the temperature (298 K).
4. Discussion
The present paper has dealt with the investigation of the properties of iron-, nickel- and zinc-phthalocyanines as possible photosensitizers. In this regard, we computed several properties, including the HOMO–LUMO gap value, dipole moment, and wavelength of maximum absorbance, with the results suggesting that all of the investigated compounds possess the requirements of good photosensitizers.
Concerning the influence of the central metal on the properties of the studied compounds, the following conclusions can be stated:
- (i).
The energies of the frontier molecular orbitals HOMO and LUMO are very slightly influenced by the type of the metal; as a result, the calculated global reactivity parameters presented in
Table 1,
Table 2,
Table 3 and
Table 4 have similar values among the same substitution type, regardless of whether the metal is iron, nickel, or zinc;
- (ii).
The central metal influences the steric parameters of the phthalocyanines, with a larger nitrogen-iron bond (1.97 Å, compared to 1.88 Å for nickel- and zinc-nitrogen) leading to different values of the Connolly Accessible Area and Connolly Solvent-Excluded Volume for iron-phthalocyanines on the one hand and nickel- and zinc-phthalocyanines on the other.
The main influence on the properties of the metal-phthalocyanines arises from the presence of the substituents. The differences are outlined in the calculated reactivity descriptors depicted in
Table 1,
Table 2,
Table 3 and
Table 4. It can be observed that the presence of the methyl groups leads to higher energies for the HOMO and LUMO orbitals, while the amino-substituted compounds are characterized by lower energies as compared to unsubstituted metal-phthalocyanines.
Another part of the present study dealt with the interactions between each metal-phthalocyanine and four fullerenes consisting of 52 atoms (one all-carbon and three fullerenes doped with isovalent nitrogen and phosphorus atoms) as possible delivery systems. The best affinities were obtained for the phthalocyanines substituted with methyl groups and where the central metal was Ni and Zn.
Based on the calculated properties, the characteristics of the ideal metal-phthalocyanines within the investigated series are those shown in
Table 19.