Abstract:For the vast majority of genes in sequenced genomes, there is limited understanding of how they are regulated. Without such knowledge, it is not possible to perform a quantitative theory-experiment dialogue on how such genes give rise to physiological and evolutionary adaptation. One category of high-throughput experiments used to understand the sequence-phenotype relationship of the transcriptome is massively parallel reporter assays (MPRAs). However, to improve the versatility and scalability of MPRA pipelines, we need a "theory of the experiment" to help us better understand the impact of various biological and experimental parameters on the interpretation of experimental data. These parameters include binding site copy number, where a large number of specific binding sites may titrate away transcription factors, as well as the presence of overlapping binding sites, which may affect analysis of the degree of mutual dependence between mutations in the regulatory region and expression levels. To that end, in this paper we create tens of thousands of synthetic single-cell gene expression outputs using both equilibrium and out-of-equilibrium models. These models make it possible to imitate the summary statistics (information footprints and expression shift matrices) used to characterize the output of MPRAs and from this summary statistic to infer the underlying regulatory architecture. Specifically, we use a more refined implementation of the so-called thermodynamic models in which the binding energies of each sequence variant are derived from energy matrices. Our simulations reveal important effects of the parameters on MPRA data and we demonstrate our ability to optimize MPRA experimental designs with the goal of generating thermodynamic models of the transcriptome with base-pair specificity. Further, this approach makes it possible to carefully examine the mapping between mutations in binding sites and their corresponding expression profiles, a tool useful not only for better designing MPRAs, but also for exploring regulatory evolution.

Thermodynamically-consistent, reduced models of gene regulatory networks

Synthesising Executable Gene Regulatory Networks from Single-cell Gene Expression Data

A bottom-up approach to gene regulation

Current approaches to gene regulatory network modelling

Theoretical and computational tools to model multistable gene regulatory networks

Deciphering regulatory architectures from synthetic single-cell expression patterns

Inferring gene regulatory networks from single-cell data: a mechanistic approach

Exactly solvable models of stochastic gene expression

A Spatial Artificial Chemistry Implementation of a Gene Regulatory Network Aimed at Generating Protein Concentration Dynamics

Modeling and Analyzing Biological Oscillations in Molecular Networks

A mathematical framework for describing and analysing gene regulatory networks

A coarse-grained bacterial cell model for resource-aware analysis and design of synthetic gene circuits

Predicting missing expression values in gene regulatory networks using a discrete logic modeling optimization guided by network stable states

Quantitative and logic modelling of molecular and gene networks

Modelling the burden caused by gene expression: an in silico investigation into the interactions between synthetic gene circuits and their chassis cell

A computational framework for optimal and Model Predictive Control of stochastic gene regulatory networks

Mathematical modeling of RNA-based architectures for closed loop control of gene expression.

Model-guided gene circuit design for engineering genetically stable cell populations in diverse applications

Deterministic and stochastic descriptions of gene expression dynamics

Reduction of dynamical biochemical reaction networks in computational biology

Optimising Boolean Synthetic Regulatory Networks to Control Cell States