Selective Control of the Apoptosis Signaling Network in Heterogeneous Cell Populations

Background Selective control in a population is the ability to control a member of the population while leaving the other members relatively unaffected. The concept of selective control is developed using cell death or apoptosis in heterogeneous cell populations as an example. Control of apoptosis is essential in a variety of therapeutic environments, including cancer where cancer cell death is a desired outcome and Alzheimer's disease where neuron survival is the desired outcome. However, in both cases these responses must occur with minimal response in other cells exposed to treatment; that is, the response must be selective. Methodology and Principal Findings Apoptosis signaling in heterogeneous cells is described by an ensemble of gene networks with identical topology but different link strengths. Selective control depends on the statistics of signaling in the ensemble of networks, and we analyze the effects of superposition, non-linearity and feedback on these statistics. Parallel pathways promote normal statistics while series pathways promote skew distributions, which in the most extreme cases become log-normal. We also show that feedback and non-linearity can produce bimodal signaling statistics, as can discreteness and non-linearity. Two methods for optimizing selective control are presented. The first is an exhaustive search method and the second is a linear programming based approach. Though control of a single gene in the signaling network yields little selectivity, control of a few genes typically yields higher levels of selectivity. The statistics of gene combinations susceptible to selective control in heterogeneous apoptosis networks is studied and is used to identify general control strategies. Conclusions and Significance We have explored two methods for the study of selectivity in cell populations. The first is an exhaustive search method limited to three node perturbations. The second is an effective linear model, based on interpolation of single node sensitivity, in which the selective combinations can be found by linear programming optimization. We found that selectivity is promoted by acting on the least sensitive nodes in the case of weak populations, while selective control of robust populations is optimized through perturbations of more sensitive nodes. High throughput experiments with heterogeneous cell lines could be designed in an analogous manner, with the further possibility of incorporating the selectivity optimization process into a closed-loop control system.

[1]  Jas. Middlemass,et al.  Archives of Neurology , 1903, British medical journal.

[2]  S. Lowen The Biophysical Journal , 1960, Nature.

[3]  E. Hill Journal of Theoretical Biology , 1961, Nature.

[4]  M. Tanner Trends in Biochemical Sciences , 1982 .

[5]  P. Sperryn,et al.  Blood. , 1989, British journal of sports medicine.

[6]  Physical Review Letters 63 , 1989 .

[7]  Ronald L. Rivest,et al.  Introduction to Algorithms , 1990 .

[8]  Stuart A. Kauffman,et al.  The origins of order , 1993 .

[9]  AC Tose Cell , 1993, Cell.

[10]  Arantxa Etxeverria The Origins of Order , 1993 .

[11]  J. Herskowitz,et al.  Proceedings of the National Academy of Sciences, USA , 1996, Current Biology.

[12]  James N. Weiss The Hill equation revisited: uses and misuses , 1997, FASEB journal : official publication of the Federation of American Societies for Experimental Biology.

[13]  S. Snyder,et al.  Proceedings of the National Academy of Sciences , 1999 .

[14]  B. Palsson,et al.  Regulation of gene expression in flux balance models of metabolism. , 2001, Journal of theoretical biology.

[15]  Scott W. Lowe,et al.  Apoptosis A Link between Cancer Genetics and Chemotherapy , 2002, Cell.

[16]  Richard Mohs,et al.  Caspase gene expression in the brain as a function of the clinical progression of Alzheimer disease. , 2003, Archives of neurology.

[17]  H. Othmer,et al.  The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster. , 2003, Journal of theoretical biology.

[18]  R. Weiss,et al.  Optimizing genetic circuits by global sensitivity analysis. , 2004, Biophysical journal.

[19]  Q. Ouyang,et al.  The yeast cell-cycle network is robustly designed. , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Herschel Rabitz,et al.  Optimal identification of biochemical reaction networks. , 2004, Biophysical journal.

[21]  Xiao-Jiang Feng,et al.  Closed-Loop Learning Control of Bio-Networks , 2004, J. Comput. Biol..

[22]  K. Lu Pinning down cell signaling, cancer and Alzheimer's disease. , 2004, Trends in biochemical sciences.

[23]  Alessandro Flammini,et al.  Characterization and modeling of protein–protein interaction networks , 2005 .

[24]  Joshua E. S. Socolar,et al.  Modeling pathways of differentiation in genetic regulatory networks with Boolean networks , 2005, Complex..

[25]  Barbara Drossel Number of attractors in random Boolean networks. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  B. Winblad,et al.  Biomarkers for apoptosis in Alzheimer's disease , 2005, International journal of geriatric psychiatry.

[27]  F. Allgöwer,et al.  Robustness properties of apoptosis models with respect to parameter variations and intrinsic noise. , 2005, Systems biology.

[28]  John Calvin Reed,et al.  Apoptosis-based therapies for hematologic malignancies. , 2005, Blood.

[29]  Nicolas E. Buchler,et al.  Nonlinear protein degradation and the function of genetic circuits. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[30]  Reinhard Lipowsky,et al.  Dynamic pattern evolution on scale-free networks. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[31]  Gabriel S. Eichler,et al.  Cell fates as high-dimensional attractor states of a complex gene regulatory network. , 2005, Physical review letters.

[32]  Ilya Shmulevich,et al.  Eukaryotic cells are dynamically ordered or critical but not chaotic. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[33]  J. Urry Complexity , 2006, Interpreting Art.

[34]  G B Ermentrout,et al.  Bistability in apoptosis: roles of bax, bcl-2, and mitochondrial permeability transition pores. , 2006, Biophysical journal.

[35]  J. Houghton,et al.  Molecular Targeting of Cell Death Signal Transduction Pathways in Cancer , 2006 .

[36]  Erwin P. Gianchandani,et al.  Flux balance analysis in the era of metabolomics , 2006, Briefings Bioinform..

[37]  B. Kholodenko Cell-signalling dynamics in time and space , 2006, Nature Reviews Molecular Cell Biology.

[38]  Bernhard O. Palsson,et al.  Iterative Reconstruction of Transcriptional Regulatory Networks: An Algorithmic Approach , 2006, PLoS Comput. Biol..

[39]  Leon Glass,et al.  Reverse Engineering the Gap Gene Network of Drosophila melanogaster , 2006, PLoS Comput. Biol..

[40]  P. Sorger,et al.  Systems biology and combination therapy in the quest for clinical efficacy , 2006, Nature chemical biology.

[41]  D. Lauffenburger,et al.  Physicochemical modelling of cell signalling pathways , 2006, Nature Cell Biology.

[42]  A. Levitzki,et al.  Signal Transduction Therapy for Cancer - Whither Now? , 2006 .

[43]  Madhukar S. Dasika,et al.  A computational framework for the topological analysis and targeted disruption of signal transduction networks. , 2006, Biophysical journal.

[44]  Monica L. Mo,et al.  Global reconstruction of the human metabolic network based on genomic and bibliomic data , 2007, Proceedings of the National Academy of Sciences.

[45]  Werner Dubitzky,et al.  Briefings in bioinformatics. , 2009, Briefings in bioinformatics.

[46]  October I Physical Review Letters , 2022 .