From molecular interaction to acute promyelocytic leukemia: Calculating leukemogenesis and remission from endogenous molecular-cellular network

Predictive modelling of acute Promyelocytic leukaemia resistance to retinoic acid therapy

Acute Promyelocytic Leukaemia (APL) arises from an aberrant chromosomal translocation involving the Retinoic Acid Receptor Alpha (RARA) gene, predominantly with the Promyelocytic Leukaemia (PML) or Promyelocytic Leukaemia Zinc Finger (PLZF) genes. The resulting oncoproteins block the haematopoietic differentiation program promoting aberrant proliferative promyelocytes. Retinoic Acid (RA) therapy is successful in most of the PML::RARA patients, while PLZF::RARA patients frequently become resistant and relapse. Recent studies pointed to various underlying molecular components, but their precise contributions remain to be deciphered. We developed a logical network model integrating signalling, transcriptional, and epigenetic regulatory mechanisms, which captures key features of the APL cell responses to RA depending on the genetic background. The explicit inclusion of the histone methyltransferase EZH2 allowed the assessment of its role in the resistance mechanism, distinguishing between its canonical and non-canonical activities. The model dynamics was thoroughly analysed using tools integrated in the public software suite maintained by the CoLoMoTo consortium (https://colomoto.github.io/). The model serves as a solid basis to assess the roles of novel regulatory mechanisms, as well as to explore novel therapeutical approaches in silico.

6-Gingerol modulates miRNAs and PODXL gene expression via methyltransferase enzymes in NB4 cells: an in silico and in vitro study

This investigation delves into the influence of predicted microRNAs on DNA methyltransferases (DNMTs) and the PODXL gene within the NB4 cell line, aiming to elucidate their roles in the pathogenesis of acute myeloid leukemia (AML). A comprehensive methodological framework was adopted to explore the therapeutic implications of 6-gingerol on DNMTs. This encompassed a suite of bioinformatics tools for protein structure prediction, docking, molecular dynamics, and ADMET profiling, alongside empirical assessments of miRNA and PODXL expression levels. Such a multifaceted strategy facilitated an in-depth understanding of 6-gingerol's potential efficacy in DNMT modulation. The findings indicate a nuanced interplay where 6-gingerol administration modulated miRNA expression levels, decreasing in DNMT1 and DNMT3A expression in NB4 cells. This alteration indirectly influenced PODXL expression, contributing to the manifestation of oncogenic phenotypes. The overexpression of DNMT1 and DNMT3A in NB4 cells may contribute to AML, which appears modulable via microRNAs such as miR-193a and miR-200c. Post-treatment with 6-gingerol, DNMT1 and DNMT3A expression alterations were observed, culminating in the upregulation of miR-193a and miR-200c. This cascade effect led to the dysregulation of tumor suppressor genes in cancer cells, including downregulation of PODXL, and the emergence of cancerous traits. These insights underscore the therapeutic promise of 6-gingerol in targeting DNMTs and microRNAs within the AML context.

Potential therapeutic targets of gastric cancer explored under endogenous network modeling of clinical data

Improvement in the survival rate of gastric cancer, a prevalent global malignancy and the leading cause of cancer-related mortality calls for more avenues in molecular therapy. This work aims to comprehend drug resistance and explore multiple-drug combinations for enhanced therapeutic treatment. An endogenous network modeling clinic data with core gastric cancer molecules, functional modules, and pathways is constructed, which is then transformed into dynamics equations for in-silicon studies. Principal component analysis, hierarchical clustering, and K-means clustering are utilized to map the attractor domains of the stochastic model to the normal and pathological phenotypes identified from the clinical data. The analyses demonstrate gastric cancer as a cluster of stable states emerging within the stochastic dynamics and elucidate the cause of resistance to anti-VEGF monotherapy in cancer treatment as the limitation of the single pathway in preventing cancer progression. The feasibility of multiple objectives of therapy targeting specified molecules and/or pathways is explored. This study verifies the rationality of the platform of endogenous network modeling, which contributes to the development of cross-functional multi-target combinations in clinical trials.

Circulating Biomarkers Associated with the Diagnosis and Prognosis of B-Cell Progenitor Acute Lymphoblastic Leukemia

Acute lymphoblastic leukemia (ALL) is a hematological disease characterized by the dysfunction of the hematopoietic system that leads to arrest at a specific stage of stem cells development, suppressing the average production of cellular hematologic components. BCP-ALL is a neoplasm of the B-cell lineage progenitor. BCP-ALL is caused and perpetuated by several mechanisms that provide the disease with its tumor potential and genetic and cytological characteristics. These pathological features are used for diagnosis and the prognostication of BCP-ALL. However, most of these paraclinical tools can only be obtained by bone marrow aspiration, which, as it is an invasive study, can delay the diagnosis and follow-up of the disease, in addition to the anesthetic risk it entails for pediatric patients. For this reason, it is crucial to find noninvasive and accessible ways to supply information concerning diagnosis, prognosis, and the monitoring of the disease, such as circulating biomarkers. In oncology, a biomarker is any measurable indicator that demonstrates the presence of malignancy, tumoral behavior, prognosis, or responses to treatments. This review summarizes circulating molecules associated with BCP-ALL with potential diagnostic value, classificatory capacity during monitoring specific clinic features of the disease, and/or capacity to identify each BCP-ALL stage regarding its evolution and outcome of the patients with BCP-ALL. In the same way, we provide and classify biomarkers that may be used in further studies focused on clinical approaches or therapeutic target identification for BCP-ALL.

Roadmap on plasticity and epigenetics in cancer

The role of plasticity and epigenetics in shaping cancer evolution and response to therapy has taken center stage with recent technological advances including single cell sequencing. This roadmap article is focused on state-of-the-art mathematical and experimental approaches to interrogate plasticity in cancer, and addresses the following themes and questions: is there a formal overarching framework that encompasses both non-genetic plasticity and mutation-driven somatic evolution? How do we measure and model the role of the microenvironment in influencing/controlling non-genetic plasticity? How can we experimentally study non-genetic plasticity? Which mathematical techniques are required or best suited? What are the clinical and practical applications and implications of these concepts?

The complex landscape of haematopoietic lineage commitments is encoded in the coarse-grained endogenous network

Haematopoietic lineage commitments are presented by a canonical roadmap in which haematopoietic stem cells or multipotent progenitors (MPPs) bifurcate into progenitors of more restricted lineages and ultimately mature to terminally differentiated cells. Although transcription factors playing significant roles in cell-fate commitments have been extensively studied, integrating such knowledge into the dynamic models to understand the underlying biological mechanism remains challenging. The hypothesis and modelling approach of the endogenous network has been developed previously and tested in various biological processes and is used in the present study of haematopoietic lineage commitments. The endogenous network is constructed based on the key transcription factors and their interactions that determine haematopoietic cell-fate decisions at each lineage branchpoint. We demonstrate that the process of haematopoietic lineage commitments can be reproduced from the landscape which orchestrates robust states of network dynamics and their transitions. Furthermore, some non-trivial characteristics are unveiled in the dynamical model. Our model also predicted previously under-represented regulatory interactions and heterogeneous MPP states by which distinct differentiation routes are intermediated. Moreover, network perturbations resulting in state transitions indicate the effects of ectopic gene expression on cellular reprogrammes. This study provides a predictive model to integrate experimental data and uncover the possible regulatory mechanism of haematopoietic lineage commitments.

Cancer Niches and Their Kikuchi Free Energy

Biological forms depend on a progressive specialization of pluripotent stem cells. The differentiation of these cells in their spatial and functional environment defines the organism itself; however, cellular mutations may disrupt the mutual balance between a cell and its niche, where cell proliferation and specialization are released from their autopoietic homeostasis. This induces the construction of cancer niches and maintains their survival. In this paper, we characterise cancer niche construction as a direct consequence of interactions between clusters of cancer and healthy cells. Explicitly, we evaluate these higher-order interactions between niches of cancer and healthy cells using Kikuchi approximations to the free energy. Kikuchi's free energy is measured in terms of changes to the sum of energies of baseline clusters of cells (or nodes) minus the energies of overcounted cluster intersections (and interactions of interactions, etc.). We posit that these changes in energy node clusters correspond to a long-term reduction in the complexity of the system conducive to cancer niche survival. We validate this formulation through numerical simulations of apoptosis, local cancer growth, and metastasis, and highlight its implications for a computational understanding of the etiopathology of cancer.

Adaptive Landscape Shaped by Core Endogenous Network Coordinates Complex Early Progenitor Fate Commitments in Embryonic Pancreas

The classical development hierarchy of pancreatic cell fate commitments describes that multipotent progenitors (MPs) first bifurcate into tip cells and trunk cells, and then these cells give rise to acinar cells and endocrine/ductal cells separately. However, lineage tracings reveal that pancreatic progenitors are highly heterogeneous in tip and trunk domains in embryonic pancreas. The progenitor fate commitments from multipotency to unipotency during early pancreas development is insufficiently characterized. In pursuing a mechanistic understanding of the complexity in progenitor fate commitments, we construct a core endogenous network for pancreatic lineage decisions based on genetic regulations and quantified its intrinsic dynamic properties using dynamic modeling. The dynamics reveal a developmental landscape with high complexity that has not been clarified. Not only well-characterized pancreatic cells are reproduced, but also previously unrecognized progenitors-tip progenitor (TiP), trunk progenitor (TrP), later endocrine progenitor (LEP), and acinar progenitors (AciP/AciP2) are predicted. Further analyses show that TrP and LEP mediate endocrine lineage maturation, while TiP, AciP, AciP2 and TrP mediate acinar and ductal lineage maturation. The predicted cell fate commitments are validated by analyzing single-cell RNA sequencing (scRNA-seq) data. Significantly, this is the first time that a redefined hierarchy with detailed early pancreatic progenitor fate commitment is obtained.

ATPR-induced differentiation and G0/G1 phase arrest in acute promyelocytic leukemia by repressing EBP50/NCF1 complex to promote the production of ROS

Our previous study has demonstrated that 4-amino-2-trifluoromethyl-phenyl Retinate (ATPR) can induce human leukemia NB4 cells differentiation and G0/G1 phase arrest, but the underlying mechanism is still unclear. In this study, we used proteomics to screen differentially expressed protein profiles in NB4 cells before and after ATPR treatment in vitro. We analyzed the peptides digested from total cellular proteins by reverse phase LC-MS/MS and then performed label-free quantitative analysis. We found 27 significantly up-regulated proteins in the ATPR group compared to the control group. NCF1 was the most significantly changed protein. Immunoprecipitation and double immunofluorescent staining showed that EBP50 bind to NCF1. We further explored the potential molecular mechanism of EBP50/NCF1 complex in ATPR-induced differentiation and G0/G1 phase arrest. The results showed that ATPR remarkably reduced the expression of EBP50 in vivo and in vitro. Interestingly, the reduction of EBP50 contributed to ROS release by modulating the subcellular localization of NCF1. The reduction of EBP50 also contributed to G0/G1 phase arrest by inhibiting CyclinD1, CyclinA2 and CDK4, as well as promoting the differentiation of NB4 cells by increasing the expression of CD11b. Furthermore, we found that the overexpression of EBP50 restrained the effects of ATPR on differentiation and G0/G1 phase arrest in NB4 cells. These results suggest that ATPR-induced differentiation and G0/G1 phase arrest in acute promyelocytic leukemia (APL) by repressing EBP50/NCF1 complex to promote the production of ROS, and the results from in vivo experiments were consistent with those from in vitro studies. Therefore, our finding results suggest that EBP50 may be a new target for ATPR in the treatment of APL.

Endogenous Molecular-Cellular Network Cancer Theory: A Systems Biology Approach

In light of ever apparent limitation of the current dominant cancer mutation theory, a quantitative hypothesis for cancer genesis and progression, endogenous molecular-cellular network hypothesis has been proposed from the systems biology perspective, now for more than 10 years. It was intended to include both the genetic and epigenetic causes to understand cancer. Its development enters the stage of meaningful interaction with experimental and clinical data and the limitation of the traditional cancer mutation theory becomes more evident. Under this endogenous network hypothesis, we established a core working network of hepatocellular carcinoma (HCC) according to the hypothesis and quantified the working network by a nonlinear dynamical system. We showed that the two stable states of the working network reproduce the main known features of normal liver and HCC at both the modular and molecular levels. Using endogenous network hypothesis and validated working network, we explored genetic mutation pattern in cancer and potential strategies to cure or relieve HCC from a totally new perspective. Patterns of genetic mutations have been traditionally analyzed by posteriori statistical association approaches in light of traditional cancer mutation theory. One may wonder the possibility of a priori determination of any mutation regularity. Here, we found that based on the endogenous network theory the features of genetic mutations in cancers may be predicted without any prior knowledge of mutation propensities. Normal hepatocyte and cancerous hepatocyte stable states, specified by distinct patterns of expressions or activities of proteins in the network, provide means to directly identify a set of most probable genetic mutations and their effects in HCC. As the key proteins and main interactions in the network are conserved through cell types in an organism, similar mutational features may also be found in other cancers. This analysis yielded straightforward and testable predictions on an accumulated and preferred mutation spectrum in normal tissue. The validation of predicted cancer state mutation patterns demonstrates the usefulness and potential of a causal dynamical framework to understand and predict genetic mutations in cancer. We also obtained the following implication related to HCC therapy, (1) specific positive feedback loops are responsible for the maintenance of normal liver and HCC; (2) inhibiting proliferation and inflammation-related positive feedback loops, and simultaneously inducing liver-specific positive feedback loop is predicated as the potential strategy to cure or relieve HCC; (3) the genesis and regression of HCC is asymmetric. In light of the characteristic property of the nonlinear dynamical system, we demonstrate that positive feedback loops must be existed as a simple and general molecular basis for the maintenance of phenotypes such as normal liver and HCC, and regulating the positive feedback loops directly or indirectly provides potential strategies to cure or relieve HCC.

Cell cycle time series gene expression data encoded as cyclic attractors in Hopfield systems

Modern time series gene expression and other omics data sets have enabled unprecedented resolution of the dynamics of cellular processes such as cell cycle and response to pharmaceutical compounds. In anticipation of the proliferation of time series data sets in the near future, we use the Hopfield model, a recurrent neural network based on spin glasses, to model the dynamics of cell cycle in HeLa (human cervical cancer) and S. cerevisiae cells. We study some of the rich dynamical properties of these cyclic Hopfield systems, including the ability of populations of simulated cells to recreate experimental expression data and the effects of noise on the dynamics. Next, we use a genetic algorithm to identify sets of genes which, when selectively inhibited by local external fields representing gene silencing compounds such as kinase inhibitors, disrupt the encoded cell cycle. We find, for example, that inhibiting the set of four kinases AURKB, NEK1, TTK, and WEE1 causes simulated HeLa cells to accumulate in the M phase. Finally, we suggest possible improvements and extensions to our model.

Cell cycle—the process in which a parent cell replicates its DNA and divides into two daughter cells—is an upregulated process in many forms of cancer. Identifying gene inhibition targets to regulate cell cycle is important to the development of effective therapies. Although modern high throughput techniques offer unprecedented resolution of the molecular details of biological processes like cell cycle, analyzing the vast quantities of the resulting experimental data and extracting actionable information remains a formidable task. Here, we create a dynamical model of the process of cell cycle using the Hopfield model (a type of recurrent neural network) and gene expression data from human cervical cancer cells and yeast cells. We find that the model recreates the oscillations observed in experimental data. Tuning the level of noise (representing the inherent randomness in gene expression and regulation) to the “edge of chaos” is crucial for the proper behavior of the system. We then use this model to identify potential gene targets for disrupting the process of cell cycle. This method could be applied to other time series data sets and used to predict the effects of untested targeted perturbations.

Beyond cancer genes: colorectal cancer as robust intrinsic states formed by molecular interactions

Colorectal cancer (CRC) has complex pathological features that defy the linear-additive reasoning prevailing in current biomedicine studies. In pursuing a mechanistic understanding behind such complexity, we constructed a core molecular–cellular interaction network underlying CRC and investigated its nonlinear dynamical properties. The hypothesis and modelling method has been developed previously and tested in various cancer studies. The network dynamics reveal a landscape of several attractive basins corresponding to both normal intestinal phenotype and robust tumour subtypes, identified by their different molecular signatures. Comparison between the modelling results and gene expression profiles from patients collected at the second affiliated hospital of Zhejiang University is presented as validation. The numerical ‘driving’ experiment suggests that CRC pathogenesis may depend on pathways involved in gastrointestinal track development and molecules associated with mesenchymal lineage differentiation, such as Stat5, BMP, retinoic acid signalling pathways, Runx and Hox transcription families. We show that the multi-faceted response to immune stimulation and therapies, as well as different carcinogenesis and metastasis routes, can be straightforwardly understood and analysed under such a framework.

Evolutionary Origins of Cancer Driver Genes and Implications for Cancer Prognosis

The cancer atavistic theory suggests that carcinogenesis is a reverse evolution process. It is thus of great interest to explore the evolutionary origins of cancer driver genes and the relevant mechanisms underlying the carcinogenesis. Moreover, the evolutionary features of cancer driver genes could be helpful in selecting cancer biomarkers from high-throughput data. In this study, through analyzing the cancer endogenous molecular networks, we revealed that the subnetwork originating from eukaryota could control the unlimited proliferation of cancer cells, and the subnetwork originating from eumetazoa could recapitulate the other hallmarks of cancer. In addition, investigations based on multiple datasets revealed that cancer driver genes were enriched in genes originating from eukaryota, opisthokonta, and eumetazoa. These results have important implications for enhancing the robustness of cancer prognosis models through selecting the gene signatures by the gene age information.

Decoding early myelopoiesis from dynamics of core endogenous network

A decade ago mainstream molecular biologists regarded it impossible or biologically ill-motivated to understand the dynamics of complex biological phenomena, such as cancer genesis and progression, from a network perspective. Indeed, there are numerical difficulties even for those who were determined to explore along this direction. Undeterred, seven years ago a group of Chinese scientists started a program aiming to obtain quantitative connections between tumors and network dynamics. Many interesting results have been obtained. In this paper we wish to test such idea from a different angle: the connection between a normal biological process and the network dynamics. We have taken early myelopoiesis as our biological model. A standard roadmap for the cell-fate diversification during hematopoiesis has already been well established experimentally, yet little was known for its underpinning dynamical mechanisms. Compounding this difficulty there were additional experimental challenges, such as the seemingly conflicting hematopoietic roadmaps and the cell-fate inter-conversion events. With early myeloid cell-fate determination in mind, we constructed a core molecular endogenous network from well-documented gene regulation and signal transduction knowledge. Turning the network into a set of dynamical equations, we found computationally several structurally robust states. Those states nicely correspond to known cell phenotypes. We also found the states connecting those stable states. They reveal the developmental routes-how one stable state would most likely turn into another stable state. Such interconnected network among stable states enabled a natural organization of cell-fates into a multi-stable state landscape. Accordingly, both the myeloid cell phenotypes and the standard roadmap were explained mechanistically in a straightforward manner. Furthermore, recent challenging observations were also explained naturally. Moreover, the landscape visually enables a prediction of a pool of additional cell states and developmental routes, including the non-sequential and cross-branch transitions, which are testable by future experiments. In summary, the endogenous network dynamics provide an integrated quantitative framework to understand the heterogeneity and lineage commitment in myeloid progenitors.

Cancer as robust intrinsic state shaped by evolution: a key issues review

Cancer is a complex disease: its pathology cannot be properly understood in terms of independent players-genes, proteins, molecular pathways, or their simple combinations. This is similar to many-body physics of a condensed phase that many important properties are not determined by a single atom or molecule. The rapidly accumulating large 'omics' data also require a new mechanistic and global underpinning to organize for rationalizing cancer complexity. A unifying and quantitative theory was proposed by some of the present authors that cancer is a robust state formed by the endogenous molecular-cellular network, which is evolutionarily built for the developmental processes and physiological functions. Cancer state is not optimized for the whole organism. The discovery of crucial players in cancer, together with their developmental and physiological roles, in turn, suggests the existence of a hierarchical structure within molecular biology systems. Such a structure enables a decision network to be constructed from experimental knowledge. By examining the nonlinear stochastic dynamics of the network, robust states corresponding to normal physiological and abnormal pathological phenotypes, including cancer, emerge naturally. The nonlinear dynamical model of the network leads to a more encompassing understanding than the prevailing linear-additive thinking in cancer research. So far, this theory has been applied to prostate, hepatocellular, gastric cancers and acute promyelocytic leukemia with initial success. It may offer an example of carrying physics inquiring spirit beyond its traditional domain: while quantitative approaches can address individual cases, however there must be general rules/laws to be discovered in biology and medicine.

Autophagy contributes to 4-Amino-2-Trifluoromethyl-Phenyl Retinate-induced differentiation in human acute promyelocytic leukemia NB4 cells

As a classic differentiation agent, all-trans retinoic acid (ATRA) has been widely used in treatment of acute promyelocytic leukemia (APL). However, clinical application of ATRA has limitations. Our previous studies suggested that 4-Amino-2-Trifluoromethyl-Phenyl Retinate (ATPR), a novel all-trans retinoic acid (ATRA) derivative designed and synthesized by our team, could induce differentiation of APL cells in vivo and in vitro. To explore the underlying mechanism of ATPR, the effect of ATPR on autophagy of APL cells was observed in the present study. The results showed that the differentiation effect of ATPR on APL cells was accompanied with autophagy induction and PML-RARα degradation via activating Notch1 signaling pathway. Moreover, inhibition of autophagy using 3-methyladenine (3-MA) or small interfering RNA (siRNA) that targets essential autophagy gene ATG5 abrogated the ATPR-induced cell differentiation. Furthermore, when pretreated with DAPT, a γ-secretase inhibitor, the Notch1 signaling pathway was blocked in APL cells, followed by the reduction of ATPR-induced autophagy and differentiation. Taken together, these results suggested that autophagy play an important role in ATPR-induced cell differentiation, which may provide a novel approach to cure APL patients.

Boolean Dynamic Modeling Approaches to Study Plant Gene Regulatory Networks: Integration, Validation, and Prediction

Mathematical models based on dynamical systems theory are well-suited tools for the integration of available molecular experimental data into coherent frameworks in order to propose hypotheses about the cooperative regulatory mechanisms driving developmental processes. Computational analysis of the proposed models using well-established methods enables testing the hypotheses by contrasting predictions with observations. Within such framework, Boolean gene regulatory network dynamical models have been extensively used in modeling plant development. Boolean models are simple and intuitively appealing, ideal tools for collaborative efforts between theorists and experimentalists. In this chapter we present protocols used in our group for the study of diverse plant developmental processes. We focus on conceptual clarity and practical implementation, providing directions to the corresponding technical literature.