Browsing by Author "Huang, Bin"

Now showing 1 - 8 of 8
  • Thumbnail Image
    Item
    Decoding biological gene regulatory networks by quantitative modeling
    (2017-04-20) Huang, Bin; Onuchic, Jose
    Gene regulatory network is essential to regulate the biological functions of cells. With the rapid development of “omics” technologies, the network can be inferred for a certain biological function. However, it still remains a challenge to understand the complex network at a systematic level. In this thesis, we utilized quantitative modeling approaches to study the nonlinear dynamics and the design principles of these biological gene regulatory networks. We aim to explain the existing experimental observations with the model, and further propose reasonable hypothesis for future experimental designs. More importantly, the understanding of the circuit’s regulatory mechanism would benefit the design of a de novo gene circuit for a new biological function. We first studied the plasticity of cell migration phenotypes during cancer metastasis, which contains two key cellular plasticity mechanisms - epithelial-tomesenchymal transition (EMT) and mesenchymal-to-amoeboid transition (MAT). In this study, we quantitatively modeled the core Rac1/RhoA gene regulatory circuit for MAT and later connected it with the core regulatory circuit for EMT. We found four different stable states, consistent with the amoeboid (A), mesenchymal (M), the hybrid amoeboid/mesenchymal (A/M), and the hybrid epithelial/mesenchymal (E/M) phenotypes that are observed in the experiment. We also explored the effects of microRNAs and EMT-inducing signals like Hepatocyte Growth Factor (HGF), and provided a new insight for the transitions among these phenotypes. To improve the traditional modeling approaches, we developed a new computational modeling method called Random Circuit Perturbation (RACIPE) to explore the dynamic behavior of gene regulatory circuits without the requirement of detailed kinetic parameters. We applied RACIPE on several gene circuits, and found the existence of robust gene expression patterns even though the model parameters are wildly perturbed. We also showed the powerful aspect of RACIPE to decipher the operating principles of the circuits. This kind of quantitative models not only works for gene regulatory network, but also is capable to be extended to study the cell-cell interactions among cancer and immune cells. The results shown the co-occurrence of three cancer states: low risk cancer with intermediate immunity (L), intermediate risk cancer with high immunity (I) and high risk cancer with low immunity state (H). We further used the model to assess the different combinations of cancer therapies.
  • Thumbnail Image
    Item
    Distinguishing mechanisms underlying EMT tristability
    (Springer International Publishing, 2017) Jia, Dongya; Jolly, Mohit K.; Tripathi, Satyendra C.; Den Hollander, Petra; Huang, Bin; Lu, Mingyang; Celiktas, Muge; Ramirez-Peña, Esmeralda; Ben-Jacob, Eshel; Onuchic, José Nelson; Hanash, Samir M.; Mani, Sendurai A.; Levine, Herbert; Bioengineering; Biosciences; Chemistry; Physics and Astronomy
    Abstract Background The Epithelial-Mesenchymal Transition (EMT) endows epithelial-looking cells with enhanced migratory ability during embryonic development and tissue repair. EMT can also be co-opted by cancer cells to acquire metastatic potential and drug-resistance. Recent research has argued that epithelial (E) cells can undergo either a partial EMT to attain a hybrid epithelial/mesenchymal (E/M) phenotype that typically displays collective migration, or a complete EMT to adopt a mesenchymal (M) phenotype that shows individual migration. The core EMT regulatory network - miR-34/SNAIL/miR-200/ZEB1 - has been identified by various studies, but how this network regulates the transitions among the E, E/M, and M phenotypes remains controversial. Two major mathematical models – ternary chimera switch (TCS) and cascading bistable switches (CBS) - that both focus on the miR-34/SNAIL/miR-200/ZEB1 network, have been proposed to elucidate the EMT dynamics, but a detailed analysis of how well either or both of these two models can capture recent experimental observations about EMT dynamics remains to be done. Results Here, via an integrated experimental and theoretical approach, we first show that both these two models can be used to understand the two-step transition of EMT - E→E/M→M, the different responses of SNAIL and ZEB1 to exogenous TGF-β and the irreversibility of complete EMT. Next, we present new experimental results that tend to discriminate between these two models. We show that ZEB1 is present at intermediate levels in the hybrid E/M H1975 cells, and that in HMLE cells, overexpression of SNAIL is not sufficient to initiate EMT in the absence of ZEB1 and FOXC2. Conclusions These experimental results argue in favor of the TCS model proposing that miR-200/ZEB1 behaves as a three-way decision-making switch enabling transitions among the E, hybrid E/M and M phenotypes.
  • Thumbnail Image
    Item
    Implications of the hybrid epithelial/mesenchymal phenotype in metastasis
    (Frontiers Media S.A., 2015) Jolly, Mohit Kumar; Boareto, Marcelo; Huang, Bin; Jia, Dongya; Lu, Mingyang; Ben-Jacob, Eshel; Onuchic, José Nelson; Levine, Herbert; Bioengineering; Chemistry; Biosciences; Center for Theoretical Biological Physics; Systems, Synthetic, and Physical Biology
    Transitions between epithelial and mesenchymal phenotypes - the epithelial to -mesenchymal transition (EMT) and its reverse the mesenchymal to epithelial transition (MET) - are hallmarks of cancer metastasis. While transitioning between the epithelial and mesenchymal phenotypes, cells can also attain a hybrid epithelial/mesenchymal (E/M) (i.e., partial or intermediate EMT) phenotype. Cells in this phenotype have mixed epithelial (e.g., adhesion) and mesenchymal (e.g., migration) properties, thereby allowing them to move collectively as clusters. If these clusters reach the bloodstream intact, they can give rise to clusters of circulating tumor cells (CTCs), as have often been seen experimentally. Here, we review the operating principles of the core regulatory network for EMT/MET that acts as a "three-way" switch giving rise to three distinct phenotypes - E, M and hybrid E/M - and present a theoretical framework that can elucidate the role of many other players in regulating epithelial plasticity. Furthermore, we highlight recent studies on partial EMT and its association with drug resistance and tumor-initiating potential; and discuss how cell-cell communication between cells in a partial EMT phenotype can enable the formation of clusters of CTCs. These clusters can be more apoptosis-resistant and have more tumor-initiating potential than singly moving CTCs with a wholly mesenchymal (complete EMT) phenotype. Also, more such clusters can be formed under inflammatory conditions that are often generated by various therapies. Finally, we discuss the multiple advantages that the partial EMT or hybrid E/M phenotype have as compared to a complete EMT phenotype and argue that these collectively migrating cells are the primary "bad actors" of metastasis.
  • Thumbnail Image
    Item
    Interrogating the topological robustness of gene regulatory circuits by randomization
    (Public Library of Science, 2017) Huang, Bin; Lu, Mingyang; Jia, Dongya; Ben-Jacob, Eshel; Levine, Herbert; Onuchic, José Nelson; Bioengineering; Biosciences; Chemistry; Physics and Astronomy; Center for Theoretical Biological Physics
    One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.
  • Thumbnail Image
    Item
    Modeling delayed processes in biological systems
    (American Physical Society, 2016) Feng, Jingchen; Sevier, Stuart A.; Huang, Bin; Jia, Dongya; Levine, Herbert; Bioengineering; Center for Theoretical Biological Physics
    Delayed processes are ubiquitous in biological systems and are often characterized by delay differential equations (DDEs) and their extension to include stochastic effects. DDEs do not explicitly incorporate intermediate states associated with a delayed process but instead use an estimated average delay time. In an effort to examine the validity of this approach, we study systems with significant delays by explicitly incorporating intermediate steps. We show that such explicit models often yield significantly different equilibrium distributions and transition times as compared to DDEs with deterministic delay values. Additionally, different explicit models with qualitatively different dynamics can give rise to the same DDEs revealing important ambiguities. We also show that DDE-based predictions of oscillatory behavior may fail for the corresponding explicit model.
  • Thumbnail Image
    Item
    Modeling the Transitions between Collective and Solitary Migration Phenotypes in Cancer Metastasis
    (Macmillan Publishers Limited, 2015) Huang, Bin; Jolly, Mohit Kumar; Lu, Mingyang; Tsarfaty, Ilan; Ben-Jacob, Eshel; Onuchic, José Nelson; Bioengineering; Biosciences; Chemistry; Physics and Astronomy; Center for Theoretical Biological Physics
    Cellular plasticity during cancer metastasis is a major clinical challenge. Two key cellular plasticity mechanisms —Epithelial-to-Mesenchymal Transition (EMT) and Mesenchymal-to-Amoeboid Transition (MAT) – have been carefully investigated individually, yet a comprehensive understanding of their interconnections remains elusive. Previously, we have modeled the dynamics of the core regulatory circuits for both EMT (miR-200/ZEB/miR-34/SNAIL) and MAT (Rac1/RhoA). We now extend our previous work to study the coupling between these two core circuits by considering the two microRNAs (miR-200 and miR-34) as external signals to the core MAT circuit. We show that this coupled circuit enables four different stable steady states (phenotypes) that correspond to hybrid epithelial/mesenchymal (E/M), mesenchymal (M), amoeboid (A) and hybrid amoeboid/mesenchymal (A/M) phenotypes. Our model recapitulates the metastasis-suppressing role of the microRNAs even in the presence of EMT-inducing signals like Hepatocyte Growth Factor (HGF). It also enables mapping the microRNA levels to the transitions among various cell migration phenotypes. Finally, it offers a mechanistic understanding for the observed phenotypic transitions among different cell migration phenotypes, specifically the Collective-to-Amoeboid Transition (CAT).
  • Thumbnail Image
    Item
    RACIPE: a computational tool for modeling gene regulatory circuits using randomization
    (Springer Nature, 2018) Huang, Bin; Jia, Dongya; Feng, Jingchen; Levine, Herbert; Onuchic, José Nelson; Lu, Mingyang; Bioengineering; Biosciences; Chemistry; Center for Theoretical Biological Physics
    BACKGROUND: One of the major challenges in traditional mathematical modeling of gene regulatory circuits is the insufficient knowledge of kinetic parameters. These parameters are often inferred from existing experimental data and/or educated guesses, which can be time-consuming and error-prone, especially for large networks. RESULTS: We present a user-friendly computational tool for the community to use our newly developed method named random circuit perturbation (RACIPE), to explore the robust dynamical features of gene regulatory circuits without the requirement of detailed kinetic parameters. Taking the network topology as the only input, RACIPE generates an ensemble of circuit models with distinct randomized parameters and uniquely identifies robust dynamical properties by statistical analysis. Here, we discuss the implementation of the software and the statistical analysis methods of RACIPE-generated data to identify robust gene expression patterns and the functions of genes and regulatory links. Finally, we apply the tool on coupled toggle-switch circuits and a published circuit of B-lymphopoiesis. CONCLUSIONS: We expect our new computational tool to contribute to a more comprehensive and unbiased understanding of mechanisms underlying gene regulatory networks. RACIPE is a free open source software distributed under (Apache 2.0) license and can be downloaded from GitHub ( https://github.com/simonhb1990/RACIPE-1.0 ).
  • Thumbnail Image
    Item
    The three-way switch operation of Rac1/RhoA GTPase-based circuit controlling amoeboid-hybrid-mesenchymal transition
    (Nature Publishing Group, 2014) Huang, Bin; Lu, Mingyang; Jolly, Mohit Kumar; Tsarfaty, Ilan; Onuchic, José Nelson; Ben-Jacob, Eshel; Bioengineering; Biosciences; Chemistry; Physics and Astronomy; Center for Theoretical Biological Physics
    Metastatic carcinoma cells exhibit at least two different phenotypes of motility and invasion - amoeboid and mesenchymal. This plasticity poses a major clinical challenge for treating metastasis, while its underlying mechanisms remain enigmatic. Transitions between these phenotypes are mediated by the Rac1/RhoA circuit that responds to external signals such as HGF/SF via c-MET pathway. Using detailed modeling of GTPase-based regulation to study the Rac1/RhoA circuit's dynamics, we found that it can operate as a three-way switch. We propose to associate the circuit's three possible states to the amoeboid, mesenchymal and amoeboid/mesenchymal hybrid phenotype. In particular, we investigated the range of existence of, and the transition between, the three states (phenotypes) in response to Grb2 and Gab1 - two downstream adaptors of c-MET. The results help to explain the regulation of metastatic cells by c-MET pathway and hence can contribute to the assessment of possible clinical interventions.