Simplification of reversible Markov chains by removal of states with low equilibrium occupancy.

We present a practical method for simplifying Markov chains on a potentially large state space when detailed balance holds. A simple and transparent technique is introduced to remove states with low equilibrium occupancy. The resulting system has fewer parameters. The resulting effective rates between the remaining nodes give dynamics identical to the original system's except on very fast timescales. This procedure amounts to using separation of timescales to neglect small capacitance nodes in a network of resistors and capacitors. We illustrate the technique by simplifying various reaction networks, including transforming an acyclic four-node network to a three-node cyclic network. For a reaction step in which a ligand binds, the law of mass action implies a forward rate proportional to ligand concentration. The effective rates in the simplified network are found to be rational functions of ligand concentration.

[1]  J. Changeux,et al.  ON THE NATURE OF ALLOSTERIC TRANSITIONS: A PLAUSIBLE MODEL. , 1965, Journal of molecular biology.

[2]  John Maindonald,et al.  Data Analysis and Graphics Using R: An Example-based Approach (Cambridge Series in Statistical and Probabilistic Mathematics) , 2003 .

[3]  D. Gillespie A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[4]  W. S. Hlavacek,et al.  On imposing detailed balance in complex reaction mechanisms. , 2006, Biophysical journal.

[5]  H. Akaike A new look at the statistical model identification , 1974 .

[6]  S. Meyn,et al.  Phase transitions and metastability in Markovian and molecular systems , 2004 .

[7]  Robust algorithm for random resistor networks using hierarchical domain structure , 2006 .

[8]  Michele Benzi Numerical Solution of Markov Chains , 2011, Numer. Linear Algebra Appl..

[9]  A. N. Shiryayev,et al.  On The Theory of Markov Chains , 1992 .

[10]  Rhjm Ralph Otten,et al.  Planarization by transformation , 1973 .

[11]  Ghanim Ullah,et al.  A data-driven model of a modal gated ion channel: The inositol 1,4,5-trisphosphate receptor in insect Sf9 cells , 2012, The Journal of general physiology.

[12]  S. Akers The Use of Wye-Delta Transformations in Network Simplification , 1960 .

[13]  David Colquhoun,et al.  Agonist‐activated ion channels , 2006, British journal of pharmacology.

[14]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[15]  L. Onsager Reciprocal Relations in Irreversible Processes. II. , 1931 .

[16]  J. Pearson,et al.  Using independent open-to-closed transitions to simplify aggregated Markov models of ion channel gating kinetics. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[17]  Sean P. Meyn,et al.  An information-theoretic framework to aggregate a Markov chain , 2009, 2009 American Control Conference.