How Retroactivity Affects the Behavior of Incoherent Feedforward Loops

Regulatory Effects of Cooperativity and Signal Profile on Adaptive Dynamics in Incoherent Feedforward Loop Networks

Cellular adaptation to external signals is essential for biological functions, and it is an important field of interest in systems biology. This study examines the impact of cooperativity on the adaptation response of the Incoherent Feedforward Loop (IFFL) network motif to various signal profiles. Through comprehensive simulations, we studied how the IFFL motif responds to constant and pulse-type signals under varying levels of cooperativity. The results of our study demonstrate that positive cooperativity generally enhances the system's ability to adapt to different signal profiles. Nevertheless, given specific signal profiles, higher levels of cooperativity may decrease the system's adaptability. On the other hand, the adaptive response breaks down for negative cooperativity. For constant signals, increased positive cooperativity leads to a response with higher amplitude, and it accelerates the response time but delays the return time required to settle back down to the pre-stimulus state. Upon signal cessation, high positive cooperativity not only slows the system's response and return times but, in some cases, can lead to a complete temporary halt in response. For the pulse-like signal, cooperativity increases the maximum amplitude of the oscillatory response. These insights highlight the delicate balance between cooperativity and signal profile in cellular adaptation mechanisms involving the IFFL network motif.

Design Principles for Perfect Adaptation in Biological Networks with Nonlinear Dynamics

Establishing a mapping between the emergent biological properties and the repository of network structures has been of great relevance in systems and synthetic biology. Adaptation is one such biological property of paramount importance that promotes regulation in the presence of environmental disturbances. This paper presents a nonlinear systems theory-driven framework to identify the design principles for perfect adaptation with respect to external disturbances of arbitrary magnitude. Based on the prior information about the network, we frame precise mathematical conditions for adaptation using nonlinear systems theory. We first deduce the mathematical conditions for perfect adaptation for constant input disturbances. Subsequently, we translate these conditions to specific necessary structural requirements for adaptation in networks of small size and then extend to argue that there exist only two classes of architectures for a network of any size that can provide local adaptation in the entire state space, namely, incoherent feed-forward (IFF) structure and negative feedback loop with buffer node (NFB). The additional positiveness constraints further narrow the admissible set of network structures. This also aids in establishing the global asymptotic stability for the steady state given a constant input disturbance. The proposed method does not assume any explicit knowledge of the underlying rate kinetics, barring some minimal assumptions. Finally, we also discuss the infeasibility of certain IFF networks in providing adaptation in the presence of downstream connections. Moreover, we propose a generic and novel algorithm based on non-linear systems theory to unravel the design principles for global adaptation. Detailed and extensive simulation studies corroborate the theoretical findings.

DOSE-L1000: unveiling the intricate landscape of compound-induced transcriptional changes

MOTIVATION

The LINCS L1000 project has collected gene expression profiles for thousands of compounds across a wide array of concentrations, cell lines, and time points. However, conventional analysis methods often fall short in capturing the rich information encapsulated within the L1000 transcriptional dose-response data.

RESULTS

We present DOSE-L1000, a database that unravels the potency and efficacy of compound-gene pairs and the intricate landscape of compound-induced transcriptional changes. Our study uses the fitting of over 140 million generalized additive models and robust linear models, spanning the complete spectrum of compounds and landmark genes within the LINCS L1000 database. This systematic approach provides quantitative insights into differential gene expression and the potency and efficacy of compound-gene pairs across diverse cellular contexts. Through examples, we showcase the application of DOSE-L1000 in tasks such as cell line and compound comparisons, along with clustering analyses and predictions of drug-target interactions. DOSE-L1000 fosters applications in drug discovery, accelerating the transition to omics-driven drug development.

AVAILABILITY AND IMPLEMENTATION

DOSE-L1000 is publicly available at https://doi.org/10.5281/zenodo.8286375.

Incoherent modulation of bi-stable dynamics orchestrates the Mushroom and Isola bifurcations

In biological networks, steady state dynamics of cell-fate regulatory genes often exhibit Mushroom and Isola kind of bifurcations. How these complex bifurcations emerge for these complex networks, and what are the minimal network structures that can generate these bifurcations, remain elusive. Herein, by employing Waddington's landscape theory and bifurcation analysis, we demonstrate that Mushroom and Isola bifurcations can be realized with four minimal network motifs that are constituted by combining a positive feedback motif with various incoherent feed-forward loops. Our study reveals that the intrinsic bi-stable dynamics originating from the positive feedback motif can be fine-tuned by altering the extent of the incoherence of these minimal networks to produce these complex bifurcations. These modeling insights will be useful in identifying the possible network motifs that may give rise to either Mushroom or Isola bifurcation in other biological systems.

Elucidating tumor-stromal metabolic crosstalk in colorectal cancer through integration of constraint-based models and LC-MS metabolomics

Colorectal cancer (CRC) is a major cause of morbidity and mortality in the United States. Tumor-stromal metabolic crosstalk in the tumor microenvironment promotes CRC development and progression, but exactly how stromal cells, in particular cancer-associated fibroblasts (CAFs), affect the metabolism of tumor cells remains unknown. Here we take a data-driven approach to investigate the metabolic interactions between CRC cells and CAFs, integrating constraint-based modeling and metabolomic profiling. Using metabolomics data, we perform unsteady-state parsimonious flux balance analysis to infer flux distributions for central carbon metabolism in CRC cells treated with or without CAF-conditioned media. We find that CAFs reprogram CRC metabolism through stimulation of glycolysis, the oxidative arm of the pentose phosphate pathway (PPP), and glutaminolysis, as well as inhibition of the tricarboxylic acid cycle. To identify potential therapeutic targets, we simulate enzyme knockouts and find that CAF-treated CRC cells are especially sensitive to inhibitions of hexokinase and glucose-6-phosphate, the rate limiting steps of glycolysis and oxidative PPP. Our work gives mechanistic insights into the metabolic interactions between CRC cells and CAFs and provides a framework for testing hypotheses towards CRC-targeted therapies.

Emergent Bistable Switches from the Incoherent Feed-Forward Signaling of a Positive Feedback Loop

Bistability is intrinsically connected to various decision making processes in living systems. The operating principles of a bistable switch, generated from a positive feedback loop, are well understood both in natural and synthetic settings. However, the fate of dynamic modularity of a positive feedback loop is unknown when it is connected to another dynamically modular signaling motif. In order to address this, here we investigate feed-forward signaling of a positive feedback loop to determine the fate of a bistable switch under such signaling. Using the potential energy based high-throughput bifurcation analysis method, we uncover that in addition to the conventional bistability the hybrid motifs generate various emergent bistable switches, namely mushroom and isola switches, which are not produced by the individual motifs. Using random parameter sampling, network perturbation, and phase plane analysis, we establish the design principles of such emergent behaviors. Incoherent feed-forward signaling of a positive feedback loop with distinct regulatory thresholds of the two arms of the feed-forward loop are the key requirements for such emergent behaviors. Our calculations show that the specific types of atypical bistable responses depend on the logic gate configuration of the signals. However, the emergent bistable behaviors of the hybrid networks do not depend on the nature of the positive feedback loop.

Design and Evaluation of Synthetic RNA-Based Incoherent Feed-Forward Loop Circuits

RNA-based regulators are promising tools for building synthetic biological systems that provide a powerful platform for achieving a complex regulation of transcription and translation. Recently, de novo-designed synthetic RNA regulators, such as the small transcriptional activating RNA (STAR), toehold switch (THS), and three-way junction (3WJ) repressor, have been utilized to construct RNA-based synthetic gene circuits in living cells. In this work, we utilized these regulators to construct type 1 incoherent feed-forward loop (IFFL) circuits in vivo and explored their dynamic behaviors. A combination of a STAR and 3WJ repressor was used to construct an RNA-only IFFL circuit. However, due to the fast kinetics of RNA–RNA interactions, there was no significant timescale difference between the direct activation and the indirect inhibition, that no pulse was observed in the experiments. These findings were confirmed with mechanistic modeling and simulation results for a wider range of conditions. To increase delay in the inhibition pathway, we introduced a protein synthesis process to the circuit and designed an RNA–protein hybrid IFFL circuit using THS and TetR protein. Simulation results indicated that pulse generation could be achieved with this RNA–protein hybrid model, and this was further verified with experimental realization in E. coli. Our findings demonstrate that while RNA-based regulators excel in speed as compared to protein-based regulators, the fast reaction kinetics of RNA-based regulators could also undermine the functionality of a circuit (e.g., lack of significant timescale difference). The agreement between experiments and simulations suggests that the mechanistic modeling can help debug issues and validate the hypothesis in designing a new circuit. Moreover, the applicability of the kinetic parameters extracted from the RNA-only circuit to the RNA–protein hybrid circuit also indicates the modularity of RNA-based regulators when used in a different context. We anticipate the findings of this work to guide the future design of gene circuits that rely heavily on the dynamics of RNA-based regulators, in terms of both modeling and experimental realization.