
==== Front
RSC Adv
RSC Adv
RA
RSCACL
RSC Advances
2046-2069
The Royal Society of Chemistry

d4ra05585c
10.1039/d4ra05585c
Chemistry
Porphyrin photosensitizer molecules as effective medicine candidates for photodynamic therapy: electronic structure information aided design†
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra05585c

Yin Wei-Huang a
Li Peng-Yuan b
Huang Hou-Hou b
Feng Lu b
Liu Shu-Hui b
Liu Xin a
https://orcid.org/0000-0001-9398-1407
Bai Fu-Quan b
a Department of Stomatology, China-Japan Union Hospital of Jilin University Changchun Jilin 130033 P. R. China liuxin99@jlu.edu.cn

b Institute of Theoretical Chemistry and College of Chemistry, Jilin University Changchun 130023 P. R. China baifq@jlu.edu.cn

16 9 2024
12 9 2024
16 9 2024
14 40 2936829383
1 8 2024
10 9 2024
This journal is © The Royal Society of Chemistry
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by-nc/3.0/ This article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given and it is not used for commercial purposes.
Traditional photosensitizers (PS) in photodynamic therapy (PDT) have restricted tissue penetrability of light and a lack of selectivity for tumor cells, which diminishes the efficiency of PDT. Our aim is to effectively screen porphyrin-based PS medication through computational simulations of large-scale design and screening of PDT candidates via a precise description of the state of the light-stimulated PS molecule. Perylene-diimide (PDI) shows an absorption band in the near-infrared region (NIR) and a great photostability. Meanwhile, the insertion of metal can enhance tumor targeting. Therefore, on the basis of the original porphyrin PS segments, a series of metalloporphyrin combined with PDI and additional allosteric Zn-porphyrin-PDI systems were designed and investigated. Geometrical structures, frontier molecular orbitals, ultraviolet-visible (UV-vis) absorption spectra, adiabatic electron affinities (AEA), especially the triplet excited states and spin–orbit coupling matrix elements (SOCME) of these expanded D–A porphyrin were studied in detail using the density functional theory (DFT) and time-dependent density functional theory (TDDFT) methods. PS candidates, conforming type I or II mechanism for PDT, have been researched carefully by molecular docking which targeted Factor-related apoptosis (Fas)/Fas ligand (Fasl) mediated signaling pathway. It was found that porphyrin-PDI, Fe2-porphyrin-PDI, Zn-porphyrin-PDI, Mg-porphyrin-PDI, Zn-porphyrin combined with PDI through single bond (compound 1), and two acetylenic bonds (compound 2) in this work would be proposed as potential PS candidates for PDT process. This study was expected to provide PS candidates for the development of novel medicines in PDT.

Novel porphyrin photosensitizer molecules, incorporating perylene diimide groups (PDI), were meticulously designed and screened using electronic structure calculations.

Department of Science and Technology of Jilin Province 10.13039/501100011789 3D5203974430 & 3D5198153430 South China University of Technology 10.13039/501100005015 2022-skllmd-09 pubstatusPaginated Article
==== Body
pmc1 Introduction

Cancer remains one of the leading causes of death worldwide, despite the notable efficacy of traditional anticancer modalities, such as surgery and chemoradiotherapy.1 However, photodynamic therapy (PDT), an alternative to contemporary cancer treatment, has emerged as a successful and clinically approved therapeutic modality.2 PDT provides several therapeutic preponderances, including selectivity, non-invasiveness, high cure rates and low side effects for the neoplastic management and non-malignant diseases.3

PDT is based on the generation of reactive oxygen species (ROS) in order to bring about the induction of death in cancerous tissues.4 The PDT process involves photosensitizers (PS), molecular oxygen, and light source.5 Mechanistically, the effect occurs when appropriate energy is absorbed optically by PS, transiting from ground state (S0) to the first singlet excited state (S1). After that, the excited triplet state (T1) is formed through Inter System Crossing (ISC). Subsequently, the T1 of PS participates in the type I electron transfer reactions to produce ROS. In the type II pathway, the PS transfers energy to triplet molecular oxygen (3O2), generating highly singlet oxygen (1O2).6,7 The production of molecular oxygen, as mentioned above in tissues and cells, initiates three main cell death models: apoptotic, necrotic, and autophagic.8–10 Significantly, Fas/Fasl mediated signaling pathway is one of the most important apoptosis-related signaling pathways.11 The abnormal expression of Fas/Fasl was found in the majority cancers, including head and neck squamous cell carcinoma.12

At present, PDT is mainly exploited in the treatment of superficial diseases, such as premalignant conditions, carcinoma in situ or superficial tumors.13 However, traditional PS in PDT has a limited excitation wavelength at the phototherapy window (620–850 nm), and other wavelengths of visible light can also damage healthy tissues.14 Meanwhile, the production of ROS is slashed by the hypoxic environment in the tumor.15 Therefore, improving the photoelectric conversion efficiency of the phototherapy window (focusing on the milder and safer wavelength of 700–850 nm) and enhancing the capacity of oxygen excitation and production are the critical measures to heighten the efficacy of PDT.

The porphyrin segment is the most prominent variety of PS in clinical trials, which is widely present in hemoglobin (iron porphyrin), hemocyanin (copper porphyrin) in animals, chlorophyll (magnesium porphyrin), and cobalt porphyrin in plants.16 Porphyrin compounds are a class of organic macrocyclic molecules with highly π-conjugated structure, with rich structural diversity and excellent photoelectric properties.17 The structure of the porphyrin molecule enables straightforward modification by substituting hydrogen atoms at the 8 β positions and the 4 meso-sites with various metal ions, substituents, or conjugated systems.18 This tunability not only facilitates synthesis and modification but also enables the optimization of their electronic and photophysical properties. As a result, various porphyrin derivatives can be produced, enhancing photosensitizing capabilities while maintaining lower toxicity and improved controllability. Paradoxically, the inadequate selectivity to neoplastic tissues, and weak absorption in the near-infrared region (NIR) of porphyrin compounds limited their further applications in clinical settings.19 Studies have shown that metalloporphyrin has a good affinity for neoplastic tissues and could selectively stay in cancer cells.20,21 Additionally, the organic chromophores were introduced into metal complexes can improve the triplet excited state lifetime.22 Perylene diimide (PDI), a super planar-conjugated system with a special energy level structure, exhibits excellent photochemical properties and tissue penetration in NIR compared to visible light.23–26 PDI is highly preferred in the field of organic solar cell research, due to its ability to be conjugated with porphyrins, allowing for adjustments in absorption, energy levels, and photovoltaic performance.27,28 This suggests that the porphyrin-PDI conjugate, forming the Donor–Acceptor (D–A) structure, holds significant promise and has the potential to enhance absorption in the near-infrared (NIR) region.29

Herein, a series of porphyrin-based derivatives in group 1 without or with metals (Fe, Zn and Mg) were designed (Scheme 1a), and the D–A metalloporphyrin complexes in group 2 were then constructed with PDI using a push–pull strategy (Scheme 1b). Density functional theory (DFT) is the most effective method to study the electronic structure of their excited states.30 During the triplet state calculation, it needs to be emphasized that among these design PS, multi-spin states are considered for Fe-porphyrin and Fe-porphyrin-PDI.31 Therefore, all calculation contents were performed in singlet and quintet states of Fe to discuss in detail.

Scheme 1 Molecular structures of (a) M-porphyrin and (b) M-porphyrin-PDI. Fe1: singlet state; Fe2: quintet state.

In brief, the aim of this study is to obtain structural information and photophysical property prediction of all designed porphyrin derivatives more quickly through theoretical simulation. Different from molecular dynamics methods for studying some medicine candidates, PDT requires a more rigorous description of its excited state electronic structure, so that its photoexcited pharmacology can be elaborated. In addition, molecular docking analysis of the PS candidates against cancer-related targeted proteins was also performed. Ultimately, effective porphyrin-based PS were screened as potential candidates for PDT. This study provides valuable information and theoretical foundation for the design and screening of porphyrin-based PS medicine.

2 Method

2.1 DFT calculations

The geometries of the ground state and excited states of all systems were optimized using the density functional theory (DFT) and time-dependent density functional theory (TDDFT) methods respectively. After literature analysis and functional testing, The M06 exchange and correlation functional was selected in this work for electronic structure and electron excited transition calculation.32,33 The 6-31G (d) basis set was used for C, N, H and O atoms.34 The LanL2DZ basis set was used for Fe, Zn and Mg.35 The implicit solvent model PCM was used to simulate the water solvent environment.36 All geometries reached the minimum point on the potential energy surfaces, which was confirmed through frequency analysis. All computations were performed without symmetry constraints. The above quantum chemical calculations were performed using with Gaussian16 program package.37

Spin–orbit coupling (SOC) was calculated through the m06 and ωb97xd method (Fe1-porphyrin and Fe2-porphyrin),38 which was coupled with the def2-TZVP basis set by Orca with TD-DFT.39 The spin–orbit couplings are defined according to the formula:1

where ĤSO represents the spin–orbit Hamiltonian. ΨSi and ΨTj stand for the wave function of ground state and triplet state, respectively.

2.2 Molecular docking

AutoDock Tools (ADT) was used for docking study.40 The three-dimensional (3D) structures of proteins were obtained from the RCSB Protein Data Bank (PDB) such as Fas [PDB ID: 1DDF] and Fasl [PDB ID: 5L19]. PDBQT files for protein and ligands preparation and grid box creation, were completed using the graphical user interface program Autodock 4.2.6. Visualization of docking results was performed by using Pymol.41 The binding free energy of complex (ΔGbind) calculated by ADT was the difference between the free energies of the complex (ΔGcomplex) and the unbound receptor (ΔGrec) and the free ligand (ΔGlig) as described by the following equations:2 ΔGbind = ΔGcomplex − ΔGreceptor − ΔGligand

3 Result and discussion

PDT has been considered as a clinical option for the treatment of solid tumors. It relies on the ingested PS, following a certain time interval, which is activated by light with a specific wavelength.42 The mechanism of PDT enforces the following key requirements of effective PS: (1) the capacity to absorb light to the first singlet excited state in the therapeutic window while oscillator strength lager than 0.01; (2) the T1–S0 energy gap larger than 0.98 eV, which is the amount of energy required to activate molecular oxygen; (3) the electron affinity of molecular oxygen should be higher than that of PS for obtaining the electron from PS; (4) the SOC values higher than 0.24 cm−1 that ensure an efficient ISC process.43–45

3.1 Molecular geometries and electronic structures

The geometrically optimized structures of the porphyrin derivatives were shown in Fig. 1a and b. All molecules consist of a porphyrin structure in the center, surrounded by four benzene rings, except for PDI. The porphyrin in the center forms a large conjugated system with a planar structure, while the four surrounding benzene rings maintain a certain dihedral angles to the porphyrin structure, and the values of these dihedral angles were listed in Table 1. A diagram of Fe1 was presented in Fig. 1c.

Fig. 1 Optimized ground-state structure of (a) M-porphyrin and (b) M-porphyrin-PDI. (c) The rotation of four benzene rings to porphyrin and the length of metal atoms to N atoms. Fe1: singlet state; Fe2: quintet state.

The dihedral angles (deg) between the porphyrin and the four benzene rings (A, B, C and D) and the bond lengths (Å) of the different metal atoms to the four N atoms (N1, N2, N3 and N4)

	Porphyrin	Fe1-porphyrin	Fe2-porphyrin	Zn-porphyrin	Mg-porphyrin	
Dihedral angles	
A-porphyrin	62.293	65.011	62.947	62.840	62.146	
B-porphyrin	62.618	64.999	62.363	62.701	61.965	
C-porphyrin	62.273	64.829	62.944	62.941	62.125	
D-porphyrin	62.625	64.998	62.365	62.843	61.965	
 	
Bond lengths	
N1–M	—	1.996	2.071	2.065	2.064	
N2–M	—	1.996	2.048	2.065	2.064	
N3–M	—	1.996	2.071	2.065	2.064	
N4–M	—	1.996	2.048	2.065	2.064	

For the porphyrin molecule, the values of the four dihedral angles were 62.29°, 62.62°, 62.27°, and 62.62°. Compared to porphyrin, the dihedral angle between the benzene ring and the porphyrin was slightly increased in Fe1-porphyrin, while it remains essentially unchanged in Fe2-porphyrin, Zn-porphyrin and Mg-porphyrin. Table 1 also lists the bond lengths of the different metal atoms with the four N atoms. In Zn-porphyrin, the length of the four N–Zn bond was essentially equal, all of which were 2.065 Å. The same was true for the N–Mg bond length of Mg-porphyrin and the N–Fe bond length of Fe1-porphyrin, which are 2.064 Å and 1.996 Å, respectively. The N–Fe bond length was slightly shorter than N–Zn and N–Mg. In Fe2-porphyrin, the two N–Fe bond lengths of 2.071 Å for B1 and B3 and 2.048 Å for B2 and B4 were significantly different from the other three molecules.

The obtained vertical excitation energies, excitation wavelengths, oscillator strengths, and the contributions of orbital transitions to electron excitation were reported in Table 2. Frontier molecular orbitals for all complexes were shown in Fig. 2. Frontier molecular orbitals are important for describing the stability using the LUMO–HOMO gap parameter.46 Among the porphyrin derivatives, the HOMO, LUMO, and LUMO+1 of all molecules was delocalized on the porphyrin structure, except for the LUMO orbital of Fe1-porphyrin, which was centrally distributed on the Fe atom. In porphyrin-PDI, Zn-porphyrin-PDI, and Mg-porphyrin-PDI complexes, the HOMO was delocalized on the porphyrin structure while the LUMO was delocalized on the PDI structure. In Fe1-porphyrin-PDI and Fe2-porphyrin-PDI, the HOMO is delocalized over the porphyrin segment. Notably, the LUMO in Fe2-porphyrin-PDI is localized exclusively on the PDI, whereas in Fe1-porphyrin-PDI, it is distributed over both the porphyrin and the PDI. The stronger orbital coupling between Fe1-porphyrin and PDI facilitates LUMO delocalization across both regions. In contrast, the weaker coupling in Fe2-porphyrin-PDI results in the LUMO being primarily confined to the PDI. This difference highlights the impact of spin states on the orbital energy and distribution within the system. Additionally, after incorporating PDI, we observed that the composition of the HOMO−1 remained largely unchanged compared to the standalone metal-porphyrin. However, the LUMO+1 orbital showed greater delocalization across both Fe1-porphyrin and PDI, indicating stronger interactions between Fe1 and the ligands. This suggests that the inclusion of PDI enhances electronic transfer between Fe1 and the surrounding ligands, potentially increasing the efficiency of ISC. Therefore, the PS could realize intramolecular electron transfer upon excited, further demonstrating its D–A structural feature. Since a smaller gap value implied electron easily to transit, our results (Fig. 3) suggested that the gap value of group 2 were smaller than group 1 (0.80, 0.74, 0.86, 0.93, 0.94 eV). It implied that the addition of PDI could reduce the HOMO–LUMO gap.

The essential vertical excitations of M-porphyrin and M-porphyrin-PDI

	Excitation energy (eV)	Wavelength (nm)	f	Main configuration	
Porphyrin	2.1715	570.97	0.0134	H → L (62%)	
Fe1-porphyrin	2.3314	531.81	0.0224	H → L+1 (57%)	
Fe2-porphyrin	2.0054	618.24	0.0106	H → L+1 (71%)	
Zn-porphyrin	2.1945	564.98	0.0501	H → L (64%)	
Mg-porphyrin	2.1735	570.43	0.0540	H → L (64%)	
PDI	2.3958	517.50	0.7816	H → L (100%)	
Porphyrin-PDI	1.6992	729.67	0.4667	H → L (96%)	
Fe1-porphyrin-PDI	1.6944	731.74	0.3899	H → L (98%)	
Fe2-porphyrin-PDI	1.1181	742.33	0.3941	H → L (77%)	
Zn-porphyrin-PDI	1.6134	768.47	0.3495	H → L (98%)	
Mg-porphyrin-PDI	1.5854	782.05	0.3578	H → L (98%)	

Fig. 2 Graphical representation of the molecular orbitals of M-porphyrin and M-porphyrin-PDI.

Fig. 3 Calculated frontier molecular orbital energies of ground states of (a) M-porphyrin and (b) M-porphyrin-PDI. (c) LUMO–HOMO energy gap.

3.2 Absorption spectra

The wavelength chosen for PDT is typically coincides with the longest-wavelength absorption band of the PS, ideally occurring between 620 and 850 nm.15 The classic UV-vis absorption spectra of porphyrin compounds can be divided into two characteristic absorption bands: a strong Soret band (350–450 nm) and several weak Q bands (500–650 nm).20 If some groups were introduced into the porphyrin substrates, the absorption wavelength and the number of peaks in the Q band would be changed.16 This was due to the symmetry of the porphyrins, variations in the conjugated system, and the degree of distortion.47 In our study on absorption spectra, group 1 (Fig. 4a) showed typical absorption characteristic of metalloporphyrin in the UV-vis region. However, the addition of metal did not produce a noticeable red-shift effect on the results. In Table 2, their maximum absorption wavelengths were mainly approximately 400 nm, and the absorption within the range of 500 to 650 nm was extremely weak, with an oscillator strength of less than 0.1. In the range of 500–650 nm, for all porphyrin derivatives, the largest contribution to the electronic excitation of all porphyrin derivatives was the HOMO → LUMO transition, except for Fe1/Fe2-porphyrin, which was a HOMO → LUMO+1 transition, and the electronic excitations of all the molecules were localized excitations on the porphyrin structure. In contrast, in porphyrin-PDI complexes, the HOMO–LUMO transitions involve charge transfer processes from porphyrin to PDI or from PDI to porphyrin, resulting in an appearance of the new absorption band at 700–800 nm and a large oscillator strength. The value of group 2 was sufficient in phototherapy window (Fig. 4b). With the insertion of PDI, it generated a “push–pull” structure with a bigger conjugated system. New peaks appeared in the spectra due to the extended charge transfer transition state, which were located in the NIR. The significant advantage of these compounds lies in their capability to absorb at substantially longer wavelengths (700–850 nm) compared to commonly used porphyrins such as Foscan® (652 nm). Group 2 were chosen to discuss in the following study. Their absorption wavelength occurred in 700–850 nm which would have a deeper tissue penetrability.

Fig. 4 Simulated absorption spectra of (a) M-porphyrin and (b) M-porphyrin-PDI.

3.3 Electron affinity

Electron Affinity (EA) included the vertical electron affinities (VEA) and adiabatic electron affinities (AEA). Considering the time of charge relaxation and charge transfer between two molecules, AEA was chosen to calculate. The region around the Zn atom was analysed, as it is the most probable site for electron loss, with a calculated value of 47.37 kcal mol−1 (2.05 eV) during the adsorption of 3O2 at a distance of 2.38 Å. Before adsorption, the overall average energies of Zn-porphyrin-PDI and 3O2 were 0.20 kcal mol−1 (0.009 eV) and 1.02 kcal mol−1 (0.044 eV), respectively. After adsorption, these values shifted to 0.16 kcal mol−1 (0.007 eV) and −31.44 kcal mol−1 (−1.36 eV), respectively. These results are used to verify our point of view. In order to verify whether the designed PS molecules could pass electrons to 3O2 to form O2− (type I), AEA values between 3O2 and PS molecules were compared (Fig. 5). An ideal PS' AEA value should be less than that of 3O2 to allow for the efficient transfer of electrons from the PS to oxygen, facilitating the reduction of oxygen.3 AEA (3O2) − AEA (PS) > 0

Fig. 5 Adiabatic electron affinities for (a) M-porphyrin and (b) M-porphyrin-PDI. The energy line of 3.73 eV for molecular oxygen was highlighted as the red dashed line.

The calculated AEA of 3O2 was 3.73 eV, hence porphyrin-PDI, Fe2-porphyrin-PDI, Zn-porphyrin-PDI and Mg-porphyrin-PDI were screened as the possible candidates for PS that followed type I PDT mechanism.

3.4 Triplet excited states

To induce irreversible destruction of neoplastic tissues, the prevailing belief is that type II mechanisms should predominate in PDT, owing to their relatively higher reactivity.48 To investigate the ability of the designed PS to excite 3O2 to 1O2, we compared the S0–T1 energy gap with the energy required to excite oxygen from 3O2 to 1O2.4 PS (T1) + 3O2 → PS (S0) + 1O2

The energy needed for the formation of singlet oxygen is 0.98 eV.44 The calculated data suggested that the value of S0–T1 energy gap was greater than 0.98 eV including porphyrin-PDI, Fe2-porphyrin-PDI (Fig. 6). It was considered that these PS molecules designed here meet this energy requirement except the Fe1-porphyrin and Fe1-porphyrin-PDI, mainly due to their non-singlet ground state and possible multi spin state characteristics. The calculated S1–T1 energy gaps (ΔEST), provided in the ESI,† indicate that all compounds possess sufficient triplet state lifetimes, increasing their chances of reacting with oxygen to generate ROS. Except for Fe1-porphyrin, the incorporation of PDI led to a slight reduction in the ΔEST across the compounds, suggesting that the presence of PDI enhances ISC efficiency by reducing the ΔEST.

Fig. 6 Triplet excitation energies for (a) M-porphyrin and (b) M-porphyrin-PDI. The energy line of 0.98 eV, required to activate molecular oxygen, was highlighted as the black dashed line.

3.5 Spin–orbit coupling matrix elements

In the process of ISC, the excited PS undergoes a transition from a singlet state to a triplet state. This triplet PS can then transfer its excitation energy to 3O2, facilitating the generation of ROS. The efficiency of the ISC hinges on the amplitudes of SOC for the S1 → Tn which effects the production of 3O2.49

The value of SOC should be greater than 0.24 cm−1. The SOC of the S1–Tn (T1–T5) was calculated to discuss whether porphyrin complexes could have the effective coupling channel from S1 to Tn (Fig. 7). Group 2 were met the requirements with 3–5 effective ISC channels. Comparing the SOC values for all the derivatives, it was observed that the insertion of PDI contributed to the enhancement of SOC's values.

Fig. 7 Spin–orbit coupling matrix elements between low-lying singlet and triplet excited states for (a) M-porphyrin and (b) M-porphyrin-PDI.

A deficiency was found through the analysis of the above data. In the calculation of triplet excited states, Zn-porphyrin-PDI (0.94 eV) and Mg-porphyrin-PDI (0.97 eV) were applicable in the type I of PDT but did not playing a role in the type II of PDT. The structure of Zn-porphyrin-PDI, performing less in type II mechanism, was selected for improvement. Based on the experience accumulated from the above work, the molecular configurations of Zn-porphyrin-PDI were expanded in group 3 by changing conjugate length and connective positions between Zn-porphyrin and PDI, or increasing the number of PDI on the basis of the Zn-porphyrin. We aimed to investigate how these structural modifications could enhance the PS's ability to generate ROS through both mechanisms. This approach was designed to better identify photosensitizers that can fulfill the requirements of both Type I and Type II mechanisms in PDT. Zn-porphyrin combined with PDI through single bond (compound 1) or two acetylenic bonds (compound 2). One PDI was added to Zn-porphyrin-PDI (compound 3). Zn-porphyrin combined with the imide position of PDI by one acetylenic bond (compound 4) or two acetylenic bonds (compound 5) (Scheme 2). It was further screened based on previous calculation contents in group 3.

Scheme 2 Molecular structures of compound 1–5.

3.6 The calculations of Zn-porphyrin-PDI derivatives

The geometrically optimized structures of group 3 were shown in Fig. 8a. The frontier molecular orbitals were shown in Fig. 8b. The LUMO and HOMO were mainly located on the PDI derivatives and the porphyrin like Zn-porphyrin-PDI, respectively. The photophysical properties together with their orbital contribution were reported in Table 3. The results of the absorption spectra, frontier molecular orbitals, AEA, triplet excited states and SOC were demonstrated in Fig. 9. In the calculation of absorption spectra, the wavelengths of compound 1–3 were in phototherapy window (700–850 nm). Compared with Zn-porphyrin-PDI, decreasing conjugate length (compound 1) would make absorption spectrum red-shift to 780.17 nm, and LUMO–HOMO gap and the triplet excited states energy grow to 2.14 eV and 1.34 eV, while AEA value decreased to 3.63 eV. Expending conjugate length (compound 2) would blue-shift the wavelength to 750.48 nm, reducing LUMO–HOMO gap to 2.06 eV, while the AEA value and triplet excited states energy were increased to 3.84 eV and 1.31 eV, respectively. Adding an extra PDI to Zn-porphyrin-PDI (compound 3) would red-shifted the absorption to 771.62 nm, shorten gap value and triplet excited states energy to 1.47 eV and 0.13 eV, and rise AEA value to 4.39 eV. After combining Zn-porphyrin with PDI's imide position (compound 4), absorption peak blue-shifted to 606.40 nm. The gap value decreased to 1.95 eV. AEA value and triplet excited states energy up to 3.77 eV and 1.26 eV. Compared with compound 4, extending conjugate length on the basis of the compound 4 (compound 5) would red-shift the absorption peak value to 613.78 nm. Gap value, AEA value and the triplet excited states energy would increase to 2.02 eV, 3.78 eV and 1.36 eV. The ΔEST of M-porphyrin and M-porphyrin-PDI are provided in Table S2 (ESI†). Compared to Zn-porphyrin-PDI, increasing the number of PDI units does not significantly alter the ΔEST. However, modifying the conjugation length between the porphyrin and PDI or changing the substitution position on PDI can notably reduce the ΔEST, thereby enhancing ISC efficiency.

Fig. 8 (a) Optimized ground-state structures of compounds 1–5 and (b) illustrations of the molecular orbitals of compounds 1–5.

The essential vertical excitations of compound 1–5

	Excitation energy (eV)	Wavelength (nm)	f	Main configuration	
Compound 1	1.5892	780.17	0.0807	H → L (99%)	
Compound 2	1.6521	750.48	0.3624	H → L (97%)	
Compound 3	1.6068	771.62	0.4452	H → L+1 (84%)	
Compound 4	2.0446	606.40	0.6136	H → L+1 (78%)	
Compound 5	2.0200	613.78	0.6529	H → L+1 (80%)	

Fig. 9 (a) Calculated absorption spectra, (b) frontier molecular orbital energies, (c) adiabatic electron affinities, (d) triplet excitation energies and (e) spin–orbit coupling matrix elements of compound 1–5.

The results showed that the shorter the conjugate distance, the longer absorption wavelength, the lower AEA value would be. By combining with imide position of PDI could narrow LUMO–HOMO gap and increase the triplet excited states energy. However, gap value was decreased by adding an PDI. These allosteric strategies may also be applied to Mg-porphyrin derivatives based on the similar electronic structures and coordination condition of Mg-porphyrin and Zn-porphyrin derivatives. Therefore, compound 1 and 2 were screened as the possible candidates in group 3 that followed type I/II and type II PDT mechanism.

3.7 Molecular docking

Based on the above calculation results, we selected PS candidates applying the type I (porphyrin-PDI, Fe2-porphyrin-PDI, Zn-porphyrin-PDI, Mg-porphyrin-PDI and compound 1) and type II (porphyrin-PDI, Fe2-porphyrin-PDI, compound 1 and 2) mechanism for molecular docking study. The molecular docking is an important method to clarify the interaction between a ligand and a target protein, which has been widely used for the structure-based drug design in pharmaceutical research.50,51 In our study, the screened PS candidates, porphyrin-PDI, Fe2-porphyrin-PDI, Zn-porphyrin-PDI, Mg-porphyrin-PDI, compound 1 and compound 2 were used as ligands. Fas and Fasl, apoptosis associated proteins, was one of the most important signal pathways of PDT research. At present, it has been verified that Aminolevulinate and chlorin e6 (they all belong to porphyrin precursors or porphyrins) can induce significantly greater Fas/Fasl which could further induce tumor cell apoptosis.52,53 Thus, they were used as targeted proteins which were docked with the six selected PS candidates. The estimated free energy of binding (EFEB) (kcal mol−1) values of the protein-ligand complexes were calculated. The EFEB of candidates which targeted Fas and Fasl were presented in Table 4.

Free energy of binding of selected compounds targeting Fas and Fasl

Ligand	Binding energy (kcal mol−1)	
Fas	Fasl	
Porphyrin-PDI	−9.8	−11.0	
Fe2-porphyrin-PDI	−10.2	−11.6	
Zn-porphyrin-PDI	−10.4	−11.6	
Mg-porphyrin-PDI	−10.3	−11.7	
Compound 1	−10.0	−9.7	
Compound 2	−10.1	−11.7	

The minimum energy docked positions and corresponding binding sites of the protein-ligand complexes were shown in Fig. 10. Meanwhile, the pink dashed lines indicated the hydrogen bonds formation between the ligand and selected targeted proteins. As can be observed from the figure, since the PDI portion possesses flexible and extended structural characteristics, the docking site initially emerges around the PDI and the final interaction is also mostly around the PDI fragment. The corresponding hydrogen bond length values between ligand molecules and the amino acid residues in the targeted proteins were also depicted in Fig. 10.

Fig. 10 Lowest energy docked positions of the PS candidates with cancer associated targeted proteins (a) Fas and (b) Fasl. The amino acid residues in the targeted proteins were shown by orange labels. Hydrogen bond interactions were represented in pink dashed lines and the corresponding hydrogen bond length values were marked.

Previous research has demonstrated that metal ion coordination can modify the electronic structure and reactivity of molecules, leading to changes in energy. Additionally, metal coordination can strengthen interactions with specific regions of proteins.54 Our results suggest that the incorporation of metals enhances the affinity between ligands and target proteins. However, the effect is relatively minor because the metal contributes to only one active site within the molecular framework, while the overall binding energy is influenced by the entire molecular structure. However, generally, more negative EFEB value was correlated with stronger binding affinity. A molecule has a good interaction with the proteins studied given by the values of EFEB lower than −6 kcal mol−1.55 All the candidates exhibited better binding capacity (from −9.8 to −11.7 kcal mol−1). These effects are generally within the realm of weak interactions, and their impacts on the spectrum are similar to those of polar solvents. These minor effects do not alter our previous judgment of being good PDT candidates. The obtained results would be significative in vitro and in vivo studies for the development of effective PS in PDT.

Finally, at the same computational level, we conducted calculations for Foscan, which is utilized in the treatment of head and neck squamous cell carcinoma. The calculated values for AEA and T1 were 2.88 eV and 1.56 eV, respectively. The S1–Tn (n = 1–5) values for spin–orbit coupling were 0.05, 0, 0.01, 0, and 0.34 cm−1, respectively. The EFEB of Fas and Fasl were −7.9 kal mol−1 and −7.4 kal mol−1, respectively. In contrast, our screened PS candidates exhibited a distinct advantage in their ability to achieve a stable T1 state and target the Fas/Fasl signaling pathway. Moreover, the porphyrin-based candidates are easily synthesized, with both porphyrin and PDI having verified feasible synthesis pathways.56–58 Although their solubility as materials is not optimal, as small-molecule PS, they are still effective at inducing the PDT process even at extremely very low concentrations.59–61 This is also the advantage of trace medication in PDT. Hence, porphyrinyl PS still possesses broad application prospects.

4 Conclusions

In summary, the conjugated porphyrin-based derivates with D–A structure were designed in this work using porphyrin without or with metals as electron donor and PDI as electron acceptor unit. Group 2 of metalloporphyrin-PDI and compound 1–3 of Zn-porphyrin-PDI PSs displayed intense absorbance in the NIR (700–850 nm) medical scope. Porphyrin-PDI, Fe2-porphyrin-PDI, Zn-porphyrin-PDI, Mg-porphyrin-PDI and compound 1 could transfer electrons to perilesional molecular oxygen to produce ROS via type I mechanism of PDT efficiently. According to type II mechanism, porphyrin-PDI, Fe2-porphyrin-PDI, compound 1 and 2 could excite 3O2 to form 1O2. All above screened candidates could also target Fas/Fasl signal pathway. Hence, the present investigations would be beneficial for the development of PS in PDT. For PDT process, in-depth understanding of the electronic state change behavior of the system induced by light excitation needs to be obtained at the chemical level in order to obtain the most effective key pharmacal efficacy information, so it is an essential research method for the development of such medicines. We believe that label extraction and machine learning based on this is an effective path for future typical pharmacologic development, and our following work is ongoing.

Data availability

All data generated or analysed during this study are included in this article.

Author contributions

Conceptualization, Fu-Quan Bai and Xin Liu; data curation, Wei-Huang Yin, Peng-Yuan Li, Hou-Hou Huang, Lu Feng and Shu-Hui Liu; funding acquisition, Fu-Quan Bai and Xin Liu; validation, Wei-Huang Yin; writing – original draft, Wei-Huang Yin and Peng-Yuan Li; writing – review & editing, Fu-Quan Bai.

Conflicts of interest

There are no conflicts to declare.

Supplementary Material

RA-014-D4RA05585C-s001

The present study was supported by a fund from Department of Science and Technology of Jilin Province, China (Grant No. 3D5203974430 & 3D5198153430) and the ‘14th Five-Year’ Science and Technology Project of the Education Department of Jilin province (Grant No. JJKH20220966KJ) and open fund of the state key laboratory of luminescent materials and devices of South China University of Technology, China (Grant No. 2022-skllmd-09).
==== Refs
Notes and references

Zaimy M. A. Saffarzadeh N. Mohammadi A. Pourghadamyari H. Izadi P. Sarli A. Moghaddam L. K. Paschepari S. R. Azizi H. Torkamandi S. Tavakkoly-Bazzaz J. New methods in the diagnosis of cancer and gene therapy of cancer based on nanoparticles Cancer Gene Therapy Cancer Gene Ther. 2017 24 233 28574057
Correia J. H. Rodrigues J. A. Pimenta S. Dong T. Yang Z. Photodynamic Therapy Review: Principles, Photosensitizers, Applications, and Future Directions Pharmaceutics 2021 13 1332 34575408
Pham T. C. Nguyen V. N. Choi Y. Lee S. Yoon J. Recent Strategies to Develop Innovative Photosensitizers for Enhanced Photodynamic Therapy Chem. Rev. 2021 121 13454 13619 34582186
Larue L. Myrzakhmetov B. Ben-Mihoub A. Moussaron A. Thomas N. Arnoux P. Baros F. Vanderesse R. Acherar S. Frochot C. Fighting Hypoxia to Improve PDT Pharmaceuticals 2019 12 163 31671658
Domka W. Bartusik-Aebisher D. Mytych W. Dynarowicz K. Aebisher D. The Use of Photodynamic Therapy for Head, Neck, and Brain Diseases Int. J. Mol. Sci. 2023 24 11867 37511625
Li X. Lovell J. F. Yoon J. Chen X. Clinical development and potential of photothermal and photodynamic therapies for cancer Nat. Rev. Clin. Oncol. 2020 17 657 674 32699309
Sarbadhikary P. George B. P. Abrahamse H. Recent Advances in Photosensitizers as Multifunctional Theranostic Agents for Imaging-Guided Photodynamic Therapy of Cancer Theranostics 2021 11 9054 9088 34522227
Alzeibak R. Mishchenko T. A. Shilyagina N. Y. Balalaeva I. V. Vedunova M. V. Krysko D. V. Targeting immunogenic cancer cell death by photodynamic therapy: past, present and future J. Immunother. Cancer 2021 9 001926
Rodrigues J. A. Correia J. H. Photodynamic Therapy for Colorectal Cancer: An Update and a Look to the Future Int. J. Mol. Sci. 2023 24 12204 37569580
Arellano L. G. Villar-Alvarez E. M. Velasco B. Domínguez-Arca V. Prieto G. Cambón A. Barbosa S. Taboada P. Light excitation of gold Nanorod-Based hybrid nanoplatforms for simultaneous bimodal phototherapy J. Mol. Liq. 2023 377 121511
Wang J. Cao Z. Wang P. Zhang X. Tang J. He Y. Huang Z. Mao X. Shi S. Kou X. Apoptotic Extracellular Vesicles Ameliorate Multiple Myeloma by Restoring Fas-Mediated Apoptosis ACS Nano 2021 15 14360 14372 34506129
Zhu J. Petit P.-F. Van den Eynde B. J. Apoptosis of tumor-infiltrating T lymphocytes: a new immune checkpoint mechanism Cancer Immunol. Immunother. 2018 68 835 847 30406374
Agostinis P. Berg K. Cengel K. A. Foster T. H. Girotti A. W. Gollnick S. O. Hahn S. M. Hamblin M. R. Juzeniene A. Kessel D. Korbelik M. Moan J. Mroz P. Nowis D. Piette J. Wilson B. C. Golab J. Photodynamic therapy of cancer: an update Ca-Cancer J. Clin. 2011 61 250 281 21617154
Gunaydin G. Gedik M. E. Ayan S. Photodynamic Therapy-Current Limitations and Novel Approaches Front. Chem. 2021 9 691697 34178948
Donohoe C. Senge M. O. Arnaut L. G. Gomes-da-Silva L. C. Cell death in photodynamic therapy: from oxidative stress to anti-tumor immunity Biochim. Biophys. Acta Rev. Canc 2019 1872 188308
Wang X. Lv H. Sun Y. Zu G. Zhang X. Song Y. Zhao F. Wang J. New porphyrin photosensitizers-Synthesis, singlet oxygen yield, photophysical properties and application in PDT Spectrochim. Acta, Part A 2022 279 121447
Jiang X. Zhao Y. Sun S. Xiang Y. Yan J. Wang J. Pei R. Research development of porphyrin-based metal-organic frameworks: targeting modalities and cancer therapeutic applications J. Mater. Chem. B 2023 11 6172 6200 37305964
Dhiman Maitra J. B. C. Elenbaas J. S. Bonkovsky H. L. Shavit J. A. Bishr Omary M. Porphyrin-Induced Protein Oxidation and Aggregation as a Mechanism of Porphyria-Associated Cell Injury Cell. Mol. Gastroenterol. Hepatol. 2019 8 536 548
Tsolekile N. Nelana S. Oluwafemi O. S. Porphyrin as Diagnostic and Therapeutic Agent Molecules 2019 24 2669 31340553
Lan M. Zhao S. Liu W. Lee C. S. Zhang W. Wang P. Photosensitizers for Photodynamic Therapy Adv. Healthc. Mater. 2019 8 e1900132 31067008
Ke J. Zhao X. Yang J. Ke K. Wang Y. Yang M. Yang W. Enhanced Ion-Selective Diffusion Achieved by Supramolecular Interaction for High Thermovoltage and Thermal Stability Energy Environ. Mater. 2023 7 e12562
Wei X. Cui W. b. Qin G.-y. Zhang X. e. Sun F. y. Li H. Guo J. f. Ren A. m. Theoretical Investigation of Ru(II) Complexes with Long Lifetime and a Large Two-Photon Absorption Cross-Section in Photodynamic Therapy J. Med. Chem. 2023 66 4167 4178 36884221
Zhao Z. Xu N. Wang Y. Ling G. Zhang P. Perylene diimide-based treatment and diagnosis of diseases J. Mater. Chem. B 2021 9 8937 8950 34657950
Li H. Yue L. Wu M. Wu F. Self-assembly of methylene violet-conjugated perylene diimide with photodynamic/photothermal properties for DNA photocleavage and cancer treatment Colloids Surf. B Biointerfaces 2020 196 111351 32916437
Lu H. Liu J. Liu Y. Xu X. Bo Z. Improving the Efficiency of Organic Solar Cells by Introducing Perylene Diimide Derivative as Third Component and Individually Dissolving Donor/Acceptor ChemSusChem 2021 14 5442 5449 34581011
Wu X. Wang H. Gong Y. Fan D. Ding P. Li Q. Qian Q. Graph neural networks for molecular and materials representation J. Mater. Inf. 2023 3 12
Pan X. Wu J. Xiao L. Yap B. Xia R. Peng X. Porphyrin Acceptors with Two Perylene Diimide Dimers for Organic Solar Cells ChemSusChem 2021 14 3614 3621 34107177
Yang F. Wu Y. Zhao J. Guo Y. Guo X. Li W. Wang J. Excited-state photophysical processes in a molecular system containing perylene bisimide and zinc porphyrin chromophores Phys. Chem. Chem. Phys. 2020 22 20891 20900 32915174
Li C. Luo Z. Yang L. Chen J. Cheng K. Xue Y. Liu G. Luo X. Wu F. Self-assembled porphyrin polymer nanoparticles with NIR-II emission and highly efficient photothermal performance in cancer therapy Mater. Today. Bio 2022 13 100198
Timothy R. A. Louis H. Adindu E. A. Gber T. E. Agwamba E. C. Offiong O. E. Pembere A. M. S. Elucidation of collagen amino acid interactions with metals (B, Ni) encapsulated graphene/PEDOT material: insight from DFT calculations and MD simulation J. Mol. Liq. 2023 390 122950
Nakashima H. Hasegawa J. Y. Nakatsuji H. On the O2 binding of Fe-porphyrin, Fe-porphycene, and Fe-corrphycene complexes J. Comput. Chem. 2006 27 1363 1372 16788910
Walker M. Harvey A. J. Sen A. Dessent C. E. Performance of M06, M06-2X, and M06-HF density functionals for conformationally flexible anionic clusters: M06 functionals perform better than B3LYP for a model system with dispersion and ionic hydrogen-bonding interactions J. Phys. Chem. A 2013 117 12590 12600 24147965
Thomas Malcomson P. R. c. á. Erhardt S. Paterson M. J. Protocols for Understanding the Redox Behavior of Copper Containing Systems ACS Omega 2022 7 45057 45066 36530299
Hehre W. J. Ditchfield R. Pople J. A. Self-Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian-Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules J. Chem. Phys. 1972 56 2257 2261
Hay P. J. Wadt W. R. Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals J. Chem. Phys. 1985 82 299 310
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 20423165
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. , Li X. , Caricato M. , Marenich A. V. , Bloino J. , Janesko B. G. , Gomperts R. , Mennucci B. , Hratchian H. P. , Ortiz J. V. , Izmaylov A. F. , Sonnenberg J. L. , Ding W. F. , Lipparini F. , Egidi F. , Goings J. , Peng B. , Petrone A. , Henderson T. , Ranasinghe D. , Zakrzewski V. G. , Gao J. , Rega N. , Zheng G. , Liang W. , Hada M. , Ehara M. , Toyota K. , Fukuda R. , Hasegawa J. , Ishida M. , Nakajima T. , Honda Y. , Kitao O. , Nakai H. , Vreven T. , Throssell K. , Montgomery Jr J. A. , Peralta J. E. , Ogliaro F. , Bearpark M. J. , Heyd J. J. , Brothers E. N. , Kudin K. N. , Staroverov V. N. , Keith T. A. , Kobayashi R. , Normand J. , Raghavachari K. , Rendell A. P. , Burant J. C. , Iyengar S. S. , Tomasi J. , Cossi M. , Millam J. M. , Klene M. , Adamo C. , Cammi R. , Ochterski J. W. , Martin R. L. , Morokuma K. , Farkas O. , Foresman J. B. and Fox D. J. , Gaussian 16 (Revision B.01), Wallingford CT, Gaussian Inc., 2016
Chai J. D. Head-Gordon M. Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections Phys. Chem. Chem. Phys. 2008 10 6615 6620 18989472
Neese F. Wennmohs F. Becker U. Riplinger C. The ORCA quantum chemistry program package J. Chem. Phys. 2020 152 224108 32534543
Morris G. M. , Huey R. and Olson A. J. , Using AutoDock for Ligand-Receptor Docking, Current protocols in bioinformatics, 2008, ch. 8, pp. 8–14
The PyMOL Molecular Graphics System, Version 2.0, Schrödinger, LLC
Ulfo L. Costantini P. E. Di Giosia M. Danielli A. Calvaresi M. EGFR-Targeted Photodynamic Therapy Pharmaceutics 2022 14 241 35213974
Soldatova A. V. Kim J. Rizzoli C. Kenney M. E. Rodgers M. A. Rosa A. Ricciardi G. Near-infrared-emitting phthalocyanines. A combined experimental and density functional theory study of the structural, optical, and photophysical properties of Pd(II) and Pt(II) alpha-butoxyphthalocyanines Inorg. Chem. 2011 50 1135 1149 21188985
Drzewiecka-Matuszek A. Rutkowska-Zbik D. Application of TD-DFT Theory to Studying Porphyrinoid-Based Photosensitizers for Photodynamic Therapy: A Review Molecules 2021 26 7176 34885763
Wang X. Bai F. Q. Liu Y. Wang Y. Zhang H. X. Lin Z. A Computational Way To Achieve More Effective Candidates for Photodynamic Therapy J. Chem. Inf. Model. 2017 57 1089 1100 28383275
Noriega L. Castro M. E. Perez-Aguilar J. M. Caballero N. A. Sanchez-Gaytan B. L. Gonzalez-Vergara E. Melendez F. J. Oxidovanadium(V) complexes as promising anticancer photosensitizers J. Inorg. Biochem. 2020 203 110862 31683130
Zu G. P. Wang J. J. Zhang Y. Chen W. B. Shi Y. Z. Guo S. W. Wang X. R. Study on Five Porphyrin-Based Photosensitizers for Singlet Oxygen Generation ChemistrySelect 2019 4 863 867
De Simone B. C. Alberto M. E. Russo N. Toscano M. Photophysical properties of heavy atom containing tetrasulfonyl phthalocyanines as possible photosensitizers in photodynamic therapy J. Comput. Chem. 2021 42 1803 1808 34236090
Spiegel M. Adamo C. Tuning the Photophysical Properties of Ru(II) Photosensitizers for PDT by Protonation and Metallation: A DFT Study J. Phys. Chem. A 2023 127 3625 3635 37040576
Valarmathi T. Premkumar R. Meera M. R. Milton A. Franklin Benial, Spectroscopic, quantum chemical and molecular docking studies on 1-amino-5-chloroanthraquinone: A targeted drug therapy for thyroid cancer Spectrochim. Acta Mol. Biomol. Spectrosc. 2021 255 119659
Ferreira L. G. Dos Santos R. N. Oliva G. Andricopulo A. D. Molecular docking and structure-based drug design strategies Molecules 2015 20 13384 13421 26205061
Salva K. A. Kim Y. H. Rahbar Z. Wood G. S. Epigenetically Enhanced PDT Induces Significantly Higher Levels of Multiple Extrinsic Pathway Apoptotic Factors than Standard PDT, Resulting in Greater Extrinsic and Overall Apoptosis of Cutaneous T-cell Lymphoma Photochem. Photobiol. 2018 94 1058 1065 29675945
Zhao L. Rao X. Zheng R. Huang C. Kong R. Yu X. Cheng H. Li S. Targeting glutamine metabolism with photodynamic immunotherapy for metastatic tumor eradication J. Controlled Release 2023 357 460 471
Abbehausen C. Zinc finger domains as therapeutic targets for metal-based compounds-an update Metallomics 2019 11 15 28 30303505
Udrea A. M. Dinache A. Staicu A. Avram S. Target Prediction of 5,10,15,20-Tetrakis(4'-Sulfonatophenyl)-Porphyrin Using Molecular Docking Pharmaceutics 2022 14 2930
Adler F. R. L. A. D. Finarelli J. D. Goldmacher J. Assour J. Korsakoff L. A simplified synthesis for meso-tetraphenylporphine J. Org. Chem. 1967 32 476
Lindsey J. S. Schreiman I. C. Hsu H. C. Kearney P. C. Marguerettaz A. M. Rothemund and adler-longo reactions revisited synthesis of tetraphenylporphyrins under equilibrium J. Org. Chem. 1987 52 827 836
Nowak-Król A. Würthner F. Progress in the synthesis of perylene bisimide dyes Org. Chem. Front. 2019 6 1272 1318
Yakavets I. Millard M. Zorin V. Lassalle H.-P. Bezdetnaya L. Current state of the nanoscale delivery systems for temoporfin-based photodynamic therapy: advanced delivery strategies J. Controlled Release 2019 304 268 287
Espitia-Almeida F. Diaz-Uribe C. Vallejo W. Gomez-Camargo D. Bohorquez A. R. R. Zarate X. Schott E. Photophysical characterization and in vitro anti-leishmanial effect of 5,10,15,20-tetrakis(4-fluorophenyl) porphyrin and the metal (Zn(II), Sn(IV), Mn(III) and V(IV)) derivatives BioMetals 2022 35 159 171 34993713
Zhang L. Geng Y. Li L. Tong X. Liu S. Liu X. Su Z. Xie Z. Zhu D. Bryce M. R. Rational design of iridium-porphyrin conjugates for novel synergistic photodynamic and photothermal therapy anticancer agents Chem. Sci. 2021 12 5918 5925 35342539
