Domains, Motifs, and Scaffolds: The Role of Modular Interactions in the Evolution and Wiring of Cell Signaling Circuits

Harnessing nature's toolbox: regulatory elements for synthetic biology

Synthetic biologists seek to engineer complex biological systems composed of modular elements. Achieving higher complexity in engineered biological organisms will require manipulating numerous systems of biological regulation: transcription; RNA interactions; protein signalling; and metabolic fluxes, among others. Exploiting the natural modularity at each level of biological regulation will promote the development of standardized tools for designing biological systems.

Protein evolution on a human signaling network

Background

The architectural structure of cellular networks provides a framework for innovations as well as constraints for protein evolution. This issue has previously been studied extensively by analyzing protein interaction networks. However, it is unclear how signaling networks influence and constrain protein evolution and conversely, how protein evolution modifies and shapes the functional consequences of signaling networks. In this study, we constructed a human signaling network containing more than 1,600 nodes and 5,000 links through manual curation of signaling pathways, and analyzed the d_N/d_S values of human-mouse orthologues on the network.

Results

We revealed that the protein d_N/d_S value decreases along the signal information flow from the extracellular space to nucleus. In the network, neighbor proteins tend to have similar d_N/d_S ratios, indicating neighbor proteins have similar evolutionary rates: co-fast or co-slow. However, different types of relationships (activating, inhibitory and neutral) between proteins have different effects on protein evolutionary rates, i.e., physically interacting protein pairs have the closest evolutionary rates. Furthermore, for directed shortest paths, the more distant two proteins are, the less chance they share similar evolutionary rates. However, such behavior was not observed for neutral shortest paths. Fast evolving signaling proteins have two modes of evolution: immunological proteins evolve more independently, while apoptotic proteins tend to form network components with other signaling proteins and share more similar evolutionary rates, possibly enhancing rapid information exchange between apoptotic and other signaling pathways.

Conclusion

Major network constraints on protein evolution in protein interaction networks previously described have been found for signaling networks. We further uncovered how network characteristics affect the evolutionary and co-evolutionary behavior of proteins and how protein evolution can modify the existing functionalities of signaling networks. These new insights provide some general principles for understanding protein evolution in the context of signaling networks.

The rules of disorder or why disorder rules

The finding that a large fraction of proteins (over 30%) in eukaryotic cells lack a unique three-dimensional structure but are functional has forced the scientific community to review its understanding of the structure-function paradigm. The involvement of many of these intrinsically unstructured proteins (IUPs) in intracellular signalling and regulatory processes as well as their central positioning (as interaction hubs) in recently mapped protein interaction networks is particularly intriguing. Here, we review the functional and structural properties of IUPs such as (i) their facilitated regulation via diverse post-translational modifications of specific amino acids (ii) scaffolding and recruitment of different binding partners in space and time via the "fly-casting" mechanism, through peptide motifs and by coupling folding with binding and (iii) conformational variability and adaptability. All of these properties allow these proteins to hold key positions in cellular organisation and regulation which in turn make them tractable as drug targets. In addition, we discuss how such properties, individually and in combination, facilitate combinatorial regulation and re-use of the same component in multiple biological processes.

From unicellular properties to multicellular behavior: bacteria quorum sensing circuitry and applications

Cell-cell communication and coordinated population-based behavior among single cell organisms have gained considerable attention in the recent years. The ability to send, receive, and process information allows unicellular organisms to act as multicellular entities and increases their chances of survival in complex environments. Quorum sensing (QS), a density-dependent cell-signaling mechanism, is one way by which bacteria 'talk' to one another. QS is commonly associated with adverse health effects such as biofilm formation, bacteria pathogenicity, and virulence. But there exists great potential to harness QS circuitry and its properties for other applications, enabling even wider societal impact. Interesting avenues are envisioned for the detection of chemicals and pathogens, the navigation of interspecies communication, the synchronization and control of cell phenotype, and the creation of novel materials based on synthetic biology. In this review, we first highlight the recent discoveries of the molecular underpinnings of QS function, with emphasis on the formation of biofilms. We then discuss how researchers have used QS circuitry to their advantage to build synthetic networks, rewire native metabolic pathways, and engineer cells for a variety of applications.

Nucleus-specific and cell cycle-regulated degradation of mitogen-activated protein kinase scaffold protein Ste5 contributes to the control of signaling competence

Saccharomyces cerevisiae cells are capable of responding to mating pheromone only prior to their exit from the G(1) phase of the cell cycle. Ste5 scaffold protein is essential for pheromone response because it couples pheromone receptor stimulation to activation of the appropriate mitogen-activated protein kinase (MAPK) cascade. In naïve cells, Ste5 resides primarily in the nucleus. Upon pheromone treatment, Ste5 is rapidly exported from the nucleus and accumulates at the tip of the mating projection via its association with multiple plasma membrane-localized molecules. We found that concomitant with its nuclear export, the rate of Ste5 turnover is markedly reduced. Preventing nuclear export destabilized Ste5, whereas preventing nuclear entry stabilized Ste5, indicating that Ste5 degradation occurs mainly in the nucleus. This degradation is dependent on ubiquitin and the proteasome. We show that Ste5 ubiquitinylation is mediated by the SCF(Cdc4) ubiquitin ligase and requires phosphorylation by the G(1) cyclin-dependent protein kinase (cdk1). The inability to efficiently degrade Ste5 resulted in pathway activation and cell cycle arrest in the absence of pheromone. These findings reveal that maintenance of this MAPK scaffold at an appropriately low level depends on its compartment-specific and cell cycle-dependent degradation. Overall, this mechanism provides a novel means for helping to prevent inadvertent stimulus-independent activation of a response and for restricting and maximizing the signaling competence of the cell to a specific cell cycle stage, which likely works hand in hand with the demonstrated role that G(1) Cdk1-dependent phosphorylation of Ste5 has in preventing its association with the plasma membrane.

An integrative approach for predicting interactions of protein regions

MOTIVATION

Protein-protein interactions are commonly mediated by the physical contact of distinct protein regions. Computational identification of interacting protein regions aids in the detailed understanding of protein networks and supports the prediction of novel protein interactions and the reconstruction of protein complexes.

RESULTS

We introduce an integrative approach for predicting protein region interactions using a probabilistic model fitted to an observed protein network. In particular, we consider globular domains, short linear motifs and coiled-coil regions as potential protein-binding regions. Possible cooperations between multiple regions within the same protein are taken into account. A.negrained confidence system allows for varying the impact of specific protein interactions and region annotations on the modeling process. We apply our prediction approach to a large training set using a maximum likelihood method, compare different scoring functions for region interactions and validate the predicted interactions against a collection of experimentally observed interactions. In addition, we analyze prediction performance with respect to the inclusion of different region types, the incorporation of confidence values for training data and the utilization of predicted protein interactions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular mechanism of the Syk activation switch

Many immune signaling pathways require activation of the Syk tyrosine kinase to link ligation of surface receptors to changes in gene expression. Despite the central role of Syk in these pathways, the Syk activation process remains poorly understood. In this work we quantitatively characterized the molecular mechanism of Syk activation in vitro using a real time fluorescence kinase assay, mutagenesis, and other biochemical techniques. We found that dephosphorylated full-length Syk demonstrates a low initial rate of substrate phosphorylation that increases during the kinase reaction due to autophosphorylation. The initial rate of Syk activity was strongly increased by either pre-autophosphorylation or binding of phosphorylated immune tyrosine activation motif peptides, and each of these factors independently fully activated Syk. Deletion mutagenesis was used to identify regions of Syk important for regulation, and residues 340-356 of the SH2 kinase linker region were identified to be important for suppression of activity before activation. Comparison of the activation processes of Syk and Zap-70 revealed that Syk is more readily activated by autophosphorylation than Zap-70, although both kinases are rapidly activated by Src family kinases. We also studied Syk activity in B cell lysates and found endogenous Syk is also activated by phosphorylation and immune tyrosine activation motif binding. Together these experiments show that Syk functions as an "OR-gate" type of molecular switch. This mechanism of switch-like activation helps explain how Syk is both rapidly activated after receptor binding but also sustains activity over time to facilitate longer term changes in gene expression.

Protein kinase A RII-like (R2D2) proteins exhibit differential localization and AKAP interaction

A-kinase anchoring proteins (AKAPs) bind to protein kinase A (PKA) via an amphipathic helix domain that interacts with a dimerization/docking domain on the regulatory (R) subunit of PKA. Four other mammalian proteins (ROPN1, ASP, SP17, and CABYR) also contain a highly conserved RII dimerization/docking (R2D2) domain, suggesting all four proteins may interact with all AKAPs in a manner similar to RII. All four of these proteins were originally detected in the flagellum of mammalian sperm. In this report, we demonstrate that all four R2D2 proteins are expressed in a wide variety of tissues and three of the proteins SP17, CABYR, and ASP are located in motile cilia of human bronchus and fallopian tubes. In addition, we detect SP17 in primary cilia. We also provide evidence that ROPN1 and ASP bind to a variety of AKAPs and this interaction can be disrupted with anchoring inhibitor peptides. The interaction of SP17 and CABYR with AKAPs appears to be much more limited. None of the R2D2 proteins appears to bind cAMP, a fundamental characteristic of the regulatory subunits of PKA. These observations suggest that R2D2 proteins utilize docking interactions with AKAPs to accomplish their function of regulating cilia and flagella. Based on location, affinity for AKAPs and lack of affinity for cAMP, it appears that each R2D2 protein has a unique role in this process.

Discordant on/off switching of gene expression in myocytes during cardiac hypertrophy in vivo

To determine whether the expression of cardiac genes changes in a graded manner or by on/off switching when cardiac myocytes change genetic programs in living animals, we have studied two indicator genes that change their expression oppositely in mouse binucleate ventricular cardiomyocytes during development and in response to cardiac hypertrophy. One is a single-copy transgene controlled by an alpha-myosin heavy chain (aMHC) promoter and coding for CFP. The other is the endogenous beta-myosin heavy chain (bMHC) gene modified to code for a YFP-bMHC fusion protein. Using high-resolution confocal microscopy, we determined the expression of the two indicator genes in individual cardiomyocytes perinatally and after inducing cardiac hypertrophy by transverse aortic constriction. Our results provide strong evidence that the cardiac genes respond by switching their expression in an on/off rather than graded manner, and that responding genes within a single cell and within the two nuclei of cardiomyocytes do not necessarily switch concordantly.

Complexity and modularity of intracellular networks: a systematic approach for modelling and simulation

Assembly of quantitative models of large complex networks brings about several challenges. One of them is the combinatorial complexity, where relatively few signalling molecules can combine to form thousands or millions of distinct chemical species. A receptor that has several separate phosphorylation sites can exist in hundreds of different states, many of which must be accounted for individually when simulating the time course of signalling. When assembly of protein complexes is being included, the number of distinct molecular species can easily increase by a few orders of magnitude. Validation, visualisation and understanding the network can become intractable. Another challenge appears when the modeller needs to recast or grow a model. Keeping track of changes and adding new elements present a significant difficulty. An approach to solve these challenges within the virtual cell (VCell) is described. Using (i) automatic extraction from pathway databases of model components (http://vcell.org/biopax) and (ii) rules of interactions that serve as reaction network generators (http://vcell.org/bionetgen), a way is provided for semi-automatic generation of quantitative mathematical models that also facilitates the reuse of model elements. In this approach, kinetic models of large, complex networks can be assembled from separately constructed modules, either directly or via rules. To implement this approach, the strength of several related technologies is combined: the BioPAX ontology, the BioNetGen rule-based description of molecular interactions and the VCell modelling and simulation framework.

System-level mapping of Escherichia coli response regulator dimerization with FRET hybrids

SUMMARY

Two-component signal transduction, featuring highly conserved histidine kinases (HKs) and response regulators (RRs), is one of the most prevalent signalling schemes in prokaryotes. RRs function as phosphorylation-activated switches to mediate diverse output responses, mostly via transcription regulation. As bacterial genomes typically encode multiple two-component proteins for distinct signalling pathways, the sequence and structural similarities of RR receiver domains create significant challenges to maintain interaction specificity. It is especially demanding for members of the OmpR/PhoB subfamily, the largest RR subfamily, which share a conserved dimerization interface for phosphorylation-mediated transcription regulation. We developed a strategy to investigate RR interaction by analysing Förster resonance energy transfer (FRET) between cyan fluorescent protein (CFP)- and yellow fluorescent protein (YFP)-fused RRs in vitro. Using the Escherichia coli RR PhoB as a model system, we were able to observe phosphorylation-dependent FRET between fluorescent protein (FP)-PhoB proteins and validated the FRET method by determining dimerization affinity and dimerization-coupled phosphorylation kinetics that recapitulated values determined by alternative methods. Further application of the FRET method to all E. coli OmpR/PhoB subfamily RRs revealed that phosphorylation-activated RR interaction is indeed a common theme for OmpR/PhoB subfamily RRs and these RRs display significant interaction specificity. Weak hetero-pair interactions were also identified between several different RRs, suggesting potential cross-regulation between distinct pathways.

Light-activated DNA binding in a designed allosteric protein

An understanding of how allostery, the conformational coupling of distant functional sites, arises in highly evolvable systems is of considerable interest in areas ranging from cell biology to protein design and signaling networks. We reasoned that the rigidity and defined geometry of an alpha-helical domain linker would make it effective as a conduit for allosteric signals. To test this idea, we rationally designed 12 fusions between the naturally photoactive LOV2 domain from Avena sativa phototropin 1 and the Escherichia coli trp repressor. When illuminated, one of the fusions selectively binds operator DNA and protects it from nuclease digestion. The ready success of our rational design strategy suggests that the helical "allosteric lever arm" is a general scheme for coupling the function of two proteins.

Estimation of the available free energy in a LOV2-J alpha photoswitch

Protein photosensors are versatile tools for studying ligand-regulated allostery and signaling. Fundamental to these processes is the amount of energy that can be provided by a photosensor to control downstream signaling events. Such regulation is exemplified by the phototropins--plant serine/threonine kinases that are activated by blue light via conserved LOV (light, oxygen and voltage) domains. The core photosensor of oat phototropin 1 is a LOV domain that interacts in a light-dependent fashion with an adjacent alpha-helix (J alpha) to control kinase activity. We used solution NMR measurements to quantify the free energy of the LOV domain-J alpha-helix binding equilibrium in the dark and lit states. These data indicate that light shifts this equilibrium by approximately 3.8 kcal mol(-1), thus quantifying the energy available through LOV-J alpha for light-driven allosteric regulation. This study provides insight into the energetics of light sensing by phototropins and benchmark values for engineering photoswitchable systems based on the LOV-J alpha interaction.

Role of the C-terminal SH3 domain and N-terminal tyrosine phosphorylation in regulation of Tim and related Dbl-family proteins

Dbl-related oncoproteins are guanine nucleotide exchange factors (GEFs) specific for Rho-family GTPases and typically possess tandem Dbl (DH) and pleckstrin homology (PH) domains that act in concert to catalyze exchange. Although the exchange potential of many Dbl-family proteins is constitutively activated by truncation, the precise mechanisms of regulation for many Dbl-family proteins are unknown. Tim and Vav are distantly related Dbl-family proteins that are similarly regulated; their Dbl homology (DH) domains interact with N-terminal helices to exclude and prevent activation of Rho GTPases. Phosphorylation, substitution, or deletion of the blocking helices relieves this autoinhibition. Here we show that two other Dbl-family proteins, Ngef and Wgef, which like Tim contain a C-terminal SH3 domain, are also activated by tyrosine phosphorylation of a blocking helix. Consequently, basal autoinhibition of DH domains by direct steric exclusion using short N-terminal helices likely represents a conserved mechanism of regulation for the large family of Dbl-related proteins. N-Terminal truncation or phosphorylation of many other Dbl-family GEFs leads to their activation; similar autoinhibition mechanisms could explain some of these events. In addition, we show that the C-terminal SH3 domain binding to a polyproline region N-terminal to the DH domain of the Tim subgroup of Dbl-family proteins provides a unique mechanism of regulated autoinhibition of exchange activity that is functionally linked to the interactions between the autoinhibitory helix and the DH domain.

Gene expansion and retention leads to a diverse tyrosine kinase superfamily in amphioxus

Tyrosine kinase (TK) proteins play a central role in cellular behavior and development of animals. The expansion of this superfamily is regarded as a key event in the evolution of the complex signaling pathways and gene networks of metazoans and is a prominent example of how shuffling of protein modules may generate molecular novelties. Using the intron/exon structure within the TK domain (TK intron code) as a complementary tool for the assignment of orthology and paralogy, we identified and studied the 118 TK proteins of the amphioxus Branchiostoma floridae genome to elucidate TK gene family evolution in metazoans and chordates in particular. Unlike all characterized metazoans to date, amphioxus has members of all known widespread TK families, with not a single loss. Putting amphioxus TKs in an evolutionary context, including new data from the cnidarian Nematostella vectensis, the echinoderm Strongylocentrotus purpuratus, and the ascidian Ciona intestinalis, we suggest new evolutionary histories for different TK families and draw a new global picture of gene loss/gain in the different phyla. Surprisingly, our survey also detected an unprecedented expansion of a group of closely related TK families, including TIE, FGFR, PDGFR, and RET, due most probably to massive gene duplication and exon shuffling. Based on their highly similar intron/exon structure at the TK domain, we suggest that this group of TK families constitute a superfamily of TK proteins, which we termed EXpanding TK, after their seemingly unique propensity to gene duplication and exon shuffling, not only in amphioxus but also across all metazoan groups. Due to this extreme tendency to both retention and expansion of TK genes, amphioxus harbors the richest and most diverse TK repertoire among all metazoans studied so far, retaining most of the gene complement of its ancestors, but having evolved its own repertoire of genetic novelties.

Specificity of the chromodomain Y chromosome family of chromodomains for lysine-methylated ARK(S/T) motifs

Previous studies have shown two homologous chromodomain modules in the HP1 and Polycomb proteins exhibit discriminatory binding to related methyllysine residues (embedded in ARKS motifs) of the histone H3 tail. Methylated ARK(S/T) motifs have recently been identified in other chromatin factors (e.g. linker histone H1.4 and lysine methyltransferase G9a). These are thought to function as peripheral docking sites for the HP1 chromodomain. In vertebrates, HP1-like chromodomains are also present in the chromodomain Y chromosome (CDY) family of proteins adjacent to a putative catalytic motif. The human genome encodes three CDY family proteins, CDY, CDYL, and CDYL2. These have putative functions ranging from establishment of histone H4 acetylation during spermiogenesis to regulation of transcription co-repressor complexes. To delineate the biochemical functions of the CDY family chromodomains, we analyzed their specificity of methyllysine recognition. We detected substantial differences among these factors. The CDY chromodomain exhibits discriminatory binding to lysine-methylated ARK(S/T) motifs, whereas the CDYL2 chromodomain binds with comparable strength to multiple ARK(S/T) motifs. Interestingly, subtle amino acid changes in the CDYL chromodomain prohibit such binding interactions in vitro and in vivo. However, point mutations can rescue binding. In support of the in vitro binding properties of the chromodomains, the full-length CDY family proteins exhibit substantial variability in chromatin localization. Our studies underscore the significance of subtle sequence differences in a conserved signaling module for diverse epigenetic regulatory pathways.

The Ser(186) phospho-acceptor site within ERK4 is essential for its ability to interact with and activate PRAK/MK5

ERK (extracellular-signal-regulated kinase) 4 [MAPK (mitogen-activated protein kinase) 4] and ERK3 (MAPK6) are atypical MAPKs. One major difference between these proteins and the classical MAPKs is substitution of the conserved T-X-Y motif within the activation loop by a single phospho-acceptor site within an S-E-G motif. In the present study we report that Ser(186) of the S-E-G motif in ERK4 is phosphorylated in vivo. Kinase-dead ERK4 is also phosphorylated on Ser(186), indicating that an ERK4 kinase, rather than autophosphorylation, is responsible. Co-expression of MK5 [MAPK-activated protein kinase 5; also known as PRAK (p38-regulated/activated kinase)], a physiological target of ERK4, increases phosphorylation of Ser(186). This is not dependent on MK5 activity, but does require interaction between ERK4 and MK5 suggesting that MK5 binding either prevents ERK4 dephosphorylation or facilitates ERK4 kinase activity. ERK4 mutants in which Ser(186) is replaced with either an alanine residue or a phospho-mimetic residue (glutamate) are unable to activate MK5 and Ser(186) is also required for cytoplasmic anchoring of MK5. Both defects seem to reflect an impaired ability of the ERK4 mutants to interact with MK5. We find that there are at least two endogenous pools of wild-type ERK4. One form exhibits reduced mobility when analysed using SDS/PAGE. This is due to MK5-dependent phosphorylation and only this retarded ERK4 species is both phosphorylated on Ser(186) and co-immunoprecipitates with wild-type MK5. We conclude that binding between ERK4 and MK5 facilitates phosphorylation of Ser(186) and stabilization of the ERK4-MK5 complex. This results in phosphorylation and activation of MK5, which in turn phosphorylates ERK4 on sites other than Ser(186) resulting in the observed mobility shift.

Proteomic approaches to the analysis of multiprotein signaling complexes

Signal transduction is one of the most active fields in modern biomedical research. Increasing evidence has shown that signaling proteins associate with each other in characteristic ways to form large signaling complexes. These diverse structures operate to boost signaling efficiency, ensure specificity and increase sensitivity of the biochemical circuitry. Traditional methods of protein analysis are inadequate to fully characterize and understand these structures, which are intricate, contain many components and are highly dynamic. Instead, proteomics technologies are currently being applied to investigate the nature and composition of multimeric signaling complexes. This review presents commonly used and potential proteomic methods of analyzing diverse protein complexes along with a discussion and a brief evaluation of alternative approaches. Challenges associated with proteomic analysis of signaling complexes are also discussed.

Biosynthetic potentials of metabolites and their hierarchical organization

A major challenge in systems biology is to understand how complex and highly connected metabolic networks are organized. The structure of these networks is investigated here by identifying sets of metabolites that have a similar biosynthetic potential. We measure the biosynthetic potential of a particular compound by determining all metabolites than can be produced from it and, following a terminology introduced previously, call this set the scope of the compound. To identify groups of compounds with similar scopes, we apply a hierarchical clustering method. We find that compounds within the same cluster often display similar chemical structures and appear in the same metabolic pathway. For each cluster we define a consensus scope by determining a set of metabolites that is most similar to all scopes within the cluster. This allows for a generalization from scopes of single compounds to scopes of a chemical family. We observe that most of the resulting consensus scopes overlap or are fully contained in others, revealing a hierarchical ordering of metabolites according to their biosynthetic potential. Our investigations show that this hierarchy is not only determined by the chemical complexity of the metabolites, but also strongly by their biological function. As a general tendency, metabolites which are necessary for essential cellular processes exhibit a larger biosynthetic potential than those involved in secondary metabolism. A central result is that chemically very similar substances with different biological functions may differ significantly in their biosynthetic potentials. Our studies provide an important step towards understanding fundamental design principles of metabolic networks determined by the structural and functional complexity of metabolites.

Using chemistry and microfluidics to understand the spatial dynamics of complex biological networks

Understanding the spatial dynamics of biochemical networks is both fundamentally important for understanding life at the systems level and also has practical implications for medicine, engineering, biology, and chemistry. Studies at the level of individual reactions provide essential information about the function, interactions, and localization of individual molecular species and reactions in a network. However, analyzing the spatial dynamics of complex biochemical networks at this level is difficult. Biochemical networks are nonequilibrium systems containing dozens to hundreds of reactions with nonlinear and time-dependent interactions, and these interactions are influenced by diffusion, flow, and the relative values of state-dependent kinetic parameters. To achieve an overall understanding of the spatial dynamics of a network and the global mechanisms that drive its function, networks must be analyzed as a whole, where all of the components and influential parameters of a network are simultaneously considered. Here, we describe chemical concepts and microfluidic tools developed for network-level investigations of the spatial dynamics of these networks. Modular approaches can be used to simplify these networks by separating them into modules, and simple experimental or computational models can be created by replacing each module with a single reaction. Microfluidics can be used to implement these models as well as to analyze and perturb the complex network itself with spatial control on the micrometer scale. We also describe the application of these network-level approaches to elucidate the mechanisms governing the spatial dynamics of two networkshemostasis (blood clotting) and early patterning of the Drosophila embryo. To investigate the dynamics of the complex network of hemostasis, we simplified the network by using a modular mechanism and created a chemical model based on this mechanism by using microfluidics. Then, we used the mechanism and the model to predict the dynamics of initiation and propagation of blood clotting and tested these predictions with human blood plasma by using microfluidics. We discovered that both initiation and propagation of clotting are regulated by a threshold response to the concentration of activators of clotting, and that clotting is sensitive to the spatial localization of stimuli. To understand the dynamics of patterning of the Drosophila embryo, we used microfluidics to perturb the environment around a developing embryo and observe the effects of this perturbation on the expression of Hunchback, a protein whose localization is essential to proper development. We found that the mechanism that is responsible for Hunchback positioning is asymmetric, time-dependent, and more complex than previously proposed by studies of individual reactions. Overall, these approaches provide strategies for simplifying, modeling, and probing complex networks without sacrificing the functionality of the network. Such network-level strategies may be most useful for understanding systems with nonlinear interactions where spatial dynamics is essential for function. In addition, microfluidics provides an opportunity to investigate the mechanisms responsible for robust functioning of complex networks. By creating nonideal, stressful, and perturbed environments, microfluidic experiments could reveal the function of pathways thought to be nonessential under ideal conditions.

Using engineered scaffold interactions to reshape MAP kinase pathway signaling dynamics

Scaffold proteins link signaling molecules into linear pathways by physically assembling them into complexes. Scaffolds may also have a higher-order role as signal-processing hubs, serving as the target of feedback loops that optimize signaling amplitude and timing. We demonstrate that the Ste5 scaffold protein can be used as a platform to systematically reshape output of the yeast mating MAP kinase pathway. We constructed synthetic positive- and negative-feedback loops by dynamically regulating recruitment of pathway modulators to an artificial binding site on Ste5. These engineered circuits yielded diverse behaviors: ultrasensitive dose response, accelerated or delayed response times, and tunable adaptation. Protein scaffolds provide a flexible platform for reprogramming cellular responses and could be exploited to engineer cells with novel therapeutic and biotechnological functions.

Membrane recognition by phospholipid-binding domains

Many different globular domains bind to the surfaces of cellular membranes, or to specific phospholipid components in these membranes, and this binding is often tightly regulated. Examples include pleckstrin homology and C2 domains, which are among the largest domain families in the human proteome. Crystal structures, binding studies and analyses of subcellular localization have provided much insight into how members of this diverse group of domains bind to membranes, what features they recognize and how binding is controlled. A full appreciation of these processes is crucial for understanding how protein localization and membrane topography and trafficking are regulated in cells.

Computational models of tandem SRC homology 2 domain interactions and application to phosphoinositide 3-kinase

Intracellular signal transduction proteins typically utilize multiple interaction domains for proper targeting, and thus a broad diversity of distinct signaling complexes may be assembled. Considering the coordination of only two such domains, as in tandem Src homology 2 (SH2) domain constructs, gives rise to a kinetic scheme that is not adequately described by simple models used routinely to interpret in vitro binding measurements. To analyze the interactions between tandem SH2 domains and bisphosphorylated peptides, we formulated detailed kinetic models and applied them to the phosphoinositide 3-kinase p85 regulatory subunit/platelet-derived growth factor beta-receptor system. Data for this system from different in vitro assay platforms, including surface plasmon resonance, competition binding, and isothermal titration calorimetry, were reconciled to estimate the magnitude of the cooperativity characterizing the sequential binding of the high and low affinity SH2 domains (C-SH2 and N-SH2, respectively). Compared with values based on an effective volume approximation, the estimated cooperativity is 3 orders of magnitude lower, indicative of significant structural constraints. Homodimerization of full-length p85 was found to be an alternative mechanism for high avidity binding to phosphorylated platelet-derived growth factor receptors, which would render the N-SH2 domain dispensable for receptor binding.

X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution

Crystallography supplies unparalleled detail on structural information critical for mechanistic analyses; however, it is restricted to describing low energy conformations of macromolecules within crystal lattices. Small angle X-ray scattering (SAXS) offers complementary information about macromolecular folding, unfolding, aggregation, extended conformations, flexibly linked domains, shape, conformation, and assembly state in solution, albeit at the lower resolution range of about 50 A to 10 A resolution, but without the size limitations inherent in NMR and electron microscopy studies. Together these techniques can allow multi-scale modeling to create complete and accurate images of macromolecules for modeling allosteric mechanisms, supramolecular complexes, and dynamic molecular machines acting in diverse processes ranging from eukaryotic DNA replication, recombination and repair to microbial membrane secretion and assembly systems. This review addresses both theoretical and practical concepts, concerns and considerations for using these techniques in conjunction with computational methods to productively combine solution scattering data with high-resolution structures. Detailed aspects of SAXS experimental results are considered with a focus on data interpretation tools suitable to model protein and nucleic acid macromolecular structures, including membrane protein, RNA, DNA, and protein-nucleic acid complexes. The methods discussed provide the basis to examine molecular interactions in solution and to study macromolecular flexibility and conformational changes that have become increasingly relevant for accurate understanding, simulation, and prediction of mechanisms in structural cell biology and nanotechnology.

In situ observation of protein phosphorylation by high-resolution NMR spectroscopy

Although the biological significance of protein phosphorylation in cellular signaling is widely appreciated, methods to directly detect these post-translational modifications in situ are lacking. Here we introduce the application of high-resolution NMR spectroscopy for observing de novo protein phosphorylation in vitro and in Xenopus laevis egg extracts and whole live oocyte cells. We found that the stepwise modification of adjacent casein kinase 2 (CK2) substrate sites within the viral SV40 large T antigen regulatory region proceeded in a defined order and through intermediate substrate release. This kinase mechanism contrasts with a more intuitive mode of CK2 action in which the kinase would remain substrate bound to perform both modification reactions without intermediate substrate release. For cellular signaling pathways, the transient availability of partially modified CK2 substrates could exert important switch-like regulatory functions.

Dynamic scaffolding in a G protein-coupled signaling system

The INAD scaffold organizes a multiprotein complex that is essential for proper visual signaling in Drosophila photoreceptor cells. Here we show that one of the INAD PDZ domains (PDZ5) exists in a redox-dependent equilibrium between two conformations--a reduced form that is similar to the structure of other PDZ domains, and an oxidized form in which the ligand-binding site is distorted through formation of a strong intramolecular disulfide bond. We demonstrate transient light-dependent formation of this disulfide bond in vivo and find that transgenic flies expressing a mutant INAD in which PDZ5 is locked in the reduced state display severe defects in termination of visual responses and visually mediated reflex behavior. These studies demonstrate a conformational switch mechanism for PDZ domain function and suggest that INAD behaves more like a dynamic machine rather than a passive scaffold, regulating signal transduction at the millisecond timescale through cycles of conformational change.

Characterization of the 4D5Flu single-chain antibody with a stimulus-responsive elastin-like peptide linker: a potential reporter of peptide linker conformation

Single-chain antibodies (scFvs) are comprised of IgG variable light and variable heavy domains tethered together by a peptide linker whose length and sequence can affect antigen binding properties. The ability to modulate antigen binding affinity through the use of environmental triggers would be of great interest for many biotechnological applications. We have characterized the antigen binding properties of an anti-fluorescein scFv, 4D5Flu, containing stimulus-responsive short elastin-like peptide linkers and nonresponsive flexible linkers. Comparison of length-matched flexible and short elastin-like peptide linkers indicates that a stimulus-responsive linker can confer stimulus-responsive control of fluorescein binding. A linker length of either six or 10 amino acids proved to have the largest thermally induced response. Similar differences in binding free energy changes indicate a common underlying mechanism of thermal responsiveness. Contrary to the thermal behavior, the effect of salt, another elastin beta-turn-inducing stimulus, stabilized antigen binding in the six- and 10-amino-acid linkers such that elastin-like linkers became less stimulus-responsive as compared with flexible linkers. Again, the thermodynamic analysis indicates a common mechanism of salt responsiveness. Characterization of the room-temperature binding affinities and evidence indicating a dimeric state of the scFvs concomitantly suggest the major contribution to the stimulus-responsive behavior derives from the perturbation of interdomain associations, rather than the linker-constrained disruption of the intramolecular association. The ability to use stimulus-responsive peptide modules to exert a novel control over protein function will likely find application in the creation of allosteric antibodies and scFv-based biosensors, and as a platform to enable the evolution of new stimulus-responsive peptides.

Building a Parallel Metabolism within the Cell

One of the key aims of synthetic biology is to engineer artificial processes inside living cells. This requires components that interact in a predictable manner, both with each other and with existing cellular systems. However, the activity of many components is constrained by their interactions with other cellular molecules and often their roles in maintaining cell health. To escape this limitation, researchers are pursuing an "orthogonal" approach, building a parallel metabolism within the cell. Components of this parallel metabolism can be sourced from evolutionarily distant species or reengineered from existing cellular molecules by using rational design and directed evolution. These approaches allow the study of basic principles in cell biology and the engineering of cells that can function as environmental sensors, simple computers, and drug factories.

Domains mediate protein-protein interactions and nucleate protein assemblies

Cell physiology is governed by an intricate mesh of physical and functional links among proteins, nucleic acids and other metabolites. The recent information flood coming from large-scale genomic and proteomic approaches allows us to foresee the possibility of compiling an exhaustive list of the molecules present within a cell, enriched with quantitative information on concentration and cellular localization. Moreover, several high-throughput experimental and computational techniques have been devised to map all the protein interactions occurring in a living cell. So far, such maps have been drawn as graphs where nodes represent proteins and edges represent interactions. However, this representation does not take into account the intrinsically modular nature of proteins and thus fails in providing an effective description of the determinants of binding. Since proteins are composed of domains that often confer on proteins their binding capabilities, a more informative description of the interaction network would detail, for each pair of interacting proteins in the network, which domains mediate the binding. Understanding how protein domains combine to mediate protein interactions would allow one to add important features to the protein interaction network, making it possible to discriminate between simultaneously occurring and mutually exclusive interactions. This objective can be achieved by experimentally characterizing domain recognition specificity or by analyzing the frequency of co-occurring domains in proteins that do interact. Such approaches allow gaining insights on the topology of complexes with unknown three-dimensional structure, thus opening the prospect of adopting a more rational strategy in developing drugs designed to selectively target specific protein interactions.

RPEL motifs link the serum response factor cofactor MAL but not myocardin to Rho signaling via actin binding

Myocardin (MC) family proteins are transcriptional coactivators for serum response factor (SRF). Each family member possesses a conserved N-terminal region containing three RPEL motifs (the "RPEL domain"). MAL/MKL1/myocardin-related transcription factor A is cytoplasmic, accumulating in the nucleus upon activation of Rho GTPase signaling, which alters interactions between G-actin and the RPEL domain. We demonstrate that MC, which is nuclear, does not shuttle through the cytoplasm and that the contrasting nucleocytoplasmic shuttling properties of MAL and MC are defined by their RPEL domains. We show that the MAL RPEL domain binds actin more avidly than that of MC and that the RPEL motif itself is an actin-binding element. RPEL1 and RPEL2 of MC bind actin weakly compared with those of MAL, while RPEL3 is of comparable and low affinity in the two proteins. Actin binding by all three motifs is required for MAL regulation. The differing behaviors of MAL and MC are specified by the RPEL1-RPEL2 unit, while RPEL3 can be exchanged between them. We propose that differential actin occupancy of multiple RPEL motifs regulates nucleocytoplasmic transport and activity of MAL.

WNK1 is a novel regulator of Munc18c-syntaxin 4 complex formation in soluble NSF attachment protein receptor (SNARE)-mediated vesicle exocytosis

Defects in soluble NSF attachment protein receptor (SNARE)-mediated granule exocytosis occur in islet beta cells, adipocytes, and/or skeletal muscle cells correlate with increased susceptibility to insulin resistance and diabetes. The serine/threonine kinase WNK1 (with no K (lysine)) has recently been implicated in exocytosis and is expressed in all three of these cell types. To search for WNK1 substrates related to exocytosis, we conducted a WNK1 two-hybrid screen, which yielded Munc18c. Munc18c is known to be a key regulator of accessibility of the target membrane (t-SNARE) protein syntaxin 4 to participate in SNARE core complex assembly, although a paucity of Munc18c-binding factors has precluded discovery of its precise functions. To validate WNK1 as a new Munc18c-interacting partner, the direct interaction between WNK1 and Munc18c was confirmed using in vitro binding analysis, and endogenous WNK1-Munc18c complexes were detected in the cytosolic and plasma membrane compartments of the islet beta cell line MIN6. This binding interaction is mediated through the N-terminal 172 residues of Munc18c and the kinase domain residues of WNK1 (residues 159-491). Expression of either of these two minimal interaction domains resulted in inhibition of glucose-stimulated insulin secretion, consistent with a functional importance for the endogenous WNK1-Munc18c complex in exocytosis. Interestingly, Munc18c failed to serve as a WNK1 substrate in kinase activity assays, suggesting that WNK1 functions in SNARE complex assembly outside its role as a kinase. Taken together, these data support a novel role for WNK1 and a new mechanism for the regulation of SNARE complex assembly by WNK1-Munc18c complexes.

Defining the specificity space of the human SRC homology 2 domain

Src homology 2 (SH2) domains are the largest family of interaction modules encoded by the human genome to recognize tyrosine-phosphorylated sequences and thereby play pivotal roles in transducing and controlling cellular signals emanating from protein-tyrosine kinases. Different SH2 domains select for distinct phosphopeptides, and the function of a given SH2 domain is often dictated by the specific motifs that it recognizes. Therefore, deciphering the phosphotyrosyl peptide motif recognized by an SH2 domain is the key to understanding its cellular function. Here we cloned all 120 SH2 domains identified in the human genome and determined the phosphotyrosyl peptide binding properties of 76 SH2 domains by screening an oriented peptide array library. Of these 76, we defined the selectivity for 43 SH2 domains and refined the binding motifs for another 33 SH2 domains. We identified a number of novel binding motifs, which are exemplified by the BRDG1 SH2 domain that selects specifically for a bulky, hydrophobic residue at P + 4 relative to the Tyr(P) residue. Based on the oriented peptide array library data, we developed scoring matrix-assisted ligand identification (or SMALI), a Web-based program for predicting binding partners for SH2-containing proteins. When applied to SH2D1A/SAP (SLAM-associated protein), a protein whose mutation or deletion underlies the X-linked lymphoproliferative syndrome, SMALI not only recapitulated known interactions but also identified a number of novel interacting proteins for this disease-associated protein. SMALI also identified a number of potential interactors for BRDG1, a protein whose function is largely unknown. Peptide in-solution binding analysis demonstrated that a SMALI score correlates well with the binding energy of a peptide to a given SH2 domain. The definition of the specificity space of the human SH2 domain provides both the necessary molecular basis and a platform for future exploration of the functions for SH2-containing proteins in cells.

Environmental variability and modularity of bacterial metabolic networks

BACKGROUND

Biological systems are often modular: they can be decomposed into nearly-independent structural units that perform specific functions. The evolutionary origin of modularity is a subject of much current interest. Recent theory suggests that modularity can be enhanced when the environment changes over time. However, this theory has not yet been tested using biological data.

RESULTS

To address this, we studied the relation between environmental variability and modularity in a natural and well-studied system, the metabolic networks of bacteria. We classified 117 bacterial species according to the degree of variability in their natural habitat. We find that metabolic networks of organisms in variable environments are significantly more modular than networks of organisms that evolved under more constant conditions.

CONCLUSION

This study supports the view that variability in the natural habitat of an organism promotes modularity in its metabolic network and perhaps in other biological systems.

Where have all the interactions gone? Estimating the coverage of two-hybrid protein interaction maps

Yeast two-hybrid screens are an important method for mapping pairwise physical interactions between proteins. The fraction of interactions detected in independent screens can be very small, and an outstanding challenge is to determine the reason for the low overlap. Low overlap can arise from either a high false-discovery rate (interaction sets have low overlap because each set is contaminated by a large number of stochastic false-positive interactions) or a high false-negative rate (interaction sets have low overlap because each misses many true interactions). We extend capture–recapture theory to provide the first unified model for false-positive and false-negative rates for two-hybrid screens. Analysis of yeast, worm, and fly data indicates that 25% to 45% of the reported interactions are likely false positives. Membrane proteins have higher false-discovery rates on average, and signal transduction proteins have lower rates. The overall false-negative rate ranges from 75% for worm to 90% for fly, which arises from a roughly 50% false-negative rate due to statistical undersampling and a 55% to 85% false-negative rate due to proteins that appear to be systematically lost from the assays. Finally, statistical model selection conclusively rejects the Erdös-Rényi network model in favor of the power law model for yeast and the truncated power law for worm and fly degree distributions. Much as genome sequencing coverage estimates were essential for planning the human genome sequencing project, the coverage estimates developed here will be valuable for guiding future proteomic screens. All software and datasets are available in Datasets S1 and S2, Figures S1–S5, and Tables S1−S6, and are also available from our Web site, http://www.baderzone.org.

The genome sequence of an organism provides a parts list of proteins, but not an instruction manual for assembling the parts into a cell. Assembly instructions now come from experiments such as two-hybrid screens that detect physical interactions between pairs of proteins. Defining the resources required for generating a full interaction map requires accurate estimates of the false-negative and false-positive rates of genome-scale screens. Two-hybrid screens often select a query protein and sample its interaction partners. True partners may be missed, and false partners may be spuriously identified. This sampling process resembles a capture–recapture experiment, except that classical capture–recapture theory assumes no false positives. Novel extensions to capture–recapture theory permit its application to proteomic screens. This new theory provides statistically grounded answers to long-standing questions: false-discovery rates of high-throughput screens (possibly over 50% per unique interaction, but probably no more than 15% per clone); the quality of different screening libraries; protein properties leading to “sticky” or “promiscuous” interactions; the global network topology; and, most importantly, the coverage of existing two-hybrid maps. Models estimate roughly 30,000 total pairwise interactions in yeast and 500,000 to 1,000,000 in metazoans. The majority of these interactions remain to be discovered.

Synthetic gene circuits: design with directed evolution

Synthetic circuits offer great promise for generating insights into nature's underlying design principles or forward engineering novel biotechnology applications. However, construction of these circuits is not straightforward. Synthetic circuits generally consist of components optimized to function in their natural context, not in the context of the synthetic circuit. Combining mathematical modeling with directed evolution offers one promising means for addressing this problem. Modeling identifies mutational targets and limits the evolutionary search space for directed evolution, which alters circuit performance without the need for detailed biophysical information. This review examines strategies for integrating modeling and directed evolution and discusses the utility and limitations of available methods.

Mechanisms of specificity in protein phosphorylation

A typical protein kinase must recognize between one and a few hundred bona fide phosphorylation sites in a background of approximately 700,000 potentially phosphorylatable residues. Multiple mechanisms have evolved that contribute to this exquisite specificity, including the structure of the catalytic site, local and distal interactions between the kinase and substrate, the formation of complexes with scaffolding and adaptor proteins that spatially regulate the kinase, systems-level competition between substrates, and error-correction mechanisms. The responsibility for the recognition of substrates by protein kinases appears to be distributed among a large number of independent, imperfect specificity mechanisms.

Function and regulation in MAPK signaling pathways: lessons learned from the yeast Saccharomyces cerevisiae

Signaling pathways that activate different mitogen-activated protein kinases (MAPKs) elicit many of the responses that are evoked in cells by changes in certain environmental conditions and upon exposure to a variety of hormonal and other stimuli. These pathways were first elucidated in the unicellular eukaryote Saccharomyces cerevisiae (budding yeast). Studies of MAPK pathways in this organism continue to be especially informative in revealing the molecular mechanisms by which MAPK cascades operate, propagate signals, modulate cellular processes, and are controlled by regulatory factors both internal to and external to the pathways. Here we highlight recent advances and new insights about MAPK-based signaling that have been made through studies in yeast, which provide lessons directly applicable to, and that enhance our understanding of, MAPK-mediated signaling in mammalian cells.

Antibody-induced oligomerization and activation of an engineered reporter enzyme

Our objective is to produce a protein biosensor (or molecular switch) that is specifically activated in solution by a monoclonal antibody. Many effector-dependent enzymes have evolved in nature, but the introduction of a novel regulatory mechanism into a normally unregulated enzyme poses a difficult design problem. We used site-saturation mutagenesis and screening to generate effector-activated variants of the reporter enzyme beta-glucuronidase (GUS). The specific activity of the purified epitope-tagged GUS variant was increased by up to approximately 500-fold by the addition of an equimolar concentration of a monoclonal antibody. This molecular switch is modular in design, so it can easily be re-engineered for the detection of other peptide-specific antibodies. Such antibody-activated reporters could someday enable point-of-care serological assays for the rapid detection of infectious diseases.

God as Genetic Engineer

The Edge of Evolution . The Search for the Limits of Darwinism. By Michael J. Behe . Free Press, New York, 2007. 331 pp. $28, C$33.99. ISBN 9780743296205.

To support belief in intelligent design, Behe argues that random mutation and natural selection cannot account for the important changes in the history of life on Earth.

Rewiring cellular morphology pathways with synthetic guanine nucleotide exchange factors

Eukaryotic cells mobilize the actin cytoskeleton to generate a remarkable diversity of morphological behaviours, including motility, phagocytosis and cytokinesis. Much of this diversity is mediated by guanine nucleotide exchange factors (GEFs) that activate Rho family GTPases-the master regulators of the actin cytoskeleton. There are over 80 Rho GEFs in the human genome (compared to only 22 genes for the Rho GTPases themselves), and the evolution of new and diverse GEFs is thought to provide a mechanism for linking the core cytoskeletal machinery to a wide range of new control inputs. Here we test this hypothesis and ask if we can systematically reprogramme cellular morphology by engineering synthetic GEF proteins. We focused on Dbl family Rho GEFs, which have a highly modular structure common to many signalling proteins: they contain a catalytic Dbl homology (DH) domain linked to diverse regulatory domains, many of which autoinhibit GEF activity. Here we show that by recombining catalytic GEF domains with new regulatory modules, we can generate synthetic GEFs that are activated by non-native inputs. We have used these synthetic GEFs to reprogramme cellular behaviour in diverse ways. The GEFs can be used to link specific cytoskeletal responses to normally unrelated upstream signalling pathways. In addition, multiple synthetic GEFs can be linked as components in series to form an artificial cascade with improved signal processing behaviour. These results show the high degree of evolutionary plasticity of this important family of modular signalling proteins, and indicate that it may be possible to use synthetic biology approaches to manipulate the complex spatio-temporal control of cell morphology.

Structure-based kinetic models of modular signaling protein function: focus on Shp2

We present here a computational, rule-based model to study the function of the SH2 domain-containing protein tyrosine phosphatase, Shp2, in intracellular signal transduction. The two SH2 domains of Shp2 differentially regulate the enzymatic activity by a well-characterized mechanism, but they also affect the targeting of Shp2 to signaling receptors in cells. Our kinetic model integrates these potentially competing effects by considering the intra- and intermolecular interactions of the Shp2 SH2 domains and catalytic site as well as the effect of Shp2 phosphorylation. Even for the isolated Shp2/receptor system, which may seem simple by certain standards, we find that the network of possible binding and phosphorylation states is composed of over 1000 members. To our knowledge, this is the first kinetic model to fully consider the modular, multifunctional structure of a signaling protein, and the computational approach should be generally applicable to other complex intermolecular interactions.

Designing biological systems

The design of artificial biological systems and the understanding of their natural counterparts are key objectives of the emerging discipline of synthetic biology. Toward both ends, research in synthetic biology has primarily focused on the construction of simple devices, such as transcription-based oscillators and switches. Construction of such devices should provide us with insight on the design of natural systems, indicating whether our understanding is complete or whether there are still gaps in our knowledge. Construction of simple biological systems may also lay the groundwork for the construction of more complex systems that have practical utility. To realize its full potential, biological systems design borrows from the allied fields of protein design and metabolic engineering. In this review, we describe the scientific accomplishments in this field, as well as its forays into biological part standardization and education of future biological designers.

Engineering modular protein interaction switches by sequence overlap

Many cellular signaling pathways contain proteins whose interactions change in response to upstream inputs, allowing for conditional activation or repression of the interaction based on the presence of the input molecule. The ability to engineer similar regulation into protein interaction elements would provide us with powerful tools for controlling cell signaling. Here we describe an approach for engineering diverse synthetic protein interaction switches. Specifically, by overlapping the sequences of pairs of protein interaction domains and peptides, we have been able to generate mutually exclusive regulation over their interactions. Thus, the hybrid protein (which is composed of the two overlapped interaction modules) can bind to either of the two respective ligands for those modules, but not to both simultaneously. We show that these synthetic switch proteins can be used to regulate specific protein-protein interactions in vivo. These switches allow us to disrupt an interaction with the addition or activation of a protein input that has no natural connection to the interaction in question. Therefore, they give us the ability to make novel connections between normally unrelated signaling pathways and to rewire the input/output relationships of cellular behaviors. Our experiments also suggest a possible mechanism by which complex regulatory proteins might have evolved from simpler components.