Changes of Cell Biochemical States Are Revealed in Protein Homomeric Complex Dynamics.

We report here a simple and global strategy to map out gene functions and target pathways of drugs, toxins, or other small molecules based on "homomer dynamics" protein-fragment complementation assays (hdPCA). hdPCA measures changes in self-association (homomerization) of over 3,500 yeast proteins in yeast grown under different conditions. hdPCA complements genetic interaction measurements while eliminating the confounding effects of gene ablation. We demonstrate that hdPCA accurately predicts the effects of two longevity and health span-affecting drugs, the immunosuppressant rapamycin and the type 2 diabetes drug metformin, on cellular pathways. We also discovered an unsuspected global cellular response to metformin that resembles iron deficiency and includes a change in protein-bound iron levels. This discovery opens a new avenue to investigate molecular mechanisms for the prevention or treatment of diabetes, cancers, and other chronic diseases of aging.

Mechanisms and Consequences of Macromolecular Phase Separation.

Over a century ago, colloidal phase separation of matter into non-membranous bodies was recognized as a fundamental organizing principal of cell "protoplasm." Recent insights into the molecular properties of such phase-separated bodies present challenges to our understanding of cellular protein interaction networks, as well as opportunities for interpreting and understanding of native and pathological genetic and molecular interactions. Here, we briefly review examples of and discuss physical principles of phase-separated cellular bodies and then reflect on how knowledge of these principles may direct future research on their functions.

Peroxisome biogenesis initiated by protein phase separation.

Peroxisomes are organelles that carry out ß-oxidation of fatty acids and amino acids. Both rare and prevalent diseases are caused by their dysfunction1. Among disease-causing variant genes are those required for protein transport into peroxisomes. The peroxisomal protein import machinery, which also shares similarities with chloroplasts2, is unique in transporting folded and large, up to 10 nm in diameter, protein complexes into peroxisomes3. Current models postulate a large pore formed by transmembrane proteins4; however, so far, no pore structure has been observed. In the budding yeast Saccharomyces cerevisiae, the minimum transport machinery includes the membrane proteins Pex13 and Pex14 and the cargo-protein-binding transport receptor, Pex5. Here we show that Pex13 undergoes liquid-liquid phase separation (LLPS) with Pex5-cargo. Intrinsically disordered regions in Pex13 and Pex5 resemble those found in nuclear pore complex proteins. Peroxisomal protein import depends on both the number and pattern of aromatic residues in these intrinsically disordered regions, consistent with their roles as 'stickers' in associative polymer models of LLPS5,6. Finally, imaging fluorescence cross-correlation spectroscopy shows that cargo import correlates with transient focusing of GFP-Pex13 and GFP-Pex14 on the peroxisome membrane. Pex13 and Pex14 form foci in distinct time frames, suggesting that they may form channels at different saturating concentrations of Pex5-cargo. Our findings lead us to suggest a model in which LLPS of Pex5-cargo with Pex13 and Pex14 results in transient protein transport channels7.

A nutrient-responsive pathway that determines M phase timing through control of B-cyclin mRNA stability.

The rate of cell-cycle progression must be tuned in response to nutrient levels to ensure that sufficient materials are synthesized to generate viable daughters. We report that accumulation of the yeast M phase B-cyclin CLB2 mRNA depends on assembly and activation of the heterogeneous nuclear RNA-binding protein (hnRNP) arginine methyltransferase Hmt1, which is promoted by the kinase Dbf2 and countered by the PP2A phosphatase Pph22. Activated Hmt1 methylates hnRNPs, which in turn stabilize CLB2 transcripts. Dbf2 activation of Hmt1 is highly cooperative, producing a sharp increase in CLB2, whereas Pph22 dephosphorylation is graded such that small changes in PP2A activity can cause large shifts in Dbf2-mediated Hmt1 activity. Starvation and rapamycin inhibition of TOR activate Pph22, causing a depletion of CLB2 and delay of M phase. We propose a general model wherein changes to Pph22 activity modulate cyclin mRNA stability to tune cell-cycle progression to environmental conditions.

Extracting insight from noisy cellular networks.

Network biologists attempt to extract meaningful relationships among genes or their products from very noisy data. We argue that what we categorize as noisy data may sometimes reflect noisy biology and therefore may shield a hidden meaning about how networks evolve and how matter is organized in the cell. We present practical solutions, based on existing evolutionary and biophysical concepts, through which our understanding of cell biology can be enormously enriched.

Endocytic proteins with prion-like domains form viscoelastic condensates that enable membrane remodeling.

Membrane invagination and vesicle formation are key steps in endocytosis and cellular trafficking. Here, we show that endocytic coat proteins with prion-like domains (PLDs) form hemispherical puncta in the budding yeast, Saccharomyces cerevisiae These puncta have the hallmarks of biomolecular condensates and organize proteins at the membrane for actin-dependent endocytosis. They also enable membrane remodeling to drive actin-independent endocytosis. The puncta, which we refer to as endocytic condensates, form and dissolve reversibly in response to changes in temperature and solution conditions. We find that endocytic condensates are organized around dynamic protein-protein interaction networks, which involve interactions among PLDs with high glutamine contents. The endocytic coat protein Sla1 is at the hub of the protein-protein interaction network. Using active rheology, we inferred the material properties of endocytic condensates. These experiments show that endocytic condensates are akin to viscoelastic materials. We use these characterizations to estimate the interfacial tension between endocytic condensates and their surroundings. We then adapt the physics of contact mechanics, specifically modifications of Hertz theory, to develop a quantitative framework for describing how interfacial tensions among condensates, the membrane, and the cytosol can deform the plasma membrane to enable actin-independent endocytosis.

Author Correction: Genome-wide SWAp-Tag yeast libraries for proteome exploration.

The version of Supplementary Table 1 originally published online with this article contained incorrect localization annotations for one plate. This error has been corrected in the online Supplementary Information.

Genetics' Piece of the PI: Inferring the Origin of Complex Traits and Diseases from Proteome-Wide Protein-Protein Interaction Dynamics.

How do common and rare genetic polymorphisms contribute to quantitative traits or disease risk and progression? Multiple human traits have been extensively characterized at the genomic level, revealing their complex genetic architecture. However, it is difficult to resolve the mechanisms by which specific variants contribute to a phenotype. Recently, analyses of variant effects on molecular traits have uncovered intermediate mechanisms that link sequence variation to phenotypic changes. Yet, these methods only capture a fraction of genetic contributions to phenotype. Here, in reviewing the field, it is proposed that complex traits can be understood by characterizing the dynamics of biochemical networks within living cells, and that the effects of genetic variation can be captured on these networks by using protein-protein interaction (PPI) methodologies. This synergy between PPI methodologies and the genetics of complex traits opens new avenues to investigate the molecular etiology of human diseases and to facilitate their prevention or treatment.

Genome-wide SWAp-Tag yeast libraries for proteome exploration.

Yeast libraries revolutionized the systematic study of cell biology. To extensively increase the number of such libraries, we used our previously devised SWAp-Tag (SWAT) approach to construct a genome-wide library of ~5,500 strains carrying the SWAT NOP1promoter-GFP module at the N terminus of proteins. In addition, we created six diverse libraries that restored the native regulation, created an overexpression library with a Cherry tag, or enabled protein complementation assays from two fragments of an enzyme or fluorophore. We developed methods utilizing these SWAT collections to systematically characterize the yeast proteome for protein abundance, localization, topology, and interactions.

Probing Inhomogeneous Diffusion in the Microenvironments of Phase-Separated Polymers under Confinement.

Biomolecular condensates formed by liquid-liquid phase separation of proteins and nucleic acids have been recently discovered to be prevalent in biology. These dynamic condensates behave like biochemical reaction vessels, but little is known about their structural organization and biophysical properties, which are likely related to condensate size. Thus, it is critical that we study them on scales found in vivo. However, previous in vitro studies of condensate assembly and physical properties have involved condensates up to 1000 times larger than those found in vivo. Here, we apply confinement microscopy to visualize condensates and control their sizes by creating appropriate confinement length scales relevant to the cell environment. We observe anomalous diffusion of probe particles embedded within confined condensates, as well as heterogeneous dynamics in condensates formed from PEG/dextran and in ribonucleoprotein complexes of RNA and the RNA-binding protein Dhh1. We propose that the observed non-Gaussian dynamics indicate a hopping diffusion mechanism inside condensates. We also observe that, for dextran-rich condensates, but not for ribonucleo condensates, probe particle diffusion depends on condensate size.

Induced and spontaneous colitis mouse models reveal complex interactions between IL-10 and IL-12/IL-23 pathways.

Crohn's disease (CD) and ulcerative colitis (UC) are the two major forms of inflammatory bowel disease (IBD). These idiopathic and chronic diseases result from inflammation of the gastrointestinal tract and are mainly mediated by the immune system. Genome wide association studies link genes of the IL-12 and IL-23 biology to both CD and UC susceptibility. IL-12 and IL-23 cytokines share a functional subunit, p40, and their respective receptors also share a functional subunit, IL-12Rß1. However, clinical trials targeting p40, and thus inhibiting both IL-12 and IL-23 pathways, provided mitigated effects on IBD, suggesting context dependent effects for each cytokine. In addition to IL-12 and IL-23, genetic deficiencies in IL-10 also result in severe IBD pathology. We generated various mouse models to determine how IL-12 or IL-23 interacts with IL-10 in IBD pathology. Whereas defects in both IL-10 and IL-12R do not impact the severity of the Dextran Sulfate Sodium (DSS)-induced colitis, combined deficiencies in both IL-10 and IL-23R aggravate the disease. In contrast to DSS-induced colitis, defects in IL-12R and IL-23R both protect from the spontaneous colitis observed in IL10-/- mice. Together, these studies exemplify the complexity of genetic and environmental interactions for identifying biological pathways predictive of pathological inflammatory processes.

Evolution of domain-peptide interactions to coadapt specificity and affinity to functional diversity.

Evolution of complexity in eukaryotic proteomes has arisen, in part, through emergence of modular independently folded domains mediating protein interactions via binding to short linear peptides in proteins. Over 30 years, structural properties and sequence preferences of these peptides have been extensively characterized. Less successful, however, were efforts to establish relationships between physicochemical properties and functions of domain-peptide interactions. To our knowledge, we have devised the first strategy to exhaustively explore the binding specificity of protein domain-peptide interactions. We applied the strategy to SH3 domains to determine the properties of their binding peptides starting from various experimental data. The strategy identified the majority (∼70%) of experimentally determined SH3 binding sites. We discovered mutual relationships among binding specificity, binding affinity, and structural properties and evolution of linear peptides. Remarkably, we found that these properties are also related to functional diversity, defined by depth of proteins within hierarchies of gene ontologies. Our results revealed that linear peptides evolved to coadapt specificity and affinity to functional diversity of domain-peptide interactions. Thus, domain-peptide interactions follow human-constructed gene ontologies, which suggest that our understanding of biological process hierarchies reflect the way chemical and thermodynamic properties of linear peptides and their interaction networks, in general, have evolved.

Exhaustive search of linear information encoding protein-peptide recognition.

High-throughput in vitro methods have been extensively applied to identify linear information that encodes peptide recognition. However, these methods are limited in number of peptides, sequence variation, and length of peptides that can be explored, and often produce solutions that are not found in the cell. Despite the large number of methods developed to attempt addressing these issues, the exhaustive search of linear information encoding protein-peptide recognition has been so far physically unfeasible. Here, we describe a strategy, called DALEL, for the exhaustive search of linear sequence information encoded in proteins that bind to a common partner. We applied DALEL to explore binding specificity of SH3 domains in the budding yeast Saccharomyces cerevisiae. Using only the polypeptide sequences of SH3 domain binding proteins, we succeeded in identifying the majority of known SH3 binding sites previously discovered either in vitro or in vivo. Moreover, we discovered a number of sites with both non-canonical sequences and distinct properties that may serve ancillary roles in peptide recognition. We compared DALEL to a variety of state-of-the-art algorithms in the blind identification of known binding sites of the human Grb2 SH3 domain. We also benchmarked DALEL on curated biological motifs derived from the ELM database to evaluate the effect of increasing/decreasing the enrichment of the motifs. Our strategy can be applied in conjunction with experimental data of proteins interacting with a common partner to identify binding sites among them. Yet, our strategy can also be applied to any group of proteins of interest to identify enriched linear motifs or to exhaustively explore the space of linear information encoded in a polypeptide sequence. Finally, we have developed a webserver located at http://michnick.bcm.umontreal.ca/dalel, offering user-friendly interface and providing different scenarios utilizing DALEL.

IL23R (Interleukin 23 Receptor) Variants Protective against Inflammatory Bowel Diseases (IBD) Display Loss of Function due to Impaired Protein Stability and Intracellular Trafficking.

Genome-wide association studies as well as murine models have shown that the interleukin 23 receptor (IL23R) pathway plays a pivotal role in chronic inflammatory diseases such as Crohn disease (CD), ulcerative colitis, psoriasis, and type 1 diabetes. Genome-wide association studies and targeted re-sequencing studies have revealed the presence of multiple potentially causal variants of the IL23R. Specifically the G149R, V362I, and R381Q IL23Rα chain variants are linked to protection against the development of Crohn disease and ulcerative colitis in humans. Moreover, the exact mechanism of action of these receptor variants has not been elucidated. We show that all three of these IL23Rα variants cause a reduction in IL23 receptor activation-mediated phosphorylation of the signal-transducing activator of transcription 3 (STAT3) and phosphorylation of signal transducing activator of transcription 4 (STAT4). The reduction in signaling is due to lower levels of cell surface receptor expression. For G149R, the receptor retention in the endoplasmic reticulum is due to an impairment of receptor maturation, whereas the R381Q and V362I variants have reduced protein stability. Finally, we demonstrate that the endogenous expression of IL23Rα protein from V362I and R381Q variants in human lymphoblastoid cell lines exhibited lower expression levels relative to susceptibility alleles. Our results suggest a convergent cause of IL23Rα variant protection against chronic inflammatory disease.

Protein Kinase C ζ Interacts with a Novel Binding Region of Gαq to Act as a Functional Effector.

Heterotrimeric G proteins play an essential role in the initiation of G protein-coupled receptor (GPCR) signaling through specific interactions with a variety of cellular effectors. We have recently reported that GPCR activation promotes a direct interaction between Gαq and protein kinase C ζ (PKCζ), leading to the stimulation of the ERK5 pathway independent of the canonical effector PLCß. We report herein that the activation-dependent Gαq/PKCζ complex involves the basic PB1-type II domain of PKCζ and a novel interaction module in Gαq different from the classical effector-binding site. Point mutations in this Gαq region completely abrogate ERK5 phosphorylation, indicating that Gαq/PKCζ association is required for the activation of the pathway. Indeed, PKCζ was demonstrated to directly bind ERK5 thus acting as a scaffold between Gαq and ERK5 upon GPCR activation. The inhibition of these protein complexes by G protein-coupled receptor kinase 2, a known Gαq modulator, led to a complete abrogation of ERK5 stimulation. Finally, we reveal that Gαq/PKCζ complexes link Gαq to apoptotic cell death pathways. Our data suggest that the interaction between this novel region in Gαq and the effector PKCζ is a key event in Gαq signaling.

An infrared reporter to detect spatiotemporal dynamics of protein-protein interactions.

We report a protein-fragment complementation assay (PCA) based on the engineered Deinococcus radiodurans infrared fluorescent protein IFP1.4. Unlike previous fluorescent protein PCAs, the IFP PCA is reversible, allowing analysis of spatiotemporal dynamics of hormone-induced signaling complexes in living yeast and mammalian cells at nanometer resolution. The inherently low background of infrared fluorescence permitted detection of subcellular reorganization of a signaling complex expressed at low abundance.

The scaffold protein Ste5 directly controls a switch-like mating decision in yeast.

Evolution has resulted in numerous innovations that allow organisms to increase their fitness by choosing particular mating partners, including secondary sexual characteristics, behavioural patterns, chemical attractants and corresponding sensory mechanisms. The haploid yeast Saccharomyces cerevisiae selects mating partners by interpreting the concentration gradient of pheromone secreted by potential mates through a network of mitogen-activated protein kinase (MAPK) signalling proteins. The mating decision in yeast is an all-or-none, or switch-like, response that allows cells to filter weak pheromone signals, thus avoiding inappropriate commitment to mating by responding only at or above critical concentrations when a mate is sufficiently close. The molecular mechanisms that govern the switch-like mating decision are poorly understood. Here we show that the switching mechanism arises from competition between the MAPK Fus3 and a phosphatase Ptc1 for control of the phosphorylation state of four sites on the scaffold protein Ste5. This competition results in a switch-like dissociation of Fus3 from Ste5 that is necessary to generate the switch-like mating response. Thus, the decision to mate is made at an early stage in the pheromone pathway and occurs rapidly, perhaps to prevent the loss of the potential mate to competitors. We argue that the architecture of the Fus3-Ste5-Ptc1 circuit generates a novel ultrasensitivity mechanism, which is robust to variations in the concentrations of these proteins. This robustness helps assure that mating can occur despite stochastic or genetic variation between individuals. The role of Ste5 as a direct modulator of a cell-fate decision expands the functional repertoire of scaffold proteins beyond providing specificity and efficiency of information processing. Similar mechanisms may govern cellular decisions in higher organisms and be disrupted in cancer.

Dissection of Cdk1-cyclin complexes in vivo.

Cyclin-dependent kinases (Cdks) are regulatory enzymes with temporal and spatial selectivity for their protein substrates that are governed by cell cycle-regulated cyclin subunits. Specific cyclin-Cdk complexes bind to and phosphorylate target proteins, coupling their activity to cell cycle states. The identification of specific cyclin-Cdk substrates is challenging and so far, has largely been achieved through indirect correlation or use of in vitro techniques. Here, we use a protein-fragment complementation assay based on the optimized yeast cytosine deaminase to systematically identify candidate substrates of budding yeast Saccharomyces cerevisiae Cdk1 and show dependency on one or more regulatory cyclins. We identified known and candidate cyclin dependencies for many predicted protein kinase Cdk1 targets and showed elusory Clb3-Cdk1-specific phosphorylation of γ-tubulin, thus establishing the timing of this event in controlling assembly of the mitotic spindle. Our strategy can be generally applied to identify substrates and accessory subunits of multisubunit protein complexes.