Linear control theory for gene network modeling

Formal reasoning about synthetic biology using higher-order-logic theorem proving

Synthetic biology is an interdisciplinary field that uses well-established engineering principles for performing the analysis of the biological systems, such as biological circuits, pathways, controllers and enzymes. Conventionally, the analysis of these biological systems is performed using paper-and-pencil proofs and computer simulation methods. However, these methods cannot ensure accurate results due to their inherent limitations. Higher-order-logic (HOL) theorem proving is proposed and used as a complementary approach for analysing linear biological systems, which is based on developing a mathematical model of the genetic circuits and the bio-controllers used in synthetic biology based on HOL and analysing it using deductive reasoning in an interactive theorem prover. The involvement of the logic, mathematics and the deductive reasoning in this method ensures the accuracy of the analysis. It is proposed to model the continuous dynamics of the genetic circuits and their associated controllers using differential equations and perform their transfer function-based analysis using the Laplace transform in a theorem prover. For illustration, the genetic circuits of activated and repressed expressions and autoactivation of protein, and phase lag and lead controllers, which are widely used in cancer-cell identifiers and multi-input receptors for precise disease detection, are formally analyzed.

Gene regulatory and signaling networks exhibit distinct topological distributions of motifs

The biological processes of cellular decision making and differentiation involve a plethora of signaling pathways and gene regulatory circuits. These networks in turn exhibit a multitude of motifs playing crucial parts in regulating network activity. Here we compare the topological placement of motifs in gene regulatory and signaling networks and observe that it suggests different evolutionary strategies in motif distribution for distinct cellular subnetworks.

Identification of control targets in Boolean molecular network models via computational algebra

BACKGROUND

Many problems in biomedicine and other areas of the life sciences can be characterized as control problems, with the goal of finding strategies to change a disease or otherwise undesirable state of a biological system into another, more desirable, state through an intervention, such as a drug or other therapeutic treatment. The identification of such strategies is typically based on a mathematical model of the process to be altered through targeted control inputs. This paper focuses on processes at the molecular level that determine the state of an individual cell, involving signaling or gene regulation. The mathematical model type considered is that of Boolean networks. The potential control targets can be represented by a set of nodes and edges that can be manipulated to produce a desired effect on the system.

RESULTS

This paper presents a method for the identification of potential intervention targets in Boolean molecular network models using algebraic techniques. The approach exploits an algebraic representation of Boolean networks to encode the control candidates in the network wiring diagram as the solutions of a system of polynomials equations, and then uses computational algebra techniques to find such controllers. The control methods in this paper are validated through the identification of combinatorial interventions in the signaling pathways of previously reported control targets in two well studied systems, a p53-mdm2 network and a blood T cell lymphocyte granular leukemia survival signaling network. Supplementary data is available online and our code in Macaulay2 and Matlab are available via http://www.ms.uky.edu/~dmu228/ControlAlg .

CONCLUSIONS

This paper presents a novel method for the identification of intervention targets in Boolean network models. The results in this paper show that the proposed methods are useful and efficient for moderately large networks.

Ultrasensitive Negative Feedback Control: A Natural Approach for the Design of Synthetic Controllers

Many of the most important potential applications of Synthetic Biology will require the ability to design and implement high performance feedback control systems that can accurately regulate the dynamics of multiple molecular species within the cell. Here, we argue that the use of design strategies based on combining ultrasensitive response dynamics with negative feedback represents a natural approach to this problem that fully exploits the strongly nonlinear nature of cellular information processing. We propose that such feedback mechanisms can explain the adaptive responses observed in one of the most widely studied biomolecular feedback systems—the yeast osmoregulatory response network. Based on our analysis of such system, we identify strong links with a well-known branch of mathematical systems theory from the field of Control Engineering, known as Sliding Mode Control. These insights allow us to develop design guidelines that can inform the construction of feedback controllers for synthetic biological systems.

Digital Signal Processing and Control for the Study of Gene Networks

Thanks to the digital revolution, digital signal processing and control has been widely used in many areas of science and engineering today. It provides practical and powerful tools to model, simulate, analyze, design, measure, and control complex and dynamic systems such as robots and aircrafts. Gene networks are also complex dynamic systems which can be studied via digital signal processing and control. Unlike conventional computational methods, this approach is capable of not only modeling but also controlling gene networks since the experimental environment is mostly digital today. The overall aim of this article is to introduce digital signal processing and control as a useful tool for the study of gene networks.

Designing Genetic Feedback Controllers

By incorporating feedback around systems we wish to manipulate, it is possible to improve their performance and robustness properties to meet pre-specified design objectives. For decades control engineers have been successfully implementing feedback controllers for complex mechanical and electrical systems such as aircraft and sports cars. Natural biological systems use feedback extensively for regulation and adaptation but apart from the most basic designs, there is no systematic framework for designing feedback controllers in Synthetic Biology. In this paper we describe how classical approaches from linear control theory can be used to close the loop. This includes the design of genetic circuits using feedback control and the presentation of a biological phase lag controller.

Modeling collective & intelligent decision making of multi-cellular populations

In the presence of unpredictable disturbances and uncertainties, cells intelligently achieve their goals by sharing information via cell-cell communication and making collective decisions, which are more reliable compared to individual decisions. Inspired by adaptive sensor network algorithms studied in communication engineering, we propose that a multi-cellular adaptive network can convert unreliable decisions by individual cells into a more reliable cell-population decision. It is demonstrated using the effector T helper (a type of immune cell) population, which plays a critical role in initiating immune reactions in response to invading foreign agents (e.g., viruses, bacteria, etc.). While each individual cell follows a simple adaptation rule, it is the combined coordination among multiple cells that leads to the manifestation of "self-organizing" decision making via cell-cell communication.

Post-translational regulation enables robust p53 regulation

Background

The tumor suppressor protein p53 plays important roles in DNA damage repair, cell cycle arrest and apoptosis. Due to its critical functions, the level of p53 is tightly regulated by a negative feedback mechanism to increase its tolerance towards fluctuations and disturbances. Interestingly, the p53 level is controlled by post-translational regulation rather than transcriptional regulation in this feedback mechanism.

Results

We analyzed the dynamics of this feedback to understand whether post-translational regulation provides any advantages over transcriptional regulation in regard to disturbance rejection. When a disturbance happens, even though negative feedback reduces the steady-state error, it can cause a system to become less stable and transiently overshoots, which may erroneously trigger downstream reactions. Therefore, the system needs to balance the trade-off between steady-state and transient errors. Feedback control and adaptive estimation theories revealed that post-translational regulation achieves a better trade-off than transcriptional regulation, contributing to a more steady level of p53 under the influence of noise and disturbances. Furthermore, post-translational regulation enables cells to respond more promptly to stress conditions with consistent amplitude. However, for better disturbance rejection, the p53- Mdm2 negative feedback has to pay a price of higher stochastic noise.

Conclusions

Our analyses suggest that the p53-Mdm2 feedback favors regulatory mechanisms that provide the optimal trade-offs for dynamic control.

Architecture-dependent robustness in a class of multiple positive feedback loops

Many types of multiple positive feedbacks with each having potentials to generate bistability exist extensively in natural, raising the question of why a particular architecture is present in a cell. In this study, the authors investigate multiple positive feedback loops across three classes: one-loop class, two-loop class and three-loop class, where each class is composed of double positive feedback loop (DPFL) or double negative feedback loop (DNFL) or both. Through large-scale sampling and robustness analysis, the authors find that for a given class, the homogeneous DPFL circuit (i.e. the coupled circuit that is composed of only DPFLs) is more robust than all the other circuits in generating bistable behaviour. In addition, stochastic simulation shows that the low stable state is more robust than the high stable state in homogeneous DPFL whereas the high-stable state is more robust than the low-stable state in homogeneous DNFL circuits. It was argued that this investigation provides insight into the relationship between robustness and network architecture.

Frequency domain analysis reveals external periodic fluctuations can generate sustained p53 oscillation

p53 is a well-known tumor suppressor protein that regulates many pathways, such as ones involved in cell cycle and apoptosis. The p53 levels are known to oscillate without damping after DNA damage, which has been a focus of many recent studies. A negative feedback loop involving p53 and MDM2 has been reported to be responsible for this oscillatory behavior, but questions remain as how the dynamics of this loop alter in order to initiate and maintain the sustained or undamped p53 oscillation. Our frequency domain analysis suggests that the sustained p53 oscillation is not completely dictated by the negative feedback loop; instead, it is likely to be also modulated by periodic DNA repair-related fluctuations that are triggered by DNA damage. According to our analysis, the p53-MDM2 feedback mechanism exhibits adaptability in different cellular contexts. It normally filters noise and fluctuations exerted on p53, but upon DNA damage, it stops performing the filtering function so that DNA repair-related oscillatory signals can modulate the p53 oscillation. Furthermore, it is shown that the p53-MDM2 feedback loop increases its damping ratio allowing p53 to oscillate at a frequency more synchronized with the other cellular efforts to repair the damaged DNA, while suppressing its inherent oscillation-generating capability. Our analysis suggests that the overexpression of MDM2, observed in many types of cancer, can disrupt the operation of this adaptive mechanism by making it less responsive to the modulating signals after DNA damage occurs.