Bei-Wen Ying

Bimodal protein solubility distribution revealed by an aggregation analysis of the entire ensemble of Escherichia coli proteins

Protein folding often competes with intermolecular aggregation, which in most cases irreversibly impairs protein function, as exemplified by the formation of inclusion bodies. Although it has been empirically determined that some proteins tend to aggregate, the relationship between the protein aggregation propensities and the primary sequences remains poorly understood. Here, we individually synthesized the entire ensemble of Escherichia coli proteins by using an in vitro reconstituted translation system and analyzed the aggregation propensities. Because the reconstituted translation system is chaperone-free, we could evaluate the inherent aggregation propensities of thousands of proteins in a translation-coupled manner. A histogram of the solubilities, based on data from 3,173 translated proteins, revealed a clear bimodal distribution, indicating that the aggregation propensities are not evenly distributed across a continuum. Instead, the proteins can be categorized into 2 groups, soluble and aggregation-prone proteins. The aggregation propensity is most prominently correlated with the structural classification of proteins, implying that the prediction of aggregation propensity requires structural information about the protein.

Comparison of sequence reads obtained from three next-generation sequencing platforms.

Next-generation sequencing technologies enable the rapid cost-effective production of sequence data. To evaluate the performance of these sequencing technologies, investigation of the quality of sequence reads obtained from these methods is important. In this study, we analyzed the quality of sequence reads and SNP detection performance using three commercially available next-generation sequencers, i.e., Roche Genome Sequencer FLX System (FLX), Illumina Genome Analyzer (GA), and Applied Biosystems SOLiD system (SOLiD). A common genomic DNA sample obtained from Escherichia coli strain DH1 was applied to these sequencers. The obtained sequence reads were aligned to the complete genome sequence of E. coli DH1, to evaluate the accuracy and sequence bias of these sequence methods. We found that the fraction of “junk” data, which could not be aligned to the reference genome, was largest in the data set of SOLiD, in which about half of reads could not be aligned. Among data sets after alignment to the reference, sequence accuracy was poorest in GA data sets, suggesting relatively low fidelity of the elongation reaction in the GA method. Furthermore, by aligning the sequence reads to the E. coli strain W3110, we screened sequence differences between two E. coli strains using data sets of three different next-generation platforms. The results revealed that the detected sequence differences were similar among these three methods, while the sequence coverage required for the detection was significantly small in the FLX data set. These results provided valuable information on the quality of short sequence reads and the performance of SNP detection in three next-generation sequencing platforms.

Noisy cell growth rate leads to fluctuating protein concentration in bacteria

The present study discusses a prime cause of fluctuating protein concentrations, which play a significant role in generating phenotypic diversity in bacteria. A genetic circuit integrated in a bacterial genome was used to evaluate the cell-to-cell variation in protein concentration. A simple dynamic model, comprising terms for synthesis and dilution, was used to elucidate the contributions of distinct noises to the fluctuation in cell protein concentration. Experimental and theoretical results demonstrated that noise in the rate of increase in cell volume (cell growth rate) serves as a source of extrinsic noise that accounts for dozens of percent of the total noise, whereas intrinsic noise in protein synthesis makes only a moderate contribution to the fluctuation in protein concentration. This suggests that such external noise in the cell growth rate has a global effect on cellular components, resulting in a large fluctuation in protein concentration in bacterial cells.

Construction of Escherichia coli gene expression level perturbation collection.

We generated 61 strains of Escherichia coli in which the expression level of a specific single gene can be changed continuously over a physiologically significant range. In each strain, one auxotrophic gene was deleted from its original position and reinserted at a specific position on the chromosome under the control of the tetA promoter. Therefore, the level of expression of the target gene can be controlled easily by altering the concentrations of inducers, e.g., anhydrotetracycline and doxycycline, in the medium. Protein and mRNA levels and changes in proliferation rate were examined in some of the strains in our collection to determine the ability to control the level of target gene expression over a physiologically significant range. These strains will be useful for extracting omics data sets and for the construction of genome-scale mathematical models, because causality between perturbations in gene expression level and their consequences can be clearly determined.

Co-translational Binding of GroEL to Nascent Polypeptides Is Followed by Post-translational Encapsulation by GroES to Mediate Protein Folding *

The eubacterial chaperonins GroEL and GroES are essential chaperones and primarily assist protein folding in the cell. Although the molecular mechanism of the GroEL system has been examined previously, the mechanism by which GroEL and GroES assist folding of nascent polypeptides during translation is still poorly understood. We previously demonstrated a co-translational involvement of the Escherichia coli GroEL in folding of newly synthesized polypeptides using a reconstituted cell-free translation system (Ying, B. W., Taguchi, H., Kondo, M., and Ueda, T. (2005) J. Biol. Chem. 280, 12035–12040). Employing the same system here, we further characterized the mechanism by which GroEL assists folding of translated proteins via encapsulation into the GroEL-GroES cavity. The stable co-translational association between GroEL and the newly synthesized polypeptide is dependent on the length of the nascent chain. Furthermore, GroES is capable of interacting with the GroEL-nascent peptide-ribosome complex, and experiments using a single-ring variant of GroEL clearly indicate that GroES association occurs only at the trans-ring, not the cis-ring, of GroEL. GroEL holds the nascent chain on the ribosome in a polypeptide length-dependent manner and post-translationally encapsulates the polypeptide using the GroES cap to accomplish the chaperonin-mediated folding process.

Multilevel comparative analysis of the contributions of genome reduction and heat shock to the Escherichia coli transcriptome

BackgroundBoth large deletions in genome and heat shock stress would lead to alterations in the gene expression profile; however, whether there is any potential linkage between these disturbances to the transcriptome have not been discovered. Here, the relationship between the genomic and environmental contributions to the transcriptome was analyzed by comparing the transcriptomes of the bacterium Escherichia coli (strain MG1655 and its extensive genomic deletion derivative, MDS42) grown in regular and transient heat shock conditions.ResultsThe transcriptome analysis showed the following: (i) there was a reorganization of the transcriptome in accordance with preferred chromosomal periodicity upon genomic or heat shock perturbation; (ii) there was a considerable overlap between the perturbed regulatory networks and the categories enriched for differentially expressed genes (DEGs) following genome reduction and heat shock; (iii) the genes sensitive to genome reduction tended to be located close to genomic scars, and some were also highly responsive to heat shock; and (iv) the genomic and environmental contributions to the transcriptome displayed not only a positive correlation but also a negatively compensated relationship (i.e., antagonistic epistasis).ConclusionThe contributions of genome reduction and heat shock to the Escherichia coli transcriptome were evaluated at multiple levels. The observations of overlapping perturbed networks, directional similarity in transcriptional changes, positive correlation and epistatic nature linked the two contributions and suggest somehow a crosstalk guiding transcriptional reorganization in response to both genetic and environmental disturbances in bacterium E. coli.

Adaptation by stochastic switching of a monostable genetic circuit in Escherichia coli

Stochastic switching is considered as a cost‐saving strategy for adaptation to environmental challenges. We show here that stochastic switching of a monostable circuit can mediate the adaptation of the engineered OSU12‐hisC Escherichia coli strain to histidine starvation. In this strain, the hisC gene was deleted from the His operon and placed under the control of a monostable foreign promoter. In response to histidine depletion, the OSU12‐hisC population shifted to a higher HisC expression level, which is beneficial under starving conditions but is not favoured by the monostable circuit. The population shift was accompanied by growth recovery and was reversible upon histidine addition. A weak directionality in stochastic switching of hisC was observed in growing microcolonies under histidine‐free conditions. Directionality and fate decision were in part dependent on the initial cellular status. Finally, microarray analysis indicated that OSU12‐hisC reorganized its transcriptome to reach the appropriate physiological state upon starvation. These findings suggest that bacteria do not necessarily need to evolve signalling mechanisms to control gene expression appropriately, even for essential genes.

Stochastic Switching Induced Adaptation in a Starved Escherichia coli Population

Population adaptation can be determined by stochastic switching in living cells. To examine how stochastic switching contributes to the fate decision for a population under severe stress, we constructed an Escherichia coli strain crucially dependent on the expression of a rewired gene. The gene essential for tryptophan biosynthesis, trpC, was removed from the native regulatory unit, the Trp operon, and placed under the extraneous control of the lactose utilisation network. Bistability of the network provided the cells two discrete phenotypes: the induced and suppressed level of trpC. The two phenotypes permitted the cells to grow or not, respectively, under conditions of tryptophan depletion. We found that stochastic switching between the two states allowed the initially suppressed cells to form a new population with induced trpC in response to tryptophan starvation. However, the frequency of the transition from suppressed to induced state dropped off dramatically in the starved population, in comparison to that in the nourished population. This reduced switching rate was compensated by increasing the initial population size, which probably provided the cell population more chances to wait for the rarely appearing fit cells from the unfit cells. Taken together, adaptation of a starved bacterial population because of stochasticity in the gene rewired from the ancient regulon was experimentally confirmed, and the nutritional status and the population size played a great role in stochastic adaptation.

Refined method for the genomic integration of complex synthetic circuits.

Genetic reconstruction of regulatory gene circuits is currently applied in systematic dynamics and structure-function studies of intact cellular networks in systems biology. We present a modified procedure for the integration of complex genetic circuits into the Escherichia coli genome, to provide an efficient synthetic approach for stochastic study and the artificial engineering of genetic networks. Linear artificial sequences of various lengths were easily integrated into the bacterial genome at one time. Comparison of the cellular concentrations of proteins encoded by genes carried on plasmids or the genome indicated that genome recombination could minimize the copy number noise in the genetic circuit, allowing precise design and interpretation of the cellular network. The refined recombination procedure allowed efficient construction of a single copy of a complex genetic circuit in cells, and the resultant reduced fluctuation in copy number led to accurate phenotypic behaviour of the genome-integrated synthetic switch corresponding to the design principle. The availability of long-fragment insertions makes the reconstruction of complex networks easy on the genome, and provides a powerful tool for precise engineering in synthetic and systems biology.

Bacterial Cells Carrying Synthetic Dual-Function Operon Survived Starvation

A synthetic dual-function operon with a bistable structure was designed and successfully integrated into the bacterial genome. Bistability was generated by the mutual inhibitory structure comprised of the promoters Ptet and Plac and the repressors LacI and TetR. Dual function essential for cell growth was introduced by replacing the genes (i.e., hisC and leuB) encoding proteins involved in the biosynthesis of histidine and leucine from their native chromosomal locations to the synthetic operon. Both colony formation and population dynamics of the cells carrying this operon showed that the cells survived starvation and the newly formed population transited between the two stable states, representing the induced hisC and leuB levels, in accordance with the nutritional status. The results strongly suggested that the synthetic design of proto-operons sensitive to external perturbations is practical and functional in native cells.