Browsing gene banks for Fe2S2 ferredoxins and structural modeling of 88 plant-type sequences: an analysis of fold and function

Predicting Kinetics and Dynamics of Spin-Dependent Processes

ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.

Fine-tuning of FeS proteins monitored via pulsed EPR redox potentiometry at Q-band

As essential electron translocating proteins in photosynthetic organisms, multiple plant-type ferredoxin (Fdx) isoforms are involved in a high number of reductive metabolic processes in the chloroplast. To allow quick cellular responses under changing environmental conditions, different plant-type Fdxs in Chlamydomonas reinhardtii were suggested to have adapted their midpoint potentials to a wide range of interaction partners. We performed pulsed electron paramagnetic resonance (EPR) monitored redox potentiometry at Q-band on three Fdx isoforms for a straightforward determination of their midpoint potentials. Additionally, site-directed mutagenesis was used to tune the midpoint potential of CrFdx1 in a range of approximately -338 to -511 mV, confirming the importance of single positions in the protein environment surrounding the [2Fe2S] cluster. Our results present a new target for future studies aiming to modify the catalytic activity of CrFdx1 that plays an essential role either as electron acceptor of photosystem I or as electron donor to hydrogenases under certain conditions. Additionally, the precisely determined redox potentials in this work using pulsed EPR demonstrate an alternative method that provides additional advantages compared with the well-established continuous wave EPR technique.

Diversification of Ferredoxins across Living Organisms

Ferredoxins, iron-sulfur (Fe-S) cluster proteins, play a key role in oxidoreduction reactions. To date, evolutionary analysis of these proteins across the domains of life have been confined to observing the abundance of Fe-S cluster types (2Fe-2S, 3Fe-4S, 4Fe-4S, 7Fe-8S (3Fe-4s and 4Fe-4S) and 2[4Fe-4S]) and the diversity of ferredoxins within these cluster types was not studied. To address this research gap, here we propose a subtype classification and nomenclature for ferredoxins based on the characteristic spacing between the cysteine amino acids of the Fe-S binding motif as a subtype signature to assess the diversity of ferredoxins across the living organisms. To test this hypothesis, comparative analysis of ferredoxins between bacterial groups, Alphaproteobacteria and Firmicutes and ferredoxins collected from species of different domains of life that are reported in the literature has been carried out. Ferredoxins were found to be highly diverse within their types. Large numbers of alphaproteobacterial species ferredoxin subtypes were found in Firmicutes species and the same ferredoxin subtypes across the species of Bacteria, Archaea, and Eukarya, suggesting shared common ancestral origin of ferredoxins between Archaea and Bacteria and lateral gene transfer of ferredoxins from prokaryotes (Archaea/Bacteria) to eukaryotes. This study opened new vistas for further analysis of diversity of ferredoxins in living organisms.

Distinct Physiological Roles of the Three Ferredoxins Encoded in the Hyperthermophilic Archaeon Thermococcus kodakarensis

High-energy electrons liberated during catabolic processes can be exploited for energy-conserving mechanisms. Maximal energy gains demand these valuable electrons be accurately shuttled from electron donor to appropriate electron acceptor. Proteinaceous electron carriers such as ferredoxins offer opportunities to exploit specific ferredoxin partnerships to ensure that electron flux to critical physiological pathways is aligned with maximal energy gains. Most species encode many ferredoxin isoforms, but very little is known about the role of individual ferredoxins in most systems. Our results detail that ferredoxin isoforms make largely unique and distinct protein interactions in vivo and that flux through one ferredoxin often cannot be recovered by flux through a different ferredoxin isoform. The results obtained more broadly suggest that ferredoxin isoforms throughout biological life have evolved not as generic electron shuttles, but rather serve as selective couriers of valuable low-potential electrons from select electron donors to desirable electron acceptors.

Control of electron flux is critical in both natural and bioengineered systems to maximize energy gains. Both small molecules and proteins shuttle high-energy, low-potential electrons liberated during catabolism through diverse metabolic landscapes. Ferredoxin (Fd) proteins—an abundant class of Fe-S-containing small proteins—are essential in many species for energy conservation and ATP production strategies. It remains difficult to model electron flow through complicated metabolisms and in systems in which multiple Fd proteins are present. The overlap of activity and/or limitations of electron flux through each Fd can limit physiology and metabolic engineering strategies. Here we establish the interplay, reactivity, and physiological role(s) of the three ferredoxin proteins in the model hyperthermophile Thermococcus kodakarensis. We demonstrate that the three loci encoding known Fds are subject to distinct regulatory mechanisms and that specific Fds are utilized to shuttle electrons to separate respiratory and energy production complexes during different physiological states. The results obtained argue that unique physiological roles have been established for each Fd and that continued use of T. kodakarensis and related hydrogen-evolving species as bioengineering platforms must account for the distinct Fd partnerships that limit flux to desired electron acceptors. Extrapolating our results more broadly, the retention of multiple Fd isoforms in most species argues that specialized Fd partnerships are likely to influence electron flux throughout biology.

Characterization of photosynthetic ferredoxin from the Antarctic alga Chlamydomonas sp. UWO241 reveals novel features of cold adaptation

The objective of this work was to characterize photosynthetic ferredoxin from the Antarctic green alga Chlamydomonas sp. UWO241, a key enzyme involved in distributing photosynthetic reducing power. We hypothesize that ferredoxin possesses characteristics typical of cold-adapted enzymes, namely increased structural flexibility and high activity at low temperatures, accompanied by low stability at moderate temperatures. To address this objective, we purified ferredoxin from UWO241 and characterized the temperature dependence of its enzymatic activity and protein conformation. The UWO241 ferredoxin protein, RNA, and DNA sequences were compared with homologous sequences from related organisms. We provide evidence for the duplication of the main ferredoxin gene in the UWO241 nuclear genome and the presence of two highly similar proteins. Ferredoxin from UWO241 has both high activity at low temperatures and high stability at moderate temperatures, representing a novel class of cold-adapted enzymes. Our study reveals novel insights into how photosynthesis functions in the cold. The presence of two distinct ferredoxin proteins in UWO241 could provide an adaptive advantage for survival at cold temperatures. The primary amino acid sequence of ferredoxin is highly conserved among photosynthetic species, and we suggest that subtle differences in sequence can lead to significant changes in activity at low temperatures.

FdC1 and Leaf-Type Ferredoxins Channel Electrons From Photosystem I to Different Downstream Electron Acceptors

Plant-type ferredoxins in Arabidopsis transfer electrons from the photosystem I to multiple redox-driven enzymes involved in the assimilation of carbon, nitrogen, and sulfur. Leaf-type ferredoxins also modulate the switch between the linear and cyclic electron routes of the photosystems. Recently, two novel ferredoxin homologs with extra C-termini were identified in the Arabidopsis genome (AtFdC1, AT4G14890; AtFdC2, AT1G32550). FdC1 was considered as an alternative electron acceptor of PSI under extreme ferredoxin-deficient conditions. Here, we showed that FdC1 could interact with some, but not all, electron acceptors of leaf-type Fds, including the ferredoxin-thioredoxin reductase (FTR), sulfite reductase (SiR), and nitrite reductase (NiR). Photoreduction assay on cytochrome c and enzyme assays confirmed its capability to receive electrons from PSI and donate electrons to the Fd-dependent SiR and NiR but not to the ferredoxin-NADP⁺ oxidoreductase (FNR). Hence, FdC1 and leaf-type Fds may play differential roles by channeling electrons from photosystem I to different downstream electron acceptors in photosynthetic tissues. In addition, the median redox potential of FdC1 may allow it to receive electrons from FNR in non-photosynthetic plastids.

Crystal structure and biochemical characterization of Chlamydomonas FDX2 reveal two residues that, when mutated, partially confer FDX2 the redox potential and catalytic properties of FDX1

The green alga Chlamydomonas reinhardtii contains six plastidic [2Fe2S]-cluster ferredoxins (FDXs), with FDX1 as the predominant isoform under photoautotrophic growth. FDX2 is highly similar to FDX1 and has been shown to interact with specific enzymes (such as nitrite reductase), as well as to share interactors with FDX1, such as the hydrogenases (HYDA), ferredoxin:NAD(P) reductase I (FNR1), and pyruvate:ferredoxin oxidoreductase (PFR1), albeit performing at low catalytic rates. Here we report the FDX2 crystal structure solved at 1.18 Å resolution. Based on differences between the Chlorella fusca FDX1 and C. reinhardtii FDX2 structures, we generated and purified point-mutated versions of the FDX2 protein and assayed them in vitro for their ability to catalyze hydrogen and NADPH photo-production. The data show that structural differences at two amino acid positions contribute to functional differences between FDX1 and FDX2, suggesting that FDX2 might have evolved from FDX1 toward a different physiological role in the cell. Moreover, we demonstrate that the mutations affect both the midpoint potentials of the FDX and kinetics of the FNR reaction, possibly due to altered binding between FDX and FNR. An effect on H2 photo-production rates was also observed, although the kinetics of the reaction were not further characterized.

Mutation of FdC2 gene encoding a ferredoxin-like protein with C-terminal extension causes yellow-green leaf phenotype in rice

Ferredoxins (Fds) are small iron-sulfur proteins that mediate electron transfer in a wide range of metabolic reactions. Besides Fds, there is a type of Fd-like proteins designated as FdC, which have conserved elements of Fds, but contain a significant C-terminal extension. So far, only two FdC genes of Arabidopsis (Arabidopsis thaliana) have been identified in higher plants and thus the functions of FdC proteins remain largely unknown. In this study, we isolated a yellow-green leaf mutant, 501ys, in rice (Oryza sativa). The mutant exhibited yellow-green leaf phenotype and reduced chlorophyll level. The phenotype of 501ys was caused by mutation of a gene on rice chromosome 3. Map-based cloning of this mutant resulted in identification of OsFdC2 gene (LOC_Os03g48040) showing high identity with Arabidopsis FdC2 gene (AT1G32550). OsFdC2 was expressed most abundantly in leaves and its encoded protein was targeted to the chloroplast. In 501ys mutant, a missense mutation was detected in DNA sequence of the gene, resulting in an amino acid change in the encoded protein. The mutant phenotype was rescued by introduction of the wild-type gene. Therefore, we successfully identified FdC2 gene via map-based cloning approach, and demonstrated that mutation of this gene caused yellow-green leaf phenotype in rice.

Hidden relationships between metalloproteins unveiled by structural comparison of their metal sites

Metalloproteins account for a substantial fraction of all proteins. They incorporate metal atoms, which are required for their structure and/or function. Here we describe a new computational protocol to systematically compare and classify metal-binding sites on the basis of their structural similarity. These sites are extracted from the MetalPDB database of minimal functional sites (MFSs) in metal-binding biological macromolecules. Structural similarity is measured by the scoring function of the available MetalS2 program. Hierarchical clustering was used to organize MFSs into clusters, for each of which a representative MFS was identified. The comparison of all representative MFSs provided a thorough structure-based classification of the sites analyzed. As examples, the application of the proposed computational protocol to all heme-binding proteins and zinc-binding proteins of known structure highlighted the existence of structural subtypes, validated known evolutionary links and shed new light on the occurrence of similar sites in systems at different evolutionary distances. The present approach thus makes available an innovative viewpoint on metalloproteins, where the functionally crucial metal sites effectively lead the discovery of structural and functional relationships in a largely protein-independent manner.

Ectopically expressed sweet pepper ferredoxin PFLP enhances disease resistance to Pectobacterium carotovorum subsp. carotovorum affected by harpin and protease-mediated hypersensitive response in Arabidopsis

Plant ferredoxin-like protein (PFLP) is a photosynthesis-type ferredoxin (Fd) found in sweet pepper. It contains an iron-sulphur cluster that receives and delivers electrons between enzymes involved in many fundamental metabolic processes. It has been demonstrated that transgenic plants overexpressing PFLP show a high resistance to many bacterial pathogens, although the mechanism remains unclear. In this investigation, the PFLP gene was transferred into Arabidopsis and its defective derivatives, such as npr1 (nonexpresser of pathogenesis-related gene 1) and eds1 (enhanced disease susceptibility 1) mutants and NAHG-transgenic plants. These transgenic plants were then infected with the soft-rot bacterial pathogen Pectobacterium carotovorum subsp. carotovorum (Erwinia carotovora ssp. carotovora, ECC) to investigate the mechanism behind PFLP-mediated resistance. The results revealed that, instead of showing soft-rot symptoms, ECC activated hypersensitive response (HR)-associated events, such as the accumulation of hydrogen peroxide (H2 O2 ), electrical conductivity leakage and expression of the HR marker genes (ATHSR2 and ATHSR3) in PFLP-transgenic Arabidopsis. This PFLP-mediated resistance could be abolished by inhibitors, such as diphenylene iodonium (DPI), 1-l-trans-epoxysuccinyl-leucylamido-(4-guanidino)-butane (E64) and benzyloxycarbonyl-Val-Ala-Asp-fluoromethylketone (z-VAD-fmk), but not by myriocin and fumonisin. The PFLP-transgenic plants were resistant to ECC, but not to its harpin mutant strain ECCAC5082. In the npr1 mutant and NAHG-transgenic Arabidopsis, but not in the eds1 mutant, overexpression of the PFLP gene increased resistance to ECC. Based on these results, we suggest that transgenic Arabidopsis contains high levels of ectopic PFLP; this may lead to the recognition of the harpin and to the activation of the HR and other resistance mechanisms, and is dependent on the protease-mediated pathway.

Identification of global ferredoxin interaction networks in Chlamydomonas reinhardtii

Ferredoxins (FDXs) can distribute electrons originating from photosynthetic water oxidation, fermentation, and other reductant-generating pathways to specific redox enzymes in different organisms. The six FDXs identified in Chlamydomonas reinhardtii are not fully characterized in terms of their biological function. In this report, we present data from the following: (a) yeast two-hybrid screens, identifying interaction partners for each Chlamydomonas FDX; (b) pairwise yeast two-hybrid assays measuring FDX interactions with proteins from selected biochemical pathways; (c) affinity pulldown assays that, in some cases, confirm and even expand the interaction network for FDX1 and FDX2; and (d) in vitro NADP(+) reduction and H2 photo-production assays mediated by each FDX that verify their role in these two pathways. Our results demonstrate new potential roles for FDX1 in redox metabolism and carbohydrate and fatty acid biosynthesis, for FDX2 in anaerobic metabolism, and possibly in state transition. Our data also suggest that FDX3 is involved in nitrogen assimilation, FDX4 in glycolysis and response to reactive oxygen species, and FDX5 in hydrogenase maturation. Finally, we provide experimental evidence that FDX1 serves as the primary electron donor to two important biological pathways, NADPH and H2 photo-production, whereas FDX2 is capable of driving these reactions at less than half the rate observed for FDX1.

Plant type ferredoxins and ferredoxin-dependent metabolism

Ferredoxin (Fd) is a small [2Fe-2S] cluster-containing protein found in all organisms performing oxygenic photosynthesis. Fd is the first soluble acceptor of electrons on the stromal side of the chloroplast electron transport chain, and as such is pivotal to determining the distribution of these electrons to different metabolic reactions. In chloroplasts, the principle sink for electrons is in the production of NADPH, which is mostly consumed during the assimilation of CO2 . In addition to this primary function in photosynthesis, Fds are also involved in a number of other essential metabolic reactions, including biosynthesis of chlorophyll, phytochrome and fatty acids, several steps in the assimilation of sulphur and nitrogen, as well as redox signalling and maintenance of redox balance via the thioredoxin system and Halliwell-Asada cycle. This makes Fds crucial determinants of the electron transfer between the thylakoid membrane and a variety of soluble enzymes dependent on these electrons. In this article, we will first describe the current knowledge on the structure and function of the various Fd isoforms present in chloroplasts of higher plants and then discuss the processes involved in oxidation of Fd, introducing the corresponding enzymes and discussing what is known about their relative interaction with Fd.

MetalPDB: a database of metal sites in biological macromolecular structures

We present here MetalPDB (freely accessible at http://metalweb.cerm.unifi.it), a novel resource aimed at conveying the information available on the three-dimensional (3D) structures of metal-binding biological macromolecules in a consistent and effective manner. This is achieved through the systematic and automated representation of metal-binding sites in proteins and nucleic acids by way of Minimal Functional Sites (MFSs). MFSs are 3D templates that describe the local environment around the metal(s) independently of the larger context of the macromolecular structure embedding the site(s), and are the central objects of MetalPDB design. MFSs are grouped into equistructural (broadly defined as sites found in corresponding positions in similar structures) and equivalent sites (equistructural sites that contain the same metals), allowing users to easily analyse similarities and variations in metal–macromolecule interactions, and to link them to functional information. The web interface of MetalPDB allows access to a comprehensive overview of metal-containing biological structures, providing a basis to investigate the basic principles governing the properties of these systems. MetalPDB is updated monthly in an automated manner.

Iron-based redox centres of reductase and oxygenase components of phenol hydroxylase from A. radioresistens: a redox chain working at highly positive redox potentials

This is the first report of the direct electrochemistry of the reductase (PHR) and oxygenase (PHO) components of phenol hydroxylase from Acinetobacter radioresistens S13 studied by cyclic and differential pulse voltammetry. The PHR contains one 2Fe2S cluster and one FAD that mediate the transfer of electrons from NAD(P)H to the non-heme diiron cluster of PHO. Cyclic and differential pulse voltammetry (CV and DPV) on glassy carbon showed two redox pairs with midpoint potentials at +131.5 ± 13 mV and -234 ± 3 mV versus normal hydrogen electrode (NHE). The first redox couple is attributed to the FeS centre, while the second one corresponds to free FAD released by the protein. DPV scans on native and guanidinium chloride treated PHR highlighted the presence of a split signal (ΔE ≈ 100 mV) attributed to heterogeneous properties of the 2Fe2S cluster interacting with the electrode, possibly due to the presence of two protein conformers and consistently with the large peak-to-peak separation and the peak broadening observed in CV. DPV experiments on gold electrodes performed on PHO confirm a consistently higher reduction potential at +396 mV vs. NHE. The positive redox potentials measured by direct electrochemistry for the FeS cluster in PHR and for the non-heme diiron cluster of PHO show that the entire phenol hydroxylase system works at higher potentials than those reported for structurally similar enzymes, for example methane monooxygenases.

Redox chemistry of the Schizosaccharomyces pombe ferredoxin electron-transfer domain and influence of Cys to Ser substitutions

Schizosaccharomyces pombe (Sp) ferredoxin contains a C-terminal electron transfer protein ferredoxin domain (etp(Fd)) that is homologous to adrenodoxin. The ferredoxin has been characterized by spectroelectrochemical methods, and Mössbauer, UV-Vis and circular dichroism spectroscopies. The Mössbauer spectrum is consistent with a standard diferric [2Fe-2S](2+) cluster. While showing sequence homology to vertebrate ferredoxins, the E°' and the reduction thermodynamics for etp(Fd) (-0.392 V) are similar to plant-type ferredoxins. Relatively stable Cys to Ser derivatives were made for each of the four bound Cys residues and variations in the visible spectrum in the 380-450 nm range were observed that are characteristic of oxygen ligated clusters, including members of the [2Fe-2S] cluster IscU/ISU scaffold proteins. Circular dichroism spectra were similar and consistent with no significant structural change accompanying these mutations. All derivatives were active in an NADPH-Fd reductase cytochrome c assay. The binding affinity of Fd to the reductase was similar, however, V(max) reflecting rate limiting electron transfer was found to decrease ~13-fold. The data are consistent with relatively minor perturbations of both the electronic properties of the cluster following substitution of the Fe-bond S atom with O, and the electronic coupling of the cluster to the protein.

Dynamics of Rhodobacter capsulatus [2FE-2S] ferredoxin VI and Aquifex aeolicus ferredoxin 5 via nuclear resonance vibrational spectroscopy (NRVS) and resonance Raman spectroscopy

We have used (57)Fe nuclear resonance vibrational spectroscopy (NRVS) to study the Fe(2)S(2)(Cys)(4) sites in oxidized and reduced [2Fe-2S] ferredoxins from Rhodobacter capsulatus (Rc FdVI) and Aquifex aeolicus (Aa Fd5). In the oxidized forms, nearly identical NRVS patterns are observed, with strong bands from Fe-S stretching modes peaking around 335 cm(-1), and additional features observed as high as the B(2u) mode at approximately 421 cm(-1). Both forms of Rc FdVI have also been investigated by resonance Raman (RR) spectroscopy. There is good correspondence between NRVS and Raman frequencies, but because of different selection rules, intensities vary dramatically between the two kinds of spectra. For example, the B(3u) mode at approximately 288 cm(-1), attributed to an asymmetric combination of the two FeS(4) breathing modes, is often the strongest resonance Raman feature. In contrast, it is nearly invisible in the NRVS, as there is almost no Fe motion in such FeS(4) breathing. NRVS and RR analysis of isotope shifts with (36)S-substituted into bridging S(2-) ions in Rc FdVI allowed quantitation of S(2-) motion in different normal modes. We observed the symmetric Fe-Fe stretching mode at approximately 190 cm(-1) in both NRVS and RR spectra. At still lower energies, the NRVS presents a complex envelope of bending, torsion, and protein modes, with a maximum at 78 cm(-1). The (57)Fe partial vibrational densities of states (PVDOS) were interpreted by normal-mode analysis with optimization of Urey-Bradley force fields. Progressively more complex D(2h) Fe(2)S(2)S'(4), C(2h) Fe(2)S(2)(SCC)(4), and C(1) Fe(2)S(2)(Cys)(4) models were optimized by comparison with the experimental spectra. After modification of the CHARMM22 all-atom force field by the addition of refined Fe-S force constants, a simulation employing the complete protein structure was used to reproduce the PVDOS, with better results in the low frequency protein mode region. This process was then repeated for analysis of data on the reduced FdVI. Finally, the degree of collectivity was used to quantitate the delocalization of the dynamic properties of the redox-active Fe site. The NRVS technique demonstrates great promise for the observation and quantitative interpretation of the dynamical properties of Fe-S proteins.

Pattern of expression and substrate specificity of chloroplast ferredoxins from Chlamydomonas reinhardtii

Ferredoxin (Fd) is the major iron-containing protein in photosynthetic organisms and is central to reductive metabolism in the chloroplast. The Chlamydomonas reinhardtii genome encodes six plant type [Fe2S2] ferredoxins, products of PETF, FDX2-FDX6. We performed the functional analysis of these ferredoxins by localizing Fd, Fdx2, Fdx3, and Fdx6 to the chloroplast by using isoform-specific antibodies and monitoring the pattern of gene expression by iron and copper nutrition, nitrogen source, and hydrogen peroxide stress. In addition, we also measured the midpoint redox potentials of Fd and Fdx2 and determined the kinetic parameters of their reactions with several ferredoxin-interacting proteins, namely nitrite reductase, Fd:NADP+ oxidoreductase, and Fd:thioredoxin reductase. We found that each of the FDX genes is differently regulated in response to changes in nutrient supply. Moreover, we show that Fdx2 (Em = -321 mV), whose expression is regulated by nitrate, is a more efficient electron donor to nitrite reductase relative to Fd. Overall, the results suggest that each ferredoxin isoform has substrate specificity and that the presence of multiple ferredoxin isoforms allows for the allocation of reducing power to specific metabolic pathways in the chloroplast under various growth conditions.

A novel, anaerobically induced ferredoxin in Chlamydomonas reinhardtii

We have found the transcript of one of at least six ferredoxin encoding genes of the green alga Chlamydomonas reinhardtii, FDX5, strongly accumulating in anaerobiosis, indicating a vital role of the encoded protein in the anaerobic metabolism of the cells. According to absorption and electron paramagnetic resonance spectroscopy, Fdx5 is a plant-type [2Fe2S]-ferredoxin with a redox potential similar to that of the ferredoxin PetF. However, although Fdx5 seems to be located in the chloroplast, it is not able to photoreduce nicotinamide adenine dinucleotide phosphate (NADP(+)) via ferredoxin-NADP-reductase, nor to be an electron donor to the plastidic [FeFe]-hydrogenase HydA1. Thus, Fdx5 seems to have a special role in a yet to be identified anaerobic pathway.

Robust recognition of zinc binding sites in proteins

Metals play a variety of roles in biological processes, and hence their presence in a protein structure can yield vital functional information. Because the residues that coordinate a metal often undergo conformational changes upon binding, detection of binding sites based on simple geometric criteria in proteins without bound metal is difficult. However, aspects of the physicochemical environment around a metal binding site are often conserved even when this structural rearrangement occurs. We have developed a Bayesian classifier using known zinc binding sites as positive training examples and nonmetal binding regions that nonetheless contain residues frequently observed in zinc sites as negative training examples. In order to allow variation in the exact positions of atoms, we average a variety of biochemical and biophysical properties in six concentric spherical shells around the site of interest. At a specificity of 99.8%, this method achieves 75.5% sensitivity in unbound proteins at a positive predictive value of 73.6%. We also test its accuracy on predicted protein structures obtained by homology modeling using templates with 30%-50% sequence identity to the target sequences. At a specificity of 99.8%, we correctly identify at least one zinc binding site in 65.5% of modeled proteins. Thus, in many cases, our model is accurate enough to identify metal binding sites in proteins of unknown structure for which no high sequence identity homologs of known structure exist. Both the source code and a Web interface are available to the public at http://feature.stanford.edu/metals.

Iron-sulfur protein folds, iron-sulfur chemistry, and evolution

An inventory of unique local protein folds around Fe-S clusters has been derived from the analysis of protein structure databases. Nearly 50 such folds have been identified, and over 90% of them harbor low-potential [2Fe-2S](2+,+) or [4Fe-4S](2+,+) clusters. In contrast, high-potential Fe-S clusters, notwithstanding their structural diversity, occur in only three different protein folds. These observations suggest that the extant population of Fe-S protein folds has to a large extent been shaped in the reducing iron- and sulfur-rich environment that is believed to have predominated on this planet until approximately two billion years ago. High-potential active sites are then surmised to be rarer because they emerged later, in a more oxidizing biosphere, in conditions where iron and sulfide had become poorly available, Fe-S clusters were less stable, and in addition faced competition from heme iron and copper active sites. Among the low-potential Fe-S active sites, protein folds hosting [4Fe-4S](2+,+) clusters outnumber those with [2Fe-2S](2+,+) ones by a factor of 3 at least. This is in keeping with the higher chemical stability and versatility of the tetranuclear clusters, compared with the binuclear ones. It is therefore suggested that, at least while novel Fe-S sites are evolving within proteins, the intrinsic chemical stability of the inorganic moiety may be more important than the stabilizing effect of the polypeptide chain. The discovery rate of novel Fe-S-containing protein folds underwent a sharp increase around 1995, and has remained stable to this day. The current trend suggests that the mapping of the Fe-S fold space is not near completion, in agreement with predictions made for protein folds in general. Altogether, the data collected and analyzed here suggest that the extant structural landscape of Fe-S proteins has been shaped to a large extent by primeval geochemical conditions on one hand, and iron-sulfur chemistry on the other.

Non-heme iron through the three domains of life

Metalloproteins are proteins capable of binding one or more metal ions, which are often required for their biological function or for regulation of their activities or for structural purposes. In high-throughput genome-level protein investigation efforts, such as Structural Genomics, the systematic experimental characterization of metal-binding properties (i.e. the investigation of the metalloproteome) is not always pursued, and remains far from trivial. In the present work we have applied a bioinformatic approach to investigate the occurrence of (putative) non-heme iron-binding proteins in 57 different organisms spanning the entire tree of life. It is found that the non-heme iron-proteome constitutes between 1% and 10% of the entire proteome of an organism. However, the iron-proteome constitutes a higher fraction of the proteome in archaea (on average 7.1% +/- 2.1%) than in bacteria (3.9% +/- 1.6%) and in eukaryota (1.1% +/- 0.4%). The analysis of the function of each putative iron-protein identified suggests that extant organisms have inherited the large majority of their iron-proteome from the last common ancestor.

Disease resistance to bacterial pathogens affected by the amount of ferredoxin-I protein in plants

SUMMARY Ferredoxin-I (Fd-I) is a fundamental protein that is involved in several metabolic pathways. The amount of Fd-I found in plants is generally regulated by environmental stress, including biotic and abiotic events. In this study, the correlation between quantity of Fd-I and plant disease resistance was investigated. Fd-I levels were increased by inoculation with Pseudomonas syringae pv. syringae but were reduced by Erwinia carotovora ssp. carotovora. Transgenic tobacco over-expressing Fd-I with the sense sweet pepper Fd-I gene (pflp) was resistant to E. carotovora ssp. carotovora and the saprophytic bacterium P. fluorescens. By contrast, transgenic tobacco with reduced total Fd-I and the antisense pflp gene was susceptible to E. carotovora ssp. carotovora and P. fluorescens. Both of these transgenic tobaccos were resistant to P. syringae pv. syringae. By contrast, the mutated E. carotovora ssp. carotovora, with a defective harpin protein, was able to invade the sense-pflp transgenic tobacco as well as the non-transgenic tobacco. An in vitro kinase assay revealed that harpin could activate unidentified kinases to phosphorylate PFLP. These results demonstrate that Fd-I plays an important role in the disease defence mechanism.

The Expression of a Plant-type Ferredoxin Redox System provides Molecular Evidence for a Plastid in the Early Dinoflagellate Perkinsus marinus

Perkinsus marinus is a parasitic protozoan with a phylogenetic positioning between Apicomplexa and dinoflagellates. It is thus of interest for reconstructing the early evolution of eukaryotes, especially with regard to the acquisition of secondary plastids in these organisms. It is also an important pathogen of oysters, and the definition of parasite-specific metabolic pathways would be beneficial for the identification of efficient treatments for infected mollusks. Although these different scientific interests have resulted in the start of a genome project for this organism, it is still unknown whether P. marinus contains a plastid or plastid-like organelle like the related dinoflagellates and Apicomplexa. Here, we show that in vitro-cultivated parasites contain transcripts of the plant-type ferredoxin and its associated reductase. Both proteins are nuclear-encoded and possess N-terminal targeting sequences similar to those characterized in dinoflagellates. Since this redox pair is exclusively found in cyanobacteria and plastid-harboring organisms its presence also in P. marinus is highly indicative of a plastid. We also provide additional evidence for such an organelle by demonstrating pharmacological sensitivity to inhibitors of plastid-localized enzymes involved in fatty acid biosynthesis (e.g. acetyl-CoA carboxylase) and by detection of genes for three enzymes of plastid-localized isoprenoid biosynthesis (1-deoxy-D-xylulose 5-phosphate reductoisomerase, (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate reductase, and (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate synthase).

FAST-NMR: functional annotation screening technology using NMR spectroscopy

An abundance of protein structures emerging from structural genomics and the Protein Structure Initiative (PSI) are not amenable to ready functional assignment because of a lack of sequence and structural homology to proteins of known function. We describe a high-throughput NMR methodology (FAST-NMR) to annotate the biological function of novel proteins through the structural and sequence analysis of protein-ligand interactions. This is based on basic tenets of biochemistry where proteins with similar functions will have similar active sites and exhibit similar ligand binding interactions, despite global differences in sequence and structure. Protein-ligand interactions are determined through a tiered NMR screen using a library composed of compounds with known biological activity. A rapid co-structure is determined by combining the experimental identification of the ligand binding site from NMR chemical shift perturbations with the protein-ligand docking program AutoDock. Our CPASS (Comparison of Protein Active Site Structures) software and database are then used to compare this active site with proteins of known function. The methodology is demonstrated using unannotated protein SAV1430 from Staphylococcus aureus.

Counting the zinc-proteins encoded in the human genome

Metalloproteins are proteins capable of binding one or more metal ions, which may be required for their biological function, or for regulation of their activities or for structural purposes. Genome sequencing projects have provided a huge number of protein primary sequences, but, even though several different elaborate analyses and annotations have been enabled by a rich and ever-increasing portfolio of bioinformatic tools, metal-binding properties remain difficult to predict as well as to investigate experimentally. Consequently, the present knowledge about metalloproteins is only partial. The present bioinformatic research proposes a strategy to answer the question of how many and which proteins encoded in the human genome may require zinc for their physiological function. This is achieved by a combination of approaches, which include: (i) searching in the proteome for the zinc-binding patterns that, on their turn, are obtained from all available X-ray data; (ii) using libraries of metal-binding protein domains based on multiple sequence alignments of known metalloproteins obtained from the Pfam database; and (iii) mining the annotations of human gene sequences, which are based on any type of information available. It is found that 1684 proteins in the human proteome are independently identified by all three approaches as zinc-proteins, 746 are identified by two, and 777 are identified by only one method. By assuming that all proteins identified by at least two approaches are truly zinc-binding and inspecting the proteins identified by a single method, it can be proposed that ca. 2800 human proteins are potentially zinc-binding in vivo, corresponding to 10% of the human proteome, with an uncertainty of 400 sequences. Available functional information suggests that the large majority of human zinc-binding proteins are involved in the regulation of gene expression. The most abundant class of zinc-binding proteins in humans is that of zinc-fingers, with Cys4 and Cys2His2 being the most common types of coordination environment.

Paramagnetic NMR spectroscopy and density functional calculations in the analysis of the geometric and electronic structures of iron-sulfur proteins

Paramagnetic NMR spectroscopy has been underutilized in the study of metalloproteins. One difficulty of the technique is that paramagnetic relaxation broadens signals from nuclei near paramagnetic centers. In systems with low electronic relaxation rates, this makes such signals difficult to observe or impossible to assign by traditional methods. We show how the challenges of detecting and assigning signals from nuclei near the metal center can be overcome through the combination of uniform and selective 2H, 13C, and 15N isotopic labeling with NMR experiments that utilize direct one-dimensional (2H, 13C, and 15N) and two-dimensional (13C-X) detection. We have developed methods for calculating NMR chemical shifts and relaxation rates by density functional theory (DFT) approaches. We use the correspondence between experimental NMR parameters and those calculated from structural models of iron-sulfur clusters derived from X-ray crystallography to validate the computational approach and to investigate how structural differences are manifested in these values. We have applied this strategy to three iron-sulfur proteins: Clostridium pasteurianum rubredoxin, Anabaena [2Fe-2S] ferredoxin, and human [2Fe-2S] ferredoxin. Provided that an accurate structural model of the iron-sulfur cluster and surrounding residues is available from diffraction data, our results show that DFT calculations can return NMR observables with excellent accuracy. This suggests that it might be possible to use calculations to refine structures or to generate structural models of active sites when crystal structures are unavailable. The approach has yielded insights into the electronic structures of these iron-sulfur proteins. In rubredoxin, the results show that substantial unpaired electron spin is delocalized across NH...S hydrogen bonds and that the reduction potential can be changed by 77 mV simply by altering the strength of one of these hydrogen bonds. In reduced [2Fe-2S] ferredoxins, hyperfine shift data have provided quantitative information on the degree of valence trapping. The approach described here for iron-sulfur proteins offers new avenues for detailed studies of these and other metalloprotein systems.

Strategy for the study of paramagnetic proteins with slow electronic relaxation rates by nmr spectroscopy: application to oxidized human [2Fe-2S] ferredoxin

NMR studies of paramagnetic proteins are hampered by the rapid relaxation of nuclei near the paramagnetic center, which prevents the application of conventional methods to investigations of the most interesting regions of such molecules. This problem is particularly acute in systems with slow electronic relaxation rates. We present a strategy that can be used with a protein with slow electronic relaxation to identify and assign resonances from nuclei near the paramagnetic center. Oxidized human [2Fe-2S] ferredoxin (adrenodoxin) was used to test the approach. The strategy involves six steps: (1) NMR signals from (1)H, (13)C, and (15)N nuclei unaffected or minimally affected by paramagnetic effects are assigned by standard multinuclear two- and three-dimensional (2D and 3D) spectroscopic methods with protein samples labeled uniformly with (13)C and (15)N. (2) The very broad, hyperfine-shifted signals from carbons in the residues that ligate the metal center are classified by amino acid and atom type by selective (13)C labeling and one-dimensional (1D) (13)C NMR spectroscopy. (3) Spin systems involving carbons near the paramagnetic center that are broadened but not hyperfine-shifted are elucidated by (13)C[(13)C] constant time correlation spectroscopy (CT-COSY). (4) Signals from amide nitrogens affected by the paramagnetic center are assigned to amino acid type by selective (15)N labeling and 1D (15)N NMR spectroscopy. (5) Sequence-specific assignments of these carbon and nitrogen signals are determined by 1D (13)C[(15)N] difference decoupling experiments. (6) Signals from (1)H nuclei in these spin systems are assigned by paramagnetic-optimized 2D and 3D (1)H[(13)C] experiments. For oxidized human ferredoxin, this strategy led to assignments (to amino acid and atom type) for 88% of the carbons in the [2Fe-2S] cluster-binding loops (residues 43-58 and 89-94). These included complete carbon spin-system assignments for eight of the 22 residues and partial assignments for each of the others. Sequence-specific assignments were determined for the backbone (15)N signals from nine of the 22 residues and ambiguous assignments for five of the others.

Miraculous catch of iron-sulfur protein sequences in the Sargasso Sea

Recent shotgun sequencing of filtered Sargasso Sea water samples has yielded data in astounding amount and diversity. Iron-sulfur proteins, which are ancient, diverse and ubiquitous, have been implemented here to further probe the sequence diversity of the Sargasso Sea database (SSDB). Sequence searches and comparisons confirm that the SSDB by and large equals in diversity the combined currently available databases. The data thus suggest that microbial diversity has so far been underestimated by orders of magnitude.

A post genomic characterization of Arabidopsis ferredoxins

In higher plant plastids, ferredoxin (Fd) is the unique soluble electron carrier protein located in the stroma. Consequently, a wide variety of essential metabolic and signaling processes depend upon reduction by Fd. The currently available plant genomes of Arabidopsis and rice (Oryza sativa) contain several genes encoding putative Fds, although little is known about the proteins themselves. To establish whether this variety represents redundancy or specialized function, we have recombinantly expressed and purified the four conventional [2Fe-2S] Fd proteins encoded in the Arabidopsis genome and analyzed their physical and functional properties. Two proteins are leaf type Fds, having relatively low redox potentials and supporting a higher photosynthetic activity. One protein is a root type Fd, being more efficiently reduced under nonphotosynthetic conditions and supporting a higher activity of sulfite reduction. A further Fd has a remarkably positive redox potential and so, although redox active, is limited in redox partners to which it can donate electrons. Immunological analysis indicates that all four proteins are expressed in mature leaves. This holistic view demonstrates how varied and essential soluble electron transfer functions in higher plants are fulfilled through a diversity of Fd proteins.

Incorporation of iron into Tritrichomonas foetus cell compartments reveals ferredoxin as a major iron-binding protein in hydrogenosomes

The intracellular transport of iron and its incorporation into organelles are poorly understood processes in eukaryotes and virtually unknown in parasitic protists. The transport of iron is of particular interest in trichomonads, which possess hydrogenosomes instead of mitochondria. The metabolic functions of hydrogenosomes, which contain a specific set of FeS proteins, entirely depend on iron acquisition. In this work the incorporation of iron into the cattle parasite Tritrichomonas foetus was monitored. Iron was efficiently taken up from (59)Fe-nitrilotriacetic acid and accumulated in the cytosol (88.9 %) and hydrogenosomes (4.7 % of the total radioactivity). Using atomic absorption spectrophotometry, an unusually high steady-state iron concentration in hydrogenosomes was determined [54.4+/-1.1 nmol Fe (mg protein)(-1)]. The concentration of iron in the cytosol was 13.4+/-0.5 nmol Fe (mg protein)(-1). Qualitative analysis of incorporated iron was performed using native gradient PAGE. The majority of the (59)Fe in the cytosol appeared as the labile-iron pool, which represents weakly bound iron associated with compounds of molecular mass ranging from 5000 to 30000 Da. Ferritin was not observed in Tt. foetus, nor in two other anaerobic protists, Entamoeba histolytica and Giardia intestinalis. Analysis of Tt. foetus hydrogenosomes showed at least nine iron-binding compounds, which were absent in metronidazole-resistant mutants. The major iron-binding compound was identified as [2Fe-2S] ferredoxin of the adrenodoxin type.

Characterization of a [2Fe-2S] protein encoded in the iron-hydrogenase operon of Thermotoga maritima

Thermotoga maritima grows optimally at 80 degrees C by fermenting carbohydrates to organic acids, CO(2), and H(2). The production of H(2) is catalyzed by a cytoplasmic, heterotrimeric (alphabetagamma) Fe-hydrogenase. This is encoded by three genes, hydC (gamma), hydB (beta) and hydA (alpha), organized within a single operon that contains five additional open reading frames (ORFs). The recombinant form of the first ORF of the operon, TM1420, was produced in Escherichia coli. It has a molecular mass of 8537+/-3 Da as determined by mass spectrometry, in agreement with the predicted amino acid sequence. Purified TM1420 is red in color, has a basic p I (8.8), and contains 1.9 Fe atoms/mol that are present as a single [2Fe-2S] cluster, as determined by UV-visible absorption and EPR spectroscopy. The protein contains five cysteine residues, but their arrangement is characteristic of a subunit or domain rather than of a ferredoxin-type protein. The reduction potential of the [2Fe-2S] cluster (-233 mV at pH 6.5 and 25 degrees C) is pH independent but decreases linearly with temperature to -296 mV (-1.15 mV/ degrees C) at 80 degrees C. TM1420 is not reduced, in vitro, by the Fe-hydrogenase nor by a pyruvate ferredoxin oxidoreductase. The protein was unstable at 70 degrees C under anaerobic conditions with a half-life of approximately 30 min. The basic nature of TM1420, its instability at the growth temperature of T. maritima, and the unusual spacing of its cysteine residues suggest that this protein does not function as a ferredoxin-type electron carrier for the Fe-hydrogenase. Instead, TM1420 is more likely part of a thermostable multi-protein complex that is involved in metal cluster assembly of the hydrogenase holoenzyme.

Bioinorganic chemistry in the postgenomic era

Genome sequencing has revolutionized all fields of life sciences. Bioinorganic chemistry is certainly not immune to this influence, which is presenting unprecedented challenges. A new goal for bioinorganic chemistry is the investigation of the linkages between inorganic elements and genomic information. This requires new advancements andor the development of new expertise in fields such as bioinformatics and genetics but also provides a driving force to push forward the exploitation of traditional analytical techniques and spectroscopic tools. The "case study" of metal homeostasis in cells is discussed to provide a flavor of the current evolution of the field.

Eukaryotic genomes contain a [2Fez.sbnd;2S] ferredoxin isoform with a conserved C-terminal sequence motif

Apicomplexan protists contain a single mitochondrial [2Fe-2S] ferredoxin sequence (mtFd) with a highly conserved C-terminal motif, VDGxxpxPH, that distinguishes it from other mtFd, which have heterogeneous C-termini. This isoform of mtFd, called 'type II ferredoxin', is widespread in eukaryotes, some species having two isoforms and others possessing only one. Because of the known modulating role of the C-terminus of type I mtFd during association with itself and other interacting proteins, the presence of a conserved C-terminus in type II mtFd suggests it evolved either as a means for optimized homodimerization or to allow interaction with a highly conserved partner(s) that is yet to be defined.

Exceptional stability of a [2Fe-2S] ferredoxin from hyperthermophilic bacterium Aquifex aeolicus

Aquifex aeolicus is the only hyperthermophile that is known to contain a plant- and mammalian-type [2Fe-2S] ferredoxin (Aae Fd1). This unique protein contains two cysteines, in addition to the four that act as ligands of the [2Fe-2S] cluster, which form a disulfide bridge. We have investigated the stability of Aae Fd1 with (wild-type) and without (C87A variant) the disulfide bond, with respect to pH, thermal and chemical perturbation, and compared the results to those for the mesophilic [2Fe-2S] ferredoxin from spinach. Unfolding reactions of all three proteins are irreversible due to cluster decomposition in the unfolded state. Wild-type and C87A Aae Fd1 proteins are extremely stable: unfolding at 20 degrees C requires high concentrations of the chemical denaturant and long incubation times. Moreover, their thermal-unfolding midpoints are 40-50 degrees higher than that for spinach ferredoxin (pH 7). The stability of the Aae Fd1 protein is significantly lower at pH 2.5 than pH 7 and 10, suggesting that ionic interactions play a role in structural integrity. Interestingly, the iron-sulfur cluster in C87A Aae Fd1 rearranges into a transient species with absorption bands at 520 and 610 nm, presumably a linear three-iron cluster, in the high-pH unfolded state.

Biogenesis of iron-sulphur clusters in amitochondriate and apicomplexan protists

During the last 4 years there has been an enormous interest in the question how iron-sulphur ([Fe-S]) clusters, which are essential building blocks for life, are synthesised and assembled into apo-proteins, both in prokaryotes and in eukaryotes. The emerging picture is that the basic mechanism of this pathway has been well conserved during evolution. In yeast and probably all other eukaryotes the mitochondrion is the place where [Fe-S] clusters are synthesised, even for extramitochondrial [Fe-S] cluster-containing proteins, and a number of proteins have been functionally characterised to a certain extent within this pathway. However, almost nothing is known about this aspect in parasitic protists, although recent studies of amitochondriate protists and on the plastid-like organelle of apicomplexan parasites, the apicoplast, have started to change this. In this article I will summarise the current view of [Fe-S] cluster biogenesis in eukaryotes and discuss its implications for amitochondriate protists and for the plastid-like organelle of apicomplexan parasites.