Juan J. Nogueira

Dynamics of CO2 Scattering off a Perfluorinated Self-Assembled Monolayer. Influence of the Incident Collision Energy, Mass Effects, and Use of Different Surface Models †

The dynamics of collisions of CO2 with a perfluorinated alkanethiol self-assembled monolayer (F-SAM) on gold were investigated by classical trajectory calculations using explicit atom (EA) and united atom (UA) models to represent the F-SAM surface. The CO2 molecule was directed perpendicularly to the surface at initial collision energies of 1.6, 4.7, 7.7, and 10.6 kcal/mol. Rotational distributions of the scattered CO2 molecules are in agreement with experimental distributions determined for collisions of CO2 with liquid surfaces of perfluoropolyether. The agreement is especially good for the EA model. The role of the mass in the efficiency of the energy transfer was investigated in separate simulations in which the mass of the F atoms was replaced by either that of hydrogen or chlorine, while keeping the potential energy function unchanged. The calculations predict the observed trend that less energy is transferred to the surface as the mass of the alkyl chains increases. Significant discrepancies were found between results obtained with the EA and UA models. The UA surface leads to an enhancement of the energy transfer efficiency in comparison with the EA surface. The reason for this is in the softer structure of the UA surface, which facilitates transfer from translation to interchain vibrational modes.

Molecular dynamics simulations of binding modes between methylene blue and DNA with alternating GC and AT sequences.

The understanding of interactions between small molecules and DNA is crucial to design new anticancer drugs targeted to DNA. Methylene blue (MB) is a phenothiazinium dye that has shown promising results in photodynamic therapy treatment. The noncovalent binding of methylene blue to DNA was experimentally and theoretically analyzed in the past, but certain features of the binding mode are still not clear. In this work, force field molecular dynamics simulations were performed to simulate the binding of methylene blue to alternating GC and AT sequences at two different ionic strengths. External, intercalative, minor groove, and major groove binding modes are discussed based on energetic and structural analyses. External and major groove complexes were found to be unstable structures, although for poly(dA-dT) the major groove binding mode cannot be discarded, especially at high ionic strengths. Minor groove and intercalative binding leads to stable adducts. The most energetically favorable orientation of the dye inside the minor groove is different for the two DNA sequences because of the different balances between the DNA deformation energy and the dye/DNA interaction energy. The intercalative binding is the most important interaction mode. The dye undergoes rotational transitions inside the intercalative pocket for both DNA sequences giving rise to three dye/DNA adducts that have different energetic and structural features. This rotational motion explains the different behavior found in experiments for the GC and AT nucleic acids at different ionic strengths.

Enhancing Intersystem Crossing in Phenotiazinium Dyes by Intercalation into DNA

Phenothiazinium dyes are used as photosensitizers in photodynamic therapy. Their mode of action is related to the generation of triplet excited states by intersystem crossing. Therefore, rationalizing the factors that influence intersystem crossing is crucial to improve the efficacy of photodynamic therapy. Here we employ quantum mechanics/molecular mechanics calculations to investigate the effect of aqueous and nucleic acid environments on the intersystem crossing mechanism in methylene blue. We find that the mechanism by which the triplet states are generated depends strongly on the environment. While intersystem crossing in water is mediated exclusively by vibronic spin-orbit coupling, it is enhanced in DNA due to a second pathway driven by electronic spin-orbit coupling. Competing charge-transfer processes, which are also possible in the presence of DNA, can therefore be suppressed by a suitable structural functionalization, thereby increasing the efficacy of photodynamic therapy.

Quenching of Charge Transfer in Nitrobenzene Induced by Vibrational Motion

Although nitrobenzene is the smallest nitro-aromatic molecule, the nature of its electronic structure is still unclear. Most notably, the lowest-energy absorption band was assessed in numerous studies providing conflicting results regarding its charge-transfer character. In this study, we employ a combination of molecular dynamics and quantum chemical methods to disentangle the nature of the lowest-energy absorption band of nitrobenzene. Surprisingly, the charge-transfer transition from the benzene moiety to the nitro group is found to be quenched by a flow of charge into the opposite direction induced by vibrational motion. Beyond clarifying the charge-transfer character of nitrobenzene, we show that the widely used approach of analyzing the ground-state minimum-energy geometry provides completely wrong conclusions.

Computational Photophysics in the Presence of an Environment

Most processes triggered by ultraviolet (UV) or visible (vis) light in nature take place in complex biological environments. The first step in these photophysical events is the excitation of the absorbing system or chromophore to an electronically excited state. Such an excitation can be monitored by the UV-vis absorption spectrum. A precise calculation of the UV-vis spectrum of a chromophore embedded in an environment is a challenging task that requires the consideration of several ingredients, besides an accurate electronic-structure method for the excited states. Two of the most important are an appropriate description of the interactions between the chromophore and the environment and accounting for the vibrational motion of the whole system. In this contribution, we review the most common theoretical methodologies to describe the environment (including quantum mechanics/continuum and quantum mechanics/molecular mechanics models) and to account for vibrational sampling (including Wigner sampling and molecular dynamics). Further, we illustrate in a series of examples how the lack of these ingredients can lead to a wrong interpretation of the electronic features behind the UV-vis absorption spectrum.

Mechanism of Ultrafast Intersystem Crossing in 2-Nitronaphthalene

Abstract Nitronaphthalene derivatives efficiently populate their electronically excited triplet states upon photoexcitation through ultrafast intersystem crossing (ISC). Despite having been studied extensively by time‐resolved spectroscopy, the reasons behind their ultrafast ISC remain unknown. Herein, we present the first ab initio nonadiabatic molecular dynamics study of a nitronaphthalene derivative, 2‐nitronaphthalene, including singlet and triplet states. We find that there are two distinct ISC reaction pathways involving different electronic states at distinct nuclear configurations. The high ISC efficiency is explained by the very small electronic and nuclear alterations that the chromophore needs to undergo during the singlet–triplet transition in the dominating ISC pathway after initial dynamics in the singlet manifold. The insights gained in this work are expected to shed new light on the photochemistry of other nitro polycyclic aromatic hydrocarbons that exhibit ultrafast intersystem crossing.

Sequential Proton-Coupled Electron Transfer Mediates Excited-State Deactivation of a Eumelanin Building Block

Skin photoprotection is commonly believed to rely on the photochemistry of 5,6-dihydroxyindole (DHI)- and 5,6-dihydroxyindole-2-carboxylic acid (DHICA)-based eumelanin building blocks. Attempts to elucidate the underlying excited-state relaxation mechanisms have been partly unsuccessful due to the marked instability to oxidation. We report a study of the excited-state deactivation of DHI using steady-state and time-resolved fluorescence accompanied by high-level quantum-chemistry calculations including solvent effects. Spectroscopic data show that deactivation of the lowest excited state of DHI in aqueous buffer proceeds on the 100 ps time scale and is 20 times faster than in methanol. Quantum-chemical calculations reveal that the excited-state decay mechanism is a sequential proton-coupled electron transfer, which involves the initial formation of a solvated electron from DHI, followed by the transfer of a proton to the solvent. This unexpected finding would prompt a revision of current notions about eumelanin photophysics and photobiology.

Challenges in Simulating Light-Induced Processes in DNA

In this contribution, we give a perspective on the main challenges in performing theoretical simulations of photoinduced phenomena within DNA and its molecular building blocks. We distinguish the different tasks that should be involved in the simulation of a complete DNA strand subject to UV irradiation: (i) stationary quantum chemical computations; (ii) the explicit description of the initial excitation of DNA with light; (iii) modeling the nonadiabatic excited state dynamics; (iv) simulation of the detected experimental observable; and (v) the subsequent analysis of the respective results. We succinctly describe the methods that are currently employed in each of these steps. While for each of them, there are different approaches with different degrees of accuracy, no feasible method exists to tackle all problems at once. Depending on the technique or combination of several ones, it can be problematic to describe the stacking of nucleobases, bond breaking and formation, quantum interferences and tunneling or even simply to characterize the involved wavefunctions. It is therefore argued that more method development and/or the combination of different techniques are urgently required. It is essential also to exercise these new developments in further studies on DNA and subsystems thereof, ideally comprising simulations of all of the different components that occur in the corresponding experiments.

Stepwise photosensitized thymine dimerization mediated by an exciton intermediate

Cyclobutane thymine dimerization is the most prominent DNA photoinduced damage. While the ultrafast mechanism that proceeds in the singlet manifold is nowadays well established, the triplet-state pathway is not completely understood. Here we report the underlying mechanism of the photosensitized dimerization process in the triplet state. Quantum chemical calculations, combined with wavefunction analysis, and nonadiabatic molecular dynamics simulations demonstrate that this is a stepwise reaction, traversing a long-lived triplet biradical intermediate, which is characterized as a Frenkel exciton with very small charge-transfer character. The low yield of the reaction is regulated by two factors: (i) a relatively large energy barrier that needs to be overcome to form the exciton intermediate, and (ii) a bifurcation of the ground-state potential-energy surface that mostly leads back to the Franck–Condon region because dimerization requires a very restricted combination of coordinates and velocities at the event of non-radiative decay to the ground state.Graphical abstract

Exciton Localization on Ru-Based Photosensitizers Induced by Binding to Lipid Membranes

The characterization of electronic properties of metal complexes embedded in membrane environments is of paramount importance to develop efficient photosensitizers in optogenetic applications. Molecular dynamics and QM/MM simulations together with quantitative wave function analysis reveal a directional electronic redistribution of the exciton formed upon excitation of [Ru(bpy)2(bpy-C17)]2+ when going from water to a lipid bilayer, despite the fact that the media influence neither the metal-to-ligand charge-transfer character nor the excitation energy of the absorption spectra. When the photosensitizer is embedded into the DOPC lipid membrane, exciton population is mainly located in the bypyridyl sites proximal to the positively charged surface of the bilayer due to electrostatic interactions. This behavior shows that the electronic structure of metal complexes can be controlled through the binding to external species, underscoring the crucial role of the environment in directing the electronic flow upon excitation and thus helping rational tuning of optogenetic agents.