Hannes Braberg

Functional Organization of the S. cerevisiae Phosphorylation Network

Reversible protein phosphorylation is a signaling mechanism involved in all cellular processes. To create a systems view of the signaling apparatus in budding yeast, we generated an epistatic miniarray profile (E-MAP) comprised of 100,000 pairwise, quantitative genetic interactions, including virtually all protein and small-molecule kinases and phosphatases as well as key cellular regulators. Quantitative genetic interaction mapping reveals factors working in compensatory pathways (negative genetic interactions) or those operating in linear pathways (positive genetic interactions). We found an enrichment of positive genetic interactions between kinases, phosphatases, and their substrates. In addition, we assembled a higher-order map from sets of three genes that display strong interactions with one another: triplets enriched for functional connectivity. The resulting network view provides insights into signaling pathway regulation and reveals a link between the cell-cycle kinase, Cak1, the Fus3 MAP kinase, and a pathway that regulates chromatin integrity during transcription by RNA polymerase II.

A Genetic Interaction Map of RNA Processing Factors Reveals Links Between Sem1/Dss1-Containing Complexes and mRNA Export and Splicing

We used a quantitative, high-density genetic interaction map, or E-MAP (Epistatic MiniArray Profile), to interrogate the relationships within and between RNA-processing pathways. Due to their complexity and the essential roles of many of the components, these pathways have been difficult to functionally dissect. Here, we report the results for 107,155 individual interactions involving 552 mutations, 166 of which are hypomorphic alleles of essential genes. Our data enabled the discovery of links between components of the mRNA export and splicing machineries and Sem1/Dss1, a component of the 19S proteasome. In particular, we demonstrate that Sem1 has a proteasome-independent role in mRNA export as a functional component of the Sac3-Thp1 complex. Sem1 also interacts with Csn12, a component of the COP9 signalosome. Finally, we show that Csn12 plays a role in pre-mRNA splicing, which is independent of other signalosome components. Thus, Sem1 is involved in three separate and functionally distinct complexes.

High-throughput, quantitative analyses of genetic interactions in E. coli.

Large-scale genetic interaction studies provide the basis for defining gene function and pathway architecture. Recent advances in the ability to generate double mutants en masse in Saccharomyces cerevisiae have dramatically accelerated the acquisition of genetic interaction information and the biological inferences that follow. Here we describe a method based on F factor–driven conjugation, which allows for high-throughput generation of double mutants in Escherichia coli. This method, termed genetic interaction analysis technology for E. coli (GIANT-coli), permits us to systematically generate and array double-mutant cells on solid media in high-density arrays. We show that colony size provides a robust and quantitative output of cellular fitness and that GIANT-coli can recapitulate known synthetic interactions and identify previously unidentified negative (synthetic sickness or lethality) and positive (suppressive or epistatic) relationships. Finally, we describe a complementary strategy for genome-wide suppressor-mutant identification. Together, these methods permit rapid, large-scale genetic interaction studies in E. coli.

A plasma-membrane E-MAP reveals links of the eisosome with sphingolipid metabolism and endosomal trafficking

The plasma membrane delimits the cell and controls material and information exchange between itself and the environment. How different plasma-membrane processes are coordinated and how the relative abundance of plasma-membrane lipids and proteins is homeostatically maintained are not yet understood. Here, we used a quantitative genetic interaction map, or E-MAP, to functionally interrogate a set of ∼400 genes involved in various aspects of plasma-membrane biology, including endocytosis, signaling, lipid metabolism and eisosome function. From this E-MAP, we derived a set of 57,799 individual interactions between genes functioning in these various processes. Using triplet genetic motif analysis, we identified a new component of the eisosome, Eis1, and linked the poorly characterized gene EMP70 to endocytic and eisosome function. Finally, we implicated Rom2, a GDP/GTP exchange factor for Rho1 and Rho2, in the regulation of sphingolipid metabolism.

From Structure to Systems: High-Resolution, Quantitative Genetic Analysis of RNA Polymerase II

RNA polymerase II (RNAPII) lies at the core of dynamic control of gene expression. Using 53 RNAPII point mutants, we generated a point mutant epistatic miniarray profile (pE-MAP) comprising ∼60,000 quantitative genetic interactions in Saccharomyces cerevisiae. This analysis enabled functional assignment of RNAPII subdomains and uncovered connections between individual regions and other protein complexes. Using splicing microarrays and mutants that alter elongation rates in vitro, we found an inverse relationship between RNAPII speed and in vivo splicing efficiency. Furthermore, the pE-MAP classified fast and slow mutants that favor upstream and downstream start site selection, respectively. The striking coordination of polymerization rate with transcription initiation and splicing suggests that transcription rate is tuned to regulate multiple gene expression steps. The pE-MAP approach provides a powerful strategy to understand other multifunctional machines at amino acid resolution.

A lipid E-MAP identifies Ubx2 as a critical regulator of lipid saturation and lipid bilayer stress.

Biological membranes are complex, and the mechanisms underlying their homeostasis are incompletely understood. Here, we present a quantitative genetic interaction map (E-MAP) focused on various aspects of lipid biology, including lipid metabolism, sorting, and trafficking. This E-MAP contains ∼250,000 negative and positive genetic interaction scores and identifies a molecular crosstalk of protein quality control pathways with lipid bilayer homeostasis. Ubx2p, a component of the endoplasmic-reticulum-associated degradation pathway, surfaces as a key upstream regulator of the essential fatty acid (FA) desaturase Ole1p. Loss of Ubx2p affects the transcriptional control of OLE1, resulting in impaired FA desaturation and a severe shift toward more saturated membrane lipids. Both the induction of the unfolded protein response and aberrant nuclear membrane morphologies observed in cells lacking UBX2 are suppressed by the supplementation of unsaturated FAs. Our results point toward the existence of dedicated bilayer stress responses for membrane homeostasis.

Protein complex compositions predicted by structural similarity.

Proteins function through interactions with other molecules. Thus, the network of physical interactions among proteins is of great interest to both experimental and computational biologists. Here we present structure-based predictions of 3387 binary and 1234 higher order protein complexes in Saccharomyces cerevisiae involving 924 and 195 proteins, respectively. To generate candidate complexes, comparative models of individual proteins were built and combined together using complexes of known structure as templates. These candidate complexes were then assessed using a statistical potential, derived from binary domain interfaces in PIBASE (). The statistical potential discriminated a benchmark set of 100 interface structures from a set of sequence-randomized negative examples with a false positive rate of 3% and a true positive rate of 97%. Moreover, the predicted complexes were also filtered using functional annotation and sub-cellular localization data. The ability of the method to select the correct binding mode among alternates is demonstrated for three camelid VHH domain—porcine α–amylase interactions. We also highlight the prediction of co-complexed domain superfamilies that are not present in template complexes. Through integration with MODBASE, the application of the method to proteomes that are less well characterized than that of S.cerevisiae will contribute to expansion of the structural and functional coverage of protein interaction space. The predicted complexes are deposited in MODBASE ().

SALIGN: a web server for alignment of multiple protein sequences and structures

SUMMARY Accurate alignment of protein sequences and/or structures is crucial for many biological analyses, including functional annotation of proteins, classifying protein sequences into families, and comparative protein structure modeling. Described here is a web interface to SALIGN, the versatile protein multiple sequence/structure alignment module of MODELLER. The web server automatically determines the best alignment procedure based on the inputs, while allowing the user to override default parameter values. Multiple alignments are guided by a dendrogram computed from a matrix of all pairwise alignment scores. When aligning sequences to structures, SALIGN uses structural environment information to place gaps optimally. If two multiple sequence alignments of related proteins are input to the server, a profile-profile alignment is performed. All features of the server have been previously optimized for accuracy, especially in the contexts of comparative modeling and identification of interacting protein partners. AVAILABILITY The SALIGN web server is freely accessible to the academic community at http://salilab.org/salign. SALIGN is a module of the MODELLER software, also freely available to academic users (http://salilab.org/modeller). CONTACT sali@salilab.org; madhusudhan@bii.a-star.edu.sg.

RNA polymerase II carboxyl-terminal domain phosphorylation regulates protein stability of the Set2 methyltransferase and histone H3 di- and trimethylation at lysine 36.

Background: Set2 catalyzes mono-, di-, and trimethylation of lysine 36 on histone H3 (H3K36). Results: Phosphorylation of the C-terminal domain of RNA polymerase II (RNAPII CTD) regulates Set2 protein levels and H3K36 methylation. Conclusion: Set2 protein is regulated to maintain balanced levels of H3K36 methylation. Significance: These results suggest that different methylation states perform distinct biological functions. Methylation of lysine 36 on histone H3 (H3K36) is catalyzed by the Set2 methyltransferase and is linked to transcriptional regulation. Previous studies have shown that trimethylation of H3K36 by Set2 is directed through its association with the phosphorylated repeats of the RNA polymerase C-terminal domain (RNAPII CTD). Here, we show that disruption of this interaction through the use of yeast mutants defective in CTD phosphorylation at serine 2 results in a destabilization of Set2 protein levels and H3K36 methylation. Consistent with this, we find that Set2 has a short half-life and is co-regulated, with RNAPII CTD phosphorylation levels, during logarithmic growth in yeast. To probe the functional consequence of uncoupling Set2-RNAPII regulation, we expressed a truncated and more stable form of Set2 that is capable of dimethylation but not trimethylation in vivo. Results of high throughput synthetic genetic analyses show that this Set2 variant has distinct genetics from either SET2 or set2Δ and is synthetically sick or lethal with a number of transcription elongation mutants. Collectively, these results provide molecular insight into the regulation of Set2 protein levels that influence H3K36 methylation states.

Adventures in time and space: Splicing efficiency and RNA polymerase II elongation rate

Control of pre-mRNA splicing is a critical part of the eukaryotic gene expression process. Extensive evidence indicates that transcription and splicing are spatiotemporally coordinated and that most splicing events occur co-transcriptionally. A kinetic coupling model has been proposed in metazoans to describe how changing RNA Polymerase II (RNAPII) elongation rate can impact which alternative splice sites are used. In Saccharomyces cerevisiae, in which most spliced genes have only a single intron and splice sites adhere to a strong consensus sequence, we recently observed that splicing efficiency was sensitive to mutations in RNAPII that increase or decrease its elongation rate. Our data revealed that RNAPII speed and splicing efficiency are generally anti-correlated: at many genes, increased elongation rate caused decreased splicing efficiency, while decreased elongation rate increased splicing efficiency. An improved splicing phenotype was also observed upon deletion of SUB1, a condition in which elongation rate is slowed. We discuss these data in the context of a growing field and expand the kinetic coupling model to apply to both alternative splicing and splicing efficiency.