Structure of [4Fe-4S] ferredoxin from Bacillus thermoproteolyticus refined at 2.3 A resolution. Structural comparisons of bacterial ferredoxins

Two local minima for structures of [4Fe-4S] clusters obtained with density functional theory methods

[4Fe-4S] clusters are essential cofactors in many proteins involved in biological redox-active processes. Density functional theory (DFT) methods are widely used to study these clusters. Previous investigations have indicated that there exist two local minima for these clusters in proteins. We perform a detailed study of these minima in five proteins and two oxidation states, using combined quantum mechanical and molecular mechanical (QM/MM) methods. We show that one local minimum (L state) has longer Fe-Fe distances than the other (S state), and that the L state is more stable for all cases studied. We also show that some DFT methods may only obtain the L state, while others may obtain both states. Our work provides new insights into the structural diversity and stability of [4Fe-4S] clusters in proteins, and highlights the importance of reliable DFT methods and geometry optimization. We recommend r2SCAN for optimizing [4Fe-4S] clusters in proteins, which gives the most accurate structures for the five proteins studied.

Structural and functional characterisation of HepTH1-5 peptide as a potential hepcidin replacement

Hepcidin is a principal regulator of iron homeostasis and its dysregulation has been recognised as a causative factor in cancers and iron disorders. The strategy of manipulating the presence of hepcidin peptide has been used for cancer treatment. However, this has demonstrated poor efficiency and has been short-lived in patients. Many studies reported using minihepcidin therapy as an alternative way to treat hepcidin dysregulation, but this was only applied to non-cancer patients. Highly conserved fish hepcidin protein, HepTH1-5, was investigated to determine its potential use in developing a hepcidin replacement for human hepcidin (Hepc25) and as a therapeutic agent by targeting the tumour suppressor protein, p53, through structure-function analysis. The authors found that HepTH1-5 is stably bound to ferroportin, compared to Hepc25, by triggering the ferroportin internalisation via Lys42 and Lys270 ubiquitination, in a similar manner to the Hepc25 activity. Moreover, the residues Ile24 and Gly24, along with copper and zinc ligands, interacted with similar residues, Lys24 and Asp1 of Hepc25, respectively, showing that those molecules are crucial to the hepcidin replacement strategy. HepTH1-5 interacts with p53 and activates its function through phosphorylation. This finding shows that HepTH1-5 might be involved in the apoptosis signalling pathway upon a DNA damage response. This study will be very helpful for understanding the mechanism of the hepcidin replacement and providing insights into the HepTH1-5 peptide as a new target for hepcidin and cancer therapeutics.Communicated by Ramaswamy H. Sarma.

Benchmark Study of Redox Potential Calculations for Iron-Sulfur Clusters in Proteins

Redox potentials have been calculated for 12 different iron–sulfur sites of 6 different types with 1–4 iron ions. Structures were optimized with combined quantum mechanical and molecular mechanical (QM/MM) methods, and the redox potentials were calculated using the QM/MM energies, single-point QM methods in a continuum solvent or by QM/MM thermodynamic cycle perturbations. We show that the best results are obtained with a large QM system (∼300 atoms, but a smaller QM system, ∼150 atoms, can be used for the QM/MM geometry optimization) and a large value of the dielectric constant (80). For absolute redox potentials, the B3LYP density functional method gives better results than TPSS, and the results are improved with a larger basis set. However, for relative redox potentials, the opposite is true. The results are insensitive to the force field (charges of the surroundings) used for the QM/MM calculations or whether the protein and solvent outside the QM system are relaxed or kept fixed at the crystal structure. With the best approach for relative potentials, mean absolute and maximum deviations of 0.17 and 0.44 V, respectively, are obtained after removing a systematic error of −0.55 V. Such an approach can be used to identify the correct oxidation states involved in a certain redox reaction.

We have studied redox potentials of 12 iron−sulfur sites of 6 types with 1−4 iron ions. Structures were optimized with combined quantum mechanical and molecular mechanical (QM/MM) methods, and the redox potentials were calculated with QM/MM, QM calculations in a continuum solvent or by QM/MM thermodynamic cycle perturbations. The best results are obtained with the second approach using ∼300 atoms in the QM model and a large dielectric constant.

Genetic variation in metronidazole metabolism and oxidative stress pathways in clinical Giardia lamblia assemblage A and B isolates

Purpose: Treatment-refractory Giardia cases have increased rapidly within the last decade. No markers of resistance nor a standardized susceptibility test have been established yet, but several enzymes and their pathways have been associated with metronidazole (MTZ) resistant Giardia. Very limited data are available regarding genetic variation in these pathways. We aimed to investigate genetic variation in metabolic pathway genes proposed to be involved in MTZ resistance in recently acquired, cultured clinical isolates.

Methods: Whole genome sequencing of 12 assemblage A2 and 8 assemblage B isolates was done, to decipher genomic variation in Giardia. Twenty-nine genes were identified in a literature search and investigated for their single nucleotide variants (SNVs) in the coding/non-coding regions of the genes, either as amino acid changing (non-synonymous SNVs) or non-changing SNVs (synonymous).

Results: In Giardia assemblage B, several genes involved in MTZ activation or oxidative stress management were found to have higher numbers of non-synonymous SNVs (thioredoxin peroxidase, nitroreductase 1, ferredoxin 2, NADH oxidase, nitroreductase 2, alcohol dehydrogenase, ferredoxin 4 and ferredoxin 1) than the average variation. For Giardia assemblage A2, the highest genetic variability was found in the ferredoxin 2, ferredoxin 6 and in nicotinamide adenine dinucleotide phosphate (NADPH) oxidoreductase putative genes. SNVs found in the ferredoxins and nitroreductases were analyzed further by alignment and homology modeling. SNVs close to the iron-sulfur cluster binding sites in nitroreductase-1 and 2 and ferredoxin 2 and 4 could potentially affect protein function. Flavohemoprotein seems to be a variable-copy gene, due to higher, but variable coverage compared to other genes investigated.

Conclusion: In clinical Giardia isolates, genetic variability is common in important genes in the MTZ metabolizing pathway and in the management of oxidative and nitrosative stress and includes high numbers of non-synonymous SNVs. Some of the identified amino acid changes could potentially affect the respective proteins important in the MTZ metabolism.

Role of the C-terminal extension stacked on the re-face of the isoalloxazine ring moiety of the flavin adenine dinucleotide prosthetic group in ferredoxin-NADP(+) oxidoreductase from Bacillus subtilis

Ferredoxin-NADP(+) oxidoreductase [EC 1.18.1.2] from Bacillus subtilis (BsFNR) is homologous to the bacterial NADPH-thioredoxin reductase, but possesses a unique C-terminal extension that covers the re-face of the isoalloxazine ring moiety of the flavin adenine dinucleotide (FAD) prosthetic group. In this report, we utilize BsFNR mutants depleted of their C-terminal residues to examine the importance of the C-terminal extension in reactions with NADPH and ferredoxin (Fd) from B. subtilis by spectroscopic and steady-state reaction analyses. The depletions of residues Y313 to K332 (whole C-terminal extension region) and S325 to K332 (His324 intact) resulted in significant increases in the catalytic efficiency with NADPH in diaphorase assay with ferricyanide, whereas Km values for ferricyanide were increased. In the cytochrome c reduction assay in the presence of B. subtilis ferredoxin, the S325-K332 depleted mutant displayed a significant decrease in the turnover rate with an Fd concentration range of 1-10 μM. The Y313-K332 depleted mutant demonstrated an increase in the rate of the direct reduction of horse heart cytochrome c in the absence of Fd. These data indicated that depletion of the C-terminal extension plays an important role in the reaction of BsFNR with ferredoxin.

The structure of a novel electron-transfer ferredoxin from Rhodopseudomonas palustris HaA2 which contains a histidine residue in its iron-sulfur cluster-binding motif

Rhodopseudomonas palustris HaA2 contains a gene, RPB3630, encoding a ferredoxin, HaPuxC, with an atypical CXXHXXC(X)nCP iron-sulfur cluster-binding motif. The ferredoxin gene is associated with a cytochrome P450 (CYP) monooxygenase-encoding gene, CYP194A3, an arrangement which is conserved in several strains of bacteria. Similar ferredoxin genes are found in other bacteria, such as Mycobacterium tuberculosis, where they are also associated with CYP genes. The crystal structure of HaPuxC has been solved at 2.3 Å resolution. The overall fold of this [3Fe-4S] cluster-containing ferredoxin is similar to other [3Fe-4S] and [4Fe-4S] species, with the loop around the iron-sulfur cluster more closely resembling those of [3Fe-4S] ferredoxins. The side chain of His17 from the cluster-binding motif in HaPuxC points away from the vacant site of the cluster and interacts with Glu61 and one of the sulfide ions of the cluster. This is the first cytochrome P450 electron-transfer partner of this type to be structurally characterized and will provide a better understanding of the electron-transfer processes between these ferredoxins and their CYP enzymes.

Heterometallic [AgFe(3)S (4)] ferredoxin variants: synthesis, characterization, and the first crystal structure of an engineered heterometallic iron-sulfur protein

Heterometallic [AgFe(3)S(4)] iron-sulfur clusters assembled in wild-type Pyrococcus furiosus ferredoxin and two variants, D14C and D14H, are characterized. The crystal structure of the [AgFe(3)S(4)] D14C variant shows that the silver(I) ion is indeed part of the cluster and is coordinated to the thiolate group of residue 14. Cyclic voltammetry shows one redox pair with a reduction potential of +220 mV versus the standard hydrogen electrode which is assigned to the [AgFe(3)S(4)](2+/+) couple. The oxidized form of the [AgFe(3)S(4)] D14C variant is stable in the presence of dioxygen, whereas the oxidized forms of the [AgFe(3)S(4)] wild type and D14H variants convert to the [Fe(3)S(4)] ferredoxin form. The monovalent d (10) silver(I) ion stabilizes the [Fe(3)S(4)](+/0) cluster fragment, as opposed to divalent d (10) metal ions, resulting in more than 0.4 V difference in reduction potentials between the silver(I) and, e.g., zinc(II) heterometallic [MFe(3)S(4)] ferredoxins. The trend in reduction potentials for the variants containing the [AgFe(3)S(4)] cluster is wild type ≤ D14C < D14H and shows the same trend as reported for the variants containing the [Fe(3)S(4)] cluster, but is different from the D14C < D14H < wild type trend reported for the [Fe(4)S(4)] ferredoxin. The similarity in the reduction potential trend for the variants containing the heterometallic [AgFe(3)S(4)] cluster and the [Fe(3)S(4)] cluster can be rationalized in terms of the electrostatic influence of the residue 14 side chains, rather than the dissociation constant of this residue, as is the case for [Fe(4)S(4)] ferredoxins. The trends in reduction potentials are in line with there being no electronic coupling between the silver(I) ion and the Fe(3)S(4) fragment.

HdrC2 from Acidithiobacillus ferrooxidans owns two iron-sulfur binding motifs but binds only one variable cluster between [4Fe-4S] and [3Fe-4S]

The heterodisulfide reductase complex HdrABC from Acidithiobacillus ferrooxidans was suggested to own novel features that act in reverse to convert the sulfane sulfur of GS( n )H species (n > 1) into sulfite in sulfur oxidation. The HdrC subunit is potentially encoded by two different highly upregulated genes sharing only 29 % identity in A. ferrooxidans grown in sulfur-containing medium, which were named as HdrC1 and HdrC2, respectively and had been confirmed to contain iron-sulfur cluster by expression and characterization, especially the HdrC1 which had been showed to bind only one [4Fe-4S] cluster by mutations. However, the mutations of the HdrC2 remain to be done and the detailed binding information of it is still unclear. Here, we report the expression, mutations, and molecular modeling of the HdrC2 from A. ferrooxidans. This HdrC2 had two identical motifs (Cx(2)Cx(2)Cx(3)C) containing total of eight cysteine residues potentially for iron-sulfur cluster binding. This purified HdrC2 was exhibited to contain one variable cluster converted between [4Fe-4S] and [3Fe-4S] according to different conditions by the UV-scanning and EPR spectra. The site-directed mutagenesis results of these eight residues further confirmed that the HdrC2 in reduction with Fe(2+) condition loaded only one [4Fe-4S](+) with spin S = 1/2 ligated by the residues of Cys73, Cys109, Cys112, and Cys115; the HdrC2 in natural aeration condition lost the Fe atom ligated by the residue of Cys73 and loaded only one [3Fe-4S](0) with spin S = 0; the HdrC2 in oxidation condition loaded only one [3Fe-4S](+) with spin S = 1/2. Molecular modeling results were also in line with the experiment results.

Expression, purification and molecular modeling of another HdrC from Acidithiobacillus ferrooxidans which binds only one [4Fe-4S] cluster

The heterodisulfide reductase complex HdrABC from Acidithiobacillus ferrooxidans was predicted to have novel features that work in reverse to catalyse the sulfane sulfur of GSnH species (n > 1) into sulfite in sulfur oxidation. There are two different highly upregulated genes potentially encoding the HdrC subunit in A. ferrooxidans grown in sulfur-containing medium. An HdrC containing iron-sulfur cluster from A. ferrooxidans corresponding to one of the genes had been expressed and biophysically characterized. Comparatively, here we report the cloning, expression, and characterization of another HdrC from A. ferrooxidans. This HdrC was expressed in inclusion bodies in all conditions tested. This purified HdrC displayed brown color and contained the [4Fe-4S] cluster confirmed by the UV-scanning and EPR spectra. This HdrC owned two identical motifs (Cx(2)Cx(2)Cx(3)C) including total of eight cysteine residues for [4Fe-4S] cluster binding. To our surprise, the site-directed mutagenesis results of these eight residues revealed that respective removal of the sulfhydryl group of Cys73, Cys76, Cys79, and Cys37 resulted in the cluster loss, but those of Cys27, Cys30, Cys33, and Cys83 had no influence, which demonstrated that this HdrC bound only one cluster, and it might be responsible for causing the HdrABC in A. ferrooxidans working in reverse. Molecular modeling results also supported the above results and showed that this cluster was ligated by Cys73, Cys76, and Cys79 in one motif and Cys37, however, in another motif.

Chloroplast DnaJ-like proteins 3 and 4 (CDJ3/4) from Chlamydomonas reinhardtii contain redox-active Fe–S clusters and interact with stromal HSP70B

In the present study we report on the identification and characterization of three novel chloroplast-targeted DnaJ-like proteins CDJ3–5, which in addition to their J-domains contain bacterial-type ferredoxin domains. In sequence databases we could identify homologues of CDJ3–5 in green algae, moss and higher plants, but not in cyanobacteria. Phylogenetic analyses allowed us to distinguish two clades containing CDJ3/4 and CDJ5 that must have diverged early in the ancestor of the ‘green lineage’ and have further diversified later on. Molecular and biochemical analysis of CDJ3 and CDJ4 from Chlamydomonas reinhardtii revealed that both proteins are weakly expressed and appear to be localized to the stroma and to thylakoid membranes respectively. The low transcript levels of the CDJ3 and CDJ4 genes declined even further in the initial phase of heat shock, but CDJ3 transcript levels strongly increased after a dark-to-light shift. Accordingly, the Arabidopsis orthologue of CDJ5 was also found to be light-inducible and to be under strong circadian control. CDJ3 and CDJ4 proteins could both be expressed in Escherichia coli and had redox-active Fe–S clusters. In vitro cross-linking studies demonstrated that CDJ3 and CDJ4 interact with chloroplast ATP-bound HSP70B (heat-shock protein 70B), presumably as dimers, and immunoprecipitation studies showed that CDJ3/4 were also in a complex with HSP70B in Chlamydomonas cell extracts. Finally, CDJ3 was found in complexes with apparent molecular masses of approx. 550–2800 kDa, which appeared to contain RNA. We speculate that the CDJ3–5 proteins might represent redox switches that act by recruiting HSP70B for the reorganization of regulatory protein complexes.

The Mycobacterium tuberculosis cytochrome P450 system

Tuberculosis remains a leading cause of human mortality. The emergence of strains of Mycobacterium tuberculosis, the causative agent, that are resistant to the major frontline antitubercular drugs increases the urgency for the development of new therapeutic agents. Sequencing of the M. tuberculosis genome revealed the existence of 20 cytochrome P450 enzymes, some of which are potential candidates for drug targeting. The recent burst of studies reporting microarray-based gene essentiality and transcriptome analyses under in vitro, ex vivo and in vivo conditions highlight the importance of selected P450 isoforms for M. tuberculosis viability and pathogenicity. Current knowledge of the structural and biochemical properties of the M. tuberculosis P450 enzymes and their putative redox partners is reviewed, with an emphasis on findings related to their physiological function(s) as well as their potential as drug targets.

The crystal structure of FdxA, a 7Fe ferredoxin from Mycobacterium smegmatis

Mycobacterium smegmatis ferredoxin FdxA, which has an orthologue ferredoxin in Mycobacterium tuberculosis, FdxC, contains both one [3Fe-4S] and one [4Fe-4S] cluster. M. smegmatis FdxA has been shown to be a preferred ferredoxin substrate of FprA [F. Fischer, D. Raimondi, A. Aliverti, G. Zanetti, Mycobacterium tuberculosis FprA, a novel bacterial NADPH-ferredoxin reductase, Eur. J. Biochem. 269 (2002) 3005-3013], an adrenodoxin reductase-like flavoprotein of M. tuberculosis, suggesting that M. tuberculosis FdxC could be the physiological partner of the enzyme in providing reducing power to the cytochromes P450. We report here the crystal structure of FdxA at 1.6A resolution (R(factor) 16.5%, R(free) 20.2%). Besides providing an insight on protein architecture for this 106-residue ferredoxin, our crystallographic investigation highlights lability of the [4Fe-4S] center, which is shown to loose a Fe atom during crystal growth. Due to their high similarity (87% sequence identity), the structure here reported can be considered a valuable model for M. tuberculosis FdxC, thus representing a step forward in the study of the complex mycobacterial redox pathways.

Multiple orientations in a physiological complex: the pyruvate-ferredoxin oxidoreductase-ferredoxin system

Ferredoxin I from Desulfovibrio africanus (Da FdI) is a small acidic [4Fe-4S] cluster protein that exchanges electrons with pyruvate-ferredoxin oxidoreductase (PFOR), a key enzyme in the energy metabolism of anaerobes. The thermodynamic properties and the electron transfer between PFOR and either native or mutated FdI have been investigated by microcalorimetry and steady-state kinetics, respectively. The association constant of the PFOR-FdI complex is 3.85 x 10(5) M(-1), and the binding affinity has been found to be highly sensitive to ionic strength, suggesting the involvement of electrostatic forces in formation of the complex. Surprisingly, the punctual or combined neutralizations of carboxylate residues surrounding the [4Fe-4S] cluster slightly affect the PFOR-FdI interaction. Furthermore, hydrophobic residues around the cluster do not seem to be crucial for the PFOR-FdI system activity; however, some of them play an important role in the stability of the FeS cluster. NMR restrained docking associated with site-directed mutagenesis studies suggested the presence of various interacting sites on Da FdI. The modification of additional acidic residues at the interacting interface, generating a FdI pentamutant, evidenced at least two distinct FdI binding sites facing the distal [4Fe-4S] cluster of the PFOR. We also used a set of various small acidic partners to investigate the specificity of PFOR toward redox partners. The remarkable flexibility of the PFOR-FdI system supports the idea that the specificity of the physiological complex has probably been "sacrificed" to improve the turnover rate and thus the efficiency of bacterial electron transfer.

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

To probe how H-bonding effects the redox potential changes in Fe-S proteins, we produced and studied a series of gaseous cubane-type analogue complexes, [Fe(4)S(4)(SEt)(3)(SC(n)H(2n+1))](2-) and [Fe(4)S(4)(SEt)(3)(SC(n)H(2n)OH)](2-) (n = 4, 6, 11; Et = C(2)H(5)). Intrinsic redox potentials for the [Fe(4)S(4)](2+/3+) redox couple involved in these complexes were measured by photoelectron spectroscopy. The oxidation energies from [Fe(4)S(4)(SEt)(3)(SC(n)H(2n)OH)](2-) to [Fe(4)S(4)(SEt)(3)(SC(n)H(2n)OH)](-) were determined directly from the photoelectron spectra to be approximately 130 meV higher than those for the corresponding [Fe(4)S(4)(SEt)(3)(SC(n)H(2n+1))](2-) systems, because of the OH...S hydrogen bond in the former. Preliminary Monte Carlo and density functional calculations showed that the H-bonding takes place between the -OH group and the S on the terminal ligand in [Fe(4)S(4)(SEt)(3)(SC(6)H(12)OH)](2-). The current data provide a direct experimental measure of a net H-bonding effect on the redox potential of [Fe(4)S(4)] clusters without the perturbation of other environmental effects.

Expression, purification and characterisation of a Bacillus subtilis ferredoxin: a potential electron transfer donor to cytochrome P450 BioI

The fer gene from Bacillus subtilis has been subcloned and overexpressed in Escherichia coli and the protein (Fer) purified to homogeneity. N-Terminal sequencing and mass spectrometry indicate that the initiator methionine is removed from the protein and that the molecular mass is 8732 Da consistent with that deduced from the gene sequence. Amino-acid sequence comparisons indicate that Fer is a ferredoxin containing a 4Fe-4S cluster. The electron paramagnetic resonance spectrum of the reduced form of Fer is typical for a [4Fe-4S](+) cluster showing rhombic signals with g values of 2.07, 1.93 and 1.88. Reduced Fer also gives rise to a magnetic circular dichroism spectrum typical of a [4Fe-4S](+) cluster. Potentiometric titrations indicate that Fer has a reduction potential of -385+/-10 mV for the [4Fe-4S](+)-[4Fe-4S](2+) redox couple, well within the normal range expected for such a ferredoxin. A proposed physiological role for Fer is as an electron donor to cytochrome P450 BioI. Studies on Fer binding to P450 BioI give rise to a K(d) value of 0.87+/-0.10 microM. Anaerobic experiments using CO-saturated buffer indicate that Fer is indeed capable of transferring electrons to this cytochrome P450 albeit at a fairly low rate.

Solution structure of the pro-hormone convertase 1 pro-domain from Mus musculus

The solution structure of the mouse pro-hormone convertase (PC) 1 pro-domain was determined using heteronuclear NMR spectroscopy and is the first structure to be obtained for any of the domains in the convertase family. The ensemble of NMR-derived structures shows a well-ordered core consisting of a four-stranded antiparallel beta-sheet with two alpha-helices packed against one side of this sheet. Sequence homology suggests that the other eukaryotic PC pro-domains will have the same overall fold and most of the residues forming the hydrophobic core of PC1 are highly conserved within the PC family. However, some of the core residues are predicted by homology to be replaced by polar amino acid residues in other PC pro-domains and this may help to explain their marginal stability. Interestingly, the folding topology observed here is also seen for the pro-domain of bacterial subtilisin despite little or no sequence homology. Both the prokaryotic and eukaryotic structures have hydrophobic residues clustered on the solvent-accessible surface of their beta-sheets although the individual residue types differ. In the bacterial case this region is buried at the binding interface with the catalytic domain and, in the eukaryotic PC family, these surface residues are conserved. We therefore propose that the hydrophobic patch in the PC1 pro-domain is involved in the binding interface with its cognate catalytic domain in a similar manner to that seen for the bacterial system. The PC1 pro-domain structure also reveals potential mechanisms for the acid-induced dissociation of the complex between pro- and catalytic domains.

Atomic resolution structures of oxidized [4Fe-4S] ferredoxin from Bacillus thermoproteolyticus in two crystal forms: systematic distortion of [4Fe-4S] cluster in the protein

Diffraction data of two crystal forms (forms I and II) of [4Fe-4S] ferredoxin from Bacillus thermoproteolyticus have been collected to 0.92 A and 1.00 A resolutions, respectively, at 100 K using synchrotron radiation. Anisotropic temperature factors were introduced for all non-hydrogen atoms in the refinement with SHELX-97, in which stereochemical restraints were applied to the protein chain but not to the [4Fe-4S] cluster. The final crystallographic R-factors are 9.8 % for 7.0-0.92 A resolution data of the form I and 11.2 % for the 13.3-1.0 A resolution data of the form II. Many hydrogen atoms as well as multiple conformations for several side-chains have been identified. The present refinement has revised the conformations of several peptide bonds and side-chains assigned previously at 2.3 A resolution; the largest correction was that the main-chain of Pro1 and the side-chain of Lys2 were changed by rotating the C(alpha)-C bond of Lys2. Although the overall structures in the two crystal forms are very similar, conformational differences are observed in the two residues at the middle (Glu29 and Asp30) and the C-terminal residues, which have large temperature factors. The [4Fe-4S] cluster is a distorted cube with non-planar rhombic faces. Slight but significant compression of the four Fe-S bonds along one direction is observed in both crystal forms, and results in the D(2d) symmetry of the cluster. The compressed direction of the cluster relative to the protein is conserved in the two crystal forms and consistent with that in one of the clusters in Clostridium acidurici ferredoxin.

Solution NMR characterization of the thermodynamics of the disulfide bond orientational isomerism and its effect of cluster electronic properties for the hyperthermostable three-iron cluster ferredoxin from the archaeon Pyrococcus furiosus

The thermodynamics and dynamics of the Cys21-Cys48 disulfide "S" if "R" conformational isomerism in the three-iron, single cubane cluster ferredoxin (Fd) from the hyperthermophilic archaeon Pyrococcus furiosus (Pf) have been characterized by (1)H NMR spectroscopy in both water and water/methanol mixed solvents. The mean interconversion rate at 25 degrees C is 3 x 10(3) s(-1) and DeltaG(298) = -0.2 kcal/mol [DeltaH = 4.0 kcal/mol; DeltaS = 14 cal/(mol.K)], with the S orientation as the more stable form at low temperature (< 0 degrees C) but the R orientation predominating at >100 degrees C, where the organism thrives. The distinct pattern of ligated Cys beta-proton contact shifts for the resolved signals and their characteristic temperature behavior for the forms of the 3Fe Fd with alternate disulfide orientations have been analyzed to determine the influences of disulfide orientation and methanol cosolvent on the topology of the inter-iron spin coupling in the 3Fe cluster. The Cys21-Cys48 disulfide orientation influences primarily the spin couplings involving the iron ligated to Cys17, whose carbonyl oxygen is a hydrogen bond acceptor to the Cys21 peptide proton. Comparison of the Cys beta-proton contact shift pattern for the alternate disulfide orientations with the pattern exhibited upon cleaving the disulfide bridge confirms an earlier [Wang, P.-L., Calzolai, L., Bren, K. L., Teng, Q., Jenney, F. E., Jr., Brereton, P. S., Howard, J. B., Adams, M. W. W., and La Mar, G. N. (1999) Biochemistry 38, 8167-8178] proposal that the structure of the same Fd with the R disulfide orientation resembles that of the Fd upon cleaving the disulfide bond.

Sequence determination of reduction potentials by cysteinyl hydrogen bonds and peptide pipoles in [4Fe-4S] ferredoxins

A sequence determinant of reduction potentials is reported for bacterial [4Fe-4S]-type ferredoxins. The residue that is four residues C-terminal to the fourth ligand of either cluster is generally an alanine or a cysteine. In five experimental ferredoxin structures, the cysteine has the same structural orientation relative to the nearest cluster, which is stabilized by the SH...S bond. Although such bonds are generally considered weak, indications that Fe-S redox site sulfurs are better hydrogen-bond acceptors than most sulfurs include the numerous amide NH...S bonds noted by Adman and our quantum mechanical calculations. Furthermore, electrostatic potential calculations of 11 experimental ferredoxin structures indicate that the extra cysteine decreases the reduction potential relative to an alanine by approximately 60 mV, in agreement with experimental mutational studies. Moreover, the decrease in potential is due to a shift in the polar backbone stabilized by the SH...S bond rather than to the slightly polar cysteinyl side chain. Thus, these cysteines can "tune" the reduction potential, which could optimize electron flow in an electron transport chain. More generally, hydrogen bonds involving sulfur can be important in protein structure/function, and mutations causing polar backbone shifts can alter electrostatics and thus affect redox properties or even enzymatic activity of a protein.

Structural and kinetic studies of the pyruvate-ferredoxin oxidoreductase/ferredoxin complex from Desulfovibrio africanus

The pyruvate-ferredoxin oxidoreductase (PFOR)/ferredoxin (Fd) system of Desulfovibrio africanus has been investigated with the aim of understanding more fully protein-protein interaction and the kinetic characteristics of electron transfer between the two redox partners. D. africanus contains three Fds (Fd I, Fd II and Fd III) able to function as electron acceptors for PFOR. The complete amino acid sequence of Fd II was determined by automatic Edman degradation. It revealed a striking similarity to that of Fd I. The protein consists of 64 residues and its amino acid sequence is in agreement with a molecular mass of 6822.5 Da as measured by electrospray MS. Fd II contains five cysteine residues of which the first four (Cys11, Cys14, Cys17 and Cys54) are likely ligands for the single [4Fe-4S] cluster. A covalently cross-linked complex between PFOR and Fd I or Fd II was obtained by using a water soluble carbodiimide. This complex exhibited a stoichiometry of one ferredoxin for one PFOR subunit and is dependent on the ionic strength. The second-order rate constants for electron transfer between PFOR and Fds determined electrochemically using cyclic voltammetry are 7 x 107 M-1.s-1 for Fd I and 2 x 107 M-1.s-1 for Fd II and Fd III. The Km values of PFOR for Fd I and Fd II measured both by the electrochemical and the spectrophotometric method have been found to be 3 microM and 5 microM, respectively. The three-dimensional modelling of Fd II and surface analysis of Fd I, Fd II and PFOR suggest that a protein-protein complex is likely to be formed between aspartic acid/glutamic acid invariant residues of Fds and lysine residues surrounding the distal [4Fe-4S] cluster of PFOR. All of these studies are indicative of the involvement of electrostatic interactions between the two redox partners.

Secondary structure extensions in Pyrococcus furiosus ferredoxin destabilize the disulfide bond relative to that in other hyperthermostable ferredoxins. Global consequences for the disulfide orientational heterogeneity

The single cubane cluster ferredoxin (Fd) from the hyperthermophilic archaeon Pyrococcus furiosus (Pf) possesses several unique properties when compared even to Fds from other hyperthermophilic archaea or bacteria. These include an equilibrium molecular heterogeneity, a six- to seven-residue increase in size, an Asp rather than the Cys as one cluster ligand, and a readily reducible disulfide bond. NMR assignments and determination of both secondary structure and tertiary contacts remote from the paramagnetic oxidized cluster of Pf 3Fe Fd with an intact disulfide bond reported previously (Teng Q., Zhou, Z. H., Smith, E. T., Busse, S. C., Howard, J. B. Adams, M. W. W., and La Mar, G. (1994) Biochemistry 33, 6316-6328) are extended here to the 4Fe oxidized cluster WT (1H and 15N) and D14C (1H only) Fds with an intact disulfide bond and to the 4Fe oxidized WT Fd (1H and 15N) with a cleaved disulfide bond. All forms are shown to possess a long (13-member) alpha-helix, two beta-sheets (one double-, one triple-stranded), and three turns outside the cluster vicinity, each with tertiary contacts among themselves as found in other Fds. While the same secondary structural elements, with similar tertiary contacts, are found in other hyperthermostable Fds, Pf Fd has two elements, the long helix and the triple-stranded beta-sheet, that exhibit extensions and form multiple tertiary contacts. All Pf Fd forms with an intact disulfide bond exhibit a dynamic equilibrium heterogeneity which is shown to modulate a hydrogen-bonding network in the hydrophobic core that radiates from the Cys21-Cys48 disulfide bond and encompasses residues Lys36, Val24, Cys21, and Cys17 and the majority of the long helix. The heterogeneity is attributed to population of the alternate S and R chiralities of the disulfide bond, each destabilized by steric interactions with the extended alpha-helix. Comparison of the chemical shifts and their temperature gradients reveals that the molecular structure of the protein with the less stable R disulfide resembles that of the Fd with a cleaved disulfide bond. Both cluster architecture (3Fe vs 4Fe) and ligand mutation (Cys for Asp14) leave the disulfide orientational heterogeneity largely unperturbed. It is concluded that the six- to seven-residue extension that results in a longer helix and larger beta-sheet in Pf Fd, relative to other hyperthermostable Fds, more likely serves to destabilize the disulfide bond, and hence make it more readily reducible, than to significantly increase protein thermostability.

Role of the Invariant Peptide Fragment Forming NH.S Hydrogen Bonds in the Active Site of Cytochrome P-450 and Chloroperoxidase: Synthesis and Properties of Cys-Containing Peptide Fe(III) and Ga(III) (Octaethylporphinato) Complexes as Models

The primary sequence of Cys-X-Gly-Y- (X, hydrophobic residue; Y, hydrophilic residue) is highly conserved in cytochrome P-450s. The amide NHs of Leu, Gly, and X are assumed to form NH.S hydrogen bonds which are also found in the active site fragment, Cys-Pro-Ala-Leu, of chloroperoxidase (CPO). [Fe(III)(OEP)(Z-cys-Leu-Gly-Leu-OMe)] (OEP, octaethylporphinato; Z, benzyloxycarbonyl) and [Fe(III)(OEP)(Z-cys-Pro-Ala-Leu-OMe)] were synthesized as P-450 and CPO model complexes containing the invariant amino acid fragment of the active site. The corresponding gallium(III) complexes were also synthesized to investigate the solution structures using two-dimensional (2D) NMR experiments because the Ga(III) ion is similar to the Fe(III) ion in the ionic radii and in the coordination geometry. The solution structures of the peptide part of the gallium complexes indicate that the invariant fragments maintain a beta I-turn-like conformation and then form NH.S hydrogen bonds between S(gamma)Cys and NH of the third and fourth amino acid residues. The hydrogen bonds have also been confirmed by the (2)H NMR spectra of N(2)H-substituted Fe(III) peptide complexes. The Fe(III)/Fe(II) redox potentials of the Fe(III) complexes indicate that the NH.S hydrogen bonds in the fragments causes a slight positive shift of the redox potential. The tri- and tetrapeptide Fe(III) complexes containing the invariant fragments of P-450 are kinetically stable at 30 degrees C in CH(2)Cl(2). In contrast, [Fe(III)(OEP)(Z-cys-Leu-OMe)] decomposed to give [Fe(II)(OEP)] (22%) and the corresponding disulfide immediately in CD(2)Cl(2) at 30 degrees C for 1 h. These results indicate that the invariant fragments involving the hydrogen bonds cause the stabilization of the high-spin Fe(III) resting state rather than the positive shift of Fe(III)/Fe(II) redox potential.

The structure of iron-sulfur proteins

Ferredoxins are a group of iron-sulfur proteins for which a wealth of structural and mutational data have recently become available. Previously unknown structures of ferredoxins which are adapted to halophilic, acidophilic or hyperthermophilic environments and new cysteine patterns for cluster ligation and non-cysteine cluster ligation have been described. Site-directed mutagenesis experiments have given insight into factors that influence the geometry, stability, redox potential, electronic properties and electron-transfer reactivity of iron-sulfur clusters.

Determination of the structure of oxidised Desulfovibrio africanus ferredoxin I by 1H NMR spectroscopy and comparison of its solution structure with its crystal structure

The solution structure of the 64 amino acid Fe4S4 ferredoxin I from Desulfovibrio africanus has been determined using two-dimensional 1H NMR spectroscopy. Sequence-specific assignments were obtained for 59 amino acid residues and the structure determined with the program DIANA on the basis of 549 nuclear Overhauser enhancement (NOE) upper distance limits, and four dihedral angle and 52 distance constraints for the Fe4S4 cluster. The NMR structure was refined using the simulated annealing and energy minimisation protocols of the program X-PLOR to yield a final family of 19 structures selected on the basis of good covalent geometry and minimal restraint violations. The r.m.s.d. values to the average structure for this family are 0.49(+/-0.07) A and 0.94(+/-0.09) A for the backbone and heavy-atoms of residues 3 to 62, respectively. The NMR structure has been compared to the previously reported X-ray structures for the two molecules within the asymmetric unit of the crystal, which have a network of seven hydrogen bonds between them. This intermolecular interface, involving residues 38, 40 to 43 and 46, has the same conformation in the solution structures showing that the crystal packing does not perturb the structure. There are three regions in which the NMR and X-ray structures differ: around the cluster, a turn involving residues 8 to 10, and a loop involving residues 29 to 32. In the family of solution structures the backbone of the loop region incorporating residues 29 to 32 is well-defined whilst in both of the X-ray molecules it is ill-defined. The small differences between the X-ray and NMR structures for the cluster environment and the turn between residues 8 to 10 probably reflects a lack of NMR constraints. The observation of relatively rapid amide NH hydrogen exchange of NH groups close to the cluster, together with rapid flipping for Phe25, which is also close to the cluster, indicates that the cluster environment is more dynamic than the corresponding regions of related Fe/S proteins.

Conformation Control of Peptides by Metal Ions. Coordination Conformation Correlation Observed in a Model for Cys-X-Y-Cys/M(2+) in Proteins

The structure of [(Boc-Cys(1)-Pro-Leu-Cys(4)-OMe)(S-tert-C(4)H(9))Hg](-) (Boc: butoxycarbonyl), 1, was studied in N,N-dimethylformamide (DMF) and compared with that of [(Boc-Cys(1)-Pro-Leu-Cys(4)-OMe)Hg], 2, in order to discuss the intrinsic structural feature of the cysteine-containing metal-binding sites of proteins: Cys(i)-X-Y-Cys(i)(+3)/M(2+). 1 was generated by the reaction of 2 with NaS-tert-C(4)H(9). The geometry of the mercury ion (Hg(2+)) in 1 was proposed to be trigonal planar by UV-vis spectroscopy and Hg L(III) edge X-ray absorption fine structure (XAFS) measurements. Extended X-ray absorption fine structure (EXAFS) calculations yielded r(Hg-S) = 2.42 Å. Analyses of the nuclear Overhauser and exchange spectroscopy (NOESY) and the rotating frame nuclear Overhauser effect spectroscopy (ROESY) spectra of 1 in DMF-d(7) gave approximate distances for the 21 (1)H-(1)H pairs of the main chain loop. These results on distance information were processed by distance geometry (DG) and restrained molecular dynamics (RMD) calculations in order to optimize the molecular structure of 1. Molecular dynamics (MD) calculations were also performed. We proposed that the trigonal planar Hg(2+) in 1 regulates the hydrogen-bonding schemes of the peptide in the same manner as the tetrahedral ions involved in the Cys(i)-X-Y-Cys(i)(+3)/M(2+) core sites in natural proteins, forming two hydrogen bonds, Cys(1) S-Leu H(N) and Cys(1) S-Cys(4) H(N). This is in contrast to 2, where the linear coordinate mercury causes another type of hydrogen-bonding scheme, Cys(1) S-Leu H(N) and Pro CO-Cys(4) H(N). Details of the effect of trigonal planar Hg(2+) on the peptide conformation were analyzed with respect to the phi, varphi, and chi torsion angles of the peptide chain. The effect of the change of the angleS-Hg-S bite angle on the conformation of Cys-Pro-Leu-Cys was also discussed on the basis of MD calculations. The distribution area of Leu (phi, varphi) in the Ramachandran plot moves from near the alpha helix region to the turn structure region as the bite angle increases from 90 to 180 degrees, accompanying the change in the hydrogen-bonding scheme. The critical bite angle is around 140 degrees. The analysis revealed that angleS-Hg-S congruent with 110 degrees, which corresponds to the tetrahedral coordination geometry of the central metal ion, allows a high flexibility of the Cys-Pro-Leu-Cys skeleton.

Topological chirality of iron-sulfur proteins

An examination of x-ray structures of single-cluster [4Fe-4S] proteins in the Protein Data Bank has revealed that all redox proteins and the glutamine 5-phosphoribosyl-l-pyrophosphate amidotransferase from Bacillus subtilis have a topological configuration arbitrarily designated as D, whereas the DNA repair enzyme endonuclease III from Escherichia coli has the opposite topological configuration L. This is the first example in which both senses of topological chirality have been observed in a class of proteins.

Novel zinc-containing ferredoxin family in thermoacidophilic archaea

The dicluster-type ferredoxins from the thermoacidophilic archaea such as Thermoplasma acidophilum and Sulfolobus sp. are known to contain an unusually long extension of unknown function in the N-terminal region. Recent x-ray structural analysis of the Sulfolobus ferredoxin has revealed the presence of a novel zinc center, which is coordinated by three histidine ligand residues in the N-terminal region and one aspartate in the ferredoxin core domain. We report here the quantitative metal analyses together with electron paramagnetic resonance and resonance Raman spectra of T. acidophilum ferredoxin, demonstrating the presence of a novel zinc center in addition to one [3Fe-4S] and one [4Fe-4S] cluster (Fe/Zn = 6.8 mol/mol). A phylogenetic tree constructed for several archaeal monocluster and dicluster type ferredoxins suggests that the zinc-containing ferredoxins of T. acidophilum and Sulfolobus sp. form an independent subgroup, which is more distantly related to the ferredoxins from the hyperthermophiles than those from the methanogenic archaea, indicating the existence of a novel group of ferredoxins, namely, a "zinc-containing ferredoxin family" in the thermoacidophilic archaea. Inspection of the N-terminal extension regions of the archaeal zinc-containing ferredoxins suggested strict conservation of three histidine and one aspartate residues as possible ligands to the novel zinc center.

Small structural changes account for the high thermostability of 1[4Fe-4S] ferredoxin from the hyperthermophilic bacterium Thermotoga maritima

BACKGROUND

The characterization of the structural features that account for the high thermostability of some proteins is of great scientific and biotechnological interest. Proteins from hyperthermophilic organisms with optimum growth temperatures of 80 degrees C and higher generally show high intrinsic stabilities. The comparison of high resolution X-ray structures of these proteins with their counterparts from mesophilic organisms has therefore helped to identify potentially stabilizing forces in a number of cases. Small monomeric proteins which comprise only a single domain, such as ferredoxins, are especially suitable for such comparisons since the search for determinants of protein stability is considerably simplified.

RESULTS

The 1.75 A crystal structure of the extremely thermostable 1[4Fe-4S] ferredoxin from Thermotoga maritima (FdTm) was determined and compared with other monocluster-containing ferredoxins with different degrees of thermostability.

CONCLUSIONS

A comparison of the three-dimensional structure of FdTm with that of ferredoxins from mesophilic organisms suggests that the very high thermostability of FdTm is unexpectedly achieved without large changes of the overall protein structure. Instead, an increased number of potentially stabilizing features is observed in FdTm, compared with mesophilic ferredoxins. These include stabilization of alpha helices, replacement of residues in strained conformation by glycines, strong docking of the N-terminal methionine and an overall increase in the number of hydrogen bonds. Most of these features stabilize several secondary structure elements and improve the overall rigidity of the polypeptide backbone. The decreased flexibility will certainly play a relevant role in shielding the iron-sulfur cluster against physiologically high temperatures and further improve the functional integrity of FdTm.

Crystal structure of the 2[4Fe-4S] ferredoxin from Chromatium vinosum: evolutionary and mechanistic inferences for [3/4Fe-4S] ferredoxins

The crystal structure of the 2[4Fe-4S] ferredoxin from Chromatium vinosum has been solved by molecular replacement using data recorded with synchrotron radiation. The crystals were hexagonal prisms that showed a strong tendency to develop into long tubes. The hexagonal prisms diffracted to 2.1 A resolution at best, and a structural model for C. vinosum ferredoxin has been built with a final R of 19.2%. The N-terminal domain coordinates the two [4Fe-4S] clusters in a fold that is almost identical to that of other known ferredoxins. However, the structure has two unique features. One is a six-residue insertion between two ligands of one cluster forming a two-turn external loop; this short loop changes the conformation of the Cys 40 ligand compared to other ferredoxins and hampers the building of one NH...S H-bond to one of the inorganic sulfurs. The other remarkable structural element is a 3.5-turn alpha-helix at the C-terminus that covers one side of the same cluster and is linked to the cluster-binding domain by a six-residue external chain segment. The charge distribution is highly asymmetric over the molecule. The structure of C. vinosum ferredoxin strongly suggests divergent evolution for bacterial [3/4Fe-4S] ferredoxins from a common ancestral cluster-binding core. The unexpected slow intramolecular electron transfer rate between the clusters in C. vinosum ferredoxin, compared to other similar proteins, may be attributed to the unusual electronic properties of one of the clusters arising from localized changes in its vicinity rather than to a global structural rearrangement.

An NMR-derived model for the solution structure of oxidized Thermotoga maritima 1[Fe4-S4] ferredoxin

The solution structure of the 60-residue 1[Fe4-S4] ferredoxin from the hyperthermophilic bacterium Thermotoga maritima was determined based on 683 distance and 35 dihedral angle restraints that were obtained from NMR data. In addition, data known from crystallographic studies of ferredoxins was used for modeling of the iron-sulfur cluster and its environment. The protein shows a globular fold very similar to the fold of the related 1[Fe4-S4] ferredoxins from Desulfovibrio gigas and Desulfovibrio africanus, and elements of regular secondary structure similar to those in other ferredoxins were found in the T. maritima protein. In particular, the T. maritima protein displayed a beta-sheet structure made up of strands located at the very NH(2) and COOH termini of the protein, and an internal alpha-helix. The internal beta-sheet observed in the D. gigas and D. africanus ferredoxins could not be confirmed in T. maritima ferredoxin and is thus suggested to be only weakly present or even absent in this protein. This result suggests that thermostability in ferredoxins is not necessarily correlated with the content of stable elements of regular secondary structure.

The prosegment-subtilisin BPN' complex: crystal structure of a specific 'foldase'

BACKGROUND

The folding of the bacterial protease subtilisin BPN' (SBT) is dependent on its 77-residue prosegment, which is then autocatalytically removed to give the mature enzyme. Mature subtilisin represents a class of proteins that lacks an efficient folding pathway. Refolding of mature SBT is extremely slow unless catalyzed by the independently expressed prosegment, leading to a bimolecular complex.

RESULTS

We report the crystal structure at 2.0 A resolution of the prosegment-SBT complex and consider its implications for prosubtilisin BPN' maturation and folding catalysis. The prosegment forms a compact domain that binds SBT through an extensive interface involving the enzyme's two parallel surface helices (residues 104-116 and 133-144), supplying negatively charged caps to the N termini of these helices. The prosegment C terminus binds in the enzyme active site in a product-like manner, with Tyr77 in the P1 binding pocket.

CONCLUSIONS

The structure of the complex supports a unimolecular mechanism for prosubtilisin cleavage, involving a 25 A rearrangement of the SBT N terminus in a late folding step. A mechanism of folding catalysis in which the two helices and their connecting beta strand form a prosegment-stabilized folding nucleus is proposed. While this putative nucleus is stabilized by prosegment binding, the N-terminal and C-terminal subdomains of SBT could fold by propagation.

Determination of the [Fe4S4]Cys4 cluster geometry of Desulfovibrio africanus ferredoxin I by 1H NMR spectroscopy

1D and 2D 1H NMR studies of the Fe4S4 cluster containing ferredoxin I from Desulfovibrio africanus have been carried out with the aim of determining the geometry of the cluster linkages with the 4 Cys side chains that bind the cluster. This required the Cys beta CH resonances of the oxidised protein to be sequence-specifically and stereo-specifically assigned, and this was accomplished by a combination of TOCSY and NOE measurements, allied to model building based on X-ray structures of related ferredoxins. An analysis of the estimated hyperfine shifts of the Cys beta CH resonances with a Karplus-type equation relating the shifts to iron-sulfur-beta carbon-beta proton dihedral angles, taken together with the relative relaxation rates of the two beta CH2 resonances, estimated from their linewidths, then allowed the iron-sulfur-beta-carbon-alpha-carbon dihedral angles to be determined. A novel representation of the NMR data is presented which shows that the cluster dihedral angles are uniquely determined by the NMR data. The analysis reveals that the dihedral angles for D. africanus ferredoxin I are similar to the corresponding angles of other ferredoxins even though there are differences in their 1H NMR spectra. The sequence-specific and stereospecific assignments have been extended by analogy to the related Fe4S4-containing D. gigas ferredoxin I, and the stereospecific assignments to the Fe4S4-containing Thermococcus litoralis ferredoxin.