1
|
Kim Y, Han Y, Hopper C, Lee J, Joo JI, Gong JR, Lee CK, Jang SH, Kang J, Kim T, Cho KH. A gray box framework that optimizes a white box logical model using a black box optimizer for simulating cellular responses to perturbations. CELL REPORTS METHODS 2024; 4:100773. [PMID: 38744288 PMCID: PMC11133856 DOI: 10.1016/j.crmeth.2024.100773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/13/2023] [Revised: 03/19/2024] [Accepted: 04/19/2024] [Indexed: 05/16/2024]
Abstract
Predicting cellular responses to perturbations requires interpretable insights into molecular regulatory dynamics to perform reliable cell fate control, despite the confounding non-linearity of the underlying interactions. There is a growing interest in developing machine learning-based perturbation response prediction models to handle the non-linearity of perturbation data, but their interpretation in terms of molecular regulatory dynamics remains a challenge. Alternatively, for meaningful biological interpretation, logical network models such as Boolean networks are widely used in systems biology to represent intracellular molecular regulation. However, determining the appropriate regulatory logic of large-scale networks remains an obstacle due to the high-dimensional and discontinuous search space. To tackle these challenges, we present a scalable derivative-free optimizer trained by meta-reinforcement learning for Boolean network models. The logical network model optimized by the trained optimizer successfully predicts anti-cancer drug responses of cancer cell lines, while simultaneously providing insight into their underlying molecular regulatory mechanisms.
Collapse
Affiliation(s)
- Yunseong Kim
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Younghyun Han
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Corbin Hopper
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Jonghoon Lee
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Jae Il Joo
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Jeong-Ryeol Gong
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Chun-Kyung Lee
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Seong-Hoon Jang
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Junsoo Kang
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Taeyoung Kim
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea
| | - Kwang-Hyun Cho
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
2
|
Hu Q, Luo M, Wang R. Identifying critical regulatory interactions in cell fate decision and transition by systematic perturbation analysis. J Theor Biol 2024; 577:111673. [PMID: 37984586 DOI: 10.1016/j.jtbi.2023.111673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2023] [Revised: 11/11/2023] [Accepted: 11/16/2023] [Indexed: 11/22/2023]
Abstract
One of the most significant challenges in biology is to elucidate the roles of various regulatory interactions in cell fate decision and transition. However, it remains to be fully clarified how they cooperate and determine fate transition. Here, a general framework based on statistical analysis and bifurcation theory is proposed to identify crucial regulatory interactions and how they play decisive roles in fate transition. More exactly, specific feedback loops determine occurrence of bifurcations by which cell fate transition can be realized. While regulatory interactions in the feedback loops determine the direction of transition. In addition, two-parameter bifurcation analysis further provides detailed understanding of how the fate transition based on statistical analysis occurs. Statistical analysis can also be used to reveal synergistic combinatorial perturbations by which fate transition can be more efficiently realized. The integrative analysis approach can be used to identify critical regulatory interactions in cell fate transition and reveal how specific cell fate transition occurs. To verify feasibility of the approach, the epithelial to mesenchymal transition (EMT) network is chosen as an illustrative example. In agreement with experimental observations, the approach reveals some critical regulatory interactions and underlying mechanisms in cell fate determination and transitions between three states. The approach can also be applied to analyze other regulatory networks related to cell fate decision and transition.
Collapse
Affiliation(s)
- Qing Hu
- Department of Mathematics, Shanghai University, Shanghai, 200444, China
| | - Min Luo
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | - Ruiqi Wang
- Department of Mathematics, Shanghai University, Shanghai, 200444, China; Newtouch Center for Mathematics of Shanghai University, Shanghai, 200444, China.
| |
Collapse
|
3
|
Zhuravleva SI, Zadorozhny AD, Shilov BV, Lagunin AA. Prediction of Amino Acid Substitutions in ABL1 Protein Leading to Tumor Drug Resistance Based on "Structure-Property" Relationship Classification Models. Life (Basel) 2023; 13:1807. [PMID: 37763211 PMCID: PMC10532460 DOI: 10.3390/life13091807] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2023] [Revised: 08/15/2023] [Accepted: 08/21/2023] [Indexed: 09/29/2023] Open
Abstract
Drug resistance to anticancer drugs is a serious complication in patients with cancer. Typically, drug resistance occurs due to amino acid substitutions (AAS) in drug target proteins. The study aimed at developing and validating a new approach to the creation of structure-property relationships (SPR) classification models to predict AASs leading to drug resistance to inhibitors of tyrosine-protein kinase ABL1. The approach was based on the representation of AASs as peptides described in terms of structural formulas. The data on drug-resistant and non-resistant variants of AAS for two isoforms of ABL1 were extracted from the COSMIC database. The given training sets (approximately 700 missense variants) were used for the creation of SPR models in MultiPASS software based on substructural atom-centric multiple neighborhoods of atom (MNA) descriptors for the description of the structural formula of protein fragments and a Bayesian-like algorithm for revealing structure-property relationships. It was found that MNA descriptors of the 6th level and peptides from 11 amino acid residues were the best combination for ABL1 isoform 1 with the prediction accuracy (AUC) of resistance to imatinib (0.897) and dasatinib (0.996). For ABL1 isoform 2 (resistance to imatinib), the best combination was MNA descriptors of the 6th level, peptides form 15 amino acids (AUC value was 0.909). The prediction of possible drug-resistant AASs was made for dbSNP and gnomAD data. The six selected most probable imatinib-resistant AASs were additionally validated by molecular modeling and docking, which confirmed the possibility of resistance for the E334V and T392I variants.
Collapse
Affiliation(s)
- Svetlana I. Zhuravleva
- Department of Bioinformatics, Pirogov Russian National Research Medical University, 117997 Moscow, Russia; (S.I.Z.); (A.D.Z.); (B.V.S.)
| | - Anton D. Zadorozhny
- Department of Bioinformatics, Pirogov Russian National Research Medical University, 117997 Moscow, Russia; (S.I.Z.); (A.D.Z.); (B.V.S.)
| | - Boris V. Shilov
- Department of Bioinformatics, Pirogov Russian National Research Medical University, 117997 Moscow, Russia; (S.I.Z.); (A.D.Z.); (B.V.S.)
| | - Alexey A. Lagunin
- Department of Bioinformatics, Pirogov Russian National Research Medical University, 117997 Moscow, Russia; (S.I.Z.); (A.D.Z.); (B.V.S.)
- Department of Bioinformatics, Institute of Biomedical Chemistry, 119121 Moscow, Russia
| |
Collapse
|
4
|
Joo JI, Park H, Cho K. Normalizing Input-Output Relationships of Cancer Networks for Reversion Therapy. ADVANCED SCIENCE (WEINHEIM, BADEN-WURTTEMBERG, GERMANY) 2023; 10:e2207322. [PMID: 37269056 PMCID: PMC10460890 DOI: 10.1002/advs.202207322] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Revised: 03/17/2023] [Indexed: 06/04/2023]
Abstract
Accumulated genetic alterations in cancer cells distort cellular stimulus-response (or input-output) relationships, resulting in uncontrolled proliferation. However, the complex molecular interaction network within a cell implicates a possibility of restoring such distorted input-output relationships by rewiring the signal flow through controlling hidden molecular switches. Here, a system framework of analyzing cellular input-output relationships in consideration of various genetic alterations and identifying possible molecular switches that can normalize the distorted relationships based on Boolean network modeling and dynamics analysis is presented. Such reversion is demonstrated by the analysis of a number of cancer molecular networks together with a focused case study on bladder cancer with in vitro experiments and patient survival data analysis. The origin of reversibility from an evolutionary point of view based on the redundancy and robustness intrinsically embedded in complex molecular regulatory networks is further discussed.
Collapse
Affiliation(s)
- Jae Il Joo
- Department of Bio and Brain EngineeringKorea Advanced Institute of Science and Technology (KAIST)Daejeon34141Republic of Korea
- Present address:
biorevert IncDaejeon34051Republic of Korea
| | - Hwa‐Jeong Park
- Department of Bio and Brain EngineeringKorea Advanced Institute of Science and Technology (KAIST)Daejeon34141Republic of Korea
- Present address:
Promega Corporationan affiliate of PromegaSouth Korea
| | - Kwang‐Hyun Cho
- Department of Bio and Brain EngineeringKorea Advanced Institute of Science and Technology (KAIST)Daejeon34141Republic of Korea
| |
Collapse
|
5
|
Kim J, Hopper C, Cho KH. Statistical control of structural networks with limited interventions to minimize cellular phenotypic diversity represented by point attractors. Sci Rep 2023; 13:6275. [PMID: 37072458 PMCID: PMC10113376 DOI: 10.1038/s41598-023-33346-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 04/12/2023] [Indexed: 05/03/2023] Open
Abstract
The underlying genetic networks of cells give rise to diverse behaviors known as phenotypes. Control of this cellular phenotypic diversity (CPD) may reveal key targets that govern differentiation during development or drug resistance in cancer. This work establishes an approach to control CPD that encompasses practical constraints, including model limitations, the number of simultaneous control targets, which targets are viable for control, and the granularity of control. Cellular networks are often limited to the structure of interactions, due to the practical difficulty of modeling interaction dynamics. However, these dynamics are essential to CPD. In response, our statistical control approach infers the CPD directly from the structure of a network, by considering an ensemble average function over all possible Boolean dynamics for each node in the network. These ensemble average functions are combined with an acyclic form of the network to infer the number of point attractors. Our approach is applied to several known biological models and shown to outperform existing approaches. Statistical control of CPD offers a new avenue to contend with systemic processes such as differentiation and cancer, despite practical limitations in the field.
Collapse
Affiliation(s)
- Jongwan Kim
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Corbin Hopper
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.
| |
Collapse
|
6
|
Shin D, Cho KH. Critical transition and reversion of tumorigenesis. Exp Mol Med 2023; 55:692-705. [PMID: 37009794 PMCID: PMC10167317 DOI: 10.1038/s12276-023-00969-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Revised: 01/10/2023] [Accepted: 01/11/2023] [Indexed: 04/04/2023] Open
Abstract
Cancer is caused by the accumulation of genetic alterations and therefore has been historically considered to be irreversible. Intriguingly, several studies have reported that cancer cells can be reversed to be normal cells under certain circumstances. Despite these experimental observations, conceptual and theoretical frameworks that explain these phenomena and enable their exploration in a systematic way are lacking. In this review, we provide an overview of cancer reversion studies and describe recent advancements in systems biological approaches based on attractor landscape analysis. We suggest that the critical transition in tumorigenesis is an important clue for achieving cancer reversion. During tumorigenesis, a critical transition may occur at a tipping point, where cells undergo abrupt changes and reach a new equilibrium state that is determined by complex intracellular regulatory events. We introduce a conceptual framework based on attractor landscapes through which we can investigate the critical transition in tumorigenesis and induce its reversion by combining intracellular molecular perturbation and extracellular signaling controls. Finally, we present a cancer reversion therapy approach that may be a paradigm-changing alternative to current cancer cell-killing therapies.
Collapse
Affiliation(s)
- Dongkwan Shin
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
- Reasearch Institute, National Cancer Center, Goyang, 10408, Republic of Korea
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.
| |
Collapse
|
7
|
Kim N, Hwang CY, Kim T, Kim H, Cho KH. A Cell-Fate Reprogramming Strategy Reverses Epithelial-to-Mesenchymal Transition of Lung Cancer Cells While Avoiding Hybrid States. Cancer Res 2023; 83:956-970. [PMID: 36710400 PMCID: PMC10015224 DOI: 10.1158/0008-5472.can-22-1559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 10/23/2022] [Accepted: 01/11/2023] [Indexed: 01/31/2023]
Abstract
The epithelial-to-mesenchymal transition (EMT) of primary cancer contributes to the acquisition of lethal properties, including metastasis and drug resistance. Blocking or reversing EMT could be an effective strategy to improve cancer treatment. However, it is still unclear how to achieve complete EMT reversal (rEMT), as cancer cells often transition to hybrid EMT states with high metastatic potential. To tackle this problem, we employed a systems biology approach and identified a core-regulatory circuit that plays the primary role in driving rEMT without hybrid properties. Perturbation of any single node was not sufficient to completely revert EMT. Inhibition of both SMAD4 and ERK signaling along with p53 activation could induce rEMT in cancer cells even with TGFβ stimulation, a primary inducer of EMT. Induction of rEMT in lung cancer cells with the triple combination approach restored chemosensitivity. This cell-fate reprogramming strategy based on attractor landscapes revealed potential therapeutic targets that can eradicate metastatic potential by subverting EMT while avoiding hybrid states. SIGNIFICANCE Network modeling unravels the highly complex and plastic process regulating epithelial and mesenchymal states in cancer cells and discovers therapeutic interventions for reversing epithelial-to-mesenchymal transition and enhancing chemosensitivity.
Collapse
Affiliation(s)
- Namhee Kim
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Chae Young Hwang
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
- Aventi Inc., Daejeon, Republic of Korea
| | - Taeyoung Kim
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Hyunjin Kim
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
| | - Kwang-Hyun Cho
- Laboratory for Systems Biology and Bio-inspired Engineering, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, Republic of Korea
- Corresponding Author: Kwang-Hyun Cho, Department of Bio and Brain Engineering, KAIST, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea. Phone: 82-42-350-4325; Fax: 82-42-350-4310; E-mail:
| |
Collapse
|
8
|
An S, Jang SY, Park SM, Lee CK, Kim HM, Cho KH. Global stabilizing control of large-scale biomolecular regulatory networks. Bioinformatics 2023; 39:6998201. [PMID: 36688702 PMCID: PMC9891247 DOI: 10.1093/bioinformatics/btad045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2022] [Revised: 01/13/2023] [Accepted: 01/20/2023] [Indexed: 01/24/2023] Open
Abstract
MOTIVATION Cellular behavior is determined by complex non-linear interactions between numerous intracellular molecules that are often represented by Boolean network models. To achieve a desired cellular behavior with minimal intervention, we need to identify optimal control targets that can drive heterogeneous cellular states to the desired phenotypic cellular state with minimal node intervention. Previous attempts to realize such global stabilization were based solely on either network structure information or simple linear dynamics. Other attempts based on non-linear dynamics are not scalable. RESULTS Here, we investigate the underlying relationship between structurally identified control targets and optimal global stabilizing control targets based on non-linear dynamics. We discovered that optimal global stabilizing control targets can be identified by analyzing the dynamics between structurally identified control targets. Utilizing these findings, we developed a scalable global stabilizing control framework using both structural and dynamic information. Our framework narrows down the search space based on strongly connected components and feedback vertex sets then identifies global stabilizing control targets based on the canalization of Boolean network dynamics. We find that the proposed global stabilizing control is superior with respect to the number of control target nodes, scalability, and computational complexity. AVAILABILITY AND IMPLEMENTATION We provide a GitHub repository that contains the DCGS framework written in Python as well as biological random Boolean network datasets (https://github.com/sugyun/DCGS). SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
| | | | | | - Chun-Kyung Lee
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hoon-Min Kim
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | | |
Collapse
|
9
|
Evaluating a therapeutic window for precision medicine by integrating genomic profiles and p53 network dynamics. Commun Biol 2022; 5:924. [PMID: 36071176 PMCID: PMC9452682 DOI: 10.1038/s42003-022-03872-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Accepted: 08/23/2022] [Indexed: 11/08/2022] Open
Abstract
The response variation to anti-cancer drugs originates from complex intracellular network dynamics of cancer. Such dynamic networks present challenges to determining optimal drug targets and stratifying cancer patients for precision medicine, although several cancer genome studies provided insights into the molecular characteristics of cancer. Here, we introduce a network dynamics-based approach based on attractor landscape analysis to evaluate the therapeutic window of a drug from cancer signaling networks combined with genomic profiles. This approach allows for effective screening of drug targets to explore potential target combinations for enhancing the therapeutic window of drug responses. We also effectively stratify patients into desired/undesired response groups using critical genomic determinants, which are network-specific origins of variability to drug response, and their dominance relationship. Our methods provide a viable and quantitative framework to connect genotype information to the phenotypes of drug response with regard to network dynamics determining the therapeutic window.
Collapse
|
10
|
Therapeutics Targeting p53-MDM2 Interaction to Induce Cancer Cell Death. Int J Mol Sci 2022; 23:ijms23095005. [PMID: 35563397 PMCID: PMC9103871 DOI: 10.3390/ijms23095005] [Citation(s) in RCA: 37] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Revised: 04/26/2022] [Accepted: 04/27/2022] [Indexed: 02/06/2023] Open
Abstract
Named as the guardian of the genome, p53 is a tumor suppressor that regulates cell function, often through many different mechanisms such as DNA repair, apoptosis, cell cycle arrest, senescence, metabolism, and autophagy. One of the genes that p53 activates is MDM2, which forms a negative feedback loop since MDM2 induces the degradation of p53. When p53 activity is inhibited, damaged cells do not undergo cell cycle arrest or apoptosis. As 50% of human cancers inactivate p53 by mutation, current research focuses on reactivating p53 by developing drugs that target the p53-MDM2 interaction, which includes the binding of MDM2 and phosphorylation of p53. The objective of this article is to provide a short list and description of p53-MDM2 antagonists that may be excellent candidates for inducing cancer cell death. Relevant articles were searched for and identified using online databases such as PubMed and ScienceDirect. Increasing p53 levels, by targeting the p53-MDM2 interaction, can help p53 play its role as a tumor suppressor and induce cancer cell death. Researchers have identified different compounds that can act as inhibitors, either by directly binding to MDM2 or by modifying p53 with phosphorylation. The results associated with the drugs demonstrate the importance of targeting such interactions to inhibit cancer cell growth, which indicates that the use of the compounds may improve cancer therapeutics.
Collapse
|
11
|
Xie J, Zhang L, Liu B, Liang X, Shi J. Single-cell analysis of p53 transitional dynamics unravels stimulus- and cell type-dependent signaling output motifs. BMC Biol 2022; 20:85. [PMID: 35410287 PMCID: PMC9004066 DOI: 10.1186/s12915-022-01290-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2021] [Accepted: 03/31/2022] [Indexed: 11/19/2022] Open
Abstract
Background To understand functional changes of complex biological networks, mathematical modeling of network topologies provides a quantitative measure of the way biological systems adapt to external stimuli. However, systemic network topology-based analysis often generates conflicting evidence depending on specific experimental conditions, leading to a limited mechanistic understanding of signaling networks and their differential dynamic outputs, an example of which is the regulation of p53 pathway responses to different stress stimuli and in variable mammalian cell types. Here, we employ a network motif approach to dissect key regulatory units of the p53 pathway and elucidate how network activities at the motif level generate context-specific dynamic responses. Results By combining single-cell imaging and mathematical modeling of dose-dependent p53 dynamics induced by three chemotherapeutics of distinct mechanism-of-actions, including Etoposide, Nutlin-3a and 5-fluorouracil, and in five cancer cell types, we uncovered novel and highly variable p53 dynamic responses, in particular p53 transitional dynamics induced at intermediate drug concentrations, and identified the functional roles of distinct positive and negative feedback motifs of the p53 pathway in modulating the central p53-Mdm2 negative feedback to generate stimulus- and cell type-specific signaling responses. The mechanistic understanding of p53 network dynamics also revealed previously unknown mediators of anticancer drug actions and phenotypic variations in cancer cells that impact drug sensitivity. Conclusions Our results demonstrate that transitional dynamics of signaling proteins such as p53, activated at intermediate stimulus levels, vary the most between the dynamic outputs of different generic network motifs and can be employed as novel quantitative readouts to uncover and elucidate the key building blocks of large signaling networks. Our findings also provide new insight on drug mediators and phenotypic heterogeneity that underlie differential drug responses. Supplementary Information The online version contains supplementary material available at 10.1186/s12915-022-01290-7.
Collapse
Affiliation(s)
- Jun Xie
- Center for Quantitative Systems Biology, Department of Physics and Department of Biology, Hong Kong Baptist University, 224 Waterloo Road, Kowloon Tong, Kowloon, Hong Kong, China
| | - Lichun Zhang
- Center for Quantitative Systems Biology, Department of Physics and Department of Biology, Hong Kong Baptist University, 224 Waterloo Road, Kowloon Tong, Kowloon, Hong Kong, China
| | - Bodong Liu
- Center for Quantitative Systems Biology, Department of Physics and Department of Biology, Hong Kong Baptist University, 224 Waterloo Road, Kowloon Tong, Kowloon, Hong Kong, China
| | - Xiao Liang
- Center for Quantitative Systems Biology, Department of Physics and Department of Biology, Hong Kong Baptist University, 224 Waterloo Road, Kowloon Tong, Kowloon, Hong Kong, China
| | - Jue Shi
- Center for Quantitative Systems Biology, Department of Physics and Department of Biology, Hong Kong Baptist University, 224 Waterloo Road, Kowloon Tong, Kowloon, Hong Kong, China.
| |
Collapse
|
12
|
Manicka S, Marques-Pita M, Rocha LM. Effective connectivity determines the critical dynamics of biochemical networks. J R Soc Interface 2022; 19:20210659. [PMID: 35042384 PMCID: PMC8767216 DOI: 10.1098/rsif.2021.0659] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Accepted: 12/02/2021] [Indexed: 11/12/2022] Open
Abstract
Living systems comprise interacting biochemical components in very large networks. Given their high connectivity, biochemical dynamics are surprisingly not chaotic but quite robust to perturbations-a feature C.H. Waddington named canalization. Because organisms are also flexible enough to evolve, they arguably operate in a critical dynamical regime between order and chaos. The established theory of criticality is based on networks of interacting automata where Boolean truth values model presence/absence of biochemical molecules. The dynamical regime is predicted using network connectivity and node bias (to be on/off) as tuning parameters. Revising this to account for canalization leads to a significant improvement in dynamical regime prediction. The revision is based on effective connectivity, a measure of dynamical redundancy that buffers automata response to some inputs. In both random and experimentally validated systems biology networks, reducing effective connectivity makes living systems operate in stable or critical regimes even though the structure of their biochemical interaction networks predicts them to be chaotic. This suggests that dynamical redundancy may be naturally selected to maintain living systems near critical dynamics, providing both robustness and evolvability. By identifying how dynamics propagates preferably via effective pathways, our approach helps to identify precise ways to design and control network models of biochemical regulation and signalling.
Collapse
Affiliation(s)
- Santosh Manicka
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing and Engineering, Indiana University, Bloomington, IN, USA
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal
| | - Manuel Marques-Pita
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing and Engineering, Indiana University, Bloomington, IN, USA
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal
- Universidade Lusófona, CICANT and COPELABS, Campo Grande 388, 1700-097 Lisbon, Portugal
| | - Luis M. Rocha
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing and Engineering, Indiana University, Bloomington, IN, USA
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal
- Binghamton University, State University of New York, Binghamton, NY, USA
| |
Collapse
|
13
|
Miranda SP, Baião FA, Fleck JL, Piccolo SR. Predicting drug sensitivity of cancer cells based on DNA methylation levels. PLoS One 2021; 16:e0238757. [PMID: 34506489 PMCID: PMC8432830 DOI: 10.1371/journal.pone.0238757] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2020] [Accepted: 06/28/2021] [Indexed: 01/22/2023] Open
Abstract
Cancer cell lines, which are cell cultures derived from tumor samples, represent one of the least expensive and most studied preclinical models for drug development. Accurately predicting drug responses for a given cell line based on molecular features may help to optimize drug-development pipelines and explain mechanisms behind treatment responses. In this study, we focus on DNA methylation profiles as one type of molecular feature that is known to drive tumorigenesis and modulate treatment responses. Using genome-wide, DNA methylation profiles from 987 cell lines in the Genomics of Drug Sensitivity in Cancer database, we used machine-learning algorithms to evaluate the potential to predict cytotoxic responses for eight anti-cancer drugs. We compared the performance of five classification algorithms and four regression algorithms representing diverse methodologies, including tree-, probability-, kernel-, ensemble-, and distance-based approaches. We artificially subsampled the data to varying degrees, aiming to understand whether training based on relatively extreme outcomes would yield improved performance. When using classification or regression algorithms to predict discrete or continuous responses, respectively, we consistently observed excellent predictive performance when the training and test sets consisted of cell-line data. Classification algorithms performed best when we trained the models using cell lines with relatively extreme drug-response values, attaining area-under-the-receiver-operating-characteristic-curve values as high as 0.97. The regression algorithms performed best when we trained the models using the full range of drug-response values, although this depended on the performance metrics we used. Finally, we used patient data from The Cancer Genome Atlas to evaluate the feasibility of classifying clinical responses for human tumors based on models derived from cell lines. Generally, the algorithms were unable to identify patterns that predicted patient responses reliably; however, predictions by the Random Forests algorithm were significantly correlated with Temozolomide responses for low-grade gliomas.
Collapse
Affiliation(s)
- Sofia P. Miranda
- Department of Industrial Engineering, Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil
| | - Fernanda A. Baião
- Department of Industrial Engineering, Pontifical Catholic University of Rio de Janeiro, Rio de Janeiro, Rio de Janeiro, Brazil
| | - Julia L. Fleck
- Mines Saint-Etienne, Univ Clermont Auvergne, CNRS, UMR 6158 LIMOS, Centre CIS, Saint-Etienne, France
| | - Stephen R. Piccolo
- Department of Biology, Brigham Young University, Provo, Utah, United States of America
| |
Collapse
|
14
|
Lee HY, Jeon Y, Kim YK, Jang JY, Cho YS, Bhak J, Cho KH. Identifying molecular targets for reverse aging using integrated network analysis of transcriptomic and epigenomic changes during aging. Sci Rep 2021; 11:12317. [PMID: 34112891 PMCID: PMC8192508 DOI: 10.1038/s41598-021-91811-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Accepted: 05/25/2021] [Indexed: 01/08/2023] Open
Abstract
Aging is associated with widespread physiological changes, including skeletal muscle weakening, neuron system degeneration, hair loss, and skin wrinkling. Previous studies have identified numerous molecular biomarkers involved in these changes, but their regulatory mechanisms and functional repercussions remain elusive. In this study, we conducted next-generation sequencing of DNA methylation and RNA sequencing of blood samples from 51 healthy adults between 20 and 74 years of age and identified aging-related epigenetic and transcriptomic biomarkers. We also identified candidate molecular targets that can reversely regulate the transcriptomic biomarkers of aging by reconstructing a gene regulatory network model and performing signal flow analysis. For validation, we screened public experimental data including gene expression profiles in response to thousands of chemical perturbagens. Despite insufficient data on the binding targets of perturbagens and their modes of action, curcumin, which reversely regulated the biomarkers in the experimental dataset, was found to bind and inhibit JUN, which was identified as a candidate target via signal flow analysis. Collectively, our results demonstrate the utility of a network model for integrative analysis of omics data, which can help elucidate inter-omics regulatory mechanisms and develop therapeutic strategies against aging.
Collapse
Affiliation(s)
- Hwang-Yeol Lee
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.,Genome Research Institute, Clinomics Inc, Ulsan, 44919, Republic of Korea
| | - Yeonsu Jeon
- Department of Biomedical Engineering, College of Information and Biotechnology, Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea.,Korea Genomics Center (KOGIC), Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea
| | - Yeon Kyung Kim
- Department of Biomedical Engineering, College of Information and Biotechnology, Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea.,Korea Genomics Center (KOGIC), Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea
| | - Jae Young Jang
- Department of Biomedical Engineering, College of Information and Biotechnology, Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea.,Korea Genomics Center (KOGIC), Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea
| | - Yun Sung Cho
- Genome Research Institute, Clinomics Inc, Ulsan, 44919, Republic of Korea
| | - Jong Bhak
- Genome Research Institute, Clinomics Inc, Ulsan, 44919, Republic of Korea. .,Department of Biomedical Engineering, College of Information and Biotechnology, Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea. .,Korea Genomics Center (KOGIC), Ulsan National Institute of Science and Technology (UNIST), Ulsan, 44919, Republic of Korea. .,Personal Genomics Institute (PGI), Genome Research Foundation (GRF), Osong, 28160, Republic of Korea.
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea.
| |
Collapse
|
15
|
Sherekar S, Viswanathan GA. Boolean dynamic modeling of cancer signaling networks: Prognosis, progression, and therapeutics. COMPUTATIONAL AND SYSTEMS ONCOLOGY 2021. [DOI: 10.1002/cso2.1017] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Shubhank Sherekar
- Department of Chemical Engineering Indian Institute of Technology Bombay, Powai Mumbai India
| | - Ganesh A. Viswanathan
- Department of Chemical Engineering Indian Institute of Technology Bombay, Powai Mumbai India
| |
Collapse
|
16
|
Gates AJ, Brattig Correia R, Wang X, Rocha LM. The effective graph reveals redundancy, canalization, and control pathways in biochemical regulation and signaling. Proc Natl Acad Sci U S A 2021; 118:e2022598118. [PMID: 33737396 PMCID: PMC8000424 DOI: 10.1073/pnas.2022598118] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The ability to map causal interactions underlying genetic control and cellular signaling has led to increasingly accurate models of the complex biochemical networks that regulate cellular function. These network models provide deep insights into the organization, dynamics, and function of biochemical systems: for example, by revealing genetic control pathways involved in disease. However, the traditional representation of biochemical networks as binary interaction graphs fails to accurately represent an important dynamical feature of these multivariate systems: some pathways propagate control signals much more effectively than do others. Such heterogeneity of interactions reflects canalization-the system is robust to dynamical interventions in redundant pathways but responsive to interventions in effective pathways. Here, we introduce the effective graph, a weighted graph that captures the nonlinear logical redundancy present in biochemical network regulation, signaling, and control. Using 78 experimentally validated models derived from systems biology, we demonstrate that 1) redundant pathways are prevalent in biological models of biochemical regulation, 2) the effective graph provides a probabilistic but precise characterization of multivariate dynamics in a causal graph form, and 3) the effective graph provides an accurate explanation of how dynamical perturbation and control signals, such as those induced by cancer drug therapies, propagate in biochemical pathways. Overall, our results indicate that the effective graph provides an enriched description of the structure and dynamics of networked multivariate causal interactions. We demonstrate that it improves explainability, prediction, and control of complex dynamical systems in general and biochemical regulation in particular.
Collapse
Affiliation(s)
- Alexander J Gates
- Network Science Institute, Northeastern University, Boston, MA 02115;
| | - Rion Brattig Correia
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal
- Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, Ministry of Education of Brazil, 70040-020 Brasília, DF, Brazil
| | - Xuan Wang
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing & Engineering, Indiana University, Bloomington, IN 47408
| | - Luis M Rocha
- Instituto Gulbenkian de Ciência, 2780-156 Oeiras, Portugal;
- Center for Social and Biomedical Complexity, Luddy School of Informatics, Computing & Engineering, Indiana University, Bloomington, IN 47408
- Department of Systems Science and Industrial Engineering, Binghamton University, Binghamton, NY 13902
| |
Collapse
|
17
|
Mannheimer JD, Prasad A, Gustafson DL. Predicting chemosensitivity using drug perturbed gene dynamics. BMC Bioinformatics 2021; 22:15. [PMID: 33413081 PMCID: PMC7789515 DOI: 10.1186/s12859-020-03947-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Accepted: 12/22/2020] [Indexed: 11/20/2022] Open
Abstract
Background One of the current directions of precision medicine is the use of computational methods to aid in the diagnosis, prognosis, and treatment of disease based on data driven approaches. For instance, in oncology, there has been a particular focus on development of algorithms and biomarkers that can be used for pre-clinical and clinical applications. In particular large-scale omics-based models to predict drug sensitivity in in vitro cancer cell line panels have been used to explore the utility and aid in the development of these models as clinical tools. Additionally, a number of web-based interfaces have been constructed for researchers to explore the potential of drug perturbed gene expression as biomarkers including the NCI Transcriptional Pharmacodynamic Workbench. In this paper we explore the influence of drug perturbed gene dynamics of the NCI Transcriptional Pharmacodynamics Workbench in computational models to predict in vitro drug sensitivity for 15 drugs on the NCI60 cell line panel. Results This work presents three main findings. First, our models show that gene expression profiles that capture changes in gene expression after 24 h of exposure to a high concentration of drug generates the most accurate predictive models compared to the expression profiles under different dosing conditions. Second, signatures of 100 genes are developed for different gene expression profiles; furthermore, when the gene signatures are applied across gene expression profiles model performance is substantially decreased when gene signatures developed using changes in gene expression are applied to non-drugged gene expression. Lastly, we show that the gene interaction networks developed on these signatures show different network topologies and can be used to inform selection of cancer relevant genes. Conclusion Our models suggest that perturbed gene signatures are predictive of drug response, but cannot be applied to predict drug response using unperturbed gene expression. Furthermore, additional drug perturbed gene expression measurements in in vitro cell lines could generate more predictive models; but, more importantly be used in conjunction with computational methods to discover important drug disease relationships.
Collapse
Affiliation(s)
- Joshua D Mannheimer
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO, USA.,Flint Animal Cancer Center, Colorado State University, Fort Collins, CO, USA
| | - Ashok Prasad
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO, USA.,Department of Chemical and Biological Engineering, Colorado State University, Fort Collins, CO, USA
| | - Daniel L Gustafson
- School of Biomedical Engineering, Colorado State University, Fort Collins, CO, USA. .,Flint Animal Cancer Center, Colorado State University, Fort Collins, CO, USA. .,Department of Clinical Sciences, Colorado State University, Fort Collins, CO, USA. .,University of Colorado, Cancer Center Developmental Therapeutics Program, University of Colorado, Aurora, CO, USA.
| |
Collapse
|
18
|
Yin Z, Guo B, Ma S, Sun Y, Mi Z, Zheng Z. DReSS: a method to quantitatively describe the influence of structural perturbations on state spaces of genetic regulatory networks. Brief Bioinform 2020; 22:6032613. [PMID: 33313791 DOI: 10.1093/bib/bbaa315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Revised: 09/23/2020] [Accepted: 10/16/2020] [Indexed: 11/14/2022] Open
Abstract
Structures of genetic regulatory networks are not fixed. These structural perturbations can cause changes to the reachability of systems' state spaces. As system structures are related to genotypes and state spaces are related to phenotypes, it is important to study the relationship between structures and state spaces. However, there is still no method can quantitively describe the reachability differences of two state spaces caused by structural perturbations. Therefore, Difference in Reachability between State Spaces (DReSS) is proposed. DReSS index family can quantitively describe differences of reachability, attractor sets between two state spaces and can help find the key structure in a system, which may influence system's state space significantly. First, basic properties of DReSS including non-negativity, symmetry and subadditivity are proved. Then, typical examples are shown to explain the meaning of DReSS and the differences between DReSS and traditional graph distance. Finally, differences of DReSS distribution between real biological regulatory networks and random networks are compared. Results show most structural perturbations in biological networks tend to affect reachability inside and between attractor basins rather than to affect attractor set itself when compared with random networks, which illustrates that most genotype differences tend to influence the proportion of different phenotypes and only a few ones can create new phenotypes. DReSS can provide researchers with a new insight to study the relation between genotypes and phenotypes.
Collapse
Affiliation(s)
- Ziqiao Yin
- Shenyuan Honors College and School of Mathematical Sciences, Beihang University, and Key Laboratory of Mathematics, Informatics and Behavioral Semantics, Ministry of Education. He currently works as a visiting scholar at Yale University
| | - Binghui Guo
- Artificial Intelligence Institute, Beijing Advanced Innovation Center for Big Data and Brain Computing, LMIB, NLSDE, School of Mathematical Sciences, Beihang University, and Peng Cheng Laboratory
| | - Shuangge Ma
- Department of Biostatistics, Yale University
| | - Yifan Sun
- School of Statistics, Renmin University of China
| | - Zhilong Mi
- Key Laboratory of Mathematics, Informatics and Behavioral Semantics, Ministry of Education, and School of Mathematical Sciences from Beihang University
| | - Zhiming Zheng
- Artificial Intelligence Institute, Beijing Advanced Innovation Center for Big Data and Brain Computing, LMIB, NLSDE, School of Mathematical Sciences, Beihang University, and Peng Cheng Laboratory
| |
Collapse
|
19
|
Network integration and modelling of dynamic drug responses at multi-omics levels. Commun Biol 2020; 3:573. [PMID: 33060801 PMCID: PMC7567116 DOI: 10.1038/s42003-020-01302-8] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2020] [Accepted: 09/14/2020] [Indexed: 12/25/2022] Open
Abstract
Uncovering cellular responses from heterogeneous genomic data is crucial for molecular medicine in particular for drug safety. This can be realized by integrating the molecular activities in networks of interacting proteins. As proof-of-concept we challenge network modeling with time-resolved proteome, transcriptome and methylome measurements in iPSC-derived human 3D cardiac microtissues to elucidate adverse mechanisms of anthracycline cardiotoxicity measured with four different drugs (doxorubicin, epirubicin, idarubicin and daunorubicin). Dynamic molecular analysis at in vivo drug exposure levels reveal a network of 175 disease-associated proteins and identify common modules of anthracycline cardiotoxicity in vitro, related to mitochondrial and sarcomere function as well as remodeling of extracellular matrix. These in vitro-identified modules are transferable and are evaluated with biopsies of cardiomyopathy patients. This to our knowledge most comprehensive study on anthracycline cardiotoxicity demonstrates a reproducible workflow for molecular medicine and serves as a template for detecting adverse drug responses from complex omics data. Using a network propagation approach with integrated multi-omic data, Selevsek et al. develop a reproducible workflow for identifying drug toxicity effects in cellular systems. This is demonstrated with the analysis of anthracycline cardiotoxicity in cardiac microtissues under the effect of multiple drugs.
Collapse
|
20
|
Joo JI, Choi M, Jang SH, Choi S, Park SM, Shin D, Cho KH. Realizing Cancer Precision Medicine by Integrating Systems Biology and Nanomaterial Engineering. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2020; 32:e1906783. [PMID: 32253807 DOI: 10.1002/adma.201906783] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 12/19/2019] [Indexed: 06/11/2023]
Abstract
Many clinical trials for cancer precision medicine have yielded unsatisfactory results due to challenges such as drug resistance and low efficacy. Drug resistance is often caused by the complex compensatory regulation within the biomolecular network in a cancer cell. Recently, systems biological studies have modeled and simulated such complex networks to unravel the hidden mechanisms of drug resistance and identify promising new drug targets or combinatorial or sequential treatments for overcoming resistance to anticancer drugs. However, many of the identified targets or treatments present major difficulties for drug development and clinical application. Nanocarriers represent a path forward for developing therapies with these "undruggable" targets or those that require precise combinatorial or sequential application, for which conventional drug delivery mechanisms are unsuitable. Conversely, a challenge in nanomedicine has been low efficacy due to heterogeneity of cancers in patients. This problem can also be resolved through systems biological approaches by identifying personalized targets for individual patients or promoting the drug responses. Therefore, integration of systems biology and nanomaterial engineering will enable the clinical application of cancer precision medicine to overcome both drug resistance of conventional treatments and low efficacy of nanomedicine due to patient heterogeneity.
Collapse
Affiliation(s)
- Jae Il Joo
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Minsoo Choi
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Seong-Hoon Jang
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Sea Choi
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Sang-Min Park
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Dongkwan Shin
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| | - Kwang-Hyun Cho
- Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Daejeon, 34141, Republic of Korea
| |
Collapse
|
21
|
Abstract
Network theory provides one of the most potent analysis tools for the study of complex systems. In this paper, we illustrate the network-based perspective in drug research and how it is coherent with the new paradigm of drug discovery. We first present data sources from which networks are built, then show some examples of how the networks can be used to investigate drug-related systems. A section is devoted to network-based inference applications, i.e., prediction methods based on interactomes, that can be used to identify putative drug-target interactions without resorting to 3D modeling. Finally, we present some aspects of Boolean networks dynamics, anticipating that it might become a very potent modeling framework to develop in silico screening protocols able to simulate phenotypic screening experiments. We conclude that network applications integrated with machine learning and 3D modeling methods will become an indispensable tool for computational drug discovery in the next years.
Collapse
Affiliation(s)
- Maurizio Recanatini
- Department of Pharmacy and
Biotechnology, Alma Mater Studiorum—University of Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| | - Chiara Cabrelle
- Department of Pharmacy and
Biotechnology, Alma Mater Studiorum—University of Bologna, Via Belmeloro 6, I-40126 Bologna, Italy
| |
Collapse
|
22
|
Abstract
Making decisions on how best to treat cancer patients requires the integration of different data sets, including genomic profiles, tumour histopathology, radiological images, proteomic analysis and more. This wealth of biological information calls for novel strategies to integrate such information in a meaningful, predictive and experimentally verifiable way. In this Perspective we explain how executable computational models meet this need. Such models provide a means for comprehensive data integration, can be experimentally validated, are readily interpreted both biologically and clinically, and have the potential to predict effective therapies for different cancer types and subtypes. We explain what executable models are and how they can be used to represent the dynamic biological behaviours inherent in cancer, and demonstrate how such models, when coupled with automated reasoning, facilitate our understanding of the mechanisms by which oncogenic signalling pathways regulate tumours. We explore how executable models have impacted the field of cancer research and argue that extending them to represent a tumour in a specific patient (that is, an avatar) will pave the way for improved personalized treatments and precision medicine. Finally, we highlight some of the ongoing challenges in developing executable models and stress that effective cross-disciplinary efforts are key to forward progress in the field.
Collapse
Affiliation(s)
- Matthew A Clarke
- Department of Biochemistry, University of Cambridge, Cambridge, UK
| | - Jasmin Fisher
- Department of Biochemistry, University of Cambridge, Cambridge, UK.
- UCL Cancer Institute, University College London, London, UK.
| |
Collapse
|
23
|
Rajapakse VN, Herrada S, Lavi O. Phenotype stability under dynamic brain-tumor environment stimuli maps glioblastoma progression in patients. SCIENCE ADVANCES 2020; 6:eaaz4125. [PMID: 32832595 PMCID: PMC7439317 DOI: 10.1126/sciadv.aaz4125] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Accepted: 03/20/2020] [Indexed: 06/11/2023]
Abstract
Although tumor invasiveness is known to drive glioblastoma (GBM) recurrence, current approaches to treatment assume a fairly simple GBM phenotype transition map. We provide new analyses to estimate the likelihood of reaching or remaining in a phenotype under dynamic, physiologically likely perturbations of stimuli ("phenotype stability"). We show that higher stability values of the motile phenotype (Go) are associated with reduced patient survival. Moreover, induced motile states are capable of driving GBM recurrence. We found that the Dormancy and Go phenotypes are equally represented in advanced GBM samples, with natural transitioning between the two. Furthermore, Go and Grow phenotype transitions are mostly driven by tumor-brain stimuli. These are difficult to regulate directly, but could be modulated by reprogramming tumor-associated cell types. Our framework provides a foundation for designing targeted perturbations of the tumor-brain environment, by assessing their impact on GBM phenotypic plasticity, and is corroborated by analyses of patient data.
Collapse
Affiliation(s)
- Vinodh N. Rajapakse
- Integrative Cancer Dynamics Unit, Laboratory of Cell Biology, CCR, NCI, NIH, Bethesda, MD, USA
| | - Sylvia Herrada
- Integrative Cancer Dynamics Unit, Laboratory of Cell Biology, CCR, NCI, NIH, Bethesda, MD, USA
| | - Orit Lavi
- Integrative Cancer Dynamics Unit, Laboratory of Cell Biology, CCR, NCI, NIH, Bethesda, MD, USA
| |
Collapse
|
24
|
Feedback analysis identifies a combination target for overcoming adaptive resistance to targeted cancer therapy. Oncogene 2020; 39:3803-3820. [PMID: 32157217 DOI: 10.1038/s41388-020-1255-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2019] [Revised: 02/26/2020] [Accepted: 02/28/2020] [Indexed: 12/19/2022]
Abstract
Targeted drugs aim to treat cancer by directly inhibiting oncogene activity or oncogenic pathways, but drug resistance frequently emerges. Due to the intricate dynamics of cancer signaling networks, which contain complex feedback regulations, cancer cells can rewire these networks to adapt to and counter the cytotoxic effects of a drug, thereby limiting the efficacy of targeted therapies. To identify a combinatorial drug target that can overcome such a limitation, we developed a Boolean network simulation and analysis framework and applied this approach to a large-scale signaling network of colorectal cancer with integrated genomic information. We discovered Src as a critical combination drug target that can overcome the adaptive resistance to the targeted inhibition of mitogen-activated protein kinase pathway by blocking the essential feedback regulation responsible for resistance. The proposed framework is generic and can be widely used to identify drug targets that can overcome adaptive resistance to targeted therapies.
Collapse
|
25
|
Fujii K, Takeishi N, Hojo M, Inaba Y, Kawahara Y. Physically-interpretable classification of biological network dynamics for complex collective motions. Sci Rep 2020; 10:3005. [PMID: 32080208 PMCID: PMC7033192 DOI: 10.1038/s41598-020-58064-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2019] [Accepted: 12/13/2019] [Indexed: 11/09/2022] Open
Abstract
Understanding biological network dynamics is a fundamental issue in various scientific and engineering fields. Network theory is capable of revealing the relationship between elements and their propagation; however, for complex collective motions, the network properties often transiently and complexly change. A fundamental question addressed here pertains to the classification of collective motion network based on physically-interpretable dynamical properties. Here we apply a data-driven spectral analysis called graph dynamic mode decomposition, which obtains the dynamical properties for collective motion classification. Using a ballgame as an example, we classified the strategic collective motions in different global behaviours and discovered that, in addition to the physical properties, the contextual node information was critical for classification. Furthermore, we discovered the label-specific stronger spectra in the relationship among the nearest agents, providing physical and semantic interpretations. Our approach contributes to the understanding of principles of biological complex network dynamics from the perspective of nonlinear dynamical systems.
Collapse
Affiliation(s)
- Keisuke Fujii
- Graduate School of Informatics, Nagoya University, Nagoya, Japan. .,RIKEN Center for Advanced Intelligence Project, Tokyo, Japan.
| | - Naoya Takeishi
- RIKEN Center for Advanced Intelligence Project, Tokyo, Japan
| | - Motokazu Hojo
- RIKEN Center for Advanced Intelligence Project, Tokyo, Japan
| | - Yuki Inaba
- Japan Institute of Sports Sciences, Tokyo, Japan
| | - Yoshinobu Kawahara
- RIKEN Center for Advanced Intelligence Project, Tokyo, Japan.,Institute of Mathematics for Industry, Kyushu University, Fukuoka, Japan
| |
Collapse
|
26
|
Choi J, Park S, Ahn J. RefDNN: a reference drug based neural network for more accurate prediction of anticancer drug resistance. Sci Rep 2020; 10:1861. [PMID: 32024872 PMCID: PMC7002431 DOI: 10.1038/s41598-020-58821-x] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2019] [Accepted: 01/20/2020] [Indexed: 11/23/2022] Open
Abstract
Cancer is one of the most difficult diseases to treat owing to the drug resistance of tumour cells. Recent studies have revealed that drug responses are closely associated with genomic alterations in cancer cells. Numerous state-of-the-art machine learning models have been developed for prediction of drug responses using various genomic data and diverse drug molecular information, but those methods are ineffective to predict drug response to untrained drugs and gene expression patterns, which is known as the cold-start problem. In this study, we present a novel deep neural network model, termed RefDNN, for improved prediction of drug resistance and identification of biomarkers related to drug response. RefDNN exploits a collection of drugs, called reference drugs, to learn representations for a high-dimensional gene expression vector and a molecular structure vector of a drug and predicts drug response labels using the reference drug-based representations. These calculations come from the observation that similar chemicals have similar effects. The proposed model not only outperformed existing computational prediction models in most comparative experiments, but also showed more robust prediction for untrained drugs and cancer types than traditional machine learning models. RefDNN exploits the ElasticNet regularization to deal with high-dimensional gene expression data, which allows identification of gene markers associated with drug resistance. Lastly, we described an application of RefDNN in exploring a new candidate drug for liver cancer. As the proposed model can guarantee good prediction of drug responses to untrained drugs for given gene expression patterns, it may be of potential benefit in drug repositioning and personalized medicine.
Collapse
Affiliation(s)
- Jonghwan Choi
- Department of Computer Science, Yonsei University, Seoul, South Korea
| | - Sanghyun Park
- Department of Computer Science, Yonsei University, Seoul, South Korea.
| | - Jaegyoon Ahn
- Department of Computer Science & Engineering, Incheon National University, Incheon, South Korea.
| |
Collapse
|
27
|
Heudobler D, Lüke F, Vogelhuber M, Klobuch S, Pukrop T, Herr W, Gerner C, Pantziarka P, Ghibelli L, Reichle A. Anakoinosis: Correcting Aberrant Homeostasis of Cancer Tissue-Going Beyond Apoptosis Induction. Front Oncol 2019; 9:1408. [PMID: 31921665 PMCID: PMC6934003 DOI: 10.3389/fonc.2019.01408] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2019] [Accepted: 11/28/2019] [Indexed: 12/16/2022] Open
Abstract
The current approach to systemic therapy for metastatic cancer is aimed predominantly at inducing apoptosis of cancer cells by blocking tumor-promoting signaling pathways or by eradicating cell compartments within the tumor. In contrast, a systems view of therapy primarily considers the communication protocols that exist at multiple levels within the tumor complex, and the role of key regulators of such systems. Such regulators may have far-reaching influence on tumor response to therapy and therefore patient survival. This implies that neoplasia may be considered as a cell non-autonomous disease. The multi-scale activity ranges from intra-tumor cell compartments, to the tumor, to the tumor-harboring organ to the organism. In contrast to molecularly targeted therapies, a systems approach that identifies the complex communications networks driving tumor growth offers the prospect of disrupting or "normalizing" such aberrant communicative behaviors and therefore attenuating tumor growth. Communicative reprogramming, a treatment strategy referred to as anakoinosis, requires novel therapeutic instruments, so-called master modifiers to deliver concerted tumor growth-attenuating action. The diversity of biological outcomes following pro-anakoinotic tumor therapy, such as differentiation, trans-differentiation, control of tumor-associated inflammation, etc. demonstrates that long-term tumor control may occur in multiple forms, inducing even continuous complete remission. Accordingly, pro-anakoinotic therapies dramatically extend the repertoire for achieving tumor control and may activate apoptosis pathways for controlling resistant metastatic tumor disease and hematologic neoplasia.
Collapse
Affiliation(s)
- Daniel Heudobler
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Florian Lüke
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Martin Vogelhuber
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Sebastian Klobuch
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Tobias Pukrop
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Wolfgang Herr
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| | - Christopher Gerner
- Institut for Analytical Chemistry, Faculty Chemistry, University Vienna, Vienna, Austria
| | - Pan Pantziarka
- The George Pantziarka TP53 Trust, London, United Kingdom
- Anticancer Fund, Brussels, Belgium
| | - Lina Ghibelli
- Department Biology, Università di Roma Tor Vergata, Rome, Italy
| | - Albrecht Reichle
- Department of Internal Medicine III, Hematology and Oncology, University Hospital Regensburg, Regensburg, Germany
| |
Collapse
|
28
|
Hou X, Li M, Jia C, Zhang X, Wang Y. Attractor - a new turning point in drug discovery. DRUG DESIGN DEVELOPMENT AND THERAPY 2019; 13:2957-2968. [PMID: 31686779 PMCID: PMC6709805 DOI: 10.2147/dddt.s216397] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Accepted: 07/28/2019] [Indexed: 11/23/2022]
Abstract
Drug discovery for complex diseases can be viewed as a challenging problem in which the influence of compounds on dynamic features of disease system should be considered, especially the strategies escaping from the disease attractors. Moreover, escaping from the disease-related attractors has been proved to be a cue for the treatment of the complex diseases. The drug discovery methodology based on the attractor theory indicates new solutions for target identification, drug discovery and drug combination design. The methodology is based on the holism level of the organism and the features of system dynamics, so it has advantages for the classification of complex diseases and drug discovery. Currently, research results of this method have increased, which expand the insight scope for drug discovery. This article introduces the major drug discovery methods in the history of pharmacy development and their characteristics, so as to illustrate the reasons and inevitability of the appearance of attractor method, its position in the history of pharmacy development, and its advantages for drug discovery and design, thereby to prove that the attractor method can indeed become the next major drug development method. In addition, it provides a comprehensive description about the concept of attractor, the pipeline of attractor analysis, the common methods of each process and its research progress, so as to provide a macroscopic framework and optional methods and tools for the follow-up researchers.
Collapse
Affiliation(s)
- Xucan Hou
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Meng Li
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Congmin Jia
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Xianbao Zhang
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| | - Yun Wang
- Department of Traditional Chinese Medicine Information Fusion and Utilization, Beijing University of Chinese Medicine, Beijing, People's Republic of China
| |
Collapse
|
29
|
Abstract
Complex disease such as cancer is often caused by genetic mutations that eventually alter the signal flow in the intra-cellular signaling network and result in different cell fate. Therefore, it is crucial to identify control targets that can most effectively block such unwanted signal flow. For this purpose, systems biological analysis provides a useful framework, but mathematical modeling of complicated signaling networks requires massive time-series measurements of signaling protein activity levels for accurate estimation of kinetic parameter values or regulatory logics. Here, we present a novel method, called SFC (Signal Flow Control), for identifying control targets without the information of kinetic parameter values or regulatory logics. Our method requires only the structural information of a signaling network and is based on the topological estimation of signal flow through the network. SFC will be particularly useful for a large-scale signaling network to which parameter estimation or inference of regulatory logics is no longer applicable in practice. The identified control targets have significant implication in drug development as they can be putative drug targets.
Collapse
|
30
|
Trosset JY, Cavé C. In Silico Target Druggability Assessment: From Structural to Systemic Approaches. Methods Mol Biol 2019; 1953:63-88. [PMID: 30912016 DOI: 10.1007/978-1-4939-9145-7_5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
This chapter will focus on today's in silico direct and indirect approaches to assess therapeutic target druggability. The direct approach tries to infer from the 3D structure the capacity of the target protein to bind small molecule in order to modulate its biological function. Algorithms to recognize and characterize the quality of the ligand interaction sites whether within buried protein cavities or within large protein-protein interface will be reviewed in the first part of the paper. In the case a ligand-binding site is already identified, indirect aspects of target druggability can be assessed. These indirect approaches focus first on target promiscuity and the potential difficulties in developing specific drugs. It is based on large-scale comparison of protein-binding sites. The second aspect concerns the capacity of the target to induce resistant pathway once it is inhibited or activated by a drug. The emergence of drug-resistant pathways can be assessed through systemic analysis of biological networks implementing metabolism and/or cell regulation signaling.
Collapse
Affiliation(s)
| | - Christian Cavé
- BioCIS UFR Pharmacie UMR CNRS 8076, Université Paris Saclay, Orsay, France
| |
Collapse
|
31
|
Ozturk K, Dow M, Carlin DE, Bejar R, Carter H. The Emerging Potential for Network Analysis to Inform Precision Cancer Medicine. J Mol Biol 2018; 430:2875-2899. [PMID: 29908887 PMCID: PMC6097914 DOI: 10.1016/j.jmb.2018.06.016] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2018] [Revised: 05/30/2018] [Accepted: 06/06/2018] [Indexed: 12/19/2022]
Abstract
Precision cancer medicine promises to tailor clinical decisions to patients using genomic information. Indeed, successes of drugs targeting genetic alterations in tumors, such as imatinib that targets BCR-ABL in chronic myelogenous leukemia, have demonstrated the power of this approach. However, biological systems are complex, and patients may differ not only by the specific genetic alterations in their tumor, but also by more subtle interactions among such alterations. Systems biology and more specifically, network analysis, provides a framework for advancing precision medicine beyond clinical actionability of individual mutations. Here we discuss applications of network analysis to study tumor biology, early methods for N-of-1 tumor genome analysis, and the path for such tools to the clinic.
Collapse
Affiliation(s)
- Kivilcim Ozturk
- Department of Medicine, Division of Medical Genetics, University of California San Diego, La Jolla, CA 92093, USA; Bioinformatics and Systems Biology Program, University of California San Diego, La Jolla, CA 92093, USA
| | - Michelle Dow
- Department of Medicine, Division of Medical Genetics, University of California San Diego, La Jolla, CA 92093, USA; Bioinformatics and Systems Biology Program, University of California San Diego, La Jolla, CA 92093, USA
| | - Daniel E Carlin
- Department of Medicine, Division of Medical Genetics, University of California San Diego, La Jolla, CA 92093, USA
| | - Rafael Bejar
- Moores Cancer Center, Division of Hematology and Oncology, University of California San Diego, La Jolla, CA 92093, USA
| | - Hannah Carter
- Department of Medicine, Division of Medical Genetics, University of California San Diego, La Jolla, CA 92093, USA; Bioinformatics and Systems Biology Program, University of California San Diego, La Jolla, CA 92093, USA; Moores Cancer Center and Institute for Genomic Medicine, University of California San Diego, La Jolla, CA 92093, USA; CIFAR, MaRS Centre, West Tower, 661 University Ave., Suite 505, Toronto, ON M5G 1M1, Canada.
| |
Collapse
|
32
|
Correia RB, Gates AJ, Wang X, Rocha LM. CANA: A Python Package for Quantifying Control and Canalization in Boolean Networks. Front Physiol 2018; 9:1046. [PMID: 30154728 PMCID: PMC6102667 DOI: 10.3389/fphys.2018.01046] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2018] [Accepted: 07/13/2018] [Indexed: 01/11/2023] Open
Abstract
Logical models offer a simple but powerful means to understand the complex dynamics of biochemical regulation, without the need to estimate kinetic parameters. However, even simple automata components can lead to collective dynamics that are computationally intractable when aggregated into networks. In previous work we demonstrated that automata network models of biochemical regulation are highly canalizing, whereby many variable states and their groupings are redundant (Marques-Pita and Rocha, 2013). The precise charting and measurement of such canalization simplifies these models, making even very large networks amenable to analysis. Moreover, canalization plays an important role in the control, robustness, modularity and criticality of Boolean network dynamics, especially those used to model biochemical regulation (Gates and Rocha, 2016; Gates et al., 2016; Manicka, 2017). Here we describe a new publicly-available Python package that provides the necessary tools to extract, measure, and visualize canalizing redundancy present in Boolean network models. It extracts the pathways most effective in controlling dynamics in these models, including their effective graph and dynamics canalizing map, as well as other tools to uncover minimum sets of control variables.
Collapse
Affiliation(s)
- Rion B. Correia
- School of Informatics, Computing, and Engineering, Indiana University, Bloomington, IN, United States
- CAPES Foundation, Ministry of Education of Brazil, Brasília, Brazil
- Instituto Gulbenkian de Ciência, Oeiras, Portugal
| | - Alexander J. Gates
- Center for Complex Networks Research, Northeastern University, Boston, MA, United States
| | - Xuan Wang
- School of Informatics, Computing, and Engineering, Indiana University, Bloomington, IN, United States
| | - Luis M. Rocha
- School of Informatics, Computing, and Engineering, Indiana University, Bloomington, IN, United States
- Instituto Gulbenkian de Ciência, Oeiras, Portugal
| |
Collapse
|
33
|
Cell signaling heterogeneity is modulated by both cell-intrinsic and -extrinsic mechanisms: An integrated approach to understanding targeted therapy. PLoS Biol 2018. [PMID: 29522507 PMCID: PMC5844524 DOI: 10.1371/journal.pbio.2002930] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
During the last decade, our understanding of cancer cell signaling networks has significantly improved, leading to the development of various targeted therapies that have elicited profound but, unfortunately, short-lived responses. This is, in part, due to the fact that these targeted therapies ignore context and average out heterogeneity. Here, we present a mathematical framework that addresses the impact of signaling heterogeneity on targeted therapy outcomes. We employ a simplified oncogenic rat sarcoma (RAS)-driven mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase-protein kinase B (PI3K-AKT) signaling pathway in lung cancer as an experimental model system and develop a network model of the pathway. We measure how inhibition of the pathway modulates protein phosphorylation as well as cell viability under different microenvironmental conditions. Training the model on this data using Monte Carlo simulation results in a suite of in silico cells whose relative protein activities and cell viability match experimental observation. The calibrated model predicts distributional responses to kinase inhibitors and suggests drug resistance mechanisms that can be exploited in drug combination strategies. The suggested combination strategies are validated using in vitro experimental data. The validated in silico cells are further interrogated through an unsupervised clustering analysis and then integrated into a mathematical model of tumor growth in a homogeneous and resource-limited microenvironment. We assess posttreatment heterogeneity and predict vast differences across treatments with similar efficacy, further emphasizing that heterogeneity should modulate treatment strategies. The signaling model is also integrated into a hybrid cellular automata (HCA) model of tumor growth in a spatially heterogeneous microenvironment. As a proof of concept, we simulate tumor responses to targeted therapies in a spatially segregated tissue structure containing tumor and stroma (derived from patient tissue) and predict complex cell signaling responses that suggest a novel combination treatment strategy. A signaling pathway is a network of molecules in a cell that is typically initiated by stimuli (e.g., microenvironmental cues) acting on receptors and internal signaling molecules to determine cell fate. Signaling pathways in cancer cells are different from those in normal cells, and this difference helps cancer cells to grow and thrive indefinitely. Drugs that target the aberrant signaling pathways in cancer cells (often referred to as targeted therapy) are promising for improving treatment outcomes of many different cancers in patients. However, most patients eventually develop resistance to these drugs. Resistance may already be present in the tumor or may emerge via mutation or via microenvironmental mediation. Tumor heterogeneity, which is characterized by subtle or dramatic differences among tumor cells, plays a key role in the development of drug resistance. Some tumor cells respond well to therapy, while others may adapt to the stress induced by the drug within the microenvironment. Moreover, removal of drug-sensitive cells may result in the competitive release of drug-resistant cells. Here, we present mathematical models to assess the impact of heterogeneity in signaling pathways within tumor cells on the outcomes of targeted therapy. We consider a simplified version of two well-known signaling pathways that modulate the growth of lung cancer cells. By using different targeted therapies, we quantify the effect of pathway inhibition on protein activity and cell viability and developed a mathematical model of the network, which is trained to reproduce these data and to develop a panel of heterogeneous in silico cells. The model predicts potential mechanisms of drug resistance and proposes combination therapies that are effective across the panel. We validate these combination therapies experimentally using the lung cancer cells and integrated the in silico cells into a computational lung tissue model that explicitly captures the microenvironment of lung cancer. Our results suggest that heterogeneity in both the tumor and microenvironment impacts treatment response in different ways and suggest a novel combination therapy for a better response.
Collapse
|