Analyses of the folding properties of ferredoxin‐like fold proteins by means of a coarse‐grained Gō model: Relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi‐state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on‐going research issues. In this study, we investigate ferredoxin‐like fold proteins which have a simple topology and multi‐state transitions. We analyze the folding processes by means of a coarse‐grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high‐free‐energy intermediate state, and ADA2h and S6, which fold into the native structure through two‐state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54‐55, are reproduced and compared with experimental observations. We show that the ferredoxin‐like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two‐state mechanism. Proteins 2014; 82:954–965. © 2013 Wiley Periodicals, Inc.     

Structural comparison of the two alternative transition states for folding of TI I27

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.     

Conformational detours during folding of a collapsed state

The protein S6 is a useful model to probe the role of partially folded states in the folding process. In the absence of salt, S6 folds from the denatured state D to the native state N without detectable intermediates. High concentrations of sodium sulfate induce the accumulation of a collapsed state C, which is off the direct folding route. However, the mutation VA85 enables S6 to fold from C directly to N through the transition state TS(C). According to the denaturant dependence of this reaction, TS(C) and C are equally compact, but the data are difficult to deconvolute. Therefore, I have measured the heat capacities (DeltaC(p)) for the D-->C and C-->TS(C) transitions. The DeltaC(p)-values suggest that C needs to increase its surface area in order to fold directly to N. This underlines that it is a misfolded state that can only fold by at least partial unfolding. In contrast to the C-state formed by S6 wildtype, the VA85 C-state is just as compact as the native state, and this may be a prerequisite for direct folding. Individual "gatekeeper" residues may thus play a disproportionately large role in guiding proteins through different folding pathways.     

Reversible folding of cysteine-rich metallothionein by an overcritical reaction path

A first-order-like state transition is considered to be involved in the restoration of the activities of a few proteins by correctly folding the protein [Phys. Rev. E 66 (2002) 021903]. In order to understand the general applicability of this mechanism, we studied a metallothionein (MT) protein with an unconventional structure, i.e., without any alpha-helix or beta-sheet. MT is a 61 amino-acid peptide. There are 6-7 Zn(2+) ions, which bind avidly to 20 conserved cysteines (Cys) of MT. These properties indicate that the structure of MT is quite different from those of the other proteins. Similar to our previous findings, the denatured MT can be folded without any aggregation via a designated stepwise quasi-static process (an over-critical reaction path). The particle size of folded MT intermediates, determined by dynamic light scattering, shrank right after the first folding stage. It is consistent with a collapse-model. In addition, results from both atomic absorption and circular dichroism (CD) indicate that the stable intermediates may fold to the native conformation but with only partial Zn(2+) binding, which in turn implies that those folding intermediates are in a molten globular state. These reversible unfolding and folding processes indicate that Cys-rich protein, MT, may also be folded by way of a first-order-like state transition mechanism. We suspect that this process may likely be involved in the reaction of the metal substitution process in metal containing enzymes.     

YbiV from Escherichia coli K12 is a HAD phosphatase

The protein YbiV from Escherichia coli K12 MG1655 is a hypothetical protein with sequence homology to the haloacid dehalogenase (HAD) superfamily of proteins. Although numerous members of this family have been identified, the functions of few are known. Using the crystal structure, sequence analysis, and biochemical assays, we have characterized YbiV as a HAD phosphatase. The crystal structure of YbiV reveals a two-domain protein, one with the characteristic HAD hydrolase fold, the other an inserted alpha/beta fold. In an effort to understand the mechanism, we also solved and report the structures of YbiV in complex with beryllofluoride (BeF3-) and aluminum trifluoride (AlF3), which have been shown to mimic the phosphorylated intermediate and transition state for hydrolysis, respectively, in analogy to other HAD phosphatases. Analysis of the structures reveals the substrate-binding cavity, which is hydrophilic in nature. Both structure and sequence homology indicate YbiV may be a sugar phosphatase, which is supported by biochemical assays that measured the release of free phosphate on a number of sugar-like substrates. We also investigated available genomic and functional data in an effort to determine the physiological substrate.     

Stabilization of the ribosomal protein S6 by trehalose is counterbalanced by the formation of a putative off-pathway species

The effect of trehalose on folding and stability of the small ribosomal protein S6 was studied. Non-disruptive point mutations distributed along the protein structure were analyzed to characterize the stabilizing effect of trehalose and map the folding pathway of S6. On average, the stability of the wild-type and S6 mutants increases by 3 kcal/mol M trehalose. Despite the non-specific thermodynamic stabilization mechanism, trehalose particularly stabilizes the less destabilized mutants. Folding/unfolding kinetics shows clearly that trehalose induces the collapse of the unfolded state to an off-pathway intermediate with non-native diffuse contacts. This state is similar to the collapsed state induced by high concentrations of stabilizing salts, as previously reported. Although it leads to the accumulation of this off-pathway intermediate, trehalose does not change the compactness of the transition state ensemble. Furthermore, the productive folding pathway of S6 is not affected by trehalose as shown by a Phi-value analysis. The unfolded state ensemble of S6 should be more compact in the presence of trehalose and therefore destabilized due to decreased conformational entropy. Increased compaction of the unfolded state ensemble might also occur for more stable mutants of S6, thus explaining the synergistic effect of trehalose and point mutations on protein stabilization.     

Mutational analysis of acylphosphatase suggests the importance of topology and contact order in protein folding.

Muscle acylphosphatase (AcP) is a small protein that folds very slowly with two-state behavior. The conformational stability and the rates of folding and unfolding have been determined for a number of mutants of AcP in order to characterize the structure of the folding transition state. The results show that the transition state is an expanded version of the native protein, where most of the native interactions are partially established. The transition state of AcP turns out to be remarkably similar in structure to that of the activation domain of procarboxypeptidase A2 (ADA2h), a protein having the same overall topology but sharing only 13% sequence identity with AcP. This suggests that transition states are conserved between proteins with the same native fold. Comparison of the rates of folding of AcP and four other proteins with the same topology, including ADA2h, supports the concept that the average distance in sequence between interacting residues (that is, the contact order) is an important determinant of the rate of protein folding.     

Different circular permutations produced different folding nuclei in proteins: a computational study

There have been many studies about the effect of circular permutation on the transition state/folding nucleus of proteins, with sometimes conflicting conclusions from different proteins and permutations. To clarify this important issue, we have studied two circular permutations of a lattice protein model with side-chains. Both permuted sequences have essentially the same native state as the original (wild-type) sequence. Circular permutant 1 cuts at the folding nucleus of the wild-type sequence. As a result, the permutant has a drastically different nucleus and folds more slowly than wild-type. In contrast, circular permutant 2 involves an incision at a site unstructured in the wild-type transition state, and the wild-type nucleus is largely retained in the permutant. In addition, permutant 2 displays both two-state and multi-state folding, with a native-like intermediate state occasionally populated. Neither the wild-type nor permutant 1 has a similar intermediate, and both fold in an apparently two-state manner. Surprisingly, permutant 2 folds at a rate identical with that of the wild-type. The intermediate in permutant 2 is stabilised by native and non-native interactions, and cannot be classified simply as on or off-pathway. So we advise caution in attributing experimental data to on or off-pathway intermediates. Finally, our work illuminates the results on alpha-spectrin SH3, chymotrypsin inhibitor 2 and beta-lactoglobulin, and supports a key assumption in the experimental efforts to locate potential nucleation sites of real proteins via circular permutations.     

Mapping the folding pathway of an immunoglobulin domain: structural detail from Phi value analysis and movement of the transition state

BACKGROUND: Do proteins that have the same structure fold by the same pathway even when they are unrelated in sequence? To address this question, we are comparing the folding of a number of different immunoglobulin-like proteins. Here, we present a detailed protein engineering phi value analysis of the folding pathway of TI I27, an immunoglobulin domain from human cardiac titin. RESULTS: TI I27 folds rapidly via a kinetic intermediate that is destabilized by most mutations. The transition state for folding is remarkably native-like in terms of solvent accessibility. We use phi value analysis to map this transition state and show that it is highly structured; only a few residues close to the N-terminal region of the protein remain completely unfolded. Interestingly, most mutations cause the transition state to become less native-like. This anti-Hammond behavior can be used as a novel means of obtaining additional structural information about the transition state. CONCLUSIONS: The residues that are involved in nucleating the folding of TI I27 are structurally equivalent to the residues that form the folding nucleus in an evolutionary unrelated fibronectin type III protein. These residues form part of the common structural core of Ig-like domains. The data support the hypothesis that interactions essential for defining the structure of these beta sandwich proteins are also important in nucleation of folding.     

Stability, unfolding, and structural changes of cofactor-free 1H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase

Stability, unfolding mechanism, spectroscopic, densimetric, and structural characteristics of the oxidatively stable C69S variant (HodC) of 1H-3-hydroxy-4-oxoquinaldine 2,4-dioxygenase (Hod) have been determined by classical and pressure modulation scanning calorimetry (DSC and PMDSC, respectively), circular dichroism (CD) spectroscopy, differential scanning densimetry (DSD), and dynamic light scattering measurements. At 25 degrees C, hexahistidine-tagged HodC has a hydrodynamic radius of 2.3 nm and is characterized by an unusually high degree of alpha-helical structure of approximately 60%, based on deconvolution of CD spectra. The percentage of beta-sheets and -turns is expected to be relatively low in view of its sequence similarity to proteins of the alpha/beta-hydrolase fold superfamily. His6HodC exhibits three-state unfolding (N <--> I <--> D) with an intermediate state I that exhibits at the transition temperature a volume larger than that of the native or denatured state. The intermediate state I is also associated with the highest isothermal expansion coefficient, alphaP, of the three states and exhibits a significantly lower percentage of alpha-helical structure than the native state. The stability difference between the native and intermediate state is rather small which makes I a potential candidate for reactions with various ligands, particularly those having a preference for the apparently preserved beta-type motifs.     

Intrinsically disordered protein from a pathogenic mesophile Mycobacterium tuberculosis adopts structured conformation at high temperature

Compared to eukaryotes, the occurrence of "intrinsically disordered" or "natively unfolded" proteins in prokaryotes has not been explored extensively. Here, we report the occurrence of an intrinsically disordered protein from the mesophilic human pathogen Mycobacterium tuberculosis. The Histidine-tagged recombinant Rv3221c biotin-binding protein is intrinsically disordered at ambient and physiological growth temperatures as revealed by circular dichroism and Fourier transform infrared (FTIR) spectroscopic studies. However, an increase in temperature induces a transition from disordered to structured state with a folding temperature of approximately 53 degrees C. Addition of a structure inducing solvent trifluoroethanol (TFE) causes the protein to fold at lower temperatures suggesting that TFE fosters hydrophobic interactions, which drives protein folding. Differential Scanning Calorimetry studies revealed that folding is endothermic and the transition from a disordered to structured state is continuous (higher-order), implying existence of intermediates during folding process. Secondary structure analysis revealed that the protein has propensity to form beta-sheets. This is in conformity with FTIR spectrum that showed an absorption peak at wave number of 1636 cm(-1), indicative of disordered beta-sheet conformation in the native state. These data suggest that although Rv3221c may be disordered under ambient or optimal growth temperature conditions, it has the potential to fold into ordered structure at high temperature driven by increased hydrophobic interactions. In contrast to the generally known behavior of other intrinsically disordered proteins folding at high temperature, Rv3221c does not appear to oligomerize or aggregate as revealed through numerous experiments including Congo red binding, Thioflavin T-binding, turbidity measurements, and examining molar ellipticity as a function of protein concentration. The amino acid composition of Rv3221c reveals that it has 24% charged and 54.9% hydrophobic amino acid residues. In this respect, this protein, although belonging to the class of intrinsically disordered proteins, has distinct features. The intrinsically disordered state and the biotin-binding feature of this protein suggest that it may participate in many biochemical processes requiring biotin as a cofactor and adopt suitable conformations upon binding other folded targets.     

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.     

Solvation of the folding-transition state in Pseudomonas aeruginosa azurin is modulated by metal: Solvation of azurin's folding nucleus

The role of water in protein folding, specifically its presence or not in the transition-state structure, is an unsolved question. There are two common classes of folding-transition states: diffuse transition states, in which almost all side chains have similar, rather low phi (phi) values, and polarized transition states, which instead display distinct substructures with very high phi-values. Apo-and zinc-forms of Pseudomonas aeruginosa azurin both fold in two-state equilibrium and kinetic reactions; while the apo-form exhibits a polarized transition state, the zinc form entails a diffuse, moving transition state. To examine the presence of water in these two types of folding-transition states, we probed the equilibrium and kinetic consequences of replacing core valines with isosteric threonines at six positions in azurin. In contrast to regular hydrophobic-to-alanine phi-value analysis, valine-to-threonine mutations do not disrupt the core packing but stabilize the unfolded state and can be used to assess the degree of solvation in the folding-transition state upon combination with regular phi-values. We find that the transition state for folding of apo-azurin appears completely dry, while that for zinc-azurin involves partially formed interactions that engage water molecules. This distinct difference between the apo-and holo-folding nuclei can be rationalized in terms of the shape of the free-energy barrier.     

Globin-like蛋白质折叠类型识别

蛋白质折叠类型识别是蛋白质结构研究的重要内容.以SCOP中的Globin-like折叠为研究对象,选择其中序列同一性小于25%的17个代表性蛋白质为训练集,采用机器和人工结合的办法进行结构比对,产生序列排比,经过训练得到了适合Globin-like折叠的概形隐马尔科夫模型(profile HMM)用于该折叠类型的识别.以Astrall.65中的68057个结构域样本进行检验,识别敏感度为99.64%,特异性100%.在折叠类型水平上,与Pfam和SUPERFAMILY单纯使用序列比对构建的HMM相比,所用模型由多于100个归为一个,仍然保持了很高的识别效果.结果表明:对序列相似度很低但具有相同折叠类型的蛋白质,可以通过引入结构比对的方法建立统一的HMM模型,实现高准确率的折叠类型识别.     

DNA toroids: stages in condensation.

The effects of polylysine (PLL) and PLL-asialoorosomucoid (AsOR) on DNA condensation have been analyzed by AFM. Different types of condensed DNA structures were observed, which show a sequence of conformational changes as circular plasmid DNA molecules condense progressively. The structures range from circular molecules with the length of the plasmid DNA to small toroids and short rods with approximately 1/6 to 1/8 the contour length of the uncondensed circular DNA. Single plasmid molecules of 6800 base pairs (bp) condense into single toroids of approximately 110 nm diameter, measured center-to-center. The results are consistent with a model for DNA condensation in which circular DNA molecules fold several times into progressively shorter rods. Structures intermediate between toroids and rods suggest that at least some toroids may form by the opening up of rods as proposed by Dunlap et al. [(1997) Nucleic Acids Res. 25, 3095]. Toroids and rods formed at lysine:nucleotide ratios of 5:1 and 6:1. This high lysine:nucleotide ratio is discussed in relation to entropic considerations and the overcharging of macroions. PLL-AsOR is much more effective than PLL alone for condensing DNA, because several PLL molecules are attached to a single AsOR molecule, resulting in an increased cation density.     

Ab initio folding of a trefoil‐fold motif reveals structural similarity with a β‐propeller blade motif

Many protein architectures exhibit evidence of internal rotational symmetry postulated to be the result of gene duplication/fusion events involving a primordial polypeptide motif. A common feature of such structures is a domain‐swapped arrangement at the interface of the N‐ and C‐termini motifs and postulated to provide cooperative interactions that promote folding and stability. De novo designed symmetric protein architectures have demonstrated an ability to accommodate circular permutation of the N‐ and C‐termini in the overall architecture; however, the folding requirement of the primordial motif is poorly understood, and tolerance to circular permutation is essentially unknown. The β‐trefoil protein fold is a threefold‐symmetric architecture where the repeating ~42‐mer “trefoil‐fold” motif assembles via a domain‐swapped arrangement. The trefoil‐fold structure in isolation exposes considerable hydrophobic area that is otherwise buried in the intact β‐trefoil trimeric assembly. The trefoil‐fold sequence is not predicted to adopt the trefoil‐fold architecture in ab initio folding studies; rather, the predicted fold is closely related to a compact “blade” motif from the β‐propeller architecture. Expression of a trefoil‐fold sequence and circular permutants shows that only the wild‐type N‐terminal motif definition yields an intact β‐trefoil trimeric assembly, while permutants yield monomers. The results elucidate the folding requirements of the primordial trefoil‐fold motif, and also suggest that this motif may sample a compact conformation that limits hydrophobic residue exposure, contains key trefoil‐fold structural features, but is more structurally homologous to a β‐propeller blade motif.     

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3?kcal?mol^?1. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.     

Soluble expression in Escherichia coli,purification and characterization of a human TF-1 cell apoptosis-related protein TFAR19

A novel human TF-1 cell apoptosis-related protein, TFAR19, cloned from a human leukemia cell line, TF-1, was first overexpressed in Escherichia coli with the sequence Met-Gly-His(6)-Gly-Thr-Asn-Gly, a hexahistidine sequence followed by a hydroxylamine cleavage site attached to its amino terminus. The resulting protein was soluble and single-step purified to homogeneity by metal chelating affinity chromatography. After cleavage of the purified His(6)-tagged TFAR19 sample with hydroxylamine, highly purified untagged TFAR19 protein was then obtained through an FPLC Resource Q column. The structural characteristics and function of the His(6)-tagged and untagged TFAR19 proteins were studied using circular dichroism, intrinsic fluorescence, and ANS-binding fluorescence spectra and apoptosis activity assay. The results show that alpha-helix is the main secondary structure of the proteins and the two forms of TFAR19 protein fold properly, which correspond well to their apoptosis activity expression. The results also indicate that the extra sequence including the His(6)-tag fused to the N-terminus of TFAR19 protein has a minimal effect on its structure and function, suggesting that the His(6)-tagged TFAR19 protein could be further used as an immobilized target for finding potential proteins which interact with TFAR19 from a cDNA library using in vitro ribosome display technique.     

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.     

Structural similarities between the N-terminal domain of Clostridium pasteurianum hydrogenase and plant-type ferredoxins

An N-terminal domain of Clostridium pasteurianum hydrogenase I, encompassing 76 residues out of the 574 composing the full-size enzyme, had previously been overproduced in Escherichia coli and shown to form a stable fold around a [2Fe-2S] cluster. This domain displays only marginal sequence similarity with [2Fe-2S] proteins of known structure, and therefore, two-dimensional 1H NMR has been implemented to elucidate features of the polypeptide fold. Despite the perturbing presence of the paramagnetic [2Fe-2S] cluster, 57 spin systems were detected in the TOCSY spectra, 52 of which were sequentially assigned through NOE connectivities. Several secondary structure elements were identified. The N terminus of the protein consists of two antiparallel beta strands followed by an alpha helix contacting both strands. Two additional antiparallel beta strands, one of them at the C terminus of the sequence, form a four-stranded beta sheet together with the two N-terminal strands. The proton resonances that can be attributed to this beta2alphabeta2 structural motif undergo no paramagnetic perturbations, suggesting that it is distant from the [2Fe-2S] cluster. In plant- and mammalian-type ferredoxins, a very similar structural pattern is found in the part of the protein farthest from the [2Fe-2S] cluster. This indicates that the N-terminal domain of C. pasteurianum hydrogenase folds in a manner very similar to those of plant- and mammalian-type ferredoxins over a significant part (ca. 50%) of its structure. Even in the vicinity of the metal site, where 1H NMR data are blurred by paramagnetic interactions, the N-terminal domains of hydrogenase and mammalian- and plant-type ferredoxins most likely display significant structural similarity, as inferred from local sequence alignments and from previously reported circular dichroism and resonance Raman spectra. These data afford structural information on a kind of [2Fe-2S] cluster-containing domain that occurs in a number of redox enzymes and complexes. In addition, together with previously published sequence alignments, they highlight the widespread distribution of the plant-type ferredoxin fold in bioenergetic systems encompassing anaerobic metabolism, photosynthesis, and aerobic respiratory chains.