A protein dynamics study of photosystem II: the effects of protein conformation on reaction center function

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Phase Transition and Phase Separation in Realistic Thylakoid Lipid Membrane of Marine Algae in All-Atom Simulations

Thylakoid membranes are specialized membranes predominantly composed of uncommon galacto- and sulfolipids, having distinct roles in photosynthesis. Large acyl chain variety and richness in polyunsaturated fatty acid (PUFA) content of thylakoid lipids further add to the compositional complexity. The function of these membrane systems is intimately dependent on the fluidity of its lipid matrix, which is strongly modulated by the lipid composition and temperature. The present work, employing extensive atomistic simulations, provides the first atomistic view of the phase transition and domain coexistence in a model membrane composed of thylakoid lipids of a commercially important red alga Gracilaria corticata between 10 and 40 °C. The growth and photosynthetic activity of marine algae are greatly influenced by the seawater temperature. So far, little is known about the molecular organization of lipids in thylakoid membranes, in particular their adaptive arrangements under temperature stress. Our simulations show that the algal thylakoid membrane undergoes a transition from a gel-like phase at a low temperature, 10-15 °C, to a homogeneous liquid-crystalline phase at a high temperature, 40 °C. Clear evidence of spontaneous phase separation into coexisting nanoscale domains is detected at intermediate temperatures nearing the optimal growth temperature range. Particularly, at 25-30 °C, we identified the formation of a stable ripple phase, where the gel-like domains rich in saturated and nearly hexagonally packed lipids were separated from fluid-like domains enriched in lipids containing PUFA chains. The phase separation is driven by the spontaneous and preferential segregation of lipids into differentially ordered domains, mainly depending on the acyl chain types. Cholesterol impairs the phase transition and the emergence of domains and induces a fairly uniform liquid-ordered phase in the membrane over the temperatures studied. This work improves the understanding of the properties and reorganization of lipids in the thylakoid membrane in response to temperature variation.

Substitution of Ca2+ and changes in the H-bond network near the oxygen-evolving complex of photosystem II

Ca²⁺, which provides binding sites for ligand water molecules W3 and W4 in the Mn₄CaO₅ cluster, is a prerequisite for O₂ evolution in photosystem II (PSII). We report structural changes in the H-bond network and the catalytic cluster itself upon the replacement of Ca²⁺ with other alkaline earth metals, using a quantum mechanical/molecular mechanical approach. The small radius of Mg²⁺ makes W3 donate an H-bond to D1-Glu189 in Mg²⁺-PSII. If an additional water molecule binds at the large surface of Ba²⁺, it donates H-bonds to D1-Glu189 and the ligand water molecule at the dangling Mn, altering the H-bond network. The potential energy profiles of the H-bond between D1-Tyr161 (TyrZ) and D1-His190 and the interconversion between the open- and closed-cubane S₂ conformations remain substantially unaltered upon the replacement of Ca²⁺. Remarkably, the O5⋯Ca²⁺ distance is shortest among all O5⋯metal distances irrespective of the radius being larger than that of Mg²⁺. Furthermore, Ca²⁺ is the only alkaline earth metal that equalizes the O5⋯metal and O2⋯metal distances and facilitates the formation of the symmetric cubane structure.

Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II

Besides absorbing light, the core antenna complex CP43 of photosystem II is of great importance in transferring excitation energy from the antenna complexes to the reaction center. Excitation energies, spectral densities, and linear absorption spectra of the complex have been evaluated by a multiscale approach. In this scheme, quantum mechanics/molecular mechanics molecular dynamics simulations are performed employing the parameterized density functional tight binding (DFTB) while the time-dependent long-range-corrected DFTB scheme is applied for the excited state calculations. The obtained average spectral density of the CP43 complex shows a very good agreement with experimental results. Moreover, the excitonic Hamiltonian of the system along with the computed site-dependent spectral densities was used to determine the linear absorption. While a Redfield-like approximation has severe shortcomings in dealing with the CP43 complex due to quasi-degenerate states, the non-Markovian full second-order cumulant expansion formalism is able to overcome the drawbacks. Linear absorption spectra were obtained, which show a good agreement with the experimental counterparts at different temperatures. This study once more emphasizes that by combining diverse techniques from the areas of molecular dynamics simulations, quantum chemistry, and open quantum systems, it is possible to obtain first-principle results for photosynthetic complexes, which are in accord with experimental findings.

Structural and dynamic insights into Mn4Ca cluster-depleted Photosystem II

In the first steps of natural oxygenic photosynthesis, sunlight is used to oxidize water molecules to protons, electrons and molecular oxygen. This reaction takes place on the Mn₄Ca cluster located in the reaction centre of Photosystem II (PSII), where the cluster is assembled and continuously repaired through a process known as photoactivation. Understanding the molecular details of such a process has important implications in different fields, in particular inspiring synthesis and repair strategies for artificial photosynthesis devices. In this regard, a detailed structural and dynamic characterization of Photosystem II lacking a Mn₄Ca cluster, namely apo PSII, is a prerequisite for the full comprehension of the photoactivation. Recently, the structure of the apo PSII was resolved at 2.55 Å resolution [Zhang et al., eLife, 2017, 6, e26933], suggesting a pre-organized structure of the protein cavity hosting the cluster. Anyway, the question of whether these findings are a feature of the method used remains open. Here, by means of classical Molecular Dynamics simulations, we characterized the structural and dynamic features of the apo PSII for different protonation states of the cluster cavity. Albeit an overall conformational stability common to all investigated systems, we found significant deviations in the conformation of the side chains of the active site with respect to the X-ray positions. Our findings suggest that not all residues acting as Mn ligands are pre-organized prior to the Mn₄Ca formation and previous local conformational changes are required in order to bind the first Mn ion in the high-affinity binding site.

Chlorophyll excitation energies and structural stability of the CP47 antenna of photosystem II: a case study in the first-principles simulation of light-harvesting complexes

Natural photosynthesis relies on light harvesting and excitation energy transfer by specialized pigment-protein complexes. Their structure and the electronic properties of the embedded chromophores define the mechanisms of energy transfer. An important example of a pigment-protein complex is CP47, one of the integral antennae of the oxygen-evolving photosystem II (PSII) that is responsible for efficient excitation energy transfer to the PSII reaction center. The charge-transfer excitation induced among coupled reaction center chromophores resolves into charge separation that initiates the electron transfer cascade driving oxygenic photosynthesis. Mapping the distribution of site energies among the 16 chlorophyll molecules of CP47 is essential for understanding excitation energy transfer and overall antenna function. In this work, we demonstrate a multiscale quantum mechanics/molecular mechanics (QM/MM) approach utilizing full time-dependent density functional theory with modern range-separated functionals to compute for the first time the excitation energies of all CP47 chlorophylls in a complete membrane-embedded cyanobacterial PSII dimer. The results quantify the electrostatic effect of the protein on the site energies of CP47 chlorophylls, providing a high-level quantum chemical excitation profile of CP47 within a complete computational model of "near-native" cyanobacterial PSII. The ranking of site energies and the identity of the most red-shifted chlorophylls (B3, followed by B1) differ from previous hypotheses in the literature and provide an alternative basis for evaluating past approaches and semiempirically fitted sets. Given that a lot of experimental studies on CP47 and other light-harvesting complexes utilize extracted samples, we employ molecular dynamics simulations of isolated CP47 to identify which parts of the polypeptide are most destabilized and which pigments are most perturbed when the antenna complex is extracted from PSII. We demonstrate that large parts of the isolated complex rapidly refold to non-native conformations and that certain pigments (such as chlorophyll B1 and β-carotene h1) are so destabilized that they are probably lost upon extraction of CP47 from PSII. The results suggest that the properties of isolated CP47 are not representative of the native complexed antenna. The insights obtained from CP47 are generalizable, with important implications for the information content of experimental studies on biological light-harvesting antenna systems.

Protein Matrix Control of Reaction Center Excitation in Photosystem II

Photosystem II (PSII) is a multisubunit pigment-protein complex that uses light-induced charge separation to power oxygenic photosynthesis. Its reaction center chromophores, where the charge transfer cascade is initiated, are arranged symmetrically along the D1 and D2 core polypeptides and comprise four chlorophyll (PD1, PD2, ChlD1, ChlD2) and two pheophytin molecules (PheoD1 and PheoD2). Evolution favored productive electron transfer only via the D1 branch, with the precise nature of primary excitation and the factors that control asymmetric charge transfer remaining under investigation. Here we present a detailed atomistic description for both. We combine large-scale simulations of membrane-embedded PSII with high-level quantum-mechanics/molecular-mechanics (QM/MM) calculations of individual and coupled reaction center chromophores to describe reaction center excited states. We employ both range-separated time-dependent density functional theory and the recently developed domain based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD), the first coupled cluster QM/MM calculations of the reaction center. We find that the protein matrix is exclusively responsible for both transverse (chlorophylls versus pheophytins) and lateral (D1 versus D2 branch) excitation asymmetry, making ChlD1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the ChlD1 → PheoD1 charge-transfer the lowest energy excitation globally within the reaction center, lower than any pigment-centered local excitation. Remarkably, no low-energy charge transfer states are located within the "special pair" PD1-PD2, which is therefore excluded as the site of initial charge separation in PSII. Finally, molecular dynamics simulations suggest that modulation of the electrostatic environment due to protein conformational flexibility enables direct excitation of low-lying charge transfer states by far-red light.

A TDDFT investigation of the Photosystem II reaction center: Insights into the precursors to charge separation

Photosystem II (PS II) captures solar energy and directs charge separation (CS) across the thylakoid membrane during photosynthesis. The highly oxidizing, charge-separated state generated within its reaction center (RC) drives water oxidation. Spectroscopic studies on PS II RCs are difficult to interpret due to large spectral congestion, necessitating modeling to elucidate key spectral features. Herein, we present results from time-dependent density functional theory (TDDFT) calculations on the largest PS II RC model reported to date. This model explicitly includes six RC chromophores and both the chlorin phytol chains and the amino acid residues <6 Å from the pigments' porphyrin ring centers. Comparing our wild-type model results with calculations on mutant D1-His-198-Ala and D2-His-197-Ala RCs, our simulated absorption-difference spectra reproduce experimentally observed shifts in known chlorophyll absorption bands, demonstrating the predictive capabilities of this model. We find that inclusion of both nearby residues and phytol chains is necessary to reproduce this behavior. Our calculations provide a unique opportunity to observe the molecular orbitals that contribute to the excited states that are precursors to CS. Strikingly, we observe two high oscillator strength, low-lying states, in which molecular orbitals are delocalized over ChlD1 and PheD1 as well as one weaker oscillator strength state with molecular orbitals delocalized over the P chlorophylls. Both these configurations are a match for previously identified exciton-charge transfer states (ChlD1+PheD1-)* and (PD2+PD1-)*. Our results demonstrate the power of TDDFT as a tool, for studies of natural photosynthesis, or indeed future studies of artificial photosynthetic complexes.

Acquirement of water-splitting ability and alteration of the charge-separation mechanism in photosynthetic reaction centers

In photosynthetic reaction centers from purple bacteria (PbRC) and the water-oxidizing enzyme, photosystem II (PSII), charge separation occurs along one of the two symmetrical electron-transfer branches. Here we report the microscopic origin of the unidirectional charge separation, fully considering electron-hole interaction, electronic coupling of the pigments, and electrostatic interaction with the polarizable entire protein environments. The electronic coupling between the pair of bacteriochlorophylls is large in PbRC, forming a delocalized excited state with the lowest excitation energy (i.e., the special pair). The charge-separated state in the active branch is stabilized by uncharged polar residues in the transmembrane region and charged residues on the cytochrome c ₂ binding surface. In contrast, the accessory chlorophyll in the D1 protein (Chl_D1) has the lowest excitation energy in PSII. The charge-separated state involves Chl_D1 ^•+ and is stabilized predominantly by charged residues near the Mn₄CaO₅ cluster and the proceeding proton-transfer pathway. It seems likely that the acquirement of water-splitting ability makes Chl_D1 the initial electron donor in PSII.

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.

The microfluidic laboratory at Synchrotron SOLEIL

A microfluidic laboratory recently opened at Synchrotron SOLEIL, dedicated to in-house research and external users. Its purpose is to provide the equipment and expertise that allow the development of microfluidic systems adapted to the beamlines of SOLEIL as well as other light sources. Such systems can be used to continuously deliver a liquid sample under a photon beam, keep a solid sample in a liquid environment or provide a means to track a chemical reaction in a time-resolved manner. The laboratory provides all the amenities required for the design and preparation of soft-lithography microfluidic chips compatible with synchrotron-based experiments. Three examples of microfluidic systems that were used on SOLEIL beamlines are presented, which allow the use of X-ray techniques to study physical, chemical or biological phenomena.

Triplet-triplet energy transfer in artificial and natural photosynthetic antennas

In photosynthetic organisms, protection against photooxidative stress due to singlet oxygen is provided by carotenoid molecules, which quench chlorophyll triplet species before they can sensitize singlet oxygen formation. In anoxygenic photosynthetic organisms, in which exposure to oxygen is low, chlorophyll-to-carotenoid triplet-triplet energy transfer (T-TET) is slow, in the tens of nanoseconds range, whereas it is ultrafast in the oxygen-rich chloroplasts of oxygen-evolving photosynthetic organisms. To better understand the structural features and resulting electronic coupling that leads to T-TET dynamics adapted to ambient oxygen activity, we have carried out experimental and theoretical studies of two isomeric carotenoporphyrin molecular dyads having different conformations and therefore different interchromophore electronic interactions. This pair of dyads reproduces the characteristics of fast and slow T-TET, including a resonance Raman-based spectroscopic marker of strong electronic coupling and fast T-TET that has been observed in photosynthesis. As identified by density functional theory (DFT) calculations, the spectroscopic marker associated with fast T-TET is due primarily to a geometrical perturbation of the carotenoid backbone in the triplet state induced by the interchromophore interaction. This is also the case for the natural systems, as demonstrated by the hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of light-harvesting proteins from oxygenic (LHCII) and anoxygenic organisms (LH2). Both DFT and electron paramagnetic resonance (EPR) analyses further indicate that, upon T-TET, the triplet wave function is localized on the carotenoid in both dyads.

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.

Interaction of the PsbH subunit with a chlorophyll bound to histidine 114 of CP47 is responsible for the red 77K fluorescence of Photosystem II

A characteristic feature of the active Photosystem II (PSII) complex is a red-shifted low temperature fluorescence emission at about 693nm. The origin of this emission has been attributed to a monomeric 'red' chlorophyll molecule located in the CP47 subunit. However, the identity and function of this chlorophyll remain uncertain. In our previous work, we could not detect the red PSII emission in a mutant of the cyanobacterium Synechocystis sp. PCC 6803 lacking PsbH, a small transmembrane subunit bound to CP47. However, it has not been clear whether the PsbH is structurally essential for the red emission or the observed effect of mutation has been indirectly caused by compromised PSII stability and function. In the present work we performed a detailed spectroscopic characterization of PSII in cells of a mutant lacking PsbH and Photosystem I and we also characterized PSII core complexes isolated from this mutant. In addition, we purified and characterized the CP47 assembly modules containing and lacking PsbH. The results clearly confirm an essential role of PsbH in the origin of the PSII red emission and also demonstrate that PsbH stabilizes the binding of one β-carotene molecule in PSII. Crystal structures of the cyanobacterial PSII show that PsbH directly interacts with a single monomeric chlorophyll ligated by the histidine 114 residue of CP47 and we conclude that this peripheral chlorophyll hydrogen-bonded to PsbH is responsible for the red fluorescence state of CP47. Given the proximity of β-carotene this state could participate in the dissipation of excessive light energy.

Characterization of thylakoid lipid membranes from cyanobacteria and higher plants by molecular dynamics simulations

The thylakoid membrane is mainly composed of non-common lipids, so called galactolipids. Despite the importance of these lipids for the function of the photosynthetic reaction centers, the molecular organization of these membranes is largely unexplored. Here we use multiscale molecular dynamics simulations to characterize the thylakoid membrane of both cyanobacteria and higher plants. We consider mixtures of up to five different galactolipids plus phosphatidylglycerol to represent these complex membranes. We find that the different lipids generally mix well, although nanoscale heterogeneities are observed especially in case of the plant membrane. The fluidity of the cyanobacterial membrane is markedly reduced compared to the plant membrane, even considering elevated temperatures at which thermophilic cyanobacteria are found. We also find that the plant membrane more readily undergoes a phase transformation to an inverted hexagonal phase. We furthermore characterized the conformation and dynamics of the cofactors plastoquinone and plastoquinol, revealing of the fast flip-flop rates for the non-reduced form. Together, our results provide a molecular view on the dynamical organization of the thylakoid membrane.

Proton-Coupled Electron Transfer During the S-State Transitions of the Oxygen-Evolving Complex of Photosystem II

The oxygen-evolving complex (OEC) of photosystem II (PSII) is a unique Mn4O5Ca cluster that catalyzes water oxidation via four photoactivated electron transfer steps. As the protein influence on the redox and protonation chemistry of the OEC remains an open question, we present a classical valence model of the OEC that allows the redox state of each Mn and the protonation state of bridging μ-oxos and terminal waters to remain in equilibrium with the PSII protein throughout the redox cycle. We find that the last bridging oxygen loses its proton during the transition from S0 to S1. Two possible S2 states are found depending on the OEC geometry: S2 has Mn4(IV) with a proton lost from a terminal water (W1) trapped by the nearby D1-D61 if O5 is closer to Mn4, or Mn1(IV), with partial deprotonation of D1-H337 and D1-E329 if O5 is closer to Mn1. In S3, the OEC is Mn4(IV) with W2 deprotonated. The estimated OEC Em's range from +0.7 to +1.3 V, enabling oxidation by P680(+), the primary electron donor in PSII. In chloride-depleted PSII, the proton release increases during the S1 to S2 transition, leaving the OEC unable to properly advance through the water-splitting cycle.

Dynamic protein conformations preferentially drive energy transfer along the active chain of the photosystem II reaction centre

One longstanding puzzle concerning photosystem II, a core component of photosynthesis, is that only one of the two symmetric branches in its reaction centre is active in electron transfer. To investigate the effect of the photosystem II environment on the preferential selection of the energy transfer pathway (a prerequisite for electron transfer), we have constructed an exciton model via extensive molecular dynamics simulations and quantum mechanics/molecular mechanics calculations based on a recent X-ray structure. Our results suggest that it is essential to take into account an ensemble of protein conformations to accurately compute the site energies. We identify the cofactor CLA606 of active chain as the most probable site for the energy excitation. We further pinpoint a number of charged protein residues that collectively lower the CLA606 site energy. Our work provides insights into the understanding of molecular mechanisms of the core machinery of the green-plant photosynthesis.

Cofactor-mediated energy and electron transfer in photosystem II occurs preferentially through one branch of the reaction centre, despite there being a symmetric path available. Here, the authors use computational methods to determine the influence of protein conformation on this selectivity.

Molecular dynamics simulations reveal highly permeable oxygen exit channels shared with water uptake channels in photosystem II

Photosystem II (PSII) catalyzes the oxidation of water in the conversion of light energy into chemical energy in photosynthesis. Water delivery and oxygen removal from the oxygen evolving complex (OEC), buried deep within PSII, are critical requirements to facilitate the reaction and minimize reactive oxygen damage. It has often been assumed that water and oxygen travel through separate channels within PSII, as demonstrated in cytochrome c oxidase. This study describes all-atom molecular dynamics simulations of PSII designed to investigate channels by fully characterizing the distribution and permeation of both water and oxygen. Interestingly, most channels found in PSII were permeable to both oxygen and water, however individual channels exhibited different energetic barriers for the two solutes. Several routes for oxygen diffusion within PSII with low energy permeation barriers were found, ensuring its fast removal from the OEC. In contrast, all routes for water showed significant energy barriers, corresponding to a much slower permeation rate for water through PSII. Two major factors were responsible for this selectivity: (1) hydrogen bonds between water and channel amino acids, and (2) steric restraints. Our results reveal the presence of a shared network of channels in PSII optimized to both facilitate the quick removal of oxygen and effectively restrict the water supply to the OEC to help stabilize and protect it from small water soluble inhibitors.

Design and application of anthracene derivative with aggregation-induced emission charateristics for visualization and monitoring of erythropoietin unfolding

Erythropoietin (EPO) is an attractive protein-unfolding/folding model because of its high degree of unfolding and folding reversibility and intermediate size. Due to its function for regulating red blood cell production by stimulating late erythroid precursor cells, EPO presents obvious values to biological research. A nonemissive anthracene derivative, that is 9,10-bis[4-(3-sulfonatopropoxyl)-styryl]anthracene sodium salt (BSPSA), with aggregation-induced emission (AIE) charateristics shows a novel phenomenon of AIE when EPO is added. The AIE biosensor for EPO shows the limit of detection is 1 × 10(-9) M. Utilizing the AIE feature of BSPSA, the unfolding process of EPO using guanidine hydrochloride is monitored, which indicates three steps for the folding structures of EPO to transform to random coil. Computational modeling suggests that the BSPSA luminogens prefer docking in the hydrophobic cavity in the EPO folding structures, and the assembly of BSPSA in this cavity makes the AIE available, making the monitoring of unfolding of EPO possible.

Extended protein/water H-bond networks in photosynthetic water oxidation

Oxidation of water molecules in the photosystem II (PSII) protein complex proceeds at the manganese-calcium complex, which is buried deeply in the lumenal part of PSII. Understanding the PSII function requires knowledge of the intricate coupling between the water-oxidation chemistry and the dynamic proton management by the PSII protein matrix. Here we assess the structural basis for long-distance proton transfer in the interior of PSII and for proton management at its surface. Using the recent high-resolution crystal structure of PSII, we investigate prominent hydrogen-bonded networks of the lumenal side of PSII. This analysis leads to the identification of clusters of polar groups and hydrogen-bonded networks consisting of amino acid residues and water molecules. We suggest that long-distance proton transfer and conformational coupling is facilitated by hydrogen-bonded networks that often involve more than one protein subunit. Proton-storing Asp/Glu dyads, such as the D1-E65/D2-E312 dyad connected to a complex water-wire network, may be particularly important for coupling protonation states to the protein conformation. Clusters of carboxylic amino acids could participate in proton management at the lumenal surface of PSII. We propose that rather than having a classical hydrophobic protein interior, the lumenal side of PSII resembles a complex polyelectrolyte with evolutionary optimized hydrogen-bonding networks. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Structure-based simulation of linear optical spectra of the CP43 core antenna of photosystem II

The linear optical spectra (absorbance, linear dichroism, circular dichroism, fluorescence) of the CP43 (PsbC) antenna of the photosystem II core complex (PSIIcc) pertaining to the S(0) → S(1) (Q(Y)) transitions of the chlorophyll (Chl) a pigments are simulated by applying a combined quantum chemical/electrostatic method to obtain excitonic couplings and local transition energies (site energies) on the basis of the 2.9 Å resolution crystal structure (Guskov et al., Nat Struct Mol Biol 16:334-342, 2009). The electrostatic calculations identify three Chls with low site energies (Chls 35, 37, and 45 in the nomenclature of Loll et al. (Nature 438:1040-1044, 2005). A refined simulation of experimental spectra of isolated CP43 suggests a modified set of site energies within 143 cm(-1) of the directly calculated values (root mean square deviation: 80 cm(-1)). In the refined set, energy sinks are at Chls 37, 43, and 45 in agreement with earlier fitting results (Raszewski and Renger, J Am Chem Soc 130:4431-4446, 2008). The present structure-based simulations reveal that a large part of the redshift of Chl 37 is due to a digalactosyldiacylglycerol lipid. This finding suggests a new role for lipids in PSIIcc, namely the tuning of optical spectra and the creation of an excitation energy funnel towards the reaction center. The analysis of electrostatic pigment-protein interactions is used to identify amino acid residues that are of potential interest for an experimental approach to an assignment of site energies and energy sinks by site-directed mutagenesis.

Calculation of chromophore excited state energy shifts in response to molecular dynamics of pigment-protein complexes

The absorption and energy transfer properties of photosynthetic pigments are strongly influenced by their local environment or "site." Local electrostatic fields vary in time with protein and chromophore molecular movement and thus transiently influence the excited state transition properties of individual chromophores. Site-specific information is experimentally inaccessible in many light-harvesting pigment-proteins due to multiple chromophores with overlapping spectra. Full quantum mechanical calculations of each chromophores excited state properties are too computationally demanding to efficiently calculate the changing excitation energies along a molecular dynamics trajectory in a pigment-protein complex. A simplified calculation of electrostatic interactions with each chromophores ground to excited state transition, the so-called charge density coupling (CDC) for site energy, CDC, has previously been developed to address this problem. We compared CDC to more rigorous quantum chemical calculations to determine its accuracy in computing excited state energy shifts and their fluctuations within a molecular dynamics simulation of the bacteriochlorophyll containing light-harvesting Fenna-Mathews-Olson (FMO) protein. In most cases CDC calculations differed from quantum mechanical (QM) calculations in predicting both excited state energy and its fluctuations. The discrepancies arose from the inability of CDC to account for the differing effects of charge on ground and excited state electron orbitals. Results of our study show that QM calculations are indispensible for site energy computations and the quantification of contributions from different parts of the system to the overall site energy shift. We suggest an extension of QM/MM methodology of site energy shift calculations capable of accounting for long-range electrostatic potential contributions from the whole system, including solvent and ions.

The role of the environment in electronic energy transfer: a molecular modeling perspective

The key role of the environment in electronic energy transfer has been underscored in recent experimental and theoretical studies. In this perspective, we provide an overview of novel quantum-mechanical methodologies aimed at describing environment effects in energy transfers. The techniques described include continuum dielectric and atomistic descriptions of the surroundings. We discuss the advantages and limitations of each technique, as well as the main insights that have emerged from their application to solvated dyads and photosynthetic pigment-protein complexes. We finally highlight the aspects that still need to be solved in order to provide a full theoretical route to the study of energy transfer phenomena in complex environments.

Structural and mechanistic investigations of photosystem II through computational methods

The advent of oxygenic photosynthesis through water oxidation by photosystem II (PSII) transformed the planet, ultimately allowing the evolution of aerobic respiration and an explosion of ecological diversity. The importance of this enzyme to life on Earth has ironically been paralleled by the elusiveness of a detailed understanding of its precise catalytic mechanism. Computational investigations have in recent years provided more and more insights into the structural and mechanistic details that underlie the workings of PSII. This review will present an overview of some of these studies, focusing on those that have aimed at elucidating the mechanism of water oxidation at the CaMn₄ cluster in PSII, and those exploring the features of the structure and dynamics of this enzyme that enable it to catalyse this energetically demanding reaction. This article is part of a Special Issue entitled: Photosystem II.

Time-dependent atomistic view on the electronic relaxation in light-harvesting system II

Aiming at a better understanding of the molecular details in light absorption during photosynthesis, spatial and temporal correlation functions as well as spectral densities have been determined. At the focus of the present study are the light-harvesting II complexes of the purple bacterium Rhodospirillum molischianum. The calculations are based on a time-dependent combination of molecular dynamics simulations and quantum chemistry methods. Using a 12 ps long trajectory, different quantum chemical methods have been compared to each other. Furthermore, several approaches to determine the couplings between the individual chromophores have been tested. Correlations between energy gap fluctuations of different individual pigments are analyzed but found to be negligible. From the energy gap fluctuations, spectral densities are extracted which serve as input for calculations of optical properties and exciton dynamics. To this end, the spectral densities are tested by determining the linear absorption of the complete two-ring system. One important difference from earlier studies is given by the severely extended length of the trajectory along which the quantum chemical calculations have been performed. Due to this extension, more accurate and reliable data have been obtained in the low frequency regime which is important in the dynamics of electronic relaxation.

Zwitterion modulation of O(2)-evolving activity of cyanobacterial photosystem II

Photosystem II (PSII) is the only enzyme in nature that can catalyze the challenging catalytic photooxidation of H(2)O into four protons, four electrons, and O(2). Slowing down turnover of the O(2)-evolving complex (OEC) is a plausible approach to gain mechanistic information on the reaction. However, modulating the kinetics of the reaction without perturbing the active site is a challenge. In this study, it is shown that the steady-state activity of cyanobacterial PSII is inhibited by small zwitterions, such as glycine betaine and β-alanine. We show that the binding of zwitterions is nondenaturing, is highly reversible, and results in the decrease of the rate of catalytic turnover by ∼50% in the presence of excess zwitterion. Control measurements of photoinduced electron transfer in O(2)-inactive PSII show that the inhibition by zwitterions is the result of a specific decrease in the rate of catalytic turnover of the OEC. Recovery of activity upon addition of an exogenous proton carrier (HCO(3)(-)) provides evidence that proton-transfer pathways, thought to be essential for the relay of protons from the OEC to the lumen, are affected. Interestingly, no inhibition is observed for spinach PSII, suggesting that zwitterions act specifically by binding to the extrinsic proteins on the lumenal side of PSII, which differ significantly between plants and cyanobacteria, to slow proton transfer on the electron donor side of PSII.

Tracking the flow of water through photosystem II using molecular dynamics and streamline tracing

The CaMn(4) cluster of the oxygen-evolving complex (OEC) of photosynthesis catalyzes the light-driven splitting of water into molecular oxygen, protons, and electrons. The OEC is buried within photosystem II (PSII), a multisubunit integral membrane protein complex, and water must find its way to the CaMn(4) cluster by moving through protein. Channels for water entrance, and proton and oxygen exit, have previously been proposed following the analysis of cavities found within X-ray structures of PSII. However, these analyses do not account for the dynamic motion of proteins and cannot track the movement of water within PSII. To study water dynamics in PSII, we performed molecular dynamics simulations and developed a novel approach for the visualization of water diffusion within protein based on a streamline tracing algorithm used in fluid dynamics and diffusion tensor imaging. We identified a system of branching pathways of water diffusion in PSII leading to the OEC that connect to a number of distinct entrance points on the lumenal surface. We observed transient changes in the connections between channels and entrance points that served to moderate both the flow of water near the OEC and the exchange of water inside and outside of the protein. Water flow was significantly altered in simulations lacking the OEC which were characterized by a simpler and wider channel with only two openings, consistent with the creation of an ion channel that allows entry of Mn(2+), Ca(2+), and Cl(-) as required for construction of the CaMn(4) cluster.

Uncovering channels in photosystem II by computer modelling: current progress, future prospects, and lessons from analogous systems

Even prior to the publication of the crystal structures for photosystem II (PSII), it had already been suggested that water, O(2) and H(+) channels exist in PSII to achieve directed transport of these molecules, and to avoid undesirable side reactions. Computational efforts to uncover these channels and investigate their properties are still at early stages, and have so far only been based on the static PSII structure. The rationale behind the proposals for such channels and the computer modelling studies thus far are reviewed here. The need to take the dynamic protein into account is then highlighted with reference to the specific issues and techniques applicable to the simulation of each of the three channels. In particular, lessons are drawn from simulation studies on other protein systems containing similar channels.

Toward understanding molecular mechanisms of light harvesting and charge separation in photosystem II

Conversion of light energy in photosynthesis is extremely fast and efficient, and understanding the nature of this complex photophysical process is challenging. This review describes current progress in understanding molecular mechanisms of light harvesting and charge separation in photosystem II (PSII). Breakthroughs in X-ray crystallography have allowed the development and testing of more detailed kinetic models than have previously been possible. However, due to the complexity of the light conversion processes, satisfactory descriptions remain elusive. Recent advances point out the importance of variations in the photochemical properties of PSII in situ in different thylakoid membrane regions as well as the advantages of combining sophisticated time-resolved spectroscopic experiments with atomic level computational modeling which includes the effects of molecular dynamics.

Toward understanding molecular mechanisms of light harvesting and charge separation in photosystem II

Conversion of light energy in photosynthesis is extremely fast and efficient, and understanding the nature of this complex photophysical process is challenging. This review describes current progress in understanding molecular mechanisms of light harvesting and charge separation in photosystem II (PSII). Breakthroughs in X-ray crystallography have allowed the development and testing of more detailed kinetic models than have previously been possible. However, due to the complexity of the light conversion processes, satisfactory descriptions remain elusive. Recent advances point out the importance of variations in the photochemical properties of PSII in situ in different thylakoid membrane regions as well as the advantages of combining sophisticated time-resolved spectroscopic experiments with atomic level computational modeling which includes the effects of molecular dynamics.

Crystal structure of cyanobacterial photosystem II at 3.0 A resolution: a closer look at the antenna system and the small membrane-intrinsic subunits

Photosystem II (PSII) is a homodimeric protein-cofactor complex embedded in the thylakoid membrane that catalyses light-driven charge separation accompanied by the water splitting reaction during oxygenic photosynthesis. In the first part of this review, we describe the current state of the crystal structure at 3.0 A resolution of cyanobacterial PSII from Thermosynechococcus elongatus [B. Loll et al., Towards complete cofactor arrangement in the 3.0 A resolution structure of photosystem II, Nature 438 (2005) 1040-1044] with emphasis on the core antenna subunits CP43 and CP47 and the small membrane-intrinsic subunits. The second part describes first the general theory of optical spectra and excitation energy transfer and how the parameters of the theory can be obtained from the structural data. Next, structure-function relationships are discussed that were identified from stationary and time-resolved experiments and simulations of optical spectra and energy transfer processes.

Spectral signatures of photosynthesis. I. Review of Earth organisms

Why do plants reflect in the green and have a "red edge" in the red, and should extrasolar photosynthesis be the same? We provide (1) a brief review of how photosynthesis works, (2) an overview of the diversity of photosynthetic organisms, their light harvesting systems, and environmental ranges, (3) a synthesis of photosynthetic surface spectral signatures, and (4) evolutionary rationales for photosynthetic surface reflectance spectra with regard to utilization of photon energy and the planetary light environment. We found the "near-infrared (NIR) end" of the red edge to trend from blue-shifted to reddest for (in order) snow algae, temperate algae, lichens, mosses, aquatic plants, and finally terrestrial vascular plants. The red edge is weak or sloping in lichens. Purple bacteria exhibit possibly a sloping edge in the NIR. More studies are needed on pigment-protein complexes, membrane composition, and measurements of bacteria before firm conclusions can be drawn about the role of the NIR reflectance. Pigment absorbance features are strongly correlated with features of atmospheric spectral transmittance: P680 in Photosystem II with the peak surface incident photon flux density at approximately 685 nm, just before an oxygen band at 687.5 nm; the NIR end of the red edge with water absorbance bands and the oxygen A-band at 761 nm; and bacteriochlorophyll reaction center wavelengths with local maxima in atmospheric and water transmittance spectra. Given the surface incident photon flux density spectrum and resonance transfer in light harvesting, we propose some rules with regard to where photosynthetic pigments will peak in absorbance: (1) the wavelength of peak incident photon flux; (2) the longest available wavelength for core antenna or reaction center pigments; and (3) the shortest wavelengths within an atmospheric window for accessory pigments. That plants absorb less green light may not be an inefficient legacy of evolutionary history, but may actually satisfy the above criteria.

Syntheses and energy transfer in multiporphyrinic arrays self-assembled with hydrogen-bonding recognition groups and comparison with covalent steroidal models

A number of new porphyrins equipped with complementary triple hydrogen-bonding groups were synthesized in good yields. Self-assembly was investigated by NMR spectroscopy, dynamic light scattering (DLS), and atomic force microscopy (AFM). These artificial antenna systems were further characterized by stationary and time-resolved fluorescence techniques to investigate several yet unsolved questions on the mechanism of excitation energy transfer (EET) in supramolecular systems. For example, the photophysics of a simple D--U[triple chemical bond]P--A dyad was studied, in which donor D and acceptor A are ZnII- metalated and free-base porphyrins, respectively, and U (uracyl) and P (2,6-diacetamidopyridyl) are complementary hydrogen-bonding groups linked by flexible spacers. In this dyad, the EET occurs with about 20 % efficiency with a lifetime of 14 ps. Reversal of the nonsymmetric triple hydrogen-bonding groups to give a A--U[triple chemical bond]P--D construct results in an EET efficiency of about 25 % and a lifetime of 19 ps. Thus, there is a slight directionality of EET mediated by these asymmetric triple hydrogen-bonding units tethered to flexible spacers. In polymeric systems of the type P-D-P[triple chemical bond]U-A-U[triple chemical bond]P-D-P, or U-D-U[triple chemical bond]P-A-P[triple chemical bond]U-D-U, the EET efficiency doubles as each donor is flanked by two acceptors. Because doubling the probability of photon capture doubles the EET efficiency, there is no energy amplification, which is consistent with the "antenna effect". For these polymeric systems, AFM images and DLS data indicate large rodlike assemblies of a few hundred nanometers, whereas the components form much smaller aggregates under the same conditions. To understand the importance of the flexible hydrogen-bonding zipper, three different covalently bridged D-B-A molecules were synthesized in which the bridge B is a rigid steroidal system and the same ester chemistry was used to link the porphyrins to each end of the steroid. The geometry inferred from molecular modeling of D-B-A indicates geometric similarities between B and some conformations of the --P[triple chemical bond]U-- supramolecular bridge. Although the EET efficiency is a factor of two greater for the steroidal systems relative to the supramolecular dyads, the rate is 50-80 times slower, but still slightly faster than that predicted by Förster-type mechanisms. Circular dichrosim (CD) spectra provide a conformational sampling of the porphyrin groups appended on the steroidal skeleton, thus allowing an estimation of the orientation factor kappa for the transition dipole moments, which significantly affects the EET rate. We conclude that the flexible hydrogen-bonded linked systems are adaptive and have variable geometries with foldamers in which the D and A groups can approach well under 1 nm. In these folded conformations, a rapid EET process occurs, probably also involving a Dexter-type exchange mechanism, thus explaining the fast EET relative to the rigid steroidal compounds. This study predicts that it is indeed possible to build large supramolecular antennas and the component design and supramolecular dynamics are essential features that dictate EET rates and efficiencies.