Origin of Low-Lying Red States in the Lhca4 Light-Harvesting Complex of Photosystem I.

The antenna complexes of Photosystem I present low-lying states visible as red-shifted and broadened absorption and fluorescence bands. Among these, Lhca4 has the most evident features of these "red" states, with a fluorescence band shifted by more than 25 nm from typical LHC emission. This signal arises from a mixing of exciton and charge-transfer (CT) states within the excitonically coupled a603-a609 chlorophyll (Chl) dimer. Here we combine molecular dynamics, multiscale quantum chemical calculations, and spectral simulations to uncover the molecular mechanism for the formation and tuning of exciton-CT interactions in Lhca4. We show that the coupling between exciton and CT states is extremely sensitive to tiny variations in the Chl dimer arrangement, explaining both the red-shifted bands and the switch between conformations with blue and red emission observed in single-molecule spectroscopy. Finally, we show that mutating the axial ligand of a603 diminishes the exciton-CT coupling, removing any red-state fingerprint.

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.

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.

Harvesting Far-Red Light with Plant Antenna Complexes Incorporating Chlorophyll d.

Increasing the absorption cross section of plants by introducing far-red absorbing chlorophylls (Chls) has been proposed as a strategy to boost crop yields. To make this strategy effective, these Chls should bind to the photosynthetic complexes without altering their functional architecture. To investigate if plant-specific antenna complexes can provide the protein scaffold to accommodate these Chls, we have reconstituted the main light-harvesting complex (LHC) of plants LHCII in vitro and in silico, with Chl d. The results demonstrate that LHCII can bind Chl d in a number of binding sites, shifting the maximum absorption ∼25 nm toward the red with respect to the wild-type complex (LHCII with Chl a and b) while maintaining the native LHC architecture. Ultrafast spectroscopic measurements show that the complex is functional in light harvesting and excitation energy transfer. Overall, we here demonstrate that it is possible to obtain plant LHCs with enhanced far-red absorption and intact functional properties.

Uncovering the interactions driving carotenoid binding in light-harvesting complexes.

Carotenoids are essential constituents of plant light-harvesting complexes (LHCs), being involved in protein stability, light harvesting, and photoprotection. Unlike chlorophylls, whose binding to LHCs is known to require coordination of the central magnesium, carotenoid binding relies on weaker intermolecular interactions (such as hydrogen bonds and van der Waals forces), whose character is far more elusive. Here we addressed the key interactions responsible for carotenoid binding to LHCs by combining molecular dynamics simulations and polarizable quantum mechanics/molecular mechanics calculations on the major LHC, LHCII. We found that carotenoid binding is mainly stabilized by van der Waals interactions with the surrounding chlorophyll macrocycles rather than by hydrogen bonds to the protein, the latter being more labile than predicted from structural data. Furthermore, the interaction network in the binding pockets is relatively insensitive to the chemical structure of the embedded carotenoid. Our results are consistent with a number of experimental data and challenge the role played by specific interactions in the assembly of pigment-protein complexes.

Understanding the Relation between Structural and Spectral Properties of Light-Harvesting Complex II.

Light-harvesting complex II (LHCII) is a pigment-protein complex present in higher plants and green algae. LHCII represents the main site of light absorption, and its role is to transfer the excitation energy toward the photosynthetic reaction centers, where primary energy conversion reactions take place. The optical properties of LHCII are known to depend on protein conformation. However, the relation between the structural and spectroscopic properties of the pigments is not fully understood yet. In this respect, previous classical molecular dynamics simulations of LHCII in a model membrane [Sci. Rep. 2015, 5, 1-10] have shown that the configuration and excitonic coupling of a chlorophyll (Chl) dimer functioning as the main terminal emitter of the complex are particularly sensitive to conformational changes. Here, we use quantum chemistry calculations to investigate in greater detail the effect of pigment-pigment interactions on the excited-state landscape. While most previous studies have used a local picture in which electrons are localized on single pigments, here we achieve a more accurate description of the Chl dimer by adopting a supramolecular picture where time-dependent density functional theory is applied to the whole system at once. Our results show that specific dimer configurations characterized by shorter inter-pigment distances can result in a sizable intensity decrease (up to 36%) of the Chl absorption bands in the visible spectral region. Such a decrease can be predicted only when accounting for Chl-Chl charge-transfer excitations, which is possible using the above-mentioned supramolecular approach. The charge-transfer character of the excitations is quantified by two types of analyses: one focusing on the composition of the excitations and the other directly on the observable total absorption intensities.

Molecular dynamics simulations in photosynthesis.

Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

Design principles of solar light harvesting in plants: Functional architecture of the monomeric antenna CP29.

In plants and green algae, light-harvesting complexes (LHCs) are a large family of chlorophyll binding proteins functioning as antennae, collecting solar photons and transferring the absorbed energy to the photosynthetic reaction centers, where light to chemical energy conversion begins. Although LHCs are all highly homologous in their structure and display a variety of common features, each complex finds a specific location and task in the energy transport. One example is CP29, which occupies a pivotal position in Photosystem II, bridging the peripheral antennae to the core. The design principles behind this specificity, however, are still unclear. Here, a synergetic approach combining steady-state and ultrafast spectroscopy, mutational analysis and structure-based exciton modeling allows uncovering the energy landscape of the chlorophylls bound to this complex. We found that, although displaying an overall highly conserved exciton structure very similar to that of other LHCs, CP29 possesses an additional terminal emitter domain. The simultaneous presence of two low energy sites facing the peripheral antennae and the core, allows CP29 to efficiently work as a conduit in the energy flux. Our results show that the LHCs share a common solid architecture but have finely tuned their structure to carry out specific functions.

Molecular Anatomy of Plant Photoprotective Switches: The Sensitivity of PsbS to the Environment, Residue by Residue.

Under strong sunlight, plants avoid photooxidation by quenching the excess absorbed energy. Quenching is triggered by PsbS, a membrane protein that is activated and deactivated by the light-dependent pH changes in the thylakoid lumen. The mechanism of action of this protein is unknown, but it was suggested that several glutamates act as pH sensors. However, the p Ka of glutamate is several pH units below the physiological values in the lumen. Thus, how can PsbS sense the pH of the lumen, and how does it respond to it? By applying a nonstandard molecular dynamics method that treats pH explicitly, we show that the lumen-exposed glutamates of PsbS have strongly shifted p Ka values and that such shifts are crucial for the pH sensitivity in physiological conditions. We also demonstrate that protonation drives a systematic unfolding of a region key for protein-protein interactions, indicating that PsbS response to pH is a functional conformational switch.

Different carotenoid conformations have distinct functions in light-harvesting regulation in plants.

To avoid photodamage plants regulate the amount of excitation energy in the membrane at the level of the light-harvesting complexes (LHCs). It has been proposed that the energy absorbed in excess is dissipated via protein conformational changes of individual LHCs. However, the exact quenching mechanism remains unclear. Here we study the mechanism of quenching in LHCs that bind a single carotenoid species and are constitutively in a dissipative conformation. Via femtosecond spectroscopy we resolve a number of carotenoid dark states, demonstrating that the carotenoid is bound to the complex in different conformations. Some of those states act as excitation energy donors for the chlorophylls, whereas others act as quenchers. Via in silico analysis we show that structural changes of carotenoids are expected in the LHC protein domains exposed to the chloroplast lumen, where acidification triggers photoprotection in vivo. We propose that structural changes of LHCs control the conformation of the carotenoids, thus permitting access to different dark states responsible for either light harvesting or photoprotection.

Coulomb couplings in solubilised light harvesting complex II (LHCII): challenging the ideal dipole approximation from TDDFT calculations.

The light harvesting complex II (LHCII), is a pigment-protein complex responsible for most of the light harvesting in plants. LHCII harvests sunlight and transfers excitation energy to the reaction centre of the photo-system, where the water oxidation process takes place. The energetics of LHCII can be modulated by means of conformational changes allowing a switch from a harvesting to a quenched state. In this state, the excitation energy is no longer transferred but converted into thermal energy to prevent photooxidation. Based on molecular dynamics simulations at the microsecond time scale, we have recently proposed that the switch between different fluorescent states can be probed by correlating shifts in the chromophore-chromophore Coulomb interactions to particular protein movements. However, these findings are based upon calculations in the ideal point dipole approximation (IDA) where the Coulomb couplings are simplified as first order dipole-dipole interactions, also assuming that the chromophore transition dipole moments lay in particular directions of space with constant moduli (FIX-IDA). In this work, we challenge this approximation using the time-dependent density functional theory (TDDFT) combined with the frozen density embedding (FDE) approach. Our aim is to establish up to which limit FIX-IDA can be applied and which chromophore types are better described under this approximation. For that purpose, we use the classical trajectories of solubilised light harvesting complex II (LHCII) we have recently reported [Liguori et al., Sci. Rep., 2015, 5, 15661] and selected three pairs of chromophores containing chlorophyll and carotenoids (Chl and Car): Chla611-Chla612, Chlb606-Chlb607 and Chla612-Lut620. Using the FDE in the Tamm-Dancoff approximation (FDEc-TDA), we show that IDA is accurate enough for predicting Chl-Chl Coulomb couplings. However, the FIX-IDA largely overestimates Chl-Car interactions mainly because the transition dipole for the Cars is not trivially oriented on the polyene chain.

Engineering a pH-Regulated Switch in the Major Light-Harvesting Complex of Plants (LHCII): Proof of Principle.

Under excess light, photosynthetic organisms employ feedback mechanisms to avoid photodamage. Photoprotection is triggered by acidification of the lumen of the photosynthetic membrane following saturation of the metabolic activity. A low pH triggers thermal dissipation of excess absorbed energy by the light-harvesting complexes (LHCs). LHCs are not able to sense pH variations, and their switch to a dissipative mode depends on stress-related proteins and allosteric cofactors. In green algae the trigger is the pigment-protein complex LHCSR3. Its C-terminus is responsible for a pH-driven conformational change from a light-harvesting to a quenched state. Here, we show that by replacing the C-terminus of the main LHC of plants with that of LHCSR3, it is possible to regulate its excited-state lifetime solely via protonation, demonstrating that the protein template of LHCs can be modified to activate reversible quenching mechanisms independent of external cofactors and triggers.

Excitation dynamics and structural implication of the stress-related complex LHCSR3 from the green alga Chlamydomonas reinhardtii.

LHCSR3 is a member of the Light-Harvesting Complexes (LHC) family, which is mainly composed of pigment-protein complexes responsible for collecting photons during the first steps of photosynthesis. Unlike related LHCs, LHCSR3 is expressed in stress conditions and has been shown to be essential for the fast component of photoprotection, non-photochemical quenching (NPQ), in the green alga Chlamydomonas reinhardtii. In plants, which do not possess LHCSR homologs, NPQ is triggered by the PSBS protein. Both PSBS and LHCSR3 possess the ability to sense pH changes but, unlike PSBS, LHCSR3 binds multiple pigments. In this work we have analyzed the properties of the pigments bound to LHCSR3 and their excited state dynamics. The data show efficient excitation energy transfer between pigments with rates similar to those observed for the other LHCs. Application of an exciton model based on a template of LHCII, the most abundant LHC, satisfactorily explains the collected steady state and time-resolved spectroscopic data, indicating that LHCSR3 has a LHC-like molecular architecture, although it probably binds less pigments. The model suggests that most of the chlorophylls have similar energy and interactions as in LHCII. The most striking difference is the localization of the lowest energy state, which is not on the Chlorophyll a (Chl a) 610-611-612 triplet as in all the LHCB antennas, but on Chl a613, which is located close to the lumen and to the pH-sensing region of the protein.

From light-harvesting to photoprotection: structural basis of the dynamic switch of the major antenna complex of plants (LHCII).

Light-Harvesting Complex II (LHCII) is largely responsible for light absorption and excitation energy transfer in plants in light-limiting conditions, while in high-light it participates in photoprotection. It is generally believed that LHCII can change its function by switching between different conformations. However, the underlying molecular picture has not been elucidated yet. The available crystal structures represent the quenched form of the complex, while solubilized LHCII has the properties of the unquenched state. To determine the structural changes involved in the switch and to identify potential quenching sites, we have explored the structural dynamics of LHCII, by performing a series of microsecond Molecular Dynamics simulations. We show that LHCII in the membrane differs substantially from the crystal and has the signatures that were experimentally associated with the light-harvesting state. Local conformational changes at the N-terminus and at the xanthophyll neoxanthin are found to strongly correlate with changes in the interactions energies of two putative quenching sites. In particular conformational disorder is observed at the terminal emitter resulting in large variations of the excitonic coupling strength of this chlorophyll pair. Our results strongly support the hypothesis that light-harvesting regulation in LHCII is coupled with structural changes.

Atomistic and Coarse Grain Topologies for the Cofactors Associated with the Photosystem II Core Complex.

Electron transfers within and between protein complexes are core processes of the electron transport chains occurring in thylakoid (chloroplast), mitochondrial, and bacterial membranes. These electron transfers involve a number of cofactors. Here we describe the derivation of molecular mechanics parameters for the cofactors associated with the function of the photosystem II core complex: plastoquinone, plastoquinol, heme b, chlorophyll A, pheophytin, and ß-carotene. Parameters were also obtained for ubiquinol and ubiquinone, related cofactors involved in the respiratory chain. Parameters were derived at both atomistic and coarse grain (CG) resolutions, compatible with the building blocks of the GROMOS united-atom and Martini CG force fields, respectively. Structural and thermodynamic properties of the cofactors were compared to experimental values when available. The topologies were further tested in molecular dynamics simulations of the cofactors in their physiological environment, e.g., either in a lipid membrane environment or in complex with the heme binding protein bacterioferritin.

Regulation of light harvesting in the green alga Chlamydomonas reinhardtii: the C-terminus of LHCSR is the knob of a dimmer switch.

Feedback mechanisms that dissipate excess photoexcitations in light-harvesting complexes (LHCs) are necessary to avoid detrimental oxidative stress in most photosynthetic eukaryotes. Here we demonstrate the unique ability of LHCSR, a stress-related LHC from the model organism Chlamydomonas reinhardtii, to sense pH variations, reversibly tuning its conformation from a light-harvesting state to a dissipative one. This conformational change is induced exclusively by the acidification of the environment, and the magnitude of quenching is correlated to the degree of acidification of the environment. We show that this ability to respond to different pH values is missing in the related major LHCII, despite high structural homology. Via mutagenesis and spectroscopic characterization, we show that LHCSR's uniqueness relies on its peculiar C-terminus subdomain, which acts as a sensor of the lumenal pH, able to tune the quenching level of the complex.

Embedding Aß42 in heterogeneous membranes depends on cholesterol asymmetries.

Using a coarse-grained lipid and peptide model, we show that the free energy stabilization of amyloid-ß in heterogeneous lipid membranes is predicted to have a dependence on asymmetric distributions of cholesterol compositions across the membrane leaflets. We find that a highly asymmetric cholesterol distribution that is depleted on the exofacial leaflet but enhanced on the cytofacial leaflet of the model lipid membrane thermodynamically favors membrane retention of a fully embedded Aß peptide. However, in the case of cholesterol redistribution that increases concentration of cholesterol on the exofacial layer, typical of aging or Alzheimer's disease, the free energy favors peptide extrusion of the highly reactive N-terminus into the extracellular space that may be vulnerable to aggregation, oligomerization, or deleterious oxidative reactivity.