Direct Measurement of the Hydrogen-Bonding Effect on the Intrinsic Redox Potentials of [4Fe−4S] Cubane Complexes

Structure of formylglycine-generating enzyme in complex with copper and a substrate reveals an acidic pocket for binding and activation of molecular oxygen

The formylglycine generating enzyme (FGE) catalyzes oxidative conversion of specific peptidyl-cysteine residues to formylglycine. FGE mediates O₂-activation and hydrogen-atom abstraction in an active site that contains Cu(i) coordinated to two cysteine residues. Similar coordination geometries are common among copper-sensing transcription factors and copper-chaperone but are unprecedented among copper-dependent oxidases. To examine the mechanism of this unusual catalyst we determined the 1.04 Å structure of FGE from Thermomonospora curvata in complex with copper and a cysteine-containing peptide substrate. This structure unveils a network of four crystallographic waters and two active site residues that form a highly acidic O₂-binding pocket juxtaposed to the trigonal planar tris-cysteine coordinated Cu(i) center. Comparison with structures of FGE in complex with Ag(i) and Cd(ii) combined with evidence from NMR spectroscopy and kinetic observations highlight several structural changes that are induced by substrate binding and prime the enzyme for O₂-binding and subsequent activation.

Fe-V sulfur clusters studied through photoelectron spectroscopy and density functional theory

Iron-vanadium sulfur cluster anions are studied by photoelectron spectroscopy (PES) at 3.492 eV (355 nm) and 4.661 eV (266 nm) photon energies, and by density functional theory (DFT) calculations. The structural properties, relative energies of different structural isomers, and the calculated first vertical detachment energies (VDEs) of different structural isomers for cluster anions FeVS1-3- and FemVnSm+n- (m + n = 3, 4; m > 0, n > 0) are investigated at a BPW91/TZVP theory level. The experimental first VDEs for these Fe-V sulfur clusters are reported. The most probable ground state structures and spin multiplicities for these clusters are tentatively assigned by comparing their theoretical and experiment first VDE values. For FeVS1-3- clusters, their first VDEs are generally observed to increase with the number of sulfur atoms from 1.45 eV to 2.86 eV. The NBO/HOMOs of the ground state of FeVS1-3- clusters are localized in a p orbital on a S atom; the partial charge distribution on the NBO/HOMO localized site of each cluster anion is responsible for the trend of their first VDEs. A less negative localized charge distribution is correlated with a higher first VDE. Structure and steric effect differences for FemVnSm+n- (m + n = 3, m > 0, n > 0) clusters are suggested to be responsible for their different first VDEs and properties. Two types of structural isomers are identified for FemVnSm+n- (m + n = 4, m > 0, n > 0) clusters: a tower structure isomer and a cubic structure isomer. The first VDEs for tower like isomers are generally higher than those for cubic like isomers of FemVnSm+n- (m + n = 4, m > 0, n > 0) clusters. Their first VDEs are can be understood through: (1) NBO/HOMO distributions, (2) structures (steric effects), and (3) partial charge numbers on the NBO/HOMO's localized sites. EBEs for excited state transitions for all Fe-V sulfur clusters are calculated employing OVGF and TDDFT approaches at the TZVP level. The OVGF approach for these Fe/V/S cluster anions is better for the higher transition energies than the TDDTF approach. The experimental and theoretical results for these Fe/V/S cluster anions are compared with their related pure iron sulfur cluster anions. Properties of the NBO/HOMO are essential for understanding and estimating the different first VDEs for Fe/V/S, and comparing them to those of the pure Fe/S cluster anions.

Photoelectron spectroscopy and density functional theory studies of (FeS)mH- (m = 2-4) cluster anions: effects of the single hydrogen

Single hydrogen containing iron hydrosulfide cluster anions (FeS)_mH^- (m = 2-4) are studied by photoelectron spectroscopy (PES) at 3.492 eV (355 nm) and 4.661 eV (266 nm) photon energies, and by Density Functional Theory (DFT) calculations. The structural properties, relative energies of different spin states and isomers, and the first calculated vertical detachment energies (VDEs) of different spin states for these (FeS)_mH^- (m = 2-4) cluster anions are investigated at various reasonable theory levels. Two types of structural isomers are found for these (FeS)_mH^- (m = 2-4) clusters: (1) the single hydrogen atom bonds to a sulfur site (SH-type); and (2) the single hydrogen atom bonds to an iron site (FeH-type). Experimental and theoretical results suggest such available different SH- and FeH-type structural isomers should be considered when evaluating the properties and behavior of these single hydrogen containing iron sulfide clusters in real chemical and biological systems. Compared to their related, respective pure iron sulfur (FeS)_m^- clusters, the first VDE trend of the diverse type (FeS)_mH_0,1^- (m = 1-4) clusters can be understood through (1) the different electron distribution properties of their highest singly occupied molecular orbital employing natural bond orbital analysis (NBO/HSOMO), and (2) the partial charge distribution on the NBO/HSOMO localized sites of each cluster anion. Generally, the properties of the NBO/HSOMOs play the principal role with regard to the physical and chemical properties of all the anions. The change of cluster VDE from low to high is associated with the change in nature of their NBO/HSOMO from a dipole bound and valence electron mixed character, to a valence p orbital on S, to a valence d orbital on Fe, and to a valence p orbital on Fe or an Fe-Fe delocalized valence bonding orbital. For clusters having the same properties for NBO/HSOMOs, the partial charge distributions at the NBO/HSOMO localized sites additionally affect their VDEs: a more negative or less positive localized charge distribution is correlated with a lower first VDE. The single hydrogen in these (FeS)_mH^- (m = 2-4) cluster anions is suggested to affect their first VDEs through the different structure types (SH- or FeH-), the nature of the NBO/HSOMOs at the local site, and the value of partial charge number at the local site of the NBO/HSOMO.

Mechanistic insights into the reductive dehydroxylation pathway for the biosynthesis of isoprenoids promoted by the IspH enzyme

Here, we report an integrated quantum mechanics/molecular mechanics (QM/MM) study of the bio-organometallic reaction pathway of the 2H+/2e- reduction of (E)-4-hydroxy-3-methylbut-2-enyl pyrophosphate (HMBPP) into the so called universal terpenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), promoted by the IspH enzyme. Our results support the viability of the bio-organometallic pathway through rotation of the OH group of HMBPP away from the [Fe4S4] cluster at the core of the catalytic site, to become engaged in a H-bond with Glu126. This rotation is synchronous with π-coordination of the C2[double bond, length as m-dash]C3 double bond of HMBPP to the apical Fe atom of the [Fe4S4] cluster. Dehydroxylation of HMBPP is triggered by a proton transfer from Glu126 to the OH group of HMBPP. The reaction pathway is completed by competitive proton transfer from the terminal phosphate group to the C2 or C4 atom of HMBPP.

Broken-Symmetry DFT Computations for the Reaction Pathway of IspH, an Iron-Sulfur Enzyme in Pathogenic Bacteria

The recently discovered methylerythritol phosphate (MEP) pathway provides new targets for the development of antibacterial and antimalarial drugs. In the final step of the MEP pathway, the [4Fe-4S] IspH protein catalyzes the 2e(-)/2H(+) reductive dehydroxylation of (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate (HMBPP) to afford the isoprenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). Recent experiments have attempted to elucidate the IspH catalytic mechanism to drive inhibitor development. Two competing mechanisms have recently emerged, differentiated by their proposed HMBPP binding modes upon 1e(-) reduction of the [4Fe-4S] cluster: (1) a Birch reduction mechanism, in which HMBPP remains bound to the [4Fe-4S] cluster through its terminal C4-OH group (ROH-bound) until the -OH is cleaved as water; and (2) an organometallic mechanism, in which the C4-OH group rotates away from the [4Fe-4S] cluster, allowing the HMBPP olefin group to form a metallacycle complex with the apical iron (η(2)-bound). We perform broken-symmetry density functional theory computations to assess the energies and reduction potentials associated with the ROH- and η(2)-bound states implicated by these competing mechanisms. Reduction potentials obtained for ROH-bound states are more negative (-1.4 to -1.0 V) than what is typically expected of [4Fe-4S] ferredoxin proteins. Instead, we find that η(2)-bound states are lower in energy than ROH-bound states when the [4Fe-4S] cluster is 1e(-) reduced. Furthermore, η(2)-bound states can already be generated in the oxidized state, yielding reduction potentials of ca. -700 mV when electron addition occurs after rotation of the HMBPP C4-OH group. We demonstrate that such η(2)-bound states are kinetically accessible both when the IspH [4Fe-4S] cluster is oxidized and 1e(-) reduced. The energetically preferred pathway gives 1e(-) reduction of the cluster after substrate conformational change, generating the 1e(-) reduced intermediate proposed in the organometallic mechanism.

Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition

With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.

Assessment of Quantum Mechanical Methods for Copper and Iron Complexes by Photoelectron Spectroscopy

Broken-symmetry density functional theory (BS-DFT) calculations are assessed for redox energetics [Cu(SCH₃)₂]^1–/0, [Cu(NCS)₂]^1–/0, [FeCl₄]^1–/0, and [Fe(SCH₃)₄]^1–/0 against vertical detachment energies (VDE) from valence photoelectron spectroscopy (PES), as a prelude to studies of metalloprotein analogs. The M06 and B3LYP hybrid functionals give VDE that agree with the PES VDE for the Fe complexes, but both underestimate it by ∼400 meV for the Cu complexes; other hybrid functionals give VDEs that are an increasing function of the amount of Hartree–Fock (HF) exchange and so cannot show good agreement for both Cu and Fe complexes. Range-separated (RS) functionals appear to give a better distribution of HF exchange since the negative HOMO energy is approximately equal to the VDEs but also give VDEs dependent on the amount of HF exchange, sometimes leading to ground states with incorrect electron configurations; the LRC-ωPBEh functional reduced to 10% HF exchange at short-range give somewhat better values for both, although still ∼150 meV too low for the Cu complexes and ∼50 meV too high for the Fe complexes. Overall, the results indicate that while HF exchange compensates for self-interaction error in DFT calculations of both Cu and Fe complexes, too much may lead to more sensitivity to nondynamical correlation in the spin-polarized Fe complexes.

Calculating standard reduction potentials of [4Fe-4S] proteins

The oxidation-reduction potentials of electron transfer proteins determine the driving forces for their electron transfer reactions. Although the type of redox site determines the intrinsic energy required to add or remove an electron, the electrostatic interaction energy between the redox site and its surrounding environment can greatly shift the redox potentials. Here, a method for calculating the reduction potential versus the standard hydrogen electrode, E°, of a metalloprotein using a combination of density functional theory and continuum electrostatics is presented. This work focuses on the methodology for the continuum electrostatics calculations, including various factors that may affect the accuracy. The calculations are demonstrated using crystal structures of six homologous HiPIPs, which give E° that are in excellent agreement with experimental results.

Understanding rubredoxin redox sites by density functional theory studies of analogues

Determining the redox energetics of redox site analogues of metalloproteins is essential in unraveling the various contributions to electron transfer properties of these proteins. Since studies of the [4Fe-4S] analogues show that the energies are dependent on the ligand dihedral angles, broken symmetry density functional theory (BS-DFT) with the B3LYP functional and double-ζ basis sets calculations of optimized geometries and electron detachment energies of [1Fe] rubredoxin analogues are compared to crystal structures and gas-phase photoelectron spectroscopy data, respectively, for [Fe(SCH(3))(4)](0/1-/2-), [Fe(S(2)-o-xyl)(2)](0/1-/2-), and Na(+)[Fe(S(2)-o-xyl)(2)](1-/2-) in different conformations. In particular, the study of Na(+)[Fe(S(2)-o-xyl)(2)](1-/2-) is the only direct comparison of calculated and experimental gas phase detachment energies for the 1-/2- couple found in the rubredoxins. These results show that variations in the inner sphere energetics by up to ∼0.4 eV can be caused by differences in the ligand dihedral angles in either or both redox states. Moreover, these results indicate that the protein stabilizes the conformation that favors reduction. In addition, the free energies and reorganization energies of oxidation and reduction as well as electrostatic potential charges are calculated, which can be used as estimates in continuum electrostatic calculations of electron transfer properties of [1Fe] proteins.

Fold versus sequence effects on the driving force for protein-mediated electron transfer

Electron transport chains composed of electron transfer reactions mainly between proteins provide fast efficient flow of energy in a variety of metabolic pathways. Reduction potentials are essential characteristics of the proteins because they determine the driving forces for the electron transfers. As both polar and charged groups from the backbone and side chains define the electrostatic environment, both the fold and the sequence will contribute. However, although the role of a specific sequence may be determined by experimental mutagenesis studies of reduction potentials, understanding the role of the fold by experiment is much more difficult. Here, continuum electrostatics and density functional theory calculations are used to analyze reduction potentials in [4Fe-4S] proteins. A key feature is that multiple homologous proteins in three different folds are compared: six high potential iron-sulfur proteins, four bacterial ferredoxins, and four nitrogenase iron proteins. Calculated absolute reduction potentials are shown to be in quantitative agreement with electrochemical reduction potentials. Calculations further demonstrate that the contribution of the backbone is larger than that of the side chains and is consistent for homologous proteins but differs for nonhomologous proteins, indicating that the fold is the major protein factor determining the reduction potential, whereas the specific amino acid sequence tunes the reduction potential for a given fold. Moreover, the fold contribution is determined mainly by the proximity of the redox site to the protein surface and the orientation of the dipoles of backbone near the redox site.

Nitrogenase structure and function relationships by density functional theory

Modern density functional theory has tremendous potential with matching popularity in metalloenzymology to reveal the unseen atomic and molecular details of structural data, spectroscopic measurements, and biochemical experiments by providing insights into unobservable structures and states, while also offering theoretical justifications for observed trends and differences. An often untapped potential of this theoretical approach is to bring together diverse experimental structural and reactivity information and allow for these to be critically evaluated at the same level. This is particularly applicable for the tantalizingly complex problem of the structure and molecular mechanism of biological nitrogen fixation. In this chapter we provide a review with extensive practical details of the compilation and evaluation of experimental data for an unbiased and systematic density functional theory analysis that can lead to remarkable new insights about the structure-function relationships of the iron-sulfur clusters of nitrogenase.

Probing ligand effects on the redox energies of [4Fe-4S] clusters using broken-symmetry density functional theory

A central issue in understanding redox properties of iron-sulfur proteins is determining the factors that tune the reduction potentials of the Fe-S clusters. Recently, Solomon and coworkers have shown that the Fe-S bond covalency of protein analogs measured by %L, the percent ligand character of the Fe 3d orbitals, from ligand K-edge X-ray absorption spectroscopy (XAS) correlates with the electrochemical redox potentials. Also, Wang and coworkers have measured electron detachment energies for iron-sulfur clusters without environmental perturbations by gas-phase photoelectron spectroscopy (PES). Here the correlations of the ligand character with redox energy and %L character are examined in [Fe(4)S(4)L(4)](2-) clusters with different ligands by broken symmetry density functional theory (BS-DFT) calculations using the B3LYP functional together with PES and XAS experimental results. These gas-phase studies assess ligand effects independently of environmental perturbations and thus provide essential information for computational studies of iron-sulfur proteins. The B3LYP oxidation energies agree well with PES data, and the %L character obtained from natural bond orbital analysis correlates with XAS values, although it systematically underestimates them because of basis set effects. The results show that stronger electron-donating terminal ligands increase %L(t), the percent ligand character from terminal ligands, but decrease %S(b), the percent ligand character from the bridging sulfurs. Because the oxidized orbital has significant Fe-L(t) antibonding character, the oxidation energy correlates well with %L(t). However, because the reduced orbital has varying contributions of both Fe-L(t) and Fe-S(b) antibonding character, the reduction energy does not correlate with either %L(t) or %S(b). Overall, BS-DFT calculations together with XAS and PES experiments can unravel the complex underlying factors in the redox energy and chemical bonding of the [4Fe-4S] clusters in iron-sulfur proteins.

Cleavage of [4Fe-4S]-type clusters: breaking the symmetry

The cleavage of [4Fe-4S]-type clusters is thought to be important in proteins such as Fe-S scaffold proteins and nitrogenase. However, most [4Fe-4S](2+) clusters in proteins have two antiferromagnetically coupled high-spin layers in which a minority spin is delocalized in each layer, thus forming a symmetric Fe(2.5+)-Fe(2.5+) pair, and how cleavage occurs between the irons is puzzling because of the shared electron. Previously, we proposed a novel mechanism for the fission of a [4Fe-4S] core into two [2Fe-2S] cores in which the minority spin localizes on one iron, thus breaking the symmetry and creating a transition state with two Fe(3+)-Fe(2+) pairs. Cleavage first through the weak Fe(2+)-S bonds lowers the activation energy. Here, we propose a test of this mechanism: break the symmetry of the cluster by changing the ligands to promote spin localization, which should enhance reactivity. The cleavage reactions for the homoligand [Fe(4)S(4)L(4)](2-) (L = SCH(3), Cl, H) and heteroligand [Fe(4)S(4)(SCH(3))(2)L(2)](2-) (L = Cl, H) clusters in the gas phase were examined via broken-symmetry density functional theory calculations. In the heteroligand clusters, the minority spin localized on the iron coordinated by the weaker electron-donor ligand, and the reaction energy and activation barrier of the cleavage were lowered, which is in accord with our proposed mechanism and consistent with photoelectron spectroscopy and collision-induced dissociation experiments. These studies suggest that proteins requiring facile fission of their [4Fe-4S] cluster in their biological function might have spin-localized [4Fe-4S] clusters.

Insight into environmental effects on bonding and redox properties of [4Fe-4S] clusters in proteins

The large differences in redox potentials between the HiPIPs and ferredoxins are generally attributed to hydrogen bonds and electrostatic effects from the protein and solvent. Recent ligand K-edge X-ray absorption studies by Solomon and co-workers show that the Fe-S covalencies of [4Fe-4S] clusters in the two proteins differ considerably apparently because of hydrogen bonds from water, indicating electronic effects may be important. However, combined density function theory (DFT) and photoelectron spectroscopy studies by our group and Wang and co-workers indicate that hydrogen bonds tune the potential of [4Fe-4S] clusters by mainly electrostatics. The DFT studies here rationalize both results, namely that the observed change in the Fe-S covalency is due to differences in ligand conformation between the two proteins rather than hydrogen bonds. Moreover, the ligand conformation affects the calculated potentials by approximately 100 mV and, thus, is a heretofore unconsidered means of tuning the potential.

Optimization of Spin-Unrestricted Density Functional Theory for Redox Properties of Rubredoxin Redox Site Analogues

Quantum chemical calculations of metal clusters in proteins for redox studies require both computational feasibility as well as accuracies of at least ∼50 mV for redox energies but only ∼0.05 Å for bond lengths. Thus, optimization of spin-unrestricted density functional theory (DFT) methods, especially the hybrid generalized gradient approximation functionals, for energies while maintaining good geometries is essential. Here, different DFT functionals with effective core potential (ECP) and full core basis sets for [Fe(SCH(3))(4)](2-/1-) and [Fe(SCH(3))(3)](1-/0), which are analogs of the iron-sulfur protein rubredoxin, are investigated in comparison to experiment as well as other more computationally intensive electron correlation methods. In particular, redox energies are calibrated against gas-phase photoelectron spectroscopy data so no approximations for the environment are needed. B3LYP gives the best balance of accuracy in energy and geometry compared B97gga1 and BHandH and is better for energies than Møller-Plesset perturbation theory series (MP2, MP3, MP4SDQ) and comparable to coupled cluster [CCSD, CCSD(T)] methods. Of the full core basis sets tested, the 6-31G** basis sets give good geometries, and addition of diffuse functions to only the sulfur significantly improves the energies. Moreover, a basis set with an ECP on only the iron gives a less accurate but still reasonable geometries and energies.

Development of a low-temperature photoelectron spectroscopy instrument using an electrospray ion source and a cryogenically controlled ion trap

The ability to control ion temperatures is critical for gas phase spectroscopy and has been a challenge in chemical physics. A low-temperature photoelectron spectroscopy instrument has been developed for the investigation of complex anions in the gas phase, including multiply charged anions, solvated species, and biological molecules. The new apparatus consists of an electrospray ionization source, a three dimensional (3D) Paul trap for ion accumulation and cooling, a time-of-flight mass spectrometer, and a magnetic-bottle photoelectron analyzer. A key feature of the new instrument is the capability to cool and tune ion temperatures from 10 to 350 K in the 3D Paul trap, which is attached to the cold head of a closed cycle helium refrigerator. Ion cooling is accomplished in the Paul trap via collisions with a background gas and has been demonstrated by observation of complete elimination of vibrational hot bands in photoelectron spectra of various anions ranging from small molecules to complex species. Further evidence of ion cooling is shown by the observation of H2-physisorbed anions at low temperatures. Cold anions result in better resolved photoelectron spectra due to the elimination of vibrational hot bands and yield more accurate energetic and spectroscopic information. Temperature-dependent studies are made possible for weakly bonded molecular and solvated clusters, allowing thermodynamic information to be obtained.

S K-edge XAS and DFT calculations on square-planar NiII-thiolate complexes: effects of active and passive H-bonding

S K-edge XAS for a low-spin NiII-thiolate complex shows a 0.2 eV shift to higher pre-edge energy but no change in Ni-S bond covalency upon H-bonding. This is different from the H-bonding effect we observed in high-spin FeIII-thiolate complexes where there is a significant decrease in Fe-S bond covalency but no change in energy due to H-bonding (Dey, A.; Okamura, T.-A.; Ueyama, N.; Hedman, B.; Hodgson, K. O.; Solomon, E. I. J. Am. Chem. Soc. 2005, 127, 12046-12053). These differences were analyzed using DFT calculations, and the results indicate that two different types of H-bonding interactions are possible in metal-thiolate systems. In the high-spin FeIII-thiolate case, the H-bonding involves a thiolate donor orbital which is also involved in bonding with the metal (active), while in the low-spin NiII-thiolate, the orbital involved in H-bonding is nonbonding with respect to the M-S bonding (passive). The contributions of active and passive H-bonds to the reduction potential and Lewis acid properties of a metal center are evaluated.

Ultrahigh-resolution study on Pyrococcus abyssi rubredoxin: II. Introduction of an O–H···Sγ–Fe hydrogen bond increased the reduction potential by 65 mV

The effect of D-H...S(gamma)-Fe hydrogen bonding on the reduction potential of rubredoxin was investigated by the introduction of an O-H...S(gamma)-Fe hydrogen bond on the surface of Pyrococcus abyssi rubredoxin. The formation of a weak hydrogen bond between Ser44-O(gamma) and Cys42-S(gamma) in mutant W4L/R5S/A44S increased the reduction potential by 56 mV. When side effects of the mutation were taken into account, the contribution of the additional cluster hydrogen bond to the reduction potential was estimated to be +65 mV. The structural analysis was based on ultrahigh-resolution structures of oxidized P. abyssi rubredoxin W4L/R5S and W4L/R5S/A44S refined to 0.69 and 0.86 A, respectively.

Observation of Cysteine Thiolate and-S···H−O Intermolecular Hydrogen Bond

The cysteine anion was produced in the gas phase by electrospray ionization and investigated by photoelectron spectroscopy at low temperature (70 K). The cysteine anion was found to exhibit the thiolate form [-SCH2CH(NH2)CO2H], rather than the expected carboxylate form [HSCH2CH(NH2)CO2-]. This observation was confirmed by two control experiments, that is, methyl cysteine [CH3SCH2CH(NH2)CO2-] and cysteine methyl ester [-SCH2CH(NH2)CO2CH3]. The electron binding energy of [-SCH2CH(NH2)CO2H] was measured to be about 0.7 eV blue-shifted relative to [-SCH2CH(NH2)CO2CH3] due to the formation of an intramolecular -S-...HO2C- hydrogen bond in the cysteine thiolate. Theoretical calculations at the CCSD(T)/6-311++G(2df,p) and B3LYP/6-311++G(2df,p) levels were carried out to estimate the strength of this intramolecular -S-...HO2C- hydrogen bond. Combining experimental measurements and theoretical calculations yielded an estimated value of 16.4 +/- 2.0 kcal/mol for the -S-...HO2C- intramolecular hydrogen-bond strength.

Direct influence of hydrogen-bonding on the reduction potential of a CuII center

Two hydrogen-bonds from geometrically constrained OH groups to coordinated oxygen donors shift the reduction potential of a Cu(II) complex by +270 mV as compared to the structurally analogous reference complex missing the OH groups.

Sulfur K-edge XAS and DFT calculations on P450 model complexes: effects of hydrogen bonding on electronic structure and redox potentials

Hydrogen bonding (H-bonding) is generally thought to play an important role in tuning the electronic structure and reactivity of metal-sulfur sites in proteins. To develop a quantitative understanding of this effect, S K-edge X-ray absorption spectroscopy (XAS) has been employed to directly probe ligand-metal bond covalency, where it has been found that protein active sites are significantly less covalent than their related model complexes. Sulfur K-edge XAS data are reported here on a series of P450 model complexes with increasing H-bonding to the ligated thiolate from its substituent. The XAS spectroscopic results show a dramatic decrease in preedge intensity. DFT calculations reproduce these effects and show that the observed changes are in fact solely due to H-bonding and not from the inductive effect of the substituent on the thiolate. These calculations also indicate that the H-bonding interaction in these systems is mainly dipolar in nature. The -2.5 kcal/mol energy of the H-bonding interaction was small relative to the large change in ligand-metal bond covalency (30%) observed in the data. A bond decomposition analysis of the total energy is developed to correlate the preedge intensity change to the change in Fe-S bonding interaction on H-bonding. This effect is greater for the reduced than the oxidized state, leading to a 260 mV increase in the redox potential. A simple model shows that E degrees should vary approximately linearly with the covalency of the Fe-S bond in the oxidized state, which can be determined directly from S K-edge XAS.

Probing the Electronic Structure of [2Fe-2S] Clusters with Three Coordinate Iron Sites by Use of Photoelectron Spectroscopy

Five series of [2Fe-2S] complexes, [Fe(2)S(2)Cl(2)(-)(x)(CN)(x)](-), [Fe(2)S(2)(SEt)(2)(-)(x)Cl(x)](-), [Fe(2)S(2)(SEt)(2)(-)(x)(CN)(x)](-), [Fe(2)S(2)Cl(2)(-)(x)(OAc)(x)](-) (OAc = acetate), and [Fe(2)S(2)(SEt)(2)(-)(x)(OPr)(x)](-) (OPr = propionate) (x = 0-2), were produced by collision-induced dissociation of the corresponding [4Fe-4S] complexes, and their electronic structures were studied by photoelectron spectroscopy. All the [2Fe-2S] complexes contain a [Fe(2)S(2)](+) core similar to that in reduced [2Fe] ferredoxins but with different coordination geometries. For the first three series, which only involve tricoordinated Fe sites, a linear relationship between the measured binding energies and the substitution number (x) was observed, revealing the independent ligand contributions to the total electron binding energies. The effect of the ligand increases in the order SEt --> Cl --> CN, conforming to their electron-withdrawing ability in the same order. The carboxylate ligands in the [Fe(2)S(2)Cl(2)(-)(x)(OAc)(x)](-) and [Fe(2)S(2)(SEt)(2)(-)(x)(OPr)(x)](-) complexes were observed to act as bidentate ligands, giving rise to tetracoordinated iron sites. This is different from their monodentate coordination behavior in the [4Fe-4S] cubane complexes, reflecting the high reactivity of the unsatisfied three-coordinate iron site in the [2Fe-2S] complexes. The [2Fe-2S] complexes with tetracoordinated iron sites exhibit lower electron binding energies, that is, higher reductive activity than the all tricoordinate planar clusters. The electronic structures of all the [2Fe-2S] complexes were shown to conform to the "inverted energy level scheme".