Effects of in situ cobalt ion addition on the activity of a gfp-oph fusion protein: the fermentation kinetics

The effects of cobalt ion addition and inducer concentration were studied in the fermentation of E. coli BL21 expressing a GFP (green fluorescent protein)-OPH (organophosphorus hydrolase) fusion protein. It was found that cobalt ion addition improved the OPH activity significantly. When 2 mM of CoCl(2) was supplied during the IPTG-induction phase, OPH activity was enhanced approximately 10-fold compared to the case without cobalt or by the addition of cobalt to the cell extracts. Results indicate, therefore, that incorporation of the cobalt during synthesis is needed for enhanced activity. Also, the maximum OPH activity was not linearly related to inducer concentration. A mathematical model was then constructed to simulate these phenomena. Model parameters were determined by constrained least-squares and optimal IPTG and cobalt addition concentrations were obtained, pinpointing the conditions for the maximum productivity. Finally, the GFP fluorescence intensity was found linear to the OPH activity in each fermentation, demonstrating the function of GFP for monitoring its fusion partner's quantity in the bioreactor.     

Environmentally-modulated changes in fluorescence distribution in cells with oscillatory genetic network dynamics

We investigated the distribution of green fluorescent protein (GFP) expression levels in a population of E. coli cells expressing an artificial genetic regulatory network, known as the “repressilator”. This network originally constructed by Elowitz and Leibler in 2000 consists of three cyclically-inhibiting promoter–repressor pairs. It is because of this architecture that the network has been known to oscillate at the single-cell level under certain conditions. A series of shake flask experiments were performed and analyzed using flow cytometry to test how cell populations carrying this system could be controlled extracellularly using the inducers anhydrotetracycline (aTc) and isopropyl-β-d-thiogalactopyranoside (IPTG). With variation of [aTc], it exhibits a novel bi-threshold behavior, such that the entire culture reaches one of three steady states at a quasi-time-invariant “reference state.” Also, there is significant hysteresis. Transiently, the middle state shows damping oscillations, while the low and high states show a stable steady state. The addition of IPTG serves to fine-tune the characteristics of the aTc-only expression, lowering the average and coefficient of variation (CV) of the distributions, and possibly perturbing the network to a different state. However, in modeling this system, the multiplicity and bi-threshold behavior are not theoretically possible according to the designed interactions. In order to explain this discrepancy, we hypothesize that one or more of the repressors have a significant nonspecific interaction with a promoter that does not contain its operator site. The new modeling results incorporating these extra interactions qualitatively match our experimental findings. After constructing plasmids to test these hypotheses, we discover that at least four of these interactions exist, which can create the low and high states and multiplicity seen experimentally. This genetic architecture has flexibility in its behavior that has not been demonstrated before, and the combination of experiment and modeling enlightened our understanding of the molecular interactions driving the network's behavior, leading us to discover the significance of nonspecific interactions.     

Detection of promoter activity by flow cytometric analysis of GFP reporter expression

Low efficiency of transfection is often the limiting factor for acquiring conclusive data in reporter assays. It is especially difficult to efficiently transfect and characterize promoters in primary human cells. To overcome this problem we have developed a system in which reporter gene expression is quantified by flow cytometry. In this system, green fluorescent protein (GFP) reporter constructs are co-transfected with a reference plasmid that codes for the mouse cell surface antigen Thy-1.1 and serves to determine transfection efficiency. Comparison of mean GFP expression of the total transfected cell population with the activity of an analogous luciferase reporter showed that the sensitivity of the two reporter systems is similar. However, because GFP expression can be analyzed at the single-cell level and in the same cells the expression of the reference plasmid can be monitored by two-color fluorescence, the GFP reporter system is in fact more sensitive, particularly in cells which can only be transfected with a low efficiency.     

A quantitative genetic competition model for sympatric speciation

I use multilocus genetics to describe assortative mating in a competition model. The intensity of competition between individuals is influenced by a quantitative character whose value is determined additively by alleles from many loci. With assortative mating based on this character, frequency- and density-dependent competition can subdivide a population with an initially unimodal character distribution. The character distribution becomes bimodal, and the subpopulations corresponding to the two modes are reproductively separated because mating is assortative. This happens if the resource distribution is unimodal, i.e. even if selection due to phenotypic carrying capacities is not disruptive. The results suggest that sympatric speciation due to frequency-dependent selection can occur in quite general ecological scenarios if mating is assortative. I also discuss the evolution of assortative mating. Since it induces bimodal phenotype distributions, assortative mating leads to a better match of the resources if their distribution is also bimodal. Moreover, in a population with a bimodal phenotype distribution, the average strength of frequency-dependent competition is lower than in a unimodal population. Therefore, assortative mating permits higher equilibrium densities than random mating even if the resource distribution is unimodal. Thus, even though it may lead to a less efficient resource use, assortative mating is favoured over random mating because it reduces frequency-dependent effects of competition.     

Evidencing the role of lactose permease in IPTG uptake by Escherichia coli in fed-batch high cell density cultures

The lac-operon and its components have been studied for decades and it is widely used as one of the common systems for recombinant protein production in Escherichia coli. However, the role of the lactose permease, encoded by the lacY gene, when using the gratuitous inducer IPTG for the overexpression of heterologous proteins, is still a matter of discussion. A lactose permease deficient strain was successfully constructed. Growing profiles and acetate production were compared with its parent strain at shake flask scale. Our results show that the lac-permease deficient strain grows slower than the parent in defined medium at shake flask scale, probably due to a downregulation of the phosphotransferase system (PTS). The distributions of IPTG in the medium and inside the cells, as well as recombinant protein production were measured by HPLC-MS and compared in substrate limiting fed-batch fermentations at different inducer concentrations. For the mutant strain, IPTG concentration in the medium depletes slower, reaching at the end of the culture higher concentration values compared with the parent strain. Final intracellular and medium concentrations of IPTG were similar for the mutant strain, while higher intracellular concentrations than in medium were found for the parent strain. Comparison of the distribution profiles of IPTG of both strains in fed-batch fermentations showed that lac-permease is crucially involved in IPTG uptake. In the absence of the transporter, apparently IPTG only diffuses, while in the presence of lac-permease, the inducer accumulates in the cytoplasm at higher rates emphasizing the significant contribution of the permease-mediated transport.     

Bistability versus bimodal distributions in gene regulatory processes from population balance

In recent times, stochastic treatments of gene regulatory processes have appeared in the literature in which a cell exposed to a signaling molecule in its environment triggers the synthesis of a specific protein through a network of intracellular reactions. The stochastic nature of this process leads to a distribution of protein levels in a population of cells as determined by a Fokker-Planck equation. Often instability occurs as a consequence of two (stable) steady state protein levels, one at the low end representing the "off" state, and the other at the high end representing the "on" state for a given concentration of the signaling molecule within a suitable range. A consequence of such bistability has been the appearance of bimodal distributions indicating two different populations, one in the "off" state and the other in the "on" state. The bimodal distribution can come about from stochastic analysis of a single cell. However, the concerted action of the population altering the extracellular concentration in the environment of individual cells and hence their behavior can only be accomplished by an appropriate population balance model which accounts for the reciprocal effects of interaction between the population and its environment. In this study, we show how to formulate a population balance model in which stochastic gene expression in individual cells is incorporated. Interestingly, the simulation of the model shows that bistability is neither sufficient nor necessary for bimodal distributions in a population. The original notion of linking bistability with bimodal distribution from single cell stochastic model is therefore only a special consequence of a population balance model.     

Fluorescence lifetime distributions of diphenylhexatriene-labeled phosphatidylcholine as a tool for the study of phospholipid-cholesterol interactions.

Fluorescence lifetimes of 1-palmitoyl-2-diphenylhexatrienylpro-pionyl-phosphatidylc hol ine in vesicles of palmitoyloleoyl phosphatidylcholine (POPC) (1:300, mol/mol) in the liquid crystalline state were determined by multifrequency phase fluorometry. On the basis of statistic criteria (chi 2red) the measured phase angles and demodulation factors were equally well fitted to unimodal Lorentzian, Gaussian, or uniform lifetime distributions. No improvement in chi 2red could be observed if the experimental data were fitted to bimodal Lorentzian distributions or a double exponential decay. The unimodal Lorentzian lifetime distribution was characterized by a lifetime center of 6.87 ns and a full width at half maximum of 0.57 ns. Increasing amounts of cholesterol in the phospholipid vesicles (0-50 mol% relative to POPC) led to a slight increase of the lifetime center (7.58 ns at 50 mol% sterol) and reduced significantly the distributional width (0.14 ns at 50 mol% sterol). Lifetime distributions of POPC-cholesterol mixtures containing greater than 20 mol% sterol were within the resolution limit and could not be distinguished from monoexponential decays on the basis of chi 2red. Cholesterol stabilizes and rigidifies phospholipid bilayers in the fluid state. Considering its effect on lifetime distributions of fluorescent phospholipids it may also act as a membrane homogenizer.     

The effect of the width of the ventilation-perfusion distribution on arterial blood oxygen content

We investigate the effect of the width of ventilation-perfusion distributions on arterial blood oxygen content. We assume that the perfusion within the alveolar volume is a continuous function of ventilation-perfusion ratio, known as the continuous ventilation-perfusion distribution, and then write down the conservation of mass equations in the lung incorporating the nonlinear relationship between oxygen concentration in the gas phase and blood oxygen content. We solve these equations for various unimodal and bimodal ventilation-perfusion distributions believed to occur in practice and calculate the arterial blood oxygen content in each case. When a subject has a unimodal ventilation-perfusion distribution we show that the fraction of cardiac output to that mode (i.e. the fraction of non-shunted blood) has a large effect on arterial oxygen blood content. However, the width of the distribution has only a negligible effect on arterial oxygen blood content. For a bimodal ventilation-perfusion distribution the location and fraction of cardiac output to each mode has a large effect on arterial oxygen blood content. Again, the width of each mode of the distribution has little effect on arterial oxygen blood content. As a result there is little point, from a clinical perspective, in developing techniques for investigating the width of modes of these distributions since all relevant clinical information is contained in the nature (i.e. unimodal or bimodal) and in the location of the modes.     

Similarity of fluorescence lifetime distributions for single tryptophan proteins in the random coil state.

The picosecond time-resolved fluorescence decay data of nine single-tryptophan (trp) proteins and two multi-trp proteins in their native and denatured states were analyzed by the maximum entropy method (MEM). In the denatured state (6 M guanidine hydrochloride) a majority of the single-trp proteins show bimodal (at 25 degrees C) and trimodal (at 85 degrees C) distributions with similar patterns and similar values for average lifetimes. In the native state of the proteins the lifetime distributions were bimodal or trimodal. These results (multimodal distributions) are contradictory to the unimodal Lorentzian distribution of lifetimes reported for some proteins in the native and denatured states. MEM analysis gives a unimodal distribution of lifetimes only when the signal-to-noise ratio is poor in the time-resolved fluorescence decay data. The unimodal distribution model is therefore not realistic for proteins in the native and denatured states. The fluorescence decay components of the bi- or trimodal distribution are associated with the rotamer structures of the indole moiety when the protein is in the random coil state.     

Exploitation of GFP fusion proteins and stress avoidance as a generic strategy for the production of high-quality recombinant proteins

A C-terminal green fluorescent protein (GFP) fusion to a model target protein, Escherichia coli CheY, was exploited both as a reporter of the accumulation of soluble recombinant protein, and to develop a generic approach to optimize protein yields. The rapid accumulation of CheY∷GFP expressed from a pET20 vector under the control of an isopropyl-β- d -thiogalactoside (IPTG)-inducible T7 RNA polymerase resulted not only in the well-documented growth arrest but also loss of culturability and overgrowth of the productive population using plasmid-deficient bacteria. The highest yields of soluble CheY∷GFP as judged from the fluorescence levels were achieved using very low concentrations of IPTG, which avoid growth arrest and loss of culturability postinduction. Optimal product yields were obtained with 8 μM IPTG, a concentration so low that insufficient T7 RNA polymerase accumulated to be detectable by Western blot analysis. The improved protocol was shown to be suitable for process scale-up and intensification. It is also applicable to the accumulation of an untagged heterologous protein, cytochrome c₂ from Neisseria gonorrhoeae , which requires both secretion and extensive post-translational modification.     

lac operon induction in Escherichia coli: Systematic comparison of IPTG and TMG induction and influence of the transacetylase LacA

Most commonly used expression systems in bacteria are based on the Escherichia coli lac promoter. Furthermore, lac operon elements are used today in systems and synthetic biology. In the majority of the cases the gratuitous inducers IPTG or TMG are used. Here we report a systematic comparison of lac promoter induction by TMG and IPTG which focuses on the aspects inducer uptake, population heterogeneity and a potential influence of the transacetylase, LacA. We provide induction curves in E. coli LJ110 and in isogenic lacY and lacA mutant strains and we show that both inducers are substrates of the lactose permease at low inducer concentrations but can also enter cells independently of lactose permease if present at higher concentrations. Using a gfp reporter strain we compared TMG and IPTG induction at single cell level and showed that bimodal induction with IPTG occurred at approximately ten-fold lower concentrations than with TMG. Furthermore, we observed that lac operon induction is influenced by the transacetylase, LacA. By comparing two Plac-gfp reporter strains with and without a lacA deletion we could show that in the lacA⁺ strain the fluorescence level decreased after few hours while the fluorescence further increased in the lacA⁻ strain. The results indicate that through the activity of LacA the IPTG concentration can be reduced below an inducing threshold concentration—an influence that should be considered if low inducer amounts are used.     

Effect of lacY expression on homogeneity of induction from the P(tac) and P(trc) promoters by natural and synthetic inducers

The role of the Escherichia coli lactose permease (LacY) in the homogeneous induction of the lactose-inducible promoters P(tac) and P(trc) by the natural inducer lactose and the synthetic inducer isopropyl-beta-D-thiogalactopyranoside (IPTG) was investigated. Lactose requires active transport by LacY, whereas IPTG can freely penetrate the cell wall. In E. coli strains lacking a functional LacY, IPTG is required for induction of P(tac) and P(trc). In E. coli strains carrying a functional LacY, induction of P(trc) and P(tac) with intermediate concentrations of lactose gave rise to two subpopulations, one fully induced and one uninduced, whereas a single, fully induced population resulted when high inducer concentrations were used. In contrast, induction with IPTG gave rise to a single population of cells at all inducer concentrations in both lacY and lacY(+) strains.     

Variation in scale shape among alternative sympatric phenotypes of Arctic charr Salvelinus alpinus from two lakes in Scotland

Landmark‐based geometric morphometric analysis was used to detect differences in scale shape between ecologically distinct phenotypes of Arctic charr Salvelinus alpinus coexisting in the same lake. Relative warp analysis and standard multivariate analyses of the partial warps, obtained after a Procrustes superimposition, showed that scale landmarks were efficient in discriminating among two closely related alternative phenotypes within each of the two lakes. In Loch Tay, S. alpinus exhibited a bimodal body size‐frequency distribution among sexually mature fish, whereas in Loch Awe, S. alpinus are unimodal in body size but segregated into two distinct spawning phenotypes. In both lakes, alternative phenotypes showed significant differences in foraging ecology, habitat use and life history. It is probable that differences in scale shape reflect differences in ecology of these forms.     

Criteria for downhill protein folding: calorimetry, chevron plot, kinetic relaxation, and single-molecule radius of gyration in chain models with subdued degrees of cooperativity

Recent investigations of possible downhill folding of small proteins such as BBL have focused on the thermodynamics of non-two-state, "barrierless" folding/denaturation transitions. Downhill folding is noncooperative and thermodynamically "one-state," a phenomenon underpinned by a unimodal conformational distribution over chain properties such as enthalpy, hydrophobic exposure, and conformational dimension. In contrast, corresponding distributions for cooperative two-state folding are bimodal with well-separated population peaks. Using simplified atomic modeling of a three-helix bundle-in a scheme that accounts for hydrophobic interactions and hydrogen bonding-and coarse-grained C(alpha) models of four real proteins with various degrees of cooperativity, we evaluate the effectiveness of several observables at defining the underlying distribution. Bimodal distributions generally lead to sharper transitions, with a higher heat capacity peak at the transition midpoint, compared with unimodal distributions. However, the observation of a sigmoidal transition is not a reliable criterion for two-state behavior, and the heat capacity baselines, used to determine the van't Hoff and calorimetric enthalpies of the transition, can introduce ambiguity. Interestingly we find that, if the distribution of the single-molecule radius of gyration were available, it would permit discrimination between unimodal and bimodal underlying distributions. We investigate kinetic implications of thermodynamic noncooperativity using Langevin dynamics. Despite substantial chevron rollovers, the relaxation of the models considered is essentially single-exponential over an extended range of native stabilities. Consistent with experiments, significant deviations from single-exponential behavior occur only under strongly folding conditions.     

Pollination,mating and reproductive fitness in a plant population with bimodal floral‐tube length

Mating patterns and natural selection play important roles in determining whether genetic polymorphisms are maintained or lost. Here, we document an atypical population of Lapeirousia anceps (Iridaceae) with a bimodal distribution of floral‐tube length and investigate the reproductive mechanisms associated with this pattern of variation. Flowers were visited exclusively by the long‐proboscid fly Moegistorhynchus longirostris (Nemestrinidae), which exhibited a unimodal distribution of proboscis length and displayed a preference for long‐tubed phenotypes. Despite being visited by a single pollinator species, allozyme markers revealed significant genetic differentiation between open‐pollinated progeny of long‐ and short‐tubed phenotypes suggesting mating barriers between them. We obtained direct evidence for mating barriers between the floral‐tube phenotypes through observations of pollinator foraging, controlled hand pollinations and measurements of pollen competition and seed set. Intermediate tube‐length phenotypes produced fewer seeds in the field than either long‐ or short‐tubed phenotypes. Although floral‐tube length bimodality may not be a stable state over long timescales, reproductive barriers to mating and low ‘hybrid’ fitness have the potential to contribute to the maintenance of this state in the short term.     

Head shape dimorphism in European glass eels (Anguilla anguilla)

The life cycle of the European eel (Anguilla anguilla) remained a mystery until the 20th century, when Schmidt discovered that the Sargasso Sea was its spawning area. However, many aspects of the eel's life cycle remain poorly understood. Among these is the bimodal distribution in head shape, with broad- and narrowheaded phenotypes reported in the yellow eel stage. Although this has been linked to dietary preferences of the yellow eels, very little is known about why, how and when this dimorphism arises during their ontogeny. To determine whether this dimorphism indeed appears in relation to trophic niche segregation, we examined head shape variation at an earlier ontogenetic stage, the glass eel stage, as at this stage eels are considered to be non-feeding. Head shape was studied in a large dataset, containing glass eels captured from the Yser river mouth, the Leopold Canal (Belgium) and from the rivers Severn, Trent and Parret (UK), by both taking measurements (head width/head length) and using an outline analysis. Our results show that there is already considerable variation in broadness and bluntness of the head at the glass eel stage. In most cases, equal support for a unimodal and bimodal head shape distribution is found, whereas some cases support head shape bimodality in glass eels, suggesting that glass eel head shape might be shifting from a unimodal to a bimodal distribution. This, in combination with the observation that variation in head width/head length ratios in non-feeding glass eels shows a similar range as in feeding yellow eels, indicates that head shape in European eel might be at least partially determined through other mechanisms than trophic segregation.     

Characterization of the AlkS/P(alkB)-expression system as an efficient tool for the production of recombinant proteins in Escherichia coli fed-batch fermentations

The availability of suitable, well-characterized, and robust expression systems remains an essential requirement for successful metabolic engineering and recombinant protein production. We investigated the suitability of the Pseudomonas putida GPo1-derived AlkS/P(alkB) expression system in strictly aqueous cultures. By applying the apolar inducer dicyclopropylketone (DCPK) to express green fluorescent protein (GFP) from this system in Escherichia coli and analyzing the resulting cultures on single-cell level by flow cytometry, we found that this expression system gives rise to a homogeneous population of cells, even though the overall system is expected to have a positive feed-back element in the expression of the regulatory gene alkS. Overexpressing E. coli's serine hydroxymethyltransferase gene glyA, we showed that the system was already fully turned on at inducer concentrations as low as 0.005% (v/v). This allows efficient mass production of recombinant enzymes even though DCPK concentrations decreased from 0.05% to 0.01% over the course of a fully aerated cultivation in aqueous medium. Therefore, we elaborated the optimum induction procedure for production of the biocatalytically promising serine hydroxymethyltransferase and found volumetric and specific productivity to increase with specific growth rate in glucose-limited fed-batch cultures. Acetate excretion as a result of recombinant protein production could be avoided in an optimized fermentation protocol by switching earlier to a linear feed. This protocol resulted in a production of a final cell dry weight (CDW) concentration of 52 g/L, producing recombinant GlyA with a maximum specific activity of 6.3 U/mg total protein.     

Framework for online optimization of recombinant protein expression in high-cell-density Escherichia coli cultures using GFP-fusion monitoring

A framework for the online optimization of protein induction using green fluorescent protein (GFP)-monitoring technology was developed for high-cell-density cultivation of Escherichia coli. A simple and unstructured mathematical model was developed that described well the dynamics of cloned chloramphenicol acetyltransferase (CAT) production in E. coli JM105 was developed. A sequential quadratic programming (SQP) optimization algorithm was used to estimate model parameter values and to solve optimal open-loop control problems for piecewise control of inducer feed rates that maximize productivity. The optimal inducer feeding profile for an arabinose induction system was different from that of an isopropyl-beta-D-thiogalactopyranoside (IPTG) induction system. Also, model-based online parameter estimation and online optimization algorithms were developed to determine optimal inducer feeding rates for eventual use of a feedback signal from a GFP fluorescence probe (direct product monitoring with 95-minute time delay). Because the numerical algorithms required minimal processing time, the potential for product-based and model-based online optimal control methodology can be realized.     

Direct detection and quantification of horizontal gene transfer by using flow cytometry and gfp as a reporter gene

A new cultivation-independent method for studying conjugal gene transfer between bacteria was evaluated. The method was based on direct detection and enumeration of donor and transconjugant bacterial cells by flow cytometry. Specific detection of transconjugants was obtained by using a conjugative plasmid tagged with a reporter gene (gfp) encoding green fluorescent protein. A chromosomal encoded repressor (lacI(ql)) repressed expression of GFP in the donor bacteria. Enumeration of the donor cells was performed after induction of GFP expression by the addition of inducer isopropyl-thio-beta-D-galactoside (IPTG). The method presented here provided simple and precise quantification of horizontal gene transfer between both Escherichia coli and Pseudomonas putida strains.     

Automated live cell imaging of green fluorescent protein degradation in individual fibroblasts.

To accurately interpret the data from fluorescent proteins as reporters of gene activation within living cells, it is important to understand the kinetics of the degradation of the reporter proteins. We examined the degradation kinetics over a large number (>1,000) of single, living cells from a clonal population of NIH3T3 fibroblasts that were stably transfected with a destabilized, enhanced green fluorescent protein (eGFP) reporter driven by the tenascin-C promoter. Data collection and quantification of the fluorescence protein within a statistically significant number of individual cells over long times (14 h) by automated microscopy was facilitated by culturing cells on micropatterned arrays that confined their migration and allowed them to be segmented using phase contrast images. To measure GFP degradation rates unambiguously, protein synthesis was inhibited with cycloheximide. Results from automated live cell microscopy and image analysis indicated a wide range of cell-to-cell variability in the GFP fluorescence within individual cells. Degradation for this reporter was analyzed as a first order rate process with a degradation half-life of 2.8 h. We found that GFP degradation rates were independent of the initial intensity of GFP fluorescence within cells. This result indicates that higher GFP abundance in some cells is likely due to higher rates of gene expression, because it is not due to systematically lower rates of protein degradation. The approach described in this study will assist the quantification and understanding of gene activity within live cells using fluorescent protein reporters.