# Porphyrins and Metalloporphyrins Combined with N-Heterocyclic Carbene (NHC) Gold(I) Complexes for Photodynamic Therapy Application: What Is the Weight of the Heavy Atom Effect?

Stefano Scoditti, Francesco Chiodo, Gloria Mazzone, Sébastien Richeter, Emilia Sicilia

► To cite this version:

Stefano Scoditti, Francesco Chiodo, Gloria Mazzone, Sébastien Richeter, Emilia Sicilia. Porphyrins and Metalloporphyrins Combined with N-Heterocyclic Carbene (NHC) Gold(I) Complexes for Photodynamic Therapy Application: What Is the Weight of the Heavy Atom Effect?. Molecules, 2022, 27 (13), pp.4046. 10.3390/molecules27134046. [hal-03735800](https://hal.umontpellier.fr/hal-03735800)

HAL Id: hal-03735800

https://hal.umontpellier.fr/hal-03735800

Submitted on 18 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Porphyrins and Metalloporphyrins Combined with N-Heterocyclic Carbene (NHC) Gold(I) Complexes for Photodynamic Therapy Application: What Is the Weight of the Heavy Atom Effect?

Stefano Scoditti <fnref n="1" /> ID, Francesco Chiodo <fnref n="1" />, Gloria Mazzone <fnref n="1" />,* ID, Sébastien Richeter <fnref n="2" /> and Emilia Sicilia <fnref n="1" />,* ID

**Abstract:** The photophysical properties of two classes of porphyrins and metalloporphyrins linked to N-heterocyclic carbene (NHC) Au(I) complexes have been investigated by means of density functional theory and its time-dependent extension for their potential application in photodynamic therapy. For this purpose, the absorption spectra, the singlet-triplet energy gaps, and the spin-orbit coupling (SOC) constants have been determined. The obtained results show that all the studied compounds possess the appropriate properties to generate cytotoxic singlet molecular oxygen, and consequently, they can be employed as photosensitizers in photodynamic therapy. Nevertheless, on the basis of the computed SOCs and the analysis of the metal contribution to the involved molecular orbitals, a different influence in terms of the heavy atom effect in promoting the intersystem crossing process has been found as a function of the identity of the metal center and its position in the center of the porphyrin core or linked to the peripheral NHC.

**Keywords:** PDT; porphyrins; metalloporphyrins; gold complexes; carbenes; DFT; TDDFT

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Scoditti, S.; Chiodo, F.; Mazzone, G.; Richeter, S.; Sicilia, E. Porphyrins and Metalloporphyrins Combined with N-Heterocyclic Carbene (NHC) Gold(I) Complexes for Photodynamic Therapy
Application: What Is the Weight of the Heavy Atom Effect? *Molecules* **2022**, *27*, 4046. https://doi.org/10.3390/molecules27134046

Academic Editors: Tiziano Marzo and Iogann Tolbatov

Received: 30 May 2022

Accepted: 21 June 2022

Published: 23 June 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

Photodynamic therapy (PDT) is a clinical therapeutic modality with considerable potential for application in a wide range of diseases, such as cancer [1] as well as the antimicrobial and environmental fields [2–7]. A PDT protocol involves the combination of three key components: light, with an appropriate wavelength, molecular oxygen, and a photosensitizer, a molecule which can produce reactive oxygen species (ROS) that are deadly to cells when excited by light of a proper wavelength in presence of oxygen. The production of ROS may occur through two pathways, named type I and type II processes. The type I mechanism involves an electron-transfer-mediated biomolecule oxidation that generates reactive oxygen species such as $O_2^{•−}$ and $HO_2^{•}$. The type II mechanism entails an energy transfer from the photosensitizer in its triplet state to molecular oxygen, promoting the formation of singlet oxygen, $^1O_2$. An efficient photosensitizer must possess a series of required chemical and photophysical properties, and some of them are of paramount importance: (i) a maximum absorption wavelength falling in the photodynamic therapeutic window (600–850 nm) to ensure a good penetration of light into living tissues [8,9]; (ii) a high intersystem spin crossing (ISC) probability between the excited singlet and triplet states; and (iii) a populated triplet state with an energy higher than 0.98 eV, the energy necessary to induce singlet oxygen formation from the ground triplet state [10–13].

Since the discovery of the potential of light-based therapies, numerous attempts have been made to treat tumors with photosensitizing agents. Porphyrins and their derivatives, such as chlorins or bacteriochlorins, surely represent one of the best classes of molecules for

<footnote n="1">¹ Department of Chemistry and Chemical Technologies, Università della Calabria, 87036 Arcavacata di Rende, Italy; stefano.scoditti@unical.it (S.S.); francescochiodo22@gmail.com (F.C.)</footnote>

<footnote n="2">² ICGM, Univ. Montpellier, CNRS, ENSCM, 34293 Montpellier, France; sebastien.richeter@umontpellier.fr</footnote>

<footnote>* Correspondence: gloria.mazzone@unical.it (G.M.); emilia.sicilia@unical.it (E.S.)</footnote>
PDT application due to their inherent chemical stability and suitable photophysical proper-
ties [14–19]. They show characteristic UV-vis spectra composed of one strong absorption
band, named the Soret or B band, in the ~380–420 nm region and weaker absorption bands,
known as Q bands [20,21]. The number and the intensity of these Q bands strongly depend
on the structure of the porphyrin core and the nature of the metal ion in the macrocycle.
Free-base porphyrins generally possess four Q bands with the lowest energy one centered
at ~650 nm. Although the molar absorption coefficients of Q bands are lower compared to
the Soret band, the excitation of the porphyrin derivatives in these Q-bands is extensively
used for PDT since they fall within the therapeutic window, enabling a deep penetration of
red light into living tissues. Despite the excellent properties, such as their low dark toxicity,
their thermodynamic stability, and their maximum absorption bands in the red region of the
visible spectrum, the application of porphyrin derivatives is limited by shortcomings such
as low water solubility and poor accumulation into cancer cells [22–26]. Over the years,
several structural changes have been proposed to improve both their PDT efficiency and
water solubility, which include the modification of the π electron architecture, the addition
of charged groups, and the inclusion of heavy atoms [27]. It is well-known, indeed, that
the presence of heavy atoms in the molecular structure of the chromophore significantly
increases the spin–orbit coupling (SOC), a relativistic effect that causes a quantum mechan-
ical mixing between states with different spin multiplicities. This phenomenon is known as
the “Heavy Atom Effect” (HAE). Therefore, as established by the Fermi golden rule [28,29],
high SOC values improve the kinetics of the transitions between these states, favoring the
population of a triplet state required for PDT action. In order to increase the quantum
yield of the triplet state, two different strategies have been proposed for porphyrinoid
dyes: (i) including a metal center within the cavity of the tetrapyrrole macrocycle and
(ii) linking a peripheral metal complex to the π system of the chromophore. Many studies,
both experimental and computational, have been carried out dealing with such an effect
for numerous metalloporphyrins [13,18,30–32]. For PDT applications, the peripheral metal
complex might increase the SOC values and consequently boost the kinetics of the ISC,
enhancing the generation of $ ^1O_2 $ through the HAE. For example, Spingler, Gasser, and
coworkers also reported that $ ^1O_2 $ quantum yields of *meso* tetra(pyridyl)porphyrins are
increased upon the coordination of Pt(II) complexes, suggesting that the Pt(II) ions promote
the formation of the T1 state through the HAE [33]. The peripheral metal complexes can
also bring additional properties such as intrinsic anticancerous properties. Thus, Brunner
and coworkers reported the synthesis of porphyrin derivatives with peripheral Pt(II) com-
plexes for dual chemo- and phototherapy because they combine within the same molecular
species the antitumor activity of Pt(II) complexes and the phototoxicity of porphyrins upon
irradiation with light [34–38].

N-heterocyclic carbenes (NHCs) are important ligands in the field of organometallic chemistry, and several molecular systems combining porphyrins and NHC-metal complexes were investigated during the last decade for different applications, including catalysis and PDT [39]. In this framework, we reported the synthesis of porphyrins with peripheral NHC-Au(I) complexes and their use as photosensitizers for PDT [39–41]. Photophysical measurements showed that peripheral Au(I) complexes improve PDT efficiency [40]. The obtained Au(I) also presents moderate cytotoxicity in the dark [39,40]. Interestingly, we and others also showed that the anticancerous properties of NHC-Au(I) complexes can be significantly improved by changing the second ancillary ligand of the linear Au(I) complex [42,43]. To have a better comprehension of the photophysical properties of porphyrins equipped with peripheral NHC-Au(I) complexes and their photodynamic properties, we report the computational evaluation of the photophysical properties of two sets of porphyrins, starting from the corresponding free-base porphyrins (Fb) named **Fused-Fb** and **Meso-Fb**. In the former case, the porphyrin and the NHCs are fused together and two neighboring β-pyrrolic carbon atoms of the porphyrin core correspond to the C4 and C5 atoms of the imidazole-2-ylidene. In the second case, the porphyrin is a wingtip group of the NHC, and both are linked together through the formation of a single $C_{meso}-N_{NHC}$
bond. The structures of these ligands are depicted in Scheme 1 together with the adopted nomenclature for the **Fused** and **Meso** species in panels a and b, respectively. The corresponding Au(I) complexes have also been studied to establish how much the presence of the Au(I) ion can enhance the ISC conversion. Additionally, we propose for both systems and their Au(I) complexes the corresponding metalloporphyrins with Zn(II) and Pd(II) as inner metal ions (see Scheme 1) based on the hypothesis that the simultaneous presence of inner and outer metal ions such as Zn(II)/Au(I) or Pd(II)/Au(I) couples could synergistically increase the SOC values and the rate of the ISC process. Density functional theory and its time-dependent formulation, TDDFT, have been employed in order to compute the photophysical properties of the two series of porphyrins and compare them with the available experimental counterparts.

<figure><img src="image_4.png" /><figcaption>**Scheme 1.** Structure of the investigated compounds.</figcaption></figure>

## 2. Results and Discussion

Optimized geometrical structures of the two compounds named **Fused-FbAu** and **Meso-FbAu** are reported in Figure 1 along with some key geometrical parameters. To allow an easier visualization and highlight the distortion of the macrocycle caused by the fused NHC ligand of the **Fused** species, only the porphyrin cores and the NHC-Au(I) complexes are reported. The addition in the *meso* position of the NHC moiety, which is perpendicular to the porphyrin plane, does not provoke any geometrical change worth mentioning. The optimized structure of the free porphyrin fused with the NHC group has been superimposed with the experimental counterpart, as reported in Figure S1 of the Supplementary Materials, which shows a very good agreement. Notwithstanding the observed deviation from the planarity of the **Fused** porphyrin compounds, the *π* conjugation in the central core is not disturbed, and the main characteristics of the spectra of the unperturbed porphyrins are conserved, as will be illustrated in the next paragraphs.

The rationalization of the UV-vis absorption spectra of porphyrin-like compounds can be conducted in terms of the classical four-orbital model first applied by Gouterman [44]. On the basis of such model, indeed, both the strong Soret band (B) in the near-UV region and the low-intensity Q-bands in the visible or near-IR region are mainly generated by the mixed electronic transitions between the two highest occupied molecular orbitals (HOMO and HOMO-1) and the two lowest unoccupied ones (LUMO and LUMO+1), here named H, H-1, L, and L+1, respectively. It has been underscored above that one of the reasons for the interest in porphyrin-like compounds in PDT resides in the Q-band absorption maxima that can be red-shifted as much as possible to fall into the so-called therapeutic window where the highest tissue permeability occurs. This is the reason that although all the most important spectral features of the whole spectrum for all the investigated compounds have been calculated, the attention has been mainly focused on the properties of the Q bands.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Optimized geometrical structure of the central core of the **Fused-FbAu** and **Meso-FbAu** compounds. Both top and side views are reported together with some key geometrical parameters. Bond lengths are in Å, and valence and dihedral angles are in degrees.</figcaption></figure>

In order to select the most appropriate protocol for the analysis of the properties of the compounds under investigation, a preliminary benchmark of the maximum absorption wavelengths ($\lambda_{max}$) of the two diagnostic porphyrin bands was carried out for both the **Fused-Fb** and **Fused-FbAu** compounds, as reported in Table S1 of the Supplementary Materials, with respect to the UV-vis experimental data [40,45]. Both compounds taken into account show an experimental UV-vis absorption spectrum characterized by a strong B band and a weaker Q band. Looking at the maximum wavelength of both bands, the best agreement with the experimental data was obtained using the M06 functional [46] that, however, gives a Q-band wavelength blue shifted with respect to the experimental value. The experimental spectra of the complexes named **Meso-Fb**, **Meso-Zn**, and **Meso-ZnAu** [41] were used to check if the M06 functional reproduces the **Meso** series well. The comparison between M06 and the experimental data can be found in Table S2 of the Supplementary Materials, which shows how the $\lambda_{max}$ of both the B and Q bands is properly computationally reproduced. Then, the M06 xc functional was adopted to perform water solvent calculations for all the considered compounds. This choice, moreover, was in agreement with previous calculations carried out by some of us on porphyrin-like systems [13,47–49].

The main features of the computed two components of the Q band, generally denoted as Qx and Qy, are collected in Tables 1 and 2 for all the examined **Fused** and **Meso** compounds, respectively. The plots of the corresponding Gouterman's molecular orbitals are reported in Figure 2 for the species **Fused-Fb**, **Fused-FbAu**, **Fused-ZnAu**, and **Fused-PdAu**, whereas for the rest of compounds, their graphical representations are shown in Figure S2.

The analysis of the data reported in Tables 1 and 2 shows that for the free-base **Fused-Fb** the low-intensity Qy component is H→L in nature (about 70%), while the Qx one involves mainly the H→L+1 (more than 60%) transition. For the free-base **Meso-Fb**, the same transitions give the main contribution, but their weight is lower. The coordination of the porphyrin to the Au(I) center through the formation of a C<sub>NHC</sub>-Au bond does not introduce any significant change to the position or to the intensity of the bands. When the Zn(II) and Pd(II) cations were included in the porphyrin cavity, only a small blue-shift of about few nanometers was observed, which was slightly more significant when Pd(II) was
used as the inner metal ion. Finally, the contemporary presence of the peripheral Au(I) and the metal cation at the center of the cavity causes a few nanometers of displacement of the two Q-band components towards the blue region of the spectrum that, once again, is more substantial for the Pd complexes. Free-base porphyrin Q bands are, obviously, of a $\pi \to \pi^*$ nature, with the MOs localized on the central porphyrin ring, as shown in Figure 2 and Figure S2, and it is evident that the gold atom as well as the carbene moiety are largely isolated from the aromatic system of the porphyrin core. The frontier MOs involved in the most important transitions are not influenced by the presence of both the peripheral Au(I) and central Zn(II) ions such as the optical properties of such compounds. An observable contribution of the metal center to the L and L+1 orbitals was noticed, instead, for the Pd(II) complexes.

<table><caption>**Table 1.** Main absorption wavelengths ($\lambda$), excitation energies ($\Delta E$), oscillator strengths ($f$), and main configuration for the examined Fused compounds in water solvent, computed at M06/6-31G** level of theory.</caption><thead><tr><th>Compound</th><th>Excited State</th><th>λ (nm)</th><th>ΔE (eV)</th><th>f</th><th>Main Configuration (%)</th></tr></thead><tbody><tr><td rowspan="2">Fused-Fb</td><td>S<sub>1</sub></td><td>605</td><td>2.05</td><td>0.084</td><td>H→L (70%), H-1→L+1 (29%)</td></tr><tr><td>S<sub>2</sub></td><td>570</td><td>2.17</td><td>0.047</td><td>H→L+1 (61%), H-1→L (38%)</td></tr><tr><td rowspan="2">Fused-Zn</td><td>S<sub>1</sub></td><td>568</td><td>2.18</td><td>0.047</td><td>H→L (62%), H-1→L+1 (37%)</td></tr><tr><td>S<sub>2</sub></td><td>567</td><td>2.18</td><td>0.022</td><td>H→L+1 (59%), H-1→L (41%)</td></tr><tr><td rowspan="2">Fused-Pd</td><td>S<sub>1</sub></td><td>543</td><td>2.28</td><td>0.019</td><td>H→L+1 (59%), H-1→L (41%)</td></tr><tr><td>S<sub>2</sub></td><td>541</td><td>2.29</td><td>0.033</td><td>H→L+1 (59%), H-1→L+1 (41%)</td></tr><tr><td rowspan="2">Fused-FbAu</td><td>S<sub>1</sub></td><td>604</td><td>2.05</td><td>0.063</td><td>H→L (67%), H-1→L+1 (32%)</td></tr><tr><td>S<sub>2</sub></td><td>571</td><td>2.17</td><td>0.055</td><td>H→L+1 (62%), H-1→LUMO (38%)</td></tr><tr><td rowspan="2">Fused-ZnAu</td><td>S<sub>1</sub></td><td>570</td><td>2.17</td><td>0.034</td><td>H→L (60%), H-1→L+1 (39%)</td></tr><tr><td>S<sub>2</sub></td><td>568</td><td>2.18</td><td>0.038</td><td>H→L+1 (60%), H-1→L (39%)</td></tr><tr><td rowspan="2">Fused-PdAu</td><td>S<sub>1</sub></td><td>544</td><td>2.28</td><td>0.028</td><td>H→L (60%), H-1→L+1 (39%)</td></tr><tr><td>S<sub>2</sub></td><td>541</td><td>2.29</td><td>0.024</td><td>H→L+1 (56%), H-1→LUMO (43%)</td></tr></tbody></table>

<table><caption>**Table 2.** Main absorption wavelengths ($\lambda$), vertical singlet electronic energies ($\Delta E$), oscillator strengths ($f$), and main configuration for the examined Meso compounds in water solvent, computed at M06/6-31G** level of theory.</caption><thead><tr><th>Compound</th><th>Excited State</th><th>λ (nm)</th><th>ΔE (eV)</th><th>f</th><th>Main Configuration</th></tr></thead><tbody><tr><td rowspan="2">Meso-Fb</td><td>S<sub>1</sub></td><td>585</td><td>2.12</td><td>0.024</td><td>H→L (45%), H-1→L+1 (24%)</td></tr><tr><td>S<sub>2</sub></td><td>553</td><td>2.24</td><td>0.031</td><td>H→L+1 (38%), H-1→L (31%),</td></tr><tr><td rowspan="2">Meso-Zn</td><td>S<sub>1</sub></td><td>557</td><td>2.23</td><td>0.027</td><td>H→L (62%), H-1→L+1 (38%)</td></tr><tr><td>S<sub>2</sub></td><td>551</td><td>2.25</td><td>0.001</td><td>H-1→L (50%), H→L+1 (50%)</td></tr><tr><td rowspan="2">Meso-Pd</td><td>S<sub>1</sub></td><td>533</td><td>2.33</td><td>0.024</td><td>H→L (61%), H-1→L+1 (38%)</td></tr><tr><td>S<sub>2</sub></td><td>528</td><td>2.35</td><td>0.000</td><td>H-1→L (53%), H→L+1 (46%)</td></tr><tr><td rowspan="2">Meso-FbAu</td><td>S<sub>1</sub></td><td>588</td><td>2.11</td><td>0.024</td><td>H→L (44%), H-1→L+1 (23%), H→L+1 (20%)</td></tr><tr><td>S<sub>2</sub></td><td>555</td><td>2.23</td><td>0.037</td><td>H→L+1 (39%), H-1→L (28%), H→L (21%)</td></tr><tr><td rowspan="2">Meso-ZnAu</td><td>S<sub>1</sub></td><td>560</td><td>2.22</td><td>0.025</td><td>H→L (62%), H-1→L+1 (38%)</td></tr><tr><td>S<sub>2</sub></td><td>554</td><td>2.24</td><td>0.008</td><td>H→L+1 (53%), H-1→LUMO (46%)</td></tr><tr><td rowspan="2">Meso-PdAu</td><td>S<sub>1</sub></td><td>533</td><td>2.32</td><td>0.020</td><td>H→L (60%), H-1→L+1 (39%)</td></tr><tr><td>S<sub>2</sub></td><td>529</td><td>2.34</td><td>0.004</td><td>H→L+1 (50%), H-1→L (49%)</td></tr></tbody></table>

The vertical excitation energies for the triplet states, together with the most significant molecular orbital contributions, are listed in Table S3. The triplet states lying below the first singlet excited state are four. From the computed vertical transition energies for all the examined compounds, it clearly appears that the energy gap between the ground singlet and the lowest-lying triplet excited state ($\Delta E_{S0-T1}$) is higher than the energy required to
promote molecular oxygen in its ¹Δg state (0.98 eV), and consequently, they are potential candidates to act as photosensitizers by the type II mechanism. The generation of cytotoxic singlet oxygen by the type II mechanism relies on the possibility of the system to trigger the radiationless ISC process between a singlet excited state (S<sub>n</sub>) and a triplet state lying below it (T<sub>m</sub>).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Four Gouterman's orbital contour plots for **Fused-Fb**, **Fused-FbAu**, **Fused-ZnAu**, and **Fused-PdAu**.</figcaption></figure>

As underlined above, the probability to populate the triplet states through an ISC depends on the spin-orbit coupling (SOC) values and the corresponding values of the singlet-triplet energy splitting ΔE<sub>Sn-Tm</sub>. Such energy gaps for all the examined compounds can be found in Figure S3 of the Supplementary Materials. The presence of heavy atoms might influence the SOC values by the HAE that, eventually, boosts the probability of ISC. As there are two low-lying singlet excited states and there are four triplets below, all the ISC processes between the S<sub>1-2</sub> and T<sub>1-4</sub> states have been considered for all the compounds under investigation. All the calculated values are collected in Table 3. The two free-base porphyrins **Fused-Fb** and **Meso-Fb** present similar SOC values, with the largest matrix element being 0.4 cm<sup>-1</sup> for the former and 0.3 cm<sup>-1</sup> for the latter. The largest values correspond to the S<sub>2</sub> → T<sub>2</sub> and S<sub>1</sub> → T<sub>1</sub> transitions for **Fused-Fb** and **Meso-Fb**, respectively. These values are similar to that computed for Foscan® (0.24 cm<sup>-1</sup> for the S<sub>1</sub> → T<sub>1</sub> transition), a photosensitizer already used in PDT. The presence of the Zn(II) in the porphyrin core does not provoke any increase in the SOC values for both the **Fused** and **Meso** compounds. The computed values show that the HAE is quite irrelevant for the **Fused** porphyrin conjugated with Au(I), whereas a slight increase was observed for the corresponding **Meso** porphyrin.
<table><caption>**Table 3.** Spin-orbit coupling constants (cm⁻¹) calculated between the two low-lying singlet excited states and the four triplets below them for all the examined systems.</caption><thead><tr><th></th><th>S<sub>1</sub>-T<sub>1</sub></th><th>S<sub>1</sub>-T<sub>2</sub></th><th>S<sub>1</sub>-T<sub>3</sub></th><th>S<sub>1</sub>-T<sub>4</sub></th><th>S<sub>2</sub>-T<sub>1</sub></th><th>S<sub>2</sub>-T<sub>2</sub></th><th>S<sub>2</sub>-T<sub>3</sub></th><th>S<sub>2</sub>-T<sub>4</sub></th></tr></thead><tbody><tr><td>Fused-Fb</td><td>0.3</td><td>0.0</td><td>0.2</td><td>0.0</td><td>0.1</td><td>0.4</td><td>0.0</td><td>0.2</td></tr><tr><td>Fused-Zn</td><td>1.5</td><td>0.1</td><td>0.0</td><td>0.0</td><td>0.1</td><td>1.6</td><td>0.5</td><td>0.0</td></tr><tr><td>Fused-Pd</td><td>17.4</td><td>1.5</td><td>0.3</td><td>5.1</td><td>0.7</td><td>17.1</td><td>5.5</td><td>0.3</td></tr><tr><td>Fused-FbAu</td><td>1.5</td><td>0.3</td><td>0.3</td><td>0.1</td><td>0.1</td><td>1.9</td><td>0.1</td><td>1.2</td></tr><tr><td>Fused-ZnAu</td><td>0.3</td><td>3.2</td><td>1.2</td><td>0.1</td><td>0.3</td><td>0.2</td><td>0.1</td><td>1.1</td></tr><tr><td>Fused-PdAu</td><td>2.2</td><td>18.0</td><td>4.5</td><td>0.1</td><td>16.9</td><td>1.0</td><td>0.1</td><td>6.1</td></tr><tr><td>Meso-Fb</td><td>0.3</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.2</td><td>0.1</td><td>0.2</td></tr><tr><td>Meso-Zn</td><td>0.3</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.1</td><td>0.2</td><td>0.1</td><td>0.2</td></tr><tr><td>Meso-Pd</td><td>0.6</td><td>18.5</td><td>1.1</td><td>0.0</td><td>17.1</td><td>0.5</td><td>0.3</td><td>8.7</td></tr><tr><td>Meso-FbAu</td><td>9.8</td><td>3.3</td><td>2.8</td><td>2.4</td><td>2.2</td><td>9.4</td><td>2.5</td><td>6.2</td></tr><tr><td>Meso-ZnAu</td><td>2.4</td><td>9.9</td><td>1.6</td><td>1.6</td><td>8.8</td><td>2.9</td><td>1.0</td><td>5.0</td></tr><tr><td>Meso-PdAu</td><td>3.4</td><td>12.1</td><td>0.9</td><td>1.8</td><td>12.2</td><td>3.9</td><td>0.6</td><td>5.1</td></tr></tbody></table>

This behavior is expected on the basis of the shapes of the frontier MO contour plots, as illustrated above, as well as that the largest SOC values are obtained for Pd(II) complexes. The largest values of the SOC elements for **Fused-Pd**, which are 17.4 and 17.1 cm⁻¹, are due to the S₁ → T₁ and S₂ → T₂ transitions, respectively. The S₁ → T₂ and S₂ → T₁ transitions correspond to the values of 18.5 and 17.1 cm⁻¹, calculated for **Meso-Pd**. According to the El-Sayed rule, as a high ISC rate can be achieved if the non-radiative transition from the singlet to triplet excited states involves different molecular orbital configurations, a change in the contributing orbitals causes a change in the transitions originating the largest SOCs. The addition of a peripheral Au(I) complex to obtain the **Fused-PdAu** does not increase the SOC amplitudes but causes a change in the nature of the contributing orbitals, and the highest SOC values become 18.0 and 16.9 cm⁻¹, corresponding to the S₁ → T₂ and S₂ → T₁ transitions. Compared to **Meso-Pd**, the amplitude of the SOCs for the corresponding Au(I) complex, **Meso-PdAu**, decreases, becoming 12.1 and 12.2 cm⁻¹, and the involved transitions do not change. In summary, the presence of the external Au(I) complex, increases the SOC values only when conjugated with the **Meso** porphyrin but not enough to have a consistent enhancement of the singlet oxygen generation, whereas no synergistic effect on the SOC is detected when Au(I) is contemporarily present with Zn(II) and Pd(II). Thus, the obtained data show that peripheral Au(I) complexes do not improve the PDT efficiency from a photophysical point of view. However, as previously suggested, the peripheral Au(I) complexes can be beneficial to combine the dark cytotoxic activity of Au(I) with the photodynamic effect of the porphyrin when activated by light of a proper wavelength. Moreover, the Au(I) center can serve as a connecting metal ion to bind ligands that can improve other useful properties, such as water solubility or the active targeting of cancer cells.

### 3. Computational Details

All the calculations were carried out using Gaussian 16 software [50] adopting the density functional theory (DFT) approach. All the geometrical structures were optimized in implicit water solvent using the integral equation formalism variant of the polarizable continuum model (IEFPCM) [51,52] at the B3LYP-D3/6-31G** level of theory [53–55]. The excitation energies were obtained by performing linear response (LR) TDDFT calculations on the optimized structures of all the investigated compounds using the M06 functional [46] and 6-31G** basis sets for all the atoms except Au and Pd, while the effective core potentials SDD [56] were used together with the split valence basis set. This functional was chosen after a preliminary TDDFT benchmark in dichloromethane, the solvent used in the experiments, in order to select the best exchange and correlation functional that is able to reproduce the recorded absorption spectra [40,45] of **Fused-Fb** and **Fused-FbAu**. For this purpose, the structures of the compounds were optimized in CH₂Cl₂, and a few exchange
and correlation functionals were taken into consideration, which were: B3LYP, PBE0 [57], CAM-B3LYP [58], M06, B3PW91 [59], and LC-wPBE [60]. The available experimental absorption spectra of **Meso-Fb**, **Meso-Zn**, and **Meso-ZnAu** were also used to check if the chosen protocol properly described the second class of porphyrins.

Subsequently, in order to mimic the cellular environment, UV spectra and SOC values were calculated in implicit water solvent. The spin-orbit matrix elements were evaluated using SOC-TD-DFT as implemented in the ORCA code [61,62]. A check was carried out to ascertain whether the outcomes of the two Orca and Gaussian codes used were superimposable when the energies of the coupled states were taken into consideration. Thus, relativistic corrections were computed by the zeroth order regular approximation (ZORA) at the ground-state-optimized geometries. Accordingly, the ZORA-DEF2-SVP and SARC-ZORA-SVP basis sets for the main and metal atoms, respectively, were employed, and the SOCs values were calculated as previously reported [13].

**4. Conclusions**

Porphyrins and their derivatives are, very likely, the most studied macrocycles for their use in photodynamic therapy (PDT) to treat different types of cancers because of their proper photophysical properties. A huge number of functionalizations have been tested, aiming at improving the properties of such compounds and enhancing the efficiency in generating singlet oxygen. Here, we have examined the photophysical properties of porphyrins containing peripheral N-heterocyclic carbene (NHC) Au(I) complexes. Both the presence of an external Au(I) ion bound to the NHC moiety and a Zn(II) or Pd(II) ion at the center of the tetrapyrrole ring were taken into consideration. Keeping in mind the underestimation of the absorption wavelengths by the TDDFT approach, the outcomes of our computational analysis show that all the compounds under examination, named **Fused** and **Meso** as a function of the way to link the NHC group to the porphyrin, are able to play their role as photosensitizers in type II PDT processes. Indeed, the maximum absorption wavelengths of the Q-band components fall in the so-called therapeutic window, allowing penetration into the tissue; the singlet-triplet energy gap is higher than that required to generate the cytotoxic singlet oxygen species; the calculated SOC elements for both the radiationless singlet-triplet couplings are similar or larger than that computed for the photosensitizer Foscan® (0.24 cm⁻¹ for the S₁ → T₁ transition), which is already used in PDT. The influence of the presence of the peripheral Au(I) center and of the Zn(II) and Pd(II) ions in the center of the porphyrin core in increasing the rate of the ISC attributable to the heavy atom effect, which causes a change in the SOC elements of two orders of magnitude only when the Pd(II) ion is involved. The presence of the external Au(I) ion, instead, causes an increase in the SOC values only when the NHC ligand is linked to the *meso* position of the porphyrin core (Meso series), but not enough to have a substantial enhancement of the singlet oxygen generation. No synergistic effect on the SOC is detected when Au(I) is contemporarily present with Zn(II) or Pd(II) ions. Plots of the four Gouterman's orbitals confirm that the Q bands are of a π → π* nature, with the MOs localized on the central porphyrin ring and neither the NHC moiety and the Au(I) atom nor the Zn(II) center give a contribution to such orbitals. Only the Pd(II) ion participates in the LUMO and LUMO+1 orbitals, provoking the observed increase in the calculated SOC elements. The outcomes of the present computational analysis show that, even if it is well-established that the introduction of heavy atoms into a photosensitizer should have a direct influence on the rate of the ISC as a consequence of the HAE, no obvious relationship exists since both the identity and the position of the introduced heavy atoms determines the final effect.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27134046/s1, Figure S1: Superimposition of the optimized (green) and crystallographic (silver) structures of **Fused-Fb**. Two selected computed and crystallographic dihedral angles in degrees are also reported.; Table S1: TD-DFT benchmark for **Fused-Fb** and **Fused-FbAu** species optimized at B3LYP/6-31G** level in CH₂Cl₂ implicit solvent; Table S2: M06 maximum wavelengths for the Q and Soret bands calculated for the **Meso-Fb**,
**Meso-Zn**, and **Meso-ZnAu** species compared with the experimental counterparts; Figure S2: Graphi-cal representation of the highest occupied (H and H-1) and the lowest unoccupied (L and L+1) MOs for **Fused-ZnAu**, **Fused-PdAu**, and all the **Meso** compounds; Table S3: Main absorption wavelengths (λ), vertical triplet electronic energies (ΔE), oscillator strengths (f), and main configuration for the examined **Fused** compounds in water solvent computed at M06/6-31G** level of theory; Table S4: Main absorption wavelengths (λ), vertical triplet electronic energies (ΔE), oscillator strengths (f), and main configuration for the examined **Meso** compounds in water solvent computed at M06/6-31G** level of theory; Figure S3: Computed singlet–triplet splitting energies (ΔE, eV) between singlet (S<sub>1</sub>-2) and triplets (T<sub>1-4</sub>) potentially involved in the ISC for (a) **Fused** and (b) **Meso** compounds.

**Author Contributions:** Conceptualization, methodology, software and formal analysis S.S.; investigation F.C.; data curation, validation, resources, supervision, G.M.; writing—review and editing, S.R. and E.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported by the project POR Calabria–FSE/FESR 2014–2020. The University of Calabria and the Calabria Region are acknowledged for financial support and CINECA is acknowledged for the computing time (Project IsC92).

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Sample Availability:** Samples of the compounds are not available from the authors.

**References**

1. Zhou, Z.; Song, J.; Nie, L.; Chen, X. Reactive Oxygen Species Generating Systems Meeting Challenges of Photodynamic Cancer Therapy. *Chem. Soc. Rev.* **2016**, *45*, 6597–6626. [CrossRef] [PubMed]

2. Han, S.-K.; Sik, R.H.; Motten, A.G.; Chignell, C.F.; Bilski, P.J. Photosensitized Oxidation of Tetrabromobisphenol A by Humic Acid in Aqueous Solution. *Photochem. Photobiol.* **2009**, *85*, 1299–1305. [CrossRef]

3. Romero, O.C.; Straub, A.P.; Kohn, T.; Nguyen, T.H. Role of Temperature and Suwannee River Natural Organic Matter on Inactivation Kinetics of Rotavirus and Bacteriophage MS2 by Solar Irradiation. *Environ. Sci. Technol* **2011**, *45*, 10385–10393. [CrossRef]

4. Kohn, T.; Nelson, K.L. Sunlight-Mediated Inactivation of MS2 Coliphage via Exogenous Singlet Oxygen Produced by Sensitizers in Natural Waters. *Environ. Sci. Technol.* **2007**, *41*, 192–197. [CrossRef]

5. Polat, E.; Kang, K. Natural Photosensitizers in Antimicrobial Photodynamic Therapy. *Biomedicines* **2021**, *9*, 584. [CrossRef] [PubMed]

6. Cieplik, F.; Deng, D.; Crielaard, W.; Buchalla, W.; Hellwig, E.; Al-Ahmad, A.; Maisch, T. Antimicrobial Photodynamic Therapy—What We Know and What We Don’t. *Crit. Rev. Microbiol.* **2018**, *44*, 571–589. [CrossRef] [PubMed]

7. Songca, S.P.; Adjei, Y. Applications of Antimicrobial Photodynamic Therapy against Bacterial Biofilms. *Int. J. Mol. Sci.* **2022**, *23*, 3209. [CrossRef]

8. Kwiatkowski, S.; Knap, B.; Przystupski, D.; Saczko, J.; Kędzierska, E.; Knap-Czop, K.; Kotlińska, J.; Michel, O.; Kotowski, K.; Kulbacka, J. Photodynamic Therapy—Mechanisms, Photosensitizers and Combinations. *Biomed. Pharmacother.* **2018**, *106*, 1098–1107. [CrossRef]

9. Dabrowski, J.M. Reactive Oxygen Species in Photodynamic Therapy: Mechanisms of Their Generation and Potentiation. In *Advances in Inorganic Chemistry*; Academic Press: Cambridge, MA, USA, 2017; Volume 70, pp. 343–394.

10. Alberto, M.E.; De Simone, B.C.; Liuzzi, S.; Marino, T.; Russo, N.; Toscano, M. Iodine Substituted Phosphorus Corrole Complexes as Possible Photosensitizers in Photodynamic Therapy: Insights from Theory. *J. Comput. Chem.* **2020**, *41*, 1395–1401. [CrossRef]

11. Zhou, C.H.; Zhao, X. Chemically Modified Fullerene Derivatives as Photosensitizers in Photodynamic Therapy: A First-Principles Study. *J. Comput. Chem.* **2012**, *33*, 861–867. [CrossRef] [PubMed]

12. Plaetzer, K.; Krammer, B.; Berlanda, J.; Berr, F.; Kiesslich, T. Photophysics and Photochemistry of Photodynamic Therapy: Fundamental Aspects. *Lasers Med. Sci.* **2009**, *24*, 259–268. [CrossRef]

13. De Simone, B.C.; Mazzone, G.; Russo, N.; Sicilia, E.; Toscano, M. Metal Atom Effect on the Photophysical Properties of Mg(II), Zn(II), Cd(II), and Pd(II) Tetraphenylporphyrin Complexes Proposed as Possible Drugs in Photodynamic Therapy. *Molecules* **2017**, *22*, 1093. [CrossRef] [PubMed]

14. Banfi, S.; Caruso, E.; Caprioli, S.; Mazzagatti, L.; Canti, G.; Ravizza, R.; Gariboldi, M.; Monti, E. Photodynamic Effects of Porphyrin and Chlorin Photosensitizers in Human Colon Adenocarcinoma Cells. *Bioorganic Med. Chem.* **2004**, *12*, 4853–4860. [CrossRef]

15. Janas, K.; Boniewska-Bernacka, E.; Dyrda, G.; Slota, R. Porphyrin and Phthalocyanine Photosensitizers Designed for Targeted Photodynamic Therapy of Colorectal Cancer. *Bioorganic Med. Chem.* **2021**, *30*, 115926. [CrossRef]
16. Amos-Tautua, B.; Songca, S.; Oluwafemi, O. Application of Porphyrins in Antibacterial Photodynamic Therapy. *Molecules* **2019**, *24*, 2456. [CrossRef]

17. Ethirajan, M.; Chen, Y.; Joshi, P.; Pandey, R.K. The Role of Porphyrin Chemistry in Tumor Imaging and Photodynamic Therapy. *Chem. Soc. Rev.* **2010**, *40*, 340–362. [CrossRef] [PubMed]

18. Wu, F.; Yang, M.; Zhang, J.; Zhu, S.; Shi, M.; Wang, K. Metalloporphyrin-Indomethacin Conjugates as New Photosensitizers for Photodynamic Therapy. *JBIC J. Biol. Inorg. Chem.* **2019**, *24*, 53–60. [CrossRef]

19. Chilakamarthi, U.; Giribabu, L. Photodynamic Therapy: Past, Present and Future. *Chem. Rec.* **2017**, *17*, 775–802. [CrossRef] [PubMed]

20. Zhao, X.; Liu, J.; Fan, J.; Chao, H.; Peng, X. Recent Progress in Photosensitizers for Overcoming the Challenges of Photodynamic Therapy: From Molecular Design to Application. *Chem. Soc. Rev.* **2021**, *50*, 4185–4219. [CrossRef]

21. Huang, H.; Song, W.; Rieffel, J.; Lovell, J.F. Emerging Applications of Porphyrins in Photomedicine. *Front. Phys.* **2015**, *3*, 23. [CrossRef]

22. Militello, M.P.; Hernández-Ramírez, R.E.; Lijanova, I.V.; Previtali, C.M.; Bertolotti, S.G.; Arbeloa, E.M. Novel PAMAM Dendrimers with Porphyrin Core as Potential Photosensitizers for PDT Applications. *J. Photochem. Photobiol. A Chem.* **2018**, *353*, 71–76. [CrossRef]

23. O'Connor, A.E.; Gallagher, W.M.; Byrne, A.T. Porphyrin and Nonporphyrin Photosensitizers in Oncology: Preclinical and Clinical Advances in Photodynamic Therapy. *Photochem. Photobiol.* **2009**, *85*, 1053–1074. [CrossRef] [PubMed]

24. Montaseri, H.; Kruger, C.A.; Abrahamse, H. Recent Advances in Porphyrin-Based Inorganic Nanoparticles for Cancer Treatment. *Int. J. Mol. Sci.* **2020**, *21*, 3358. [CrossRef] [PubMed]

25. Tsolekile, N.; Nelana, S.; Oluwafemi, O.S. Porphyrin as Diagnostic and Therapeutic Agent. *Molecules* **2019**, *24*, 2669. [CrossRef]

26. Lammer, A.D.; Cook, M.E.; Sessler, J.L. Synthesis and Anti-Cancer Activities of a Water Soluble Gold(III) Porphyrin. *J. Porphyr. Phthalocyanines* **2015**, *19*, 398–403. [CrossRef]

27. Martinez De Pinillos Bayona, A.; Mroz, P.; Thunshelle, C.; Hamblin, M.R. Design Features for Optimization of Tetrapyrrole Macrocycles as Antimicrobial and Anticancer Photosensitizers. *Chem. Biol. Drug Des.* **2017**, *89*, 192–206. [CrossRef]

28. Robinson, G.W.; Frosch, R.P. Electronic Excitation Transfer and Relaxation. *J. Chem. Phys.* **2004**, *38*, 1187. [CrossRef]

29. Abell, G.C.; Mozumder, A. Theory of Intersystem Crossing in Aromatic Hydrocarbons. *J. Chem. Phys.* **2003**, *56*, 4058. [CrossRef]

30. Horiuchi, H.; Terashima, K.; Sakai, A.; Suda, D.; Yoshihara, T.; Kobayashi, A.; Tobita, S.; Okutsu, T. The Effect of Central Metal on the Photodynamic Properties of Silylated Tetraphenylporphyrin Derivative. *J. Photochem. Photobiol. A Chem.* **2016**, *321*, 72–78. [CrossRef]

31. Deng, J.; Li, H.; Yang, M.; Wu, F. Palladium Porphyrin Complexes for Photodynamic Cancer Therapy: Effect of Porphyrin Units and Metal. *Photochem. Photobiol. Sci.* **2020**, *19*, 905–912. [CrossRef]

32. Toubia, I.; Nguyen, C.; Diring, S.; Ali, L.M.A.; Larue, J.; Aoun, R.; Frochot, C.; Gary-Bobo, M.; Kobeissi, M.; Odobel, F. Synthesis and Anticancer Activity of Gold Porphyrin Linked to Malonate Diamine Platinum Complexes. *Inorg. Chem.* **2019**, *58*, 12395–12406. [CrossRef] [PubMed]

33. Naik, A.; Rubbiani, R.; Gasser, G.; Spingler, B. Visible-Light-Induced Annihilation of Tumor Cells with Platinum-Porphyrin Conjugates. *Angew. Chem.* **2014**, *126*, 7058–7061. [CrossRef]

34. Lottner, C.; Knuechel, R.; Bernhardt, G.; Brunner, H. Distribution and Subcellular Localization of a Water-Soluble Hematoporphyrin–Platinum(II) Complex in Human Bladder Cancer Cells. *Cancer Lett.* **2004**, *215*, 167–177. [CrossRef]

35. Lottner, C.; Knuechel, R.; Bernhardt, G.; Brunner, H. Combined Chemotherapeutic and Photodynamic Treatment on Human Bladder Cells by Hematoporphyrin–Platinum(II) Conjugates. *Cancer Lett.* **2004**, *203*, 171–180. [CrossRef]

36. Brunner, H.; Schellerer, K.-M. New Porphyrin Platinum Conjugates for the Cytostatic and Photodynamic Tumor Therapy. *Inorg. Chim. Acta* **2003**, *350*, 39–48. [CrossRef]

37. Lottner, C.; Bart, K.-C.; Bernhardt, G.; Brunner, H. Hematoporphyrin-Derived Soluble Porphyrin–Platinum Conjugates with Combined Cytotoxic and Phototoxic Antitumor Activity. *J. Med. Chem.* **2002**, *45*, 2064–2078. [CrossRef] [PubMed]

38. Brunner, H.; Schellerer, K.-M.; Treittinger, B. Synthesis and in Vitro Testing of Hematoporphyrin Type Ligands in Platinum(II) Complexes as Potent Cytostatic and Phototoxic Antitumor Agents. *Inorg. Chim. Acta* **1997**, *264*, 67–79. [CrossRef]

39. Longevial, J.; Rose, C.; Poyac, L.; Clément, S.; Richeter, S. Molecular Systems Combining Porphyrinoids and N-Heterocyclic Carbenes. *Eur. J. Inorg. Chem.* **2021**, *2021*, 776–791. [CrossRef]

40. Longevial, J.F.; Langlois, A.; Buisson, A.; Devillers, C.H.; Clément, S.; van der Lee, A.; Harvey, P.D.; Richeter, S. Synthesis, Characterization, and Electronic Properties of Porphyrins Conjugated with N-Heterocyclic Carbene (NHC)–Gold(I) Complexes. *Organometallics* **2016**, *35*, 663–672. [CrossRef]

41. Rose, C.; Lichon, L.; Daurat, M.; Clément, S.; Gary-Bobo, M.; Richeter, S. In Vitro Toxicity and Photodynamic Properties of Porphyrinoids Bearing Imidazolium Salts and N-Heterocyclic Carbene Gold(I) Complexes. *Comptes Rendus. Chim.* **2021**, *24*, 83–99. [CrossRef]

42. Longevial, J.-F.; El Cheikh, K.; Aggad, D.; Lebrun, A.; van der Lee, A.; Tielens, F.; Clément, S.; Morère, A.; Garcia, M.; Gary-Bobo, M.; et al. Porphyrins Conjugated with Peripheral Thiolato Gold(I) Complexes for Enhanced Photodynamic Therapy. *Chem.-Eur. J.* **2017**, *23*, 14017–14026. [CrossRef] [PubMed]

43. Scattolin, T.; Lippmann, P.; Beliš, M.; Hecke, K.; Ott, I.; Nolan, S.P. A Simple Synthetic Entryway into (N-heterocyclic Carbene)Gold-steroidyl Complexes and Their Anticancer Activity. *Appl. Organomet. Chem.* **2022**, e6624. [CrossRef]
44. Gouterman, M. Spectra of Porphyrins. *J. Mol. Spectrosc.* **1961**, *6*, 138–163. [CrossRef]

45. Lefebvre, J.F.; Leclercq, D.; Gisselbrecht, J.P.; Richeter, S. Synthesis, Characterization, and Electronic Properties of Metalloporphyrins Annulated to Exocyclic Imidazole and Imidazolium Rings. *Eur. J. Org. Chem.* **2010**, *2010*, 1912–1920. [CrossRef]

46. Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. *Theor. Chem. Acc.* **2008**, *120*, 215–241. [CrossRef]

47. Alberto, M.E.; Mazzone, G.; Regina, C.; Russo, N.; Sicilia, E. Theoretical Exploration of the Photophysical Properties of Two-Component Ru II–Porphyrin Dyes as Promising Assemblies for a Combined Antitumor Effect. *Dalton Trans.* **2020**, *49*, 12653–12661. [CrossRef]

48. Alberto, M.E.; De Simone, B.C.; Mazzone, G.; Marino, T.; Russo, N. Photophysical Properties of Free and Metallated Meso-Substituted Tetrabenzotriazaporphyrin from Density Functional Theory Investigation. *Dye. Pigment.* **2015**, *120*, 335–339. [CrossRef]

49. Mazzone, G.; Russo, N.; Sicilia, E. Theoretical Investigation of the Absorption Spectra and Singlet-Triplet Energy Gap of Positively Charged Tetraphenylporphyrins as Potential Photodynamic Therapy Photosensitizers. *Can. J. Chem.* **2013**, *91*, 902–906. [CrossRef]

50. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. *G16_C01*, *Gaussian 16*, *Revision C.01*; Gaussian, Inc.: Wallingford, CT, USA, 2016.

51. Miertuš, S.; Tomasi, J. Approximate Evaluations of the Electrostatic Free Energy and Internal Energy Changes in Solution Processes. *Chem. Phys.* **1982**, *65*, 239–245. [CrossRef]

52. Miertuš, S.; Scrocco, E.; Tomasi, J. Electrostatic Interaction of a Solute with a Continuum. A Direct Utilizaion of AB Initio Molecular Potentials for the Prevision of Solvent Effects. *Chem. Phys.* **1981**, *55*, 117–129. [CrossRef]

53. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. *J. Chem. Phys.* **2010**, *132*, 154104. [CrossRef]

54. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. *Phys. Rev. B* **1988**, *37*, 785. [CrossRef]

55. Becke, A.D. Density-functional Thermochemistry. III. The Role of Exact Exchange. *J. Chem. Phys.* **1998**, *98*, 5648. [CrossRef]

56. Andrae, D.; Häußermann, U.; Dolg, M.; Stoll, H.; Preuß, H. Energy-Adjustedab Initio Pseudopotentials for the Second and Third Row Transition Elements. *Theor. Chim. Acta* **1990**, *77*, 123–141. [CrossRef]

57. Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. *J. Chem. Phys.* **1999**, *110*, 6158. [CrossRef]

58. Yanai, T.; Tew, D.P.; Handy, N.C. A New Hybrid Exchange-Correlation Functional Using the Coulomb-Attenuating Method (CAM-B3LYP). *Chem. Phys. Lett.* **2004**, *393*, 51–57. [CrossRef]

59. Perdew, J.P.; Wang, Y. Accurate and Simple Analytic Representation of the Electron-Gas Correlation Energy. *Phys. Rev. B* **1992**, *45*, 13244. [CrossRef]

60. Henderson, T.M.; Izmaylov, A.F.; Scalmani, G.; Scuseria, G.E. Can Short-Range Hybrids Describe Long-Range-Dependent Properties? *J. Chem. Phys.* **2009**, *131*, 044108. [CrossRef]

61. Neese, F. Software Update: The ORCA Program System—Version 5.0. *Wiley Interdiscip. Rev. Comput. Mol. Sci.* **2022**, e1606. [CrossRef]

62. Neese, F.; Wiley, J. The ORCA Program System. *Wiley Interdiscip. Rev. Comput. Mol. Sci.* **2012**, *2*, 73–78. [CrossRef]