Negatively Charged Residues and Hydrogen Bonds Tune the Ligand Histidine pKa Values of Rieske Iron−Sulfur Proteins

Rieske Oxygenases and Other Ferredoxin-Dependent Enzymes: Electron Transfer Principles and Catalytic Capabilities

Enzymes that depend on sophisticated electron transfer via ferredoxins (Fds) exhibit outstanding catalytic capabilities, but despite decades of research, many of them are still not well understood or exploited for synthetic applications. This review aims to provide a general overview of the most important Fd-dependent enzymes and the electron transfer processes involved. While several examples are discussed, we focus in particular on the family of Rieske non-heme iron-dependent oxygenases (ROs). In addition to illustrating their electron transfer principles and catalytic potential, the current state of knowledge on structure-function relationships and the mode of interaction between the redox partner proteins is reviewed. Moreover, we highlight several key catalyzed transformations, but also take a deeper dive into their engineerability for biocatalytic applications. The overall findings from these case studies highlight the catalytic capabilities of these biocatalysts and could stimulate future interest in developing additional Fd-dependent enzyme classes for synthetic applications.

Tools for analyzing protonation states and for tracing proton transfer pathways with examples from the Rb. sphaeroides photosynthetic reaction centers

Protons participate in many reactions. In proteins, protons need paths to move in and out of buried active sites. The vectorial movement of protons coupled to electron transfer reactions establishes the transmembrane electrochemical gradient used for many reactions, including ATP synthesis. Protons move through hydrogen bonded chains of waters and hydroxy side chains via the Grotthuss mechanism and by proton binding and release from acidic and basic residues. MCCE analysis shows that proteins exist in a large number of protonation states. Knowledge of the equilibrium ensemble can provide a rational basis for setting protonation states in simulations that fix them, such as molecular dynamics (MD). The proton path into the Q_B site in the bacterial reaction centers (RCs) of Rb. sphaeroides is analyzed by MD to provide an example of the benefits of using protonation states found by the MCCE program. A tangled web of side chains and waters link the cytoplasm to Q_B. MCCE analysis of snapshots from multiple trajectories shows that changing the input protonation state of a residue in MD biases the trajectory shifting the proton affinity of that residue. However, the proton affinity of some residues is more sensitive to the input structure. The proton transfer networks derived from different trajectories are quite robust. There are some changes in connectivity that are largely restricted to the specific residues whose protonation state is changed. Trajectories with Q_B^•- are compared with earlier results obtained with Q_B [Wei et. al Photosynthesis Research volume 152, pages153-165 (2022)] showing only modest changes. While introducing new methods the study highlights the difficulty of establishing the connections between protein conformation.

Histidine Ligated Iron-Sulfur Peptides

Iron‐sulfur clusters are thought to be ancient cofactors that could have played a role in early protometabolic systems. Thus far, redox active, prebiotically plausible iron‐sulfur clusters have always contained cysteine ligands to the cluster. However, extant iron‐sulfur proteins can be found to exploit other modes of binding, including ligation by histidine residues, as seen with [2Fe‐2S] Rieske and MitoNEET proteins. Here, we investigated the ability of cysteine‐ and histidine‐containing peptides to coordinate a mononuclear Fe²⁺ center and a [2Fe‐2S] cluster and compare their properties with purified iron‐sulfur proteins. The iron‐sulfur peptides were characterized by UV‐vis, circular dichroism, and paramagnetic NMR spectroscopies and cyclic voltammetry. Small (≤6 amino acids) peptides can coordinate [2Fe‐2S] clusters through a combination of cysteine and histidine residues with similar reduction potentials as their corresponding proteins. Such complexes may have been important for early cell‐like systems.

A site-differentiated [4Fe-4S] cluster controls electron transfer reactivity of Clostridium acetobutylicum [FeFe]-hydrogenase I

One of the many functions of reduction–oxidation (redox) cofactors is to mediate electron transfer in biological enzymes catalyzing redox-based chemical transformation reactions. There are numerous examples of enzymes that utilize redox cofactors to form electron transfer relays to connect catalytic sites to external electron donors and acceptors. The compositions of relays are diverse and tune transfer thermodynamics and kinetics towards the chemical reactivity of the enzyme. Diversity in relay design is exemplified among different members of hydrogenases, enzymes which catalyze reversible H₂ activation, which also couple to diverse types of donor and acceptor molecules. The [FeFe]-hydrogenase I from Clostridium acetobutylicum (CaI) is a member of a large family of structurally related enzymes where interfacial electron transfer is mediated by a terminal, non-canonical, His-coordinated, [4Fe–4S] cluster. The function of His coordination was examined by comparing the biophysical properties and reactivity to a Cys substituted variant of CaI. This demonstrated that His coordination strongly affected the distal [4Fe–4S] cluster spin state, spin pairing, and spatial orientations of molecular orbitals, with a minor effect on reduction potential. The deviations in these properties by substituting His for Cys in CaI, correlated with pronounced changes in electron transfer and reactivity with the native electron donor–acceptor ferredoxin. The results demonstrate that differential coordination of the surface localized [4Fe–4S]His cluster in CaI is utilized to control intermolecular and intramolecular electron transfer where His coordination creates a physical and electronic environment that enables facile electron exchange between electron carrier molecules and the iron–sulfur cluster relay for coupling to reversible H₂ activation at the catalytic site.

Histidine coordination of the distal [4Fe–4S] cluster in [FeFe]-hydrogenase was demonstrated to tune the cluster spin-states, spin-pairing and surrounding molecular orbitals to enable more facile electron transfer compared to cysteine coordination.

The cytochrome b lysine 329 residue is critical for ubihydroquinone oxidation and proton release at the Qo site of bacterial cytochrome bc1

The ubihydroquinone:cytochrome (cyt) c oxidoreductase (or cyt bc₁) is an important enzyme for photosynthesis and respiration. In bacteria like Rhodobacter capsulatus, this membrane complex has three subunits, the iron‑sulfur protein (ISP) with its Fe₂S₂ cluster, cyt c₁ and cyt b, forming two catalytic domains, the Q_o (hydroquinone (QH₂) oxidation) and Q_i (quinone (Q) reduction) sites. At the Q_o site, the electron transfer pathways originating from QH₂ oxidation are known, but their associated proton release routes are less well defined. Earlier, we demonstrated that the His291 of cyt b is important for this latter process. In this work, using the bacterial cyt bc₁ and site directed mutagenesis, we show that Lys329 of cyt b is also critical for electron and proton transfer at the Q_o site. Of the mutants examined, Lys329Arg was photosynthesis proficient and had quasi-wild type cyt bc₁ activity. In contrast, the Lys329Ala and Lys329Asp were photosynthesis-impaired and contained defective but assembled cyt bc₁. In particular, the bifurcated electron transfer and associated proton(s) release reactions occurring during QH₂ oxidation were drastically impaired in Lys329Asp mutant. Furthermore, in silico docking studies showed that in this mutant the location and the H-bonding network around the Fe₂S₂ cluster of ISP on cyt b surface was different than the wild type enzyme. Based on these experimental findings and theoretical considerations, we propose that the presence of a positive charge at position 329 of cyt b is critical for efficient electron transfer and proton release for QH₂ oxidation at the Q_o site of cyt bc₁.

Importance of Hydrogen Bonding in Fine Tuning the [2Fe-2S] Cluster Redox Potential of HydC from Thermotoga maritima

Iron-sulfur clusters form one of the largest and most diverse classes of enzyme cofactors in nature. They may serve as structural factors, form electron transfer chains between active sites and external redox partners, or form components of enzyme active sites. Their specific role is a consequence of the cluster type and the surrounding protein environment. The relative effects of these factors are not completely understood, and it is not yet possible to predict the properties of iron-sulfur clusters based on amino acid sequences or rationally tune their properties to generate proteins with new desirable functions. Here, we generate mutations in a [2Fe-2S] cluster protein, the TmHydC subunit of the trimeric [FeFe]-hydrogenase from Thermotoga maritima, to study the factors that affect its redox potential. Saturation mutagenesis of Val131 was used to tune the redox potential over a 135 mV range and revealed that cluster redox potential and electronic properties correlate with amino acid hydrophobicity and the ability to form hydrogen bonds to the cluster. Proline scanning mutagenesis between pairs of ligating cysteines was used to remove backbone amide hydrogen bonds to the cluster and decrease the redox potential by up to 132 mV, without large structural changes in most cases. However, substitution of Gly83 with proline caused a change of HydC to a [4Fe-4S] cluster protein with a redox potential of -526 mV. Together, these results confirm the importance of hydrogen bonding in tuning cluster redox potentials and demonstrate the versatility of iron-sulfur cluster protein folds at binding different types of clusters.

The reduction rates of DEPC-modified mutant Thermus thermophilus Rieske proteins differ when there is a negative charge proximal to the cluster

Rieske and Rieske-type proteins are electron transport proteins involved in key biological processes such as respiration, photosynthesis, and detoxification. They have a [2Fe-2S] cluster ligated by two cysteines and two histidines. A series of mutations, L135E, L135R, L135A, and Y158F, of the Rieske protein from Thermus thermophilus has been produced which probe the effects of the neighboring residues, in the second sphere, on the dynamics of cluster reduction and the reactivity of the ligating histidines. These properties were probed using titrations and modifications with diethyl pyrocarbonate (DEPC) at various pH values monitored using UV-Visible and circular dichroism spectrophotometry. These results, along with results from EPR studies, provide information on ligating histidine modification and rate of reduction of each of the mutant proteins. L135R, L135A, and Y158F react with DEPC similarly to wild type, resulting in modified protein with a reduced [2Fe-2S] cluster in <90 min, whereas L135E requires >15 h under the same conditions. Thus, the negative charge slows down the rate of reduction and provides an explanation as to why negatively charged residues are rarely, if ever, found in the equivalent position of other Rieske and Rieske-type proteins.

Fast proton-coupled electron transfer observed for a high-fidelity structural and functional [2Fe-2S] Rieske model

Rieske cofactors have a [2Fe–2S] cluster with unique {His₂Cys₂} ligation and distinct Fe subsites. The histidine ligands are functionally relevant, since they allow for coupling of electron and proton transfer (PCET) during quinol oxidation in respiratory and photosynthetic ET chains. Here we present the highest fidelity synthetic analogue for the Rieske [2Fe–2S] cluster reported so far. This synthetic analogue 5^x– emulates the heteroleptic {His₂Cys₂} ligation of the [2Fe–2S] core, and it also serves as a functional model that undergoes fast concerted proton and electron transfer (CPET) upon reaction of the mixed-valent (ferrous/ferric) protonated 5H^2– with TEMPO. The thermodynamics of the PCET square scheme for 5^x– have been determined, and three species (diferric 5^2–, protonated diferric 5H^–, and mixed-valent 5^3–) have been characterized by X-ray diffraction. pK_a values for 5H^– and 5H^2– differ by about 4 units, and the reduction potential of 5H^– is shifted anodically by about +230 mV compared to that of 5^2–. While the N–H bond dissociation free energy of 5H^2– (60.2 ± 0.5 kcal mol^–1) and the free energy, ΔG°_CPET, of its reaction with TEMPO (−6.3 kcal mol^–1) are similar to values recently reported for a homoleptic {N₂/N₂}-coordinated [2Fe–2S] cluster, CPET is significantly faster for 5H^2– with biomimetic {N₂/S₂} ligation (k = (9.5 ± 1.2) × 10⁴ M^–1 s^–1, ΔH^‡ = 8.7 ± 1.0 kJ mol^–1, ΔS^‡ = −120 ± 40 J mol^–1 K^–1, and ΔG^‡ = 43.8 ± 0.3 kJ mol^–1 at 293 K). These parameters, and the comparison with homoleptic analogues, provide important information and new perspectives for the mechanistic understanding of the biological Rieske cofactor.

Dissection of hydrogen bond interaction network around an iron-sulfur cluster by site-specific isotope labeling of hyperthermophilic archaeal Rieske-type ferredoxin

The electronic structure and geometry of redox-active metal cofactors in proteins are tuned by the pattern of hydrogen bonding with the backbone peptide matrix. In this study we developed a method for selective amino acid labeling of a hyperthermophilic archaeal metalloprotein with engineered Escherichia coli auxotroph strains, and we applied this to resolve the hydrogen bond interactions with the reduced Rieske-type [2Fe-2S] cluster by two-dimensional pulsed electron spin resonance technique. Because deep electron spin-echo envelope modulation of two histidine (14)N(δ) ligands of the cluster decreased non-coordinating (15)N signal intensities via the cross-suppression effect, an inverse labeling strategy was employed in which (14)N amino acid-labeled archaeal Rieske-type ferredoxin samples were examined in an (15)N-protein background. This has directly identified Lys45 N(α) as providing the major pathway for the transfer of unpaired electron spin density from the reduced cluster by a "through-bond" mechanism. All other backbone peptide nitrogens interact more weakly with the reduced cluster. The extension of this approach will allow visualizing the three-dimensional landscape of preferred pathways for the transfer of unpaired spin density from a paramagnetic metal center onto the protein frame, and will discriminate specific interactions by a "through-bond" mechanism from interactions which are "through-space" in various metalloproteins.

Buffering capacity and size of siRNA polyplexes influence cytokine levels

Induction of cytokines by small interfering RNA (siRNA) polyplexes has been a significant concern of researchers attempting to minimize the toxicity of this promising therapy. Although cationic carriers of siRNA are known to increase cytokine levels, few systematic studies have been done to determine what properties of the carrier are important to modulate cytokines. Because branched histidine-lysine (HK) peptides are effective carriers of siRNA and their sequence can be readily modified, we selected this class of carrier to determine which sequences of the peptide were important for cytokine induction. With the use of peripheral blood mononuclear cells (PBMCs), the HK peptide with a higher number of histidines (H3K(+H)4b) in complex with siRNA induced lower levels of cytokines compared with other HK (e.g., H2K4b, H3K4b, H3K(+N)4b) siRNA nanoplexes. Notably, these peptides' siRNA polyplexes showed a similar pattern of cytokine induction when injected intravenously in a mouse model, i.e., the HK with higher content of histidines induced cytokines the least. As indicated by the pH-sensitive dye within acidic endosomes, the greater pH-buffering capacity of H3K(+H)4b compared with other HK peptides may explain why cytokine levels were reduced. In addition to buffering capacity, the size of HK polyplexes markedly influenced cytokine production.

Reduction potentials and acidity constants of Mn superoxide dismutase calculated by QM/MM free-energy methods

We used two theoretical methods to estimate reduction potentials and acidity constants in Mn superoxide dismutase (MnSOD), namely combined quantum mechanical and molecular mechanics (QM/MM) thermodynamic cycle perturbation (QTCP) and the QM/MM-PBSA approach. In the latter, QM/MM energies are combined with continuum solvation energies calculated by solving the Poisson-Boltzmann equation (PB) or by the generalised Born approach (GB) and non-polar solvation energies calculated from the solvent-exposed surface area. We show that using the QTCP method, we can obtain accurate and precise estimates of the proton-coupled reduction potential for MnSOD, 0.30±0.01 V, which compares favourably with experimental estimates of 0.26-0.40 V. However, the calculated potentials depend strongly on the DFT functional used: The B3LYP functional gives 0.6 V more positive potentials than the PBE functional. The QM/MM-PBSA approach leads to somewhat too high reduction potentials for the coupled reaction and the results depend on the solvation model used. For reactions involving a change in the net charge of the metal site, the corresponding results differ by up to 1.3 V or 24 pK(a) units, rendering the QM/MM-PBSA method useless to determine absolute potentials. However, it may still be useful to estimate relative shifts, although the QTCP method is expected to be more accurate.

A combined far-infrared spectroscopic and electrochemical approach for the study of iron-sulfur proteins

Herein, we present the development of a far-infrared spectroscopic approach for studying metalloenzyme active sites in a redox-dependent manner. An electrochemical cell with 5 mm path and based on silicon windows was found to be appropriate for the measurement of aqueous solutions down to 200 cm(-1) . The cell was probed with the infrared redox signature of the metal-ligand vibrations of different iron-sulfur proteins. Each Fe-S cluster type was found to show a specific spectral signature. As a common feature, a downshift of the frequency of the Fe-S vibrations was seen upon reduction, in line with the increase of the Fe-S bond. This downshift was found to be fully reversible. Electrochemically induced FTIR difference spectroscopy in the far infrared is now possible, opening new perspectives on the understanding of metalloproteins in function of the redox state.

Engineering the redox potential over a wide range within a new class of FeS proteins

MitoNEET is a newly discovered mitochondrial protein and a target of the TZD class of antidiabetes drugs. MitoNEET is homodimeric with each protomer binding a [2Fe-2S] center through a rare 3-Cys and 1-His coordination geometry. Both the fold and the coordination of the [2Fe-2S] centers suggest that it could have novel properties compared to other known [2Fe-2S] proteins. We tested the robustness of mitoNEET to mutation and the range over which the redox potential (E(M)) could be tuned. We found that the protein could tolerate an array of mutations that modified the E(M) of the [2Fe-2S] center over a range of ∼700 mV, which is the largest E(M) range engineered in an FeS protein and, importantly, spans the cellular redox range (+200 to -300 mV). These properties make mitoNEET potentially useful for both physiological studies and industrial applications as a stable, water-soluble, redox agent.

Redox potential of the Rieske iron-sulfur protein quantum-chemical and electrostatic study

Quantum-chemical study of structures, energies, and effective partial charge distribution for several models of the Rieske protein redox center is performed in terms of the B3LYP density functional method in combination with the broken symmetry approach using three different atomic basis sets. The structure of the redox complex optimized in vacuum differs markedly from that inside the protein. This means that the protein matrix imposes some stress on the active site resulting in distortion of its structure. The redox potentials calculated for the real active site structure are in a substantially better agreement with the experiment than those calculated for the idealized structure. This shows an important role of the active site distortion in tuning its redox potential. The reference absolute electrode potential of the standard hydrogen electrode is used that accounts for the correction caused by the water surface potential. Electrostatic calculations are performed in the framework of the polarizable solute model. Two dielectric permittivities of the protein are employed: the optical permittivity for calculation of the intraprotein electric field, and the static permittivity for calculation of the dielectric response energy. Only this approach results in a reasonable agreement of the calculated and experimental redox potentials.

Modifications of protein environment of the [2Fe-2S] cluster of the bc1 complex: effects on the biophysical properties of the rieske iron-sulfur protein and on the kinetics of the complex

The rate-determining step in the overall turnover of the bc(1) complex is electron transfer from ubiquinol to the Rieske iron-sulfur protein (ISP) at the Q(o)-site. Structures of the ISP from Rhodobacter sphaeroides show that serine 154 and tyrosine 156 form H-bonds to S-1 of the [2Fe-2S] cluster and to the sulfur atom of the cysteine liganding Fe-1 of the cluster, respectively. These are responsible in part for the high potential (E(m)(,7) approximately 300 mV) and low pK(a) (7.6) of the ISP, which determine the overall reaction rate of the bc(1) complex. We have made site-directed mutations at these residues, measured thermodynamic properties using protein film voltammetry to evaluate the E(m) and pK(a) values of ISPs, explored the local proton environment through two-dimensional electron spin echo envelope modulation, and characterized function in strains S154T, S154C, S154A, Y156F, and Y156W. Alterations in reaction rate were investigated under conditions in which concentration of one substrate (ubiquinol or ISP(ox)) was saturating and the other was varied, allowing calculation of kinetic terms and relative affinities. These studies confirm that H-bonds to the cluster or its ligands are important determinants of the electrochemical characteristics of the ISP, likely through electron affinity of the interacting atom and the geometry of the H-bonding neighborhood. The calculated parameters were used in a detailed Marcus-Brønsted analysis of the dependence of rate on driving force and pH. The proton-first-then-electron model proposed accounts naturally for the effects of mutation on the overall reaction.

Two-dimensional pulsed electron spin resonance characterization of 15N-labeled archaeal Rieske-type ferredoxin

Two-dimensional electron spin-echo envelope modulation (ESEEM) analysis of the uniformly (15)N-labeled archaeal Rieske-type [2Fe-2S] ferredoxin (ARF) from Sulfolobus solfataricus P1 has been conducted in comparison with the previously characterized high-potential protein homologs. Major differences among these proteins were found in the hyperfine sublevel correlation (HYSCORE) lineshapes and intensities of the signals in the (++) quadrant, which are contributed from weakly coupled (non-coordinated) peptide nitrogens near the reduced clusters. They are less pronounced in the HYSCORE spectra of ARF than those of the high-potential protein homologs, and may account for the tuning of Rieske-type clusters in various redox systems.

M234Glu is a component of the proton sponge in the reaction center from photosynthetic bacteria

Bacterial reaction centers use light energy to couple the uptake of protons to the successive semi-reduction of two quinones, namely Q(A) and Q(B). These molecules are situated symmetrically in regard to a non-heme iron atom. Four histidines and one glutamic acid, M234Glu, constitute the five ligands of this atom. By flash-induced absorption spectroscopy and delayed fluorescence we have studied in the M234EH and M234EL variants the role played by this acidic residue on the energetic balance between the two quinones as well as in proton uptake. Delayed fluorescence from the P(+)Q(A)(-) state (P is the primary electron donor) and temperature dependence of the rate of P(+)Q(A)(-) charge recombination that are in good agreement show that in the two RC variants, both Q(A)(-) and Q(B)(-) are destabilized by about the same free energy amount: respectively approximately 100 +/- 5 meV and 90 +/- 5 meV for the M234EH and M234EL variants, as compared to the WT. Importantly, in the M234EH and M234EL variants we observe a collapse of the high pH band (present in the wild-type reaction center) of the proton uptake amplitudes associated with formation of Q(A)(-) and Q(B)(-). This band has recently been shown to be a signature of a collective behaviour of an extended, multi-entry, proton uptake network. M234Glu seems to play a central role in the proton sponge-like system formed by the RC protein.

Resonance Raman studies of the (His)(Cys)3 2Fe-2S cluster of MitoNEET: comparison to the (Cys)4 mutant and implications of the effects of pH on the labile metal center

MitoNEET is a 2Fe-2S outer mitochondrial membrane protein that was initially identified as a target for anti-diabetic drugs. It exhibits a novel protein fold, and in contrast to other 2Fe-2S proteins such as Rieske proteins and ferredoxins, the metal clusters in the mitoNEET homodimer are each coordinated by one histidine residue and three cysteine residues. The interaction of the ligating His87 residue with the 2Fe-2S moiety is especially significant because previous studies have shown that replacement with Cys in the H87C mutant stabilizes the cluster against release. Here, we report the resonance Raman spectra of this naturally occurring Fe(2)S(2)(His)(Cys)(3) protein to assess local structural changes associated with cluster lability. Comparison of mitoNEET to its ferredoxin-like H87C mutant indicates that Raman peaks in the approximately 250-300 cm(-1) region of mitoNEET are influenced by the Fe-His87 moiety. Systematic pH-dependent resonance Raman spectral changes were observed in this spectral region for native mitoNEET but not the H87C mutant. The approximately 250-300 cm(-1) region of native mitoNEET is also sensitive to phosphate buffer. Thus, conditions that influence cluster release are shown here to concomitantly affect the resonance Raman spectrum in the region with Fe-His contribution. These results support the hypothesis that the Fe-N(His87) interaction is modulated within the physiological pH range, and this modulation may be critical to the function of mitoNEET.

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.

Which factors determine the acidity of the phytochromobilin chromophore of plant phytochrome?

Quantum chemical calculations aimed at identifying the factors controlling the acidity of phytochromobilin, the tetrapyrrole chromophore of the plant photoreceptor phytochrome, are reported. Phytochrome is converted from an inactive (Pr) to an active form (Pfr) through a series of events initiated by a Z --> E photoisomerization of phytochromobilin, forming the Lumi-R intermediate, and much controversy exists as to whether the protonation state of the chromophore (cationic in Pr with all nitrogens protonated) changes during the photoactivation. Here, relative ground (S0) and excited-state (S1) pKa s of all four pyrrole moieties of phytochromobilin in all 64 possible configurations with respect to the three methine bridges are calculated in a protein-like environment, using a recently benchmarked level of theory. Accordingly, the relationships between acidity and chromophore geometry and charge distribution, hydrogen bonding, and light absorption are investigated in some detail, and discussed in terms of possible mechanisms making a proton transfer reaction more probable along the Pr --> Pfr reaction than in the parent cationic Pr state. It is found that charge distribution in the cationic species, intra-molecular hydrogen bonding in the neutral, and hydrogen bonding with two highly conserved aspartate and histidine residues have a significant effect on the acidity, while overall chromophore geometry and electronic state are less important factors. Furthermore, based on the calculations, two processes that may facilitate a proton transfer by substantially lowering the pKa s relative to their Pr values are identified: (i) a thermal Z,anti --> Z,syn isomerization at C5, occurring after formation of Lumi-R; (ii) a perturbation of the hydrogen bonding network which in Pr comprises the nitrogens of pyrroles A, B and C and the two aspartate and histidine residues.

Relative Ground and Excited-State pKa Values of Phytochromobilin in the Photoactivation of Phytochrome: A Computational Study

The conversion of the plant photoreceptor phytochrome from an inactive (Pr) to an active form (Pfr) is accomplished by a red-light induced Z --> E photoisomerization of its phytochromobilin chromophore. In recent years, the question whether the photoactivation involves a change in chromophore protonation state has been the subject of many experimental studies. Here, we have used quantum chemical methods to calculate relative ground and excited-state pKa values of the different pyrrole moieties of phytochromobilin in a protein-like environment. Assuming (based on experimental data) a Pr ZaZsZa chromophore and considering isomerizations at C15 and C5, it is found that moieties B and C are the strongest acids both in the ground state and in the bright first singlet excited state, which is rationalized in simple geometric and electronic terms. It is also shown that neither light absorption nor isomerization increases the acidity of phytochromobilin relative to the reference Pr state with all pyrrolenic nitrogens protonated. Hence, provided that the subset of chromophore geometries under investigation is biologically relevant, there appears to be no intrinsic driving force for a proton-transfer event. In a series of benchmark calculations, the performance of ab initio and time-dependent density functional theory methods for excited-state studies of phytochromobilin is evaluated in light of available experimental data.

Rieske protein from Thermus thermophilus: 15N NMR titration study demonstrates the role of iron-ligated histidines in the pH dependence of the reduction potential

A unique feature of Rieske proteins is the pH dependence of their reduction potentials. It has been proposed that protonation of the Nepsilon2 atoms of the two histidine rings ligated to the iron-sulfur cluster is coupled with cluster reduction (electron transfer). We have incorporated [15Ndelta1, 15Nepsilon2]-histidine into the Rieske protein from Thermus thermophilis and have used 15N NMR spectroscopy to determine the pKa values of the histidine residues in the oxidized state of the protein. As expected from studies of a Rieske-type ferredoxin, the signals from the 15Ndelta1 atoms directly bound to iron were too broad to be detected, but broad signals could be detected from the 15Nepsilon2 atom of each of the ligated histidine rings. We measured the chemical shifts of these signals as a function of pH between pH 6 and pH 12 and fitted them to theoretical titration curves. The results yielded well-separated pKa values for the two histidines (7.46 and 9.24), with Hill coefficients close to unity. The pKa values are in excellent agreement with values predicted from the pH dependence of the reduction potentials (7.85 and 9.65).

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.

Development, validation, and application of adapted PEOE charges to estimate pKa values of functional groups in protein-ligand complexes

For routine pK(a) calculations of protein-ligand complexes in drug design, the PEOE method to compute partial charges was modified. The new method is applicable to a large scope of proteins and ligands. The adapted charges were parameterized using experimental free energies of solvation of amino acids and small organic ligands. For a data set of 80 small organic molecules, a correlation coefficient of r(2) = 0.78 between calculated and experimental solvation free energies was obtained. Continuum electrostatics pK(a) calculations based on the Poisson-Boltzmann equation were carried out on a validation set of nine proteins for which 132 experimental pK(a) values are known. In total, an overall RMSD of 0.88 log units between calculated and experimentally determined data is achieved. In particular, the predictions of significantly shifted pK(a) values are satisfactory, and reasonable estimates of protonation states in the active sites of lysozyme and xylanase could be obtained. Application of the charge-assignment and pK(a)-calculation procedure to protein-ligand complexes provides clear structural interpretations of experimentally observed changes of protonation states of functional groups upon complex formation. This information is essential for the interpretation of thermodynamic data of protein-ligand complex formation and provides the basis for the reliable factorization of the free energy of binding in enthalpic and entropic contributions. The modified charge-assignment procedure forms the basis for future automated pK(a) calculations of protein-ligand complexes.

Identification of Hydrogen Bonds to the Rieske Cluster through the Weakly Coupled Nitrogens Detected by Electron Spin Echo Envelope Modulation Spectroscopy

The interaction of the reduced[2Fe-2S] cluster of isolated Rieske fragment from the bc1 complex of Rhodobacter sphaeroides with nitrogens (14N and 15N) from the local protein environment has been studied by X- and S-band pulsed EPR spectroscopy. The two-dimensional electron spin echo envelope modulation spectra of uniformly 15N-labeled protein show two well resolved cross-peaks with weak couplings of approximately 0.3-0.4 and 1.1 MHz in addition to couplings in the range of 6-8 MHz from two coordinating Ndelta of histidine ligands. The quadrupole coupling constants for weakly coupled nitrogens determined from S-band electron spin echo envelope modulation spectra identify them as Nepsilon of histidine ligands and peptide nitrogen (Np), respectively. Analysis of the line intensities in orientation-selected S-band spectra indicated that Np is the backbone N-atom of Leu-132 residue. The hyperfine couplings from Nepsilon and Np demonstrate the predominantly isotropic character resulting from the transfer of unpaired spin density onto the 2s orbitals of the nitrogens. Spectra also show that other peptide nitrogens in the protein environment must carry a 5-10 times smaller amount of spin density than the Np of Leu-132 residue. The appearance of the excess unpaired spin density on the Np of Leu-132 residue indicates its involvement in hydrogen bond formation with the bridging sulfur of the Rieske cluster. The configuration of the hydrogen bond therefore provides a preferred path for spin density transfer. Observation of similar splittings in the 15N spectra of other Rieske-type proteins and [2Fe-2S] ferredoxins suggests that a hydrogen bond between the bridging sulfur and peptide nitrogen is a common structural feature of [2Fe-2S] clusters.

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.

2-Oxoquinoline 8-monooxygenase oxygenase component: active site modulation by Rieske-[2Fe-2S] center oxidation/reduction

2-Oxoquinoline 8-monooxygenase is a Rieske non-heme iron oxygenase that catalyzes the NADH-dependent oxidation of the N-heterocyclic aromatic compound 2-oxoquinoline to 8-hydroxy-2-oxoquinoline in the soil bacterium Pseudomonas putida 86. The crystal structure of the oxygenase component of 2-oxoquinoline 8-monooxygenase shows a ring-shaped, C3-symmetric arrangement in which the mononuclear Fe(II) ion active site of one monomer is at a distance of 13 A from the Rieske-[2Fe-2S] center of a second monomer. Structural analyses of oxidized, reduced, and substrate bound states reveal the molecular bases for a new function of Fe-S clusters. Reduction of the Rieske center modulates the mononuclear Fe through a chain of conformational changes across the subunit interface, resulting in the displacement of Fe and its histidine ligand away from the substrate binding site. This creates an additional coordination site at the mononuclear Fe(II) ion and can open a pathway for dioxygen to bind in the substrate-containing active site.

A Rieske ferredoxin typifying a subtype within Rieske proteins: spectroscopic, biochemical and stability studies

A new subtype of archaeal Rieske ferredoxin (RFd) has been identified in the genome of the thermoacidophilic archaeon Acidianus ambivalens. The gene is inserted in an atypical genomic context in a gene cluster encoding a NiFe hydrogenase. Sequence and phyletic analysis showed that the protein is related to bacterial RFd but not to any of the known archaeal Rieske proteins. The recombinant 14 kDa protein isolated from Escherichia coli behaved as a dimer in solution. It contained approximately 2 Fe/mol and all visible and EPR spectroscopic features typical of Rieske centre-containing proteins. However, its redox potential (+170 mV) was significantly higher than those of canonical RFd. This difference is rationalized in terms of the protein structure environment, as discrete amino acid substitutions in key positions around the metal centre account for the higher potential.