Abstract
The parameter explosion problem is a crucial bottleneck in modelling gene regulatory networks (GRNs), limiting the size of models that can be optimised to experimental data. By discretising state, but not time, Boolean delay equations (BDEs) provide a significant reduction in parameter numbers, whilst still providing dynamical complexity comparable to more biochemically detailed models, such as those based on differential equations. Here, we explore several approaches to optimising BDEs to timeseries data, using a simple circadian clock model as a case study. We compare the effectiveness of two optimisers on our problem: a genetic algorithm (GA) and an elite accumulative sampling (EAS) algorithm that provides robustness to data discretisation. Our results show that both methods are able to distinguish effectively between alternative architectures, yielding excellent fits to data. We also perform a landscape analysis, providing insights into the properties that determine optimiser performance (e.g. number of local optima and basin sizes). Our results provide a promising platform for the analysis of more complex GRNs, and suggest the possibility of leveraging cost landscapes to devise more efficient optimisation schemes.