Untangling the Hairball: Fitness-Based Asymptotic Reduction of Biological Networks.

Complex mathematical models of interaction networks are routinely used for prediction in systems biology. However, it is difficult to reconcile network complexities with a formal understanding of their behavior. Here, we propose a simple procedure (called ϕ¯) to reduce biological models to functional submodules, using statistical mechanics of complex systems combined with a fitness-based approach inspired by in silico evolution. The ϕ¯ algorithm works by putting parameters or combination of parameters to some asymptotic limit, while keeping (or slightly improving) the model performance, and requires parameter symmetry breaking for more complex models. We illustrate ϕ¯ on biochemical adaptation and on different models of immune recognition by T cells. An intractable model of immune recognition with close to a hundred individual transition rates is reduced to a simple two-parameter model. The ϕ¯ algorithm extracts three different mechanisms for early immune recognition, and automatically discovers similar functional modules in different models of the same process, allowing for model classification and comparison. Our procedure can be applied to biological networks based on rate equations using a fitness function that quantifies phenotypic performance.

[1]  H. Metzger,et al.  An unusual mechanism for ligand antagonism. , 1998, Science.

[2]  Paul François,et al.  Phenotypic model for early T-cell activation displaying sensitivity, specificity, and antagonism , 2013, Proceedings of the National Academy of Sciences.

[3]  James R Faeder,et al.  Stochastic effects and bistability in T cell receptor signaling. , 2008, Journal of theoretical biology.

[4]  Jonathan R. Karr,et al.  A Whole-Cell Computational Model Predicts Phenotype from Genotype , 2012, Cell.

[5]  T. McKeithan,et al.  Kinetic proofreading in T-cell receptor signal transduction. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[6]  Konrad Paul Kording,et al.  Could a Neuroscientist Understand a Microprocessor? , 2016, bioRxiv.

[7]  C. Janeway,et al.  Cross-antagonism of a T cell clone expressing two distinct T cell receptors. , 1999, Immunity.

[8]  Paul François,et al.  Principles of adaptive sorting revealed by in silico evolution. , 2013, Physical review letters.

[9]  L. F. Abbott,et al.  Generating Coherent Patterns of Activity from Chaotic Neural Networks , 2009, Neuron.

[10]  J. Sethna,et al.  Parameter Space Compression Underlies Emergent Theories and Predictive Models , 2013, Science.

[11]  Bryan C. Daniels,et al.  Perspective: Sloppiness and emergent theories in physics, biology, and beyond. , 2015, The Journal of chemical physics.

[12]  J. Doyle,et al.  Robust perfect adaptation in bacterial chemotaxis through integral feedback control. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[13]  Jeremy Gunawardena,et al.  Models in biology: ‘accurate descriptions of our pathetic thinking’ , 2014, BMC Biology.

[14]  Mark K Transtrum,et al.  Model reduction by manifold boundaries. , 2014, Physical review letters.

[15]  Ronald N Germain,et al.  Modeling T Cell Antigen Discrimination Based on Feedback Control of Digital ERK Responses , 2005, PLoS biology.

[16]  R. Germain,et al.  Variability and Robustness in T Cell Activation from Regulated Heterogeneity in Protein Levels , 2008, Science.

[17]  Ofer Feinerman,et al.  Quantitative challenges in understanding ligand discrimination by αβ T cells , 2008 .

[18]  Christopher R. Myers,et al.  Universally Sloppy Parameter Sensitivities in Systems Biology Models , 2007, PLoS Comput. Biol..

[19]  Peng Qiu,et al.  Bridging Mechanistic and Phenomenological Models of Complex Biological Systems , 2015, PLoS Comput. Biol..

[20]  Ronald D. Vale,et al.  A DNA-Based T Cell Receptor Reveals a Role for Receptor Clustering in Ligand Discrimination , 2016, Cell.

[21]  James Sharpe,et al.  An atlas of gene regulatory networks reveals multiple three-gene mechanisms for interpreting morphogen gradients , 2010, Molecular systems biology.

[22]  James G. King,et al.  Reconstruction and Simulation of Neocortical Microcircuitry , 2015, Cell.

[23]  Paul François,et al.  Evolving phenotypic networks in silico. , 2014, Seminars in cell & developmental biology.

[24]  Omer Dushek,et al.  Phenotypic models of T cell activation , 2014, Nature Reviews Immunology.

[25]  Paul François,et al.  The Case for Absolute Ligand Discrimination: Modeling Information Processing and Decision by Immune T Cells , 2015, bioRxiv.

[26]  S. Shen-Orr,et al.  Network motifs: simple building blocks of complex networks. , 2002, Science.

[27]  Olivier Tassy,et al.  Evolutionary plasticity of segmentation clock networks , 2011, Development.

[28]  K. S. Brown,et al.  Statistical mechanical approaches to models with many poorly known parameters. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[29]  Arthur D Lander,et al.  The edges of understanding , 2010, BMC Biology.

[30]  Bryan C. Daniels,et al.  Automated adaptive inference of phenomenological dynamical models , 2015, Nature Communications.

[31]  Mark K. Transtrum,et al.  Information topology identifies emergent model classes , 2014 .

[32]  Paul François,et al.  Phenotypic spandrel: absolute discrimination and ligand antagonism , 2015, bioRxiv.

[33]  K. H. Lee,et al.  The statistical mechanics of complex signaling networks: nerve growth factor signaling , 2004, Physical biology.

[34]  F. Kondrashov,et al.  The evolution of gene duplications: classifying and distinguishing between models , 2010, Nature Reviews Genetics.

[35]  W. Lim,et al.  Defining Network Topologies that Can Achieve Biochemical Adaptation , 2009, Cell.

[36]  Eduardo Sontag,et al.  Multiple steady states and the form of response functions to antigen in a model for the initiation of T-cell activation , 2017, Royal Society Open Science.

[37]  E. Siggia,et al.  Predicting embryonic patterning using mutual entropy fitness and in silico evolution , 2010, Development.

[38]  Naomi S. Altman,et al.  Points of Significance: Model selection and overfitting , 2016, Nature Methods.

[39]  Paul François,et al.  A case study of evolutionary computation of biochemical adaptation , 2008, Physical biology.

[40]  Francis Corson,et al.  Geometry, epistasis, and developmental patterning , 2012, Proceedings of the National Academy of Sciences.

[41]  Jonathon Howard,et al.  Drawing an elephant with four complex parameters , 2010 .

[42]  Bryan C. Daniels,et al.  Efficient Inference of Parsimonious Phenomenological Models of Cellular Dynamics Using S-Systems and Alternating Regression , 2014, PloS one.