Toggling Between Two Limit Cycles in a Molecular Ecosystem

In this article, we introduce an oscillatory molecular system with two independent limit cycles, dependent on initial conditions. The system is composed of two crossed Predator-Prey oscillators, in which the prey on one side is consumed by its predator to produce the predator of the other species. After providing a mathematical model for the system, we numerically analyze its behavior over a large range of enzymatic concentration conditions. Finally, we provide a surrogate model to evaluate the bistability region of the system to efficiently find optimal conditions.