,
Salvador Buse
,
Erik Winfree
Creative Commons Attribution 4.0 International license
Many molecular systems are best understood in terms of prototypical species and reactions. The central dogma and related biochemistry are rife with examples: gene i is transcribed into RNA i, which is translated into protein i; kinase n phosphorylates substrate m; protein p dimerizes with protein q. Engineered nucleic acid systems also often have this form: oligonucleotide i hybridizes to complementary oligonucleotide j; signal strand n displaces the output of seesaw gate m; hairpin p triggers the opening of target q. When there are many variants of a small number of prototypes, it can be conceptually cleaner and computationally more efficient to represent the full system in terms of indexed species (e.g. for dimerization, M_p, D_pq) and indexed reactions (M_p + M_q → D_pq). Here, we formalize the Indexed Chemical Reaction Network (ICRN) model and describe a Python software package designed to simulate such systems in the well-mixed and reaction-diffusion settings, using a differentiable programming framework originally developed for large-scale neural network models, taking advantage of GPU acceleration when available. Notably, this framework makes it straightforward to train the models’ initial conditions and rate constants to optimize a target behavior, such as matching experimental data, performing a computation, or exhibiting spatial pattern formation. The natural map of indexed chemical reaction networks onto neural network formalisms provides a tangible yet general perspective for translating concepts and techniques from the theory and practice of neural computation into the design of biomolecular systems.
@InProceedings{lee_et_al:LIPIcs.DNA.31.4,
author = {Lee, Inhoo and Buse, Salvador and Winfree, Erik},
title = {{Differentiable Programming of Indexed Chemical Reaction Networks and Reaction-Diffusion Systems}},
booktitle = {31st International Conference on DNA Computing and Molecular Programming (DNA 31)},
pages = {4:1--4:23},
series = {Leibniz International Proceedings in Informatics (LIPIcs)},
ISBN = {978-3-95977-399-7},
ISSN = {1868-8969},
year = {2025},
volume = {347},
editor = {Schaeffer, Josie and Zhang, Fei},
publisher = {Schloss Dagstuhl -- Leibniz-Zentrum f{\"u}r Informatik},
address = {Dagstuhl, Germany},
URL = {https://drops.dagstuhl.de/entities/document/10.4230/LIPIcs.DNA.31.4},
URN = {urn:nbn:de:0030-drops-238534},
doi = {10.4230/LIPIcs.DNA.31.4},
annote = {Keywords: Differentiable Programming, Chemical Reaction Networks, Reaction-Diffusion Systems}
}