MobsPy: A programming language for biochemical reaction networks

Biochemical Reaction Networks (BCRNs) model species and their interactions via reactions. They have been extensively used in chemistry and extended to biological settings by generalizing the reactions' kinetics. However, detailed models of biochemical processes tend to result in complex BCRN models. We present the Meta-species Oriented Biosystem Syntax (MobsPy), a language designed to simplify the modeling process using the concept of meta-species. Meta-species are constructed using a bottom-up approach from base species, which represent elementary, simple characteristics. These characteristics are then combined to create meta-species with all their complex behavior. The combined species have characteristics that are the Cartesian product of the base species' characteristics and feature inheritance of reactions involving the base species. New reactions can involve all the states of a meta-species or only a subset that is selected via a query. In particular, reactions of meta-species can express a state change of one of the reactants. MobsPy is deployed as a Python package. We showcase its modeling capabilities by building concise models for biochemical systems from the literature.