Search

Recent neutron diffraction studies of photoactive yellow protein (PYP) proposed that the H bond between protonated Glu46 and the chromophore [ionized p-coumaric acid (pCA)] was a low-barrier H bond (LBHB). Using the atomic coordinates of the high-resolution crystal structure, we analyzed the energetics of the short H bond by two independent methods: electrostatic pK(a) calculations and a quantum mechanical/molecular mechanical (QM/MM) approach. (i) In the QM/MM optimized geometry, we reproduced the two short H-bond distances of the crystal structure: Tyr42-pCA (2.50 Å) and Glu46-pCA (2.57 Å). However, the H atoms obviously belonged to the Tyr or Glu moieties, and were not near the midpoint of the donor and acceptor atoms. (ii) The potential-energy curves of the two H bonds resembled those of standard asymmetric double-well potentials, which differ from those of LBHB. (iii) The calculated pK(a) values for Glu46 and pCA were 8.6 and 5.4, respectively. The pK(a) difference was unlikely to satisfy the prerequisite for LBHB. (iv) The LBHB in PYP was originally proposed to stabilize the ionized pCA because deprotonated Arg52 cannot stabilize it. However, the calculated pK(a) of Arg52 and QM/MM optimized geometry suggested that Arg52 was protonated on the protein surface. The short H bond between Glu46 and ionized pCA in the PYP ground state could be simply explained by electrostatic stabilization without invoking LBHB.

The crystal structure of photosystem II (PSII) analyzed at a resolution of 1.9 Å revealed a remarkably short H-bond between redox-active tyrosine Y(Z) and D1-His190 (2.46 Å donor-acceptor distance). Using large-scale quantum mechanical/molecular mechanical (QM/MM) calculations with the explicit PSII protein environment, we were able to reproduce this remarkably short H-bond in the original geometry of the crystal structure in the neutral [Y(Z)O···H···N(ε)-His-N(δ)H···O═Asn] state, but not in the oxidized states, indicating that the neutral state was the one observed in the crystal structure. In addition to the appropriate redox/protonation state of Y(Z) and D1-His190, we found that the presence of a cluster of water molecules played a key role in shortening the distance between Y(Z) and D1-His190. The orientations of the water molecules in the cluster were energetically stabilized by the highly polarized PSII protein environment, where the Ca ion of the oxygen-evolving complex (OEC) and the OEC ligand D1-Glu189 were also involved.

In photosynthesis, light is captured by antenna proteins. These proteins transfer the excitation energy with almost 100% quantum efficiency to the reaction centers, where charge separation takes place. The time scale and pathways of this transfer are controlled by the protein scaffold, which holds the pigments at optimal geometry and tunes their excitation energies (site energies). The detailed understanding of the tuning of site energies by the protein has been an unsolved problem since the first high-resolution crystal structure of a light-harvesting antenna appeared >30 years ago [Fenna RE, Matthews BW (1975) Nature 258:573-577]. Here, we present a combined quantum chemical/electrostatic approach to compute site energies that considers the whole protein in atomic detail and provides the missing link between crystallography and spectroscopy. The calculation of site energies of the Fenna-Matthews-Olson protein results in optical spectra that are in quantitative agreement with experiment and reveals an unexpectedly strong influence of the backbone of two alpha-helices. The electric field from the latter defines the direction of excitation energy flow in the Fenna-Matthews-Olson protein, whereas the effects of amino acid side chains, hitherto thought to be crucial, largely compensate each other. This result challenges the current view of how energy flow is regulated in pigment-protein complexes and demonstrates that attention has to be paid to the backbone architecture.

Cd(2+) binding at the bacterial photosynthetic reaction center (bRC) from Rhodobacter sphaeroides is known to inhibit proton transfer (PT) from bulk solvent to the secondary quinone Q(B). To elucidate this mechanism, we calculated the pK(a) for residues along the water channels connecting Q(B) with the stromal side based on the crystal structures of WT-bRC and Cd(2+)-bound bRC. Upon Cd(2+) binding, we observed the release of two protons from His-H126/128 at the Cd(2+) binding site and significant pK(a) shifts for residues along the PT pathways. Remarkably, Asp-L213 near Q(B), which is proposed to play a significant role in PT, resulted in a decrease in pK(a) upon Cd(2+) binding. The direct electrostatic influence of the Cd(2+)-positive charge on these pK(a) shifts was small. Instead, conformational changes of amino acid side chains induced electrostatically by Cd(2+) binding were the main mechanism for these pK(a) shifts. The long-range electrostatic influence over approximately 12 A between Cd(2+) and Asp-L213 is likely to originate from a set of Cd(2+)-induced successive reorientations of side chains (Asp-H124, His-H126, His-H128, Asp-H170, Glu-H173, Asp-M17, and Asp-L210), which propagate along the PT pathways as a "domino" effect.