Prospects for Declarative Mathematical Modeling of Complex Biological Systems

Declarative modeling uses symbolic expressions to represent models. With such expressions, one can formalize high-level mathematical computations on models that would be difficult or impossible to perform directly on a lower-level simulation program, in a general-purpose programming language. Examples of such computations on models include model analysis, relatively general-purpose model reduction maps, and the initial phases of model implementation, all of which should preserve or approximate the mathematical semantics of a complex biological model. The potential advantages are particularly relevant in the case of developmental modeling, wherein complex spatial structures exhibit dynamics at molecular, cellular, and organogenic levels to relate genotype to multicellular phenotype. Multiscale modeling can benefit from both the expressive power of declarative modeling languages and the application of model reduction methods to link models across scale. Based on previous work, here we define declarative modeling of complex biological systems by defining the operator algebra semantics of an increasingly powerful series of declarative modeling languages including reaction-like dynamics of parameterized and extended objects; we define semantics-preserving implementation and semantics-approximating model reduction transformations; and we outline a “meta-hierarchy” for organizing declarative models and the mathematical methods that can fruitfully manipulate them.

[1]  Ben Milder Smith, et al , 1992 .

[2]  Anders Logg,et al.  Unified form language: A domain-specific language for weak formulations of partial differential equations , 2012, TOMS.

[3]  Przemyslaw Prusinkiewicz,et al.  Animation of plant development , 1993, SIGGRAPH.

[4]  G. Venugopala Reddy,et al.  Cell-Based Model of the Generation and Maintenance of the Shape and Structure of the Multilayered Shoot Apical Meristem of Arabidopsis thaliana , 2018, Bulletin of Mathematical Biology.

[5]  Bruce E. Shapiro,et al.  Pycellerator: an arrow-based reaction-like modelling language for biological simulations , 2016, Bioinform..

[6]  M. Maggioni,et al.  Manifold parametrizations by eigenfunctions of the Laplacian and heat kernels , 2008, Proceedings of the National Academy of Sciences.

[7]  Klaus Gärtner,et al.  Meshing Piecewise Linear Complexes by Constrained Delaunay Tetrahedralizations , 2005, IMR.

[8]  W. Ebeling Stochastic Processes in Physics and Chemistry , 1995 .

[9]  Max Welling,et al.  Semi-Supervised Classification with Graph Convolutional Networks , 2016, ICLR.

[10]  Mark Kirkpatrick,et al.  Compensatory Drift and the Evolutionary Dynamics of Dosage-Sensitive Duplicate Genes , 2015, Genetics.

[11]  Nikolaus Grigorieff,et al.  Severing enzymes amplify microtubule arrays through lattice GTP-tubulin incorporation , 2018, Science.

[12]  David Joyner,et al.  Open source computer algebra systems: SymPy , 2012, ACCA.

[13]  Terry Winograd,et al.  FRAME REPRESENTATIONS AND THE DECLARATIVE/PROCEDURAL CONTROVERSY , 1975 .

[14]  Bard Ermentrout,et al.  Simplifying and reducing complex models , 2002 .

[15]  Pierre Vandergheynst,et al.  Wavelets on Graphs via Spectral Graph Theory , 2009, ArXiv.

[16]  Eric Mjolsness,et al.  Towards Measurable Types for Dynamical Process Modeling Languages , 2010, MFPS.

[17]  E. Mjolsness,et al.  Analysis of cell division patterns in the Arabidopsis shoot apical meristem , 2015, Proceedings of the National Academy of Sciences.

[18]  Bruce E. Shapiro,et al.  Modeling the organization of the WUSCHEL expression domain in the shoot apical meristem , 2005, ISMB.

[19]  John W. Backus,et al.  Can programming be liberated from the von Neumann style?: a functional style and its algebra of programs , 1978, CACM.

[20]  Eric Mjolsness,et al.  Dependency diagrams and graph-constrained correlation dynamics: new systems for probabilistic graphical modeling , 2012 .

[21]  Darren Engwirda,et al.  Conforming restricted Delaunay mesh generation for piecewise smooth complexes , 2016, ArXiv.

[22]  William S. Hlavacek,et al.  BioNetGen: software for rule-based modeling of signal transduction based on the interactions of molecular domains , 2004, Bioinform..

[23]  Bruce E. Shapiro,et al.  An enzyme mechanism language for the mathematical modeling of metabolic pathways , 2005, Bioinform..

[24]  D. Mattis,et al.  The uses of quantum field theory in diffusion-limited reactions , 1998 .

[25]  E. Deinum,et al.  Efficient event-driven simulations shed new light on microtubule organization in the plant cortical array , 2014, Front. Physics.

[26]  Xiaohui Xie,et al.  Parameter inference for discretely observed stochastic kinetic models using stochastic gradient descent , 2010, BMC Systems Biology.

[27]  Roeland M. H. Merks,et al.  Adapting a Plant Tissue Model to Animal Development: Introducing Cell Sliding into VirtualLeaf , 2018, Bulletin of Mathematical Biology.

[28]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[29]  Eric Mjolsness,et al.  Learning dynamic Boltzmann distributions as reduced models of spatial chemical kinetics. , 2018, The Journal of chemical physics.

[30]  C. Schwechheimer,et al.  DiSUMO-LIKE Interacts with RNA-Binding Proteins and Affects Cell-Cycle Progression during Maize Embryogenesis , 2018, Current Biology.

[31]  Luca Cardelli,et al.  On process rate semantics , 2008, Theor. Comput. Sci..

[32]  Eric Mjolsness,et al.  Mathematics of small stochastic reaction networks: a boundary layer theory for eigenstate analysis. , 2013, The Journal of chemical physics.

[33]  L. Peliti Path integral approach to birth-death processes on a lattice , 1985 .

[34]  Brendan J. Frey,et al.  Extending Factor Graphs so as to Unify Directed and Undirected Graphical Models , 2002, UAI.

[35]  L. Wolpert Positional information and the spatial pattern of cellular differentiation. , 1969, Journal of theoretical biology.

[36]  M. Doi Second quantization representation for classical many-particle system , 1976 .

[37]  Luidnel Maignan,et al.  Global Graph Transformations , 2015, GCM@ICGT.

[38]  David M. Mount,et al.  A point-placement strategy for conforming Delaunay tetrahedralization , 2000, SODA '00.

[39]  Justin B. Kinney,et al.  Modeling multi-particle complexes in stochastic chemical systems , 2016 .

[40]  F. Nédélec,et al.  Collective Langevin dynamics of flexible cytoskeletal fibers , 2007, 0903.5178.

[41]  J. Glimm,et al.  Quantum Physics: A Functional Integral Point of View , 1981 .

[42]  Anil N. Hirani,et al.  Discrete exterior calculus , 2005, math/0508341.

[43]  J S Smith,et al.  ANI-1: an extensible neural network potential with DFT accuracy at force field computational cost , 2016, Chemical science.

[44]  Alexander Rand,et al.  Collars and Intestines: Practical Conforming Delaunay Refinement , 2009, IMR.

[45]  E. Mjolsness,et al.  An auxin-driven polarized transport model for phyllotaxis , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[46]  Eric Mjolsness,et al.  Stochastic parameterized grammars: formalization, inference and modeling applications , 2009 .

[47]  Raymond Wightman,et al.  Severing at sites of microtubule crossover contributes to microtubule alignment in cortical arrays. , 2007, The Plant journal : for cell and molecular biology.

[48]  Bandan Chakrabortty,et al.  A Plausible Microtubule-B ased Mechanism for Cell Division Orientation in Plant Embryogenesis Graphical , 2018 .

[49]  David H. Sharp,et al.  A connectionist model of development. , 1991, Journal of theoretical biology.

[50]  Gerald J. Sussman,et al.  Structure and Interpretation of Computer Programs, Second Edition , 1996 .

[51]  D. Arnold,et al.  Finite element exterior calculus: From hodge theory to numerical stability , 2009, 0906.4325.

[52]  Eric Mjolsness,et al.  Stochastic process semantics for dynamical grammars , 2006, Annals of Mathematics and Artificial Intelligence.

[53]  Bruce E. Shapiro,et al.  Cellerator: extending a computer algebra system to include biochemical arrows for signal transduction simulations , 2003, Bioinform..

[54]  Bruce E. Shapiro,et al.  Using cellzilla for plant growth simulations at the cellular level , 2013, Front. Plant Sci..

[55]  M. Doi Stochastic theory of diffusion-controlled reaction , 1976 .

[56]  Eric Mjolsness,et al.  A Hierarchical Exact Accelerated Stochastic Simulation Algorithm , 2012, The Journal of chemical physics.

[57]  E. Mjolsness,et al.  Topological object types for morphodynamic modeling languages , 2012, 2012 IEEE 4th International Symposium on Plant Growth Modeling, Simulation, Visualization and Applications.

[58]  Uri Stav,et al.  Non-Linear Index Coding Outperforming the Linear Optimum , 2007, FOCS.

[59]  Jean-Louis Giavitto,et al.  Topological rewriting and the geometrization of programming , 2008 .

[60]  Eric Mjolsness,et al.  Time-ordered product expansions for computational stochastic system biology , 2012, Physical biology.

[61]  Eric Mjolsness,et al.  Model reduction for stochastic CaMKII reaction kinetics in synapses by graph-constrained correlation dynamics. , 2015, Physical biology.

[62]  Eric Mjolsness,et al.  Stochastic Process Semantics for Dynamical Grammar Syntax: An Overview , 2005, AI&M.

[63]  François Fages,et al.  Abstract interpretation and types for systems biology , 2008, Theor. Comput. Sci..

[64]  Richard C. Moore,et al.  The early stages of duplicate gene evolution , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[65]  Olivier Michel,et al.  Declarative modeling of a neurulation-like process , 2007, Biosyst..

[66]  V. A. Likhoshvai,et al.  Combined in silico/in vivo analysis of mechanisms providing for root apical meristem self-organization and maintenance , 2012, Annals of botany.

[67]  J. Lambek,et al.  Introduction to higher order categorical logic , 1986 .

[68]  Gordon D. Plotkin,et al.  The origins of structural operational semantics , 2004, J. Log. Algebraic Methods Program..

[69]  A. Ludlow,et al.  Atomic clock performance enabling geodesy below the centimetre level , 2018, Nature.

[70]  T Poggio,et al.  Regularization Algorithms for Learning That Are Equivalent to Multilayer Networks , 1990, Science.

[71]  Christopher D. Thompson-Walsh,et al.  Graphs, Rewriting and Pathway Reconstruction for Rule-Based Models , 2012, FSTTCS.

[72]  Jean Vaillant,et al.  Differential analysis on stratified spaces , 2003 .

[73]  Hartmut Ehrig,et al.  Fundamentals of Algebraic Graph Transformation , 2006, Monographs in Theoretical Computer Science. An EATCS Series.

[74]  Mathieu Coppey,et al.  Stochastic theory of diffusion-controlled reactions , 2003 .

[75]  Alberto Martelli,et al.  An Efficient Unification Algorithm , 1982, TOPL.

[76]  Sidney L. Shaw,et al.  Sustained Microtubule Treadmilling in Arabidopsis Cortical Arrays , 2003, Science.

[77]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[78]  Vincent Danos,et al.  Stochastic mechanics of graph rewriting , 2016, 2016 31st Annual ACM/IEEE Symposium on Logic in Computer Science (LICS).

[79]  Jean-Louis Giavitto,et al.  MGS: a Programming Language for the Transformations of Topological Collections , 2011 .

[80]  Adrienne Roeder,et al.  Tessellations and Pattern Formation in Plant Growth and Development , 2012, 1209.2937.

[81]  Alain Pumir,et al.  Strain- or Stress-Sensing in Mechanochemical Patterning by the Phytohormone Auxin , 2019, bioRxiv.

[82]  Kei-Hoi Cheung,et al.  BioPAX – A community standard for pathway data sharing , 2010, Nature Biotechnology.

[83]  Darren J Wilkinson,et al.  Bayesian parameter inference for stochastic biochemical network models using particle Markov chain Monte Carlo , 2011, Interface Focus.

[84]  Max Gunzburger,et al.  POD and CVT-based reduced-order modeling of Navier-Stokes flows , 2006 .

[85]  Eric Mjolsness,et al.  Multilevel Artificial Neural Network Training for Spatially Correlated Learning , 2018, SIAM J. Sci. Comput..

[86]  Petros Koumoutsakos,et al.  An exact accelerated stochastic simulation algorithm. , 2009, The Journal of chemical physics.

[87]  David Cohen-Steiner,et al.  Inferring Local Homology from Sampled Stratified Spaces , 2007, 48th Annual IEEE Symposium on Foundations of Computer Science (FOCS'07).

[88]  Vincent Danos,et al.  Rule-Based Modelling of Cellular Signalling , 2007, CONCUR.

[89]  Peter I. Booth,et al.  Monoidal closed, Cartesian closed and convenient categories of topological spaces , 1980 .

[90]  X. Trepat,et al.  Supracellular contraction at the rear of neural crest cell groups drives collective chemotaxis , 2018, Science.

[91]  Stéphane Lafon,et al.  Diffusion maps , 2006 .

[92]  M. Kac A stochastic model related to the telegrapher's equation , 1974 .