Computational and biochemical methods to measure the activity of carboxysomes and protein organelles in vivo.

Cyanobacteria are photosynthetic microorganisms that play important ecological roles as major contributors to global nutrient cycles. Cyanobacteria are highly efficient in carrying out oxygenic photosynthesis because they possess carboxysomes, a class of bacterial microcompartments (BMC) in which a polyhedral protein shell encapsulates the enzymes ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) and carbonic anhydrase and functions as the key component of the cyanobacterial CO2-concentrating mechanism (CCM). Elevated CO2 levels within the carboxysome shell as a result of carbonic anhydrase activity increase the efficiency of RuBisCO. Yet, there remain many questions regarding the flux or exclusion of metabolites across the shell and how the activity of BMCs varies over time. These questions have been difficult to address using traditional ensemble techniques due to the heterogeneity of BMCs extracted from their native hosts or with heterologous expression. In this chapter, we describe a method to film and extract quantitative information about carboxysome activity using molecular biology and live cell, timelapse microscopy. In our method, the production of carboxysomes is first controlled by deleting the native genes required for carboxysome assembly and then re-introducing them under the control of an inducible promoter. This system enables carboxysomes to be tracked through multiple generations of cells and provides a way to quantify the total biomass accumulation attributed to a single carboxysome. While the method presented here was developed specifically for carboxysomes, it could be modified to track and quantify the activity of bacterial microcompartments in general.

Role of carboxysomes in cyanobacterial CO2 assimilation: CO2 concentrating mechanisms and metabolon implications.

Many carbon-fixing organisms have evolved CO2 concentrating mechanisms (CCMs) to enhance the delivery of CO2 to RuBisCO, while minimizing reactions with the competitive inhibitor, molecular O2 . These distinct types of CCMs have been extensively studied using genetics, biochemistry, cell imaging, mass spectrometry, and metabolic flux analysis. Highlighted in this paper, the cyanobacterial CCM features a bacterial microcompartment (BMC) called 'carboxysome' in which RuBisCO is co-encapsulated with the enzyme carbonic anhydrase (CA) within a semi-permeable protein shell. The cyanobacterial CCM is capable of increasing CO2 around RuBisCO, leading to one of the most efficient processes known for fixing ambient CO2 . The carboxysome life cycle is dynamic and creates a unique subcellular environment that promotes activity of the Calvin-Benson (CB) cycle. The carboxysome may function within a larger cellular metabolon, physical association of functionally coupled proteins, to enhance metabolite channelling and carbon flux. In light of CCMs, synthetic biology approaches have been used to improve enzyme complex for CO2 fixations. Research on CCM-associated metabolons has also inspired biologists to engineer multi-step pathways by providing anchoring points for enzyme cascades to channel intermediate metabolites towards valuable products.

Asymmetric survival in single-cell lineages of cyanobacteria in response to photodamage.

Oxygenic photosynthesis is driven by the coupled action of the light-dependent pigment-protein complexes, photosystem I and II, located within the internal thylakoid membrane system. However, photosystem II is known to be prone to photooxidative damage. Thus, photosynthetic organisms have evolved a repair cycle to continuously replace the damaged proteins in photosystem II. However, it has remained difficult to deconvolute the damage and repair processes using traditional ensemble approaches. Here, we demonstrate an automated approach using time-lapse fluorescence microscopy and computational image analysis to study the dynamics and effects of photodamage in single cells at subcellular resolution in cyanobacteria. By growing cells in a two-dimensional layer, we avoid shading effects, thereby generating uniform and reproducible growth conditions. Using this platform, we analyzed the growth and physiology of multiple strains simultaneously under defined photoinhibitory conditions stimulated by UV-A light. Our results reveal an asymmetric cellular response to photodamage between sibling cells and the generation of an elusive subcellular structure, here named a 'photoendosome,' derived from the thylakoid which could indicate the presence of a previously unknown photoprotective mechanism. We anticipate these results to be a starting point for further studies to better understand photodamage and repair at the single-cell level.

Dynamic and single cell characterization of a CRISPR-interference toolset in Pseudomonas putida KT2440 for ß-ketoadipate production from p-coumarate.

Pseudomonas putida KT2440 is a well-studied bacterium for the conversion of lignin-derived aromatic compounds to bioproducts. The development of advanced genetic tools in P. putida has reduced the turnaround time for hypothesis testing and enabled the construction of strains capable of producing various products of interest. Here, we evaluate an inducible CRISPR-interference (CRISPRi) toolset on fluorescent, essential, and metabolic targets. Nuclease-deficient Cas9 (dCas9) expressed with the arabinose (8K)-inducible promoter was shown to be tightly regulated across various media conditions and when targeting essential genes. In addition to bulk growth data, single cell time lapse microscopy was conducted, which revealed intrinsic heterogeneity in knockdown rate within an isoclonal population. The dynamics of knockdown were studied across genomic targets in exponentially-growing cells, revealing a universal 1.75 ± 0.38 h quiescent phase after induction where 1.5 ± 0.35 doublings occur before a phenotypic response is observed. To demonstrate application of this CRISPRi toolset, ß-ketoadipate, a monomer for performance-advantaged nylon, was produced at a 4.39 ± 0.5 g/L and yield of 0.76 ± 0.10 mol/mol from p-coumarate, a hydroxycinnamic acid that can be derived from grasses. These cultivation metrics were achieved by using the higher strength IPTG (1K)-inducible promoter to knockdown the pcaIJ operon in the ßKA pathway during early exponential phase. This allowed the majority of the carbon to be shunted into the desired product while eliminating the need for a supplemental carbon and energy source to support growth and maintenance.

Zam Is a Redox-Regulated Member of the RNB-Family Required for Optimal Photosynthesis in Cyanobacteria.

The zam gene mediating resistance to acetazolamide in cyanobacteria was discovered thirty years ago during a drug tolerance screen. We use phylogenetics to show that Zam proteins are distributed across cyanobacteria and that they form their own unique clade of the ribonuclease II/R (RNB) family. Despite being RNB family members, multiple sequence alignments reveal that Zam proteins lack conservation and exhibit extreme degeneracy in the canonical active site-raising questions about their cellular function(s). Several known phenotypes arise from the deletion of zam, including drug resistance, slower growth, and altered pigmentation. Using room-temperature and low-temperature fluorescence and absorption spectroscopy, we show that deletion of zam results in decreased phycocyanin synthesis rates, altered PSI:PSII ratios, and an increase in coupling between the phycobilisome and PSII. Conserved cysteines within Zam are identified and assayed for function using in vitro and in vivo methods. We show that these cysteines are essential for Zam function, with mutation of either residue to serine causing phenotypes identical to the deletion of Zam. Redox regulation of Zam activity based on the reversible oxidation-reduction of a disulfide bond involving these cysteine residues could provide a mechanism to integrate the 'central dogma' with photosynthesis in cyanobacteria.

Computational modeling and evolutionary implications of biochemical reactions in bacterial microcompartments.

Bacterial microcompartments (BMCs) are protein-encapsulated compartments found across at least 23 bacterial phyla. BMCs contain a variety of metabolic processes that share the commonality of toxic or volatile intermediates, oxygen-sensitive enzymes and cofactors, or increased substrate concentration for magnified reaction rates. These compartmentalized reactions have been computationally modeled to explore the encapsulated dynamics, ask evolutionary-based questions, and develop a more systematic understanding required for the engineering of novel BMCs. Many crucial aspects of these systems remain unknown or unmeasured, such as substrate permeabilities across the protein shell, feasibility of pH gradients, and transport rates of associated substrates into the cell. This review explores existing BMC models, dominated in the literature by cyanobacterial carboxysomes, and highlights potentially important areas for exploration.

ANALYTICAL SINGULAR VALUE DECOMPOSITION FOR A CLASS OF STOICHIOMETRY MATRICES.

We present the analytical singular value decomposition of the stoichiometry matrix for a spatially discrete reaction-diffusion system. The motivation for this work is to develop a matrix decomposition that can reveal hidden spatial flux patterns of chemical reactions. We consider a 1D domain with two subregions sharing a single common boundary. Each of the subregions is further partitioned into a finite number of compartments. Chemical reactions can occur within a compartment, whereas diffusion is represented as movement between adjacent compartments. Inspired by biology, we study both (1) the case where the reactions on each side of the boundary are different and only certain species diffuse across the boundary and (2) the case where reactions and diffusion are spatially homogeneous. We write the stoichiometry matrix for these two classes of systems using a Kronecker product formulation. For the first scenario, we apply linear perturbation theory to derive an approximate singular value decomposition in the limit as diffusion becomes much faster than reactions. For the second scenario, we derive an exact analytical singular value decomposition for all relative diffusion and reaction time scales. By writing the stoichiometry matrix using Kronecker products, we show that the singular vectors and values can also be written concisely using Kronecker products. Ultimately, we find that the singular value decomposition of the reaction-diffusion stoichiometry matrix depends on the singular value decompositions of smaller matrices. These smaller matrices represent modified versions of the reaction-only stoichiometry matrices and the analytically known diffusion-only stoichiometry matrix. Lastly, we present the singular value decomposition of the model for the Calvin cycle in cyanobacteria and demonstrate the accuracy of our formulation. The MATLAB code, available at www.github.com/MathBioCU/ReacDiffStoicSVD, provides routines for efficiently calculating the SVD for a given reaction network on a 1D spatial domain.

Carbon isotope evidence for the global physiology of Proterozoic cyanobacteria.

Ancestral cyanobacteria are assumed to be prominent primary producers after the Great Oxidation Event [≈2.4 to 2.0 billion years (Ga) ago], but carbon isotope fractionation by extant marine cyanobacteria (α-cyanobacteria) is inconsistent with isotopic records of carbon fixation by primary producers in the mid-Proterozoic eon (1.8 to 1.0 Ga ago). To resolve this disagreement, we quantified carbon isotope fractionation by a wild-type planktic ß-cyanobacterium (Synechococcus sp. PCC 7002), an engineered Proterozoic analog lacking a CO2-concentrating mechanism, and cyanobacterial mats. At mid-Proterozoic pH and pCO2 values, carbon isotope fractionation by the wild-type ß-cyanobacterium is fully consistent with the Proterozoic carbon isotope record, suggesting that cyanobacteria with CO2-concentrating mechanisms were apparently the major primary producers in the pelagic Proterozoic ocean, despite atmospheric CO2 levels up to 100 times modern. The selectively permeable microcompartments central to cyanobacterial CO2-concentrating mechanisms ("carboxysomes") likely emerged to shield rubisco from O2 during the Great Oxidation Event.

Engineering living building materials for enhanced bacterial viability and mechanical properties.

Living building materials (LBMs) utilize microorganisms to produce construction materials that exhibit mechanical and biological properties. A hydrogel-based LBM containing bacteria capable of microbially induced calcium carbonate precipitation (MICP) was recently developed. Here, LBM design factors, i.e., gel/sand ratio, inclusion of trehalose, and MICP pathways, are evaluated. The results show that non-saturated LBM (gel/sand = 0.13) and gel-saturated LBM (gel/sand = 0.30) underwent distinct failure modes. The inclusion of trehalose maintains bacterial viability under ambient conditions with low relative humidity, without affecting mechanical properties of the LBM. Comparison of biotic and abiotic LBM shows that MICP efficiency in this material is subject to the pathway selected: the LBM with heterotrophic ureolytic Escherichia coli demonstrated the most mechanical enhancement from the abiotic controls, compared with either ureolytic or CO2-concentrating mechanisms from Synechococcus. The study shows that tailoring of LBM properties can be accomplished in a manner that considers both LBM microstructure and MICP pathways.

Effect of pH on the activity of ice-binding protein from Marinomonas primoryensis.

The ability of an ice-binding protein (IBP) from Marinomonas primoryensis (MpIBP) to influence ice crystal growth and structure in nonphysiological pH environments was investigated in this work. The ability for MpIBP to retain ice interactivity under stressed environmental conditions was determined via (1) a modified splat assay to determine ice recrystallization inhibition (IRI) of polycrystalline ice and (2) nanoliter osmometry to evaluate the ability of MpIBP to dynamically shape the morphology of a single ice crystal. Circular dichroism (CD) was used to relate the IRI and DIS activity of MpIBP to secondary structure. The results illustrate that MpIBP secondary structure was stable between pH 6 and pH 10. It was found that MpIBP did not interact with ice at pH ≤ 4 or pH ≥ 13. At 6 ≤ pH ≥ 12 MpIBP exhibited a reduction in grain size of ice crystals as compared to control solutions and demonstrated dynamic ice shaping at 6 ≤ pH ≥ 10. The results substantiate that MpIBP retains some secondary structure and function in non-neutral pH environments; thereby, enabling its potential utility in nonphysiological materials science and engineering applications.

Proximity-based proteomics reveals the thylakoid lumen proteome in the cyanobacterium Synechococcus sp. PCC 7002.

Cyanobacteria possess unique intracellular organization. Many proteomic studies have examined different features of cyanobacteria to learn about the intracellular structures and their respective functions. While these studies have made great progress in understanding cyanobacterial physiology, the conventional fractionation methods used to purify cellular structures have limitations; specifically, certain regions of cells cannot be purified with existing fractionation methods. Proximity-based proteomics techniques were developed to overcome the limitations of biochemical fractionation for proteomics. Proximity-based proteomics relies on spatiotemporal protein labeling followed by mass spectrometry of the labeled proteins to determine the proteome of the region of interest. We performed proximity-based proteomics in the cyanobacterium Synechococcus sp. PCC 7002 with the APEX2 enzyme, an engineered ascorbate peroxidase. We determined the proteome of the thylakoid lumen, a region of the cell that has remained challenging to study with existing methods, using a translational fusion between APEX2 and PsbU, a lumenal subunit of photosystem II. Our results demonstrate the power of APEX2 as a tool to study the cell biology of intracellular features and processes, including photosystem II assembly in cyanobacteria, with enhanced spatiotemporal resolution.

Genome-Wide Analysis of RNA Decay in the Cyanobacterium Synechococcus sp. Strain PCC 7002.

RNA degradation is an important process that influences the ultimate concentration of individual proteins inside cells. While the main enzymes that facilitate this process have been identified, global maps of RNA turnover are available for only a few species. Even in these cases, there are few sequence elements that are known to enhance or destabilize a native transcript; even fewer confer the same effect when added to a heterologous transcript. To address this knowledge gap, we assayed genome-wide RNA degradation in the cyanobacterium Synechococcus sp. strain PCC 7002 by collecting total RNA samples after stopping nascent transcription with rifampin. We quantified the abundance of each position in the transcriptome as a function of time using RNA-sequencing data and later analyzed the global mRNA decay map using machine learning principles. Half-lives, calculated on a per-ORF (open reading frame) basis, were extremely short, with a median half-life of only 0.97 min. Despite extremely rapid turnover of most mRNA, transcripts encoding proteins involved in photosynthesis were both highly expressed and highly stable. Upon inspection of these stable transcripts, we identified an enriched motif in the 3' untranslated region (UTR) that had similarity to Rho-independent terminators. We built statistical models for half-life prediction and used them to systematically identify sequence motifs in both 5' and 3' UTRs that correlate with stabilized transcripts. We found that transcripts linked to a terminator containing a poly(U) tract had a longer half-life than both those without a poly(U) tract and those without a terminator.IMPORTANCE RNA degradation is an important process that affects the final concentration of individual mRNAs, affecting protein expression and cellular physiology. Studies of how RNA is degraded increase our knowledge of this fundamental process as well as enable the creation of genetic tools to manipulate RNA stability. By studying global transcript turnover, we searched for sequence elements that correlated with transcript (in)stability and used these sequences to guide tool design. This study probes global RNA turnover in a cyanobacterium, Synechococcus sp. strain PCC 7002, that both has a unique array of RNases that facilitate RNA degradation and is an industrially relevant strain that could be used to convert CO2 and sunlight into useful products.

Life cycle of a cyanobacterial carboxysome.

Carboxysomes, prototypical bacterial microcompartments (BMCs) found in cyanobacteria, are large (~1 GDa) and essential protein complexes that enhance CO2 fixation. While carboxysome biogenesis has been elucidated, the activity dynamics, lifetime, and degradation of these structures have not been investigated, owing to the inability of tracking individual BMCs over time in vivo. We have developed a fluorescence-imaging platform to simultaneously measure carboxysome number, position, and activity over time in a growing cyanobacterial population, allowing individual carboxysomes to be clustered on the basis of activity and spatial dynamics. We have demonstrated both BMC degradation, characterized by abrupt activity loss followed by polar recruitment of the deactivated complex, and a subclass of ultraproductive carboxysomes. Together, our results reveal the BMC life cycle after biogenesis and describe the first method for measuring activity of single BMCs in vivo.

Mechanical regulation of photosynthesis in cyanobacteria.

Photosynthetic organisms regulate their responses to many diverse stimuli in an effort to balance light harvesting with utilizable light energy for carbon fixation and growth (source-sink regulation). This balance is critical to prevent the formation of reactive oxygen species that can lead to cell death. However, investigating the molecular mechanisms that underlie the regulation of photosynthesis in cyanobacteria using ensemble-based measurements remains a challenge due to population heterogeneity. Here, to address this problem, we used long-term quantitative time-lapse fluorescence microscopy, transmission electron microscopy, mathematical modelling and genetic manipulation to visualize and analyse the growth and subcellular dynamics of individual wild-type and mutant cyanobacterial cells over multiple generations. We reveal that mechanical confinement of actively growing Synechococcus sp. PCC 7002 cells leads to the physical disassociation of phycobilisomes and energetic decoupling from the photosynthetic reaction centres. We suggest that the mechanical regulation of photosynthesis is a critical failsafe that prevents cell expansion when light and nutrients are plentiful, but when space is limiting. These results imply that cyanobacteria must convert a fraction of the available light energy into mechanical energy to overcome frictional forces in the environment, providing insight into the regulation of photosynthesis and how microorganisms navigate their physical environment.

Development of both type I-B and type II CRISPR/Cas genome editing systems in the cellulolytic bacterium Clostridium thermocellum.

The robust lignocellulose-solubilizing activity of C. thermocellum makes it a top candidate for consolidated bioprocessing for biofuel production. Genetic techniques for C. thermocellum have lagged behind model organisms thus limiting attempts to improve biofuel production. To improve our ability to engineer C. thermocellum, we characterized a native Type I-B and heterologous Type II Clustered Regularly-Interspaced Short Palindromic Repeat (CRISPR)/cas (CRISPR associated) systems. We repurposed the native Type I-B system for genome editing. We tested three thermophilic Cas9 variants (Type II) and found that GeoCas9, isolated from Geobacillus stearothermophilus, is active in C. thermocellum. We employed CRISPR-mediated homology directed repair to introduce a nonsense mutation into pyrF. For both editing systems, homologous recombination between the repair template and the genome appeared to be the limiting step. To overcome this limitation, we tested three novel thermophilic recombinases and demonstrated that exo/beta homologs, isolated from Acidithiobacillus caldus, are functional in C. thermocellum. For the Type I-B system an engineered strain, termed LL1586, yielded 40% genome editing efficiency at the pyrF locus and when recombineering machinery was expressed this increased to 71%. For the Type II GeoCas9 system, 12.5% genome editing efficiency was observed and when recombineering machinery was expressed, this increased to 94%. By combining the thermophilic CRISPR system (either Type I-B or Type II) with the recombinases, we developed a new tool that allows for efficient CRISPR editing. We are now poised to enable CRISPR technologies to better engineer C. thermocellum for both increased lignocellulose degradation and biofuel production.

Genome engineering of E. coli for improved styrene production.

Microbial production of exogenous organic compounds is challenging as biosynthetic pathways are often complex and produce metabolites that are toxic to the hosts. Biogenic styrene is an example of this problem, which if addressed could result in a more sustainable supply of this important component of the plastics industry. In this study, we engineered Escherichia coli for the production of styrene. We systematically optimized the production capability by first screening different pathway expression levels in E. coli strains. We then further designed and constructed a transcription regulator library targeting 54 genes with 85,420 mutations, and tested this library for increased styrene resistance and production. A series of tolerant mutants not only exhibited improved styrene tolerance but also produced higher styrene concentrations compared to the parent strain. The best producing mutant, ST05 LexA_E45I, produced a 3.45-fold increase in styrene compared to the parent strain. The produced styrene was extracted via gas stripping into dodecane and used in a direct free radical synthesis of polystyrene.

Engineered Ureolytic Microorganisms Can Tailor the Morphology and Nanomechanical Properties of Microbial-Precipitated Calcium Carbonate.

We demonstrate for the first time that the morphology and nanomechanical properties of calcium carbonate (CaCO3) can be tailored by modulating the precipitation kinetics of ureolytic microorganisms through genetic engineering. Many engineering applications employ microorganisms to produce CaCO3. However, control over bacterial calcite morphology and material properties has not been demonstrated. We hypothesized that microorganisms genetically engineered for low urease activity would achieve larger calcite crystals with higher moduli. We compared precipitation kinetics, morphology, and nanomechanical properties for biogenic CaCO3 produced by two Escherichia coli (E. coli) strains that were engineered to display either high or low urease activity and the native producer Sporosarcina pasteurii. While all three microorganisms produced calcite, lower urease activity was associated with both slower initial calcium depletion rate and increased average calcite crystal size. Both calcite crystal size and nanoindentation moduli were also significantly higher for the low-urease activity E. coli compared with the high-urease activity E. coli. The relative resistance to inelastic deformation, measured via the ratio of nanoindentation hardness to modulus, was similar across microorganisms. These findings may enable design of novel advanced engineering materials where modulus is tailored to the application while resistance to irreversible deformation is not compromised.

Cyanobacterial carboxysome mutant analysis reveals the influence of enzyme compartmentalization on cellular metabolism and metabolic network rigidity.

Cyanobacterial carboxysomes encapsulate carbonic anhydrase and ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). Genetic deletion of the major structural proteins encoded within the ccm operon in Synechococcus sp. PCC 7002 (ΔccmKLMN) disrupts carboxysome formation and significantly affects cellular physiology. Here we employed both metabolite pool size analysis and isotopically nonstationary metabolic flux analysis (INST-MFA) to examine metabolic regulation in cells lacking carboxysomes. Under high CO2 environments (1%), the ΔccmKLMN mutant could recover growth and had a similar central flux distribution as the control strain, with the exceptions of moderately decreased photosynthesis and elevated biomass protein content and photorespiration activity. Metabolite analyses indicated that the ΔccmKLMN strain had significantly larger pool sizes of pyruvate (>18 folds), UDPG (uridine diphosphate glucose), and aspartate as well as higher levels of secreted organic acids (e.g., malate and succinate). These results suggest that the ΔccmKLMN mutant is able to largely maintain a fluxome similar to the control strain by changing in intracellular metabolite concentrations and metabolite overflows under optimal growth conditions. When CO2 was insufficient (0.2%), provision of acetate moderately promoted mutant growth. Interestingly, the removal of microcompartments may loosen the flux network and promote RuBisCO side-reactions, facilitating redirection of central metabolites to competing pathways (i.e., pyruvate to heterologous lactate production). This study provides important insights into metabolic regulation via enzyme compartmentation and cyanobacterial compensatory responses.

Impact of overexpression of cytosolic isoform of O-acetylserine sulfhydrylase on soybean nodulation and nodule metabolome.

Nitrogen-fixing nodules, which are also major sites of sulfur assimilation, contribute significantly to the sulfur needs of whole soybean plants. Nodules are the predominant sites for cysteine accumulation and the activity of O-acetylserine(thiol)lyase (OASS) is central to the sulfur assimilation process in plants. Here, we examined the impact of overexpressing OASS on soybean nodulation and nodule metabolome. Overexpression of OASS did not affect the nodule number, but negatively impacted plant growth. HPLC measurement of antioxidant metabolites demonstrated that levels of cysteine, glutathione, and homoglutathione nearly doubled in OASS overexpressing nodules when compared to control nodules. Metabolite profiling by LC-MS and GC-MS demonstrated that several metabolites related to serine, aspartate, glutamate, and branched-chain amino acid pathways were significantly elevated in OASS overexpressing nodules. Striking differences were also observed in the flavonoid levels between the OASS overexpressing and control soybean nodules. Our results suggest that OASS overexpressing plants compensate for the increase in carbon requirement for sulfur assimilation by reducing the biosynthesis of some amino acids, and by replenishing the TCA cycle through fatty acid hydrolysis. These data may indicate that in OASS overexpressing soybean nodules there is a moderate decease in the supply of energy metabolites to the nodule, which is then compensated by the degradation of cellular components to meet the needs of the nodule energy metabolism.