Probing the Structure of [NiFeSe] Hydrogenase with QM/MM Computations

The geometry and vibrational behavior of selenocysteine [NiFeSe] hydrogenase isolated from Desulfovibrio vulgaris Hildenborough have been investigated using a hybrid quantum mechanical (QM)/ molecular mechanical (MM) approach. Structural models have been built based on the three conformers identified in the recent crystal structure resolved at 1.3 Å from X-ray crystallography. In the models, a diamagnetic Ni2+ atom was modeled in combination with both Fe2+ and Fe3+ to investigate the effect of iron oxidation on geometry and vibrational frequency of the nonproteic ligands, CO and CN-, coordinated to the Fe atom. Overall, the QM/MM optimized geometries are in good agreement with the experimentally resolved geometries. Analysis of computed vibrational frequencies, in comparison with experimental Fourier-transform infrared (FTIR) frequencies, suggests that a mixture of conformers as well as Fe2+ and Fe3+ oxidation states may be responsible for the acquired vibrational spectra.DFG, 390540038, EXC 2008: UniSysCatDFG, 414044773, Open Access Publizieren 2019 - 2020 / Technische Universität BerlinEC/H2020/810856/EU/Twin to Illuminate Metals in Biology and Biocatalysis through Biospectroscopy/TIMB

Structural and electronic properties of the active site of [ZnFe] SulE

The function of the recently isolated sulerythrin (SulE) has been investigated using a combination of structural and electronic analyses based on quantum mechanical calculations. In the SulE structure of Fushinobu et al. (2003), isolated from a strictly aerobic archaeon, Sulfolobus tokadaii, a dioxygen-containing species was tentatively included at the active site during crystallographic refinement although the substrate specificity of SulE remains unclear. Studies have suggested that a structurally related enzyme, rubrerythrin, functions as a hydrogen peroxide reductase. Since SulE is a truncated version of rubrerythrin, the enzymes are hypothesized to function similarly. Hence, using available X-ray crystallography data (1.7 Å), we constructed various models of SulE containing a ZnII–Fe active site, differing in the nature of the substrate specificity (O2, H2O2), the oxidation level and the spin state of the iron ion, and the protonation states of the coordinating glutamate residues. Also, the substrate H2O2 is modeled in two possible configurations, differing in the orientation of the hydrogen atoms. Overall, the optimized geometries with an O2 substrate do not show good agreement with the experimentally resolved geometry. In contrast, excellent agreement between crystal structure arrangement and optimized geometries is achieved considering a H2O2 substrate and FeII in both spin states, when Glu92 is protonated. These results suggest that the dioxo species detected at the [ZnFe] active site of sulerythrin is H2O2, rather than an O2 molecule in agreement with experimental data indicating that only the diferrous oxidation state of the dimetal site in rubrerythrin reacts rapidly with H2O2. Based on our computations, we proposed a possible reaction pathway for substrate binding at the ZnFeII site of SulE with a H2O2 substrate. In this reaction pathway, Fe or another electron donor, such as NAD(P)H, catalyzes the reduction of H2O2 to water at the zinc–iron site

Insights into the structure of the active site of the O-2-tolerant membrane bound [NiFe] hydrogenase of R. eutropha H16 by molecular modelling

Dieser Beitrag ist mit Zustimmung des Rechteinhabers aufgrund einer (DFG geförderten) Allianz- bzw. Nationallizenz frei zugänglich.This publication is with permission of the rights owner freely accessible due to an Alliance licence and a national licence (funded by the DFG, German Research Foundation) respectively.Structural models for the Ni-B state of the wild-type and C81S protein variant of the membrane-bound [NiFe] hydrogenase from Ralstonia eutrophaH16 were derived by applying the homology model technique combined with molecular simulations and a hybrid quantum mechanical/molecular mechanical approach. The active site structure was assessed by comparing calculated and experimental IR spectra, confirming the view that the active site structure is very similar to those of anaerobic standard hydrogenases. In addition, the data suggest the presence of a water molecule in the second coordination sphere of the active centre.DFG, EXC 314, Unifying Concepts in Catalysi

Ultrafast proton release reaction and primary photochemistry of phycocyanobilin in solution observed with fs-time-resolved mid-IR and UV/Vis spectroscopy

Deactivation processes of photoexcited (λ$_{ex}$ = 580 nm) phycocyanobilin (PCB) in methanol were investigated by means of UV/Vis and mid-IR femtosecond (fs) transient absorption (TA) as well as static fluorescence spectroscopy, supported by density-functional-theory calculations of three relevant ground state conformers, PCB$_{A}$, PCB$_{B}$ and PCB$_{C}$, their relative electronic state energies and normal mode vibrational analysis. UV/Vis fs-TA reveals time constants of 2.0, 18 and 67 ps, describing decay of PCB$_{B}$*, of PCB$_{A}$* and thermal re-equilibration of PCB$_{A}$, PCB$_{B}$ and PCB$_{C}$, respectively, in line with the model by Dietzek et al. (Chem Phys Lett 515:163, 2011) and predecessors. Significant substantiation and extension of this model is achieved first via mid-IR fs-TA, i.e. identification of molecular structures and their dynamics, with time constants of 2.6, 21 and 40 ps, respectively. Second, transient IR continuum absorption (CA) is observed in the region above 1755 cm$^{-1}$ (CA1) and between 1550 and 1450 cm$^{-1}$ (CA2), indicative for the IR absorption of highly polarizable protons in hydrogen bonding networks (X–H…Y). This allows to characterize chromophore protonation/deprotonation processes, associated with the electronic and structural dynamics, on a molecular level. The PCB photocycle is suggested to be closed via a long living (> 1 ns), PCBC-like (i.e. deprotonated), fluorescent species

Ultrafast protein response in the Pfr state of Cph1 phytochrome

Photoisomerization is a fundamental process in several classes of photoreceptors. Phytochromes sense red and far-red light in their Pr and Pfr states, respectively. Upon light absorption, these states react via individual photoreactions to the other state. Cph1 phytochrome shows a photoisomerization of its phycocyanobilin (PCB) chromophore in the Pfr state with a time constant of 0.7 ps. The dynamics of the PCB chromophore has been described, but whether or not the apoprotein exhibits an ultrafast response too, is not known. Here, we compare the photoreaction of 13C/15N labeled apoprotein with unlabeled apoprotein to unravel ultrafast apoprotein dynamics in Cph1. In the spectral range from 1750 to 1620 cm−1 we assigned several signals due to ultrafast apoprotein dynamics. A bleaching signal at 1724 cm−1 is tentatively assigned to deprotonation of a carboxylic acid, probably Asp207, and signals around 1670 cm−1 are assigned to amide I vibrations of the capping helix close to the chromophore. These signals remain after photoisomerization. The apoprotein dynamics appear upon photoexcitation or concomitant with chromophore isomerization. Thus, apoprotein dynamics occur prior to and after photoisomerization on an ultrafast time-scale. We discuss the origin of the ultrafast apoprotein response with the ‘Coulomb hammer’ mechanism, i.e. an impulsive change of electric field and Coulombic force around the chromophore upon excitation

Understanding flavin electronic structure and spectra

Flavins have emerged as central to electron bifurcation, signaling, and countless enzymatic reactions. In bifurcation, two electrons acquired as a pair are separated in coupled transfers wherein the energy of both is concentrated on one of the two. This enables organisms to drive demanding reactions based on abundant low-grade chemical fuel. To enable incorporation of this and other flavin capabilities into designed materials and devices, it is essential to understand fundamental principles of flavin electronic structure that make flavins so reactive and tunable by interactions with protein. Emerging computational tools can now replicate spectra of flavins and are gaining capacity to explain reactivity at atomistic resolution, based on electronic structures. Such fundamental understanding can moreover be transferrable to other chemical systems. A variety of computational innovations have been critical in reproducing experimental properties of flavins including their electronic spectra, vibrational signatures, and nuclear magnetic resonance (NMR) chemical shifts. A computational toolbox for understanding flavin reactivity moreover must be able to treat all five oxidation and protonation states, in addition to excited states that participate in flavoprotein's light-driven reactions. Therefore, we compare emerging hybrid strategies and their successes in replicating effects of hydrogen bonding, the surrounding dielectric, and local electrostatics. These contribute to the protein's ability to modulate flavin reactivity, so we conclude with a survey of methods for incorporating the effects of the protein residues explicitly, as well as local dynamics. Computation is poised to elucidate the factors that affect a bound flavin's ability to mediate stunningly diverse reactions, and make life possible.DFG, 390540038, EXC 2008: Unifying Systems in Catalysis "UniSysCat

Insights from surface enhanced resonance Raman spectroscopy and QM/MM calculations

Understanding the coupling between heme reduction and proton translocation in cytochrome c oxidase (CcO) is still an open problem. The propionic acids of heme a3 have been proposed to act as a proton loading site (PLS) in the proton pumping pathway, yet this proposal could not be verified by experimental data so far. We have set up an experiment where the redox states of the two hemes in CcO can be controlled via external electrical potential. Surface enhanced resonance Raman (SERR) spectroscopy was applied to simultaneously monitor the redox state of the hemes and the protonation state of the heme propionates. Simulated spectra based on QM/MM calculations were used to assign the resonant enhanced CH2 twisting modes of the propionates to the protonation state of the individual heme a and heme a3 propionates respectively. The comparison between calculated and measured H2OD2O difference spectra allowed a sound band assignment. In the fully reduced enzyme at least three of the four heme propionates were found to be protonated whereas in the presence of a reduced heme a and an oxidized heme a3 only protonation of one heme a3 propionates was observed. Our data supports the postulated scenario where the heme a3 propionates are involved in the proton pathway

Characterisation of the Cyanate Inhibited State of Cytochrome c Oxidase

Heme-copper oxygen reductases are terminal respiratory enzymes, catalyzing the reduction of dioxygen to water and the translocation of protons across the membrane. Oxygen consumption is inhibited by various substances. Here we tested the relatively unknown inhibition of cytochrome c oxidase (CcO) with isocyanate. In contrast to other more common inhibitors like cyanide, inhibition with cyanate was accompanied with the rise of a metal to ligand charge transfer (MLCT) band around 638nm. Increasing the cyanate concentration furthermore caused selective reduction of heme a. The presence of the CT band allowed for the first time to directly monitor the nature of the ligand via surface-enhanced resonance Raman (SERR) spectroscopy. Analysis of isotope sensitive SERR spectra in comparison with Density Functional Theory (DFT) calculations identified not only the cyanate monomer as an inhibiting ligand but suggested also presence of an uretdion ligand formed upon dimerization of two cyanate ions. It is therefore proposed that under high cyanate concentrations the catalytic site of CcO promotes cyanate dimerization. The two excess electrons that are supplied from the uretdion ligand lead to the observed physiologically inverse electron transfer from heme a(3) to heme a

Domain motions and electron transfer dynamics in 2Fe-superoxide reductase

Superoxide reductases are non-heme iron enzymes that represent valuable model systems for the reductive detoxification of reactive oxygen species. In the present study, we applied different theoretical methods to study the structural dynamics of a prototypical 2Fe-superoxide reductase and its influence on electron transfer towards the active site. Using normal mode and essential dynamics analyses, we could show that enzymes of this type are capable of well-defined, electrostatically triggered domain movements, which may allow conformational proofreading for cellular redox partners involved in intermolecular electron transfer. Moreover, these global modes of motion were found to enable access to molecular configurations with decreased tunnelling distances between the active site and the enzyme's second iron centre. Using all-atom classical molecular dynamics simulations and the tunnelling pathway model, however, we found that electron transfer between the two metal sites is not accelerated under these conditions. This unexpected finding suggests that the unperturbed enzymatic structure is optimized for intramolecular electron transfer, which provides an indirect indication of the biological relevance of such a mechanism. Consistently, efficient electron transfer was found to depend on a distinct route, which is accessible via the equilibrium geometry and characterized by a quasi conserved tyrosine that could enable multistep-tunnelling (hopping). Besides these explicit findings, the present study demonstrates the importance of considering both global and local protein dynamics, and a generalized approach for the functional analysis of these aspects is provided

Halogenation of tyrosine perturbs large-scale protein self-organization

Protein halogenation is a common non-enzymatic post-translational modification contributing to aging, oxidative stress-related diseases and cancer. Here, we report a genetically encodable halogenation of tyrosine residues in a reconstituted prokaryotic filamentous cell-division protein (FtsZ) as a platform to elucidate the implications of halogenation that can be extrapolated to living systems of much higher complexity. We show how single halogenations can fine-tune protein structures and dynamics of FtsZ with subtle perturbations collectively amplified by the process of FtsZ self-organization. Based on experiments and theories, we have gained valuable insights into the mechanism of halogen influence. The bending of FtsZ structures occurs by affecting surface charges and internal domain distances and is reflected in the decline of GTPase activities by reducing GTP binding energy during polymerization. Our results point to a better understanding of the physiological and pathological effects of protein halogenation and may contribute to the development of potential diagnostic tools