Constraint and Contingency in Multifunctional Gene Regulatory Circuits
Gene regulatory circuits are at the heart of many fundamental biological processes , ranging from developmental patterning in multicellular organisms 1 to chemotaxis in bacteria 2 ., Regulatory circuits are usually multifunctional ., This means that they can form different metastable gene expression states under different physiological conditions , in different tissues , or in different stages of embryonic development ., The segment polarity network of Drosophila melanogaster offers an example , where the same regulatory circuit affects several developmental processes , including embryonic segmentation and the development of the flys wing 3 ., Similarly , in the vertebrate neural tube , a single circuit is responsible for interpreting a morphogen gradient to produce three spatially distinct ventral progenitor domains 4 ., Other notable examples include the bistable competence control circuit of Bacillus subtilis 5 and the lysis-lysogeny switch of bacteriophage lambda 6 ., Multifunctional regulatory circuits are also relevant to synthetic biology , where artificial oscillators 7 , toggle switches 8 , and logic gates 9 are engineered to control biological processes ., The functions of gene regulatory circuits are embodied in their gene expression patterns ., An important property of natural circuits , and a design goal of synthetic circuits , is that these patterns should be robust to perturbations ., Such perturbations include nongenetic perturbations , such as stochastic fluctuations in protein concentrations and environmental change ., Much attention has focused on understanding 1 , 2 , 4 , 10 , 11 and engineering 12–14 circuits that are robust to nongenetic perturbations ., Equally important is the robustness of circuit functions to genetic perturbations , such as those caused by point mutation or recombination ., Multiple studies have asked what renders biological circuitry robust to such genetic changes 15–20 ., With few exceptions 21 , 22 , these studies have focused on circuits with one function , embodied in their gene expression pattern ., Such monofunctional circuits tend to have several properties ., First , many circuits exist that have the same gene expression pattern 17–19 , 23–28 ., Second , these circuits can vary greatly in their robustness 16 , 18 , 29 ., And third , they can often be reached from one another via a series of function-preserving mutational events 18 , 19 , 30 ., Taken together , these observations suggest that the robustness of the many circuits with a given regulatory function can be tuned via incremental mutational change ., Most circuits have multiple functions , but how these observations translate to such multifunctional circuits is largely unknown ., In a given space of possible circuits , how many circuits exist that have a given number of k specific functions ( expression patterns ) ?, What is the relationship between this number of functions and the robustness of each function ?, Do circuits with any combination of functions exist , or are some combinations “prohibited ? ”, Pertinent earlier work showed that there are indeed fewer multifunctional circuits than monofunctional circuits 21 , but this investigation had two main limitations ., First , it considered circuits so large that the space of circuits and their functions could not be exhaustively explored , and restricted itself to mostly bifunctional circuits ., Second , it included only topological circuit variants ( i . e . , who interacts with whom ) , and ignored variations in the signal-integration logic of cis-regulatory regions ., These regions encode regulatory programs , which specify the input-output mapping of regulatory signals ( input ) to gene expression pattern ( output ) 31–33 ., Variations in cis-regulatory regions 34 , such as mutations that change the spacing between transcription factor binding sites 35 , are known to impact circuit function 36 , 37 , and their inclusion in a computational model of regulatory circuits is thus important ., Here , we overcome these limitations by focusing on regulatory circuits that are sufficiently small that an entire space of circuits can be exhaustively explored ., Specifically , we focus on circuits that comprise only three genes and all possible regulatory interactions between them ., Small circuits like this play an important role in some biological processes ., Examples include the kaiABC gene cluster in Cyanobacteria , which is responsible for circadian oscillations 38 , the gap gene system in Dropsophila , which is responsible for the interpretation of morphogen gradients during embryogenesis 19 , and the krox-otx-gatae feedback loop in starfish , which is necessary for endoderm specification 39 ., Additionally , theoretical studies of small regulatory circuits have provided several general insights into the features of circuit design and function ., Examples include biochemical adaptation in feedback loops 40 and response delays in feed-forward loops 41 , among others 16 , 19 , 23 , 42–45 ., Lastly , there is a substantial body of evidence suggesting that small regulatory circuits form the building blocks of larger regulatory networks 34 , 46–48 , further warranting their study ., For two reasons , we chose Boolean logic circuits 49 as our modeling framework ., First , they allow us not only to vary circuit topology 45 , but also a circuits all-important signal-integration logic 44 ., Second , Boolean circuits have been successful in explaining properties of biological circuits ., For example , they have been used to explain the dynamics of gene expression in the segment polarity genes of Drosophila melanogaster 50 , the development of primordial floral organ cells of Arabidopsis thaliana 51 , gene expression cascades after gene knockout in Saccharomyces cerevisiae 52 , and the temporal and spatial expression dynamics of the genes responsible for endomesoderm specification in the sea urchin embryo 53 ., We consider a specific gene expression pattern as the function of a circuit like this , because it is this pattern that ultimately drives embryonic pattern formation and physiological processes ., Multifunctional circuits are circuits with multiple gene expression patterns , and here we study the constraints that multifunctionality imposes on the robustness and other properties of regulatory circuits ., The questions we ask include the following:, ( i ) How many circuits have a given number k of functions ?, ( ii ) What is the relationship between multifunctionality and robustness to genetic perturbation ?, ( iii ) Are some multifunctional circuits more robust than others ?, ( iv ) Is it possible to change one multifunctional circuit into another through a series of small genetic changes that do not jeopardize circuit function ?, We consider circuits of genes ( Fig . 1A ) ., We choose a compact representation of a circuits genotype G that allows us to represent both a circuits signal-integration logic and its architecture by a single binary vector of length ( Fig . 1B ) ., Changes to this vector can be caused by mutations in the cis-regulatory regions of DNA ., Such mutations may alter the binding affinity of a transcription factor to its binding site , thereby creating or removing a regulatory interaction 34 ., Alternatively , they may affect the distance of a transcription factor binding site from the transcription start site , changing its rotational position on the DNA helix ., In turn , this may alter the regulatory effect of the transcription factor 54 , and change the downstream genes signal-integration logic ., Lastly , such mutations may change the distance between adjacent transcription factor binding sites , enabling or disabling a functional interaction between proximally bound transcription factors 35 ., We note that mutations in G could also be conceptualized as changes in the DNA binding domain of a transcription factor ., However , evolutionary evidence from microbes suggest that alterations in the structure and logic of regulatory circuits occurs preferentially via changes in cis-regulatory regions , rather than via changes in the transcription factors that bind these regions 55 ., The dynamics of the expression states of a circuits N genes begin with a prespecified initial state , which represents regulatory influences outside or upstream of the circuit , such as transcription factors that are not part of the circuit but can influence its expression state ., The initial state reflects the fact that small circuits are typically embedded in larger regulatory networks 34 , 46–48 , which provide the circuit with different regulatory inputs under different environmental or tissue-specific conditions ., Through the regulatory interactions specified in the circuits genotype , the circuits gene expression state changes from this initial state , until it may reach a stable ( i . e . , fixed-point ) equilibrium state ., We consider a circuits function to be a mapping from an initial expression state to an equilibrium expression state ( Fig . 1C ) ., In the main text , we consider only circuit functions that involve fixed point equilibria , but we consider periodic equilibrium states in the Supporting Online Material ., A circuit could in principle have as many as functions , as long as the initial expression states are all different from one another , and the equilibrium expression states are all different from one another ( Material and Methods ) ., The circuits we study may map multiple initial states to the same equilibrium state , but our definition of function ignores all but one of these initial states ., While a definition of function that includes many-to-one mappings between initial and equilibrium states can be biologically sensible , our intent is to investigate specific pairs of inputs ( i . e . , ) and outputs ( i . e . , ) , as is typical for circuits in development and physiology 56–58 ., We emphasize that a circuit can express its k functions individually , or in various combinations , such that the same circuit could be said to have between one and k functions ., For brevity , we refer to a specific set of k functions as a multifunction or a k-function and to circuits that have at least one function as viable ., The space of circuits we explore here contains possible genotypes ., We exhaustively determine the equilibrium expression states of each genotype for all initial states , thereby providing a complete genotype-to-phenotype ( function ) map ., We use this map to partition the space of genotypes into genotype networks 17–19 , 21 ., A genotype network consists of a single connected set of genotypes ( circuits ) that have identical functions , and where two circuits are connected neighbors if their corresponding genotypes differ by a single element ( Fig . 1D ) ., Note that such single mutations may correspond to larger mutational changes in the cis-regulatory regions of DNA ., For example , mutations that change the distance between binding sites , or between a binding site and a transcription start site , may involve the addition or deletion of large segments of DNA 26 , 59–62 ., We first asked how the number of genotypes that have k functions depends on k ., Fig . 2 shows that this number decreases exponentially , implying that multifunctionality constrains the number of viable genotypes severely ., For instance , increasing k from 1 to 2 decreases the number of viable genotypes by 34%; further increasing k from 2 to 3 leads to an additional 39% decrease ., However , there is always at least one genotype with a given number k of functions , for any ., In other words , even in these small circuits , multiple genotypes exist that have many functions ., Thus far , we have determined the number of genotypes with a given number k of functions , but we did not distinguish between the actual functions that these genotypes can have ., For example , there are 64 variants of function , since there are potential initial states and potential equilibrium states ( ) ., Analogously , simple combinatorics ( Text S1 ) shows that there are 1204 variants of functions , and the number of variants increases dramatically with greater k , up to a maximum of variants of functions ., This is possible because individual functions can occur in different possible combinations in multifunctional circuits ( Material and Methods ) ., The solid line in the inset of Fig . 2 indicates how this number of possible different functions scales with k ., We next asked whether there exist circuits ( genotypes ) for each of these possible combinations of functions , or whether some multifunctions are prohibited ., The open circles in the inset of Fig . 2 show the answer: These circles lie exactly on the solid line that indicates the number of possible combinations of functions for each value of k ( Text S1 ) ., This means that no multifunction is prohibited ., In other words , even though multifunctionality constrains the number of viable genotypes , there is always at least one genotype with k functions , and in any possible combination ., As gene regulatory circuits are often involved in crucial biological processes , their functions should be robust to perturbation ., We therefore asked whether the constraints imposed by multifunctionality also impact the robustness of circuits and their functions ., In studying robustness , we differentiate between the robustness of a genotype ( circuit ) and the robustness of a k-function ., We assess the robustness of a genotype as the proportion of all possible single-mutants that have the same k-function , and the robustness of a k-function as the average robustness of all genotypes with that k-function 17 , 18 , 51 , 63 ( Material and Methods ) ., We refer to the collection of genotypes with a given k-function as a genotype set , which may comprise one or more genotype networks ., We emphasize that a genotype may be part of several different genotype sets , because genotypes typically have more than one k-function ., Fig . 3A shows that the robustness of a k-function decreases approximately linearly as k increases , indicating a trade-off between multifunctionality and robustness ., However , some degree of robustness is maintained so long as ., For larger k , some functions exist that have zero robustness ( Text S1 ) , that is , none of the circuits with these functions can tolerate a change in their regulatory genotype ., The inset of Fig . 3A reveals a similar inverse relationship between the size of a genotype set and the number of functions k , implying that multifunctions become increasingly less “designable” 64 — fewer circuits have them — as k increases ( Text S1 ) ., For example , for as few as functions , the genotype set may comprise a single genotype , reducing the corresponding robustness of the k-function to zero ., For each value of k , the maximum proportion of genotypes with a given k-function is equal to the square of the maximum proportion of genotypes with a function , explaining the triangular shape of the data in the inset ., This triangular shape indicates that the genotype set of a given k-function is always smaller than the union of the k constituent genotypes sets ., Additionally , we find that the robustness of a k-function and the size of its genotype set are strongly correlated ( Fig . S1 ) , indicating that the genotypes of larger genotype sets are , on average , more robust than those of smaller genotype sets ., This result is not trivial because the structure of a genotype set may change with its size ., For example , large genotype sets may comprise many isolated genotypes , or their genotype networks might be structured as long linear chains ., In either case , the robustness of a k-function would decrease as the size of its genotype set increased ., We have so far focused on the properties of the genotype sets of k-functions , but have not considered the properties of the genotype networks that make up these sets ., Therefore , we next asked how genotypic robustness varies across the genotype networks of k-functions ., In Figs ., 3B–D , we show the distributions of genotypic robustness for representative genotype networks with functions ., These distributions highlight the inherent variability in genotypic robustness that is present in the genotype networks of multifunctions , indicating that genotypic robustness is an evolvable property of multifunctional circuits ., Indeed , in Fig . S2 , we show the results of random walks on these genotype networks , which confirm that it is almost always possible to increase genotypic robustness through a series of mutational steps that preserve the k-function ., In Fig . S3 , we show in which dynamic regimes ( Material and Methods ) the circuits in these same genotype networks lie ., We have shown that the genotype set of any k-function is non-empty ( Fig . 2 ) , meaning that there are no “prohibited” k-functions ., We now ask how the genotypes with a given k-function are organized in genotype space ., More specifically , is it possible to connect any two circuits with the same k-function through a sequence of small genotypic changes where each change in the sequence preserves this k-function ?, In other words , are all genotypes with a given k-function part of the same genotype network , or do such genotypes occur on multiple disconnected genotype networks ?, Fig . 4 shows the relationship between the number of genotype networks in a genotype set and the number of circuit functions k ., For monofunctional circuits ( ) , the genotype set always consists of a single , connected genotype network ., This implies that any genotype in the genotype set can be reached from any other via a series of function-preserving mutational events ., In contrast , for circuits with functions , the genotype set often fragments into several isolated genotype networks , indicating that some regions of the genotype set cannot be reached from some others without jeopardizing circuit function ., The most extreme fragmentation occurs for functions , where some genotype sets break up into more than 20 isolated genotype networks ., Fig . S4 provides a schematic illustration of how fragmentation can occur in a k-functions genotype set , despite the fact that the genotype sets of the k constituent monofunctions consist of genotype networks that are themselves connected ., Fig . S5 provides a concrete example of fragmentation , depicting one genotype from each of the several genotype networks of a bifunctions genotype set ., The proportion of k-functions with genotype sets that comprise a single genotype network is shown in the inset of Fig . 4 ., This proportion decreases dramatically as the number of functions increases from to , such that only 16% of genotype sets comprise a single genotype network when ., Figs ., 4B–D show that the distributions of the number of genotype networks per genotype set are typically left-skewed ., This implies that when fragmentation occurs , the genotype set usually fragments into only a few genotype networks ., However , the distribution of genotype network sizes across all genotype sets is heavy-tailed and often spans several orders of magnitude ( Fig . S6 ) ., This means that the number of genotypes per genotype network is highly variable ., We next ask whether the number of genotypes in the genotype set of a k-function can be predicted from the number of genotypes in the genotype sets of the k constituent monofunctions ., To address this question , we define the fractional size of a genotype set as the number of genotypes in the set , divided by the number of genotypes in genotype space ., We first observe that the maximum fractional size of a genotype set of a k-function is equal to ( Fig . S6 ) , which is the maximum fractional size of a genotype set for monofunctional circuits 44 raised to the kth power ., In general , we find that the fractional size of a genotype set of a k-function can be approximated with reasonable accuracy by the product of the fractional sizes of the genotype sets of the k constituent monofunctions , but that the accuracy of this approximation decreases as k increases ( Fig . S7 ) ., While these fractional genotype set sizes may be quite small , we note that their absolute sizes are still fairly large , even in the tiny circuits considered here ., For example , for functions the maximum genotype set size is 262 , 144 ., For functions , the maximum is 32 , 768 ., In evolution , a circuit may acquire a new regulatory function while preserving its pre-existing functions ., An example is the highly-conserved hedgehog regulatory circuit , which patterns the insect wing blade ., In butterflies , this regulatory circuit has acquired a new function ., It helps form the wings eyespots , an antipredatory adaptation that arose after the insect body plan 65 ., This example illustrates that a regulatory circuit may acquire additional functions incrementally via gradual genetic change ., The order in which the mutations leading to a new function arise and go to fixation can have a profound impact upon the evolution of such phenotypes 66 ., In particular , early mutations have the potential to influence the phenotypic effects of later mutations , which can lead to a phenomenon known as historical contingency ., We next ask whether it is possible for a circuit to incrementally evolve regulatory functions in any order , or whether this evolutionary process is susceptible to historical contingency ., In other words , is it possible that some sequence of genetic changes that lead a circuit to have k functions also preclude it from gaining an additional function ?, The genotype space framework allows us to address this question in a systematic way , because it permits us to see contingency as a result of genotype set fragmentation ., Specifically , contingency means that , as a result of fragmentation , the genotype network of a new function may become inaccessible from at least one of the genotype networks of a k-functions genotype set ., To ask whether this occurs in our model regulatory circuits , we considered all permutations of every k-function ., These permutations reflect every possible order in which a circuit may acquire a specific combination of k functions through a sequence of genetic changes ., To determine the frequency with which historical contingency occurs , we calculate the number of genotype networks per genotype set , as the k functions are incrementally added ., This procedure is outlined in Fig . S4 and detailed in the Material and Methods section ., We note that historical contingency is not possible when because all monofunctions comprise genotype sets with a single connected genotype network ., Historical contingency is also not possible when , because there is only one genotype that yields this combination ( Fig . 2 ) ., In Fig . 5 , we show the relationship between the proportion of k-functions that exhibit historical contingency and the number of functions k ., For as few as functions , 43% of all k-functions exhibit historical contingency ., This percentage is highest for , where 94% of combinations are contingent ., The inset of Fig . 5 shows the proportion of the permutations of a k-function in which genotype set fragmentation may preclude the evolution of the k-function ., Again , this proportion is highest for functions ., These results highlight an additional constraint of multifunctionality ., Not only does the number of genotypes with k functions decrease as k increases , but the dependence upon the temporal order in which these functions evolve tends to increase ., In the Supporting Online Material , we repeat the above calculations to show how our results scale to equilibrium expression states with period ( For the sake of computational tractability , we restrict our attention to the case where all equilibrium expression states have the same period P ) ., We show that the exponential decrease in the number of circuits with k functions also holds for periodic equilibrium expression states , but that the maximum number of functions per circuit decreases with increasing ( Fig . S8 ) ., So long as , it is possible for a circuit to have more than one function ., In this case , the inverse relationship between robustness to genetic perturbation and the number of functions k also holds ( Fig . S9 ) ., Similarly , the results pertaining to genotype set fragmentation hold so long as ( Fig . S10 ) ., Lastly , the results pertaining to historical contingency only hold when ., This is because it is not possible for a circuit with an equilibrium expression pattern of period to have more than functions , which is a prerequisite for historical contingency ( Material and Methods ) ., Taken together , these additional observations show that the results obtained for fixed-point equilibrium expression states can also apply to periodic equilibrium expression states , so long as is not too large ., We have used a Boolean model of gene regulatory circuits to exhaustively characterize the functions of all possible combinations of circuit topologies and signal-integration functions in three-gene circuits ., The most basic question we have addressed is whether multifunctionality is easy or difficult to attain in regulatory circuits ., Our results show that while the number of circuits with k functions decreases sharply as k increases , there are generally thousands of circuits with k functions , so long as k is not exceedingly large ., Thus , multifunctionality is relatively easy to attain , even in the tiny circuits examined here ., It is worth considering how this result might translate to larger circuits ., In a related model of gene regulatory circuits with genes , the genotype sets of bifunctions comprised an average of circuits 21 , which is over an order of magnitude more circuits per bifunction than observed here ( Fig . 3 , inset ) ., For a greater number of functions k , we expect the number of circuits per k-function to increase as the number of genes N in the regulatory circuit increases ., This is because the maximum number of circuits with a given k-function is , which is the total number of circuits with N genes ( ) multiplied by the maximum proportion of circuits per multifunction ( ) ., For a given number of functions k , this quotient will increase hyper-exponentially as N increases , indicating a dramatic increase in the maximum number of circuits per k-function ., More generally , because the fractional size of a k-functions genotype set can be approximated as the product of the fractional sizes of the genotype sets of its k constituent monofunctions ( Fig . S7 ) and because the total number of circuits increases exponentially with N , our observation that there are many circuits with k functions is expected to scale to larger circuits ., The next question we asked is whether there is a tradeoff between the robustness of a k-function and the number of functions k ., We found that the robustness of a k-function decreases as k increases ., However , some degree of robustness is generally maintained , so long as k is not too large ., These observations suggest that the number of circuit functions generally does not impose severe constraints on the evolution of circuit genotypes , unless the number of functions is very large ., Our current knowledge of biological circuits is too limited to allow us to count the number of functions per circuit ., However , we can ask whether the functional “burden” on biological circuits is very high ., If so , we would expect that the genes that form these circuits and their regulatory regions cannot tolerate genetic perturbations , and that they have thus accumulated few or no genetic changes in their evolutionary history ., However , this is not the case ., The biochemical activities and regulatory regions of circuit genes can diverge extensively without affecting circuit function 55 , 59 , 61 , 67 , and the very different circuit architectures of distantly related species can have identical function 24 , 28 ., Further , circuits are highly robust to the experimental perturbation of their architecture , such as the rewiring of regulatory interactions 20 ., More indirect evidence comes from the study of genes with multiple functions , identified through gene ontology annotations ., The rate of evolution of these genes is significantly but only weakly correlated with the number of known functions 68 ., Thus , the functional burden on biological genes and circuits is not sufficiently high to preclude evolutionary change ., Previous studies of monofunctional regulatory circuits have revealed broad distributions of circuit robustness to genetic perturbation 16 , 18 , 29 ., We therefore asked if this is also the case for multifunctional circuits ., We found that circuit robustness was indeed variable , but that the mean and variance of the distributions of circuit robustness decreased as the number of functions k increased ., Thus , variation in circuit robustness persists in multifunctional circuits , so long as k is not too large ., This provides further evidence that robustness to mutational change may be considered the rule , rather than the exception , in biological networks 1 , 18 , 20 , 29 ., However , to make the claim that robustness to genetic perturbation is an evolvable property in multifunctional regulatory circuits requires not only variability in circuit robustness , but also the ability to change one circuit into another via a series of mutations that do not affect any of the circuits functions ., We therefore asked whether it is possible to interconvert any two circuits with the same function via a series of function-preserving mutational changes ., We showed that this is always possible for monofunctions , but not necessarily for multifunctions , because these often comprise fragmented genotype sets ., Genotype set fragmentation has also been observed at lower levels of biological organization , such as the mapping from RNA sequence to secondary structure 69 ., Such fragmentation has two evolutionary implications , as has recently been discussed for RNA phenotypes 70 ., First , the mutational robustness of a phenotype ( function ) depends upon which genotype network its sequences inhabit , as we have also shown for regulatory circuits ( Fig . S11 ) ., Second , it can lead to historical contingency , where the phenotypic effects of future mutations depend upon the current genetic background ., Such contingency indeed occurs in our circuits , because the specific genotype network that a circuit ( genotype ) occupies may be influenced by the temporal order in which a circuits functions ( phenotypes ) have evolved ., This order in turn may affect a circuits ability to evolve new functions ., These observations hinge on the assumption that the space between two ( disconnected ) parts of a fragmented genotype set is not easily traversed ., For example , in RNA it is well known that pairs of so-called compensatory mutations can allow transitions between genotype networks 71 , thus alleviating the historical contingency caused by fragmentation ., To assess whether an analogous phenomenon might exist for regulatory circuits , we calculated the average distance between all pairs of genotypes on distinct genotype networks for circuits with the same k-function ., We found that this distance decreases as the number of functions k increases , indicating an increased proximity between genotype networks ( Fig . S12 ) ., However , those pairs of genotypes in any two different genotype networks that had the minimal distance of two mutations never exceeded 1% of all pairs of genotypes on these networks , and was as low as 0 . 03% for functions ( Fig . S12A , inset ) ., This means that transitions between genotype networks through few mutations are not usually possible in these model regulatory circuits ., Thus , the multiple genotype networks of a genotype set can indee
Gene regulatory circuits drive the development , physiology , and behavior of organisms from bacteria to humans ., The phenotypes or functions of such circuits are embodied in the gene expression patterns they form ., Regulatory circuits are typically multifunctional , forming distinct gene expression patterns in different embryonic stages , tissues , or physiological states ., Any one circuit with a single function can be realized by many different regulatory genotypes ., Multifunctionality presumably constrains this number , but we do not know to what extent ., We here exhaustively characterize a genotype space harboring millions of model regulatory circuits and all their possible functions ., As a circuits number of functions increases , the number of genotypes with a given number of functions decreases exponentially but can remain very large for a modest number of functions ., However , the sets of circuits that can form any one set of functions becomes increasingly fragmented ., As a result , historical contingency becomes widespread in circuits with many functions ., Whether a circuit can acquire an additional function in the course of its evolution becomes increasingly dependent on the function it already has ., Circuits with many functions also become increasingly brittle and sensitive to mutation ., These observations are generic properties of a broad class of circuits and independent of any one circuit genotype or phenotype .
Many essential biological processes , ranging from embryonic patterning to circadian rhythms , are driven by gene regulatory circuits , which comprise small sets of genes that turn each other on or off to form a distinct pattern of gene expression ., Gene regulatory circuits often have multiple functions ., This means that they can form different gene expression patterns at different times or in different tissues ., We know little about multifunctional gene regulatory circuits ., For example , we do not know how multifunctionality constrains the evolution of such circuits , how many circuits exist that have a given number of functions , and whether tradeoffs exist between multifunctionality and the robustness of a circuit to mutation ., Because it is not currently possible to answer these questions experimentally , we use a computational model to exhaustively enumerate millions of regulatory circuits and all their possible functions , thereby providing the first comprehensive study of multifunctionality in model regulatory circuits ., Our results highlight limits of circuit designability that are relevant to both systems biologists and synthetic biologists .
systems biology, computer science, genetics, biology, computational biology, computerized simulations, gene networks