Molecular stripping, targets and decoys as modulators of oscillations in the NF-κB/IκBα/DNA genetic network

Phosphorylation of DNA-binding domains of CLOCK-BMAL1 complex for PER-dependent inhibition in circadian clock of mammalian cells

In mammals, CLOCK and BMAL1 proteins form a heterodimer that binds to E-box sequences and activates transcription of target genes, including Period (Per). Translated PER proteins then bind to the CLOCK-BMAL1 complex to inhibit its transcriptional activity. However, the molecular mechanism and the impact of this PER-dependent inhibition on the circadian clock oscillation remain elusive. We previously identified Ser38 and Ser42 in a DNA-binding domain of CLOCK as phosphorylation sites at the PER-dependent inhibition phase. In this study, knockout rescue experiments showed that nonphosphorylatable (Ala) mutations at these sites shortened circadian period, whereas their constitutive-phospho-mimetic (Asp) mutations completely abolished the circadian rhythms. Similarly, we found that nonphosphorylatable (Ala) and constitutive-phospho-mimetic (Glu) mutations at Ser78 in a DNA-binding domain of BMAL1 also shortened the circadian period and abolished the rhythms, respectively. The mathematical modeling predicted that these constitutive-phospho-mimetic mutations weaken the DNA binding of the CLOCK-BMAL1 complex and that the nonphosphorylatable mutations inhibit the PER-dependent displacement (reduction of DNA-binding ability) of the CLOCK-BMAL1 complex from DNA. Biochemical experiments supported the importance of these phosphorylation sites for displacement of the complex in the PER2-dependent inhibition. Our results provide direct evidence that phosphorylation of CLOCK-Ser38/Ser42 and BMAL1-Ser78 plays a crucial role in the PER-dependent inhibition and the determination of the circadian period.

Negative feedback systems for modelling NF-κB transcription factor oscillatory activity

Low-dimensional negative feedback systems (NFSs) were developed within a signal flow model to describe the oscillatory activities of NF-κB caused by interactions with its inhibitor IκBα. The NFSs were established as 3rd- and 4th-order linear systems containing unperturbed and perturbed negative feedback (NF) loops with constant or time-varying NF strengths and a feed-forward loop. NF-related analytical solutions to the NFSs representing the time courses of NF-κB and IκBα were determined and their exact mathematical relationship was found. The NFS's parameters were determined to fit the experimental time courses of NF-κB in TNF-α-stimulated embryonic fibroblasts, rela-/- embryonic fibroblasts reconstituted with RelA, C9L cells, GFP-p65 knock-in embryonic fibroblasts and embryogenic fibroblasts lacking Iκβ and IκBε, LPS-stimulated IC-21 macrophages treated or not with DCPA, and anti-IgM-stimulated DT40 B-lymphocytes. The unperturbed and perturbed NFSs describing the above biosystems generated isochronous and non-isochronous solutions, depending on a constant or time-varying NF strength, respectively. The oscillation period of the NF-coupled solutions, the phase difference between them and the time delays in the appearance of cytoplasmic IκBα after stimulation of NF-κB were determined. A significant divergence between the IκBα solutions to the NFSs and the IκBα experimental courses led to a rejection of the NF coupling between NF-κB and IκBα in the above biosystems. It was shown that neither the linearity nor the low dimensionality of the NFSs altered the NF relationship and the divergence between the IκBα solutions to the NFS and IκBα experimental time courses. Although the NF relationship between IκBα and NF-κB was not confirmed in all the experimental data analyzed, delayed negative feedback was found in some cases.

Extracellular Signal-Regulated Kinases: One Pathway, Multiple Fates

This comprehensive review delves into the multifaceted aspects of ERK signaling and the intricate mechanisms underlying distinct cellular fates. ERK1 and ERK2 (ERK) govern proliferation, transformation, epithelial-mesenchymal transition, differentiation, senescence, or cell death, contingent upon activation strength, duration, and context. The biochemical mechanisms underlying these outcomes are inadequately understood, shaped by signaling feedback and the spatial localization of ERK activation. Generally, ERK activation aligns with the Goldilocks principle in cell fate determination. Inadequate or excessive ERK activity hinders cell proliferation, while balanced activation promotes both cell proliferation and survival. Unraveling the intricacies of how the degree of ERK activation dictates cell fate requires deciphering mechanisms encompassing protein stability, transcription factors downstream of ERK, and the chromatin landscape.

The multifunctional role of intrinsic disorder in NF-κB signaling

The nuclear factor-κB (NF-κB) transcription activation system involves disordered regions of both the NF-κB dimers and their inhibitors, the IκBs. The system is well-studied both at the cellular and biophysical levels affording a unique opportunity to compare and contrast the conclusions from both types of experiments. Through a combination of both experiments and theory, we have discovered that the RelA/p50 heterodimer and its inhibitor IκBα operate under kinetic control. Intrinsically disordered parts of both proteins are directly involved in temporal control and their folding and unfolding determines the rates of various processes. In this review, we show how the dynamic state of the intrinsically disordered sequences define the rates of intracellular processes.

Design of Oscillatory Networks through Post-Translational Control of Network Components

Many essential functions in biological systems, including cell cycle progression and circadian rhythm regulation, are governed by the periodic behaviors of specific molecules. These periodic behaviors arise from the precise arrangement of components in biomolecular networks that generate oscillatory output signals. The dynamic properties of individual components of these networks, such as maturation delays and degradation rates, often play a key role in determining the network's oscillatory behavior. In this study, we explored the post-translational modulation of network components as a means to generate genetic circuits with oscillatory behaviors and perturb the oscillation features. Specifically, we used the NanoDeg platform-A bifunctional molecule consisting of a target-specific nanobody and a degron tag-to control the degradation rates of the circuit's components and predicted the effect of NanoDeg-mediated post-translational depletion of a key circuit component on the behavior of a series of proto-oscillating network topologies. We modeled the behavior of two main classes of oscillators, namely relaxation oscillator topologies (the activator-repressor and the Goodwin oscillator) and ring oscillator topologies (repressilators). We identified two main mechanisms by which non-oscillating networks could be induced to oscillate through post-translational modulation of network components: an increase in the separation of timescales of network components and mitigation of the leaky expression of network components. These results are in agreement with previous findings describing the effect of timescale separation and mitigation of leaky expression on oscillatory behaviors. This work thus validates the use of tools to control protein degradation rates as a strategy to modulate existing oscillatory signals and construct oscillatory networks. In addition, this study provides the design rules to implement such an approach based on the control of protein degradation rates using the NanoDeg platform, which does not require genetic manipulation of the network components and can be adapted to virtually any cellular protein. This work also establishes a framework to explore the use of tools for post-translational perturbations of biomolecular networks and generates desired behaviors of the network output.

Combined multiple transcriptional repression mechanisms generate ultrasensitivity and oscillations

Transcriptional repression can occur via various mechanisms, such as blocking, sequestration and displacement. For instance, the repressors can hold the activators to prevent binding with DNA or can bind to the DNA-bound activators to block their transcriptional activity. Although the transcription can be completely suppressed with a single mechanism, multiple repression mechanisms are used together to inhibit transcriptional activators in many systems, such as circadian clocks and NF-κB oscillators. This raises the question of what advantages arise if seemingly redundant repression mechanisms are combined. Here, by deriving equations describing the multiple repression mechanisms, we find that their combination can synergistically generate a sharply ultrasensitive transcription response and thus strong oscillations. This rationalizes why the multiple repression mechanisms are used together in various biological oscillators. The critical role of such combined transcriptional repression for strong oscillations is further supported by our analysis of formerly identified mutations disrupting the transcriptional repression of the mammalian circadian clock. The hitherto unrecognized source of the ultrasensitivity, the combined transcriptional repressions, can lead to robust synthetic oscillators with a previously unachievable simple design.

Diverse role of decoys on emergence and precision of oscillations in a biomolecular clock

Biomolecular clocks are key drivers of oscillatory dynamics in diverse biological processes including cell-cycle regulation, circadian rhythms, and pattern formation during development. A minimal clock implementation is based on the classical Goodwin oscillator, in which a repressor protein inhibits its own synthesis via time-delayed negative feedback. Clock motifs, however, do not exist in isolation; its components are open to interacting with the complex environment inside cells. For example, there are ubiquitous high-affinity binding sites along the genome, known as decoys, where transcription factors such as repressor proteins can potentially interact. This binding affects the availability of transcription factors and has often been ignored in theoretical studies. How does such genomic decoy binding impact the clock's robustness and precision? To address this question, we systematically analyze deterministic and stochastic models of the Goodwin oscillator in the presence of reversible binding of the repressor to a finite number of decoy sites. Our analysis reveals that the relative stability of decoy-bound repressors compared to the free repressor plays distinct roles on the emergence and precision of oscillations. Interestingly, active degradation of the bound repressor can induce sustained oscillations that are otherwise absent without decoys. In contrast, decoy abundances can kill oscillation dynamics if the bound repressor is protected from degradation. Taking into account low copy-number fluctuations in clock components, we show that the degradation of the bound repressors enhances precision by attenuating noise in both the amplitude and period of oscillations. Overall, these results highlight the versatile role of otherwise hidden decoys in shaping the stochastic dynamics of biological clocks and emphasize the importance of synthetic decoys in designing robust clocks.

Therapeutic potential of beta-caryophyllene against aflatoxin B1-Induced liver toxicity: biochemical and molecular insights in rats

Aflatoxin B1 (AFB1) is a mycotoxin highly toxic and carcinogenic to humans due to its potential to induce oxidative stress. The Beta-caryophyllene (BCP) have been highlighted for its broad spectrum of pharmacological effects. The present study aimed to investigate the beneficial effects of BCP against the susceptibility of hepatic and renal tissues to AFB1 toxicity, in biochemical parameters to assess organ function, tissue oxidation, and the immunocontent of oxidative and inflammatory proteins. Male Wistar rats was exposed to AFB1 (250 μg/kg, i.g.) and/or BCP (100 mg/kg, i.p.) for 14 successive days. It was found that exposure to AFB1 did not change the measured renal toxicity parameters. Also, AFB1 increased liver injury biomarkers (gamma glutamyl transferase and alkaline phosphatase) and reduced levels of non-enzymatic antioxidant defenses (ascorbic acid and non-protein thiol), however did not cause changes in the lipid peroxidation levels. Moreover, AFB1 interfered in oxidative pathway regulated by Kelch-like ECH-associated protein (Keap1)/nuclear factor (erythroid-derived 2)-like 2 (Nrf2), overacting Glutathione-S-Transferase (GST) activity. Lastly, a main effect of AFB1 on the total interleukin 1 beta (IL-1β) was observed. Remarkably, the associated treatment of AFB1 + BCP improved altered liver parameters. In addition, BCP and AFB1 + BCP groups showed an increase in the levels of inhibitor of nuclear factor kappa-B kinase subunit beta (IKKβ). Thus, these results indicated that BCP has potential protective effect against AFB1 induced hepatotoxicity.

How Retroactivity Affects the Behavior of Incoherent Feedforward Loops

An incoherent feedforward loop (IFFL) is a network motif known for its ability to accelerate responses and generate pulses. It remains an open question to understand the behavior of IFFLs in contexts with high levels of retroactivity, where an upstream transcription factor binds to numerous downstream binding sites. Here we study the behavior of IFFLs by simulating and comparing ODE models with different levels of retroactivity. We find that increasing retroactivity in an IFFL can increase, decrease, or keep the network's response time and pulse amplitude constant. This suggests that increasing retroactivity, traditionally considered an impediment to designing robust synthetic systems, could be exploited to improve the performance of IFFLs. In contrast, we find that increasing retroactivity in a negative autoregulated circuit can only slow the response. The ability of an IFFL to flexibly handle retroactivity may have contributed to its significant abundance in both bacterial and eukaryotic regulatory networks.

Polyelectrolyte interactions enable rapid association and dissociation in high-affinity disordered protein complexes

Highly charged intrinsically disordered proteins can form complexes with very high affinity in which both binding partners fully retain their disorder and dynamics, exemplified by the positively charged linker histone H1.0 and its chaperone, the negatively charged prothymosin α. Their interaction exhibits another surprising feature: The association/dissociation kinetics switch from slow two-state-like exchange at low protein concentrations to fast exchange at higher, physiologically relevant concentrations. Here we show that this change in mechanism can be explained by the formation of transient ternary complexes favored at high protein concentrations that accelerate the exchange between bound and unbound populations by orders of magnitude. Molecular simulations show how the extreme disorder in such polyelectrolyte complexes facilitates (i) diffusion-limited binding, (ii) transient ternary complex formation, and (iii) fast exchange of monomers by competitive substitution, which together enable rapid kinetics. Biological polyelectrolytes thus have the potential to keep regulatory networks highly responsive even for interactions with extremely high affinities.

Enhancement of gene expression noise from transcription factor binding to genomic decoy sites

The genome contains several high-affinity non-functional binding sites for transcription factors (TFs) creating a hidden and unexplored layer of gene regulation. We investigate the role of such “decoy sites” in controlling noise (random fluctuations) in the level of a TF that is synthesized in stochastic bursts. Prior studies have assumed that decoy-bound TFs are protected from degradation, and in this case decoys function to buffer noise. Relaxing this assumption to consider arbitrary degradation rates for both bound/unbound TF states, we find rich noise behaviors. For low-affinity decoys, noise in the level of unbound TF always monotonically decreases to the Poisson limit with increasing decoy numbers. In contrast, for high-affinity decoys, noise levels first increase with increasing decoy numbers, before decreasing back to the Poisson limit. Interestingly, while protection of bound TFs from degradation slows the time-scale of fluctuations in the unbound TF levels, the decay of bound TFs leads to faster fluctuations and smaller noise propagation to downstream target proteins. In summary, our analysis reveals stochastic dynamics emerging from nonspecific binding of TFs and highlights the dual role of decoys as attenuators or amplifiers of gene expression noise depending on their binding affinity and stability of the bound TF.

Force-Dependent Facilitated Dissociation Can Generate Protein-DNA Catch Bonds

Cellular structures are continually subjected to forces, which may serve as mechanical signals for cells through their effects on biomolecule interaction kinetics. Typically, molecular complexes interact via "slip bonds," so applied forces accelerate off rates by reducing transition energy barriers. However, biomolecules with multiple dissociation pathways may have considerably more complicated force dependencies. This is the case for DNA-binding proteins that undergo "facilitated dissociation," in which competitor biomolecules from solution enhance molecular dissociation in a concentration-dependent manner. Using simulations and theory, we develop a generic model that shows that proteins undergoing facilitated dissociation can form an alternative type of molecular bond, known as a "catch bond," for which applied forces suppress protein dissociation. This occurs because the binding by protein competitors responsible for the facilitated dissociation pathway can be inhibited by applied forces. Within the model, we explore how the force dependence of dissociation is regulated by intrinsic factors, including molecular sensitivity to force and binding geometry and the extrinsic factor of competitor protein concentration. We find that catch bonds generically emerge when the force dependence of the facilitated unbinding pathway is stronger than that of the spontaneous unbinding pathway. The sharpness of the transition between slip- and catch-bond kinetics depends on the degree to which the protein bends its DNA substrate. This force-dependent kinetics is broadly regulated by the concentration of competitor biomolecules in solution. Thus, the observed catch bond is mechanistically distinct from other known physiological catch bonds because it requires an extrinsic factor-competitor proteins-rather than a specific intrinsic molecular structure. We hypothesize that this mechanism for regulating force-dependent protein dissociation may be used by cells to modulate protein exchange, regulate transcription, and facilitate diffusive search processes.

Efficient analysis of stochastic gene dynamics in the non-adiabatic regime using piecewise deterministic Markov processes

Single-cell experiments show that gene expression is stochastic and bursty, a feature that can emerge from slow switching between promoter states with different activities. In addition to slow chromatin and/or DNA looping dynamics, one source of long-lived promoter states is the slow binding and unbinding kinetics of transcription factors to promoters, i.e. the non-adiabatic binding regime. Here, we introduce a simple analytical framework, known as a piecewise deterministic Markov process (PDMP), that accurately describes the stochastic dynamics of gene expression in the non-adiabatic regime. We illustrate the utility of the PDMP on a non-trivial dynamical system by analysing the properties of a titration-based oscillator in the non-adiabatic limit. We first show how to transform the underlying chemical master equation into a PDMP where the slow transitions between promoter states are stochastic, but whose rates depend upon the faster deterministic dynamics of the transcription factors regulated by these promoters. We show that the PDMP accurately describes the observed periods of stochastic cycles in activator and repressor-based titration oscillators. We then generalize our PDMP analysis to more complicated versions of titration-based oscillators to explain how multiple binding sites lengthen the period and improve coherence. Last, we show how noise-induced oscillation previously observed in a titration-based oscillator arises from non-adiabatic and discrete binding events at the promoter site.

Modeling the therapeutic efficacy of NFκB synthetic decoy oligodeoxynucleotides (ODNs)

BACKGROUND

Transfection of NF κB synthetic decoy Oligodeoxynucleotides (ODNs) has been proposed as a promising therapeutic strategy for a variety of diseases arising from constitutive activation of the eukaryotic transcription factor NF κB. The decoy approach faces some limitations under physiological conditions notably nuclease-induced degradation.

RESULTS

In this work, we show how a systems pharmacology model of NF κB regulatory networks displaying oscillatory temporal dynamics, can be used to predict quantitatively the dependence of therapeutic efficacy of NF κB synthetic decoy ODNs on dose, unbinding kinetic rates and nuclease-induced degradation rates. Both deterministic mass action simulations and stochastic simulations of the systems biology model show that the therapeutic efficacy of synthetic decoy ODNs is inversely correlated with unbinding kinetic rates, nuclease-induced degradation rates and molecular stripping rates, but is positively correlated with dose. We show that the temporal coherence of the stochastic dynamics of NF κB regulatory networks is most sensitive to adding NF κB synthetic decoy ODNs having unbinding time-scales that are in-resonance with the time-scale of the limit cycle of the network.

CONCLUSIONS

The pharmacokinetics/pharmacodynamics (PK/PD) predicted by the systems-level model should provide quantitative guidance for in-depth translational research of optimizing the thermodynamics/kinetic properties of synthetic decoy ODNs.

Potential role of DNA methylation as a facilitator of target search processes for transcription factors through interplay with methyl-CpG-binding proteins

Eukaryotic genomes contain numerous non-functional high-affinity sequences for transcription factors. These sequences potentially serve as natural decoys that sequester transcription factors. We have previously shown that the presence of sequences similar to the target sequence could substantially impede association of the transcription factor Egr-1 with its targets. In this study, using a stopped-flow fluorescence method, we examined the kinetic impact of DNA methylation of decoys on the search process of the Egr-1 zinc-finger protein. We analyzed its association with an unmethylated target site on fluorescence-labeled DNA in the presence of competitor DNA duplexes, including Egr-1 decoys. DNA methylation of decoys alone did not affect target search kinetics. In the presence of the MeCP2 methyl-CpG-binding domain (MBD), however, DNA methylation of decoys substantially (∼10-30-fold) accelerated the target search process of the Egr-1 zinc-finger protein. This acceleration did not occur when the target was also methylated. These results suggest that when decoys are methylated, MBD proteins can block them and thereby allow Egr-1 to avoid sequestration in non-functional locations. This effect may occur in vivo for DNA methylation outside CpG islands (CGIs) and could facilitate localization of some transcription factors within regulatory CGIs, where DNA methylation is rare.

Interrogating the topological robustness of gene regulatory circuits by randomization

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Cells are able to robustly carry out their essential biological functions, possibly because of multiple layers of tight regulation via complex, yet well-designed, gene regulatory networks involving a substantial number of genes. State-of-the-art genomics technology has enabled the mapping of these large gene networks, yet it remains a tremendous challenge to elucidate their design principles and the regulatory mechanisms underlying their biological functions such as signal processing and decision-making. One of the key barriers is the absence of accurate kinetics for the regulatory interactions, especially from in vivo experiments. To this end, we have developed a new computational modeling method, Random Circuit Perturbation (RACIPE), to explore the dynamic behaviors of gene regulatory circuits without the requirement of detailed kinetic parameters. RACIPE takes a network topology as the input, and generates an unbiased ensemble of models with varying kinetic parameters. Each model is subjected to simulation, followed by statistical analysis for the ensemble. We tested RACIPE on several gene circuits, and found that the predicted gene expression patterns from all of the models converge to experimentally observed gene state clusters. We expect RACIPE to be a powerful method to identify the role of network topology in determining network operating principles.

Evaluation of 8-hydroxy-2-deoxyguanosine and NFkB activation, oxidative stress response, acetylcholinesterase activity, and histopathological changes in rainbow trout brain exposed to linuron

Linuron is a widely used herbicide to control grasses and annual broad leaf weeds. It is known that linuron has toxic effects on different organisms. However, the toxic effects of linuron on aquatic organisms, especially fish, is completely unknown. Thus, we aimed to investigate changes in 8-hydroxy-2-deoxyguanosine (8-OHdG) and nuclear factor kappa B (NFkB) activity, histopathological changes, antioxidant responses and acetylcholinesterase (AChE) activity in rainbow trout brain after exposure to linuron. Fish were exposed to 30μg/L, 120μg/L and 240μg/L concentrations of linuron for twenty-one days. Brain tissues were taken from fish for 8-OHdG and NFkB activity, histopathological examination and determination of superoxide dismutase (SOD), catalase (CAT) enzyme activity, lipid peroxidation (LPO), and reduced glutathione (GSH) levels. Our data indicated that high linuron concentrations caused a decrease in GSH levels, SOD and CAT activities in brain tissues (p<0.05). LPO levels were significantly increased by 240μg/L linuron. All concentrations caused a significant inhibition in brain AChE enzyme activity (p<0.05). Immunopositivity was detected for 8-OHdG and NFkB, and linuron exposure caused histopathological damage to the brain tissues. The results of this study can provide useful information for understanding of linuron-induced toxicity.

A comparative review of computational methods for pathway perturbation analysis: dynamical and topological perspectives

Stem cells offer great promise within the field of regenerative medicine but despite encouraging results, the large scale use of stem cells for therapeutic applications still faces challenges when it comes to controlling signaling pathway responses with respect to environmental perturbations.