Context-dependent Binding Research Articles

FuzPred: a web server for the sequence-based prediction of the context-dependent binding modes of proteins.

Proteins form complex interactions in the cellular environment to carry out their functions. They exhibit a wide range of binding modes depending on the cellular conditions, which result in a variety of ordered or disordered assemblies. To help rationalise the binding behavior of proteins, the FuzPred server predicts their sequence-based binding modes without specifying their binding partners. The binding mode defines whether the bound state is formed through a disorder-to-order transition resulting in a well-defined conformation, or through a disorder-to-disorder transition where the binding partners remain conformationally heterogeneous. To account for the context-dependent nature of the binding modes, the FuzPred method also estimates the multiplicity of binding modes, the likelihood of sampling multiple binding modes. Protein regions with a high multiplicity of binding modes may serve as regulatory sites or hot-spots for structural transitions in the assembly. To facilitate the interpretation of the predictions, protein regions with different interaction behaviors can be visualised on protein structures generated by AlphaFold. The FuzPred web server (https://fuzpred.bio.unipd.it) thus offers insights into the structural and dynamical changes of proteins upon interactions and contributes to development of structure-function relationships under a variety of cellular conditions.

A direct comparison between AML1-ETO and ETO2-GLIS2 leukemia fusion proteins reveals context-dependent binding and regulation of target genes and opposite functions in cell differentiation.

The ETO-family transcriptional corepressors, including ETO, ETO2, and MTGR1, are all involved in leukemia-causing chromosomal translocations. In every case, an ETO-family corepressor acquires a DNA-binding domain (DBD) to form a typical transcription factor—the DBD binds to DNA, while the ETO moiety manifests transcriptional activity. A directly comparative study of these “homologous” fusion transcription factors may clarify their similarities and differences in regulating transcription and leukemogenesis. Here, we performed a side-by-side comparison between AML1-ETO and ETO2-GLIS2, the most common fusion proteins in M2-and M7-subtypes of acute myeloid leukemia, respectively, by inducible expression of them in U937 leukemia cells. We found that, although AML1-ETO and ETO2-GLIS2 can use their own DBDs to bind DNA, they share a large proportion of genome-wide binding regions dependent on other cooperative transcription factors, including the ETS-, bZIP- and bHLH-family proteins. AML1-ETO acts as either transcriptional repressor or activator, whereas ETO2-GLIS2 mainly acts as activator. The repressor-versus-activator functions of AML1-ETO might be determined by the abundance of cooperative transcription factors/cofactors on the target genes. Importantly, AML1-ETO and ETO2-GLIS2 differentially regulate key transcription factors in myeloid differentiation including PU.1 and C/EBPβ. Consequently, AML1-ETO inhibits, but ETO2-GLIS2 facilitates, myeloid differentiation of U937 cells. This function of ETO2-GLIS2 is reminiscent of a similar effect of MLL-AF9 as previously reported. Taken together, this directly comparative study between AML1-ETO and ETO2-GLIS2 in the same cellular context provides insights into context-dependent transcription regulatory mechanisms that may underlie how these seemingly “homologous” fusion transcription factors exert distinct functions to drive different subtypes of leukemia.

HnRNP A1/A2 Proteins Assemble onto 7SK snRNA via Context Dependent Interactions

7SK small nuclear RNA (snRNA) is an abundant and ubiquitously expressed noncoding RNA that functions to modulate the activity of RNA Polymerase II (RNAPII) in part by stabilizing distinct pools of 7SK-protein complexes. Prevailing models suggest that the secondary structure of 7SK is dynamically remodeled within its alternative RNA-protein pools such that its architecture differentially regulates the exchange of cognate binding partners. The nuclear hnRNP A1/A2 proteins influence the biology of 7SK snRNA via processes that require an intact stem loop (SL) 3 domain; however, the molecular details by which hnRNPs assemble onto 7SK snRNA are yet to be described. Here, we have taken an integrated approach to present a detailed description of the 7SK-hnRNP A1 complex. We show that unbound 7SK snRNA adopts at least two major conformations in solution, with significant structural differences localizing to the SL2-3 linker and the base of SL3. Phylogenetic analysis indicates that this same region is the least genetically conserved feature of 7SK snRNA. By performing DMS modifications with the presence of excess protein, we reveal that hnRNP A1 binds with selectivity to SL3 through mechanisms that increase the flexibility of the RNA adjacent to putative binding sites. Calorimetric titrations further validate that hnRNP A1-SL3 assembly is complex with the affinity of discrete binding events modulated by the surrounding RNA structure. To interpret this context-dependent binding phenomenon, we determined a 3D model of SL3 to show that it folds to position minimal hnRNP A1/A2 binding sites (5′-Y/RAG-3′) within different local environments. SL3-protein complexes resolved by SEC-MALS-SAXS confirm that up to four hnRNP A1 proteins bind along the entire surface of SL3 via interactions that preserve the overall structural integrity of this domain. In sum, the collective results presented here reveal a specific role for a folded SL3 domain to scaffold hnRNP A1/A2-7SK assembly via mechanisms modulated by the surrounding RNA structure.

Sequence-based prediction of protein binding mode landscapes

Interactions between disordered proteins involve a wide range of changes in the structure and dynamics of the partners involved. These changes can be classified in terms of binding modes, which include disorder-to-order (DO) transitions, when proteins fold upon binding, as well as disorder-to-disorder (DD) transitions, when the conformational heterogeneity is maintained in the bound states. Furthermore, systematic studies of these interactions are revealing that proteins may exhibit different binding modes with different partners. Proteins that exhibit this context-dependent binding can be referred to as fuzzy proteins. Here we investigate amino acid code for fuzzy binding in terms of the entropy of the probability distribution of transitions towards decreasing order. We implement these entropy calculations into the FuzPred (http://protdyn-fuzpred.org) algorithm to predict the range of context-dependent binding modes of proteins from their amino acid sequences. As we illustrate through a variety of examples, this method identifies those binding sites that are sensitive to the cellular context or post-translational modifications, and may serve as regulatory points of cellular pathways.

Sequence-Based Prediction of Fuzzy Protein Interactions

It is becoming increasingly recognised that disordered proteins may be fuzzy, in that they can exhibit a wide variety of binding modes. In addition to the well-known process of folding upon binding (disorder-to-order transition), many examples are emerging of interacting proteins that remain disordered in their bound states (disorder-to-disorder transitions). Furthermore, disordered proteins may populate ordered and disordered states to different extents depending on their partners (context-dependent binding). Here we assemble three datasets comprising disorder-to-order, context-dependent, and disorder-to-disorder transitions of 828 protein regions represented in 2157 complexes and elucidate the sequence-determinants of the different interaction modes. We found that fuzzy interactions originate from local sequence compositions that promote the sampling of a wide range of different structures. Based on this observation, we developed the FuzPred method (http://protdyn-fuzpred.org) of predicting the binding modes of disordered proteins based on their amino acid sequences, without specifying their partners. We thus illustrate how the amino acid sequences of proteins can encode a wide range of conformational changes upon binding, including transitions from disordered to ordered and from disordered to disordered states.

Conformational Dynamics and Cooperativity Drive the Specificity of a Protein-Ligand Interaction

Molecular recognition is critical for the fidelity of signal transduction in biology. Conversely, the disruption of protein-protein interactions can lead to disease. Thus, comprehension of the molecular determinants of specificity is essential for understanding normal biological signaling processes and for the development of precise therapeutics. Although high-resolution structures have provided atomic details of molecular interactions, much less is known about the influence of cooperativity and conformational dynamics. Here, we used the Tiam2 PSD-95/Dlg/ZO-1 (PDZ) domain and a quadruple mutant (QM), engineered by swapping the identity of four residues important for specificity in the Tiam1 PDZ into the Tiam2 PDZ domain, as a model system to investigate the role of cooperativity and dynamics in PDZ ligand specificity. Surprisingly, equilibrium binding experiments found that the ligand specificity of the Tiam2 QM was switched to that of the Tiam1 PDZ. NMR-based studies indicated that Tiam2 QM PDZ, but not other mutants, had extensive microsecond to millisecond motions distributed throughout the entire domain suggesting structural cooperativity between the mutated residues. Thermodynamic analyses revealed energetic cooperativity between residues in distinct specificity subpockets that was dependent upon the identity of the ligand, indicating a context-dependent binding mechanism. Finally, isothermal titration calorimetry experiments showed distinct entropic signatures along the mutational trajectory from the Tiam2 wild-type to the QM PDZ domain. Collectively, our studies provide unique insights into how structure, conformational dynamics, and thermodynamics combine to modulate ligand-binding specificity and have implications for the evolution, regulation, and design of protein-ligand interactions.

Chemoenzymatic Synthesis of N-glycan Positional Isomers and Evidence for Branch Selective Binding by Monoclonal Antibodies and Human C-type Lectin Receptors.

Here, we describe a strategy for the rapid preparation of pure positional isomers of complex N-glycans to complement an existing array comprising a larger number of N-glycans and smaller glycan structures. The expanded array was then employed to study context-dependent binding of structural glycan fragments by monoclonal antibodies and C-type lectins. A partial enzymatic elongation of semiprotected core structures was combined with the protecting-group-aided separation of positional isomers by preparative HPLC. This methodology, which avoids the laborious chemical differentiation of antennae, was employed for the preparation of eight biantennary N-glycans with Galβ1,4GlcNAc (LN), GalNAcβ1,4GlcNAc (LDN), and GalNAcβ1,4[Fucα1,3]GlcNAc (LDNF) motifs presented on either one or both antennae. Screening of the binding specificities of three anti-LeX monoclonal IgM antibodies raised against S. mansoni glycans and three C-type lectin receptors of the innate immune system, namely DC-SIGN, DC-SIGNR, and LSECtin, revealed a surprising context-dependent fine specificity for the recognition of the glycan motifs. Moreover, we observed a striking selection of one individual positional isomer over the other by the C-type lectins tested, underscoring the biological relevance of the structural context of glycan elements in molecular recognition.

Neuronal oscillations: unavoidable and useful?

Neuronal systems have a high propensity to engage in oscillatory activity because both the properties of individual neurons and canonical circuit motifs favour rhythmic activity. In addition, coupled oscillators can engage in a large variety of dynamical regimes, ranging from synchronization with different phase offsets to chaotic behaviour. Which regime prevails depends on differences between preferred oscillation frequencies, coupling strength and coupling delays. The ability of delay coupled oscillator networks to generate a rich repertoire of temporally structured activation sequences is exploited by central pattern generator networks for the control of movements. However, it is less clear whether temporal patterning of neuronal discharges also plays a role in cognitive processes. Here, it will be argued that the temporal patterning of neuronal discharges emerging from delay coupled oscillator networks plays a pivotal role in all instances in which selective relations have to be established between the responses of distributed assemblies of neurons. Examples are the dynamic formation of functional networks, the selective routing of activity in densely interconnected networks, the attention-dependent selection of sensory signals, the fast and context-dependent binding of responses for further joint processing in pattern recognition and the formation of associations by learning. Special consideration is given to arguments that challenge a functional role of oscillations and synchrony in cognition because of the volatile nature of these phenomena and recent evidence will be reviewed suggesting that this volatility is functionally advantageous.

Abstract 391: A comparative genomics approach to understanding the control of cell context dependent P53 binding

Abstract Cancer is considered to be a genetic disease characterized by sequential accumulation of mutations. Recent studies have provided evidence that epigenetic changes and regulatory sequence mutations could also dysregulate oncogenes and tumor suppressors. As a transcription factor, P53 is activated in response to oncogenic stress and exerts distinct anti-proliferative functions based on the stressor and cell type. Though a number of ChIP-Seq studies have identified thousands of P53 binding sites in mammalian genomes, the functionality of these binding sites remains to be established. In addition, we know little about what controls stress and cell context specific binding profile of P53. Traditionally, mutations in the coding regions of p53 have been extensively studied to gain insights on its role in cancer and to identify strategies to restore the functions of p53 in cancer cells. We hypothesize that mutations in or epigenetic silencing of functionally important P53 binding sites play an important role in tumorigenesis as well. Since functional regulatory regions tend to be more evolutionarily conserved, in this project we propose using a comparative genomics approach to identify functional P53 binding sites and determine if these regions are involved in tumorigenesis. In order to study how P53 binding following DNA damage differs between cell types, we aim to perform comparison between P53 ChIP-Seq data generated in our lab from Drosophila embryos at different developmental stages and a cell line (Kc167) as well as published datasets from mammalian stem cells and differentiated cells. Differential expression analysis using RNA-Seq exhibited that at an early stem cell-like stage there is P53-dependent induction of pro-apoptotic genes in response to DNA damage but not in the differentiated stages. We are seeking to establish functionally significant P53 binding sites by combining ChIP-Seq data and gene expression data obtained using RNA-Seq. The importance of these binding sites will be verified by CRISPR-Cas9 -mediated genome editing. We will also perform comparison studies of DNA damage-induced P53 binding in human, mouse and Drosophila to identify analogous patterns. Knowledge gained from this study will help us to understand the role of non-coding regulatory regions in tumorigenesis, and predict patient response to apoptosis inducing therapeutic agents. It may also lead to novel strategies to restore cellular sensitivity to chemotherapy or radiation. Citation Format: Varsha Sundaresan, Ying Li, Benedetto DiCiaccio, Victor T. Lin, Lei Zhou. A comparative genomics approach to understanding the control of cell context dependent P53 binding [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 391. doi:10.1158/1538-7445.AM2017-391

The expanding role of the Ehmt2/G9a complex in neurodevelopment

Epigenetic regulators play a crucial role in neurodevelopment. One such epigenetic complex, Ehmt1/2 (G9a/GLP), is essential for repressing gene transcription by methylating H3K9 in a highly tissue- and temporal-specific manner. Recently, data has emerged suggesting that this complex plays additional roles in regulating the activity of numerous other non-histone proteins. While much is known about the downstream effects of Ehmt1/2 function, evidence is only beginning to come to light suggesting the control of Ehmt1/2 function may be, at least in part, due to context-dependent binding partners. Here we review emerging roles for the Ehmt1/2 complex suggesting that it may play a much larger role than previously recognized, and discuss binding partners that we and others have recently characterized which act to coordinate its activity during early neurodevelopment.

Co-regulation in embryonic stem cells via context-dependent binding of transcription factors

With the accumulation of genome-wide binding data for many transcription factors (TFs) in the same cell type or cellular condition, it is of great current interest to systematically infer the complex regulatory logic among multiple TFs. In particular, ChIP-Seq data have been generated for 14 core TFs critical to the maintenance and reprogramming of mouse embryonic stem cells (ESCs). This provides a great opportunity to study the regulatory collaboration and interaction among these TFs and with other unknown co-regulators. In combination with liquid association among gene expression profiles, we develop a computational method to predict context-dependent (CD) co-egulators of these core TFs in ESCs from pairwise binding datasets. That is, co-occupancy between a core TF and a predicted co-regulator depends on the presence or absence of binding sites of another core TF, which is regarded as a binding context. Unbiased external validation confirms that the predicted CD binding of a co-regulator is reliable. Our results reveal a detailed CD co-regulation network among the 14 core TFs and provide many other potential co-regulators showing strong agreement with the literature. See www.stat.ucla.edu/~zhou/CMF for software and source code.

클라이언트 매시업의 편리한 이동 사용자 인터페이스 개발

Client based web mashups have become one of the most important architecture in web application development. Although there are well known methods to generate presentation view codes for web services, navigations between views and service requests are still developed manually in most web mashups. In this paper, we propose the concept of deterministic data binding from output data to input parameters of another method. Using binding relations, we can model navigation menus for service method requests. For a given set of data mappings between services, we investigate context dependent binding conditions and discuss the generation of views and navigation menus. The proposed approach provides UI for users to navigate services of client mashup page using simple and convenient interface even when the number of services and the size of the mashup page grows. In order to show the usability of the proposed approach, we present a historic tourism service.

The Specificity and Topology of Chromatin Interaction Pathways in Yeast

Packaging of DNA into chromatin has a profound impact on gene expression. To understand how changes in chromatin influence transcription, we analyzed 165 mutants of chromatin machinery components in Saccharomyces cerevisiae. mRNA expression patterns change in 80% of mutants, always with specific effects, even for loss of widespread histone marks. The data are assembled into a network of chromatin interaction pathways. The network is function based, has a branched, interconnected topology, and lacks strict one-to-one relationships between complexes. Chromatin pathways are not separate entities for different gene sets, but share many components. The study evaluates which interactions are important for which genes and predicts additional interactions, for example between Paf1C and Set3C, as well as a role for Mediator in subtelomeric silencing. The results indicate the presence of gene-dependent effects that go beyond context-dependent binding of chromatin factors and provide a framework for understanding how specificity is achieved through regulating chromatin.

Sox17 modulates Wnt3A/β-catenin-mediated transcriptional activation of the Lef-1 promoter

Wnt/β-catenin-dependent activation of lymphoid enhancer factor 1 (Lef-1) plays an important role in numerous developmental processes. In this context, transcription of the Lef-1 gene is increased by Wnt-mediated TCF4/β-catenin activation on the Lef-1 promoter through mechanisms that remain poorly defined. In mouse airway submucosal gland progenitor cells, Wnt3A transiently induces Lef-1 gene expression, and this process is required for epithelial cell proliferation and glandular morphogenesis. In the present study, we sought to identify additional candidate transcriptional regulators of the Lef-1 gene during glandular morphogenesis. To this end, we found that Sox17 expression is dramatically downregulated in early glandular progenitor cells that induce Lef-1 expression. Wnt stimulation of undifferentiated primary airway epithelial cells induced similar changes in Sox17 and Lef-1 expression. Reporter assays revealed that ectopic expression of Sox17 suppresses Wnt3A/β-catenin activation of the Lef-1 promoter in cell lines. EMSA and ChIP analyses defined several Sox17- and TCF4-binding sites that collaborate in transcriptional control of the Lef-1 promoter. More specifically, Sox17 bound to four sites in the Lef-1 promoter, either directly or indirectly through TCF complexes. The DNA- or β-catenin-binding domains of Sox17 controlled context-specific binding of Sox17/TCF complexes on the Lef-1 promoter. Combinatorial site-directed mutagenesis of Sox17- or TCF-binding sites in the Lef-1 promoter demonstrated that these sites control Wnt/β-catenin-mediated induction and/or repression. These findings demonstrate for the first time that Sox17 can directly regulate Wnt/β-catenin-dependent transcription of the Lef-1 promoter and reveal new context-dependent binding sites in the Lef-1 promoter that facilitate protein-protein interactions between Sox17 and TCF4.

Specificity landscapes of DNA binding molecules elucidate biological function

Evaluating the specificity spectra of DNA binding molecules is a nontrivial challenge that hinders the ability to decipher gene regulatory networks or engineer molecules that act on genomes. Here we compare the DNA sequence specificities for different classes of proteins and engineered DNA binding molecules across the entire sequence space. These high-content data are visualized and interpreted using an interactive "specificity landscape" which simultaneously displays the affinity and specificity of a million-plus DNA sequences. Contrary to expectation, specificity landscapes reveal that synthetic DNA ligands match, and often surpass, the specificities of eukaryotic DNA binding proteins. The landscapes also identify differential specificity constraints imposed by diverse structural folds of natural and synthetic DNA binders. Importantly, the sequence context of a binding site significantly influences binding energetics, and utilizing the full contextual information permits greater accuracy in annotating regulatory elements within a given genome. Assigning such context-dependent binding values to every DNA sequence across the genome yields predictive genome-wide binding landscapes (genomescapes). A genomescape of a synthetic DNA binding molecule provided insight into its differential regulatory activity in cultured cells. The approach we describe will accelerate the creation of precision-tailored DNA therapeutics and uncover principles that govern sequence-specificity of DNA binding molecules.

Distributed processing and temporal codes in neuronal networks

The cerebral cortex presents itself as a distributed dynamical system with the characteristics of a small world network. The neuronal correlates of cognitive and executive processes often appear to consist of the coordinated activity of large assemblies of widely distributed neurons. These features require mechanisms for the selective routing of signals across densely interconnected networks, the flexible and context dependent binding of neuronal groups into functionally coherent assemblies and the task and attention dependent integration of subsystems. In order to implement these mechanisms, it is proposed that neuronal responses should convey two orthogonal messages in parallel. They should indicate (1) the presence of the feature to which they are tuned and (2) with which other neurons (specific target cells or members of a coherent assembly) they are communicating. The first message is encoded in the discharge frequency of the neurons (rate code) and it is proposed that the second message is contained in the precise timing relationships between individual spikes of distributed neurons (temporal code). It is further proposed that these precise timing relations are established either by the timing of external events (stimulus locking) or by internal timing mechanisms. The latter are assumed to consist of an oscillatory modulation of neuronal responses in different frequency bands that cover a broad frequency range from <2 Hz (delta) to >40 Hz (gamma) and ripples. These oscillations limit the communication of cells to short temporal windows whereby the duration of these windows decreases with oscillation frequency. Thus, by varying the phase relationship between oscillating groups, networks of functionally cooperating neurons can be flexibly configurated within hard wired networks. Moreover, by synchronizing the spikes emitted by neuronal populations, the saliency of their responses can be enhanced due to the coincidence sensitivity of receiving neurons in very much the same way as can be achieved by increasing the discharge rate. Experimental evidence will be reviewed in support of the coexistence of rate and temporal codes. Evidence will also be provided that disturbances of temporal coding mechanisms are likely to be one of the pathophysiological mechanisms in schizophrenia.

Structural Basis of Affinity Maturation and Intramolecular Cooperativity in a Protein-Protein Interaction

Although protein-protein interactions are involved in nearly all cellular processes, general rules for describing affinity and selectivity in protein-protein complexes are lacking, primarily because correlations between changes in protein structure and binding energetics have not been well determined. Here, we establish the structural basis of affinity maturation for a protein-protein interaction system that we had previously characterized energetically. This model system exhibits a 1500-fold affinity increase. Also, its affinity maturation is restricted by negative intramolecular cooperativity. With three complex and six unliganded variant X-ray crystal structures, we provide molecular snapshots of protein interface remodeling events that span the breadth of the affinity maturation process and present a comprehensive structural view of affinity maturation. Correlating crystallographically observed structural changes with measured energetic changes reveals molecular bases for affinity maturation, intramolecular cooperativity, and context-dependent binding.

Pax3 target gene recognition occurs through distinct modes that are differentially affected by disease‐associated mutations

The paired box protein Pax3 is an essential regulator of muscle and neural crest-derived cell types, including melanocytes. Within this lineage, Pax3 has been shown to regulate the genes encoding microphthalmia-associated transcription factor (Mitf) and tyrosinase-related protein-1 (Trp-1), despite each having dissimilar Pax3 recognition sequences. We have, therefore, examined the structural requirements for Pax3 binding to the MITF and TRP-1 promoter elements, focusing on the contribution of the paired domain and homeodomain to Pax3 target site recognition. Unexpectedly, although the MITF element is characterized by suboptimal recognition motifs for the paired domain and homeodomain, it sustains a higher level of Pax3 binding than TRP-1, which contains a canonical paired domain site. The basis for this difference involves a context-dependent cooperative binding event requiring both the paired domain and homeodomain, while the paired domain alone is sufficient for TRP-1 recognition. Significantly, the analysis of Waardenburg syndrome mutations reveals marked disparity in their effects on MITF and TRP-1 binding that further underscores mechanistic differences in their interaction with Pax3. Importantly, these mutations also exert distinct effects on the ability of Pax3 to regulate reporter genes fused to either the MITF or TRP-1 promoters. Our results, therefore, establish that Pax3 can regulate target genes through alternate modes of DNA recognition that are differentially impacted by disease-causing mutations, which together have important implications for understanding Pax3-regulated gene networks.

On the role of biophysical properties of cortical neurons in binding and segmentation of visual scenes.

Neuroscience is progressing vigorously, and knowledge at different levels of description is rapidly accumulating. To establish relationships between results found at these different levels is one of the central challenges. In this simulation study, we demonstrate how microscopic cellular properties, taking the example of the action of modulatory substances onto the membrane leakage current, can provide the basis for the perceptual functions reflected in the macroscopic behavior of a cortical network. In the first part, the action of the modulatory system on cortical dynamics is investigated. First, it is demonstrated that the inclusion of these biophysical properties in a model of the primary visual cortex leads to the dynamic formation of synchronously active neuronal assemblies reflecting a context-dependent binding and segmentation of image components. Second, it is shown that the differential regulation of the leakage current can be used to bias the interactions of multiple cortical modules. This allows the flexible use of different feature domains for scene segmentation. Third, we demonstrate how, within the proposed architecture, the mapping of a moving stimulus onto the spatial dimension of the network results in an increased speed of synchronization. In the second part, we demonstrate how the differential regulation of neuromodulatory activity can be achieved in a self-consistent system. Three different mechanisms are described and investigated. This study thus demonstrates how a modulatory system, affecting the biophysical properties of single cells, can be used to achieve context-dependent processing at the system level.