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

Evolutionary synthetic biology

Signaling networks process vast amounts of environmental information to generate specific cellular responses. As cellular environments change, signaling networks adapt accordingly. Here, I will discuss how the integration of synthetic biology and directed evolution approaches is shedding light on the molecular mechanisms that guide the evolution of signaling networks. In particular, I will review studies that demonstrate how different types of mutations, from the replacement of individual amino acids to the shuffling of modular domains, lead to markedly different evolutionary trajectories and consequently to diverse network rewiring. Moreover, I will argue that intrinsic evolutionary properties of signaling proteins, such as the robustness of wild type functions, the promiscuous nature of evolutionary intermediates, and the modular decoupling between binding and catalysis, play important roles in the evolution of signaling networks. Finally, I will argue that rapid advances in our ability to synthesize DNA will radically alter how we study signaling network evolution at the genome-wide level.

Tomato TFT1 is required for PAMP-triggered immunity and mutations that prevent T3S effector XopN from binding to TFT1 attenuate Xanthomonas virulence

XopN is a type III effector protein from Xanthomonas campestris pathovar vesicatoria that suppresses PAMP-triggered immunity (PTI) in tomato. Previous work reported that XopN interacts with the tomato 14-3-3 isoform TFT1; however, TFT1's role in PTI and/or XopN virulence was not determined. Here we show that TFT1 functions in PTI and is a XopN virulence target. Virus-induced gene silencing of TFT1 mRNA in tomato leaves resulted in increased growth of Xcv ΔxopN and Xcv ΔhrpF demonstrating that TFT1 is required to inhibit Xcv multiplication. TFT1 expression was required for Xcv-induced accumulation of PTI5, GRAS4, WRKY28, and LRR22 mRNAs, four PTI marker genes in tomato. Deletion analysis revealed that the XopN C-terminal domain (amino acids 344–733) is sufficient to bind TFT1. Removal of amino acids 605–733 disrupts XopN binding to TFT1 in plant extracts and inhibits XopN-dependent virulence in tomato, demonstrating that these residues are necessary for the XopN/TFT1 interaction. Phos-tag gel analysis and mass spectrometry showed that XopN is phosphorylated in plant extracts at serine 688 in a putative 14-3-3 recognition motif. Mutation of S688 reduced XopN's phosphorylation state but was not sufficient to inhibit binding to TFT1 or reduce XopN virulence. Mutation of S688 and two leucines (L64,L65) in XopN, however, eliminated XopN binding to TFT1 in plant extracts and XopN virulence. L64 and L65 are required for XopN to bind TARK1, a tomato atypical receptor kinase required for PTI. This suggested that TFT1 binding to XopN's C-terminal domain might be stabilized via TARK1/XopN interaction. Pull-down and BiFC analyses show that XopN promotes TARK1/TFT1 complex formation in vitro and in planta by functioning as a molecular scaffold. This is the first report showing that a type III effector targets a host 14-3-3 involved in PTI to promote bacterial pathogenesis.

Bacterial pathogens of plants and animals employ the type III secretion system to secrete and translocate effector proteins into host cells to suppress defense responses. Biochemical analyses have revealed that several effector proteins mimic host enzyme activities to directly interfere with pathogen perception, defense signal transduction, and/or secretion of antimicrobial compounds. The XopN effector from Xanthomonas campestris pathovar vesicatoria, the causal agent of bacterial spot on tomato and pepper plants, suppresses PAMP-triggered immunity (PTI) but structural modeling predicts it encodes a unique scaffolding-like protein with no obvious enzymatic fold. We postulate that XopN mediates its virulence function by associating with host proteins regulating defense. XopN was previously shown to physically interact with the tomato 14-3-3 protein TFT1 in planta. The significance of TFT1 in plant immune signaling and the relevance of the XopN/TFT1 interaction were not determined. Here we show that TFT1 is a positive regulator of PTI in tomato required to inhibit Xcv growth. Moreover, we provide evidence that TFT1 is a bona fide target of XopN because mutations that disrupt XopN/TFT1 binding also eliminate XopN-dependent virulence in tomato. This is the first example of a bacterial effector targeting a 14-3-3 associated with host immunity.

NEMO ensures signaling specificity of the pleiotropic IKKβ by directing its kinase activity toward IκBα

Besides activating NFκB by phosphorylating IκBs, IKKα/IKKβ kinases are also involved in regulating metabolic insulin signaling, the mTOR pathway, Wnt signaling, and autophagy. How IKKβ enzymatic activity is targeted to stimulus-specific substrates has remained unclear. We show here that NEMO, known to be essential for IKKβ activation by inflammatory stimuli, is also a specificity factor that directs IKKβ activity toward IκBα. Physical interaction and functional competition studies with mutant NEMO and IκB proteins indicate that NEMO functions as a scaffold to recruit IκBα to IKKβ. Interestingly, expression of NEMO mutants that allow for IKKβ activation by the cytokine IL-1, but fail to recruit IκBs, results in hyperphosphorylation of alternative IKKβ substrates. Furthermore IKK's function in autophagy, which is independent of NFκB, is significantly enhanced without NEMO as IκB scaffold. Our work establishes a role for scaffolds such as NEMO in determining stimulus-specific signal transduction via the pleiotropic signaling hub IKK.

Natural strategies for the spatial optimization of metabolism in synthetic biology

Metabolism is a highly interconnected web of chemical reactions that power life. Though the stoichiometry of metabolism is well understood, the multidimensional aspects of metabolic regulation in time and space remain difficult to define, model and engineer. Complex metabolic conversions can be performed by multiple species working cooperatively and exchanging metabolites via structured networks of organisms and resources. Within cells, metabolism is spatially regulated via sequestration in subcellular compartments and through the assembly of multienzyme complexes. Metabolic engineering and synthetic biology have had success in engineering metabolism in the first and second dimensions, designing linear metabolic pathways and channeling metabolic flux. More recently, engineering of the third dimension has improved output of engineered pathways through isolation and organization of multicell and multienzyme complexes. This review highlights natural and synthetic examples of three-dimensional metabolism both inter- and intracellularly, offering tools and perspectives for biological design.

Shoc2 is targeted to late endosomes and required for Erk1/2 activation in EGF-stimulated cells

Shoc2 is the putative scaffold protein that interacts with RAS and RAF, and positively regulates signaling to extracellular signal-regulated protein kinases 1 and 2 (ERK1/2). To elucidate the mechanism by which Shoc2 regulates ERK1/2 activation by the epidermal growth factor (EGF) receptor (EGFR), we studied subcellular localization of Shoc2. Upon EGFR activation, endogenous Shoc2 and red fluorescent protein tagged Shoc2 were translocated from the cytosol to a subset of late endosomes containing Rab7. The endosomal recruitment of Shoc2 was blocked by overexpression of a GDP-bound H-RAS (N17S) mutant and RNAi knockdown of clathrin, suggesting the requirement of RAS activity and clathrin-dependent endocytosis. RNAi depletion of Shoc2 strongly inhibited activation of ERK1/2 by low, physiological EGF concentrations, which was rescued by expression of wild-type recombinant Shoc2. In contrast, the Shoc2 (S2G) mutant, that is myristoylated and found in patients with the Noonan-like syndrome, did not rescue ERK1/2 activation in Shoc2-depleted cells. Shoc2 (S2G) was not located in late endosomes but was present on the plasma membrane and early endosomes. These data suggest that targeting of Shoc2 to late endosomes may facilitate EGFR-induced ERK activation under physiological conditions of cell stimulation by EGF, and therefore, may be involved in the spatiotemporal regulation of signaling through the RAS-RAF module.

Light-dependent, dark-promoted interaction between Arabidopsis cryptochrome 1 and phytochrome B proteins

Plant photoreceptors transduce environmental light cues to downstream signaling pathways, regulating a wide array of processes during growth and development. Two major plant photoreceptors with critical roles in photomorphogenesis are phytochrome B (phyB), a red/far-red absorbing photoreceptor, and cryptochrome 1 (CRY1), a UV-A/blue photoreceptor. Despite substantial genetic evidence for cross-talk between phyB and CRY1 pathways, a direct interaction between these proteins has not been observed. Here, we report that Arabidopsis phyB interacts directly with CRY1 in a light-dependent interaction. Surprisingly, the interaction is light-dissociated; CRY1 interacts specifically with the dark/far-red (Pr) state of phyB, but not with the red light-activated (Pfr) or the chromophore unconjugated form of the enzyme. The interaction is also regulated by light activation of CRY1; phyB Pr interacts only with the unstimulated form of CRY1 but not with the photostimulated protein. Further studies reveal that a small domain extending from the photolyase homology region (PHR) of CRY1 regulates the specificity of the interaction with different conformational states of phyB. We hypothesize that in plants, the phyB/CRY1 interaction may mediate cross-talk between the red/far-red- and blue/UV-sensing pathways, enabling fine-tuning of light responses to different spectral inputs.

Domain-mediated protein interaction prediction: From genome to network

Protein-protein interactions (PPIs), involved in many biological processes such as cellular signaling, are ultimately encoded in the genome. Solving the problem of predicting protein interactions from the genome sequence will lead to increased understanding of complex networks, evolution and human disease. We can learn the relationship between genomes and networks by focusing on an easily approachable subset of high-resolution protein interactions that are mediated by peptide recognition modules (PRMs) such as PDZ, WW and SH3 domains. This review focuses on computational prediction and analysis of PRM-mediated networks and discusses sequence- and structure-based interaction predictors, techniques and datasets for identifying physiologically relevant PPIs, and interpreting high-resolution interaction networks in the context of evolution and human disease.

Structural mechanism of disulphide bond-mediated redox switches

The oxidation of cysteine sulphydryl in proteins produces sulphenic acid that can form a reversible disulphide bond with another cysteine. The disulphide bond formation often triggers switches in protein structure and activity, especially when the distance between the two cysteine sulphur atoms is longer than the resulting disulphide bond distance. As an early example for the reversible disulphide bond-mediated functional switches, the reduced and oxidized forms of the bacterial transcription factor OxyR were characterized by X-ray crystallography. Recently, the Drosophila vision signalling protein, the association of inactivation-no-afterpotential D (INAD) was analysed by structural and functional methods. The two conserved cysteines of INAD were found to cycle between reduced and oxidized states during the light signal processing in Drosophila eyes, which was achieved by conformation dependent modulation of the disulphide bond redox potential. The production of the hypertension control peptide angiotensins was also shown to be controlled by the reversible disulphide bond in the precursor protein angiotensinogen. The crystal structure of the complex of angiotensiongen with its processing enzyme renin elucidated the role of the disulphide bond in stabilizing the precursor-enzyme complex facilitating the production of angiotensins. The increasing importance of the disulphide bond-mediated redox switches in normal and diseased states has implications in the development of novel antioxidant-based therapeutic approaches.

Beyond 'furballs' and 'dumpling soups' - towards a molecular architecture of signaling complexes and networks

The molecular architectures of intracellular signaling networks are largely unknown. Understanding their design principles and mechanisms of processing information is essential to grasp the molecular basis of virtually all biological processes. This is particularly challenging for human pathologies like cancers, as essentially each tumor is a unique disease with vastly deranged signaling networks. However, even in normal cells we know almost nothing. A few 'signalosomes', like the COP9 and the TCR signaling complexes have been described, but detailed structural information on their architectures is largely lacking. Similarly, many growth factor receptors, for example EGF receptor, insulin receptor and c-Met, signal via huge protein complexes built on large platform proteins (Gab, Irs/Dok, p130Cas[BCAR1], Frs families etc.), which are structurally not well understood. Subsequent higher order processing events remain even more enigmatic. We discuss here methods that can be employed to study signaling architectures, and the importance of too often neglected features like macromolecular crowding, intrinsic disorder in proteins and the sophisticated cellular infrastructures, which need to be carefully considered in order to develop a more mature understanding of cellular signal processing.

Structure of an enzyme-derived phosphoprotein recognition domain

Membrane Associated Guanylate Kinases (MAGUKs) contain a protein interaction domain (GK^dom) derived from the enzyme Guanylate Kinase (GK^enz). Here we show that GK^dom from the MAGUK Discs large (Dlg) is a phosphoprotein recognition domain, specifically recognizing the phosphorylated form of the mitotic spindle orientation protein Partner of Inscuteable (Pins). We determined the structure of the Dlg-Pins complex to understand the dramatic transition from nucleotide kinase to phosphoprotein recognition domain. The structure reveals that the region of the GK^dom that once served as the GMP binding domain (GBD) has been co-opted for protein interaction. Pins makes significantly more contact with the GBD than does GMP, but primarily with residues that are conserved between enzyme and domain revealing the versatility of the GBD as a platform for nucleotide and protein interactions. Mutational analysis reveals that the GBD is also used to bind the GK ligand MAP1a, suggesting that this is a common mode of MAGUK complex assembly. The GK^enz undergoes a dramatic closing reaction upon GMP binding but the protein-bound GK^dom remains in the ‘open’ conformation indicating that the dramatic conformational change has been lost in the conversion from nucleotide kinase to phosphoprotein recognition domain.

The emerging contribution of sequence context to the specificity of protein interactions mediated by PDZ domains

The canonical binding mode of PDZ domains to target motifs involves a small interface, unlikely to fully account for PDZ-target interaction specificities. Here, we review recent work on sequence context, defined as the regions surrounding not only the PDZ domains but also their target motifs. We also address the theoretical problem of defining the core of PDZ domains and the practical issue of designing PDZ constructs. Sequence context is found to introduce structural diversity, to impact the stability and solubility of constructs, and to deeply influence binding affinity and specificity, thereby increasing the difficulty of predicting PDZ-motif interactions. We expect that sequence context will have similar importance for other protein interactions mediated by globular domains binding to short linear motifs.

Motif switches: decision-making in cell regulation

Tight regulation of gene products from transcription to protein degradation is required for reliable and robust control of eukaryotic cell physiology. Many of the mechanisms directing cell regulation rely on proteins detecting the state of the cell through context-dependent, tuneable interactions. These interactions underlie the ability of proteins to make decisions by combining regulatory information encoded in a protein's expression level, localisation and modification state. This raises the question, how do proteins integrate available information to correctly make decisions? Over the past decade pioneering work on the nature and function of intrinsically disordered protein regions has revealed many elegant switching mechanisms that underlie cell signalling and regulation, prompting a reevaluation of their role in cooperative decision-making.

Genome-wide functional annotation of dual-specificity protein- and lipid-binding modules that regulate protein interactions

Emerging evidence indicates that membrane lipids regulate protein networking by directly interacting with protein-interaction domains (PIDs). As a pilot study to identify and functionally annodate lipid-binding PIDs on a genomic scale, we performed experimental and computational studies of PDZ domains. Characterization of 70 PDZ domains showed that ~40% had submicromolar membrane affinity. Using a computational model built from these data, we predicted the membrane-binding properties of 2,000 PDZ domains from 20 species. The accuracy of the prediction was experimentally validated for 26 PDZ domains. We also subdivided lipid-binding PDZ domains into three classes based on the interplay between membrane- and protein-binding sites. For different classes of PDZ domains, lipid binding regulates their protein interactions by different mechanisms. Functional studies of a PDZ domain protein, rhophilin 2, suggest that all classes of lipid-binding PDZ domains serve as genuine dual-specificity modules regulating protein interactions at the membrane under physiological conditions.

Rewiring kinase specificity with a synthetic adaptor protein

Signaling cascades are managed in time and space by interactions between and among proteins. These interactions are often aided by adaptor proteins, which guide enzyme-substrate pairs into proximity. Miniature proteins are a class of small, well-folded protein domains possessing engineered binding properties. Here we made use of two miniature proteins with complementary binding properties to create a synthetic adaptor protein that effectively redirects a ubiquitous signaling event: tyrosine phosphorylation. We report that miniature-protein-based adaptor 3 uses templated catalysis to redirect the Src family kinase Hck to phosphorylate hDM2, a negative regulator of the p53 tumor suppressor and a poor Hck substrate. Phosphorylation occurs with multiple turnover and at a single site targeted by c-Abl kinase in the cell.

Distinct peptide binding specificities of Src homology 3 (SH3) protein domains can be determined by modulation of local energetics across the binding interface

The yeast Nbp2p SH3 and Bem1p SH3b domains bind certain target peptides with similar high affinities, yet display vastly different affinities for other targets. To investigate this unusual behavior, we have solved the structure of the Nbp2p SH3-Ste20 peptide complex and compared it with the previously determined structure of the Bem1p SH3b bound to the same peptide. Although the Ste20 peptide interacts with both domains in a structurally similar manner, extensive in vitro studies with domain and peptide mutants revealed large variations in interaction strength across the binding interface of the two complexes. Whereas the Nbp2p SH3 made stronger contacts with the peptide core RXXPXXP motif, the Bem1p SH3b domain made stronger contacts with residues flanking the core motif. Remarkably, this modulation of local binding energetics can explain the distinct and highly nuanced binding specificities of these two domains.

The plasma membrane potential and the organization of the actin cytoskeleton of epithelial cells

The establishment and maintenance of the polarized epithelial phenotype require a characteristic organization of the cytoskeletal components. There are many cellular effectors involved in the regulation of the cytoskeleton of epithelial cells. Recently, modifications in the plasma membrane potential (PMP) have been suggested to participate in the modulation of the cytoskeletal organization of epithelia. Here, we review evidence showing that changes in the PMP of diverse epithelial cells promote characteristic modifications in the cytoskeletal organization, with a focus on the actin cytoskeleton. The molecular paths mediating these effects may include voltage-sensitive integral membrane proteins and/or peripheral proteins sensitive to surface potentials. The voltage dependence of the cytoskeletal organization seems to have implications in several physiological processes, including epithelial wound healing and apoptosis.

Structural determinants of the specificity of a membrane binding domain of the scaffold protein Ste5 of budding yeast: implications in signaling by the scaffold protein in MAPK pathway

In the mitogen activated protein kinase (MAPK) cascades of budding yeast, the scaffold protein Ste5 is recruited to the plasma membrane to transmit pheromone induced signal. A region or domain of Ste5 i.e. residues P44-R67, referred here as Ste5PM24, has been known to be involved in direct interactions with the membrane. In order to gain structural insights into membrane interactions of Ste5, here, we have investigated structures and interactions of two synthetic peptide fragments of Ste5, Ste5PM24, and a hyperactive mutant, Ste5PM24LM, by NMR, ITC, and fluorescence spectroscopy, with lipid membranes. We observed that Ste5PM24 predominantly interacted only with the anionic lipid vesicles. By contrast, Ste5PM24LM exhibited binding with negatively charged as well as zwitterionic or mixed lipid vesicles. Binding of Ste5 peptides with the negatively charged lipid vesicles were primarily driven by hydrophobic interactions. NMR studies revealed that Ste5PM24 assumes dynamic or transient conformations in zwitterionic dodecylphosphocholine (DPC) micelles. By contrast, NMR structure, obtained in anionic sodium dodecyl sulphate (SDS), demonstrated amphipathic helical conformations for the central segment of Ste5PM24. The hydrophobic surface of the helix was found to be buried inside the micelles. Taken together, these results provide important insights toward the structure and specificity determinants of the scaffold protein interactions with the plasma membrane.

Phosphosite mapping of P-type plasma membrane H+-ATPase in homologous and heterologous environments

Phosphorylation is an important posttranslational modification of proteins in living cells and primarily serves regulatory purposes. Several methods were employed for isolating phosphopeptides from proteolytically digested plasma membranes of Arabidopsis thaliana. After a mass spectrometric analysis of the resulting peptides we could identify 10 different phosphorylation sites in plasma membrane H(+)-ATPases AHA1, AHA2, AHA3, and AHA4/11, five of which have not been reported before, bringing the total number of phosphosites up to 11, which is substantially higher than reported so far for any other P-type ATPase. Phosphosites were almost exclusively (9 of 10) in the terminal regulatory domains of the pumps. The AHA2 isoform was subsequently expressed in the yeast Saccharomyces cerevisiae. The plant protein was phosphorylated at multiple sites in yeast, and surprisingly, seven of nine of the phosphosites identified in AHA2 were identical in the plant and fungal systems even though none of the target sequences in AHA2 show homology to proteins of the fungal host. These findings suggest an unexpected accessibility of the terminal regulatory domain of plasma membrane H(+)-ATPase to protein kinase action.

Unraveling the design principle for motif organization in signaling networks

Cellular signaling networks display complex architecture. Defining the design principle of this architecture is crucial for our understanding of various biological processes. Using a mathematical model for three-node feed-forward loops, we identify that the organization of motifs in specific manner within the network serves as an important regulator of signal processing. Further, incorporating a systemic stochastic perturbation to the model we could propose a possible design principle, for higher-order organization of motifs into larger networks in order to achieve specific biological output. The design principle was then verified in a large, complex human cancer signaling network. Further analysis permitted us to classify signaling nodes of the network into robust and vulnerable nodes as a result of higher order motif organization. We show that distribution of these nodes within the network at strategic locations then provides for the range of features displayed by the signaling network.

The Levinthal paradox of the interactome

The central biological question of the 21st century is: how does a viable cell emerge from the bewildering combinatorial complexity of its molecular components? Here, we estimate the combinatorics of self-assembling the protein constituents of a yeast cell, a number so vast that the functional interactome could only have emerged by iterative hierarchic assembly of its component sub-assemblies. A protein can undergo both reversible denaturation and hierarchic self-assembly spontaneously, but a functioning interactome must expend energy to achieve viability. Consequently, it is implausible that a completely "denatured" cell could be reversibly renatured spontaneously, like a protein. Instead, new cells are generated by the division of pre-existing cells, an unbroken chain of renewal tracking back through contingent conditions and evolving responses to the origin of life on the prebiotic earth. We surmise that this non-deterministic temporal continuum could not be reconstructed de novo under present conditions.

Evidence for alteration of calpain/calpastatin system in PBMC of cystic fibrosis patients

We are here reporting that in peripheral blood mononuclear cells (PBMC) of patients homozygous for F508del-CFTR the calpain-calpastatin system undergoes a profound alteration. In fact, calpain basal activity, almost undetectable in control PBMC, becomes measurable at a significant extent in cells from cystic fibrosis (CF) patients, also due to a 40-60% decrease in both calpastatin protein and inhibitory activity. Constitutive protease activation in CF patients' cells induces a large accumulation of the mutated cystic fibrosis transmembrane conductance regulator (CFTR) in the 100kD+70kD split forms as well as a degradation of proteins associated to the CFTR complex. Specifically, the scaffolding protein Na(+)/H(+) exchanger 3 regulatory factor-1 (NHERF-1) is converted in two distinct fragments showing masses of 35kD and 20kD, being however the latter form the most represented one, thereby indicating that in CF-PBMC the CFTR complex undergoes a large disorganization. In conclusion, our observations are providing new information on the role of calpain in the regulation of plasma membrane ion conductance and provide additional evidence on the transition of this protease activity from a physiological to a pathological function.

Characterization of PDZ domain-peptide interaction interface based on energetic patterns

PDZ domain is one of the abundant modular domains that recognize short peptide sequences to mediate protein-protein interactions. To decipher the binding specificity of PDZ domain, we analyzed the interactions between 11 mouse PDZ domains and 217 [corrected] peptides using a method called MIEC-SVM, which energetically characterizes the domain-peptide interaction using molecular interaction energy components (MIECs) and predicts binding specificity using support vector machine (SVM). Cross-validation and leave-one-domain-out test showed that the MIEC-SVM using all 44 PDZ-peptide residue pairs at the interaction interface outperformed the sequence-based methods in the literature. A further feature (residue pair) selection procedure illustrated that 16 residue pairs were uninformative to the binding specificity, even though they contributed significantly (~50%) to the binding energy. If only using the 28 informative residue pairs, the performance of the MIEC-SVM on predicting the PDZ binding specificity was significantly improved. This analysis suggests that the informative and uninformative residue interactions between the PDZ domain and the peptide may represent those contributing to binding specificity and affinity, respectively. We performed additional structural and energetic analyses to shed light on understanding how the PDZ-peptide recognition is established. The success of the MIEC-SVM method on PDZ domains in this study and SH3 domains in our previous studies illustrates its generality on characterizing protein-peptide interactions and understanding protein recognition from a structural and energetic viewpoint.

Evolution at the subgene level: domain rearrangements in the Drosophila phylogeny

Although the possibility of gene evolution by domain rearrangements has long been appreciated, current methods for reconstructing and systematically analyzing gene family evolution are limited to events such as duplication, loss, and sometimes, horizontal transfer. However, within the Drosophila clade, we find domain rearrangements occur in 35.9% of gene families, and thus, any comprehensive study of gene evolution in these species will need to account for such events. Here, we present a new computational model and algorithm for reconstructing gene evolution at the domain level. We develop a method for detecting homologous domains between genes and present a phylogenetic algorithm for reconstructing maximum parsimony evolutionary histories that include domain generation, duplication, loss, merge (fusion), and split (fission) events. Using this method, we find that genes involved in fusion and fission are enriched in signaling and development, suggesting that domain rearrangements and reuse may be crucial in these processes. We also find that fusion is more abundant than fission, and that fusion and fission events occur predominantly alongside duplication, with 92.5% and 34.3% of fusion and fission events retaining ancestral architectures in the duplicated copies. We provide a catalog of ∼9,000 genes that undergo domain rearrangement across nine sequenced species, along with possible mechanisms for their formation. These results dramatically expand on evolution at the subgene level and offer several insights into how new genes and functions arise between species.

Binding free energy landscape of domain-peptide interactions

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.

The complex biological processes occurring in living organisms are enabled by numerous networks of interacting proteins. It is therefore of great interest to understand the physical interplay between proteins and, in particular, how this process gives rise to highly specific network connectivities. For a long time, the dominant molecular view of protein-protein interactions was the docking of more or less static folded structures, with specificity obtained from a complementarity in shape and charge distributions. Lately it has been realized that many of the links in protein networks are mediated by interactions between folded domains, on the one hand, and disordered polypeptide segments, on the other. We use an all-atom Monte Carlo based approach which attempts to capture this domain-peptide binding process in full and apply it to representative members of a common domain family. This allows us to examine and compare detailed aspects of the binding free energy landscapes which underlie specificity and affinity. Being able to model domain-peptide binding in a physically sound, yet computationally tractable way is essential for identifying molecular binding mechanisms and opens up possibilities for modifying interaction networks in a controlled way.

Loss of EBP50 stimulates EGFR activity to induce EMT phenotypic features in biliary cancer cells

Scaffold proteins form multiprotein complexes that are central to the regulation of intracellular signaling. The scaffold protein ezrin-radixin-moesin-binding phosphoprotein 50 (EBP50) is highly expressed at the plasma membrane of normal biliary epithelial cells and binds epidermal growth factor receptor (EGFR), a tyrosine kinase receptor with oncogenic properties. This study investigated EBP50-EGFR interplay in biliary cancer. We report that in a collection of 106 cholangiocarcinomas, EBP50 was delocalized to the cytoplasm of tumor cells in 66% of the cases. Ectopic expression of EBP50 was correlated with the presence of satellite nodules and with the expression of EGFR, which was at the plasma membrane, implying a loss of interaction with EBP50 in these cases. In vitro, loss of interaction between EBP50 and EGFR was mimicked by EBP50 depletion using a small interfering RNA approach in human biliary carcinoma cells co-expressing the two proteins at their plasma membrane, and in which interaction between EBP50 and EGFR was validated. EBP50 depletion caused an increase in EGFR expression at their surface, and a sustained activation of the receptor and of its downstream effectors (extracellular signal-regulated kinase 1/2, signal transducer and activator of transcription 3) in both basal and EGF-stimulated conditions. Cells lacking EBP50 showed epithelial-to-mesenchymal transition-associated features, including reduction in E-cadherin and cytokeratin-19 expression, induction of S100A4 and of the E-cadherin transcriptional repressor, Slug, and loss of cell polarity. Accordingly, depletion of EBP50 induced the disruption of adherens junctional complexes, the development of lamellipodia structures and the subsequent acquisition of motility properties. All these phenotypic changes were prevented upon inhibition of EGFR tyrosine kinase by gefitinib. These findings indicate that loss of EBP50 at the plasma membrane in tumor cells may contribute to biliary carcinogenesis through EGFR activation.

Modularity of the bacterial cell cycle enables independent spatial and temporal control of DNA replication

BACKGROUND

Complex regulatory circuits in biology are often built of simpler subcircuits or modules. In most cases, the functional consequences and evolutionary origins of modularity remain poorly defined.

RESULTS

Here, by combining single-cell microscopy with genetic approaches, we demonstrate that two separable modules independently govern the temporal and spatial control of DNA replication in the asymmetrically dividing bacterium Caulobacter crescentus. DNA replication control involves DnaA, which promotes initiation, and CtrA, which silences initiation. We show that oscillations in DnaA activity dictate the periodicity of replication while CtrA governs the asymmetric replicative fates of daughter cells. Importantly, we demonstrate that DnaA activity oscillates independently of CtrA.

CONCLUSIONS

The genetic separability of spatial and temporal control modules in Caulobacter reflects their evolutionary history. DnaA is the central component of an ancient and phylogenetically widespread circuit that governs replication periodicity in Caulobacter and most other bacteria. By contrast, CtrA, which is found only in the asymmetrically dividing α-proteobacteria, was integrated later in evolution to enforce replicative asymmetry on daughter cells.

Domain orientation in the N-Terminal PDZ tandem from PSD-95 is maintained in the full-length protein

Tandem PDZ domains have been suggested to form structurally independent supramodules. However, dissimilarity between crystallography and NMR models emphasize their malleable conformation. Studies in full-length scaffold proteins are needed to examine the effect of tertiary interactions within their native context. Using single-molecule fluorescence to characterize the N-terminal PDZ tandem in PSD-95, we provide the first direct evidence that PDZ tandems can be structurally independent within a full-length scaffold protein. Molecular refinement using our data converged on a single structure with an antiparallel alignment of the ligand-binding sites. Devoid of interaction partners, single-molecule conditions captured PSD-95 in its unbound, ground state. Interactions between PDZ domains could not be detected while fluctuation correlation spectroscopy showed that other conformations are dynamically sampled. We conclude that ultra-weak interactions stabilize the conformation providing a "low-relief" energy landscape that allows the domain orientation to be flipped by environmental interactions.

Sequence-specific long range networks in PSD-95/discs large/ZO-1 (PDZ) domains tune their binding selectivity

Protein-protein interactions mediated by modular protein domains are critical for cell scaffolding, differentiation, signaling, and ultimately, evolution. Given the vast number of ligands competing for binding to a limited number of domain families, it is often puzzling how specificity can be achieved. Selectivity may be modulated by intradomain allostery, whereby a remote residue is energetically connected to the functional binding site via side chain or backbone interactions. Whereas several energetic pathways, which could mediate intradomain allostery, have been predicted in modular protein domains, there is a paucity of experimental data to validate their existence and roles. Here, we have identified such functional energetic networks in one of the most common protein-protein interaction modules, the PDZ domain. We used double mutant cycles involving site-directed mutagenesis of both the PDZ domain and the peptide ligand, in conjunction with kinetics to capture the fine energetic details of the networks involved in peptide recognition. We performed the analysis on two homologous PDZ-ligand complexes and found that the energetically coupled residues differ for these two complexes. This result demonstrates that amino acid sequence rather than topology dictates the allosteric pathways. Furthermore, our data support a mechanism whereby the whole domain and not only the binding pocket is optimized for a specific ligand. Such cross-talk between binding sites and remote residues may be used to fine tune target selectivity.

Structure collisions between interacting proteins

Protein-protein interactions take place at defined binding interfaces. One protein may bind two or more proteins at different interfaces at the same time. So far it has been commonly accepted that non-overlapping interfaces allow a given protein to bind other proteins simultaneously while no collisions occur between the binding protein structures. To test this assumption, we performed a comprehensive analysis of structural protein interactions to detect potential collisions. Our results did not indicate cases of biologically relevant collisions in the Protein Data Bank of protein structures. However, we discovered a number of collisions that originate from alternative protein conformations or quaternary structures due to different experimental conditions.

Recruitment interactions can override catalytic interactions in determining the functional identity of a protein kinase

The yeast Saccharomyces cerevisae has four distinct mitogen-activated protein kinase kinases (MAPKKs), each of which has a distinct functional identity characterized by communication with specific upstream and downstream partners to form distinct functional pathways. These four kinases belong to one family, sharing closely related catalytic domains. How have these four related kinases diverged to take on four distinct functional roles? The specificity of an enzyme for a particular substrate is often thought to reside in differences in the catalytic domain. However, many kinases, including MAPKKs, have modular interaction domains and motifs that have been shown to play an important role in determining the specificity of kinases through recruitment to specific partners and complexes. Here we probe the relative importance of catalytic domain interactions versus recruitment interactions in defining the functional identity of MAPKKs by asking whether we can use recruitment interactions to force other MAPKK catalytic domains to play the functional role of the mating MAPKK, Ste7. We find that two alternative MAPKKs, Pbs2 and Mkk2, can be forced to functionally replace the mating MAPKK Ste7, but only if the proper set of recruitment interactions are grafted onto their catalytic domains. These results show that within a family of kinases, recruitment interactions can play a dominant role in defining functional identity, and is consistent with a model in which new kinase functions can arise through recombination of existing catalytic domains with new interaction modules.

Protein structural modularity and robustness are associated with evolvability

Theory suggests that biological modularity and robustness allow for maintenance of fitness under mutational change, and when this change is adaptive, for evolvability. Empirical demonstrations that these traits promote evolvability in nature remain scant however. This is in part because modularity, robustness, and evolvability are difficult to define and measure in real biological systems. Here, we address whether structural modularity and/or robustness confer evolvability at the level of proteins by looking for associations between indices of protein structural modularity, structural robustness, and evolvability. We propose a novel index for protein structural modularity: the number of regular secondary structure elements (helices and strands) divided by the number of residues in the structure. We index protein evolvability as the proportion of sites with evidence of being under positive selection multiplied by the average rate of adaptive evolution at these sites, and we measure this as an average over a phylogeny of 25 mammalian species. We use contact density as an index of protein designability, and thus, structural robustness. We find that protein evolvability is positively associated with structural modularity as well as structural robustness and that the effect of structural modularity on evolvability is independent of the structural robustness index. We interpret these associations to be the result of reduced constraints on amino acid substitutions in highly modular and robust protein structures, which results in faster adaptation through natural selection.

Regioselective hydroformylation of allylic alcohols

A highly regioselective hydroformylation of allylic alcohols is reported toward the synthesis of β-hydroxy-acid and aldehyde products. The selectivity is achieved through the use of a ligand that reversibly binds to alcohols in situ, allowing for a directed hydroformylation to occur. The application to trisubstituted olefins was also demonstrated, which yields a single diastereomer product consistent with a stereospecific addition of CO and hydrogen.

Towards genome-scale signalling network reconstructions

Biological signalling networks allow living organisms to issue an integrated response to current conditions and make limited predictions about future environmental changes. Small-scale dynamic models of signalling cascades, including mitogen-activated protein kinase cascades, have been developed to generate hypotheses about signal transduction. Owing to technical limitations, these models and the hypotheses they generate have focused on a limited subset of signalling molecules. Now that we can simultaneously measure a substantial portion of the molecular components of a cell, we can begin to develop and test systems-level models of cellular signalling and regulatory processes, therefore gaining insights into the 'thought' processes of a cell.

Scaffold proteins: hubs for controlling the flow of cellular information

The spatial and temporal organization of molecules within a cell is critical for coordinating the many distinct activities carried out by the cell. In an increasing number of biological signaling processes, scaffold proteins have been found to play a central role in physically assembling the relevant molecular components. Although most scaffolds use a simple tethering mechanism to increase the efficiency of interaction between individual partner molecules, these proteins can also exert complex allosteric control over their partners and are themselves the target of regulation. Scaffold proteins offer a simple, flexible strategy for regulating selectivity in pathways, shaping output behaviors, and achieving new responses from preexisting signaling components. As a result, scaffold proteins have been exploited by evolution, pathogens, and cellular engineers to reshape cellular behavior.

Microbial linguistics: perspectives and applications of microbial cell-to-cell communication

Inter-cellular communication via diffusible small molecules is a defining character not only of multicellular forms of life but also of single-celled organisms. A large number of bacterial genes are regulated by the change of chemical milieu mediated by the local population density of its own species or others. The cell density-dependent "autoinducer" molecules regulate the expression of those genes involved in genetic competence, biofilm formation and persistence, virulence, sporulation, bioluminescence, antibiotic production, and many others. Recent innovations in recombinant DNA technology and micro-/nano-fluidics systems render the genetic circuitry responsible for cell-to-cell communication feasible to and malleable via synthetic biological approaches. Here we review the current understanding of the molecular biology of bacterial intercellular communication and the novel experimental protocols and platforms used to investigate this phenomenon. A particular emphasis is given to the genetic regulatory circuits that provide the standard building blocks which constitute the syntax of the biochemical communication network. Thus, this review gives focus to the engineering principles necessary for rewiring bacterial chemo-communication for various applications, ranging from population-level gene expression control to the study of host-pathogen interactions.

RNA templating the epigenome: long noncoding RNAs as molecular scaffolds

Cellular pathways must be synergized, controlled and organized to manage homeostasis. To achieve high selectivity within the crowded cellular milieu the cell utilizes scaffolding complexes whose role is to bring molecules in proximity thereby controlling and enhancing intermolecular interactions and signaling events. To date, scaffolds have been shown to be composed of proteinaceous units; however, recent evidence has supported the idea that non-coding RNAs may also play a similar role. In this point of view article we discuss recent data on ncRNA scaffolds, with particular focus on ncRNA HOTAIR. Using our current knowledge of signaling networks we discuss the role that RNA may play in writing and regulating histone modifications and the information needed for correct gene expression. Further, we speculate on additional, yet undiscovered roles that ncRNAs may be playing as molecular scaffolds.

A Conserved residue in the yeast Bem1p SH3 domain maintains the high level of binding specificity required for function

The yeast Bem1p SH3b and Nbp2p SH3 domains are unusual because they bind to peptides containing the same consensus sequence, yet they perform different functions and display low sequence similarity. In this work, by analyzing the interactions of these domains with six biologically relevant peptides containing the consensus sequence, they are shown to possess finely tuned and distinct binding specificities. We also identify a residue in the Bem1p SH3b domain that inhibits binding, yet is highly conserved for the purpose of preventing nonspecific interactions. Substitution of this residue results in a marked reduction of in vivo function that is caused by titration of the domain away from its proper targets through nonspecific interactions with other proteins. This work provides a clear illustration of the importance of intrinsic binding specificity for the function of protein-protein interaction modules, and the key role of "negative" interactions in determining the specificity of a domain.

The emergence of modularity in biological systems

In this review, we discuss modularity and hierarchy in biological systems. We review examples from protein structure, genetics, and biological networks of modular partitioning of the geometry of biological space. We review theories to explain modular organization of biology, with a focus on explaining how biology may spontaneously organize to a structured form. That is, we seek to explain how biology nucleated from among the many possibilities in chemistry. The emergence of modular organization of biological structure will be described as a symmetry-breaking phase transition, with modularity as the order parameter. Experimental support for this description will be reviewed. Examples will be presented from pathogen structure, metabolic networks, gene networks, and protein-protein interaction networks. Additional examples will be presented from ecological food networks, developmental pathways, physiology, and social networks.

Evolution of complex gene regulatory circuits by addition of refinements

How do complex gene regulatory circuits evolve? These circuits involve many interacting components, which work together to specify patterns of gene expression. They typically include many subtle mechanistic features, but in most cases it is unclear whether these features are essential for the circuit to work at all, or if instead they make a functional circuit work better. In the latter case, such a feature is here termed 'dispensable', and it is plausible that the feature has been added at a late stage in the evolution of the circuit. This review describes experimental tests of this question, using the phage λ gene regulatory circuit. Several features of this circuit are found to be dispensable, in the sense that the circuitry works without these features, though not as well as the wild type. In some cases, second-site suppressor mutations are needed to confer near-normal behavior in the absence of such a feature. These findings are discussed here in the context of a two-stage model for evolution of gene regulatory circuits. In this model, a circuit evolves by assembly of a primitive or basic form, followed by adjustment of parameters and addition of qualitatively new features. Pathways are suggested for the addition of such features to a more basic form. Selected examples in other systems are described. Some of the dispensable features of phage λ may be evolutionary refinements. Finding that a feature is dispensable, however, does not prove that it is a late addition - it is possible that it was essential early in evolution, and became dispensable as the circuit evolved. Conversely, a late addition might have become essential. As ongoing work provides additional examples of dispensable features, it may become clearer how often they represent refinements.

Evolutionary reshaping of fungal mating pathway scaffold proteins

Scaffold proteins play central roles in the function of many signaling pathways. Among the best-studied examples are the Ste5 and Far1 proteins of the yeast Saccharomyces cerevisiae. These proteins contain three conserved modules, the RING and PH domains, characteristic of some ubiquitin-ligating enzymes, and a vWA domain implicated in protein-protein interactions. In yeast, Ste5p regulates the mating pathway kinases while Far1p coordinates the cellular polarity machinery. Within the fungal lineage, the Basidiomycetes and the Pezizomycetes contain a single Far1-like protein, while several Saccharomycotina species, belonging to the CTG (Candida) clade, contain both a classic Far1-like protein and a Ste5-like protein that lacks the vWA domain. We analyzed the function of C. albicans Ste5p (Cst5p), a member of this class of structurally distinct Ste5 proteins. CST5 is essential for mating and still coordinates the mitogen-activated protein (MAP) kinase (MAPK) cascade elements in the absence of the vWA domain; Cst5p interacts with the MEK kinase (MEKK) C. albicans Ste11p (CaSte11p) and the MAPK Cek1 as well as with the MEK Hst7 in a vWA domain-independent manner. Cst5p can homodimerize, similar to Ste5p, but can also heterodimerize with Far1p, potentially forming heteromeric signaling scaffolds. We found direct binding between the MEKK CaSte11p and the MEK Hst7p that depends on a mobile acidic loop absent from S. cerevisiae Ste11p but related to the Ste7-binding region within the vWA domain of Ste5p. Thus, the fungal lineage has restructured specific scaffolding modules to coordinate the proteins required to direct the gene expression, polarity, and cell cycle regulation essential for mating.

Evolving specificity from variability for protein interaction domains

An important question in modular domain-peptide interactions, which play crucial roles in many biological processes, is how the diverse specificities exhibited by different members of a domain family are encoded in a common scaffold. Analysis of the Src homology (SH) 2 family has revealed that its specificity is determined, in large part, by the configuration of surface loops that regulate ligand access to binding pockets. In a distinct manner, SH3 domains employ loops for ligand recognition. The PDZ domain, in contrast, achieves specificity by co-evolution of binding-site residues. Thus, the conformational and sequence variability afforded by surface loops and binding sites provides a general mechanism by which to encode the wide spectrum of specificities observed for modular protein interaction domains.

Rule-based modelling and simulation of biochemical systems with molecular finite automata

The authors propose a theoretical formalism, molecular finite automata (MFA), to describe individual proteins as rule-based computing machines. The MFA formalism provides a framework for modelling individual protein behaviours and systems-level dynamics via construction of programmable and executable machines. Models specified within this formalism explicitly represent the context-sensitive dynamics of individual proteins driven by external inputs and represent protein-protein interactions as synchronised machine reconfigurations. Both deterministic and stochastic simulations can be applied to quantitatively compute the dynamics of MFA models. They apply the MFA formalism to model and simulate a simple example of a signal-transduction system that involves an MAP kinase cascade and a scaffold protein.