Investigating the mechanisms of photosynthetic proteins using continuum electrostatics.

Computational methods based on continuum electrostatics are widely used in theoretical biochemistry to analyze the function of proteins. Continuum electrostatic methods in combination with quantum chemical and molecular mechanical methods can help to analyze even very complex biochemical systems. In this article, applications of these methods to proteins involved in photosynthesis are reviewed. After giving a short introduction to the basic concepts of the continuum electrostatic model based on the Poisson-Boltzmann equation, we describe the application of this approach to the docking of electron transfer proteins, to the comparison of isofunctional proteins, to the tuning of absorption spectra, to the analysis of the coupling of electron and proton transfer, to the analysis of the effect of membrane potentials on the energetics of membrane proteins, and to the kinetics of charge transfer reactions. Simulations as those reviewed in this article help to analyze molecular mechanisms on the basis of the structure of the protein, guide new experiments, and provide a better and deeper understanding of protein functions.

pH modulates the quinone position in the photosynthetic reaction center from Rhodobacter sphaeroides in the neutral and charge separated states.

The structure of the photosynthetic reaction-center from Rhodobacter sphaeroides has been determined at four different pH values (6.5, 8.0, 9.0, 10.0) in the neutral and in charge separated states. At pH 8.0, in the neutral state, we obtain a resolution of 1.87 A, which is the best ever reported for the bacterial reaction center protein. Our crystallographic data confirm the existence of two different binding positions of the secondary quinone (QB). We observe a new orientation of QB in its distal position, which shows no ring-flip compared to the orientation in the proximal position. Datasets collected for the different pH values show a pH-dependence of the population of the proximal position. The new orientation of QB in the distal position and the pH-dependence could be confirmed by continuum electrostatics calculations. Our calculations are in agreement with the experimentally observed proton uptake upon charge separation. The high resolution of our crystallographic data allows us to identify new water molecules and external residues being involved in two previously described hydrogen bond proton channels. These extended proton-transfer pathways, ending at either of the two oxo-groups of QB in its proximal position, provide additional evidence that ring-flipping is not required for complete protonation of QB upon reduction.

Redox-linked protonation state changes in cytochrome bc1 identified by Poisson-Boltzmann electrostatics calculations.

Cytochrome bc(1) is a major component of biological energy conversion that exploits an energetically favourable redox reaction to generate a transmembrane proton gradient. Since the mechanistic details of the coupling of redox and protonation reactions in the active sites are largely unresolved, we have identified residues that undergo redox-linked protonation state changes. Structure-based Poisson-Boltzmann/Monte Carlo titration calculations have been performed for completely reduced and completely oxidised cytochrome bc(1). Different crystallographically observed conformations of Glu272 and surrounding residues of the cytochrome b subunit in cytochrome bc(1) from Saccharomyces cerevisiae have been considered in the calculations. Coenzyme Q (CoQ) has been modelled into the CoQ oxidation site (Q(o)-site). Our results indicate that both conformational and protonation state changes of Glu272 of cytochrome b may contribute to the postulated gating of CoQ oxidation. The Rieske iron-sulphur cluster could be shown to undergo redox-linked protonation state changes of its histidine ligands in the structural context of the CoQ-bound Q(o)-site. The proton acceptor role of the CoQ ligands in the CoQ reduction site (Q(i)-site) is supported by our results. A modified path for proton uptake towards the Q(i)-site features a cluster of conserved lysine residues in the cytochrome b (Lys228) and cytochrome c(1) subunits (Lys288, Lys289, Lys296). The cardiolipin molecule bound close to the Q(i)-site stabilises protons in this cluster of lysine residues.

Theoretical investigation of the behavior of titratable groups in proteins.

This paper presents a theoretical analysis of the titration behavior of strongly interacting titratable residues in proteins. Strongly interacting titratable residues exist in many proteins such as for instance bacteriorhodopsin, cytochrome c oxidase, cytochrome bc(1), or the photosynthetic reaction center. Strong interaction between titratable groups can lead to irregular titration behavior. We analyze under which circumstances titration curves can become irregular. We demonstrate that conformational flexibility alone can not lead to irregular titration behavior. Strong interaction between titratable groups is a necessary, but not sufficient condition for irregular titration curves. In addition, the two interacting groups also need to titrate in the same pH-range. These two conditions together lead to irregular titration curves. The mutation of a single residue within a cluster of interacting titratable residues can influence the titration behavior of the other titratable residues in the cluster. We demonstrate this effect on a cluster of four interacting residues. This example underlines that mutational studies directed at identifying the role of a certain titratable residue in a cluster of interacting residues should always be accompanied by an analysis of the effect of the mutation on the titration behavior of the other residues.

Negatively charged residues and hydrogen bonds tune the ligand histidine pKa values of Rieske iron-sulfur proteins.

Rieske proteins carry a redox-active iron-sulfur cluster, which is bound by two histidine and two cysteine side chains. The reduction potential of Rieske proteins depends on pH. This pH dependence can be described by two pK(a) values, which have been assigned to the two iron-coordinating histidines. Rieske proteins are commonly grouped into two major classes: Rieske proteins from quinol-oxidizing cytochrome bc complexes, in which the ligand histidines titrate in the physiological pH range, and bacterial ferredoxin Rieske proteins, in which the ligand histidines are protonated at physiological pH. In the study presented here, we have calculated pK(a) values of the cluster ligand histidines using a combined density functional theory/continuum electrostatics approach. Experimental pK(a) values for a bc-type and a ferredoxin Rieske protein could be reproduced. We could identify functionally important differences between the two proteins: hydrogen bonds toward the cluster, which are present in bc-type Rieske proteins, and negatively charged residues, which are present in ferredoxin Rieske proteins. We removed these differences by mutating the proteins in our calculations. The Rieske centers in the mutated proteins have very similar pK(a) values. We thus conclude that the studied structural differences are the main reason for the different pH-titration behavior of the proteins. Interestingly, the shift caused by neutralizing the negative charges in ferredoxin Rieske proteins is larger than the shift caused by removing the hydrogen bonds toward the cluster in bc-type Rieske proteins.