Spatial and Electronic Correlations in the PE545 Light-Harvesting Complex

Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation

Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.

Fluorescence enhancement of a ligand-activated fluorescent protein induced by collective noncovalent interactions

Fluorescent proteins contain an internal chromophore constituted of amino acids or an external chromophore covalently bonded to the protein. To increase their fluorescence intensities, many research groups have attempted to mutate amino acids within or near the chromophore. Recently, a new type of fluorescent protein, called UnaG, in which the ligand binds to the protein through many noncovalent interactions was discovered. Later, a series of mutants of the UnaG protein were introduced, which include eUnaG with valine 2 mutated to leucine emitting significantly stronger fluorescence than the wild type and V2T mutant, in which valine 2 is mutated to threonine, emitting weaker fluorescence than the wild type. Interestingly, the single mutation sites of both eUnaG and V2T mutants are distant from the fluorophore, bilirubin, which renders the mechanism of such fluorescence enhancement or reduction unclear. To elucidate the origin of fluorescence intensity changes induced by the single mutations, we carried out extensive analyses on MD simulations for the original UnaG, eUnaG and V2T, and found that the bilirubin ligand bound to eUnaG is conformationally more rigid than the wild-type, particularly in the skeletal dihedral angles, possibly resulting in the increase of quantum yield through a reduction of non-radiative decay. On the other hand, the bilirubin bound to the V2T appears to be flexible than that in the UnaG. Furthermore, examining the structural correlations between the ligand and proteins, we found evidence that the bilirubin ligand is encapsulated in different environments composed of protein residues and water molecules that increase or decrease the stability of the ligand. The changed protein stability affects the mobility and confinement of water molecules captured between bilirubin and the protein. Since the flexible ligand contains multiple hydrogen bond (H-bond) donors and acceptors, the H-bonding structure and dynamics of bound water molecules are highly correlated with the rigidity of the bound ligand. Our results suggest that, to understand the fluorescence properties of protein mutants, especially the ones with noncovalently bound fluorophores with internal rotations, the interaction network among protein residues, ligand, and water molecules within the binding cavity should be investigated rather than focusing on the local structure near the fluorescing moiety. Our in-depth simulation study may offer a foundation for the design principles for engineering this new class of fluorescent proteins.

Determination of the protonation preferences of bilin pigments in cryptophyte antenna complexes

The light-harvesting mechanisms of cryptophyte antenna complexes have attracted considerable attention due to their ability to exhibit maximal photosynthetic activity under very low-light conditions and to display several colors, as well as the observation of vibronic coherent features in their two-dimensional electronic spectra. However, detailed investigations on the interplay between the protein environment and their light-harvesting properties are hampered by the uncertainty related to the protonation state of the underlying bilin pigments. Here we study the protonation preferences of four types of bilin pigments including 15,16-dihydrobiliverdin (DBV), phycoerythrobilin (PEB), phycocyanobilin (PCB) and mesobiliverdin (MBV), which are found in phycoerythrin PE545 and phycocyanin PC577, PC612, PC630 and PC645 complexes. We apply quantum chemical calculations coupled to continuum solvation calculations to predict the intrinsic acidity of bilins in aqueous solution, and then combine molecular dynamics simulations with empirical pKa estimates to investigate the impact of the local protein environment on the acidity of the pigments. We also report measurements of the absorption spectra of the five complexes in a wide range of pH in order to validate our simulations and investigate possible changes in the light harvesting properties of the complexes in the range of physiological pH found in the lumen (pH ∼ 5-7). The results suggest a pKa > 7 for DBV and MBV pigments in the α polypeptide chains of PE545 and PC630/PC645 complexes, which are not coordinated to a negatively charged amino acid. For the other PEB, DBV and PCB pigments, which interact with a Glu or Asp side chain, higher pKa values (pKa > 8) are estimated. Overall, the results support a preferential population of the fully protonated state for bilins in cryptophyte complexes under physiological conditions regardless of the specific type of pigment and local protein environment.

Exciton transport in the PE545 complex: insight from atomistic QM/MM-based quantum master equations and elastic network models

In this paper, we work out a parameterization of environmental noise within the Haken-Strobl-Reinenker (HSR) model for the PE545 light-harvesting complex, based on atomic-level quantum mechanics/molecular mechanics (QM/MM) simulations. We use this approach to investigate the role of various auto- and cross-correlations in the HSR noise tensor, confirming that site-energy autocorrelations (pure dephasing) terms dominate the noise-induced exciton mobility enhancement, followed by site energy-coupling cross-correlations for specific triplets of pigments. Interestingly, several cross-correlations of the latter kind, together with coupling-coupling cross-correlations, display clear low-frequency signatures in their spectral densities in the 30-70 [Formula: see text] region. These slow components lie at the limits of validity of the HSR approach, which requires that environmental fluctuations be faster than typical exciton transfer time scales. We show that a simple coarse-grained elastic-network-model (ENM) analysis of the PE545 protein naturally spotlights collective normal modes in this frequency range that represent specific concerted motions of the subnetwork of cysteines covalenty linked to the pigments. This analysis strongly suggests that protein scaffolds in light-harvesting complexes are able to express specific collective, low-frequency normal modes providing a fold-rooted blueprint of exciton transport pathways. We speculate that ENM-based mixed quantum classical methods, such as Ehrenfest dynamics, might be promising tools to disentangle the fundamental designing principles of these dynamical processes in natural and artificial light-harvesting structures.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Towards an ab initio description of the optical spectra of light-harvesting antennae: application to the CP29 complex of photosystem II

Light-harvesting pigment-protein complexes (PPC) represent the fundamental units through which the photosynthetic organisms absorb sunlight and funnel the energy to the reaction centre for carrying out the primary energy conversion reactions of photosynthesis. Here we apply a multiscale computational strategy to a specific PPC present in the photosystem II of plants and algae (CP29) to investigate in what detail should the environment effects due to protein and membrane/solvent be included for an accurate description of optical spectra. We find that a refinement of the crystal structure is needed before any meaningful quantum chemical calculations of pigment transition energies can be performed. For this purpose we apply classical molecular dynamics simulations of the PPC within its natural environment and we perform ab initio computations of the exciton Hamiltonian of the complex, including the environment either implicitly by the polarizable continuum model (PCM) or explicitly using the polarizable QM/MM methodology (MMPol). However, PCM essentially leads to an unspecific redshift of all transition energies, and MMPol is able to reveal site-specific changes in the optical properties of the pigments. Based on the latter and the excitonic couplings obtained within a polarizable QM/MM methodology, optical spectra are calculated, which are in good qualitative agreement with experimental data. A weakness of the approach is however found in the overestimation of the fluctuations of the excitonic parameters of the pigments along the MD trajectory. An explanation for such a finding in terms of the limits of the force fields commonly used for protein cofactors is presented and discussed.

Combining classical molecular dynamics and quantum mechanical methods for the description of electronic excitations: The case of carotenoids

Carotenoids are important actors both in light-harvesting (LH) and in photoprotection functions of photosynthetic pigment-protein complexes. A deep theoretical investigation of this multiple role is still missing owing to the difficulty of describing the delicate interplay between electronic and nuclear degrees of freedom. A possible strategy is to combine accurate quantum mechanical (QM) methods with classical molecular dynamics. To do this, however, accurate force-fields (FF) are necessary. This article presents a new FF for the different carotenoids present in LH complexes of plants. The results show that all the important structural properties described by the new FF are in very good agreement with QM reference values. This increased accuracy in the simulation of the structural fluctuations is also reflected in the description of excited states. Both the energy order and the different nature of the lowest singlet states are preserved during the dynamics when the new FF is used, whereas an unphysical mixing is found when a standard FF is used.

Developing accurate molecular mechanics force fields for conjugated molecular systems

A rapid method to parameterize the intramolecular component of classical force fields is proposed and applied to a molecular semiconductor, oligomers of conjugated polymers and a biological chromophore.

Limits and potentials of quantum chemical methods in modelling photosynthetic antennae

A critical overview of quantum chemical approaches to simulate the light-harvesting process in photosynthetic antennae is presented together with a perspective on the developments that need to be introduced to reach a quantitative predictive power.

Simulation of two-dimensional electronic spectra of phycoerythrin 545 at ambient temperature

By using a hierarchical equations-of-motion approach, we reproduce the two-dimensional electronic spectra of phycoerythrin 545 from Rhodomonas CS24 at ambient temperature (294 K). The simulated spectra are in agreement with the experimental results reported in Wong et al. (Nat. Chem. 2012, 4, 396). The evolutions of cross peaks for rephasing spectra and diagonal peaks for nonrephasing spectra have also been plotted. The peaks oscillate with the population times, with frequencies, phases, and amplitudes of the oscillating curves also being qualitatively consistent with the experimental results.

Different Types of Vibrations Interacting with Electronic Excitations in Phycoerythrin 545 and Fenna-Matthews-Olson Antenna Systems

The interest in the phycoerythrin 545 (PE545) photosynthetic antenna system of marine algae and the Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria has drastically increased since long-lived quantum coherences were reported for these complexes. For the PE545 complex, this phenomenon is clearly visible even at ambient temperatures, while for the FMO system it is more prominent at lower temperatures. The key to elucidate the role of the environment in these long-lived quantum effects is the spectral density. Here, we employ molecular dynamics simulations combined with quantum chemistry calculations to study the coupling between the biological environment and the vertical excitation energies of the bilin pigment molecules in PE545 and compare them to prior calculations on the FMO complex. It is found that the overall strength of the resulting spectral densities for the PE545 system is similar to the experiment-based counterpart but also to those in the FMO complex. Molecular analysis, however, reveals that the origin for the spectral densities in the low frequency range, which is most important for excitonic transitions, is entirely different. In the case of FMO, this part of the spectral density is due to environmental fluctuations, while, in case of PE545, it is essentially only due to internal modes of the bilin molecules. This finding sheds new light on possible explanations of the long-lived quantum coherences and that the reasons might actually be different in dissimilar systems.

Molecular basis of the exciton–phonon interactions in the PE545 light-harvesting complex

A fully polarizable QM/MM approach is used in combination with classical MD simulations to predict the pigment-dependent spectral densities of the PE545 antenna complex and account for their effects on the exciton dynamics.