Structure-Based Exciton Hamiltonian and Dynamics for the Reconstituted Wild-type CP29 Protein Antenna Complex of the Photosystem II

Probing the Effect of Mutations on Light Harvesting in CP29 by Transient Absorption and First-Principles Simulations

Natural light harvesting is exceptionally efficient thanks to the local energy funnel created within light-harvesting complexes (LHCs). To understand the design principles underlying energy transport in LHCs, ultrafast spectroscopy is often complemented by mutational studies that introduce perturbations into the excitonic structure of the natural complexes. However, such studies may fall short of identifying all excitation energy transfer (EET) pathways and their changes upon mutation. Here, we show that a synergistic combination of first-principles calculations and ultrafast spectroscopy can give unprecedented insight into the EET pathways occurring within LHCs. We measured the transient absorption spectra of the minor CP29 complex of plants and of two mutants, systematically mapping the kinetic components seen in experiments to the simulated exciton dynamics. With our combined strategy, we show that EET in CP29 is surprisingly robust to the changes in the exciton states induced by mutations, explaining the versatility of plant LHCs.

Protein Effects on the Excitation Energies and Exciton Dynamics of the CP24 Antenna Complex

In this study, the site energy fluctuations, energy transfer dynamics, and some spectroscopic properties of the minor light-harvesting complex CP24 in a membrane environment were determined. For this purpose, a 3 μs-long classical molecular dynamics simulation was performed for the CP24 complex. Furthermore, using the density functional tight binding/molecular mechanics molecular dynamics (DFTB/MM MD) approach, we performed excited state calculations for the chlorophyll a and chlorophyll b molecules in the complex starting from five different positions of the MD trajectory. During the extended simulations, we observed variations in the site energies of the different sets as a result of the fluctuating protein environment. In particular, a water coordination to Chl-b 608 occurred only after about 1 μs in the simulations, demonstrating dynamic changes in the environment of this pigment. From the classical and the DFTB/MM MD simulations, spectral densities and the (time-dependent) Hamiltonian of the complex were determined. Based on these results, three independent strongly coupled chlorophyll clusters were revealed within the complex. In addition, absorption and fluorescence spectra were determined together with the exciton relaxation dynamics, which reasonably well agrees with experimental time scales.

First-principles simulation of excitation energy transfer and transient absorption spectroscopy in the CP29 light-harvesting complex

The dynamics of delocalized excitons in light-harvesting complexes (LHCs) can be investigated using different experimental techniques, and transient absorption (TA) spectroscopy is one of the most valuable methods for this purpose. A careful interpretation of TA spectra is essential for the clarification of excitation energy transfer (EET) processes occurring during light-harvesting. However, even in the simplest LHCs, a physical model is needed to interpret transient spectra as the number of EET processes occurring at the same time is very large to be disentangled from measurements alone. Physical EET models are commonly built by fittings of the microscopic exciton Hamiltonians and exciton-vibrational parameters, an approach that can lead to biases. Here, we present a first-principles strategy to simulate EET and transient absorption spectra in LHCs, combining molecular dynamics and accurate multiscale quantum chemical calculations to obtain an independent estimate of the excitonic structure of the complex. The microscopic parameters thus obtained are then used in EET simulations to obtain the population dynamics and the related spectroscopic signature. We apply this approach to the CP29 minor antenna complex of plants for which we follow the EET dynamics and transient spectra after excitation in the chlorophyll b region. Our calculations reproduce all the main features observed in the transient absorption spectra and provide independent insight on the excited-state dynamics of CP29. The approach presented here lays the groundwork for the accurate simulation of EET and unbiased interpretation of transient spectra in multichromophoric systems.

The origin of pigment-binding differences in CP29 and LHCII: the role of protein structure and dynamics

The first step of photosynthesis in plants is performed by the light-harvesting complexes (LHC), a large family of pigment-binding proteins embedded in the photosynthetic membranes. These complexes are conserved across species, suggesting that each has a distinct role. However, they display a high degree of sequence homology and their static structures are almost identical. What are then the structural features that determine their different properties? In this work, we compared the two best-characterized LHCs of plants: LHCII and CP29. Using molecular dynamics simulations, we could rationalize the difference between them in terms of pigment-binding properties. The data also show that while the loops between the helices are very flexible, the structure of the transmembrane regions remains very similar in the crystal and the membranes. However, the small structural differences significantly affect the excitonic coupling between some pigment pairs. Finally, we analyzed in detail the structure of the long N-terminus of CP29, showing that it is structurally stable and it remains on top of the membrane even in the absence of other proteins. Although the structural changes upon phosphorylation are minor, they can explain the differences in the absorption properties of the pigments observed experimentally.

Frequency-Domain Spectroscopic Study of the Photosystem I Supercomplexes, Isolated IsiA Monomers, and the Intact IsiA Ring

The PSI₃-IsiA₁₈ supercomplex is one of the largest and most complicated assemblies in photosynthesis. The IsiA ring, composed of 18 IsiA monomers (IsiA₁₈) surrounding the PSI trimer (PSI₃), forms under iron-deficient conditions in cyanobacteria and acts as a peripheral antenna. Based on the supercomplex structure recently determined via cryo-EM imaging, we model various optical spectra of the IsiA monomers and IsiA₁₈ ring. Comparison of the absorption and emission spectra of the isolated IsiA monomers and the full ring reveals that about 2.7 chlorophylls (Chls) are lost in the isolated IsiA monomers. The best fits for isolated monomers spectra are obtained assuming the absence of Chl 508 and Chl 517 and 70% loss of Chl 511. The best model describing all three hexamers and the entire ring suggests that the lowest energy pigments are Chls 511, 514, and 517. Based on the modeling results presented in this work, we conclude that there are most likely three entry points for EET from the IsiA₆ hexamer to the PSI core monomer, with two of these entry points likely being located next to each other (i.e., nine entry points from IsiA₁₈ to the PSI₃ trimer). Finally, we show that excitation energy transfer inside individual monomers is fast (<2 ps at T = 5 K) and at least 20 times faster than intermonomer energy transfer.

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 pigments are moderately coupled (∼15 cm-1), but ctR spectra may differ significantly from exact spectra when strong interpigment coupling (&gt;100 cm-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.

Time-dependent atomistic simulations of the CP29 light-harvesting complex

Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.

Exciton Lifetime Distributions and Population Dynamics in the FMO Protein Complex from Prosthecochloris aestuarii

Significant protein rearrangement upon excitation and energy transfer in Fenna-Matthews-Olson protein of Prosthecochloris aestuarii results in a modified energy landscape, which induces more changes in pigment site energies than predicted by the "standard" hole-burning theory. The energy changes are elucidated by simulations while investigating the effects of site-dependent disorder, both static (site-energy distribution widths) and dynamic (spectral density shapes). The resulting optimized site energies and their fluctuations are consistent with relative differences observed in inhomogeneous widths calculated by recent molecular dynamic simulations. Two sets of different spectral densities reveal how their shapes affect the population dynamics and distribution of exciton lifetimes. Calculations revealed the wavelength-dependent distributions of exciton lifetimes (T ₁) in the femtosecond to picosecond time frame. We suggest that the calculated multimodal and asymmetric wavelength-dependent T ₁ distributions offer more insight into the interpretation of resonant hole-burned (HB) spectra, kinetic traces in two-dimensional (2D) electronic spectroscopy experiments, and widely used global analyses in fitting data from transient absorption experiments.

On wavelength-dependent exciton lifetime distributions in reconstituted CP29 antenna of the photosystem II and its site-directed mutants

To provide more insight into the excitonic structure and exciton lifetimes of the wild type (WT) CP29 complex of photosystem II, we measured high-resolution (low temperature) absorption, emission, and hole burned spectra for the A2 and B3 mutants, which lack chlorophylls a612 and b614 (Chls), respectively. Experimental and modeling results obtained for the WT CP29 and A2/B3 mutants provide new insight on the mutation-induced changes at the molecular level and shed more light on energy transfer dynamics. Simulations of the A2 and B3 optical spectra, using the second-order non-Markovian theory, and comparison with improved fits of WT CP29 optical spectra provide more insight into their excitonic structure, mutation induced changes, and frequency-dependent distributions of exciton lifetimes (T₁). A new Hamiltonian obtained for WT CP29 reveals that deletion of Chls a612 or b614 induces changes in the site energies of all remaining Chls. Hamiltonians obtained for A2 and B3 mutants are discussed in the context of the energy landscape of chlorophylls, excitonic structure, and transfer kinetics. Our data suggest that the lowest exciton states in A2 and B3 mutants are contributed by a611(57%), a610(17%), a615(15%) and a615(58%), a611(20%), a612(15%) Chls, respectively, although other compositions of lowest energy states are also discussed. Finally, we argue that the calculated exciton decay times are consistent with both the hole-burning and recent transient absorption measurements. Wavelength-dependent T₁ distributions offer more insight into the interpretation of kinetic traces commonly described by discrete exponentials in global analysis/global fitting of transient absorption experiments.

Structural and functional differentiation of the light-harvesting protein Lhcb4 during land plant diversification

About 475 million years ago, plants originated from an ancestral green alga and evolved first as non-vascular and later as vascular plants, becoming the primary producers of biomass on lands. During that time, the light-harvesting complex II (LHCII), responsible for sunlight absorption and excitation energy transfer to the photosystem II (PSII) core, underwent extensive differentiation. Lhcb4 is an ancestral LHCII that, in flowering plants, differentiated into up to three isoforms, Lhcb4.1, Lhcb4.2 and Lhcb4.3. The pivotal position of Lhcb4 in the PSII-LHCII supercomplex (PSII-LHCIIsc) allows functioning as linker for either S- or M-trimers of LHCII to the PSII core. The increased accumulation of Lhcb4.3 observed in PSII-LHCIIsc of plants acclimated to moderate and high light intensities induced us to investigate, whether this isoform has a preferential localization in a specific PSII-LHCIIsc conformation that might explain its light-dependent accumulation. In this work, by combining an improved method for separation of different forms of PSII-LHCIIsc from thylakoids of Pisum sativum L. grown at increasing irradiances with quantitative proteomics, we assessed that Lhcb4.3 is abundant in PSII-LHCIIsc of type C₂ S₂ , and, interestingly, similar results were found for the PsbR subunit. Phylogenetic comparative analysis on different taxa of the Viridiplantae lineage and structural modeling further pointed out to an effect of the evolution of different Lhcb4 isoforms on the light-dependent modulation of the PSII-LHCIIsc organization. This information provides new insight on the properties of the Lhcb4 and its isoforms and their role on the structure, function and regulation of PSII.

Chlorophyll a/b binding-specificity in water-soluble chlorophyll protein

We altered the chlorophyll (Chl) binding sites in various versions of water-soluble chlorophyll protein (WSCP) by amino acid exchanges to alter their preferences for either Chl a or Chl b. WSCP is ideally suited for this mutational analysis since it forms a tetrameric complex with only four identical Chl binding sites. A loop of 4-6 amino acids is responsible for Chl a versus Chl b selectivity. We show that a single amino acid exchange within this loop changes the relative Chl a/b affinities by a factor of 40. We obtained crystal structures of this WSCP variant binding either Chl a or Chl b. The Chl binding sites in these structures were compared with those in the major light-harvesting complex (LHCII) of the photosynthetic apparatus in plants to search for similar structural features involved in Chl a/b binding specificity.