Magnus Carlquist

Experimental methods and modeling techniques for description of cell population heterogeneity

With the continuous development, in the last decades, of analytical techniques providing complex information at single cell level, the study of cell heterogeneity has been the focus of several research projects within analytical biotechnology. Nonetheless, the complex interplay between environmental changes and cellular responses is yet not fully understood, and the integration of this new knowledge into the strategies for design, operation and control of bioprocesses is far from being an established reality. Indeed, the impact of cell heterogeneity on productivity of large scale cultivations is acknowledged but seldom accounted for. In order to include population heterogeneity mechanisms in the development of novel bioprocess control strategies, a reliable mathematical description of such phenomena has to be developed. With this review, we search to summarize the potential of currently available methods for monitoring cell population heterogeneity as well as model frameworks suitable for describing dynamic heterogeneous cell populations. We will furthermore underline the highly important coordination between experimental and modeling efforts necessary to attain a reliable quantitative description of cell heterogeneity, which is a necessity if such models are to contribute to the development of improved control of bioprocesses.

Physiological heterogeneities in microbial populations and implications for physical stress tolerance.

BackgroundTraditionally average values of the whole population are considered when analysing microbial cell cultivations. However, a typical microbial population in a bioreactor is heterogeneous in most phenotypes measurable at a single-cell level. There are indications that such heterogeneity may be unfavourable on the one hand (reduces yields and productivities), but also beneficial on the other hand (facilitates quick adaptation to new conditions - i.e. increases the robustness of the fermentation process). Understanding and control of microbial population heterogeneity is thus of major importance for improving microbial cell factory processes.ResultsIn this work, a dual reporter system was developed and applied to map growth and cell fitness heterogeneities within budding yeast populations during aerobic cultivation in well-mixed bioreactors. The reporter strain, which was based on the expression of green fluorescent protein (GFP) under the control of the ribosomal protein RPL22a promoter, made it possible to distinguish cell growth phases by the level of fluorescence intensity. Furthermore, by exploiting the strong correlation of intracellular GFP level and cell membrane integrity it was possible to distinguish subpopulations with high and low cell membrane robustness and hence ability to withstand freeze-thaw stress. A strong inverse correlation between growth and cell membrane robustness was observed, which further supports the hypothesis that cellular resources are limited and need to be distributed as a trade-off between two functions: growth and robustness. In addition, the trade-off was shown to vary within the population, and the occurrence of two distinct subpopulations shifting between these two antagonistic modes of cell operation could be distinguished.ConclusionsThe reporter strain enabled mapping of population heterogeneities in growth and cell membrane robustness towards freeze-thaw stress at different phases of cell cultivation. The described reporter system is a valuable tool for understanding the effect of environmental conditions on population heterogeneity of microbial cells and thereby to understand cell responses during industrial process-like conditions. It may be applied to identify more robust subpopulations, and for developing novel strategies for strain improvement and process design for more effective bioprocessing.

Flavonoids as inhibitors of human carbonyl reductase 1

Human carbonyl reductase 1 (CBR1), that is one of the enzymes responsible for the reduced efficiency of treatments by the antineoplastic agents anthracyclines, was functionally expressed in Saccharomyces cerevisiae. CBR1 was purified and kinetically characterised using daunorubicin as substrate. CBR1-catalysed reduction of daunorubicin followed an apparent Michaelis-Menten kinetics with K(M)=85.2+/-26.7microM and V(max)=3490+/-220micromol/(mingprotein). The type of inhibition for the flavonoid compound rutin was determined by studying initial reaction rates in the presence of rutin. The inhibition kinetics was found to follow an apparent mixed inhibition with K(ic)=1.8+/-1.2microM and K(iu)=2.8+/-1.6microM. IC50-values were also determined for a set of flavonoids in order to identify essential structure for inhibition activity. Computational docking experiments of the four best inhibitors to the catalytic site of CBR1 showed that the flavonoid skeleton structure was the binding part of the molecule. The presence of a sugar moiety in 1 and 2, or a sugar mimicking part in 9, directed the orientation of the flavonoid so that the sugars were pointing outwards, giving rise to a stabilising effect to the binding. Finally, additional binding epitopes that interacted with various parts of the flavonoid ligand were identified and could potentially be targeted for further improvement of inhibition activity. These included; hydrogen-binding sites surrounding Ser139 and Cys226, Met234 and Tyr193 or Trp229; aromatic-aromatic interaction with Tyr193, Trp229 or NADPH; van der Waals interactions with Ile140.

Cell mass and cell cycle dynamics of an asynchronous budding yeast population: experimental observations, flow cytometry data analysis, and multi-scale modeling.

Despite traditionally regarded as identical, cells in a microbial cultivation present a distribution of phenotypic traits, forming a heterogeneous cell population. Moreover, the degree of heterogeneity is notably enhanced by changes in micro‐environmental conditions. A major development in experimental single‐cell studies has taken place in the last decades. It has however not been fully accompanied by similar contributions within data analysis and mathematical modeling. Indeed, literature reporting, for example, quantitative analyses of experimental single‐cell observations and validation of model predictions for cell property distributions against experimental data is scarce. This study focuses on the experimental and mathematical description of the dynamics of cell size and cell cycle position distributions, of a population of Saccharomyces cerevisiae, in response to the substrate consumption observed during batch cultivation. The good agreement between the proposed multi‐scale model (a population balance model [PBM] coupled to an unstructured model) and experimental data (both the overall physiology and cell size and cell cycle distributions) indicates that a mechanistic model is a suitable tool for describing the microbial population dynamics in a bioreactor. This study therefore contributes towards the understanding of the development of heterogeneous populations during microbial cultivations. More generally, it consists of a step towards a paradigm change in the study and description of cell cultivations, where average cell behaviors observed experimentally now are interpreted as a potential joint result of various co‐existing single‐cell behaviors, rather than a unique response common to all cells in the cultivation. Biotechnol. Bioeng. 2013; 110: 812–826.

NADH-dependent biosensor in Saccharomyces cerevisiae: principle and validation at the single cell level.

A reporter system was constructed to measure perturbations in the NADH/NAD+ co-factor balance in yeast, by using the green fluorescent protein gene under the control of the GPD2 promoter that is induced under conditions of excess of NADH. High fluorescence levels were obtained in a glycerol 3-phosphate dehydrogenase double deletion strain (gpd1Δgpd2Δ), which is deficient in the ability to regenerate NAD+ via glycerol formation. The responsiveness of the reporter system to externally induced perturbations in NADH oxidation was also evaluated in the gpd1Δgpd2Δ strain background by addition of acetoin, as well as by introduction of a set of heterologous xylose reductases (XRs) having different selectivities for NADH. Addition of acetoin during cell proliferation under oxygen-limited conditions resulted in a more than 2-fold decrease in mean fluorescence intensity as compared to the control experiment. Strains carrying XRs with different selectivities for NADH could be distinguished at the single cell level, so that the XR with the highest selectivity for NADH displayed the lowest fluorescence. In conclusion, the designed system successfully allowed for monitoring perturbations in the cellular redox metabolism caused by environmental changes, or by heterologous gene expression. The reporter system displayed high resolution in distinguishing cytosolic NADH oxidation capacity and hence has potential to be used for high-throughput screening based on the fluorescence of single cells.

Exploiting cell metabolism for biocatalytic whole-cell transamination by recombinant Saccharomyces cerevisiae.

The potential of Saccharomyces cerevisiae for biocatalytic whole-cell transamination was investigated using the kinetic resolution of racemic 1-phenylethylamine (1-PEA) to (R)-1-PEA as a model reaction. As native yeast do not possess any ω-transaminase activity for the reaction, a recombinant yeast biocatalyst was constructed by overexpressing the gene coding for vanillin aminotransferase from Capsicum chinense. The yeast-based biocatalyst could use glucose as the sole co-substrate for the supply of amine acceptor via cell metabolism. In addition, the biocatalyst was functional without addition of the co-factor pyridoxal-5′-phosphate (PLP), which can be explained by a high inherent cellular capacity to sustain PLP-dependent reactions in living cells. In contrast, external PLP supplementation was required when cell viability was low, as it was the case when using pyruvate as a co-substrate. Overall, the results indicate a potential for engineered S. cerevisiae as a biocatalyst for whole-cell transamination and with glucose as the only co-substrate for the supply of amine acceptor and PLP.

Biocatalytic potential of vanillin aminotransferase from Capsicum chinense

BackgroundThe conversion of vanillin to vanillylamine is a key step in the biosynthetic route towards capsaicinoids in pungent cultivars of Capsicum sp. The reaction has previously been annotated to be catalysed by PAMT (putative aminotransferase; [GenBank: AAC78480.1, Swiss-Prot: O82521]), however, the enzyme has previously not been biochemically characterised in vitro.ResultsThe biochemical activity of the transaminase was confirmed by direct measurement of the reaction with purified recombinant enzyme. The enzyme accepted pyruvate, and oxaloacetate but not 2-oxoglutarate as co-substrate, which is in accordance with other characterised transaminases from the plant kingdom. The enzyme was also able to convert (S)-1-phenylethylamine into acetophenone with high stereo-selectivity. Additionally, it was shown to be active at a broad pH range.ConclusionsWe suggest PAMT to be renamed to VAMT (vanillin aminotransferase, abbreviation used in this study) as formation of vanillin from vanillylamine could be demonstrated. Furthermore, due to high stereoselectivity and activity at physiological pH, VAMT is a suitable candidate for biocatalytic transamination in a recombinant whole-cell system.

Applying Mechanistic Models in Bioprocess Development

The available knowledge on the mechanisms of a bioprocess system is central to process analytical technology. In this respect, mechanistic modeling has gained renewed attention, since a mechanistic model can provide an excellent summary of available process knowledge. Such a model therefore incorporates process-relevant input (critical process variables)-output (product concentration and product quality attributes) relations. The model therefore has great value in planning experiments, or in determining which critical process variables need to be monitored and controlled tightly. Mechanistic models should be combined with proper model analysis tools, such as uncertainty and sensitivity analysis. When assuming distributed inputs, the resulting uncertainty in the model outputs can be decomposed using sensitivity analysis to determine which input parameters are responsible for the major part of the output uncertainty. Such information can be used as guidance for experimental work; i.e., only parameters with a significant influence on model outputs need to be determined experimentally. The use of mechanistic models and model analysis tools is demonstrated in this chapter. As a practical case study, experimental data from Saccharomyces cerevisiae fermentations are used. The data are described with the well-known model of Sonnleitner and Käppeli (Biotechnol Bioeng 28:927-937, 1986) and the model is analyzed further. The methods used are generic, and can be transferred easily to other, more complex case studies as well.

Comparison of engineered Saccharomyces cerevisiae and engineered Escherichia coli for the production of an optically pure keto alcohol.

In this study, the production of enantiomerically pure (1R,4S,6S)-6-hydroxy-bicyclo[2.2.2]octane-2-one ((−)-2) through stereoselective bioreduction was used as a model reaction for the comparison of engineered Saccharomyces cerevisiae and engineered Escherichia coli as biocatalysts. For both microorganisms, over-expression of the gene encoding the NADPH-dependent aldo-keto reductase YPR1 resulted in high purity of the keto alcohol (−)-2 (>99% ee, 97–98% de). E. coli had three times higher initial reduction rate but S. cerevisiae continued the reduction reaction for a longer time period, thus reaching a higher conversion of the substrate (95%). S. cerevisiae was also more robust than E. coli, as demonstrated by higher viability during bioreduction. It was also investigated whether the NADPH regeneration rate was sufficient to supply the over-expressed reductase with NADPH. Five strains of each microorganism with varied carbon flux through the NADPH regenerating pentose phosphate pathway were genetically constructed and compared. S. cerevisiae required an increased NADPH regeneration rate to supply YPR1 with co-enzyme while the native NADPH regeneration rate was sufficient for E. coli.

Engineered baker’s yeast as whole-cell biocatalyst for one-pot stereo-selective conversion of amines to alcohols

BackgroundOne-pot multi-step biocatalysis is advantageous over step-by-step synthesis as it reduces the number of process operation units, leading to significant process intensification. Whole-cell biocatalysis with metabolically active cells is especially valuable since all enzymes can be co-expressed in the cell whose metabolism can be exploited for supply of co-substrates and co-factors.ResultsIn this study, a heterologous enzymatic system consisting of ω-transaminase and ketone reductase was introduced in Saccharomyces cerevisiae, and evaluated for one-pot stereo-selective conversion of amines to alcohols. The system was applied for simultaneous kinetic resolution of racemic 1-phenylethylamine to (R)-1-phenylethylamine and reduction of the ketone intermediate to (R)-1-phenylethanol. Glucose was used as sole co-substrate for both the supply of amine acceptor and the regeneration of NADPH in the reduction step.ConclusionsThe whole-cell biocatalyst was shown to sustain transaminase-reductase-catalyzed enantioselective conversion of amines to alcohols with glucose as co-substrate. The transamination catalyzed by recombinant vanillin aminotransferase from Capsicum chinense proved to be the rate-limiting step as a three-fold increase in transaminase gene copy number led to a two-fold increased conversion. The (R)-selective NADPH-dependent alcohol dehydrogenase from Lactobacillus kefir proved to be efficient in catalyzing the reduction of the acetophenone generated in the transamination reaction.