A Data-driven Optimization Method for Coarse-graining Gene Regulatory Networks

Overview

Journal iScience

Publisher Cell Press

Date 2023 Jan 26

PMID 36698721

Authors

Cristian Caranica

Mingyang Lu

Affiliations

Soon will be listed here.

Abstract

One major challenge in systems biology is to understand how various genes in a gene regulatory network (GRN) collectively perform their functions and control network dynamics. This task becomes extremely hard to tackle in the case of large networks with hundreds of genes and edges, many of which have redundant regulatory roles and functions. The existing methods for model reduction usually require the detailed mathematical description of dynamical systems and their corresponding kinetic parameters, which are often not available. Here, we present a data-driven method for coarse-graining large GRNs, named SacoGraci, using ensemble-based mathematical modeling, dimensionality reduction, and gene circuit optimization by Markov Chain Monte Carlo methods. SacoGraci requires network topology as the only input and is robust against errors in GRNs. We benchmark and demonstrate its usage with synthetic, literature-based, and bioinformatics-derived GRNs. We hope SacoGraci will enhance our ability to model the gene regulation of complex biological systems.

References

Thiele I, Palsson B . A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010; 5(1):93-121. PMC: 3125167. DOI: 10.1038/nprot.2009.203. View

Lourenco D, Brauner R, Rybczynska M, Nihoul-Fekete C, McElreavey K, Bashamboo A . Loss-of-function mutation in GATA4 causes anomalies of human testicular development. Proc Natl Acad Sci U S A. 2011; 108(4):1597-602. PMC: 3029689. DOI: 10.1073/pnas.1010257108. View

Jia D, Jolly M, Tripathi S, den Hollander P, Huang B, Lu M . Distinguishing mechanisms underlying EMT tristability. Cancer Converg. 2018; 1(1):2. PMC: 5876698. DOI: 10.1186/s41236-017-0005-8. View

Yang Y, Workman S, Wilson M . The molecular pathways underlying early gonadal development. J Mol Endocrinol. 2018; . DOI: 10.1530/JME-17-0314. View

Prescott T, Papachristodoulou A . Layered decomposition for the model order reduction of timescale separated biochemical reaction networks. J Theor Biol. 2014; 356:113-22. DOI: 10.1016/j.jtbi.2014.04.007. View

Lu M, Poon B, Ma J . A New Method for Coarse-Grained Elastic Normal-Mode Analysis. J Chem Theory Comput. 2011; 2(3):464-471. PMC: 3133770. DOI: 10.1021/ct050307u. View

Kulkarni V, Arastoo R, Bhat A, Subramanian K, Kothare M, Riedel M . Gene regulatory network modeling using literature curated and high throughput data. Syst Synth Biol. 2013; 6(3-4):69-77. PMC: 3528886. DOI: 10.1007/s11693-012-9100-4. View

Nieto M, Huang R, Jackson R, Thiery J . EMT: 2016. Cell. 2016; 166(1):21-45. DOI: 10.1016/j.cell.2016.06.028. View

Ohnesorg T, Vilain E, Sinclair A . The genetics of disorders of sex development in humans. Sex Dev. 2014; 8(5):262-72. DOI: 10.1159/000357956. View

10.

Feist A, Herrgard M, Thiele I, Reed J, Palsson B . Reconstruction of biochemical networks in microorganisms. Nat Rev Microbiol. 2009; 7(2):129-43. PMC: 3119670. DOI: 10.1038/nrmicro1949. View

11.

Tian X, Zhang H, Xing J . Coupled reversible and irreversible bistable switches underlying TGFβ-induced epithelial to mesenchymal transition. Biophys J. 2013; 105(4):1079-89. PMC: 3752104. DOI: 10.1016/j.bpj.2013.07.011. View

12.

Pranzatelli T, Michael D, Chiorini J . ATAC2GRN: optimized ATAC-seq and DNase1-seq pipelines for rapid and accurate genome regulatory network inference. BMC Genomics. 2018; 19(1):563. PMC: 6069842. DOI: 10.1186/s12864-018-4943-z. View

13.

Chauhan L, Ram U, Hari K, Jolly M . Topological signatures in regulatory network enable phenotypic heterogeneity in small cell lung cancer. Elife. 2021; 10. PMC: 8012062. DOI: 10.7554/eLife.64522. View

14.

Matsumoto H, Kiryu H, Furusawa C, Ko M, Ko S, Gouda N . SCODE: an efficient regulatory network inference algorithm from single-cell RNA-Seq during differentiation. Bioinformatics. 2017; 33(15):2314-2321. PMC: 5860123. DOI: 10.1093/bioinformatics/btx194. View

15.

Gribble S, Kim H, Bonner J, Wang X, Dorsky R . Tcf3 inhibits spinal cord neurogenesis by regulating sox4a expression. Development. 2009; 136(5):781-9. PMC: 2685945. DOI: 10.1242/dev.027995. View

16.

Atilgan A, Durell S, Jernigan R, Demirel M, Keskin O, Bahar I . Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys J. 2001; 80(1):505-15. PMC: 1301252. DOI: 10.1016/S0006-3495(01)76033-X. View

17.

Zhang J, Tian X, Zhang H, Teng Y, Li R, Bai F . TGF-β-induced epithelial-to-mesenchymal transition proceeds through stepwise activation of multiple feedback loops. Sci Signal. 2014; 7(345):ra91. DOI: 10.1126/scisignal.2005304. View

18.

Kmiecik S, Gront D, Kolinski M, Wieteska L, Dawid A, Kolinski A . Coarse-Grained Protein Models and Their Applications. Chem Rev. 2016; 116(14):7898-936. DOI: 10.1021/acs.chemrev.6b00163. View

19.

Karacosta L, Anchang B, Ignatiadis N, Kimmey S, Benson J, Shrager J . Mapping lung cancer epithelial-mesenchymal transition states and trajectories with single-cell resolution. Nat Commun. 2019; 10(1):5587. PMC: 6898514. DOI: 10.1038/s41467-019-13441-6. View

20.

Margolin A, Nemenman I, Basso K, Wiggins C, Stolovitzky G, Dalla Favera R . ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinformatics. 2006; 7 Suppl 1:S7. PMC: 1810318. DOI: 10.1186/1471-2105-7-S1-S7. View