Analysis of Cytokine Profiles in Pediatric Myocarditis Multicenter Study.

Acute myocarditis (AM) is an inflammatory disease of the heart muscle that can progress to fulminant myocarditis (FM), a severe and life-threatening condition. The cytokine profile of myocarditis in children, especially in relation to fulminant myocarditis, is not well understood. This study aims to evaluate the cytokine profiles of acute and fulminant myocarditis in children. Pediatric patients diagnosed with myocarditis were included in the study. Cytokine levels were measured using a multiplexed fluorescent bead-based immunoassay. Statistical analysis was performed to compare patient characteristics and cytokine levels between FM, AM, and healthy control (HC) groups. Principal component analysis (PCA) was applied to cytokine groups that were independent among the FM, AM, and HC groups. The study included 22 patients with FM and 14 with AM patients. We identified four cytokines that were significantly higher in the FM group compared to the AM group: IL1-RA (p = 0.002), IL-8 (p = 0.005), IL-10 (p = 0.011), and IL-15 (p = 0.005). IL-4 was significantly higher in the AM group compared to FM and HC groups (p = 0.006 and 0.0015). PDGF-AA, and VEGF-A were significantly lower in the FM group than in the AM group (p = 0.013 and <0.001). Similar results were obtained in PCA. Cytokine profiles might be used to differentiate pediatric FM from AM, stratify severity, and predict prognosis. The targeted therapy that works individual cytokines might provide a potential treatment for reducing the onset of the FM and calming the condition, and further studies are needed.

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.

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.

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.

Processive phosphorylation of ERK MAP kinase in mammalian cells.

The mitogen-activated protein (MAP) kinase pathway is comprised of a three-tiered kinase cascade. The distributive kinetic mechanism of two-site MAP kinase phosphorylation inherently generates a nonlinear switch-like response. However, a linear graded response of MAP kinase has also been observed in mammalian cells, and its molecular mechanism remains unclear. To dissect these input-output behaviors, we quantitatively measured the kinetic parameters involved in the MEK (MAPK/ERK kinase)-ERK MAP kinase signaling module in HeLa cells. Using a numerical analysis based on experimentally determined parameters, we predicted in silico and validated in vivo that ERK is processively phosphorylated in HeLa cells. Finally, we identified molecular crowding as a critical factor that converts distributive phosphorylation into processive phosphorylation. We proposed the term quasi-processive phosphorylation to describe this mode of ERK phosphorylation that is operated under the physiological condition of molecular crowding. The generality of this phenomenon may provide a new paradigm for a diverse set of biochemical reactions including multiple posttranslational modifications.

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.

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.

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.

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.

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.

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.

Distinct predictive performance of Rac1 and Cdc42 in cell migration.

We propose a new computation-based approach for elucidating how signaling molecules are decoded in cell migration. In this approach, we performed FRET time-lapse imaging of Rac1 and Cdc42, members of Rho GTPases which are responsible for cell motility, and quantitatively identified the response functions that describe the conversion from the molecular activities to the morphological changes. Based on the identified response functions, we clarified the profiles of how the morphology spatiotemporally changes in response to local and transient activation of Rac1 and Cdc42, and found that Rac1 and Cdc42 activation triggers laterally propagating membrane protrusion. The response functions were also endowed with property of differentiator, which is beneficial for maintaining sensitivity under adaptation to the mean level of input. Using the response function, we could predict the morphological change from molecular activity, and its predictive performance provides a new quantitative measure of how much the Rho GTPases participate in the cell migration. Interestingly, we discovered distinct predictive performance of Rac1 and Cdc42 depending on the migration modes, indicating that Rac1 and Cdc42 contribute to persistent and random migration, respectively. Thus, our proposed predictive approach enabled us to uncover the hidden information processing rules of Rho GTPases in the cell migration.

An optogenetic system for interrogating the temporal dynamics of Akt.

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.

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.