Team:Duke/Modeling

Introduction

Synthetic biologists, in their goal of designing biological systems with novel functions, are particularly interested in design techniques for synthetic genetic networks, given the many emergent properties observed in natural gene circuits. The simplest of these gene circuits is the bistable toggle switch, in which mutually repressing genes create a system that can be tuned to make either gene dominant at will. Previous studies have shown that an essential criterion for the bistability of these systems is the cooperative action of the repressors. Several toggle switch systems, including the first switch pioneered by Collins, have been proposed to introduce cooperativity in the form of the cooperative binding of multiple repressor molecules. In this study, we consider alternative methods for generating cooperativity though molecular titration.

Analogous to chemical titration, molecular titration generates a nonlinear “all-or-nothing” response with a “sink” of buffer molecules that sequester and inactivate repressor molecules at low concentration of repressor. A naturally occurring example of this is protein sequestration, and a previous study has constructed a synthetic circuit with proteins using this principle and demonstrated an ultrasensitive response in yeast. We propose to look at novel forms of molecular titration that involve the widely used dCas9 system in an attempt to construct new building blocks for synthetic circuits. In particular, instead of having a protein and its sequestering protein as a buffer, we wish to use small RNA molecules (anti gRNAs) and decoy binding sites on a plasmid that are complementary to the gRNA component of the dcas9 system to suppress the activity of the cas9 complex and act as a buffer against repression. In this study, we will model the properties of this system by considering how its equilibrium conditions vary with respect to total gRNA concentration and we hope to demonstrate this system’s ability to exhibit an ultrasensitive response under biologically reasonable conditions and to characterize its behavior in general over a wide range of parameters.

Methods

The approach we took in modeling our systems is based on Buchler’s approach in modeling molecular titration in a system consisting of a generic repressor and its buffering molecule (1). Their model considers the equilibrium concentration of the active repressor. Taking the total concentration of the other molecules in the system as constants, it is possible to calculate this concentration as a function of the total concentration of repressor. This function is obtained as a solution to a system of equations involving the equilibrium expression for the sequestering reaction in addition to the mass balance equations for each component of the system. The resulting system of equations can be solved as a quadratic polynomial in terms of the concentration of active repressor, giving an explicit expression for this equilibrium concentration in terms of the dissociation constants and the total concentrations of each component in the system.

Figure 1: Biochemical reactions involved in the (a) gRNA:anti-gRNA system, and the (b) decoy binding sites system

The system that we want to model using Cas9 is more complex in that neither Cas9 nor the gRNA will actively repress transcription on its own. Because we’re interested in the concentration of the active Cas9:gRNA complex, the equilibrium involving the formation of the active complex needs to be considered, in addition to the equilibrium of the sequestering reaction. Please see the appendix for details of how the resulting polynomial was solved for. This additional equation in the system of equations has the effect of increasing the degree of this polynomial so that it is no longer practical to explicitly solve for the concentration of Cas9:gRNA. To obtain solutions to these polynomials, we used the MATLAB vpasolve function to calculate the numerical solutions for the value of Cas9:gRNA concentration over a range of values for total gRNA concentration. Manual inspection to eliminate the other physically impossible solutions to the polynomial was necessary to obtain the physically relevant graph of active complex versus total gRNA concentration.

A challenge in modeling these systems was finding the equilibrium constants that describe these reactions, as the biochemistry of the Cas9/gRNA system is still a fairly unexplored field of research. Where available, these values were taken directly from published sources; but for some of the parameters, it was necessary to make approximations based on the closest available information.

As can be seen from Figure 1, there are 3 distinct dissociation constants that are relevant to our systems.

Table 1: Values of dissociation constants used in modeling

The first, the dissociation constant of the gRNA:Cas9 complex, was estimated using a dissociation constant obtained through a gel-mobility shift assay of a RNA binding protein that targets small 28bp RNA, which is roughly the same length as the gRNAs we are interested in (2). For the dissociation constant of the gRNA:Cas9:decoyDNA complex, we used the dissociation constant for the complex as reported by Jennifer Doudna (3). Lastly, for the dissociation constants for the formation of the gRNA:anti-gRNA duplex, they were computed from the Gibbs free energy change for the binding of gRNA:anti-gRNA duplexes of similar length. Unfortunately, this approach yielded a very large range in values for the dissociation constant, so we decided to run the model for both extremes reported (4),(5).