Natural photoreceptors as a source of fluorescent proteins, biosensors, and optogenetic tools.

Genetically encoded optical tools have revolutionized modern biology by allowing detection and control of biological processes with exceptional spatiotemporal precision and sensitivity. Natural photoreceptors provide researchers with a vast source of molecular templates for engineering of fluorescent proteins, biosensors, and optogenetic tools. Here, we give a brief overview of natural photoreceptors and their mechanisms of action. We then discuss fluorescent proteins and biosensors developed from light-oxygen-voltage-sensing (LOV) domains and phytochromes, as well as their properties and applications. These fluorescent tools possess unique characteristics not achievable with green fluorescent protein-like probes, including near-infrared fluorescence, independence of oxygen, small size, and photosensitizer activity. We next provide an overview of available optogenetic tools of various origins, such as LOV and BLUF (blue-light-utilizing flavin adenine dinucleotide) domains, cryptochromes, and phytochromes, enabling control of versatile cellular processes. We analyze the principles of their function and practical requirements for use. We focus mainly on optical tools with demonstrated use beyond bacteria, with a specific emphasis on their applications in mammalian cells.

Unnatural Direct Interspecies Electron Transfer Enabled by Living Cell-Cell Click Chemistry.

Direct interspecies electron transfer (DIET) is essential for maintaining the function and stability of anaerobic microbial consortia. However, only limited natural DIET modes have been identified and DIET engineering remains highly challenging. In this study, an unnatural DIET between Shewanella oneidensis MR-1 (SO, electron donating partner) and Rhodopseudomonas palustris (RP, electron accepting partner) was artificially established by a facile living cell-cell click chemistry strategy. By introducing alkyne- or azide-modified monosaccharides onto the cell outer surface of the target species, precise covalent connections between different species in high proximity were realized through a fast click chemistry reaction. Remarkably, upon covalent connection, outer cell surface C-type cytochromes mediated DIET between SO and RP was achieved and identified, although this was never realized naturally. Moreover, this connection directly shifted the natural H2 mediated interspecies electron transfer (MIET) to DIET between SO and RP, which delivered superior interspecies electron exchange efficiency. Therefore, this work demonstrated a naturally unachievable DIET and an unprecedented MIET shift to DIET accomplished by cell-cell distance engineering, offering an efficient and versatile solution for DIET engineering, which extends our understanding of DIET and opens up new avenues for DIET exploration and applications.

Spectroscopic evidence for direct flavin-flavin contact in a bifurcating electron transfer flavoprotein.

A remarkable charge transfer (CT) band is described in the bifurcating electron transfer flavoprotein (Bf-ETF) from Rhodopseudomonas palustris (RpaETF). RpaETF contains two FADs that play contrasting roles in electron bifurcation. The Bf-FAD accepts electrons pairwise from NADH, directs one to a lower-reduction midpoint potential (E°) carrier, and the other to the higher-E° electron transfer FAD (ET-FAD). Previous work noted that a CT band at 726 nm formed when ET-FAD was reduced and Bf-FAD was oxidized, suggesting that both flavins participate. However, existing crystal structures place them too far apart to interact directly. We present biochemical experiments addressing this conundrum and elucidating the nature of this CT species. We observed that RpaETF missing either FAD lacked the 726 nm band. Site-directed mutagenesis near either FAD produced altered yields of the CT species, supporting involvement of both flavins. The residue substitutions did not alter the absorption maximum of the signal, ruling out contributions from residue orbitals. Instead, we propose that the residue identities modulate the population of a protein conformation that brings the ET-flavin and Bf-flavin into direct contact, explaining the 726 nm band based on a CT complex of reduced ET-FAD and oxidized Bf-FAD. This is corroborated by persistence of the 726 nm species during gentle protein denaturation and simple density functional theory calculations of flavin dimers. Although such a CT complex has been demonstrated for free flavins, this is the first observation of such, to our knowledge, in an enzyme. Thus, Bf-ETFs may optimize electron transfer efficiency by enabling direct flavin-flavin contact.

Immobilization of Rhodopseudomonas palustris P1 on glass pumice to improve the removal of NH4 + -N and NO2 - -N from aquaculture pond water.

We conducted this research in order to investigate the potential of a new material called glass pumice for use as a microorganism immobilization carrier to improve aquaculture pond water quality. The pH adjustment capacity and the Rhodopseudomonas palustris P1 cell adsorption capacity of glass pumice were measured. The immobilized Rps. palustris P1 and the free sample were compared to determine which had an enhanced NH4 + -N and NO2 - -N removal efficiency. The results showed that glass pumice significantly affected the pH of the acid solution (P < 0.05); the pH increased from 3.0 ± 0.08 to 7.21 ± 0.13 in 12 H. Rps. palustris P1 adsorption to glass pumice was rapid and reached equilibrium within 60 Min. The Langmuir adsorption parameter data showed that glass pumice had a higher affinity for Rps. palustris P1 than SiO2 powder, with an adsorption capacity of 4.02 × 108  cells g-1 . The maximum NH4 + -N and NO2 - -N removal rates by immobilized Rps. palustris P1 were 134.82 ± 0.67% and 93.68 ± 0.14% higher than those of nonimmobilized P1, respectively. Based on the above results, we propose that glass pumice is potential as a microorganism carrier material in aquaculture water treatment.

Structural Perturbations of Rhodopseudomonas palustris Form II RuBisCO Mutant Enzymes That Affect CO2 Fixation.

The enzyme ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) and its central role in capturing atmospheric CO2 via the Calvin-Benson-Bassham (CBB) cycle have been well-studied. Previously, a form II RuBisCO from Rhodopseudomonas palustris, a facultative anaerobic bacterium, was shown to assemble into a hexameric holoenzyme. Unlike previous studies with form II RuBisCO, the R. palustris enzyme could be crystallized in the presence of the transition state analogue 2-carboxyarabinitol 1,5-bisphosphate (CABP), greatly facilitating the structure-function studies reported here. Structural analysis of mutant enzymes with substitutions in form II-specific residues (Ile165 and Met331) and other conserved and semiconserved residues near the enzyme's active site identified subtle structural interactions that may account for functional differences between divergent RuBisCO enzymes. In addition, using a distantly related aerobic bacterial host, further selection of a suppressor mutant enzyme that overcomes negative enzymatic functions was accomplished. Structure-function analyses with negative and suppressor mutant RuBisCOs highlighted the importance of interactions involving different parts of the enzyme's quaternary structure that influenced partial reactions that constitute RuBisCO's carboxylation mechanism. In particular, structural perturbations in an intersubunit interface appear to affect CO2 addition but not the previous step in the enzymatic mechanism, i.e., the enolization of substrate ribulose 1,5-bisphosphate (RuBP). This was further substantiated by the ability of a subset of carboxylation negative mutants to support a previously described sulfur-salvage function, one that appears to rely solely on the enzyme's ability to catalyze the enolization of a substrate analogous to RuBP.

Structural basis of transcriptional regulation by CouR, a repressor of coumarate catabolism, in Rhodopseudomonas palustris.

The MarR family transcriptional regulator CouR, from the soil bacterium Rhodopseudomonas palustris CGA009, has recently been shown to negatively regulate a p-coumarate catabolic operon. Unlike most characterized MarR repressors that respond to small metabolites at concentrations in the millimolar range, repression by CouR is alleviated by the 800-Da ligand p-coumaroyl-CoA with high affinity and specificity. Here we report the crystal structures of ligand-free CouR as well as the complex with p-coumaroyl-CoA, each to 2.1-Å resolution, and the 2.85-Å resolution cocrystal structure of CouR bound to an oligonucleotide bearing the cognate DNA operator sequence. In combination with binding experiments that uncover specific residues important for ligand and DNA recognition, these structures provide glimpses of a MarR family repressor in all possible states, providing an understanding of the molecular basis of DNA binding and the conformation alterations that accompany ligand-induced dissociation for activation of the operon.

Ultrafast energy relaxation in single light-harvesting complexes.

Energy relaxation in light-harvesting complexes has been extensively studied by various ultrafast spectroscopic techniques, the fastest processes being in the sub-100-fs range. At the same time, much slower dynamics have been observed in individual complexes by single-molecule fluorescence spectroscopy (SMS). In this work, we use a pump-probe-type SMS technique to observe the ultrafast energy relaxation in single light-harvesting complexes LH2 of purple bacteria. After excitation at 800 nm, the measured relaxation time distribution of multiple complexes has a peak at 95 fs and is asymmetric, with a tail at slower relaxation times. When tuning the excitation wavelength, the distribution changes in both its shape and position. The observed behavior agrees with what is to be expected from the LH2 excited states structure. As we show by a Redfield theory calculation of the relaxation times, the distribution shape corresponds to the expected effect of Gaussian disorder of the pigment transition energies. By repeatedly measuring few individual complexes for minutes, we find that complexes sample the relaxation time distribution on a timescale of seconds. Furthermore, by comparing the distribution from a single long-lived complex with the whole ensemble, we demonstrate that, regarding the relaxation times, the ensemble can be considered ergodic. Our findings thus agree with the commonly used notion of an ensemble of identical LH2 complexes experiencing slow random fluctuations.

Postnovo: Postprocessing Enables Accurate and FDR-Controlled de Novo Peptide Sequencing.

De novo sequencing offers an alternative to database search methods for peptide identification from mass spectra. Since it does not rely on a predetermined database of expected or potential sequences in the sample, de novo sequencing is particularly appropriate for samples lacking a well-defined or comprehensive reference database. However, the low accuracy of many de novo sequence predictions has prevented the widespread use of the variety of sequencing tools currently available. Here, we present a new open-source tool, Postnovo, that postprocesses de novo sequence predictions to find high-accuracy results. Postnovo uses a predictive model to rescore and rerank candidate sequences in a manner akin to database search postprocessing tools such as Percolator. Postnovo leverages the output from multiple de novo sequencing tools in its own analyses, producing many times the length of amino acid sequence information (including both full- and partial-length peptide sequences) at an equivalent false discovery rate (FDR) compared to any individual tool. We present a methodology to reliably screen the sequence predictions to a desired FDR given the Postnovo sequence score. We validate Postnovo with multiple data sets and demonstrate its ability to identify proteins that are missed by database search even in samples with paired reference databases.

Apoprotein heterogeneity increases spectral disorder and a step-wise modification of the B850 fluorescence peak position.

It has already been established that the quaternary structure of the main light-harvesting complex (LH2) from the photosynthetic bacterium Rhodopseudomonas palustris is a nonameric 'ring' of PucAB heterodimers and under low-light culturing conditions an increased diversity of PucB synthesis occurs. In this work, single molecule fluorescence emission studies show that different classes of LH2 'rings' are present in "low-light" adapted cells and that an unknown chaperon process creates multiple sub-types of 'rings' with more conformational sub-states and configurations. This increase in spectral disorder significantly augments the cross-section for photon absorption and subsequent energy flow to the reaction centre trap when photon availability is a limiting factor. This work highlights yet another variant used by phototrophs to gather energy for cellular development.

Identification of protein W, the elusive sixth subunit of the Rhodopseudomonas palustris reaction center-light harvesting 1 core complex.

The X-ray crystal structure of the Rhodopseudomonas (Rps.) palustris reaction center-light harvesting 1 (RC-LH1) core complex revealed the presence of a sixth protein component, variably referred to in the literature as helix W, subunit W or protein W. The position of this protein prevents closure of the LH1 ring, possibly to allow diffusion of ubiquinone/ubiquinol between the RC and the cytochrome bc1 complex in analogous fashion to the well-studied PufX protein from Rhodobacter sphaeroides. The identity and function of helix W have remained unknown for over 13years; here we use a combination of biochemistry, mass spectrometry, molecular genetics and electron microscopy to identify this protein as RPA4402 in Rps. palustris CGA009. Protein W shares key conserved sequence features with PufX homologs, and although a deletion mutant was able to grow under photosynthetic conditions with no discernible phenotype, we show that a tagged version of protein W pulls down the RC-LH1 complex. Protein W is not encoded in the photosynthesis gene cluster and our data indicate that only approximately 10% of wild-type Rps. palustris core complexes contain this non-essential subunit; functional and evolutionary consequences of this observation are discussed. The ability to purify uniform RC-LH1 and RC-LH1-protein W preparations will also be beneficial for future structural studies of these bacterial core complexes.

Membrane-Bound Bacteriophytochrome-Like Complex of Phototrophic Purple Non-Sulfur Bacterium Rhodopseudomonas palustris.

A pigment-protein complex of yellow color with absorption maxima at 682 and 776 nm, characteristic for bacteriophytochromes, was isolated from the photosynthetic membranes of the purple bacterium Rhodopseudomonas palustris. Zinc-induced fluorescence of the complex indicated the presence of the biliverdin chromophore covalently bound to the protein. The parameters of low-temperature fluorescence (λ excitation at 680 nm, λ emission at 695 nm) indicated the ability of the complex to undergo photoconversion. These data, as well as the kinetics of accumulation of the red (Pr)-form on far red light, allowed the complex to be classified as a bacteriophytochrome-like complex with its localization in the photosynthetic membranes of Rps. palustris.

The NMR solution structure and function of RPA3313: a putative ribosomal transport protein from Rhodopseudomonas palustris.

Protein function elucidation often relies heavily on amino acid sequence analysis and other bioinformatics approaches. The reliance is extended to structure homology modeling for ligand docking and protein-protein interaction mapping. However, sequence analysis of RPA3313 exposes a large, unannotated class of hypothetical proteins mostly from the Rhizobiales order. In the absence of sequence and structure information, further functional elucidation of this class of proteins has been significantly hindered. A high quality NMR structure of RPA3313 reveals that the protein forms a novel split ßßαß fold with a conserved ligand binding pocket between the first ß-strand and the N-terminus of the α-helix. Conserved residue analysis and protein-protein interaction prediction analyses reveal multiple protein binding sites and conserved functional residues. Results of a mass spectrometry proteomic analysis strongly point toward interaction with the ribosome and its subunits. The combined structural and proteomic analyses suggest that RPA3313 by itself or in a larger complex may assist in the transportation of substrates to or from the ribosome for further processing. Proteins 2016; 85:93-102. © 2016 Wiley Periodicals, Inc.

The Lipid†A from Rhodopseudomonas palustris Strain BisA53 LPS Possesses a Unique Structure and Low Immunostimulant Properties.

The search for novel lipid A analogues from any biological source that can act as antagonists, displaying inhibitory activity towards the production of pro-inflammatory cytokines, or as immunomodulators in mammals, is a very topical issue. To this aim, the structure and immunological properties of the lipopolysaccharide lipid A from the purple nonsulfur bacterium Rhodopseudomonas palustris strain BisA53 have been determined. This lipid A displays a unique structural feature, with a non-phosphorylated skeleton made up of the tetrasaccharide Manp-α-(1→4)-GlcpN3N-ß-1→6-GlcpN3N-α-(1→1)-α-GalpA, and four primary amide-linked 14:0(3-OH) and, as secondary O-acyl substituents, a 16:0 and the very long-chain fatty acid 26:0(25-OAc), appended on the GlcpN3N units. This lipid A architecture is definitely rare, so far identified only in the genus Bradyrhizobium. Immunological tests on both murine bone-marrow-derived and human monocyte-derived macrophages revealed an extremely low immunostimulant capability of this LPS lipid A.

Mutation design of a thermophilic Rubisco based on three-dimensional structure enhances its activity at ambient temperature.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) plays a central role in carbon dioxide fixation on our planet. Rubisco from a hyperthermophilic archaeon Thermococcus kodakarensis (Tk-Rubisco) shows approximately twenty times the activity of spinach Rubisco at high temperature, but only one-eighth the activity at ambient temperature. We have tried to improve the activity of Tk-Rubisco at ambient temperature, and have successfully constructed several mutants which showed higher activities than the wild-type enzyme both in vitro and in vivo. Here, we designed new Tk-Rubisco mutants based on its three-dimensional structure and a sequence comparison of thermophilic and mesophilic plant Rubiscos. Four mutations were introduced to generate new mutants based on this strategy, and one of the four mutants, T289D, showed significantly improved activity compared to that of the wild-type enzyme. The crystal structure of the Tk-Rubisco T289D mutant suggested that the increase in activity was due to mechanisms distinct from those involved in the improvement in activity of Tk-Rubisco SP8, a mutant protein previously reported to show the highest activity at ambient temperature. Combining the mutations of T289D and SP8 successfully generated a mutant protein (SP8-T289D) with the highest activity to date both in vitro and in vivo. The improvement was particularly pronounced for the in vivo activity of SP8-T289D when introduced into the mesophilic, photosynthetic bacterium Rhodopseudomonas palustris, which resulted in a strain with nearly two-fold higher specific growth rates compared to that of a strain harboring the wild-type enzyme at ambient temperature. Proteins 2016; 84:1339-1346. © 2016 Wiley Periodicals, Inc.

DNA-directed spatial assembly of photosynthetic light-harvesting proteins.

This manuscript describes the surface immobilization of a light-harvesting complex to prescribed locations directed by the sequence-selective recognition of duplex DNA. An engineered light-harvesting complex (RC-LH1) derived from Rhodopseudomonas (Rps.) palustris containing the zinc finger (ZF) domain zif268 was prepared. The zif268 domain directed the binding of zfRC-LH1 to target double-stranded DNA sequences both in solution and when immobilized on lithographically defined micro-patterns. Excitation energy transfer from the carotenoids to the bacteriochlorophyll pigments within zfRC-LH1 confirmed that the functional and structural integrity of the complex is retained after surface immobilization.

Hydrogen production under salt stress conditions by a freshwater Rhodopseudomonas palustris strain.

Hydrogen represents a possible alternative energy carrier to face the growing request for energy and the shortage of fossil fuels. Photofermentation for the production of H2 constitutes a promising way for integrating the production of energy with waste treatments. Many wastes are characterized by high salinity, and polluted seawater can as well be considered as a substrate. Moreover, the application of seawater for bacterial culturing is considered cost-effective. The aims of this study were to assess the capability of the metabolically versatile freshwater Rhodopseudomonas palustris 42OL of producing hydrogen on salt-containing substrates and to investigate its salt stress response strategy, never described before. R. palustris 42OL was able to produce hydrogen in media containing up to 3 % added salt concentration and to grow in media containing up to 4.5 % salinity without the addition of exogenous osmoprotectants. While the hydrogen production performances in absence of sea salts were higher than in their presence, there was no significant difference in performances between 1 and 2 % of added sea salts. Nitrogenase expression levels indicated that the enzyme was not directly inhibited during salt stress, but a regulation of its expression may have occurred in response to salt concentration increase. During cell growth and hydrogen production in the presence of salts, trehalose was accumulated as a compatible solute; it protected the enzymatic functionality against salt stress, thus allowing hydrogen production. The possibility of producing hydrogen on salt-containing substrates widens the range of wastes that can be efficiently used in production processes.

Single-molecule spectroscopy reveals photosynthetic LH2 complexes switch between emissive states.

Photosynthetic organisms flourish under low light intensities by converting photoenergy to chemical energy with near unity quantum efficiency and under high light intensities by safely dissipating excess photoenergy and deleterious photoproducts. The molecular mechanisms balancing these two functions remain incompletely described. One critical barrier to characterizing the mechanisms responsible for these processes is that they occur within proteins whose excited-state properties vary drastically among individual proteins and even within a single protein over time. In ensemble measurements, these excited-state properties appear only as the average value. To overcome this averaging, we investigate the purple bacterial antenna protein light harvesting complex 2 (LH2) from Rhodopseudomonas acidophila at the single-protein level. We use a room-temperature, single-molecule technique, the anti-Brownian electrokinetic trap, to study LH2 in a solution-phase (nonperturbative) environment. By performing simultaneous measurements of fluorescence intensity, lifetime, and spectra of single LH2 complexes, we identify three distinct states and observe transitions occurring among them on a timescale of seconds. Our results reveal that LH2 complexes undergo photoactivated switching to a quenched state, likely by a conformational change, and thermally revert to the ground state. This is a previously unobserved, reversible quenching pathway, and is one mechanism through which photosynthetic organisms can adapt to changes in light intensities.

DEVELOPMENT OF SINGLE-DOMAIN NEAR-INFRARED FLUORESCENT PROTEIN GAF-FP BASED ON BACTERIAL PHYTOCHROME.

Fluorescent proteins (FPs) are widely used as genetically encoded markers for noninvasive and quantitative study of biological processes. Development of biomarkers that fluoresce in the near-infrared spectral range allows the study of animals at a deeper level due to high permeability of tissues to light in this wavelength range, compared to the visible light. For widespread use of FPs, such properties as low molecular weight and the monomer become important. In this paper, we developed a FP called the GAF-FP and based on the chromophore- binding domain of bacterial phytochrome from Rhodopseudomonas palustris (RpBphP1). GAF-FP has a molecular mass of ~ 19 kDa, 2 times lower than that of other FP based on BphPs and 1.4 times less than the commonly used GFP-like proteins. Unlike most other near-infrared FP, GAF-FP is a monomer, has high photostability and its structure can withstand the introduction of small peptide inserts. Moreover, GAF-FP can covalently bind two different tetrapyrrole chromophores: phycocyanobilin (PCB) and biliverdin (BV), which is found in mammalian tissues. GAF-FP with BV as a chromophore (GAF-FP­BV) has a main absorption band with a maximum at 635 nm and fluorescence maximum at 670 nm, whereby GAF-FP has a high signal to background ratio even if localized at a depth of several mm below the tissue surface. Apart from the near-infrared absorption band, GAF-FP­BV also has also an absorption band in the violet spectral range with a maximum at 378 nm. This property has been used by us to create a chimeric protein consisting of a modified luciferase from Renilla reniformis (RLuc8) and GAF-FP. We have shown that the chimeric protein is capable of resonance energy transfer from the substrate, which is oxidized by luciferase, to chromophore of GAF-FP­BV. In the absence of energy acceptor, RLuc8 catalyzes the cleavage of the substrate with light radiation having a peak of 400 nm. At the same time, as a part of GAF-FP­RLuc8 chimeric protein, the energy from the substrate is transferred to the chromophore of FP and then emitted in the near-infrared spectral range corresponding to GAF-FP fluorescence. These results open the way for the creation of new small near-infrared FPs based on various natural BphPs with a prospect of their wider use in cell and molecular biology.

Kinetically and Crystallographically Guided Mutations of a Benzoate CoA Ligase (BadA) Elucidate Mechanism and Expand Substrate Permissivity.

A benzoate CoA ligase (BadA), isolated from the bacterium Rhodopseudomonas palustris, catalyzes the conversion of benzoate to benzoyl CoA on the catabolic pathway of aromatic carboxylic acids. Herein, apparent Michaelis constants K(app)cat and K(app)M were determined for an expanded array of 31 substrates chosen to systematically probe the active site architecture of the enzyme and provide a baseline for expansion of wild-type substrate specificity. Acyl CoA products were observed for 25 of the 31 substrates; in general, BadA converted ortho-substituted substrates better than the corresponding meta and para regioisomers, and the turnover number was more affected by steric rather than electronic effects. The kinetic data are interpreted in relation to six crystal structures of BadA in complex with several substrates and a benzoyl-AMP reaction intermediate. In contrast to other known natural substrate-bound benzoate ligase structures, all substrate-bound BadA structures adopted the thiolation conformation instead of the adenylation conformation. We also observed all the aryl carboxylates to be uniquely oriented within the active site, relative to other structures. Together, the kinetics and structural data suggested a mechanism that involves substrate binding in the thiolation conformation, followed by substrate rotation to an active orientation upon the transition to the adenylation conformation. On the basis of this hypothesis and the structural data, sterically demanding active site residues were mutated, and the substrate specificity was expanded substantially versus that of BadA. Novel activities were seen for substrates with larger substituents, including phenyl acetate. Additionally, the mutant Lys427Ala identified this nonconserved residue as essential for the thiolation step of BadA, but not adenylation. These variously acylated CoAs can serve as novel substrates of acyl CoA-dependent acyltransferases in coupled enzyme assays to produce analogues of bioactive natural products.

Charge transfer through a cytochrome multiheme chain: theory and simulation.

We study sequential charge transfer within a chain of four heme cofactors located in the c-type cytochrome subunit of the photoreaction center of Rhodopseudomonas viridis from a theoretical perspective. Molecular dynamics simulations of the thermodynamic integration type are used to compute two key energies of Marcus' theory of charge transfer, the driving force ∆G and the reorganization energy λ. Due to the small exposure of the cofactors to the solvent and to charged amino acids, the outer sphere contribution to the reorganization energy almost vanishes. Interheme effective electronic couplings are estimated using ab initio wave functions and a well-parameterized semiempirical scheme for long-range interactions. From the resulting charge transfer rates, we conclude that at most the two heme molecules closest to the membrane participate in a fast recharging of the photoreaction center, whereas the remaining hemes are likely to have a different function, such as intermediate electron storage. Finally, we suggest means to verify or falsify this hypothesis.