Next Article in Journal
Nitrergic Enteric Neurons in Health and Disease—Focus on Animal Models
Previous Article in Journal
Temperature Stress and Responses in Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rational Design of Modified Oxobacteriochlorins as Potential Photodynamic Therapy Photosensitizers

Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria, 87036 Rende, Italy
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(8), 2002; https://doi.org/10.3390/ijms20082002
Submission received: 4 April 2019 / Revised: 18 April 2019 / Accepted: 20 April 2019 / Published: 24 April 2019
(This article belongs to the Section Molecular Pathology, Diagnostics, and Therapeutics)

Abstract

:
The modulation of the photophysical properties of a series of recently synthetized oxobacteriochlorins with the introduction of heavy atoms in the macrocycles, was investigated at density functional level of theory and by means of the time-dependent TDDFT formulation. Absorption frequencies, singlet-triplet energy gaps and spin-orbit coupling (SOC) constants values were computed for all the investigated compounds. Results show how the sulfur- selenium- and iodine-substituted compounds possess improved properties that make them suitable for application in photodynamic therapy (PDT).

Graphical Abstract

1. Introduction

Photodynamic therapy (PDT) is used in cancer treatment due to its mini-invasive nature, limited side effects, easy procedure and short treatment time in respect to surgery or chemotherapy [1,2,3,4,5]. PDT is based upon the so-called photodynamic effect, that can be defined as the tissue destruction or damage by the combined action of a photosensitizer, a visible or near visible light, and molecular oxygen. The key cytotoxic agent is the highly reactive singlet exited oxygen (O2, 1Δg) produced in a photophysical cycle (type II photoexcitation process) which can be briefly described as it follows [6]: Absorption of light with the appropriate wavelength by a photosensitizer in its ground state (S0) leads to the excitation of one electron into a higher-energy orbital. The excited singlet (S1) that is formed may undergo an ‘intersystem crossing’ (ISC) to generate a more stable excited triplet state (T1). The long lifetime of the triplet state allows to transfer its energy to the molecular oxygen (O2), which passes from its triplet ground state (3Σg) to the singlet excited state (1Δg). The efficiency of the ISC process strictly depends on the amplitude of the spin-orbit coupling singlet constants (SOC) [7] and it is possible only if the lowest triplet excited state (T1) lies above the energy needed to produce the molecular singlet oxygen (0.98 eV). In order to better penetrate into the tissue, the photosensitizer absorption wavelength should fall in the so called PDT absorption window that is between 500 to 800 nm. To improve the penetration of the light to the tissue, many scientists work on excitation with longer wavelength at near-IR or even longer wavenumbers [5]. In this spectroscopic region, even if the light deeply penetrates into the tissues it does not produce side effects due to the high frequency vibrations. In addition, the photosensitizer must be non-toxic in the dark, not subjected to intermolecular aggregation phenomena (which would diminish its phototherapeutic efficacy), and soluble in aqueous environment.
Although research in this area has increased significantly in recent years [4,5,6], the number of molecules, that have passed clinical trials to be used in medical protocols, are few and essentially based on porphyrin-like structures. Among these, we mention Foscan® (Temoporfin) [8], Photofrin® (porfimer sodium) [9], and Laserphyrin® (Talaporfin) [10]. Second generation photosensitizers such as Lutrin® (Lutetium texaphyrin) are in different clinical phases [11]. The success of the porphyrin-like systems as PDT agents is due to the fact that, as it is required, they are not toxic, easily soluble in physiological solutions and able to absorb into the spectral region of the therapeutic window. Furthermore, their long-wavelength absorption Q band can be tuned easily introducing different chemical groups at distinct sites around the perimeter of the hydroporphyrin skeleton. In particular, the presence of auxochromes positioned along the bisection axis of pyrrole rings causes batochromic shifts with respect to the unsubstituted chorines or bacteriochlorins.
Recently, a series of oxo- and thioxo- substituted bacteriochlorins [12,13] have been synthetized and their photophysical properties determined. Moreover, other photosensitizers with different chemical structures, squaraine [14], BODIPY [15], aza-BODIPY [16] and BOIMPY [17] have been proposed for their use in PDT on the basis of cytotoxicity tests. Some of these works underline the possibility to redshift the absorption Q band as a function of the substituents position in the macrocyclic skeleton, being this, one of the most important requirements to satisfy in designing new PDT drugs.
In this scenario, modern methods of theoretical and computational chemistry can offer significant aids for a reliable and rational design of targeted drugs by shortening the usually long chain from chemical synthesis to clinical approval. In particular, density functional based methods (DFT) are able to predict the structural, energetic and spectroscopic properties of a great number of complex chemical systems conjugating reasonable computational effort with chemical accuracy.
In this work we will show how, by using DFT, it is possible not only to correctly reproduce a series of photophysical properties, but also predict their change as a result of a targeted substitution of atoms or functional groups. Starting from some recently synthesized oxobacteriochlorins [12], the calculation of their absorption spectra, singlet-triplet energy gaps and spin-orbit couplings constants allowed to suggest chemical modifications that fulfill the criteria necessary to propose them as photosensitizers in photodynamic therapy.

2. Results and Discussion

The ground and excited singlet and triplet states optimized structures (see Table S1 for Cartesian coordinates) present a planar conformation of the hydroporphyrin macrocycle in all studied systems (see Scheme 1).
In order to evaluate the effects of the heavy atom on the photophysical properties of the considered systems, vertical excitation energies at TDDFT level of theory have been calculated (see Table 1 and Figure S1). From the obtained values it results a batochromic shift in going from compounds 1, with an oxo group (599 nm), to 2 (623 nm) and 3 (678 nm) in which S and a Se atoms, respectively substitute the oxygen one.
The dioxo molecule (4) shows a Q band at 602 nm meaning that the addition of a second oxo group does not change the absorption behavior significantly. The same situation is found in systems 7 and 9, where two phenyl groups are present in the phorphyrin-like skeleton. In all systems, the transitions involve the HOMO and LUMO orbitals and are π-π in nature. This can be inferred from the molecular orbital pictures of all the three compounds, in which a π delocalization on the entire macrocyclic ring (see Figure 1) is present. Several triplet states have been localized above the low-lying excited singlet one. In particular, three triplet states have been found in unsubstituted molecules (1, 4–7, 9) while four ones appeared when a sulfur or a selenium atom replaces the oxo group (2, 3, 8). T1 derives from a HOMO-LUMO transition while T2 and T3 are originated by an electron transfer involving different orbitals (see Table 1). T4, present in the thioxi and selenoxi containing molecules, is almost degenerate with the S1 state. For all systems, the energy gap between the ground state and the low-lying triplet (T1) assumes values higher than that required to activate molecular oxygen from the 3Σg electronic state to 1Δg one (0.98 eV) and therefore satisfies one of the criteria necessary for an efficient PDT photosensitizer (see Table 1).
As previously mentioned, the photosensitizers active in PDT must possess an efficient singlet-triplet intersystem crossing whose kinetic constant can be derived from the Fermi golden rule throughout the following expression [7]:
k I S C ( S 1 T j ) = 2 π S 1 | H S O | T j 2 4 π λ k B T n S k n n ! e x p ( S k ) e x p [ ( Δ E S T j + n ω k + λ s ) 2 4 λ k B T ]
where S 1 | H S O | T j are the spin-orbit matrix elements, ωk the high-frequency vibrational modes and ΔES-T is energy gap between the S1 and Tj states at their equilibrium geometries. In this formulation, the Franck–Condon-weighted density of states is estimated by the Marcus theory for the non-radiative transition from S1 to each triplet excited state (Tj) lying below.
From the reported equation, it can be inferred how the ISC efficiency depends mainly on the values of the spin orbit couplings other than the singlet-triplet energy gaps for the involved ISC. For all the investigated compounds, we have decided to compute the SOCs concerning the S1Tj transitions. The same procedure has been adopted for the selenium containing molecules 3 and 8 for which the S2 state is populated (see Table 1). This should mean that, in agreement with the Kasha rules [18], after the excitation, which populates the S2 state, this state decays on S1 through a fast internal conversion process (IC).
The computed S 1 | H S O | T j , together with the relative ΔES1-Tj energy gaps are collected in Table 2.
From that Table it clearly appears that the presence of heavy atoms in the porphyrin-like skeleton greatly influences the amplitude of the spin orbit coupling constants. In fact, in going from system 1 (oxo), to 2 (thioxo), to 3 (selenoxo), the S 1 | H S O | T 1 increases significantly assuming values of 114.7 and 685.8 cm−1 for 2 and 3, respectively. The same behavior is found for compounds 7 (0.2 cm−1) and 8 (675.4 cm−1). It is interesting to note that for these two last compounds, having a T4 state with energy similar to that of S1, the SOCs values for the S1T4 transitions have quite high values (52.2 and 239.5 cm−1 for 2 and 3, respectively).
The presence of halogen atoms also contributes to increase the S 1 | H S O | T 1 matrix element, especially in the case of iodine (17.4 cm−1). While for the other examined systems the values of the other S1–Tj SOCs are negligible, for the iodine-containing system (6) the S 1 | H S O | T 2 and S 1 | H S O | T 3 matrix elements increase considerably compared to those of the unsubstituted system (see Table 2).
These trends can be rationalized using the El Sayed rules [19] that establish that SOC values increase significantly if a symmetry change of the molecular orbitals involved in the transitions, occurs. The orbital pictures of singlet (S1) and triplets (T1–T4) for systems 1 and 3 are reported in Figure 1.
For the unsubstituted oxobacteriochlorin (see Table 1) S1 and T1 transition involves the HOMO-LUMO orbitals that are both π-π * in nature. Although for T1 and T2 the transitions are due to HOMO-1→LUMO and HOMO→LUMO+1, respectively, Figure 1 shows that all these orbitals are of π type. The presence of selenium in 3 affects significantly the ISC in the S1→T1 and S1→T4 transitions since the orbital compositions change becoming nπ→ππ with different percentage (Table 1). In accordance with the El Sayed rules [19], the SOCs values should be high. In other considered intersystem crossings, the transition involves orbitals with almost the same symmetry in agreement with the low SOC values obtained by our calculations.
Same correlations between orbital nature and SOC amplitudes have been previously found for S, Se and Te-substituted deoxyguanosines [20] and thymidines [21]. Furthermore, the heavy atom effects have been explored in the case of halogen substituted aza-BODIPY [22], BOIMPY [23] and metal containing tetraphenylporphyrin [24].
Considering the computed spin-orbit coupling constants, we can suggest the following possible faster deactivation paths for the considered ISC processes:
  • S 0 A S 2 IC S 1 ISC T 1 or S 0 A S 2 IC S 1 ISC T 4   IC T 1 for compound 3;
  • S 0 A S 1 ISC T 1 or S 0 A S 1 ISC T 4   IC T 1 for compounds 2 and 8;
  • S 0 A S 1 ISC T 1 or S 0 A S 1 ISC T 2   IC T 1 for compound 6;
  • S 0 A S 1 ISC T 1 for all other systems.
An accurate prediction of k I S C ( S 1 T j ) requires that the adiabatic singlet-triplet energy gaps are taken into consideration as well as the calculation of the Franck–Condon weighted density of states. Work is in progress in this direction.

3. Materials and Methods

All the computations have been carried out by using the Gaussian 09 [25] code employing DFT and its time dependent formulation TDDFT [26]. The hybrid B3LYP [27,28] exchange–correlation functional has been used in conjunction with the 6-31+G* basis set for all the atoms and the SDD pseudopotential [29] for iodine. No symmetry constraints have been imposed during the geometry optimizations.
Computations have been performed in toluene to mimic the experimental environment, by using a dielectric constant of ε = 2.37. The integral equation formalism for the polarizable continuum model (IEFPCM) [30], which corresponds to a linear response in nonequilibrium solvation, has been employed. Spin–orbit matrix elements S 1 | H S O | T j have been computed with the DALTON code [31] by using the spin–orbit coupling operators for effective core potentials with an effective nuclear charge [32] for iodine atoms containing systems. The atomic mean field approximation [33] has been used in the other cases. For this purpose, the B3LYP functional has been employed in conjunction with the cc-pVDZ basis set for all the atoms except for iodine for which the corresponding pseudopotential has been adopted.
The choice of the computational protocol has been suggested by our previous benchmarks in computing different photophysical properties for several compounds [34,35,36,37,38,39].

4. Conclusions

Stimulated by the recent synthesis of a series of oxo- and dioxobacteriochlorine that should produce singlet oxygen under UV-Vis irradiation, we have computed the main photophysical properties of such compounds relevant in photodynamic therapy. In addition, we have demonstrated that the inclusion of a heavy atom (S, Se, I) in the skeleton of these systems strongly improves the spin-orbit coupling constant values, and consequently the efficiency of the intersystem crossing, ensuring high yield of the cytotoxic singlet molecular oxygen agent.
All the studied systems show absorption wavelengths that fall in the so-called photodynamic window and possess a singlet-triplet energy gap higher than that necessary to generate 1Δg O2 (0.98 eV).
The computed SOCs have values higher than that of Foscan® already used as photosensitizer in the PDT clinical protocols.
We hope that our work can stimulate the synthesis and the characterization of new possible PDT drugs.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/8/2002/s1.

Author Contributions

M.E.A., B.C.D.S., E.S., M.T. and N.R. made equal contributions to the study and the publication of this work.

Funding

This research received no external funding.

Acknowledgments

We thank the Dipartimento di Chimica e Tecnologie Chimiche, Università della Calabria for financial support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Agostinis, P.; Berg, K.; Cengel, K.A.; Foster, T.H.; Girotti, A.W.; Gollnick, S.O.; Hahn, S.M.; Hamblin, M.R.; Juze-niene, A.; Kessel, D.; et al. Photodynamic Therapy of Cancer: An Update. CA Cancer J. Clin. 2011, 61, 250–281. [Google Scholar] [CrossRef] [PubMed]
  2. Dolmans, D.E.; Fukumura, D.; Jain, R.K. Photodynamic therapy for cancer. Nat. Rev. Cancer 2003, 3, 380–387. [Google Scholar] [CrossRef]
  3. Allison, R.R.; Bagnato, V.S.; Sibata, C.H. Future of oncologic photodynamic therapy. Future Oncol. 2010, 6, 929–940. [Google Scholar] [CrossRef] [Green Version]
  4. Yano, S.; Hirohara, S.; Obata, M.; Hagiya, Y.; Ogura, S.; Ikeda, I.; Kataoka, H.; Tanaka, M.; Joh, T. Current states and future views in photodynamic therapy. J. Photochem. Photobiol. C 2011, 12, 46–67. [Google Scholar] [CrossRef]
  5. Dabrowski, J.M.; Arnaut, L.G. Photodynamic therapy (PDT) of cancer: From local to systemic treatment. Photochem Photobiol Sci. 2015, 14, 1765–1780. [Google Scholar] [CrossRef] [PubMed]
  6. Zhu, J.; Dominijanni, A.; Rodríguez-Corrales, J.Á.; Prussin, R.; Zhao, Z.; Li, T.; Robertson, J.L.; Brewer, K.J. Visible light-induced cytotoxicity of Ru,Os–polyazine complexes towards rat malignant glioma. Inorg. Chim. Acta 2017, 454, 155–161. [Google Scholar] [CrossRef]
  7. Marian, C.M. Spin-orbit coupling and intersystem crossing in molecules. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 187–203. [Google Scholar] [CrossRef]
  8. 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. Bioorg. Med. Chem. 2004, 12, 4853–4860. [Google Scholar] [CrossRef]
  9. García-Díaz, M.; Sánchez-García, D.; Soriano, J.; Sagrista, L.; Mora, M.; Villanueva, A.; Stockert, J.; Cañete, M.; Nonell, S. Temocene: The porphycene analogue of temoporfin (Foscan®). MedChemComm 2011, 2, 616–619. [Google Scholar]
  10. Nanashima, A.; Abo, T.; Nonaka, T.; Nonaka, Y.; Morisaki, T.; Uehara, R.; Ohnita, K.; Fukuda, D.; Murakami, G.; Tou, K.; et al. Photodynamic therapy using talaporfin sodium (Laserphyrin®) for bile duct carcinoma: A preliminary clinical trial. Anticancer. Res. 2012, 32, 4931–4938. [Google Scholar] [PubMed]
  11. 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. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Liu, M.; Chen, C.-Y.; Hood, D.; Taniguchi, M.; Diers, J.R.; Bocian, D.F.; Holten, D.; Lindsey, J.S. Synthesis, photophysics and electronic structure of oxobacteriochlorins. New J. Chem. 2017, 41, 3732–3744. [Google Scholar] [CrossRef]
  13. Tamiaki, H.; Xu, M.; Kinoshita, Y. Synthesis of oxo-, thioxo- and methylene-substituted bacteriochlorins by modifying chlorophyll-a and their electronic absorption in visible and near-infrared regions. J. Photochem. Photobiol. A Chem. 2013, 252, 60–68. [Google Scholar] [CrossRef]
  14. Kamkaew, A.; Lim, S.H.; Lee, H.B.; Kiew, L.V.; Chung, L.Y.; Burgess, K. BODIPY dyes in photodynamic therapy. Chem. Soc. Rev. 2013, 42, 77–88. [Google Scholar] [CrossRef] [PubMed]
  15. Ramaiah, D.; Eckert, I.; Arun, K.T.; Weidenfeller, L.; Epe, B. Squaraine dyes for photodynamic therapy: Study of their cytotoxicity and genotoxicity in bacteria and mammalian cells. Photochem. Photobiol. 2002, 76, 672–677. [Google Scholar] [CrossRef]
  16. Loudet, A.; Burgess, K. BODIPY dyes and their derivatives: Syntheses and spectroscopic properties. Chem. Rev. 2007, 107, 4891–4932. [Google Scholar] [CrossRef]
  17. Patalag, L.J.; Jones, P.G.; Werz, D.B. BOIMPYs: Rapid Access to a Family of Red-Emissive Fluorophores and NIR Dyes. Angew. Chem., Int. Ed. 2016, 55, 13340–13344. [Google Scholar] [CrossRef] [PubMed]
  18. Kasha, M. Characterization of electronic transitions in complex molecules. Discuss Faraday Soc. 1950, 9, 14–19. [Google Scholar] [CrossRef]
  19. El-Sayed, M.A. Triplet state. Its radiative and nonradiative properties. Acc. Chem. Res. 1968, 1, 8–16. [Google Scholar] [CrossRef]
  20. Pirillo, J.; Mazzone, G.; Russo, N.; Bertini, L. Photophysical properties of S, Se and te-substituted deoxyguanosines: Insight into their ability to act as chemotherapeutic agents. J. Chem. Inf. Model 2017, 57, 234–242. [Google Scholar] [CrossRef]
  21. Pirillo, J.; De Simone, B.C.; Russo, N. Photophysical properties prediction of seleniumand tellurium-substituted thymidine as potential UVA chemotherapeutic agents. Theor. Chem. Acc. 2016, 135, 8. [Google Scholar] [CrossRef]
  22. De Simone, B.C.; Mazzone, G.; Pirillo, J.; Russo, N.; Sicilia, E. Halogen atom effect on the photophysical properties of substituted aza-BODIPY derivatives. Phys. Chem. Chem. Phys. 2017, 19, 2530–2536. [Google Scholar] [CrossRef]
  23. De Simone, B.C.; Mazzone, G.; Russo, N.; Sicilia, E.; Toscano, M. Excitation energies, singlet–triplet energy gaps, spin–orbit matrix elements and heavy atom effects in BOIMPYs as possible photosensitizers for photodynamic therapy: A computational investigation. Phys. Chem. Chem. Phys. 2018, 20, 2656–2661. [Google Scholar] [CrossRef] [PubMed]
  24. 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. [Google Scholar] [CrossRef] [PubMed]
  25. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, USA, 2010; Available online: http://gaussian.com/ (accessed on 23 April 2019).
  26. Casida, M.E. Recent Advances in Density Functional Methods, Part I; Chong, D.P., Ed.; World Scientific: Singapore, 1995; pp. 155–192. [Google Scholar]
  27. Becke, A.D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  28. 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 Condens. Matter Mater. Phys. 1988, 37, 785–789. [Google Scholar] [CrossRef]
  29. Andrae, D.; Haussermann, U.; Dolg, M.; Stoll, H.; Preuss, H. Energy-adjusted ab initio pseudopotentials for the second and third row transition elements. Theor. Chim. Acta 1990, 77, 123–141. [Google Scholar] [CrossRef]
  30. Cossi, M.; Barone, V. Solvent effect on vertical electronic transitions by the polarizable continuum model. J. Chem. Phys. 2000, 112, 2427–2435. [Google Scholar] [CrossRef]
  31. Dalton. A Molecular Electronic Structure Program. Release Dalton. 2011. Available online: http://daltonpr ogram.org/ (accessed on 24 February 2016).
  32. Ågren, H.; Vahtras, O.; Minaev, B. Response theory and calculations of spin-orbit coupling phenomena in molecules. Adv. Quantum Chem. 1996, 27, 71–162. [Google Scholar]
  33. Koseki, S.; Schmidt, M.W.; Gordon, M.S. Effective nuclear charges for the first- through third-row transition metal elements in spin–orbit calculations. J. Phys. Chem. A 1998, 102, 10430–10436. [Google Scholar] [CrossRef]
  34. De Simone, B.C.; Mazzone, G.; Russo, N.; Sicilia, E.; Toscano, M. Computational investigation of the influence of halogen atoms on the photophysical properties of tetraphenylporphyrin and its Zinc(II) complexes. J. Phys. Chem. A 2018, 122, 2809–2815. [Google Scholar] [CrossRef]
  35. Alberto, M.E.; De Simone, B.C.; Mazzone, G.; Marino, T.; Russo, N. Photophysical properties of free and metallated mesosubstituted tetrabenzotriazaporphyrin from density functional theory investigation. Dyes Pigm. 2015, 120, 335–339. [Google Scholar] [CrossRef]
  36. De Simone, B.C.; Mazzone, G.; Sang-aroon, W.; Marino, T.; Russo, N.; Sicilia, E. Theoretical insight into joint photodynamic action of a gold(I) complex and a BODIPY chromophore for singlet oxygen generation. Phys. Chem. Chem. Phys. 2019, 21, 3446–3452. [Google Scholar] [CrossRef]
  37. Alberto, M.E.; De Simone, B.C.; Mazzone, G.; Sicilia, E.; Russo, N. The heavy atom effect on Zn(II) phthalocyanine derivatives: A theoretical exploration of the photophysical properties. Phys. Chem. Chem. Phys. 2015, 17, 23595–23601. [Google Scholar] [CrossRef]
  38. Alberto, M.E.; Iuga, C.; Quartarolo, A.D.; Russo, N. Bisanthracene Bis(dicarboxylic imide)s as potential photosensitizers in Photodynamic therapy. A theoretical investigation. J. Chem. Inf. Model. 2013, 53, 2334–2340. [Google Scholar] [CrossRef] [PubMed]
  39. Alberto, M.E.; De Simone, B.C.; Russo, N.; Sicilia, E.; Toscano, M. Can BODIPY Dimers Act as Photosensitizers in Photodynamic Therapy? A Theoretical Prediction. Front. Phys. 2018, 6, 143. [Google Scholar] [CrossRef]
Scheme 1. Schematic representation of compounds investigated in this work.
Scheme 1. Schematic representation of compounds investigated in this work.
Ijms 20 02002 sch001
Figure 1. Molecular orbital composition for the excited states involved in the singlet and triplet excitation for molecules for 1 and 3 molecules.
Figure 1. Molecular orbital composition for the excited states involved in the singlet and triplet excitation for molecules for 1 and 3 molecules.
Ijms 20 02002 g001
Table 1. B3LYP/6-31+G* vertical excitation energies, ΔE in eV and (nm), oscillator strengths f and main orbital transitions for the studied compounds (in toluene). H = HOMO; L = LUMO. Experimental values are taken from ref. [12].
Table 1. B3LYP/6-31+G* vertical excitation energies, ΔE in eV and (nm), oscillator strengths f and main orbital transitions for the studied compounds (in toluene). H = HOMO; L = LUMO. Experimental values are taken from ref. [12].
MoleculeStateΔE fTransitionsλexp
1S12.07(599)0.247H→L, 87%690
T11.31(948) H→L, 93%
T21.36(910) H−1→L, 84%
T32.01(617) H→L+1, 87%
2S11.99(623)0.244H→L, 80%
T11.22(1013) H→L, 81%
T21.38(898) H-1→L, 59%
T31.67(741) H→L+1, 75%
T41.82(682) H-2→L, 57%
3S11.83(678)0.000H-1→L, 77%
S21.95(636)0.279H→L, 79%
T11.15(1076) H→L, 85%
T21.38(898) H-2→L, 55%
T31.54(806) H-1→L, 70%
T41.62(764) H→L+1, 70%
4S1206(602)0.200H→L, 82%680
T11.29(963) H→L, 66%
T21.44(863) H-1→L, 58%
T31.93(641) H→L+1, 84%
5S12.02(613)0.333H→L, 89%703
T11.29(964) H→L, 98%
T21.33(929) H-1→L, 87%
T32.00(620) H→L+1, 87%
6S12.02(614)0.362H→L, 89%
T11.29(964) H→L, 100%
T21.33(933) H-1→L, 90%
T32.00(619) H→L+1, 87%
7S11.99(624)0.457H→L, 86%709
T11.26(988) H→L, 88%
T21.38(896) H-1→L, 80%
T31.98(627) H→L+1, 86%
8S11.83(679)0.000H-2→L, 60%
S21.89(655)0.373H→L, 82%
T11.12(1109) H→L, 87%
T21.40(883) H-1→L, 47%
T31.54(806) H-2→L, 53%
T41.59(781) H→L+1, 72%
9S12.01(616)0.355H→L, 84%695
T11.23(1007) H→L, 66%
T21.45(858) H-1→L, 63%
T31.95(635) H→L+1, 85%
Table 2. Spin–orbit matrix elements (cm−1) and singlet-triplet energy gaps (eV) calculated at the B3LYP/6-31G(d) level of theory.
Table 2. Spin–orbit matrix elements (cm−1) and singlet-triplet energy gaps (eV) calculated at the B3LYP/6-31G(d) level of theory.
Cmpd | Ψ S 1 | H ^ s o | Ψ T 1 | ΔES1-T1 | Ψ S 1 | H ^ s o | Ψ T 2 | ΔES1-T2 | Ψ S 1 | H ^ s o | Ψ T 3 | ΔES1-T3 | Ψ S 1 | H ^ s o | Ψ T 4 | ΔES1-T4
12.60.81.90.70.40.1
2114.70.82.20.60.60.352.20.2
3685.80.78.40.51.50.3239.50.2
42.10.83.20.61.00.1
52.80.732.10.690.90.0
617.40.771.10.6935.20.0
70.20.80.20.610.10.0
8675.40.79.20.433.90.3258.60.2
90.10.80.10.560.20.1

Share and Cite

MDPI and ACS Style

Alberto, M.E.; De Simone, B.C.; Sicilia, E.; Toscano, M.; Russo, N. Rational Design of Modified Oxobacteriochlorins as Potential Photodynamic Therapy Photosensitizers. Int. J. Mol. Sci. 2019, 20, 2002. https://doi.org/10.3390/ijms20082002

AMA Style

Alberto ME, De Simone BC, Sicilia E, Toscano M, Russo N. Rational Design of Modified Oxobacteriochlorins as Potential Photodynamic Therapy Photosensitizers. International Journal of Molecular Sciences. 2019; 20(8):2002. https://doi.org/10.3390/ijms20082002

Chicago/Turabian Style

Alberto, Marta Erminia, Bruna Clara De Simone, Emilia Sicilia, Marirosa Toscano, and Nino Russo. 2019. "Rational Design of Modified Oxobacteriochlorins as Potential Photodynamic Therapy Photosensitizers" International Journal of Molecular Sciences 20, no. 8: 2002. https://doi.org/10.3390/ijms20082002

APA Style

Alberto, M. E., De Simone, B. C., Sicilia, E., Toscano, M., & Russo, N. (2019). Rational Design of Modified Oxobacteriochlorins as Potential Photodynamic Therapy Photosensitizers. International Journal of Molecular Sciences, 20(8), 2002. https://doi.org/10.3390/ijms20082002

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop