Serguei Vassiliev

Structural-functional role of chloride in photosystem II.

Chloride binding in photosystem II (PSII) is essential for photosynthetic water oxidation. However, the functional roles of chloride and possible binding sites, during oxygen evolution, remain controversial. This paper examines the functions of chloride based on its binding site revealed in the X-ray crystal structure of PSII at 1.9 Å resolution. We find that chloride depletion induces formation of a salt bridge between D2-K317 and D1-D61 that could suppress the transfer of protons to the lumen.

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.

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.

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 Ems 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.

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.

Factors Controlling the Redox Potential of ZnCe6 in an Engineered Bacterioferritin Photochemical 'Reaction Centre'

Photosystem II (PSII) of photosynthesis has the unique ability to photochemically oxidize water. Recently an engineered bacterioferritin photochemical ‘reaction centre’ (BFR-RC) using a zinc chlorin pigment (ZnCe6) in place of its native heme has been shown to photo-oxidize bound manganese ions through a tyrosine residue, thus mimicking two of the key reactions on the electron donor side of PSII. To understand the mechanism of tyrosine oxidation in BFR-RCs, and explore the possibility of water oxidation in such a system we have built an atomic-level model of the BFR-RC using ONIOM methodology. We studied the influence of axial ligands and carboxyl groups on the oxidation potential of ZnCe6 using DFT theory, and finally calculated the shift of the redox potential of ZnCe6 in the BFR-RC protein using the multi-conformational molecular mechanics–Poisson-Boltzmann approach. According to our calculations, the redox potential for the first oxidation of ZnCe6 in the BRF-RC protein is only 0.57 V, too low to oxidize tyrosine. We suggest that the observed tyrosine oxidation in BRF-RC could be driven by the ZnCe6 di-cation. In order to increase the efficiency of tyrosine oxidation, and ultimately oxidize water, the first potential of ZnCe6 would have to attain a value in excess of 0.8 V. We discuss the possibilities for modifying the BFR-RC to achieve this goal.

Water and Oxygen Diffusion Pathways Within Photosystem II. Computational Studies of Controlled Substrate Access and Product Release

The problem of finding and characterizing channels responsible for controlling the movement of small molecules within protein complexes is addressed, using photosystem II (PSII) of photosynthesis as an example. The shortcomings of traditional methods of searching for channels in static crystallographic structures are discussed. Basic concepts and advantages of molecular dynamics (MD) simulations and existing computational approaches based on MD structures to finding water channels are introduced. The use of MD to observe water diffusion within proteins, visualize water flow, and localize water channels is described. This is followed by more detailed characterization and analysis of the found channels, with respect to their permitivity for water. The second part of the chapter focuses on oxygen channels within PSII. The experimental and computational techniques for characterization of oxygen channels in proteins are reviewed. A powerful computational technique, based on analysis of 3D energy maps, is introduced. The construction of oxygen 3D free energy maps for PSII and their analysis using a wavefront propagation technique is described to find oxygen diffusion pathways and characterize the rate of oxygen diffusion within PSII.