Mechanistic insights into ferredoxin-NADP(H) reductase catalysis involving the conserved glutamate in the active site.

Plant-type ferredoxin-NADP(H) reductases (FNRs) are flavoenzymes harboring one molecule of noncovalently bound flavin adenine dinucleotide that catalyze reversible reactions between obligatory one-electron carriers and obligatory two-electron carriers. A glutamate next to the C-terminus is strictly conserved in FNR and has been proposed to function as proton donor/acceptor during catalysis. However, experimental studies of this proposed function led to contradicting conclusions about the role of this glutamate in the catalytic mechanism. In the present work, we study the titration behavior of the glutamate in the active site of FNR using theoretical methods. Protonation probabilities for maize FNR were computed for the reaction intermediates of the catalytic cycle by Poisson-Boltzmann electrostatic calculations and Metropolis Monte Carlo titration. The titration behavior of the highly conserved glutamate was found to vary depending on the bound substrates NADP(H) and ferredoxin and also on the redox states of these substrates and the flavin adenine dinucleotide. Our results support the involvement of the glutamate in the FNR catalytic mechanism not only as a proton donor but also as a key residue for stabilizing and destabilizing reaction intermediates. On the basis of our findings, we propose a model rationalizing the function of the glutamate in the reaction cycle, which allows reinterpretation of previous experimental results.

Simulating the proton transfer in gramicidin A by a sequential dynamical Monte Carlo method.

The large interest in long-range proton transfer in biomolecules is triggered by its importance for many biochemical processes such as biological energy transduction and drug detoxification. Since long-range proton transfer occurs on a microsecond time scale, simulating this process on a molecular level is still a challenging task and not possible with standard simulation methods. In general, the dynamics of a reactive system can be described by a master equation. A natural way to describe long-range charge transfer in biomolecules is to decompose the process into elementary steps which are transitions between microstates. Each microstate has a defined protonation pattern. Although such a master equation can in principle be solved analytically, it is often too demanding to solve this equation because of the large number of microstates. In this paper, we describe a new method which solves the master equation by a sequential dynamical Monte Carlo algorithm. Starting from one microstate, the evolution of the system is simulated as a stochastic process. The energetic parameters required for these simulations are determined by continuum electrostatic calculations. We apply this method to simulate the proton transfer through gramicidin A, a transmembrane proton channel, in dependence on the applied membrane potential and the pH value of the solution. As elementary steps in our reaction, we consider proton uptake and release, proton transfer along a hydrogen bond, and rotations of water molecules that constitute a proton wire through the channel. A simulation of 8 mus length took about 5 min on an Intel Pentium 4 CPU with 3.2 GHz. We obtained good agreement with experimental data for the proton flux through gramicidin A over a wide range of pH values and membrane potentials. We find that proton desolvation as well as water rotations are equally important for the proton transfer through gramicidin A at physiological membrane potentials. Our method allows to simulate long-range charge transfer in biological systems at time scales, which are not accessible by other methods.

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.

Effects of histidine protonation and phosphorylation on histidine-containing phosphocarrier protein structure, dynamics, and physicochemical properties.

Previous structural studies of the histidine-containing phosphocarrier protein (HPr) have shown that active site residue His15 can adopt two distinct conformations which were termed OPEN and CLOSED. Using molecular dynamics simulations and protonation probability calculations, we were able to show that these two conformations correspond to different protonation forms of the histidine ring. The CLOSED-to-OPEN transition requires His15 to adopt a conformation with higher energy, which is compensated by the favorable energetic consequences of protonation. Calculations of the conformational energy of His15 show that HPr exists mainly in the CLOSED form at pH 7. The very low apparent pKa value (3.2-4.5) of the CLOSED conformation and the fact that the imidazole ring of residue 15 is primarily unprotonated at Ndelta1 at neutral pH ensure that His15 is ideally primed to be specifically phosphorylated at Ndelta1. In contrast to unphosphorylated HPr, the phosphorylated form exhibits no conformational transitions, and the CLOSED state is stable even for the protonated imidazole ring due to favorable interactions between the phosphate group and the backbone of Ala16 and Arg17. These observations from MD simulations are confirmed by a simple four-microstate model which can explain both the pH-dependent conformational change of unphosphorylated HPr and the conformational rigidity of phosphorylated HPr. Our study suggests that the predominant CLOSED conformation is relevant for HPr function in the phosphotransfer reaction, while the OPEN form of unphosphorylated HPr might be important for its additional regulatory function, in which an OPEN conformation of His15 is recognized by the transcriptional regulator CcpA.

Effect of HPr phosphorylation on structure, dynamics, and interactions in the course of transcriptional control.

The serine46-phosphorylated form of the bacterial protein HPr fulfils an essential function in carbon catabolite repression (CCR). Using molecular dynamics (MD) we studied the effect of Ser46 phosphorylation on the molecular properties of HPr and its capability to act as the co-repressor of carbon catabolite protein A (CcpA). The calculated pK (a) values for a representative set of HPr(Ser46P) structures indicate that the phosphate group of HPr(Ser46P) exists predominantly in the unprotonated form under neutral conditions. A hydrogen bond detected in HPr(Ser46P) between one phosphate-group oxygen and a side-chain hydrogen of Asn43-an amino acid conserved in all HPr proteins of Gram-positive bacteria that regulate their carbon consumption by CCR-might fulfil an important role in CcpA-HPr(Ser46P) complex formation. MD simulations show that the Ser46P-Asn43 hydrogen bond present in the unbound structure is replaced by intermolecular interactions upon complex formation. The degree to which amino acids in the CcpA-HPr(Ser46P) interface contribute to cofactor binding was analyzed by in silico alanine scanning. Lys307, Arg303, Asp296, Val300, and Tyr295 of CcpA were identified as important amino acids for the CcpA-HPr(Ser46P) interaction. Three of these residues are directly involved in sensing the correct phosphorylation state at His15(HPr) and Ser46(HPr). A substitution of interface residues Val319, Val314, Ser316, Leu321 and Gln320 by alanine showed that these amino acids, which contact helix alpha2 of HPr(Ser46P), play a less prominent role for complex formation.

Formation and characterization of an all-ferrous Rieske cluster and stabilization of the [2Fe-2S]0 core by protonation.

The all-ferrous Rieske cluster, [2Fe-2S](0), has been produced in solution and characterized by protein-film voltammetry and UV-visible, EPR, and Mössbauer spectroscopies. The [2Fe-2S](0) cluster, in the overexpressed soluble domain of the Rieske protein from the bovine cytochrome bc(1) complex, is formed at -0.73 V at pH 7. Therefore, at pH 7, the [2Fe-2S](1+/0) couple is 1.0 V below the [2Fe-2S](2+/1+) couple. The two cluster-bound ferrous irons are both high spin (S = 2), and they are coupled antiferromagnetically (-J > or = 30 cm(-1), H =-2JS1.S2) to give a diamagnetic (S = 0) ground state. The ability of the Rieske cluster to exist in three oxidation states (2+, 1+, and 0) without an accompanying coupled reaction, such as a conformational change or protonation, is highly unusual. However, uncoupled reduction to the [2Fe-2S](0) state occurs at pH > 9.8 only, and at high pH the intact cluster persists in solution for <1 min. At pH < 9.8, the all-ferrous cluster is stabilized significantly by protonation. A combination of experimental data and calculations based on density functional theory suggests strongly that the proton binds to one of the cluster mu(2)-sulfides, consistent with observations that reduced [3Fe-4S] clusters are protonated also. The implications for our understanding of coupled reactions at iron-sulfur clusters and of the factors that determine the relative stabilities of their different oxidation states are discussed.