Photosynthetic Light-harvesting Complexes Research Articles

Mpsqd: A matrix product state based Python package to simulate closed and open system quantum dynamics.

We introduce a Python package based on matrix product states (MPS) to simulate both the time-dependent Schrödinger equation (TDSE) and the hierarchical equations of motion (HEOM). The wave function in the TDSE or the reduced density operator/auxiliary density operators in the HEOM are represented using MPS. A matrix product operator (MPO) is then constructed to represent the Hamiltonian in the TDSE or the generalized Liouvillian in the HEOM. The fourth-order Runge-Kutta method and the time-dependent variational principle are used to propagate the MPS. Several examples, including the nonadiabatic interconversion dynamics of the pyrazine molecule, excitation energy transfer dynamics in molecular aggregates and photosynthetic light-harvesting complexes, the spin-boson model, a laser driven two-state model, the Holstein model, and charge transport in the Anderson impurity model, are presented to demonstrate the capability of the package.

The Unusual Functional Role of Protein Flexibility in Photosynthetic Light Harvesting: Protein Dynamics Studied Using Neutron Scattering

In addition to investigations of the three-dimensional protein structure, information on the dynamical properties of proteins is indispensable for an understanding of protein function in general. Correlations between protein dynamics and function are typically anticipated when both molecular mobility and function are concurrently affected under specific temperatures or hydration conditions. In contrast, excitation energy transfer within the major photosynthetic light-harvesting complex II (LHC II) presents an atypical case, as it remains fully operational even at cryogenic temperatures, primarily depending on the interactions between electronic states and involving harmonic protein vibrations only. This review summarizes recent work on vibrational and conformational protein dynamics of LHC II and directly relates these findings to its light-harvesting function. In addition, we give a comprehensive introduction into the use of neutron spectroscopy and molecular dynamics simulations to investigate the protein dynamics of photosynthetic protein complexes in solution, which is information complementary to that obtained by protein crystallography.

Simulation of the Pump-Probe Spectra and Excitation Energy Relaxation of the B850 Band of the LH2 Complex in Purple Bacteria.

Ultrafast spectroscopic techniques have been vital in studying excitation energy transfer (EET) in photosynthetic light harvesting complexes. In this paper, we simulate the pump-probe spectra of the B850 band of the light harvesting complex 2 (LH2) of purple bacteria, by using the hierarchical equation of motion method and the optical response function approach. The ground state bleach, stimulated emission, and excited state absorption components of the pump-probe spectra are analyzed in detail. The laser pulse-induced population dynamics are also simulated to help understand the main features of the pump-probe spectra and the EET process. It is shown that the excitation energy relaxation is an ultrafast process with multiple time scales. The first 40 fs of the pump-probe spectra is dominated by the relaxation of the k = ±1 states to both the k = 0 and higher energy states. Dynamics on a longer time scale around 200 fs reflects the relaxation of higher energy states to the k = 0 state.

Direct observation of triplet energy transfer between chlorophylls and carotenoids in the core antenna of photosystem I from Thermosynechococcus elongatus

Quenching of chlorophyll triplet states by carotenoids is an essential photoprotective process, which prevents formation of reactive singlet oxygen in photosynthetic light-harvesting complexes. The process is usually very efficient in oxygenic organisms under physiological conditions, thus preventing any observable accumulation of chlorophyll triplets. However, it subsequently prevents also the determination of the triplet transfer rate. Here we report results of nanosecond transient absorption spectroscopy on photosystem I core complexes, where a major part of chlorophyll a triplet states (~60 %) accumulates on a nanosecond time scale at ambient temperature. As a consequence, the triplet energy transfer could be resolved and the transfer time was determined to be about 24 ns. A smaller fraction of chlorophyll a triplet states (~40 %) is quenched with a faster rate, which could not be determined. Our analysis indicates that these chlorophylls are in direct contact with carotenoids. The overall chlorophyll triplet yield in the core antenna was estimated to be ~0.3 %, which is a value two orders of magnitude smaller than in most other photosynthetic light-harvesting complexes. This explains why slower quenching of chlorophyll triplet states is sufficient for photoprotection of photosystem I. Nevertheless, the core antenna of photosystem I represents one of only few photosynthetic complexes of oxygenic organisms in which the quenching rate of the majority of chlorophyll triplets can be directly monitored under physiological temperature.

Geometry-Independent Ultrafast Energy Transfer in Bioinspired Arrays Containing Electronically Coupled BODIPY Dimers as Energy Donors.

In photosynthetic light-harvesting complexes, strong interaction between chromophores enables efficient absorption of solar radiation and has been suggested to enable ultrafast energy funneling to the reaction center. To examine whether similar effects can be realized in synthetic systems, and to determine the mechanisms of energy transfer, we synthesized and characterized a series of bioinspired arrays containing strongly-coupled BODIPY dimers as energy donors and chlorin derivatives as energy acceptors. The BODIPY dimers feature broad absorption in the range of 500-600 nm, complementing the chlorin absorption to provide absorption across the entire visible spectrum. Ultrafast (~10 ps) energy transfer was observed from photoexcited BODIPY dyads to chlorin subunits. Surprisingly, the energy-transfer rate is nearly independent of the position where the BODIPY dimer is attached to the chlorin and of the type of connecting linker. In addition, the energy-transfer rate from BODIPY dimers to chlorin is slower than the corresponding rate in arrays containing BODIPY monomers. The lower rate, corresponding to less efficient through-bond transfer, is most likely due to weaker electronic coupling between the ground state of the chlorin acceptor and the delocalized electronic state of the BODIPY dimer, compared to the localized state of a BODIPY monomer.

Origin of Energy Dissipation in the Oligomeric Fucoxanthin-Chlorophyll a/c Binding Proteins.

Fucoxanthin-chlorophyll proteins (FCPs) are a family of photosynthetic light-harvesting complex (LHC) proteins found in diatoms. They efficiently capture photons and regulate their functions, ensuring diatom survival in highly fluctuating light. FCPs are present in different oligomeric states in vivo, but functional differences among these FCP oligomers are not yet fully understood. Here we characterized two types of antenna complexes (FCP-B/C dimers and FCP-A tetramers) that coexist in the marine centric diatom Chaetoceros gracilis using both time-resolved fluorescence and transient absorption spectroscopy. We found that the FCP-B/C complex did not show fluorescence quenching, whereas FCP-A was severely quenched, via an ultrafast excitation energy transfer (EET) pathway from Chl a Qy to the fucoxanthin S1/ICT state. These results highlight the functional differences between FCP dimers and tetramers and indicate that the EET pathway from Chl a to carotenoids is an energy dissipation mechanism conserved in a variety of photosynthetic organisms.

Resonant vibrations produce quantum bridge over high-energy states in heterogeneous antenna.

Photosynthetic light-harvesting complexes usually contain several pools of molecules with a big difference in transition energies, for example, chlorophylls a and b in plant antennas. Some pathways of the excitation energy transfer may include pigments from the low-energy pool separated by a site occupied by a high-energy molecule. We demonstrate that such pathways may be functional if high-frequency intramolecular vibrations fall in resonance with the energy gap between the neighboring molecules belonging to different pools. In this case, a vibration-assisted mixing of the excited states can produce delocalized vibronic states playing a role of 'quantum bridge' that facilitates a passage over high-energy barrier. We perform calculations of the excitation dynamics in the model three-state system with the parameters emerging from our previous studies of real antennas. Simulation of the dynamics in an explicit electron-vibrational basis demonstrates that the rate of transfer between the two chlorophylls a through the chlorophyll b intermediate is increased by a factor of 1.7-2 in the presence of resonant vibration. A possible influence of energetic disorder and other (non-resonant) vibrations on this effect is discussed.

Light-harvesting proteins intermolecular order in the Langmuir-Blodgett (LB) nanofilms – methods and applications.

Light-harvesting proteins can be used as versatile building blocks for the realization of functional biomaterials. A pathway to the hierarchical arrangement of protein molecules into functional biomaterials is provided by mesoscale (>50 nm) self-assembly. This nanotechnological approach opens new avenue in applications of light-harvesting proteins in materials science and bioelectronics (e.g., biosensors, energy-conversion devices).Phycocyanin (PC), a phycobiliprotein found in the majority of cyanobacteria and red algae’s photosynthetic light-harvesting phycobilisome complex, has the potential to be used in a wide range of nanobiotechnological applications, such as biological photovoltaic devices and as a natural photosensitizer for dye-sensitized solar cells. However, such applications frequently need for certain operating conditions, with thermal stability as one of the most crucial. Here we compare three advanced methods: Cryo-EM, X-ray nanodifraction and XFEL, previously used for probing local temperature-induced self-assembly in thermally annealed PC LB multilayered nanofilms.

Protein-Based Model for Energy Transfer between Photosynthetic Light-Harvesting Complexes Is Constructed Using a Direct Protein-Protein Conjugation Strategy.

Photosynthetic organisms utilize dynamic and complex networks of pigments bound within light-harvesting complexes to transfer solar energy from antenna complexes to reaction centers. Understanding the principles underlying the efficiency of these energy transfer processes, and how they may be incorporated into artificial light-harvesting systems, is facilitated by the construction of easily tunable model systems. We describe a protein-based model to mimic directional energy transfer between light-harvesting complexes using a circular permutant of the tobacco mosaic virus coat protein (cpTMV), which self-assembles into a 34-monomer hollow disk. Two populations of cpTMV assemblies, one labeled with donor chromophores and another labeled with acceptor chromophores, were coupled using a direct protein-protein bioconjugation method. Using potassium ferricyanide as an oxidant, assemblies containing o-aminotyrosine were activated toward the addition of assemblies containing p-aminophenylalanine. Both of these noncanonical amino acids were introduced into the cpTMV monomers through amber codon suppression. This coupling strategy has the advantages of directly, irreversibly, and site-selectively coupling donor with acceptor protein assemblies and avoids cross-reactivity with native amino acids and undesired donor-donor or acceptor-acceptor combinations. The coupled donor-acceptor model was shown to transfer energy from an antenna disk containing donor chromophores to a downstream disk containing acceptor chromophores. This model ultimately provides a controllable and modifiable platform for understanding photosynthetic interassembly energy transfer and may lead to the design of more efficient functional light-harvesting materials.

At the origin of the selectivity of the chlorophyll-binding sites in Light Harvesting Complex II (LHCII)

The photosynthetic light-harvesting complexes (LHCs) are responsible for light absorption due to their pigment-binding properties. These pigments are primarily Chlorophyll (Chl) molecules of type a and b, which ensure an excellent coverage of the visible light spectrum. To date, it is unclear which factors drive the selective binding of different Chl types in the LHC binding pockets. To gain insights into this, we employed molecular dynamics simulations on LHCII binding different Chl types. From the resulting trajectories, we have calculated the binding affinities per each Chl-binding pocket using the Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) model. To further examine the importance of the nature of the axial ligand in tuning the Chl selectivity of the binding sites, we used Density Functional Theory (DFT) calculations. The results indicate that some binding pockets have a clear Chl selectivity, and the factors governing these selectivities are identified. Other binding pockets are promiscuous, which is consistent with previous in vitro reconstitution studies. DFT calculations show that the nature of the axial ligand is not a major factor in determining the Chl binding pocket selectivity, which is instead probably controlled by the folding process.

β,β-Directly Linked Porphyrin Rings: Synthesis, Photophysical Properties, and Fullerene Binding.

Cyclic porphyrin oligomers have been studied as models for photosynthetic light-harvesting antenna complexes and as potential receptors for supramolecular chemistry. Here, we report the synthesis of unprecedented β,β-directly linked cyclic zinc porphyrin oligomers, the trimer (CP3) and tetramer (CP4), by Yamamoto coupling of a 2,3-dibromoporphyrin precursor. Their three-dimensional structures were confirmed by nuclear magnetic resonance (NMR) spectroscopy, mass spectrometry, and single-crystal X-ray diffraction analyses. The minimum-energy geometries of CP3 and CP4 have propeller and saddle shapes, respectively, as calculated using density functional theory. Their different geometries result in distinct photophysical and electrochemical properties. The smaller dihedral angles between the porphyrin units in CP3, compared with CP4, result in stronger π-conjugation, splitting the ultraviolet-vis absorption bands and shifting them to longer wavelengths. Analysis of the crystallographic bond lengths indicates that the central benzene ring of the CP3 is partially aromatic [harmonic oscillator model of aromaticity (HOMA) 0.52], whereas the central cyclooctatetraene ring of the CP4 is non-aromatic (HOMA -0.02). The saddle-shaped structure of CP4 makes it a ditopic receptor for fullerenes, with affinity constants of (1.1 ± 0.4) × 105 M-1 for C70 and (2.2 ± 0.1) × 104 M-1 for C60, respectively, in toluene solution at 298 K. The formation of a 1:2 complex with C60 is confirmed by NMR titration and single-crystal X-ray diffraction.

Cyanobacterial Phycobilisome Allostery as Revealed by Quantitative Mass Spectrometry.

Phycobilisomes (PBSs) are the major photosynthetic light-harvesting complexes in cyanobacteria and red algae. PBS, a multisubunit protein complex, has two major interfaces that comprise intrinsically disordered regions (IDRs): rod-core and core-membrane. IDRs do not form regular, three-dimensional structures on their own. Their presence in the photosynthetic pigment-protein complexes portends their structural and functional importance. A recent model suggests that PB-loop, an IDR located on the PBS subunit ApcE and C-terminal extension (CTE) of the PBS subunit ApcG, forms a structural protrusion on the PBS core-membrane side, facing the thylakoid membrane. Here, the structural synergy between the rod-core region and the core-membrane region was investigated using quantitative mass spectrometry (MS). The AlphaFold-predicted CpcG-CTE structure was first modeled onto the PBS rod-core region, guided and justified by the isotopically encoded structural MS data. Quantitative cross-linking MS analysis revealed that the structural proximity of the PB-loop in ApcE and ApcG-CTE is significantly disturbed in the absence of six PBS rods, which are attached to PBS via CpcG-CTE, indicative of drastic conformational changes and decreased structural integrity. These results suggest that CpcG-rod attachment on the PBS rod-core side is essentially required for the PBS core-membrane structural assembly. The hypothesized long-range synergy between the rod-core interface (where the orange carotenoid protein also functions) and the terminal energy emitter of PBS must have important regulatory roles in PBS core assembly, light-harvesting, and excitation energy transmission. These data also lend strategies that genetic truncation of the light-harvesting antennas aimed for improved photosynthetic productivity must rely on an in-depth understanding of their global structural integrity.

Charge transfer as a mechanism for chlorophyll fluorescence concentration quenching

Highly concentrated solutions of chlorophyll display rapid fluorescence quenching. The same devastating energy loss is not seen in photosynthetic light-harvesting antenna complexes, despite the need for chromophores to be in close proximity to facilitate energy transfer. A promising, though unconfirmed mechanism for the observed quenching is energy transfer from an excited chlorophyll monomer to a closely associated chlorophyll pair that subsequently undergoes rapid nonradiative decay to the ground state via a short-lived intermediate charge-transfer state. In this work, we make use of newly emerging fast methods in quantum chemistry to assess the feasibility of this proposed mechanism. We calculate rate constants for the initial charge separation, based on Marcus free-energy surfaces extracted from molecular dynamics simulations of solvated chlorophyll pairs, demonstrating that this pathway will compete with fluorescence (i.e., drive quenching) at experimentally measured quenching concentrations. We show that the rate of charge separation is highly sensitive to interchlorophyll distance and the relative orientations of chromophores within a quenching pair. We discuss possible solvent effects on the rate of charge separation (and consequently the degree of quenching), using the light-harvesting complex II (LH2) protein from rps. acidophila as a specific example of how this process might be controlled in a protein environment. Crucially, we reveal that the LH2 antenna protein prevents quenching, even at the high chlorophyll concentrations required for efficient energy transfer, by restricting the range of orientations that neighboring chlorophyll pairs can adopt.

Excitation dynamics in photosynthetic light-harvesting complex B850: exact solution versus Redfield and Förster limits.

We model the excitation dynamics in bacterial B850 antenna using the standard Redfield, modified Redfield, and Förster approaches and comparing them with the exact solution obtained with hierarchical equation of motion (HEOM). We have found that the modified Redfield is capable of reproducing the dynamics associated with downhill relaxation from higher exciton levels, but fails to explain the migration of quasi-equilibrated excitation over the B850 ring (and its spectral signatures like transient anisotropy decay). Neglecting the population-to-coherence transfers leads to a quick decoherence between the exciton states resulting in unrealistically fast delocalization looking like instantaneous transfer around the ring. The standard (non-secular) Redfield gives a more satisfactory picture of this kind of migration, but in some cases the results can be corrupted by artifacts emerging from the one-phonon character of this theory.

Cryo-electron microscopy structure of the intact photosynthetic light-harvesting antenna-reaction center complex from a green sulfur bacterium.

The photosynthetic reaction center complex (RCC) of green sulfur bacteria (GSB) consists of the membrane-imbedded RC core and the peripheric energy transmitting proteins called Fenna-Matthews-Olson (FMO). Functionally, FMO transfers the absorbed energy from a huge peripheral light-harvesting antenna named chlorosome to the RC core where charge separation occurs. In vivo, one RC was found to bind two FMOs, however, the intact structure of RCC as well as the energy transfer mechanism within RCC remain to be clarified. Here we report a structure of intact RCC which contains a RC core and two FMO trimers from a thermophilic green sulfur bacterium Chlorobaculum tepidum at 2.9 Å resolution by cryo-electron microscopy. The second FMO trimer is attached at the cytoplasmic side asymmetrically relative to the first FMO trimer reported previously. We also observed two new subunits (PscE and PscF) and the N-terminal transmembrane domain of a cytochrome-containing subunit (PscC) in the structure. These two novel subunits possibly function to facilitate the binding of FMOs to RC core and to stabilize the whole complex. A new bacteriochlorophyll (numbered as 816) was identified at the interspace between PscF and PscA-1, causing an asymmetrical energy transfer from the two FMO trimers to RC core. Based on the structure, we propose an energy transfer network within this photosynthetic apparatus.

Tight inner ring architecture and quantum motion of nuclei enable efficient energy transfer in bacterial light harvesting.

The efficient, directional transfer of absorbed solar energy between photosynthetic light-harvesting complexes continues to pose intriguing questions. In this work, we identify the pathways of energy flow between the B800 and B850 rings in the LH2 complex of Rhodopseudomonas molischianum using fully quantum mechanical path integral methods to simulate the excited-state dynamics of the 24 bacteriochlorophyll molecules and their coupling to 50 normal mode vibrations in each chromophore. While all pigments are identical, the tighter packing of the inner B850 ring is responsible for the thermodynamic stabilization of the inner ring. Molecular vibrations enable the 1-ps flow of energy to the B850 states, which would otherwise be kinetically inaccessible. A classical treatment of the vibrations leads to uniform equilibrium distribution of the excitation, with only 67% transferred to the inner ring. However, spontaneous fluctuations associated with the quantum motion of the nuclei increase the transfer efficiency to 90%.

Static Disorder has Dynamic Impact on Energy Transport in Biomimetic Light-Harvesting Complexes.

Despite extensive studies, many questions remain about what structural and energetic factors give rise to the remarkable energy transport efficiency of photosynthetic light-harvesting protein complexes, owing largely to the inability to synthetically control such factors in these natural systems. Herein, we demonstrate energy transfer within a biomimetic light-harvesting complex consisting of identical chromophores attached in a circular array to a protein scaffold derived from the tobacco mosaic virus coat protein. We confirm the capability of energy transport by observing ultrafast depolarization in transient absorption anisotropy measurements and a redshift in time-resolved emission spectra in these complexes. Modeling the system with kinetic Monte Carlo simulations recapitulates the observed anisotropy decays, suggesting an inter-site hopping rate as high as 1.6 ps–1. With these simulations, we identify static disorder in orientation, site energy, and degree of coupling as key remaining factors to control to achieve long-range energy transfer in these systems. We thereby establish this system as a highly promising, bottom-up model for studying long-range energy transfer in light-harvesting protein complexes.

Accuracy of approximate methods for the calculation of absorption-type linear spectra with a complex system-bath coupling.

The accuracy of approximate methods for calculating linear optical spectra depends on many variables. In this study, we fix most of these parameters to typical values found in photosynthetic light-harvesting complexes of plants and determine the accuracy of approximate spectra with respect to exact calculation as a function of the energy gap and interpigment coupling in a pigment dimer. We use a spectral density with the first eight intramolecular modes of chlorophyll a and include inhomogeneous disorder for the calculation of spectra. We compare the accuracy of absorption, linear dichroism, and circular dichroism spectra calculated using the Full Cumulant Expansion (FCE), coherent time-dependent Redfield (ctR), and time-independent Redfield and modified Redfield methods. As a reference, we use spectra calculated with the exact stochastic path integral evaluation method. We find the FCE method to be the most accurate for the calculation of all spectra. The ctR method performs well for the qualitative calculation of absorption and linear dichroism spectra when the pigments are moderately coupled (∼15cm-1), but ctR spectra may differ significantly from exact spectra when strong interpigment coupling (>100cm-1) is present. The dependence of the quality of Redfield and modified Redfield spectra on molecular parameters is similar, and these methods almost always perform worse than ctR, especially when the interpigment coupling is strong or the excitonic energy gap is small (for a given coupling). The accuracy of approximate spectra is not affected by resonance with intramolecular modes for typical system-bath coupling and disorder values found in plant light-harvesting complexes.

Machine learning for excitation energy transfer dynamics

A well-known approach to describe the dynamics of an open quantum system is to compute the master equation evolving the reduced density matrix of the system. This approach plays an important role in describing excitation transfer through photosynthetic light harvesting complexes (LHCs). The hierarchical equations of motion (HEOM) was adapted by Ishizaki and Fleming (J. Chem. Phys., 2009) to simulate open quantum dynamics in the biological regime. We generate a set of time dependent observables that depict the coherent propagation of electronic excitations through the LHCs by solving the HEOM. We solve the inverse problem using classical machine learning (ML) models as this is a computationally intractable problem. The objective here is to determine whether a trained ML model can perform Hamiltonian tomography by using the time dependence of the observables as inputs. We demonstrate the capability of convolutional neural networks to tackle this research problem. The models developed here can predict Hamiltonian parameters such as excited state energies and inter-site couplings of a system up to 99.28\% accuracy and mean-squared error as low as 0.65.

Structural dynamics of light harvesting proteins, photosynthetic membranes, and cells observed by spectral editing solid-state NMR.

Photosynthetic light-harvesting complexes have a remarkable capacity to perform robust photo-physics at ambient temperatures and in fluctuating environments. Protein conformational dynamics and membrane mobility are processes that contribute to the light-harvesting efficiencies and control photoprotective responses. This short review describes the application of magic angle spinning nuclear magnetic resonance (NMR) spectroscopy for characterizing the structural dynamics of pigment, protein, and thylakoid membrane components related to light harvesting and photoprotection. I will discuss the use of dynamics-based spectral editing solid-state NMR for distinguishing rigid and mobile components and assessing protein, pigment, and lipid dynamics on sub-nanosecond to millisecond timescales. Dynamic spectral editing NMR has been applied to investigate light-harvesting complex II protein conformational dynamics inside lipid bilayers and in native membranes. Furthermore, we used the NMR approach to assess thylakoid membrane dynamics. Finally, it is shown that dynamics-based spectral editing NMR for reducing spectral complexity by filtering motion-dependent signals enabled us to follow processes in live photosynthetic cells.