Introduction

Biohydrogen production offers an appealing alternative to fossil fuel burning. It does not involve the production of toxic side products and is CO2-neutral. The different hydrogenases (Fe-, FeFe-, and NiFe-) taking part in biological hydrogen turnover evolved into extremely efficient catalysts since the primordial earth1,2,3. Many different options for biohydrogen evolution are available and are studied such as fermentative hydrogen production4, direct biophotolysis of water5,6, and exploitation of microbial consortia7. Direct biophotolysis is the most energy efficient process since it is converting sunlight and water directly to O2 and H2 without organic intermediates. However, it faces the challenge that the FeFe-hydrogenases, which are considered the most efficient H2-producing enzymes, are working specifically under anaerobic conditions and are inactivated by oxygen.

Cyanobacteria perform oxygenic photosynthesis and harbor two types of NiFe-hydrogenases8. Protein film voltammetry experiments showed that the cyanobacterial bidirectional NiFe-hydrogenase has a bias to hydrogen production9, which was considered uncommon for NiFe-enzymes. In addition, this hydrogenase is reduced by the low potential electron donors ferredoxin and flavodoxin10 supporting its role as H2-producer. Since this enzyme remains partially active in the presence of low concentrations of oxygen9 it is a promising candidate for biotechnological photohydrogen production. One of the major questions to be solved in this respect is if the hydrogen turnover of the cyanobacterial NiFe-hydrogenase is able to cope with the rate of photosynthetic electron transport. However, the rate of hydrogen formation at its active site is not known. We therefore developed an approach to first quantify the number of hydrogenase molecules in the cells. On this basis the turnover rates for hydrogen production of the enzyme could be determined.

The assembly of metalloenzymes such as hydrogenases requires a complex array of additional proteins11. Accessory proteins ensure that the correct metal is inserted into the active site. This is owing to the so-called Irving-Williams series of the affinity of the transition metals to whatever ligand. Predominantly due to their decreasing ionic radii and the increasing electrostatic effect the affinity increases in the order Mg2+  < Mn2+  < Fe2+  <Co2+  <Ni2+  <Cu2+ > Zn2+ 12. Cells are thus facing a problem. They need to avoid that a wrong metal is inserted, which might be bound with a higher affinity but does not confer catalytic activity. They tackle this challenge by expressing a number of metal-binding accessory proteins that channel the correct metal into the proper active site not allowing these metals to escape from these routes13,14. As soon as the respective metal-ion is part of the folded holoenzyme the chance to be exchanged becomes negligible at least inside the cell. In addition, metal ion homeostasis inside the cells is tightly controlled according to their binding affinities so that they are not able to supersede metals in inappropriate active sites. In brief, cells should contain the number of metal ions that correspond to the number of appropriate binding sites including enzymes and storage proteins.

Ni-enzymes are an intriguing example of these mechanisms and recent years saw a considerable increase in our understanding of how they acquire nickel11,15. Especially, active NiFe-hydrogenases are the product of the sophisticated interaction of a minimum of 8 accessory proteins and a nickel transporter16,17.

In this study we used the model cyanobacterium Synechocystis sp. PCC 6803. It harbors a single NiFe-hydrogenase that produces hydrogen under fermentative conditions as well as photohydrogen when cells are shifted from anaerobic conditions into the light10,18,19,20. This NiFe-hydrogenase and the urease are the only known nickel enzymes in Synechocystis21.

Since the nickel transporter necessary for the expression of an active hydrogenase remained elusive we deleted the putative transporter nikKLMQOnik) to unravel its specificity. ICP-MS was used to determine the number of nickel ions per cell in wild type cells, the hydrogenase deletion strain (ΔhoxH) and the hydrogenase overexpression strain (oehox). On this basis we were able to determine for the first time the number of hydrogenase molecules per cell and to deduce the in-vivo turnover frequency of the enzyme in different strains. This allowed us to compare the efficiency of cyanobacterial in-vivo photohydrogen production with optimized in-vitro systems of FeFe- and NiFe-hydrogenases from heterotrophic organisms that were tethered to PSI22,23. Knowledge about the enzymatic characteristics of the cyanobacterial NiFe-hydrogenase is basic to evaluate the feasibility of biotechnological approaches with this enzyme.

Results

Characterization of Δnik mutant

The putative nickel transporter NikKLMQO (sll0381–sll0385) was deleted from the genome of Synechocystis and complete segregation was checked by PCR and Southern blotting (s. Fig. S1). Wild type cells and Δnik were grown in media with and without NTA and with supplements of Co2+ or Ni2+ (Fig. 1). NTA is a very efficient chelator of nickel and cobalt and therefore lowers the concentration of the free metal ions in the growth medium. Consequently, a high affinity transporter is necessary for their uptake. If NTA is added in combination with high concentrations (300 µM) of nickel or cobalt it masks these metals and it is thus possible to increase their concentration above otherwise toxic levels. This large supply enables the cells to unspecifically take up these metals via other metal transporters.

Figure 1
figure 1

Growth curves of wild type cells and Δnik in BG-11 and BG-11 with 1 mM NTA alone or 1 mM NTA and either 300 µM Co2+ or 300 µM Ni2+. Growth of Δnik can be rescued upon addition of NTA and 300 µM Ni2+, which allows for unspecific uptake of nickel. In each case three independent cultures were measured in triplicates. The standard deviations are indicated and are smaller than the data points if not visible. Please note that the curves of Δnik with 1 mM NTA and Δnik with 1 mM and 300 µM Co2+ are on top of each other.

The Δnik mutant grows significantly slower than wild type cells under normal conditions. If the availability of nickel is further decreased by addition of NTA (+1 mM NTA and +1 mM NTA + 300 µM Co2+) Δnik cells loose their ability to grow. However, the growth phenotype of Δnik can be rescued by addition of 1 mM NTA and 300 µM Ni2+, which facilitates unspecific nickel uptake. Addition of 300 µM Co2+ did not have any effect. This strongly suggests that nikKLMQO (sll0381, sll0382, sll0383, sll0384, and sll0385) encodes a nickel uptake system in Synechocystis.

Hydrogenase and urease activity

Synechocystis possesses two known nickel-containing enzymes: a bidirectional NiFe-hydrogenase and a urease. In order to further characterize the Δnik mutant and to investigate the nickel dependency of both enzymes, samples were taken from wild type cells, ΔhoxH, oehox, and Δnik grown under normal conditions and subjected to measurements of hydrogenase and urease activity (Fig. 2) Hydrogenase activity was determined by adding the artificial electron donor methylviologen, which allows the quantification of functional hydrogenase in these cells. Δnik shows only about 1/10th of the hydrogenase activity of the wild type cells (Fig. 2A) and its urease activity is not detectable under the assay conditions used (Fig. 2B). These results further support the assumption that nikKLMQO codes for a nickel transporter.

Figure 2
figure 2

Hydrogenase (A) and urease (B) activities of the different strains used in this study. The activities of three independent cultures were measured in triplicates. The standard deviation is indicated.

Metal content of different Synechocystis strains

In order to quantify the metal content of wild type cells, ΔhoxH, oehox, and Δnik their cells were harvested, washed with EDTA, oxidatively broken down and subjected to ICP-MS. Table 1 shows the number of metal atoms per cell. The obtained values were in good agreement with published data as described in the following in detail.

Table 1 The number of metal atoms detected in the cells of the different strains used in this study.

Of the metals measured magnesium is the most abundant. Since every chlorophyll molecule contains one magnesium atom we plotted the number of magnesium atoms per cell versus the number of chlorophyll molecules (s. Fig. S2). This yields a linear relationship with a y-axis intercept at about 4.36 × 107 atoms per cell. This is the number of magnesium atoms not bound to chlorophyll and corresponds in the wild type with a cell size of 3.1 µm3 to a concentration of 21 mM. This is well in the range of magnesium concentration of 15 to 25 mM found in cells24.

Most of the phosphorus in the cells resides in RNA that makes up about 20% of their dry weight25. There are 3.5 × 108 P-atoms in a wild type cell (Table 1). 20% of their dry weight would correspond to about 3.3 × 108 P-atoms in RNA leaving sufficient P-atoms for the DNA and the different nucleotides (ATP, GTP, ADP, GDP etc.).

The sulfur content of a cell gives a direct measure of its protein content since 5% of all amino acids are methionine and cysteine in an average proteome26. A plot of the S-atoms per cell versus the P-atoms per cell shows a linear relationship (s. Fig. S3) suggesting a constant ratio of protein to RNA for the different strains as also described by the Redfield ratio for phytoplankton27.

Previous investigations found 4.5 × 105 Mn-atoms per cell in Mn-limited cultures of Synechocystis28 and 4.2 × 106 Fe-atoms in normal cultures29. In our case these values are in a similar range with 3.6 × 105 Mn and 6.2 × 106 Fe. Probably small variations in the preparation of the medium and the washing procedure of the cells cause these deviations.

It is obvious from the table that the cells of the Δnik mutant have a considerably lower nickel content (4,000 per cell) compared to the other strains (15,000 in wild type cells, 17,000 in oehox and 13,000 in ΔhoxH) as expected. All other metal contents are higher in the Δnik indicating an increased uptake possibly triggered by the lack of nickel.

From inspection of the cells at the light microscope and measurements at the Coulter Counter it was obvious that the cells of the ΔhoxH mutant had a larger volume (7.7 µm3 compared to 3.1 µm3 in the wild type cells). On the basis of a volume ratio of 2.5 the data obtained from the ICP-MS was corrected in ΔhoxH.

If comparing the corrected metal content per cell of the ΔhoxH (Table 1) it is obvious that most transition metal contents are close to those of the other strains except cobalt, which is significantly higher and which might indicate a higher uptake activity caused by the lack of the hydrogenase. Taken together the ICP-MS measurements yielded convincing numbers of metal atoms per cell that were all well in range with previously published data.

Photosystem stoichiometry and number of photosystems per cell

As we were interested to determine the ratio between hydrogenase molecules to PSI and PSII, the number photosystems per cell was determined. The chlorophyll content of the different strains (Table 2) and the 77 K chlorophyll fluorescence spectra (Fig. 3) were used to estimate the total number of photosystems per cell. According to established protocols a ratio of the photosystems can be approximated by deconvoluting the areas of their fluorescence maxima in the 77 K spectra30. Using this method we determined the ratio of PSI/PSII in wild type cells to 5.6 in the oehox to 3.6 and in the ΔhoxH as 3.8.

Table 2 Number of chlorophyll molecules and calculated PSII dimers and PSI monomers per cell.
Figure 3
figure 3

77 K spectra of wild type cells, oehox, and ΔhoxH. The peak at 720 nm originates from PSI and the two peaks at 685 and 695 nm correspond to CP47 and the reaction center of PSII, respectively. All the curves were normalized to the same fluorescence emission at 720 nm.

The number of chlorophyll molecules per photosystem II reaction center is 3531 and 96 per photosystem I reaction center32. About 25% of the chlorophyll is bound to small CABs (chlorophyll a binding protein) in Synechocystis33. Based on this data it is possible to estimate the number of photosystems according to the formulas given in the materials and methods section. Wild type cells contain 23,600 PSII dimers (47,200 PSII reaction centers) and 268,000 PSI reaction centers (Table 2). This corresponds to 188,800 Mn atoms and 3.21 × 106 Fe atoms. In all strains the metal atoms of the reaction centers amount to about half of the respective metal measured in the different strains (compare Tables 1 and 2).

Number of NiFe-hydrogenase molecules per cell

The hydrogenase activity assay used with 5 mM methylviologen and 10 mM dithionite allows activity measurements with whole cells. As shown in Fig. 4 the activity is linearly dependent on the number of cells added into the assay mixture over a wide range of cell numbers and could not be saturated even with about 6 × 109 cells/ml. This correlation convincingly shows that the assay directly mirrors the amount of hydrogenase enzymes in the cells.

Figure 4
figure 4

Hydrogenase acitivity as measured in the presence of 5 mM methylviologen and 10 mM dithionite of different amounts of cells. Cells from three independent wild type cultures were measured.

When nickel content and hydrogenase activity per cell of wild type cells, oehox, and ΔhoxH were plotted against each other a linear relationship was obtained (Fig. 5). Since the correct insertion of the different metals into the different metalloenzymes requires a tight control of their content inside the cells13,14 we can thus assume that the different strains contain just the number of nickel atoms necessary for the number of active sites of their hydrogenase molecules besides a constant amount of other nickel dependent enzymes (e.g. urease) and nickel storage proteins (e.g. HypB1). It was possible to quantify the nickel pool that is not bound to the NiFe-hydrogenase by including ΔhoxH (corrected for its larger cell size) in our measurements. These Ni-atoms are not bound to hydrogenase molecules but to all other Ni-dependent enzymes. Since there is a single Ni-atom bound in the active site of the NiFe-hydrogenase we used this correlation to count the number of NiFe-hydrogenases in the cells according to their nickel content (corrected by the number of nickel atoms that are also present in ΔhoxH and therefore not bound to the NiFe-hydrogenase).

Figure 5
figure 5

The number of Ni-atoms per cell plotted against the hydrogenase activity of ΔhoxH (open squares), wild type cells (filled circles) and the hydrogenase overexpression strain (oehox)(red triangles). Of each strain four independent cultures were measured at least in triplicates. Standard deviations are indicated.

Ni-content and hydrogenase measurements correlated in a linear manner with a good statistical support with a correlation coefficient of 0.87. Using the derived equation we could therefore calculate that there are 1,100 ± 94 hydrogenase molecules per wild type cell and 3,790 ± 485 molecules in the cells of the overexpression strain.

To determine the variability of the number of hydrogenase molecules per cell and its ratio to the photosynthetic reaction centers we measured wild type cells from different culturing conditions and cell densities. Thereby, we obtained values of 400 to 2000 hydrogenase molecules per cell and ratios of one hydrogenase per 50 to 240 PSI centers and 40 to 55 PSII centers.

Turnover frequency of the NiFe-hydrogenase

Based on the knowledge of the precise number of hydrogenase molecules per cell it is possible to determine the H2 turnover frequency per active site. Under the assay conditions with methylviologen as artificial electron donor and dithionite as reducing agent this turnover frequency is 1,340 ± 115 H2/s. This was calculated on the basis of the evolved H2 (Fig. 4) and the number of hydrogenase molecules in the wild type cells. Photohydrogen production is measured under physiological conditions without artificial electron donors but instead by illuminating anaerobic, dark-adapted cells. We determined turnover frequencies of NiFe-hydrogenases of different cyanobacterial strains when producing photohydrogen under in–vivo conditions using previously published data20. We first quantified the number of NiFe-hydrogenases on the basis of the available methylviologen and dithionite dependent measurements by means of the equation shown in Fig. 5 and then calculated the turnover frequencies based on the measured rate of photohydrogen production.

These calculations represent turnover frequencies of the NiFe-hydrogenase under physiological in-vivo conditions. Numbers are ranging from 18 H2/s for wild type cells to about 120 H2/s in case of the ΔctaIΔcyd mutant in normal BG-11 medium with nitrate (Table 3). The ΔctaIΔcyd mutant lacks the two respiratory terminal oxidases, the cytochrome c oxidase (ctaI) and the cytochrome bd oxidase (cyd) localized in the thylakoid membrane. These oxidases are able to withdraw part of the electrons from the light reaction even under anaerobic conditions due to the oxygen produced by PSII activity. Consequently the deletion mutant shows a higher photohydrogen production rate at the onset of light20.

Table 3 Turnover frequency for hydrogen production in the light in different strains.

In wild type cells the turnover frequency of the hydrogenase was dependent on the culturing conditions. The highest values of photohydrogen production of wild type cells were found at low cell densities (which corresponds to high photosynthetic activity) if the cells were shifted to nitrate free medium and directly measured. Under these conditions a maximum turnover frequency of the wild type cells of 62 H2/s was found.

For comparison the M55 mutant (45 H2/s), the mutant lacking all terminal respiratory oxidases (ΔctaIΔcydΔctaII, 75 H2/s) and another double mutant (ΔcydΔctaII, 69 H2/s) are also shown in Table 3. Since the M55 has no NDH-1 its cyclic electron transfer is very slow. It has a low hydrogenase activity and produces photohydrogen for several minutes19,20.

Discussion

Cyanobacteria are known to harbor a number of different nickel dependent enzymes like the NiFe-hydrogenases, urease and SOD. The expression of these enzymes requires high affinity nickel uptake. Here we show on the basis of growth curves (Fig. 1) and activity measurements of hydrogenase and urease (Fig. 2) that the sole high affinity nickel transporter in Synechocystis belongs to a new type of energy coupling factor transporters and is encoded by nikKLMQO (sll0381 to sll0385). This result confirms predictions based on bioinformatics approaches34. In contrast to previous assumptions HupE is not able to complement the lack of this transporter35. Determination of nickel by ICP-MS of the hydrogenase deletion strain, wild type cells and a hydrogenase overexpression strain allowed the correlation of their nickel content and hydrogenase activity (Fig. 5). On the basis of this data we calculated about 1,100 hydrogenase molecules per wild type cell. Since we found 268,000 photosystem I and 47,200 photosystem II centers per cell on the basis of their chlorophyll content and 77 K spectra there is about 1 hydrogenase molecule per 240 PSI and 43 PSII reaction centers.

The high number of about 13,200 nickel atoms per cell not inserted into the active site of the hydrogenase (Fig. 5) suggests that the number of nickel dependent enzymes in Synechocystis might be higher than predicted on the basis of its genome sequence. Only the bidirectional hydrogenase, the urease and some other nickel binding proteins like the hydrogenase accessory genes HypB1, HypB2 and the urease accessory genes are known until now. Therefore, the actual metalloproteome of Synechocystis is probably underestimated as also shown for Pyrococcus furiosus36.

From this point on it is possible to also calculate the turnover frequencies of this NiFe-hydrogenase under different conditions and in different strains. In the presence of the artificial electron donor methylviologen we found the highest kcat of about 1340 H2/s. This is in a similar range as the published value of the purified soluble hydrogenase (SH) of Ralstonia eutropha. The SH is closely related to the cyanobacterial bidirectional enzyme and has a kcat of 372 H2/s when reducing ferricyanide with hydrogen37.

Highest rates of hydrogen turnover in Synechcoystis are achieved during photohydrogen evolution. The turnover frequency of wild type cells of up to 62 H2/s in the absence of nitrate still seems to be low since an oxidase mutant (ΔctaIΔcyd) even in the presence of nitrate shows about 120 H2/s (Table 3), which is up to now the strain with the highest rate achieved19,20. These results show that the hydrogen turnover is mostly dependent on the electron supply to the enzyme in-vivo and not limited by its characteristics.

The turnover frequencies of a variety of NiFe-hydrogenases has been measured. In this respect especially the SH of R. eutropha is interesting. For this enzyme a kcat of 342 H2/s has been determined for the reduction of NAD+ with H238. However, it has been shown by protein film voltammetry that the lower limit of the NiFe-hydrogenase of Allochromatium vinosum is 6,000 H2/s for the oxidation of H239. These latter results suggest that the catalytic turnover of NiFe-hydrogenases can be significantly enhanced when it is possible to accelerate electron transfer rates to the enzyme.

It is very interesting to compare the turnover frequencies obtained in this study with those obtained by others using in-vitro constructs for light driven hydrogen production. Especially FeFe-hydrogenases are considered as highly efficient catalysts. Turnover frequencies of these enzymes of up to 21,000 ± 12,000 H2/s have been measured40. A FeFe-hydrogenase of Clostridium acetobutylicum tethered to PSI via a sulfhydryl linker showed a turnover rate of 105 ± 22 e/PSI/s22 corresponding to about 53 H2/s. In another study it was shown that the membrane bound hydrogenase (MBH) of R. eutropha fused to the PsaE subunit of photosystem I yielded a rate of 75 ± 19 H2/s per complex when assembled on a gold electrode23. Thus, the rates of about 62 H2/s in wild type cells and 120 H2/s of the oxidase mutant measured here are comparable to and even higher than these systems working under optimal conditions.

The measured in-vivo turnover frequency of the bidirectional NiFe-hydrogenase of 120 H2/s in the oxidase mutant (ΔctaIΔcyd) is also above the in-vivo photosynthetic electron transfer rate of 47 e/s estimated for PSI22 (equivalent to 24 H2/s) and even the rates measured in-vitro of 230 e/s for isolated PSI (equivalent to 115 H2/s) saturated with ascorbate and ferredoxin22. Thus, the catalytic efficiency of the cyanobacterial NiFe-hydrogenase should not limit photohydrogen production when redirecting the cellular electron flow to this enzyme.

Our investigations expand previous electrochemical investigations9 and show that the bidirectional NiFe-hydrogenase of Synechocystis is a highly efficient catalyst for hydrogen evolution in-vivo. They show that metabolic engineering especially to direct photosynthetic electron transfer to this hydrogenase holds a great potential for high rates of biotechnological H2 production. The catalytic rates found here outcompete in-vitro systems poised at optimal redox potentials and working in the presence of plenty of electron donors. Since the cyanobacterial bidirectional NiFe-hydrogenase does not suffer from irreversible inhibition by oxygen and is partially active under low oxygen levels9 it holds an additional advantage compared to FeFe-hydrogenases. Moreover, in-vivo production systems would benefit from the self-repair inherent to living organisms.

Material and Methods

Growth conditions

Synechocystis sp. PCC 6803 wild type cells and the mutant strains were grown in BG-1141 supplemented with 5 mM TES pH 8. This medium contains a trace metal mix including Mn2+, Cu2+, Co2+, Zn2+ and MoO4. The iron concentration in the medium is 22.9 µM. The mutant strains were segregated on BG-11 plates containing 50 µg/ml kanamycin. For physiological measurements cultures were grown in 200 ml BG-11 medium bubbled with air at a light intensity of 50 µE/m2/s. To deplete the cells of Co2+ and Ni2+ 1 mM nitrilotriacetic acid was added in some experiments as previously described35. At an OD750 between 3 and 4 the cells were harvested and subjected to all further measurements.

Construction of mutants

To create a deletion mutant (Δnik) of the putative nickel transporter genes PCR products containing the up- and downstream regions of the nikKLMQO (sll0381, sll0382, sll0383, sll0384 and sll0385) gene cluster and the kanamycin resistance cassette were prepared by PCR fusion as described previously35. The following primers were used P1nik CCGTGGGCAAAATCTACCCT and P2nik CTTTCTGGCTGGATGATGGGGCGATTAAAGGCGATCAGCAAACTGTGGG for the upstream region and P3nik TGTTGGAATTTAATCGCGGCCTCGAACCCACTTTTAAATGGCGAATCGC and P4nik ATCGTCAGGATCGTCTGCGG for the downstream region. The resulting PCR product was ligated into the pGEM-T vector (Promega, Madison, USA), checked by sequencing and transformed into Synechocystis sp, PCC 6803 wild type cells using established protocols42.

The construction of the deletion strain of the gene of the large hydrogenase subunit ΔhoxH was described previously18.

The strain overproducing the NiFe-hydrogenase was constructed by inserting the psbAII-promoter upstream of the hox-operon and the hypFCDEAB-operon form Nostoc sp. PCC 7120 under the control of the same promoter as described43. This strain retains the original structure of the hox-operon.

Hydrogenase measurements

Hydrogenase was determined directly in whole cells as previously described using a Clark-type electrode pretreated for H2 measurements20. The electrode was connected to a lab-made control box maintaining a potential of −600 mV. An aliquot of the cells corresponding to 5 µg chlorophyll was given into the electrode chamber that was kept at a constant temperature of 28 °C. The chamber contained a total volume of 1 ml with the cells, 5 mM methylviologen and 10 mM sodium dithionite.

Urease measurements

To measure the urease activity cells were pelleted by centrifugation and resuspended in a buffer with 100 mM Tris pH 8 and 150 mM NaCl. The cell suspension was vigorously vortexed with glass beads (diameter 0.17–0.18 mm) 3 × 3 min at 6 °C with intermittent cooling on ice. After centrifuging the suspension at 800 × g for 1 min at 4 °C the liquid phase was removed and centrifuged at 1,300 × g for 10 min to remove residual glass beads. An aliquot of this cell homogenate was used to determine the urease activity according to Kaltwasser and Schlegel44. The assay mixture contained in 1 ml 0.8 mM α-ketoglutarate, 0.03 M urea, 0.25 mM NADH and 9 U glutamate dehydrogenase and was buffered with 0.04 M Tris pH 8.0. The activity was measured by following the consumption of NADH at 366 nm with a spectrophotometer (Shimadzu UV-2501PC, Kyoto, Japan) at 20 °C.

Determination of the cell number

The cell number and diameter was determined at least in triplicates with a Multisizer 3 (Beckman Coulter, Krefeld, Germany) with a capillary opening of 20 µm.

Chlorophyll determination

A defined volume of cell suspension was pelleted by centrifugation for 5 min at 15,000 × g. The supernatant was removed. To extract the chlorophyll the cell pellet was vigorously vortexed with 1 ml methanol for 5 min. After an additional centrifugation at 15,000 × g for 10 min the supernatant was measured at a spectrophotometer at the wavelengths 665, 665.5, 666, 666.5 and 750 nm. The amount of chlorophyll was determined by taking the highest absorbance between 665 and 666.5 nm (Abshigh) according to the formula [Chl a] = (Abshigh − Abs750)/0.0809 µg/ml45.

77 K fluorescence measurements

Chlorophyll fluorescence spectra at 77 K were recorded with a Hitachi F-4500 fluorescence spectrometer (Hitachi, Tokyo, Japan). The samples were diluted to 2 µg chlorophyll/ml frozen in liquid nitrogen and excited at 440 nm. The spectrum was recorded between 660 and 800 nm. By deconvolution of the spectrum it is possible to determine the ratio of PSI/PSII by using the areas of the different component bands30.

Estimation of photosystems per cell

To estimate the number of photosystems per cell we used the total amount of chlorophyll measured in the cells, the cell number, the ratio of the two photosystems as determined from the 77 K spectra and the number of chlorophyll molecules per reaction center as given from the most recent X-ray structures31,32 according to the following formulas:

$$Number\,of\,PSI\,centers\,per\,cell=\frac{number\,of\,chlorophylls\,\times \,0.75}{cell\,number}\,\times \,\frac{\frac{PSI}{PSII}}{\frac{PSI}{PSII}\,\times \,96+35}$$
(1)
$$Number\,of\,PSII\,centers\,per\,cell=\frac{number\,of\,chlorophylls\,\times \,0.75}{cell\,number}\,\times \,\frac{1}{\frac{PSI}{PSII}\,\times \,96+35}$$
(2)

PSI/PSII is the ratio of the photosystems as determined from the 77 K spectra, where PSII was set to 1. The total amount of chlorophyll was reduced by 25% since this proportion of chlorophyll is bound to CABs (chlorophyll a binding proteins)33.

Treatment of the cells for ICP-MS measurements

To analyze the different strains for their metal content their cell pellets of 150 ml culture were washed twice in 3 ml 0.2 M ultrapure EDTA pH 8. The pellet was resuspended in 3 ml distilled water and 700 µl of this suspension were pelleted again. This pellet was dried at 65 °C overnight and dissolved in 1 ml 69% aristar grade HNO3 and 500 µl ultrapure 30% H2O2. After three days in daylight the cell constituents were completely broken down by oxidation reactions. All steps and the harvesting were done in Greiner Bio one polypropylene tubes. An equal volume of BG-11 medium treated the same way was used as a control. For all cultures at least three replicates were prepared.

ICP-MS

Samples were analyzed by inductively coupled-mass-spectrometry (ICP-MS) using a high resolution mass spectrometer (Thermo Element XR). All selected isotopes (25Mg, 31P, 32S, 44Ca, 55Mn, 56Fe, 59Co, 60Ni, 63Cu, 66Zn and 115In) were measured in medium ass resolution (4,000 RP) and carefully monitored for interferences. Prior to analysis samples were diluted 100fold and spiked with 2.5 ng/l In for internal standardization. Analytical results are blank-subtracted averages of three runs. Accuracy of the results was monitored by analyzing CRM’s NIST 1643e, NIST 1640 a and LGC6019 as unknowns. Average precision of results was estimated from replicate analyses of 12 samples and found to be better than <0.5% RSD for all elements except Ni with 2.5% RSD. Additional control measurements also including 98Mo were performed using quadrupole-based ICP-MS (Agilent 7500 cs) with lower sample dilution of 50fold, now giving a reproducibility of 0.8 and 0.3% RSD for Ni and Mo, respectively. A spiking experiment with 3 samples resulted in recoveries of >95% for all analytes.

Data availability

All data generated or analysed during this study are included in this published article (and its Supplementary Information files).