A hybrid in silico/in-cell controller for microbial bioprocesses with process-model mismatch

Bioprocess optimization using mathematical models is prevalent, yet the discrepancy between model predictions and actual processes, known as process-model mismatch (PMM), remains a significant challenge. This study proposes a novel hybrid control system called the hybrid in silico/in-cell controller (HISICC) to address PMM by combining model-based optimization (in silico feedforward controller) with feedback controllers utilizing synthetic genetic circuits integrated into cells (in-cell feedback controller). We demonstrated the efficacy of HISICC using two engineered Escherichia coli strains, TA1415 and TA2445, previously developed for isopropanol (IPA) production. TA1415 contains a metabolic toggle switch (MTS) to manage the competition between cell growth and IPA production for intracellular acetyl-CoA by responding to external input of isopropyl β-D-1-thiogalactopyranoside (IPTG). TA2445, in addition to the MTS, has a genetic circuit that detects cell density to autonomously activate MTS. The combination of TA2445 with an in silico controller exemplifies HISICC implementation. We constructed mathematical models to optimize IPTG input values for both strains based on the two-compartment model and validated these models using experimental data of the IPA production process. Using these models, we evaluated the robustness of HISICC against PMM by comparing IPA yields with two strains in simulations assuming various magnitudes of PMM in cell growth rates. The results indicate that the in-cell feedback controller in TA2445 effectively compensates for PMM by modifying MTS activation timing. In conclusion, the HISICC system presents a promising solution to the PMM problem in bioprocess engineering, paving the way for more efficient and reliable optimization of microbial bioprocesses.

Decoding cellular deformation from pseudo-simultaneously observed Rho GTPase activities

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.

Single-Cell Information Analysis Reveals That Skeletal Muscles Incorporate Cell-to-Cell Variability as Information Not Noise

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.

Trans-omic Analysis Reveals ROS-Dependent Pentose Phosphate Pathway Activation after High-Frequency Electrical Stimulation in C2C12 Myotubes

Skeletal muscle adaptation is mediated by cooperative regulation of metabolism, signal transduction, and gene expression. However, the global regulatory mechanism remains unclear. To address this issue, we performed electrical pulse stimulation (EPS) in differentiated C2C12 myotubes at low and high frequency, carried out metabolome and transcriptome analyses, and investigated phosphorylation status of signaling molecules. EPS triggered extensive and specific changes in metabolites, signaling phosphorylation, and gene expression during and after EPS in a frequency-dependent manner. We constructed trans-omic network by integrating these data and found selective activation of the pentose phosphate pathway including metabolites, upstream signaling molecules, and gene expression of metabolic enzymes after high-frequency EPS. We experimentally validated that activation of these molecules after high-frequency EPS was dependent on reactive oxygen species (ROS). Thus, the trans-omic analysis revealed ROS-dependent activation in signal transduction, metabolome, and transcriptome after high-frequency EPS in C2C12 myotubes, shedding light on possible mechanisms of muscle adaptation.

•

We performed electrical pulse stimulation in differentiated C2C12 myotubes

•

We constructed trans-omic network after high-frequency electrical pulse stimulation

•

Trans-omic network integrates metabolome, transcriptome, and signaling molecules

•

We identified ROS-dependent pentose phosphate pathway activation

Automatic Quantitative Segmentation of Myotubes Reveals Single-cell Dynamics of S6 Kinase Activation

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.

System identification of signaling dependent gene expression with different time-scale data

Cells decode information of signaling activation at a scale of tens of minutes by downstream gene expression with a scale of hours to days, leading to cell fate decisions such as cell differentiation. However, no system identification method with such different time scales exists. Here we used compressed sensing technology and developed a system identification method using data of different time scales by recovering signals of missing time points. We measured phosphorylation of ERK and CREB, immediate early gene expression products, and mRNAs of decoder genes for neurite elongation in PC12 cell differentiation and performed system identification, revealing the input-output relationships between signaling and gene expression with sensitivity such as graded or switch-like response and with time delay and gain, representing signal transfer efficiency. We predicted and validated the identified system using pharmacological perturbation. Thus, we provide a versatile method for system identification using data with different time scales.

Selective control of up-regulated and down-regulated genes by temporal patterns and doses of insulin

Secretion of insulin transiently increases after eating, resulting in a high circulating concentration. Fasting limits insulin secretion, resulting in a low concentration of insulin in the circulation. We analyzed transcriptional responses to different temporal patterns and doses of insulin in the hepatoma FAO cells and identified 13 up-regulated and 16 down-regulated insulin-responsive genes (IRGs). The up-regulated IRGs responded more rapidly than did the down-regulated IRGs to transient stepwise or pulsatile increases in insulin concentration, whereas the down-regulated IRGs were repressed at lower concentrations of insulin than those required to stimulate the up-regulated IRGs. Mathematical modeling of the insulin response as two stages-(i) insulin signaling to transcription and (ii)transcription and mRNA stability-indicated that the first stage was the more rapid stage for the down-regulated IRGs, whereas the second stage of transcription was the more rapid stage for the up-regulated IRGs. A subset of the IRGs that were up-regulated or down-regulated in the FAO cells was similarly regulated in the livers of rats injected with a single dose of insulin. Thus, not only can cells respond to insulin but they can also interpret the intensity and pattern of signal to produce distinct transcriptional responses. These results provide insight that may be useful in treating obesity and type 2 diabetes associated with aberrant insulin production or tissue responsiveness.

Robustness of MEK-ERK Dynamics and Origins of Cell-to-Cell Variability in MAPK Signaling

Cellular signaling processes can exhibit pronounced cell-to-cell variability in genetically identical cells. This affects how individual cells respond differentially to the same environmental stimulus. However, the origins of cell-to-cell variability in cellular signaling systems remain poorly understood. Here, we measure the dynamics of phosphorylated MEK and ERK across cell populations and quantify the levels of population heterogeneity over time using high-throughput image cytometry. We use a statistical modeling framework to show that extrinsic noise, particularly that from upstream MEK, is the dominant factor causing cell-to-cell variability in ERK phosphorylation, rather than stochasticity in the phosphorylation/dephosphorylation of ERK. We furthermore show that without extrinsic noise in the core module, variable (including noisy) signals would be faithfully reproduced downstream, but the within-module extrinsic variability distorts these signals and leads to a drastic reduction in the mutual information between incoming signal and ERK activity.

•

Active MEK and ERK levels differ profoundly among genetically identical cells

•

A statistical framework is developed to identify the causes of this variability

•

Analysis shows that extrinsic noise upstream MEK-ERK module causes cell variability

•

Within-module extrinsic variability distorts signals

Reconstruction of insulin signal flow from phosphoproteome and metabolome data

Cellular homeostasis is regulated by signals through multiple molecular networks that include protein phosphorylation and metabolites. However, where and when the signal flows through a network and regulates homeostasis has not been explored. We have developed a reconstruction method for the signal flow based on time-course phosphoproteome and metabolome data, using multiple databases, and have applied it to acute action of insulin, an important hormone for metabolic homeostasis. An insulin signal flows through a network, through signaling pathways that involve 13 protein kinases, 26 phosphorylated metabolic enzymes, and 35 allosteric effectors, resulting in quantitative changes in 44 metabolites. Analysis of the network reveals that insulin induces phosphorylation and activation of liver-type phosphofructokinase 1, thereby controlling a key reaction in glycolysis. We thus provide a versatile method of reconstruction of signal flow through the network using phosphoproteome and metabolome data.

FRET imaging and statistical signal processing reveal positive and negative feedback loops regulating the morphology of randomly migrating HT-1080 cells

Cell migration plays an important role in many physiological processes. Rho GTPases (Rac1, Cdc42, RhoA) and phosphatidylinositols have been extensively studied in directional cell migration. However, it remains unclear how Rho GTPases and phosphatidylinositols regulate random cell migration in space and time. We have attempted to address this issue using fluorescence resonance energy transfer (FRET) imaging and statistical signal processing. First, we acquired time-lapse images of random migration of HT-1080 fibrosarcoma cells expressing FRET biosensors of Rho GTPases and phosphatidyl inositols. We developed an image-processing algorithm to extract FRET values and velocities at the leading edge of migrating cells. Auto- and cross-correlation analysis suggested the involvement of feedback regulations among Rac1, phosphatidyl inositols and membrane protrusions. To verify the feedback regulations, we employed an acute inhibition of the signaling pathway with pharmaceutical inhibitors. The inhibition of actin polymerization decreased Rac1 activity, indicating the presence of positive feedback from actin polymerization to Rac1. Furthermore, treatment with PI3-kinase inhibitor induced an adaptation of Rac1 activity, i.e. a transient reduction of Rac1 activity followed by recovery to the basal level. In silico modeling that reproduced the adaptation predicted the existence of a negative feedback loop from Rac1 to actin polymerization. Finally, we identified MLCK as the probable controlling factor in the negative feedback. These findings quantitatively demonstrate positive and negative feedback loops that involve actin, Rac1 and MLCK, and account for the ordered patterns of membrane dynamics observed in randomly migrating cells.

SH3BP1, an exocyst-associated RhoGAP, inactivates Rac1 at the front to drive cell motility

The coordination of the several pathways involved in cell motility is poorly understood. Here, we identify SH3BP1, belonging to the RhoGAP family, as a partner of the exocyst complex and establish a physical and functional link between two motility-driving pathways, the Ral/exocyst and Rac signaling pathways. We show that SH3BP1 localizes together with the exocyst to the leading edge of motile cells and that SH3BP1 regulates cell migration via its GAP activity upon Rac1. SH3BP1 loss of function induces abnormally high Rac1 activity at the front, as visualized by in vivo biosensors, and disorganized and instable protrusions, as revealed by cell morphodynamics analysis. Consistently, constitutively active Rac1 mimics the phenotype of SH3BP1 depletion: slow migration and aberrant cell morphodynamics. Our finding that SH3BP1 downregulates Rac1 at the motile-cell front indicates that Rac1 inactivation in this location, as well as its activation by GEF proteins, is a fundamental requirement for cell motility.