Katsuyuki KUNIDA

The optimization of bioprocess inputs using mathematical models is widely practiced. However, the mismatch between model prediction and the actual process [called process-model mismatch (PMM)] is problematic; when a large PMM exists, the process inputs optimized using the mathematical model in advance are no longer optimal for the actual process. In this study, we propose a hybrid control system that combines model-based optimization (in silico feedforward controller) and feedback controllers using synthetic genetic circuits integrated into cells (in-cell feedback controller) - which we named the hybrid in silico/in-cell controller (HISICC) - as a solution to this PMM issue. As a proof of concept for HISICC, we constructed a mathematical model of an engineered Escherichia coli strain for the isopropanol production process that was previously developed. This strain contains an in-cell feedback controller, and its combination with an in silico controller can be regarded as an example of HISICC. We demonstrated the robustness of HISICC against PMM by comparing the strain with another strain with no in-cell feedback controller in simulations assuming PMM of various magnitudes.

Limitations in simultaneously observing the activity of multiple molecules in live cells prevent researchers from elucidating how these molecules coordinate the dynamic regulation of cellular functions. Here, we propose the motion-triggered average (MTA) algorithm to characterize pseudo-simultaneous dynamic changes in arbitrary cellular deformation and molecular activities. Using MTA, we successfully extract a pseudo-simultaneous time series from individually observed activities of three Rho GTPases: Cdc42, Rac1, and RhoA. To verify that this time series encoded information on cell-edge movement, we use a mathematical regression model to predict the edge velocity from the activities of the three molecules. The model accurately predicts the unknown edge velocity, providing numerical evidence that these Rho GTPases regulate edge movement. Data preprocessing using MTA combined with mathematical regression provides an effective strategy for reusing numerous individual observations of molecular activities.

Cell-to-cell variability in signal transduction in biological systems is often considered noise. However, intercellular variation (i.e., cell-to-cell variability) has the potential to enable individual cells to encode different information. Here, we show that intercellular variation increases information transmission of skeletal muscle. We analyze the responses of multiple cultured myotubes or isolated skeletal muscle fibers as a multiple-cell channel composed of single-cell channels. We find that the multiple-cell channel, which incorporates intercellular variation as information, not noise, transmitted more information in the presence of intercellular variation than in the absence according to the "response diversity effect," increasing in the gradualness of dose response by summing the cell-to-cell variable dose responses. We quantify the information transmission of human facial muscle contraction during intraoperative neurophysiological monitoring and find that information transmission of muscle contraction is comparable to that of a multiple-cell channel. Thus, our data indicate that intercellular variation can increase the information capacity of tissues.

Automatic cell segmentation is a powerful method for quantifying signaling dynamics at single-cell resolution in live cell fluorescence imaging. Segmentation methods for mononuclear and round shape cells have been developed extensively. However, a segmentation method for elongated polynuclear cells, such as differentiated C2C12 myotubes, has yet to be developed. In addition, myotubes are surrounded by undifferentiated reserve cells, making it difficult to identify background regions and subsequent quantification. Here we developed an automatic quantitative segmentation method for myotubes using watershed segmentation of summed binary images and a two-component Gaussian mixture model. We used time-lapse fluorescence images of differentiated C2C12 cells stably expressing Eevee-S6K, a fluorescence resonance energy transfer (FRET) biosensor of S6 kinase (S6K). Summation of binary images enhanced the contrast between myotubes and reserve cells, permitting detection of a myotube and a myotube center. Using a myotube center instead of a nucleus, individual myotubes could be detected automatically by watershed segmentation. In addition, a background correction using the two-component Gaussian mixture model permitted automatic signal intensity quantification in individual myotubes. Thus, we provide an automatic quantitative segmentation method by combining automatic myotube detection and background correction. Furthermore, this method allowed us to quantify S6K activity in individual myotubes, demonstrating that some of the temporal properties of S6K activity such as peak time and half-life of adaptation show different dose-dependent changes of insulin between cell population and individuals.Key words: time lapse images, cell segmentation, fluorescence resonance energy transfer, C2C12, myotube.

The dynamic activity of the serine/threonine kinase Akt is crucial for the regulation of diverse cellular functions, but the precise spatiotemporal control of its activity remains a critical issue. Herein, we present a photo-activatable Akt (PA-Akt) system based on a light-inducible protein interaction module of Arabidopsis thaliana cryptochrome2 (CRY2) and CIB1. Akt fused to CRY2phr, which is a minimal light sensitive domain of CRY2 (CRY2-Akt), is reversibly activated by light illumination in several minutes within a physiological dynamic range and specifically regulates downstream molecules and inducible biological functions. We have generated a computational model of CRY2-Akt activation that allows us to use PA-Akt to control the activity quantitatively. The system provides evidence that the temporal patterns of Akt activity are crucial for generating one of the downstream functions of the Akt-FoxO pathway; the expression of a key gene involved in muscle atrophy (Atrogin-1). The use of an optical module with computational modeling represents a general framework for interrogating the temporal dynamics of biomolecules by predictive manipulation of optogenetic modules.