Nonlinearity arising from noncooperative transcription factor binding enhances negative feedback and promotes genetic oscillations

Simplifications and approximations in a single-gene circuit modeling

The absence of detailed knowledge about regulatory interactions makes the use of phenomenological assumptions mandatory in cell biology modeling. Furthermore, the challenges associated with the analysis of these models compel the implementation of mathematical approximations. However, the constraints these methods introduce to biological interpretation are sometimes neglected. Consequently, understanding these restrictions is a very important task for systems biology modeling. In this article, we examine the impact of such simplifications, taking the case of a single-gene autoinhibitory circuit; however, our conclusions are not limited solely to this instance. We demonstrate that models grounded in the same biological assumptions but described at varying levels of detail can lead to different outcomes, that is, different and contradictory phenotypes or behaviors. Indeed, incorporating specific molecular processes like translation and elongation into the model can introduce instabilities and oscillations not seen when these processes are assumed to be instantaneous. Furthermore, incorporating a detailed description of promoter dynamics, usually described by a phenomenological regulatory function, can lead to instability, depending on the cooperative binding mechanism that is acting. Consequently, although the use of a regulating function facilitates model analysis, it may mask relevant aspects of the system's behavior. In particular, we observe that the two cooperative binding mechanisms, both compatible with the same sigmoidal function, can lead to different phenotypes, such as transcriptional oscillations with different oscillation frequencies.

Self-consistent theory of transcriptional control in complex regulatory architectures

Individual regulatory proteins are typically charged with the simultaneous regulation of a battery of different genes. As a result, when one of these proteins is limiting, competitive effects have a significant impact on the transcriptional response of the regulated genes. Here we present a general framework for the analysis of any generic regulatory architecture that accounts for the competitive effects of the regulatory environment by isolating these effects into an effective concentration parameter. These predictions are formulated using the grand-canonical ensemble of statistical mechanics and the fold-change in gene expression is predicted as a function of the number of transcription factors, the strength of interactions between the transcription factors and their DNA binding sites, and the effective concentration of the transcription factor. The effective concentration is set by the transcription factor interactions with competing binding sites within the cell and is determined self-consistently. Using this approach, we analyze regulatory architectures in the grand-canonical ensemble ranging from simple repression and simple activation to scenarios that include repression mediated by DNA looping of distal regulatory sites. It is demonstrated that all the canonical expressions previously derived in the case of an isolated, non-competing gene, can be generalised by a simple substitution to their grand canonical counterpart, which allows for simple intuitive incorporation of the influence of multiple competing transcription factor binding sites. As an example of the strength of this approach, we build on these results to present an analytical description of transcriptional regulation of the lac operon.

Feedback loops interlocked at competitive binding sites amplify and facilitate genetic oscillations

Positive and negative feedback loops are often present in regulatory networks for genetic oscillations. Relative time scales and integration of these feedback loops are key to robust oscillations in expression levels. Using examples from the circadian clock and synthetic genetic oscillators, we study positive and negative feedback loops interlocked at competitive binding sites. In the mammalian circadian clock, a key clock gene Bmal1 is regulated by the activator ROR and the repressor REV-ERB. Conversely, Bmal1 activates both of them, forming interlocked feedback loops. Previous experiments indicate that the activator and repressor compete for the same binding sites in the Bmal1 promoter. Transcription patterns predict that ROR peaks later than REV-ERB and, moreover, the peak phase difference between them is small. Using mathematical modeling we reveal an optimal ratio of dissociation constants of an activator and a repressor for the competitive binding sites to enhance the amplitude of Bmal1 oscillations. This optimal ratio arises only when the amplitude of the repressor is larger than that of the activator. Secondly, we reveal that the preference of binding sites for an activator and a repressor depends on their relative time scales. A previous study demonstrated that noncompetitive binding sites are preferable for synthetic genetic oscillators that comprise a fast activator and a slow repressor with a large time scale separation. Here we show that when their time scales are similar, competitive binding sites are more likely to generate oscillation than noncompetitive sites. In contrast, for a slow activator and a fast repressor with a small phase difference as in Bmal1 regulation, noncompetitive binding sites are advantageous for amplifying oscillations. Our results, therefore, predict that additional mechanisms are necessary to compensate the disadvantage of the Bmal1 promoter and further facilitate amplification under the regulation by ROR and REV-ERB.

Multiple binding sites for transcriptional repressors can produce regular bursting and enhance noise suppression

Cells may control fluctuations in protein levels by means of negative autoregulation, where transcription factors bind DNA sites to repress their own production. Theoretical studies have assumed a single binding site for the repressor, while in most species it is found that multiple binding sites are arranged in clusters. We study a stochastic description of negative autoregulation with multiple binding sites for the repressor. We find that increasing the number of binding sites induces regular bursting of gene products. By tuning the threshold for repression, we show that multiple binding sites can also suppress fluctuations. Our results highlight possible roles for the presence of multiple binding sites of negative autoregulators.

Promoters Architecture-Based Mechanism for Noise-Induced Oscillations in a Single-Gene Circuit

It is well known that single-gene circuits with negative feedback loop can lead to oscillatory gene expression when they operate with time delay. In order to generate these oscillations many processes can contribute to properly timing such delay. Here we show that the time delay coming from the transitions between internal states of the cis-regulatory system (CRS) can drive sustained oscillations in an auto-repressive single-gene circuit operating in a small volume like a cell. We found that the cooperative binding of repressor molecules is not mandatory for a oscillatory behavior if there are enough binding sites in the CRS. These oscillations depend on an adequate balance between the CRS kinetic, and the synthesis/degradation rates of repressor molecules. This finding suggest that the multi-site CRS architecture can play a key role for oscillatory behavior of gene expression. Finally, our results can also help to synthetic biologists on the design of the promoters architecture for new genetic oscillatory circuits.

Role of DNA binding sites and slow unbinding kinetics in titration-based oscillators

Genetic oscillators, such as circadian clocks, are constantly perturbed by molecular noise arising from the small number of molecules involved in gene regulation. One of the strongest sources of stochasticity is the binary noise that arises from the binding of a regulatory protein to a promoter in the chromosomal DNA. In this study, we focus on two minimal oscillators based on activator titration and repressor titration to understand the key parameters that are important for oscillations and for overcoming binary noise. We show that the rate of unbinding from the DNA, despite traditionally being considered a fast parameter, needs to be slow to broaden the space of oscillatory solutions. The addition of multiple, independent DNA binding sites further expands the oscillatory parameter space for the repressor-titration oscillator and lengthens the period of both oscillators. This effect is a combination of increased effective delay of the unbinding kinetics due to multiple binding sites and increased promoter ultrasensitivity that is specific for repression. We then use stochastic simulation to show that multiple binding sites increase the coherence of oscillations by mitigating the binary noise. Slow values of DNA unbinding rate are also effective in alleviating molecular noise due to the increased distance from the bifurcation point. Our work demonstrates how the number of DNA binding sites and slow unbinding kinetics, which are often omitted in biophysical models of gene circuits, can have a significant impact on the temporal and stochastic dynamics of genetic oscillators.