DNA Reaction System That Acquires Classical Conditioning.

Biochemical reaction networks can exhibit plastic adaptation to alter their functions in response to environmental changes. This capability is derived from the structure and dynamics of the reaction networks and the functionality of the biomolecule. This plastic adaptation in biochemical reaction systems is essentially related to memory and learning capabilities, which have been studied in DNA computing applications for the past decade. However, designing DNA reaction systems with memory and learning capabilities using the dynamic properties of biochemical reactions remains challenging. In this study, we propose a basic DNA reaction system design that acquires classical conditioning, a phenomenon underlying memory and learning, as a typical learning task. Our design is based on a simple mechanism of five DNA strand displacement reactions and two degradative reactions. The proposed DNA circuit can acquire or lose a new function under specific conditions, depending on the input history formed by repetitive stimuli, by exploiting the dynamic properties of biochemical reactions induced by different input timings.

Renewable DNA Proportional-Integral Controller with Photoresponsive Molecules.

A molecular robot is an intelligent molecular system. A typical control problem of molecular robots is to maintain the concentration of a specific DNA strand at the desired level, which is typically attained by a molecular feedback control mechanism. A molecular feedback system can be constructed in a bottom-up method by transforming a nonlinear chemical reaction system into a pseudo-linear system. This method enables the implementation of a molecular proportional-integral (PI) controller on a DNA reaction system. However, a DNA reaction system is driven by fuel DNA strand consumption, and without a sufficient amount of fuel strands, the molecular PI controller cannot perform normal operations as a concentration regulator. In this study, we developed a design method for a molecular PI control system to regenerate fuel strands by introducing photoresponsive reaction control. To this end, we employed a photoresponsive molecule, azobenzene, to guide the reaction direction forward or backward using light irradiation. We validated our renewable design of the PI controller by numerical simulations based on the reaction kinetics. We also confirmed the proof-of-principle of our renewable design by conducting experiments using a basic DNA circuit.

Individual hematopoietic stem cells in human bone marrow of patients with aplastic anemia or myelodysplastic syndrome stably give rise to limited cell lineages.

Mutation of the phosphatidylinositol N-acetylglucosaminyltransferase subunit A (PIG-A) gene in hematopoietic stem cells (HSCs) results in the loss of glycosylphosphatidylinositol-anchored proteins (GPI-APs) on HSCs, but minimally affects their development, and thus can be used as a clonal maker of HSCs. We analyzed GPI-APs expression on six major lineage cells in a total of 574 patients with bone marrow (BM) failure in which microenvironment itself is thought to be unaffected, including aplastic anemia (AA) or myelodysplastic syndrome (MDS). GPI-APs-deficient (GPI-APs(-) ) cells were detected in 250 patients. Whereas the GPI-APs(-) cells were seen in all six lineages in a majority of patients who had higher proportion ([dbmtequ]3%) of GPI-APs(-) cells, they were detected in only limited lineages in 92.9% of cases in the lower proportion (<3%) group. In all 250 cases, the same lineages of GPI-APs(-) cells were detected even after 6-18-month intervals, indicating that the GPI-APs(-) cells reflect hematopoiesis maintained by a self-renewing HSC in most of cases. The frequency of clones with limited lineages seen in mild cases of AA was similar to that in severe cases, and clones with limited lineages were seen even in two health volunteer cases. These results strongly suggest most individual HSCs produce only restricted lineages even in a steady state. While this restriction could reflect heterogeneity in the developmental potential of HSCs, we propose an alternative model in which the BM microenvironment is mosaic in supporting commitment of progenitors toward distinct lineages. Our computer simulation based on this model successfully recapitulated the observed clinical data.

Epidermal growth factor receptor mutation in combination with expression of MIG6 alters gefitinib sensitivity.

BACKGROUND: Epidermal growth factor receptor (EGFR) signaling plays an important role in the regulation of cell proliferation, survival, metastasis, and invasion in various tumors. Earlier studies showed that the EGFR is frequently overexpressed in non-small-cell lung cancer (NSCLC) and EGFR mutations at specific amino acid residues in the kinase domain induce altered responsiveness to gefitinib, a small molecule EGFR tyrosine kinase inhibitor. However, the mechanism underlying the drug response modulated by EGFR mutation is still largely unknown. To elucidate drug response in EGFR signal transduction pathway in which complex dynamics of multiple molecules involved, a systematic approach is necessary. In this paper, we performed experimental and computational analyses to clarify the underlying mechanism of EGFR signaling and cell-specific gefitinib responsiveness in three H1299-derived NSCLC cell lines; H1299 wild type (H1299WT), H1299 with an overexpressed wild type EGFR (H1299EGFR-WT), and H1299 with an overexpressed mutant EGFR L858R (H1299L858R; gefitinib sensitive mutant). RESULTS: We predicted and experimentally verified that Mig6, which is a known negative regulator of EGFR and specifically expressed in H1299L858R cells, synergized with gefitinib to suppress cellular growth. Computational analyses indicated that this inhibitory effect is amplified at the phosphorylation/dephosphorylation steps of MEK and ERK. CONCLUSIONS: Thus, we showed that L858R receptor mutation in combination with expression of its negative regulator, Mig6, alters signaling outcomes and results in variable drug sensitivity.

Computation of restoration of ligand response in the random kinetics of a prostate cancer cell signaling pathway.

Mutation and/or dysfunction of signaling proteins in the mitogen activated protein kinase (MAPK) signal transduction pathway are frequently observed in various kinds of human cancer. Consistent with this fact, in the present study, we experimentally observe that the epidermal growth factor (EGF) induced activation profile of MAP kinase signaling is not straightforward dose-dependent in the PC3 prostate cancer cells. To find out what parameters and reactions in the pathway are involved in this departure from the normal dose-dependency, a model-based pathway analysis is performed. The pathway is mathematically modeled with 28 rate equations yielding those many ordinary differential equations (ODE) with kinetic rate constants that have been reported to take random values in the existing literature. This has led to us treating the ODE model of the pathways kinetics as a random differential equations (RDE) system in which the parameters are random variables. We show that our RDE model captures the uncertainty in the kinetic rate constants as seen in the behavior of the experimental data and more importantly, upon simulation, exhibits the abnormal EGF dose-dependency of the activation profile of MAP kinase signaling in PC3 prostate cancer cells. The most likely set of values of the kinetic rate constants obtained from fitting the RDE model into the experimental data is then used in a direct transcription based dynamic optimization method for computing the changes needed in these kinetic rate constant values for the restoration of the normal EGF dose response. The last computation identifies the parameters, i.e., the kinetic rate constants in the RDE model, that are the most sensitive to the change in the EGF dose response behavior in the PC3 prostate cancer cells. The reactions in which these most sensitive parameters participate emerge as candidate drug targets on the signaling pathway.

Ligand-specific c-Fos expression emerges from the spatiotemporal control of ErbB network dynamics.

Activation of ErbB receptors by epidermal growth factor (EGF) or heregulin (HRG) determines distinct cell-fate decisions, although signals propagate through shared pathways. Using mathematical modeling and experimental approaches, we unravel how HRG and EGF generate distinct, all-or-none responses of the phosphorylated transcription factor c-Fos. In the cytosol, EGF induces transient and HRG induces sustained ERK activation. In the nucleus, however, ERK activity and c-fos mRNA expression are transient for both ligands. Knockdown of dual-specificity phosphatases extends HRG-stimulated nuclear ERK activation, but not c-fos mRNA expression, implying the existence of a HRG-induced repressor of c-fos transcription. Further experiments confirmed that this repressor is mainly induced by HRG, but not EGF, and requires new protein synthesis. We show how a spatially distributed, signaling-transcription cascade robustly discriminates between transient and sustained ERK activities at the c-Fos system level. The proposed control mechanisms are general and operate in different cell types, stimulated by various ligands.

Topological analysis of MAPK cascade for kinetic ErbB signaling.

Ligand-induced homo- and hetero-dimer formation of ErbB receptors results in different biological outcomes irrespective of recruitment and activation of similar effector proteins. Earlier experimental research indicated that cells expressing both EGFR (epidermal growth factor receptor) and the ErbB4 receptor (E1/4 cells) induced E1/4 cell-specific B-Raf activation and higher extracellular signal-regulated kinase (ERK) activation, followed by cellular transformation, than cells solely expressing EGFR (E1 cells) in Chinese hamster ovary (CHO) cells. Since our experimental data revealed the presence of positive feedback by ERK on upstream pathways, it was estimated that the cross-talk/feedback pathway structure of the Raf-MEK-ERK cascade might affect ERK activation dynamics in our cell system. To uncover the regulatory mechanism concerning the ERK dynamics, we used topological models and performed parameter estimation for all candidate structures that possessed ERK-mediated positive feedback regulation of Raf. The structure that reliably reproduced a series of experimental data regarding signal amplitude and duration of the signaling molecules was selected as a solution. We found that the pathway structure is characterized by ERK-mediated positive feedback regulation of B-Raf and B-Raf-mediated negative regulation of Raf-1. Steady-state analysis of the estimated structure indicated that the amplitude of Ras activity might critically affect ERK activity through ERK-B-Raf positive feedback coordination with sustained B-Raf activation in E1/4 cells. However, Rap1 that positively regulates B-Raf activity might be less effective concerning ERK and B-Raf activity. Furthermore, we investigated how such Ras activity in E1/4 cells can be regulated by EGFR/ErbB4 heterodimer-mediated signaling. From a sensitivity analysis of the detailed upstream model for Ras activation, we concluded that Ras activation dynamics is dominated by heterodimer-mediated signaling coordination with a large initial speed of dimerization when the concentration of the ErbB4 receptor is considerably high. Such characteristics of the signaling cause the preferential binding of the Grb2-SOS complex to heterodimer-mediated signaling molecules.

Quantitative transcriptional control of ErbB receptor signaling undergoes graded to biphasic response for cell differentiation.

ErbB receptor ligands, epidermal growth factor (EGF) and heregulin (HRG), induce dose-dependent transient and sustained intracellular signaling, proliferation, and differentiation of MCF-7 breast cancer cells, respectively. In an effort to delineate the ligand-specific cell determination mechanism, we investigated time course gene expressions induced by EGF and HRG that induce distinct cellular phenotypes in MCF-7 cells. To analyze independently the effects of ligand dosage and time for gene expression, we developed a statistical method for estimating the two effects. Our results indicated that signal transduction pathways convey quantitative properties of the dose-dependent activation of ErbB receptor to early transcription. The results also implied that moderate changes in the expression levels of a number of genes, not the predominant regulation of a few specific genes, might cooperatively work at the early stage of the transcription for determining cell fate. However, the EGF- and HRG-induced distinct signal durations resulted in the ligand-oriented biphasic induction of proteins after 20 min. The selected gene list and HRG-induced prolonged signaling suggested that transcriptional feedback to the intracellular signaling results in a graded to biphasic response in the cell determination process and that each ErbB receptor is inextricably responsible for the control of amplitude and duration of cellular biochemical reactions.