Ultrafast infrared observation of exciton equilibration from oriented single crystals of photosystem II

Excitation landscape of the CP43 photosynthetic antenna complex from multiscale simulations

Photosystem II (PSII), the principal enzyme of oxygenic photosynthesis, contains two integral light harvesting proteins (CP43 and CP47) that bind chlorophylls and carotenoids. The two intrinsic antennae play crucial roles in excitation energy transfer and photoprotection. CP43 interacts most closely with the reaction center of PSII, specifically with the branch of the reaction center (D1) that is responsible for primary charge separation and electron transfer. Deciphering the function of CP43 requires detailed atomic-level insights into the properties of the embedded pigments. To advance this goal, we employ a range of multiscale computational approaches to determine the site energies and excitonic profile of CP43 chlorophylls, using large all-atom models of a membrane-bound PSII monomer. In addition to time-dependent density functional theory (TD-DFT) used in the context of a quantum-mechanics/molecular-mechanics setup (QM/MM), we present a thorough analysis using the perturbed matrix method (PMM), which enables us to utilize information from long-timescale molecular dynamics simulations of native PSII-complexed CP43. The excited state energetics and excitonic couplings have both similarities and differences compared with previous experimental fits and theoretical calculations. Both static TD-DFT and dynamic PMM results indicate a layered distribution of site energies and reveal specific groups of chlorophylls that have shared contributions to low-energy excitations. Importantly, the contribution to the lowest energy exciton does not arise from the same chlorophylls at each system configuration, but rather changes as a function of conformational dynamics. An unexpected finding is the identification of a low-energy charge-transfer excited state within CP43 that involves a lumenal (C2) and the central (C10) chlorophyll of the complex. The results provide a refined basis for structure-based interpretation of spectroscopic observations and for further deciphering excitation energy transfer in oxygenic photosynthesis.

How can Phycobilisome, the unique light harvesting system in certain algae working highly efficiently: The connection in between structures and functions

Algae, which are ubiquitous in ecosystems, have evolved a variety of light-harvesting complexes to better adapt to diverse habitats. Phycobilisomes/phycobiliproteins, unique to cyanobacteria, red algae, and certain cryptomonads, compensate for the lack of chlorophyll absorption, allowing algae to capture and efficiently transfer light energy in aquatic environments. With the advancement of microscopy and spectroscopy, the structure and energy transfer processes of increasingly complex phycobilisomes have been elucidated, providing us with a vivid portrait of the dynamic adaptation of their structures to the light environment in which algae thrive: 1) Cyanobacteria living on the surface of the water use short, small phycobilisomes to absorb red-orange light and reduce the damage from blue-violet light via multiple methods; 2) Large red algae inhabiting the depths of the ocean have evolved long and dense phycobilisomes containing phycoerythrin to capture the feeble blue-green light; 3) In far-red light environments such as caves, algae use special allophycocyanin cores to optimally utilize the far-red light; 4) When the environment shifts, algae can adjust the length, composition and density of their rods to better adapt; 5) By carefully designing the position of the pigments, phycobilisomes can transfer light energy to the reaction center with nearly 100% efficiency via three energy transfer processes.

Four-dimensional microED of conformational dynamics in protein microcrystals on the femto-to-microsecond timescales

As structural determination of protein complexes approaches atomic resolution, there is an increasing focus on conformational dynamics. Here we conceptualize the combination of two techniques which have become established in recent years: microcrystal electron diffraction and ultrafast electron microscopy. We show that the extremely low dose of pulsed photoemission still enables microED due to the strength of the electron bunching from diffraction of the protein crystals. Indeed, ultrafast electron diffraction experiments on protein crystals have already been demonstrated to be effective in measuring intermolecular forces in protein microcrystals. We discuss difficulties that may arise in the acquisition and processing of data and the overall feasibility of the experiment, paying specific attention to dose and signal-to-noise ratio. In doing so, we outline a detailed workflow that may be effective in minimizing the dose on the specimen. A series of model systems that would be good candidates for initial experiments is provided.

Alternative Fast and Slow Primary Charge-Separation Pathways in Photosystem II

Photosystem-II (PSII) is a multi-subunit protein complex that harvests sunlight to perform oxygenic photosynthesis. Initial light-activated charge separation takes place at a reaction centre consisting of four chlorophylls and two pheophytins. Understanding the processes following light excitation remains elusive due to spectral congestion, the ultrafast nature, and multi-component behaviour of the charge-separation process. Here, using advanced computational multiscale approaches which take into account the large-scale configurational flexibility of the system, we identify two possible primary pathways to radical-pair formation that differ by three orders of magnitude in their kinetics. The fast (short-range) pathway is dominant, but the existence of an alternative slow (long-range) charge-separation pathway hints at the evolution of redundancy that may serve other purposes, adaptive or protective, related to formation of the unique oxidative species that drives water oxidation in PSII.

Light-induced reversible reorganizations in closed Type II reaction centre complexes: physiological roles and physical mechanisms

The purpose of this review is to outline our understanding of the nature, mechanism and physiological significance of light-induced reversible reorganizations in closed Type II reaction centre (RC) complexes. In the so-called 'closed' state, purple bacterial RC (bRC) and photosystem II (PSII) RC complexes are incapable of generating additional stable charge separation. Yet, upon continued excitation they display well-discernible changes in their photophysical and photochemical parameters. Substantial stabilization of their charge-separated states has been thoroughly documented-uncovering light-induced reorganizations in closed RCs and revealing their physiological importance in gradually optimizing the operation of the photosynthetic machinery during the dark-to-light transition. A range of subtle light-induced conformational changes has indeed been detected experimentally in different laboratories using different bRC and PSII-containing preparations. In general, the presently available data strongly suggest similar structural dynamics of closed bRC and PSII RC complexes, and similar physical mechanisms, in which dielectric relaxation processes and structural memory effects of proteins are proposed to play important roles.

Exciton quenching by oxidized chlorophyll Z across the two adjacent monomers in a photosystem II core dimer

This study aimed to clarify (1) which pigment in a photosystem II (PSII) core complex is responsible for the 695-nm emission at 77 K and (2) the molecular basis for the oxidation-induced fluorescence quenching in PSII. Picosecond time-resolved fluorescence dynamics was compared between the dimeric and monomeric PSII with and without addition of an oxidant. The results indicated that the excitation-energy flow to the 695-nm-emitting chlorophyll (Chl) at 36 K and 77 K was hindered upon monomerization, clearly demonstrating significant exciton migration from the Chls on one monomer to the 695-nm-emitting pigment on the adjacent monomer. Oxidation of the redox-active Chl, which is named Chl_Z caused almost equal quenching of the 684-nm and 695-nm emission bands in the dimer, and lower quenching of the 695-nm band in the monomer. These results suggested two possible scenarios responsible for the 695-nm emission band: (A) Chl11-13 pair and the oxidized Chl_ZD1 work as the 695-nm emitting Chl and the quenching site, respectively, and (B) Chl29 and the oxidized Chl_ZD2 work as the 695-nm emitting Chl and the quenching site, respectively.

Observation of Cation Chromophore Photoisomerization of a Fluorescent Protein Using Millisecond Synchrotron Serial Crystallography and Infrared Vibrational and Visible Spectroscopy

The chromophores of reversibly switchable fluorescent proteins (rsFPs) undergo photoisomerization of both the trans and cis forms. Concurrent with cis/trans photoisomerisation, rsFPs typically become protonated on the phenolic oxygen resulting in a blue shift of the absorption. A synthetic rsFP referred to as rsEospa, derived from EosFP family, displays the same spectroscopic behavior as the GFP-like rsFP Dronpa at pH 8.4 and involves the photoconversion between nonfluorescent neutral and fluorescent anionic chromophore states. Millisecond time-resolved synchrotron serial crystallography of rsEospa at pH 8.4 shows that photoisomerization is accompanied by rearrangements of the same three residues as seen in Dronpa. However, at pH 5.5 we observe that the OFF state is identified as the cationic chromophore with additional protonation of the imidazolinone nitrogen which is concurrent with a newly formed hydrogen bond with the Glu212 carboxylate side chain. FTIR spectroscopy resolves the characteristic up-shifted carbonyl stretching frequency at 1713 cm^-1 for the cationic species. Electronic spectroscopy furthermore distinguishes the cationic absorption band at 397 nm from the neutral species at pH 8.4 seen at 387 nm. The observation of photoisomerization of the cationic chromophore state demonstrates the conical intersection for the electronic configuration, where previously fluorescence was proposed to be the main decay route for states containing imidazolinone nitrogen protonation. We present the full time-resolved room-temperature X-ray crystallographic, FTIR, and UV/vis assignment and photoconversion modeling of rsEospa.

From antenna to reaction center: Pathways of ultrafast energy and charge transfer in photosystem II

The photosystem II core complex (PSII-CC) is a photosynthetic complex that contains antenna proteins, which collect energy from sunlight, and a reaction center, which converts the collected energy to redox potential. Understanding the interplay between the antenna proteins and the reaction center will facilitate the development of more efficient solar energy conversion technologies. Here, we study the sub-100-ps dynamics of PSII-CC with two-dimensional electronic-vibrational spectroscopy, which connects energy flows with physical space, allowing a direct mapping of energy transfer pathways. Our results reveal a complex dynamical scheme which includes a specific pathway that connects CP43 to the reaction center. Resolving this pathway experimentally provides insights into the energy conversion processes in natural photosynthesis.

The photosystem II core complex (PSII-CC) is the smallest subunit of the oxygenic photosynthetic apparatus that contains core antennas and a reaction center, which together allow for rapid energy transfer and charge separation, ultimately leading to efficient solar energy conversion. However, there is a lack of consensus on the interplay between the energy transfer and charge separation dynamics of the core complex. Here, we report the application of two-dimensional electronic-vibrational (2DEV) spectroscopy to the spinach PSII-CC at 77 K. The simultaneous temporal and spectral resolution afforded by 2DEV spectroscopy facilitates the separation and direct assignment of coexisting dynamical processes. Our results show that the dominant dynamics of the PSII-CC are distinct in different excitation energy regions. By separating the excitation regions, we are able to distinguish the intraprotein dynamics and interprotein energy transfer. Additionally, with the improved resolution, we are able to identify the key pigments involved in the pathways, allowing for a direct connection between dynamical and structural information. Specifically, we show that C505 in CP43 and the peripheral chlorophyll Chlz_D1 in the reaction center are most likely responsible for energy transfer from CP43 to the reaction center.

Ultrafast excitation quenching by the oxidized photosystem II reaction center

Photosystem II (PSII) is the pigment–protein complex driving the photoinduced oxidation of water and reduction of plastoquinone in all oxygenic photosynthetic organisms. Excitations in the antenna chlorophylls are photochemically trapped in the reaction center (RC) producing the chlorophyll–pheophytin radical ion pair P+ Pheo−. When electron donation from water is inhibited, the oxidized RC chlorophyll P+ acts as an excitation quencher, but knowledge on the kinetics of quenching is limited. Here, we used femtosecond transient absorption spectroscopy to compare the excitation dynamics of PSII with neutral and oxidized RC (P+). We find that equilibration in the core antenna has a major lifetime of about 300 fs, irrespective of the RC redox state. Two-dimensional electronic spectroscopy revealed additional slower energy equilibration occurring on timescales of 3–5 ps, concurrent with excitation trapping. The kinetics of PSII with open RC can be described well with previously proposed models according to which the radical pair P+ Pheo− is populated with a main lifetime of about 40 ps, which is primarily determined by energy transfer between the core antenna and the RC chlorophylls. Yet, in PSII with oxidized RC (P+), fast excitation quenching was observed with decay lifetimes as short as 3 ps and an average decay lifetime of about 90 ps, which is shorter than the excited-state lifetime of PSII with open RC. The underlying mechanism of this extremely fast quenching prompts further investigation.

Independent Trajectory Mixed Quantum-Classical Approaches Based on the Exact Factorization

Mixed quantum-classical dynamics based on the exact factorization exploits the "derived" electron-nuclear correlation (ENC) term aiming for the description of quantum coherences. The ENC contains interactions between the phase of electronic states and nuclear quantum momenta which depend on the spatial shape of the nuclear density.The original surface hopping based on the exact factorization (SHXF) [\textit{J. Phys. Chem. Lett.} \textbf{2018}, \textit{9}, 1097] exploits frozen Gaussian functions to construct the nuclear density in the ENC term while the phase of electronic states is approximated as a fictitious nuclear momentum change.However, in reality, the width of nuclear wave packets varies in time depending on the shape of potential energy surfaces.In this work, we present a modified SHXF approach and a newly-developed Ehrenfest dynamics based on the exact factorization (EhXF) with time-dependent Gaussian functions and phases by enforcing total energy conservation.We perform numerical tests for various one-dimensional two-state model Hamiltonians.Overall, the time-dependent width of Gaussian functions and the energy conserving phase show a reliable decoherence compared to the original frozen Gaussian-based SHXF and the exact quantum mechanical calculation.Especially, the energy conserving phase is crucial for EhXF to reproduce the correct quantum dynamics.

Coupled- and Independent-Trajectory Approaches Based on the Exact Factorization Using the PyUNIxMD Package

We present mixed quantum-classical approaches based on the exact factorization framework. The electron-nuclear correlation term in the exact factorization enables us to deal with quantum coherences by accounting for electronic and nuclear nonadiabatic couplings effectively within classical nuclei approximation. We compare coupled- and independent-trajectory approximations with each other to understand algorithms in description of the bifurcation of nuclear wave packets and the correct spatial distribution of electronic wave functions along with nuclear trajectories. Finally, we show numerical results for comparisons of coupled- and independent-trajectory approaches for the photoisomerization of a protonated Schiff base from excited state molecular dynamics (ESMD) simulations with the recently developed Python-based ESMD code, namely, the PyUNIxMD program.

Coherent Two-Dimensional and Broadband Electronic Spectroscopies

Over the past few decades, coherent broadband spectroscopy has been widely used to improve our understanding of ultrafast processes (e.g., photoinduced electron transfer, proton transfer, and proton-coupled electron transfer reactions) at femtosecond resolution. The advances in femtosecond laser technology along with the development of nonlinear multidimensional spectroscopy enabled further insights into ultrafast energy transfer and carrier relaxation processes in complex biological and material systems. New discoveries and interpretations have led to improved design principles for optimizing the photophysical properties of various artificial systems. In this review, we first provide a detailed theoretical framework of both coherent broadband and two-dimensional electronic spectroscopy (2DES). We then discuss a selection of experimental approaches and considerations of 2DES along with best practices for data processing and analysis. Finally, we review several examples where coherent broadband and 2DES were employed to reveal mechanisms of photoinitiated ultrafast processes in molecular, biological, and material systems. We end the review with a brief perspective on the future of the experimental techniques themselves and their potential to answer an even greater range of scientific questions.

Femtosecond anisotropy excitation spectroscopy to disentangle the Qx and Qy absorption in chlorophyll a

Femtosecond anisotropy excitation spectroscopy is a powerful tool for unraveling contributions of electronic transitions. Here, we applied it to chlorophyll a, identifying the contribution of Qx and Qy transitions within its absorption spectrum.

Impact of energy limitations on function and resilience in long-wavelength Photosystem II

Photosystem II (PSII) uses the energy from red light to split water and reduce quinone, an energy-demanding process based on chlorophyll a (Chl-a) photochemistry. Two types of cyanobacterial PSII can use chlorophyll d (Chl-d) and chlorophyll f (Chl-f) to perform the same reactions using lower energy, far-red light. PSII from Acaryochloris marina has Chl-d replacing all but one of its 35 Chl-a, while PSII from Chroococcidiopsis thermalis, a facultative far-red species, has just 4 Chl-f and 1 Chl-d and 30 Chl-a. From bioenergetic considerations, the far-red PSII were predicted to lose photochemical efficiency and/or resilience to photodamage. Here, we compare enzyme turnover efficiency, forward electron transfer, back-reactions and photodamage in Chl-f-PSII, Chl-d-PSII, and Chl-a-PSII. We show that: (i) all types of PSII have a comparable efficiency in enzyme turnover; (ii) the modified energy gaps on the acceptor side of Chl-d-PSII favour recombination via P_D1⁺Phe^- repopulation, leading to increased singlet oxygen production and greater sensitivity to high-light damage compared to Chl-a-PSII and Chl-f-PSII; (iii) the acceptor-side energy gaps in Chl-f-PSII are tuned to avoid harmful back reactions, favouring resilience to photodamage over efficiency of light usage. The results are explained by the differences in the redox tuning of the electron transfer cofactors Phe and Q_A and in the number and layout of the chlorophylls that share the excitation energy with the primary electron donor. PSII has adapted to lower energy in two distinct ways, each appropriate for its specific environment but with different functional penalties.

Hydroquinone redox mediator enhances the photovoltaic performances of chlorophyll-based bio-inspired solar cells

Chlorophyll (Chl) derivatives have recently been proposed as photoactive materials in next-generation bio-inspired solar cells, because of their natural abundance, environmental friendliness, excellent photoelectric performance, and biodegradability. However, the intrinsic excitation dynamics of Chl derivatives remain unclear. Here, we show sub-nanosecond pump-probe time-resolved absorption spectroscopy of Chl derivatives both in solution and solid film states. We observe the formation of triplet-excited states of Chl derivatives both in deoxygenated solutions and in film samples by adding all-trans-β-carotene as a triplet scavenger. In addition, radical species of the Chl derivatives in solution were identified by adding hydroquinone as a cation radical scavenger and/or anion radical donor. These radical species (either cations or anions) can become carriers in Chl-derivative-based solar cells. Remarkably, the introduction of hydroquinone to the film samples enhanced the carrier lifetimes and the power conversion efficiency of Chl-based solar cells by 20% (from pristine 1.29% to 1.55%). This enhancement is due to a charge recombination process of Chl-A⁺/Chl-D^-, which is based on the natural Z-scheme process of photosynthesis.

PyUNIxMD: A Python-based excited state molecular dynamics package

Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.

Machine Learning-Assisted Excited State Molecular Dynamics with the State-Interaction State-Averaged Spin-Restricted Ensemble-Referenced Kohn-Sham Approach

We present a machine learning-assisted excited state molecular dynamics (ML-ESMD) based on the ensemble density functional theory framework. Since we represent a diabatic Hamiltonian in terms of generalized valence bond ansatz within the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS) method, we can avoid singularities near conical intersections, which are crucial in excited state molecular dynamics simulations. We train the diabatic Hamiltonian elements and their analytical gradients with the SchNet architecture to construct machine learning models, while the phase freedom of off-diagonal elements of the Hamiltonian is cured by introducing the phase-less loss function. Our machine learning models show reasonable accuracy with mean absolute errors of ∼0.1 kcal/mol and ∼0.5 kcal/mol/Å for the diabatic Hamiltonian elements and their gradients, respectively, for penta-2,4-dieniminium cation. Moreover, by exploiting the diabatic representation, our models can predict correct conical intersection structures and their topologies. In addition, our ML-ESMD simulations give almost identical result with a direct dynamics at the same level of theory.

Structural basis of light-harvesting in the photosystem II core complex

Photosystem II (PSII) is a membrane-spanning, multi-subunit pigment-protein complex responsible for the oxidation of water and the reduction of plastoquinone in oxygenic photosynthesis. In the present review, the recent explosive increase in available structural information about the PSII core complex based on X-ray crystallography and cryo-electron microscopy is described at a level of detail that is suitable for a future structure-based analysis of light-harvesting processes. This description includes a proposal for a consistent numbering scheme of protein-bound pigment cofactors across species. The structural survey is complemented by an overview of the state of affairs in structure-based modeling of excitation energy transfer in the PSII core complex with emphasis on electrostatic computations, optical properties of the reaction center, the assignment of long-wavelength chlorophylls, and energy trapping mechanisms.

Advances and opportunities in ultrafast X-ray crystallography and ultrafast structural optical crystallography of nuclear and electronic protein dynamics

Both nuclear and electronic dynamics contribute to protein function and need multiple and complementary techniques to reveal their ultrafast structural dynamics response. Real-space information obtained from the measurement of electron density dynamics by X-ray crystallography provides aspects of both, while the molecular physics of coherence parameters and frequency-frequency correlation needs spectroscopy methods. Ultrafast pump-probe applications of protein dynamics in crystals provide real-space information through direct X-ray crystallographic structure analysis or through structural optical crystallographic analysis. A discussion of methods of analysis using ultrafast macromolecular X-ray crystallography and ultrafast nonlinear structural optical crystallography is presented. The current and future high repetition rate capabilities provided by X-ray free electron lasers for ultrafast diffraction studies provide opportunities for optical control and optical selection of nuclear coherence which may develop to access higher frequency dynamics through improvements of sensitivity and time resolution to reveal coherence directly. Specific selection of electronic coherence requires optical probes, which can provide real-space structural information through photoselection of oriented samples and specifically in birefringent crystals. Ultrafast structural optical crystallography of photosynthetic energy transfer has been demonstrated, and the theory of two-dimensional structural optical crystallography has shown a method for accessing the structural selection of electronic coherence.

Coherent two-dimensional electronic and infrared crystallography

The two-dimensional electronic and infrared spectroscopy of oriented single crystals is sensitive to structure and point group symmetry. The third order response of crystals is generally different from measurements of isotropic solutions because each coherence path that contributes to the measured field scales to the ensemble average of the four-point correlation functions of the four field-dipole interactions involved in the respective Feynman paths. An analytical evaluation of 2D optical crystallography which depends on the crystal symmetry, laboratory orientation, and the orientation in the crystallographic frame is presented. Applying a symmetry operator in the basis of the allowed polarised radiation modes provides a method for evaluation of non-zero fourth rank tensor elements alternative to direct inspection methods. Uniaxial and biaxial systems are distinguished and the contributions to the rephasing and non-rephasing directions are evaluated for isolated and coupled oscillators. By exploiting coordinate analysis, the extension of non-linear electronic and infrared crystallography for coupled oscillators demonstrates the structural, directional, and symmetry dependent selection of coherences to the four-wave mixing signal.

Femtosecond infrared spectroscopy of chlorophyll f-containing photosystem I

Femtosecond time resolved infrared spectroscopy of far-red light grown photosystem I shows chlorophyll f contributions in light harvesting and charge separation.

Directed polaron propagation in linear polypeptides induced by intramolecular vibrations and external electric pulses

We study the propagation of α-helix polarons in a model describing the nonadiabatic interaction between an electron and a lattice of quantum mechanical oscillators at physiological temperature. We show that when excited by a subpicosecond electric pulse, as induced by experimentally observed subpicosecond charge separation, the polaron is displaced by up to hundreds of lattice sites before the electron becomes delocalised. We discuss biophysical implications of our results.

Surface Hopping Dynamics beyond Nonadiabatic Couplings for Quantum Coherence

Description of correct electron-nuclear couplings is crucial in modeling of nonadiabatic dynamics. Within traditional semiclassical or mixed quantum-classical dynamics, the coupling between quantum electronic states and classical nuclear trajectories is governed by nonadiabatic coupling vectors coupled to classical nuclear momenta. This enables us to develop a very powerful nonadiabatic dynamics algorithm, namely, surface hopping dynamics, which can describe the splitting of nuclear wave packets and detailed balance. Despite its efficiency and practicality, it suffers from the lack of quantum decoherence due to incorrect accounts for the electron-nuclear coupling. Here we present a new surface hopping algorithm based on the exact electron-nuclear correlation from the exact factorization of molecular wave functions. This algorithm demands comparable computational costs to existing surface hopping methods. Numerical simulations with two-state models and a multidimensional multistate realistic molecule show that the electron-nuclear coupling beyond the nonadiabatic coupling terms can describe the quantum coherence properly.