Insights into the anticancer photodynamic activity of Ir(III) and Ru(II) polypyridyl complexes bearing ß-carboline ligands.

Ir(III) and Ru(II) polypyridyl complexes are promising photosensitizers (PSs) for photodynamic therapy (PDT) due to their outstanding photophysical properties. Herein, one series of cyclometallated Ir(III) complexes and two series of Ru(II) polypyridyl derivatives bearing three different thiazolyl-ß-carboline N^N' ligands have been synthesized, aiming to evaluate the impact of the different metal fragments ([Ir(C^N)2]+ or [Ru(N^N)2]2+) and N^N' ligands on the photophysical and biological properties. All the compounds exhibit remarkable photostability under blue-light irradiation and are emissive (605 < λem < 720 nm), with the Ru(II) derivatives displaying higher photoluminescence quantum yields and longer excited state lifetimes. The Ir PSs display pKa values between 5.9 and 7.9, whereas their Ru counterparts are less acidic (pKa > 9.3). The presence of the deprotonated form in the Ir-PSs favours the generation of reactive oxygen species (ROS) since, according to theoretical calculations, it features a low-lying ligand-centered triplet excited state (T1 = 3LC) with a long lifetime. All compounds have demonstrated anticancer activity. Ir(III) complexes 1-3 exhibit the highest cytotoxicity in dark conditions, comparable to cisplatin. Their activity is notably enhanced by blue-light irradiation, resulting in nanomolar IC50 values and phototoxicity indexes (PIs) between 70 and 201 in different cancer cell lines. The Ir(III) PSs are also activated by green (with PI between 16 and 19.2) and red light in the case of complex 3 (PI = 8.5). Their antitumor efficacy is confirmed by clonogenic assays and using spheroid models. The Ir(III) complexes rapidly enter cells, accumulating in mitochondria and lysosomes. Upon photoactivation, they generate ROS, leading to mitochondrial dysfunction and lysosomal damage and ultimately cell apoptosis. Additionally, they inhibit cancer cell migration, a crucial step in metastasis. In contrast, Ru(II) complex 6 exhibits moderate mitochondrial activity. Overall, Ir(III) complexes 1-3 show potential for selective light-controlled cancer treatment, providing an alternative mechanism to chemotherapy and the ability to inhibit lethal cancer cell dissemination.

XMS-CASPT2//XMS-CASPT2 and XMS-CASPT2//CASSCF at comparison: The impact of dynamic correlation in the excited state optimization of nitronaphthalene.

Here, analytical extended multi-state complete active space second-order perturbation method (XMS-CASPT2) gradients are used to rationalize the decreasing triplet quantum yield trend in 2-nitronaphthalene, 1-nitronaphthalene, and 2-methyl-1-nitronaphthalene, a series of nitro-substituted aromatic compounds. Comparison with the XMS-CASPT2//CASSCF (where CASSCF stands for complete active space self-consistent field method) results highlights the importance of dynamic correlation in geometry optimization and challenges the validity of an XMS-CASPT2//CASSCF approach: XMS-CASPT2 S1 minima leads to planar structures, while CASSCF optimizations trigger a pyramidalization of the nitro group. The XMS-CASPT2 results correlate the reported decreasing triplet quantum yield trend in these species to a decrease in S1 to T2 population transfer and an increase in S1-S0 decay, while no such correlation is observed when using XMS-CASPT2//CASSCF data.

Quantum-Chemistry Study of the Photophysical Properties of 4-Thiouracil and Comparisons with 2-Thiouracil.

DNA in living beings is constantly damaged by exogenous and endogenous agents. However, in some cases, DNA photodamage can have interesting applications, as it happens in photodynamic therapy. In this work, the current knowledge on the photophysics of 4-thiouracil has been extended by further quantum-chemistry studies to improve the agreement between theory and experiments, to better understand the differences with 2-thiouracil, and, last but not least, to verify its usefulness as a photosensitizer for photodynamic therapy. This study has been carried out by determining the most favorable deactivation paths of UV-vis photoexcited 4-thiouracil by means of the photochemical reaction path approach and an efficient combination of the complete-active-space second-order perturbation theory//complete-active-space self-consistent field (CASPT2//CASSCF), (CASPT2//CASPT2), time-dependent density functional theory (TDDFT), and spin-flip TDDFT (SF-TDDFT) methodologies. By comparing the data computed herein for both 4-thiouracil and 2-thiouracil, a rationale is provided on the relatively higher yields of intersystem crossing, triplet lifetime and singlet oxygen production of 4-thiouracil, and the relatively higher yield of phosphorescence of 2-thiouracil.

A First Proposal on the Nitrobenzene Photorelease Mechanism of NO2 and Its Relation to NO Formation through a Roaming Mechanism.

Despite the fact that NO2 is considered to be the main photoproduct of nitrobenzene photochemistry, no mechanism has ever been proposed to rationalize its formation. NO photorelease is instead a more studied process, probably due to its application in the drug delivery sector and the study of roaming mechanisms. In this contribution, a photoinduced mechanism accounting for the formation of NO2 in nitrobenzene is theorized based on CASPT2, CASSCF, and DFT electronic structure calculations and CASSCF classical dynamics. A triplet nπ* state is shown to evolve toward C-NO2 dissociation, being, in fact, the only low-lying excited state favoring such a deformation. Along the triplet dissociation path, the possibility to decay to the singlet ground state results in the frustration of the dissociation and in the recombination of the fragments, either back to the nitro or the nitrite isomer. The thermal decomposition of the latter to NO constitutes globally a roaming mechanism of NO formation.

Governing the emissive properties of 4-aminobiphenyl-2-pyrimidine push-pull systems via the restricted torsion of N,N-disubstituted amino groups.

Donor-acceptor-substituted biphenyl derivatives are particularly interesting model compounds, which exhibit intramolecular charge transfer because of the extent of charge transfer between both substituents. The connection of a 4-[1,1'-biphenyl]-4-yl-2-pyrimidinyl) moiety to differently disubstituted amino groups at the biphenyl terminal can offer push-pull compounds with distinctive photophysical properties. Herein, we report a comprehensive study of the influence of the torsion angle of the disubstituted amino group on the emissive properties of two pull-push systems: 4-[4-(4-N,N-dimethylaminophenyl)phenyl]-2,6-diphenylpyrimidine (D1) and 4-[4-(4-N,N-diphenylaminophenyl)phenyl]-2,6-diphenylpyrimidine (D2). The torsion angle of the disubstituted amino group, either N,N-dimethyl-amine or N,N-diphenyl-amine, at the biphenyl end governs their emissive properties. A drastic fluorescence quenching occurs in D1 as the solvent polarity increases, whereas D2 maintains its emission independently of the solvent polarity. Theoretical calculations on D1 support the presence of a twisted geometry for the lowest energy, charge-transfer excited state (S1,90), which corresponds to the minimum energy structure in polar solvents and presents a small energy barrier to move from the excited to the ground state, thereby favoring the non-radiative pathway and reducing the fluorescence efficiency. In contrast, this twisted structure is absent in D2 due to the steric hindrance of the phenyl groups attached to the amine group, making the non-radiative decay less favorable. Our findings provide insights into the crucial role of the substituent in the donor moiety of donor-acceptor systems on both the singlet excited state and the intramolecular charge-transfer process.

On the importance of equatorial metal-centered excited states in the photophysics of cyclometallated Ir(III) complexes.

In the present contribution, the following three cyclometallated Ir(III) complexes were theoretically investigated using density functional theory calculations to explain their different photophysical properties: [Ir(ppy)2(bpy)]+, where Hppy is 2-phenylpyridine and bpy is 2,2'-bipyridine, [Ir(ppy)2(pbpy)]+, where pbpy is 6-phenyl-2,2'-bipyridine, and [Ir(ppy)2(dpbpy)]+, where dpbpy is 6,6'-diphenyl-2,2'-bipyridine. Despite sharing the same molecular skeleton, with the only difference being the addition of one or two phenyl groups attached to the ancillary bpy ligand, the complexes show different emission quantum yields in CH2Cl2 solution (0.196, 0.049 and 0.036, respectively). Such a behavior was previously justified as a consequence of a different ability to non-radiatively decay through an axial metal-centered (MC) triplet state. In the present contribution, a new non-radiative decay path has been characterized to be mediated by the so-called equatorial MC states, in which an Ir-Nbpy bond is elongated instead of an Ir-Nppy bond as observed in the axial MC states. The decay path involving the equatorial MC states is more favorable than that associated with the axial MC states, and the different ability to decay through the former better explains the photoemission properties exhibited by the three complexes.

On the photorelease of nitric oxide by nitrobenzene derivatives: A CASPT2//CASSCF model.

Nitroaromatic compounds can photorelease nitric oxide after UV absorption. The efficiency of the photoreaction depends on the molecular structure, and two features have been pointed out as particularly important for the yield of the process: the presence of methyl groups at the ortho position with respect to the nitro group and the degree of conjugation of the molecule. In this paper, we provide a theoretical characterization at the CASPT2//CASSCF (complete active space second-order perturbation theory//complete active space self-consistent field) level of theory of the photorelease of NO for four molecules derived from nitrobenzene through the addition of ortho methyl groups and/or the elongation of the conjugation. Our previously described mechanism obtained for the photorelease of NO in nitrobenzene has been adopted as a model for the process. According to this model, the process proceeds through a reactive singlet-triplet crossing (STC) region that the system can reach from the triplet 3(πOπ*) minimum. The energy barrier that must be surmounted in order to populate the reactive STC can be associated with the efficiency of the photoreaction. Here, the obtained results display clear differences in the efficiency of the photoreaction in the studied systems and can be correlated with experimental results. Thus, the model proves its ability to highlight the differences in the photoreaction efficiency for the nitroaromatic compounds studied here.

On the Importance of Ligand-Centered Excited States in the Emission of Cyclometalated Ir(III) Complexes.

The photophysical behavior of the cyclometalating Ir(III) complexes [Ir(ppy)2(bpy)]+, where Hppy is 2-phenylpyridine and bpy is 2,2'-bipyridine (complex 1), and [Ir(diFppy)2(dtb-bpy)]+, where diFppy is 2-(2,4-difluorophenyl)pyridine and dtb-bpy is 4,4'-di-tert-butyl-2,2'-bipyridine (complex 2), has been theoretically investigated by performing density functional theory calculations. The two complexes share the same molecular skeleton, complex 2 being derived from complex 1 through the addition of fluoro and tert-butyl substituents, but present notable differences in their photophysical properties. The remarkable difference in their emission quantum yields (0.196 for complex 1 in dichloromethane and 0.71 for complex 2 in acetonitrile) has been evaluated by characterizing both radiative and nonradiative decay paths. It has emerged that the probability of decaying through the nonradiative triplet metal-centered state, normally associated with the loss of the emission quantum yield, does not appear to be the reason behind the reported substantially different emission efficiency. A more critical factor appears to be the ability of complex 2 to emit from both the usual metal-to-ligand charge-transfer state and from two additional ligand-centered states, as supported by the fact that the respective minima belong to the potential energy surface of the lowest triplet T1 state and that their phosphorescence lifetimes are in the same order of magnitude. In contrast, the emission of complex 1 can be originated only from the metal-to-ligand charge-transfer state, being the only emissive T1 minimum. The results constitute a significant case in which the emission from ligand-centered states is the key for determining the high emission quantum yield of a complex.

Building a Functionalizable, Potent Chemiluminescent Agent: A Rational Design Study on 6,8-Substituted Luminol Derivatives.

Luminol is a prominent chemiluminescent (CL) agent, finding applications across numerous fields, including forensics, immunoassays, and imaging. Different substitution patterns on the aromatic ring can enhance or decrease its CL efficiency. We herein report a systematic study on the synthesis and photophysics of all possible 6,8-disubstituted luminol derivatives bearing H, Ph, and/or Me substituents. Their CL responses are monitored at three pH values (8, 10, and 12), thus revealing the architecture with the optimum CL efficiency. The most efficient pattern is used for the synthesis of a strongly CL luminol derivative, bearing a functional group for further, straightforward derivatization. This adduct exhibits an unprecedented increase in chemiluminescence efficiency at pH = 12, pH = 10, and especially at pH = 8 (closer to the biologically relevant conditions) compared to luminol. Complementary work on the fluorescence of the emissive species as well as quantum chemistry computations are employed for the rationalization of the observed results.

On the chemiluminescence emission of luminol: protic and aprotic solvents and encapsulation to improve the properties in aqueous solution.

Luminol is a popular molecule that is currently gaining further interest due to its potential role for non-invasive cancer treatments. Design of more efficient derivatives in this context would benefit from a clear knowledge on the origin of the distinct intensity and spectroscopic properties in protic and aprotic solvents observed experimentally, which are still not rationalized. By efficiently combining molecular dynamics, quantum methodologies based on density functional theory and multiconfigurational quantum chemistry and hybrid approaches, and developing herein a computational approach for accurately determining "molar negative extinction (or gain) coefficients of emission", we firstly demonstrate that the amino and imino forms of the 3-aminophthalate dianion are responsible for the chemiluminescence in protic and aprotic medium, respectively. Secondly, we show that the coupling between the adjacent amino and carboxylate groups of luminol existing in aprotic solvents must be kept in aqueous solution to increase the chemiexcitation and emission intensity. Thirdly, modifications of luminol are proposed and simulated showing improved performances as compared to the parent molecule (stronger emission electronic transition and longer emission wavelengths) under the physiological conditions of relevance in biological and medical applications.

Regiochemical memory in the adiabatic photolysis of thymine-derived oxetanes. A combined ultrafast spectroscopic and CASSCF/CASPT2 computational study.

The photoinduced cycloreversion of oxetanes has been thoroughly investigated in connection with the photorepair of the well-known DNA (6-4) photoproducts. In the present work, the direct photolysis of the two regioisomers arising from the irradiation of benzophenone (BP) and 1,3-dimethylthymine (DMT), namely the head-to-head (HH-1) and head-to-tail (HT-1) oxetane adducts, has been investigated by combining ultrafast spectroscopy and theoretical multiconfigurational quantum chemistry analysis. Both the experimental and computational results agree with the involvement of an excited triplet exciplex 3[BPDMT]* for the photoinduced oxetane cleavage to generate 3BP* and DMT through an adiabatic photochemical reaction. The experimental signature of 3[BPDMT]* is the appearance of an absorption band at ca. 400 nm, detected by femtosecond transient absorption spectroscopy. Its formation is markedly regioselective, as it is more efficient and proceeds faster for HH-1 (∼2.8 ps) than for HT-1 (∼6.3 ps). This is in line with the theoretical analysis, which predicts an energy barrier to reach the triplet exciplex for HT-1, in contrast with a less hindered profile for HH-1. Finally, the more favorable adiabatic cycloreversion of HH-1 compared to that of HT-1 is explained by its lower probability to reach the intersystem crossing with the ground state, which would induce a radiationless deactivation process leading either to a starting adduct or to a dissociated BP and DMT.

On the Intrinsically Low Quantum Yields of Pyrimidine DNA Photodamages: Evaluating the Reactivity of the Corresponding Minimum Energy Crossing Points.

The low quantum yield of photoformation of cyclobutane pyrimidine dimers and pyrimidine-pyrimidone (6-4) adducts in DNA bases is usually associated with the presence of more favorable nonreactive decay paths and with the unlikeliness of exciting the system in a favorable conformation. Here, we prove that the ability of the reactive conical intersection to bring the system either back to the absorbing conformation or to the photoproduct must be considered as a fundamental factor in the low quantum yields of the mentioned photodamage. In support of the proposed model, the one order of magnitude difference in the quantum yield of formation of the cyclobutane thymine dimer with respect to the thymine-thymine (6-4) adduct is rationalized here by comparing the reactive ability of the seam of intersections leading respectively to the cyclobutane thymine dimer and the oxetane precursor of the thymine-thymine (6-4) adduct at the CASPT2 level of theory.

How important is roaming in the photodegradation of nitrobenzene?

At low excitation energies nitrobenzene photoreleases NO with low translational and rotational energy, while at higher excitation energies NO is photoreleased with both low and high translational and rotational energy. The fast products are formed through a singlet-triplet crossing (STC) region featuring an oxaziridine ring, while a ground state roaming mechanism was suggested to produce the slow molecules. Computing translational and rotational energies performing CASSCF classical dynamics, we here prove how the same oxaziridine STC can account for both fast and slow photoproducts, depending on the region of the seam through which the ground state is populated. A roaming-type STC/CI has also been characterized, from which slow NO molecules can also be formed through a roaming photodegradation mechanism, here in the excited state. The higher accessibility of the oxaziridine STC mechanism, 1.53 eV lower in energy than the roaming path, questions the contribution of roaming in nitrobenzene NO photoproduction.

Induced Night Vision by Singlet-Oxygen-Mediated Activation of Rhodopsin.

In humans, vision is limited to a small fraction of the whole electromagnetic spectrum. One possible strategy for enhancing vision in deep-red or poor-light conditions consists of recruiting chlorophyll derivatives in the rod photoreceptor cells of the eye, as suggested in the case of some deep-sea fish. Here, we employ all-atom molecular simulations and high-level quantum chemistry calculations to rationalize how chlorin e6 (Ce6), widely used in photodynamic therapy although accompanied by enhanced visual sensitivity, mediates vision in the dark, shining light on a fascinating but largely unknown molecular mechanism. First, we identify persistent interaction sites between Ce6 and the extracellular loops of rhodopsin, the transmembrane photoreceptor protein responsible for the first steps in vision. Triggered by Ce6 deep-red light absorption, the retinal within rhodopsin can be isomerized thus starting the visual phototransduction cascade. Our data largely exclude previously hypothesized energy-transfer mechanisms while clearly lending credence to a retinal isomerization indirectly triggered by singlet oxygen, proposing an alternative mechanism to rationalize photosensitizer-mediated night vision.

OpenMolcas: From Source Code to Insight.

In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.

Similar chemical structures, dissimilar triplet quantum yields: a CASPT2 model rationalizing the trend of triplet quantum yields in nitroaromatic systems.

The photophysics of nitroaromatic compounds is characterized by an ultrafast decay into the triplet manifold and by significant triplet quantum yields. The latter quantity changes drastically depending on the system, as shown for 2-nitronaphthalene, 1-nitronaphthalene, and 2-methyl-1-nitronaphthalene, whose triplet quantum yields have been previously measured to be 0.93 ± 0.15, 0.64 ± 0.12, and 0.33 ± 0.05, respectively (J. Phys. Chem. A, 2013, 117, 14100). In this study, we rationalize the reported trend of the triplet quantum yield on the basis of the different abilities of the excited S1 state to reach a previously unreported conical intersection with the ground state. This path is in competition with the path leading to the triplet state, which appears to be equally favorable in the three systems. The energy barriers from the S1 CASPT2//CASSCF minima to a CASPT2 minimum-energy-crossing-point of the S1/S0 conical intersection have been computed to follow the same trend as the triplet quantum yields of the nitroaromatic systems under analysis. The path has also been characterized for nitrobenzene; an energy barrier was obtained that nicely fits the derived model and is in agreement with its triplet quantum yield value (>0.8). The ability of the present model to not only rationalize the experimental data of a single molecule but also to reproduce a trend for four slightly different systems demonstrates its reliability.

Pyrene, a Test Case for Deep-Ultraviolet Molecular Photophysics.

We determined the complete relaxation dynamics of pyrene in ethanol from the second bright state, employing experimental and theoretical broadband heterodyne detected transient grating and two-dimensional photon echo (2DPE) spectroscopy, using pulses with duration of 6 fs and covering a spectral range spanning from 250 to 300 nm. Multiple lifetimes are assigned to conical intersections through a cascade of electronic states, eventually leading to a rapid population of the lowest long-living excited state and subsequent slow vibrational cooling. The lineshapes in the 2DPE spectra indicate that the efficiency of the population transfer depends on the kinetic energy deposited into modes required to reach a sloped conical intersection, which mediates the decay to the lowest electronic state. The presented experimental-theoretical protocol paves the way for studies on deep-ultraviolet-absorbing biochromophores ubiquitous in genomic and proteic systems.

Azole-containing cationic bis-cyclometallated iridium(iii) isocyanide complexes: a theoretical insight into the emission energy and emission efficiency.

Using a density functional theory approach, we explore the emission properties of a family of bis-cyclometallated cationic iridium(iii) complexes of general formula [Ir(C^N)2(CN-tert-Bu)2]+ that have tert-butyl isocyanides as neutral auxiliary ligands. Taking the [Ir(ppy)2(CN-tert-Bu)2]+ complex (Hppy = 2-phenylpyridine) as a reference, the effect of replacing the pyridine ring in the cyclometallating ppy ligand by a five-membered azole ring has been examined. To this end, two series of complexes differing by the nature of the atom (either nitrogen or carbon) linking the azole to the phenyl ring of the cyclometallating ligand have been designed. Each series is composed of three molecules having an increasing number of nitrogen atoms (2 to 4) in the azole ring. The emission energies computed for the azole-containing [Ir(C^N)2(CN-tert-Bu)2]+ complexes show a generalized blue-shift compared to [Ir(ppy)2(CN-tert-Bu)2]+, in agreement with the experimental data available for two of the six complexes designed here. The electronic nature of the lowest-lying triplet (T1) is clearly established as a ligand-centred (3LC) state associated with the cyclometallating ligands, and cannot be described as a simple HOMO → LUMO promotion. Therefore, no clear trend based on the sole use of molecular orbital energies can be inferred to predict the emission properties. The significant oscillation in the emission quantum yield (ranging from 0.1% to 52%) experimentally reported is rationalized by the energy gap between the emitting T1 state and a non-radiative triplet state having metal-centred (3MC) d-d* nature. On the basis of such a model, two of the here proposed systems are expected to display significant emission quantum yields in the blue region of the visible spectra, which make them good candidates for electroluminescent applications.

Molecular Basis of the Chemiluminescence Mechanism of Luminol.

Light emission from luminol is probably one of the most popular chemiluminescence reactions due to its use in forensic science, and has recently displayed promising applications for the treatment of cancer in deep tissues. The mechanism is, however, very complex and distinct possibilities have been proposed. By efficiently combining DFT and CASPT2 methodologies, the chemiluminescence mechanism has been studied in three steps: 1) luminol oxygenation to generate the chemiluminophore, 2) a chemiexcitation step, and 3) generation of the light emitter. The findings demonstrate that the luminol double-deprotonated dianion activates molecular oxygen, diazaquinone is not formed, and the chemiluminophore is formed through the concerted addition of oxygen and concerted elimination of nitrogen. The peroxide bond, in comparison to other isoelectronic chemical functionalities (-NH-NH-, -N- -N- -, and -S-S-), is found to have the best chemiexcitation efficiency, which allows the oxygenation requirement to be rationalized and establishes general design principles for the chemiluminescence efficiency. Electron transfer from the aniline ring to the OO bond promotes the excitation process to create an excited state that is not the chemiluminescent species. To produce the light emitter, proton transfer between the amino and carbonyl groups must occur; this requires highly localized vibrational energy during chemiexcitation.