Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations

Coupled Nuclear and Electron Dynamics in Chlorophyll Unraveled by XMS-CASPT2 X-ray Absorption Spectra

Attosecond spectroscopy, especially time-resolved X-ray absorption spectra (XAS), enables direct observation of ultrafast molecular dynamics. The complementary and even preceding development of theoretical simulations can offer the necessary guidance and stimulate new experiments. In this work, we simulated high-level XAS for the magnesium and nitrogen K-edge of chlorophyll a. In our previous work on the ultrafast relaxation process in the Q-band, our quantum dynamics simulations found the Qx and Qy states to be energetically close and therefore strongly coupled. Here, we analyze the strong coupling between Qx and Qy via XAS, indicating promising possibilities for experimental observation. The excited-state energies, potential energy surfaces, and XAS are computed at the XMS-CASPT2 level of theory to capture the complex multireference character of chlorophyll excitations. In our simulated spectra, we could follow the ultrafast population transfer between Qx and Qy and thus draw conclusions about the strong vibrational coupling between them.

The Q-Band Energetics and Relaxation of Chlorophylls a and b as Revealed by Visible-to-Near Infrared Time-Resolved Absorption Spectroscopy

Chlorophyll (Chl) is the most abundant light-harvesting pigment of oxygenic photosynthetic organisms; however, the Q-band energetics and relaxation dynamics remain unclear. In this work, we have applied femtosecond time-resolved (fs-TA) absorption spectroscopy in 430-1,700 nm to Chls a and b in diluted pyridine solutions under selective optical excitation within their Q-bands. The results revealed distinct near-infrared absorption features of the Bx,y ← Qy and Bx,y ← Qx transitions in 930-1,700 nm, which together with the steady-state absorption in 400-700 nm unveiled the Qx(0,0)-state energy that lies 1,000 ± 400 and 600 ± 400 cm-1 above the Qy(0,0)-state for Chls a and b, respectively. In addition, the Qx-to-Qy internal conversion time constants are estimated to be less than 80 fs for Chls a and b. These findings may shed light on understanding the roles of the Chls in the primary excitation energy transfer reactions of photosynthesis.

Complementing Adiabatic and Nonadiabatic Methods To Understand Internal Conversion Dynamics in Porphyrin Derivatives

We analyze the internal conversion dynamics within the Qy and Qx excited states of both bare and functionalized porphyrins, which are known to exhibit significantly different time constants experimentally. Through the integration of two complementary approaches, static calculation of per-mode reorganization energies and nonadiabatic molecular dynamics, we achieve a comprehensive understanding of the factors determining the different behavior of the two molecules. We identify the key normal and essential modes responsible for the population transfer between excited states and discuss the efficacy of different statistical and nonstatistical analyses in providing a full physics-based description of the phenomenon.

Computational Exploration of the Mechanism of Action of a Sorafenib-Containing Ruthenium Complex as an Anticancer Agent for Photoactivated Chemotherapy

Ruthenium(II) polypyridyl complexes are being tested as potential anticancer agents in different therapies, which include conventional chemotherapy and light-activated approaches. A mechanistic study on a recently synthesized dual-action Ru(II) complex [Ru(bpy)2(sora)Cl]+ is described here. It is characterized by two mono-dentate leaving ligands, namely, chloride and sorafenib ligands, which make it possible to form a di-aquo complex able to bind DNA. At the same time, while the released sorafenib can induce ferroptosis, the complex is also able to act as a photosensitizer according to type II photodynamic therapy processes, thus generating one of the most harmful cytotoxic species, 1O2. In order to clarify the mechanism of action of the drug, computational strategies based on density functional theory are exploited. The photophysical properties of the complex, which include the absorption spectrum, the kinetics of ISC, and the character of all the excited states potentially involved in 1O2 generation, as well as the pathway providing the di-aquo complex, are fully explored. Interestingly, the outcomes show that light is needed to form the mono-aquo complex, after releasing both chloride and sorafenib ligands, while the second solvent molecule enters the coordination sphere of the metal once the system has come back to the ground-state potential energy surface. In order to simulate the interaction with canonical DNA, the di-aquo complex interaction with a guanine nucleobase as a model has also been studied. The whole study aims to elucidate the intricate details of the photodissociation process, which could help with designing tailored metal complexes as potential anticancer agents.

Molecular ordering-enhanced circularly polarized luminescence of chiral 1,10-phenanthroline derivatives

2,9-Bis((1R,2S,5R)-2-isopropyl-5-methylcyclohexanoxy)-1,10-phenanthroline [2,9-di-L-menthoxy-1,10-phenanthroline] (Men2Phen) and 2,9-bis(2-(S)-methylbutoxy)-1,10-phenanthroline (MB2Phen) were synthesized as chiral derivatives of 1,10-phenanthroline (Phen). Differences in rigidity and bulkiness of the chiral substituents at the 2- and 9-positions of the Phen backbone led to distinctive molecular dissymmetry in the ground state resulting in remarkable differences in circular dichroism. Men2Phen exhibited efficient circularly polarized luminescence (CPLm) at an anisotropy factor of 10-2 in the solid state based on molecular ordering disclosed by X-ray crystal analysis, while it showed much lower anisotropy factor in solution. MB2Phen, which was rather amorphous and did not afford good crystals, showed only negligible CPLm both in the solid state and in solution.

Unraveling Copper Exchange in the Atox1-Cu(I)-Mnk1 Heterodimer: A Simulation Approach

Copper, an essential metal for various cellular processes, requires tight regulation to prevent cytotoxicity. Intracellular pathways crucial for maintaining optimal copper levels involve soluble and membrane transporters, namely, metallochaperones and P-type ATPases, respectively. In this study, we used a simulation workflow based on free-energy perturbation (FEP) theory and parallel bias metadynamics (PBMetaD) to predict the Cu(I) exchange mechanism between the human Cu(I) chaperone, Atox1, and one of its two physiological partners, ATP7A. ATP7A, also known as the Menkes disease protein, is a transmembrane protein and one of the main copper-transporting ATPases. It pumps copper into the trans-Golgi network for the maturation of cuproenzymes and is also essential for the efflux of excess copper across the plasma membrane. In this analysis, we utilized the nuclear magnetic resonance (NMR) structure of the Cu(I)-mediated complex between Atox1 and the first soluble domain of the Menkes protein (Mnk1) as a starting point. Independent free-energy simulations were conducted to investigate the dissociation of both Atox1 and Mnk1. The calculations revealed that the two dissociations require free energy values of 6.3 and 6.2 kcal/mol, respectively, following a stepwise dissociation mechanism.

Excited-state normal-mode analysis: The case of porphyrins

We systematically applied excited-state normal mode analysis to investigate and compare the relaxation and internal conversion dynamics of a free-base porphyrin (BP) with those of a novel functional porphyrin (FP) derivative. We discuss the strengths and limitations of this method and employ it to predict very different dynamical behaviors of the two compounds and to clarify the role of high reorganization energy modes in driving the system toward critical regions of the potential energy landscape. We identify the modes of vibrations along which the energy gap between two excited-state potential energy surfaces within the Q band manifold may vanish and find that the excess energy to reach this "touching" region is significantly reduced in the case of FP (0.16 eV) as compared to the one calculated for BP (0.92 eV). Our findings establish a link between the chemical functionalization and the electronic and vibrational structure that can be exploited to control the internal conversion pathways in a systematic way.

Simulation workflows to predict the circular dichroism and circularly polarized luminescence of chiral materials

Chiral materials are attracting considerable interest in various fields in view of their unique properties and optical activity. Indeed, the peculiar features of chiral materials to absorb and emit circularly polarized light enable their use in an extensive range of applications. Motivated by the interest in boosting the development of chiral materials characterized by enhanced chiroptical properties such as circular dichroism (CD) and circular polarized luminescence (CPL), we herein illustrate in this tutorial how theoretical simulations can be used for the predictions and interpretations of chiroptical data and for the identification of chiral geometries. We are focusing on computational frameworks that can be used to investigate the theoretical aspects of chiral materials' photophysical and conformational characteristics. We will then illustrate ab initio methods based on density functional theory (DFT) and its time-dependent extension (TD-DFT) to simulate CD and CPL signals, and we will exemplify a variety of enhanced sampling techniques useful for an adequate sampling of the configurational space for chiral systems.

Signatures of intramolecular vibrational and vibronic Q x -Q y coupling effects in absorption and CD spectra of chlorophyll dimers

An electron-vibrational coupling model that includes the vibronic (non-adiabatic) coupling between the Qy and Qx transitions of chlorophyll (Chl), created by Reimers and coworkers (Scientific Rep. 3, 2761, 2013) is extended here to chlorophyll dimers with interchlorophyll excitonic coupling. The model is applied to a Chl a dimer of the water-soluble chlorophyll binding protein (WSCP). As for isolated chlorophyll, the vibronic coupling is found to have a strong influence on the high-frequency vibrational sideband in the absorption spectrum, giving rise to a band splitting. In contrast, in the CD spectrum the interplay of vibronic coupling and static disorder leads to a strong suppression of the vibrational sideband in excellent agreement with the experimental data. The conservative nature of the CD spectrum in the low-energy region is found to be caused by a delicate balance of the intermonomer excitonic coupling between the purely electronic Qy transition and the Qy transition involving intramolecular vibrational excitations on one hand and the coupling to higher-energy electronic transitions on the other hand.

Thermal site energy fluctuations in photosystem I: new insights from MD/QM/MM calculations

Cyanobacterial photosystem I (PSI) is one of the most efficient photosynthetic machineries found in nature. Due to the large scale and complexity of the system, the energy transfer mechanism from the antenna complex to the reaction center is still not fully understood. A central element is the accurate evaluation of the individual chlorophyll excitation energies (site energies). Such an evaluation must include a detailed treatment of site specific environmental influences on structural and electrostatic properties, but also their evolution in the temporal domain, because of the dynamic nature of the energy transfer process. In this work, we calculate the site energies of all 96 chlorophylls in a membrane-embedded model of PSI. The employed hybrid QM/MM approach using the multireference DFT/MRCI method in the QM region allows to obtain accurate site energies under explicit consideration of the natural environment. We identify energy traps and barriers in the antenna complex and discuss their implications for energy transfer to the reaction center. Going beyond previous studies, our model also accounts for the molecular dynamics of the full trimeric PSI complex. Via statistical analysis we show that the thermal fluctuations of single chlorophylls prevent the formation of a single prominent energy funnel within the antenna complex. These findings are also supported by a dipole exciton model. We conclude that energy transfer pathways may form only transiently at physiological temperatures, as thermal fluctuations overcome energy barriers. The set of site energies provided in this work sets the stage for theoretical and experimental studies on the highly efficient energy transfer mechanisms in PSI.

Q-Band relaxation in chlorophyll: new insights from multireference quantum dynamics

The ultrafast relaxation within the Q-bands of chlorophyll plays a crucial role in photosynthetic light-harvesting. We investigate this process via nuclear and electronic quantum dynamics on multireference potential energy surfaces.