Exploring the energetics of water permeation in photosystem II by multiple steered molecular dynamics simulations

The Mn(4)Ca 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 Mn(4)Ca cluster by moving through protein. Molecular dynamics simulations were used to determine the energetic barriers for water permeation though PSII extrinsic proteins. Potentials of mean force (PMFs) for water were derived by using the technique of multiple steered molecular dynamics (MSMD). Calculation of free energy profiles for water permeation allowed us to characterize previously identified water channels, and discover new pathways for water movement toward the Mn(4)Ca cluster. Our results identify the main constriction sites in these pathways which may serve as selectivity filters that restrict both the access of solutes detrimental to the water oxidation reaction and loss of Ca(2+) and Cl(-) from the active site.     

Not only enthalpy: large entropy contribution to ion permeation barriers in single-file channels

The effect of channel length on the barrier for potassium ion permeation through single-file channels has been studied by means of all-atom molecular dynamics simulations. Using series of peptidic gramicidin-like and simplified ring-structured channels, both embedded in model membranes, we obtained two distinct types of behavior: saturation of the central free energy barriers for peptidic channels and a linear increase in simplified ring-structured channels with increasing channel length. The saturation of the central free energy barrier for the peptidic channels occurs at relatively short lengths, and it is correlated with the desolvation from the bulk water. Remarkably, decomposition of free energy barriers into enthalpic and entropic terms reveals an entropic cost for ion permeation. Furthermore, this entropic cost dominates the ion permeation free energy barrier, since the corresponding free energy contribution is higher than the enthalpic barrier. We conclude that the length dependence of the free energy is enthalpy-dominated, but the entropy is the major contribution to the permeation barrier. The decrease in rotational water motion and the reduction of channel mobility are putative origins for the overall entropic penalty.     

Structural characteristics of channels and pathways in photosystem II including the identification of an oxygen channel

Here we use crystal structures to investigate and review channels and pathways for the transfer of substrates (water, plastoquinone (PQ)) and products (electrons, protons, oxygen, reduced PQ (PQH(2))) to, and from, the redox active catalytic sites of photosystem II (PSII). A putative oxygen channel has been identified which is about 21A in length, leading from the water splitting site to the lumen. This channel follows a path along the lumenal surface of CP43, passing across the interface of the large extrinsic loop which joins the fifth and sixth transmembrane helices of this chlorophyll binding protein. In so doing it seems to minimise interactions with the excited states of chlorophylls bound within the PSII complex, especially those that constitute the primary electron donor, P680. Two additional channels leading from the water splitting site, and also exiting at the lumen, were also identified. Their hydrophilic nature suggests that they probably facilitate the delivery of water to, and protons from, the catalytic site. Also discussed are unique features in the electron transfer pathway of PSII, as compared with those of purple photosynthetic bacteria, and structural implications of the PSII Q(B)-site in terms of PQ protonation and PQ/PQH(2) diffusion.     

QM/MM computational studies of substrate water binding to the oxygen-evolving centre of photosystem II

This paper reports computational studies of substrate water binding to the oxygen-evolving centre (OEC) of photosystem II (PSII), completely ligated by amino acid residues, water, hydroxide and chloride. The calculations are based on quantum mechanics/molecular mechanics hybrid models of the OEC of PSII, recently developed in conjunction with the X-ray crystal structure of PSII from the cyanobacterium Thermosynechococcus elongatus. The model OEC involves a cuboidal Mn3CaO4Mn metal cluster with three closely associated manganese ions linked to a single mu4-oxo-ligated Mn ion, often called the 'dangling manganese'. Two water molecules bound to calcium and the dangling manganese are postulated to be substrate molecules, responsible for dioxygen formation. It is found that the energy barriers for the Mn(4)-bound water agree nicely with those of model complexes. However, the barriers for Ca-bound waters are substantially larger. Water binding is not simply correlated to the formal oxidation states of the metal centres but rather to their corresponding electrostatic potential atomic charges as modulated by charge-transfer interactions. The calculations of structural rearrangements during water exchange provide support for the experimental finding that the exchange rates with bulk 18 O-labelled water should be smaller for water molecules coordinated to calcium than for water molecules attached to the dangling manganese. The models also predict that the S1-->S2 transition should produce opposite effects on the two water-exchange rates.     

Photosynthetic generation of oxygen

The oxygen in the atmosphere is derived from light-driven oxidation of water at a catalytic centre contained within a multi-subunit enzyme known as photosystem II (PSII). PSII is located in the photosynthetic membranes of plants, algae and cyanobacteria and its oxygen-evolving centre (OEC) consists of four manganese ions and a calcium ion surrounded by a highly conserved protein environment. Recently, the structure of PSII was elucidated by X-ray crystallography thus revealing details of the molecular architecture of the OEC. This structural information, coupled with an extensive knowledge base derived from a wide range of biophysical, biochemical and molecular biological studies, has provided a framework for understanding the chemistry of photosynthetic oxygen generation as well as opening up debate about its evolutionary origin.     

Water and ion permeation in bAQP1 and GlpF channels: a kinetic Monte Carlo study

The kinetic Monte Carlo reaction-path-following technique is applied to determine the lowest-energy water pathway and the coordinating amino acids in bAQP1 and GlpF channels, both treated as rigid. In bAQP1, water molecules pass through the pore between the asparagine-proline-alanine (NPA) and selectivity filter (SF) sites one at a time. The water chain is interrupted at the SF where one water forms three stable hydrogen bonds with protein atoms. In this SF, water's conformation depends on the protonation locus of H182. In GlpF, two water molecules bond simultaneously to the NPA asparagines and pass through the SF in zigzag fashion. No water single-file forms in rigid GlpF. To accommodate a single file of waters requires narrowing the GlpF pore. Our results reveal that in both proteins a proposed bipolar water arrangement is thermally disrupted in the NPA region, especially in the cytoplasmic part of the pore. The equilibrium hydrogen-bonded chain is occasionally interrupted in the hydrophobic zones adjacent to the NPA motifs. The permeation of alkali cations through bAQP1 and GlpF is barred due to a large free-energy barrier in the NPA region as well as a large energy barrier blocking entry from the cytoplasm. Permeation of halides is prevented due to two large energy barriers in the cytoplasmic and periplasmic pores as well as a large free-energy barrier barring entry from the periplasm. Our results, based on modeling charge permeation, support an electrostatic rather than orientational basis for proton exclusion. Binding within the aquaporin pore cannot compensate sufficiently for dehydration of the protonic charge; there is also an electrostatic barrier in the NPA region blocking proton transport. The highly ordered single file of waters, which is drastically interrupted at the SF of bAQP1, may also contribute to proton block.     

Neurophilin-1 is a downstream target of transcription factor Ets-1 in human umbilical vein endothelial cells

Photosystem II complex (PSII) of thylakoid membranes uses light energy to oxidise extremely stable water and produce oxygen (2H(2)O-->O(2)+4H(+)+4e(-)). PSII is compared with cytochrome c oxidase that catalyses the opposite reaction coupled to proton translocation. Cytochrome c oxidase has proton and water channels, and a tentative oxygen channel. I propose that functional PSII complexes also need a specific oxygen channel to direct O(2) from the water molecules bound to specific Mn atoms of the Mn cluster within PSII out to the membrane surface. The function of this channel will be to prevent oxygen being accessible to the radical pair P680(+)Pheo(-), thereby preventing singlet oxygen generation from the triplet P680 state in functional PSII. The important role of singlet oxygen in structurally perturbed non-functional photosystem II is also discussed.     

Perturbations at the chloride site during the photosynthetic oxygen-evolving cycle

Photosystem II (PSII) catalyzes the oxidation of water to O₂ at the manganese-containing, oxygen-evolving complex (OEC). Photoexcitation of PSII results in the oxidation of the OEC; four sequential oxidation reactions are required for the generation and release of molecular oxygen. Therefore, with flash illumination, the OEC cycles among five S _n states. Chloride depletion inhibits O₂ evolution. However, the binding site of chloride in the OEC is not known, and the role of chloride in oxygen evolution has not as yet been elucidated. We have employed reaction-induced FT-IR spectroscopy and selective flash excitation, which cycles PSII samples through the S state transitions. On the time scale employed, these FT-IR difference spectra reflect long-lived structural changes in the OEC. Bromide substitution supports oxygen evolution and was used to identify vibrational bands arising from structural changes at the chloride-binding site. Contributions to the vibrational spectrum from bromide-sensitive bands were observed on each flash. Sulfate treatment led to an elimination of oxygen evolution activity and of the FT-IR spectra assigned to the S₃ to S₀ (third flash) and S₀ to S₁ transitions (fourth flash). However, sulfate treatment changed, but did not eliminate, the FT-IR spectra obtained with the first and second flashes. Solvent isotope exchange in chloride-exchanged samples suggests flash-dependent structural changes, which alter protein dynamics during the S state cycle. Supported by NSF MCB 03-55421.     

Access channels and methanol binding site to the CaMn4 cluster in Photosystem II based on solvent accessibility simulations, with implications for substrate water access

Given the tightly packed environment of Photosystem II (PSII), channels are expected to exist within the protein to allow the movement of small molecules to and from the oxygen evolving centre. In this report, we calculate solvent contact surfaces from the PSII crystal structures to identify such access channels for methanol and water molecules. In a previous study of the effects of methanol on the EPR split S1-, S3-, and S0-signals [Su et al. (2006) Biochemistry 45, 7617-7627], we proposed that methanol binds to one and the same Mn ion in all S-states. We find here that while channels of methanol dimensions were able to make contact with the CaMn4 cluster, only 3Mn and 4Mn were accessible to methanol. Combining this observation with spectroscopic data in the literature, we propose that 3Mn is the ion to which methanol binds. Furthermore, by calculating solvent contact surfaces for water, we found analogous and more extensive water accessible channels within PSII. On the basis of their structure, orientation, and electrostatic properties, we propose functional assignments of these channels as passages for substrate water access to the CaMn4 cluster, and for the exit of O2 and H+ that are released during water oxidation. Finally, we discuss the possible existence of a gating mechanism for the control of substrate water access to the CaMn4 cluster, based on the observation of a gap within the channel system that is formed by Ca2+ and several mechanistically very significant residues in the vicinity of the cluster.     

Fluorescence kinetics of PSII crystals containing Ca2 + or Sr2 + in the oxygen evolving complex

Photosystem II (PSII) is the pigment–protein complex which converts sunlight energy into chemical energy by catalysing the process of light-driven oxidation of water into reducing equivalents in the form of protons and electrons. Three-dimensional structures from x-ray crystallography have been used extensively to model these processes. However, the crystal structures are not necessarily identical to those of the solubilised complexes. Here we compared picosecond fluorescence of solubilised and crystallised PSII core particles isolated from the thermophilic cyanobacterium Thermosynechococcus elongatus. The fluorescence of the crystals is sensitive to the presence of artificial electron acceptors (K₃Fe(CN)₃) and electron transport inhibitors (DCMU). In PSII with reaction centres in the open state, the picosecond fluorescence of PSII crystals and solubilised PSII is indistinguishable. Additionally we compared picosecond fluorescence of native PSII with PSII in which Ca² in the oxygen evolving complex (OEC) is biosynthetically replaced by Sr^2 +. With the Sr^2 + replaced OEC the average fluorescence decay slows down slightly (81 ps to 85 ps), and reaction centres are less readily closed, indicating that both energy transfer/trapping and electron transfer are affected by the replacement.     

Carbonic anhydrase activity of the photosystem II OEC33 protein from pea

The purpose of this study was to identify the location of one of the two sources of carbonic anhydrase (CA) activity associated with the PSII complex in chloroplast membranes. We tested the hypothesis that the extrinsic 33 kDa protein, OEC33, associated with the oxygen-evolving complex (OEC), is one source of CA activity. We found that precursor OEC33 expressed in Escherichia coli exhibits CA activity, but the expressed precursors of OEC24 or OEC17 do not. The CA activity of OEC33 remained after treatment at 90 degrees C for 15 min. Additional biochemical evidence supports the hypothesis. Only those wash treatments that remove the OEC33 from PSII also remove CA activity. Both immunoblot and CA activity show that the CA tracks the OEC33, in parallel, when PSII undergoes washing at different CaCl2 concentrations. The OEC33 protein purified by HiTrap Q anion exchange chromatography has CA activity that is inhibited by an antibody against OEC33. PSII membranes washed with 1 M CaCl2 to remove OEC33 can be reconstituted either with extracted, purified, OEC33 or with the E. coli-expressed precursor OEC33. Reconstitution partially restores both oxygen evolution and CA activity. For maximal CA activity, OEC33 requires manganese as a cofactor.     

Inhibition of photosynthesis by viral infection: Effect on PSII structure and function

The effect of viral infection on photosynthesis was investigated in Nicotiana benthamiana Gray plants infected with different strains of pepper and paprika mild mottle viruses (PMMoV and PaMMoV) and chimeric viral genomes derived from them. In both symptomatic and asymptomatic leaves of virus-infected plants, photosynthetic electron transport in photosystem II (PSII) was reduced. In all cases analyzed, viral infection affected the polypeptide pattern of the oxygen-evolving complex (OEC) in thylakoid membranes. The levels of both the 24 and 16 kDa proteins were reduced to a differing extent when compared with the levels in healthy control. This loss of the OEC extrinsic proteins affected the oxygen evolution rates of thylakoid membranes and leaves from infected plants. Additionally, viral coat protein (CP) was found associated with the chloroplasts and the thylakoid membranes of the infected plants. The CP accumulation level was dependent upon both the post-infection time and the virus analyzed, but independent of the CP itself since hybrid viruses did not behave as their parental viruses with the same CP, with respect to PSII inhibition, CP accumulation rates and OEC protein levels. Modulated chlorophyll (Chl) fluorescence and oxygen evolution measurements carried out in both types of leaves showed that the quantum yield of PSII electron transport was diminished in infected plants with respect to those of control plants. The decrease in electron transport efficiency was mainly caused by a reduction in the fraction of open reaction centers. The infected plants also showed a reduction in the efficiency of excitation capture in PSII by photoprotective thermal dissipation of excess excitation energy.     

Ion permeation in normal and batrachotoxin-modified Na+ channels in the squid giant axon

Na+ permeation through normal and batrachotoxin (BTX)-modified squid axon Na+ channels was characterized. Unmodified and toxin-modified Na+ channels were studied simultaneously in outside-out membrane patches using the cut-open axon technique. Current-voltage relations for both normal and BTX-modified channels were measured over a wide range of Na+ concentrations and voltages. Channel conductance as a function of Na+ concentration curves showed that within the range 0.015-1 M Na+ the normal channel conductance is 1.7-2.5-fold larger than the BTX-modified conductance. These relations cannot be fitted by a simple Langmuir isotherm. Channel conductance at low concentrations was larger than expected from a Michaelis-Menten behavior. The deviations from the simple case were accounted for by fixed negative charges located in the vicinity of the channel entrances. Fixed negative charges near the pore mouths would have the effect of increasing the local Na+ concentration. The results are discussed in terms of energy profiles with three barriers and two sites, taking into consideration the effect of the fixed negative charges. Either single- or multi-ion pore models can account for all the permeation data obtained in both symmetric and asymmetric conditions. In a temperature range of 5-15 degrees C, the estimated Q10 for the conductance of the BTX-modified Na+ channel was 1.53. BTX appears not to change the Na+ channel ion selectively (for the conditions used) or the surface charge located near the channel entrances.     

Quantifying the ion selectivity of the Ca2+ site in photosystem II: evidence for direct involvement of Ca2+ in O2 formation.

Calcium is an essential cofactor in the oxygen-evolving complex (OEC) of photosystem II (PSII). The removal of Ca2+ or its substitution by any metal ion except Sr2+ inhibits oxygen evolution. We used steady-state enzyme kinetics to measure the rate of O2 evolution in PSII samples treated with an extensive series of mono-, di-, and trivalent metal ions in order to determine the basis for the affinity of metal ions for the Ca2+-binding site. Our results show that the Ca2+-binding site in PSII behaves very similarly to the Ca2+-binding sites in other proteins, and we discuss the implications this has for the structure of the site in PSII. Activity measurements as a function of time show that the binding site achieves equilibrium in 4 h for all of the PSII samples investigated. The binding affinities of the metal ions are modulated by the 17 and 23 kDa extrinsic polypeptides; their removal decreases the free energy of binding of the metal ions by 2.5 kcal/mol, but does not significantly change the time required to reach equilibrium. Monovalent ions are effectively excluded from the Ca2+-binding site, exhibiting no inhibition of O2 evolution. Di- and trivalent metal ions with ionic radii similar to that of Ca2+ (0.99 A) bind competitively with Ca2+ and have the highest binding affinity, while smaller metal ions bind more weakly and much larger ones do not bind competitively. This is consistent with a size-selective Ca2+-binding site that has a rigid array of coordinating ligands. Despite the large number of metal ions that competitively replace Ca2+ in the OEC, only Sr2+ is capable of partially restoring activity. Comparing the physical characteristics of the metal ions studied, we identify the pK(a) of the aqua ion as the factor that determines the functional competence of the metal ion. This suggests that Ca2+ is directly involved in the chemistry of water oxidation and is not only a structural cofactor in the OEC. We propose that the role of Ca2+ is to act as a Lewis acid, binding a substrate water molecule and tuning its reactivity.     

Evidence against bicarbonate bound in the O2-evolving complex of photosystem II

The oxidation of water to molecular oxygen by photosystem II (PSII) is inhibited in bicarbonate-depleted media. One contribution to the inhibition is the binding of bicarbonate to the non-heme iron, which is required for efficient electron transfer on the electron-acceptor side of PSII. There are also proposals that bicarbonate is required for formation of O 2 by the manganese-containing O 2-evolving complex (OEC). Previous work indicates that a bicarbonate ion does not bind reversibly close to the OEC, but it remains possible that bicarbonate is bound sufficiently tightly to the OEC that it cannot readily exchange with bicarbonate in solution. In this study, we have used NH 2OH to destroy the OEC, which would release any tightly bound bicarbonate ions from the active site, and mass spectrometry to detect any released bicarbonate as CO 2. The amount of CO 2 per PSII released by the NH 2OH treatment is observed to be comparable to the background level, although N 2O, a product of the reaction of NH 2OH with the OEC, is detected in good yield. These results strongly argue against tightly bound bicarbonate ions in the OEC.     

Azide as a probe of proton transfer reactions in photosynthetic oxygen evolution

In oxygenic photosynthesis, photosystem II (PSII) is the multisubunit membrane protein responsible for the oxidation of water to O₂ and the reduction of plastoquinone to plastoquinol. One electron charge separation in the PSII reaction center is coupled to sequential oxidation reactions at the oxygen-evolving complex (OEC), which is composed of four manganese ions and one calcium ion. The sequentially oxidized forms of the OEC are referred to as the S_n states. S₁ is the dark-adapted state of the OEC. Flash-induced oxygen production oscillates with period four and occurs during the S₃ to S₀ transition. Chloride plays an important, but poorly understood role in photosynthetic water oxidation. Chloride removal is known to block manganese oxidation during the S₂ to S₃ transition. In this work, we have used azide as a probe of proton transfer reactions in PSII. PSII was sulfate-treated to deplete chloride and then treated with azide. Steady state oxygen evolution measurements demonstrate that azide inhibits oxygen evolution in a chloride-dependent manner and that azide is a mixed or noncompetitive inhibitor. This result is consistent with two azide binding sites, one at which azide competes with chloride and one at which azide and chloride do not compete. At pH 7.5, the K_i for the competing site was estimated as 1 mM, and the K_i′ for the uncompetitive site was estimated as 8 mM. Vibrational spectroscopy was then used to monitor perturbations in the frequency and amplitude of the azide antisymmetric stretching band. These changes were induced by laser-induced charge separation in the PSII reaction center. The results suggest that azide is involved in proton transfer reactions, which occur before manganese oxidation, on the donor side of chloride-depleted PSII.     

Damage to the oxygen-evolving complex by superoxide anion, hydrogen peroxide, and hydroxyl radical in photoinhibition of photosystem II

Under strong illumination of a photosystem II (PSII) membrane, endogenous superoxide anion, hydrogen peroxide, and hydroxyl radical were successively produced. These compounds then cooperatively resulted in a release of manganese from the oxygen-evolving complex (OEC) and an inhibition of oxygen evolution activity. The OEC inactivation was initiated by an acceptor-side generated superoxide anion, and hydrogen peroxide was most probably responsible for the transportation of reactive oxygen species (ROS) across the PSII membrane from the acceptor-side to the donor-side. Besides ROS being generated in the acceptor-side induced manganese loss; there may also be a ROS-independent manganese loss in the OEC of PSII. Both superoxide anion and hydroxyl radical located inside the PSII membrane were directly identified by a spin trapping-electron spin resonance (ESR) method in combination with a lipophilic spin trap, 5-(diethoxyphosphoryl)-5-phenethyl-1-pyrroline N-oxide (DEPPEPO). The endogenous hydrogen peroxide production was examined by oxidation of thiobenzamide.     

The evolution of Photosystem II: insights into the past and future

This article attempts to address the molecular origin of Photosystem II (PSII), the central component in oxygenic photosynthesis. It discusses the possible evolution of the relevant cofactors needed for splitting water into molecular O₂ with respect to the following functional domains in PSII: the reaction center (RC), the oxygen evolving complex (OEC), and the manganese stabilizing protein (MSP). Possible ancestral sources of the relevant cofactors are considered, as are scenarios of how these components may have been brought together to produce the intermediate steps in the evolution of PSII. Most importantly, the driving forces that maintained these intermediates for continued adaptation are considered. We then apply our understanding of the evolution of PSII to the bioengineering of a water oxidizing catalyst for utilization of solar energy.     

Charge delocalization in proton channels, I: the aquaporin channels and proton blockage

The explicit contribution to the free energy barrier and proton conductance from the delocalized nature of the excess proton is examined in aquaporin channels using an accurate all-atom molecular dynamics computer simulation model. In particular, the channel permeation free energy profiles are calculated and compared for both a delocalized (fully Grotthuss shuttling) proton and a classical (nonshuttling) hydronium ion along two aquaporin channels, Aqp1 and GlpF. To elucidate the effects of the bipolar field thought to arise from two alpha-helical macrodipoles on proton blockage, free energy profiles were also calculated for computational mutants of the two channels where the bipolar field was eliminated by artificially discharging the backbone atoms. Comparison of the free energy profiles between the proton and hydronium cases indicates that the magnitude of the free energy barrier and position of the barrier peak for the fully delocalized and shuttling proton are somewhat different from the case of the (localized) classical hydronium. The proton conductance through the two aquaporin channels is also estimated using Poisson-Nernst-Planck theory for both the Grotthuss shuttling excess proton and the classical hydronium cation.     

Binding and functions of the two chloride ions in the oxygen-evolving center of photosystem II

Photosynthesis Research - Light-driven water oxidation in photosynthesis occurs at the oxygen-evolving center (OEC) of photosystem II (PSII). Chloride ions (Cl?) are essential for oxygen...