Systems analysis of the CO2 concentrating mechanism in cyanobacteria

  1. Niall M Mangan  Is a corresponding author
  2. Michael P Brenner  Is a corresponding author
  1. Harvard University, United States
6 figures, 4 tables and 1 additional file

Figures

Schematic of the CCM in cyanobacteria.

Outer and cell membranes (in black), as well as, thylakoid membranes where the light reactions take place (in green) are treated together. Carboxysomes are shown as four hexagons evenly spaced along the centerline of the cell. The model treats a spherically symmetric cell, with one carboxysome at the center. Active HCO3 transport into the cell is indicated (in light blue), as well as active conversion from CO2 to HCO3, sometimes called ‘facilitated uptake’ or ‘scavenging’, at membranes (in orange). Both CO2 and HCO3 can leak in and out of the cell, with CO2 leaking out much more readily. Both species passively diffuse across the carboxysome shell. Carbonic anhydrase (red) and RuBisCO (blue) are contained in the carboxysomes and facilitate reactions as shown.

https://doi.org/10.7554/eLife.02043.003
Figure 2 with 2 supplements
Phase space for HCO3 transport, jc, and carboxysome permeability kc.

Plotted are the parameter values at which the CO2 concentration reaches some critical value. The left most line (dark blue) indicates for what values of jc and kc the CO2 concentration in the carboxysome would half-saturate RuBisCO (Km). The middle line (light blue) indicates the parameter values which would result in a CO2 concentration where 99% of all RuBisCO reactions are carboxylation reactions and only 1% are oxygenation reactions when O2 concentration is 260 μM. The right most (red) line indicates the parameter values which result in carbonic anyhdrase saturating. Here α = 0, so there is no CO2 scavenging or facilitated uptake. The dotted line (grey) shows the kc and jc values, where the HCO3 concentration in the cytosol is 30 mM. The HCO3 concentration in the cytosol does not vary appreciably with kc in this parameter regime, and reaches 30 mM at jc0.6cms. All other parameters, such as reaction rates are held fixed and the value can be found in the Table 1 and Table 2.

https://doi.org/10.7554/eLife.02043.007
Figure 2—figure supplement 1
Phase space for HCO3 transport and carboxysome permeability.

Solid lines show lines of constant CO2 concentration in the carboxysome for Dc=1e5cm2s, or the diffusion constant of small molecule in water. Dashed lines show the same lines of constant CO2 concentration, bur for Dc=1e7cm2s, or the diffusion constant of a small molecule in a 60% sucrose solution.

https://doi.org/10.7554/eLife.02043.008
Figure 2—figure supplement 2
Phase space for HCO3 transport, jc, and carboxysome permeability, kc.

Plotted are the parameter values at which CO2 concentration reaches some critical value. The left most line (dark blue) indicates for what values of jc and kc the CO2 concentration in the carboxysome would saturate RuBisCO. The middle line (light blue) indicates the parameter values which would result in a CO2 concentration where 99% of all RuBisCO reactions are carboxylation reactions and only 1% are oxygenation reactions when O2 concentration is 260 µM. The right most (red) line indicates the parameter values which result in carbonic anyhdrase saturating. Here α=0cms (solid lines) and α=0cms (dashed line), showing the effect of CO2 scavenging or facilitated uptake on the phase space. All other parameters, such as reaction rates are held fixed and the value can be found Table 1.

https://doi.org/10.7554/eLife.02043.009
Figure 3 with 1 supplement
Numerical solution (diamonds and circles) and analytic solutions (carbonic anhydrase unsaturated, solid lines, and saturated, dashed lines) correspond well.

HCO3 transport is varied, and all other system parameters are held constant. The CO2 concentration above which RuBisCO is saturated is Km (grey dashed line). The CO2 concentration where the oxygen reaction error rate will be 1% is C99% (grey dash-dotted line). The transition between carbonic anyhdrase being unsatruated and saturated happens where the two analytic solutions meet (where the dashed and solid red lines meet). Below a critical value of transport, jc103cms the level of transport is lower than the HCO3 leaking through the cell membrane. A value of kc=103 cm/s for the carboxysome permeability was used for these calculations.

https://doi.org/10.7554/eLife.02043.010
Figure 3—figure supplement 1
No effect of localizing carbonic anhydrase to the shell of the carboxysome.

We assume the same amount of carbonic anhydrase and RuBisCO activity for each simulation and compare the case with the enzymes evenly distributed throughout the carboxysome to the case where the carbonic anhydrase is localized to the inner carboxysome shell. The (-.-) lines are for no organization and (x) for localization with Dc=105cm2s. The (…) lines are for no organization and (o) are for localization with Dc=107cm2s.

https://doi.org/10.7554/eLife.02043.011
Concentration of CO2 in the carboxysome with varying carboxysome permeability (A).

Numerical solution (diamonds and circles) and analytic solutions (carbonic anhydrase unsaturated, solid lines, and saturated, dashed lines) correspond well. On all plots CO2 (red circle) < HCO3 (blue diamond). Concentration in the cell along the radius, r, with carboxysome permeability kc=105cms (B), kc=103cms (C), kc=1cms (D). Grey dotted lines in (B), (C), (D) indicate location of the carboxysome shell boundary. The transition from low CO2 at high permeability (D) to maximum CO2 concentration at optimal permeability (C) occurs at kc=DRc=2cms. At low carboxysome permeability (B) HCO3 diffusion into the carboxysome is slower than consumption. For all subplots α=0cms and jc=0.6cms. Qualitative results remain the same with varying jc, increasing α will increase the gradient of CO2 across the cell as CO2 is converted to HCO3 at the cell membrane.

https://doi.org/10.7554/eLife.02043.012
Size of the HCO3 flux in one cell from varying sources, as the proportion of CO2 to HCO3 outside the cell changes changes.

We show results for three carboxysome permeabilities, kc, and only the scavenging is effected. Total external inorganic carbon is 15μM, jc=1cms and αKα=1cms. Scavenging is negligibly small for all values of kc shown. Unless there is very little HCO3 in the environment, HCO3 transport seems to be more efficient than CO2 facilitated uptake.

https://doi.org/10.7554/eLife.02043.015
Concentration of CO2 achieved through various cellular organizations of enzymes, where we have selected the HCO3 transport level such that the HCO3 concentration in the cytosol is 30 mM.

O2 concentration is 260 μM. The oxygenation error rates, as a percent of total RuBisCO reactions are indicated on the concentration bars. The cellular organizations investigated are RuBisCO and carbonic anhydrase distributed throughout the entire cytosol, co-localizing RuBisCO and carbonic anhydrase on a scaffold at the center of the cell without a carboxysome shell, RuBisCO and carbonic anhydrase encapsulated in a carboxysome with high permeability at the center of the cell, and RuBisCO and carbonic anhydrase encapsulated in a carboxysome with optimal permeability at the center of the cell.

https://doi.org/10.7554/eLife.02043.004

Tables

Table 1

Parameter values chosen for main set of simulations, unless otherwise indicated

https://doi.org/10.7554/eLife.02043.005
ParameterDefinitionValueReference
Houtconcentration of bicarbonate outside the cell14 μM*(Price et al., 2008)
Coutconcentration of carbon dioxide outside of cell0.14 μM*(Price et al., 2008)
Ddiffusion constant of small molecules, CO2 and HCO310−5 cm2s(Fridlyand et al., 1996)
kmcpermeability of cell membrane to CO20.3 cms(Missner et al., 2008; Gutknecht et al., 1977)
kmHpermeability of cell membrane to HCO33×104cms(Missner et al., 2008; Gutknecht et al., 1977)
Rcradius of carboxysome5×10−6 cm(Cheng et al., 2008; Schmid et al., 2006)
Rbradius of bacteria5 × 10−5 cm(Savage et al., 2010)
jcHCO3 transport rate resulting in 30mM cytosolic HCO3 pool0.6 cms*calculated here
kcoptimal carboxysome permeability10−3 cms*calculated here
Vcellcell volume5.2×1010μLcalculated
SAcellcell surface area3×108cm2calculated
  1. *

    these parameters are varied in the text, but these values are use unless noted otherwise.

Table 2
https://doi.org/10.7554/eLife.02043.006
Enzyme reactionactive siteskcat [1s]Vmax in ‘cell’ [μMs]Vmax in carboxysome [μMs]K1/2 [μM]
carbonic anhydrase hydration808 × 1048.8 × 1031.5 × 1073.2 × 103
carbonic anhydrase dehydration804.6 × 1041.5 × 1048.8 × 1069.3 × 103
RuBisCO carboxylation2160261781.8 × 105270
  1. Vmax in cell and carboxysome refer to the volumetric reaction rate assuming the enzymes are distributed throughout the entire cell or only carboxysome. Vba (Vmax for carbonic anhydrase dehydration) is estimated by assuming Keq = 5 and using parameters found in (Heinhorst et al., 2006). Vca is Vmax for carbonic anhydrase hydration. Similarly, Kba, and Kca are K1/2 for dehydration and hydration respectively.

Table 3

Fate of carbon brought into the cell for jc = 0.6cm/s and kc = 10–3cm/s

https://doi.org/10.7554/eLife.02043.013
formula[picomoles(cells)]% of flux
HCO3 transportjcHout3.26 × 10−4
HCO3 leakagekmH(HoutHcytosol(Rb))3.2 × 10−498.6%
CO2 leakagekmC(CoutCcytosol(Rb))4.5 × 10−41.4 %
carboxylationVmaxCC+Km(1+OKO)8.2 × 10−80.03 %
oxygenationVmaxOCC+KmO(1+CKC)6.7 × 10−102 × 10−4 %
Table 4

Fate of carbon brought into the cell for jc = 0.06 cm/s and kc = 10–3cm/s

https://doi.org/10.7554/eLife.02043.014
formula[picomoles(cells)]% of flux
HCO3 transportjcHout2.8 × 10−5
HCO3 leakagekmH(HoutHcytosol(Rb))2.7 × 10−596.6 %
CO2 leakagekmC(CoutCcytosol(Rb))8.8 × 10−73.2 %
carboxylationVmaxCC+Km(1+OKO)5.4 × 10−80.2 %
oxygenationVmaxOCC+KmO(1+CKC)2.3 × 10−98 × 10−3 %

Additional files

Supplementary file 1

Mathematical derivation appendix. Mathematical derivations of analytic solutions for a spherical cell with reactions organized in a variety of ways. We present analytic solutions for the concentration of CO2 and HCO3 a carboxysome located at the center of the cell. We derive analytic solutions assuming a number of different cases for the enzymatic rates in the carboxysome: RuBisCO reaction rate negligible with carbonic anhydrase saturated and unsaturated, RuBisCO reaction rate non-negligible with carbonic anhydrase unsaturated. Additionally we derive analytic solutions for the enzymatic reactions throughout the cell and localized to a scaffold without a carboxysome shell.

https://doi.org/10.7554/eLife.02043.016

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Niall M Mangan
  2. Michael P Brenner
(2014)
Systems analysis of the CO2 concentrating mechanism in cyanobacteria
eLife 3:e02043.
https://doi.org/10.7554/eLife.02043